Phylogenomics of T4 cyanophages: lateral gene transfer in the 'core' and origins of host genes

Discarded diversity: Novel megaphages, auxiliary metabolic genes, and virally encoded CRISPR-Cas systems in landfills

Background

Viruses are the most abundant microbial entity on the planet, impacting microbial community structure and ecosystem services. Despite outnumbering Bacteria and Archaea by an order of magnitude, viruses have been comparatively underrepresented in reference databases. Metagenomic examinations have illustrated that viruses of Bacteria and Archaea have been specifically understudied in engineered environments. Here we employed metagenomic and computational biology methods to examine the diversity, host interactions, and genetic systems of viruses predicted from 27 samples taken from three municipal landfills across North America.

Results

We identified numerous viruses that are not represented in reference databases, including the third largest bacteriophage genome identified to date (~678 kbp), and note a cosmopolitan diversity of viruses in landfills that are distinct from viromes in other systems. Host-virus interactions were examined via host CRISPR spacer to viral protospacer mapping which captured hyper-targeted viral populations and six viral populations predicted to infect across multiple phyla. Virally-encoded auxiliary metabolic genes (AMGs) were identified with the potential to augment hosts' methane, sulfur, and contaminant degradation metabolisms, including AMGs not previously reported in literature. CRISPR arrays and CRISPR-Cas systems were identified from predicted viral genomes, including the two largest bacteriophage genomes to contain these genetic features. Some virally encoded Cas effector proteins appear distinct relative to previously reported Cas systems and are interesting targets for potential genome editing tools.

Conclusions

Our observations indicate landfills, as heterogeneous contaminated sites with unique selective pressures, are key locations for diverse viruses and atypical virus-host dynamics.

A novel long-tailed myovirus represents a new T4-like cyanophage cluster

Cyanophages affect the abundance, diversity, metabolism, and evolution of picocyanobacteria in marine ecosystems. Here we report an estuarine Synechococcus phage, S-CREM2, which represents a novel viral genus and leads to the establishment of a new T4-like cyanophage clade named cluster C. S-CREM2 possesses the longest tail (~418 nm) among isolated cyanomyoviruses and encodes six tail-related proteins that are exclusively homologous to those predicted in the cluster C cyanophages. Furthermore, S-CREM2 may carry three regulatory proteins in the virion, which may play a crucial role in optimizing the host intracellular environment for viral replication at the initial stage of infection. The cluster C cyanophages lack auxiliary metabolic genes (AMGs) that are commonly found in cyanophages of the T4-like clusters A and B and encode unique AMGs like an S-type phycobilin lyase gene. A variation in the composition of tRNA and cis-regulatory RNA genes was observed between the marine and freshwater phage strains in cluster C, reflecting their different modes of coping with hosts and habitats. The cluster C cyanophages are widespread in estuarine and coastal regions and exhibit equivalent or even higher relative abundance compared to those of clusters A and B cyanophages in certain estuarine regions. The isolation of cyanophage S-CREM2 provides new insights into the phage-host interactions mediated by both newly discovered AMGs and virion-associated proteins and emphasizes the ecological significance of cluster C cyanophages in estuarine environments.

Cryo-EM structure of cyanophage P-SCSP1u offers insights into DNA gating and evolution of T7-like viruses

Cyanophages, together with their host cyanobacteria, play important roles in marine biogeochemical cycles and control of marine food webs. The recently identified MPP-C (Marine Picocyanobacteria Podovirus clade C) cyanophages, belonging to the T7-like podoviruses, contain the smallest genomes among cyanopodoviruses and exhibit distinct infection kinetics. However, understanding of the MPP-C cyanophage infection process is hindered by the lack of high-resolution structural information. Here, we report the cryo-EM structure of the cyanophage P-SCSP1u, a representative member of the MPP-C phages, in its native form at near-atomic resolution, which reveals the assembly mechanism of the capsid and molecular interaction of the portal-tail complex. Structural comparison of the capsid proteins of P-SCSP1u and other podoviruses with known structures provides insights into the evolution of T7-like viruses. Furthermore, our study provides the near-atomic resolution structure of portal-tail complex for T7-like viruses. On the basis of previously reported structures of phage T7, we identify an additional valve and gate to explain the DNA gating mechanism for the T7-like viruses.

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.

The Impact of Viral Infection on the Chemistries of the Earth's Most Abundant Photosynthesizes: Metabolically Talented Aquatic Cyanobacteria

Cyanobacteria are the most abundant photosynthesizers on earth, and as such, they play a central role in marine metabolite generation, ocean nutrient cycling, and the control of planetary oxygen generation. Cyanobacteriophage infection exerts control on all of these critical processes of the planet, with the phage-ported homologs of genes linked to photosynthesis, catabolism, and secondary metabolism (marine metabolite generation). Here, we analyze the 153 fully sequenced cyanophages from the National Center for Biotechnology Information (NCBI) database and the 45 auxiliary metabolic genes (AMGs) that they deliver into their hosts. Most of these AMGs are homologs of those found within cyanobacteria and play a key role in cyanobacterial metabolism-encoding proteins involved in photosynthesis, central carbon metabolism, phosphate metabolism, methylation, and cellular regulation. A greater understanding of cyanobacteriophage infection will pave the way to a better understanding of carbon fixation and nutrient cycling, as well as provide new tools for synthetic biology and alternative approaches for the use of cyanobacteria in biotechnology and sustainable manufacturing.

Characterization and genomic analysis of a novel Synechococcus phage S-H9-2 belonging to Bristolvirus genus isolated from the Yellow Sea

Cyanophages are known to influence the population dynamics and community structure of cyanobacteria and thus play an important role in biogeochemical cycles in aquatic ecosystems. In this study, a novel Synechococcus phage S-H9-2 infecting Synechococcus sp. WH 8102 was isolated from the coastal water of the Yellow Sea. Synechococcus phage S-H9-2 contains a 187,320 bp genome of double-stranded DNA with a G + C content of 40.3%, 202 potential open reading frames (ORFs), and 15 tRNAs. Phylogenetic analysis and nucleotide-based intergenomic similarity suggest that Synechococcus phage S-H9-2 belongs to the Bristolvirus genus under the family Kyanoviridae. Homologs of the S-H9-2 open reading frame can be found in a variety of marine environments, as shown by the results of mapping the genome sequence of S-H9-2 to the Global Ocean Viromes 2.0 dataset. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as phylogenetic relationships based on complete genome sequences, reflect the mechanism of phage-host interaction and host-specific strategies for adaptation to environmental conditions. This study enriches the current genomic database of cyanophage and contributed to our understanding of the virus-host interactions and their adaption to the environment.

CRISPR-resolved virus-host interactions in a municipal landfill include non-specific viruses, hyper-targeted viral populations, and interviral conflicts

Viruses are the most abundant microbial guild on the planet, impacting microbial community structure and ecosystem services. Viruses are specifically understudied in engineered environments, including examinations of their host interactions. We examined host-virus interactions via host CRISPR spacer to viral protospacer mapping in a municipal landfill across two years. Viruses comprised ~ 4% of both the unassembled reads and assembled basepairs. A total of 458 unique virus-host connections captured hyper-targeted viral populations and host CRISPR array adaptation over time. Four viruses were predicted to infect across multiple phyla, suggesting that some viruses are far less host-specific than is currently understood. We detected 161 viral elements that encode CRISPR arrays, including one with 187 spacers, the longest virally-encoded CRISPR array described to date. Virally-encoded CRISPR arrays targeted other viral elements in interviral conflicts. CRISPR-encoding proviruses integrated into host chromosomes were latent examples of CRISPR-immunity-based superinfection exclusion. The bulk of the observed virus-host interactions fit the one-virus-one-host paradigm, but with limited geographic specificity. Our networks highlight rare and previously undescribed complex interactions influencing the ecology of this dynamic engineered system. Our observations indicate landfills, as heterogeneous contaminated sites with unique selective pressures, are key locations for atypical virus-host dynamics.

Spatial and temporal metagenomics of river compartments reveals viral community dynamics in an urban impacted stream

Although river ecosystems comprise less than 1% of Earth's total non-glaciated area, they are critical modulators of microbially and virally orchestrated global biogeochemical cycles. However, most studies either use data that is not spatially resolved or is collected at timepoints that do not reflect the short life cycles of microorganisms. As a result, the relevance of microbiome interactions and the impacts they have over time on biogeochemical cycles are poorly understood. To assess how viral and microbial communities change over time, we sampled surface water and pore water compartments of the wastewater-impacted River Erpe in Germany every 3 hours over a 48-hour period resulting in 32 metagenomes paired to geochemical and metabolite measurements. We reconstructed 6,500 viral and 1,033 microbial genomes and found distinct communities associated with each river compartment. We show that 17% of our vMAGs clustered to viruses from other ecosystems like wastewater treatment plants and rivers. Our results also indicated that 70% of the viral community was persistent in surface waters, whereas only 13% were persistent in the pore waters taken from the hyporheic zone. Finally, we predicted linkages between 73 viral genomes and 38 microbial genomes. These putatively linked hosts included members of the Competibacteraceae, which we suggest are potential contributors to carbon and nitrogen cycling. Together, these findings demonstrate that microbial and viral communities in surface waters of this urban river can exist as stable communities along a flowing river; and raise important considerations for ecosystem models attempting to constrain dynamics of river biogeochemical cycles.

Spatiotemporal Shift of T4-Like Phage Community Structure in the Three Largest Estuaries of China

Estuaries are one of the most highly productive and economically important ecosystems at the continent-ocean interface. Estuary productivity is largely determined by the microbial community structure and activity. Viruses are major agents of microbial mortality and are key drivers of global geochemical cycles. However, the taxonomic diversity of viral communities and their spatial-temporal distribution in estuarine ecosystems have been poorly studied. In this study, we investigated the T4-like viral community composition at three major Chinese estuaries in winter and in summer. Diverse T4-like viruses, which were divided into three main clusters (Clusters I to III), were revealed. The Marine Group of Cluster III, with seven identified subgroups, was the most dominant (averaging 76.5% of the total sequences) in the Chinese estuarine ecosystems. Significant variations of T4-like viral community composition were observed among estuaries and seasons, with higher diversity occurring in winter. Among various environmental variables, temperature was a main driver of the viral communities. This study demonstrates viral assemblage diversification and seasonality in Chinese estuarine ecosystems. IMPORTANCE Viruses are ubiquitous but largely uncharacterized members of aquatic environments that cause significant mortality in microbial communities. Recent large-scale oceanic projects have greatly advanced our understanding of viral ecology in marine environments, but those studies mostly focused on oceanic regions. There have yet to be spatiotemporal studies of viral communities in estuarine ecosystems, which are unique habitats that play a significant role in global ecology and biogeochemistry. This work is the first comprehensive study that provides a detailed picture of the spatial and seasonal variation of viral communities (specifically, T4-like viral communities) in three major estuarine ecosystems in China. These findings provide much-needed knowledge regarding estuarine viral ecosystems, which currently lags in oceanic ecosystem research.

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.

Two Novel Yersinia pestis Bacteriophages with a Broad Host Range: Potential as Biocontrol Agents in Plague Natural Foci

Bacteriophages (phages) have been successfully used as disinfectors to kill bacteria in food and the environment and have been used medically for curing human diseases. The objective of this research was to elucidate the morphological and genomic characteristics of two novel Yersinia pestis phages, vB_YpeM_ MHS112 (MHS112) and vB_YpeM_GMS130 (GMS130), belonging to the genus Gaprivervirus, subfamily Tevenvirinae, family Myoviridae. Genome sequencing showed that the sizes of MHS112 and GMS130 were 170507 and 168552 bp, respectively. A total of 303 and 292 open reading frames with 2 tRNA and 3 tRNA were predicted in MHS112 and GMS130, respectively. The phylogenetic relationships were analysed among the two novel Y. pestis phages, phages in the genus Gaprivervirus, and several T4-like phages infecting the Yersinia genus. The bacteriophage MHS112 and GMS130 exhibited a wider lytic host spectrum and exhibited comparative temperature and pH stability. Such features signify that these phages do not need to rely on Y. pestis as their host bacteria in the ecological environment, while they could be based on more massive Enterobacteriales species to propagate and form ecological barriers against Y. pestis pathogens colonised in plague foci. Such characteristics indicated that the two phages have potential as biocontrol agents for eliminating the endemics of animal plague in natural plague foci.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml^-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml^-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

A Unique Set of Auxiliary Metabolic Genes Found in an Isolated Cyanophage Sheds New Light on Marine Phage-Host Interactions

Cyanophages, viruses that infect cyanobacteria, are abundant and widely distributed in aquatic ecosystems, playing important roles in regulating the abundance, activity, diversity, and evolution of cyanobacteria. A T4-like cyanophage, S-SCSM1, infecting Synechococcus and Prochlorococcus strains of different ecotypes, was isolated from the South China Sea in this study. For the first time, a mannose-6-phosphate isomerase (MPI) gene was identified in the cultured cyanophage. At least 11 phylogenetic clusters of cyanophage MPIs were retrieved and identified from the marine metagenomic data sets, indicating that cyanophage MPIs in the marine environment are extremely diverse. The existence of 24 genes encoding 2-oxoglutarate (2OG)-Fe(II) oxygenase superfamily proteins in the S-SCSM1 genome emphasizes their potential importance and diverse functions in reprogramming host metabolism during phage infection. Novel cell wall synthesis and modification genes found in the S-SCSM1 genome indicate that diverse phenotypic modifications imposed by phages on cyanobacterial hosts remain to be discovered. Two noncoding RNAs of cis-regulatory elements in the S-SCSM1 genome were predicted to be associated with host exopolysaccharide metabolism and photosynthesis. The isolation and genomic characterization of cyanophage S-SCSM1 provide more information on the genetic diversity of cyanophages and phage-host interactions in the marine environment. IMPORTANCE Cyanophages play important ecological roles in aquatic ecosystems. Genomic and proteomic characterizations of the T4-like cyanophage S-SCSM1 indicate that novel and diverse viral genes and phage-host interactions in the marine environment remain unexplored. The first identified mannose-6-phosphate isomerase (MPI) gene from a cultured cyanophage was found in the S-SCSM1 genome, although MPIs were previously found in viral metagenomes at high frequencies similar to those of the cyanophage photosynthetic gene psbA. The presence of 24 genes encoding 2-oxoglutarate (2OG)-Fe(II) oxygenase superfamily proteins, novel cell wall synthesis and modification genes, a nonbleaching protein A gene, and 2 noncoding RNAs of cis-regulatory elements in the S-SCSM1 genome as well as the presence of a virion-associated regulatory protein indicate the diverse functions that cyanophages have in reprogramming the metabolism and modifying the phenotypes of hosts during infection.

Genomes from Uncultivated Pelagiphages Reveal Multiple Phylogenetic Clades Exhibiting Extensive Auxiliary Metabolic Genes and Cross-Family Multigene Transfers

For the abundant marine Alphaproteobacterium Pelagibacter (SAR11), and other bacteria, phages are powerful forces of mortality. However, little is known about the most abundant Pelagiphages in nature, such as the widespread HTVC023P-type, which is currently represented by two cultured phages. Using viral metagenomic data sets and fluorescence-activated cell sorting, we recovered 80 complete, undescribed Podoviridae genomes that form 10 phylogenomically distinct clades (herein, named Clades I to X) related to the HTVC023P-type. These expanded the HTVC023P-type pan-genome by 15-fold and revealed 41 previously unknown auxiliary metabolic genes (AMGs) in this viral lineage. Numerous instances of partner-AMGs (colocated and involved in related functions) were observed, including partners in nucleotide metabolism, DNA hypermodification, and Curli biogenesis. The Type VIII secretion system (T8SS) responsible for Curli biogenesis was identified in nine genomes and expanded the repertoire of T8SS proteins reported thus far in viruses. Additionally, the identified T8SS gene cluster contained an iron-dependent regulator (FecR), as well as a histidine kinase and adenylate cyclase that can be implicated in T8SS function but are not within T8SS operons in bacteria. While T8SS are lacking in known Pelagibacter, they contribute to aggregation and biofilm formation in other bacteria. Phylogenetic reconstructions of partner-AMGs indicate derivation from cellular lineages with a more recent transfer between viral families. For example, homologs of all T8SS genes are present in syntenic regions of distant Myoviridae Pelagiphages, and they appear to have alphaproteobacterial origins with a later transfer between viral families. The results point to an unprecedented multipartner-AMG transfer between marine Myoviridae and Podoviridae. Together with the expansion of known metabolic functions, our studies provide new prospects for understanding the ecology and evolution of marine phages and their hosts. IMPORTANCE One of the most abundant and diverse marine bacterial groups is Pelagibacter. Phages have roles in shaping Pelagibacter ecology; however, several Pelagiphage lineages are represented by only a few genomes. This paucity of data from even the most widespread lineages has imposed limits on the understanding of the diversity of Pelagiphages and their impacts on hosts. Here, we report 80 complete genomes, assembled directly from environmental data, which are from undescribed Pelagiphages and render new insights into the manipulation of host metabolism during infection. Notably, the viruses have functionally related partner genes that appear to be transferred between distant viruses, including a suite that encode a secretion system which both brings a new functional capability to the host and is abundant in phages across the ocean. Together, these functions have important implications for phage evolution and for how Pelagiphage infection influences host biology in manners extending beyond canonical viral lysis and mortality.

A Novel Freshwater Cyanophage Mae-Yong1326-1 Infecting Bloom-Forming Cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is a major harmful cyanobacterium causing water bloom worldwide. Cyanophage has been proposed as a promising tool for cyanobacterial bloom. In this study, M. aeruginosa FACHB-1326 was used as an indicator host to isolate cyanophage from Lake Taihu. The isolated Microcystis cyanophage Mae-Yong1326-1 has an elliptical head of about 47 nm in diameter and a slender flexible tail of about 340 nm in length. Mae-Yong1326-1 could lyse cyanobacterial strains across three orders (Chroococcales, Nostocales, and Oscillatoriales) in the host range experiments. Mae-Yong1326-1 was stable in stability tests, maintaining high titers at 0–40 °C and at a wide pH range of 3–12. Mae-Yong 1326-1 has a burst size of 329 PFU/cell, which is much larger than the reported Microcystis cyanophages so far. The complete genome of Mae-Yong1326-1 is a double-stranded DNA of 48, 822 bp, with a G + C content of 71.80% and long direct terminal repeats (DTR) of 366 bp, containing 57 predicted ORFs. No Mae-Yong1326-1 ORF was found to be associated with virulence factor or antibiotic resistance. PASC scanning illustrated that the highest nucleotide sequence similarity between Mae-Yong1326-1 and all known phages in databases was only 17.75%, less than 70% (the threshold to define a genus), which indicates that Mae-Yong1326-1 belongs to an unknown new genus. In the proteomic tree based on genome-wide sequence similarities, Mae-Yong1326-1 distantly clusters with three unclassified Microcystis cyanophages (MinS1, Mwe-Yong1112-1, and Mwes-Yong2). These four Microcystis cyanophages form a monophyletic clade, which separates at a node from the other clade formed by two independent families (Zierdtviridae and Orlajensenviridae) of Caudoviricetes class. We propose to establish a new family to harbor the Microcystis cyanophages Mae-Yong1326-1, MinS1, Mwe-Yong1112-1, and Mwes-Yong2. This study enriched the understanding of freshwater cyanophages.

A Novel Alteromonas Phage Lineage with a Broad Host Range and Small Burst Size

Alteromonas is an opportunistic marine bacterium that persists in the global ocean and has important ecological significance. However, current knowledge about the diversity and ecology of alterophages (phages that infect Alteromonas) is lacking. Here, three similar phages infecting Alteromonas macleodii ATCC 27126^T were isolated and physiologically characterized. Transmission electron microscopy revealed Siphoviridae morphology, with an oblate icosahedral head and a long noncontractile tail. Notably, these members displayed a small burst size (15–19 plaque-forming units/cell) yet an extensively broad host spectrum when tested on 175 Alteromonas strains. Such unique infection kinetics are potentially associated with discrepancies in codon usage bias from the host tRNA inventory. Phylogenetic analysis indicated that the three phages are closely evolutionarily related; they clustered at the species level and represent a novel genus. Three auxiliary metabolic genes with roles in nucleotide metabolism and putative biofilm dispersal were found in these phage genomes, which revealed important biogeochemical significance of these alterophages in marine ecosystems. Our isolation and characterization of these novel phages expand the current understanding of alterophage diversity, evolution, and phage–host interactions.

IMPORTANCE The marine bacterium Alteromonas is prevalent in the global ocean with crucial ecological significance; however, little is known about the diversity and evolution of its bacteriophages that profoundly affect the bacterial communities. Our study characterized a novel genus of three newly isolated Alteromonas phages that exhibited a distinct infection strategy of broad host spectrum and small burst size. This strategy is likely a consequence of the viral trade-off between virulence and lysis profiles during phage–host coevolution, and our work provides new insight into viral evolution and infection strategies.

Auxiliary Metabolic Gene Functions in Pelagic and Benthic Viruses of the Baltic Sea

Marine microbial communities are facing various ecosystem fluctuations (e.g., temperature, organic matter concentration, salinity, or redox regimes) and thus have to be highly adaptive. This might be supported by the acquisition of auxiliary metabolic genes (AMGs) originating from virus infections. Marine bacteriophages frequently contain AMGs, which allow them to augment their host’s metabolism or enhance virus fitness. These genes encode proteins for the same metabolic functions as their highly similar host homologs. In the present study, we analyzed the diversity, distribution, and composition of marine viruses, focusing on AMGs to identify their putative ecologic role. We analyzed viruses and assemblies of 212 publicly available metagenomes obtained from sediment and water samples across the Baltic Sea. In general, the virus composition in both compartments differed compositionally. While the predominant viral lifestyle was found to be lytic, lysogeny was more prevalent in sediments than in the pelagic samples. The highest proportion of AMGs was identified in the genomes of Myoviridae. Overall, the most abundantly occurring AMGs are encoded for functions that protect viruses from degradation by their hosts, such as methylases. Additionally, some detected AMGs are known to be involved in photosynthesis, 7-cyano-7-deazaguanine synthesis, and cobalamin biosynthesis among other functions. Several AMGs that were identified in this study were previously detected in a large-scale analysis including metagenomes from various origins, i.e., different marine sites, wastewater, and the human gut. This supports the theory of globally conserved core AMGs that are spread over virus genomes, regardless of host or environment.

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.

Viral community analysis in a marine oxygen minimum zone indicates increased potential for viral manipulation of microbial physiological state

Microbial communities in oxygen minimum zones (OMZs) are known to have significant impacts on global biogeochemical cycles, but viral influence on microbial processes in these regions are much less studied. Here we provide baseline ecological patterns using microscopy and viral metagenomics from the Eastern Tropical North Pacific (ETNP) OMZ region that enhance our understanding of viruses in these climate-critical systems. While extracellular viral abundance decreased below the oxycline, viral diversity and lytic infection frequency remained high within the OMZ, demonstrating that viral influences on microbial communities were still substantial without the detectable presence of oxygen. Viral community composition was strongly related to oxygen concentration, with viral populations in low-oxygen portions of the water column being distinct from their surface layer counterparts. However, this divergence was not accompanied by the expected differences in viral-encoded auxiliary metabolic genes (AMGs) relating to nitrogen and sulfur metabolisms that are known to be performed by microbial communities in these low-oxygen and anoxic regions. Instead, several abundant AMGs were identified in the oxycline and OMZ that may modulate host responses to low-oxygen stress. We hypothesize that this is due to selection for viral-encoded genes that influence host survivability rather than modulating host metabolic reactions within the ETNP OMZ. Together, this study shows that viruses are not only diverse throughout the water column in the ETNP, including the OMZ, but their infection of microorganisms has the potential to alter host physiological state within these biogeochemically important regions of the ocean.

Record-Breaking Rain Event Altered Estuarine Viral Assemblages

Viruses are the dominant biological entity in the ocean, play a vital role in biogeochemical cycles, and provide their hosts with novel metabolic capabilities through auxiliary metabolic genes (AMGs). Hurricane Harvey was a category 4 hurricane that made landfall on the Texas coast in 2017 and lashed the Houston area with 1.4–1.7 × 1010 m3 of rainfall. In this paper, we aim to characterize how the changes in abiotic conditions brought by Hurricane Harvey altered the viral assemblages of Galveston Bay at the taxonomic level and determine how viral ecosystem functions were altered. Metagenomes of the viruses and their hosts were sequenced from a transect in Galveston Bay over the five weeks following the storm. Our results show that the viral assemblages of Galveston Bay dramatically changed following Hurricane Harvey’s landfall. Of the abiotic parameters measured, salinity had the strongest effect on shaping the viral assemblages. In the five weeks following Hurricane Harvey, there was a steady increase of metabolic genes and putative viral infections. Our study provides the first in-depth look at how marine viral assemblages respond and recover from extreme rainfall events, which models predict will become more frequent and intense with climate change.

phoH-carrying virus communities responded to multiple factors and their correlation network with prokaryotes in sediments along Bohai Sea, Yellow Sea, and East China Sea in China

Viruses carrying phoH genes are an important functional group that may boost phosphate metabolism of their prokaryote hosts and affect phosphorus cycle in the ocean. However, at present, very little is known about the phoH-carrying viruses' community structure and diversity in marine sediments, as well as their correlation network with prokaryotes and environment. Here, via a large spatial scale investigation along the Bohai Sea, Yellow Sea, and East China Sea, for the first time, diverse unknown benthic phoH-carrying viruses were uncovered, which were mainly affiliated to three clusters. Interestingly, these viruses presented a very distinct community structure compared to those in seawaters. Correlation network analysis implied that these viruses might mainly infect the prokaryotes of Gamm-/Delta-proteobacteria, Thaumarchaeota, and Cyanobacteria in sediments. Distinct virus-prokaryote correlation network modules were shown in different sea areas. These modules' highly nested feature implied their coevolution with prokaryotes during long-term arms race. Their distribution in sediments was influenced by multiple factors including geographic separation and the key environmental variables of total organic carbon and total phosphorus, and responded to terrestrial inputs and coastal aquaculture activities. The results of this study provide novel insights into the benthic virus communities potentially participating in phosphorus cycling in the ocean.

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.

Identification of the tail assembly chaperone genes of T4-Like phages suggests a mechanism other than translational frameshifting for biogenesis of their encoded proteins

Tape measure (TM) proteins are essential for the formation of long-tailed phages. TM protein assembly into tails requires the action of tail assembly chaperones (TACs). TACs (e.g. gpG and gpT of E. coli phage lambda) are usually produced in a short (TAC-N) and long form (TAC-NC) with the latter comprised of TAC-N with an additional C-terminal domain (TAC-C). TAC-NC is generally synthesized through a ribosomal frameshifting mechanism. TAC encoding genes have never been identified in the intensively studied Escherichia coli phage T4, or any related phages. Here, we have bioinformatically identified putative TAC encoding genes in diverse T4-like phage genomes. The frameshifting mechanism for producing TAC-NC appears to be conserved in several T4-like phage groups. However, the group including phage T4 itself likely employs a different strategy whereby TAC-N and TAC-NC are encoded by separate genes (26 and 51 in phage T4).

Perspective on taxonomic classification of uncultivated viruses

Historically, virus taxonomy has been limited to describing viruses that were readily cultivated in the laboratory or emerging in natural biomes. Metagenomic analyses, single-particle sequencing, and database mining efforts have yielded new sequence data on an astounding number of previously unknown viruses. As metagenomes are relatively free of biases, these data provide an unprecedented insight into the vastness of the virosphere, but to properly value the extent of this diversity it is critical that the viruses are taxonomically classified. Inclusion of uncultivated viruses has already improved the process as well as the understanding of the taxa, viruses, and their evolutionary relationships. The continuous development and testing of computational tools will be required to maintain a dynamic virus taxonomy that can accommodate the new discoveries.

Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class Caudoviricetes, represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus-host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

A Novel Phage Infecting Alteromonas Represents a Distinct Group of Siphophages Infecting Diverse Aquatic Copiotrophs

Bacteriophages play critical roles in impacting microbial community succession both ecologically and evolutionarily. Although the majority of phage genetic diversity has been increasingly unveiled, phages infecting members of the ecologically important genus Alteromonas remain poorly understood. Here, we present a comprehensive analysis of a newly isolated alterophage, vB_AcoS-R7M (R7M), to characterize its life cycle traits, genomic features, and putative evolutionary origin. R7M harbors abundant genes identified as host-like auxiliary metabolic genes facilitating viral propagation. Genomic analysis suggested that R7M is distinct from currently known alterophages. Interestingly, R7M was found to share a set of similar characteristics with a number of siphophages infecting diverse aquatic opportunistic copiotrophs. We therefore proposed the creation of one new subfamily (Queuovirinae) to group with these evolutionarily related phages. Notably, tail genes were less likely to be shared among them, and baseplate-related genes varied the most. In-depth analyses indicated that R7M has replaced its distal tail with a Rhodobacter capsulatus gene transfer agent (RcGTA)-like baseplate and further acquired a putative receptor interaction site targeting Alteromonas. These findings suggest that horizontal exchanges of viral tail adsorption apparatuses are widespread and vital for phages to hunt new hosts and to adapt to new niches. IMPORTANCE The evolution and ecology of phages infecting members of Alteromonas, a marine opportunistic genus that is widely distributed and of great ecological significance, remain poorly understood. The present study integrates physiological and genomic evidence to characterize the properties and putative phage-host interactions of a newly isolated Alteromonas phage, vB_AcoS-R7M (R7M). A taxonomic study reveals close evolutionary relationships among R7M and a number of siphophages infecting various aquatic copiotrophs. Their similar head morphology and overall genetic framework suggest their putative common ancestry and the grouping of a new viral subfamily. However, their major difference lies in the viral tail adsorption apparatuses and the horizontal exchanges of which possibly account for variations in host specificity. These findings outline an evolutionary scenario for the emergence of diverse viral lineages of a shared genetic pool and give insights into the genetics and ecology of viral host jumps.

Individuals at risk for rheumatoid arthritis harbor differential intestinal bacteriophage communities with distinct metabolic potential

Rheumatoid arthritis (RA) is an autoimmune disease characterized in seropositive individuals by the presence of anti-cyclic citrullinated protein (CCP) antibodies. RA is linked to the intestinal microbiota, yet the association of microbes with CCP serology and their contribution to RA is unclear. We describe intestinal phage communities of individuals at risk for developing RA, with or without anti-CCP antibodies, whose first-degree relatives have been diagnosed with RA. We show that at-risk individuals harbor intestinal phage compositions that diverge based on CCP serology, are dominated by Streptococcaceae, Bacteroidaceae, and Lachnospiraceae phages, and may originate from disparate ecosystems. These phages encode unique repertoires of auxiliary metabolic genes, which associate with anti-CCP status, suggesting that these phages directly influence the metabolic and immunomodulatory capability of the microbiota. This work sets the stage for the use of phages as preclinical biomarkers and provides insight into a possible microbial-based causation of RA disease development.

Influence of Non-canonical DNA Bases on the Genomic Diversity of Tevenvirinae

The Tevenvirinae viruses are some of the most common viruses on Earth. Representatives of this subfamily have long been used in the molecular biology studies as model organisms – since the emergence of the discipline. Tevenvirinae are promising agents for phage therapy in animals and humans, since their representatives have only lytic life cycle and many of their host bacteria are pathogens. As confirmed experimentally, some Tevenvirinae have non-canonical DNA bases. Non-canonical bases can play an essential role in the diversification of closely related viruses. The article performs a comparative and evolutionary analysis of Tevenvirinae genomes and components of Tevenvirinae genomes. A comparative analysis of these genomes and the genes associated with the synthesis of non-canonical bases allows us to conclude that non-canonical bases have a major influence on the divergence of Tevenvirinae viruses within the same habitats. Supposedly, Tevenvirinae developed a strategy for changing HGT frequency in individual populations, which was based on the accumulation of proteins for the synthesis of non-canonical bases and proteins that used those bases as substrates. Owing to this strategy, ancestors of Tevenvirinae with the highest frequency of HGT acquired genes that allowed them to exist in a certain niche, and ancestors with the lowest HGT frequency preserved the most adaptive of those genes. Given the origin and characteristics of genes associated with the synthesis of non-canonical bases in Tevenvirinae, one can assume that other phages may have similar strategies. The article demonstrates the dependence of genomic diversity of closely related Tevenvirinae on non-canonical bases.

The genomic content and context of auxiliary metabolic genes in roseophages

Marine bacteriophages frequently possess auxiliary metabolic genes (AMGs) that accelerate host metabolism during phage infection. The significance of AMGs in phage infecting the ecologically important Roseobacter clade, found predominantly in marine environments, remains to be determined. Here, we analysed the distribution and genomic context of 180 AMGs, annotated into 20 types, across 50 roseophage genomes. Roseophages share seven high-frequency AMGs (trx, grx, RNR, thyX, DCD, phoH, and mazG)， most of them involved in the nucleotide biosynthesis pathway that represent conserved intra and inter operational taxonomic units (OTUs), and share ≥97% full-length DNA sequence similarity. Sporadic AMGs (dUTPase, lexA, degS, Que, NAPRT, AHL, pcnB, ctrA, RTX, RNR-nrdA, RNR-nrdE, wclP, and flgJ), present in only one or two OTUs, show high functional diversity. The roseophage AMG repertoire weakly correlates with environmental factors, while host range partially explains the sporadic AMG distribution. Locally co-linear blocks distribution index (LDI) analysis indicated that high-frequency roseopodovirus AMGs are restricted to particular genomic islands, possibly originating from limited historical acquisition events. Low-frequency roseopodovirus AMGs and all roseosiphovirus AMGs have high LDI values, implying multiple historical acquisition events. In summary, roseophages have acquired a range of AMGs through horizontal gene transfer, and the forces shaping the evolution of roseophages are described.

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.

Viral Ecogenomics of Arctic Cryopeg Brine and Sea Ice

Arctic regions, which are changing rapidly as they warm 2 to 3 times faster than the global average, still retain microbial habitats that serve as natural laboratories for understanding mechanisms of microbial adaptation to extreme conditions. Seawater-derived brines within both sea ice (sea-ice brine) and ancient layers of permafrost (cryopeg brine) support diverse microbes adapted to subzero temperatures and high salinities, yet little is known about viruses in these extreme environments, which, if analogous to other systems, could play important evolutionary and ecosystem roles. Here, we characterized viral communities and their functions in samples of cryopeg brine, sea-ice brine, and melted sea ice. Viral abundance was high in cryopeg brine (1.2 × 108 ml-1) and much lower in sea-ice brine (1.3 × 105 to 2.1 × 105 ml-1), which roughly paralleled the differences in cell concentrations in these samples. Five low-input, quantitative viral metagenomes were sequenced to yield 476 viral populations (i.e., species level; ≥10 kb), only 12% of which could be assigned taxonomy by traditional database approaches, indicating a high degree of novelty. Additional analyses revealed that these viruses: (i) formed communities that differed between sample type and vertically with sea-ice depth; (ii) infected hosts that dominated these extreme ecosystems, including Marinobacter, Glaciecola, and Colwellia; and (iii) encoded fatty acid desaturase (FAD) genes that likely helped their hosts overcome cold and salt stress during infection, as well as mediated horizontal gene transfer of FAD genes between microbes. Together, these findings contribute to understanding viral abundances and communities and how viruses impact their microbial hosts in subzero brines and sea ice.IMPORTANCE This study explores viral community structure and function in remote and extreme Arctic environments, including subzero brines within marine layers of permafrost and sea ice, using a modern viral ecogenomics toolkit for the first time. In addition to providing foundational data sets for these climate-threatened habitats, we found evidence that the viruses had habitat specificity, infected dominant microbial hosts, encoded host-derived metabolic genes, and mediated horizontal gene transfer among hosts. These results advance our understanding of the virosphere and how viruses influence extreme ecosystems. More broadly, the evidence that virally mediated gene transfers may be limited by host range in these extreme habitats contributes to a mechanistic understanding of genetic exchange among microbes under stressful conditions in other systems.

Isolation and complete genome sequence of a novel cyanophage, S-B05, infecting an estuarine Synechococcus strain: insights into environmental adaptation

A new cyanophage, S-B05, infecting a phycoerythrin-enriched (PE-type) Synechococcus strain was isolated by the liquid infection method, and its morphology and genetic features were examined. Phylogenetic analysis and morphological observation confirmed that S-B05 belongs to the family Myoviridae of the order Caudovirales. Its genome was fully sequenced, and found to be 208,857 bp in length with a G + C content of 39.9%. It contained 280 potential open reading frames and 123 conserved domains. Ninety-eight functional genes responsible for cyanophage structuring and packaging, DNA replication and regulation, and photosynthesis were identified, as well as genes encoding 172 hypothetical proteins. The genome of S-B05 is most similar to that of Prochlorococcus phage P-TIM68. Homologues of open reading frames of S-B05 can be found in various marine environments, as revealed by comparison of the S-B05 genome sequence to sequences in marine viral metagenomic databases. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as the phylogenetic relationships based on AMGs and the complete genome sequence, reflect the phage-host interaction mechanism or the specific adaptation strategy of the host to environmental conditions. The genome sequence information reported here will provide an important basis for further study of the adaptive evolution and ecological role of cyanophages and their hosts in the marine environment.

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.

Genomic and Seasonal Variations among Aquatic Phages Infecting the Baltic Sea Gammaproteobacterium Rheinheimera sp. Strain BAL341

Phages are important in aquatic ecosystems as they influence their microbial hosts through lysis, gene transfer, transcriptional regulation, and expression of phage metabolic genes. Still, there is limited knowledge of how phages interact with their hosts, especially at fine scales. Here, a Rheinheimera phage-host system constituting highly similar phages infecting one host strain is presented. This relatively limited diversity has previously been seen only when smaller numbers of phages have been isolated and points toward ecological constraints affecting the Rheinheimera phage diversity. The variation of metabolic genes among the species points toward various fitness advantages, opening up possibilities for future hypothesis testing. Phage-host dynamics monitored over several years point toward recurring “kill-the-winner” oscillations and an ecological niche fulfilled by this system in the Baltic Sea. Identifying and quantifying ecological dynamics of such phage-host model systems in situ allow us to understand and study the influence of phages on aquatic ecosystems.

Knowledge in aquatic virology has been greatly improved by culture-independent methods, yet there is still a critical need for isolating novel phages to identify the large proportion of “unknowns” that dominate metagenomes and for detailed analyses of phage-host interactions. Here, 54 phages infecting Rheinheimera sp. strain BAL341 (Gammaproteobacteria) were isolated from Baltic Sea seawater and characterized through genome content analysis and comparative genomics. The phages showed a myovirus-like morphology and belonged to a novel genus, for which we propose the name Barbavirus. All phages had similar genome sizes and numbers of genes (80 to 84 kb; 134 to 145 genes), and based on average nucleotide identity and genome BLAST distance phylogeny, the phages were divided into five species. The phages possessed several genes involved in metabolic processes and host signaling, such as genes encoding ribonucleotide reductase and thymidylate synthase, phoH, and mazG. One species had additional metabolic genes involved in pyridine nucleotide salvage, possibly providing a fitness advantage by further increasing the phages’ replication efficiency. Recruitment of viral metagenomic reads (25 Baltic Sea viral metagenomes from 2012 to 2015) to the phage genomes showed pronounced seasonal variations, with increased relative abundances of barba phages in August and September synchronized with peaks in host abundances, as shown by 16S rRNA gene amplicon sequencing. Overall, this study provides detailed information regarding genetic diversity, phage-host interactions, and temporal dynamics of an ecologically important aquatic phage-host system.

IMPORTANCE Phages are important in aquatic ecosystems as they influence their microbial hosts through lysis, gene transfer, transcriptional regulation, and expression of phage metabolic genes. Still, there is limited knowledge of how phages interact with their hosts, especially at fine scales. Here, a Rheinheimera phage-host system constituting highly similar phages infecting one host strain is presented. This relatively limited diversity has previously been seen only when smaller numbers of phages have been isolated and points toward ecological constraints affecting the Rheinheimera phage diversity. The variation of metabolic genes among the species points toward various fitness advantages, opening up possibilities for future hypothesis testing. Phage-host dynamics monitored over several years point toward recurring “kill-the-winner” oscillations and an ecological niche fulfilled by this system in the Baltic Sea. Identifying and quantifying ecological dynamics of such phage-host model systems in situ allow us to understand and study the influence of phages on aquatic ecosystems.

Rates of Molecular Evolution in a Marine Synechococcus Phage Lineage

Cyanophages are characterized by vast genomic diversity and the formation of stable ecotypes over time. The evolution of phage diversity includes vertical processes, such as mutation, and horizontal processes, such as recombination and gene transfer. Here, we study the contribution of vertical and horizontal processes to short-term evolution of marine cyanophages. Analyzing time series data of Synechococcus-infecting Myoviridae ecotypes spanning up to 17 years, we found a high contribution of recombination relative to mutation (r/m) in all ecotypes. Additionally, we found a molecular clock of substitution and recombination in one ecotype, RIM8. The estimated RIM8 evolutionary rates are 2.2 genome-wide substitutions per year (1.275 × 10⁻⁵ substitutions/site/year) and 29 genome-wide nucleotide alterations due to recombination per year. We found 26 variable protein families, of which only two families have a predicted functional annotation, suggesting that they are auxiliary metabolic genes with bacterial homologs. A comparison of our rate estimates to other phage evolutionary rate estimates in the literature reveals a negative correlation of phage substitution rates with their genome size. A comparison to evolutionary rates in bacterial organisms further shows that phages have high rates of mutation and recombination compared to their bacterial hosts. We conclude that the increased recombination rate in phages likely contributes to their vast genomic diversity.

Expanding the Diversity of Myoviridae Phages Infecting Lactobacillus plantarum-A Novel Lineage of Lactobacillus Phages Comprising Five New Members

Lactobacillus plantarum is a bacterium with probiotic properties and promising applications in the food industry and agriculture. So far, bacteriophages of this bacterium have been moderately addressed. We examined the diversity of five new L. plantarum phages via whole genome shotgun sequencing and in silico protein predictions. Moreover, we looked into their phylogeny and their potential genomic similarities to other complete phage genome records through extensive nucleotide and protein comparisons. These analyses revealed a high degree of similarity among the five phages, which extended to the vast majority of predicted virion-associated proteins. Based on these, we selected one of the phages as a representative and performed transmission electron microscopy and structural protein sequencing tests. Overall, the results suggested that the five phages belong to the family Myoviridae, they have a long genome of 137,973-141,344 bp, a G/C content of 36.3-36.6% that is quite distinct from their host's, and surprisingly, 7 to 15 tRNAs. Only an average 41/174 of their predicted genes were assigned a function. The comparative analyses unraveled considerable genetic diversity for the five L. plantarum phages in this study. Hence, the new genus "Semelevirus" was proposed, comprising exclusively of the five phages. This novel lineage of Lactobacillus phages provides further insight into the genetic heterogeneity of phages infecting Lactobacillus sp. The five new Lactobacillus phages have potential value for the development of more robust starters through, for example, the selection of mutants insensitive to phage infections. The five phages could also form part of phage cocktails, which producers would apply in different stages of L. plantarum fermentations in order to create a range of organoleptic outputs.

Adaptation to sub-optimal hosts is a driver of viral diversification in the ocean

Cyanophages of the Myoviridae family include generalist viruses capable of infecting a wide range of hosts including those from different cyanobacterial genera. While the influence of phages on host evolution has been studied previously, it is not known how the infection of distinct hosts influences the evolution of cyanophage populations. Here, using an experimental evolution approach, we investigated the adaptation of multiple cyanophage populations to distinct cyanobacterial hosts. We show that when infecting an "optimal" host, whose infection is the most efficient, phage populations accumulated only a few mutations. However, when infecting "sub-optimal" hosts, different mutations spread in the phage populations, leading to rapid diversification into distinct subpopulations. Based on our results, we propose a model demonstrating how shifts in microbial abundance, which lead to infection of "sub-optimal" hosts, act as a driver for rapid diversification of viral populations.

Proteorhodopsins dominate the expression of phototrophic mechanisms in seasonal and dynamic marine picoplankton communities

The most abundant and ubiquitous microbes in the surface ocean use light as an energy source, capturing it via complex chlorophyll-based photosystems or simple retinal-based rhodopsins. Studies in various ocean regimes compared the abundance of these mechanisms, but few investigated their expression. Here we present the first full seasonal study of abundance and expression of light-harvesting mechanisms (proteorhodopsin, PR; aerobic anoxygenic photosynthesis, AAnP; and oxygenic photosynthesis, PSI) from deep-sequenced metagenomes and metatranscriptomes of marine picoplankton (<1 µm) at three coastal stations of the San Pedro Channel in the Pacific Ocean. We show that, regardless of season or sampling location, the most common phototrophic mechanism in metagenomes of this dynamic region was PR (present in 65–104% of the genomes as estimated by single-copy recA), followed by PSI (5–104%) and AAnP (5–32%). Furthermore, the normalized expression (RNA to DNA ratio) of PR genes was higher than that of oxygenic photosynthesis (average ± standard deviation 26.2 ± 8.4 vs. 11 ± 9.7), and the expression of the AAnP marker gene was significantly lower than both mechanisms (0.013 ± 0.02). We demonstrate that PR expression was dominated by the SAR11-cluster year-round, followed by other Alphaproteobacteria, unknown-environmental clusters and Gammaproteobacteria. This highly dynamic system further allowed us to identify a trend for PR spectral tuning, in which blue-absorbing PR genes dominate in areas with low chlorophyll-a concentrations (<0.25 µgL⁻¹). This suggests that PR phototrophy is not an accessory function but instead a central mechanism that can regulate photoheterotrophic population dynamics.

Discovery of several novel, widespread, and ecologically distinct marine Thaumarchaeota viruses that encode amoC nitrification genes

Much of the diversity of prokaryotic viruses has yet to be described. In particular, there are no viral isolates that infect abundant, globally significant marine archaea including the phylum Thaumarchaeota. This phylum oxidizes ammonia, fixes inorganic carbon, and thus contributes to globally significant nitrogen and carbon cycles in the oceans. Metagenomics provides an alternative to culture-dependent means for identifying and characterizing viral diversity. Some viruses carry auxiliary metabolic genes (AMGs) that are acquired via horizontal gene transfer from their host(s), allowing inference of what host a virus infects. Here we present the discovery of 15 new genomically and ecologically distinct Thaumarchaeota virus populations, identified as contigs that encode viral capsid and thaumarchaeal ammonia monooxygenase genes (amoC). These viruses exhibit depth and latitude partitioning and are distributed globally in various marine habitats including pelagic waters, estuarine habitats, and hydrothermal plume water and sediments. We found evidence of viral amoC expression and that viral amoC AMGs sometimes comprise up to half of total amoC DNA copies in cellular fraction metagenomes, highlighting the potential impact of these viruses on N cycling in the oceans. Phylogenetics suggest they are potentially tailed viruses and share a common ancestor with related marine Euryarchaeota viruses. This work significantly expands our view of viruses of globally important marine Thaumarchaeota.

High Diversity of Myocyanophage in Various Aquatic Environments Revealed by High-Throughput Sequencing of Major Capsid Protein Gene With a New Set of Primers

Myocyanophages, a group of viruses infecting cyanobacteria, are abundant and play important roles in elemental cycling. Here we investigated the particle-associated viral communities retained on 0.2 μm filters and in sediment samples (representing ancient cyanophage communities) from four ocean and three lake locations, using high-throughput sequencing and a newly designed primer pair targeting a gene fragment (∼145-bp in length) encoding the cyanophage gp23 major capsid protein (MCP). Diverse viral communities were detected in all samples. The fragments of 142-, 145-, and 148-bp in length were most abundant in the amplicons, and most sequences (>92%) belonged to cyanophages. Additionally, different sequencing depths resulted in different diversity estimates of the viral community. Operational taxonomic units obtained from deep sequencing of the MCP gene covered the majority of those obtained from shallow sequencing, suggesting that deep sequencing exhibited a more complete picture of cyanophage community than shallow sequencing. Our results also revealed a wide geographic distribution of marine myocyanophages, i.e., higher dissimilarities of the myocyanophage communities corresponded with the larger distances between the sampling sites. Collectively, this study suggests that the newly designed primer pair can be effectively used to study the community and diversity of myocyanophage from different environments, and the high-throughput sequencing represents a good method to understand viral diversity.

The future of genomics in polar and alpine cyanobacteria

In recent years, genomic analyses have arisen as an exciting way of investigating the functional capacity and environmental adaptations of numerous micro-organisms of global relevance, including cyanobacteria. In the extreme cold of Arctic, Antarctic and alpine environments, cyanobacteria are of fundamental ecological importance as primary producers and ecosystem engineers. While their role in biogeochemical cycles is well appreciated, little is known about the genomic makeup of polar and alpine cyanobacteria. In this article, we present ways that genomic techniques might be used to further our understanding of cyanobacteria in cold environments in terms of their evolution and ecology. Existing examples from other environments (e.g. marine/hot springs) are used to discuss how methods developed there might be used to investigate specific questions in the cryosphere. Phylogenomics, comparative genomics and population genomics are identified as methods for understanding the evolution and biogeography of polar and alpine cyanobacteria. Transcriptomics will allow us to investigate gene expression under extreme environmental conditions, and metagenomics can be used to complement tradition amplicon-based methods of community profiling. Finally, new techniques such as single cell genomics and metagenome assembled genomes will also help to expand our understanding of polar and alpine cyanobacteria that cannot readily be cultured.

Co-existence of freshwater and marine T4-like myoviruses in a typical subtropical estuary

Viruses are the most abundant biological entities on Earth and play an important role in microbial community dynamics and biogeochemical cycling, yet their ecological characteristics in estuarine ecosystems are unclear. Here, virioplankton communities in a typical subtropical estuary, the Jiulong River estuary (JRE) in China, were investigated. The abundance of virioplankton ranged from 1.01 ± 0.05 × 107 to 1.62 ± 0.09 × 107 particles mL-1 in JRE, and the population size of viruses was correlated with temperature and nutrient levels. Three tailed viral morphotypes (myovirus, siphovirus and podovirus) were observed. Phylogenetic analysis showed that most of the g23 sequences in the JRE fell into three previously established groups (Marine, Paddy and Lake Groups) and two potential Estuary Groups. This demonstrates the co-existence of typical freshwater and marine T4-like myoviruses in the estuarine ecosystem, suggesting the movement of viruses and their hosts among biomes. Additionally, the spatial variation of g23 sequences suggests a geographic distribution pattern of T4-like myoviruses in the JRE, which might be shaped by the environmental gradient and/or their host distribution. These results provide valuable insights into the abundance, diversity and distribution patterns of virioplankton, as well as the factors influencing them, in subtropical estuarine ecosystems.

The Gut Microbiota Facilitates Drifts in the Genetic Diversity and Infectivity of Bacterial Viruses

The intestinal microbiota and human health are intimately linked, but interactions between bacteria and bacteriophages in the context of the mammalian intestine remain largely unexplored. We used comparative population genomics to study a tripartite network consisting of a virulent bacteriophage, its bacterial host, and a phage-insensitive bacterial strain both in vitro and within the murine gut. The bacteriophage adapted to infect the insensitive strain when the three partners co-existed in the gut of conventional mice, but not in dixenic mice or in planktonic cultures. The molecular changes associated with modifications in the bacteriophage host spectrum included single amino acid substitutions and an unusual homologous intragenomic recombination event within the genome of the bacteriophage. An intermediate bacterial host isolated from the murine microbiota mediated bacteriophage adaptation. Our data indicate that by offering access to new hosts, the microbiota shifts the genetic diversity of bacteriophages, thereby promoting long-term persistence of bacteriophage populations.

Cyanophage-encoded lipid desaturases: oceanic distribution, diversity and function

Cyanobacteria are among the most abundant photosynthetic organisms in the oceans; viruses infecting cyanobacteria (cyanophages) can alter cyanobacterial populations, and therefore affect the local food web and global biochemical cycles. These phages carry auxiliary metabolic genes (AMGs), which rewire various metabolic pathways in the infected host cell, resulting in increased phage fitness. Coping with stress resulting from photodamage appears to be a central necessity of cyanophages, yet the overall mechanism is poorly understood. Here we report a novel, widespread cyanophage AMG, encoding a fatty acid desaturase (FAD), found in two genotypes with distinct geographical distribution. FADs are capable of modulating the fluidity of the host’s membrane, a fundamental stress response in living cells. We show that both viral FAD (vFAD) families are Δ9 lipid desaturases, catalyzing the desaturation at carbon 9 in C16 fatty acid chains. In addition, we present a comprehensive fatty acid profiling for marine cyanobacteria, which suggests a unique desaturation pathway of medium- to long-chain fatty acids no longer than C16, in accordance with the vFAD activity. Our findings suggest that cyanophages are capable of fiddling with the infected host’s membranes, possibly leading to increased photoprotection and potentially enhancing viral-encoded photosynthetic proteins, resulting in a new viral metabolic network.