RNA editing in the human ENCODE RNA-seq data
Eddie Park, Brian Williams, Barbara J Wold, Ali Mortazavi, Eddie Park, Brian Williams, Barbara J Wold, Ali Mortazavi
Abstract
RNA-seq data can be mined for sequence differences relative to the reference genome to identify both genomic SNPs and RNA editing events. We analyzed the long, polyA-selected, unstranded, deeply sequenced RNA-seq data from the ENCODE Project across 14 human cell lines for candidate RNA editing events. On average, 43% of the RNA sequencing variants that are not in dbSNP and are within gene boundaries are A-to-G(I) RNA editing candidates. The vast majority of A-to-G(I) edits are located in introns and 3' UTRs, with only 123 located in protein-coding sequence. In contrast, the majority of non-A-to-G variants (60%-80%) map near exon boundaries and have the characteristics of splice-mapping artifacts. After filtering out all candidates with evidence of private genomic variation using genome resequencing or ChIP-seq data, we find that up to 85% of the high-confidence RNA variants are A-to-G(I) editing candidates. Genes with A-to-G(I) edits are enriched in Gene Ontology terms involving cell division, viral defense, and translation. The distribution and character of the remaining non-A-to-G variants closely resemble known SNPs. We find no reproducible A-to-G(I) edits that result in nonsynonymous substitutions in all three lymphoblastoid cell lines in our study, unlike RNA editing in the brain. Given that only a fraction of sites are reproducibly edited in multiple cell lines and that we find a stronger association of editing and specific genes suggests that the editing of the transcript is more important than the editing of any individual site.
Figures
References
- The 1000 Genomes Project Consortium 2010. A map of human genome variation from population-scale sequencing. Nature 467: 1061–1073
- Agranat L, Raitskin O, Sperling J, Sperling R 2008. The editing enzyme ADAR and the mRNA surveillance protein hUpf1 interact in the cell nucleus. Proc Natl Acad Sci 105: 5028–5033
- Athanasiadis A, Rich A, Maas S 2004. Widespread A-to-I RNA editing of Alu-containing mRNAs in the human transcriptome. PLoS Biol 2: e391 doi: 10.1371/journal.pbio.0020391
- Bahn JH, Lee JH, Li G, Greer C, Peng G, Xiao X 2012. Accurate identification of A-to-I RNA editing in human by transcriptome sequencing. Genome Res 22: 142–150
- Burns CM, Chu H, Rueter SM, Hutchinson LK, Canton H, Sanders-Bush E, Emeson RB 1997. Regulation of serotonin-2C receptor G-protein coupling by RNA editing. Nature 387: 303–308
- Cenci C, Barzotti R, Galeano F, Corbelli S, Rota R, Massimi L, Di Rocco C, O'Connell MA, Gallo A 2008. Down-regulation of RNA editing in pediatric astrocytomas: ADARB1 editing activity inhibits cell migration and proliferation. J Biol Chem 283: 7251–7260
- de Hoon MJL, Imoto S, Nolan J, Miyano S 2004. Open Source Clustering Software. Bioinformatics 20: 1453–1454
- Doria M, Neri F, Gallo A, Farace MG, Michienzi A 2009. Editing of HIV-1 RNA by the double-stranded RNA deaminase ADAR stimulates viral infection. Nucleic Acids Res 37: 5848–5858
- Gerber AP, Keller W 2001. RNA editing by base deamination: more enzymes, more targets, new mysteries. Trends Biochem Sci 26: 376–384
- Iizasa H, Nishikura K 2009. A new function for the RNA-editing enzyme ADAR. Nat Immunol 10: 16–18
- Kim U, Wang Y, Sanford T, Zeng Y, Nishikura K 1994. Molecular cloning of cDNA for double-stranded RNA adenosine deaminase, a candidate enzyme for nuclear RNA editing. Proc Natl Acad Sci 91: 11457–11461
- Kim DD, Kim TT, Walsh T, Kobayashi Y, Matise TC, Buyske S, Gabriel A 2004. Widespread RNA editing of embedded Alu elements in the human transcriptome. Genome Res 14: 1719–1725
- Kiran A, Baranov PV 2010. DARNED: A DAtabase of RNa EDiting in humans. Bioinformatics 26: 772–776
- Kleinman CL, Majewski J 2012. Comment on “Widespread RNA and DNA sequence differences in the human transcriptome.” Science 335: 1302 doi: 10.1126/science.1209658
- Kumar M, Carmichael GC 1997. Nuclear antisense RNA induces extensive adenosine modifications and nuclear retention of target transcripts. Proc Natl Acad Sci 94: 3542–3547
- Langmead B, Trapnell C, Pop M, Salzberg SL 2009. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25 doi: 10.1186/gb-2009-10-3-r25
- Li M, Wang IX, Li Y, Bruzel A, Richards AL, Toung JM, Cheung VG 2011. Widespread RNA and DNA sequence differences in the human transcriptome. Science 333: 53–58
- Lin W, Piskol R, Tan MH, Li JB 2012. Comment on “Widespread RNA and DNA sequence differences in the human transcriptome.” Science 335: 1302 doi: 10.1126/science.1210624
- Luciano DJ, Mirsky H, Vendetti NJ, Maas S 2004. RNA editing of a miRNA precursor. RNA 10: 1174–1177
- Maas S, Kawahara Y, Tamburro KM, Nishikura K 2006. A-to-I RNA editing and human disease. RNA Biol 3: 1–9
- McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G 2010. GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol 28: 495–501
- Montgomery SB, Sammeth M, Gutierrez-Arcelus M, Larch RP, Ingle C, Nisbett J, Guigo R, Dermitzakis ET 2010. Transcriptome genetics using second generation sequencing in a Caucasian population. Nature 464: 773–777
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B 2008. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5: 621–628
- Nishikura K 2006. Editor meets silencer: Crosstalk between RNA editing and RNA interference. Nat Rev Mol Cell Biol 7: 919–931
- Nishikura K 2010. Functions and regulation of RNA editing by ADAR deaminases. Annu Rev Biochem 79: 321–349
- Nishikura K, Yoo C, Kim U, Murray JM, Estes PA, Cash FE, Liebhaber SA 1991. Substrate specificity of the dsRNA unwinding/modifying activity. EMBO J 10: 3523–3532
- Page RD 2002. Visualizing phylogenetic trees using TreeView. Curr Protoc Bioinformatics 6.2.1–6.2.15.
- Paz N, Levanon EY, Amariglio N, Heimberger AB, Ram Z, Constantini S, Barbash ZS, Adamsky K, Safran M, Hirschberg A, et al. 2007. Altered adenosine-to-inosine RNA editing in human cancer. Genome Res 17: 1586–1595
- Peng PL, Zhong X, Tu W, Soundarapandian MM, Molner P, Zhu D, Lau L, Liu S, Liu F, Lu Y 2006. ADARB1-dependent RNA editing of AMPA receptor subunit GluR2 determines vulnerability of neurons in forebrain ischemia. Neuron 49: 719–733
- Peng Z, Cheng Y, Tan BC, Kang L, Tian Z, Zhu Y, Zhang W, Liang Y, Hu X, Tan X, et al. 2012. Comprehensive analysis of RNA-Seq data reveals extensive RNA editing in a human transcriptome. Nat Biotechnol. doi: 10.1038/nbt.2122
- Pickrell JK, Gilad Y, Pritchard JK 2012. Comment on “Widespread RNA and DNA sequence differences in the human transcriptome.” Science 335: 1302 doi: 10.1126/science.1210484
- Polson AG, Bass BL 1994. Preferential selection of adenosines for modification by double-stranded RNA adenosine deaminase. EMBO J 13: 5701–5711
- Prasanth KV, Prasanth SG, Xuan Z, Hearn S, Freier SM, Bennett CF, Zhang MQ, Spector DL 2005. Regulating gene expression through RNA nuclear retention. Cell 123: 249–263
- Schrider DR, Gout JF, Hahn MW 2011. Very few RNA and DNA sequence differences in the human transcriptome. PLoS ONE 6: e25842 doi: 10.1371/journal.pone.0025842
- Seeburg PH, Higuchi M, Sprengel R 1998. RNA editing of brain glutamate receptor channels: Mechanism and physiology. Brain Res Brain Res Rev 26: 217–229
- Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K 2001. dbSNP: The NCBI database of genetic variation. Nucleic Acids Res 29: 308–311
- Taylor DR, Puig M, Darnell ME, Mihalik K, Feinstone SM 2005. New antiviral pathway that mediates hepatitis C virus replicon interferon sensitivity through ADAR. J Virol 79: 6291–6298
- Wagner RW, Smith JE, Cooperman BS, Nishikura K 1989. A double-stranded RNA unwinding activity introduces structural alterations by means of adenosine to inosine conversions in mammalian cells and Xenopus eggs. Proc Natl Acad Sci 86: 2647–2651
- Wang Q, Miyakoda M, Yang W, Khillan J, Stachura D, Weiss M, Nishikura K 2004. Stress-induced apoptosis associated with null mutation of ADAR RNA editing deaminase gene. J Biol Chem 279: 4952–4961
- Wang Q, Zhang Z, Blackwell K, Carmichael GG 2005. Vigilins bind to promiscuously A-to-I-edited RNAs and are involved in the formation of heterochromatin. Curr Biol 15: 384–391
- Wang K, Li M, Hakonarson H 2010. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 38: e164 doi: 10.1093/nar/gkq603
- Zheng H, Fu TB, Lazinski D, Taylor J 1992. Editing on the genomic RNA of human hepatitis delta virus. J Virol 66: 4693–4697
Source: PubMed