Estimating the prevalence of functional exonic splice regulatory information

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.

The Power of Clinical Diagnosis for Deciphering Complex Genetic Mechanisms in Rare Diseases

Complex genetic disease mechanisms, such as structural or non-coding variants, currently pose a substantial difficulty in frontline diagnostic tests. They thus may account for most unsolved rare disease patients regardless of the clinical phenotype. However, the clinical diagnosis can narrow the genetic focus to just a couple of genes for patients with well-established syndromes defined by prominent physical and/or unique biochemical phenotypes, allowing deeper analyses to consider complex genetic origin. Then, clinical-diagnosis-driven genome sequencing strategies may expedite the development of testing and analytical methods to account for complex disease mechanisms as well as to advance functional assays for the confirmation of complex variants, clinical management, and the development of new therapies.

Diagnostic potential of the amniotic fluid cells transcriptome in deciphering mendelian disease: a proof-of-concept

RNA sequencing (RNA-seq) is emerging in genetic diagnoses as it provides functional support for the interpretation of variants of uncertain significance. However, the use of amniotic fluid (AF) cells for RNA-seq has not yet been explored. Here, we examined the expression of clinically relevant genes in AF cells (n = 48) compared with whole blood and fibroblasts. The number of well-expressed genes in AF cells was comparable to that in fibroblasts and much higher than that in blood across different disease categories. We found AF cells RNA-seq feasible and beneficial in prenatal diagnosis (n = 4) as transcriptomic data elucidated the molecular consequence leading to the pathogenicity upgrade of variants in CHD7 and COL1A2 and revising the in silico prediction of a variant in MYRF. AF cells RNA-seq could become a reasonable choice for postnatal patients with advantages over fibroblasts and blood as it prevents invasive procedures.

Modulation of pre-mRNA structure by hnRNP proteins regulates alternative splicing of MALT1

Alternative splicing plays key roles for cell type-specific regulation of protein function. It is controlled by cis-regulatory RNA elements that are recognized by RNA binding proteins (RBPs). The MALT1 paracaspase is a key factor of signaling pathways that mediate innate and adaptive immune responses. Alternative splicing of MALT1 is critical for controlling optimal T cell activation. We demonstrate that MALT1 splicing depends on RNA structural elements that sequester the splice sites of the alternatively spliced exon7. The RBPs hnRNP U and hnRNP L bind competitively to stem-loop RNA structures that involve the 5' and 3' splice sites flanking exon7. While hnRNP U stabilizes RNA stem-loop conformations that maintain exon7 skipping, hnRNP L disrupts these RNA elements to facilitate recruitment of the essential splicing factor U2AF2, thereby promoting exon7 inclusion. Our data represent a paradigm for the control of splice site selection by differential RBP binding and modulation of pre-mRNA structure.

Effects of 14 F9 synonymous codon variants on hemophilia B expression: Alteration of splicing along with protein expression

There is growing evidence that synonymous codon variants (SCVs) can cause disease through the disruption of different processes of protein production. The aim of the study is to investigate whether the 14 SCVs reported in the F9 variant database were the pathogenic causes of hemophilia B. The impacts of SCVs on splicing and protein expression were detected using a combination of in silico prediction, in vitro minigene splicing assay and cell expression detection. The splicing transcripts were identified and quantified by co-amplification fluorescent PCR. The mechanism of splicing was verified by a modified pU1snRNA and pU7snRNA approach. Aberrant splicing patterns were found in eight SCVs. Five of the 8 SCVs produced almost all aberrant splicing isoforms, which were expected to truncate protein, three of them presented a partial defect on both splicing and protein secretion, the overall effects were consistent with the residual Factor IX activity of the affected cases. Neither the pre-messenger RNA (mRNA) splicing process nor the protein function was impaired in the rest six SCVs. In conclusion, our study firstly revealed the pathogenic mechanism of the 14 F9 SCVs and highlighted the importance of performing mRNA splicing analysis and protein expression studies of SCVs in inherited disorders.

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

Transcripts containing premature termination codons (PTCs) can be subject to nonsense-associated alternative splicing (NAS). Two models have been evoked to explain this, scanning and splice motif disruption. The latter postulates that exonic cis motifs, such as exonic splice enhancers (ESEs), are disrupted by nonsense mutations. We employ genome-wide transcriptomic and k-mer enrichment methods to scrutinize this model. First, we show that ESEs are prone to disruptive nonsense mutations owing to their purine richness and paucity of TGA, TAA and TAG. The motif model correctly predicts that NAS rates should be low (we estimate 5–30%) and approximately in line with estimates for the rate at which random point mutations disrupt splicing (8–20%). Further, we find that, as expected, NAS-associated PTCs are predictable from nucleotide-based machine learning approaches to predict splice disruption and, at least for pathogenic variants, are enriched in ESEs. Finally, we find that both in and out of frame mutations to TAA, TGA or TAG are associated with exon skipping. While a higher relative frequency of such skip-inducing mutations in-frame than out of frame lends some credence to the scanning model, these results reinforce the importance of considering splice motif modulation to understand the etiology of PTC-associated disease.

Splicing in the Diagnosis of Rare Disease: Advances and Challenges

Mutations which affect splicing are significant contributors to rare disease, but are frequently overlooked by diagnostic sequencing pipelines. Greater ascertainment of pathogenic splicing variants will increase diagnostic yields, ending the diagnostic odyssey for patients and families affected by rare disorders, and improving treatment and care strategies. Advances in sequencing technologies, predictive modeling, and understanding of the mechanisms of splicing in recent years pave the way for improved detection and interpretation of splice affecting variants, yet several limitations still prohibit their routine ascertainment in diagnostic testing. This review explores some of these advances in the context of clinical application and discusses challenges to be overcome before these variants are comprehensively and routinely recognized in diagnostics.

Large-scale comparative evaluation of user-friendly tools for predicting variant-induced alterations of splicing regulatory elements

Discriminating which nucleotide variants cause disease or contribute to phenotypic traits remains a major challenge in human genetics. In theory, any intragenic variant can potentially affect RNA splicing by altering splicing regulatory elements (SREs). However, these alterations are often ignored mainly because pioneer SRE predictors have proved inefficient. Here, we report the first large-scale comparative evaluation of four user-friendly SRE-dedicated algorithms (QUEPASA, HEXplorer, SPANR, and HAL) tested both as standalone tools and in multiple combined ways based on two independent benchmark datasets adding up to >1,300 exonic variants studied at the messenger RNA level and mapping to 89 different disease-causing genes. These methods display good predictive power, based on decision thresholds derived from the receiver operating characteristics curve analyses, with QUEPASA and HAL having the best accuracies either as standalone or in combination. Still, overall there was a tight race between the four predictors, suggesting that all methods may be of use. Additionally, QUEPASA and HEXplorer may be beneficial as well for predicting variant-induced creation of pseudoexons deep within introns. Our study highlights the potential of SRE predictors as filtering tools for identifying disease-causing candidates among the plethora of variants detected by high-throughput DNA sequencing and provides guidance for their use in genomic medicine settings.

Splicing impact of deep exonic missense variants in CAPN3 explored systematically by minigene functional assay

Improving the accuracy of variant interpretation during diagnostic sequencing is a major goal for genomic medicine. To explore an often-overlooked splicing effect of missense variants, we developed the functional assay ("minigene") for the majority of exons of CAPN3, the gene responsible for limb girdle muscular dystrophy. By systematically screening 21 missense variants distributed along the gene, we found that eight clinically relevant missense variants located at a certain distance from the exon-intron borders (deep exonic missense variants) disrupted normal splicing of CAPN3 exons. Several recent machine learning-based computational tools failed to predict splicing impact for the majority of these deep exonic missense variants, highlighting the importance of including variants of this type in the training sets during the future algorithm development. Overall, 24 variants in CAPN3 gene were explored, leading to the change in the American College of Medical Genetics and Genomics classification of seven of them when results of the "minigene" functional assay were considered. Our findings reveal previously unknown splicing impact of several clinically important variants in CAPN3 and draw attention to the existence of deep exonic variants with a disruptive effect on gene splicing that could be overlooked by the current approaches in clinical genetics.

Mutational bias and the protein code shape the evolution of splicing enhancers

Exonic splicing enhancers (ESEs) are enriched in exons relative to introns and bind splicing activators. This study considers a fundamental question of co-evolution: How did ESE motifs become enriched in exons prior to the evolution of ESE recognition? We hypothesize that the high exon to intron motif ratios necessary for ESE function were created by mutational bias coupled with purifying selection on the protein code. These two forces retain certain coding motifs in exons while passively depleting them from introns. Through the use of simulations, genomic analyses, and high throughput splicing assays, we confirm the key predictions of this hypothesis, including an overlap between protein and splicing information in ESEs. We discuss the implications of mutational bias as an evolutionary driver in other cis-regulatory systems.

Splicing is regulated by cis-acting elements in pre-mRNAs such as exonic or intronic splicing enhancers and silencers. Here the authors show that exonic splicing enhancers are enriched in exons compared to introns due to mutational bias coupled with purifying selection on the protein code.

An intronic splice site alteration in combination with a large deletion affecting VPS13B (COH1) causes Cohen syndrome

Cohen syndrome (CS) is a rare, autosomal recessive disorder characterized by intellectual disability, postnatal microcephaly, facial abnormalities, abnormal truncal fat distribution, myopia, and pigmentary retinopathy. It is often considered an underdiagnosed condition, especially in children with developmental delay and intellectual disability. Here we report on four individuals from a large Jordanian family clinically diagnosed with CS. Using Trio Exome Sequencing (Trio-WES) and MLPA analyses we identified a maternally inherited novel intronic nucleotide substitution c.3446-23T>G leading to the activation of a cryptic splice site and a paternally inherited multi-exon deletion in VPS13B (previously termed COH1) in the index patient. Expression analysis showed a strong decrease of VPS13B mRNA levels and direct sequencing of cDNA confirmed splicing at a cryptic upstream splice acceptor site, resulting in the inclusion of 22 intronic bases. This extension results in a frameshift and a premature stop of translation (p.Gly1149Valfs*9). Segregation analysis revealed that three affected maternal cousins were homozygous for the intronic splice site variant. Our data show causality of both alterations and strongly suggest the expansion of the diagnostic strategy to search for intronic splice variants in molecularly unconfirmed patients affected by CS.

Genetics' Piece of the PI: Inferring the Origin of Complex Traits and Diseases from Proteome-Wide Protein-Protein Interaction Dynamics

How do common and rare genetic polymorphisms contribute to quantitative traits or disease risk and progression? Multiple human traits have been extensively characterized at the genomic level, revealing their complex genetic architecture. However, it is difficult to resolve the mechanisms by which specific variants contribute to a phenotype. Recently, analyses of variant effects on molecular traits have uncovered intermediate mechanisms that link sequence variation to phenotypic changes. Yet, these methods only capture a fraction of genetic contributions to phenotype. Here, in reviewing the field, it is proposed that complex traits can be understood by characterizing the dynamics of biochemical networks within living cells, and that the effects of genetic variation can be captured on these networks by using protein-protein interaction (PPI) methodologies. This synergy between PPI methodologies and the genetics of complex traits opens new avenues to investigate the molecular etiology of human diseases and to facilitate their prevention or treatment.

Interplay between coding and exonic splicing regulatory sequences

The inclusion of exons during the splicing process depends on the binding of splicing factors to short low-complexity regulatory sequences. The relationship between exonic splicing regulatory sequences and coding sequences is still poorly understood. We demonstrate that exons that are coregulated by any given splicing factor share a similar nucleotide composition bias and preferentially code for amino acids with similar physicochemical properties because of the nonrandomness of the genetic code. Indeed, amino acids sharing similar physicochemical properties correspond to codons that have the same nucleotide composition bias. In particular, we uncover that the TRA2A and TRA2B splicing factors that bind to adenine-rich motifs promote the inclusion of adenine-rich exons coding preferentially for hydrophilic amino acids that correspond to adenine-rich codons. SRSF2 that binds guanine/cytosine-rich motifs promotes the inclusion of GC-rich exons coding preferentially for small amino acids, whereas SRSF3 that binds cytosine-rich motifs promotes the inclusion of exons coding preferentially for uncharged amino acids, like serine and threonine that can be phosphorylated. Finally, coregulated exons encoding amino acids with similar physicochemical properties correspond to specific protein features. In conclusion, the regulation of an exon by a splicing factor that relies on the affinity of this factor for specific nucleotide(s) is tightly interconnected with the exon-encoded physicochemical properties. We therefore uncover an unanticipated bidirectional interplay between the splicing regulatory process and its biological functional outcome.

Specific combinations of biallelic POLR3A variants cause Wiedemann-Rautenstrauch syndrome

BACKGROUND

Wiedemann-Rautenstrauch syndrome (WRS) is a form of segmental progeria presenting neonatally, characterised by growth retardation, sparse scalp hair, generalised lipodystrophy with characteristic local fatty tissue accumulations and unusual face. We aimed to understand its molecular cause.

METHODS

We performed exome sequencing in two families, targeted sequencing in 10 other families and performed in silico modelling studies and transcript processing analyses to explore the structural and functional consequences of the identified variants.

RESULTS

Biallelic POLR3A variants were identified in eight affected individuals and monoallelic variants of the same gene in four other individuals. In the latter, lack of genetic material precluded further analyses. Multiple variants were found to affect POLR3A transcript processing and were mostly located in deep intronic regions, making clinical suspicion fundamental to detection. While biallelic POLR3A variants have been previously reported in 4H syndrome and adolescent-onset progressive spastic ataxia, recurrent haplotypes specifically occurring in individuals with WRS were detected. All WRS-associated POLR3A amino acid changes were predicted to perturb substantially POLR3A structure/function.

CONCLUSION

Biallelic mutations in POLR3A, which encodes for the largest subunit of the DNA-dependent RNA polymerase III, underlie WRS. No isolated functional sites in POLR3A explain the phenotype variability in POLR3A-related disorders. We suggest that specific combinations of compound heterozygous variants must be present to cause the WRS phenotype. Our findings expand the molecular mechanisms contributing to progeroid disorders.

Exonic splice regulation imposes strong selection at synonymous sites

What proportion of coding sequence nucleotides have roles in splicing, and how strong is the selection that maintains them? Despite a large body of research into exonic splice regulatory signals, these questions have not been answered. This is because, to our knowledge, previous investigations have not explicitly disentangled the frequency of splice regulatory elements from the strength of the evolutionary constraint under which they evolve. Current data are consistent both with a scenario of weak and diffuse constraint, enveloping large swaths of sequence, as well as with well-defined pockets of strong purifying selection. In the former case, natural selection on exonic splice enhancers (ESEs) might primarily act as a slight modifier of codon usage bias. In the latter, mutations that disrupt ESEs are likely to have large fitness and, potentially, clinical effects. To distinguish between these scenarios, we used several different methods to determine the distribution of selection coefficients for new mutations within ESEs. The analyses converged to suggest that ∼15%-20% of fourfold degenerate sites are part of functional ESEs. Most of these sites are under strong evolutionary constraint. Therefore, exonic splice regulation does not simply impose a weak bias that gently nudges coding sequence evolution in a particular direction. Rather, the selection to preserve these motifs is a strong force that severely constrains the evolution of a substantial proportion of coding nucleotides. Thus synonymous mutations that disrupt ESEs should be considered as a potentially common cause of single-locus genetic disorders.

Decoding a cancer-relevant splicing decision in the RON proto-oncogene using high-throughput mutagenesis

Mutations causing aberrant splicing are frequently implicated in human diseases including cancer. Here, we establish a high-throughput screen of randomly mutated minigenes to decode the cis-regulatory landscape that determines alternative splicing of exon 11 in the proto-oncogene MST1R (RON). Mathematical modelling of splicing kinetics enables us to identify more than 1000 mutations affecting RON exon 11 skipping, which corresponds to the pathological isoform RON∆165. Importantly, the effects correlate with RON alternative splicing in cancer patients bearing the same mutations. Moreover, we highlight heterogeneous nuclear ribonucleoprotein H (HNRNPH) as a key regulator of RON splicing in healthy tissues and cancer. Using iCLIP and synergy analysis, we pinpoint the functionally most relevant HNRNPH binding sites and demonstrate how cooperative HNRNPH binding facilitates a splicing switch of RON exon 11. Our results thereby offer insights into splicing regulation and the impact of mutations on alternative splicing in cancer.

RES complex is associated with intron definition and required for zebrafish early embryogenesis

Pre-mRNA splicing is a critical step of gene expression in eukaryotes. Transcriptome-wide splicing patterns are complex and primarily regulated by a diverse set of recognition elements and associated RNA-binding proteins. The retention and splicing (RES) complex is formed by three different proteins (Bud13p, Pml1p and Snu17p) and is involved in splicing in yeast. However, the importance of the RES complex for vertebrate splicing, the intronic features associated with its activity, and its role in development are unknown. In this study, we have generated loss-of-function mutants for the three components of the RES complex in zebrafish and showed that they are required during early development. The mutants showed a marked neural phenotype with increased cell death in the brain and a decrease in differentiated neurons. Transcriptomic analysis of bud13, snip1 (pml1) and rbmx2 (snu17) mutants revealed a global defect in intron splicing, with strong mis-splicing of a subset of introns. We found these RES-dependent introns were short, rich in GC and flanked by GC depleted exons, all of which are features associated with intron definition. Using these features, we developed and validated a predictive model that classifies RES dependent introns. Altogether, our study uncovers the essential role of the RES complex during vertebrate development and provides new insights into its function during splicing.

Exonic mutations and exon skipping: Lessons learned from DFNA5

Dysregulation of splicing is a common factor underlying many inherited diseases including deafness. For one deafness-associated gene, DFNA5, perturbation of exon 8 splicing results in a constitutively active truncated protein. To date, only intronic mutations have been reported to cause exon 8 skipping in patients with DFNA5-related deafness. In five families with postlingual progressive autosomal dominant non-syndromic hearing loss, we employed two next-generation sequencing platforms-OtoSCOPE and whole exome sequencing-followed by variant filtering and prioritization based on both minor allele frequency and functional consequence using a customized bioinformatics pipeline to identify three novel and two recurrent mutations in DFNA5 that segregated with hearing loss in these families. The three novel mutations are all missense variants within exon 8 that are predicted computationally to decrease splicing efficiency or abolish it completely. We confirmed their functional impact in vitro using mini-genes carrying each mutant DFNA5 exon 8. In so doing, we present the first exonic mutations in DFNA5 to cause deafness, expand the mutational spectrum of DFNA5-related hearing loss, and highlight the importance of assessing the effect of coding variants on splicing.

Distinct functional consequences of ECEL1/DINE missense mutations in the pathogenesis of congenital contracture disorders

Endothelin-converting enzyme-like 1 (ECEL1, also termed DINE in rodents), a membrane-bound metalloprotease, has been identified as a gene responsible for distal arthrogryposis (DA). ECEL1-mutated DA is generally characterized by ocular phenotypes in addition to the congenital limb contractures that are common to all DA subtypes. Until now, the consequences of the identified pathogenic mutations have remained incompletely understood because of a lack of detailed phenotypic analyses in relevant mouse models. In this study, we generated a new knock-in mouse strain that carries an ECEL1/DINE pathogenic G607S missense mutation, based on a previous study reporting atypical DA hindlimb phenotypes in two siblings with the mutation. We compared the morphological phenotypes of G607S knock-in mice with C760R knock-in mice that we previously established. Both C760R and G607S knock-in mouse embryos showed similar axonal arborization defects with normal trajectory patterns from the spinal cord to the target hindlimb muscles, as well as axon guidance defects of the abducens nerves. Intriguingly, distinct phenotypes in DINE protein localization and mRNA expression were identified in these knock-in mouse lines. For G607S, DINE mRNA and protein expression was decreased or almost absent in motor neurons. In the C760R mutant mice DINE was expressed and localized in the somata of motor neurons but not in axons. Our mutant mouse data suggest that ECEL1/DINE G607S and C760R mutations both lead to motor innervation defects as primary causes in ECEL1-mutated congenital contracture disorders. However, the functional consequences of the two mutations are distinct, with loss of axonal transport of ECEL1/DINE in C760R mutants and mRNA expression deficits in G607S mutants.