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

Rapid evolutionary change, constraints and the maintenance of polymorphism in natural populations of Drosophila melanogaster

Allele frequencies can shift rapidly within natural populations. Under certain conditions, repeated rapid allele frequency shifts can lead to the long-term maintenance of polymorphism. In recent years, studies of the model insect Drosophila melanogaster have suggested that this phenomenon is more common than previously believed and is often driven by some form of balancing selection, such as temporally fluctuating or sexually antagonistic selection. Here we discuss some of the general insights into rapid evolutionary change revealed by large-scale population genomic studies, as well as the functional and mechanistic causes of rapid adaptation uncovered by single-gene studies. As an example of the latter, we consider a regulatory polymorphism of the D. melanogaster fezzik gene. Polymorphism at this site has been maintained at intermediate frequency over an extended period of time. Regular observations from a single population over a period of 7 years revealed significant differences in the frequency of the derived allele and its variance across collections between the sexes. These patterns are highly unlikely to arise from genetic drift alone or from the action of sexually antagonistic or temporally fluctuating selection individually. Instead, the joint action of sexually antagonistic and temporally fluctuating selection can best explain the observed rapid and repeated allele frequency shifts. Temporal studies such as those reviewed here further our understanding of how rapid changes in selection can lead to the long-term maintenance of polymorphism as well as improve our knowledge of the forces driving and limiting adaptation in nature.

From sub-Saharan Africa to China: Evolutionary history and adaptation of Drosophila melanogaster revealed by population genomics

Drosophila melanogaster is a widely used model organism for studying environmental adaptation. However, the genetic diversity of populations in Asia is poorly understood, leaving a notable gap in our knowledge of the global evolution and adaptation of this species. We sequenced genomes of 292 D. melanogaster strains from various ecological settings in China and analyzed them along with previously published genome sequences. We have identified six global genetic ancestry groups, despite the presence of widespread genetic admixture. The strains from China represent a unique ancestry group, although detectable differentiation exists among populations within China. We deciphered the global migration and demography of D. melanogaster, and identified widespread signals of adaptation, including genetic changes in response to insecticides. We validated the effects of insecticide resistance variants using population cage trials and deep sequencing. This work highlights the importance of population genomics in understanding the genetic underpinnings of adaptation, an effort that is particularly relevant given the deterioration of ecosystems.

The pseudoentropy of allele frequency trajectories, the persistence of variation, and the effective population size

To concisely describe how genetic variation, at individual loci or across whole genomes, changes over time, and to follow transitory allelic changes, we introduce a quantity related to entropy, that we term pseudoentropy. This quantity emerges in a diffusion analysis of the mean time a mutation segregates in a population. For a neutral locus with an arbitrary number of alleles, the mean time of segregation is generally proportional to the pseudoentropy of initial allele frequencies. After the initial time point, pseudoentropy generally decreases, but other behaviours are possible, depending on the genetic diversity and selective forces present. For a biallelic locus, pseudoentropy and entropy coincide, but they are distinct quantities with more than two alleles. Thus for populations with multiple biallelic loci, the language of entropy suffices. Then entropy, combined across loci, serves as a concise description of genetic variation. We used individual based simulations to explore how this entropy behaves under different evolutionary scenarios. In agreement with predictions, the entropy associated with unlinked neutral loci decreases over time. However, deviations from free recombination and neutrality have clear and informative effects on the entropy's behaviour over time. Analysis of publicly available data of a natural D. melanogaster population, that had been sampled over seven years, using a sliding-window approach, yielded considerable variation in entropy trajectories of different genomic regions. These mostly follow a pattern that suggests a substantial effective population size and a limited effect of positive selection on genome-wide diversity over short time scales.

The utility of Drosophila melanogaster as a fungal infection model

Invasive fungal diseases have profound effects upon human health and are on increase globally. The World Health Organization (WHO) in 2022 published the fungal priority list calling for improved public health interventions and advance research. Drosophila melanogaster presents an excellent model system to dissect host-pathogen interactions and has been proved valuable to study immunopathogenesis of fungal diseases. In this review we highlight the recent advances in fungal-Drosophila interplay with an emphasis on the recently published WHO's fungal priority list and we focus on available tools and technologies.

Cis-regulatory polymorphism at fiz ecdysone oxidase contributes to polygenic evolutionary response to malnutrition in Drosophila

We investigate the contribution of a candidate gene, fiz (fezzik), to complex polygenic adaptation to juvenile malnutrition in Drosophila melanogaster. Experimental populations maintained for >250 generations of experimental evolution to a nutritionally poor larval diet (Selected populations) evolved several-fold lower fiz expression compared to unselected Control populations. Here we show that this divergence in fiz expression is mediated by a cis-regulatory polymorphism. This polymorphism, originally sampled from a natural population in Switzerland, is distinct from a second cis-regulatory SNP previously identified in non-African D. melanogaster populations, implying that two independent cis-regulatory variants promoting high fiz expression segregate in non-African populations. Enzymatic analyses of Fiz protein expressed in E. coli demonstrate that it has ecdysone oxidase activity acting on both ecdysone and 20-hydroxyecdysone. Four of five fiz paralogs annotated to ecdysteroid metabolism also show reduced expression in Selected larvae, implying that malnutrition-driven selection favored general downregulation of ecdysone oxidases. Finally, as an independent test of the role of fiz in poor diet adaptation, we show that fiz knockdown by RNAi results in faster larval growth on the poor diet, but at the cost of greatly reduced survival. These results imply that downregulation of fiz in Selected populations was favored by selection on the nutritionally poor diet because of its role in suppressing growth in response to nutrient shortage. However, they suggest that fiz downregulation is only adaptive in combination with other changes evolved by Selected populations, which ensure that the organism can sustain the faster growth promoted by fiz downregulation.

Genetic basis of variation in thermal developmental plasticity for Drosophila melanogaster body pigmentation

Seasonal differences in insect pigmentation are attributed to the influence of ambient temperature on pigmentation development. This thermal plasticity is adaptive and heritable, and thereby capable of evolving. However, the specific genes contributing to the variation in plasticity that can drive its evolution remain largely unknown. To address this, we analysed pigmentation and pigmentation plasticity in Drosophila melanogaster. We measured two components of pigmentation in the thorax and abdomen: overall darkness and the proportion of length covered by darker pattern elements (a trident in the thorax and bands in the abdomen) in females from two developmental temperatures (17 or 28°C) and 191 genotypes. Using a GWAS approach to identify the genetic basis of variation in pigmentation and its response to temperature, we identified numerous dispersed QTLs, including some mapping to melanogenesis genes (yellow, ebony, and tan). Remarkably, we observed limited overlap between QTLs for variation within specific temperatures and those influencing thermal plasticity, as well as minimal overlap between plasticity QTLs across pigmentation components and across body parts. For most traits, consistent with selection favouring the retention of plasticity, we found that lower plasticity alleles were often at lower frequencies. The functional analysis of selected candidate QTLs and pigmentation genes largely confirmed their contributions to variation in pigmentation and/or pigmentation plasticity. Overall, our study reveals the existence and underlying basis of extensive and trait-specific genetic variation for pigmentation and pigmentation plasticity, offering a rich reservoir of raw material for natural selection to shape the evolution of these traits independently.

Geographic and seasonal variation of the for gene reveal signatures of local adaptation in Drosophila melanogaster

In the early 1980s, the observation that Drosophila melanogaster larvae differed in their foraging behaviour laid the foundation for the work that would later lead to the discovery of the foraging gene (for) and its associated foraging phenotypes, rover and sitter. Since then, the molecular characterization of the for gene and our understanding of the mechanisms that maintain its phenotypic variants in the laboratory have progressed enormously. However, the significance and dynamics of such variation are yet to be investigated in nature. With the advent of next-generation sequencing, it is now possible to identify loci underlying the adaptation of populations in response to environmental variation. Here, I present the results of a genotype-environment association analysis that quantifies variation at the for gene among samples of D. melanogaster structured across space and time. These samples consist of published genomes of adult flies collected worldwide, and at least twice per site of collection (during spring and fall). Both an analysis of genetic differentiation based on Fs⁢t values and an analysis of population structure revealed an east-west gradient in allele frequency. This gradient may be the result of spatially varying selection driven by the seasonality of precipitation. These results support the hypothesis that different patterns of gene flow as expected under models of isolation by distance and potentially isolation by environment are driving genetic differentiation among populations. Overall, this study is essential for understanding the mechanisms underlying the evolution of foraging behaviour in D. melanogaster.

A cosmopolitan inversion facilitates seasonal adaptation in overwintering Drosophila

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

Quantification and modeling of turnover dynamics of de novo transcripts in Drosophila melanogaster

Most of the transcribed eukaryotic genomes are composed of non-coding transcripts. Among these transcripts, some are newly transcribed when compared to outgroups and are referred to as de novo transcripts. De novo transcripts have been shown to play a major role in genomic innovations. However, little is known about the rates at which de novo transcripts are gained and lost in individuals of the same species. Here, we address this gap and estimate the de novo transcript turnover rate with an evolutionary model. We use DNA long reads and RNA short reads from seven geographically remote samples of inbred individuals of Drosophila melanogaster to detect de novo transcripts that are gained on a short evolutionary time scale. Overall, each sampled individual contains around 2500 unspliced de novo transcripts, with most of them being sample specific. We estimate that around 0.15 transcripts are gained per year, and that each gained transcript is lost at a rate around 5× 10-5 per year. This high turnover of transcripts suggests frequent exploration of new genomic sequences within species. These rate estimates are essential to comprehend the process and timescale of de novo gene birth.

Pervasive G × E interactions shape adaptive trajectories and the exploration of the phenotypic space in artificial selection experiments

Quantitative genetics models have shown that long-term selection responses depend on initial variance and mutational influx. Understanding limits of selection requires quantifying the role of mutational variance. However, correlative responses to selection on nonfocal traits can perturb the selection response on the focal trait; and generations are often confounded with selection environments so that genotype by environment (G×E) interactions are ignored. The Saclay divergent selection experiments (DSEs) on maize flowering time were used to track the fate of individual mutations combining genotyping data and phenotyping data from yearly measurements (DSEYM) and common garden experiments (DSECG) with four objectives: (1) to quantify the relative contribution of standing and mutational variance to the selection response, (2) to estimate genotypic mutation effects, (3) to study the impact of G×E interactions in the selection response, and (4) to analyze how trait correlations modulate the exploration of the phenotypic space. We validated experimentally the expected enrichment of fixed beneficial mutations with an average effect of +0.278 and +0.299 days to flowering, depending on the genetic background. Fixation of unfavorable mutations reached up to 25% of incoming mutations, a genetic load possibly due to antagonistic pleiotropy, whereby mutations fixed in the selection environment (DSEYM) turned to be unfavorable in the evaluation environment (DSECG). Global patterns of trait correlations were conserved across genetic backgrounds but exhibited temporal patterns. Traits weakly or uncorrelated with flowering time triggered stochastic exploration of the phenotypic space, owing to microenvironment-specific fixation of standing variants and pleiotropic mutational input.

Evolutionary Insights from a Large-Scale Survey of Population-Genomic Variation

The field of genomics has ushered in new methods for studying molecular-genetic variation in natural populations. However, most population-genomic studies still rely on small sample sizes (typically, <100 individuals) from single time points, leaving considerable uncertainties with respect to the behavior of relatively young (and rare) alleles and, owing to the large sampling variance of measures of variation, to the specific gene targets of unusually strong selection. Genomic sequences of ∼1,700 haplotypes distributed over a 10-year period from a natural population of the microcrustacean Daphnia pulex reveal evolutionary-genomic features at a refined scale, including previously hidden information on the behavior of rare alleles predicted by recent theory. Background selection, resulting from the recurrent introduction of deleterious alleles, appears to strongly influence the dynamics of neutral alleles, inducing indirect negative selection on rare variants and positive selection on common variants. Temporally fluctuating selection increases the persistence of nonsynonymous alleles with intermediate frequencies, while reducing standing levels of variation at linked silent sites. Combined with the results from an equally large metapopulation survey of the study species, classes of genes that are under strong positive selection can now be confidently identified in this key model organism. Most notable among rapidly evolving Daphnia genes are those associated with ribosomes, mitochondrial functions, sensory systems, and lifespan determination.

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

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

Geographic Variation in Genomic Signals of Admixture Between Two Closely Related European Sepsid Fly Species

The extent of interspecific gene flow and its consequences for the initiation, maintenance, and breakdown of species barriers in natural systems remain poorly understood. Interspecific gene flow by hybridization may weaken adaptive divergence, but can be overcome by selection against hybrids, which may ultimately promote reinforcement. An informative step towards understanding the role of gene flow during speciation is to describe patterns of past gene flow among extant species. We investigate signals of admixture between allopatric and sympatric populations of the two closely related European dung fly species Sepsis cynipsea and S. neocynipsea (Diptera: Sepsidae). Based on microsatellite genotypes, we first inferred a baseline demographic history using Approximate Bayesian Computation. We then used genomic data from pooled DNA of natural and laboratory populations to test for past interspecific gene flow based on allelic configurations discordant with the inferred population tree (ABBA-BABA test with D-statistic). Comparing the detected signals of gene flow with the contemporary geographic relationship among interspecific pairs of populations (sympatric vs. allopatric), we made two contrasting observations. At one site in the French Cevennes, we detected an excess of past interspecific gene flow, while at two sites in Switzerland we observed lower signals of past microsatellite genotypes gene flow among populations in sympatry compared to allopatric populations. These results suggest that the species boundaries between these two species depend on the past and/or present eco-geographic context in Europe, which indicates that there is no uniform link between contemporary geographic proximity and past interspecific gene flow in natural populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-023-09612-5.

Timesweeper: accurately identifying selective sweeps using population genomic time series

Despite decades of research, identifying selective sweeps, the genomic footprints of positive selection, remains a core problem in population genetics. Of the myriad methods that have been developed to tackle this task, few are designed to leverage the potential of genomic time-series data. This is because in most population genetic studies of natural populations, only a single period of time can be sampled. Recent advancements in sequencing technology, including improvements in extracting and sequencing ancient DNA, have made repeated samplings of a population possible, allowing for more direct analysis of recent evolutionary dynamics. Serial sampling of organisms with shorter generation times has also become more feasible due to improvements in the cost and throughput of sequencing. With these advances in mind, here we present Timesweeper, a fast and accurate convolutional neural network-based tool for identifying selective sweeps in data consisting of multiple genomic samplings of a population over time. Timesweeper analyzes population genomic time-series data by first simulating training data under a demographic model appropriate for the data of interest, training a one-dimensional convolutional neural network on said simulations, and inferring which polymorphisms in this serialized data set were the direct target of a completed or ongoing selective sweep. We show that Timesweeper is accurate under multiple simulated demographic and sampling scenarios, identifies selected variants with high resolution, and estimates selection coefficients more accurately than existing methods. In sum, we show that more accurate inferences about natural selection are possible when genomic time-series data are available; such data will continue to proliferate in coming years due to both the sequencing of ancient samples and repeated samplings of extant populations with faster generation times, as well as experimentally evolved populations where time-series data are often generated. Methodological advances such as Timesweeper thus have the potential to help resolve the controversy over the role of positive selection in the genome. We provide Timesweeper as a Python package for use by the community.

Genetic differentiation following recent domestication events: A study of farmed Nile tilapia (Oreochromis niloticus) populations

Nile tilapia (Oreochromis niloticus) is among the most farmed finfish worldwide, distributed across different environmental conditions. Its wide distribution has mainly been facilitated by several breeding programs and widespread dissemination of genetically improved strains. In the first Nile tilapia study exploiting a whole-genome pooled sequencing (Poolseq) approach, we identified the genetic structure and signatures of selection in diverse, farmed Nile tilapia populations, with a particular focus on the GIFT strain, developed in the 1980s, and currently managed by WorldFish (GIFTw). We also investigated important farmed strains from The Philippines and Africa. Using both SNP array data and Poolseq SNPs, we characterized the population structure of these samples. We observed the greatest separation between the Asian and African populations and greater admixture in the Asian populations than in the African ones. We also established that the SNP array data were able to successfully resolve relationships between these diverse Nile tilapia populations. The Poolseq data identified genomic regions with high levels of differentiation (F _ST) between GIFTw and the other populations. Gene ontology terms associated with mesoderm development were significantly enriched in the genes located in these regions. A region on chromosome Oni06 was genetically differentiated in pairwise comparisons between GIFTw and all other populations. This region contains genes associated with muscle-related traits and overlaps with a previously published QTL for fillet yield, suggesting that these traits may have been direct targets for selection on GIFT. A nearby region was also identified using XP-EHH to detect genomic differentiation using the SNP array data. Genomic regions with high or extended homozygosity within each population were also identified. This study provides putative genomic landmarks associated with the recent domestication process in several Nile tilapia populations, which could help to inform their genetic management and improvement.

An Ancestral Balanced Inversion Polymorphism Confers Global Adaptation

Since the pioneering work of Dobzhansky in the 1930s and 1940s, many chromosomal inversions have been identified, but how they contribute to adaptation remains poorly understood. In Drosophila melanogaster, the widespread inversion polymorphism In(3R)Payne underpins latitudinal clines in fitness traits on multiple continents. Here, we use single-individual whole-genome sequencing, transcriptomics, and published sequencing data to study the population genomics of this inversion on four continents: in its ancestral African range and in derived populations in Europe, North America, and Australia. Our results confirm that this inversion originated in sub-Saharan Africa and subsequently became cosmopolitan; we observe marked monophyletic divergence of inverted and noninverted karyotypes, with some substructure among inverted chromosomes between continents. Despite divergent evolution of this inversion since its out-of-Africa migration, derived non-African populations exhibit similar patterns of long-range linkage disequilibrium between the inversion breakpoints and major peaks of divergence in its center, consistent with balancing selection and suggesting that the inversion harbors alleles that are maintained by selection on several continents. Using RNA-sequencing, we identify overlap between inversion-linked single-nucleotide polymorphisms and loci that are differentially expressed between inverted and noninverted chromosomes. Expression levels are higher for inverted chromosomes at low temperature, suggesting loss of buffering or compensatory plasticity and consistent with higher inversion frequency in warm climates. Our results suggest that this ancestrally tropical balanced polymorphism spread around the world and became latitudinally assorted along similar but independent climatic gradients, always being frequent in subtropical/tropical areas but rare or absent in temperate climates.

Population genomics reveals mechanisms and dynamics of de novo expressed open reading frame emergence in Drosophila melanogaster

Novel genes are essential for evolutionary innovations and differ substantially even between closely related species. Recently, multiple studies across many taxa showed that some novel genes arise de novo, that is, from previously noncoding DNA. To characterize the underlying mutations that allowed de novo gene emergence and their order of occurrence, homologous regions must be detected within noncoding sequences in closely related sister genomes. So far, most studies do not detect noncoding homologs of de novo genes because of incomplete assemblies and annotations, and long evolutionary distances separating genomes. Here, we overcome these issues by searching for de novo expressed open reading frames (neORFs), the not-yet fixed precursors of de novo genes that emerged within a single species. We sequenced and assembled genomes with long-read technology and the corresponding transcriptomes from inbred lines of Drosophila melanogaster, derived from seven geographically diverse populations. We found line-specific neORFs in abundance but few neORFs shared by lines, suggesting a rapid turnover. Gain and loss of transcription is more frequent than the creation of ORFs, for example, by forming new start and stop codons. Consequently, the gain of ORFs becomes rate limiting and is frequently the initial step in neORFs emergence. Furthermore, transposable elements (TEs) are major drivers for intragenomic duplications of neORFs, yet TE insertions are less important for the emergence of neORFs. However, highly mutable genomic regions around TEs provide new features that enable gene birth. In conclusion, neORFs have a high birth-death rate, are rapidly purged, but surviving neORFs spread neutrally through populations and within genomes.

Evolutionary Insights from a Large-scale Survey of Population-genomic Variation

Results from data on > 1000 haplotypes distributed over a nine-year period from a natural population of the microcrustacean Daphnia pulex reveal evolutionary-genomic features at a refined scale, including key population-genetic properties that are obscured in studies with smaller sample sizes. Background selection, resulting from the recurrent introduction of deleterious alleles, appears to strongly influence the dynamics of neutral alleles, inducing indirect negative selection on rare variants and positive selection on common variants. Fluctuating selection increases the persistence of nonsynonymous alleles with intermediate frequencies, while reducing standing levels of variation at linked silent sites. Combined with the results from an equally large metapopulation survey of the study species, regions of gene structure that are under strong purifying selection and classes of genes that are under strong positive selection in this key species can be confidently identified. Most notable among rapidly evolving Daphnia genes are those associated with ribosomes, mitochondrial functions, sensory systems, and lifespan determination.

A Single Nucleotide Variant in the PPARγ-homolog Eip75B Affects Fecundity in Drosophila

Single nucleotide polymorphisms are the most common type of genetic variation, but how these variants contribute to the adaptation of complex phenotypes is largely unknown. Experimental evolution and genome-wide association studies have demonstrated that variation in the PPARγ-homolog Eip75B has associated with longevity and life-history differences in the fruit fly Drosophila melanogaster. Using RNAi knockdown, we first demonstrate that reduced expression of Eip75B in adult flies affects lifespan, egg-laying rate, and egg volume. We then tested the effects of a naturally occurring SNP within a cis-regulatory domain of Eip75B by applying two complementary approaches: a Mendelian randomization approach using lines of the Drosophila Genetic Reference Panel, and allelic replacement using precise CRISPR/Cas9-induced genome editing. Our experiments reveal that this natural polymorphism has a significant pleiotropic effect on fecundity and egg-to-adult viability, but not on longevity or other life-history traits. Our results provide a rare functional validation at the nucleotide level and identify a natural allelic variant affecting fitness and life-history adaptation.

Natural alleles at the Doa locus underpin evolutionary changes in Drosophila lifespan and fecundity

'Evolve and resequence' (E&R) studies in Drosophila melanogaster have identified many candidate loci underlying the evolution of ageing and life history, but experiments that validate the effects of such candidates remain rare. In a recent E&R study we have identified several alleles of the LAMMER kinase Darkener of apricot (Doa) as candidates for evolutionary changes in lifespan and fecundity. Here, we use two complementary approaches to confirm a functional role of Doa in life-history evolution. First, we used transgenic RNAi to study the effects of Doa at the whole-gene level. Ubiquitous silencing of expression in adult flies reduced both lifespan and fecundity, indicating pleiotropic effects. Second, to characterize segregating variation at Doa, we examined four candidate single nucleotide polymorphisms (SNPs; Doa-1, -2, -3, -4) using a genetic association approach. Three candidate SNPs had effects that were qualitatively consistent with expectations based on our E&R study: Doa-2 pleiotropically affected both lifespan and late-life fecundity; Doa-1 affected lifespan (but not fecundity); and Doa-4 affected late-life fecundity (but not lifespan). Finally, the last candidate allele (Doa-3) also affected lifespan, but in the opposite direction from predicted.

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.

Genetic diversity loss in the Anthropocene

Anthropogenic habitat loss and climate change are reducing species' geographic ranges, increasing extinction risk and losses of species' genetic diversity. Although preserving genetic diversity is key to maintaining species' adaptability, we lack predictive tools and global estimates of genetic diversity loss across ecosystems. We introduce a mathematical framework that bridges biodiversity theory and population genetics to understand the loss of naturally occurring DNA mutations with decreasing habitat. By analyzing genomic variation of 10,095 georeferenced individuals from 20 plant and animal species, we show that genome-wide diversity follows a mutations-area relationship power law with geographic area, which can predict genetic diversity loss from local population extinctions. We estimate that more than 10% of genetic diversity may already be lost for many threatened and nonthreatened species, surpassing the United Nations' post-2020 targets for genetic preservation.

OUP accepted manuscript

There are many problems in biology and related disciplines involving stochasticity, where a signal can only be detected when it lies above a threshold level, while signals lying below threshold are simply not detected. A consequence is that the detected signal is conditioned to lie above threshold, and is not representative of the actual signal. In this work, we present some general results for the conditioning that occurs due to the existence of such an observational threshold. We show that this conditioning is relevant, for example, to gene-frequency trajectories, where many loci in the genome are simultaneously measured in a given generation. Such a threshold can lead to severe biases of allele frequency estimates under purifying selection. In the analysis presented, within the context of Markov chains such as the Wright–Fisher model, we address two key questions: (1) “What is a natural measure of the strength of the conditioning associated with an observation threshold?” (2) “What is a principled way to correct for the effects of the conditioning?”. We answer the first question in terms of a proportion. Starting with a large number of trajectories, the relevant quantity is the proportion of these trajectories that are above threshold at a later time and hence are detected. The smaller the value of this proportion, the stronger the effects of conditioning. We provide an approximate analytical answer to the second question, that corrects the bias produced by an observation threshold, and performs to reasonable accuracy in the Wright–Fisher model for biologically plausible parameter values.

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.