-
Bernstein duality revisited: frequency-dependent selection, coordinated mutation and opposing environments
Authors:
Fernando Cordero,
Sebastian Hummel,
Grégoire Véchambre
Abstract:
Population models usually come in pairs; one process describes forward evolution (e.g. type composition) and the other describes backward evolution (e.g. lines of descent). These processes are often linked by a formal relationship known as duality. Ideally, one of the two processes is easier to analyze, and the duality relation is so simple that properties of the more involved process can be infer…
▽ More
Population models usually come in pairs; one process describes forward evolution (e.g. type composition) and the other describes backward evolution (e.g. lines of descent). These processes are often linked by a formal relationship known as duality. Ideally, one of the two processes is easier to analyze, and the duality relation is so simple that properties of the more involved process can be inferred from the simpler one. This is the case when the forward process admits a moment dual. Unfortunately, moment duality seems to be the exception rather than the rule. Various approaches have been used to analyze models in the absence of a moment dual, one of them is based on Bernstein duality and another one on Siegmund duality. As a rule of thumb, the first approach seems to work well whenever the ancestral processes are positive recurrent; the second one, in contrast, works well in situations where the ancestral structures can grow to infinity (in size). The second approach was recently used to provide a full characterization of the long-term behavior of a broad class of $Λ$-Wright--Fisher processes subject to frequency-dependent selection and opposing environments. In this paper, we use the first approach to complete the picture, i.e. we describe the long-term behavior of a different class of $Λ$-Wright--Fisher processes, which covers many of the cases that were not covered by the aforementioned result (the two classes intersect, but none is a proper subset of the other one). Moreover, we extend the notion of Bernstein duality to cases with (single and coordinated) mutations {and environmental selection}, and we use it to show ergodic properties of the process.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Constructing Maximal Germ Couplings of Brownian Motions with Drift
Authors:
Sebastian Hummel,
Adam Quinn Jaffe
Abstract:
Consider all the possible ways of coupling together two Brownian motions with the same starting position but with different drifts onto the same probability space. It is known that there exist couplings which make these processes agree for some random, positive, maximal initial length of time. Presently, we provide an explicit, elementary construction of such couplings.
Consider all the possible ways of coupling together two Brownian motions with the same starting position but with different drifts onto the same probability space. It is known that there exist couplings which make these processes agree for some random, positive, maximal initial length of time. Presently, we provide an explicit, elementary construction of such couplings.
△ Less
Submitted 28 September, 2023;
originally announced September 2023.
-
Absorption and stationary times for the $Λ$-Wright-Fisher process
Authors:
Airam Blancas,
Adrián González Casanova,
Sebastian Hummel,
Sandra Palau
Abstract:
We derive stationary and fixation times for the multi-type $Λ$-Wright-Fisher process with and without the classic linear drift that models mutations. Our method relies on a grand coupling of the process realized through the so-called lookdown-construction. A well-known process embedded in this construction is the fixation line. We generalise the process to our setup and make use of the associated…
▽ More
We derive stationary and fixation times for the multi-type $Λ$-Wright-Fisher process with and without the classic linear drift that models mutations. Our method relies on a grand coupling of the process realized through the so-called lookdown-construction. A well-known process embedded in this construction is the fixation line. We generalise the process to our setup and make use of the associated explosion times to obtain a representation of the fixation and stationary times in terms of the waiting time in a coupon collector problem.
△ Less
Submitted 17 August, 2023;
originally announced August 2023.
-
Mean-field interacting multi-type birth-death processes with a view to applications in phylodynamics
Authors:
William S. DeWitt,
Steven N. Evans,
Ella Hiesmayr,
Sebastian Hummel
Abstract:
Multi-type birth-death processes underlie approaches for inferring evolutionary dynamics from phylogenetic trees across biological scales, ranging from deep-time species macroevolution to rapid viral evolution and somatic cellular proliferation. A limitation of current phylogenetic birth-death models is that they require restrictive linearity assumptions that yield tractable message-passing likeli…
▽ More
Multi-type birth-death processes underlie approaches for inferring evolutionary dynamics from phylogenetic trees across biological scales, ranging from deep-time species macroevolution to rapid viral evolution and somatic cellular proliferation. A limitation of current phylogenetic birth-death models is that they require restrictive linearity assumptions that yield tractable message-passing likelihoods, but that also preclude interactions between individuals. Many fundamental evolutionary processes -- such as environmental carrying capacity or frequency-dependent selection -- entail interactions, and may strongly influence the dynamics in some systems. Here, we introduce a multi-type birth-death process in mean-field interaction with an ensemble of replicas of the focal process. We prove that, under quite general conditions, the ensemble's stochastically evolving interaction field converges to a deterministic trajectory in the limit of an infinite ensemble. In this limit, the replicas effectively decouple, and self-consistent interactions appear as nonlinearities in the infinitesimal generator of the focal process. We investigate a special case that is rich enough to model both carrying capacity and frequency-dependent selection while yielding tractable message-passing likelihoods in the context of a phylogenetic birth-death model.
△ Less
Submitted 31 March, 2024; v1 submitted 12 July, 2023;
originally announced July 2023.
-
$Λ$-Wright--Fisher processes with general selection and opposing environmental effects: fixation and coexistence
Authors:
Fernando Cordero,
Sebastian Hummel,
Grégoire Véchambre
Abstract:
Our results characterize the long-term behavior for a broad class of $Λ$-Wright--Fisher processes with frequency-dependent and environmental selection. In particular, we reveal a rich variety of parameter-dependent behaviors and provide explicit criteria to discriminate between them. That includes the situation in which the (entire) boundary is repelling -- a new phenomenon in this context. This h…
▽ More
Our results characterize the long-term behavior for a broad class of $Λ$-Wright--Fisher processes with frequency-dependent and environmental selection. In particular, we reveal a rich variety of parameter-dependent behaviors and provide explicit criteria to discriminate between them. That includes the situation in which the (entire) boundary is repelling -- a new phenomenon in this context. This has significant biological implications, because it means that selection alone can maintain coexistence. If a boundary point is attractive, we derive polynomial/exponential decay rates for the probability of not being polynomially/exponentially close to that boundary, depending on some weak/strong integrability conditions. Moreover, we provide a handy representation of the fixation probability. In our proofs we make use of Siegmund duality. The dual process can be sandwiched near the boundaries in between transformed Lévy processes. In this way we relate the boundary behavior of the dual process to fluctuation properties of these Lévy processes and shed new light on previously established conditions for attractive/repelling boundary points. Our method allows us to treat models that so far could not be analyzed by means of moment or Bernstein duality. This closes an existing gap in the literature.
△ Less
Submitted 26 February, 2024; v1 submitted 20 December, 2021;
originally announced December 2021.
-
Lines of descent in a Moran model with frequency-dependent selection and mutation
Authors:
Ellen Baake,
Luigi Esercito,
Sebastian Hummel
Abstract:
We study ancestral structures for the two-type Moran model with mutation and frequency-dependent selection under the nonlinear dominance or fittest-type-wins scheme. Under appropriate conditions, both lead, in distribution, to the same type-frequency process. Reasoning through the mutations on the ancestral selection graph (ASG), we develop the corresponding killed and pruned lookdown ASG and use…
▽ More
We study ancestral structures for the two-type Moran model with mutation and frequency-dependent selection under the nonlinear dominance or fittest-type-wins scheme. Under appropriate conditions, both lead, in distribution, to the same type-frequency process. Reasoning through the mutations on the ancestral selection graph (ASG), we develop the corresponding killed and pruned lookdown ASG and use them to determine the present and ancestral type distributions. To this end, we establish factorial moment dualities to the Moran model and a relative. We extend the results to the diffusion limit and present applications for finite population size as well as moderate and weak selection.
△ Less
Submitted 21 October, 2022; v1 submitted 17 November, 2020;
originally announced November 2020.
-
General selection models: Bernstein duality and minimal ancestral structures
Authors:
Fernando Cordero,
Sebastian Hummel,
Emmanuel Schertzer
Abstract:
$Λ$-Wright--Fisher processes provide a robust framework to describe the type-frequency evolution of an infinite neutral population. We add a polynomial drift to the corresponding stochastic differential equation to incorporate frequency-dependent selection. A decomposition of the drift allows us to approximate the solution of the stochastic differential equation by a sequence of Moran models. The…
▽ More
$Λ$-Wright--Fisher processes provide a robust framework to describe the type-frequency evolution of an infinite neutral population. We add a polynomial drift to the corresponding stochastic differential equation to incorporate frequency-dependent selection. A decomposition of the drift allows us to approximate the solution of the stochastic differential equation by a sequence of Moran models. The genealogical structure underlying the Moran model leads in the large population limit to a generalisation of the ancestral selection graph of Krone and Neuhauser. Building on this object, we construct a continuous-time Markov chain and relate it to the forward process via a new form of duality, which we call Bernstein duality. We adapt classical methods based on the moment duality to determine the time to absorption and criteria for the accessibility of the boundaries; this extends a recent result by González Casanova and Spanò. An intriguing feature of the construction is that the same forward process is compatible with multiple backward models. In this context we introduce suitable notions for minimality among the ancestral processes and characterise the corresponding parameter sets. In this way we recover classic ancestral structures as minimal ones.
△ Less
Submitted 10 May, 2021; v1 submitted 15 March, 2019;
originally announced March 2019.
-
Lines of descent in the deterministic mutation-selection model with pairwise interaction
Authors:
Ellen Baake,
Fernando Cordero,
Sebastian Hummel
Abstract:
We consider the mutation--selection differential equation with pairwise interaction (or, equivalently, the diploid mutation--selection equation) and establish the corresponding ancestral process, which is a random tree and a variant of the ancestral selection graph. The formal relation to the forward model is given via duality. To make the tree tractable, we prune branches upon mutations, thus red…
▽ More
We consider the mutation--selection differential equation with pairwise interaction (or, equivalently, the diploid mutation--selection equation) and establish the corresponding ancestral process, which is a random tree and a variant of the ancestral selection graph. The formal relation to the forward model is given via duality. To make the tree tractable, we prune branches upon mutations, thus reducing it to its informative parts. The hierarchies inherent in the tree are encoded systematically via tripod trees with weighted leaves; this leads to the stratified ancestral selection graph. The latter also satisfies a duality relation with the mutation--selection equation. Each of the dualities provides a stochastic representation of the solution of the differential equation. This allows us to connect the equilibria and their bifurcations to the long-term behaviour of the ancestral process. Furthermore, with the help of the stratified ancestral selection graph, we obtain explicit results about the ancestral type distribution in the case of unidirectional mutation.
△ Less
Submitted 5 May, 2021; v1 submitted 3 December, 2018;
originally announced December 2018.
-
A probabilistic view on the deterministic mutation-selection equation: dynamics, equilibria, and ancestry via individual lines of descent
Authors:
Ellen Baake,
Fernando Cordero,
Sebastian Hummel
Abstract:
We reconsider the deterministic haploid mutation-selection equation with two types. This is an ordinary differential equation that describes the type distribution (forward in time) in a population of infinite size. This paper establishes ancestral (random) structures inherent in this deterministic model. In a first step, we obtain a representation of the deterministic equation's solution (and, in…
▽ More
We reconsider the deterministic haploid mutation-selection equation with two types. This is an ordinary differential equation that describes the type distribution (forward in time) in a population of infinite size. This paper establishes ancestral (random) structures inherent in this deterministic model. In a first step, we obtain a representation of the deterministic equation's solution (and, in particular, of its equilibrium) in terms of an ancestral process called the killed ancestral selection graph. This representation allows one to understand the bifurcations related to the error threshold phenomenon from a genealogical point of view. Next, we characterise the ancestral type distribution by means of the pruned lookdown ancestral selection graph and study its properties at equilibrium. We also provide an alternative characterisation in terms of a piecewise-deterministic Markov process. Throughout, emphasis is on the underlying dualities as well as on explicit results.
△ Less
Submitted 28 March, 2018; v1 submitted 12 October, 2017;
originally announced October 2017.
-
On the Borel Inseparability of Game Tree Languages
Authors:
Szczepan Hummel,
Henryk Michalewski,
Damian Niwinski
Abstract:
The game tree languages can be viewed as an automata-theoretic counterpart of parity games on graphs. They witness the strictness of the index hierarchy of alternating tree automata, as well as the fixed-point hierarchy over binary trees. We consider a game tree language of the first non-trivial level, where Eve can force that 0 repeats from some moment on, and its dual, where Adam can force tha…
▽ More
The game tree languages can be viewed as an automata-theoretic counterpart of parity games on graphs. They witness the strictness of the index hierarchy of alternating tree automata, as well as the fixed-point hierarchy over binary trees. We consider a game tree language of the first non-trivial level, where Eve can force that 0 repeats from some moment on, and its dual, where Adam can force that 1 repeats from some moment on. Both these sets (which amount to one up to an obvious renaming) are complete in the class of co-analytic sets. We show that they cannot be separated by any Borel set, hence {\em a fortiori} by any weakly definable set of trees. This settles a case left open by L.Santocanale and A.Arnold, who have thoroughly investigated the separation property within the $μ$-calculus and the automata index hierarchies. They showed that separability fails in general for non-deterministic automata of type $Σ^μ_{n} $, starting from level $n=3$, while our result settles the missing case $n=2$.
△ Less
Submitted 10 February, 2009;
originally announced February 2009.