Fixation in haploid populations exhibiting density dependence I: The non-neutral case

We extend the one-locus two allele Moran model of fixation in a haploid population to the case where the total size of the population is not fixed. The model is defined as a two-dimensional birth-and-death process for allele number. Changes in allele number occur through density-independent death events and birth events whose per capita rate decreases linearly with the total population density. Uniquely for models of this type, the latter is determined by these same birth-and-death events. This provides a framework for investigating both the effects of fluctuation in total population number through demographic stochasticity, and deterministic density-dependent changes in mean density, on allele fixation. We analyze this model using a combination of asymptotic analytic approximations supported by numerics. We find that for advantageous mutants demographic stochasticity of the resident population does not affect the fixation probability, but that deterministic changes in total density do. In contrast, for deleterious mutants, the fixation probability increases with increasing resident population fluctuation size, but is relatively insensitive to initial density. These phenomena cannot be described by simply using a harmonic mean effective population size.     

Absorption and fixation times for neutral and quasi-neutral populations with density dependence

We study a generalisation of Moran’s population-genetic model that incorporates density dependence. Rather than assuming fixed population size, we allow the number of individuals to vary stochastically with the same events that change allele number, according to a logistic growth process with density dependent mortality. We analyse the expected time to absorption and fixation in the ‘quasi-neutral’ case: both types have the same carrying capacity, achieved through a trade-off of birth and death rates. Such types would be competitively neutral in a classical, fixed-population Wright-Fisher model. Nonetheless, we find that absorption times are skewed compared to the Wright-Fisher model. The absorption time is longer than the Wright-Fisher prediction when the initial proportion of the type with higher birth rate is large, and shorter when it is small. By contrast, demographic stochasticity has no effect on the fixation or absorption times of truly neutral alleles in a large population. Our calculations provide the first analytic results on hitting times in a two-allele model, when the population size varies stochastically.     

Absorption and fixation times for neutral and quasi-neutral populations with density dependence

We study a generalisation of Moran’s population-genetic model that incorporates density dependence. Rather than assuming fixed population size, we allow the number of individuals to vary stochastically with the same events that change allele number, according to a logistic growth process with density dependent mortality. We analyse the expected time to absorption and fixation in the ‘quasi-neutral’ case: both types have the same carrying capacity, achieved through a trade-off of birth and death rates. Such types would be competitively neutral in a classical, fixed-population Wright–Fisher model. Nonetheless, we find that absorption times are skewed compared to the Wright–Fisher model. The absorption time is longer than the Wright–Fisher prediction when the initial proportion of the type with higher birth rate is large, and shorter when it is small. By contrast, demographic stochasticity has no effect on the fixation or absorption times of truly neutral alleles in a large population. Our calculations provide the first analytic results on hitting times in a two-allele model, when the population size varies stochastically.     

Fixation probabilities for the Moran process with three or more strategies: general and coupling results

We study fixation probabilities for the Moran stochastic process for the evolution of a population with three or more types of individuals and frequency-dependent fitnesses. Contrary to the case of populations with two types of individuals, in which fixation probabilities may be calculated by an exact formula, here we must solve a large system of linear equations. We first show that this system always has a unique solution. Other results are upper and lower bounds for the fixation probabilities obtained by coupling the Moran process with three strategies with birth–death processes with only two strategies. We also apply our bounds to the problem of evolution of cooperation in a population with three types of individuals already studied in a deterministic setting by Núñez Rodríguez and Neves (J Math Biol 73:1665–1690, 2016). We argue that cooperators will be fixated in the population with probability arbitrarily close to 1 for a large region of initial conditions and large enough population sizes.

     

How should we test for the role of behaviourin population dynamics?

The behaviour of an individual affects the probability that it will find food or a mate, and whether it will avoid becoming food itself. Because the birth and death of individuals are the constituents of the birth and death rates of populations, it seems likely that population dynamics is affected by variation in the behaviour of individuals. While intuitively appealing, there are few data to support this contention. I suggest that this lack of support stems from the failure to make use of manipulative experiments to test the hypothesis. Manipulative experiments are not the only approach to testing these hypotheses, but this powerful tool has not been used as effectively as it might. The behaviour of an animal can be manipulated with information-carrying chemicals in its environment, and pharmacologically. Variation in behaviour among individuals that is genetically based can also provide experimental material to test hypotheses about the role of behaviour in population dynamics. Manipulative experiments have the advantages of increased statistical power, and the elimination of unmeasured covariates. They have the disadvantage that they can introduce artifacts into the study system. It seems unlikely, however, that different kinds of manipulations would produce the same kinds of artifacts.     

Steady states and outbreaks of two-phase nonlinear age-structured model of population dynamics with discrete time delay

In this paper we study the nonlinear age-structured model of a polycyclic two-phase population dynamics including delayed effect of population density growth on the mortality. Both phases are modelled as a system of initial boundary values problem for semi-linear transport equation with delay and initial problem for nonlinear delay ODE. The obtained system is studied both theoretically and numerically. Three different regimes of population dynamics for asymptotically stable states of autonomous systems are obtained in numerical experiments for the different initial values of population density. The quasi-periodical travelling wave solutions are studied numerically for the autonomous system with the different values of time delays and for the system with oscillating death rate and birth modulus. In both cases it is observed three types of travelling wave solutions: harmonic oscillations, pulse sequence and single pulse.     

Continuous-discrete models of the dynamics of an isolated population and of 2 competing species

The paper presents the analysis of various mathematical models for dynamics of isolated population and for competition between two species. It is assumed that mortality is continuous and birth of individuals of new generations takes place in certain fixed moments. Influence of winter upon the population dynamics and conditions of classic discrete model "deduction" of population dynamics (in particular, Moran-Ricker and Hassel's models) are investigated. Dynamic regimes of models under various assumptions about the birth and death rates upon the population states are also examined. Analysis of models of isolated population dynamics with nonoverlapping generations showed the density changes regularly if the birth rate is constant. Moreover, there exists a unique global stable level and population size stabilizes asymptotically at this equilibrium, i.e. cycle and chaotic regimes in various discrete models depend on correlation between individual productivity and population state in previous time. When the correlation is exponential upon mean population size the discrete Hassel model is realized. Modification of basis model, based on the assumption that during winter survival/death changes are constant, showed that population size at global level is stable. Generally, the dependence of population rate upon "winter parameters" has nonlinear character. Nonparametric models of competition between two species does not vary if the individual productivity is constant. In a phase space there are several stable stationary states and population stabilizes at one or other level asymptotically. So, in discrete models of competition between two species oscillation can be explained by dependence of population growth rate on the population size at previous times.     

Red environmental noise and the appearance of delayed density dependence in age-structured populations

Previous work suggests that red environmental noise can lead to the spurious appearance of delayed density dependence (DDD) in unstructured populations regulated only by direct density dependence. We analysed the effect of noise reddening on the pattern of spurious DDD in several variants of the density-dependent age-structured population model. We found patterns of spurious DDD in structured populations with either density-dependent fertility or density-dependent survival of the first age class, inconsistent with predictions from unstructured population models. Moreover, we found that nonspurious negative DDD always emerges in populations with deterministic chaotic dynamics, regardless of population structure or the type of environmental noise. The effect of noise reddening in generating spurious DDD is often negligible in the chaotic region of population deterministic dynamics. Our findings suggest that differences in species' life histories may exhibit different patterns of spurious DDD (owing to noise reddening) than predicted by unstructured models.     

Obtaining birth and mortality patterns from structured population trajectories

"A method is presented for unravelling the demographic equation for structured populations. A solution to the McKendrick-von Foerster equation is constructed using spline functions and this is fitted to stage-structured population data in such a way that the solution is smooth, positive, and does not imply negative death rates. The smoothness of the surface, and hence the complexity of the population model, is determined in a statistically optimum manner using cross validation. Time- and age-dependent death rates can be obtained as well as time-dependent birth rates. Confidence intervals are obtained for population size and death rates that give a 95% probability that the true population dynamics are within the intervals. Practical application of the method is demonstrated, and comparison made with three alternative methods."     

Demographic heterogeneity impacts density-dependent population dynamics

Among-individual variation in vital parameters such as birth and death rates that is unrelated to age, stage, sex, or environmental fluctuations is referred to as demographic heterogeneity. This kind of heterogeneity is prevalent in ecological populations, but is almost always left out of models. Demographic heterogeneity has been shown to affect demographic stochasticity in small populations and to increase growth rates for density-independent populations. The latter is due to ??cohort selection,?? where the most frail individuals die out first, lowering the cohort??s average mortality as it ages. The importance of cohort selection to population dynamics has only recently been recognized. We use a continuous-time model with density dependence, based on the logistic equation, to study the effects of demographic heterogeneity in mortality and reproduction. Reproductive heterogeneity is introduced in three ways: parent fertility, offspring viability, and parent?Coffspring correlation. We find that both the low-density growth rate and the equilibrium population size increase as the magnitude of mortality heterogeneity increases or as parent?Coffspring phenotypic correlation increases. Population dynamics are affected by complex interactions among the different types of heterogeneity, and trade-off scenarios are examined which can sometimes reverse the effect of increased heterogeneity. We show that there are a number of different homogeneous approximations to heterogeneous models, but all fail to capture important parts of the dynamics of the full model.     

Genetic instability and clonal expansion

Inactivation of tumor suppressor genes can lead to clonal expansion. We study the evolutionary dynamics of this process and calculate the probability that inactivation of a tumor suppressor gene is preceded by mutations in genes that confer genetic instability. Unstable cells might have a slower rate of clonal expansion than stable cells because of an increased probability of generating lethal mutations or inducing apoptosis. We show that the different growth rates of genetically stable and unstable cells during clonal expansion represent, in general, only a small disadvantage for genetic instability. The intuitive reason for this conclusion is that robust clonal expansion, where cellular birth rates are significantly greater than death rates, occurs on a much faster time scale than waiting for those mutations that allow clonal expansion. Moreover, in special cases where clonal expansion is very slow, genetically unstable cells have a higher probability to accumulate additional mutations during clonal expansion that confer a selective advantage. Clonal expansion represents a major disadvantage for genetic instability only when inactivation of the tumor suppressor gene leads to a very small increase of the cellular reproductive rate that is cancelled by the increased mortality of unstable cells.     

Comparing the consequences of induced and constitutive plant resistance for herbivore population dynamics

Although it has been suggested that induced and constitutive plant resistance should have different effects on insect herbivore population dynamics, there is little experimental evidence that plant resistance can influence herbivore populations longer than one season. We used a density-manipulation experiment and model fitting to examine the effects of constitutive and induced resistance on herbivore dynamics over both the short and long term. We used likelihood methods to fit population dynamic models to recruitment data for populations of Mexican bean beetles on soybean varieties with no resistance, constitutive resistance, or induced resistance. We compared model configurations that fit parameters for resistance types separately to models that did not account for resistance type. Models representing the hypothesis that the three resistance types differed in their effects on beetle dynamics received the most support. Induced resistance resulted in lower population growth rates and stronger density dependence than no resistance. Constitutive resistance resulted in lower population growth rates and stronger density dependence than induced resistance. Constitutive resistance had a stronger effect on both short-term beetle recruitment and predicted beetle population dynamics than induced resistance. The results of this study suggest that induced and constitutive resistance can differ in their effects on herbivore populations even in a relatively complex system.     

Duality,ancestral and diffusion processes in models with selection

The ancestral selection graph in population genetics was introduced by Krone and Neuhauser [Krone, S.M., Neuhauser, C., 1997. Ancestral process with selection. Theor. Popul. Biol. 51, 210–237] as an analogue of the coalescent genealogy of a sample of genes from a neutrally evolving population. The number of particles in this graph, followed backwards in time, is a birth and death process with quadratic death and linear birth rates. In this paper an explicit form of the probability distribution of the number of particles is obtained by using the density of the allele frequency in the corresponding diffusion model obtained by Kimura [Kimura, M., 1955. Stochastic process and distribution of gene frequencies under natural selection. Cold Spring Harbor Symposia on Quantitative Biology 20, 33–53]. It is shown that the process of fixation of the allele in the diffusion model corresponds to convergence of the ancestral process to its stationary measure. The time to fixation of the allele conditional on fixation is studied in terms of the ancestral process.     

The effect of self-fertilization, inbreeding depression, and population size on autopolyploid establishment

The minority cytotype exclusion principle describes how random mating between diploid and autotetraploid cytotypes hinders establishment of the rare cytotype. We present deterministic and stochastic models to ascertain how selfing, inbreeding depression, unreduced gamete production, and finite population size affect minority cytotype exclusion and the establishment of autotetraploids. Results demonstrate that higher selfing rates and lower inbreeding depression in autotetraploids facilitate establishment of autotetraploid populations. Stochastic effects due to finite population size increase the probability of polyploid establishment and decrease the mean time to tetraploid fixation. Our results extend the minority cytotype exclusion principle to include important features of plant reproduction and demonstrate that variation in mating system parameters significantly influences the conditions necessary for polyploid establishment.     

Establishing a beachhead: a stochastic population model with an Allee effect applied to species invasion

We formulated a spatially explicit stochastic population model with an Allee effect in order to explore how invasive species may become established. In our model, we varied the degree of migration between local populations and used an Allee effect with variable birth and death rates. Because of the stochastic component, population sizes below the Allee effect threshold may still have a positive probability for successful invasion. The larger the network of populations, the greater the probability of an invasion occurring when initial population sizes are close to or above the Allee threshold. Furthermore, if migration rates are low, one or more than one patch may be successfully invaded, while if migration rates are high all patches are invaded.     

Competition and evolutionary stability of plants in a spatially structured habitat

We modelled the population dynamics of two types of plants with limited dispersal living in a lattice structured habitat. Each site of the square lattice model was either occupied by an individual or vacant. Each individual reproduced to its neighbors. We derived a criterion for the invasion of a rare type into a population composed of a resident type based on a pair-approximation method, in which the dynamics of both average densities and the nearest neighbor correlations were considered. Based on this invasibility criterion, we showed that, when there is a tradeoff between birth and death rates, the evolutionarily stable type is the one that has the highest ratio of birth rate to mortality. If these types are different species, they form segregated spatial patterns in the lattice model in which intraspecific competitive interactions occur more frequently than interspecific interactions. However, stable coexistence is not possible in the lattice model contrary to results from completely mixed population models. This clearly shows that the casual conclusion, based on traditional well mixed population models, that different species can coexist if intraspecific competition is stronger than interspecific competition, does not hold for spatially structured population models.     

Tiller demography in seashore populations of Agrostis stolonifera,Festuca rubra and Poa irrigata

Abstract. The dynamics of tillers in natural populations of three cohabiting perennial grass species, Agrostis stolonifera, Festuca rubra and Poa irrigata (= Poa pratensis ssp. irrigata) were studied for five years in a Baltic seashore meadow. The process of tiller population maintenance was very dynamic. Both birth and death rates of tillers were high, particularly in A stolonifera, and the turnover rate of the populations was high. Recruitment was mainly by vegetative tillers, produced continuously throughout the growing season. The proportion of flowering tillers was low, but varied between years. Considerable year-to-year variation was also found in birth and death rates. Despite this between-year variation and the differences found between species in flowering frequency, pattern of survivorship and tiller longevity, population sizes of the species remained relatively constant.     

Evolutionary tracking is determined by differential selection on demographic rates and density dependence

Recent ecological forecasts predict that ~25% of species worldwide will go extinct by 2050. However, these estimates are primarily based on environmental changes alone and fail to incorporate important biological mechanisms such as genetic adaptation via evolution. Thus, environmental change can affect population dynamics in ways that classical frameworks can neither describe nor predict. Furthermore, often due to a lack of data, forecasting models commonly describe changes in population demography by summarizing changes in fecundity and survival concurrently with the intrinsic growth rate (r). This has been shown to be an oversimplification as the environment may impose selective pressure on specific demographic rates (birth and death) rather than directly on r (the difference between the birth and death rates). This differential pressure may alter population response to density, in each demographic rate, further diluting the information combined to produce r. Thus, when we consider the potential for persistence via adaptive evolution, populations with the same r can have different abilities to persist amidst environmental change. Therefore, we cannot adequately forecast population response to climate change without accounting for demography and selection on density dependence. Using a continuous‐time Markov chain model to describe the stochastic dynamics of the logistic model of population growth and allow for trait evolution via mutations arising during birth events, we find persistence via evolutionary tracking more likely when environmental change alters birth rather than the death rate. Furthermore, species that evolve responses to changes in the strength of density dependence due to environmental change are less vulnerable to extinction than species that undergo selection independent of population density. By incorporating these key demographic considerations into our predictive models, we can better understand how species will respond to climate change.