Evidence in disease and non-disease contexts that nonsense mutations cause altered splicing via motif disruption

Selection on synonymous sites: the unwanted transcript hypothesis

Although translational selection to favour codons that match the most abundant tRNAs is not readily observed in humans, there is nonetheless selection in humans on synonymous mutations. We hypothesize that much of this synonymous site selection can be explained in terms of protection against unwanted RNAs - spurious transcripts, mis-spliced forms or RNAs derived from transposable elements or viruses. We propose not only that selection on synonymous sites functions to reduce the rate of creation of unwanted transcripts (for example, through selection on exonic splice enhancers and cryptic splice sites) but also that high-GC content (but low-CpG content), together with intron presence and position, is both particular to functional native mRNAs and used to recognize transcripts as native. In support of this hypothesis, transcription, nuclear export, liquid phase condensation and RNA degradation have all recently been shown to promote GC-rich transcripts and suppress AU/CpG-rich ones. With such 'traps' being set against AU/CpG-rich transcripts, the codon usage of native genes has, in turn, evolved to avoid such suppression. That parallel filters against AU/CpG-rich transcripts also affect the endosomal import of RNAs further supports the unwanted transcript hypothesis of synonymous site selection and explains the similar design rules that have enabled the successful use of transgenes and RNA vaccines.

Novel mutation leading to splice donor loss in a conserved site of DMD gene causes Duchenne muscular dystrophy with cryptorchidism

BACKGROUND

As one of the most common congenital abnormalities in male births, cryptorchidism has been found to have a polygenic aetiology according to previous studies of common variants. However, little is known about genetic predisposition of rare variants for cryptorchidism, since rare variants have larger effective size on diseases than common variants.

METHODS

In this study, a cohort of 115 Chinese probands with cryptorchidism was analysed using whole-genome sequencing, alongside 19 parental controls and 2136 unaffected men. Additionally, CRISPR-Cas9 editing of a conserved variant was performed in a mouse model, with MRI screening used to observe the phenotype.

RESULTS

In 30 of 115 patients (26.1%), we identified four novel genes (ARSH, DMD, MAGEA4 and SHROOM2) affecting at least five unrelated patients and four known genes (USP9Y, UBA1, BCORL1 and KDM6A) with the candidate rare pathogenic variants affecting at least two cases. Burden tests of rare variants revealed the genome-wide significances for newly identified genes (p<2.5×10-6) under the Bonferroni correction. Surprisingly, novel and known genes were mainly found on X chromosome (seven on X and one on Y) and all rare X-chromosomal segregating variants exhibited a maternal inheritance rather than de novo origin. CRISPR-Cas9 mouse modelling of a splice donor loss variant in DMD (NC_000023.11:g.32454661C>G), which resides in a conserved site across vertebrates, replicated bilateral cryptorchidism phenotypes, confirmed by MRI at 4 and 10 weeks. The movement tests further revealed symptoms of Duchenne muscular dystrophy (DMD) in transgenic mice.

CONCLUSION

Our results revealed the role of the DMD gene mutation in causing cryptorchidism. The results also suggest that maternal-X inheritance of pathogenic defects could have a predominant role in the development of cryptorchidism.

Assessing the reliability of point mutation as data augmentation for deep learning with genomic data

BACKGROUND

Deep neural networks (DNNs) have the potential to revolutionize our understanding and treatment of genetic diseases. An inherent limitation of deep neural networks, however, is their high demand for data during training. To overcome this challenge, other fields, such as computer vision, use various data augmentation techniques to artificially increase the available training data for DNNs. Unfortunately, most data augmentation techniques used in other domains do not transfer well to genomic data.

RESULTS

Most genomic data possesses peculiar properties and data augmentations may significantly alter the intrinsic properties of the data. In this work, we propose a novel data augmentation technique for genomic data inspired by biology: point mutations. By employing point mutations as substitutes for codons, we demonstrate that our newly proposed data augmentation technique enhances the performance of DNNs across various genomic tasks that involve coding regions, such as translation initiation and splice site detection.

CONCLUSION

Silent and missense mutations are found to positively influence effectiveness, while nonsense mutations and random mutations in non-coding regions generally lead to degradation. Overall, point mutation-based augmentations in genomic datasets present valuable opportunities for improving the accuracy and reliability of predictive models for DNA sequences.

Truncated variants of thyroid hormone receptor beta display disease-inflicting malfunctioning at cellular level

Thyroid hormone receptor β (THRβ) is a member of the nuclear receptor superfamily of ligand-modulated transcription factors. Upon ligand binding, THRβ sequentially recruits the components of transcriptional machinery to modulate target gene expression. In addition to regulating diverse physiological processes, THRβ plays a crucial role in hypothalamus-pituitary-thyroid axis feedback regulation. Anomalies in THRβ gene/protein structure are associated with onset of diverse disease states. In this study, we investigated disease-inflicting truncated variants of THRβ using in-silico analysis and cell-based assays. We examined the THRβ truncated variants on multiple test parameters, including subcellular localization, ligand-receptor interactions, transcriptional functions, interaction with heterodimeric partner RXR, and receptor-chromatin interactions. Moreover, molecular dynamic simulation approaches predicted that shortened THRβ-LBD due to point mutations contributes proportionally to the loss of structural integrity and receptor stability. Deviant subcellular localization and compromised transcriptional function were apparent with these truncated variants. Present study shows that 'mitotic bookmarking' property of some THRβ variants is also affected. The study highlights that structural and conformational attributes of THRβ are necessary for normal receptor functioning, and any deviations may contribute to the underlying cause of the inflicted diseases. We anticipate that insights derived herein may contribute to improved mechanistic understanding to assess disease predisposition.

Escaping from CRISPR-Cas-mediated knockout: the facts, mechanisms, and applications

Clustered regularly interspaced short palindromic repeats and associated Cas protein (CRISPR-Cas), a powerful genome editing tool, has revolutionized gene function investigation and exhibits huge potential for clinical applications. CRISPR-Cas-mediated gene knockout has already become a routine method in research laboratories. However, in the last few years, accumulating evidences have demonstrated that genes knocked out by CRISPR-Cas may not be truly silenced. Functional residual proteins could be generated in such knockout organisms to compensate the putative loss of function, termed herein knockout escaping. In line with this, several CRISPR-Cas-mediated knockout screenings have discovered much less abnormal phenotypes than expected. How does knockout escaping happen and how often does it happen have not been systematically reviewed yet. Without knowing this, knockout results could easily be misinterpreted. In this review, we summarize these evidences and propose two main mechanisms allowing knockout escaping. To avoid the confusion caused by knockout escaping, several strategies are discussed as well as their advantages and disadvantages. On the other hand, knockout escaping also provides convenient tools for studying essential genes and treating monogenic disorders such as Duchenne muscular dystrophy, which are discussed in the end.

Novel biallelic mutations in TMEM126B cause splicing defects and lead to Leigh-like syndrome with severe complex I deficiency

Leigh syndrome (LS)/Leigh-like syndrome (LLS) is one of the most common mitochondrial disease subtypes, caused by mutations in either the nuclear or mitochondrial genomes. Here, we identified a novel intronic mutation (c.82-2 A > G) and a novel exonic insertion mutation (c.290dupT) in TMEM126B from a Chinese patient with clinical manifestations of LLS. In silico predictions, minigene splicing assays and patients’ RNA analyses determined that the c.82-2 A > G mutation resulted in complete exon 2 skipping, and the c.290dupT mutation provoked partial and complete exon 3 skipping, leading to translational frameshifts and premature termination. Functional analysis revealed the impaired mitochondrial function in patient-derived lymphocytes due to severe complex I content and assembly defect. Altogether, this is the first report of LLS in a patient carrying mutations in TMEM126B. Our data uncovers the functional effect and the molecular mechanism of the pathogenic variants c.82-2 A > G and c.290dupT, which expands the gene mutation spectrum of LLS and clinical spectrum caused by TMEM126B mutations, and thus help to clinical diagnosis of TMEM126B mutation‐related mitochondrial diseases.

Predictive functional, statistical and structural analysis of CSNK2A1 and CSNK2B variants linked to neurodevelopmental diseases

Okur-Chung Neurodevelopmental Syndrome (OCNDS) and Poirier-Bienvenu Neurodevelopmental Syndrome (POBINDS) were recently identified as rare neurodevelopmental disorders. OCNDS and POBINDS are associated with heterozygous mutations in the CSNK2A1 and CSNK2B genes which encode CK2α, a serine/threonine protein kinase, and CK2β, a regulatory protein, respectively, which together can form a tetrameric enzyme called protein kinase CK2. A challenge in OCNDS and POBINDS is to understand the genetic basis of these diseases and the effect of the various CK2⍺ and CK2β mutations. In this study we have collected all variants available to date in CSNK2A1 and CSNK2B, and identified hotspots. We have investigated CK2⍺ and CK2β missense mutations through prediction programs which consider the evolutionary conservation, functionality and structure or these two proteins, compared these results with published experimental data on CK2α and CK2β mutants, and suggested prediction programs that could help predict changes in functionality of CK2α mutants. We also investigated the potential effect of CK2α and CK2β mutations on the 3D structure of the proteins and in their binding to each other. These results indicate that there are functional and structural consequences of mutation of CK2α and CK2β, and provide a rationale for further study of OCNDS and POBINDS-associated mutations. These data contribute to understanding the genetic and functional basis of these diseases, which is needed to identify their underlying mechanisms.

Twelve exonic variants in the SLC12A1 and CLCNKB genes alter RNA splicing in a minigene assay

Background: Bartter syndrome (BS) is a rare renal tubular disease caused by gene variants in SLC12A1, KCNJ1, CLCNKA, CLCNKB, BSND or MAGED2 genes. There is growing evidence that many exonic mutations can affect the pre-mRNA normal splicing and induce exon skipping by altering various splicing regulatory signals. Therefore, the aim of this study was to gain new insights into the consequences of exonic mutations associated with BS on pre-mRNA splicing.

Methods: We analyzed all the missense, nonsense and synonymous variants described in six pathogenic genes by bioinformatics programs and identified candidate mutations that may promote exon skipping through a minigene system.

Results: Results of the study showed that 12 of 14 candidate variants distributed in SLC12A1 (c.728G&gt;A, C.735C&gt;G, c.904C&gt;T, c.905G&gt;A, c.1304C&gt;T, c.1493C&gt;T, c.2221A&gt;T) and CLCNKB (c.226C&gt;T, c.228A&gt;C, c.229G&gt;A, c.229G&gt;C, c.1979C&gt;A) were identified to induce splicing alterations. These variants may not only disrupt exonic splicing enhancers (ESEs) but also generate new exonic splicing silencers (ESSs), or disturb the classic splicing sites.

Conclusion: To our knowledge, this is a comprehensive study regarding alterations in pre-mRNA of exonic variants in BS pathogenic genes. Our results reinforce the necessity of assessing the consequences of exonic variants at the mRNA level.

T-G-A Deficiency Pattern in Protein-Coding Genes and Its Potential Reason

If a stop codon appears within one gene, then its translation will be terminated earlier than expected. False folding of premature protein will be adverse to the host; hence, all functional genes would tend to avoid the intragenic stop codons. Therefore, we hypothesize that there will be less frequency of nucleotides corresponding to stop codons at each codon position of genes. Here, we validate this inference by investigating the nucleotide frequency at a large scale and results from 19,911 prokaryote genomes revealed that nucleotides coinciding with stop codons indeed have the lowest frequency in most genomes. Interestingly, genes with three types of stop codons all tend to follow a T-G-A deficiency pattern, suggesting that the property of avoiding intragenic termination pressure is the same and the major stop codon TGA plays a dominant role in this effect. Finally, a positive correlation between the TGA deficiency extent and the base length was observed in start-experimentally verified genes of Escherichia coli (E. coli). This strengthens the proof of our hypothesis. The T-G-A deficiency pattern observed would help to understand the evolution of codon usage tactics in extant organisms.

In rice splice variants that restore the reading frame after frameshifting indel introduction are common, often induced by the indels and sometimes lead to organism-level rescue

The introduction of frameshifting non-3n indels enables the identification of gene-trait associations. However, it has been hypothesised that recovery of the original reading frame owing to usage of non-canonical splice forms could cause rescue. To date there is very little evidence for organism-level rescue by such a mechanism and it is unknown how commonly indels induce, or are otherwise associated with, frame-restoring splice forms. We perform CRISPR/Cas9 editing of randomly selected loci in rice to investigate these issues. We find that the majority of loci have a frame-restoring isoform. Importantly, three quarters of these isoforms are not seen in the absence of the indels, consistent with indels commonly inducing novel isoforms. This is supported by analysis in the context of NMD knockdowns. We consider in detail the two top rescue candidates, in wax deficient anther 1 (wda1) and brittle culm (bc10), finding that organismal-level rescue in both cases is strong but owing to different splice modification routes. More generally, however, as frame-restoring isoforms are low abundance and possibly too disruptive, such rescue we suggest to be the rare exception, not the rule. Nonetheless, assuming that indels commonly induce frame-restoring isoforms, these results emphasize the need to examine RNA level effects of non-3n indels and suggest that multiple non-3n indels in any given gene are advisable to probe a gene’s trait associations.

As protein coding genes are read in units of three (codons), insertions or deletions (indels) that are not a multiple of three long (non 3n) are expected to be especially harmful. Whether they are is important both for interpreting the results of non-3n indel experiments to probe a gene’s functional importance and for diagnostics. Particularly enigmatic are incidences where some non-3n changes in a gene compromise phenotypes while other seemingly comparable ones do not. One explanation for the latter is that a non-3n indel might be rescued via a frame-restoring splice form. Here we examine this hypothesis by inducing non-3n indels in many genes in rice and find that many non-3n indels are associated with a splice form that restores the reading frame. In the majority of these cases the indel appears to induce the potential rescuing splice form. We examine two top hit cases in detail and show functional rescue by splice modification. More generally, the frame-restoring forms are, however, low abundance and probably result in compromised proteins. We conclude then that splice mediated rescue is possible, but probably uncommon. Nonetheless it should not be overlooked in experimental design and interpretation.