Beneficial fitness effects are not exponential for two viruses

The distribution of fitness effects for beneficial mutations is of paramount importance in determining the outcome of adaptation. It is generally assumed that fitness effects of beneficial mutations follow an exponential distribution, for example, in theoretical treatments of quantitative genetics, clonal interference, experimental evolution, and the adaptation of DNA sequences. This assumption has been justified by the statistical theory of extreme values, because the fitnesses conferred by beneficial mutations should represent samples from the extreme right tail of the fitness distribution. Yet in extreme value theory, there are three different limiting forms for right tails of distributions, and the exponential describes only those of distributions in the Gumbel domain of attraction. Using beneficial mutations from two viruses, we show for the first time that the Gumbel domain can be rejected in favor of a distribution with a right-truncated tail, thus providing evidence for an upper bound on fitness effects. Our data also violate the common assumption that small-effect beneficial mutations greatly outnumber those of large effect, as they are consistent with a uniform distribution of beneficial effects.     

A general extreme value theory model for the adaptation of DNA sequences under strong selection and weak mutation

Recent theoretical studies of the adaptation of DNA sequences assume that the distribution of fitness effects among new beneficial mutations is exponential. This has been justified by using extreme value theory and, in particular, by assuming that the distribution of fitnesses belongs to the Gumbel domain of attraction. However, extreme value theory shows that two other domains of attraction are also possible: the Fréchet and Weibull domains. Distributions in the Fréchet domain have right tails that are heavier than exponential, while distributions in the Weibull domain have right tails that are truncated. To explore the consequences of relaxing the Gumbel assumption, we generalize previous adaptation theory to allow all three domains. We find that many of the previously derived Gumbel-based predictions about the first step of adaptation are fairly robust for some moderate forms of right tails in the Weibull and Fréchet domains, but significant departures are possible, especially for predictions concerning multiple steps in adaptation.     

Distributions of beneficial fitness effects in RNA

Beneficial mutations are the driving force of evolution by natural selection. Yet, relatively little is known about the distribution of the fitness effects of beneficial mutations in populations. Recent work of Gillespie and Orr suggested some of the first generalizations for the distributions of beneficial fitness effects and, surprisingly, they depend only weakly on biological details. In particular, the theory suggests that beneficial mutations obey an exponential distribution of fitness effects, with the same exponential parameter across different regions of genotype space, provided only that few possible beneficial mutations are available to that genotype. Here we tested this hypothesis with a quasi-empirical model of RNA evolution in which fitness is based on the secondary structures of molecules and their thermodynamic stabilities. The fitnesses of randomly selected genotypes appeared to follow a Gumbel-type distribution and thus conform to a basic assumption of adaptation theory. However, the observed distributions of beneficial fitness effects conflict with specific predictions of the theory. In particular, the distributions of beneficial fitness effects appeared exponential only when the vast majority of small-effect beneficial mutations were ignored. Additionally, the distribution of beneficial fitness effects varied with the fitness of the parent genotype. We believe that correlation of the fitness values among similar genotypes is likely the cause of the departure from the predictions of recent adaptation theory. Although in conflict with the current theory, these results suggest that more complex statistical generalizations about beneficial mutations may be possible.     

The distribution of fitness effects among beneficial mutations in Fisher's geometric model of adaptation

Recent models of adaptation at the DNA sequence level assume that the fitness effects of new mutations show certain statistical properties. In particular, these models assume that the distribution of fitness effects among new mutations is in the domain of attraction of the so-called Gumbel-type extreme value distribution. This assumption has not, however, been justified on any biological or theoretical grounds. In this note, I study random mutation in one of the simplest models of mutation and adaptation-Fisher's geometric model. I show that random mutation in this model yields a distribution of mutational effects that belongs to the Gumbel type. I also show that the distribution of fitness effects among rare beneficial mutations in Fisher's model is asymptotically exponential. I confirm these analytic findings with exact computer simulations. These results provide some support for the use of Gumbel-type extreme value theory in studies of adaptation and point to a surprising connection between recent phenotypic- and sequence-based models of adaptation: in both, the distribution of fitness effects among rare beneficial mutations is approximately exponential.     

The distribution of beneficial and fixed mutation fitness effects close to an optimum

The distribution of the selection coefficients of beneficial mutations is pivotal to the study of the adaptive process, both at the organismal level (theories of adaptation) and at the gene level (molecular evolution). A now famous result of extreme value theory states that this distribution is an exponential, at least when considering a well-adapted wild type. However, this prediction could be inaccurate under selection for an optimum (because fitness effect distributions have a finite right tail in this case). In this article, we derive the distribution of beneficial mutation effects under a general model of stabilizing selection, with arbitrary selective and mutational covariance between a finite set of traits. We assume a well-adapted wild type, thus taking advantage of the robustness of tail behaviors, as in extreme value theory. We show that, under these general conditions, both beneficial mutation effects and fixed effects (mutations escaping drift loss) are beta distributed. In both cases, the parameters have explicit biological meaning and are empirically measurable; their variation through time can also be predicted. We retrieve the classic exponential distribution as a subcase of the beta when there are a moderate to large number of weakly correlated traits under selection. In this case too, we provide an explicit biological interpretation of the parameters of the distribution. We show by simulations that these conclusions are fairly robust to a lower adaptation of the wild type and discuss the relevance of our findings in the context of adaptation theories and experimental evolution.     

Quantifying the adaptive potential of an antibiotic resistance enzyme

For a quantitative understanding of the process of adaptation, we need to understand its "raw material," that is, the frequency and fitness effects of beneficial mutations. At present, most empirical evidence suggests an exponential distribution of fitness effects of beneficial mutations, as predicted for Gumbel-domain distributions by extreme value theory. Here, we study the distribution of mutation effects on cefotaxime (Ctx) resistance and fitness of 48 unique beneficial mutations in the bacterial enzyme TEM-1 β-lactamase, which were obtained by screening the products of random mutagenesis for increased Ctx resistance. Our contributions are threefold. First, based on the frequency of unique mutations among more than 300 sequenced isolates and correcting for mutation bias, we conservatively estimate that the total number of first-step mutations that increase Ctx resistance in this enzyme is 87 [95% CI 75-189], or 3.4% of all 2,583 possible base-pair substitutions. Of the 48 mutations, 10 are synonymous and the majority of the 38 non-synonymous mutations occur in the pocket surrounding the catalytic site. Second, we estimate the effects of the mutations on Ctx resistance by determining survival at various Ctx concentrations, and we derive their fitness effects by modeling reproduction and survival as a branching process. Third, we find that the distribution of both measures follows a Fréchet-type distribution characterized by a broad tail of a few exceptionally fit mutants. Such distributions have fundamental evolutionary implications, including an increased predictability of evolution, and may provide a partial explanation for recent observations of striking parallel evolution of antibiotic resistance.     

The distribution of fitness effects among beneficial mutations

We know little about the distribution of fitness effects among new beneficial mutations, a problem that partly reflects the rarity of these changes. Surprisingly, though, population genetic theory allows us to predict what this distribution should look like under fairly general assumptions. Using extreme value theory, I derive this distribution and show that it has two unexpected properties. First, the distribution of beneficial fitness effects at a gene is exponential. Second, the distribution of beneficial effects at a gene has the same mean regardless of the fitness of the present wild-type allele. Adaptation from new mutations is thus characterized by a kind of invariance: natural selection chooses from the same spectrum of beneficial effects at a locus independent of the fitness rank of the present wild type. I show that these findings are reasonably robust to deviations from several assumptions. I further show that one can back calculate the mean size of new beneficial mutations from the observed mean size of fixed beneficial mutations.     

Exploiting the Adaptation Dynamics to Predict the Distribution of Beneficial Fitness Effects

Adaptation of asexual populations is driven by beneficial mutations and therefore the dynamics of this process, besides other factors, depends on the distribution of beneficial fitness effects. It is known that on uncorrelated fitness landscapes, this distribution can only be of three types: truncated, exponential and power law. We performed extensive stochastic simulations to study the adaptation dynamics on rugged fitness landscapes, and identified two quantities that can be used to distinguish the underlying distribution of beneficial fitness effects. The first quantity studied here is the fitness difference between successive mutations that spread in the population, which is found to decrease in the case of truncated distributions, remains nearly a constant for exponentially decaying distributions and increases when the fitness distribution decays as a power law. The second quantity of interest, namely, the rate of change of fitness with time also shows quantitatively different behaviour for different beneficial fitness distributions. The patterns displayed by the two aforementioned quantities are found to hold good for both low and high mutation rates. We discuss how these patterns can be exploited to determine the distribution of beneficial fitness effects in microbial experiments.     

The Distribution of Fitness Effects of Beneficial Mutations in Pseudomonas aeruginosa

Understanding how beneficial mutations affect fitness is crucial to our understanding of adaptation by natural selection. Here, using adaptation to the antibiotic rifampicin in the opportunistic pathogen Pseudomonas aeruginosa as a model system, we investigate the underlying distribution of fitness effects of beneficial mutations on which natural selection acts. Consistent with theory, the effects of beneficial mutations are exponentially distributed where the fitness of the wild type is moderate to high. However, when the fitness of the wild type is low, the data no longer follow an exponential distribution, because many beneficial mutations have large effects on fitness. There is no existing population genetic theory to explain this bias towards mutations of large effects, but it can be readily explained by the underlying biochemistry of rifampicin–RNA polymerase interactions. These results demonstrate the limitations of current population genetic theory for predicting adaptation to severe sources of stress, such as antibiotics, and they highlight the utility of integrating statistical and biophysical approaches to adaptation.     

ADAPTIVE WALKS AND DISTRIBUTION OF BENEFICIAL FITNESS EFFECTS

We study the adaptation dynamics of a maladapted asexual population on rugged fitness landscapes with many local fitness peaks. The distribution of beneficial fitness effects is assumed to belong to one of the three extreme value domains, viz. Weibull, Gumbel, and Fréchet. We work in the strong selection‐weak mutation regime in which beneficial mutations fix sequentially, and the population performs an uphill walk on the fitness landscape until a local fitness peak is reached. A striking prediction of our analysis is that the fitness difference between successive steps follows a pattern of diminishing returns in the Weibull domain and accelerating returns in the Fréchet domain, as the initial fitness of the population is increased. These trends are found to be robust with respect to fitness correlations. We believe that this result can be exploited in experiments to determine the extreme value domain of the distribution of beneficial fitness effects. Our work here differs significantly from the previous ones that assume the selection coefficient to be small. On taking large effect mutations into account, we find that the length of the walk shows different qualitative trends from those derived using small selection coefficient approximation.     

Distributions of dimeric tandem repeats in non-coding and coding DNA sequences

We study the length distribution functions for the 16 possible distinct dimeric tandem repeats in DNA sequences of diverse taxonomic partitions of GenBank (known human and mouse genomes, and complete genomes of Caenorhabditis elegans and yeast). For coding DNA, we find that all 16 distribution functions are exponential. For non-coding DNA, the distribution functions for most of the dimeric repeats have surprisingly long tails, that fit a power-law function. We hypothesize that: (i) the exponential distributions of dimeric repeats in protein coding sequences indicate strong evolutionary pressure against tandem repeat expansion in coding DNA sequences; and (ii) long tails in the distributions of dimers in non-coding DNA may be a result of various mutational mechanisms. These long, non-exponential tails in the distribution of dimeric repeats in non-coding DNA are hypothesized to be due to the higher tolerance of non-coding DNA to mutations. By comparing genomes of various phylogenetic types of organisms, we find that the shapes of the distributions are not universal, but rather depend on the specific class of species and the type of a dimer.     

The Properties of Adaptive Walks in Evolving Populations of Fungus

The rarity of beneficial mutations has frustrated efforts to develop a quantitative theory of adaptation. Recent models of adaptive walks, the sequential substitution of beneficial mutations by selection, make two compelling predictions: adaptive walks should be short, and fitness increases should become exponentially smaller as successive mutations fix. We estimated the number and fitness effects of beneficial mutations in each of 118 replicate lineages of Aspergillus nidulans evolving for approximately 800 generations at two population sizes using a novel maximum likelihood framework, the results of which were confirmed experimentally using sexual crosses. We find that adaptive walks do indeed tend to be short, and fitness increases become smaller as successive mutations fix. Moreover, we show that these patterns are associated with a decreasing supply of beneficial mutations as the population adapts. We also provide empirical distributions of fitness effects among mutations fixed at each step. Our results provide a first glimpse into the properties of multiple steps in an adaptive walk in asexual populations and lend empirical support to models of adaptation involving selection towards a single optimum phenotype. In practical terms, our results suggest that the bulk of adaptation is likely to be accomplished within the first few steps.     

Genetic variation and the fate of beneficial mutations in asexual populations

The fate of a newly arising beneficial mutation depends on many factors, such as the population size and the availability and fitness effects of other mutations that accumulate in the population. It has proved difficult to understand how these factors influence the trajectories of particular mutations, since experiments have primarily focused on characterizing successful clones emerging from a small number of evolving populations. Here, we present the results of a massively parallel experiment designed to measure the full spectrum of possible fates of new beneficial mutations in hundreds of experimental yeast populations, whether these mutations are ultimately successful or not. Using strains in which a particular class of beneficial mutation is detectable by fluorescence, we followed the trajectories of these beneficial mutations across 592 independent populations for 1000 generations. We find that the fitness advantage provided by individual mutations plays a surprisingly small role. Rather, underlying "background" genetic variation is quickly generated in our initially clonal populations and plays a crucial role in determining the fate of each individual beneficial mutation in the evolving population.     

Clonal interference, multiple mutations and adaptation in large asexual populations

Two important problems affect the ability of asexual populations to accumulate beneficial mutations and hence to adapt. First, clonal interference causes some beneficial mutations to be outcompeted by more-fit mutations that occur in the same genetic background. Second, multiple mutations occur in some individuals, so even mutations of large effect can be outcompeted unless they occur in a good genetic background that contains other beneficial mutations. In this article, we use a Monte Carlo simulation to study how these two factors influence the adaptation of asexual populations. We find that the results depend qualitatively on the shape of the distribution of the fitness effects of possible beneficial mutations. When this distribution falls off slower than exponentially, clonal interference alone reasonably describes which mutations dominate the adaptation, although it gives a misleading picture of the evolutionary dynamics. When the distribution falls off faster than exponentially, an analysis based on multiple mutations is more appropriate. Using our simulations, we are able to explore the limits of validity of both of these approaches, and we explore the complex dynamics in the regimes where neither one is fully applicable.     

Mutational Spectrum Drives the Rise of Mutator Bacteria

Understanding how mutator strains emerge in bacterial populations is relevant both to evolutionary theory and to reduce the threat they pose in clinical settings. The rise of mutator alleles is understood as a result of their hitchhiking with linked beneficial mutations, although the factors that govern this process remain unclear. A prominent but underappreciated fact is that each mutator allele increases only a specific spectrum of mutational changes. This spectrum has been speculated to alter the distribution of fitness effects of beneficial mutations, potentially affecting hitchhiking. To study this possibility, we analyzed the fitness distribution of beneficial mutations generated from different mutator and wild-type Escherichia coli strains. Using antibiotic resistance as a model system, we show that mutational spectra can alter these distributions substantially, ultimately determining the competitive ability of each strain across environments. Computer simulation showed that the effect of mutational spectrum on hitchhiking dynamics follows a non-linear function, implying that even slight spectrum-dependent fitness differences are sufficient to alter mutator success frequency by several orders of magnitude. These results indicate an unanticipated central role for the mutational spectrum in the evolution of bacterial mutation rates. At a practical level, this study indicates that knowledge of the molecular details of resistance determinants is crucial for minimizing mutator evolution during antibiotic therapy.     

The probability of parallel evolution

Abstract How often will natural selection drive parallel evolution at the DNA sequence level? More precisely, what is the probability that selection will cause two populations that live in identical environments to substitute the same beneficial mutation? Here I show that, under fairly general conditions, the answer is simple: if a wild‐type sequence can mutate to n different beneficial mutations, replicate populations will on average fix the same mutation with probability P= 2/(n + 1). This probability, which is derived using extreme value theory, is independent of most biological details, including the length of the gene in question and the precise distribution of fitness effects among alleles. I conclude that the probability of parallel evolution under natural selection is nearly twice as large as that under neutrality.     

Inferring the distribution of mutational effects on fitness in Drosophila

The properties of the distribution of deleterious mutational effects on fitness (DDME) are of fundamental importance for evolutionary genetics. Since it is extremely difficult to determine the nature of this distribution, several methods using various assumptions about the DDME have been developed, for the purpose of parameter estimation. We apply a newly developed method to DNA sequence polymorphism data from two Drosophila species and compare estimates of the parameters of the distribution of the heterozygous fitness effects of amino acid mutations for several different distribution functions. The results exclude normal and gamma distributions, since these predict too few effectively lethal mutations and power-law distributions as a result of predicting too many lethals. Only the lognormal distribution appears to fit both the diversity data and the frequency of lethals. This DDME arises naturally in complex systems when independent factors contribute multiplicatively to an increase in fitness-reducing damage. Several important parameters, such as the fraction of effectively neutral non-synonymous mutations and the harmonic mean of non-neutral selection coefficients, are robust to the form of the DDME. Our results suggest that the majority of non-synonymous mutations in Drosophila are under effective purifying selection.     

The Evolutionarily Stable Distribution of Fitness Effects

The distribution of fitness effects (DFE) of new mutations is a key parameter in determining the course of evolution. This fact has motivated extensive efforts to measure the DFE or to predict it from first principles. However, just as the DFE determines the course of evolution, the evolutionary process itself constrains the DFE. Here, we analyze a simple model of genome evolution in a constant environment in which natural selection drives the population toward a dynamic steady state where beneficial and deleterious substitutions balance. The distribution of fitness effects at this steady state is stable under further evolution and provides a natural null expectation for the DFE in a population that has evolved in a constant environment for a long time. We calculate how the shape of the evolutionarily stable DFE depends on the underlying population genetic parameters. We show that, in the absence of epistasis, the ratio of beneficial to deleterious mutations of a given fitness effect obeys a simple relationship independent of population genetic details. Finally, we analyze how the stable DFE changes in the presence of a simple form of diminishing-returns epistasis.     

Note on the computation of critical effective population sizes.

This work extends the work of Whitlock in examining the critical effective population sizes from the fixation of both deleterious and beneficial mutations under drift and selection to prevent mutation breakdown of the population. The validity of approximations for the probability of fixation depends on the nature of the assumed distribution for the fitness effect of both types of mutations. Using no approximation for the probability of fixation and assuming a heavy tailed fitness effect distribution, the current model indicates that the coefficients of variation for the fitness effect distributions of both types of mutations and the fitness effect distribution mean for the beneficial mutations are important predictors of the critical effective population size. The current model further predicts that very small populations can be sustained if the fitness effect variances for both types of mutations and the mean for beneficial mutations are large.