Parasite diversity drives rapid host dynamics and evolution of resistance in a bacteria‐phage system

Host–parasite evolutionary interactions are typically considered in a pairwise species framework. However, natural infections frequently involve multiple parasites. Altering parasite diversity alters ecological and evolutionary dynamics as parasites compete and hosts resist multiple infection. We investigated the effects of parasite diversity on host–parasite population dynamics and evolution using the pathogen Pseudomonas aeruginosa and five lytic bacteriophage parasites. To manipulate parasite diversity, bacterial populations were exposed for 24 hours to either phage monocultures or diverse communities containing up to five phages. Phage communities suppressed host populations more rapidly but also showed reduced phage density, likely due to interphage competition. The evolution of resistance allowed rapid bacterial recovery that was greater in magnitude with increases in phage diversity. We observed no difference in the extent of resistance with increased parasite diversity, but there was a profound impact on the specificity of resistance; specialized resistance evolved to monocultures through mutations in a diverse set of genes. In summary, we demonstrate that parasite diversity has rapid effects on host–parasite population dynamics and evolution by selecting for different resistance mutations and affecting the magnitude of bacterial suppression and recovery. Finally, we discuss the implications of phage diversity for their use as biological control agents.     

Do‐or‐die life cycles and diverse post‐infection resistance mechanisms limit the evolution of parasite host ranges

In light of the dynamic nature of parasite host ranges and documented potential for rapid host shifts, the observed high host specificity of most parasites remains an ecological paradox. Different variants of host‐use trade‐offs have become a mainstay of theoretical explanations of the prevalence of host specialism, but empirical evidence for such trade‐offs is rare. We propose an alternative theory based on basic features of the parasite life cycle: host selection and subsequent intrahost replication. We introduce a new concept of effective burst size that accounts for the fact that successful host selection does not guarantee intrahost replication. Our theory makes a general prediction that a parasite will expand its host range if its effective burst size is positive. An in silico model of bacteria‐phage coevolution verifies our predictions and demonstrates that the tendency for relatively narrow host ranges in parasites can be explained even in the absence of trade‐offs.     

Conflicting selection alters the trajectory of molecular evolution in a tripartite bacteria–plasmid–phage interaction

Bacteria engage in a complex network of ecological interactions, which includes mobile genetic elements (MGEs) such as phages and plasmids. These elements play a key role in microbial communities as vectors of horizontal gene transfer but can also be important sources of selection for their bacterial hosts. In natural communities, bacteria are likely to encounter multiple MGEs simultaneously and conflicting selection among MGEs could alter the bacterial evolutionary response to each MGE. Here, we test the effect of interactions with multiple MGEs on bacterial molecular evolution in the tripartite interaction between the bacterium, Pseudomonas fluorescens, the lytic bacteriophage, SBW25φ2, and conjugative plasmid, pQBR103, using genome sequencing of experimentally evolved bacteria. We show that individually, both plasmids and phages impose selection leading to bacterial evolutionary responses that are distinct from bacterial populations evolving without MGEs, but that together, plasmids and phages impose conflicting selection on bacteria, constraining the evolutionary responses observed in pairwise interactions. Our findings highlight the likely difficulties of predicting evolutionary responses to multiple selective pressures from the observed evolutionary responses to each selective pressure alone. Understanding evolution in complex microbial communities comprising many species and MGEs will require that we go beyond studies of pairwise interactions.     

The cost of phage resistance in a plant pathogenic bacterium is context‐dependent

Parasites are ubiquitous features of living systems and many parasites severely reduce the fecundity or longevity of their hosts. This parasite‐imposed selection on host populations should strongly favor the evolution of host resistance, but hosts typically face a trade‐off between investment in reproductive fitness and investment in defense against parasites. The magnitude of such a trade‐off is likely to be context‐dependent, and accordingly costs that are key in shaping evolution in nature may not be easily observable in an artificial environment. We set out to assess the costs of phage resistance for a plant pathogenic bacterium in its natural plant host versus in a nutrient‐rich, artificial medium. We demonstrate that mutants of Pseudomonas syringae that have evolved resistance via a single mutational step pay a substantial cost for this resistance when grown on their tomato plant hosts, but do not realize any measurable growth rate costs in nutrient‐rich media. This work demonstrates that resistance to phage can significantly alter bacterial growth within plant hosts, and therefore that phage‐mediated selection in nature is likely to be an important component of bacterial pathogenicity.     

Turnover in local parasite populations temporarily favors host outcrossing over self‐fertilization during experimental evolution

The ubiquity of outcrossing in plants and animals is difficult to explain given its costs relative to self‐fertilization. Despite these costs, exposure to changing environmental conditions can temporarily favor outcrossing over selfing. Therefore, recurring episodes of environmental change are predicted to favor the maintenance of outcrossing. Studies of host–parasite coevolution have provided strong support for this hypothesis. However, it is unclear whether multiple exposures to novel parasite genotypes in the absence of coevolution are sufficient to favor outcrossing. Using the nematode Caenorhabditis elegans and the bacterial parasite Serratia marcescens, we studied host responses to parasite turnover. We passaged several replicates of a host population that was well‐adapted to the S. marcescens strain Sm2170 with either Sm2170 or one of three novel S. marcescens strains, each derived from Sm2170, for 18 generations. We found that hosts exposed to novel parasites maintained higher outcrossing rates than hosts exposed to Sm2170. Nonetheless, host outcrossing rates declined over time against all but the most virulent novel parasite strain. Hosts exposed to the most virulent novel strain exhibited increased outcrossing rates for approximately 12 generations, but did not maintain elevated levels of outcrossing throughout the experiment. Thus, parasite turnover can transiently increase host outcrossing. These results suggest that recurring episodes of parasite turnover have the potential to favor the maintenance of host outcrossing. However, such maintenance may require frequent exposure to novel virulent parasites, rapid rates of parasite turnover, and substantial host gene flow.     

Coevolutionary dynamics shape the structure of bacteria‐phage infection networks

Coevolution—reciprocal evolutionary change among interacting species driven by natural selection—is thought to be an important force in shaping biodiversity. This ongoing process takes place within tangled networks of species interactions. In microbial communities, evolutionary change between hosts and parasites occurs at the same time scale as ecological change. Yet, we still lack experimental evidence of the role of coevolution in driving changes in the structure of such species interaction networks. Filling this gap is important because network structure influences community persistence through indirect effects. Here, we quantified experimentally to what extent coevolutionary dynamics lead to contrasting patterns in the architecture of bacteria–phage infection networks. Specifically, we look at the tendency of these networks to be organized in a nested pattern by which the more specialist phages tend to infect only a proper subset of those bacteria infected by the most generalist phages. We found that interactions between coevolving bacteria and phages become less nested over time under fluctuating dynamics, and more nested under arms race dynamics. Moreover, when coevolution results in high average infectivity, phages and bacteria differ more from each other over time under arms race dynamics than under fluctuating dynamics. The tradeoff between the fitness benefits of evolving resistance/infectivity traits and the costs of maintaining them might explain these differences in network structure. Our study shows that the interaction pattern between bacteria and phages at the community level depends on the way coevolution unfolds.     

Pre‐adapting parasitic phages to a pathogen leads to increased pathogen clearance and lowered resistance evolution with Pseudomonas aeruginosa cystic fibrosis bacterial isolates

Recent years have seen renewed interest in phage therapy – the use of viruses to specifically kill disease‐causing bacteria – because of the alarming rise in antibiotic resistance. However, a major limitation of phage therapy is the ease at with bacteria can evolve resistance to phages. Here, we determined whether in vitro experimental coevolution can increase the efficiency of phage therapy by limiting the resistance evolution of intermittent and chronic cystic fibrosis Pseudomonas aeruginosa lung isolates to four different phages. We first pre‐adapted all phage strains against all bacterial strains and then compared the efficacy of pre‐adapted and nonadapted phages against ancestral bacterial strains. We found that evolved phages were more efficient in reducing bacterial densities than ancestral phages. This was primarily because only 50% of bacterial strains were able to evolve resistance to evolved phages, whereas all bacteria were able to evolve some level of resistance to ancestral phages. Although the rate of resistance evolution did not differ between intermittent and chronic isolates, it incurred a relatively higher growth cost for chronic isolates when measured in the absence of phages. This is likely to explain why evolved phages were more effective in reducing the densities of chronic isolates. Our data show that pathogen genotypes respond differently to phage pre‐adaptation, and as a result, phage therapies might need to be individually adjusted for different patients.     

No effect of host–parasite co‐evolution on host range expansion

Antagonistic co‐evolution between hosts and parasites (reciprocal selection for resistance and infectivity) is hypothesized to play an important role in host range expansion by selecting for novel infectivity alleles, but tests are lacking. Here, we determine whether experimental co‐evolution between a bacterium (Pseudomonas fluorescens SBW25) and a phage (SBW25Φ2) affects interstrain host range: the ability to infect different strains of P. fluorescens other than SBW25. We identified and tested a genetically and phenotypically diverse suite of co‐evolved phage variants of SBW25Φ2 against both sympatric and allopatric co‐evolving hosts (P. fluorescens SBW25) and a large set of other P. fluorescens strains. Although all co‐evolved phage had a greater host range than the ancestral phage and could differentially infect co‐evolved variants of P. fluorescens SBW25, none could infect any of the alternative P. fluorescens strains. Thus, parasite generalism at one genetic scale does not appear to affect generalism at other scales, suggesting fundamental genetic constraints on parasite adaptation for this virus.     

BREX is a novel phage resistance system widespread in microbial genomes

The perpetual arms race between bacteria and phage has resulted in the evolution of efficient resistance systems that protect bacteria from phage infection. Such systems, which include the CRISPR‐Cas and restriction‐modification systems, have proven to be invaluable in the biotechnology and dairy industries. Here, we report on a six‐gene cassette in Bacillus cereus which, when integrated into the Bacillus subtilis genome, confers resistance to a broad range of phages, including both virulent and temperate ones. This cassette includes a putative Lon‐like protease, an alkaline phosphatase domain protein, a putative RNA‐binding protein, a DNA methylase, an ATPase‐domain protein, and a protein of unknown function. We denote this novel defense system BREX (Bacteriophage Exclusion) and show that it allows phage adsorption but blocks phage DNA replication. Furthermore, our results suggest that methylation on non‐palindromic TAGGAG motifs in the bacterial genome guides self/non‐self discrimination and is essential for the defensive function of the BREX system. However, unlike restriction‐modification systems, phage DNA does not appear to be cleaved or degraded by BREX, suggesting a novel mechanism of defense. Pan genomic analysis revealed that BREX and BREX‐like systems, including the distantly related Pgl system described in Streptomyces coelicolor, are widely distributed in ~10% of all sequenced microbial genomes and can be divided into six coherent subtypes in which the gene composition and order is conserved. Finally, we detected a phage family that evades the BREX defense, implying that anti‐BREX mechanisms may have evolved in some phages as part of their arms race with bacteria.     

Experimental evolution with a multicellular host causes diversification within and between microbial parasite populations—Differences in emerging phenotypes of two different parasite strains

Host‐parasite coevolution is predicted to have complex evolutionary consequences, potentially leading to the emergence of genetic and phenotypic diversity for both antagonists. However, little is known about variation in phenotypic responses to coevolution between different parasite strains exposed to the same experimental conditions. We infected Caenorhabditis elegans with one of two strains of Bacillus thuringiensis and either allowed the host and the parasite to experimentally coevolve (coevolution treatment) or allowed only the parasite to adapt to the host (one‐sided parasite adaptation). By isolating single parasite clones from evolved populations, we found phenotypic diversification of the ancestral strain into distinct clones, which varied in virulence toward ancestral hosts and competitive ability against other parasite genotypes. Parasite phenotypes differed remarkably not only between the two strains, but also between and within different replicate populations, indicating diversification of the clonal population caused by selection. This study highlights that the evolutionary selection pressure mediated by a multicellular host causes phenotypic diversification, but not necessarily with the same phenotypic outcome for different parasite strains.     

Parasite host range and the evolution of host resistance

Parasite host range plays a pivotal role in the evolution and ecology of hosts and the emergence of infectious disease. Although the factors that promote host range and the epidemiological consequences of variation in host range are relatively well characterized, the effect of parasite host range on host resistance evolution is less well understood. In this study, we tested the impact of parasite host range on host resistance evolution. To do so, we used the host bacterium Pseudomonas fluorescens SBW25 and a diverse suite of coevolved viral parasites (lytic bacteriophage Φ2) with variable host ranges (defined here as the number of host genotypes that can be infected) as our experimental model organisms. Our results show that resistance evolution to coevolved phages occurred at a much lower rate than to ancestral phage (approximately 50% vs. 100%), but the host range of coevolved phages did not influence the likelihood of resistance evolution. We also show that the host range of both single parasites and populations of parasites does not affect the breadth of the resulting resistance range in a naïve host but that hosts that evolve resistance to single parasites are more likely to resist other (genetically) more closely related parasites as a correlated response. These findings have important implications for our understanding of resistance evolution in natural populations of bacteria and viruses and other host–parasite combinations with similar underlying infection genetics, as well as the development of phage therapy.     

Genomic evidence for population‐specific responses to co‐evolving parasites in a New Zealand freshwater snail

Reciprocal co‐evolving interactions between hosts and parasites are a primary source of strong selection that can promote rapid and often population‐ or genotype‐specific evolutionary change. These host–parasite interactions are also a major source of disease. Despite their importance, very little is known about the genomic basis of co‐evolving host–parasite interactions in natural populations, especially in animals. Here, we use gene expression and sequence evolution approaches to take critical steps towards characterizing the genomic basis of interactions between the freshwater snail Potamopyrgus antipodarum and its co‐evolving sterilizing trematode parasite, Microphallus sp., a textbook example of natural coevolution. We found that Microphallus‐infected P. antipodarum exhibit systematic downregulation of genes relative to uninfected P. antipodarum. The specific genes involved in parasite response differ markedly across lakes, consistent with a scenario where population‐level co‐evolution is leading to population‐specific host–parasite interactions and evolutionary trajectories. We also used an F_ST‐based approach to identify a set of loci that represent promising candidates for targets of parasite‐mediated selection across lakes as well as within each lake population. These results constitute the first genomic evidence for population‐specific responses to co‐evolving infection in the P. antipodarum‐Microphallus interaction and provide new insights into the genomic basis of co‐evolutionary interactions in nature.     

Inter‐ and intraspecific variation of spider mite susceptibility to fungal infections: Implications for the long‐term success of biological control

Spider mites are severe pests of several annual and perennial crops worldwide, often causing important economic damages. As rapid evolution of pesticide resistance in this group hampers the efficiency of chemical control, alternative control strategies, such as the use of entomopathogenic fungi, are being developed. However, while several studies have focused on the evaluation of the control potential of different fungal species and/or isolates as well as their compatibility with other control methods (e.g., predators or chemical pesticides), knowledge on the extent of inter‐ and intraspecific variation in spider mite susceptibility to fungal infection is as yet incipient. Here, we measured the mortality induced by two generalist fungi, Beauveria bassiana and Metarhizium brunneum, in 12 spider mite populations belonging to different Tetranychus species: T. evansi, T. ludeni, and T. urticae (green and red form), within a full factorial experiment. We found that spider mite species differed in their susceptibility to infection by both fungal species. Moreover, we also found important intraspecific variation for this trait. These results draw caution on the development of single strains as biocontrol agents. Indeed, the high level of intraspecific variation suggests that (a) the one‐size‐fits‐all strategy may fail to control spider mite populations and (b) hosts resistance to infection may evolve at a rapid pace. Finally, we propose future directions to better understand this system and improve the long‐term success of spider mite control strategies based on entomopathogenic fungi.     

Gene copy number variations as signatures of adaptive evolution in the parthenogenetic,plant‐parasitic nematode Meloidogyne incognita

Adaptation to changing environmental conditions represents a challenge to parthenogenetic organisms, and until now, how phenotypic variants are generated in clones in response to the selection pressure of their environment remains poorly known. The obligatory parthenogenetic root‐knot nematode species Meloidogyne incognita has a worldwide distribution and is the most devastating plant‐parasitic nematode. Despite its asexual reproduction, this species exhibits an unexpected capacity of adaptation to environmental constraints, for example, resistant hosts. Here, we used a genomewide comparative hybridization strategy to evaluate variations in gene copy numbers between genotypes of M. incognita resulting from two parallel experimental evolution assays on a susceptible vs. resistant host plant. We detected gene copy number variations (CNVs) associated with the ability of the nematodes to overcome resistance of the host plant, and this genetic variation may reflect an adaptive response to host resistance in this parthenogenetic species. The CNV distribution throughout the nematode genome is not random and suggests the occurrence of genomic regions more prone to undergo duplications and losses in response to the selection pressure of the host resistance. Furthermore, our analysis revealed an outstanding level of gene loss events in nematode genotypes that have overcome the resistance. Overall, our results support the view that gene loss could be a common class of adaptive genetic mechanism in response to a challenging new biotic environment in clonal animals.     

A plant‐specific in vitro ubiquitination analysis system

Protein ubiquitination requires the concerted action of three enzymes: ubiquitin‐activating enzyme (E1), ubiquitin‐conjugating enzyme (E2) and ubiquitin ligase (E3). These ubiquitination enzymes belong to an abundant protein family that is encoded in all eukaryotic genomes. Describing their biochemical characteristics is an important part of their functional analysis. It has been recognized that various E2/E3 specificities exist, and that detection of E3 ubiquitination activity in vitro may depend on the recruitment of E2s. Here, we describe the development of an in vitro ubiquitination system based on proteins encoded by genes from Arabidopsis. It includes most varieties of Arabidopsis E2 proteins, which are tested with several RING‐finger type E3 ligases. This system permits determination of E3 activity in combination with most of the E2 sub‐groups that have been identified in the Arabidopsis genome. At the same time, E2/E3 specificities have also been explored. The components used in this system are all from plants, particularly Arabidopsis, making it very suitable for ubiquitination assays of plant proteins. Some E2 proteins that are not easily expressed in Escherichia coli were transiently expressed and purified from plants before use in ubiquitination assays. This system is also adaptable to proteins of species other than plants. In this system, we also analyzed two mutated forms of ubiquitin, K48R and K63R, to detect various types of ubiquitin conjugation.     

Diversity in CRISPR‐based immunity protects susceptible genotypes by restricting phage spread and evolution

Diversity in host resistance often associates with reduced pathogen spread. This may result from ecological and evolutionary processes, likely with feedback between them. Theory and experiments on bacteria–phage interactions have shown that genetic diversity of the bacterial adaptive immune system can limit phage evolution to overcome resistance. Using the CRISPR–Cas bacterial immune system and lytic phage, we engineered a host–pathogen system where each bacterial host genotype could be infected by only one phage genotype. With this model system, we explored how CRISPR diversity impacts the spread of phage when they can overcome a resistance allele, how immune diversity affects the evolution of the phage to increase its host range and if there was feedback between these processes. We show that increasing CRISPR diversity benefits susceptible bacteria via a dilution effect, which limits the spread of the phage. We suggest that this ecological effect impacts the evolution of novel phage genotypes, which then feeds back into phage population dynamics.     

Evolution of bacterial life‐history traits is sensitive to community structure

Very few studies have experimentally assessed the evolutionary effects of species interactions within the same trophic level. Here we show that when Serratia marcescens evolve in multispecies communities, their growth rate exceeds the growth rate of the bacteria that evolved alone, whereas the biomass yield gets lower. In addition to the community effects per se, we found that few species in the communities caused strong effects on S. marcescens evolution. The results indicate that evolutionary responses (of a focal species) are different in communities, compared to species evolving alone. Moreover, selection can lead to very different outcomes depending on the community structure. Such context dependencies cast doubt on our ability to predict the course of evolution in the wild, where species often inhabit very different kinds of communities.     

The experimental evolution of herbicide resistance in Chlamydomonas reinhardtii results in a positive correlation between fitness in the presence and absence of herbicides

Pleiotropic fitness trade-offs will be key determinants of the evolutionary dynamics of selection for pesticide resistance. However, for herbicide resistance, empirical support for a fitness cost of resistance is mixed, and it is therefore also questionable what further ecological trade-offs can be assumed to apply to herbicide resistance. Here, we test the existence of trade-offs by experimentally evolving herbicide resistance in Chlamydomonas reinhardtii. Although fitness costs are detected for all herbicides, we find that, counterintuitively, the most resistant populations also have the lowest fitness costs as measured by growth rate in the ancestral environment. Furthermore, after controlling for differences in the evolutionary dynamics of resistance to different herbicides, we also detect significant positive correlations between resistance, fitness in the ancestral environment and cross-resistance to other herbicides. We attribute this to the highest levels of nontarget-site resistance being achieved by fixing mutations that more broadly affect cellular physiology, which results in both more cross-resistance and less overall antagonistic pleiotropy on maximum growth rate. Consequently, the lack of classical ecological trade-offs could present a major challenge for herbicide resistance management.