Splicing of Nascent RNA Coincides with Intron Exit from RNA Polymerase II

Co-transcriptional splicing efficiencies differ within genes and between cell types

Pre-mRNA splicing is carried out by the spliceosome and involves splice site recognition, removal of introns, and ligation of exons. Components of the spliceosome have been shown to interact with the elongating RNA polymerase II (RNAPII) which is thought to allow splicing to occur concurrently with transcription. However, little is known about the regulation and efficiency of co-transcriptional splicing in human cells. In this study, we used Bru-seq and BruChase-seq to determine the co-transcriptional splicing efficiencies of 17,000 introns expressed across 6 human cell lines. We found that less than half of all introns across these 6 cell lines were co-transcriptionally spliced. Splicing efficiencies for individual introns showed variations across cell lines, suggesting that splicing may be regulated in a cell-type specific manner. Moreover, the splicing efficiency of introns varied within genes. The efficiency of co-transcriptional splicing did not correlate with gene length, intron position, splice site strengths, or the intron/neighboring exons GC content. However, we identified binding signals from multiple RNA binding proteins (RBPs) that correlated with splicing efficiency, including core spliceosomal machinery components-such as SF3B4, U2AF1 and U2AF2 showing higher binding signals in poorly spliced introns. In addition, multiple RBPs, such as BUD13, PUM1 and SND1, showed preferential binding in exons that flank introns with high splicing efficiencies. The nascent RNA splicing patterns presented here across multiple cell types add to our understanding of the complexity in RNA splicing, wherein RNA-binding proteins may play important roles in determining splicing outcomes in a cell type- and intron-specific manner.

Dynamic imaging of nascent RNA reveals general principles of transcription dynamics and stochastic splice site selection

The activities of RNA polymerase and the spliceosome are responsible for the heterogeneity in the abundance and isoform composition of mRNA in human cells. However, the dynamics of these megadalton enzymatic complexes working in concert on endogenous genes have not been described. Here, we establish a quasi-genome-scale platform for observing synthesis and processing kinetics of single nascent RNA molecules in real time. We find that all observed genes show transcriptional bursting. We also observe large kinetic variation in intron removal for single introns in single cells, which is inconsistent with deterministic splice site selection. Transcriptome-wide footprinting of the U2AF complex, nascent RNA profiling, long-read sequencing, and lariat sequencing further reveal widespread stochastic recursive splicing within introns. We propose and validate a unified theoretical model to explain the general features of transcription and pervasive stochastic splice site selection.

Alternative Splicing Dynamics of the Hypothalamus-Pituitary-Ovary Axis During Pubertal Transition in Gilts

The timing of puberty in mammals marks the point at which reproduction becomes possible. Abnormalities in the timing of puberty may exert a series of negative effects on subsequent health outcomes. Alternative splicing (AS) has not only emerged as a significant factor in the transcription of genes but it is also reported to play a role in the timing of puberty. However, to date, the changes and dynamics of AS during the onset of puberty is extremely seldom explored. In the present study, we used gilts as a research model to investigated the dynamics of AS and differentially expressed AS (DEAS) events within the hypothalamus-pituitary-ovary (HPO) axis across pre-, in-, and post-puberty. We detected 3,390, 6,098, and 9,085 DEAS events in the hypothalamus, pituitary, and ovary when compared across pre-, in-, and post-pubertal stages, respectively. Within the entire HPO axis, we also identified 22,889, 22,857, and 21,055 DEAS events in the pre-, in-, and post-pubertal stages, respectively. Further analysis revealed that the differentially spliced genes (DSGs) associated with staged DEAS events were likely to be enriched in the oxytocin signaling pathway, thyroid hormone signaling pathway, GnRH signaling pathway, and oocyte meiosis signaling pathway. The DSGs associated with DEAS events across the entire HPO axis were enriched in endocytosis signaling pathway, the MAPK signaling pathway, and the Rap1 signaling pathway. Moreover. the ASs of TAC1, TACR3, CYP19A1, ESR1, ESRRA, and FSHR were likely to regulate the functions of the certain HPO tissues during the onset of puberty. Collectively, the AS dynamics and DEAS events were comprehensively profiled in hypothalamus, pituitary, and ovary across the pre-, in-, and post-pubertal stages in pigs. These findings may enhance our knowledge of how puberty is regulated by AS and shed new light on the molecular mechanisms underlying the timing of puberty in mammals.

Gene Architecture Facilitates Intron-Mediated Enhancement of Transcription

Introns impact several vital aspects of eukaryotic organisms like proteomic plasticity, genomic stability, stress response and gene expression. A role for introns in the regulation of gene expression at the level of transcription has been known for more than thirty years. The molecular basis underlying the phenomenon, however, is still not entirely clear. An important clue came from studies performed in budding yeast that indicate that the presence of an intron within a gene results in formation of a multi-looped gene architecture. When looping is defective, these interactions are abolished, and there is no enhancement of transcription despite normal splicing. In this review, we highlight several potential mechanisms through which looping interactions may enhance transcription. The promoter-5′ splice site interaction can facilitate initiation of transcription, the terminator-3′ splice site interaction can enable efficient termination of transcription, while the promoter-terminator interaction can enhance promoter directionality and expedite reinitiation of transcription. Like yeast, mammalian genes also exhibit an intragenic interaction of the promoter with the gene body, especially exons. Such promoter-exon interactions may be responsible for splicing-dependent transcriptional regulation. Thus, the splicing-facilitated changes in gene architecture may play a critical role in regulation of transcription in yeast as well as in higher eukaryotes.

Long non-coding RNAs and splicing

In this review I focus on the role of splicing in long non-coding RNA (lncRNA) life. First, I summarize differences between the splicing efficiency of protein-coding genes and lncRNAs and discuss why non-coding RNAs are spliced less efficiently. In the second half of the review, I speculate why splice sites are the most conserved sequences in lncRNAs and what additional roles could splicing play in lncRNA metabolism. I discuss the hypothesis that the splicing machinery can, besides its dominant role in intron removal and exon joining, protect cells from undesired transcripts.

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.

POINT technology illuminates the processing of polymerase-associated intact nascent transcripts

Mammalian chromatin is the site of both RNA polymerase II (Pol II) transcription and coupled RNA processing. However, molecular details of such co-transcriptional mechanisms remain obscure, partly because of technical limitations in purifying authentic nascent transcripts. We present a new approach to characterize nascent RNA, called polymerase intact nascent transcript (POINT) technology. This three-pronged methodology maps nascent RNA 5′ ends (POINT-5), establishes the kinetics of co-transcriptional splicing patterns (POINT-nano), and profiles whole transcription units (POINT-seq). In particular, we show by depletion of the nuclear exonuclease Xrn2 that this activity acts selectively on cleaved 5′ P-RNA at polyadenylation sites. Furthermore, POINT-nano reveals that co-transcriptional splicing either occurs immediately after splice site transcription or is delayed until Pol II transcribes downstream sequences. Finally, we connect RNA cleavage and splicing with either premature or full-length transcript termination. We anticipate that POINT technology will afford full dissection of the complexity of co-transcriptional RNA processing.

•

POINT methodology dissects intact nascent RNA processing

•

Specificity of Xrn2 exonuclease in co-transcriptional RNA degradation

•

Splicing suppresses Xrn2-dependent premature termination

•

Different kinetic classes of co-transcriptional splicing in human genes

Revealing nascent RNA processing dynamics with nano-COP

During maturation, eukaryotic precursor RNAs undergo processing events including intron splicing, 3'-end cleavage, and polyadenylation. Here we describe nanopore analysis of co-transcriptional processing (nano-COP), a method for probing the timing and patterns of RNA processing. An extension of native elongating transcript sequencing, which quantifies transcription genome-wide through short-read sequencing of nascent RNA 3' ends, nano-COP uses long-read nascent RNA sequencing to observe global patterns of RNA processing. First, nascent RNA is stringently purified through a combination of 4-thiouridine metabolic labeling and cellular fractionation. In contrast to cDNA or short-read-based approaches relying on reverse transcription or amplification, the sample is sequenced directly through nanopores to reveal the native context of nascent RNA. nano-COP identifies both active transcription sites and splice isoforms of single RNA molecules during synthesis, providing insight into patterns of intron removal and the physical coupling between transcription and splicing. The nano-COP protocol yields data within 3 d.

2passtools: two-pass alignment using machine-learning-filtered splice junctions increases the accuracy of intron detection in long-read RNA sequencing

Transcription of eukaryotic genomes involves complex alternative processing of RNAs. Sequencing of full-length RNAs using long reads reveals the true complexity of processing. However, the relatively high error rates of long-read sequencing technologies can reduce the accuracy of intron identification. Here we apply alignment metrics and machine-learning-derived sequence information to filter spurious splice junctions from long-read alignments and use the remaining junctions to guide realignment in a two-pass approach. This method, available in the software package 2passtools (https://github.com/bartongroup/2passtools), improves the accuracy of spliced alignment and transcriptome assembly for species both with and without existing high-quality annotations.

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.

Fast and furious: insights of back splicing regulation during nascent RNA synthesis

Alternative splicing of eukaryotic precursor (messenger) RNAs in the nucleus not only increases transcriptomic complexity, but also expands proteomic and functional diversity. In addition to basic types of alternative splicing, recent transcriptome-wide analyses have also suggested other new types of non-canonical splicing, such as back splicing and recursive splicing, and their widespread expression across species Increasing lines of evidence have suggested mechanisms for back splicing, including insights from analyses of nascent RNA sequencing. In this review, we discuss our current understanding of back splicing regulation, and highlight its distinct characteristics in processing during nascent RNA synthesis by taking advantage of metabolic tagging nascent RNA sequencing. Features of recursive splicing are also discussed in the perspective of nascent RNA sequencing.

Co-transcriptional splicing regulates 3' end cleavage during mammalian erythropoiesis

Pre-mRNA processing steps are tightly coordinated with transcription in many organisms. To determine how co-transcriptional splicing is integrated with transcription elongation and 3' end formation in mammalian cells, we performed long-read sequencing of individual nascent RNAs and precision run-on sequencing (PRO-seq) during mouse erythropoiesis. Splicing was not accompanied by transcriptional pausing and was detected when RNA polymerase II (Pol II) was within 75-300 nucleotides of 3' splice sites (3'SSs), often during transcription of the downstream exon. Interestingly, several hundred introns displayed abundant splicing intermediates, suggesting that splicing delays can take place between the two catalytic steps. Overall, splicing efficiencies were correlated among introns within the same transcript, and intron retention was associated with inefficient 3' end cleavage. Remarkably, a thalassemia patient-derived mutation introducing a cryptic 3'SS improved both splicing and 3' end cleavage of individual β-globin transcripts, demonstrating functional coupling between the two co-transcriptional processes as a determinant of productive gene output.

Cotranscriptional and Posttranscriptional Features of the Transcriptome in Soybean Shoot Apex and Leaf

Transcription is the first step of central dogma, in which the genetic information stored in DNA is copied into RNA. In addition to mature RNA sequencing (RNA-seq), high-throughput nascent RNA assays have been established and applied to provide detailed transcriptional information. Here, we present the profiling of nascent RNA from trifoliate leaves and shoot apices of soybean. In combination with nascent RNA (chromatin-bound RNA, CB RNA) and RNA-seq, we found that introns were largely spliced cotranscriptionally. Although alternative splicing (AS) was mainly determined at nascent RNA biogenesis, differential AS between the leaf and shoot apex at the mature RNA level did not correlate well with cotranscriptional differential AS. Overall, RNA abundance was moderately correlated between nascent RNA and mature RNA within each tissue, but the fold changes between the leaf and shoot apex were highly correlated. Thousands of novel transcripts (mainly non-coding RNA) were detected by CB RNA-seq, including the overlap of natural antisense RNA with two important genes controlling soybean reproductive development, FT2a and Dt1. Taken together, we demonstrated the adoption of CB RNA-seq in soybean, which may shed light on gene expression regulation of important agronomic traits in leguminous crops.

Processing of coding and non-coding RNAs in plant development and environmental responses

Precursor RNAs undergo extensive processing to become mature RNAs. RNA transcripts are subjected to 5′ capping, 3′-end processing, splicing, and modification; they also form dynamic secondary structures during co-transcriptional and post-transcriptional processing. Like coding RNAs, non-coding RNAs (ncRNAs) undergo extensive processing. For example, secondary small interfering RNA (siRNA) transcripts undergo RNA processing, followed by further cleavage to become mature siRNAs. Transcriptome studies have revealed roles for co-transcriptional and post-transcriptional RNA processing in the regulation of gene expression and the coordination of plant development and plant–environment interactions. In this review, we present the latest progress on RNA processing in gene expression and discuss phased siRNAs (phasiRNAs), a kind of germ cell-specific secondary small RNA (sRNA), focusing on their functions in plant development and environmental responses.

Preparation of Mammalian Nascent RNA for Long Read Sequencing

Long read sequencing technologies now allow high-quality sequencing of RNAs (or their cDNAs) that are hundreds to thousands of nucleotides long. Long read sequences of nascent RNA provide single-nucleotide-resolution information about co-transcriptional RNA processing events-e.g., splicing, folding, and base modifications. Here, we describe how to isolate nascent RNA from mammalian cells through subcellular fractionation of chromatin-associated RNA, as well as how to deplete poly(A)+ RNA and rRNA, and, finally, how to generate a full-length cDNA library for use on long read sequencing platforms. This approach allows for an understanding of coordinated splicing status across multi-intron transcripts by revealing patterns of splicing or other RNA processing events that cannot be gained from traditional short read RNA sequencing. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Subcellular fractionation Basic Protocol 2: Nascent RNA isolation and adapter ligation Basic Protocol 3: cDNA amplicon preparation.

Widespread Transcriptional Readthrough Caused by Nab2 Depletion Leads to Chimeric Transcripts with Retained Introns

Nascent RNA sequencing has revealed that pre-mRNA splicing can occur shortly after introns emerge from RNA polymerase II (RNA Pol II). Differences in co-transcriptional splicing profiles suggest regulation by cis- and/or trans-acting factors. Here, we use single-molecule intron tracking (SMIT) to identify a cohort of regulators by machine learning in budding yeast. Of these, Nab2 displays reduced co-transcriptional splicing when depleted. Unexpectedly, these splicing defects are attributable to aberrant "intrusive" transcriptional readthrough from upstream genes, as revealed by long-read sequencing. Transcripts that originate from the intron-containing gene's own transcription start site (TSS) are efficiently spliced, indicating no direct role of Nab2 in splicing per se. This work highlights the coupling between transcription, splicing, and 3' end formation in the context of gene organization along chromosomes. We conclude that Nab2 is required for proper 3' end processing, which ensures gene-specific control of co-transcriptional RNA processing.

The Bromodomain Protein 4 Contributes to the Regulation of Alternative Splicing

The bromodomain protein 4 (BRD4) is an atypical kinase and histone acetyl transferase (HAT) that binds to acetylated histones and contributes to chromatin remodeling and early transcriptional elongation. During transcription, BRD4 travels with the elongation complex. Since most alternative splicing events take place co-transcriptionally, we asked if BRD4 plays a role in regulating alternative splicing. We report that distinct patterns of alternative splicing are associated with a conditional deletion of BRD4 during thymocyte differentiation in vivo. Similarly, the depletion of BRD4 in T cell acute lymphoblastic leukemia (T-ALL) cells alters patterns of splicing. Most alternatively spliced events affected by BRD4 are exon skipping. Importantly, BRD4 interacts with components of the splicing machinery, as assessed by both immunoprecipitation (IP) and proximity ligation assays (PLAs), and co-localizes on chromatin with the splicing regulator, FUS. We propose that BRD4 contributes to patterns of alternative splicing through its interaction with the splicing machinery during transcription elongation.

Macrophage development and activation involve coordinated intron retention in key inflammatory regulators

Monocytes and macrophages are essential components of the innate immune system. Herein, we report that intron retention (IR) plays an important role in the development and function of these cells. Using Illumina mRNA sequencing, Nanopore direct cDNA sequencing and proteomics analysis, we identify IR events that affect the expression of key genes/proteins involved in macrophage development and function. We demonstrate that decreased IR in nuclear-detained mRNA is coupled with increased expression of genes encoding regulators of macrophage transcription, phagocytosis and inflammatory signalling, including ID2, IRF7, ENG and LAT. We further show that this dynamic IR program persists during the polarisation of resting macrophages into activated macrophages. In the presence of proinflammatory stimuli, intron-retaining CXCL2 and NFKBIZ transcripts are rapidly spliced, enabling timely expression of these key inflammatory regulators by macrophages. Our study provides novel insights into the molecular factors controlling vital regulators of the innate immune response.

Pathogenic impact of transcript isoform switching in 1,209 cancer samples covering 27 cancer types using an isoform-specific interaction network

Under normal conditions, cells of almost all tissue types express the same predominant canonical transcript isoform at each gene locus. In cancer, however, splicing regulation is often disturbed, leading to cancer-specific switches in the most dominant transcripts (MDT). To address the pathogenic impact of these switches, we have analyzed isoform-specific protein-protein interaction disruptions in 1,209 cancer samples covering 27 different cancer types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) project of the International Cancer Genomics Consortium (ICGC). Our study revealed large variations in the number of cancer-specific MDT (cMDT) with the highest frequency in cancers of female reproductive organs. Interestingly, in contrast to the mutational load, cancers arising from the same primary tissue had a similar number of cMDT. Some cMDT were found in 100% of all samples in a cancer type, making them candidates for diagnostic biomarkers. cMDT tend to be located at densely populated network regions where they disrupted protein interactions in the proximity of pathogenic cancer genes. A gene ontology enrichment analysis showed that these disruptions occurred mostly in protein translation and RNA splicing pathways. Interestingly, samples with mutations in the spliceosomal complex tend to have higher number of cMDT, while other transcript expressions correlated with mutations in non-coding splice-site and promoter regions of their genes. This work demonstrates for the first time the large extent of cancer-specific alterations in alternative splicing for 27 different cancer types. It highlights distinct and common patterns of cMDT and suggests novel pathogenic transcripts and markers that induce large network disruptions in cancers.

Revisiting the window of opportunity for cotranscriptional splicing in budding yeast

We reported previously that, in budding yeast, transcription rate affects both the efficiency and fidelity of pre-mRNA splicing, especially of ribosomal protein transcripts. Here, we report that the majority of ribosomal protein transcripts with non-consensus 5' splice sites are spliced less efficiently when transcription is faster, and more efficiently with slower transcription. These results support the "window of opportunity" model, and we suggest a possible mechanism to explain these findings.

Transcription and splicing: A two-way street

RNA synthesis by RNA polymerase II and RNA processing are closely coupled during the transcription cycle of protein-coding genes. This coupling affords opportunities for quality control and regulation of gene expression and the effects can go in both directions. For example, polymerase speed can affect splice site selection and splicing can increase transcription and affect the chromatin landscape. Here we review the many ways that transcription and splicing influence one another, including how splicing "talks back" to transcription. We will also place the connections between transcription and splicing in the context of other RNA processing events that define the exons that will make up the final mRNA. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Linking transcription, RNA polymerase II elongation and alternative splicing

Gene expression is an intricately regulated process that is at the basis of cell differentiation, the maintenance of cell identity and the cellular responses to environmental changes. Alternative splicing, the process by which multiple functionally distinct transcripts are generated from a single gene, is one of the main mechanisms that contribute to expand the coding capacity of genomes and help explain the level of complexity achieved by higher organisms. Eukaryotic transcription is subject to multiple layers of regulation both intrinsic — such as promoter structure — and dynamic, allowing the cell to respond to internal and external signals. Similarly, alternative splicing choices are affected by all of these aspects, mainly through the regulation of transcription elongation, making it a regulatory knob on a par with the regulation of gene expression levels. This review aims to recapitulate some of the history and stepping-stones that led to the paradigms held today about transcription and splicing regulation, with major focus on transcription elongation and its effect on alternative splicing.

Light in the transcription landscape: chromatin, RNA polymerase II and splicing throughout Arabidopsis thaliana's life cycle

Plants have a high level of developmental plasticity that allows them to respond and adapt to changes in the environment. Among the environmental cues, light controls almost every aspect of A. thaliana's life cycle, including seed maturation, seed germination, seedling de-etiolation and flowering time. Light signals induce massive reprogramming of gene expression, producing changes in RNA polymerase II transcription, alternative splicing, and chromatin state. Since splicing reactions occur mainly while transcription takes place, the regulation of RNAPII transcription has repercussions in the splicing outcomes. This cotranscriptional nature allows a functional coupling between transcription and splicing, in which properties of the splicing reactions are affected by the transcriptional process. Chromatin landscapes influence both transcription and splicing. In this review, we highlight, summarize and discuss recent progress in the field to gain a comprehensive insight on the cross-regulation between chromatin state, RNAPII transcription and splicing decisions in plants, with a special focus on light-triggered responses. We also introduce several examples of transcription and splicing factors that could be acting as coupling factors in plants. Unravelling how these connected regulatory networks operate, can help in the design of better crops with higher productivity and tolerance.

H3K36 Methylation and the Chromodomain Protein Eaf3 Are Required for Proper Cotranscriptional Spliceosome Assembly

In the eukaryotic cell, spliceosomes assemble onto pre-mRNA cotranscriptionally. Spliceosome assembly takes place in the context of the chromatin environment, suggesting that the state of the chromatin may affect splicing. The molecular details and mechanisms through which chromatin affects splicing, however, are still unclear. Here, we show a role for the histone methyltransferase Set2 and its histone modification, H3K36 methylation, in pre-mRNA splicing through high-throughput sequencing. Moreover, the effect of H3K36 methylation on pre-mRNA splicing is mediated through the chromodomain protein Eaf3. We find that Eaf3 is recruited to intron-containing genes and that Eaf3 interacts with the splicing factor Prp45. Eaf3 acts with Prp45 and Prp19 after formation of the precatalytic B complex around the time of splicing activation, thus revealing the step in splicing that is regulated by H3K36 methylation. These studies support a model whereby H3K36 facilitates recruitment of an "adapter protein" to support efficient, constitutive splicing.

Methodologies for Transcript Profiling Using Long-Read Technologies

RNA sequencing using next-generation sequencing technologies (NGS) is currently the standard approach for gene expression profiling, particularly for large-scale high-throughput studies. NGS technologies comprise high throughput, cost efficient short-read RNA-Seq, while emerging single molecule, long-read RNA-Seq technologies have enabled new approaches to study the transcriptome and its function. The emerging single molecule, long-read technologies are currently commercially available by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), while new methodologies based on short-read sequencing approaches are also being developed in order to provide long range single molecule level information-for example, the ones represented by the 10x Genomics linked read methodology. The shift toward long-read sequencing technologies for transcriptome characterization is based on current increases in throughput and decreases in cost, making these attractive for de novo transcriptome assembly, isoform expression quantification, and in-depth RNA species analysis. These types of analyses were challenging with standard short sequencing approaches, due to the complex nature of the transcriptome, which consists of variable lengths of transcripts and multiple alternatively spliced isoforms for most genes, as well as the high sequence similarity of highly abundant species of RNA, such as rRNAs. Here we aim to focus on single molecule level sequencing technologies and single-cell technologies that, combined with perturbation tools, allow the analysis of complete RNA species, whether short or long, at high resolution. In parallel, these tools have opened new ways in understanding gene functions at the tissue, network, and pathway levels, as well as their detailed functional characterization. Analysis of the epi-transcriptome, including RNA methylation and modification and the effects of such modifications on biological systems is now enabled through direct RNA sequencing instead of classical indirect approaches. However, many difficulties and challenges remain, such as methodologies to generate full-length RNA or cDNA libraries from all different species of RNAs, not only poly-A containing transcripts, and the identification of allele-specific transcripts due to current error rates of single molecule technologies, while the bioinformatics analysis on long-read data for accurate identification of 5' and 3' UTRs is still in development.

Post-transcriptional splicing of nascent RNA contributes to widespread intron retention in plants

In eukaryotes, genes are transcribed by RNA polymerase-II (Pol-II) and introns are removed by the spliceosome largely cotranscriptionally^1-3; analysis using long-read sequencing revealed that splicing occurs immediately after Pol-II passes introns in yeast^4,5. Here, we developed a Nanopore-based method to profile chromatin-bound RNA that enables the simultaneous detection of splicing status, Pol-II position and polyadenylation at the genome-wide scale in Arabidopsis. We found that more than half of the introns remain unspliced after Pol-II transcribes 1 kb past the 3' splice site, which is much slower than the rate of splicing reported in yeast^4,5. Many of the full-length chromatin-bound RNA molecules are polyadenylated, yet still contain unspliced introns at specific positions. These introns are nearly absent in the cytoplasm and are resistant to nonsense-mediated decay, suggesting that they are post-transcriptionally spliced before the transcripts are released into the cytoplasm; we therefore termed these introns post-transcriptionally spliced introns (pts introns). Analysis of around 6,500 public RNA-sequencing libraries found that the splicing of pts introns requires the function of splicing-related proteins such as PRMT5 and SKIP, and is also influenced by various environmental signals. The majority of the intron retention events in Arabidopsis are at pts introns, suggesting that chromatin-tethered post-transcriptional splicing is a major contributor to the widespread intron retention that is observed in plants, and could be a mechanism to produce fully spliced functional mRNAs for rapid response.

Global transcriptome analysis reveals circadian control of splicing events in Arabidopsis thaliana

The circadian clock of Arabidopsis thaliana controls many physiological and molecular processes, allowing plants to anticipate daily changes in their environment. However, developing a detailed understanding of how oscillations in mRNA levels are connected to oscillations in co/post-transcriptional processes, such as splicing, has remained a challenge. Here we applied a combined approach using deep transcriptome sequencing and bioinformatics tools to identify novel circadian-regulated genes and splicing events. Using a stringent approach, we identified 300 intron retention, eight exon skipping, 79 alternative 3' splice site usage, 48 alternative 5' splice site usage, and 350 multiple (more than one event type) annotated events under circadian regulation. We also found seven and 721 novel alternative exonic and intronic events. Depletion of the circadian-regulated splicing factor AtSPF30 homologue resulted in the disruption of a subset of clock-controlled splicing events. Altogether, our global circadian RNA-seq coupled with an in silico, event-centred, splicing analysis tool offers a new approach for studying the interplay between the circadian clock and the splicing machinery at a global scale. The identification of many circadian-regulated splicing events broadens our current understanding of the level of control that the circadian clock has over this co/post-transcriptional regulatory layer.

Epigenome-based splicing prediction using a recurrent neural network

Alternative RNA splicing provides an important means to expand metazoan transcriptome diversity. Contrary to what was accepted previously, splicing is now thought to predominantly take place during transcription. Motivated by emerging data showing the physical proximity of the spliceosome to Pol II, we surveyed the effect of epigenetic context on co-transcriptional splicing. In particular, we observed that splicing factors were not necessarily enriched at exon junctions and that most epigenetic signatures had a distinctly asymmetric profile around known splice sites. Given this, we tried to build an interpretable model that mimics the physical layout of splicing regulation where the chromatin context progressively changes as the Pol II moves along the guide DNA. We used a recurrent-neural-network architecture to predict the inclusion of a spliced exon based on adjacent epigenetic signals, and we showed that distinct spatio-temporal features of these signals were key determinants of model outcome, in addition to the actual nucleotide sequence of the guide DNA strand. After the model had been trained and tested (with >80% precision-recall curve metric), we explored the derived weights of the latent factors, finding they highlight the importance of the asymmetric time-direction of chromatin context during transcription.

The changing paradigm of intron retention: regulation, ramifications and recipes

Intron retention (IR) is a form of alternative splicing that has long been neglected in mammalian systems although it has been studied for decades in non-mammalian species such as plants, fungi, insects and viruses. It was generally assumed that mis-splicing, leading to the retention of introns, would have no physiological consequence other than reducing gene expression by nonsense-mediated decay. Relatively recent landmark discoveries have highlighted the pivotal role that IR serves in normal and disease-related human biology. Significant technical hurdles have been overcome, thereby enabling the robust detection and quantification of IR. Still, relatively little is known about the cis- and trans-acting modulators controlling this phenomenon. The fate of an intron to be, or not to be, retained in the mature transcript is the direct result of the influence exerted by numerous intrinsic and extrinsic factors at multiple levels of regulation. These factors have altered current biological paradigms and provided unexpected insights into the transcriptional landscape. In this review, we discuss the regulators of IR and methods to identify them. Our focus is primarily on mammals, however, we broaden the scope to non-mammalian organisms in which IR has been shown to be biologically relevant.

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

Exon Definition Facilitates Reliable Control of Alternative Splicing in the RON Proto-Oncogene

Alternative splicing is a key step in eukaryotic gene expression that allows for the production of multiple transcript and protein isoforms from the same gene. Even though splicing is perturbed in many diseases, we currently lack insights into regulatory mechanisms promoting its precision and efficiency. We analyze high-throughput mutagenesis data obtained for an alternatively spliced exon in the proto-oncogene RON and determine the functional units that control this splicing event. Using mathematical modeling of distinct splicing mechanisms, we show that alternative splicing is based in RON on a so-called "exon definition" mechanism. Here, the recognition of the adjacent exons by the spliceosome is required for removal of an intron. We use our model to analyze the differences between the exon and intron definition scenarios and find that exon definition prevents the accumulation of deleterious, partially spliced retention products during alternative splicing regulation. Furthermore, it modularizes splicing control, as multiple regulatory inputs are integrated into a common net input, irrespective of the location and nature of the corresponding cis-regulatory elements in the pre-messenger RNA. Our analysis suggests that exon definition promotes robust and reliable splicing outcomes in RON splicing.

The Features and Regulation of Co-transcriptional Splicing in Arabidopsis

Precursor mRNA (pre-mRNA) splicing is essential for gene expression in most eukaryotic organisms. Previous studies from mammals, Drosophila, and yeast show that the majority of splicing events occurs co-transcriptionally. In plants, however, the features of co-transcriptional splicing (CTS) and its regulation still remain largely unknown. Here, we used chromatin-bound RNA sequencing to study CTS in Arabidopsis thaliana. We found that CTS is widespread in Arabidopsis seedlings, with a large proportion of alternative splicing events determined co-transcriptionally. CTS efficiency correlated with gene expression level, the chromatin landscape and, most surprisingly, the number of introns and exons of individual genes, but is independent of gene length. In combination with enhanced crosslinking and immunoprecipitation sequencing analysis, we further showed that the hnRNP-like proteins RZ-1B and RZ-1C promote efficient CTS globally through direct binding, frequently to exonic sequences. Notably, this general effect of RZ-1B/1C on splicing promotion is mainly observed at the chromatin level, not at the mRNA level. RZ-1C promotes CTS of multiple-exon genes in association with its binding to regions both proximal and distal to the regulated introns. We propose that RZ-1C promotes efficient CTS of genes with multiple exons through cooperative interactions with many exons, introns, and splicing factors. Our work thus reveals important features of CTS in plants and provides methodologies for the investigation of CTS and RNA-binding proteins in plants.

Alternative Splicing of Circadian Clock Genes Correlates With Temperature in Field-Grown Sugarcane

Alternative Splicing (AS) is a mechanism that generates different mature transcripts from precursor mRNAs (pre-mRNAs) of the same gene. In plants, a wide range of physiological and metabolic events are related to AS, as well as fast responses to changes in temperature. AS is present in around 60% of intron-containing genes in Arabidopsis, 46% in rice, and 38% in maize and it is widespread among the circadian clock genes. Little is known about how AS influences the circadian clock of C4 plants, like commercial sugarcane, a C4 crop with a complex hybrid genome. This work aims to test if the daily dynamics of AS forms of circadian clock genes are regulated by environmental factors, such as temperature, in the field. A systematic search for AS in five sugarcane clock genes, ScLHY, ScPRR37, ScPRR73, ScPRR95, and ScTOC1 using different organs of sugarcane sampled during winter, with 4 months old plants, and during summer, with 9 months old plants, revealed temperature- and organ-dependent expression of at least one alternatively spliced isoform in all genes. Expression of AS isoforms varied according to the season. Our results suggest that AS events in circadian clock genes are correlated with temperature.

Dido3-dependent SFPQ recruitment maintains efficiency in mammalian alternative splicing

Alternative splicing is facilitated by accessory proteins that guide spliceosome subunits to the primary transcript. Many of these splicing factors recognize the RNA polymerase II tail, but SFPQ is a notable exception even though essential for mammalian RNA processing. This study reveals a novel role for Dido3, one of three Dido gene products, in alternative splicing. Binding of the Dido3 amino terminus to histones and to the polymerase jaw domain was previously reported, and here we show interaction between its carboxy terminus and SFPQ. We generated a mutant that eliminates Dido3 but preserves other Dido gene products, mimicking reduced Dido3 levels in myeloid neoplasms. Dido mutation suppressed SFPQ binding to RNA and increased skipping for a large group of exons. Exons bearing recognition sequences for alternative splicing factors were nonetheless included more efficiently. Reduced SFPQ recruitment may thus account for increased skipping of SFPQ-dependent exons, but could also generate a splicing factor surplus that becomes available to competing splice sites. Taken together, our data indicate that Dido3 is an adaptor that controls SFPQ utilization in RNA splicing. Distributing splicing factor recruitment over parallel pathways provides mammals with a simple mechanism to regulate exon usage while maintaining RNA splicing efficiency.

Nascent RNA analyses: tracking transcription and its regulation

The programmes that direct an organism's development and maintenance are encoded in its genome. Decoding of this information begins with regulated transcription of genomic DNA into RNA. Although transcription and its control can be tracked indirectly by measuring stable RNAs, it is only by directly measuring nascent RNAs that the immediate regulatory changes in response to developmental, environmental, disease and metabolic signals are revealed. Multiple complementary methods have been developed to quantitatively track nascent transcription genome-wide at nucleotide resolution, all of which have contributed novel insights into the mechanisms of gene regulation and transcription-coupled RNA processing. Here we critically evaluate the array of strategies used for investigating nascent transcription and discuss the recent conceptual advances they have provided.

Does co-transcriptional regulation of alternative splicing mediate plant stress responses?

Plants display exquisite control over gene expression to elicit appropriate responses under normal and stress conditions. Alternative splicing (AS) of pre-mRNAs, a process that generates two or more transcripts from multi-exon genes, adds another layer of regulation to fine-tune condition-specific gene expression in animals and plants. However, exactly how plants control splice isoform ratios and the timing of this regulation in response to environmental signals remains elusive. In mammals, recent evidence indicate that epigenetic and epitranscriptome changes, such as DNA methylation, chromatin modifications and RNA methylation, regulate RNA polymerase II processivity, co-transcriptional splicing, and stability and translation efficiency of splice isoforms. In plants, the role of epigenetic modifications in regulating transcription rate and mRNA abundance under stress is beginning to emerge. However, the mechanisms by which epigenetic and epitranscriptomic modifications regulate AS and translation efficiency require further research. Dynamic changes in the chromatin landscape in response to stress may provide a scaffold around which gene expression, AS and translation are orchestrated. Finally, we discuss CRISPR/Cas-based strategies for engineering chromatin architecture to manipulate AS patterns (or splice isoforms levels) to obtain insight into the epigenetic regulation of AS.

Nuclear sorting of RNA

The majority of the mammalian genome is transcribed by RNA polymerase II, yielding a vast amount of noncoding RNA (ncRNA) in addition to the standard production of mRNA. The typical nuclear biogenesis of mRNA relies on the tightly controlled coupling of co- and post-transcriptional processing events, which ultimately results in the export of transcripts into the cytoplasm. These processes are subject to surveillance by nuclear RNA decay pathways to prevent the export of aberrant, or otherwise "non-optimal," transcripts. However, unlike mRNA, many long ncRNAs are nuclear retained and those that maintain enduring functions must employ precautions to evade decay. Proper sorting and localization of RNA is therefore an essential activity in eukaryotic cells and the formation of ribonucleoprotein complexes during early stages of RNA synthesis is central to deciding such transcript fate. This review details our current understanding of the pathways and factors that direct RNAs towards a particular destiny and how transcripts combat the adverse conditions of the nucleus. This article is categorized under: RNA Export and Localization > Nuclear Export/Import RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Spt5 modulates cotranscriptional spliceosome assembly in Saccharomyces cerevisiae

There is increasing evidence from yeast to humans that pre-mRNA splicing occurs mainly cotranscriptionally, such that splicing and transcription are functionally coupled. Currently, there is little insight into the contribution of the core transcription elongation machinery to cotranscriptional spliceosome assembly and pre-mRNA splicing. Spt5 is a member of the core transcription elongation machinery and an essential protein, whose absence in budding yeast causes defects in pre-mRNA splicing. To determine how Spt5 affects pre-mRNA splicing, we used the auxin-inducible degron system to conditionally deplete Spt5 in Saccharomyces cerevisiae and assayed effects on cotranscriptional spliceosome assembly and splicing. We show that Spt5 is needed for efficient splicing and for the accumulation of U5 snRNPs at intron-containing genes, and therefore for stable cotranscriptional assembly of spliceosomes. The defect in cotranscriptional spliceosome assembly can explain the relatively mild splicing defect as being a consequence of the failure of cotranscriptional splicing. Coimmunoprecipitation of Spt5 with core spliceosomal proteins and all spliceosomal snRNAs suggests a model whereby Spt5 promotes cotranscriptional pre-mRNA splicing by stabilizing the association of U5 snRNP with spliceosome complexes as they assemble on the nascent transcript. If this phenomenon is conserved in higher eukaryotes, it has the potential to be important for cotranscriptional regulation of alternative splicing.

MicroRNA miR-1002 Enhances NMNAT-Mediated Stress Response by Modulating Alternative Splicing

Understanding endogenous regulation of stress resistance and homeostasis maintenance is critical to developing neuroprotective therapies. Nicotinamide mononucleotide adenylyltransferase (NMNAT) is a conserved essential enzyme that confers extraordinary protection and stress resistance in many neurodegenerative disease models. Drosophila Nmnat is alternatively spliced to two mRNA variants, RA and RB. RB translates to protein isoform PD with robust protective activity and is upregulated upon stress to confer enhanced neuroprotection. The mechanisms regulating the alternative splicing and stress response of NMNAT remain unclear. We have discovered a Drosophila microRNA, dme-miR-1002, which promotes the splicing of NMNAT pre-mRNA to RB by disrupting a pre-mRNA stem-loop structure. NMNAT pre-mRNA is preferentially spliced to RA in basal conditions, whereas miR-1002 enhances NMNAT PD-mediated stress protection by binding via RISC component Argonaute1 to the pre-mRNA, facilitating the splicing switch to RB. These results outline a new process for microRNAs in regulating alternative splicing and modulating stress resistance.

Getting the Entire Message: Progress in Isoform Sequencing

The advent of second-generation sequencing and its application to RNA sequencing have revolutionized the field of genomics by allowing quantification of gene expression, as well as the definition of transcription start/end sites, exons, splice sites and RNA editing sites. However, due to the sequencing of fragments of cDNAs, these methods have not given a reliable picture of complete RNA isoforms. Third-generation sequencing has filled this gap and allows end-to-end sequencing of entire RNA/cDNA molecules. This approach to transcriptomics has been a "niche" technology for a couple of years but now is becoming mainstream with many different applications. Here, we review the background and progress made to date in this rapidly growing field. We start by reviewing the progressive realization that alternative splicing is omnipresent. We then focus on long-noncoding RNA isoforms and the distinct combination patterns of exons in noncoding and coding genes. We consider the implications of the recent technologies of direct RNA sequencing and single-cell isoform RNA sequencing. Finally, we discuss the parameters that define the success of long-read RNA sequencing experiments and strategies commonly used to make the most of such data.

A novel enhancer RNA, Hmrhl, positively regulates its host gene, phkb, in chronic myelogenous leukemia

Noncoding RNAs are increasingly being accredited with key roles in gene regulation during development and disease. Here we report the discovery and characterization of a novel long noncoding RNA, Hmrhl, which shares synteny and partial sequence similarity with the mouse lncRNA, Mrhl. The human homolog, Hmrhl, transcribed from intron 14 of phkb gene, is 5.5 kb in size, expressed in all tissues examined and is associated with chromatin. Analysis of Hmrhl locus using ENCODE database revealed that it exhibits hallmarks of enhancers like the open chromatin configuration, binding of transcription factors, enhancer specific histone signature etc. in the K562 Chronic Myelogenous Leukemia (CML) cells. We compared the expression of Hmrhl in the normal lymphoblast cell line, GM12878, with that of K562 cells and lymphoma samples and show that it is highly upregulated in leukemia as well as several cases of lymphoma. Further, we validated the enhancer properties of Hmrhl locus in K562 cells with the help of ChIP-qPCR and Luciferase assay. Moreover, siRNA mediated down-regulation of Hmrhl in K562 cells leads to a concomitant down regulation of its parent gene, phkb, showing that Hmrhl functions as an enhancer RNA and positively regulates its host gene, phkb, in chronic myelogenous leukemia. This study is significant in view of the fact that a better understanding of mechanism of gene regulation under normal conditions and its perturbation in cancer could in turn help in its therapeutic intervention through molecular medicine/RNA based drug discovery.

Nascent RNA and the Coordination of Splicing with Transcription

At each active protein-encoding gene, nascent RNA is tethered to the DNA axis by elongating RNA polymerase II (Pol II) and is continuously altered by splicing and other processing events during its synthesis. This review discusses the development of three major methods that enable us to track the conversion of precursor messenger RNA (pre-mRNA) to messenger RNA (mRNA) products in vivo: live-cell imaging, metabolic labeling of RNA, and RNA-seq of purified nascent RNA. These approaches are complementary, addressing distinct issues of transcription rates and intron lifetimes alongside spatial information regarding the gene position of Pol II at which spliceosomes act. The findings will be placed in the context of active transcription units, each of which-because of the presence of nascent RNA, Pol II, and features of the chromatin environment-will recruit a potentially gene-specific constellation of RNA binding proteins and processing machineries.

Conserved sequences in the Drosophila mod(mdg4) intron promote poly(A)-independent transcription termination and trans-splicing

Alternative splicing (AS) is a regulatory mechanism of gene expression that greatly expands the coding capacities of genomes by allowing the generation of multiple mRNAs from a single gene. In Drosophila, the mod(mdg4) locus is an extreme example of AS that produces more than 30 different mRNAs via trans-splicing that joins together the common exons and the 3′ variable exons generated from alternative promoters. To map the regions required for trans-splicing, we have developed an assay for measuring trans-splicing events and identified a 73-bp region in the last common intron that is critical for trans-splicing of three pre-mRNAs synthesized from different DNA strands. We have also found that conserved sequences in the distal part of the last common intron induce polyadenylation-independent transcription termination and are enriched by paused RNA polymerase II (RNAP II). These results suggest that all mod(mdg4) mRNAs are formed by joining in trans the 5′ splice site in the last common exon with the 3′ splice site in one of the alternative exons.

Aberrant splicing in B-cell acute lymphoblastic leukemia

Aberrant splicing is a hallmark of leukemias with mutations in splicing factor (SF)-encoding genes. Here we investigated its prevalence in pediatric B-cell acute lymphoblastic leukemias (B-ALL), where SFs are not mutated. By comparing these samples to normal pro-B cells, we found thousands of aberrant local splice variations (LSVs) per sample, with 279 LSVs in 241 genes present in every comparison. These genes were enriched in RNA processing pathways and encoded ∼100 SFs, e.g. hnRNPA1. HNRNPA1 3'UTR was most pervasively mis-spliced, yielding the transcript subject to nonsense-mediated decay. To mimic this event, we knocked it down in B-lymphoblastoid cells and identified 213 hnRNPA1-regulated exon usage events comprising the hnRNPA1 splicing signature in pediatric leukemia. Some of its elements were LSVs in DICER1 and NT5C2, known cancer drivers. We searched for LSVs in other leukemia and lymphoma drivers and discovered 81 LSVs in 41 additional genes. Seventy-seven LSVs out of 81 were confirmed using two large independent B-ALL RNA-seq datasets, and the twenty most common B-ALL drivers, including NT5C2, showed higher prevalence of aberrant splicing than of somatic mutations. Thus, post-transcriptional deregulation of SF can drive widespread changes in B-ALL splicing and likely contributes to disease pathogenesis.

Noncoding AUG circRNAs constitute an abundant and conserved subclass of circles

Circular RNAs (circRNAs) are a subset of noncoding RNAs previously considered as products of missplicing. Now, circRNAs are considered functional molecules, although to date, only few functions have been experimentally validated. Here, based on RNA sequencing from the ENCODE consortium, we identify and characterize a subset of circRNAs, coined AUG circRNAs, encompassing the annotated translational start codon from the protein-coding host genes. AUG circRNAs are more abundantly expressed and conserved than other groups of circRNAs, and they display flanking sequences that suggest an Alu-independent mechanism of biogenesis. The AUG circRNAs contain part of bona fide open reading frame, and in the recent years, several studies have reported cases of circRNA translation. However, using thorough cross-species analysis, extensive ribosome profiling, proteomics analyses, and experimental data on a selected panel of AUG circRNAs, we observe no indications of translation of AUG circRNAs or any other circRNAs. Our data provide a comprehensive classification of circRNAs and, collectively, the data suggest that the AUG circRNAs constitute an abundant subclass of circRNAs produced independently of primate-specific Alu elements.

Circular exonic RNAs: When RNA structure meets topology

Although RNA circularization was first documented in the 1990s, the extent to which it occurs was not known until recent advances in high-throughput sequencing enabled the widespread identification of circular RNAs (circRNAs). Despite this, many aspects of circRNA biogenesis, structure, and function yet remain obscure. This review focuses on circular exonic RNAs, a subclass of circRNAs that are generated through backsplicing. Here, I hypothesize that RNA secondary structure can be the common factor that promotes both exon skipping and spliceosomal RNA circularization, and that backsplicing of double-stranded regions could generate topologically linked circRNA molecules. CircRNAs manifest themselves by the presence of tail-to-head exon junctions, which were previously attributed to post-transcriptional exon permutation and repetition. I revisit these observations and argue that backsplicing does not automatically imply RNA circularization because tail-to-head exon junctions give only local information about transcript architecture and, therefore, they are in principle insufficient to determine globally circular topology. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Global donor and acceptor splicing site kinetics in human cells

RNA splicing is an essential part of eukaryotic gene expression. Although the mechanism of splicing has been extensively studied in vitro, in vivo kinetics for the two-step splicing reaction remain poorly understood. Here, we combine transient transcriptome sequencing (TT-seq) and mathematical modeling to quantify RNA metabolic rates at donor and acceptor splice sites across the human genome. Splicing occurs in the range of minutes and is limited by the speed of RNA polymerase elongation. Splicing kinetics strongly depends on the position and nature of nucleotides flanking splice sites, and on structural interactions between unspliced RNA and small nuclear RNAs in spliceosomal intermediates. Finally, we introduce the 'yield' of splicing as the efficiency of converting unspliced to spliced RNA and show that it is highest for mRNAs and independent of splicing kinetics. These results lead to quantitative models describing how splicing rates and yield are encoded in the human genome.

RNA Polymerase II Phosphorylated on CTD Serine 5 Interacts with the Spliceosome during Co-transcriptional Splicing

The highly intronic nature of protein coding genes in mammals necessitates a co-transcriptional splicing mechanism as revealed by mNET-seq analysis. Immunoprecipitation of MNase-digested chromatin with antibodies against RNA polymerase II (Pol II) shows that active spliceosomes (both snRNA and proteins) are complexed to Pol II S5P CTD during elongation and co-transcriptional splicing. Notably, elongating Pol II-spliceosome complexes form strong interactions with nascent transcripts, resulting in footprints of approximately 60 nucleotides. Also, splicing intermediates formed by cleavage at the 5′ splice site are associated with nearly all spliced exons. These spliceosome-bound intermediates are frequently ligated to upstream exons, implying a sequential, constitutive, and U12-dependent splicing process. Finally, lack of detectable spliced products connected to the Pol II active site in human HeLa or murine lymphoid cells suggests that splicing does not occur immediately following 3′ splice site synthesis. Our results imply that most mammalian splicing requires exon definition for completion.

•

S5P CTD Pol II associates with the catalytic spliceosome

•

Elongating Pol II complexes protect about 60 newly synthesized nucleotides

•

Co-transcriptional splicing associated with dominant 5′ ss intermediate

•

U12-dependent introns are sequentially spliced in association with Pol II