Strong impact of thermal environment on the quantitative genetic basis of a key stress tolerance trait

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.

Within-population variation in body size plasticity in response to combined nutritional and thermal stress is partially independent from variation in development time

Ongoing climate change has forced animals to face changing thermal and nutritional environments. Animals can adjust to such combinations of stressors via plasticity. Body size is a key trait influencing organismal fitness, and plasticity in this trait in response to nutritional and thermal conditions varies among genetically diverse, locally adapted populations. The standing genetic variation within a population can also influence the extent of body size plasticity. We generated near-isogenic lines from a newly collected population of Drosophila melanogaster at the mid-point of east coast Australia and assayed body size for all lines in combinations of thermal and nutritional stress. We found that isogenic lines showed distinct underlying patterns of body size plasticity in response to temperature and nutrition that were often different from the overall population response. We then tested whether plasticity in development time could explain, and therefore regulate, variation in body size to these combinations of environmental conditions. We selected five genotypes that showed the greatest variation in response to combined thermal and nutritional stress and assessed the correlation between response of developmental time and body size. While we found significant genetic variation in development time plasticity, it was a poor predictor of body size among genotypes. Our results therefore suggest that multiple developmental pathways could generate genetic variation in body size plasticity. Our study emphasizes the need to better understand genetic variation in plasticity within a population, which will help determine the potential for populations to adapt to ongoing environmental change.

Plasticity's role in adaptive evolution depends on environmental change components

To forecast extinction risks of natural populations under climate change and direct human impacts, an integrative understanding of both phenotypic plasticity and adaptive evolution is essential. To date, the evidence for whether, when, and how much plasticity facilitates adaptive responses in changing environments is contradictory. We argue that explicitly considering three key environmental change components - rate of change, variance, and temporal autocorrelation - affords a unifying framework of the impact of plasticity on adaptive evolution. These environmental components each distinctively effect evolutionary and ecological processes underpinning population viability. Using this framework, we develop expectations regarding the interplay between plasticity and adaptive evolution in natural populations. This framework has the potential to improve predictions of population viability in a changing world.

Are Signals of Local Environmental Adaptation Diluted by Laboratory Culture?

Insects have the ability to readily adapt to changes in environmental conditions, however the strength of local environmental adaptation signals under divergent conditions and the occurrence of trait inertia after relaxation of selection, remains poorly understood, especially for traits of climate stress resistance (CSR) and their phenotypic plasticity. The strength of environmental adaptation signals depend on several selection pressures present in the local environment, while trait inertia often occurs when there is a weakening or removal of a source of selection. Here, using Drosophila melanogaster, we asked whether signals of adaptation in CSR traits (critical thermal limits, heat and chill survival and, desiccation and starvation resistance) persist after exposure to laboratory culture for different durations (two vs. ten generations) across four climatically distinct populations. We show that culture duration has large effects on CSR traits and can both amplify or dilute signals of local adaptation. Effects were however dependent upon interactions between the source population, acclimation (adult acclimation at either 18 °C, 23 °C or 28 °C) conditions and the sex of the flies. Trait plasticity is markedly affected by the interaction between the source population, the specific acclimation conditions employed, and the duration in the laboratory. Therefore, a complex matrix of dynamic CSR trait responses is shown in space and time. Given these strong interaction effects, 'snapshot' estimates of environmental adaptation can result in misleading conclusions about the fitness consequences of climate variability.

Genotype-by-environment interactions modulate the rate of microevolution in reproductive timing in humans

Evidence from natural populations shows that changes in environmental conditions can cause rapid modifications in the evolutionary potential of phenotypes, partly through genotype-by-environment interactions (G×E). Therefore, the overall rate of microevolution should depend on fluctuations in environmental conditions, even when directional selection is sustained over several generations. We tested this hypothesis in a preindustrial human population that experienced a microevolutionary change in age at first reproduction (AFR) of mothers, using the annual infant mortality rate (IMR) as an indicator of environmental conditions during their early life. Using quantitative genetics analyses, we found that G×Es explained a nonnegligible fraction of the additive genetic variance in AFR and in relative fitness, as well as of the genetic covariance between AFR and fitness (i.e., the Robertson-Price covariance). The covariance was stronger for individuals exposed to unfavorable early-life environmental conditions. Our results unravel the presence of G×Es in an important life history trait and its impact on the rate of microevolution, which appears to have been sensitive to short-term fluctuations in local environmental conditions.

Conservation genetics as a management tool: The five best-supported paradigms to assist the management of threatened species

About 50 y ago, Crow and Kimura [An Introduction to Population Genetics Theory (1970)] and Ohta and Kimura [Genet. Res. 22, 201-204 (1973)] laid the foundations of conservation genetics by predicting the relationship between population size and genetic marker diversity. This work sparked an enormous research effort investigating the importance of population dynamics, in particular small population size, for population mean performance, population viability, and evolutionary potential. In light of a recent perspective [J. C. Teixeira, C. D. Huber, Proc. Natl. Acad. Sci. U.S.A. 118, 10 (2021)] that challenges some fundamental assumptions in conservation genetics, it is timely to summarize what the field has achieved, what robust patterns have emerged, and worthwhile future research directions. We consider theory and methodological breakthroughs that have helped management, and we outline some fundamental and applied challenges for conservation genetics.

Genetic variation for upper thermal tolerance diminishes within and between populations with increasing acclimation temperature in Atlantic salmon

Populations may counteract lasting temperature changes or recurrent extremes through plasticity or adaptation. However, it remains underexplored how outbreeding, either naturally, unintentionally, or facilitated, may modify a local response potential and whether genotype-by-environment interactions or between-trait correlations can restrict this potential. We quantified population differences and outbreeding effects, within-population genetic variation, and plasticity of these, for thermal performance proxy traits using 32 pedigreed wild, domesticated, and wild-domesticated Atlantic salmon families reared under common-garden conditions. Following exposure to ambient cold (11.6 °C) or ~4° and ~8° warmer summer temperatures, populations differed notably for body length and critical thermal maximum (CT_max) and for thermal plasticity of length, condition, and CT_max, but not for haematocrit. Line-cross analysis suggested mostly additive and some dominant outbreeding effects on means and solely additive outbreeding effects on plasticity. Heritability was detected for all traits. However, with increasing acclimation temperature, differences in CT_max between populations and CT_max heritability diminished, and CT_max breeding values re-ranked. Furthermore, CT_max and body size were negatively correlated at the genetic and phenotypic levels, and there was indirect evidence for a positive correlation between growth potential and thermal performance breadth for growth. Thus, population differences (including those between wild and domesticated populations) in thermal performance and plasticity may present a genetic resource in addition to the within-population genetic variance to facilitate, or impede, thermal adaptation. However, unfavourable genotype-by-environment interactions and negative between-trait correlations may generally hamper joint evolution in response to an increase in average temperature and temporary extremes.

Alternative models for transgenerational epigenetic inheritance: Molecular psychiatry beyond mice and man

Mental illness remains the greatest chronic health burden globally with few in-roads having been made despite significant advances in genomic knowledge in recent decades. The field of psychiatry is constantly challenged to bring new approaches and tools to address and treat the needs of vulnerable individuals and subpopulations, and that has to be supported by a continuous growth in knowledge. The majority of neuropsychiatric symptoms reflect complex gene-environment interactions, with epigenetics bridging the gap between genetic susceptibility and environmental stressors that trigger disease onset and drive the advancement of symptoms. It has more recently been demonstrated in preclinical models that epigenetics underpins the transgenerational inheritance of stress-related behavioural phenotypes in both paternal and maternal lineages, providing further supporting evidence for heritability in humans. However, unbiased prospective studies of this nature are practically impossible to conduct in humans so preclinical models remain our best option for researching the molecular pathophysiologies underlying many neuropsychiatric conditions. While rodents will remain the dominant model system for preclinical studies (especially for addressing complex behavioural phenotypes), there is scope to expand current research of the molecular and epigenetic pathologies by using invertebrate models. Here, we will discuss the utility and advantages of two alternative model organisms–Caenorhabditis elegans and Drosophila melanogaster–and summarise the compelling insights of the epigenetic regulation of transgenerational inheritance that are potentially relevant to human psychiatry.

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

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

Reconstructing hotspots of genetic diversity from glacial refugia and subsequent dispersal in Italian common toads (Bufo bufo)

Genetic diversity feeds the evolutionary process and allows populations to adapt to environmental changes. However, we still lack a thorough understanding of why hotspots of genetic diversity are so 'hot'. Here, we analysed the relative contribution of bioclimatic stability and genetic admixture between divergent lineages in shaping spatial patterns of genetic diversity in the common toad Bufo bufo along the Italian peninsula. We combined population genetic, phylogeographic and species distribution modelling (SDM) approaches to map ancestral areas, glacial refugia, and secondary contact zones. We consistently identified three phylogeographic lineages, distributed in northern, central and southern Italy. These lineages expanded from their ancestral areas and established secondary contact zones, before the last interglacial. SDM identified widespread glacial refugia in peninsular Italy, sometimes located under the present-day sea-level. Generalized linear models indicated genetic admixture as the only significant predictor of the levels of population genetic diversity. Our results show that glacial refugia contributed to preserving both levels and patterns of genetic diversity across glacial-interglacial cycles, but not to their formation, and highlight a general principle emerging in Mediterranean species: higher levels of genetic diversity mark populations with substantial contributions from multiple genetic lineages, irrespective of the location of glacial refugia.

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.

Distinct cold tolerance traits independently vary across genotypes in Drosophila melanogaster

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

Intraspecific variation in thermal acclimation and tolerance between populations of the winter ant, Prenolepis imparis

Thermal phenotypic plasticity, otherwise known as acclimation, plays an essential role in how organisms respond to short-term temperature changes. Plasticity buffers the impact of harmful temperature changes; therefore, understanding variation in plasticity in natural populations is crucial for understanding how species will respond to the changing climate. However, very few studies have examined patterns of phenotypic plasticity among populations, especially among ant populations. Considering that this intraspecies variation can provide insight into adaptive variation in populations, the goal of this study was to quantify the short-term acclimation ability and thermal tolerance of several populations of the winter ant, Prenolepis imparis. We tested for correlations between thermal plasticity and thermal tolerance, elevation, and body size. We characterized the thermal environment both above and below ground for several populations distributed across different elevations within California, USA. In addition, we measured the short-term acclimation ability and thermal tolerance of those populations. To measure thermal tolerance, we used chill-coma recovery time (CCRT) and knockdown time as indicators of cold and heat tolerance, respectively. Short-term phenotypic plasticity was assessed by calculating acclimation capacity using CCRT and knockdown time after exposure to both high and low temperatures. We found that several populations displayed different chill-coma recovery times and a few displayed different heat knockdown times, and that the acclimation capacities of cold and heat tolerance differed among most populations. The high-elevation populations displayed increased tolerance to the cold (faster CCRT) and greater plasticity. For high-temperature tolerance, we found heat tolerance was not associated with altitude; instead, greater tolerance to the heat was correlated with increased plasticity at higher temperatures. These current findings provide insight into thermal adaptation and factors that contribute to phenotypic diversity by revealing physiological variance among populations.

Life-History Evolution and the Genetics of Fitness Components in Drosophila melanogaster

Life-history traits or "fitness components"-such as age and size at maturity, fecundity and fertility, age-specific rates of survival, and life span-are the major phenotypic determinants of Darwinian fitness. Analyzing the evolution and genetics of these phenotypic targets of selection is central to our understanding of adaptation. Due to its simple and rapid life cycle, cosmopolitan distribution, ease of maintenance in the laboratory, well-understood evolutionary genetics, and its versatile genetic toolbox, the "vinegar fly" Drosophila melanogaster is one of the most powerful, experimentally tractable model systems for studying "life-history evolution." Here, I review what has been learned about the evolution and genetics of life-history variation in D. melanogaster by drawing on numerous sources spanning population and quantitative genetics, genomics, experimental evolution, evolutionary ecology, and physiology. This body of work has contributed greatly to our knowledge of several fundamental problems in evolutionary biology, including the amount and maintenance of genetic variation, the evolution of body size, clines and climate adaptation, the evolution of senescence, phenotypic plasticity, the nature of life-history trade-offs, and so forth. While major progress has been made, important facets of these and other questions remain open, and the D. melanogaster system will undoubtedly continue to deliver key insights into central issues of life-history evolution and the genetics of adaptation.

Using natural laboratories to study evolution to global warming: contrasting altitudinal, latitudinal, and urbanization gradients

Demonstrating the likelihood of evolution in response to global warming is important, yet challenging. We discuss how three spatial thermal gradients (latitudinal, altitudinal, and urbanization) can be used as natural laboratories to inform about the gradual thermal evolution of populations by applying a space-for-time substitution (SFTS) approach. We compare thermal variables and confounding non-thermal abiotic variables, methodological approaches and evolutionary aspects associated with each type of gradient. On the basis of an overview of recent insect studies, we show that a key assumption of SFTS, local thermal adaptation along these gradients, is often but not always met, requiring explicit validation. To increase realism when applying SFTS, we highlight the importance of integrating daily temperature fluctuations, multiple stressors and multiple interacting species. Finally, comparative studies, especially across gradient types, are important to provide more robust inferences of evolution under gradual global warming. Integrating these research directions will further strengthen the still underused, yet powerful SFTS approach to infer gradual evolution under global warming.

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.

Adaptation to environmental stress at different timescales

Environments are changing rapidly, and to cope with these changes, organisms have to adapt. Adaptation can take many shapes and occur at different speeds, depending on the type of response, the trait, the population, and the environmental conditions. The biodiversity crisis that we are currently facing illustrates that numerous species and populations are not capable of adapting with sufficient speed to ongoing environmental changes. Here, we discuss current knowledge on the ability of animals and plants to adapt to environmental stress on different timescales, mainly focusing on thermal stress and ectotherms. We discuss within-generation responses that can be fast and induced within minutes or hours, evolutionary adaptations that are often slow and take several generations, and mechanisms that lay somewhere in between and that include epigenetic transgenerational effects. To understand and predict the impacts of environmental change and stress on biodiversity, we suggest that future studies should (1) have an increased focus on understanding the type and speed of responses to fast environmental changes; (2) focus on the importance of environmental fluctuations and the predictability of environmental conditions on adaptive capabilities, preferably in field studies encompassing several fitness components; and (3) look at ecosystem responses to environmental stress and their resilience when disturbed.