Transcription Elongation Rate Research Articles

Beyond Transcription: Roles of Transcription Factors in Pre-mRNA Splicing.

Whereas individual steps of protein-coding gene expression in eukaryotes can be studied in isolation in vitro, it has become clear that these steps are intimately connected within cells. Connections not only ensure quality control but also fine-tune the gene expression process, which must adapt to environmental changes while remaining robust. In this review, we systematically present proven and potential mechanisms by which sequence-specific DNA-binding transcription factors can alter gene expression beyond transcription initiation and regulate pre-mRNA splicing, and thereby mRNA isoform production, by (i) influencing transcription elongation rates, (ii) binding to pre-mRNA to recruit splicing factors, and/or (iii) blocking the association of splicing factors with pre-mRNA. We propose various mechanistic models throughout the review, in some cases without explicit supportive evidence, in hopes of providing fertile ground for future studies.

Functions for fission yeast splicing factors SpSlu7 and SpPrp18 in alternative splice-site choice and stress-specific regulated splicing.

Budding yeast spliceosomal factors ScSlu7 and ScPrp18 interact and mediate intron 3’ss choice during second step pre-mRNA splicing. The fission yeast genome with abundant multi-intronic transcripts, degenerate splice signals and SR proteins is an apt unicellular fungal model to deduce roles for core spliceosomal factors in alternative splice-site choice, intron retention and to study the cellular implications of regulated splicing. From our custom microarray data we deduce a stringent reproducible subset of S. pombe alternative events. We examined the role of factors SpSlu7 or SpPrp18 for these splice events and investigated the relationship to growth phase and stress. Wild-type log and stationary phase cells showed ats1+ exon 3 skipped and intron 3 retained transcripts. Interestingly the non-consensus 5’ss in ats1+ intron 3 caused SpSlu7 and SpPrp18 dependent intron retention. We validated the use of an alternative 5’ss in dtd1+ intron 1 and of an upstream alternative 3’ss in DUF3074 intron 1. The dtd1+ intron 1 non-canonical 5’ss yielded an alternative mRNA whose levels increased in stationary phase. Utilization of dtd1+ intron 1 sub-optimal 5’ ss required functional SpPrp18 and SpSlu7 while compromise in SpSlu7 function alone hampered the selection of the DUF3074 intron 1 non canonical 3’ss. We analysed the relative abundance of these splice isoforms during mild thermal, oxidative and heavy metal stress and found stress-specific splice patterns for ats1+ and DUF3074 intron 1 some of which were SpSlu7 and SpPrp18 dependent. By studying ats1+ splice isoforms during compromised transcription elongation rates in wild-type, spslu7-2 and spprp18-5 mutant cells we found dynamic and intron context-specific effects in splice-site choice. Our work thus shows the combinatorial effects of splice site strength, core splicing factor functions and transcription elongation kinetics to dictate alternative splice patterns which in turn serve as an additional recourse of gene regulation in fission yeast.

A code for transcription elongation speed

ABSTRACTThe two major steps of gene expression are transcription and translation. While hundreds of studies regarding the effect of sequence features on the translation elongation process have been published, very few connect sequence features to the transcription elongation rate. We suggest, for the first time, that short transcript sub-sequences have a typical effect on RNA polymerase (RNAP) speed: we show that nucleotide 5-mers tend to have typical RNAP speed (or transcription rate), which is consistent along different parts of genes and among different groups of genes with high correlation. We also demonstrate that relative RNAP speed correlates with mRNA levels of endogenous and heterologous genes. Furthermore, we show that the estimated transcription and translation elongation rates correlate in endogenous genes. Finally, we demonstrate that our results are consistent for different high resolution experimental measurements of RNAP densities. These results suggest for the first time that transcription elongation is partly encoded in the transcript, affected by the codon-usage, and optimized by evolution with a significant effect on gene expression and organismal fitness.

Determinants of Histone H3K4 Methylation Patterns

Various factors differentially recognize trimethylated histone H3 lysine 4 (H3K4me3) near promoters, H3K4me2 just downstream, and promoter-distal H3K4me1 to modulate gene expression. This methylation "gradient" is thought to result from preferential binding of the H3K4 methyltransferase Set1/complex associated with Set1 (COMPASS) to promoter-proximalRNA polymerase II. However, other studies havesuggested that location-specific cues allosterically activate Set1. Chromatin immunoprecipitationsequencing (ChIP-seq) experiments show that H3K4 methylation patterns on active genes are not universal or fixed and change in response to both transcription elongation rate and frequency as well as reduced COMPASS activity. Fusing Set1 to RNA polymerase II results in H3K4me2 throughout transcribed regions and similarly extended H3K4me3 on highly transcribed genes. Tethered Set1 still requires histone H2B ubiquitylation for activity. These results show that higher-level methylations reflect not only Set1/COMPASS recruitment but also multiple rounds of transcription. This model provides a simple explanation for non-canonical methylationpatterns at some loci or in certain COMPASS mutants.

Transcription elongation rate has a tissue-specific impact on alternative cleavage and polyadenylation in Drosophila melanogaster.

Alternative polyadenylation (APA) is a mechanism that generates multiple mRNA isoforms with different 3′UTRs and/or coding sequences from a single gene. Here, using 3′ region extraction and deep sequencing (3′READS), we have systematically mapped cleavage and polyadenylation sites (PASs) in Drosophila melanogaster, expanding the total repertoire of PASs previously identified for the species, especially those located in A-rich genomic sequences. Cis-element analysis revealed distinct sequence motifs around fly PASs when compared to mammalian ones, including the greater enrichment of upstream UAUA elements and the less prominent presence of downstream UGUG elements. We found that over 75% of mRNA genes in Drosophila melanogaster undergo APA. The head tissue tends to use distal PASs when compared to the body, leading to preferential expression of APA isoforms with long 3′UTRs as well as with distal terminal exons. The distance between the APA sites and intron location of PAS are important parameters for APA difference between body and head, suggesting distinct PAS selection contexts. APA analysis of the RpII215C4 mutant strain, which harbors a mutant RNA polymerase II (RNAPII) with a slower elongation rate, revealed that a 50% decrease in transcriptional elongation rate leads to a mild trend of more usage of proximal, weaker PASs, both in 3′UTRs and in introns, consistent with the “first come, first served” model of APA regulation. However, this trend was not observed in the head, suggesting a different regulatory context in neuronal cells. Together, our data expand the PAS collection for Drosophila melanogaster and reveal a tissue-specific effect of APA regulation by RNAPII elongation rate.

Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment

While intron retention (IR) is considered a widely conserved and distinct mechanism of gene expression control, its regulation is poorly understood. Here we show that DNA methylation directly regulates IR. We also find reduced occupancy of MeCP2 near the splice junctions of retained introns, mirroring the reduced DNA methylation at these sites. Accordingly, MeCP2 depletion in tissues and cells enhances IR. By analysing the MeCP2 interactome using mass spectrometry and RNA co-precipitation, we demonstrate that decreased MeCP2 binding near splice junctions facilitates IR via reduced recruitment of splicing factors, including Tra2b, and increased RNA polymerase II stalling. These results suggest an association between IR and a slower rate of transcription elongation, which reflects inefficient splicing factor recruitment. In summary, our results reinforce the interdependency between alternative splicing involving IR and epigenetic controls of gene expression.

A saga of cancer epigenetics: linking epigenetics to alternative splicing

The discovery of an increasing number of alternative splicing events in the human genome highlighted that ∼94% of genes generate alternatively spliced transcripts that may produce different protein isoforms with diverse functions. It is now well known that several diseases are a direct and indirect consequence of aberrant splicing events in humans. In addition to the conventional mode of alternative splicing regulation by 'cis' RNA-binding sites and 'trans' RNA-binding proteins, recent literature provides enormous evidence for epigenetic regulation of alternative splicing. The epigenetic modifications may regulate alternative splicing by either influencing the transcription elongation rate of RNA polymerase II or by recruiting a specific splicing regulator via different chromatin adaptors. The epigenetic alterations and aberrant alternative splicing are known to be associated with various diseases individually, but this review discusses/highlights the latest literature on the role of epigenetic alterations in the regulation of alternative splicing and thereby cancer progression. This review also points out the need for further studies to understand the interplay between epigenetic modifications and aberrant alternative splicing in cancer progression.

Rapid Sampling of Escherichia coli After Changing Oxygen Conditions Reveals Transcriptional Dynamics.

Escherichia coli is able to shift between anaerobic and aerobic metabolism by adapting its gene expression, e.g., of metabolic genes, to the new environment. The dynamics of gene expression that result from environmental shifts are limited, amongst others, by the time needed for regulation and transcription elongation. In this study, we examined gene expression dynamics after an anaerobic-to-aerobic shift on a short time scale (0.5, 1, 2, 5, and 10 min) by RNA sequencing with emphasis on delay times and transcriptional elongation rates (TER). Transient expression patterns and timing of differential expression, characterized by delay and elongation, were identified as key features of the dataset. Gene ontology enrichment analysis revealed early upregulation of respiratory and iron-related gene sets. We inferred specific TERs of 89 operons with a mean TER of 42.0 nt/s and mean delay time of 22.4 s. TERs correlate with sequence features, such as codon bias, whereas delay times correlate with the involvement of regulators. The presented data illustrate that at very short times after a shift in oxygenation, extensional changes of the transcriptome, such as temporary responses, can be observed. Besides regulation, TERs contribute to the dynamics of gene expression.

Measuring In Vivo Supercoil Dynamics and Transcription Elongation Rates in Bacterial Chromosomes.

DNA gyrase is the only topoisomerase that can catalytically introduce negative supercoils into covalently closed DNA. The enzyme plays a critical role in many phases of DNA biochemistry. There are only a few methods that allow one to measure supercoiling in chromosomal DNA and analyze the role of gyrase in transcription and its interaction with the other three bacterial topoisomerases. Here, we provide molecular tools for measuring supercoil density in the chromosome and for connecting the dots between transcription and DNA topology.

Sub1 contacts the RNA polymerase II stalk to modulate mRNA synthesis.

Biogenesis of messenger RNA is critically influenced by the phosphorylation state of the carboxy-terminal domain (CTD) in the largest RNA polymerase II (RNAPII) subunit. Several kinases and phosphatases are required to maintain proper CTD phosphorylation levels and, additionally, several other proteins modulate them, including Rpb4/7 and Sub1. The Rpb4/7 heterodimer, constituting the RNAPII stalk, promote phosphatase functions and Sub1 globally influences CTD phosphorylation, though its mechanism remains mostly unknown. Here, we show that Sub1 physically interacts with the RNAPII stalk domain, Rpb4/7, likely through its C-terminal region, and associates with Fcp1. While Rpb4 is not required for Sub1 interaction with RNAPII complex, a fully functional heterodimer is required for Sub1 association to promoters. We also demonstrate that a complete CTD is necessary for proper association of Sub1 to chromatin and to the RNAPII. Finally, genetic data show a functional relationship between Sub1 and the RNAPII clamp domain. Altogether, our results indicate that Sub1, Rpb4/7 and Fcp1 interaction modulates CTD phosphorylation. In addition, Sub1 interaction with Rpb4/7 can also modulate transcription start site selection and transcription elongation rate likely by influencing the clamp function.

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.

Super elongation complex contains a TFIIF-related subcomplex

ABSTRACTSuper elongation complex (SEC) belongs to a family of RNA polymerase II (Pol II) elongation factors that has similar properties as TFIIF, a general transcription factor that increases the transcription elongation rate by reducing pausing. Although SEC has TFIIF-like functional properties, it apparently lacks sequence and structural homology. Using HHpred, we find that SEC contains an evolutionarily related TFIIF-like subcomplex. We show that the SEC subunit ELL interacts with the Pol II Rbp2 subunit, as expected for a TFIIF-like factor. These findings suggest a new model for how SEC functions as a Pol II elongation factor and how it suppresses Pol II pausing.

Quantitative regulation of FLC via coordinated transcriptional initiation and elongation

The basis of quantitative regulation of gene expression is still poorly understood. In Arabidopsis thaliana, quantitative variation in expression of FLOWERING LOCUS C (FLC) influences the timing of flowering. In ambient temperatures, FLC expression is quantitatively modulated by a chromatin silencing mechanism involving alternative polyadenylation of antisense transcripts. Investigation of this mechanism unexpectedly showed that RNA polymerase II (Pol II) occupancy changes at FLC did not reflect RNA fold changes. Mathematical modeling of these transcriptional dynamics predicted a tight coordination of transcriptional initiation and elongation. This prediction was validated by detailed measurements of total and chromatin-bound FLC intronic RNA, a methodology appropriate for analyzing elongation rate changes in a range of organisms. Transcription initiation was found to vary ∼ 25-fold with elongation rate varying ∼ 8- to 12-fold. Premature sense transcript termination contributed very little to expression differences. This quantitative variation in transcription was coincident with variation in H3K36me3 and H3K4me2 over the FLC gene body. We propose different chromatin states coordinately influence transcriptional initiation and elongation rates and that this coordination is likely to be a general feature of quantitative gene regulation in a chromatin context.

Effects of Transcription Elongation Rate and Xrn2 Exonuclease Activity on RNA Polymerase II Termination Suggest Widespread Kinetic Competition

The torpedo model of transcription termination asserts that the exonuclease Xrn2 attacks the 5'PO4-end exposed by nascent RNA cleavage and chases down the RNA polymerase. We tested this mechanism using a dominant-negative human Xrn2 mutant and found that it delayed termination genome-wide. Xrn2 nuclease inactivation caused strong termination defects downstream of most poly(A) sites and modest delays at some histone and U snRNA genes, suggesting that the torpedo mechanism is not limited to poly(A) site-dependent termination. A central untested feature of the torpedo model is that there is kinetic competition between the exonuclease and the pol II elongation complex. Using pol II rate mutants, we found that slow transcription robustly shifts termination upstream, and fast elongation extends the zone of termination further downstream. These results suggest that kinetic competition between elongating pol II and the Xrn2 exonuclease is integral to termination of transcription on most human genes.

Expanded GAA repeats impede transcription elongation through the FXN gene and induce transcriptional silencing that is restricted to the FXN locus.

Friedreich's ataxia (FRDA) is a severe neurodegenerative disease caused by homozygous expansion of the guanine-adenine-adenine (GAA) repeats in intron 1 of the FXN gene leading to transcriptional repression of frataxin expression. Post-translational histone modifications that typify heterochromatin are enriched in the vicinity of the repeats, whereas active chromatin marks in this region are underrepresented in FRDA samples. Yet, the immediate effect of the expanded repeats on transcription progression through FXN and their long-range effect on the surrounding genomic context are two critical questions that remain unanswered in the molecular pathogenesis of FRDA. To address these questions, we conducted next-generation RNA sequencing of a large cohort of FRDA and control primary fibroblasts. This comprehensive analysis revealed that the GAA-induced silencing effect does not influence expression of neighboring genes upstream or downstream of FXN. Furthermore, no long-range silencing effects were detected across a large portion of chromosome 9. Additionally, results of chromatin immunoprecipitation studies confirmed that histone modifications associated with repressed transcription are confined to the FXN locus. Finally, deep sequencing of FXN pre-mRNA molecules revealed a pronounced defect in the transcription elongation rate in FRDA cells when compared with controls. These results indicate that approaches aimed to reactivate frataxin expression should simultaneously address deficits in transcription initiation and elongation at the FXN locus.

Cwf16p Associating with the Nineteen Complex Ensures Ordered Exon Joining in Constitutive Pre-mRNA Splicing in Fission Yeast.

Exons are ligated in an ordered manner without the skipping of exons in the constitutive splicing of pre-mRNAs with multiple introns. To identify factors ensuring ordered exon joining in constitutive pre-mRNA splicing, we previously screened for exon skipping mutants in Schizosaccharomyces pombe using a reporter plasmid, and characterized three exon skipping mutants named ods1 (ordered splicing 1), ods2, and ods3, the responsible genes of which encode Prp2/U2AF59, U2AF23, and SF1, respectively. They form an SF1-U2AF59-U2AF23 complex involved in recognition of the branch and 3′ splice sites in pre-mRNA. In the present study, we identified a fourth ods mutant, ods4, which was isolated in an exon-skipping screen. The ods4 + gene encodes Cwf16p, which interacts with the NineTeen Complex (NTC), a complex thought to be involved in the first catalytic step of the splicing reaction. We isolated two multi-copy suppressors for the ods4-1 mutation, Srp2p, an SR protein essential for pre-mRNA splicing, and Tif213p, a translation initiation factor, in S. pombe. The overexpression of Srp2p suppressed the exon-skipping phenotype of all ods mutants, whereas Tif213p suppressed only ods4-1, which has a mutation in the translational start codon of the cwf16 gene. We also showed that the decrease in the transcriptional elongation rate induced by drug treatment suppressed exon skipping in ods4-1. We propose that Cwf16p/NTC participates in the early recognition of the branch and 3′ splice sites and cooperates with the SF1-U2AF59-U2AF23 complex to maintain ordered exon joining.

Bacteriophage Xp10 anti-termination factor p7 induces forward translocation by host RNA polymerase.

Regulation of transcription elongation is based on response of RNA polymerase (RNAP) to various pause signals and is modulated by various accessory factors. Here we report that a 7 kDa protein p7 encoded by bacteriophage Xp10 acts as an elongation processivity factor of RNAP of host bacterium Xanthomonas oryzae, a major rice pathogen. Our data suggest that p7 stabilizes the upstream DNA duplex of the elongation complex thus disfavouring backtracking and promoting forward translocated states of the elongation complex. The p7-induced ‘pushing’ of RNAP and modification of RNAP contacts with the upstream edge of the transcription bubble lead to read-through of various types of pauses and termination signals and generally increase transcription processivity and elongation rate, contributing for transcription of an extremely long late genes operon of Xp10. Forward translocation was observed earlier upon the binding of unrelated bacterial elongation factor NusG, suggesting that this may be a general pathway of regulation of transcription elongation.

Coupling of Transcription with mRNA Processing in time and Space

The cellular transcription and mRNA processing machineries appear to have co-evolved to allow coupling of the reactions they perform in space and time. Coupling in space is achieved by “recruitment” of processing factors to the nascent transcript, often by binding to the pol II CTD. Much less is known about how coupling in time is achieved, that is to say the coordination of rates of elongation and processing. Slow elongation is predicted to lengthen the window of opportunity for an event to occur at an upstream site on the nascent transcript before it must compete with an RNA sequence element further downstream. Consistent with this prediction, in a few cases slow elongation is associated with enhanced inclusion of alternative exons (ref. 1) through processing at upstream splice sites and processing at upstream alternative poly(A) sites (ref. 2). We investigated kinetic coupling in human cell lines that express pol II large subunits with mutations in funnel and trigger loop domains that slow down and speed up elongation rates genome-wide. Slow and fast elongation affected constitutive and alternative splicing, frequently altering exon inclusion and intron retention in ways not predicted by the “window of opportunity” model. Unexpectedly, slow and fast elongation often both increased or both decreased inclusion of a particular exon or retained intron. Together our results suggest that an optimal rate of transcriptional elongation is required for normal co-transcriptional pre-mRNA splicing.

Simultaneous measurement of genome-wide transcription elongation speeds and rates of RNA polymerase II transition into active elongation with 4sUDRB-seq.

4sUDRB-seq separately measures, on a genomic scale, the distinct contributions of transcription elongation speed and rate of RNA polymerase II (Pol II) transition into active elongation (TAE) to the overall mRNA production rate. It uses reversible inhibition of transcription elongation with 5,6-dichloro-1-β-D-ribofuranosylbenzimidazole (DRB), combined with a pulse of 4-thiouridine (4sU), to tag newly transcribed RNA. After DRB removal, cells are collected at several time points, and tagged RNA is biotinylated, captured on streptavidin beads and sequenced. 4sUDRB-seq enables the comparison of elongation speeds between different developmental stages or different cell types, and it allows the impact of specific transcription factors on transcription elongation speed versus TAE to be studied. RNA preparation takes ∼4 d to complete, with deep sequencing requiring an additional ∼4-11 d plus 1-3 d for bioinformatics analysis. The experimental protocol requires basic molecular biology skills, whereas data analysis requires knowledge in bioinformatics, particularly MATLAB and the Linux environment.

RNA folding retrospective: lessons from ribozymes big and small.

RNA folding retrospective: lessons from ribozymes big and small.