Massively parallel reporter assays discover de novo exonic splicing mutants in paralogs of Autism genes.

To determine the contribution of defective splicing in Autism Spectrum Disorders (ASD), the most common neurodevelopmental disorder, a high throughput Massively Parallel Splicing Assay (MaPSY) was employed and identified 42 exonic splicing mutants out of 725 coding de novo variants discovered in the sequencing of ASD families. A redesign of the minigene constructs in MaPSY revealed that upstream exons with strong 5' splice sites increase the magnitude of skipping phenotypes observed in downstream exons. Select hits were validated by RT-PCR and amplicon sequencing in patient cell lines. Exonic splicing mutants were enriched in probands relative to unaffected siblings -especially synonymous variants (7.5% vs 3.5%, respectively). Of the 26 genes disrupted by exonic splicing mutations, 6 were in known ASD genes and 3 were in paralogs of known ASD genes. Of particular interest was a synonymous variant in TNRC6C - an ASD gene paralog with interactions with other ASD genes. Clinical records of 3 ASD patients with TNRC6C variant revealed respiratory issues consistent with phenotypes observed in TNRC6 depleted mice. Overall, this study highlights the need for splicing analysis in determining variant pathogenicity, especially as it relates to ASD.

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.

Future directions for high-throughput splicing assays in precision medicine.

Classification of variants of unknown significance is a challenging technical problem in clinical genetics. As up to one-third of disease-causing mutations are thought to affect pre-mRNA splicing, it is important to accurately classify splicing mutations in patient sequencing data. Several consortia and healthcare systems have conducted large-scale patient sequencing studies, which discover novel variants faster than they can be classified. Here, we compare the advantages and limitations of several high-throughput splicing assays aimed at mitigating this bottleneck, and describe a data set of ~5,000 variants that we analyzed using our Massively Parallel Splicing Assay (MaPSy). The Critical Assessment of Genome Interpretation group (CAGI) organized a challenge, in which participants submitted machine learning models to predict the splicing effects of variants in this data set. We discuss the winning submission of the challenge (MMSplice) which outperformed existing software. Finally, we highlight methods to overcome the limitations of MaPSy and similar assays, such as tissue-specific splicing, the effect of surrounding sequence context, classifying intronic variants, synthesizing large exons, and amplifying complex libraries of minigene species. Further development of these assays will greatly benefit the field of clinical genetics, which lack high-throughput methods for variant interpretation.

Changes in the process of alternative RNA splicing results in soluble B and T lymphocyte attenuator with biological and clinical implications in critical illness.

BACKGROUND: Critically ill patients with sepsis and acute respiratory distress syndrome have severely altered physiology and immune system modifications. RNA splicing is a basic molecular mechanism influenced by physiologic alterations. Immune checkpoint inhibitors, such as B and T Lymphocyte Attenuator (BTLA) have previously been shown to influence outcomes in critical illness. We hypothesize altered physiology in critical illness results in alternative RNA splicing of the immune checkpoint protein, BTLA, resulting in a soluble form with biologic and clinical significance. METHODS: Samples were collected from critically ill humans and mice. Levels soluble BTLA (sBTLA) were measured. Ex vivo experiments assessing for cellular proliferation and cytokine production were done using splenocytes from critically ill mice cultured with sBTLA. Deep RNA sequencing was done to look for alternative splicing of BTLA. sBTLA levels were fitted to models to predict sepsis diagnosis. RESULTS: sBTLA is increased in the blood of critically ill humans and mice and can predict a sepsis diagnosis on hospital day 0 in humans. Alternative RNA splicing results in a premature stop codon that results in the soluble form. sBTLA has a clinically relevant impact as splenocytes from mice with critical illness cultured with soluble BTLA have increased cellular proliferation. CONCLUSION: sBTLA is produced as a result of alternative RNA splicing. This isoform of BTLA has biological significance through changes in cellular proliferation and can predict the diagnosis of sepsis.

Hereditary cancer genes are highly susceptible to splicing mutations.

Substitutions that disrupt pre-mRNA splicing are a common cause of genetic disease. On average, 13.4% of all hereditary disease alleles are classified as splicing mutations mapping to the canonical 5' and 3' splice sites. However, splicing mutations present in exons and deeper intronic positions are vastly underreported. A recent re-analysis of coding mutations in exon 10 of the Lynch Syndrome gene, MLH1, revealed an extremely high rate (77%) of mutations that lead to defective splicing. This finding is confirmed by extending the sampling to five other exons in the MLH1 gene. Further analysis suggests a more general phenomenon of defective splicing driving Lynch Syndrome. Of the 36 mutations tested, 11 disrupted splicing. Furthermore, analyzing past reports suggest that MLH1 mutations in canonical splice sites also occupy a much higher fraction (36%) of total mutations than expected. When performing a comprehensive analysis of splicing mutations in human disease genes, we found that three main causal genes of Lynch Syndrome, MLH1, MSH2, and PMS2, belonged to a class of 86 disease genes which are enriched for splicing mutations. Other cancer genes were also enriched in the 86 susceptible genes. The enrichment of splicing mutations in hereditary cancers strongly argues for additional priority in interpreting clinical sequencing data in relation to cancer and splicing.

Defective splicing of the RB1 transcript is the dominant cause of retinoblastomas.

Defective splicing is a common cause of genetic diseases. On average, 13.4% of all hereditary disease alleles are classified as splicing mutations with most mapping to the critical GT or AG nucleotides within the 5' and 3' splice sites. However, splicing mutations are underreported and the fraction of splicing mutations that compose all disease alleles varies greatly across disease gene. For example, there is a great excess (46%; ~threefold) of hereditary disease alleles that map to splice sites in RB1 that cause retinoblastoma. Furthermore, mutations in the exons and deeper intronic position may also affect splicing. We recently developed a high-throughput method that assays reported disease mutations for their ability to disrupt pre-mRNA splicing. Surprisingly, 27% of RB1-coding mutations tested also disrupt splicing. High-throughput in vitro spliceosomal assembly assay reveals heterogeneity in which stage of spliceosomal assembly is affected by splicing mutations. 58% of exonic splicing mutations were primarily blocked at the A complex in transition to the B complex and 33% were blocked at the B complex. Several mutants appear to reduce more than one step in the assembly. As RB1 splicing mutants are enriched in retinoblastoma disease alleles, additional priority should be allocated to this class of allele while interpreting clinical sequencing experiments. Analysis of the spectrum of RB1 variants observed in 60,706 exomes identifies 197 variants that have enough potential to disrupt splicing to warrant further consideration.

The effects of structure on pre-mRNA processing and stability.

Pre-mRNA molecules can form a variety of structures, and both secondary and tertiary structures have important effects on processing, function and stability of these molecules. The prediction of RNA secondary structure is a challenging problem and various algorithms that use minimum free energy, maximum expected accuracy and comparative evolutionary based methods have been developed to predict secondary structures. However, these tools are not perfect, and this remains an active area of research. The secondary structure of pre-mRNA molecules can have an enhancing or inhibitory effect on pre-mRNA splicing. An example of enhancing structure can be found in a novel class of introns in zebrafish. About 10% of zebrafish genes contain a structured intron that forms a bridging hairpin that enforces correct splice site pairing. Negative examples of splicing include local structures around splice sites that decrease splicing efficiency and potentially cause mis-splicing leading to disease. Splicing mutations are a frequent cause of hereditary disease. The transcripts of disease genes are significantly more structured around the splice sites, and point mutations that increase the local structure often cause splicing disruptions. Post-splicing, RNA secondary structure can also affect the stability of the spliced intron and regulatory RNA interference pathway intermediates, such as pre-microRNAs. Additionally, RNA secondary structure has important roles in the innate immune defense against viruses. Finally, tertiary structure can also play a large role in pre-mRNA splicing. One example is the G-quadruplex structure, which, similar to secondary structure, can either enhance or inhibit splicing through mechanisms such as creating or obscuring RNA binding protein sites.

Pathogenic variants that alter protein code often disrupt splicing.

The lack of tools to identify causative variants from sequencing data greatly limits the promise of precision medicine. Previous studies suggest that one-third of disease-associated alleles alter splicing. We discovered that the alleles causing splicing defects cluster in disease-associated genes (for example, haploinsufficient genes). We analyzed 4,964 published disease-causing exonic mutations using a massively parallel splicing assay (MaPSy), which showed an 81% concordance rate with splicing in patient tissue. Approximately 10% of exonic mutations altered splicing, mostly by disrupting multiple stages of spliceosome assembly. We present a large-scale characterization of exonic splicing mutations using a new technology that facilitates variant classification and keeps pace with variant discovery.

Residue-by-Residue View of In Vitro FUS Granules that Bind the C-Terminal Domain of RNA Polymerase II.

Phase-separated states of proteins underlie ribonucleoprotein (RNP) granules and nuclear RNA-binding protein assemblies that may nucleate protein inclusions associated with neurodegenerative diseases. We report that the N-terminal low-complexity domain of the RNA-binding protein Fused in Sarcoma (FUS LC) is structurally disordered and forms a liquid-like phase-separated state resembling RNP granules. This state directly binds the C-terminal domain of RNA polymerase II. Phase-separated FUS lacks static structures as probed by fluorescence microscopy, indicating they are distinct from both protein inclusions and hydrogels. We use solution nuclear magnetic resonance spectroscopy to directly probe the dynamic architecture within FUS liquid phase-separated assemblies. Importantly, we find that FUS LC retains disordered secondary structure even in the liquid phase-separated state. Therefore, we propose that disordered protein granules, even those made of aggregation-prone prion-like domains, are dynamic and disordered molecular assemblies with transiently formed protein-protein contacts.