Adaptive walks toward a moving optimum

We investigate how the dynamics and outcomes of adaptation by natural selection are affected by environmental stability by simulating adaptive walks in response to an environmental change of fixed magnitude but variable speed. Here we consider monomorphic lineages that adapt by the sequential fixation of beneficial mutations. This is modeled by selecting short RNA sequences for folding stability and secondary structure conservation at increasing temperatures. Using short RNA sequences allows us to describe adaptive outcomes in terms of genotype (sequence) and phenotype (secondary structure) and to follow the dynamics of fitness increase. We find that slower rates of environmental change affect the dynamics of adaptive walks by reducing the fitness effect of fixed beneficial mutations, as well as by increasing the range of time in which the substitutions of largest effect are likely to occur. In addition, adaptation to slower rates of environmental change results in fitter endpoints with fewer possible end phenotypes relative to lineages that adapt to a sudden change. This suggests that care should be taken when experiments using sudden environmental changes are used to make predictions about adaptive responses to gradual change.     

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 Adaptive Potential of the Middle Domain of Yeast Hsp90

The distribution of fitness effects (DFEs) of new mutations across different environments quantifies the potential for adaptation in a given environment and its cost in others. So far, results regarding the cost of adaptation across environments have been mixed, and most studies have sampled random mutations across different genes. Here, we quantify systematically how costs of adaptation vary along a large stretch of protein sequence by studying the distribution of fitness effects of the same ≈2,300 amino-acid changing mutations obtained from deep mutational scanning of 119 amino acids in the middle domain of the heat shock protein Hsp90 in five environments. This region is known to be important for client binding, stabilization of the Hsp90 dimer, stabilization of the N-terminal-Middle and Middle-C-terminal interdomains, and regulation of ATPase–chaperone activity. Interestingly, we find that fitness correlates well across diverse stressful environments, with the exception of one environment, diamide. Consistent with this result, we find little cost of adaptation; on average only one in seven beneficial mutations is deleterious in another environment. We identify a hotspot of beneficial mutations in a region of the protein that is located within an allosteric center. The identified protein regions that are enriched in beneficial, deleterious, and costly mutations coincide with residues that are involved in the stabilization of Hsp90 interdomains and stabilization of client-binding interfaces, or residues that are involved in ATPase–chaperone activity of Hsp90. Thus, our study yields information regarding the role and adaptive potential of a protein sequence that complements and extends known structural information.     

FISHER'S GEOMETRIC MODEL WITH A MOVING OPTIMUM

Fisher's geometric model has been widely used to study the effects of pleiotropy and organismic complexity on phenotypic adaptation. Here, we study a version of Fisher's model in which a population adapts to a gradually moving optimum. Key parameters are the rate of environmental change, the dimensionality of phenotype space, and the patterns of mutational and selectional correlations. We focus on the distribution of adaptive substitutions, that is, the multivariate distribution of the phenotypic effects of fixed beneficial mutations. Our main results are based on an “adaptive‐walk approximation,” which is checked against individual‐based simulations. We find that (1) the distribution of adaptive substitutions is strongly affected by the ecological dynamics and largely depends on a single composite parameter γ, which scales the rate of environmental change by the “adaptive potential” of the population; (2) the distribution of adaptive substitution reflects the shape of the fitness landscape if the environment changes slowly, whereas it mirrors the distribution of new mutations if the environment changes fast; (3) in contrast to classical models of adaptation assuming a constant optimum, with a moving optimum, more complex organisms evolve via larger adaptive steps.     

Environmental heterogeneity enhances clonal interference

Clonal interference (CI) is a phenomenon that may be important in several asexual microbes. It occurs when population sizes are large and mutation rates to new beneficial alleles are of significant magnitude. Here we explore the role of gene flow and spatial heterogeneity in selection strength in the adaptation of asexuals. We consider a subdivided population of individuals that are adapting, through new beneficial mutations, and that migrate between different patches. The fitness effect of each mutation depends on the patch and all mutations considered are assumed to be unconditionally beneficial. We find that spatial variation in selection pressure affects the rate of adaptive evolution and its qualitative effects depend on the level of gene flow. In particular, we find that both low migration and high levels of heterogeneity lead to enhanced CI. In contrast, for high levels of migration the rate of fixation of adaptive mutations is higher when environmental heterogeneity is present. In addition, we observe that the level of fitness variation is higher and simultaneous fixation of multiple mutations tends to occur in the regime of low migration rates and high heterogeneity.     

The length of adaptive walks is insensitive to starting fitness in Aspergillus nidulans

Adaptation involves the successive substitution of beneficial mutations by selection, a process known as an adaptive walk. Gradualist models of adaptation, which assume that all mutations are small relative to the distance to a fitness optimum, predict that adaptive walks should be longer when the founding genotype is less well adapted. More recent work modeling adaptation as a sequence of moves in phenotype or genotype space predicts, by contrast, much shorter adaptive walks irrespective of the fitness of the founding genotype. Here, we provide what is, to the best of our knowledge, the first direct test of these alternative models, measuring the length of adaptive walks in evolving lineages of fungus that differ initially in fitness. Contrary to the gradualist view, we show that the length of adaptive walks in the fungus Aspergillus nidulans is insensitive to starting fitness and involves just two mutations on average. This arises because poorly adapted populations tend to fix mutations of larger average effect than those of better-adapted populations. Our results suggest that the length of adaptive walks may be independent of the fitness of the founding genotype and, moreover, that poorly adapted populations can quickly adapt to novel environments.     

Properties of selected mutations and genotypic landscapes under Fisher's geometric model

The fitness landscape—the mapping between genotypes and fitness—determines properties of the process of adaptation. Several small genotypic fitness landscapes have recently been built by selecting a handful of beneficial mutations and measuring fitness of all combinations of these mutations. Here, we generate several testable predictions for the properties of these small genotypic landscapes under Fisher's geometric model of adaptation. When the ancestral strain is far from the fitness optimum, we analytically compute the fitness effect of selected mutations and their epistatic interactions. Epistasis may be negative or positive on average depending on the distance of the ancestral genotype to the optimum and whether mutations were independently selected, or coselected in an adaptive walk. Simulations show that genotypic landscapes built from Fisher's model are very close to an additive landscape when the ancestral strain is far from the optimum. However, when it is close to the optimum, a large diversity of landscape with substantial roughness and sign epistasis emerged. Strikingly, small genotypic landscapes built from several replicate adaptive walks on the same underlying landscape were highly variable, suggesting that several realizations of small genotypic landscapes are needed to gain information about the underlying architecture of the fitness landscape.     

Diverse phenotypic and genetic responses to short‐term selection in evolving Escherichia coli populations

Beneficial mutations fuel adaptation by altering phenotypes that enhance the fit of organisms to their environment. However, the phenotypic effects of mutations often depend on ecological context, making the distribution of effects across multiple environments essential to understanding the true nature of beneficial mutations. Studies that address both the genetic basis and ecological consequences of adaptive mutations remain rare. Here, we characterize the direct and pleiotropic fitness effects of a collection of 21 first‐step beneficial mutants derived from naïve and adapted genotypes used in a long‐term experimental evolution of Escherichia coli. Whole‐genome sequencing was able to identify the majority of beneficial mutations. In contrast to previous studies, we find diverse fitness effects of mutations selected in a simple environment and few cases of genetic parallelism. The pleiotropic effects of these mutations were predominantly positive but some mutants were highly antagonistic in alternative environments. Further, the fitness effects of mutations derived from the adapted genotypes were dramatically reduced in nearly all environments. These findings suggest that many beneficial variants are accessible from a single point on the fitness landscape, and the fixation of alternative beneficial mutations may have dramatic consequences for niche breadth reduction via metabolic erosion.     

Genomewide Patterns of Substitution in Adaptively Evolving Populations of the RNA Bacteriophage MS2

Experimental evolution of bacteriophage provides a powerful means of studying the genetics of adaptation, as every substitution contributing to adaptation can be identified and characterized. Here, I use experimental evolution of MS2, an RNA bacteriophage, to study its adaptive response to a novel environment. To this end, three lines of MS2 were adapted to rapid growth and lysis at cold temperature for a minimum of 50 phage generations and subjected to whole-genome sequencing. Using this system, I identified adaptive substitutions, monitored changes in frequency of adaptive mutations through the course of the experiment, and measured the effect on phage growth rate of each substitution. All three lines showed a substantial increase in fitness (a two- to threefold increase in growth rate) due to a modest number of substitutions (three to four). The data show some evidence that the substitutions occurring early in the experiment have larger beneficial effects than later ones, in accordance with the expected diminishing returns relationship between the fitness effects of a mutation and its order of substitution. Patterns of molecular evolution seen here—primarily a paucity of hitchhiking mutations—suggest an abundant supply of beneficial mutations in this system. Nevertheless, some beneficial mutations appear to have been lost, possibly due to accumulation of beneficial mutations on other genetic backgrounds, clonal interference, and negatively epistatic interactions with other beneficial mutations.     

Reciprocal sign epistasis between frequently experimentally evolved adaptive mutations causes a rugged fitness landscape

The fitness landscape captures the relationship between genotype and evolutionary fitness and is a pervasive metaphor used to describe the possible evolutionary trajectories of adaptation. However, little is known about the actual shape of fitness landscapes, including whether valleys of low fitness create local fitness optima, acting as barriers to adaptive change. Here we provide evidence of a rugged molecular fitness landscape arising during an evolution experiment in an asexual population of Saccharomyces cerevisiae. We identify the mutations that arose during the evolution using whole-genome sequencing and use competitive fitness assays to describe the mutations individually responsible for adaptation. In addition, we find that a fitness valley between two adaptive mutations in the genes MTH1 and HXT6/HXT7 is caused by reciprocal sign epistasis, where the fitness cost of the double mutant prohibits the two mutations from being selected in the same genetic background. The constraint enforced by reciprocal sign epistasis causes the mutations to remain mutually exclusive during the experiment, even though adaptive mutations in these two genes occur several times in independent lineages during the experiment. Our results show that epistasis plays a key role during adaptation and that inter-genic interactions can act as barriers between adaptive solutions. These results also provide a new interpretation on the classic Dobzhansky-Muller model of reproductive isolation and display some surprising parallels with mutations in genes often associated with tumors.     

Disentangling genetic and epigenetic determinants of ultrafast adaptation

A major rationale for the advocacy of epigenetically mediated adaptive responses is that they facilitate faster adaptation to environmental challenges. This motivated us to develop a theoretical–experimental framework for disclosing the presence of such adaptation‐speeding mechanisms in an experimental evolution setting circumventing the need for pursuing costly mutation–accumulation experiments. To this end, we exposed clonal populations of budding yeast to a whole range of stressors. By growth phenotyping, we found that almost complete adaptation to arsenic emerged after a few mitotic cell divisions without involving any phenotypic plasticity. Causative mutations were identified by deep sequencing of the arsenic‐adapted populations and reconstructed for validation. Mutation effects on growth phenotypes, and the associated mutational target sizes were quantified and embedded in data‐driven individual‐based evolutionary population models. We found that the experimentally observed homogeneity of adaptation speed and heterogeneity of molecular solutions could only be accounted for if the mutation rate had been near estimates of the basal mutation rate. The ultrafast adaptation could be fully explained by extensive positive pleiotropy such that all beneficial mutations dramatically enhanced multiple fitness components in concert. As our approach can be exploited across a range of model organisms exposed to a variety of environmental challenges, it may be used for determining the importance of epigenetic adaptation‐speeding mechanisms in general.     

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.     

Allele surfing promotes microbial adaptation from standing variation

The coupling of ecology and evolution during range expansions enables mutations to establish at expanding range margins and reach high frequencies. This phenomenon, called allele surfing, is thought to have caused revolutions in the gene pool of many species, most evidently in microbial communities. It has remained unclear, however, under which conditions allele surfing promotes or hinders adaptation. Here, using microbial experiments and simulations, we show that, starting with standing adaptive variation, range expansions generate a larger increase in mean fitness than spatially uniform population expansions. The adaptation gain results from ‘soft’ selective sweeps emerging from surfing beneficial mutations. The rate of these surfing events is shown to sensitively depend on the strength of genetic drift, which varies among strains and environmental conditions. More generally, allele surfing promotes the rate of adaptation per biomass produced, which could help developing biofilms and other resource‐limited populations to cope with environmental challenges.     

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.     

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.     

Long-term adaptation of epistatic genetic networks

Gene networks are likely to govern most traits in nature. Mutations at these genes often show functional epistatic interactions that lead to complex genetic architectures and variable fitness effects in different genetic backgrounds. Understanding how epistatic genetic systems evolve in nature remains one of the great challenges in evolutionary biology. Here we combine an analytical framework with individual-based simulations to generate novel predictions about long-term adaptation of epistatic networks. We find that relative to traits governed by independently evolving genes, adaptation with epistatic gene networks is often characterized by longer waiting times to selective sweeps, lower standing genetic variation, and larger fitness effects of adaptive mutations. This may cause epistatic networks to either adapt more slowly or more quickly relative to a nonepistatic system. Interestingly, epistatic networks may adapt faster even when epistatic effects of mutations are on average deleterious. Further, we study the evolution of epistatic properties of adaptive mutations in gene networks. Our results show that adaptive mutations with small fitness effects typically evolve positive synergistic interactions, whereas adaptive mutations with large fitness effects evolve positive synergistic and negative antagonistic interactions at approximately equal frequencies. These results provide testable predictions for adaptation of traits governed by epistatic networks and the evolution of epistasis within networks.     

Pleiotropic effects of beneficial mutations in Escherichia coli

Micromutational models of adaptation have placed considerable weight on antagonistic pleiotropy as a mechanism that prevents mutations of large effect from achieving fixation. However, there are few empirical studies of the distribution of pleiotropic effects, and no studies that have examined this distribution for a large number of adaptive mutations. Here we examine the form and extent of pleiotropy associated with beneficial mutations in Escherichia coli. To do so, we used a collection of independently evolved genotypes, each of which contains a beneficial mutation that confers increased fitness in a glucose-limited environment. To determine the pleiotropic effects of these mutations, we examined the fitnesses of the mutants in five novel resource environments. Our results show that the majority of mutations had significant fitness effects in alternative resources, such that pleiotropy was common. The predominant form of this pleiotropy was positive--that is, most mutations that conferred increased fitness in glucose also conferred increased fitness in novel resources. We did detect some deleterious pleiotropic effects, but they were primarily limited to one of the five resources, and within this resource, to only a subset of mutants. Although pleiotropic effects were generally positive, fitness levels were lower and more variable on resources that differed most in their mechanisms of uptake and catabolism from that of glucose. Positive pleiotropic effects were strongly correlated in magnitude with their direct effects, but no such correlation was found among mutants with deleterious pleiotropic effects. Whereas previous studies of populations evolved on glucose for longer periods of time showed consistent declines on some of the resources used here, our results suggest that deleterious pleiotropic effects were limited to only a subset of the beneficial mutations available.     

Coalescent under the evolution of coadaptation

Adaptation to novel environments arises either from new beneficial mutations or by utilizing pre‐existing genetic variation. When standing variation is used as the source of new adaptation, fitness effects of alleles may be altered through an environmental change. Alternatively, changes in epistatic genetic backgrounds may convert formerly neutral mutations into beneficial alleles in the new genetic background. By extending the coalescent theory to describe the genealogical histories of two interacting loci, I here investigated the hitchhiking effect of epistatic selection on the amount and pattern of sequence diversity at the linked neutral regions. Assuming a specific form of epistasis between two new mutations that are independently neutral, but together form a coadapted haplotype, I demonstrate that the footprints of epistatic selection differ markedly between the interacting loci depending on the order and relative timing of the two mutational events, even though both mutations are equally essential for the formation of an adaptive gene combination. Our results imply that even when neutrality tests could detect just a single instance of adaptive substitution, there may, in fact, be numerous other hidden mutations that are left undetected, but still play indispensable roles in the evolution of a new adaptation. We expect that the integration of the coalescent framework into the general theory of polygenic inheritance would clarify the connection between factors driving phenotypic evolution and their consequences on underlying DNA sequence changes, which should further illuminate the evolutionary foundation of coadapted systems.     

The genetic basis of parallel and divergent phenotypic responses in evolving populations of Escherichia coli

Pleiotropy plays a central role in theories of adaptation, but little is known about the distribution of pleiotropic effects associated with different adaptive mutations. Previously, we described the phenotypic effects of a collection of independently arising beneficial mutations in Escherichia coli. We quantified their fitness effects in the glucose environment in which they evolved and their pleiotropic effects in five novel resource environments. Here we use a candidate gene approach to associate the phenotypic effects of the mutations with the underlying genetic changes. Among our collection of 27 adaptive mutants, we identified a total of 21 mutations (18 of which were unique) encompassing five different loci or gene regions. There was limited resolution to distinguish among loci based on their fitness effects in the glucose environment, demonstrating widespread parallelism in the direct response to selection. However, substantial heterogeneity in mutant effects was revealed when we examined their pleiotropic effects on fitness in the five novel environments. Substitutions in the same locus clustered together phenotypically, indicating concordance between molecular and phenotypic measures of divergence.