Fluctuations in fitness distributions and the effects of weak linked selection on sequence evolution

Limits to selection on standing variation in an asexual population

We consider how a population of N haploid individuals responds to directional selection on standing variation, with no new variation from recombination or mutation. Individuals have trait values z1,…,zN, which are drawn from a distribution ψ; the fitness of individual i is proportional to [Formula: see text] . For illustration, we consider the Laplace and Gaussian distributions, which are parametrised only by the variance V0, and show that for large N, there is a scaling limit which depends on a single parameter NV0. When selection is weak relative to drift (NV0≪1), the variance decreases exponentially at rate 1/N, and the expected ultimate gain in log fitness (scaled by V0), is just NV0, which is the same as Robertson's (1960) prediction for a sexual population. In contrast, when selection is strong relative to drift (NV0≫1), the ultimate gain can be found by approximating the establishment of alleles by a branching process in which each allele competes independently with the population mean and the fittest allele to establish is certain to fix. Then, if the probability of survival to time t∼1/V0 of an allele with value z is P(z), with mean P¯, the winning allele is the fittest of NP¯ survivors drawn from a distribution ψP/P¯. The expected ultimate change is ∼2log(1.15NV0) for a Gaussian distribution, and ∼-12log0.36NV0-log-log0.36NV0 for a Laplace distribution. This approach also predicts the variability of the process, and its dynamics; we show that in the strong selection regime, the expected genetic variance decreases as ∼t-3 at large times. We discuss how these results may be related to selection on standing variation that is spread along a linear chromosome.

A quantitative genetic model of background selection in humans

Across the human genome, there are large-scale fluctuations in genetic diversity caused by the indirect effects of selection. This "linked selection signal" reflects the impact of selection according to the physical placement of functional regions and recombination rates along chromosomes. Previous work has shown that purifying selection acting against the steady influx of new deleterious mutations at functional portions of the genome shapes patterns of genomic variation. To date, statistical efforts to estimate purifying selection parameters from linked selection models have relied on classic Background Selection theory, which is only applicable when new mutations are so deleterious that they cannot fix in the population. Here, we develop a statistical method based on a quantitative genetics view of linked selection, that models how polygenic additive fitness variance distributed along the genome increases the rate of stochastic allele frequency change. By jointly predicting the equilibrium fitness variance and substitution rate due to both strong and weakly deleterious mutations, we estimate the distribution of fitness effects (DFE) and mutation rate across three geographically distinct human samples. While our model can accommodate weaker selection, we find evidence of strong selection operating similarly across all human samples. Although our quantitative genetic model of linked selection fits better than previous models, substitution rates of the most constrained sites disagree with observed divergence levels. We find that a model incorporating selective interference better predicts observed divergence in conserved regions, but overall our results suggest uncertainty remains about the processes generating fitness variation in humans.

A field theoretic approach to non-equilibrium population genetics in the strong selection regime

Natural populations are virtually never observed in equilibrium, yet equilibrium approximations comprise the majority of our understanding of population genetics. Using standard tools from statistical physics, a formalism is presented that re-expresses the stochastic equations describing allelic evolution as a partition functional over all possible allelic trajectories ('paths') governed by selection, mutation, and drift. A perturbative field theory is developed for strong additive selection, relevant to disease variation, that facilitates the straightforward computation of closed-form approximations for time-dependent moments of the allele frequency distribution across a wide range of non-equilibrium scenarios; examples are presented for constant population size, exponential growth, bottlenecks, and oscillatory size, all of which align well to simulations and break down just above the drift barrier. Equilibration times are computed and, even for static population size, generically extend beyond the order 1/s timescale associated with exponential frequency decay. Though the mutation load is largely robust to variable population size, perturbative drift-based corrections to the deterministic trajectory are readily computed. Under strong selection, the variance of a new mutation's frequency (related to homozygosity) is dominated by drift-driven dynamics and a transient increase in variance often occurs prior to equilibrating. The excess kurtosis over skew squared is roughly constant (i.e., independent of selection, provided 2Ns ≳ 5) for static population size, and thus potentially sensitive to deviation from equilibrium. These insights highlight the value of such closed-form approximations, naturally generated from Feynman diagrams in a perturbative field theory, which can simply and accurately capture the parameter dependences describing a variety of non-equilibrium population genetic phenomena of interest.

Population genetics of polymorphism and divergence in rapidly evolving populations

In rapidly evolving populations, numerous beneficial and deleterious mutations can arise and segregate within a population at the same time. In this regime, evolutionary dynamics cannot be analyzed using traditional population genetic approaches that assume that sites evolve independently. Instead, the dynamics of many loci must be analyzed simultaneously. Recent work has made progress by first analyzing the fitness variation within a population, and then studying how individual lineages interact with this traveling fitness wave. However, these "traveling wave" models have previously been restricted to extreme cases where selection on individual mutations is either much faster or much slower than the typical coalescent timescale Tc. In this work, we show how the traveling wave framework can be extended to intermediate regimes in which the scaled fitness effects of mutations (Tcs) are neither large nor small compared to one. This enables us to describe the dynamics of populations subject to a wide range of fitness effects, and in particular, in cases where it is not immediately clear which mutations are most important in shaping the dynamics and statistics of genetic diversity. We use this approach to derive new expressions for the fixation probabilities and site frequency spectra of mutations as a function of their scaled fitness effects, along with related results for the coalescent timescale Tc and the rate of adaptation or Muller's ratchet. We find that competition between linked mutations can have a dramatic impact on the proportions of neutral and selected polymorphisms, which is not simply summarized by the scaled selection coefficient Tcs. We conclude by discussing the implications of these results for population genetic inferences.

Linkage disequilibrium between rare mutations

The statistical associations between mutations, collectively known as linkage disequilibrium, encode important information about the evolutionary forces acting within a population. Yet in contrast to single-site analogues like the site frequency spectrum, our theoretical understanding of linkage disequilibrium remains limited. In particular, little is currently known about how mutations with different ages and fitness costs contribute to expected patterns of linkage disequilibrium, even in simple settings where recombination and genetic drift are the major evolutionary forces. Here, I introduce a forward-time framework for predicting linkage disequilibrium between pairs of neutral and deleterious mutations as a function of their present-day frequencies. I show that the dynamics of linkage disequilibrium become much simpler in the limit that mutations are rare, where they admit a simple heuristic picture based on the trajectories of the underlying lineages. I use this approach to derive analytical expressions for a family of frequency-weighted linkage disequilibrium statistics as a function of the recombination rate, the frequency scale, and the additive and epistatic fitness costs of the mutations. I find that the frequency scale can have a dramatic impact on the shapes of the resulting linkage disequilibrium curves, reflecting the broad range of time scales over which these correlations arise. I also show that the differences between neutral and deleterious linkage disequilibrium are not purely driven by differences in their mutation frequencies and can instead display qualitative features that are reminiscent of epistasis. I conclude by discussing the implications of these results for recent linkage disequilibrium measurements in bacteria. This forward-time approach may provide a useful framework for predicting linkage disequilibrium across a range of evolutionary scenarios.

Within-host evolutionary dynamics of seasonal and pandemic human influenza A viruses in young children

The evolution of influenza viruses is fundamentally shaped by within-host processes. However, the within-host evolutionary dynamics of influenza viruses remain incompletely understood, in part because most studies have focused on infections in healthy adults based on single timepoint data. Here, we analyzed the within-host evolution of 82 longitudinally sampled individuals, mostly young children, infected with A/H1N1pdm09 or A/H3N2 viruses between 2007 and 2009. For A/H1N1pdm09 infections during the 2009 pandemic, nonsynonymous minority variants were more prevalent than synonymous ones. For A/H3N2 viruses in young children, early infection was dominated by purifying selection. As these infections progressed, nonsynonymous variants typically increased in frequency even when within-host virus titers decreased. Unlike the short-lived infections of adults where de novo within-host variants are rare, longer infections in young children allow for the maintenance of virus diversity via mutation-selection balance creating potentially important opportunities for within-host virus evolution.

Cryptic selection forces and dynamic heritability in generalized phenotypic evolution

Individuals with different phenotypes can have widely-varying responses to natural selection, yet many classical approaches to evolutionary dynamics emphasize only how a population's average phenotype increases in fitness over time. However, recent experimental results have produced examples of populations that have multiple fitness peaks, or that experience frequency-dependence that affects the direction and strength of selection on certain individuals. Here, we extend classical fitness gradient formulations of natural selection in order to describe the dynamics of a phenotype distribution in terms of its moments-such as the mean, variance, and skewness. The number of governing equations in our model can be adjusted in order to capture different degrees of detail about the population. We compare our simplified model to direct Wright-Fisher simulations of evolution in several canonical fitness landscapes, and we find that our model provides a low-dimensional description of complex dynamics not typically explained by classical theory, such as cryptic selection forces due to selection on trait ranges, time-variation of the heritability, and nonlinear responses to stabilizing or disruptive selection due to asymmetric trait distributions. In addition to providing a framework for extending general understanding of common qualitative concepts in phenotypic evolution - such as fitness gradients, selection pressures, and heritability - our approach has practical importance for studying evolution in contexts in which genetic analysis is infeasible.

Effective models and the search for quantitative principles in microbial evolution

Microbes evolve rapidly. Yet they do so in idiosyncratic ways, which depend on the specific mutations that are beneficial or deleterious in a given situation. At the same time, some population-level patterns of adaptation are strikingly similar across different microbial systems, suggesting that there may also be simple, quantitative principles that unite these diverse scenarios. We review the search for simple principles in microbial evolution, ranging from the biophysical level to emergent evolutionary dynamics. A key theme has been the use of effective models, which coarse-grain over molecular and cellular details to obtain a simpler description in terms of a few effective parameters. Collectively, these theoretical approaches provide a set of quantitative principles that facilitate understanding, prediction, and potentially control of evolutionary phenomena, though formidable challenges remain due to the ecological complexity of natural populations.

The Effect of Strong Purifying Selection on Genetic Diversity

Purifying selection reduces genetic diversity, both at sites under direct selection and at linked neutral sites. This process, known as background selection, is thought to play an important role in shaping genomic diversity in natural populations. Yet despite its importance, the effects of background selection are not fully understood. Previous theoretical analyses of this process have taken a backward-time approach based on the structured coalescent. While they provide some insight, these methods are either limited to very small samples or are computationally prohibitive. Here, we present a new forward-time analysis of the trajectories of both neutral and deleterious mutations at a nonrecombining locus. We find that strong purifying selection leads to remarkably rich dynamics: neutral mutations can exhibit sweep-like behavior, and deleterious mutations can reach substantial frequencies even when they are guaranteed to eventually go extinct. Our analysis of these dynamics allows us to calculate analytical expressions for the full site frequency spectrum. We find that whenever background selection is strong enough to lead to a reduction in genetic diversity, it also results in substantial distortions to the site frequency spectrum, which can mimic the effects of population expansions or positive selection. Because these distortions are most pronounced in the low and high frequency ends of the spectrum, they become particularly important in larger samples, but may have small effects in smaller samples. We also apply our forward-time framework to calculate other quantities, such as the ultimate fates of polymorphisms or the fitnesses of their ancestral backgrounds.

Beneficial mutation-selection dynamics in finite asexual populations: a free boundary approach

Using a free boundary approach based on an analogy with ice melting models, we propose a deterministic PDE framework to describe the dynamics of fitness distributions in the presence of beneficial mutations with non-epistatic effects on fitness. Contrarily to most approaches based on deterministic models, our framework does not rely on an infinite population size assumption, and successfully captures the transient as well as the long time dynamics of fitness distributions. In particular, consistently with stochastic individual-based approaches or stochastic PDE approaches, it leads to a constant asymptotic rate of adaptation at large times, that most deterministic approaches failed to describe. We derive analytic formulas for the asymptotic rate of adaptation and the full asymptotic distribution of fitness. These formulas depend explicitly on the population size, and are shown to be accurate for a wide range of population sizes and mutation rates, compared to individual-based simulations. Although we were not able to derive an analytic description for the transient dynamics, numerical computations lead to accurate predictions and are computationally efficient compared to stochastic simulations. These computations show that the fitness distribution converges towards a travelling wave with constant speed, and whose profile can be computed analytically.

The Nonstationary Dynamics of Fitness Distributions: Asexual Model with Epistasis and Standing Variation

Various models describe asexual evolution by mutation, selection, and drift. Some focus directly on fitness, typically modeling drift but ignoring or simplifying both epistasis and the distribution of mutation effects (traveling wave models). Others follow the dynamics of quantitative traits determining fitness (Fisher's geometric model), imposing a complex but fixed form of mutation effects and epistasis, and often ignoring drift. In all cases, predictions are typically obtained in high or low mutation rate limits and for long-term stationary regimes, thus losing information on transient behaviors and the effect of initial conditions. Here, we connect fitness-based and trait-based models into a single framework, and seek explicit solutions even away from stationarity. The expected fitness distribution is followed over time via its cumulant generating function, using a deterministic approximation that neglects drift. In several cases, explicit trajectories for the full fitness distribution are obtained for arbitrary mutation rates and standing variance. For nonepistatic mutations, especially with beneficial mutations, this approximation fails over the long term but captures the early dynamics, thus complementing stationary stochastic predictions. The approximation also handles several diminishing returns epistasis models (e.g., with an optimal genotype); it can be applied at and away from equilibrium. General results arise at equilibrium, where fitness distributions display a "phase transition" with mutation rate. Beyond this phase transition, in Fisher's geometric model, the full trajectory of fitness and trait distributions takes a simple form; robust to the details of the mutant phenotype distribution. Analytical arguments are explored regarding why and when the deterministic approximation applies.

Collective Fluctuations in the Dynamics of Adaptation and Other Traveling Waves

The dynamics of adaptation are difficult to predict because it is highly stochastic even in large populations. The uncertainty emerges from random genetic drift arising in a vanguard of particularly fit individuals of the population. Several approaches have been developed to analyze the crucial role of genetic drift on the expected dynamics of adaptation, including the mean fitness of the entire population, or the fate of newly arising beneficial deleterious mutations. However, little is known about how genetic drift causes fluctuations to emerge on the population level, where it becomes palpable as variations in the adaptation speed and the fitness distribution. Yet these phenomena control the decay of genetic diversity and variability in evolution experiments and are key to a truly predictive understanding of evolutionary processes. Here, we show that correlations induced by these emergent fluctuations can be computed at any arbitrary order by a suitable choice of a dynamical constraint. The resulting linear equations exhibit fluctuation-induced terms that amplify short-distance correlations and suppress long-distance ones. These terms, which are in general not small, control the decay of genetic diversity and, for wave-tip dominated ("pulled") waves, lead to anticorrelations between the tip of the wave and the lagging bulk of the population. While it is natural to consider the process of adaptation as a branching random walk in fitness space subject to a constraint (due to finite resources), we show that other traveling wave phenomena in ecology and evolution likewise fall into this class of constrained branching random walks. Our methods, therefore, provide a systematic approach toward analyzing fluctuations in a wide range of population biological processes, such as adaptation, genetic meltdown, species invasions, or epidemics.

Deleterious passengers in adapting populations

Most new mutations are deleterious and are eventually eliminated by natural selection. But in an adapting population, the rapid amplification of beneficial mutations can hinder the removal of deleterious variants in nearby regions of the genome, altering the patterns of sequence evolution. Here, we analyze the interactions between beneficial "driver" mutations and linked deleterious "passengers" during the course of adaptation. We derive analytical expressions for the substitution rate of a deleterious mutation as a function of its fitness cost, as well as the reduction in the beneficial substitution rate due to the genetic load of the passengers. We find that the fate of each deleterious mutation varies dramatically with the rate and spectrum of beneficial mutations and the deleterious substitution rate depends nonmonotonically on the population size and the rate of adaptation. By quantifying this dependence, our results allow us to estimate which deleterious mutations will be likely to fix and how many of these mutations must arise before the progress of adaptation is significantly reduced.

Genetic diversity in the interference selection limit

Pervasive natural selection can strongly influence observed patterns of genetic variation, but these effects remain poorly understood when multiple selected variants segregate in nearby regions of the genome. Classical population genetics fails to account for interference between linked mutations, which grows increasingly severe as the density of selected polymorphisms increases. Here, we describe a simple limit that emerges when interference is common, in which the fitness effects of individual mutations play a relatively minor role. Instead, similar to models of quantitative genetics, molecular evolution is determined by the variance in fitness within the population, defined over an effectively asexual segment of the genome (a "linkage block"). We exploit this insensitivity in a new "coarse-grained" coalescent framework, which approximates the effects of many weakly selected mutations with a smaller number of strongly selected mutations that create the same variance in fitness. This approximation generates accurate and efficient predictions for silent site variability when interference is common. However, these results suggest that there is reduced power to resolve individual selection pressures when interference is sufficiently widespread, since a broad range of parameters possess nearly identical patterns of silent site variability.

Population genomics of rapid adaptation by soft selective sweeps

Organisms can often adapt surprisingly quickly to evolutionary challenges, such as the application of pesticides or antibiotics, suggesting an abundant supply of adaptive genetic variation. In these situations, adaptation should commonly produce 'soft' selective sweeps, where multiple adaptive alleles sweep through the population at the same time, either because the alleles were already present as standing genetic variation or arose independently by recurrent de novo mutations. Most well-known examples of rapid molecular adaptation indeed show signatures of such soft selective sweeps. Here, we review the current understanding of the mechanisms that produce soft sweeps and the approaches used for their identification in population genomic data. We argue that soft sweeps might be the dominant mode of adaptation in many species.