Nascent-seq indicates widespread cotranscriptional pre-mRNA splicing in Drosophila

RNA structure in alternative splicing regulation: from mechanism to therapy

Alternative splicing is a highly intricate process that plays a crucial role in post-transcriptional regulation and significantly expands the functional proteome of a limited number of coding genes in eukaryotes. Its regulation is multifactorial, with RNA structure exerting a significant impact. Aberrant RNA conformations lead to dysregulation of splicing patterns, which directly affects the manifestation of disease symptoms. In this review, the molecular mechanisms of RNA secondary structure-mediated splicing regulation are summarized, with a focus on the complex interplay between aberrant RNA conformations and disease phenotypes resulted from splicing defects. This study also explores additional factors that reshape structural conformations, enriching our understanding of the mechanistic network underlying structure-mediated splicing regulation. In addition, an emphasis has been placed on the clinical role of targeting aberrant splicing corrections in human diseases. The principal mechanisms of action behind this phenomenon are described, followed by a discussion of prospective development strategies and pertinent challenges.

Analysis of a detailed multi-stage model of stochastic gene expression using queueing theory and model reduction

We introduce a biologically detailed, stochastic model of gene expression describing the multiple rate-limiting steps of transcription, nuclear pre-mRNA processing, nuclear mRNA export, cytoplasmic mRNA degradation and translation of mRNA into protein. The processes in sub-cellular compartments are described by an arbitrary number of processing stages, thus accounting for a significantly finer molecular description of gene expression than conventional models such as the telegraph, two-stage and three-stage models of gene expression. We use two distinct tools, queueing theory and model reduction using the slow-scale linear-noise approximation, to derive exact or approximate analytic expressions for the moments or distributions of nuclear mRNA, cytoplasmic mRNA and protein fluctuations, as well as lower bounds for their Fano factors in steady-state conditions. We use these to study the phase diagram of the stochastic model; in particular we derive parametric conditions determining three types of transitions in the properties of mRNA fluctuations: from sub-Poissonian to super-Poissonian noise, from high noise in the nucleus to high noise in the cytoplasm, and from a monotonic increase to a monotonic decrease of the Fano factor with the number of processing stages. In contrast, protein fluctuations are always super-Poissonian and show weak dependence on the number of mRNA processing stages. Our results delineate the region of parameter space where conventional models give qualitatively incorrect results and provide insight into how the number of processing stages, e.g. the number of rate-limiting steps in initiation, splicing and mRNA degradation, shape stochastic gene expression by modulation of molecular memory.

Timing is everything: advances in quantifying splicing kinetics

Splicing is a highly regulated process critical for proper pre-mRNA maturation and the maintenance of a healthy cellular environment. Splicing events are impacted by ongoing transcription, neighboring splicing events, and cis and trans regulatory factors on the respective pre-mRNA transcript. Within this complex regulatory environment, splicing kinetics have the potential to influence splicing outcomes but have historically been challenging to study in vivo. In this review, we highlight recent technological advancements that have enabled measurements of global splicing kinetics and of the variability of splicing kinetics at single introns. We demonstrate how identifying features that are correlated with splicing kinetics has increased our ability to form potential models for how splicing kinetics may be regulated in vivo.

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

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

RNA Pol II-dependent transcription efficiency fine-tunes A-to-I editing levels

A-to-I RNA editing is a widespread epitranscriptomic phenomenon leading to the conversion of adenosines to inosines, which are primarily interpreted as guanosines by cellular machines. Consequently, A-to-I editing can alter splicing or lead to recoding of transcripts. As misregulation of editing can cause a variety of human diseases, A-to-I editing requires tight regulation of the extent of deamination, particularly in protein-coding regions. The bulk of A-to-I editing occurs cotranscriptionally. Thus, we studied A-to-I editing regulation in the context of transcription and pre-mRNA processing. We show that stimulation of transcription impacts editing levels. Activation of the transcription factor MYC leads to an up-regulation of A-to-I editing, particularly in transcripts that are suppressed upon MYC activation. Moreover, low pre-mRNA synthesis rates and low pre-mRNA expression levels support high levels of editing. We also show that editing levels greatly differ between nascent pre-mRNA and mRNA in a cellular system, as well as in mouse tissues. Editing levels can increase or decrease from pre-mRNA to mRNA and can vary across editing targets and across tissues, showing that pre-mRNA processing is an important layer of editing regulation. Several lines of evidence suggest that the differences emerge during pre-mRNA splicing. Moreover, actinomycin D treatment of primary neuronal cells and editing level analysis suggests that regulation of editing levels also depends on transcription.

Truncating the spliceosomal 'rope protein' Prp45 results in Htz1 dependent phenotypes

Spliceosome assembly contributes an important but incompletely understood aspect of splicing regulation. Prp45 is a yeast splicing factor which runs as an extended fold through the spliceosome, and which may be important for bringing its components together. We performed a whole genome analysis of the genetic interaction network of the truncated allele of PRP45 (prp45(1-169)) using synthetic genetic array technology and found chromatin remodellers and modifiers as an enriched category. In agreement with related studies, H2A.Z-encoding HTZ1, and the components of SWR1, INO80, and SAGA complexes represented prominent interactors, with htz1 conferring the strongest growth defect. Because the truncation of Prp45 disproportionately affected low copy number transcripts of intron-containing genes, we prepared strains carrying intronless versions of SRB2, VPS75, or HRB1, the most affected cases with transcription-related function. Intron removal from SRB2, but not from the other genes, partly repaired some but not all the growth phenotypes identified in the genetic screen. The interaction of prp45(1-169) and htz1Δ was detectable even in cells with SRB2 intron deleted (srb2Δi). The less truncated variant, prp45(1-330), had a synthetic growth defect with htz1Δ at 16°C, which also persisted in the srb2Δi background. Moreover, htz1Δ enhanced prp45(1-330) dependent pre-mRNA hyper-accumulation of both high and low efficiency splicers, genes ECM33 and COF1, respectively. We conclude that while the expression defects of low expression intron-containing genes contribute to the genetic interactome of prp45(1-169), the genetic interactions between prp45 and htz1 alleles demonstrate the sensitivity of spliceosome assembly, delayed in prp45(1-169), to the chromatin environment.

Pre-mRNA splicing and its cotranscriptional connections

Transcription of eukaryotic genes by RNA polymerase II (Pol II) yields RNA precursors containing introns that must be spliced out and the flanking exons ligated together. Splicing is catalyzed by a dynamic ribonucleoprotein complex called the spliceosome. Recent evidence has shown that a large fraction of splicing occurs cotranscriptionally as the RNA chain is extruded from Pol II at speeds of up to 5 kb/minute. Splicing is more efficient when it is tethered to the transcription elongation complex, and this linkage permits functional coupling of splicing with transcription. We discuss recent progress that has uncovered a network of connections that link splicing to transcript elongation and other cotranscriptional RNA processing events.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

Insplico: effective computational tool for studying splicing order of adjacent introns genome-wide with short and long RNA-seq reads

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Cotranscriptional RNA processing and modification in plants

The activities of RNA polymerases shape the epigenetic landscape of genomes with profound consequences for genome integrity and gene expression. A fundamental event during the regulation of eukaryotic gene expression is the coordination between transcription and RNA processing. Most primary RNAs mature through various RNA processing and modification events to become fully functional. While pioneering results positioned RNA maturation steps after transcription ends, the coupling between the maturation of diverse RNA species and their transcription is becoming increasingly evident in plants. In this review, we discuss recent advances in our understanding of the crosstalk between RNA Polymerase II, IV, and V transcription and nascent RNA processing of both coding and noncoding RNAs.

3D genome organization around nuclear speckles drives mRNA splicing efficiency

The nucleus is highly organized such that factors involved in transcription and processing of distinct classes of RNA are organized within specific nuclear bodies. One such nuclear body is the nuclear speckle, which is defined by high concentrations of protein and non-coding RNA regulators of pre-mRNA splicing. What functional role, if any, speckles might play in the process of mRNA splicing remains unknown. Here we show that genes localized near nuclear speckles display higher spliceosome concentrations, increased spliceosome binding to their pre-mRNAs, and higher co-transcriptional splicing levels relative to genes that are located farther from nuclear speckles. We show that directed recruitment of a pre-mRNA to nuclear speckles is sufficient to drive increased mRNA splicing levels. Finally, we show that gene organization around nuclear speckles is highly dynamic with differential localization between cell types corresponding to differences in Pol II occupancy. Together, our results integrate the longstanding observations of nuclear speckles with the biochemistry of mRNA splicing and demonstrate a critical role for dynamic 3D spatial organization of genomic DNA in driving spliceosome concentrations and controlling the efficiency of mRNA splicing.

A simple and robust method for isolating and analyzing chromatin-bound RNAs in Arabidopsis

BACKGROUND

Chromatin-bound RNAs are the primary product of transcription that undergo on-chromatin processing such as capping, splicing, and polyadenylation. These processing steps then determine the fate of RNAs. Albeit its vital importance, a simple and robust method for isolating different fractions of chromatin-bound RNAs is missing in plants.

RESULT

Here, we describe our updated method and the associated step-by-step protocol for chromatin-bound RNAs isolation in A. thaliana. The chromatin-bound RNAs isolation is based on the 1 M UREA wash that removes the majority of non-chromatin-associated proteins from the nucleus, as previously developed in mammalian cells. On-demand, the isolated chromatin-bound RNAs can be either used directly for gene-specific analysis or subject to further rRNA removal and also the optional polyadenylated RNA removal, followed by high-throughput sequencing. Detailed protocols for these procedures are also provided. Comparison of sequencing results of chromatin-bound RNAs with and without polyadenylated RNA removal revealed that a small fraction of CB-RNAs is polyadenylated but not yet fully spliced, representing RNA-processing intermediate on-chromatin.

CONCLUSION

This optimized chromatin-bound RNAs purification method is simple and robust and can be used to study transcription and its-coupled RNA processing in plants.

Non-coding antisense transcripts: fine regulation of gene expression in cancer

Natural antisense transcripts (NATs) are coding or non-coding RNA sequences transcribed on the opposite direction from the same genomic locus. NATs are widely distributed throughout the human genome and seem to play crucial roles in physiological and pathological processes, through newly described and targeted mechanisms. NATs represent the intricate complexity of the genome organization and constitute another layer of potential targets in disease. Here, we focus on the interesting and unique role of non-coding NATs in cancer, paying particular attention to those acting as miRNA sponges.

Transcriptome-Wide Detection of Intron/Exon Definition in the Endogenous Pre-mRNA Transcripts of Mammalian Cells and Its Regulation by Depolarization

Pairing of splice sites across an intron or exon is the central point of intron or exon definition in pre-mRNA splicing with the latter mode proposed for most mammalian exons. However, transcriptome-wide pairing within endogenous transcripts has not been examined for the prevalence of each mode in mammalian cells. Here we report such pairings in rat GH3 pituitary cells by measuring the relative abundance of nuclear RNA-Seq reads at the intron start or end (RISE). Interestingly, RISE indexes are positively correlated between 5′ and 3′ splice sites specifically across introns or exons but inversely correlated with the usage of adjacent exons. Moreover, the ratios between the paired indexes were globally modulated by depolarization, which was disruptible by 5-aza-Cytidine. The nucleotide matrices of the RISE-positive splice sites deviate significantly from the rat consensus, and short introns or exons are enriched with the cross-intron or -exon RISE pairs, respectively. Functionally, the RISE-positive genes cluster for basic cellular processes including RNA binding/splicing, or more specifically, hormone production if regulated by depolarization. Together, the RISE analysis identified the transcriptome-wide regulation of either intron or exon definition between weak splice sites of short introns/exons in mammalian cells. The analysis also provides a way to further track the splicing intermediates and intron/exon definition during the dynamic regulation of alternative splicing by extracellular factors.

Comprehensive analysis of the circadian nuclear and cytoplasmic transcriptome in mouse liver

In eukaryotes, RNA is synthesised in the nucleus, spliced, and exported to the cytoplasm where it is translated and finally degraded. Any of these steps could be subject to temporal regulation during the circadian cycle, resulting in daily fluctuations of RNA accumulation and affecting the distribution of transcripts in different subcellular compartments. Our study analysed the nuclear and cytoplasmic, poly(A) and total transcriptomes of mouse livers collected over the course of a day. These data provide a genome-wide temporal inventory of enrichment in subcellular RNA, and revealed specific signatures of splicing, nuclear export and cytoplasmic mRNA stability related to transcript and gene lengths. Combined with a mathematical model describing rhythmic RNA profiles, we could test the rhythmicity of export rates and cytoplasmic degradation rates of approximately 1400 genes. With nuclear export times usually much shorter than cytoplasmic half-lives, we found that nuclear export contributes to the modulation and generation of rhythmic profiles of 10% of the cycling nuclear mRNAs. This study contributes to a better understanding of the dynamic regulation of the transcriptome during the day-night cycle.

circMbl functions in cis and in trans to regulate gene expression and physiology in a tissue-specific fashion

Muscleblind (mbl) is an essential muscle and neuronal splicing regulator. Mbl hosts multiple circular RNAs (circRNAs), including circMbl, which is conserved from flies to humans. Here, we show that mbl-derived circRNAs are key regulators of MBL by cis- and trans-acting mechanisms. By generating fly lines to specifically modulate the levels of all mbl RNA isoforms, including circMbl, we demonstrate that the two major mbl protein isoforms, MBL-O/P and MBL-C, buffer their own levels by producing different types of circRNA isoforms in the eye and fly brain, respectively. Moreover, we show that circMbl has unique functions in trans, as knockdown of circMbl results in specific morphological and physiological phenotypes. In addition, depletion of MBL-C or circMbl results in opposite behavioral phenotypes, showing that they also regulate each other in trans. Together, our results illuminate key aspects of mbl regulation and uncover cis and trans functions of circMbl in vivo.

Circular RNAs in Pregnancy and the Placenta

The emerging field of circular RNAs (circRNAs) has identified their novel roles in the development and function of many cancers and inspired the interest of many researchers. circRNAs are also found throughout the healthy body, as well as in other pathological states, but while research into the function and abundance of circRNAs has progressed, our overall understanding of these molecules remains primitive. Importantly, recent studies are elucidating new roles for circRNAs in pregnancy, particularly in the placenta. Given that many of the genes responsible for circRNA production in cancer are also highly expressed in the placenta, it is likely that the same genes act in the production of circRNAs in the placenta. Furthermore, placental development can be referred to as ‘controlled cancer’, as it shares many key signalling pathways and hallmarks with tumour growth and metastasis. Hence, the roles of circRNAs in this field are important to study with respect to pregnancy success but also may provide novel insights for cancer progression. This review illuminates the known roles of circRNAs in pregnancy and the placenta, as well as demonstrating differential placental expressions of circRNAs between complicated and uncomplicated pregnancies.

Co-transcriptional splicing efficiency is a gene-specific feature that can be regulated by TGFβ

Differential splicing efficiency of specific introns is a mechanism that dramatically increases protein diversity, based on selection of alternative exons for the final mature mRNA. However, it is unclear whether splicing efficiency of introns within the same gene is coordinated and eventually regulated as a mechanism to control mature mRNA levels. Based on nascent chromatin-associated RNA-sequencing data, we now find that co-transcriptional splicing (CTS) efficiency tends to be similar between the different introns of a gene. We establish that two well-differentiated strategies for CTS efficiency exist, at the extremes of a gradient: short genes that produce high levels of pre-mRNA undergo inefficient splicing, while long genes with relatively low levels of pre-mRNA have an efficient splicing. Notably, we observe that genes with efficient CTS display a higher level of mature mRNA relative to their pre-mRNA levels. Further, we show that the TGFβ signal transduction pathway regulates the general CTS efficiency, causing changes in mature mRNA levels. Taken together, our data indicate that CTS efficiency is a gene-specific characteristic that can be regulated to control gene expression.

Co-transcriptional splicing efficiency is a gene-specific characteristic that can be regulated by TGFβ to modulate gene expression.

Longitudinal intronic RNA-Seq analysis of Parkinson's disease patients reveals disease-specific nascent transcription

Transcriptomic studies usually focus on either gene or exon-based annotations, and only limited experiments have reported changes in reads mapping to introns. The analysis of intronic reads allows the detection of nascent transcription that is not influenced by steady-state RNA levels and provides information on actively transcribed genes. Here, we describe substantial intronic transcriptional changes in Parkinson's disease (PD) patients compared to healthy controls (CO) at two different timepoints; at the time of diagnosis (BL) and three years later (V08). We used blood RNA-Seq data from the Parkinson's Progression Markers Initiative (PPMI) cohort and identified significantly changed transcription of intronic reads only in PD patients during this follow-up period. In CO subjects, only nine transcripts demonstrated differentially expressed introns between visits. However, in PD patients, 4873 transcripts had differentially expressed introns at visit V08 compared to BL, many of them in genes previously associated with neurodegenerative diseases, such as LRRK2, C9orf72, LGALS3, KANSL1AS1, and ALS2. In addition, at the time of diagnosis (BL visit), we identified 836 transcripts (e.g. SNCA, DNAJC19, PRRG4) and at visit V08, 2184 transcripts (e.g. PINK1, GBA, ALS2, PLEKHM1) with differential intronic expression specific to PD patients. In contrast, reads mapping to exonic regions demonstrated little variation indicating highly specific changes only in intronic transcription. Our study demonstrated that PD is characterized by substantial changes in the nascent transcription, and description of these changes could help to understand the molecular pathology underpinning this disease.

CMTr cap-adjacent 2'-O-ribose mRNA methyltransferases are required for reward learning and mRNA localization to synapses

Cap-adjacent nucleotides of animal, protist and viral mRNAs can be O-methylated at the 2‘ position of the ribose (cOMe). The functions of cOMe in animals, however, remain largely unknown. Here we show that the two cap methyltransferases (CMTr1 and CMTr2) of Drosophila can methylate the ribose of the first nucleotide in mRNA. Double-mutant flies lack cOMe but are viable. Consistent with prominent neuronal expression, they have a reward learning defect that can be rescued by conditional expression in mushroom body neurons before training. Among CMTr targets are cell adhesion and signaling molecules. Many are relevant for learning, and are also targets of Fragile X Mental Retardation Protein (FMRP). Like FMRP, cOMe is required for localization of untranslated mRNAs to synapses and enhances binding of the cap binding complex in the nucleus. Hence, our study reveals a mechanism to co-transcriptionally prime mRNAs by cOMe for localized protein synthesis at synapses.

The two cap methyltransferases (CMTrs) redundantly methylate riboses of first cap adjacent nucleotides in messenger RNAs in Drosophila. Here, CMTrs are required for reward learning and localization of untranslated messenger RNAs to synapses.

Pseudouridine synthases modify human pre-mRNA co-transcriptionally and affect pre-mRNA processing

Pseudouridine is a modified nucleotide that is prevalent in human mRNAs and is dynamically regulated. Here, we investigate when in their life cycle mRNAs become pseudouridylated to illuminate the potential regulatory functions of endogenous mRNA pseudouridylation. Using single-nucleotide resolution pseudouridine profiling on chromatin-associated RNA from human cells, we identified pseudouridines in nascent pre-mRNA at locations associated with alternatively spliced regions, enriched near splice sites, and overlapping hundreds of binding sites for RNA-binding proteins. In vitro splicing assays establish a direct effect of individual endogenous pre-mRNA pseudouridines on splicing efficiency. We validate hundreds of pre-mRNA sites as direct targets of distinct pseudouridine synthases and show that PUS1, PUS7, and RPUSD4-three pre-mRNA-modifying pseudouridine synthases with tissue-specific expression-control widespread changes in alternative pre-mRNA splicing and 3' end processing. Our results establish a vast potential for cotranscriptional pre-mRNA pseudouridylation to regulate human gene expression via alternative pre-mRNA processing.

Transcription and splicing dynamics during early Drosophila development

Widespread cotranscriptional splicing has been demonstrated from yeast to human. However, most studies to date addressing the kinetics of splicing relative to transcription used either Saccharomyces cerevisiae or metazoan cultured cell lines. Here, we adapted native elongating transcript sequencing technology (NET-seq) to measure cotranscriptional splicing dynamics during the early developmental stages of Drosophila melanogaster embryos. Our results reveal the position of RNA polymerase II (Pol II) when both canonical and recursive splicing occur. We found heterogeneity in splicing dynamics, with some RNAs spliced immediately after intron transcription, whereas for other transcripts no splicing was observed over the first 100 nt of the downstream exon. Introns that show splicing completion before Pol II has reached the end of the downstream exon are necessarily intron-defined. We studied the splicing dynamics of both nascent pre-mRNAs transcribed in the early embryo, which have few and short introns, as well as pre-mRNAs transcribed later in embryonic development, which contain multiple long introns. As expected, we found a relationship between the proportion of spliced reads and intron size. However, intron definition was observed at all intron sizes. We further observed that genes transcribed in the early embryo tend to be isolated in the genome whereas genes transcribed later are often overlapped by a neighboring convergent gene. In isolated genes, transcription termination occurred soon after the polyadenylation site, while in overlapped genes, Pol II persisted associated with the DNA template after cleavage and polyadenylation of the nascent transcript. Taken together, our data unravel novel dynamic features of Pol II transcription and splicing in the developing Drosophila embryo.

The Hox transcription factor Ultrabithorax binds RNA and regulates co-transcriptional splicing through an interplay with RNA polymerase II

Transcription factors (TFs) play a pivotal role in cell fate decision by coordinating gene expression programs. Although most TFs act at the DNA layer, few TFs bind RNA and modulate splicing. Yet, the mechanistic cues underlying TFs activity in splicing remain elusive. Focusing on the Drosophila Hox TF Ultrabithorax (Ubx), our work shed light on a novel layer of Ubx function at the RNA level. Transcriptome and genome-wide binding profiles in embryonic mesoderm and Drosophila cells indicate that Ubx regulates mRNA expression and splicing to promote distinct outcomes in defined cellular contexts. Our results demonstrate a new RNA-binding ability of Ubx. We find that the N51 amino acid of the DNA-binding Homeodomain is non-essential for RNA interaction in vitro, but is required for RNA interaction in vivo and Ubx splicing activity. Moreover, mutation of the N51 amino acid weakens the interaction between Ubx and active RNA Polymerase II (Pol II). Our results reveal that Ubx regulates elongation-coupled splicing, which could be coordinated by a dynamic interplay with active Pol II on chromatin. Overall, our work uncovered a novel role of the Hox TFs at the mRNA regulatory layer. This could be an essential function for other classes of TFs to control cell diversity.

Using Drosophila to uncover molecular and physiological functions of circRNAs

Circular RNAs (circRNAs) are a class of covalently closed RNA molecules generated by backsplicing. circRNAs are expressed in a tissue-specific manner, accumulate with age in neural tissues, and are highly stable. In many cases, circRNAs are generated at the expense of a linear transcript as back-splicing competes with linear splicing. Some circRNAs regulate gene expression in cis, and some circRNAs can be translated into protein. The advent of deep sequencing and new bioinformatic tools has allowed detection of thousands of circRNAs in eukaryotes. Studying the functions of circRNAs is done using a combination of molecular and genetic methods. The unique genetic tools that can be used in studies of Drosophila melanogaster are ideal for deciphering the functions of circRNAs in vivo. These tools include the GAL4-UAS system, which can be used to manipulate the levels of circRNAs with exquisite temporal and spatial control, and genetic interaction screening, which could be used to identify pathways regulated by circRNAs. Research performed in Drosophila has revealed circRNAs production mechanisms, details of their translation, and their physiological functions. Due to their short lifecycle and the existence of excellent neurodegeneration models, Drosophila can also be used to study the role of circRNAs in aging and age-related disorders. Here, we review molecular and genetic tools and methods for detecting, manipulating, and studying circRNAs in Drosophila.

CircRNAs and their regulatory roles in cancers

Circular RNAs (circRNAs), a novel type of non-coding RNAs (ncRNAs), have a covalently closed circular structure resulting from pre-mRNA back splicing via spliceosome and ribozymes. They can be classified differently in accordance with different criteria. As circRNAs are abundant, conserved, and stable, they can be used as diagnostic markers in various diseases and targets to develop new therapies. There are various functions of circRNAs, including sponge for miR/proteins, role of scaffolds, templates for translation, and regulators of mRNA translation and stability. Without m7G cap and poly-A tail, circRNAs can still be degraded in several ways, including RNase L, Ago-dependent, and Ago-independent degradation. Increasing evidence indicates that circRNAs can be modified by N-6 methylation (m6A) in many aspects such as biogenesis, nuclear export, translation, and degradation. In addition, they have been proved to play a regulatory role in the progression of various cancers. Recently, methods of detecting circRNAs with high sensitivity and specificity have also been reported. This review presents a detailed overview of circRNAs regarding biogenesis, biomarker, functions, degradation, and dynamic modification as well as their regulatory roles in various cancers. It’s particularly summarized in detail in the biogenesis of circRNAs, regulation of circRNAs by m6A modification and mechanisms by which circRNAs affect tumor progression respectively. Moreover, existing circRNA detection methods and their characteristics are also mentioned.

RNA polymerase II speed: a key player in controlling and adapting transcriptome composition

RNA polymerase II (RNA Pol II) speed or elongation rate, i.e., the number of nucleotides synthesized per unit of time, is a major determinant of transcriptome composition. It controls co-transcriptional processes such as splicing, polyadenylation, and transcription termination, thus regulating the production of alternative splice variants, circular RNAs, alternatively polyadenylated transcripts, or read-through transcripts. RNA Pol II speed itself is regulated in response to intra- and extra-cellular stimuli and can in turn affect the transcriptome composition in response to these stimuli. Evidence points to a potentially important role of transcriptome composition modification through RNA Pol II speed regulation for adaptation of cells to a changing environment, thus pointing to a function of RNA Pol II speed regulation in cellular physiology. Analyzing RNA Pol II speed dynamics may therefore be central to fully understand the regulation of physiological processes, such as the development of multicellular organisms. Recent findings also raise the possibility that RNA Pol II speed deregulation can be detrimental and participate in disease progression. Here, we review initial and current approaches to measure RNA Pol II speed, as well as providing an overview of the factors controlling speed and the co-transcriptional processes which are affected. Finally, we discuss the role of RNA Pol II speed regulation in cell physiology.

Challenges and Opportunities in Understanding Genetics of Fungal Diseases: Towards a Functional Genomics Approach

Infectious diseases are a leading cause of morbidity and mortality worldwide, and human pathogens have long been recognized as one of the main sources of evolutionary pressure, resulting in a high variable genetic background in immune-related genes. The study of the genetic contribution to infectious diseases has undergone tremendous advances over the last decades. Here, focusing on genetic predisposition to fungal diseases, we provide an overview of the available approaches for studying human genetic susceptibility to infections, reviewing current methodological and practical limitations. We describe how the classical methods available, such as family-based studies and candidate gene studies, have contributed to the discovery of crucial susceptibility factors for fungal infections. We will also discuss the contribution of novel unbiased approaches to the field, highlighting their success but also their limitations for the fungal immunology field. Finally, we show how a systems genomics approach can overcome those limitations and can lead to efficient prioritization and identification of genes and pathways with a critical role in susceptibility to fungal diseases. This knowledge will help to stratify at-risk patient groups and, subsequently, develop early appropriate prophylactic and treatment strategies.

SPLICE-q: a Python tool for genome-wide quantification of splicing efficiency

BACKGROUND

Introns are generally removed from primary transcripts to form mature RNA molecules in a post-transcriptional process called splicing. An efficient splicing of primary transcripts is an essential step in gene expression and its misregulation is related to numerous human diseases. Thus, to better understand the dynamics of this process and the perturbations that might be caused by aberrant transcript processing it is important to quantify splicing efficiency.

RESULTS

Here, we introduce SPLICE-q, a fast and user-friendly Python tool for genome-wide SPLICing Efficiency quantification. It supports studies focusing on the implications of splicing efficiency in transcript processing dynamics. SPLICE-q uses aligned reads from strand-specific RNA-seq to quantify splicing efficiency for each intron individually and allows the user to select different levels of restrictiveness concerning the introns' overlap with other genomic elements such as exons of other genes. We applied SPLICE-q to globally assess the dynamics of intron excision in yeast and human nascent RNA-seq. We also show its application using total RNA-seq from a patient-matched prostate cancer sample.

CONCLUSIONS

Our analyses illustrate that SPLICE-q is suitable to detect a progressive increase of splicing efficiency throughout a time course of nascent RNA-seq and it might be useful when it comes to understanding cancer progression beyond mere gene expression levels. SPLICE-q is available at: https://github.com/vrmelo/SPLICE-q.

Modifying splice site usage with ModCon: Maintaining the genetic code while changing the underlying mRNP code

Codon degeneracy of amino acid sequences permits an additional “mRNP code” layer underlying the genetic code that is related to RNA processing. In pre-mRNA splicing, splice site usage is determined by both intrinsic strength and sequence context providing RNA binding sites for splicing regulatory proteins. In this study, we systematically examined modification of splicing regulatory properties in the neighborhood of a GT site, i.e. potential splice site, without altering the encoded amino acids.

We quantified the splicing regulatory properties of the neighborhood around a potential splice site by its Splice Site HEXplorer Weight (SSHW) based on the HEXplorer score algorithm. To systematically modify GT site neighborhoods, either minimizing or maximizing their SSHW, we designed the novel stochastic optimization algorithm ModCon that applies a genetic algorithm with stochastic crossover, insertion and random mutation elements supplemented by a heuristic sliding window approach.

To assess the achievable range in SSHW in human splice donors without altering the encoded amino acids, we applied ModCon to a set of 1000 randomly selected Ensembl annotated human splice donor sites, achieving substantial and accurate changes in SSHW. Using ModCon optimization, we successfully switched splice donor usage in a splice site competition reporter containing coding sequences from FANCA, FANCB or BRCA2, while retaining their amino acid coding information.

The ModCon algorithm and its R package implementation can assist in reporter design by either introducing novel splice sites, silencing accidental, undesired splice sites, and by generally modifying the entire mRNP code while maintaining the genetic code.

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.

The Exon Junction Complex and intron removal prevent re-splicing of mRNA

Accurate splice site selection is critical for fruitful gene expression. Recently, the mammalian EJC was shown to repress competing, cryptic, splice sites (SS). However, the evolutionary generality of this remains unclear. Here, we demonstrate the Drosophila EJC suppresses hundreds of functional cryptic SS, even though most bear weak splicing motifs and are seemingly incompetent. Mechanistically, the EJC directly conceals cryptic splicing elements by virtue of its position-specific recruitment, preventing aberrant SS definition. Unexpectedly, we discover the EJC inhibits scores of regenerated 5' and 3' recursive SS on segments that have already undergone splicing, and that loss of EJC regulation triggers faulty resplicing of mRNA. An important corollary is that certain intronless cDNA constructs yield unanticipated, truncated transcripts generated by resplicing. We conclude the EJC has conserved roles to defend transcriptome fidelity by (1) repressing illegitimate splice sites on pre-mRNAs, and (2) preventing inadvertent activation of such sites on spliced segments.

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.

Circular sisRNA identification and characterisation

Stable Intronic Sequence RNA (sisRNA) is a relatively new class of non-coding RNA. Found in many organisms, these sisRNA produced from their host genes are generally involved in regulatory roles, controlling gene expression at multiple levels through active involvement in regulatory feedback loops. Large scale identification of sisRNA via genome-wide RNA sequencing has been difficult, largely in part due to its low abundance. Done on its own, RNA sequencing often yields a large mass of information that is ironically uninformative; the potential sisRNA reads being masked by other highly abundant RNA species like ribosomal RNA and messenger RNA. In this review, we present a practical workflow for the enrichment of circular sisRNA through the use of transcriptionally quiescent systems, rRNA-depletion, and RNase R treatment prior to deep sequencing. This workflow allows circular sisRNA to be reliably detected. We also present various methods to experimentally validate the circularity and stability of the circular sisRNA identified, as well as a few methods for further functional characterisation.

ELAV/Hu RNA binding proteins determine multiple programs of neural alternative splicing

ELAV/Hu factors are conserved RNA binding proteins (RBPs) that play diverse roles in mRNA processing and regulation. The founding member, Drosophila Elav, was recognized as a vital neural factor 35 years ago. Nevertheless, little was known about its impacts on the transcriptome, and potential functional overlap with its paralogs. Building on our recent findings that neural-specific lengthened 3' UTR isoforms are co-determined by ELAV/Hu factors, we address their impacts on splicing. While only a few splicing targets of Drosophila are known, ectopic expression of each of the three family members (Elav, Fne and Rbp9) alters hundreds of cassette exon and alternative last exon (ALE) splicing choices. Reciprocally, double mutants of elav/fne, but not elav alone, exhibit opposite effects on both classes of regulated mRNA processing events in larval CNS. While manipulation of Drosophila ELAV/Hu RBPs induces both exon skipping and inclusion, characteristic ELAV/Hu motifs are enriched only within introns flanking exons that are suppressed by ELAV/Hu factors. Moreover, the roles of ELAV/Hu factors in global promotion of distal ALE splicing are mechanistically linked to terminal 3' UTR extensions in neurons, since both processes involve bypass of proximal polyadenylation signals linked to ELAV/Hu motifs downstream of cleavage sites. We corroborate the direct action of Elav in diverse modes of mRNA processing using RRM-dependent Elav-CLIP data from S2 cells. Finally, we provide evidence for conservation in mammalian neurons, which undergo broad programs of distal ALE and APA lengthening, linked to ELAV/Hu motifs downstream of regulated polyadenylation sites. Overall, ELAV/Hu RBPs orchestrate multiple broad programs of neuronal mRNA processing and isoform diversification in Drosophila and mammalian neurons.

Intron exon boundary junctions in human genome have in-built unique structural and energetic signals

Precise identification of correct exon–intron boundaries is a prerequisite to analyze the location and structure of genes. The existing framework for genomic signals, delineating exon and introns in a genomic segment, seems insufficient, predominantly due to poor sequence consensus as well as limitations of training on available experimental data sets. We present here a novel concept for characterizing exon–intron boundaries in genomic segments on the basis of structural and energetic properties. We analyzed boundary junctions on both sides of all the exons (3 28 368) of protein coding genes from human genome (GENCODE database) using 28 structural and three energy parameters. Study of sequence conservation at these sites shows very poor consensus. It is observed that DNA adopts a unique structural and energy state at the boundary junctions. Also, signals are somewhat different for housekeeping and tissue specific genes. Clustering of 31 parameters into four derived vectors gives some additional insights into the physical mechanisms involved in this biological process. Sites of structural and energy signals correlate well to the positions playing important roles in pre-mRNA splicing.

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.

Detection of aberrant splicing events in RNA-seq data using FRASER

Aberrant splicing is a major cause of rare diseases.  However, its prediction from genome sequence alone remains in most cases inconclusive. Recently, RNA sequencing has proven to be an effective complementary avenue to detect aberrant splicing. Here, we develop FRASER, an algorithm to detect aberrant splicing from RNA sequencing data. Unlike existing methods, FRASER captures not only alternative splicing but also intron retention events. This typically doubles the number of detected aberrant events and identified a pathogenic intron retention in MCOLN1 causing mucolipidosis. FRASER automatically controls for latent confounders, which are widespread and affect sensitivity substantially. Moreover, FRASER is based on a count distribution and multiple testing correction, thus reducing the number of calls by two orders of magnitude over commonly applied z score cutoffs, with a minor loss of sensitivity. Applying FRASER to rare disease diagnostics is demonstrated by reprioritizing a pathogenic aberrant exon truncation in TAZ from a published dataset. FRASER is easy to use and freely available.

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.

NineTeen Complex-subunit Salsa is required for efficient splicing of a subset of introns and dorsal-ventral patterning

The NineTeen Complex (NTC), also known as pre-mRNA-processing factor 19 (Prp19) complex, regulates distinct spliceosome conformational changes necessary for splicing. During Drosophila midblastula transition, splicing is particularly sensitive to mutations in NTC-subunit Fandango, which suggests differential requirements of NTC during development. We show that NTC-subunit Salsa, the Drosophila ortholog of human RNA helicase Aquarius, is rate-limiting for splicing of a subset of small first introns during oogenesis, including the first intron of gurken Germline depletion of Salsa and splice site mutations within gurken first intron impair both adult female fertility and oocyte dorsal-ventral patterning, due to an abnormal expression of Gurken. Supporting causality, the fertility and dorsal-ventral patterning defects observed after Salsa depletion could be suppressed by the expression of a gurken construct without its first intron. Altogether, our results suggest that one of the key rate-limiting functions of Salsa during oogenesis is to ensure the correct expression and efficient splicing of the first intron of gurken mRNA. Retention of gurken first intron compromises the function of this gene most likely because it undermines the correct structure and function of the transcript 5'UTR.

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.

Ribo-Pop: simple, cost-effective, and widely applicable ribosomal RNA depletion

The measurement of RNA abundance derived from massively parallel sequencing experiments is an essential technique. Methods that reduce ribosomal RNA levels are usually required prior to sequencing library construction because ribosomal RNA typically comprises the vast majority of a total RNA sample. For some experiments, ribosomal RNA depletion is favored over poly(A) selection because it offers a more inclusive representation of the transcriptome. However, methods to deplete ribosomal RNA are generally proprietary, complex, inefficient, applicable to only specific species, or compatible with only a narrow range of RNA input levels. Here, we describe Ribo-Pop (ribosomal RNA depletion for popular use), a simple workflow and antisense oligo design strategy that we demonstrate works over a wide input range and can be easily adapted to any organism with a sequenced genome. We provide a computational pipeline for probe selection, a streamlined 20-min protocol, and ready-to-use oligo sequences for several organisms. We anticipate that our simple and generalizable "open source" design strategy would enable virtually any laboratory to pursue full transcriptome sequencing in their organism of interest with minimal time and resource investment.

Intron retention and its impact on gene expression and protein diversity: A review and a practical guide

Intron retention (IR) occurs when a complete and unspliced intron remains in mature mRNA. An increasing body of literature has demonstrated a major role for IR in numerous biological functions, including several that impact human health and disease. Although experimental technologies used to study other forms of mRNA splicing can also be used to investigate IR, a specialized downstream computational analysis is optimal for IR discovery and analysis. Here we provide a review of IR and its biological implications, as well as a practical guide for how to detect and analyze it. Several methods, including long read third generation direct RNA sequencing, are described. We have developed an R package, FakIR, to facilitate the execution of the bioinformatic tasks recommended in this review and a tutorial on how to fit them to users aims. Additionally, we provide guidelines and experimental protocols to validate IR discovery and to evaluate the potential impact of IR on gene expression and protein output. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Splicing Regulation/Alternative Splicing RNA Methods > RNA Analyses in vitro and In Silico.

Overlapping Activities of ELAV/Hu Family RNA Binding Proteins Specify the Extended Neuronal 3' UTR Landscape in Drosophila

The tissue-specific deployment of highly extended neural 3' UTR isoforms, generated by alternative polyadenylation (APA), is a broad and conserved feature of metazoan genomes. However, the factors and mechanisms that control neural APA isoforms are not well understood. Here, we show that three ELAV/Hu RNA binding proteins (Elav, Rbp9, and Fne) have similar capacities to induce a lengthened 3' UTR landscape in an ectopic setting. These factors promote accumulation of chromatin-associated, 3' UTR-extended, nascent transcripts, through inhibition of proximal polyadenylation site (PAS) usage. Notably, Elav represses an unannotated splice isoform of fne, switching the normally cytoplasmic Fne toward the nucleus in elav mutants. We use genomic profiling to reveal strong and broad loss of neural APA in elav/fne double mutant CNS, the first genetic background to largely abrogate this distinct APA signature. Overall, we demonstrate how regulatory interplay and functionally overlapping activities of neural ELAV/Hu RBPs drives the neural APA landscape.

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.

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.

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.

Challenges in detecting and quantifying intron retention from next generation sequencing data

Intron retention (IR) occurs when an intron is transcribed into pre-mRNA and remains in the final mRNA. An increasing body of literature has demonstrated a major role for IR in numerous biological functions and in disease. Here we give an overview of the different computational approaches for detecting IR events from sequencing data. We show that these are based on different biological and computational assumptions that may lead to dramatically different results. We describe the various approaches for mitigating errors in detecting intron retention and for discovering IR signatures between different conditions.

Global Co-transcriptional Splicing in Arabidopsis and the Correlation with Splicing Regulation in Mature RNAs

RNA splicing and spliceosome assembly in eukaryotes occur mainly during transcription. However, co-transcriptional splicing has not yet been explored in plants. Here, we built transcriptomes of nascent chromatin RNAs in Arabidopsis thaliana and showed that nearly all introns undergo co-transcriptional splicing, which occurs with higher efficiency for introns in protein-coding genes than for those in noncoding RNAs. Total intron number and intron position are two predominant features that correlate with co-transcriptional splicing efficiency, and introns with alternative 5' or 3' splice sites are less efficiently spliced. Furthermore, we found that mutations in genes encoding trans-acting proteins lead to more introns with increased splicing defects in nascent RNAs than in mature RNAs, and that introns with increased splicing defects in mature RNAs are inefficiently spliced at the co-transcriptional level. Collectively, our results not only uncovered widespread co-transcriptional splicing in Arabidopsis but also identified features that may affect or be affected by co-transcriptional splicing efficiency.