Regional and local scale metapopulation dynamics in the interaction between Callosobruchus maculatus and Anisopteromalus calandrae

Habitat structure increases the persistence of many extinction‐prone resource–consumer interactions. Metapopulation theory is one of the leading approaches currently used to explain why local, ephemeral populations persist at a regional scale. Central to the metapopulation concept is the amount of dispersal occurring between patches, too much or too little can result in regional extinction. In this study, the role of dispersal on the metapopulation dynamics of an over‐exploitative host–parasitoid interaction is assessed. In the absence of the parasitoid the highly vagile bruchid, Callosobruchus maculatus, can maintain a similar population size regardless of the permeability of the inter‐patch matrix and exhibits strong negative density‐dependence. After the introduction of the parasitoid the size of the bruchid population decreases with a corresponding increase in the occurrence of empty patches. In this case, limiting the dispersal of both species decouples the interaction to a greater extent and results in larger regional bruchid populations. Given the disparity between the dispersal rates of the two species, it is proposed that the more dispersive host benefits from the reduction in landscape permeability by increasing the opportunity to colonise empty patches and rescue extinction prone populations. Associated with the introduction of the parasitoid is a shift in the strength of density‐dependence as the population moves from bottom–up towards top–down regulation. The importance of local and regional scale measurements is apparent when the role of individual patches on regional dynamics is considered. By only taking regional dynamics into account the importance of dispersal regime on local dynamics is overlooked. Similarly, when local dynamics were examined, patches were found to have different influences on regional dynamics depending on dispersal regime and patch location.     

Invasion and persistence of infectious agents in fragmented host populations

One of the important questions in understanding infectious diseases and their prevention and control is how infectious agents can invade and become endemic in a host population. A ubiquitous feature of natural populations is that they are spatially fragmented, resulting in relatively homogeneous local populations inhabiting patches connected by the migration of hosts. Such fragmented population structures are studied extensively with metapopulation models. Being able to define and calculate an indicator for the success of invasion and persistence of an infectious agent is essential for obtaining general qualitative insights into infection dynamics, for the comparison of prevention and control scenarios, and for quantitative insights into specific systems. For homogeneous populations, the basic reproduction ratio R(0) plays this role. For metapopulations, defining such an 'invasion indicator' is not straightforward. Some indicators have been defined for specific situations, e.g., the household reproduction number R*. However, these existing indicators often fail to account for host demography and especially host migration. Here we show how to calculate a more broadly applicable indicator R(m) for the invasion and persistence of infectious agents in a host metapopulation of equally connected patches, for a wide range of possible epidemiological models. A strong feature of our method is that it explicitly accounts for host demography and host migration. Using a simple compartmental system as an example, we illustrate how R(m) can be calculated and expressed in terms of the key determinants of epidemiological dynamics.     

Host–parasitoid persistence over variable spatio‐temporally susceptible habitats: bottom–up effects of ephemeral resources

We experimentally and theoretically investigated the persistence of hosts and parasitoids interacting in a metapopulation structure consisting of ephemeral local patches (MELPs). We used a host–parasitoid system consisting of necrophagous Diptera species and their pupal parasitoids. The basal resources used by the host species were assumed to be ephemeral, supporting only one generation of individuals before completely disappearing from the environment. We experimentally measured the host–parasitoid persistence and the effects of local demographic processes in two scenarios: 1) constant occurrence of basal resources at a single site (no dispersion or colonization of other sites) and 2) variable occurrence of basal resources between two sites (colonization of a new patch requiring species dispersal). The experimental setup and findings were then formalized into a mathematical model describing the interaction dynamics in a MELP structure. We evaluated the contribution of several factors to the host–parasitoid coexistence, such as resource allocation probability (probability of resource appearance in a site), variation in resource size and number of sites available to receive resources in the MELP. We found that demographic fluctuations and environmental stochasticity affected the density of migrants, patch habitat connectivity, persistence and spatial distribution of interacting species.     

Population structure of a large blue butterfly and its specialist parasitoid in a fragmented landscape

Habitat fragmentation may interrupt trophic interactions if herbivores and their specific parasitoids respond differently to decreasing connectivity of populations. Theoretical models predict that species at higher trophic levels are more negatively affected by isolation than lower trophic level species. By combining ecological data with genetic information from microsatellite markers we tested this hypothesis on the butterfly Maculinea nausithous and its specialist hymenopteran parasitoid Neotypus melanocephalus. We assessed the susceptibility of both species to habitat fragmentation by measuring population density, rate of parasitism, overall genetic differentiation (theta(ST)) and allelic richness in a large metapopulation. We also simulated the dynamics of genetic differentiation among local populations to asses the relative effects of migration rate, population size, and haplodiploid (parasitoid) and diploid (host) inheritance on metapopulation persistence. We show that parasitism by N. melanocephalus is less frequent at larger distances to the nearest neighbouring population of M. nausithous hosts, but that host density itself is not affected by isolation. Allelic richness was independent of isolation, but the mean genetic differentiation among local parasitoid populations increased with the distance between these populations. Overall, genetic differentiation in the parasitoid wasp was much greater than in the butterfly host and our simulations indicate that this difference is due to a combination of haplodiploidy and small local population sizes. Our results thus support the hypothesis that Neotypus parasitoid wasps are more sensitive to habitat fragmentation than their Maculinea butterfly hosts.     

Metapopulation dynamics in an aphid-parasitoid system

Metapopulation theory makes a number of predictions concerning the effects of dispersal on the persistence of predator-prey or host-parasitoid systems. While the stabilising effects of dispersal have been shown in a number of laboratory studies, evidence from field studies remains scarce due to a lack of suitable model systems. I describe a host-parasitoid system that shows a classical metapopulation structure with frequent extinctions and colonisations consisting of the aphidiid Lysiphlebus hirticornisand the aphid Metopeurum fuscoviride. Both the parasitoid and the aphid are specialists on their respective hosts. I followed the dynamics of host and parasitoid on individually marked tansy (Tanacetum vulgare) plants, the host of M. fuscoviride. Dynamics of host and parasitoid populations were characterized by frequent extinctions and colonisations. Mean longevity of aphid colonies was only 3.1 weeks. Parasitism by L. hirticorniswas a main cause of extinction for the aphid as rates of parasitism often reached 100%, in particular towards the end of the field season. Patchiness in this system occurs at two spatial scales. Aphid colonies form on single tansy ramets = shoots but movements of aphid individuals among ramets within a particular tansy genet are frequent. Because aphids can persist on a genet for a large numer of generations, it is argued that local populations form on genets rather than ramets. The number of host and parasitoid extinctions described in this study exceeds the number of extinctions usually observed in field studies of host-parasitoid metapopulations. It is suggested that aphid-parasitoid systems such as the one studied in this paper may be good models to test the predictions of metapopulation theory.     

Dances with fire: Tracking metapopulation dynamics of<Emphasis Type="Italic">polygonella basiramia</Emphasis> in Florida scrub (USA)

The regional persistence of species subject to local population colonization and extinction necessarily depends on how landscape features and disturbance affect metapopulation dynamics. Here, we characterize the metapopulation structure and short-term dynamics ofPolygonella basiramia. This rare, short-lived perennial herb is endemic to Florida scrublands and lacks a seed bank. Fires create the open sand gaps within a shrub matrix that support this species but also kill established plants. Thus, persistence depends on frequent colonization of unoccupied gaps. We are monitoring population dynamics within and among 1204 gaps distributed among 19 shrub patches. Considerable subpopulation turnover is evident at the gap level with rates of gap extinction exceeding rates of colonization in the first year. Whether declines in overall abundance continue is likely to depend on patterns of disturbance and regional stochasticity in this dynamic landscape.Polygonella is more likely to occupy larger and less isolated gaps, demonstrating that landscape features and disturbance strongly affect metapopulation dynamics. BecausePolygonella basiramia displays characteristics, occupancy patterns, and turnover dynamics consistent with metapopulation theory, it represents a model system for studying plant metapopulations.     

Temporal/spatial structure and the dynamical property of laboratory host-parasitoid systems

The effects of spatial structure in terms of local capacity, or the maximum number of larvae surviving competition at resource patches, and temporal structure in terms of the period vulnerable to parasitoid attack in host populations on the persistence of host-parasitoid systems were quantitatively evaluated by laboratory experiments and well-parameterized model analyses. One of two bruchid beetles,Callosobruchus maculatus andC. phaseoli, were used as a host with Heterospilus prosopidis used as the parasitoid.C. maculatus, in which few larvae survive competition to become adults in each bean, andC. phaseoli, in which many larvae become adults in each bean, along with two kinds of beans, the mung and the azuki, were combined to construct four (2×2) resource-herbivorous host-parasitoid systems that differed in local capacity and vulnerable period. The mung-C. maculatus system with the parasitoid was the most persistent, i.e., took the longest time for extinction of either the host or parasitoid to occur. Since this resource-herbivorous host combination exhibited the lowest local capacity and the shortest vulnerable period, these two conditions possibly promoted the persistence of the system. A model incorporating the host population structure supported the observed persistence. Furthermore, the possible contribution of the timing of density-dependent competition of the host on the host-parasitoid persistence is predicted.     

Does colonization asymmetry matter in metapopulations?

Despite the considerable evidence showing that dispersal between habitat patches is often asymmetric, most of the metapopulation models assume symmetric dispersal. In this paper, we develop a Monte Carlo simulation model to quantify the effect of asymmetric dispersal on metapopulation persistence. Our results suggest that metapopulation extinctions are more likely when dispersal is asymmetric. Metapopulation viability in systems with symmetric dispersal mirrors results from a mean field approximation, where the system persists if the expected per patch colonization probability exceeds the expected per patch local extinction rate. For asymmetric cases, the mean field approximation underestimates the number of patches necessary for maintaining population persistence. If we use a model assuming symmetric dispersal when dispersal is actually asymmetric, the estimation of metapopulation persistence is wrong in more than 50% of the cases. Metapopulation viability depends on patch connectivity in symmetric systems, whereas in the asymmetric case the number of patches is more important. These results have important implications for managing spatially structured populations, when asymmetric dispersal may occur. Future metapopulation models should account for asymmetric dispersal, while empirical work is needed to quantify the patterns and the consequences of asymmetric dispersal in natural metapopulations.     

Host-parasitoid extinction and colonization in a fragmented prairie landscape

Few field studies of natural populations have examined the factors influencing local extinctions and colonization of empty habitat patches for a prey species and its predator. In this study, I carried out a census of planthopper (Prokelisia crocea; Hemiptera: Delphacidae) and egg parasitoid (Anagrus columbi; Hymenoptera: Mymaridae) incidence and densities in 147 host-plant patches (Spartina pectinata; Poaceae) over seven planthopper generations in a tall-grass prairie landscape. For both species, the likelihood of going extinct in a patch was related to a number of patch-specific variables: density, temporal variability in density, proportion of hosts parasitized (planthopper only), host-plant density, patch size, patch isolation, and composition of the surrounding matrix. Colonization likelihood was related only to the physical attributes of the patch. There was high patch turnover in this prairie landscape. On average, planthoppers went extinct in 23% of the patches and A. columbi went extinct in 51% of the patches in each generation. For the planthopper, extinction likelihood increased with a decrease in patch size and the proportion of the matrix composed of mudflat. Parasitism of eggs had no effect on the extinction likelihood of local P. crocea populations, suggesting that A. columbi may not play a major role in the patch dynamics of its host. The likelihood of extinction for A. columbi was dependent on factors that spanned three trophic levels. An increase in plant density, decrease in host density and decrease in parasitoid density all increased the likelihood of A. columbi extinction within a patch. The dependency on multiple trophic levels may explain the higher extinction risk for the parasitoid than its host. A. columbi extinction was also affected by the matrix habitat surrounding the patch—the effect was the opposite of that for P. crocea. Finally, vacant patches were colonized at rates of 53% and 34% per generation for the planthopper and parasitoid, respectively. For both species, colonization probabilities decreased with an increase in patch isolation. High host densities in a patch also favored high rates of colonization by A. columbi. I discuss how anthropogenic changes to the prairie landscape can affect the metapopulation dynamics and persistence time of this host-parasitoid interaction.     

Effects of colonization asymmetries on metapopulation persistence

Ocean currents, prevailing winds, and the hierarchical structures of river networks are known to create asymmetries in re-colonization between habitat patches. The impacts of such asymmetries on metapopulation persistence are seldom considered, especially rarely in theoretical studies. Considering three classical models (the island, the stepping stone and the distance-dependent model), we explore how metapopulation persistence is affected by (i) asymmetry in dispersal strength, in which the colonization rate between two patches differs in direction, and (ii) asymmetry in connectivity, in which the overall colonization pattern displays asymmetry (circulating or dendritic networks). Viability can be drastically reduced when directional bias in dispersal strength is higher than 25%. Re-colonization patterns that allow for strong local connectivity provide the highest persistence compared to systems that allow circulation. Finally, asymmetry has relatively weak effects when metapopulations maintain strong general connectivity.     

Local and spatial dynamics of a host–parasitoid system in a field experiment

Local extinction and colonisation rates are key factors in host–parasitoid metapopulation theory, but experimental evidence from the field is scarce. We studied the host–parasitoid system consisting of the aphid Metopeurum fuscoviride and its specialist parasitoid Lysiphlebus hirticornis. This system is characterised by a patchy distribution of the host plants (Tanacetum vulgare) and by frequent extinctions of local aphid populations. In a first field experiment, we found that the presence of the parasitoid increases the likelihood of extinction of local host populations (=all aphids living on one plant). In a second field experiment, we manipulated the distance between local host populations. Parasitoid colonisation rate strongly decreased with increasing distance between local host populations. Thus, our results show the importance of parasitoids for local host populations extinction and of distance between local host populations for parasitoid colonisation rate, suggesting the importance of spatial processes for host–parasitoid systems in the field.     

The persistence of vertically transmitted fungi in grass metapopulations

Theory predicts that (i) vertical transmission of parasites (i.e. when they are passed directly from a host to its offspring) selects for benign association with the host and that (ii) vertically transmitted parasites that lower their hosts' fitness cannot persist if they are not able to infect horizontally (i.e. contagiously) other host individuals in the population. In this paper, we develop a mathematical model to examine whether mutualism is a prerequisite for persistence of exclusively vertically transmitted (from maternal plant to offspring via seeds) fungal endophytes in structured grass metapopulations. Interestingly, endophyte survival does not require plant mutualism, even in a metapopulation consisting of qualitatively identical patches, if vertical transmission of the fungus is perfect, i.e. if all established seedlings in offspring of the endophyte-infected plant are infected. In more realistic situations, when the metapopulation consists of qualitatively different patches, endophyte-infected plants may persist at the metapopulation level even if the vertical transmission is imperfect (due to hyphae inviability or failure to grow into all seeds) and the endophyte decreases the host grass fitness in certain environments. These results have biological importance because they (i) question the requirement of a mutualistic nature in exclusively vertically transmitted symbionts and (ii) emphasize the importance of habitat diversity in relation to symbiont success in vertical transmission.     

How patch configuration affects the impact of disturbances on metapopulation persistence

Disturbances affect metapopulations directly through reductions in population size and indirectly through habitat modification. We consider how metapopulation persistence is affected by different disturbance regimes and the way in which disturbances spread, when metapopulations are compact or elongated, using a stochastic spatially explicit model which includes metapopulation and habitat dynamics. We discover that the risk of population extinction is larger for spatially aggregated disturbances than for spatially random disturbances. By changing the spatial configuration of the patches in the system--leading to different proportions of edge and interior patches--we demonstrate that the probability of metapopulation extinction is smaller when the metapopulation is more compact. Both of these results become more pronounced when colonization connectivity decreases. Our results have important management implication as edge patches, which are invariably considered to be less important, may play an important role as disturbance refugia.     

Metapopulation Persistence in Random Fragmented Landscapes

Habitat destruction and land use change are making the world in which natural populations live increasingly fragmented, often leading to local extinctions. Although local populations might undergo extinction, a metapopulation may still be viable as long as patches of suitable habitat are connected by dispersal, so that empty patches can be recolonized. Thus far, metapopulations models have either taken a mean-field approach, or have modeled empirically-based, realistic landscapes. Here we show that an intermediate level of complexity between these two extremes is to consider random landscapes, in which the patches of suitable habitat are randomly arranged in an area (or volume). Using methods borrowed from the mathematics of Random Geometric Graphs and Euclidean Random Matrices, we derive a simple, analytic criterion for the persistence of the metapopulation in random fragmented landscapes. Our results show how the density of patches, the variability in their value, the shape of the dispersal kernel, and the dimensionality of the landscape all contribute to determining the fate of the metapopulation. Using this framework, we derive sufficient conditions for the population to be spatially localized, such that spatially confined clusters of patches act as a source of dispersal for the whole landscape. Finally, we show that a regular arrangement of the patches is always detrimental for persistence, compared to the random arrangement of the patches. Given the strong parallel between metapopulation models and contact processes, our results are also applicable to models of disease spread on spatial networks.     

Spatial ecology of host–parasitoid interactions: a threatened butterfly and its specialised parasitoid

Habitat conservation for threatened temperate insect species is often guided by one of two paradigms: a metapopulation approach focusing on patch area, isolation and number; or a habitat approach focusing on maintaining high quality habitat for the focal species. Recent research has identified the additive and interacting importance of both approaches for maintaining populations of threatened butterflies. For specialised host-parasitoid interactions, understanding the consequences of habitat characteristics for the interacting species is important, because (1) specialised parasitoids are particularly vulnerable to the consequences of fragmentation, and (2) altered interaction frequencies resulting from changes to habitat management or the spatial configuration of habitat are likely to have consequences for host dynamics. The spatial ecology of Cotesia bignellii, a specialist parasitoid of the threatened butterfly Euphydryas aurinia, was investigated at two spatial scales: within habitat patches (at the scale of individual aggregations of larvae, or ‘webs’) and among habitat patches (the scale of local populations). Parasitism rates were investigated in relation to larval web size, vegetation sward height and host density. Within patches, the probability of a larval webs being parasitized increased significantly with increasing number of larvae in the web, and parasitism rates increased significantly with increasing web isolation. The proportion of webs parasitized was significantly and negatively correlated with cluster density. Among habitat patches the proportion of parasitized webs decreased as cluster density increased. Clusters with a high proportion of larval webs parasitized tended to have lower parasitism rates per larval web. These results support the call for relatively large and continuous habitat patches to maintain stable parasitoid and host populations. Conservation efforts directed towards maintenance of high host plant density could allow E. aurinia to reduce parasitism risk, while providing C. bignellii with sufficient larval webs to allow population persistence.     

Metapopulation viability analysis of the bog fritillary butterfly using RAMAS/GIS

Many species living in man-shaped landscapes are restricted to small natural habitat patches and form metapopulations; predicting their future is a central issue in applied ecology. We examined the viability of the bog fritillary butterfly Proclossiana eunomia Esper, a specialist glacial relict species, in a highly fragmented landscape (<1% of suitable habitat in 10 km²), by way of population viability analysis. We used comprehensive data from a long-term study in which a patchy population was monitored during ten consecutive years to parameterise a spatially structured metapopulation model using commercially available platform RAMAS/GIS 3.0. Population growth rate was density-dependent and modulated by various climatic variables acting on different developmental stages of the butterfly. Density dependence was probably related to larval parasitism by a specific parasitoid. Population size was negatively affected by an increase in the mean temperature. Dispersal was modelled as the observed proportion of movements between patches, taking into account the probability of emigration out of a given patch. Our model provided results close to the picture of the system drawn from the field data and was considered as useful in making predictions about the metapopulation. Demographic parameters proved to have a far higher impact on metapopulation persistence than dispersal or correlation of local dynamics. Scenarios simulating both global warming and management of habitat patches by rustic herbivore grazing indicated a decrease in the viability of the metapopulation. Our results prompted the regional nature conservation agency to modify the planned management regime. We urge conservation biologists to use structured population models including local population dynamics for viability analysis targeted to such threatened metapopulations in highly fragmented landscapes.     

Population dynamics and sex ratio of a parasitoid altered by fungal-infected diet of host butterfly

Variation of host quality affects population dynamics of parasitoids, even at the landscape scale. What causes host quality to vary and the subsequent mechanisms by which parasitoid population dynamics are affected can be complex. Here, we examine the indirect interaction of a plant pathogen with a parasitoid wasp. Under laboratory conditions, parasitoids from hosts fed fungus-infected plants weighed less than those from hosts fed uninfected plants, indicating that the fungus causes the hosts to be of poor quality. However, parasitoids reared from hosts fed fungal-infected diet also tended to be female, a characteristic associated with high host quality. The pathogen, herbivore and parasitoid persist regionally as metapopulations in a shared landscape in Aland, Finland. In an analysis of the metapopulation dynamics of the parasitoid over 6 years, the probability of colonization of a host population increased by more than twofold in patches occupied by the plant pathogen. While we cannot determine that the relationship is causal, a compelling explanation is that the plant pathogen facilitates the establishment by the parasitoid by increasing the fraction of female offspring. This is a novel mechanism of spatial multi-trophic level interactions.     

Temporal autocorrelation can enhance the persistence and abundance of metapopulations comprised of coupled sinks

In spatially heterogeneous landscapes, some habitats may be persistent sources, providing immigrants to sustain populations in unfavorable sink habitats (where extinction is inevitable without immigration). Recent theoretical and empirical studies of source-sink systems demonstrate that temporally variable local growth rates in sinks can substantially increase average abundance of a persisting population, provided that the variation is positively autocorrelated--in effect, temporal variation inflates average abundance. Here we extend these results to a metapopulation in which all habitat patches are sinks. Using numerical studies of a population with discrete generations (buttressed by analytic results), we show that temporal variation and moderate dispersal can jointly permit indefinite persistence of the metapopulation and that positive autocorrelation both lowers the magnitude of variation required for persistence and increases the average abundance of persisting metapopulations. These effects are weakened--but not destroyed--if variation in local growth rates is spatially synchronized and dispersal is localized. We show that the inflationary effect is robust to a number of extensions of the basic model, including demographic stochasticity and density dependence. Because ecological and environmental processes contributing to temporally variable growth rates in natural populations are typically autocorrelated, these observations may have important implications for species persistence.     

Metapopulation persistence in fragmented landscapes: significant interactions between genetic and demographic processes

We formulated a mathematical model in order to study the joint influence of demographic and genetic processes on metapopulation viability. Moreover, we explored the influence of habitat structure, matrix quality and disturbance on the interplay of these processes. We showed that the conditions that allow metapopulation persistence under the synergistic action of genetic and demographic processes depart significantly from predictions based on a mere superposition of the effects of each process separately. Moreover, an optimal dispersal rate exists that maximizes the range of survival rates of dispersers under which metapopulation persists and at the same time allows the largest sustainable patch removal and patch‐size reduction. The relative impact of patch removal and patch‐size reduction depends both on matrix quality and the dispersal strategy of the species: metapopulation persistence is more affected by patch‐size reduction (patch removal) for low (high)‐dispersing species, in presence of a low (high) quality matrix. Avoidance of inbreeding, through increased dispersal when the rate of inbreeding in a population is large, has positive effects on low‐dispersing species, but impairs the persistence of high‐dispersing species. Finally, size heterogeneity between patches largely influences metapopulation dynamics; the presence of large patches, even at the expense of other patches being smaller, can have positive effects on persistence in particular for species of low dispersing ability.     

Coevolutionary fine-tuning: evidence for genetic tracking between a specialist wasp parasitoid and its aphid host in a dual metapopulation interaction

In the interaction between two ecologically-associated species, the population structure of one species may affect the population structure of the other. Here, we examine the population structures of the aphid Metopeurum fuscoviride, a specialist on tansy Tanacetum vulgare, and its specialist primary hymenopterous parasitoid Lysiphlebus hirticornis, both of which are characterized by multivoltine life histories and a classic metapopulation structure. Samples of the aphid host and the parasitoid were collected from eight sites in and around Jena, Germany, where both insect species co-occur, and then were genotyped using suites of polymorphic microsatellite markers. The host aphid was greatly differentiated in terms of its spatial population genetic patterning, while the parasitoid was, in comparison, only moderately differentiated. There was a positive Mantel test correlation between pairwise shared allele distance (DAS) of the host and parasitoid, i.e. if host subpopulation samples were more similar between two particular sites, so were the parasitoid subpopulation samples. We argue that while the differences in the levels of genetic differentiation are due to the differences in the biology of the species, the correlations between host and parasitoid are indicative of dependence of the parasitoid population structure on that of its aphid host. The parasitoid is genetically tracking behind the aphid host, as can be expected in a classic metapopulation structure where host persistence depends on a delay between host and parasitoid colonization of the patch. The results may also have relevance to the Red Queen hypothesis, whereupon in the 'arms race' between parasitoid and its host, the latter 'attempts' to evolve away from the former.