RNA structure and the mechanisms of alternative splicing

DDX RNA helicases: key players in cellular homeostasis and innate antiviral immunity

RNA helicases are integral in RNA metabolism, performing important roles in cellular homeostasis and stress responses. In particular, the DExD/H-box (DDX) helicase family possesses a conserved catalytic core that binds structural features rather than specific sequences in RNA targets. DDXs have critical roles in all aspects of RNA metabolism including ribosome biogenesis, translation, RNA export, and RNA stability. Importantly, functional specialization within this family arises from divergent N and C termini and is driven at least in part by gene duplications with 18 of the 42 human helicases having paralogs. In addition to their key roles in the homeostatic control of cellular RNA, these factors have critical roles in RNA virus infection. The canonical RIG-I-like receptors (RLRs) play pivotal roles in cytoplasmic sensing of viral RNA structures, inducing antiviral gene expression. Additional RNA helicases function as viral sensors or regulators, further diversifying the innate immune defense arsenal. Moreover, some of these helicases have been coopted by viruses to facilitate their replication. Altogether, DDX helicases exhibit functional specificity, playing intricate roles in RNA metabolism and host defense. This review will discuss the mechanisms by which these RNA helicases recognize diverse RNA structures in cellular and viral RNAs, and how this impacts RNA processing and innate immune responses.

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We found frequent splicing order differences between alleles and uncovered significant single nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This included SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also had a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, showed the most allele-specific splicing orders, which frequently co-occurred with allele-specific AS, APA or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Splicing factor hnRNPA1 regulates alternative splicing of LOXL2 to enhance the production of LOXL2Δ13

Lysyl oxidase-like 2 (LOXL2) is a member of the lysyl oxidase family and has the ability to catalyze the cross-linking of extracellular matrix collagen and elastin. High expression of LOXL2 is related to tumor cell proliferation, invasion, and metastasis. LOXL2 contains 14 exons. Previous studies have found that LOXL2 has abnormal alternative splicing and exon skipping in a variety of tissues and cells, resulting in a new alternatively spliced isoform denoted LOXL2Δ13. LOXL2Δ13 lacks LOXL2WT exon 13, but its encoded protein has greater ability to induce tumor cell proliferation, invasion, and metastasis. However, the molecular events that produce LOXL2Δ13 are still unclear. In this study, we found that overexpression of the splicing factor hnRNPA1 in cells can regulate the alternative splicing of LOXL2 and increase the expression of LOXL2Δ13. The exonic splicing silencer exists at the 3' splice site and 5' splice site of LOXL2 exon 13. HnRNPA1 can bind to the exonic splicing silencer and inhibit the inclusion of exon 13. The RRM domain of hnRNPA1 and phosphorylation of hnRNPA1 at S91 and S95 are important for the regulation of LOXL2 alternative splicing. These results show that hnRNPA1 is a splicing factor that enhances the production of LOXL2Δ13.

Identification of RNA structures and their roles in RNA functions

The development of high-throughput RNA structure profiling methods in the past decade has greatly facilitated our ability to map and characterize different aspects of RNA structures transcriptome-wide in cell populations, single cells and single molecules. The resulting high-resolution data have provided insights into the static and dynamic nature of RNA structures, revealing their complexity as they perform their respective functions in the cell. In this Review, we discuss recent technical advances in the determination of RNA structures, and the roles of RNA structures in RNA biogenesis and functions, including in transcription, processing, translation, degradation, localization and RNA structure-dependent condensates. We also discuss the current understanding of how RNA structures could guide drug design for treating genetic diseases and battling pathogenic viruses, and highlight existing challenges and future directions in RNA structure research.

LncRNA-mediated orchestrations of alternative splicing in the landscape of breast cancer

Alternative splicing (AS) is a fundamental post-transcriptional process in eukaryotes, enabling a single gene to generate diverse mRNA transcripts, thereby enhancing protein variability. This process involves the excision of introns and the joining of exons in pre-mRNA(s) to form mature mRNA. The resulting mature mRNAs exhibit various combinations of exons, contributing to functional diversity. Dysregulation of AS can substantially modulate protein functions, impacting the onset and progression of numerous diseases, including cancer. Non-coding RNAs (ncRNAs) are distinct from protein-coding RNAs and consist of short and long types. Long non-coding RNAs (lncRNAs) play an important role in regulating several cellular processes, particularly alternative splicing, according to new research. This review provides insight into the latest discoveries concerning how lncRNAs influence alternative splicing within the realm of breast cancer. Additionally, it explores potential therapeutic strategies focused on targeting lncRNAs.

Unearthing a novel function of SRSF1 in binding and unfolding of RNA G-quadruplexes

SRSF1 governs splicing of over 1500 mRNA transcripts. SRSF1 contains two RNA-recognition motifs (RRMs) and a C-terminal Arg/Ser-rich region (RS). It has been thought that SRSF1 RRMs exclusively recognize single-stranded exonic splicing enhancers, while RS lacks RNA-binding specificity. With our success in solving the insolubility problem of SRSF1, we can explore the unknown RNA-binding landscape of SRSF1. We find that SRSF1 RS prefers purine over pyrimidine. Moreover, SRSF1 binds to the G-quadruplex (GQ) from the ARPC2 mRNA, with both RRMs and RS being crucial. Our binding assays show that the traditional RNA-binding sites on the RRM tandem and the Arg in RS are responsible for GQ binding. Interestingly, our FRET and circular dichroism data reveal that SRSF1 unfolds the ARPC2 GQ, with RS leading unfolding and RRMs aiding. Our saturation transfer difference NMR results discover that Arg residues in SRSF1 RS interact with the guanine base but not other nucleobases, underscoring the uniqueness of the Arg/guanine interaction. Our luciferase assays confirm that SRSF1 can alleviate the inhibitory effect of GQ on gene expression in the cell. Given the prevalence of RNA GQ and SR proteins, our findings unveil unexplored SR protein functions with broad implications in RNA splicing and translation.

Identification of MYC intron 2 regions that modulate expression

MYC pre-mRNA is spliced with high fidelity to produce the transcription factor known to regulate cellular differentiation, proliferation, apoptosis, and alternative splicing. The mechanisms underpinning the pre-mRNA splicing of MYC, however, remain mostly unexplored. In this study, we examined the interaction of heterogeneous nuclear ribonucleoprotein C (HNRNPC) with MYC intron 2. Building off published eCLIP studies, we confirmed this interaction with poly(U) regions in intron 2 of MYC and found that full binding is correlated with optimal protein production. The interaction appears to be compensatory, as mutational disruption of all three poly(U) regions was required to reduce both HNRNPC binding capacity and fidelity of either splicing or translation. Poly(U) sequences in MYC intron 2 were relatively conserved across sequences from several different species. Lastly, we identified a short sequence just upstream of an HNRNPC binding region that when removed enhances MYC translation.

Unearthing SRSF1's Novel Function in Binding and Unfolding of RNA G-Quadruplexes

SRSF1 governs splicing of over 1,500 mRNA transcripts. SRSF1 contains two RNA-recognition motifs (RRMs) and a C-terminal Arg/Ser-rich region (RS). It has been thought that SRSF1 RRMs exclusively recognize single-stranded exonic splicing enhancers, while RS lacks RNA-binding specificity. With our success in solving the insolubility problem of SRSF1, we can explore the unknown RNA-binding landscape of SRSF1. We find that SRSF1 RS prefers purine over pyrimidine. Moreover, SRSF1 binds to the G-quadruplex (GQ) from the ARPC2 mRNA, with both RRMs and RS being crucial. Our binding assays show that the traditional RNA-binding sites on the RRM tandem and the Arg in RS are responsible for GQ binding. Interestingly, our FRET and circular dichroism data reveal that SRSF1 unfolds the ARPC2 GQ, with RS leading unfolding and RRMs aiding. Our saturation transfer difference NMR results discover that Arg residues in SRSF1 RS interact with the guanine base but other nucleobases, underscoring the uniqueness of the Arg/guanine interaction. Our luciferase assays confirm that SRSF1 can alleviate the inhibitory effect of GQ on gene expression in the cell. Given the prevalence of RNA GQ and SR proteins, our findings unveil unexplored SR protein functions with broad implications in RNA splicing and translation.

Genes, genomes, and developmental process

The view advanced by Madole & Harden falls back on the dogma of a gene as a DNA sequence that codes for a fixed product with an invariant function regardless of temporal and spatial contexts. This outdated perspective entrenches the metaphor of genes as static units of information and glosses over developmental complexities.

HCS-Splice: A High-Content Screening Method to Advance the Discovery of RNA Splicing-Modulating Therapeutics

Nucleic acid therapeutics have demonstrated an impressive acceleration in recent years. They work through multiple mechanisms of action, including the downregulation of gene expression and the modulation of RNA splicing. While several drugs based on the former mechanism have been approved, few target the latter, despite the promise of RNA splicing modulation. To improve our ability to discover novel RNA splicing-modulating therapies, we developed HCS-Splice, a robust cell-based High-Content Screening (HCS) assay. By implementing the use of a two-colour (GFP/RFP) fluorescent splicing reporter plasmid, we developed a versatile, effective, rapid, and robust high-throughput strategy for the identification of potent splicing-modulating molecules. The HCS-Splice strategy can also be used to functionally confirm splicing mutations in human genetic disorders or to screen drug candidates. As a proof-of-concept, we introduced a dementia-related splice-switching mutation in the Microtubule-Associated Protein Tau (MAPT) exon 10 splicing reporter. We applied HCS-Splice to the wild-type and mutant reporters and measured the functional change in exon 10 inclusion. To demonstrate the applicability of the method in cell-based drug discovery, HCS-Splice was used to evaluate the efficacy of an exon 10-targeting siRNA, which was able to restore the correct alternative splicing balance.

A network of DZF proteins controls alternative splicing regulation and fidelity

Proteins containing DZF (domain associated with zinc fingers) modules play important roles throughout gene expression, from transcription to translation. Derived from nucleotidyltransferases but lacking catalytic residues, DZF domains serve as heterodimerization surfaces between DZF protein pairs. Three DZF proteins are widely expressed in mammalian tissues, ILF2, ILF3 and ZFR, which form mutually exclusive ILF2-ILF3 and ILF2-ZFR heterodimers. Using eCLIP-Seq, we find that ZFR binds across broad intronic regions to regulate the alternative splicing of cassette and mutually exclusive exons. ZFR preferentially binds dsRNA in vitro and is enriched on introns containing conserved dsRNA elements in cells. Many splicing events are similarly altered upon depletion of any of the three DZF proteins; however, we also identify independent and opposing roles for ZFR and ILF3 in alternative splicing regulation. Along with widespread involvement in cassette exon splicing, the DZF proteins control the fidelity and regulation of over a dozen highly validated mutually exclusive splicing events. Our findings indicate that the DZF proteins form a complex regulatory network that leverages dsRNA binding by ILF3 and ZFR to modulate splicing regulation and fidelity.

New insights into the centrosome-associated spliceosome components as regulators of ciliogenesis and tissue identity

Biomolecular condensates are membrane-less assemblies of proteins and nucleic acids. Centrosomes are biomolecular condensates that play a crucial role in nuclear division, cytoskeletal remodeling, and cilia formation in animal cells. Spatial omics technology is providing new insights into the dynamic exchange of spliceosome components between the nucleus and the centrosome/cilium. Intriguingly, centrosomes are emerging as cytoplasmic sites for information storage, enriched with RNA molecules and RNA-processing proteins. Furthermore, growing evidence supports the view that nuclear spliceosome components assembled at the centrosome function as potential coordinators of splicing subprograms, pluripotency, and cell differentiation. In this article, we first discuss the current understanding of the centrosome/cilium complex, which controls both stem cell differentiation and pluripotency. We next explore the molecular mechanisms that govern splicing factor assembly and disassembly around the centrosome and examine how RNA processing pathways contribute to ciliogenesis. Finally, we discuss numerous unresolved compelling questions regarding the centrosome-associated spliceosome components and transcript variants within the cytoplasm as sources of RNA-based secondary messages in the regulation of cell identity and cell fate determination. This article is categorized under: RNA-Based Catalysis > RNA Catalysis in Splicing and Translation RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Processing > Splicing Regulation/Alternative Splicing RNA Processing > RNA Processing.

An automated framework for evaluation of deep learning models for splice site predictions

A novel framework for the automated evaluation of various deep learning-based splice site detectors is presented. The framework eliminates time-consuming development and experimenting activities for different codebases, architectures, and configurations to obtain the best models for a given RNA splice site dataset. RNA splicing is a cellular process in which pre-mRNAs are processed into mature mRNAs and used to produce multiple mRNA transcripts from a single gene sequence. Since the advancement of sequencing technologies, many splice site variants have been identified and associated with the diseases. So, RNA splice site prediction is essential for gene finding, genome annotation, disease-causing variants, and identification of potential biomarkers. Recently, deep learning models performed highly accurately for classifying genomic signals. Convolutional Neural Network (CNN), Long Short-Term Memory (LSTM) and its bidirectional version (BLSTM), Gated Recurrent Unit (GRU), and its bidirectional version (BGRU) are promising models. During genomic data analysis, CNN's locality feature helps where each nucleotide correlates with other bases in its vicinity. In contrast, BLSTM can be trained bidirectionally, allowing sequential data to be processed from forward and reverse directions. Therefore, it can process 1-D encoded genomic data effectively. Even though both methods have been used in the literature, a performance comparison was missing. To compare selected models under similar conditions, we have created a blueprint for a series of networks with five different levels. As a case study, we compared CNN and BLSTM models' learning capabilities as building blocks for RNA splice site prediction in two different datasets. Overall, CNN performed better with [Formula: see text] accuracy ([Formula: see text] improvement), [Formula: see text] F1 score ([Formula: see text] improvement), and [Formula: see text] AUC-PR ([Formula: see text] improvement) in human splice site prediction. Likewise, an outperforming performance with [Formula: see text] accuracy ([Formula: see text] improvement), [Formula: see text] F1 score ([Formula: see text] improvement), and [Formula: see text] AUC-PR ([Formula: see text] improvement) is achieved in C. elegans splice site prediction. Overall, our results showed that CNN learns faster than BLSTM and BGRU. Moreover, CNN performs better at extracting sequence patterns than BLSTM and BGRU. To our knowledge, no other framework is developed explicitly for evaluating splice detection models to decide the best possible model in an automated manner. So, the proposed framework and the blueprint would help selecting different deep learning models, such as CNN vs. BLSTM and BGRU, for splice site analysis or similar classification tasks and in different problems.

Investigating the NRAS 5' UTR as a target for small molecules

Neuroblastoma RAS (NRAS) is an oncogene that is deregulated and highly mutated in cancers including melanomas and acute myeloid leukemias. The 5' untranslated region (UTR) (5' UTR) of the NRAS mRNA contains a G-quadruplex (G4) that regulates translation. Here we report a novel class of small molecule that binds to the G4 structure located in the 5' UTR of the NRAS mRNA. We used a small molecule microarray screen to identify molecules that selectively bind to the NRAS-G4 with submicromolar affinity. One compound inhibits the translation of NRAS in vitro but showed only moderate effects on the NRAS levels in cellulo. Rapid Amplification of cDNA Ends and RT-PCR analysis revealed that the predominant NRAS transcript does not possess the G4 structure. Thus, although NRAS transcripts lack a G4 in many cell lines the concept of targeting folded regions within 5' UTRs to control translation remains a highly attractive strategy.

LncRNA BC promotes lung adenocarcinoma progression by modulating IMPAD1 alternative splicing

BACKGROUND

The therapeutic value of targeted therapies in patients with lung cancer is reduced when tumours acquire secondary resistance after an initial period of successful treatment. However, the molecular events behind the resistance to targeted therapies in lung cancer remain largely unknown.

AIMS

To discover the important role and mechanism of lncRNA BC in promoting tumor metastasis and influencing clinical prognosis of LUAD.

MATERIALS & METHODS

Microarrays were used to screen a comprehensive set of lncRNAs with differential expression profiles in lung cancer cells. The functional role and mechanism of lncRNA were further investigated by gain- and loss-of-function assays. RNA pull-down, protein assays, and mass spectrometry were used to identify proteins that interacted with lncRNA. TaqMan PCR was used to measure lncRNA in lung adenocarcinoma and adjacent nontumor tissues from 428 patients. The clinical significance of lncRNA identified was statistically confirmed in this cohort of patients.

RESULTS

In this study, we show that the long non-coding RNA BC009639 (BC) is involved in acquired resistance to EGFR-targeted therapies. Among the 235 long non-coding RNAs that were differentially expressed in lung cancer cell lines, with different metastatic potentials, BC promoted growth, invasion, metastasis, and resistance to EGFR-tyrosine kinase inhibitors (EGFR-TKIs), both in vitro and in vivo. BC was highly expressed in 428 patients with lung adenocarcinoma (LUAD) and high BC expression correlated with reduced efficacy of EGFR-TKI therapy. To uncover the molecular mechanism of BC-mediated EGFR-TKI resistance in lung cancer, we screened and identified nucleolin and hnRNPK that interact with BC. BC formed the splicing complex with nucleolin and hnRNPK to facilitate the production of a non-protein-coding inositol monophosphatase domain containing 1 (IMPAD1) splice variant, instead of the protein-coding variant. The BC-mediated alternative splicing (AS) of IMPAD1 resulted in the induction of the epithelial-mesenchymal transition and resistance to EGFR-TKI in lung cancer. High BC expression correlated with clinical progress and poor survival among 402 patients with LUAD.

DISSCUSSION

Through alternative splicing, BC boosted the non-coding IMPAD1-203 transcript variant while suppressing the IMPAD1-201 variant. In order to control the processing of pre-mRNA, BC not only attracted RNA binding proteins (NCL, IGF2BP1) or splicing factors (hnRNPK), but also controlled the formation of the splicing-regulator complex by creating RNA-RNA-duplexes.

CONCLUSION

Our results reveal an important role for BC in mediating resistance to EGFR-targeted therapy in LUAD through IMPAD1 AS and in implication for the targeted therapy resistance.

Modulation of pre-mRNA structure by hnRNP proteins regulates alternative splicing of MALT1

Alternative splicing plays key roles for cell type-specific regulation of protein function. It is controlled by cis-regulatory RNA elements that are recognized by RNA binding proteins (RBPs). The MALT1 paracaspase is a key factor of signaling pathways that mediate innate and adaptive immune responses. Alternative splicing of MALT1 is critical for controlling optimal T cell activation. We demonstrate that MALT1 splicing depends on RNA structural elements that sequester the splice sites of the alternatively spliced exon7. The RBPs hnRNP U and hnRNP L bind competitively to stem-loop RNA structures that involve the 5' and 3' splice sites flanking exon7. While hnRNP U stabilizes RNA stem-loop conformations that maintain exon7 skipping, hnRNP L disrupts these RNA elements to facilitate recruitment of the essential splicing factor U2AF2, thereby promoting exon7 inclusion. Our data represent a paradigm for the control of splice site selection by differential RBP binding and modulation of pre-mRNA structure.

Structural basis of mRNA maturation: Time to put it together

In eukaryotes, the expression of genetic information begins in the cell nucleus with precursor messenger RNA (pre-mRNA) transcription and processing into mature mRNA. The mRNA is subsequently recognized and packaged by proteins into an mRNA ribonucleoprotein complex (mRNP) and exported to the cytoplasm for translation. Each of the nuclear mRNA maturation steps is carried out by a dedicated molecular machine. Here, we highlight recent structural and mechanistic insights into how these machines function, including the capping enzyme, the spliceosome, the 3'-end processing machinery, and the transcription-export complex. While we increasingly understand individual steps of nuclear gene expression, many questions remain. For example, we are only beginning to reveal how mature mRNAs are recognized and packaged for nuclear export and how mRNA maturation events are coupled to transcription and to each other. Advances in the preparation of recombinant and endogenous protein-nucleic acid complexes, cryo-electron microscopy, and machine learning promise exciting insights into the mechanisms of nuclear gene expression and its spatial organization.

Cooperative engagement and subsequent selective displacement of SR proteins define the pre-mRNA 3D structural scaffold for early spliceosome assembly

We recently reported that serine–arginine-rich (SR) protein-mediated pre-mRNA structural remodeling generates a pre-mRNA 3D structural scaffold that is stably recognized by the early spliceosomal components. However, the intermediate steps between the free pre-mRNA and the assembled early spliceosome are not yet characterized. By probing the early spliceosomal complexes in vitro and RNA-protein interactions in vivo, we show that the SR proteins bind the pre-mRNAs cooperatively generating a substrate that recruits U1 snRNP and U2AF65 in a splice signal-independent manner. Excess U1 snRNP selectively displaces some of the SR protein molecules from the pre-mRNA generating the substrate for splice signal-specific, sequential recognition by U1 snRNP, U2AF65 and U2AF35. Our work thus identifies a novel function of U1 snRNP in mammalian splicing substrate definition, explains the need for excess U1 snRNP compared to other U snRNPs in vivo, demonstrates how excess SR proteins could inhibit splicing, and provides a conceptual basis to examine if this mechanism of splicing substrate definition is employed by other splicing regulatory proteins.

Regulation of alternative splicing of PaFT and PaFDL1, the FT and FD homologs in Platanus acerifolia

Alternative splicing (AS) selects different alternative splice sites and produces a variety of transcripts with different exon/intron combinations, which may result in multiple protein isoforms. The splicing signals include cis-elements and RNA structures; however, the mechanisms of AS regulation in plants have yet to be elucidated. Previous studies have shown that in Platanus acerifolia, the FLOWERING LOCUS T (FT) homolog PaFT has a unique and complex AS pattern, in which most of the splice forms of PaFT involve the first and/or second intron, and the FD homolog PaFDL1 produces two transcripts via AS, whereas the other FT homolog PaFTL is not regulated by AS. In this study, the regulatory mechanism of the AS of PaFT was demonstrated to be conserved in different plant species. To define the distribution of the AS regulatory signals, the intron-swap, site-directed mutagenesis of alternative splice sites, and deletion experiment were performed. For the PaFT gene, all the signals that regulate the AS of the first intron were located within this intron, while the usage of the first alternative splice site in the second intron was determined by the first intron. Meanwhile, the AS of PaFDL1 might be co-regulated by exons and the first intron. Additionally, the first alternative splice site and adjacent region in PaFT intron 1 might contain cis-elements and/or RNA structures that affect the use of the other sites. This study had provided a deeper insight into the distribution of AS signals in plants, namely the AS signals of different splice sites might exist in the intron where the sites were present, and might also be distributed in exons or other introns.

Dynamics and functional roles of splicing factor autoregulation

Non-core spliceosome components are essential, conserved regulators of alternative splicing. They provide concentration-dependent control of diverse pre-mRNAs. Many splicing factors direct unproductive splicing of their own pre-mRNAs through negative autoregulation. However, the impact of such feedback loops on splicing dynamics at the single-cell level remains unclear. Here, we developed a system to quantitatively analyze negative autoregulatory splicing dynamics by splicing factor SRSF1 in response to perturbations in single HEK293 cells. We show that negative autoregulatory splicing provides critical functions for gene regulation, establishing a ceiling of SRSF1 protein concentration, reducing cell-cell heterogeneity in SRSF1 levels, and buffering variation in transcription. Most important, it adapts SRSF1 splicing activity to variations in demand from other pre-mRNA substrates. A minimal mathematical model of autoregulatory splicing explains these experimentally observed features and provides values for effective biochemical parameters. These results reveal the unique functional roles that splicing negative autoregulation plays in homeostatically regulating transcriptional programs.

Predicting functional riboSNitches in the context of alternative splicing

RNAs are the major regulators of gene expression, and their secondary structures play crucial roles at different levels. RiboSNitches are disease-associated SNPs that cause changes in the pre-mRNA secondary structural ensemble. Several riboSNitches have been detected in the 5' and 3' untranslated regions and lncRNA. Although cases of secondary structural elements playing a regulatory role in alternative splicing are known, regions specific to splicing events, such as splice junctions have not received much attention. We tested splice-site mutations for their efficiency in disrupting the secondary structure and hypothesized that these could play a crucial role in alternative splicing. Multiple riboSNitch prediction methods were applied to obtain overlapping results that are potentially more reliable. Putative riboSNitches were identified from aberrant 5' and 3' splice site mutations, cancer-causing somatic mutations, and genes that harbor the regulatory RNA secondary structural elements. Our workflow for predicting riboSNitches associated with alternative splicing is novel and paves the way for subsequent experimental validation.

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

Splicing is highly regulated and is modulated by numerous factors. Quantitative predictions for how a mutation will affect precursor mRNA (pre-mRNA) structure and downstream function are particularly challenging. Here, we use a novel chemical probing strategy to visualize endogenous precursor and mature MAPT mRNA structures in cells. We used these data to estimate Boltzmann suboptimal structural ensembles, which were then analyzed to predict consequences of mutations on pre-mRNA structure. Further analysis of recent cryo-EM structures of the spliceosome at different stages of the splicing cycle revealed that the footprint of the B^act complex with pre-mRNA best predicted alternative splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene, achieving 74% accuracy. We further developed a β-regression weighting framework that incorporates splice site strength, RNA structure, and exonic/intronic splicing regulatory elements capable of predicting, with 90% accuracy, the effects of 47 known and 6 newly discovered mutations on inclusion of exon 10 of MAPT. This combined experimental and computational framework represents a path forward for accurate prediction of splicing-related disease-causing variants.

Secondary structures in RNA synthesis, splicing and translation

Even though the functional role of mRNA molecules is primarily decided by the nucleotide sequence, several properties are determined by secondary structure conformations. Examples of secondary structures include long range interactions, hairpins, R-loops and G-quadruplexes and they are formed through interactions of non-adjacent nucleotides. Here, we discuss advances in our understanding of how secondary structures can impact RNA synthesis, splicing, translation and mRNA half-life. During RNA synthesis, secondary structures determine RNA polymerase II (RNAPII) speed, thereby influencing splicing. Splicing is also determined by RNA binding proteins and their binding rates are modulated by secondary structures. For the initiation of translation, secondary structures can control the choice of translation start site. Here, we highlight the mechanisms by which secondary structures modulate these processes, discuss advances in technologies to detect and study them systematically, and consider the roles of RNA secondary structures in disease.

Identification of RBMX as a splicing regulator in Parkinsonian mimetic induced alternative splicing of α-synuclein

α-Synuclein (α-syn) plays a precipitating role in Parkinson's disease (PD) due to its tendency to form oligomers and fibrils. The presence of smaller isoforms of α-syn was widely noticed in the affected brain regions of PD patients. 112-synuclein (112-syn) which lacks exon-5, possess enhanced aggregation propensity and forms intracellular inclusions. However, the factors responsible for the skipping of exon-5 are not completely understood. In this context, we aimed to identity the cis & trans-acting elements governing alternative splicing (AS) events by the Parkinsonian agent (MPP⁺) using minigene constructs. Minigene-I and -II were constructed by pruning the intron-4 and -5 regions respectively without altering the branch point adenosine to preserve splicing machinery. Also, chimeric minigenes were engineered by replacing either 5' (Mini-III) or 3' (Mini-IV) flanking intronic regions of exon-5 with other intronic regions (intron-3 and -2) that are not responsive to MPP⁺ induced splicing. While all the above minigenes exhibited MPP⁺-induced skipping of exon-5, Minigene-III did not generate the spliced product indicating that the 5' flanking intronic region (316 bp) of exon-5 possess cis-acting elements responsible for oxidant-induced alternative splicing. RNA-Binding Protein Database (RBDP) analysis revealed the presence of four putative RNA binding proteins (RBPs), namely, RBMX, MBNL1, KHDRBS3 and SFRS1 that may bind to the 316 bp region of intron-4and their expression was substantially reduced following MPP⁺ treatment. Further, overexpression of RBMX mitigated MPP⁺-induced generation of 112-syn and also reduced intracellular α-syn aggregates. Overall, our study identified the pivotal role of the splicing regulator, RBMX, in the pathophysiology of PD.

Micro-dissection and integration of long and short reads to create a robust catalog of kidney compartment-specific isoforms

Studying isoform expression at the microscopic level has always been a challenging task. A classical example is kidney, where glomerular and tubulo-interstitial compartments carry out drastically different physiological functions and thus presumably their isoform expression also differs. We aim at developing an experimental and computational pipeline for identifying isoforms at microscopic structure-level. We microdissected glomerular and tubulo-interstitial compartments from healthy human kidney tissues from two cohorts. The two compartments were separately sequenced with the PacBio RS II platform. These transcripts were then validated using transcripts of the same samples by the traditional Illumina RNA-Seq protocol, distinct Illumina RNA-Seq short reads from European Renal cDNA Bank (ERCB) samples, and annotated GENCODE transcript list, thus identifying novel transcripts. We identified 14,739 and 14,259 annotated transcripts, and 17,268 and 13,118 potentially novel transcripts in the glomerular and tubulo-interstitial compartments, respectively. Of note, relying solely on either short or long reads would have resulted in many erroneous identifications. We identified distinct pathways involved in glomerular and tubulo-interstitial compartments at the isoform level, creating an important experimental and computational resource for the kidney research community.

Modulation of RNA Splicing by Oligonucleotides: Mechanisms of Action and Therapeutic Implications

Dysregulation of RNA splicing causes many diseases and disorders. Several therapeutic approaches have been developed to correct aberrant alternative splicing events for the treatment of cancers and hereditary diseases, including gene therapy and redirecting splicing, using small molecules or splice switching oligonucleotides (SSO). Significant advances in the chemistry and pharmacology of nucleic acid have led to the development of clinically approved SSO drugs for the treatment of spinal muscular dystrophy and Duchenne muscular dystrophy (DMD). In this review, we discuss the mechanisms of SSO action with emphasis on "less common" approaches to modulate alternative splicing, including bipartite and bifunctional SSO, oligonucleotide decoys for splice factors and SSO-mediated mRNA degradation via AS-NMD and NGD pathways. We briefly discuss the current progress and future perspectives of SSO therapy for rare and ultrarare diseases.

A simple guide to de novo transcriptome assembly and annotation

A transcriptome constructed from short-read RNA sequencing (RNA-seq) is an easily attainable proxy catalog of protein-coding genes when genome assembly is unnecessary, expensive or difficult. In the absence of a sequenced genome to guide the reconstruction process, the transcriptome must be assembled de novo using only the information available in the RNA-seq reads. Subsequently, the sequences must be annotated in order to identify sequence-intrinsic and evolutionary features in them (for example, protein-coding regions). Although straightforward at first glance, de novo transcriptome assembly and annotation can quickly prove to be challenging undertakings. In addition to familiarizing themselves with the conceptual and technical intricacies of the tasks at hand and the numerous pre- and post-processing steps involved, those interested must also grapple with an overwhelmingly large choice of tools. The lack of standardized workflows, fast pace of development of new tools and techniques and paucity of authoritative literature have served to exacerbate the difficulty of the task even further. Here, we present a comprehensive overview of de novo transcriptome assembly and annotation. We discuss the procedures involved, including pre- and post-processing steps, and present a compendium of corresponding tools.

Genome-wide identification and functional analysis of circRNAs in Trichophyton rubrum conidial and mycelial stages

Background

Circular RNAs (circRNAs) are a group of noncoding RNAs that participate in gene expression regulation in various pathways. The essential roles of circRNAs have been revealed in many species. However, knowledge of circRNAs in fungi is still not comprehensive.

Results

Trichophyton rubrum (T. rubrum) is considered a model organism of human pathogenic filamentous fungi and dermatophytes. In this study, we performed a genome-wide investigation of circRNAs in T. rubrum based on high-throughput sequencing and ultimately identified 4254 circRNAs. Most of these circRNAs were specific to the conidial or mycelial stage, revealing a developmental stage-specific expression pattern. In addition, 940 circRNAs were significantly differentially expressed between the conidial and mycelial stages. PCR experiments conducted on seven randomly selected differentially expressed (DE-) circRNAs confirmed the circularized structures and relative expression levels of these circRNAs. Based on their genome locations, most circRNAs originated from intergenic regions, unlike those in plants and animals. Furthermore, we constructed circRNA-miRNA-mRNA regulatory networks that included 661 DE-circRNAs targeting 140 miRNAs and further regulating 2753 mRNAs. The relative expression levels of two randomly selected circRNA-miRNA-mRNA axes were investigated by qRT-PCR, and the competing endogenous RNA (ceRNA) network theory was validated. Functional enrichment analysis of the target genes suggested that they were significantly involved in posttranscriptional processes and protein synthesis as well as some small-molecule metabolism processes. CircRNAs are relatively more conserved in closely related dermatophytes but rarely conserved in distantly related species. Tru_circ07138_001 is a highly conserved circRNA that was conserved in all ten dermatophytes analyzed in our study and three distantly related species. Its host gene TERG_07138 was also highly conserved in two of these distantly related species Gallus gallus and Caenorhabditis elegans. The specific role of this circRNA deserves further exploration.

Conclusions

Our study is the first to provide a global profile of circRNAs in T. rubrum as well as dermatophytes. These results could serve as valuable resources for research on circRNA regulatory mechanisms in fungi and reveal new insights for further investigation of the physical characteristics of these significant human fungal pathogens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08184-y.

Alternative Splicing Mechanisms Underlying Opioid-Induced Hyperalgesia

Prolonged use of opioids can cause opioid-induced hyperalgesia (OIH). The impact of alternative splicing on OIH remains partially characterized. A study of the absolute and relative modes of action of alternative splicing further the understanding of the molecular mechanisms underlying OIH. Differential absolute and relative isoform profiles were detected in the trigeminal ganglia and nucleus accumbens of mice presenting OIH behaviors elicited by chronic morphine administration relative to control mice. Genes that participate in glutamatergic synapse (e.g., Grip1, Grin1, Wnk3), myelin protein processes (e.g., Mbp, Mpz), and axon guidance presented absolute and relative splicing associated with OIH. Splicing of genes in the gonadotropin-releasing hormone receptor pathway was detected in the nucleus accumbens while splicing in the vascular endothelial growth factor, endogenous cannabinoid signaling, circadian clock system, and metabotropic glutamate receptor pathways was detected in the trigeminal ganglia. A notable finding was the prevalence of alternatively spliced transcription factors and regulators (e.g., Ciart, Ablim2, Pbx1, Arntl2) in the trigeminal ganglia. Insights into the nociceptive and antinociceptive modulatory action of Hnrnpk were gained. The results from our study highlight the impact of alternative splicing and transcriptional regulators on OIH and expose the need for isoform-level research to advance the understanding of morphine-associated hyperalgesia.

SON drives oncogenic RNA splicing in glioblastoma by regulating PTBP1/PTBP2 switching and RBFOX2 activity

While dysregulation of RNA splicing has been recognized as an emerging target for cancer therapy, the functional significance of RNA splicing and individual splicing factors in brain tumors is poorly understood. Here, we identify SON as a master regulator that activates PTBP1-mediated oncogenic splicing while suppressing RBFOX2-mediated non-oncogenic neuronal splicing in glioblastoma multiforme (GBM). SON is overexpressed in GBM patients and SON knockdown causes failure in intron removal from the PTBP1 transcript, resulting in PTBP1 downregulation and inhibition of its downstream oncogenic splicing. Furthermore, SON forms a complex with hnRNP A2B1 and antagonizes RBFOX2, which leads to skipping of RBFOX2-targeted cassette exons, including the PTBP2 neuronal exon. SON knockdown inhibits proliferation and clonogenicity of GBM cells in vitro and significantly suppresses tumor growth in orthotopic xenografts in vivo. Collectively, our study reveals that SON-mediated RNA splicing is a GBM vulnerability, implicating SON as a potential therapeutic target in brain tumors.

Discovery of a pre-mRNA structural scaffold as a contributor to the mammalian splicing code

The specific recognition of splice signals at or near exon-intron junctions is not explained by their weak conservation and instead is postulated to require a multitude of features embedded in the pre-mRNA strand. We explored the possibility of 3D structural scaffold of AdML-a model pre-mRNA substrate-guiding early spliceosomal components to the splice signal sequences. We find that mutations in the non-cognate splice signal sequences impede recruitment of early spliceosomal components due to disruption of the global structure of the pre-mRNA. We further find that the pre-mRNA segments potentially interacting with the early spliceosomal component U1 snRNP are distributed across the intron, that there is a spatial proximity of 5' and 3' splice sites within the pre-mRNA scaffold, and that an interplay exists between the structural scaffold and splicing regulatory elements in recruiting early spliceosomal components. These results suggest that early spliceosomal components can recognize a 3D structural scaffold beyond the short splice signal sequences, and that in our model pre-mRNA, this scaffold is formed across the intron involving the major splice signals. This provides a conceptual basis to analyze the contribution of recognizable 3D structural scaffolds to the splicing code across the mammalian transcriptome.

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures

RNA structures play a fundamental role in nearly every aspect of cellular physiology and pathology. Gaining insights into the functions of RNA molecules requires accurate predictions of RNA secondary structures. However, the existing thermodynamic folding models remain less accurate than desired, even when chemical probing data, such as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivities, are used as restraints. Unlike most SHAPE-directed algorithms that only consider SHAPE restraints for base pairing, we extract two-dimensional structural features encoded in SHAPE data and establish robust relationships between characteristic SHAPE patterns and loop motifs of various types (hairpin, internal, and bulge) and lengths (2-11 nucleotides). Such characteristic SHAPE patterns are closely related to the sugar pucker conformations of loop residues. Based on these patterns, we propose a computational method, SHAPELoop, which refines the predicted results of the existing methods, thereby further improving their prediction accuracy. In addition, SHAPELoop can provide information about local or global structural rearrangements (including pseudoknots) and help researchers to easily test their hypothesized secondary structures.

RNA architecture influences plant biology

The majority of the genome is transcribed to RNA in living organisms. RNA transcripts can form astonishing arrays of secondary and tertiary structures via Watson-Crick, Hoogsteen, or wobble base pairing. In vivo, RNA folding is not a simple thermodynamic event of minimizing free energy. Instead, the process is constrained by transcription, RNA-binding proteins, steric factors, and the microenvironment. RNA secondary structure (RSS) plays myriad roles in numerous biological processes, such as RNA processing, stability, transportation, and translation in prokaryotes and eukaryotes. Emerging evidence has also implicated RSS in RNA trafficking, liquid-liquid phase separation, and plant responses to environmental variations such as temperature and salinity. At molecular level, RSS is correlated with splicing, polyadenylation, protein synthesis, and miRNA biogenesis and functions. In this review, we summarize newly reported methods for probing RSS in vivo and functions and mechanisms of RSS in plant physiology.

A prometastatic splicing program regulated by SNRPA1 interactions with structured RNA elements

Aberrant alternative splicing is a hallmark of cancer, yet the underlying regulatory programs that control this process remain largely unknown. Here, we report a systematic effort to decipher the RNA structural code that shapes pathological splicing during breast cancer metastasis. We discovered a previously unknown structural splicing enhancer that is enriched near cassette exons with increased inclusion in highly metastatic cells. We show that the spliceosomal protein small nuclear ribonucleoprotein polypeptide A' (SNRPA1) interacts with these enhancers to promote cassette exon inclusion. This interaction enhances metastatic lung colonization and cancer cell invasion, in part through SNRPA1-mediated regulation of PLEC alternative splicing, which can be counteracted by splicing modulating morpholinos. Our findings establish a noncanonical regulatory role for SNRPA1 as a prometastatic splicing enhancer in breast cancer.

Unannotated splicing regulatory elements in deep intron space

Deep intron space harbors a diverse array of splicing regulatory elements that cooperate with better-known exon-proximal elements to enforce proper tissue-specific and development-specific pre-mRNA processing. Many deep intron elements have been highly conserved through vertebrate evolution, yet remain poorly annotated in the human genome. Recursive splicing exons (RS-exons) and intraexons promote noncanonical, multistep resplicing pathways in long introns, involving transient intermediate structures that are greatly underrepresented in RNA-seq datasets. Decoy splice sites and decoy exons act at a distance to inhibit splicing catalysis at annotated splice sites, with functional consequences such as exon skipping and intron retention. RNA:RNA bridges can juxtapose distant sequences within or across introns to activate deep intron splicing enhancers and silencers, to loop out exons to be skipped, or to select one member of a mutually exclusive set of exons. Similarly, protein bridges mediated by interactions among transcript-bound RNA binding proteins (RBPs) can modulate splicing outcomes. Experimental disruption of deep intron elements serving any of these functions can abrogate normal splicing, strongly suggesting that natural mutations of deep intron elements can do likewise to cause human disease. Understanding noncanonical splicing pathways and discovering deep intron regulatory signals, many of which map hundreds to many thousands of nucleotides from annotated splice junctions, is of great academic interest for basic scientists studying alternative splicing mechanisms. Hopefully, this knowledge coupled with increased analysis of deep intron sequences will also have important medical applications, as better interpretation of deep intron mutations may reveal new disease mechanisms and suggest new therapies. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing.

Conserved long-range base pairings are associated with pre-mRNA processing of human genes

The ability of nucleic acids to form double-stranded structures is essential for all living systems on Earth. Current knowledge on functional RNA structures is focused on locally-occurring base pairs. However, crosslinking and proximity ligation experiments demonstrated that long-range RNA structures are highly abundant. Here, we present the most complete to-date catalog of conserved complementary regions (PCCRs) in human protein-coding genes. PCCRs tend to occur within introns, suppress intervening exons, and obstruct cryptic and inactive splice sites. Double-stranded structure of PCCRs is supported by decreased icSHAPE nucleotide accessibility, high abundance of RNA editing sites, and frequent occurrence of forked eCLIP peaks. Introns with PCCRs show a distinct splicing pattern in response to RNAPII slowdown suggesting that splicing is widely affected by co-transcriptional RNA folding. The enrichment of 3'-ends within PCCRs raises the intriguing hypothesis that coupling between RNA folding and splicing could mediate co-transcriptional suppression of premature pre-mRNA cleavage and polyadenylation.

U7 snRNA: A tool for gene therapy

Most U-rich small nuclear ribonucleoproteins (snRNPs) are complexes that mediate the splicing of pre-mRNAs. U7 snRNP is an exception in that it is not involved in splicing but is a key factor in the unique 3' end processing of replication-dependent histone mRNAs. However, by introducing controlled changes in the U7 snRNA histone binding sequence and in the Sm motif, it can be used as an effective tool for gene therapy. The modified U7 snRNP (U7 Sm OPT) is thus not involved in the processing of replication-dependent histone pre-mRNA but targets splicing by inducing efficient skipping or inclusion of selected exons. U7 Sm OPT is of therapeutic importance in diseases that are an outcome of splicing defects, such as myotonic dystrophy, Duchenne muscular dystrophy, amyotrophic lateral sclerosis, β-thalassemia, HIV-1 infection and spinal muscular atrophy. The benefits of using U7 Sm OPT for gene therapy are its compact size, ability to accumulate in the nucleus without causing any toxic effects in the cells, and no immunoreactivity. The risk of transgene misregulation by using U7 Sm OPT is also low because it is involved in correcting the expression of an endogenous gene controlled by its own regulatory elements. Altogether, using U7 Sm OPT as a tool in gene therapy can ensure lifelong treatment, whereas an oligonucleotide or other drug/compound would require repeated administration. It would thus be strategic to harness these unique properties of U7 snRNP and deploy it as a tool in gene therapy.

Transcription Regulation Through Nascent RNA Folding

Folding of RNA into secondary structures through intramolecular base pairing determines an RNA's three-dimensional architecture and associated function. Simple RNA structures like stem loops can provide specialized functions independent of coding capacity, such as protein binding, regulation of RNA processing and stability, stimulation or inhibition of translation. RNA catalysis is dependent on tertiary structures found in the ribosome, tRNAs and group I and II introns. While the extent to which non-coding RNAs contribute to cellular maintenance is generally appreciated, the fact that both non-coding and coding RNA can assume relevant structural states has only recently gained attention. In particular, the co-transcriptional folding of nascent RNA of all classes has the potential to regulate co-transcriptional processing, RNP (ribonucleoprotein particle) formation, and transcription itself. Riboswitches are established examples of co-transcriptionally folded coding RNAs that directly regulate transcription, mainly in prokaryotes. Here we discuss recent studies in both prokaryotes and eukaryotes showing that structure formation may carry a more widespread regulatory logic during RNA synthesis. Local structures forming close to the catalytic center of RNA polymerases have the potential to regulate transcription by reducing backtracking. In addition, stem loops or more complex structures may alter co-transcriptional RNA processing or its efficiency. Several examples of functional structures have been identified to date, and this review provides an overview of physiologically distinct processes where co-transcriptionally folded RNA plays a role. Experimental approaches such as single-molecule FRET and in vivo structural probing to further advance our insight into the significance of co-transcriptional structure formation are discussed.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

The oncogenic and tumor suppressive roles of RNA-binding proteins in human cancers

Posttranscriptional regulation is a mechanism for the cells to control gene regulation at the RNA level. In this process, RNA-binding proteins (RBPs) play central roles and orchestrate the function of RNA molecules in multiple steps. Accumulating evidence has shown that the aberrant regulation of RBPs makes  contributions to the initiation and progression of tumorigenesis via numerous mechanisms such as genetic changes, epigenetic alterations, and noncoding RNA-mediated regulations. In this article, we review the effects caused by RBPs and their functional diversity in the malignant transformation of cancer cells that occurs through the involvement of these proteins in various stages of RNA regulation including alternative splicing, stability, polyadenylation, localization, and translation. Besides this, we review the various interactions between RBPs and other crucial posttranscriptional regulators such as microRNAs and long noncoding RNAs in the pathogenesis of cancer. Finally, we discuss the potential approaches for targeting RBPs in human cancers.

The Chloroplast Trans-Splicing RNA-Protein Supercomplex from the Green Alga Chlamydomonas reinhardtii

In eukaryotes, RNA trans-splicing is a significant RNA modification process for the end-to-end ligation of exons from separately transcribed primary transcripts to generate mature mRNA. So far, three different categories of RNA trans-splicing have been found in organisms within a diverse range. Here, we review trans-splicing of discontinuous group II introns, which occurs in chloroplasts and mitochondria of lower eukaryotes and plants. We discuss the origin of intronic sequences and the evolutionary relationship between chloroplast ribonucleoprotein complexes and the nuclear spliceosome. Finally, we focus on the ribonucleoprotein supercomplex involved in trans-splicing of chloroplast group II introns from the green alga Chlamydomonas reinhardtii. This complex has been well characterized genetically and biochemically, resulting in a detailed picture of the chloroplast ribonucleoprotein supercomplex. This information contributes substantially to our understanding of the function of RNA-processing machineries and might provide a blueprint for other splicing complexes involved in trans- as well as cis-splicing of organellar intron RNAs.

Intronic tRNAs of mitochondrial origin regulate constitutive and alternative splicing

BACKGROUND

The presence of nuclear mitochondrial DNA (numtDNA) has been reported within several nuclear genomes. Next to mitochondrial protein-coding genes, numtDNA sequences also encode for mitochondrial tRNA genes. However, the biological roles of numtDNA remain elusive.

RESULTS

Employing in silico analysis, we identify 281 mitochondrial tRNA homologs in the human genome, which we term nimtRNAs (nuclear intronic mitochondrial-derived tRNAs), being contained within introns of 76 nuclear host genes. Despite base changes in nimtRNAs when compared to their mtRNA homologs, a canonical tRNA cloverleaf structure is maintained. To address potential functions of intronic nimtRNAs, we insert them into introns of constitutive and alternative splicing reporters and demonstrate that nimtRNAs promote pre-mRNA splicing, dependent on the number and positioning of nimtRNA genes and splice site recognition efficiency. A mutational analysis reveals that the nimtRNA cloverleaf structure is required for the observed splicing increase. Utilizing a CRISPR/Cas9 approach, we show that a partial deletion of a single endogenous nimtRNALys within intron 28 of the PPFIBP1 gene decreases inclusion of the downstream-located exon 29 of the PPFIBP1 mRNA. By employing a pull-down approach followed by mass spectrometry, a 3'-splice site-associated protein network is identified, including KHDRBS1, which we show directly interacts with nimtRNATyr by an electrophoretic mobility shift assay.

CONCLUSIONS

We propose that nimtRNAs, along with associated protein factors, can act as a novel class of intronic splicing regulatory elements in the human genome by participating in the regulation of splicing.

Allosteric regulation of U1 snRNP by splicing regulatory proteins controls spliceosomal assembly

Alternative splicing is responsible for much of the transcriptomic and proteomic diversity observed in eukaryotes and involves combinatorial regulation by many cis-acting elements and trans-acting factors. SR and hnRNP splicing regulatory proteins often have opposing effects on splicing efficiency depending on where they bind the pre-mRNA relative to the splice site. Position-dependent splicing repression occurs at spliceosomal E-complex, suggesting that U1 snRNP binds but cannot facilitate higher order spliceosomal assembly. To test the hypothesis that the structure of U1 snRNA changes during activation or repression, we developed a method to structure-probe native U1 snRNP in enriched conformations that mimic activated or repressed spliceosomal E-complexes. While the core of U1 snRNA is highly structured, the 5' end of U1 snRNA shows different SHAPE reactivities and psoralen crosslinking efficiencies depending on where splicing regulatory elements are located relative to the 5' splice site. A motif within the 5' splice site binding region of U1 snRNA is more reactive toward SHAPE electrophiles when repressors are bound, suggesting U1 snRNA is bound, but less base-paired. These observations demonstrate that splicing regulators modulate splice site selection allosterically.

The mechanism and detection of alternative splicing events in circular RNAs

Circular RNAs (circRNAs) are considered as functional biomolecules with tissue/development-specific expression patterns. Generally, a single gene may generate multiple circRNA variants by alternative splicing, which contain different combinations of exons and/or introns. Due to the low abundance of circRNAs as well as overlapped with their linear counterparts, circRNA enrichment protocol is needed prior to sequencing. Compared with numerous algorithms, which use back-splicing reads for detection and functional characterization of circRNAs, original bioinformatic analyzing tools have been developed to large-scale determination of full-length circRNAs and accurate quantification. This review provides insights into the complexity of circRNA biogenesis and surveys the recent progresses in the experimental and bioinformatic methodologies that focus on accurately full-length circRNAs identification.

The role of hypoxia on Alzheimer's disease-related APP and Tau mRNA formation

The removal of introns from mRNA precursors (pre-mRNAs) is an essential step in eukaryotic gene expression. The splicing machinery heavily contributes to biological complexity and especially to the ability of cells to adapt to altered cellular conditions. Hypoxia also plays a key role in the pathophysiology of many diseases, including Alzheimer's disease (AD). In the presented study, we have examined the influence of cellular hypoxia on mRNA splice variant formation from Alzheimer's disease-related Tau and APP genes in brain cells. We have shown that the hypoxic microenvironment influenced the formation of Tau mRNA splice variants, but had no effect on APP mRNA splice variant formation. Additionally, our presented results indicate that splicing factor SRSF1 but not SRSF5 alters the formation of Tau cellular mRNA splice variants in hypoxic cells. Obtained results have also shown that hypoxic brain cells possess enhanced CLK1-4 kinase mRNA levels. This study underlines that cellular hypoxia can influence disease development through changing pre-mRNA splicing.

RNA structures in alternative splicing and back-splicing

Alternative splicing greatly expands the transcriptomic and proteomic diversities related to physiological and developmental processes in higher eukaryotes. Splicing of long noncoding RNAs, and back- and trans- splicing further expanded the regulatory repertoire of alternative splicing. RNA structures were shown to play an important role in regulating alternative splicing and back-splicing. Application of novel sequencing technologies made it possible to identify genome-wide RNA structures and interaction networks, which might provide new insights into RNA splicing regulation in vitro to in vivo. The emerging transcription-folding-splicing paradigm is changing our understanding of RNA alternative splicing regulation. Here, we review the insights into the roles and mechanisms of RNA structures in alternative splicing and back-splicing, as well as how disruption of these structures affects alternative splicing and then leads to human diseases. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

A transient α-helix in the N-terminal RNA recognition motif of polypyrimidine tract binding protein senses RNA secondary structure

The polypyrimidine tract binding protein (PTB) is a multi-domain protein involved in alternative splicing, mRNA localization, stabilization, polyadenylation and translation initiation from internal ribosome entry sites (IRES). In this latter process, PTB promotes viral translation by interacting extensively with complex structured regions in the 5′-untranslated regions of viral RNAs at pyrimidine-rich targets located in single strand and hairpin regions. To better understand how PTB recognizes structured elements in RNA targets, we solved the solution structure of the N-terminal RNA recognition motif (RRM) in complex with an RNA hairpin embedding the loop sequence UCUUU, which is frequently found in IRESs of the picornovirus family. Surprisingly, a new three-turn α3 helix C-terminal to the RRM, folds upon binding the RNA hairpin. Although α3 does not mediate any contacts to the RNA, it acts as a sensor of RNA secondary structure, suggesting a role for RRM1 in detecting pyrimidine tracts in the context of structured RNA. Moreover, the degree of helix formation depends on the RNA loop sequence. Finally, we show that the α3 helix region, which is highly conserved in vertebrates, is crucial for PTB function in enhancing Encephalomyocarditis virus IRES activity.

Congenital posterior cervical spine malformation due to biallelic c.240-4T>G RIPPLY2 variant: A discrete entity

The clinical and radiological spectrum of spondylocostal dysostosis syndromes encompasses distinctive costo-vertebral anomalies. RIPPLY2 biallelic pathogenic variants were described in two distinct cervical spine malformation syndromes: Klippel-Feil syndrome and posterior cervical spine malformation. RIPPLY2 is involved in the determination of rostro-caudal polarity and somite patterning during development. To date, only four cases have been reported. The current report aims at further delineating the posterior malformation in three new patients. Three patients from two unrelated families underwent clinical and radiological examination through X-ray, 3D computed tomography and brain magnetic resonance imaging. After informed consent was obtained, family-based whole exome sequencing (WES) was performed. Complex vertebral segmentation defects in the cervico-thoracic spine were observed in all patients. WES led to the identification of the homozygous splicing variant c.240-4T>G in all subjects. This variant is predicted to result in aberrant splicing of Exon 4. The current report highlights a subtype of cervical spine malformation with major atlo-axoidal malformation compromising spinal cord integrity. This distinctive mutation-specific pattern of malformation differs from Klippel-Feil syndrome and broadens the current classification, defining a sub-type of RIPPLY2-related skeletal disorder. Of note, the phenotype of one patient overlaps with oculo-auriculo-vertebral spectrum disorder.

Accessing a New Dimension in TP53 Biology: Multiplex Long Amplicon Digital PCR to Specifically Detect and Quantitate Individual TP53 Transcripts

TP53, the most commonly-mutated gene in cancer, undergoes complex alternative splicing. Different TP53 transcripts play different biological roles, both in normal function and in the progression of diseases such as cancer. The study of TP53’s alternative RNA splice forms and their use as clinical biomarkers has been hampered by limited specificity and quantitative accuracy of current methods. TP53 RNA splice variants differ at both 5’ and 3’ ends, but because they have a common central region of 618 bp, the individual TP53 transcripts are impossible to specifically detect and precisely quantitate using standard PCR-based methods or short-read RNA sequencing. Therefore, we devised multiplex probe-based long amplicon droplet digital PCR (ddPCR) assays, which for the first time allow precise end-to-end quantitation of the seven major TP53 transcripts, with amplicons ranging from 0.85 to 1.85 kb. Multiple modifications to standard ddPCR assay procedures were required to enable specific co-amplification of these long transcripts and to overcome issues with secondary structure. Using these assays, we show that several TP53 transcripts are co-expressed in breast cancers, and illustrate the potential for this method to identify novel TP53 transcripts in tumour cells. This capability will facilitate a new level of biological and clinical understanding of the alternatively-spliced TP53 isoforms.