QTL for the thermotolerance effect of heat hardening, knockdown resistance to heat and chill-coma recovery in an intercontinental set of recombinant inbred lines of Drosophila melanogaster

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.

Mating success at elevated temperature is associated to thermal adaptation in a set of recombinant inbred lines of Drosophila melanogaster

In insects, mating ability at elevated temperature can be relevant for adaptation to heat-stressed environments and global warming. Here, we examined copulation latency (T1), copulation duration (T2), and mating frequency (T3, an index of mating success) in two related sets of recombinant inbred lines (RIL) in Drosophila melanogaster at both elevated (33 °C) and benign (25 °C) temperatures. One of these RIL sets (RIL-SH2) was shown to be consistently more resistant in both heat knockdown and heat-shock survival assays than its related set (RIL-D48) in previous studies. Negative correlations across RILs were found between T1 and T3 in this study. Flies from the heat-resistant set of RIL (RIL-SH2) were better able to mate at elevated temperature than flies from the heat-susceptible set (RIL-D48). Quantitative trait locus (QTL) mapping identified temperature-dependent QTLs for all traits (T1, T2 and T3) on all the three major chromosomes. Mating success at elevated temperature was found to be influenced by multiple QTLs. At elevated temperature, several QTLs for mating traits co-localized with QTLs that were previously associated with thermotolerance. The genetic basis for T1, T2 and T3 at the elevated temperature was found to be largely different from the genetic basis controlling the variation for mating success at benign temperature, as there was only a very low (or even null) number of QTLs overlapping across 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.

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.

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.

Negative genetic correlation between longevity and its hormetic extension by dietary restriction in Drosophila melanogaster

Longevity is a highly malleable trait which is influenced by many genetic and environmental factors including nutrition. Mild stress of dietary restriction (DR) is often beneficial by extending longevity in many organisms. Here, DR-induced effects on longevity were tested for genetic variation in a set of recombinant inbred lines (RIL) in D. melanogaster. Genetic variability was significant in the longevity response following a DR-treatment across RIL, with detrimental effects in several RIL but beneficial effects in other RIL. One quantitative trait locus (QTL) was consistently significant in the middle of chromosome 2 for DR-induced changes in longevity, including hormesis (an increase in longevity by DR). Another QTL co-localized with a previously found QTL for starvation resistance in females. Several other QTL were also significant on most chromosomal arms. Longevity in controls was negatively correlated to DR effects across RIL for longevity in females, the sex showing higher DR-induced hormesis. This negative genetic correlation highlights the importance to further investigate the effects of genetic variation in the strength of DR-induced hormesis in longevity and its sex-specificity.

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.

Are model organisms representative for climate change research? Testing thermal tolerance in wild and laboratory zebrafish populations

Model organisms can be useful for studying climate change impacts, but it is unclear whether domestication to laboratory conditions has altered their thermal tolerance and therefore how representative of wild populations they are. Zebrafish in the wild live in fluctuating thermal environments that potentially reach harmful temperatures. In the laboratory, zebrafish have gone through four decades of domestication and adaptation to stable optimal temperatures with few thermal extremes. If maintaining thermal tolerance is costly or if genetic traits promoting laboratory fitness at optimal temperature differ from genetic traits for high thermal tolerance, the thermal tolerance of laboratory zebrafish could be hypothesized to be lower than that of wild zebrafish. Furthermore, very little is known about the thermal environment of wild zebrafish and how close to their thermal limits they live. Here, we compared the acute upper thermal tolerance (critical thermal maxima; CT_max) of wild zebrafish measured on-site in West Bengal, India, to zebrafish at three laboratory acclimation/domestication levels: wild-caught, F₁ generation wild-caught and domesticated laboratory AB-WT line. We found that in the wild, CT_max increased with increasing site temperature. Yet at the warmest site, zebrafish lived very close to their thermal limit, suggesting that they may currently encounter lethal temperatures. In the laboratory, acclimation temperature appeared to have a stronger effect on CT_max than it did in the wild. The fish in the wild also had a 0.85-1.01°C lower CT_max compared to all laboratory populations. This difference between laboratory-held and wild populations shows that environmental conditions can affect zebrafish's thermal tolerance. However, there was no difference in CT_max between the laboratory-held populations regardless of the domestication duration. This suggests that thermal tolerance is maintained during domestication and highlights that experiments using domesticated laboratory-reared model species can be appropriate for addressing certain questions on thermal tolerance and global warming impacts.

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.

Hawaiian picture-winged Drosophila exhibit adaptive population divergence along a narrow climatic gradient on Hawaii Island

Anthropogenic influences on global processes and climatic conditions are increasingly affecting ecosystems throughout the world.Hawaii Island's native ecosystems are well studied and local long-term climatic trends well documented, making these ecosystems ideal for evaluating how native taxa may respond to a warming environment.This study documents adaptive divergence of populations of a Hawaiian picture-winged Drosophila, D. sproati, that are separated by only 7 km and 365 m in elevation.Representative laboratory populations show divergent behavioral and physiological responses to an experimental low-intensity increase in ambient temperature during maturation. The significant interaction of source population by temperature treatment for behavioral and physiological measurements indicates differential adaptation to temperature for the two populations.Significant differences in gene expression among males were mostly explained by the source population, with eleven genes in males also showing a significant interaction of source population by temperature treatment.The combined behavior, physiology, and gene expression differences between populations illustrate the potential for local adaptation to occur over a fine spatial scale and exemplify nuanced response to climate change.

A transcriptomics assessment of oxygen-temperature interactions reveals novel candidate genes underlying variation in thermal tolerance and survival

While single stress responses are fairly well researched, multiple, interactive stress responses are not-despite the obvious importance thereof. Here, using D. melanogaster, we investigated the effects of simultaneous exposures to low O₂ (hypoxia) and varying thermal conditions on mortality rates, estimates of thermal tolerance and the transcriptome. We used combinations of 21 (normoxia), 10 or 5kPa O₂ with control (23°C), cold (4°C) or hot (31°C) temperature exposures before assaying chill coma recovery time (CCRT) and heat knock down time (HKDT) as measures of cold and heat tolerance respectively. We found that mortality was significantly affected by temperature, oxygen partial pressure (PO₂) and the interaction between the two. Cold treatments resulted in low mortality (<5%), regardless of PO₂ treatment; while hot treatments resulted in higher mortality (∼20%), especially at 5kPa O₂ which was lethal for most flies (∼80%). Both CCRT and HKDT were significantly affected by temperature, but not PO₂, of the treatments, and the interaction of temperature and PO₂ was non-significant. Hot treatments led to significantly longer CCRT, and shorter HKDT in comparison to cold treatments. Global gene expression profiling provided the first transcriptome level response to the combined stress of PO₂ and temperature, showing that stressful treatments resulted in higher mortality and induced transcripts that were associated with protein kinases, catabolic processes (proteases, hydrolases, peptidases) and membrane function. Several genes and pathways that may be responsible for the protective effects of combined PO₂ and cold treatments were identified. We found that urate oxidase was upregulated in all three cold treatments, regardless of the PO₂. Small heat shock proteins Hsp22 and Hsp23 were upregulated after both 10 and 21kPa O₂-hot treatments. Collectively, the data from PO₂-hot treatments suggests that hypoxia does exacerbate heat stress, through an as yet unidentified mechanism. Hsp70B and an unannotated transcript (CG6733) were significantly differentially expressed after 5kPa O₂-cold and 10kPa O₂-hot treatments relative to their controls. Downregulation of these transcripts was correlated with reduced thermal tolerance (longer CCRT and shorter HKDT), suggesting that these genes may be important candidates for future research.

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.

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.

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.

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.

What Can Plasticity Contribute to Insect Responses to Climate Change?

Plastic responses figure prominently in discussions on insect adaptation to climate change. Here we review the different types of plastic responses and whether they contribute much to adaptation. Under climate change, plastic responses involving diapause are often critical for population persistence, but key diapause responses under dry and hot conditions remain poorly understood. Climate variability can impose large fitness costs on insects showing diapause and other life cycle responses, threatening population persistence. In response to stressful climatic conditions, insects also undergo ontogenetic changes including hardening and acclimation. Environmental conditions experienced across developmental stages or by prior generations can influence hardening and acclimation, although evidence for the latter remains weak. Costs and constraints influence patterns of plasticity across insect clades, but they are poorly understood within field contexts. Plastic responses and their evolution should be considered when predicting vulnerability to climate change-but meaningful empirical data lag behind theory.

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.

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.

Fine-mapping and selective sweep analysis of QTL for cold tolerance in Drosophila melanogaster

There is a growing interest in investigating the relationship between genes with signatures of natural selection and genes identified in QTL mapping studies using combined population and quantitative genetics approaches. We dissected an X-linked interval of 6.2 Mb, which contains two QTL underlying variation in chill coma recovery time (CCRT) in Drosophila melanogaster from temperate (European) and tropical (African) regions. This resulted in two relatively small regions of 131 kb and 124 kb. The latter one co-localizes with a very strong selective sweep in the European population. We examined the genes within and near the sweep region individually using gene expression analysis and P-element insertion lines. Of the genes overlapping with the sweep, none appears to be related to CCRT. However, we have identified a new candidate gene of CCRT, brinker, which is located just outside the sweep region and is inducible by cold stress. We discuss these results in light of recent population genetics theories on quantitative traits.

Developmental thermal plasticity among Drosophila melanogaster populations

Many biotic and abiotic variables influence the dispersal and distribution of organisms. Temperature has a major role in determining these patterns because it changes daily, seasonally and spatially, and these fluctuations have a significant impact on an organism's behaviour and fitness. Most ecologically relevant phenotypes that are adaptive are also complex and thus they are influenced by many underlying loci that interact with the environment. In this study, we quantified the degree of thermal phenotypic plasticity within and among populations by measuring chill-coma recovery times of lines reared from egg to adult at two different environmental temperatures. We used sixty genotypes from six natural populations of Drosophila melanogaster sampled along a latitudinal gradient in South America. We found significant variation in thermal plasticity both within and among populations. All populations exhibit a cold acclimation response, with flies reared at lower temperatures having increased resistance to cold. We tested a series of environmental parameters against the variation in population mean thermal plasticity and discovered the mean thermal plasticity was significantly correlated with altitude of origin of the population. Pairing our data with previous experiments on viability fitness assays in the same populations in fixed and variable environments suggests an adaptive role of this thermal plasticity in variable laboratory environments. Altogether, these data demonstrate abundant variation in adaptive thermal plasticity within and among populations.

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

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

Heat stress survival in the pre-adult stage of the life cycle in an intercontinental set of recombinant inbred lines of Drosophila melanogaster

In insects, pre-adult stages of the life cycle are exposed to variation in temperature that may differ from that in adults. However, the genetic basis for adaptation to environmental temperature could be similar between the pre-adult and the adult stages of the life cycle. Here, we tested quantitative trait loci (QTL) for heat-stress survival in larvae of D. melanogaster, with and without a mild-heat-stress pre-treatment. Two sets of recombinant inbred lines derived from lines artificially selected for high and low levels of knockdown resistance to high temperature in young flies were used as mapping population. There was no apparent increase in heat-shock survival between heat-pretreated and non-pretreated larvae. There was a positive correlation between the two experimental conditions of heat-shock survival (with and without a heat pre-treatment) except for males from one set of lines. Several QTL were identified involving all three major chromosomes. Many QTL for larval thermotolerance overlapped with thermotolerance-QTL identified in previous studies for adults. One new thermotolerance-QTL was found but these QTL explained only a small fraction of the phenotypic variance and were only significant in larvae that received no heat pre-treatment. Several candidate genes mapped within QTL ranges. We discuss an overall co-localization for thermotolerance-QTL between the adult fly in previous studies and the pre-adult stage of the life cycle in this study.

Quantitative trait loci for response to ethanol in an intercontinental set of recombinant inbred lines of Drosophila melanogaster

Alcohol, a drug widely abused, impacts the central nervous system functioning of diverse organisms. The behavioral responses to acute alcohol exposure are remarkably similar among humans and fruit flies. In its natural environment, rich in fermentation products, the fruit fly Drosophila melanogaster encounters relatively high levels of ethanol. The effects of ethanol and its metabolites on Drosophila have been studied for decades, as a model for adaptive evolution. Although extensive work has been done for elucidating patterns of genetic variation, substantially less is known about the genomic regions or genes that underlie the genetic variation of this important trait. To identify regions containing genes involved in the responses to ethanol, we used a mapping population of recombinant inbred (RIL) lines to map quantitative trait loci (QTL) that affect variation in resistance and recovery from ethanol sedation in adults and ethanol resistance in larvae. We mapped fourteen QTL affecting the response to ethanol on the three chromosomes. Seven of the QTL influence the resistance to ethanol in adults, two QTL are related to ethanol-coma recovery in adults and five affect the survival to ethanol in larvae. Most of the QTL were trait specific, suggesting that overlapping but generally unique genetic architectures underlie each trait. Each QTL explained up to 16.8% of the genetic variance among lines. Potential candidate loci contained within our QTL regions were identified and analyzed.

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

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

Genetics of climate change adaptation

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

QTL for survival to UV-C radiation in Drosophila melanogaster

PURPOSE

The aim of this study was to investigate tolerance to UV-C (ultraviolet C, 280-100 nm) radiation in Drosophila melanogaster, implementing a quantitative trait locus (QTL) mapping approach. This is of interest to test for genetic variation in survival to UV (ultraviolet) radiation.

MATERIALS AND METHODS

We performed a QTL scan in D. melanogaster recombinant inbred lines (RIL) constructed from parental stocks derived from a crossing between northern and southern hemisphere populations that segregated substantial genetic variation in thermal resistance in a previous study. Here, two experimental treatments were implemented: Continuous and cyclic UV-C radiation.

RESULTS

Significant QTL were detected on all three major chromosomes. Among these, multiple trait composite interval mapping revealed a significant QTL in the pericentromeric region of chromosome 2, a genome region consistently implicated in thermotolerance in previous studies.

CONCLUSIONS

This study shows substantial genetic variation for UV-C radiation resistance in D. melanogaster, with QTL for survival to UV-C radiation generally overlapping with major thermotolerance QTL. The genetic architecture of UV-C radiation resistance appears to be more complex in continuously irradiated individuals.

Survival of heat stress with and without heat hardening in Drosophila melanogaster: interactions with larval density

Survival of a potentially lethal high temperature stress is a genetically variable thermal adaptation trait in many organisms. Organisms cope with heat stress by basal or induced thermoresistance. Here, we tested quantitative trait loci (QTL) for heat stress survival (HSS) in Drosophila melanogaster, with and without a cyclic heat-hardening pre-treatment, for flies that were reared at low (LD) or high (HD) density. Mapping populations were two panels of recombinant inbred lines (RIL), which were previously constructed from heat stress-selected stocks: RIL-D48 and RIL-SH2, derived from backcrosses to stocks of low and high heat resistance, respectively. HSS increased with heat hardening in both LD and HD flies. In addition, HSS increased consistently with density in non-hardened flies. There was a significant interaction between heat hardening and density effects in RIL-D48. Several QTL were significant for both density and hardening treatments. Many QTL overlapped with thermotolerance QTL identified for other traits in previous studies based on LD cultures only. However, three new QTL were found in HD only (cytological ranges: 12E–16F6; 30A3–34C2; 49C–50C). Previously found thermotolerance QTL were also significant for flies from HD cultures.

Cold tolerance and cold-induced modulation of gene expression in two Drosophila virilis group species with different distributions

The importance of high and low temperature tolerance in adaptation to changing environmental conditions has evoked new interest in modulations in gene expression and metabolism linked with stress tolerance. We investigated the effects of rapid cold hardening and cold acclimatization on the chill coma recovery times of two Drosophila virilis group species, Drosophila montana and D. virilis, with different distributions and utilized a candidate gene approach to trace changes in their gene expression during and after the cold treatments. The study showed that cold acclimatization clearly decreases chill coma recovery times in both species, whereas rapid cold hardening did not have a significant effect. Microarray analysis revealed several genes showing expression changes during different stages of cold response. Amongst the 219 genes studied, two genes showed rather consistent expression changes: hsr-omega, which was up-regulated in both study species during cold acclimatization, and Eip71CD, which was down-regulated in nearly all of the cold treatments. In addition, 29 genes showed expression changes that were more treatment- and/or species specific. Overall, different stages of cold response elicited changes mainly in genes involved in heat shock response, circadian rhythm and metabolism.

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

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

Consistent effects of a major QTL for thermal resistance in field-released Drosophila melanogaster

Molecular genetic markers can be used to identify quantitative trait loci (QTL) for thermal resistance and this has allowed characterization of a major QTL for knockdown resistance to high temperature in Drosophila melanogaster. The QTL showed trade-off associations with cold resistance under laboratory conditions. However, assays of thermal tolerance conducted in the laboratory may not necessarily reflect performance at varying temperatures in the field. Here we tested if lines with different genotypes in this QTL show different thermal performance under high and low temperatures in the field using a release recapture assay. We found that lines carrying the QTL genotype for high thermal tolerance were significantly better at locating resources in the field releases under hot temperatures while the QTL line carrying the contrasting genotype were superior at cold temperatures. Further, we studied copulatory success between the different QTL genotypes at different temperatures. We found higher copulatory success in males of the high tolerance QTL genotype under hot temperature conditions, while there was no difference in females at cold temperatures. The results allow relating components of field fitness at different environmental temperatures with genotypic variation in a QTL for thermal tolerance.

Quantitative trait loci for longevity in heat-stressed Drosophila melanogaster

Longevity is a typical quantitative trait which is influenced by multiple genes. Here we explore the genetic variation in longevity of Drosophila melanogaster in both mildly heat-stressed and control flies. Quantitative trait loci (QTL) analysis for longevity was performed in a single-sex environment at 25°C with and without a mild heat-stress pre-treatment, using a previously reported set of recombinant inbred lines (RIL). QTL regions for longevity in heat-stressed flies overlapped with QTL for longevity in control flies. All longevity QTL co-localized with QTL for longevity identified in previous studies using very different sets of RIL in mixed sex environments, though the genome is nearly saturated with QTL for longevity when considering all previous studies. Heat stress decreased the number of significant QTL for longevity if compared to the control environment. Our mild heat-stress pre-treatment had a beneficial effect (hormesis) more often in shorter-lived than in longer-lived RIL.

Genome-wide deficiency screen for the genomic regions responsible for heat resistance in Drosophila melanogaster

Background

Temperature adaptation is one of the most important determinants of distribution and population size of organisms in nature. Recently, quantitative trait loci (QTL) mapping and gene expression profiling approaches have been used for detecting candidate genes for heat resistance. However, the resolution of QTL mapping is not high enough to examine the individual effects of various genes in each QTL. Heat stress-responsive genes, characterized by gene expression profiling studies, are not necessarily responsible for heat resistance. Some of these genes may be regulated in association with the heat stress response of other genes.

Results

To evaluate which heat-responsive genes are potential candidates for heat resistance with higher resolution than previous QTL mapping studies, we performed genome-wide deficiency screen for QTL for heat resistance. We screened 439 isogenic deficiency strains from the DrosDel project, covering 65.6% of the Drosophila melanogaster genome in order to map QTL for thermal resistance. As a result, we found 19 QTL for heat resistance, including 3 novel QTL outside the QTL found in previous studies.

Conclusion

The QTL found in this study encompassed 19 heat-responsive genes found in the previous gene expression profiling studies, suggesting that they were strong candidates for heat resistance. This result provides new insights into the genetic architecture of heat resistance. It also emphasizes the advantages of genome-wide deficiency screen using isogenic deficiency libraries.

Identification of X-linked quantitative trait loci affecting cold tolerance in Drosophila melanogaster and fine mapping by selective sweep analysis

Drosophila melanogaster is a cosmopolitan species that colonizes a great variety of environments. One trait that shows abundant evidence for naturally segregating genetic variance in different populations of D. melanogaster is cold tolerance. Previous work has found quantitative trait loci (QTL) exclusively on the second and the third chromosomes. To gain insight into the genetic architecture of cold tolerance on the X chromosome and to compare the results with our analyses of selective sweeps, a mapping population was derived from a cross between substitution lines that solely differed in the origin of their X chromosome: one originates from a European inbred line and the other one from an African inbred line. We found a total of six QTL for cold tolerance factors on the X chromosome of D. melanogaster. Although the composite interval mapping revealed slightly different QTL profiles between sexes, a coherent model suggests that most QTL overlapped between sexes, and each explained around 5-14% of the genetic variance (which may be slightly overestimated). The allelic effects were largely additive, but we also detected two significant interactions. Taken together, this provides evidence for multiple QTL that are spread along the entire X chromosome and whose effects range from low to intermediate. One detected transgressive QTL influences cold tolerance in different ways for the two sexes. While females benefit from the European allele increasing their cold tolerance, males tend to do better with the African allele. Finally, using selective sweep mapping, the candidate gene CG16700 for cold tolerance colocalizing with a QTL was identified.

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

A comprehensive understanding of the genetic basis of phenotypic adaptation in nature requires the identification of the functional allelic variation underlying adaptive phenotypes. The manner in which organisms respond to temperature extremes is an adaptation in many species. In the current study, we investigate the role of molecular variation in senescence marker protein-30 (Smp-30) on natural phenotypic variation in cold tolerance in Drosophila melanogaster. Smp-30 encodes a product that is thought to be involved in the regulation of Ca2+ ion homeostasis and has been shown previously to be differentially expressed in response to cold stress. Thus, we sought to assess whether molecular variation in Smp-30 was associated with natural phenotypic variation in cold tolerance in a panel of naturally derived inbred lines from a population in Raleigh, North Carolina. We identified four non-coding polymorphisms that were strongly associated with natural phenotypic variation in cold tolerance. Interestingly, two polymorphisms that were in close proximity to one another (2 bp apart) exhibited opposite phenotypic effects. Consistent with the maintenance of a pair of antagonistically acting polymorphisms, tests of molecular evolution identified a significant excess of maintained variation in this region, suggesting balancing selection is acting to maintain this variation. These results suggest that multiple mutations in non-coding regions can have significant effects on phenotypic variation in adaptive traits within natural populations, and that balancing selection can maintain polymorphisms with opposite effects on phenotypic variation.

Three selections are better than one: clinal variation of thermal QTL from independent selection experiments in Drosophila

We report the results of two independent selection experiments that have exposed distinct populations of Drosophila melanogaster to different forms of thermal selection. A recombinant population derived from Arvin California and Zimbabwe isofemale lines was exposed to laboratory natural selection at two temperatures (T(AZ): 18°C and 28°C). Microsatellite mapping identified quantitative trait loci (QTL) on the X-chromosome between the replicate "Hot" and "Cold" populations. In a separate experiment, disruptive selection was imposed on an outbred California population for the "knockdown" temperature (T(KD)) in a thermal column. Microsatellite mapping of the "High" and "Low" populations also uncovered primarily X-linked QTL. Notably, a marker in the shaggy locus at band 3A was significantly differentiated in both experiments. Finer scale mapping of the 3A region has narrowed the QTL to the shaggy gene region, which contains several candidate genes that function in circadian rhythms. The same allele that was increased in frequency in the High T(KD) populations is significantly clinal in North America and is more common at the warm end of the cline (Florida vs. Maine; however, the cline was not apparent in Australia). Together, these studies show that independent selection experiments can uncover the same target of selection and that evolution in the laboratory can recapitulate putatively adaptive clinal variation in nature.

Combined expression patterns of QTL-linked candidate genes best predict thermotolerance in Drosophila melanogaster

Knockdown resistance to high temperature (KRHT) is a thermal adaptation trait in Drosophila melanogaster. Here we used quantitative real-time PCR (qRT-PCR) to test for possible associations between KRHT and the expression of candidate genes within quantitative trait loci (QTL) in eight recombinant inbred lines (RIL). hsp60 and hsc70-3 map within an X-linked QTL, while CG10383, catsup, ddc, trap1, and cyp6a13 are linked in a KRHT-QTL on chromosome 2. hsc70-3 expression increased by heat-hardening. Principal Components analysis revealed that catsup, ddc and trap1 were either co-expressed or combined in their expression levels. This composite expression variable (e-PC1) was positively associated to KRHT in non-hardened RIL. In heat-hardened flies, hsp60 was negatively related to hsc70-3 on e-PC2, with effects on KRHT. These results are consistent with the notion that QTL can be shaped by expression variation in combined candidate loci. We found composite variables of gene expression (e-PCs) that best correlated to KRHT. Network effects with other untested linked loci are apparent because, in spite of their associations with KRHT phenotypes, e-PCs were sometimes uncorrelated with their QTL genotype.

Variations in morphological and life-history traits under extreme temperatures in Drosophila ananassae

Using half-sib analysis, we analysed the consequences of extreme rearing temperatures on genetic and phenotypic variations in the morphological and life-history traits of Drosophila ananassae. Paternal half-sib covariance contains a relatively small proportion of the epistatic variance and lacks the dominance variance and variance due to maternal effect, which provides more reliable estimates of additive genetic variance. Experiments were performed on a mass culture population of D. ananassae collected from Kanniyakumari (India). Two extremely stressful temperatures (18 degree C and 32 degree C) and one standard temperature (25 degree C) were used to examine the effect of stressful and non-stressful environments on the morphological and life-history traits in males and females. Mean values of various morphological traits differed signifi cantly among different temperature regimens in both males and females. Rearing at 18 degree C and 32 degree C resulted in decreased thorax length, wing-to-thorax (w/t) ratio, sternopleural bristle number, ovariole number, sex comb-tooth number and testis length. Phenotypic variances increased under stressful temperatures in comparison with non-stressful temperatures. Heritability and evolvability based on among-sires (males), among-dams (females), and the sum of the two components (sire + dam) showed higher values at both the stressful temperatures than at the non-stressful temperature. These differences reflect changes in additive genetic variance. Viability was greater at the high than the low extreme temperature. As viability is an indicator of stress, we can assume that stress was greater at 18 degree C than at 32 degree C in D. ananassae. The genetic variations for all the quantitative and life-history traits were higher at low temperature. Variation in sexual traits was more pronounced as compared with other morphometric traits, which shows that sexual traits are more prone to thermal stress. Our results agree with the hypothesis that genetic variation is increased in stressful environments.