Quantitative trait loci for thermotolerance phenotypes in Drosophila melanogaster

Transcriptomic analysis of the response of Spodoptera frugiperda (Lepidoptera: Noctuidae) to short-term low-temperature stress

Spodoptera frugiperda is a major invasive pest that poses a serious threat to crops worldwide. Low temperature is a key factor limiting the survival and reproduction for this pest. To study the responses of S. frugiperda to low-temperature stress, high-throughput sequencing was used to perform transcriptomic analysis on the 6th instar larvae under low-temperature stress at 5 °C and 10 °C, along with 25 °C as a control. As a result, 215 differentially expressed genes (DEGs) were identified under different low-temperature stresses. Upon functional annotation of the DEGs in KEGG and GO databases, the number of DEGs annotated in control vs. LT10 comparison was the largest (n = 150), whereas fewer DEGs (n = 89) were annotated in control vs. LT5 comparison. This discrepancy suggested that S. frugiperda might adopt different strategies to cope with low-temperature stress. The DEGs in the GO database were particularly associated with cell catalytic activity, cell anatomical entity process, cell apoptosis, and cell binding channel. KEGG annotation analysis of the different low-temperature stresses showed that most of the enriched pathways were related to carbon metabolism, oxidative phosphorylation, and lipid metabolism. The results will be the basis for mastering the cold tolerant mechanism of S. frugiperda, and is of great significance for its prevention.

Genetic Architecture of the Thermal Tolerance Landscape in Drosophila melanogaster

Increased environmental temperatures associated with global warming strongly impact natural populations of ectothermic species. Therefore, it is crucial to understand the genetic basis and evolutionary potential of heat tolerance. However, heat tolerance and its genetic components depend on the methodology, making it difficult to predict the adaptive responses to global warming. Here, we measured the knockdown time for 100 lines from the Drosophila Genetic Reference Panel (DGRP) at four different static temperatures, and we estimated their thermal-death-time (TDT) curves, which incorporate the magnitude and the time of exposure to thermal stress, to determine the genetic basis of the thermal tolerance landscape. Through quantitative genetic analyses, the knockdown time showed a significant heritability at different temperatures and that its genetic correlations decreased as temperatures differences increased. Significant genotype-by-sex and genotype-by-environment interactions were noted for heat tolerance. We also discovered genetic variability for the two parameters of TDT: CTmax and thermal sensitivity. Taking advantage of the DGRP, we performed a GWAS and identified multiple variants associated with the TDT parameters, which mapped to genes related to signalling and developmental functions. We performed functional validations for some candidate genes using RNAi, which revealed that genes such as mam, KNCQ, or robo3 affect the knockdown time at a specific temperature but are not associated with the TDT parameters. In conlusion, the thermal tolerance landscape display genetic variation and plastic responses, which may facilitate the adaptation of Drosophila populations to a changing world.

Geographic variation in larval cold tolerance and exposure across the invasion front of a widely established forest insect

Under global climate change, high and low temperature extremes can drive shifts in species distributions. Across the range of a species, thermal tolerance is based on acclimatization, plasticity, and may undergo selection, shaping resilience to temperature stress. In this study, we measured variation in cold temperature tolerance of early instar larvae of an invasive forest insect, Lymantria dispar dispar L. (Lepidoptera: Erebidae), using populations sourced from a range of climates within the current introduced range in the Eastern United States. We tested for population differences in chill coma recovery (CCR) by measuring recovery time following a period of exposure to a nonlethal cold temperature in 2 cold exposure experiments. A 3rd experiment quantified growth responses after CCR to evaluate sublethal effects. Our results indicate that cold tolerance is linked to regional climate, with individuals from populations sourced from colder climates recovering faster from chill coma. While this geographic gradient is seen in many species, detecting this pattern is notable for an introduced species founded from a single point-source introduction. We demonstrate that the cold temperatures used in our experiments occur in nature during cold spells after spring egg hatch, but impacts to growth and survival appear low. We expect that population differences in cold temperature performance manifest more from differences in temperature-dependent growth than acute exposure. Evaluating intraspecific variation in cold tolerance increases our understanding of the role of climatic gradients on the physiology of an invasive species, and contributes to tools for predicting further expansion.

Naturally segregating genetic variants contribute to thermal tolerance in a Drosophila melanogaster model system

Thermal tolerance is a fundamental physiological complex trait for survival in many species. For example, everyday tasks such as foraging, finding a mate, and avoiding predation are highly dependent on how well an organism can tolerate extreme temperatures. Understanding the general architecture of the natural variants within the genes that control this trait is of high importance if we want to better comprehend thermal physiology. Here, we take a multipronged approach to further dissect the genetic architecture that controls thermal tolerance in natural populations using the Drosophila Synthetic Population Resource as a model system. First, we used quantitative genetics and Quantitative Trait Loci mapping to identify major effect regions within the genome that influences thermal tolerance, then integrated RNA-sequencing to identify differences in gene expression, and lastly, we used the RNAi system to (1) alter tissue-specific gene expression and (2) functionally validate our findings. This powerful integration of approaches not only allows for the identification of the genetic basis of thermal tolerance but also the physiology of thermal tolerance in a natural population, which ultimately elucidates thermal tolerance through a fitness-associated lens.

Naturally segregating genetic variants contribute to thermal tolerance in a D. melanogaster model system

Thermal tolerance is a fundamental physiological complex trait for survival in many species. For example, everyday tasks such as foraging, finding a mate, and avoiding predation, are highly dependent on how well an organism can tolerate extreme temperatures. Understanding the general architecture of the natural variants of the genes that control this trait is of high importance if we want to better comprehend how this trait evolves in natural populations. Here, we take a multipronged approach to further dissect the genetic architecture that controls thermal tolerance in natural populations using the Drosophila Synthetic Population Resource (DSPR) as a model system. First, we used quantitative genetics and Quantitative Trait Loci (QTL) mapping to identify major effect regions within the genome that influences thermal tolerance, then integrated RNA-sequencing to identify differences in gene expression, and lastly, we used the RNAi system to 1) alter tissue-specific gene expression and 2) functionally validate our findings. This powerful integration of approaches not only allows for the identification of the genetic basis of thermal tolerance but also the physiology of thermal tolerance in a natural population, which ultimately elucidates thermal tolerance through a fitness-associated lens.

Pleiotropic fitness effects of the lncRNA Uhg4 in Drosophila melanogaster

BACKGROUND

Long noncoding RNAs (lncRNAs) are a diverse class of RNAs that are critical for gene regulation, DNA repair, and splicing, and have been implicated in development, stress response, and cancer. However, the functions of many lncRNAs remain unknown. In Drosophila melanogaster, U snoRNA host gene 4 (Uhg4) encodes an antisense long noncoding RNA that is host to seven small nucleolar RNAs (snoRNAs). Uhg4 is expressed ubiquitously during development and in all adult tissues, with maximal expression in ovaries; however, it has no annotated function(s).

RESULTS

We used CRISPR-Cas9 germline gene editing to generate multiple deletions spanning the promoter region and first exon of Uhg4. Females showed arrested egg development and both males and females were sterile. In addition, Uhg4 deletion mutants showed delayed development and decreased viability, and changes in sleep and responses to stress. Whole-genome RNA sequencing of Uhg4 deletion flies and their controls identified co-regulated genes and genetic interaction networks associated with Uhg4. Gene ontology analyses highlighted a broad spectrum of biological processes, including regulation of transcription and translation, morphogenesis, and stress response.

CONCLUSION

Uhg4 is a lncRNA essential for reproduction with pleiotropic effects on multiple fitness traits.

Urban evolution of thermal physiology in a range-expanding, mycophagous fruit fly, Drosophila tripunctata

In Drosophila spp., their often high number of annual generations, large population sizes and large amounts of standing genetic variation should predispose them to undergo contemporary adaptation to climatic warming. However, a number of laboratory experimental evolution studies in this group of organisms suggest strong limits on the rate and magnitude of contemporary thermal adaptation. Here, we explore this discrepancy by examining the potential for rapid evolutionary divergence between wild populations of Drosophila tripunctata Loew, 1862 from rural and urban sites. We performed a multi-generation common garden study and found evidence for the evolution of higher heat tolerance (critical thermal maximum) in flies from urban populations. We also detected evolutionary divergence in cold resistance (chill coma recovery time), with diminished cold resistance in flies from urban populations, although the effect was weaker than the shift in heat tolerance. Our study provides evidence of contemporary urban thermal adaptation, although the magnitude of phenotypic change lagged the magnitude of environmental temperature change across the urbanization gradient, suggesting potential limits on the evolution of urban thermal physiology.

Sexual dimorphism in phenotypic plasticity and persistence under environmental change: An extension of theory and meta-analysis of current data

Populations must adapt to environmental changes to remain viable. Both evolution and phenotypic plasticity contribute to adaptation, with plasticity possibly being more important for coping with rapid change. Adaptation is complex in species with separate sexes, as the sexes can differ in the strength or direction of natural selection, the genetic basis of trait variation, and phenotypic plasticity. Many species show sex differences in plasticity, yet how these differences influence extinction susceptibility remains unclear. We first extend theoretical models of population persistence in changing environments and show that persistence is affected by sexual dimorphism for phenotypic plasticity, trait genetic architecture, and sex-specific selection. Our models predict that female-biased adaptive plasticity-particularly in traits with modest-to-low cross-sex genetic correlations-typically promotes persistence, though we also identify conditions where sexually monomorphic or male-biased plasticity promotes persistence. We then perform a meta-analysis of sex-specific plasticity under manipulated thermal conditions. Although examples of sexually dimorphic plasticity are widely observed, systematic sex differences are rare. An exception-cold resistance-is systematically female-biased and represents a trait wherein sexually dimorphic plasticity might elevate population viability in changing environments. We discuss our results in light of debates about the roles of evolution and plasticity in extinction susceptibility.

Differential Transcriptome Analysis Reveals Genes Related to Low- and High-Temperature Stress in the Fall Armyworm, Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is regarded as one of the world’s most harmful plant pests. This research examines the molecular response processes of FAW to low temperature (4°C) and high temperature (40°C) when gene expression is compared to controls (25°C). A total of 211,967 unigenes were collected, at least 14,338 of which were annotated with gene descriptions, gene ontology terms, and metabolic pathways. There were 50 Gene Ontology (GO) functional sub-groups and 21 EggNOG words as a result. Differentially expresses genes (DEGs) with log₂FC ≥ 2 were identified and compared at various temperatures. In comparison to the 25°C treated group, we discovered 199 and 1,248 individual DEGs co-regulated at 4 and 40°C, respectively. Comparing transcriptome profiles for differential gene expression revealed a number of DEGs, including cytochrome P450, odorant binding proteins (OBPs), and immune system genes previously implicated in cold and high temperature stresses. The enrichment pathways were identified using Kyoto Encyclopedia of Genes and Genomics (KEGG) analysis, and heatmaps of similar unigenes from both treatment groups (T4 and T40) were plotted. We used quantitative reverse transcription PCR (RT-qPCR) to confirm the RNA-seq data on 10 up- and down-regulated DEGs. These findings provide a foundation for future understanding of FAW adaptation mechanisms and the underlying basis underlying the response to low and high temperatures.

Genetic variation in the heat-stress survival of embryos is largely decoupled from adult thermotolerance in an intercontinental set of recombinant lines of Drosophila melanogaster

In insects, thermal adaptation works on the genetic variation for thermotolerance of not only larvae and adults but also of the immobile stages of the life cycle including eggs. In contrast to adults and larvae, the genetic basis for thermal adaptation in embryos (eggs) remains to be tested in the model insect Drosophila melanogaster. Quantitative-trait loci (QTL) for heat-stress resistance in embryos could largely differ from previously identified QTL for larvae and adults. Here we used an intercontinental set of recombinant inbred lines (RIL), which were previously used to identify thermotolerance-QTLs in adults and larvae because of their high variation segregating for adult thermotolerance. Eggs appeared to be more heat resistant than larvae and adults from previous studies on these RIL, though different heat-shock assays were used in previous studies. We found that variation in thermotolerance in embryos can be, at least partially, genetically decoupled from thermotolerance in the adult insect. Some RIL that are heat resistant in the adult and larvae can be heat susceptible in embryos. Only one small-effect QTL out of five autosomal QTL co-localized between embryo and other ontogenetic stages. These results suggest that selection for thermal adaptation in adult flies and larvae is predicted to have only a small impact on embryo thermotolerance. In addition, heat-stress tolerance of insects can be measured across ontogenetic stages including embryos in order to better predict thermal adaptive limits of populations and species.

Cardiac performance in heat-stressed flies of heat-susceptible and heat-resistant Drosophila melanogaster

Thermotolerance is a complex trait that can greatly differ between heat-susceptible (HS) and heat-adapted populations of small insects including Drosophila, with short-term effects after a sub-lethal level of heat stress on many physiological functions. Cardiac performance could accordingly be more robust in heat-resistant (HR) than in HS individuals under heat stress. Here, we tested heart performance under heat-stress effects in two recombinant inbred lines (RIL) of Drosophila melanogaster that dramatically differ in heat knockdown resistance. Heart rate did not strongly differ between heat-susceptible and heat-tolerant flies after a sub-lethal heat stress. Instead, heat-susceptible flies showed a much higher arrhythmia incidence, a longer duration of each heartbeat, and a larger amount of bradycardia than heat-tolerant flies. The highly conserved cardiac proteins SERCA, RyR and NCX that participate in the excitation/contraction coupling, did not differ in activity level between HR and HS flies. Available information for both RIL suggests that heart performance under heat stress may be linked, at least partially, to candidate genes of previously identified quantitative trait loci (QTL) for thermotolerance. This study indicates that HR flies can be genetically more robust in their heart performance than HS flies under even sub-lethal levels of heat stress.

The genetic architecture of temperature adaptation is shaped by population ancestry and not by selection regime

BACKGROUND

Understanding the genetic architecture of temperature adaptation is key for characterizing and predicting the effect of climate change on natural populations. One particularly promising approach is Evolve and Resequence, which combines advantages of experimental evolution such as time series, replicate populations, and controlled environmental conditions, with whole genome sequencing. Recent analysis of replicate populations from two different Drosophila simulans founder populations, which were adapting to the same novel hot environment, uncovered very different architectures-either many selection targets with large heterogeneity among replicates or fewer selection targets with a consistent response among replicates.

RESULTS

Here, we expose the founder population from Portugal to a cold temperature regime. Although almost no selection targets are shared between the hot and cold selection regime, the adaptive architecture was similar. We identify a moderate number of targets under strong selection (19 selection targets, mean selection coefficient = 0.072) and parallel responses in the cold evolved replicates. This similarity across different environments indicates that the adaptive architecture depends more on the ancestry of the founder population than the specific selection regime.

CONCLUSIONS

These observations will have broad implications for the correct interpretation of the genomic responses to a changing climate in natural populations.

Adaptation of Drosophila subobscura chromosomal inversions to climatic variables: the Balkan natural population of Avala

The adaptive value of chromosomal inversions continues raising relevant questions in evolutionary biology. In many species of the Drosophila genus, different inversions have been recognized to be related to thermal adaptation, but it is necessary to determine to which specific climatic variables the inversions are adaptive. With this aim, the behavior of thermal adapted inversions of Drosophila subobscura regarding climatic variables was studied in the natural population of Avala (Serbia) during the 2014-2017 period. The results obtained were compared with those previously reported in the Font Groga (Barcelona, Spain) population, which presents different climatic and environmental conditions. In both populations, it was observed that most thermal adapted inversions were significantly associated with the first, second or both principal components, which were related with maximum, minimum and mean temperatures. Moreover, a significant increase over years (2004-2017) for the minimum temperature was detected. In parallel, a significant variation over time in Avala was only observed for the frequencies of 'warm' and 'non-thermal' adapted inversions of the U chromosome. However, stability in the chromosomal inversion polymorphism was observed for the 2014-2017 period which might result from the temporal span of the study and/or selective process acting on the population.

False and true positives in arthropod thermal adaptation candidate gene lists

Genome-wide studies are prone to false positives due to inherently low priors and statistical power. One approach to ameliorate this problem is to seek validation of reported candidate genes across independent studies: genes with repeatedly discovered effects are less likely to be false positives. Inversely, genes reported only as many times as expected by chance alone, while possibly representing novel discoveries, are also more likely to be false positives. We show that, across over 30 genome-wide studies that reported Drosophila and Daphnia genes with possible roles in thermal adaptation, the combined lists of candidate genes and orthologous groups are rapidly approaching the total number of genes and orthologous groups in the respective genomes. This is consistent with the expectation of high frequency of false positives. The majority of these spurious candidates have been identified by one or a few studies, as expected by chance alone. In contrast, a noticeable minority of genes have been identified by numerous studies with the probabilities of such discoveries occurring by chance alone being exceedingly small. For this subset of genes, different studies are in agreement with each other despite differences in the ecological settings, genomic tools and methodology, and reporting thresholds. We provide a reference set of presumed true positives among Drosophila candidate genes and orthologous groups involved in response to changes in temperature, suitable for cross-validation purposes. Despite this approach being prone to false negatives, this list of presumed true positives includes several hundred genes, consistent with the "omnigenic" concept of genetic architecture of complex traits.

Toxicological Study and Genetic Basis of BTEX Susceptibility in Drosophila melanogaster

Benzene, toluene, ethylbenzene and xylene, also known as BTEX, are released into environmental media by petroleum product exploratory and exploitative activities and are harmful to humans and animals. Testing the effects of these chemicals on a significantly large scale requires an inexpensive, rapidly developing model organism such as Drosophila melanogaster. In this study, the toxicological profile of benzene, toluene, ethylbenzene, p-xylene, m-xylene, and o-xylene in D. melanogaster was evaluated. Adult animals were monitored for acute toxicity effects. Similarly, first instar larvae reared separately on the same compounds were monitored for the ability to develop into adult flies (eclosion). Further, the impact of fixed concentrations of benzene and xylene on apoptosis and mitosis were investigated in adult progenitor tissues found in third instar larvae. Toluene is the most toxic to adult flies with an LC₅₀ of 0.166 mM, while a significant and dose-dependent decrease in fly eclosion was observed with benzene, p-xylene, and o-xylene. An increase in apoptosis and mitosis was also observed in animals exposed to benzene and p-xylene. Through Genome Wide Association Screening (GWAS), 38 regions of the D. melanogaster genome were identified as critical for responses to p-xylene. This study reveals the strength of D. Melanogaster genetics as an accessible approach to study BTEX compounds.

Measuring the evolutionary potential of a winter-active parasitic wasp to climate change

In temperate climates, as a consequence of warming winters, an increasing number of ectothermic species are remaining active throughout winter months instead of diapausing, rendering them increasingly vulnerable to unpredictable cold events. One species displaying a shift in overwintering strategy is the parasitoid wasp and biological control agent Aphidius avenae. The current study aimed to better understand the consequence of a changing overwintering strategy on the evolutionary potential of an insect population to adapt to the cold stress events, set to increase in frequency, even during milder winters. Using a parental half-sibling breeding design, narrow-sense heritability of the cold tolerance, morphology and longevity of A. avenae was estimated. The heritability of cold tolerance was estimated at 0.07 (CI_95% = [0.00; 0.25]) for the Critical Thermal Minima (CT_min) and 0.11 (CI_95% = [0.00; 0.34]) for chill coma temperature; estimates much lower than those obtained for morphological traits (tibia length 0.20 (CI_95% = [0.03; 0.37]); head width 0.23 (CI_95% = [0.09; 0.39]); wing surface area 0.28 (CI_95% = [0.11; 0.47])), although comparable with the heritability estimate of 0.12 obtained for longevity (CI_95% = [0.00; 0.25]). The heritability estimates obtained thus suggest that A. avenae possesses low adaptive potential against cold stress. If such estimates are indicative of the evolutionary potential of A. avenae cold tolerance, more emphasis may be placed on adaptive phenotypic plasticity at the individual level to persist in a changing climate, with potential implications for the biological control function they provide.

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.

Heat knockdown resistance and chill-coma recovery as correlated responses to selection on mating success at high temperature in Drosophila buzzatii

Reproduction and related traits such as mating success are strongly affected by thermal stress. We tested direct and correlated responses to artificial selection in replicated lines of Drosophila buzzatii that were selected for mating success at high temperature. Knockdown resistance at high temperature (KRHT) and chill-coma recovery (CCR) were tested as correlated selection responses. Virgin flies were allowed to mate for four hours at 33°C in three replicated lines (S lines) to obtain the selected flies and then returned at 25°C to lay eggs. Other three replicated lines were maintained at 25°C without any selection as control (C lines). After 15 selection generations, KRHT and CCR were measured. Both traits were assessed in flies that did not receive any hardening pretreatments as well as in flies that were either heat or cold hardened. Thermotolerance traits showed significant correlated responses with higher KRHT in S than in C lines, both with a heat-hardening pretreatment and without a heat-hardening pretreatment. CCR time was longer in S than in C lines both with a cold-hardening pretreatment and without a cold-hardening pretreatment. Hardening treatments improved both KRHT and CCR in all cases excepting KRHT in C lines. Overall, KRHT and CCR showed an antagonistic pattern of correlated responses to our selection regime, suggesting either pleiotropy or tightly linked trait-specific genes partially affecting KRHT and CCR.

Three Quantitative Trait Loci Explain More than 60% of Variation for Chill Coma Recovery Time in a Natural Population of Drosophila ananassae

Ectothermic species such as insects are particularly vulnerable to climatic fluctuations. Nevertheless, many insects that evolved and diversified in the tropics have successfully colonized temperate regions all over the globe. To shed light on the genetic basis of cold tolerance in such species, we conducted a quantitative trait locus (QTL) mapping experiment for chill coma recovery time (CCRT) in Drosophila ananassae, a cosmopolitan species that has expanded its range from tropical to temperate regions. We created a mapping population of recombinant inbred advanced intercross lines (RIAILs) from two founder strains with diverging CCRT phenotypes. The RIAILs were phenotyped for their CCRT and, together with the founder strains, genotyped for polymorphic markers with double-digest restriction site-associated DNA (ddRAD) sequencing. Using a hierarchical mapping approach that combined standard interval mapping and a multiple-QTL model, we mapped three QTL which altogether explained 64% of the phenotypic variance. For two of the identified QTL, we found evidence of epistasis. To narrow down the list of cold tolerance candidate genes, we cross-referenced the QTL intervals with genes that we previously identified as differentially expressed in response to cold in D. ananassae, and with thermotolerance candidate genes of D. melanogaster. Among the 58 differentially expressed genes that were contained within the QTL, GF15058 showed a significant interaction of the CCRT phenotype and gene expression. Further, we identified the orthologs of four D. melanogaster thermotolerance candidate genes, MtnA, klarsicht, CG5246 (D.ana/GF17132) and CG10383 (D.ana/GF14829) as candidates for cold tolerance in D. ananassae.

Stage-specific genotype-by-environment interactions for cold and heat hardiness in Drosophila melanogaster

Environments often vary across a life cycle, imposing fluctuating natural selection across development. Such fluctuating selection can favor different phenotypes in different life stages, but stage-specific evolutionary responses will depend on genetic variance, covariance, and their interaction across development and across environments. Thus, quantifying how genetic architecture varies with plastic responses to the environment and across development is vital to predict whether stage-specific adaptation will occur in nature. Additionally, the interaction of genetic variation and environmental plasticity (GxE) may be stage-specific, leading to a three-way interaction between genotype, environment, and development or GxDxE. To test for these patterns, we exposed larvae and adults of Drosophila melanogaster isogenic lines derived from a natural population to extreme heat and cold stress after developmental acclimation to cool (18 °C) and warm (25 °C) conditions and measured genetic variance for thermal hardiness. We detected significant GxE that was specific to larvae and adults for cold and heat hardiness (GxDxE), but no significant genetic correlation across development for either trait at either acclimation temperature. However, cross-development phenotypic correlations for acclimation responses suggest that plasticity itself may be developmentally constrained, though rigorously testing this hypothesis requires more experimentation. These results illustrate the potential for stage-specific adaptation within a complex life cycle and demonstrate the importance of measuring traits at appropriate developmental stages and environmental conditions when predicting evolutionary responses to changing climates.

Genetic Decoupling of Thermal Hardiness across Metamorphosis in Drosophila melanogaster

As organisms age the environment fluctuates, exerting differential selection across ontogeny. In particular, highly seasonal environments expose life stages to often drastically different thermal environments. This developmental variation is particularly striking in organisms with complex life cycles, wherein life history stages also exhibit distinct morphologies, physiologies, and behaviors. Genes acting pleiotropically on thermal responses may produce genetic correlations across ontogeny, constraining the independent evolution of each life stage to their respective thermal environments. To investigate whether developmental genetic correlations constrain the evolution thermal hardiness of the fly Drosophila melanogaster, we applied quantitative genetic analyses to cold hardiness measured in both larvae and adults from isogenic lines of the Drosophila Genetic Reference Panel (DGRP), using survival at stressful low temperatures as the phenotypic metric. Using full genome resequencing data for the DGRP, we also implemented genome-wide association (GWA) analysis using Bayesian Sparse Linear Mixed Models (BSLMMs) to estimate associations between naturally segregating variation and cold hardiness for both larvae and adults. Quantitative genetic analyses revealed no significant genetic correlation for cold hardiness between life stages, suggesting complete genetic decoupling of thermal hardiness across the metamorphic boundary. Both quantitative genetic and GWA analyses suggested that polygenic variation underlies cold hardiness in both stages, and that associated loci largely affected one stage or the other, but not both. However, reciprocal enrichment tests and correlations between BSLMM parameters for each life stage support some shared physiological mechanisms that may reflect common cellular thermal response pathways. Overall, these results suggest no developmental genetic constraints on cold hardiness across metamorphosis in D. melanogaster, an important consideration in evolutionary models of responses to changing climates. Genetic correlations for environmental sensitivity across ontogeny remains largely unexplored in other organisms, thus assessing the generality of genetic decoupling will require further quantitative or population genetic analysis in additional species.

Costs of cold acclimation on survival and reproductive behavior in Drosophila melanogaster

Fitness is determined by the ability of an organism to both survive and reproduce; however, the mechanisms that lead to increased survival may not have the same effect on reproductive success. We used nineteen natural Drosophila melanogaster genotypes from the Drosophila Genetic Reference Panel to determine if adaptive plasticity following short-term acclimation through rapid cold-hardening (RCH) affects mating behavior and mating success. We confirmed that exposure to the acclimation temperature is beneficial to survival following cold stress; however, we found that this same acclimation temperature exposure led to less efficient male courtship and a significant decrease in the likelihood of mating. Cold tolerance and the capacity to respond plastically to cold stress were not correlated with mating behavior following acclimation, suggesting that the genetic control of the physiological effects of the cold temperature exposure likely differ between survival and behavioral responses. We also tested whether the exposure of males to the acclimation temperature influenced courtship song. This exposure again significantly increased courtship duration; however, courtship song was unchanged. These results illustrate costs of short-term acclimation on survival and reproductive components of fitness and demonstrate the pronounced effect that short-term thermal environment shifts can have on reproductive success.

Environmental and genetic control of cold tolerance in the Glanville fritillary butterfly

Thermal tolerance has a major effect on individual fitness and species distributions and can be determined by genetic variation and phenotypic plasticity. We investigate the effects of developmental and adult thermal conditions on cold tolerance, measured as chill coma recovery (CCR) time, during the early and late adult stage in the Glanville fritillary butterfly. We also investigate the genetic basis of cold tolerance by associating CCR variation with polymorphisms in candidate genes that have a known role in insect physiology. Our results demonstrate that a cooler developmental temperature leads to reduced cold tolerance in the early adult stage, whereas cooler conditions during the adult stage lead to increased cold tolerance. This suggests that adult acclimation, but not developmental plasticity, of adult cold tolerance is adaptive. This could be explained by the ecological conditions the Glanville fritillary experiences in the field, where temperature during early summer, but not spring, is predictive of thermal conditions during the butterfly's flight season. In addition, an amino acid polymorphism (Ala-Glu) in the gene flightin, which has a known function in insect flight and locomotion, was associated with CCR. These amino acids have distinct biochemical properties and may thus affect protein function and/or structure. To our knowledge, our study is the first to link genetic variation in flightin to cold tolerance, or thermal adaptation in general.

Chromosomal Thermal Index: a comprehensive way to integrate the thermal adaptation of Drosophila subobscura whole karyotype

Drosophila has demonstrated to be an excellent model to study the adaptation of organisms to global warming, with inversion chromosomal polymorphism having a key role in this adaptation. Here, we introduce a new index (Chromosomal Thermal Index or CTI) to quantify the thermal adaptation of a population according to its composition of "warm" and "cold" adapted inversions. This index is intuitive, has good statistical properties, and can be used to hypothesis on the effect of global warming on natural populations. We show the usefulness of CTI using data from European populations of D. subobscura, sampled in different years. Out of 15 comparisons over time, nine showed significant increase of CTI, in accordance with global warming expectations. Although large regions of the genome outside inversions contain thermal adaptation genes, our results show that the total amount of warm or cold inversions in populations seems to be directly involved in thermal adaptation, whereas the interactions between the inversions content of homologous and non-homologous chromosomes are not relevant.

Seasonal variation in basal and plastic cold tolerance: Adaptation is influenced by both long- and short-term phenotypic plasticity

Understanding how thermal selection affects phenotypic distributions across different time scales will allow us to predict the effect of climate change on the fitness of ectotherms. We tested how seasonal temperature variation affects basal levels of cold tolerance and two types of phenotypic plasticity in Drosophila melanogaster. Developmental acclimation occurs as developmental stages of an organism are exposed to seasonal changes in temperature and its effect is irreversible, while reversible short-term acclimation occurs daily in response to diurnal changes in temperature. We collected wild flies from a temperate population across seasons and measured two cold tolerance metrics (chill-coma recovery and cold stress survival) and their responses to developmental and short-term acclimation. Chill-coma recovery responded to seasonal shifts in temperature, and phenotypic plasticity following both short-term and developmental acclimation improved cold tolerance. This improvement indicated that both types of plasticity are adaptive, and that plasticity can compensate for genetic variation in basal cold tolerance during warmer parts of the season when flies tend to be less cold tolerant. We also observed a significantly stronger trade-off between basal cold tolerance and short-term acclimation during warmer months. For the longer-term developmental acclimation, a trade-off persisted regardless of season. A relationship between the two types of plasticity may provide additional insight into why some measures of thermal tolerance are more sensitive to seasonal variation than others.

Multiple Trait Covariance Association Test Identifies Gene Ontology Categories Associated with Chill Coma Recovery Time in Drosophila melanogaster

The genomic best linear unbiased prediction (GBLUP) model has proven to be useful for prediction of complex traits as well as estimation of population genetic parameters. Improved inference and prediction accuracy of GBLUP may be achieved by identifying genomic regions enriched for causal genetic variants. We aimed at searching for patterns in GBLUP-derived single-marker statistics, by including them in genetic marker set tests, that could reveal associations between a set of genetic markers (genomic feature) and a complex trait. GBLUP-derived set tests proved to be powerful for detecting genomic features, here defined by gene ontology (GO) terms, enriched for causal variants affecting a quantitative trait in a population with low degree of relatedness. Different set test approaches were compared using simulated data illustrating the impact of trait- and genomic feature-specific factors on detection power. We extended the most powerful single trait set test, covariance association test (CVAT), to a multiple trait setting. The multiple trait CVAT (MT-CVAT) identified functionally relevant GO categories associated with the quantitative trait, chill coma recovery time, in the unrelated, sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel.

CRISPR-induced null alleles show that Frost protects Drosophila melanogaster reproduction after cold exposure

The ability to survive and reproduce after cold exposure is important in all kingdoms of life. However, even in a sophisticated genetic model system like Drosophila melanogaster, few genes have been identified as functioning in cold tolerance. The accumulation of the Frost (Fst) gene transcript increases after cold exposure, making it a good candidate for a gene that has a role in cold tolerance. However, despite extensive RNAi knockdown analysis, no role in cold tolerance has been assigned to Fst. CRISPR is an effective technique for completely knocking down genes, and less likely to produce off-target effects than GAL4-UAS RNAi systems. We have used CRISPR-mediated homologous recombination to generate Fst null alleles, and these Fst alleles uncovered a requirement for FST protein in maintaining female fecundity following cold exposure. However, FST does not have a direct role in survival following cold exposure. FST mRNA accumulates in the Malpighian tubules, and the FST protein is a highly disordered protein with a putative signal peptide for export from the cell. Future work is needed to determine whether FST is exported from the Malpighian tubules and directly interacts with female reproductive tissues post-cold exposure, or if it is required for other repair/recovery functions that indirectly alter energy allocation to reproduction.

Quantitative Trait Loci and Antagonistic Associations for Two Developmentally Related Traits in the Drosophila Head

In insects, some developmentally related traits are negatively correlated. Here, we mapped Quantitative Trait Loci (QTL) for traits of eye size and head capsule, in an intercontinental set of recombinant inbred lines (RILs) of Drosophila melanogaster Composite interval mapping identified QTL on all major chromosomes. Two negatively correlated traits (size of eyes and between-eyes distance) were influenced by one QTL that appeared to be antagonistic between the traits (QTL cytological range is 25F5-30A6), consistent with a negative genetic correlation between these traits of the head capsule. Comparisons of QTL across traits indicated a nonrandom distribution over the genome, with a considerable overlap between some QTL across traits. Developmentally-related traits were influenced by QTL in a pattern that is consistent both with 1) the sign of the genetic correlation between the traits and 2) a constraint in the micro-evolutionary differentiation in the traits.

Multiple Independent Retroelement Insertions in the Promoter of a Stress Response Gene Have Variable Molecular and Functional Effects in Drosophila

Promoters are structurally and functionally diverse gene regulatory regions. The presence or absence of sequence motifs and the spacing between the motifs defines the properties of promoters. Recent alternative promoter usage analyses in Drosophila melanogaster revealed that transposable elements significantly contribute to promote diversity. In this work, we analyzed in detail one of the transposable element insertions, named FBti0019985, that has been co-opted to drive expression of CG18446, a candidate stress response gene. We analyzed strains from different natural populations and we found that besides FBti0019985, there are another eight independent transposable elements inserted in the proximal promoter region of CG18446. All nine insertions are solo-LTRs that belong to the roo family. We analyzed the sequence of the nine roo insertions and we investigated whether the different insertions were functionally equivalent by performing 5'-RACE, gene expression, and cold-stress survival experiments. We found that different insertions have different molecular and functional consequences. The exact position where the transposable elements are inserted matters, as they all showed highly conserved sequences but only two of the analyzed insertions provided alternative transcription start sites, and only the FBti0019985 insertion consistently affects CG18446 expression. The phenotypic consequences of the different insertions also vary: only FBti0019985 was associated with cold-stress tolerance. Interestingly, the only previous report of transposable elements inserting repeatedly and independently in a promoter region in D. melanogaster, were also located upstream of a stress response gene. Our results suggest that functional validation of individual structural variants is needed to resolve the complexity of insertion clusters.

Adaptation to Low Temperature Exposure Increases Metabolic Rates Independently of Growth Rates

Metabolic cold adaptation is a pattern where ectotherms from cold, high-latitude, or -altitude habitats have higher metabolic rates than ectotherms from warmer habitats. When found, metabolic cold adaptation is often attributed to countergradient selection, wherein short, cool growing seasons select for a compensatory increase in growth rates and development times of ectotherms. Yet, ectotherms in high-latitude and -altitude environments face many challenges in addition to thermal and time constraints on lifecycles. In addition to short, cool growing seasons, high-latitude and - altitude environments are characterized by regular exposure to extreme low temperatures, which cause ectotherms to enter a transient state of immobility termed chill coma. The ability to resume activity quickly after chill coma increases with latitude and altitude in patterns consistent with local adaptation to cold conditions. We show that artificial selection for fast and slow chill coma recovery among lines of the fly Drosophila melanogaster also affects rates of respiratory metabolism. Cold-hardy fly lines, with fast recovery from chill coma, had higher respiratory metabolic rates than control lines, with cold-susceptible slow-recovering lines having the lowest metabolic rates. Fast chill coma recovery was also associated with higher respiratory metabolism in a set of lines derived from a natural population. Although their metabolic rates were higher than control lines, fast-recovering cold-hardy lines did not have faster growth rates or development times than control lines. This suggests that raised metabolic rates in high-latitude and -altitude species may be driven by adaptation to extreme low temperatures, illustrating the importance of moving "Beyond the Mean".

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.

Genomic Prediction for Quantitative Traits Is Improved by Mapping Variants to Gene Ontology Categories in Drosophila melanogaster

Predicting individual quantitative trait phenotypes from high-resolution genomic polymorphism data is important for personalized medicine in humans, plant and animal breeding, and adaptive evolution. However, this is difficult for populations of unrelated individuals when the number of causal variants is low relative to the total number of polymorphisms and causal variants individually have small effects on the traits. We hypothesized that mapping molecular polymorphisms to genomic features such as genes and their gene ontology categories could increase the accuracy of genomic prediction models. We developed a genomic feature best linear unbiased prediction (GFBLUP) model that implements this strategy and applied it to three quantitative traits (startle response, starvation resistance, and chill coma recovery) in the unrelated, sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel. Our results indicate that subsetting markers based on genomic features increases the predictive ability relative to the standard genomic best linear unbiased prediction (GBLUP) model. Both models use all markers, but GFBLUP allows differential weighting of the individual genetic marker relationships, whereas GBLUP weighs the genetic marker relationships equally. Simulation studies show that it is possible to further increase the accuracy of genomic prediction for complex traits using this model, provided the genomic features are enriched for causal variants. Our GFBLUP model using prior information on genomic features enriched for causal variants can increase the accuracy of genomic predictions in populations of unrelated individuals and provides a formal statistical framework for leveraging and evaluating information across multiple experimental studies to provide novel insights into the genetic architecture of complex traits.

Does Genetic Variation Maintained by Environmental Heterogeneity Facilitate Adaptation to Novel Selection?

Environmental heterogeneity helps maintain genetic variation in fitness. Therefore, one might predict that populations living in heterogeneous environments have higher adaptive potential than populations living in homogeneous environments. Such a prediction could be useful in guiding conservation priorities without requiring detailed genetic studies. However, this prediction will be true only if the additional genetic variation maintained by environmental heterogeneity can be used to respond to novel selection. Here we examine the effect of environmental heterogeneity on future adaptability using replicated experimental Drosophila melanogaster populations that had previously evolved for ∼100 generations under one of four selective regimes: constant salt-enriched larvae medium, constant cadmium-enriched larvae medium, and two heterogeneous regimes that vary either temporally or spatially between the two media. Replicates of these experimental populations were subjected to a novel heat stress while being maintained in their original larval diet selection regimes. Adaptation to increased temperature was measured with respect to female productivity and male siring success after ∼20 generations. For female productivity, there was evidence of adaptation overall and heterogeneous populations had a larger adaptive response than homogeneous populations. There was less evidence of adaptation overall for male siring success and no support for faster adaptation in heterogeneous populations.

Artificial selection on chill-coma recovery time in Drosophila melanogaster: Direct and correlated responses to selection

Artificial selection can be used to create populations with extreme phenotypic responses to environmental stressors. When artificial selection is applied to a single component of a stress response, this selection may result in correlated responses in other stress responses, a phenomenon called cross-tolerance, which is ultimately controlled by the genetic correlations among traits. We selected for extreme responses to cold tolerance by selecting for chill-coma recovery time from a single temperate population of Drosophila melanogaster. Chill-coma recovery time is a common metric of low, but non-lethal, cold temperature tolerance. Replicated divergent artificial selection was applied to a genetically variable base population for 31 generations, resulting in two cold resistant, two cold susceptible, and two unselected control lines. To quantify the relationship between selection on chill-coma recovery and other metrics of thermal performance, we also measured survivorship after acute cold exposure, survivorship after chronic cold exposure, survivorship after cold exposure following a pre-treatment period (rapid cold hardening), starvation tolerance, and heat tolerance. We find that chill-coma recovery time is heritable within this population and that there is an asymmetric response to increased and decreased chill-coma recovery time. Surprisingly, we found no cross-tolerances between selection on chill-coma recovery time and the other environmental stress response traits. These results suggest that although artificial selection has dramatically altered chill-coma recovery time, the correlated response to selection on other stress response phenotypes has been negligible. The lack of a correlated response suggests that chill-coma recovery time in these selection lines is likely genetically independent from measures of cold survivorship tested here.

Thermal reactionomes reveal divergent responses to thermal extremes in warm and cool-climate ant species

Background

The distributions of species and their responses to climate change are in part determined by their thermal tolerances. However, little is known about how thermal tolerance evolves. To test whether evolutionary extension of thermal limits is accomplished through enhanced cellular stress response (enhanced response), constitutively elevated expression of protective genes (genetic assimilation) or a shift from damage resistance to passive mechanisms of thermal stability (tolerance), we conducted an analysis of the reactionome: the reaction norm for all genes in an organism’s transcriptome measured across an experimental gradient. We characterized thermal reactionomes of two common ant species in the eastern U.S, the northern cool-climate Aphaenogaster picea and the southern warm-climate Aphaenogaster carolinensis, across 12 temperatures that spanned their entire thermal breadth.

Results

We found that at least 2 % of all genes changed expression with temperature. The majority of upregulation was specific to exposure to low temperatures. The cool-adapted A. picea induced expression of more genes in response to extreme temperatures than did A. carolinensis, consistent with the enhanced response hypothesis. In contrast, under high temperatures the warm-adapted A. carolinensis downregulated many of the genes upregulated in A. picea, and required more extreme temperatures to induce down-regulation in gene expression, consistent with the tolerance hypothesis. We found no evidence for a trade-off between constitutive and inducible gene expression as predicted by the genetic assimilation hypothesis.

Conclusions

These results suggest that increases in upper thermal limits may require an evolutionary shift in response mechanism away from damage repair toward tolerance and prevention.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2466-z) contains supplementary material, which is available to authorized users.

Novel SNP identification in exon 3 of HSP90AA1 gene and their association with heat tolerance traits in Karan Fries (Bos taurus × Bos indicus) cows under tropical climatic condition

Heat shock proteins (HSPs) act as molecular chaperones those are preferentially transcribed in respose to heat stress and the polymorphism in HSP genes associated with heat tolerance traits in cows. HSP90AA1 gene has been mapped on Bos taurus autosome 21 (BTA-21) and spans nearly 5368 bp comprising of 11 exons out of which the first exon does not translate. The present study was done on Karan Fries (5/8 HF × 3/8 Tharparkar) cows reared in tropical climate with the objectives of identifying single-nucleotide polymorphisms (SNPs) in targeted regions (exon 3) of HSP90AA1 gene and analyzing their association with heat tolerance traits in Karan Fries cows. Respiration rate (RR) and rectal temperature (RT) were recorded once daily for four consecutive days during probable extreme hours in different seasons or temperature humidity index (THI), viz., winter, spring, and summer. For detecting single-nucleotide polymorphisms, sequence data were analyzed using BioEdit software (version 7.2). Comparative sequence analysis of HSP90AA1 gene showed point mutation, viz., g.1209A>G (exon 3) as compared to Bos taurus (NCBI Ref Seq: AC_000178.1). Association analysis indicated that THI was influenced (P < 0.01) by RR, RT, and HTC. Similarly, SNPs at locus g.1209A>G were categorized into three genotypes, i.e., AA, AG, and GG, and the least squares means (LSMEANS) of RR, RT, and HTC for GG (homozygous) genotype were significantly lower (P < 0.01) than AA (homozygous) and AG (heterozygous) genotypes. These findings may partly suggest that cows with GG genotypes were favored for heat tolerance trait, which can be used as an aid to selection for thermo-tolerance Karan Fries cows for better adaptation in subtropical and tropical hot climate.

Population genetic evidence for cold adaptation in European Drosophila melanogaster populations

We studied Drosophila melanogaster populations from Europe (the Netherlands and France) and Africa (Rwanda and Zambia) to uncover genetic evidence of adaptation to cold. We present here four lines of evidence for genes involved in cold adaptation from four perspectives: (i) the frequency of SNPs at genes previously known to be associated with chill-coma recovery time (CCRT), startle reflex (SR) and resistance to starvation stress (RSS) vary along environmental gradients and therefore among populations; (ii) SNPs of genes that correlate significantly with latitude and altitude in African and European populations overlap with SNPs that correlate with a latitudinal cline from North America; (iii) at the genomewide level, the top candidate genes are enriched in gene ontology (GO) terms that are related to cold tolerance; (iv) GO enriched terms from North American clinal genes overlap significantly with those from Africa and Europe. Each SNP was tested in 10 independent runs of Bayenv2, using the median Bayes factors to ascertain candidate genes. None of the candidate genes were found close to the breakpoints of cosmopolitan inversions, and only four candidate genes were linked to QTLs related to CCRT. To overcome the limitation that we used only four populations to test correlations with environmental gradients, we performed simulations to estimate the power of our approach for detecting selection. Based on our results, we propose a novel network of genes that is involved in cold adaptation.

Transcriptome responses to heat- and cold-stress in ladybirds (Cryptolaemus montrouzieri Mulasnt) analyzed by deep-sequencing

Background

Changed temperature not only threaten agricultural production, but they also affect individual biological behavior, population and community of many insects, and consequently reduce the stability of our ecosystem. Insect’s ability to respond to temperature stress evolved through a complex adaptive process, thus resulting in varied temperature tolerance among different insects. Both high and low extreme temperatures are detrimental to insect development since they constitute an important abiotic stress capable of inducing abnormal biological responses. Many studies on heat or cold tolerance of ladybirds have focused on measurements of physiological and biochemical indexes such as supercooling point, higher/lower lethal temperatures, survival rate, dry body weight, water content, and developmental duration. And studies of the molecular mechanisms of ladybird responses to heat or cold stress have focused on single genes, such as those encoding heat shock proteins, but has not been analyzed by transcriptome profiling.

Results

In this study, we report the use of Digital Gene Expression (DGE) tag profiling to gain insight into transcriptional events associated with heat- and cold-stress in C. montrouzieri. About 6 million tags (49 bp in length) were sequenced in a heat stress group, a cold stress group and a negative control group. We obtained 687 and 573 genes that showed significantly altered expression levels following heat and cold shock treatments, respectively. Analysis of the global gene expression pattern suggested that 42 enzyme-encoding genes mapped to many Gene Ontology terms are associated with insect’s response to heat- and cold-stress.

Conclusions

These results provide a global assessment of genes and molecular mechanisms involved in heat and cold tolerance.

Patterns of longevity and fecundity at two temperatures in a set of heat-selected recombinant inbred lines of Drosophila melanogaster

Quantitative trait loci (QTL) were mapped for longevity and fecundity at two temperatures, 20 and 30 °C, in two sets of recombinant inbred lines (RIL) highly differing in thermotolerance. Early fecundity (EF) and longevity showed a negative association between temperatures. For instance, longevity was higher and fecundity was lower in the RIL panel showing higher life span at 30 °C. One X-linked QTL (7B3-12E) co-localized for longevity and EF at 20 °C, with one QTL allele showing a positive additive effect on longevity and a negative effect on EF. The across-RIL genetic correlation between longevity and EF was not significant within each temperature, and most QTL that affect life span have no effect on EF at each temperature. EF and longevity can mostly be genetically uncoupled in the thermotolerance-divergent RIL within each temperature as opposed to between temperatures. QTL were mostly temperature specific, although some trait-specific QTL showed possible antagonistic effects between temperatures.

The Genomic Basis of Postponed Senescence in Drosophila melanogaster

Natural populations harbor considerable genetic variation for lifespan. While evolutionary theory provides general explanations for the existence of this variation, our knowledge of the genes harboring naturally occurring polymorphisms affecting lifespan is limited. Here, we assessed the genetic divergence between five Drosophila melanogaster lines selected for postponed senescence for over 170 generations (O lines) and five lines from the same base population maintained at a two week generation interval for over 850 generations (B lines). On average, O lines live 70% longer than B lines, are more productive at all ages, and have delayed senescence for other traits than reproduction. We performed population sequencing of pools of individuals from all B and O lines and identified 6,394 genetically divergent variants in or near 1,928 genes at a false discovery rate of 0.068. A 2.6 Mb region at the tip of the X chromosome contained many variants fixed for alternative alleles in the two populations, suggestive of a hard selective sweep. We also assessed genome wide gene expression of O and B lines at one and five weeks of age using RNA sequencing and identified genes with significant (false discovery rate < 0.05) effects on gene expression with age, population and the age by population interaction, separately for each sex. We identified transcripts that exhibited the transcriptional signature of postponed senescence and integrated the gene expression and genetic divergence data to identify 98 (175) top candidate genes in females (males) affecting postponed senescence and increased lifespan. While several of these genes have been previously associated with Drosophila lifespan, most are novel and constitute a rich resource for future functional validation.

Genetics of Intraspecies Variation in Avoidance Behavior Induced by a Thermal Stimulus in Caenorhabditis elegans

Individuals within a species vary in their responses to a wide range of stimuli, partly as a result of differences in their genetic makeup. Relatively little is known about the genetic and neuronal mechanisms contributing to diversity of behavior in natural populations. By studying intraspecies variation in innate avoidance behavior to thermal stimuli in the nematode Caenorhabditis elegans, we uncovered genetic principles of how different components of a behavioral response can be altered in nature to generate behavioral diversity. Using a thermal pulse assay, we uncovered heritable variation in responses to a transient temperature increase. Quantitative trait locus mapping revealed that separate components of this response were controlled by distinct genomic loci. The loci we identified contributed to variation in components of thermal pulse avoidance behavior in an additive fashion. Our results show that the escape behavior induced by thermal stimuli is composed of simpler behavioral components that are influenced by at least six distinct genetic loci. The loci that decouple components of the escape behavior reveal a genetic system that allows independent modification of behavioral parameters. Our work sets the foundation for future studies of evolution of innate behaviors at the molecular and neuronal level.

Localization of quantitative trait loci for diapause and other photoperiodically regulated life history traits important in adaptation to seasonally varying environments

Seasonally changing environments at high latitudes present great challenges for the reproduction and survival of insects, and photoperiodic cues play an important role in helping them to synchronize their life cycle with prevalent and forthcoming conditions. We have mapped quantitative trait loci (QTL) responsible for the photoperiodic regulation of four life history traits, female reproductive diapause, cold tolerance, egg-to-eclosion development time and juvenile body weight in Drosophila montana strains from different latitudes in Canada and Finland. The F2 progeny of the cross was reared under a single photoperiod (LD cycle 16:8), which the flies from the Canadian population interpret as early summer and the flies from the Finnish population as late summer. The analysis revealed a unique QTL for diapause induction on the X chromosome and several QTL for this and the other measured traits on the 4th chromosome. Flies' cold tolerance, egg-to-eclosion development time and juvenile body weight had several QTL also on the 2nd, 3rd and 5th chromosome, some of the peaks overlapping with each other. These results suggest that while the downstream output of females' photoperiodic diapause response is partly under a different genetic control from that of the other traits in the given day length, all traits also share some QTL, possibly involving genes with pleiotropic effects and/or multiple tightly linked genes. Nonoverlapping QTL detected for some of the traits also suggest that the traits are potentially capable of independent evolution, even though this may be restricted by epistatic interactions and/or correlations and trade-offs between the traits.

Accounting for genetic architecture improves sequence based genomic prediction for a Drosophila fitness trait

The ability to predict quantitative trait phenotypes from molecular polymorphism data will revolutionize evolutionary biology, medicine and human biology, and animal and plant breeding. Efforts to map quantitative trait loci have yielded novel insights into the biology of quantitative traits, but the combination of individually significant quantitative trait loci typically has low predictive ability. Utilizing all segregating variants can give good predictive ability in plant and animal breeding populations, but gives little insight into trait biology. Here, we used the Drosophila Genetic Reference Panel to perform both a genome wide association analysis and genomic prediction for the fitness-related trait chill coma recovery time. We found substantial total genetic variation for chill coma recovery time, with a genetic architecture that differs between males and females, a small number of molecular variants with large main effects, and evidence for epistasis. Although the top additive variants explained 36% (17%) of the genetic variance among lines in females (males), the predictive ability using genomic best linear unbiased prediction and a relationship matrix using all common segregating variants was very low for females and zero for males. We hypothesized that the low predictive ability was due to the mismatch between the infinitesimal genetic architecture assumed by the genomic best linear unbiased prediction model and the true genetic architecture of chill coma recovery time. Indeed, we found that the predictive ability of the genomic best linear unbiased prediction model is markedly improved when we combine quantitative trait locus mapping with genomic prediction by only including the top variants associated with main and epistatic effects in the relationship matrix. This trait-associated prediction approach has the advantage that it yields biologically interpretable prediction models.

Genetic Architecture of Micro-Environmental Plasticity in Drosophila melanogaster

Individuals of the same genotype do not have the same phenotype for quantitative traits when reared under common macro-environmental conditions, a phenomenon called micro-environmental plasticity. Genetic variation in micro-environmental plasticity is assumed in models of the evolution of phenotypic variance, and is important in applied breeding and personalized medicine. Here, we quantified genetic variation for micro-environmental plasticity for three quantitative traits in the inbred, sequenced lines of the Drosophila melanogaster Genetic Reference Panel. We found substantial genetic variation for micro-environmental plasticity for all traits, with broad sense heritabilities of the same magnitude or greater than those of trait means. Micro-environmental plasticity is not correlated with residual segregating variation, is trait-specific, and has genetic correlations with trait means ranging from zero to near unity. We identified several candidate genes associated with micro-environmental plasticity of startle response, including Drosophila Hsp90, setting the stage for future genetic dissection of this phenomenon.

The polymorphisms in the promoter of HSP90 gene and their association with heat tolerance of bay scallop

The heat shock protein 90 (HSP90) is a highly abundant and ubiquitous molecular chaperone which plays essential roles in many cellular processes. In the present study, the messenger RNA (mRNA) expressions of HSP90 after acute heat stress were investigated in two bay scallop populations (Argopecten irradians irradians and Argopecten irradians concentricus). The heat-resistant scallop A. i. concentricus, which is distributed in Zhanjiang, China, exhibited significantly higher induction of HSP90 compared with that of the heat-sensitive scallop A. i. irradians, which is distributed in Qinhuangdao, China. The promoter sequence of HSP90 gene from bay scallop (AiHSP90) was cloned, and the polymorphisms within this region were investigated by sequencing to analyze their association with heat tolerance. A total of six single nucleotide polymorphisms (SNPs), including -1167 T-C, -1023 A-C, -799 C-T, -774 A-G, -686 C-T, and -682 A-C, were identified in the amplified promoter region, and most of them affected the putative transcription factor binding sites except for locus -1167. All the six SNP sites were found to be associated with heat tolerance after Hardy-Weinberg equilibrium (HWE) and association analysis. Moreover, haplotypes CACACC and TCTATC were also found to be associated with heat tolerance based on the result of linkage disequilibrium and association analysis. The results provided insights into the molecular mechanisms underlying the thermal adaptation of different congener endemic bay scallops, which suggested that the increased heat tolerance of A. i. concentricus (compared with A. i. irradians) was associated with the higher expression of AiHSP90. Meanwhile, the six genotypes (-1167 TT, -1023 CC, -799 TT, -774 GG, -686 CC, and -682 AA) and two haplotypes (CACACC and TCTATC) could be used as potential markers for scallop selection breeding with higher heat tolerance.