Fixation probability in a two-locus model by the ancestral recombination-selection graph

We use the ancestral influence graph (AIG) for a two-locus, two-allele selection model in the limit of a large population size to obtain an analytic approximation for the probability of ultimate fixation of a single mutant allele A. We assume that this new mutant is introduced at a given locus into a finite population in which a previous mutant allele B is already segregating with a wild type at another linked locus. We deduce that the fixation probability increases as the recombination rate increases if allele A is either in positive epistatic interaction with B and allele B is beneficial or in no epistatic interaction with B and then allele A itself is beneficial. This holds at least as long as the recombination fraction and the selection intensity are small enough and the population size is large enough. In particular this confirms the Hill-Robertson effect, which predicts that recombination renders more likely the ultimate fixation of beneficial mutants at different loci in a population in the presence of random genetic drift even in the absence of epistasis. More importantly, we show that this is true from weak negative epistasis to positive epistasis, at least under weak selection. In the case of deleterious mutants, the fixation probability decreases as the recombination rate increases. This supports Muller's ratchet mechanism to explain the accumulation of deleterious mutants in a population lacking recombination.     

Selection in a subdivided population with local extinction and recolonization

In a subdivided population, local extinction and subsequent recolonization affect the fate of alleles. Of particular interest is the interaction of this force with natural selection. The effect of selection can be weakened by this additional source of stochastic change in allele frequency. The behavior of a selected allele in such a population is shown to be equivalent to that of an allele with a different selection coefficient in an unstructured population with a different size. This equivalence allows use of established results for panmictic populations to predict such quantities as fixation probabilities and mean times to fixation. The magnitude of the quantity N(e)s(e), which determines fixation probability, is decreased by extinction and recolonization. Thus deleterious alleles are more likely to fix, and advantageous alleles less likely to do so, in the presence of extinction and recolonization. Computer simulations confirm that the theoretical predictions of both fixation probabilities and mean times to fixation are good approximations.     

A diffusion approximation for selection and drift in a subdivided population

The population-genetic consequences of population structure are of great interest and have been studied extensively. An area of particular interest is the interaction among population structure, natural selection, and genetic drift. At first glance, different results in this area give very different impressions of the effect of population subdivision on effective population size (N(e)), suggesting that no single value of N(e) can completely characterize a structured population. Results presented here show that a population conforming to Wright's island model of subdivision with genic selection can be related to an idealized panmictic population (a Wright-Fisher population). This equivalent panmictic population has a larger size than the actual population; i.e., N(e) is larger than the actual population size, as expected from many results for this type of population structure. The selection coefficient in the equivalent panmictic population, referred to here as the effective selection coefficient (s(e)), is smaller than the actual selection coefficient (s). This explains how the fixation probability of a selected allele can be unaffected by population subdivision despite the fact that subdivision increases N(e), for the product N(e)s(e) is not altered by subdivision.     

The effect of deleterious alleles on adaptation in asexual populations

We calculate the fixation probability of a beneficial allele that arises as the result of a unique mutation in an asexual population that is subject to recurrent deleterious mutation at rate U. Our analysis is an extension of previous works, which make a biologically restrictive assumption that selection against deleterious alleles is stronger than that on the beneficial allele of interest. We show that when selection against deleterious alleles is weak, beneficial alleles that confer a selective advantage that is small relative to U have greatly reduced probabilities of fixation. We discuss the consequences of this effect for the distribution of effects of alleles fixed during adaptation. We show that a selective sweep will increase the fixation probabilities of other beneficial mutations arising during some short interval afterward. We use the calculated fixation probabilities to estimate the expected rate of fitness improvement in an asexual population when beneficial alleles arise continually at some low rate proportional to U. We estimate the rate of mutation that is optimal in the sense that it maximizes this rate of fitness improvement. Again, this analysis relaxes the assumption made previously that selection against deleterious alleles is stronger than on beneficial alleles.     

Linkage disequilibrium under genetic hitchhiking in finite populations

The model of genetic hitchhiking predicts a reduction in sequence diversity at a neutral locus closely linked to a beneficial allele. In addition, it has been shown that the same process results in a specific pattern of correlations (linkage disequilibrium) between neutral polymorphisms along the chromosome at the time of fixation of the beneficial allele. During the hitchhiking event, linkage disequilibrium on either side of the beneficial allele is built up whereas it is destroyed across the selected site. We derive explicit formulas for the expectation of the covariance measure D and standardized linkage disequilibrium sigma 2D between a pair of polymorphic sites. For our analysis we use the approximation of a star-like genealogy at the selected site. The resulting expressions are approximately correct in the limit of large selection coefficients. Using simulations we show that the resulting pattern of linkage disequilibrium is quickly-i.e., in <0.1N generations-destroyed after the fixation of the beneficial allele for moderately distant neutral loci, where N is the diploid population size.     

Fixation probabilities in a spatially heterogeneous environment

 We consider a simple model of a one-locus, two-allele population inhibiting a two-patch system and experiencing spatially heterogeneous viability selection. The populaton size is finite. We use a diffusion approximation and singular perturbation techniques to find the probability of fixation of a mutant allele. We focus on situations in which each allele is advantageous in one patch and deleterious in the other patch. Our theoretical results support the previous conclusions that, under certain conditions, small populations respond faster to selection than do large populations. We emphasize that knowledge of the dependence of migration rates on population size is crucial in evaluating the effects of population size on the rate of evolution.     

Adaptive evolution of asexual populations under Muller's ratchet

We study the population genetics of adaptation in nonequilibrium haploid asexual populations. We find that the accumulation of deleterious mutations, due to the operation of Muller's ratchet, can considerably reduce the rate of fixation of advantageous alleles. Such reduction can be approximated reasonably well by a reduction in the effective population size. In the absence of Muller's ratchet, a beneficial mutation can only become fixed if it creates the best possible genotype; if Muller's ratchet operates, however, mutations initially arising in a nonoptimal genotype can also become fixed in the population, since the loss of the least-loaded class implies that an initially nonoptimal background can become optimal. We show that, while the rate at which adaptive mutations become fixed is reduced, the rate of fixation of deleterious mutations due to the ratchet is not changed by the presence of beneficial mutations as long as the rate of their occurrence is low and the deleterious effects of mutations (s(d)) are higher than the beneficial effects (s(a)). When s(a) > s(d), the advantage of a beneficial mutation can outweigh the deleterious effects of associated mutations. Under these conditions, a beneficial allele can drag to fixation deleterious mutations initially associated with it at a higher rate than in the absence of advantageous alleles. We propose analytical approximations for the rates of accumulation of deleterious and beneficial mutations. Furthermore, when allowing for the possible occurrence of interference between beneficial alleles, we find that the presence of deleterious mutations of either very weak or very strong effect can marginally increase the rate of accumulation of beneficial mutations over that observed in the absence of such deleterious mutations.     

Fixation probability for a beneficial allele and a mutant strategy in a linear game under weak selection in a finite island model

The effect of population structure on the probability of fixation of a newly introduced mutant under weak selection is studied using a coalescent approach. Wright's island model in a framework of a finite number of demes is assumed and two selection regimes are considered: a beneficial allele model and a linear game among offspring. A first-order approximation of the fixation probability for a single mutant with respect to the intensity of selection is deduced. The approximation requires the calculation of expected coalescence times, under neutrality, for lineages starting from two or three sampled individuals. The results are obtained in a general setting without assumptions on the number of demes, the deme size or the migration rate, which allows for simultaneous coalescence or migration events in the genealogy of the sampled individuals. Comparisons are made with limit cases as the deme size or the number of demes goes to infinity or the migration rate goes to zero for which a diffusion approximation approach is possible. Conditions for selection to favor a mutant strategy replacing a resident strategy in the context of a linear game in a finite island population are addressed.     

Soft sweeps: molecular population genetics of adaptation from standing genetic variation

A population can adapt to a rapid environmental change or habitat expansion in two ways. It may adapt either through new beneficial mutations that subsequently sweep through the population or by using alleles from the standing genetic variation. We use diffusion theory to calculate the probabilities for selective adaptations and find a large increase in the fixation probability for weak substitutions, if alleles originate from the standing genetic variation. We then determine the parameter regions where each scenario-standing variation vs. new mutations-is more likely. Adaptations from the standing genetic variation are favored if either the selective advantage is weak or the selection coefficient and the mutation rate are both high. Finally, we analyze the probability of "soft sweeps," where multiple copies of the selected allele contribute to a substitution, and discuss the consequences for the footprint of selection on linked neutral variation. We find that soft sweeps with weaker selective footprints are likely under both scenarios if the mutation rate and/or the selection coefficient is high.     

The Yule Approximation for the Site Frequency Spectrum after a Selective Sweep

In the area of evolutionary theory, a key question is which portions of the genome of a species are targets of natural selection. Genetic hitchhiking is a theoretical concept that has helped to identify various such targets in natural populations. In the presence of recombination, a severe reduction in sequence diversity is expected around a strongly beneficial allele. The site frequency spectrum is an important tool in genome scans for selection and is composed of the numbers , where is the number of single nucleotide polymorphisms (SNPs) present in from individuals. Previous work has shown that both the number of low- and high-frequency variants are elevated relative to neutral evolution when a strongly beneficial allele fixes. Here, we follow a recent investigation of genetic hitchhiking using a marked Yule process to obtain an analytical prediction of the site frequency spectrum in a panmictic population at the time of fixation of a highly beneficial mutation. We combine standard results from the neutral case with the effects of a selective sweep. As simulations show, the resulting formula produces predictions that are more accurate than previous approaches for the whole frequency spectrum. In particular, the formula correctly predicts the elevation of low- and high-frequency variants and is significantly more accurate than previously derived formulas for intermediate frequency variants.     

Fixation dynamics of beneficial alleles in prokaryotic polyploid chromosomes and plasmids

Theoretical population genetics has been mostly developed for sexually reproducing diploid and for monoploid (haploid) organisms, focusing on eukaryotes. The evolution of bacteria and archaea is often studied by models for the allele dynamics in monoploid populations. However, many prokaryotic organisms harbor multicopy replicons—chromosomes and plasmids—and theory for the allele dynamics in populations of polyploid prokaryotes remains lacking. Here, we present a population genetics model for replicons with multiple copies in the cell. Using this model, we characterize the fixation process of a dominant beneficial mutation at 2 levels: the phenotype and the genotype. Our results show that depending on the mode of replication and segregation, the fixation of the mutant phenotype may precede genotypic fixation by many generations; we term this time interval the heterozygosity window. We furthermore derive concise analytical expressions for the occurrence and length of the heterozygosity window, showing that it emerges if the copy number is high and selection strong. Within the heterozygosity window, the population is phenotypically adapted, while both alleles persist in the population. Replicon ploidy thus allows for the maintenance of genetic variation following phenotypic adaptation and consequently for reversibility in adaptation to fluctuating environmental conditions.     

Contrasting dynamics of a mutator allele in asexual populations of differing size

Mutators have been shown to hitchhike in asexual populations when the anticipated beneficial mutation supply rate of the mutator subpopulation, NU(b) (for subpopulation of size N and beneficial mutation rate U(b)) exceeds that of the wild-type subpopulation. Here, we examine the effect of total population size on mutator dynamics in asexual experimental populations of Saccharomyces cerevisiae. Although mutators quickly hitchhike to fixation in smaller populations, mutator fixation requires more and more time as population size increases; this observed delay in mutator hitchhiking is consistent with the expected effect of clonal interference. Interestingly, despite their higher beneficial mutation supply rate, mutators are supplanted by the wild type in very large populations. We postulate that this striking reversal in mutator dynamics is caused by an interaction between clonal interference, the fitness cost of the mutator allele, and infrequent large-effect beneficial mutations in our experimental populations. Our work thus identifies a potential set of circumstances under which mutator hitchhiking can be inhibited in natural asexual populations, despite recent theoretical predictions that such populations should have a net tendency to evolve ever-higher genomic mutation rates.     

AN ADAPTIVE WALK BY HUMAN IMMUNODEFICIENCY VIRUS TYPE 1 THROUGH A FLUCTUATING FITNESS LANDSCAPE

The mutational landscape model of adaptive sequence evolution has been used to explain an unexpected strong positive linear relationship between marginal fitness and mean site‐specific amino acid frequency in the functionally important HIV‐1 gp120 V3 protein region. The model predicts a positive linear relationship between the probability that a particular beneficial allele, among several, is the next to spread to fixation during an adaptive walk, its transition probability, and the allele's selection coefficient. Here, stochastic simulation is used to confirm the intuition that the linear relationship between transition probabilities and selection coefficients, predicted by the model, should, under fluctuating selection, produce a linear relationship between allele frequency, averaged across populations, and fitness. In addition, these relationships hold for the effective population size and mutation rate of HIV‐1 and for the moderately strong selection observed for V3. A survey of the strength of mutation for diverse organisms suggests that these relationships may be widely applicable.     

Selection in a subdivided population with dominance or local frequency dependence

The interplay between population structure and natural selection is an area of great interest. It is known that certain types of population subdivision do not alter fixation probabilities of selected alleles under genic, frequency-independent selection. In the presence of dominance for fitness or frequency-dependent selection these same types of subdivision can have large effects on fixation probabilities. For example, the barrier to fixation of a fitter allele due to underdominance is reduced by subdivision. Analytic results presented here relate a subdivided population that conforms to a finite island model to an approximately equivalent panmictic population. The size of this equivalent population is different from (larger than) the actual size of the subdivided population. Selection parameters are also different in the hypothetical equivalent population. As expected, the degree of dominance is lower in the equivalent population. The results are not limited to dominance but cover any form of polynomial frequency dependence.     

The fixation probability of a beneficial allele in a population dividing by binary fission

We derive formulae for the fixation probability, P, of a rare benefical allele segregating in a population of fixed size which reproduces by binary fission, in terms of the selection coefficient for the beneficial allele, s. We find that an earlier result P 4s does not depend on the assumption of binary fission, but depends on an assumption about the ordering of events in the life cycle. We find that P 2s for mutations occurring during chromosome replication and P 2.8s for mutations occurring at random times between replication events.     

The Hill-Robertson effect and the evolution of recombination

In finite populations, genetic drift generates interference between selected loci, causing advantageous alleles to be found more often on different chromosomes than on the same chromosome, which reduces the rate of adaptation. This "Hill-Robertson effect" generates indirect selection to increase recombination rates. We present a new method to quantify the strength of this selection. Our model represents a new beneficial allele (A) entering a population as a single copy, while another beneficial allele (B) is sweeping at another locus. A third locus affects the recombination rate between selected loci. Using a branching process model, we calculate the probability distribution of the number of copies of A on the different genetic backgrounds, after it is established but while it is still rare. Then, we use a deterministic model to express the change in frequency of the recombination modifier, due to hitchhiking, as A goes to fixation. We show that this method can give good estimates of selection for recombination. Moreover, it shows that recombination is selected through two different effects: it increases the fixation probability of new alleles, and it accelerates selective sweeps. The relative importance of these two effects depends on the relative times of occurrence of the beneficial alleles.     

Selective sweeps for recessive alleles and for other modes of dominance

A selective sweep describes the reduction of linked genetic variation due to strong positive selection. If s is the fitness advantage of a homozygote for the beneficial allele and h its dominance coefficient, it is usually assumed that h = 1/2, i.e. the beneficial allele is co-dominant. We complement existing theory for selective sweeps by assuming that h is any value in [0, 1]. We show that genetic diversity patterns under selective sweeps with strength s and dominance 0 < h < 1 are similar to co-dominant sweeps with selection strength 2hs. Moreover, we focus on the case h = 0 of a completely recessive beneficial allele. We find that the length of the sweep, i.e. the time from occurrence until fixation of the beneficial allele, is of the order of \({\sqrt{N/s}}\) generations, if N is the population size. Simulations as well as our results show that genetic diversity patterns in the recessive case h = 0 greatly differ from all other cases.     

Impact of demography on extinction/fixation events

In this article we consider diffusion processes modeling the dynamics of multiple allelic proportions (with fixed and varying population size). We are interested in the way alleles extinctions and fixations occur. We first prove that for the Wright–Fisher diffusion process with selection, alleles get extinct successively (and not simultaneously), until the fixation of one last allele. Then we introduce a very general model with selection, competition and Mendelian reproduction, derived from the rescaling of a discrete individual-based dynamics. This multi-dimensional diffusion process describes the dynamics of the population size as well as the proportion of each type in the population. We prove first that alleles extinctions occur successively and second that depending on population size dynamics near extinction, fixation can occur either before extinction almost surely, or not. The proofs of these different results rely on stochastic time changes, integrability of one-dimensional diffusion processes paths and multi-dimensional Girsanov’s tranform.

     

Probability of fixation in a heterogeneous environment

We investigate the probability of fixation of a new mutation arising in a metapopulation that ranges over a heterogeneous selective environment. Using simulations, we test the performance of several approximations of this probability, including a new analytical approximation based on separation of the timescales of selection and migration. We extend all approximations to multideme metapopulations with arbitrary population structure. Our simulations show that no single approximation produces accurate predictions of fixation probabilities for all cases of potential interest. At the limits of low and high migration, previously published approximations are found to be highly accurate. The new separation-of-timescales approach provides the best approximations for intermediate rates of migration among habitats, provided selection is not too intense. For nonzero migration and relatively strong selection, all approximations perform poorly. However, the probability of fixation is bounded above and below by the approximations based on low and high migration limits. Surprisingly, in our simulations with symmetric migration, heterogeneous selection in a metapopulation never decreased-and sometimes substantially increased-the probability of fixation of a new allele compared to metapopulations experiencing homogeneous selection with the same mean selection intensity.