PAT-seq: a method to study the integration of 3′-UTR dynamics with gene expression in the eukaryotic transcriptome

A major objective of systems biology is to quantitatively integrate multiple parameters from genome-wide measurements. To integrate gene expression with dynamics in poly(A) tail length and adenylation site, we developed a targeted next-generation sequencing approach, Poly(A)-Test RNA-sequencing. PAT-seq returns (i) digital gene expression, (ii) polyadenylation site/s, and (iii) the polyadenylation-state within and between eukaryotic transcriptomes. PAT-seq differs from previous 3' focused RNA-seq methods in that it depends strictly on 3' adenylation within total RNA samples and that the full-native poly(A) tail is included in the sequencing libraries. Here, total RNA samples from budding yeast cells were analyzed to identify the intersect between adenylation state and gene expression in response to loss of the major cytoplasmic deadenylase Ccr4. Furthermore, concordant changes to gene expression and adenylation-state were demonstrated in the classic Crabtree-Warburg metabolic shift. Because all polyadenylated RNA is interrogated by the approach, alternative adenylation sites, noncoding RNA and RNA-decay intermediates were also identified. Most important, the PAT-seq approach uses standard sequencing procedures, supports significant multiplexing, and thus replication and rigorous statistical analyses can for the first time be brought to the measure of 3'-UTR dynamics genome wide.

[1]  M. Hengartner,et al.  Analysis of C. elegans intestinal gene expression and polyadenylation by fluorescence-activated nuclei sorting and 3′-end-seq , 2012, Nucleic acids research.

[2]  R Parker,et al.  A turnover pathway for both stable and unstable mRNAs in yeast: evidence for a requirement for deadenylation. , 1993, Genes & development.

[3]  K. Nishida,et al.  Mechanisms and consequences of alternative polyadenylation. , 2011, Molecules and Cells.

[4]  Helga Thorvaldsdóttir,et al.  Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration , 2012, Briefings Bioinform..

[5]  Job Harms,et al.  THE LANDSCAPE OF , 2010 .

[6]  John B. VanCuylenberg,et al.  ePAT: a simple method to tag adenylated RNA to measure poly(A)-tail length and other 3' RACE applications. , 2012, RNA.

[7]  P. Kapranov,et al.  Comprehensive Polyadenylation Site Maps in Yeast and Human Reveal Pervasive Alternative Polyadenylation , 2010, Cell.

[8]  J. Yates,et al.  Mammalian miRNA RISC recruits CAF1 and PABP to affect PABP-dependent deadenylation. , 2009, Molecular cell.

[9]  Clifford A. Meyer,et al.  Gene expression profiling of human breast tissue samples using SAGE-Seq. , 2010, Genome research.

[10]  B. Tian,et al.  Reprogramming of 3′ Untranslated Regions of mRNAs by Alternative Polyadenylation in Generation of Pluripotent Stem Cells from Different Cell Types , 2009, PloS one.

[11]  I. Nookaew,et al.  A comprehensive comparison of RNA-Seq-based transcriptome analysis from reads to differential gene expression and cross-comparison with microarrays: a case study in Saccharomyces cerevisiae , 2012, Nucleic acids research.

[12]  C. Mayr,et al.  Widespread Shortening of 3′UTRs by Alternative Cleavage and Polyadenylation Activates Oncogenes in Cancer Cells , 2009, Cell.

[13]  M. Robinson,et al.  A scaling normalization method for differential expression analysis of RNA-seq data , 2010, Genome Biology.

[14]  T. Beilharz,et al.  Using Klenow-mediated extension to measure poly(A)-tail length and position in the transcriptome. , 2014, Methods in molecular biology.

[15]  P. Sharp,et al.  Proliferating Cells Express mRNAs with Shortened 3' Untranslated Regions and Fewer MicroRNA Target Sites , 2008, Science.

[16]  L. Steinmetz,et al.  Functional consequences of bidirectional promoters. , 2011, Trends in genetics : TIG.

[17]  Peter J. Shepard,et al.  Complex and dynamic landscape of RNA polyadenylation revealed by PAS-Seq. , 2011, RNA.

[18]  Seth M. Kelly,et al.  Mutation of the conserved polyadenosine RNA binding protein, ZC3H14/dNab2, impairs neural function in Drosophila and humans , 2011, Proceedings of the National Academy of Sciences.

[19]  Ji Huang,et al.  [Serial analysis of gene expression]. , 2002, Yi chuan = Hereditas.

[20]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[21]  M. Swanson,et al.  Control of cleavage site selection during mRNA 3′ end formation by a yeast hnRNP , 1998, The EMBO journal.

[22]  J. Mata,et al.  A Network of Multiple Regulatory Layers Shapes Gene Expression in Fission Yeast , 2007, Molecular cell.

[23]  P. Cramer,et al.  Global analysis of eukaryotic mRNA degradation reveals Xrn1-dependent buffering of transcript levels. , 2013, Molecular cell.

[24]  A. Sachs,et al.  Poly(A) Tail Length Control in Saccharomyces cerevisiae Occurs by Message-Specific Deadenylation , 1998, Molecular and Cellular Biology.

[25]  D. Bartel,et al.  Extensive alternative polyadenylation during zebrafish development , 2012, Genome research.

[26]  D. Bartel,et al.  Formation, Regulation and Evolution of Caenorhabditis elegans 3′UTRs , 2010, Nature.

[27]  David I. K. Martin,et al.  microRNA-Mediated Messenger RNA Deadenylation Contributes to Translational Repression in Mammalian Cells , 2009, PloS one.

[28]  T. Babak,et al.  A quantitative atlas of polyadenylation in five mammals , 2012, Genome research.

[29]  Vicent Pelechano,et al.  An efficient method for genome-wide polyadenylation site mapping and RNA quantification , 2013, Nucleic acids research.

[30]  D. Riddle,et al.  Deep SAGE analysis of the Caenorhabditis elegans transcriptome , 2010, Nucleic acids research.

[31]  Elisa Izaurralde,et al.  Deadenylation is a widespread effect of miRNA regulation. , 2008, RNA.

[32]  R. Parker,et al.  The Transcription Factor Associated Ccr4 and Caf1 Proteins Are Components of the Major Cytoplasmic mRNA Deadenylase in Saccharomyces cerevisiae , 2001, Cell.

[33]  F. Real,et al.  Key contribution of CPEB4-mediated translational control to cancer progression , 2011, Nature Medicine.

[34]  T. Hughes,et al.  Genome-Wide Analysis of mRNA Stability Using Transcription Inhibitors and Microarrays Reveals Posttranscriptional Control of Ribosome Biogenesis Factors , 2004, Molecular and Cellular Biology.

[35]  M. Gerstein,et al.  The Transcriptional Landscape of the Yeast Genome Defined by RNA Sequencing , 2008, Science.

[36]  Frederick P. Roth,et al.  Next generation software for functional trend analysis , 2009, Bioinform..

[37]  Sylvain Foissac,et al.  New class of gene-termini-associated human RNAs suggests a novel RNA copying mechanism , 2010, Nature.

[38]  T. Preiss,et al.  Transcriptome-wide measurement of mRNA polyadenylation state. , 2009, Methods.

[39]  E. Heinz,et al.  The cellular roles of Ccr4-NOT in model and pathogenic fungi—implications for fungal virulence , 2013, Front. Genet..

[40]  Andrew H. Beck,et al.  3′-End Sequencing for Expression Quantification (3SEQ) from Archival Tumor Samples , 2010, PloS one.

[41]  M. Rigoulet,et al.  The Warburg and Crabtree effects: On the origin of cancer cell energy metabolism and of yeast glucose repression. , 2011, Biochimica et biophysica acta.

[42]  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.

[43]  Shimyn Slomovic,et al.  Addition of poly(A) and poly(A)-rich tails during RNA degradation in the cytoplasm of human cells , 2010, Proceedings of the National Academy of Sciences.

[44]  G. Barsh,et al.  Digital gene expression for non-model organisms. , 2011, Genome research.

[45]  D. Bartel,et al.  Poly(A)-tail profiling reveals an embryonic switch in translational control , 2014, Nature.

[46]  Sebastian D. Mackowiak,et al.  The Landscape of C. elegans 3′UTRs , 2010, Science.

[47]  Chong-Jian Chen,et al.  Differential genome-wide profiling of tandem 3' UTRs among human breast cancer and normal cells by high-throughput sequencing. , 2011, Genome research.

[48]  T. Preiss,et al.  Widespread use of poly(A) tail length control to accentuate expression of the yeast transcriptome. , 2007, RNA.

[49]  Gordon K Smyth,et al.  Statistical Applications in Genetics and Molecular Biology Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments , 2011 .

[50]  R. Brem,et al.  Noncanonical transcript forms in yeast and their regulation during environmental stress. , 2010, RNA.

[51]  Douglas M. Hawkins,et al.  A variance-stabilizing transformation for gene-expression microarray data , 2002, ISMB.