Sequence determinants of polyadenylation-mediated regulation
Open accessAbstract
The cleavage and polyadenylation reaction is a crucial step in transcription termination and pre-mRNA maturation in human cells. Despite extensive research, the encoding of polyadenylation-mediated regulation of gene expression within the DNA sequence is not well understood. Here, we utilized a massively parallel reporter assay to inspect the effect of over 12,000 rationally designed polyadenylation sequences (PASs) on reporter gene expression and cleavage efficiency. We find that the PAS sequence can modulate gene expression by over five orders of magnitude. By using a uniquely designed scanning mutagenesis data set, we gain mechanistic insight into various modes of action by which the cleavage efficiency affects the sensitivity or robustness of the PAS to mutation. Furthermore, we employ motif discovery to identify both known and novel sequence motifs associated with PAS-mediated regulation. By leveraging the large scale of our data, we train a deep learning model for the highly accurate prediction of RNA levels from DNA sequence alone (R = 0.83). Moreover, we devise unique approaches for predicting exact cleavage sites for our reporter constructs and for endogenous transcripts. Taken together, our results expand our understanding of PAS-mediated regulation, and provide an unprecedented resource for analyzing and predicting PAS for regulatory genomics applications.
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Integrity signals
Registry
Cites
Statistics
Claims the paper makes
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Claims other papers make about this one
Each is checked against this paper’s own text: does it say what the citing paper says it does?
Sequence-based models of post-transcriptional regulation rely on genome annotations and transformed measurements extracted from RNA-seq to isolate each regulatory mechanism.
“Similarly, sequence-based models of post-transcriptional regulation rely on genome annotations and transformed measurements extracted from RNA-seq to isolate each regulatory mechanism (for example, percent spliced-in for splicing)14–23.”
From Predicting RNA-seq coverage from DNA sequence as a unifying model of gene regulation (2025) · cites this paper as [18]
The Discussion describes sequence-based models of polyadenylation-mediated regulation, and Methods > Computational analyses describes a transformed RNA-level measurement. But Results > High-throughput measurement of expression levels and cleavage maps says the models use measurements from a designed reporter library, not genome annotations and RNA-seq-derived measurements to isolate regulatory mechanisms. The paper does not address percent spliced-in or splicing. (LLM, conf=0.85)
Experimental — accuracy not yet measured
Self-check
Percentages, sums and reported differences re-computed from the paper’s own numbers.
The self-check needs the paper’s full text.
Cites
No references could be resolved yet.
Research works printed but not linked: 70
- Akhtar MN, Bukhari SA, Fazal Z, Qamar R, Shahmuradov IA. 2010 POLYAR, a new computer program for prediction of poly(A) sites in human sequences. BMC Genomics 11: 646 10.1186/1471-2164-11-646 21092114 PMC3053588 · DOI
- Alipanahi B, Delong A, Weirauch MT, Frey BJ. 2015 Predicting the sequence specificities of DNA- and RNA-binding proteins by deep learning. Nat Biotechnol 33: 831– 838. 10.1038/nbt.3300 26213851 · DOI · Paper page
- Angermueller C, Pärnamaa T, Parts L, Stegle O. 2016 Deep learning for computational biology. Mol Syst Biol 12: 878 10.15252/msb.20156651 27474269 PMC4965871 · DOI · Paper page
- Antoniou M, Geraghty F, Hurst J, Grosveld F. 1998 Efficient 3′-end formation of human β-globin mRNA in vivo requires sequences within the last intron but occurs independently of the splicing reaction. Nucleic Acids Res 26: 721– 729. 10.1093/nar/26.3.721 9443963 PMC147328 · DOI
- Bagga PS, Ford LP, Chen F, Wilusz J. 1995 The G-rich auxiliary downstream element has distinct sequence and position requirements and mediates efficient 3′ end pre-mRNA processing through a trans-acting factor. Nucleic Acids Res 23: 1625– 1631. 10.1093/nar/23.9.1625 7784220 PMC306907 · DOI
- Bailey TL. 2011 DREME: motif discovery in transcription factor ChIP-seq data. Bioinformatics 27: 1653– 1659. 10.1093/bioinformatics/btr261 21543442 PMC3106199 · DOI
- Bailey TL, Machanick P. 2012 Inferring direct DNA binding from ChIP-seq. Nucleic Acids Res 40: e128 10.1093/nar/gks433 22610855 PMC3458523 · DOI
- Bogard N, Linder J, Rosenberg AB, Seelig G. 2019 A deep neural network for predicting and engineering alternative polyadenylation. Cell 178: 91–106.e23. 10.1016/j.cell.2019.04.046 PMC6599575 31178116 · DOI · Paper page
- Bohnlein S, Hauber J, Cullen BR. 1989 Identification of a U5-specific sequence required for efficient polyadenylation within the human immunodeficiency virus long terminal repeat. J Virol 63: 421– 424. 2908926 PMC247699
- Brister JR, Ako-Adjei D, Bao Y, Blinkova O. 2015 NCBI Viral Genomes Resource. Nucleic Acids Res 43: D571– D577. 10.1093/nar/gku1207 25428358 PMC4383986 · DOI
- Chang T-H, Wu L-C, Chen Y-T, Huang H-D, Liu B-J, Cheng K-F, Horng J-T. 2011 Characterization and prediction of mRNA polyadenylation sites in human genes. Med Biol Eng Comput 49: 463– 472. 10.1007/s11517-011-0732-4 21286831 · DOI
- Chen T, Guestrin C. 2016 XGBoost: a scalable tree boosting system. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, pp. 785– 794. New York.
- Cheng Y, Miura RM, Tian B. 2006 Prediction of mRNA polyadenylation sites by support vector machine. Bioinformatics 22: 2320– 2325. 10.1093/bioinformatics/btl394 16870936 · DOI
- Connelly S, Manley JL. 1988 A functional mRNA polyadenylation signal is required for transcription termination by RNA polymerase II. Genes Dev 2: 440– 452. 10.1101/gad.2.4.440 2836265 · DOI
- Deng Z, Zhang S, Gu S, Ni X, Zeng W, Li X. 2018 Useful bicistronic reporter system for studying poly(A) site-defining cis elements and regulation of alternative polyadenylation. Int J Mol Sci 19: 279 10.3390/ijms19010279 PMC5796225 29342112 · DOI
- Gao X, Zhang J, Wei Z, Hakonarson H. 2018 DeepPolyA: a convolutional neural network approach for polyadenylation site prediction. IEEE Access 6: 24340– 24349. 10.1109/ACCESS.2018.2825996 · DOI
- Goodman DB, Church GM, Kosuri S. 2013 Causes and effects of N-terminal codon bias in bacterial genes. Science 342: 475– 479. 10.1126/science.1241934 24072823 · DOI · Paper page
- Goodwin EC, Rottman FM. 1992 The 3′-flanking sequence of the bovine growth hormone gene contains novel elements required for efficient and accurate polyadenylation. J Biol Chem 267: 16330– 16334. 1644817
- Graveley BR, Gilmartin GM. 1996 A common mechanism for the enhancement of mRNA 3′ processing by U3 sequences in two distantly related lentiviruses. J Virol 70: 1612– 1617. 8627681 PMC189984
- Hans H, Alwine JC. 2000 Functionally significant secondary structure of the simian virus 40 late polyadenylation signal. Mol Cell Biol 20: 2926– 2932. 10.1128/MCB.20.8.2926-2932.2000 10733596 PMC85533 · DOI
- Harrow J, Frankish A, Gonzalez JM, Tapanari E, Diekhans M, Kokocinski F, Aken BL, Barrell D, Zadissa A, Searle S, 2012 GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res 22: 1760– 1774. 10.1101/gr.135350.111 22955987 PMC3431492 · DOI · Paper page
- Hart RP, McDevitt MA, Ali H, Nevins JR. 1985 Definition of essential sequences and functional equivalence of elements downstream of the adenovirus E2A and the early simian virus 40 polyadenylation sites. Mol Cell Biol 5: 2975– 2983. 10.1128/MCB.5.11.2975 3018490 PMC369109 · DOI
- Hu, J, Lutz, CS, Wilusz, J, Tian, B. 2005 Bioinformatic identification of candidate cis-regulatory elements involved in human mRNA polyadenylation. RNA 11: 1485– 1493. 10.1261/rna.2107305 16131587 PMC1370832 · DOI
- Hulo C, de Castro E, Masson P, Bougueleret L, Bairoch A, Xenarios I, Le Mercier P. 2011 ViralZone: a knowledge resource to understand virus diversity. Nucleic Acids Res 39: D576– D582. 10.1093/nar/gkq901 20947564 PMC3013774 · DOI · Paper page
- Kelley DR, Snoek J, Rinn JL. 2016 Basset: learning the regulatory code of the accessible genome with deep convolutional neural networks. Genome Res 26: 990– 999. 10.1101/gr.200535.115 27197224 PMC4937568 · DOI · Paper page
- Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. 2007 The role of site accessibility in microRNA target recognition. Nat Genet 39: 1278– 1284. 10.1038/ng2135 17893677 · DOI · Paper page
- Kheradpour P, Ernst J, Melnikov A, Rogov P, Wang L, Zhang X, Alston J, Mikkelsen TS, Kellis M. 2013 Systematic dissection of regulatory motifs in 2000 predicted human enhancers using a massively parallel reporter assay. Genome Res 23: 800– 811. 10.1101/gr.144899.112 23512712 PMC3638136 · DOI · Paper page
- Legendre M, Gautheret D. 2003 Sequence determinants in human polyadenylation site selection. BMC Genomics 4: 7 10.1186/1471-2164-4-7 12600277 PMC151664 · DOI
- Leung MKK, Delong A, Frey BJ. 2018 Inference of the human polyadenylation code. Bioinformatics 34: 2889– 2898. 10.1093/bioinformatics/bty211 29648582 PMC6129302 · DOI · Paper page
- Levitt N, Briggs D, Gil A, Proudfoot NJ. 1989 Definition of an efficient synthetic poly(A) site. Genes Dev 3: 1019– 1025. 10.1101/gad.3.7.1019 2570734 · DOI
- Lin Y, Li Z, Ozsolak F, Kim SW, Arango-Argoty G, Liu TT, Tenenbaum SA, Bailey T, Monaghan AP, Milos PM, 2012 An in-depth map of polyadenylation sites in cancer. Nucleic Acids Res 40: 8460– 8471. 10.1093/nar/gks637 22753024 PMC3458571 · DOI
- Lorenz R, Bernhart SH, Höner zu Siederdissen C, Tafer H, Flamm C, Stadler PF, Hofacker IL. 2011 ViennaRNA Package 2.0. Algorithms Mol Biol 6: 26 10.1186/1748-7188-6-26 22115189 PMC3319429 · DOI · Paper page
- Lubliner S, Regev I, Lotan-Pompan M, Edelheit S, Weinberger A, Segal E. 2015 Core promoter sequence in yeast is a major determinant of expression level. Genome Res 25: 1008– 1117. 10.1101/gr.188193.114 25969468 PMC4484384 · DOI
- Lutz CS, Moreira A. 2011 Alternative mRNA polyadenylation in eukaryotes: an effective regulator of gene expression. Wiley Interdiscip Rev RNA 2: 22– 31. 10.1002/wrna.47 21956967 · DOI
- Magana-Mora A, Kalkatawi M, Bajic VB. 2017 Omni-PolyA: a method and tool for accurate recognition of Poly(A) signals in human genomic DNA. BMC Genomics 18: 620 10.1186/s12864-017-4033-7 28810905 PMC5558757 · DOI
- Marín RM, Vaníček J. 2011 Efficient use of accessibility in microRNA target prediction. Nucleic Acids Res 39: 19– 29. 10.1093/nar/gkq768 20805242 PMC3017612 · DOI
- Matoulkova E, Michalova E, Vojtesek B, Hrstka R. 2012 The role of the 3′ untranslated region in post-transcriptional regulation of protein expression in mammalian cells. RNA Biol 9: 563– 576. 10.4161/rna.20231 22614827 · DOI · Paper page
- McDevitt MA, Hart RP, Wong WW, Nevins JR. 1986 Sequences capable of restoring poly(A) site function define two distinct downstream elements. EMBO J 5: 2907– 2913. 10.1002/j.1460-2075.1986.tb04586.x 3024967 PMC1167241 · DOI
- McLeay RC, Bailey TL. 2010 Motif Enrichment Analysis: a unified framework and an evaluation on ChIP data. BMC Bioinformatics 11: 165 10.1186/1471-2105-11-165 20356413 PMC2868005 · DOI
- Mogno I, Kwasnieski JC, Cohen BA. 2013 Massively parallel synthetic promoter assays reveal the in vivo effects of binding site variants. Genome Res 23: 1908– 1915. 10.1101/gr.157891.113 23921661 PMC3814890 · DOI
- Moreira A, Wollerton M, Monks J, Proudfoot NJ. 1995 Upstream sequence elements enhance poly(A) site efficiency of the C2 complement gene and are phylogenetically conserved. EMBO J 14: 3809– 3819. 10.1002/j.1460-2075.1995.tb00050.x 7641699 PMC394455 · DOI
- Moreira A, Takagaki Y, Brackenridge S, Wollerton M, Manley JL, Proudfoot NJ. 1998 The upstream sequence element of the C2 complement poly(A) signal activates mRNA 3′ end formation by two distinct mechanisms. Genes Dev 12: 2522– 2534. 10.1101/gad.12.16.2522 9716405 PMC317083 · DOI
- Muerdter F, Boryń ŁM, Arnold CD. 2015 STARR-seq—Principles and applications. Genomics 106: 145– 150. 10.1016/j.ygeno.2015.06.001 26072434 · DOI
- Natalizio BJ. 2002 Upstream elements present in the 3′-untranslated region of collagen genes influence the processing efficiency of overlapping polyadenylation signals. J Biol Chem 277: 42733– 42740. 10.1074/jbc.M208070200 12200454 · DOI
- Noderer WL, Flockhart RJ, Bhaduri A, Diaz de Arce AJ, Zhang J, Khavari PA, Wang CL. 2014 Quantitative analysis of mammalian translation initiation sites by facs-seq. Mol Syst Biol 10: 748 10.15252/msb.20145136 25170020 PMC4299517 · DOI · Paper page
- Nunes NM, Li W, Tian B, Furger A. 2010 A functional human Poly(A) site requires only a potent DSE and an A-rich upstream sequence. EMBO J 29: 1523– 1536. 10.1038/emboj.2010.42 20339349 PMC2876958 · DOI
- Quang D, Xie X. 2016 DanQ: a hybrid convolutional and recurrent deep neural network for quantifying the function of DNA sequences. Nucleic Acids Res 44: e107 10.1093/nar/gkw226 27084946 PMC4914104 · DOI
- Rabani M, Pieper L, Chew G-L, Schier AF. 2017 A massively parallel reporter assay of 3′ UTR sequences identifies in vivo rules for mRNA degradation. Mol Cell 68: 1083– 1094.e5. 10.1016/j.molcel.2017.11.014 29225039 PMC5994907 · DOI
- Rosenberg AB, Patwardhan RP, Shendure J, Seelig G. 2015 Learning the sequence determinants of alternative splicing from millions of random sequences. Cell 163: 698– 711. 10.1016/j.cell.2015.09.054 26496609 · DOI
- Sadofsky M, Connelly S, Manley JL, Alwine JC. 1985 Identification of a sequence element on the 3′ side of AAUAAA which is necessary for simian virus 40 late mRNA 3′-end processing. Mol Cell Biol 5: 2713– 2719. 10.1128/MCB.5.10.2713 3016512 PMC367009 · DOI
- Schek N, Cooke C, Alwine JC. 1992 Definition of the upstream efficiency element of the simian virus 40 late polyadenylation signal by using in vitro analyses. Mol Cell Biol 12: 5386– 5393. 10.1128/MCB.12.12.5386 1333042 PMC360476 · DOI · Paper page
- Shalem O, Sharon E, Lubliner S, Regev I, Lotan-Pompan M, Yakhini Z, Segal E. 2015 Systematic dissection of the sequence determinants of gene 3′ end mediated expression control. PLoS Genet 11: e1005147 10.1371/journal.pgen.1005147 25875337 PMC4398552 · DOI
- Sharon E, Kalma Y, Sharp A, Raveh-Sadka T, Levo M, Zeevi D, Keren L, Yakhini Z, Weinberger A, Segal E. 2012 Inferring gene regulatory logic from high-throughput measurements of thousands of systematically designed promoters. Nat Biotechnol 30: 521– 530. 10.1038/nbt.2205 22609971 PMC3374032 · DOI · Paper page
- Sharon E, van Dijk D, Kalma Y, Keren L, Manor O, Yakhini Z, Segal E. 2014 Probing the effect of promoters on noise in gene expression using thousands of designed sequences. Genome Res 24: 1698– 1706. 10.1101/gr.168773.113 25030889 PMC4199362 · DOI
- Sheets MD, Ogg SC, Wickens MP. 1990 Point mutations in AAUAAA and the poly (A) addition site: effects on the accuracy and efficiency of cleavage and polyadenylation in vitro. Nucleic Acids Res 18: 5799– 5805. 10.1093/nar/18.19.5799 2170946 PMC332317 · DOI
- Shrikumar A, Greenside P, Kundaje A. 2017 Learning important features through propagating activation differences. arXiv: 1704.02685 [cs.CV]. · arXiv
- Sittler A, Gallinaro H, Jacob M. 1994 Upstream and downstream cis-acting elements for cleavage at the L4 polyadenylation site of adenovirus-2. Nucleic Acids Res 22: 222– 231. 10.1093/nar/22.2.222 8121807 PMC307775 · DOI
- Smith RP, Taher L, Patwardhan RP, Kim MJ, Inoue F, Shendure J, Ovcharenko I, Ahituv N. 2013 Massively parallel decoding of mammalian regulatory sequences supports a flexible organizational model. Nat Genet 45: 1021– 1028. 10.1038/ng.2713 23892608 PMC3775494 · DOI · Paper page
- Thomas LF, Saetrom P. 2012 Single nucleotide polymorphisms can create alternative polyadenylation signals and affect gene expression through loss of microRNA-regulation. PLoS Comput Biol 8: e1002621 10.1371/journal.pcbi.1002621 22915998 PMC3420919 · DOI
- Vainberg Slutskin I, Weingarten-Gabbay S, Nir R, Weinberger A, Segal E. 2018 Unraveling the determinants of microRNA mediated regulation using a massively parallel reporter assay. Nat Commun 9: 529 10.1038/s41467-018-02980-z 29410437 PMC5802814 · DOI
- Valsamakis A, Zeichner S, Carswell S, Alwine JC. 1991 The human immunodeficiency virus type 1 polyadenylylation signal: a 3′ long terminal repeat element upstream of the AAUAAA necessary for efficient polyadenylylation. Proc Natl Acad Sci 88: 2108– 2112. 10.1073/pnas.88.6.2108 1848693 PMC51178 · DOI
- Wang R, Zheng D, Yehia G, Tian B. 2018 A compendium of conserved cleavage and polyadenylation events in mammalian genes. Genome Res 28: 1427– 1441. 10.1101/gr.237826.118 30143597 PMC6169888 · DOI
- Weingarten-Gabbay S, Elias-Kirma S, Nir R, Gritsenko AA, Stern-Ginossar N, Yakhini Z, Weinberger A, Segal E. 2016 Systematic discovery of cap-independent translation sequences in human and viral genomes. Science 351: aad4939 10.1126/science.aad4939 26816383 · DOI · Paper page
- Wu X, Bartel DP. 2017 Widespread influence of 3′-end structures on mammalian mRNA processing and stability. Cell 169: 905– 917.e11. 10.1016/j.cell.2017.04.036 28525757 PMC5546744 · DOI
- Yoon OK, Hsu TY, Im JH, Brem RB. 2012 Genetics and regulatory impact of alternative polyadenylation in human B-lymphoblastoid cells. PLoS Genet 8: e1002882 10.1371/journal.pgen.1002882 22916029 PMC3420953 · DOI
- Zarudnaya MI, Kolomiets IM, Potyahaylo AL, Hovorun DM. 2003 Downstream elements of mammalian pre-mRNA polyadenylation signals: primary, secondary and higher-order structures. Nucleic Acids Res 31: 1375– 1386. 10.1093/nar/gkg241 12595544 PMC149834 · DOI
- Zhang F, Cole CN. 1987 Identification of a complex associated with processing and polyadenylation in vitro of herpes simplex virus type 1 thymidine kinase precursor RNA. Mol Cell Biol 7: 3277– 3286. 10.1128/MCB.7.9.3277 2823124 PMC367965 · DOI
- Zhang F, Denome RM, Cole CN. 1986 Fine-structure analysis of the processing and polyadenylation region of the herpes simplex virus type 1 thymidine kinase gene by using linker scanning, internal deletion, and insertion mutations. Mol Cell Biol 6: 4611– 4623. 10.1128/MCB.6.12.4611 2879221 PMC367246 · DOI
- Zhao J, Hyman L, Moore C. 1999 Formation of mRNA 3′ ends in eukaryotes: mechanism, regulation, and interrelationships with other steps in mRNA synthesis. Microbiol Mol Biol Rev 63: 405– 445. 10357856 PMC98971
- Zhou J, Troyanskaya OG. 2015 Predicting effects of noncoding variants with deep learning–based sequence model. Nat Methods 12: 931– 934. 10.1038/nmeth.3547 26301843 PMC4768299 · DOI · Paper page
Cited by
Sorted by citations
Predicting RNA-seq coverage from DNA sequence as a unifying model of gene regulation
Linder, J., Srivastava, D., Yuan, H., Agarwal, V., & Kelley, D. R. (2025) · Nature Genetics · doi:10.1038/s41588-024-02053-6 · Cited by 265 · see citing works
- Authors
- Vainberg Slutskin, I., Weinberger, A., & Segal, E.
- Affiliations
- Weizmann Institute of Science (IL)
- Corresponding author
- Segal, Eran
- Venue
- Genome Research
- Publisher
- Cold Spring Harbor Laboratory
- Published
- 2019-09-17
- Volume
- 29
- Issue
- 10
- Pages
- 1635-1647
- ISSN
- 1088-9051 (ISSN-L)
- 1549-5469
- PMC
- PMC6771402
- OpenAlex
- W2974101135
- Cited by
- 27
- References
- 0
- Citations by year
- 2026: 2 · 2025: 5 · 2024: 4 · 2023: 4 · 2022: 5 · 2021: 7
- Type
- posted content
- Language
- en
- Open access
- bronze
- doi.org
- Locations
- Genome Research (publishedVersion)
- Europe PMC (PubMed Central) (publishedVersion, cc-by-nc)
- PubMed Central (publishedVersion, cc-by-nc)
- PubMed (publishedVersion)
- HighWire Press Open Archive (submittedVersion)
- Funders
- European Research Council
- Israel Science Foundation
- Topics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA Research and Splicing
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Chromatin Dynamics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA and protein synthesis mechanisms
- Keywords
- polyadenylation, cleavage and polyadenylation, transcription termination, gene expression regulation, massively parallel reporter assay, scanning mutagenesis, sequence motifs, cleavage site prediction, regulatory genomics
- MeSH terms
- Deep Learning, Gene Expression Regulation, Humans, RNA Precursors, RNA, Messenger, Mutagenesis, Genes, Reporter, 3' Untranslated Regions, Genomics, Polyadenylation, High-Throughput Nucleotide Sequencing, Transcription Termination, Genetic
About these checks
Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.
Statistics and arithmetic are recomputed, not judged.
Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.
What counts as a claim ›How we measure confidence ›How we check this ›
Are you an author, or do you see a problem with this page? Request a review