Evolution of Caenorhabditis elegans host defense under selection by the bacterial parasite Serratia marcescens

Abstract

Parasites can impose strong selection on hosts. In response, some host populations have adapted via the evolution of defenses that prevent or impede infection by parasites. However, host populations have also evolved life history shifts that maximize host fitness despite infection. Outcrossing and self-fertilization can have contrasting effects on evolutionary trajectories of host populations. While selfing and outcrossing are known to affect the rate at which host populations adapt in response to parasites, these mating systems may also influence the specific traits that underlie adaptation to parasites. Here, we determined the role of evolved host defense versus altered life history,in mixed mating (selfing and outcrossing) and obligately outcrossing C. elegans host populations after experimental evolution with the bacterial parasite, S. marcescens. Similar to previous studies, we found that both mixed mating and obligately outcrossing host populations adapted to S. marcescens exposure, and that the obligately outcrossing populations exhibited the greatest rates of adaptation. Regardless of the host population mating system, exposure to parasites did not significantly alter reproductive timing or total fecundity over the course of experimental evolution. However, both mixed mating and obligately outcrossing host populations exhibited significantly reduced mortality rates in the presence of the parasite after experimental evolution. Therefore, adaptation in both the mixed mating and obligately outcrossing populations was driven, at least in part, by the evolution of increased host defense and not changes in host life history. Thus, the host mating system altered the rate of adaptation, but not the nature of adaptive change in the host populations.

About the Authors:

McKenna J. Penley

Affiliation: Department of Biology, Emory University, Atlanta, Georgia, United States of America

Giang T. Ha

Affiliation: Department of Biology, Emory University, Atlanta, Georgia, United States of America

Levi T. Morran

* E-mail: [email protected]

Affiliation: Department of Biology, Emory University, Atlanta, Georgia, United States of America

Introduction

Parasites are prevalent in nature and capable of imposing significant fitness costs on their hosts. Consequently, host populations have evolved numerous strategies for mitigating or preventing infection by parasites [1, 2]. Hosts commonly adapt via the evolution of novel or enhanced forms of defense in response to parasites. The specific traits and mechanisms that underlie host defense can vary greatly, from heighted cell-mediated immunity to recognition and avoidance of the parasite upon exposure [3]. As opposed to an inducible defense strategy, some host populations have instead evolved altered life histories [4–6]. Indeed, after evolving under frequent threat of infection, some host populations exhibit reduced times to sexual maturity and earlier reproduction as a means of maximizing fitness [7, 8]. Thus, the nature of host adaptation to parasites can vary greatly, as hosts employ many different mechanisms of combating fitness loss due to parasite infection.

Many factors may influence the evolutionary trajectories of host populations in response to selection imposed by parasites [9]. The host mating system can determine whether host populations respond sufficiently to selection from coevolving parasites or are driven to extinction [10]. Further, self-fertilization and outcrossing, different forms of sexual reproduction, can have very different effects on the evolutionary trajectories of populations in response to natural selection [11–19]. In general, outcrossing can facilitate more rapid adaptation to selection than self-fertilization by accelerating the breakdown of linkage disequilibrium within populations [20–22]. Further, breaking genetic linkage can permit the assembly of novel beneficial allelic combinations [23, 24] and disassociate beneficial mutations from deleterious mutations [21, 22, 25–27]. Additionally, the propensity of selfing to generate and maintain linkage disequilibrium can slow adaptation due to selective interference from deleterious mutations at loci in close proximity to beneficial alleles. Overall, adaptation via outcrossing may result in the fixation of fewer deleterious alleles and incorporate a greater number of loci relative to self-fertilization. By increasing the efficacy of selection, outcrossing may also permit selection to act on a greater number of traits simultaneously, whereas adaptation in selfing populations may often be restricted to the trait with the greatest contribution to fitness [25–28]. Therefore, the host mating system may play a prominent role in determining the mechanisms underlying host adaptation to parasites. Nonetheless, it is generally unclear how selfing and outcrossing may affect the mechanisms and specific traits that contribute to host adaptation in response to selection imposed by parasites

Caenorhabditis elegans are androdioecious nematodes that reproduce via self-fertilization and outcrossing (mixed mating) [29]. However, C. elegans populations harboring the fog-2(q7) mutant allele, which disables sperm production in hermaphrodites, are obligately outcrossing [30]. Conversely, C. elegans populations carrying the xol-1(y9) mutation, which is lethal in males, are obligately selfing [31]. Serratia marcescens is a virulent bacterial parasite that can kill C. elegans hosts when consumed [32]. Previous host-parasite experimental evolution studies, in which hosts were consistently exposed to the same genotype of S. marcescens each generation, found that both obligately outcrossing and mixed mating C. elegans populations rapidly adapted to S. marcescens. However, obligately selfing populations did not adapt to the parasites [10, 20, 33–35], and the obligately outcrossing populations exhibited greater increases in fitness relative to the mixed mating populations [20, 34]. Therefore, high levels of outcrossing facilitated greater rates of adaptation to the parasites. Indeed, greater rates of outcrossing evolved in mixed mating host populations, indicating that outcrossing was favored by selection imposed by the parasite [10, 20, 33].

Apart from different rates of adaptive change, it is unclear if mixed mating can alter the evolutionary trajectories of host populations relative to obligate outcrossing. In a previous experiment, hosts were required to survive exposure to S. marcescens strain Sm2170, locate and consume their Escherichia coli food source, then reproduce before experimental passage [10]. Within the context of the experimental system, host fitness was determined by both survival in the presence of the parasite and fecundity after parasite exposure. The host populations may have adapted to the parasite and the selection regime through increased host defense and/or the evolution of greater or earlier reproduction. Here, we tested for adaptation in these mixed mating and obligately outcrossing host populations that were passaged either in the presence of live S. marcescens Sm2170 or heat-killed S. marcescens, as a control, for thirty generations [10]. We then investigated the evolutionary trajectories of the host populations by linking changes in fitness with changes in host defense versus shifts in reproductive timing.

Materials and methods

Experimental host populations

Experimental host populations consisted of 5 replicates of mixed mating C. elegans strain PX382 and five replicates of obligately outcrossing C. elegans strain PX386, both from an overall CB4856 background [20], which was originally obtained from the Caenorhabditis Genetics Center. Strain PX386 harbors the fog-2(q71) mutation which prevents hermaphrodites from producing sperm, therefore creating an obligately outcrossing population. All experimental populations were independently mutagenized with Ethyl Methanesulfonate (EMS) to introduce novel genetic variation prior to selection [10]. The ancestor to each experimental population was frozen and stored at -80°C prior to experimental evolution, while the experimental populations were frozen after 30 generations of experimental evolution. Each of the ten replicate populations (five mixed mating and five obligately outcrossing) was then replicated across each treatment, such that a sample of each replicate population underwent each treatment.

Experimental evolution

All populations underwent 30 generations of passage on Serratia Selection Plates (SSP) as described in in a previous study [10]. Briefly, S. marcescens strain Sm2170, obtained from the lab of Curtis Lively at Indiana University, was grown on approximately 1/3 of a Petri dish, OP50 E. coli on the opposite 1/3 of the plate, with ampicillin (200 μg/mL) spread in the middle 1/3 of the plate. C. elegans are plated on the S. marcescens lawn and must survive exposure to the parasite and reproduce on the OP50 lawn to have measurable fitness during the experiment. Five genetically independent mixed mating and five genetically independent obligately outcrossing populations were exposed to 2 different treatments, (Evolution) exposure to a fixed strain of S. marcescens Sm2170, and (Control) exposure to heat-killed Sm2170. SSPs were constructed using NGM Lite (US Biological, Swampscott, MA) and stored at 20°C throughout the experiment. C. elegans host populations were transferred every generation following methods outlined in the previous study [10].

Change in fitness

Change in fitness after 30 generations of experimental evolution was measured using Competitive Fitness Assays (CFA) [20]. CFAs were conducted for the ancestral populations and experimental populations, all of which were thawed after storage at -80°C, using methods from [10]. Briefly, 100 experimental or ancestral individuals and 100 GFP labeled tester strain individuals (strain JK2735) were placed on SSPs with live S. marcescens. Individuals were permitted eat, grow, and reproduce. Then, we measured the ratio of experimental or ancestral population offspring to tester strain offspring after one generation [20]. At least two replicate CFAs were conducted for each host population. Importantly, the GFP marker is dominant, therefore the offspring of outcrossing between experimental males and GFP-labeled hermaphrodites (all tester strain individuals were hermaphrodites) appear as tester strain offspring. Given that the obligately outcrossing host populations harbor significantly greater numbers of males, relative to mixed mating populations, the competitive fitness measurements of the obligately outcrossing populations are more conservative estimates. We performed a Shapiro-Wilk test and Levene’s test, and found that the data did not violate assumptions of normality nor homogeneity of variance. Then we performed an ANOVA in JMP 12 (SAS Institute, Cary, NC) testing the effects of host mating system (mixed mating versus obligately outcrossing), treatment (Ancestor, Control, Evolution), host replicate population (treated as a random factor nested within mating system), and the interaction of mating system and treatment on the mean ratio of experimental to tester strain offspring. Further, we conducted least squares mean linear contrast tests within the mixed mating and obligately outcrossing populations contrasting the means of Evolution treatment with the means of Control and Ancestral populations. Lastly, we performed a Kruskal Wallis test to analyze the effects of host mating system on the change in fitness relative to the ancestral population. The change in host fitness relative to the ancestral population was calculated as described in [20]. Bonferroni correction was applied to this data set due to multiple statistical tests (P < 0.025).

Host defense

Host defense was measured using Survival Assays (SA), wherein 200 C. elegans from the either an ancestral or experimental (after 30 generations of selection) population were exposed to live Sm2170 on SSPs and scored for survival after 48 hours. We calculated the mean mortality rates for each experimental and ancestral population, running at least 3 replicate assays for each population. We chose a 48 hour time period because Sm2170 rapidly kills infected C. elegans, often within 24 hours of consumption. C. elegans hosts move off the Sm2170 lawn and onto OP50 within the first 24 hours of exposure to the SSP. Further, C. elegans hosts must reproduce within the first 48 hours of exposure on the SSP to ensure that their offspring are passaged to the next generation. So, measuring rates of mortality after 48 hours of exposure is an indicator of the number of individuals that survived to reproductive age in the context of our experiment. We performed Shapiro-Wilk and Levene’s tests, and found that the data did not violate assumptions of ANOVA. We then performed an ANOVA in JMP 12 testing the effects of host mating system, treatment, host replicate population (treated as a random factor nested within mating system), and the interaction of mating system and treatment on the mean mortality rates calculated for each ancestral and experimental population. Further, we conducted least squares mean linear contrast tests within treatment to test the effects of individual treatments (Evolution versus Control and Evolution versus Ancestral populations).

Host fecundity

Host fecundity was calculated for both the ancestral and experimental populations (after 30 generations of selection) by measuring reproduction over seven days, in the absence of the parasite, beginning at the L4 lifestage of the hermaphrodite or female. For obligately outcrossing populations, L4 female offspring were picked onto 35mm plates seeded with 20μl of OP50, along with two adult males from the female’s respective replicate population. All three worms were passaged together throughout the assay to permit further reproduction. For mixed mating populations, one hermaphrodite per plate, starting at the L4 lifestage, was passaged throughout the assay. Importantly, fecundity was assessed using the predominant mating system employed by the populations after 30 generations of experimental evolution. Specifically, the obligately outcrossing populations reproduced via obligate outcrossing, whereas most offspring (~80%) in the mixed mating populations were produced via selfing [10]. C. elegans were transferred to a new 35mm plate at 48 hours after the initial plating. Then focal individuals were again transferred after a subsequent 36 hours and 120 hours from the beginning of the assay to clearly delineate between focal parents and offspring, and to assess reproductive timing. We conducted this assay for five replicate females or hermaphrodites from each replicate population. The quantity of L1 offspring produced was measured after each passage and on day seven, and the mean lifetime fecundity and 48-hour fecundity was calculated for each replicate population. For both the mean lifetime and mean 48-hour fecundity data sets, we found that the data did not violate the assumptions of ANOVA using Shapiro-Wilk and Levene’s tests. Then, we performed an ANOVA in JMP 12 testing the effects of host mating system, treatment, host replicate population (treated as a random factor nested within mating system), and the interaction of mating system and treatment on the mean fecundity calculated for each ancestral and experimental population. We also conducted least squares mean linear contrast tests within treatment to test the effects of individual treatments (Control versus Evolution and Control versus Ancestral populations). Bonferroni correction was applied to this data due to multiple statistical tests (P < 0.025).

Results

Change in fitness

After thirty generations of exposure to the same strain of S. marcescens in the Evolution treatment, both the mixed mating and obligately outcrossing populations exhibited increased competitive fitness in the presence of the parasite, relative to the Control and Ancestral populations (Fig 1A and Table 1; Mixed mating: Evolution vs Control; F_1,16 = 8.25, P = 0. 011, Evolution vs Ancestor; F_1,16 = 4.70, P = 0. 045) (Fig 1B and Table 1; Obligate outcrossing: Evolution vs Control; F_1,16 = 18.14, P = 0.0006, Evolution vs Ancestor; F_1,16 = 59.64, P < 0. 0001). Therefore, both the mixed mating and obligately outcrossing populations adapted to S. marcescens exposure. However, the rate of adaptation differed between mating systems. The obligately outcrossing populations adapted to S. marcescens at a greater rate than the mixed mating populations. Specifically, the obligately outcrossing host populations in the Evolution treatment exhibited a significantly greater increase in fitness relative to their ancestral populations, than did the mixed mating hosts populations (Obligate outcrossing increase in fitness relative to ancestor = 131%, mixed mating increase in fitness relative to ancestor = 24%; χ²₁ = 5.771, P = 0.0163). However, it should be noted that the competitive fitness of the obligately outcrossing ancestral populations was lower than the fitness of the ancestral mixed mating populations (Fig 1), which contributed to the difference in the rates of adaptation between mating systems.

[Figure omitted. See PDF.]

Fig 1. Competitive fitness during exposure to S. marcescens.

Populations exposed to the Control and Evolution treatments, as well their Ancestral populations, were competed against a tester strain in the presence of S. marcescens. Experimental and tester strain individuals were exposed to the parasite, starting at a 1:1 ratio, for one generation. Greater mean proportions of experimental individuals indicate greater competitive fitness. (A) Mixed mating populations exposed to the Evolution treatment exhibit greater competitive fitness than the Ancestral or Control populations. (B) Obligately outcrossing populations exposed to the Evolution treatment exhibit greater competitive fitness than their respective Ancestral and Control populations.

https://doi.org/10.1371/journal.pone.0181913.g001

[Figure omitted. See PDF.]

Table 1. ANOVA table for competitive fitness.

https://doi.org/10.1371/journal.pone.0181913.t001

Host defense

The increased fitness in the presence of S. marcescens exhibited by C. elegans populations in the Evolution treatment could be the result of increased levels of defense against the parasite, a life history shift resulting in increased host fecundity, or a combination of both changes. We assessed changes in host defense by exposing Control, Evolution, and Ancestral populations to S. marcescens on SSPs and we measured rates of host mortality over 48 hours. Overall, we found a significant treatment effect (Table 2), in the mixed mating and obligately outcrossing populations from the Evolution treatment exhibiting significantly reduced mortality rates relative to Control (Fig 2 and Table 2; F_1,16 = 20.69 P = 0. 0003) and Ancestral populations (Fig 2 and Table 2; F_1,16 = 17.20, P = 0.0008). Therefore, adaptation in both the mixed mating and obligately outcrossing populations was at least partially driven by increased host defense.

[Figure omitted. See PDF.]

Fig 2. Host defense during exposure to S. marcescens.

Populations evolved under the Control and Evolution treatments, as well as Ancestral populations, were exposed to S. marcescens on SSPs for 48 hours and then host mortality rates were assessed. (A) Mixed mating populations from the Evolution treatment exhibited reduced mortality rates. (B) Obligately outcrossing populations from the Evolution treatment exhibited reduced mortality rates. Overall, persistent exposure to the parasite selected for the evolution of increased host defense regardless of the host mating system.

https://doi.org/10.1371/journal.pone.0181913.g002

[Figure omitted. See PDF.]

Table 2. ANOVA table for mean host mortality rates.

https://doi.org/10.1371/journal.pone.0181913.t002

Host fecundity

We assessed changes in host fecundity over the course of experimental evolution as a potential mechanism contributing to adaptation. Specifically, we measured the lifetime fecundity, in the absence of the parasite, for populations that evolved under the Control and Evolution treatments, as well as the Ancestral populations. Overall, we did not observe a treatment effect in the mixed mating populations as the Control and Evolution populations did not exhibit significantly different levels of fecundity (Fig 3A and Table 3; F_1,14 = 0.077, P = 0.785) However, we found that the mean lifetime fecundity increased in the obligately outcrossing Control populations relative to the Evolution (Fig 3B and Table 3; F_1,14 = 18.49, P = 0.0007) and Ancestral (Fig 3B and Table 3; F_1,14 = 26.11, P = 0.0002) populations. Nonetheless, the Control populations did not adapt to S. marcescens exposure (Fig 1), whereas the populations exposed to the Evolution treatment exhibited adaptive change. Therefore the adaptive change in the Evolution populations cannot be attributed to increases in host fecundity in the absence of the parasite.

[Figure omitted. See PDF.]

Fig 3. Host lifetime fecundity.

Lifetime fecundity was assessed in the absence of S. marcescens for Control, Evolution, and Ancestral populations in the (A) mixed mating populations and the (B) obligately outcrossing populations. Control populations exhibited greater levels of fecundity than the Ancestral and Evolution populations in obligately outcrossing populations, but no treatment differences were detected in the mixed mating populations.

https://doi.org/10.1371/journal.pone.0181913.g003

[Figure omitted. See PDF.]

Table 3. ANOVA table for mean lifetime fecundity.

https://doi.org/10.1371/journal.pone.0181913.t003

Within the context of our experiment, fecundity within the first 48 hours of exposure to the SSP is likely a better measure of fitness than the host’s total lifetime fecundity. During experimental evolution, larval hosts were transferred to the SSP and then their offspring were transferred to a new SSP after 4 days. Early reproduction likely disproportionately contributed to fitness relative to lifetime fecundity, given the short time window for reproduction before subsequent transfer. Further, S. marcescens strain Sm2170 can rapidly kill infected hosts, often in less than 48 hours (10). Therefore, rapid reproduction may have been under selection in our experiment to permit reproduction despite infection. We measured the mean 48-hour fecundity of the Control, Evolution, and Ancestral populations in the absence of the parasite. Fecundity was measured in the absence of the parasite to determine the role of shifts in reproductive timing as a potential mechanism underlying changes in fitness exhibited by the experimental populations. Overall, we found significantly greater levels of mean 48-hour fecundity in the Control populations (Fig 4 and Table 4). Additionally, we found that the mixed mating populations exhibited greater 48-hour fecundity than the obligately outcrossing populations (Fig 4 and Table 4). Importantly, these early fecundity differences may be due to mate limitation as a result of reliance on male mating within the obligately outcrossing populations, whereas the mixed mating hermaphrodites self-fertilized. The observed shift in reproductive timing within the Control populations may be the result of adaptation to the experimental evolution regime. Nonetheless, whether lifetime or 48-hour fecundity, the result of increased fecundity in the Control populations does not sufficiently explain the adaptive change observed in the Evolution populations relative to the Controls.

[Figure omitted. See PDF.]

Fig 4. Host 48-hour fecundity.

C. elegans fecundity was measured over a 48 hour period in the absence of S. marcescens. Control populations exhibited increased mean 48-hour fecundity in both the (A) mixed mating and (B) obligately outcrossing populations. Additionally, the mixed mating populations exhibited greater mean 48-hour fecundity than the obligately outcrossing populations.

https://doi.org/10.1371/journal.pone.0181913.g004

[Figure omitted. See PDF.]

Table 4. ANOVA table for mean 48-hour fecundity.

https://doi.org/10.1371/journal.pone.0181913.t004

Discussion

Overall, our results reinforce the results presented by Morran et al [20]. Obligately outcrossing and mixed mating populations evolved in the presence of live S. marcescens both adapted to parasite exposure over 30 generations (Fig 1). Further, the obligately outcrossing populations exhibited a greater increase in fitness relative to their ancestral populations. Nonetheless, the rate and magnitude of adaptation are only two aspects of a population’s evolutionary trajectory, and outcrossing and selfing have the potential to influence many aspects of adaptive change in populations. Given the contrasting genetic consequences of outcrossing and selfing, both the traits evolving under selection and the genetic architecture underlying adaptation may differ between mixed mating and obligately outcrossing populations.

Here, we assessed the role of increased host defense against parasites versus shifts in life history as means of characterizing adaptive evolution in the obligately outcrossing and mixed mating populations. We assayed survival in the presence of S. marcescens as a means of testing host defense, while we measured fecundity in the absence of the parasite to test for shifts in life history. Importantly, we conducted fecundity assays in the absence of the parasite, because parasite-induced fecundity is a form of host defense and we sought to assess general changes in host life history. Regardless of the host mating system, we found that host populations passaged in the presence of the live parasite evolved increased defense to the parasite, resulting in decreased mortality in the presence of the parasite in comparison to ancestral populations and populations passaged with heat-killed parasite (Fig 2). With regards to changes in fecundity, the populations passaged with heat-killed parasites evolved increased 48-hour fecundity (Fig 4). However, evolution under selection imposed by live parasites did not result in increased host fecundity or shifts in host reproductive timing, relative to the Control treatment, for either mixed mating or obligately outcrossing populations (Figs 3 and 4). Therefore, adaptive evolution was, at least in part, due to the evolution of greater defense to the parasite in both the obligately outcrossing and mixed mating populations and not the result of changes in fecundity. Thus, the mixed mating and obligately outcrossing populations followed similar evolutionary trajectories in that exposure to parasite resulted in increased host defense as opposed to general shifts in life history. Yet, host defense is a broad category and there may be multiple mechanisms driving the evolution of defense within the obligately outcrossing and mixed mating populations.

S. marcescens Sm2170 is a highly virulent parasite that can rapidly kill up to 80% of C. elegans hosts upon ingestion in naive host strains [32]. Therefore, it is most likely that selection for defense was stronger than selection for changes in fecundity in host populations. A shift in reproductive timing likely did not yield a sufficient advantage to supersede the evolution of defense in host populations exposed to such a virulent parasite. However, Control obligately outcrossing populations evolved greater 48-hour and total fecundity (Figs 3B and 4B) over the course of the experiment, while the mixed mating populations evolved greater 48-hour fecundity in the Control treatment (Fig 4A). Yet these shifts in fecundity did not result in increased fitness in the presence of the parasite for mixed mating populations and only moderately increased fitness for obligately outcrossing populations (Fig 1). Therefore, it seems that life history changes were favored by the experimental passage protocol, but not in the presence of the parasite. Rather, selection for defense may have superseded selection for changes in fecundity and reproductive timing in the presence of the parasite, regardless of the host mating system.

Different mating systems can have vastly different effects on the evolutionary trajectories of populations [36, 37]. Here, we found that obligate outcrossing increased the rate of adaptation under selection from parasites, relative to rates of adaptation in mixed mating host populations. Greater rates of adaptation in the obligately outcrossing populations may have been driven by a more rapid breakdown of linkage between beneficial and deleterious mutations [36, 37], or by more rapidly uniting beneficial mutations from different lineages [23, 24], relative to the mixed mating populations. It is possible that sexual selection was stronger on males in the obligately outcrossing populations than in the mixed mating populations, which may have also contributed to the different rates of adaptation [38–40]. However, outcrossing rates and male frequencies were maintained at high levels for several generations in the mixed mating populations exposed to parasites [10], so sexual selection is unlikely to account for the differences in the obligately outcrossing and mixed mating populations in the Evolution treatment. Conversely, male frequencies were maintained at much greater levels in the Control obligately outcrossing populations compared to the Control mixed mating populations. Therefore, sexual selection may have contributed to the increased total fecundity exhibited by the obligately outcrossing populations evolved in the absence of parasites (Fig 3B).

Overall, the host mating system did not alter the general response of the host populations. It may be that selection imposed by the parasite favoring defense was strong enough to overcome any constraints specific to mixed mating or obligate outcrossing. Additionally, the magnitude of outcrossing rates in the mixed mating populations may have been sufficient to facilitate adaptation to S. marcescens in a manner very similar to the obligately outcrossing populations. Perhaps only moderate amounts of outcrossing were necessary to facilitate the evolution of increased host defense. Indeed, our results support the prediction that mixed mating populations can gain at least some of the benefits associated with obligate outcrossing [11, 41, 42].

Although the obligately outcrossing and mixed mating populations evolved greater levels of host defense in general, many different mechanisms may contribute to these changes [1]. Resistance, tolerance, and avoidance are just three potential mechanisms that could underlie the increases in host defense. However multiple traits and many different genetic architectures may contribute to these mechanisms. Further, increased defense may have evolved through different traits in separate populations. Both behavioral [30–33] and cellular [43–45] traits are known to facilitate defense in C. elegans. Additionally, many loci are known to contribute to these traits [44–46]. Therefore, the specific traits underlying defense and the loci that contributed to the evolution of defense may have been significantly influenced by the host mating system. Nonetheless, strong selection, like that imposed by virulent parasites, may often supersede the direct effects of the host mating system on a host population’s evolutionary trajectory and promote parallel evolution. Further work is necessary to identify the components of defense in our experimental populations and more precisely determine the influence of the host mating system on their evolutionary trajectories.

Supporting information

[Figure omitted. See PDF.]

S1 File. Fitness data.

Raw competitive fitness data from the ancestral and experimental populations.

https://doi.org/10.1371/journal.pone.0181913.s001

(XLS)

S2 File. Host mortality data.

Raw host mortality assay data from the ancestral and experimental population.

https://doi.org/10.1371/journal.pone.0181913.s002

(XLS)

S3 File. Host fecundity data.

Raw host fecundity data over multiple days of reproduction presented as daily and total sums of offspring from ancestral and experimental populations.

https://doi.org/10.1371/journal.pone.0181913.s003

(XLS)

Acknowledgments

We thank Sarah Hesse for logistical assistance. We also thank the members of the Morran lab, the GerMo HideR labs, and two anonymous reviewers for comments that improved this work. Some nematode strains used in this work were provided by the Caenorhabditis Genetics Center, which is funded by the NIH National Center for Research Resources (NCRR).

Author Contributions

1. Conceptualization: MJP LTM.

2. Data curation: MJP.

3. Formal analysis: LTM.

4. Funding acquisition: LTM.

5. Investigation: MJP GTH.

6. Methodology: MJP GTH LTM.

7. Project administration: MJP LTM.

8. Resources: MJP LTM.

9. Supervision: LTM.

10. Validation: MJP LTM.

11. Visualization: MJP LTM.

12. Writing – original draft: MJP LTM.

13. Writing – review & editing: MJP GTH LTM.

Citation: Penley MJ, Ha GT, Morran LT (2017) Evolution of Caenorhabditis elegans host defense under selection by the bacterial parasite Serratia marcescens. PLoS ONE 12(8): e0181913. https://doi.org/10.1371/journal.pone.0181913

References

1. Parker BJ, Barribeau SM, Laughton AM, de Roode JC, Gerardo NM. Non-immunological defense in an evolutionary framework. Trends Ecol Evol. 2011;26(5):242–8. pmid:21435735.

2. de Roode JC, Lefevre T, Hunter MD. Ecology. Self-medication in animals. Science. 2013;340(6129):150–1. pmid:23580516.

3. Schmid-Hempel P. Immune defence, parasite evasion strategies and their relevance for 'macroscopic phenomena' such as virulence. Philos Trans R Soc Lond B Biol Sci. 2009;364(1513):85–98. Epub 2008/10/22. pmid:18930879;

4. Hily JM, Poulicard N, Mora MA, Pagan I, Garcia-Arenal F. Environment and host genotype determine the outcome of a plant-virus interaction: from antagonism to mutualism. New Phytol. 2016;209(2):812–22. pmid:26365599.

5. Agnew P, CK J, Michalakis Y. Host life history responses to parasitism. Microbes Infect. 2000;2(8):891–6. pmid:10962272.

6. Michalakis Y, Hochberg ME. Parasitic effects on host life-history traits: a review of recent studies. Parasite. 1994;1(4):291–4. pmid:9140497.

7. Chadwick W, Little TJ. A parasite-mediated life-history shift in Daphnia magna. Proc Biol Sci. 2005;272(1562):505–9. pmid:15799946;

8. Polak M, Starmer WT. Parasite-induced risk of mortality elevates reproductive effort in male Drosophila. Proc Biol Sci. 1998;265(1411):2197–201. pmid:9872007;

9. Schmid-Hempel P. Evolutionary ecology of insect immune defenses. Annu Rev Entomol. 2005;50:529–51. Epub 2004/10/09. pmid:15471530.

10. Morran LT, Schmidt OG, Gelarden IA, Parrish RC 2nd, Lively CM. Running with the Red Queen: host-parasite coevolution selects for biparental sex. Science. 2011;333(6039):216–8. Epub 2011/07/09. pmid:21737739.

11. Charlesworth D, Charlesworth B. Inbreeding depression and its evolutionary consequences. Annu Rev Ecol Syst. 1987;18:237–68. ISI:A1987K958800011.

12. Charlesworth D, Morgan MT, Charlesworth B. Mutation accumulation in finite, outbreeding, and inbreeding populations. Genetics Research. 1993;61:39–56.

13. Lande R, Schemske DW. The evolution of self-fertilization and inbreeding depression in plants .1. genetic models. Evolution 1985;39(1):24–40. ISI:A1985ACA4300003. pmid:28563655

14. Goldberg EE, Kohn JR, Lande R, Robertson KA, Smith SA, Igic B. Species selection maintains self-incompatibility. Science. 2010;330:493–5. pmid:20966249

15. Heller J, Maynard Smith J. Does Muller's Ratchet work with selfing? Genetics Research. 1979;8:269–94.

16. Stebbins GL. Self-fertilization and population variation in higher plants. American Naturalist. 1957;91:337–54.

17. Lynch M, Conery J, Burger R. Mutational meltdowns in sexual populations. Evolution. 1995;47:1744–57.

18. Teotonio H, Carvalho S, Manoel D, Roque M, Chelo IM. Evolution of outcrossing in experimental populations of Caenorhabditis elegans. PLoS ONE. 2012;7(4):e35811. pmid:22540006

19. Loewe L, Cutter AD. On the potential for extinction by Muller's ratchet in Caenorhabditis elegans. BMC Evol Biol. 2008;8:125. Epub 2008/05/02. pmid:18447910.

20. Morran LT, Parmenter MD, Phillips PC. Mutation load and rapid adaptation favour outcrossing over self-fertilization. Nature. 2009;462(7271):350–2. Epub 2009/10/23. pmid:19847164.

21. Hill WG, Robertson A. The effect of linkage on the limits to artificial selection. Genet Res. 1966;8:269–94. pmid:5980116

22. Felsenstein J. The evolutionary advantage of recombination. Genetics. 1974;78:737–56. pmid:4448362

23. Fisher RA. The Genetical Theory of Natural Selection. Oxford: Clarendon Press; 1930.

24. Muller HJ. Some genetic aspects of sex. American Naturalist. 1932;66:118–38.

25. Barton NH. Linkage and the limits to natural selection. Genetics. 1995;140:821–41. pmid:7498757

26. Keightley PD, Otto SP. Interference among deleterious mutations favours sex and recombination in finite populations. Nature. 2006;443(7107):89–92. pmid:16957730.

27. Gandon S, Otto SP. The evolution of sex and recombination in response to abiotic or coevolutionary fluctuations in epistasis. Genetics. 2007;175(4):1835–53. pmid:17277371

28. Roze D, Barton NH. The Hill-Robertson effect and the evolution of recombination. Genetics. 2006;173:1793–811. pmid:16702422

29. Brenner S. The genetics of Caenorhabditis elegans. Genetics. 1974;77(1):71–94. pmid:4366476.

30. Schedl T, Kimble J. fog-2, a germ-line-specific sex determination gene required for hermaphrodite spermatogenesis in Caenorhabditis elegans. Genetics. 1988;119(1):43–61. pmid:3396865.

31. Miller LM, Plenefisch JD, Casson LP, Meyer BJ. xol-1—a gene that controls the male modes of both sex determination and X-chromosome dosage compensation in C. elegans. Cell. 1988;55(1):167–83. ISI:A1988Q510600018. pmid:3167975

32. Schulenburg H, Ewbank JJ. Diversity and specificity in the interaction between Caenorhabditis elegans and the pathogen Serratia marcescens. BMC Evol Biol. 2004;4(1):49. Epub 2004/11/24. pmid:15555070;

33. Morran LT, Parrish RC 2nd, Gelarden IA, Lively CM. Temporal Dynamics of Outcrossing and Host Mortality Rates in Host-Pathogen Experimental Coevolution. Evolution. 2013;67:1860–8. pmid:23815644

34. Morran LT, Parrish RC 2nd, Gelarden IA, Allen MB, Lively CM. Experimental coevolution: rapid local adaptation by parasites depends on host mating system Am Nat. 2014;184 Suppl 1:S91–S100. Epub 2014/07/26. pmid:25061681.

35. Parrish RC 2nd, Penley MJ, Morran LT. The Integral Role of Genetic Variation in the Evolution of Outcrossing in the Caenorhabditis elegans-Serratia marcescens Host-Parasite System. PLoS One. 2016;11(4):e0154463. pmid:27119159;

36. Lively CM, Morran LT. The ecology of sexual reproduction. J Evol Biol. 2014;27(7):1292–303. Epub 2014/03/13. pmid:24617324;

37. Hartfield M, Keightley PD. Current hypotheses for the evolution of sex and recombination. Integr Zool. 2012;7(2):192–209. Epub 2012/06/14. pmid:22691203.

38. Long TA, Agrawal AF, Rowe L. The effect of sexual selection on offspring fitness depends on the nature of genetic variation. Curr Biol. 2012;22(3):204–8. pmid:22226747.

39. Whitlock MC, Agrawal AF. Purging the genome with sexual selection: reducing mutation load through selection on males. Evolution. 2009;63(3):569–82. pmid:19154364.

40. Maclellan K, Whitlock MC, Rundle HD. Sexual selection against deleterious mutations via variable male search success. Biol Lett. 2009;5(6):795–7. pmid:19625301;

41. Goodwillie C, Kalisz S, Ekert C. The evolutionary enigma of mixed mating systems in plants: Occurrence, theoretical explanations, and empirical evidence. Annu Rev Ecol Evol Syst. 2005;36:47–79.

42. Uyenoyama MK, Waller DM. Coevolution of self-fertilization and inbreeding depression. I. Mutation-selection balance at one and two loci. Theor Popul Biol. 1991;40(1):14–46. pmid:1948770.

43. Mallo GV, Kurz CL, Couillault C, Pujol N, Granjeaud S, Kohara Y, et al. Inducible antibacterial defense system in C. elegans. Current Biology. 2002;12(14):1209–14. Epub 2002/08/15. pmid:12176330.

44. Wong D, Bazopoulou D, Pujol N, Tavernarakis N, Ewbank J. Genome-wide investigation reveal pathogen-specific and shared signatures in response to Caenorhabditis elegans to infection. Genome Biology. 2007;8(9):R194. R194. pmid:17875205

45. Englemann I, Pujol N. Innate Immunity in C. elegans. In: Söderhäll K, editor. Invertebrate Immunity. Boston, MA: Landes Bioscience and Spring Science+Business Media; 2010.

46. Glater EE, Rockman MV, Bargmann CI. Multigenic natural variation underlies Caenorhabditis elegans olfactory preference for the bacterial pathogen Serratia marcescens. G3. 2014;4(2):265–76. Epub 2013/12/19. pmid:24347628;