Minigene splicing assays reveal new insights into exonic variants of the SLC12A3 gene in Gitelman Syndrome.

Gitelman syndrome (GS) is a kind of salt-losing tubular disease, most of which is caused by SLC12A3 gene variants, and missense variants account for the majority. Recently, the phenomenon of exon skipping, in which exonic variants disrupt normal pre-mRNA splicing, has been related to a variety of diseases. The purpose of this study was to identify the effect of previously presumed missense SLC12A3 variants on pre-mRNA splicing using bioinformatics tools and minigenes. The results revealed that, among ten candidate variants, six variants (c.602G>A, c.602G>T, c.1667C>T, c.1925G>A, c.2548G>C and c.2549G>C) led to complete or incomplete exon skipping by affecting exonic splicing regulatory elements and/or disturbing canonical splice sites. It is worth mentioning that this is the largest study on pre-mRNA splicing of SLC12A3 exonic variants. In addition, our study emphasizes the importance of detecting splicing function at the mRNA level in GS and indicates that minigene analysis is a valuable tool for splicing functional assays of variants in vitro.

Quantitative prediction of variant effects on alternative splicing using endogenous pre-messenger RNA structure probing

Splicing is a highly regulated process that depends on numerous factors. It is particularly challenging to quantitatively predict how a mutation will affect precursor messenger RNA (mRNA) structure and the subsequent functional consequences. Here we use a novel Mutational Profiling (-MaP) methodology to obtain highly reproducible endogenous precursor and mature mRNA structural probing data in vivo. We use these data to estimate Boltzmann suboptimal ensembles, and predict the structural consequences of mutations on precursor mRNA structure. Together with a structural analysis of recent cryo-EM spliceosome structures at different stages of the splicing cycle, we determined that the footprint of the Bact complex on precursor mRNA is best able to predict splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene. However, structure alone only achieves 74% accuracy. We therefore developed a β-regression weighting framework that incorporates splice site strength, structure and exonic/intronic splicing regulatory elements which together achieves 90% accuracy for 47 known and six newly discovered splice-altering variants. This combined experimental/computational framework represents a path forward for accurate prediction of splicing related disease-causing variants.

Role of Splicing Regulatory Elements and In Silico Tools Usage in the Identification of Deep Intronic Splicing Variants in Hereditary Breast/Ovarian Cancer Genes

The contribution of deep intronic splice-altering variants to hereditary breast and ovarian cancer (HBOC) is unknown. Current computational in silico tools to predict spliceogenic variants leading to pseudoexons have limited efficiency. We assessed the performance of the SpliceAI tool combined with ESRseq scores to identify spliceogenic deep intronic variants by affecting cryptic sites or splicing regulatory elements (SREs) using literature and experimental datasets. Our results with 233 published deep intronic variants showed that SpliceAI, with a 0.05 threshold, predicts spliceogenic deep intronic variants affecting cryptic splice sites, but is less effective in detecting those affecting SREs. Next, we characterized the SRE profiles using ESRseq, showing that pseudoexons are significantly enriched in SRE-enhancers compared to adjacent intronic regions. Although the combination of SpliceAI with ESRseq scores (considering ∆ESRseq and SRE landscape) showed higher sensitivity, the global performance did not improve because of the higher number of false positives. The combination of both tools was tested in a tumor RNA dataset with 207 intronic variants disrupting splicing, showing a sensitivity of 86%. Following the pipeline, five spliceogenic deep intronic variants were experimentally identified from 33 variants in HBOC genes. Overall, our results provide a framework to detect deep intronic variants disrupting splicing.

RAD51D Aberrant Splicing in Breast Cancer: Identification of Splicing Regulatory Elements and Minigene-Based Evaluation of 53 DNA Variants

RAD51D loss-of-function variants increase lifetime risk of breast and ovarian cancer. Splicing disruption is a frequent pathogenic mechanism associated with variants in susceptibility genes. Herein, we have assessed the splicing and clinical impact of splice-site and exonic splicing enhancer (ESE) variants identified through the study of ~113,000 women of the BRIDGES cohort. A RAD51D minigene with exons 2–9 was constructed in splicing vector pSAD. Eleven BRIDGES splice-site variants (selected by MaxEntScan) were introduced into the minigene by site-directed mutagenesis and tested in MCF-7 cells. The 11 variants disrupted splicing, collectively generating 25 different aberrant transcripts. All variants but one produced negligible levels (<3.4%) of the full-length (FL) transcript. In addition, ESE elements of the alternative exon 3 were mapped by testing four overlapping exonic microdeletions (≥30-bp), revealing an ESE-rich interval (c.202_235del) with critical sequences for exon 3 recognition that might have been affected by germline variants. Next, 26 BRIDGES variants and 16 artificial exon 3 single-nucleotide substitutions were also assayed. Thirty variants impaired splicing with variable amounts (0–65.1%) of the FL transcript, although only c.202G > A demonstrated a complete aberrant splicing pattern without the FL transcript. On the other hand, c.214T > C increased efficiency of exon 3 recognition, so only the FL transcript was detected (100%). In conclusion, 41 RAD51D spliceogenic variants (28 of which were from the BRIDGES cohort) were identified by minigene assays. We show that minigene-based mapping of ESEs is a powerful approach for identifying ESE hotspots and ESE-disrupting variants. Finally, we have classified nine variants as likely pathogenic according to ACMG/AMP-based guidelines, highlighting the complex relationship between splicing alterations and variant interpretation.

Systematic characterization of short intronic splicing-regulatory elements

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and MAPT minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.

Functional characterization of splicing regulatory elements

Background: RNA binding protein-RNA interactions mediate a variety of processes including pre-mRNA splicing, translation, decay, polyadenylation and many others. Previous high-throughput studies have characterized general sequence features associated with increased and decreased splicing of certain exons, but these studies are limited by not knowing the mechanisms, and in particular, the mediating RNA binding proteins, underlying these associations. Results: Here we utilize ENCODE data from diverse data modalities to identify functional splicing regulatory elements and their associated RNA binding proteins. We identify features which make splicing events more sensitive to depletion of RNA binding proteins, as well as which RNA binding proteins act as splicing regulators sensitive to depletion. To analyze the sequence determinants underlying RBP-RNA interactions impacting splicing, we assay tens of thousands of sequence variants in a high-throughput splicing reporter called Vex-seq and confirm a small subset in their endogenous loci using CRISPR base editors. Finally, we leverage other large transcriptomic datasets to confirm the importance of RNA binding proteins which we designed experiments around and identify additional RBPs which may act as additional splicing regulators of the exons studied. Conclusions: This study identifies sequence and other features underlying splicing regulation mediated specific RNA binding proteins, as well as validates and identifies other potentially important regulators of splicing in other large transcriptomic datasets.

A Comprehensive Analysis and Splicing Characterization of Naturally Occurring Synonymous Variants in the ATP7B Gene

Next-generation sequencing is effective for the molecular diagnosis of genetic diseases. However, the identification of the clinical significance of synonymous variants remains a challenge. Our previous study showed that some synonymous variants in ATP7B gene produced splicing disruptions, leading to Wilson disease (WD). To test the hypothesis that synonymous variants of ATP7B cause abnormal splicing by disrupting authentic splice sites or splicing regulatory elements, we used computational tools and minigene assays to characterize 253 naturally occurring ATP7B gene synonymous variants in this study. Human Splicing Finder (HSF) and ESE Finder 3.0 were used to predict the impact of these rare synonymous variants on pre-mRNA splicing. Then, we cloned 14 different wild-type Minigene_ATP7B_ex constructs for in vitro minigene assay, including 16 exons of ATP7B gene. After computational prediction, 85 candidate variants were selected to be introduced into the corresponding Minigene_ATP7B_ex constructs for splicing assays. Using this two-step procedure, we demonstrated that 11 synonymous variants in ExAc database (c.1620C&gt;T, c.3888C&gt;T, c.1554C&gt;T, c.1677C&gt;T, c.1830G&gt;A, c.1875T&gt;A, c.2826C&gt;A, c.4098G&gt;A, c.2994C&gt;T, c.3243G&gt;A, and c.3747G&gt;A) disrupted RNA splicing in vitro, and two (c.1620C&gt;T and c.3243G&gt;A) of these caused a complete exon skipping. The results not only provided a reliable experimental basis for the genetic diagnosis of WD patients but also offered some new insights into the pathogenicity of synonymous variants in genetic diseases.

Identification of eight exonic variants in the SLC4A1, ATP6V1B1 and ATP6V0A4 gene 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 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, 8 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*; c.1571C>T, p.Pro524Leu and c.1572G>A, p.Pro524Pro) were identified to result in whole or part of exon skipping by either disruption of ESEs and generation of ESSs, 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