Demographic and genetic collapses in spatially structured populations: insights from a long‐term survey in wild fish metapopulations

Unraveling the relationship between demographic declines and genetic changes over time is of critical importance to predict the persistence of at‐risk populations and to propose efficient conservation plans. This is particularly relevant in spatially structured populations (i.e. metapopulations) in which the spatial arrangement of local populations can modulate both demographic and genetic changes. We used ten‐year demo‐genetic monitoring to test 1) whether demographic declines were associated with genetic diversity declines and 2) whether the spatial structure of a metapopulation can weaken or reinforce these demographic and genetic temporal trends. We continuously surveyed, over time and across their entire range, two metapopulations of an endemic freshwater fish species Leuciscus burdigalensis: one metapopulation that had experienced a recent demographic decline and a second metapopulation that was stable over time. In the declining metapopulation, the number of alleles rapidly decreased, the inbreeding coefficient increased, and a genetic bottleneck emerged over time. In contrast, genetic indices were constant over time in the stable metapopulation. We further show that, in the declining metapopulation, demographic and genetic declines were not homogeneously distributed across the metapopulation. We notably identify one local population situated downstream as a ‘reservoir’ of individuals and genetic variability that dampens both the demographic and genetic declines measured at the metapopulation level. We demonstrate the usefulness of long‐term monitoring that combines both genetic and demographic parameters to understand and predict temporal population fluctuations of at‐risk species living in a metapopulation context.     

Environment,but not migration rate,influences extinction risk in experimental metapopulations

Ecological theory suggests that several demographic factors influence metapopulation extinction risk, including synchrony in population size between subpopulations, metapopulation size and the magnitude of fluctuations in population size. Theoretically, each of these is influenced by the rate of migration between subpopulations. Here we report on an experiment where we manipulated migration rate within metapopulations of the freshwater zooplankton Daphnia magna to examine how migration influenced each of these demographic variables, and subsequent effects on metapopulation extinction. In addition, our experimental procedures introduced unplanned but controlled differences between metapopulations in light intensity, enabling us to examine the relative influences of environmental and demographic factors. We found that increasing migration rate increased subpopulation synchrony. We failed to detect effects of migration on population size and fluctuations in population size at the metapopulation or subpopulation level, however. In contrast, light intensity did not influence synchrony, but was positively correlated with population size and negatively correlated with population fluctuation. Finally, synchrony did not influence time to extinction, while population size and the magnitude of fluctuations did. We conclude that environmental factors had a greater influence on extinction risk than demographic factors, and that metapopulation size and fluctuation were more important to extinction risk than metapopulation synchrony.     

Allee-like effects in metapopulation dynamics

The existences of the Allee effect at the local population level and of the Allee-like effect at the metapopulation level are important for both ecology and conservation. Although there have been a great many papers on the Allee effect, they have mainly referred to only local populations and have not dealt with the relationship between the two. In this paper, we begin with local population dynamics and then construct a model including both local population and metapopulation dynamics. Then we simulate with computer at these two levels. The results indicate that the Allee-like effect in a metapopulation may emerge from the imposed Allee effect at the local population level. This threshold fraction of occupied patches below which the metapopulation goes extinct is seriously affected by the per capita migration rate, the survival rate during migration and the initial population size on the occupied patches. We also find that severe demographic stochasticity may compound the metapopulation extinction risk posed by the Allee effect. These conclusions are helpful for nature conservation, especially for the preservation of rare species.     

Evolution of reproductive effort in a metapopulation with local extinctions and ecological succession

Abstract Using a metapopulation model, we study how local extinctions, limited population life span, and local demographic disequilibrium affect the evolution of the reproductive effort in a species with overlapping generations but no senescence. We show that in a metapopulation with saturation of all sites and an infinite deme maximal life span (no succession), local extinctions simply constitute an additional source of extrinsic mortality. When either the hypothesis of an infinite deme maximal life span or the saturation hypothesis is relaxed, nontrivial predictions arise. in particular, we find interactions between the evolutionarily stable reproductive effort strategy and the demographic dynamics in the metapopulation. We predict that larger reproductive effort may be selected for in habitats of poorer productivity, contrary to what would be predicted in a single population. Also, we predict that higher dispersal rates should favor selection for lower reproductive efforts. However, metapopulation parameters that favor high dispersal rates also favor larger reproductive efforts. Conflicting selection pressures in the metapopulation also allow maintaining evolutionarily stable polymorphism between a low and high reproductive effort for particular trade-offs between survival and fecundity.     

局域种群的Allee效应和集合种群的同步性

从包含Allee效应的局域种群出发,建立了耦合映像格子模型,即集合种群模型.通过分析和计算机模拟表明:(1)当局域种群受到Allee效应强度较大时,集合种群同步灭绝;(2)而当Allee效应强度相对较弱时,通过稳定局域种群动态(减少混沌)使得集合种群发生同步波动,而这种同步波动能够增加集合种群的灭绝风险;(3)斑块间的连接程度对集合种群同步波动的发生有很大的影响,适当的破碎化有利于集合种群的续存.全局迁移和Allee效应结合起来增加了集合种群同步波动的可能,从而增加集合种群的灭绝风险.这些结果对理解同步性的机理、利用同步机理来制定物种保护策略和害虫防治都有重要的意义.     

On valuing patches: estimating contributions to metapopulation growth with reverse-time capture-recapture modelling

Metapopulation ecology has historically been rich in theory, yet analytical approaches for inferring demographic relationships among local populations have been few. We show how reverse-time multi-state capture-recapture models can be used to estimate the importance of local recruitment and interpopulation dispersal to metapopulation growth. We use 'contribution metrics' to infer demographic connectedness among eight local populations of banner-tailed kangaroo rats, to assess their demographic closure, and to investigate sources of variation in these contributions. Using a 7 year dataset, we show that: (i) local populations are relatively independent demographically, and contributions to local population growth via dispersal within the system decline with distance; (ii) growth contributions via local survival and recruitment are greater for adults than juveniles, while contributions involving dispersal are greater for juveniles; (iii) central populations rely more on local recruitment and survival than peripheral populations; (iv) contributions involving dispersal are not clearly related to overall metapopulation density; and (v) estimated contributions from outside the system are unexpectedly large. Our analytical framework can classify metapopulations on a continuum between demographic independence and panmixia, detect hidden population growth contributions, and make inference about other population linkage forms, including rescue effects and source-sink structures. Finally, we discuss differences between demographic and genetic population linkage patterns for our system.     

Evolution of Complex Density-Dependent Dispersal Strategies

The question of how dispersal behavior is adaptive and how it responds to changes in selection pressure is more relevant than ever, as anthropogenic habitat alteration and climate change accelerate around the world. In metapopulation models where local populations are large, and thus local population size is measured in densities, density-dependent dispersal is expected to evolve to a single-threshold strategy, in which individuals stay in patches with local population density smaller than a threshold value and move immediately away from patches with local population density larger than the threshold. Fragmentation tends to convert continuous populations into metapopulations and also to decrease local population sizes. Therefore we analyze a metapopulation model, where each patch can support only a relatively small local population and thus experience demographic stochasticity. We investigated the evolution of density-dependent dispersal, emigration and immigration, in two scenarios: adult and natal dispersal. We show that density-dependent emigration can also evolve to a nonmonotone, “triple-threshold” strategy. This interesting phenomenon results from an interplay between the direct and indirect benefits of dispersal and the costs of dispersal. We also found that, compared to juveniles, dispersing adults may benefit more from density-dependent vs. density-independent dispersal strategies.     

Spatial connectedness imposes local‐ and metapopulation‐level selection on life history through feedbacks on demography

Dispersal evolution impacts the fluxes of individuals and hence, connectivity in metapopulations. Connectivity is therefore decoupled from the structural connectedness of the patches within the spatial network. Because of demographic feedbacks, local selection also drives the evolution of other life history traits. We investigated how different levels of connectedness affect trait evolution in experimental metapopulations of the two‐spotted spider mite. We separated local‐ and metapopulation‐level selection and linked trait divergence to population dynamics. With lower connectedness, an increased starvation resistance and delayed dispersal evolved. Reproductive performance evolved locally by transgenerational plasticity or epigenetic processes. Costs of dispersal, but also changes in local densities and temporal fluctuations herein are found to be putative drivers. In addition to dispersal, demographic traits are able to evolve in response to metapopulation connectedness at both the local and metapopulation level by genetic and/or non‐genetic inheritance. These trait changes impact the persistence of spatially structured populations.     

Determinants of Genetic Structure in a Nonequilibrium Metapopulation of the Plant Silene latifolia

Population genetic differentiation will be influenced by the demographic history of populations, opportunities for migration among neighboring demes and founder effects associated with repeated extinction and recolonization. In natural populations, these factors are expected to interact with each other and their magnitudes will vary depending on the spatial distribution and age structure of local demes. Although each of these effects has been individually identified as important in structuring genetic variance, their relative magnitude is seldom estimated in nature. We conducted a population genetic analysis in a metapopulation of the angiosperm, Silene latifolia, from which we had more than 20 years of data on the spatial distribution, demographic history, and extinction and colonization of demes. We used hierarchical Bayesian methods to disentangle which features of the populations contributed to among population variation in allele frequencies, including the magnitude and direction of their effects. We show that population age, long-term size and degree of connectivity all combine to affect the distribution of genetic variance; small, recently-founded, isolated populations contributed most to increase F _ST in the metapopulation. However, the effects of population size and population age are best understood as being modulated through the effects of connectivity to other extant populations, i.e. F _ST diminishes as populations age, but at a rate that depends how isolated the population is. These spatial and temporal correlates of population structure give insight into how migration, founder effect and within-deme genetic drift have combined to enhance and restrict genetic divergence in a natural metapopulation.     

Fitting parameters of stochastic birth-death models to metapopulation data

Populations that are structured into small local patches are a common feature of ecological and epidemiological systems. Models describing this structure are often referred to as metapopulation models in ecology or household models in epidemiology. Small local populations are subject to demographic stochasticity. Theoretical studies of household disease models without resistant stages (SIS models) have shown that local stochasticity can be ignored for between patch disease transmission if the number of connected patches is large. In that case the distribution of the number of infected individuals per household reaches a stationary distribution described by a birth-death process with a constant immigration term. Here we show how this result, in conjunction with the balancing condition for birth-death processes, provides a framework to estimate demographic parameters from a frequency distribution of local population sizes. The parameter estimation framework is applicable to estimate parameters of disease transmission models as well as metapopulation models.     

Demographic and genetic factors shaping contemporary metapopulation effective size and its empirical estimation in salmonid fish

The preservation of biodiversity requires an understanding of the maintenance of its components, including genetic diversity. Effective population size determines the amount of genetic variance maintained in populations, but its estimation can be complex, especially when populations are interconnected in a metapopulation. Theory predicts that the effective size of a metapopulation (meta-N(e)) can be decreased or increased by population subdivision, but little empirical work has evaluated these predictions. Here, we use neutral genetic markers and simulations to estimate the effective size of a putative metapopulation in Atlantic salmon (Salmo salar). For a weakly structured set of rivers, we find that meta-N(e) is similar to the sum of local deme sizes, whereas higher genetic differentiation among demes dramatically reduces meta-N(e) estimates. Interdemic demographic processes, such as asymmetrical gene flow, may explain this pattern. However, simulations also suggest that unrecognized population subdivision can also introduce downward bias into empirical estimation, emphasizing the importance of identifying the proper scale of distinct demographic and genetic processes. Under natural patterns of connectivity, evolutionary potential may generally be maintained at higher levels than the local population, with implications for conservation given ongoing species declines and habitat fragmentation.     

Metapopulation dynamics: Does it help to have more of the same?

With the interest in conservation biology shifting towards processes from patterns, and to populations from communities, the theory of metapopulation dynamics is replacing the equilibrium theory of island biogeography as the population ecology paradigm in conservation biology. The simplest models of metapopulation dynamics make predictions about the effects of habitat fragmentation - size and isolation of habitat patches - on metapopulation persistence. The simple models may be enriched by considerations of the effects of demographic and environmental stochasticity on the size and extinction probability of local populations. Environmental stochasticity affects populations at two levels: it makes local extinctions more probable, and it also decreases metapopulation persistence time by increasing the correlation of extinction events across populations. Some controversy has arisen over the significance of correlated extinctions, and how they may affect the optimal subdivision of metapopulations to maximize their persistence time.     

Comparison of reintroduction and enhancement effects on metapopulation viability

Metapopulation viability depends upon a balance of extinction and colonization of local habitats by a species. Mechanisms that can affect this balance include physical characteristics related to natural processes (e.g. succession) as well as anthropogenic actions. Plant restorations can help to produce favorable metapopulation dynamics and consequently increase viability; however, to date no studies confirm this is true. Population viability analysis (PVA) allows for the use of empirical data to generate theoretical future projections in the form of median time to extinction and probability of extinction. In turn, PVAs can inform and aid the development of conservation, recovery, and management plans. Pitcher's thistle (Cirsium pitcheri) is a dune endemic that exhibited metapopulation dynamics. We projected viability of three natural and two restored populations with demographic data spanning 15–23 years to determine the degree the addition of reintroduced population affects metapopulation viability. The models were validated by comparing observed and projected abundances and adjusting parameters associated with demographic and environmental stochasticity to improve model performance. Our chosen model correctly predicted yearly population abundance for 60% of the population‐years. Using that model, 50‐year projections showed that the addition of reintroductions increases metapopulation viability. The reintroduction that simulated population performance in early‐successional habitats had the maximum benefit. In situ enhancements of existing populations proved to be equally effective. This study shows that restorations can facilitate and improve metapopulation viability of species dependent on metapopulation dynamics for survival with long‐term persistence of C. pitcheri in Indiana likely to depend on continued active management.     

Climate, social factors and research disturbance influence population dynamics in a declining sociable weaver metapopulation

Population trends are determined by gains through reproduction and immigration, and losses through mortality and emigration. These demographic quantities and resulting population dynamics are affected by different external and internal drivers. We examined how these demographic quantities were affected by weather, research-induced disturbance, local density, colony site and year in a metapopulation of 17 sociable weaver (Philetairus socius) colonies over 17 years of study (4 years for reproduction). Most colonies declined, but at different rates. The four demographic quantities were related to different drivers. Survival strongly varied among years and colonies and was positively related to rainfall and negatively related to extreme temperature (together explaining 30 % of variation) and disturbance (measured as number of captures conducted at a colony; 7 %). There was a trend for a positive relationship between reproduction and rainfall (50 %). Movement was mainly related to local density: individuals were more likely to emigrate from small to large colonies and from colonies that were either well below or above their long-term mean. They were more likely to immigrate into colonies that were nearby, and below their mean size. We then quantified the effects of these relationships on metapopulation dynamics using a multi-site matrix projection model. Rainfall was potentially a strong driver of metapopulation dynamics. In addition, field-work disturbance might have contributed to the decline of this metapopulation but could not explain its full magnitude. Hence, through a combination of analytical methods we were able to obtain information on the main drivers affecting dynamics in a declining metapopulation.     

Host–parasitoid spatial models: the interplay of demographic stochasticity and dynamics

Host–parasitoid metapopulation models have typically been deterministic models formulated with population numbers as a continuous variable. Spatial heterogeneity in local population abundance is a typical (and often essential) feature of these models and means that, even when average population density is high, some patches have small population sizes. In addition, large temporal population fluctuations are characteristic of many of these models, and this also results in periodically small local population sizes. Whenever population abundances are small, demographic stochasticity can become important in several ways. To investigate this problem, we have reformulated a deterministic, host–parasitoid metapopulation as an integer-based model in which encounters between hosts and parasitoids, and the fecundity of individuals are modelled as stochastic processes. This has a number of important consequences: (1) stochastic fluctuations at small population sizes tend to be amplified by the dynamics to cause massive population variability, i.e. the demographic stochasticity has a destabilizing effect; (2) the spatial patterns of local abundance observed in the deterministic counterpart are largely maintained (although the area of ''spatial chaos'' is extended); (3) at small population sizes, dispersal by discrete individuals leads to a smaller fraction of new patches being colonized, so that parasitoids with small dispersal rates have a greater tendency for extinction and higher dispersal rates have a larger competitive advantage; and (4) competing parasitoids that could coexist in the deterministic model due to spatial segregation cannot now coexist for any combination of parameters.     

Incorporating movement into models of grey seal population dynamics

1. One of the most difficult problems in developing spatially explicit models of population dynamics is the validation and parameterization of the movement process. We show how movement models derived from capture-recapture analysis can be improved by incorporating them into a spatially explicit metapopulation model that is fitted to a time series of abundance data. 2. We applied multisite capture-recapture analysis techniques to photo-identification data collected from female grey seals at the four main breeding colonies in the North Sea between 1999 and 2001. The best-fitting movement models were then incorporated into state-space metapopulation models that explicitly accounted for demographic and observational stochasticity. 3. These metapopulation models were fitted to a 20-year time series of pup production data for each colony using a Bayesian approach. The best-fitting model, based on the Akaike Information Criterion (AIC), had only a single movement parameter, whose confidence interval was 82% less than that obtained from the capture-recapture study, but there was some support for a model that included an effect of distance between colonies. 4. The state-space modelling provided improved estimates of other demographic parameters. 5. The incorporation of movement, and the way in which it was modelled, affected both local and regional dynamics. These differences were most evident as colonies approached their carrying capacities, suggesting that our ability to discriminate between models should improve as the length of the grey seal time series increases.     

The applicability of metapopulation theory to large mammals

Metapopulation theory has become a common framework in conservation biology and it is sometimes suggested that a metapopulation approach should be used for management of large mammals. However, it has also been suggested that metapopulation theory would not be applicable to species with long generations compared to those with short ones. In this paper, we review how and on what empirical ground metapopulation terminology has been applied to insects, small mammals and large mammals. The review showed that the metapopulation term sometimes was used for population networks which only fulfilled the broadest possible definition of a metapopulation, i.e. they were subpopulations connected by migrating individuals. We argue that the metapopulation concept should be reserved for networks that also show some kind of metapopulation dynamics. Otherwise it applies to almost all populations and loses its substance. We found much empirical support for metapopulation dynamics in both insects and small mammals, but not in large mammals. A possible reason is the methods used to confirm the existence of metapopulation dynamics. For insects and small mammals, the common approach is to study population turnover through patch occupancy data. Such data is difficult to obtain for large mammals, since longer temporal scales need to be covered to record extinctions and colonizations. Still, many populations of large mammals are exposed to habitat fragmentation and the resulting subpopulations sometimes have high risks of extinction. If there is migration between the subpopulations, the metapopulation framework could provide valuable information on their population dynamics. We suggest that a metapopulation approach can be interesting for populations of large mammals, when there are discrete breeding subpopulations and when these subpopulations have different growth rates and demographic fates. Thus, a comparison of the subpopulations’ demographic fates, rather than subpopulation turnover, can be a feasible alternative for studies of metapopulation dynamics in large mammals.     

Consequences of metapopulation collapse: comparison of genetic attributes between two Allegheny woodrat metapopulations

Disruptions in metapopulation connectivity due to demographic pressure can leave genetically isolated subpopulations susceptible to genetic drift, accumulation of deleterious alleles, and inbreeding depression. Such a scenario may be playing out within Allegheny woodrat (Neotoma magister) metapopulations as a series of synergistic extrinsic pressures have contributed to the rangewide decline of the species over the last 40 years. Our goal was to elucidate the effects of demographic collapse on metapopulation function by using 11 microsatellites markers to quantify differences in patterns of connectivity and genetic diversity between a demographically stable metapopulation and one in severe demographic decline. The demographically diminished metapopulation had lower levels of genetic diversity than the stable metapopulation at all levels evaluated (metapopulation-, subpopulation-, and individual-scales). In contrast to patterns of connectivity observed within the stable metapopulation, peripheral subpopulations in the diminished metapopulation had become completely isolated and were drifting toward genetic fixation, likely as a result of the extirpation of stepping-stone subpopulations. The declining genetic parameters observed within these isolated peripheral subpopulations suggest that inbreeding depression may be contributing significantly to their demographic decline. Allegheny woodrats readily express the genetic consequences of metapopulation decline due to the low effective population sizes of subpopulations and the species’ limited dispersal capacity. Differences in genetic parameters observed between demographically stable and diminished Allegheny woodrat metapopulations emphasize the risks posed to metapopulation function and associated genetic processes introduced with demographic decline.     

Allee threshold and extinction threshold for spatially explicit metapopulation dynamics with Allee effects

In this paper, we investigate a spatially explicit metapopulation model with Allee effects. We refer to the patch occupancy model introduced by Levins (Bull Entomol Soc Am 15:237–240, 1969) as a spatially implicit metapopulation model, i.e., each local patch is either occupied or vacant and a vacant patch can be recolonized by a randomly chosen occupied patch from anywhere in the metapopulation. When we transform the model into a spatially explicit one by using a lattice model, the obtained model becomes theoretically equivalent to a “lattice logistic model” or a “basic contact process”. One of the most popular or standard metapopulation models with Allee effects, developed by Amarasekare (Am Nat 152:298–302, 1998), supposes that those effects are introduced formally by means of a logistic equation. However, it is easier to understand the ecological meaning of associating Allee effects with this model if we suppose that only the logistic colonization term directly suffers from Allee effects. The resulting model is also well defined, and therefore we can naturally examine it by Monte Carlo simulation and by doublet and triplet decoupling approximation. We then obtain the following specific features of one-dimensional lattice space: (1) the metapopulation as a whole does not have an Allee threshold for initial population size even when each local population follows the Allee effects; and (2) a metapopulation goes extinct when the extinction rate of a local population is lower than that in the spatially implicit model. The real ecological metapopulation lies between two extremes: completely mixing interactions between patches on the one hand and, on the other, nearest neighboring interactions with only two nearest neighbors. Thus, it is important to identify the metapopulation structure when we consider the problems of invasion species such as establishment or the speed of expansion.