An efficient method for genome-wide polyadenylation site mapping and RNA quantification

The use of alternative poly(A) sites is common and affects the post-transcriptional fate of mRNA, including its stability, subcellular localization and translation. Here, we present a method to identify poly(A) sites in a genome-wide and strand-specific manner. This method, termed 3′T-fill, initially fills in the poly(A) stretch with unlabeled dTTPs, allowing sequencing to start directly after the poly(A) tail into the 3′-untranslated regions (UTR). Our comparative analysis demonstrates that it outperforms existing protocols in quality and throughput and accurately quantifies RNA levels as only one read is produced from each transcript. We use this method to characterize the diversity of polyadenylation in Saccharomyces cerevisiae, showing that alternative RNA molecules are present even in a genetically identical cell population. Finally, we observe that overlap of convergent 3′-UTRs is frequent but sharply limited by coding regions, suggesting factors that restrict compression of the yeast genome.

[1]  R. Elkon,et al.  E2F mediates enhanced alternative polyadenylation in proliferation , 2012, Genome Biology.

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

[3]  Lira Mamanova,et al.  Low-bias, strand-specific transcriptome Illumina sequencing by on-flowcell reverse transcription (FRT-seq) , 2011, Nature Protocols.

[4]  Xiang-Dong Fu,et al.  A multiplex RNA-seq strategy to profile poly(A+) RNA: application to analysis of transcription response and 3' end formation. , 2011, Genomics.

[5]  C. Lutz,et al.  Alternative mRNA polyadenylation in eukaryotes: an effective regulator of gene expression , 2011, Wiley interdisciplinary reviews. RNA.

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

[7]  P. Brown,et al.  Identification of RNA recognition elements in the Saccharomyces cerevisiae transcriptome , 2010, Nucleic acids research.

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

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

[10]  Helga Thorvaldsdóttir,et al.  Integrative Genomics Viewer , 2011, Nature Biotechnology.

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

[12]  Serban Nacu,et al.  Fast and SNP-tolerant detection of complex variants and splicing in short reads , 2010, Bioinform..

[13]  B. Tian,et al.  Progressive lengthening of 3′ untranslated regions of mRNAs by alternative polyadenylation during mouse embryonic development , 2009, Proceedings of the National Academy of Sciences.

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

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

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

[17]  Tyson A. Clark,et al.  HITS-CLIP yields genome-wide insights into brain alternative RNA processing , 2008, Nature.

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

[19]  W. Huber,et al.  which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. MAnorm: a robust model for quantitative comparison of ChIP-Seq data sets , 2011 .

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

[21]  T. Hashimshony,et al.  CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification. , 2012, Cell reports.

[22]  S. Kuersten,et al.  The power of the 3′ UTR: translational control and development , 2003, Nature Reviews Genetics.

[23]  Vicent Pelechano,et al.  Genome-wide polyadenylation site mapping. , 2012, Methods in enzymology.

[24]  G. Azzam,et al.  Developmental RNA processing of 3′UTRs in Hox mRNAs as a context-dependent mechanism modulating visibility to microRNAs , 2010, Development.

[25]  N. Woo,et al.  Distinct Role of Long 3′ UTR BDNF mRNA in Spine Morphology and Synaptic Plasticity in Hippocampal Neurons , 2008, Cell.

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

[27]  T. Borodina,et al.  Transcriptome analysis by strand-specific sequencing of complementary DNA , 2009, Nucleic acids research.

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

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

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

[31]  W. Huber,et al.  Detecting differential usage of exons from RNA-seq data , 2012, Genome research.

[32]  Bin Tian,et al.  A large-scale analysis of mRNA polyadenylation of human and mouse genes , 2005, Nucleic acids research.

[33]  J. Corden,et al.  Termination of cryptic unstable transcripts is directed by yeast RNA-binding proteins Nrd1 and Nab3. , 2006, Molecular cell.