Post-transcriptional splicing can occur in a slow-moving zone around the gene.

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

Self-supervised learning on millions of primary RNA sequences from 72 vertebrates improves sequence-based RNA splicing prediction.

Language models pretrained by self-supervised learning (SSL) have been widely utilized to study protein sequences, while few models were developed for genomic sequences and were limited to single species. Due to the lack of genomes from different species, these models cannot effectively leverage evolutionary information. In this study, we have developed SpliceBERT, a language model pretrained on primary ribonucleic acids (RNA) sequences from 72 vertebrates by masked language modeling, and applied it to sequence-based modeling of RNA splicing. Pretraining SpliceBERT on diverse species enables effective identification of evolutionarily conserved elements. Meanwhile, the learned hidden states and attention weights can characterize the biological properties of splice sites. As a result, SpliceBERT was shown effective on several downstream tasks: zero-shot prediction of variant effects on splicing, prediction of branchpoints in humans, and cross-species prediction of splice sites. Our study highlighted the importance of pretraining genomic language models on a diverse range of species and suggested that SSL is a promising approach to enhance our understanding of the regulatory logic underlying genomic sequences.

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.

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.

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.

RNA-binding proteins in pain.

RNAs are meticulously controlled by proteins. Through direct and indirect associations, every facet in the brief life of an mRNA is subject to regulation. RNA-binding proteins (RBPs) permeate biology. Here, we focus on their roles in pain. Chronic pain is among the largest challenges facing medicine and requires new strategies. Mounting pharmacologic and genetic evidence obtained in pre-clinical models suggests fundamental roles for a broad array of RBPs. We describe their diverse roles that span RNA modification, splicing, stability, translation, and decay. Finally, we highlight opportunities to expand our understanding of regulatory interactions that contribute to pain signaling. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Regulation RNA in Disease and Development > RNA in Disease.

On a kneading theory for gene-splicing.

Two well-known facets in protein synthesis in eukaryotic cells are transcription of DNA to pre-RNA in the nucleus and the translation of messenger-RNA (mRNA) to proteins in the cytoplasm. A critical intermediate step is the removal of segments (introns) containing ∼97% of the nucleic-acid sites in pre-RNA and sequential alignment of the retained segments (exons) to form mRNA through a process referred to as splicing. Alternative forms of splicing enrich the proteome while abnormal splicing can enhance the likelihood of a cell developing cancer or other diseases. Mechanisms for splicing and origins of splicing errors are only partially deciphered. Our goal is to determine if rules on splicing can be inferred from data analytics on nucleic-acid sequences. Toward that end, we represent a nucleic-acid site as a point in a plane defined in terms of the anterior and posterior sub-sequences of the site. The "point-set" representation expands analytical approaches, including the use of statistical tools, to characterize genome sequences. It is found that point-sets for exons and introns are visually different, and that the differences can be quantified using a family of generalized moments. We design a machine-learning algorithm that can recognize individual exons or introns with 91% accuracy. Point-set distributions and generalized moments are found to differ between organisms.

Structural basis of U12-type intron engagement by the fully assembled human minor spliceosome.

The minor spliceosome, which is responsible for the splicing of U12-type introns, comprises five small nuclear RNAs (snRNAs), of which only one is shared with the major spliceosome. In this work, we report the 3.3-angstrom cryo-electron microscopy structure of the fully assembled human minor spliceosome pre-B complex. The atomic model includes U11 small nuclear ribonucleoprotein (snRNP), U12 snRNP, and U4atac/U6atac.U5 tri-snRNP. U11 snRNA is recognized by five U11-specific proteins (20K, 25K, 35K, 48K, and 59K) and the heptameric Sm ring. The 3' half of the 5'-splice site forms a duplex with U11 snRNA; the 5' half is recognized by U11-35K, U11-48K, and U11 snRNA. Two proteins, CENATAC and DIM2/TXNL4B, specifically associate with the minor tri-snRNP. A structural analysis uncovered how two conformationally similar tri-snRNPs are differentiated by the minor and major prespliceosomes for assembly.

Structural and mechanistic insights into activation of the human RNA ligase RTCB by Archease.

RNA ligases of the RTCB-type play an essential role in tRNA splicing, the unfolded protein response and RNA repair. RTCB is the catalytic subunit of the pentameric human tRNA ligase complex. RNA ligation by the tRNA ligase complex requires GTP-dependent activation of RTCB. This active site guanylylation reaction relies on the activation factor Archease. The mechanistic interplay between both proteins has remained unknown. Here, we report a biochemical and structural analysis of the human RTCB-Archease complex in the pre- and post-activation state. Archease reaches into the active site of RTCB and promotes the formation of a covalent RTCB-GMP intermediate through coordination of GTP and metal ions. During the activation reaction, Archease prevents futile RNA substrate binding to RTCB. Moreover, monomer structures of Archease and RTCB reveal additional states within the RNA ligation mechanism. Taken together, we present structural snapshots along the reaction cycle of the human tRNA ligase.

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 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.

Successful skipping of abnormal pseudoexon by antisense oligonucleotides in vitro for a patient with beta-propeller protein-associated neurodegeneration.

Pathogenic variants in WDR45 on chromosome Xp11 cause neurodegenerative disorder beta-propeller protein-associated neurodegeneration (BPAN). Currently, there is no effective therapy for BPAN. Here we report a 17-year-old female patient with BPAN and show that antisense oligonucleotide (ASO) was effective in vitro. The patient had developmental delay and later showed extrapyramidal signs since the age of 15 years. MRI findings showed iron deposition in the globus pallidus and substantia nigra on T2 MRI. Whole genome sequencing and RNA sequencing revealed generation of pseudoexon due to inclusion of intronic sequences triggered by an intronic variant that is remote from the exon-intron junction: WDR45 (OMIM #300526) chrX(GRCh37):g.48935143G > C, (NM_007075.4:c.235 + 159C > G). We recapitulated the exonization of intron sequences by a mini-gene assay and further sought antisense oligonucleotide that induce pseudoexon skipping using our recently developed, a dual fluorescent splicing reporter system that encodes two fluorescent proteins, mCherry, a transfection marker designed to facilitate evaluation of exon skipping and split eGFP, a splicing reaction marker. The results showed that the 24-base ASO was the strongest inducer of pseudoexon skipping. Our data presented here have provided supportive evidence for in vivo preclinical studies.

RNA aptamer reveals nuclear TDP-43 pathology is an early aggregation event that coincides with STMN-2 cryptic splicing and precedes clinical manifestation in ALS.

TDP-43 is an aggregation-prone protein which accumulates in the hallmark pathological inclusions of amyotrophic lateral sclerosis (ALS). However, the analysis of deeply phenotyped human post-mortem samples has shown that TDP-43 aggregation, revealed by standard antibody methods, correlates poorly with symptom manifestation. Recent identification of cryptic-splicing events, such as the detection of Stathmin-2 (STMN-2) cryptic exons, are providing evidence implicating TDP-43 loss-of-function as a potential driving pathomechanism but the temporal nature of TDP-43 loss and its relation to the disease process and clinical phenotype is not known. To address these outstanding questions, we used a novel RNA aptamer, TDP-43APT, to detect TDP-43 pathology and used single molecule in situ hybridization to sensitively reveal TDP-43 loss-of-function and applied these in a deeply phenotyped human post-mortem tissue cohort. We demonstrate that TDP-43APT identifies pathological TDP-43, detecting aggregation events that cannot be detected by classical antibody stains. We show that nuclear TDP-43 pathology is an early event, occurring prior to cytoplasmic accumulation and is associated with loss-of-function measured by coincident STMN-2 cryptic splicing pathology. Crucially, we show that these pathological features of TDP-43 loss-of-function precede the clinical inflection point and are not required for region specific clinical manifestation. Furthermore, we demonstrate that gain-of-function in the form of extensive cytoplasmic accumulation, but not loss-of-function, is the primary molecular correlate of clinical manifestation. Taken together, our findings demonstrate implications for early diagnostics as the presence of STMN-2 cryptic exons and early TDP-43 aggregation events could be detected prior to symptom onset, holding promise for early intervention in ALS.

Detecting and Validating MAPT Mutations in Neurodegeneration Patients and Analysis of Exon Splicing Consequences.

Mutation of MAPT has been observed in patients with parkinsonism, progressive supranuclear palsy, and corticobasal degeneration and is a significant cause of frontotemporal dementia. In this chapter, we discuss considerations for next-generation sequencing analysis to identify MAPT mutations in patient genomic DNA and describe the validation of these mutations by Sanger sequencing. One of the most common effects of MAPT mutations is differential splicing of exon 10, which leads to an imbalance in the proportion of 3-repeat and 4-repeat tau isoforms. We describe how to investigate the effect of novel DNA variants on the splicing efficiency of this exon in vitro using the exon-trapping technique, also known as the splicing reporter minigene assay.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans.

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that 'complex' organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this 'drift barrier' theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

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.