Consistent levels of A-to-I RNA editing across individuals in coding sequences and non-conserved Alu repeats
- Shoshana Greenberger†1, 2,
- Erez Y Levanon†3,
- Nurit Paz-Yaacov4,
- Aviv Barzilai2,
- Michal Safran4,
- Sivan Osenberg4,
- Ninette Amariglio4,
- Gideon Rechavi4 and
- Eli Eisenberg5Email author
© Greenberger et al; licensee BioMed Central Ltd. 2010
Received: 28 April 2010
Accepted: 28 October 2010
Published: 28 October 2010
Adenosine to inosine (A-to-I) RNA-editing is an essential post-transcriptional mechanism that occurs in numerous sites in the human transcriptome, mainly within Alu repeats. It has been shown to have consistent levels of editing across individuals in a few targets in the human brain and altered in several human pathologies. However, the variability across human individuals of editing levels in other tissues has not been studied so far.
Here, we analyzed 32 skin samples, looking at A-to-I editing level in three genes within coding sequences and in the Alu repeats of six different genes. We observed highly consistent editing levels across different individuals as well as across tissues, not only in coding targets but, surprisingly, also in the non evolutionary conserved Alu repeats.
Our findings suggest that A-to-I RNA-editing of Alu elements is a tightly regulated process and, as such, might have been recruited in the course of primate evolution for post-transcriptional regulatory mechanisms.
Site-selective adenosine to inosine (A-to-I) RNA-editing is an essential post-transcriptional mechanism for expanding the proteomic repertoire. It is carried out by members of the double-stranded RNA-specific ADAR family predominantly acting on precursor messenger RNAs . As inosines in mRNA are recognized as guanosines (G) by the ribosome in the course of translation, RNA-editing may lead to a codon exchange resulting in a modified protein. ADAR-mediated RNA editing is essential for the development and normal life of both invertebrates and vertebrates [2–5]. Additionally, altered editing patterns have been found to be associated with various diseases including inflammation , SLE , epilepsy , depression , ALS [10, 11] and malignant brain tumors such as gliomas and astrocytomas [12–14]. A-to-I editing affects numerous sites in the human transcriptome, most of which are located in Alu elements within untranslated regions [15–19]. Alu is a retrotransposon, about 280 bp long, belonging to the class of Short Interspersed Nuclear Elements (SINEs). More than one million copies, are present in the human genome, comprising 10% of the whole genome mass . The precise role of RNA editing in Alu repeats is yet a mystery. However, it might affect gene expression through a number of mechanisms : As inosine pairs with cytosine, editing influence the stability of RNA molecules by creating and disrupting secondary structures. At another level, since inosine is recognized as guanosine by the splicing machinery A-I editing can lead to modification of splice sites in introns, inducing premature termination, frame-shift, or new exon formation [22, 23].
Consistent levels of A-to-I RNA-editing across human individuals were previously observed in a few recoding sites, mainly glutamate receptors [10, 12, 24–26]. However, the variability of the editing level among different human individuals has so far not been studied for neither most recoding editing targets nor the abundant Alu editing. Low variance of the editing level indicates a tight regulation and might attest to the functional importance of the specific editing event.
To test this, we studied the variability in editing levels across different human individuals and a variety of tissues, genes and sites.
Results and Discussion
Editing levels are consistent among skin samples of human individuals both in coding and non-coding sequences
In healthy brain tissue, editing levels for the recoding sites within the glutamate receptor are highly uniform across individuals [10, 12, 24]. As a first step, we tested whether the same is true for other editing sites occurring within coding sequences expressed in other tissues.
First, we investigated 32 skin samples. Although mice ADAR1 or ADAR2 knockout die in-utero or shortly after birth, RNA editing is implicated as relevant to the skin in humans by the observation that mutations in ADAR1 lead to Dyschromatosis symmetrica hereditaria , a pigmenting genodermatosis with an autosomal dominant inheritance reported predominantly in Japanese and Chinese individuals . It is interesting to note that the effect of these mutations on ADAR1 function or editing pattern has not been found.
Next, we tested the variability in editing levels of non-recoding sites. In particular, we focused on editing targets within the primate-specific Alu repeats which constitute the vast majority of editing sites in the human genome. In this category we looked at a single editing site in 3 targets selected from the RNA editing database : FYN (last intron, chr6:112,094,677; (hg18)), CARD11 (the 14th intron, chr7:2,942,082), and KLRD1 (3' UTR, chr12:10,359,728). Using Sequenom analysis of the above 32 skin samples, we found low average editing levels in these sites, 19, 2, and 3% (with standard deviations of 11, 3 and 3) for FYN, KLRD1 and CARD11, respectively (Figure 1B, C). In order to quantify the degree of consistency among individuals, we compared the individual editing levels of the 6 different sites in the FLNA, CYFIP2, BLCAP, KLRD1, CARD11 and FYN genes. A statistically significant difference (p < 0.05 by Mann-Whitney test) between editing levels of the different sites was observed for 13 out of the 15 pairs of editing sites tested (see supplementary Table 1). We therefore conclude that the variation of results across individuals is significantly lower than the difference in editing efficiency between different sites (F-ratio 41.8 by ANOVA analysis; p-value 1.1E-28).
Comparison of editing levels in different tissues
It is established that editing levels are particularly high in some tissues, including the brain . This can be attributed to the elevated level of ADAR2 in the brain. Here we aim at comparing the editing levels in a number of non-brain tissues. For this purpose we analyzed editing data  for clusters of editing sites in the 3'UTR Alu sequences of 3 different genes: MDM4, NRIP3 and THOC5 (see Additional file 1) in the following tissues: lung, kidney, prostate, uterus, liver and glioma tumor (not including normal brain tissue). We found low variability among different tissues and a clear distinction between different editing sites of the 3 genes. That is, given the six values of editing levels in the six non-brain tissues, for two different sites within the same Alu repeat, the tissue-to-tissue variance is low enough to allow one to tell (using a standard statistical test, e.g. t-test or Mann-Whitney) that the two sets of six measurements describe two sites differing in their editing efficiency. The tissue-to-tissue similarity is high enough to allow pair-wise distinction between different sites in a cluster in 51/91, 148/300 and 50/105 of the pairs, in MDM4, NRIP3 and THOC5, respectively (p < 0.05; Mann-Whitney test. supplementary tables 3-5). ANOVA analysis yielded F-ratios 12.71, 18.20 and 5.29 (p-values 7.5E-16, 1.4E-26 and 4.7E-07) for MDM4, NRIP3 and THOC5, respectively, again showing that the tissue-to-tissue variability is an order of magnitude lower than the site-to-site variability within the same cluster. The standard-deviation to mean ratio is 0.34, 0.49 and 0.62 for MDM4, NRIP3 and THOC5 (averaged over all sites in the same gene).
Our findings demonstrate that editing levels display low variability among different human individuals not only in coding but also in the non-coding, non-conserved, Alu sequences.
RNA Editing of several evolutionary-conserved recoding sites is known to be of critical importance to proper cell development and function. Altered editing patterns in these sites are associated with severe phenotypes [4, 8, 35]. It is therefore not surprising to find that editing levels in these sites are uniform among individuals, suggesting tight regulation as one would expect. The three newly discovered recoding sites considered in this work are also extremely well evolutionary conserved, and it is thus reasonable to believe that they also have a critical role, yet to be explored. Consistently, we found that they also all have low variability of editing levels.
Virtually all A-to-I editing events in primates occur in the primate-specific Alu repeats, and are therefore not conserved among many species. The biological importance of these sites, if any, is yet to be shown. We therefore tested the consistency of editing levels in these sites in order to obtain a hint as to their potential role. Editing events of sites with highly variable levels among individuals are unlikely to serve as an important link in a chain of events being part of a biological pathway. The other side of the coin is that sites whose editing exhibit a consistent pattern in terms of editing levels might have been recruited in the course of evolution to serve a functional role. Surprisingly, we found that editing levels in sites with Alu repeats do present us with such consistent patterns. The variability among sites, including neighboring ones, is very high: where some sites show an extremely low level of a few percent, neighboring sites could be edited up to 60-80%. Yet, these seemingly fluctuating patterns are actually consistent among individuals - strong sites are strong in all samples and weak site are weak in all samples.
As all sites are edited by the same ADAR enzymes, what could be the mechanism providing the wide range of efficiencies on one hand, together with significant consistency on the other? We propose that A to I editing is subjected to two levels of control - ADAR expression and structural patterns of the dsRNA. The essential ADARs enzymes are expected to be tightly regulated to have consistent levels among different individuals. Indeed, it was found that editing enzymes are tightly regulated during development [25, 36]. Here too, we found ADAR1 levels to be consistent, in most cases, between different individuals, supporting this assumption (see Additional file 1). Therefore, although the relation between editing levels and ADARs expression is probably not a simple linear one [12, 25, 31], it is likely that ADAR expression controls editing level. However, given a certain level of the ADAR enzymes, it seems likely that sequence and structural differences between sites determine their editing efficiency relative to other sites. The sequence and the resulting dsRNA structure formed by Alu, vary significantly from site to site, but are shared by all samples. Sequence analysis [37, 38]of editing sites revealed a number of weak motifs. However, these alone cannot account for the observed tightly-regulated editing profiles. It is therefore plausible that structural motifs may take a role and should be analyzed as well.
The massive expansion of the Alu repeats in the primate genome has increased by order of magnitudes the amount of A-to-I editing in their transcriptome. The above analysis suggests that this phenomenon provided the primates with thousands of well-controlled and consistent transcriptomic "switches" that can be utilized for biological functions. It is not clear how many of these were actually adopted. One distinct example is the adoption of editing site in the NARF gene to create a new editing-assisted splicing event, resulting in a whole-new primate-specific alternatively-spliced exon . The surprising consistency shown in editing of arbitrarily chosen sites within Alu's tantalizes us to wonder whether this mechanism of increasing diversification by creating new editing-assisted splicing events, could be more widespread. Of particular interest are miRNA target sites within Alu repeats [39, 40], where one can clearly see how regulated control over single nucleotides might result in an efficient mechanism affecting the translation rate of the harboring gene.
In summary, we show that editing events within Alu sequences exhibit a consistent pattern of editing levels across individuals. This might be attributed to sequence and structural motifs controlling the editing efficiency. Therefore, Alu intervention in the genome provides the primates with thousands of well controlled binary transcriptomic switches , available for use as additional regulatory mechanisms. Evidence for sporadic use of these switches already exists, but it is yet unclear how widespread this phenomenon is. As A-to-I editing is most abundant in the brain, the fascinating question than arises whether the above-mentioned mechanism might have played a role in primates' brain evolution.
Human skin tissues
The study was approved by the Institutional Helsinki Committee at Sheba Medical Center, Tel Hashomer, Israel and informed consent was properly obtained by all participants. Human skin tissues were frozen by liquid nitrogen after their removal at surgery and kept at -70°C until further use. Thirty two skin samples were tested; 20 inflammatory skin lesions with the following clinical and pathological diagnoses: Atopic/nummular Dermatitis n = 7, Drug eruption n = 2, Psoriasis n = 1, Allergic contact dermatitis n = 1 and Cutaneus T-cell Lymphoma (Mycosis Fungoides) n = 9. 12 Normal skin samples were collected at the Chaim Sheba Medical Center.
RNA purification, reverse transcription (RT) and A-to-I RNA editing reading
Total RNA was isolated using TRIzol reagent (Invitrogen) according to the manufacturer's instructions. Random-primed cDNA synthesis was done on 2 μg of total RNA using M-MLV reverse transcriptase (Invitrogen) according to the manufacturer's instructions. For analyzing editing levels we used Sequenom (San Diego, CA) MassARRAY Compact Analyzer and MassARRAY Assay Design 2 software, as described before . The primer sequences and reaction conditions are available as supplementary data (Additional file 1).
Direct sequencing and A-to-I editing reading of Alu
In order to compare editing levels of different sites in a highly edited region within an Alu repeat, we direct sequenced the 3' UTR of FYN transcripts. PCR reaction for cDNA products was carried out using 100 ng (1 λ) of cDNA obtained as described above, at an annealing temperature of 60°C. The PCR product was separated by 1.5% agarose gel electrophoresis, and extracted using QiaQuick gel extraction kit (QIAGEN) and according to the manufacturer's instructions. cDNA sequencing was carried out in genetic analyzer 3100 (Applied Biosystems/Hitachi; Foster City, CA), and according to the manufacturer's sequencing protocol. High-quality sequences were analyzed. Sequencing results were read using Sequencher 4.2 software (1991-2004 Gene Codes Corp., Ann Arbor, MI); editing quantification was carried out using DS gene 1.5 software (Accelerys Inc. Discovery Studio 2003; San Diego, CA).
Reference Human Genome version used in the research
Mar. 2006 (hg18) assembly
We thank Ofra Maydan Sherf for performing the sequenome analysis. S.G was supported by Talpiot Medical Leadership Program; Sheba Medical Center, Israel. E.Y.L was supported by the Israel Science Foundation (Legacy Heritage Science Initiative). We thank the Flight Attendant Medical Research Institute (FAMRI) for their support. The work of E.E. was supported by the Israel Science Foundation [grant number 365/06] and the Israel Ministry for Science and Technology (Scientific Infrastructure Program).
- Bass BL: RNA editing by adenosine deaminases that act on RNA. Annu Rev Biochem. 2002, 71: 817-846. 10.1146/annurev.biochem.71.110601.135501.PubMed CentralPubMedView Article
- Chen CX, Cho DS, Wang Q, Lai F, Carter KC, Nishikura K: A third member of the RNA-specific adenosine deaminase gene family, ADAR3, contains both single- and double-stranded RNA binding domains. Rna. 2000, 6 (5): 755-767. 10.1017/S1355838200000170.PubMed CentralPubMedView Article
- Palladino MJ, Keegan LP, O'Connell MA, Reenan RA: A-to-I pre-mRNA editing in Drosophila is primarily involved in adult nervous system function and integrity. Cell. 2000, 102 (4): 437-449. 10.1016/S0092-8674(00)00049-0.PubMedView Article
- Higuchi M, Maas S, Single FN, Hartner J, Rozov A, Burnashev N, Feldmeyer D, Sprengel R, Seeburg PH: Point mutation in an AMPA receptor gene rescues lethality in mice deficient in the RNA-editing enzyme ADAR2. Nature. 2000, 406 (6791): 78-81. 10.1038/35017558.PubMedView Article
- Hartner JC, Walkley CR, Lu J, Orkin SH: ADAR1 is essential for the maintenance of hematopoiesis and suppression of interferon signaling. Nature immunology. 2009, 10 (1): 109-115. 10.1038/ni.1680.PubMed CentralPubMedView Article
- Patterson JB, Samuel CE: Expression and regulation by interferon of a double-stranded-RNA-specific adenosine deaminase from human cells: evidence for two forms of the deaminase. Mol Cell Biol. 1995, 15 (10): 5376-5388.PubMed CentralPubMedView Article
- Laxminarayana D, O'Rourke KS, Maas S, Olorenshaw I: Altered editing in RNA editing adenosine deaminase ADAR2 gene transcripts of systemic lupus erythematosus T lymphocytes. Immunology. 2007, 121 (3): 359-369. 10.1111/j.1365-2567.2007.02582.x.PubMed CentralPubMedView Article
- Brusa R, Zimmermann F, Koh DS, Feldmeyer D, Gass P, Seeburg PH, Sprengel R: Early-onset epilepsy and postnatal lethality associated with an editing-deficient GluR-B allele in mice. Science. 1995, 270 (5242): 1677-1680. 10.1126/science.270.5242.1677.PubMedView Article
- Gurevich I, Tamir H, Arango V, Dwork AJ, Mann JJ, Schmauss C: Altered editing of serotonin 2C receptor pre-mRNA in the prefrontal cortex of depressed suicide victims. Neuron. 2002, 34 (3): 349-356. 10.1016/S0896-6273(02)00660-8.PubMedView Article
- Kawahara Y, Ito K, Sun H, Aizawa H, Kanazawa I, Kwak S: Glutamate receptors: RNA editing and death of motor neurons. Nature. 2004, 427 (6977): 801-10.1038/427801a.PubMedView Article
- Kwak S, Kawahara Y: Deficient RNA editing of GluR2 and neuronal death in amyotropic lateral sclerosis. Journal of molecular medicine (Berlin, Germany). 2005, 83 (2): 110-120.View Article
- Maas S, Patt S, Schrey M, Rich A: Underediting of glutamate receptor GluR-B mRNA in malignant gliomas. Proc Natl Acad Sci USA. 2001, 98 (25): 14687-14692. 10.1073/pnas.251531398.PubMed CentralPubMedView Article
- Paz N, Levanon EY, Amariglio N, Heimberger AB, Ram Z, Constantini S, Barbash ZS, Adamsky K, Safran M, Hirschberg A, et al: Altered adenosine-to-inosine RNA editing in human cancer. Genome research. 2007, 17 (11): 1586-1595. 10.1101/gr.6493107.PubMed CentralPubMedView Article
- Cenci C, Barzotti R, Galeano F, Corbelli S, Rota R, Massimi L, Di Rocco C, O'Connell MA, Gallo A: Down-regulation of RNA editing in pediatric astrocytomas: ADAR2 editing activity inhibits cell migration and proliferation. The Journal of biological chemistry. 2008
- Athanasiadis A, Rich A, Maas S: Widespread A-to-I RNA Editing of Alu-Containing mRNAs in the Human Transcriptome. PLoS Biol. 2004, 2 (12): e391-10.1371/journal.pbio.0020391.PubMed CentralPubMedView Article
- Blow M, Futreal PA, Wooster R, Stratton MR: A survey of RNA editing in human brain. Genome research. 2004, 14 (12): 2379-2387. 10.1101/gr.2951204.PubMed CentralPubMedView Article
- Kim DD, Kim TT, Walsh T, Kobayashi Y, Matise TC, Buyske S, Gabriel A: Widespread RNA editing of embedded alu elements in the human transcriptome. Genome research. 2004, 14 (9): 1719-1725. 10.1101/gr.2855504.PubMed CentralPubMedView Article
- Levanon EY, Eisenberg E, Yelin R, Nemzer S, Hallegger M, Shemesh R, Fligelman ZY, Shoshan A, Pollock SR, Sztybel D, et al: Systematic identification of abundant A-to-I editing sites in the human transcriptome. Nat Biotechnol. 2004, 22 (8): 1001-1005. 10.1038/nbt996.PubMedView Article
- Morse DP, Aruscavage PJ, Bass BL: RNA hairpins in noncoding regions of human brain and Caenorhabditis elegans mRNA are edited by adenosine deaminases that act on RNA. Proc Natl Acad Sci USA. 2002, 99 (12): 7906-7911. 10.1073/pnas.112704299.PubMed CentralPubMedView Article
- Batzer MA, Deininger PL: Alu repeats and human genomic diversity. Nat Rev Genet. 2002, 3 (5): 370-379. 10.1038/nrg798.PubMedView Article
- Chen LL, DeCerbo JN, Carmichael GG: Alu element-mediated gene silencing. Embo J. 2008, 27 (12): 1694-1705. 10.1038/emboj.2008.94.PubMed CentralPubMedView Article
- Rueter SM, Dawson TR, Emeson RB: Regulation of alternative splicing by RNA editing. Nature. 1999, 399 (6731): 75-80. 10.1038/19992.PubMedView Article
- Lev-Maor G, Sorek R, Shomron N, Ast G: The birth of an alternatively spliced exon: 3' splice-site selection in Alu exons. Science. 2003, 300 (5623): 1288-1291. 10.1126/science.1082588.PubMedView Article
- Paschen W, Hedreen JC, Ross CA: RNA editing of the glutamate receptor subunits GluR2 and GluR6 in human brain tissue. Journal of neurochemistry. 1994, 63 (5): 1596-1602. 10.1046/j.1471-4159.1994.63051596.x.PubMedView Article
- Wahlstedt H, Daniel C, Enstero M, Ohman M: Large-scale mRNA sequencing determines global regulation of RNA editing during brain development. Genome research. 2009, 19 (6): 978-986. 10.1101/gr.089409.108.PubMed CentralPubMedView Article
- Keegan LP, Brindle J, Gallo A, Leroy A, Reenan RA, O'Connell MA: Tuning of RNA editing by ADAR is required in Drosophila. Embo J. 2005, 24 (12): 2183-2193. 10.1038/sj.emboj.7600691.PubMed CentralPubMedView Article
- Miyamura Y, Suzuki T, Kono M, Inagaki K, Ito S, Suzuki N, Tomita Y: Mutations of the RNA-specific adenosine deaminase gene (DSRAD) are involved in dyschromatosis symmetrica hereditaria. American journal of human genetics. 2003, 73 (3): 693-699. 10.1086/378209.PubMed CentralPubMedView Article
- Zhang XJ, He PP, Li M, He CD, Yan KL, Cui Y, Yang S, Zhang KY, Gao M, Chen JJ, et al: Seven novel mutations of the ADAR gene in Chinese families and sporadic patients with dyschromatosis symmetrica hereditaria (DSH). Hum Mutat. 2004, 23 (6): 629-630. 10.1002/humu.9246.PubMedView Article
- Clutterbuck DR, Leroy A, O'Connell MA, Semple CA: A bioinformatic screen for novel A-I RNA editing sites reveals recoding editing in BC10. Bioinformatics. 2005, 21 (11): 2590-2595. 10.1093/bioinformatics/bti411.PubMedView Article
- Levanon EY, Hallegger M, Kinar Y, Shemesh R, Djinovic-Carugo K, Rechavi G, Jantsch MF, Eisenberg E: Evolutionarily conserved human targets of adenosine to inosine RNA editing. Nucleic acids research. 2005, 33 (4): 1162-1168. 10.1093/nar/gki239.PubMed CentralPubMedView Article
- Riedmann EM, Schopoff S, Hartner JC, Jantsch MF: Specificity of ADAR-mediated RNA editing in newly identified targets. Rna. 2008, 14 (6): 1110-1118. 10.1261/rna.923308.PubMed CentralPubMedView Article
- Zilberman DE, Safran M, Paz N, Amariglio N, Simon A, Fridman E, Kleinmann N, Ramon J, Rechavi G: Does RNA editing play a role in the development of urinary bladder cancer?. Urologic oncology. 2009
- Paul MS, Bass BL: Inosine exists in mRNA at tissue-specific levels and is most abundant in brain mRNA. Embo J. 1998, 17 (4): 1120-1127. 10.1093/emboj/17.4.1120.PubMed CentralPubMedView Article
- Lev-Maor G, Sorek R, Levanon EY, Paz N, Eisenberg E, Ast G: RNA-editing-mediated exon evolution. Genome Biol. 2007, 8 (2): R29-10.1186/gb-2007-8-2-r29.PubMed CentralPubMedView Article
- Singh M, Kesterson RA, Jacobs MM, Joers JM, Gore JC, Emeson RB: Hyperphagia-mediated obesity in transgenic mice misexpressing the RNA-editing enzyme ADAR2. The Journal of biological chemistry. 2007, 282 (31): 22448-22459. 10.1074/jbc.M700265200.PubMedView Article
- Palladino MJ, Keegan LP, O'Connell MA, Reenan RA: dADAR, a Drosophila double-stranded RNA-specific adenosine deaminase is highly developmentally regulated and is itself a target for RNA editing. Rna. 2000, 6 (7): 1004-1018. 10.1017/S1355838200000248.PubMed CentralPubMedView Article
- Polson AG, Bass BL: Preferential selection of adenosines for modification by double-stranded RNA adenosine deaminase. Embo J. 1994, 13 (23): 5701-5711.PubMed CentralPubMed
- Kleinberger Y, Eisenberg E: Large-scale analysis of structural, sequence and thermodynamic characteristics of A-to-I RNA editing sites in human Alu repeats. BMC Genomics. 2010, 11: 453-10.1186/1471-2164-11-453.PubMed CentralPubMedView Article
- Liang H, Landweber LF: Hypothesis: RNA editing of microRNA target sites in humans?. Rna. 2007
- Smalheiser NR, Torvik VI: Alu elements within human mRNAs are probable microRNA targets. Trends Genet. 2006, 22 (10): 532-536. 10.1016/j.tig.2006.08.007.PubMedView Article
- Paz-Yaacov N, Levanon EY, Nevo E, Kinar Y, Harmelin A, Jacob-Hirsch J, Amariglio N, Eisenberg E, Rechavi G: Adenosine-to-inosine RNA editing shapes transcriptome diversity in primates. Proc Natl Acad Sci USA. 2010, 107 (27): 12174-12179. 10.1073/pnas.1006183107.PubMed CentralPubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.