Spatiotemporal dynamics and genome-wide association genome-wide association analysis of desiccation tolerance in Drosophila melanogaster

Elemental stoichiometry and insect chill tolerance: evolved and plastic changes in organismal Na+ and K+ content in Drosophila

Acclimation and evolutionary adaptation can produce phenotypic changes that allow organisms to cope with challenges. Determining the relative contributions and the underlying mechanisms driving phenotypic shifts from acclimation and adaptation is of central importance to understanding animal responses to change. Rates of evolution have traditionally been considered slow relative to ecological processes that shape biodiversity. Many organisms nonetheless show patterns of genetic variation that suggest that adaptation may act sufficiently fast to allow continuous change in phenotypes in response to environmental change (called 'adaptive tracking'). In Drosophila, both plastic and evolved differences in chill tolerance are associated with ionoregulation. Here, we combine an acclimation experiment, field collections along a well-characterized latitudinal cline, and a replicated field experiment to assess the concordance in the direction, magnitude, and potential mechanisms of acclimation and adaptation on chill coma recovery and elemental (Na and K) stoichiometry in both sexes of Drosophila melanogaster. Acclimation strongly shaped chill coma recovery, spatial adaptation produced comparatively modest effects, and temporal adaptation had no significant effect. Leveraging knowledge on the mechanisms underlying variation in chill tolerance traits, we find that relationships between elemental stoichiometry and chill coma recovery in the context of acclimation may differ from those that are associated with spatial adaptive change.

Drosophila melanogaster pigmentation demonstrates adaptive phenotypic parallelism but genomic unpredictability over multiple timescales

Populations are capable of responding to environmental change over ecological timescales via adaptive tracking. However, the translation from patterns of allele frequency change to rapid adaptation of complex traits remains unresolved. We used abdominal pigmentation in Drosophila melanogaster as a model phenotype to address the nature, genetic architecture, and repeatability of rapid adaptation in the field. We show that D. melanogaster pigmentation evolves as a highly parallel and deterministic response to shared environmental gradients across latitude and season in natural North American populations. We then experimentally evolved replicate, genetically diverse fly populations in field mesocosms to remove any confounding effects of demography and/or cryptic structure that may drive patterns in wild populations; we show that pigmentation rapidly responds, in parallel, in fewer than ten generations. Thus, pigmentation evolves concordantly in response to spatial and temporal climatic gradients. We next examined whether phenotypic differentiation was associated with allele frequency change at loci with established links to genetic variance in pigmentation in natural populations. We found that across all spatial and temporal scales, phenotypic patterns were associated with variation at pigmentation-related loci, and the sets of genes we identified in each context were largely nonoverlapping. Therefore, our findings suggest that parallel phenotypic evolution is associated with an unpredictable genomic response, with distinct components of the polygenic architecture shifting across each environmental gradient to produce redundant adaptive patterns.

Significance Statement

Shifts in global climate conditions have heightened our need to understand the dynamics and pace of adaptation in natural populations. In order to anticipate the population-level response to rapidly changing environmental conditions, we need to understand whether trait evolution is predictable over short timescales, and whether the genetic basis of adaptation is shared or distinct across multiple timescales. Here, we explored parallelism in the adaptive response of a complex phenotype, D. melanogaster pigmentation, to shared conditions that varied over multiple spatiotemporal scales. Our results demonstrate that while phenotypic adaptation proceeds as a predictable response to environmental gradients, even over short timescales, the genetic basis of the adaptive response is variable and nuanced across spatial and temporal contexts.

Clinal variation in natural populations of Drosophila melanogaster: An old debate about natural selection and neutral processes

Distinguishing between environmental adaptations and neutral processes poses a challenge in population genetics and evolutionary studies, particularly when phenomena can be explained by both processes. Clines are genotypic or phenotypic characters correlated with environmental variables, because of that correlation, they are used as examples of spatially varying selection. At the same time, many genotypic clines can be explained by demographic history, like isolation by distance or secondary contact zones. Clines have been extensively studied in Drosophila melanogaster, especially in North America and Australia, where they are attributed to both differential selection and various demographic processes. This review explores existing literature supporting this conclusion and suggests new approaches to better understand the influence of these processes on clines. These innovative approaches aim to shed light on the longstanding debate regarding the importance of natural selection versus neutral processes in maintaining variation in natural populations.

A cosmopolitan inversion facilitates seasonal adaptation in overwintering Drosophila

Fluctuations in the strength and direction of natural selection through time are a ubiquitous feature of life on Earth. One evolutionary outcome of such fluctuations is adaptive tracking, wherein populations rapidly adapt from standing genetic variation. In certain circumstances, adaptive tracking can lead to the long-term maintenance of functional polymorphism despite allele frequency change due to selection. Although adaptive tracking is likely a common process, we still have a limited understanding of aspects of its genetic architecture and its strength relative to other evolutionary forces such as drift. Drosophila melanogaster living in temperate regions evolve to track seasonal fluctuations and are an excellent system to tackle these gaps in knowledge. By sequencing orchard populations collected across multiple years, we characterized the genomic signal of seasonal demography and identified that the cosmopolitan inversion In(2L)t facilitates seasonal adaptive tracking and shows molecular footprints of selection. A meta-analysis of phenotypic studies shows that seasonal loci within In(2L)t are associated with behavior, life history, physiology, and morphological traits. We identify candidate loci and experimentally link them to phenotype. Our work contributes to our general understanding of fluctuating selection and highlights the evolutionary outcome and dynamics of contemporary selection on inversions.

Genetic basis and repeatability for desiccation resistance in Drosophila melanogaster (Diptera: Drosophilidae)

Dehydration is a stress factor for organisms inhabiting natural habitats where water is scarce. Thus, it may be expected that species facing arid environments will develop mechanisms that maximize resistance to desiccation. Insects are excellent models for studying the effects of dehydration as well as the mechanisms and processes that prevent water loss since the effect of desiccation is greater due to the higher area/volume ratio than larger animals. Even though physiological and behavioral mechanisms to cope with desiccation are being understood, the genetic basis underlying the mechanisms related to variation in desiccation resistance and the context-dependent effect remain unsolved. Here we analyze the genetic bases of desiccation resistance in Drosophila melanogaster and identify candidate genes that underlie trait variation. Our quantitative genetic analysis of desiccation resistance revealed sexual dimorphism and extensive genetic variation. The phenotype-genotype association analyses (GWAS) identified 71 candidate genes responsible for total phenotypic variation in desiccation resistance. Half of these candidate genes were sex-specific suggesting that the genetic architecture underlying this adaptive trait differs between males and females. Moreover, the public availability of desiccation data analyzed on the same lines but in a different lab allows us to investigate the reliability and repeatability of results obtained in independent screens. Our survey indicates a pervasive micro-environment lab-dependent effect since we did not detect overlap in the sets of genes affecting desiccation resistance identified between labs.

Gene expression differences consistent with water loss reduction underlie desiccation tolerance of natural Drosophila populations

BACKGROUND

Climate change is one of the main factors shaping the distribution and biodiversity of organisms, among others by greatly altering water availability, thus exposing species and ecosystems to harsh desiccation conditions. However, most of the studies so far have focused on the effects of increased temperature. Integrating transcriptomics and physiology is key to advancing our knowledge on how species cope with desiccation stress, and these studies are still best accomplished in model organisms.

RESULTS

Here, we characterized the natural variation of European D. melanogaster populations across climate zones and found that strains from arid regions were similar or more tolerant to desiccation compared with strains from temperate regions. Tolerant and sensitive strains differed not only in their transcriptomic response to stress but also in their basal expression levels. We further showed that gene expression changes in tolerant strains correlated with their physiological response to desiccation stress and with their cuticular hydrocarbon composition, and functionally validated three of the candidate genes identified. Transposable elements, which are known to influence stress response across organisms, were not found to be enriched nearby differentially expressed genes. Finally, we identified several tRNA-derived small RNA fragments that differentially targeted genes in response to desiccation stress.

CONCLUSIONS

Overall, our results showed that basal gene expression differences across individuals should be analyzed if we are to understand the genetic basis of differential stress survival. Moreover, tRNA-derived small RNA fragments appear to be relevant across stress responses and allow for the identification of stress-response genes not detected at the transcriptional level.

Direct observation of adaptive tracking on ecological time scales in Drosophila

Direct observation of evolution in response to natural environmental change can resolve fundamental questions about adaptation, including its pace, temporal dynamics, and underlying phenotypic and genomic architecture. We tracked the evolution of fitness-associated phenotypes and allele frequencies genome-wide in 10 replicate field populations of Drosophila melanogaster over 10 generations from summer to late fall. Adaptation was evident over each sampling interval (one to four generations), with exceptionally rapid phenotypic adaptation and large allele frequency shifts at many independent loci. The direction and basis of the adaptive response shifted repeatedly over time, consistent with the action of strong and rapidly fluctuating selection. Overall, we found clear phenotypic and genomic evidence of adaptive tracking occurring contemporaneously with environmental change, thus demonstrating the temporally dynamic nature of adaptation.

Refining the evolutionary time machine: An assessment of whole genome amplification using single historical Daphnia eggs

Whole genome sequencing is instrumental for the study of genome variation in natural populations, delivering important knowledge on genomic modifications and potential targets of natural selection at the population level. Large dormant eggbanks of aquatic invertebrates such as the keystone herbivore Daphnia, a microcrustacean widespread in freshwater ecosystems, provide detailed sedimentary archives to study genomic processes over centuries. To overcome the problem of limited DNA amounts in single Daphnia dormant eggs, we developed an optimized workflow for whole genome amplification (WGA), yielding sufficient amounts of DNA for downstream whole genome sequencing of individual historical eggs, including polyploid lineages. We compare two WGA kits, applied to recently produced Daphnia magna dormant eggs from laboratory cultures, and to historical dormant eggs of Daphnia pulicaria collected from Arctic lake sediment between 10 and 300 years old. Resulting genome coverage breadth in most samples was ~70%, including those from >100-year-old isolates. Sequence read distribution was highly correlated among samples amplified with the same kit, but less correlated between kits. Despite this, a high percentage of genomic positions with single nucleotide polymorphisms in one or more samples (maximum of 74% between kits, and 97% within kits) were recovered at a depth required for genotyping. As a by-product of sequencing we obtained 100% coverage of the mitochondrial genomes even from the oldest isolates (~300 years). The mitochondrial DNA provides an additional source for evolutionary studies of these populations. We provide an optimized workflow for WGA followed by whole genome sequencing including steps to minimize exogenous DNA.

Genomic Responses to Climate Change: Making the Most of the Drosophila Model

It is pressing to understand how animal populations evolve in response to climate change. We argue that new sequencing technologies and the use of historical samples are opening unprecedented opportunities to investigate genome-wide responses to changing environments. However, there are important challenges in interpreting the emerging findings. First, it is essential to differentiate genetic adaptation from phenotypic plasticity. Second, it is extremely difficult to map genotype, phenotype, and fitness. Third, neutral demographic processes and natural selection affect genetic variation in similar ways. We argue that Drosophila melanogaster, a classical model organism with decades of climate adaptation research, is uniquely suited to overcome most of these challenges. In the near future, long-term time series genome-wide datasets of D. melanogaster natural populations will provide exciting opportunities to study adaptation to recent climate change and will lay the groundwork for related research in non-model systems.

Broad geographic sampling reveals the shared basis and environmental correlates of seasonal adaptation in Drosophila

To advance our understanding of adaptation to temporally varying selection pressures, we identified signatures of seasonal adaptation occurring in parallel among Drosophila melanogaster populations. Specifically, we estimated allele frequencies genome-wide from flies sampled early and late in the growing season from 20 widely dispersed populations. We identified parallel seasonal allele frequency shifts across North America and Europe, demonstrating that seasonal adaptation is a general phenomenon of temperate fly populations. Seasonally fluctuating polymorphisms are enriched in large chromosomal inversions, and we find a broad concordance between seasonal and spatial allele frequency change. The direction of allele frequency change at seasonally variable polymorphisms can be predicted by weather conditions in the weeks prior to sampling, linking the environment and the genomic response to selection. Our results suggest that fluctuating selection is an important evolutionary force affecting patterns of genetic variation in Drosophila.

No water, no mating: Connecting dots from behaviour to pathways

Insects hold considerable ecological and agricultural importance making it vital to understand the factors impacting their reproductive output. Environmental stressors are examples of such factors which have a substantial and significant influence on insect reproductive fitness. Insects are also ectothermic and small in size which makes them even more susceptible to environmental stresses. The present study assesses the consequence of desiccation on the mating latency and copulations duration in tropical Drosophila melanogaster. We tested flies for these reproductive behavioral parameters at varying body water levels and with whole metabolome analysis in order to gain a further understanding of the physiological response to desiccation. Our results showed that the duration of desiccation is positively correlated with mating latency and mating failure, while having no influence on the copulation duration. The metabolomic analysis revealed three biological pathways highly affected by desiccation: starch and sucrose metabolism, galactose metabolism, and phenylalanine, tyrosine and tryptophan biosynthesis. These results are consistent with carbohydrate metabolism providing an energy source in desiccated flies and also suggests that the phenylalanine biosynthesis pathway plays a role in the reproductive fitness of the flies. Desiccation is a common issue with smaller insects, like Drosophila and other tropical insects, and our findings indicate that this lack of ambient water can immediately and drastically affect the insect reproductive behaviour, which becomes more crucial because of unpredictable and dynamic weather conditions.

Allelic polymorphism at foxo contributes to local adaptation in Drosophila melanogaster

The insulin/insulin-like growth factor signalling pathway has been hypothesized as a major determinant of life-history profiles that vary adaptively in natural populations. In Drosophila melanogaster, multiple components of this pathway vary predictably with latitude; this includes foxo, a conserved gene that regulates insulin signalling and has pleiotropic effects on a variety of fitness-associated traits. We hypothesized that allelic variation at foxo contributes to genetic variance for size-related traits that vary adaptively with latitude. We first examined patterns of variation among natural populations along a latitudinal transect in the eastern United States and show that thorax length, wing area, wing loading, and starvation tolerance exhibit significant latitudinal clines for both males and females but that development time does not vary predictably with latitude. We then generated recombinant outbred populations and show that naturally occurring allelic variation at foxo, which exhibits stronger clinality than expected, is associated with the same traits that vary with latitude in the natural populations. Our results suggest that allelic variation at foxo contributes to adaptive patterns of life-history variation in natural populations of this genetic model.

Genomic signatures of thermal adaptation are associated with clinal shifts of life history in a broadly distributed frog

Temperature is a critical driver of ectotherm life‐history strategies, whereby a warmer environment is associated with increased growth, reduced longevity and accelerated senescence. Increasing evidence indicates that thermal adaptation may underlie such life‐history shifts in wild populations. Single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) can help uncover the molecular mechanisms of temperature‐driven variation in growth, longevity and senescence. However, our understanding of these mechanisms is still limited, which reduces our ability to predict the response of non‐model ectotherms to global temperature change.

In this study, we examined the potential role of thermal adaptation in clinal shifts of life‐history traits (i.e. life span, senescence rate and recruitment) in the Columbia spotted frog Rana luteiventris along a broad temperature gradient in the western United States.

We took advantage of extensive capture–recapture datasets of 20,033 marked individuals from eight populations surveyed annually for 14–18 years to examine how mean annual temperature and precipitation influenced demographic parameters (i.e. adult survival, life span, senescence rate, recruitment and population growth). After showing that temperature was the main climatic predictor influencing demography, we used RAD‐seq data (50,829 SNPs and 6,599 putative CNVs) generated for 352 individuals from 31 breeding sites to identify the genomic signatures of thermal adaptation.

Our results showed that temperature was negatively associated with annual adult survival and reproductive life span and positively associated with senescence rate. By contrast, recruitment increased with temperature, promoting the long‐term viability of most populations. These temperature‐dependent demographic changes were associated with strong genomic signatures of thermal adaptation. We identified 148 SNP candidates associated with temperature including three SNPs located within protein‐coding genes regulating resistance to cold and hypoxia, immunity and reproduction in ranids. We also identified 39 CNV candidates (including within 38 transposable elements) for which normalized read depth was associated with temperature.

Our study indicates that both SNPs and structural variants are associated with temperature and could eventually be found to play a functional role in clinal shifts in senescence rate and life‐history strategies in R. luteiventris. These results highlight the potential role of different sources of molecular variation in the response of ectotherms to environmental temperature variation in the context of global warming.

Genetic basis of susceptibility to low-dose paraquat and variation between the sexes in Drosophila melanogaster

Toxicant resistance is a complex trait, affected both by genetics and the environment. Like most complex traits, it can exhibit sexual dimorphism, yet sex is often overlooked as a factor in studies of toxicant resistance. Paraquat, one such toxicant, is a commonly used herbicide and is known to produce mitochondrial oxidative stress, decrease dopaminergic neurons and dopamine (DA) levels, and decrease motor ability. While the main effects of paraquat are well-characterized, less is known about the naturally occurring variation in paraquat susceptibility. The purpose of this study was to map the genes contributing to low-dose paraquat susceptibility in Drosophila melanogaster, and to determine if susceptibility differs between the sexes. One hundred of the Drosophila Genetic Reference Panel (DGRP) lines were scored for susceptibility via climbing ability and used in a genome-wide association study (GWAS). Variation in seventeen genes in females and thirty-five genes in males associated with paraquat susceptibility. Only two candidate genes overlapped between the sexes despite a significant positive correlation between male and female susceptibilities. Many associated polymorphisms had significant interactions with sex, with most having conditionally neutral effects. Conditional neutrality between the sexes probably stems from sex-biased expression which may result from partial resolution of sexual conflict. Candidate genes were verified with RNAi knockdowns, gene expression analyses, and DA quantification. Several of these genes are novel associations with paraquat susceptibility. This research highlights the importance of assessing both sexes when studying toxicant susceptibility.

Eco-genetics of desiccation resistance in Drosophila

Climate change globally perturbs water circulation thereby influencing ecosystems including cultivated land. Both harmful and beneficial species of insects are likely to be vulnerable to such changes in climate. As small animals with a disadvantageous surface area to body mass ratio, they face a risk of desiccation. A number of behavioural, physiological and genetic strategies are deployed to solve these problems during adaptation in various Drosophila species. Over 100 desiccation-related genes have been identified in laboratory and wild populations of the cosmopolitan fruit fly Drosophila melanogaster and its sister species in large-scale and single-gene approaches. These genes are involved in water sensing and homeostasis, and barrier formation and function via the production and composition of surface lipids and via pigmentation. Interestingly, the genetic strategy implemented in a given population appears to be unpredictable. In part, this may be due to different experimental approaches in different studies. The observed variability may also reflect a rich standing genetic variation in Drosophila allowing a quasi-random choice of response strategies through soft-sweep events, although further studies are needed to unravel any underlying principles. These findings underline that D. melanogaster is a robust species well adapted to resist climate change-related desiccation. The rich data obtained in Drosophila research provide a framework to address and understand desiccation resistance in other insects. Through the application of powerful genetic tools in the model organism D. melanogaster, the functions of desiccation-related genes revealed by correlative studies can be tested and the underlying molecular mechanisms of desiccation tolerance understood. The combination of the wealth of available data and its genetic accessibility makes Drosophila an ideal bioindicator. Accumulation of data on desiccation resistance in Drosophila may allow us to create a world map of genetic evolution in response to climate change in an insect genome. Ultimately these efforts may provide guidelines for dealing with the effects of climate-related perturbations on insect population dynamics in the future.

Temperature, rainfall and wind variables underlie environmental adaptation in natural populations of Drosophila melanogaster

While several studies in a diverse set of species have shed light on the genes underlying adaptation, our knowledge on the selective pressures that explain the observed patterns lags behind. Drosophila melanogaster is a valuable organism to study environmental adaptation because this species originated in Southern Africa and has recently expanded worldwide, and also because it has a functionally well-annotated genome. In this study, we aimed to decipher which environmental variables are relevant for adaptation of D. melanogaster natural populations in Europe and North America. We analysed 36 whole-genome pool-seq samples of D. melanogaster natural populations collected in 20 European and 11 North American locations. We used the BayPass software to identify single nucleotide polymorphisms (SNPs) and transposable elements (TEs) showing signature of adaptive differentiation across populations, as well as significant associations with 59 environmental variables related to temperature, rainfall, evaporation, solar radiation, wind, daylight hours, and soil type. We found that in addition to temperature and rainfall, wind related variables are also relevant for D. melanogaster environmental adaptation. Interestingly, 23%-51% of the genes that showed significant associations with environmental variables were not found overly differentiated across populations. In addition to SNPs, we also identified 10 reference transposable element insertions associated with environmental variables. Our results showed that genome-environment association analysis can identify adaptive genetic variants that are undetected by population differentiation analysis while also allowing the identification of candidate environmental drivers of adaptation.

Phenotypic Responses to and Genetic Architecture of Sterility Following Exposure to Sub-Lethal Temperature During Development

Thermal tolerance range, based on temperatures that result in incapacitating effects, influences species’ distributions and has been used to predict species’ response to increasing temperature. Reproductive performance may also be negatively affected at less extreme temperatures, but such sublethal heat-induced sterility has been relatively ignored in studies addressing the potential effects of, and ability of species’ to respond to, predicted climate warming. The few studies examining the link between increased temperature and reproductive performance typically focus on adults, although effects can vary between life history stages. Here we assessed how sublethal heat stress during development impacted subsequent adult fertility and its plasticity, both of which can provide the raw material for evolutionary responses to increased temperature. We quantified phenotypic and genetic variation in fertility of Drosophila melanogaster reared at standardized densities in three temperatures (25, 27, and 29°C) from a set of lines of the Drosophila Genetic Reference Panel (DGRP). We found little phenotypic variation at the two lower temperatures with more variation at the highest temperature and for plasticity. Males were more affected than females. Despite reasonably large broad-sense heritabilities, a genome-wide association study found little evidence for additive genetic variance and no genetic variants were robustly linked with reproductive performance at specific temperatures or for phenotypic plasticity. We compared results on heat-induced male sterility with other DGRP results on relevant fitness traits measured after abiotic stress and found an association between male susceptibility to sterility and male lifespan reduction following oxidative stress. Our results suggest that sublethal stress during development has profound negative consequences on male adult reproduction, but despite phenotypic variation in a population for this response, there is limited evolutionary potential, either through adaptation to a specific developmental temperature or plasticity in response to developmental heat-induced sterility.

Combining Experimental Evolution and Genomics to Understand How Seed Beetles Adapt to a Marginal Host Plant

Genes that affect adaptive traits have been identified, but our knowledge of the genetic basis of adaptation in a more general sense (across multiple traits) remains limited. We combined population-genomic analyses of evolve-and-resequence experiments, genome-wide association mapping of performance traits, and analyses of gene expression to fill this knowledge gap and shed light on the genomics of adaptation to a marginal host (lentil) by the seed beetle Callosobruchus maculatus. Using population-genomic approaches, we detected modest parallelism in allele frequency change across replicate lines during adaptation to lentil. Mapping populations derived from each lentil-adapted line revealed a polygenic basis for two host-specific performance traits (weight and development time), which had low to modest heritabilities. We found less evidence of parallelism in genotype-phenotype associations across these lines than in allele frequency changes during the experiments. Differential gene expression caused by differences in recent evolutionary history exceeded that caused by immediate rearing host. Together, the three genomic datasets suggest that genes affecting traits other than weight and development time are likely to be the main causes of parallel evolution and that detoxification genes (especially cytochrome P450s and beta-glucosidase) could be especially important for colonization of lentil by C. maculatus.

Inference of natural selection from ancient DNA

Evolutionary processes, including selection, can be indirectly inferred based on patterns of genomic variation among contemporary populations or species. However, this often requires unrealistic assumptions of ancestral demography and selective regimes. Sequencing ancient DNA from temporally spaced samples can inform about past selection processes, as time series data allow direct quantification of population parameters collected before, during, and after genetic changes driven by selection. In this Comment and Opinion, we advocate for the inclusion of temporal sampling and the generation of paleogenomic datasets in evolutionary biology, and highlight some of the recent advances that have yet to be broadly applied by evolutionary biologists. In doing so, we consider the expected signatures of balancing, purifying, and positive selection in time series data, and detail how this can advance our understanding of the chronology and tempo of genomic change driven by selection. However, we also recognize the limitations of such data, which can suffer from postmortem damage, fragmentation, low coverage, and typically low sample size. We therefore highlight the many assumptions and considerations associated with analyzing paleogenomic data and the assumptions associated with analytical methods.

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.

The Linked Selection Signature of Rapid Adaptation in Temporal Genomic Data

The majority of empirical population genetic studies have tried to understand the evolutionary processes that have shaped genetic variation in a single sample taken from a present-day population. However, genomic data collected over tens of generations in both natural and laboratory populations are increasingly used to find selected loci underpinning adaptation over these short timescales. Although these studies have been quite successful in detecting selection on large-effect loci, the fitness differences between individuals are often polygenic, such that selection leads to allele frequency changes that are difficult to distinguish from genetic drift. However, one promising signal comes from polygenic selection's effect on neutral sites that become stochastically associated with the genetic backgrounds that lead to fitness differences between individuals. Previous theoretical work has established that the random associations between a neutral allele and heritable fitness backgrounds act to reduce the effective population size experienced by this neutral allele. These associations perturb neutral allele frequency trajectories, creating autocovariance in the allele frequency changes across generations. Here, we show how temporal genomic data allow us to measure the temporal autocovariance in allele frequency changes and characterize the genome-wide impact of polygenic selection. We develop expressions for these temporal autocovariances, showing that their magnitude is determined by the level of additive genetic variation, recombination, and linkage disequilibria in a region. Furthermore, by using analytic expressions for the temporal variances and autocovariances in allele frequency, we demonstrate that one can estimate the additive genetic variation for fitness and the drift-effective population size from temporal genomic data. We also show how the proportion of total variation in allele frequency change due to linked selection can be estimated from temporal data. Overall, we demonstrate that temporal genomic data offer opportunities to identify the role of linked selection on genome-wide diversity over short timescales, and can help bridge population genetic and quantitative genetic studies of adaptation.

Microbiome composition shapes rapid genomic adaptation of Drosophila melanogaster

Population genomic data has revealed patterns of genetic variation associated with adaptation in many taxa. Yet understanding the adaptive process that drives such patterns is challenging; it requires disentangling the ecological agents of selection, determining the relevant timescales over which evolution occurs, and elucidating the genetic architecture of adaptation. Doing so for the adaptation of hosts to their microbiome is of particular interest with growing recognition of the importance and complexity of host-microbe interactions. Here, we track the pace and genomic architecture of adaptation to an experimental microbiome manipulation in replicate populations of Drosophila melanogaster in field mesocosms. Shifts in microbiome composition altered population dynamics and led to divergence between treatments in allele frequencies, with regions showing strong divergence found on all chromosomes. Moreover, at divergent loci previously associated with adaptation across natural populations, we found that the more common allele in fly populations experimentally enriched for a certain microbial group was also more common in natural populations with high relative abundance of that microbial group. These results suggest that microbiomes may be an agent of selection that shapes the pattern and process of adaptation and, more broadly, that variation in a single ecological factor within a complex environment can drive rapid, polygenic adaptation over short timescales.

A clinal polymorphism in the insulin signaling transcription factor foxo contributes to life-history adaptation in Drosophila

A fundamental aim of adaptation genomics is to identify polymorphisms that underpin variation in fitness traits. In Drosophila melanogaster, latitudinal life-history clines exist on multiple continents and make an excellent system for dissecting the genetics of adaptation. We have previously identified numerous clinal single-nucleotide polymorphism in insulin/insulin-like growth factor signaling (IIS), a pathway known from mutant studies to affect life history. However, the effects of natural variants in this pathway remain poorly understood. Here we investigate how two clinal alternative alleles at foxo, a transcriptional effector of IIS, affect fitness components (viability, size, starvation resistance, fat content). We assessed this polymorphism from the North American cline by reconstituting outbred populations, fixed for either the low- or high-latitude allele, from inbred DGRP lines. Because diet and temperature modulate IIS, we phenotyped alleles across two temperatures (18°C, 25°C) and two diets differing in sugar source and content. Consistent with clinal expectations, the high-latitude allele conferred larger body size and reduced wing loading. Alleles also differed in starvation resistance and expression of insulin-like receptor, a transcriptional target of FOXO. Allelic reaction norms were mostly parallel, with few GxE interactions. Together, our results suggest that variation in IIS makes a major contribution to clinal life-history adaptation.

Desiccation-induced changes in recombination rate and crossover interference in Drosophila melanogaster: evidence for fitness-dependent plasticity

Meiotic recombination is evolutionarily ambiguous, as being associated with both benefits and costs to its bearers, with the resultant dependent on a variety of conditions. While existing theoretical models explain the emergence and maintenance of recombination, some of its essential features remain underexplored. Here we focus on one such feature, recombination plasticity, and test whether recombination response to stress is fitness-dependent. We compare desiccation stress effects on recombination rate and crossover interference in chromosome 3 between desiccation-sensitive and desiccation-tolerant Drosophila lines. We show that relative to desiccation-tolerant genotypes, desiccation-sensitive genotypes exhibit a significant segment-specific increase in single- and double-crossover frequencies across the pericentromeric region of chromosome 3. Significant changes (relaxation) in crossover interference were found for the interval pairs flanking the centromere and extending to the left arm of the chromosome. These results indicate that desiccation is a recombinogenic factor and that desiccation-induced changes in both recombination rate and crossover interference are fitness-dependent, with a tendency of less fitted individuals to produce more variable progeny. Such dependence may play an important role in the regulation of genetic variation in populations experiencing environmental challenges.

Dissection of Complex, Fitness-Related Traits in Multiple Drosophila Mapping Populations Offers Insight into the Genetic Control of Stress Resistance

We leverage two complementary Drosophila melanogaster mapping panels to genetically dissect starvation resistance-an important fitness trait. Using >1600 genotypes from the multiparental Drosophila Synthetic Population Resource (DSPR), we map numerous starvation stress QTL that collectively explain a substantial fraction of trait heritability. Mapped QTL effects allowed us to estimate DSPR founder phenotypes, predictions that were correlated with the actual phenotypes of these lines. We observe a modest phenotypic correlation between starvation resistance and triglyceride level, traits that have been linked in previous studies. However, overlap among QTL identified for each trait is low. Since we also show that DSPR strains with extreme starvation phenotypes differ in desiccation resistance and activity level, our data imply multiple physiological mechanisms contribute to starvation variability. We additionally exploited the Drosophila Genetic Reference Panel (DGRP) to identify sequence variants associated with starvation resistance. Consistent with prior work these sites rarely fall within QTL intervals mapped in the DSPR. We were offered a unique opportunity to directly compare association mapping results across laboratories since two other groups previously measured starvation resistance in the DGRP. We found strong phenotypic correlations among studies, but extremely low overlap in the sets of genomewide significant sites. Despite this, our analyses revealed that the most highly associated variants from each study typically showed the same additive effect sign in independent studies, in contrast to otherwise equivalent sets of random variants. This consistency provides evidence for reproducible trait-associated sites in a widely used mapping panel, and highlights the polygenic nature of starvation resistance.