Nascent RNA structure modulates the transcriptional dynamics of RNA polymerases

Collective Variables and Facilitated Conformational Opening during Translocation of Human Mitochondrial RNA Polymerase (POLRMT) from Atomic Simulations

Collective variable (CV) identification is challenging in complex dynamical systems. To study the translocation of a single-subunit RNA polymerase (RNAP) during human mitochondrial transcription, we employed all-atom molecular dynamics (MD) as a vehicle to illustrate CV refinement in conformational samplings and dimension reduction analyses. RNAP translocation is an essential mechanical step of transcription elongation that dictates gene expression. The translocation generally follows from polymerization product release and proceeds to initial binding or preinsertion of incoming nucleotides. The human mitochondrial DNA-dependent RNAP (or POLRMT) plays a critical role in cellular metabolism and can be a key molecular off-target in the design of nucleotide analogue antiviral and antitumor drugs due to its structural similarities with many viral RNAPs or RNA-dependent RNA polymerases (RdRps). While POLRMT shares particularly high structural similarity with bacteriophage T7 RNAP, previous experimental studies and our current simulations suggest that POLRMT's mechanochemical coupling mechanisms may be distinct. In the current work, we modeled POLRMT elongation complexes and performed equilibrium MD simulations on the pre- and post-translocation models, with extensive samplings around two potential translocation paths (with or without coupling to the fingers subdomain conformational change). We then compared time-lagged independent component analysis (tICA) and the neural network implementation of the variational approach for Markov processes (VAMPnets) as dimensional reduction methods on selected atomic coordinate sets to best represent the sampled features from the MD simulations. Our results indicate that POLRMT translocation is likely coupled with NTP binding to enable fingers subdomain opening at post-translocation which would otherwise be nonstabilized, or the translocations may proceed futilely without the fingers opening for incoming NTP initial binding or incorporation. The time scale of the coupled translocation reaches over hundreds of microseconds, as predicted by the VAMPnets analyses. Such a time scale seems to match a last postcatalytic kinetic step suggested for the POLRMT elongation cycle by previous experimental measurements. Our MD simulation studies combining atomic coordinate refinements and dimension reduction analyses on top of extensive conformational samplings thus suggest a variation of Brownian ratcheting in POLRMT translocation, as if the Brownian motions of translocation are coupled with NTP binding, which captures transient fingers subdomain opening to couple the translocation with a sustained fingers opening.

Precursor RNA structural patterns at SF3B1 mutation sensitive cryptic 3' splice sites

SF3B1 is a core component of the spliceosome involved in branch point recognition and 3' splice site selection. SF3B1 mutation is common in myelodysplastic syndrome and other blood disorders. The most common mutation in SF3B1 is K700E, a lysine to glutamic acid change within the pre-mRNA interacting heat repeat domain. A hallmark of SF3B1 mutation is an increased use of cryptic 3' splice sites; however, the properties distinguishing SF3B1-sensitive splice junctions from other alternatively spliced junctions are unknown. We identify a subset of 192 core splice junctions that are mis-spliced with SF3B1 K700E mutation. We use our core set to test whether SF3B1-sensitive splice sites are different from control cryptic 3' splice sites via RNA structural accessibility. As a comparison, we define a set of SF3B1-resistant splice junctions with cryptic splice site use that does not change with SF3B1 K700E mutation. We find sequence differences between SF3B1-sensitive and SF3B1-resistant junctions, particularly at the cryptic sites. SF3B1-sensitive cryptic 3' splice sites are within an extended polypyrimidine tract and have lower splice site strength scores. We develop experimental RNA structure data for 83 SF3B1-sensitive junctions and 39 SF3B1-resistant junctions. We find that the pattern of structural accessibility at the NAG splicing motif in cryptic and canonical 3' splice sites is similar. In addition, this pattern can be found in both SF3B1-resistant and SF3B1-sensitive junctions. However, SF3B1-sensitive junctions have cryptic splice sites that are less structurally distinct from the canonical splice sites. In addition, SF3B1-sensitive splice junctions are overall more flexible than SF3B1-resistant junctions. Our results suggest that the SF3B1-sensitive splice junctions have unique structure and sequence properties, containing poorly differentiated, weak splice sites that lead to altered 3' splice site recognition in the presence of SF3B1 mutation.

Probabilistic and machine-learning methods for predicting local rates of transcription elongation from nascent RNA sequencing data

Rates of transcription elongation vary within and across eukaryotic gene bodies. Here, we introduce new methods for predicting elongation rates from nascent RNA sequencing data. First, we devise a probabilistic model that predicts nucleotide-specific elongation rates as a generalized linear function of nearby genomic and epigenomic features. We validate this model with simulations and apply it to public PRO-seq (Precision Run-On Sequencing) and epigenomic data for four cell types, finding that reductions in local elongation rate are associated with cytosine nucleotides, DNA methylation, splice sites, RNA stem-loops, CTCF (CCCTC-binding factor) binding sites, and several histone marks, including H3K36me3 and H4K20me1. By contrast, increases in local elongation rate are associated with thymines, A+T-rich and low-complexity sequences, and H3K79me2 marks. We then introduce a convolutional neural network that improves our local rate predictions. Our analysis is the first to permit genome-wide predictions of relative nucleotide-specific elongation rates.

The role of RNA structure in 3' end processing in eukaryotes

Maturation of pre-mRNA into fully functional mRNA involves a series of highly coordinated steps that are essential for eukaryotic gene expression. RNA structure has been found to play regulatory roles in many of these steps, including cleavage, polyadenylation, and termination. Recent advances in structure probing techniques have been instrumental in revealing how nascent transcript conformation contributes to these dynamic, co-transcriptional processes. In this review, we present examples where RNA structure affects accessibility and/or function of key processing enzymes, thereby influencing the efficiency and precision of 3' end processing machinery. We also discuss emerging technologies that could further enhance our understanding of RNA structure mediated regulation of 3' end processing.

Conformational Dynamics of Influenza A Virus Ribonucleoprotein Complexes during RNA Synthesis

Viral ribonucleoproteins (vRNPs) are the cornerstones of viral proliferation, as they form the macromolecular complexes that are responsible for the transcription and replication of most single-stranded RNA viruses. The influenza A virus (IAV) polymerase catalyzes RNA synthesis within the context of vRNPs where genomic viral RNA (vRNA) is packaged by the viral nucleoprotein (NP). We used high-speed atomic force microscopy and electron microscopy to study the conformational dynamics of individual IAV recombinant RNPs (rRNPs) during RNA synthesis. The rRNPs present an annular organization that allows for the real-time tracking of conformational changes in the NP-vRNA template caused by the advancing polymerase. We demonstrate that the rRNPs undergo a well-defined conformational cycle during RNA synthesis, which can be interpreted in light of previous transcription models. We also present initial estimations of the average RNA synthesis rate in the rRNP and its dependence on the nucleotide concentration and stability of the nascent RNA secondary structures. Furthermore, we provide evidence that rRNPs can perform consecutive cycles of RNA synthesis, accounting for their ability to recycle and generate multiple copies of RNA.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

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

Analytic delay distributions for a family of gene transcription models

Models intended to describe the time evolution of a gene network must somehow include transcription, the DNA-templated synthesis of RNA, and translation, the RNA-templated synthesis of proteins. In eukaryotes, the DNA template for transcription can be very long, often consisting of tens of thousands of nucleotides, and lengthy pauses may punctuate this process. Accordingly, transcription can last for many minutes, in some cases hours. There is a long history of introducing delays in gene expression models to take the transcription and translation times into account. Here we study a family of detailed transcription models that includes initiation, elongation, and termination reactions. We establish a framework for computing the distribution of transcription times, and work out these distributions for some typical cases. For elongation, a fixed delay is a good model provided elongation is fast compared to initiation and termination, and there are no sites where long pauses occur. The initiation and termination phases of the model then generate a nontrivial delay distribution, and elongation shifts this distribution by an amount corresponding to the elongation delay. When initiation and termination are relatively fast, the distribution of elongation times can be approximated by a Gaussian. A convolution of this Gaussian with the initiation and termination time distributions gives another analytic approximation to the transcription time distribution. If there are long pauses during elongation, because of the modularity of the family of models considered, the elongation phase can be partitioned into reactions generating a simple delay (elongation through regions where there are no long pauses), and reactions whose distribution of waiting times must be considered explicitly (initiation, termination, and motion through regions where long pauses are likely). In these cases, the distribution of transcription times again involves a nontrivial part and a shift due to fast elongation processes.

Kinetic pathway of HIV-1 TAR cotranscriptional folding

The Trans-Activator Receptor (TAR) RNA, located at the 5'-end untranslated region (5' UTR) of the human immunodeficiency virus type 1 (HIV-1), is pivotal in the virus's life cycle. As the initial functional domain, it folds during the transcription of viral mRNA. Although TAR's role in recruiting the Tat protein for trans-activation is established, the detailed kinetic mechanisms at play during early transcription, especially at points of temporary transcriptional pausing, remain elusive. Moreover, the precise physical processes of transcriptional pause and subsequent escape are not fully elucidated. This study focuses on the folding kinetics of TAR and the biological implications by integrating computer simulations of RNA folding during transcription with nuclear magnetic resonance (NMR) spectroscopy data. The findings reveal insights into the folding mechanism of a non-native intermediate that triggers transcriptional pause, along with different folding pathways leading to transcriptional pause and readthrough. The profiling of the cotranscriptional folding pathway and identification of kinetic structural intermediates reveal a novel mechanism for viral transcriptional regulation, which could pave the way for new antiviral drug designs targeting kinetic cotranscriptional folding pathways in viral RNAs.

Reciprocating RNA Polymerase batters through roadblocks

RNA polymerases must transit through protein roadblocks to produce full-length transcripts. Here we report real-time measurements of Escherichia coli RNA polymerase passing through different barriers. As intuitively expected, assisting forces facilitated, and opposing forces hindered, RNA polymerase passage through lac repressor protein bound to natural binding sites. Force-dependent differences were significant at magnitudes as low as 0.2 pN and were abolished in the presence of the transcript cleavage factor GreA, which rescues backtracked RNA polymerase. In stark contrast, opposing forces promoted passage when the rate of RNA polymerase backtracking was comparable to, or faster than the rate of dissociation of the roadblock, particularly in the presence of GreA. Our experiments and simulations indicate that RNA polymerase may transit after roadblocks dissociate, or undergo cycles of backtracking, recovery, and ramming into roadblocks to pass through. We propose that such reciprocating motion also enables RNA polymerase to break protein-DNA contacts that hold RNA polymerase back during promoter escape and RNA chain elongation. This may facilitate productive transcription in vivo.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

Introduction to Optical Tweezers: Background, System Designs, and Applications

Optical tweezers are a means to manipulate objects with light. With the technique, microscopically small objects can be held and steered, allowing for accurate measurement of the forces applied to these objects. Optical tweezers can typically obtain a nanometer spatial resolution, a picoNewton force resolution, and a millisecond time resolution, which makes the technique well suited for the study of biological processes from the single-cell down to the single-molecule level. In this chapter, we aim to provide an introduction to the use of optical tweezers for single-molecule analyses. We start from the basic principles and methodology involved in optical trapping, force calibration, and force measurements. Next, we describe the components of an optical tweezers setup and their experimental relevance. Finally, we will provide an overview of the broad applications in context of biological research, with the emphasis on the measurement modes, experimental assays, and possible combinations with fluorescence microscopy techniques.

DNA-sequence and epigenomic determinants of local rates of transcription elongation

Across all branches of life, transcription elongation is a crucial, regulated phase in gene expression. Many recent studies in eukaryotes have focused on the regulation of promoter-proximal pausing of RNA Polymerase II (Pol II), but rates of productive elongation also vary substantially throughout the gene body, both within and across genes. Here, we introduce a probabilistic model for systematically evaluating potential determinants of the local elongation rate based on nascent RNA sequencing (NRS) data. Our model is derived from a unified model for both the kinetics of Pol II movement along the DNA template and the generation of NRS read counts at steady state. It allows for a continuously variable elongation rate along the gene body, with the rate at each nucleotide defined by a generalized linear relationship with nearby genomic and epigenomic features. High-dimensional feature vectors are accommodated through a sparse-regression extension. We show with simulations that the model allows accurate detection of associated features and accurate prediction of local elongation rates. In an analysis of public PRO-seq and epigenomic data, we identify several features that are strongly associated with reductions in the local elongation rate, including DNA methylation, splice sites, RNA stem-loops, CTCF binding sites, and several histone marks, including H3K36me3 and H4K20me1. By contrast, low-complexity sequences and H3K79me2 marks are associated with increases in elongation rate. In an analysis of DNA k -mers, we find that cytosine nucleotides are strongly associated with reductions in local elongation rate, particularly when preceded by guanines and followed by adenines or thymines. Increases in elongation rate are associated with thymines and A+T-rich k -mers. These associations are generally shared across cell types, and by considering them our model is effective at predicting features of held-out PRO-seq data. Overall, our analysis is the first to permit genome-wide predictions of relative nucleotide-specific elongation rates based on complex sets of genomic and epigenomic covariates. We have made predictions available for the K562, CD14+, MCF-7, and HeLa-S3 cell types in a UCSC Genome Browser track.

Genome-wide probing of eukaryotic nascent RNA structure elucidates cotranscriptional folding and its antimutagenic effect

The transcriptional intermediates of RNAs fold into secondary structures with multiple regulatory roles, yet the details of such cotranscriptional RNA folding are largely unresolved in eukaryotes. Here, we present eSPET-seq (Structural Probing of Elongating Transcripts in eukaryotes), a method to assess the cotranscriptional RNA folding in Saccharomyces cerevisiae. Our study reveals pervasive structural transitions during cotranscriptional folding and overall structural similarities between nascent and mature RNAs. Furthermore, a combined analysis with genome-wide R-loop and mutation rate approximations provides quantitative evidence for the antimutator effect of nascent RNA folding through competitive inhibition of the R-loops, known to facilitate transcription-associated mutagenesis. Taken together, we present an experimental evaluation of cotranscriptional folding in eukaryotes and demonstrate the antimutator effect of nascent RNA folding. These results suggest genome-wide coupling between the processing and transmission of genetic information through RNA folding.

U1 snRNP increases RNA Pol II elongation rate to enable synthesis of long genes

The expansion of introns within mammalian genomes poses a challenge for the production of full-length messenger RNAs (mRNAs), with increasing evidence that these long AT-rich sequences present obstacles to transcription. Here, we investigate RNA polymerase II (RNAPII) elongation at high resolution in mammalian cells and demonstrate that RNAPII transcribes faster across introns. Moreover, we find that this acceleration requires the association of U1 snRNP (U1) with the elongation complex at 5' splice sites. The role of U1 to stimulate elongation rate through introns reduces the frequency of both premature termination and transcriptional arrest, thereby dramatically increasing RNA production. We further show that changes in RNAPII elongation rate due to AT content and U1 binding explain previous reports of pausing or termination at splice junctions and the edge of CpG islands. We propose that U1-mediated acceleration of elongation has evolved to mitigate the risks that long AT-rich introns pose to transcript completion.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

Uncovering the mechanisms of transcription elongation by eukaryotic RNA polymerases I, II, and III

Eukaryotes express three nuclear RNA polymerases (Pols I, II, and III) that are essential for cell survival. Despite extensive investigation of the three Pols, significant knowledge gaps regarding their biochemical properties remain because each Pol has been evaluated independently under disparate experimental conditions and methodologies. To advance our understanding of the Pols, we employed identical in vitro transcription assays for direct comparison of their elongation rates, elongation complex (EC) stabilities, and fidelities. Pol I is the fastest, most likely to misincorporate, forms the least stable EC, and is most sensitive to alterations in reaction buffers. Pol II is the slowest of the Pols, forms the most stable EC, and negligibly misincorporated an incorrect nucleotide. The enzymatic properties of Pol III were intermediate between Pols I and II in all assays examined. These results reveal unique enzymatic characteristics of the Pols that provide new insights into their evolutionary divergence.

The pausing zone and control of RNA polymerase II elongation by Spt5: Implications for the pause-release model

The pause-release model of transcription proposes that 40-100 bases from the start site RNA Pol II pauses, followed by release into productive elongation. Pause release is facilitated by the PTEFb phosphorylation of the RNA Pol II elongation factor, Spt5. We mapped paused polymerases by eNET-seq and found frequent pausing in zones that extend ∼0.3-3 kb into genes even when PTEFb is inhibited. The fraction of paused polymerases or pausing propensity declines gradually over several kb and not abruptly as predicted for a discrete pause-release event. Spt5 depletion extends pausing zones, suggesting that it promotes the maturation of elongation complexes to a low-pausing state. The expression of mutants after Spt5 depletion showed that phosphomimetic substitutions in the CTR1 domain diminished pausing throughout genes. By contrast, mutants that prevent the phosphorylation of the Spt5 RNA-binding domain strengthened pausing. Thus, distinct Spt5 phospho-isoforms set the balance between pausing and elongation.

RNA helicase-dependent gene looping impacts messenger RNA processing

DDX5 and DDX17 are DEAD-box RNA helicase paralogs which regulate several aspects of gene expression, especially transcription and splicing, through incompletely understood mechanisms. A transcriptome analysis of DDX5/DDX17-depleted human cells confirmed the large impact of these RNA helicases on splicing and revealed a widespread deregulation of 3' end processing. In silico analyses and experiments in cultured cells showed the binding and functional contribution of the genome organizing factor CTCF to chromatin sites at or near a subset of DDX5/DDX17-dependent exons that are characterized by a high GC content and a high density of RNA Polymerase II. We propose the existence of an RNA helicase-dependent relationship between CTCF and the dynamics of transcription across DNA and/or RNA structured regions, that contributes to the processing of internal and terminal exons. Moreover, local DDX5/DDX17-dependent chromatin loops spatially connect RNA helicase-regulated exons with their cognate promoter, and we provide the first direct evidence that de novo gene looping modifies alternative splicing and polyadenylation. Overall our findings uncover the impact of DDX5/DDX17-dependent chromatin folding on pre-messenger RNA processing.

Alkaline pH has an unexpected effect on transcriptional pausing during synthesis of the E. coli pH-responsive riboswitch

Riboswitches are 5′-untranslated regions of mRNA that change their conformation in response to ligand binding, allowing post-transcriptional gene regulation. This ligand-based model of riboswitch function has been expanded with the discovery of a “pH-responsive element” (PRE) riboswitch in Escherichia coli. At neutral pH, the PRE folds into a translationally inactive structure with an occluded ribosome-binding sequence, whereas at alkaline pH, the PRE adopts a translationally active structure. This unique riboswitch does not rely on ligand binding in a traditional sense to modulate its alternative folding outcomes. Rather, pH controls riboswitch folding by two possible modes that are yet to be distinguished; pH either regulates the transcription rate of RNA polymerase (RNAP) or acts on the RNA itself. Previous work suggested that RNAP pausing is prolonged by alkaline pH at two sites, stimulating PRE folding into the active structure. To date, there has been no rigorous exploration into how pH influences RNAP pausing kinetics during PRE synthesis. To provide that understanding and distinguish between pH acting on RNAP versus RNA, we investigated RNAP pausing kinetics at key sites for PRE folding under different pH conditions. We find that pH influences RNAP pausing but not in the manner proposed previously. Rather, alkaline pH either decreases or has no effect on RNAP pause longevity, suggesting that the modulation of RNAP pausing is not the sole mechanism by which pH affects PRE folding. These findings invite the possibility that the RNA itself actively participates in the sensing of pH.

Screening thousands of transcribed coding and non-coding regions reveals sequence determinants of RNA polymerase II elongation potential

Precise regulation of transcription by RNA polymerase II (RNAPII) is critical for organismal growth and development. However, what determines whether an engaged RNAPII will synthesize a full-length transcript or terminate prematurely is poorly understood. Notably, RNAPII is far more susceptible to termination when transcribing non-coding RNAs than when synthesizing protein-coding mRNAs, but the mechanisms underlying this are unclear. To investigate the impact of transcribed sequence on elongation potential, we developed a method to screen the effects of thousands of INtegrated Sequences on Expression of RNA and Translation using high-throughput sequencing (INSERT-seq). We found that higher AT content in non-coding RNAs, rather than specific sequence motifs, drives RNAPII termination. Further, we demonstrate that 5' splice sites autonomously stimulate processive transcription, even in the absence of polyadenylation signals. Our results reveal a potent role for the transcribed sequence in dictating gene output and demonstrate the power of INSERT-seq toward illuminating these contributions.

Landscape Zooming toward the Prediction of RNA Cotranscriptional Folding

RNA molecules fold as they are transcribed. Cotranscriptional folding of RNA plays a critical role in RNA functions in vivo. Present computational strategies focus on simulations where large structural changes may not be completely sampled. Here, we describe an alternative approach to predicting cotranscriptional RNA folding by zooming in and out of the RNA folding energy landscape. By classifying the RNA structural ensemble into "partitions" based on long, stable helices, we zoom out of the landscape and predict the overall slow folding kinetics from the interpartition kinetic network, and for each interpartition transition, we zoom in on the landscape to simulate the kinetics. Applications of the model to the 117-nucleotide E. coli SRP RNA and the 59-nucleotide HIV-1 TAR RNA show agreements with the experimental data and new structural and kinetic insights into biologically significant conformational switches and pathways for these important systems. This approach, by zooming in/out of an RNA folding landscape at different resolutions, might allow us to treat large RNAs in vivo with transcriptional pause, transcription speed, and other in vivo effects.

Transcriptional Pausing as a Mediator of Bacterial Gene Regulation

Cellular life depends on transcription of DNA by RNA polymerase to express genetic information. RNA polymerase has evolved not just to read information from DNA and write it to RNA but also to sense and process information from the cellular and extracellular environments. Much of this information processing occurs during transcript elongation, when transcriptional pausing enables regulatory decisions. Transcriptional pauses halt RNA polymerase in response to DNA and RNA sequences and structures at locations and times that help coordinate interactions with small molecules and transcription factors important for regulation. Four classes of transcriptional pause signals are now evident after decades of study: elemental pauses, backtrack pauses, hairpin-stabilized pauses, and regulator-stabilized pauses. In this review, I describe current understanding of the molecular mechanisms of these four classes of pause signals, remaining questions about how RNA polymerase responds to pause signals, and the many exciting directions now open to understand pausing and the regulation of transcript elongation on a genome-wide scale. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Targeting RNA structures in diseases with small molecules

RNA is crucial for gene expression and regulation. Recent advances in understanding of RNA biochemistry, structure and molecular biology have revealed the importance of RNA structure in cellular processes and diseases. Various approaches to discovering drug-like small molecules that target RNA structure have been developed. This review provides a brief introduction to RNA structural biology and how RNA structures function as disease regulators. We summarize approaches to targeting RNA with small molecules and highlight their advantages, shortcomings and therapeutic potential.

A Codon-Pair Bias Associated With Network Interactions in Influenza A, B, and C Genomes

A new codon-pair bias present in the genomes of different types of influenza virus is described. Codons with fewer network interactions are more frequency paired together than other codon-pairs in influenza A, B, and C genomes. A shared feature among three different influenza types suggests an evolutionary bias. Codon-pair preference can affect both speed of protein translation and RNA structure. This newly identified bias may provide insight into drivers of virus evolution.

Dissecting nucleotide selectivity in viral RNA polymerases

Designing antiviral therapeutics is of great concern per current pandemics caused by novel coronavirus or SARS-CoV-2. The core polymerase enzyme in the viral replication/transcription machinery is generally conserved and serves well for drug target. In this work we briefly review structural biology and computational clues on representative single-subunit viral polymerases that are more or less connected with SARS-CoV-2 RNA dependent RNA polymerase (RdRp), in particular, to elucidate how nucleotide substrates and potential drug analogs are selected in the viral genome synthesis. To do that, we first survey two well studied RdRps from Polio virus and hepatitis C virus in regard to structural motifs and key residues that have been identified for the nucleotide selectivity. Then we focus on related structural and biochemical characteristics discovered for the SARS-CoV-2 RdRp. To further compare, we summarize what we have learned computationally from phage T7 RNA polymerase (RNAP) on its stepwise nucleotide selectivity, and extend discussion to a structurally similar human mitochondria RNAP, which deserves special attention as it cannot be adversely affected by antiviral treatments. We also include viral phi29 DNA polymerase for comparison, which has both helicase and proofreading activities on top of nucleotide selectivity for replication fidelity control. The helicase and proofreading functions are achieved by protein components in addition to RdRp in the coronavirus replication-transcription machine, with the proofreading strategy important for the fidelity control in synthesizing a comparatively large viral genome.

Transcription Regulation Through Nascent RNA Folding

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

Residence time analysis of RNA polymerase transcription dynamics: A Bayesian sticky HMM approach

The time spent by a single RNA polymerase (RNAP) at specific locations along the DNA, termed "residence time," reports on the initiation, elongation, and termination stages of transcription. At the single-molecule level, this information can be obtained from dual ultrastable optical trapping experiments, revealing a transcriptional elongation of RNAP interspersed with residence times of variable duration. Successfully discriminating between long and short residence times was used by previous approaches to learn about RNAP's transcription elongation dynamics. Here, we propose an approach based on the Bayesian sticky hidden Markov model that treats all residence times for an Escherichia coli RNAP on an equal footing without a priori discriminating between long and short residence times. Furthermore, our method has two additional advantages: we provide full distributions around key point statistics and directly treat the sequence dependence of RNAP's elongation rate. By applying our approach to experimental data, we find assigned relative probabilities on long versus short residence times, force-dependent average residence time transcription elongation dynamics, ∼10% drop in the average backtracking durations in the presence of GreB, and ∼20% drop in the average residence time as a function of applied force in the presence of RNaseA.

Basic mechanisms and kinetics of pause-interspersed transcript elongation

RNA polymerase pausing during elongation is an important mechanism in the regulation of gene expression. Pausing along DNA templates is thought to be induced by distinct signals encoded in the nucleic acid sequence and halt elongation complexes to allow time for necessary co-transcriptional events. Pausing signals have been classified as those producing short-lived elemental, long-lived backtracked, or hairpin-stabilized pauses. In recent years, structural microbiology and single-molecule studies have significantly advanced our understanding of the paused states, but the dynamics of these states are still uncertain, although several models have been proposed to explain the experimentally observed pausing behaviors. This review summarizes present knowledge about the paused states, discusses key discrepancies among the kinetic models and their basic assumptions, and highlights the importance and challenges in constructing theoretical models that may further our biochemical understanding of transcriptional pausing.

Causes and consequences of RNA polymerase II stalling during transcript elongation

The journey of RNA polymerase II (Pol II) as it transcribes a gene is anything but a smooth ride. Transcript elongation is discontinuous and can be perturbed by intrinsic regulatory barriers, such as promoter-proximal pausing, nucleosomes, RNA secondary structures and the underlying DNA sequence. More substantial blocking of Pol II translocation can be caused by other physiological circumstances and extrinsic obstacles, including other transcribing polymerases, the replication machinery and several types of DNA damage, such as bulky lesions and DNA double-strand breaks. Although numerous different obstacles cause Pol II stalling or arrest, the cell somehow distinguishes between them and invokes different mechanisms to resolve each roadblock. Resolution of Pol II blocking can be as straightforward as temporary backtracking and transcription elongation factor S-II (TFIIS)-dependent RNA cleavage, or as drastic as premature transcription termination or degradation of polyubiquitylated Pol II and its associated nascent RNA. In this Review, we discuss the current knowledge of how these different Pol II stalling contexts are distinguished by the cell, how they overlap with each other, how they are resolved and how, when unresolved, they can cause genome instability.

RNA polymerase II-binding aptamers in human ACRO1 satellites disrupt transcription in cis

Transcription elongation is a highly regulated process affected by many proteins, RNAs and the underlying DNA. Here we show that the nascent RNA can interfere with transcription in human cells, extending our previous findings from bacteria and yeast. We identified a variety of Pol II-binding aptamers (RAPs), prominent in repeat elements such as ACRO1 satellites, LINE1 retrotransposons and CA simple repeats, and also in several protein-coding genes. ACRO1 repeat, when translated in silico, exhibits ~50% identity with the Pol II CTD sequence. Taken together with a recent proposal that proteins in general tend to interact with RNAs similar to their cognate mRNAs, this suggests a mechanism for RAP binding. Using a reporter construct, we show that ACRO1 potently inhibits Pol II elongation in cis. We propose a novel mode of transcriptional regulation in humans, in which the nascent RNA binds Pol II to silence its own expression.

So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes

Transcription of DNA into mRNA and translation of mRNA into proteins are two major processes underlying gene expression. Due to the distinct molecular mechanisms, timings, and locales of action, these processes are mainly considered to be independent. During the last two decades, however, multiple factors and elements were shown to coordinate transcription and translation, suggesting an intricate level of synchronization. This review discusses the molecular mechanisms that impact both processes in eukaryotic cells of different origins. The emerging global picture suggests evolutionarily conserved regulation and coordination between transcription and mRNA translation, indicating the importance of this phenomenon for the fine-tuning of gene expression and the adjustment to constantly changing conditions.

Computational discovery and modeling of novel gene expression rules encoded in the mRNA

The transcript is populated with numerous overlapping codes that regulate all steps of gene expression. Deciphering these codes is very challenging due to the large number of variables involved, the non-modular nature of the codes, biases and limitations in current experimental approaches, our limited knowledge in gene expression regulation across the tree of life, and other factors. In recent years, it has been shown that computational modeling and algorithms can significantly accelerate the discovery of novel gene expression codes. Here, we briefly summarize the latest developments and different approaches in the field.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Nascent Transcript Folding Plays a Major Role in Determining RNA Polymerase Elongation Rates

Transcription elongation rates influence RNA processing, but sequence-specific regulation is poorly understood. We addressed this in vivo, analyzing RNAPI in S. cerevisiae. Mapping RNAPI by Miller chromatin spreads or UV crosslinking revealed 5' enrichment and strikingly uneven local polymerase occupancy along the rDNA, indicating substantial variation in transcription speed. Two features of the nascent transcript correlated with RNAPI distribution: folding energy and GC content in the transcription bubble. In vitro experiments confirmed that strong RNA structures close to the polymerase promote forward translocation and limit backtracking, whereas high GC in the transcription bubble slows elongation. A mathematical model for RNAPI elongation confirmed the importance of nascent RNA folding in transcription. RNAPI from S. pombe was similarly sensitive to transcript folding, as were S. cerevisiae RNAPII and RNAPIII. For RNAPII, unstructured RNA, which favors slowed elongation, was associated with faster cotranscriptional splicing and proximal splice site use, indicating regulatory significance for transcript folding.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

High-resolution modeling of the selection on local mRNA folding strength in coding sequences across the tree of life

BACKGROUND

mRNA can form local secondary structure within the protein-coding sequence, and the strength of this structure is thought to influence gene expression regulation. Previous studies suggest that secondary structure strength may be maintained under selection, but the details of this phenomenon are not well understood.

RESULTS

We perform a comprehensive study of the selection on local mRNA folding strengths considering variation between species across the tree of life. We show for the first time that local folding strength selection tends to follow a conserved characteristic profile in most phyla, with selection for weak folding at the two ends of the coding region and for strong folding elsewhere in the coding sequence, with an additional peak of selection for strong folding located downstream of the start codon. The strength of this pattern varies between species and organism groups, and we highlight contradicting cases. To better understand the underlying evolutionary process, we show that selection strengths in the different regions are strongly correlated, and report four factors which have a clear predictive effect on local mRNA folding selection within the coding sequence in different species.

CONCLUSIONS

The correlations observed between selection for local secondary structure strength in the different regions and with the four genomic and environmental factors suggest that they are shaped by the same evolutionary process throughout the coding sequence, and might be maintained under direct selection related to optimization of gene expression and specifically translation regulation.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

Enhancer RNAs in cancer: regulation, mechanisms and therapeutic potential

Enhancers are distal genomic elements critical for gene regulation and cell identify control during development and diseases. Many human cancers were found to associate with enhancer malfunction, due to genetic and epigenetic alterations, which in some cases directly drive tumour growth. Conventionally, enhancers are known to provide DNA binding motifs to recruit transcription factors (TFs) and to control target genes. However, recent progress found that most, if not all, active enhancers pervasively transcribe noncoding RNAs that are referred to as enhancer RNAs (eRNAs). Increasing evidence points to functional roles of at least a subset of eRNAs in gene regulation in both normal and cancer cells, adding new insights into the action mechanisms of enhancers. eRNA expression was observed to be widespread but also specific to tumour types and individual patients, serving as opportunities to exploit them as potential diagnosis markers or therapeutic targets. In this review, we discuss the brief history of eRNA research, their functional mechanisms and importance in cancer gene regulation, as well as their therapeutic and diagnostic values in cancer. We propose that further studies of eRNAs in cancer will offer a promising 'eRNA targeted therapy' for human cancer intervention.

Oxazinomycin arrests RNA polymerase at the polythymidine sequences

Oxazinomycin is a C-nucleoside antibiotic that is produced by Streptomyces hygroscopicus and closely resembles uridine. Here, we show that the oxazinomycin triphosphate is a good substrate for bacterial and eukaryotic RNA polymerases (RNAPs) and that a single incorporated oxazinomycin is rapidly extended by the next nucleotide. However, the incorporation of several successive oxazinomycins or a single oxazinomycin in a certain sequence context arrested a fraction of the transcribing RNAP. The addition of Gre RNA cleavage factors eliminated the transcriptional arrest at a single oxazinomycin and shortened the nascent RNAs arrested at the polythymidine sequences suggesting that the transcriptional arrest was caused by backtracking of RNAP along the DNA template. We further demonstrate that the ubiquitous C-nucleoside pseudouridine is also a good substrate for RNA polymerases in a triphosphorylated form but does not inhibit transcription of the polythymidine sequences. Our results collectively suggest that oxazinomycin functions as a Trojan horse substrate and its inhibitory effect is attributable to the oxygen atom in the position corresponding to carbon five of the uracil ring.

Single-nucleotide control of tRNA folding cooperativity under near-cellular conditions

RNA folding is often studied by renaturing full-length RNA in vitro and tracking folding transitions. However, the intracellular transcript folds as it emerges from the RNA polymerase. Here, we investigate the folding pathways and stability of numerous late-transcriptional intermediates of yeast and Escherichia coli transfer RNAs (tRNAs). Transfer RNA is a highly regulated functional RNA that undergoes multiple steps of posttranscriptional processing and is found in very different lengths during its lifetime in the cell. The precursor transcript is extended on both the 5' and 3' ends of the cloverleaf core, and these extensions get trimmed before addition of the 3'-CCA and aminoacylation. We studied the thermodynamics and structures of the precursor tRNA and of late-transcriptional intermediates of the cloverleaf structure. We examined RNA folding at both the secondary and tertiary structural levels using multiple biochemical and biophysical approaches. Our findings suggest that perhaps nature has selected for a single-base addition to control folding to the functional 3D structure. In near-cellular conditions, yeast tRNA^Phe and E. coli tRNA^Ala transcripts fold in a single, cooperative transition only when nearly all of the nucleotides in the cloverleaf are transcribed by indirectly enhancing folding cooperativity. Furthermore, native extensions on the 5' and 3' ends do not interfere with cooperative core folding. This highly controlled cooperative folding has implications for recognition of tRNA by processing and modification enzymes and quality control of tRNA in cells.

C, Da LT, Yu J. A Viral T7 RNA Polymerase Ratcheting Along DNA With Fidelity Control

RNA polymerase (RNAP) from bacteriophage T7 is a representative single-subunit viral RNAP that can transcribe with high promoter activities without assistances from transcription factors. We accordingly studied this small transcription machine computationally as a model system to understand underlying mechanisms of mechano-chemical coupling and fidelity control in the RNAP transcription elongation. Here we summarize our computational work from several recent publications to demonstrate first how T7 RNAP translocates via Brownian alike motions along DNA right after the catalytic product release. Then we show how the backward translocation motions are prevented at post-translocation upon successful nucleotide incorporation, which is also subject to stepwise nucleotide selection and acts as a pawl for "selective ratcheting". The structural dynamics and energetics features revealed from our atomistic molecular dynamics (MD) simulations and related analyses on the single-subunit T7 RNAP thus provided detailed and quantitative characterizations on the Brownian-ratchet working scenario of a prototypical transcription machine with sophisticated nucleotide selectivity for fidelity control. The presented mechanisms can be more or less general for structurally similar viral or mitochondrial RNAPs and some of DNA polymerases, or even for the RNAP engine of the more complicated transcription machinery in higher organisms.

The fitness consequences of synonymous mutations in Escherichia coli: Experimental evidence for a pleiotropic effect of translational selection

Codon usage bias (CUB) is a universal feature of genomes, and in most species CUB of protein coding genes is positively correlated with expression level and degree of evolutionary conservation. There is mounting experimental evidence that CUB is due in part to selection for translational efficiency and/or accuracy, i.e., translational selection. However, there is a paucity of experimental data on whether and how CUB acts in trans - does the usage of preferred codons in a highly expressed gene affect the translation of other genes by freeing up more ribosomes, thereby increasing their availability to translate all mRNA transcripts in the cell? We investigated this question by creating two extreme versions of the highly expressed Escherichia coli β-lactamase (bla) gene, one comprised almost entirely of unpreferred codons, and a second comprised almost entirely of preferred codons. We monitored the fitness effects of these synonymous mutations over hundreds of generations in two selective environments that allowed us to disentangle translational effects acting in cis from those acting in trans. In a selective environment for maximizing translational efficiency in trans of a gene (tetA) encoding a tetracycline resistance protein, unpreferred synonymous mutations had a negative impact on long-term fitness, whereas preferred mutations had a positive impact on long-term fitness, providing strong experimental evidence for a pleiotropic effect of translational selection.

Transcription factor regulation of RNA polymerase's torque generation capacity

During transcription, RNA polymerase (RNAP) supercoils DNA as it translocates. The resulting torsional stress in DNA can accumulate and, in the absence of regulatory mechanisms, becomes a barrier to RNAP elongation, causing RNAP stalling, backtracking, and transcriptional arrest. Here we investigate whether and how a transcription factor may regulate both torque-induced Escherichia coli RNAP stalling and the torque generation capacity of RNAP. Using a unique real-time angular optical trapping assay, we found that RNAP working against a resisting torque was highly prone to extensive backtracking. We then investigated transcription in the presence of GreB, a transcription factor known to rescue RNAP from the backtracked state. We found that GreB greatly suppressed RNAP backtracking and remarkably increased the torque that RNAP was able to generate by 65%, from 11.2 pN⋅nm to 18.5 pN·nm. Variance analysis of the real-time positional trajectories of RNAP after a stall revealed the kinetic parameters of backtracking and GreB rescue. These results demonstrate that backtracking is the primary mechanism by which torsional stress limits transcription and that the transcription factor GreB effectively enhances the torsional capacity of RNAP. These findings suggest a broader role for transcription factors in regulating RNAP functionality and elongation.

TEFM Enhances Transcription Elongation by Modifying mtRNAP Pausing Dynamics

Regulation of transcription elongation is one of the key mechanisms employed to control gene expression. The single-subunit mitochondrial RNA polymerase (mtRNAP) transcribes mitochondrial genes, such as those related to ATP synthesis. We investigated how mitochondrial transcription elongation factor (TEFM) enhances mtRNAP transcription elongation using a single-molecule optical-tweezers transcription assay, which follows transcription dynamics in real time and allows the separation of pause-free elongation from transcriptional pauses. We found that TEFM enhances the stall force of mtRNAP. Although TEFM does not change the pause-free elongation rate, it enhances mtRNAP transcription elongation by reducing the frequency of long-lived pauses and shortening their durations. Furthermore, we demonstrate how mtRNAP passes through the conserved sequence block II, which is the key sequence for the switch between DNA replication and transcription in mitochondria. Our findings elucidate how both TEFM and mitochondrial genomic DNA sequences directly control the transcription elongation dynamics of mtRNAP.

Dynamic turnover of paused Pol II complexes at human promoters

Paused RNA polymerase II (Pol II) that piles up near most human promoters is the target of mechanisms that control entry into productive elongation. Whether paused Pol II is a stable or dynamic target remains unresolved. We report that most 5' paused Pol II throughout the genome is turned over within 2 min. This process is revealed under hypertonic conditions that prevent Pol II recruitment to promoters. This turnover requires cell viability but is not prevented by inhibiting transcription elongation, suggesting that it is mediated at the level of termination. When initiation was prevented by triptolide during recovery from high salt, a novel preinitiated state of Pol II lacking the pausing factor Spt5 accumulated at transcription start sites. We propose that Pol II occupancy near 5' ends is governed by a cycle of ongoing assembly of preinitiated complexes that transition to pause sites followed by eviction from the DNA template. This model suggests that mechanisms regulating the transition to productive elongation at pause sites operate on a dynamic population of Pol II that is turning over at rates far higher than previously suspected. We suggest that a plausible alternative to elongation control via escape from a stable pause is by escape from premature termination.

Synonymous nucleotide modification of the KCNH2 gene affects both mRNA characteristics and translation of the encoded hERG ion channel

Synonymous nucleotide variation is increasingly recognized as a factor than can affect protein expression, but the underlying mechanisms are incompletely understood. Here, we investigated whether synonymous changes could affect expression of the potassium voltage-gated channel subfamily H member 2 (KCNH2) gene, encoding the human ether-a-go-go-related gene (hERG) ion channel, which is linked to hereditary cardiac arrhythmia. We examined a previously described synthetic version (hERG-codon modified (CM)) with synonymous substitutions designed to reduce GC content, rare codons, and mRNA secondary structure relative to the native construct (hERG-NT). hERG-CM exhibited lower protein expression than hERG-NT in HEK293T cells. We found that the steady-state abundance of hERG-NT mRNA was greater than hERG-CM because of an enhanced transcription rate and increased mRNA stability for hERG-NT. Translation of hERG-CM was independently reduced, contributing to the overall greater synthesis of hERG-NT channel protein. This was partially offset, however, by a higher aggregation of a newly synthesized hERG-NT channel, resulting in nonfunctional protein. Regional mRNA analyses of chimeras of hERG-NT and hERG-CM revealed that synonymous changes in the 5' segments of the coding region had the greatest influence on hERG synthesis at both the mRNA and protein levels. Taken together, these results indicate that synonymous nucleotide variations within the coding region, particularly in the 5' region of the hERG mRNA, can affect both transcription and translation. These findings support the notion that greater attention should be given to the effects of synonymous genetic variation when analyzing hERG DNA sequences in the study of hereditary cardiac disease.

Pause sequences facilitate entry into long-lived paused states by reducing RNA polymerase transcription rates

Transcription by RNA polymerase (RNAP) is interspersed with sequence-dependent pausing. The processes through which paused states are accessed and stabilized occur at spatiotemporal scales beyond the resolution of previous methods, and are poorly understood. Here, we combine high-resolution optical trapping with improved data analysis methods to investigate the formation of paused states at enhanced temporal resolution. We find that pause sites reduce the forward transcription rate of nearly all RNAP molecules, rather than just affecting the subset of molecules that enter long-lived pauses. We propose that the reduced rates at pause sites allow time for the elongation complex to undergo conformational changes required to enter long-lived pauses. We also find that backtracking occurs stepwise, with states backtracked by at most one base pair forming quickly, and further backtracking occurring slowly. Finally, we find that nascent RNA structures act as modulators that either enhance or attenuate pausing, depending on the sequence context.

Codon-Resolution Analysis Reveals a Direct and Context-Dependent Impact of Individual Synonymous Mutations on mRNA Level

Codon usage bias (CUB) refers to the observation that synonymous codons are not used equally frequently in a genome. CUB is stronger in more highly expressed genes, a phenomenon commonly explained by stronger natural selection on translational accuracy and/or efficiency among these genes. Nevertheless, this phenomenon could also occur if CUB regulates gene expression at the mRNA level, a hypothesis that has not been tested until recently. Here, we attempt to quantify the impact of synonymous mutations on mRNA level in yeast using 3,556 synonymous variants of a heterologous gene encoding green fluorescent protein (GFP) and 523 synonymous variants of an endogenous gene TDH3. We found that mRNA level was positively correlated with CUB among these synonymous variants, demonstrating a direct role of CUB in regulating transcript concentration, likely via regulating mRNA degradation rate, as our additional experiments suggested. More importantly, we quantified the effects of individual synonymous mutations on mRNA level and found them dependent on 1) CUB and 2) mRNA secondary structure, both in proximal sequence contexts. Our study reveals the pleiotropic effects of synonymous codon usage and provides an additional explanation for the well-known correlation between CUB and gene expression level.