Insertions and deletions (indels) are a major source of genetic variation within species and may result in functional changes to coding or regulatory sequences. In this study we report that an indel polymorphism in the 3’ untranslated region (UTR) of the metallothionein gene MtnA is associated with gene expression variation in natural populations of Drosophila melanogaster. A derived allele of MtnA with a 49-bp deletion in the 3' UTR segregates at high frequency in populations outside of sub-Saharan Africa. The frequency of the deletion increases with latitude across multiple continents and approaches 100% in northern Europe. Flies with the deletion have more than 4-fold higher MtnA expression than flies with the ancestral sequence. Using reporter gene constructs in transgenic flies, we show that the 3' UTR deletion significantly contributes to the observed expression difference. Population genetic analyses uncovered signatures of a selective sweep in the MtnA region within populations from northern Europe. We also find that the 3’ UTR deletion is associated with increased oxidative stress tolerance. These results suggest that the 3' UTR deletion has been a target of selection for its ability to confer increased levels of MtnA expression in northern European populations, likely due to a local adaptive advantage of increased oxidative stress tolerance. Although molecular variation is abundant in natural populations, understanding how this variation affects organismal phenotypes that are subject to natural selection remains a major challenge in the field of evolutionary genetics. Here we show that a deletion mutation in a noncoding region of the Drosophila melanogaster Metallothionein A gene leads to a significant increase in gene expression and increases survival under oxidative stress. The deletion is in high frequency in three distinct geographic regions: in northern European populations, in northern populations along the east coast of North America, and in southern populations along the east coast of Australia. In northern European populations the deletion shows population genetic signatures of recent positive selection. Thus, we provide evidence for a regulatory polymorphism that underlies local adaptation in natural populations.

During the evolution of heteromorphic sex chromosomes, the sex-specific Y chromosome degenerates, while the X chromosome evolves new mechanisms of regulation. Using bioinformatic and experimental approaches, we investigate the expression of the X chromosome in Drosophila melanogaster. We observe nearly complete X chromosome dosage compensation in male somatic tissues, but not in testis. The X chromosome contains disproportionately fewer genes with high expression in testis than the autosomes, even after accounting for the lack of dosage compensation, which suggests that another mechanism suppresses their expression in the male germline. This is consistent with studies of reporter genes and transposed genes, which find that the same gene has higher expression when autosomal than when X-linked. Using a new reporter gene that is expressed in both testis and somatic tissues, we find that the suppression of X-linked gene expression is limited to genes with high expression in testis and that the extent of the suppression is positively correlated with expression level.

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatio-temporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster . Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterise variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

A particularly well-studied evolutionary model is the vinegar fly Drosophila melanogaster, a cosmopolitan insect of ancestral southern-central African origin. Recent work suggests that it expanded out of Africa ∼9,000 years ago, and spread from the Middle East into Europe ∼1,800 years ago. During its global expansion, this human commensal adapted to novel climate zones and habitats. Despite much work on phenotypic differentiation and adaptation on several continents (especially North America and Australia), typically in the context of latitudinal clines, little is known about phenotypic divergence among European populations. Here, we sought to provide a continent-wide study of phenotypic differentiation among European populations of D. melanogaster. In a consortium-wide phenomics effort, we assayed 16 fitness-related traits on a panel of 173 isofemale lines from 9 European populations, with the majority of traits measured by several groups using semi-standardized protocols. For most fitness-related traits, we found significant differentiation among populations on a continental scale. Despite inevitable differences in assay conditions among labs, the reproducibility and hence robustness of our measurements were overall remarkably good. Several fitness components (e.g., viability, development time) exhibited significant latitudinal or longitudinal clines, and populations differed markedly in multivariate trait structure. Notably, populations experiencing higher humidity/rainfall and lower maximum temperature showed higher viability, fertility, starvation resistance, and lifespan at the expense of lower heat-shock survival, suggesting a pattern of local adaptation. Our results indicate that derived populations of this tropical fly have been shaped by pervasive spatially varying multivariate selection and adaptation to different climates on the European continent.

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome datasets from natural populations of this species have been published over the last 20 years. A major challenge is the integration of these disparate datasets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution and population structure 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 population samples from over 100 locations in >20 countries on four continents based on a combination of 121 unpublished and 150 previously published genomic datasets. Several of these locations have been sampled at different seasons across multiple years. This dataset, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental meta-data. A web-based genome browser and web portal provide easy access to the SNP dataset. 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 dataset. Our resource will enable population geneticists to analyze spatio-temporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.