Joy Bose
Rebecca D. Schulte- Department of Behavioral Biology, University of Osnabrueck, Osnabrueck, Germany
Host–parasite interactions represent one of the strongest selection pressures in nature. They are often governed by genotype-specific (GxG) interactions resulting in host genotypes that differ in resistance and parasite genotypes that differ in virulence depending on the antagonist’s genotype. Another type of GxG interactions, which is often neglected but which certainly influences host–parasite interactions, are those between coinfecting parasite genotypes. Mechanistically, within-host parasite interactions may range from competition for limited host resources to cooperation for more efficient host exploitation. The exact type of interaction, i.e., whether competitive or cooperative, is known to affect life-history traits such as virulence. However, the latter has been shown for chosen genotype combinations only, not considering whether the specific genotype combination per se may influence the interaction (i.e., GxG interactions). Here, we want to test for the presence of GxG interactions between coinfections of the bacterium Bacillus thuringiensis infecting the nematode Caenorhabditis elegans by combining two non-pathogenic and five pathogenic strains in all possible ways. Furthermore, we evaluate whether the type of interaction, reflected by the direction of virulence change of multiple compared to single infections, is genotype-specific. Generally, we found no indication for GxG interactions between non-pathogenic and pathogenic bacterial strains, indicating that virulence of pathogenic strains is equally affected by both non-pathogenic strains. Specific genotype combinations, however, differ in the strength of virulence change, indicating that the interaction type between coinfecting parasite strains and thus the virulence mechanism is specific for different genotype combinations. Such interactions are expected to influence host–parasite interactions and to have strong implications for coevolution.
Introduction
Genotype by genotype (GxG) interactions are interactions between species which depend on the species’ genotypes. Thus, the exact phenotype of the interaction depends on the genotypes that are involved (Hamilton, 1980; Thompson, 2005). Such interactions can be mutualistic like cooperation (Hoeksema and Thompson, 2007; Heath, 2010; Gorton et al., 2012) or antagonistic like parasitism (Webster and Woolhouse, 1998; Carius et al., 2001; Schmid-Hempel and Ebert, 2003; Schulenburg and Ewbank, 2004; Grech et al., 2006; Lambrechts et al., 2006, 2009; Rauch et al., 2006; Salvaudon et al., 2007; de Roode and Altizer, 2010; Luijckx et al., 2011). They are one of the fundamental requirements for coevolution, as shown for the coevolution between hosts and parasites (Hamilton et al., 1990; Agrawal and Lively, 2002).
Host–parasite GxG interactions are not static. They can be influenced by abiotic factors like temperature or nutrition (Thomas and Blanford, 2003; Lazzaro and Little, 2009; Wolinska and King, 2009; Sadd, 2011; Tollenaere and Laine, 2012; Gsell et al., 2013; Mboup et al., 2013), but also biotic factors, namely a third species. Facultative symbiotic microbes have been shown to influence or even mediate GxG interactions between parasites and hosts (Oliver et al., 2005; Tetard-Jones et al., 2007; Koch and Schmid-Hempel, 2012; Rouchet and Vorburger, 2012). Recently, Seppälä et al. (2012) revealed GxG interactions between coinfecting parasite species that are furthermore influenced by the environmental factor dose. Coinfections are likely to also influence GxG interactions between host and parasite, resulting in rather complex interactions (i.e., GxGxGxE). Next to coinfections of different parasite species, coinfections of different genotypes of the same parasite species are frequently observed in nature. GxG interactions for the latter remain to be shown.
Although knowledge on the role of GxG interactions between coinfecting parasites is scarce, the consequences of two parasites simultaneously infecting the same host (i.e., coinfections) have been studied in detail. Such coinfections may either be antagonistic, like resource competition (faster but not optimal host exploitation is selected at the within-host level) or spiteful competition (competitors are directly fought even though it is costly for the actor), or they may be mutualistic like public good cooperation (produced goods can be used by all individuals causing a fitness advantage for both the producer and the recipient; Buckling and Brockhurst, 2008; Mideo, 2009). It has furthermore been described that the host immune response may interfere with the interaction between competing parasites. It may do so by either preventing infection of one species if another is present or by facilitating its infection (Cox, 2001; Read and Taylor, 2001; Mideo, 2009). Such scenarios are likely to reveal GxGxG interactions between the host and the coinfecting parasites: depending on the host genotype, infection by certain parasite genotypes may be selectively favored over other parasite genotypes.
Studies on multiple infections usually assume that one of the potential interaction mechanisms is at work. We argue that depending on the exact parasite genotypes coinfecting a host, the mechanism might be different. Most theoretical and empirical studies on multiple infections focused on parasite virulence (e.g., Alizon et al., 2013), and virulence is expected to increase (e.g., resource competition) or decrease (e.g., spiteful competition) in multiple compared to single infections depending on the exact type of interaction. Therefore, we consider the presence of multiple infections with lower and with higher virulence compared to the mean of the corresponding single infections, according to the exact parasite genotype combination tested, as an indication for such mechanisms.
Here, we used the bacterial microparasite Bacillus thuringiensis (Bt) and its nematode host Caenorhabditis elegans. This model system has been well described (Wei et al., 2003; Schulenburg and Müller, 2004; Schulte et al., 2010). Bt causes persistent gut infections in C. elegans that potentially lead to host death (Borgonie et al., 1995, 1996; Schulenburg and Müller, 2004). The infection is caused by oral uptake of bacterial spores during feeding. The spores are associated with crystal toxins (Cry-toxins), which destroy intestinal cells. Cry-toxins are thought to be the prime determinant for the infection (Griffitts and Aroian, 2005), but also other virulence factors like phospholipase C, proteinases and hemolysins have been described (George and Crickmore, 2012). Once host resources are made available, the spores germinate and proliferate vegetatively until cells sporulate. Both antagonists show potential for specific interactions: Bt strains show high specificity against nematodes, including C. elegans (Wei et al., 2003; Schulenburg and Ewbank, 2004; Schulenburg and Müller, 2004; Schulte et al., 2010), and C. elegans expresses specific immune reactions toward different pathogens (Alper et al., 2007; Wong et al., 2007; Schulenburg et al., 2008).
Bacillus thuringiensis has the potential for within-host interaction since different strains differ in their growth rates and produce different public Cry-toxins to exploit the host (Payne, 1992; Payne et al., 1993). Bt is furthermore capable of bacteriocin production, substances harming other Bt-strains (Abriouel et al., 2011), which therefore can be classified as spiteful behavior. Thus, different strains may interact differently.
Since we were especially interested in the within-host interactions, we used the same outbred host population and all possible combinations of single and double infections of two non-pathogenic and five pathogenic Bt strains. As phenotypic proxies for virulence we measured host survival and host reproduction. If parasite genotypes do not interact, multiple infections should take an intermediate value of the two corresponding single infections. Interaction however may result in reduced or decreased virulence, depending on the exact type of interaction.
Our aims for this study were (i) to test for GxG interactions between coinfecting parasite strains, i.e., whether the virulence in a coinfection is influenced by the infecting genotypes, (ii) to elucidate whether there is a general difference between single and double infections which would indicate that most interactions between coinfecting parasites are of the same type, (iii) to test whether the exact change in virulence between single and mixed infections depends on the coinfecting genotypes, indicating that the exact interaction mechanism between coinfecting genotypes is genotype specific, and (iv) to find indication for the exact interaction mechanism for specific genotype combinations.
Materials and Methods
Bacterial Strains and Nematodes
In total, we used seven different Bt strains. The nematocidal strains B-18243, B-18245, B-18246, B-18247, and B-18679 were provided by the Agricultural Research Service Patent Culture Collection (United States Department of Agriculture). Different Bt-strains differ in their genotype and Cry-toxin production (Payne, 1992; Payne et al., 1993; Schnepf et al., 2001), as do the strains we used here (Schulte et al., 2010, H. Schulenburg, personal communication). The other two strains were non-nematocidal, thus non-pathogenic toward nematodes, namely DSM-350 (German Collection of Microorganisms and Cell Cultures) and 407 Cry- (kindly provided by Christina Nielsen-LeRoux; Lereclus et al., 1989). Prior to the experiments, all Bt strains were cultured in large quantities and stored in aliquots at -20°C (Borgonie et al., 1995).
As nematode host we used a genetically diverse and outcrossed population of C. elegans to simulate natural conditions. Thus, our results are likely to be valid for C. elegans as a species and not only for a specific genotype. The population was originally prepared by Henrique Teotónio by consecutive crosses among 16 natural isolates (Teotónio et al., 2012). At least 13 of these isolates are genetically diverse and the outcrossed population is more diverse than natural populations (Rockman and Kruglyak, 2009; Teotónio et al., 2012). This population was adapted over 10 generations in 40 replicates to our laboratory conditions. Afterward the replicates were mixed, aliquoted and cryopreserved at -80°C (Stiernagle, 2006) for later usage. Otherwise, worms were maintained following standard procedures (Stiernagle, 2006).
Experimental Design and Protocol
To compare single and double infections, we tested all Bt strains in all possible combinations (Figure 1). Thus, we had five different treatments: non-pathogenic single Bt, non-pathogenic mixed Bt, non-pathogenic-pathogenic mixed Bt, pathogenic single Bt and pathogenic mixed Bt. In mixed infection treatments, both strains were mixed in equal proportions. Importantly, the total Bt concentration was identical in all treatments. The whole experiment was replicated four times. One replicate of the combination between 407 Cry- and B-18247 for survival was lost during experimental procedure.
FIGURE 1. Schematic overview of the experimental design showing different treatments and different strain combinations within each treatment. Every strain combination was replicated four times, except for survival of the combination of 407 Cry- and B-18247 (three times).
The experiment was performed at 18°C and at 70% humidity. 20 age synchronized hermaphroditic worms of the last larval stage (L4) were transferred into a “worm-ball” (Sicard et al., 2007) containing peptone-free nematode-growth medium. Each worm-ball was inoculated with 3.5 × 107 Bt particles and ad libitum Escherichia coli OP50 as a food source to prevent any effects caused by starvation. Beforehand, Bt was grown for 6 days on nematode growth medium. After 3 days, nematode survival and population size were estimated as proxies for virulence. Survival rate was measured as the number of surviving worms divided by the sum of surviving and dead worms. Population size was estimated as the number of worms per ball by washing them off with sterile water, counting twice a subsample of 20 μl and extrapolating the total number of worms. For statistical analysis of population size, we used the logarithm to the base 10.
Statistical Analysis
All statistical analyses were performed using IBM SPSS Statistics Version 20.
To test for GxG interactions between coinfecting non-pathogenic and pathogenic Bt strains, we used a generalized linear model with “non-pathogenic strain,” “pathogenic strain” and the interaction term as fixed factors and survival rate and population size as dependent variables. GxG interactions are revealed by a significant interaction term. Since such a model based analysis is based on a full-factorial design, this analysis is not possible for the mixed pathogenic treatment (the five pathogenic strains cannot be grouped unequivocally into two factors which needed to be tested in all combinations; Figure 1).
To elucidate whether there is a general difference between single and double infections, we analyzed a general linear mixed model with “treatment” (non-pathogenic-pathogenic mixed Bt, pathogenic single Bt and pathogenic mixed Bt) as fixed factor, “combination (treatment)” as random factor and survival rate and population size as dependent variables. The exact difference between the three treatment categories was tested using LSD (Fisher’s Least Significance Difference) as post hoc test.
Our main question is, however, whether the change in virulence between single and multiple infections and thus the type of interaction between coinfecting strains depends on their exact genotypes. This should become visible by plotting single and multiple infections of the mixed treatments. For each mixed infection, we calculated the mean from the corresponding single infections. For this, we created random pairs of the four replicates (e.g., the mean of single infection A and B of replicate 1 was calculated). We estimated the influence of the factors multitude of infection (mean of single versus mixed infections) and strain combination (10 combinations for the non-pathogenic-pathogenic mixed and the pathogenic mixed treatment, each) and the interaction between both as fixed factors in a generalized linear model. A significant interaction term reveals in how far the difference between the mixed infection and mean single infections depended on the genotype combination tested. This analysis was performed for the mixed non-pathogenic-pathogenic and the mixed pathogenic treatments separately.
Finally, we were interested in the exact type of interaction between coinfecting strains. Depending on the interaction type, mixed infections should be more or less virulent than the mean of the corresponding single infections. Thus, we tested whether the slopes between the two infection types of each strain combination differed from zero using a 1-sample t-test.
Results
We did not find any indications for GxG interactions in coinfections of non-pathogenic and pathogenic strains as indicated by non-significant interaction terms (Table 1; Figure 2). However, non-pathogenic strains differed in their effect on host survival and pathogenic strains tend to differ.
TABLE 1. Statistical results of the generalized linear model testing for GxG interactions in the non-pathogenic-pathogenic mixed treatment using non-pathogenic Bt strain, pathogenic Bt strain and the interaction between both as fixed factors.
FIGURE 2. Line graphs of (A) survival rate and (B) logarithmic population size of the non-pathogenic-pathogenic mixed treatment. Vertical lines represent the one-fold standard error.
Generally, infection treatments differed from each other (Table 2; Figure 3). The mixed non-pathogenic-pathogenic treatment differed from the pathogenic treatments for both variables, while the two pathogenic treatments marginally differed for population size, but not for survival rate.
TABLE 2. Statistical results of the comparison of the three pathogenic treatments (mixed non-pathogenic-pathogenic Bt, pure pathogenic Bt, mixed pathogenic Bt) using a general linear mixed model with treatment as fixed factor and strain combination(treatment) as random factor and a post hoc test (LSD).
FIGURE 3. Boxplots of (A) mean host survival rate and (B) mean logarithmic population size per strain combination of different treatments. The horizontal line indicates the median, the box the 25% quartile above and below the median, and the whiskers the data range. Circles represent outliers and stars represent extreme values.
Furthermore, we found indications that the interaction mechanism depends on the coinfecting parasite genotypes (Table 3). In how far the mixed treatments differ from the mean of the single values for survival depends on the exact genotype combination tested, as indicated by a trend for the interaction term for the non-pathogenic-pathogenic mixed treatment (Figure 4A) and a significant interaction term for the pathogenic mixed treatment (Figure 4B). This was not the case for population size (Figures 4C,D).
TABLE 3. Statistical results of the generalized linear model testing for a genotype-specific difference between single and multiple infections using multitude of infection (mean of single infections versus mixed infection), strain combination and the interaction between both as fixed factors.
FIGURE 4. Line graphs of mixed treatments and the mean of the corresponding two single infections. (A) Survival rate of the non-pathogenic-pathogenic mixed treatment. (B) Survival rate of the pathogenic mixed treatment. (C) Logarithmic population size of the non-pathogenic-pathogenic mixed treatment. (D) Logarithmic population size of the pathogenic mixed treatment. Dashed lines represent combinations, for which the slope differs significantly from zero, dashed-dotted lines represent combinations, for which the slope shows a tendency to differ from zero. Vertical lines represent the onefold standard error.
The significant interaction terms are not only caused by deviations in one direction. We found that for survival, the slopes of the strain combinations deviate in both directions from zero. For the non-pathogenic-pathogenic mixed treatment, we found one case in which the mix is more virulent than the mean of the single infections (407 Cry- and B-18246: 1-sample t-test, t = 8.653, df = 3, p = 0.003; Figure 4A) and one combination in which it is less virulent than the mean of the two single infections (407 Cry- and B-18243: 1-sample t-test, t = -15.236, df = 3, p = 0.001; Figure 4A). For the pathogenic mixed treatment, the mix is more virulent in one combination (B-18243 and B-18247: 1-sample t-test, t = 5.644, df = 3, p = 0.011; Figure 4B) and less in two combinations (B-18246 and B-18247: 1-sample t-test, t = -9.2724, df = 3, p = 0.003; B-18246 and B-18679: 1-sample t-test, t = -3.342, df = 3, p = 0.044; Figure 4B). For host population size, one combination of the non-pathogenic-pathogenic mixed treatment performed worse in mixed infections (407 Cry- and B-18247: 1-sample t-test, t = -3.669, df = 3, p = 0.035; Figure 4C). One combination of the pathogenic mixed treatment did show higher virulence in the mixed compared to the single treatments (B-18245 and B-18247: 1-sample t-test, t = 4.948, df = 3, p = 0.016; Figure 4D), two showed a trend (B-18243 and B-18245: 1-sample t-test, t = 2.950, df = 3, p = 0.060; B-18245 and B-18246: 1-sample t-test, t = 2.759, df = 3, p = 0.070; Figure 4D).
Discussion
Theoretical models on multiple infections usually assume only one mechanism of parasite interaction to be present, although different interaction types like resource competition or public good cooperation have been described (Buckling and Brockhurst, 2008; Mideo, 2009). Furthermore, GxG interactions between coinfecting parasites occur (Seppälä et al., 2012), which may influence the parasite interaction. Here, our aim was to test whether the interaction mechanism can be determined by GxG interactions between coinfecting parasites.
We could not show that the virulence of coinfections depends on GxG interactions between both infecting strains. Yet, we cannot exclude their presence since we were only able to test for such interactions in the non-pathogenic-pathogenic mixed treatment, and not in the pathogenic mixed treatment. The latter might have been more suitable for such an analysis since both and not only one strain are pathogenic and reduce fitness. We also only tested a rather small set of non-pathogenic strains in a low sample size, which might have obscured GxG interactions.
However, we provide evidence that within-host interactions of Bt genotypes follow a general pattern: Over all tested genotype combinations, we found that mixed non-pathogenic-pathogenic combinations are less virulent than pathogenic combinations (Figure 3). This may be explained by the infection characteristics of this model system (Schnepf et al., 1998): Bt spores, which are taken up by the host are already associated with Cry-toxins. These proteins are of prime importance for the infection, since they cause pore formation in the gut epithelium. They are released and activated in the gut, thus representing a public good (Raymond et al., 2007). In mixed non-pathogenic-pathogenic combinations, only one strain produces toxins and therefore virulence is reduced. What is known about Cry-toxin genes is therefore in line with the theory on public good cooperation and non-cooperating cheaters (Brown et al., 2002; Buckling and Brockhurst, 2008; Alizon and Lion, 2011).
Although we were not able to show that coinfections are characterized by GxG interactions, the virulence difference between multiple and the corresponding single infections depends on the exact genotype combination that coinfects a host (Figure 4). For some genotype combinations, virulence is higher, for others it is lower in double infections compared to the mean single infections. This indicates that different interaction mechanisms may be involved in the Bt-C. elegans model system. For example, next to the production of public Cry-toxins, Bt may compete by producing bacteriocins (Abriouel et al., 2011). The exact interaction mechanisms, which are involved in this model system remain to be shown. Importantly, we show that the virulence difference between double and the corresponding single infections is characterized by GxG interactions within the host. These interactions were not revealed by the interaction term between the genotypes of coinfecting strains (see Table 1).
Surprisingly, we did not find similar trends for population size as we revealed for host survival. For example, population size is not affected by the strain combination tested but survival rate is (Table 3). Since we took our measurements on a certain day, we cannot draw any conclusions about lifetime fitness of hosts and about potential costs. However, C. elegans is known to produce the majority of offspring during the first days of adulthood (Huang et al., 2004). One explanation for our results might be that reproduction is not affected by the genotypes of Bt strains, and that hosts produce their offspring before they die. However, since they die earlier when they are infected by certain Bt genotype combinations but population size remains the same, reproduction could also be shifted toward early reproduction. Alternatively, this could be caused by host genotype effects. Since we used a genetically diverse worm population, it is possible that the host genotype influences the interaction between different bacterial genotypes.
Thus, host–parasite interactions can be influenced by genotype specific within-host interactions. Since host–parasite GxG interactions have been revealed for many model systems including the one we used here (Schulte et al., 2010), it is likely that genotype specific within-host interactions are influenced by the host genotype and vice versa, resulting in GxGxG interactions. A potential mechanism for GxGxG interactions could be immune mediated interaction between parasites: depending on the host genotype, the immune-system could be more effective against some parasite genotypes and thus facilitate infections by others (Cox, 2001; Read and Taylor, 2001; Mideo, 2009). The general idea that a third party can influence GxG interactions is not new. For example, host–parasite GxG interactions are known to be influenced or even mediated by the genotypes of host endosymbionts (Oliver et al., 2005; Tetard-Jones et al., 2007; Koch and Schmid-Hempel, 2012; Rouchet and Vorburger, 2012). The situation can even be complicated if the interactions are influences by environmental factors (Seppälä et al., 2012).
That host and parasite genotypes do interact is one of the key assumptions for mathematical models on host–parasite coevolution (Hamilton et al., 1990; Agrawal and Lively, 2002). Thus, changes of allele frequencies and changes of genotypes over time require GxG interactions to occur. Here, we reveal additionally genotype specific within-host interactions of different parasite genotypes. The resulting GxGxG interactions should result in even higher evolutionary rates. If the outcome of host–parasite GxG interactions depends on within-host parasite GxG interactions, the frequency of multiple infections should influence host–parasite coevolution. It has indeed been shown that the frequency of multiple infections has strong implications for host–parasite interactions (e.g., Turner and Chao, 1998, 1999), but also the initial frequency of each strain is of importance (Taylor et al., 1997; Harrison et al., 2006; Rumbaugh et al., 2009; Zwart et al., 2009). Thus, many factors influence host–parasite coevolution under natural condition, making the outcome and epidemiology even more difficult to predict (e.g., de Roode et al., 2005; Alizon and Lion, 2011). This may explain why examples for predicted evolutionary dynamics are rarely found in nature (but see Burdon and Thrall, 2000; Jokela et al., 2003; Ebert, 2008; Dupas et al., 2009; Laine, 2009). Within-host GxG interaction is another factor that should be considered in the study of host–parasite interactions.
Our study provides an intriguing experimental evidence for the potential importance of GxG interactions among distinct pathogen genotypes within their host. We reveal that GxG within-host interactions can occur between genotypes of the same species. Since genotypes may be more or less virulent in combination than the corresponding single infections, GxG within-host interactions are furthermore likely include different interaction mechanisms and to influence the interaction between host and parasites, which is also genotype specific in this model system. What remains to be shown in future experiments is whether the exact outcome of the interaction between different parasite genotypes also depends on the host genotype. Taken together, GxG interactions between coinfecting parasites may have strong implications for epidemiology and virulence evolution.
Author Contributions
Rebecca D. Schulte designed the experiment. Rebecca D. Schulte and Joy Bose acquired, analyzed and interpreted the data and wrote the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Michael Keuler, Michaela Kloesener, and Verena Langer for experimental help, Thomas D’Souza and Hinrich Schulenburg for statistical advice and Michaela Kloesener, Nadine Timmermeyer and three reviewers for helpful comments on the manuscript. Thanks to the members of the Department of Behavioral Biology, University of Osnabrueck and members of the Bacillus-Invertebrate Cluster of the Priority Programme SPP1399 funded by the DFG for discussion of the data. This research was funded by the DFG grant SCHU2522/1-1 of the German Science foundation.
References
Abriouel, H., Franz, C. M. A. P., Ben Omar, N., and Galvez, A. (2011). Diversity and applications of Bacillus bacteriocins. FEMS Microbiol. Rev. 35, 201–232. doi: 10.1111/j.1574-6976.2010.00244.x
Agrawal, A., and Lively, C. M. (2002). Infection genetics: gene-for-gene versus matching-alleles models and all points in between. Evol. Ecol. Res. 4, 79–90.
Alizon, S., De Roode, J. C., and Michalakis, Y. (2013). Multiple infections and the evolution of virulence. Ecol. Lett. 16, 556–567. doi: 10.1111/ele.12076
Alizon, S., and Lion, S. (2011). Within-host parasite cooperation and the evolution of virulence. Proc. Biol. Sci. 278, 3738–3747. doi: 10.1098/rspb.2011.0471
Alper, S., Mcbride, S. J., Lackford, B., Freedman, J. H., and Schwartz, D. A. (2007). Specificity and complexity of the Caenorhabditis elegans innate immune response. Mol. Cell. Biol. 27, 5544–5553. doi: 10.1128/MCB.02070-06
Borgonie, G., Claeys, M., Leyns, F., Arnaut, G., De Waele, D., and Coomans, A. (1996). Effect of nematicidal Bacillus thuringiensis strains on free-living nematodes. 2. Ultrastructural analysis of the intoxication process in Caenorhabditis elegans. Fundam Appl. Nematol. 19, 407–414.
Borgonie, G., Van Driessche, R., Leyns, F., Arnaut, G., De Waele, D., and Coomans, A. (1995). Germination of Bacillus thuringiensis spores in Bacteriophagous Nematodes (Nematoda: Rhabditida). J. Invertebr. Pathol. 65, 61–67. doi: 10.1006/jipa.1995.1008
Brown, S. P., Hochberg, M. E., and Grenfell, B. T. (2002). Does multiple infection select for raised virulence? Trends Microbiol. 10, 401–405. doi: 10.1016/S0966-842X(02)02413-7
Buckling, A., and Brockhurst, M. A. (2008). Kin selection and the evolution of virulence. Heredity 100, 484–488. doi: 10.1038/sj.hdy.6801093
Burdon, J. J., and Thrall, P. H. (2000). Coevolution at multiple spatial scales: Linum marginale-Melampsora lini – from the individual to the species. Evol. Ecol. 14, 261–281. doi: 10.1023/a:1011042721274
Carius, H. J., Little, T. J., and Ebert, D. (2001). Genetic variation in a host–parasite association: potential for coevolution and frequency-dependent selection. Evolution 55, 1136–1145. doi: 10.1111/j.0014-3820.2001.tb00633.x
Cox, F. E. G. (2001). Concomitant infections, parasites and immune responses. Parasitology 122, S23–S38. doi: 10.1017/s003118200001698x
de Roode, J. C., and Altizer, S. (2010). Host–parasite genetic interactions and virulence-transmission relationships in natural populations of Monarch butterflies. Evolution 64, 502–514. doi: 10.1111/j.1558-5646.2009.00845.x
de Roode, J. C., Pansini, R., Cheesman, S. J., Helinski, M. E. H., Huijben, S., Wargo, A. R., et al. (2005). Virulence and competitive ability in genetically diverse malaria infections. Proc. Natl. Acad. Sci. U.S.A. 102, 7624–7628. doi: 10.1073/pnas.0500078102
Dupas, S., Dubuffet, A., Carton, Y., and Poirie, M. (2009). “Local, geographic and phylogenetic scales of coevolution in Drosophila-parasitoid interactions,” in Advances in Parasitology, ed. G. Prevost (San Diego, USA: Elsevier Academic Press Inc.), 281–295.
Ebert, D. (2008). Host–parasite coevolution: insights from the Daphnia-parasite model system. Curr. Opin. Microbiol. 11, 290–301. doi: 10.1016/j.mib.2008.05.012
George, Z., and Crickmore, N. (2012). “Bacillus thuringiensis applications in agriculture,” in Bacillus thuringiensis Biotechnology, ed. E. Sansinenea (Berlin: Springer Science+Business Media B.V.), 19–39.
Gorton, A. J., Heath, K. D., Pilet-Nayel, M. -L., Baranger, A., and Stinchcombe, J. R. (2012). Mapping the Genetic Basis of Symbiotic Variation in Legume-Rhizobium Interactions in Medicago truncatula. G3 (Bethesda) 2, 1291–1303. doi: 10.1534/g3.112.003269
Grech, K., Watt, K., and Read, A. F. (2006). Host–parasite interactions for virulence and resistance in a malaria model system. J. Evol. Biol. 19, 1620–1630. doi: 10.1111/j.1420-9101.2006.01116.x
Griffitts, J. S., and Aroian, R. V. (2005). Many roads to resistance: how invertebrates adapt to Bt toxins. Bioessays 27, 614–624. doi: 10.1002/bies.20239
Gsell, A. S., Domis, L. N. D. S., Van Donk, E., and Ibelings, B. W. (2013). Temperature alters host genotype-specific susceptibility to chytrid infection. PLoS ONE 8:e71737. doi: 10.1371/journal.pone.0071737
Hamilton, W. D., Axelrod, R., and Tanese, R. (1990). Sexual reproduction as an adaptation to resist parasites (A review). Proc. Natl. Acad. Sci. U.S.A. 87, 3566–3573. doi: 10.1073/pnas.87.9.3566
Harrison, F., Browning, L. E., Vos, M., and Buckling, A. (2006). Cooperation and virulence in acute Pseudomonas aeruginosa infections. BMC Biol. 4:21. doi: 10.1186/1741-7007-4-21
Heath, K. D. (2010). Intergenomic epistasis and coevolutionary constraints in plants and rhizobia. Evolution 64, 1446–1458. doi: 10.1111/j.1558-5646.2009.00913.x
Hoeksema, J. D., and Thompson, J. N. (2007). Geographic structure in a widespread plant-mycorrhizal interaction: pines and false truffles. J. Evol. Biol. 20, 1148–1163. doi: 10.1111/j.1420-9101.2006.01287.x
Huang, C., Xiong, C. J., and Kornfeld, K. (2004). Measurements of age-related changes of physiological processes that predict lifespan of Caenorhabditis elegans. Proc. Natl. Acad. Sci. U.S.A. 101, 8084–8089. doi: 10.1073/pnas.0400848101
Jokela, J., Lively, C. M., Dybdahl, M. F., and Fox, J. A. (2003). Genetic variation in sexual and clonal lineages of a freshwater snail. Biol. J. Linn. Soc. 79, 165–181. doi: 10.1046/j.1095-8312.2003.00181.x
Koch, H., and Schmid-Hempel, P. (2012). Gut microbiota instead of host genotype drive the specificity in the interaction of a natural host–parasite system. Ecol. Lett. 15, 1095–1103. doi: 10.1111/j.1461-0248.2012.01831.x
Laine, A. L. (2009). Role of coevolution in generating biological diversity: spatially divergent selection trajectories. J. Exp. Bot. 60, 2957–2970. doi: 10.1093/jxb/erp168
Lambrechts, L., Chevillon, C., Albright, R. G., Thaisomboonsuk, B., Richardson, J. H., Jarman, R. G., et al. (2009). Genetic specificity and potential for local adaptation between dengue viruses and mosquito vectors. BMC Evol. Biol. 9:160. doi: 10.1186/1471-2148-9-160
Lambrechts, L., Fellous, S., and Koella, J. C. (2006). Coevolutionary interactions between host and parasite genotypes. Trends Parasitol. 22, 12–16. doi: 10.1016/j.pt.2005.11.008
Lazzaro, B. P., and Little, T. J. (2009). Immunity in a variable world. Philos. Trans. R. Soc. Lond. B Biol. Sci. 364, 15–26. doi: 10.1098/rstb.2008.0141
Lereclus, D., Arantes, O., Chaufaux, J., and Lecadet, M. M. (1989). Transformation and expression of a clones delta-endotoxin gene in Bacillus thuringiensis. FEMS Microbiol. Lett. 60, 211–218. doi: 10.1111/j.1574-6968.1989.tb03448.x
Luijckx, P., Ben-Ami, F., Mouton, L., Du Pasquier, L., and Ebert, D. (2011). Cloning of the unculturable parasite Pasteuria ramosa and its Daphnia host reveals extreme genotype-genotype interactions. Ecol. Lett. 14, 125–131. doi: 10.1111/j.1461-0248.2010.01561.x
Mboup, M., Bahri, B., Leconte, M., De Vallavieille-Pope, C., Kaltz, O., and Enjalbert, J. (2013). Genetic structure and local adaptation of European wheat yellow rust populations: the role of temperature-specific adaptation. Evol. Appl. 5, 341–352. doi: 10.1111/j.1752-4571.2011.00228.x
Mideo, N. (2009). Parasite adaptations to within-host competition. Trends Parasitol. 25, 261–268. doi: 10.1016/j.pt.2009.03.001
Oliver, K. M., Moran, N. A., and Hunter, M. S. (2005). Variation in resistance to parasitism in aphids is due to symbionts not host genotype. Proc. Natl. Acad. Sci. U.S.A. 102, 12795–12800. doi: 10.1073/pnas.08.007
Payne, J. M. (1992). Isolates of Bacillus thuringiensis that are active against nematodes. US patent 5151363, US patent 5424410.
Payne, J. M., Cannon, R. J. C., and Bagley, A. L. (1993). Bacillus thuringiensis isolates for controlling acarides. US Patent 5424410 A.
Rauch, G., Kalbe, M., and Reusch, T. B. H. (2006). One day is enough: rapid and specific host–parasite interactions in a stickleback-trematode system. Biol. Lett. 2, 382–384. doi: 10.1098/rsbl.2006.0462
Raymond, B., Davis, D., and Bonsall, M. B. (2007). Competition and reproduction in mixed infections of pathogenic and non-pathogenic Bacillus spp. J. Invertebr. Pathol. 96, 151–155. doi: 10.1016/j.jip.2007.03.003
Read, A. F., and Taylor, L. H. (2001). The ecology of genetically diverse infections. Science 292, 1099–1102. doi: 10.1126/science.1059410
Rockman, M. V., and Kruglyak, L. (2009). Recombinational landscape and population genomics of Caenorhabditis elegans. PLoS Genet. 5:3. doi: 10.1371/journal.pgen.1000419
Rouchet, R., and Vorburger, C. (2012). Strong specificity in the interaction between parasitoids and symbiont-protected hosts. J. Evol. Biol. 25, 2369–2375. doi: 10.1111/j.1420-9101.2012.02608.x
Rumbaugh, K. P., Diggle, S. P., Watters, C. M., Ross-Gillespie, A., Griffin, A. S., and West, S. A. (2009). Quorum sensing and the social evolution of bacterial virulence. Curr. Biol. 19, 341–345. doi: 10.1016/j.cub.2009.01.050
Sadd, B. M. (2011). Food-environment mediates the outcome of specific interactions between a bumblebee and its trypanosome parasite. Evolution 65, 2995–3001. doi: 10.1111/j.1558-5646.2011.01345.x
Salvaudon, L., Heraudet, V., and Shykoff, J. A. (2007). Genotype-specific interactions and the trade-off between host and parasite fitness. BMC Evol. Biol. 7:189. doi: 10.1186/1471-2148-7-189
Schmid-Hempel, P., and Ebert, D. (2003). On the evolutionary ecology of specific immune defence. Trends Ecol. Evol. 18, 27–32. doi: 10.1016/S0169-5347(02)00013-7
Schnepf, E., Crickmore, N., Van Rie, J., Lereclus, D., Baum, J., Feitelson, J., et al. (1998). Bacillus thuringiensis and its pesticidal crystal proteins. Microbiol. Mol. Biol. Rev. 62, 775–806.
Schnepf, H. E., Schwab, G. E., Payne, J. M., Narva, K. E., and Foncerrada, L. (2001). Nematicidal proteins. US patent 20010010932.
Schulenburg, H., and Ewbank, J. J. (2004). Diversity and specificity in the interaction between Caenorhabditis elegans and the pathogen Serratia marcescens. BMC Evol. Biol. 4:49. doi: 10.1186/1471-2148-4-49
Schulenburg, H., Hoeppner, M. P., Weiner, J., and Bornberg-Bauer, E. (2008). Specificity of the innate immune system and diversity of C-type lectin domain (CTLD) proteins in the nematode Caenorhabditis elegans. Immunobiology 213, 237–250. doi: 10.1016/j.imbio.2007.12.004
Schulenburg, H., and Müller, S. (2004). Natural variation in the response of Caenorhabditis elegans towards Bacillus thuringiensis. Parasitology 128, 433–443. doi: 10.1017/S003118200300461X
Schulte, R. D., Makus, C., Hasert, B., Michiels, N. K., and Schulenburg, H. (2010). Multiple reciprocal adaptations and rapid genetic change upon experimental coevolution of an animal host and its microbial parasite. Proc. Natl. Acad. Sci. U.S.A. 107, 7359–7364. doi: 10.1073/pnas.1003113107
Seppälä, O., Karvonen, A., Rellstab, C., Louhi, K. -R., and Jokela, J. (2012). Reciprocal interaction matrix reveals complex genetic and dose-dependent specificity among coinfecting parasites. Am. Nat. 180, 306–315. doi: 10.1086/666985
Sicard, M., Hering, S., Schulte, R., Gaudriault, S., and Schulenburg, H. (2007). The effect of Photorhabdus luminescens (Enterobacteriaceae) on the survival, development, reproduction and behaviour of Caenorhabditis elegans (Nematoda: Rhabditidae). Environ. Microbiol. 9, 12–25. doi: 10.1111/j.1462-2920.2006.01099.x
Stiernagle, T. (2006). “Maintenance of C. elegans,” in The C. elegans Research Community (http://www.wormbook.org), ed. WormBook. doi: 10.1895/wormbook.1.101.1
Taylor, L. H., Walliker, D., and Read, A. F. (1997). Mixed-genotype infections of the rodent malaria Plasmodium chabaudi are more infectious to mosquitoes than single infections. Parasitology 115, 121–132. doi: 10.1017/S0031182097001145
Teotónio, H., Carvalho, S., Manoel, D., Roque, M., and Chelo, I. M. (2012). Evolution of outcrossing in experimental populations of Caenorhabditis elegans. PLoS ONE 7:e35811. doi: 10.1371/journal.pone.0035811
Tetard-Jones, C., Kertesz, M. A., Gallois, P., and Preziosi, R. F. (2007). Genotype-by-genotype interactions modified by a third species in a plant-insect system. Am. Nat. 170, 492–499. doi: 10.1086/520115
Thomas, M. B., and Blanford, S. (2003). Thermal biology in insect-parasite interactions. Trends Ecol. Evol. 18, 344–350. doi: 10.1016/s0169-5347(03)00069-7
Thompson, J. N. (2005). The Geographic Mosaic of Coevolution. Chicago, IL: University of Chicago Press.
Tollenaere, C., and Laine, A. L. (2012). Investigating the production of sexual resting structures in a plant pathogen reveals unexpected self-fertility and genotype-by-environment effects. J. Evol. Biol. 26, 1716–1726. doi: 10.1111/jeb.12169
Turner, P. E., and Chao, L. (1998). Sex and the evolution of intrahost competition in RNA virus phi 6. Genetics 150, 523–532.
Turner, P. E., and Chao, L. (1999). Prisoner’s dilemma in an RNA virus. Nature 398, 441–443. doi: 10.1038/18913
Webster, J. P., and Woolhouse, M. E. J. (1998). Selection and strain specificity of compatibility between snail intermediate hosts and their parasitic schistosomes. Evolution 52, 1627–1634. doi: 10.2307/2411336
Wei, J. Z., Hale, K., Carta, L., Platzer, E., Wong, C., Fang, S. C., et al. (2003). Bacillus thuringiensis crystal proteins that target nematodes. Proc. Natl. Acad. Sci. U.S.A. 100, 2760–2765. doi: 10.1073/pnas.0538072100
Wolinska, J., and King, K. C. (2009). Environment can alter selection in host–parasite interactions. Trends Parasitol. 25, 236–244. doi: 10.1016/j.pt.2009.02.004
Wong, D., Bazopoulou, D., Pujol, N., Tavernarakis, N., and Ewbank, J. J. (2007). Genome-wide investigation reveals pathogen-specific and shared signatures in the response of Caenorhabditis elegans to infection. Genome Biol. 8, R194. doi: 10.1186/gb-2007-8-9-r194
Keywords: GxG, parasite, host, multiple infection, within-host interaction, Caenorhabditis elegans, Bacillus thuringiensis, virulence
Citation: Bose J and Schulte RD (2014) Testing GxG interactions between coinfecting microbial parasite genotypes within hosts. Front. Genet. 5:124. doi: 10.3389/fgene.2014.00124
Received: 05 December 2013; Accepted: 21 April 2014;
Published online: 14 May 2014.
Edited by:
Louis Lambrechts, Institut Pasteur, FranceReviewed by:
Suleyman Yildirim, Walter Reed Army Institute of Research, USASimon Fellous, Institut National de la Recherche Agronomique, France
Copyright © 2014 Bose and Schulte. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rebecca D. Schulte, Department of Behavioral Biology, University of Osnabrueck, Barbarastraße 11, 49076 Osnabrueck, Germany e-mail: rebecca.schulte@biologie.uni-osnabrueck.de