An individual-based model for the Lenski experiment, and the deceleration of the relative fitness

The ancestral selection graph for a Λ-asymmetric Moran model

Motivated by the question of the impact of selective advantage in populations with skewed reproduction mechanisms, we study a Moran model with selection. We assume that there are two types of individuals, where the reproductive success of one type is larger than the other. The higher reproductive success may stem from either more frequent reproduction, or from larger numbers of offspring, and is encoded in a measure Λ for each of the two types. Λ-reproduction here means that a whole fraction of the population is replaced at a reproductive event. Our approach consists of constructing a Λ-asymmetric Moran model in which individuals of the two populations compete, rather than considering a Moran model for each population. Provided the measure are ordered stochastically, we can couple them. This allows us to construct the central object of this paper, the Λ-asymmetric ancestral selection graph, leading to a pathwise duality of the forward in time Λ-asymmetric Moran model with its ancestral process. We apply the ancestral selection graph in order to obtain scaling limits of the forward and backward processes, and note that the frequency process converges to the solution of an SDE with discontinuous paths. Finally, we derive a Griffiths representation for the generator of the SDE and use it to find a semi-explicit formula for the probability of fixation of the less beneficial of the two types.

'Evolutionary poker': an agent-based model of interactome emergence and epistasis tested against Lenski's long-term E. coli experiments

A simple agent-based model is presented that produces results matching the experimental data found by Lenski's group for ≤50,000 generations of Escherichia coli bacteria under continuous selective pressure. Although various mathematical models have been devised previously to model the Lenski data, the present model has advantages in terms of overall simplicity and conceptual accessibility. The model also clearly illustrates a number of features of the evolutionary process that are otherwise not obvious, such as the roles of epistasis and historical contingency in adaptation and why evolution is time irreversible ('Dollo's law'). The reason for this irreversibility is that genomes become increasingly integrated or organized, and this organization becomes a novel selective factor itself, against which future generations must compete. Selection for integrated or synergistic networks, systems or sets of mutations or traits, not for individual mutations, confers the main adaptive advantage. The result is a punctuated form of evolution that follows a logarithmic occurrence probability, in which evolution proceeds very quickly when interactomes begin to form but which slows as interactomes become more robust and the difficulty of integrating new mutations increases. Sufficient parameters exist in the game to suggest not only how equilibrium or stasis is reached but also the conditions in which it will be punctuated, the factors governing the rate at which genomic organization occurs and novel traits appear, and how population size, genome size and gene variability affect these.

Segregational instability of multicopy plasmids: A population genetics approach

Plasmids are extra-chromosomal genetic elements that encode a wide variety of phenotypes and can be maintained in bacterial populations through vertical and horizontal transmission, thus increasing bacterial adaptation to hostile environmental conditions like those imposed by antimicrobial substances. To circumvent the segregational instability resulting from randomly distributing plasmids between daughter cells upon division, nontransmissible plasmids tend to be carried in multiple copies per cell, with the added benefit of exhibiting increased gene dosage and resistance levels. But carrying multiple copies also results in a high metabolic burden to the bacterial host, therefore reducing the overall fitness of the population. This trade-off poses an existential question for plasmids: What is the optimal plasmid copy number? In this manuscript, we address this question by postulating and analyzing a population genetics model to evaluate the interaction between selective pressure, the number of plasmid copies carried by each cell, and the metabolic burden associated with plasmid bearing in the absence of selection for plasmid-encoded traits. Parameter values of the model were estimated experimentally using Escherichia coli K12 carrying a multicopy plasmid encoding for a fluorescent protein and bla TEM-1, a gene conferring resistance to β-lactam antibiotics. By numerically determining the optimal plasmid copy number for constant and fluctuating selection regimes, we show that plasmid copy number is a highly optimized evolutionary trait that depends on the rate of environmental fluctuation and balances the benefit between increased stability in the absence of selection with the burden associated with carrying multiple copies of the plasmid.

Kingman’s model with random mutation probabilities: convergence and condensation I

For a one-locus haploid infinite population with discrete generations, the celebrated model of Kingman describes the evolution of fitness distributions under the competition of selection and mutation, with a constant mutation probability. This paper generalises Kingman’s model by using independent and identically distributed random mutation probabilities, to reflect the influence of a random environment. The weak convergence of fitness distributions to the globally stable equilibrium is proved. Condensation occurs when almost surely a positive proportion of the population travels to and condenses at the largest fitness value. Condensation may occur when selection is favoured over mutation. A criterion for the occurrence of condensation is given.

Haldane's formula in Cannings models: the case of moderately strong selection

For a class of Cannings models we prove Haldane’s formula, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi (s_N) \sim \frac{2s_N}{\rho ^2}$$\end{document}π(sN)∼2sNρ2, for the fixation probability of a single beneficial mutant in the limit of large population size N and in the regime of moderately strong selection, i.e. for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s_N \sim N^{-b}$$\end{document}sN∼N-b and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0< b<1/2$$\end{document}0<b<1/2. Here, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s_N$$\end{document}sN is the selective advantage of an individual carrying the beneficial type, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho ^2$$\end{document}ρ2 is the (asymptotic) offspring variance. Our assumptions on the reproduction mechanism allow for a coupling of the beneficial allele’s frequency process with slightly supercritical Galton–Watson processes in the early phase of fixation.

The Wright-Fisher model with efficiency

In this article, we propose a Wright-Fisher model with two types of individuals: the inefficient individuals, those who need more resources to reproduce and can have a higher growth rate, and the efficient individuals. In this model, the total amount of resource N is fixed, and the population size varies randomly depending on the number of efficient individuals. We show that, as N increases, the frequency process of efficient individuals converges to a diffusion which is a generalization of the Wright-Fisher diffusion with selection. The genealogy of this model is given by a branching-coalescing process that we call the Ancestral Selection/Efficiency Graph, and that is an extension of the Ancestral Selection Graph (Krone and Neuhauser, 1997a,b). The main contribution of this paper is that, in evolving populations, inefficiency can arise as a promoter of selective advantage and not necessarily as a trade-off.