CRISPR-dCas13d-based deep screening of proximal and distal splicing-regulatory elements.

Pre-mRNA splicing, a key process in gene expression, can be therapeutically modulated using various drug modalities, including antisense oligonucleotides (ASOs). However, determining promising targets is hampered by the challenge of systematically mapping splicing-regulatory elements (SREs) in their native sequence context. Here, we use the catalytically inactive CRISPR-RfxCas13d RNA-targeting system (dCas13d/gRNA) as a programmable platform to bind SREs and modulate splicing by competing against endogenous splicing factors. SpliceRUSH, a high-throughput screening method, was developed to map SREs in any gene of interest using a lentivirus gRNA library that tiles the genetic region, including distal intronic sequences. When applied to SMN2, a therapeutic target for spinal muscular atrophy, SpliceRUSH robustly identifies not only known SREs but also a previously unknown distal intronic SRE, which can be targeted to alter exon 7 splicing using either dCas13d/gRNA or ASOs. This technology enables a deeper understanding of splicing regulation with applications for RNA-based drug discovery.

PRPF8-mediated dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects.

The carboxy-terminus of the spliceosomal protein PRPF8, which regulates the RNA helicase Brr2, is a hotspot for mutations causing retinitis pigmentosa-type 13, with unclear role in human splicing and tissue-specificity mechanism. We used patient induced pluripotent stem cells-derived cells, carrying the heterozygous PRPF8 c.6926 A > C (p.H2309P) mutation to demonstrate retinal-specific endophenotypes comprising photoreceptor loss, apical-basal polarity and ciliary defects. Comprehensive molecular, transcriptomic, and proteomic analyses revealed a role of the PRPF8/Brr2 regulation in 5'-splice site (5'SS) selection by spliceosomes, for which disruption impaired alternative splicing and weak/suboptimal 5'SS selection, and enhanced cryptic splicing, predominantly in ciliary and retinal-specific transcripts. Altered splicing efficiency, nuclear speckles organisation, and PRPF8 interaction with U6 snRNA, caused accumulation of active spliceosomes and poly(A)+ mRNAs in unique splicing clusters located at the nuclear periphery of photoreceptors. Collectively these elucidate the role of PRPF8/Brr2 regulatory mechanisms in splicing and the molecular basis of retinal disease, informing therapeutic approaches.

Loss-of-function mutation in PRMT9 causes abnormal synapse development by dysregulation of RNA alternative splicing.

Protein arginine methyltransferase 9 (PRMT9) is a recently identified member of the PRMT family, yet its biological function remains largely unknown. Here, by characterizing an intellectual disability associated PRMT9 mutation (G189R) and establishing a Prmt9 conditional knockout (cKO) mouse model, we uncover an important function of PRMT9 in neuronal development. The G189R mutation abolishes PRMT9 methyltransferase activity and reduces its protein stability. Knockout of Prmt9 in hippocampal neurons causes alternative splicing of ~1900 genes, which likely accounts for the aberrant synapse development and impaired learning and memory in the Prmt9 cKO mice. Mechanistically, we discover a methylation-sensitive protein-RNA interaction between the arginine 508 (R508) of the splicing factor 3B subunit 2 (SF3B2), the site that is exclusively methylated by PRMT9, and the pre-mRNA anchoring site, a cis-regulatory element that is critical for RNA splicing. Additionally, using human and mouse cell lines, as well as an SF3B2 arginine methylation-deficient mouse model, we provide strong evidence that SF3B2 is the primary methylation substrate of PRMT9, thus highlighting the conserved function of the PRMT9/SF3B2 axis in regulating pre-mRNA splicing.

Phenotypic Description of A Patient with ODLURO Syndrome and Functional Characterization of the Pathogenetic Role of A Synonymous Variant c.186G>A in KMT2E Gene.

O'Donnell-Luria-Rodan (ODLURO) syndrome is an autosomal dominant disorder caused by mutations in the KMT2E gene. The clinical phonotype of the affected individuals is typically characterized by global developmental delay, autism, epilepsy, hypotonia, macrocephaly, and very mild dysmorphic facial features. In this report, we describe the case of a 6-year-old boy with ODLURO syndrome who is a carrier of the synonymous mutation c.186G>A (p.Ala62=) in the KMT2E gene, predicted to alter splicing by in silico tools. Given the lack of functional studies on the c.186G>A variant, in order to assess its potential functional effect, we sequenced the patient's cDNA demonstrating its impact on the mechanism of splicing. To the best of our knowledge, our patient is the second to date reported carrying this synonymous mutation, but he is the first whose functional investigation has confirmed the deleterious consequence of the variant, resulting in exon 4 skipping. Additionally, we suggest a potential etiological mechanism that could be responsible for the aberrant splicing mechanism in KMT2E.

Identification of a novel intronic variant of ATP6V0A2 in a Han-Chinese family with cutis laxa.

BACKGROUND: Cutis laxa is a connective tissue disease caused by abnormal synthesis or secretion of skin elastic fibers, leading to skin flabby and saggy in various body parts. It can be divided into congenital cutis laxa and acquired cutis laxa, and inherited cutis laxa syndromes is more common in clinic. METHODS: In this study, we reported a case of a Han-Chinese male newborn with ATP6V0A2 gene variant leading to cutis laxa. The proband was identified by whole-exome sequencing to determine the novel variant, and their parents were verified by Sanger sequencing. Bioinformatics analysis and minigene assay were used to verify the effect of this variant on splicing function. RESULTS: The main manifestations of the proband are skin laxity, abnormal facial features, and enlargement of the anterior fontanelle. Whole-exome sequencing showed that the newborn carried a non-canonical splicing-site variant c.117 + 5G > T, p. (?) in ATP6V0A2 gene. Sanger sequencing showed that both parents of the proband carried the heterozygous variant. The results of bioinformatics analysis and minigene assay displayed that the variant site affected the splicing function of pre-mRNA of the ATP6V0A2 gene. CONCLUSIONS: In this study, it was identified that ATP6V0A2 gene c. 117 + 5G > T may be the cause of the disease. The non-canonical splicing variants of ATP6V0A2 gene were rarely reported in the past, and this variant expanded the variants spectrum of the gene. The functional study of minigene assay plays a certain role in improving the level of evidence for the pathogenicity of splicing variants, which lays a foundation for prenatal counseling and follow-up gene therapy.

Utility of oligonucleotide in upregulating circular RNA production in a cellular model.

Circular RNAs (circRNAs), are a covalently closed, single-stranded RNA without 5'- and 3'-termini, commonly stem from the exons of precursor mRNAs (pre-mRNAs). They have recently garnered interest, with studies uncovering their pivotal roles in regulating various aspects of cell functions and disease progressions. A notable feature of circRNA lies in the mechanism of its biogenesis involving a specialized form of splicing: back-splicing. A splicing process that relies on interactions between introns flanking the circularizing exon to bring the up and downstream splice sites in proximity through the formation of a prerequisite hairpin structure, allowing the spliceosomes to join the two splice sites together to produce a circular RNA molecule. Based on this mechanism, we explored the feasibility of facilitating the formation of such a prerequisite hairpin structure by utilizing a newly designed oligonucleotide, CircuLarIzation Promoting OligoNucleotide (CLIP-ON), to promote the production of circRNA in cells. CLIP-ON was designed to hybridize with and physically bridge two distal sequences in the flanking introns of the circularizing exons. The feasibility of CLIP-ON was confirmed in HeLa cells using a model pre-mRNA, demonstrating the applicability of CLIP-ON as a trans-acting modulator to upregulate the production of circRNAs in a cellular environment.

Bruno 1/CELF regulates splicing and cytoskeleton dynamics to ensure correct sarcomere assembly in Drosophila flight muscles.

Muscles undergo developmental transitions in gene expression and alternative splicing that are necessary to refine sarcomere structure and contractility. CUG-BP and ETR-3-like (CELF) family RNA-binding proteins are important regulators of RNA processing during myogenesis that are misregulated in diseases such as Myotonic Dystrophy Type I (DM1). Here, we report a conserved function for Bruno 1 (Bru1, Arrest), a CELF1/2 family homolog in Drosophila, during early muscle myogenesis. Loss of Bru1 in flight muscles results in disorganization of the actin cytoskeleton leading to aberrant myofiber compaction and defects in pre-myofibril formation. Temporally restricted rescue and RNAi knockdown demonstrate that early cytoskeletal defects interfere with subsequent steps in sarcomere growth and maturation. Early defects are distinct from a later requirement for bru1 to regulate sarcomere assembly dynamics during myofiber maturation. We identify an imbalance in growth in sarcomere length and width during later stages of development as the mechanism driving abnormal radial growth, myofibril fusion, and the formation of hollow myofibrils in bru1 mutant muscle. Molecularly, we characterize a genome-wide transition from immature to mature sarcomere gene isoform expression in flight muscle development that is blocked in bru1 mutants. We further demonstrate that temporally restricted Bru1 rescue can partially alleviate hypercontraction in late pupal and adult stages, but it cannot restore myofiber function or correct structural deficits. Our results reveal the conserved nature of CELF function in regulating cytoskeletal dynamics in muscle development and demonstrate that defective RNA processing due to misexpression of CELF proteins causes wide-reaching structural defects and progressive malfunction of affected muscles that cannot be rescued by late-stage gene replacement.

Transcription elongation defects link oncogenic SF3B1 mutations to targetable alterations in chromatin landscape.

Transcription and splicing of pre-messenger RNA are closely coordinated, but how this functional coupling is disrupted in human diseases remains unexplored. Using isogenic cell lines, patient samples, and a mutant mouse model, we investigated how cancer-associated mutations in SF3B1 alter transcription. We found that these mutations reduce the elongation rate of RNA polymerase II (RNAPII) along gene bodies and its density at promoters. The elongation defect results from disrupted pre-spliceosome assembly due to impaired protein-protein interactions of mutant SF3B1. The decreased promoter-proximal RNAPII density reduces both chromatin accessibility and H3K4me3 marks at promoters. Through an unbiased screen, we identified epigenetic factors in the Sin3/HDAC/H3K4me pathway, which, when modulated, reverse both transcription and chromatin changes. Our findings reveal how splicing factor mutant states behave functionally as epigenetic disorders through impaired transcription-related changes to the chromatin landscape. We also present a rationale for targeting the Sin3/HDAC complex as a therapeutic strategy.

Highly Reactive Group I Introns Ubiquitous in Pathogenic Fungi.

Systemic fungal infections are a growing public health threat, and yet viable antifungal drug targets are limited as fungi share a similar proteome with humans. However, features of RNA metabolism and the noncoding transcriptomes in fungi are distinctive. For example, fungi harbor highly structured RNA elements that humans lack, such as self-splicing introns within key housekeeping genes in the mitochondria. However, the location and function of these mitochondrial riboregulatory elements has largely eluded characterization. Here we used an RNA-structure-based bioinformatics pipeline to identify the group I introns interrupting key mitochondrial genes in medically relevant fungi, revealing their fixation within a handful of genetic hotspots and their ubiquitous presence across divergent phylogenies of fungi, including all highest priority pathogens such as Candida albicans, Candida auris, Aspergillus fumigatus and Cryptococcus neoformans. We then biochemically characterized two representative introns from C. albicans and C. auris, demonstrating their exceptionally efficient splicing catalysis relative to previously-characterized group I introns. Indeed, the C. albicans mitochondrial intron displays extremely rapid catalytic turnover, even at ambient temperatures and physiological magnesium ion concentrations. Our results unmask a significant new set of players in the RNA metabolism of pathogenic fungi, suggesting a promising new type of antifungal drug target.

The impact of genetically controlled splicing on exon inclusion and protein structure.

Common variants affecting mRNA splicing are typically identified though splicing quantitative trait locus (sQTL) mapping and have been shown to be enriched for GWAS signals by a similar degree to eQTLs. However, the specific splicing changes induced by these variants have been difficult to characterize, making it more complicated to analyze the effect size and direction of sQTLs, and to determine downstream splicing effects on protein structure. In this study, we catalogue sQTLs using exon percent spliced in (PSI) scores as a quantitative phenotype. PSI is an interpretable metric for identifying exon skipping events and has some advantages over other methods for quantifying splicing from short read RNA sequencing. In our set of sQTL variants, we find evidence of selective effects based on splicing effect size and effect direction, as well as exon symmetry. Additionally, we utilize AlphaFold2 to predict changes in protein structure associated with sQTLs overlapping GWAS traits, highlighting a potential new use-case for this technology for interpreting genetic effects on traits and disorders.

The spectrum of pre-mRNA splicing in autism.

Disruptions in spatiotemporal gene expression can result in atypical brain function. Specifically, autism spectrum disorder (ASD) is characterized by abnormalities in pre-mRNA splicing. Abnormal splicing patterns have been identified in the brains of individuals with ASD, and mutations in splicing factors have been found to contribute to neurodevelopmental delays associated with ASD. Here we review studies that shed light on the importance of splicing observed in ASD and that explored the intricate relationship between splicing factors and ASD, revealing how disruptions in pre-mRNA splicing may underlie ASD pathogenesis. We provide an overview of the research regarding all splicing factors associated with ASD and place a special emphasis on five specific splicing factors-HNRNPH2, NOVA2, WBP4, SRRM2, and RBFOX1-known to impact the splicing of ASD-related genes. In the discussion of the molecular mechanisms influenced by these splicing factors, we lay the groundwork for a deeper understanding of ASD's complex etiology. Finally, we discuss the potential benefit of unraveling the connection between splicing and ASD for the development of more precise diagnostic tools and targeted therapeutic interventions. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Evolution and Genomics > Computational Analyses of RNA RNA-Based Catalysis > RNA Catalysis in Splicing and Translation.

The exon junction complex is required for DMD gene splicing fidelity and myogenic differentiation.

Deposition of the exon junction complex (EJC) upstream of exon-exon junctions helps maintain transcriptome integrity by preventing spurious re-splicing events in already spliced mRNAs. Here we investigate the importance of EJC for the correct splicing of the 2.2-megabase-long human DMD pre-mRNA, which encodes dystrophin, an essential protein involved in cytoskeletal organization and cell signaling. Using targeted RNA-seq, we show that knock-down of the eIF4A3 and Y14 core components of EJC in a human muscle cell line causes an accumulation of mis-splicing events clustered towards the 3' end of the DMD transcript (Dp427m). This deregulation is conserved in the short Dp71 isoform expressed ubiquitously except in adult skeletal muscle and is rescued with wild-type eIF4A3 and Y14 proteins but not with an EJC assembly-defective mutant eIF4A3. MLN51 protein and EJC-associated ASAP/PSAP complexes independently modulate the inclusion of the regulated exons 71 and 78. Our data confirm the protective role of EJC in maintaining splicing fidelity, which in the DMD gene is necessary to preserve the function of the critical C-terminal protein-protein interaction domain of dystrophin present in all tissue-specific isoforms. Given the role of the EJC in maintaining the integrity of dystrophin, we asked whether the EJC could also be involved in the regulation of a mechanism as complex as skeletal muscle differentiation. We found that eIF4A3 knockdown impairs myogenic differentiation by blocking myotube formation. Collectively, our data provide new insights into the functional roles of EJC in human skeletal muscle.

Comprehensive transcriptome analysis reveals altered mRNA splicing and post-transcriptional changes in the aged mouse brain.

A comprehensive understanding of molecular changes during brain aging is essential to mitigate cognitive decline and delay neurodegenerative diseases. The interpretation of mRNA alterations during brain aging is influenced by the health and age of the animal cohorts studied. Here, we carefully consider these factors and provide an in-depth investigation of mRNA splicing and dynamics in the aging mouse brain, combining short- and long-read sequencing technologies with extensive bioinformatic analyses. Our findings encompass a spectrum of age-related changes, including differences in isoform usage, decreased mRNA dynamics and a module showing increased expression of neuronal genes. Notably, our results indicate a reduced abundance of mRNA isoforms leading to nonsense-mediated RNA decay and suggest a regulatory role for RNA-binding proteins, indicating that their regulation may be altered leading to the reshaping of the aged brain transcriptome. Collectively, our study highlights the importance of studying mRNA splicing events during brain aging.

Exitron splicing of odor receptor genes in Drosophila.

Proper expression of odor receptor genes is critical for the function of olfactory systems. In this study, we identified exitrons (exonic introns) in four of the 39 Odorant receptor (Or) genes expressed in the Drosophila antenna. Exitrons are sequences that can be spliced out from within a protein-coding exon, thereby altering the encoded protein. We focused on Or88a, which encodes a pheromone receptor, and found that exitron splicing of Or88a is conserved across five Drosophila species over 20 My of evolution. The exitron was spliced out in 15% of Or88a transcripts. Removal of this exitron creates a non-coding RNA rather than an RNA that encodes a stable protein. Our results suggest the hypothesis that in the case of Or88a, exitron splicing could act in neuronal modulation by decreasing the level of functional Or transcripts. Activation of Or88a-expressing olfactory receptor neurons via either optogenetics or pheromone stimulation increased the level of exitron-spliced transcripts, with optogenetic activation leading to a 14-fold increase. A fifth Or can also undergo an alternative splicing event that eliminates most of the canonical open reading frame. Besides these cases of alternative splicing, we found alternative polyadenylation of four Ors, and exposure of Or67c to its ligand ethyl lactate in the antenna downregulated all of its 3' isoforms. Our study reveals mechanisms by which neuronal activity could be modulated via regulation of the levels of Or isoforms.

Novel Splicing Variants in the ARR3 Gene Cause the Female-Limited Early-Onset High Myopia.

Purpose: Variants in the ARR3 gene have been linked to early-onset high myopia (eoHM) with a unique X-linked female-limited inheritance. However, the clinical validity of this gene-disease association has not been systematically evaluated. Methods: We identified two Chinese families with novel ARR3 splicing variants associated with eoHM. Minigene constructs were generated to assess the effects of the variants on splicing. We integrated previous evidence to curate the clinical validity of ARR3 and eoHM using the ClinGen framework. Results: The variants c.39+1G>A and c.100+4A>G were identified in the two families. Minigene analysis showed both variants resulted in abnormal splicing and introduction of premature termination codons. Based on genetic and experimental evidence, the ARR3-eoHM relationship was classified as "definitive." Conclusions: Our study identified two novel splicing variants of the ARR3 gene linked to eoHM and confirmed their functional validity via minigene assay. This research expanded the mutational spectrum of ARR3 and confirmed the minigene assay technique as an effective tool for understanding variant effects on splicing mechanisms.

RNA therapeutics for ß-thalassemia.

ß-thalassemia is an autosomal recessive disease, caused by one or more mutations in the ß-globin gene that reduces or abolishes ß-globin chain synthesis causing an imbalance in the ratio of α- and ß-globin chain. Therefore, the ability to target mutations will provide a good result in the treatment of ß-thalassemia. RNA therapeutics represents a promising class of drugs inclusive antisense oligonucleotides (ASO), small interfering RNA (siRNA), microRNA (miRNA) and APTAMER have investigated in clinical trials for treatment of human diseases as ß-thalassemia; Especially, ASO therapeutics can completely treat ß-thalassemia patients by the way of making ASO infiltrating through erythrocyte progenitor cells, migrating to the nucleus and hybridizing with abnormal splicing sites to suppress an abnormal splicing pattern of ß-globin pre-mRNA. As a result, the exactly splicing process is restored to increase the expression of ß-globin which increases the amount of mature hemoglobin of red blood cells of ß-thalassemia patients. Furthermore, current study demonstrates that RNA-based therapeutics get lots of good results for ß-thalassemia patients. Then, this chapter focuses on current advances of RNA-based therapeutics and addresses current challenges with their development and application for treatment of ß-thalassemia patients.

A common cellular response to broad splicing perturbations is characterized by metabolic transcript downregulation driven by the Mdm2-p53 axis.

Disruptions in core cellular processes elicit stress responses that drive cell-state changes leading to organismal phenotypes. Perturbations in the splicing machinery cause widespread mis-splicing, resulting in p53-dependent cell-state changes that give rise to cell-type-specific phenotypes and disease. However, a unified framework for how cells respond to splicing perturbations, and how this response manifests itself in nuanced disease phenotypes, has yet to be established. Here, we show that a p53-stabilizing Mdm2 alternative splicing event and the resulting widespread downregulation of metabolic transcripts are common events that arise in response to various splicing perturbations in both cellular and organismal models. Together, our results classify a common cellular response to splicing perturbations, put forth a new mechanism behind the cell-type-specific phenotypes that arise when splicing is broadly disrupted, and lend insight into the pleiotropic nature of the effects of p53 stabilization in disease.

Interrogations of single-cell RNA splicing landscapes with SCASL define new cell identities with physiological relevance.

RNA splicing shapes the gene regulatory programs that underlie various physiological and disease processes. Here, we present the SCASL (single-cell clustering based on alternative splicing landscapes) method for interrogating the heterogeneity of RNA splicing with single-cell RNA-seq data. SCASL resolves the issue of biased and sparse data coverage on single-cell RNA splicing and provides a new scheme for classifications of cell identities. With previously published datasets as examples, SCASL identifies new cell clusters indicating potentially precancerous and early-tumor stages in triple-negative breast cancer, illustrates cell lineages of embryonic liver development, and provides fine clusters of highly heterogeneous tumor-associated CD4 and CD8 T cells with functional and physiological relevance. Most of these findings are not readily available via conventional cell clustering based on single-cell gene expression data. Our study shows the potential of SCASL in revealing the intrinsic RNA splicing heterogeneity and generating biological insights into the dynamic and functional cell landscapes in complex tissues.

Detecting and understanding meaningful cancerous mutations based on computational models of mRNA splicing.

Cancer research has long relied on non-silent mutations. Yet, it has become overwhelmingly clear that silent mutations can affect gene expression and cancer cell fitness. One fundamental mechanism that apparently silent mutations can severely disrupt is alternative splicing. Here we introduce Oncosplice, a tool that scores mutations based on models of proteomes generated using aberrant splicing predictions. Oncosplice leverages a highly accurate neural network that predicts splice sites within arbitrary mRNA sequences, a greedy transcript constructor that considers alternate arrangements of splicing blueprints, and an algorithm that grades the functional divergence between proteins based on evolutionary conservation. By applying this tool to 12M somatic mutations we identify 8K deleterious variants that are significantly depleted within the healthy population; we demonstrate the tool's ability to identify clinically validated pathogenic variants with a positive predictive value of 94%; we show strong enrichment of predicted deleterious mutations across pan-cancer drivers. We also achieve improved patient survival estimation using a proposed set of novel cancer-involved genes. Ultimately, this pipeline enables accelerated insight-gathering of sequence-specific consequences for a class of understudied mutations and provides an efficient way of filtering through massive variant datasets - functionalities with immediate experimental and clinical applications.