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

Genetic Architecture of the Thermal Tolerance Landscape in Drosophila melanogaster

Increased environmental temperatures associated with global warming strongly impact natural populations of ectothermic species. Therefore, it is crucial to understand the genetic basis and evolutionary potential of heat tolerance. However, heat tolerance and its genetic components depend on the methodology, making it difficult to predict the adaptive responses to global warming. Here, we measured the knockdown time for 100 lines from the Drosophila Genetic Reference Panel (DGRP) at four different static temperatures, and we estimated their thermal-death-time (TDT) curves, which incorporate the magnitude and the time of exposure to thermal stress, to determine the genetic basis of the thermal tolerance landscape. Through quantitative genetic analyses, the knockdown time showed a significant heritability at different temperatures and that its genetic correlations decreased as temperatures differences increased. Significant genotype-by-sex and genotype-by-environment interactions were noted for heat tolerance. We also discovered genetic variability for the two parameters of TDT: CTmax and thermal sensitivity. Taking advantage of the DGRP, we performed a GWAS and identified multiple variants associated with the TDT parameters, which mapped to genes related to signalling and developmental functions. We performed functional validations for some candidate genes using RNAi, which revealed that genes such as mam, KNCQ, or robo3 affect the knockdown time at a specific temperature but are not associated with the TDT parameters. In conlusion, the thermal tolerance landscape display genetic variation and plastic responses, which may facilitate the adaptation of Drosophila populations to a changing world.

(Limited) Predictability of thermal adaptation in invertebrates

Evolutionary genomic approaches provide powerful tools to understand variation in and evolution of physiological processes. Untargeted genomic or transcriptomic screens can identify functionally annotated candidate genes linked to specific physiological processes, in turn suggesting evolutionary roles for these processes. Such studies often aim to inform modeling of the potential of natural populations to adapt to climate change, but these models are most accurate when evolutionary responses are repeatable, and thus predictable. Here, we synthesize the evolutionary genetic and comparative transcriptomic literature on terrestrial and marine invertebrates to assess whether evolutionary responses to temperature are repeatable within populations, across populations and across species. There is compelling evidence for repeatability, sometimes even across species. However, responses to laboratory selection and geographic variation across thermal gradients appear to be highly idiosyncratic. We also survey whether genetic/transcriptomic studies repeatedly identify candidate genes in three functional groups previously associated with the response to thermal stress: heat shock protein (Hsp) genes, proteolysis genes and immunity genes. Multiple studies across terrestrial and marine species identify candidates included in these gene sets. Yet, each of the gene sets are identified in only a minority of studies. Together, these patterns suggest that there is limited predictability of evolutionary responses to natural selection, including across studies within species. We discuss specific patterns for the candidate gene sets, implications for predictive modeling, and other potential applications of evolutionary genetics in elucidating physiology and gene function. Finally, we discuss limitations of inferences from available evolutionary genetic studies and directions for future research.

Rare variations within the serine/arginine-rich splicing factor PtoRSZ21 modulate stomatal size to determine drought tolerance in Populus

Rare variants contribute significantly to the 'missing heritability' of quantitative traits. The genome-wide characteristics of rare variants and their roles in environmental adaptation of woody plants remain unexplored. Utilizing genome-wide rare variant association study (RVAS), expression quantitative trait loci (eQTL) mapping, genetic transformation, and molecular experiments, we explored the impact of rare variants on stomatal morphology and drought adaptation in Populus. Through comparative analysis of five world-wide Populus species, we observed the influence of mutational bias and adaptive selection on the distribution of rare variants. RVAS identified 75 candidate genes correlated with stomatal size (SS)/stomatal density (SD), and a rare haplotype in the promoter of serine/arginine-rich splicing factor PtoRSZ21 emerged as the foremost association signal governing SS. As a positive regulator of drought tolerance, PtoRSZ21 can recruit the core splicing factor PtoU1-70K to regulate alternative splicing (AS) of PtoATG2b (autophagy-related 2). The rare haplotype PtoRSZ21hap2 weakens binding affinity to PtoMYB61, consequently affecting PtoRSZ21 expression and SS, ultimately resulting in differential distribution of Populus accessions in arid and humid climates. This study enhances the understanding of regulatory mechanisms that underlie AS induced by rare variants and might provide targets for drought-tolerant varieties breeding in Populus.

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.

Single-fly genome assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1 Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

Single-fly assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

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.

Rapid evolutionary diversification of the flamenco locus across simulans clade Drosophila species

Suppression of transposable elements (TEs) is paramount to maintain genomic integrity and organismal fitness. In D. melanogaster, the flamenco locus is a master suppressor of TEs, preventing the mobilization of certain endogenous retrovirus-like TEs from somatic ovarian support cells to the germline. It is transcribed by Pol II as a long (100s of kb), single-stranded, primary transcript, and metabolized into ~24-32 nt Piwi-interacting RNAs (piRNAs) that target active TEs via antisense complementarity. flamenco is thought to operate as a trap, owing to its high content of recent horizontally transferred TEs that are enriched in antisense orientation. Using newly-generated long read genome data, which is critical for accurate assembly of repetitive sequences, we find that flamenco has undergone radical transformations in sequence content and even copy number across simulans clade Drosophilid species. Drosophila simulans flamenco has duplicated and diverged, and neither copy exhibits synteny with D. melanogaster beyond the core promoter. Moreover, flamenco organization is highly variable across D. simulans individuals. Next, we find that D. simulans and D. mauritiana flamenco display signatures of a dual-stranded cluster, with ping-pong signals in the testis and/or embryo. This is accompanied by increased copy numbers of germline TEs, consistent with these regions operating as functional dual-stranded clusters. Overall, the physical and functional diversity of flamenco orthologs is testament to the extremely dynamic consequences of TE arms races on genome organization, not only amongst highly related species, but even amongst individuals.

Unveiling Subtle Geographical Clines: Phenotypic Effects and Dynamics of Circadian Clock Gene Polymorphisms

Our understanding of the gene regulatory network that constitutes the circadian clock has greatly increased in recent decades, notably due to the use of Drosophila as a model system. In contrast, the analysis of natural genetic variation that enables the robust function of the clock under a broad range of environments has developed more slowly. In the current study, we analyzed comprehensive genome sequencing data from wild European populations of Drosophila, which were densely sampled through time and space. We identified hundreds of single nucleotide polymorphisms (SNPs) in nine genes associated with the clock, 276 of which exhibited a latitudinal cline in their allele frequencies. While the effect sizes of these clinal patterns were small, indicating subtle adaptations driven by natural selection, they provided important insights into the genetic dynamics of circadian rhythms in natural populations. We selected nine SNPs in different genes and assessed their impact on circadian and seasonal phenotypes by reconstructing outbred populations fixed for either of the SNP alleles, from inbred DGRP strains. The circadian free-running period of the locomotor activity rhythm was affected by an SNP in doubletime (dbt) and eyes absent (Eya). The SNPs in Clock (Clk), Shaggy (Sgg), period (per), and timeless (tim) affected the acrophase. The alleles of the SNP in Eya conferred different levels of diapause and the chill coma recovery response.

Population admixtures in medaka inferred by multiple arbitrary amplicon sequencing

Cost-effective genotyping can be achieved by sequencing PCR amplicons. Short 3-10 base primers can arbitrarily amplify thousands of loci using only a few primers. To improve the sequencing efficiency of the multiple arbitrary amplicon sequencing (MAAS) approach, we designed new primers and examined their efficiency in sequencing and genotyping. To demonstrate the effectiveness of our method, we applied it to examining the population structure of the small freshwater fish, medaka (Oryzias latipes). We obtained 2987 informative SNVs with no missing genotype calls for 67 individuals from 15 wild populations and three artificial strains. The estimated phylogenic and population genetic structures of the wild populations were consistent with previous studies, corroborating the accuracy of our genotyping method. We also attempted to reconstruct the genetic backgrounds of a commercial orange mutant strain, Himedaka, which has caused a genetic disturbance in wild populations. Our admixture analysis focusing on Himedaka showed that at least two wild populations had genetically been contributed to the nuclear genome of this mutant strain. Our genotyping methods and results will be useful in quantitative assessments of genetic disturbance by this commercially available strain.

Patterns of Population Structure and Introgression Among Recently Differentiated Drosophila melanogaster Populations

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

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.

2b-RAD Genotyping of the Seagrass Cymodocea nodosa Along a Latitudinal Cline Identifies Candidate Genes for Environmental Adaptation

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria—Canary Islands, Faro—Portugal, and Ebro Delta—Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, “regulation of transcription” and “signalling” were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

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.

A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality—encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality—most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.

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.

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

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.

Clinal and seasonal changes are correlated in Drosophila melanogaster natural populations

Spatial and seasonal variations in the environment are ubiquitous. Environmental heterogeneity can affect natural populations and lead to covariation between environment and allele frequencies. Drosophila melanogaster is known to harbor polymorphisms that change both with latitude and seasons. Identifying the role of selection in driving these changes is not trivial, because nonadaptive processes can cause similar patterns. Given the environment changes in similar ways across seasons and along the latitudinal gradient, one promising approach may be to look for parallelism between clinal and seasonal changes. Here, we test whether there is a genome-wide correlation between clinal and seasonal changes, and whether the pattern is consistent with selection. Allele frequency estimates were obtained from pooled samples from seven different locations along the east coast of the United States, and across seasons within Pennsylvania. We show that there is a genome-wide correlation between clinal and seasonal variations, which cannot be explained by linked selection alone. This pattern is stronger in genomic regions with higher functional content, consistent with natural selection. We derive a way to biologically interpret these correlations and show that around 3.7% of the common, autosomal variants could be under parallel seasonal and spatial selection. Our results highlight the contribution of natural selection in driving fluctuations in allele frequencies in natural fly populations and point to a shared genomic basis to climate adaptation that happens over space and time in D. melanogaster.

Population Size, Sex and Purifying Selection: Comparative Genomics of Two Sister Taxa of the Wild Yeast Saccharomyces paradoxus

This study uses population genomic data to estimate demographic and selection parameters in two sister lineages of the wild yeast Saccharomyces paradoxus and compare their evolution. We first estimate nucleotide and recombinational diversities in each of the two lineages to infer their population size and frequency of sex and then analyze the rate of mutation accumulation since divergence from their inferred common ancestor to estimate the generation time and efficacy of selection. We find that one of the lineages has significantly higher silent nucleotide diversity and lower linkage disequilibrium, indicating a larger population with more frequent sexual generations. The same lineage also shows shorter generation time and higher efficacy of purifying selection, the latter consistent with the finding of larger population size and more frequent sex. Similar analyses are also performed on the ancestries of individual strains within lineages and we find significant differences between strains implying variation in rates of mitotic cell divisions. Our sample includes some strains originating in the Chernobyl nuclear-accident exclusion zone, which has been subjected to high levels of radiation for nearly 30 years now. We find no evidence, however, for increased rates of mutation. Finally, there is a positive correlation between rates of mutation accumulation and length of growing period, as measured by latitude of the place of origin of strains. Our study illustrates the power of genomic analyses in estimating population and life history parameters and testing predictions based on population genetic theory.

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.

Cold adaptation drives population genomic divergence in the ecological specialist, Drosophila montana

Detecting signatures of ecological adaptation in comparative genomics is challenging, but analysing population samples with characterised geographic distributions, such as clinal variation, can help identify genes showing covariation with important ecological variation. Here, we analysed patterns of geographic variation in the cold-adapted species Drosophila montana across phenotypes, genotypes and environmental conditions and tested for signatures of cold adaptation in population genomic divergence. We first derived the climatic variables associated with the geographic distribution of 24 populations across two continents to trace the scale of environmental variation experienced by the species, and measured variation in the cold tolerance of the flies of six populations from different geographic contexts. We then performed pooled whole genome sequencing of these six populations, and used Bayesian methods to identify SNPs where genetic differentiation is associated with both climatic variables and the population phenotypic measurements, while controlling for effects of demography and population structure. The top candidate SNPs were enriched on the X and fourth chromosomes, and they also lay near genes implicated in other studies of cold tolerance and population divergence in this species and its close relatives. We conclude that ecological adaptation has contributed to the divergence of D. montana populations throughout the genome and in particular on the X and fourth chromosomes, which also showed highest interpopulation F_ST . This study demonstrates that ecological selection can drive genomic divergence at different scales, from candidate genes to chromosome-wide effects.

Experimental evolution supports signatures of sexual selection in genomic divergence

Comparative genomics has contributed to the growing evidence that sexual selection is an important component of evolutionary divergence and speciation. Divergence by sexual selection is implicated in faster rates of divergence of the X chromosome and of genes thought to underlie sexually selected traits, including genes that are sex biased in expression. However, accurately inferring the relative importance of complex and interacting forms of natural selection, demography, and neutral processes that occurred in the evolutionary past is challenging. Experimental evolution provides an opportunity to apply controlled treatments for multiple generations and examine the consequent genomic divergence. Here, we altered sexual selection intensity, elevating sexual selection in polyandrous lines and eliminating it in monogamous lines, and examined patterns of allele frequency divergence in the genome of Drosophila pseudoobscura after more than 160 generations of experimental evolution. Divergence is not uniform across the genome but concentrated in "islands," many of which contain candidate genes implicated in mating behaviors and other sexually selected phenotypes. These are more often seen on the X chromosome, which also shows greater divergence in F ST than neutral expectations. There are characteristic signatures of selection seen in these regions, with lower diversity on the X chromosome than the autosomes, and differences in diversity on the autosomes between selection regimes. Reduced Tajima's D within some of the divergent regions may imply that selective sweeps have occurred, despite considerable recombination. These changes are associated with both differential gene expression between the lines and sex-biased gene expression within the lines. Our results are very similar to those thought to implicate sexual selection in divergence between species and natural populations, and hence provide experimental support for the likely role of sexual selection in driving such types of genetic divergence, but also illustrate how variable outcomes can be for different genomic regions.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Does geographic variation in thermal tolerance in Daphnia represent trade-offs or conditional neutrality?

Geographic variation in thermal tolerance in Daphnia seems to represent genetic load at the loci specifically responsible for heat tolerance resulting from conditional neutrality. We see no evidence of trade-offs between fitness-related traits at 25 °C vs. 10 °C or between two algal diets across Daphnia magna clones from a variety of locations representing the opposite ends of the distribution of long-term heat tolerance. Likewise, we found no evidence of within-environment trade-offs between heat tolerance and fitness-related traits in any of the environments. Neither short-term and long-term heat tolerance shows any consistent relationship with lipid fluorescence polarization and lipid peroxidation across clones or environments. Pervasive positive correlations between fitness-related traits indicate differences in genetic load rather than trade-off based local adaptation or thermal specialization. For heat tolerance such differences may be caused by either relaxation of stabilizing selection due to lower exposure to high temperature extremes, i.e., conditional neutrality, or by small effective population size followed by the recent range expansion.

The genomics of rapid climatic adaptation and parallel evolution in North American house mice

Parallel changes in genotype and phenotype in response to similar selection pressures in different populations provide compelling evidence of adaptation. House mice (Mus musculus domesticus) have recently colonized North America and are found in a wide range of environments. Here we measure phenotypic and genotypic differentiation among house mice from five populations sampled across 21° of latitude in western North America, and we compare our results to a parallel latitudinal cline in eastern North America. First, we show that mice are genetically differentiated between transects, indicating that they have independently colonized similar environments in eastern and western North America. Next, we find genetically-based differences in body weight and nest building behavior between mice from the ends of the western transect which mirror differences seen in the eastern transect, demonstrating parallel phenotypic change. We then conduct genome-wide scans for selection and a genome-wide association study to identify targets of selection and candidate genes for body weight. We find some genomic signatures that are unique to each transect, indicating population-specific responses to selection. However, there is significant overlap between genes under selection in eastern and western house mouse transects, providing evidence of parallel genetic evolution in response to similar selection pressures across North America.

Population Genomics on the Fly: Recent Advances in Drosophila

Drosophila melanogaster, a small dipteran of African origin, represents one of the best-studied model organisms. Early work in this system has uniquely shed light on the basic principles of genetics and resulted in a versatile collection of genetic tools that allow to uncover mechanistic links between genotype and phenotype. Moreover, given its worldwide distribution in diverse habitats and its moderate genome-size, Drosophila has proven very powerful for population genetics inference and was one of the first eukaryotes whose genome was fully sequenced. In this book chapter, we provide a brief historical overview of research in Drosophila and then focus on recent advances during the genomic era. After describing different types and sources of genomic data, we discuss mechanisms of neutral evolution including the demographic history of Drosophila and the effects of recombination and biased gene conversion. Then, we review recent advances in detecting genome-wide signals of selection, such as soft and hard selective sweeps. We further provide a brief introduction to background selection, selection of noncoding DNA and codon usage and focus on the role of structural variants, such as transposable elements and chromosomal inversions, during the adaptive process. Finally, we discuss how genomic data helps to dissect neutral and adaptive evolutionary mechanisms that shape genetic and phenotypic variation in natural populations along environmental gradients. In summary, this book chapter serves as a starting point to Drosophila population genomics and provides an introduction to the system and an overview to data sources, important population genetic concepts and recent advances in the field.

Adaptation, ancestral variation and gene flow in a 'Sky Island' Drosophila species

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south-western USA separated by large expanses of desert. Using population genomic data from more than 300 wild-caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30-100 KYA) and have recently expanded even further (0.2-2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long-term recurrent selection in these genes. In contrast, we find evidence of long-term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll-regulated antimicrobial peptides.

Demographic analyses of a new sample of haploid genomes from a Swedish population of Drosophila melanogaster

European and African natural populations of Drosophila melanogaster have been the focus of several studies aiming at inferring demographic and adaptive processes based on genetic variation data. However, in these analyses little attention has been given to gene flow between African and European samples. Here we present a dataset consisting of 14 fully sequenced haploid genomes sampled from a natural population from the northern species range (Umeå, Sweden). We co-analyzed this new data with an African population to compare the likelihood of several competing demographic scenarios for European and African populations and show that gene flow improves the fit of demographic models to data.

The Evolution of Phenotypic Plasticity in Response to Temperature Stress

Phenotypic plasticity is the ability of a single genotype to produce different phenotypes in response to environmental variation. The importance of phenotypic plasticity in natural populations and its contribution to phenotypic evolution during rapid environmental change is widely debated. Here, we show that thermal plasticity of gene expression in natural populations is a key component of its adaptation: evolution to novel thermal environments increases ancestral plasticity rather than mean genetic expression. We determined the evolution of plasticity in gene expression by conducting laboratory natural selection on a Drosophila simulans population in hot and cold environments. After more than 60 generations in the hot environment, 325 genes evolved a change in plasticity relative to the natural ancestral population. Plasticity increased in 75% of these genes, which were strongly enriched for several well-defined functional categories (e.g., chitin metabolism, glycolysis, and oxidative phosphorylation). Furthermore, we show that plasticity in gene expression of populations exposed to different temperatures is rather similar across species. We conclude that most of the ancestral plasticity can evolve further in more extreme environments.

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

Organisms living in seasonally variable environments utilize cues such as light and temperature to induce plastic responses, enabling them to exploit favorable seasons and avoid unfavorable ones. Local adapation can result in variation in seasonal responses, but the genetic basis and evolutionary history of this variation remains elusive. Many insects, including Drosophila melanogaster, are able to undergo an arrest of reproductive development (diapause) in response to unfavorable conditions. In D. melanogaster, the ability to diapause is more common in high latitude populations, where flies endure harsher winters, and in the spring, reflecting differential survivorship of overwintering populations. Using a novel hybrid swarm-based genome wide association study, we examined the genetic basis and evolutionary history of ovarian diapause. We exposed outbred females to different temperatures and day lengths, characterized ovarian development for over 2800 flies, and reconstructed their complete, phased genomes. We found that diapause, scored at two different developmental cutoffs, has modest heritability, and we identified hundreds of SNPs associated with each of the two phenotypes. Alleles associated with one of the diapause phenotypes tend to be more common at higher latitudes, but these alleles do not show predictable seasonal variation. The collective signal of many small-effect, clinally varying SNPs can plausibly explain latitudinal variation in diapause seen in North America. Alleles associated with diapause are segregating in Zambia, suggesting that variation in diapause relies on ancestral polymorphisms, and both pro- and anti-diapause alleles have experienced selection in North America. Finally, we utilized outdoor mesocosms to track diapause under natural conditions. We found that hybrid swarms reared outdoors evolved increased propensity for diapause in late fall, whereas indoor control populations experienced no such change. Our results indicate that diapause is a complex, quantitative trait with different evolutionary patterns across time and space.

Phenotypic coupling of sleep and starvation resistance evolves in D. melanogaster

BACKGROUND

One hypothesis for the function of sleep is that it serves as a mechanism to conserve energy. Recent studies have suggested that increased sleep can be an adaptive mechanism to improve survival under food deprivation in Drosophila melanogaster. To test the generality of this hypothesis, we compared sleep and its plastic response to starvation in a temperate and tropical population of Drosophila melanogaster.

RESULTS

We found that flies from the temperate population were more starvation resistant, and hypothesized that they would engage in behaviors that are considered to conserve energy, including increased sleep and reduced movement. Surprisingly, temperate flies slept less and moved more when they were awake compared to tropical flies, both under fed and starved conditions, therefore sleep did not correlate with population-level differences in starvation resistance. In contrast, total sleep and percent change in sleep when starved were strongly positively correlated with starvation resistance within the tropical population, but not within the temperate population. Thus, we observe unexpectedly complex relationships between starvation and sleep that vary both within and across populations. These observations falsify the simple hypothesis of a straightforward relationship between sleep and energy conservation. We also tested the hypothesis that starvation is correlated with metabolic phenotypes by investigating stored lipid and carbohydrate levels, and found that stored metabolites partially contributed towards variation starvation resistance.

CONCLUSIONS

Our findings demonstrate that the function of sleep under starvation can rapidly evolve on short timescales and raise new questions about the physiological correlates of sleep and the extent to which variation in sleep is shaped by natural selection.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal 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 characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Light dependent courtship behavior in Drosophila simulans and D. melanogaster

Differences in courtship signals and perception are well-known among Drosophila species. One such described difference is the dependency on light, and thus presumably vision, for copulation success. Many studies have described a difference in light-dependent copulation success between D. melanogaster and D. simulans, identifying D. simulans as a light-dependent species, and D. melanogaster as a light-independent one. However, many of these studies use assays of varying design and few strains to represent the entire species. Here, we attempt to better characterize this purported difference using 11 strains of each species, paired by collection location, in behavioral assays conducted at two different exposure times. We show that, while there is a species-wide difference in magnitude of light-dependent copulation success, D. melanogaster copulation success is, on average, still impaired in the dark at both exposure times we measured. Additionally, there is significant variation in strain-specific ability to copulate in the dark in both species across two different exposure times. We find that this variation correlates strongly with longitude in D. melanogaster, but not in D. simulans. We hypothesize that differences in species history and demography may explain behavioral variation. Finally, we use courtship assays to show that light-dependent copulation success in one D. simulans strain is driven in part by both males and females. We discuss potential differences in courtship signals and/or signal importance between these species and potential for further comparative studies for functional characterization.

Neuronal Function and Dopamine Signaling Evolve at High Temperature in Drosophila

Neuronal activity is temperature sensitive and affects behavioral traits important for individual fitness, such as locomotion and courtship. Yet, we do not know enough about the evolutionary response of neuronal phenotypes in new temperature environments. Here, we use long-term experimental evolution of Drosophila simulans populations exposed to novel temperature regimes. Here, we demonstrate a direct relationship between thermal selective pressure and the evolution of neuronally expressed molecular and behavioral phenotypes. Several essential neuronal genes evolve lower expression at high temperatures and higher expression at low temperatures, with dopaminergic neurons standing out by displaying the most consistent expression change across independent replicates. We functionally validate the link between evolved gene expression and behavioral changes by pharmacological intervention in the experimentally evolved D. simulans populations as well as by genetically triggered expression changes of key genes in D. melanogaster. As natural temperature clines confirm our results for Drosophila and Anopheles populations, we conclude that neuronal dopamine evolution is a key factor for temperature adaptation.

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.

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.

Population genomics perspectives on convergent adaptation

Convergent adaptation is the independent evolution of similar traits conferring a fitness advantage in two or more lineages. Cases of convergent adaptation inform our ideas about the ecological and molecular basis of adaptation. In judging the degree to which putative cases of convergent adaptation provide an independent replication of the process of adaptation, it is necessary to establish the degree to which the evolutionary change is unexpected under null models and to show that selection has repeatedly, independently driven these changes. Here, we discuss the issues that arise from these questions particularly for closely related populations, where gene flow and standing variation add additional layers of complexity. We outline a conceptual framework to guide intuition as to the extent to which evolutionary change represents the independent gain of information owing to selection and show that this is a measure of how surprised we should be by convergence. Additionally, we summarize the ways population and quantitative genetics and genomics may help us address questions related to convergent adaptation, as well as open new questions and avenues of research. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Genomic divergence and adaptive convergence in Drosophila simulans from Evolution Canyon, Israel

Biodiversity refugia formed by unique features of the Mediterranean arid landscape, such as the dramatic ecological contrast of "Evolution Canyon," provide a natural laboratory in which local adaptations to divergent microclimate conditions can be investigated. Significant insights have been provided by studies of Drosophila melanogaster diversifying along the thermal gradient in Evolution Canyon, but a comparative framework to survey adaptive convergence across sister species at the site has been lacking. To fill this void, we present an analysis of genomic polymorphism and evolutionary divergence of Drosophila simulans, a close relative of Drosophila melanogaster with which it co-occurs on both slopes of the canyon. Our results show even deeper interslope divergence in D. simulans than in D. melanogaster, with extensive signatures of selective sweeps present in flies from both slopes but enhanced in the population from the hotter and drier south-facing slope. Interslope divergence was enriched for genes related to electrochemical balance and transmembrane transport, likely in response to increased selection for dehydration resistance on the hotter slope. Both species shared genomic regions that underwent major selective sweeps, but the overall level of adaptive convergence was low, demonstrating no shortage of alternative genomic solutions to cope with the challenges of the microclimate contrast. Mobile elements were a major source of genetic polymorphism and divergence, affecting all parts of the genome, including coding sequences of mating behavior-related genes.

Interactions Between Biotic and Abiotic Factors Affect Survival in Overwintering Drosophila suzukii (Diptera: Drosophilidae)

Drosophila suzukii Matsumura is an invasive species affecting berry crops and cherries throughout North America, South America, and Europe. Previous research suggests that in temperate climates, the overwintering success of D. suzukii is likely dependent on access to food, shelter, and adequate cold hardening. We performed a multi-state study under field conditions for two winters to determine whether D. suzukii sex, phenotype (summer-morphotype, winter-morphotype), and life stage (adults, pupae) affected survival over time while recording naturally-occurring spatial and temporal variation in temperature. Access to food was provided and the flies were buried under leaf litter. Baited traps were deployed to determine whether local populations of D. suzukii were active throughout the winter season. The duration of exposure, mean daily temperature, and cumulative time below freezing significantly affected survival. Below freezing, D. suzukii survival was significantly reduced, particularly in northern locations. In contrast, we observed sustained survival up to 10 wk in southern locations among adults and pupae. Biotic factors also significantly affected survival outcomes: female survival was greater than male survival, winter-morphotype survival was greater than summer-morphotype survival, and adult survival was greater than pupal survival. In the north, wild D. suzukii were captured only in early winter, while in the south they were found throughout the winter. These data suggest that although adult D. suzukii may overwinter in sheltered microclimates, this ability may be limited in regions where the ground temperature, or site of overwintering, falls below freezing for extended durations.

Pervasive Linked Selection and Intermediate-Frequency Alleles Are Implicated in an Evolve-and-Resequencing Experiment of Drosophila simulans

We develop analytical and simulation tools for evolve-and-resequencing experiments and apply them to a new study of rapid evolution in Drosophila simulans Likelihood test statistics applied to pooled population sequencing data suggest parallel evolution of 138 SNPs across the genome. This number is reduced by orders of magnitude from previous studies (thousands or tens of thousands), owing to differences in both experimental design and statistical analysis. Whole genome simulations calibrated from Drosophila genetic data sets indicate that major features of the genome-wide response could be explained by as few as 30 loci under strong directional selection with a corresponding hitchhiking effect. Smaller effect loci are likely also responding, but are below the detection limit of the experiment. Finally, SNPs showing strong parallel evolution in the experiment are intermediate in frequency in the natural population (usually 30-70%) indicative of balancing selection in nature. These loci also exhibit elevated differentiation among natural populations of D. simulans, suggesting environmental heterogeneity as a potential balancing mechanism.

Population Genetic and Functional Analysis of a cis-Regulatory Polymorphism in the DrosophilamelanogasterMetallothionein A gene

Although gene expression can vary extensively within and among populations, the genetic basis of this variation and the evolutionary forces that maintain it are largely unknown. In Drosophilamelanogaster, a 49-bp insertion/deletion (indel) polymorphism in the Metallothionein A (MtnA) gene is associated with variation in MtnA expression and oxidative stress tolerance. To better understand the functional and evolutionary significance of this polymorphism, we investigated it in several worldwide populations. In a German population, the deletion was present at a high and stable frequency over multiple seasons and years, and was associated with increased MtnA expression. There was, however, no evidence that the polymorphism was maintained by overdominant, seasonally fluctuating, or sexually antagonistic selection. The deletion was rare in a population from the species' ancestral range in sub-Saharan Africa and is likely the result of non-African admixture, suggesting that it spread to high frequency following the species' out-of-Africa expansion. Using data from a North American population, we found that the deletion was associated with MtnA expression and tolerance to oxidative stress induced by menadione sodium bisulfite. Our results are consistent with the deletion being selectively favored in temperate populations due to the increased MtnA expression and oxidative stress tolerance that it confers.

Evolutionary insights from large scale resequencing datasets in Drosophila melanogaster

Drosophila melanogaster has long been used as an evolutionary model system. Its small genome size, well-annotated genome, and ease of sampling, also makes it a choice species for genome resequencing studies. Hundreds of genomic samples from populations worldwide are available and are currently being used to tackle a wide range of evolutionary questions. In this review, we focused on three insights that have increased our understanding of the evolutionary history of this species, and that have implications for the study of evolutionary processes in other species as well. Because of technical limitations, most of the studies so far have focused on SNP variants. However, long-read sequencing techniques should allow us in the near future to include other type of genomic variants that also influence genome evolution.

Functional Analysis of a Putative Target of Spatially Varying Selection in the Menin1 Gene of Drosophila melanogaster

While significant effort has been devoted to investigating the potential influence of spatially varying selection on genomic variation, relatively little effort has been devoted to experimental analysis of putative variants or genes experiencing such selection. Previous population genetic work identified an amino acid polymorphism in the Mnn1 gene as one of the most strongly latitudinally differentiated SNPs in the genome of Drosophila melanogaster in the United States and Australia. Here we report the results of our transgenic analysis of this amino acid polymorphism. Genotypes carrying alternative Mnn1 alleles differed in multiple phenotypes in a direction generally consistent with phenotypic differences previously observed along latitudinal clines. These results support inferences from earlier population genomic work that this variant influences fitness, and support the idea that the alleles exhibiting clines may be likely to have pleiotropic effects that are correlated along the axes favored by natural selection.

Spatially varying cis-regulatory divergence in Drosophila embryos elucidates cis-regulatory logic

Spatial patterning of gene expression is a key process in development, yet how it evolves is still poorly understood. Both cis- and trans-acting changes could participate in complex interactions, so to isolate the cis-regulatory component of patterning evolution, we measured allele-specific spatial gene expression patterns in D. melanogaster × simulans hybrid embryos. RNA-seq of cryo-sectioned slices revealed 66 genes with strong spatially varying allele-specific expression. We found that hunchback, a major regulator of developmental patterning, had reduced expression of the D. simulans allele specifically in the anterior tip of hybrid embryos. Mathematical modeling of hunchback cis-regulation suggested a candidate transcription factor binding site variant, which we verified as causal using CRISPR-Cas9 genome editing. In sum, even comparing morphologically near-identical species we identified surprisingly extensive spatial variation in gene expression, suggesting not only that development is robust to many such changes, but also that natural selection may have ample raw material for evolving new body plans via changes in spatial patterning.