Phage-bacterial evolutionary interactions: experimental models and complications

Phage treatment of bacterial infections has offered some hope even as the crisis of antimicrobial resistance continues to be on the rise. However, bacterial resistance to phage is another looming challenge capable of undermining the effectiveness of phage therapy. Moreover, the consideration of including phage therapy in modern medicine calls for more careful research around every aspect of phage study. In an attempt to adequately prepare for the events of phage resistance, many studies have attempted to experimentally evolve phage resistance in different bacterial strains, as well as train phages to evolve counter-infectivity of resistant bacterial mutants, in view of answering such questions as coevolutionary dynamics between phage and bacteria, mechanisms of phage resistance, fitness costs of phage resistance on bacteria, etc. In this review, we summarised many such studies and by careful examination, highlighted critical issues to the outcome of phage therapy. We also discuss the insufficiency of many of these in vitro studies to represent actual disease conditions during phage application, alongside other complications that exist in phage-bacterial evolutionary interactions. Conclusively, we present the exploitation of phage-bacterial interactions for successful infection managements, as well as some future perspectives to direct phage research.

Phenotypic flux: The role of physiology in explaining the conundrum of bacterial persistence amid phage attack

Bacteriophages, the viruses of bacteria, have been studied for over a century. They were not only instrumental in laying the foundations of molecular biology, but they are also likely to play crucial roles in shaping our biosphere and may offer a solution to the control of drug-resistant bacterial infections. However, it remains challenging to predict the conditions for bacterial eradication by phage predation, sometimes even under well-defined laboratory conditions, and, most curiously, if the majority of surviving cells are genetically phage-susceptible. Here, I propose that even clonal phage and bacterial populations are generally in a state of continuous 'phenotypic flux', which is caused by transient and nongenetic variation in phage and bacterial physiology. Phenotypic flux can shape phage infection dynamics by reducing the force of infection to an extent that allows for coexistence between phages and susceptible bacteria. Understanding the mechanisms and impact of phenotypic flux may be key to providing a complete picture of phage-bacteria coexistence. I review the empirical evidence for phenotypic variation in phage and bacterial physiology together with the ways they have been modeled and discuss the potential implications of phenotypic flux for ecological and evolutionary dynamics between phages and bacteria, as well as for phage therapy.

Resistance of Dickeya solani strain IPO 2222 to lytic bacteriophage ΦD5 results in fitness tradeoffs for the bacterium during infection

Resistance to bacteriophage infections protects bacteria in phage-replete environments, enabling them to survive and multiply in the presence of their viral predators. However, such resistance may confer costs for strains, reducing their ecological fitness as expressed as competitiveness for resources or virulence or both. There is limited knowledge about such costs paid by phage-resistant plant pathogenic bacteria in their natural habitats. This study analyzed the costs of phage resistance paid by the phytopathogenic pectinolytic bacterium Dickeya solani both in vitro and in potato (Solanum tuberosum L.) plants. Thirteen Tn5 mutants of D. solani IPO 2222 were identified that exhibited resistance to infection by lytic bacteriophage vB_Dsol_D5 (ΦD5). The genes disrupted in these mutants encoded proteins involved in the synthesis of bacterial envelope components (viz. LPS, EPS and capsule). Although phage resistance did not affect most of the phenotypes of ΦD5-resistant D. solani such as growth rate, production of effectors, swimming and swarming motility, use of various carbon and nitrogen sources and biofilm formation evaluated in vitro, all phage resistant mutants were significantly compromised in their ability to survive on leaf surfaces as well as to grow within and cause disease symptoms in potato plants.

Further Considerations on How to Improve Phage Therapy Experimentation, Practice, and Reporting: Pharmacodynamics Perspectives

Phage therapy uses bacterial viruses (bacteriophages) to infect and kill targeted pathogens. Approximately one decade ago, I started publishing on how possibly to improve upon phage therapy experimentation, practice, and reporting. Here, I gather and expand upon some of those suggestions. The issues emphasized are (1) that using ratios of antibacterial agents to bacteria is not how dosing is accomplished in the real world, (2) that it can be helpful to not ignore Poisson distributions as a means of either anticipating or characterizing phage therapy success, and (3) how to calculate a concept of 'inundative phage densities.' Together, these are issues of phage therapy pharmacodynamics, meaning they are ways of thinking about the potential for phage therapy treatments to be efficacious mostly independent of the details of delivery of phages to targeted bacteria. Much emphasis is placed on working with Poisson distributions to better align phage therapy with other antimicrobial treatments.

Gauging genetic diversity of generalists: A test of genetic and ecological generalism with RNA virus experimental evolution

Generalist viruses, those with a comparatively larger host range, are considered more likely to emerge on new hosts. The potential to emerge in new hosts has been linked to viral genetic diversity, a measure of evolvability. However, there is no consensus on whether infecting a larger number of hosts leads to higher genetic diversity, or whether diversity is better maintained in a homogeneous environment, similar to the lifestyle of a specialist virus. Using experimental evolution with the RNA bacteriophage phi6, we directly tested whether genetic generalism (carrying an expanded host range mutation) or environmental generalism (growing on heterogeneous hosts) leads to viral populations with more genetic variation. Sixteen evolved viral lineages were deep sequenced to provide genetic evidence for population diversity. When evolved on a single host, specialist and generalist genotypes both maintained the same level of diversity (measured by the number of single nucleotide polymorphisms (SNPs) above 1%, P = 0.81). However, the generalist genotype evolved on a single host had higher SNP levels than generalist lineages under two heterogeneous host passaging schemes (P = 0.001, P < 0.001). RNA viruses’ response to selection in alternating hosts reduces standing genetic diversity compared to those evolving in a single host to which the virus is already well-adapted.

Antagonistic pleiotropy for carbon use is rare in new mutations

Pleiotropic effects of mutations underlie diverse biological phenomena such as ageing and specialization. In particular, antagonistic pleiotropy ("AP": when a mutation has opposite fitness effects in different environments) generates tradeoffs, which may constrain adaptation. Models of adaptation typically assume that AP is common - especially among large-effect mutations - and that pleiotropic effect sizes are positively correlated. Empirical tests of these assumptions have focused on de novo beneficial mutations arising under strong selection. However, most mutations are actually deleterious or neutral, and may contribute to standing genetic variation that can subsequently drive adaptation. We quantified the incidence, nature, and effect size of pleiotropy for carbon utilization across 80 single mutations in Escherichia coli that arose under mutation accumulation (i.e., weak selection). Although ∼46% of the mutations were pleiotropic, only 11% showed AP; among beneficial mutations, only ∼4% showed AP. In some environments, AP was more common in large-effect mutations; and AP effect sizes across environments were often negatively correlated. Thus, AP for carbon use is generally rare (especially among beneficial mutations); is not consistently enriched in large-effect mutations; and often involves weakly deleterious antagonistic effects. Our unbiased quantification of mutational effects therefore suggests that antagonistic pleiotropy may be unlikely to cause maladaptive tradeoffs.

Leaky resistance and the conditions for the existence of lytic bacteriophage

In experimental cultures, when bacteria are mixed with lytic (virulent) bacteriophage, bacterial cells resistant to the phage commonly emerge and become the dominant population of bacteria. Following the ascent of resistant mutants, the densities of bacteria in these simple communities become limited by resources rather than the phage. Despite the evolution of resistant hosts, upon which the phage cannot replicate, the lytic phage population is most commonly maintained in an apparently stable state with the resistant bacteria. Several mechanisms have been put forward to account for this result. Here we report the results of population dynamic/evolution experiments with a virulent mutant of phage Lambda, λ^VIR, and Escherichia coli in serial transfer cultures. We show that, following the ascent of λ^VIR-resistant bacteria, λ^VIR is maintained in the majority of cases in maltose-limited minimal media and in all cases in nutrient-rich broth. Using mathematical models and experiments, we show that the dominant mechanism responsible for maintenance of λ^VIR in these resource-limited populations dominated by resistant E. coli is a high rate of either phenotypic or genetic transition from resistance to susceptibility—a hitherto undemonstrated mechanism we term "leaky resistance." We discuss the implications of leaky resistance to our understanding of the conditions for the maintenance of phage in populations of bacteria—their “existence conditions.”

While it is clear that bacteriophage abound in bacterial communities, their role in the ecology and evolution of these communities remains poorly understood. Fundamental questions remain unanswered, such as, are phage regulating the population densities of their host bacteria? And how are virulent phage maintained in bacterial communities, following the seemingly inevitable evolution of resistant bacteria? Here we present a theoretical and experimental investigation to provide evidence for a new mechanism for maintaining phage in populations dominated by resistant bacteria. This mechanism, which we term “leaky resistance,” is based on a high rate of either phenotypic or genetic transition from resistance to susceptibility.

Genome Variation in the Model Halophilic Bacterium Salinibacter ruber

The halophilic bacterium Salinibacter ruber is an abundant and ecologically important member of halophilic communities worldwide. Given its broad distribution and high intraspecific genetic diversity, S. ruber is considered one of the main models for ecological and evolutionary studies of bacterial adaptation to hypersaline environments. However, current insights on the genomic diversity of this species is limited to the comparison of the genomes of two co-isolated strains. Here, we present a comparative genomic analysis of eight S. ruber strains isolated at two different time points in each of two different Mediterranean solar salterns. Our results show an open pangenome with contrasting evolutionary patterns in the core and accessory genomes. We found that the core genome is shaped by extensive homologous recombination (HR), which results in limited sequence variation within population clusters. In contrast, the accessory genome is modulated by horizontal gene transfer (HGT), with genomic islands and plasmids acting as gateways to the rest of the genome. In addition, both types of genetic exchange are modulated by restriction and modification (RM) or CRISPR-Cas systems. Finally, genes differentially impacted by such processes reveal functional processes potentially relevant for environmental interactions and adaptation to extremophilic conditions. Altogether, our results support scenarios that conciliate “Neutral” and “Constant Diversity” models of bacterial evolution.

Emerging Interaction Patterns in the Emiliania huxleyi-EhV System

Viruses are thought to be fundamental in driving microbial diversity in the oceanic planktonic realm. That role and associated emerging infection patterns remain particularly elusive for eukaryotic phytoplankton and their viruses. Here we used a vast number of strains from the model system Emiliania huxleyi/Emiliania huxleyi Virus to quantify parameters such as growth rate (µ), resistance (R), and viral production (Vp) capacities. Algal and viral abundances were monitored by flow cytometry during 72-h incubation experiments. The results pointed out higher viral production capacity in generalist EhV strains, and the virus-host infection network showed a strong co-evolution pattern between E. huxleyi and EhV populations. The existence of a trade-off between resistance and growth capacities was not confirmed.

Adaptation of the pathogen, Pseudomonas syringae, during experimental evolution on a native vs. alternative host plant

The specialization and distribution of pathogens among species has substantial impact on disease spread, especially when reservoir hosts can maintain high pathogen densities or select for increased pathogen virulence. Theory predicts that optimal within‐host growth rate will vary among host genotypes/species and therefore that pathogens infecting multiple hosts should experience different selection pressures depending on the host environment in which they are found. This should be true for pathogens with broad host ranges, but also those experiencing opportunistic infections on novel hosts or that spill over among host populations. There is very little empirical data, however, regarding how adaptation to one host might directly influence infectivity and growth on another. We took an experimental evolution approach to examine short‐term adaptation of the plant pathogen, Pseudomonas syringae pathovar tomato, to its native tomato host compared with an alternative host, Arabidopsis, in either the presence or absence of bacteriophages. After four serial passages (20 days of selection in planta), we measured bacterial growth of selected lines in leaves of either the focal or alternative host. We found that passage through Arabidopsis led to greater within‐host bacterial densities in both hosts than did passage through tomato. Whole genome resequencing of evolved isolates identified numerous single nucleotide polymorphisms based on our novel draft assembly for strain PT23. However, there was no clear pattern of clustering among plant selection lines at the genetic level despite the phenotypic differences observed. Together, the results emphasize that previous host associations can influence the within‐host growth rate of pathogens.

Evolutionary Ecology of Prokaryotic Immune Mechanisms

Bacteria have a range of distinct immune strategies that provide protection against bacteriophage (phage) infections. While much has been learned about the mechanism of action of these defense strategies, it is less clear why such diversity in defense strategies has evolved. In this review, we discuss the short- and long-term costs and benefits of the different resistance strategies and, hence, the ecological conditions that are likely to favor the different strategies alone and in combination. Finally, we discuss some of the broader consequences, beyond resistance to phage and other genetic elements, resulting from the operation of different immune strategies.

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.

Genomic and Gene-Expression Comparisons among Phage-Resistant Type-IV Pilus Mutants of Pseudomonas syringae pathovar phaseolicola

Pseudomonas syringae pv. phaseolicola (Pph) is a significant bacterial pathogen of agricultural crops, and phage Φ6 and other members of the dsRNA virus family Cystoviridae undergo lytic (virulent) infection of Pph, using the type IV pilus as the initial site of cellular attachment. Despite the popularity of Pph/phage Φ6 as a model system in evolutionary biology, Pph resistance to phage Φ6 remains poorly characterized. To investigate differences between phage Φ6 resistant Pph strains, we examined genomic and gene expression variation among three bacterial genotypes that differ in the number of type IV pili expressed per cell: ordinary (wild-type), non-piliated, and super-piliated. Genome sequencing of non-piliated and super-piliated Pph identified few mutations that separate these genotypes from wild type Pph--and none present in genes known to be directly involved in type IV pilus expression. Expression analysis revealed that 81.1% of gene ontology (GO) terms up-regulated in the non-piliated strain were down-regulated in the super-piliated strain. This differential expression is particularly prevalent in genes associated with respiration--specifically genes in the tricarboxylic acid cycle (TCA) cycle, aerobic respiration, and acetyl-CoA metabolism. The expression patterns of the TCA pathway appear to be generally up and down-regulated, in non-piliated and super-piliated Pph respectively. As pilus retraction is mediated by an ATP motor, loss of retraction ability might lead to a lower energy draw on the bacterial cell, leading to a different energy balance than wild type. The lower metabolic rate of the super-piliated strain is potentially a result of its loss of ability to retract.

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.

Convergent evolution toward an improved growth rate and a reduced resistance range in Prochlorococcus strains resistant to phage

Prochlorococcus is an abundant marine cyanobacterium that grows rapidly in the environment and contributes significantly to global primary production. This cyanobacterium coexists with many cyanophages in the oceans, likely aided by resistance to numerous co-occurring phages. Spontaneous resistance occurs frequently in Prochlorococcus and is often accompanied by a pleiotropic fitness cost manifested as either a reduced growth rate or enhanced infection by other phages. Here, we assessed the fate of a number of phage-resistant Prochlorococcus strains, focusing on those with a high fitness cost. We found that phage-resistant strains continued evolving toward an improved growth rate and a narrower resistance range, resulting in lineages with phenotypes intermediate between those of ancestral susceptible wild-type and initial resistant substrains. Changes in growth rate and resistance range often occurred in independent events, leading to a decoupling of the selection pressures acting on these phenotypes. These changes were largely the result of additional, compensatory mutations in noncore genes located in genomic islands, although genetic reversions were also observed. Additionally, a mutator strain was identified. The similarity of the evolutionary pathway followed by multiple independent resistant cultures and clones suggests they undergo a predictable evolutionary pathway. This process serves to increase both genetic diversity and infection permutations in Prochlorococcus populations, further augmenting the complexity of the interaction network between Prochlorococcus and its phages in nature. Last, our findings provide an explanation for the apparent paradox of a multitude of resistant Prochlorococcus cells in nature that are growing close to their maximal intrinsic growth rates.

Characterization of a single mutation in TraQ in a strain of Escherichia coli partially resistant to Qβ infection

Bacteria and virulent bacteriophages are in a prey–predator relationship. Experimental models under simplified conditions with the presence of bacteria and bacteriophages have been used to elucidate the mechanisms that have enabled both prey and predator to coexist over long periods. In experimental coevolution conducted with Escherichia coli and the virulent RNA bacteriophage Qβ in serial transfer, both coexisted for at least for 54 days, during which time they continued to change genetically and phenotypically. By day 16, an E. coli strain partially resistant to Qβ appeared and caused an approximately 10⁴-fold decrease in Qβ amplification. Whole-genome analysis of this strain suggested that a single mutation in TraQ was responsible for the partially resistant phenotype. TraQ interacts with propilin, encoded by the traA gene and a precursor of pilin, which is a component of the F pilus. The present study was performed to elucidate the mechanism underlying the coexistence of E. coli and Qβ by investigating how a mutation in TraQ altered the physiological state of E. coli, and thus the amplification of Qβ. Overexpression of wild-type TraQ in the partially resistant E. coli strain resulted in recovery of both TraA protein content, including propilin and pilin, and Qβ amplification to levels comparable to those observed in the susceptible strain. Intriguingly, overexpression of the mutant TraQ in the partially resistant strains also increased the levels of TraA protein and Qβ amplification, but these increases were smaller than those observed in the wild-type strain or the partially resistant strain expressing wild-type TraQ. The results of this study represent an example of how E. coli can become partially resistant to RNA bacteriophage infection via changes in a protein involved in maturation of a receptor rather than in the receptor itself and of how E. coli can stably coexist with virulent RNA bacteriophages.

Phage therapy: Should bacterial resistance to phages be a concern, even in the long run?

Bacteriophage therapy, the use of viruses that infect bacteria as antimicrobials, has been championed as a promising alternative to conventional antibiotics. Although in the laboratory bacterial resistance against phages arises rapidly, resistance so far has been an only minor problem for the effectiveness of phage therapy. Resistance to antibiotics, however, has become a major issue after decades of extensive use. Should we expect similar problems after long-term use of phages as antimicrobials? Like antibiotics, phages are often noted to be drivers of bacterial evolution. Should we expect phage-treated pathogens to develop a general resistance to phages over time, a resistance against which only, for example, hypothetically co-evolved phages might be infective? Here we argue that the global infection patterns of phages suggest that this is not necessarily a concern as environmental phages often can infect bacteria with which those phages lack any recent co-evolutionary history.

Co-evolutionary dynamics of the bacteria Vibrio sp. CV1 and phages V1G, V1P1, and V1P2: implications for phage therapy

Bacterial infections are the second largest cause of mortality in shrimp hatcheries. Among them, bacteria from the genus Vibrio constitute a major threat. As the use of antibiotics may be ineffective and banned from the food sector, alternatives are required. Historically, phage therapy, which is the use of bacteriophages, is thought to be a promising option to fight against bacterial infections. However, as for antibiotics, resistance can be rapidly developed. Since the emergence of resistance is highly undesirable, a formal characterization of the dynamics of its acquisition is mandatory. Here, we explored the co-evolutionary dynamics of resistance between the bacteria Vibrio sp. CV1 and the phages V1G, V1P1, and V1P2. Single-phage treatments as well as a cocktail composed of the three phages were considered. We found that in the presence of a single phage, bacteria rapidly evolved resistance, and the phages decreased their infectivity, suggesting that monotherapy may be an inefficient treatment to fight against Vibrio infections in shrimp hatcheries. On the contrary, the use of a phage cocktail considerably delayed the evolution of resistance and sustained phage infectivity for periods in which shrimp larvae are most susceptible to bacterial infections, suggesting the simultaneous use of multiple phages as a serious strategy for the control of vibriosis. These findings are very promising in terms of their consequences to different industrial and medical scenarios where bacterial infections are present.

Viral host-adaptation: insights from evolution experiments with phages

Phages, viral parasites of bacteria, share fundamental features of pathogenic animal and plant viruses and represent a highly tractable empirical model system to understand viral evolution and in particular viral host-adaptation. Phage adaptation to a particular host genotype often results in improved fitness by way of parallel evolution whereby independent lineages hit upon identical adaptive solutions. By contrast, phage adaptation to an evolving host population leads to the evolution of increasing host-range over time and correlated phenotypic and genetic divergence between populations. Phage host-range expansion frequently occurs by a process of stepwise evolution of multiple mutations, and host-shifts are often constrained by mutational availability, pleiotropic costs or ecological conditions.

Predation in homogeneous and heterogeneous phage environments affects virulence determinants of Pseudomonas aeruginosa

The rise of bacterial variants in the presence of lytic phages has been one of the basic grounds for evolution studies. However, there are incongruent results among different studies investigating the effect of phage resistance acquisition on bacterial fitness and virulence. We used experimental evolution to generate three classes of Pseudomonas aeruginosa variants under selective pressure from two different homogeneous phage environments and one heterogeneous phage environment. The fitness and virulence determinants of the variants, such as growth, motility, biofilm formation, resistance to oxidative stress, and the production of siderophores and chromophores, changed significantly compared to the control. Variants with similar colony morphology that were developed through different phage treatments have different phenotypic traits. Also, mRNA transcription for genes associated with certain phenotypic traits changed significantly; however, sequencing did not reveal any point mutations in selected gene loci. Furthermore, the appearance of small colony variants and melanogenic variants and the increase in pyocyanin and pyoverdin production for some variants are believed to affect the virulence of the population. The knowledge gained from this study will fundamentally contribute to our understanding of the evolutionary dynamics of bacteria under phage selective pressure which is crucial to the efficient utilization of bacteriophages in medical contexts.

What Can Phages Tell Us about Host-Pathogen Coevolution?

The outcomes of host-parasite interactions depend on the coevolutionary forces acting upon them, but because every host-parasite relation is enmeshed in a web of biotic and abiotic interactions across a heterogeneous landscape, host-parasite coevolution has proven difficult to study. Simple laboratory phage-bacteria microcosms can ameliorate this difficulty by allowing controlled, well-replicated experiments with a limited number of interactors. Genetic, population, and life history data obtained from these studies permit a closer examination of the fundamental correlates of host-parasite coevolution. In this paper, I describe the results of phage-bacteria coevolutionary studies and their implications for the study of host-parasite coevolution. Recent experimental studies have confirmed phage-host coevolutionary dynamics in the laboratory and have shown that coevolution can increase parasite virulence, specialization, adaptation, and diversity. Genetically, coevolution frequently proceeds in a manner best described by the Gene for Gene model, typified by arms race dynamics, but certain contexts can result in Red Queen dynamics according to the Matching Alleles model. Although some features appear to apply only to phage-bacteria systems, other results are broadly generalizable and apply to all instances of antagonistic coevolution. With laboratory host-parasite coevolutionary studies, we can better understand the perplexing array of interactions that characterize organismal diversity in the wild.

Bacteria-virus coevolution

Phages, viruses of bacteria, are ubiquitous. Many phages require host cell death to successfully complete their life cycle, resulting in reciprocal evolution of bacterial resistance and phage infectivity (antagonistic coevolution). Such coevolution can have profound consequences at all levels of biological organisation. Here, we review genetic and ecological factors that contribute to determining coevolutionary dynamics between bacteria and phages. We also consider some of the consequences of bacteria-phage coevolution, such as determining rates of molecular evolution and structuring communities, and how these in turn feedback into driving coevolutionary dynamics.

Genomic island variability facilitates Prochlorococcus-virus coexistence

Prochlorococcus cyanobacteria are extremely abundant in the oceans, as are the viruses that infect them. How hosts and viruses coexist in nature remains unclear, although the presence of both susceptible and resistant cells may allow this coexistence. Combined whole-genome sequencing and PCR screening technology now enables us to investigate the effect of resistance on genome evolution and the genomic mechanisms behind the long-term coexistence of Prochlorococcus and their viruses. Here we present a genome analysis of 77 substrains selected for resistance to ten viruses, revealing mutations primarily in non-conserved, horizontally transferred genes that localize to a single hypervariable genomic island. Mutations affected viral attachment to the cell surface and imposed a fitness cost to the host, manifested by significantly lower growth rates or a previously unknown mechanism of more rapid infection by other viruses. The mutant genes are generally uncommon in nature yet some carry polymorphisms matching those found experimentally. These data are empirical evidence indicating that viral-attachment genes are preferentially located in genomic islands and that viruses are a selective pressure enhancing the diversity of both island genes and island gene content. This diversity emerges as a genomic mechanism that reduces the effective host population size for infection by a given virus, thus facilitating long-term coexistence between viruses and their hosts in nature.

Using experimental evolution to explore natural patterns between bacterial motility and resistance to bacteriophages

Resistance of bacteria to phages may be gained by alteration of surface proteins to which phages bind, a mechanism that is likely to be costly as these molecules typically have critical functions such as movement or nutrient uptake. To address this potential trade-off, we combine a systematic study of natural bacteria and phage populations with an experimental evolution approach. We compare motility, growth rate and susceptibility to local phages for 80 bacteria isolated from horse chestnut leaves and, contrary to expectation, find no negative association between resistance to phages and bacterial motility or growth rate. However, because correlational patterns (and their absence) are open to numerous interpretations, we test for any causal association between resistance to phages and bacterial motility using experimental evolution of a subset of bacteria in both the presence and absence of naturally associated phages. Again, we find no clear link between the acquisition of resistance and bacterial motility, suggesting that for these natural bacterial populations, phage-mediated selection is unlikely to shape bacterial motility, a key fitness trait for many bacteria in the phyllosphere. The agreement between the observed natural pattern and the experimental evolution results presented here demonstrates the power of this combined approach for testing evolutionary trade-offs.

Parallel changes in host resistance to viral infection during 45,000 generations of relaxed selection

The dynamics of host susceptibility to parasites are often influenced by trade-offs between the costs and benefits of resistance. We assayed changes in the resistance to three viruses in six lines of Escherichia coli that had been evolving for almost 45,000 generations in their absence. The common ancestor of these lines was completely resistant to T6, partially resistant to T6* (a mutant of T6 with altered host range), and sensitive to λ. None of the populations changed with respect to resistance to T6, whereas all six evolved increased susceptibility to T6*, probably ameliorating a cost of resistance. More surprisingly, however, the majority of lines evolved complete resistance to λ, despite not encountering that virus during this period. By coupling our results with previous work, we infer that resistance to λ evolved as a pleiotropic effect of a beneficial mutation that downregulated an unused metabolic pathway. The strong parallelism between the lines implies that selection had almost deterministic effects on the evolution of these patterns of host resistance. The opposite outcomes for resistance to T6* and λ demonstrate that the evolution of host resistance under relaxed selection cannot be fully predicted by simple trade-off models.

An experimental study of the population and evolutionary dynamics of Vibrio cholerae O1 and the bacteriophage JSF4

Studies of Vibrio cholerae in the environment and infected patients suggest that the waning of cholera outbreaks is associated with rise in the density of lytic bacteriophage. In accordance with mathematical models, there are seemingly realistic conditions where phage predation could be responsible for declines in the incidence of cholera. Here, we present the results of experiments with the El Tor strain of V. cholerae (N16961) and a naturally occurring lytic phage (JSF4), exploring the validity of the main premise of this model: that phage predation limits the density of V. cholerae populations. At one level, the results of our experiments are inconsistent with this hypothesis. JSF4-resistant V. cholerae evolve within a short time following their confrontation with these viruses and their populations become limited by resources rather than phage predation. At a larger scale, however, the results of our experiments are not inconsistent with the hypothesis that bacteriophage modulate outbreaks of cholera. We postulate that the resistant bacteria that evolved play an insignificant role in the ecology or pathogenicity of V. cholerae. Relative to the phage-sensitive cells from whence they are derived, the evolved JSF4-resistant V. cholerae have fitness costs and other characters that are likely to impair their ability to compete with the sensitive cells in their natural habitat and may be avirulent in human hosts. The results of this in vitro study make predictions that can be tested in natural populations of V. cholerae and cholera-infected patients.

Adaptation rates of lytic viruses depend critically on whether host cells survive the bottleneck

We use a branching process approach to estimate the substitution rate, the rate at which beneficial mutations occur and fix, in populations of lytic viruses whose growth is controlled by periodic population bottlenecks. Our model predicts that substitution rates, and by extension adaptation rates, are profoundly affected by the survival of infected host cells at the bottleneck. In particular, we find that direct transfer (or environmental) bottlenecks, in which some fraction of both free virus and host cells are preserved, are associated with relatively slow adaptation rates for the virus. In contrast, viruses can adapt much more quickly when only free virus is transferred to a new host population, as is typical in an epidemiological setting. Finally, when premature lysis of the host-cell population is induced at the bottleneck, we predict that adaptation rates for the virus will, in general, be faster still. These results hold irrespective of the life-history trait affected by the beneficial mutation. The substitution rates in the presence of environmental bottlenecks are predicted to be as much as an order of magnitude lower than in the other two cases.

Mean-square stability of a stochastic model for bacteriophage infection with time delays

We consider the stability properties of the positive equilibrium of a stochastic model for bacteriophage infection with discrete time delay. Conditions for mean-square stability of the trivial solution of the linearized system around the equilibrium are given by the construction of suitable Lyapunov functionals. The numerical simulations of the strong solutions of the arising stochastic delay differential system suggest that, even for the original non-linear model, the longer the incubation time the more the phage and bacteria populations can coexist on a stable equilibrium in a noisy environment for very long time.

Host density impacts relative fitness of bacteriophage Phi6 genotypes in structured habitats

Spatially structured environments may impact evolution by restricting population sizes, limiting opportunities for genetic mixis, or weakening selection against deleterious genotypes. When habitat structure impedes dispersal, low-productivity (less virulent) infectious parasites may benefit from their prudent exploitation of local hosts. Here we explored the combined ability for habitat structure and host density to dictate the relative reproductive success of differentially productive parasites. To do so, we allowed two RNA bacteriophage Phi6 genotypes to compete in structured and unstructured (semi-solid versus liquid) habitats while manipulating the density of Pseudomonas hosts. In the unstructured habitats, the more-productive phage strain experienced a relatively constant fitness advantage regardless of starting host density. By contrast, in structured habitats, restricted phage dispersal may have magnified the importance of local productivity, thus allowing the relative fitness of the less-productive virus to improve as host density increased. Further data suggested that latent period (duration of cellular infection) and especially burst size (viral progeny produced per cell) were the phage "life-history" traits most responsible for our results. We discuss the relevance of our findings for selection occurring in natural phage populations and for the general evolutionary epidemiology of infectious parasites.

Is there a cost of virus resistance in marine cyanobacteria?

Owing to their abundance and diversity, it is generally perceived that viruses are important for structuring microbial communities and regulating biogeochemical cycles. The ecological impact of viruses on microbial food webs, however, may be influenced by evolutionary processes, including the ability of bacteria to evolve resistance to viruses and the theoretical prediction that this resistance should be accompanied by a fitness cost. We conducted experiments using phylogenetically distinct strains of marine Synechococcus (Cyanobacteria) to test for a cost of resistance (COR) to viral isolates collected from Mount Hope Bay, Rhode Island. In addition, we examined whether fitness costs (1) increased proportionally with 'total resistance', the number of viruses for which a strain had evolved resistance, or (2) were determined more by 'compositional resistance', the identity of the viruses to which it evolved resistance. A COR was only found in half of our experiments, which may be attributed to compensatory mutations or the inability to detect a small COR. When detected, the COR resulted in a approximately 20% reduction in relative fitness compared to ancestral strains. The COR was unaffected by total resistance, suggesting a pleiotropic fitness response. Under competitive conditions, however, the COR was dependent on compositional resistance, suggesting that fitness costs were associated with the identity of a few particular viruses. Our study provides the first evidence for a COR in marine bacteria, and suggests that Synechococcus production may be influenced by the composition of co-occurring viruses.

Virus population extinction via ecological traps

Populations are at risk of extinction when unsuitable or when sink habitat exceeds a threshold frequency in the environment. Sinks that present cues associated with high-quality habitats, termed ecological traps, have especially detrimental effects on net population growth at metapopulation scales. Ecological traps for viruses arise naturally, or can be engineered, via the expression of viral-binding sites on cells that preclude viral reproduction. We present a model for virus population growth in a heterogeneous host community, parameterized with data from populations of the RNA bacteriophage Phi6 presented with mixtures of suitable host bacteria and either neutral or trap cells. We demonstrate that viruses can sustain high rates of population growth in the presence of neutral non-hosts as long as some host cells are present, whereas trap cells dramatically reduce viral fitness. In addition, we demonstrate that the efficacy of traps for viral elimination is frequency dependent in spatially structured environments such that population viability is a nonlinear function of habitat loss in dispersal-limited virus populations. We conclude that the ecological concepts applied to species conservation in altered landscapes can also contribute to the development of trap cell therapies for infectious human viruses.

Viral ecology and the maintenance of novel host use

Viruses can occasionally emerge by infecting new host species. However, the early phases of emergence can hinge upon ecological sustainability of the virus population, which is a product of both within-host population growth and between-host transmission. Insufficient growth or transmission can force virus extinction before the latter phases of emergence, where genetic adaptations that improve host use may occur. We examined the early phase of emergence by studying the population dynamics of RNA phages in replicated laboratory environments containing native and novel host bacteria. To predict the breadth of transmission rates allowing viral persistence on each species, we developed a simple model based on in vitro data for phage growth rate over a range of initial population densities on both hosts. Validation of these predictions using serial passage experiments revealed a range of transmission rates for which the native host was a source and the novel host was a sink. In this critical range of transmission rates, periodic exposure to the native host was sufficient for the maintenance of the viral population on the novel host. We argue that this effect should facilitate adaptation by the virus to utilize the novel host--often crucial in subsequent phases of emergence.

Experimental coevolution with bacteria and phage. The Pseudomonas fluorescens--Phi2 model system

Parasites are ubiquitous in biological systems and antagonistic coevolution between hosts and parasites is thought be a major ecological and evolutionary force. Recent experiments using laboratory populations of bacteria and their parasitic viruses, phage, have provided the first direct empirical evidence of antagonistic coevolution in action. In this article we describe this model system and synthesise recent findings that address the causes and consequences of antagonistic coevolution.

Reconciling empirical ecology with neutral community models

Neutral community models embody the idea that individuals are ecologically equivalent, having equal fitness over all environmental conditions, and describe how the spatial dynamics and speciation of such individuals can produce a wide range of patterns of distribution, diversity, and abundance. Neutral models have been controversial, provoking a rush of tests and comments. The debate has been spurred by the suggestion that we should test mechanisms. However, the mechanisms and the spatial scales of interest have never clearly been described, and consequently, the tests have often been only peripherally relevant. At least two mechanisms are present in spatially structured neutral models. Dispersal limitation causes clumping of a species, which increases the strength of intraspecific competition and reduces the strength of interspecific competition. This may prolong coexistence and enhance local and regional diversity. Speciation is present in some neutral models and gives a donor-controlled input of new species, many of which remain rare or are short lived, but which directly add to species diversity. Spatial scale is an important consideration in neutral models. Ecological equivalence and equal fitness have implicit spatial scales because dispersal limitation and its emergent effects operate at population levels, and populations and communities are defined at a chosen spatial scale in recent neutral models; equality is measured relative to a metacommunity, and this necessitates defining the spatial scale of that metacommunity. Furthermore, dispersal has its own scales. Thorough empirical tests of neutral models will require both tests of mechanisms and pattern-producing ability, and will involve coupling theoretical models and experiments.

Spatial heterogeneity and the stability of host-parasite coexistence

Spatially heterogeneous environments can theoretically promote more stable coexistence of hosts and parasites by reducing the risk of parasite attack either through providing permanent spatial refuges or through providing ephemeral refuges by reducing dispersal. In experimental populations of Pseudomonas aeruginosa and the bacteriophage PP7, spatial heterogeneity promoted stable coexistence of host and parasite, while coexistence was significantly less stable in the homogeneous environment. Phage populations were found to be persisting on subpopulations of sensitive bacteria. Transferring populations to fresh microcosms every 24 h prevented the development of permanent spatial refuges. However, the lower dispersal rates in the heterogeneous environment were found to reduce parasite transmission thereby creating ephemeral refuges from phage attack. These results suggest that spatial heterogeneity can stabilize an otherwise unstable host-parasite interaction even in the absence of permanent spatial refuges.

The effect of a bacteriophage on diversification of the opportunistic bacterial pathogen, Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen that colonizes the lungs of cystic fibrosis (CF) patients. CF lungs often contain a diverse range of P. aeruginosa phenotypes, some of which are likely to contribute to the persistence of infection, yet the causes of diversity are unclear. While the ecological heterogeneity of the lung environment and therapeutic regimes are probable factors, a role for parasitic bacteriophage cannot be ruled out. Parasites have been implicated as a key ecological variable driving the evolution of diversity in host populations. PP7 drove cycles of morphological diversification in host populations of P. aeruginosa due to the de novo evolution of small-rough colony variants that coexisted with large diffuse colony morph bacteria. In the absence of phage, bacteria only displayed the large diffuse colony morphology of the wild-type. Further assays revealed there to be two distinct types of resistant bacteria; these had very different ecological phenotypes, yet each carried a cost of resistance.