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.

Rapid seasonal evolution in innate immunity of wild Drosophila melanogaster.

Understanding the rate of evolutionary change and the genetic architecture that facilitates rapid adaptation is a current challenge in evolutionary biology. Comparative studies show that genes with immune function are among the most rapidly evolving genes across a range of taxa. Here, we use immune defence in natural populations of Drosophila melanogaster to understand the rate of evolution in natural populations and the genetics underlying rapid change. We probed the immune system using the natural pathogens Enterococcus faecalis and Providencia rettgeri to measure post-infection survival and bacterial load of wild D. melanogaster populations collected across seasonal time along a latitudinal transect along eastern North America (Massachusetts, Pennsylvania and Virginia). There are pronounced and repeatable changes in the immune response over the approximately 10 generations between spring and autumn collections, with a significant but less distinct difference observed among geographical locations. Genes with known immune function are not enriched among alleles that cycle with seasonal time, but the immune function of a subset of seasonally cycling alleles in immune genes was tested using reconstructed outbred populations. We find that flies containing seasonal alleles in Thioester-containing protein 3 (Tep3) have different functional responses to infection and that epistatic interactions among seasonal Tep3 and Drosomycin-like 6 (Dro6) alleles underlie the immune phenotypes observed in natural populations. This rapid, cyclic response to seasonal environmental pressure broadens our understanding of the complex ecological and genetic interactions determining the evolution of immune defence in natural populations.

Global Transcriptional Profiling of Diapause and Climatic Adaptation in Drosophila melanogaster.

Wild populations of the model organism Drosophila melanogaster experience highly heterogeneous environments over broad geographical ranges as well as over seasonal and annual timescales. Diapause is a primary adaptation to environmental heterogeneity, and in D. melanogaster the propensity to enter diapause varies predictably with latitude and season. Here we performed global transcriptomic profiling of naturally occurring variation in diapause expression elicited by short day photoperiod and moderately low temperature in two tissue types associated with neuroendocrine and endocrine signaling, heads, and ovaries. We show that diapause in D. melanogaster is an actively regulated phenotype at the transcriptional level, suggesting that diapause is not a simple physiological or reproductive quiescence. Differentially expressed genes and pathways are highly distinct in heads and ovaries, demonstrating that the diapause response is not uniform throughout the soma and suggesting that it may be comprised of functional modules associated with specific tissues. Genes downregulated in heads of diapausing flies are significantly enriched for clinally varying single nucleotide polymorphism (SNPs) and seasonally oscillating SNPs, consistent with the hypothesis that diapause is a driving phenotype of climatic adaptation. We also show that chromosome location-based coregulation of gene expression is present in the transcriptional regulation of diapause. Taken together, these results demonstrate that diapause is a complex phenotype actively regulated in multiple tissues, and support the hypothesis that natural variation in diapause propensity underlies adaptation to spatially and temporally varying selective pressures.

Genomic evidence of rapid and stable adaptive oscillations over seasonal time scales in Drosophila.

In many species, genomic data have revealed pervasive adaptive evolution indicated by the fixation of beneficial alleles. However, when selection pressures are highly variable along a species' range or through time adaptive alleles may persist at intermediate frequencies for long periods. So called "balanced polymorphisms" have long been understood to be an important component of standing genetic variation, yet direct evidence of the strength of balancing selection and the stability and prevalence of balanced polymorphisms has remained elusive. We hypothesized that environmental fluctuations among seasons in a North American orchard would impose temporally variable selection on Drosophila melanogaster that would drive repeatable adaptive oscillations at balanced polymorphisms. We identified hundreds of polymorphisms whose frequency oscillates among seasons and argue that these loci are subject to strong, temporally variable selection. We show that these polymorphisms respond to acute and persistent changes in climate and are associated in predictable ways with seasonally variable phenotypes. In addition, our results suggest that adaptively oscillating polymorphisms are likely millions of years old, with some possibly predating the divergence between D. melanogaster and D. simulans. Taken together, our results are consistent with a model of balancing selection wherein rapid temporal fluctuations in climate over generational time promotes adaptive genetic diversity at loci underlying polygenic variation in fitness related phenotypes.

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster.

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.

A small system--high-resolution study of metabolic adaptation in the central metabolic pathway to temperate climates in Drosophila melanogaster.

In this article, we couple the geographic variation in 127 single-nucleotide polymorphism (SNP) frequencies in genes of 46 enzymes of central metabolism with their associated cis-expression variation to predict latitudinal or climatic-driven gene expression changes in the metabolic architecture of Drosophila melanogaster. Forty-two percent of the SNPs in 65% of the genes show statistically significant clines in frequency with latitude across the 20 local population samples collected from southern Florida to Ontario. A number of SNPs in the screened genes are also associated with significant expression variation within the Raleigh population from North Carolina. A principal component analysis of the full variance-covariance matrix of latitudinal changes in SNP-associated standardized gene expression allows us to identify those major genes in the pathway and its associated branches that are likely targets of natural selection. When embedded in a central metabolic context, we show that these apparent targets are concentrated in the genes of the upper glycolytic pathway and pentose shunt, those controlling glycerol shuttle activity, and finally those enzymes associated with the utilization of glutamate and pyruvate. These metabolites possess high connectivity and thus may be the points where flux balance can be best shifted. We also propose that these points are conserved points associated with coupling energy homeostasis and energy sensing in mammals. We speculate that the modulation of gene expression at specific points in central metabolism that are associated with shifting flux balance or possibly energy-state sensing plays a role in adaptation to climatic variation.

Variation in Drosophila melanogaster central metabolic genes appears driven by natural selection both within and between populations.

In this report, we examine the hypothesis that the drivers of latitudinal selection observed in the eastern US Drosophila melanogaster populations are reiterated within seasons in a temperate orchard population in Pennsylvania, USA. Specifically, we ask whether alleles that are apparently favoured in northern populations are also favoured early in the spring, and decrease in frequency from the spring to autumn with the population expansion. We use SNP data collected for 46 metabolic genes and 128 SNPs representing the central metabolic pathway and examine for the aggregate SNP allele frequencies whether the association of allele change with latitude and that with increasing days of spring-autumn season are reversed. Testing by random permutation, we observe a highly significant negative correlation between these associations that is consistent with this expectation. This correlation is stronger when we confine our analysis to only those alleles that show significant latitudinal changes. This pattern is not caused by association with chromosomal inversions. When data are resampled using SNPs for amino acid change the relationship is not significant but is supported when SNPs associated with cis-expression are only considered. Our results suggest that climate factors driving latitudinal molecular variation in a metabolic pathway are related to those operating on a seasonal level within populations.

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.

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.

On the Long-term Stability of Clines in Some Metabolic Genes in Drosophila melanogaster.

Very little information exists for long-term changes in genetic variation in natural populations. Here we take the unique opportunity to compare a set of data for SNPs in 15 metabolic genes from eastern US collections of Drosophila melanogaster that span a large latitudinal range and represent two collections separated by 12 to 13 years. We also expand this to a 22-year interval for the Adh gene and approximately 30 years for the G6pd and Pgd genes. During these intervals, five genes showed a statistically significant change in average SNP allele frequency corrected for latitude. While much remains unchanged, we see five genes where latitudinal clines have been lost or gained and two where the slope significantly changes. The long-term frequency shift towards a southern favored Adh S allele reported in Australia populations is not observed in the eastern US over a period of 21 years. There is no general pattern of southern-favored or northern-favored alleles increasing in frequency across the genes. This observation points to the fluid nature of some allelic variation over this time period and the action of selective responses or migration that may be more regional than uniformly imposed across the cline.

A highly pleiotropic amino acid polymorphism in the Drosophila insulin receptor contributes to life-history adaptation.

Finding the specific nucleotides that underlie adaptive variation is a major goal in evolutionary biology, but polygenic traits pose a challenge because the complex genotype-phenotype relationship can obscure the effects of individual alleles. However, natural selection working in large wild populations can shift allele frequencies and indicate functional regions of the genome. Previously, we showed that the two most common alleles of a complex amino acid insertion-deletion polymorphism in the Drosophila insulin receptor show independent, parallel clines in frequency across the North American and Australian continents. Here, we report that the cline is stable over at least a five-year period and that the polymorphism also demonstrates temporal shifts in allele frequency concurrent with seasonal change. We tested the alleles for effects on levels of insulin signaling, fecundity, development time, body size, stress tolerance, and life span. We find that the alleles are associated with predictable differences in these traits, consistent with patterns of Drosophila life-history variation across geography that likely reflect adaptation to the heterogeneous climatic environment. These results implicate insulin signaling as a major mediator of life-history adaptation in Drosophila, and suggest that life-history trade-offs can be explained by extensive pleiotropy at a single locus.

The intensity of selection acting on the couch potato gene--spatial-temporal variation in a diapause cline.

Cosmopolitan populations of Drosophila melanogaster have co-opted a form of reproductive diapause to overwinter in northern populations. Polymorphism in the couch potato gene has been implicated in genetic variation for this diapause trait. Using a collection of 20 populations from Florida to Canada and 11 collections from 3 years in a Pennsylvania orchard, we estimated the allele frequencies for 15 single nucleotide polymorphisms (SNPs) in the couch potato gene. These include the specific polymorphism associated with diapause inducability. We find that the SNP polymorphism, 48034(A/T), is correlated with latitude and its frequencies are predicted by the incidence of diapause trait. We find that the clinal patterns for cpo SNPs sampled in 1997 are similar to the same SNPs sampled in 2009-2010. SNPs that show apparent associations with cpo expression are also clinal with the low-expression allele increasing in frequency, as would be predicted from functional knockout studies of cpo. Finally, we see a significant pattern where the frequency of the diapause-causing allele drops in frequency during the summer season, consistent with the drop in the incidence of the diapause trait. The selection required to drive this response is large, roughly 24% to 59% per generation depending on the degree of dominance.

What determines adult cognitive skills? Influences of pre-school, school, and post-school experiences in Guatemala.

Most empirical investigations of the effects of cognitive skills assume that they are produced by schooling. Drawing on longitudinal data to estimate production functions for adult verbal and nonverbal cognitive skills, we find that: (1) School attainment has a significant and substantial effect on adult verbal cognitive skills but not on adult nonverbal cognitive skills; and (2) Pre-school and post-school experiences also have substantial positive significant effects on adult cognitive skills. Pre-school experiences captured by height for age at 6 years substantially and significantly increase adult nonverbal cognitive skills, even after controlling for school attainment. Post-school tenure in skilled jobs has significant positive effects on both types of cognitive skills. The findings (1) reinforce the importance of early life investments; (2) support the importance of childhood nutrition ("Flynn effect") and work complexity in explaining increases in nonverbal cognitive skills; (3) call into question interpretations of studies reporting productivity impacts of cognitive skills that do not control for endogeneity; and (4) point to limitations in using adult school attainment alone to represent human capital.

Hunchback-like protein is expressed in cleavage blastomeres, gastrula epithelium, and ciliary structures in gastropods.

We report the expression of Hunchback (Hb)-like protein during embryonic and larval development in two caenogastropods, Crepidula fornicata and Ilyanassa obsoleta. During the cleavage stages of these species, Hb-like protein is uniformly expressed in micromere and macromere nuclei. At gastrulation, gastropod Hb-like protein is expressed in the surface epithelium that undergoes epiboly. During organogenesis, gastropod Hb-like protein is expressed in the developing ciliated structures associated with feeding and locomotion. We find no detectable gradient or regionalization of Hb-like protein in gastropod embryos or larvae that resembles the graded Hb pattern of expression observed in dipteran insect embryos. Rather we found that the spatiotemporal expression profile of gastropod Hb-like protein is nearly identical to the Hb-like patterns obtained from the polychaete Capitella sp. I and is highly similar to those reported for clitellate annelids. Based upon the comparative data collected from both ecdysozoans and lophotrochozoan lineages, our results support the hypothesis that the role of Hb in anteroposterior patterning is a derived trait specific to arthropods, and that the ancestral function of lophotrochozoan Hb-like protein played a role in the formation of the cleavage-stage blastomeres and the gastrula epithelium and in structures associated with larval feeding and locomotion.