Cell-based splicing of minigenes

Deciphering RNA splicing logic with interpretable machine learning

Machine learning methods, particularly neural networks trained on large datasets, are transforming how scientists approach scientific discovery and experimental design. However, current state-of-the-art neural networks are limited by their uninterpretability: Despite their excellent accuracy, they cannot describe how they arrived at their predictions. Here, using an "interpretable-by-design" approach, we present a neural network model that provides insights into RNA splicing, a fundamental process in the transfer of genomic information into functional biochemical products. Although we designed our model to emphasize interpretability, its predictive accuracy is on par with state-of-the-art models. To demonstrate the model's interpretability, we introduce a visualization that, for any given exon, allows us to trace and quantify the entire decision process from input sequence to output splicing prediction. Importantly, the model revealed uncharacterized components of the splicing logic, which we experimentally validated. This study highlights how interpretable machine learning can advance scientific discovery.

Hypoparathyroidism, deafness and renal dysplasia syndrome caused by a GATA3 splice site mutation leading to the activation of a cryptic splice site

The HDR syndrome is a rare autosomal dominant disorder characterised by Hypoparathyroidism, Deafness, and Renal dysplasia, and is caused by inactivating heterozygous germline mutations in the GATA3 gene. We report an 11-year-old girl with HDR syndrome caused by a heterozygous mutation located at the splice acceptor site of exon 5 of the GATA3 gene (NM_001002295.2: c.925-1G>T). Functional studies using a minigene assay showed that this splice site mutation abolished the normal splicing of the GATA3 pre-mRNA and led to the use of a cryptic splice acceptor site, resulting in the loss of the first seven nucleotides (TCTGCAG) of exon 5 in the GATA3 mRNA. These findings increase the understanding of the mechanisms by which GATA3 splicing mutations can cause HDR syndrome.

Alternative splicing of HDAC7 regulates its interaction with 14-3-3 proteins to alter histone marks and target gene expression

Chromatin regulation and alternative splicing are both critical mechanisms guiding gene expression. Studies have demonstrated that histone modifications can influence alternative splicing decisions, but less is known about how alternative splicing may impact chromatin. Here, we demonstrate that several genes encoding histone-modifying enzymes are alternatively spliced downstream of T cell signaling pathways, including HDAC7, a gene previously implicated in controlling gene expression and differentiation in T cells. Using CRISPR-Cas9 gene editing and cDNA expression, we show that differential inclusion of HDAC7 exon 9 controls the interaction of HDAC7 with protein chaperones, resulting in changes to histone modifications and gene expression. Notably, the long isoform, which is induced by the RNA-binding protein CELF2, promotes expression of several critical T cell surface proteins including CD3, CD28, and CD69. Thus, we demonstrate that alternative splicing of HDAC7 has a global impact on histone modification and gene expression that contributes to T cell development.

Identification of seven exonic variants in the SLC4A1, ATP6V1B1, and ATP6V0A4 genes that alter RNA splicing by minigene assay

Primary distal renal tubular acidosis (dRTA) is a rare tubular disease associated with variants in SLC4A1, ATP6V0A4, ATP6V1B1, FOXⅠ1, or WDR72 genes. Currently, there is growing evidence that all types of exonic variants can alter splicing regulatory elements, affecting the precursor messenger RNA (pre-mRNA) splicing process. This study was to determine the consequences of variants associated with dRTA on pre-mRNA splicing combined with predictive bioinformatics tools and minigene assay. As a result, among the 15 candidate variants, 7 variants distributed in SLC4A1 (c.1765C>T, p.Arg589Cys), ATP6V1B1 (c.368G>T, p.Gly123Val; c.370C>T, p.Arg124Trp; c.484G>T, p.Glu162* and c.1102G>A, p.Glu368Lys) and ATP6V0A4 genes (c.322C>T, p.Gln108* and c.1572G>A, p.Pro524Pro) were identified to result in complete or incomplete exon skipping by either disruption of exonic splicing enhancers (ESEs) and generation of exonic splicing silencers, or interference with the recognition of the classic splicing site, or both. To our knowledge, this is the first study on pre-mRNA splicing of exonic variants in the dRTA-related genes. These results highlight the importance of assessing the effects of exonic variants at the mRNA level and suggest that minigene analysis is an effective tool for evaluating the effects of splicing on variants in vitro.

Reversion mutation of cDNA CA8-204 minigene construct produces a truncated functional peptide that regulates calcium release in vitro and produces profound analgesia in vivo

Intracellular calcium is critical in orchestrating neuronal excitability and analgesia. Carbonic anhydrase-8 (CA8) regulates intracellular calcium signaling through allosteric inhibition of neuronal inositol trisphosphate receptor 1 (ITPR1) to produce profound analgesia. Recently, we reported the "G" allele at rs6471859 represents cis-eQTL regulating alternative splicing of a 1697 bp transcript (CA8-204^G) with a retained intron, alternative polyadenylation site and a new stop codon producing a functional 26 kDa peptide with an extended exon 3. In this study we show the reversion mutation (G to C) at rs6471859 within the CA8-204^G expression vector also produced a stable 1697 bp transcript (CA8-204^C) coding for a smaller peptide (~ 22 kDa) containing only the first three CA8 exons. Surprisingly, this peptide inhibited ITPR1 (pITPR1) activation, ITPR1-mediated calcium release in vitro; and produced profound analgesia in vivo. This is the first report showing CA8-204^C codes for a functional peptide sufficient to regulate calcium signaling and produce profound analgesia.

Deep profiling and custom databases improve detection of proteoforms generated by alternative splicing

Alternative pre-mRNA splicing has long been proposed to contribute greatly to proteome complexity. However, the extent to which mature mRNA isoforms are successfully translated into protein remains controversial. Here, we used high-throughput RNA sequencing and mass spectrometry (MS)–based proteomics to better evaluate the translation of alternatively spliced mRNAs. To increase proteome coverage and improve protein quantitation, we optimized cell fractionation and sample processing steps at both the protein and peptide level. Furthermore, we generated a custom peptide database trained on analysis of RNA-seq data with MAJIQ, an algorithm optimized to detect and quantify differential and unannotated splice junction usage. We matched tandem mass spectra acquired by data-dependent acquisition (DDA) against our custom RNA-seq based database, as well as SWISS-PROT and RefSeq databases to improve identification of splicing-derived proteoforms by 28% compared with use of the SWISS-PROT database alone. Altogether, we identified peptide evidence for 554 alternate proteoforms corresponding to 274 genes. Our increased depth and detection of proteins also allowed us to track changes in the transcriptome and proteome induced by T-cell stimulation, as well as fluctuations in protein subcellular localization. In sum, our data here confirm that use of generic databases in proteomic studies underestimates the number of spliced mRNA isoforms that are translated into protein and provides a workflow that improves isoform detection in large-scale proteomic experiments.

A Multiplexed Assay for Exon Recognition Reveals that an Unappreciated Fraction of Rare Genetic Variants Cause Large-Effect Splicing Disruptions

Mutations that lead to splicing defects can have severe consequences on gene function and cause disease. Here, we explore how human genetic variation affects exon recognition by developing a multiplexed functional assay of splicing using Sort-seq (MFASS). We assayed 27,733 variants in the Exome Aggregation Consortium (ExAC) within or adjacent to 2,198 human exons in the MFASS minigene reporter and found that 3.8% (1,050) of variants, most of which are extremely rare, led to large-effect splice-disrupting variants (SDVs). Importantly, we find that 83% of SDVs are located outside of canonical splice sites, are distributed evenly across distinct exonic and intronic regions, and are difficult to predict a priori. Our results indicate extant, rare genetic variants can have large functional effects on splicing at appreciable rates, even outside the context of disease, and MFASS enables their empirical assessment at scale.

MAJIQ-SPEL: web-tool to interrogate classical and complex splicing variations from RNA-Seq data

Summary

Analysis of RNA sequencing (RNA-Seq) data have highlighted the fact that most genes undergo alternative splicing (AS) and that these patterns are tightly regulated. Many of these events are complex, resulting in numerous possible isoforms that quickly become difficult to visualize, interpret and experimentally validate. To address these challenges we developed MAJIQ-SPEL, a web-tool that takes as input local splicing variations (LSVs) quantified from RNA-Seq data and provides users with visualization and quantification of gene isoforms associated with those. Importantly, MAJIQ-SPEL is able to handle both classical (binary) and complex, non-binary, splicing variations. Using a matching primer design algorithm it also suggests to users possible primers for experimental validation by RT-PCR and displays those, along with the matching protein domains affected by the LSV, on UCSC Genome Browser for further downstream analysis.

Availability and implementation

Program and code will be available athttp://majiq.biociphers.org/majiq-spel.

Supplementary information

Supplementary data are available atBioinformatics online.

Double-target Antisense U1snRNAs Correct Mis-splicing Due to c.639+861C>T and c.639+919G>A GLA Deep Intronic Mutations

Fabry disease is a rare X-linked lysosomal storage disorder caused by deficiency of the α-galactosidase A (α-Gal A) enzyme, which is encoded by the GLA gene. GLA transcription in humans produces a major mRNA encoding α-Gal A and a minor mRNA of unknown function, which retains a 57-nucleotide-long cryptic exon between exons 4 and 5, bearing a premature termination codon. NM_000169.2:c.639+861C>T and NM_000169.2:c.639+919G>A GLA deep intronic mutations have been described to cause Fabry disease by inducing overexpression of the alternatively spliced mRNA, along with a dramatic decrease in the major one. Here, we built a wild-type GLA minigene and two minigenes that carry mutations c.639+861C>T and c.639+919G>A. Once transfected into cells, the minigenes recapitulate the molecular patterns observed in patients, at the mRNA, protein, and enzymatic level. We constructed a set of specific double-target U1asRNAs to correct c.639+861C>T and c.639+919G>A GLA mutations. Efficacy of U1asRNAs in inducing the skipping of the cryptic exon was evaluated upon their transient co-transfection with the minigenes in COS-1 cells, by real-time polymerase chain reaction (PCR), western blot analysis, and α-Gal A enzyme assay. We identified a set of U1asRNAs that efficiently restored α-Gal A enzyme activity and the correct splicing pathways in reporter minigenes. We also identified a unique U1asRNA correcting both mutations as efficently as the mutation-specific U1asRNAs. Our study proves that an exon skipping-based approach recovering α-Gal A activity in the c.639+861C>T and c.639+919G>A GLA mutations is active.

The Dengue Virus NS5 Protein Intrudes in the Cellular Spliceosome and Modulates Splicing

Dengue virus NS5 protein plays multiple functions in the cytoplasm of infected cells, enabling viral RNA replication and counteracting host antiviral responses. Here, we demonstrate a novel function of NS5 in the nucleus where it interferes with cellular splicing. Using global proteomic analysis of infected cells together with functional studies, we found that NS5 binds spliceosome complexes and modulates endogenous splicing as well as minigene-derived alternative splicing patterns. In particular, we show that NS5 alone, or in the context of viral infection, interacts with core components of the U5 snRNP particle, CD2BP2 and DDX23, alters the inclusion/exclusion ratio of alternative splicing events, and changes mRNA isoform abundance of known antiviral factors. Interestingly, a genome wide transcriptome analysis, using recently developed bioinformatics tools, revealed an increase of intron retention upon dengue virus infection, and viral replication was improved by silencing specific U5 components. Different mechanistic studies indicate that binding of NS5 to the spliceosome reduces the efficiency of pre-mRNA processing, independently of NS5 enzymatic activities. We propose that NS5 binding to U5 snRNP proteins hijacks the splicing machinery resulting in a less restrictive environment for viral replication.

Splicing Regulation of Pro-Inflammatory Cytokines and Chemokines: At the Interface of the Neuroendocrine and Immune Systems

Alternative splicing plays a key role in posttranscriptional regulation of gene expression, allowing a single gene to encode multiple protein isoforms. As such, alternative splicing amplifies the coding capacity of the genome enormously, generates protein diversity, and alters protein function. More than 90% of human genes undergo alternative splicing, and alternative splicing is especially prevalent in the nervous and immune systems, tissues where cells need to react swiftly and adapt to changes in the environment through carefully regulated mechanisms of cell differentiation, migration, targeting, and activation. Given its prevalence and complexity, this highly regulated mode of gene expression is prone to be affected by disease. In the following review, we look at how alternative splicing of signaling molecules—cytokines and their receptors—changes in different pathological conditions, from chronic inflammation to neurologic disorders, providing means of functional interaction between the immune and neuroendocrine systems. Switches in alternative splicing patterns can be very dynamic and can produce signaling molecules with distinct or antagonistic functions and localization to different subcellular compartments. This newly discovered link expands our understanding of the biology of immune and neuroendocrine cells, and has the potential to open new windows of opportunity for treatment of neurodegenerative disorders.