Source author record

Anton Wakolbinger

Anton Wakolbinger 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

22works
4topics
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

22 published item(s)

preprint2025arXiv

The tournament ratchet's clicktime process, and metastability in a Moran model

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. In the classical variant, an individual's selective advantage is proportional to the difference between the population average and the individual's mutation load, whereas in the ratchet with {\em tournament selection} only the signs of the differences of the individual mutation loads matter. In a parameter regime which leads to slow clicking (i.e. to a loss of the currently fittest class at a rate $\ll 1/N$) we prove that the rescaled process of click times of the tournament ratchet converges as $N\to \infty$ to a Poisson process. Central ingredients in the proof are a thorough analysis of the metastable behaviour of a two-type Moran model with selection and deleterious mutation (which describes the size of the fittest class up to its extinction time) and a lower estimate on the size of the new fittest class at a clicktime.

preprint2023arXiv

Asymptotic Gaussianity via coalescence probabilities in the Hammond-Sheffield urn

For the renormalised sums of the random $\pm 1$-colouring of the connected components of $\mathbb Z$ generated by the coalescing renewal processes in the "power law Pólya's urn" of Hammond and Sheffield we prove functional convergence towards fractional Brownian motion, closing a gap in the tightness argument of their paper. In addition, in the regime of the strong renewal theorem we gain insights into the coalescing renewal processes in the Hammond-Sheffield urn (such as the asymptotic depth of most recent common ancestors) and are able to control the coalescence probabilities of two, three and four individuals that are randomly sampled from $[n]$. This allows us to obtain a new, conceptual proof of the asymptotic Gaussianity (including the functional convergence) of the renormalised sums of more general colourings, which can be seen as an invariance principle beyond the main result of Hammond and Sheffield. In this proof, a key ingredient of independent interest is a sufficient criterion for the asymptotic Gaussianity of the renormalised sums in randomly coloured random partitions of $[n]$, based on Stein's method. Along the way we also prove a statement on the asymptotics of the coalescence probabilities in the long-range seedbank model of Blath, González Casanova, Kurt, and Spanò.

preprint2022arXiv

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

For a class of Cannings models we prove Haldane's formula, $π(s_N) \sim \frac{2s_N}{ρ^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 $s_N \sim N^{-b}$ and $0< b<1/2$. Here, $s_N$ is the selective advantage of an individual carrying the beneficial type, and $ρ^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.

preprint2020arXiv

Doob--Martin boundary of Rémy's tree growth chain

Rémy's algorithm is a Markov chain that iteratively generates a sequence of random trees in such a way that the $n^{\mathrm{th}}$ tree is uniformly distributed over the set of rooted, planar, binary trees with $2n+1$ vertices. We obtain a concrete characterization of the Doob--Martin boundary of this transient Markov chain and thereby delineate all the ways in which, loosely speaking, this process can be conditioned to "go to infinity" at large times. A (deterministic) sequence of finite rooted, planar, binary trees converges to a point in the boundary if for each $m$ the random rooted, planar, binary tree spanned by $m+1$ leaves chosen uniformly at random from the $n^{\mathrm{th}}$ tree in the sequence converges in distribution as $n$ tends to infinity -- a notion of convergence that is analogous to one that appears in the recently developed theory of graph limits. We show that a point in the Doob--Martin boundary may be identified with the following ensemble of objects: a complete separable $\mathbb{R}$-tree that is rooted and binary in a suitable sense, a diffuse probability measure on the $\mathbb{R}$-tree that allows us to make sense of sampling points from it, and a kernel on the $\mathbb{R}$-tree that describes the probability that the first of a given pair of points is below and to the left of their most recent common ancestor while the second is below and to the right. The Doob--Martin boundary corresponds bijectively to the set of extreme points of the closed convex set of normalized nonnegative harmonic functions, in other words, the minimal and full Doob--Martin boundaries coincide. These results are in the spirit of the identification of graphons as limit objects in the theory of graph limits.

preprint2020arXiv

Haldane's formula in Cannings models: The case of moderately weak selection

We introduce a Cannings model with directional selection via a paintbox construction and establish a strong duality with the line counting process of a new \emph{Cannings ancestral selection graph} in discrete time. This duality also yields a formula for the fixation probability of the beneficial type. Haldane's formula states that for a single selectively advantageous individual in a population of haploid individuals of size $N$ the prob\-ability of fixation is asymptotically (as $N\to \infty$) equal to the selective advantage of haploids $s_N$ divided by half of the offspring variance. For a class of offspring distributions within Kingman attraction we prove this asymptotics for sequences $s_N$ obeying $N^{-1} \ll s_N \ll N^{-1/2} $, which is a regime of "moderately weak selection". It turns out that for $ s_N \ll N^{-2/3} $ the Cannings ancestral selection graph is so close to the ancestral selection graph of a Moran model that a suitable coupling argument allows to play the problem back asymptotically to the fixation probability in the Moran model, which can be computed explicitly.

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.

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

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

The Lenski experiment investigates the long-term evolution of bacterial populations. Its design allows the direct comparison of the reproductive fitness of an evolved strain with its founder ancestor. It was observed by Wiser et al. (2013) that the relative fitness over time increases sublinearly, a behaviour which is commonly attributed to effects like clonal interference or epistasis. In this paper we present an individual-based probabilistic model that captures essential features of the design of the Lenski experiment. We assume that each beneficial mutation increases the individual reproduction rate by a fixed amount, which corresponds to the absence of epistasis in the continuous-time (intraday) part of the model, but leads to an epistatic effect in the discrete-time (interday) part of the model. Using an approximation by near-critical Galton-Watson processes, we prove that under some assumptions on the model parameters which exclude clonal interference, the relative fitness process converges, after suitable rescaling, in the large population limit to a power law function.

preprint2016arXiv

Radix sort trees in the large

The trie-based radix sort algorithm stores pairwise different infinite binary strings in the leaves of a binary tree in a way that the Ulam-Harris coding of each leaf equals a prefix (that is, an initial segment) of the corresponding string, with the prefixes being of minimal length so that they are pairwise different. We investigate the {\em radix sort tree chains} -- the tree-valued Markov chains that arise when successively storing infinite binary strings $Z_1,\ldots, Z_n$, $n=1,2,\ldots$ according to the trie-based radix sort algorithm, where the source strings $Z_1, Z_2,\ldots$ are independent and identically distributed. We establish a bijective correspondence between the full Doob--Martin boundary of the radix sort tree chain with a {\em symmetric Bernoulli source} (that is, each $Z_k$ is a fair coin-tossing sequence) and the family of radix sort tree chains for which the common distribution of the $Z_k$ is a diffuse probability measure on $\{0,1\}^\infty$. In essence, our result characterizes all the ways that it is possible to condition such a chain of radix sort trees consistently on its behavior "in the large".

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.

preprint2016arXiv

The fixation time of a strongly beneficial allele in a structured population

For a beneficial allele which enters a large unstructured population and eventually goes to fixation, it is known that the time to fixation is approximately $2\log(α)/α$ for a large selection coefficient $α$. For a population that is distributed over finitely many colonies, with migration between these colonies, we detect various regimes of the migration rate $μ$ for which the fixation times have different asymptotics as $α\to \infty$. If $μ$ is of order $α$, the allele fixes (as in the spatially unstructured case) in time $\sim 2\log(α)/α$. If $μ$ is of order $α^γ, 0\leq γ\leq 1$, the fixation time is $\sim (2 + (1-γ)Δ) \log(α)/α$, where $Δ$ is the number of migration steps that are needed to reach all other colonies starting from the colony where the beneficial allele appeared. If $μ= 1/\log(α)$, the fixation time is $\sim (2+S)\log(α)/α$, where $S$ is a random time in a simple epidemic model. The main idea for our analysis is to combine a new moment dual for the process conditioned to fixation with the time reversal in equilibrium of a spatial version of Neuhauser and Krone's ancestral selection graph.

preprint2015arXiv

A path-valued Markov process indexed by the ancestral mass

A family of Feller branching diffusions $Z^x$, $x \ge 0$, with nonlinear drift and initial value $x$ can, with a suitable coupling over the {\em ancestral masses} $x$, be viewed as a path-valued process indexed by $x$. For a coupling due to Dawson and Li, which in case of a linear drift describes the corresponding Feller branching diffusion, and in our case makes the path-valued process Markovian, we find an SDE solved by $Z$, which is driven by a random point measure on excursion space. In this way we are able to identify the infinitesimal generator of the path-valued process. We also establish path properties of $x\mapsto Z^x$ using various couplings of $Z$ with classical Feller branching diffusions.

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 Kingman tree length process has infinite quadratic variation

In the case of neutral populations of fixed sizes in equilibrium whose genealogies are described by the Kingman $N$-coalescent back from time $t$ consider the associated processes of total tree length as $t$ increases. We show that the (càdlàg) process to which the sequence of compensated tree length processes converges as $N$ tends to infinity is a process of infinite quadratic variation; therefore this process cannot be a semimartingale. This answers a question posed in Pfaffelhuber et al. (2011).

preprint2014arXiv

The evolving beta coalescent

In mathematical population genetics, it is well known that one can represent the genealogy of a population by a tree, which indicates how the ancestral lines of individuals in the population coalesce as they are traced back in time. As the population evolves over time, the tree that represents the genealogy of the population also changes, leading to a tree-valued stochastic process known as the evolving coalescent. Here we will consider the evolving coalescent for populations whose genealogy can be described by a beta coalescent, which is known to give the genealogy of populations with very large family sizes. We show that as the size of the population tends to infinity, the evolution of certain functionals of the beta coalescent, such as the total number of mergers, the total branch length, and the total length of external branches, converges to a stationary stable process. Our methods also lead to new proofs of known asymptotic results for certain functionals of the non-evolving beta coalescent.

preprint2013arXiv

From Brownian motion with a local time drift to Feller's branching diffusion with logistic growth

We give a new proof for a Ray-Knight representation of Feller's branching diffusion with logistic growth in terms of the local times of a reflected Brownian motion $H$ with a drift that is affine linear in the local time accumulated by $H$ at its current level. In \cite{LPW}, such a representation was obtained by an approximation through Harris paths that code the genealogies of particle systems. The present proof is purely in terms of stochastic analysis, and is inspired by previous work of Norris, Rogers and Williams \cite{NRW}.

preprint2013arXiv

The spatial Lambda-Fleming-Viot process: an event-based construction and a lookdown representation

We construct a measure-valued equivalent to the spatial Lambda-Fleming-Viot process (SLFV) introduced in [Eth08]. In contrast with the construction carried out in [Eth08], we fix the realization of the sequence of reproduction events and obtain a quenched evolution of the local genetic diversities. To this end, we use a particle representation which highlights the role of the genealogies in the attribution of types (or alleles) to the individuals of the population. This construction also enables us to clarify the state-space of the SLFV and to derive several path properties of the measure-valued process as well as of the labeled trees describing the genealogical relations between a sample of individuals. We complement it with a look-down construction which provides a particle system whose empirical distribution at time t, seen as a process in t, has the law of the quenched SLFV. In all these results, the facts that we work with a fixed configuration of events and that reproduction occurs only locally in space introduce serious technical issues that are overcome by controlling the number of events occurring and of particles present in a given area over macroscopic time intervals.

preprint2013arXiv

Trees under attack: a Ray-Knight representation of Feller's branching diffusion with logistic growth

We obtain a representation of Feller's branching diffusion with logistic growth in terms of the local times of a reflected Brownian motion $H$ with a drift that is affine linear in the local time accumulated by $H$ at its current level. As in the classical Ray-Knight representation, the excursions of $H$ are the exploration paths of the trees of descendants of the ancestors at time $t=0$, and the local time of $H$ at height $t$ measures the population size at time $t$ (see e.g. \cite{LG4}). We cope with the dependence in the reproduction by introducing a pecking order of individuals: an individual explored at time $s$ and living at time $t=H_s$ is prone to be killed by any of its contemporaneans that have been explored so far. The proof of our main result relies on approximating $H$ with a sequence of Harris paths $H^N$ which figure in a Ray-Knight representation of the total mass of a branching particle system. We obtain a suitable joint convergence of $H^N$ together with its local times {\em and} with the Girsanov densities that introduce the dependence in the reproduction.

preprint2012arXiv

The total external branch length of Beta-coalescents

For $1<α<2$ we derive the asymptotic distribution of the total length of {\em external} branches of a Beta$(2-α, α)$-coalescent as the number $n$ of leaves becomes large. It turns out the fluctuations of the external branch length follow those of $τ_n^{2-α}$ over the entire parameter regime, where $τ_n$ denotes the random number of coalescences that bring the $n$ lineages down to one. This is in contrast to the fluctuation behavior of the total branch length, which exhibits a transition at $α_0 = (1+\sqrt 5)/2$.

preprint2011arXiv

Trickle-down processes and their boundaries

It is possible to represent each of a number of Markov chains as an evolving sequence of connected subsets of a directed acyclic graph that grow in the following way: initially, all vertices of the graph are unoccupied, particles are fed in one-by-one at a distinguished source vertex, successive particles proceed along directed edges according to an appropriate stochastic mechanism, and each particle comes to rest once it encounters an unoccupied vertex. Examples include the binary and digital search tree processes, the random recursive tree process and generalizations of it arising from nested instances of Pitman's two-parameter Chinese restaurant process, tree-growth models associated with Mallows' phi model of random permutations and with Schuetzenberger's non-commutative q-binomial theorem, and a construction due to Luczak and Winkler that grows uniform random binary trees in a Markovian manner. We introduce a framework that encompasses such Markov chains, and we characterize their asymptotic behavior by analyzing in detail their Doob-Martin compactifications, Poisson boundaries and tail sigma-fields.

preprint2010arXiv

The tree length of an evolving coalescent

A well-established model for the genealogy of a large population in equilibrium is Kingman's coalescent. For the population together with its genealogy evolving in time, this gives rise to a time-stationary tree-valued process. We study the sum of the branch lengths, briefly denoted as tree length, and prove that the (suitably compensated) sequence of tree length processes converges, as the population size tends to infinity, to a limit process with cadlag paths, infinite infinitesimal variance, and a Gumbel distribution as its equilibrium.