1
|
Igelbrink JL, González Casanova A, Smadi C, Wakolbinger A. Muller's ratchet in a near-critical regime: Tournament versus fitness proportional selection. Theor Popul Biol 2024; 158:121-138. [PMID: 38844263 DOI: 10.1016/j.tpb.2024.06.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 01/31/2024] [Accepted: 06/01/2024] [Indexed: 06/22/2024]
Abstract
Muller's ratchet, in its prototype version, models a haploid, asexual population whose size N is constant over the generations. Slightly deleterious mutations are acquired along the lineages at a constant rate, and individuals carrying less mutations have a selective advantage. The classical variant considers fitness proportional selection, but other fitness schemes are conceivable as well. Inspired by the work of Etheridge et al. (2009) we propose a parameter scaling which fits well to the "near-critical" regime that was in the focus of Etheridge et al. (2009) (and in which the mutation-selection ratio diverges logarithmically as N→∞). Using a Moran model, we investigate the"rule of thumb" given in Etheridge et al. (2009) for the click rate of the "classical ratchet" by putting it into the context of new results on the long-time evolution of the size of the best class of the ratchet with (binary) tournament selection. This variant of Muller's ratchet was introduced in González Casanova et al. (2023), and was analysed there in a subcritical parameter regime. Other than that of the classical ratchet, the size of the best class of the tournament ratchet follows an autonomous dynamics up to the time of its extinction. It turns out that, under a suitable correspondence of the model parameters, this dynamics coincides with the so called Poisson profile approximation of the dynamics of the best class of the classical ratchet.
Collapse
Affiliation(s)
- J L Igelbrink
- Goethe-Universität, Institut für Mathematik, Frankfurt am Main, 60629, Germany; Johannes Gutenberg-Universität, Institut für Mathematik, Staudingerweg 9, Mainz, 55128, Germany.
| | - A González Casanova
- Universidad Nacional Autónoma de México (UNAM), Instituto de Matemáticas, Circuito exterior, Ciudad Universitaria, 04510, Mexico; Department of Statistics, 367 Evans Hall, University of California, Berkeley, CA 94720-386, United States of America
| | - C Smadi
- Univ. Grenoble Alpes, INRAE, LESSEM, Grenoble, 38000, France; Univ. Grenoble Alpes, CNRS, Institut Fourier, Grenoble, 38000, France
| | - A Wakolbinger
- Goethe-Universität, Institut für Mathematik, Frankfurt am Main, 60629, Germany
| |
Collapse
|
2
|
Olofsson P, Chipkin L, Daileda RC, Azevedo RBR. Mutational meltdown in asexual populations doomed to extinction. J Math Biol 2023; 87:88. [PMID: 37994999 DOI: 10.1007/s00285-023-02019-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2023] [Revised: 07/03/2023] [Accepted: 10/20/2023] [Indexed: 11/24/2023]
Abstract
Asexual populations are expected to accumulate deleterious mutations through a process known as Muller's ratchet. Lynch and colleagues proposed that the ratchet eventually results in a vicious cycle of mutation accumulation and population decline that drives populations to extinction. They called this phenomenon mutational meltdown. Here, we analyze mutational meltdown using a multi-type branching process model where, in the presence of mutation, populations are doomed to extinction. We analyse the change in size and composition of the population and the time of extinction under this model.
Collapse
Affiliation(s)
- Peter Olofsson
- Department of Mathematics, Trinity University, San Antonio, TX, 78212, USA
- Department of Mathematics, Physics and Chemical Engineering, Jönköping University, 551 11, Jönköping, Sweden
| | - Logan Chipkin
- Department of Biology and Biochemistry, University of Houston, Houston, TX, 77204, USA
| | - Ryan C Daileda
- Department of Mathematics, Trinity University, San Antonio, TX, 78212, USA
| | - Ricardo B R Azevedo
- Department of Biology and Biochemistry, University of Houston, Houston, TX, 77204, USA.
| |
Collapse
|
3
|
Roitershtein A, Rastegar R, Chapkin RS, Ivanov I. Extinction scenarios in evolutionary processes: a multinomial Wright-Fisher approach. J Math Biol 2023; 87:63. [PMID: 37751048 PMCID: PMC10586398 DOI: 10.1007/s00285-023-01993-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2019] [Revised: 08/16/2023] [Accepted: 08/31/2023] [Indexed: 09/27/2023]
Abstract
We study a discrete-time multi-type Wright-Fisher population process. The mean-field dynamics of the stochastic process is induced by a general replicator difference equation. We prove several results regarding the asymptotic behavior of the model, focusing on the impact of the mean-field dynamics on it. One of the results is a limit theorem that describes sufficient conditions for an almost certain path to extinction, first eliminating the type which is the least fit at the mean-field equilibrium. The effect is explained by the metastability of the stochastic system, which under the conditions of the theorem spends almost all time before the extinction event in a neighborhood of the equilibrium. In addition to the limit theorems, we propose a maximization principle for a general deterministic replicator dynamics and study its implications for the stochastic model.
Collapse
Affiliation(s)
| | - Reza Rastegar
- Occidental Petroleum Corporation, Houston, TX, 77046, USA
| | - Robert S Chapkin
- Department of Nutrition - Program in Integrative Nutrition & Complex Diseases, Texas A &M University, College Station, TX, 77843, USA
| | - Ivan Ivanov
- Department of Veterinary Physiology and Pharmacology, Texas A &M University, College Station, TX, 77843, USA.
| |
Collapse
|
4
|
Bräutigam C, Smerlak M. Diffusion approximations in population genetics and the rate of Muller's ratchet. J Theor Biol 2022; 550:111236. [PMID: 35926567 DOI: 10.1016/j.jtbi.2022.111236] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2021] [Revised: 07/13/2022] [Accepted: 07/25/2022] [Indexed: 10/16/2022]
Abstract
The Wright-Fisher binomial model of allele frequency change is often approximated by a scaling limit in which selection, mutation and drift all decrease at the same 1/N rate. This construction restricts the applicability of the resulting 'Wright-Fisher diffusion equation' to the weak selection, weak mutation regime of evolution. We argue that diffusion approximations of the Wright-Fisher model can be used more generally, for instance in cases where genetic drift is much weaker than selection. One important example of this regime is Muller's ratchet phenomenon, whereby deleterious mutations slowly but irreversibly accumulate through rare stochastic fluctuations. Using a modified diffusion equation we derive improved analytical estimates for the mean click time of the ratchet.
Collapse
Affiliation(s)
- Camila Bräutigam
- Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany
| | - Matteo Smerlak
- Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany.
| |
Collapse
|
5
|
Sakamoto T, Innan H. Muller's ratchet of the Y chromosome with gene conversion. Genetics 2022; 220:iyab204. [PMID: 34791206 PMCID: PMC8733426 DOI: 10.1093/genetics/iyab204] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Accepted: 10/28/2021] [Indexed: 11/13/2022] Open
Abstract
Muller's ratchet is a process in which deleterious mutations are fixed irreversibly in the absence of recombination. The degeneration of the Y chromosome, and the gradual loss of its genes, can be explained by Muller's ratchet. However, most theories consider single-copy genes, and may not be applicable to Y chromosomes, which have a number of duplicated genes in many species, which are probably undergoing concerted evolution by gene conversion. We developed a model of Muller's ratchet to explore the evolution of the Y chromosome. The model assumes a nonrecombining chromosome with both single-copy and duplicated genes. We used analytical and simulation approaches to obtain the rate of gene loss in this model, with special attention to the role of gene conversion. Homogenization by gene conversion makes both duplicated copies either mutated or intact. The former promotes the ratchet, and the latter retards, and we ask which of these counteracting forces dominates under which conditions. We found that the effect of gene conversion is complex, and depends upon the fitness effect of gene duplication. When duplication has no effect on fitness, gene conversion accelerates the ratchet of both single-copy and duplicated genes. If duplication has an additive fitness effect, the ratchet of single-copy genes is accelerated by gene duplication, regardless of the gene conversion rate, whereas gene conversion slows the degeneration of duplicated genes. Our results suggest that the evolution of the Y chromosome involves several parameters, including the fitness effect of gene duplication by increasing dosage and gene conversion rate.
Collapse
Affiliation(s)
- Takahiro Sakamoto
- Department of Evolutionary Studies of Biosystems, SOKENDAI, The Graduate University for Advanced Studies, Hayama, Kanagawa 240-0193, Japan
| | - Hideki Innan
- Department of Evolutionary Studies of Biosystems, SOKENDAI, The Graduate University for Advanced Studies, Hayama, Kanagawa 240-0193, Japan
| |
Collapse
|
6
|
Zhao L, Lascoux M, Waxman D. An informational transition in conditioned Markov chains: Applied to genetics and evolution. J Theor Biol 2016; 402:158-70. [PMID: 27105672 DOI: 10.1016/j.jtbi.2016.04.021] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2015] [Revised: 02/19/2016] [Accepted: 04/17/2016] [Indexed: 11/18/2022]
Abstract
In this work we assume that we have some knowledge about the state of a population at two known times, when the dynamics is governed by a Markov chain such as a Wright-Fisher model. Such knowledge could be obtained, for example, from observations made on ancient and contemporary DNA, or during laboratory experiments involving long term evolution. A natural assumption is that the behaviour of the population, between observations, is related to (or constrained by) what was actually observed. The present work shows that this assumption has limited validity. When the time interval between observations is larger than a characteristic value, which is a property of the population under consideration, there is a range of intermediate times where the behaviour of the population has reduced or no dependence on what was observed and an equilibrium-like distribution applies. Thus, for example, if the frequency of an allele is observed at two different times, then for a large enough time interval between observations, the population has reduced or no dependence on the two observed frequencies for a range of intermediate times. Given observations of a population at two times, we provide a general theoretical analysis of the behaviour of the population at all intermediate times, and determine an expression for the characteristic time interval, beyond which the observations do not constrain the population's behaviour over a range of intermediate times. The findings of this work relate to what can be meaningfully inferred about a population at intermediate times, given knowledge of terminal states.
Collapse
Affiliation(s)
- Lei Zhao
- Centre for Computational Systems Biology, Fudan University, 220 Handan Road, Shanghai 200433, PR China
| | - Martin Lascoux
- Centre for Computational Systems Biology, Fudan University, 220 Handan Road, Shanghai 200433, PR China; Evolutionary Biology Center, Department of Ecology and Genetics, Uppsala University, Uppsala 75236, Sweden
| | - David Waxman
- Centre for Computational Systems Biology, Fudan University, 220 Handan Road, Shanghai 200433, PR China.
| |
Collapse
|
7
|
The influence of genetic drift on the formation and stability of polymorphisms arising from negative frequency-dependent selection. J Theor Biol 2016; 391:51-64. [DOI: 10.1016/j.jtbi.2015.11.011] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2015] [Revised: 11/13/2015] [Accepted: 11/17/2015] [Indexed: 11/20/2022]
|
8
|
Effect of drift, selection and recombination on the equilibrium frequency of deleterious mutations. J Theor Biol 2015; 365:238-46. [PMID: 25451760 DOI: 10.1016/j.jtbi.2014.10.023] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2014] [Revised: 08/30/2014] [Accepted: 10/22/2014] [Indexed: 11/22/2022]
Abstract
We study the stationary state of a population evolving under the action of random genetic drift, selection and recombination in which both deleterious and reverse beneficial mutations can occur. We find that the equilibrium fraction of deleterious mutations decreases as the population size is increased. We calculate exactly the steady state frequency in a nonrecombining population when population size is infinite and for a neutral finite population, and obtain bounds on the fraction of deleterious mutations. We also find that for small and very large populations, the number of deleterious mutations depends weakly on recombination, but for moderately large populations, recombination alleviates the effect of deleterious mutations. An analytical argument shows that recombination decreases disadvantageous mutations appreciably when beneficial mutations are rare as is the case in adapting microbial populations, whereas it has a moderate effect on codon bias where the mutation rates between the preferred and unpreferred codons are comparable.
Collapse
|
9
|
Abstract
Competition between independently arising beneficial mutations is enhanced in spatial populations due to the linear rather than exponential growth of clones. Recent theoretical studies have pointed out that the resulting fitness dynamics is analogous to a surface growth process, where new layers nucleate and spread stochastically, leading to the build up of scale-invariant roughness. This scenario differs qualitatively from the standard view of adaptation in that the speed of adaptation becomes independent of population size while the fitness variance does not. Here we exploit recent progress in the understanding of surface growth processes to obtain precise predictions for the universal, non-Gaussian shape of the fitness distribution for one-dimensional habitats, which are verified by simulations. When the mutations are deleterious rather than beneficial the problem becomes a spatial version of Muller's ratchet. In contrast to the case of well-mixed populations, the rate of fitness decline remains finite even in the limit of an infinite habitat, provided the ratio [Formula: see text] between the deleterious mutation rate and the square of the (negative) selection coefficient is sufficiently large. Using, again, an analogy to surface growth models we show that the transition between the stationary and the moving state of the ratchet is governed by directed percolation.
Collapse
Affiliation(s)
- Jakub Otwinowski
- Emory University, Physics Department Atlanta, Georgia, USA. University of Pennsylvania, Biology Department, Philadelphia, Pennsylvania, USA
| | | |
Collapse
|
10
|
Metzger JJ, Eule S. Distribution of the fittest individuals and the rate of Muller's ratchet in a model with overlapping generations. PLoS Comput Biol 2013; 9:e1003303. [PMID: 24244123 PMCID: PMC3820511 DOI: 10.1371/journal.pcbi.1003303] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2013] [Accepted: 09/10/2013] [Indexed: 11/17/2022] Open
Abstract
Muller's ratchet is a paradigmatic model for the accumulation of deleterious mutations in a population of finite size. A click of the ratchet occurs when all individuals with the least number of deleterious mutations are lost irreversibly due to a stochastic fluctuation. In spite of the simplicity of the model, a quantitative understanding of the process remains an open challenge. In contrast to previous works, we here study a Moran model of the ratchet with overlapping generations. Employing an approximation which describes the fittest individuals as one class and the rest as a second class, we obtain closed analytical expressions of the ratchet rate in the rare clicking regime. As a click in this regime is caused by a rare, large fluctuation from a metastable state, we do not resort to a diffusion approximation but apply an approximation scheme which is especially well suited to describe extinction events from metastable states. This method also allows for a derivation of expressions for the quasi-stationary distribution of the fittest class. Additionally, we confirm numerically that the formulation with overlapping generations leads to the same results as the diffusion approximation and the corresponding Wright-Fisher model with non-overlapping generations. Muller's ratchet is a paradigmatic model in population genetics which describes the fixation of a deleterious mutation in a population of finite size due to an unfortunate stochastic fluctuation. Obtaining quantitative predictions of the ratchet rate, i.e. the frequency with which such a mutation fixes, is believed to be important for understanding a broad range of effects ranging from the degeneration of the Y-chromosome to the evolution of sex as a means of avoiding the fixation of deleterious mutations. To obtain a better understanding of how Muller's ratchet operates, we have considered a model with overlapping generations, which allows for the application of methods specifically tailored for the analysis of rare stochastic fluctuations which drive the ratchet. We obtain concise and accurate results for the rate of Muller's ratchet. Additionally, we are able to predict the full distribution of the frequency of the fittest individuals, a quantity of central interest in understanding the ratchet rate and possibly experimentally much more accessible than the rate, in particular when the ratchet rate is very large.
Collapse
Affiliation(s)
- Jakob J Metzger
- Max Planck Institute for Dynamics and Self-Organization (MPIDS), Göttingen, Germany ; Institute for Nonlinear Dynamics, Department of Physics, University of Göttingen, Göttingen, Germany
| | | |
Collapse
|
11
|
Abstract
Adaptive landscape, proposed by Sewall Wright, has provided a conceptual framework to describe dynamical behaviours. However, it is still a challenge to explicitly construct such a landscape, and apply it to quantify interesting evolutionary processes. This is particularly true for neutral evolution. In this work, the authors study one-dimensional Wright Fisher process, and analytically obtain an adaptive landscape as a potential function. They provide the complete characterisation for dynamical behaviours of all possible mutation rates under the influence of mutation and random drift. This same analysis has been applied to situations with additive selection and random drift for all possible selection rates. The critical state dividing the basins of two stable states is directly obtained by the landscape. In addition, the landscape is able to handle situations with pure random drift, which would be non-normalisable for its stationary distribution. The nature of non-normalisation is from the singularity of adaptive landscape. In addition, they propose a new type of neutral evolution. It has the same probability for all possible states. The new type of neutral evolution describes the non-neutral alleles with 0%. They take the equal effect of mutation and random drift as an example.
Collapse
|
12
|
Good BH, Desai MM. Fluctuations in fitness distributions and the effects of weak linked selection on sequence evolution. Theor Popul Biol 2013; 85:86-102. [PMID: 23337315 DOI: 10.1016/j.tpb.2013.01.005] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2012] [Revised: 01/02/2013] [Accepted: 01/11/2013] [Indexed: 02/02/2023]
Abstract
Evolutionary dynamics and patterns of molecular evolution are strongly influenced by selection on linked regions of the genome, but our quantitative understanding of these effects remains incomplete. Recent work has focused on predicting the distribution of fitness within an evolving population, and this forms the basis for several methods that leverage the fitness distribution to predict the patterns of genetic diversity when selection is strong. However, in weakly selected populations random fluctuations due to genetic drift are more severe, and neither the distribution of fitness nor the sequence diversity within the population are well understood. Here, we briefly review the motivations behind the fitness-distribution picture, and summarize the general approaches that have been used to analyze this distribution in the strong-selection regime. We then extend these approaches to the case of weak selection, by outlining a perturbative treatment of selection at a large number of linked sites. This allows us to quantify the stochastic behavior of the fitness distribution and yields exact analytical predictions for the sequence diversity and substitution rate in the limit that selection is weak.
Collapse
Affiliation(s)
- Benjamin H Good
- Department of Organismic and Evolutionary Biology, Department of Physics, and FAS Center for Systems Biology, Harvard University, United States
| | | |
Collapse
|
13
|
Pfaffelhuber P, Staab PR, Wakolbinger A. Muller’s ratchet with compensatory mutations. ANN APPL PROBAB 2012. [DOI: 10.1214/11-aap836] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
14
|
Absorbing phenomena and escaping time for Muller's ratchet in adaptive landscape. BMC SYSTEMS BIOLOGY 2012; 6 Suppl 1:S10. [PMID: 23046686 PMCID: PMC3403021 DOI: 10.1186/1752-0509-6-s1-s10] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
BACKGROUND The accumulation of deleterious mutations of a population directly contributes to the fate as to how long the population would exist, a process often described as Muller's ratchet with the absorbing phenomenon. The key to understand this absorbing phenomenon is to characterize the decaying time of the fittest class of the population. Adaptive landscape introduced by Wright, a re-emerging powerful concept in systems biology, is used as a tool to describe biological processes. To our knowledge, the dynamical behaviors for Muller's ratchet over the full parameter regimes are not studied from the point of the adaptive landscape. And the characterization of the absorbing phenomenon is not yet quantitatively obtained without extraneous assumptions as well. METHODS We describe how Muller's ratchet can be mapped to the classical Wright-Fisher process in both discrete and continuous manners. Furthermore, we construct the adaptive landscape for the system analytically from the general diffusion equation. The constructed adaptive landscape is independent of the existence and normalization of the stationary distribution. We derive the formula of the single click time in finite and infinite potential barrier for all parameters regimes by mean first passage time. RESULTS We describe the dynamical behavior of the population exposed to Muller's ratchet in all parameters regimes by adaptive landscape. The adaptive landscape has rich structures such as finite and infinite potential, real and imaginary fixed points. We give the formula about the single click time with finite and infinite potential. And we find the single click time increases with selection rates and population size increasing, decreases with mutation rates increasing. These results provide a new understanding of infinite potential. We analytically demonstrate the adaptive and unadaptive states for the whole parameters regimes. Interesting issues about the parameters regions with the imaginary fixed points is demonstrated. Most importantly, we find that the absorbing phenomenon is characterized by the adaptive landscape and the single click time without any extraneous assumptions. These results suggest a graphical and quantitative framework to study the absorbing phenomenon.
Collapse
|
15
|
Abstract
The accumulation of deleterious mutations is driven by rare fluctuations that lead to the loss of all mutation free individuals, a process known as Muller’s ratchet. Even though Muller’s ratchet is a paradigmatic process in population genetics, a quantitative understanding of its rate is still lacking. The difficulty lies in the nontrivial nature of fluctuations in the fitness distribution, which control the rate of extinction of the fittest genotype. We address this problem using the simple but classic model of mutation selection balance with deleterious mutations all having the same effect on fitness. We show analytically how fluctuations among the fittest individuals propagate to individuals of lower fitness and have dramatically amplified effects on the bulk of the population at a later time. If a reduction in the size of the fittest class reduces the mean fitness only after a delay, selection opposing this reduction is also delayed. This delayed restoring force speeds up Muller’s ratchet. We show how the delayed response can be accounted for using a path-integral formulation of the stochastic dynamics and provide an expression for the rate of the ratchet that is accurate across a broad range of parameters.
Collapse
|
16
|
A unified treatment of the probability of fixation when population size and the strength of selection change over time. Genetics 2011; 188:907-13. [PMID: 21527780 DOI: 10.1534/genetics.111.129288] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The fixation probability is determined when population size and selection change over time and differs from Kimura's result, with long-term implications for a population. It is found that changes in population size are not equivalent to the corresponding changes in selection and can result in less drift than anticipated.
Collapse
|