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 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.     

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.     

Size distribution of the quasi-circular vegetation patches in the Yellow River Delta,China

Recent empirical studies have suggested that the patch-size distribution of vegetation can be fitted by a power law, truncated power law, or lognormal model to provide explanatory mechanisms for vegetation pattern formation in arid and semiarid regions. However, contradictory results have been reported. Therefore, additional empirical studies are necessary to test the patch-size distribution of vegetation over several regions before it can be considered as an indicator for assessing the discontinuous transition of ecosystems and understanding the mechanisms of vegetation pattern formation. Analogous to arid and semiarid regions of the world, vegetation patterns are characterized by a two-phase mosaic composed of dense vegetation patches interspersed with areas of bare soil, referred to as quasi-circular vegetation patches (QVPs), in the Yellow River Delta (YRD), China. However, research on the patch-size distribution of the QVPs reflecting vegetation patterns and ecosystem functioning is lacking. To fill this gap, for the first time, we examined the patch-size distribution of the QVPs using the fused IKONOS high-spatial-resolution image and evaluated the statistical distributions that better fit the patch size data of the QVPs in the YRD. We found that a power law, truncated power law, or lognormal distribution was not supported in the study area, whereas gamma distribution reasonably fits the size data of QVPs, implying that micro-depressions, combined with the water-limited and salinization environments had considerable effects on vegetation pattern formation. Our results provide helpful insights and suggest that further studies are needed to classify different types of QVPs. Additionally, more efficient approaches need to be used to fit the statistical distributions for elucidating the spatial vegetation patterns in the YRD.     

Testing the extreme value domain of attraction for distributions of beneficial fitness effects

In modeling evolutionary genetics, it is often assumed that mutational effects are assigned according to a continuous probability distribution, and multiple distributions have been used with varying degrees of justification. For mutations with beneficial effects, the distribution currently favored is the exponential distribution, in part because it can be justified in terms of extreme value theory, since beneficial mutations should have fitnesses in the extreme right tail of the fitness distribution. While the appeal to extreme value theory seems justified, the exponential distribution is but one of three possible limiting forms for tail distributions, with the other two loosely corresponding to distributions with right-truncated tails and those with heavy tails. We describe a likelihood-ratio framework for analyzing the fitness effects of beneficial mutations, focusing on testing the null hypothesis that the distribution is exponential. We also describe how to account for missing the smallest-effect mutations, which are often difficult to identify experimentally. This technique makes it possible to apply the test to gain-of-function mutations, where the ancestral genotype is unable to grow under the selective conditions. We also describe how to pool data across experiments, since we expect few possible beneficial mutations in any particular experiment.     

Multiple adaptive substitutions during evolution in novel environments

We consider an asexual population under strong selection-weak mutation conditions evolving on rugged fitness landscapes with many local fitness peaks. Unlike the previous studies in which the initial fitness of the population is assumed to be high, here we start the adaptation process with a low fitness corresponding to a population in a stressful novel environment. For generic fitness distributions, using an analytic argument we find that the average number of steps to a local optimum varies logarithmically with the genotype sequence length and increases as the correlations among genotypic fitnesses increase. When the fitnesses are exponentially or uniformly distributed, using an evolution equation for the distribution of population fitness, we analytically calculate the fitness distribution of fixed beneficial mutations and the walk length distribution.     

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.     

Evolutionary dynamics, epistatic interactions, and biological information

We investigate a definition of biological information that connects population genetics with the tools of information theory by focusing on the distribution of genotypes found in a population. Previous research has treated loci as non-interacting by making specific approximations in the calculation of information-theoretic quantities. We expand earlier mathematical forms to include epistasis, or interactions between mutations at all pairs of loci. Application of our improved measure of biological information to evolution on two-locus, two-allele fitness landscapes demonstrates that mutual information between loci reflects epistatic interaction of mutations. Finally, we consider four-locus, two-allele fitness landscapes with modular structure. As modular interactions are inherently epistatic, we demonstrate that our refined approximation provides insight into the underlying structure of these non-trivial fitness landscapes.     

Fitness landscapes emerging from pharmacodynamic functions in the evolution of multidrug resistance

Adaptive evolution often involves beneficial mutations at more than one locus. In this case, the trajectory and rate of adaptation is determined by the underlying fitness landscape, that is, the fitness values and mutational connectivity of all genotypes under consideration. Drug resistance, especially resistance to multiple drugs simultaneously, is also often conferred by mutations at several loci so that the concept of fitness landscapes becomes important. However, fitness landscapes underlying drug resistance are not static but dependent on drug concentrations, which means they are influenced by the pharmacodynamics of the drugs administered. Here, I present a mathematical framework for fitness landscapes of multidrug resistance based on Hill functions describing how drug concentrations affect fitness. I demonstrate that these ‘pharmacodynamic fitness landscapes’ are characterized by pervasive epistasis that arises through (i) fitness costs of resistance (even when these costs are additive), (ii) nonspecificity of resistance mutations to drugs, in particular cross‐resistance, and (iii) drug interactions (both synergistic and antagonistic). In the latter case, reciprocal drug suppression may even lead to reciprocal sign epistasis, so that the doubly resistant genotype occupies a local fitness peak that may be difficult to access by evolution. Simulations exploring the evolutionary dynamics on some pharmacodynamic fitness landscapes with both constant and changing drug concentrations confirm the crucial role of epistasis in determining the rate of multidrug resistance 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.     

Environment Determines Epistatic Patterns for a ssDNA Virus

Despite the accumulation of substantial quantities of information about epistatic interactions among both deleterious and beneficial mutations in a wide array of experimental systems, neither consistent patterns nor causal explanations for these interactions have yet emerged. Furthermore, the effects of mutations depend on the environment in which they are characterized, implying that the environment may also influence epistatic interactions. Recent work with beneficial mutations for the single-stranded DNA bacteriophage ID11 demonstrated that interactions between pairs of mutations could be understood by means of a simple model that assumes that mutations have additive phenotypic effects and that epistasis arises through a nonlinear phenotype–fitness map with a single intermediate optimum. To determine whether such a model could also explain changes in epistatic patterns associated with changes in environment, we measured epistatic interactions for these same mutations under conditions for which we expected to find the wild-type ID11 at different distances from its phenotypic optimum by assaying fitnesses at three different temperatures: 33°, 37°, and 41°. Epistasis was present and negative under all conditions, but became more pronounced as temperature increased. We found that the additive-phenotypes model explained these patterns as changes in the parameters of the phenotype–fitness map, but that a model that additionally allows the phenotypes to vary across temperatures performed significantly better. Our results show that ostensibly complex patterns of fitness effects and epistasis across environments can be explained by assuming a simple structure for the genotype–phenotype relationship.     

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.     

Coevolution to the edge of chaos: coupled fitness landscapes, poised states, and coevolutionary avalanches

We introduce a broadened framework to study aspects of coevolution based on the NK class of statistical models of rugged fitness landscapes. In these models the fitness contribution of each of N genes in a genotype depends epistatically on K other genes. Increasing epistatic interactions increases the rugged multipeaked character of the fitness landscape. Coevolution is thought of, at the lowest level, as a coupling of landscapes such that adaptive moves by one player deform the landscapes of its immediate partners. In these models we are able to tune the ruggedness of landscapes, how richly intercoupled any two landscapes are, and how many other players interact with each player. All these properties profoundly alter the character of the coevolutionary dynamics. In particular, these parameters govern how readily coevolving ecosystems achieve Nash equilibria, how stable to perturbations such equilibria are, and the sustained mean fitness of coevolving partners. In turn, this raises the possibility that an evolutionary metadynamics due to natural selection may sculpt landscapes and their couplings to achieve coevolutionary systems able to coadapt well. The results suggest that sustained fitness is optimized when landscape ruggedness relative to couplings between landscapes is tuned such that Nash equilibria just tenuously form across the ecosystem. In this poised state, coevolutionary avalanches appear to propagate on all length scales in a power law distribution. Such avalanches may be related to the distribution of small and large extinction events in the record.     

Population Size Affects Adaptation in Complex Ways: Simulations on Empirical Adaptive Landscapes

Do large populations always outcompete smaller ones? Does increasing the mutation rate have a similar effect to increasing the population size, with respect to the adaptation of a population? How important are substitutions in determining the adaptation rate? In this study, we ask how population size and mutation rate interact to affect adaptation on empirical adaptive landscapes. Using such landscapes, we do not need to make many ad hoc assumption about landscape topography, such as about epistatic interactions among mutations or about the distribution of fitness effects. Moreover, we have a better understanding of all the mutations that occur in a population and their effects on the average fitness of the population than we can know in experimental studies. Our results show that the evolutionary dynamics of a population cannot be fully explained by the population mutation rate \(N\mu\); even at constant \(N\mu\), there can be dramatic differences in the adaptation of populations of different sizes. Moreover, the substitution rate of mutations is not always equivalent to the adaptation rate, because we observed populations adapting to high adaptive peaks without fixing any mutations. Finally, in contrast to some theoretical predictions, even on the most rugged landscapes we study, small population size is never an advantage over larger population size. These result show that complex interactions among multiple factors can affect the evolutionary dynamics of populations, and simple models should be taken with caution.     

Exploring the complexity of the HIV-1 fitness landscape

Although fitness landscapes are central to evolutionary theory, so far no biologically realistic examples for large-scale fitness landscapes have been described. Most currently available biological examples are restricted to very few loci or alleles and therefore do not capture the high dimensionality characteristic of real fitness landscapes. Here we analyze large-scale fitness landscapes that are based on predictive models for in vitro replicative fitness of HIV-1. We find that these landscapes are characterized by large correlation lengths, considerable neutrality, and high ruggedness and that these properties depend only weakly on whether fitness is measured in the absence or presence of different antiretrovirals. Accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. While the relative extent to which mutations affect fitness on their own (main effects) or in combination with other mutations (epistasis) is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged, less neutral, and more correlated than expected from the distribution of main effects and epistatic interactions alone. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place.     

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.     

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.     

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.     

Effects of genotype on rates of substitution during experimental evolution

Rates of molecular evolution may vary widely between populations, yet the causes of this variation are still incompletely understood. Genetic differences between populations may make an important contribution to variation in rates of evolution, owing to differences in fitness, population size, mutation rates, or in the distribution of fitness effects (DFEs) of available beneficial mutations. By whole genome sequencing of Escherichia coli populations experimentally evolved in the presence of a quinolone antibiotic, we found that rates of substitution varied by genotype, with evidence for a contribution from a genotype's starting fitness. Subsequent targeted sequencing showed that genotypes with high average substitution rates were more likely to undergo the simultaneous fixation of several mutations, consistent with theoretical models of multiple mutation dynamics. Moreover, patterns of substitution were indicative of epistatic relationships between known resistance mutations.     

Statistical analyses support power law distributions found in neuronal avalanches

The size distribution of neuronal avalanches in cortical networks has been reported to follow a power law distribution with exponent close to -1.5, which is a reflection of long-range spatial correlations in spontaneous neuronal activity. However, identifying power law scaling in empirical data can be difficult and sometimes controversial. In the present study, we tested the power law hypothesis for neuronal avalanches by using more stringent statistical analyses. In particular, we performed the following steps: (i) analysis of finite-size scaling to identify scale-free dynamics in neuronal avalanches, (ii) model parameter estimation to determine the specific exponent of the power law, and (iii) comparison of the power law to alternative model distributions. Consistent with critical state dynamics, avalanche size distributions exhibited robust scaling behavior in which the maximum avalanche size was limited only by the spatial extent of sampling ("finite size" effect). This scale-free dynamics suggests the power law as a model for the distribution of avalanche sizes. Using both the Kolmogorov-Smirnov statistic and a maximum likelihood approach, we found the slope to be close to -1.5, which is in line with previous reports. Finally, the power law model for neuronal avalanches was compared to the exponential and to various heavy-tail distributions based on the Kolmogorov-Smirnov distance and by using a log-likelihood ratio test. Both the power law distribution without and with exponential cut-off provided significantly better fits to the cluster size distributions in neuronal avalanches than the exponential, the lognormal and the gamma distribution. In summary, our findings strongly support the power law scaling in neuronal avalanches, providing further evidence for critical state dynamics in superficial layers of cortex.