Periodicity, persistence, and collapse in host-parasitoid systems with egg limitation

There is an emerging consensus that parasitoids are limited by the number of eggs which they can lay as well as the amount of time they can search for their hosts. Since egg limitation tends to destabilize host-parasitoid dynamics, successful control of insect pests by parasitoids requires additional stabilizing mechanisms such as heterogeneity in the distribution of parasitoid attacks and host density-dependence. To better understand how egg limitation, search limitation, heterogeneity in parasitoid attacks, and host density-dependence influence host-parasitoid dynamics, discrete time models accounting for these factors are analyzed. When parasitoids are purely egg-limited, a complete anaylsis of the host-parasitoid dynamics are possible. The analysis implies that the parasitoid can invade the host system only if the parasitoid's intrinsic fitness exceeds the host's intrinsic fitness. When the parasitoid can invade, there is a critical threshold, CV*>1, of the coefficient of variation (CV) of the distribution of parasitoid attacks that determines that outcome of the invasion. If parasitoid attacks sufficiently aggregated (i.e., CV>CV*), then the host and parasitoid coexist. Typically (in a topological sense), this coexistence is shown to occur about a periodic attractor or a stable equilibrium. If the parasitoid attacks are sufficiently random (i.e. CV1. When CV<1, the parasitoid exhibits highly oscillatory dynamics. Alternatively, when parasitoid attacks are sufficiently aggregated but not overly aggregated (i.e. CV>1 but close to 1), the host and parasitoid coexist about a stable equilibrium with low host densities. The implications of these results for classical biological control are discussed.     

Host-Parasitoid Dynamics and the Success of Biological Control When Parasitoids Are Prone to Allee Effects

In sexual organisms, low population density can result in mating failures and subsequently yields a low population growth rate and high chance of extinction. For species that are in tight interaction, as in host-parasitoid systems, population dynamics are primarily constrained by demographic interdependences, so that mating failures may have much more intricate consequences. Our main objective is to study the demographic consequences of parasitoid mating failures at low density and its consequences on the success of biological control. For this, we developed a deterministic host-parasitoid model with a mate-finding Allee effect, allowing to tackle interactions between the Allee effect and key determinants of host-parasitoid demography such as the distribution of parasitoid attacks and host competition. Our study shows that parasitoid mating failures at low density result in an extinction threshold and increase the domain of parasitoid deterministic extinction. When proned to mate finding difficulties, parasitoids with cyclic dynamics or low searching efficiency go extinct; parasitoids with high searching efficiency may either persist or go extinct, depending on host intraspecific competition. We show that parasitoids suitable as biocontrol agents for their ability to reduce host populations are particularly likely to suffer from mate-finding Allee effects. This study highlights novel perspectives for understanding of the dynamics observed in natural host-parasitoid systems and improving the success of parasitoid introductions.     

Differential cannibalism and population dynamics in a host-parasitoid system

The effects of host cannibalism on a host-parasitoid system were explored through experiment and modelling. In individual encounters between parasitized and unparasitized Plodia interpunctella larvae, parasitized larvae were more likely to be cannibalized. Inclusion of this differential cannibalism into a simple Lotka-Volterra-type model of host-parasitoid population dynamics generates alternative stable states-including stable coexistence and extinction of the parasitoid — which depend on starting conditions. Possible mechanisms for differential cannibalism, and its implications for studies of host-parasitoid populations and biological control programmes are discussed.     

Sex ratio genetics and the competitiveness of parasitic wasps

In the first part of the paper we analyse dynamics of the genetic mechanisms responsible for maintaining biased sex ratios in host-parasitoid interactions. We begin by reviewing recent results relating to the maintenance of sibmating in haplo-diploid populations. We then investigate the evolutionary stable sex ratio in populations in which all or some of the females mate with their brothers. In particular, we derive a diallelic one-locus model for studying evolutionary stable sex ratios in partially sibmating haplo-diploid populations. In the second part of the paper we review the impact of sex ratio on host-parasitoid populations. We then analyse how the sex ratio strategy of one parasitoid species may affect its interaction with another parasitoid species competing for the same host. In particular we show that, although a female biased sex ratio may enhance the inherent competitiveness of one species, it may also destabilize the ecological interaction of the three species so that all become extinct.     

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.     

Periodicity,persistence, and collapse in host–parasitoid systems with egg limitation

There is an emerging consensus that parasitoids are limited by the number of eggs which they can lay as well as the amount of time they can search for their hosts. Since egg limitation tends to destabilize host–parasitoid dynamics, successful control of insect pests by parasitoids requires additional stabilizing mechanisms such as heterogeneity in the distribution of parasitoid attacks and host density-dependence. To better understand how egg limitation, search limitation, heterogeneity in parasitoid attacks, and host density-dependence influence host–parasitoid dynamics, discrete time models accounting for these factors are analyzed. When parasitoids are purely egg-limited, a complete anaylsis of the host–parasitoid dynamics are possible. The analysis implies that the parasitoid can invade the host system only if the parasitoid’s intrinsic fitness exceeds the host’s intrinsic fitness. When the parasitoid can invade, there is a critical threshold, CV ^*>1, of the coefficient of variation (CV) of the distribution of parasitoid attacks that determines that outcome of the invasion. If parasitoid attacks sufficiently aggregated (i.e., CV>CV ^*), then the host and parasitoid coexist. Typically (in a topological sense), this coexistence is shown to occur about a periodic attractor or a stable equilibrium. If the parasitoid attacks are sufficiently random (i.e. CV<CV ^*), then the parasitoid drives the host to extinction. When parasitoids are weakly search-limited as well as egg-limited, coexistence about a global attractor occurs even if CV<CV ^*. However, numerical simulations suggest that the nature of this attractor depends critically on whether CV<1 or CV>1. When CV<1, the parasitoid exhibits highly oscillatory dynamics. Alternatively, when parasitoid attacks are sufficiently aggregated but not overly aggregated (i.e. CV>1 but close to 1), the host and parasitoid coexist about a stable equilibrium with low host densities. The implications of these results for classical biological control are discussed.     

Causes and consequences of small population size for a specialist parasitoid wasp

The parasitoid wasp Cotesia melitaearum lives in extremely small extinction-prone populations in the Åland islands of southwest Finland. Intensive observational data from two generations, a laboratory competition experiment, and 8 years of survey data were used to measure the causes, extent and consequences of small population size for this parasitoid. In the spring generations of 1999 and of 2000 we observed 21 out of 23 and 26 populations respectively, ranging in size from 2 to 103 parasitoid cocoons. Within these populations the fraction of individuals surviving to adulthood decreased with increasing parasitoid population size. The largest source of mortality was predation (44%) followed by parasitism (20%) and unknown causes (10%). In the field about 30% of the host butterfly larvae are parasitized by a competing parasitoid, Hyposoter horticola. A laboratory competition experiment showed that C. melitaearum eggs died when laid in post-diapause host larvae occupied by H. horticola. Consequently one-third of the progeny of the over-wintering generation of C. melitaearum from the field die as a result of larval competition. The survey of host and parasitoid population dynamics over 8 years showed that extinction of local host butterfly populations occupied by the parasitoid was not associated with current parasitoid population size. Over the same period small parasitoid populations were more likely to become extinct than large populations. However, parasitoid population size was not related to parasitoid extinction when the host also became extinct. These data suggest that the parasitoid populations are kept small through the action of natural enemies and competitors, some of which are density dependent. Local populations are so small that they become extinct frequently and rarely measurably affect the population dynamics of their host. It is likely that this parasitoid persists in Åland because of the spatial asynchrony of local population dynamics.     

Coexistence of multiple parasitoids on a single host due to differences in parasitoid phenology

There are many well-documented cases in which multiple parasitoids can coexist on a single host species. We examine a theoretical framework to assess whether parasitoid coexistence can be explained through differences in timing of parasitoid oviposition and parasitoid emergence. This study explicitly includes the phenology of host and parasitoid development and explores how this mechanism affects the population dynamics. Coexistence of the host with two parasitoids requires a balance between parasitoid fecundity and survival and occurs most readily if one parasitoid attacks earlier but emerges later than the other parasitoid. The host density can either be decreased or increased when a second coexisting parasitoid is introduced into the system. However, there always exists a single parasitoid type that is most effective at depressing the host density, although this type may not be successful due to parasitoid competition. The coexistence of multiple parasitoids also affects the population dynamics. For instance, population oscillations can be removed by the introduction of a second parasitoid. In general, subtle differences in parasitoid phenology can give rise to different outcomes in a host–multi-parasitoid system, and this may offer some insight into why establishing criteria for the ‘ideal’ biological control agent has been so challenging.     

Persistence of host and parasite populations subject to experimental size-selective removal

Predators have the potential to limit the spread of pathogens not only by selecting infected prey but also by shaping prey demographics. We tested this idea with an epidemiological experiment in which we simulated variable levels of size-selective predation on zooplankton hosts and monitored the persistence of host and parasite populations. In the absence of simulated predation, the virulent protozoan Caullerya mesnili frequently drove its host Daphnia galeata to extinction. Uninfected control populations showed lower extinction rates and higher average densities than infected populations in the absence of simulated predation (all of the latter went extinct or remained infected). With a weak removal rate of the largest hosts, the proportion of populations in which the parasite drove the host to extinction decreased, while the number of populations in which the host persisted and the parasite went extinct increased. Host-parasite coexistence was also observed in some cases. With intermediate levels of removal, most of the parasite populations went extinct, while the host populations persisted. With an even higher removal rate, Daphnia were driven to extinction as well. Thus, variation in one factor, size-selective mortality, resulted in four different patterns of population dynamics. Our results highlight the potential role of predation in shaping the epidemiology and community structure of host-parasite systems.     

The effect of evolution on host-parasitoid systems

It is well known that a simple first-order difference equation can exhibit complex population dynamics, such as sustained oscillations and chaos. An interesting problem is whether such oscillatory dynamics are expected to occur in real populations. This paper assumes that the resident system is composed of 1-host and 1-parasitoid and that only the host is allowed to evolve, but not the parasitoid. Based on the invasibility of a host to host-parasitoid systems, we investigate the dynamics of the host-parasitoid system favored by natural selection. We consider two cases. In the first case, the host's evolution involving both the intrinsic growth rate and the sensitivity to density is considered. In the second case, the host's evolution involving both the intrinsic growth rate and the vulnerability to the parasitoid is considered. In both cases, we see that the dynamics with a stable equilibrium will not be favored by natural selection without the trade-off between the host's traits which are allowed to evolve. The host-parasitoid system with a stable equilibrium will be eventually invaded by a host type that develops an unstable equilibrium with the parasitoid. If there is a trade-off between the host's traits which are allowed to evolve, a host-parasitoid system with a stable equilibrium can be favored by natural selection.     

Multiple attractors in host-parasitoid interactions: coexistence and extinction

This paper considers the dynamics of a two-dimensional discrete-time model for host-parasitoid interactions, and shows that the model has two attractors: the fixed point where two species coexist and a boundary cycle where the parasitoid is absent. The analysis with the Liapunov exponent confirms that this kind of bistability is common in this model. The generality of this phenomenon in host-parasitoid interactions is also discussed.     

Parasitoid foraging decisions mediated by artificial vibrations

In any host-parasitoid interaction, the detection of the host in its microhabitat is a crucial step in the interaction that has a major influence on the parasitoid's fitness. We used bioassays to investigate whether the parasitoid Sympiesis sericeicornis (Hymenoptera: Eulophidae) uses the vibrations of its host Phyllonorycter sp. gr. blancardella (Lepidoptera: Gracillariidae) to detect it in the leaf mine. From the large variety of signals produced by the host, we identified a unitary signal that could be reproduced artificially. We then sent this signal into emptied mines to analyse the reaction of the parasitoid. The wasps reacted by increasing both the time spent foraging on the mine and the number of insertions of the ovipositor to detect the ‘virtual’ host. This study is the first to show that parasitoids detect host vibrations.     

How populations persist when asexuality requires sex: the spatial dynamics of coping with sperm parasites

The twofold cost of sex implies that sexual and asexual reproduction do not coexist easily. Asexual forms tend to outcompete sexuals but may eventually suffer higher extinction rates, creating tension between short- and long-term advantages of different reproductive modes. The 'short-sightedness' of asexual reproduction takes a particularly intriguing form in gynogenetic species complexes, in which an asexual species requires sperm from a related sexual host species to trigger embryogenesis. Asexuals are then predicted to outcompete their host, after which neither species can persist. We examine whether spatial structure can explain continued coexistence of the species complex, and assess the evidence based on data on the Amazon molly (Poecilia formosa). A modification of the Levins metapopulation model creates two regions of good prospects for coexistence, connected by a region of poorer patch occupancy levels. In the first case, mate discrimination and/or niche differentiation keep local extinction rates low, and most patches contain both species; the other possibility resembles host-parasite dynamics where parasites frequently drive the host locally extinct. Several dynamical features are counterintuitive and relate to the parasitic nature of interactions in the species complex: for example, high local extinction rates of the asexual species can be beneficial for its own persistence. This creates a link from the evolution of sexual reproduction to that of prudent predation.     

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.     

Dynamical consequences of optimal host feeding on host-parasitoid population dynamics

This study examines the influence of various host-feeding patterns on host-parasitoid population dynamics. The following types of host-feeding patterns are considered: concurrent and non-destructive, non-concurrent and non-destructive, and non-concurrent and destructive. The host-parasitoid population dynamics is described by the Lotka-Volterra continuous-time model. This study shows that when parasitoids behave optimally, i.e. they maximize their fitness measured by the instantaneous per capita growth rate, the non-destructive type of host feeding stabilizes host-parasitoid dynamics. Other types of host feeding, i.e. destructive, concurrent, or non-concurrent, do not qualitatively change the neutral stability of the Lotka-Volterra model. Moreover, it is shown that the pattern of host feeding which maximizes parasitoid fitness is either non-concurrent and destructive, or concurrent and non-destructive host feeding, depending on the host abundance and parameters of the model. The effects of the adaptive choice of host-feeding patterns on host-parasitoid population dynamics are discussed.     

Prolonged diapause and the stability of host-parasitoid interactions

We investigated the effect on host-parasitoid dynamics of prolonged diapause, a feature of the life history of many animals living in unpredictable environments, by modifying the classical May (J. Anim. Ecol. 47 (1978) 833) host-parasitoid model. We considered three patterns of development of host and parasitoid: (a) prolonged parasitoid diapause controlled by host physiology, (b) parasitoid interference in host development, preventing parasitized hosts from prolonging diapause, and (c) host diapause independent of parasitoid attack. We found that single-year prolonged diapause shifted the boundaries of the May model towards a slight increase in stability. Longer periods of diapause prolongation had a stronger influence, but this influence remained modest if we considered realistic parameter values. In contrast to other recent studies, our results suggest that prolonged diapause does not necessarily compensate for the destabilizing effects of time lags on the influence of parasitoids on population dynamics.     

Niche modification and stability of competitive systems. II. Persistence of interspecific competitive systems with parasitoid wasps

Effects of niche shift in ecological time scale on the population dynamics of competing species were studied in the experimental populations of two parasitoid wasp species, Anisopteromalus calandrae and Heterospilus prosopidis (both are solitary parasites), on a host, the azuki bean weevil, Callosobruchus chinensis. Four resource conditions were set up with combination of kind of bean (azuki or black eye), and host distribution (uniform or clumped). In each resource condition, four developmental stages of hosts were provided as a resource spectrum for parasitoid wasps. Population dynamics of the two wasp populations were investigated in each resource condition in Multi-Generation Competitive Systems (MGCS), in which fresh hosts of four developmental stages were periodically introduced and were parasitized competitively by the two wasp species. Competitive coexistence of both wasps occurred in the azuki-clumped condition, where the peaks of the resource utilization curves separated in the two species; pupae in A. calandrae and the early or late fourth instar in H. prosopidis, A. calandrae was eliminated in the azuki-uniform condition and H. prosopidis went extinct in two black eye conditions irrespective of host distributions. The degrees of overlap of the resource utilization patterns of the two wasp species during population dynamics were not significantly different among resource conditions irrespective of the results of coexistence or extinction. Even in the azuki-clumped condition, however, extinction of A. calandrae was observed when resource partitioning could not be realized with only the late fourth instar larvae available to wasps. Further analytical experiments showed that parasitizing ability of A. calandrae increased with host density per bean with azuki beans, but A. calandrae could express higher parasitizing ability with black eye beans than H. prosopidis irrespective of host density per bean. The flexibility in parasitizing ability by A. calandrae for various host stages under different resource conditions was thought to be the major factor in determining the competitive coexistence or the extinction of either species under different resource conditions. The present experiments also suggested that different second-best host stages between competitors could be a major contributing factor to competitive coexistence.     

Spatial dynamics in a host–multiparasitoid community

1. The harlequin bug, a herbivore on bladderpod, is attacked by two specialist egg parasitoids Trissolcus murgantiae and Ooencyrtus johnsonii . Ooencyrtus can out-compete Trissolcus in the laboratory, but coexistence is the norm in field populations. Despite the heavy mortality inflicted by the two parasitoids, the host–parasitoid interaction is persistent in all sites that have been studied in southern California.
2. I manipulated inter-patch distances in a field experiment to determine whether spatial processes drive parasitoid coexistence and/or host–parasitoid dynamics. I first tested the hypothesis that the parasitoids coexist via a dispersal–competition trade-off. Both parasitoid species took significantly longer to colonize isolated patches than well-connected patches, suggesting that they have comparable dispersal abilities. Ooencyrtus did not exclude Trissolcus even when inter-patch distances were reduced to 25–30% of those observed in natural populations. These data suggest that parasitoid coexistence can occur in the absence of a dispersal advantage to the inferior competitor.
3. Since the treatments with isolated vs. well-connected patches did not differ in parasitoid composition, I next asked whether isolation would destabilize, or drive extinct, the host–multiparasitoid interaction. No local extinctions of bugs or parasitoids were observed during the 18-month study period. Bug populations in the isolated patches were no more variable than those in the well-connected patches. In fact, temporal variability in the experimentally isolated patches was comparable to that observed in highly isolated natural populations.
4. These data argue against a strong effect of spatial processes on host–parasitoid dynamics. Local processes may mediate both parasitoid coexistence as well as the host–parasitoid interaction.     

Consequences for host-parasitoid interactions of grazing-dependent habitat heterogeneity

1. Environmental heterogeneity can produce effects that cascade up to higher trophic levels and affect species interactions. We hypothesized that grazing-dependent habitat heterogeneity and grazing-independent host plant heterogeneity would influence directly and indirectly a host-parasitoid interaction in a woodland habitat. 2. Thistles were planted randomly in 20 birch woodlands, half of which are grazed by cattle. The abundances of two species of seed herbivore and their shared parasitoid were measured, and related to habitat and host-plant heterogeneity. 3. The presence of cattle grazing created a structurally and compositionally distinct plant assemblage from the ungrazed seminatural situation. Grazing did not affect the number or dispersion of the host plant underpinning the host-parasitoid interaction. 4. The density of one insect herbivore, Tephritis conura, and its parasitoid Pteromalus elevatus was significantly increased by the presence of cattle; but another herbivore, Xyphosia miliaria, was unaffected. The percentage of parasitism of T. conura was increased in grazed habitat occurring at twice the rate found in ungrazed habitat. 5. The increase in T. conura abundance was correlated with increased species richness and cover of forbs in grazed sites. This effect of grazing-dependent habitat variation on host insect density cascaded up to parasitoid density and percentage of parasitism. Habitat heterogeneity had a further direct, positive effect on parasitoid density and percentage of parasitism after controlling for host-insect density. 6. Independent of grazing, heterogeneity in host-plant flowering, architecture and stature further affected T. conura and its parasitoid's densities. Parasitoid density was also affected by the dispersion of the host plant. 7. A combination of habitat and host-plant scale environmental heterogeneity influenced a host-parasitoid interaction indirectly and directly, providing a rare example of an anthropogenic disturbance positively affecting a tertiary trophic level. This finding highlights the need to consider not only the importance of bottom-up effects for top-down processes, but also the role of environmental heterogeneity arising from anthropogenic disturbance for trophic interactions such as parasitism.     

Host-parasitoid association and diffuse coevolution: when to be a generalist?

In host-parasitoid communities, hosts are subjected to selective pressures from numerous parasitoid species, and parasitoids may attack several host species. The specificity of host resistance and parasitoid virulence is thus a key factor in host-parasitoid coevolution. A continuum of strategies exists, from strict specificity to a generalist strategy. The optimal level of specificity may differ in host and parasitoid. I investigated the optimal level of resistance specificity using a model in which the host could be attacked by two parasitoid species, with variable levels of defense specificity. The fitness of a parasitoid attacking two host species with different levels of virulence specificity was also modeled. Finally, a fluctuating environment was simulated by introducing variable probabilities of encounters between antagonistic species over several generations. If the frequency of encounters with the antagonistic species is fixed, then both host and parasitoid gain from a strategy of exclusive specialization toward the most frequent antagonist. If the frequency of encounters fluctuates between generations, generalist host resistance and partially specialist parasitoid virulence are favored. Generalist host resistance may be considered to be a bet-hedging response to an unpredictable environment. This asymmetry in host-parasitoid coevolution may account for some of the genetic structures observed in the field for host-parasitoid associations.