- Open Access
Comparative analysis of A-to-I editing in human and non-human primate brains reveals conserved patterns and context-dependent regulation of RNA editing
© The Author(s). 2017
- Received: 1 November 2016
- Accepted: 15 March 2017
- Published: 6 April 2017
A-to-I RNA editing is an important process for generating molecular diversity in the brain through modification of transcripts encoding several proteins important for neuronal signaling. We investigated the relationships between the extent of editing at multiple substrate transcripts (5HT2C, MGLUR4, CADPS, GLUR2, GLUR4, and GABRA3) in brain tissue obtained from adult humans and rhesus macaques. Several patterns emerged from these studies revealing conservation of editing across primate species. Additionally, variability in the human population allows us to make novel inferences about the co-regulation of editing at different editing sites and even across different brain regions.
- 5HT2C Receptor
- Editing Site
- Adenosine Deaminase Acting
- Editing Efficiency
- Discrete Brain Region
Modification of messenger RNAs by adenosine-to-inosine (A-to-I) editing is frequently observed in metazoan transcripts and represents a potentially important mechanism for facilitating advantageous molecular diversification in mammals. This kind of RNA editing has been postulated to play a prominent role in nervous system function by re-coding transcripts such that they encode functionally distinct proteins. A-to-I editing involves targeted conversion of specific adenosines on substrate transcripts to inosines by hydrolytic deamination, a reaction carried out by a family of enzymes known as Adenosine Deaminases Acting on RNA (ADARs) which interact with target substrates by binding regions of double stranded RNA. A-to-I editing effects RNA function because inosine base pairs with cytosine with similar Watson-Crick geometry as guanosine. As such, it behaves like guanosine in most cellular contexts including translation, alternative splicing, and RNA induced gene silencing  . Importantly, alterations in protein function as a result of editing have been observed for several neurotransmitter receptors ranging from alterations in intercellular signaling for 5HT2C receptors  to changes in biophysical properties of ion channel subunits GLUR2 , GLUR4 , and GABRA3 . While A-to-I editing is recognized as an important and ubiquitous phenomenon in mammals, little is known about how the process is regulated at individual sites and in distinct tissues. Knockout studies in mice     indicate that the two ADAR enzymes (ADAR1 and ADAR2) have overlapping substrate specificities in certain cases and more distinct specificities in others. Previous work has indicated that many of the sites examined in this study are edited by both Adar1 and Adar2 in mice (5HT2C A site, B site, and C sites; GLUR4-Flip/Flop R/G sites, CADPS, and Gabra3), while the GLUR2 Q/R site and the D site on the 5HT2C receptor appear to be primarily targeted by Adar2. Conversely, the substrate specificity of ADAR enzymes at MGLUR4 While these studies do not clearly delineate ADAR substrate specificity in all relevant contexts they suggest that the unique mechanisms regulate A-to-I editing at each respective substrate. The following analyses were carried out to test the hypothesis that A-to-I editing is regulated by conserved substrate-specific and brain region-specific mechanisms. We focused on several A-to-I editing events in open reading frames which alter amino acid codons in transcripts important for synaptic communication in the CNS. This was accomplished by performing ‘head-to-head’ comparisons of A-to-I editing profiles in brain tissue samples harvested from striatum and cortex of both human and rhesus macaque looking specifically for correlations in the extent of editing across and between substrates as well as with ADAR mRNA levels.
For this study, we simultaneously quantify the extent of editing at several substrates in cortical and striatal tissue samples by targeted multiplex transcript analysis . Interestingly, we observed similar patterns of editing for many of these substrates in human and rhesus macaque brains supporting the notions that the precise extent of editing at each of these substrates is conserved, tightly regulated, and functionally important. These results provide evidence that editing at certain substrates may be co-regulated within certain anatomical contexts and regulated by divergent mechanisms in others. When comparing ADAR mRNA expression levels with editing profiles we observe no direct relationship between mRNA expression and the extent of editing at any site analyzed in this study, supporting the hypothesis that A-to-I editing is regulated downstream of ADAR expression. Interestingly, we did observe a direct relationship between the mRNA expression levels of ADAR1 and ADAR2 mRNA transcripts over a large dynamic range in both species. Furthermore, observations in the human population indicate that certain individuals harbor multi-substrate deficiencies in A-to-I editing that generalize across each of the analyzed substrates while others display apparent deficiencies in editing at only a few substrates in certain neuroanatomical contexts.
Covariant analysis of A-to-I editing
It should be noted that the human tissue samples represent a semi-matched cohort originally obtained to test the hypothesis that RNA processing of 5HT2C transcripts was aberrant in individuals with a genetic disorder known as Prader-Willi Syndrome (PWS). After observing that 5HT2C editing was not altered in the PWS we decided to leverage the remaining available tissue to investigate other aspects of the regulation of RNA editing presented in this study. After we analyzed editing in each tissue sample, we first verified that PWS diagnosis, sex, age, and postmortem interval (PMI) did not have an effect on patterns of editing or RNA expression. Correlation analysis (ANCOVA p > .05) did not reveal any evidence that the parameters discussed here are significantly influenced by diagnosis, sex, age or PMI (Additional file 1: Figure S1).
Overall patterns of editing in the brain
While the mean extent of editing at most sites was found to be similar in monkey and human, the extent of editing at GLUR4-Flip and GLUR4-Flop transcripts is higher in both brain regions in the monkey (Cortex: Flip +25 p < .0001, Flop +11 p < .0001; Striatum: Flip +44 p < .0001, Flop +6 p < .01,). Additional interspecies differences were observed for the extent of editing at various substrates in specific brain regions, these are summarized in Fig. 1. Divergence in the extent of editing at certain substrates while other substrates retain similar patterns across species supports the hypothesis that dynamic regulation of editing occurs in a substrate specific rather than universal manner .
Co-regulation of RNA editing
When profiling the extent of editing in the human cohort we noticed higher levels of variability than we typically observe in rodent cohorts . This was largely driven by a few individuals with editing deficiencies in editing at specific substrates, or in some cases deficiency in editing at all substrates. We took advantage of these deviations to examine the possibility that editing efficiency at certain substrates may be linked with each other, or if deficiencies in editing manifest across multiple brain regions.
We also observed direct correlation in the extent of editing at GLUR4-Flip and GLUR4-Flop transcripts in the cortex (Pearson correlation p > .001 r2 = .62). Thus, the extent of editing at GLUR4-Flip is a strong predictor of the extent of editing at GLUR4-Flop. This direct relationship could be related to the fact that these transcripts are generated by alternative splicing of the same pre-mRNA , and thus are subject to similar regulatory mechanisms. This effect was not observed in striatum further supporting a model where context-specific mechanisms play a role in regulating A-to-I RNA editing in the brain.
Similar regulation of editing across brain regions
GLUR2 also displays a high level of correlation between brain regions but these results are difficult to interpret due to the limited range; as ~98-99% of the transcripts were edited in all tissues from all individuals with the exception of one patient. Interestingly, this individual had the lowest levels of GLUR2 Q/R site editing in both striatum and cortex implying some generalized deficiency in editing at this particular substrate.
Attempts to correlate editing across brain regions and substrates in the rhesus macaque cohort were precluded by the relatively low levels of inter-animal variability. Interestingly, despite the expected genetic diversity of this population, the extent of editing varied little between individuals (Fig. 1). These animals lived in highly controlled environments compared to the real-world experienced by the individuals in the human cohort. This observation supports the notion that some environmental factors may influence the extent of editing, although we cannot rule out the contribution of significant genetic variation within the human cohort.
Our analysis of A-to-I editing in dissected brain regions from humans and non-human primates allowed us to make several observations and revealed some patterns that may be useful for understanding the regulation of A-to-I editing. One of the most interesting observations from the human cohort was the remarkable inter-individual variability. This contrasts with observations made in rodents and primates reared in a laboratory setting which display more similar editing profiles  . We decided to leverage the inter-individual variability in editing within this human cohort to investigate inter-substrate and cross- anatomical regulation of A-to-I editing. This revealed that editing at certain sites might predict the extent of editing at other sites as observed for GABRA3 and 5HT2C. It also revealed that the extent of editing at 5HT2C transcripts generalizes across at least two distinct brain regions analyzed in these studies. Many of the individuals in this cohort expressed ‘normal’ editing profiles for many substrates while displaying deficiencies in editing at only a few sites.The fact that these deficiencies were generally not observed in more than one brain region in most subjects suggests that the deficit does not result from an inherent polymorphism at the substrate transcript effecting secondary structure as such a polymorphism would be expected to manifest across different anatomical contexts. The individuals with low editing at certain substrates in discrete brain regions could represent examples of region specific plasticity in editing efficiency, however we cannot rule out the possibility that editing varies across sub-regions of cortex and striatum such that variations in precise dissection location could contribute to the observed variation. Interestingly, One subject had low rates of editing for all substrates in each brain region suggesting that this individual may harbor some trait locus that results in a global editing deficiency (although we cannot rule out a unique environmental insult). Regardless, the existence of this subject also provides compelling evidence for some upstream mechanism capable of affecting editing globally at all substrates across brain regions. Another subject harbored a deficiency in editing only at 5HT2C transcripts, which generalized across both brain regions. However, editing at each of the other sites was ‘normal’ compared to the rest of the cohort. We tested the hypothesis that a sequence polymorphism within the 5HT2C gene of this individual might affect a component of the secondary RNA structure required for A-to-I editing. However, analysis of the genomic sequence encoding hee predicted 5HT2C transcript RNA duplex did not reveal any polymorphisms. This implies that other factors or possibly polymorphisms at more remote locations could influence editing at this substrate in humans.
It is interesting to consider the prospect that the processes regulating editing are influenced by the environment and can be dynamically tuned throughout life providing adaptive plasticity. Alternatively, the differences could be relatively non-malleable and dictated by a molecular balance coded uniquely in each individual genome. In a further effort to discriminate these models, we looked at the patterns of editing in matched brain regions in 12 rhesus monkeys that where raised in a controlled laboratory environment. The variability in editing was generally much lower at each site in the monkey cohort compared to the humans suggesting that similar environmental conditions may foster the manifestation of similar editing profiles. Despite differences in variability, monkeys and humans displayed similar anatomical patterns of editing at some sites as highlighted in Fig. 1. The observation that these patterns are retained across primate species suggests that regulated expression of specific editing profiles in discrete brain regions confers some utility advantage and has been conserved through primate evolution.
The results of these studies indicate that there are likely global mechanisms responsible for regulating editing at all of these substrates typified by the case of global deficiency in A-to-I editing manifesting in each substrate and across both brain regions analyzed. Moreover, several examples of context specific regulation of editing were observed as some individuals demonstrated deficiencies in editing only at specific substrates and in specific brain regions. Taken together, these results imply that several distinct levels of regulation exist which can affect the efficiency of editing either globally or only in specific contexts. It is tempting to speculate that editing at different substrates is independently regulated by specific mechanisms inherent to each unique transcript. However, we cannot rule out contributions made by independent cell populations; the observed differences could result from altered regulation of editing in specific cells rather than at specific substrate transcripts.
To test the hypothesis that editing is effected by ADAR expression, we measured mRNA expression levels of ADAR1 and ADAR2 in each tissue sample and compared these levels to the extent of editing at each substrate. We observed a range of ADAR mRNA expression levels across tissue samples however none of the relative expression levels appeared to positively correlate with the extent of editing at any substrate analyzed in this study. Of note however, we did observe characteristic ratios of ADAR1:ADAR2 mRNA expression specific to striatum and cortex (Fig. 3). These fixed ratios may have important implications for the regulation of editing through potential competition between ADAR enzymes at certain editing sites.
It has been suggested by some studies that the 5HT2C editing could play a role in the etiology of the human genetic disorder PWS . However, comparison of matched tissue samples from PWS patients and normal controls indicates similar editing patterns regardless of PWS diagnosis. The work presented here does not rule out the possibility that 5HT2C editing could be altered in PWS patients within certain discrete brain regions such as those regulating hunger and satiety. Unfortunately, attempts to procure precisely dissected brain tissue from ventromedial and lateral hypothalamus were unsuccessful due to the rarity of the disease and the limited availability of these relatively small brain structures.
This work clearly demonstrates that mechanisms for maintaining precise, context-specific editing profiles for certain edited transcripts are conserved across primate species. Divergent patterns observed in some human individuals suggest that while editing often falls within characteristic ranges, it can vary widely in certain circumstances. The fact that evolution has conserved editing within specific ranges across species suggests that the differences in editing observed in some human individuals from these studies are likely to have consequences for CNS function. It will be important to identify the neurological manifestations of altered editing profiles by studying larger human cohorts with detailed psychiatric and neurological histories. Identification of factors capable of influencing the editing profile will facilitate more complete understanding of feedback processes that influence the composition of edited transcripts in the brain.
Tissue collection and processing
All human tissue was obtained from the NICHD Brain and Tissue Bank for Developmental Disorders (University of Maryland in Baltimore, MD) some of the samples were initially selected for a different study because they were obtained from individuals with Prader-Willi Syndrome. Covariant analysis indicated that Prader-Willi diagnosis did not significantly influence any of the parameters discussed here. Tissue was processed by crushing into fine powder while frozen using a ‘Cryo-cup’ grinder and pestle (Biospec products, Tulsa, OK) and mechanically homogenized in Tri-reagent (Ambion, Austin, TX). RNA was isolated from the Tri-reagent mixture according to manufacturer’s instructions and stored at−80C until further analysis. Rhesus macaque RNA samples were obtained as a gift from Karoly Mirnics (Vanderbilt University, Nashville, TN) and prepared as described in . All monkeys were female Indian rhesus monkeys (macaca mulatta) from the colony at Oregon National Primate Research Center and were representatives of two groups divided by spontaneous physical activity as indicated in Mitchel et al. . Covariant analysis did not indicate an effect of physical activity on editing at any of the substrates analyzed (Additional file 4: Figure S4).
Tissue quality control
Choroid plexus is an ependymal tissue which expresses high levels of 5HT2C transcript bearing a unique editing profile. For this reason it was important to exclude any samples containing choroid plexus tissue from any 5HT2C editing analysis. This was accomplished by measuring the levels of clotting factor V (F5), a choroid plexus specific marker, in each tissue sample by Real-Time RT-PCR. Any samples containing detectable F5 were excluded from analysis of 5HT2C editing.
Analysis of A-to-I editing
Editing profiles were determined by high-throughput multiplexed transcript (HTMTA) as described in  . Reactions amplifying each respective substrate transcript were performed in parallel using one of 24 unique barcoded primer sets allowing for 24 sample multiplexing.
The FASTX-Toolkit was used to split the reads according each of the 24 different barcodes. For each read, we then used a script for pairwise alignments against all of the reference sequences allowing for only one mismatch (at the editing site) and assigned each match to its corresponding candidate gene. We then counted the nucleotide composition at the editing position and considered any adenosine at the corresponding site “not edited” and any guanosine “edited”, other nucleotides were detected at a very low rate presumably due to sequencing or polymerase errors and were not counted.
Measuring mRNA expression
F5 was detected using a Fam labeled Taqman probe and primer set (Mm00484202_m1) and human ADAR expression was detected using probe and primer sets ADAR1 (Hs00241666_m1) and ADAR2 (Hs00953724_m1), and rhesus ADAR expression was detected with probes for ADAR1 (Hs00241666_m1) and ADAR2 (Rh00955199_m1). For qRT-PCR, first strand cDNA was synthesized using random primers using the High-Capacity Reverse Transcription Kit and expression was analyzed according to manufacturer’s instructions for Taqman real-time PCR (Applied Biosystems, Forster City, CA) on a CFX96 platform (Biorad, Hercules, CA). Parallel amplification efficiencies were verified by PCR miner algorithm  and data was analyzed by CFX manager (Biorad, Hercules, CA).
Analysis of covariance was performed in human samples by plotting the extent of editing at each respective site by Prader-Willi diagnosis, age, sex, and PMI. Covariance was determined by T-test for Prader-Willi diagnosis and sex; or linear regression for the continuous variables ADAR expression, age, and PMI. No significant effect was identified for any of these variables (p > .05). Standard statistical analysis (t-test) was used for pair wise comparisons. When analyzing correlations we used Pearson’s product moment to determine p-values.
Thank you to Dr. Karoly Mirnics for generously providing rhesus RNA samples.
This work was funded by NIH grant P50 MH096972 05 awarded to RE. The funders had no part in the study design, data collection, or analysis.
Availability of data and material
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
RTO Designed the study, carried out molecular analysis, and drafted the manuscript. XW provided informatics support and sequencing analysis. MVM and RBE Assisted in study design. All authors read and approved the final manuscript.
The authors have no competing interests to disclose.
Consent for publication
Human tissue was obtained from the Maryland Brain and Tissue Bank and approved by The University of Maryland Institutional Review Board (HM-HP-00042077). All rhesus macaque tissue was obtained with approval of the Institutional Animal Care and Use Committee of the University of Oregon-Health Sciences.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Li JB, Church GM. Deciphering the functions and regulation of brain-enriched A-to-I RNA editing. Nat Neurosci. 2013;16(11):1518–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Mallela A, Nishikura K. A-to-I editing of protein coding and noncoding RNAs. Crit Rev Biochem Mol Biol. 2012;47(6):493–501.View ArticlePubMedGoogle Scholar
- Burns CM, et al. Regulation of serotonin-2C receptor G-protein coupling by RNA editing. Nature. 1997;387(6630):303–8.View ArticlePubMedGoogle Scholar
- Berg KA, et al. RNA-editing of the 5-HT(2C) receptor alters agonist-receptor-effector coupling specificity. Br J Pharmacol. 2001;134(2):386–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Sommer B, et al. RNA editing in brain controls a determinant of ion flow in glutamate-gated channels. Cell. 1991;67(1):11–9.View ArticlePubMedGoogle Scholar
- Puchalski RB, et al. Selective RNA editing and subunit assembly of native glutamate receptors. Neuron. 1994;13(1):131–47.View ArticlePubMedGoogle Scholar
- Ohlson J, et al. Editing modifies the GABA(A) receptor subunit alpha3. RNA. 2007;13(5):698–703.View ArticlePubMedPubMed CentralGoogle Scholar
- Higuchi M, et al. Point mutation in an AMPA receptor gene rescues lethality in mice deficient in the RNA-editing enzyme ADAR2. Nature. 2000;406(6791):78–81.View ArticlePubMedGoogle Scholar
- Hartner JC, et al. Liver disintegration in the mouse embryo caused by deficiency in the RNA-editing enzyme ADAR1. J Biol Chem. 2004;279(6):4894–902.View ArticlePubMedGoogle Scholar
- Kubota-Sakashita M, et al. A role of ADAR2 and RNA editing of glutamate receptors in mood disorders and schizophrenia. Mol Brain. 2014;7:5.View ArticlePubMedPubMed CentralGoogle Scholar
- Horsch M, et al. Requirement of the RNA-editing enzyme ADAR2 for normal physiology in mice. J Biol Chem. 2011;286(21):18614–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Morabito MV, et al. High-throughput multiplexed transcript analysis yields enhanced resolution of 5-hydroxytryptamine 2C receptor mRNA editing profiles. Mol Pharmacol. 2010;77(6):895–902.View ArticlePubMedPubMed CentralGoogle Scholar
- Jepson JE, et al. Engineered alterations in RNA editing modulate complex behavior in Drosophila: regulatory diversity of adenosine deaminase acting on RNA (ADAR) targets. J Biol Chem. 2011;286(10):8325–37.View ArticlePubMedGoogle Scholar
- Carmel L, Koonin EV, Dracheva S. Dependencies among editing sites in serotonin 2C receptor mRNA. PLoS Comput Biol. 2012;8(9):e1002663.View ArticlePubMedPubMed CentralGoogle Scholar
- Sommer B, et al. Flip and flop: a cell-specific functional switch in glutamate-operated channels of the CNS. Science. 1990;249(4976):1580–5.View ArticlePubMedGoogle Scholar
- Liu Y, et al. Functionally distinct double-stranded RNA-binding domains associated with alternative splice site variants of the interferon-inducible double-stranded RNA-specific adenosine deaminase. J Biol Chem. 1997;272(7):4419–28.View ArticlePubMedGoogle Scholar
- George CX, Samuel CE. Characterization of the 5’-flanking region of the human RNA-specific adenosine deaminase ADAR1 gene and identification of an interferon-inducible ADAR1 promoter. Gene. 1999;229(1-2):203–13.View ArticlePubMedGoogle Scholar
- Mittaz L, et al. Cloning of a human RNA editing deaminase (ADARB1) of glutamate receptors that maps to chromosome 21q22.3. Genomics. 1997;41(2):210–7.View ArticlePubMedGoogle Scholar
- O’Neil RT, Emeson RB. Quantitative analysis of 5HT(2C) receptor RNA editing patterns in psychiatric disorders. Neurobiol Dis. 2012;45(1):8–13.View ArticlePubMedGoogle Scholar
- Doe CM, et al. Loss of the imprinted snoRNA mbii-52 leads to increased 5htr2c pre-RNA editing and altered 5HT2CR-mediated behaviour. Hum Mol Genet. 2009;18(12):2140–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Mitchell AC, et al. Gene transcripts associated with BMI in the motor cortex and caudate nucleus of calorie restricted rhesus monkeys. Genomics. 2012;99(3):144–51.View ArticlePubMedGoogle Scholar
- Hawrylycz MJ, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature. 2012;489(7416):391–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hood JL, et al. Reovirus-mediated induction of ADAR1 (p150) minimally alters RNA editing patterns in discrete brain regions. Mol Cell Neurosci. 2014;61:97–109.View ArticlePubMedGoogle Scholar
- Zhao S, Fernald RD. Comprehensive algorithm for quantitative real-time polymerase chain reaction. J Comput Biol. 2005;12(8):1047–64.View ArticlePubMedPubMed CentralGoogle Scholar