Source author record

Ellen Baake

Ellen Baake appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2019arXiv

Modelling and simulating Lenski's long-term evolution experiment

We revisit the model by Wiser, Ribeck, and Lenski (Science \textbf{342} (2013), 1364--1367), which describes how the mean fitness increases over time due to beneficial mutations in Lenski's long-term evolution experiment. We develop the model further both conceptually and mathematically. Conceptually, we describe the experiment with the help of a Cannings model with mutation and selection, where the latter includes diminishing returns epistasis. The analysis sheds light on the growth dynamics within every single day and reveals a runtime effect, that is, the shortening of the daily growth period with increasing fitness; and it allows to clarify the contribution of epistasis to the mean fitness curve. Mathematically, we explain rigorous results in terms of a law of large numbers (in the limit of infinite population size and for a certain asymptotic parameter regime), and present approximations based on heuristics and supported by simulations for finite populations.

preprint2018arXiv

Fragmentation process, pruning poset for rooted forests, and Möbius inversion

We consider a discrete-time Markov chain, called fragmentation process, that describes a specific way of successively removing objects from a linear arrangement. The process arises in population genetics and describes the ancestry of the genetic material of individuals in a population experiencing recombination. We aim at the law of the process over time. To this end, we investigate sets of realisations of this process that agree with respect to a specific order of events and represent each such set by a rooted (binary) tree. The probability of each tree is, in turn, obtained by Möbius inversion on a suitable poset of all rooted forests that can be obtained from the tree by edge deletion; we call this poset the \textit{pruning poset}. Dependencies within the fragments make it difficult to obtain explicit expressions for the probabilities of the trees. We therefore construct an auxiliary process for every given tree, which is i.i.d. over time, and which allows to give a pathwise construction of realisations that match the tree.

preprint2017arXiv

Lines of descent under selection

We review recent progress on ancestral processes related to mutation-selection models, both in the deterministic and the stochastic setting. We mainly rely on two concepts, namely, the killed ancestral selection graph and the pruned lookdown ancestral selection graph. The killed ancestral selection graph gives a representation of the type of a random individual from a stationary population, based upon the individual's potential ancestry back until the mutations that define the individual's type. The pruned lookdown ancestral selection graph allows one to trace the ancestry of individuals from a stationary distribution back into the distant past, thus leading to the stationary distribution of ancestral types. We illustrate the results by applying them to a prototype model for the error threshold phenomenon.

preprint2016arXiv

Haldane linearisation done right: Solving the nonlinear recombination equation the easy way

The nonlinear recombination equation from population genetics has a long history and is notoriously difficult to solve, both in continuous and in discrete time. This is particularly so if one aims at full generality, thus also including degenerate parameter cases. Due to recent progress for the continuous time case via the identification of an underlying stochastic fragmentation process, it became clear that a direct general solution at the level of the corresponding ODE itself should also be possible. This paper shows how to do it, and how to extend the approach to the discrete-time case as well.

preprint2016arXiv

The common ancestor type distribution of a $Λ$-Wright-Fisher process with selection and mutation

Using graphical methods based on a `lookdown' and pruned version of the {\em ancestral selection graph}, we obtain a representation of the type distribution of the ancestor in a two-type Wright-Fisher population with mutation and selection, conditional on the overall type frequency in the old population. This extends results from Lenz, Kluth, Baake, and Wakolbinger (Theor. Pop. Biol., 103 (2015), 27-37) to the case of heavy-tailed offspring, directed by a reproduction measure $Λ$. The representation is in terms of the equilibrium tail probabilities of the line-counting process $L$ of the graph. We identify a strong pathwise Siegmund dual of $L$, and characterise the equilibrium tail probabilities of $L$ in terms of hitting probabilities of the dual process.

preprint2015arXiv

Feller's Contributions to Mathematical Biology

This is a review of William Feller's important contributions to mathematical biology. The seminal paper [Feller1951] "Diffusion processes in genetics" was particularly influential on the development of stochastic processes at the interface to evolutionary biology, and interesting ideas in this direction (including a first characterization of what is nowadays known as "Feller's branching diffusion") already shaped up in the paper [Feller 1939] (written in German) "The foundations of a probabistic treatment of Volterra's theory of the struggle for life". Feller's article "On fitness and the cost of natural selection" [Feller 1967] contains a critical analysis of the concept of "genetic load".

preprint2015arXiv

Looking down in the ancestral selection graph: A probabilistic approach to the common ancestor type distribution

In a (two-type) Wright-Fisher diffusion with directional selection and two-way mutation, let $x$ denote today's frequency of the beneficial type, and given $x$, let $h(x)$ be the probability that, among all individuals of today's population, the individual whose progeny will eventually take over in the population is of the beneficial type. Fearnhead [Fearnhead, P., 2002. The common ancestor at a nonneutral locus. J. Appl. Probab. 39, 38-54] and Taylor [Taylor, J. E., 2007. The common ancestor process for a Wright-Fisher diffusion. Electron. J. Probab. 12, 808-847] obtained a series representation for $h(x)$. We develop a construction that contains elements of both the ancestral selection graph and the lookdown construction and includes pruning of certain lines upon mutation. Besides being interesting in its own right, this construction allows a transparent derivation of the series coefficients of $h(x)$ and gives them a probabilistic meaning.

preprint2015arXiv

The general recombination equation in continuous time and its solution

The process of recombination in population genetics, in its deterministic limit, leads to a nonlinear ODE in the Banach space of finite measures on a locally compact product space. It has an embedding into a larger family of nonlinear ODEs that permits a systematic analysis with lattice-theoretic methods for general partitions of finite sets. We discuss this type of system, reduce it to an equivalent finite-dimensional nonlinear problem, and establish a connection with an ancestral partitioning process, backward in time. We solve the finite-dimensional problem recursively for generic sets of parameters and briefly discuss the singular cases, and how to extend the solution to this situation.

preprint2014arXiv

Rare event simulation in immune biology: Models of negative selection in T-cell maturation

We present a probabilistic T-cell model that includes negative selection and takes contrasting models of tissue-restricted antigen (TRA) expression in the thymus into account. We start from the basic model of van den Berg, Rand, and Burroughs (2001) and include negative selection via individual-based T-cell modelling, in which each T-cell is defined by its stimulation rates to all relevant self antigens. We present a simulation approach based on partial tilting of the stimulation rates recognized by a single T-cell. We investigate the effects of negative selection for diverging modes of thymic antigen presentation, namely arbitrary TRA presentation, and more or less strict emulation of tissue-specific cell lines. We observe that negative selection leads to truncation of the tail of the distribution of the stimulation rates mature T-cells receive from self antigens, i.e., the self background is reduced. This increases the activation probabilities of single T-cells in the presence of non-self antigens.

preprint2013arXiv

Single--crossover recombination and ancestral recombination trees

We consider the Wright-Fisher model for a population of $N$ individuals, each identified with a sequence of a finite number of sites, and single-crossover recombination between them. We trace back the ancestry of single individuals from the present population. In the $N \to \infty$ limit without rescaling of parameters or time, this ancestral process is described by a random tree, whose branching events correspond to the splitting of the sequence due to recombination. With the help of a decomposition of the trees into subtrees, we calculate the probabilities of the topologies of the ancestral trees. At the same time, these probabilities lead to a semi-explicit solution of the deterministic single-crossover equation. The latter is a discrete-time dynamical system that emerges from the Wright-Fisher model via a law of large numbers and has been waiting for a solution for many decades.

preprint2013arXiv

The common ancestor process revisited

We consider the Moran model in continuous time with two types, mutation, and selection. We concentrate on the ancestral line and its stationary type distribution. Building on work by Fearnhead (J. Appl. Prob. 39 (2002), 38-54) and Taylor (Electron. J. Probab. 12 (2007), 808-847), we characterise this distribution via the fixation probability of the offspring of all individuals of favourable type (regardless of the offspring's types). We concentrate on a finite population and stay with the resulting discrete setting all the way through. This way, we extend previous results and gain new insight into the underlying particle picture.

preprint2013arXiv

The Moran model with selection: Fixation probabilities, ancestral lines, and an alternative particle representation

We reconsider the Moran model in continuous time with population size $N$, two allelic types, and selection. We introduce a new particle representation, which we call the labelled Moran model, and which has the same distribution of type frequencies as the original Moran model, provided the initial values are chosen appropriately. In the new model, individuals are labelled $1,2, \dots, N$; neutral resampling events may take place between arbitrary labels, whereas selective events only occur in the direction of increasing labels. With the help of elementary methods only, we not only recover fixation probabilities, but also obtain detailed insight into the number and nature of the selective events that play a role in the fixation process forward in time.

preprint2011arXiv

Deterministic and stochastic aspects of single-crossover recombination

This contribution is concerned with mathematical models for the dynamics of the genetic composition of populations evolving under recombination. Recombination is the genetic mechanism by which two parent individuals create the mixed type of their offspring during sexual reproduction. The corresponding models are large, nonlinear dynamical systems (for the deterministic treatment that applies in the infinite-population limit), or interacting particle systems (for the stochastic treatment required for finite populations). We review recent progress on these difficult problems. In particular, we present a closed solution of the deterministic continuous-time system, for the important special case of single crossovers; we extract an underlying linearity; we analyse how this carries over to the corresponding stochastic setting; and we provide a solution of the analogous deterministic discrete-time dynamics, in terms of its generalised eigenvalues and a simple recursion for the corresponding coefficients.

preprint2011arXiv

Moment closure in a Moran model with recombination

We extend the Moran model with single-crossover recombination to include general recombination and mutation. We show that, in the case without resampling, the expectations of products of marginal processes defined via partitions of sites form a closed hierarchy, which is exhaustively described by a finite system of differential equations. One thus has the exceptional situation of moment closure in a nonlinear system. Surprisingly, this property is lost when resampling (i.e., genetic drift) is included.

preprint2009arXiv

Single--crossover recombination in discrete time

Modelling the process of recombination leads to a large coupled nonlinear dynamical system. Here, we consider a particular case of recombination in {\em discrete} time, allowing only for {\em single crossovers}. While the analogous dynamics in {\em continuous} time admits a closed solution, this no longer works for discrete time. A more general model (i.e. without the restriction to single crossovers) has been studied before and was solved algorithmically by means of Haldane linearisation. Using the special formalism introduced by Baake and Baake (2003), we obtain further insight into the single-crossover dynamics and the particular difficulties that arise in discrete time. We then transform the equations to a solvable system in a two-step procedure: linearisation followed by diagonalisation. Still, the coefficients of the second step must be determined in a recursive manner, but once this is done for a given system, they allow for an explicit solution valid for all times.

preprint1999arXiv

Biological evolution through mutation, selection, and drift: An introductory review

Motivated by present activities in (statistical) physics directed towards biological evolution, we review the interplay of three evolutionary forces: mutation, selection, and genetic drift. The review addresses itself to physicists and intends to bridge the gap between the biological and the physical literature. We first clarify the terminology and recapitulate the basic models of population genetics, which describe the evolution of the composition of a population under the joint action of the various evolutionary forces. Building on these foundations, we specify the ingredients explicitly, namely, the various mutation models and fitness landscapes. We then review recent developments concerning models of mutational degradation. These predict upper limits for the mutation rate above which mutation can no longer be controlled by selection, the most important phenomena being error thresholds, Muller's ratchet, and mutational meltdowns. Error thresholds are deterministic phenomena, whereas Muller's ratchet requires the stochastic component brought about by finite population size. Mutational meltdowns additionally rely on an explicit model of population dynamics, and describe the extinction of populations. Special emphasis is put on the mutual relationship between these phenomena. Finally, a few connections with the process of molecular evolution are established.