Dynamic response of RNA editing to temperature in Drosophila
© Reenan et al.; licensee BioMed Central. 2015
Received: 15 December 2014
Accepted: 15 December 2014
Published: 3 January 2015
Adenosine-to-inosine RNA editing is a highly conserved process that post-transcriptionally modifies mRNA, generating proteomic diversity, particularly within the nervous system of metazoans. Transcripts encoding proteins involved in neurotransmission predominate as targets of such modifications. Previous reports suggest that RNA editing is responsive to environmental inputs in the form of temperature alterations. However, the molecular determinants underlying temperature-dependent RNA editing responses are not well understood.
Using the poikilotherm Drosophila, we show that acute temperature alterations within a normal physiological range result in substantial changes in RNA editing levels. Our examination of particular sites reveals diversity in the patterns with which editing responds to temperature, and these patterns are conserved across five species of Drosophilidae representing over 10 million years of divergence. In addition, we show that expression of the editing enzyme, ADAR (adenosine deaminase acting on RNA), is dramatically decreased at elevated temperatures, partially, but not fully, explaining some target responses to temperature. Interestingly, this reduction in editing enzyme levels at elevated temperature is only partially reversed by a return to lower temperatures. Lastly, we show that engineered structural variants of the most temperature-sensitive editing site, in a sodium channel transcript, perturb thermal responsiveness in RNA editing profile for a particular RNA structure.
Our results suggest that the RNA editing process responds to temperature alterations via two distinct molecular mechanisms: through intrinsic thermo-sensitivity of the RNA structures that direct editing, and due to temperature sensitive expression or stability of the RNA editing enzyme. Environmental cues, in this case temperature, rapidly reprogram the Drosophila transcriptome through RNA editing, presumably resulting in altered proteomic ratios of edited and unedited proteins.
KeywordsADAR RNA editing RNA thermometers Temperature Editosome Poikilotherm Inosine A-to-I editing
Natural DNAs are usually limited to double-stranded helical shapes; however, RNA is different — the repertoire of possible RNA secondary and tertiary structures appears limitless. RNA secondary structure is strongly correlated with function, and both the structure and corresponding thermodynamic stability of an RNA molecule contribute to functional regulation . Dynamic RNA structures are acutely responsive to input in the form of molecular and environmental factors; it is this mutability of RNA structure that allows RNA to act as a sensor and elicit rapid cellular responses .
RNA structure is fundamentally sensitive to abiotic factors, such as temperature and metal ion concentration. Bacterial RNA thermometers, riboswitches sensitive to temperature, are responsive regulatory elements that control translation of heat-shock, cold-shock, and virulence genes [3,4]. Yet, there is no direct evidence of translational RNA thermometers in eukaryotes. With the addition of large expanses of intronic sequence, eukaryotic RNA thermometers could be considerably less conserved than those found in bacteria, confounding detection. Indeed, the regulation of alternative splicing by the eukaryotic thymidine pyrophosphate riboswitch depends on complex long-distance nucleotide interactions . Therefore, temperature-sensitive structures found in eukaryotic mRNA could, in theory, act anywhere in the transcript to alter processing, translation, transport, degradation or protein binding.
Adenosine-to-inosine (A-to-I) RNA editing is a post-transcriptional modification known to be directed by secondary  and tertiary RNA structures, including those that act over a distance of up to several thousand nucleotides [7,8]. We reasoned that there might be a class of eukaryotic RNA thermometer-like structures that, instead of controlling translation, rather exert their effects on RNA editing levels. A-to-I RNA editing involves the hydrolytic deamination of adenosine into inosine, which is read as guanosine by the protein synthesis machinery. Editing, therefore, has the ability to recode the transcriptome at select sites  and diversify the proteome.
ADARs (adenosine deaminases acting on RNA), the highly conserved proteins responsible for A-to-I editing in all metazoa, are found localized to both cytoplasm and nucleus. While there are multiple ADAR proteins in mammals, the mammalian ADAR2 and single Drosophila ADAR (dADAR) appear to function primarily in the neuronal nucleus . This is consistent with the observation that ADAR’s target transcripts encode proteins involved in chemical and electrical neurotransmission. Phenotypes in model organisms with editing deficiencies range from embryonic lethality  and seizures  (mouse) to defects in motor control  (Drosophila) and chemotaxis  (Caenorhabditis elegans). Such phenotypic consequences of ADAR deficiency indicate that editing plays an integral role in organismal behavior and viability.
Some evidence suggests that environmental factors, specifically temperature, affect editing of select transcripts [15-17], and that RNA structure may regulate this relationship [18,19]. For example, Garrett and Rosenthal showed that editing in an octopus delayed rectifier potassium channel transcript correlates with ambient water temperature for different species, and even suggest that editing differences contribute to octopus adaptation . Additionally, auto-editing of the Drosophila adar (dadar) transcript decreases as temperature increases. Given that the ratio of edited to unedited dADAR isoforms fine-tunes the global ‘editosome’  as well as complex organismal behavior, these data point toward an intriguing role of thermal control in global RNA editing .
However, little research has delved into the widespread effect of temperature on editing, especially in a genetically tractable organism. Because ADARs recognize duplex RNA structures, the information directing editing is carried in RNA secondary and tertiary structures, and changes in editing should partially reflect thermal alteration in RNA shapes . To observe temperature effects on RNA editing we used Drosophila, a poikilotherm with robust RNA editing machinery and a well-characterized editosome.
Although poikilotherms experience environmental temperature undiminished, they have evolved physiological processes that allow them to thrive within the natural temperature ranges of their habitats . It is possible that differential RNA editing is one process that allows poikilothermic animals such as Drosophila to rapidly adapt to varying environmental temperatures. Approximately 28% of ADAR-mediated RNA modifications occur in coding regions in the Drosophila transcriptome . More specifically, A-to-I RNA editing in coding regions leads mostly to non-synonymous amino acid substitutions within proteins involved in neurotransmission [24,25]. As nervous system function is highly sensitive to temperature fluctuations , this suggests that editing could provide a mechanism of rapid temperature responsiveness within the nervous system.
The above observations and hypotheses led us to undertake a comprehensive survey of 54 Drosophila editing sites across a 20°C biologically relevant temperature range . We created a high throughput batch-processing method of editing analysis to test the thermo-sensitivity of editing, which revealed diverse patterns of temperature response between different editing sites. Part of the overall pattern seen was due to altered dADAR levels and largely consistent with observations from parallel heat-shock studies . Certain alterations in editing due to elevated temperature are reversible, although individual editing sites respond differently. To test the contribution of RNA structure to temperature-sensitive editing changes, we used genetically engineered Drosophila strains in which stabilizing and destabilizing structural mutations have been introduced into the endogenous locus of a natural editing substrate , revealing altered temperature sensitivity of editing due to engineered RNA structural changes. Our results suggest that RNA editing is acutely sensitive to temperature, and that this response is partially affected by the thermo-sensitive secondary and tertiary RNA structures that direct editing.
Widespread RNA editing decreases at elevated temperature
We grew engineered D. melanogaster at 25°C before transferring newly-ecclosed male animals to equivalent incubators held at 10°C, 20°C and 30°C. We used only males to control for possible sex-specific differences. Drosophila eggs are most viable at 20°C ; therefore, the 20°C temperature range chosen is biologically relevant. After 72 hours of temperature acclimation, the animals were snap frozen at −80°C. We then performed RNA editing analysis using previously published methods .
Thermo-sensitive editing at individual sites occupies a range of possible patterns. For example, editing at sites 2 to 4 in the calcium sensor synaptotagmin-1 transcript is universally insensitive to temperature (Figure 1B). However, editing of paralytic at sites 1 to 3 decreases significantly at 30°C (Figure 1C). Contrary to the overall trend (Figure S2A), editing of site 6 in the potassium channel shab transcript is potentiated by temperature, significantly increasing in editing level over the temperature range (Figure 1D). Notably, shab site 7, located adjacent to site 6, is edited at nearly 100%, and is temperature-insensitive within the 20°C range studied (Figure 1D).
As there are several examples of multiple editing sites within a single transcript that respond differently to temperature (see Additional file 4: Figure S3), it is unlikely that the observed changes in editing are due to fluctuating transcript levels. Thus, while the general trend is for many sites to decrease editing level at 30°C, there are notable examples that are temperature invariant, as well as those whose temperature profile is counter to the overall trend.
Thermo-sensitive editing patterns are largely conserved
The quite diverse thermo-responsive editing patterns at paralytic sites 1 to 3 (Figure 2B), shab site 6 (Figure 2C), sites 1 to 3 in the complexin SNARE-binding protein transcript (see Additional file 5: Figure S4A), and dadar auto-editing (see Additional file 5: Figure S4B) are generally highly conserved, although absolute editing varies slightly and slope also varies significantly between species (see Additional file 6: Table S2). However, editing at the single site in uncoordinated-13, which encodes a protein involved in calmodulin binding in the larval central nervous system, is not conserved, either in absolute levels or in response to temperature (see Additional file 5: Figure S4C, Additional file 6: Table S2).
It is possible that the slight but significant changes in thermo-responsiveness of editing between species are a result of altered RNA structures, perhaps due to single nucleotide changes or polymorphisms within the genes of the Drosophila species studied here. As most structures that direct editing are formed between exonic and intronic sequences, changes in primary sequence often occur in intronic cis elements, leading to alterations in RNA structure and a corresponding change in editing. However, it is notable that editing at both dadar  (see Additional file 5: Figure S4B) and unc-13 (RA Reenan, personal observation, Additional file 5: Figure S4C) is directed by an entirely exonic secondary structure, whereas synaptotagmin-1  and paralytic  are sites directed by structures comprised of paired exon and intron sequences. One would, therefore, expect the structures of dadar and unc-13 to be under higher sequence—and, therefore, structural—conservation, leading to highly similar editing at these sites across Drosophilidae. While auto-editing of the dadar transcript satisfies this prediction, editing in unc-13 is highly variable (see Additional file 5: Figure S4B and C), suggesting that either the unc-13 RNA structure or the specificity of editing machinery has evolved in the species studied, leading to the observed variability.
Our observations suggest that, in general, patterns of editing responsiveness to temperature are highly conserved and gene-specific in nature. The fact that insensitivity to temperature, temperature-sensitivity and temperature-potentiated site patterns are conserved is significant, and further suggests that while RNA structures that direct editing at certain sites are highly thermodynamically stable, others may be inherently more temperature labile, allowing editing levels to track with temperature based upon an as yet unknown interaction between temperature and RNA structure.
ADAR protein expression level is temperature-sensitive
Changes in editing levels may manifest due to effects at the RNA level (RNA folding, interactions with transcription/splicing machinery, negative influence of RNA chaperones) and/or directly via changes in dADAR protein concentration. We assayed dADAR protein in extracts from the heads of male flies reared at the experimental temperatures, as described above. We used engineered animals in which Jepson et al. tagged the endogenous dadar gene with an HA (hemagglutinin) epitope sequence using homologous recombination (HR) . The appropriate control for mutations generated via HR is an allele in which no targeted mutations have been introduced through the engineering process, but which contains a loxP sequence, the only remnant from the HR procedure, within an intron . The intermediate allele generated during the HR process contains a 5.5 kb mini-white gene, inserted within the same intron in the dadar locus, which produces a viable dADAR hypomorph with approximately 20% of the protein found in the loxP ADAR+ control .
Temperature-dependent editing patterns in loxP control and dadar hypomorphic animals are different and diverse even within a transcript (see Additional file 8: Figure S6), suggesting that dADAR levels are not fully responsible for temperature responsiveness. Finally, when editing in the hypomorph is compared to editing at different temperatures (see Additional file 9: Figure S7), it is clear that the dADAR level accounts for some, but not all, of the observed temperature dependent editing patterns (R = 0.48). Previously, we defined editing sites as ‘low-efficiency’ or ‘high-efficiency’, depending on their responsiveness to the hypomorphic dadar allele. Low-efficiency sites are sensitive to decreased dADAR concentrations, while high-efficiency sites are largely insensitive and are edited at a fairly constant level regardless of dADAR concentration . In the present study, we manipulate the environment (30°C) to mimic the genetic hypomorph. However, the pattern of low-efficiency and high-efficiency sites is different between environmental (30°C) and genetic (hypomorph) lowered dADAR states (Figure 3 and Additional file 9: Figure S7), suggesting that the dADAR concentration is not the only factor driving editing levels and another temperature-sensitive factor, such as RNA folding, also directs editing.
Thermo-sensitive editing response is partially reversible
Edited dADAR isoforms are not differentially sensitive to temperature
The dADAR protein auto-edits its own transcript, causing an amino acid substitution at position 458 in the deaminase domain. Unedited dadar transcripts produce a peptide encoding a serine at position 458 (dADARS), while edited transcripts produce the glycine isoform (dADARG). We previously reported that the dADARS and dADARG isoforms target the same editing sites, but edit them to slightly different levels, with the dADARG isoform displaying less activity in vivo . Auto-editing is stable between 10°C and 20°C at about 55%, although dadar editing decreases to 37% at 30°C (see Additional file 2: Figure S2B, Additional file 3: Table S1), increasing the proportion of the unedited dADARS isoform. We, therefore, suspected that our observed overall temperature-dependent editing patterns (Figure 1A) might be a result of differential stability of the dADARS and dADARG isoforms.
Structural RNA mutations affect temperature-dependent patterns
If RNA structures directing editing are partially responsible for the thermo-sensitivity of editing at specific sites, we reasoned that as temperature changes, the specific structures required for RNA editing may be more or less common leading to altered editing levels.
The ‘DCS delete’ mutation, which is predicted to decrease secondary structure (Figure 6B), decreases editing at all three paralytic editing sites, while the ‘DCS zip’ mutation (Figure 6C), which increases secondary structure, is known to increase editing at all three sites . The thermo-responsive editing profile of these structural mutations appears to be a product of both temperature and structure. The ‘DCS delete’ mutation displays decreased editing at all three sites compared to the loxP control, as expected, yet the thermo-responsive editing pattern, measured as slope, is significantly different at all three sites than that seen in the loxP control (Figure 6D, Additional file 12: Table S3). Similarly, the ‘DCS zip’ mutation shows decreased thermo-responsive editing patterns at all three sites compared to the loxP control (Figure 6D), in addition to increasing editing in all conditions, as expected . The altered thermo-sensitivity of these mutations, defined here as slope, is very subtle, yet is still significant (see Additional file 12: Table S3).
We also tested the thermo-responsiveness of structural mutations in the paralytic tertiary pseudoknot structure (see Additional file 12: Table S3, Additional file 13: Figure S10 A-D). The individual pseudoknot mutations, ‘Loop > α’ and ‘Dock > α’, both disrupt the tertiary RNA interaction, and selectively abolish editing at site 1, while preserving editing at sites 2 and 3. The double mutation, ‘Loop/Dock > α/α’, combines these mutations and is predicted to restore the tertiary pseudoknot structure, rescues editing at site 1 and increases editing at all three sites . These mutations reveal temperature sensitivity patterns (slope) distinct from that seen in the loxP control (see Additional file 12: Table S3, Additional file 13: Figure S10E), suggesting that engineered RNA structures are affected differently by temperature, which manifests in the thermo-sensitive editing response of paralytic editing sites.
The effects of acute environmental changes on the process of RNA editing have not previously been assessed in great detail. Here, we demonstrate the effect of temperature on RNA editing patterns using multiple known target RNAs. Our analysis reveals that editing generally decreases at 30°C, but that there are exceptions to this trend and that the temperature profile of each editing site is distinct (Figure 1A). As D. melanogaster exhibit fairly constant viability and fertility between approximately 12°C and 30°C , we suggest that the changes we observe are normal and represent an adaptive cellular response rather than a breakdown in editing efficiency. Indeed, the potentiation of certain editing sites by elevated temperature, in opposition to the overall trend, argues for a complex underlying hierarchy of responses to temperature (Figure 1D, Additional file 2: Figure S2B).
Congruent with these data is our previous observation that certain sites are more sensitive to dADAR auto-editing and concentration of dADAR protein . In addition, all known editing sites have varying degrees of duplex RNA structure, although other than double-strandedness, no common sequence or structural property has been identified. This suggests that editing sites may behave differently in response to temperature because variable RNA structures act as mechanistically distinct ‘thermometers’, depending on the thermostability of the final RNA structure that directs editing at that particular adenosine, or the thermal influence on the folding pathway of that structure.
This hypothesis is supported by the observation that neighboring sites within a transcript, which are predicted to be involved in the same governing RNA structure, often show similar temperature profiles (Figure 1B-C). However, there are other transcripts (Figure 1D, Additional file 4: Figure S3) in which adjacent sites behave very differently in response to temperature. This could indicate that these sites are ‘more functionally important’ than labile sites. Indeed, with a few noted exceptions, the top 25% of sites ranked by editing level in Figure 1A are largely temperature insensitive.
We suggest that sites that are edited close to 100% are edited at this level regardless of temperature, perhaps because these structures are the most thermodynamically stable, and possible alternative structures are much less likely to form within the investigated temperature range. In support, Tian et al. observed that while editing of a wild type Gabra3 I/M site is insensitive to temperature, synthetic constructs containing structural mutations are temperature-sensitive. This work suggests that thermodynamically less stable structures carry an intrinsic temperature sensitivity, adding that there must be an evolutionary advantage to temperature-insensitivity at this particular site .
The five Drosophilidae species studied here encompass over 10 million years of divergent evolution (Figure 2D) and represent diverse geographical (and temperature) ranges. While D. melanogaster and D. simulans are now cosmopolitan species, D. yakuba is found in savannah climes , D. erecta is native to west central Africa , and D. sechellia was previously confined to the Seychelles islands . Our data show that, generally, editing versus temperature profiles is conserved across species (Figure 2, Additional file 5: Figure S4, and Additional file 6: Table S2). However, editing sites that are temperature-sensitive are also more variable across species, suggesting that labile RNA structures are under less selective pressure.
One mechanism through which temperature affects RNA editing is through dADAR concentration. We found that dADAR concentration is stable between 10°C and 20°C, but significantly decreases at 30°C (Figure 3), as does RNA editing for many sites (see Additional file 2: Figure S2A). However, levels of the auto-editing isoforms of the dADAR protein do not show differential sensitivity to temperature (Figure 5). The editing enzyme is likely to play an important role in adaptation to temperature; Ma et al. discovered that the hypnos-2 mutant allele of dadar resulted in a loss of editing and increased susceptibility to heat-shock . It is possible that decreased dADAR expression at 30°C is adaptive outside of the direct impact on widespread or specific editing. For example, ADAR is also implicated in the RNAi pathway. ADAR is known to edit RNAi precursors and also competes with the RNAi machinery for double-stranded RNA substrates [39,40]. Further, ADARs affect gene expression through heterochromatic gene silencing , providing a possible evolutionary significance for the sensitivity of dADAR expression to temperature. Thus, temperature alterations can change gene expression through dADAR concentration.
To test our overall hypothesis, that the responsiveness of RNA editing to temperature is mediated in part by RNA structure, we tested several knock-in mutations in the paralytic locus . These mutants were designed to probe and confirm secondary and tertiary RNA base pairing interactions directing the editing of three particular adenosines. In fact, these three adenosines in this unique and complex RNA structure are among the most temperature sensitive ADAR editing sites we observed (Figure 1B, Additional file 3: Table S1). Our mutations represent bidirectional changes to both secondary (DCS, Figure 6) and tertiary (pseudoknot, Additional file 13: Figure S10) structural elements.
Interestingly, in this particular transcript, editing level appears to be responsive to both structure and temperature in an additive manner (Figure 6D, Additional file 13: Figure S10E), suggesting levels of regulation at both the molecular and abiotic levels. Combining mutations designed to restore the pseudoknot structure does not restore the thermo-sensitivity of editing to loxP control levels. This is likely due to altered stability of the RNA tertiary structure: we previously showed that combining the Loop > α and Dock > α’ mutations rescues editing at paralytic site 1, but also results in increased editing at all three sites compared to loxP . This mutational pair (α-α’) recreate a kissing loop interaction that is substantially stabilized by additional hydrogen bonds and potentially increased stacking energies, suggesting that slightly altered stability of the structure alters both absolute editing level as well as thermo-sensitivity of the whole multi-duplex structural complex encompassing the three paralytic sites.
The most conclusive results from these mutants suggest the stabilization of editing levels against temperature by mutations that increase editing, such as the DCS zip and Loop/Dock > α/α’ mutations. Both of these mutations increase editing at all three sites and also result in less thermo-sensitive editing profiles for paralytic sites (Figures 6 and Additional file 13: Figure S10). These mutations likely increase the thermo-stability of secondary and tertiary structures around these sites by adding base pairing interactions or stabilizing base pairs. However, this observation is specific to these particular sites in paralytic, in the most complex RNA tertiary structure yet described for a dADAR substrate, and is not necessarily representative of global trends.
RNA editing is intricately linked to splicing . Splicing often occurs co-transcriptionally  and many of the RNA structures that direct editing involve the participation of intronic cis elements. ADAR must, therefore, edit before RNA structures are resolved and removed to allow transcript splicing. For example, the paralytic HR mutation DCS zip causes the splice donor to be sequestered in an additional secondary structure. This results in increased editing at all three sites in paralytic, but also results in a splicing defect and a temperature-sensitive phenotype . In the present study, we chose PCR primers in exons adjacent to the editing sites of interest and did not observe any temperature-dependent alternative spliceoforms after RT-PCR (see Additional file 14: Figure S11). These data reflect splicing decisions local to the editing sites of interest, and long-range splicing may be affected. Additionally, temperature-sensitive ADAR recruitment may compete with splicing machinery, especially since the structures that direct editing are often formed between exon and downstream intron, and many editing sites are located near exon-intron boundaries .
The temperature range studied here is biologically relevant to Drosophila, and an individual animal is likely to experience the full 20°C temperature range many times during a one to two month lifetime on many temporal scales. Therefore, poikilothermic Drosophila require mechanisms of rapid and reversible cellular adjustment, especially at the level of the critical neuronal proteome, considering the intrinsic temperature sensitivity of neuronal signaling. Although the present study does not address whether or not the response of editing to temperature is adaptive in Drosophila, editing in a number of the potassium channel transcript sites in this study is known to affect channel function.
For example, five highly edited sites in the shab potassium channel transcript affect the voltage-dependence and kinetics of the resulting channels . Characterization of nine Shaker potassium channel isoforms resulting from editing at four sites in the shaker transcript revealed functional epistasis of edited residues on channel inactivation , and gating currents of the Shaker channel are known to be highly temperature sensitive . Finally, editing of four adenosines in the ether-a-go-go potassium channel transcript results in altered channel activation and inactivation kinetics, and editing at a single site alters channel sensitivity to extracellular Mg2+ concentrations . Based on these data, it seems likely that flies living at a different temperature for more than a day or two would be expressing channels altered in their editing profile, and altered in their neuronal signaling properties. We envision that such environmental modulation of channel properties plays a significant and possibly adaptive role in tuning channel kinetics in response to temperature.
Additional intriguing evidence for the role of editing in temperature adaptation exists in the literature. The timing of Drosophila embryogenesis is affected by temperature , and dADAR is expressed in the developing embryo in multiple tissues, including the developing nervous system . Although editing is not known to occur early in development at substantial levels, dADAR may play an additional cellular role in temperature-dependent developmental timing. For example, the relative timing of developmental events scales across temperatures and is constant between many species .
Garrett and Rosenthal have suggested that, due to its unique molecular properties, RNA editing may act as a mechanism for cold adaptation and acclimation because the conversion of an adenosine into an inosine often results in the coding shift to a smaller amino acid side chain . A prevailing hypothesis in the field suggests that enzyme isoforms with increased flexibility around an active site require a lower enthalpy of the rate-limiting step of catalysis, leading to a decrease in cold sensitivity. Because editing often results in amino acid substitutions, changes at the RNA level via editing could affect temperature sensitivity at the protein level.
Our survey of 54 editing sites represented a broad sample of the amino acid changes possible through RNA editing ; the current analysis includes examples of all 17 amino acids created or destroyed via editing. However, we did not see any overall evidence of edited amino acid directionality due to temperature (see Additional file 15: Figure S12). While our data do not support the idea that editing is an adaptive response to temperature by increasing active site lability, as suggested by Garrett and Rosenthal, they also do not preclude this possibility. Further studies in diverse organisms and ADAR targets are necessary.
RNA structures possess the ability to respond to biotic signals, including molecular ligands, and abiotic signals, such as temperature, regulating diverse cellular tasks such as transcriptional regulation and protein synthesis. These structures are inexorably linked to RNA processing events, including RNA editing. Here, we present a survey of the response of RNA editing to temperature alterations in Drosophila. Our results indicate multiple molecular mechanisms of temperature response, including RNA structural stability and dADAR concentration. We found that individual editing sites behave very differently in response to a range of temperatures, suggesting that editing of some adenosines is predominantly influenced by RNA structure, while others are governed by additive mechanisms including enzyme concentration. These data present the first evidence supporting the existence of RNA editing ‘thermometers’, which directly regulate editing in response to environmental temperature. It will be interesting to determine the effects of temperature and other environmental cues on the other functions of ADAR as a regulator of heterochromatic gene silencing .
We raised wild type (Canton-S) Drosophila at 25°C in humidity-controlled incubators with 12-hour light/dark cycles. As age is known to affect editing , we collected animals daily as they ecclosed and shifted newly emerged animals to one of three temperatures: 10°C, 20°C and 30°C, using temperature, humidity and light-controlled (12-hour light/dark) incubators. From experience we know that flies can live within this temperature range for 72 hours, although they are most often subjected to natural temperatures between 20°C and 30°C. We allowed these animals to remain in the incubators for 72 hours, and then immediately froze flies from this population for long-term storage at −80°C.
Additionally, we raised and collected D. erecta, D. sechellia, D. simulans and D. yakuba in the same manner as D. melanogaster above. We did the same with mutant D. melanogaster in which the endogenous dadar gene had been HA-tagged using homologous recombination, as well as the dADAR hypomorph , the dADARS and dADARG engineered animals  and engineered Drosophila with mutations in the paralytic locus .
RNA editing analyses
For each unique species/genotype/temperature sample, we extracted RNA from male fly heads (N = 15 to 20) using TriReagent (Molecular Research Center, Inc., Cincinnati, OH). We chose to use only males to control for sex-related editing changes. For each sample, three separately extracted RNA biological replicates were used. We amplified cDNAs via RT-PCR using random primers, performed PCR using target-specific primers (two per technical PCR replicate per RNA sample) and electrophoresed samples on agarose gels. We cleaned the PCR products using Wizard Gel and PCR Cleanup Kits (Promega, Madison, WI) and detected editing by Sanger sequencing (University of Wisconsin Biotechnology Center, Madison, WI, USA). For species-specific editing, we designed primers to regions of high sequence conservation.
In silico editing quantification
We wrote a program to calculate the ratio of mixed peak areas in chromatograms generated by Sanger sequencing. This program replicates and replaces the procedure that was previously done by hand , generating more accurate, consistent and timely results that correlate highly with those generated by hand. For each mixed chromatogram peak, the program truncates the curves to mitigate the influence from adjacent curves. It then uses linear interpolation to complete parts of the curve missing due to graphical artifacts and linear extrapolation to reconstruct the truncated tails of the curves; higher-moment extrapolation showed little effect on the results in practice. Finally, it computes the area underneath each curve to calculate the mixed peak ratio of the curves, which is a number between 0 and 1 corresponding to the editing level G/(G + A).
Western blot analysis
We prepared whole-head lysates from 30 male heads per 100 μl radioimmunoprecipitation assay (RIPA) buffer (150 mM NaCl, 1.0% NP-40, 0.5% sodium deoxycholate, 0.1% SDS, 50 mM Tris–HCl, pH 8.0), including 1 mM dithiothreitol (DTT), 1 mM phenylmethanesulfonylfluoride (PMSF) and protease inhibitor mixed tablet (Roche). We mixed approximately 25 μg protein per sample with 6 X sample buffer containing SDS and β-mercaptoethanol and loaded samples on 12.5% SDS-PAGE gels, which we then transferred to a 0.2 μm pore size polyvinylidene difluoride membrane (Bio-Rad). Primary antibodies included anti-actin (Millipore; MAB1501; 1:40,000), and anti-HA (Covance; MMS-101P; 1:750). The secondary antibody used was horseradish peroxidase (HRP)-conjugated goat anti-mouse (Abcam; ab5930; 1:6,000). Both technical replicates and biological replicates were performed for each analysis and are presented as supplemental figures. We performed densitometry using ImageJ analysis software (US National Institutes of Health, Bethesda, MD, USA).
For editing analyses, we performed one-way analyses of variance ANOVAs (α = 0.05) followed by Dunnett post-hoc tests in which editing at 10°C and 30°C was compared to editing at 20°C. Exact P-values are noted in supplemental tables (see Additional file 3: Table S1, Additional file 6: Table S2, Additional file 12: Table S3) and on graphs where possible. To compare curve shapes, a population of slopes was calculated from biological and statistical replicates. The populations of slopes between 10°C and 20°C, and 20°C and 30°C for all mutations/species were compared to the appropriate controls (D. melanogaster for Figures 2 and Additional file 4: Figure S3, and loxP for Figures 6 and Additional file 13: Figure S10) using one-way ANOVAs (α = 0.05) followed by Dunnett post-hoc tests.
We thank the laboratory of Dr. Stephen Helfand for use of incubators and Dr. James Jepson for advice and the dadar hypomorph editing data. We thank Patricia Santos for preliminary studies and Dr. Nadia Singh for additional research support. We thank Dr. Erica Larschan, Brian Jones and Kasia Sierzputowska for manuscript edits.
- Bonetti A, Carninci P. Structure by sequencing: discovery of hidden messages in the noncoding RNA fraction. Mol Cell. 2012;48:153–5.PubMedView ArticleGoogle Scholar
- Wan Y, Kertesz M, Spitale RC, Segal E, Chang HY. Understanding the transcriptome through RNA structure. Nat Rev Genet. 2011;12:641–55.PubMedView ArticleGoogle Scholar
- Narberhaus F, Waldminghaus T, Chowdhury S. RNA thermometers. FEMS Microbiol Rev. 2006;30:3–16.PubMedView ArticleGoogle Scholar
- Krajewski SS, Narberhaus F. Temperature-driven differential gene expression by RNA thermosensors. Biochim Biophys Acta. 1839;2014:978–88.Google Scholar
- Li S, Breaker RR. Eukaryotic TPP riboswitch regulation of alternative splicing involving long-distance base pairing. Nucleic Acids Res. 2013;41:3022–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Rieder LE, Reenan RA. The intricate relationship between RNA structure, editing, and splicing. Semin Cell Dev Biol. 2012;23:281–8.PubMedView ArticleGoogle Scholar
- Reenan RA. Molecular determinants and guided evolution of species-specific RNA editing. Nature. 2005;434:409–13.PubMedView ArticleGoogle Scholar
- Rieder LE, Staber CJ, Hoopengardner B, Reenan RA. Tertiary structural elements determine extent and specificity of messenger RNA editing. Nat Commun. 2013;4:2232.PubMedView ArticleGoogle Scholar
- Savva YA, Jepson JE, Sahin A, Sugden AU, Dorsky JS, Alpert L, et al. Auto-regulatory RNA editing fine-tunes mRNA re-coding and complex behavioufigure 1r in Drosophila. Nat Commun. 2012;3:790.PubMed CentralPubMedView ArticleGoogle Scholar
- Savva YA, Rieder LE, Reenan RA. The ADAR protein family. Genome Biol. 2012;13:252.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Q, Khillan J, Gadue P, Nishikura K. Requirement of the RNA editing deaminase ADAR1 gene for embryonic erythropoiesis. Science. 2000;290:1765–8.PubMedView 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
- 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–49.PubMedView ArticleGoogle Scholar
- Tonkin LA, Saccomanno L, Morse DP, Brodigan T, Krause M, Bass BL. RNA editing by ADARs is important for normal behavior in Caenorhabditis elegans. EMBO J. 2002;21:6025–35.PubMed CentralPubMedView ArticleGoogle Scholar
- Kurihara-Yonemoto S, Kubo T. Increased accumulation of intron-containing transcripts in rice mitochondria caused by low temperature: is cold-sensitive RNA editing implicated? Curr Genet. 2010;56:529–41.PubMedView ArticleGoogle Scholar
- Karcher D, Bock R. Temperature sensitivity of RNA editing and intron splicing reactions in the plastid ndhB transcript. Curr Genet. 2002;41:48–52.PubMedView ArticleGoogle Scholar
- Stocker J, Huang HW, Wang HM, Chang HW, Chiu CC, Cho CL, et al. Reduction of RNA A-to-I editing in Drosophila acclimated to heat shock. Kaohsiung J Med Sci. 2013;29:478–83.PubMedView ArticleGoogle Scholar
- Tian N, Yang Y, Sachsenmaier N, Muggenhumer D, Bi J, Waldsich C, et al. A structural determinant required for RNA editing. Nucleic Acids Res. 2011;39:5669–81.PubMed CentralPubMedView ArticleGoogle Scholar
- Dahabieh MS, Samanta D, Brodovitch JC, Frech C, O’Neill MA, Pinto BM. Sequence-dependent structural dynamics of primate adenosine-to-inosine editing substrates. Chembiochem. 2012;13:2714–21.PubMedView ArticleGoogle Scholar
- Garrett S, Rosenthal JJ. RNA editing underlies temperature adaptation in K+ channels from polar octopuses. Science. 2012;335:848–51.PubMed CentralPubMedView ArticleGoogle Scholar
- Pinto Y, Cohen HY, Levanon EY. Mammalian conserved ADAR targets comprise only a small fragment of the human editosome. Genome Biol. 2014;15:R5.PubMed CentralPubMedView ArticleGoogle Scholar
- Economos AC, Lints FA. Developmental temperature and life span in Drosophila melanogaster. II Oscillating temperature. Gerontology. 1986;32:28–36.PubMedView ArticleGoogle Scholar
- Carmi S, Borukhov I, Levanon EY. Identification of widespread ultra-edited human RNAs. PLoS Genet. 2011;7:e1002317.PubMed CentralPubMedView ArticleGoogle Scholar
- Hood JL, Emeson RB. Editing of neurotransmitter receptor and ion channel RNAs in the nervous system. Curr Top Microbiol Immunol. 2012;353:61–90.PubMed CentralPubMedGoogle Scholar
- Schmauss C, Howe JR. RNA editing of neurotransmitter receptors in the mammalian brain. Sci STKE. 2002;2002:pe26.PubMedGoogle Scholar
- Montgomery JC, Macdonald JA. Effects of temperature on nervous-system - implications for behavioral performance. Am J Physiol. 1990;259:R191–6.PubMedGoogle Scholar
- Dillon ME, Wang G, Garrity PA, Huey RB. Review: thermal preference in Drosophila. J Therm Biol. 2009;34:109–19.PubMed CentralPubMedView ArticleGoogle Scholar
- Ryan MY, Maloney R, Fineberg JD, Reenan RA, Horn R. RNA editing in eag potassium channels: biophysical consequences of editing a conserved S6 residue. Channels (Austin). 2012;6:443–52.View ArticleGoogle Scholar
- St Laurent G, Tackett MR, Nechkin S, Shtokalo D, Antonets D, Savva YA, et al. Genome-wide analysis of A-to-I RNA editing by single-molecule sequencing in Drosophila. Nat Struct Mol Biol. 2013;20:1333–9.PubMedView ArticleGoogle Scholar
- Tamura K, Subramanian S, Kumar S. Temporal patterns of fruit fly (Drosophila) evolution revealed by mutation clocks. Mol Biol Evol. 2004;21:36–44.PubMedView ArticleGoogle Scholar
- 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:2183–93.PubMed CentralPubMedView 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.Google Scholar
- Staber CJ, Gell S, Jepson JE, Reenan RA. Perturbing A-to-I RNA editing using genetics and homologous recombination. Methods Mol Biol. 2011;718:41–73.PubMedView ArticleGoogle Scholar
- Hoffmann AA. Physiological climatic limits in Drosophila: patterns and implications. J Exp Biol. 2010;213:870–80.PubMedView ArticleGoogle Scholar
- Lachaise D, Silvain JF. How two Afrotropical endemics made two cosmopolitan human commensals: the Drosophila melanogaster-D. simulans palaeogeographic riddle. Genetica. 2004;120:17–39.PubMedView ArticleGoogle Scholar
- Rio B, Couturier G, Lemeunier F, Lachaise D. Evolution of a seasonal specialization in Drosophila-Erecta (Dipt-Drosophilidae). Annales De La Societe Entomologique De France. 1983;19:235–48.Google Scholar
- Lachaise D, Cariou ML, David JR, Lemeunier F, Tsacas L, Ashburner M. Historical biogeography of the Drosophila-Melanogaster species subgroup. Evol Biol. 1988;22:159–225.Google Scholar
- Ma EB, Gu XQ, Wu XH, Xu T, Haddad GG. Mutation in pre-mRNA adenosine deaminase markedly attenuates neuronal tolerance to O-2 deprivation in Drosophila melanogaster (vol 107, pg 685, 2001). J Clin Invest. 2001;107:1203–1203.View ArticleGoogle Scholar
- Wu D, Lamm AT, Fire AZ. Competition between ADAR and RNAi pathways for an extensive class of RNA targets. Nat Struct Mol Biol. 2011;18:1094–101.PubMed CentralPubMedView ArticleGoogle Scholar
- Heale BS, Keegan LP, McGurk L, Michlewski G, Brindle J, Stanton CM, et al. Editing independent effects of ADARs on the miRNA/siRNA pathways. EMBO J. 2009;28:3145–56.PubMed CentralPubMedView ArticleGoogle Scholar
- Savva YA, Jepson JE, Chang YJ, Whitaker R, Jones BC, St Laurent G, et al. RNA editing regulates transposon-mediated heterochromatic gene silencing. Nat Commun. 2013;4:2745.PubMed CentralPubMedView ArticleGoogle Scholar
- Khodor YL, Rodriguez J, Abruzzi KC, Tang CH, Marr 2nd MT, Rosbash M. Nascent-seq indicates widespread cotranscriptional pre-mRNA splicing in Drosophila. Genes Dev. 2011;25:2502–12.PubMed CentralPubMedView ArticleGoogle Scholar
- Ryan MY, Maloney R, Reenan R, Horn R. Characterization of five RNA editing sites in Shab potassium channels. Channels (Austin). 2008;2:202–9.View ArticleGoogle Scholar
- Ingleby L, Maloney R, Jepson J, Horn R, Reenan R. Regulated RNA editing and functional epistasis in Shaker potassium channels. J Gen Physiol. 2009;133:17–27.PubMed CentralPubMedView ArticleGoogle Scholar
- Rodriguez BM, Sigg D, Bezanilla F. Voltage gating of Shaker K+ channels. The effect of temperature on ionic and gating currents. J Gen Physiol. 1998;112:223–42.PubMed CentralPubMedView ArticleGoogle Scholar
- Powsner L. The effects of temperature on the durations of the developmental stages of drosophila melanogaster. Physiol Zool. 1935;8:474–520.Google Scholar
- Chen J, Lakshmi GG, Hays DL, McDowell KM, Ma E, Vaughn JC. Spatial and temporal expression of dADAR mRNA and protein isoforms during embryogenesis in Drosophila melanogaster. Differentiation. 2009;78:312–20.PubMedView ArticleGoogle Scholar
- Kuntz SG, Eisen MB. Drosophila embryogenesis scales uniformly across temperature in developmentally diverse species. PLoS Genet. 2014;10:e1004293.PubMed CentralPubMedView ArticleGoogle Scholar
- Garrett SC, Rosenthal JJ. A role for A-to-I RNA editing in temperature adaptation. Physiology (Bethesda). 2012;27:362–9.View ArticleGoogle Scholar
- Ma E, Tucker MC, Chen Q, Haddad GG. Developmental expression and enzymatic activity of pre-mRNA deaminase in Drosophila melanogaster. Brain Res Mol Brain Res. 2002;102:100–4.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.