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

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.

Epigenomic Changes in Ostrinia Moths Under Elevated Pupal and Adult Temperature

Epigenetic changes in the methylation of DNA may occur in response to environmental stressors, including warming climates. DNA methylation may also play an important role in regulating gene expression during both male and female reproduction in many insect species. However, it is currently unknown how DNA methylation shifts when individuals are reproducing under warmer temperatures. We exposed European corn borer moths (Ostrinia nubilalis) to heat during the pupal and adult life stages then investigated changes in DNA methylation across the genome using enzymatic methyl-seq (EM-seq). We compared methylation patterns in reproductive males and females exposed to heat (28°C) to those that experienced an ambient temperature (23°C). We found that heat exposure led to a small but significant increase in the percentage of methylated CpG sites throughout the genome in both sexes. However, DNA methylation rates were higher in females and differential methylation following heat exposure localised to unique regions in each sex. In males, methylation shifted within genes belonging to pathways including Hippo signalling, ubiquitin-mediated proteolysis, DNA damage repair and spermatogenesis. In females, differential methylation occurred in genes related to histone modification and oogenesis. Our results suggest that DNA methylation patterns respond to moderate heat exposure in Lepidoptera and provide insight into epigenetic responses to heatwaves, suggesting novel pathways that may be involved in responding to heat stress during metamorphosis and reproduction.

(Limited) Predictability of thermal adaptation in invertebrates

Evolutionary genomic approaches provide powerful tools to understand variation in and evolution of physiological processes. Untargeted genomic or transcriptomic screens can identify functionally annotated candidate genes linked to specific physiological processes, in turn suggesting evolutionary roles for these processes. Such studies often aim to inform modeling of the potential of natural populations to adapt to climate change, but these models are most accurate when evolutionary responses are repeatable, and thus predictable. Here, we synthesize the evolutionary genetic and comparative transcriptomic literature on terrestrial and marine invertebrates to assess whether evolutionary responses to temperature are repeatable within populations, across populations and across species. There is compelling evidence for repeatability, sometimes even across species. However, responses to laboratory selection and geographic variation across thermal gradients appear to be highly idiosyncratic. We also survey whether genetic/transcriptomic studies repeatedly identify candidate genes in three functional groups previously associated with the response to thermal stress: heat shock protein (Hsp) genes, proteolysis genes and immunity genes. Multiple studies across terrestrial and marine species identify candidates included in these gene sets. Yet, each of the gene sets are identified in only a minority of studies. Together, these patterns suggest that there is limited predictability of evolutionary responses to natural selection, including across studies within species. We discuss specific patterns for the candidate gene sets, implications for predictive modeling, and other potential applications of evolutionary genetics in elucidating physiology and gene function. Finally, we discuss limitations of inferences from available evolutionary genetic studies and directions for future research.

The Role of Genetic Variation in Shaping Phenotypic Responses to Diet in Aging Drosophila melanogaster

Nutrition plays a central role in healthy living, however, extensive variability in individual responses to dietary interventions complicates our understanding of its effects. Here we present a comprehensive study utilizing the Drosophila Genetic Reference Panel (DGRP), investigating how genetic variation influences responses to diet and aging. Quantitative genetic analyses of the impact of dietary restriction on lifespan, locomotor activity, dry weight, and heat knockdown time were performed. Locomotor activity, dry weight and heat knockdown time were measured on the same individual flies. We found significant genotype-by-diet interaction (GDI) and genotype-by-age interaction (GAI) for all traits. Therefore, environmental factors play a crucial role in shaping trait variation at different ages and diets, and/or distinct genetic variation influences these traits at different ages and diets. Our genome wide association study also identified a quantitative trait locus for age-dependent dietary response. The observed GDI and GAI indicates that susceptibility to environmental influences changes as organisms age, which could have significant implications for dietary recommendations and interventions aimed at promoting healthy aging in humans. The identification of associations between DNA sequence variation and age-dependent dietary responses opens new avenues for research into the genetic mechanisms underlying these interactions.

Evolution and Plasticity of Gene Expression Under Progressive Warming in Drosophila subobscura

Understanding the molecular mechanisms of thermal adaptation is crucial to predict the impacts of global warming. However, there is still a lack of research on the effects of rising temperatures over time and of studies involving different populations from the same species. The present study focuses on these two aspects, which are of great importance in understanding how organisms cope and adapt to ongoing changes in their environment. This study investigates the impact of global warming on the gene expression patterns of Drosophila subobscura populations from two different latitudinal locations after 23 generations of evolution. Our results indicate that evolutionary changes depend on the genetic background of the populations, with different starting points for thermal evolution, and that high-latitude populations show more pronounced evolutionary changes, with some evidence of convergence towards low-latitude populations. We found an interplay between plasticity and selection, with the high-latitude population showing fewer initial plastic genes and lower levels of adaptive plasticity, but a greater magnitude of change in both plastic and selective responses during evolution under warming conditions compared with its low-latitude counterpart. A substantial proportion of the transcriptome was observed to be evolving, despite the lack of observable response at higher-order phenotypic traits. The interplay between plasticity and selection may prove to be an essential component in shaping species' evolutionary responses to climate change. Furthermore, the value of conducting studies on multiple populations of the same species is emphasised, given the identification of differences between populations with different backgrounds in several contexts.

Experimental Evolution in a Warming World: The Omics Era

A comprehensive understanding of the genetic mechanisms that shape species responses to thermal variation is essential for more accurate predictions of the impacts of climate change on biodiversity. Experimental evolution with high-throughput resequencing approaches (evolve and resequence) is a highly effective tool that has been increasingly employed to elucidate the genetic basis of adaptation. The number of thermal evolve and resequence studies is rising, yet there is a dearth of efforts to integrate this new wealth of knowledge. Here, we review this literature showing how these studies have contributed to increase our understanding on the genetic basis of thermal adaptation. We identify two major trends: highly polygenic basis of thermal adaptation and general lack of consistency in candidate targets of selection between studies. These findings indicate that the adaptive responses to specific environments are rather independent. A review of the literature reveals several gaps in the existing research. Firstly, there is a paucity of studies done with organisms of diverse taxa. Secondly, there is a need to apply more dynamic and ecologically relevant thermal environments. Thirdly, there is a lack of studies that integrate genomic changes with changes in life history and behavioral traits. Addressing these issues would allow a more in-depth understanding of the relationship between genotype and phenotype. We highlight key methodological aspects that can address some of the limitations and omissions identified. These include the need for greater standardization of methodologies and the utilization of new technologies focusing on the integration of genomic and phenotypic variation in the context of thermal adaptation.

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.

"A comparison of thermal stress response between Drosophila melanogaster and Drosophila pseudoobscura reveals differences between species and sexes"

The environment is changing faster than anticipated due to climate change, making species more vulnerable to its impacts. The level of vulnerability of species is influenced by factors such as the degree and duration of exposure, as well as the physiological sensitivity of organisms to changes in their environments, which has been shown to vary among species, populations, and individuals. Here, we compared physiological changes in fecundity, critical thermalmaximum (CTmax), respiratory quotient (RQ), and DNA damage in ovaries in response to temperature stress in two species of fruit fly, Drosophila melanogaster (25 vs. 29.5 °C) and Drosophila pseudoobscura (20.5 vs. 25 °C). The fecundity of D. melanogaster was more affected by high temperatures when exposed during egg through adult development, while D. pseudoobscura was most significantly affected when exposed to high temperatures exclusively during egg through pupal development. Additionally, D. melanogaster males exhibited a decrease of CTmax under high temperatures, while females showed an increase of CTmax when exposed to high temperatures during egg through adult development. while D. pseudoobscura females and males showed an increased CTmax only when reared at high temperatures during egg through pupae development. Moreover, both species showed an acceleration in oogenesis and an increase in apoptosis due to heat stress. These changes can likely be attributed to key differences in the geographic range, thermal range, development time, and other different factors between these two systems. Through this comparison of variation in physiology and developmental response to thermal stress, we found important differences between species and sexes that suggest future work needs to account for these factors separately in understanding the effects of constant increased temperatures.

Scoring thermal limits in small insects using open-source, computer-assisted motion detection

Scoring thermal tolerance traits live or with recorded video can be time consuming and susceptible to observer bias, and as with many physiological measurements, there can be trade-offs between accuracy and throughput. Recent studies show that automated particle tracking is a viable alternative to manually scoring videos, although some of the software options are proprietary and costly. In this study, we present a novel strategy for automated scoring of thermal tolerance videos by inferring motor activity with motion detection using an open-source Python command line application called DIME (detector of insect motion endpoint). We apply our strategy to both dynamic and static thermal tolerance assays, and our results indicate that DIME can accurately measure thermal acclimation responses, generally agrees with visual estimates of thermal limits, and can significantly increase throughput over manual methods.

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

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

A predominant role of genotypic variation in both expression of sperm competition genes and paternity success in Drosophila melanogaster

Sperm competition is a crucial aspect of male reproductive success in many species, including Drosophila melanogaster, and seminal fluid proteins (Sfps) can influence sperm competitiveness. However, the combined effect of environmental and genotypic variation on sperm competition gene expression remains poorly understood. Here, we used Drosophila Genetic Reference Panel (DGRP) inbred lines and manipulated developmental population density (i.e. larval density) to test the effects of genotype, environment and genotype-by-environment interactions (GEI) on the expression of the known sperm competition genes Sex Peptide, Acp36DE and CG9997. High larval density resulted in reduced adult body size, but expression of sperm competition genes remained unaffected. Furthermore, we found no significant GEI but genotypic effects in the expression of SP and Acp36DE. Our results also revealed GEI for relative competitive paternity success (second male paternity; P2), with genes' expression positively correlated with P2. Given the effect of genotype on the expression of genes, we conducted a genome-wide association study (GWAS) and identified polymorphisms in putative cis-regulatory elements as predominant factors regulating the expression of SP and Acp36DE. The association of genotypic variation with sperm competition outcomes, and the resilience of sperm competition genes' expression against environmental challenges, demonstrates the importance of genome variation background in reproductive fitness.

To RNA-binding and beyond: Emerging facets of the role of Rbfox proteins in development and disease

The RNA-binding Fox-1 homologue (Rbfox) proteins represent an ancient family of splicing factors, conserved through evolution. All members share an RNA recognition motif (RRM), and a particular affinity for the GCAUG signature in target RNA molecules. The role of Rbfox, as a splice factor, deciding the tissue-specific inclusion/exclusion of an exon, depending on its binding position on the flanking introns, is well known. Rbfox often acts in concert with other splicing factors, and forms splicing regulatory networks. Apart from this canonical role, recent studies show that Rbfox can also function as a transcription co-factor, and affects mRNA stability and translation. The repertoire of Rbfox targets is vast, including genes involved in the development of tissue lineages, such as neurogenesis, myogenesis, and erythropoeiesis, and molecular processes, including cytoskeletal dynamics, and calcium handling. A second layer of complexity is added by the fact that Rbfox expression itself is regulated by multiple mechanisms, and, in vertebrates, exhibits tissue-specific expression. The optimum dosage of Rbfox is critical, and its misexpression is etiological to various disease conditions. In this review, we discuss the contextual roles played by Rbfox as a tissue-specific regulator for the expression of many important genes with diverse functions, through the lens of the emerging data which highlights its involvement in many human diseases. Furthermore, we explore the mechanistic details provided by studies in model organisms, with emphasis on the work with Drosophila. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > Splicing Regulation/Alternative Splicing.

Increased pupal temperature has reversible effects on thermal performance and irreversible effects on immune system and fecundity in adult ladybirds

The environmental conditions an organism encounters during development vary in their lasting impact on adult phenotypes. In the context of ongoing climate change, it is particularly relevant to understand how high developmental temperatures can impact adult traits, and whether these effects persist or diminish during adulthood. Here, we assessed the effects of pupal temperature (17 °C - normal temperature, 26 °C - increased temperature, or 35 °C - heat wave) on adult Harmonia axyridis thermal stress tolerance, immune function, starvation resistance, and fecundity. The temperature during pupation significantly affected all investigated traits in fresh adults. Heat acclimation decreased adult haemocyte concentration, cold tolerance, and total egg production, and had a positive effect on heat tolerance and starvation resistance. The negative effects of heat acclimation on cold tolerance diminished after seven days. In contrast, heat acclimation had a lasting positive effect on adult heat tolerance. Our results provide a broad assessment of the effects of developmental thermal acclimation on H. axyridis adult phenotypes. The relative plasticity of several adult traits after thermal acclimation may be consequential for the future geographic distribution and local performance of various insect species.

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.

Genetic variation of morphological scaling in Drosophila melanogaster

Morphological scaling relationships between the sizes of individual traits and the body captures the characteristic shape of a species, and their evolution is the primary mechanism of morphological diversification. However, we have almost no knowledge of the genetic variation of scaling, which is critical if we are to understand how scaling evolves. Here we explore the genetics of population scaling relationships (scaling relationships fit to multiple genetically-distinct individuals in a population) by describing the distribution of individual scaling relationships (genotype-specific scaling relationships that are unseen or cryptic). These individual scaling relationships harbor the genetic variation in the developmental mechanisms that regulate trait growth relative to body growth, and theoretical studies suggest that their distribution dictates how the population scaling relationship will respond to selection. Using variation in nutrition to generate size variation within 197 isogenic lineages of Drosophila melanogaster, we reveal extensive variation in the slopes of the wing-body and leg-body individual scaling relationships among genotypes. This variation reflects variation in the nutritionally-induced size plasticity of the wing, leg, and body. Surprisingly, we find that variation in the slope of individual scaling relationships primarily results from variation in nutritionally-induced plasticity of body size, not leg or wing size. These data allow us to predict how different selection regimes affect scaling in Drosophila, and is the first step in identifying the genetic targets of such selection. More generally, our approach provides a framework for understanding the genetic variation of scaling, an important prerequisite to explaining how selection changes scaling and morphology.

A Drosophila screen identifies a role for histone methylation in ER stress preconditioning

Stress preconditioning occurs when transient, sublethal stress events impact an organism's ability to counter future stresses. Although preconditioning effects are often noted in the literature, very little is known about the underlying mechanisms. To model preconditioning, we exposed a panel of genetically diverse Drosophila melanogaster to a sublethal heat shock and measured how well the flies survived subsequent exposure to endoplasmic reticulum (ER) stress. The impact of preconditioning varied with genetic background, ranging from dying half as fast to four and a half times faster with preconditioning compared to no preconditioning. Subsequent association and transcriptional analyses revealed that histone methylation, transcriptional regulation, and immune status are all candidate preconditioning modifier pathways. Strikingly, almost all subunits (7/8) in the Set1/COMPASS complex were identified as candidate modifiers of preconditioning. Functional analysis of Set1 knockdown flies demonstrated that loss of Set1 led to the transcriptional dysregulation of canonical ER stress genes during preconditioning. Based on these analyses, we propose a model of preconditioning in which Set1 helps to establish an interim transcriptional 'memory' of previous stress events, resulting in a preconditioned response to subsequent stress.

Signaling through the dystrophin glycoprotein complex affects the stress-dependent transcriptome in Drosophila

Deficiencies in the human dystrophin glycoprotein complex (DGC), which links the extracellular matrix with the intracellular cytoskeleton, cause muscular dystrophies, a group of incurable disorders associated with heterogeneous muscle, brain and eye abnormalities. Stresses such as nutrient deprivation and aging cause muscle wasting, which can be exacerbated by reduced levels of the DGC in membranes, the integrity of which is vital for muscle health and function. Moreover, the DGC operates in multiple signaling pathways, demonstrating an important function in gene expression regulation. To advance disease diagnostics and treatment strategies, we strive to understand the genetic pathways that are perturbed by DGC mutations. Here, we utilized a Drosophila model to investigate the transcriptomic changes in mutants of four DGC components under temperature and metabolic stress. We identified DGC-dependent genes, stress-dependent genes and genes dependent on the DGC for a proper stress response, confirming a novel function of the DGC in stress-response signaling. This perspective yields new insights into the etiology of muscular dystrophy symptoms, possible treatment directions and a better understanding of DGC signaling and regulation under normal and stress conditions.

Molecular Mechanisms of Winter Survival

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

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

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

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

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

Fluctuating heat stress during development exposes reproductive costs and putative benefits

Temperature and thermal variability are increasing worldwide, with well-known survival consequences. However, effects on other potentially more thermally sensitive reproductive traits are less understood, especially when considering thermal variation. Studying the consequences of male reproduction in the context of climate warming and ability to adapt is becoming increasingly relevant. Our goals were to test how exposure to different average temperatures that either fluctuated or remained constant impacts different male reproductive performance traits and to assess adaptive potential to future heat stress. We took advantage of a set of Drosophila melanogaster isogenic lines of different genotypes, exposing them to four different thermal conditions. These conditions represented a benign and a stressful mean temperature, applied either constantly or fluctuating around the mean and experienced during development when heat stress avoidance is hindered because of restricted mobility. We measured subsequent male reproductive performance for mating success, fertility, number of offspring produced and offspring sex ratio, and calculated the influence of thermal stress on estimated heritability and evolvability of these reproductive traits. Both costs and benefits to different thermal conditions on reproductive performance were found, with some responses varying between genotypes. Mating success improved under fluctuating benign temperature conditions and declined as temperature stress increased regardless of genotype. Fertility and productivity were severely reduced at fluctuating mean high temperature for all genotypes, but some genotypes were unaffected at constant high mean temperature. These more thermally robust genotypes showed a slight increase in productivity under the fluctuating benign condition compared to constant high temperature, despite both thermal conditions sharing the same temperature for 6 hr daily. Increasing thermal stress resulted in higher heritability and evolvability. Overall, the effects of temperature on reproductive performance depended on the trait and genotype; performance of some traits slightly increased when high temperatures were experienced for short periods but decreased substantially even when experiencing a benign temperature for a portion of each day. While thermal stress increased genetic variation that could provide adaptive potential against climate warming, this is unlikely to compensate for the overall severe negative effect on reproductive performance as mean temperature and variance increase.

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

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

Pronounced Plastic and Evolutionary Responses to Unpredictable Thermal Fluctuations in Drosophila simulans

Organisms are exposed to temperatures that vary, for example on diurnal and seasonal time scales. Thus, the ability to behaviorally and/or physiologically respond to variation in temperatures is a fundamental requirement for long-term persistence. Studies on thermal biology in ectotherms are typically performed under constant laboratory conditions, which differ markedly from the variation in temperature across time and space in nature. Here, we investigate evolutionary adaptation and environmentally induced plastic responses of Drosophila simulans to no fluctuations (constant), predictable fluctuations or unpredictable fluctuations in temperature. We whole-genome sequenced populations exposed to 20 generations of experimental evolution under the three thermal regimes and examined the proteome after short-term exposure to the same three regimes. We find that unpredictable fluctuations cause the strongest response at both genome and proteome levels. The loci showing evolutionary responses were generally unique to each thermal regime, but a minor overlap suggests either common laboratory adaptation or that some loci were involved in the adaptation to multiple thermal regimes. The evolutionary response, i.e., loci under selection, did not coincide with induced responses of the proteome. Thus, genes under selection in fluctuating thermal environments are distinct from genes important for the adaptive plastic response observed within a generation. This information is key to obtain a better understanding and prediction of the effects of future increases in both mean and variability of temperatures.

Predicting Thermal Adaptation by Looking Into Populations' Genomic Past

Molecular evolution offers an insightful theory to interpret the genomic consequences of thermal adaptation to previous events of climate change beyond range shifts. However, disentangling often mixed footprints of selective and demographic processes from those due to lineage sorting, recombination rate variation, and genomic constrains is not trivial. Therefore, here we condense current and historical population genomic tools to study thermal adaptation and outline key developments (genomic prediction, machine learning) that might assist their utilization for improving forecasts of populations' responses to thermal variation. We start by summarizing how recent thermal-driven selective and demographic responses can be inferred by coalescent methods and in turn how quantitative genetic theory offers suitable multi-trait predictions over a few generations via the breeder's equation. We later assume that enough generations have passed as to display genomic signatures of divergent selection to thermal variation and describe how these footprints can be reconstructed using genome-wide association and selection scans or, alternatively, may be used for forward prediction over multiple generations under an infinitesimal genomic prediction model. Finally, we move deeper in time to comprehend the genomic consequences of thermal shifts at an evolutionary time scale by relying on phylogeographic approaches that allow for reticulate evolution and ecological parapatric speciation, and end by envisioning the potential of modern machine learning techniques to better inform long-term predictions. We conclude that foreseeing future thermal adaptive responses requires bridging the multiple spatial scales of historical and predictive environmental change research under modern cohesive approaches such as genomic prediction and machine learning frameworks.