Antisense sequencing improves the accuracy and precision of A-to-I editing measurements using the peak height ratio method
© Rinkevich et al; BioMed Central Ltd. 2012
Received: 22 September 2011
Accepted: 24 January 2012
Published: 24 January 2012
A-to-I RNA editing is found in all phyla of animals and contributes to transcript diversity that may have profound impacts on behavior and physiology. Many transcripts of genes involved in axonal conductance, synaptic transmission and modulation are the targets of A-to-I RNA editing. There are a number of methods to measure the extent of A-to-I RNA editing, but they are generally costly and time consuming. One way to determine the frequency of A-to-I RNA editing is the peak height ratio method, which compares the size of peaks on electropherograms that represent unedited and edited sites.
Sequencing of 4 editing sites of the Dα6 nicotinic acetylcholine receptor subunit with an antisense primer (which uses T/C peaks to measure unedited and edited sites, respectively) showed very accurate and precise measurements of A-to-I RNA editing. The accuracy and precision were excellent for all editing sites, including those edited with high or low frequencies. The frequency of A-to-I RNA editing was comparable to the editing frequency as measured by clone counting from the same sample. Sequencing these same sites with the sense primer (which uses A/G peaks) yielded inaccurate and imprecise measurements.
We have validated and improved the accuracy and precision of the peak height ratio method to measure the frequency of A-to-I RNA editing, and shown that results are primer specific. Thus, the correct sequencing primer must be utilized for the most dependable data. When compared to other methods used to measure the frequency of A-to-I RNA editing, the major benefits of the peak height ratio are that this method is inexpensive, fast, non-labor intensive and easily adaptable to many laboratory and field settings.
A-to-I RNA editing is catalyzed by adenosine deaminases that act on RNA (ADARs) that bind to double stranded pre-mRNAs and convert adenosine (A) to inosine (I) which is recognized by the ribosome as a guanosine . ADARs are found in all animals, but are absent from protists, plants and fungi . There are three ADAR genes in vertebrates. ADAR1 and ADAR2 both have catalytic activity, whereas ADAR3 lacks activity, although the functional domains of ADAR3 are conserved. ADAR3 is likely is a duplicate of ADAR2. A single ADAR gene exists in insects. Drosophila melanogaster dAdar is homologous to vertebrate ADAR2 .
A-to-I RNA editing regulates behavior and life history traits in many phyla of animals. The widespread conservation of this pathway is thought to be a viral defense mechanism [3, 4]. A-to-I RNA editing occurs in protein coding and non-coding sequences, transposable elements, introns, 5' and 3' untranslated regions of the pre-mRNA that may result in changes in the amino acid sequences, splice sites or levels of transcripts [4, 5]. RNA editing frequently results in non-synonymous substitutions that can be critically important for proper function or tissue distribution. For example, RNA editing results in a Q/R amino acid substitution in the pore loop domain of the AMPA receptor GluR-B subunit that makes the channel impermeable to Ca2+. Unedited GluR-B transcripts lead to neuronal death that causes seizures and premature death in editing deficient mice .
The transcripts of many ligand-gated or voltage-sensitive ion channels and G-protein coupled receptors are targets of A-to-I RNA editing [4, 7–9]. Genome wide studies in D. melanogaster have shown wide-spread editing of these genes . RNA editing of these genes is most common in regions that code for functionally important amino acids in the protein. In voltage-gated K+, Na+ and Ca2+ channels, residues involved in channel gating or inactivation are edited. Editing sites on nAChR or GABA receptor subunits occur in crucial areas in the ligand-binding domain and TM2 that forms the channel pore [9–12].
There are a number of dAdar transcripts produced and dAdar expression is important for a number of functions in the fruit fly. Four common transcripts of dAdar result from alternative splicing and these transcripts are themselves subject to RNA editing . D. melanogaster dAdar null mutants show significant deficiencies in motor control and mating that grow progressively worse with age. Nervous system morphology is greatly affected by dAdar null mutants . Behavioral deficits were also seen in D. melanogaster that had reduced dAdar expression. These mutants did not fly, exhibited diurnal activity patterns and displayed temperature sensitive paralysis . Males deficient in dAdar took longer to engage in courtship behaviors and had a dramatically altered courtship song .
There are a variety of methods for measuring the extent of A-to-I RNA editing [17, 18] and these are important for the burgeoning study of A-to-I RNA editing (since 2007, there have been more than 1080 journal articles, book chapters and meeting presentations on A-to-I RNA editing). However, current methods to measure the extent of RNA editing tend to be costly, time consuming and use radioactive materials or generate hazardous waste. Using the Dα6 nicotinic acetylcholine receptor subunit, we demonstrate a detailed and improved method for a highly accurate and precise estimate of the frequency of editing using peak height ratios from Sanger sequencing electropherograms and compare its costs and benefits to other methods for measuring RNA editing frequency.
Sequences of the primers used
Determination of editing with Dα6IR2
Estimates of A-to-I RNA editing frequency between Dα6IR2 and Dα6285F sequencing primers
Expected Editing Rate
0.94 ± 0.02
0.96 ± 0.01
0.92 ± 0.03
0.92 ± 0.03
0.88 ± 0.01*
0.90 ± 0.01
0.87 ± 0.00*
0.88 ± 0.01
0.76 ± 0.01
0.74 ± 0.01
0.75 ± 0.02
0.76 ± 0.01
0.51 ± 0.01*
0.50 ± 0.03
0.50 ± 0.02
0.52 ± 0.01*
0.25 ± 0.03
0.24 ± 0.02
0.23 ± 0.01
0.25 ± 0.02
0.11 ± 0.03
0.09 ± 0.01
0.09 ± 0.02
0.09 ± 0.03
0.07 ± 0.03
0.07 ± 0.02
0.08 ± 0.05
0.07 ± 0.05
0.85 ± 0.12
0.88 ± 0.17
0.97 ± 0.03
0.97 ± 0.02
0.86 ± 0.05
0.86 ± 0.05
0.95 ± 0.01*
0.93 ± 0.00*
0.66 ± 0.04*
0.79 ± 0.04
0.84 ± 0.03*
0.81 ± 0.03*
0.38 ± 0.02*
0.58 ± 0.02*
0.67 ± 0.01*
0.60 ± 0.05*
0.19 ± 0.03*
0.33 ± 0.04*
0.43 ± 0.04*
0.39 ± 0.13
0.06 ± 0.03
0.17 ± 0.02*
0.33 ± 0.06*
0.23 ± 0.10*
0.04 ± 0.02
0.13 ± 0.04
0.25 ± 0.04*
0.23 ± 0.06*
Comparison of the reliability of estimating A-to-I RNA editing of the Dα6 subunit using the peak height ratio method between Dα6IR2 and Dα6285F primers at different rates of expected editing
0.98 ± 0.03
1.00 ± 0.03
0.97 ± 0.05
0.97 ± 0.05
0.94 ± 0.10
0.88 ± 0.10*
0.80 ± 0.07*
0.85 ± 0.04*
Determination of editing with Dα6285F
Determination of Dα6 editing estimates was unreliable (inaccurate and imprecise) when sequenced with Dα6285F at all four editing sites (Figure 1B, Table 2). The observed frequency of editing was significantly different than expected in 17 of 28 comparisons and varied by as much as 23% (Table 2). In all, 20 estimates were significantly different from the expected editing rate by more than 5%. The slope of the line from a plot of observed values of editing vs. those expected was significantly different from 1.0 for three of the four editing sites and the r2 value was less than 1.0 in 3 of the 4 cases (Table 3). This indicates use of the sense Dα6285F primer is not a reliable method for determining the frequency of editing.
Validation of the peak height ratio method using a known sample
We have validated and verified a very accurate, precise, fast and cost-effective method for estimating editing rates. The accuracy and precision of the editing estimates were primer specific (Figures 1 and 2, Tables 2 and 3). Results with primer Dα6IR2 were excellent, while those with primer Dα 6285F were much less reliable (Figure 1). When sequenced with Dα 6IR2, unedited and edited sites appeared as a mix of T/C signals, respectively. In contrast, results with Dα 6285F showed unedited and edited sites as A or G, respectively (Figure 1B, Table 2). Therefore, it is imperative that the proper primer (Dα 6IR2 in this case) be used to accurately estimate the frequency of A-to-I RNA editing.
The disparity in editing estimation between primers may be a consequence of the sequencing chemistry commonly used for dye-terminator cycle sequencing prior to analysis on automated, capillary DNA analyzers. Current nucleotide formulations substitute dITP for dGTP to prevent compression of poly-G tracts in the resulting electropherogram. However, the incorporation of dITP is less efficient than dGTP .
Comparison of methods used to measure A-to-I RNA editing
3 Analyses a
Peak Height Ratio
6 (RNA isolation, RT, PCR, PCR purification, sequencing, analysis)
Poisoned Primer Extension
7 (RNA isolation, RT, primer labeling, PCR, gel, imaging, analysis)
$120.63 F b $123.40 32P
$133.24 F b $157.84 32P
15.0 F b
9 (RNA isolation, RT, PCR, purification, RE digestion, purification, gel, imaging, analysis)
$10.37 F b
$31.10 F b
3.50 F b
Ultra High Throughput Sequencing
9 (RNA isolation, RT, PCR(×2), purification(×2), hybridization, sequencing, analysis)
One to Several c
10 (RNA isolation, RT, PCR, PCR purification, cloning, transformation, colony screening, colony growth, plasmid purification, sequencing, analysis)
Besides these savings of time and money, there are many other technical advantages of using peak height ratios over other methods. RNA editing may introduce or eliminate a restriction enzyme site. Poisoned primer extension utilizes ddGTP as a reaction terminator in edited transcripts. A larger band will result from unedited transcripts because dATP would be incorporated into the product. These methods utilize fluorescent or phosphoimaging systems for quantification.
The major drawbacks to poisoned primer extension and restriction digests are that they only allow for the quantification of a single editing site, they use radiolabeled or fluorescently labeled primers, and are somewhat labor intensive. The high cost of the restriction digest with radioisotopes and poisoned primer extension is mainly due to the high cost of fluorescent dyes and 32P. These two protocols are cost effective only when larger numbers of replications can be performed in a short period of time. In the case of restriction enzyme digestion, the edited site may alter a restriction enzyme recognition sequence for which an inexpensive or widely used restriction enzyme may not be readily available. In the case of adjacent or nearby editing sites, multiple enzymes would be needed to be used to account for recognition site variation. Complications may arise in cases when a restriction enzyme is not available to recognize the change in recognition sequence.
A major limitation of the poisoned primer extension method is that is can only assess one editing site at a time. In order to obtain the same amount of data across four editing sites as used in this experiment, four unique reactions would need to be run on each sample. This would require extensive sample planning and management. Also, using sense primers to assess editing at adjacent editing sites would be particularly troublesome, as in the case of editing sites 415 and 416 of the Dα6 nAChR subunit used in this experiment. Two sense primers would have to be utilized to account for transcripts with edited and unedited versions at site 415 in order to accurately assess editing at site 416 as mismatches at the 3' end of primer with template can lead to reduced amplification efficiency. Conversely, an antisense primer could be used to assess editing at site 416 without regard for editing status at site 415. However, this would require the use of ddCTP as the reaction terminator. The cost to measure the same 4 editing sites of Dα6 by poisoned primer extension would be 15 to 17 fold higher than by the peak height method (Table 4).
Ultra high throughput sequencing (UHTS) is extremely cost prohibitive and better suited for experiments that may not require comparing many biological samples. Next generation UHTS is the most accurate method to measure the frequency of RNA editing. It can even detect rare transcripts that are missed by clone counting methods . The major disadvantages are cost and the short reads generated may only be useful for multiple editing sites that are nearby if the user needs to know what sites are edited on a specific transcript. However, the undeniable major advantage of UHTS is that many editing sites on many transcripts of many genes can be evaluated.
Clone counting can be performed to quantify the frequency of editing. While there are many advantages to this method, the major drawback is that a large number of clones need to be sequenced to ensure an accurate reflection of editing frequency. This process may take a few days to complete and the cost of sequencing a large number of clones may be substantial. Screening colonies for positive inserts usually requires screening many more colonies than will actually go for sequencing. The waste generated by growing colonies on plates and in liquid media needs to properly disposal by autoclaving or incinerating at an approved facility.
The peak height ratio method with the antisense primer utilizes the different intensities in the T/C signal of eletropherograms. The disadvantage of this method is that it does not allow for the identification of which editing sites are edited on each transcript. However, the cost, labor, turnaround time is unequivocal. The peak height ratio method is very advantageous in that there are few steps required to complete an analysis from biological sample to data point.
Applications of this method
The technical efficiency of the peak height ratio method makes it an ideal method for demonstrating the extent of A-to-I RNA editing in high school and introductory biology classes, or in research laboratories as a significant cost reduction method. Additionally, a mobile lab unit could be assembled for rapid sample processing in the field and other remote areas. The RNA extraction, reverse transcription and PCR could be done in the field and processed samples could be mailed to a sequencing facility. Processing field collected animals would help overcome potential changes in allele frequency that may result from genetic drift if they were returned to a lab and reared for a second generation. This would also reduce the risk of sample degradation and shipping nucleic acids across international borders is typically less burdensome than shipping organisms or tissues.
The peak height method could also be used to estimate allele frequencies within populations. Using pools of animals, it is possible to simultaneously evaluate allele frequencies from as many as 10 diploid or 20 haploid individuals based on the upper (0.95) and lower (0.05) detection rates that we used in this analysis. This would allow for the detection of a single allele out of 20 potential alleles. It is likely that this method could be validated for even lower detection rates. This method would be extremely valuable in our lab for evaluating the frequency of insecticide resistance alleles from field collected populations.
The accuracy and precision of the estimate of A-to-I RNA editing using the peak height ratio method with sequences of Dα6 from the antisense Dα 6IR2 primer is in very good agreement with expected values and is comparable to the quantitative clone counting method. It is also very cost effective and fast compared to other current methods, especially when evaluating editing at multiple sites. Because of these many advantages, it is likely that this method will prove to be a powerful and useful tool in evaluating the extent of A-to-I RNA editing in this rapidly growing field of study, and has other uses in other fields (such as population genetics) as well.
This study was funded by Dow Agrosciences (FDR and JGS) with no competing interests.
- Nishikura K: Functions and regulation of RNA editing by ADAR deaminases. Annu Rev Biochem. 2010, 79: 321-349. 10.1146/annurev-biochem-060208-105251.PubMedPubMed CentralView ArticleGoogle Scholar
- Jin Y, Zhang W, Li Q: Origins and evolution of ADAR-mediated RNA editing. IUBMB Life. 2009, 61 (6): 572-578. 10.1002/iub.207.PubMedView ArticleGoogle Scholar
- Bass BL, Weintraub H: An unwinding activity that covalently modifies its double-stranded RNA substrate. Cell. 1988, 55 (6): 1089-1098. 10.1016/0092-8674(88)90253-X.PubMedView ArticleGoogle Scholar
- Bass BL: RNA editing by adenosine deaminases that act on RNA. Ann Rev Biochem. 2002, 71: 817-846. 10.1146/annurev.biochem.71.110601.135501.PubMedPubMed CentralView ArticleGoogle Scholar
- Fukui T, Itoh M: RNA editing in P transposable element read-through transcripts in Drosophila melanogaster. Genetica. 2010, 138: 1119-1126. 10.1007/s10709-010-9499-z.PubMedView ArticleGoogle Scholar
- Higuchi M, Maas S, Single FN, Hartner J, Rozov A, Burnashev N, Feldmeyer D, Sprengel R, Seeburg PH: Point mutation in an AMPA recepetor gene rescues lethality in mice deficient in the RNA-editing enzyme ADAR2. Nature. 2000, 406: 78-81. 10.1038/35017558.PubMedView ArticleGoogle Scholar
- O'Connell MA: RNA editing: rewriting receptors. Curr Biol. 1997, 7: R437-R439. 10.1016/S0960-9822(06)00212-0.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Hoopengardner B, Bhalla T, Staber D, Reenan R: Nervous system targets of RNA editing identified by comparative genomics. Science. 2003, 301: 832-836. 10.1126/science.1086763.PubMedView ArticleGoogle Scholar
- Jones AK, Buckingham SD, Papadaki M, Yokota M, Sattelle BM, Matsuda K, Sattelle DB: Splice-variant and stage-specific RNA editing of the Drosophila GABA receptor modulates agonist potency. J Neurosci. 2009, 29: 4287-4292. 10.1523/JNEUROSCI.5251-08.2009.PubMedView ArticleGoogle Scholar
- Yao X, Song F, Zhang Y, Shao Y, Li J, Liu Z: Nicotinic acetylcholine receptor beta-1 subunit from the brown planthopper, Nilaparvata lugens: A-to-I RNA editing and its possible roles in neonicotinoid sensitivity. Insect Biochem Mol Biol. 2009, 39: 348-354. 10.1016/j.ibmb.2009.02.001.PubMedView ArticleGoogle Scholar
- Song W, Liu Z, Tan J, Nomura Y, Dong K: RNA editing generates tissue-specific sodium channels with distinct gating properties. J Biol Chem. 2004, 279 (31): 32554-32561. 10.1074/jbc.M402392200.PubMedPubMed CentralView ArticleGoogle Scholar
- 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: 1004-1018. 10.1017/S1355838200000248.PubMedPubMed CentralView ArticleGoogle Scholar
- 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: 437-449. 10.1016/S0092-8674(00)00049-0.PubMedView ArticleGoogle Scholar
- Jepson JEC, Reenan RA: Adenosine-to-inosine genetic recoding is required in the adult stage nervous system for coordinated behavior in Drosophila. J Biol Chem. 2009, 284 (45): 31391-31400. 10.1074/jbc.M109.035048.PubMedPubMed CentralView ArticleGoogle Scholar
- Jepson J, Savya Y, Yokose C, Sugden A, Sahin A, Reenan R: Engineered alterations in RNA editing modulate complex behavior in Drosophila: Regulatory diversity of adenosine deaminase acting on RNA (ADAR) targets. J Biol Chem. 2010, 286 (10): 8325-8337.PubMedPubMed CentralView ArticleGoogle Scholar
- Jepson JEC, Reenan RA: Genetic approaches to studying adenosine-to-inosine RNA editing. Meth Enzymol. 2007, 424: 265-287.PubMedView ArticleGoogle Scholar
- Abbas AI, Urban DJ, Jensen NH, Farrell MS, Kroeze WK, Mieczkowski P, Wang Z, Roth BL: Assessing serotonin receptor mRNA editing frequency by a novel ultra high-throughput sequencing method. Nucl Acids Res. 2010, 38 (10): e118-10.1093/nar/gkq107.PubMedPubMed CentralView ArticleGoogle Scholar
- Innis MA, Myambo KB, Gelfand DH, Brow MA: DNA sequencing with Thermus aquaticus DNA polymerase and direct sequencing of polymerase chain reaction-amplified DNA. PNAS. 1988, 85: 9436-9440. 10.1073/pnas.85.24.9436.PubMedPubMed CentralView ArticleGoogle Scholar
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.