# Fixation probability of a beneficial mutation conferring decreased generation time in changing environments

- Fangshu Cui
^{1}Email author and - Bo Yuan
^{1}

**12(Suppl 4)**:48

https://doi.org/10.1186/s12918-018-0575-9

© The Author(s). 2018

**Published: **24 April 2018

## Abstract

### Background

One central building block of population genetics is the fixation probability. It is a probabilistic understanding of the eventual fate of new mutations. Moreover, the fixation probability of new beneficial mutations plays an important effect on the adaptation of populations to environmental challenges. Great progress has been made in the study of the beneficial mutations that increases offspring number. However, the fixation probability of beneficial mutations with a shorter generation time under various genetic and ecological conditions has not been explored.

### Results

Here we extend the classical result of the fixation probability of beneficial mutations obtained by Haldane, and estimate the fixation probability of a beneficial mutation with a reduced generation time in a changing environment. Assuming that the selective advantage is very small, we concentrate all the changing factors of environment on a single quantity: effective selective advantage. Using a time-dependent branching process, we get the analytic approximation for the fixation probability of beneficial mutations that decrease the generation time. Then, we apply this approximation to four interesting biological cases.

### Conclusions

In these instances, we show the comparison of the approximation with the accurate values. We find that they are consistent, demonstrating the effectiveness of our result for the fixation probability of beneficial mutations conferring a reduced replication time.

## Keywords

## Background

One cornerstone of population genetics is the fixation probability, i.e. the probability that mutations survive loss. This is a probabilistic comprehension of the eventual fate of the beneficial, neutral, or deleterious mutations. Furthermore, the fixation probability of new beneficial mutations has a significant impact on the rate of adaptation of populations [1, 2]. When a new beneficial mutation goes into a population, to finish an adaptive step, it has to escape random loss due to genetic drift, rise to enough copies, and eventually get fixed. Actually, the frequency of a beneficial mutation fluctuates over time. When this frequency is low, the genetic drift is likely to lead to the disappearance of the beneficial mutation, which needs a stochastic process. Once this frequency is large enough, a deterministic model can be used to well approximate the further increment of this frequency. Usually, the calculation for the frequency of a beneficial mutation is equivalent to evaluating the probability that the beneficial mutation survives an earliest stage of strong genetic drift.

Since 1920s, interest in the computation of fixation probabilities has been maintained for nearly one century and considerable progress has been made in this problem. Generally, there are three methods to estimate the fixation probabilities: Markov chain, branching process and diffusion approximation. When the individuals and genotypes in a population can be enumerated, the Markov chain method can obtain the fixation probability precisely. Therefore, this method is characteristically practicable only when the population size is very small [3, 4]. Once the population size becomes large, the discrete branching processes are in wide use [5–13]. Since the branching process method presumes that the population is large enough that the destiny of each mutation is independent of all others, it gives an approximation to the real fixation probability. When the selection is weak in a large population, the diffusion approximation approaches are usually used [9, 14–16]. Furthermore, many literatures have tried to integrate and reconcile the discrete and continuous methods [8, 9, 17, 18].

Classically, if a native wild-type individual has on average one offspring per generation, a beneficial mutant has on average 1 + *s* offspring per generation, where the parameter *s* (*s* > 0) is the selective advantage. This mechanism of the selective advantage is defined as fecundity, which is fundamental in a large amount of literature in population genetics [19–22]. Assuming a Poisson offspring distribution and a small, constant selection coefficient in a population of constant size, Haldane [20] gives the well-known result that the fixation probability is approximately 2*s*, for a mutation that increases fecundity.

Nevertheless, the mutants in many organisms may produce the same number of offspring as the wild-type in a shorter generation time: so-called “generation time” mutants. For example, in the bacteria population, a mutant that has the antibiotic resistance completes the cell cycle and produces two offspring faster than the drug-sensitive individuals. In this case, a reduced replication time is obviously a more suitable mechanism for the selective advantage. Wahl and DeHaan [12] have firstly demonstrated that the classic approximation 2*s* for the fixation probability of a beneficial mutation does not hold for this mutation conferring a decreased development time. Using a model with a Poisson-distributed offspring with mean 2 and a weak, constant selective advantage, they have approximated the fixation probability of this “generation time” mutant as *s*/ln(2) for a population of constant size. Therefore, if all mutations are assumed to increase the offspring number, it leads to an overestimate of the order 2 ln(2) for the fixation probability of the mutation that reduces the replication time.

The study of fixation probability under diverse genetic factors and ecological scenarios has been explored [23]. A series of articles have estimated the survival probability of beneficial mutations when the population size changes [3, 4, 8, 9, 18, 24]. Ewens [24] derived the fixation probability of a beneficial mutation in two cases of changing population sizes: a cyclic sequence of population sizes and the population size that first increases and then remains constant. Otto and Whitlock [8] studied the survival probability of beneficial mutants under several demographic models of population size change, including a single change, exponential growth or decline, logistic growth or decline, and fluctuating size. Wahl and Gerrish [9] examined the influence of population bottlenecks on the fixation probability. Lambert [18] and Parsons and Quince [3, 4] developed the fixation probability of beneficial mutations when the population size changed dynamically. Since Pollak [25] first studied the fixation probability in a subdivided population, great process has been made in this probability in spatially heterogeneous populations [16, 26–28]. The influence of linked loci on the fixation probability of an advantageous allele has been widely investigated [29–32]. Studies on time-dependent selection mainly focus on random fluctuations of selection coefficients [33–36]. Recently, Waxman [37] and Uecker and Hermisson [38] addressed the question of the establishment of new beneficial mutants when the change of selection coefficients and population sizes follows an explicit trend. Peischl and Kirkpatrick [39] derived analytical approximations for the fixation probability of favorable mutations in arbitrarily changing environment that used a novel approach assuming small environmental fluctuations. In these studies, the beneficial mutations are assumed to increase the average number of offspring. However, when a mutant confers an advantage in generation time, the effects of various genetic and ecological conditions on the fixation probability have not received enough attention.

In this article, we use a time-dependent branching process to study the fixation probability of “generation time” mutants in changing environments. Assuming weak selection, we centralize all the environmental changes into a single parameter: effective selective advantage, and deduce an analytical approximation for this fixation probability. We apply our result to four absorbing biological cases, including the monotone increase and the periodic change of the selective advantage in a population of constant size, the changing population size, and the stochastic fluctuations in selection. In these conditions, our approximation compares well with the numerical calculation, which demonstrates the effectiveness of our result.

## Methods

### Branching process

In probability theory, the branching process is a mathematical object known as stochastic process. It is used to model reproduction, that is, to model a population in which each individual produces stochastic number of offspring in the next generation. It can also be used to model some other similar dynamics, for example, the dispersion of surnames in genealogy, the spread of neutrons in a nuclear reactor and so on.

*k*(

*k*= 0, 1, 2, ...) offspring with the assigned probability

*f*

_{ k }in the next generation. Then, these offspring have the same reproductive capacity with the ancestor, that is, each of them produce

*k*offspring with the probability

*f*

_{ k }, and so on. The crucial hypothesis of the branching process model is that the offspring distribution of each individual is independent and identical. However, only when the mutant lineage is a small part of the population, this hypothesis is true. If the population size is constant and the mutant lineage becomes a large part of the population, the density of each individual must rely on the others.

### Fixation probability of beneficial mutations in a constant environment

*r*offspring on average per generation, this individual will have

*r*

^{ t }offspring and a “fecundity” mutant having selective advantage

*s*will have (

*r*(1 +

*s*))

^{ t }offspring in

*t*generations. In the “growth” stage of a population, a lot of offspring are reproduced. Due to the assumption of constant population size, only part of offspring can survive in the next generation, which is called the “sampling” stage. We use

*y*(

*x*) to be the probability generating function that represent the number of offspring reproduced by a “fecundity” mutant lineage in one generation, that is, the probability generating function in the “growth” stage. The

*y*(

*x*) is described as follows: if the probability that an individual produces

*i*offspring is

*p*

_{ i },

*y*(

*x*) =

*p*

_{0}+

*p*

_{1}

*x*+

*p*

_{2}

*x*

^{2}+ ... [40]. Assuming that the offspring obey the Poisson distribution and the average number of offspring of a wild-type individual is

*r*= 2, a mutant with a fecundity advantage

*s*has 2(1 +

*s*) offspring on average per generation. Hence,

*N*to denote the population size, which involves both the wild type and the mutant. Because the population size is constant, only

*N*of the newly produced offspring will survive in every generation. Owing to the supposition that the average number of offspring of each wild-type individual is

*r*, the probability that each offspring survives is 1/

*r*. Conversely, the probability that each offspring dies is 1 − 1/

*r*. Therefore, the probability generating function in the “sampling” stage is

*z*(

*x*) = 1 − 1/

*r*+ (1/

*r*)

*x*. According to the property of the probability generating function, the overall probability generating function of “growth” stage followed by “sampling” stage is

*ϕ*(

*x*) =

*y*(

*z*(

*x*)) =

*y*∘

*z*(

*x*), which is the total probability generating function for the number of offspring in a “fecundity” mutant lineage after one generation. Consequently, the whole probability generating function for the “fecundity” mutant lineage after

*n*generation is

*x*= 0, we can eliminate the higher terms and obtain the probability that the mutant lineage eventually disappear. Hence, the extinction probability of a “fecundity” mutation is \( q=\underset{n\to \infty }{\lim }{\phi}_n(0) \). And, the fixation probability of the “fecundity” mutation is [20].

*τ*generations, where

*τ*is constant. The “fecundity” mutant lineage will experience

*τ*sequential “growth” stages and one “sampling” stage, so the total probability generating function is

Here, the *ϕ*(*x*) can be simply denoted as *ϕ*(*x*) = *y*_{
τ
} ∘ *z*(*x*). Particularly, the situation that the bottlenecks occur every *τ* = 1 generations is equal to a constant population size.

*t*

_{ g }(

*t*

_{ g }< 1) be the generation time of the “generation time” mutant lineage. In

*t*wild-type generations, a “generation time” mutant will produce \( {r}^{t/{t}_g} \) offspring on average. Assuming that the whole growth rate of both the “fecundity” mutant lineage and the “generation time” mutant lineage is the same, we have \( {r}^{t/{t}_g}={\left(r\left(1+s\right)\right)}^t \) and obtain

Let \( \tilde{s}={\log}_r\left(1+s\right) \), then \( {t}_g=1/\left(1+\tilde{s}\right) \). When *s* is very small, we have \( \tilde{s}\approx s/\ln (r) \).

*τ*+ 1 generations between two population bottlenecks. If a mutation of this type firstly appears at the start of a “growth” stage, the situation above will firstly happen before

*n*

_{1}population bottlenecks. Therefore, we have

Extending the Eq. (6), the situation that the “generation time” mutant lineage undergo an extra generation occurs before *n*_{
i
} consecutive population bottlenecks, we find \( \left({n}_i\tau +i\right)\left[1/\left(1+\tilde{s}\right)\right]\le {n}_i\tau \) and have \( {n}_i\ge i/\left(\tilde{s}\tau \right) \).

*ϕ*

^{+}(

*x*) =

*y*

_{τ + 1}∘

*z*(

*x*). For all the other “growth” and “sampling” stages, the probability generating function is

*ϕ*(

*x*) =

*y*

_{ τ }∘

*z*(

*x*). For instance, if

*τ*= 5, \( \tilde{s}=1/20 \), we obtain

*n*

_{1}≥ 4. Thus, the probability generating function of this mutant lineage after eight population bottlenecks is

Accordingly, the fixation probability of the “generation time” mutation can be computed by the Eq. (3).

### Fixation probability of beneficial mutations in a variable environment

*s*

_{ k }(

*s*

_{ k }> 0, and

*s*

_{ k }≪ 1,

*k*= 1, 2, 3,...) changes in time in a changing environment. Extending the classical result for the fixation probability of beneficial mutations in a constant environment, Peischl and Kirkpatrick [39] have provided this probability in a time-dependent branching process:

*s*

_{ e }is the effective selective advantage.

*k*in a variable environment as:

*ε*

_{ k },

*k*= 1, 2, 3, ... is the disturbance function in generation

*k*.

*ε*

_{ k }is a smooth and bounded function that maps [0,1] to [− 1,1], and

*ε*

_{ k }(1) = 0,

*k*= 1, 2, 3, .... Assuming that the offspring obey the Poisson distribution and its instantaneous variation is very small in the changing environment, we have max

_{x, k}[

*ε*

_{ k }(

*x*)] ≪ 1, \( {\max}_{x,k}\left[{\varepsilon}_k^{\prime }(x)\right]\ll 1 \). We substitute the Eq. (11) into the Eq. (9) and expand in a Taylor series

Here, \( \sigma ={\overline{\phi}}^{\prime}\left(1-\overline{p}\right) \), and *O*(*ε*^{2}) denotes the order of max_{x, k}[*ε*_{
k
}^{2}(*x*)] and \( {\max}_{x,k,l}\left[{\varepsilon}_k(x){\varepsilon}_l^{\prime }(x)\right] \). The average number of offspring in the reference environment is more than 1, so we have 0 < *σ* < 1. Because *ε*_{
k
} is bounded, the \( {\sum}_{k=0}^{\infty }{\varepsilon}_k\left(1-\overline{p}\right){\sigma}^k \) in the Eq. (12) is convergent for every series of environments.

*s*

_{ k }of the “generation time” mutation in variable environments is very small, i.e.

*s*

_{ k }≪ 1,

*k*= 1, 2, 3, ..., and \( \overline{s}\ll 1 \), we have

*x*− 1), we obtain an approximation of the disturbance function

*ε*

_{ k }

*ω*

_{ k }}can be understood as a probability distribution. And, the value of

*ω*

_{ k }diminished over time. Here we define the effective selective advantage

*s*

_{ e }as:

*s*

_{ e }can be regarded as a weighted average. Consequently, in a variable environment where the selective advantage changes over time, an analytical approximation for the fixation probability of “generation time” mutations is given by

## Results

We apply our analytic approximation to four interesting biological instances of changing environments. In these cases, we compare our approximation for the fixation probability of “generation time” mutations (Eq. (18)) to the exact value acquired by numerical iteration of Eq. (9).

### Monotonously increasing selection

*s*

_{0}to

*s*

_{∞}in a population of constant size:

*s*

_{∞}≪ 1, we have \( {e}^{-{s}_{\infty }}\approx 1-{s}_{\infty } \). The effective selective advantage acquired by the Eq. (17) is

### Cyclically changing selection

*s*

_{ mean }is the mean of selective advantage,

*Δs*is the amplitude of fluctuations,

*ρ*is a parameter that controls the length of a circle of fluctuations, and

*θ*is used to decide the initial selective advantage.

### Varying population size

*s*is constant in a population in which the population size changes over time. The average number of offspring of each wild-type individual in generation

*n*is

*N*

_{n + 1}/

*N*

_{ n }, and that of each mutant is

*N*

_{ n }is the population size of the wild-type individuals in generation

*n*. We use the solution of the Beverton-Holt equation as the demographic dynamics of the wild-type individuals:

Here, *N*_{0} denotes the initial population size, *g* denotes the growth rate of the population, and *C* denotes the carrying capacity.

### Stochastic fluctuations in selection

*s*

_{ mean }is the mean of selective advantage,

*λ*∈ [0, 1] is the correlation coefficient between

*s*

_{k − 1}and

*s*

_{ k }, and

*ξ*

_{ k }is a white noise that has mean 0 and variance

*σ*

^{2}.

*s*

_{0}, it becomes time dependent and the expected selective advantage in generation

*k*is

## Discussion

In the nature, if an organism carries a beneficial mutation, it will propagate more quickly than the wild-type individuals. In a large number of literature, this means that the beneficial mutant can reproduce more surviving offspring in each generation, i.e. “fecundity” mutant. Nevertheless, a great deal of organisms may produce the same number of offspring as the wild-type in a shorter generation time, i.e. “generation time” mutant. In many cases, the mechanism of the beneficial mutations is very important. Therefore, we studied the fixation probability of “generation time” mutations in variable environments that has not been explored.

The adaptation of populations to environmental challenges relies on the fixation probability of new beneficial mutations. There are a variety of origins of environmental change, and variations can happen in all time scales, from transient variations to transitions in the geological time scales. However, a majority of researches of adaptation depend on the detachment of time scales in the evolution and ecology. This is maybe a serious simplification in many situations. Thus, we explored the establishment probability of beneficial mutations in changing environments.

In a population, individual alleles can undergo changing selection pressures even though the outside environment is constant. This is because multiple selected alleles separate and interfere owing to the linkage or the epistasis [41]. There are many examples, such as the evolution of the compensatory mutation, clonal interference, adaptive genes that flow across a genetic barrier, and so on.

In population genetics, the deduction of the fixation probability from branching processes is a common method. Using branching processes, it is easy to obtain simple analytic results, which are precise when the population size is large enough that *Ns* ≫ 1. The shortcoming of this method is that it is only applicable to beneficial and originally rare mutations. So far, branching processes have been also applied to some other natural phenomenon, including the spread of communicable diseases, the increment of tumor cells and so on [42, 43].

To show how our analytic approximation can be applied to special instances of variable environments, we used it to a few biological cases. They are the monotonously increasing selection, periodically changing selection, changing population size and random fluctuations of selection. In these examples, in order to guarantee the positive result of Eq. (9), we focused on the scenarios that the selective advantage *s*_{
k
} > 0. Generally, our approximations match well with the accurate values (see Figs. 2-5). Nevertheless, our result underestimates the fixation probability of “generation time” mutations if \( {s}_k\ll \overline{s} \) in initial generations (see Fig. 5).

The theoretical framework can be tested by the experiments, such as the recent evolution experiments by Bell and Gonzalez [44]. The results from these experiments compare well with the analytic forecasts on evolutionary rescue [45]. For the microbial population in these experiments, the environment they live in can be manipulated by an automatic liquid processing system. Therefore, this system can be applied to verify the theory of adaptation in variable environments. For the future work, we hope that the fixation probability can be further explored both theoretically and experimentally.

## Conclusions

In this article, we expanded the classical result of the fixation probability of beneficial mutations acquired by Haldane, and calculated an approximation for the fixation probability of “generation time” mutations in a variable environment. When the selective advantage is weak, all the environmental changes are condensed into a unitary quantity: effective selective advantage. This parameter is a weighted mean across the selective advantage per generation, and the weights diminish monotonously over time. Consequently, this fixation probability in changing environments is decided by the environments that the population experienced and the average influence.

We employed our result to four attractive biological cases, which are the monotone increase of selection, the periodic change of selection, varying population size, and stochastic fluctuations in selection. In these situations, our approximation is in good accordance with the precise value, which certifies the effectiveness of our result.

## Declarations

### Acknowledgements

Not applicable.

### Funding

The publication cost of this article was funded by an NNSFC grant (#14Z103010221).

### Availability of data and materials

Not applicable.

### About this supplement

This article has been published as part of *BMC Systems Biology* Volume 12 Supplement 4, 2018: Selected papers from the 11th International Conference on Systems Biology (ISB 2017). The full contents of the supplement are available online at https://bmcsystbiol.biomedcentral.com/articles/supplements/volume-12-supplement-4.

### Authors’ contributions

FC carried the research and writing. BY oversaw the total project, and participated in research and writing. Both authors read and approved the final manuscript.

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Gillespie JH. Genetic drift in an infinite population: the pseudohitchhiking model. Genetics. 2000;155(2):909–19.PubMedPubMed CentralGoogle Scholar
- Orr HA. The rate of adaptation in asexuals. Genetics. 2000;155(2):961–8.PubMedPubMed CentralGoogle Scholar
- Parsons TL, Quince C. Fixation in haploid populations exhibiting density dependence I: the non-neutral case. Theor Popul Biol. 2007;72(1):121–35.View ArticlePubMedGoogle Scholar
- Parsons TL, Quince C. Fixation in haploid populations exhibiting density dependence II: the quasi-neutral case. Theor Popul Biol. 2007;72(4):468–79.View ArticlePubMedGoogle Scholar
- Athreya KB. Rates of decay for the survival probability of a mutant gene. J Math Biol. 1992;30(6):577–81.View ArticlePubMedGoogle Scholar
- Haccou P, Iwasa Y. Establishment probability in fluctuating environments: a branching process model. Theor Popul Biol. 1996;50(3):254–80.View ArticlePubMedGoogle Scholar
- Lange K, Fan RZ. Branching process models for mutant genes in nonstationary populations. Theor Popul Biol. 1997;51(2):118–33.View ArticlePubMedGoogle Scholar
- Otto SP, Whitlock MC. The probability of fixation in populations of changing size. Genetics. 1997;146(2):723–33.PubMedPubMed CentralGoogle Scholar
- Wahl LM, Gerrish PJ. The probability that beneficial mutations are lost in populations with periodic bottlenecks. Evolution. 2001;55(12):2606–10.View ArticlePubMedGoogle Scholar
- Johnson T, Gerrish PJ. The fixation probability of a beneficial allele in a population dividing by binary fission. Genetica. 2002;115(3):283–7.View ArticlePubMedGoogle Scholar
- De Oliveira VM, Campos PRA. Dynamics of fixation of advantageous mutations. Physica A. 2004;337(3–4):546–54.View ArticleGoogle Scholar
- Wahl LM, DeHaan CS. Fixation probability favors increased fecundity over reduced generation time. Genetics. 2004;168(2):1009–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Champagnat N, Lambert A. Evolution of discrete populations and the canonical diffusion of adaptive dynamics. Ann Appl Probab. 2007;17(1):102–55.View ArticleGoogle Scholar
- Yamazaki T. The effects of overdominance of linkage in a multilocus system. Genetics. 1977;86(1):227–36.PubMedPubMed CentralGoogle Scholar
- Gavrilets S, Gibson N. Fixation probabilities in a spatially heterogeneous environment. Popul Ecol. 2002;44(2):51–8.View ArticleGoogle Scholar
- Whitlock MC. Fixation probability and time in subdivided populations. Genetics. 2003;164(2):767–79.PubMedPubMed CentralGoogle Scholar
- Kimura M, Ohta T. Probability of fixation of a mutant gene in a finite population when selective advantage decreases with time. Genetics. 1970;65(3):525–34.PubMedPubMed CentralGoogle Scholar
- Lambert A. Probability of fixation under weak selection: a branching process unifying approach. Theor Popul Biol. 2006;69(4):419–41.View ArticlePubMedGoogle Scholar
- Haldane JBS. The causes of evolution. New York: Harper and Brothers; 1932.Google Scholar
- Haldane JBS. A mathematical theory of natural and artificial selection. Math Proc Camb Philos Soc. 1927;3:607–15.View ArticleGoogle Scholar
- Kimura M. Some problems of stochastic processes in genetics. Ann Math Stat. 1957;28(4):882–901.View ArticleGoogle Scholar
- Kimura M. On the probability of fixation of mutant genes in a population. Genetics. 1962;47(6):713–9.PubMedPubMed CentralGoogle Scholar
- Patwa Z, Wahl LM. The fixation probability of beneficial mutations. J R Soc Interface. 2008;5(28):1279–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Ewens WJ. The probability of survival of a new mutant in a fluctuating environment. Heredity. 1967;22(3):438–43.View ArticleGoogle Scholar
- Pollak E. On the survival of a gene in a subdivided population. J Appl Prob. 1966;3(1):142–55.View ArticleGoogle Scholar
- Whitlock MC, Gomulkiewicz R. Probability of fixation in a heterogeneous environment. Genetics. 2005;171(3):1407–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Gordo I, Campos PRA. Adaptive evolution in a spatially structured asexual population. Genetica. 2006;127(1-3):217–29.View ArticlePubMedGoogle Scholar
- Perfeito L, Pereira MI, Campos PRA, Gordo I. The effect of spatial structure on adaptation in Escherichia coli. Biol Lett. 2008;4(1):57–9.View ArticlePubMedGoogle Scholar
- Barton NH. The reduction in fixation probability caused by substitutions at linked loci. Genet Res. 1994;64(03):199–208.View ArticleGoogle Scholar
- Barton NH. Linkage and the limits to natural selection. Genetics. 1995;140(2):821–41.PubMedPubMed CentralGoogle Scholar
- Johnson T, Barton NH. The effect of deleterious alleles on adaptation in asexual populations. Genetics. 2002;162(1):395–411.PubMedPubMed CentralGoogle Scholar
- Takahasi KR, Tajima F. Evolution of coadaptation in a two-locus epistatic system. Evolution. 2005;59(11):2324–32.View ArticlePubMedGoogle Scholar
- Jensen L. Random selective advantages of genes and their probabilities of fixation. Genet Res. 1973;21(3):215–9.View ArticlePubMedGoogle Scholar
- Karlin S, Levikson B. Temporal fluctuations in selection intensities: case of small population size. Theor Popul Biol. 1974;6(3):383–412.View ArticleGoogle Scholar
- Takahata N, Ishii K, Matsuda H. Effect of temporal fluctuation of selection coefficient on gene frequency in a population. Proc Natl Acad Sci U S A. 1975;72(11):4541–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Huillet T. On the Karlin–Kimura approaches to the Wright–fisher diffusion with fluctuating selection. J Stat Mech. 2011;2011(2):P02016.View ArticleGoogle Scholar
- Waxman D. A unified treatment of the probability of fixation when population size and the strength of selection change over time. Genetics. 2011;188(4):907–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Uecker H, Hermisson J. On the fixation process of a beneficial mutation in a variable environment. Genetics. 2011;188(4):915–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Peischl S, Kirkpatrick M. Establishment of new mutations in changing environments. Genetics. 2012;191(3):895–906.View ArticlePubMedPubMed CentralGoogle Scholar
- Harris TE. The theory of branching processes. Berlin: Springer-Verlag; 1963.View ArticleGoogle Scholar
- Hartfield M, Otto SP. Recombination and hitchhiking of deleterious alleles. Evolution. 2011;65(9):2421–34.View ArticlePubMedGoogle Scholar
- Lloyd-Smith JO, Schreiber SJ, Kopp PE, Getz WM. Superspreading and the effect of individual variation on disease emergence. Nature. 2005;438(7066):355–9.View ArticlePubMedGoogle Scholar
- Bozic I, Antal T, Ohtsuki H, Carter H, Kim D, Chen S, Karchin R, Kinzler KW, Vogelstein B, Nowak MA. Accumulation of driver and passenger mutations during tumor progression. Proc Natl Acad Sci U S A. 2010;107(43):18545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Bell G, Gonzalez A. Evolutionary rescue can prevent extinction following environmental change. Ecol Lett. 2009;12(9):942–8.View ArticlePubMedGoogle Scholar
- Orr HA, Unckless RL. Population extinction and the genetics of adaptation. Am Nat. 2008;172(2):160–9.View ArticlePubMedGoogle Scholar