Source author record

Jason Schweinsberg

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

ResearcherUnclaimed source record

math.PR Populations and Evolution math-ph math.MP Tissues and Organs

Catalog footprint

What is connected

17works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Particle configurations for branching Brownian motion with an inhomogeneous branching rate

Aiming to understand the distribution of fitness levels of individuals in a large population undergoing selection, we study the particle configurations of branching Brownian motion where each particle independently moves as Brownian motion with negative drift, particles can die or undergo dyadic fission, and the difference between the birth rate and the death rate is proportional to the particle's location. Under some assumptions, we obtain the limit in probability of the number of particles in any given interval and an explicit formula for the asymptotic empirical density of the fitness distribution. We show that after a sufficiently long time, the fitness distribution from the lowest to the highest fitness levels approximately evolves as a traveling wave with a profile which is asymptotically related the the Airy function. Our work complements the results in Roberts and Schweinsberg (2021), giving a fuller picture of the fitness distribution.

preprint2021arXiv

A Gaussian particle distribution for branching Brownian motion with an inhomogeneous branching rate

Motivated by the goal of understanding the evolution of populations undergoing selection, we consider branching Brownian motion in which particles independently move according to one-dimensional Brownian motion with drift, each particle may either split into two or die, and the difference between the birth and death rates is a linear function of the position of the particle. We show that, under certain assumptions, after a sufficiently long time, the empirical distribution of the positions of the particles is approximately Gaussian. This provides mathematically rigorous justification for results in the biology literature indicating that the distribution of the fitness levels of individuals in a population over time evolves like a Gaussian traveling wave.

preprint2020arXiv

$Λ$-coalescents arising in populations with dormancy

Consider a population evolving from year to year through three seasons: spring, summer and winter. Every spring starts with $N$ dormant individuals waking up independently of each other according to a given distribution. Once an individual is awake, it starts reproducing at a constant rate. By the end of spring, all individuals are awake and continue reproducing independently as Yule processes during the whole summer. In the winter, $N$ individuals chosen uniformly at random go to sleep until the next spring, and the other individuals die. We show that because an individual that wakes up unusually early can have a large number of surviving descendants, for some choices of model parameters the genealogy of the population will be described by a $Λ$-coalescent. In particular, the beta coalescent can describe the genealogy when the rate at which individuals wake up increases exponentially over time. We also characterize the set of all $Λ$-coalescents that can arise in this framework.

preprint2020arXiv

Mutation timing in a spatial model of evolution

Motivated by models of cancer formation in which cells need to acquire $k$ mutations to become cancerous, we consider a spatial population model in which the population is represented by the $d$-dimensional torus of side length $L$. Initially, no sites have mutations, but sites with $i-1$ mutations acquire an $i$th mutation at rate $μ_i$ per unit area. Mutations spread to neighboring sites at rate $α$, so that $t$ time units after a mutation, the region of individuals that have acquired the mutation will be a ball of radius $αt$. We calculate, for some ranges of the parameter values, the asymptotic distribution of the time required for some individual to acquire $k$ mutations. Our results, which build on previous work of Durrett, Foo, and Leder, are essentially complete when $k = 2$ and when $μ_i = μ$ for all $i$.

preprint2015arXiv

Rigorous results for a population model with selection I: evolution of the fitness distribution

We consider a model of a population of fixed size $N$ undergoing selection. Each individual acquires beneficial mutations at rate $μ_N$, and each beneficial mutation increases the individual's fitness by $s_N$. Each individual dies at rate one, and when a death occurs, an individual is chosen with probability proportional to the individual's fitness to give birth. Under certain conditions on the parameters $μ_N$ and $s_N$, we obtain rigorous results for the rate at which mutations accumulate in the population and the distribution of the fitnesses of individuals in the population at a given time. Our results confirm predictions of Desai and Fisher (2007).

preprint2015arXiv

Rigorous results for a population model with selection II: genealogy of the population

We consider a model of a population of fixed size $N$ undergoing selection. Each individual acquires beneficial mutations at rate $μ_N$, and each beneficial mutation increases the individual's fitness by $s_N$. Each individual dies at rate one, and when a death occurs, an individual is chosen with probability proportional to the individual's fitness to give birth. Under certain conditions on the parameters $μ_N$ and $s_N$, we show that the genealogy of the population can be described by the Bolthausen-Sznitman coalescent. This result confirms predictions of Desai, Walczak, and Fisher (2013), and Neher and Hallatschek (2013).

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

Critical branching Brownian motion with absorption: particle configurations

We consider critical branching Brownian motion with absorption, in which there is initially a single particle at $x > 0$, particles move according to independent one-dimensional Brownian motions with the critical drift of $-\sqrt{2}$, and particles are absorbed when they reach zero. Here we obtain asymptotic results concerning the behavior of the process before the extinction time, as the position $x$ of the initial particle tends to infinity. We estimate the number of particles in the system at a given time and the position of the right-most particle. We also obtain asymptotic results for the configuration of particles at a typical time.

preprint2013arXiv

The genealogy of branching Brownian motion with absorption

We consider a system of particles which perform branching Brownian motion with negative drift and are killed upon reaching zero, in the near-critical regime where the total population stays roughly constant with approximately N particles. We show that the characteristic time scale for the evolution of this population is of order $(\log N)^3$, in the sense that when time is measured in these units, the scaled number of particles converges to a variant of Neveu's continuous-state branching process. Furthermore, the genealogy of the particles is then governed by a coalescent process known as the Bolthausen-Sznitman coalescent. This validates the nonrigorous predictions by Brunet, Derrida, Muller and Munier for a closely related model.

preprint2012arXiv

Critical branching Brownian motion with absorption: survival probability

We consider branching Brownian motion on the real line with absorption at zero, in which particles move according to independent Brownian motions with the critical drift of $-\sqrt{2}$. Kesten (1978) showed that almost surely this process eventually dies out. Here we obtain upper and lower bounds on the probability that the process survives until some large time $t$. These bounds improve upon results of Kesten (1978), and partially confirm nonrigorous predictions of Derrida and Simon (2007).

preprint2011arXiv

Consensus in the two-state Axelrod model

The Axelrod model is a spatial stochastic model for the dynamics of cultures which, similarly to the voter model, includes social influence, but differs from the latter by also accounting for another social factor called homophily, the tendency to interact more frequently with individuals who are more similar. Each individual is characterized by its opinions about a finite number of cultural features, each of which can assume the same finite number of states. Pairs of adjacent individuals interact at a rate equal to the fraction of features they have in common, thus modeling homophily, which results in the interacting pair having one more cultural feature in common, thus modeling social influence. It has been conjectured based on numerical simulations that the one-dimensional Axelrod model clusters when the number of features exceeds the number of states per feature. In this article, we prove this conjecture for the two-state model with an arbitrary number of features.

preprint2011arXiv

Dynamics of the evolving Bolthausen-Sznitman coalescent

Consider a population of fixed size that evolves over time. At each time, the genealogical structure of the population can be described by a coalescent tree whose branches are traced back to the most recent common ancestor of the population. As time goes forward, the genealogy of the population evolves, leading to what is known as an evolving coalescent. We will study the evolving coalescent for populations whose genealogy can be described by the Bolthausen-Sznitman coalescent. We obtain the limiting behavior of the evolution of the time back to the most recent common ancestor and the total length of the branches in the tree. By similar methods, we also obtain a new result concerning the number of blocks in the Bolthausen-Sznitman coalescent.

preprint2010arXiv

Survival of near-critical branching Brownian motion

Consider a system of particles performing branching Brownian motion with negative drift $μ= \sqrt{2 - ε}$ and killed upon hitting zero. Initially there is one particle at $x>0$. Kesten showed that the process survives with positive probability if and only if $ε>0$. Here we are interested in the asymptotics as $\eps\to 0$ of the survival probability $Q_μ(x)$. It is proved that if $L= π/\sqrtε$ then for all $x \in \R$, $\lim_{ε\to 0} Q_μ(L+x) = θ(x) \in (0,1)$ exists and is a travelling wave solution of the Fisher-KPP equation. Furthermore, we obtain sharp asymptotics of the survival probability when $x<L$ and $L-x \to \infty$. The proofs rely on probabilistic methods developed by the authors in a previous work. This completes earlier work by Harris, Harris and Kyprianou and confirms predictions made by Derrida and Simon, which were obtained using nonrigorous PDE methods.

preprint2010arXiv

The number of small blocks in exchangeable random partitions

Suppose $Π$ is an exchangeable random partition of the positive integers and $Π_n$ is its restriction to $\{1, ..., n\}$. Let $K_n$ denote the number of blocks of $Π_n$, and let $K_{n,r}$ denote the number of blocks of $Π_n$ containing $r$ integers. We show that if $0 < α< 1$ and $K_n/(n^α \ell(n))$ converges in probability to $Γ(1-α)$, where $\ell$ is a slowly varying function, then $K_{n,r}/(n^α \ell(n))$ converges in probability to $αΓ(r - α)/r!$. This result was previously known when the convergence of $K_n/(n^α \ell(n))$ holds almost surely, but the result under the hypothesis of convergence in probability has significant implications for coalescent theory. We also show that a related conjecture for the case when $K_n$ grows only slightly slower than $n$ fails to be true.

preprint2009arXiv

A waiting time problem arising from the study of multi-stage carcinogenesis

We consider the population genetics problem: how long does it take before some member of the population has $m$ specified mutations? The case $m=2$ is relevant to onset of cancer due to the inactivation of both copies of a tumor suppressor gene. Models for larger $m$ are needed for colon cancer and other diseases where a sequence of mutations leads to cells with uncontrolled growth.

preprint2008arXiv

The waiting time for m mutations

We consider a model of a population of fixed size N in which each individual gets replaced at rate one and each individual experiences a mutation at rate μ. We calculate the asymptotic distribution of the time that it takes before there is an individual in the population with m mutations. Several different behaviors are possible, depending on how μchanges with N. These results have applications to the problem of determining the waiting time for regulatory sequences to appear and to models of cancer development.

preprint2007arXiv

Beta-coalescents and continuous stable random trees

Coalescents with multiple collisions, also known as $Λ$-coalescents, were introduced by Pitman and Sagitov in 1999. These processes describe the evolution of particles that undergo stochastic coagulation in such a way that several blocks can merge at the same time to form a single block. In the case that the measure $Λ$ is the $\operatorname {Beta}(2-α,α)$ distribution, they are also known to describe the genealogies of large populations where a single individual can produce a large number of offspring. Here, we use a recent result of Birkner et al. to prove that Beta-coalescents can be embedded in continuous stable random trees, about which much is known due to the recent progress of Duquesne and Le Gall. Our proof is based on a construction of the Donnelly--Kurtz lookdown process using continuous random trees, which is of independent interest. This produces a number of results concerning the small-time behavior of Beta-coalescents. Most notably, we recover an almost sure limit theorem of the present authors for the number of blocks at small times and give the multifractal spectrum corresponding to the emergence of blocks with atypical size. Also, we are able to find exact asymptotics for sampling formulae corresponding to the site frequency spectrum and the allele frequency spectrum associated with mutations in the context of population genetics.