A branching process, its application in biology: influence of demographic parameters on the social structure in mammal groups

Branching processes are widely used in biology. This theoretical tool is used in cell dynamics, epidemics and population dynamics. In population dynamics, branching processes are mainly used to access extinction probabilities of populations, groups or families, with the Galton-Watson branching process. Many mammal species live in socially-structured groups, and the smallest units of these groups are lineages (or families) of kin-related individuals. In many primate species, these lineages are matrilines, as females remain in their natal groups most of the time, whereas males generally disperse. Lineage parameters, such as numbers of matrilines, size of each matriline and average degree of relatedness, could strongly influence the genetic composition of groups. Evidence indicates that division along matrilines could induce substantial differentiation among fission groups. Here, we develop a novel mathematical model based on the branching process theory describing demographic dynamics of groups. The main result of this model is an explicit analytical expression of the joint distribution of numbers of lineages and sizes of socially-structured groups. We investigated the influence of parameters such as natality and mortality on the outcome of the process, including extinction probability. Finally, we discuss this theoretical result with respect to biological significance.     

A note on modelling the dynamics of budding yeast populations using branching processes

A multitype branching process is proposed as a model for the behaviour of populations of the budding yeast Saccharomyces Cerevisiae. Using the idea of branching processes counted by random characteristics, we are able to obtain explicit expressions describing different aspects of the asymptotic composition of such populations. The main purpose of this note is to show that the branching process approach is an alternative to deterministic population models based on differential equation methods.Supported by the Swedish Natural Science Research Council     

Population Genetic Consequences of the Allee Effect and the Role of Offspring-Number Variation

A strong demographic Allee effect in which the expected population growth rate is negative below a certain critical population size can cause high extinction probabilities in small introduced populations. But many species are repeatedly introduced to the same location and eventually one population may overcome the Allee effect by chance. With the help of stochastic models, we investigate how much genetic diversity such successful populations harbor on average and how this depends on offspring-number variation, an important source of stochastic variability in population size. We find that with increasing variability, the Allee effect increasingly promotes genetic diversity in successful populations. Successful Allee-effect populations with highly variable population dynamics escape rapidly from the region of small population sizes and do not linger around the critical population size. Therefore, they are exposed to relatively little genetic drift. It is also conceivable, however, that an Allee effect itself leads to an increase in offspring-number variation. In this case, successful populations with an Allee effect can exhibit less genetic diversity despite growing faster at small population sizes. Unlike in many classical population genetics models, the role of offspring-number variation for the population genetic consequences of the Allee effect cannot be accounted for by an effective-population-size correction. Thus, our results highlight the importance of detailed biological knowledge, in this case on the probability distribution of family sizes, when predicting the evolutionary potential of newly founded populations or when using genetic data to reconstruct their demographic history.     

Evolutionary Branching in a Finite Population: Deterministic Branching vs. Stochastic Branching

Adaptive dynamics formalism demonstrates that, in a constant environment, a continuous trait may first converge to a singular point followed by spontaneous transition from a unimodal trait distribution into a bimodal one, which is called “evolutionary branching.” Most previous analyses of evolutionary branching have been conducted in an infinitely large population. Here, we study the effect of stochasticity caused by the finiteness of the population size on evolutionary branching. By analyzing the dynamics of trait variance, we obtain the condition for evolutionary branching as the one under which trait variance explodes. Genetic drift reduces the trait variance and causes stochastic fluctuation. In a very small population, evolutionary branching does not occur. In larger populations, evolutionary branching may occur, but it occurs in two different manners: in deterministic branching, branching occurs quickly when the population reaches the singular point, while in stochastic branching, the population stays at singularity for a period before branching out. The conditions for these cases and the mean branching-out times are calculated in terms of population size, mutational effects, and selection intensity and are confirmed by direct computer simulations of the individual-based model.     

EVOLUTION OF DISPERSAL RATES IN METAPOPULATION MODELS: BRANCHING AND CYCLIC DYNAMICS IN PHENOTYPE SPACE

We study the evolution of dispersal rates in a two patch metapopulation model. The local dynamics in each patch are given by difference equations, which, together with the rate of dispersal between the patches, determine the ecological dynamics of the metapopulation. We assume that phenotypes are given by their dispersal rate. The evolutionary dynamics in phenotype space are determined by invasion exponents, which describe whether a mutant can invade a given resident population. If the resident metapopulation is at a stable equilibrium, then selection on dispersal rates is neutral if the population sizes in the two patches are the same, while selection drives dispersal rates to zero if the local abundances are different. With non-equilibrium metapopulation dynamics, non-zero dispersal rates can be maintained by selection. In this case, and if the patches are ecologically identical, dispersal rates always evolve to values which induce synchronized metapopulation dynamics. If the patches are ecologically different, evolutionary branching into two coexisting dispersal phenotypes can be observed. Such branching can happen repeatedly, leading to polymorphisms with more than two phenotypes. If there is a cost to dispersal, evolutionary cycling in phenotype space can occur due to the dependence of selection pressures on the ecological attractor of the resident population, or because phenotypic branching alternates with the extinction of one of the branches. Our results extend those of Holt and McPeek (1996), and suggest that phenotypic branching is an important evolutionary process. This process may be relevant for sympatric speciation.     

A tale of two cycles – distinguishing quasi-cycles and limit cycles in finite predator–prey populations

Periodic predator – prey dynamics in constant environments are usually taken as indicative of deterministic limit cycles. It is known, however, that demographic stochasticity in finite populations can also give rise to regular population cycles, even when the corresponding deterministic models predict a stable equilibrium. Specifically, such quasi-cycles are expected in stochastic versions of deterministic models exhibiting equilibrium dynamics with weakly damped oscillations. The existence of quasi-cycles substantially expands the scope for natural patterns of periodic population oscillations caused by ecological interactions, thereby complicating the conclusive interpretation of such patterns. Here we show how to distinguish between quasi-cycles and noisy limit cycles based on observing changing population sizes in predator – prey populations. We start by confirming that both types of cycle can occur in the individual-based version of a widely used class of deterministic predator – prey model. We then show that it is feasible and straightforward to accurately distinguish between the two types of cycle through the combined analysis of autocorrelations and marginal distributions of population sizes. Finally, by confronting these results with real ecological time series, we demonstrate that by using our methods even short and imperfect time series allow quasi-cycles and limit cycles to be distinguished reliably.     

The role of weak selection and high mutation rates in nearly neutral evolution

Neutral dynamics occur in evolution if all types are ‘effectively equal’ in their reproductive success, where the definition of ‘effectively equal’ depends on the population size and the details of mutations. Empirically observed neutral genetic evolution in extremely large clonal populations can only be explained under current models if selection is completely absent. Such models typically consider the case where population dynamics occurs on a different timescale to evolution. However, this assumption is invalid when mutations are not rare in a whole population. We show that this has important consequences for the occurrence of neutral evolution in clonal populations. In highly connected type spaces, neutral dynamics can occur for all population sizes despite significant selective differences, via the forming of effectively neutral networks connecting rare neutral types. Biological implications include an explanation for the high diversity of rare types that survive in large clonal populations, and a theoretical justification for the use of neutral null models.     

On the importance of evolving phenotype distributions on evolutionary diversification

Evolutionary branching occurs when a population with a unimodal phenotype distribution diversifies into a multimodally distributed population consisting of two or more strains. Branching results from frequency-dependent selection, which is caused by interactions between individuals. For example, a population performing a social task may diversify into a cooperator strain and a defector strain. Branching can also occur in multi-dimensional phenotype spaces, such as when two tasks are performed simultaneously. In such cases, the strains may diverge in different directions: possible outcomes include division of labor (with each population performing one of the tasks) or the diversification into a strain that performs both tasks and another that performs neither. Here we show that the shape of the population’s phenotypic distribution plays a role in determining the direction of branching. Furthermore, we show that the shape of the distribution is, in turn, contingent on the direction of approach to the evolutionary branching point. This results in a distribution–selection feedback that is not captured in analytical models of evolutionary branching, which assume monomorphic populations. Finally, we show that this feedback can influence long-term evolutionary dynamics and promote the evolution of division of labor.     

Local environment and density-dependent feedbacks determine population growth in a forest herb

Linking spatial variation in environmental factors to variation in demographic rates is essential for a mechanistic understanding of the dynamics of populations. However, we still know relatively little about such links, partly because feedbacks via intraspecific density make them difficult to observe in natural populations. We conducted a detailed field study and investigated simultaneous effects of environmental factors and the intraspecific density of individuals on the demography of the herb Lathyrus vernus. In regression models of vital rates we identified effects associated with spring shade on survival and growth, while density was negatively correlated with these vital rates. Density was also negatively correlated with average individual size in the study plots, which is consistent with self-thinning. In addition, average plant sizes were larger than predicted by density in plots that were less shaded by the tree canopy, indicating an environmentally determined carrying capacity. A size-structured integral projection model based on the vital rate regressions revealed that the identified effects of shade and density were strong enough to produce differences in stable population sizes similar to those observed in the field. The results illustrate how the local environment can determine dynamics of populations and that intraspecific density may have to be more carefully considered in studies of plant demography and population viability analyses of threatened species. We conclude that demographic approaches incorporating information about both density and key environmental factors are powerful tools for understanding the processes that interact to determine population dynamics and abundances.     

Phylodynamic Inference for Structured Epidemiological Models

Coalescent theory is routinely used to estimate past population dynamics and demographic parameters from genealogies. While early work in coalescent theory only considered simple demographic models, advances in theory have allowed for increasingly complex demographic scenarios to be considered. The success of this approach has lead to coalescent-based inference methods being applied to populations with rapidly changing population dynamics, including pathogens like RNA viruses. However, fitting epidemiological models to genealogies via coalescent models remains a challenging task, because pathogen populations often exhibit complex, nonlinear dynamics and are structured by multiple factors. Moreover, it often becomes necessary to consider stochastic variation in population dynamics when fitting such complex models to real data. Using recently developed structured coalescent models that accommodate complex population dynamics and population structure, we develop a statistical framework for fitting stochastic epidemiological models to genealogies. By combining particle filtering methods with Bayesian Markov chain Monte Carlo methods, we are able to fit a wide class of stochastic, nonlinear epidemiological models with different forms of population structure to genealogies. We demonstrate our framework using two structured epidemiological models: a model with disease progression between multiple stages of infection and a two-population model reflecting spatial structure. We apply the multi-stage model to HIV genealogies and show that the proposed method can be used to estimate the stage-specific transmission rates and prevalence of HIV. Finally, using the two-population model we explore how much information about population structure is contained in genealogies and what sample sizes are necessary to reliably infer parameters like migration rates.     

Demographic response of plant populations to habitat fragmentation and temporal environmental variability

Many plant species currently exist in fragmented populations of different sizes, while they also experience unpredictable climatic fluctuation over time. However, we still understand little about how plant demography responds to such spatial and temporal environmental variability. We studied population dynamics of an understory perennial herb Trillium camschatcense in the Tokachi plain of Hokkaido, Japan, where a significant effect of forest fragmentation on seedling recruitment was previously reported. Four populations across a range of fragment sizes were studied for 6 years, and the data were analyzed using matrix population models. Per capita fecundity (the number of recruits per plant) varied greatly among populations, but the variation in population growth rates (λ) was mainly driven by the variation in stasis and growth rates, suggesting that the general trend of reduced fecundity in fragmented populations may not be readily translated into subsequent dynamics. Temporal variation in λ among years was more than 2 times larger than spatial variation among populations, and this result was likely attributable to the contrasting response of correlation structures among demographic rates. The among-population variation in λ was dampened by negative covariation between matrix elements possibly due to density-dependent regulation as well as an inherent constraint that some elements are not independent, whereas positive covariation between matrix elements resulted in large temporal variation in λ. Our results show that population dynamics responded differently to habitat fragmentation and temporal variability of the environment, emphasizing the need to discriminate these spatial and temporal variations in demographic models. Although no populations were projected to be declining in stochastic simulations, correlation between current habitat size and plant density implies historical λ is positively related to habitat size.     

Occupancy versus colonization–extinction models for projecting population trends at different spatial scales

Understanding spatiotemporal population trends and their drivers is a key aim in population ecology. We further need to be able to predict how the dynamics and sizes of populations are affected in the long term by changing landscapes and climate. However, predictions of future population trends are sensitive to a range of modeling assumptions. Deadwood‐dependent fungi are an excellent system for testing the performance of different predictive models of sessile species as these species have different rarity and spatial population dynamics, the populations are structured at different spatial scales, and they utilize distinct substrates. We tested how the projected large‐scale occupancies of species with differing landscape‐scale occupancies are affected over the coming century by different modeling assumptions. We compared projections based on occupancy models against colonization–extinction models, conducting the modeling at alternative spatial scales and using fine‐ or coarse‐resolution deadwood data. We also tested effects of key explanatory variables on species occurrence and colonization–extinction dynamics. The hierarchical Bayesian models applied were fitted to an extensive repeated survey of deadwood and fungi at 174 patches. We projected higher occurrence probabilities and more positive trends using the occupancy models compared to the colonization–extinction models, with greater difference for the species with lower occupancy, colonization rate, and colonization:extinction ratio than for the species with higher estimates of these statistics. The magnitude of future increase in occupancy depended strongly on the spatial modeling scale and resource resolution. We encourage using colonization–extinction models over occupancy models, modeling the process at the finest resource‐unit resolution that is utilizable by the species, and conducting projections for the same spatial scale and resource resolution at which the model fitting is conducted. Further, the models applied should include key variables driving the metapopulation dynamics, such as the availability of suitable resource units, habitat quality, and spatial connectivity.     

High host density favors greater virulence: a model of parasite–host dynamics based on multi-type branching processes

We use a multitype continuous time Markov branching process model to describe the dynamics of the spread of parasites of two types that can mutate into each other in a common host population. While most mathematical models for the virulence of infectious diseases focus on the interplay between the dynamics of host populations and the optimal characteristics for the success of the pathogen, our model focuses on how pathogen characteristics may change at the start of an epidemic, before the density of susceptible hosts decline. We envisage animal husbandry situations where hosts are at very high density and epidemics are curtailed before host densities are much reduced. The use of three pathogen characteristics: lethality, transmissibility and mutability allows us to investigate the interplay of these in relation to host density. We provide some numerical illustrations and discuss the effects of the size of the enclosure containing the host population on the encounter rate in our model that plays the key role in determining what pathogen type will eventually prevail. We also present a multistage extension of the model to situations where there are several populations and parasites can be transmitted from one of them to another. We conclude that animal husbandry situations with high stock densities will lead to very rapid increases in virulence, where virulent strains are either more transmissible or favoured by mutation. Further the process is affected by the nature of the farm enclosures.     

Olive Fruit Fly (Bactrocera oleae) Population Dynamics in the Eastern Mediterranean: Influence of Exogenous Uncertainty on a Monophagous Frugivorous Insect

Despite of the economic importance of the olive fly (Bactrocera oleae) and the large amount of biological and ecological studies on the insect, the factors driving its population dynamics (i.e., population persistence and regulation) had not been analytically investigated until the present study. Specifically, our study investigated the autoregressive process of the olive fly populations, and the joint role of intrinsic and extrinsic factors molding the population dynamics of the insect. Accounting for endogenous dynamics and the influences of exogenous factors such as olive grove temperature, the North Atlantic Oscillation and the presence of potential host fruit, we modeled olive fly populations in five locations in the Eastern Mediterranean region. Our models indicate that the rate of population change is mainly shaped by first and higher order non-monotonic, endogenous dynamics (i.e., density-dependent population feedback). The olive grove temperature was the main exogenous driver, while the North Atlantic Oscillation and fruit availability acted as significant exogenous factors in one of the five populations. Seasonal influences were also relevant for three of the populations. In spite of exogenous effects, the rate of population change was fairly stable along time. We propose that a special reproductive mechanism, such as reproductive quiescence, allows populations of monophagous fruit flies such as the olive fly to remain stable. Further, we discuss how weather factors could impinge constraints on the population dynamics at the local level. Particularly, local temperature dynamics could provide forecasting cues for management guidelines. Jointly, our results advocate for establishing monitoring programs and for a major focus of research on the relationship between life history traits and populations dynamics.     

Effect of dispersal at range edges on the structure of species ranges

Range edges are of particular interest to ecology because they hold key insights into the limits of the realized niche and associated population dynamics. A recent feature of Oikos summarized the state of the art on range edge ecology. While the typical question is what causes range edges, another important question is how range edges influence the distribution of abundances across a species geographic range when dispersal is present. We used a single species population dynamics model on a coupled-lattice to determine the effects of dispersal on peripheral populations as compared to populations at the core of the range. In the absence of resource gradients, the reduced neighborhood and thus lower connectivity or higher isolation among populations at the range edge alone led to significantly lower population sizes in the periphery of the range than in the core. Lower population sizes mean higher extinction risks and lower adaptability at the range edge, which could inhibit or slow range expansions, and thus effectively stabilize range edges. The strength of this effect depended on the potential population growth rate and the maximum dispersal distance. Lower potential population growth rates led to a stronger effect of dispersal resulting in a higher difference in population sizes between the two areas. The differential effect of dispersal on population sizes at the core and periphery of the range in the absence of resource gradients implies that traditional, habitat-based distribution models could result in misleading conclusions about the habitat quality in the periphery. Lower population sizes at the periphery are also relevant to conservation, because habitat removal not only eliminates populations but also creates new edges. Populations bordering these new edges may experience declines, due to their increased isolation.     

Estimating spatial coupling in epidemiological systems: a mechanistic approach

In recent years, ecologists and epidemiologists have paid increasing attention to the influence of spatial structure in shaping the dynamics and determining the persistence of populations. This is fundamentally affected by the concept of 'coupling'—the flux of individuals moving between separate populations. In this paper, we contrast how coupling is typically implemented in epidemic models with more detailed approaches. Our aim is to link the popular phenomenological formulations with the results of mechanistic models. By concentrating on the behaviour of simple epidemic systems, we relate explicit movement patterns with observed levels of coupling, validating the standard formulation. The analysis is then extended to include a brief study of how the correlation between stochastic populations is affected by coupling, the underlying deterministic dynamics and the relative population sizes.     

Metapopulation models for historical inference

The genealogical process for a sample from a metapopulation, in which local populations are connected by migration and can undergo extinction and subsequent recolonization, is shown to have a relatively simple structure in the limit as the number of populations in the metapopulation approaches infinity. The result, which is an approximation to the ancestral behaviour of samples from a metapopulation with a large number of populations, is the same as that previously described for other metapopulation models, namely that the genealogical process is closely related to Kingman's unstructured coalescent. The present work considers a more general class of models that includes two kinds of extinction and recolonization, and the possibility that gamete production precedes extinction. In addition, following other recent work, this result for a metapopulation divided into many populations is shown to hold both for finite population sizes and in the usual diffusion limit, which assumes that population sizes are large. Examples illustrate when the usual diffusion limit is appropriate and when it is not. Some shortcomings and extensions of the model are considered, and the relevance of such models to understanding human history is discussed.     

Dynamics of a feline retrovirus (FeLV) in host populations with variable spatial structure.

The predictions of epidemic models are remarkably affected by the underlying assumptions concerning host population dynamics and the relation between host density and disease transmission. Furthermore, hypotheses underlying distinct models are rarely tested. Domestic cats (Felis catus) can be used to compare models and test their predictions, because cat populations show variable spatial structure that probably results in variability in the relation between density and disease transmission. Cat populations also exhibit various dynamics. We compare four epidemiological models of Feline Leukaemia Virus (FeLV). We use two different incidence terms, i.e. proportionate mixing and pseudo-mass action. Population dynamics are modelled as logistic or exponential growth. Compared with proportionate mixing, mass action incidence with logistic growth results in a threshold population size under which the virus cannot persist in the population. Exponential growth of host populations results in systems where FeLV persistence at a steady prevalence and depression of host population growth are biologically unlikely to occur. Predictions of our models account for presently available data on FeLV dynamics in various populations of cats. Thus, host population dynamics and spatial structure can be determinant parameters in parasite transmission, host population depression, and disease control.     

Life history and population size variability in a relict plant. Different routes towards long-term persistence

A central tenet of conservation biology is that population size affects the persistence of populations. However, many narrow endemic species combine small population ranges and sizes with long persistence, thereby challenging this tenet. I examined the performance of three different-sized populations of Petrocoptis pseudoviscosa (Caryophyllaceae), a palaeoendemic rupicolous herb distributed along a small valley in the Spanish Pyrenees. Reproductive and demographic parameters were recorded over 6 years, and deterministic and stochastic matrix models were constructed to explore population dynamics and extinction risk. Populations differed greatly in structure, fecundity, recruitment, survival rate, and life span. Strong differentiation in life-history parameters and their temporal variability resulted in differential population vulnerability under current conditions and simulated global changes such as habitat fragmentation or higher climatic fluctuations. This study provides insights into the capacity of narrow endemics to survive both at extreme environmental conditions and at small population sizes. When dealing with species conservation, the population size–extinction risk relationship may be too simplistic for ancient, ecologically restricted organisms, and some knowledge of life history may be most important to assess their future.     

MAXIMUM-LIKELIHOOD ESTIMATION OF POPULATION DIVERGENCE TIMES AND POPULATION PHYLOGENY IN MODELS WITHOUT MUTATION

In this paper we present a method for estimating population divergence times by maximum likelihood in models without mutation. The maximum-likelihood estimator is compared to a commonly applied estimator based on Wright's F_ST statistic. Simulations suggest that the maximum-likelihood estimator is less biased and has a lower variance than the F_ST-based estimator. The maximum-likelihood estimator provides a statistical framework for the analysis of population history given genetic data. We demonstrate how maximum-likelihood estimates of the branching pattern of divergence of multiple populations may be obtained. We also describe how the method may be applied to test hypotheses such as whether populations have maintained equal population sizes. We illustrate the method by applying it to two previously published sets of human restriction fragment length polymorphism (RFLP) data.