The genomic content and context of auxiliary metabolic genes in marine cyanomyoviruses

Genomic insights into redox-driven microbial processes for carbon decomposition in thawing Arctic soils and permafrost

Climate change is rapidly transforming Arctic landscapes where increasing soil temperatures speed up permafrost thaw. This exposes large carbon stocks to microbial decomposition, possibly worsening climate change by releasing more greenhouse gases. Understanding how microbes break down soil carbon, especially under the anaerobic conditions of thawing permafrost, is important to determine future changes. Here, we studied the microbial community dynamics and soil carbon decomposition potential in permafrost and active layer soils under anaerobic laboratory conditions that simulated an Arctic summer thaw. The microbial and viral compositions in the samples were analyzed based on metagenomes, metagenome-assembled genomes, and metagenomic viral contigs (mVCs). Following the thawing of permafrost, there was a notable shift in microbial community structure, with fermentative Firmicutes and Bacteroidota taking over from Actinobacteria and Proteobacteria over the 60-day incubation period. The increase in iron and sulfate-reducing microbes had a significant role in limiting methane production from thawed permafrost, underscoring the competition within microbial communities. We explored the growth strategies of microbial communities and found that slow growth was the major strategy in both the active layer and permafrost. Our findings challenge the assumption that fast-growing microbes mainly respond to environmental changes like permafrost thaw. Instead, they indicate a common strategy of slow growth among microbial communities, likely due to the thermodynamic constraints of soil substrates and electron acceptors, and the need for microbes to adjust to post-thaw conditions. The mVCs harbored a wide range of auxiliary metabolic genes that may support cell protection from ice formation in virus-infected cells.

IMPORTANCE

As the Arctic warms, thawing permafrost unlocks carbon, potentially accelerating climate change by releasing greenhouse gases. Our research delves into the underlying biogeochemical processes likely mediated by the soil microbial community in response to the wet and anaerobic conditions, akin to an Arctic summer thaw. We observed a significant shift in the microbial community post-thaw, with fermentative bacteria like Firmicutes and Bacteroidota taking over and switching to different fermentation pathways. The dominance of iron and sulfate-reducing bacteria likely constrained methane production in the thawing permafrost. Slow-growing microbes outweighed fast-growing ones, even after thaw, upending the expectation that rapid microbial responses to dominate after permafrost thaws. This research highlights the nuanced and complex interactions within Arctic soil microbial communities and underscores the challenges in predicting microbial response to environmental change.

Novel phages of Pseudomonas syringae unveil numerous potential auxiliary metabolic genes

Relatively few phages that infect plant pathogens have been isolated and investigated. The Pseudomonas syringae species complex is present in various environments, including plants. It can cause major crop diseases, such as bacterial canker on apricot trees. This study presents a collection of 25 unique phage genomes that infect P. syringae. These phages were isolated from apricot orchards with bacterial canker symptoms after enrichment with 21 strains of P. syringae. This collection comprises mostly virulent phages, with only three being temperate. They belong to 14 genera, 11 of which are newly discovered, and 18 new species, revealing great genetic diversity within this collection. Novel DNA packaging systems have been identified bioinformatically in one of the new phage species, but experimental confirmation is required to define the precise mechanism. Additionally, many phage genomes contain numerous potential auxiliary metabolic genes with diversified putative functions. At least three phages encode genes involved in bacterial tellurite resistance, a toxic metalloid. This suggests that viruses could play a role in bacterial stress tolerance. This research emphasizes the significance of continuing the search for new phages in the agricultural ecosystem to unravel novel ecological diversity and new gene functions. This work contributes to the foundation for future fundamental and applied research on phages infecting phytopathogenic bacteria.

Benchmarking bioinformatic virus identification tools using real-world metagenomic data across biomes

BACKGROUND

As most viruses remain uncultivated, metagenomics is currently the main method for virus discovery. Detecting viruses in metagenomic data is not trivial. In the past few years, many bioinformatic virus identification tools have been developed for this task, making it challenging to choose the right tools, parameters, and cutoffs. As all these tools measure different biological signals, and use different algorithms and training and reference databases, it is imperative to conduct an independent benchmarking to give users objective guidance.

RESULTS

We compare the performance of nine state-of-the-art virus identification tools in thirteen modes on eight paired viral and microbial datasets from three distinct biomes, including a new complex dataset from Antarctic coastal waters. The tools have highly variable true positive rates (0-97%) and false positive rates (0-30%). PPR-Meta best distinguishes viral from microbial contigs, followed by DeepVirFinder, VirSorter2, and VIBRANT. Different tools identify different subsets of the benchmarking data and all tools, except for Sourmash, find unique viral contigs. Performance of tools improved with adjusted parameter cutoffs, indicating that adjustment of parameter cutoffs before usage should be considered.

CONCLUSIONS

Together, our independent benchmarking facilitates selecting choices of bioinformatic virus identification tools and gives suggestions for parameter adjustments to viromics researchers.

The insights into the phage communities of fermented foods in the age of viral metagenomics

Phages play a critical role in the assembly and regulation of fermented food microbiome through lysis and lysogenic lifestyle, which in turn affects the yield and quality of fermented foods. Therefore, it is important to investigate and characterize the diversity and function of phages under complex microbial communities and nutrient substrate conditions to provide novel insights into the regulation of traditional spontaneous fermentation. Viral metagenomics has gradually garnered increasing attention in fermented food research to elucidate phage functions and characterize the interactions between phages and the microbial community. Advances in this technology have uncovered a wide range of phages associated with the production of traditional fermented foods and beverages. This paper reviews the common methods of viral metagenomics applied in fermented food research, and summarizes the ecological functions of phages in traditional fermented foods. In the future, combining viral metagenomics with culturable methods and metagenomics will broaden the scope of research on fermented food systems, revealing the complex role of phages and intricate phage-bacterium interactions.

Coordinated transcriptional response to environmental stress by a Synechococcus virus

Viruses are a major control on populations of microbes. Often, their virulence is examined in controlled laboratory conditions. Yet, in nature, environmental conditions lead to changes in host physiology and fitness that may impart both costs and benefits on viral success. Phosphorus (P) is a major abiotic control on the marine cyanobacterium Synechococcus. Some viruses infecting Synechococcus have acquired, from their host, a gene encoding a P substrate binding protein (PstS), thought to improve virus replication under phosphate starvation. Yet, pstS is uncommon among cyanobacterial viruses. Thus, we asked how infections with viruses lacking PstS are affected by P scarcity. We show that the production of infectious virus particles of such viruses is reduced in low P conditions. However, this reduction in progeny is not caused by impaired phage genome replication, thought to be a major sink for cellular phosphate. Instead, transcriptomic analysis showed that under low P conditions, a PstS-lacking cyanophage increased the expression of a specific gene set that included mazG, hli2, and gp43 encoding a pyrophosphatase, a high-light inducible protein and DNA polymerase, respectively. Moreover, several of the upregulated genes were controlled by the host's phoBR two-component system. We hypothesize that recycling and polymerization of nucleotides liberates free phosphate and thus allows viral morphogenesis, albeit at lower rates than when phosphate is replete or when phages encode pstS. Altogether, our data show how phage genomes, lacking obvious P-stress-related genes, have evolved to exploit their host's environmental sensing mechanisms to coordinate their own gene expression in response to resource limitation.

Viruses under the Antarctic Ice Shelf are active and potentially involved in global nutrient cycles

Viruses play an important role in the marine ecosystem. However, our comprehension of viruses inhabiting the dark ocean, and in particular, under the Antarctic Ice Shelves, remains limited. Here, we mine single-cell genomic, transcriptomic, and metagenomic data to uncover the viral diversity, biogeography, activity, and their role as metabolic facilitators of microbes beneath the Ross Ice Shelf. This is the largest Antarctic ice shelf with a major impact on global carbon cycle. The viral community found in the cavity under the ice shelf mainly comprises endemic viruses adapted to polar and mesopelagic environments. The low abundance of genes related to lysogenic lifestyle (<3%) does not support a predominance of the Piggyback-the-Winner hypothesis, consistent with a low-productivity habitat. Our results indicate a viral community actively infecting key ammonium and sulfur-oxidizing chemolithoautotrophs (e.g. Nitrosopumilus spp, Thioglobus spp.), supporting a "kill-the-winner" dynamic. Based on genome analysis, these viruses carry specific auxiliary metabolic genes potentially involved in nitrogen, sulfur, and phosphorus acquisition. Altogether, the viruses under Antarctic ice shelves are putatively involved in programming the metabolism of ecologically relevant microbes that maintain primary production in these chemosynthetically-driven ecosystems, which have a major role in global nutrient cycles.

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.

Depth-driven patterns in lytic viral diversity, auxiliary metabolic gene content, and productivity in offshore oligotrophic waters

Introduction

Marine viruses regulate microbial population dynamics and biogeochemical cycling in the oceans. The ability of viruses to manipulate hosts' metabolism through the expression of viral auxiliary metabolic genes (AMGs) was recently highlighted, having important implications in energy production and flow in various aquatic environments. Up to now, the presence and diversity of viral AMGs is studied using -omics data, and rarely using quantitative measures of viral activity alongside.

Methods

In the present study, four depth layers (5, 50, 75, and 1,000 m) with discrete hydrographic features were sampled in the Eastern Mediterranean Sea; we studied lytic viral community composition and AMG content through metagenomics, and lytic production rates through the viral reduction approach in the ultra-oligotrophic Levantine basin where knowledge regarding viral actions is rather limited.

Results and Discussion

Our results demonstrate depth-dependent patterns in viral diversity and AMG content, related to differences in temperature, nutrients availability, and host bacterial productivity and abundance. Although lytic viral production rates were similar along the water column, the virus-to-bacteria ratio was higher and the particular set of AMGs was more diverse in the bathypelagic (1,000 m) than the shallow epipelagic (5, 50, and 75 m) layers, revealing that the quantitative effect of viruses on their hosts may be the same along the water column through the intervention of different AMGs. In the resource- and energy-limited bathypelagic waters of the Eastern Mediterranean, the detected AMGs could divert hosts' metabolism toward energy production, through a boost in gluconeogenesis, fatty-acid and glycan biosynthesis and metabolism, and sulfur relay. Near the deep-chlorophyll maximum depth, an exceptionally high percentage of AMGs related to photosynthesis was noticed. Taken together our findings suggest that the roles of viruses in the deep sea might be even more important than previously thought as they seem to orchestrate energy acquisition and microbial community dynamics, and thus, biogeochemical turnover in the oceans.

Narrow host range phages infect essential bacteria for water purification reactions in groundwater-fed rapid sand filters

Biofiltration is used worldwide to provide safe potable water due to its low energy demand and excellent treatment performance. For instance, in Denmark, over 95% of drinking water is supplied through groundwater-fed rapid sand filters (RSF). Bacteriophages, viruses that infect bacteria, have been shown to shape the taxonomic and functional composition of microbial communities across a range of natural and engineering systems. However, phages in the biofiltration systems are rarely studied, despite the central role microbes play in water purification. To probe this, metagenomic data from surface water, groundwater and mixed source water biofiltration units (n = 26 from China, Europe and USA) for drinking water production were analysed to characterize prokaryotic viruses and to identify their potential microbial hosts. The source water type and geographical location are found to exert influence on the composition of the phageome in biofilters. Although the viral abundance (71,676 ± 17,841 RPKM) in biofilters is only 14.4% and 17.0% lower than those of the nutrient-rich wastewater treatment plants and fresh surface waters, the richness (1,441 ± 1,046) and diversity (Inverse Simpson: 91 ± 61) in biofiltration units are significantly less by a factor of 2-5 and 3-4, respectively. In depth analysis of data from 24 groundwater-fed RSFs in Denmark revealed a core phageome shared by most RSFs, which was consistently linked to dominant microbial hosts involved in key biological reactions for water purification. Finally, the high number of specific links detected between phages and bacterial species and the large proportion of lytic phages (77%) led to the conjecture that phages regulate bacterial populations through predation, preventing the proliferation of dominant species and contributing to the established functional redundancy among the dominant microbial groups. In conclusion, bacteriophages are likely to play a significant role in water treatment within biofilters, particularly through interactions with key bacterial species.

Unveiling the hidden role of aquatic viruses in hydrocarbon pollution bioremediation

Hydrocarbon pollution poses substantial environmental risks to water and soil. Bioremediation, which utilizes microorganisms to manage pollutants, offers a cost-effective solution. However, the role of viruses, particularly bacteriophages (phages), in bioremediation remains unexplored. This study examines the diversity and activity of hydrocarbon-degradation genes encoded by environmental viruses, focusing on phages, within public databases. We identified 57 high-quality phage-encoded auxiliary metabolic genes (AMGs) related to hydrocarbon degradation, which we refer to as virus-encoded hydrocarbon degradation genes (vHYDEGs). These genes are encoded by taxonomically diverse aquatic phages and highlight the under-characterized global virosphere. Six protein families involved in the initial alkane hydroxylation steps were identified. Phylogenetic analyses revealed the diverse evolutionary trajectories of vHYDEGs across habitats, revealing previously unknown biodegraders linked evolutionarily with vHYDEGs. Our findings suggest phage AMGs may contribute to alkane and aromatic hydrocarbon degradation, participating in the initial, rate-limiting hydroxylation steps, thereby aiding hydrocarbon pollution bioremediation and promoting their propagation. To support future research, we developed vHyDeg, a database containing identified vHYDEGs with comprehensive annotations, facilitating the screening of hydrocarbon degradation AMGs and encouraging their bioremediation applications.

Ubiquitous, B12-dependent virioplankton utilizing ribonucleotide-triphosphate reductase demonstrate interseasonal dynamics and associate with a diverse range of bacterial hosts in the pelagic ocean

Through infection and lysis of their coexisting bacterial hosts, viruses impact the biogeochemical cycles sustaining globally significant pelagic oceanic ecosystems. Currently, little is known of the ecological interactions between lytic viruses and their bacterial hosts underlying these biogeochemical impacts at ecosystem scales. This study focused on populations of lytic viruses carrying the B12-dependent Class II monomeric ribonucleotide reductase (RNR) gene, ribonucleotide-triphosphate reductase (Class II RTPR), documenting seasonal changes in pelagic virioplankton and bacterioplankton using amplicon sequences of Class II RTPR and the 16S rRNA gene, respectively. Amplicon sequence libraries were analyzed using compositional data analysis tools that account for the compositional nature of these data. Both virio- and bacterioplankton communities responded to environmental changes typically seen across seasonal cycles as well as shorter term upwelling-downwelling events. Defining Class II RTPR-carrying viral populations according to major phylogenetic clades proved a more robust means of exploring virioplankton ecology than operational taxonomic units defined by percent sequence homology. Virioplankton Class II RTPR populations showed positive associations with a broad phylogenetic diversity of bacterioplankton including dominant taxa within pelagic oceanic ecosystems such as Prochlorococcus and SAR11. Temporal changes in Class II RTPR virioplankton, occurring as both free viruses and within infected cells, indicated possible viral-host pairs undergoing sustained infection and lysis cycles throughout the seasonal study. Phylogenetic relationships inferred from Class II RTPR sequences mirrored ecological patterns in virio- and bacterioplankton populations demonstrating possible genome to phenome associations for an essential viral replication gene.

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.

Functional Profiling and Evolutionary Analysis of a Marine Microalgal Virus Pangenome

Phycodnaviridae are large double-stranded DNA viruses, which facilitate studies of host-virus interactions and co-evolution due to their prominence in algal infection and their role in the life cycle of algal blooms. However, the genomic interpretation of these viruses is hampered by a lack of functional information, stemming from the surprising number of hypothetical genes of unknown function. It is also unclear how many of these genes are widely shared within the clade. Using one of the most extensively characterized genera, Coccolithovirus, as a case study, we combined pangenome analysis, multiple functional annotation tools, AlphaFold structural modeling, and literature analysis to compare the core and accessory pangenome and assess support for novel functional predictions. We determined that the Coccolithovirus pangenome shares 30% of its genes with all 14 strains, making up the core. Notably, 34% of its genes were found in at most three strains. Core genes were enriched in early expression based on a transcriptomic dataset of Coccolithovirus EhV-201 algal infection, were more likely to be similar to host proteins than the non-core set, and were more likely to be involved in vital functions such as replication, recombination, and repair. In addition, we generated and collated annotations for the EhV representative EhV-86 from 12 different annotation sources, building up information for 142 previously hypothetical and putative membrane proteins. AlphaFold was further able to predict structures for 204 EhV-86 proteins with a modelling accuracy of good-high. These functional clues, combined with generated AlphaFold structures, provide a foundational framework for the future characterization of this model genus (and other giant viruses) and a further look into the evolution of the Coccolithovirus proteome.

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.

Genetic diversity of virus auxiliary metabolism genes associated with phosphorus metabolism in Napahai plateau wetland

Viruses play important roles in ecosystems by interfering with the central metabolic pathways of the host during infection via the expression of auxiliary metabolic genes (AMGs), altering the productivity of ecosystems and thus affecting geochemical cycling. In this study, the genetic diversity of phosphorus metabolism AMGs phoH, phoU and pstS was investigated by phylogenetic analysis, PCoA analysis, and alpha diversity analysis based on metagenomic data. It was found that the majority of the sequences were unique to Napahai plateau wetland. It was shown that the genetic diversity of phoH, phoU and pstS genes was independent of both habitats and host origins. In addition, the metabolic pathway of AMGs associated with the phosphorus cycling was identified based on metagenomic data. When phosphorus is deficient, virus utilizes AMGs to affect the metabolic pathway, contributing to higher phosphorus levels in the host and facilitating virus survival, replication, and propagation in the host cell.

Diversity, evolution and life strategies of CbK-like phages

Caulobacter phage CbK has been extensively studied as a model system in virology and bacteriology. Lysogeny-related genes have been found in each CbK-like isolate, suggesting a life strategy of both lytic and lysogenic cycles. However, whether CbK-related phages can enter lysogeny is still undetermined. This study identified new CbK-like sequences and expanded the collection of CbK-related phages. A common ancestry with a temperate lifestyle was predicted for the group, however, which subsequently evolved into two clades of different genome sizes and host associations. Through the examination of phage recombinase genes, alignment of attachment sites on the phage and bacterial genomes (attP-attB pairing), and the experimental validation, different lifestyles were found among the different members. A majority of clade II members retain a lysogenic lifestyle, whereas all clade I members have evolved into an obligate lytic lifestyle via a loss of the gene encoding Cre-like recombinase and the coupled attP fragment. We postulated that the loss of lysogeny may be a by-product of the increase in phage genome size, and vice versa. Clade I is likely to overcome the costs through maintaining more auxiliary metabolic genes (AMGs), particularly for those involved in protein metabolism, to strengthen host takeover and further benefit virion production.

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 and ecological importance of viruses in the marine environment]

The ocean is the largest reservoir of viruses on the planet with estimates of up to several billions per liter. These viruses represent a major driving force not only for the evolution and for structuring the microbial world, but also for the functioning and the balance of marine ecosystems. With the advances in high throughput sequencing techniques, we are beginning to uncover the diversity and the complexity of this marine virosphere. This review synthesizes milestones in the field of marine viral ecology, including the diversity of these fascinating microorganisms, their impact on microbial mortality and cycling of nutrients and energy in the ocean.

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.

Viral community-wide auxiliary metabolic genes differ by lifestyles, habitats, and hosts

BACKGROUND

Viral-encoded auxiliary metabolic genes (AMGs) are important toolkits for modulating their hosts' metabolisms and the microbial-driven biogeochemical cycles. Although the functions of AMGs have been extensively reported in numerous environments, we still know little about the drivers that shape the viral community-wide AMG compositions in natural ecosystems. Exploring the drivers of viral community-wide AMG compositions is critical for a deeper understanding of the complex interplays among viruses, hosts, and the environments.

RESULTS

Here, we investigated the impact of viral lifestyles (i.e., lytic and lysogenic), habitats (i.e., water, particle, and sediment), and prokaryotic hosts on viral AMG profiles by utilizing metagenomic and metatranscriptomic techniques. We found that viral lifestyles were the most important drivers, followed by habitats and host identities. Specifically, irrespective of what habitats viruses came from, lytic viruses exhibited greater AMG diversity and tended to encode AMGs for chaperone biosynthesis, signaling proteins, and lipid metabolism, which could boost progeny reproduction, whereas temperate viruses were apt to encode AMGs for host survivability. Moreover, the lytic and temperate viral communities tended to mediate the microbial-driven biogeochemical cycles, especially nitrogen metabolism, in different manners via AMGs. When focusing on each lifestyle, we further found clear dissimilarity in AMG compositions between water and sediment, as well the divergent AMGs encoded by viruses infecting different host orders.

CONCLUSIONS

Overall, our study provides a first systematic characterization of the drivers of viral community-wide AMG compositions and further expands our knowledge of the distinct interactions of lytic and temperate viruses with their prokaryotic hosts from an AMG perspective, which is critical for understanding virus-host-environment interactions in natural conditions. Video Abstract.

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.

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.

Structural characterization of a soil viral auxiliary metabolic gene product - a functional chitosanase

Metagenomics is unearthing the previously hidden world of soil viruses. Many soil viral sequences in metagenomes contain putative auxiliary metabolic genes (AMGs) that are not associated with viral replication. Here, we establish that AMGs on soil viruses actually produce functional, active proteins. We focus on AMGs that potentially encode chitosanase enzymes that metabolize chitin - a common carbon polymer. We express and functionally screen several chitosanase genes identified from environmental metagenomes. One expressed protein showing endo-chitosanase activity (V-Csn) is crystalized and structurally characterized at ultra-high resolution, thus representing the structure of a soil viral AMG product. This structure provides details about the active site, and together with structure models determined using AlphaFold, facilitates understanding of substrate specificity and enzyme mechanism. Our findings support the hypothesis that soil viruses contribute auxiliary functions to their hosts.

Climate-Endangered Arctic Epishelf Lake Harbors Viral Assemblages with Distinct Genetic Repertoires

Milne Fiord, located on the coastal margin of the Last Ice Area (LIA) in the High Arctic (82°N, Canada), harbors an epishelf lake, a rare type of ice-dependent ecosystem in which a layer of freshwater overlies marine water connected to the open ocean. This microbe-dominated ecosystem faces catastrophic change due to the deterioration of its ice environment related to warming temperatures. We produced the first assessment of viral abundance, diversity, and distribution in this vulnerable ecosystem and explored the niches available for viral taxa and the functional genes underlying their distribution. We found that the viral community in the freshwater layer was distinct from, and more diverse than, the community in the underlying seawater and contained a different set of putative auxiliary metabolic genes, including the sulfur starvation-linked gene tauD and the gene coding for patatin-like phospholipase. The halocline community resembled the freshwater more than the marine community, but harbored viral taxa unique to this layer. We observed distinct viral assemblages immediately below the halocline, at a depth that was associated with a peak of prasinophyte algae and the viral family Phycodnaviridae. We also assembled 15 complete circular genomes, including a putative Pelagibacter phage with a marine distribution. It appears that despite its isolated and precarious situation, the varied niches in this epishelf lake support a diverse viral community, highlighting the importance of characterizing underexplored microbiota in the Last Ice Area before these ecosystems undergo irreversible change.

IMPORTANCE Viruses are key to understanding polar aquatic ecosystems, which are dominated by microorganisms. However, studies of viral communities are challenging to interpret because the vast majority of viruses are known only from sequence fragments, and their taxonomy, hosts, and genetic repertoires are unknown. Our study establishes a basis for comparison that will advance understanding of viral ecology in diverse global environments, particularly in the High Arctic. Rising temperatures in this region mean that researchers have limited time remaining to understand the biodiversity and biogeochemical cycles of ice-dependent environments and the consequences of these rapid, irreversible changes. The case of the Milne Fiord epishelf lake has special urgency because of the rarity of this type of “floating lake” ecosystem and its location in the Last Ice Area, a region of thick sea ice with global importance for conservation efforts.

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.

Viral diversity and biogeochemical potential revealed in different prawn-culture sediments by virus-enriched metagenome analysis

As the most numerous biological entities on Earth, viruses affect the microbial dynamics, metabolism and biogeochemical cycles in the aquatic ecosystems. Viral diversity and functions in ocean have been relatively well studied, but our understanding of viruses in mariculture systems is limited. To fill this knowledge gap, we studied viral diversity and potential biogeochemical impacts of sediments from four different prawn-mariculture ecosystems (mono-culture of prawn and poly-culture of prawn with jellyfish, sea cucumber, and clam) using a metagenomic approach with prior virus-like particles (VLPs) separation. We found that the order Caudovirales was the predominant viral category and accounted for the most volume (78.39% of classified viruses). Sediment viruses were verified to have a high diversity by using the construct phylogenetic tree of terL gene, with three potential novel clades being identified. Meanwhile, compared with viruses inhabiting other ecosystems based on gene-sharing network, our results revealed that mariculture sediments harbored considerable unexplored viral diversity and that maricultural species were potentially important drivers of the viral community structure. Notably, viral auxiliary metabolic genes were identified and suggested that viruses influence carbon and sulfur cycling, as well as cofactors/vitamins and amino acid metabolism, which indirectly participate in biogeochemical cycling. Overall, our findings revealed the genomic diversity and ecological function of viral communities in prawn mariculture sediments, and suggested the role of viruses in microbial ecology and biogeochemistry.

Droplet Microfluidics for High-Resolution Virology

Microfluidics has enabled a new era of cellular and molecular assays due to the small length scales, parallelization, and the modularity of various analysis and actuation functions. Droplet microfluidics, in particular, has been instrumental in providing new tools for biology with its ability to quickly and reproducibly generate drops that act as individual reactors. A notable beneficiary of this technology has been single-cell RNA sequencing, which has revealed new heterogeneities and interactions for the fundamental unit of life. However, viruses far surpass the diversity of cellular life, affect the dynamics of all ecosystems, and are a chronic source of global health crises. Despite their impact on the world, high-throughput and high-resolution viral profiling has been difficult, with conventional methods being limited to population-level averaging, large sample volumes, and few cultivable hosts. Consequently, most viruses have not been identified and studied. Droplet microfluidics holds the potential to address many of these limitations and offers new levels of sensitivity and throughput for virology. This Feature highlights recent efforts that have applied droplet microfluidics to the detection and study of viruses, including for diagnostics, virus-host interactions, and cell-independent virus assays. In combination with traditional virology methods, droplet microfluidics should prove a potent tool toward achieving a better understanding of the most abundant biological species on Earth.

Global Composition of the Bacteriophage Community in Honey Bees

The microbial communities in animal digestive systems are critical for host development and health. They stimulate the immune system during development, synthesize important chemical compounds like hormones, aid in digestion, competitively exclude pathogens, etc. Compared to the bacterial and fungal components of the microbiome, we know little about the temporal and spatial dynamics of bacteriophage communities in animal digestive systems. Recently, the bacteriophages of the honey bee gut were characterized in two European bee populations. Most of the bacteriophages described in these two reports were novel, harbored many metabolic genes in their genomes, and had a community structure that suggests coevolution with their bacterial hosts. To describe the conservation of bacteriophages in bees and begin to understand their role in the bee microbiome, we sequenced the virome of Apis mellifera from Austin, TX, and compared bacteriophage compositions among three locations around the world. We found that most bacteriophages from Austin are novel, sharing no sequence similarity with anything in public repositories. However, many bacteriophages are shared among the three bee viromes, indicating specialization of bacteriophages in the bee gut. Our study, along with the two previous bee virome studies, shows that the bee gut bacteriophage community is simple compared to that of many animals, consisting of several hundred types of bacteriophages that primarily infect four of the dominant bacterial phylotypes in the bee gut. IMPORTANCE Viruses that infect bacteria (bacteriophages) are abundant in the microbial communities that live on and in plants and animals. However, our knowledge of the structure, dynamics, and function of these viral communities lags far behind our knowledge of their bacterial hosts. We sequenced the first bacteriophage community of honey bees from the United States and compared the U.S. honey bee bacteriophage community to those of samples from Europe. Our work is an important characterization of an economically critical insect species and shows how bacteriophage communities can contain highly conserved individuals and be highly variable in composition across a wide geographic range.

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.

Viruses affect picocyanobacterial abundance and biogeography in the North Pacific Ocean

The photosynthetic picocyanobacteria Prochlorococcus and Synechococcus are models for dissecting how ecological niches are defined by environmental conditions, but how interactions with bacteriophages affect picocyanobacterial biogeography in open ocean biomes has rarely been assessed. We applied single-virus and single-cell infection approaches to quantify cyanophage abundance and infected picocyanobacteria in 87 surface water samples from five transects that traversed approximately 2,200 km in the North Pacific Ocean on three cruises, with a duration of 2-4 weeks, between 2015 and 2017. We detected a 550-km-wide hotspot of cyanophages and virus-infected picocyanobacteria in the transition zone between the North Pacific Subtropical and Subpolar gyres that was present in each transect. Notably, the hotspot occurred at a consistent temperature and displayed distinct cyanophage-lineage composition on all transects. On two of these transects, the levels of infection in the hotspot were estimated to be sufficient to substantially limit the geographical range of Prochlorococcus. Coincident with the detection of high levels of virally infected picocyanobacteria, we measured an increase of 10-100-fold in the Synechococcus populations in samples that are usually dominated by Prochlorococcus. We developed a multiple regression model of cyanophages, temperature and chlorophyll concentrations that inferred that the hotspot extended across the North Pacific Ocean, creating a biological boundary between gyres, with the potential to release organic matter comparable to that of the sevenfold-larger North Pacific Subtropical Gyre. Our results highlight the probable impact of viruses on large-scale phytoplankton biogeography and biogeochemistry in distinct regions of the oceans.

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.

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.

Comparative genomic analysis of dwarf Vibrio myoviruses defines a conserved gene cluster for successful phage infection

Tailed bacteriophages have been at the center of attention, not only for their ability to infect and kill pathogenic bacteria but also due to their peculiar and intriguing complex contractile tail structure. Tailed bacteriophages with contractile tails are known to have a Myoviridae morphotype and are members of the order Caudovirales. Large bacteriophages with a genome larger than 150 kbp have been studied for their ability to use multiple infection and lysis strategies to replicate more efficiently. On the other hand, smaller bacteriophages with fewer genes are represented in the GenBank database in greater numbers, and have several genes with unknown function. Isolation and molecular characterization of a newly reported bacteriophage named Athena1 revealed that it is a strongly lytic bacteriophage with a genome size of 39,826 bp. This prompted us to perform a comparative genomic analysis of Vibrio myoviruses with a genome size of no more than 50 kbp. The results revealed a pattern of genomic organization that includes sets of genes responsible for virion morphogenesis, replication/recombination of DNA, and lysis/lysogeny switching. By studying phylogenetic gene markers, we were able to draw conclusions about evolutionary events that shaped the genomic mosaicism of these phages, pinpointing the importance of a conserved organization of the genomic region encoding the baseplate protein for successful infection of Gram-negative bacteria. In addition, we propose the creation of new genera for dwarf Vibrio myoviruses. Comparative genomics of phages infecting aquatic bacteria could provide information that is useful for combating fish pathogens in aquaculture, using novel strategies.

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.

Genomic Characterization of Two Novel RCA Phages Reveals New Insights into the Diversity and Evolution of Marine Viruses

Viruses are the most abundant living entities in marine ecosystems, playing critical roles in altering the structure and function of microbial communities and driving ocean biogeochemistry. Phages that infect Roseobacter clade-affiliated (RCA) cluster strains are an important component of marine viral communities. Here, we characterize the genome sequences of two new RCA phages, CRP-9 and CRP-13, which infect RCA strain FZCC0023. Genomic analysis reveals that CRP-9 and CRP-13 represent a novel evolutionary lineage of marine phages. They both have a DNA replication module most similar to those in Cobavirus group phages. In contrast, their morphogenesis and packaging modules are distinct from those in cobaviruses but homologous to those in HMO-2011-type phages. The genomic architecture of CRP-9 and CRP-13 suggests a genomic recombination event between distinct phage groups. Metagenomic data sets were examined for metagenome-assembled viral genomes (MAVGs) with similar recombinant genome architectures. Fifteen CRP-9-type MAVGs were identified from marine viromes. Additionally, 158 MAVGs were identified containing HMO-2011-type morphogenesis and packaging modules with other types of DNA replication genes, providing more evidence that recombination between different phage groups is a major driver of phage evolution. Altogether, this study significantly expands the understanding of diversity and evolution of marine roseophages. Meanwhile, the analysis of these novel RCA phages and MAVGs highlights the critical role of recombination in shaping phage diversity. These phage sequences are valuable resources for inferring the evolutionary connection of distinct phage groups. IMPORTANCE Diversity and evolution of phages that infect the relatively slow-growing but dominant Roseobacter lineages are largely unknown. In this study, RCA phages CRP-9 and CRP-13 have been isolated on a Roseobacter RCA strain and shown to have a unique genomic architecture, which appears to be the result of a recombination event. CRP-9 and CRP-13 have a DNA replication module most similar to those in Cobavirus group phages and morphogenesis and packaging modules most similar to those in HMO-2011-type phages. HMO-2011-type morphogenesis and packaging modules are found in combination with distinct types of DNA replication genes, suggesting compatibility with various DNA replication modules. Altogether, this study contributes toward a better understanding of marine viral diversity and evolution.

Viruses and Their Interactions With Bacteria and Archaea of Hypersaline Great Salt Lake

Viruses play vital biogeochemical and ecological roles by (a) expressing auxiliary metabolic genes during infection, (b) enhancing the lateral transfer of host genes, and (c) inducing host mortality. Even in harsh and extreme environments, viruses are major players in carbon and nutrient recycling from organic matter. However, there is much that we do not yet understand about viruses and the processes mediated by them in the extreme environments such as hypersaline habitats. The Great Salt Lake (GSL) in Utah, United States is a hypersaline ecosystem where the biogeochemical role of viruses is poorly understood. This study elucidates the diversity of viruses and describes virus–host interactions in GSL sediments along a salinity gradient. The GSL sediment virosphere consisted of Haloviruses (32.07 ± 19.33%) and members of families Siphoviridae (39.12 ± 19.8%), Myoviridae (13.7 ± 6.6%), and Podoviridae (5.43 ± 0.64%). Our results demonstrate that salinity alongside the concentration of organic carbon and inorganic nutrients (nitrogen and phosphorus) governs the viral, bacteria, and archaeal diversity in this habitat. Computational host predictions for the GSL viruses revealed a wide host range with a dominance of viruses that infect Proteobacteria, Actinobacteria, and Firmicutes. Identification of auxiliary metabolic genes for photosynthesis (psbA), carbon fixation (rbcL, cbbL), formaldehyde assimilation (SHMT), and nitric oxide reduction (NorQ) shed light on the roles played by GSL viruses in biogeochemical cycles of global relevance.

Incomplete tricarboxylic acid cycle and proton gradient in Pandoravirus massiliensis: is it still a virus?

The discovery of Acanthamoeba polyphaga Mimivirus, the first isolated giant virus of amoeba, challenged the historical hallmarks defining a virus. Giant virion sizes are known to reach up to 2.3 µm, making them visible by optical microscopy. Their large genome sizes of up to 2.5 Mb can encode proteins involved in the translation apparatus. We have investigated possible energy production in Pandoravirus massiliensis. Mitochondrial membrane markers allowed for the detection of a membrane potential in purified virions and this was enhanced by a regulator of the tricarboxylic acid cycle but abolished by the use of a depolarizing agent. Bioinformatics was employed to identify enzymes involved in virion proton gradient generation and this approach revealed that eight putative P. massiliensis proteins exhibited low sequence identities with known cellular enzymes involved in the universal tricarboxylic acid cycle. Further, all eight viral genes were transcribed during replication. The product of one of these genes, ORF132, was cloned and expressed in Escherichia coli, and shown to function as an isocitrate dehydrogenase, a key enzyme of the tricarboxylic acid cycle. Our findings show for the first time that a membrane potential can exist in Pandoraviruses, and this may be related to tricarboxylic acid cycle. The presence of a proton gradient in P. massiliensis makes this virus a form of life for which it is legitimate to ask the question “what is a virus?”.

Characterization and genomic analysis of the first Oceanospirillum phage, vB_OliS_GJ44, representing a novel siphoviral cluster

Background

Marine bacteriophages play key roles in the community structure of microorganisms, biogeochemical cycles, and the mediation of genetic diversity through horizontal gene transfer. Recently, traditional isolation methods, complemented by high-throughput sequencing metagenomics technology, have greatly increased our understanding of the diversity of bacteriophages. Oceanospirillum, within the order Oceanospirillales, are important symbiotic marine bacteria associated with hydrocarbon degradation and algal blooms, especially in polar regions. However, until now there has been no isolate of an Oceanospirillum bacteriophage, and so details of their metagenome has remained unknown.

Results

Here, we reported the first Oceanospirillum phage, vB_OliS_GJ44, which was assembled into a 33,786 bp linear dsDNA genome, which includes abundant tail-related and recombinant proteins. The recombinant module was highly adapted to the host, according to the tetranucleotides correlations. Genomic and morphological analyses identified vB_OliS_GJ44 as a siphovirus, however, due to the distant evolutionary relationship with any other known siphovirus, it is proposed that this virus could be classified as the type phage of a new Oceanospirivirus genus within the Siphoviridae family. vB_OliS_GJ44 showed synteny with six uncultured phages, which supports its representation in uncultured environmental viral contigs from metagenomics. Homologs of several vB_OliS_GJ44 genes have mostly been found in marine metagenomes, suggesting the prevalence of this phage genus in the oceans.

Conclusions

These results describe the first Oceanospirillum phage, vB_OliS_GJ44, that represents a novel viral cluster and exhibits interesting genetic features related to phage–host interactions and evolution. Thus, we propose a new viral genus Oceanospirivirus within the Siphoviridae family to reconcile this cluster, with vB_OliS_GJ44 as a representative member.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07978-4.

Differences in Gut Virome Related to Barrett Esophagus and Esophageal Adenocarcinoma

The relationship between viruses (dominated by bacteriophages or phages) and lower gastrointestinal (GI) tract diseases has been investigated, whereas the relationship between gut bacteriophages and upper GI tract diseases, such as esophageal diseases, which mainly include Barrett’s esophagus (BE) and esophageal adenocarcinoma (EAC), remains poorly described. This study aimed to reveal the gut bacteriophage community and their behavior in the progression of esophageal diseases. In total, we analyzed the gut phage community of sixteen samples from patients with esophageal diseases (six BE patients and four EAC patients) as well as six healthy controls. Differences were found in the community composition of abundant and rare bacteriophages among three groups. In addition, the auxiliary metabolic genes (AMGs) related to bacterial exotoxin and virulence factors such as lipopolysaccharides (LPS) biosynthesis proteins were found to be more abundant in the genome of rare phages from BE and EAC samples compared to the controls. These results suggest that the community composition of gut phages and functional traits encoded by them were different in two stages of esophageal diseases. However, the findings from this study need to be validated with larger sample sizes in the future.

Exploring Viral Diversity in a Gypsum Karst Lake Ecosystem Using Targeted Single-Cell Genomics

Little is known about the diversity and distribution of viruses infecting green sulfur bacteria (GSB) thriving in euxinic (sulfuric and anoxic) habitats, including gypsum karst lake ecosystems. In this study, we used targeted cell sorting combined with single-cell sequencing to gain insights into the gene content and genomic potential of viruses infecting sulfur-oxidizing bacteria Chlorobium clathratiforme, obtained from water samples collected during summer stratification in gypsum karst Lake Kirkilai (Lithuania). In total, 82 viral contigs were bioinformatically identified in 62 single amplified genomes (SAGs) of C. clathratiforme. The majority of viral gene and protein sequences showed little to no similarity with phage sequences in public databases, uncovering the vast diversity of previously undescribed GSB viruses. We observed a high level of lysogenization in the C. clathratiforme population, as 87% SAGs contained intact prophages. Among the thirty identified auxiliary metabolic genes (AMGs), two, thiosulfate sulfurtransferase (TST) and thioredoxin-dependent phosphoadenosine phosphosulfate (PAPS) reductase (cysH), were found to be involved in the oxidation of inorganic sulfur compounds, suggesting that viruses can influence the metabolism and cycling of this essential element. Finally, the analysis of CRISPR spacers retrieved from the consensus C. clathratiforme genome imply persistent and active virus–host interactions for several putative phages prevalent among C. clathratiforme SAGs. Overall, this study provides a glimpse into the diversity of phages associated with naturally occurring and highly abundant sulfur-oxidizing bacteria.

Ecophysiological Features Shape the Distribution of Prophages and CRISPR in Sulfate Reducing Prokaryotes

Sulfate reducing prokaryotes (SRP) are a phylogenetically and physiologically diverse group of microorganisms that use sulfate as an electron acceptor. SRP have long been recognized as key players of the carbon and sulfur cycles, and more recently, they have been identified to play a relevant role as part of syntrophic and symbiotic relations and the human microbiome. Despite their environmental relevance, there is a poor understanding about the prevalence of prophages and CRISPR arrays and how their distribution and dynamic affect the ecological role of SRP. We addressed this question by analyzing the results of a comprehensive survey of prophages and CRISPR in a total of 91 genomes of SRP with several genotypic, phenotypic, and physiological traits, including genome size, cell volume, minimum doubling time, cell wall, and habitat, among others. Our analysis discovered 81 prophages in 51 strains, representing the 56% of the total evaluated strains. Prophages are non-uniformly distributed across the SRP phylogeny, where prophage-rich lineages belonged to Desulfovibrionaceae and Peptococcaceae. Furthermore, our study found 160 CRISPR arrays in 71 SRP, which is more abundant and widely spread than previously expected. Although there is no correlation between presence and abundance of prophages and CRISPR arrays at the strain level, our analysis showed that there is a directly proportional relation between cellular volumes and number of prophages per cell. This result suggests that there is an additional selective pressure for strains with smaller cells to get rid of foreign DNA, such as prophages, but not CRISPR, due to less availability of cellular resources. Analysis of the prophage genes encoding viral structural proteins reported that 44% of SRP prophages are classified as Myoviridae, and comparative analysis showed high level of homology, but not synteny, among prophages belonging to the Family Desulfovibrionaceae. We further recovered viral-like particles and structures that resemble outer membrane vesicles from D. vulgaris str. Hildenborough. The results of this study improved the current understanding of dynamic interactions between prophages and CRISPR with their hosts in both cultured and hitherto-uncultured SRP strains, and how their distribution affects the microbial community dynamics in several sulfidogenic natural and engineered environments.

Extended Evaluation of Viral Diversity in Lake Baikal through Metagenomics

Lake Baikal is a unique oligotrophic freshwater lake with unusually cold conditions and amazing biological diversity. Studies of the lake's viral communities have begun recently, and their full diversity is not elucidated yet. Here, we performed DNA viral metagenomic analysis on integral samples from four different deep-water and shallow stations of the southern and central basins of the lake. There was a strict distinction of viral communities in areas with different environmental conditions. Comparative analysis with other freshwater lakes revealed the highest similarity of Baikal viromes with those of the Asian lakes Soyang and Biwa. Analysis of new data, together with previously published data allowed us to get a deeper insight into the diversity and functional potential of Baikal viruses; however, the true diversity of Baikal viruses in the lake ecosystem remains still unknown. The new metaviromic data will be useful for future studies of viral composition, distribution, and the dynamics associated with global climatic and anthropogenic impacts on this ecosystem.

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.

The Microbial Composition in Circumneutral Thermal Springs from Chignahuapan, Puebla, Mexico Reveals the Presence of Particular Sulfur-Oxidizing Bacterial and Viral Communities

Terrestrial thermal springs are widely distributed globally, and these springs harbor a broad diversity of organisms of biotechnological interest. In Mexico, few studies exploring this kind of environment have been described. In this work, we explore the microbial community in Chignahuapan hot springs, which provides clues to understand these ecosystems' diversity. We assessed the diversity of the microorganism communities in a hot spring environment with a metagenomic shotgun approach. Besides identifying similarities and differences with other ecosystems, we achieved a systematic comparison against 11 metagenomic samples from diverse localities. The Chignahuapan hot springs show a particular prevalence of sulfur-oxidizing bacteria from the genera Rhodococcus, Thermomonas, Thiomonas, Acinetobacter, Sulfurovum, and Bacillus, highlighting those that are different from other recovered bacterial populations in circumneutral hot springs environments around the world. The co-occurrence analysis of the bacteria and viruses in these environments revealed that within the Rhodococcus, Thiomonas, Thermonas, and Bacillus genera, the Chignahuapan samples have specific species of bacteria with a particular abundance, such as Rhodococcus erytropholis. The viruses in the circumneutral hot springs present bacteriophages within the order Caudovirales (Siphoviridae, Myoviridae, and Podoviridae), but the family of Herelleviridae was the most abundant in Chignahuapan samples. Furthermore, viral auxiliary metabolic genes were identified, many of which contribute mainly to the metabolism of cofactors and vitamins as well as carbohydrate metabolism. Nevertheless, the viruses and bacteria present in the circumneutral environments contribute to the sulfur cycle. This work represents an exhaustive characterization of a community structure in samples collected from hot springs in Mexico and opens opportunities to identify organisms of biotechnological interest.

Single-virus genomics and beyond

Viruses are extremely diverse and modulate important biological and ecological processes globally. However, much of viral diversity remains uncultured and yet to be discovered. Several powerful culture-independent tools, in particular metagenomics, have substantially advanced virus discovery. Among those tools is single-virus genomics, which yields sequenced reference genomes from individual sorted virus particles without the need for cultivation. This new method complements virus culturing and metagenomic approaches and its advantages include targeted investigation of specific virus groups and investigation of genomic microdiversity within viral populations. In this Review, we provide a brief history of single-virus genomics, outline how this emergent method has facilitated advances in virus ecology and discuss its current limitations and future potential. Finally, we address how this method may synergistically intersect with other single-virus and single-cell approaches.

Viral ecogenomics across the Porifera

BACKGROUND

Viruses directly affect the most important biological processes in the ocean via their regulation of prokaryotic and eukaryotic populations. Marine sponges form stable symbiotic partnerships with a wide diversity of microorganisms and this high symbiont complexity makes them an ideal model for studying viral ecology. Here, we used morphological and molecular approaches to illuminate the diversity and function of viruses inhabiting nine sponge species from the Great Barrier Reef and seven from the Red Sea.

RESULTS

Viromic sequencing revealed host-specific and site-specific patterns in the viral assemblages, with all sponge species dominated by the bacteriophage order Caudovirales but also containing variable representation from the nucleocytoplasmic large DNA virus families Mimiviridae, Marseilleviridae, Phycodnaviridae, Ascoviridae, Iridoviridae, Asfarviridae and Poxviridae. Whilst core viral functions related to replication, infection and structure were largely consistent across the sponge viromes, functional profiles varied significantly between species and sites largely due to differential representation of putative auxiliary metabolic genes (AMGs) and accessory genes, including those associated with herbicide resistance, heavy metal resistance and nylon degradation. Furthermore, putative AMGs varied with the composition and abundance of the sponge-associated microbiome. For instance, genes associated with antimicrobial activity were enriched in low microbial abundance sponges, genes associated with nitrogen metabolism were enriched in high microbial abundance sponges and genes related to cellulose biosynthesis were enriched in species that host photosynthetic symbionts.

CONCLUSIONS

Our results highlight the diverse functional roles that viruses can play in marine sponges and are consistent with our current understanding of sponge ecology. Differential representation of putative viral AMGs and accessory genes across sponge species illustrate the diverse suite of beneficial roles viruses can play in the functional ecology of these complex reef holobionts. Video Abstract.

Genomic analysis of Synechococcus phage S-B43 and its adaption to the coastal environment

Synechococcus dominate picocyanobacterial communities in coastal environments. However, only a few Synechococcus phages have been described from the coastal seas of the Northwest Pacific Ocean. Here a new Synechococcus phage, S-B43 was isolated from the Bohai Sea, a semi-closed coastal sea of the Northwest Pacific Ocean. S-B43 is a member of Myoviridae, containing 275 predicted open reading frames. Fourteen auxiliary metabolic genes (AMG) were identified from the genome of S-B43, including five photosynthetic associated genes and several AMGs related to its adaption to the high turbidity and eutrophic coastal environment with a low ratio of phosphorus to nitrogen (HNLP). The occurrences of 31 AMGs among 34 cyanophage genomes indicates that AMGs zwf, gnd, speD, petF and those coding for FECH and thioredoxin were more common in coastal areas than in the open ocean and AMGs pebS and ho1 were more prevalent in the open ocean. The occurrence of cyanophage AMGs in different environments might be a reflection of the environmental adaption of their hosts. This study contributes to our understanding of the interactions between cyanobacteria and cyanophages and their environmental adaption to the coastal environment.