Transition Densities and Sample Frequency Spectra of Diffusion Processes with Selection and Variable Population Size

Advances in empirical population genetics have made apparent the need for models that simultaneously account for selection and demography. To address this need, we here study the Wright–Fisher diffusion under selection and variable effective population size. In the case of genic selection and piecewise-constant effective population sizes, we obtain the transition density by extending a recently developed method for computing an accurate spectral representation for a constant population size. Utilizing this extension, we show how to compute the sample frequency spectrum in the presence of genic selection and an arbitrary number of instantaneous changes in the effective population size. We also develop an alternate, efficient algorithm for computing the sample frequency spectrum using a moment-based approach. We apply these methods to answer the following questions: If neutrality is incorrectly assumed when there is selection, what effects does it have on demographic parameter estimation? Can the impact of negative selection be observed in populations that undergo strong exponential growth?     

Numerical solution of structured population models

Numerical methods are presented for a general age-structured population model with demographic rates depending on age and the total population size. The accuracy of these methods is established by solving problems for which alternate solution techniques are available and are used for comparison. The methods reliably solve test problems with a variety of dynamic behavior. Simulations of a blowfly population exhibit cyclic fluctuations, whereas a simulated squirrel population reaches a stable age distribution and stable equilibrium population size. Life-history attributes are easily studied from the computed solutions, and are discussed for these examples. Recovery of a stressed population back to equilibrium is examined by computing the transition in age structure, and the transient behavior of other properties of the population such as the per capita growth rate, the average age, and the generation length.     

Population growth of a gorgonian coral: equilibrium and non-equilibrium sensitivity to changes in life history variables

Summary A size dependent model of population growth of the Caribbean gorgonian Plexaura A is developed based on observed rates of survival, growth and colony fragmentation at a site in the San Blas Islands, Panama. Sensitivity and elasticity analyses indicate that the fate of large colonies has the greatest effect on population growth. Variables which directly affect the generation of large colonies have the next greatest effect on population growth. These variables include the recruitment of large fragments, and the survivorship of colonies in the next smaller size class. Sexual reproduction has an extremely limited ability to affect population growth. Vegetative reproduction has a greater potential effect on growth rates. Environmental conditions regularly change the matrix of transition probabilities which predicts population growth. This keeps the population from approaching its stable size class distribution. Deviations from the stable size class distribution alter sensitivity and elasticity and in this case have the effect of increasing the importance of survivorship of the smallest colonies. Nonequilibrium conditions alter sensitivity analyses and it is important to assess whether populations are at equilibrium and to determine the effects of such deviations on the sensitivity analysis.     

Population and prehistory II: space-limited human populations in constant environments

We present a population model to examine the forces that determined the quality and quantity of human life in early agricultural societies where cultivable area is limited. The model is driven by the non-linear and interdependent relationships between the age distribution of a population, its behavior and technology, and the nature of its environment. The common currency in the model is the production of food, on which age-specific rates of birth and death depend. There is a single non-trivial equilibrium population at which productivity balances caloric needs. One of the most powerful controls on equilibrium hunger level is fertility control. Gains against hunger are accompanied by decreases in population size. Increasing worker productivity does increase equilibrium population size but does not improve welfare at equilibrium. As a case study we apply the model to the population of a Polynesian valley before European contact.     

Unimodal Tree Size Distributions Possibly Result from Relatively Strong Conservatism in Intermediate Size Classes

Tree size distributions have long been of interest to ecologists and foresters because they reflect fundamental demographic processes. Previous studies have assumed that size distributions are often associated with population trends or with the degree of shade tolerance. We tested these associations for 31 tree species in a 20 ha plot in a Dinghushan south subtropical forest in China. These species varied widely in growth form and shade-tolerance. We used 2005 and 2010 census data from that plot. We found that 23 species had reversed J shaped size distributions, and eight species had unimodal size distributions in 2005. On average, modal species had lower recruitment rates than reversed J species, while showing no significant difference in mortality rates, per capita population growth rates or shade-tolerance. We compared the observed size distributions with the equilibrium distributions projected from observed size-dependent growth and mortality. We found that observed distributions generally had the same shape as predicted equilibrium distributions in both unimodal and reversed J species, but there were statistically significant, important quantitative differences between observed and projected equilibrium size distributions in most species, suggesting that these populations are not at equilibrium and that this forest is changing over time. Almost all modal species had U-shaped size-dependent mortality and/or growth functions, with turning points of both mortality and growth at intermediate size classes close to the peak in the size distribution. These results show that modal size distributions do not necessarily indicate either population decline or shade-intolerance. Instead, the modal species in our study were characterized by a life history strategy of relatively strong conservatism in an intermediate size class, leading to very low growth and mortality in that size class, and thus to a peak in the size distribution at intermediate sizes.     

Stability of an age specific population with density dependent fertility

A nonlinear version of the Lotka-Sharpe model of population growth is considered in which the age specific fertility is a function of the population size. The stability of an equilibrium population distribution is investigated with respect to both global and local perturbations. Sufficient conditions for such stability are presented, as are estimates for the rate of return of the population to the equilibrium configuration. Particular attention is paid to those situations in which the age dependent stability criteria coincide with those of age independent models.     

Population dynamics in two-patch environments: Some anomalous consequences of an optimal habitat distribution

The effect of dispersal on population size and stability is explored for a population that disperses passively between two discrete habitat patches. Two basic models are considered. In the first model, a single population experiences density-dependent growth in both patches. A graphical construction is presented which allows one to determine the spatial pattern of abundance at equilibrium for most reasonable growth models and rates of dispersal. It is shown under rather general conditions that this equilibrium is unique and globally stable. In the second model, the dispersing population is a food-limited predator that occurs in both a source habitat (which contains a prey population) and a sink habitat (which does not). Passive dispersal between source and sink habitats can stabilize an otherwise unstable predator-prey interaction. The conditions allowing this are explored in some detail. The theory of optimal habitat selection predicts the evolutionarily stable distribution of a population, given that individuals can freely move among habitats so as to maximize individual fitness. This theory is used to develop a heuristic argument for why passive dispersal should always be selectively disadvantageous (ignoring kin effects) in a spatially heterogeneous but temporally constant environment. For both the models considered here, passive dispersal may lead to a greater number of individuals in both habitats combined than if there were no dispersal. This implies that the evolution of an optimal habitat distribution may lead to a reduction in population size; in the case of the predator-prey model, it may have the additional effect of destabilizing the interaction. The paper concludes with a discussion of the disparate effects habitat selection might have on the geographical range occupied by a species.     

Estimating population birth rates of zooplankton when rates of egg deposition and hatching are periodic

Summary I present a general method of computing finite birth and death rates of natural zooplankton populations from changes in the age distribution of eggs and changes in population size. The method is applicable to cases in which eggs hatch periodically owing to variable rates of oviposition. When morphological criteria are used to determine the age distribution of eggs at the beginning and end of a sampling interval, egg mortality can be incorporated in estimates of population birth rate. I raised laboratory populations of Asplanchna priodonta, a common planktonic rotifer, in semicontinuous culture to evaluate my method of computing finite birth rate. The Asplanchna population became synchronized to a daily addition of food but grew by the same amount each day once steady state was achieved. The steady-state rate of growth, which can be computed from the volume-specific dilution rate of the culture, was consistent with the finite birth rate predicted from the population's egg ratio and egg age distribution.     

Population growth makes waves in the distribution of pairwise genetic differences.

Episodes of population growth and decline leave characteristic signatures in the distribution of nucleotide (or restriction) site differences between pairs of individuals. These signatures appear in histograms showing the relative frequencies of pairs of individuals who differ by i sites, where i = 0, 1, .... In this distribution an episode of growth generates a wave that travels to the right, traversing 1 unit of the horizontal axis in each 1/2u generations, where u is the mutation rate. The smaller the initial population, the steeper will be the leading face of the wave. The larger the increase in population size, the smaller will be the distribution's vertical intercept. The implications of continued exponential growth are indistinguishable from those of a sudden burst of population growth Bottlenecks in population size also generate waves similar to those produced by a sudden expansion, but with elevated uppertail probabilities. Reductions in population size initially generate L-shaped distributions with high probability of identity, but these converge rapidly to a new equilibrium. In equilibrium populations the theoretical curves are free of waves. However, computer simulations of such populations generate empirical distributions with many peaks and little resemblance to the theory. On the other hand, agreement is better in the transient (nonequilibrium) case, where simulated empirical distributions typically exhibit waves very similar to those predicted by theory. Thus, waves in empirical distributions may be rich in information about the history of population dynamics.     

The mid-depth method and HIV-1: a practical approach for testing hypotheses of viral epidemic history.

We introduce the mid-depth method, a practical approach for testing hypotheses of demographic history using genealogies reconstructed from sequence data. The relative positions of internal nodes within a genealogy contain information about past population dynamics. We explain how this information can be used to (1) test the null hypothesis of constant population size and (2) estimate the growth rate and current population size of an exponentially growing population. Simulation tests indicate that, as expected, estimates of exponential growth rates are sometimes biased. The mid-depth method is computationally rapid and does not require knowledge of the sample's mutation rate. However, it does assume that the reconstructed genealogy is correct and is therefore best suited to the analysis of variation-rich viral data sets. When applied to HIV-1 sequence data, the mid-depth method provides phylogenetic evidence of different exponential growth rates for subtypes A and B. We posit that this difference in growth rate reflects the different transmission routes and epidemiological histories of the two subtypes.     

Ontogenetic symmetry and asymmetry in energetics

Body size ( $\equiv $ biomass) is the dominant determinant of population dynamical processes such as giving birth or dying in almost all species, with often drastically different behaviour occurring in different parts of the growth trajectory, while the latter is largely determined by food availability at the different life stages. This leads to the question under what conditions unstructured population models, formulated in terms of total population biomass, still do a fair job. To contribute to answering this question we first analyze the conditions under which a size-structured model collapses to a dynamically equivalent unstructured one in terms of total biomass. The only biologically meaningful case where this occurs is when body size does not affect any of the population dynamic processes, this is the case if and only if the mass-specific ingestion rate, the mass-specific biomass production and the mortality rate of the individuals are independent of size, a condition to which we refer as “ontogenetic symmetry”. Intriguingly, under ontogenetic symmetry the equilibrium biomass-body size spectrum is proportional to 1/size, a form that has been conjectured for marine size spectra and subsequently has been used as prior assumption in theoretical papers dealing with the latter. As a next step we consider an archetypical class of models in which reproduction takes over from growth upon reaching an adult body size, in order to determine how quickly discrepancies from ontogenetic symmetry lead to relevant novel population dynamical phenomena. The phenomena considered are biomass overcompensation, when additional imposed mortality leads, rather unexpectedly, to an increase in the equilibrium biomass of either the juveniles or the adults (a phenomenon with potentially big consequences for predators of the species), and the occurrence of two types of size-structure driven oscillations, juvenile-driven cycles with separated extended cohorts, and adult-driven cycles in which periodically a front of relatively steeply decreasing frequencies moves up the size distribution. A small discrepancy from symmetry can already lead to biomass overcompensation; size-structure driven cycles only occur for somewhat larger discrepancies.     

Dynamics of an age-structured population with Allee effects and harvesting

We propose a discrete-time, age-structured population model to study the impact of Allee effects and harvesting. It is assumed that survival probabilities from one age class to the next are constants and fertility rate is a function of weighted total population size. Global extinction is certain if the maximal growth rate of the population is less than one. The model can have multiple attractors and the asymptotic dynamics of the population depend on its initial distribution if the maximal growth rate is larger than one. An Allee threshold depending on the components of the unstable interior equilibrium is derived when only the last age class can reproduce. The population becomes extinct if its initial population distribution is below the threshold. Harvesting on any particular age class can decrease the magnitude of the possible stable interior equilibrium and increase the magnitude of the unstable interior equilibrium simultaneously.     

Dynamics of an age-structured population with Allee effects and harvesting

We propose a discrete-time, age-structured population model to study the impact of Allee effects and harvesting. It is assumed that survival probabilities from one age class to the next are constants and fertility rate is a function of weighted total population size. Global extinction is certain if the maximal growth rate of the population is less than one. The model can have multiple attractors and the asymptotic dynamics of the population depend on its initial distribution if the maximal growth rate is larger than one. An Allee threshold depending on the components of the unstable interior equilibrium is derived when only the last age class can reproduce. The population becomes extinct if its initial population distribution is below the threshold. Harvesting on any particular age class can decrease the magnitude of the possible stable interior equilibrium and increase the magnitude of the unstable interior equilibrium simultaneously.     

How can we model selectively neutral density dependence in evolutionary games

The problem of density dependence appears in all approaches to the modelling of population dynamics. It is pertinent to classic models (i.e., Lotka-Volterra's), and also population genetics and game theoretical models related to the replicator dynamics. There is no density dependence in the classic formulation of replicator dynamics, which means that population size may grow to infinity. Therefore the question arises: How is unlimited population growth suppressed in frequency-dependent models? Two categories of solutions can be found in the literature. In the first, replicator dynamics is independent of background fitness. In the second type of solution, a multiplicative suppression coefficient is used, as in a logistic equation. Both approaches have disadvantages. The first one is incompatible with the methods of life history theory and basic probabilistic intuitions. The logistic type of suppression of per capita growth rate stops trajectories of selection when population size reaches the maximal value (carrying capacity); hence this method does not satisfy selective neutrality. To overcome these difficulties, we must explicitly consider turn-over of individuals dependent on mortality rate. This new approach leads to two interesting predictions. First, the equilibrium value of population size is lower than carrying capacity and depends on the mortality rate. Second, although the phase portrait of selection trajectories is the same as in density-independent replicator dynamics, pace of selection slows down when population size approaches equilibrium, and then remains constant and dependent on the rate of turn-over of individuals.     

Polymorphism in sexual versus non-sexual disease transmission

Pathogens causing sexually transmitted diseases (STDs) often consist of related strains that cause non-sexually transmitted, or ''ordinary infectious'', diseases (OIDs). We use differential equation models of single populations to derive conditions under which a genetic variant with one (e.g. sexual) transmission mode can invade and successfully displace a genetic variant with a different (e.g. non-sexual) transmission mode. Invasion by an STD is easier if the equilibrium population size in the presence of an OID is smaller; conversely an OID can invade more easily if the equilibrium size of the population with the STD is larger. Invasion of an STD does not depend on the degree of sterility caused by the infection, but does depend on the added mortality caused by a resident OID. In contrast, the ability of an OID to invade a population at equilibrium with an STD decreases as the degree of sterility caused by the STD increases. When equilibrium population sizes for a population infected with an STD are above the point at which non-sexual contacts exceed sexual contacts (the sexual–social crossover point) and when equilibrium population sizes for an OID are below this point, there can be a stable genetic polymorphism for transmission mode. This is most likely when the STD is mildly sterilizing, and the OID causes low or intermediate levels of added mortality. Because we assume the strains are competitively equivalent and there are no heterogeneities associated with the transmission process, the polymorphism is maintained by density-dependent selection brought about by pathogen effects on population size.     

Defining the fire trap: Extension of the persistence equilibrium model in mesic savannas

Mesic savannas are dominated by trees that are strong resprouters caught in a frequent fire trap. Persistence within this fire trap has been described by a resprout curve of SizeNext ~ f(Pre‐fire size), defined by the Michaelis‐Menten function. A key feature of this resprout curve is a stable persistence equilibrium that represents the size of individual plants upon which a population will converge over successive inter‐fire time steps under a given fire regime. Here, we contend that such a resprout curve does not adequately describe resprout tree dynamics in frequently burnt mesic savannas because it is constrained to an asymptote. We propose a new framework for modelling the resprout curve, which recognizes that local environmental stochasticity and growth patterns can interact to change the growth response function entirely, and thus more readily reflect the range of feasible resprout responses. Importantly, we define an unstable equilibrium representing the size above which individuals have escaped the fire trap and explore mechanisms that can shift an individual from persistence to escape. Through a case study from northern Australia, we confirm that our framework provides a simple yet practical approach to defining these critical aspects of savanna tree growth dynamics: persistence and escape.     

Reconstructing the dynamics of ancient human populations from radiocarbon dates: 10 000 years of population growth in Australia

Measuring trends in the size of prehistoric populations is fundamental to our understanding of the demography of ancient people and their responses to environmental change. Archaeologists commonly use the temporal distribution of radiocarbon dates to reconstruct population trends, but this can give a false picture of population growth because of the loss of evidence from older sites. We demonstrate a method for quantifying this bias, and we use it to test for population growth through the Holocene of Australia. We used model simulations to show how turnover of site occupation across an archaeological landscape, interacting with erasure of evidence at abandoned sites, can create an increase in apparent site occupation towards the present when occupation density is actually constant. By estimating the probabilities of abandonment and erasure from archaeological data, we then used the model to show that this effect does not account for the observed increase in occupation through the Holocene in Australia. This is best explained by population growth, which was low for the first part of the Holocene but accelerated about 5000 years ago. Our results provide new evidence for the dynamism of non-agricultural populations through the Holocene.     

How to compute the effective size of spatiotemporally structured populations using separation of time scales

Calculations to derive effective population size become highly complicated when complex population structure is considered. We provide an easy method of computing the effective size of a subdivided population with overlapping generations (a spatiotemporally structured population) using an approximation based on separation of time scales. We also numerically compute the effective size to verify the accuracy of the derived formula. Various interesting quantities, including moments of coalescent time, are readily derived using this approach.     

On stochastic logistic population growth models with immigration and multiple births

This paper develops a stochastic logistic population growth model with immigration and multiple births. The differential equations for the low-order cumulant functions (i.e., mean, variance, and skewness) of the single birth model are reviewed, and the corresponding equations for the multiple birth model are derived. Accurate approximate solutions for the cumulant functions are obtained using moment closure methods for two families of model parameterizations, one for badger and the other for fox population growth. For both model families, the equilibrium size distribution may be approximated well using the Normal approximation, and even more accurately using the saddlepoint approximation. It is shown that in comparison with the corresponding single birth model, the multiple birth mechanism increases the skewness and the variance of the equilibrium distribution, but slightly reduces its mean. Moreover, the type of density-dependent population control is shown to influence the sign of the skewness and the size of the variance.     

Effects of diffusion on total biomass in heterogeneous continuous and discrete-patch systems

Theoretical models of populations on a system of two connected patches previously have shown that when the two patches differ in maximum growth rate and carrying capacity, and in the limit of high diffusion, conditions exist for which the total population size at equilibrium exceeds that of the ideal free distribution, which predicts that the total population would equal the total carrying capacity of the two patches. However, this result has only been shown for the Pearl-Verhulst growth function on two patches and for a single-parameter growth function in continuous space. Here, we provide a general criterion for total population size to exceed total carrying capacity for three commonly used population growth rates for both heterogeneous continuous and multi-patch heterogeneous landscapes with high population diffusion. We show that a sufficient condition for this situation is that there is a convex positive relationship between the maximum growth rate and the parameter that, by itself or together with the maximum growth rate, determines the carrying capacity, as both vary across a spatial region. This relationship occurs in some biological populations, though not in others, so the result has ecological implications.