Fixation of new alleles and the extinction of small populations: drift load, beneficial alleles, and sexual selection

With a small effective population size, random genetic drift is more important than selection in determining the fate of new alleles. Small populations therefore accumulate deleterious mutations. Left unchecked, the effect of these fixed alleles is to reduce the reproductive capacity of a species, eventually to the point of extinction. New beneficial mutations, if fixed by selection, can restore some of this lost fitness. This paper derives the overall change in fitness due to fixation of new deleterious and beneficial alleles, as a function of the distribution of effects of new mutations and the effective population size. There is a critical effective size below which a population will on average decline in fitness, but above which beneficial mutations allow the population to persist. With reasonable estimates of the relevant parameters, this critical effective size is likely to be a few hundred. Furthermore, sexual selection can act to reduce the fixation probability of deleterious new mutations and increase the probability of fixing new beneficial mutations. Sexual selection can therefore reduce the risk of extinction of small populations.     

Fixation probability and time in subdivided populations

New alleles arising in a population by mutation ultimately are either fixed or lost. Either is possible, for both beneficial and deleterious alleles, because of stochastic changes in allele frequency due to genetic drift. Spatially structured populations differ from unstructured populations in the probability of fixation and the time that this fixation takes. Previous results have generally made many assumptions: that all demes contribute to the next generation in exact proportion to their current sizes, that new mutations are beneficial, and that new alleles have additive effects. In this article these assumptions are relaxed, allowing for an arbitrary distribution among demes of reproductive success, both beneficial and deleterious effects, and arbitrary dominance. The effects of population structure can be expressed with two summary statistics: the effective population size and a variant of Wright's F(ST). In general, the probability of fixation is strongly affected by population structure, as is the expected time to fixation or loss. Population structure changes the effective size of the species, often strongly downward; smaller effective size increases the probability of fixing deleterious alleles and decreases the probability of fixing beneficial alleles. On the other hand, population structure causes an increase in the homozygosity of alleles, which increases the probability of fixing beneficial alleles but somewhat decreases the probability of fixing deleterious alleles. The probability of fixing new beneficial alleles can be simply described by 2hs(1 - F(ST))N(e)/N(tot), where hs is the change in fitness of heterozygotes relative to the ancestral homozygote, F(ST) is a weighted version of Wright's measure of population subdivision, and N(e) and N(tot) are the effective and census sizes, respectively. These results are verified by simulation for a broad range of population structures, including the island model, the stepping-stone model, and a model with extinction and recolonization.     

Mutational effects on the clonal interference phenomenon

We study the process of fixation of beneficial mutations in an asexual population by means of a theoretical model. Particularly, we wish to investigate how the supply of deleterious and beneficial mutations influences the dynamics of the adaptive process of an evolving population. It is well known that the deleterious mutations drastically affect the fate of beneficial mutations. In addition, an increasing supply of favorable mutations, to compensate the decay of the fitness due to the accumulation of deleterious mutations, produces the clonal interference phenomenon where advantageous mutations in distinct lineages compete to reach fixation. This competition imposes a limit to the speed of adaptation of the population. Intuitively, we would expect that the interplay of the two mechanisms would conspire to ensure fixation of only large-effect beneficial mutations. Our results, however, show that beneficial mutations of small effect have an increased probability of fixation when both beneficial and deleterious mutations rates are increased.     

From bad to good: Fitness reversals and the ascent of deleterious mutations

Deleterious mutations are considered a major impediment to adaptation, and there are straightforward expectations for the rate at which they accumulate as a function of population size and mutation rate. In a simulation model of an evolving population of asexually replicating RNA molecules, initially deleterious mutations accumulated at rates nearly equal to that of initially beneficial mutations, without impeding evolutionary progress. As the mutation rate was increased within a moderate range, deleterious mutation accumulation and mean fitness improvement both increased. The fixation rates were higher than predicted by many population-genetic models. This seemingly paradoxical result was resolved in part by the observation that, during the time to fixation, the selection coefficient (s) of initially deleterious mutations reversed to confer a selective advantage. Significantly, more than half of the fixations of initially deleterious mutations involved fitness reversals. These fitness reversals had a substantial effect on the total fitness of the genome and thus contributed to its success in the population. Despite the relative importance of fitness reversals, however, the probabilities of fixation for both initially beneficial and initially deleterious mutations were exceedingly small (on the order of 10⁻⁵ of all mutations).     

The speed of adaptation in large asexual populations

In large asexual populations, beneficial mutations have to compete with each other for fixation. Here, I derive explicit analytic expressions for the rate of substitution and the mean beneficial effect of fixed mutations, under the assumptions that the population size N is large, that the mean effect of new beneficial mutations is smaller than the mean effect of new deleterious mutations, and that new beneficial mutations are exponentially distributed. As N increases, the rate of substitution approaches a constant, which is equal to the mean effect of new beneficial mutations. The mean effect of fixed mutations continues to grow logarithmically with N. The speed of adaptation, measured as the change of log fitness over time, also grows logarithmically with N for moderately large N, and it grows double-logarithmically for extremely large N. Moreover, I derive a simple formula that determines whether at given N beneficial mutations are expected to compete with each other or go to fixation independently. Finally, I verify all results with numerical simulations.     

Fixation probability for lytic viruses: the attachment-lysis model

The fixation probability of a beneficial mutation is extremely sensitive to assumptions regarding the organism's life history. In this article we compute the fixation probability using a life-history model for lytic viruses, a key model organism in experimental studies of adaptation. The model assumes that attachment times are exponentially distributed, but that the lysis time, the time between attachment and host cell lysis, is constant. We assume that the growth of the wild-type viral population is controlled by periodic sampling (population bottlenecks) and also include the possibility that clearance may occur at a constant rate, for example, through washout in a chemostat. We then compute the fixation probability for mutations that increase the attachment rate, decrease the lysis time, increase the burst size, or reduce the probability of clearance. The fixation probability of these four types of beneficial mutations can be vastly different and depends critically on the time between population bottlenecks. We also explore mutations that affect lysis time, assuming that the burst size is constrained by the lysis time, for experimental protocols that sample either free phage or free phage and artificially lysed infected cells. In all cases we predict that the fixation probability of beneficial alleles is remarkably sensitive to the time between population bottlenecks.     

Fixation probabilities depend on life history: fecundity, generation time and survival in a burst-death model

The burst-death model has been developed to describe the life history of organisms with variable generation times and a burst of a fixed number of offspring. The model also includes an optional constant clearance rate, such as washout from a chemostat, and the possibility of sustained periods of population growth followed by severe bottlenecks, as in serial passaging. In this model, a beneficial mutation can either increase the burst rate or the burst size, or reduce the clearance rate, thus increasing survival. In this article we examine the effects of these three possible mechanisms on both the Malthusian fitness and the fixation probability of the lineage. We find that equivalent relative increases in the burst rate or burst size confer equivalent increases in the Malthusian fitness of a lineage, whereas increasing survival typically has a more moderate effect on Malthusian fitness. In contrast, for beneficial mutations that confer the same increase in fitness, mutations that increase survival are the most likely to fix, followed by mutations that increase the burst rate. Mutations that increase the burst size are the least likely to fix. These results imply that mutant lineages with the highest Malthusian fitness are not, in many cases, the most likely to escape extinction.     

Model and test in a fungus of the probability that beneficial mutations survive drift

Determining the probability of fixation of beneficial mutations is critically important for building predictive models of adaptive evolution. Despite considerable theoretical work, models of fixation probability have stood untested for nearly a century. However, recent advances in experimental and theoretical techniques permit the development of models with testable predictions. We developed a new model for the probability of surviving genetic drift, a major component of fixation probability, for novel beneficial mutations in the fungus Aspergillus nidulans, based on the life-history characteristics of its colony growth on a solid surface. We tested the model by measuring the probability of surviving drift in 11 adapted strains introduced into wild-type populations of different densities. We found that the probability of surviving drift increased with mutant invasion fitness, and decreased with wild-type density, as expected. The model accurately predicted the survival probability for the majority of mutants, yielding one of the first direct tests of the extinction probability of beneficial mutations.     

Hitchhiking of Deleterious Alleles and the Cost of Adaptation in Partially Selfing Species

Self-fertilization is generally seen to be disadvantageous in the long term. It increases genetic drift, which subsequently reduces polymorphism and the efficiency of selection, which also challenges adaptation. However, high selfing rates can increase the fixation probability of recessive beneficial mutations, but existing theory has generally not accounted for the effect of linked sites. Here, we analyze a model for the fixation probability of deleterious mutants that hitchhike with selective sweeps in diploid, partially selfing populations. Approximate analytical solutions show that, conditional on the sweep not being lost by drift, higher inbreeding rates increase the fixation probability of the deleterious allele, due to the resulting reduction in polymorphism and effective recombination. When extending the analysis to consider a distribution of deleterious alleles, as well as the average fitness increase after a sweep, we find that beneficial alleles generally need to be more recessive than the previously assumed dominance threshold (h < 1/2) for selfing to be beneficial from one-locus theory. Our results highlight that recombination aiding the efficiency of selection on multiple loci amplifies the fitness benefits of outcrossing over selfing, compared to results obtained from one-locus theory. This effect additionally increases the parameter range under which obligate outcrossing is beneficial over partial selfing.     

The impact of male-killing bacteria on host evolutionary processes

Male-killing bacteria are maternally inherited endosymbionts that selectively kill male offspring of their arthropod hosts. Using both analytical techniques and computer simulations, we studied the impact of these bacteria on the population genetics of their hosts. In particular, we derived and corroborated formulas for the fixation probability of mutant alleles, mean times to fixation and fixation or extinction, and heterozygosity for varying male-killer prevalence. Our results demonstrate that infections with male-killing bacteria impede the spread of beneficial alleles, facilitate the spread of deleterious alleles, and reduce genetic variation. The reason for this lies in the strongly reduced fitness of infected females combined with no or very limited gene flow from infected females to uninfected individuals. These two properties of male-killer-infected populations reduce the population size relevant for the initial emergence and spread of mutations. In contrast, use of Wright's equation relating sex ratio to effective population size produces misleading predictions. We discuss the relationship to the similar effect of background selection, the impact of other sex-ratio-distorting endosymbionts, and how our results affect the interpretation of empirical data on genetic variation in male-killer-infected populations.     

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.     

Rapid evolutionary escape by large populations from local fitness peaks is likely in nature

Fitness interactions between loci in the genome, or epistasis, can result in mutations that are individually deleterious but jointly beneficial. Such epistasis gives rise to multiple peaks on the genotypic fitness landscape. The problem of evolutionary escape from such local peaks has been a central problem of evolutionary genetics for at least 75 years. Much attention has focused on models of small populations, in which the sequential fixation of valley genotypes carrying individually deleterious mutations operates most quickly owing to genetic drift. However, valley genotypes can also be subject to mutation while transiently segregating, giving rise to copies of the high fitness escape genotype carrying the jointly beneficial mutations. In the absence of genetic recombination, these mutations may then fix simultaneously. The time for this process declines sharply with increasing population size, and it eventually comes to dominate evolutionary behavior. Here we develop an analytic expression for N(crit), the critical population size that defines the boundary between these regimes, which shows that both are likely to operate in nature. Frequent recombination may disrupt high-fitness escape genotypes produced in populations larger than N(crit) before they reach fixation, defining a third regime whose rate again slows with increasing population size. We develop a novel expression for this critical recombination rate, which shows that in large populations the simultaneous fixation of mutations that are beneficial only jointly is unlikely to be disrupted by genetic recombination if their map distance is on the order of the size of single genes. Thus, counterintuitively, mass selection alone offers a biologically realistic resolution to the problem of evolutionary escape from local fitness peaks in natural populations.     

The effect of population bottlenecks on mutation rate evolution in asexual populations

In the absence of recombination, a mutator allele can spread through a population by hitchhiking with beneficial mutations that appear in its genetic background. Theoretical studies over the past decade have shown that the survival and fixation probability of beneficial mutations can be severely reduced by population size bottlenecks. Here, we use computational modelling and evolution experiments with the yeast S. cerevisiae to examine whether population bottlenecks can affect mutator dynamics in adapting asexual populations. In simulation, we show that population bottlenecks can inhibit mutator hitchhiking with beneficial mutations and are most effective at lower beneficial mutation supply rates. We then subjected experimental populations of yeast propagated at the same effective population size to three different bottleneck regimes and observed that the speed of mutator hitchhiking was significantly slower at smaller bottlenecks, consistent with our theoretical expectations. Our results, thus, suggest that bottlenecks can be an important factor in mutation rate evolution and can in certain circumstances act to stabilize or, at least, delay the progressive elevation of mutation rates in asexual populations. Additionally, our findings provide the first experimental support for the theoretically postulated effect of population bottlenecks on beneficial mutations and demonstrate the usefulness of studying mutator frequency dynamics for understanding the underlying dynamics of fitness‐affecting mutations.     

Transposable elements and fitness of bacteria

A stochastic model was designed to describe the evolution of bacterial cultures during 10,000 generations. It is based on a decreasing law for the generation of beneficial mutations as they become fixed in the genomes. Seven beneficial mutations on average were necessary to improve the relative fitness from 1.0 to 1.43 and the model was consistent with the population biology and the genetic data of 12 experimental lines. In one bacterial line, comparison between the model and the data suggests that pivotal mutations mediated by insertion sequences account for a large part of bacterial adaptation. In a more detailed analysis of one simulation, it was shown that only 0.01% of the mutations generated by a population over 10,000 generations can go to fixation as a consequence of their improved fitness. However in the model, the probability of being better fit than its parent should be set initially at ca. 10% to promote an evolution similar to the observed data.     

Compensating for our load of mutations: freezing the meltdown of small populations

We have investigated the reduction of fitness caused by the fixation of new deleterious mutations in small populations within the framework of Fisher's geometrical model of adaptation. In Fisher's model, a population evolves in an n-dimensional character space with an adaptive optimum at the origin. The model allows us to investigate compensatory mutations, which restore fitness losses incurred by other mutations, in a context-dependent manner. We have conducted a moment analysis of the model, supplemented by the numerical results of computer simulations. The mean reduction of fitness (i.e., expected load) scaled to one is approximately n/(n+2Ne), where Ne is the effective population size. The reciprocal relationship between the load and Ne implies that the fixation of deleterious mutations is unlikely to cause extinction when there is a broad scope for compensatory mutations, except in very small populations. Furthermore, the dependence of load on n implies that pleiotropy plays a large role in determining the extinction risk of small populations. Differences and similarities between our results and those of a previous study on the effects of Ne and n are explored. That the predictions of this model are qualitatively different from studies ignoring compensatory mutations implies that we must be cautious in predicting the evolutionary fate of small populations and that additional data on the nature of mutations is of critical importance.     

Beneficial mutation selection balance and the effect of linkage on positive selection

When beneficial mutations are rare, they accumulate by a series of selective sweeps. But when they are common, many beneficial mutations will occur before any can fix, so there will be many different mutant lineages in the population concurrently. In an asexual population, these different mutant lineages interfere and not all can fix simultaneously. In addition, further beneficial mutations can accumulate in mutant lineages while these are still a minority of the population. In this article, we analyze the dynamics of such multiple mutations and the interplay between multiple mutations and interference between clones. These result in substantial variation in fitness accumulating within a single asexual population. The amount of variation is determined by a balance between selection, which destroys variation, and beneficial mutations, which create more. The behavior depends in a subtle way on the population parameters: the population size, the beneficial mutation rate, and the distribution of the fitness increments of the potential beneficial mutations. The mutation-selection balance leads to a continually evolving population with a steady-state fitness variation. This variation increases logarithmically with both population size and mutation rate and sets the rate at which the population accumulates beneficial mutations, which thus also grows only logarithmically with population size and mutation rate. These results imply that mutator phenotypes are less effective in larger asexual populations. They also have consequences for the advantages (or disadvantages) of sex via the Fisher-Muller effect; these are discussed briefly.     

The Probability of Fixation in Populations of Changing Size

The rate of adaptive evolution of a population ultimately depends on the rate of incorporation of beneficial mutations. Even beneficial mutations may, however, be lost from a population since mutant individuals may, by chance, fail to reproduce. In this paper, we calculate the probability of fixation of beneficial mutations that occur in populations of changing size. We examine a number of demographic models, including a population whose size changes once, a population experiencing exponential growth or decline, one that is experiencing logistic growth or decline, and a population that fluctuates in size. The results are based on a branching process model but are shown to be approximate solutions to the diffusion equation describing changes in the probability of fixation over time. Using the diffusion equation, the probability of fixation of deleterious alleles can also be determined for populations that are changing in size. The results developed in this paper can be used to estimate the fixation flux, defined as the rate at which beneficial alleles fix within a population. The fixation flux measures the rate of adaptive evolution of a population and, as we shall see, depends strongly on changes that occur in population size.     

Fitness effects of fixed beneficial mutations in microbial populations

Beneficial mutations are intuitively relevant to understanding adaptation, yet not all beneficial mutations are of consequence to the long-term evolutionary outcome of adaptation. Many beneficial mutations-mostly those of small effect-are lost due either to (1) genetic drift or to (2) competition among clones carrying different beneficial mutations, a phenomenon called the "Hill-Robertson effect" for sexual populations and "clonal interference" for asexual populations. Competition among clones becomes more prevalent with increasing genetic linkage and increasing population size, and it is thus generally characteristic of microbial populations. Together, these two phenomena suggest that only those beneficial mutations of large fitness effect should achieve fixation, despite the fact that most beneficial mutations produced are predicted to have very small fitness effects. Here, we confirm this prediction-both empirically and theoretically-by showing that fitness effects of fixed beneficial mutations follow a distribution whose mode is positive.     

Fixation probabilities when generation times are variable: the burst death model

Estimating the fixation probability of a beneficial mutation has a rich history in theoretical population genetics. Typically, to attain mathematical tractability, we assume that generation times are fixed, while the number of offspring per individual is stochastic. However, fixation probabilities are extremely sensitive to these assumptions regarding life history. In this article, we compute the fixation probability for a "burst-death" life-history model. The model assumes that generation times are exponentially distributed, but the number of offspring per individual is constant. We estimate the fixation probability for populations of constant size and for populations that grow exponentially between periodic population bottlenecks. We find that the fixation probability is, in general, substantially lower in the burst-death model than in classical models. We also note striking qualitative differences between the fates of beneficial mutations that increase burst size and mutations that increase the burst rate. In particular, once the burst size is sufficiently large relative to the wild type, the burst-death model predicts that fixation probability depends only on burst rate.     

Fixation of clonal lineages under Muller's ratchet

For clonal lineages of finite size that differ in their deleterious mutational effects, the probability of fixation is investigated by mathematical theory and Monte Carlo simulations. If these fitness effects are sufficiently small in one or both lineages, then the lineage with the less deleterious effects will become fixed with high probability. If, however, in both lineages the deleterious effects are larger than a threshold s(c), then the probability of fixation is independent of the fitness effects and depends only on the initial frequencies of the lineages. This threshold decreases with decreasing genomic mutation rate U and increases with population size N. (For N = 10(5), we have s(c) approximately = 0.1 if U = 1, and s(c) approximately = 0.015 if U = 0.1). Above the threshold, the competition is not driven by the ratio of mean fitnesses of the lineages, but by the relative sizes of the zero-mutation classes, which are independent of the fitness effects of the mutations. After the loss of the zero-mutation class of a lineage, the other lineage will spread to fixation with high probability and within a short time span. If the mutation rates of the lineages differ substantially, the lineage with the lower mutation rate is fixed with very high probability unless the lineage with the larger mutation rate has very slightly deleterious mutational effects. If the mutation rates differ by not more than a few percent, then the lineage with the higher mutation rate and the more deleterious effects can become fixed with appreciable probability for a certain range of parameters. The independence of the fixation probability on the fitness effects in a single population leads to dramatic effects in metapopulations: lineages with more deleterious effects have a much higher fixation probability. The critical value s(c), above which this phenomenon occurs, decreases as the migration rate between the subpopulations decreases.