Transcriptome dynamics of a broad host-range cyanophage and its hosts

Expression and characterization of the complete cyanophage genome PP in the heterologous host Synechococcus elongatus PCC 7942

In this study, we successfully integrated the full-length genome of the cyanophage PP into the non-host cyanobacterium Synechococcus elongatus PCC 7942, facilitated by conjugation via Escherichia coli. To address the challenge posed by the toxic open reading frames (ORFs) of PP in E. coli, we first identified and characterized three toxic ORFs. The PP genome was subsequently rearranged, and the expression of these toxic ORFs was controlled using a tandem-induction switch system. The full-length PP genome was then successfully integrated into the genome of S. elongatus PCC 7942. Interestingly, the integration of the PP genome led to a reduction in photosynthesis and carbon fixation in S. elongatus PCC 7942, resembling the effects typically associated with cyanophage infection. Transcriptomic analysis showed that 32 of the 41 ORFs in the PP genome were actively transcribed in S. elongatus PCC 7942, significantly affecting energy metabolism and carbon fixation pathways. These effects were further confirmed by metabolomic analysis.

Infection and Genomic Properties of Single- and Double-Stranded DNA Cellulophaga Phages

Bacterial viruses (phages) are abundant and ecologically impactful, but laboratory-based experimental model systems vastly under-represent known phage diversity, particularly for ssDNA phages. Here, we characterize the genomes and infection properties of two unrelated marine flavophages-ssDNA generalist phage phi18:4 (6.5 Kbp) and dsDNA specialist phage phi18:1 (39.2 Kbp)-when infecting the same Cellulophaga baltica strain #18 (Cba18), of the class Flavobacteriia. Phage phi18:4 belongs to a new family of ssDNA phages, has an internal lipid membrane, and its genome encodes primarily structural proteins, as well as a DNA replication protein common to ssDNA phages and a unique lysis protein. Phage phi18:1 is a siphovirus that encodes several virulence genes, despite not having a known temperate lifestyle, a CAZy enzyme likely for regulatory purposes, and four DNA methyltransferases dispersed throughout the genome that suggest both host modulation and phage DNA protection against host restriction. Physiologically, ssDNA phage phi18:4 has a shorter latent period and smaller burst size than dsDNA phage phi18:1, and both phages efficiently infect this host. These results help augment the diversity of characterized environmental phage-host model systems by studying infections of genomically diverse phages (ssDNA vs. dsDNA) on the same host.

Combined metabolomic and genomic analyses reveal phage-specific and infection stage-specific alterations to marine Roseobacter metabolism

Phages can reshape the metabolic network of hosts to support specific requirements for replication during infection. However, metabolomic profiling of phage-elicited host global metabolic alterations and the linkage of phage-encoded auxiliary metabolic genes to these alterations are understudied. In this study, the dynamics of intracellular metabolites of Dinoroseobacter shibae DFL12, a member of marine environmentally and biogeochemically relevant Roseobacter clade, in response to four distinct lytic roseophage infections were investigated. Metabolomic profiling indicated that roseophage infections significantly altered host metabolism in a phage-specific manner. Pathway enrichment analyses showed that the central carbon pathway and DNA, amino acid, and coenzyme metabolism were commonly altered by roseophages, revealing a central role of these pathways in phage replication. Furthermore, clear infection stage-specific host responses were observed, corresponding to different metabolic demands of phage replication in the early and late infection stages. Interestingly, the content of host vitamin B1, which is the essential nutrient provided by D. shibae to its symbiotic microalgae, increased in the early infection stage for most roseophages, implying that phage infection may impact the symbiosis of D. shibae with microalgae. Finally, combined metabolomic and phage genomics analyses showed that roseophages adopt different strategies to expand the host pyrimidine pool (recycling or de novo synthesis of pyrimidine nucleotides), and this difference was likely related to variation in the GC content between phage and host genomes. Collectively, these results highlight the potential importance of phage-specific and infection stage-specific host metabolic reprogramming in marine phage-host interactions, bacteria-microalgae symbiosis, and biogeochemical cycles.

Adaptive loss of tRNA gene expression leads to phage resistance in a marine Synechococcus cyanobacterium

Synechococcus is a significant primary producer in the oceans, coexisting with cyanophages, which are important agents of mortality. Bacterial resistance against phage infection is a topic of significant interest, yet little is known for ecologically relevant systems. Here we use exogenous gene expression and gene disruption to investigate mechanisms underlying intracellular resistance of marine Synechococcus WH5701 to the Syn9 cyanophage. The restriction-modification and Gabija defence systems possessed by Synechococcus WH5701 did not contribute to resistance. Instead, resistance was primarily driven by insufficient levels of LeuTAA tRNA, preventing translation of key phage genes in a passive, intracellular mode of resistance. Restoring cellular tRNA expression rendered the cyanobacterium sensitive to infection. We propose an evolutionary scenario whereby changes in cell codon usage, acquisition of tRNAs by the phage and loss of cell and phage tRNA expression resulted in an effective means of resistance, highlighting the dynamic interplay between bacteria and phages in shaping their co-evolutionary trajectories.

Cyanobacteria-cyanophage interactions between freshwater and marine ecosystems based on large-scale cyanophage genomic analysis

The disparities in harmful algal blooms dynamics are largely attributed to variations in cyanobacteria populations within aquatic ecosystems. However, cyanobacteria-cyanophage interactions and their role in shaping cyanobacterial populations has been previously underappreciated. To address this knowledge gap, we isolated and sequenced 42 cyanophages from diverse water sources in China, with the majority (n = 35) originating from freshwater sources. We designated these sequences as the "Novel Cyanophage Genome sequence Collection" (NCGC). NCGC displayed notable genetic variations, with 95 % (40/42) of the sequences representing previously unidentified taxonomic ranks. By integrating NCGC with public data of cyanophages and cyanobacteria, we found evidence for more frequent historical cyanobacteria-cyanophage interactions in freshwater ecosystems. This was evidenced by a higher prevalence of prophage integrase-related genes in freshwater cyanophages (37.97 %) than marine cyanophages (7.42 %). In addition, freshwater cyanophages could infect a broader range of cyanobacteria orders (n = 4) than marine ones (n = 0). Correspondingly, freshwater cyanobacteria harbored more defense systems per million base pairs in their genomes, indicating more frequent phage infections. Evolutionary and cyanophage epidemiological studies suggest that interactions between cyanobacteria and cyanophages in freshwater and marine ecosystems are interconnected, and that brackish water can act as a transitional zone for freshwater and marine cyanophages. In conclusion, our research significantly expands the genetic information database of cyanophage, offering a wider selection of cyanophages to control harmful cyanobacterial blooms. Additionally, we represent a pioneering large-scale and comprehensive analysis of cyanobacteria and cyanophage sequencing data, and it provides theoretical guidance for the application of cyanophages in different environments.

Identification of Shemin pathway genes for tetrapyrrole biosynthesis in bacteriophage sequences from aquatic environments

Tetrapyrroles such as heme, chlorophyll, and vitamin B12 are essential for various metabolic pathways. They derive from 5-aminolevulinic acid (5-ALA), which can be synthesized by a single enzyme (5-ALA synthase or AlaS, Shemin pathway) or by a two-enzyme pathway. The genomes of some bacteriophages from aquatic environments carry various tetrapyrrole biosynthesis genes. Here, we analyze available metagenomic datasets and identify alaS homologs (viral alaS, or valaS) in sequences corresponding to marine and freshwater phages. The genes are found individually or as part of complete or truncated three-gene loci encoding heme-catabolizing enzymes. Amino-acid sequence alignments and three-dimensional structure prediction support that the valaS sequences likely encode functional enzymes. Indeed, we demonstrate that is the case for a freshwater phage valaS sequence, as it can complement an Escherichia coli 5-ALA auxotroph, and an E. coli strain overexpressing the gene converts the typical AlaS substrates glycine and succinyl-CoA into 5-ALA. Thus, our work identifies valaS as an auxiliary metabolic gene in phage sequences from aquatic environments, further supporting the importance of tetrapyrrole metabolism in bacteriophage biology.

Ploidy levels in diverse picocyanobacteria from the Baltic Sea

In nature, the number of genome or chromosome copies within cells (ploidy) can vary between species and environmental conditions, potentially influencing how organisms adapt to changing environments. Although ploidy levels cannot be easily determined by standard genome sequencing, understanding ploidy is crucial for the quantitative interpretation of molecular data. Cyanobacteria are known to contain haploid, oligoploid, and polyploid species. The smallest cyanobacteria, picocyanobacteria (less than 2 μm in diameter), have a widespread distribution ranging from marine to freshwater environments, contributing significantly to global primary production. In this study, we determined the ploidy level of genetically and physiologically diverse brackish picocyanobacteria isolated from the Baltic Sea using a qPCR assay targeting the rbcL gene. The strains contained one to four genome copies per cell. The ploidy level was not linked with phylogeny based on the identity of the 16S rRNA gene. The variation of ploidy among the brackish strains was lower compared to what has been reported for freshwater strains and was more similar to what has been reported for marine strains. The potential ecological advantage of polyploidy among picocyanobacteria has yet to be described. Our study highlights the importance of considering ploidy to interpret the abundance and adaptation of brackish picocyanobacteria.

Accounting for cellular-level variation in lysis: implications for virus-host dynamics

Viral impacts on microbial populations depend on interaction phenotypes-including viral traits spanning the adsorption rate, latent period, and burst size. The latent period is a key viral trait in lytic infections. Defined as the time from viral adsorption to viral progeny release, the latent period of bacteriophage is conventionally inferred via one-step growth curves in which the accumulation of free virus is measured over time in a population of infected cells. Developed more than 80 years ago, one-step growth curves do not account for cellular-level variability in the timing of lysis, potentially biasing inference of viral traits. Here, we use nonlinear dynamical models to understand how individual-level variation of the latent period impacts virus-host dynamics. Our modeling approach shows that inference of the latent period via one-step growth curves is systematically biased-generating estimates of shorter latent periods than the underlying population-level mean. The bias arises because variability in lysis timing at the cellular level leads to a fraction of early burst events, which are interpreted, artefactually, as an earlier mean time of viral release. We develop a computational framework to estimate latent period variability from joint measurements of host and free virus populations. Our computational framework recovers both the mean and variance of the latent period within simulated infections including realistic measurement noise. This work suggests that reframing the latent period as a distribution to account for variability in the population will improve the study of viral traits and their role in shaping microbial populations.IMPORTANCEQuantifying viral traits-including the adsorption rate, burst size, and latent period-is critical to characterize viral infection dynamics and develop predictive models of viral impacts across scales from cells to ecosystems. Here, we revisit the gold standard of viral trait estimation-the one-step growth curve-to assess the extent to which assumptions at the core of viral infection dynamics lead to ongoing and systematic biases in inferences of viral traits. We show that latent period estimates obtained via one-step growth curves systematically underestimate the mean latent period and, in turn, overestimate the rate of viral killing at population scales. By explicitly incorporating trait variability into a dynamical inference framework that leverages both virus and host time series, we provide a practical route to improve estimates of the mean and variance of viral traits across diverse virus-microbe systems.

Profiling the interplay and coevolution of Microcystis aeruginosa and cyanosiphophage Mic1

The cyanosiphophage Mic1 specifically infects the bloom-forming Microcystis aeruginosa FACHB 1339 from Lake Chaohu, China. Previous genomic analysis showed that its 92,627 bp double-stranded DNA genome consists of 98 putative open reading frames, 63% of which are of unknown function. Here, we investigated the transcriptome dynamics of Mic1 and its host using RNA sequencing. In the early, middle, and late phases of the 10 h lytic cycle, the Mic1 genes are sequentially expressed and could be further temporally grouped into two distinct clusters in each phase. Notably, six early genes, including gp49 that encodes a TnpB-like transposase, immediately reach the highest transcriptional level in half an hour, representing a pioneer cluster that rapidly regulates and redirects host metabolism toward the phage. An in-depth analysis of the host transcriptomic profile in response to Mic1 infection revealed significant upregulation of a polyketide synthase pathway and a type III-B CRISPR system, accompanied by moderate downregulation of the photosynthesis and key metabolism pathways. The constant increase of phage transcripts and relatively low replacement rate over the host transcripts indicated that Mic1 utilizes a unique strategy to gradually take over a small portion of host metabolism pathways after infection. In addition, genomic analysis of a less-infective Mic1 and a Mic1-resistant host strain further confirmed their dynamic interplay and coevolution via the frequent horizontal gene transfer. These findings provide insights into the mutual benefit and symbiosis of the highly polymorphic cyanobacteria M. aeruginosa and cyanophages.

IMPORTANCE

The highly polymorphic Microcystis aeruginosa is one of the predominant bloom-forming cyanobacteria in eutrophic freshwater bodies and is infected by diverse and abundant cyanophages. The presence of a large number of defense systems in M. aeruginosa genome suggests a dynamic interplay and coevolution with the cyanophages. In this study, we investigated the temporal gene expression pattern of Mic1 after infection and the corresponding transcriptional responses of its host. Moreover, the identification of a less-infective Mic1 and a Mic1-resistant host strain provided the evolved genes in the phage-host coevolution during the multiple-generation cultivation in the laboratory. Our findings enrich the knowledge on the interplay and coevolution of M. aeruginosa and its cyanophages and lay the foundation for the future application of cyanophage as a potential eco-friendly and bio-safe agent in controlling the succession of harmful cyanobacterial blooms.

Exploring the transcriptional landscape of phage-host interactions using novel high-throughput approaches

In the last decade, powerful high-throughput sequencing approaches have emerged to analyse microbial transcriptomes at a global scale. However, to date, applications of these approaches to microbial viruses such as phages remain scarce. Tailoring these techniques to virus-infected bacteria promises to obtain a detailed picture of the underexplored RNA biology and molecular processes during infection. In addition, transcriptome study of stress and perturbations induced by phages in their infected bacterial hosts is likely to reveal new fundamental mechanisms of bacterial metabolism and gene regulation. Here, we provide references and blueprints to implement emerging transcriptomic approaches towards addressing transcriptome architecture, RNA-RNA and RNA-protein interactions, RNA modifications, structures and heterogeneity of transcription profiles in infected cells that will provide guides for future directions in phage-centric therapeutic applications and microbial synthetic biology.

Environment-specific virocell metabolic reprogramming

Viruses impact microbial systems through killing hosts, horizontal gene transfer, and altering cellular metabolism, consequently impacting nutrient cycles. A virus-infected cell, a "virocell," is distinct from its uninfected sister cell as the virus commandeers cellular machinery to produce viruses rather than replicate cells. Problematically, virocell responses to the nutrient-limited conditions that abound in nature are poorly understood. Here we used a systems biology approach to investigate virocell metabolic reprogramming under nutrient limitation. Using transcriptomics, proteomics, lipidomics, and endo- and exo-metabolomics, we assessed how low phosphate (low-P) conditions impacted virocells of a marine Pseudoalteromonas host when independently infected by two unrelated phages (HP1 and HS2). With the combined stresses of infection and nutrient limitation, a set of nested responses were observed. First, low-P imposed common cellular responses on all cells (virocells and uninfected cells), including activating the canonical P-stress response, and decreasing transcription, translation, and extracellular organic matter consumption. Second, low-P imposed infection-specific responses (for both virocells), including enhancing nitrogen assimilation and fatty acid degradation, and decreasing extracellular lipid relative abundance. Third, low-P suggested virocell-specific strategies. Specifically, HS2-virocells regulated gene expression by increasing transcription and ribosomal protein production, whereas HP1-virocells accumulated host proteins, decreased extracellular peptide relative abundance, and invested in broader energy and resource acquisition. These results suggest that although environmental conditions shape metabolism in common ways regardless of infection, virocell-specific strategies exist to support viral replication during nutrient limitation, and a framework now exists for identifying metabolic strategies of nutrient-limited virocells in nature.

Increased cyanophage infection at the bottom of the euphotic zone, especially in the fall

Picocyanobacteria contribute greatly to offshore primary production with cells extending through the deep euphotic zone. Literature indicates high viral infection of cyanobacteria in ocean transition zones. We postulate that the bottom of the euphotic zone is a transition zone, where communities transition from phototrophic to aphotic processes. We use single-copy core genes to examine cyanophage to cyanobacteria ratios in cellular metagenomes in the subtropical North Atlantic and Pacific. Cyanophage to cyanobacteria terL/rpoB ratios generally increase to >10 in the deep euphotic zone. As light levels decrease in the fall, Prochlorococcus in the deep euphotic zone experience reduced light levels. We find clear differences between spring (Geotraces GA02) and fall (GA03) in the North Atlantic, with terL/rpoB ratios increasing to >40 in the fall. When examining 23 months of the North Pacific Hawaii Ocean Timeseries, the depth of elevated cyanophage to cyanobacteria ratios in cellular metagenomes negatively correlated with surface photosynthetic radiation (PAR), particularly with the change in PAR, which reflected the season. In fall, all picocyanobacteria ecotypes were found at depths enriched with viruses, while in summer, only low light ecotypes were affected. Thus, we find high cyanophage infection both in the deep euphotic zone and during seasonal transitions.

Top-heavy trophic structure within benthic viral dark matter

A better understanding of system-specific viral ecology in diverse environments is needed to predict patterns of virus-host trophic structure in the Anthropocene. This study characterised viral-host trophic structure within coral reef benthic cyanobacterial mats-a globally proliferating cause and consequence of coral reef degradation. We employed deep longitudinal multi-omic sequencing to characterise the viral assemblage (ssDNA, dsDNA, and dsRNA viruses) and profile lineage-specific host-virus interactions within benthic cyanobacterial mats sampled from Bonaire, Caribbean Netherlands. We recovered 11,012 unique viral populations spanning at least 10 viral families across the orders Caudovirales, Petitvirales, and Mindivirales. Gene-sharing network analyses provided evidence for extensive genomic novelty of mat viruses from reference and environmental viral sequences. Analysis of coverage ratios of viral sequences and computationally predicted hosts spanning 15 phyla and 21 classes revealed virus-host abundance (from DNA) and activity (from RNA) ratios consistently exceeding 1:1, suggesting a top-heavy intra-mat trophic structure with respect to virus-host interactions. Overall, our article contributes a curated database of viral sequences found in Caribbean coral reef benthic cyanobacterial mats (vMAT database) and provides multiple lines of field-based evidence demonstrating that viruses are active members of mat communities, with broader implications for mat functional ecology and demography.

Characterization and genomic analysis of an oceanic cyanophage infecting marine Synechococcus reveal a novel genus

Cyanophages play a crucial role in the biogeochemical cycles of aquatic ecosystems by affecting the population dynamics and community structure of cyanobacteria. In this study, a novel cyanophage, Nanhaivirus ms29, that infects Synechococcus sp. MW02 was isolated from the ocean basin in the South China Sea. It was identified as a T4-like phage using transmission electron microscopy. Phylogenetic analysis demonstrated that this cyanophage is distinct from other known T4-like cyanophage, belonging to a novel genus named Nanhaivirus within the family Kyanoviridae, according to the most recent classification proposed by the International Committee on Taxonomy of Viruses (ICTV). The genome of this novel cyanophage is composed of 178,866 bp of double-stranded DNA with a G + C content of 42.5%. It contains 217 potential open reading frames (ORFs) and 6 tRNAs. As many as 30 auxiliary metabolic genes (AMGs) were identified in the genome, which related to photosynthesis, carbon metabolism, nutrient uptake and stress tolerance, possibly reflecting a genomic adaption to the oligotrophic environment. Read-mapping analysis showed that Nanhaivirus ms29 mainly distributed in temperate and tropical epipelagic waters. This study enriches of the virus gene database of cyanophages and provides valuable insights into the phylogeny of cyanophages and their interactions with their hosts.

From natural to artificial cyanophages: Current progress and application prospects

The over proliferation of harmful cyanobacteria and their cyanotoxins resulted in damaged aquatic ecosystem, polluted drinking water and threatened human health. Cyanophages are a kind of viruses that exclusively infect cyanobacteria, which is considered as a potential strategy to deal with cyanobacterial blooms. Nevertheless, the infecting host range and/or lysis efficiency of natural cyanophages is limited, rising the necessity of constructing non-natural cyanophages via artificial modification, design and synthesis to expand their host range and/or efficiency. The paper firstly reviewed representative cyanophages such as P60 with a short latent period of 1.5 h and S-CBS1 having a burst size up to 200 PFU/cell. To explore the in-silico design principles, we critically summarized the interactions between cyanophages and the hosts, indicating modifying the recognized receptors, enhancing the adsorption ability, changing the lysogeny and excluding the defense of hosts are important for artificial cyanophages. The research progress of synthesizing artificial cyanophages were summarized subsequently, raising the importance of developing genetic manipulation technologies and their rescue strategies in the future. Meanwhile, Large-scale preparation of cyanophages for bloom control is a big challenge. The application prospects of artificial cyanophages besides cyanobacteria bloom control like adaptive evolution and phage therapy were discussed at last. The review will promote the design, synthesis and application of cyanophages for cyanobacteria blooms, which may provide new insights for the related water pollution control and ensuring hydrosphere security.

Cyanophage technology in removal of cyanobacteria mediated harmful algal blooms: A novel and eco-friendly method

Cyanophages are highly abundant specific viruses that infect cyanobacterial cells. In recent years, the cyanophages and cyanobacteria interactions drew attention to environmental restoration due to their discovery in marine and freshwater systems. Cyanobacterial harmful algal blooms (cyanoHABs) are increasing throughout the world and contaminating aquatic ecosystems. The blooms cause severe environmental problems including unpleasant odors and cyanotoxin production. Cyanotoxins have been reported to be lethal agents for living beings and can harm animals, people, aquatic species, recreational activities, and drinking water reservoirs. Biological remediation of cyanoHABs in aquatic systems is a sustainable and eco-friendly approach to increasing surface water quality. Therefore, this study compiles the fragmented information with the solution of removal of cyanoHABs using cyanophage therapy techniques. To date, scant information exists in terms of bloom formation, cyanophage occurrence, and mode of action to remediate cyanoHABs. Overall, this study illustrates cyanobacterial toxin production and its impacts on the environment, the mechanisms involved in the cyanophage-cyanobacteria interaction, and the application of cyanophages for the removal of toxic cyanobacterial blooms.

An Estuarine Cyanophage S-CREM1 Encodes Three Distinct Antitoxin Genes and a Large Number of Non-Coding RNA Genes

Cyanophages play important roles in regulating the population dynamics, community structure, metabolism, and evolution of cyanobacteria in aquatic ecosystems. Here, we report the genomic analysis of an estuarine cyanophage, S-CREM1, which represents a new genus of T4-like cyanomyovirus and exhibits new genetic characteristics. S-CREM1 is a lytic phage which infects estuarine Synechococcus sp. CB0101. In contrast to many cyanomyoviruses that usually have a broad host range, S-CREM1 only infected the original host strain. In addition to cyanophage-featured auxiliary metabolic genes (AMGs), S-CREM1 also contains unique AMGs, including three antitoxin genes, a MoxR family ATPase gene, and a pyrimidine dimer DNA glycosylase gene. The finding of three antitoxin genes in S-CREM1 implies a possible phage control of host cells during infection. One small RNA (sRNA) gene and three cis-regulatory RNA genes in the S-CREM1 genome suggest potential molecular regulations of host metabolism by the phage. In addition, S-CREM1 contains a large number of tRNA genes which may reflect a genomic adaption to the nutrient-rich environment. Our study suggests that we are still far from understanding the viral diversity in nature, and the complicated virus-host interactions remain to be discovered. The isolation and characterization of S-CREM1 further our understanding of the gene diversity of cyanophages and phage-host interactions in the estuarine environment.

Associations between picocyanobacterial ecotypes and cyanophage host genes across ocean basins and depth

Background

Cyanophages, viruses that infect cyanobacteria, are globally abundant in the ocean's euphotic zone and are a potentially important cause of mortality for marine picocyanobacteria. Viral host genes are thought to increase viral fitness by either increasing numbers of genes for synthesizing nucleotides for virus replication, or by mitigating direct stresses imposed by the environment. The encoding of host genes in viral genomes through horizontal gene transfer is a form of evolution that links viruses, hosts, and the environment. We previously examined depth profiles of the proportion of cyanophage containing various host genes in the Eastern Tropical North Pacific Oxygen Deficient Zone (ODZ) and at the subtropical North Atlantic (BATS). However, cyanophage host genes have not been previously examined in environmental depth profiles across the oceans.

Methodology

We examined geographical and depth distributions of picocyanobacterial ecotypes, cyanophage, and their viral-host genes across ocean basins including the North Atlantic, Mediterranean Sea, North Pacific, South Pacific, and Eastern Tropical North and South Pacific ODZs using phylogenetic metagenomic read placement. We determined the proportion of myo and podo-cyanophage containing a range of host genes by comparing to cyanophage single copy core gene terminase (terL). With this large dataset (22 stations), network analysis identified statistical links between 12 of the 14 cyanophage host genes examined here with their picocyanobacteria host ecotypes.

Results

Picyanobacterial ecotypes, and the composition and proportion of cyanophage host genes, shifted dramatically and predictably with depth. For most of the cyanophage host genes examined here, we found that the composition of host ecotypes predicted the proportion of viral host genes harbored by the cyanophage community. Terminase is too conserved to illuminate the myo-cyanophage community structure. Cyanophage cobS was present in almost all myo-cyanophage and did not vary in proportion with depth. We used the composition of cobS phylotypes to track changes in myo-cyanophage composition.

Conclusions

Picocyanobacteria ecotypes shift with changes in light, temperature, and oxygen and many common cyanophage host genes shift concomitantly. However, cyanophage phosphate transporter gene pstS appeared to instead vary with ocean basin and was most abundant in low phosphate regions. Abundances of cyanophage host genes related to nutrient acquisition may diverge from host ecotype constraints as the same host can live in varying nutrient concentrations. Myo-cyanophage community in the anoxic ODZ had reduced diversity. By comparison to the oxic ocean, we can see which cyanophage host genes are especially abundant (nirA, nirC, and purS) or not abundant (myo psbA) in ODZs, highlighting both the stability of conditions in the ODZ and the importance of nitrite as an N source to ODZ endemic LLV Prochlorococcus.

Diversity in Infection Specificity between the Bloom-forming Microalga Heterosigma akashiwo and Its dsDNA Virus, Heterosigma akashiwo Virus

Heterosigma akashiwo virus (HaV) is a dsDNA virus that infects the bloom-forming raphidoflagellate Heterosigma akashiwo. Both the host and its virus are phenotypically diverse in terms of infection specificity. Their relationships have been examined based on the occurrence or absence of algal lysis following virus inoculation; however, variations in the strain-level host-virus relationship regarding infectivity and lysis rates remain unclear. Therefore, we performed a series of cross-infectivity tests using 60 H. akashiwo and 22 HaV strains isolated from the coastal waters of western Japan. The host strains were divided into 5 different groups and viruses into 4 groups. Using a representative strain from each group, algal lysis was observed in 14 of the (5×4=) 20 host-virus combinations; the concentration of infectious units in each HaV suspension was then assessed using the most probable number (MPN) assay on the five host strains. Virus titers ranged between 1.1×10¹ and 2.1×10⁷ infectious units mL^-1; the titer of each viral lysate was differently estimated using distinct H. akashiwo strains as hosts. These results suggest that (1) a clonal viral lysate comprises virions with different intraspecific infection specificities and/or (2) the efficiency and error rates of each intracellular replication process vary in each host-virus combination.

A Novel Inovirus Reprograms Metabolism and Motility of Marine Alteromonas

Members from the Inoviridae family with striking features are widespread, highly diverse, and ecologically pervasive across multiple hosts and environments. However, a small number of inoviruses have been isolated and studied. Here, a filamentous phage infecting Alteromonas abrolhosensis, designated ϕAFP1, was isolated from the South China Sea and represented a novel genus of Inoviridae. ϕAFP1 consisted of a single-stranded DNA genome (5986 bp), encoding eight putative ORFs. Comparative analyses revealed ϕAFP1 could be regarded as genetic mosaics having homologous sequences with Ralstonia and Stenotrophomonas phages. The temporal transcriptome analysis of A. abrolhosensis to ϕAFP1 infection revealed that 7.78% of the host genes were differentially expressed. The genes involved in translation processes, ribosome pathways, and degradation of multiple amino acid pathways at the plateau period were upregulated, while host material catabolic and bacterial motility-related genes were downregulated, indicating that ϕAFP1 might hijack the energy of the host for the synthesis of phage proteins. ϕAFP1 exerted step-by-step control on host genes through the appropriate level of utilizing host resources. Our study provided novel information for a better understanding of filamentous phage characteristics and phage-host interactions. IMPORTANCE Alteromonas is widely distributed and plays a vital role in biogeochemical in marine environments. However, little information about Alteromonas phages is available. Here, we isolated and characterized the biological characteristics and genome sequence of a novel inovirus infecting Alteromonas abrolhosensis, designated ϕAFP1, representing a novel viral genus of Inoviridae. We then presented a comprehensive view of the ϕAFP1 phage-Alteromonas abrolhosensis interactions, elucidating reprogramed host metabolism and motility. Our study provided novel information for better comprehension of filamentous phage characteristics and phage-host interactions.

Protist impacts on marine cyanovirocell metabolism

The fate of oceanic carbon and nutrients depends on interactions between viruses, prokaryotes, and unicellular eukaryotes (protists) in a highly interconnected planktonic food web. To date, few controlled mechanistic studies of these interactions exist, and where they do, they are largely pairwise, focusing either on viral infection (i.e., virocells) or protist predation. Here we studied population-level responses of Synechococcus cyanobacterial virocells (i.e., cyanovirocells) to the protist Oxyrrhis marina using transcriptomics, endo- and exo-metabolomics, photosynthetic efficiency measurements, and microscopy. Protist presence had no measurable impact on Synechococcus transcripts or endometabolites. The cyanovirocells alone had a smaller intracellular transcriptional and metabolic response than cyanovirocells co-cultured with protists, displaying known patterns of virus-mediated metabolic reprogramming while releasing diverse exometabolites during infection. When protists were added, several exometabolites disappeared, suggesting microbial consumption. In addition, the intracellular cyanovirocell impact was largest, with 4.5- and 10-fold more host transcripts and endometabolites, respectively, responding to protists, especially those involved in resource and energy production. Physiologically, photosynthetic efficiency also increased, and together with the transcriptomics and metabolomics findings suggest that cyanovirocell metabolic demand is highest when protists are present. These data illustrate cyanovirocell responses to protist presence that are not yet considered when linking microbial physiology to global-scale biogeochemical processes.

Cyanophage infections reduce photosynthetic activity and expression of CO2 fixation genes in the freshwater bloom-forming cyanobacterium Aphanizomenon flos-aquae

Cyanobacteria play a significant role in ecosystem functioning as photosynthetic and CO₂ fixing microorganisms. Whether and to what extent cyanophages alter these carbon and energy cycles in their cyanobacterial hosts is still poorly understood. In this study, we investigated changes in photosynthetic activity (PSII), expression of genes associated with the light phase of photosynthesis (psbA, petA, ndhK) and carbon metabolism (rbcL, zwf) as well as intracellular ATP and NADHP concentrations in freshwater bloom-forming filamentous cyanobacterium Aphanizomenon flos-aquae infected by cyanophage vB_AphaS-CL131. We found that PSII activity and expression level of rbcL genes, indicating potential for CO₂ fixation, had decreased in response to cyanophage adsorption and DNA injection. During the period of viral DNA replication and assembly, PSII performance and gene expression remained at this decreased level and did not change significantly, indicating lack of transcriptional shutdown by the cyanophage. Combined, these observations suggest that although there is little to no interference between cyanophage DNA replication, host transcription and cellular metabolism, A. flos-aquae underwent a physiological state-shift toward lower efficiency of carbon and energy cycling. This further suggest potential cascading effect for co-occurring non-infected members of the microbial community.

Global Transcriptomic Analysis of Bacteriophage-Host Interactions between a Kayvirus Therapeutic Phage and Staphylococcus aureus

Kayviruses are polyvalent broad host range staphylococcal phages with a potential to combat staphylococcal infections. However, the implementation of rational phage therapy in medicine requires a thorough understanding of the interactions between bacteriophages and pathogens at omics level. To evaluate the effect of a phage used in therapy on its host bacterium, we performed differential transcriptomic analysis by RNA-Seq from bacteriophage K of genus Kayvirus infecting two Staphylococcus aureus strains, prophage-less strain SH1000 and quadruple lysogenic strain Newman. The temporal transcriptional profile of phage K was comparable in both strains except for a few loci encoding hypothetical proteins. Stranded sequencing revealed transcription of phage noncoding RNAs that may play a role in the regulation of phage and host gene expression. The transcriptional response of S. aureus to phage K infection resembles a general stress response with differential expression of genes involved in a DNA damage response. The host transcriptional changes involved upregulation of nucleotide, amino acid and energy synthesis and transporter genes and downregulation of host transcription factors. The interaction of phage K with variable genetic elements of the host showed slight upregulation of gene expression of prophage integrases and antirepressors. The virulence genes involved in adhesion and immune evasion were only marginally affected, making phage K suitable for therapy. IMPORTANCE Bacterium Staphylococcus aureus is a common human and veterinary pathogen that causes mild to life-threatening infections. As strains of S. aureus are becoming increasingly resistant to multiple antibiotics, the need to search for new therapeutics is urgent. A promising alternative to antibiotic treatment of staphylococcal infections is a phage therapy using lytic phages from the genus Kayvirus. Here, we present a comprehensive view on the phage-bacterium interactions on transcriptomic level that improves the knowledge of molecular mechanisms underlying the Kayvirus lytic action. The results will ensure safer usage of the phage therapeutics and may also serve as a basis for the development of new antibacterial strategies.

Picocyanobacterial Synechococcus in marine ecosystem: Insights from genetic diversity, global distribution, and potential function

Marine Synechococcus, a main group of picocyanobacteria, has been ubiquitously observed across the global oceans. Synechococcus exhibits high phylogenetical and phenotypical diversity, and horizontal gene transfer makes its genetic evolution much more intricate. With the development of measurement technologies and analysis methods, the genomic information and niche partition of each Synechococcus lineage tend to be precisely described, but the global analysis is still lacking. Therefore, it is necessary to summarize existing studies and integrate published data to gain a comprehensive understanding of Synechococcus on genetic variation, niche division, and potential functions. In this review, the maximum likelihood trees are constructed based on existing sequence data, including both phylogenetic and pigmentary gene markers. The global distribution characteristics of abundance, lineages, and pigment types are concluded through pooled analysis of more than 700 samples obtained from approximately 50 scientific research cruises. The potential functions of Synechococcus are explored in element cycles and biological interactions. Future work on Synechococcus is suggested to focus on not only elucidating the nature of Synechococcus biodiversity but also demonstrating its interactions with the ecosystem by combining bioinformatics and macroscopic isotope-labeled environmental parameters.

Nutrient driven transcriptional changes during phage infection in an aquatic Gammaproteobacterium

Phages modulate bacterial metabolism during infection by regulating gene expression, which influences aquatic nutrient cycling. However, the effects of shifting nutrient regimes are less understood. Here, we analyzed transcriptomes of an ecologically relevant Gammaproteobacterium and its lytic phage in high (HNM) and low (LNM) nutrient medium. Despite different infection characteristics, including reduced burst size and longer latent period in LNM, the phage had a fixed expression profile. Bacterial transcription was instead different depending on nutrient regime, with HNM bacteria focusing on growth while LNM bacteria focused on motility and membrane transport. Additionally, phage infection had a larger effect on bacterial gene expression in LNM compared to HNM, e.g. suppressing increased iron uptake and altering expression of phosphorus uptake genes. Overall, phage infection influenced host metabolism more in LNM, which was more similar to natural conditions, emphasizing the importance of considering natural conditions to understand phage and host ecology.

Isolation and Characterization of a Novel Cyanophage Encoding Multiple Auxiliary Metabolic Genes

As significant drivers of cyanobacteria mortality, cyanophages have been known to regulate the population dynamics, metabolic activities, and community structure of this most important marine autotrophic picoplankton and, therefore, influence the global primary production and biogeochemical cycle in aquatic ecosystems. In the present study, a lytic Synechococcus phage, namely S-SZBM1, was isolated and identified. Cyanophage S-SZBM1 has a double-stranded DNA genome of 177,834 bp with a G+C content of 43.31% and contains a total of 218 predicted ORFs and six tRNA genes. Phylogenetic analysis and nucleotide-based intergenomic similarity suggested that cyanophage S-SZBM1 belongs to a new genus under the family Kyanoviridae. A variety of auxiliary metabolic genes (AMGs) that have been proved or speculated to relate to photosynthesis, carbon metabolism, nucleotide synthesis and metabolism, cell protection, and other cell metabolism were identified in cyanophage S-SZBM1 genome and may affect host processes during infection. In addition, 24 of 32 predicted structural proteins were identified by a high-throughput proteome analysis which were potentially involved in the assembly processes of virion. The genomic and proteomic analysis features of cyanophage S-SZBM1 offer a valuable insight into the interactions between cyanophages and their hosts during infection.

A Novel and Ubiquitous Marine Methylophage Provides Insights into Viral-Host Coevolution and Possible Host-Range Expansion in Streamlined Marine Heterotrophic Bacteria

The methylotrophic OM43 clade are Gammaproteobacteria that comprise some of the smallest free-living cells known and have highly streamlined genomes. OM43 represents an important microbial link between marine primary production and remineralization of carbon back to the atmosphere. Bacteriophages shape microbial communities and are major drivers of mortality and global marine biogeochemistry. Recent cultivation efforts have brought the first viruses infecting members of the OM43 clade into culture. Here, we characterize a novel myophage infecting OM43 called Melnitz. Melnitz was isolated independently from water samples from a subtropical ocean gyre (Sargasso Sea) and temperate coastal (Western English Channel) systems. Metagenomic recruitment from global ocean viromes confirmed that Melnitz is globally ubiquitous, congruent with patterns of host abundance. Bacteria with streamlined genomes such as OM43 and the globally dominant SAR11 clade use riboswitches as an efficient method to regulate metabolism. Melnitz encodes a two-piece tmRNA (ssrA), controlled by a glutamine riboswitch, providing evidence that riboswitch use also occurs for regulation during phage infection of streamlined heterotrophs. Virally encoded tRNAs and ssrA found in Melnitz were phylogenetically more closely related to those found within the alphaproteobacterial SAR11 clade and their associated myophages than those within their gammaproteobacterial hosts. This suggests the possibility of an ancestral host transition event between SAR11 and OM43. Melnitz and a related myophage that infects SAR11 were unable to infect hosts of the SAR11 and OM43, respectively, suggesting host transition rather than a broadening of host range. IMPORTANCE Isolation and cultivation of viruses are the foundations on which the mechanistic understanding of virus-host interactions and parameterization of bioinformatic tools for viral ecology are based. This study isolated and characterized the first myophage known to infect the OM43 clade, expanding our knowledge of this understudied group of microbes. The nearly identical genomes of four strains of Melnitz isolated from different marine provinces and the global abundance estimations from metagenomic data suggest that this viral population is globally ubiquitous. Genome analysis revealed several unusual features in Melnitz and related genomes recovered from viromes, such as a curli operon and virally encoded tmRNA controlled by a glutamine riboswitch, neither of which are found in the host. Further phylogenetic analysis of shared genes indicates that this group of viruses infecting the gammaproteobacterial OM43 shares a recent common ancestor with viruses infecting the abundant alphaproteobacterial SAR11 clade. Host ranges are affected by compatible cell surface receptors, successful circumvention of superinfection exclusion systems, and the presence of required accessory proteins, which typically limits phages to singular narrow groups of closely related bacterial hosts. This study provides intriguing evidence that for streamlined heterotrophic bacteria, virus-host transitioning may not be necessarily restricted to phylogenetically related hosts but is a function of shared physical and biochemical properties of the cell.

Genetic engineering of marine cyanophages reveals integration but not lysogeny in T7-like cyanophages

Marine cyanobacteria of the genera Synechococcus and Prochlorococcus are the most abundant photosynthetic organisms on earth, spanning vast regions of the oceans and contributing significantly to global primary production. Their viruses (cyanophages) greatly influence cyanobacterial ecology and evolution. Although many cyanophage genomes have been sequenced, insight into the functional role of cyanophage genes is limited by the lack of a cyanophage genetic engineering system. Here, we describe a simple, generalizable method for genetic engineering of cyanophages from multiple families, that we named REEP for REcombination, Enrichment and PCR screening. This method enables direct investigation of key cyanophage genes, and its simplicity makes it adaptable to other ecologically relevant host-virus systems. T7-like cyanophages often carry integrase genes and attachment sites, yet exhibit lytic infection dynamics. Here, using REEP, we investigated their ability to integrate and maintain a lysogenic life cycle. We found that these cyanophages integrate into the host genome and that the integrase and attachment site are required for integration. However, stable lysogens did not form. The frequency of integration was found to be low in both lab cultures and the oceans. These findings suggest that T7-like cyanophage integration is transient and is not part of a classical lysogenic cycle.

Type I CRISPR-Cas provides robust immunity but incomplete attenuation of phage-induced cellular stress

During infection, phages manipulate bacteria to redirect metabolism towards viral proliferation. To counteract phages, some bacteria employ CRISPR-Cas systems that provide adaptive immunity. While CRISPR-Cas mechanisms have been studied extensively, their effects on both the phage and the host during phage infection remains poorly understood. Here, we analysed the infection of Serratia by a siphovirus (JS26) and the transcriptomic response with, or without type I-E or I-F CRISPR-Cas immunity. In non-immune Serratia, phage infection altered bacterial metabolism by upregulating anaerobic respiration and amino acid biosynthesis genes, while flagella production was suppressed. Furthermore, phage proliferation required a late-expressed viral Cas4 homologue, which did not influence CRISPR adaptation. While type I-E and I-F immunity provided robust defence against phage infection, phage development still impacted the bacterial host. Moreover, DNA repair and SOS response pathways were upregulated during type I immunity. We also discovered that the type I-F system is controlled by a positive autoregulatory feedback loop that is activated upon phage targeting during type I-F immunity, leading to a controlled anti-phage response. Overall, our results provide new insight into phage-host dynamics and the impact of CRISPR immunity within the infected cell.

Confocal microscopy reveals alterations of thylakoids in Limnospira fusiformis during prophage induction

The alkaliphilic cyanobacterium Limnospira fusiformis is an integral part in food webs of tropical soda lakes. Recently, sudden breakdowns of Limnospira sp. blooms in their natural environment have been linked to cyanophage infections. We studied ultrastructural details and prophage components in the laboratory by means of confocal laser scanning microscopy (CLSM) and transmission electron microscopy (TEM). For a comparison at the subcellular level, we included transmission electron microscopy (TEM) material of infected cells collected during a field survey. Compared to TEM, CLSM has the advantage to rapidly providing results for whole, intact cells. Moreover, many cells can be studied at once. We chemically induced lysogenic cyanophages by means of mitomycin C (MMC) treatments and studied the ultrastructural alterations of host cells. In parallel, the number of cyanophages was obtained by flow cytometry. After treatment of the culture with MMC, flow cytometry showed a strong increase in viral counts, i.e., prophage induction. CLSM reflected the re-organization of L. fusiformis with remarkable alterations of thylakoid arrangements after prophage induction. Our study provides a first step towards 3D visualization of ultrastructure of cyanobacteria and showed the high potential of CLSM to investigate viral-mediated modifications in these groups.

Hostile Takeover: How Viruses Reprogram Prokaryotic Metabolism

To reproduce, prokaryotic viruses must hijack the cellular machinery of their hosts and redirect it toward the production of viral particles. While takeover of the host replication and protein synthesis apparatus has long been considered an essential feature of infection, recent studies indicate that extensive reprogramming of host primary metabolism is a widespread phenomenon among prokaryotic viruses that is required to fulfill the biosynthetic needs of virion production. In this review we provide an overview of the most significant recent findings regarding virus-induced reprogramming of prokaryotic metabolism and suggest how quantitative systems biology approaches may be used to provide a holistic understanding of metabolic remodeling during lytic viral infection. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Frequency of mispackaging of Prochlorococcus DNA by cyanophage

Prochlorococcus cells are the numerically dominant phototrophs in the open ocean. Cyanophages that infect them are a notable fraction of the total viral population in the euphotic zone, and, as vehicles of horizontal gene transfer, appear to drive their evolution. Here we examine the propensity of three cyanophages-a podovirus, a siphovirus, and a myovirus-to mispackage host DNA in their capsids while infecting Prochlorococcus, the first step in phage-mediated horizontal gene transfer. We find the mispackaging frequencies are distinctly different among the three phages. Myoviruses mispackage host DNA at low and seemingly fixed frequencies, while podo- and siphoviruses vary in their mispackaging frequencies by orders of magnitude depending on growth light intensity. We link this difference to the concentration of intracellular reactive oxygen species and protein synthesis rates, both parameters increasing in response to higher light intensity. Based on our findings, we propose a model of mispackaging frequency determined by the imbalance between the production of capsids and the number of phage genome copies during infection: when protein synthesis rate increase to levels that the phage cannot regulate, they lead to an accumulation of empty capsids, in turn triggering more frequent host DNA mispackaging errors.

Temporal transcriptomes of a marine cyanopodovirus and its Synechococcus host during infection

Marine picocyanobacteria belonging to genera Synechococcus and Prochlorococcus are genetically diverged and distributed into distinct biogeographical patterns, and both are infected by genetically closely related cyanopodoviruses. Previous studies have not fully explored whether the two virus–host systems share similar gene expression patterns during infection. Whole‐genome expression dynamics of T7‐like cyanopodovirus P‐SSP7 and its host Prochlorococcus strain MED4 have already been reported. Here, we conducted genomic and transcriptomic analyses on T7‐like cyanopodovirus S‐SBP1 during its infection on Synechococcus strain WH7803. S‐SBP1 has a latent period of 8 h and phage DNA production of 30 copies per cell. In terms of whole‐genome phylogenetic relationships and average nucleotide identity, S‐SBP1 was most similar to cyanopodovirus S‐RIP2, which also infects Synechococcus WH7803. Three hypervariable genomic islands were identified when comparing the genomes of S‐SBP1 and S‐RIP2. Single nucleotide variants were also observed in three S‐SBP1 genes, which were located within the island regions. Based on RNA‐seq analysis, S‐SBP1 genes clustered into three temporal expression classes, whose gene content was similar to that of P‐SSP7. Thirty‐two host genes were upregulated during phage infection, including those involved in carbon metabolism, ribosome components, and stress response. These upregulated genes were similar to those upregulated by Prochlorococcus MED4 in response to infection by P‐SSP7. Our study demonstrates a programmed temporal expression pattern of cyanopodoviruses and hosts during infection.

Phage gene expression and host responses lead to infection-dependent costs of CRISPR immunity

CRISPR-Cas immune systems are widespread in bacteria and archaea, but not ubiquitous. Previous work has demonstrated that CRISPR immunity is associated with an infection-induced fitness cost, which may help explain the patchy distribution observed. However, the mechanistic basis of this cost has remained unclear. Using Pseudomonas aeruginosa PA14 and its phage DMS3vir as a model, we perform a 30-day evolution experiment under phage mediated selection. We demonstrate that although CRISPR is initially selected for, bacteria carrying mutations in the phage receptor rapidly invade the population following subsequent reinfections. We then test three potential mechanisms for the observed cost of CRISPR: (1) autoimmunity from the acquisition of self-targeting spacers, (2) immunopathology or energetic costs from increased cas gene expression and (3) toxicity caused by phage gene expression prior to CRISPR-mediated cleavage. We find that phages can express genes before the immune system clears the infection and that expression of these genes can have a negative effect on host fitness. While infection does not lead to increased expression of cas genes, it does cause differential expression of multiple other host processes that may further contribute to the cost of CRISPR immunity. In contrast, we found little support for infection-induced autoimmunological and immunopathological effects. Phage gene expression prior to cleavage of the genome by the CRISPR-Cas immune system is therefore the most parsimonious explanation for the observed phage-induced fitness cost.

Viruses of freshwater bloom-forming cyanobacteria: genomic features, infection strategies and coexistence with the host

Freshwater bloom-forming cyanobacteria densely grow in the aquatic environments, leading to an increase in the viral-contact rate. They possess numerous antiviral genes, as well as cell differentiation- and physiological performance-related genes, owing to genome expansion. Their genomic features and unique lifestyles suggest that they coexist with cyanoviruses in ways different from marine cyanobacteria. Furthermore, genome contents of isolated freshwater bloom-forming cyanobacterial viruses have little in common with those of marine cyanoviruses studied to date. They lack the marine cyanoviral hallmark genes that sustain photosynthetic activity and redirect host metabolism to viral reproduction; therefore, they are predicted to share metabolisms and precursor pools with host cyanobacteria to ensure efficient viral reproduction and avoid nutrient deficiencies and antiviral response. Additionally, cyanovirus-cyanobacteria coexistence strategies may change as bloom density increases. Diverse genotypic populations of cyanoviruses and hosts coexist and fluctuate under high viral-contact rate conditions, leading to their rapid coevolution through antiviral responses. The ancestral and newly evolved genotypes coexist, thereby expanding the diversity levels of host and viral populations. Bottleneck events occurring due to season-related decreases in bloom-forming species abundance provide each genotype within cyanobacterial population an equal chance to increase in prevalence during the next bloom and enhance further diversification.

Lysogeny in the oceans: Lessons from cultivated model systems and a reanalysis of its prevalence

In the oceans, viruses that infect bacteria (phages) influence a variety of microbially mediated processes that drive global biogeochemical cycles. The nature of their influence is dependent upon infection mode, be it lytic or lysogenic. Temperate phages are predicted to be prevalent in marine systems where they are expected to execute both types of infection modes. Understanding the range and outcomes of temperate phage-host interactions is fundamental for evaluating their ecological impact. Here, we (i) review phage-mediated rewiring of host metabolism, with a focus on marine systems, (ii) consider the range and nature of temperate phage-host interactions, and (iii) draw on studies of cultivated model systems to examine the consequences of lysogeny among several dominant marine bacterial lineages. We also readdress the prevalence of lysogeny among marine bacteria by probing a collection of 1239 publicly available bacterial genomes, representing cultured and uncultivated strains, for evidence of complete prophages. Our conservative analysis, anticipated to underestimate true prevalence, predicts 18% of the genomes examined contain at least one prophage, the majority (97%) were found within genomes of cultured isolates. These results highlight the need for cultivation of additional model systems to better capture the diversity of temperate phage-host interactions in the oceans.

Nitrogen Flow in Diazotrophic Cyanobacterium Aphanizomenon flos-aquae Is Altered by Cyanophage Infection

Viruses can significantly influence cyanobacteria population dynamics and activity, and through this the biogeochemical cycling of major nutrients. However, surprisingly little attention has been given to understand how viral infections alter the ability of diazotrophic cyanobacteria for atmospheric nitrogen fixation and its release to the environment. This study addressed the importance of cyanophages for net ¹⁵N₂ assimilation rate, expression of nitrogenase reductase gene (nifH) and changes in nitrogen enrichment (¹⁵N/¹⁴N) in the diazotrophic cyanobacterium Aphanizomenon flos-aquae during infection by the cyanophage vB_AphaS-CL131. We found that while the growth of A. flos-aquae was inhibited by cyanophage addition (decreased from 0.02 h^–1 to 0.002 h^–1), there were no significant differences in nitrogen fixation rates (control: 22.7 × 10^–7 nmol N heterocyte^–1; infected: 23.9 × 10^–7 nmol N heterocyte^–1) and nifH expression level (control: 0.6–1.6 transcripts heterocyte^–1; infected: 0.7–1.1 transcripts heterocyte^–1) between the infected and control A. flos-aquae cultures. This implies that cyanophage genome replication and progeny production within the vegetative cells does not interfere with the N₂ fixation reactions in the heterocytes of these cyanobacteria. However, higher ¹⁵N enrichment at the poles of heterocytes of the infected A. flos-aquae, revealed by NanoSIMS analysis indicates the accumulation of fixed nitrogen in response to cyanophage addition. This suggests reduced nitrogen transport to vegetative cells and the alterations in the flow of fixed nitrogen within the filaments. In addition, we found that cyanophage lysis resulted in a substantial release of ammonium into culture medium. Cyanophage infection seems to substantially redirect N flow from cyanobacterial biomass to the production of N storage compounds and N release.

Predetermined clockwork microbial worlds: Current understanding of aquatic microbial diel response from model systems to complex environments

In the photic zone of aquatic ecosystems, microorganisms with different metabolisms and their viruses form complex interactions and food webs. Within these interactions, phototrophic microorganisms such as eukaryotic microalgae and cyanobacteria interact directly with sunlight, and thereby generate circadian rhythms. Diel cycling originally generated in microbial phototrophs is directly transmitted toward heterotrophic microorganisms utilizing the photosynthetic products as they are excreted or exuded. Such diel cycling seems to be indirectly propagated toward heterotrophs as a result of complex biotic interactions. For example, cell death of phototrophic microorganisms induced by viral lysis and protistan grazing provides additional resources of dissolved organic matter to the microbial community, and so generates diel cycling in other heterotrophs with different nutrient dependencies. Likewise, differences in the diel transmitting pathway via complex interactions among heterotrophs, and between heterotrophs and their viruses, may also generate higher variation and time lag diel rhythms in different heterotrophic taxa. Thus, sunlight and photosynthesis not only contribute energy and carbon supply, but also directly or indirectly control diel cycling of the microbial community through complex interactions in the photic zone of aquatic ecosystems.

Cyanophage host-derived genes reflect contrasting selective pressures with depth in the oxic and anoxic water column of the Eastern Tropical North Pacific

Cyanophages encode host-derived genes that may increase their fitness. We examined the relative abundance of 18 host-derived cyanophages genes in metagenomes and viromes along depth profiles from the Eastern Tropical North Pacific Oxygen Deficient Zone (ETNP ODZ) where Prochlorococcus dominates a secondary chlorophyll maximum within the ODZ. Cyanophages at the oxic primary chlorophyll maximum encoded genes related to light and phosphate stress (psbA, psbD and pstS in T4-like and psbA in T7-like), but the proportion of cyanophage with these genes decreased with depth. The proportion of cyanophage with purine biosynthesis genes increased with depth in T4-like, but not T7-like cyanophages. No additional host-derived genes were found in deep T7-like cyanophages, suggesting that T4-like and T7-like cyanophages have different host-derived gene acquisition strategies, possibly linked to their different genome packaging mechanisms. In contrast to the ETNP, in the oxic North Atlantic T4-like cyanophages encoded psbA and pstS throughout the euphotic zone. Differences in pstS between the ETNP and the North Atlantic stations were consistent with differences in phosphate concentrations in those regimes. We suggest that the low proportion of cyanophage with psbA within the ODZ reflects the stably stratified low-light conditions occupied by their hosts, a Prochlorococcus ecotype endemic to ODZs.

Rethinking Phage Ecology by Rooting it Within an Established Plant Framework

Despite the abundance and significance of bacteriophages to microbial ecosystems, no broad ecological frameworks exist within which to determine "bacteriophage types" that reflect their ecological strategies and ways in which they interact with bacterial cells. To address this, we repurposed the well-established Grime's triangular CSR framework, which classifies plants according to three axes: competitiveness (C), ability to tolerate stress (S), and capacity to cope with disturbance (R). This framework is distinguished from other accepted schemes, as it seeks to identify individual characteristics of plants to understand their biological strategies and roles within an ecosystem. Our repurposing of the CSR triangle is based on phage transcription and the observation that typically phages have three major distinguishable transcription phases: early, middle, and late. We hypothesize that the proportion of genes expressed in these phases reflects key information about the phage "ecological strategy," namely the C, S, and R strategies, allowing us to examine phages in a similar way to how plants are projected onto the triangle. In the "phage version" of this scheme, we suggest: (1) that some phages prioritize the early phase of transcription that shuts off host defense mechanisms, which reflects competitiveness; (2) other phages prioritize tuning resource management mechanisms in the cell such as nucleotide metabolism during their "mid" expression profile to tolerate stress; and (3) a further subset of phages (termed Ruderals) survive disturbance by investing significant resources into regeneration so they express a higher proportion of their genes during late infection. We examined 42 published phage transcriptomes and show that they fall into discrete CSR categories according to their expression profiles. We discuss these positions in the context of their biology, which is largely consistent with our predictions of specific phage characteristics. In this opinion article, we suggest a starting point to ascribe phages into different functional types and thus understand them in an ecological framework. We suggest that this may have far-reaching implications for the application of phages in therapy and their exploitation to manipulate bacterial communities. We invite further use of this framework via our online tool; www.PhageCSR.ml.

Quantification of T4-Like and T7-Like Cyanophages Using the Polony Method Show They Are Significant Members of the Virioplankton in the North Pacific Subtropical Gyre

The North Pacific Subtropical Gyre (NPSG) is one of the largest biomes on Earth, with the cyanobacterium Prochlorococcus being the most abundant primary producer year-round. Viruses that infect cyanobacteria (cyanophages) influence cyanobacterial mortality, diversity and evolution. Two major cyanophage families are the T4-like cyanomyoviruses and T7-like cyanopodoviruses, yet their abundances and distribution patterns remain unknown due to difficulty in quantifying their populations. To address this limitation, we previously adapted the polony method (for PCR colony) to quantify T7-like cyanophages and applied it to spring populations in the Red Sea. Here, we further adapted the method for the quantification of T4-like cyanophages and analyzed the abundances of T4-like and T7-like cyanophage populations in the photic zone of the NPSG in summer 2015 and spring 2016. Combined, the peak abundances of these two cyanophage families reached 2.8 × 10⁶ and 1.1 × 10⁶ cyanophages ⋅ ml^–1 in the summer and spring, respectively. They constituted between 3 and 16% of total virus-like particles (VLPs), comprising a substantial component of the virioplankton in the NPSG. While both cyanophage families were highly abundant, the T4-like cyanophages were generally 1.3–4.4 fold more so. In summer, cyanophages had similar and reproducible distribution patterns with depth. Abundances were relatively low in the upper mixed layer and increased to form a pronounced subsurface peak at 100 m (1.9 × 10⁶ and 9.1 × 10⁵ phages ⋅ ml^–1 for the T4-like and T7-like cyanophages, respectively), coincident with the maximum in Prochlorococcus populations. Less vertical structure in cyanophage abundances was apparent in the spring profile, despite a subsurface peak in Prochlorococcus numbers. In the summer upper mixed layer, cyanophages constituted a smaller proportion of VLPs than below it and cyanophage to cyanobacteria ratios were considerably lower (1.3–2.8) than those of VLPs to bacteria (8.1–21.2). Differences in abundances between the two families and their contribution to VLPs with depth suggest differences in cyanophage production and/or decay processes relative to other members of the virioplankton in the upper mixed layer. These findings highlight the importance of quantifying distinct populations within the virioplankton to gain accurate understanding of their distribution patterns.

Temporal transcriptional patterns of cyanophage genes suggest synchronized infection of cyanobacteria in the oceans

BACKGROUND

Based on the peak expression times during infection, early, middle, and late genes have been characterized in viruses (cyanophages) that infect the unicellular cyanobacterium Prochlorococcus. Laboratory experiments show that some cyanophages can only replicate in the light and thus exhibit diurnal infection rhythms under light-dark cycles. Field evidence also suggests synchronized infection of Prochlorococcus by cyanophages in the oceans, which should result in progressive expression of cyanophage early, middle, and late genes. However, distinct temporal expression patterns have not been observed in cyanophage field populations.

RESULTS

In this study, we reanalyzed a previous metatranscriptomic dataset collected in the North Pacific Subtropical Gyre. In this dataset, it was previously shown that aggregate transcripts from cyanophage scaffolds display diurnal transcriptional rhythms with transcript abundances decreasing at night. By mapping metatranscriptomic reads to individual viral genes, we identified periodically expressed genes from putative viruses infecting the cyanobacteria Prochlorococcus and Synechococcus, heterotrophic bacteria, and algae. Of the 41 cyanophage genes, 35 were from cyanomyoviruses. We grouped the periodically expressed cyanomyovirus genes into early, middle, and late genes based on the conserved temporal expression patterns of their orthologs in cyanomyovirus laboratory cultures. We found that the peak expression times of late genes in cyanophage field populations were significantly later than those of early and middle genes, which were similar to the temporal expression patterns of synchronized cyanophage laboratory cultures.

CONCLUSIONS

The significantly later peak expression times of late genes in cyanomyovirus field populations suggest that cyanophage infection of Prochlorococcus is synchronized in the North Pacific Subtropical Gyre. The night-time peak expression of late genes also suggests synchronized lysis of Prochlorococcus at night, which might result in synchronized release of dissolved organic matter to the marine food web. Video abstract.

Isolation and Characterization of the Novel Phage JD032 and Global Transcriptomic Response during JD032 Infection of Clostridioides difficile Ribotype 078

C. difficile is one of the most clinically significant intestinal pathogens. Although phages have been shown to effectively control C. difficile infection, the host responses to phage predation have not been fully studied. In this study, we reported the isolation and characterization of a new phage, JD032, and analyzed the global transcriptomic changes in the hypervirulent RT078 C. difficile strain, TW11, during phage JD032 infection. We found that bacterial host mRNA was progressively replaced with phage transcripts, three temporal categories of JD032 gene expression, the extensive interplay between phage-bacterium, antiphage-like responses of the host and phage evasion, and decreased expression of sporulation- and virulence-related genes of the host after phage infection. These findings confirmed the complexity of interactions between C. difficile and phages and suggest that phages undergoing a lytic cycle may also cause different phenotypes in hosts, similar to prophages, which may inspire phage therapy for the control of C. difficile.

Insights into the interaction between phages and their bacterial hosts are crucial for the development of phage therapy. However, only one study has investigated global gene expression of Clostridioides (formerly Clostridium) difficile carrying prophage, and transcriptional reprogramming during lytic infection has not been studied. Here, we presented the isolation, propagation, and characterization of a newly discovered 35,109-bp phage, JD032, and investigated the global transcriptomes of both JD032 and C. difficile ribotype 078 (RT078) strain TW11 during JD032 infection. Transcriptome sequencing (RNA-seq) revealed the progressive replacement of bacterial host mRNA with phage transcripts. The expressed genes of JD032 were clustered into early, middle, and late temporal categories that were functionally similar. Specifically, a gene (JD032_orf016) involved in the lysis-lysogeny decision was identified as an early expression gene. Only 17.7% (668/3,781) of the host genes were differentially expressed, and more genes were downregulated than upregulated. The expression of genes involved in host macromolecular synthesis (DNA/RNA/proteins) was altered by JD032 at the level of transcription. In particular, the expression of the ropA operon was downregulated. Most noteworthy is that the gene expression of some antiphage systems, including CRISPR-Cas, restriction-modification, and toxin-antitoxin systems, was suppressed by JD032 during infection. In addition, bacterial sporulation, adhesion, and virulence factor genes were significantly downregulated. This study provides the first description of the interaction between anaerobic spore-forming bacteria and phages during lytic infection and highlights new aspects of C. difficile phage-host interactions.

IMPORTANCEC. difficile is one of the most clinically significant intestinal pathogens. Although phages have been shown to effectively control C. difficile infection, the host responses to phage predation have not been fully studied. In this study, we reported the isolation and characterization of a new phage, JD032, and analyzed the global transcriptomic changes in the hypervirulent RT078 C. difficile strain, TW11, during phage JD032 infection. We found that bacterial host mRNA was progressively replaced with phage transcripts, three temporal categories of JD032 gene expression, the extensive interplay between phage-bacterium, antiphage-like responses of the host and phage evasion, and decreased expression of sporulation- and virulence-related genes of the host after phage infection. These findings confirmed the complexity of interactions between C. difficile and phages and suggest that phages undergoing a lytic cycle may also cause different phenotypes in hosts, similar to prophages, which may inspire phage therapy for the control of C. difficile.

Phage-specific metabolic reprogramming of virocells

Ocean viruses are abundant and infect 20–40% of surface microbes. Infected cells, termed virocells, are thus a predominant microbial state. Yet, virocells and their ecosystem impacts are understudied, thus precluding their incorporation into ecosystem models. Here we investigated how unrelated bacterial viruses (phages) reprogram one host into contrasting virocells with different potential ecosystem footprints. We independently infected the marine Pseudoalteromonas bacterium with siphovirus PSA-HS2 and podovirus PSA-HP1. Time-resolved multi-omics unveiled drastically different metabolic reprogramming and resource requirements by each virocell, which were related to phage–host genomic complementarity and viral fitness. Namely, HS2 was more complementary to the host in nucleotides and amino acids, and fitter during infection than HP1. Functionally, HS2 virocells hardly differed from uninfected cells, with minimal host metabolism impacts. HS2 virocells repressed energy-consuming metabolisms, including motility and translation. Contrastingly, HP1 virocells substantially differed from uninfected cells. They repressed host transcription, responded to infection continuously, and drastically reprogrammed resource acquisition, central carbon and energy metabolisms. Ecologically, this work suggests that one cell, infected versus uninfected, can have immensely different metabolisms that affect the ecosystem differently. Finally, we relate phage–host genome complementarity, virocell metabolic reprogramming, and viral fitness in a conceptual model to guide incorporating viruses into ecosystem models.

Metabolic and biogeochemical consequences of viral infection in aquatic ecosystems

Ecosystems are controlled by 'bottom-up' (resources) and 'top-down' (predation) forces. Viral infection is now recognized as a ubiquitous top-down control of microbial growth across ecosystems but, at the same time, cell death by viral predation influences, and is influenced by, resource availability. In this Review, we discuss recent advances in understanding the biogeochemical impact of viruses, focusing on how metabolic reprogramming of host cells during lytic viral infection alters the flow of energy and nutrients in aquatic ecosystems. Our synthesis revealed several emerging themes. First, viral infection transforms host metabolism, in part through virus-encoded metabolic genes; the functions performed by these genes appear to alleviate energetic and biosynthetic bottlenecks to viral production. Second, viral infection depends on the physiological state of the host cell and on environmental conditions, which are challenging to replicate in the laboratory. Last, metabolic reprogramming of infected cells and viral lysis alter nutrient cycling and carbon export in the oceans, although the net impacts remain uncertain. This Review highlights the need for understanding viral infection dynamics in realistic physiological and environmental contexts to better predict their biogeochemical consequences.

Cooccurrence of Broad- and Narrow-Host-Range Viruses Infecting the Bloom-Forming Toxic Cyanobacterium Microcystis aeruginosa

Viruses play important roles in regulating the abundance and composition of bacterial populations in aquatic ecosystems. The bloom-forming toxic cyanobacterium Microcystis aeruginosa is predicted to interact with diverse cyanoviruses, resulting in Microcystis population diversification. However, current knowledge of the genomes from these viruses and their infection programs is limited to those of Microcystis virus Ma-LMM01. Here, we performed a time series sampling at a small pond in Japan during a Microcystis bloom and then investigated the genomic information and transcriptional dynamics of Microcystis-interacting viruses using metagenomic and metatranscriptomic approaches. We identified 15 viral genomic fragments classified into three groups, groups I (including Ma-LMM01), II (high abundance and transcriptional activity), and III (new lineages). According to the phylogenetic distribution of Microcystis strains possessing spacers against each viral group, the group II-original viruses interacted with all three phylogenetically distinct Microcystis population types (phylotypes), whereas the groups I and III-original viruses interacted with only one or two phylotypes, indicating the cooccurrence of broad- (group II) and narrow (groups I and III)-host-range viruses in the bloom. These viral fragments showed the highest transcriptional levels during daytime regardless of their genomic differences. Interestingly, M. aeruginosa expressed antiviral defense genes in the environment, unlike what was seen with an Ma-LMM01 infection in a previous culture experiment. Given that broad-host-range viruses often induce antiviral responses within alternative hosts, our findings suggest that such antiviral responses might inhibit viral multiplication, mainly that of broad-host-range viruses like those in group II.IMPORTANCE The bloom-forming toxic cyanobacterium Microcystis aeruginosa is thought to have diversified its population through the interactions between host and viruses in antiviral defense systems. However, current knowledge of viral genomes and infection programs is limited to those of Microcystis virus Ma-LMM01, which was a narrow host range in which it can escape from the highly abundant host defense systems. Our metagenomic approaches unveiled the cooccurrence of narrow- and broad-host-range Microcystis viruses, which included fifteen viral genomic fragments from Microcystis blooms that were classified into three groups. Interestingly, Microcystis antiviral defense genes were expressed against viral infection in the environment, unlike what was seen in a culture experiment with Ma-LMM01. Given that viruses with a broad host range often induce antiviral responses within alternative hosts, our findings suggest that antiviral responses inhibit viral reproduction, especially that of broad-range viruses like those in group II. This paper augments our understanding of the interactions between M. aeruginosa and its viruses and fills an important knowledge gap.

Resistance in marine cyanobacteria differs against specialist and generalist cyanophages

Long-term coexistence between unicellular cyanobacteria and their lytic viruses (cyanophages) in the oceans is thought to be due to the presence of sensitive cells in which cyanophages reproduce, ultimately killing the cell, while other cyanobacteria survive due to resistance to infection. Here, we investigated resistance in marine cyanobacteria from the genera Synechococcus and Prochlorococcus and compared modes of resistance against specialist and generalist cyanophages belonging to the T7-like and T4-like cyanophage families. Resistance was extracellular in most interactions against specialist cyanophages irrespective of the phage family, preventing entry into the cell. In contrast, resistance was intracellular in practically all interactions against generalist T4-like cyanophages. The stage of intracellular arrest was interaction-specific, halting at various stages of the infection cycle. Incomplete infection cycles proceeded to various degrees of phage genome transcription and translation as well as phage genome replication in numerous interactions. In a particularly intriguing case, intracellular capsid assembly was observed, but the phage genome was not packaged. The cyanobacteria survived the encounter despite late-stage infection and partial genome degradation. We hypothesize that this is tolerated due to genome polyploidy, which we found for certain strains of both Synechococcus and Prochlorococcus Our findings unveil a heavy cost of promiscuous entry of generalist phages into nonhost cells that is rarely paid by specialist phages and suggests the presence of unknown mechanisms of intracellular resistance in the marine unicellular cyanobacteria. Furthermore, these findings indicate that the range for virus-mediated horizontal gene transfer extends beyond hosts to nonhost cyanobacterial cells.

Nitrogen sourcing during viral infection of marine cyanobacteria

The building blocks of a virus derived from de novo biosynthesis during infection and/or catabolism of preexisting host cell biomass, and the relative contribution of these 2 sources has important consequences for understanding viral biogeochemistry. We determined the uptake of extracellular nitrogen (N) and its biosynthetic incorporation into both virus and host proteins using an isotope-labeling proteomics approach in a model marine cyanobacterium Synechococcus WH8102 infected by a lytic cyanophage S-SM1. By supplying dissolved N as ¹⁵N postinfection, we found that proteins in progeny phage particles were composed of up to 41% extracellularly derived N, while proteins of the infected host cell showed almost no isotope incorporation, demonstrating that de novo amino acid synthesis continues during infection and contributes specifically and substantially to phage replication. The source of N for phage protein synthesis shifted over the course of infection from mostly host derived in the early stages to more medium derived later on. We show that the photosystem II reaction center proteins D1 and D2, which are auxiliary metabolic genes (AMGs) in the S-SM1 genome, are made de novo during infection in an apparently light-dependent manner. We also identified a small set of host proteins that continue to be produced during infection; the majority are homologs of AMGs in S-SM1 or other viruses, suggesting selective continuation of host protein production during infection. The continued acquisition of nutrients by the infected cell and their utilization for phage replication are significant for both evolution and biogeochemical impact of viruses.

Cyanophage MazG is a pyrophosphohydrolase but unable to hydrolyse magic spot nucleotides

Bacteriophage possess a variety of auxiliary metabolic genes of bacterial origin. These proteins enable them to maximize infection efficiency, subverting bacterial metabolic processes for the purpose of viral genome replication and synthesis of the next generation of virion progeny. Here, we examined the enzymatic activity of a cyanophage MazG protein - a putative pyrophosphohydrolase previously implicated in regulation of the stringent response via reducing levels of the central alarmone molecule (p)ppGpp. We demonstrate, however, that the purified viral MazG shows no binding or hydrolysis activity against (p)ppGpp. Instead, dGTP and dCTP appear to be the preferred substrates of this protein, consistent with a role preferentially hydrolysing deoxyribonucleotides from the high GC content host Synechococcus genome. This showcases a new example of the fine-tuned nature of viral metabolic processes.

Host-hijacking and planktonic piracy: how phages command the microbial high seas

Microbial communities living in the oceans are major drivers of global biogeochemical cycles. With nutrients limited across vast swathes of the ocean, marine microbes eke out a living under constant assault from predatory viruses. Viral concentrations exceed those of their bacterial prey by an order of magnitude in surface water, making these obligate parasites the most abundant biological entities in the ocean. Like the pirates of the 17th and 18th centuries that hounded ships plying major trade and exploration routes, viruses have evolved mechanisms to hijack microbial cells and repurpose their cargo and indeed the vessels themselves to maximise viral propagation. Phenotypic reconfiguration of the host is often achieved through Auxiliary Metabolic Genes - genes originally derived from host genomes but maintained and adapted in viral genomes to redirect energy and substrates towards viral synthesis. In this review, we critically evaluate the literature describing the mechanisms used by bacteriophages to reconfigure host metabolism and to plunder intracellular resources to optimise viral production. We also highlight the mechanisms used when, in challenging environments, a 'batten down the hatches' strategy supersedes that of 'plunder and pillage'. Here, the infecting virus increases host fitness through phenotypic augmentation in order to ride out the metaphorical storm, with a concomitant impact on host substrate uptake and metabolism, and ultimately, their interactions with their wider microbial community. Thus, the traditional view of the virus-host relationship as predator and prey does not fully characterise the variety or significance of the interactions observed. Recent advances in viral metagenomics have provided a tantalising glimpse of novel mechanisms of viral metabolic reprogramming in global oceans. Incorporation of these new findings into global biogeochemical models requires experimental evidence from model systems and major improvements in our ability to accurately predict protein function from sequence data.