Complex regulation of ADAR-mediated RNA-editing across tissues
- Melanie A. Huntley†1Email author,
- Melanie Lou†1,
- Leonard D. Goldstein1,
- Michael Lawrence1,
- Gerrit J.P. Dijkgraaf2,
- Joshua S. Kaminker1 and
- Robert Gentleman1
© Huntley et al. 2015
Received: 4 September 2015
Accepted: 11 December 2015
Published: 15 January 2016
RNA-editing is a tightly regulated, and essential cellular process for a properly functioning brain. Dysfunction of A-to-I RNA editing can have catastrophic effects, particularly in the central nervous system. Thus, understanding how the process of RNA-editing is regulated has important implications for human health. However, at present, very little is known about the regulation of editing across tissues, and individuals.
Here we present an analysis of RNA-editing patterns from 9 different tissues harvested from a single mouse. For comparison, we also analyzed data for 5 of these tissues harvested from 15 additional animals. We find that tissue specificity of editing largely reflects differential expression of substrate transcripts across tissues. We identified a surprising enrichment of editing in intronic regions of brain transcripts, that could account for previously reported higher levels of editing in brain. There exists a small but remarkable amount of editing which is tissue-specific, despite comparable expression levels of the edit site across multiple tissues. Expression levels of editing enzymes and their isoforms can explain some, but not all of this variation.
Together, these data suggest a complex regulation of the RNA-editing process beyond transcript expression levels.
RNA-editing is an enzyme-mediated post-transcriptional modification whereby individual nucleotide residues are altered in the nascent RNA transcript molecule . In contrast to the traditional dogma of DNA makes RNA makes protein, RNA-editing is a non-heritable change to the RNA, which can result in both increased diversity of the transcriptional repertoire, and of the resulting peptide products.
The most common RNA-editing alteration is enzymatic deamination of adenosine nucleotides to inosine, by the family of ADAR enzymes [2, 3]. These are referred to as canonical A-to-I editing events. Since much of the cellular machinery treats inosine as a guanine nucleotide, A-to-I events are typically discovered and referred to as A-to-G in sequencing data.
There are 3 members of the ADAR family in vertebrates: Adar, Adarb1, and Adarb2 (also called Adar1, Adar2, and Adar3, respectively). A number of substrates are reported to be edited by Adar and Adarb1 [4–6]. However, there is no evidence to support the enzymatic activity of Adarb2 .
The gene for Adar produces two main protein isoforms. Expression of the longer Adar isoform, p150, is inducible by interferon stimulation  and has both nuclear and cytoplasmic localization [9, 10]. The shorter isoform of Adar, p110, is more constitutively expressed and nuclear localized . It has been suggested that the nuclear Adar proteins edit pre-mRNAs prior to splicing, while the cytoplasmic p150 isoform may edit viral RNAs or microRNAs in the cytoplasm [4, 11–14].
Owing to the role of interferon in defense against microbial infection, the p150 isoform of Adar has been thought to function in the editing of double stranded viral RNA. In humans, double stranded RNA structures often form in the long 3’ UTR regions of genes when juxtaposed Alu repeat elements are present. These inverted transposable elements, which are remnants of ancient viral insertions, thus form the secondary structure substrate for hyper-editing .
In contrast to the hyper-editing of many adenosine residues in close proximity, site-selective editing tends to occur consistently at specific residues, and even within the coding sequence of transcripts. Aside from the propensity of RNA substrate to form hairpin structures, little is known about other mechanisms regulating this more selective form of editing. The importance of site-selective editing events is highlighted by the necessity of an Adarb1 mediated edit in the neurotransmitter receptor protein GluR-B for developmental viability in mice .
The ADAR gene family is highly conserved across metazoans, and highly expressed in the nervous system [4, 17, 18]. This pattern of evolutionary constraint suggests an important role of RNA-editing in the development and maintenance of the metazoan nervous system. RNA-editing deficiencies are observed in epilepsy, amyotrophic lateral sclerosis (ALS), Aicardi-Goutieres syndrome (AGS), schizophrenia, suicidal depression, and other neurodegenerative diseases [19–21].
RNA-editing has also been implicated in functions and disease outside the nervous system. Site-specific editing of AZIN1 is associated with more aggressive tumor-initiating potential in hepatocellular carcinoma , while an overall decrease in hyper-editing of Alu repeats has been found in brain, prostate, lung, kidney and testis tumors . Analysis of RNA-seq data from The Cancer Genome Atlas (TCGA) demonstrated alterations in editing levels between tumors and matched normal tissues from the same patients [24, 25]. In some cancers, the editing level correlated with patient survival . Cell viability assays for a small number of editing events that result in protein alterations were observed to significantly effect cell survival .
It is clear from these examples that RNA-editing is functionally important, and must require tight regulation to prevent dysfunction. However very little is known about how this regulation occurs. We do not have a thorough understanding of the variability in editing between individuals, nor the intra-individual variability across tissues to properly address this question. Furthermore, there has not yet been a genome-wide study to address, in an unbiased manner, the scope of functionally relevant RNA-editing in healthy individuals.
Previous high throughput RNA-editing surveys have been limited in the range of tissues studied, with a bias towards central nervous system (CNS) tissues, or in their genomic breadth. Transcriptome-wide RNA-seq analyses of RNA-editing in mice have largely been limited to whole brain and cerebral cortex [26, 27].
A valuable study of RNA-editing was performed across seven distinct tissues from a single human individual. However, the study used specific padlock probes, not transcriptome wide RNA-seq, and resulted in only 569 putative RNA-editing sites .
In order to investigate the mechanisms underlying the regulation of RNA-editing, we present an analysis of a novel transcriptome wide RNA-seq dataset from nine tissues in a single individual mouse. Additionally we present a comparison and combined analysis of 20 previously published RNA-seq samples from mice of the same strain across a variety of tissues. To the best of our knowledge, this is the only dataset for which we can address in a genome-wide fashion, questions both on the variability of RNA-editing across tissues in a single individual and within specific tissues across individuals. We then use these data to investigate the underlying mechanisms regulating the important process of RNA-editing.
Collecting high confidence RNA-DNA differences
Genic context of editing sites
We next considered whether certain classes of genes were preferentially edited within our samples as this might reveal a regulatory aspect of the editing process. Within each tissue, we collected genes containing protein altering RDDs and tested for enrichment of any functional category amongst the genes normally expressed within that tissue. We performed a Gene Ontology analysis for each tissue individually, where the universe of genes was simply the set of genes with an RPKM expression level ≥ 10 within that tissue. Additional file 3: Table S1 reports the results of this analysis, filtering for odds ratios > 10, and at least 2 edited genes per category. Only the edited genes within brain and within testis show significant enrichment for any gene ontology. We find that genes involved in ionotropic glutamate receptor signalling are enriched for non-synonymous RDDs amongst all genes expressed within brain tissue. This result is observed in both the high and low coverage brain tissue samples. These sites are highly important for normal brain functioning and often have editing frequencies near 100 %. However, the same six edited genes (Gria2, Gria4, Grik2, Grik5, Fcho2, and Cadps) are found to recur in most of the significant GO categories for the brain samples. Likewise for the testis results, three edited genes (Nup153, Ddx25, and Plekha8) account for the results. Thus there are few ontologies enriched for protein altering RNA-editing, and these are driven by a very small number of edited genes.
Variation in RNA-editing within and across individuals
We next looked at the inter-individual variability of calling RDDs. If RDDs are functionally important, we assume that their presence will be consistent across individuals, when looking in the same type of tissue. Here we included the previously published RNA-seq samples (excluding the male samples from Brawand et al. (2011) due to insufficient sequencing depth), to serve as biological replicates for four of our tissues. We had 7 adult brain samples, and 2 samples each for the heart, kidney, and liver tissues.
Within the 7 brain samples we found a much higher proportion of RDD sites that were shared (by any two individuals, and by all individuals) than across tissues within a single individual (Fig. 4 b). Thus there is more similarity in the editing of brain transcripts from different individuals than from different tissues within the same individual. We also observe an increase in the proportion of brain RDDs shared between individuals when we consider coding and then non-synonymous sites. Thus shared RDDs within brain transcripts, which we presume to be functionally relevant, appear enriched for non-synonymous sites.
We next wanted to see if this pattern held true in the other tissues. Since we only had 2 samples per tissue for the heart, kidney, and liver, for consistency we subset our brain samples to include only 2 samples (one from our generated data, and one from Brawand et al. (2011)). The brain samples had the highest consistency in RDD recurrence within coding and non-synonymous sites (Additional file 4: Figure S4). The other tissues had higher consistency than brain overall, but not for coding and non-synonymous sites. Additionally, except for kidney, there was more similarity in editing within a tissue from different individuals than across tissues. Thus, for the presumed most functionally relevant RDDs, we observe more differences between tissues than between individuals.
Tissue specific RNA-editing
ADAR expression and isoform usage
Interestingly, we did find some examples of RDDs for which the transcript was sufficiently expressed, but where the editing frequency varied across tissues. To examine what might be causing this variability in editing frequency across tissues, we used the RNA-seq data to investigate the expression and isoform usage of the ADAR enzymes. This provided us with the unique opportunity to study the quantitative effect of ADAR expression on RNA-editing.
Adar is most highly expressed in the brain and spleen samples, which is consistent with brain and spleen having higher number of observed RDDs after accounting for sequencing depth (see Fig. 2 a). Adarb1 has the highest expression in brain, and seems to be more brain specific in its expression than Adar. However both Adar and Adarb1 have relatively modest expression in brain, with RPKMs lower than 50 in all samples. Adarb2 is exclusively expressed in the brain, however it is expressed at a very low level, and to date, there is no evidence supporting the editing activity of this enzyme . Thus, Adarb1 and Adarb2 have nearly exclusive expression within brain tissue, while Adar has more appreciable expression across multiple tissues.
From the literature, it is expected that there are differences in Adar isoform usage across tissues, with the shorter p110 isoform more highly expressed in CNS tissues. Therefore, it is important to study gene expression at the isoform level as well. When we examined the data across our samples, we confirmed the elevated expression of the p110 isoform within brain samples (Fig. 6 b). The longer p150 isoform was observed more frequently within the non-brain samples.
While spleen and thymus have Adar expression levels similar to those found in brain, the difference in isoform usage may result in differential regulation of editing between these tissues. There are some reported motifs prefered by Adar and Adarb1, but little is known about the specificity differences between the p110 and p150 isoforms of Adar. Both isoforms contain a z-alpha DNA binding domain at their N-terminal end, but the alternative 5’ UTR and splice of the p150 isoform encodes a second z-alpha domain upstream of the first. This extra z-alpha domain could potentially direct the editing enzymes to prefer a subset of nascent transcripts with proximity to a particular DNA sequence that interact with the z-alpha domains. Thus differences in the resulting protein isoforms may provide a mechanism for different selectivity of editing substrates.
We also discovered a correlation of the p110 (J1 junction) with an internal in-frame alternative splicing event (J9) that results in the addition of 26 amino acids between the second double stranded RNA-binding domain and the deamination domain. This region includes four potential phosphorylation sites. The functional effect of this change, in altering the relative position of the deamination domain to the bound RNA transcript, could play a role in the sequence context of RNA-editing events. Thus isoform usage may influence and correlate with specific motifs surrounding editing sites within a sample.
To test this hypothesis, we analyzed the neighboring nucleotides surrounding edited sites across tissues to look for tissue specific patterns that might elucidate editing specificity of the Adar isoforms (Additional file 6: Figure S6). We did not observe a strong difference in motifs surrounding edit sites across tissues. However, we did find that the genic context of the site appeared to have an effect, with sites in 3’ UTRs having the most consistent neighboring nucleotides across tissues.
Regulation of differential editing
To test the hypothesis that editing enzyme expression levels might account for differences in editing frequency per site across tissues, we first needed a collection of sites that were sufficiently expressed across multiple tissues, and varibly edited. For this we required sites that had a minimum of 10X coverage in at least 2 samples, and for an RDD to have been called at the site within at least 1 sample.
We observed an RDD site in the gene Tmem63b which is expressed in all 9 tissues, but edited only within the brain samples. This site is found within the fifth coding exon upstream of the 3’ UTR, and results in a non-synonymous change from a glutamine to an arginine. The editing frequency at this site within brain samples ranges from 41 to 74 %. We also observed this brain specific pattern of editing in the human Tmem63b gene, using the human Illumina Body Map 2 RNA-seq dataset (NCBI GEO accession GSE30611, Additional file 7: Figure S7).
The gene Flnb contains a previously reported non-synonymous RDD . The authors investigated RNA-editing within mouse brain tissue, and these samples are included in our analysis. We also observe RNA-editing at this site within our brain samples, but additionally observe a wide range of editing frequencies across tissues. In our samples, this site is most highly edited within brain and heart tissue, and this pattern is also found within the human Illumina Body Map 2 dataset (Additional file 7: Figure S7).
We also observed a non-synonymous RDD in the gene Copa. This site was previously unreported in mouse, but has been studied in human. This RDD has a much lower editing frequency, with the highest being 27 %, in one of the heart samples. This variably edited RDD site was one of the very few for which brain did not have the highest editing frequency. Curiously, the pattern of editing frequency across our mouse tissues was not consistent with that observed in the human Illumina Body Map 2 dataset (Additional file 7: Figure S7).
Having identified a collection of variably edited sites, we then tested the hypothesis that editing enzyme expression levels could explain the variations in editing across samples. To do this we calculated RPKM expression levels for the Adar genes (across coding exons), and RPKM estimates for the mutually exclusive 5‘ exons that distinguish the p110 and p150 isoforms of Adar. These expression levels were used to look for correlations with per site editing frequencies across samples.
Previous studies suggest that RNA-editing is prevalent within transcripts of the central nervous system. In fact, the most well known A-to-I editing sites occur within coding regions of neuronal signalling genes. Among these, the transcripts of several genes involved with ionotropic glutamate signalling are typically edited at 100 % frequency. This begs the question of why these sites are not simply changed in the genome rather than edited in the transcript. Several studies report the importance of these editing events in modulating receptor activity in a developmental manner [11, 32–36] and provide a rationale for why these changes are not simply genomically encoded instead. This highlights an advantage of regulating transcriptome diversity via RNA-editing.
Here we have quantified that the number of sites which undergo A-to-I editing is remarkably higher in transcripts derived from the brain. While the ability to detect editing events increases with sequencing depth [22, 37], we continue to observe the brain transcriptome as an outlier even after accounting for this difference. We also confirm previous findings [38, 39] that the average editing frequency per site is elevated in brain tissues. These findings are now consistent across humans, non-human primates, and mice, suggesting a conserved role for RNA-editing in the neural transcriptome.
Interestingly, in our study, the vast majority of edited sites within brain samples occurred in intronic regions of genes. This enrichment for intronic editing accounted for the abnormally high amount of editing sites in brain tissues. The abundance of editing events within intronic sequence is intriguing as we did not expect high levels of sequencing coverage over introns. The RNA-sequencing library prep used in this study was a standard poly-d(T) protocol which should enrich for fully mature, poly-adenlyated mRNA transcripts that have already undergone splicing. However, intron retention appears to be a more common phenomenon within neural transcripts . The recently reported phenomenon of recursive splicing in vertebrates  may result in slower processing of introns within the brain, and this may explain the abundance of intronic sequence we observe within our brain tissue data.
Previous studies have found that RNA editing in mice occurs less frequently than in humans . This is thought to be related to hyper editing of Alu repeat elements within the 3’ UTRs of human transcripts. Alu repeats are an active primate specific SINE retrotransposon and are abundant within the human genome. When Alu sequences are juxtaposed within transcripts, they can form a double-stranded RNA hairpin loop structure which is an ideal substrate for ADAR enzymes. Mice do not have Alu repeats, but they do have four related active SINEs (B1, B2, 1D, and B4). The higher divergence between these mouse repeats is thought to result in lower editing levels than those observed within the Alu elements of humans, as fewer double stranded hairpin loops are expected to form within mouse transcripts. However, we do still observe a high proportion of editing sites within the 3’ UTR of transcripts. In spleen and thymus, the majority of the editing occurs within these exons. Therefore, despite lower levels of editing in mouse samples, the distribution of edited sites is largely consistent between humans and mice.
RNA-editing patterns varied across all tissues examined and sites edited in multiple tissues were not common. We demonstrate that tissue specific editing is largely a consequence of the tissue specific expression of the edited transcript. While our data do not argue against the existence of tissue specific editing cofactors, it appears that a major regulator of editing is the expression of the edited transcripts.
However, we did find evidence for a small number of transcripts expressed in multiple tissues but with variable levels of editing. In several of these cases a change in the expression or isoform usage of the editing enzymes correlates with the observed difference in editing. This suggests that the editing frequency is regulated by expression of a particular editing enzyme. In fact, the editing sites in both Flnb and Copa have been previously studied in the context of human hepatocellular carcinoma . Using both tumor and cell line over-expression and knock down experimental data the authors present a compelling case for Adarb1 levels regulating editing in both genes. While they did not investigate the p150 isoform of Adar, they found that the p110 isoform enhanced editing in Flnb, but not in Copa.
While the gene Tmem63b is expressed in multiple tissues, editing is only observed within brain. This simpler, binary editing status is consistent with the brain having the highest expression of both Adarb1 and nearly exclusive expression of the Adar p110 isoform. However, Adars can dimerize to form homo- and heterodimers [44–49]. This combinatorial effect can add yet another layer of complexity to the regulation of editing, making even the simpler cases of differential editing, such as Tmem63b, less straightforward to explain. After accounting for substrate expression, we believe that the editing enzyme expression and isoform usage plays an important role in regulating the observed differences in the amount and location of editing across tissues and cell types.
This study within mouse tissues is largely consistent with genome wide patterns of editing described in human datasets , and further explains how RNA-editing differences can be differentially regulated. The recurrence in editing patterns within tissues, across individuals, and even species, reinforces RNA-editing as functionally important for organismal development, survival, and evolution. Future studies testing the specificity of Adar isoforms, homo- and hetero-dimer combinations, and enzyme or substrate saturation are needed to further elucidate the regulatory mechanisms of this important cellular process.
Isolation of DNA and RNA from mouse tissues
An 8 week old C57BL/6J male mouse was obtained from The Jackson Laboratory and housed according to the animal use guidelines of Genentech, conforming to California State legal and ethical practices. This animal was euthanized by exsanguination with phosphate buffered saline (PBS) and tissues were collected, flash frozen and stored at –80 °C until further processing. DNA and RNA were extracted with either the All Prep DNA/RNA mini kit (kidney, liver, lung, spleen, testis and thymus), the RNeasy Fibrous tissue mini kit (heart and skin) or the RNeasy Lipid tissue mini kit (brain) according to manufacturers protocol (Qiagen).
We generated 100 bp, paired-end standard RNA-seq libraries from the 9 tissues (brain, skin, heart, liver, thymus, testis, spleen, lung and kidney) of the single adult male mouse described above. Each RNA-seq library resulted in 28 million to 77 million uniquely mapped reads.
Additionally, from this same mouse, we extracted genomic DNA and generated 75 bp paired end whole genome sequence libraries from the brain and spleen tissue. These produced 586 million and 445 million uniquely mapped reads, respectively.
For biological replicates of our dataset from a single mouse, we analyzed previously published RNA-seq data from 3 relevant datasets (GSE30352, GSE39866, and ERP000614), totalling 20 additional samples from 15 animals. The first dataset included 76 bp single-end RNA-seq data from 6 tissues (brain, cerebellum, heart, kidney, liver, and testis) of adult C57BL/6 mice . This experiment used 2 separate mice for each tissue (one female, one male), with the exception of the testis tissue which only had the single male sample. In this dataset, only the female samples were sequenced to moderate depth (11 million to 20 million uniquely aligned reads), while the male samples had low coverage (1 million to 5 million uniquely aligned reads).
The second dataset included 80 bp single-end RNA-seq data from the cerebral cortex of embryonic (E17) and adult (3 - 4 months old) C57BL/6J female mice . This experiment included 4 embryonic samples (resulting in 25 million to 29 million uniquely aligned reads) and 3 adult samples (with 26 million to 29 million uniquely mapped reads).
Finally, the third dataset included whole brain derived 76 bp paired-end RNA-seq data from two adult (ages 2 – 3 months) C57BL/6NJ female mice . These data resulted in 15 million and 18 million uniquely mapped reads.
Read processing and mapping
Our nine RNA-seq samples derived from a single individual mouse, and two whole genome DNA samples were sequenced on an Illumina HiSeq machine. The fastq sequence files for all RNA-seq samples, including the 20 previously published samples, were then filtered for read quality (keeping reads where at least 70 % of the cycles had Phred scores ≥ 23), and ribosomal RNA contamination. The remaining reads were then aligned to the mouse reference genome (mm9) using the GSNAP alignment tool . Alignments were produced using the following GSNAP parameters: “-M 2 -n 10 -B 2 -i 1 -N 1 -w 200000 -E 1 –pairmax-rna =200000 –clip-overlap”. These steps, and the downstream processing of the resulting alignments to obtain read counts and RPKMs per gene (over coding exons of RefSeq gene models) are implemented in the Bioconductor package, HTSeqGenie (v 3.12.0) . Only uniquely mapped reads were used for further analysis.
Our two DNA samples were processed similarly for quality, and aligned with the following GSNAP parameters specific for whole genome sequence: “-M 2 -n 10 -B 2 -i 1 –pairmax-dna =1000 –terminal-threshold =1000 –gmap-mode =none –clip-overlap”.
Calling RNA variants
Single nucleotide variants (SNVs) were first called using the Bioconductor package VariantTools . To increase our confidence in identifying true positive RNA-editing events, we required a mimimum of 10X coverage, at least two unique reads supporting the edited variant and base quality scores of 23 or higher. We discarded variants found solely within the 10 bp at either end of reads, due to the higher proportion of mismatches which tend to occur within these regions.
To infer the strand of a transcript, we considered only RNA-editing events that occurred within annotated genes, and discarded SNVs that were located in more than one gene on opposite strands.
To exclude spurious variants arising from regions of the genome that are problematic for short read alignment, we also discarded SNVs within homopolymer tracts, UCSC repeat masker track elements (low complexity, simple repeat, and satellites), UCSC multi-mappability regions, and self-chain regions (Additional file 8: Figure S1).
For added stringency, we also removed any intronic SNV that was within 5 bp of an exon-intron boundary. Although this was a very small fraction of the total number of SNVs, many were spurious alignments of reads into the intron with only a few basepairs originating from the downstream exon.
Identification of RNA-DNA differences
A collection of high confidence RNA-DNA differences (RDDs) was then produced by stringent filtering of the SNVs from RNA using the whole genome sequence data. After all the filtering above, an SNV in the RNA was considered an RDD if the genomic data showed zero reads supporting the edited base, and had > 20X total coverage at the site in question. Although highly stringent, this gave us confidence in calling the remaining RDD sites RNA-editing events (Additional file 9: Figure S9).
We then further selected RDDs that represented the ADAR mediated A-to-I editing events (A-to-G or T-to-C changes for genes on the positive or negative strand, respectively).
RNA-seq data generated for this study are available from NCBI’s Gene Expression Omnibus (GSE74747). Whole genome sequence data generated for this study are available from the European Nucleotide Archive (ERP010840).
We thank Alison Bruce, Jens Reeder, Thomas Sandmann, Haiyin Chen, Gerard Manning and Noelyn Kljavin for their discussion and contributions supporting this manuscript.
Open Access This 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.
- Gott JM, Emeson RB. Functions and mechanisms of RNA editing. Annu Rev Genet. 2000; 34:499–531.PubMedView ArticleGoogle Scholar
- Polson AG, Crain PF, Pomerantz SC, McCloskey JA, Bass BL. The mechanism of adenosine to inosine conversion by the double-stranded RNA unwinding/modifying activity: a high-performance liquid chromatography-mass spectrometry analysis. Biochemistry. 1991; 30:11507–14.PubMedView ArticleGoogle Scholar
- Savva YA, Rieder LE, Reenan RA. The ADAR protein family. Genome Biol. 2012; 13:252–2012.PubMedPubMed CentralView ArticleGoogle Scholar
- Bass BL. RNA editing by adenosine deaminases that act on RNA. Annu Rev Biochem. 2002; 71:817–46.PubMedPubMed CentralView ArticleGoogle Scholar
- Nishikura K. Functions and regulation of RNA editing by ADAR deaminases. Annu Rev Biochem. 2010; 79:321–49.PubMedPubMed CentralView ArticleGoogle Scholar
- Sommer B, Kohler M, Sprengel R, Seeburg PH. RNA editing in brain controls a determinant of ion flow in glutamate-gated channels. Cell. 1991; 67:11–9.PubMedView ArticleGoogle Scholar
- Chen CX, Cho DS, Wang Q, Lai F, Carter KC, Nishikura KA. third member of the RNA-specific adenosine deaminase gene family, ADAR3, contains both single- and double-stranded RNA binding domains. RNA. 2000; 6:755–67.PubMedPubMed CentralView ArticleGoogle Scholar
- George CX, Samuel CE. Human RNA-specific adenosine deaminase ADAR1 transcripts possess alternative exon 1 structures that initiate from different promoters, one constitutively active and the other interferon inducible. Proc Natl Acad Sci U S A. 1999; 96:4621–6.PubMedPubMed CentralView ArticleGoogle Scholar
- 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:5376–88.PubMedPubMed CentralView ArticleGoogle Scholar
- Poulsen H, Nilsson J, Damgaard CK, Egebjerg J, Kjems J. CRM1 mediates the export of ADAR1 through a nuclear export signal within the Z-DNA binding domain. Mol Cell Biol. 2001; 21:7862–71.PubMedPubMed CentralView ArticleGoogle Scholar
- Seeburg PH, Hartner J. Regulation of ion channel/neurotransmitter receptor function by RNA editing. Curr Opin Neurobiol. 2003; 13:279–83.PubMedView ArticleGoogle Scholar
- Valente L, Nishikura K. ADAR gene family and A-to-I RNA editing: diverse roles in posttranscriptional gene regulation. Prog Nucleic Acid Res Mol Biol. 2005; 79:299–338.PubMedView ArticleGoogle Scholar
- Toth AM, Zhang P, Das S, George CX, Samuel CE. Interferon action and the double-stranded RNA-dependent enzymes ADAR1 adenosine deaminase and PKR protein kinase. Prog Nucleic Acid Res Mol Biol. 2006; 81:369–434.PubMedView ArticleGoogle Scholar
- George CX, Gan Z, Liu Y, Samuel CE. Adenosine deaminases acting on RNA, RNA editing, and interferon action. J Interferon Cytokine Res. 2011; 31:99–117.PubMedPubMed CentralView ArticleGoogle Scholar
- Osenberg S, Dominissini D, Rechavi G, Eisenberg E. Widespread cleavage of A-to-I hyperediting substrates. RNA. 2009; 15:1632–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Higuchi M, Maas S, Single FN, Hartner J, Rozov A, Burnashev N, et al. Point mutation in an AMPA receptor gene rescues lethality in mice deficient in the RNA-editing enzyme ADAR2. Nature. 2000; 406:78–81.PubMedView ArticleGoogle Scholar
- Jin Y, Zhang W, Li Q. Origins and evolution of ADAR-mediated RNA editing. IUBMB Life. 2009; 61:572–8.PubMedView ArticleGoogle Scholar
- Jepson JE, Savva YA, Yokose C, Sugden AU, Sahin A, Reenan RA. 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:8325–37.PubMedPubMed CentralView ArticleGoogle Scholar
- Maas S, Kawahara Y, Tamburro KM, Nishikura K. A-to-I RNA editing and human disease. RNA Biol. 2006; 3:1–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Miyamura Y, Suzuki T, Kono M, Inagaki K, Ito S, Suzuki N, et al. Mutations of the RNA-specific adenosine deaminase gene (DSRAD) are involved in dyschromatosis symmetrica hereditaria. Am J Hum Genet. 2003; 73:693–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Rice GI, Kasher PR, Forte GM, Mannion NM, Greenwood SM, Szynkiewicz M, et al. Mutations in ADAR1 cause Aicardi-Goutieres syndrome associated with a type I interferon signature. Nat Genet. 2012; 44:1243–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen L. Characterization and comparison of human nuclear and cytosolic editomes. Proc Natl Acad Sci U S A. 2013; 110:E2741—7.PubMedGoogle Scholar
- Paz N, Levanon EY, Amariglio N, Heimberger AB, Ram Z, Constantini S, et al. Altered adenosine-to-inosine RNA editing in human cancer. Genome Res. 2007; 17:1586–95.PubMedPubMed CentralView ArticleGoogle Scholar
- Paz-Yaacov N, Bazak L, Buchumenski I, Porath HT, Danan-Gotthold M, Knisbacher BA, et al. Elevated RNA editing activity is a major contributor to transcriptomic diversity in tumors. Cell Rep. 2015; 13:267–76.PubMedView ArticleGoogle Scholar
- Han L, Diao L, Yu S, Xu X, Li J, Zhang R, et al. The genomic landscape and clinical relevance of A-to-I RNA editing in human cancers. Cancer Cell. 2015; 28:515–28.PubMedView ArticleGoogle Scholar
- Danecek P, Nellaker C, McIntyre RE, Buendia-Buendia JE, Bumpstead S, Ponting CP, et al. High levels of RNA-editing site conservation amongst 15 laboratory mouse strains. Genome Biol. 2012; 13:26–2012.PubMedView ArticleGoogle Scholar
- Dillman AA, Hauser DN, Gibbs JR, Nalls MA, McCoy MK, Rudenko IN, et al. mRNA expression, splicing and editing in the embryonic and adult mouse cerebral cortex. Nat Neurosci. 2013; 16:499–506.PubMedPubMed CentralView ArticleGoogle Scholar
- Li JB, Levanon EY, Yoon JK, Aach J, Xie B, Leproust E, et al. Genome-wide identification of human RNA editing sites by parallel DNA capturing and sequencing. Science. 2009; 324:1210–3.PubMedView ArticleGoogle Scholar
- Kiran AM, O’Mahony JJ, Sanjeev K, Baranov PV. Darned in 2013: inclusion of model organisms and linking with Wikipedia. Nucleic Acids Res. 2013; 41:D258—61.PubMedView ArticleGoogle Scholar
- Ramaswami G, Li JB. RADAR: a rigorously annotated database of A-to-I RNA editing. Nucleic Acids Res. 2014; 42:D109—13.PubMedView ArticleGoogle Scholar
- Hoopengardner B, Bhalla T, Staber C, Reenan R. Nervous system targets of RNA editing identified by comparative genomics. Science. 2003; 301:832–6.PubMedView ArticleGoogle Scholar
- Burns CM, Chu H, Rueter SM, Hutchinson LK, Canton H, Sanders-Bush E, et al. Regulation of serotonin-2C receptor G-protein coupling by RNA editing. Nature. 1997; 387:303–8.PubMedView ArticleGoogle Scholar
- Ohlson J, Pedersen JS, Haussler D, Ohman M. Editing modifies the GABA(A) receptor subunit alpha3. RNA. 2007; 13:698–703.PubMedPubMed CentralView ArticleGoogle Scholar
- Rula EY, Lagrange AH, Jacobs MM, Hu N, Macdonald RL, Emeson RB. Developmental modulation of GABA(A) receptor function by RNA editing. J Neurosci. 2008; 28:6196–201.PubMedPubMed CentralView ArticleGoogle Scholar
- Whitney NP, Peng H, Erdmann NB, Tian C, Monaghan DT, Zheng JC. Calcium-permeable AMPA receptors containing Q/R-unedited GluR2 direct human neural progenitor cell differentiation to neurons. FASEB J. 2008; 22:2888–900.PubMedPubMed CentralView ArticleGoogle Scholar
- Wahlstedt H, Daniel C, Enstero M, Ohman M. Large-scale mRNA sequencing determines global regulation of RNA editing during brain development. Genome Res. 2009; 19:978–86.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee JH, Ang JK, Xiao X. Analysis and design of RNA sequencing experiments for identifying RNA editing and other single-nucleotide variants. RNA. 2013; 19:725–32.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Q, Xiao X. Genome sequence-independent identification of RNA editing sites. Nat Methods. 2015; 12:347–50.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen JY, Peng Z, Zhang R, Yang XZ, Tan BC, Fang H, et al. RNA editome in rhesus macaque shaped by purifying selection. PLoS Genet. 2014:;10.Google Scholar
- Braunschweig U, Barbosa-Morais NL, Pan Q, Nachman EN, Alipanahi B, Gonatopoulos-Pournatzis T, et al. Widespread intron retention in mammals functionally tunes transcriptomes. Genome Res. 2014; 24:1774–86.PubMedPubMed CentralView ArticleGoogle Scholar
- Sibley CR, Emmett W, Blazquez L, Faro A, Haberman N, Briese M, et al. Recursive splicing in long vertebrate genes. Nature. 2015; 521:371–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Neeman Y, Levanon EY, Jantsch MF, Eisenberg E. RNA editing level in the mouse is determined by the genomic repeat repertoire. RNA. 2006; 12:1802–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Chan TH, Lin CH, Qi L, Fei J, Li Y, Yong KJ, et al. A disrupted RNA editing balance mediated by ADARs (Adenosine DeAminases that act on RNA) in human hepatocellular carcinoma. Gut. 2014; 63:832–43.PubMedPubMed CentralView ArticleGoogle Scholar
- Jaikaran DC, Collins CH, MacMillan AM. Adenosine to inosine editing by ADAR2 requires formation of a ternary complex on the GluR-B R/G site. J Biol Chem. 2002; 277:37624–9.PubMedView ArticleGoogle Scholar
- Cho DS, Yang W, Lee JT, Shiekhattar R, Murray JM, Nishikura K. Requirement of dimerization for RNA editing activity of adenosine deaminases acting on RNA. J Biol Chem. 2003; 278:17093–102.PubMedView ArticleGoogle Scholar
- Gallo A, Keegan LP, Ring GM, O’Connell MA. An ADAR that edits transcripts encoding ion channel subunits functions as a dimer. EMBO J. 2003; 22:3421–30.PubMedPubMed CentralView ArticleGoogle Scholar
- Chilibeck KA, Wu T, Liang C, Schellenberg MJ, Gesner EM, Lynch JM, et al.FRET analysis of in vivo dimerization by RNA-editing enzymes. J Biol Chem. 2006; 281:16530–5.PubMedView ArticleGoogle Scholar
- Poulsen H, Jorgensen R, Heding A, Nielsen FC, Bonven B, Egebjerg J. Dimerization of ADAR2 is mediated by the double-stranded RNA binding domain. RNA. 2006; 12:1350–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Valente L, Nishikura K. RNA binding-independent dimerization of adenosine deaminases acting on RNA and dominant negative effects of nonfunctional subunits on dimer functions. J Biol Chem. 2007; 282:16054–61.PubMedPubMed CentralView ArticleGoogle Scholar
- Brawand D, Soumillon M, Necsulea A, Julien P, Csardi G, Harrigan P, et al. The evolution of gene expression levels in mammalian organs. Nature. 2011; 478:343–8.PubMedView ArticleGoogle Scholar
- Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010; 26:873–81.PubMedPubMed CentralView ArticleGoogle Scholar
- Pau G, Reeder J. HTSeqGenie: A NGS analysis pipeline. R package. 2013.Google Scholar
- Lawrence M, Degenhardt J, Gentleman R. VariantTools: Tools for Working with Genetic Variants. R package. 2013.Google Scholar