Alternative polyadenylation (APA) is widely present in the human genome and plays a key role in carcinogenesis. We conducted a comprehensive analysis of the APA products in glioblastoma multiforme (GBM, one of the most lethal brain tumors) and normal brain tissues and further developed a computational pipeline, RNAelements (http://sysbio.zju.edu.cn/RNAelements/), using covariance model from known RNA binding protein (RBP) targets acquired by RNA Immunoprecipitation (RIP) analysis. We identified 4530 APA isoforms for 2733 genes in GBM, and found that 182 APA isoforms from 148 genes showed significant differential expression between normal and GBM brain tissues. We then focused on three genes with long and short APA isoforms that show inconsistent expression changes between normal and GBM brain tissues. These were myocyte enhancer factor 2D, heat shock factor binding protein 1, and polyhomeotic homolog 1 (Drosophila). Using the RNAelements program, we found that RBP binding sites were enriched in the alternative regions between the first and the last polyadenylation sites, which would result in the short APA forms escaping regulation from those RNA binding proteins. To the best of our knowledge, this report is the first comprehensive APA isoform dataset for GBM and normal brain tissues. Additionally, we demonstrated a putative novel APA-mediated mechanism for controlling RNA stability and translation for APA isoforms. These observations collectively lay a foundation for novel diagnostics and molecular mechanisms that can inform future therapeutic interventions for GBM.
Get full access to this article
View all access options for this article.
References
1.
AbdelmohsenK, GorospeM. 2010. Posttranscriptional regulation of cancer traits by HuR. Wiley interdisciplinary reviews. RNA, 1:214–1229.
2.
AkhtarMN, BukhariSA, FazalZ, QamarR, ShahmuradovIA. 2010. POLYAR, a new computer program for prediction of poly(A) sites in human sequences. BMC Genom, 11:646.
3.
AnJJ, GharamiK, LiaoGYet al.2008. Distinct role of long 3' UTR BDNF mRNA in spine morphology and synaptic plasticity in hippocampal neurons. Cell, 134:175–187.
AndreassiC, RiccioA. 2009. To localize or not to localize: mRNA fate is in 3'UTR ends. Trends Cell Biol, 19:465–474.
6.
BaileyTL, WilliamsN, MislehC, LiWW. 2006. MEME: Discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res, 34:W369–373.
7.
BandyopadhyayS, ChiangCY, SrivastavaJet al.2010. A human MAP kinase interactome. Nat Methods, 7:801–805.
8.
BaroniTE, ChitturSV, GeorgeAD, TenenbaumSA. 2008. Advances in RIP-chip analysis: RNA-binding protein immunoprecipitation-microarray profiling. Methods Mol Biol, 419:93–108.
9.
BrennerS, JohnsonM, BridghamJet al.2000. Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays. Nat Biotechnol, 18:630–634.
10.
BubenikJL, LaddAN, GerberCA, BudimanME, DriscollDM. 2009. Known turnover and translation regulatory RNA-binding proteins interact with the 3' UTR of SECIS-binding protein 2. RNA Biol, 6:73–83.
11.
ChabanonH, MickleburghI, HeskethJ. 2004. Zipcodes and postage stamps: mRNA localisation signals and their trans-acting binding proteins. Brief Funct Genom Proteom, 3:240–256.
12.
ChenYT, ScanlanMJ, VendittiCAet al.2005. Identification of cancer/testis-antigen genes by massively parallel signature sequencing. Proc Natl Acad Sci USA, 102:7940–7945.
13.
ChengY, MiuraRM, TianB. 2006. Prediction of mRNA polyadenylation sites by support vector machine. Bioinformatics, 22:2320–2325.
14.
De Sousa AbreuR, Sanchez-DiazPC, VogelCet al.2009. Genomic analyses of musashi1 downstream targets show a strong association with cancer-related processes. J Biol Chem, 284:12125–12135.
15.
EiringAM, NevianiP, SanthanamRet al.2008. Identification of novel posttranscriptional targets of the BCR/ABL oncoprotein by ribonomics: Requirement of E2F3 for BCR/ABL leukemogenesis. Blood, 111:816–828.
16.
FanJ, IshmaelFT, FangXet al.2011. Chemokine transcripts as targets of the RNA-binding protein HuR in human airway epithelium. J Immunol, 186:2482–2494.
FuY, SunY, LiY, LiJ, RaoX, ChenC, XuA. 2011. Differential genome-wide profiling of tandem 3' UTRs among human breast cancer and normal cells by high-throughput sequencing. Genome Res, 21:741–747.
19.
GalganoA, ForrerM, JaskiewiczL, KanitzA, ZavolanM, GerberAP. 2008. Comparative analysis of mRNA targets for human PUF-family proteins suggests extensive interaction with the miRNA regulatory system. PLoS ONE, 3:e3164.
20.
GoutsiasJ, LeeNH. 2007. Computational and experimental approaches for modeling gene regulatory networks. Curr Pharm Des, 13:1415–1436.
21.
GruberAR, NeubockR, HofackerIL, WashietlS. 2007. The RNAz web server: Prediction of thermodynamically stable and evolutionarily conserved RNA structures. Nucleic Acids Res, 35:W335–338.
22.
HafnerM, LandthalerM, BurgerLet al.2010. Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell, 141:129–141.
23.
HalbeisenRE, GalganoA, ScherrerT, GerberAP. 2008. Post-transcriptional gene regulation: From genome-wide studies to principles. Cell Mol Life Sci, 65:798–813.
24.
HavgaardJH, LyngsoRB, GorodkinJ. 2005. The FOLDALIGN web server for pairwise structural RNA alignment and mutual motif search. Nucleic Acids Res, 33:W650–653.
25.
HieronymusH, SilverPA. 2004. A systems view of mRNP biology. Genes Dev, 18:2845–2860.
26.
HillerM, PudimatR, BuschA, BackofenR. 2006. Using RNA secondary structures to guide sequence motif finding towards single-stranded regions. Nucleic Acids Res, 34:e117.
HofackerIL, PriwitzerB, StadlerPF. 2004. Prediction of locally stable RNA secondary structures for genome-wide surveys. Bioinformatics, 20:186–190.
30.
IrierHA, QuanY, YooJ, DingledineR. 2009a. Control of glutamate receptor 2 (GluR2) translational initiation by its alternative 3' untranslated regions. Mol Pharmacol, 76:1145–1149.
31.
IrierHA, ShawR, LauA, FengY, DingledineR.2009b. Translational regulation of GluR2 mRNAs in rat hippocampus by alternative 3' untranslated regions. J Neurochem, 109:584–594.
32.
JenalM, ElkonR, Loayza-PuchFet al.2012. The poly(A)-binding protein nuclear 1 suppresses alternative cleavage and polyadenylation sites. Cell, 149:538–553.
33.
JiY, XuX, StormoGD. 2004. A graph theoretical approach for predicting common RNA secondary structure motifs including pseudoknots in unaligned sequences. Bioinformatics, 20:1591–1602.
34.
JongeneelCV, DelorenziM, IseliCet al.2005. An atlas of human gene expression from massively parallel signature sequencing (MPSS)Genome Res, 15:1007–1014.
35.
JongeneelCV, IseliC, StevensonBJet al.2003. Comprehensive sampling of gene expression in human cell lines with massively parallel signature sequencing. Proc Natl Acad Sci USA, 100:4702–4705.
36.
KatoY, ZhaoM, MorikawaAet al.2000. Big mitogen-activated kinase regulates multiple members of the MEF2 protein family. J Biol Chem, 275:18534–18540.
37.
KazanH, RaD, ChanET, HughesTR, MorrisQ. 2010. RNAcontext: A new method for learning the sequence and structure binding preferences of RNA-binding proteins. PLoS Comput Biol, 6:e1000832.
KeeneJD, KomisarowJM, FriedersdorfMB. 2006. RIP-Chip: The isolation and identification of mRNAs, microRNAs and protein components of ribonucleoprotein complexes from cell extracts. Nat Protoc, 1:302–307.
KishoreS, Jaskiewiczl, BurgerL, HausserJ, KhorshidM, ZavolanM. 2011. A quantitative analysis of CLIP methods for identifying binding sites of RNA-binding proteins. Nat Methods, 8:559–564.
42.
KonigJ, ZarnackK, RotGet al.2010. iCLIP reveals the function of hnRNP particles in splicing at individual nucleotide resolution. Nat Struct Mol Biol, 17:909–915.
43.
LambertA, FontaineJF, LegendreMet al.2004. The ERPIN server: An interface to profile-based RNA motif identification. Nucleic Acids Res, 32:W160–165.
44.
LauAG, IrierHA, GuJet al.2010. Distinct 3'UTRs differentially regulate activity-dependent translation of brain-derived neurotrophic factor (BDNF)Proc Natl Acad Sci USA.
45.
LebedevaS, JensM, TheilKet al.2011. Transcriptome-wide analysis of regulatory interactions of the RNA-binding protein HuR. Mol Cell, 43:340–352.
46.
LinB, MadanA, YoonJGet al.2010. Massively parallel signature sequencing and bioinformatics analysis identifies up-regulation of TGFBI and SOX4 in human glioblastoma. PLoS ONE, 5:e10210.
47.
LinCH, PatelDJ. 1997. Structural basis of DNA folding and recognition in an AMP-DNA aptamer complex: Distinct architectures but common recognition motifs for DNA and RNA aptamers complexed to AMP. Chem Biol, 4:817–832.
48.
Lopez de SilanesI, ZhanM, LalA, YangX, GorospeM. 2004. Identification of a target RNA motif for RNA-binding protein HuR. Proc Natl Acad Sci USA, 101:2987–2992.
49.
MackeTJ, EckerDJ, GutellRR, GautheretD, CaseDA, SampathR. 2001. RNAMotif, an RNA secondary structure definition and search algorithm. Nucleic Acids Res, 29:4724–4735.
50.
MayrC, BartelDP. 2009. Widespread shortening of 3'UTRs by alternative cleavage and polyadenylation activates oncogenes in cancer cells. Cell, 138:673–684.
51.
Mazan-MamczarzK, KuwanoY, ZhanMet al.2009. Identification of a signature motif in target mRNAs of RNA-binding protein AUF1. Nucleic Acids Res, 37:204–214.
52.
MooreMJ. 2005. From birth to death: The complex lives of eukaryotic mRNAs. Science, 309:1514–1518.
53.
MukherjeeN, CorcoranDL, NusbaumJDet al.2011. Integrative regulatory mapping indicates that the RNA-binding protein HuR couples pre-mRNA processing and mRNA stability. Mol Cell, 43:327–339.
54.
MukherjeeN, LagerPJ, FriedersdorfMB, ThompsonMA, KeeneJD. 2009. Coordinated posttranscriptional mRNA population dynamics during T-cell activation. Mol Syst Biol, 5:288.
OtakeY, SoundararajanS, SenguptaTKet al.2007. Overexpression of nucleolin in chronic lymphocytic leukemia cells induces stabilization of bcl2 mRNA. Blood, 109:3069–3075.
57.
PavesiG, MauriG, StefaniM, PesoleG. 2004. RNAProfile: An algorithm for finding conserved secondary structure motifs in unaligned RNA sequences. Nucleic Acids Res, 32:3258–3269.
58.
RabaniM, KerteszM, SegalE. 2008. Computational prediction of RNA structural motifs involved in posttranscriptional regulatory processes. Proc Natl Acad Sci USA, 105:14885–14890.
59.
ReederJ, GiegerichR. 2007. Locomotif: From graphical motif description to RNA motif search. Bioinformatics, 23:i392–400.
60.
ReinartzJ, BruynsE, LinJZet al.2002. Massively parallel signature sequencing (MPSS) as a tool for in-depth quantitative gene expression profiling in all organisms. Brief Funct Genom Proteom, 1:95–104.
61.
Sanchez-BeatoM, SanchezE, Gonzalez-CarreroJet al.2006. Variability in the expression of polycomb proteins in different normal and tumoral tissues. A pilot study using tissue microarrays. Mod Pathol, 19:684–694.
62.
SandbergR, NeilsonJR, SarmaA, SharpPA, BurgeCB. 2008. Proliferating cells express mRNAs with shortened 3' untranslated regions and fewer microRNA target sites. Science, 320:1643–1647.
63.
SatyalSH, ChenD, FoxSG, KramerJM, MorimotoRI. 1998. Negative regulation of the heat shock transcriptional response by HSBP1. Genes Dev, 12:1962–1974.
64.
Scotto-LavinoR, DuG, FrohmanMA. 2006. 3' end cDNA amplification using classic RACE. Nat Protoc, 1:2742–2745.
65.
SiddharthanR, SiggiaED, Van NimwegenE. 2005. PhyloGibbs: A Gibbs sampling motif finder that incorporates phylogeny. PLoS Comput Biol, 1:e67.
66.
SinghP, AlleyTL, WrightSMet al.2009. Global changes in processing of mRNA 3' untranslated regions characterize clinically distinct cancer subtypes. Cancer Res, 69:9422–9430.
67.
SommerG, DittmannJ, KuehnertJet al.2011. The RNA-binding protein La contributes to cell proliferation and CCND1 expression. Oncogene, 30:434–444.
68.
TianB, HuJ, ZhangH, LutzCS. 2005. A large-scale analysis of mRNA polyadenylation of human and mouse genes. Nucleic Acids Res, 33:201–212.
69.
TollerveyJR, CurkT, RogelJBet al.2011. Characterizing the RNA targets and position-dependent splicing regulation by TDP-43. Nat Neurosci, 14:452–458.
70.
UleJ, JensenKB, RuggiuM, MeleA, UleA, DarnellRB. 2003. CLIP identifies Nova-regulated RNA networks in the brain. Science, 302:1212–1215.
71.
WangZ, GersteinM, SnyderM. 2009a. RNA-Seq: A revolutionary tool for transcriptomics. Nat Rev Genet, 10:57–63.
72.
WangZ, KayikciM, BrieseMet al.2010. iCLIP predicts the dual splicing effects of TIA-RNA interactions. PLoS Biol, 8:e1000530.
73.
WangZ, TollerveyJ, BrieseM, TurnerF, UleJ. 2009b. CLIP: Construction of cDNA libraries for high-throughput sequencing from RNAs cross-linked to proteins in vivo. Methods, 48:287–293.
74.
WeinbergZ, BreakerRR. 2011. R2R–software to speed the depiction of aesthetic consensus RNA secondary structures. BMC Bioinform, 12:3.
75.
XuX, JiY, StormoGD. 2007. RNA Sampler: A new sampling based algorithm for common RNA secondary structure prediction and structural alignment. Bioinformatics, 23:1883–1891.
76.
YangC, MaiguelDA, CarrierF. 2002. Identification of nucleolin and nucleophosmin as genotoxic stress-responsive RNA-binding proteins. Nucleic Acids Res, 30:2251–2260.
77.
YangCC, OrnatskyOI, McDermottJC, CruzTF, ProdyCA. 1998. Interaction of myocyte enhancer factor 2 (MEF2) with a mitogen-activated protein kinase, ERK5/BMK1. Nucleic Acids Res, 26:4771–4777.
78.
YaoZ, WeinbergZ, RuzzoWL. 2006. CMfinder— A covariance model based RNA motif finding algorithm. Bioinformatics, 22:445–452.
79.
ZhangH, HuJ, RecceM, TianB. 2005a. PolyA_DB: A database for mammalian mRNA polyadenylation. Nucleic Acids Res, 33:D116–120.
80.
ZhangH, LeeJY, TianB. 2005b. Biased alternative polyadenylation in human tissues. Genome Biol, 6:R100.
81.
ZhangJ, TsaprailisG, BowdenGT. 2008. Nucleolin stabilizes Bcl-X L messenger RNA in response to UVA irradiation. Cancer Res, 68:1046–1054.
82.
ZlotorynskiE, AgamiR. 2008. A PASport to cellular proliferation. Cell, 134:208–210.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
