Unique genetic signatures of local adaptation over space and time for diapause, an ecologically relevant complex trait, in Drosophila melanogaster

AMPK regulates HIF-1α to induce pupal diapause in the cotton bollworm, Helicoverpa armigera

Diapause is an adaptive strategy employed by insects to endure adverse environmental conditions and is characterized by reduced metabolic activity, primarily due to a decreased respiratory rate. AMP-activated protein kinase (AMPK) serves as an intracellular energy regulator, modulating energy metabolism in response to metabolic fluctuations. However, its role in pupal diapause of the cotton bollworm, Helicoverpa armigera, remains unclear. In this study, we found that AMPK and its active form, P-AMPK, are highly expressed in diapause-destined pupae. Furthermore, activation of AMPK delayed the development of nondiapause-destined pupae, suggesting a critical role for AMPK in the regulation of pupal diapause in H. armigera. Manipulating AMPK activity in H. armigera epidermal (HaEpi) cells and pupae significantly influenced the expression of hypoxia-inducible factor-1α (HIF-1α), which our laboratory previously reported as a key inducer of pupal diapause through the reduction of mitochondrial activity in H. armigera. Histone deacetylase 4 (HDAC4), a shuttle protein phosphorylated by AMPK which translocates between the cytoplasm and the nucleus, was found to exhibit significantly higher expression in diapause-destined pupal brains compared to their nondiapause counterparts. AMPK in both HaEpi cells and pupae positively regulated the protein levels of P-HDAC4 by binding to the HDAC4 promoter. Additionally, HDAC4 was shown to enhance HIF-1α expression in diapause-destined individuals. HDAC4 binds to and deacetylates heat shock protein 70 (HSP70), and reduced acetylation of HSP70 was found to significantly elevate HIF-1α protein levels. The AMPK-HIF-1α signaling pathway appears to play a pivotal role in reducing mitochondrial activity and facilitating diapause induction in H. armigera pupae.

Trans-Specific Polymorphisms Between Cryptic Daphnia Species Affect Fitness and Behavior

Shared polymorphisms, loci with identical alleles across species, are of unique interest in evolutionary biology as they may represent cases of selection maintaining ancient genetic variation post-speciation, or contemporary selection promoting convergent evolution. In this study, we investigate the abundance of shared polymorphism between two members of the Daphnia pulex species complex. We test whether the presence of shared mutations is consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting or convergent evolution. We analyzed over 2,000 genomes from six taxa in the D. pulex species group and examined the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that North American and European D. pulex diverged over 10 million years ago, yet retained tens of thousands of shared polymorphisms. We suggest that the number of shared polymorphisms between North American and European D. pulex cannot be fully explained by hybridization or incomplete lineage sorting alone. We show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting convergent and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild.

Rapid wing size evolution in African fig flies (Zaprionus indianus) following temperate colonization

Invasive species often encounter novel selective pressures in their invaded range, and understanding their potential for rapid evolution is critical for developing effective management strategies. Zaprionus indianus is an invasive drosophilid native to Africa that reached Florida in 2005 and likely re-establishes temperate North American populations each year. We addressed two evolutionary questions in this system: first, do populations evolve phenotypic changes in the generations immediately following colonization of temperate environments? Second, does Z. indianus evolve directional phenotypic changes along a latitudinal cline? We established isofemale lines from wild collections across space and time and measured twelve ecologically relevant phenotypes, using a reference population as a control. Z. indianus evolved smaller wings following colonization, suggesting early colonizers have larger wings, but smaller wings are favorable after colonization. No other phenotypes changed significantly following colonization or across latitudes, but we did see significant post-colonization changes in principal components of all phenotypes. We documented substantial laboratory evolution and effects of the laboratory environment across multiple phenotypes, emphasizing the importance of controlling for both possibilities when conducting common garden studies. Our results demonstrate the potential for rapid adaptation in Z. indianus, which could contribute to its success and expansion throughout invaded ecosystems.

Limited population structure but signals of recent selection in introduced African Fig Fly (Zaprionus indianus) in North America

Invasive species have devastating consequences for human health, food security, and the environment. Many invasive species adapt to new ecological niches following invasion, but little is known about the early steps of adaptation. Here we examine population genomics of a recently introduced drosophilid in North America, the African Fig Fly, Zaprionus indianus. This species is likely intolerant of subfreezing temperatures and recolonizes temperate environments yearly. We generated a new chromosome-level genome assembly for Z. indianus. Using resequencing of over 200 North American individuals collected over four years in temperate Virginia, plus a single collection from subtropical Florida, we tested for signatures of recolonization, population structure, and adaptation within invasive populations. We show founding populations are sometimes small and contain close genetic relatives, yet temporal population structure and differentiation of populations is mostly absent across recurrent recolonization events. Although we find limited signals of genome-wide spatial or temporal population structure, we identify haplotypes on the X chromosome that are repeatedly differentiated between Virginia and Florida populations. These haplotypes show signatures of natural selection and are not found in African populations. We also find evidence for several large structural polymorphisms segregating within North America populations and show X chromosome evolution in invasive populations is strikingly different from the autosomes. These results show that despite limited population structure, populations may rapidly evolve genetic differences early in an invasion. Further uncovering how these genomic regions influence invasive potential and success in new environments will advance our understanding of how organisms evolve in changing environments.

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.

Anticipating change: The impact of simulated seasonal heterogeneity on heat tolerances along a latitudinal cline

An understanding of thermal limits and variation across geographic regions is central to predicting how any population may respond to global change. Latitudinal clines, in particular, have been used to demonstrate that populations can be locally adapted to their own thermal environment and, as a result, not all populations will be equally impacted by an increase in temperature. But how robust are these signals of thermal adaptation to the other ecological challenges that animals commonly face in the wild? Seasonal changes in population density, food availability, or photoperiod are common ecological challenges that could disrupt patterns of thermal tolerance along a cline if each population differentially used these signals to anticipate future temperatures and adjust their thermal tolerances accordingly. In this study, we aimed to test the robustness of a cline in thermal tolerance to simulated signals of seasonal heterogeneity. Experimental animals were derived from clones of the Australian water flea, Daphnia carinata, sampled from nine distinct populations along a latitudinal transect in Eastern Australia. We then factorially combined summer (18 h light, 6 h dark) and winter (6 h light, 18 h dark) photoperiods with high (5 million algal cells individual-1 day-1) and low (1 million algal cells individual-1 day-1) food availabilities, before performing static heat shock assays to measure thermal tolerance. We found that the thermal tolerances of the clonal populations were sensitive to both measures of seasonal change. In general, higher food availability led to an increase in thermal tolerances, with the magnitude of the increase varying by clone. In contrast, a switch in photoperiod led to rank-order changes in thermal tolerances, with heat resistance increasing for some clones, and decreasing for others. Heat resistance, however, still declined with increasing latitude, irrespective of the manipulation of seasonal signals, with clones from northern populations always showing greater thermal resistance, most likely driven by adaptation to winter thermal conditions. While photoperiod and food availability can clearly shape thermal tolerances for specific populations, they are unlikely to overwhelm overarching signals of thermal adaptation, and thus, observed clines in heat resistance will likely have remained robust to these forms of seasonal heterogeneity.

Across two continents: The genomic basis of environmental adaptation in house mice (Mus musculus domesticus) from the Americas

Replicated clines across environmental gradients can be strong evidence of adaptation. House mice (Mus musculus domesticus) were introduced to the Americas by European colonizers and are now widely distributed from Tierra del Fuego to Alaska. Multiple aspects of climate, such as temperature, vary predictably across latitude in the Americas. Past studies of North American populations across latitudinal gradients provided evidence of environmental adaptation in traits related to body size, metabolism, and behavior and identified candidate genes using selection scans. Here, we investigate genomic signals of environmental adaptation on a second continent, South America, and ask whether there is evidence of parallel adaptation across multiple latitudinal transects in the Americas. We first identified loci across the genome showing signatures of selection related to climatic variation in mice sampled across a latitudinal transect in South America, accounting for neutral population structure. Consistent with previous results, most candidate SNPs were in putatively regulatory regions. Genes that contained the most extreme outliers relate to traits such as body weight or size, metabolism, immunity, fat, eye function, and the cardiovascular system. We then compared these results with the results of analyses of published data from two transects in North America. While most candidate genes were unique to individual transects, we found significant overlap among candidate genes identified independently in the three transects. These genes are diverse, with functions relating to metabolism, immunity, cardiac function, and circadian rhythm, among others. We also found parallel shifts in allele frequency in candidate genes across latitudinal gradients. Finally, combining data from all three transects, we identified several genes associated with variation in body weight. Overall, our results provide strong evidence of shared responses to selection and identify genes that likely underlie recent environmental adaptation in house mice across North and South America.

Balancing selection and the functional effects of shared polymorphism in cryptic Daphnia species

The patterns of genetic variation within and between related taxa represent the genetic history of a species. Shared polymorphisms, loci with identical alleles across species, are of unique interest as they may represent cases of ancient selection maintaining functional variation post-speciation. In this study, we investigate the abundance of shared polymorphism in the Daphnia pulex species complex. We test whether shared mutations are consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting, or convergent evolution. We analyzed over 2,000 genomes from North American and European D. pulex and several outgroup species to examine the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that while North American and European D. pulex diverged over ten million years ago, they retained tens of thousands of shared alleles. We found that the number of shared polymorphisms between North American and European D. pulex cannot be explained by hybridization or incomplete lineage sorting alone. Instead, we show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting young and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild. Ultimately, our findings provide insights into the genetic basis of adaptation and the maintenance of genetic diversity between species.

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.

The limits of stress-tolerance for zooplankton resting stages in freshwater ponds

In seasonal environments, many organisms evolve strategies such as diapause to survive stressful periods. Understanding the link between habitat stability and diapause strategy can help predict a population's survival in a changing world. Indeed, resting stages may be an important way freshwater organisms can survive periods of drought or freezing, and as the frequency and extent of drought or freezing vary strongly among habitats and are predicted to change with climate change, it raises questions about how organisms cope with, and survive, environmental stress. Using Daphnia magna as a model system, we tested the ability of resting stages from different populations to cope with stress during diapause. The combination of elevated temperatures and wet conditions during diapause shows to prevent hatching altogether. In contrast, hatching is relatively higher after a dry and warm diapause, but declines with rising temperatures, while time to hatch increases. Resting stages produced by populations from summer-dry habitats perform slightly, but consistently, better at higher temperatures and dryness, supporting the local adaptation hypothesis. A higher trehalose content in resting eggs from summer-dry habitat might explain such pattern. Considering that temperatures and summer droughts are projected to increase in upcoming years, it is fundamental to know how resting stages resist stressful conditions so as to predict and protect the ecological functioning of freshwater ecosystems.

Genomes from historical Drosophila melanogaster specimens illuminate adaptive and demographic changes across more than 200 years of evolution

The ability to perform genomic sequencing on long-dead organisms is opening new frontiers in evolutionary research. These opportunities are especially notable in the case of museum collections, from which countless documented specimens may now be suitable for genomic analysis-if data of sufficient quality can be obtained. Here, we report 25 newly sequenced genomes from museum specimens of the model organism Drosophila melanogaster, including the oldest extant specimens of this species. By comparing historical samples ranging from the early 1800s to 1933 against modern-day genomes, we document evolution across thousands of generations, including time periods that encompass the species' initial occupation of northern Europe and an era of rapidly increasing human activity. We also find that the Lund, Sweden population underwent local genetic differentiation during the early 1800s to 1933 interval (potentially due to drift in a small population) but then became more similar to other European populations thereafter (potentially due to increased migration). Within each century-scale time period, our temporal sampling allows us to document compelling candidates for recent natural selection. In some cases, we gain insights regarding previously implicated selection candidates, such as ChKov1, for which our inferred timing of selection favors the hypothesis of antiviral resistance over insecticide resistance. Other candidates are novel, such as the circadian-related gene Ahcy, which yields a selection signal that rivals that of the DDT resistance gene Cyp6g1. These insights deepen our understanding of recent evolution in a model system, and highlight the potential of future museomic studies.

Histone methylation regulates reproductive diapause in Drosophila melanogaster

Fluctuating environments threaten fertility and viability. To better match the immediate, local environment, many organisms adopt alternative phenotypic states, a phenomenon called "phenotypic plasticity." Natural populations that predictably encounter fluctuating environments tend to be more plastic than conspecific populations that encounter a constant environment, suggesting that phenotypic plasticity can be adaptive. Despite pervasive evidence of such "adaptive phenotypic plasticity," gene regulatory mechanisms underlying plasticity remains poorly understood. Here we test the hypothesis that environment-dependent phenotypic plasticity is mediated by epigenetic factors. To test this hypothesis, we exploit the adaptive reproductive arrest of Drosophila melanogaster females, called diapause. Using an inbred line from a natural population with high diapause plasticity, we demonstrate that diapause is determined epigenetically: only a subset of genetically identical individuals enter diapause and this diapause plasticity is epigenetically transmitted for at least three generations. Upon screening a suite of epigenetic marks, we discovered that the active histone marks H3K4me3 and H3K36me1 are depleted in diapausing ovaries. Using ovary-specific knockdown of histone mark writers and erasers, we demonstrate that H3K4me3 and H3K36me1 depletion promotes diapause. Given that diapause is highly polygenic, that is, distinct suites of alleles mediate diapause plasticity across distinct genotypes, we also investigated the potential for genetic variation in diapause-determining epigenetic marks. Specifically, we asked if these histone marks were similarly depleted in diapause of a genotypically distinct line. We found evidence of divergence in both the gene expression program and histone mark abundance. This study reveals chromatin determinants of phenotypic plasticity and suggests that these determinants may be genotype-dependent, offering new insight into how organisms may exploit and evolve epigenetic mechanisms to persist in fluctuating environments.

Genome-wide association studies of human and rat BMI converge on synapse, epigenome, and hormone signaling networks

A vexing observation in genome-wide association studies (GWASs) is that parallel analyses in different species may not identify orthologous genes. Here, we demonstrate that cross-species translation of GWASs can be greatly improved by an analysis of co-localization within molecular networks. Using body mass index (BMI) as an example, we show that the genes associated with BMI in humans lack significant agreement with those identified in rats. However, the networks interconnecting these genes show substantial overlap, highlighting common mechanisms including synaptic signaling, epigenetic modification, and hormonal regulation. Genetic perturbations within these networks cause abnormal BMI phenotypes in mice, too, supporting their broad conservation across mammals. Other mechanisms appear species specific, including carbohydrate biosynthesis (humans) and glycerolipid metabolism (rodents). Finally, network co-localization also identifies cross-species convergence for height/body length. This study advances a general paradigm for determining whether and how phenotypes measured in model species recapitulate human biology.

Selection and geography shape male reproductive tract transcriptomes in Drosophila melanogaster

Transcriptome analysis of several animal clades suggests that male reproductive tract gene expression evolves quickly. However, the factors influencing the abundance and distribution of within-species variation, the ultimate source of interspecific divergence, are poorly known. Drosophila melanogaster, an ancestrally African species that has recently spread throughout the world and colonized the Americas in the last roughly 100 years, exhibits phenotypic and genetic latitudinal clines on multiple continents, consistent with a role for spatially varying selection in shaping its biology. Nevertheless, geographic expression variation in the Americas is poorly described, as is its relationship to African expression variation. Here, we investigate these issues through the analysis of two male reproductive tissue transcriptomes [testis and accessory gland (AG)] in samples from Maine (USA), Panama, and Zambia. We find dramatic differences between these tissues in differential expression between Maine and Panama, with the accessory glands exhibiting abundant expression differentiation and the testis exhibiting very little. Latitudinal expression differentiation appears to be influenced by the selection of Panama expression phenotypes. While the testis shows little latitudinal expression differentiation, it exhibits much greater differentiation than the accessory gland in Zambia vs American population comparisons. Expression differentiation for both tissues is non-randomly distributed across the genome on a chromosome arm scale. Interspecific expression divergence between D. melanogaster and D. simulans is discordant with rates of differentiation between D. melanogaster populations. Strongly heterogeneous expression differentiation across tissues and timescales suggests a complex evolutionary process involving major temporal changes in the way selection influences expression evolution in these organs.

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.

Distinct signals of clinal and seasonal allele frequency change at eQTLs in Drosophila melanogaster

Populations of short-lived organisms can respond to spatial and temporal environmental heterogeneity through local adaptation. Local adaptation can be reflected on both phenotypic and genetic levels, and it has been documented in many organisms. Although complex fitness-related phenotypes have been shown to vary across latitudinal clines and seasons in similar ways in Drosophila melanogaster populations, the comparative signals of local adaptation across space and time remain poorly understood. Here, we examined patterns of allele frequency change across a latitudinal cline and between seasons at previously reported expression quantitative trait loci (eQTLs). We divided eQTLs into groups by using differential expression profiles of fly populations collected across latitudinal clines or exposed to different environmental conditions. In general, we find that eQTLs are enriched for clinally varying polymorphisms, and that these eQTLs change in frequency in concordant ways across the cline and in response to starvation and chill-coma. The enrichment of eQTLs among seasonally varying polymorphisms is more subtle, and the direction of allele frequency change at eQTLs appears to be somewhat idiosyncratic. Taken together, we suggest that clinal adaptation at eQTLs is at least partially distinct from seasonal adaptation.

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.

Powerful, efficient QTL mapping in Drosophila melanogaster using bulked phenotyping and pooled sequencing

Despite the value of recombinant inbred lines for the dissection of complex traits, large panels can be difficult to maintain, distribute, and phenotype. An attractive alternative to recombinant inbred lines for many traits leverages selecting phenotypically extreme individuals from a segregating population, and subjecting pools of selected and control individuals to sequencing. Under a bulked or extreme segregant analysis paradigm, genomic regions contributing to trait variation are revealed as frequency differences between pools. Here, we describe such an extreme quantitative trait locus, or extreme quantitative trait loci, mapping strategy that builds on an existing multiparental population, the Drosophila Synthetic Population Resource, and involves phenotyping and genotyping a population derived by mixing hundreds of Drosophila Synthetic Population Resource recombinant inbred lines. Simulations demonstrate that challenging, yet experimentally tractable extreme quantitative trait loci designs (≥4 replicates, ≥5,000 individuals/replicate, and selecting the 5-10% most extreme animals) yield at least the same power as traditional recombinant inbred line-based quantitative trait loci mapping and can localize variants with sub-centimorgan resolution. We empirically demonstrate the effectiveness of the approach using a 4-fold replicated extreme quantitative trait loci experiment that identifies 7 quantitative trait loci for caffeine resistance. Two mapped extreme quantitative trait loci factors replicate loci previously identified in recombinant inbred lines, 6/7 are associated with excellent candidate genes, and RNAi knock-downs support the involvement of 4 genes in the genetic control of trait variation. For many traits of interest to drosophilists, a bulked phenotyping/genotyping extreme quantitative trait loci design has considerable advantages.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

Natural variation in the transcriptional response of Drosophila melanogaster to oxidative stress

Broadly distributed species must cope with diverse and changing environmental conditions, including various forms of stress. Cosmopolitan populations of Drosophila melanogaster are more tolerant to oxidative stress than those from the species’ ancestral range in sub-Saharan Africa, and the degree of tolerance is associated with an insertion/deletion polymorphism in the 3′ untranslated region of the Metallothionein A (MtnA) gene that varies clinally in frequency. We examined oxidative stress tolerance and the transcriptional response to oxidative stress in cosmopolitan and sub-Saharan African populations of D. melanogaster, including paired samples with allelic differences at the MtnA locus. We found that the effect of the MtnA polymorphism on oxidative stress tolerance was dependent on the genomic background, with the deletion allele increasing tolerance only in a northern, temperate population. Genes that were differentially expressed under oxidative stress included MtnA and other metallothioneins, as well as those involved in glutathione metabolism and other genes known to be part of the oxidative stress response or the general stress response. A gene coexpression analysis revealed further genes and pathways that respond to oxidative stress including those involved in additional metabolic processes, autophagy, and apoptosis. There was a significant overlap among the genes induced by oxidative and cold stress, which suggests a shared response pathway to these two stresses. Interestingly, the MtnA deletion was associated with consistent changes in the expression of many genes across all genomic backgrounds, regardless of the expression level of the MtnA gene itself. We hypothesize that this is an indirect effect driven by the loss of microRNA binding sites within the MtnA 3′ untranslated region.

Neuroendocrine Regulation of Reproductive Dormancy in the Fruit Fly Drosophila melanogaster: A Review of Juvenile Hormone-Dependent Regulation

Animals can adjust their physiology, helping them survive and reproduce under a wide range of environmental conditions. One of the strategies to endure unfavorable environmental conditions such as low temperature and limited food supplies is dormancy. In some insect species, this may manifest as reproductive dormancy, which causes their reproductive organs to be severely depleted under conditions unsuitable for reproduction. Reproductive dormancy in insects is induced by a reduction in juvenile hormones synthesized in the corpus allatum (pl. corpora allata; CA) in response to winter-specific environmental cues, such as low temperatures and short-day length. In recent years, significant progress has been made in the study of dormancy-inducing conditions dependent on CA control mechanisms in Drosophila melanogaster. This review summarizes dormancy control mechanisms in D. melanogaster and discusses the implications for future studies of insect dormancy, particularly focusing on juvenile hormone-dependent regulation.

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.

Accurate, ultra-low coverage genome reconstruction and association studies in Hybrid Swarm mapping populations

Genetic association studies seek to uncover the link between genotype and phenotype, and often utilize inbred reference panels as a replicable source of genetic variation. However, inbred reference panels can differ substantially from wild populations in their genotypic distribution, patterns of linkage-disequilibrium, and nucleotide diversity. As a result, associations discovered using inbred reference panels may not reflect the genetic basis of phenotypic variation in natural populations. To address this problem, we evaluated a mapping population design where dozens to hundreds of inbred lines are outbred for few generations, which we call the Hybrid Swarm. The Hybrid Swarm approach has likely remained underutilized relative to pre-sequenced inbred lines due to the costs of genome-wide genotyping. To reduce sequencing costs and make the Hybrid Swarm approach feasible, we developed a computational pipeline that reconstructs accurate whole genomes from ultra-low-coverage (0.05X) sequence data in Hybrid Swarm populations derived from ancestors with phased haplotypes. We evaluate reconstructions using genetic variation from the Drosophila Genetic Reference Panel as well as variation from neutral simulations. We compared the power and precision of Genome-Wide Association Studies using the Hybrid Swarm, inbred lines, recombinant inbred lines (RILs), and highly outbred populations across a range of allele frequencies, effect sizes, and genetic architectures. Our simulations show that these different mapping panels vary in their power and precision, largely depending on the architecture of the trait. The Hybrid Swam and RILs outperform inbred lines for quantitative traits, but not for monogenic ones. Taken together, our results demonstrate the feasibility of the Hybrid Swarm as a cost-effective method of fine-scale genetic mapping.

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.