Contrasting patterns of genetic diversity at three different genetic markers in a marine mammal metapopulation

Many studies use genetic markers to explore population structure and variability within species. However, only a minority use more than one type of marker and, despite increasing evidence of a link between heterozygosity and individual fitness, few ask whether diversity correlates with population trajectory. To address these issues, we analyzed data from the Steller’s sea lion, Eumetiopias jubatus, where three stocks are distributed over a vast geographical range and where both genetic samples and detailed demographic data have been collected from many diverse breeding colonies. To previously published mitochondrial DNA(mtDNA) and microsatellite data sets,we have added new data for amplified fragment length polymorphism (AFLP) markers, comprising 238 loci scored in 285 sea lions sampled from 23 natal rookeries. Genotypic diversity was low relative to most vertebrates, with only 37 loci (15.5%) being polymorphic. Moreover, contrasting geographical patterns of genetic diversity were found at the three markers, with Nei’s gene diversity tending to be higher for AFLPs and microsatellites in rookeries of the western and Asian stocks, while the highest mtDNA values were found in the eastern stock. Overall, and despite strongly contrasting demographic histories, after applying phylogenetic correction we found little correlation between genetic diversity and either colony size or demography. In contrast, we were able to show a highly significant positive relationship between AFLP diversity and current population size across a range of pinniped species, even though equivalent analyses did not reveal significant trends for either microsatellites or mtDNA.

Metapopulation extinction caused by mutation accumulation

Theory suggests that the risk of extinction by mutation accumulation can be comparable to that by environmental stochasticity for an isolated population smaller than a few thousand individuals. Here we show that metapopulation structure, habitat loss or fragmentation, and environmental stochasticity can be expected to greatly accelerate the accumulation of mildly deleterious mutations, lowering the genetic effective size to such a degree that even large metapopulations may be at risk of extinction. Because of mutation accumulation, viable metapopulations may need to be far larger and better connected than would be required under just stochastic demography.

Genetic drift and within-host metapopulation dynamics of HIV-1 infection

Most HIV replication occurs in solid lymphoid tissue, which has prominent architecture at the histological level, which separates groups of productively infected CD4+ cells. Nevertheless, current population models of HIV assume panmixis within lymphoid tissue. We present a simple “metapopulation” model of HIV replication, where the population of infected cells is comprised of a large number of small populations, each of which is established by a few founder viruses and undergoes turnover. To test this model, we analyzed viral genetic variation of infected cell subpopulations within the spleen and demonstrated the action of founder effects as well as significant variation in the extent of genetic differentiation between subpopulations among patients. The combination of founder effects and subpopulation turnover can result in an effective population size much lower than the actual population size and may contribute to the importance of genetic drift in HIV evolution despite a large number of infected cells.

Green turtle somatic growth dynamics in a spatially disjunct Great Barrier Reef metapopulation

The growth dynamics of green sea turtles resident in four separate foraging grounds of the southern Great Barrier Reef genetic stock were assessed using a nonparametric regression modeling approach. Juveniles recruit to these grounds at the same size, but grow at foraging-ground-dependent rates that result in significant differences in expected size- or age-at-maturity. Mean age-at-maturity was estimated to vary from 25-50 years depending on the ground. This stock comprises mainly the same mtDNA haplotype, so geographic variability might be due to local environmental conditions rather than genetic factors, although the variability was not a function of latitudinal variation in environmental conditions or whether the food stock was seagrass or algae. Temporal variability in growth rates was evident in response to local environmental stochasticity, so geographic variability might be due to local food stock dynamics. Despite such variability, the expected size-specific growth rate function at all grounds displayed a similar nonmonotonic growth pattern with a juvenile growth spurt at 60-70 cm curved carapace length, (CCL) or 15-20 years of age. Sex-specific growth differences were also evident with females tending to grow faster than similar-sized males after the Juvenile growth spurt. It is clear that slow sex-specific growth displaying both spatial and temporal variability and a juvenile growth spurt are distinct growth behaviors of green turtles from this stock.

A stochastic metapopulation model accounting for habitat dynamics

A stochastic metapopulation model accounting for habitat dynamics is presented. This is the stochastic SIS logistic model with the novel aspect that it incorporates varying carrying capacity. We present results of Kurtz and Barbour, that provide deterministic and diffusion approximations for a wide class of stochastic models, in a form that most easily allows their direct application to population models. These results are used to show that a suitably scaled version of the metapopulation model converges, uniformly in probability over finite time intervals, to a deterministic model previously studied in the ecological literature. Additionally, they allow us to establish a bivariate normal approximation to the quasi-stationary distribution of the process. This allows us to consider the effects of habitat dynamics on metapopulation modelling through a comparison with the stochastic SIS logistic model and provides an effective means for modelling metapopulations inhabiting dynamic landscapes.

The role of habitat disturbance and recovery in metapopulation persistence

Classical metapopulation theory assumes a static landscape. However, empirical evidence indicates many metapopulations are driven by habitat succession and disturbance. We develop a stochastic metapopulation model, incorporating habitat disturbance and recovery, coupled with patch colonization and extinction, to investigate the effect of habitat dynamics on persistence. We discover that habitat dynamics play a fundamental role in metapopulation dynamics. The mean number of suitable habitat patches is not adequate for characterizing the dynamics of the metapopulation. For a fixed mean number of suitable patches, we discover that the details of how disturbance affects patches and how patches recover influences metapopulation dynamics in a fundamental way. Moreover, metapopulation persistence is dependent not only oil the average lifetime of a patch, but also on the variance in patch lifetime and the synchrony in patch dynamics that results from disturbance. Finally, there is an interaction between the habitat and metapopulation dynamics, for instance declining metapopulations react differently to habitat dynamics than expanding metapopulations. We close, emphasizing the importance of using performance measures appropriate to stochastic systems when evaluating their behavior, such as the probability distribution of the state of the. metapopulation, conditional on it being extant (i.e., the quasistationary distribution).

An R package for simulating metapopulation dynamics and range expansion under environmental change

The metapopulation paradigm is central in ecology and conservation biology to understand the dynamics of spatially-structured populations in fragmented landscapes. Metapopulations are often studied using simulation modelling, and there is an increasing demand of user-friendly software tools to simulate metapopulation responses to environmental change. Here we describe the MetaLandSim R package, mwhich integrates ideas from metapopulation and graph theories to simulate the dynamics of real and virtual metapopulations. The package offers tools to (i) estimate metapopulation parameters from empirical data, (ii) to predict variation in patch occupancy over time in static and dynamic landscapes, either real or virtual, and (iii) to quantify the patterns and speed of metapopulation expansion into empty landscapes. MetaLandSim thus provides detailed information on metapopulation processes, which can be easily combined with land use and climate change scenarios to predict metapopulation dynamics and range expansion for a variety of taxa and ecological systems.

Marine protected areas for spatially structured exploited stocks

Many harvested marine and terrestrial populations have segments of their range protected in areas free from exploitation. Reasons for areas being protected from harvesting include conservation, tourism, research, protection of breeding grounds, stock recovery, harvest regulation, or habitat that is uneconomical to exploit. In this paper we consider the problem of optimally exploiting a single species local population that is connected by dispersing larvae to an unharvested local population. We define a spatially-explicit population dynamics model and apply dynamic optimization techniques to determine policies for harvesting the exploited patch. We then consider how reservation affects yield and spawning stock abundance when compared to policies that have not recognised the spatial structure of the metapopulation. Comparisons of harvest strategies between an exploited metapopulation with and without a harvest refuge are also made. Results show that in a 2 local population metapopulation with unidirectional larval transfer, the optimal exploitation of the harvested population should be conducted as if it were independent of the reserved population. Numerical examples suggest that relative source populations should be exploited if the objective is to maximise spawning stock abundance within a harvested metapopulation that includes a protected local population. However, this strategy can markedly reduce yield over a sink harvested reserve system and may require strict regulation for conservation goals to be realised. If exchange rates are high, results indicate that spawning stock abundance can be less in a reserve system than in a fully exploited metapopulation. In order to maximise economic gain in the reserve system, results indicate that relative sink populations should be harvested. Depending on transfer levels, loss in harvest through reservation can be minimal, and is likely to be compensated by the potential environmental and economic benefits of the reserve.

Optimal release strategies for biological control agents: an application of stochastic dynamic programming to population management

1. Establishing biological control agents in the field is a major step in any classical biocontrol programme, yet there are few general guidelines to help the practitioner decide what factors might enhance the establishment of such agents. 2. A stochastic dynamic programming (SDP) approach, linked to a metapopulation model, was used to find optimal release strategies (number and size of releases), given constraints on time and the number of biocontrol agents available. By modelling within a decision-making framework we derived rules of thumb that will enable biocontrol workers to choose between management options, depending on the current state of the system. 3. When there are few well-established sites, making a few large releases is the optimal strategy. For other states of the system, the optimal strategy ranges from a few large releases, through a mixed strategy (a variety of release sizes), to many small releases, as the probability of establishment of smaller inocula increases. 4. Given that the probability of establishment is rarely a known entity, we also strongly recommend a mixed strategy in the early stages of a release programme, to accelerate learning and improve the chances of finding the optimal approach.

A landscape-scale test of the predictive ability of a spatially explicit model for population viability analysis

1. Although population viability analysis (PVA) is widely employed, forecasts from PVA models are rarely tested. This study in a fragmented forest in southern Australia contrasted field data on patch occupancy and abundance for the arboreal marsupial greater glider Petauroides volans with predictions from a generic spatially explicit PVA model. This work represents one of the first landscape-scale tests of its type. 2. Initially we contrasted field data from a set of eucalypt forest patches totalling 437 ha with a naive null model in which forecasts of patch occupancy were made, assuming no fragmentation effects and based simply on remnant area and measured densities derived from nearby unfragmented forest. The naive null model predicted an average total of approximately 170 greater gliders, considerably greater than the true count (n = 81). 3. Congruence was examined between field data and predictions from PVA under several metapopulation modelling scenarios. The metapopulation models performed better than the naive null model. Logistic regression showed highly significant positive relationships between predicted and actual patch occupancy for the four scenarios (P = 0.001-0.006). When the model-derived probability of patch occupancy was high (0.50-0.75, 0.75-1.00), there was greater congruence between actual patch occupancy and the predicted probability of occupancy. 4. For many patches, probability distribution functions indicated that model predictions for animal abundance in a given patch were not outside those expected by chance. However, for some patches the model either substantially over-predicted or under-predicted actual abundance. Some important processes, such as inter-patch dispersal, that influence the distribution and abundance of the greater glider may not have been adequately modelled. 5. Additional landscape-scale tests of PVA models, on a wider range of species, are required to assess further predictions made using these tools. This will help determine those taxa for which predictions are and are not accurate and give insights for improving models for applied conservation management.

Detecting environmental impacts on metapopulations of mound spring invertebrates - Assessing an incidence function model

We use a stochastic patch occupancy model of invertebrates in the Mound Springs ecosystem of South Australia to assess the ability of incidence function models to detect environmental impacts on metapopulations. We assume that the probability of colonisation decreases with increasing isolation and the probability of extinction is constant across spring vents. We run the models to quasi-equilibrium, and then impose an impact by increasing the local extinction probability. We sample the output at various times pre- and postimpact, and examine the probability of detecting a significant change in population parameters. The incidence function model approach turns out to have little power to detect environmental impacts on metapopulations with small numbers of patches. (C) 2001 Elsevier Science Ltd. All rights reserved.

Disease, habitat fragmentation and conservation

Habitat loss and the resultant fragmentation of remaining habitat is the primary cause of loss of biological diversity. How do these processes affect the dynamics of parasites and pathogens? Hess has provided some important insights into this problem using metapopulation models for pathogens that exhibit 'S-I' dynamics; for example, pathogens such as rabies in which the host population may be divided into susceptible and infected individuals. A major assumption of Hess's models is that infected patches become extinct, rather than recovering and becoming resistant to future infections. In this paper, we build upon this framework in two different ways: first, we examine the consequences of including patches that are resistant to infection; second, we examine the consequences of including a second species of host that can act as a reservoir for the pathogen. Both of these effects are likely to be important from a conservation perspective. The results of both sets of analysis indicate that the benefits of corridors and other connections that allow species to disperse through the landscape far outweigh the possible risks of increased pathogen transmission. Even in the commonest case, where harmful pathogens are maintained by a common reservoir host, increased landscape connectance still allows greater coexistence and persistence of a threatened or endangered host.

The use of stochastic dynamic programming in optimal landscape reconstruction for metapopulations

A decision theory framework can be a powerful technique to derive optimal management decisions for endangered species. We built a spatially realistic stochastic metapopulation model for the Mount Lofty Ranges Southern Emu-wren (Stipiturus malachurus intermedius), a critically endangered Australian bird. Using diserete-time Markov,chains to describe the dynamics of a metapopulation and stochastic dynamic programming (SDP) to find optimal solutions, we evaluated the following different management decisions: enlarging existing patches, linking patches via corridors, and creating a new patch. This is the first application of SDP to optimal landscape reconstruction and one of the few times that landscape reconstruction dynamics have been integrated with population dynamics. SDP is a powerful tool that has advantages over standard Monte Carlo simulation methods because it can give the exact optimal strategy for every landscape configuration (combination of patch areas and presence of corridors) and pattern of metapopulation occupancy, as well as a trajectory of strategies. It is useful when a sequence of management actions can be performed over a given time horizon, as is the case for many endangered species recovery programs, where only fixed amounts of resources are available in each time step. However, it is generally limited by computational constraints to rather small networks of patches. The model shows that optimal metapopulation, management decisions depend greatly on the current state of the metapopulation,. and there is no strategy that is universally the best. The extinction probability over 30 yr for the optimal state-dependent management actions is 50-80% better than no management, whereas the best fixed state-independent sets of strategies are only 30% better than no management. This highlights the advantages of using a decision theory tool to investigate conservation strategies for metapopulations. It is clear from these results that the sequence of management actions is critical, and this can only be effectively derived from stochastic dynamic programming. The model illustrates the underlying difficulty in determining simple rules of thumb for the sequence of management actions for a metapopulation. This use of a decision theory framework extends the capacity of population viability analysis (PVA) to manage threatened species.

SPAnDX: a process-based population dynamics model to explore management and climate change impacts on an invasive alien plant, Acacia nilotica

This paper describes a process-based metapopulation dynamics and phenology model of prickly acacia, Acacia nilotica, an invasive alien species in Australia. The model, SPAnDX, describes the interactions between riparian and upland sub-populations of A. nilotica within livestock paddocks, including the effects of extrinsic factors such as temperature, soil moisture availability and atmospheric concentrations of carbon dioxide. The model includes the effects of management events such as changing the livestock species or stocking rate, applying fire, and herbicide application. The predicted population behaviour of A. nilotica was sensitive to climate. Using 35 years daily weather datasets for five representative sites spanning the range of conditions that A. nilotica is found in Australia, the model predicted biomass levels that closely accord with expected values at each site. SPAnDX can be used as a decision-support tool in integrated weed management, and to explore the sensitivity of cultural management practices to climate change throughout the range of A. nilotica. The cohort-based DYMEX modelling package used to build and run SPAnDX provided several advantages over more traditional population modelling approaches (e.g. an appropriate specific formalism (discrete time, cohort-based, process-oriented), user-friendly graphical environment, extensible library of reusable components, and useful and flexible input/output support framework). (C) 2003 Published by Elsevier Science B.V.