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.     

Temporal genetic structuring of a specialist parasitoid,Lysiphlebus hirticornis Mackauer (Hymenoptera: Braconidae) attacking a specialist aphid on tansy

In insect species characterized by inbreeding, limited dispersal, and a metapopulation structure, high genetic differentiation and reduced genetic diversity within local populations are expected. Using the model system Lysiphlebus hirticornis Mackauer, a specialist parasitoid of the tansy aphid, Metopeurum fuscoviride Stroyan (Hemiptera: Aphididae), we examined within‐site temporal population dynamics and genetics, including molecular variation at the tansy plant level. Aphid‐parasitoid dynamics were surveyed and parasitoids sampled from 72 tansy plants at 11 sites in and around Jena, Germany, over one growing season. Thereafter, parasitoid samples were genotyped at 11 polymorphic microsatellite loci. Colonization, extinction, and recolonization events occurred during the season. Allele numbers and identities were highly variable over time. When samples from all sites were pooled, allele number over all loci showed a decreasing trend with time. At the level of sites, temporal changes in genetic diversity were more variable. Analysis of molecular variance revealed that samples at the plant level explained the highest variance compared to at site level. We conclude that the genetic structuring of this insect is very fine grained (i.e. at the tansy plant level) and the temporal genetic diversity is explained by a combination of extinction and recolonization events, as well as inbreeding. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 737–749.     

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.     

Metapopulation structure of the specialized herbivore Macrosiphoniella tanacetaria (Homoptera,Aphididae)

We investigated population dynamics, genetic diversity and spatial structure in the aphid species Macrosiphoniella tanacetaria, a specialist herbivore feeding on tansy, Tanacetum vulgare. Tansy plants (genets) consist of many shoots (ramets), and genets are grouped in sites. Thus, aphids feeding on tansy can cluster at the level of ramets, genets and sites. We studied aphid population dynamics in 1997 and 2001 and found that within sites: (i). at any time, aphids used only a fraction of the available ramets and genets; (ii). at the level of ramets, most aphid colonies survived only one week; (iii). at the level of genets, mean survival time was less than 4 weeks; and (iv). colonization and extinction events occurred throughout the season. We sampled aphids in seven sites in the Alsace region, France (4-45 km apart) and two sites in Germany in 1999 to study genetic structure within and between populations. Genetic analyses using nine microsatellite loci showed that: (i). genotypic variability was high, (ii). none of the populations was in Hardy-Weinberg equilibrium, (iii). heterozygote deficits and linkage disequilibria were frequent, and (iv). all populations were genetically differentiated, even at a small geographical scale. Renewed sampling of the Alsace sites in 2001 showed that three populations had become extinct and significant genetic changes had occurred in the remaining four populations. The frequencies of extinction and colonization events at several spatial scales suggest a hierarchical metapopulation structure for M. tanacetaria. Frequent population turnover and drift are likely causes for the genetic differentiation of M. tanacetaria populations.     

One step ahead: a parasitoid disperses farther and forms a wider geographic population than its fig wasp host

The structure of populations across landscapes influences the dynamics of their interactions with other species. Understanding the geographic structure of populations can thus shed light on the potential for interacting species to co‐evolve. Host–parasitoid interactions are widespread in nature and also represent a significant force in the evolution of plant–insect interactions. However, there have been few comparisons of population structure between an insect host and its parasitoid. We used microsatellite markers to analyse the population genetic structure of Pleistodontes imperialis sp. 1, a fig‐pollinating wasp of Port Jackson fig (Ficus rubiginosa), and its main parasitoid, Sycoscapter sp. A, in eastern Australia. Besides exploring this host–parasitoid system, our study also constitutes, to our knowledge, the first study of population structure in a nonpollinating fig wasp species. We collected matched samples of pollinators and parasitoids at several sites in two regions separated by up to 2000 km. We found that pollinators occupying the two regions represent distinct populations, but, in contrast, parasitoids formed a single population across the wide geographic range sampled. We observed genetic isolation by distance for each species, but found consistently lower F_ST and R_ST values between sites for parasitoids compared with pollinators. Previous studies have indicated that pollinators of monoecious figs can disperse over very long distances, and we provide the first genetic evidence that their parasitoids may disperse as far, if not farther. The contrasting geographic population structures of host and parasitoid highlight the potential for geographic mosaics in this important symbiotic system.     

Geographic structure with no evidence for host-associated lineages in European populations of Lysiphlebus testaceipes,an introduced biological control agent

Lysiphlebus testaceipes (Cress.) is an aphidiine parasitoid originally introduced to Europe as a biological control agent of citrus aphids in the Mediterranean. It has rapidly become widespread in coastal areas continuing gradually to expand inland. L. testaceipes exploited a large number of aphids in Europe, including new hosts and significantly changed the relative abundance of the native parasitoids. This behavior may reflect a broad oligophagy of the introduced parasitoid or it may require the evolution of host specialization that results in genetically differentiated subpopulations on different hosts. To address this issue we used the mitochondrial cytochrome oxidase subunit I and seven microsatellite loci to analyze the structure of genetic variation for L. testaceipes samples collected from 12 different aphid hosts across seven European countries, as well as some samples from Benin, Costa Rica, USA, Algeria and Libya for comparison. Only five COI haplotypes with moderate divergence were identified overall. There was no evidence for the association of haplotypes with different aphid hosts in the European samples, but there was geographic structuring in this variation. Haplotype diversity was highest in France, where L. testaceipes was introduced, but only a single haplotype was detected in areas of south-eastern Europe that were invaded subsequently. The analysis of microsatellite variation confirmed the lack of host-associated genetic structure, as well as differentiation between populations from south-western and south-eastern Europe. The parasitoid L. testaceipes in Europe is thus an opportunistic oligophagous species with a population structure shaped by the processes of introduction and expansion rather than by host exploitation.     

Strong dispersal in a parasitoid wasp overwhelms habitat fragmentation and host population dynamics

The population dynamics of a parasite depend on species traits, host dynamics and the environment. Those dynamics are reflected in the genetic structure of the population. Habitat fragmentation has a greater impact on parasites than on their hosts because resource distribution is increasingly fragmented for species at higher trophic levels. This could lead to either more or less genetic structure than the host, depending on the relative dispersal rates of species. We examined the spatial genetic structure of the parasitoid wasp Hyposoter horticola, and how it was influenced by dispersal, host population dynamics and habitat fragmentation. The host, the Glanville fritillary butterfly, lives as a metapopulation in a fragmented landscape in the Åland Islands, Finland. We collected wasps throughout the 50 by 70 km archipelago and determined the genetic diversity, spatial population structure and genetic differentiation using 14 neutral DNA microsatellite loci. We compared the genetic structure of the wasp with that of the host butterfly using published genetic data collected over the shared landscape. Using maternity assignment, we also identified full‐siblings among the sampled parasitoids to estimate the dispersal range of individual females. We found that because the parasitoid is dispersive, it has low genetic structure, is not very sensitive to habitat fragmentation and has less spatial genetic structure than its butterfly host. The wasp is sensitive to regional rather than local host dynamics, and there is a geographic mosaic landscape for antagonistic co‐evolution of host resistance and parasite virulence.     

EVOLUTION OF AN APHID-PARASITOID INTERACTION: VARIATION IN RESISTANCE TO PARASITISM AMONG APHID POPULATIONS SPECIALIZED ON DIFFERENT PLANTS

The evolution of associations between herbivorous insects and their parasitoids is likely to be influenced by the relationship between the herbivore and its host plants. If populations of specialized herbivorous insects are structured by their host plants such that populations on different hosts are genetically differentiated, then the traits affecting insect-parasitoid interactions may exhibit an associated structure. The pea aphid (Acyrthosiphon pisum) is a herbivorous insect species comprised of genetically distinct groups that are specialized on different host plants (Via 1991a, 1994). Here, we examine how the genetic differentiation of pea aphid populations on different host plants affects their interaction with a parasitoid wasp, Aphidius ervi. We performed four experiments. (1) By exposing pea aphids from both alfalfa and clover to parasitoids from both crops, we demonstrate that pea aphid populations that are specialized on alfalfa are successfully parasitized less often than are populations specialized on clover. This difference in parasitism rate does not depend upon whether the wasps were collected from alfalfa or clover fields. (2) When we controlled for potential differences in aphid and parasitoid behavior between the two host plants and ensured that aphids were attacked, we found that pea aphids from alfalfa were still parasitized less often than pea aphids from clover. Thus, the difference in parasitism rates is not due to behavior of either aphids or wasps, but appears to be a physiologically based difference in resistance to parasitism. (3) Replicates of pea aphid clones reared on their own host plant and on a common host plant, fava bean, exhibited the same pattern of resistance as above. Thus, there do not appear to be nutritional or secondary chemical effects on the level of physiological resistance in the aphids due to feeding on clover or alfalfa, and therefore the difference in resistance on the two crops appears to be genetically based. (4) We assayed for genetic variation in resistance among individual pea aphid clones collected from clover fields and found no detectable genetic variation for resistance to parasitism within two populations sampled from clover. This is in contrast to Henter and Via's (1995) report of abundant genetic variation in resistance to this parasitoid within a pea aphid population on alfalfa. Low levels of genetic variation may be one factor that constrains the evolution of resistance to parasitism in the populations of pea aphids from clover, leading them to remain more susceptible than populations of the same species from alfalfa.     

Asymmetry in host and parasitoid diffuse coevolution: when the red queen has to keep a finger in more than one pie

BACKGROUND: Coevolution between pairs of antagonistic species is generally considered an endless "arms race" between attack and defense traits to counteract the adaptive responses of the other species. PRESENTATION OF THE HYPOTHESIS: When more than two species are involved, diffuse coevolution of hosts and parasitoids could be asymmetric because consumers can choose their prey whereas preys do not choose their predator. This asymmetry may lead to differences in the rate of evolution of the antagonistic species in response to selection. The more long-standing the coevolution of a given pair of antagonistic populations, the higher should be the fitness advantage for the consumer. Therefore, the main prediction of the hypothesis is that the consumer trophic level is more likely to win the coevolution race. TESTING THE HYPOTHESIS: We propose testing the asymmetry hypothesis by focusing on the tritrophic system plant/aphid/aphid parasitoid. The analysis of the genetic variability in the virulence of several parasitoid populations and in the defenses of several aphid species or several clones of the same aphid species could be compared. Moreover, the analysis of the neutral population genetic structure of the parasitoid as a function of the aphid host, the plant host and geographic isolation may complement the detection of differences between host and parasitoid trophic specialization. IMPLICATIONS OF THE HYPOTHESIS: Genetic structures induced by the arms race between antagonistic species may be disturbed by asymmetry in coevolution, producing neither rare genotype advantages nor coevolutionary hotspots. Thus this hypothesis profoundly changes our understanding of coevolution and may have important implications in terms of pest management.     

Evolution of microparasites in spatially and genetically structured host populations: The example of RHDV infecting rabbits

Several studies have shown that classical results of microparasite evolution could not extend to the case where the host species shows an important spatial structure. Rabbit haemorrhagic disease virus (RHDV), responsible for rabbit haemorrhagic disease (RHD), which recently emerged in rabbits, has strains within a wide range of virulence, thus providing an interesting example of competition between strains infecting a host species with a metapopulation structure. In addition, rabbits may show a genetic diversity regarding RHDV susceptibility. In the present paper we use the example of the rabbit-RHDV interaction to study the competition between strains of a same microparasite in a host population that is both spatially and genetically structured. Using metapopulation models we show that the evolution of the microparasite is guided by a trade-off between its capacity to invade subpopulations potentially infected by other strains and its capacity to persist within the subpopulation. In such a context, host genetic diversity acts by reducing the number of hosts susceptible to each strain, often favouring more persistent—and generally less virulent—strains. We also show that even in a stochastic context where host genes regularly go locally extinct, the microparasite pressure helps maintain the genetic diversity in the long term while reinforcing gene loss risk in the short term. Finally, we study how different demographic and epidemiologic parameters affect the coevolution between the rabbit and RHDV.     

Honeydew feeding increased the longevity of two egg parasitoids of the pine processionary moth

The longevity of a generalist (Ooencyrtus pityocampae) and a specialist (Baryscapus servadeii) egg parasitoid of the pine processionary moth (Thau‐metopoea pityocampa) was compared under laboratory feeding conditions including water and honeydew from aphid species growing on maritime pine (Pinus pinaster) or pedunculate oak (Quercus robur). The longevity of both parasitoid species increased when specimens were fed with honeydew. This increase was larger for the generalist (3.7, 32.0 and 38.0 days) than for the specialist (3.0, 23.3 and 21.5 days) parasitoid species when fed with water, oak and pine aphid honeydew respectively. The phenology of the specialist species B. servadeii is well‐adapted with its host availability with or without food supply. The generalist species O. pityocampae could overlap its host emergence curve during 14.0–20.0 days when fed with oak and pine aphid honeydew respectively, vs. no overlap when no supplementary food was provided. Analysis of honeydew composition indicated that sugars and amino acids may have distinct effects on parasitoid longevity.     

Does variation in host plant association and symbiont infection of pea aphid populations induce genetic and behaviour differentiation of its main parasitoid, Aphidius ervi?

The host-associated differentiation (HAD) hypothesis states that higher trophic levels in parasitic associations should exhibit similar divergence in case of host sympatric speciation. We tested HAD on populations of Aphidius ervi the main parasitoid of the pea aphid Acyrthosiphon pisum, emerging from host populations specialized on either alfalfa or red clover. Host and parasitoid populations were assessed for genetic variation and structure, while considering geography, host plant and host aphid protective symbionts Regiella insecticola and Hamiltonella defensa as potential covariables. Cluster and hierarchical analyses were used to assess the contribution of these variables to population structure, based on genotyping pea aphids and associated A. ervi with microsatellites, and host aphid facultative symbionts with 16S rDNA markers. Pea aphid genotypes were clearly distributed in two groups closely corresponding with their plant origins, confirming strong plant associated differentiation of this aphid in North America. Overall parasitism by A. ervi averaged 21.5 % across samples, and many parasitized aphids producing a wasp hosted defensive bacteria, indicating partial or ineffective protective efficacy of these symbionts in the field. The A. ervi population genetic data failed to support differentiation according to the host plant association of their pea aphid host. Potential for parasitoid specialization was also explored in experiments where wasps from alfalfa and clover aphids were reciprocally transplanted on alternate hosts, the hypothesis being that wasp behaviour and parasitic stages should be most adapted to their host of origin. Results revealed higher probability of oviposition on the alfalfa aphids, but higher adult emergence success on red clover aphids, with no interaction as expected under HAD. We conclude that our study provides no support for the HAD in this system. We discuss factors that might impair A. ervi specialization on its divergent aphid hosts on alfalfa and clover.     

Genotype matching in a parasitoid–host genotypic food web: an approach for measuring effects of environmental change

Food webs typically quantify interactions between species, whereas evolution operates through the success of alleles within populations of a single species. To bridge this gap, we quantify genotypic interaction networks among individuals of a single specialized parasitoid species and its obligate to cyclically parthenogenetic aphid host along a climatic gradient. As a case study for the kinds of questions genotype food webs could be used to answer, we show that genetically similar parasitoids became more likely to attack genetically similar hosts in warmer sites (i.e. there was network‐wide congruence between the within‐species shared allelic distance of the parasitoid and that of its host). Narrowing of host‐genotype‐niche breadth by parasitoids could reduce resilience of the network to changes in host genetic structure or invasion by novel host genotypes and inhibit biological control. Thus, our approach can be easily used to detect changes to sub‐species‐level food webs, which may have important ecological and evolutionary implications, such as promoting host‐race specialization or the accelerated loss of functional diversity following extinctions of closely related genotypes.     

THE POTENTIAL FOR COEVOLUTION IN A HOST-PARASITOID SYSTEM. II. GENETIC VARIATION WITHIN A POPULATION OF WASPS IN THE ABILITY TO PARASITIZE AN APHID HOST

Much of the study of coevolution has focused on the adaptations that have resulted from interactions between species. For reciprocal evolution to occur, there must be genetic variation in each species for traits that directly affect their interaction. Here I report evidence of significant additive genetic variance within a population of parasitic wasps in the ability to successfully parasitize an aphid host. These data, combined with companion work documenting clonal variation in a population of aphids from the same site, provide evidence that within the same population both a host and its parasitoid have the potential for specific and reciprocal genetic interactions.     

Food web structure of aphids and their parasitoids in Belgian fruit agroecosystems

Community structures of aphids and their parasitoids were studied in fruit crop habitats of eastern Belgium in 2014 and 2015. Quantitative food webs of these insects were constructed separately for each year, and divided into subwebs on three host‐plant categories, fruit crop plants, non‐crop woody and shrub plants and non‐crop herbaceous plants. The webs were analyzed using the standard food web statistics designed for binary data. During the whole study period, 78 plant species were recorded as host plants of 71 aphid species, from which 48 parasitoid species emerged. The community structure, aphid / parasitoid species‐richness ratio and trophic link number varied between the two years, whereas the realized connectance between parasitoids and aphids was relatively constant. A new plant–aphid–parasitoid association for Europe was recorded. Dominant parasitoid species in the study sites were Ephedrus persicae, Binodoxys angelicae and Praon volucre: the first species was frequently observed on non‐crop trees and shrubs, but the other two on non‐crop herbaceous plants. The potential influence, through indirect interactions, of parasitoids on aphid communities was assessed with quantitative parasitoid‐overlap diagrams. Symmetrical links were uncommon, and abundant aphid species seemed to have large indirect effects on less abundant species. These results show that trophic indirect interactions through parasitoids may govern aphid populations in fruit crop habitats with various non‐crop plants, implying the importance for landscape management and biological control of aphid pests in fruit agroecosystems.     

Different genetic structures revealed resident populations of a specialist parasitoid wasp in contrast to its migratory host

Genetic comparisons of parasitoids and their hosts are expected to reflect ecological and evolutionary processes that influence the interactions between species. The parasitoid wasp, Cotesia vestalis, and its host diamondback moth (DBM), Plutella xylostella, provide opportunities to test whether the specialist natural enemy migrates seasonally with its host or occurs as resident population. We genotyped 17 microsatellite loci and two mitochondrial genes for 158 female adults of C. vestalis collected from 12 geographical populations, as well as nine microsatellite loci for 127 DBM larvae from six separate sites. The samplings covered both the likely source (southern) and immigrant (northern) areas of DBM from China. Populations of C. vestalis fell into three groups, pointing to isolation in northwestern and southwestern China and strong genetic differentiation of these populations from others in central and eastern China. In contrast, DBM showed much weaker genetic differentiation and high rates of gene flow. TESS analysis identified the immigrant populations of DBM as showing admixture in northern China. Genetic disconnect between C. vestalis and its host suggests that the parasitoid did not migrate yearly with its host but likely consisted of resident populations in places where its host could not survive in winter.     

Genetic structure and gene flow in a metapopulation of an endangered plant species,Silene tatarica

We investigated the distribution of genetic variation within and between seven subpopulations in a riparian population of Silene tatarica in northern Finland by using amplified fragment length polymorphism (AFLP) markers. A Bayesian approach-based clustering program indicated that the marker data contained not only one panmictic population, but consisted of seven clusters, and that each original sample site seems to consist of a distinct subpopulation. A coalescent-based simulation approach shows recurrent gene flow between subpopulations. Relative high FST values indicated a clear subpopulation differentiation. However, amova analysis and UPGMA-dendrogram did not suggest any hierarchical regional structuring among the subpopulations. There was no correlation between geographical and genetic distances among the subpopulations, nor any correlation between the subpopulation census size and amount of genetic variation. Estimates of gene flow suggested a low level of gene flow between the subpopulations, and the assignment tests proposed a few long-distance bidirectional dispersal events between the subpopulations. No apparent difference was found in within-subpopulation genetic diversity among upper, middle and lower regions along the river. Relative high amounts of linkage disequilibrium at subpopulation level indicated recent population bottlenecks or admixture, and at metapopulation levels a high subpopulation turnover rate. The overall pattern of genetic variation within and between subpopulations also suggested a 'classical' metapopulation structure of the species suggested by the ecological surveys.     

Effects of habitat reduction on the persistence of Ichneumon eumerus (Hymenoptera: Ichneumonidae), the specialist parasitoid of Maculinea rebeli (Lepidoptera: Lycaenidae)

The lycaenid butterfly, Maculinea rebeli, and its specialist parasitoid, Ichneumon eumerus, live in small, closed populations. Given the threatened status of the butterfly, it is reasonable to assume that its specialist parasitoid is even more vulnerable to local extinction than the butterfly host. Based on a mechanistic model recently developed for the tightly-woven community surrounding M. rebeli at a site in the Spanish Pyrenees, we investigate how the removal of habitat, and more particularly, specific habitat promoting the persistence of the butterfly, affects the population persistence of the parasitoid. Because of the relatively small impact of the parasitoid on the butterfly population in the Spanish Pyrenees, guidelines for conserving the parasitoid are only slightly more restrictive than those for its host. It is argued that at sites of more marginal quality for the butterfly than the reference site, achieving the dual aims of conserving both species will be more problematic. © Rapid Science Ltd. 1998     

Ecological specialization of the aphid Aphis gossypii Glover on cultivated host plants

Many plant-feeding insect species considered to be polyphagous are in fact composed of genetically differentiated sympatric populations that use different hosts and between which gene flow still exists. We studied the population genetic structure of the cotton-melon aphid Aphis gossypii that is considered as one of the most polyphagous aphid species. We used eight microsatellites to analyse the genetic diversity of numerous samples of A. gossypii collected over several years at a large geographical scale on annual crops from different plant families. The number of multilocus genotypes detected was extremely low and the genotypes were found to be associated with host plants. Five host races were unambiguously identified (Cucurbitaceae, cotton, eggplant, potato and chili- or sweet pepper). These host races were dominated by asexual clones. Plant transfer experiments using several specialized clones further confirmed the existence of host-associated trade-offs. Finally, both genetic and experimental data suggested that plants of the genus Hibiscus may be used as refuge for the specialized clones. Resource abundance is discussed as a key factor involved in the process of ecological specialization in A. gossypii.     

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.