Genetic variation in the shape of cold‐survival curves in a single fly population suggests potential for selection from climate variability

Impact of Heat Stress on Transposable Element Expression and Derived Small RNAs in Drosophila subobscura

Global warming is forcing insect populations to move and adapt, triggering adaptive genetic responses. Thermal stress is known to alter gene expression, repressing the transcription of active genes, and inducing others, such as those encoding heat shock proteins. It has also been related to the activation of some specific transposable element (TE) families. However, the actual magnitude of this stress on the whole genome and the factors involved in these genomic changes are still unclear. We studied mRNAs and small RNAs in gonads of two Drosophila subobscura populations, considered a good model to study adaptation to temperature changes. In control conditions, we found that a few genes and TE families were differentially expressed between populations, pointing out their putative involvement in the adaptation of populations to their different environments. Under heat stress, sex-specific changes in gene expression together with a trend toward overexpression, mainly of heat shock response-related genes, were observed. We did not observe large changes of TE expression nor small RNA production due to stress. Only population and sex-specific expression changes of some TE families (mainly retrotransposons), or the amounts of siRNAs and piRNAs, derived from specific TE families were observed, as well as the piRNA production from some piRNA clusters. Changes in small RNA amounts and TE expression could not be clearly correlated, indicating that other factors as chromatin modulation could also be involved. This work provides the first whole transcriptomic study including genes, TEs, and small RNAs after a heat stress in D. subobscura.

Molecular Mechanisms of Winter Survival

Winter provides many challenges for insects, including direct injury to tissues and energy drain due to low food availability. As a result, the geographic distribution of many species is tightly coupled to their ability to survive winter. In this review, we summarize molecular processes associated with winter survival, with a particular focus on coping with cold injury and energetic challenges. Anticipatory processes such as cold acclimation and diapause cause wholesale transcriptional reorganization that increases cold resistance and promotes cryoprotectant production and energy storage. Molecular responses to low temperature are also dynamic and include signaling events during and after a cold stressor to prevent and repair cold injury. In addition, we highlight mechanisms that are subject to selection as insects evolve to variable winter conditions. Based on current knowledge, despite common threads, molecular mechanisms of winter survival vary considerably across species, and taxonomic biases must be addressed to fully appreciate the mechanistic basis of winter survival across the insect phylogeny.

Analysis of heat and cold tolerance of a freeze-tolerant soil invertebrate distributed from temperate to Arctic regions: evidence of selection for extreme cold tolerance

Tolerance to thermal extremes is critical for the geographic distributions of ectotherm species, many of which are probably going to be modified by future climatic changes. To predict species distributions it is important to understand the potential of species to adapt to changing thermal conditions. Here, we tested whether the thermal tolerance traits of a common freeze-tolerant potworm were correlated with climatic conditions and if adaptation to extreme cold constrains the evolutionary potential for high temperature tolerance. Further, we tested if evolution of thermal tolerance traits is associated with costs in other fitness traits (body size and reproduction). Lastly, we tested if slopes of temperature-survival curves (i.e., the sensitivity distribution) are related to tolerance itself. Using 24 populations of the potworm, Enchytraeus albidus Henle (Enchytraeidae), collected from a wide range of climatic conditions, we established a common garden experiment in which we determined high and low temperature tolerance (using survival as endpoint), average reproductive output and adult body size. Heat tolerance was not related to environmental temperatures whereas lower lethal temperature was about 10 °C lower in Arctic populations than in populations from temperate regions. Reproduction was not related to environmental temperature, but was negatively correlated with cold tolerance. One explanation for the trade-off between cold tolerance and reproduction could be that the more cold-hardy populations need to channel energy to large glycogen reserves at the expense of less energy expenditure for reproduction. Adult body size was negatively related to environmental temperature. Finally, the slopes of temperature-survival curves were significantly correlated with critical temperature limits for heat and cold tolerance; i.e., slopes increased with thermal tolerance. Our results suggest that relatively heat-sensitive populations possess genetic variation, leaving room for improved heat tolerance through evolutionary processes, which may alleviate the effects of a warmer future climate in the Arctic. On the other hand, we observed relatively narrow sensitivity distributions (i.e., less variation) in the most heat tolerant populations. Taken together, our results suggest that both cold and heat tolerance can only be selected for (and improved) until a certain limit has been reached.

Transcriptomic and functional genetic evidence for distinct ecophysiological responses across complex life cycle stages

Organisms with complex life cycles demonstrate a remarkable ability to change their phenotypes across development, presumably as an evolutionary adaptation to developmentally variable environments. Developmental variation in environmentally sensitive performance, and thermal sensitivity in particular, has been well documented in holometabolous insects. For example, thermal performance in adults and juvenile stages exhibit little genetic correlation (genetic decoupling) and can evolve independently, resulting in divergent thermal responses. Yet, we understand very little about how this genetic decoupling occurs. We tested the hypothesis that genetic decoupling of thermal physiology is driven by fundamental differences in physiology between life stages, despite a potentially conserved Cellular Stress Response. We used RNAseq to compare transcript expression in response to a cold stressor in Drosophila melanogaster larvae and adults and used RNAi (RNA interference) to test whether knocking down nine target genes differentially affected larval and adult cold tolerance. Transcriptomic responses of whole larvae and adults during and following exposure to -5°C were largely unique both in identity of responding transcripts and in temporal dynamics. Further, we analyzed the tissue-specificity of differentially-expressed transcripts from FlyAtlas 2 data, and concluded that stage-specific differences in transcription were not simply driven by differences in tissue composition. In addition, RNAi of target genes resulted in largely stage-specific and sometimes sex-specific effects on cold tolerance. The combined evidence suggests that thermal physiology is largely stage-specific at the level of gene expression, and thus natural selection may be acting on different loci during the independent thermal adaptation of different life stages.

Pool-GWAS on reproductive dormancy in Drosophila simulans suggests a polygenic architecture

The genetic basis of adaptation to different environments has been of long-standing interest to evolutionary biologists. Dormancy is a well-studied adaptation to facilitate overwintering. In Drosophila melanogaster, a moderate number of genes with large effects have been described, which suggests a simple genetic basis of dormancy. On the other hand, genome-wide scans for dormancy suggest a polygenic architecture in insects. In D. melanogaster, the analysis of the genetic architecture of dormancy is complicated by the presence of cosmopolitan inversions. Here, we performed a genome-wide scan to characterize the genetic basis of this ecologically extremely important trait in the sibling species of D. melanogaster, D. simulans that lacks cosmopolitan inversions. We performed Pool-GWAS in a South African D. simulans population for dormancy incidence at 2 temperature regimes (10 and 12°C, LD 10:14). We identified several genes with SNPs that showed a significant association with dormancy (P-value < 1e-13), but the overall modest response suggests that dormancy is a polygenic trait with many loci of small effect. Our results shed light on controversies on reproductive dormancy in Drosophila and have important implications for the characterization of the genetic basis of this trait.

Local thermal environment and warming influence supercooling and drive widespread shifts in the metabolome of diapausing Pieris rapae butterflies

Global climate change has the potential to negatively impact biological systems as organisms are exposed to novel temperature regimes. Increases in annual mean temperature have been accompanied by disproportionate rates of change in temperature across seasons, and winter is the season warming most rapidly. Yet, we know relatively little about how warming will alter the physiology of overwintering organisms. Here, we simulated future warming conditions by comparing diapausing Pieris rapae butterfly pupae collected from disparate thermal environments and by exposing P. rapae pupae to acute and chronic increases in temperature. First, we compared internal freezing temperatures (supercooling points) of diapausing pupae that were developed in common-garden conditions but whose parents were collected from northern Vermont, USA, or North Carolina, USA. Matching the warmer winter climate of North Carolina, North Carolina pupae had significantly higher supercooling points than Vermont pupae. Next, we measured the effects of acute and chronic warming exposure in Vermont pupae and found that warming induced higher supercooling points. We further characterized the effects of chronic warming by profiling the metabolomes of Vermont pupae via untargeted LC-MS metabolomics. Warming caused significant changes in abundance of hundreds of metabolites across the metabolome. Notably, there were warming-induced shifts in key biochemical pathways, such as pyruvate metabolism, fructose and mannose metabolism, and β-alanine metabolism, suggesting shifts in energy metabolism and cryoprotection. These results suggest that warming affects various aspects of overwintering physiology in P. rapae and may be detrimental depending on the frequency and variation of winter warming events. Further research is needed to ascertain the extent to which the effects of warming are felt among a broader set of populations of P. rapae, and among other species, in order to better predict how insects may respond to changes in winter thermal environments.

Natural variation in the transcriptional response of Drosophila melanogaster to oxidative stress

Broadly distributed species must cope with diverse and changing environmental conditions, including various forms of stress. Cosmopolitan populations of Drosophila melanogaster are more tolerant to oxidative stress than those from the species’ ancestral range in sub-Saharan Africa, and the degree of tolerance is associated with an insertion/deletion polymorphism in the 3′ untranslated region of the Metallothionein A (MtnA) gene that varies clinally in frequency. We examined oxidative stress tolerance and the transcriptional response to oxidative stress in cosmopolitan and sub-Saharan African populations of D. melanogaster, including paired samples with allelic differences at the MtnA locus. We found that the effect of the MtnA polymorphism on oxidative stress tolerance was dependent on the genomic background, with the deletion allele increasing tolerance only in a northern, temperate population. Genes that were differentially expressed under oxidative stress included MtnA and other metallothioneins, as well as those involved in glutathione metabolism and other genes known to be part of the oxidative stress response or the general stress response. A gene coexpression analysis revealed further genes and pathways that respond to oxidative stress including those involved in additional metabolic processes, autophagy, and apoptosis. There was a significant overlap among the genes induced by oxidative and cold stress, which suggests a shared response pathway to these two stresses. Interestingly, the MtnA deletion was associated with consistent changes in the expression of many genes across all genomic backgrounds, regardless of the expression level of the MtnA gene itself. We hypothesize that this is an indirect effect driven by the loss of microRNA binding sites within the MtnA 3′ untranslated region.

Evidence for genetic isolation and local adaptation in the field cricket Gryllus campestris

Understanding how species can thrive in a range of environments is a central challenge for evolutionary ecology. There is strong evidence for local adaptation along large-scale ecological clines in insects. However, potential adaptation among neighbouring populations differing in their environment has been studied much less. We used RAD sequencing to quantify genetic divergence and clustering of ten populations of the field cricket Gryllus campestris in the Cantabrian Mountains of northern Spain, and an outgroup on the inland plain. Our populations were chosen to represent replicate high and low altitude habitats. We identified genetic clusters that include both high and low altitude populations indicating that the two habitat types do not hold ancestrally distinct lineages. Using common-garden rearing experiments to remove environmental effects, we found evidence for differences between high and low altitude populations in physiological and life-history traits. As predicted by the local adaptation hypothesis, crickets with parents from cooler (high altitude) populations recovered from periods of extreme cooling more rapidly than those with parents from warmer (low altitude) populations. Growth rates also differed between offspring from high and low altitude populations. However, contrary to our prediction that crickets from high altitudes would grow faster, the most striking difference was that at high temperatures, growth was fastest in individuals from low altitudes. Our findings reveal that populations a few tens of kilometres apart have independently evolved adaptations to their environment. This suggests that local adaptation in a range of traits may be commonplace even in mobile invertebrates at scales of a small fraction of species' distributions.

Chill coma onset and recovery fail to reveal true variation in thermal performance among populations of Drosophila melanogaster

Species from colder climates tend to be more chill tolerant regardless of the chill tolerance trait measured, but for Drosophila melanogaster, population-level differences in chill tolerance among populations are not always found when a single trait is measured in the laboratory. We measured chill coma onset temperature, chill coma recovery time, and survival after chronic cold exposure in replicate lines derived from multiple paired African and European D. melanogaster populations. The populations in our study were previously found to differ in chronic cold survival ability, which is believed to have evolved independently in each population pair; however, they did not differ in chill coma onset temperature and chill coma recovery time in a manner that reflected their geographic origins, even though these traits are known to vary with origin latitude among Drosophila species and are among the most common metrics of thermal tolerance in insects. While it is common practice to measure only one chill tolerance trait when comparing chill tolerance among insect populations, our results emphasise the importance of measuring more than one thermal tolerance trait to minimize the risk of missing real adaptive variation in insect thermal tolerance.

Laboratory diet influences cold tolerance in a genotype-dependent manner in Drosophila melanogaster

Cold stress can reduce insect fitness and is an important determinant of species distributions and responses to climate change. Cold tolerance is influenced by genotype and environmental conditions, with factors such as day length and temperature having a particularly strong influence. Recent studies also indicate that diet impacts cold tolerance, but it is unclear whether diet-mediated shifts in cold tolerance are consistent across distinct genotypes. The goal of this study was to determine the extent to which commonly used artificial diets influence cold tolerance in Drosophila melanogaster, and whether these effects are consistent across genetically distinct lines. Specifically, we tested the impact of different fly diets on 1) ability to survive cold stress, 2) critical thermal minimum (CT_min), and 3) the ability to maintain reproduction after cold stress. Experiments were conducted across six isogenic lines from the Drosophila Genetic Reference Panel, and these lines were reared on different fly diets. Cold shock survival, CT_min, and reproductive output pre- and post-cold exposure varied considerably across diet and genotype combinations, suggesting strong genotype by environment interactions shape nutritionally mediated changes in cold tolerance. For example, in some lines cold shock survival remained consistently high or low across diets, while in others cold shock survival ranged from 5% to 75% depending on diet. Ultimately, these results add to a growing literature that cold tolerance is shaped by complex interactions between genotype and environment and inform practical considerations when selecting a laboratory diet for thermal tolerance experiments in Drosophila.

Evolutionary impacts of winter climate change on insects

Overwintering is a serious challenge for insects, and winters are rapidly changing as climate shifts. The capacity for phenotypic plasticity and evolutionary adaptation will determine which species profit or suffer from these changes. Here we discuss current knowledge on the potential and evidence for evolution in winter-relevant traits among insect species and populations. We conclude that the best evidence for evolutionary shifts in response to changing winters remain those related to changes in phenology, but all evidence points to cold hardiness as also having the potential to evolve in response to climate change. Predicting future population sizes and ranges relies on understanding to what extent evolution in winter-related traits is possible, and remains a serious challenge.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

Distinct cold tolerance traits independently vary across genotypes in Drosophila melanogaster

Cold tolerance, the ability to cope with low temperature stress, is a critical adaptation in thermally variable environments. An individual's cold tolerance comprises several traits including minimum temperatures for growth and activity, ability to survive severe cold, and ability to resume normal function after cold subsides. Across species, these traits are correlated, suggesting they were shaped by shared evolutionary processes or possibly share physiological mechanisms. However, the extent to which cold tolerance traits and their associated mechanisms covary within populations has not been assessed. We measured five cold tolerance traits-critical thermal minimum, chill coma recovery, short- and long-term cold tolerance, and cold-induced changes in locomotor behavior-along with cold-induced expression of two genes with possible roles in cold tolerance (heat shock protein 70 and frost)-across 12 lines of Drosophila melanogaster derived from a single population. We observed significant genetic variation in all traits, but few were correlated across genotypes, and these correlations were sex-specific. Further, cold-induced gene expression varied by genotype, but there was no evidence supporting our hypothesis that cold-hardy lines would have either higher baseline expression or induction of stress genes. These results suggest cold tolerance traits possess unique mechanisms and have the capacity to evolve independently.

Temperature Tolerance and Thermal Environment of European Seed Bugs

Heteroptera, or true bugs populate many climate zones, coping with different environmental conditions. The aim of this study was the evaluation of their thermal limits and derived traits, as well as climatological parameters which might influence their distribution. We assessed the thermal limits (critical thermal maxima, CT_max, and minima, CT_min) of eight seed bug species (Lygaeidae, Pyrrhocoridae) distributed over four Köppen–Geiger climate classification types (KCC), approximately 6° of latitude, and four European countries (Austria, Italy, Croatia, Bulgaria). In test tubes, a temperature ramp was driven down to −5 °C for CT_min and up to 50 °C for CT_max (0.25 °C/min) until the bugs’ voluntary, coordinated movement stopped. In contrast to CT_min, CT_max depended significantly on KCC, species, and body mass. CT_max showed high correlation with bioclimatic parameters such as annual mean temperature and mean maximum temperature of warmest month (BIO5), as well as three parameters representing temperature variability. CT_min correlated with mean annual temperature, mean minimum temperature of coldest month (BIO6), and two parameters representing variability. Although the derived trait cold tolerance (TC = BIO6 − CT_min) depended on several bioclimatic variables, heat tolerance (TH = CT_max − BIO5) showed no correlation. Seed bugs seem to have potential for further range shifts in the face of global warming.

Rapid cold hardening: ecological relevance, physiological mechanisms and new perspectives

Rapid cold hardening (RCH) is a type of phenotypic plasticity that allows ectotherms to quickly enhance cold tolerance in response to brief chilling (lasting minutes to hours). In this Review, we summarize the current state of knowledge of this important phenotype and provide new directions for research. As one of the fastest adaptive responses to temperature known, RCH allows ectotherms to cope with sudden cold snaps and to optimize their performance during diurnal cooling cycles. RCH and similar phenotypes have been observed across a diversity of ectotherms, including crustaceans, terrestrial arthropods, amphibians, reptiles, and fish. In addition to its well-defined role in enhancing survival to extreme cold, RCH also protects against nonlethal cold injury by preserving essential functions following cold stress, such as locomotion, reproduction, and energy balance. The capacity for RCH varies across species and across genotypes of the same species, indicating that RCH can be shaped by selection and is likely favored in thermally variable environments. Mechanistically, RCH is distinct from other rapid stress responses in that it typically does not involve synthesis of new gene products; rather, the existing cellular machinery regulates RCH through post-translational signaling mechanisms. However, the protective mechanisms that enhance cold hardiness are largely unknown. We provide evidence that RCH can be induced by multiple triggers in addition to low temperature, and that rapidly induced tolerance and cross-tolerance to a variety of environmental stressors may be a general feature of stress responses that requires further investigation.

Evidence for a rapid cold hardening response in cultured Drosophila S2 cells

The ability to quickly respond to changes in environmental temperature is critical for organisms living in thermally variable environments. To cope with sudden drops in temperature, insects and other ectotherms are capable of rapid cold hardening (RCH), in which mild chilling significantly enhances cold tolerance within minutes. While the ecological significance of RCH is well established, the mechanisms underlying RCH are still poorly understood. Previous work has demonstrated that RCH is regulated at the cellular level by post-translational signaling mechanisms, and here we tested the hypothesis that cultured cells are capable of RCH. A 2 h cold shock at -8°C significantly reduced the metabolic viability of Drosophila S2 cells, but pre-treatment with RCH at 4°C for 2 h prevented this decrease in viability. Thus, S2 cells are capable of RCH in a similar manner to whole insects and provide a new system for investigating the cell biology of RCH.

Testing the thermal limits: Non-linear reaction norms drive disparate thermal acclimation responses in Drosophila melanogaster

Critical thermal limits are important ecological parameters for studying thermal biology and for modelling species' distributions under current and changing climatic conditions (including predicting the risk of extinction for species from future warming). However, estimates of the critical thermal limits are biased by the choice of assay and assay conditions, which differ among studies. Furthermore, estimates of the potential for phenotypic plasticity (thermal acclimation) to buffer against thermal variability are usually based on single assay conditions and (usually linear) extrapolation from a few acclimation temperatures. We produced high resolution estimates of adult acclimation capacity for upper tolerance limits at different assay conditions (ramping rates and knock-down temperatures) using CTmax (dynamic) and knock-down (static) thermal assays in the model species Drosophila melanogaster. We found the reaction norms to be highly dependent on assay conditions. We confirmed that progressively lower ramping rates or higher knock-down temperatures led to overall lower tolerance estimates. More surprisingly, extended assays (lower ramping rates or lower knock-down temperatures) also led to increasingly non-linear reaction norms for upper thermal tolerance across adult acclimation temperatures. Our results suggest that the magnitude (capacity) and direction (beneficial or detrimental) of acclimation responses are highly sensitive to assay conditions. The results offer a framework for comparison of acclimation responses between different assay conditions and a potential for explaining disparate acclimation capacity theories. We advocate cautious interpretation of acclimation capacities and careful consideration of assay conditions, which should represent realistic environmental conditions based on species' ecological niches.