Genetic variation in senescence marker protein-30 is associated with natural variation in cold tolerance in Drosophila

Metabolomes of bumble bees reared in common garden conditions suggest constitutive differences in energy and toxin metabolism across populations

Cold tolerance of ectotherms can vary strikingly among species and populations. Variation in cold tolerance can reflect differences in genomes and transcriptomes that confer cellular-level protection from cold; additionally, shifts in protein function and abundance can be altered by other cellular constituents as cold-exposed insects often have shifts in their metabolomes. Even without a cold challenge, insects from different populations may vary in cellular composition that could alter cold tolerance, but investigations of constitutive differences in metabolomes across wild populations remain rare. To address this gap, we reared Bombus vosnesenskii queens collected from Oregon and California (USA) that differ in cold tolerance (CTmin = -6 °C and 0 °C, respectively) in common garden conditions, and measured offspring metabolomes using untargeted LC-MS/MS. Oregon bees had higher levels of metabolites associated with carbohydrate (sorbitol, lactitol, maltitol, and sorbitol-6-phosphate) and amino acid (hydroxyproline, ornithine, and histamine) metabolism. Exogenous metabolites, likely derived from the diet, also varied between Oregon and California bees, suggesting population-level differences in toxin metabolism. Overall, our results reveal constitutive differences in metabolomes for bumble bees reared in common garden conditions from queens collected in different locations despite no previous cold exposure.

Transcriptome analysis of Drosophila melanogaster laboratory strains of different geographical origin after long-term laboratory maintenance

Positive selection may be the main factor of the between-population divergence in gene expression. Expression profiles of two Drosophila melanogaster laboratory strains of different geographical origin and long-term laboratory maintenance were analyzed using microchip arrays encompassing probes for 18,500 transcripts. The Russian strain D18 and the North American strain Canton-S were compared. A set of 223 known or putative genes demonstrated significant changes in expression levels between these strains. Differentially expressed genes (DEG) were enriched in response to DDT (p = .0014), proteolysis (p = 2.285E-5), transmembrane transport (p = 1.03E-4), carbohydrate metabolic process (p = .0317), protein homotetramerization (p = .0444), and antibacterial humoral response (p = 425E-4). The expression in subset of genes from different categories was verified by qRT-PCR. Analysis of transcript abundance between Canton-S and D18 strains allowed to select several genes to estimate their participation in latitude adaptation. Expression of selected genes was analyzed in five D. melanogaster lines of different geographic origins by qRT-PCR, and we found two candidate genes that may be associated with latitude adaptation in adult flies-smp-30 and Cda9. Quite possible that several alleles of these genes may be important for insect survival in the environments of global warming. It is interesting that the number of genes involved in local adaptation demonstrates expression level appropriate to their geographical origin even after decades of laboratory maintenance.

Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel

Insecticide resistance is considered a classic model of microevolution, where a strong selective agent is applied to a large natural population, resulting in a change in frequency of alleles that confer resistance. While many insecticide resistance variants have been characterized at the gene level, they are typically single genes of large effect identified in highly resistant pest species. In contrast, multiple variants have been implicated in DDT resistance in Drosophila melanogaster; however, only the Cyp6g1 locus has previously been shown to be relevant to field populations. Here we use genome-wide association studies (GWAS) to identify DDT-associated polygenes and use selective sweep analyses to assess their adaptive significance. We identify and verify two candidate DDT resistance loci. A largely uncharacterized gene, CG10737, has a function in muscles that ameliorates the effects of DDT, while a putative detoxifying P450, Cyp6w1, shows compelling evidence of positive selection.

Application of Response Surface Methods To Determine Conditions for Optimal Genomic Prediction

An epistatic genetic architecture can have a significant impact on prediction accuracies of genomic prediction (GP) methods. Machine learning methods predict traits comprised of epistatic genetic architectures more accurately than statistical methods based on additive mixed linear models. The differences between these types of GP methods suggest a diagnostic for revealing genetic architectures underlying traits of interest. In addition to genetic architecture, the performance of GP methods may be influenced by the sample size of the training population, the number of QTL, and the proportion of phenotypic variability due to genotypic variability (heritability). Possible values for these factors and the number of combinations of the factor levels that influence the performance of GP methods can be large. Thus, efficient methods for identifying combinations of factor levels that produce most accurate GPs is needed. Herein, we employ response surface methods (RSMs) to find the experimental conditions that produce the most accurate GPs. We illustrate RSM with an example of simulated doubled haploid populations and identify the combination of factors that maximize the difference between prediction accuracies of best linear unbiased prediction (BLUP) and support vector machine (SVM) GP methods. The greatest impact on the response is due to the genetic architecture of the population, heritability of the trait, and the sample size. When epistasis is responsible for all of the genotypic variance and heritability is equal to one and the sample size of the training population is large, the advantage of using the SVM method vs. the BLUP method is greatest. However, except for values close to the maximum, most of the response surface shows little difference between the methods. We also determined that the conditions resulting in the greatest prediction accuracy for BLUP occurred when genetic architecture consists solely of additive effects, and heritability is equal to one.

Cold acclimation wholly reorganizes the Drosophila melanogaster transcriptome and metabolome

Cold tolerance is a key determinant of insect distribution and abundance, and thermal acclimation can strongly influence organismal stress tolerance phenotypes, particularly in small ectotherms like Drosophila. However, there is limited understanding of the molecular and biochemical mechanisms that confer such impressive plasticity. Here, we use high-throughput mRNA sequencing (RNA-seq) and liquid chromatography – mass spectrometry (LC-MS) to compare the transcriptomes and metabolomes of D. melanogaster acclimated as adults to warm (rearing) (21.5 °C) or cold conditions (6 °C). Cold acclimation improved cold tolerance and led to extensive biological reorganization: almost one third of the transcriptome and nearly half of the metabolome were differentially regulated. There was overlap in the metabolic pathways identified via transcriptomics and metabolomics, with proline and glutathione metabolism being the most strongly-supported metabolic pathways associated with increased cold tolerance. We discuss several new targets in the study of insect cold tolerance (e.g. dopamine signaling and Na⁺-driven transport), but many previously identified candidate genes and pathways (e.g. heat shock proteins, Ca²⁺ signaling, and ROS detoxification) were also identified in the present study, and our results are thus consistent with and extend the current understanding of the mechanisms of insect chilling tolerance.

Natural Populations of Drosophila melanogaster Reveal Features of an Uncharacterized Circadian Property: The Lower Temperature Limit of Rhythmicity

Most cyclic biological processes are under control of a circadian molecular timing system that synchronizes these phenomena to the 24-h day. One generic property of circadian-controlled processes is that they operate within a specific temperature range, below which the manifestation of rhythm ceases. Little is known about the evolutionary relevance of the lower temperature limit of rhythmicity or about the mechanism underlying the loss of overt circadian behavior below this lower limit, especially in one model organism of chronobiology, Drosophila melanogaster. Natural populations of Drosophila are evolving under divergent selection pressures and so provide a source of diversity necessary to address these issues. Using lines derived from African populations, we find that there is natural variation in the expression of rhythmic behavior under low-temperature conditions. We found evidence that this variability is evolutionarily relevant at extremely low temperature (12 °C) because high-altitude populations exhibit selection for locally adapted genomes that contribute to rhythmic behavior. Lines resistant to 15 °C show an additional layer of diversity in their response to temperature extremes because some lines are resistant to low temperature (15 °C) only, whereas others are cross-resistant to high and low temperature (15 °C and 30 °C). Genetic analysis of one cold-resistant circadian line at 15 °C reveals that the phenotype maps to the X-chromosome but not to the core clock genes, per and sgg. Analysis of the central clock cells of this line reveals that maintenance of rhythm is associated with robust clock function, which is compromised in a standard laboratory strain. These data indicate that the cold-resistant circadian phenotype is clock based. This study highlights the importance of using natural populations to inform us of the basic features of circadian traits, especially those that might be under temperature-based selection.

Increased abundance of frost mRNA during recovery from cold stress is not essential for cold tolerance in adult Drosophila melanogaster

Frost (Fst) is a candidate gene associated with the response to cold in Drosophila melanogaster because Fst mRNA accumulation increases during recovery from low temperature exposure. We investigated the contribution of Fst expression to chill-coma recovery time, acute cold tolerance and rapid cold hardening (RCH) in adult D. melanogaster by knocking down Fst mRNA expression using GAL4/UAS-mediated RNA interference. In this experiment, four UAS-Fst and one tubulin-GAL4 lines were used. We predicted that if Fst is essential for cold tolerance phenotypes, flies with low Fst mRNA levels should be less cold tolerant than flies with normal levels of cold-induced Fst mRNA. Cold-induced Fst abundance and recovery time from chill-coma were not negatively correlated in male or female flies. Survival of 2 h exposures to sub-zero temperatures in Fst knockdown lines was not lower than that in a control line. Moreover, a low temperature pretreatment increased survival of severe cold exposure in flies regardless of Fst abundance level during recovery from cold stress, suggesting that Fst expression is not essential for RCH. Thus, cold-induced Fst accumulation is not essential for cold tolerance measured as chill-coma recovery time, survival to acute cold stress and RCH response in adult D. melanogaster.

Proteomic profiling of thermal acclimation in Drosophila melanogaster

Thermal acclimation drastically alters thermotolerance of ectotherms, but the mechanisms determining this plastic response are not fully understood. The present study investigates the proteomic response (2D-DIGE) of adult Drosophila melanogaster acclimated at 11, 25 or 31 °C. As expected 11 °C-acclimation improved cold tolerance and 31 °C-acclimation improved heat tolerance. We hypothesized that the marked organismal responses to acclimation could be detected at the proteomic level assuming that changes in the abundance of specific proteins are linked to the physiological changes underlying the phenotypic response. The 31 °C-acclimated flies displayed a particular divergent proteomic profile where molecular chaperones made up a large number of the proteins that were modulated during heat acclimation. Many other proteins showed significant modulation during acclimation including proteins involved in iron ion and cell redox homeostasis, carbohydrate and energy metabolism, chromatin remodeling and translation, and contractile machinery. Interestingly the changes in protein abundance were often unrelated to transcriptional activity of the genes coding for the proteins, except for the most strongly expressed proteins (e.g. Hsp70). The 11 °C-acclimation evoked weak proteomic response despite the marked effect on the organismal phenotype. Thus the acquired cold tolerance observed here may involve regulatory process such as posttranslational regulation rather than de novo protein synthesis.

Reestablishment of ion homeostasis during chill-coma recovery in the cricket Gryllus pennsylvanicus

The time required to recover from cold-induced paralysis (chill-coma) is a common measure of insect cold tolerance used to test central questions in thermal biology and predict the effects of climate change on insect populations. The onset of chill-coma in the fall field cricket (Gryllus pennsylvanicus, Orthoptera: Gryllidae) is accompanied by a progressive drift of Na(+) and water from the hemolymph to the gut, but the physiological mechanisms underlying recovery from chill-coma are not understood for any insect. Using a combination of gravimetric methods and atomic absorption spectroscopy, we demonstrate that recovery from chill-coma involves a reestablishment of hemolymph ion content and volume driven by removal of Na(+) and water from the gut. Recovery is associated with a transient elevation of metabolic rate, the time span of which increases with increasing cold exposure duration and closely matches the duration of complete osmotic recovery. Thus, complete recovery from chill-coma is metabolically costly and encompasses a longer period than is required for the recovery of muscle potentials and movement. These findings provide evidence that physiological mechanisms of hemolymph ion content and volume regulation, such as ion-motive ATPase activity, are instrumental in chill-coma recovery and may underlie natural variation in insect cold tolerance.

Genomic variation in natural populations of Drosophila melanogaster

This report of independent genome sequences of two natural populations of Drosophila melanogaster (37 from North America and 6 from Africa) provides unique insight into forces shaping genomic polymorphism and divergence. Evidence of interactions between natural selection and genetic linkage is abundant not only in centromere- and telomere-proximal regions, but also throughout the euchromatic arms. Linkage disequilibrium, which decays within 1 kbp, exhibits a strong bias toward coupling of the more frequent alleles and provides a high-resolution map of recombination rate. The juxtaposition of population genetics statistics in small genomic windows with gene structures and chromatin states yields a rich, high-resolution annotation, including the following: (1) 5'- and 3'-UTRs are enriched for regions of reduced polymorphism relative to lineage-specific divergence; (2) exons overlap with windows of excess relative polymorphism; (3) epigenetic marks associated with active transcription initiation sites overlap with regions of reduced relative polymorphism and relatively reduced estimates of the rate of recombination; (4) the rate of adaptive nonsynonymous fixation increases with the rate of crossing over per base pair; and (5) both duplications and deletions are enriched near origins of replication and their density correlates negatively with the rate of crossing over. Available demographic models of X and autosome descent cannot account for the increased divergence on the X and loss of diversity associated with the out-of-Africa migration. Comparison of the variation among these genomes to variation among genomes from D. simulans suggests that many targets of directional selection are shared between these species.

Genetics of climate change adaptation

The rapid rate of current global climate change is having strong effects on many species and, at least in some cases, is driving evolution, particularly when changes in conditions alter patterns of selection. Climate change thus provides an opportunity for the study of the genetic basis of adaptation. Such studies include a variety of observational and experimental approaches, such as sampling across clines, artificial evolution experiments, and resurrection studies. These approaches can be combined with a number of techniques in genetics and genomics, including association and mapping analyses, genome scans, and transcription profiling. Recent research has revealed a number of candidate genes potentially involved in climate change adaptation and has also illustrated that genetic regulatory networks and epigenetic effects may be particularly relevant for evolution driven by climate change. Although genetic and genomic data are rapidly accumulating, we still have much to learn about the genetic architecture of climate change adaptation.

Photoperiodic regulation of cold tolerance and expression levels of regucalcin gene in Drosophila montana

Temperature-induced plasticity of cold tolerance has been reported in many insect species, but cold tolerance can also be affected by changes in day (or night) length. In the present study we elucidate the direct and indirect effects of photoperiod on the cold tolerance of females of two Drosophila montana strains--one which possesses a robust photoperiodic diapause and another which does not. In the diapause-strain the time needed for recovery from chill coma showed a positive correlation with day length, but diapause itself played only a minor role in photoperiodic acclimation. The strain that was not able to enter to diapause as a response to day length also lacked photoperiodic cold acclimation ability indicating that this strain has deficiencies in its photoperiodic time measurement system. In the diapause-strain, the expression level of regucalcin gene was more than two times higher in diapausing than in non-diapausing females maintained in a single photoperiod, but day length per se did not cause significant changes in expression levels of this gene in either of the strains. In the non-diapausing strain this gene showed no expression changes in any comparison. Overall, the study shows that a decrease in day length can induce cold acclimation in D. montana, while changes in regucalcin expression are linked with photoperiodic diapause.

A comparative study of the short term cold resistance response in distantly related Drosophila species: the role of regucalcin and frost

The molecular basis of short term cold resistance (indexed as chill-coma recovery time) has been mostly addressed in D. melanogaster, where candidate genes (Dca (also known as smp-30) and Frost (Fst)) have been identified. Nevertheless, in Drosophila, the ability to tolerate short term exposure to low temperatures evolved several times independently. Therefore, it is unclear whether variation in the same candidate genes is also responsible for short term cold resistance in distantly related Drosophila species. It should be noted that Dca is a candidate gene for cold resistance in the Sophophora subgenus only, since there is no orthologous gene copy in the Drosophila subgenus. Here we show that, in D. americana (Drosophila subgenus), there is a north-south gradient for a variant at the 5′ non-coding region of regucalcin (a Dca-like gene; in D. melanogaster the proteins encoded by the two genes share 71.9% amino acid identities) but in our D. americana F2 association experiment there is no association between this polymorphism and chill-coma recovery times. Moreover, we found no convincing evidence that this gene is up-regulated after cold shock in both D. americana and D. melanogaster. Size variation in the Fst PEST domain (putatively involved in rapid protein degradation) is observed when comparing distantly related Drosophila species, and is associated with short term cold resistance differences in D. americana. Nevertheless, this effect is likely through body size variation. Moreover, we show that, even at two hours after cold shock, when up-regulation of this gene is maximal in D. melanogaster (about 48 fold expression change), in D. americana this gene is only moderately up-regulated (about 3 fold expression change). Our work thus shows that there are important differences regarding the molecular basis of cold resistance in distantly related Drosophila species.

Molecular basis of adaptive shift in body size in Drosophila melanogaster: functional and sequence analyses of the Dca gene

Latitudinal body size clines in animals conforming to Bergmann's rule occur on many continents but isolating their underlying genetic basis remains a challenge. In Drosophila melanogaster, the gene Dca accounts for approximately 5-10% of the natural wing size variation (McKechnie SW, Blacket MJ, Song SV, Rako L, Carroll X, Johnson TK, Jensen LT, Lee SF, Wee CW, Hoffmann AA. 2010. A clinally varying promoter polymorphism associated with adaptive variation in wing size in Drosophila. Mol Ecol. 19:775-784). We present here functional evidence that Dca is a negative regulator of wing size. A significant negative latitudinal cline of Dca gene expression was detected in synchronized third instar larvae. In addition, we clarified the evolutionary history of the three most common Dca promoter alleles (Dca237-1, Dca237-2, and Dca247) and showed that the insertion allele (Dca247), whose frequency increases with latitude, is associated with larger wing centroid size and higher average cell number in male flies. Finally, we showed that the overall linkage disequilibrium (LD) was low in the Dca promoter and that the insertion/deletion polymorphism that defines the Dca alleles was in strong LD with two other upstream sites. Our results provide strong support that Dca is a candidate for climatic adaptation in D. melanogaster.

Rapid growth reduces cold resistance: evidence from latitudinal variation in growth rate, cold resistance and stress proteins

Background

Physiological costs of rapid growth may contribute to the observation that organisms typically grow at submaximal rates. Although, it has been hypothesized that faster growing individuals would do worse in dealing with suboptimal temperatures, this type of cost has never been explored empirically. Furthermore, the mechanistic basis of the physiological costs of rapid growth is largely unexplored.

Methodology/Principal Finding

Larvae of the damselfly Ischnura elegans from two univoltine northern and two multivoltine southern populations were reared at three temperatures and after emergence given a cold shock. Cold resistance, measured by chill coma recovery times in the adult stage, was lower in the southern populations. The faster larval growth rates in the southern populations contributed to this latitudinal pattern in cold resistance. In accordance with their assumed role in cold resistance, Hsp70 levels were lower in the southern populations, and faster growing larvae had lower Hsp70 levels. Yet, individual variation in Hsp70 levels did not explain variation in cold resistance.

Conclusions/Significance

We provide evidence for a novel cost of rapid growth: reduced cold resistance. Our results indicate that the reduced cold resistance in southern populations of animals that change voltinism along the latitudinal gradient may not entirely be explained by thermal selection per se but also by the costs of time constraint-induced higher growth rates. This also illustrates that stressors imposed in the larval stage may carry over and shape fitness in the adult stage and highlights the importance of physiological costs in the evolution of life-histories at macro-scales.