Theoretical Population Biology最新文献

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.

In this paper we study invasion probabilities and invasion times of cooperative parasites spreading in spatially structured host populations. The spatial structure of the host population is given by a random geometric graph on ${[0,1]}^n$, $n\in N$, with a Poisson$\left(N\right)$-distributed number of vertices and in which vertices are connected over an edge when they have a distance of at most $r_N$ with $r_N$ of order $N^{(\beta -1)/n}$ for some $0<\beta <1$. At a host infection many parasites are generated and parasites move along edges to neighbouring hosts. We assume that parasites have to cooperate to infect hosts, in the sense that at least two parasites need to attack a host simultaneously. We find lower and upper bounds on the invasion probability of the parasites in terms of survival probabilities of branching processes with cooperation. Furthermore, we characterize the asymptotic invasion time.

An important ingredient of the proofs is a comparison with infection dynamics of cooperative parasites in host populations structured according to a complete graph, i.e. in well-mixed host populations. For these infection processes we can show that invasion probabilities are asymptotically equal to survival probabilities of branching processes with cooperation. Furthermore, we build on proof techniques developed in Brouard and Pokalyuk (2022), where an analogous invasion process has been studied for host populations structured according to a configuration model.

We substantiate our results with simulations.

Understanding the conditions that promote the evolution of migration is important in ecology and evolution. When environments are fixed and there is one most favorable site, migration to other sites lowers overall growth rate and is not favored. Here we ask, can environmental variability favor migration when there is one best site on average? Previous work suggests that the answer is yes, but a general and precise answer remained elusive. Here we establish new, rigorous inequalities to show (and use simulations to illustrate) how stochastic growth rate can increase with migration when fitness (dis)advantages fluctuate over time across sites. The effect of migration between sites on the overall stochastic growth rate depends on the difference in expected growth rates and the variance of the fluctuating difference in growth rates. When fluctuations (variance) are large, a population can benefit from bursts of higher growth in sites that are worse on average. Such bursts become more probable as the between-site variance increases. Our results apply to many ($\ge$ 2) sites, and reveal an interplay between the length of paths between sites, the average differences in site-specific growth rates, and the size of fluctuations. Our findings have implications for evolutionary biology as they provide conditions for departure from the reduction principle, and for ecological dynamics: even when there are superior sites in a sea of poor habitats, variability and habitat quality across space determine the importance of migration.

The host microbiome can be considered an ecological community of microbes present inside a complex and dynamic host environment. The host is under selective pressure to ensure that its microbiome remains beneficial. The host can impose a range of ecological filters including the immune response that can influence the assembly and composition of the microbial community. How the host immune response interacts with the within-microbiome community dynamics to affect the assembly of the microbiome has been largely unexplored. We present here a mathematical framework to elucidate the role of host immune response and its interaction with the balance of ecological interactions types within the microbiome community. We find that highly mutualistic microbial communities characteristic of high community density are most susceptible to changes in immune control and become invasion prone as host immune control strength is increased. Whereas highly competitive communities remain relatively stable in resisting invasion to changing host immune control. Our model reveals that the host immune control can interact in unexpected ways with a microbial community depending on the prevalent ecological interactions types for that community. We stress the need to incorporate the role of host-control mechanisms to better understand microbiome community assembly and stability.

Sexual selection plays a crucial role in modern evolutionary theory, offering valuable insight into evolutionary patterns and species diversity. Recently, a comprehensive definition of sexual selection has been proposed, defining it as any selection that arises from fitness differences associated with nonrandom success in the competition for access to gametes for fertilization. Previous research on discrete traits demonstrated that non-random mating can be effectively quantified using Jeffreys (or symmetrized Kullback-Leibler) divergence, capturing information acquired through mating influenced by mutual mating propensities instead of random occurrences. This novel theoretical framework allows for detecting and assessing the strength of sexual selection and assortative mating.

In this study, we aim to achieve two primary objectives. Firstly, we demonstrate the seamless alignment of the previous theoretical development, rooted in information theory and mutual mating propensity, with the aforementioned definition of sexual selection. Secondly, we extend the theory to encompass quantitative traits. Our findings reveal that sexual selection and assortative mating can be quantified effectively for quantitative traits by measuring the information gain relative to the random mating pattern. The connection of the information indices of sexual selection with the classical measures of sexual selection is established.

Additionally, if mating traits are normally distributed, the measure capturing the underlying information of assortative mating is a function of the square of the correlation coefficient, taking values within the non-negative real number set [0, +∞).

It is worth noting that the same divergence measure captures information acquired through mating for both discrete and quantitative traits. This is interesting as it provides a common context and can help simplify the study of sexual selection patterns.

The coalescent is a stochastic process representing ancestral lineages in a population undergoing neutral genetic drift. Originally defined for a well-mixed population, the coalescent has been adapted in various ways to accommodate spatial, age, and class structure, along with other features of real-world populations. To further extend the range of population structures to which coalescent theory applies, we formulate a coalescent process for a broad class of neutral drift models with arbitrary – but fixed – spatial, age, sex, and class structure, haploid or diploid genetics, and any fixed mating pattern. Here, the coalescent is represented as a random sequence of mappings $C={\left(C_t\right)}_{t=0}^{\infty }$ from a finite set $G$ to itself. The set $G$ represents the “sites” (in individuals, in particular locations and/or classes) at which these alleles can live. The state of the coalescent, $C_t:G\to G$, maps each site $g\in G$ to the site containing $g$’s ancestor, $t$ time-steps into the past. Using this representation, we define and analyze coalescence time, coalescence branch length, mutations prior to coalescence, and stationary probabilities of identity-by-descent and identity-by-state. For low mutation, we provide a recipe for computing identity-by-descent and identity-by-state probabilities via the coalescent. Applying our results to a diploid population with arbitrary sex ratio $r$, we find that measures of genetic dissimilarity, among any set of sites, are scaled by $4r(1-r)$ relative to the even sex ratio case.

The introduction of the spatial Lambda-Fleming–Viot model ($\Lambda$V) in population genetics was mainly driven by the pioneering work of Alison Etheridge, in collaboration with Nick Barton and Amandine Véber about ten years ago (Barton et al., 2010; Barton et al., 2013). The $\Lambda$V model provides a sound mathematical framework for describing the evolution of a population of related individuals along a spatial continuum. It alleviates the “pain in the torus” issue with Wright and Malécot’s isolation by distance model and is sampling consistent, making it a tool of choice for statistical inference. Yet, little is known about the potential connections between the $\Lambda$V and other stochastic processes generating trees and the spatial coordinates along the corresponding lineages. This work focuses on a version of the $\Lambda$V whereby lineages move rapidly over small distances. Using simulations, we show that the induced $\Lambda$V tree-generating process is well approximated by a birth–death model. Our results also indicate that Brownian motions modelling the movements of lines of descent along birth–death trees do not generally provide a good approximation of the $\Lambda$V due to habitat boundaries effects that play an increasingly important role in the long run. Accounting for habitat boundaries through reflected Brownian motions considerably increases the similarity to the $\Lambda$V model however. Finally, we describe efficient algorithms for fast simulation of the backward and forward in time versions of the $\Lambda$V model.

Muller’s ratchet, in its prototype version, models a haploid, asexual population whose size $N$ is constant over the generations. Slightly deleterious mutations are acquired along the lineages at a constant rate, and individuals carrying less mutations have a selective advantage. The classical variant considers fitness proportional selection, but other fitness schemes are conceivable as well. Inspired by the work of Etheridge et al. (2009) we propose a parameter scaling which fits well to the “near-critical” regime that was in the focus of Etheridge et al. (2009) (and in which the mutation–selection ratio diverges logarithmically as $N\to \infty$). Using a Moran model, we investigate the“rule of thumb” given in Etheridge et al. (2009) for the click rate of the “classical ratchet” by putting it into the context of new results on the long-time evolution of the size of the best class of the ratchet with (binary) tournament selection. This variant of Muller’s ratchet was introduced in González Casanova et al. (2023), and was analysed there in a subcritical parameter regime. Other than that of the classical ratchet, the size of the best class of the tournament ratchet follows an autonomous dynamics up to the time of its extinction. It turns out that, under a suitable correspondence of the model parameters, this dynamics coincides with the so called Poisson profile approximation of the dynamics of the best class of the classical ratchet.

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.