Tradeoffs between phage resistance and nitrogen fixation drive the evolution of genes essential for cyanobacterial heterocyst functionality.

Harmful blooms caused by diazotrophic (nitrogen-fixing) Cyanobacteria are becoming increasingly frequent and negatively impact aquatic environments worldwide. Cyanophages (viruses infecting Cyanobacteria) can potentially regulate cyanobacterial blooms, yet Cyanobacteria can rapidly acquire mutations that provide protection against phage infection. Here, we provide novel insights into cyanophage:Cyanobacteria interactions by characterizing the resistance to phages in two species of diazotrophic Cyanobacteria: Nostoc sp. and Cylindrospermopsis raciborskii. Our results demonstrate that phage resistance is associated with a fitness tradeoff by which resistant Cyanobacteria have reduced ability to fix nitrogen and/or to survive nitrogen starvation. Furthermore, we use whole-genome sequence analysis of 58 Nostoc-resistant strains to identify several mutations associated with phage resistance, including in cell surface-related genes and regulatory genes involved in the development and function of heterocysts (cells specialized in nitrogen fixation). Finally, we employ phylogenetic analyses to show that most of these resistance genes are accessory genes whose evolution is impacted by lateral gene transfer events. Together, these results further our understanding of the interplay between diazotrophic Cyanobacteria and their phages and suggest that a tradeoff between phage resistance and nitrogen fixation affects the evolution of cell surface-related genes and of genes involved in heterocyst differentiation and nitrogen fixation.

Sequencing the genomes of LPP-1, the first isolated cyanophage, and its relative LPP-2 reveal different integration mechanisms in closely related phages.

In the early 1960s, the first cyanophage was isolated. The description of this phage, named LPP-1, led to the extensive investigation of various cyanophages and to the study of their interactions with their cyanobacterial hosts towards controlling blooms. Here, the genomes of LPP-1 and its putative relative, LPP-2 were sequenced. Sequencing these genomes revealed that LPP-1 and LPP-2 are members of a group of short-tailed cyanophages, which are distantly related to the T7-like cyanophages. Most of the phages in this group have the ability to lysogenize their hosts. Their ability to switch between lytic and lysogenic infection may explain the formation of cyanobacterial blooms despite the persistence of their phages. This lysogenic capacity of the LPP-1-like phages occurs despite the lack of an obvious integrase gene within their genomes. Interestingly, we show that LPP-2 integrates into the host genome through an integration site in high proximity to a recombination endonuclease that may have integrase activity. Further understanding of cyanobacterial-phage relationships may provide insight into their population dynamics and suggest novel approaches for control of destructive cyanobacterial blooms.

Isolation and characterization of a novel Lambda-like phage infecting the bloom-forming cyanobacteria Cylindrospermopsis raciborskii.

Cylindrospermopsis raciborskii is a central bloom-forming cyanobacteria. However, despite its ecological significance, little is known of its interactions with the phages that infect it. Currently, only a single sequenced genome of a Cylindrospermopsis-infecting phage is publicly available. Here we describe the isolation and characterization of Cr-LKS3, a second phage infecting Cylindrospermopsis. Cr-LKS3 is a siphovirus with a higher genome similarity to prophages within heterotrophic bacteria genomes than to any other cyanophage/cyano-prophage, suggesting that it represents a novel cyanophage group. The function, order and orientation of the 72 genes in the Cr-LKS3 genome are highly similar to those of Escherichia virus Lambda (hereafter Lambda), despite the very low sequence similarity between these phages, showing high evolutionary convergence despite the substantial difference in host characteristics. Similarly to Lambda, the genome of Cr-LKS3 contains various genes that are known to be central to lysogeny, suggesting it can enter a lysogenic cycle. Cr-LKS3 has a unique ability to infect a host with a dramatically different GC content, without carrying any tRNA genes to compensate for this difference. This ability, together with its potential lysogenic lifestyle shed light on the complex interactions between C. raciborskii and its phages.

Coevolutionary phage training leads to greater bacterial suppression and delays the evolution of phage resistance.

The evolution of antibiotic-resistant bacteria threatens to become the leading cause of worldwide mortality. This crisis has renewed interest in the practice of phage therapy. Yet, bacteria's capacity to evolve resistance may debilitate this therapy as well. To combat the evolution of phage resistance and improve treatment outcomes, many suggest leveraging phages' ability to counter resistance by evolving phages on target hosts before using them in therapy (phage training). We found that in vitro, λtrn, a phage trained for 28 d, suppressed bacteria ∼1,000-fold for three to eight times longer than its untrained ancestor. Prolonged suppression was due to a delay in the evolution of resistance caused by several factors. Mutations that confer resistance to λtrn are ∼100× less common, and while the target bacterium can evolve complete resistance to the untrained phage in a single step, multiple mutations are required to evolve complete resistance to λtrn. Mutations that confer resistance to λtrn are more costly than mutations for untrained phage resistance. Furthermore, when resistance does evolve, λtrn is better able to suppress these forms of resistance. One way that λtrn improved was through recombination with a gene in a defunct prophage in the host genome, which doubled phage fitness. This transfer of information from the host genome is an unexpected but highly efficient mode of training phage. Lastly, we found that many other independently trained λ phages were able to suppress bacterial populations, supporting the important role training could play during phage therapeutic development.

Dynamics of Adaptation During Three Years of Evolution Under Long-Term Stationary Phase.

Many bacterial species that cannot sporulate, such as the model bacterium Escherichia coli, can nevertheless survive for years, following exhaustion of external resources, in a state termed long-term stationary phase (LTSP). Here we describe the dynamics of E. coli adaptation during the first three years spent under LTSP. We show that during this time, E. coli continuously adapts genetically through the accumulation of mutations. For nonmutator clones, the majority of mutations accumulated appear to be adaptive under LTSP, reflected in an extremely convergent pattern of mutation accumulation. Despite the rapid and convergent manner in which populations adapt under LTSP, they continue to harbor extensive genetic variation. The dynamics of evolution of mutation rates under LTSP are particularly interesting. The emergence of mutators affects overall mutation accumulation rates as well as the mutational spectra and the ultimate spectrum of adaptive alleles acquired under LTSP. With time, mutators can evolve even higher mutation rates through the acquisition of additional mutation rate-enhancing mutations. Different mutator and nonmutator clones within a single population and time point can display extreme variation in their mutation rates, resulting in differences in both the dynamics of adaptation and their associated deleterious burdens. Despite these differences, clones that vary greatly in their mutation rates tend to coexist within their populations for many years, under LTSP.

Culture Volume Influences the Dynamics of Adaptation under Long-Term Stationary Phase.

Escherichia coli and many other bacterial species, which are incapable of sporulation, can nevertheless survive within resource exhausted media by entering a state termed long-term stationary phase (LTSP). We have previously shown that E. coli populations adapt genetically under LTSP in an extremely convergent manner. Here, we examine how the dynamics of LTSP genetic adaptation are influenced by varying a single parameter of the experiment-culture volume. We find that culture volume affects survival under LTSP, with viable counts decreasing as volumes increase. Across all volumes, mutations accumulate with time, and the majority of mutations accumulated demonstrate signals of being adaptive. However, positive selection appears to affect mutation accumulation more strongly at higher, compared with lower volumes. Finally, we find that several similar genes are likely involved in adaptation across volumes. However, the specific mutations within these genes that contribute to adaptation can vary in a consistent manner. Combined, our results demonstrate how varying a single parameter of an evolutionary experiment can substantially influence the dynamics of observed adaptation.

Adaptations Accumulated under Prolonged Resource Exhaustion Are Highly Transient.

Many nonsporulating bacterial species can survive for years within exhausted growth media in a state termed long-term stationary phase (LTSP). We have been carrying out evolutionary experiments aimed at elucidating the dynamics of genetic adaptation under LTSP. We showed that Escherichia coli adapts to prolonged resource exhaustion through the highly convergent acquisition of mutations. In the most striking example of such convergent adaptation, we observed that across all independently evolving LTSP populations, over 90% of E. coli cells carry mutations to one of three specific sites of the RNA polymerase core enzyme (RNAPC). These LTSP adaptations reduce the ability of the cells carrying them to grow once fresh resources are again provided. Here, we examine how LTSP populations recover from costs associated with their adaptation once resources are again provided to them. We demonstrate that due to the ability of LTSP populations to maintain high levels of standing genetic variation during adaptation, costly adaptations are very rapidly purged from the population once they are provided with fresh resources. We further demonstrate that recovery from costs acquired during adaptation under LTSP occurs more rapidly than would be possible if LTSP adaptations had fixed during the time populations spent under resource exhaustion. Finally, we previously reported that under LTSP, some clones develop a mutator phenotype, greatly increasing their mutation accumulation rates. Here, we show that the mechanisms by which populations recover from costs associated with fixed adaptations may depend on mutator status.IMPORTANCE Many bacterial species can survive for decades under starvation, following the exhaustion of external growth resources. We have previously shown that bacteria genetically adapt under these conditions in a manner that reduces their ability to grow once resources again become available. Here, we study how populations that have been subject to very prolonged resource exhaustion recover from costs associated with their adaptation. We demonstrate that rapid adaptations acquired under prolonged starvation tend to be highly transient, rapidly reducing in frequency once bacteria are no longer starved. Our results shed light on the longer-term consequences of bacterial survival under prolonged starvation. More generally, these results may also be applicable to understanding longer-term consequences of rapid adaptation to additional conditions as well.

Rapid Genetic Adaptation during the First Four Months of Survival under Resource Exhaustion.

Many bacteria, including the model bacterium Escherichia coli can survive for years within spent media, following resource exhaustion. We carried out evolutionary experiments, followed by whole genome sequencing of hundreds of evolved clones to study the dynamics by which E. coli adapts during the first 4 months of survival under resource exhaustion. Our results reveal that bacteria evolving under resource exhaustion are subject to intense selection, manifesting in rapid mutation accumulation, enrichment in functional mutation categories and extremely convergent adaptation. In the most striking example of convergent adaptation, we found that across five independent populations adaptation to conditions of resource exhaustion occurs through mutations to the three same specific positions of the RNA polymerase core enzyme. Mutations to these three sites are strongly antagonistically pleiotropic, in that they sharply reduce exponential growth rates in fresh media. Such antagonistically pleiotropic mutations, combined with the accumulation of additional mutations, severely reduce the ability of bacteria surviving under resource exhaustion to grow exponentially in fresh media. We further demonstrate that the three positions at which these resource exhaustion mutations occur are conserved for the ancestral E. coli allele, across bacterial phyla, with the exception of nonculturable bacteria that carry the resource exhaustion allele at one of these positions, at very high frequencies. Finally, our results demonstrate that adaptation to resource exhaustion is not limited by mutational input and that bacteria are able to rapidly adapt under resource exhaustion in a temporally precise manner through allele frequency fluctuations.

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.

Diversity and evolutionary relationships of T7-like podoviruses infecting marine cyanobacteria.

Phages are extremely abundant in the oceans, influencing the population dynamics, diversity and evolution of their hosts. Here we assessed the diversity and phylogenetic relationships among T7-like cyanophages using DNA polymerase (replication), major capsid (structural) and photosynthesis psbA (host-derived) genes from isolated phages. DNA polymerase and major capsid phylogeny divided them into two discrete clades with no evidence for gene exchange between clades. Clade A phages primarily infect Synechococcus while clade B phages infect either Synechococcus or Prochlorococcus. The major capsid gene of one of the phages from clade B carries a putative intron. Nearly all clade B phages encode psbA whereas clade A phages do not. This suggests an ancient separation between cyanophages from these two clades, with the acquisition or loss of psbA occurring around the time of their divergence. A mix and match of clustering patterns was found for the replication and structural genes within each major clade, even among phages infecting different host genera. This is suggestive of numerous gene exchanges within each major clade and indicates that core phage functions have not coevolved with specific hosts. In contrast, clustering of phage psbA broadly tracks that of the host genus. These findings suggest that T7-like cyanophages evolve through clade-limited gene exchanges and that different genes are subjected to vastly different selection pressures.

Virus-host swinging party in the oceans: Incorporating biological complexity into paradigms of antagonistic coexistence.

Bacteria and their viruses (phages) are antagonists, yet have coexisted in nature for billions of years. Models proposed to explain the paradox of antagonistic coexistence generally reach two types of solutions: Arms race-like dynamics that lead to hosts and viruses with increasing resistance and infection ranges; and population fluctuations between diverse host and viral types due to a metabolic cost of resistance. Recently, we found that populations of the marine cyanobacterium, Prochlorococcus, consist of cells with extreme hypervariability in gene sequence and gene content in a viral susceptibility region of the genome. Furthermore, we found a novel cost of resistance where resistance to one set of viruses is accompanied by changes in infection dynamics by other viruses. In this combined mini-review and commentary paper we discuss these findings in the context of existing ecological, evolutionary and genetic models of host-virus coexistence. We suggest that this coexistence is governed mainly by fluctuations between microbial subpopulations with differing viral susceptibility regions and that these fluctuations are driven by both metabolic and enhanced infection costs of resistance. Furthermore, we suggest that enhanced infection leads to passive host-switching by viruses, preventing the development of hosts with universal resistance. These findings highlight the vital importance of community complexity for host-virus coexistence.

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.