Simultaneous Measurement of Transcriptional and Post-transcriptional Parameters by 3′ End RNA-Seq

Summary Cellular RNA levels are determined by transcription and decay rates, which are fundamental in understanding gene expression regulation. Measurement of these two parameters is usually performed independently, complicating analysis as well as introducing methodological biases and batch effects that hamper direct comparison. Here, we present a simple approach of concurrent sequencing of S. cerevisiae poly(A)+ and poly(A)− RNA 3′ ends to simultaneously estimate total RNA levels, transcription, and decay rates from the same RNA sample. The transcription data generated correlate well with reported estimates and also reveal local RNA polymerase stalling and termination sites with high precision. Although the method by design uses brief metabolic labeling of newly synthesized RNA with 4-thiouracil, the results demonstrate that transcription estimates can also be gained from unlabeled RNA samples. These findings underscore the potential of the approach, which should be generally applicable to study a range of biological questions in diverse organisms.

[1]  S. Preibisch,et al.  Global analysis of nascent RNA reveals transcriptional pausing in terminal exons. , 2010, Molecular cell.

[2]  Sarah J. Wheelan,et al.  Transcriptome-Wide Binding Sites for Components of the Saccharomyces cerevisiae Non-Poly(A) Termination Pathway: Nrd1, Nab3, and Sen1 , 2011, PLoS genetics.

[3]  Paul Theodor Pyl,et al.  HTSeq—a Python framework to work with high-throughput sequencing data , 2014, bioRxiv.

[4]  L. Steinmetz,et al.  Bidirectional promoters generate pervasive transcription in yeast , 2009, Nature.

[5]  L. Steinmetz,et al.  Roadblock termination by reb1p restricts cryptic and readthrough transcription. , 2014, Molecular cell.

[6]  B. Séraphin,et al.  Futile cycle of transcription initiation and termination modulates the response to nucleotide shortage in S. cerevisiae. , 2008, Molecular cell.

[7]  Achim Tresch,et al.  Comparative dynamic transcriptome analysis (cDTA) reveals mutual feedback between mRNA synthesis and degradation , 2012, Genome research.

[8]  Fidel Ramírez,et al.  deepTools: a flexible platform for exploring deep-sequencing data , 2014, Nucleic Acids Res..

[9]  N. Friedman,et al.  High-Resolution Sequencing and Modeling Identifies Distinct Dynamic RNA Regulatory Strategies , 2014, Cell.

[10]  Leighton J. Core,et al.  Nascent RNA Sequencing Reveals Widespread Pausing and Divergent Initiation at Human Promoters , 2008, Science.

[11]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[12]  B. Pugh,et al.  Comprehensive Genome-wide Protein-DNA Interactions Detected at Single-Nucleotide Resolution , 2011, Cell.

[13]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[14]  J. Weissman,et al.  Nascent transcript sequencing visualizes transcription at nucleotide resolution , 2011, Nature.

[15]  L. Steinmetz,et al.  The Nuclear PolyA-Binding Protein Nab2p Is Essential for mRNA Production. , 2015, Cell reports.

[16]  A. Bohm,et al.  Mechanism of poly(A) polymerase: structure of the enzyme-MgATP-RNA ternary complex and kinetic analysis. , 2007, Structure.

[17]  T. Jensen,et al.  A Nuclear Export Block Triggers the Decay of Newly Synthesized Polyadenylated RNA , 2018, Cell reports.

[18]  Nathan Morris,et al.  Codon Optimality Is a Major Determinant of mRNA Stability , 2015, Cell.

[19]  J. Corden,et al.  Regulation of yeast NRD1 expression by premature transcription termination. , 2006, Molecular cell.

[20]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[21]  G. Sanguinetti,et al.  Transcriptome-wide RNA processing kinetics revealed using extremely short 4tU labeling , 2015, Genome Biology.

[22]  Daniel Schulz,et al.  Transcriptome Surveillance by Selective Termination of Noncoding RNA Synthesis , 2013, Cell.

[23]  Richard J Maraia,et al.  Transcription termination by the eukaryotic RNA polymerase III. , 2013, Biochimica et biophysica acta.

[24]  J. Gagneur,et al.  TT-seq maps the human transient transcriptome , 2016, Science.

[25]  Leighton J. Core,et al.  Precise Maps of RNA Polymerase Reveal How Promoters Direct Initiation and Pausing , 2013, Science.

[26]  G. Chanfreau,et al.  Common genomic elements promote transcriptional and DNA replication roadblocks , 2016, Genome research.

[27]  Alex P. Reynolds,et al.  Native Elongating Transcript Sequencing Reveals Human Transcriptional Activity at Nucleotide Resolution , 2015, Cell.

[28]  Johannes Söding,et al.  Uniform transitions of the general RNA polymerase II transcription complex , 2010, Nature Structural &Molecular Biology.

[29]  Galt P. Barber,et al.  BigWig and BigBed: enabling browsing of large distributed datasets , 2010, Bioinform..

[30]  D. Gresham,et al.  Determination of in vivo RNA kinetics using RATE-seq , 2014, RNA.

[31]  B. Séraphin,et al.  Cryptic Pol II Transcripts Are Degraded by a Nuclear Quality Control Pathway Involving a New Poly(A) Polymerase , 2005, Cell.

[32]  David Tollervey,et al.  A Transcriptome-wide Atlas of RNP Composition Reveals Diverse Classes of mRNAs and lncRNAs , 2013, Cell.

[33]  G. Sanguinetti,et al.  Strand-specific , high-resolution mapping of modified RNA polymerase II , 2016 .

[34]  Achim Tresch,et al.  Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast , 2011, Molecular systems biology.

[35]  D. Gautheret,et al.  Nonsense-Mediated Decay Restricts LncRNA Levels in Yeast Unless Blocked by Double-Stranded RNA Structure , 2016, Molecular cell.

[36]  Srinivas Ramachandran,et al.  Transcription of Nearly All Yeast RNA Polymerase II-Transcribed Genes Is Dependent on Transcription Factor TFIID. , 2017, Molecular cell.

[37]  D. Libri,et al.  Transcription termination and the control of the transcriptome: why, where and how to stop , 2015, Nature Reviews Molecular Cell Biology.

[38]  M. V. Nogués,et al.  Bovine pancreatic ribonuclease: fifty years of the first enzymatic reaction mechanism. , 2011, Biochemistry.

[39]  Cizhong Jiang,et al.  A compiled and systematic reference map of nucleosome positions across the Saccharomyces cerevisiae genome , 2009, Genome Biology.

[40]  T. Jensen,et al.  Dealing with pervasive transcription. , 2013, Molecular cell.

[41]  U. K. Laemmli,et al.  The anchor-away technique: rapid, conditional establishment of yeast mutant phenotypes. , 2008, Molecular cell.

[42]  Roy Parker,et al.  RNA Degradation in Saccharomyces cerevisae , 2012, Genetics.

[43]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .

[44]  K. Helin,et al.  Gene silencing triggers polycomb repressive complex 2 recruitment to CpG islands genome wide. , 2014, Molecular cell.

[45]  R. Zimmer,et al.  High-resolution gene expression profiling for simultaneous kinetic parameter analysis of RNA synthesis and decay. , 2008, RNA.