No Arabic abstract
Consider a two-type Moran population of size $N$ subject to selection and mutation, which is immersed in a varying environment. The population is susceptible to exceptional changes in the environment, which accentuate the selective advantage of the fit individuals. In this setting, we show that the type-composition in the population is continuous with respect to the environment. This allows us to replace the deterministic environment by a random one, which is driven by a subordinator. Assuming that selection, mutation and the environment are weak in relation to $N$, we show that the type-frequency process, with time speed up by $N$, converges as $Ntoinfty$ to a Wright--Fisher-type SDE with a jump term modeling the effect of the environment. Next, we study the asymptotic behavior of the limiting model in the far future and in the distant past, both in the annealed and in the quenched setting. Our approach builds on the genealogical picture behind the model. The latter is described by means of an extension of the ancestral selection graph (ASG). The formal relation between forward and backward objects is given in the form of a moment duality between the type-frequency process and the line-counting process of a pruned version of the ASG. This relation yields characterizations of the annealed and the quenched moments of the asymptotic type distribution. A more involved pruning of the ASG allows us to obtain annealed and quenched results for the ancestral type distribution. In the absence of mutations, one of the types fixates and our results yield expressions for the fixation probabilities.
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 $Lambda$. 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.
We study ancestral structures for the two-type Moran model with two-way mutation and frequency-dependent selection that follows the nonlinear dominance or fittest-type-wins scheme. Both schemes lead, in distribution, to the same type-frequency process. Reasoning through the mutation structure on the ancestral selection graph (ASG), we derive processes suitable to determine the type distribution of the present and ancestral population, leading to, respectively, the killed and pruned lookdown ASG. To this end, we establish factorial moment dualities to the Moran model and a relative thereof, respectively. Finally, we extend the results to the diffusion limit.
The stationary distribution of the diffusion limit of the 2-island, 2-allele Wright-Fisher with small but otherwise arbitrary mutation and migration rates is investigated. Following a method developed by Burden and Tang (2016, 2017) for approximating the forward Kolmogorov equation, the stationary distribution is obtained to leading order as a set of line densities on the edges of the sample space, corresponding to states for which one island is bi-allelic and the other island is non-segregating, and a set of point masses at the corners of the sample space, corresponding to states for which both islands are simultaneously non-segregating. Analytic results for the corner probabilities and line densities are verified independently using the backward generator and for the corner probabilities using the coalescent.
The stationary distribution of a sample taken from a Wright-Fisher diffusion with general small mutation rates is found using a coalescent approach. The approximation is equivalent to having at most one mutation in the coalescent tree to the first order in the rates. The sample probabilities characterize an approximation for the stationary distribution from the Wright-Fisher diffusion. The approach is different from Burden and Tang (2016,2017) who use a probability flux argument to obtain the same results from a forward diffusion generator equation. The solution has interest because the solution is not known when rates are not small. An analogous solution is found for the configuration of alleles in a general exchangeable binary coalescent tree. In particular an explicit solution is found for a pure birth process tree when individuals reproduce at rate lambda.
The transition distribution of a sample taken from a Wright-Fisher diffusion with general small mutation rates is found using a coalescent approach. The approximation is equivalent to having at most one mutation in the coalescent tree of the sample up to the most recent common ancestor with additional mutations occurring on the lineage from the most recent common ancestor to the time origin if complete coalescence occurs before the origin. The sampling distribution leads to an approximation for the transition density in the diffusion with small mutation rates. This new solution has interest because the transition density in a Wright-Fisher diffusion with general mutation rates is not known.