Human branch point consensus sequence is yUnAy

Invention of an oral medication for cardiac Fabry disease caused by RNA mis-splicing

Pathogenic RNA splicing variants have emerged as promising therapeutic targets due to their role in disease while preserving coding sequences. In this study, we developed RECTAS-2.0, a small molecule designed to correct RNA mis-splicing caused by the GLA c.639+919G>A mutation, which leads to the inclusion of a 57-nucleotide poison exon, resulting in later-onset Fabry disease, particularly prevalent in East Asia. RECTAS-2.0 restored normal GLA mRNA splicing and α-galactosidase activity in patient-derived B-lymphoblastoid cell lines and induced pluripotent stem cell-derived cardiomyocytes. Furthermore, oral administration of RECTAS-2.0 effectively corrected splicing in a transgenic mouse model, demonstrating its substantial splice-switching activity and safety for clinical application. RECTAS-2.0 demonstrated potential applicability to other genetic disorders that involve similar exon competition. These findings underscore the therapeutic potential of RECTAS-2.0 for Fabry disease and highlight its broader implications for RNA splicing-targeted therapies in genetic disorders.

Intron-Derived Lariat RNAs Go Stable

During pre-mRNA splicing, introns are removed by the spliceosome, and the flanking exons are ligated to form mature mRNA, which is subsequently translated into protein. Traditionally, intronic RNAs have been regarded as "junk", presumed to be degraded for nucleotide turnover. Notably, after debranching, some linearized lariat RNAs can be further processed into snoRNAs, miRNAs, and other long non-coding RNAs. However, recent studies have shown that many intron-derived lariat RNAs can escape degradation and remain stable across various eukaryotic organisms, indicating they may play significant roles in cellular processes. Moreover, these naturally retained lariat RNAs are frequently observed in circular forms in vivo, suggesting that their linear tails are highly susceptible to degradation. This highlights lariat RNAs as an important source of circular RNAs. Furthermore, many lariat-derived circRNAs have been detected in the cytoplasm, implying active nuclear export and potential roles in cytoplasmic processes. In this review, we provide an overview of the life cycle of intron-derived lariat RNAs, focusing on their biogenesis, degradation, and retention. We also discuss the mechanisms that enable their resistance to degradation and the biological functions of stable lariat RNAs, shedding light on these seemingly "nonsense" yet inevitably produced non-coding intronic RNAs.

Salvage pathway of vitamin B12 absorption in chickens with mutant tumor virus a receptor

The tumor virus A receptor (TVA), a member of the low-density lipoprotein receptor (LDLR) family, serves as an entry receptor for Avian Leukosis Virus (ALV) subgroups A and K, as well as a receptor for vitamin B12 bound to transcobalamin. Naturally occurring genetic variants in the TVA gene determine susceptibility or resistance to ALV-A and -K, but the effects of these mutated TVA on vitamin B12 uptake have not been investigated systemically. We found four TVA variants comprising the wild type (TVAWT), a single nucleotide polymorphism variant (TVASNP), and two partial deletions in the splicing branch point region (TVAR). This study investigates the relationship between the various genotypes of TVA alleles and uptake of vitamin B12 in chickens. A protein interaction model suggested that mutant TVAs (i.e., TVASNP, TVAR) may have reduced ability to take up vitamin B12 due to a disrupted LDL-A domain, a pivotal region involved in vitamin B12 uptake; however, we found no significant difference in absorption of vitamin B12 in TVAWT and TVASNP chickens, or levels of its metabolite in serum. Notably, TVAR chickens had significantly higher levels of vitamin B12 than TVAWT chickens, a finding contrary to the predicted lower uptake. Expression of vitamin B12 carrier related genes (i.e., AMN, GIF, and TCN2) in chickens showed a stepwise increase: TVAWT > TVASNP > TVAR. These results suggest a mechanism by which mutant TVA chickens with a disrupted TVA protein acquire natural resistance to ALV-A -K, with no impairment of vitamin B12 metabolism.

Translation of the downstream ORF from bicistronic mRNAs by human cells: Impact of codon usage and splicing in the upstream ORF

Biochemistry textbooks describe eukaryotic mRNAs as monocistronic. However, increasing evidence reveals the widespread presence and translation of upstream open reading frames preceding the "main" ORF. DNA and RNA viruses infecting eukaryotes often produce polycistronic mRNAs and viruses have evolved multiple ways of manipulating the host's translation machinery. Here, we introduce an experimental model to study gene expression regulation from virus-like bicistronic mRNAs in human cells. The model consists of a short upstream ORF and a reporter downstream ORF encoding a fluorescent protein. We have engineered synonymous variants of the upstream ORF to explore large parameter space, including codon usage preferences, mRNA folding features, and splicing propensity. We show that human translation machinery can translate the downstream ORF from bicistronic mRNAs, albeit reporter protein levels are thousand times lower than those from the upstream ORF. Furthermore, synonymous recoding of the upstream ORF exclusively during elongation significantly influences its own translation efficiency, reveals cryptic splice signals, and modulates the probability of downstream ORF translation. Our results are consistent with a leaky scanning mechanism facilitating downstream ORF translation from bicistronic mRNAs in human cells, offering new insights into the role of upstream ORFs in translation regulation.

Hyperactivation of NF-κB signaling in splicing factor mutant myelodysplastic syndromes and therapeutic approaches

The transcription factor NF-κB plays a critical role in the control of innate and adaptive immunity and inflammation. Several recent studies have demonstrated that the mutation of different splicing factor genes, including SF3B1, SRSF2 and U2AF1, in myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML) result in hyperactive NF-κB signaling through the aberrant splicing of different target genes. The presence of U2AF1 and SF3B1 mutations in the bone marrow cells of MDS and AML patients induces oncogenic isoforms of the target gene IRAK4, leading to hyperactivation of NF-κB signaling and an increase in the fitness of leukemic stem and progenitor cells (LSPCs). The potent IRAK4 inhibitor CA-4948 has shown efficacy in both pre-clinical studies and MDS clinical trials, with splicing factor mutant patients showing the higher response rates. Emerging data has, however, revealed that co-targeting of IRAK4 and its paralog IRAK1 is required to maximally suppress LSPC function in vitro and in vivo by inducing cellular differentiation. These findings provide a link between the presence of the commonly mutated splicing factor genes and activation of innate immune signaling pathways in myeloid malignancies and have important implications for targeted therapy in these disorders.

Two novel deep intronic variants cause Duchenne muscular dystrophy by splice-altering mechanism

Duchenne muscular dystrophy (DMD) is a genetic disorder characterized by progressive muscle degeneration and weakness, due to mutations in the DMD gene, which encodes the dystrophin protein. While mutations within the coding regions of DMD have been extensively studied, recent focus has shifted to deep intronic variants for their potential impact on disease severity. Here, we characterize two deep intronic variants, c.8669-19_8669-24del and c.6439-1016_6439-3376del, in unrelated DMD patients. These variants were identified using targeted long-read sequencing on patients' DNA. RNA sequencing/reverse transcription polymerase chain reaction on RNA extracted from muscle biopsies revealed the presence of a pseudoexon or retention of part of the intron in the transcript, resulting in the introduction of premature termination codons. This study enhances our understanding of pseudoexon activation mechanisms in DMD and underscores the diverse genetic abnormalities contributing to the disease's complexity.

Somatic Recombination Between an Ancient and a Recent NOTCH2 Gene Variant Is Associated with the NOTCH2 Gain-of-Function Phenotype in Chronic Lymphocytic Leukemia

Constitutively active NOTCH2 signaling is a hallmark in chronic lymphocytic leukemia (CLL). The precise underlying defect remains obscure. Here we show that the mRNA sequence coding for the NOTCH2 negative regulatory region (NRR) is consistently deleted in CLL cells. The most common NOTCH2ΔNRR-DEL2 deletion is associated with two intronic single nucleotide variations (SNVs) which either create (CTTAT, G>A for rs2453058) or destroy (CTCGT, A>G for rs5025718) a putative splicing branch point sequence (BPS). Phylogenetic analysis demonstrates that rs2453058 is part of an ancient NOTCH2 gene variant (*1A01) which is associated with type 2 diabetes mellitus (T2DM) and is two times more frequent in Europeans than in East Asians, resembling the differences in CLL incidence. In contrast, rs5025718 belongs to a recent NOTCH2 variant (*1a4) that dominates the world outside Africa. Nanopore sequencing indicates that somatic reciprocal crossing over between rs2453058 (*1A01) and rs5025718 (*1a4) leads to recombined NOTCH2 alleles with altered BPS patterns in NOTCH2*1A01/*1a4 CLL cases. This would explain the loss of the NRR domain by aberrant pre-mRNA splicing and consequently the NOTCH2 gain-of-function phenotype. Together, our findings suggest that somatic recombination of inherited NOTCH2 variants might be relevant to CLL etiology and may at least partly explain its geographical clustering.

Genetic excision of the regulatory cardiac troponin I extension in high-heart rate mammal clades

Mammalian cardiac troponin I (cTnI) contains a highly conserved amino-terminal extension harboring protein kinase A targets [serine-23 and -24 (Ser23/24)] that are phosphorylated during β-adrenergic stimulation to defend diastolic filling by means of an increased cardiomyocyte relaxation rate. In this work, we show that the Ser23/24-encoding exon 3 of TNNI3 was pseudoexonized multiple times in shrews and moles to mimic Ser23/24 phosphorylation without adrenergic stimulation, facilitating the evolution of exceptionally high resting heart rates (~1000 beats per minute). We further reveal alternative exon 3 splicing in distantly related bat families and confirm that both cTnI splice variants are incorporated into cardiac myofibrils. Because exon 3 of human TNNI3 exhibits a relatively low splice strength score, our findings offer an evolutionarily informed strategy to excise this exon to improve diastolic function during heart failure.

Intronic RNA secondary structural information captured for the human MYC pre-mRNA

To address the lack of intronic reads in secondary structure probing data for the human MYC pre-mRNA, we developed a method that combines spliceosomal inhibition with RNA probing and sequencing. Here, the SIRP-seq method was applied to study the secondary structure of human MYC RNAs by chemically probing HeLa cells with dimethyl sulfate in the presence of the small molecule spliceosome inhibitor pladienolide B. Pladienolide B binds to the SF3B complex of the spliceosome to inhibit intron removal during splicing, resulting in retained intronic sequences. This method was used to increase the read coverage over intronic regions of MYC. The purpose for increasing coverage across introns was to generate complete reactivity profiles for intronic sequences via the DMS-MaPseq approach. Notably, depth was sufficient for analysis by the program DRACO, which was able to deduce distinct reactivity profiles and predict multiple secondary structural conformations as well as their suggested stoichiometric abundances. The results presented here provide a new method for intronic RNA secondary structural analyses, as well as specific structural insights relevant to MYC RNA splicing regulation and therapeutic targeting.

Human cells contain myriad excised linear intron RNAs with links to gene regulation and potential utility as biomarkers

A previous study using Thermostable Group II Intron Reverse Transcriptase sequencing (TGIRT-seq) found human plasma contains short (≤300 nt) structured full-length excised linear intron (FLEXI) RNAs with potential to serve as blood-based biomarkers. Here, TGIRT-seq identified >9,000 different FLEXI RNAs in human cell lines, including relatively abundant FLEXIs with cell-type-specific expression patterns. Analysis of public CLIP-seq datasets identified 126 RNA-binding proteins (RBPs) that have binding sites within the region corresponding to the FLEXI or overlapping FLEXI splice sites in pre-mRNAs, including 53 RBPs with binding sites for ≥30 different FLEXIs. These included splicing factors, transcription factors, a chromatin remodeling protein, cellular growth regulators, and proteins with cytoplasmic functions. Analysis of ENCODE datasets identified subsets of these RBPs whose knockdown impacted FLEXI host gene mRNA levels or proximate alternative splicing, indicating functional interactions. Hierarchical clustering identified six subsets of RBPs whose FLEXI binding sites were co-enriched in six subsets of functionally related host genes: AGO1-4 and DICER, including but not limited to agotrons or mirtron pre-miRNAs; DKC1, NOLC1, SMNDC1, and AATF (Apoptosis Antagonizing Transcription Factor), including but not limited to snoRNA-encoding FLEXIs; two subsets of alternative splicing factors; and two subsets that included RBPs with cytoplasmic functions (e.g., LARP4, PABPC4, METAP2, and ZNF622) together with regulatory proteins. Cell fractionation experiments showed cytoplasmic enrichment of FLEXI RNAs with binding sites for RBPs with cytoplasmic functions. The subsets of host genes encoding FLEXIs with binding sites for different subsets of RBPs were co-enriched with non-FLEXI other short and long introns with binding sites for the same RBPs, suggesting overarching mechanisms for coordinately regulating expression of functionally related genes. Our findings identify FLEXIs as a previously unrecognized large class of cellular RNAs and provide a comprehensive roadmap for further analyzing their biological functions and the relationship of their RBPs to cellular regulatory mechanisms.

Taxonomy of introns and the evolution of minor introns

Classification of introns, which is crucial to understanding their evolution and splicing, has historically been binary and has resulted in the naming of major and minor introns that are spliced by their namesake spliceosome. However, a broad range of intron consensus sequences exist, leading us to here reclassify introns as minor, minor-like, hybrid, major-like, major and non-canonical introns in 263 species across six eukaryotic supergroups. Through intron orthology analysis, we discovered that minor-like introns are a transitory node for intron conversion across evolution. Despite close resemblance of their consensus sequences to minor introns, these introns possess an AG dinucleotide at the -1 and -2 position of the 5' splice site, a salient feature of major introns. Through combined analysis of CoLa-seq, CLIP-seq for major and minor spliceosome components, and RNAseq from samples in which the minor spliceosome is inhibited we found that minor-like introns are also an intermediate class from a splicing mechanism perspective. Importantly, this analysis has provided insight into the sequence elements that have evolved to make minor-like introns amenable to recognition by both minor and major spliceosome components. We hope that this revised intron classification provides a new framework to study intron evolution and splicing.

Refining clinically relevant parameters for mis-splicing risk in shortened introns with donor-to-branchpoint space constraint

Intronic deletions that critically shorten donor-to-branchpoint (D-BP) distance of a precursor mRNA impose biophysical space constraint on assembly of the U1/U2 spliceosomal complex, leading to canonical splicing failure. Here we use a series of β-globin (HBB) gene constructs with intron 1 deletions to define D-BP lengths that present low/no risk of mis-splicing and lengths which are critically short and likely elicit clinically relevant mis-splicing. We extend our previous observation in EMD intron 5 of 46 nt as the minimal productive D-BP length, demonstrating spliceosome assembly constraint persists at D-BP lengths of 47-56 nt. We exploit the common HBB exon 1 β-thalassemia variant that strengthens a cryptic donor (NM_000518.5(HBB):c.79G > A) to provide a simple barometer for the earliest signs of space constraint, via cryptic donor activation. For clinical evaluation of intronic deletions, we assert D-BP lengths > 60 nt present low mis-splicing risk while space constraint increases exponentially with D-BP lengths < 55 nt, with critical risk and profound splicing abnormalities with D-BP lengths < 50 nt.

Conserved intronic secondary structures with concealed branch sites regulate alternative splicing of poison exons

Alternative splicing (AS) generates multiple RNA isoforms and increases the complexities of transcriptomes and proteomes. However, it remains unclear how RNA structures contribute to AS regulation. Here, we systematically search transcriptomes for secondary structures with concealed branch sites (BSs) in the alternatively spliced introns and predict thousands of them from six organisms, of which many are evolutionarily conserved. Intriguingly, a highly conserved stem-loop structure with concealed BSs is found in animal SF3B3 genes and colocalizes with a downstream poison exon (PE). Destabilization of this structure allows increased usage of the BSs and results in enhanced PE inclusion in human and Drosophila cells, leading to decreased expression of SF3B3. This structure is experimentally validated using an in-cell SHAPE-MaP assay. Through RNA interference screens of 28 RNA-binding proteins, we find that this stem-loop structure is sensitive to U2 factors. Furthermore, we find that SF3B3 also facilitates DNA repair and protects genome stability by enhancing interaction between ERCC6/CSB and arrested RNA polymerase II. Importantly, both Drosophila and human cells with the secondary structure mutated by genome editing exhibit altered DNA repair in vivo. This study provides a novel and common mechanism for AS regulation of PEs and reveals a physiological function of SF3B3 in DNA repair.

CaTCH: Calculating transcript complexity of human genes

The findings based on whole transcriptome sequencing suggest that alternative splicing occurs in approximately 95% of human multi-exon genes, thus, playing a crucial role in promoting proteome diversity. According to the latest GENCODE annotations, most genes have less than four transcripts, positively correlating with the number of exons. Thus, it is more accurate to measure the splice variant efficiency of a gene with respect to the number of exons, which is a measure of Transcript Complexity (TC). In addition to that, the theoretical number of transcripts is substantially higher than the actual number of transcripts produced by Alternative Splicing Events, and the features restricting this phenomenon need to be explored. In this method, we have extracted the data of various features contributing to TC from different databases. Linear regression is used to identify the determinant features and to train and test the model of TC. The results indicate that exon length is the determining feature of TC, followed by coding potential, presence of chromatin signature, and 5' splice site dinucleotide, all of which negatively affect a gene's TC, except exon length. To further classify the genes based on TC, random forest is used to identify the determinant features.•The splicing efficiency of a gene can be inferred by the transcript complexity, which is the number of transcripts per exon.•CaTCH is a linear regression-based model to calculate the transcript complexity of human genes, which can be calculated from the exon length, coding potentiality, presence of chromatin signature/s, and 5' splice site dinucleotide.

A systematic screen identifies Saf5 as a link between splicing and transcription in fission yeast

Splicing is an important step of gene expression regulation in eukaryotes, as there are many mRNA precursors that can be alternatively spliced in different tissues, at different cell cycle phases or under different external stimuli. We have developed several integrated fluorescence-based in vivo splicing reporter constructs that allow the quantification of fission yeast splicing in vivo on intact cells, and we have compared their splicing efficiency in a wild type strain and in a prp2-1 (U2AF65) genetic background, showing a clear dependency between Prp2 and a consensus signal at 5' splicing site (5'SS). To isolate novel genes involved in regulated splicing, we have crossed the reporter showing more intron retention with the Schizosaccharomyces pombe knock out collection. Among the candidate genes involved in the regulation of splicing, we have detected strong splicing defects in two of the mutants -Δcwf12, a member of the NineTeen Complex (NTC) and Δsaf5, a methylosome subunit that acts together with the survival motor neuron (SMN) complex in small nuclear ribonucleoproteins (snRNP) biogenesis. We have identified that strains with mutations in cwf12 have inefficient splicing, mainly when the 5'SS differs from the consensus. However, although Δsaf5 cells also have some dependency on 5'SS sequence, we noticed that when one intron of a given pre-mRNA was affected, the rest of the introns of the same pre-mRNA had high probabilities of being also affected. This observation points Saf5 as a link between transcription rate and splicing.

GPATCH8 modulates mutant SF3B1 mis-splicing and pathogenicity in hematologic malignancies

Mutations in the RNA splicing factor gene SF3B1 are common across hematologic and solid cancers and result in widespread alterations in splicing, yet there is currently no therapeutic means to correct this mis-splicing. Here, we utilize synthetic introns uniquely responsive to mutant SF3B1 to identify trans factors required for aberrant mutant SF3B1 splicing activity. This revealed the G-patch domain-containing protein GPATCH8 as required for mutant SF3B1-induced splicing alterations and impaired hematopoiesis. GPATCH8 is involved in quality control of branchpoint selection, interacts with the RNA helicase DHX15, and functionally opposes SURP and G-patch domain containing 1 (SUGP1), a G-patch protein recently implicated in SF3B1-mutant diseases. Silencing of GPATCH8 corrected one-third of mutant SF3B1-dependent splicing defects and was sufficient to improve dysfunctional hematopoiesis in SF3B1-mutant mice and primary human progenitors. These data identify GPATCH8 as a novel splicing factor required for mis-splicing by mutant SF3B1 and highlight the therapeutic impact of correcting aberrant splicing in SF3B1-mutant cancers.

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.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Programmable RNA writing with trans-splicing

RNA editing offers the opportunity to introduce either stable or transient modifications to nucleic acid sequence without permanent off-target effects, but installation of arbitrary edits into the transcriptome is currently infeasible. Here, we describe Programmable RNA Editing & Cleavage for Insertion, Substitution, and Erasure (PRECISE), a versatile RNA editing method for writing RNA of arbitrary length and sequence into existing pre-mRNAs via 5' or 3' trans-splicing. In trans-splicing, an exogenous template is introduced to compete with the endogenous pre-mRNA, allowing for replacement of upstream or downstream exon sequence. Using Cas7-11 cleavage of pre-mRNAs to bias towards editing outcomes, we boost the efficiency of RNA trans-splicing by 10-100 fold, achieving editing rates between 5-50% and 85% on endogenous and reporter transcripts, respectively, while maintaining high-fidelity. We demonstrate PRECISE editing across 11 distinct endogenous transcripts of widely varying expression levels, showcasing more than 50 types of edits, including all 12 possible transversions and transitions, insertions ranging from 1 to 1,863 nucleotides, and deletions. We show high efficiency replacement of exon 4 of MECP2, addressing most mutations that drive the Rett Syndrome; editing of SHANK3 transcripts, a gene involved in Autism; and replacement of exon 1 of HTT, removing the hallmark repeat expansions of Huntington's disease. Whole transcriptome sequencing reveals the high precision of PRECISE editing and lack of off-target trans-splicing activity. Furthermore, we combine payload engineering and ribozymes for protein-free, high-efficiency trans-splicing, with demonstrated efficiency in editing HTT exon 1 via AAV delivery. We show that the high activity of PRECISE editing enables editing in non-dividing neurons and patient-derived Huntington's disease fibroblasts. PRECISE editing markedly broadens the scope of genetic editing, is straightforward to deliver over existing gene editing tools like prime editing, lacks permanent off-targets, and can enable any type of genetic edit large or small, including edits not otherwise possible with existing RNA base editors, widening the spectrum of addressable diseases.

Structural basis of branching during RNA splicing

Branching is a critical step in RNA splicing that is essential for 5' splice site selection. Recent spliceosome structures have led to competing models for the recognition of the invariant adenosine at the branch point. However, there are no structures of any splicing complex with the adenosine nucleophile docked in the active site and positioned to attack the 5' splice site. Thus we lack a mechanistic understanding of adenosine selection and splice site recognition during RNA splicing. Here we present a cryo-electron microscopy structure of a group II intron that reveals that active site dynamics are coupled to the formation of a base triple within the branch-site helix that positions the 2'-OH of the adenosine for nucleophilic attack on the 5' scissile phosphate. This structure, complemented with biochemistry and comparative analyses to splicing complexes, supports a base triple model of adenosine recognition for branching within group II introns and the evolutionarily related spliceosome.

Generation of Rare Human NMDA Receptor Variants in Mice

The analysis of rare NMDAR gene variants in mice, coupled with a fundamental understanding of NMDAR function, plays a crucial role in achieving therapeutic success when addressing NMDAR dysfunctions in human patients. For the generation of such NMDAR mouse models, a basic knowledge of receptor structure, along with skills in database sequence analysis, cloning in E. coli, genetic manipulation of embryonic stem (ES) cells, and ultimately the genetic modification of mouse embryos, is essential. Primarily, this chapter will focus on the design and synthesis of NMDAR gene-targeting vectors that can be used successfully for the genetic manipulation of mice. We will outline the core principles of the design and synthesis of a gene targeting vector that facilitates the introduction of single-point mutations in NMDAR-encoding genes in mice. The transformation of ES cells, selection of positive ES cell colonies, manipulation of mouse embryos, and genotyping strategies will be described briefly.

Exploiting public databases of genomic variation to quantify evolutionary constraint on the branch point sequence in 30 plant and animal species

The branch point sequence is a degenerate intronic heptamer required for the assembly of the spliceosome during pre-mRNA splicing. Disruption of this motif may promote alternative splicing and eventually cause phenotype variation. Despite its functional relevance, the branch point sequence is not included in most genome annotations. Here, we predict branch point sequences in 30 plant and animal species and attempt to quantify their evolutionary constraints using public variant databases. We find an implausible variant distribution in the databases from 16 of 30 examined species. Comparative analysis of variants from whole-genome sequencing shows that variants submitted from exome sequencing or false positive variants are widespread in public databases and cause these irregularities. We then investigate evolutionary constraint with largely unbiased public variant databases in 14 species and find that the fourth and sixth position of the branch point sequence are more constrained than coding nucleotides. Our findings show that public variant databases should be scrutinized for possible biases before they qualify to analyze evolutionary constraint.

Trichothiodystrophy-associated MPLKIP maintains DBR1 levels for proper lariat debranching and ectodermal differentiation

The brittle hair syndrome Trichothiodystrophy (TTD) is characterized by variable clinical features, including photosensitivity, ichthyosis, growth retardation, microcephaly, intellectual disability, hypogonadism, and anaemia. TTD-associated mutations typically cause unstable mutant proteins involved in various steps of gene expression, severely reducing steady-state mutant protein levels. However, to date, no such link to instability of gene-expression factors for TTD-associated mutations in MPLKIP/TTDN1 has been established. Here, we present seven additional TTD individuals with MPLKIP mutations from five consanguineous families, with a newly identified MPLKIP variant in one family. By mass spectrometry-based interaction proteomics, we demonstrate that MPLKIP interacts with core splicing factors and the lariat debranching protein DBR1. MPLKIP-deficient primary fibroblasts have reduced steady-state DBR1 protein levels. Using Human Skin Equivalents (HSEs), we observed impaired keratinocyte differentiation associated with compromised splicing and eventually, an imbalanced proteome affecting skin development and, interestingly, also the immune system. Our data show that MPLKIP, through its DBR1 stabilizing role, is implicated in mRNA splicing, which is of particular importance in highly differentiated tissue.

Shell protein composition specified by the lncRNA NEAT1 domains dictates the formation of paraspeckles as distinct membraneless organelles

Many membraneless organelles (MLOs) formed through phase separation play crucial roles in various cellular processes. Although these MLOs co-exist in cells, how they maintain their independence without coalescence or engulfment remains largely unknown. Here, we investigated the molecular mechanism by which paraspeckles with core-shell architecture scaffolded by NEAT1_2 long noncoding RNAs exist as distinct MLOs. We identified NEAT1 deletion mutants that assemble paraspeckles that are incorporated into nuclear speckles. Several paraspeckle proteins, including SFPQ, HNRNPF and BRG1, prevent this incorporation and thus contribute to the segregation of paraspeckles from nuclear speckles. Shell localization of these proteins in the paraspeckles, which is determined by NEAT1_2 long noncoding RNA domains, is required for this segregation process. Conversely, U2-related spliceosomal proteins are involved in internalizing the paraspeckles into nuclear speckles. This study shows that the paraspeckle shell composition dictates the independence of MLOs in the nucleus, providing insights into the importance of the shell in defining features and functions of MLOs.

A novel splice variant of the human MSTN gene encodes a myostatin-specific myostatin inhibitor

BACKGROUND

Myostatin, encoded by the MSTN gene comprising 3 exons, is a potent negative regulator of skeletal muscle growth. Although a variety of myostatin inhibitors have been invented for increasing muscle mass in muscle wasting diseases, no effective inhibitor is currently available for clinical use. Myostatin isoforms in several animals have been reported to inhibit myostatin, but an isoform has never been identified for the human MSTN gene, a conserved gene among animals. Here, a splice variant of the human MSTN gene was explored.

METHODS

Transcripts and proteins were analysed by reverse transcription-PCR amplification and western blotting, respectively. Proteins were expressed from expression plasmid. Myostatin signalling was assayed by the SMAD-responsive luciferase activity. Cell proliferation was assayed by the Cell Counting Kit-8 (CCK-8) assay and cell counting. Cell cycle was analysed by the FastFUCCI system.

RESULTS

Reverse transcription-PCR amplification of the full-length MSTN transcript in CRL-2061 rhabdomyosarcoma cells revealed two bands consisting of a thick expected-size product and a thin additional small-size product. Sequencing of the small-size product showed a 963-bp deletion in the 5' end of exon 3, creating exon 3s, which contained unusual splice acceptor TG dinucleotides. The novel variant was identified in other human cell lines, although it was not identified in skeletal muscle. The 251-amino acid isoform encoded by the novel variant (myostatin-b) was identified in CRL-2061 rhabdomyosarcoma cells. Transfection of a myostatin-b expression plasmid into CRL-2061 and myoblast cells inhibited endogenous myostatin signalling (44%, P < 0.001 and 63%, P < 0.001, respectively). Furthermore, myostatin-b inhibited myostatin signalling induced by recombinant myostatin (68.8%, P < 0.001). In remarkable contrast, myostatin-b did not inhibit the myostatin signalling induced by recombinant growth differentiation factor 11 (9.2%, P = 0.70), transforming growth factor β (+3.1%, P = 0.83) or activin A (+1.1%, P = 0.96). These results indicate the myostatin-specific inhibitory effect of myostatin-b. Notably, the expression of myostatin-b in myoblasts significantly enhanced cell proliferation higher than the mock-transfected cells by the CCK-8 and direct cell counting assays (60%, P < 0.05 and 39%, P < 0.05, respectively). Myostatin-b increased the percentage of S-phase cells significantly higher than that of the mock-transfected cells (53% vs. 80%, P < 0.05).

CONCLUSIONS

We cloned a novel human MSTN variant produced by unorthodox splicing. The variant encoded a novel myostatin isoform, myostatin-b, that inhibited myostatin signalling by myostatin-specific manner and enhanced myoblast proliferation by shifting cell cycle. Myostatin-b, which has myostatin-specific inhibitory activity, could be developed as a natural myostatin inhibitor.

FexSplice: A LightGBM-Based Model for Predicting the Splicing Effect of a Single Nucleotide Variant Affecting the First Nucleotide G of an Exon

Single nucleotide variants (SNVs) affecting the first nucleotide G of an exon (Fex-SNVs) identified in various diseases are mostly recognized as missense or nonsense variants. Their effect on pre-mRNA splicing has been seldom analyzed, and no curated database is available. We previously reported that Fex-SNVs affect splicing when the length of the polypyrimidine tract is short or degenerate. However, we cannot readily predict the splicing effects of Fex-SNVs. We here scrutinized the available literature and identified 106 splicing-affecting Fex-SNVs based on experimental evidence. We similarly identified 106 neutral Fex-SNVs in the dbSNP database with a global minor allele frequency (MAF) of more than 0.01 and less than 0.50. We extracted 115 features representing the strength of splicing cis-elements and developed machine-learning models with support vector machine, random forest, and gradient boosting to discriminate splicing-affecting and neutral Fex-SNVs. Gradient boosting-based LightGBM outperformed the other two models, and the length and nucleotide compositions of the polypyrimidine tract played critical roles in the discrimination. Recursive feature elimination showed that the LightGBM model using 15 features achieved the best performance with an accuracy of 0.80 ± 0.12 (mean and SD), a Matthews Correlation Coefficient (MCC) of 0.57 ± 0.15, an area under the curve of the receiver operating characteristics curve (AUROC) of 0.86 ± 0.08, and an area under the curve of the precision-recall curve (AUPRC) of 0.87 ± 0.09 using a 10-fold cross-validation. We developed a web service program, named FexSplice that accepts a genomic coordinate either on GRCh37/hg19 or GRCh38/hg38 and returns a predicted probability of aberrant splicing of A, C, and T variants.

Circular stable intronic RNAs possess distinct biological features and are deregulated in bladder cancer

Until recently, intronic lariats were regarded as short-lasting splicing byproducts with no apparent function; however, increasing evidence of stable derivatives suggests regulatory roles. Yet little is known about their characteristics, functions, distribution, and expression in healthy and tumor tissue. Here, we profiled and characterized circular stable intronic sequence RNAs (sisRNAs) using total RNA-Seq data from bladder cancer (BC; n = 457, UROMOL cohort), healthy tissue (n = 46), and fractionated cell lines (n = 5). We found that the recently-discovered full-length intronic circles and the stable lariats formed distinct subclasses, with a surprisingly high intronic circle fraction in BC (∼45%) compared to healthy tissues (0-20%). The stable lariats and their host introns were characterized by small transcript sizes, highly conserved BP regions, enriched BP motifs, and localization in multiple cell fractions. Additionally, circular sisRNAs showed tissue-specific expression patterns. We found nine circular sisRNAs as differentially expressed across early-stage BC patients with different prognoses, and sisHNRNPK expression correlated with progression-free survival. In conclusion, we identify distinguishing biological features of circular sisRNAs and point to specific candidates (incl. sisHNRNPK, sisWDR13 and sisMBNL1) that were highly expressed, had evolutionary conserved sequences, or had clinical correlations, which may facilitate future studies and further insights into their functional roles.

Early Splicing Complexes and Human Disease

Over the last decade, our understanding of spliceosome structure and function has significantly improved, refining the study of the impact of dysregulated splicing on human disease. As a result, targeted splicing therapeutics have been developed, treating various diseases including spinal muscular atrophy and Duchenne muscular dystrophy. These advancements are very promising and emphasize the critical role of proper splicing in maintaining human health. Herein, we provide an overview of the current information on the composition and assembly of early splicing complexes-commitment complex and pre-spliceosome-and their association with human disease.

Recursive splicing discovery using lariats in total RNA sequencing

Recursive splicing is a non-canonical splicing mechanism in which an intron is removed in segments via multiple splicing reactions. Relatively few recursive splice sites have been identified with high confidence in human introns, and more comprehensive analyses are needed to better characterize where recursive splicing happens and whether or not it has a regulatory function. In this study, we use an unbiased approach using intron lariats to search for recursive splice sites in constitutive introns and alternative exons in the human transcriptome. We find evidence for recursive splicing in a broader range of intron sizes than previously reported and detail a new location for recursive splicing at the distal ends of cassette exons. In addition, we identify evidence for the conservation of these recursive splice sites among higher vertebrates and the use of these sites to influence alternative exon exclusion. Together, our data demonstrate the prevalence of recursive splicing and its potential influence on gene expression through alternatively spliced isoforms.

Structural basis of catalytic activation in human splicing

Pre-mRNA splicing follows a pathway driven by ATP-dependent RNA helicases. A crucial event of the splicing pathway is the catalytic activation, which takes place at the transition between the activated Bact and the branching-competent B* spliceosomes. Catalytic activation occurs through an ATP-dependent remodelling mediated by the helicase PRP2 (also known as DHX16)1-3. However, because PRP2 is observed only at the periphery of spliceosomes3-5, its function has remained elusive. Here we show that catalytic activation occurs in two ATP-dependent stages driven by two helicases: PRP2 and Aquarius. The role of Aquarius in splicing has been enigmatic6,7. Here the inactivation of Aquarius leads to the stalling of a spliceosome intermediate-the BAQR complex-found halfway through the catalytic activation process. The cryogenic electron microscopy structure of BAQR reveals how PRP2 and Aquarius remodel Bact and BAQR, respectively. Notably, PRP2 translocates along the intron while it strips away the RES complex, opens the SF3B1 clamp and unfastens the branch helix. Translocation terminates six nucleotides downstream of the branch site through an assembly of PPIL4, SKIP and the amino-terminal domain of PRP2. Finally, Aquarius enables the dissociation of PRP2, plus the SF3A and SF3B complexes, which promotes the relocation of the branch duplex for catalysis. This work elucidates catalytic activation in human splicing, reveals how a DEAH helicase operates and provides a paradigm for how helicases can coordinate their activities.

IARA: a complete and curated atlas of the biogenesis of spliceosome machinery during RNA splicing

Splicing is one of the most important post-transcriptional processing systems and is responsible for the generation of transcriptome diversity in all living eukaryotes. Splicing is regulated by the spliceosome machinery, which is responsible for each step of primary RNA processing. However, current molecules and stages involved in RNA splicing are still spread over different studies. Thus, a curated atlas of spliceosome-related molecules and all involved stages during RNA processing can provide all researchers with a reliable resource to better investigate this important mechanism. Here, we present IARA (website access: https://pucpr-bioinformatics.github.io/atlas/), an extensively curated and constantly updated catalog of molecules involved in spliceosome machinery. IARA has a map of the steps involved in the human splicing mechanism, and it allows a detailed overview of the molecules involved throughout the distinct steps of splicing.

A Trem2R47H mouse model without cryptic splicing drives age- and disease-dependent tissue damage and synaptic loss in response to plaques

BACKGROUND

The TREM2 R47H variant is one of the strongest genetic risk factors for late-onset Alzheimer's Disease (AD). Unfortunately, many current Trem2 R47H mouse models are associated with cryptic mRNA splicing of the mutant allele that produces a confounding reduction in protein product. To overcome this issue, we developed the Trem2R47H NSS (Normal Splice Site) mouse model in which the Trem2 allele is expressed at a similar level to the wild-type Trem2 allele without evidence of cryptic splicing products.

METHODS

Trem2R47H NSS mice were treated with the demyelinating agent cuprizone, or crossed with the 5xFAD mouse model of amyloidosis, to explore the impact of the TREM2 R47H variant on inflammatory responses to demyelination, plaque development, and the brain's response to plaques.

RESULTS

Trem2R47H NSS mice display an appropriate inflammatory response to cuprizone challenge, and do not recapitulate the null allele in terms of impeded inflammatory responses to demyelination. Utilizing the 5xFAD mouse model, we report age- and disease-dependent changes in Trem2R47H NSS mice in response to development of AD-like pathology. At an early (4-month-old) disease stage, hemizygous 5xFAD/homozygous Trem2R47H NSS (5xFAD/Trem2R47H NSS) mice have reduced size and number of microglia that display impaired interaction with plaques compared to microglia in age-matched 5xFAD hemizygous controls. This is associated with a suppressed inflammatory response but increased dystrophic neurites and axonal damage as measured by plasma neurofilament light chain (NfL) level. Homozygosity for Trem2R47H NSS suppressed LTP deficits and loss of presynaptic puncta caused by the 5xFAD transgene array in 4-month-old mice. At a more advanced (12-month-old) disease stage 5xFAD/Trem2R47H NSS mice no longer display impaired plaque-microglia interaction or suppressed inflammatory gene expression, although NfL levels remain elevated, and a unique interferon-related gene expression signature is seen. Twelve-month old Trem2R47H NSS mice also display LTP deficits and postsynaptic loss.

CONCLUSIONS

The Trem2R47H NSS mouse is a valuable model that can be used to investigate age-dependent effects of the AD-risk R47H mutation on TREM2 and microglial function including its effects on plaque development, microglial-plaque interaction, production of a unique interferon signature and associated tissue damage.

Long Non-Coding RNA Generated from CDKN1A Gene by Alternative Polyadenylation Regulates p21 Expression during DNA Damage Response

Alternative Polyadenylation (APA) is an emerging mechanism for dynamic changes in gene expression. Previously, we described widespread APA occurrence in introns during the DNA damage response (DDR). Here, we show that a DNA damage activated APA event occurs in the first intron of CDKN1A , inducing an alternate last exon (ALE)-containing lncRNA. We named this lncRNA SPUD (Selective Polyadenylation Upon Damage). SPUD localizes to polysomes in the cytoplasm and is detectable as multiple isoforms in available high throughput studies. SPUD has low abundance compared to the CDKN1A full-length isoform and is induced in cancer and normal cells under a variety of DNA damaging conditions in part through p53 transcriptional activation. RNA binding protein (RBP) HuR and the transcriptional repressor CTCF regulate SPUD levels. SPUD induction increases p21 protein, but not CDKN1A full-length levels, affecting p21 functions in cell-cycle, CDK2 expression, and cell viability. Like CDKN1A full-length isoform, SPUD can bind two competitive p21 translational regulators, the inhibitor calreticulin and the activator CUGBP1; SPUD can change their association with CDKN1A full-length in a DDR-dependent manner. Together, these results show a new regulatory mechanism by which a lncRNA controls p21 expression post-transcriptionally, highlighting lncRNA relevance in DDR progression and cellcycle.

The evolution of pre-mRNA splicing and its machinery revealed by reduced extremophilic red algae

The Cyanidiales are a group of mostly thermophilic and acidophilic red algae that thrive near volcanic vents. Despite their phylogenetic relationship, the reduced genomes of Cyanidioschyzon merolae and Galdieria sulphuraria are strikingly different with respect to pre-mRNA splicing, a ubiquitous eukaryotic feature. Introns are rare and spliceosomal machinery is extremely reduced in C. merolae, in contrast to G. sulphuraria. Previous studies also revealed divergent spliceosomes in the mesophilic red alga Porphyridium purpureum and the red algal derived plastid of Guillardia theta (Cryptophyta), along with unusually high levels of unspliced transcripts. To further examine the evolution of splicing in red algae, we compared C. merolae and G. sulphuraria, investigating splicing levels, intron position, intron sequence features, and the composition of the spliceosome. In addition to identifying 11 additional introns in C. merolae, our transcriptomic analysis also revealed typical eukaryotic splicing in G. sulphuraria, whereas most transcripts in C. merolae remain unspliced. The distribution of intron positions within their host genes was examined to provide insight into patterns of intron loss in red algae. We observed increasing variability of 5' splice sites and branch donor regions with increasing intron richness. We also found these relationships to be connected to reductions in and losses of corresponding parts of the spliceosome. Our findings highlight patterns of intron and spliceosome evolution in related red algae under the pressures of genome reduction.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

Alternative pre-mRNA splicing as a mechanism for terminating Toll-like Receptor signaling

While inflammation induced by Toll-like receptor (TLR) signaling is required to combat infection, persistent inflammation can damage host tissues and contribute to a myriad of acute and chronic inflammatory disorders. Thus, it is essential not only that TLR signaling be activated in the presence of pathogens but that TLR signaling is ultimately terminated. One mechanism that limits persistent TLR signaling is alternative pre-mRNA splicing. In addition to encoding the canonical mRNAs that produce proteins that promote inflammation, many genes in the TLR signaling pathway also encode alternative mRNAs that produce proteins that are dominant negative inhibitors of signaling. Many of these negative regulators are induced by immune challenge, so production of these alternative isoforms represents a negative feedback loop that limits persistent inflammation. While these alternative splicing events have been investigated on a gene by gene basis, there has been limited systemic analysis of this mechanism that terminates TLR signaling. Here we review what is known about the production of negatively acting alternative isoforms in the TLR signaling pathway including how these inhibitors function, how they are produced, and what role they may play in inflammatory disease.

Mechanism and modeling of human disease-associated near-exon intronic variants that perturb RNA splicing

It is estimated that 10%-30% of disease-associated genetic variants affect splicing. Splicing variants may generate deleteriously altered gene product and are potential therapeutic targets. However, systematic diagnosis or prediction of splicing variants is yet to be established, especially for the near-exon intronic splice region. The major challenge lies in the redundant and ill-defined branch sites and other splicing motifs therein. Here, we carried out unbiased massively parallel splicing assays on 5,307 disease-associated variants that overlapped with branch sites and collected 5,884 variants across the 5' splice region. We found that strong splice sites and exonic features preserve splicing from intronic sequence variation. Whereas the splice-altering mechanism of the 3' intronic variants is complex, that of the 5' is mainly splice-site destruction. Statistical learning combined with these molecular features allows precise prediction of altered splicing from an intronic variant. This statistical model provides the identity and ranking of biological features that determine splicing, which serves as transferable knowledge and out-performs the benchmarking predictive tool. Moreover, we demonstrated that intronic splicing variants may associate with disease risks in the human population. Our study elucidates the mechanism of splicing response of intronic variants, which classify disease-associated splicing variants for the promise of precision medicine.

Genome-wide detection of human variants that disrupt intronic branchpoints

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter.

Structural Context of a Critical Exon of Spinal Muscular Atrophy Gene

Humans contain two nearly identical copies of Survival Motor Neuron genes, SMN1 and SMN2. Deletion or mutation of SMN1 causes spinal muscular atrophy (SMA), one of the leading genetic diseases associated with infant mortality. SMN2 is unable to compensate for the loss of SMN1 due to predominant exon 7 skipping, leading to the production of a truncated protein. Antisense oligonucleotide and small molecule-based strategies aimed at the restoration of SMN2 exon 7 inclusion are approved therapies of SMA. Many cis-elements and transacting factors have been implicated in regulation of SMN exon 7 splicing. Also, several structural elements, including those formed by a long-distance interaction, have been implicated in the modulation of SMN exon 7 splicing. Several of these structures have been confirmed by enzymatic and chemical structure-probing methods. Additional structures formed by inter-intronic interactions have been predicted by computational algorithms. SMN genes generate a vast repertoire of circular RNAs through inter-intronic secondary structures formed by inverted Alu repeats present in large number in SMN genes. Here, we review the structural context of the exonic and intronic cis-elements that promote or prevent exon 7 recognition. We discuss how structural rearrangements triggered by single nucleotide substitutions could bring drastic changes in SMN2 exon 7 splicing. We also propose potential mechanisms by which inter-intronic structures might impact the splicing outcomes.

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

From the comparative study of a circRNA originating from an mammalian ATXN2L intron to understanding the genesis of intron lariat-derived circRNAs

Circular intronic RNAs (ciRNAs) are still unexplored regarding mechanisms for their emergence. We considered the ATXN2L intron lariat-derived circular RNA (ciRNA-ATXN2L) as an opportunity to conduct a cross-species examination of ciRNA genesis. To this end, we investigated 207 datasets from 4 tissues and from 13 mammalian species. While in eight species, ciRNA-ATXN2L was never detected, in pigs and rabbits, ciRNA-ATXN2L was expressed in all tissues and sometimes at very high levels. Bovine tissues were an intermediate case and in macaques and cats, only ciRNA-ATXN2L traces were detected. The pattern of ciRNA-ATXN2L restricted to only five species is not related to a particular evolution of intronic sequences. To empower our analysis, we considered 221 additional introns including 80 introns where a lariat-derived ciRNA was previously described. The primary driver of micro-ciRNA genesis (< 155 nt as ciRNA-ATXN2L) appears to be the absence of a canonical "A" (i.e. a "tnA" located in the usual branching region) to build the lariat around this adenosine. The balance between available "non canonical-A" (no ciRNA genesis) and "non-A" (ciRNA genesis) for use as a branch point to build the lariat could modify the expression level of ciRNA-ATXN2L. In addition, the rare localization of the 2'-5' bond in an open RNA secondary structure could also negatively affect the lifetime of ciRNAs (macaque ciRNA-ATXN2L). Our analyses suggest that ciRNA-ATXN2L is likely a functionless splice remnant. This study provides a better understanding of the ciRNAs origin, especially drivers for micro ciRNA genesis.

Somatic Mutations in Core Spliceosome Components Promote Tumorigenesis and Generate an Exploitable Vulnerability in Human Cancer

Simple Summary

High throughput exome sequencing approaches have disclosed recurrent cancer-associated mutations in spliceosomal components, which drive aberrant pre-mRNA processing events and support the tumor phenotype. At the same time, mutations in spliceosome genes and aberrant splicing regulation establish a selective vulnerability of cancer cells to splicing-targeting approaches, which could be exploited therapeutically. It is conceivable that a better understanding of the mechanisms and roles of abnormal splicing in tumor metabolism will facilitate the development of a novel generation of tumor-targeting drugs. In this review, we describe recent advances in the elucidation of the biological impact and biochemical effects of somatic mutations in core spliceosome components on splicing choices and their associated targetable vulnerabilities.

Abstract

Alternative pre-mRNA processing enables the production of distinct mRNA and protein isoforms from a single gene, thus greatly expanding the coding potential of eukaryotic genomes and fine-tuning gene expression programs. Splicing is carried out by the spliceosome, a complex molecular machinery which assembles step-wise on mRNA precursors in the nucleus of eukaryotic cells. In the last decade, exome sequencing technologies have allowed the identification of point mutations in genes encoding splicing factors as a recurrent hallmark of human cancers, with higher incidence in hematological malignancies. These mutations lead to production of splicing factors that reduce the fidelity of the splicing process and yield splicing variants that are often advantageous for cancer cells. However, at the same time, these mutations increase the sensitivity of transformed cells to splicing inhibitors, thus offering a therapeutic opportunity for novel targeted strategies. Herein, we review the recent literature documenting cancer-associated mutations in components of the early spliceosome complex and discuss novel therapeutic strategies based on small-molecule spliceosome inhibitors that exhibit strong anti-tumor effects, particularly against cancer cells harboring mutations in spliceosomal components.

ABCA4 c.859-25A>G, a Frequent Palestinian Founder Mutation Affecting the Intron 7 Branchpoint, Is Associated With Early-Onset Stargardt Disease

Purpose

The effect of noncoding variants is often unknown in the absence of functional assays. Here, we characterized an ABCA4 intron 7 variant, c.859-25A>G, identified in Palestinian probands with Stargardt disease (STGD) or cone-rod dystrophy (CRD). We investigated the effect of this variant on the ABCA4 mRNA and retinal phenotype, and its prevalence in Palestine.

Methods

The ABCA4 gene was sequenced completely or partially in 1998 cases with STGD or CRD. The effect of c.859-25A>G on splicing was investigated in silico using SpliceAI and in vitro using splice assays. Homozygosity mapping was performed for 16 affected individuals homozygous for c.859-25A>G. The clinical phenotype was assessed using functional and structural analyses including visual acuity, full-field electroretinography, and multimodal imaging.

Results

The smMIPs-based ABCA4 sequencing revealed c.859-25A>G in 10 Palestinian probands from Hebron and Jerusalem. SpliceAI predicted a significant effect of this putative branchpoint-inactivating variant on the nearby intron 7 splice acceptor site. Splice assays revealed exon 8 skipping and two partial inclusions of intron 7, each having a deleterious effect. Additional genotyping revealed another 46 affected homozygous or compound heterozygous individuals carrying variant c.859-25A>G. Homozygotes shared a genomic segment of 59.6 to 87.9 kb and showed severe retinal defects on ophthalmoscopic evaluation.

Conclusions

The ABCA4 variant c.859-25A>G disrupts a predicted branchpoint, resulting in protein truncation because of different splice defects, and is associated with early-onset STGD1 when present in homozygosity. This variant was found in 25/525 Palestinian inherited retinal dystrophy probands, representing one of the most frequent inherited retinal disease-causing variants in West-Bank Palestine.

Non-canonical splice junction processing increases the diversity of RBFOX2 splicing isoforms

The underlying mechanisms of splicing regulation through non-canonical splice junction processing remain largely unknown. Here, we identified two RBFOX2 splicing isoforms by alternative 3' splice site selection of exon 9; the non-canonical splice junction processed RBFOX2 transcript (RBFOX2-N.C.) was expressed by the selection of the 3' splice GG acceptor sequence. The cytoplasmic localization of RBFOX2-C., a canonical splice junction-processed RBFOX2 transcript, was different from that of RBFOX2-N.C., which showed nuclear localization. In addition, we confirmed that RBFOX2-C. showed a significantly stronger localization into stress granules than RBFOX2-N.C. upon sodium arsenite treatment. Next, we investigated the importance of non-canonical 3' splice GG sequence selection of specific cis-regulatory elements using minigene constructs of the RBFOX2 gene. We found that the non-canonical 3' splice GG sequence and suboptimal branch point site adjacent region were critical for RBFOX2-N.C. expression through a non-canonical 3' splice selection. Our results suggest a regulatory mechanism for the non-canonical 3' splice selection in the RBFOX2 gene, providing a basis for studies related to the regulation of alternative pre-mRNA splicing through non-canonical splice junction processing.

The splicing effect of variants at branchpoint elements in cancer genes

PURPOSE

Branchpoint elements are required for intron removal, and variants at these elements can result in aberrant splicing. We aimed to assess the value of branchpoint annotations generated from recent large-scale studies to select branchpoint-abrogating variants, using hereditary cancer genes as model.

METHODS

We identified branchpoint elements in 119 genes associated with hereditary cancer from 3 genome-wide experimentally-inferred and 2 predicted branchpoint data sets. We then identified variants that occur within branchpoint elements from public databases. We compared conservation, unique variant observations, and population frequencies at different nucleotides within branchpoint motifs. Finally, selected minigene assays were performed to assess the splicing effect of variants at branchpoint elements within mismatch repair genes.

RESULTS

There was poor overlap between predicted and experimentally-inferred branchpoints. Our analysis of cancer genes suggested that variants at -2 nucleotide, -1 nucleotide, and branchpoint positions in experimentally-inferred canonical motifs are more likely to be clinically relevant. Minigene assay data showed the -2 nucleotide to be more important to branchpoint motif integrity but also showed fluidity in branchpoint usage.

CONCLUSION

Data from cancer gene analysis suggest that there are few high-risk alleles that severely impact function via branchpoint abrogation. Results of this study inform a general scheme to prioritize branchpoint motif variants for further study.

Ribosomal DNA promoter recognition is determined in vivo by cooperation between UBTF1 and SL1 and is compromised in the UBTF-E210K neuroregression syndrome

Transcription of the ~200 mouse and human ribosomal RNA genes (rDNA) by RNA Polymerase I (RPI/PolR1) accounts for 80% of total cellular RNA, around 35% of all nuclear RNA synthesis, and determines the cytoplasmic ribosome complement. It is therefore a major factor controlling cell growth and its misfunction has been implicated in hypertrophic and developmental disorders. Activation of each rDNA repeat requires nucleosome replacement by the architectural multi-HMGbox factor UBTF to create a 15.7 kbp nucleosome free region (NFR). Formation of this NFR is also essential for recruitment of the TBP-TAFI factor SL1 and for preinitiation complex (PIC) formation at the gene and enhancer-associated promoters of the rDNA. However, these promoters show little sequence commonality and neither UBTF nor SL1 display significant DNA sequence binding specificity, making what drives PIC formation a mystery. Here we show that cooperation between SL1 and the longer UBTF1 splice variant generates the specificity required for rDNA promoter recognition in cell. We find that conditional deletion of the TAF1B subunit of SL1 causes a striking depletion of UBTF at both rDNA promoters but not elsewhere across the rDNA. We also find that while both UBTF1 and -2 variants bind throughout the rDNA NFR, only UBTF1 is present with SL1 at the promoters. The data strongly suggest an induced-fit model of RPI promoter recognition in which UBTF1 plays an architectural role. Interestingly, a recurrent UBTF-E210K mutation and the cause of a pediatric neurodegeneration syndrome provides indirect support for this model. E210K knock-in cells show enhanced levels of the UBTF1 splice variant and a concomitant increase in active rDNA copies. In contrast, they also display reduced rDNA transcription and promoter recruitment of SL1. We suggest the underlying cause of the UBTF-E210K syndrome is therefore a reduction in cooperative UBTF1-SL1 promoter recruitment that may be partially compensated by enhanced rDNA activation.

Differential fates of introns in gene expression due to global alternative splicing

The discovery of introns over four decades ago revealed a new vision of genes and their interrupted arrangement. Throughout the years, it has appeared that introns play essential roles in the regulation of gene expression. Unique processing of excised introns through the formation of lariats suggests a widespread role for these molecules in the structure and function of cells. In addition to rapid destruction, these lariats may linger on in the nucleus or may even be exported to the cytoplasm, where they remain stable circular RNAs (circRNAs). Alternative splicing (AS) is a source of diversity in mature transcripts harboring retained introns (RI-mRNAs). Such RNAs may contain one or more entire retained intron(s) (RIs), but they may also have intron fragments resulting from sequential excision of smaller subfragments via recursive splicing (RS), which is characteristic of long introns. There are many potential fates of RI-mRNAs, including their downregulation via nuclear and cytoplasmic surveillance systems and the generation of new protein isoforms with potentially different functions. Various reports have linked the presence of such unprocessed transcripts in mammals to important roles in normal development and in disease-related conditions. In certain human neurological-neuromuscular disorders, including myotonic dystrophy type 2 (DM2), frontotemporal dementia/amyotrophic lateral sclerosis (FTD/ALS) and Duchenne muscular dystrophy (DMD), peculiar processing of long introns has been identified and is associated with their pathogenic effects. In this review, we discuss different mechanisms involved in the processing of introns during AS and the functions of these large sections of the genome in our biology.

Phylogeny and conservation of plant U2A/U2A', a core splicing component in U2 spliceosomal complex

This study systematically identifies 112 U2A genes from 80 plant species by combinatory bioinformatics analysis, which is important for understanding their phylogenetic history, expression profiles and for predicting specific functions. In eukaryotes, a pre-mRNA can generate multiple transcripts by removing certain introns and joining corresponding exons, thus greatly expanding the transcriptome and proteome diversity. The spliceosome is a mega-Dalton ribonucleoprotein (RNP) complex that is essential for the process of splicing. In spliceosome components, the U2 small nuclear ribonucleoprotein (U2 snRNP) forms the pre-spliceosome by association with the branch site. An essential component that promotes U2 snRNP assembly, named U2A, has been extensively identified in humans, yeast and nematodes. However, studies examining U2A genes in plants are scarce. In this study, we performed a comprehensive analysis and identified a total of 112 U2A genes from 80 plant species representing dicots, monocots, mosses and algae. Comparisons of the gene structures, protein domains, and expression patterns of 112 U2A genes indicated that the conserved functions were likely retained by plant U2A genes and important for responses to internal and external stimuli. In addition, analysis of alternative transcripts and splice sites of U2A genes indicated that the fifth intron contained a conserved alternative splicing event that might be important for its molecular function. Our work provides a general understanding of this splicing factor family in terms of genes and proteins, and it will serve as a fundamental resource that will contribute to further mechanistic characterization in plants.

Spliceator: multi-species splice site prediction using convolutional neural networks

Background

Ab initio prediction of splice sites is an essential step in eukaryotic genome annotation. Recent predictors have exploited Deep Learning algorithms and reliable gene structures from model organisms. However, Deep Learning methods for non-model organisms are lacking.

Results

We developed Spliceator to predict splice sites in a wide range of species, including model and non-model organisms. Spliceator uses a convolutional neural network and is trained on carefully validated data from over 100 organisms. We show that Spliceator achieves consistently high accuracy (89–92%) compared to existing methods on independent benchmarks from human, fish, fly, worm, plant and protist organisms.

Conclusions

Spliceator is a new Deep Learning method trained on high-quality data, which can be used to predict splice sites in diverse organisms, ranging from human to protists, with consistently high accuracy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04471-3.

Evolution of the Early Spliceosomal Complex-From Constitutive to Regulated Splicing

Pre-mRNA splicing is a major process in the regulated expression of genes in eukaryotes, and alternative splicing is used to generate different proteins from the same coding gene. Splicing is a catalytic process that removes introns and ligates exons to create the RNA sequence that codifies the final protein. While this is achieved in an autocatalytic process in ancestral group II introns in prokaryotes, the spliceosome has evolved during eukaryogenesis to assist in this process and to finally provide the opportunity for intron-specific splicing. In the early stage of splicing, the RNA 5' and 3' splice sites must be brought within proximity to correctly assemble the active spliceosome and perform the excision and ligation reactions. The assembly of this first complex, termed E-complex, is currently the least understood process. We focused in this review on the formation of the E-complex and compared its composition and function in three different organisms. We highlight the common ancestral mechanisms in S. cerevisiae, S. pombe, and mammals and conclude with a unifying model for intron definition in constitutive and regulated co-transcriptional splicing.