What can genetics tell us about population connectivity?

Genetic data are often used to assess ‘population connectivity’ because it is difficult to measure dispersal directly at large spatial scales. Genetic connectivity, however, depends primarily on the absolute number of dispersers among populations, whereas demographic connectivity depends on the relative contributions to population growth rates of dispersal vs. local recruitment (i.e. survival and reproduction of residents). Although many questions are best answered with data on genetic connectivity, genetic data alone provide little information on demographic connectivity. The importance of demographic connectivity is clear when the elimination of immigration results in a shift from stable or positive population growth to negative population growth. Otherwise, the amount of dispersal required for demographic connectivity depends on the context (e.g. conservation or harvest management), and even high dispersal rates may not indicate demographic interdependence. Therefore, it is risky to infer the importance of demographic connectivity without information on local demographic rates and how those rates vary over time. Genetic methods can provide insight on demographic connectivity when combined with these local demographic rates, data on movement behaviour, or estimates of reproductive success of immigrants and residents. We also consider the strengths and limitations of genetic measures of connectivity and discuss three concepts of genetic connectivity that depend upon the evolutionary criteria of interest: inbreeding connectivity, drift connectivity, and adaptive connectivity. To conclude, we describe alternative approaches for assessing population connectivity, highlighting the value of combining genetic data with capture‐mark‐recapture methods or other direct measures of movement to elucidate the complex role of dispersal in natural populations.     

Stability and distribution of predator–prey systems: local and regional mechanisms and patterns

Explaining the coexistence and distribution of species in time and space remains a fundamental challenge. While species coexistence depends on both local and regional mechanisms, it is sometimes unclear which role each mechanism takes in a given ecosystem. Consequently, it is very hard to predict the response of the ecosystem to environmental changes. Here, we develop a model to study spatial patterns of coexistence, focusing on predator–prey and host–parasite populations. We show, both theoretically and empirically, that these systems may exhibit both local and regional patterns and mechanisms of coexistence. Changes in environmental parameters, such as spatial connectivity, may lead to a transition from regional to local coexistence or it may lead directly to extinction, depending on demographic parameters. This demonstrates the importance of simultaneously analysing interacting mechanisms that act at different spatial scales to understand the response of ecosystems to environmental changes.     

Unifying population and landscape ecology with spatial capture–recapture

Spatial heterogeneity in the environment induces variation in population demographic rates and dispersal patterns, which result in spatio‐temporal variation in density and gene flow. Unfortunately, applying theory to learn about the role of spatial structure on populations has been hindered by the lack of mechanistic spatial models and inability to make precise observations of population state and structure. Spatial capture–recapture (SCR) represents an individual‐based analytic framework for overcoming this fundamental obstacle that has limited the utility of ecological theory. SCR methods make explicit use of spatial encounter information on individuals in order to model density and other spatial aspects of animal population structure, and they have been widely adopted in the last decade. We review the historical context and emerging developments in SCR models that enable the integration of explicit ecological hypotheses about landscape connectivity, movement, resource selection, and spatial variation in density, directly with individual encounter history data obtained by new technologies (e.g. camera trapping, non‐invasive DNA sampling). We describe ways in which SCR methods stand to advance the study of animal population ecology.     

Differentiating Migration and Dispersal Processes for Pond-Breeding Amphibians

Abstract Understanding the movement of animals is critical to many aspects of conservation such as spread of emerging disease, proliferation of invasive species, changes in land-use patterns, and responses to global climate change. Movement processes are especially important for amphibian management and conservation as species declines and extinctions worldwide become ever more apparent. To better integrate behavioral and ecological data on amphibian movements with our use of spatially explicit demographic models and guide effective conservation solutions, I present 1) a synopsis of the literature regarding behavior, ecology, and evolution of movement in pond-breeding amphibians possessing biphasic life cycles to distinguish between migration and dispersal processes, 2) a working hypothesis of juvenile-based dispersal, and 3) a discussion of conservation issues that follow from distinguishing the spatial and temporal movements of amphibians at different scales. I define amphibian migration as intrapopulational, round-trip movements toward and away from aquatic breeding sites. Population-level management, in general, can be focused on spatial scales of <1.0 km with attention focused on adult population and juveniles that remain near the natal wetland. I define amphibian dispersal as interpopulational, unidirectional movements from natal sites to other breeding sites. Metapopulation- or landscape-level management can be focused on movements among populations at spatial scales >1.0–10.0 km and on importance of terrestrial connectivity. The ultimate goal of conservation for amphibians should be long-term regional persistence by addressing management issues at both local and metapopulation scales.     

Spatial–temporal analysis of species range expansion: the case of the mountain pine beetle, Dendroctonus ponderosae

Aim The spatial extent of western Canada’s current epidemic of mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae, Scolytinae), is increasing. The roles of the various dispersal processes acting as drivers of range expansion are poorly understood for most species. The aim of this paper is to characterize the movement patterns of the mountain pine beetle in areas where range expansion is occurring, in order to describe the fine‐scale spatial dynamics of processes associated with mountain pine beetle range expansion. Location Three regions of Canada’s Rocky Mountains: Kicking Horse Pass, Yellowhead Pass and Pine Pass. Methods Data on locations of mountain pine beetle‐attacked trees of predominantly lodgepole pine (Pinus contorta var. latifolia) were obtained from annual fixed‐wing aircraft surveys of forest health and helicopter‐based GPS surveys of mountain pine beetle‐damaged areas in British Columbia and Alberta. The annual (1999–2005) spatial extents of outbreak ranges were delineated from these data. Spatial analysis was conducted using the spatial–temporal analysis of moving polygons (STAMP), a recently developed pattern‐based approach. Results We found that distant dispersal patterns (spot infestations) were most often associated with marginal increases in the areal size of mountain pine beetle range polygons. When the mountain pine beetle range size increased rapidly relative to the years examined, local dispersal patterns (adjacent infestation) were more common. In Pine Pass, long‐range dispersal (> 2 km) markedly extended the north‐east border of the mountain pine beetle range. In Yellowhead Pass and Kicking Horse Pass, the extension of the range occurred incrementally via ground‐based spread. Main conclusions Dispersal of mountain pine beetle varies with geography as well as with host and beetle population dynamics. Although colonization is mediated by habitat connectivity, during periods of low overall habitat expansion, dispersal to new distant locations is common, whereas during periods of rapid invasion, locally connected spread is the dominant mode of dispersal. The propensity for long‐range transport to establish new beetle populations, and thus to be considered a driver of range expansion, is likely to be determined by regional weather patterns, and influenced by local topography. We conclude that STAMP appears to be a useful approach for examining changes in biogeograpical ranges, with the potential to reveal both fine‐ and large‐scale patterns.     

Spatiotemporal Model of Barley and Cereal Yellow Dwarf Virus Transmission Dynamics with Seasonality and Plant Competition

Many generalist pathogens are influenced by the spatial distributions and relative abundances of susceptible host species. The spatial structure of host populations can influence patterns of infection incidence (or disease outbreaks), and the effects of a generalist pathogen on host community dynamics in a spatially heterogeneous community may differ from predictions derived via simple models. In this paper, we model the transmission of a generalist pathogen within a patch framework that incorporates the movement of vectors between discrete host patches to investigate the effects of local host community composition and vector movement rates on disease dynamics.     

School’s Out: Seasonal Variation in the Movement Patterns of School Children

School children are core groups in the transmission of many common infectious diseases, and are likely to play a key role in the spatial dispersal of disease across multiple scales. However, there is currently little detailed information about the spatial movements of this epidemiologically important age group. To address this knowledge gap, we collaborated with eight secondary schools to conduct a survey of movement patterns of school pupils in primary and secondary schools in the United Kingdom. We found evidence of a significant change in behaviour between term time and holidays, with term time weekdays characterised by predominately local movements, and holidays seeing much broader variation in travel patterns. Studies that use mathematical models to examine epidemic transmission and control often use adult commuting data as a proxy for population movements. We show that while these data share some features with the movement patterns reported by school children, there are some crucial differences between the movements of children and adult commuters during both term-time and holidays.     

Larval dispersal and movement patterns of coral reef fishes,and implications for marine reserve network design

Well‐designed and effectively managed networks of marine reserves can be effective tools for both fisheries management and biodiversity conservation. Connectivity, the demographic linking of local populations through the dispersal of individuals as larvae, juveniles or adults, is a key ecological factor to consider in marine reserve design, since it has important implications for the persistence of metapopulations and their recovery from disturbance. For marine reserves to protect biodiversity and enhance populations of species in fished areas, they must be able to sustain focal species (particularly fishery species) within their boundaries, and be spaced such that they can function as mutually replenishing networks whilst providing recruitment subsidies to fished areas. Thus the configuration (size, spacing and location) of individual reserves within a network should be informed by larval dispersal and movement patterns of the species for which protection is required. In the past, empirical data regarding larval dispersal and movement patterns of adults and juveniles of many tropical marine species have been unavailable or inaccessible to practitioners responsible for marine reserve design. Recent empirical studies using new technologies have also provided fresh insights into movement patterns of many species and redefined our understanding of connectivity among populations through larval dispersal. Our review of movement patterns of 34 families (210 species) of coral reef fishes demonstrates that movement patterns (home ranges, ontogenetic shifts and spawning migrations) vary among and within species, and are influenced by a range of factors (e.g. size, sex, behaviour, density, habitat characteristics, season, tide and time of day). Some species move <0.1–0.5 km (e.g. damselfishes, butterflyfishes and angelfishes), <0.5–3 km (e.g. most parrotfishes, goatfishes and surgeonfishes) or 3–10 km (e.g. large parrotfishes and wrasses), while others move tens to hundreds (e.g. some groupers, emperors, snappers and jacks) or thousands of kilometres (e.g. some sharks and tuna). Larval dispersal distances tend to be <5–15 km, and self‐recruitment is common. Synthesising this information allows us, for the first time, to provide species, specific advice on the size, spacing and location of marine reserves in tropical marine ecosystems to maximise benefits for conservation and fisheries management for a range of taxa. We recommend that: (i) marine reserves should be more than twice the size of the home range of focal species (in all directions), thus marine reserves of various sizes will be required depending on which species require protection, how far they move, and if other effective protection is in place outside reserves; (ii) reserve spacing should be <15 km, with smaller reserves spaced more closely; and (iii) marine reserves should include habitats that are critical to the life history of focal species (e.g. home ranges, nursery grounds, migration corridors and spawning aggregations), and be located to accommodate movement patterns among these. We also provide practical advice for practitioners on how to use this information to design, evaluate and monitor the effectiveness of marine reserve networks within broader ecological, socioeconomic and management contexts.     

From Shelf to Shelf: Assessing Historical and Contemporary Genetic Differentiation and Connectivity across the Gulf of Mexico in Gag,Mycteroperca microlepis

Describing patterns of connectivity among populations of species with widespread distributions is particularly important in understanding the ecology and evolution of marine species. In this study, we examined patterns of population differentiation, migration, and historical population dynamics using microsatellite and mitochondrial loci to test whether populations of the epinephelid fish, Gag, Mycteroperca microlepis, an important fishery species, are genetically connected across the Gulf of Mexico and if so, whether that connectivity is attributable to either contemporary or historical processes. Populations of Gag on the Campeche Bank and the West Florida Shelf show significant, but low magnitude, differentiation. Time since divergence/expansion estimates associated with historical population dynamics indicate that any population or spatial expansions indicated by population genetics would have likely occurred in the late Pleistocene. Using coalescent-based approaches, we find that the best model for explaining observed spatial patterns of contemporary genetic variation is one of asymmetric gene flow, with movement from Campeche Bank to the West Florida Shelf. Both estimated migration rates and ecological data support the hypothesis that Gag populations throughout the Gulf of Mexico are connected via present day larval dispersal. Demonstrating this greatly expanded scale of connectivity for Gag highlights the influence of “ghost” populations (sensu Beerli) on genetic patterns and presents a critical consideration for both fisheries management and conservation of this and other species with similar genetic patterns.     

The Role of Human Transportation Networks in Mediating the Genetic Structure of Seasonal Influenza in the United States

Recent studies have demonstrated the importance of accounting for human mobility networks when modeling epidemics in order to accurately predict spatial dynamics. However, little is known about the impact these movement networks have on the genetic structure of pathogen populations and whether these effects are scale-dependent. We investigated how human movement along the aviation and commuter networks contributed to intra-seasonal genetic structure of influenza A epidemics in the continental United States using spatially-referenced hemagglutinin nucleotide sequences collected from 2003–2013 for both the H3N2 and H1N1 subtypes. Comparative analysis of these transportation networks revealed that the commuter network is highly spatially-organized and more heavily traveled than the aviation network, which instead is characterized by high connectivity between all state pairs. We found that genetic distance between sequences often correlated with distance based on interstate commuter network connectivity for the H1N1 subtype, and that this correlation was not as prevalent when geographic distance or aviation network connectivity distance was assessed against genetic distance. However, these patterns were not as apparent for the H3N2 subtype at the scale of the continental United States. Finally, although sequences were spatially referenced at the level of the US state of collection, a community analysis based on county to county commuter connections revealed that commuting communities did not consistently align with state geographic boundaries, emphasizing the need for the greater availability of more specific sequence location data. Our results highlight the importance of utilizing host movement data in characterizing the underlying genetic structure of pathogen populations and demonstrate a need for a greater understanding of the differential effects of host movement networks on pathogen transmission at various spatial scales.     

Contrasting demographic and genetic estimates of dispersal in the endangered Coahuilan box turtle: a contemporary approach to conservation

The evolutionary viability of an endangered species depends upon gene flow among subpopulations and the degree of habitat patch connectivity. Contrasting population connectivity over ecological and evolutionary timescales may provide novel insight into what maintains genetic diversity within threatened species. We employed this integrative approach to evaluating dispersal in the critically endangered Coahuilan box turtle (Terrapene coahuila) that inhabits isolated wetlands in the desert‐spring ecosystem of Cuatro Ciénegas, Mexico. Recent wetland habitat loss has altered the spatial distribution and connectivity of habitat patches; and we therefore predicted that T. coahuila would exhibit limited movement relative to estimates of historic gene flow. To evaluate contemporary dispersal patterns, we employed mark–recapture techniques at both local (wetland complex) and regional (intercomplex) spatial scales. Gene flow estimates were obtained by surveying genetic variation at nine microsatellite loci in seven subpopulations located across the species’ geographical range. The mark–recapture results at the local spatial scale reveal frequent movement among wetlands that was unaffected by interwetland distance. At the regional spatial scale, dispersal events were relatively less frequent between wetland complexes. The complementary analysis of population genetic substructure indicates strong historic gene flow (global F_ST = 0.01). However, a relationship of genetic isolation by distance across the geographical range suggests that dispersal limitation exists at the regional scale. Our approach of contrasting direct and indirect estimates of dispersal at multiple spatial scales in T. coahuila conveys a sustainable evolutionary trajectory of the species pending preservation of threatened wetland habitats and a range‐wide network of corridors.     

Parasites as indicators of movement and population connectivity of a non-diadromous, tropical estuarine teleost: king threadfin Polydactylus macrochir

Temporal and spatial patterns in parasite assemblages were examined to evaluate the degree of movement and connectivity of post-recruitment life-history stages of a large, non-diadromous tropical estuarine teleost, king threadfin Polydactylus macrochir, collected from 18 locations across northern Australia. Ten parasites types (juvenile stages of two nematodes and seven cestodes, and adults of an acanthocephalan) were deemed to be suitable for use as biological tags, in that they were considered to have a long residence time in the fish, were relatively easy to find and were morphologically very different to each other which aided discrimination. Univariate and discriminant function analysis of these parasites revealed little difference in temporal replicates collected from five locations, suggesting that the parasite communities were stable over the timeframes explored. Univariate, discriminant function, and Bray-Curtis similarity analyses indicated significant spatial heterogeneity, with Bray-Curtis classification accuracies ranging from 55 to 100% for locations in north-western and northern Australia, 24 to 88% in the Gulf of Carpentaria, and 39 to 88% on the east coast of Queensland. Few differences were observed among locations separated by <200 km. The observed patterns of parasite infection are in agreement with concurrent studies of movement and connectivity of P. macrochir in that they indicate a complex population structure across northern Australia. These results should be considered when reviewing the management arrangements for this species.     

A gravity model for the spread of a pollinator-borne plant pathogen

Many pathogens of plants are transmitted by arthropod vectors whose movement between individual hosts is influenced by foraging behavior. Insect foraging has been shown to depend on both the quality of hosts and the distances between hosts. Given the spatial distribution of host plants and individual variation in quality, vector foraging patterns may therefore produce predictable variation in exposure to pathogens. We develop a "gravity" model to describe the spatial spread of a vector-borne plant pathogen from underlying models of insect foraging in response to host quality using the pollinator-borne smut fungus Microbotryum violaceum as a case study. We fit the model to spatially explicit time series of M. violaceum transmission in replicate experimental plots of the white campion Silene latifolia. The gravity model provides a better fit than a mean field model or a model with only distance-dependent transmission. The results highlight the importance of active vector foraging in generating spatial patterns of disease incidence and for pathogen-mediated selection for floral traits.     

Patterns of connectivity among populations of a coral reef fish

Knowledge of the patterns and scale of connectivity among populations is essential for the effective management of species, but our understanding is still poor for marine species. We used otolith microchemistry of newly settled bicolor damselfish (Stegastes partitus) in the Mesoamerican Reef System (MRS), Western Caribbean, to investigate patterns of connectivity among populations over 2 years. First, we assessed spatial and temporal variability in trace elemental concentrations from the otolith edge to make a ‘chemical map’ of potential source reef(s) in the region. Significant otolith chemical differences were detected at three spatial scales (within-atoll, between-atolls, and region-wide), such that individuals were classified to locations with moderate (52 % jackknife classification) to high (99 %) accuracy. Most sites at Turneffe Atoll, Belize showed significant temporal variability in otolith concentrations on the scale of 1–2 months. Using a maximum likelihood approach, we estimated the natal source of larvae recruiting to reefs across the MRS by comparing ‘natal’ chemical signatures from the otolith of recruits to the ‘chemical map’ of potential source reef(s). Our results indicated that populations at both Turneffe Atoll and Banco Chinchorro supply a substantial amount of individuals to their own reefs (i.e., self-recruitment) and thus emphasize that marine conservation and management in the MRS region would benefit from localized management efforts as well as international cooperation.     

From one host to another: tracking vector movements using microsatellite markers

In principle, the solution to stopping the spread of any vectorborne pathogen is a simple one – just stop infected vectors from biting new hosts and the pathogen cannot spread. Importantly, this does not necessarily require killing all vectors, or protecting all hosts. Transmission only occurs when an infected vector moves to a new host, and so knowing how vectors move between hosts in nature and how they choose hosts is crucial to understanding transmission. For example, the infection status of a potential vector or that of a potential host would have a huge influence on pathogen transmission if it affected vector movement or host choice. Remarkably little is known about how vectors move between and choose hosts in nature, in part because of the logistical difficulties of tracking vector movement. This is why the article by Levin and Parker ( 2014 ) in this issue of Molecular Ecology is so exciting.     

Connectivity, sustainability, and yield: bridging the gap between conventional fisheries management and marine protected areas

A substantial shift toward use of marine protected areas (MPAs) for conservation and fisheries management is currently underway. This shift to explicit spatial management presents new challenges and uncertainties for ecologists and resource managers. In particular, the potential for MPAs to change population sustainability, fishery yield, and ecosystem properties depends on the poorly understood consequences of three critical forms of connectivity over space: larval dispersal, juvenile and adult swimming, and movement of fishermen. Conventional fishery management describes the dynamics and current status of fish populations, with increasing recent emphasis on sustainability, often through reference points that reflect individual replacement. These compare lifetime egg production (LEP) to a critical replacement threshold (CRT) whose value is uncertain. Sustainability of spatially distributed populations also depends on individual replacement, but through all possible paths created by larval dispersal and LEP at each location. Model calculations of spatial replacement considering larval connectivity alone indicate sustainability and yield depend on species dispersal distance and the distribution of LEP created by species habitat distribution and fishing mortality. Adding MPAs creates areas with high LEP, increasing sustainability, but not necessarily yield. Generally, short distance dispersers will persist in almost all MPAs, while sustainability of long distance dispersers requires a specific density of MPAs along the coast. The value of that density also depends on the uncertain CRT, as well as fishing rate. MPAs can increase yield in areas with previously low LEP but for short distance dispersers, high yields will require many small MPAs. The paucity of information on larval dispersal distances, especially in cases with strong advection, renders these projections uncertain. Adding juvenile and adult movement to these calculations reduces LEP near the edges in MPAs, if movement is within a home-range, but more broadly over space if movement is diffusive. Adding movement of fishermen shifts effort on the basis of anticipated revenues and fishing costs, leading to lower LEP near ports, for example. Our evolving understanding of connectivity in spatial management could form the basis for a new, spatially oriented replacement reference point for sustainability, with associated new uncertainties.     

Vector-Borne Pathogen and Host Evolution in a Structured Immuno-Epidemiological System

Vector-borne disease transmission is a common dissemination mode used by many pathogens to spread in a host population. Similar to directly transmitted diseases, the within-host interaction of a vector-borne pathogen and a host’s immune system influences the pathogen’s transmission potential between hosts via vectors. Yet there are few theoretical studies on virulence–transmission trade-offs and evolution in vector-borne pathogen–host systems. Here, we consider an immuno-epidemiological model that links the within-host dynamics to between-host circulation of a vector-borne disease. On the immunological scale, the model mimics antibody-pathogen dynamics for arbovirus diseases, such as Rift Valley fever and West Nile virus. The within-host dynamics govern transmission and host mortality and recovery in an age-since-infection structured host-vector-borne pathogen epidemic model. By considering multiple pathogen strains and multiple competing host populations differing in their within-host replication rate and immune response parameters, respectively, we derive evolutionary optimization principles for both pathogen and host. Invasion analysis shows that the \({\mathcal {R}}_0\) maximization principle holds for the vector-borne pathogen. For the host, we prove that evolution favors minimizing case fatality ratio (CFR). These results are utilized to compute host and pathogen evolutionary trajectories and to determine how model parameters affect evolution outcomes. We find that increasing the vector inoculum size increases the pathogen \({\mathcal {R}}_0\), but can either increase or decrease the pathogen virulence (the host CFR), suggesting that vector inoculum size can contribute to virulence of vector-borne diseases in distinct ways.     

The spatial spread and final size of models for the deterministic host-vector epidemic

A disease is considered which is transferred between two populations, termed hosts and vectors. The disease is transmitted solely from infected vector to uninfected host and from infected host to uninfected vector. Two models are formulated in which infectious individuals are introduced at time t = 0 into the populations of susceptibles, thus triggering an epidemic through those populations. Conditions are established for a major epidemic to occur, and the final size of the epidemic is obtained for these models when no spatial aspect is considered. When a spatial aspect is included in the models, again the condition for a major epidemic is obtained. The pandemic theorem is proved rigorously, giving a lower bound for the proportion of each population, at each point, who eventually suffer the epidemic. The behavior a long way from the initial focus of infection is also rigorously obtained.     

An empirical test of a diffusion model: predicting clouded apollo movements in a novel environment

Functional connectivity is a fundamental concept in conservation biology because it sets the level of migration and gene flow among local populations. However, functional connectivity is difficult to measure, largely because it is hard to acquire and analyze movement data from heterogeneous landscapes. Here we apply a Bayesian state-space framework to parameterize a diffusion-based movement model using capture-recapture data on the endangered clouded apollo butterfly. We test whether the model is able to disentangle the inherent movement behavior of the species from landscape structure and sampling artifacts, which is a necessity if the model is to be used to examine how movements depend on landscape structure. We show that this is the case by demonstrating that the model, parameterized with data from a reference landscape, correctly predicts movements in a structurally different landscape. In particular, the model helps to explain why a movement corridor that was constructed as a management measure failed to increase movement among local populations. We illustrate how the parameterized model can be used to derive biologically relevant measures of functional connectivity, thus linking movement data with models of spatial population 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.