The mode of host-parasite interaction shapes coevolutionary dynamics and the fate of host cooperation

Eggshell colour differences in a classic example of coevolved eggshell mimicry

Avian brood parasitism is a model system for understanding coevolutionary arms races, and the great reed warbler (Acrocephalus arundinaceus, hereafter 'warbler') and its parasite the common cuckoo (Cuculus canorus, hereafter 'cuckoo') are prime examples of this coevolutionary struggle. Here, warblers select for egg colour mimicry by rejecting poorly matched cuckoo eggs. Contrary to long-held assumptions, recent work showed that warblers tend to reject lighter and browner eggs but tended to accept darker and bluer eggs rather than basing rejection decisions solely on perceived colour differences (i.e. the degree of mimicry). This counterintuitive, colour-biased rejection behaviour would select for bluer and darker cuckoo eggs, but would only be adaptive if cuckoos were consistently lighter and browner than warbler eggs. Therefore, we tested whether warbler eggs were consistently bluer and darker than cuckoo eggs. To do so, we re-analysed eggshell reflectance spectra of warblers and the cuckoos that parasitized them in the Czech Republic. As expected, we found that warbler eggs were significantly bluer and darker than the cuckoo eggs at the population level. Thus, we demonstrate imperfect mimicry in a long-coevolved cuckoo host-race and provide insights for exploring the coevolutionary interactions among hosts and their brood parasites.

On the difficult evolutionary transition from the free-living lifestyle to obligate symbiosis

Obligate symbiosis evolved from free-living individuals most likely via the intermediate stage of facultative symbiosis. However, why should facultative symbionts, who can live independently but also benefit from their partners if these are available, give up this best of both worlds? Using the adaptive dynamics approach, we analyse a simple model, focusing on one partner of the symbiosis, to gain more insight into the selective forces that make individuals forgo the ability to reproduce in the free-living state. Our results suggest that, similar to the parasitism-mutualism continuum, the free-living way of life and obligate symbiosis are two extremes of a continuum of the ability to reproduce independently of a partner. More importantly, facultative symbiosis should be the rule as for many parameter combinations completely giving up independent reproduction or adopting a pure free-living strategy is not so easy. We also show that if host encounter comes at a cost, individuals that put more effort into increasing the chances to meet with their partners are more likely to give up the ability to reproduce independently. Finally, our model does not specify the ecological interactions between hosts and symbionts but we discuss briefly how the ecological nature of an interaction can influence the transition from facultative to obligate symbiosis.

Insights into the Natural Defenses of a Coral Reef Fish Against Gill Ectoparasites: Integrated Metabolome and Microbiome Approach

Understanding natural defense mechanisms against parasites can be a valuable tool for the development of innovative therapies. We have previously identified a butterflyfish species (Chaetodon lunulatus) that avoids gill monogenean parasites while living amongst closely related parasitized species. The metabolome and microbiome of several sympatric butterflyfish species from the island of Moorea (French Polynesia) were previously described. In this study, we used the previously generated datasets in an attempt to identify metabolites and bacteria potentially involved in parasite defense mechanisms. We investigated the interplay between the gill mucus metabolome and microbiome of the non-susceptible C. lunulatus versus sympatric butterflyfish species that were always found parasitized in the Central and Eastern Indo-Pacific. After observing significant differences between the metabolome and bacteria of susceptible versus non-susceptible fish, we obtained the discriminant metabolites and operational taxonomic units (OTUs) using a supervised analysis. Some of the most important discriminant metabolites were identified as peptides, and three new peptides derived from β-subunit hemoglobin from C. lunulatus (CLHbβ-1, CLHbβ-2, and CLHbβ-3) were purified, characterized and synthesized to confirm their structures. We also identified specific bacterial families and OTUs typical from low-oxygen habitats in C. lunulatus gill mucus. By using a correlation network between the two datasets, we found a Fusobacteriaceae strain exclusively present in C. lunulatus and highly correlated to the peptides. Finally, we discuss the possible involvement of these peptides and Fusobacteriaceae in monogenean avoidance by this fish species.

CRISPR evolution and bacteriophage persistence in the context of population bottlenecks

Population bottlenecks often cause strong reductions in genetic diversity and alter population structure. In the context of host-parasite interactions, bottlenecks could in theory benefit either the host or the pathogen. We predicted that bottlenecking of bacterial populations that evolve CRISPR immunity against bacteriophages (phage) would benefit the pathogen, because CRISPR spacer diversity can rapidly drive phages extinct. To test this, we bottlenecked populations of bacteria and phage, tracking phage persistence and the evolution of bacterial resistance mechanisms. Contrary to our prediction, bottlenecking worked in the advantage of the host. With some possible exceptions, this effect was not caused by CRISPR immunity. This host benefit is consistent with a dilution effect disproportionately affecting phage. This study provides further insight into how bottlenecking influences bacteria-phage dynamics, the role of dilution in bacteria-phage interactions, and the evolution of host immune systems.

Genomic characterization of three novel Basilisk-like phages infecting Bacillus anthracis

Background

In the present study, we sequenced the complete genomes of three novel bacteriophages v_B-Bak1, v_B-Bak6, v_B-Bak10 previously isolated from historical anthrax burial sites in the South Caucasus country of Georgia. We report here major trends in the molecular evolution of these phages, which we designate as “Basilisk-Like-Phages” (BLPs), and illustrate patterns in their evolution, genomic plasticity and core genome architecture.

Results

Comparative whole genome sequence analysis revealed a close evolutionary relationship between our phages and two unclassified Bacillus cereus group phages, phage Basilisk, a broad host range phage (Grose JH et al., J Vir. 2014;88(20):11846-11860) and phage PBC4, a highly host-restricted phage and close relative of Basilisk (Na H. et al. FEMS Microbiol. letters. 2016;363(12)). Genome comparisons of phages v_B-Bak1, v_B-Bak6, and v_B-Bak10 revealed significant similarity in sequence, gene content, and synteny with both Basilisk and PBC4. Transmission electron microscopy (TEM) confirmed the three phages belong to the Siphoviridae family. In contrast to the broad host range of phage Basilisk and the single-strain specificity of PBC4, our three phages displayed host specificity for Bacillus anthracis. Bacillus species including Bacillus cereus, Bacillus subtilis, Bacillus anthracoides, and Bacillus megaterium were refractory to infection.

Conclusions

Data reported here provide further insight into the shared genomic architecture, host range specificity, and molecular evolution of these rare B. cereus group phages. To date, the three phages represent the only known close relatives of the Basilisk and PBC4 phages and their shared genetic attributes and unique host specificity for B. anthracis provides additional insight into candidate host range determinants.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5056-4) contains supplementary material, which is available to authorized users.

Host-parasite fluctuating selection in the absence of specificity

Fluctuating selection driven by coevolution between hosts and parasites is important for the generation of host and parasite diversity across space and time. Theory has focused primarily on infection genetics, with highly specific ‘matching-allele’ frameworks more likely to generate fluctuating selection dynamics (FSD) than ‘gene-for-gene’ (generalist–specialist) frameworks. However, the environment, ecological feedbacks and life-history characteristics may all play a role in determining when FSD occurs. Here, we develop eco-evolutionary models with explicit ecological dynamics to explore the ecological, epidemiological and host life-history drivers of FSD. Our key result is to demonstrate for the first time, to our knowledge, that specificity between hosts and parasites is not required to generate FSD. Furthermore, highly specific host–parasite interactions produce unstable, less robust stochastic fluctuations in contrast to interactions that lack specificity altogether or those that vary from generalist to specialist, which produce predictable limit cycles. Given the ubiquity of ecological feedbacks and the variation in the nature of specificity in host–parasite interactions, our work emphasizes the underestimated potential for host–parasite coevolution to generate fluctuating selection.

Emerging Schistosomiasis in Europe: A Need to Quantify the Risks

The recent recurrent outbreaks of urogenital schistosomiasis in the south of Europe were unanticipated and caught scientists and health authorities unprepared. It is now time to learn lessons from these outbreaks and to implement concrete procedures in order to better quantify the risks and prevent future outbreaks of schistosomiasis in Europe. In this context, we propose a reflection on the factors that currently hamper our ability to quantify these risks and argue that we are incapable of predicting future outbreaks. We base our reflexion on an ecological two-step filter concept that drives host-parasite interactions, namely the encounter and the compatibility filters.

Antibiotic stress selects against cooperation in the pathogenic bacterium Pseudomonas aeruginosa

Cheats are a pervasive threat to public goods production in natural and human communities, as they benefit from the commons without contributing to it. Although ecological antagonisms such as predation, parasitism, competition, and abiotic environmental stress play key roles in shaping population biology, it is unknown how such stresses generally affect the ability of cheats to undermine cooperation. We used theory and experiments to address this question in the pathogenic bacterium, Pseudomonas aeruginosa Although public goods producers were selected against in all populations, our competition experiments showed that antibiotics significantly increased the advantage of nonproducers. Moreover, the dominance of nonproducers in mixed cultures was associated with higher resistance to antibiotics than in either monoculture. Mathematical modeling indicates that accentuated costs to producer phenotypes underlie the observed patterns. Mathematical analysis further shows how these patterns should generalize to other taxa with public goods behaviors. Our findings suggest that explaining the maintenance of cooperative public goods behaviors in certain natural systems will be more challenging than previously thought. Our results also have specific implications for the control of pathogenic bacteria using antibiotics and for understanding natural bacterial ecosystems, where subinhibitory concentrations of antimicrobials frequently occur.

Host and Parasite Evolution in a Tangled Bank

Most hosts and parasites exist in diverse communities wherein they interact with other species, spanning the parasite-mutualist continuum. These additional interactions have the potential to impose selection on hosts and parasites and influence the patterns and processes of their evolution. Yet, host-parasite interactions are almost exclusively studied in species pairs. A wave of new research has incorporated a multispecies community context, showing that additional ecological interactions can alter components of host and parasite fitness, as well as interaction specificity and virulence. Here, we synthesize these findings to assess the effects of increased species diversity on the patterns and processes of host and parasite evolution. We argue that our understanding of host-parasite interactions would benefit from a richer biotic perspective.

Gill monogenean communities (Platyhelminthes, Monogenea, Dactylogyridae) of butterflyfishes from tropical Indo-West Pacific Islands

We studied the monogenean communities of 34 species of butterflyfish from the tropical Indo-West Pacific, identifying 13 dactylogyrid species (including two species that are presently undescribed). Monogenean assemblages differed significantly between host species in terms of taxonomic structure, intensity and prevalence. Parasite richness ranged from 0 (Chaetodon lunulatus) to 11 (C. auriga, C. citrinellus and C. lunula). Host specificity varied between the dactylogyrids species, being found on 2–29 of the 34 chaetodontid species examined. Sympatric butterflyfish species were typically parasitized by different combinations of dactylogyrid species, suggesting the existence of complex host–parasite interactions. We identified six clusters of butterflyfish species based on the similarities of their dactylogyrid communities. Dactylogyrid richness and diversity were not related to host size, diet specialization, depth range or phylogeny of butterflyfish species. However, there was a weak positive correlation between monogenean richness and diversity and host geographical range. Most communities of dactylogyrids were dominated by Haliotrema aurigae and H. angelopterum, indicating the importance of the genus Haliotrema in shaping monogenean communities of butterflyfishes. This study casts light on the structure of the monogenean communities of butterflyfishes, suggesting that the diversity and complexity of community structures arises from a combination of host species-specific parameters.

Eco-evolutionary dynamics of social dilemmas

Social dilemmas are an integral part of social interactions. Cooperative actions, ranging from secreting extra-cellular products in microbial populations to donating blood in humans, are costly to the actor and hence create an incentive to shirk and avoid the costs. Nevertheless, cooperation is ubiquitous in nature. Both costs and benefits often depend non-linearly on the number and types of individuals involved-as captured by idioms such as 'too many cooks spoil the broth' where additional contributions are discounted, or 'two heads are better than one' where cooperators synergistically enhance the group benefit. Interaction group sizes may depend on the size of the population and hence on ecological processes. This results in feedback mechanisms between ecological and evolutionary processes, which jointly affect and determine the evolutionary trajectory. Only recently combined eco-evolutionary processes became experimentally tractable in microbial social dilemmas. Here we analyse the evolutionary dynamics of non-linear social dilemmas in settings where the population fluctuates in size and the environment changes over time. In particular, cooperation is often supported and maintained at high densities through ecological fluctuations. Moreover, we find that the combination of the two processes routinely reveals highly complex dynamics, which suggests common occurrence in nature.

Host-parasite coevolution: why changing population size matters

Host-parasite coevolution is widely assumed to have a major influence on biological evolution, especially as these interactions impose high selective pressure on the reciprocally interacting antagonists. The exact nature of the underlying dynamics is yet under debate and may be determined by recurrent selective sweeps (i.e., arms race dynamics), negative frequency-dependent selection (i.e., Red Queen dynamics), or a combination thereof. These interactions are often associated with reciprocally induced changes in population size, which, in turn, should have a strong impact on co-adaptation processes, yet are neglected in most current work on the topic. Here, we discuss potential consequences of temporal variations in population size on host-parasite coevolution. The limited empirical data available and the current theoretical literature in this field highlight that the consideration of such interaction-dependent population size changes is likely key for the full understanding of the coevolutionary dynamics, and, thus, a more realistic view on the complex nature of species interactions.

Epidemiological and Evolutionary Outcomes in Gene-for-Gene and Matching Allele Models

Gene-for-gene (GFG) and matching-allele (MA) models are qualitatively different paradigms for describing the outcome of genetic interactions between hosts and pathogens. The GFG paradigm was largely built on the foundations of Flor's early work on the flax-flax rust interaction and is based on the concept of genetic recognition leading to incompatible disease outcomes, typical of host immune recognition. In contrast, the MA model is based on the assumption that genetic recognition leads to compatible interactions, which can result when pathogens require specific host factors to cause infection. Results from classical MA and GFG models have led to important predictions regarding various coevolutionary phenomena, including the role of fitness costs associated with resistance and infectivity, the distribution of resistance genes in wild populations, patterns of local adaptation and the evolution and maintenance of sexual reproduction. Empirical evidence (which we review briefly here), particularly from recent molecular advances in understanding of the mechanisms that determine the outcome of host-pathogen encounters, suggests considerable variation in specific details of the functioning of interactions between hosts and pathogens, which may contain elements of both models. In this regard, GFG and MA scenarios likely represent endpoints of a continuum of potentially more complex interactions that occur in nature. Increasingly, this has been recognized in theoretical studies of coevolutionary processes in plant host-pathogen and animal host-parasite associations (e.g., departures from strict GFG/MA assumptions, diploid genetics, multi-step infection processes). However, few studies have explored how different genetic assumptions about host resistance and pathogen infectivity might impact on disease epidemiology or pathogen persistence within and among populations. Here, we use spatially explicit simulations of the basic MA and GFG scenarios to highlight qualitative differences between these scenarios with regard to patterns of disease and impacts on host demography. Given that such impacts drive evolutionary trajectories, future theoretical advances that aim to capture more complex genetic scenarios should explicitly address the interaction between epidemiology and different models of host-pathogen interaction genetics.

Host population bottlenecks drive parasite extinction during antagonistic coevolution

Host-parasite interactions are often characterized by large fluctuations in host population size, and we investigated how such host bottlenecks affected coevolution between a bacterium and a virus. Previous theory suggests that host bottlenecks should provide parasites with an evolutionary advantage, but instead we found that phages were rapidly driven to extinction when coevolving with hosts exposed to large genetic bottlenecks. This was caused by the stochastic loss of sensitive bacteria, which are required for phage persistence and infectivity evolution. Our findings emphasize the importance of feedbacks between ecological and coevolutionary dynamics, and how this feedback can qualitatively alter coevolutionary dynamics.

Host-parasite coevolution in populations of constant and variable size

BACKGROUND

The matching-allele and gene-for-gene models are widely used in mathematical approaches that study the dynamics of host-parasite interactions. Agrawal and Lively (Evolutionary Ecology Research 4:79-90, 2002) captured these two models in a single framework and numerically explored the associated time discrete dynamics of allele frequencies.

RESULTS

Here, we present a detailed analytical investigation of this unifying framework in continuous time and provide a generalization. We extend the model to take into account changing population sizes, which result from the antagonistic nature of the interaction and follow the Lotka-Volterra equations. Under this extension, the population dynamics become most complex as the model moves away from pure matching-allele and becomes more gene-for-gene-like. While the population densities oscillate with a single oscillation frequency in the pure matching-allele model, a second oscillation frequency arises under gene-for-gene-like conditions. These observations hold for general interaction parameters and allow to infer generic patterns of the dynamics.

CONCLUSION

Our results suggest that experimentally inferred dynamical patterns of host-parasite coevolution should typically be much more complex than the popular illustrations of Red Queen dynamics. A single parasite that infects more than one host can substantially alter the cyclic dynamics.

Enhancing and initiating phage-based therapies

Drug development has typically been a primary foundation of strategy for systematic, long-range management of pathogenic cells. However, drug development is limited in speed and flexibility when response is needed to changes in pathogenic cells, especially changes that produce drug-resistance. The high replication speed and high diversity of phages are potentially useful for increasing both response speed and response flexibility when changes occur in either drug resistance or other aspects of pathogenic cells. We present strategy, with some empirical details, for (1) using modern molecular biology and biophysics to access these advantages during the phage therapy of bacterial infections, and (2) initiating use of phage capsid-based drug delivery vehicles (DDVs) with procedures that potentially overcome both drug resistance and other present limitations in the use of DDVs for the therapy of neoplasms. The discussion of phage therapy includes (a) historical considerations, (b) changes that appear to be needed in clinical tests if use of phage therapy is to be expanded, (c) recent work on novel phages and its potential use for expanding the capabilities of phage therapy and (d) an outline for a strategy that encompasses both theory and practice for expanding the applications of phage therapy. The discussion of DDVs starts by reviewing current work on DDVs, including work on both liposomal and viral DDVs. The discussion concludes with some details of the potential use of permeability constrained phage capsids as DDVs.

Coevolutionary diversification creates nested-modular structure in phage-bacteria interaction networks

Phage and their bacterial hosts are the most diverse and abundant biological entities in the oceans, where their interactions have a major impact on marine ecology and ecosystem function. The structure of interaction networks for natural phage-bacteria communities offers insight into their coevolutionary origin. At small phylogenetic scales, observed communities typically show a nested structure, in which both hosts and phages can be ranked by their range of resistance and infectivity, respectively. A qualitatively different multi-scale structure is seen at larger phylogenetic scales; a natural assemblage sampled from the Atlantic Ocean displays large-scale modularity and local nestedness within each module. Here, we show that such 'nested-modular' interaction networks can be produced by a simple model of host-phage coevolution in which infection depends on genetic matching. Negative frequency-dependent selection causes diversification of hosts (to escape phages) and phages (to track their evolving hosts). This creates a diverse community of bacteria and phage, maintained by kill-the-winner ecological dynamics. When the resulting communities are visualized as bipartite networks of who infects whom, they show the nested-modular structure characteristic of the Atlantic sample. The statistical significance and strength of this observation varies depending on whether the interaction networks take into account the density of the interacting strains, with implications for interpretation of interaction networks constructed by different methods. Our results suggest that the apparently complex community structures associated with marine bacteria and phage may arise from relatively simple coevolutionary origins.

The evolution of bacterial mutation rates under simultaneous selection by interspecific and social parasitism

Many bacterial populations harbour substantial numbers of hypermutable bacteria, in spite of hypermutation being associated with deleterious mutations. One reason for the persistence of hypermutators is the provision of novel mutations, enabling rapid adaptation to continually changing environments, for example coevolving virulent parasites. However, hypermutation also increases the rate at which intraspecific parasites (social cheats) are generated. Interspecific and intraspecific parasitism are therefore likely to impose conflicting selection pressure on mutation rate. Here, we combine theory and experiments to investigate how simultaneous selection from inter- and intraspecific parasitism affects the evolution of bacterial mutation rates in the plant-colonizing bacterium Pseudomonas fluorescens. Both our theoretical and experimental results suggest that phage presence increases and selection for public goods cooperation (the production of iron-scavenging siderophores) decreases selection for mutator bacteria. Moreover, phages imposed a much greater growth cost than social cheating, and when both selection pressures were imposed simultaneously, selection for cooperation did not affect mutation rate evolution. Given the ubiquity of infectious phages in the natural environment and clinical infections, our results suggest that phages are likely to be more important than social interactions in determining mutation rate evolution.

Lotka-Volterra dynamics kills the Red Queen: population size fluctuations and associated stochasticity dramatically change host-parasite coevolution

BACKGROUND

Host-parasite coevolution is generally believed to follow Red Queen dynamics consisting of ongoing oscillations in the frequencies of interacting host and parasite alleles. This belief is founded on previous theoretical work, which assumes infinite or constant population size. To what extent are such sustained oscillations realistic?

RESULTS

Here, we use a related mathematical modeling approach to demonstrate that ongoing Red Queen dynamics is unlikely. In fact, they collapse rapidly when two critical pieces of realism are acknowledged: (i) population size fluctuations, caused by the antagonism of the interaction in concordance with the Lotka-Volterra relationship; and (ii) stochasticity, acting in any finite population. Together, these two factors cause fast allele fixation. Fixation is not restricted to common alleles, as expected from drift, but also seen for originally rare alleles under a wide parameter space, potentially facilitating spread of novel variants.

CONCLUSION

Our results call for a paradigm shift in our understanding of host-parasite coevolution, strongly suggesting that these are driven by recurrent selective sweeps rather than continuous allele oscillations.