Interplay between Pleiotropy and Secondary Selection Determines Rise and Fall of Mutators in Stress Response

Mutators are clones whose mutation rate is about two to three orders of magnitude higher than the rate of wild-type clones and their roles in adaptive evolution of asexual populations have been controversial. Here we address this problem by using an ab initio microscopic model of living cells, which combines population genetics with a physically realistic presentation of protein stability and protein-protein interactions. The genome of model organisms encodes replication controlling genes (RCGs) and genes modeling the mismatch repair (MMR) complexes. The genotype-phenotype relationship posits that the replication rate of an organism is proportional to protein copy numbers of RCGs in their functional form and there is a production cost penalty for protein overexpression. The mutation rate depends linearly on the concentration of homodimers of MMR proteins. By simulating multiple runs of evolution of populations under various environmental stresses—stationary phase, starvation or temperature-jump—we find that adaptation most often occurs through transient fixation of a mutator phenotype, regardless of the nature of stress. By contrast, the fixation mechanism does depend on the nature of stress. In temperature jump stress, mutators take over the population due to loss of stability of MMR complexes. In contrast, in starvation and stationary phase stresses, a small number of mutators are supplied to the population via epigenetic stochastic noise in production of MMR proteins (a pleiotropic effect), and their net supply is higher due to reduced genetic drift in slowly growing populations under stressful environments. Subsequently, mutators in stationary phase or starvation hitchhike to fixation with a beneficial mutation in the RCGs, (second order selection) and finally a mutation stabilizing the MMR complex arrives, returning the population to a non-mutator phenotype. Our results provide microscopic insights into the rise and fall of mutators in adapting finite asexual populations.     

The Genetic Basis of Phenotypic Adaptation I: Fixation of Beneficial Mutations in the Moving Optimum Model

We study the genetic basis of adaptation in a moving optimum model, in which the optimal value for a quantitative trait increases over time at a constant rate. We first analyze a one-locus two-allele model with recurrent mutation, for which we derive accurate analytical approximations for (i) the time at which a previously deleterious allele becomes beneficial, (ii) the waiting time for a successful new mutation, and (iii) the time the mutant allele needs to reach fixation. On the basis of these results, we show that the shortest total time to fixation is for alleles with intermediate phenotypic effect. We derive an approximation for this “optimal” effect, and we show that it depends in a simple way on a composite parameter, which integrates the ecological parameters and the genetic architecture of the trait. In a second step, we use stochastic computer simulations of a multilocus model to study the order in which mutant alleles with different effects go to fixation. In agreement with the one-locus results, alleles with intermediate effect tend to become fixed earlier than those with either small or large effects. However, the effect size of the fastest mutations differs from the one predicted in the one-locus model. We show how these differences can be explained by two specific effects of multilocus genetics. Finally, we discuss our results in the light of three relevant timescales acting in the system—the environmental, mutation, and fixation timescales—which define three parameter regimes leading to qualitative differences in the adaptive substitution pattern.     

A Coalescent Model for a Sweep of a Unique Standing Variant

The use of genetic polymorphism data to understand the dynamics of adaptation and identify the loci that are involved has become a major pursuit of modern evolutionary genetics. In addition to the classical “hard sweep” hitchhiking model, recent research has drawn attention to the fact that the dynamics of adaptation can play out in a variety of different ways and that the specific signatures left behind in population genetic data may depend somewhat strongly on these dynamics. One particular model for which a large number of empirical examples are already known is that in which a single derived mutation arises and drifts to some low frequency before an environmental change causes the allele to become beneficial and sweeps to fixation. Here, we pursue an analytical investigation of this model, bolstered and extended via simulation study. We use coalescent theory to develop an analytical approximation for the effect of a sweep from standing variation on the genealogy at the locus of the selected allele and sites tightly linked to it. We show that the distribution of haplotypes that the selected allele is present on at the time of the environmental change can be approximated by considering recombinant haplotypes as alleles in the infinite-alleles model. We show that this approximation can be leveraged to make accurate predictions regarding patterns of genetic polymorphism following such a sweep. We then use simulations to highlight which sources of haplotypic information are likely to be most useful in distinguishing this model from neutrality, as well as from other sweep models, such as the classic hard sweep and multiple-mutation soft sweeps. We find that in general, adaptation from a unique standing variant will likely be difficult to detect on the basis of genetic polymorphism data from a single population time point alone, and when it can be detected, it will be difficult to distinguish from other varieties of selective sweeps. Samples from multiple populations and/or time points have the potential to ease this difficulty.     

Natural selection and the genetics of adaptation in threespine stickleback

Growing knowledge of the molecular basis of adaptation in wild populations is expanding the study of natural selection. We summarize ongoing efforts to infer three aspects of natural selection—mechanism, form and history—from the genetics of adaptive evolution in threespine stickleback that colonized freshwater after the last ice age. We tested a mechanism of selection for reduced bony armour in freshwater by tracking genotype and allele frequency changes at an underlying major locus (Ectodysplasin) in transplanted stickleback populations. We inferred disruptive selection on genotypes at the same locus in a population polymorphic for bony armour. Finally, we compared the distribution of phenotypic effect sizes of genes underlying changes in body shape with that predicted by models of adaptive peak shifts following colonization of freshwater. Studies of the effects of selection on genes complement efforts to identify the molecular basis of adaptive differences, and improve our understanding of phenotypic evolution.     

Detecting Selection on Temporal and Spatial Scales: A Genomic Time-Series Assessment of Selective Responses to Devil Facial Tumor Disease

Detecting loci under selection is an important task in evolutionary biology. In conservation genetics detecting selection is key to investigating adaptation to the spread of infectious disease. Loci under selection can be detected on a spatial scale, accounting for differences in demographic history among populations, or on a temporal scale, tracing changes in allele frequencies over time. Here we use these two approaches to investigate selective responses to the spread of an infectious cancer—devil facial tumor disease (DFTD)—that since 1996 has ravaged the Tasmanian devil (Sarcophilus harrisii). Using time-series ‘restriction site associated DNA’ (RAD) markers from populations pre- and post DFTD arrival, and DFTD free populations, we infer loci under selection due to DFTD and investigate signatures of selection that are incongruent among methods, populations, and times. The lack of congruence among populations influenced by DFTD with respect to inferred loci under selection, and the direction of that selection, fail to implicate a consistent selective role for DFTD. Instead genetic drift is more likely driving the observed allele frequency changes over time. Our study illustrates the importance of applying methods with different performance optima e.g. accounting for population structure and background selection, and assessing congruence of the results.     

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.     

A General Model of Negative Frequency Dependent Selection Explains Global Patterns of Human ABO Polymorphism

The ABO locus in humans is characterized by elevated heterozygosity and very similar allele frequencies among populations scattered across the globe. Using knowledge of ABO protein function, we generated a simple model of asymmetric negative frequency dependent selection and genetic drift to explain the maintenance of ABO polymorphism and its loss in human populations. In our models, regardless of the strength of selection, models with large effective population sizes result in ABO allele frequencies that closely match those observed in most continental populations. Populations must be moderately small to fall out of equilibrium and lose either the A or B allele (N_e ≤ 50) and much smaller (N_e ≤ 25) for the complete loss of diversity, which nearly always involved the fixation of the O allele. A pattern of low heterozygosity at the ABO locus where loss of polymorphism occurs in our model is consistent with small populations, such as Native American populations. This study provides a general evolutionary model to explain the observed global patterns of polymorphism at the ABO locus and the pattern of allele loss in small populations. Moreover, these results inform the range of population sizes associated with the recent human colonization of the Americas.     

Key questions in the genetics and genomics of eco-evolutionary dynamics

Increasing acceptance that evolution can be ‘rapid'' (or ‘contemporary'') has generated growing interest in the consequences for ecology. The genetics and genomics of these ‘eco-evolutionary dynamics'' will be—to a large extent—the genetics and genomics of organismal phenotypes. In the hope of stimulating research in this area, I review empirical data from natural populations and draw the following conclusions. (1) Considerable additive genetic variance is present for most traits in most populations. (2) Trait correlations do not consistently oppose selection. (3) Adaptive differences between populations often involve dominance and epistasis. (4) Most adaptation is the result of genes of small-to-modest effect, although (5) some genes certainly have larger effects than the others. (6) Adaptation by independent lineages to similar environments is mostly driven by different alleles/genes. (7) Adaptation to new environments is mostly driven by standing genetic variation, although new mutations can be important in some instances. (8) Adaptation is driven by both structural and regulatory genetic variation, with recent studies emphasizing the latter. (9) The ecological effects of organisms, considered as extended phenotypes, are often heritable. Overall, the study of eco-evolutionary dynamics will benefit from perspectives and approaches that emphasize standing genetic variation in many genes of small-to-modest effect acting across multiple traits and that analyze overall adaptation or ‘fitness''. In addition, increasing attention should be paid to dominance, epistasis and regulatory variation.     

Generalized population models and the nature of genetic drift

The Wright–Fisher model of allele dynamics forms the basis for most theoretical and applied research in population genetics. Our understanding of genetic drift, and its role in suppressing the deterministic forces of Darwinian selection has relied on the specific form of sampling inherent to the Wright–Fisher model and its diffusion limit. Here we introduce and analyze a broad class of forward-time population models that share the same mean and variance as the Wright–Fisher model, but may otherwise differ. The proposed class unifies and further generalizes a number of population-genetic processes of recent interest, including the Λ and Cannings processes. Even though these models all have the same variance effective population size, they encode a rich diversity of alternative forms of genetic drift, with significant consequences for allele dynamics. We characterize in detail the behavior of standard population-genetic quantities across this family of generalized models. Some quantities, such as heterozygosity, remain unchanged; but others, such as neutral absorption times and fixation probabilities under selection, deviate by orders of magnitude from the Wright–Fisher model. We show that generalized population models can produce startling phenomena that differ qualitatively from classical behavior — such as assured fixation of a new mutant despite the presence of genetic drift. We derive the forward-time continuum limits of the generalized processes, analogous to Kimura’s diffusion limit of the Wright–Fisher process, and we discuss their relationships to the Kingman and non-Kingman coalescents. Finally, we demonstrate that some non-diffusive, generalized models are more likely, in certain respects, than the Wright–Fisher model itself, given empirical data from Drosophila populations.     

Evidence that Adaptation in Drosophila Is Not Limited by Mutation at Single Sites

Adaptation in eukaryotes is generally assumed to be mutation-limited because of small effective population sizes. This view is difficult to reconcile, however, with the observation that adaptation to anthropogenic changes, such as the introduction of pesticides, can occur very rapidly. Here we investigate adaptation at a key insecticide resistance locus (Ace) in Drosophila melanogaster and show that multiple simple and complex resistance alleles evolved quickly and repeatedly within individual populations. Our results imply that the current effective population size of modern D. melanogaster populations is likely to be substantially larger (≥100-fold) than commonly believed. This discrepancy arises because estimates of the effective population size are generally derived from levels of standing variation and thus reveal long-term population dynamics dominated by sharp—even if infrequent—bottlenecks. The short-term effective population sizes relevant for strong adaptation, on the other hand, might be much closer to census population sizes. Adaptation in Drosophila may therefore not be limited by waiting for mutations at single sites, and complex adaptive alleles can be generated quickly without fixation of intermediate states. Adaptive events should also commonly involve the simultaneous rise in frequency of independently generated adaptive mutations. These so-called soft sweeps have very distinct effects on the linked neutral polymorphisms compared to the standard hard sweeps in mutation-limited scenarios. Methods for the mapping of adaptive mutations or association mapping of evolutionarily relevant mutations may thus need to be reconsidered.     

Detecting Signatures of Selection Through Haplotype Differentiation Among Hierarchically Structured Populations

The detection of molecular signatures of selection is one of the major concerns of modern population genetics. A widely used strategy in this context is to compare samples from several populations and to look for genomic regions with outstanding genetic differentiation between these populations. Genetic differentiation is generally based on allele frequency differences between populations, which are measured by F_ST or related statistics. Here we introduce a new statistic, denoted hapFLK, which focuses instead on the differences of haplotype frequencies between populations. In contrast to most existing statistics, hapFLK accounts for the hierarchical structure of the sampled populations. Using computer simulations, we show that each of these two features—the use of haplotype information and of the hierarchical structure of populations—significantly improves the detection power of selected loci and that combining them in the hapFLK statistic provides even greater power. We also show that hapFLK is robust with respect to bottlenecks and migration and improves over existing approaches in many situations. Finally, we apply hapFLK to a set of six sheep breeds from Northern Europe and identify seven regions under selection, which include already reported regions but also several new ones. We propose a method to help identifying the population(s) under selection in a detected region, which reveals that in many of these regions selection most likely occurred in more than one population. Furthermore, several of the detected regions correspond to incomplete sweeps, where the favorable haplotype is only at intermediate frequency in the population(s) under selection.     

The effects of migration and drift on local adaptation to a heterogeneous environment

Local adaptation experiments are widely used to quantify the levels of adaptation within a heterogeneous environment. However, theoretical studies generally focus on the probability of fixation of alleles or the mean fitness of populations, rather than local adaptation as it is commonly measured experimentally or in field studies. Here, we develop mathematical models and use them to generate analytical predictions for the level of local adaptation as a function of selection, migration and genetic drift. First, we contrast mean fitness and local adaptation measures and show that the latter can be expressed in a simple and general way as a function of the spatial covariance between population mean phenotype and local environmental conditions. Second, we develop several approximations of a population genetics model to show that the system exhibits different behaviours depending on the rate of migration. The main insights are the following: with intermediate migration, both genetic drift and migration decrease local adaptation; with low migration, drift decreases local adaptation but migration speeds up adaptation; with high migration, genetic drift has no effect on local adaptation. Third, we extend this analysis to cases where the trait under selection is continuous using classical quantitative genetics theory. Finally, we discuss these results in the light of recent experimental work on local adaptation.     

Robust Identification of Local Adaptation from Allele Frequencies

Comparing allele frequencies among populations that differ in environment has long been a tool for detecting loci involved in local adaptation. However, such analyses are complicated by an imperfect knowledge of population allele frequencies and neutral correlations of allele frequencies among populations due to shared population history and gene flow. Here we develop a set of methods to robustly test for unusual allele frequency patterns and correlations between environmental variables and allele frequencies while accounting for these complications based on a Bayesian model previously implemented in the software Bayenv. Using this model, we calculate a set of “standardized allele frequencies” that allows investigators to apply tests of their choice to multiple populations while accounting for sampling and covariance due to population history. We illustrate this first by showing that these standardized frequencies can be used to detect nonparametric correlations with environmental variables; these correlations are also less prone to spurious results due to outlier populations. We then demonstrate how these standardized allele frequencies can be used to construct a test to detect SNPs that deviate strongly from neutral population structure. This test is conceptually related to F_ST and is shown to be more powerful, as we account for population history. We also extend the model to next-generation sequencing of population pools—a cost-efficient way to estimate population allele frequencies, but one that introduces an additional level of sampling noise. The utility of these methods is demonstrated in simulations and by reanalyzing human SNP data from the Human Genome Diversity Panel populations and pooled next-generation sequencing data from Atlantic herring. An implementation of our method is available from http://gcbias.org.     

The Role of Geography in Human Adaptation

Various observations argue for a role of adaptation in recent human evolution, including results from genome-wide studies and analyses of selection signals at candidate genes. Here, we use genome-wide SNP data from the HapMap and CEPH-Human Genome Diversity Panel samples to study the geographic distributions of putatively selected alleles at a range of geographic scales. We find that the average allele frequency divergence is highly predictive of the most extreme F_ST values across the whole genome. On a broad scale, the geographic distribution of putatively selected alleles almost invariably conforms to population clusters identified using randomly chosen genetic markers. Given this structure, there are surprisingly few fixed or nearly fixed differences between human populations. Among the nearly fixed differences that do exist, nearly all are due to fixation events that occurred outside of Africa, and most appear in East Asia. These patterns suggest that selection is often weak enough that neutral processes—especially population history, migration, and drift—exert powerful influences over the fate and geographic distribution of selected alleles.     

A coastal cline in sodium accumulation in Arabidopsis thaliana is driven by natural variation of the sodium transporter AtHKT1;1

The genetic model plant Arabidopsis thaliana, like many plant species, experiences a range of edaphic conditions across its natural habitat. Such heterogeneity may drive local adaptation, though the molecular genetic basis remains elusive. Here, we describe a study in which we used genome-wide association mapping, genetic complementation, and gene expression studies to identify cis-regulatory expression level polymorphisms at the AtHKT1;1 locus, encoding a known sodium (Na(+)) transporter, as being a major factor controlling natural variation in leaf Na(+) accumulation capacity across the global A. thaliana population. A weak allele of AtHKT1;1 that drives elevated leaf Na(+) in this population has been previously linked to elevated salinity tolerance. Inspection of the geographical distribution of this allele revealed its significant enrichment in populations associated with the coast and saline soils in Europe. The fixation of this weak AtHKT1;1 allele in these populations is genetic evidence supporting local adaptation to these potentially saline impacted environments.     

The Repertoire and Dynamics of Evolutionary Adaptations to Controlled Nutrient-Limited Environments in Yeast

The experimental evolution of laboratory populations of microbes provides an opportunity to observe the evolutionary dynamics of adaptation in real time. Until very recently, however, such studies have been limited by our inability to systematically find mutations in evolved organisms. We overcome this limitation by using a variety of DNA microarray-based techniques to characterize genetic changes—including point mutations, structural changes, and insertion variation—that resulted from the experimental adaptation of 24 haploid and diploid cultures of Saccharomyces cerevisiae to growth in either glucose, sulfate, or phosphate-limited chemostats for ∼200 generations. We identified frequent genomic amplifications and rearrangements as well as novel retrotransposition events associated with adaptation. Global nucleotide variation detection in ten clonal isolates identified 32 point mutations. On the basis of mutation frequencies, we infer that these mutations and the subsequent dynamics of adaptation are determined by the batch phase of growth prior to initiation of the continuous phase in the chemostat. We relate these genotypic changes to phenotypic outcomes, namely global patterns of gene expression, and to increases in fitness by 5–50%. We found that the spectrum of available mutations in glucose- or phosphate-limited environments combined with the batch phase population dynamics early in our experiments allowed several distinct genotypic and phenotypic evolutionary pathways in response to these nutrient limitations. By contrast, sulfate-limited populations were much more constrained in both genotypic and phenotypic outcomes. Thus, the reproducibility of evolution varies with specific selective pressures, reflecting the constraints inherent in the system-level organization of metabolic processes in the cell. We were able to relate some of the observed adaptive mutations (e.g., transporter gene amplifications) to known features of the relevant metabolic pathways, but many of the mutations pointed to genes not previously associated with the relevant physiology. Thus, in addition to answering basic mechanistic questions about evolutionary mechanisms, our work suggests that experimental evolution can also shed light on the function and regulation of individual metabolic pathways.     

Ethanol resistance in Drosophila melanogaster has increased in parallel cold‐adapted populations and shows a variable genetic architecture within and between populations

Understanding the genetic properties of adaptive trait evolution is a fundamental crux of biological inquiry that links molecular processes to biological diversity. Important uncertainties persist regarding the genetic predictability of adaptive trait change, the role of standing variation, and whether adaptation tends to result in the fixation of favored variants. Here, we use the recurrent evolution of enhanced ethanol resistance in Drosophila melanogaster during this species’ worldwide expansion as a promising system to add to our understanding of the genetics of adaptation. We find that elevated ethanol resistance has evolved at least three times in different cooler regions of the species’ modern range—not only at high latitude but also in two African high‐altitude regions. Applying a bulk segregant mapping framework, we find that the genetic architecture of ethanol resistance evolution differs substantially not only between our three resistant populations, but also between two crosses involving the same European population. We then apply population genetic scans for local adaptation within our quantitative trait locus regions, and we find potential contributions of genes with annotated roles in spindle localization, membrane composition, sterol and alcohol metabolism, and other processes. We also apply simulation‐based analyses that confirm the variable genetic basis of ethanol resistance and hint at a moderately polygenic architecture. However, these simulations indicate that larger‐scale studies will be needed to more clearly quantify the genetic architecture of adaptive evolution and to firmly connect trait evolution to specific causative loci.     

The contribution of statistical physics to evolutionary biology

Evolutionary biology shares many concepts with statistical physics: both deal with populations, whether of molecules or organisms, and both seek to simplify evolution in very many dimensions. Often, methodologies have undergone parallel and independent development, as with stochastic methods in population genetics. Here, we discuss aspects of population genetics that have embraced methods from physics: non-equilibrium statistical mechanics, travelling waves and Monte-Carlo methods, among others, have been used to study polygenic evolution, rates of adaptation and range expansions. These applications indicate that evolutionary biology can further benefit from interactions with other areas of statistical physics; for example, by following the distribution of paths taken by a population through time.     

Designing gene drives to limit spillover to non-target populations

The prospect of utilizing CRISPR-based gene-drive technology for controlling populations has generated much excitement. However, the potential for spillovers of gene-drive alleles from the target population to non-target populations has raised concerns. Here, using mathematical models, we investigate the possibility of limiting spillovers to non-target populations by designing differential-targeting gene drives, in which the expected equilibrium gene-drive allele frequencies are high in the target population but low in the non-target population. We find that achieving differential targeting is possible with certain configurations of gene-drive parameters, but, in most cases, only under relatively low migration rates between populations. Under high migration, differential targeting is possible only in a narrow region of the parameter space. Because fixation of the gene drive in the non-target population could severely disrupt ecosystems, we outline possible ways to avoid this outcome. We apply our model to two potential applications of gene drives—field trials for malaria-vector gene drives and control of invasive species on islands. We discuss theoretical predictions of key requirements for differential targeting and their practical implications.     

Segregational Drift Constrains the Evolutionary Rate of Prokaryotic Plasmids

Plasmids are extrachromosomal genetic elements in prokaryotes that have been recognized as important drivers of microbial ecology and evolution. Plasmids are found in multiple copies inside their host cell where independent emergence of mutations may lead to intracellular genetic heterogeneity. The intracellular plasmid diversity is thus subject to changes upon cell division. However, the effect of plasmid segregation on plasmid evolution remains understudied. Here, we show that genetic drift during cell division—segregational drift—leads to the rapid extinction of novel plasmid alleles. We established a novel experimental approach to control plasmid allele frequency at the levels of a single cell and the whole population. Following the dynamics of plasmid alleles in an evolution experiment, we find that the mode of plasmid inheritance—random or clustered—is an important determinant of plasmid allele dynamics. Phylogenetic reconstruction of our model plasmid in clinical isolates furthermore reveals a slow evolutionary rate of plasmid-encoded genes in comparison to chromosomal genes. Our study provides empirical evidence that genetic drift in plasmid evolution occurs at multiple levels: the host cell and the population of hosts. Segregational drift has implications for the evolutionary rate heterogeneity of extrachromosomal genetic elements.