Research Article  Open  Published:
The competition between simple and complex evolutionary trajectories in asexual populations
BMC Evolutionary Biologyvolume 15, Article number: 55 (2015)
Abstract
Background
On rugged fitness landscapes where sign epistasis is common, adaptation can often involve either individually beneficial “uphill” mutations or more complex mutational trajectories involving fitness valleys or plateaus. The dynamics of the evolutionary process determine the probability that evolution will take any specific path among a variety of competing possible trajectories. Understanding this evolutionary choice is essential if we are to understand the outcomes and predictability of adaptation on rugged landscapes.
Results
We present a simple model to analyze the probability that evolution will eschew immediately uphill paths in favor of crossing fitness valleys or plateaus that lead to higher fitness but less accessible genotypes. We calculate how this probability depends on the population size, mutation rates, and relevant selection pressures, and compare our analytical results to WrightFisher simulations.
Conclusion
We find that the probability of valley crossing depends nonmonotonically on population size: intermediate size populations are most likely to follow a “greedy” strategy of acquiring immediately beneficial mutations even if they lead to evolutionary dead ends, while larger and smaller populations are more likely to cross fitness valleys to reach distant advantageous genotypes. We explicitly identify the boundaries between these different regimes in terms of the relevant evolutionary parameters. Above a certain threshold population size, we show that the probability that the population finds the more distant peak depends only on a single simple combination of the relevant parameters.
Background
In an adapting population, evolution often has the potential to follow many distinct mutational trajectories. In order to predict how the population will adapt, we must understand how evolution chooses among these possibilities. Many experimental and theoretical studies have analyzed this question, focusing primarily on the simple case where epistasis is absent, so that each mutation has some fixed fitness effect [16]. This work can explain the probability that a given mutation will fix as a population adapts, as a function of its fitness effect, the population size, mutation rate, distribution of fitness effects of other mutations, and other parameters of the evolutionary process.
However, the fitness effect of a mutation often depends on the genetic background in which it occurs. A particularly interesting form of this phenomenon, sign epistasis, occurs when several mutations are individually neutral or deleterious but their combination is beneficial [7]. Sign epistasis has been observed repeatedly in experiments [813], and plays a central role in the evolution of complex phenotypes that involve multiple interacting components. When sign epistasis is present, adaptation can involve passing through genotypes of lower fitness — i.e. a population may have to cross a fitness valley or plateau. Thus the fate of a mutation depends not only on its fitness, but also on its adaptive potential [14].
Several recent theoretical studies have analyzed the evolutionary dynamics of fitness valley crossing [1520]. This work has focused on calculating the rate at which adapting populations cross a valley or plateau, in the absence of any other possible mutational trajectories. However, individually beneficial mutations may often compete with more complex evolutionary trajectories. We must then ask how likely evolution is to eschew the immediately uphill paths, and instead cross valleys or plateaus to reach better but less accessible genotypes. In other words, when the fitness landscape is rugged, we wish to understand whether evolution will take the more “farsighted” path to reach distant advantageous genotypes, rather than a “greedy” trajectory that fixes immediately beneficial mutations regardless of whether these may lead to evolutionary dead ends.
In this article, we analyze this evolutionary choice between immediately beneficial mutations and more complex mutational trajectories that ultimately lead to higher fitness. We calculate the probability that an adapting population will follow each type of competing trajectory, as a function of the population size, mutation rates, and selection pressures. We focus on asexual populations, where the only way for a population to acquire a complex adaptation is for a single lineage to acquire each mutation in turn. Our analysis is similar in spirit to earlier work which also considered the tradeoff between shortterm and longterm fitness advantages [2124]. However, these earlier studies dealt with competition between different strictly uphill or neutral paths, and considered the case where the less beneficial initial mutation led to better longterm evolutionary opportunities. In contrast, our analysis describes the competition between uphill mutations and more complex trajectories. While these two cases can be qualitatively similar in very small populations, they lead to very different dynamics in larger populations where the sign of the effect of the intermediate mutation can play a crucial role.
Our results show that population size has a crucial impact on how “farsighted” evolution can be. This dependence is not monotonic: evolution at intermediate population sizes is most “greedy”, while both larger and smaller populations are more likely to eschew uphill paths in favor of complex trajectories. In large populations, our results show that a single parameter reliably predicts the extent of this evolutionary “foresight” across a wide range of parameters. Finally, we describe how our analysis can be generalized to predict how evolution will choose among even more complex trajectories, such as broad fitness valleys with multiple intermediate genotypes, and we discuss evolution in genotype spaces with many possible evolutionary paths.
Methods
We are interested in how a population makes an evolutionary choice when confronted with multiple possible mutational trajectories. Specifically, we focus on the extent to which adaptation proceeds by crossing fitness valleys rather than acquiring immediately beneficial (uphill) mutations. Of course, the relative frequency of valley crossing will depend on the number of available fitness valleys, their depth, and the fitness advantage of the multiplemutants, as well as the distribution of fitness effects (DFE) of the uphill mutants. Our goal is to understand how the prevalence of valley crossing depends on these factors.
Model
Throughout most of this article, we consider the simplest context in which we can address this question: the choice between a single uphill path and a single fitness valley. Specifically, we consider a haploid asexual population of constant size N which can either acquire an uphill mutation (u) that confers an immediate fitness advantage s _{ u }, or alternatively acquire a deleterious fitness valley intermediate (i) with fitness deficit δ _{ i } on which background a doublemutant (v) with fitness s _{ v }>s _{ u } can arise. This scenario is illustrated in Figure 1. We also consider the case of a fitness plateau, where δ _{ i }=0.
Because we are interested in the evolutionary choice between competing mutational trajectories, we assume that these two trajectories are mutually exclusive (i.e. the mixed genotypes ui and uv are strongly deleterious), so that only one genotype (either u or v) can eventually fix in the population. As a measure of evolutionary foresight, we analyze the probability that the doublemutant v fixes as a function of the relevant mutation rates, selection coefficients, and population size. In some situations, we could imagine that after either genotype u or v fixes, another set of competing potential trajectories become available. In this case, our analysis predicts the longterm relative ratio of fixed uphill versus valleycrossing mutations. In the Discussion, we consider how this model can be extended to the situation where there are many different competing uphill paths and valleys, and to broader fitness valleys involving multiple intermediate genotypes.
Simulations
In addition, we compare our analytical predictions for valley crossing probability to WrightFisher simulations. Each simulated population was evolved until either the uphill genotype or valleycrossing genotype fixed. Valley crossing probabilities were then inferred from the number of trials in which the valleycrossing genotype fixed, out of 1000 trials per parameter set.
Results
In the absence of the uphill genotype, fitness valley crossing can be modeled as a homogeneous Poisson process with rates as calculated by [17]. In small populations, the primary role of the uphill genotype is to introduce an effective time limit on this process: once an uphill mutation destined to survive drift first occurs, it very quickly fixes, leading to the extinction of the wildtype. The probability of valleycrossing can thus be calculated as the probability that the intermediate i fixes before the uphill genotype u. An example of this is shown in Figure 1b.
In larger populations, the dynamics are more complex, as illustrated in Figure 1c. Rather than leading to a single cutoff time for valleycrossing to occur, the singlemutant occurs and gradually increases in frequency. This leads to a decline in the size of the wildtype background on which intermediate and valleycrossing mutants can arise, and a corresponding increase in the mean fitness of the population (Figure 1c). These effects gradually reduce the rate at which intermediates are produced, and make these intermediates effectively more deleterious relative to the mean fitness. These factors reduce the rate of the valleycrossing process. Thus valleycrossing becomes an inhomogenous Poisson process, with a rate that depends on the random appearance time T _{ u } of the uphill mutant.
In general, these effects of interference and tunneling are complex. However, the analysis becomes simpler in two specific regimes. When the expected drift time of the intermediate genotype is short, we can neglect the changing background fitness due to the uphill mutant during this drift time (Figure 1d). Alternatively, for very large populations (N μ>1), the Poisson process approximation breaks down and both uphill and intermediate mutations can be treated deterministically (Figure 1e), and only the valleycrossing genotype must be treated stochastically.
These various regimes are illustrated in Figure 2. We now analyze each in turn, assuming weak selection (s _{ j }<1) for all genotypes throughout. Taken together, this provides a complete picture of the probability that evolution will eschew the immediately uphill path in favor of the more complex adaptation.
Small populations
When the population size is small enough that the probability of stochastic tunneling is very low, the population is generally clonal or nearly clonal, and moves in Markovian jumps between neighboring genotypes. The transition between genotypes i and j occurs at rate
where π _{ i j } is the probability that a single j mutant will give rise to a lineage that fixes, given by the standard formula,
We refer to this as the sequential fixation regime. Because we are considering neutral and weakly deleterious intermediates, we account for the possibility of backmutation to the wild type if the intermediate fixes. Therefore, the process can be modeled as an absorbing states Markov chain, where the wild type and intermediates act as transient states, and the uphill and doublemutant genotypes act as absorbing states. From elementary Markov chain theory, we find
As the population size increases, π _{ wu }→2s _{ u } and π _{ wi }→0, so $\frac {\pi _{\textit {wu}}}{\pi _{\textit {wi}}} \rightarrow \infty $ , and P _{cross}→0. Thus we find that within the sequential fixation regime, larger population sizes are less likely to cross fitness valleys.
Stochastic tunneling
For large populations, the probability that deleterious intermediates will fix declines drastically, and successful double mutants will instead arise on the unfixed singlemutant background in a process known as stochastic tunneling [16]. This transition occurs when
where p _{ v } is the probability that the intermediate lineage survives drift long enough to give rise to an ultimately successful double mutant lineage (we will explicitly calculate this probability below). We can then model the appearance of an intermediate mutant lineage destined to give rise to a double mutant lineage as a Poisson process. The rate λ _{ v } at which these lineages appear is given by the rate at which intermediate mutations arise times the probability of success of the lineage, integrated over the drift time t _{ d } after appearance of the singlemutant intermediate:
Here N _{wt} is the wildtype population size, and p _{ i }(t _{ d }) is the cumulative probability that a singlemutant lineage will give rise to a successful doublemutant lineage by time t _{ d } after it appears. This probability is given by [17]:
where
and δ _{ i,eff} and s _{ v,eff} are the fitnesses of the intermediate and valleycrossing genotypes relative to the (timedependent) mean population fitness. These effective fitnesses and N _{wt} depend on the background at time t _{ d } in a way we must now consider. The background in turn is determined by the frequency f _{ u }(T _{ v }+t _{ d }) of the uphill genotype at time t _{ d } after the appearance of the first singlemutant intermediate lineage destined for success at time T _{ v } (Figure 1c). Thus the uphill genotype frequency sets the fitness background on which valleycrossing probabilities are determined.
To calculate the relevant effective parameters, we note that the appearance of uphill lineages destined for success can be modeled as a Poisson process. Moreover, because the valleycrossing genotypes make up a tiny fraction of the population (unless the doublemutant has already established), we can treat the genetic background on which these uphill lineages appear as essentially fixed. Therefore, the first uphill lineage destined to survive genetic drift will appear at time T _{ u }, distributed exponentially with rate
Once a successful uphill lineage appears, we assume it establishes in time τ _{est}=γ _{ e }/(2s _{ u }), where γ _{ e }≈.577 is the EulerMascheroni constant [25], and then sweeps deterministically according to
where $\hat {t}\equiv t{T_{u}}\tau _{\text {est}}$ is the time after establishment of the uphill mutant.
Conditioning on the appearance time T _{ u }, we can thus work out our effective parameters
These effective parameters encompass the two main effects of the sweeping uphill mutation on the valley crossing probability: the first represents the declining wildtype background on which new mutations can arise, and the remaining two represent the decreasing relative fitness of the valleycrossing lineage.
We are interested in the probability that a doublemutant lineage destined for success appears before the uphill genotype fixes. Integrating over all possible appearance times T _{ u }, this is given by:
This integral is a complete solution for the probability of valleycrossing in the stochastic tunneling regime, provided that the population size is small enough that the Poisson process approximation above holds. Although it does not have a simple closedform solution, we can easily evaluate the integral numerically. Alternatively, there is a simple and relevant parameter regime in which background fitnesses change slowly. We now consider this case, and show that it allows us to evaluate our expression for the valleycrossing probability explicitly. In a later section below, we turn to the alternative case where the Poisson process approximation breaks down, and we can instead treat all singlemutants deterministically; the valleycrossing probability also simplifies considerably in this very large population regime.
Slowly changing background fitness
One of the main complications of Equation (13) is the integral over possible drift times t _{ d }, which reflects the increasing effective deleteriousness of the intermediate genotype as the uphill genotype sweeps to fixation and increases the mean fitness of the background population. However, when s _{ u } is small or δ _{ i } is large, this background fitness changes slowly compared to the intermediate drift time. In this case, we can treat the background during intermediate drift as effectively constant (Figure 1d). This eliminates the need for an integral over t _{ d }, since the timedependent probability of success of a singlemutant at time t is fully determined by f _{ u }(t+t _{ d })≈f _{ u }(t).
We can further simplify the analysis if we treat the probability of crossing the valley as a function of two probabilities: the probability P _{ v,1} that the first successful valleycrossing lineage appears before the first successful uphill lineage establishes, and the probability P _{ v,2} that a successful valleycrossing lineage appears after the uphill mutant establishes:
The calculation of P _{ v,1} takes place on a purely wildtype background, so we can use
where p _{ v }, the probability that the intermediate lineage survives drift long enough to give rise to an ultimately successful double mutant lineage, is given by [17]:
Meanwhile, λ _{ u } is unchanged from the original analysis. P _{ v,1} is determined by a race between these two exponential random variables. Using basic properties of the exponential, the probability that the doublemutant appears before the uphill genotype establishes is therefore
where Δ τ represents the difference between the mean drift time τ _{drift} of the valleycrossing lineage and the mean establishment time τ _{est} of the uphill lineage,
Here τ _{est} is as given above, and from [17], we can approximate the drift time as
If the successful uphill lineage establishes, the first successful valleycrossing lineage still has a chance to appear and outcompete it, albeit on a declining wildtype background. Thus for P _{ v,2} we get a similar integral as in the original analysis. However, since we are approximating p _{ i } as constant, the rate of successful lineage generation simplifies to
Integrating and assuming mutation rates are small compared to selection pressures, we find:
Combining these results, we find
We expect this result to be valid provided that δ _{ i } is effectively constant over the expected drift time. This will hold when
Semideterministic tunneling
We now consider the case where N μ _{ i }>1 and N μ _{ u }>1, and hence the Poisson process approximation used to derive Equation (13) breaks down. Fortunately, in this regime the number of singlemutant intermediates and uphill mutants in the population are well approximated by their deterministic expectation (Figure 1e). Thus the only random variable is the appearance time of the first successful doublemutant lineage, which occurs with rate
Because intermediates never make up a large portion of the population, f _{ u } is unaffected by f _{ i }, and hence is still given by Equation (9). We can then approximate the frequency f _{ i } of the singlemutant intermediates using mutationselection balance with a declining wildtype population:
where f _{ i } ^{∗} gives the independent deterministic dynamics (mutationselection balance) of single mutants on the wildtype background, and (1−f _{ u }) is the size of the wildtype background. It is useful to transform this into an integral in the frequency domain:
We note from equation 9 that
We further approximate
Combining these expressions and assuming μ _{ u }≪s _{ u }, we find
where
This gives
where we have defined the useful quantity
Thus we see that in very large populations, the probability of valleycrossing depends in a simple way on the single composite parameter Γ. The form of this composite parameter depends crucially on whether the fitness cost of the intermediate genotype is large or small, as defined in Equation (31).
Discussion
Our results have shown that the size of a population strongly influences how “farsighted” it can be. In small populations, genetic drift is strong relative to selection, so the evolutionary dynamics proceeds by sequential fixation. Since fixation of a deleterious intermediate becomes less likely in larger populations, this means that increasing the population size initially decreases the relative influence of fitnessvalley crossing. However, as the population size increases further, beneficial mutations take longer to fix, maintaining diversity in the population and allowing double mutants to stochastically tunnel on the declining wildtype background [1517]. Together, these effects lead to a nonmonotonic relationship between population size and the probability that evolution will favor the complex adaptation over the directly uphill path, as illustrated in Figure 3a. This nonmonotonic dependence on population size is similar in spirit to the results of earlier work analyzing evolution on epistatic landscapes in the absence of fitness valleys [2123].
It is interesting to note that P _{cross} does not immediately begin to rise with the onset of tunneling. Instead, the dependence is more complex, as a consequence of the tradeoff between increasing mutation rates and fixation times. Nevertheless, for populations in which the transition to valleycrossing behavior occurs in the semideterministic regime, we can derive a simple expression for the threshold size at which the population will tend to cross valleys with probability P _{cross}. A straightforward inversion of (32) gives
valid for large population sizes in the semideterministic regime. Thus in this regime the threshold size above which a population exhibits a given degree of foresight (i.e. has a particular P _{cross}) depends only on Γ. To illustrate this, in Figure 3b we show P _{cross} as a function of N Γ for a variety of simulations across different values of μ _{ u },μ _{ i },μ _{ v },s _{ u },δ _{ i } and s _{ v }. It is clear that even across a wide parameter range, N Γ is a reliable predictor of valley crossing probability in the semideterministic limit.
Multiple intermediates and evolutionary predictability
Recently, Szendro et al. [26] simulated evolution across a wide range of population sizes on an experimentallyderived epistatic fitness landscape, finding that “evolutionary entropy” (i.e. unpredictability over all possible outcomes, given by $S =  \sum p_{j} \log p_{j}$ ) varied nonmonotonically with population size. Specifically, these authors found that entropy initially decreased with N above a characteristic population size N∝1/μ, before increasing again above a second characteristic size N∝1/μ ^{2}; they argued that these points were related to the supply rate of single and double mutants respectively. Our analysis is consistent with these results. For example, the increase in entropy at N∝1/μ ^{2} found by [26] corresponds to our result that valleycrossing begins to significantly influence evolution when N∼1/Γ: this is approximately proportional to 1/μ ^{2}, albeit with an additional log dependence on μ from the γ factor that would be harder to observe experimentally.
We find related behavior if we extend our analysis to valleys with more intermediates. As a simple example, we consider a fitness landscape where we add deleterious intermediates u _{0} and v _{0} (each occurring at rate μ _{0} and leading to fitness cost δ _{ i }) to the uphill and valleycrossing branches respectively, so that we now have competition between a singleintermediate valley and a twointermediate valley. In a large enough population, mutationselection balance will ensure that
If we assume that these subpopulations are large enough that doublemutants behave deterministically, then we find the crossing probability obeys
Thus our analysis of valleys with multiple intermediates suggests that the N μ ^{2} entropy peak is not unique: as the population grows larger, there should be entropy peaks corresponding to foresight across valleys with increasing numbers of intermediates. The emergence of such second peaks has been observed in simulations [23], and our model offers a quantitative outline of where such peaks should occur given the relevant evolutionary parameters. In the example above, for instance, there would be an entropy peak approximately proportional to N μ ^{3}, and in general, we could expect entropy peaks at points approximately proportional to N μ ^{n} for (n−1)intermediate valleys. In practice, however, the semideterministic approximation will break down for any sizable number of intermediates, unless the population size is unrealistically large.
Many paths
Throughout this paper, we have assumed the presence of a single uphill mutation and fitness valley. We now consider how our analysis can be extended to predict how evolution chooses among many such possible mutational trajectories.
In small populations that are in the sequential fixation regime, we simply add additional transient transition matrix elements representing different mutations, with the uphill mutations transitioning to the uphill absorbing state, and similarly for the valleycrossing mutations. When stochastic tunneling is important, we must instead add the rates of single valleycrossing mutants to get a total rate of
and similarly find the total rate of uphill mutants that are destined to survive drift,
where in the last equality we have replaced a discrete collection of uphill mutations with a continuous distribution of uphill fitness effects ρ(s), and a discrete collection of mutation rates μ _{ u } with a total beneficial mutation rate U _{ b }. This is valid as long as Λ _{ u }≪1. Once an uphill mutation destined to survive drift occurs, the probability that it has fitness s is given by the ratio between its partial rate and the total rate; formally, the probability density is given by:
Using these expressions, we can integrate our results from the analysis over all possible trajectories. However, we note that if there are a large number of weakly beneficial mutations, it is possible the first successful lineage to appear will be outcompeted by a stronger uphill mutation that arises later but fixes first. Our analysis applies provided that we consider only uphill mutations that reach a significant portion k of the population before a new, more fit uphill mutant is expected to be produced: i.e.
where τ _{ k } is the expected time for a singlemutant destined for success to make up frequency k of the population. This is consistent with our intuition that as the population size grows larger, we increasingly expect the mutations of largest effect to dominate the dynamics.
Conclusions
Using a simple threelocus fitness landscape model (Figure 1a), we identified several regimes of valleycrossing with qualitatively different behavior (Figure 2). By examining the behavior in each of these regimes in turn, we found that the probability of valleycrossing has a complex, nonmonotonic dependence on population size (Figure 3a), and identified a parameter Γ that reliably predicts the population size at which valley crossing becomes preferred (Figure 3b). Finally, we showed how these results can be extended to fitness valleys with more intermediates, and to fitness landscapes with many possible evolutionary trajectories, as is the case in most naturally occurring populations.
Abbreviations
 DFE:

Distribution of fitness effects
References
 1
Perfeito L, Fernandes L, Mota C, Gordo I. Adaptive mutations in bacteria: High rate and small effects. Science. 2007; 317(5839):813–5.
 2
Tenaillon O, RodríguezVerdugo A, Gaut RL, McDonald P, Bennett AF, Long AD, et al.The molecular diversity of adaptive convergence. Science. 2012; 335(6067):457–61.
 3
Lang GI, Rice DP, Hickman MJ, Sodergren E, Weinstock GM, Botstein D, Desai MM. Pervasive genetic hitchhiking and clonal interference in forty evolving yeast populations. Nature. 2013; 500(7464):571–4.
 4
Gerrish P, Lenski R. The fate of competing beneficial mutations in an asexual population. Genetica. 1998; 102:127–44.
 5
Schiffels S, Szollosi GJ, Mustonen V, Lassig M. Emergent neutrality in adaptive asexual evolution. Genetics. 2011; 189(4):1361–75.
 6
Good BH, Rouzine IM, Balick DJ, Hallatschek O, Desai MM. Distribution of fixed beneficial mutations and the rate of adaptation in asexual populations. Proc Nat Acad Sci. 2012; 109:4950–5.
 7
Weinreich DM, Watson RA, Chao L. Perspective: Sign epistasis and genetic constraint on evolutionary trajectories. Evolution. 2005; 59(6):1165–74.
 8
Weinreich DM, Delaney NF, DePristo MA, Hartl DL. Darwinian evolution can follow only very few mutational paths to fitter proteins. Science. 2006; 312(5770):111–4.
 9
Poelwijk FJ, Kiviet DJ, Weinreich DM, Tans SJ. Empirical fitness landscapes reveal accessible evolutionary paths. Nature. 2007; 445:383–6.
 10
Kvitek DJ, Sherlock G. Reciprocal sign epistasis between frequently experimentally evolved adaptive mutations causes a rugged fitness landscape. PLoS Genet. 2011; 7(4):1002056.
 11
Silva RF, Mendonça SC, Carvalho LM, Reis AM, Gordo I, Trindade S, et al.Pervasive sign epistasis between conjugative plasmids and drugresistance chromosomal mutations. PLoS Genet. 2011; 7(7):1002181.
 12
Dawid A, Kiviet DJ, Kogenaru M, de Vos M, Tans SJ. Multiple peaks and reciprocal sign epistasis in an empirically determined genotypephenotype landscape. Chaos: Interdisciplinary J Nonlinear Sci. 2010; 20(2):026105.
 13
Salverda ML, Dellus E, Gorter FA, Debets AJ, Van Der Oost J, Hoekstra RF, et al.Initial mutations direct alternative pathways of protein evolution. PLoS Genet. 2011; 7(3):1001321.
 14
Woods RJ, Barrick JE, Cooper TF, Shrestha U, Kauth MR, Lenski RE. Secondorder selection for evolvability in a large escherichia coli population. Science. 2011; 331(6023):1433–6.
 15
Weinreich DM, Chao L. Rapid evolutionary escape by large populations from local fitness peaks is likely in nature. Evolution. 2005; 59:1175–82.
 16
Iwasa Y, Michor F, Nowak MA. Stochastic tunnels in evolutionary dynamics. Genetics. 2004; 166(3):1571–9.
 17
Weissman DB, Desai MM, Fisher DS, Feldman MW. The rate at which asexual populations cross fitness valleys. Theor Population Biol. 2009; 75(4):286–300.
 18
Gokhale CS, Iwasa Y, Nowak MA, Traulsen A. The pace of evolution across fitness valleys. J Theor Biol. 2009; 259(3):613–20.
 19
Altland A, Fischer A, Krug J, Szendro IG. Rare events in population genetics: stochastic tunneling in a twolocus model with recombination. Phys Rev Lett. 2011; 106(8):088101.
 20
Weissman DB, Feldman MW, Fisher DS. The rate of fitnessvalley crossing in sexual populations. Genetics. 2010; 186(4):1389–410.
 21
Rozen DE, Habets MG, Handel A, de Visser JAG. Heterogeneous adaptive trajectories of small populations on complex fitness landscapes. PLoS One. 2008; 3(3):1715.
 22
Handel A, Rozen DE. The impact of population size on the evolution of asexual microbes on smooth versus rugged fitness landscapes. BMC Evolutionary Biol. 2009; 9(1):236.
 23
Jain K, Krug J, Park SC. Evolutionary advantage of small populations on complex fitness landscapes. Evolution. 2011; 65(7):1945–55.
 24
Van Nimwegen E, Crutchfield JP. Metastable evolutionary dynamics: crossing fitness barriers or escaping via neutral paths?. Bull Math Biol. 2000; 62(5):799–848.
 25
Desai MM, Fisher DS. Beneficial mutation–selection balance and the effect of linkage on positive selection. Genetics. 2007; 176(3):1759–8.
 26
Szendro IG, Franke J, de Visser JAG, Krug J. Predictability of evolution depends nonmonotonically on population size. Proc Nat Acad Sci. 2013; 110(2):571–6.
Acknowledgements
We thank Benjamin Good, Sergey Kryazhimskiy, Elizabeth Jerison, and other members of the Desai lab for useful discussions, and Katya Kosheleva for help with the figures. This work was supported by the Harvard HCRP, Harvard Herchel Smith, and Harvard PRISE programs (I.E.O.), and the James S. McDonnell Foundation, the Alfred P. Sloan Foundation, the Harvard Milton Fund, grant PHY 1313638 from the NSF, and grant GM104239 from the NIH (M.M.D.). Simulations in this paper were performed on the Odyssey cluster supported by the Research Computing Group at Harvard University.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
IEO and MMD designed the research, conducted the analysis, and wrote the paper. Both authors read and approved the final manuscript.
Rights and permissions
About this article
Received
Accepted
Published
DOI
Keywords
 Epistasis
 Rugged fitness landscape
 Fitness valley