Skip to main content
  • Research article
  • Open access
  • Published:

Complex adaptive responses during antagonistic coevolution between Tribolium castaneum and its natural parasite Nosema whitei revealed by multiple fitness components



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 [1]. Parasite species might represent more than half of the known biodiversity [2] and this ubiquity paired with detrimental fitness effects, can affect host population dynamics [3, 4], genetic diversity [5], biodiversity, ecosystem functioning and community structure [6]. 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 [7], polyandry [8], and the evolution of sex [9], 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 [10]. 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 [13] or a mosaic of patterns, with "no adaptation" being the most common finding [14]. 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 [17] the dynamic nature of coevolutionary interactions should temporarily lead to local adaptation of one or both antagonists [10, 1824]. As a consequence of local adaptation of one antagonist, local maladaptation of the other (Table 1) is also commonly observed [2527] 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 [14], as ideally, identical but naive host and/or parasite populations are needed to unequivocally interpret results obtained with cross infection studies (Table 1).

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.

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 [17]. 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 [32]. 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 [35], this added complexity could explain why a substantial proportion of studies fail to demonstrate clear-cut patterns of local adaptation [10]. 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 [36].

We have already tested several predictions of the coevolutionary theory using experimental coevolution of Tribolium castaneum and its natural, obligately killing parasite, Nosema whitei [37] 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 [3840], 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 [31]: (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 [38] 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).

Figure 1
figure 1

Heatmap of mortality for all host-parasite combinations for A) the control and B) coevolved hosts separately. Shades of red indicate the observed mortality, with darker shades corresponding to higher mortality (see legend). C) Barplot of mortality when exposed to own parasites (dark grey bars) and foreign parasites (light grey bars). Corresponding statistical details can be found in Table 1. Error bars denote ± 1 S.E.

Table 2 Results of generalized linear model of host mortality after exposure to N.whitei using binomial error distribution. A)

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).

Figure 2
figure 2

Correlations between fitness measures. A) Correlation between mean mortality in the eight randomly sampled beetles in each experimental block, and mean spore load as a response variable. The solid lines shows the best fitting model (R2 = 0.80, F2,57 = 117.01, P < 0.001). Overall, there was a significant correlation between spore load and sampled mortality (Spearman rank correlation, r = 0.89, P < 0.001). B) Association of host (virulence) and parasite (failure due to resistance or infectivity) mortalities. The pairwise differences between host populations (virulence) are plotted against the respective pairwise differences in the parasite populations. Parasite population differences correlated with coevolved host differences, but only for coevolved (broken line; Mantel test, Coevolved hosts: r = 0.74, P = 0.009) and not for control hosts (r = 0.16, P = 0.308). Shown is the best fitting linear model for the coevolved hosts. s

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).

Table 3 Results of a generalized linear model of spore load in a randomly collected subsample of eight beetles per experimental block.
Figure 3
figure 3

Heatmap of residual spore load controlled for mortality for all host-parasite combinations for A) the control and B) coevolved hosts separately. Shades of red indicate residual spore load, with lighter shades corresponding to higher residual spore load (see legend). C) Bar charts of residual spore load when exposed to own parasites (dark grey bars) and foreign parasites (light grey bars). Corresponding statistical details can be found in Table 2. Error bars denote ± 1 S.E.

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. [14] 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 [15] 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 [17]. However, more complex genetic architecture with epistatic interactions between many genes do not necessarily lead to such easily interpretable outcomes [45]. 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, [31]) 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 [46] 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 [10]. 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 [49]. Although such oscillations initially occur in this system, they dampened rapidly over time [38]. 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 [50] 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 [36].

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 [39], 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 [13]. 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 [5256]. 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 [39]. 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.


Selection experiment

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. [39]. 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. [39]. 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% [39]. 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 [39]. 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 [39]. 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 [31]. 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.

Statistical analyses

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 [61].


  1. Schmid-Hempel P: Evolutionary parasitology: the integrated study of infections, immunology, ecology and genetics. 2011, Oxford: Oxford University Press

    Google Scholar 

  2. Toft CA: Communities of parasites with parasitic life-styles. Community ecology. Edited by: Diamond JM. 1986, Case TJ: Harper & Row, 445-463.

    Google Scholar 

  3. Little TJ: The evolutionary significance of parasitism: do parasite-driven genetic dynamics occur ex silico?. Journal of Evolutionary Biology. 2002, 15 (1): 1-9. 10.1046/j.1420-9101.2002.00366.x.

    Article  Google Scholar 

  4. Hudson PJ, Dobson AP, Newborn D: Prevention of population cycles by parasite removal. Science. 1998, 282 (5397): 2256-2258.

    Article  PubMed  CAS  Google Scholar 

  5. Haldane JBS: Disease and evolution. La Ricerca Scientifica. 1949, 19 (suppl): 68-76.

    Google Scholar 

  6. 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. 10.1016/j.tree.2006.04.007.

    Article  Google Scholar 

  7. Hamilton WD, Zuk M: Heritable true fitness and bright birds: a role for parasites. Science. 1982, 218 (4570): 384-387. 10.1126/science.7123238.

    Article  PubMed  CAS  Google Scholar 

  8. Baer B, Schmid-Hempel P: Experimental variation in polyandry affects parasite loads and fitness in a bumble-bee. Nature. 1999, 397 (6715): 151-154. 10.1038/16451.

    Article  CAS  Google Scholar 

  9. Hamilton WD: Sex versus non-sex versus parasite. Oikos. 1980, 35 (2): 282-290. 10.2307/3544435.

    Article  Google Scholar 

  10. Greischar MA, Koskella B: A synthesis of experimental work on parasite local adaptation. Ecology Letters. 2007, 10 (5): 418-434. 10.1111/j.1461-0248.2007.01028.x.

    Article  PubMed  Google Scholar 

  11. Legros M, Koella JC: Experimental evolution of specialization by a microsporidian parasite. Bmc Evolutionary Biology. 2010, 10:

    Google Scholar 

  12. Ebert D, Hamilton WD: Sex against virulence: the coevolution of parasitic diseases. Trends in Ecology & Evolution. 1996, 11 (2): 79-82. 10.1016/0169-5347(96)81047-0.

    Article  CAS  Google Scholar 

  13. Adiba S, Huet M, Kaltz O: Experimental evolution of local parasite maladaptation. Journal of Evolutionary Biology. 2010, 23 (6): 1195-1205. 10.1111/j.1420-9101.2010.01985.x.

    Article  PubMed  CAS  Google Scholar 

  14. Schulte RD, Makus C, Hasert B, Michiels NK, Schulenburg H: Host-parasite local adaptation after experimental coevolution of Caenorhabditis elegans and its microparasite Bacillus thuringiensis. Proceedings of the Royal Society B: Biological Sciences. 2011

    Google Scholar 

  15. Thompson JN: Specific hypotheses on the geographic mosaic of coevolution. American Naturalist. 1999, 153: S1-S14. 10.1086/303208.

    Article  Google Scholar 

  16. Buckling A, Rainey PB: The role of parasites in sympatric and allopatric host diversification. Nature. 2002, 420 (6915): 496-499. 10.1038/nature01164.

    Article  PubMed  CAS  Google Scholar 

  17. 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. 10.1111/j.1461-0248.2010.01561.x.

    Article  PubMed  Google Scholar 

  18. Ebert D: Virulence and local adaptation of a horizontally transmitted parasite. Science. 1994, 265 (5175): 1084-1086. 10.1126/science.265.5175.1084.

    Article  PubMed  CAS  Google Scholar 

  19. Kaltz O, Shykoff JA: Local adaptation in host-parasite systems. Heredity. 1998, 81: 361-370. 10.1046/j.1365-2540.1998.00435.x.

    Article  Google Scholar 

  20. Lively CM, Dybdahl MF: Parasite adaptation to locally common host genotypes. Nature. 2000, 405 (6787): 679-681. 10.1038/35015069.

    Article  PubMed  CAS  Google Scholar 

  21. Hatcher MJ, Hogg JC, Dunn AM: Local adaptation and enhanced virulence of Nosema granulosis artificially introduced into novel populations of its crustacean host, Gammarus duebeni. International Journal for Parasitology. 2005, 35 (3): 265-274. 10.1016/j.ijpara.2004.12.004.

    Article  PubMed  Google Scholar 

  22. Laine AL: Spatial scale of local adaptation in a plant-pathogen metapopulation. Journal of Evolutionary Biology. 2005, 18 (4): 930-938. 10.1111/j.1420-9101.2005.00933.x.

    Article  PubMed  Google Scholar 

  23. Refardt D, Ebert D: Inference of parasite local adaptation using two different fitness components. Journal of Evolutionary Biology. 2007, 20 (3): 921-929. 10.1111/j.1420-9101.2007.01307.x.

    Article  PubMed  CAS  Google Scholar 

  24. 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. 10.1038/nature03913.

    Article  PubMed  CAS  Google Scholar 

  25. Kaltz O, Gandon S, Michalakis Y, Shykoff JA: Local maladaptation in the anther-smut fungus Microbotryum violaceum to its host plant Silene latifolia: Evidence from a cross-inoculation experiment. Evolution. 1999, 53 (2): 395-407. 10.2307/2640776.

    Article  Google Scholar 

  26. 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. 10.1098/rspb.1996.0019.

    Article  CAS  Google Scholar 

  27. Thompson JN, Nuismer SL, Gomulkiewicz R: Coevolution and maladaptation. Integr Comp Biol. 2002, 42 (2): 381-387. 10.1093/icb/42.2.381.

    Article  PubMed  Google Scholar 

  28. Thompson JN: The coevolutionary process. 1994, The University of Chicago Press

    Chapter  Google Scholar 

  29. 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. 10.1098/rspb.1999.0679.

    Article  Google Scholar 

  30. Wilfert L, Schmid-Hempel P: The genetic architecture of susceptibility to parasites. Bmc Evolutionary Biology. 2008, 8:

    Google Scholar 

  31. Wegner KM, Berenos C, Schmid-Hempel P: Nonadditive genetic components in resistance of the red flour beetle Tribolium castanaeum against parasite infection. Evolution. 2008, 62 (9): 2381-2392. 10.1111/j.1558-5646.2008.00444.x.

    Article  PubMed  CAS  Google Scholar 

  32. Read AF, Taylor LH: The ecology of genetically diverse infections. Science. 2001, 292 (5519): 1099-1102. 10.1126/science.1059410.

    Article  PubMed  CAS  Google Scholar 

  33. Gibson G: Hints of hidden heritability in GWAS. Nature Genetics. 2010, 42 (7): 558-560. 10.1038/ng0710-558.

    Article  PubMed  CAS  Google Scholar 

  34. 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. 10.1038/nature08494.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  35. Ridenhour BJ, Nuismer SL: Polygenic traits and parasite local adaptation. Evolution. 2007, 61 (2): 368-376. 10.1111/j.1558-5646.2007.00029.x.

    Article  PubMed  Google Scholar 

  36. Kopp M, Gavrilets S: Multilocus genetics and the coevolution of quantitative traits. Evolution. 2006, 60 (7): 1321-1336.

    Article  PubMed  CAS  Google Scholar 

  37. Burges HD, Weiser J: Occurrence of pathogens of the flour beetle, Tribolium castaneum. Journal of Invertebrate Pathology. 1973, 22 (3): 464-466. 10.1016/0022-2011(73)90178-X.

    Article  Google Scholar 

  38. 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. 10.1111/j.1420-9101.2011.02306.x.

    Article  PubMed  CAS  Google Scholar 

  39. 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. 10.1111/j.1420-9101.2009.01821.x.

    Article  PubMed  CAS  Google Scholar 

  40. 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. 10.1098/rspb.2010.1211.

    Article  PubMed Central  Google Scholar 

  41. Blaser M, Schmid-Hempel P: Determinants of virulence for the parasite Nosema whitei in its host Tribolium castaneum. Journal of Invertebrate Pathology. 2005, 89 (3): 251-257. 10.1016/j.jip.2005.04.004.

    Article  PubMed  Google Scholar 

  42. Milner RJ: Nosema whitei, a Microsporidan pathogen of some species of Tribolium :1: morphology, life-cycle, and generation time. Journal of Invertebrate Pathology. 1972, 19 (2): 231-238. 10.1016/0022-2011(72)90214-5.

    Article  Google Scholar 

  43. Peters AD, Lively CM: The Red Queen and fluctuating epistasis: A population genetic analysis of antagonistic coevolution. American Naturalist. 1999, 154 (4): 393-405. 10.1086/303247.

    Article  PubMed  Google Scholar 

  44. 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. 10.1038/nature06291.

    Article  PubMed  CAS  Google Scholar 

  45. Lambrechts L: Dissecting the genetic architecture of host-pathogen specificity. PLoS Pathogens. 2010, 6 (8): e1001019-10.1371/journal.ppat.1001019.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 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  Google Scholar 

  47. Imhoof B, Schmid-Hempel P: Patterns of local adaptation of a protozoan parasite to its bumblebee host. Oikos. 1998, 82 (1): 59-65. 10.2307/3546917.

    Article  Google Scholar 

  48. 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. 10.1016/j.ijpara.2004.12.004.

    Article  PubMed  Google Scholar 

  49. Koskella B, Lively CM: Advice of the rose: Experimental coevolution of a trematode parasite and its snail host. Evolution. 2007, 61 (1): 152-159. 10.1111/j.1558-5646.2007.00012.x.

    Article  PubMed  Google Scholar 

  50. Hoeksema JD, Forde SE: A meta-analysis of factors affecting local adaptation between interacting species. American Naturalist. 2008, 171 (3): 275-290. 10.1086/527496.

    Article  PubMed  Google Scholar 

  51. Wegner KM, Berenos C, Schmid-Hempel P: Host genetic architecture in single and multiple infections. Journal of Evolutionary Biology. 2009, 22 (2): 396-404. 10.1111/j.1420-9101.2008.01657.x.

    Article  PubMed  CAS  Google Scholar 

  52. 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. 10.1073/pnas.1003113107.

    Article  CAS  Google Scholar 

  53. Duffy MA, Sivars-Becker L: Rapid evolution and ecological host-parasite dynamics. Ecology Letters. 2007, 10 (1): 44-53. 10.1111/j.1461-0248.2006.00995.x.

    Article  PubMed  Google Scholar 

  54. Duncan AB, Mitchell SE, Little TJ: Parasite-mediated selection and the role of sex and diapause in Daphnia. Journal of Evolutionary Biology. 2006, 19 (4): 1183-1189. 10.1111/j.1420-9101.2006.01085.x.

    Article  PubMed  CAS  Google Scholar 

  55. 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. 10.1098/rspb.2001.1945.

    Article  Google Scholar 

  56. Lohse K, Gutierrez A, Kaltz O: Experimental evolution of resistance in Paramecium caudatum against the bacterial parasite Holospora undulata. Evolution. 2006, 60 (6): 1177-1186.

    Article  PubMed  Google Scholar 

  57. Mackinnon MJ, Read AF: Immunity promotes virulence evolution in a malaria model. Plos Biology. 2004, 2 (9): 1286-1292.

    Article  CAS  Google Scholar 

  58. 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. 10.1098/rspb.2000.1100.

    Article  CAS  Google Scholar 

  59. 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. 10.1098/rsbl.2005.0369.

    Article  PubMed  PubMed Central  Google Scholar 

  60. 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.

    Article  PubMed  Google Scholar 

  61. R_Development_Core_Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing Vienna. 2010

    Google Scholar 

Download references


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.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Camillo Bérénos.

Additional information

Authors' contributions

CB designed and performed the experiment, analyzed the data and wrote the manuscript. PSH and KMW designed the experiment and wrote the manuscript. All authors have approved the final version.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Bérénos, C., Schmid-Hempel, P. & Wegner, K.M. Complex adaptive responses during antagonistic coevolution between Tribolium castaneum and its natural parasite Nosema whitei revealed by multiple fitness components. BMC Evol Biol 12, 11 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: