Evolution of spatially structured host–parasite interactions

Spatial structure has dramatic effects on the demography and the evolution of species. A large variety of theoretical models have attempted to understand how local dispersal may shape the coevolution of interacting species such as host–parasite interactions. The lack of a unifying framework is a serious impediment for anyone willing to understand current theory. Here, we review previous theoretical studies in the light of a single epidemiological model that allows us to explore the effects of both host and parasite migration rates on the evolution and coevolution of various life‐history traits. We discuss the impact of local dispersal on parasite virulence, various host defence strategies and local adaptation. Our analysis shows that evolutionary and coevolutionary outcomes crucially depend on the details of the host–parasite life cycle and on which life‐history trait is involved in the interaction. We also discuss experimental studies that support the effects of spatial structure on the evolution of host–parasite interactions. This review highlights major similarities between some theoretical results, but it also reveals an important gap between evolutionary and coevolutionary models. We discuss possible ways to bridge this gap within a more unified framework that would reconcile spatial epidemiology, evolution and coevolution.     

Parasite establishment and host extinction in model communities

Studies of host–parasite dynamics usually consider one, or at most two, host species, neglecting the possible effects of other species on the focal hosts and vice versa. To explore the interaction of community structure with host–parasite dynamics, we model the invasion of stable communities of varying size by a parasite. The communities are generated with random interaction coefficients and connectance 0.5. Each community is invaded by parasites with different values of virulence (disease-induced host mortality rate), specificity and transmission rate. The result of each invasion is determined by numerically simulating the dynamics of the community. We classify the outcomes by whether the parasite successfully establishes in the focal host population(s), and, if so, by the proportion of host and non-host species that go extinct as a result of the parasite's introduction. We discuss how the structure of the community and the interaction between hosts and other species affect several important processes of disease ecology: the density threshold for parasite invasion, extinction cascades caused by the parasite, and the frequency of extinctions of hosts and non-hosts. In our simulated communities, non-host species went extinct more frequently than hosts, suggesting the importance of the community context of disease. In some cases, the parasite's invasion induced regular population cycles in the previously stable community.     

Relative migration rates and local adaptation in a mosquito-protozoan interaction

Theory predicts that the direction of local adaptation depends on the relative migration rates of hosts and parasites. Here we measured relative migration rates and tested for local adaptation in the interaction between a tree hole mosquito (Ochlerotatus sierrensis) and a protozoan parasite (Lambornella clarki). We found strong support for the hypothesis that the host migrates more than its parasite. Hosts colonized artificial tree holes in the field at a much higher rate than the parasite. Field releases of the parasite demonstrated that it colonizes and persists in natural tree holes where it was previously absent, suggesting that parasite distribution is limited by its migratory ability. Although the host migrates more than its parasite, we found no evidence for local adaptation by hosts and some evidence for local adaptation by parasites. Other life history traits of the host and parasite may also influence patterns in local adaptation, particularly parasite virulence and host dormancy.     

Experimental evolution of local parasite maladaptation

Sign and magnitude of local adaptation in host–parasite systems may vary with ecological, epidemiological or genetic parameters. To investigate the role of host genetic background, we established long‐term experimental populations of different genotypes of the protozoan Paramecium caudatum, infected with the bacterial parasite Holospora undulata. We observed the evolution of an overall pattern of parasite local maladaptation for infectivity, indicating a general coevolutionary disadvantage of this parasite. Maladaptation extended to host populations with the same genetic background, similar to extending from the local to a higher regional level in natural populations. Patterns for virulence were qualitatively similar, but with less statistical support. A nonsignificant correlation with levels of (mal)adaptation for infectivity suggests independent evolution of these traits. Our results indicate similar (co)evolutionary trajectories in populations with different genetic backgrounds. Nonetheless, the correlated clines of genetic distance and parasite performance illustrate how genetic background can shape spatial gradients of local adaptation.     

Coevolution between hosts and parasites with partially overlapping geographic ranges

Many host species interact with a specific parasite within only a fraction of their geographical range. Where host and parasite overlap geographically, selection may be reciprocal constituting a coevolutionary hot spot. Host evolution, however, may be driven primarily by selection imposed by alternative biotic or abiotic factors that occur outside such hot spots. To evaluate the importance of coevolutionary hot spots for host and parasite evolution, we analyse a spatially explicit genetic model for a host that overlaps with a parasite in only part of its geographical range. Our results show that there is a critical amount of overlap beyond which reciprocal selection leads to a coevolutionary response in the host. This critical amount of overlap depends upon the explicit spatial configuration of hot spots. When the amount of overlap exceeds this first critical level, host-parasite coevolution commonly generates stable allele frequency clines rather than oscillations. It is within this region that one of the primary predictions of the geographic mosaic theory is realized, and local maladaptation is prevalent in both species. Past a further threshold of overlap between the species oscillations do evolve, but allele frequencies in both species are spatially synchronous and local maladaptation is absent in both species. A consequence of such transitions between coevolutionary dynamics is that parasite adaptation is inversely proportional to the fraction of its host's range that it occupies. Hence, as the geographical range of a parasite increases, it becomes increasingly maladapted to the host. This suggests a novel mechanism through which the geographical range of parasites may be limited.     

Temperature-mediated patterns of local adaptation in a natural plant-pathogen metapopulation

There have been numerous investigations of parasite local adaptation, a phenomenon important from the perspectives of both basic and applied evolutionary ecology. Recent work has demonstrated that temperature has striking effects on parasite performance by mediating trade-offs in parasite life history and through genotype × environment interactions. To test whether parasite local adaptation is mediated by temperature, I measured the performance of sympatric populations against allopatric populations of a fungal pathogen, Podosphaera plantaginis , on its host Plantago lanceolata , across a temperature gradient. I used data on parasite life history and epidemiology to derive fitness estimates to measure local adaptation. The results demonstrate unambiguously that trajectories of host–parasite co-evolution are tightly coupled with parasite adaptation to the abiotic habitat, as the strength, and even direction, of local adaptation varied with temperature. Patterns of local adaptation further depended on how parasite fitness was estimated, highlighting the importance of choosing relevant fitness measures in studies of local adaptation.     

Coevolution of daily activity timing in a host–parasite system

Coevolutionary theories applied in the study of host–parasite systems indicate that lineages exhibit progressive trends in response to reciprocal selective pressures. Avian brood parasites have generated intense interest as models for coevolutionary processes. Similar to avian cuckoos, Polistes wasp social parasites usurp a nest and exploit the parental care of a congeneric species to rear their own brood. In the present study, we show a coevolutionary arms race in the daily activity pattern in a Polistes host–parasite pair. We measured the daily activity rate, in constant laboratory conditions, of both host and parasite females during the period in which nest usurpations occur. The parasites showed a hyperkinesis in the middle of the day. As the field observations suggested, this mid-day activity is used to perform host nest usurpation attempts. Timing the usurpations allows the parasite to maximize its usurpation attempts during daytime when the host defence is lower. A field comparison of host presence on the nest in two populations with different parasitism rates showed that populations under strong parasitic pressure exhibit timing counteradaptations to optimize nest defence. This study provides the first example of a mutual coadaptation in timing activity in a parasite–host system.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 399–405.     

Heterogeneous selection promotes maintenance of polymorphism in host–parasite interactions

Polymorphism in loci affecting host resistance and parasite virulence is characteristic for nearly all species and this genetic variation is considered to have profound consequences for the patterns of disease incidence, prevalence and evolution. The gene-for-gene (GFG) system is a well-characterized genetic interaction of host recognition and parasite antigenic loci for a wide range of plant-parasite interactions. Long-term maintenance of polymorphism in GFG systems has remained puzzling for both theoreticians and empiricists. Traditionally this diversity has been explained by tradeoffs with other life-history traits closely linked with fitness, yet empirical evidence for such costs has remained mixed. Here we argue that incorporating simple ecological reality – spatial structuring and gradient of environmental conditions – into host–parasite research will help us understand how polymorphism is maintained. While environmental conditions (biotic and abiotic factors) have been studied in depth in plant pathology for their influence on disease severity and plant yield, they have been rarely set into an evolutionary framework. We briefly review recent data on natural plant–parasite metapopulations and theoretical models moving from single population models towards metapopulation theory to reveal in just how many ways spatial structuring may affect the coevolutionary process. We clarify also how spatially heterogeneous selection, through G×E (or G×G×E) interactions, may be particularly important for natural host–parasite interactions and suggest that this provides the unifying ground upon which future theoretical and empirical work should be build on.     

Detecting local adaptation in a natural plant-pathogen metapopulation: a laboratory vs. field transplant approach

Antagonistic coevolution between hosts and parasites in spatially structured populations can result in local adaptation of parasites. Traditionally parasite local adaptation has been investigated in field transplant experiments or in the laboratory under a constant environment. Despite the conceptual importance of local adaptation in studies of (co)evolution, to date no study has provided a comparative analysis of these two methods. Here, using information on pathogen population dynamics, I tested local adaptation of the specialist phytopathogen, Podosphaera plantaginis, to its host, Plantago lanceolata at three different spatial scales: sympatric host population, sympatric host metapopulation and allopatric host metapopulations. The experiment was carried out as a field transplant experiment with greenhouse-reared host plants from these three different origins introduced into four pathogen populations. In contrast to results of an earlier study performed with these same host and parasite populations under laboratory conditions, I did not find any evidence for parasite local adaptation. For interactions governed by strain-specific resistance, field studies may not be sensitive enough to detect mean parasite population virulence. Given that parasite transmission potential may be mediated by the abiotic environment and genotype-by-environment interactions, I suggest that relevant environmental variation should be incorporated into laboratory studies of parasite local adaptation.     

A synthesis of experimental work on parasite local adaptation

The study of parasite local adaptation, whereby parasites perform better on sympatric hosts than on allopatric hosts and/or better on their own host population than do other parasites, is of great importance to both basic and applied biology. Theoretical examination of host-parasite coevolution predicts that parasite migration rate, generation time and virulence all contribute to the pattern of parasite local adaptation, such that parasites with greater dispersal ability, more frequent reproduction and/or high virulence ought to exhibit increased infectivity on local hosts. Here, we present a meta-analysis of experimental work from 57 host-parasite systems across 54 local adaptation studies to directly test theoretical predictions concerning the effect of each attribute on parasite adaptation. As expected, we find that studies of parasites with higher migration rates than their hosts report local adaptation, as measured by infection success, significantly more often than studies of parasites with relatively low migration rates. Furthermore, this synthesis serves to identify biases in the current body of work and highlight areas with the greatest need for further study. We emphasize the importance of unifying the field with regard to experimental methods, local adaptation definitions and reported statistics for cross-infection studies.     

The impact of migration from parasite-free patches on antagonistic host-parasite coevolution

Natural populations of hosts and parasites are often subdivided and patchily distributed such that some regions of a host species' range will be free from a given parasite. Host migration from parasite-free to parasite-containing patches is expected to alter coevolutionary dynamics by changing the evolutionary potential of antagonists. Specifically, host immigration can favor parasites by increasing transmission opportunities, or hosts by introducing genetic variation. We tested these predictions in coevolving populations of Pseudomonas fluorescens and phage Phi2 that received immigrants from phage-free populations. We observed a negative quadratic relationship between sympatric resistance to phage and host immigration rate (highest at intermediate immigration) but a positive quadratic relationship between coevolution rate and host immigration rate (lowest at intermediate immigration). These results indicate that for a wide range of rates, host immigration from parasite-free patches can increase the evolutionary potential of parasites, and increase the coevolutionary rate if parasite adaptation is limiting in the absence of immigration.     

Superinfection and the coevolution of parasite virulence and host recovery

Parasite strategies of host exploitation may be affected by host defence strategies and multiple infections. In particular, within‐host competition between multiple parasite strains has been shown to select for higher virulence. However, little is known on how multiple infections could affect the coevolution between host recovery and parasite virulence. Here, we extend a coevolutionary model introduced by van Baalen (Proc. R. Soc. B, 265, 1998, 317) to account for superinfection. When the susceptibility to superinfection is low, we recover van Baalen's results and show that there are two potential evolutionary endpoints: one with avirulent parasites and poorly defended hosts, and another one with high virulence and high recovery. However, when the susceptibility to superinfection is above a threshold, the only possible evolutionary outcome is one with high virulence and high investment into defence. We also show that within‐host competition may select for lower host recovery, as a consequence of selection for more virulent strains. We discuss how different parasite and host strategies (superinfection facilitation, competitive exclusion) as well as demographic and environmental parameters, such as host fecundity or various costs of defence, may affect the interplay between multiple infections and host–parasite coevolution. Our model shows the interplay between coevolutionary dynamics and multiple infections may be affected by crucial mechanistic or ecological details.     

The evolution of sex and recombination in response to abiotic or coevolutionary fluctuations in epistasis

Evolutionary biologists have identified several factors that could explain the widespread phenomena of sex and recombination. One hypothesis is that host-parasite interactions favor sex and recombination because they favor the production of rare genotypes. A problem with many of the early models of this so-called Red Queen hypothesis is that several factors are acting together: directional selection, fluctuating epistasis, and drift. It is thus difficult to identify what exactly is selecting for sex in these models. Is one factor more important than the others or is it the synergistic action of these different factors that really matters? Here we focus on the analysis of a simple model with a single mechanism that might select for sex: fluctuating epistasis. We first analyze the evolution of sex and recombination when the temporal fluctuations are driven by the abiotic environment. We then analyze the evolution of sex and recombination in a two-species coevolutionary model, where directional selection is absent (allele frequencies remain fixed) and temporal variation in epistasis is induced by coevolution with the antagonist species. In both cases we contrast situations with weak and strong selection and derive the evolutionarily stable (ES) recombination rate. The ES recombination rate is most sensitive to the period of the cycles, which in turn depends on the strength of epistasis. In particular, more virulent parasites cause more rapid cycles and consequently increase the ES recombination rate of the host. Although the ES strategy is maximized at an intermediate period, some recombination is favored even when fluctuations are very slow. By contrast, the amplitude of the cycles has no effect on the ES level of sex and recombination, unless sex and recombination are costly, in which case higher-amplitude cycles allow the evolution of higher rates of sex and recombination. In the coevolutionary model, the amount of recombination in the interacting species also has a large effect on the ES, with evolution favoring higher rates of sex and recombination than in the interacting species. In general, the ES recombination rate is less than or equal to the recombination rate that would maximize mean fitness. We also discuss the effect of migration when sex and recombination evolve in a metapopulation. We find that intermediate parasite migration rates maximize the degree of local adaptation of the parasite and lead to a higher ES recombination rate in the host.     

Host-dependent genetic structure of parasite populations: differential dispersal of seabird tick host races

Abstract Despite the fact that parasite dispersal is likely to be one of the most important processes influencing the dynamics and coevolution of host-parasite interactions, little information is available on the factors that affect it. In most cases, opportunities for parasite dispersal should be closely linked to host biology. Here we use microsatellite genetic markers to compare the population structure and dispersal of two host races of the seabird tick Ixodes uriae at the scale of the North Atlantic. Interestingly, tick populations showed high within-population genetic variation and relatively low population differentiation. However, gene flow at different spatial scales seemed to depend on the host species exploited. The black-legged kittiwake ( Rissa tridactyla ) had structured tick populations showing patterns of isolation by distance, whereas tick populations of the Atlantic puffin ( Fratercula arctica ) were only weakly structured at the largest scale considered. Host-dependent rates of tick dispersal between colonies will alter infestation probabilities and local dynamics and may thus modify the adaptation potential of ticks to local hosts. Moreover, as I. uriae is a vector of the Lyme disease agent Borrelia burgdorferi sensu lato in both hemispheres, the large-scale movements of birds and the subsequent dispersal of ticks will have important consequences for the dynamics and coevolutionary interactions of this microparasite with its different vertebrate and invertebrate hosts.     

LOCAL MALADAPTATION IN THE ANTHER-SMUT FUNGUS MICROBOTRYUM VIOLACEUM TO ITS HOST PLANT SILENE LATIFOLIA: EVIDENCE FROM A CROSS-INOCULATION EXPERIMENT

Conventional wisdom holds that parasites evolve more rapidly than their hosts and are therefore locally adapted, that is, better at exploiting sympatric than allopatric hosts. We studied local adaptation in the insect-transmitted fungal pathogen Microbotryum violaceum and its host plant Silene latifolia. Infection success was tested in sympatric (local) and allopatric (foreign) combinations of pathogen and host from 14 natural populations from a metapopulation. Seedlings from up to 10 seed families from each population were exposed to sporidial suspensions from each of four fungal strains derived from the same population, from a near-by population (< 10 km distance), and from two populations at an intermediate (< 30 km) and remote (< 170 km) distance, respectively. We obtained significant pathogen X plant interactions in infection success (proportion of diseased plants) at both fungal population and strain level. There was an overall pattern of local maladaptation of this pathogen: average fungal infection success was significantly lower on sympatric hosts (mean proportion of diseased plants = 0.32 ± 0.03 SE) than on allopatric hosts (0.40 ± 0.02). Five of the 14 fungal populations showed no strong reduction in infection success on sympatric hosts, and three even tended to perform better on sympatric hosts. This pattern is consistent with models of time-lagged cycles predicting patterns of local adaptation in host-parasite systems to emerge only on average. Several factors may restrict the evolutionary potential of this pathogen relative to that of its host. First, a predominantly selfing breeding system may limit its ability to generate new virulence types by sexual recombination, whereas the obligately outcrossing host 5. latifolia may profit from rearrangement of resistance alleles by random mating. Second, populations often harbor only a few infected individuals, so virulence variation may be further reduced by drift. Third, migration rates among host plant populations are much higher than among pathogen populations, possibly because pollinators prefer healthy over diseased plants. Migration among partly isolated populations may therefore introduce novel host plant resistance variants more often than novel parasite virulence variants. That migration contributes to the coevolutionary dynamics in this system is supported by the geographic pattern of infectivity. Infection success increased over the first 10–km range of host-pathogen population distances, which is likely the natural range of gene exchange.     

THE IMPACT OF SPATIAL SCALE AND HABITAT CONFIGURATION ON PATTERNS OF TRAIT VARIATION AND LOCAL ADAPTATION IN A WILD PLANT PARASITE

Theory indicates that spatial scale and habitat configuration are fundamental for coevolutionary dynamics and how diversity is maintained in host–pathogen interactions. Yet, we lack empirical data to translate the theory to natural host–parasite systems. In this study, we conduct a multiscale cross‐inoculation study using the specialist wild plant pathogen Podosphaera plantaginis on its host plant Plantago lanceolata. We apply the same sampling scheme to a region with highly fragmented (Åland) and continuous (Saaremaa) host populations. Although theory predicts higher parasite virulence in continuous regions, we did not detect differences in traits conferring virulence among the regions. Patterns of adaptation were highly scale dependent. We detected parasite maladaptation among regions, and among populations separated by intermediate distances (6.0–40.0 km) within the fragmented region. In contrast, parasite performance did not vary significantly according to host origin in the continuous landscape. For both regions, differentiation among populations was much larger for genetic variation than for phenotypic variation, indicating balancing selection maintaining phenotypic variation within populations. Our findings illustrate the critical role of spatial scale and habitat configuration in driving host–parasite coevolution. The absence of more aggressive strains in the continuous landscape, in contrast to theoretical predictions, has major implications for long‐term decision making in conservation, agriculture, and public health.     

Differential virulence in a multiple‐host parasite of bumble bees: resolving the paradox of parasite survival?

Parasite virulence determines both the impact that parasites have on their hosts and parasite fitness. While most studies of virulence have involved single-species host–parasite interactions, the majority of parasites are likely to use multiple concurrent host species. Our understanding of how this impacts on parasite epidemiology and virulence is limited. Using the bumble bee Bombus lucorum , which exists in sympatry with B. terrestris in multi-species assemblages, and their generalist micosporidian parasite Nosema bombi , we tested whether the apparent paradox of parasite maintenance due to parasite virulence in a single host, B. terrestris , could be resolved through understanding the parasite's virulence in this sympatric host species. Nosema bombi significantly impacted colony growth, individual longevity, and individual development in B. lucorum . However, these effects were different both qualitatively and quantitatively to the parasite's impact in B. terrestris . Infected colonies of B. lucorum successfully produced both male and female reproductives, and infected female reproductives were capable of successful mating. Variation in life-history across host species may explain differences in the virulence, or impact of the parasite in B. terrestris and B. lucorum , with species with shorter life-cycles being more likely to transfer the parasite from one annual generation to the next. These results suggest that to understand the virulence and epidemiology of multi-host parasites we need to examine their ecological interactions across their various host species.     

The role of defensive symbionts in host–parasite coevolution

Understanding the coevolution of hosts and parasites is a long‐standing goal of evolutionary biology. There is a well‐developed theoretical framework to describe the evolution of host–parasite interactions under the assumption of direct, two‐species interactions, which can result in arms race dynamics or sustained genotype fluctuations driven by negative frequency dependence (Red Queen dynamics). However, many hosts rely on symbionts for defence against parasites. Whilst the ubiquity of defensive symbionts and their potential importance for disease control are increasingly recognized, there is still a gap in our understanding of how symbionts mediate or possibly take part in host–parasite coevolution. Herein we address this question by synthesizing information already available from theoretical and empirical studies. First, we briefly introduce current hypotheses on how defensive mutualisms evolved from more parasitic relationships and highlight exciting new experimental evidence showing that this can occur very rapidly. We go on to show that defensive symbionts influence virtually all important determinants of coevolutionary dynamics, namely the variation in host resistance available to selection by parasites, the specificity of host resistance, and the trade‐off structure between host resistance and other components of fitness. In light of these findings, we turn to the limited theory and experiments available for such three‐species interactions to assess the role of defensive symbionts in host–parasite coevolution. Specifically, we discuss under which conditions the defensive symbiont may take over from the host the reciprocal adaptation with parasites and undergo its own selection dynamics, thereby altering or relaxing selection on the hosts' own immune defences. Finally, we address potential effects of defensive symbionts on the evolution of parasite virulence. This is an important problem for which there is no single, clear‐cut prediction. The selection on parasite virulence resulting from the presence of defensive symbionts in their hosts will depend on the underlying mechanism of defence. We identify the evolutionary predictions for different functional categories of symbiont‐conferred resistance and we evaluate the empirical literature for supporting evidence. We end this review with outstanding questions and promising avenues for future research to improve our understanding of symbiont‐mediated coevolution between hosts and parasites.     

Rapid evolution and ecological host-parasite dynamics

Traditionally, the termination of parasite epidemics has been attributed to ecological causes: namely, the depletion of susceptible hosts as a result of mortality or acquired immunity. Here, we suggest that epidemics can also end because of rapid host evolution. Focusing on a particular host–parasite system, Daphnia dentifera and its parasite Metschnikowia bicuspidata , we show that Daphnia from lakes with recent epidemics were more resistant to infection and had less variance in susceptibility than Daphnia from lakes without recent epidemics. However, our studies revealed little evidence for genetic variation in infectivity or virulence in Metschnikowia . Incorporating the observed genetic variation in host susceptibility into an epidemiological model parameterized for this system reveals that rapid evolution can explain the termination of epidemics on time scales matching what occurs in lake populations. Thus, not only does our study provide rare evidence for parasite-mediated selection in natural populations, it also suggests that rapid evolution has important effects on short-term host–parasite dynamics.     

A meta-analysis of factors affecting local adaptation between interacting species

Adaptive divergence among populations can result in local adaptation, whereby genotypes in native environments exhibit greater fitness than genotypes in novel environments. A body of theory has developed that predicts how different species traits, such as rates of gene flow and generation times, influence local adaptation in coevolutionary species interactions. We used a meta-analysis of local-adaptation studies across a broad range of host-parasite interactions to evaluate predictions about the effect of species traits on local adaptation. We also evaluated how experimental design influences the outcome of local adaptation experiments. In reciprocally designed experiments, the relative gene flow rate of hosts versus parasites was the strongest predictor of local adaptation, with significant parasite local adaptation only in the studies in which parasites had greater gene flow rates than their hosts. When nonreciprocal studies were included in analyses, species traits did not explain significant variation in local adaptation, although the overall level of local adaptation observed was lower in the nonreciprocal than in the reciprocal studies. This formal meta-analysis across a diversity of host-parasite systems lends insight into the role of both biology (species traits) and biologists (experimental design) in detecting local adaptation in coevolving species interactions.