Adaptation in Tunably Rugged Fitness Landscapes: The Rough Mount Fuji Model

Much of the current theory of adaptation is based on Gillespie’s mutational landscape model (MLM), which assumes that the fitness values of genotypes linked by single mutational steps are independent random variables. On the other hand, a growing body of empirical evidence shows that real fitness landscapes, while possessing a considerable amount of ruggedness, are smoother than predicted by the MLM. In the present article we propose and analyze a simple fitness landscape model with tunable ruggedness based on the rough Mount Fuji (RMF) model originally introduced by Aita et al. in the context of protein evolution. We provide a comprehensive collection of results pertaining to the topographical structure of RMF landscapes, including explicit formulas for the expected number of local fitness maxima, the location of the global peak, and the fitness correlation function. The statistics of single and multiple adaptive steps on the RMF landscape are explored mainly through simulations, and the results are compared to the known behavior in the MLM model. Finally, we show that the RMF model can explain the large number of second-step mutations observed on a highly fit first-step background in a recent evolution experiment with a microvirid bacteriophage.     

Mutational fitness landscape and drug resistance

Robust technology has been developed to systematically quantify fitness landscapes that provide valuable opportunities to improve our understanding of drug resistance and define new avenues to develop drugs with reduced resistance susceptibility. We outline the critical importance of drug resistance studies and the potential for fitness landscape approaches to contribute to this effort. We describe the major technical advancements in mutational scanning, which is the primary approach used to quantify protein fitness landscapes. There are many complex steps to consider in planning and executing mutational scanning projects including developing a selection scheme, generating mutant libraries, tracking the frequency of variants using next-generation sequencing, and processing and interpreting the data. Key experimental parameters impacting each of these steps are discussed to aid in planning fitness landscape studies. There is a strong need for improved understanding of drug resistance, and fitness landscapes provide a promising new approach.     

The Effect of Bacterial Recombination on Adaptation on Fitness Landscapes with Limited Peak Accessibility

There is ample empirical evidence revealing that fitness landscapes are often complex: the fitness effect of a newly arisen mutation can depend strongly on the allelic state at other loci. However, little is known about the effects of recombination on adaptation on such fitness landscapes. Here, we investigate how recombination influences the rate of adaptation on a special type of complex fitness landscapes. On these landscapes, the mutational trajectories from the least to the most fit genotype are interrupted by genotypes with low relative fitness. We study the dynamics of adapting populations on landscapes with different compositions and numbers of low fitness genotypes, with and without recombination. Our results of the deterministic model (assuming an infinite population size) show that recombination generally decelerates adaptation on these landscapes. However, in finite populations, this deceleration is outweighed by the accelerating Fisher-Muller effect under certain conditions. We conclude that recombination has complex effects on adaptation that are highly dependent on the particular fitness landscape, population size and recombination rate.     

The genetics of speciation: Insights from Fisher's geometric model

Research in speciation genetics has uncovered many robust patterns in intrinsic reproductive isolation, and fitness landscape models have been useful in interpreting these patterns. Here, we examine fitness landscapes based on Fisher's geometric model. Such landscapes are analogous to models of optimizing selection acting on quantitative traits, and have been widely used to study adaptation and the distribution of mutational effects. We show that, with a few modifications, Fisher's model can generate all of the major findings of introgression studies (including “speciation genes” with strong deleterious effects, complex epistasis and asymmetry), and the major patterns in overall hybrid fitnesses (including Haldane's Rule, the speciation clock, heterosis, hybrid breakdown, and male–female asymmetry in the F1). We compare our approach to alternative modeling frameworks that assign fitnesses to genotypes by identifying combinations of incompatible alleles. In some cases, the predictions are importantly different. For example, Fisher's model can explain conflicting empirical results about the rate at which incompatibilities accumulate with genetic divergence. In other cases, the predictions are identical. For example, the quality of reproductive isolation is little affected by the manner in which populations diverge.     

A sharp minimum on the mean number of steps taken in adaptive walks

It was recently conjectured by H.A. Orr that from a random initial point on a random fitness landscape of alphabetic sequences with one-mutation adjacency, chosen from a larger class of landscapes, no adaptive algorithm can arrive at a local optimum in fewer than on average e-1 steps. Here, using an example in which the mean number of steps to a local optimum equals (A-1)/A, where A is the number of distinct "letters" in the "alphabet" from which sequences are constructed, it is shown that as originally stated, the conjecture does not hold. It is also demonstrated that (A-1)/A is a sharp minimum on the mean number of steps taken in adaptive walks on fitness landscapes of alphabetic sequences with one-mutation adjacency. As the example that achieves the new lower bound has properties that are not often considered as potential attributes for fitness landscapes-non-identically distributed fitnesses and negative fitness correlations for adjacent points-a weaker set of conditions characteristic of more commonly studied fitness landscapes is proposed under which the lower bound on the mean length of adaptive walks is conjectured to equal e-1.     

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.     

Towards a general theory of adaptive walks on rugged landscapes

Adaptive evolution, to a large extent, is a complex combinatorial optimization process. In this article we take beginning steps towards developing a general theory of adaptive "walks" via fitter variants in such optimization processes. We introduce the basic idea of a space of entities, each a 1-mutant neighbor of many other entities in the space, and the idea of a fitness ascribed to each entity. Adaptive walks proceed from an initial entity, via fitter neighbors, to locally or globally optimal entities that are fitter than their neighbors. We develop a general theory for the number of local optima, lengths of adaptive walks, and the number of alternative local optima accessible from any given initial entity, for the baseline case of an uncorrelated fitness landscape. Most fitness landscapes are correlated, however. Therefore we develop parts of a universal theory of adaptation on correlated landscapes by adaptive processes that have sufficient numbers of mutations per individual to "jump beyond" the correlation lengths in the underlying landscape. In addition, we explore the statistical character of adaptive walks in two independent complex combinatorial optimization problems, that of evolving a specific cell type in model genetic networks, and that of finding good solutions to the traveling salesman problem. Surprisingly, both show similar statistical features, encouraging the hope that a general theory for adaptive walks on correlated and uncorrelated landscapes can be found. In the final section we explore two limits to the efficacy of selection. The first is new, and surprising: for a wide class of systems, as the complexity of the entities under selection increases, the local optima that are attainable fall progressively closer to the mean properties of the underlying space of entities. This may imply that complex biological systems, such as genetic regulatory systems, are "close" to the mean properties of the ensemble of genomic regulatory systems explored by evolution. The second limit shows that with increasing complexity and a fixed mutation rate, selection often becomes unable to pull an adapting population to those local optima to which connected adaptive walks via fitter variants exist. These beginning steps in theory development are applied to maturation of the immune response, and to the problem of radiation and stasis. Despite the limitations of the adaptive landscape metaphor, we believe that further development along the lines begun here will prove useful.     

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.     

Evolutionary accessibility of mutational pathways

Functional effects of different mutations are known to combine to the total effect in highly nontrivial ways. For the trait under evolutionary selection ('fitness'), measured values over all possible combinations of a set of mutations yield a fitness landscape that determines which mutational states can be reached from a given initial genotype. Understanding the accessibility properties of fitness landscapes is conceptually important in answering questions about the predictability and repeatability of evolutionary adaptation. Here we theoretically investigate accessibility of the globally optimal state on a wide variety of model landscapes, including landscapes with tunable ruggedness as well as neutral 'holey' landscapes. We define a mutational pathway to be accessible if it contains the minimal number of mutations required to reach the target genotype, and if fitness increases in each mutational step. Under this definition accessibility is high, in the sense that at least one accessible pathway exists with a substantial probability that approaches unity as the dimensionality of the fitness landscape (set by the number of mutational loci) becomes large. At the same time the number of alternative accessible pathways grows without bounds. We test the model predictions against an empirical 8-locus fitness landscape obtained for the filamentous fungus Aspergillus niger. By analyzing subgraphs of the full landscape containing different subsets of mutations, we are able to probe the mutational distance scale in the empirical data. The predicted effect of high accessibility is supported by the empirical data and is very robust, which we argue reflects the generic topology of sequence spaces. Together with the restrictive assumptions that lie in our definition of accessibility, this implies that the globally optimal configuration should be accessible to genome wide evolution, but the repeatability of evolutionary trajectories is limited owing to the presence of a large number of alternative mutational pathways.     

植物景观遗传学研究进展

植物景观遗传学是新兴的景观遗传学交叉学科的一个重要研究方向。目前植物景观遗传学的研究虽落后于动物,但其在生物多样性保护方面具有的巨大潜力不可忽视。从景观特征对遗传结构、环境因素对适应性遗传变异影响两个方面,系统综述了近十年来国际上植物景观遗传学的研究焦点和研究进展,比较了植物景观遗传学与动物景观遗传学研究在研究设计和研究方法上的异同,并基于将来植物景观遗传学由对空间遗传结构的描述发展为对景观遗传效应的量化分析及预测的发展框架,具体针对目前景观特征与遗传结构研究设计的系统性差、遗传结构与景观格局在时间上的误配、适应性位点与环境变量的模糊匹配、中性遗传变异与适应性遗传变异研究的分隔、景观与遗传关系分析方法的局限等五个方面提出了研究对策。     

Adaptive evolutionary walks require neutral intermediates in RNA fitness landscapes

In RNA fitness landscapes with interconnected networks of neutral mutations, neutral precursor mutations can play an important role in facilitating the accessibility of epistatic adaptive mutant combinations. I use an exhaustively surveyed fitness landscape model based on short sequence RNA genotypes (and their secondary structure phenotypes) to calculate the minimum rate at which mutants initially appearing as neutral are incorporated into an adaptive evolutionary walk. I show first, that incorporating neutral mutations significantly increases the number of point mutations in a given evolutionary walk when compared to estimates from previous adaptive walk models. Second, that incorporating neutral mutants into such a walk significantly increases the final fitness encountered on that walk - indeed evolutionary walks including neutral steps often reach the global optimum in this model. Third, and perhaps most importantly, evolutionary paths of this kind are often extremely winding in their nature and have the potential to undergo multiple mutations at a given sequence position within a single walk; the potential of these winding paths to mislead phylogenetic reconstruction is briefly considered.     

THE INEVITABILITY OF UNCONDITIONALLY DELETERIOUS SUBSTITUTIONS DURING ADAPTATION

Studies on the genetics of adaptation from new mutations typically neglect the possibility that a deleterious mutation might fix. Nonetheless, here we show that, in many regimes, the first mutation to fix is most often deleterious, even when fitness is expected to increase in the long term. In particular, we prove that this phenomenon occurs under weak mutation for any house‐of‐cards model with an equilibrium distribution. We find that the same qualitative results hold under Fisher's geometric model. We also provide a simple intuition for the surprising prevalence of unconditionally deleterious substitutions during early adaptation. Importantly, the phenomenon we describe occurs on fitness landscapes without any local maxima and is therefore distinct from “valley crossing.” Our results imply that the common practice of ignoring deleterious substitutions leads to qualitatively incorrect predictions in many regimes. Our results also have implications for the substitution process at equilibrium and for the response to a sudden decrease in population size.     

Thermoregulation, Flight, and the Evolution of Wing Pattern in Pierid Butterflies: The Topography of Adaptive Landscapes

SYNOPSIS. This paper describes a case study of adaptation, constraint,and evolutionary innovation in pierid butterflies. I developa framework for discussing these issues that focuses on thequestions: What is the form of the adaptive landscape relatingfitness to phenotypic characters? How do such landscapes differfor evolutionarily related groups? I examine the evolution ofwing pigment patterns and thermoregulatory behavior for butterfliesin two subfamilies in the family Pieridae, with three principalresults. First, I show that thermoregulation can be an importantcomponent of fitness in pierids, and that wing color and thermoregulatorybehavior are important phenotypic characters determining thermoregulatoryperformance and the adaptive landscape. Second, I show how limitson possible variation in wing color and behavior constrain evolutionwithin one subfamily of pierids, and how these constraints areset by the physical and biochemical mechanisms of adaptation.Third, I show how evolutionary innovation may have resultedfrom the addition of a new, behavioral dimension to the landscape,and how this addition has altered the functional interrelationsamong various elements of the wing color pattern. I suggestthat comparative analyses of the form and determinants of theadaptive landscape may be useful in identifying evolutionaryinnovations, and complement theoretical analyses of evolutionarydynamics on such fitness surfaces.     

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.     

The Time Scale of Evolutionary Innovation

A fundamental question in biology is the following: what is the time scale that is needed for evolutionary innovations? There are many results that characterize single steps in terms of the fixation time of new mutants arising in populations of certain size and structure. But here we ask a different question, which is concerned with the much longer time scale of evolutionary trajectories: how long does it take for a population exploring a fitness landscape to find target sequences that encode new biological functions? Our key variable is the length, of the genetic sequence that undergoes adaptation. In computer science there is a crucial distinction between problems that require algorithms which take polynomial or exponential time. The latter are considered to be intractable. Here we develop a theoretical approach that allows us to estimate the time of evolution as function of We show that adaptation on many fitness landscapes takes time that is exponential in even if there are broad selection gradients and many targets uniformly distributed in sequence space. These negative results lead us to search for specific mechanisms that allow evolution to work on polynomial time scales. We study a regeneration process and show that it enables evolution to work in polynomial time.     

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.     

Context dependence in complex adaptive landscapes: frequency and trait‐dependent selection surfaces within an adaptive radiation of Caribbean pupfishes

The adaptive landscape provides the foundational bridge between micro‐ and macroevolution. One well‐known caveat to this perspective is that fitness surfaces depend on ecological context, including competitor frequency, traits measured, and resource abundance. However, this view is based largely on intraspecific studies. It is still unknown how context‐dependence affects the larger features of peaks and valleys on the landscape which ultimately drive speciation and adaptive radiation. Here, I explore this question using one of the most complex fitness landscapes measured in the wild in a sympatric pupfish radiation endemic to San Salvador Island, Bahamas by tracking survival and growth of laboratory‐reared F2 hybrids. I present new analyses of the effects of competitor frequency, dietary isotopes, and trait subsets on this fitness landscape. Contrary to expectations, decreasing competitor frequency increased survival only among very common phenotypes, whereas less common phenotypes rarely survived despite few competitors, suggesting that performance, not competitor frequency, shapes large‐scale features of the fitness landscape. Dietary isotopes were weakly correlated with phenotype and growth, but did not explain additional survival variation. Nonlinear fitness surfaces varied substantially among trait subsets, revealing one‐, two‐, and three‐peak landscapes, demonstrating the complexity of selection in the wild, even among similar functional traits.     

The rank ordering of genotypic fitness values predicts genetic constraint on natural selection on landscapes lacking sign epistasis

Sewall Wright's genotypic fitness landscape makes explicit one mechanism by which epistasis for fitness can constrain evolution by natural selection. Wright distinguished between landscapes possessing multiple fitness peaks and those with only a single peak and emphasized that the former class imposes substantially greater constraint on natural selection. Here I present novel formalism that more finely partitions the universe of possible fitness landscapes on the basis of the rank ordering of their genotypic fitness values. In this report I focus on fitness landscapes lacking sign epistasis (i.e., landscapes that lack mutations the sign of whose fitness effect varies epistatically), which constitute a subset of Wright's single peaked landscapes. More than one fitness rank ordering lacking sign epistasis exists for L > 2 (where L is the number of interacting loci), and I find that a highly statistically significant effect exists between landscape membership in fitness rank-ordering partition and two different proxies for genetic constraint, even within this subset of landscapes. This statistical association is robust to population size, permitting general inferences about some of the characteristics of fitness rank orderings responsible for genetic constraint on natural selection.     

Evolution in alternating environments with tunable interlandscape correlations

Natural populations are often exposed to temporally varying environments. Evolutionary dynamics in varying environments have been extensively studied, although understanding the effects of varying selection pressures remains challenging. Here, we investigate how cycling between a pair of statistically related fitness landscapes affects the evolved fitness of an asexually reproducing population. We construct pairs of fitness landscapes that share global fitness features but are correlated with one another in a tunable way, resulting in landscape pairs with specific correlations. We find that switching between these landscape pairs, depending on the ruggedness of the landscape and the interlandscape correlation, can either increase or decrease steady‐state fitness relative to evolution in single environments. In addition, we show that switching between rugged landscapes often selects for increased fitness in both landscapes, even in situations where the landscapes themselves are anticorrelated. We demonstrate that positively correlated landscapes often possess a shared maximum in both landscapes that allows the population to step through sub‐optimal local fitness maxima that often trap single landscape evolution trajectories. Finally, we demonstrate that switching between anticorrelated paired landscapes leads to ergodic‐like dynamics where each genotype is populated with nonzero probability, dramatically lowering the steady‐state fitness in comparison to single landscape evolution.     

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.