Theoretical Population Biology最新文献

Distributional properties of tree shape statistics under random phylogenetic tree models play an important role in investigating the evolutionary forces underlying the observed phylogenies. In this paper, we study two subtree counting statistics, the number of cherries and that of pitchforks for the Ford model, the alpha model introduced by Daniel Ford. It is a one-parameter family of random phylogenetic tree models which includes the proportional to distinguishable arrangement (PDA) and the Yule models, two tree models commonly used in phylogenetics. Based on a non-uniform version of the extended Pólya urn models in which negative entries are permitted for their replacement matrices, we obtain the strong law of large numbers and the central limit theorem for the joint distribution of these two statistics for the Ford model. Furthermore, we derive a recursive formula for computing the exact joint distribution of these two statistics. This leads to exact formulas for their means and higher order asymptotic expansions of their second moments, which allows us to identify a critical parameter value for the correlation between these two statistics. That is, when the number of tree leaves is sufficiently large, they are negatively correlated for $0\le \alpha \le 1/2$ and positively correlated for $1/2<\alpha <1$.

Predicting temporal dynamics of genetic diversity is important for assessing long-term population persistence. In stage-structured populations, especially in perennial plant species, genetic diversity is often compared among life history stages, such as seedlings, juveniles, and flowerings, using neutral genetic markers. The comparison among stages is sometimes referred to as demographic genetic structure, which has been regarded as a proxy of potential genetic changes because individuals in mature stages will die and be replaced by those in more immature stages over the course of time. However, due to the lack of theoretical examination, the basic property of the stage-wise genetic diversity remained unclear. We developed a matrix model which was made up of difference equations of the probability of non-identical-by-descent of each life history stage at a neutral locus to describe the dynamics and the inter-stage differences of genetic diversity in stage-structured plant populations. Based on the model, we formulated demographic genetic structure as well as the annual change rate of the probability of non-identical-by-descent (denoted as $\eta$). We checked if theoretical expectations on demographic genetic structure and $\eta$ obtained from our model agreed with computational results of stochastic simulation using randomly generated 3,000 life histories. We then examined the relationships of demographic genetic structure with effective population size $N_e$, which is the determinants of diversity loss per generation time. Theoretical expectations on $\eta$ and demographic genetic structure fitted well to the results of stochastic simulation, supporting the validity of our model. Demographic genetic structure varied independently of $N_e$ and $\eta$, while having a strong correlation with stable stage distribution: genetic diversity was lower in stages with fewer individuals. Our results indicate that demographic genetic structure strongly reflects stable stage distribution, rather than temporal genetic dynamics, and that inferring future genetic diversity solely from demographic genetic structure would be misleading. Instead of demographic genetic structure, we propose $\eta$ as an useful tool to predict genetic diversity at the same time scale as population dynamics (i.e., per year), facilitating evaluation on population viability from a genetic point of view.

The concept of individual admixture (IA) assumes that the genome of individuals is composed of alleles inherited from $K$ ancestral populations. Each copy of each allele has the same chance $q_k$ to originate from population $k$, and together with the allele frequencies $p$ in all populations at all $M$ markers, comprises the admixture model. Here, we assume a supervised scheme, i.e. allele frequencies $p$ are given through a reference database of size $N$, and $q$ is estimated via maximum likelihood for a single sample. We study laws of large numbers and central limit theorems describing effects of finiteness of both, $M$ and $N$, on the estimate of $q$. We recall results for the effect of finite $M$, and provide a central limit theorem for the effect of finite $N$, introduce a new way to express the uncertainty in estimates in standard barplots, give simulation results, and discuss applications in forensic genetics.

Ancestral state reconstruction is one of the most important tasks in evolutionary biology. Conditions under which we can reliably reconstruct the ancestral state have been studied for both discrete and continuous traits. However, the connection between these results is unclear, and it seems that each model needs different conditions. In this work, we provide a unifying theory on the consistency of ancestral state reconstruction for various types of trait evolution models. Notably, we show that for a sequence of nested trees with bounded heights, the necessary and sufficient conditions for the existence of a consistent ancestral state reconstruction method under discrete models, the Brownian motion model, and the threshold model are equivalent. When tree heights are unbounded, we provide a simple counter-example to show that this equivalence is no longer valid.

This study deals with the problem of the population shrinking in habitats affected by aging and excessive migration outflows. First, a control-oriented population dynamics model was proposed that catches the effect of depopulation. The model also includes the effect of spatial interaction-driven migration flows on population size. The resulting model is a non-homogeneous ordinary differential equation. It includes such phenomena that are important from the control point of view, such as the influence of migration costs on population dynamics, the impact of aging on population size, or the effect of the habitats’ carrying capacity on migration flows. Based on the model, controllability conditions are formulated and a control strategy was developed that is meant to avoid the depopulation of the habitat. The control method acts on the migration costs to achieve the control goal and requires only population size measurements. Simulation measurements are presented in the paper to show the effectiveness of the proposed modeling and control methods.

Numerous traits under migration–selection balance are shown to exhibit complex patterns of genetic architecture with large variance in effect sizes. However, the conditions under which such genetic architectures are stable have yet to be investigated, because studying the influence of a large number of small allelic effects on the maintenance of spatial polymorphism is mathematically challenging, due to the high complexity of the systems that arise. In particular, in the most simple case of a haploid population in a two-patch environment, while it is known from population genetics that polymorphism at a single major-effect locus is stable in the symmetric case, there exist no analytical predictions on how this polymorphism holds when a polygenic background also contributes to the trait. Here we propose to answer this question by introducing a new eco-evo methodology that allows us to take into account the combined contributions of a major-effect locus and of a quantitative background resulting from small-effect loci, where inheritance is encoded according to an extension to the infinitesimal model. In a regime of small variance contributed by the quantitative loci, we justify that traits are concentrated around the major alleles, according to a normal distribution, using new convex analysis arguments. This allows a reduction in the complexity of the system using a separation of time scales approach. We predict an undocumented phenomenon of loss of polymorphism at the major-effect locus despite strong selection for local adaptation, because the quantitative background slowly disrupts the rapidly established polymorphism at the major-effect locus, which is confirmed by individual-based simulations. Our study highlights how segregation of a quantitative background can greatly impact the dynamics of major-effect loci by provoking migrational meltdowns. We also provide a comprehensive toolbox designed to describe how to apply our method to more complex population genetic models.

The concept of the average effect of an allele pervades much of evolutionary population genetics. In this context the average effect of an allele is often considered as the main component of the “fitness” of that allele. It is widely believed that, if this fitness component for an allele is positive, then the frequency of this allele will increase, at least for one generation in discrete-time models. In this note we show that this is not necessarily the case since the average effect of an allele on fitness may be different from its marginal additive fitness even in a one-locus setting in non-random-mating populations.

The Gini coefficient of the life table is a concentration index that provides information on lifespan variation. Originally proposed by economists to measure income and wealth inequalities, it has been widely used in population studies to investigate variation in ages at death. We focus on the complement of the Gini coefficient, Drewnowski’s index, which is a measure of equality. We study its mathematical properties and analyze how changes over time relate to changes in life expectancy. Further, we identify the threshold age below which mortality improvements are translated into decreasing lifespan variation and above which these improvements translate into increasing lifespan inequality. We illustrate our theoretical findings simulating scenarios of mortality improvement in the Gompertz model, and showing an example of application to Swedish life table data. Our experiments demonstrate how Drewnowski’s index can serve as an indicator of the shape of mortality patterns. These properties, along with our analytical findings, support studying lifespan variation alongside life expectancy trends in multiple species.

Principal component analysis (PCA) is one of the most frequently-used approach to describe population structure from multilocus genotype data. Regarding geographic range expansions of modern humans, interpretations of PCA have, however, been questioned, as there is uncertainty about the wave-like patterns that have been observed in principal components. It has indeed been argued that wave-like patterns are mathematical artifacts that arise generally when PCA is applied to data in which genetic differentiation increases with geographic distance. Here, we present an alternative theory for the observation of wave-like patterns in PCA. We study a coalescent model – the umbrella model – for the diffusion of genetic variants. The model is based on genetic drift without any particular geographical structure. In the umbrella model, splits from an ancestral population occur almost continuously in time, giving birth to small daughter populations at a regular pace. Our results provide detailed mathematical descriptions of eigenvalues and eigenvectors for the PCA of sampled genomic sequences under the model. When variants uniquely represented in the sample are removed, the PCA eigenvectors are defined as cosine functions of increasing periodicity, reproducing wave-like patterns observed in equilibrium isolation-by-distance models. Including singleton variants in the analysis, the eigenvectors corresponding to the largest eigenvalues exhibit complex wave shapes. The accuracy of our predictions is further investigated with coalescent simulations. Our analysis supports the hypothesis that highly structured wave-like patterns could arise from genetic drift only, and may not always be artificial outcomes of spatially structured data. Genomic data related to the peopling of the Americas are reanalyzed in the light of our new theory.

Chan, Durrett, and Lanchier introduced a multitype contact process with temporal heterogeneity involving two species competing for space on the d-dimensional integer lattice. Time is divided into two seasons. They proved that there is an open set of the parameters for which both species can coexist when their dispersal range is sufficiently large. Numerical simulations suggested that three species can coexist in the presence of two seasons. The main point of this paper is to prove that this conjecture is incorrect. To do this we prove results for a more general ODE model and contrast its behavior with other related systems that have been studied in order to understand the competitive exclusion principle.