AMS eBook CollectionsOne of the world's most respected mathematical collections, available in digital format for your library or institution
Critical population and error threshold on the sharp peak landscape for a Moran model
About this Title
Raphaël Cerf, Université Paris Sud, Mathématique, Bâtiment $425$, 91405 Orsay Cedex–France
Publication: Memoirs of the American Mathematical Society
Publication Year:
2015; Volume 233, Number 1096
ISBNs: 978-1-4704-0967-8 (print); 978-1-4704-1964-6 (online)
DOI: https://doi.org/10.1090/memo/1096
Published electronically: May 19, 2014
Keywords: Moran,
quasispecies,
error threshold
MSC: Primary 60F10; Secondary 92D25
View full volume PDF
Table of Contents
Chapters
- 1. Introduction
- 2. The Model
- 3. Main Results
- 4. Coupling
- 5. Normalized Model
- 6. Lumping
- 7. Monotonicity
- 8. Stochastic Bounds
- 9. Birth and Death Processes
- 10. The Neutral Phase
Abstract
The goal of this work is to propose a finite population counterpart to Eigen’s model, which incorporates stochastic effects. We consider a Moran model describing the evolution of a population of size $m$ of chromosomes of length $\ell$ over an alphabet of cardinality $\kappa$. The mutation probability per locus is $q$. We deal only with the sharp peak landscape: the replication rate is $\sigma >1$ for the master sequence and $1$ for the other sequences. We study the equilibrium distribution of the process in the regime where \[
\displaylines { \ell \to +\infty \,,\qquad m\to +\infty \,,\qquad q\to 0\,,\cr {\ell q} \to a\in ]0,+\infty [\,, \qquad \frac {m}{\ell }\to \alpha \in [0,+\infty ]\,.}\] We obtain an equation $\alpha \,\phi (a)=\ln \kappa$ in the parameter space $(a,\alpha )$ separating the regime where the equilibrium population is totally random from the regime where a quasispecies is formed. We observe the existence of a critical population size necessary for a quasispecies to emerge and we recover the finite population counterpart of the error threshold. Moreover, in the limit of very small mutations, we obtain a lower bound on the population size allowing the emergence of a quasispecies: if $\alpha < \ln \kappa /\ln \sigma$ then the equilibrium population is totally random, and a quasispecies can be formed only when $\alpha \geq \ln \kappa /\ln \sigma$. Finally, in the limit of very large populations, we recover an error catastrophe reminiscent of Eigen’s model: if $\sigma \exp (-a)\leq 1$ then the equilibrium population is totally random, and a quasispecies can be formed only when $\sigma \exp (-a)>1$. These results are supported by computer simulations.