Theoretical Population Biology最新文献

Social behavior is divided into four types: altruism, spite, mutualism, and selfishness. The former two are costly to the actor; therefore, from the perspective of natural selection, their existence can be regarded as mysterious. One potential setup which encourages the evolution of altruism and spite is repeated interaction. Players can behave conditionally based on their opponent's previous actions in the repeated interaction. On the one hand, the retaliatory strategy (who behaves altruistically when their opponent behaved altruistically and behaves non-altruistically when the opponent player behaved non-altruistically) is likely to evolve when players choose altruistic or selfish behavior in each round. On the other hand, the anti-retaliatory strategy (who is spiteful when the opponent was not spiteful and is not spiteful when the opponent player was spiteful) is likely to evolve when players opt for spiteful or mutualistic behavior in each round. These successful conditional behaviors can be favored by natural selection. Here, we notice that information on opponent players’ actions is not always available. When there is no such information, players cannot determine their behavior according to their opponent's action. By investigating the case of altruism, a previous study (Kurokawa, 2017, Mathematical Biosciences, 286, 94–103) found that persistent altruistic strategies, which choose the same action as the own previous action, are favored by natural selection. How, then, should a spiteful conditional strategy behave when the player does not know what their opponent did? By studying the repeated game, we find that persistent spiteful strategies, which choose the same action as the own previous action, are favored by natural selection. Altruism and spite differ concerning whether retaliatory or anti-retaliatory strategies are favored by natural selection; however, they are identical concerning whether persistent strategies are favored by natural selection.

In Athreya et al. (2021), models from population genetics were used to define stochastic dynamics in the space of graphons arising as continuum limits of dense graphs. In the present paper we exhibit an example of a simple neutral population genetics model for which this dynamics is a Markovian diffusion that can be characterized as the solution of a martingale problem. In particular, we consider a Markov chain in the space of finite graphs that resembles a Moran model with resampling and mutation. We encode the finite graphs as graphemes, which can be represented as a triple consisting of a vertex set (or more generally, a topological space), an adjacency matrix, and a sampling (Borel) measure. We equip the space of graphons with convergence of sample subgraph densities and show that the grapheme-valued Markov chain converges to a grapheme-valued diffusion as the number of vertices goes to infinity. We show that the grapheme-valued diffusion has a stationary distribution that is linked to the Griffiths–Engen–McCloskey (GEM) distribution. In a companion paper (Greven et al. 2023), we build up a general theory for obtaining grapheme-valued diffusions via genealogies of models in population genetics.

We consider a single genetic locus with two alleles $A_1$ and $A_2$ in a large haploid population. The locus is subject to selection and two-way, or recurrent, mutation. Assuming the allele frequencies follow a Wright–Fisher diffusion and have reached stationarity, we describe the asymptotic behaviors of the conditional gene genealogy and the latent mutations of a sample with known allele counts, when the count $n_1$ of allele $A_1$ is fixed, and when either or both the sample size $n$ and the selection strength $\left|\alpha \right|$ tend to infinity. Our study extends previous work under neutrality to the case of non-neutral rare alleles, asserting that when selection is not too strong relative to the sample size, even if it is strongly positive or strongly negative in the usual sense ($\alpha \to -\infty$ or $\alpha \to +\infty$), the number of latent mutations of the $n_1$ copies of allele $A_1$ follows the same distribution as the number of alleles in the Ewens sampling formula. On the other hand, very strong positive selection relative to the sample size leads to neutral gene genealogies with a single ancient latent mutation. We also demonstrate robustness of our asymptotic results against changing population sizes, when one of $\left|\alpha \right|$ or $n$ is large.

We study the response of a quantitative trait to exponential directional selection in a finite haploid population, both at the genetic and the phenotypic level. We assume an infinite sites model, in which the number of new mutations per generation in the population follows a Poisson distribution (with mean $\Theta$) and each mutation occurs at a new, previously monomorphic site. Mutation effects are beneficial and drawn from a distribution. Sites are unlinked and contribute additively to the trait. Assuming that selection is stronger than random genetic drift, we model the initial phase of the dynamics by a supercritical Galton–Watson process. This enables us to obtain time-dependent results. We show that the copy-number distribution of the mutant in generation $n$, conditioned on non-extinction until $n$, is described accurately by the deterministic increase from an initial distribution with mean 1. This distribution is related to the absolutely continuous part $W^{+}$ of the random variable, typically denoted $W$, that characterizes the stochasticity accumulating during the mutant’s sweep. A suitable transformation yields the approximate dynamics of the mutant frequency distribution in a Wright–Fisher population of size $N$. Our expression provides a very accurate approximation except when mutant frequencies are close to 1. On this basis, we derive explicitly the (approximate) time dependence of the expected mean and variance of the trait and of the expected number of segregating sites. Unexpectedly, we obtain highly accurate approximations for all times, even for the quasi-stationary phase when the expected per-generation response and the trait variance have equilibrated. The latter refine classical results. In addition, we find that $\Theta$ is the main determinant of the pattern of adaptation at the genetic level, i.e., whether the initial allele-frequency dynamics are best described by sweep-like patterns at few loci or small allele-frequency shifts at many. The number of segregating sites is an appropriate indicator for these patterns. The selection strength determines primarily the rate of adaptation. The accuracy of our results is tested by comprehensive simulations in a Wright–Fisher framework. We argue that our results apply to more complex forms of directional selection.

We consider the Moran model of population genetics with two types, mutation, and selection, and investigate the line of descent of a randomly-sampled individual from a contemporary population. We trace this ancestral line back into the distant past, far beyond the most recent common ancestor of the population (thus connecting population genetics to phylogeny), and analyse the mutation process along this line.

To this end, we use the pruned lookdown ancestral selection graph (Lenz et al., 2015), which consists of a set of potential ancestors of the sampled individual at any given time. Relative to the neutral case (that is, without selection), we obtain a general bias towards the beneficial type, an increase in the beneficial mutation rate, and a decrease in the deleterious mutation rate. This sheds new light on previous analytical results. We discuss our findings in the light of a well-known observation at the interface of phylogeny and population genetics, namely, the difference in the mutation rates (or, more precisely, mutation fluxes) estimated via phylogenetic methods relative to those observed in pedigree studies.

Infectious disease agents can influence each other’s dynamics in shared host populations. We consider such influence for two mosquito-borne infections where one pathogen is endemic at the time that a second pathogen invades. We regard a setting where the vector has a bias towards biting host individuals infected with the endemic pathogen and where there is a cost to co-infected hosts. As a motivating case study, we regard Plasmodium spp., that cause avian malaria, as the endemic pathogen, and Usutu virus (USUV) as the invading pathogen. Hosts with malaria attract more mosquitoes compared to susceptible hosts, a phenomenon named vector bias. The possible trade-off between the vector-bias effect and the co-infection mortality is studied using a compartmental epidemic model. We focus first on the basic reproduction number $R_0$ for Usutu virus invading into a malaria-endemic population, and then explore the long-term dynamics of both pathogens once Usutu virus has become established. We find that the vector bias facilitates the introduction of malaria into a susceptible population, as well as the introduction of Usutu in a malaria-endemic population. In the long term, however, both a vector bias and co-infection mortality lead to a decrease in the number of individuals infected with either pathogen, suggesting that avian malaria is unlikely to be a promoter of Usutu invasion. This proposed approach is general and allows for new insights into other negative associations between endemic and invading vector-borne pathogens.

We approach the questions, what part of evolutionary change results from selection, and what is the adaptive information flow into a population undergoing selection, as a problem of quantifying the divergence of typical trajectories realized under selection from the expected dynamics of their counterparts under a null stochastic-process model representing the absence of selection. This approach starts with a formulation of adaptation in terms of information and from that identifies selection from the genetic parameters that generate information flow; it is the reverse of a historical approach that defines selection in terms of fitness, and then identifies adaptive characters as those amplified in relative frequency by fitness. Adaptive information is a relative entropy on distributions of histories computed directly from the generators of stochastic evolutionary population processes, which in large population limits can be approximated by its leading exponential dependence as a large-deviation function. We study a particular class of generators that represent the genetic dependence of explicit transitions around reproductive cycles in terms of stoichiometry, familiar from chemical reaction networks. Following Smith (2023), which showed that partitioning evolutionary events among genetically distinct realizations of lifecycles yields a more consistent causal analysis through the Price equation than the construction from units of selection and fitness, here we show that it likewise yields more complete evolutionary information measures.

In this article, discrete and stochastic changes in (effective) population size are incorporated into the spectral representation of a biallelic diffusion process for drift and small mutation rates. A forward algorithm inspired by Hidden-Markov-Model (HMM) literature is used to compute exact sample allele frequency spectra for three demographic scenarios: single changes in (effective) population size, boom-bust dynamics, and stochastic fluctuations in (effective) population size. An approach for fully agnostic demographic inference from these sample allele spectra is explored, and sufficient statistics for stepwise changes in population size are found. Further, convergence behaviours of the polymorphic sample spectra for population size changes on different time scales are examined and discussed within the context of inference of the effective population size. Joint visual assessment of the sample spectra and the temporal coefficients of the spectral decomposition of the forward diffusion process is found to be important in determining departure from equilibrium. Stochastic changes in (effective) population size are shown to shape sample spectra particularly strongly.