Host-parasite coevolution can lead to local adaptation of either parasite or host if there is specificity (GxG interactions) and asymmetric evolutionary potential between host and parasite. This has been demonstrated both experimentally and in field studies, but a substantial proportion of studies fail to detect such clear-cut patterns. One explanation for this is that adaptation can be masked by counter-adaptation by the antagonist. Additionally, genetic architecture underlying the interaction is often highly complex thus preventing specific adaptive responses. Here, we have employed a reciprocal cross-infection experiment to unravel the adaptive responses of two components of fitness affecting both parties with different complexities of the underlying genetic architecture (i.e. mortality and spore load). Furthermore, our experimental coevolution of hosts (Tribolium castaneum) and parasites (Nosema whitei) included paired replicates of naive hosts from identical genetic backgrounds to allow separation between host- and parasite-specific responses.
In hosts, coevolution led to higher resistance and altered resistance profiles compared to paired control lines. Host genotype × parasite genotype interactions (GH × GP) were observed for spore load (the trait of lower genetic complexity), but not for mortality. Overall parasite performance correlated with resistance of its matching host coevolution background reflecting a directional and unspecific response to strength of selection during coevolution. Despite high selective pressures exerted by the obligatory killing parasite, and host- and parasite-specific mortality profiles, no general pattern of local adaptation was observed, but one case of parasite maladaptation was consistently observed on both coevolved and control host populations. In addition, the use of replicate control host populations in the assay revealed one case of host maladaptation and one case of parasite adaptation that was masked by host counter-adaptation, suggesting the presence of complex and probably dynamically changing fitness landscapes.
Our results demonstrate that the use of replicate naive populations can be a useful tool to differentiate between host and parasite adaptation in complex and dynamic fitness landscapes. The absence of clear local adaptation patterns during coevolution with a sexual host showing a complex genetic architecture for resistance suggests that directional selection for generality may be more important attributes of host-parasite coevolution than commonly assumed.
The antagonistic arms race between hosts and parasites is one of the most evolutionarily significant biotic interaction . Parasite species might represent more than half of the known biodiversity  and this ubiquity paired with detrimental fitness effects, can affect host population dynamics [3, 4], genetic diversity , biodiversity, ecosystem functioning and community structure . Consequently, parasites play a key role in various evolutionary and ecological theories, and are thought to be a major factor explaining the existence of male display in birds , polyandry , and the evolution of sex , despite the evolutionary costs that are associated with these phenomena.
Generally speaking, it is assumed that parasites have a larger evolutionary potential, due to their larger population size, shorter generation time, and higher mutation rates . Yet, sexual hosts with complex genetic architecture underlying resistance can prevent adaptation of (asexual) parasites by producing more heterogeneous offspring, which in turn could lead to the evolution of generalist parasites [11, 12]. Indeed, studies reporting experimental evolution using (facultatively) sexual hosts have shown parasite maladaptation  or a mosaic of patterns, with "no adaptation" being the most common finding . This suggests that sexual hosts might stay ahead in the coevolutionary game by producing genetically diverse offspring, especially when genetic architecture of resistance traits are complex.
In spatially structured populations, sub-populations may experience different and asynchronous coevolutionary trajectories, which could lead to rapid between-population divergence [15, 16]. Given a genetic basis of infectivity and resistance  the dynamic nature of coevolutionary interactions should temporarily lead to local adaptation of one or both antagonists [10, 18–24]. As a consequence of local adaptation of one antagonist, local maladaptation of the other (Table 1) is also commonly observed [25–27] and is indeed predicted by the geographic mosaic theory of coevolution [28, 29]. Few studies using multicellular hosts have been able to separate the two , as ideally, identical but naive host and/or parasite populations are needed to unequivocally interpret results obtained with cross infection studies (Table 1).
Possible outcomes (rows, columns) of a cross-infection experiment where parasite performance is assayed, both, on coevolved host populations, and on control host populations paired for the same original host lines.
Control host populations
Mortality in matching combinations > non-matching combinations
Mortality in matching combinations < non-matching combinations
No difference between matching and non-matching combinations
Mortality in matching combinations > non-matching combinations
1: Parasite adaptation
2: Parasite maladaptation < host maladaptation
3: Host maladaptation
Coevolved host populations
Mortality in matching combinations < non-matching combinations
4: Parasite adaptation < host adaptation
5: Parasite maladaptation
6: Host adaptation
No difference between matching and non-matching combinations
7: Parasite adaptation = host adaptation
8: Parasite maladaptation = host maladaptation
9: No adaptation
The numbered entries refer to the interpretation for each possible combination of outcomes a. A matching combination refers to hosts and parasite being from the same experimental coevolution replicate; non-matching otherwise.
a Differentiating between parasite adaptation and host maladaptation can be facilitated if parasite fitness is assayed on replicate host populations of identical genetic background as their coevolving antagonists. In this table we have used host mortality as a measure of both parasite and host performance, as these traits are very closely linked to fitness in both antagonists in our study system [39, 41]. For other study systems, relevant traits may differ. This experimental approach has, to our knowledge, never been used, but could most readily be performed in a laboratory setting.
One requirement for local adaptation in such cross-infection studies is that the outcome of exposure relies on strong host genotype (GH) × parasite genotype (GP) interactions (GH xGP). These interactions are typically suggestive of only a few underlying genes and low genetic diversity in infecting inocula . However, evidence for a more complex genetic architecture underlying resistance is accumulating [30, 31] and in natural populations, multiple infections seem to be the rule rather than the exception . Both processes are expected to lead to less noticeable GH × GP interactions. The importance of genetic architecture is clearly illustrated by genome-wide association studies which repeatedly demonstrate that seemingly highly heritable human diseases can insufficiently be explained by genetic markers, even when taking into account the joint effects of many tens of loci [33, 34]. As theory predicts that the extent to which antagonists are able to adapt decreases with increasing number of loci involved in the interaction , this added complexity could explain why a substantial proportion of studies fail to demonstrate clear-cut patterns of local adaptation . Consequently theory predicts that in a population where host genetic variation in resistance is large, parasites may be trapped in the middle of the host phenotypic distribution, thereby leading to a reduction in parasite variance, and a generalist strategy .
We have already tested several predictions of the coevolutionary theory using experimental coevolution of Tribolium castaneum and its natural, obligately killing parasite, Nosema whitei  demonstrating coevolutionary change of both parties within local demes by time shift experiments [38, 39], which represents a crucial prerequisite for local adaptation. Here, by using hosts and parasites from the same long-term study as described in references [38–40], we now tried to assess in more detail whether the reciprocal phenotypic changes led to local adaptation/specificity of either antagonist. In detail we had the following objectives: First of all, we determined whether coevolution actually leads to global and/or population-specific local adaptation patterns by performing a fully reciprocal cross infection experiment. Then, to determine the influence of trait specific genetic architecture we focused on two fitness measures with different degrees of complexity in genetic architecture : (i) Host mortality (55 days post exposure) is a proxy both for host fitness and parasite fitness, as parasite transmission is only possible after host death [41, 42], and shows a complex genetic architecture; and (ii) spore load, which is a measure of both transmission potential of the parasite and host resistance [39, 41] and shows a rather simple genetic architecture. By cross-infecting replicate host lines from both the coevolved treatment and their paired control lines from the same genetic background kept under parasite free conditions, we increase the power to discriminate between adaptive responses of the parasite and the host, respectively. For example, if adaptation of a parasite to its own host line is neutralized by counter-adaptation of its matched coevolved host line, the parasite's adaptive response could still be detected when exposed to their matched naïve host lines, which were kept in identical but parasite-free conditions (see Table 1 for all the possible outcomes of a cross-infection experiment as performed in this paper).
If selection on parasite infectivity leads to the evolution of generality in host use, similarity in performance between parasite isolates should reflect the respective selective background and correlate with similarity in resistance between matching host lines, because resistant host lines should select for more infective/virulent parasites. This would be compatible with a form of coevolution where genetic variance of sexually reproducing hosts in combination with the complexity of host genetic architecture can prevent specific adaptation of the parasite. Complementing the longitudinal trajectory of host-parasite coevolution in this system  with the spatial dimension of local adaptation will provide a better understanding of the evolution of specificity as well as the selective landscape of host-parasite coevolution.
During the experiment, larvae raised in medium with parasites had significantly higher mortality rates than those raised in control (parasite-free) medium (29% vs 3%, Treatment: F6,1531 = 16.7, P < 0.001), regardless of their selection regime (Control: 38% vs 3%; contrast for control vs parasite sources: z-value = -5.04, P < 0.001; Coevolved: 20% vs 3%; contrast for control vs parasite sources: z-value = -3.67, P < 0.001). Furthermore, mortality induced by parasites that had coevolved with any of the host lines (i.e. averaged over parasites from all coevolved lines pooled) was lower than mortality induced by ancestral parasites (25.4 ± 10.2% S.E. vs. 44.9 ± 17.7; contrast analysis, z-value = 6.27, P < 0.001). Similarly, ancestral parasites achieved higher spore loads than evolved parasites (20'040 ± 3'554 S.E. vs. 8'436 ± 990 S.E., contrast analysis: z-value = 4.473, P < 0.001).
Main effects, GHOST x GPARASITE interactions and the response to selection
Analyzing the cross-infection experiment for contemporary combinations (i.e. from the same generation) separately, strong main effects of host line and parasite isolate on host mortality were found (Figure 1, Table 2). Coevolved host lines were more resistant to N. whitei, but evolved resistance depended on the parasite isolate used (see Parasite × Selection interaction term in Table 2). Host × parasite interactions differed between coevolved and control selection regime, which is an expected result of antagonistic coevolution (Figure 1, three-way interaction in Table 2), but no significant overall host × parasite interactions were found (Figure 1, Table 2).
Results of generalized linear model of host mortality after exposure to N.whitei using binomial error distribution. A)
Results of post hoc contrast analysis of mortality
Within host (coevolved lines)
Within parasite (coevolved lines)
Within host (control lines)
Within parasite (control lines)
Levels of significance for the GLM-model fits were tested using analysis of deviance with chi-square distribution. Post hoc test results for contrasts between sympatric and allopatric combinations are shown below. Positive Z values indicate combinations where local antagonists show higher mortality than foreign antagonists.
Mortality of the beetles collected (in blocks) for measuring spore loads correlated with mortality in the survival experiment (Spearman rank correlation, r = 0.81, P < 0.001). Not surprisingly, both mortality in the experiment and sampled mortality correlated with spore load of the randomly collected subsample (Spearman rank correlation, r = 0.74, P < 0.001, and r = 0.89, P < 0.001 respectively, Figure 2a).
When analyzing the experiment using the coevolved parasites only, spore load differed mainly between beetles that had already died or were still alive at time of collection (Table 3). Additionally, spore load differed between host lines and parasite isolates, and was generally lower in coevolved hosts, but there was also a significant host line × parasite isolate interaction (Figure 3a, Figure 3b, Table 3). The difference in spore load between coevolved and control host lines depended on the parasite isolate, but also on different host line × parasite isolate interactions under each selection regime (Figures 3a and 3b, Table 3).
Results of a generalized linear model of spore load in a randomly collected subsample of eight beetles per experimental block.
Individual mortality (Dead/Alive)
Results of post hoc contrast analysis of spore load
Within host (coevolved lines)
Within parasite (coevolved lines)
Within host (control lines)
Within parasite (control lines)
Levels of significance of GLM model fits were tested using analysis of deviance with chi-square distribution. Post hoc test results for contrasts between sympatric and allopatric combinations are shown below. Positive Z values indicate combinations where were local antagonists show higher spore load than foreign antagonists.
Any evidence for local adaptation?
When analyzing mortality, no global local parasite adaptation (in the sense ofhigher mortality on "own" hosts) was observed in either the coevolved (contrast analysis, z = 0.33, P = 0.932) or control regime (z = -0.029, P = 0.999, Figure 1c). When testing for population-specific host or parasite adaptation, the majority of cases revealed no significant adaptation (8 out of 10 comparisons in both the coevolved and control selection regime, Table 2). Within-parasite line tests for adaptation revealed no significant adaptation while within-host line tests revealed that 4 out of 10 combinations showed significant adaptation or maladaptation (Table 2). Both parasite adaptation (better with "own" host, positive Z-value, Table 2) and maladaptation (better with "foreign" hosts, negative Z-value, Table 2) were observed when assayed on host lines from either selection regime (Table 2).
Although sporeload measurements suggested local adaptation of the parasite (Figure 3c), we found no overall significant difference between own and foreign combinations in either the coevolved (z = 0.27, P = 0.954) or control treatment (z = -0.23, P = 0.967, Figure 3c). However, unlike host mortality, parasite adaptation with respect to spore load was only observed when parasites were exposed to control host lines (parasite isolate 4, Table 3). Of the nine possible scenarios of local adaptation (Table 1), we found evidence for scenario 7 (parasite adaptation = host adaptation) in line 4 when looking at mortality and spore load (Table 2 Table 3), scenario 5 (parasite maladaptation) in line 1 when looking at mortality, and scenario 3 (host maladaptation) in line 6 for mortality. In the remaining parasite and host lines no form of local adaptation could be detected.
Correlation between host susceptibility and matching parasite virulence
Phenotypic differences in virulence (defined as parasite-exposed host mortality) between parasite isolates correlated with differentiation in resistance of its matching coevolved hosts in the coevolved hosts (Mantel test, r = 0.74, P = 0.009, Figure 2b). Such a correlation was absent when parasites were assayed on the control host lines (Mantel test, r = 0.16, P = 0.308). Parasite population differentiation in spore load did not correlate with host population differentiation in spore load (Mantel test, coevolved regime: r = -0.53, P = 0.976; control regime: r = -0.11, P = 0.575), possibly reflecting the stronger underlying host × parasite interactions.
Variation in mortality among host lines upon exposure to parasite isolates did not differ significantly between control and coevolution treatment (mean coefficient of variance: control: 63.8, coevolved: 68.6, pairwise t-test t = 1.33, df = 4, P = 0.253). Similarly, variation in induced host mortality among parasite isolates did not differ between selection regimes (mean coefficient of variance: control: 47.4, coevolved: 55.9; pairwise t-test, t = 0.65, df = 4, P = 0.551). Nevertheless, host and parasite population differentiation in performance were smaller in the coevolved selection regime than under control conditions (pairwise t-test on pairwise distances, parasite performance: t = -3.26, df = 9, P = 0.009; host resistance: t = -3.57, df = 9, P = 0.005) indicating a flattening of the fitness landscape during coevolution.
A general lack of local adaptation after experimental coevolution
As proposed by Schulte et al.  a distinction needs to be made between "mosaic adaptation" (where only some host or parasite populations show local adaptation, while other populations show no adaptation or even maladaptation) and "local adaptation" (where hosts or parasite populations show local adaptation across a range of tested populations). Our results provide evidence for mosaic adaptation, indicating that adaptation and coevolution occurs, but we currently have no evidence for local adaptation in our experimental lines. Note that while parasite maladaptation can often be ascribed to host (counter-) adaptation, our experimental set up provides the means to separate these two explanatory mechanisms. Based on the predicted possible scenarios, one case where local combinations showed lower mortality could be ascribed to local parasite maladaptation, and one case of higher local mortality could be ascribed to host maladaptation (Table 2). One parasite isolate showed consistent adaptation to its host genetic background for both measured traits, which was not observed when looking at the coevolved combinations separately (Table 2, Table 3). In summary, our results are in concordance with the geographic mosaic theory of coevolution  which predicts that local adaptation may vary between demes, both spatially and temporarily, leading to a lack of an overall pattern of local adaptation. Additionally, in two out of five populations, no adaptation was detected, suggesting that these populations may be so-called "cold spots" [15, 28].
The effect of complexity of genetic architecture on GHOST x GPARASITE interactions
A substantial body of theoretical and empirical coevolutionary work is based on either the implicit or explicit assumption of high specificity in host-parasite interactions, which should "ideally" lead to negative frequency-dependent selection and spatio-temporal adaptation of antagonists [20, 43, 44]. Typically host × parasite (GH × GP) interactions occur in systems where only a few major genes and low genetic diversity in infecting inocula are involved in the compatibility . However, more complex genetic architecture with epistatic interactions between many genes do not necessarily lead to such easily interpretable outcomes . Despite the strong positive correlation between mortality and spore load, significant GHOST x GPARASITE interactions were indeed found for spore load but no such interactions were detected for mortality (Table 2, Table 3). As genetic architecture underlying mortality is more complex than that of spore load, (as host mortality is most likely a composite trait, with possible contributions of both resistance and tolerance, ) this can potentially explain the difference in GH × GP interactions between our fitness measures.
Relative importance of specific and unspecific adaptation
Our results support other laboratory-based coevolutionary experiments that demonstrate very little evidence of specific local adaptation [13, 14, 46]. Pouillan et al  showed that if phage Φ2 was allowed to adapt to Pseudomonas fluorescens SBW25 specific adaptation to host genetic background was observed, while coevolution did not lead to adaptation, but rather to a broader infectivity range. This may indicate that the increased genotypic variance created by rapid reciprocal coevolutionary changes may preclude adaptation, and that local adaptation detected by field observations may be the result of directional selection for generality rather than evidence of underlying coevolution [20, 47, 48].
There are four additional factors that may explain the absence of local adaptation in this experiment. First, it has been shown that an increase in the number of sympatric units dramatically decreases the chance to detect local adaptation, but increases the probability to find host or parasite main effects . And indeed, a substantial amount of variation in outcome of exposure to N. whitei was explained by such main effects in both parties (Table 1). Second, temporal oscillations in parasite infectivity potentially mask local adaptation, especially when these oscillations are not synchronized between populations . Although such oscillations initially occur in this system, they dampened rapidly over time . Therefore this seems not a crucial factor explaining the lack of local adaptation at the time point when this experiment was carried out. Third, a recent meta-analyses  showed that forces generating diversity such as mutation, migration are essential for local adaptation to occur. In the absence of such forces, erosion of standing genetic variation can decrease evolutionary potential of antagonists. In our experiment, the sexual and recombining hosts may be at an advantage compared with the asexual parasites. As no migration was allowed in our experiment, lineage sorting could have led to depletion of clonal diversity, thus precluding parasite adaptation. Finally, and most importantly, theoretical work suggests that if resistance is polygenic, such as the case for the T.castaneum - N. whitei interaction [31, 51], selection may favour generalist parasites. Under such a scenario, within-host genetic variance is expected to increase during coevolution, while within-parasite genetic variance decreases, due to its dependency on the mean of the host phenotype distribution .
A loss of parasite performance due to relaxed selection?
The predominant lack of (overall and population specific) local adaptation cannot be interpreted as a lack of coevolutionary responses of the parasites. Even though at the start of the experiment all lines were inoculated with the same spore cocktail , parasites rapidly differentiated in performance, while on average inducing lower host mortality upon exposure than the ancestral parasites. The null hypothesis is that diversification can be a result of genetic drift due to repeated bottlenecking. But diversification seems to be driven by the selective environment of their coevolving host genetic background, as parasite performance correlated strongly with host resistance (Figure 2A) meaning that coevolution on a low resistance background may be equal to relaxed selection. Our results thus demonstrate that susceptibility of a given host line can select for general attenuation of the parasite, a pattern that differs from parasite maladaptation . Therefore we can show that coevolution may lead to adaptive changes and divergence between parasite isolates without leading to a pattern of general local adaptation. This leads to two important evolutionary implications. First, the genetic background of the host and its evolutionary potential seems crucial in determining the evolutionary trajectory of the host-parasite system. It also shows the importance of a directional selection component for generality during antagonistic coevolution giving support to previous studies in other systems [52–56]. Second, the observed loss of virulence under relaxed selection (i.e. a low resistance background) indicates that virulence may be costly in this system, and thus can be selected against if selection on growth rate is higher than selection on killing rate . That differences in resistance can affect the evolution of virulence, additionally, confirms both theoretical and empirical work [57, 58]. However, the correlation between resistance and virulence only manifested in parasites infecting the coevolved host lines. This different pattern of mortality caused by parasite exposure in the two selection regimes is likely a result of the weaker host × parasite interactions in coevolved hosts (see significant GH × GP x Selection interaction, Table 1).
We presented evidence for rapid differentiation in infection profiles between parasite populations during antagonistic coevolution, which is likely to be driven by differences in selection coefficients imposed by differing resistance levels of coevolving host lines. At the same time, coevolution has led to increased host resistance and smaller differences in resistance profile between host lines. We did, however, not find general local adaptation in either host or parasite. By comparing parasite performance on coevolved and paired control populations we could nevertheless differentiate between adaptation and maladaptation of the antagonist and could observe likely cases of parasite adaptation, parasite maladaptation, and host maladaptation. Our results may suggest that, in this system, evolution towards generalist exploitation of hosts with complex genetic architectures is more likely than the evolution towards increased specificity, thereby explaining the lack of clear-cut patterns of local adaptation.
All hosts and parasites used in this study originated from a coevolution experiment that has been continuously running for over three years prior to the infection experiment presented here, using Tribolium castaneum and its natural, microsporidian parasite Nosema whitei, of which the protocol is described in more detail in ref. . A total of five random populations were chosen as our experimental lines for this cross-infection experiment; in particular, we used lines nr.1, 3, 4, 5 and 6 from the study of ref. . All lines had a different genetic background, and thus by implication differing levels of initial resistance. At the start of the selection experiment (before generation 1), every F1 hybrid line was divided in half, and each half was subjected to one of two selection regimes. In the "coevolution" regime, lines were subjected to coevolution with the Nosema whitei, which led to mortalities of up to 40% . In the "control" regime lines of identical origin and genetic background were maintained in the absence of parasites. In this way, half of the beetles of line 3, for example, were assigned to coevolution and the other half of line 3 to control - and so on for every line (representing the different genetic backgrounds). Population size was kept constant at 500 adult beetles in both selection regimes, by always collecting this number as breeders to initiate each following generation. The ancestral parasite used to inoculate the host lines at the start of the experiment was the same for all host lines, and consisted of a mixture (in equal proportions) of eight Nosema whitei isolates . The regimes were maintained for a total of 16 generations since the start of the selection experiment, which has been described in more detail in . After 16 generations we relaxed selection for one generation to start the cross-infection assay. Relaxed selection was achieved by rearing all lines from both treatments under control conditions (i.e. with standard flour and environmental conditions, and in the absence of parasites) to avoid potential trans-generational effects of exposure to parasites on resistance of assayed individuals in the infection experiment [59, 60].
Infection and survival experiment
In the following generation 17, after one generation of relaxation, we collected 500 unsexed adult beetles from each of these experimental lines to produce offspring that were later used in the cross-infection experiment. The resulting larval offspring were assigned randomly to the following treatments: (i) Exposure to the ancestral parasites (the same mixture of spores that was used to inoculate the coevolution experiment in generation 0); (ii) Exposure to spores from all 5 coevolved host lines collected after 17 generations of coevolution. All five host lines were subjected to all five parasite isolates in a full factorial design, meaning that for each host line there was one "sympatric" combination and four "allopatric" combinations, while each individual larvae still was only exposed to one parasite source; (iii) Controls that were not exposed to parasite spores. Freshly hatched larvae (1-2 days old) were randomly collected from a single jar for each host line, and were subsequently placed individually into glass vials (13 × 40 mm), containing 0.1 g of either parasite-inoculated flour (5 × 104 spores/gram) or parasite-free medium. A total of 1543 larvae were successfully distributed (5 replicate lines × 7 treatments × 2 selection regimes × ca. 25 larvae each). After their assignment to a treatment the larvae were kept under standard environmental conditions (24 h dark, 33°C, 70% humidity), and vials were checked for survival when the experiment was terminated (55 d after the assignment of larvae).
Spore load measurement
From each experimental block (replicate host line × selection regime × infection treatment) we collected 8 random beetles for N. whitei spore load analysis. Spore load was measured using quantitative real time PCR of a 220 bp product of N. whitei 16sRNA using methods as described in . DNA was extracted using 96-well plate extraction kits (Qiagen), quantified using a Nanodrop 8000 (Thermo scientific), and diluted to 5 ng μl-1. Spore load was quantified in twofold for each sample, and quantification of samples with a difference of more than 1 CT between replicates were repeated, but excluded from further analysis if the difference remained higher than 1 CT. To quantify spores we used a duplicated fourfold dilution series of the same ten standard samples in every run. The highest dilution representing the detection threshold was set to the equivalent of one N. whitei spore. Negative controls, using ddH20 instead of host DNA, were additionally used in every run.
To test if exposure to parasites caused higher host mortality as compared to control conditions, we used a Generalized linear mixed model (glmmPQL from MASS package in R), with individual mortality as response, treatment (exposure/control) as fixed factor, and host line and selection regime as random factors. To test whether host mortality (binary trait: alive/dead) when exposed to coevolved parasites differed from exposure to ancestral parasites we used contrast analysis on the results of the model above. Individual survival (binary trait: alive/dead) was subsequently analyzed for all contemporary combinations using a GLM, with the factors host line, parasite source and selection regime, but all possible interactions were kept in the model. To test for directional selection on infectivity, we used Mantel tests to analyze the association between coevolved host resistance (which dictate the selection pressures on the parasites) and virulence of matching parasites when assayed on both coevolved and host lines.
To test if parasite-induced host mortality is predictive of transmission potential (spore load), a Spearman's rank correlation on mean values of each host-parasite combination was used. Spore load was subsequently analyzed for all contemporary combinations, using a GLM-model with spore load as response (quasipoisson error distribution), and individual mortality (dead/alive) as a fixed factor. Additional fixed factors were host line, parasite isolate, selection regime, with all interactions kept in the model. To compare means of total mortality and spore load of interest, we used multiple contrast analysis. Levels of significance for the GLM models were tested with analysis of deviance following a chi-square distribution. All statistical analyses were conducted with the R statistical package .
The authors thank Daniel Trujillo-Vallegas, Natasha Rossel and Miguel Hon for help with cleaning and autoclaving everything needed for this experiment, Tilmann Silber for his RT-PCR work, and EO and TB groups for helpful feedback during discussions. We also thank Andrew Morgan, Oliver Otti, Michael Brockhurst and three anonymous referees for constructive feedback on an earlier draft of the manuscript. Supported by the Genetic Diversity Center of ETH Zurich (GDC) and CCES. Financially supported by SNF grant 31-120451 to KMW and ETH grant nr. TH-09 60-1 to PSH.
Experimental Ecology, Institute of Integrative Biology
Institute of Evolutionary Biology, University of Edinburgh
Evolutionary Ecology of Marine Fishes, Leibniz Institute for Marine Sciences (IfM-Geomar)
Alfred Wegener Institute for Polar and Marine Science, Wadden Sea Station Sylt
Schmid-Hempel P: Evolutionary parasitology: the integrated study of infections, immunology, ecology and genetics. Oxford: Oxford University Press; 2011.
Toft CA: Communities of parasites with parasitic life-styles. In Community ecology. Edited by: Diamond JM. Case TJ: Harper & Row; 1986:445–463.
Little TJ: The evolutionary significance of parasitism: do parasite-driven genetic dynamics occur ex silico?Journal of Evolutionary Biology 2002,15(1):1–9.View Article
Hudson PJ, Dobson AP, Newborn D: Prevention of population cycles by parasite removal.Science 1998,282(5397):2256–2258.PubMedView Article
Haldane JBS: Disease and evolution.La Ricerca Scientifica 1949,19(suppl):68–76.
Hudson PJ, Dobson AP, Lafferty KD: Is a healthy ecosystem one that is rich in parasites?Trends in Ecology & Evolution 2006,21(7):381–385.View Article
Hamilton WD, Zuk M: Heritable true fitness and bright birds: a role for parasites.Science 1982,218(4570):384–387.PubMedView Article
Baer B, Schmid-Hempel P: Experimental variation in polyandry affects parasite loads and fitness in a bumble-bee.Nature 1999,397(6715):151–154.View Article
Hamilton WD: Sex versus non-sex versus parasite.Oikos 1980,35(2):282–290.View Article
Greischar MA, Koskella B: A synthesis of experimental work on parasite local adaptation.Ecology Letters 2007,10(5):418–434.PubMedView Article
Legros M, Koella JC: Experimental evolution of specialization by a microsporidian parasite.Bmc Evolutionary Biology 2010., 10:
Ebert D, Hamilton WD: Sex against virulence: the coevolution of parasitic diseases.Trends in Ecology & Evolution 1996,11(2):79–82.View Article
Adiba S, Huet M, Kaltz O: Experimental evolution of local parasite maladaptation.Journal of Evolutionary Biology 2010,23(6):1195–1205.PubMedView Article
Schulte RD, Makus C, Hasert B, Michiels NK, Schulenburg H: Host-parasite local adaptation after experimental coevolution ofCaenorhabditis elegansand its microparasiteBacillus thuringiensis.Proceedings of the Royal Society B: Biological Sciences 2011.
Thompson JN: Specific hypotheses on the geographic mosaic of coevolution.American Naturalist 1999, 153:S1-S14.View Article
Buckling A, Rainey PB: The role of parasites in sympatric and allopatric host diversification.Nature 2002,420(6915):496–499.PubMedView Article
Luijckx P, Ben-Ami F, Mouton L, Du Pasquier L, Ebert D: Cloning of the unculturable parasite Pasteuria ramosa and its Daphnia host reveals extreme genotype-genotype interactions.Ecology Letters 2011,14(2):125–131.PubMedView Article
Ebert D: Virulence and local adaptation of a horizontally transmitted parasite.Science 1994,265(5175):1084–1086.PubMedView Article
Kaltz O, Shykoff JA: Local adaptation in host-parasite systems.Heredity 1998, 81:361–370.View Article
Lively CM, Dybdahl MF: Parasite adaptation to locally common host genotypes.Nature 2000,405(6787):679–681.PubMedView Article
Hatcher MJ, Hogg JC, Dunn AM: Local adaptation and enhanced virulence ofNosema granulosisartificially introduced into novel populations of its crustacean host,Gammarus duebeni.International Journal for Parasitology 2005,35(3):265–274.PubMedView Article
Laine AL: Spatial scale of local adaptation in a plant-pathogen metapopulation.Journal of Evolutionary Biology 2005,18(4):930–938.PubMedView Article
Refardt D, Ebert D: Inference of parasite local adaptation using two different fitness components.Journal of Evolutionary Biology 2007,20(3):921–929.PubMedView Article
Morgan AD, Gandon S, Buckling A: The effect of migration on local adaptation in a coevolving host-parasite system.Nature 2005,437(7056):253–256.PubMedView Article
Kaltz O, Gandon S, Michalakis Y, Shykoff JA: Local maladaptation in the anther-smut fungusMicrobotryum violaceumto its host plantSilene latifolia: Evidence from a cross-inoculation experiment.Evolution 1999,53(2):395–407.View Article
Morand S, Manning SD, Woolhouse MEJ: Parasite-host coevolution and geographic patterns of parasite infectivity and host susceptibility.Proceedings of the Royal Society of London Series B-Biological Sciences 1996,263(1366):119–128.View Article
Thompson JN: The coevolutionary process. The University of Chicago Press; 1994.
Nuismer SL, Thompson JN, Gomulkiewicz R: Gene flow and geographically structured coevolution.Proceedings of the Royal Society of London Series B-Biological Sciences 1999,266(1419):605–609.View Article
Wilfert L, Schmid-Hempel P: The genetic architecture of susceptibility to parasites.Bmc Evolutionary Biology 2008., 8:
Wegner KM, Berenos C, Schmid-Hempel P: Nonadditive genetic components in resistance of the red flour beetleTribolium castanaeumagainst parasite infection.Evolution 2008,62(9):2381–2392.PubMedView Article
Read AF, Taylor LH: The ecology of genetically diverse infections.Science 2001,292(5519):1099–1102.PubMedView Article
Gibson G: Hints of hidden heritability in GWAS.Nature Genetics 2010,42(7):558–560.PubMedView Article
Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al.: Finding the missing heritability of complex diseases.Nature 2009,461(7265):747–753.PubMedView Article
Ridenhour BJ, Nuismer SL: Polygenic traits and parasite local adaptation.Evolution 2007,61(2):368–376.PubMedView Article
Kopp M, Gavrilets S: Multilocus genetics and the coevolution of quantitative traits.Evolution 2006,60(7):1321–1336.PubMed
Burges HD, Weiser J: Occurrence of pathogens of the flour beetle,Tribolium castaneum.Journal of Invertebrate Pathology 1973,22(3):464–466.View Article
Berenos C, Schmid-Hempel P, Wegner KM: Experimental coevolution leads to a decrease in parasite-induced host mortality.Journal of Evolutionary Biology 2011,24(8):1777–1782.PubMedView Article
Berenos C, Schmid-Hempel P, Wegner KM: Evolution of host resistance and trade-offs between virulence and transmission potential in an obligately killing parasite.Journal of Evolutionary Biology 2009,22(10):2049–2056.PubMedView Article
Berenos C, Wegner KM, Schmid-Hempel P: Antagonistic coevolution with parasites maintains host genetic diversity: an experimental test.Proceedings of the Royal Society B-Biological Sciences 2011,278(1703):218–224.View Article
Blaser M, Schmid-Hempel P: Determinants of virulence for the parasiteNosema whiteiin its hostTribolium castaneum.Journal of Invertebrate Pathology 2005,89(3):251–257.PubMedView Article
Milner RJ: Nosema whitei, a Microsporidan pathogen of some species ofTribolium:1: morphology, life-cycle, and generation time.Journal of Invertebrate Pathology 1972,19(2):231–238.View Article
Peters AD, Lively CM: The Red Queen and fluctuating epistasis: A population genetic analysis of antagonistic coevolution.American Naturalist 1999,154(4):393–405.PubMedView Article
Decaestecker E, Gaba S, Raeymaekers JAM, Stoks R, Van Kerckhoven L, Ebert D, De Meester L: Host-parasite 'Red Queen' dynamics archived in pond sediment.Nature 2007, 450:870-U816.PubMedView Article
Lambrechts L: Dissecting the genetic architecture of host-pathogen specificity.PLoS Pathogens 2010,6(8):e1001019.PubMedView Article
Poullain V, Gandon S, Brockhurst MA, Buckling A, Hochberg ME: The evolution of specificity in evolving and coevolving antagonistic interactions between a bacteria and its phage.Evolution 2008,62(1):1–11.PubMed
Imhoof B, Schmid-Hempel P: Patterns of local adaptation of a protozoan parasite to its bumblebee host.Oikos 1998,82(1):59–65.View Article
Hatcher M, Hogg J, Dunn A: Local adaptation and enhanced virulence of artificially introduced into novel populations of its crustacean host.International Journal for Parasitology 2005,35(3):265–274.PubMedView Article
Koskella B, Lively CM: Advice of the rose: Experimental coevolution of a trematode parasite and its snail host.Evolution 2007,61(1):152–159.PubMedView Article
Hoeksema JD, Forde SE: A meta-analysis of factors affecting local adaptation between interacting species.American Naturalist 2008,171(3):275–290.PubMedView Article
Wegner KM, Berenos C, Schmid-Hempel P: Host genetic architecture in single and multiple infections.Journal of Evolutionary Biology 2009,22(2):396–404.PubMedView Article
Schulte RD, Makus C, Hasert B, Michiels NK, Schulenburg H: Multiple reciprocal adaptations and rapid genetic change upon experimental coevolution of an animal host and its microbial parasite.Proceedings of the National Academy of Sciences 2010,107(16):7359–7364.View Article
Duncan AB, Mitchell SE, Little TJ: Parasite-mediated selection and the role of sex and diapause inDaphnia.Journal of Evolutionary Biology 2006,19(4):1183–1189.PubMedView Article
Buckling A, Rainey PB: Antagonistic coevolution between a bacterium and a bacteriophage.Proceedings of the Royal Society of London Series B-Biological Sciences 2002,269(1494):931–936.View Article
Lohse K, Gutierrez A, Kaltz O: Experimental evolution of resistance inParamecium caudatumagainst the bacterial parasiteHolospora undulata.Evolution 2006,60(6):1177–1186.PubMed
Mackinnon MJ, Read AF: Immunity promotes virulence evolution in a malaria model.Plos Biology 2004,2(9):1286–1292.View Article
Gandon S, Michalakis Y: Evolution of parasite virulence against qualitative or quantitative host resistance.Proceedings of the Royal Society of London Series B-Biological Sciences 2000,267(1447):985–990.View Article
Sadd BM, Kleinlogel Y, Schmid-Hempel R, Schmid-Hempel P: Trans-generational immune priming in a social insect.Biology Letters 2005,1(4):386–388.PubMedView Article
Roth O, Joop G, Eggert H, Hilbert J, Daniel J, Schmid-Hempel P, Kurtz J: Paternally derived immune priming for offspring in the red flour beetle,Tribolium castaneum.Journal of Animal Ecology 2009,79(2):403–413.PubMedView Article
R_Development_Core_Team: R: A language and environment for statistical computing.R Foundation for Statistical Computing Vienna 2010.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.