Cannings models, population size changes and multiple-merger coalescents

Multiple-merger coalescents, e.g. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document}Λ-n-coalescents, have been proposed as models of the genealogy of n sampled individuals for a range of populations whose genealogical structures are not captured well by Kingman’s n-coalescent. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document}Λ-n-coalescents can be seen as the limit process of the discrete genealogies of Cannings models with fixed population size, when time is rescaled and population size \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\rightarrow \infty $$\end{document}N→∞. As established for Kingman’s n-coalescent, moderate population size fluctuations in the discrete population model should be reflected by a time-change of the limit coalescent. For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document}Λ-n-coalescents, this has been explicitly shown for only a limited subclass of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document}Λ-n-coalescents and exponentially growing populations. This article gives a more general construction of time-changed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document}Λ-n-coalescents as limits of specific Cannings models with rather arbitrary time changes.


Introduction
The genealogies of samples from populations with highly variant offspring numbers, for instance due to sweepstake reproduction or rapid selection, are not well modelled by Kingman's n-coalescent. As a more realistic alternative, multiple-merger coalescents, especially Λ-coalescents have been proposed, as reviewed in Tellier and Lemaire (2014), Irwin et al. (2016) and Eldon et al. (2016). Λ-n-coalescents, introduced by Pitman (1999), Sagitov (1999), Donnelly and Kurtz (1999), are Markovian processes (Π t ) t≥0 , which describe the genealogy of a set of individuals {1, . . . , n}. This is done by representing the ancestral lineages present at time t of these individuals by the sets of offspring of each ancestral lineage in the sample. Thus, (Π t ) t≥0 can be defined as a random process with states in the set of partitions of {1, . . . , n} and transitions via merging of blocks (i.e. merging of ancestral lineages to a common ancestor). For a Λ-n-coalescent, the infinitesimal rates of any merger of k of b present lineages is given by λ b,k := where Λ is a finite measure on [0, 1]. This includes Kingman's n-coalescent if Λ is the Dirac measure in 0. As in the case of Kingman's n-coalescent being the limit genealogy from samples taken from a discrete Wright-Fisher or Moran model, Λ-n-coalescents can be constructed as the (weak) limit of genealogies from samples of size n taken from Cannings models. The limit is reached as population size N goes to infinity and time is rescaled, see Möhle and Sagitov (2001). Time is rescaled by using [c −1 N ] generations in the discrete model as one unit of evolutionary (coalescent) time in the limit, where c N is the probability that two individuals picked in a generation have the same parent one generation before. In the discrete models, the population size N is fixed across all generations. Only populations in an equilibrium state are described well by models with fixed population sizes. This idealized condition often does not apply to natural populations. In particular, due to fluctuating environmental conditions population sizes are expected to fluctuate likewise. Two standard models of population size changes are timespans of exponential growth or decline, as well as population bottlenecks, where population size drops to a fixed size smaller than N for a timespan on the evolutionary (coalescent) timescale. Such changes are featured in coalescent simulators as ms (Hudson 2002) or msprime (Kelleher et al. 2016). The latter changes are also the model of population size changes in PSMC (Li and Durbin 2011) or similar approaches as SMC++ (Terhorst et al. 2017). For the Wright-Fisher model, which converges to Kingman's n-coalescent if population size N is fixed for all generations, the same scaling c −1 N from discrete genealogy to limit is valid for population size changes which maintain a population size of order N at all times, see Griffiths and Tavare (1994) or Kaj and Krone (2003). The resulting limit process is Kingman's n-coalescent, whose timescale is (non-linearly) transformed. However, size changes too extreme can yield a non-bifurcating (multiple merger) genealogy, see Birkner et al. (2009, Sect. 6.1). For Λ-n-coalescents, the link between fluctuating population sizes in the discrete models and the time-change in the coalescent limit is somewhat less established. While conditions for convergence of the discrete genealogies to a limit process are given in Möhle (2002), no explicit construction of haploid Cannings models leading to an analogous limit, a Λ-n-coalescent with changed time scale, is given. For a specific case, the Dirac n-coalescent for an exponentially growing population, such a construction has been given in Matuszewski et al. (2017), based on the fixed-N Cannings model (modified Moran model) from Eldon and Wakeley (2006). However, also other Λ-n-coalescents (or Cannings models which should converge to these) with changed time scale have been recently discussed and applied as models of genealogies, see Spence et al. (2016), Kato et al. (2017), Alter and Louzoun (2016) and Hoscheit and Pybus (2019). This leads to the goal of this article, which is to extend the approach in Matuszewski et al. (2017) to explicitly give a construction of timechanged Λ-coalescents as limits of Cannings models with fluctuating population sizes. The Cannings models used are modified Moran models, see e.g. Huillet and Möhle (2013), and the Cannings models introduced in Schweinsberg (2003). The main tool to establish the convergence to the time-changed Λ-n-coalescent is, as in Matuszewski et al. (2017), applying Möhle (2002, Thm. 2.2). For diploid Cannings models, the umbrella model from Koskela and Wilke Berenguer (2019) gives a general framework to add population size changes, selection, recombination and population structure to the fixed-N -model. There, if one only considers population size changes, the limit is a time-changed Ξ -n-coalescent, a coalescent process with simultaneous multiple mergers. The focus in the present paper is slightly different though, the aim is to explicitly construct Cannings models that converge, after linear time scaling, to a time-changed Λ-n-coalescent, while Koskela and Wilke Berenguer (2019) concentrates on the convergence itself.

Models and main results
Cannings models (Cannings 1974(Cannings , 1975 describe the probabilistic structure of the pedigree (offspring-parent relations) of a finite population in generations and offspring sizes are exchangeable, i.e. . . . , ν (v) σ (N v ) ) for any permutation σ ∈ S N v . The offspring generation v + 1 then consists of these individuals in arbitrary order (independent of the parents). The case N v = N for all v is denoted as the fixed-N case. From now on, look at the genealogy of the population in generation 0. For convenience, denote the generations in reverse order by r = −v, i.e. if one looks i generations back, this is denoted by r = i. The population sizes N r are defined relative to a reference size N , in a way that if N → ∞, also N r → ∞. From now on, use N = N 0 . The goal is to establish a limit process of the discrete genealogies as N → ∞. The discrete genealogy of a sample of size n in generation 0 is a random process (R (N ) r ) r ∈N 0 with values in the partitions of {1, . . . , n}, where i, j are in the same block of R (N ) r iff they share the same ancestor in generation r . The terminology from Möhle (2002) is used with slight adaptations. Let c N ,r be the probability that two arbitrary individuals in generation r − 1 have the same ancestor in generation r in the model with reference population size N . To clarify, c N ,r is the coalescence probability for individuals in generation r − 1 if population sizes are variable, while c N denotes the coalescence probability in the fixed-N case. Define r =1 c N ,r and let be its shifted pseudo-inverse. For l and a 1 , . . . , a l ≥ 1, set as the probability that in generation r − 1, from l i=1 a i individuals sampled from the Cannings model, specific sets of a 1 ≥ · · · ≥ a l individuals each find a common ancestor one generation before (generation r ), where ancestors of different sets are different. For l = 1, a 1 = 2, c N ,r = Φ (N ) 1 (r ; 2). See Möhle (1998) for details. Consider a sequence of fixed-N Cannings models for each N → ∞ with c N → 0 as N → ∞ and transition probabilities Φ (N ) l (a 1 , . . . , a l ) for a merger of a 1 , . . . , a l ≥ 1 individuals , converging to a Λ-n-coalescent (Π t ) t≥0 with infinitesimal transition rates φ l (a 1 , . . . , a l ) for a 1 ≥ · · · ≥ a l ≥ 2 as N → ∞, see (Möhle and Sagitov 2001, Thm. 2.1). We will establish a variant of Möhle (2002, Corollary 2.4) to show convergence of a variety of Cannings models with variable population sizes to (time-changed) Λ-n-coalescents. For this, we need some assumptions. Most importantly, an asymptotically infinite sum needs to be controlled. For this, we introduce a concept of o-terms: if summing the first O(x −1 N ) members of y N still vanishes as N → ∞. This will be denoted by y N For instance, this condition is satisfied if y N (r )/x N have the same null sequence majorant as N → ∞ for all r ∈ N. Fix t > 0. We assume for all 0 ≤ r ≤ [c −1 N t]: -Population size changes of order N leading to a well-defined population size profile in coalescent time, i.e.
for positive and finite functions c 1 , c 2 , ν : The first class of Cannings models used to construct time-changed Λ-n-coalescents are modified Moran models. In a modified Moran model, only a single individual has more than one offspring (and may have many offspring). Following Huillet and Möhle (2013), define the modified (haploid) Moran model with fixed population size N . Let (U N (z)) z∈N be i.i.d. random variables with values in {2, . . . , N }, let U N be a r.v. with their common distribution. In each generation z ∈ Z, -One randomly chosen individual has U N (z) offspring, -U N (z) − 1 randomly chosen individuals have no offspring, -The other N − U N (z) have one offspring each, Specific modified Moran models leading to Dirac n-coalescents as genealogy limits have been introduced as population models with skewed offspring distributions, see Eldon and Wakeley (2006) and Matuszewski et al. (2017), for fixed and variable population sizes. For any Λ-n-coalescent with Λ([0, 1]) = 1 (denoted by Λ ∈ M[0, 1]), Möhle (2013, Prop. 3.4) show that there always exist fixed-N modified Moran models such that their rescaled genealogies converge to the Λ-n-coalescent. These can be constructed via a random variable U N , that is distributed like the merger size of the first merger in a Λ-N -coalescent. As shown in (Huillet and Möhle 2013, Eq. 9), this means where λ N is the total transition rate of the Λ-N -coalescent and X has distribution Λ.
To add population size changes from generation to generation to the modified Moran model, the relationship between offspring and parent generation needs to be defined. This will be done by adjusting the fixed-N model: If in generation r , there are N r individuals, first run a fixed-N r modified Moran model, producing N r (potential) offspring. Let U n,r denote the number of offspring in generation r −1 of the multiplying parent in generation r . If population size declines from generation r to r − 1, sample N r −1 individuals randomly (without replacement) from the N r potential offspring consisting of U N r offspring of the multiplying parent and N r − U N r single offspring. If individuals are added to the population (population growth), i.e.
to still end up with a modified Moran model one has two options. Additional individuals can be added as further offspring of the already multiplying parent. A second option is to add individuals as offspring of the non-reproducing individuals from generation r in the fixed-N r model. Each originally non-reproducing parent can have one offspring, so this allows one to add U N r − 1 ≥ 1 individuals. The number of additional individuals can be divided between these two options, let A n,r denote the individuals added as Fig. 1 Example of allocation of individuals when population size is increasing. Left: start with a fixed size Moran model with U 6 = 4. Right: population increases by d 6,r = 3, from which A n,r = 1 individual is allocated to the multiplying parent from generation r in the fixed size model (and 2 to non-reproducing individuals from the fixed size model in generation r ). This results in U 6,r = 5 offspring to the multiplying parent (which means that d N ,r − A N ,r individuals are added as offspring of non-reproducing parents from the fixed-N r model). Expressed differently, there are N r −1 = N r + d N ,r offspring from which U n,r = U N r + A N ,r share the same parent, while all other offspring are single offspring of other parents (who all differ). See Fig. 1 for an example.
While some care has to be taken to not change coalescence probabilities (see Remark 5 for an example), there will be different possibilities to choose A N ,r . For Dirac-n-coalescents with exponential growth (on the coalescent time scale), Matuszewski et al. (2017) used A n,r = d N ,r . A reasonable approach may also be to set A n,r (close to) proportional to the fraction U N r /N r of offspring coming from the multiplying parent: Each of the d N ,r added individuals are added to the multiplying parent with probability U N r /N r (with the obvious constraint that after U N r − 1 individuals are added as offspring of non-reproducing parents, all further individuals need to be added to the multiplying parent). As for the fixed-size models, we consider the genealogy of a sample of n individuals, which is denoted by (R (N ) r ) r ∈N 0 The main results of the present paper show that the two allocation schemes allow one to construct Λ-n-coalescent limits of the genealogies of these modified Moran models if population sizes vary in the discrete models in ways described by Eq. (4).
Define a modified Moran model for fixed N by

In the discrete model, additional individuals can be added in any way so that the resulting model is still a modified Moran model.
For Λ not covered by Theorem 1, one can choose slightly different modified Moran models that converge to a Λ-n-coalescent limit for an arbitrary population size profile on the coalescent time scale.
In the discrete model, additional individuals are added solely as offspring of non-reproducing parents from the fixed-N r model, unless E((U N ) 2 ) → ∞ as N → ∞. In that case, they can be added any way that preserves that the model is still a modified Moran model.

Remark 1
The condition of Λ([0, 1]) = 1 in both theorems is not very important: If one scales by c 2 c N instead of c N for any c 2 > 0, the rescaled discrete genealogies converge to the c 2 Λ-n-coalescent.
To get convergence to a time-changed Λ-n-coalescent in Theorem 2, i.e.

R (N )
in the Skorohod-sense as N → ∞, one needs that exists for all t ≥ 0. The following corollary shows that, at least for certain measures Λ, this condition is met.
The specific models used in each of the two theorems are not the only possibilities of modified Moran models with variable population sizes to converge to Λ-n-coalescents. For instance, if one only allows certain population size changes, one can also use the modified Moran model used in Theorem 2 for some Λ covered by Theorem 1.
Consider an exponentially growing modified Moran model population on the coalescent time scale, i.e. ν(t) = exp(−ρt) Then, there exist population sizes satisfying Eq. (4) for ν so that the genealogies (R In the discrete model, additional individuals can be added in any way so that the resulting model is still a modified Moran model.

Individuals are added only as offspring of non-reproducing parents (in the fixed-N r model) if the population size increases.
For Beta(2 − α, α)-n-coalescents, 1 ≤ α < 2, genealogies sampled from the fixed-N Cannings models introduced in Schweinsberg (2003) also converge weakly to these Beta coalescent processes (after rescaling of time) as N → ∞. This model (for fixed population size N ) lets each individual in any generation r produce a number of (potential) offspring X (r ) i , i.i.d. across individuals and generations, distributed as a tail-heavy random variable X with E(X ) = μ > 1, i.e.
where C > 0 is a constant and 1 < α < 2. Then, N offspring are chosen to form the next generation. If less than N offspring are produced, the missing next generation individuals are arbitrarily associated with parents. Here, this is done by randomly choosing a parent, which preserves exchangeability and makes the model a Cannings model. The genealogies of a sample of size n converge as N → ∞ and time rescaled by c −1 N to the Beta(2 − α,α)-n-coalescent, see Schweinsberg (2003, Thm. 4). This model can very easily extended to variable population sizes by just sampling from the potential offspring. The tail-heavy distributions used produce, asymptotically as N → ∞, enough potential offspring to cover growing population sizes of order N as allowed by Eq. (4).
This gives us an alternative Cannings model with variable population sizes to define time-changed Beta coalescents as the limit of its discrete genealogies.

Corollary 3 For a population size profiles of exponential growth (on the coalescent scale) with growth rate ρ and for c N
This implies that the waiting time between coalescent events are Gompertz distributed with parameters a = λ b e ργ t 0 and b = ργ , i.e. the waiting time T for the next coalescence event, given the last coalescence at t 0 into b lineages, fulfills Remark 2 It is well-known that for Kingman's n-coalescent with exponential growth, waiting times for coalescence events follow a Gompertz distribution, e.g. see Slatkin and Hudson (1991, Eq. 5), Polanski et al. (2003). For time-changed Dirac coalescents appearing as limits of modified Moran models with ν(t) = ex p(−ρt), Eq. (10) appeared in Matuszewski et al. (2017).

Discussion
As for the Wright-Fisher model, genealogies of samples taken from (haploid) modified Moran and other Cannings models can be approximated by a time-change of their limit coalescent process, when the population sizes of the discrete models are fluctuating, but are always of the same order of size. As for models with fixed population size, time intervals of [c −1 N t] generations in the discrete model correspond to a time interval of length t in the continuous time limit. The approach of this study was to build on existing Cannings models that converge for fixed population size to the Λ-n-coalescent and just change the population sizes gradually from generation to generation, which includes adjusting parent-offspring allocation between generations. This raises the question whether the used Cannings models and the adjustment of ancestral relationships have biological interpretations and are a reasonable model for at least some real populations.

Interpretation of the Cannings models and allocation schemes used
The modified Moran models used to construct a time-changed Λ-n-coalescent with Λ([0, 1]) = 1 [defined via Eq. (6), introduced in Huillet and Möhle (2013)] can be described as follows (for fixed N ): On top of a standard Moran model choice of one parent M with two offspring and one individual in the parent generation with no offpring, there is a random probability X for each other individual in the parent generation to not have offspring in the next generation. X is drawn from Λ, potentially only activated in a given generation with a low probability N −γ , γ ∈ (1, 2). From the individuals that have offspring, all but M reproduce once, and M replaces itself and all non-reproducing individuals by its offspring. These models capture the concept of sweepstake reproduction (Hedgecock and Pudovkin 2011), though the assumption of a single individual with more than one offspring is rather artificial. For a non-random X and large families appearing occasionally at rate of order N −γ , this model is very similar to the discrete modified Moran model from Eldon and Wakeley (2006) used to describe sweepstake reproduction (and that was used in Matuszewski et al. (2017) as a basis to construct a time-changed Dirac n-coalescent). Both models lead to the same Dirac coalescent limit and have the same time rescaling order c −1 N . In Eldon and Wakeley (2006), instead of randomly choosing individuals to not reproduce with probability X , a fixed number of ≈ N X − 2 individuals are chosen at random to not reproduce on top of the Moran choice (again with a small probability in each generation for this to happen). For random X , similar models also appear in Hartmann and Huillet (2018) and Eldon (2012). The other class of Cannings models used to capture skewed offspring distributions, defined via Eq. (9), lead to the specific class of Beta(2−α, α)-n-coalescents. They have been proposed as a model of type-III survivorship, where all individuals produce many offspring with a high juvenile mortality, see e.g. Steinrücken et al. (2013, Sect. 2.3), also leading to sweepstake-like phenomena. While both classes of Cannings models allow the Bolthausen-Sznitman n-coalescent (Λ = Beta(1, 1)) as a possible limit model, the discrete models used to explicitly construct it are not based on modelling a directed selection process due to selective advantages of certain ancestral lineages. Thus, the results do not answer whether adding population size changes to a model of rapid selection or genetic draft as in Desai et al. (2013), Neher and Hallatschek (2013), Schweinsberg (2017) also leads to its rescaled genealogies being described by a time-changed Bolthausen-Sznitman n-coalescent. To construct time-changed Λ-n-coalescents as limits of genealogies in modified Moran models, the approach here is to adjust fixed-N modified Moran models for growing or decreasing population sizes. Sampling the next generation from the fixed-N offspring when there is population decline maintains on average the ratio between the large family U N and the rest off the individuals. This means that the population decrease, e.g. due to less resources available, has the same chance to affect each offspring of the fixedsize model. Additional individuals can be added to the family of the multiplying parent or by allowing parents with no offspring from the fixed-N allocation scheme to have exactly one offspring. For some sequences of modified Moran models, any partition of additional individuals to these two allocation forms is possible, e.g. allocate them randomly to the multiplying parent (with U N r offspring) from the fixed-size model with probability U N r /N r (with the constraint that we cannot add more than U N r − 1 individuals to non-reproducing parents). The merit of this random allocation is that it is trying to maintain the ratio U N r /N r from the fixed-size model. As for sampling a smaller number of individuals, this describes that population size increase, e.g. due to more resources available, follows (approximately and on average) the sweepstake pattern of the fixed-N model. From a biological viewpoint, other allocation schemes can also be interpreted: Adding the additional offspring completely to the largest family, as done in Matuszewski et al. (2017), could describe a scenario where new resources become available and only the multiple-offspring parent can claim them for its offspring. In contrast, adding individuals as single offspring of non-reproducing parents from the fixed-size model relaxes the (viability) "selection" pressure of the modified Moran model by allowing more non-multiplying parents (resp. their offspring) to survive, e.g. due to the additional resources. For the models covered in Theorem 3 from Schweinsberg (2003), population size changes in either direction are modelled by sampling from a pool of more individuals than the current population size, thus additional or decreasing resources affect the offspring of different parents in the same way.

Influence of the choice of Cannings model on the limit
Many results in the present paper allow us to scale the time in the discrete models with c −1 N as in the fixed N case so that the scaled genealogies converge to a time-changed Λn-coalescent (Π G(t) ) t≥0 . This time-change G(t) depends both on the population size profile ν on the coalescent time scale from Eq. (4) and the (asymptotic properties of) the coalescence probabilities c N , i.e. how many discrete generation correspond to one unit of coalescent time. For instance, consider an exponentially growing population (on the coalescent time scale, ν(t) = exp(−ρt) for ρ > 0) and two different models leading to a time-changed Beta(2 − α, α)-n-coalescent (α ∈ (1, 2)): the ones from Corollary 1 and Theorem 3. From Eq. (10), we see that G depends on the product γρ. For the model from Corollary 2, γ = α and for the one from Proposition 3, it is γ = α − 1. Thus, the exact same time-changed Λ-n-coalescent can appear as limit model for genealogies with different population size profiles on the coalescent time scale. As already discussed in (Matuszewski et al. 2017) in the case of timechanged Dirac-n-coalescents, this poses a problem for inference: If one wants to infer ρ directly (instead of the compound parameter γρ), γ has to be known. This means that specifying/identifying the Cannings model leading to the limit process would be necessary to directly estimate ρ. This is very similar to the effect that e.g. Watterson's estimator only estimates the mutation rate on the coalescent time scale, and not the mutation rate in one generation, see e.g. Eldon andWakeley (2006, p. 2627). Another example for different ν leading to the same time-scaled coalescent limit for different Cannings models is given by the genealogy limit from the Wright-Fisher model and the (usual) Moran models. It is well known, see e.g. Griffiths and Tavare (1994), that the rescaled genealogy of a sample from a Wright-Fisher model with population size profile ν converges to Kingman's n-coalescent with time change G as in Eq. (22) with γ = 1. However, for the classic Moran model, Proposition 1 shows that Eq. (22) holds with γ = 2. For families of Cannings models, if the coalescence probability c N is of order log(N ) −1 , a curious phenomenon appears: Population size changes of order N do not even alter the limit genealogy. An example is the model from Proposition 3 for the Bolthausen-Sznitman n-coalescent (Λ = Beta(1, 1)). One can interpret this for a population described by the model as follows: Even instantaneous bottlenecks or expansions do not influence the effect that a very large family appearing in a generation has on the genealogy. How the population reproduces, i.e. how the offspring distributions compare between different parents, is thus fully controlling the genealogy, regardless of changes that alter the population sizes overall, e.g. changes in range and/or resources.

Proofs
This section contains the proof of the presented statements as well as some further remarks.

Converging to a time-changed coalescent: sufficient conditions
First, recall this special case of Möhle (2002, Thm. 2

bution (Skorohod-sense) to a continuous-time Markov chain with transition function exp(Qt)
, where Q is a transition rate matrix with entries q a 1 ,...,a l , a 1 ≥ · · · ≥ a l ≥ 2 (so diagonal entries are the negative row sums of the other entries).
Remark 3 When compared to the original formulation of Möhle (2002, Thm 2.2), the limit here can be described as a homogeneous Markov chain with rate matrix Q instead of the more complicated original description of the transition probabilities as a product integral of matrix-valued measures. This directly follows from the stronger condition (12) Möhle (2002, p. 209), see also Eq. (24) therein. Then, the form of the transition function is described on Möhle (2002, p. 203). Now, recall the conditions (4), (5). Additionally, consider the following control condition for the fluctuations of c N as N → ∞: for all r ≤ [tc −1 N ]. For instance, when (4) holds, this condition is satisfied if c N = f (N ), where f is regularly varying (at ∞). If Eqs. (5) and (13) hold, choosing l = 1 in Eq. (5) yields as N → ∞. Now, we establish an easy-to-verify variant of Möhle (2002, Corollary 2.4).

Lemma 2 Consider a sequence of Cannings models with reference size N
as N → ∞ Proof Size changes of order N satisfy the first part of Condition (11). Its second part is then satisfied by (14), which in turn is satisfied due to (5) and (13). Also due to (14), as N → ∞ and [s] ≤ c −1 N t and thus its pseudo-inverse by with an appropriate t ≥ t. This implies that the time change function G −1 N for the discrete models in Corollary 4 is of order c −1 N . Knowing this, we compute as N → ∞ The second equation is valid due to the uniform convergence of Φ (N r ) l (a 1 , . . . , a l )c −1 N r in r for N → ∞ (N r is bounded from below on the timescale used). This allows us to pull out φ l (a 1 , . . . , a l ). This shows that condition 12 is satisfied and thus establishes the convergence of (R ) t≥0 to the same Λ-n-coalescent as the fixed-size model. Eq. (15) follows as described in Möhle (1998, Sec. 4).

Remark 4
The condition for Eq. (15) to hold is a weak condition, since G −1 N (t) is of order c −1 N . Additionally, the linear scaling in (15) makes it easy to introduce a mutation structure. Let mutation be introduced in the discrete model by allowing mutations from parent to offspring with a rate μ N . If μ N c −1 N → θ as N → ∞, the mutations on the time-scaled Λ-n-coalescent are given by a Poisson point process with homogeneous intensity θ . The next step is to establish a special case of Lemma 2 which only considers modified Moran models with changing population sizes.
Remark 5 Depending on the magnitude of a population size increase, adding individuals as further offspring of the multiplying parent from the fixed-size modified Moran model can strongly increase coalescence probabilities. For instance, for a population expansion of size N m, if one just expands by adding d N ,r = N m to the offspring number of the individual with multiple offspring in a single generation, the coalescence probability for this generation is dominated by the population size change. Then U N ,r ≥ N m, leading to c N ,r Thus, from generation r − 1 to r , coalescence is still happening with positive probability as N → ∞, which shows that a potential limit coalescent cannot just be a (non-degenerately) time-changed Λ-n-coalescent, a continuous-time (inhomogeneous) Markovian process. This has an implication for modelling of real populations: The genealogy of a sudden population expansion, happening at a specific generation, where a single genotype/individual is responsible for the population growth, is not given by a time-changed continuous-time Λ-n-coalescent.
We recall some properties of fixed-N modified Moran models.
for all k ≥ 2.

coalescent limit for the fixed-N modified Moran model.
Proof This is shown by applying Lemma 2. All conditions but Eq. (5) of it are clearly fulfilled under the assumptions of the proposition currently proven, see also Lemma 3. To show (5), first assume d N ,r ≥ 0. Then, where s(n, k) are Stirling numbers of the first kind. From Eq. (17), we see that uniformly in r as N → ∞ (convergence of the first summand is at least as fast as for N − (t)), while we will now show that the sum following this term in Eq. (19) vanishes asymptotically. We will give an upper bound for E(U l Nr A k−l N ,r ) (N r ) a 1 c N , independent from r . For this, we need to recall several technical points: For an upper bound, we can always omit terms of the form which follows from the fact that x → x a 1 (x) a 1 decreases for x ≥ a 1 . With all this, we can observe that, for 0 ≤ r ≤ c −1 N t and 0 ≤ l < k < a 1 , as N → ∞. We can thus establish convergence, uniform in r , in Eq. (19), since we have just shown the uniform convergence of all its summands: as N → ∞.
Regardless of the allocation of the new individuals, the population model is a modified Moran model with a single multiplying parent. Thus, to show Eq. (5) one only needs to show Φ (N ) Equation ( * ) follows from Eq. (21) and, for the first factor, from N −1 d N ,r ≤ d N being a null sequence. Now, consider d N ,r < 0. Then, we get the offspring population by sampling N r −1 individuals out of N r , from which U N r share one common parent. Thus, this is again a modified Moran model, where U n,r is conditionally hypergeometrically distributed . Using the factorial moment of the hypergeometric distribution leads to
Finally, the following lemma provides sufficient conditions for shifting the timechange G −1 N from pre-limit to limit. Lemma 4 Assume that for a Cannings model with variable population sizes (N r ) r ∈N , the discrete n-coalescents satisfy (R (N ) (1).

Further assume that Eq. (4) is satisfied for a positive real function ν and that c N
since for a sequence of functions, the inverses converge iff the original functions converge and since c f has the inverse t → f −1 (c −1 t). The shift by -1 does not alter the limit here, since its effect vanishes as N → ∞ due to the multiplication with c N . It is important to note here that any terms of order o (c N ) can be omitted when computing F N . Thus, we can replace c N ,r by c N r and even by c * N r = f (N r ) for a constant c 2 . Since N −1 N tc −1 N → ν(t), analogous to Griffiths and Tavare (1994), we can show, for f (x) = cx −γ , where for convergence, observe that there is pointwise convergence inside the integral and that bounded convergence is applicable since Eq. (4) ensures that the integrand is in Remark 7 -The integral representation of the time change is a deterministic version of the coalescent intensity from Kaj and Krone (2003, Sect. 1.3), just applied to Cannings models leading to non-Kingman Λ-n-coalescents. -As described in Möhle (1998, Section 4), the time-changed Λ-n-coalescent limit can also be expressed by its infinitesimal rates for a merger of k of n present lineages. This is also the form in which the limit process of the diploid umbrella model from Koskela and Wilke Berenguer (2019) is given. -Conditioned that the limit coalescent (Π G(t) ) t≥0 has at time T 0 = t 0 coalesced into a state with b blocks, what is the distribution of the waiting time T for the next coalescence event? If T = t, this means that in the non-rescaled Λ-n-coalescent (Π t ) t≥0 , we wait G(t) − G(t 0 ) for the next coalescence. This waiting time T in (Π t ) t≥0 is exponentially distributed with parameter λ b (total rate of coalescence). Thus,

Proof of Corollary 3
The form of G is a direct consequence of Lemma 4, since ν(t) = ex p(−ρt). Then, plugging G into Eq. (23) yields the distribution for the next coales-cence event, the Gompertz distribution parameters as e.g. described in Lenart (2014, Eq. 3) can be read off.

Proofs of convergence to a time-changed coalescent: modified Moran models
The modified Moran models used in Theorems 1 and 2 were introduced in Huillet and Möhle (2013, Prop. 4), the latter model with a small modification to ensure that there is always a parent with at least two offspring, see also Huillet and Möhle (2013, Example 4.1).

Proof of Theorem 1
Assume that E((U N ) 2 )(N − 1) −1 0 as N → ∞ also holds for any subsequence. If not, restrict to a subsequence for which this is true and define the limit only along this subsequence. First, we verify that c N = N −γ + o (N −γ ), thus converges to 0 and that the fixed-N model converges to the Λ-n-coalescent. Let c N be the coalescence probability in a fixed-N modified Moran model with U N as the number of offspring of the multiplying parent. This ensures 2((N ) 2 ) −1 ≤ c N ≤ 1. Moreover, the assumptions made ensure that (N c N ) N ∈N has a lower bound > 0, so we can define A n s.t. P(A N )c N = N −γ for any γ ∈ (1, 2). Then, the following is satisfied as N → ∞ and X where Eq. ( * ) uses c N = λ −1 N , which follows from Eq. (18) with k = 2. This establishes the convergence to the Λ-n-coalescent in the fixed-N case. Now we assume variable population sizes (N r ) r ∈N . First, observe that, since c N = O(N −γ ), it is enough to add occasionally a single individual from generation r to r − 1 to generate any population size changes allowed in Eq. (4) including bottlenecks which are instantaneous on the coalescent time scale. This single individual can then be added as offspring of a non-multiplying parent from the fixed-N r model or as an offspring of the already multiplying parent. To see the latter, observe that E((U N ) 2 ) = c N (N ) 2 ∼ N −γ N 2 → ∞. Then, as in Huillet and Möhle (2013, third remark p. 8), one has .
Thus, Proposition 1 allows one to add the one individual also to the already multiplying parent.
Thus, we have verified all conditions but Eq. (13) to apply Proposition 2. However, this follows from c N regularly varying. Since c N = N −γ + o(N −γ ), we can also shift the non-linear time-change to the limit due to Lemma 4.
For the proof of the next theorem, we use the following Lemma 5 Let Λ ∈ M[0, 1] and let U N be distributed as in Eq. (6) for any N ∈ N. Let (N r ) r ∈N satisfy Eq. (4). Then, Eq. (13) is satisfied.
Proof Equation (18) shows c N = λ −1 N , where λ N is the total transition rate for the first jump of a Λ-N -coalescent. Without restriction, assume N r ≥ N (for N < N r , an analogous proof provides bounds for c N c Nr ). Further assume that the N -coalescent is just the restriction of the N r -coalescent on individuals {1, . . . , N }. Any merger in the N -coalescent is then also a merger in the N r -coalescent, which shows λ N ≤ λ N r . In contrast, the first merger in the N r -coalescent is only a merger in the N -coalescent if it features at least two individuals from {1, . . . , N }. The probability of this is bounded from below by the probability N (N −1) N r (N r −1) that the first two of the blocks merged in the N are from {1, . . . , N }. This implies Proof of Theorem 2 In the fixed-N case, Eq. (6) implies that E((U N ) 2 )(N −1) −1 → 0 necessarily needs that λ N → ∞ as N → ∞. This is equivalent to 1 0 x −2 Λ(dx) = ∞, see (Pitman 1999, Eq. 7). Thus, convergence to the Λ-n-coalescent is shown in (Huillet and Möhle 2013, Prop. 3.4). Now, switch to variable population sizes (N r ) r ∈N . Since N c N = E((U N ) 2 )(N − 1) −1 → 0 as N → ∞, it is enough to add one individual per generation to cover any population growth profile covered by Eq. (4). This can always be done by letting a parent not reproducing in the fixed-size model reproduce (once). To add as further offspring of the multiplying parent, assume E((U N ) 2 ) → ∞ as N → ∞. Then, E(U N ) E((U N ) 2 ) → 0 as shown in the proof of Theorem 1. Thus, Proposition 1 provides that at most A N ,r ≤ c 4 E(U N ) for arbitrary c 4 > 0 individuals can be added per generation to the already multiplying parent. This allows for adding up to any fixed number k ∈ N individuals per generation. Additionally, Eq. (13) is satisfied due to Lemma 5. We can thus apply Proposition 2, with an arbitrary allocation of additional individuals that yields a modified Moran model.
To prove the Corollaries 1 and 2, we collect some properties of the modified Moran models with U N = U N with U N given by Eq. (6) leading to Beta-(a, b)-coalescents for a ∈ (0, 2], b > 0. From Huillet and Möhle (2013, Eq. (10) + Corollary A.1), Proof of Corollary 1 From Eq. (25), it follows that for a ∈ (0, 1), Thus, such Λ-n-coalescents are covered by Theorem 2. Additionally from Eq. (25), c N has a form that is covered by Lemma 4, which allows us to shift the time-change G in Theorem 2 to the limit coalescent and also shows the form of G in Corollary 1.

Proof of Corollary 2
We reiterate the proof of Theorem 2. Let Λ = Beta(a, b) for a ∈ (1, 2), which satisfies x −2 Λ(dx) = ∞. Thus, in the fixed-N case, again (Huillet and Möhle 2013, Prop. 3.4) ensures the convergence of the discrete genealogies to the Λ-n-coalescent when properly rescaled as N → ∞. Furthermore, Lemma 5 shows that Eq. (13) is satisfied Since ν(t) = ex p(−ρt), we can use N r = N (1 − ρc N ) to satisfy (4). Thus, we only need to show that the population size increase per generation does not violate the conditions of Proposition 2. Indeed, individuals at most have to be added. These can be added as A n,r additional offspring of the multiplying parent from the fixed-N r model, if the condition to apply it from Lemma 2 are met. For Λ considered here, one has E((U N ) 2 )) ∼ N 2 c N ∼ N a → ∞, see Eq. (25). From Lemma 2 we see that then we are allowed to add O(E(U N )) individuals. Huillet and Möhle (2013, Remark p. 9) shows E(U N ) = c 5 N c N for a constant c 5 > 0, so such growth is indeed covered (and we can then still use A N ,r < d N r and add the other individuals to non-reproducing parents from the fixed-N r model). Thus, we can establish convergence using Proposition 2 and shift the time-change G to the limit using Lemma 4, since c N is essentially a negative power of N .

Proofs of converging to a time-changed coalescent: model from Schweinsberg (2003)
Proof of Lemma 1 It suffices to reiterate the proof of Schweinsberg (2003, Lemma 5) briefly. For u ∈ [0, 1], consider the generating function f (u) := E(u X ). Let d n,r = d N ,r /N r . Then, S N ,r := N r i=1 X (r ) i fulfills P(S N ,r ≤ N r −1 ) ≤ u −N r (1+d N ,r ) E(u S N ,r ) = (u −(1+d N ,r ) f (u)) N r .

Proof of Theorem 3
Recall that, for 1 < α < 2, the coalescence probability in this (fixed-N ) model satisfies c N ∼ Cα B(2 − α, α)E(X ) −α N 1−α , where B is the Beta function, see Schweinsberg (2003, Lemma 13). For α = 1, instead c N ∼ (log N ) −1 , see Schweinsberg (2003, Lemma 16). Check the conditions necessary to apply Lemma 2: the model and the assumptions above satisfy c N → 0 as N → ∞, (4) and, since c N is regularly varying, also (13). The changes of population sizes from generation to generation are enough to cover instantanous population size changes on the coalescent time scale (and these are the most extreme changes allowed in Eq. (4)): for a (coalescent time) instantaneous change of size m N , one can set |d N ,r | = m √ c N =: d N → 0 as N → ∞ for ( √ c N ) −1 generations. Thus, only (5) needs to be verified. Since N r −1 offspring are sampled from N r i=1 X (r ) i potential offspring, the transition probabilities of the discrete coalescent can be formulated analogously to Eq. (2) as Φ (N ) (r ; a 1 , . . . , a l ) = This means one just needs to show that Φ (N r ) (a 1 , . . . , a l ) which follows from c −1 N Φ (N r ) (a 1 , . . . , a l ) To show the latter, uniform convergence, proceed as follows. First, observe that, since (13) holds, c −1 N Φ (N r ) (a 1 , . . . , a l ) = c −1 N r Φ (N r ) (a 1 , . . . , a l ) →φ(a 1 ,...,a l ) c N r c N is uniformly bounded (again, since N r is bounded from below by N − (t), there is uniform convergence in r of the first factor as N → ∞). Thus, we only need to show − 1 → 0. For this, observe that the function x → a 1 −x a 2 −x for x < a 2 is strictly increasing (decreasing) if a 1 − a 2 > 0 (if a 1 − a 2 < 0). This shows that there are b 1 , b 2 ∈ N 0 so that This implies that it is sufficient to show Since d N → 0, this vanishes uniformly. Thus, Lemma 2 can be applied, establishing convergence of (R . Lemma 4 then ensures that the timechange G can be shifted to the limit, since c N is either essentially a negative power or a logarithm of N .