Adventures in time and space: splicing efficiency and RNA polymerase II elongation rate

Alternative splicing is coupled to gene expression in a subset of variably expressed genes

Numerous factors regulate alternative splicing of human genes at a co-transcriptional level. However, how alternative splicing depends on the regulation of gene expression is poorly understood. We leveraged data from the Genotype-Tissue Expression (GTEx) project to show a significant association of gene expression and splicing for 6874 (4.9%) of 141,043 exons in 1106 (13.3%) of 8314 genes with substantially variable expression in nine GTEx tissues. About half of these exons demonstrate higher inclusion with higher gene expression, and half demonstrate higher exclusion, with the observed direction of coupling being highly consistent across different tissues and in external datasets. The exons differ with respect to multiple characteristics and are enriched for hundreds of isoform-specific Gene Ontology annotations suggesting an important regulatory mechanism. Notably, splicing-expression coupling of exons with roles in JUN and MAP kinase signalling could play an important role during cell division.

Splicing-specific transcriptome-wide association uncovers genetic mechanisms for schizophrenia

Recent studies have highlighted the essential role of RNA splicing, a key mechanism of alternative RNA processing, in establishing connections between genetic variations and disease. Genetic loci influencing RNA splicing variations show considerable influence on complex traits, possibly surpassing those affecting total gene expression. Dysregulated RNA splicing has emerged as a major potential contributor to neurological and psychiatric disorders, likely due to the exceptionally high prevalence of alternatively spliced genes in the human brain. Nevertheless, establishing direct associations between genetically altered splicing and complex traits has remained an enduring challenge. We introduce Spliced-Transcriptome-Wide Associations (SpliTWAS) to integrate alternative splicing information with genome-wide association studies to pinpoint genes linked to traits through exon splicing events. We applied SpliTWAS to two schizophrenia (SCZ) RNA-sequencing datasets, BrainGVEX and CommonMind, revealing 137 and 88 trait-associated exons (in 84 and 67 genes), respectively. Enriched biological functions in the associated gene sets converged on neuronal function and development, immune cell activation, and cellular transport, which are highly relevant to SCZ. SpliTWAS variants impacted RNA-binding protein binding sites, revealing potential disruption of RNA-protein interactions affecting splicing. We extended the probabilistic fine-mapping method FOCUS to the exon level, identifying 36 genes and 48 exons as putatively causal for SCZ. We highlight VPS45 and APOPT1, where splicing of specific exons was associated with disease risk, eluding detection by conventional gene expression analysis. Collectively, this study supports the substantial role of alternative splicing in shaping the genetic basis of SCZ, providing a valuable approach for future investigations in this area.

On a kneading theory for gene-splicing

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

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.

Alternative splicing is coupled to gene expression in a subset of variably expressed genes

Numerous factors regulate alternative splicing of human genes at a co-transcriptional level. However, how alternative splicing depends on the regulation of gene expression is poorly understood. We leveraged data from the Genotype-Tissue Expression (GTEx) project to show a significant association of gene expression and splicing for 6874 (4.9%) of 141,043 exons in 1106 (13.3%) of 8314 genes with substantially variable expression in ten GTEx tissues. About half of these exons demonstrate higher inclusion with higher gene expression, and half demonstrate higher exclusion, with the observed direction of coupling being highly consistent across different tissues and in external datasets. The exons differ with respect to sequence characteristics, enriched sequence motifs, RNA polymerase II binding, and inferred transcription rate of downstream introns. The exons were enriched for hundreds of isoform-specific Gene Ontology annotations, suggesting that the coupling of expression and alternative splicing described here may provide an important gene regulatory mechanism that might be used in a variety of biological contexts. In particular, higher inclusion exons could play an important role during cell division.

PARP1 Regulates Circular RNA Biogenesis though Control of Transcriptional Dynamics

Circular RNAs (circRNAs) are a recently discovered class of RNAs derived from protein-coding genes that have important biological and pathological roles. They are formed through backsplicing during co-transcriptional alternative splicing; however, the unified mechanism that accounts for backsplicing decisions remains unclear. Factors that regulate the transcriptional timing and spatial organization of pre-mRNA, including RNAPII kinetics, the availability of splicing factors, and features of gene architecture, have been shown to influence backsplicing decisions. Poly (ADP-ribose) polymerase I (PARP1) regulates alternative splicing through both its presence on chromatin as well as its PARylation activity. However, no studies have investigated PARP1's possible role in regulating circRNA biogenesis. Here, we hypothesized that PARP1's role in splicing extends to circRNA biogenesis. Our results identify many unique circRNAs in PARP1 depletion and PARylation-inhibited conditions compared to the wild type. We found that while all genes producing circRNAs share gene architecture features common to circRNA host genes, genes producing circRNAs in PARP1 knockdown conditions had longer upstream introns than downstream introns, whereas flanking introns in wild type host genes were symmetrical. Interestingly, we found that the behavior of PARP1 in regulating RNAPII pausing is distinct between these two classes of host genes. We conclude that the PARP1 pausing of RNAPII works within the context of gene architecture to regulate transcriptional kinetics, and therefore circRNA biogenesis. Furthermore, this regulation of PARP1 within host genes acts to fine tune their transcriptional output with implications in gene function.

Full-Length Transcriptome of the Great Himalayan Leaf-Nosed Bats (Hipposideros armiger) Optimized Genome Annotation and Revealed the Expression of Novel Genes

The Great Himalayan Leaf-nosed bat (Hipposideros armiger) is one of the most representative species of all echolocating bats and is an ideal model for studying the echolocation system of bats. An incomplete reference genome and limited availability of full-length cDNAs have hindered the identification of alternatively spliced transcripts, which slowed down related basic studies on bats' echolocation and evolution. In this study, we analyzed five organs from H. armiger for the first time using PacBio single-molecule real-time sequencing (SMRT). There were 120 GB of subreads generated, including 1,472,058 full-length non-chimeric (FLNC) sequences. A total of 34,611 alternative splicing (AS) events and 66,010 Alternative Polyadenylation (APA) sites were detected by transcriptome structural analysis. Moreover, a total of 110,611 isoforms were identified, consisting of 52% new isoforms of known genes and 5% of novel gene loci, as well as 2112 novel genes that have not been annotated before in the current reference genome of H. armiger. Furthermore, several key novel genes, including Pol, RAS, NFKB1, and CAMK4, were identified as being associated with nervous, signal transduction, and immune system processes, which may be involved in regulating the auditory nervous perception and immune system that helps bats to regulate in echolocation. In conclusion, the full-length transcriptome results optimized and replenished existing H. armiger genome annotation in multiple ways and offer advantages for newly discovered or previously unrecognized protein-coding genes and isoforms, which can be used as a reference resource.

Rate of transcription elongation and sequence-specific pausing by RNA polymerase I directly influence rRNA processing

One of the first steps in ribosome biogenesis is transcription of the ribosomal DNA by RNA polymerase I (Pol I). Processing of the resultant rRNA begins cotranscriptionally, and perturbation of Pol I transcription elongation results in defective rRNA processing. Mechanistic insight regarding the link between transcription elongation and ribosome assembly is lacking because of limited in vivo methods to assay Pol I transcription. Here, we use native elongating transcript sequencing (NET-Seq) with a strain of Saccharomyces cerevisiae containing a point mutation in Pol I, rpa190-F1205H, which results in impaired rRNA processing and ribosome assembly. We previously demonstrated that this mutation caused a mild reduction in the transcription elongation rate of Pol I in vitro; however, transcription elongation by the mutant has not been characterized in vivo. Here, our findings demonstrate that the mutant Pol I has an increased pause propensity during processive transcription elongation both in vitro and in vivo. NET-Seq reveals that rpa190-F1205H Pol I displays alternative pause site preferences in vivo. Specifically, the mutant is sensitized to A/G residues in the RNA:DNA hybrid and at the last incorporated nucleotide position. Furthermore, both NET-Seq and EM analysis of Miller chromatin spreads reveal pileups of rpa190-F1205H Pol I throughout the ribosomal DNA, particularly at the 5' end of the 35S gene. This combination of in vitro and in vivo analyses of a Pol I mutant provides novel insights into Pol I elongation properties and indicates how these properties are crucial for efficient cotranscriptional rRNA processing and ribosome assembly.

The Use of circRNAs as Biomarkers of Cancer

CircRNAs are a subclass of lncRNAs that have been found to be abundantly present in a wide range of species, including humans. CircRNAs are generally produced by a noncanonical splicing event called backsplicing that is dependent on the canonical splicing machinery, giving rise to circRNAs classified into three main categories: exonic circRNA, circular intronic RNA, and exon-intron circular RNA. Notably, circRNAs possess functional importance and display their functions through different mechanisms of action including sponging miRNAs, or even being translated into functional proteins. In addition, circRNAs also have great potential as biomarkers, particularly in cancer, thanks to their high stability, tissue type and developmental stage specificity, and their presence in biological fluids, which make them promising candidates as noninvasive biomarkers. In this chapter, we describe the most commonly used techniques for the study of circRNAs as cancer biomarkers, including high-throughput techniques such as RNA-Seq and microarrays, and other methods to analyze the presence of specific circRNAs in patient samples.

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.

The Splicing Factor hnRNP M Is a Critical Regulator of Innate Immune Gene Expression in Macrophages

While transcriptional control of innate immune gene expression is well characterized, almost nothing is known about how pre-mRNA splicing decisions influence, or are influenced by, macrophage activation. Here, we demonstrate that the splicing factor hnRNP M is a critical repressor of innate immune gene expression and that its function is regulated by pathogen sensing cascades. Loss of hnRNP M led to hyperinduction of a unique regulon of inflammatory and antimicrobial genes following diverse innate immune stimuli. While mutating specific serines on hnRNP M had little effect on its ability to control pre-mRNA splicing or transcript levels of housekeeping genes in resting macrophages, it greatly impacted the protein's ability to dampen induction of specific innate immune transcripts following pathogen sensing. These data reveal a previously unappreciated role for pattern recognition receptor signaling in controlling splicing factor phosphorylation and establish pre-mRNA splicing as a critical regulatory node in defining innate immune outcomes.

Single-molecule FRET method to investigate the dynamics of transcription elongation through the nucleosome by RNA polymerase II

Transcription elongation through the nucleosome is a precisely coordinated activity to ensure timely production of RNA and accurate regulation of co-transcriptional histone modifications. Nucleosomes actively participate in transcription regulation at various levels and impose physical barriers to RNA polymerase II (RNAPII) during transcription elongation. Despite its high significance, the detailed dynamics of how RNAPII translocates along nucleosomal DNA during transcription elongation and how the nucleosome structure dynamically conforms to the changes necessary for RNAPII progression remain poorly understood. Transcription elongation through the nucleosome is a complex process and investigating the changes of the nucleosome structure during this process by ensemble measurements is daunting. This is because it is nearly impossible to synchronize elongation complexes within a nucleosome or a sub-nucleosome to a designated location at a high enough efficiency for desired sample homogeneity. Here we review our recently developed single-molecule FRET experimental system and method that has fulfilled this deficiency. With our method, one can follow the changes in the structure of individual nucleosomes during transcription elongation. We demonstrated that this method enables the detailed measurements of the kinetics of transcription elongation through the nucleosome and its regulation by a transcription factor, which can be easily extended to investigations of the roles of environmental variables and histone post-translational modifications in regulating transcription elongation.

Insights into the biogenesis and potential functions of exonic circular RNA

Circular RNAs (circRNAs) exhibit unique properties due to their covalently closed nature. Models of circRNAs synthesis and function are emerging but much remains undefined about this surprisingly prevalent class of RNA. Here, we identified exonic circRNAs from human and mouse RNA-sequencing datasets, documenting multiple new examples. Addressing function, we found that many circRNAs co-sediment with ribosomes, indicative of their translation potential. By contrast, circRNAs with potential to act as microRNA sponges were scarce, with some support for a collective sponge function by groups of circRNAs. Addressing circRNA biogenesis, we delineated several features commonly associated with circRNA occurrence. CircRNA-producing genes tend to be longer and to contain more exons than average. Back-splice acceptor exons are strongly enriched at ordinal position 2 within genes, and circRNAs typically have a short exon span with two exons being the most prevalent. The flanking introns either side of circRNA loci are exceptionally long. Of note also, single-exon circRNAs derive from unusually long exons while multi-exon circRNAs are mostly generated from exons of regular length. These findings independently validate and extend similar observations made in a number of prior studies. Furthermore, we analysed high-resolution RNA polymerase II occupancy data from two separate human cell lines to reveal distinctive transcription dynamics at circRNA-producing genes. Specifically, RNA polymerase II traverses the introns of these genes at above average speed concomitant with an accentuated slow-down at exons. Collectively, these features indicate how a perturbed balance between transcription and linear splicing creates important preconditions for circRNA production. We speculate that these preconditions need to be in place so that looping interactions between flanking introns can promote back-splicing to raise circRNA production to appreciable levels.

Promoter-proximal pausing mediated by the exon junction complex regulates splicing

Promoter-proximal pausing of RNA polymerase II (Pol II) is a widespread transcriptional regulatory step across metazoans. Here we find that the nuclear exon junction complex (pre-EJC) is a critical and conserved regulator of this process. Depletion of pre-EJC subunits leads to a global decrease in Pol II pausing and to premature entry into elongation. This effect occurs, at least in part, via non-canonical recruitment of pre-EJC components at promoters. Failure to recruit the pre-EJC at promoters results in increased binding of the positive transcription elongation complex (P-TEFb) and in enhanced Pol II release. Notably, restoring pausing is sufficient to rescue exon skipping and the photoreceptor differentiation defect associated with depletion of pre-EJC components in vivo. We propose that the pre-EJC serves as an early transcriptional checkpoint to prevent premature entry into elongation, ensuring proper recruitment of RNA processing components that are necessary for exon definition.

Splicing analysis of 26 F8 nucleotide variations using a minigene assay

BACKGROUND

Classically, the study of splicing impact of variation located near the splice site is performed by both in silico and mRNA analysis. However, RNA sample was rarely available.

OBJECTIVE

To characterize a panel of putative haemophilia A splicing variations.

MATERIALS AND METHODS

Twenty-six F8 variations identified from a cohort of 2075 haemophilia A families were studied using both bioinformatic tools and in vitro minigene assays in HeLa and Huh7 cells.

RESULTS

An aberrant splicing was demonstrated for 21/26 tested sequence variations. A good correlation between in silico and in vitro analysis was obtained for variations affecting donor splice site (12/14) and for the synonymous variations located inside an exon (6/6). Conversely, no concordant results were observed for the six variations affecting acceptor splice sites. The variations resulted more frequently in exon skipping (n = 13) than in activation of nearby cryptic splice sites (n = 5), in use of a de novo splice site (n = 2) or in insertion of large intronic sequences (n = 1). This study allowed to reclassify 5 synonymous substitutions c.1167A>G (p.Gln389Gln), c.1569G>T (p.Leu523Leu), c.1752G>A (p.Gln584Gln), c.5586G>A (p.Leu1862Leu) and c.6066C>T (p.Gly2022Gly) as splicing variations. The pathological significance of five variations remained unclear (c.222G>A [p.Thr74Thr], c.237C>T [p.Asn79Asn], c.240C>T [p.Ile80Ile], c.2113+5_2113+8del and c.2113+5G>A).

DISCUSSION

The minigene assay herein gave additional evidences for the clinical significance of 21/26 F8 putative splice site mutations. Such investigation should be performed for each F8 putative splice site variation for which no mRNA sample is available, notably to greatly improve the genetic counselling given to female carriers.

Adaptive transcription-splicing resynchronization upon losing an essential splicing factor

Essential genes form the core of a genome and are therefore thought to be indispensable for cellular viability. However, recent findings have challenged this notion in that cells may survive in the absence of some essential genes provided that relevant genetic modifiers are in existence. We therefore hypothesized that the loss of an essential gene may not always be fatefully detrimental; instead, it may pave the way towards genome evolution. We experimentally tested this hypothesis in the context of pre-messenger RNA splicing by evolving yeast cells harbouring a permanent loss of the essential splicing factor Prp28 in the presence of a genetic modifier. Here, we show that cellular fitness can be restored by compensatory mutations that alter either the splicing machinery per se or the Spt-Ada-Gcn5 acetyltransferase transcription co-activator complex in the cells with no Prp28. Biochemical and genetic analysis revealed that slowing down transcription compensates for splicing deficiency, which in turn boosts cellular fitness. In addition, we found that inefficient splicing also conversely decreases nascent RNA production. Taken together, our data suggest that transcription-splicing synchronization contributes to robustness in the gene-expression pathway and argue that the intrinsic interconnectivity within a biological system can be exploited for compensatory evolution and system re-optimization.

A Defective mRNA Cleavage and Polyadenylation Complex Facilitates Expansions of Transcribed (GAA)n Repeats Associated with Friedreich's Ataxia

Expansions of microsatellite repeats are responsible for numerous hereditary diseases in humans, including myotonic dystrophy and Friedreich's ataxia. Whereas the length of an expandable repeat is the main factor determining disease inheritance, recent data point to genomic trans modifiers that can impact the likelihood of expansions and disease progression. Detection of these modifiers may lead to understanding and treating repeat expansion diseases. Here, we describe a method for the rapid, genome-wide identification of trans modifiers for repeat expansion in a yeast experimental system. Using this method, we found that missense mutations in the endoribonuclease subunit (Ysh1) of the mRNA cleavage and polyadenylation complex dramatically increase the rate of (GAA)_n repeat expansions but only when they are actively transcribed. These expansions correlate with slower transcription elongation caused by the ysh1 mutation. These results reveal an interplay between RNA processing and repeat-mediated genome instability, confirming the validity of our approach.

Messenger RNA Life-Cycle in Cancer Cells: Emerging Role of Conventional and Non-Conventional RNA-Binding Proteins?

Functional specialization of cells and tissues in metazoans require specific gene expression patterns. Biological processes, thus, need precise temporal and spatial coordination of gene activity. Regulation of the fate of messenger RNA plays a crucial role in this context. In the present review, the current knowledge related to the role of RNA-binding proteins in the whole mRNA life-cycle is summarized. This field opens up a new angle for understanding the importance of the post-transcriptional control of gene expression in cancer cells. The emerging role of non-classic RNA-binding proteins is highlighted. The goal of this review is to encourage readers to view, through the mRNA life-cycle, novel aspects of the molecular basis of cancer and the potential to develop RNA-based therapies.

Alternative mRNA processing sites decrease genetic variability while increasing functional diversity

Recent large-scale RNA sequencing efforts have revealed the extensive diversity of mRNA molecules produced from most eukaryotic coding genes, which arises from the usage of alternative, cryptic or non-canonical splicing and intronic polyadenylation sites. The prevailing view regarding the tremendous diversity of coding gene transcripts is that mRNA processing is a flexible and more-or-less noisy process leading to a diversity of proteins on which natural selection can act depending on protein-mediated cellular functions. However, this concept raises two main questions. First, do alternative mRNA processing pathways have a role other than generating mRNA and protein diversity? Second, is the cellular function of mRNA variants restricted to the biogenesis of functional protein isoforms? Here, I propose that the co-transcriptional use of alternative mRNA processing sites allows first, the resolution of co-transcriptional biophysical constraints that may otherwise result in DNA instability, and second, increases the diversity of cellular functions of mRNAs in a manner that is not restricted to protein synthesis.

Quantifying the influence of 5'-RNA modifications on RNA polymerase I activity

For ensemble and single-molecule analyses of transcription, the use of synthetic transcription elongation complexes has been a versatile and powerful tool. However, structural analyses demonstrate that short RNA substrates, often employed in these assays, would occupy space within the RNA polymerase. Most commercial RNA oligonucleotides do not carry a 5'-triphosphate as would be present on a natural, de novo synthesized RNA. To examine the effects of 5'-moities on transcription kinetics, we measured nucleotide addition and 3'-dinucleotide cleavage by eukaryotic RNA polymerase I using 5'-hydroxyl and 5'-triphosphate RNA substrates. We found that 5' modifications had no discernable effect on the kinetics of nucleotide addition; however, we observed clear, but modest, effects on the rate of backtracking and/or dinucleotide cleavage. These data suggest that the 5'-end may influence RNA polymerase translocation, consistent with previous prokaryotic studies, and these findings may have implications on kinetic barriers that confront RNA polymerases during the transition from initiation to elongation.

Splicing and transcription touch base: co-transcriptional spliceosome assembly and function

Several macromolecular machines collaborate to produce eukaryotic messenger RNA. RNA polymerase II (Pol II) translocates along genes that are up to millions of base pairs in length and generates a flexible RNA copy of the DNA template. This nascent RNA harbours introns that are removed by the spliceosome, which is a megadalton ribonucleoprotein complex that positions the distant ends of the intron into its catalytic centre. Emerging evidence that the catalytic spliceosome is physically close to Pol II in vivo implies that transcription and splicing occur on similar timescales and that the transcription and splicing machineries may be spatially constrained. In this Review, we discuss aspects of spliceosome assembly, transcription elongation and other co-transcriptional events that allow the temporal coordination of co-transcriptional splicing.

RNA Binding by Histone Methyltransferases Set1 and Set2

Histone methylation at H3K4 and H3K36 is commonly associated with genes actively transcribed by RNA polymerase II (RNAPII) and is catalyzed by Saccharomyces cerevisiae Set1 and Set2, respectively. Here we report that both methyltransferases can be UV cross-linked to RNA in vivo High-throughput sequencing of the bound RNAs revealed strong Set1 enrichment near the transcription start site, whereas Set2 was distributed along pre-mRNAs. A subset of transcripts showed notably high enrichment for Set1 or Set2 binding relative to RNAPII, suggesting functional posttranscriptional interactions. In particular, Set1 was strongly bound to the SET1 mRNA, Ty1 retrotransposons, and noncoding RNAs from the ribosomal DNA (rDNA) intergenic spacers, consistent with its previously reported silencing roles. Set1 lacking RNA recognition motif 2 (RRM2) showed reduced in vivo cross-linking to RNA and reduced chromatin occupancy. In addition, levels of H3K4 trimethylation were decreased, whereas levels of dimethylation were increased. We conclude that RNA binding by Set1 contributes to both chromatin association and methyltransferase activity.

Splicing Factor Mutations in Myelodysplasias: Insights from Spliceosome Structures

Somatic mutations of pre-mRNA splicing factors recur among patients with myelodysplastic syndrome (MDS) and related malignancies. Although these MDS-relevant mutations alter the splicing of a subset of transcripts, the mechanisms by which these single amino acid substitutions change gene expression remain controversial. New structures of spliceosome intermediates and associated protein complexes shed light on the molecular interactions mediated by 'hotspots' of the SF3B1 and U2AF1 pre-mRNA splicing factors. The frequently mutated SF3B1 residues contact the pre-mRNA splice site. Based on structural homology with other spliceosome subunits, and recent findings of altered RNA binding by mutant U2AF1 proteins, we suggest that affected U2AF1 residues also contact pre-mRNA. Altered pre-mRNA recognition emerges as a molecular theme among MDS-relevant mutations of pre-mRNA splicing factors.

Antisense Oligonucleotides Promote Exon Inclusion and Correct the Common c.-32-13T>G GAA Splicing Variant in Pompe Disease

The most common variant causing Pompe disease is c.-32-13T>G (IVS1) in the acid α-glucosidase (GAA) gene, which weakens the splice acceptor of GAA exon 2 and induces partial and complete exon 2 skipping. It also allows a low level of leaky wild-type splicing, leading to a childhood/adult phenotype. We hypothesized that cis-acting splicing motifs may exist that could be blocked using antisense oligonucleotides (AONs) to promote exon inclusion. To test this, a screen was performed in patient-derived primary fibroblasts using a tiling array of U7 small nuclear RNA (snRNA)-based AONs. This resulted in the identification of a splicing regulatory element in GAA intron 1. We designed phosphorodiamidate morpholino oligomer-based AONs to this element, and these promoted exon 2 inclusion and enhanced GAA enzyme activity to levels above the disease threshold. These results indicate that the common IVS1 GAA splicing variant in Pompe disease is subject to negative regulation, and inhibition of a splicing regulatory element using AONs is able to restore canonical GAA splicing and endogenous GAA enzyme activity.

Perfect timing: splicing and transcription rates in living cells

An important step toward understanding gene regulation is the elucidation of the time necessary for the completion of individual steps. Measurement of reaction rates can reveal potential nodes for regulation. For example, measurements of in vivo transcription elongation rates reveal regulation by DNA sequence, gene architecture, and chromatin. Pre-mRNA splicing is regulated by transcription elongation rates and vice versa, yet the rates of RNA processing reactions remain largely elusive. Since the 1980s, numerous model systems and approaches have been used to determine the precise timing of splicing in vivo. Because splicing can be co-transcriptional, the position of Pol II when splicing is detected has been used as a proxy for time by some investigators. In addition to these 'distance-based' measurements, 'time-based' measurements have been possible through live cell imaging, metabolic labeling of RNA, and gene induction. Yet splicing rates can be convolved by the time it takes for transcription, spliceosome assembly and spliceosome disassembly. The variety of assays and systems used has, perhaps not surprisingly, led to reports of widely differing splicing rates in vivo. Recently, single molecule RNA-seq has indicated that splicing occurs more quickly than previously deduced. Here we comprehensively review these findings and discuss evidence that splicing and transcription rates are closely coordinated, facilitating the efficiency of gene expression. On the other hand, introduction of splicing delays through as yet unknown mechanisms provide opportunity for regulation. More work is needed to understand how cells optimize the rates of gene expression for a range of biological conditions. WIREs RNA 2017, 8:e1401. doi: 10.1002/wrna.1401 For further resources related to this article, please visit the WIREs website.

A Conserved Nuclear Cyclophilin Is Required for Both RNA Polymerase II Elongation and Co-transcriptional Splicing in Caenorhabditis elegans

The elongation phase of transcription by RNA Polymerase II (Pol II) involves numerous events that are tightly coordinated, including RNA processing, histone modification, and chromatin remodeling. RNA splicing factors are associated with elongating Pol II, and the interdependent coupling of splicing and elongation has been documented in several systems. Here we identify a conserved, multi-domain cyclophilin family member, SIG-7, as an essential factor for both normal transcription elongation and co-transcriptional splicing. In embryos depleted for SIG-7, RNA levels for over a thousand zygotically expressed genes are substantially reduced, Pol II becomes significantly reduced at the 3' end of genes, marks of transcription elongation are reduced, and unspliced mRNAs accumulate. Our findings suggest that SIG-7 plays a central role in both Pol II elongation and co-transcriptional splicing and may provide an important link for their coordination and regulation.

Strand-specific, high-resolution mapping of modified RNA polymerase II

Reversible modification of the RNAPII C‐terminal domain links transcription with RNA processing and surveillance activities. To better understand this, we mapped the location of RNAPII carrying the five types of CTD phosphorylation on the RNA transcript, providing strand‐specific, nucleotide‐resolution information, and we used a machine learning‐based approach to define RNAPII states. This revealed enrichment of Ser5P, and depletion of Tyr1P, Ser2P, Thr4P, and Ser7P in the transcription start site (TSS) proximal ~150 nt of most genes, with depletion of all modifications close to the poly(A) site. The TSS region also showed elevated RNAPII relative to regions further 3′, with high recruitment of RNA surveillance and termination factors, and correlated with the previously mapped 3′ ends of short, unstable ncRNA transcripts. A hidden Markov model identified distinct modification states associated with initiating, early elongating and later elongating RNAPII. The initiation state was enriched near the TSS of protein‐coding genes and persisted throughout exon 1 of intron‐containing genes. Notably, unstable ncRNAs apparently failed to transition into the elongation states seen on protein‐coding genes.

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

Protein-coding genes in eukaryotes are transcribed by RNA polymerase II (Pol II) and introns are removed from pre-mRNA by the spliceosome. Understanding the time lag between Pol II progression and splicing could provide mechanistic insights into the regulation of gene expression. Here, we present two single-molecule nascent RNA sequencing methods that directly determine the progress of splicing catalysis as a function of Pol II position. Endogenous genes were analyzed on a global scale in budding yeast. We show that splicing is 50% complete when Pol II is only 45 nt downstream of introns, with the first spliced products observed as introns emerge from Pol II. Perturbations that slow the rate of spliceosome assembly or speed up the rate of transcription caused splicing delays, showing that regulation of both processes determines in vivo splicing profiles. We propose that matched rates streamline the gene expression pathway, while allowing regulation through kinetic competition.

Bidirectional coupling of splicing and ATM signaling in response to transcription-blocking DNA damage

In response to DNA damage cells activate intricate protein networks to ensure genomic fidelity and tissue homeostasis. DNA damage response signaling pathways coordinate these networks and determine cellular fates, in part, by modulating RNA metabolism. Here we discuss a replication-independent pathway activated by transcription-blocking DNA lesions, which utilizes the ATM signaling kinase to regulate spliceosome function in a reciprocal manner. We present a model according to which, displacement of co-transcriptional spliceosomes from lesion-arrested RNA polymerases, culminates in R-loop formation and non-canonical ATM activation. ATM signals in a feed-forward fashion to further impede spliceosome organization and regulates UV-induced gene expression and alternative splicing genome-wide. This reciprocal coupling between ATM and the spliceosome highlights the importance of ATM signaling in the cellular response to transcription-blocking lesions and supports a key role of the splicing machinery in this process.

An intronic RNA structure modulates expression of the mRNA biogenesis factor Sus1

Sus1 is a conserved protein involved in chromatin remodeling and mRNA biogenesis. Unlike most yeast genes, the SUS1 pre-mRNA of Saccharomyces cerevisiae contains two introns and is alternatively spliced, retaining one or both introns in response to changes in environmental conditions. SUS1 splicing may allow the cell to control Sus1 expression, but the mechanisms that regulate this process remain unknown. Using in silico analyses together with NMR spectroscopy, gel electrophoresis, and UV thermal denaturation experiments, we show that the downstream intron (I2) of SUS1 forms a weakly stable, 37-nucleotide stem-loop structure containing the branch site near its apical loop and the 3' splice site after the stem terminus. A cellular assay revealed that two of four mutants containing altered I2 structures had significantly impaired SUS1 expression. Semiquantitative RT-PCR experiments indicated that all mutants accumulated unspliced SUS1 pre-mRNA and/or induced distorted levels of fully spliced mRNA relative to wild type. Concomitantly, Sus1 cellular functions in histone H2B deubiquitination and mRNA export were affected in I2 hairpin mutants that inhibited splicing. This work demonstrates that I2 structure is relevant for SUS1 expression, and that this effect is likely exerted through modulation of splicing.

LEDGF/p75 interacts with mRNA splicing factors and targets HIV-1 integration to highly spliced genes

The host chromatin-binding factor LEDGF/p75 interacts with HIV-1 integrase and directs integration to active transcription units. To understand how LEDGF/p75 recognizes transcription units, we sequenced 1 million HIV-1 integration sites isolated from cultured HEK293T cells. Analysis of integration sites showed that cancer genes were preferentially targeted, raising concerns about using lentivirus vectors for gene therapy. Additional analysis led to the discovery that introns and alternative splicing contributed significantly to integration site selection. These correlations were independent of transcription levels, size of transcription units, and length of the introns. Multivariate analysis with five parameters previously found to predict integration sites showed that intron density is the strongest predictor of integration density in transcription units. Analysis of previously published HIV-1 integration site data showed that integration density in transcription units in mouse embryonic fibroblasts also correlated strongly with intron number, and this correlation was absent in cells lacking LEDGF. Affinity purification showed that LEDGF/p75 is associated with a number of splicing factors, and RNA sequencing (RNA-seq) analysis of HEK293T cells lacking LEDGF/p75 or the LEDGF/p75 integrase-binding domain (IBD) showed that LEDGF/p75 contributes to splicing patterns in half of the transcription units that have alternative isoforms. Thus, LEDGF/p75 interacts with splicing factors, contributes to exon choice, and directs HIV-1 integration to transcription units that are highly spliced.

Genetic interaction mapping reveals a role for the SWI/SNF nucleosome remodeler in spliceosome activation in fission yeast

Although numerous regulatory connections between pre-mRNA splicing and chromatin have been demonstrated, the precise mechanisms by which chromatin factors influence spliceosome assembly and/or catalysis remain unclear. To probe the genetic network of pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe, we constructed an epistatic mini-array profile (E-MAP) and discovered many new connections between chromatin and splicing. Notably, the nucleosome remodeler SWI/SNF had strong genetic interactions with components of the U2 snRNP SF3 complex. Overexpression of SF3 components in ΔSWI/SNF cells led to inefficient splicing of many fission yeast introns, predominantly those with non-consensus splice sites. Deletion of SWI/SNF decreased recruitment of the splicing ATPase Prp2, suggesting that SWI/SNF promotes co-transcriptional spliceosome assembly prior to first step catalysis. Importantly, defects in SWI/SNF as well as SF3 overexpression each altered nucleosome occupancy along intron-containing genes, illustrating that the chromatin landscape both affects—and is affected by—co-transcriptional splicing.

It has recently become apparent that most introns are removed from pre-mRNA while the transcript is still engaged with RNA polymerase II (RNAPII). To gain insight into possible roles for chromatin in co-transcriptional splicing, we generated a genome-wide genetic interaction map in fission yeast and uncovered numerous connections between splicing and chromatin. The SWI/SNF remodeling complex is typically thought to activate gene expression by relieving barriers to polymerase elongation imposed by nucleosomes. Here we show that this remodeler is important for an early step in splicing in which Prp2, an RNA-dependent ATPase, is recruited to the assembling spliceosome to promote catalytic activation. Interestingly, introns with sub-optimal splice sites are particularly dependent on SWI/SNF, suggesting the impact of nucleosome dynamics on the kinetics of spliceosome assembly and catalysis. By monitoring nucleosome occupancy, we show significant alterations in nucleosome density in particular splicing and chromatin mutants, which generally paralleled the levels of RNAPII. Taken together, our findings challenge the notion that nucleosomes simply act as barriers to elongation; rather, we suggest that polymerase pausing at nucleosomes can activate gene expression by allowing more time for co-transcriptional splicing.

Getting up to speed with transcription elongation by RNA polymerase II

Recent advances in sequencing techniques that measure nascent transcripts and that reveal the positioning of RNA polymerase II (Pol II) have shown that the pausing of Pol II in promoter-proximal regions and its release to initiate a phase of productive elongation are key steps in transcription regulation. Moreover, after the release of Pol II from the promoter-proximal region, elongation rates are highly dynamic throughout the transcription of a gene, and vary on a gene-by-gene basis. Interestingly, Pol II elongation rates affect co-transcriptional processes such as splicing, termination and genome stability. Increasing numbers of factors and regulatory mechanisms have been associated with the steps of transcription elongation by Pol II, revealing that elongation is a highly complex process. Elongation is thus now recognized as a key phase in the regulation of transcription by Pol II.

Environmental stresses modulate abundance and timing of alternatively spliced circadian transcripts in Arabidopsis

Environmental stresses profoundly altered accumulation of nonsense mRNAs including intron-retaining (IR) transcripts in Arabidopsis. Temporal patterns of stress-induced IR mRNAs were dissected using both oscillating and non-oscillating transcripts. Broad-range thermal cycles triggered a sharp increase in the long IR CCA1 isoforms and altered their phasing to different times of day. Both abiotic and biotic stresses such as drought or Pseudomonas syringae infection induced a similar increase. Thermal stress induced a time delay in accumulation of CCA1 I4Rb transcripts, whereas functional mRNA showed steady oscillations. Our data favor a hypothesis that stress-induced instabilities of the central oscillator can be in part compensated through fluctuations in abundance and out-of-phase oscillations of CCA1 IR transcripts. Taken together, our results support a concept that mRNA abundance can be modulated through altering ratios between functional and nonsense/IR transcripts. SR45 protein specifically bound to the retained CCA1 intron in vitro, suggesting that this splicing factor could be involved in regulation of intron retention. Transcriptomes of nonsense-mediated mRNA decay (NMD)-impaired and heat-stressed plants shared a set of retained introns associated with stress- and defense-inducible transcripts. Constitutive activation of certain stress response networks in an NMD mutant could be linked to disequilibrium between functional and nonsense mRNAs.