Discovery of physiological and cancer-related regulators of 3′ UTR processing with KAPAC

3′ Untranslated regions (3' UTRs) length is regulated in relation to cellular state. To uncover key regulators of poly(A) site use in specific conditions, we have developed PAQR, a method for quantifying poly(A) site use from RNA sequencing data and KAPAC, an approach that infers activities of oligomeric sequence motifs on poly(A) site choice. Application of PAQR and KAPAC to RNA sequencing data from normal and tumor tissue samples uncovers motifs that can explain changes in cleavage and polyadenylation in specific cancers. In particular, our analysis points to polypyrimidine tract binding protein 1 as a regulator of poly(A) site choice in glioblastoma.

[1]  M. Zavolan,et al.  3' End Sequencing Library Preparation with A-seq2 , 2017, Journal of visualized experiments : JoVE.

[2]  J. Ule,et al.  High-Resolution RNA Maps Suggest Common Principles of Splicing and Polyadenylation Regulation by TDP-43 , 2017, Cell reports.

[3]  B. Tian,et al.  Alternative polyadenylation of mRNA precursors , 2016, Nature Reviews Molecular Cell Biology.

[4]  Gene W. Yeo,et al.  Robust transcriptome-wide discovery of RNA binding protein binding sites with enhanced CLIP (eCLIP) , 2016, Nature Methods.

[5]  Liewei Wang,et al.  Measure transcript integrity using RNA-seq data , 2016, BMC Bioinformatics.

[6]  Ralf Schmidt,et al.  A comprehensive analysis of 3′ end sequencing data sets reveals novel polyadenylation signals and the repressive role of heterogeneous ribonucleoprotein C on cleavage and polyadenylation , 2015, bioRxiv.

[7]  Wei Li,et al.  NUDT21-spanning CNVs lead to neuropsychiatric disease and altered MeCP2 abundance via alternative polyadenylation , 2015, eLife.

[8]  Anne-Claude Gingras,et al.  An alternative splicing event amplifies evolutionary differences between vertebrates , 2015, Science.

[9]  Auinash Kalsotra,et al.  Systematic Profiling of Poly(A)+ Transcripts Modulated by Core 3’ End Processing and Splicing Factors Reveals Regulatory Rules of Alternative Cleavage and Polyadenylation , 2015, PLoS genetics.

[10]  Hyunmin Kim,et al.  Pre-mRNA splicing is facilitated by an optimal RNA polymerase II elongation rate , 2014, Genes & development.

[11]  Wei Li,et al.  Dynamic analyses of alternative polyadenylation from RNA-seq reveal a 3′-UTR landscape across seven tumour types , 2014, Nature Communications.

[12]  John R Yates,et al.  CPSF30 and Wdr33 directly bind to AAUAAA in mammalian mRNA 3′ processing , 2014, Genes & development.

[13]  M. Zavolan,et al.  Reconstitution of CPSF active in polyadenylation: recognition of the polyadenylation signal by WDR33 , 2014, Genes & development.

[14]  Paul Theodor Pyl,et al.  HTSeq – A Python framework to work with high-throughput sequencing data , 2014, bioRxiv.

[15]  Wei Li,et al.  CFIm25 links Alternative Polyadenylation to Glioblastoma Tumor Suppression , 2014, Nature.

[16]  Piotr J. Balwierz,et al.  ISMARA: automated modeling of genomic signals as a democracy of regulatory motifs , 2014, Genome research.

[17]  J. Chen,et al.  Alternative splicing in cancer: implications for biology and therapy , 2014, Oncogene.

[18]  M. Zavolan,et al.  Means to an end: mechanisms of alternative polyadenylation of messenger RNA precursors , 2013, Wiley interdisciplinary reviews. RNA.

[19]  Julie L. Yang,et al.  Ubiquitously transcribed genes use alternative polyadenylation to achieve tissue-specific expression , 2013, Genes & development.

[20]  S. Liebhaber,et al.  αCP Poly(C) Binding Proteins Act as Global Regulators of Alternative Polyadenylation , 2013, Molecular and Cellular Biology.

[21]  Julian König,et al.  Direct Competition between hnRNP C and U2AF65 Protects the Transcriptome from the Exonization of Alu Elements , 2013, Cell.

[22]  Andreas R. Gruber,et al.  Cleavage factor Im is a key regulator of 3′ UTR length , 2012, RNA biology.

[23]  Ji Wan,et al.  Transcriptome-wide analyses of CstF64–RNA interactions in global regulation of mRNA alternative polyadenylation , 2012, Proceedings of the National Academy of Sciences.

[24]  N. Goshima,et al.  Alternative 3′‐end processing of long noncoding RNA initiates construction of nuclear paraspeckles , 2012, The EMBO journal.

[25]  Bronwen L. Aken,et al.  GENCODE: The reference human genome annotation for The ENCODE Project , 2012, Genome research.

[26]  Wei Li,et al.  RSeQC: quality control of RNA-seq experiments , 2012, Bioinform..

[27]  Raymond K. Auerbach,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[28]  Larry N. Singh,et al.  U1 snRNP Determines mRNA Length and Regulates Isoform Expression , 2012, Cell.

[29]  Mihaela Zavolan,et al.  Genome-wide analysis of pre-mRNA 3' end processing reveals a decisive role of human cleavage factor I in the regulation of 3' UTR length. , 2012, Cell reports.

[30]  Douglas L Black,et al.  Neuronal regulation of pre-mRNA splicing by polypyrimidine tract binding proteins, PTBP1 and PTBP2 , 2012, Critical reviews in biochemistry and molecular biology.

[31]  R. Elkon,et al.  The Poly(A)-Binding Protein Nuclear 1 Suppresses Alternative Cleavage and Polyadenylation Sites , 2012, Cell.

[32]  N. Proudfoot.  Ending the message: poly(A) signals then and now. , 2011, Genes & development.

[33]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[34]  Larry N. Singh,et al.  U1 snRNP protects pre-mRNAs from premature cleavage and polyadenylation , 2010, Nature.

[35]  J. Ule,et al.  iCLIP reveals the function of hnRNP particles in splicing at individual nucleotide resolution , 2010, Nature Structural &Molecular Biology.

[36]  Peter F. Stadler,et al.  Fast Mapping of Short Sequences with Mismatches, Insertions and Deletions Using Index Structures , 2009, PLoS Comput. Biol..

[37]  Wen Zhu,et al.  Splicing factors PTBP1 and PTBP2 promote proliferation and migration of glioma cell lines. , 2009, Brain : a journal of neurology.

[38]  B. Tian,et al.  Progressive lengthening of 3′ untranslated regions of mRNAs by alternative polyadenylation during mouse embryonic development (Proceedings of the National Academy of Sciences of the United States of America (2009) 106, 17, (7028-7033) DOI 10.1073/pnas.0900028106) , 2009 .

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

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

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

[42]  G. Varani,et al.  RNA-binding proteins: modular design for efficient function , 2007, Nature Reviews Molecular Cell Biology.

[43]  H. Handa,et al.  Knock-down of 25 kDa subunit of cleavage factor Im in Hela cells alters alternative polyadenylation within 3′-UTRs , 2006, Nucleic acids research.

[44]  B. Blencowe,et al.  An RNA map predicting Nova-dependent splicing regulation , 2006, Nature.

[45]  T. Hampton,et al.  The Cancer Genome Atlas , 2020, Indian Journal of Medical and Paediatric Oncology.

[46]  W. Keller,et al.  An interaction between U2AF 65 and CF Im links the splicing and 3′ end processing machineries , 2006, The EMBO journal.

[47]  Alexandre V. Morozov,et al.  Statistical mechanical modeling of genome-wide transcription factor occupancy data by MatrixREDUCE , 2006, ISMB.

[48]  M. Wollerton,et al.  Polypyrimidine Tract Binding Protein Modulates Efficiency of Polyadenylation , 2004, Molecular and Cellular Biology.

[49]  E. Wahle,et al.  3'-End processing of pre-mRNA in eukaryotes. , 1999, FEMS microbiology reviews.

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

[51]  Sven Rahmann,et al.  Genome analysis , 2022 .

[52]  K. Baggerly,et al.  Splicing factors PTBP 1 and PTBP 2 promote proliferation and migration of glioma cell lines , 2009 .

[53]  Thomas Ragg,et al.  The RIN: an RNA integrity number for assigning integrity values to RNA measurements , 2006, BMC Molecular Biology.

[54]  Chin H. Li,et al.  The Ucsc Genome Browser Database: Update 2011 , 2022 .