Depth-stratified functional and taxonomic niche specialization in the 'core' and 'flexible' Pacific Ocean Virome

A global atlas of soil viruses reveals unexplored biodiversity and potential biogeochemical impacts

Historically neglected by microbial ecologists, soil viruses are now thought to be critical to global biogeochemical cycles. However, our understanding of their global distribution, activities and interactions with the soil microbiome remains limited. Here we present the Global Soil Virus Atlas, a comprehensive dataset compiled from 2,953 previously sequenced soil metagenomes and composed of 616,935 uncultivated viral genomes and 38,508 unique viral operational taxonomic units. Rarefaction curves from the Global Soil Virus Atlas indicate that most soil viral diversity remains unexplored, further underscored by high spatial turnover and low rates of shared viral operational taxonomic units across samples. By examining genes associated with biogeochemical functions, we also demonstrate the viral potential to impact soil carbon and nutrient cycling. This study represents an extensive characterization of soil viral diversity and provides a foundation for developing testable hypotheses regarding the role of the virosphere in the soil microbiome and global biogeochemistry.

Environmental viromes reveal the global distribution signatures of deep-sea DNA viruses

INTRODUCTION

Viruses are abundant and ecologically significant in marine ecosystems. However, the virome of deep-sea sediments is not extensively investigated.

OBJECTIVES

To explore the distribution pattern of deep-sea viruses on a global scale, the viromes of DNA viruses isolated from 138 sediments of 5 deep-sea ecosystems were characterized.

METHODS

The viral particles were purified from each sediment sample. Then the viral DNAs were extracted and subjected to viral metagenomic analysis.

RESULTS

Here, we constructed a global deep-sea environmental virome dataset by analyzing the viral DNA of 138 sediment samples. A total of 347,737 viral operational taxonomic units (vOTUs) were identified, of which 84.94% were hitherto unknown, indicating that deep sea was a reservoir of novel DNA viruses. Furthermore, circular viral genome analysis revealed 98,581 complete genomes. The classified vOTUs included eukaryotic (44.55%) and prokaryotic (25.75%) viruses, and were taxonomically assigned to 63 viral families. The composition and abundance of the deep-sea sediment viromes were dependent on the deep-sea ecosystem as opposed to geographical region. Further analysis revealed that the viral community differentiation in different deep-sea ecosystems was driven by the virus-mediated energy metabolism.

CONCLUSION

Our findings showed that deep-sea ecosystems are a reservoir of novel DNA viruses and the viral community is shaped by the environmental characteristics of deep-sea ecosystems, thus presenting critical information for determining the ecological significance of viruses in global deep-sea ecosystems.

Viral potential to modulate microbial methane metabolism varies by habitat

Methane is a potent greenhouse gas contributing to global warming. Microorganisms largely drive the biogeochemical cycling of methane, yet little is known about viral contributions to methane metabolism (MM). We analyzed 982 publicly available metagenomes from host-associated and environmental habitats containing microbial MM genes, expanding the known MM auxiliary metabolic genes (AMGs) from three to 24, including seven genes exclusive to MM pathways. These AMGs are recovered on 911 viral contigs predicted to infect 14 prokaryotic phyla including Halobacteriota, Methanobacteriota, and Thermoproteota. Of those 24, most were encoded by viruses from rumen (16/24), with substantially fewer by viruses from environmental habitats (0-7/24). To search for additional MM AMGs from an environmental habitat, we generate metagenomes from methane-rich sediments in Vrana Lake, Croatia. Therein, we find diverse viral communities, with most viruses predicted to infect methanogens and methanotrophs and some encoding 13 AMGs that can modulate host metabolisms. However, none of these AMGs directly participate in MM pathways. Together these findings suggest that the extent to which viruses use AMGs to modulate host metabolic processes (e.g., MM) varies depending on the ecological properties of the habitat in which they dwell and is not always predictable by habitat biogeochemical properties.

Eco-evolutionary dynamics of gut phageome in wild gibbons (Hoolock tianxing) with seasonal diet variations

It has been extensively studied that the gut microbiome provides animals flexibility to adapt to food variability. Yet, how gut phageome responds to diet variation of wild animals remains unexplored. Here, we analyze the eco-evolutionary dynamics of gut phageome in six wild gibbons (Hoolock tianxing) by collecting individually-resolved fresh fecal samples and parallel feeding behavior data for 15 consecutive months. Application of complementary viral and microbial metagenomics recovers 39,198 virulent and temperate phage genomes from the feces. Hierarchical cluster analyses show remarkable seasonal diet variations in gibbons. From high-fruit to high-leaf feeding period, the abundances of phage populations are seasonally fluctuated, especially driven by the increased abundance of virulent phages that kill the Lachnospiraceae hosts, and a decreased abundance of temperate phages that piggyback the Bacteroidaceae hosts. Functional profiling reveals an enrichment through horizontal gene transfers of toxin-antitoxin genes on temperate phage genomes in high-leaf season, potentially conferring benefits to their prokaryotic hosts. The phage-host ecological dynamics are driven by the coevolutionary processes which select for tail fiber and DNA primase genes on virulent and temperate phage genomes, respectively. Our results highlight complex phageome-microbiome interactions as a key feature of the gibbon gut microbial ecosystem responding to the seasonal diet.

High niche specificity and host genetic diversity of groundwater viruses

Viruses are key members of microbial communities that exert control over host abundance and metabolism, thereby influencing ecosystem processes and biogeochemical cycles. Aquifers are known to host taxonomically diverse microbial life, yet little is known about viruses infecting groundwater microbial communities. Here, we analysed 16 metagenomes from a broad range of groundwater physicochemistries. We recovered 1571 viral genomes that clustered into 468 high-quality viral operational taxonomic units. At least 15% were observed to be transcriptionally active, although lysis was likely constrained by the resource-limited groundwater environment. Most were unclassified (95%), and the remaining 5% were Caudoviricetes. Comparisons with viruses inhabiting other aquifers revealed no shared species, indicating substantial unexplored viral diversity. In silico predictions linked 22.4% of the viruses to microbial host populations, including to ultra-small prokaryotes, such as Patescibacteria and Nanoarchaeota. Many predicted hosts were associated with the biogeochemical cycling of carbon, nitrogen, and sulfur. Metabolic predictions revealed the presence of 205 putative auxiliary metabolic genes, involved in diverse processes associated with the utilization of the host's intracellular resources for biosynthesis and transformation reactions, including those involved in nucleotide sugar, glycan, cofactor, and vitamin metabolism. Viruses, prokaryotes overall, and predicted prokaryotic hosts exhibited narrow spatial distributions, and relative abundance correlations with the same groundwater parameters (e.g. dissolved oxygen, nitrate, and iron), consistent with host control over viral distributions. Results provide insights into underexplored groundwater viruses, and indicate the large extent to which viruses may manipulate microbial communities and biogeochemistry in the terrestrial subsurface.

Genomic and transcriptomic insights into complex virus-prokaryote interactions in marine biofilms

Marine biofilms are complex communities of microorganisms that play a crucial ecological role in oceans. Although prokaryotes are the dominant members of these biofilms, little is known about their interactions with viruses. By analysing publicly available and newly sequenced metagenomic data, we identified 2446 virus-prokaryote connections in 84 marine biofilms. Most of these connections were between the bacteriophages in the Uroviricota phylum and the bacteria of Proteobacteria, Cyanobacteria and Bacteroidota. The network of virus-host pairs is complex; a single virus can infect multiple prokaryotic populations or a single prokaryote is susceptible to several viral populations. Analysis of genomes of paired prokaryotes and viruses revealed the presence of 425 putative auxiliary metabolic genes (AMGs), 239 viral genes related to restriction-modification (RM) systems and 38,538 prokaryotic anti-viral defence-related genes involved in 15 defence systems. Transcriptomic evidence from newly established biofilms revealed the expression of viral genes, including AMGs and RM, and prokaryotic defence systems, indicating the active interplay between viruses and prokaryotes. A comparison between biofilms and seawater showed that biofilm prokaryotes have more abundant defence genes than seawater prokaryotes, and the defence gene composition differs between biofilms and the surrounding seawater. Overall, our study unveiled active viruses in natural biofilms and their complex interplay with prokaryotes, which may result in the blooming of defence strategists in biofilms. The detachment of bloomed defence strategists may reduce the infectivity of viruses in seawater and result in the emergence of a novel role of marine biofilms.

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.

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.

Virome reveals effect of Ulva prolifera green tide on the structural and functional profiles of virus communities in coastal environments

Viruses are widely distributed in marine environments, where they influence the transformation of matter and energy by modulating host metabolism. Driven by eutrophication, green tides are a rising concern in Chinese coastal areas, and are a serious ecological disaster that negatively affects coastal ecosystems and disrupts biogeochemical cycles. Although the composition of bacterial communities in green algae has been investigated, the diversity and roles of viruses in green algal blooms are largely unexplored. Therefore, the diversity, abundance, lifestyle, and metabolic potential of viruses in a natural bloom in Qingdao coastal area were investigated at three different stages (pre-bloom, during-bloom, and post-bloom) by metagenomics analysis. The dsDNA viruses, Siphoviridae, Myoviridae, Podoviridae, and Phycodnaviridae, were found to dominate the viral community. The viral dynamics exhibited distinct temporal patterns across different stages. The composition of the viral community varied during the bloom, especially in populations with low abundance. The lytic cycle was most predominant, and the abundance of lytic viruses increased slightly in the post-bloom stage. The diversity and richness of the viral communities varied distinctly during the green tide, and the post-bloom stage favored viral diversity and richness. The total organic carbon, dissolved oxygen, NO^3-, NO^2-, PO₄^3-, chlorophyll-a contents, and temperature variably co-influenced the viral communities. The primary hosts included bacteria, algae, and other microplankton. Network analysis revealed the closer links between the viral communities as the bloom progressed. Functional prediction revealed that the viruses possibly influenced the biodegradation of microbial hydrocarbons and carbon by metabolic augmentation via auxiliary metabolic genes. The composition, structure, metabolic potential, and interaction taxonomy of the viromes differed significantly across the different stages of the green tide. The study demonstrated that the ecological event shaped the viral communities during algal bloom, and the viral communities played a significant role in phycospheric microecology.

Integrative omics framework for characterization of coral reef ecosystems from the Tara Pacific expedition

Coral reef science is a fast-growing field propelled by the need to better understand coral health and resilience to devise strategies to slow reef loss resulting from environmental stresses. Key to coral resilience are the symbiotic interactions established within a complex holobiont, i.e. the multipartite assemblages comprising the coral host organism, endosymbiotic dinoflagellates, bacteria, archaea, fungi, and viruses. Tara Pacific is an ambitious project built upon the experience of previous Tara Oceans expeditions, and leveraging state-of-the-art sequencing technologies and analyses to dissect the biodiversity and biocomplexity of the coral holobiont screened across most archipelagos spread throughout the entire Pacific Ocean. Here we detail the Tara Pacific workflow for multi-omics data generation, from sample handling to nucleotide sequence data generation and deposition. This unique multidimensional framework also includes a large amount of concomitant metadata collected side-by-side that provide new assessments of coral reef biodiversity including micro-biodiversity and shape future investigations of coral reef dynamics and their fate in the Anthropocene.

Water mass age structures the auxiliary metabolic gene content of free-living and particle-attached deep ocean viral communities

BACKGROUND

Viruses play important roles in the ocean's biogeochemical cycles. Yet, deep ocean viruses are one of the most under-explored fractions of the global biosphere. Little is known about the environmental factors that control the composition and functioning of their communities or how they interact with their free-living or particle-attached microbial hosts.

RESULTS

We analysed 58 viral communities associated with size-fractionated free-living (0.2-0.8 μm) and particle-attached (0.8-20 μm) cellular metagenomes from bathypelagic (2150-4018 m deep) microbiomes obtained during the Malaspina expedition. These metagenomes yielded 6631 viral sequences, 91% of which were novel, and 67 represented high-quality genomes. Taxonomic classification assigned 53% of the viral sequences to families of tailed viruses from the order Caudovirales. Computational host prediction associated 886 viral sequences to dominant members of the deep ocean microbiome, such as Alphaproteobacteria (284), Gammaproteobacteria (241), SAR324 (23), Marinisomatota (39), and Chloroflexota (61). Free-living and particle-attached viral communities had markedly distinct taxonomic composition, host prevalence, and auxiliary metabolic gene content, which led to the discovery of novel viral-encoded metabolic genes involved in the folate and nucleotide metabolisms. Water mass age emerged as an important factor driving viral community composition. We postulated this was due to changes in quality and concentration of dissolved organic matter acting on the host communities, leading to an increase of viral auxiliary metabolic genes associated with energy metabolism among older water masses.

CONCLUSIONS

These results shed light on the mechanisms by which environmental gradients of deep ocean ecosystems structure the composition and functioning of free-living and particle-attached viral communities. Video Abstract.

Viruses of sulfur oxidizing phototrophs encode genes for pigment, carbon, and sulfur metabolisms

Viral infections modulate bacterial metabolism and ecology. Here, we investigated the hypothesis that viruses influence the ecology of purple and green sulfur bacteria in anoxic and sulfidic lakes, analogs of euxinic oceans in the geologic past. By screening metagenomes from lake sediments and water column, in addition to publicly-available genomes of cultured purple and green sulfur bacteria, we identified almost 300 high and medium-quality viral genomes. Viruses carrying the gene psbA, encoding the small subunit of photosystem II protein D1, were ubiquitous, suggesting viral interference with the light reactions of sulfur oxidizing autotrophs. Viruses predicted to infect these autotrophs also encoded auxiliary metabolic genes for reductive sulfur assimilation as cysteine, pigment production, and carbon fixation. These observations show that viruses have the genomic potential to modulate the production of metabolic markers of phototrophic sulfur bacteria that are used to identify photic zone euxinia in the geologic past.

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.

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.

Potential Auxiliary Metabolic Capabilities and Activities Reveal Biochemical Impacts of Viruses in Municipal Wastewater Treatment Plants

Viruses influence biogeochemical cycles in oceans, freshwater, soil, and human gut through infection and by modulating virocell metabolism through virus-encoded auxiliary metabolic genes (vAMGs). However, the geographical distribution, potential metabolic function, and engineering significance of vAMGs in wastewater treatment plants (WWTPs) remain to be explored. Here, 752 single-contig viral genomes with high confidence, 510 of which belonged to Caudovirales, were recovered from the activated sludge metagenomes of 32 geographically distributed WWTPs. A total of 101 vAMGs involved in various metabolic pathways were identified, the most common of which were the queuosine biosynthesis genes folE, queD, and queE and the sulfur metabolism gene cysH. Phylogenetic analysis and virus-host relationship prediction revealed the probable evolutionary histories of vAMGs involved in carbon (acpP and prsA), nitrogen (amoC), sulfur (cysH), and phosphate (phoH) metabolism, which potentially mediate microbial carbon and nutrient cycling. Notably, 11 of the 38 (28.3%) vAMGs identified in the metagenomes with corresponding metatranscriptomes were transcriptionally expressed, implying an active functional state. This meta-analysis provides the first broad catalog of vAMGs in municipal WWTPs and how they may assist in the basic physiological reactions of their microbial hosts or nutrient cycling in the WWTPs, and therefore, may have important effects on the engineering of wastewater treatment processes.

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.

Microbial and phage communities as well as their interaction in PO saponification wastewater treatment systems

Viruses or phages were considered affecting microbial community composition, metabolic process, and biogeochemical cycles. However, phage communities and their potential associations with microbial community are not well understood in the activated sludge (AS) of wastewater treatment plants (WWTPs). In this study, we explored the interactions between phages and microbial community by using propylene oxide (PO) saponification WWTPs as an example. Bacterial, eukaryal and archaeal communities were investigated and 34 phage contigs (>10 kb) were recovered from PO saponification WWTPs. At least 3 complete phage genomes were assembled. In all 34 phages, 21 of them have been predicted to their host. The association network analysis showed that abundant phages were associated with abundant microorganisms. This result conformed to Kill-the-Winner model. Notably, 45 auxiliary metabolic genes (AMGs) were identified from phage genomes (including small contig fragments). They influenced bacterial metabolism through facilitating phages replication and avoiding host death. Collectively, our results suggested that phage community affect microbial community and metabolic pathways by killing their hosts and AMGs transfer in AS of PO saponification WWTPs.

Over two decades of research on the marine RNA virosphere

RNA viruses (realm: Riboviria), including RNA phages and eukaryote-infecting RNA viruses, are essential components of marine ecosystems. A large number of marine RNA viruses have been discovered in the last two decades because of the rapid development of next-generation sequencing (NGS) technology. Indeed, the combination of NGS and state-of-the-art meta-omics methods (viromics, the study of all viruses in a specific environment) has led to a fundamental understanding of the taxonomy and genetic diversity of RNA viruses in the sea, suggesting the complex ecological roles played by RNA viruses in this complex ecosystem. Furthermore, comparisons of viromes in the context of highly variable marine niches reveal the biogeographic patterns and ecological impact of marine RNA viruses, whose role in global ecology is becoming increasingly clearer. In this review, we summarize the characteristics of the global marine RNA virosphere and outline the taxonomic hierarchy of RNA viruses with a specific focus on their ancient evolutionary history. We also review the development of methodology and the major progress resulting from its applications in RNA viromics. The aim of this review is not only to provide an in-depth understanding of multifaceted aspects of marine RNA viruses, but to offer future perspectives on developing a better methodology for discovery, and exploring the evolutionary origin and major ecological significance of marine RNA virosphere.

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.

The microbiome and its association with antibiotic resistance genes in the hadal biosphere at the Yap Trench

The hadal biosphere, the deepest part of the ocean, is known as the least-explored aquatic environment and hosts taxonomically diverse microbial communities. However, the microbiome and its association with antibiotic resistance genes (ARGs) in the hadal ecosystem remain unknown. Here, we profiled the microbiome diversity and ARG occurrence in seawater and sediments of the Yap Trench (YT) using metagenomic sequencing. Within the prokaryote (bacteria and archaea) lineages, the main components of bacteria were Gammaproteobacteria (77.76 %), Firmicutes (8.36 %), and Alphaproteobacteria (2.25 %), whereas the major components of archaea were Nitrososphaeria (6.51 %), Nanoarchaeia (0.42 %), and Thermoplasmata (0.25 %), respectively. Taxonomy of viral contigs showed that the classified viral communities in YT seawater and sediments were dominated by Podoviridae (45.96 %), Siphoviridae (29.41 %), and Myoviridae (24.63 %). A large majority of viral contigs remained uncharacterized and exhibited endemicity. A total of 48 ARGs encoding resistance to 12 antibiotic classes were identified and their hosts were bacteria and viruses. Novel ARG subtypes mexF_YTV-1, mexF_YTV-2, mexF_YTV-3, vanR_YTV-1, vanS_YTV-1 (carried by unclassified viruses), and bacA_YTB-1 (carried by phylum Firmicutes) were detected in seawater samples. Overall, our findings imply that the hadal environment of the YT is a repository of viral and ARG diversity.

Viruses in Subsurface Environments

Over the past 20 years, our knowledge of virus diversity and abundance in subsurface environments has expanded dramatically through application of quantitative metagenomic approaches. In most subsurface environments, viral diversity and abundance rival viral diversity and abundance observed in surface environments. Most of these viruses are uncharacterized in terms of their hosts and replication cycles. Analysis of accessory metabolic genes encoded by subsurface viruses indicates that they evolved to replicate within the unique features of their environments. The key question remains: What role do these viruses play in the ecology and evolution of the environments in which they replicate? Undoubtedly, as more virologists examine the role of viruses in subsurface environments, new insights will emerge.

Into the Dark: Exploring the Deep Ocean with Single-Virus Genomics

Single-virus genomics (SVGs) has been successfully applied to ocean surface samples allowing the discovery of widespread dominant viruses overlooked for years by metagenomics, such as the uncultured virus vSAG 37-F6 infecting the ubiquitous Pelagibacter spp. In SVGs, one uncultured virus at a time is sorted from the environmental sample, whole-genome amplified, and sequenced. Here, we have applied SVGs to deep-ocean samples (200–4000 m depth) from global Malaspina and MEDIMAX expeditions, demonstrating the feasibility of this method in deep-ocean samples. A total of 1328 virus-like particles were sorted from the North Atlantic Ocean, the deep Mediterranean Sea, and the Pacific Ocean oxygen minimum zone (OMZ). For this proof of concept, sixty single viruses were selected at random for sequencing. Genome annotation identified 27 of these genomes as bona fide viruses, and detected three auxiliary metabolic genes involved in nucleotide biosynthesis and sugar metabolism. Massive protein profile analysis confirmed that these viruses represented novel viral groups not present in databases. Although they were not previously assembled by viromics, global fragment recruitment analysis showed a conserved profile of relative abundance of these viruses in all analyzed samples spanning different oceans. Altogether, these results reveal the feasibility in using SVGs in this vast environment to unveil the genomes of relevant viruses.

Conceptualizing the 3D niche and vertical space use

Spatial partitioning in ecological communities has predominantly been described in two dimensions, yet habitat is complex and 3D. Complex space use mediates community structure and interaction strength by expanding spatial, temporal, and dietary dimensions. Vertical stratification of resources provides opportunities for novel specializations, creating a 3D niche. Competition and predation are mediated by 3D space use, as individuals use the vertical axis to access prey, flee predators, or avoid competitors. The 3D niche is important for long-term conservation strategies as species must navigate tradeoffs in habitat use between strata-specific threats and suboptimal habitat patches. Ultimately, elucidating the 3D niche has implications for protected area management and corridor design that directly influence species persistence and ecosystem function in a rapidly changing world.

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.

Infection strategy and biogeography distinguish cosmopolitan groups of marine jumbo bacteriophages

Recent research has underscored the immense diversity and key biogeochemical roles of large DNA viruses in the ocean. Although they are important constituents of marine ecosystems, it is sometimes difficult to detect these viruses due to their large size and complex genomes. This is true for "jumbo" bacteriophages, which have genome sizes >200 kbp and large capsids reaching up to 0.45 µm in diameter. In this study, we sought to assess the genomic diversity and distribution of these bacteriophages in the ocean by generating and analyzing jumbo phage genomes from metagenomes. We recover 85 marine jumbo phages that ranged in size from 201 to 498 kilobases, and we examine their genetic similarities and biogeography together with a reference database of marine jumbo phage genomes. By analyzing Tara Oceans metagenomic data, we show that although most jumbo phages can be detected in a range of different size fractions, 17 of our bins tend to be found in those greater than 0.22 µm, potentially due to their large size. Our network-based analysis of gene-sharing patterns reveals that jumbo bacteriophages belong to five genome clusters that are typified by diverse replication strategies, genomic repertoires, and potential host ranges. Our analysis of jumbo phage distributions in the ocean reveals that depth is a major factor shaping their biogeography, with some phage genome clusters occurring preferentially in either surface or mesopelagic waters, respectively. Taken together, our findings indicate that jumbo phages are widespread community members in the ocean with complex genomic repertoires and ecological impacts that warrant further targeted investigation.

Organochlorine contamination enriches virus-encoded metabolism and pesticide degradation associated auxiliary genes in soil microbiomes

Viruses significantly influence local and global biogeochemical cycles and help bacteria to survive in different environments by encoding various auxiliary metabolic genes (AMGs) associated with energy acquisition, stress tolerance and degradation of xenobiotics. Here we studied whether bacterial (dsDNA) virus encoded AMGs are enriched in organochlorine pesticide (OCP) contaminated soil in China and if viral AMGs include genes linked to OCP biodegradation. Using metagenomics, we found that OCP-contaminated soils displayed a lower bacterial, but higher diversity of viruses that harbored a higher relative abundance of AMGs linked to pesticide degradation and metabolism. Furthermore, the diversity and relative abundance of AMGs significantly increased along with the severity of pesticide contamination, and several biodegradation genes were identified bioinformatically in viral metagenomes. Functional assays were conducted to experimentally demonstrate that virus-encoded L-2-haloacid dehalogenase gene (L-DEX) is responsible for the degradation of L-2-haloacid pesticide precursors, improving bacterial growth at sub-inhibitory pesticide concentrations. Taken together, these results demonstrate that virus-encoded AMGs are linked to bacterial metabolism and biodegradation, being more abundant and diverse in soils contaminated with pesticides. Moreover, our findings highlight the importance of virus-encoded accessory genes for bacterial ecology in stressful environments, providing a novel avenue for using viruses in the bioremediation of contaminated soils.

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.

Diversity of Pseudomonas aeruginosa Temperate Phages

Modern sequencing technologies have provided insight into the genetic diversity of numerous species, including the human pathogen Pseudomonas aeruginosa. Bacterial genomes often harbor bacteriophage genomes (prophages), which can account for upwards of 20% of the genome. Prior studies have found P. aeruginosa prophages that contribute to their host’s pathogenicity and fitness. These advantages come in many different forms, including the production of toxins, promotion of biofilm formation, and displacement of other P. aeruginosa strains. While several different genera and species of P. aeruginosa prophages have been studied, there has not been a comprehensive study of the overall diversity of P. aeruginosa-infecting prophages. Here, we present the results of just such an analysis. A total of 6,852 high-confidence prophages were identified from 5,383 P. aeruginosa genomes from strains isolated from the human body and other environments. In total, 3,201 unique prophage sequences were identified. While 53.1% of these prophage sequences displayed sequence similarity to publicly available phage genomes, novel and highly mosaic prophages were discovered. Among these prophages, there is extensive diversity, including diversity within the functionally conserved integrase and C repressor coding regions, two genes responsible for prophage entering and persisting through the lysogenic life cycle. Analysis of integrase, C repressor, and terminase coding regions revealed extensive reassortment among P. aeruginosa prophages. This catalog of P. aeruginosa prophages provides a resource for future studies into the evolution of the species.

IMPORTANCE Prophages play a critical role in the evolution of their host species and can also contribute to the virulence and fitness of pathogenic species. Here, we conducted a comprehensive investigation of prophage sequences from 5,383 publicly available Pseudomonas aeruginosa genomes from human as well as environmental isolates. We identified a diverse population of prophages, including tailed phages, inoviruses, and microviruses; 46.9% of the prophage sequences found share no significant sequence similarity with characterized phages, representing a vast array of novel P. aeruginosa-infecting phages. Our investigation into these prophages found substantial evidence of reassortment. In producing this, the first catalog of P. aeruginosa prophages, we uncovered both novel prophages as well as genetic content that have yet to be explored.

Cyanophage Distribution Across European Lakes of the Temperate-Humid Continental Climate Zone Assessed Using PCR-Based Genetic Markers

Studies of the diversity and distribution of freshwater cyanophages are generally limited to the small geographical areas, in many cases including only one or few lakes. Data from dozens of various lakes distributed at a larger distance are necessary to understand their spatial distribution and sensitivity to biotic and abiotic factors. Thus, the objective of this study was to analyze the diversity and distribution of cyanophages within the infected cells using marker genes (psbA, nblA, and g91) in 21 Polish and Lithuanian lakes. Physicochemical factors that might be related to them were also analyzed. The results demonstrated that genetic markers representing cyanophages were observed in most lakes studied. The frequently detected gene was psbA with 88% of cyanophage-positive samples, while nblA and g91 were found in approximately 50% of lakes. The DNA sequence analyses for each gene demonstrated low variability between them, although the psbA sequences branched within the larger cluster of marine Synechoccocuss counterparts. The principal component analysis allowed to identify significant variation between the lakes that presented high and low cyanobacterial biomass. The lakes with high cyanobacterial biomass were further separated by country and the different diversity of cyanobacteria species, particularly Planktothrix agardhii, was dominant in the Polish lakes and Planktolyngbya limnetica in the Lithuanian lakes. The total phosphorous and the presence of cyanophage genes psbA and nblA were the most important factors that allowed differentiation for the Polish lakes, while the pH and the genes g91 and nblA for the Lithuanian lakes.

Flavobacterium flabelliforme sp. nov. and Flavobacterium geliluteum sp. nov., Two Multidrug-Resistant Psychrotrophic Species Isolated From Antarctica

Despite unfavorable Antarctic conditions, such as cold temperatures, freeze-thaw cycles, high ultraviolet radiation, dryness and lack of nutrients, microorganisms were able to adapt and surprisingly thrive in this environment. In this study, eight cold-adapted Flavobacterium strains isolated from a remote Antarctic island, James Ross Island, were studied using a polyphasic taxonomic approach to determine their taxonomic position. Phylogenetic analyses based on the 16S rRNA gene and 92 core genes clearly showed that these strains formed two distinct phylogenetic clusters comprising three and five strains, with average nucleotide identities significantly below 90% between both proposed species as well as between their closest phylogenetic relatives. Phenotyping revealed a unique pattern of biochemical and physiological characteristics enabling differentiation from the closest phylogenetically related Flavobacterium spp. Chemotaxonomic analyses showed that type strains P4023^T and P7388^T were characterized by the major polyamine sym-homospermidine and a quinone system containing predominantly menaquinone MK-6. In the polar lipid profile phosphatidylethanolamine, an ornithine lipid and two unidentified lipids lacking a functional group were detected as major lipids. These characteristics along with fatty acid profiles confirmed that these species belong to the genus Flavobacterium. Thorough genomic analysis revealed the presence of numerous cold-inducible or cold-adaptation associated genes, such as cold-shock proteins, proteorhodopsin, carotenoid biosynthetic genes or oxidative-stress response genes. Genomes of type strains surprisingly harbored multiple prophages, with many of them predicted to be active. Genome-mining identified biosynthetic gene clusters in type strain genomes with a majority not matching any known clusters which supports further exploratory research possibilities involving these psychrotrophic bacteria. Antibiotic susceptibility testing revealed a pattern of multidrug-resistant phenotypes that were correlated with in silico antibiotic resistance prediction. Interestingly, while typical resistance finder tools failed to detect genes responsible for antibiotic resistance, genomic prediction confirmed a multidrug-resistant profile and suggested even broader resistance than tested. Results of this study confirmed and thoroughly characterized two novel psychrotrophic Flavobacterium species, for which the names Flavobacterium flabelliforme sp. nov. and Flavobacterium geliluteum sp. nov. are proposed.

Diversity and distribution of viruses inhabiting the deepest ocean on Earth

As the most abundant biological entities on the planet, viruses significantly influence the overall functioning of marine ecosystems. The abundance, distribution, and biodiversity of viral communities in the upper ocean have been relatively well studied, but our understanding of viruses in the hadal biosphere remains poor. Here, we established the oceanic trench viral genome dataset (OTVGD) by analysing 19 microbial metagenomes derived from seawater and sediment samples of the Mariana, Yap, and Kermadec Trenches. The trench viral communities harbored remarkably high novelty, and they were predicted to infect ecologically important microbial clades, including Thaumarchaeota and Oleibacter. Significant inter-trench and intra-trench exchange of viral communities was proposed. Moreover, viral communities in different habitats (seawater/sediment and depth-stratified ocean zones) exhibited distinct niche-dependent distribution patterns and genomic properties. Notably, microbes and viruses in the hadopelagic seawater seemed to preferably adopt lysogenic lifestyles compared to those in the upper ocean. Furthermore, niche-specific auxiliary metabolic genes were identified in the hadal viral genomes, and a novel viral D-amino acid oxidase was functionally and phylogenetically characterized, suggesting the contribution of these genes in the utilization of refractory organic matter. Together, these findings highlight the genomic novelty, dynamic movement, and environment-driven diversification of viral communities in oceanic trenches, and suggest that viruses may influence the hadal ecosystem by reprogramming the metabolism of their hosts and modulating the community of keystone microbes.

Deep sea sediments associated with cold seeps are a subsurface reservoir of viral diversity

In marine ecosystems, viruses exert control on the composition and metabolism of microbial communities, influencing overall biogeochemical cycling. Deep sea sediments associated with cold seeps are known to host taxonomically diverse microbial communities, but little is known about viruses infecting these microorganisms. Here, we probed metagenomes from seven geographically diverse cold seeps across global oceans to assess viral diversity, virus-host interaction, and virus-encoded auxiliary metabolic genes (AMGs). Gene-sharing network comparisons with viruses inhabiting other ecosystems reveal that cold seep sediments harbour considerable unexplored viral diversity. Most cold seep viruses display high degrees of endemism with seep fluid flux being one of the main drivers of viral community composition. In silico predictions linked 14.2% of the viruses to microbial host populations with many belonging to poorly understood candidate bacterial and archaeal phyla. Lysis was predicted to be a predominant viral lifestyle based on lineage-specific virus/host abundance ratios. Metabolic predictions of prokaryotic host genomes and viral AMGs suggest that viruses influence microbial hydrocarbon biodegradation at cold seeps, as well as other carbon, sulfur and nitrogen cycling via virus-induced mortality and/or metabolic augmentation. Overall, these findings reveal the global diversity and biogeography of cold seep viruses and indicate how viruses may manipulate seep microbial ecology and biogeochemistry.

Virus-associated organosulfur metabolism in human and environmental systems

Viruses influence the fate of nutrients and human health by killing microorganisms and altering metabolic processes. Organosulfur metabolism and biologically derived hydrogen sulfide play dynamic roles in manifestation of diseases, infrastructure degradation, and essential biological processes. Although microbial organosulfur metabolism is well studied, the role of viruses in organosulfur metabolism is unknown. Here, we report the discovery of 39 gene families involved in organosulfur metabolism encoded by 3,749 viruses from diverse ecosystems, including human microbiomes. The viruses infect organisms from all three domains of life. Six gene families encode for enzymes that degrade organosulfur compounds into sulfide, whereas others manipulate organosulfur compounds and may influence sulfide production. We show that viral metabolic genes encode key enzymatic domains, are translated into protein, and are maintained after recombination, and sulfide provides a fitness advantage to viruses. Our results reveal viruses as drivers of organosulfur metabolism with important implications for human and environmental health.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Expanding standards in viromics: in silico evaluation of dsDNA viral genome identification, classification, and auxiliary metabolic gene curation

BACKGROUND

Viruses influence global patterns of microbial diversity and nutrient cycles. Though viral metagenomics (viromics), specifically targeting dsDNA viruses, has been critical for revealing viral roles across diverse ecosystems, its analyses differ in many ways from those used for microbes. To date, viromics benchmarking has covered read pre-processing, assembly, relative abundance, read mapping thresholds and diversity estimation, but other steps would benefit from benchmarking and standardization. Here we use in silico-generated datasets and an extensive literature survey to evaluate and highlight how dataset composition (i.e., viromes vs bulk metagenomes) and assembly fragmentation impact (i) viral contig identification tool, (ii) virus taxonomic classification, and (iii) identification and curation of auxiliary metabolic genes (AMGs).

RESULTS

The in silico benchmarking of five commonly used virus identification tools show that gene-content-based tools consistently performed well for long (≥3 kbp) contigs, while k-mer- and blast-based tools were uniquely able to detect viruses from short (≤3 kbp) contigs. Notably, however, the performance increase of k-mer- and blast-based tools for short contigs was obtained at the cost of increased false positives (sometimes up to ∼5% for virome and ∼75% bulk samples), particularly when eukaryotic or mobile genetic element sequences were included in the test datasets. For viral classification, variously sized genome fragments were assessed using gene-sharing network analytics to quantify drop-offs in taxonomic assignments, which revealed correct assignations ranging from ∼95% (whole genomes) down to ∼80% (3 kbp sized genome fragments). A similar trend was also observed for other viral classification tools such as VPF-class, ViPTree and VIRIDIC, suggesting that caution is warranted when classifying short genome fragments and not full genomes. Finally, we highlight how fragmented assemblies can lead to erroneous identification of AMGs and outline a best-practices workflow to curate candidate AMGs in viral genomes assembled from metagenomes.

CONCLUSION

Together, these benchmarking experiments and annotation guidelines should aid researchers seeking to best detect, classify, and characterize the myriad viruses 'hidden' in diverse sequence datasets.

Efficient dilution-to-extinction isolation of novel virus-host model systems for fastidious heterotrophic bacteria

Microbes and their associated viruses are key drivers of biogeochemical processes in marine and soil biomes. While viruses of phototrophic cyanobacteria are well-represented in model systems, challenges of isolating marine microbial heterotrophs and their viruses have hampered experimental approaches to quantify the importance of viruses in nutrient recycling. A resurgence in cultivation efforts has improved the availability of fastidious bacteria for hypothesis testing, but this has not been matched by similar efforts to cultivate their associated bacteriophages. Here, we describe a high-throughput method for isolating important virus-host systems for fastidious heterotrophic bacteria that couples advances in culturing of hosts with sequential enrichment and isolation of associated phages. Applied to six monthly samples from the Western English Channel, we first isolated one new member of the globally dominant bacterial SAR11 clade and three new members of the methylotrophic bacterial clade OM43. We used these as bait to isolate 117 new phages, including the first known siphophage-infecting SAR11, and the first isolated phage for OM43. Genomic analyses of 13 novel viruses revealed representatives of three new viral genera, and infection assays showed that the viruses infecting SAR11 have ecotype-specific host ranges. Similar to the abundant human-associated phage ɸCrAss001, infection dynamics within the majority of isolates suggested either prevalent lysogeny or chronic infection, despite a lack of associated genes, or host phenotypic bistability with lysis putatively maintained within a susceptible subpopulation. Broader representation of important virus-host systems in culture collections and genomic databases will improve both our understanding of virus-host interactions, and accuracy of computational approaches to evaluate ecological patterns from metagenomic data.

Ecogenomics and metabolic potential of the South Atlantic Ocean microbiome

The unique combination of depth, salinity, and water masses make the South Atlantic Ocean an ecosystem of special relevance within the global ocean. Yet, the microbiome of this ecosystem has received less attention than other regions of the global Ocean. This has hampered our understanding of the diversity and metabolic potential of the microorganisms that dwell in this habitat. To fill this knowledge gap, we analyzed a collection of 31 metagenomes from the Atlantic Ocean that spanned the epipelagic, mesopelagic and bathypelagic zones (surface to 4000 m). Read-centric and gene-centric analysis revealed the unique taxonomic and functional composition of metagenomes from each depth zone, which was driven by differences in physical and chemical parameters. In parallel, a total of 40 metagenome-assembled genomes were obtained, which recovered one third of the total community. Phylogenomic reconstruction revealed that many of these genomes are derived from poorly characterized taxa of Bacteria and Archaea. Genomes derived from heterotrophic bacteria of the aphotic zone displayed a large apparatus of genes suited for the utilization of recalcitrant organic compounds such as cellulose, chitin and alkanes. In addition, we found genomic evidence suggesting that mixotrophic bacteria from the bathypelagic zone could perform carbon fixation through the Calvin-Benson-Bassham cycle, fueled by sulfur oxidation. Finally, we found that the viral communities shifted throughout the water column regarding their targeted hosts and virus-to-microbe ratio, in response to shifts in the composition and functioning their microbial counterparts. Our findings shed light on the microbial and viral drivers of important biogeochemical processes that take place in the South Atlantic Ocean.

Hybrid assembly of an agricultural slurry virome reveals a diverse and stable community with the potential to alter the metabolism and virulence of veterinary pathogens

BACKGROUND

Viruses are the most abundant biological entities on Earth, known to be crucial components of microbial ecosystems. However, there is little information on the viral community within agricultural waste. There are currently ~ 2.7 million dairy cattle in the UK producing 7-8% of their own bodyweight in manure daily, and 28 million tonnes annually. To avoid pollution of UK freshwaters, manure must be stored and spread in accordance with guidelines set by DEFRA. Manures are used as fertiliser, and widely spread over crop fields, yet little is known about their microbial composition. We analysed the virome of agricultural slurry over a 5-month period using short and long-read sequencing.

RESULTS

Hybrid sequencing uncovered more high-quality viral genomes than long or short-reads alone; yielding 7682 vOTUs, 174 of which were complete viral genomes. The slurry virome was highly diverse and dominated by lytic bacteriophage, the majority of which represent novel genera (~ 98%). Despite constant influx and efflux of slurry, the composition and diversity of the slurry virome was extremely stable over time, with 55% of vOTUs detected in all samples over a 5-month period. Functional annotation revealed a diverse and abundant range of auxiliary metabolic genes and novel features present in the community, including the agriculturally relevant virulence factor VapE, which was widely distributed across different phage genera that were predicted to infect several hosts. Furthermore, we identified an abundance of phage-encoded diversity-generating retroelements, which were previously thought to be rare on lytic viral genomes. Additionally, we identified a group of crAssphages, including lineages that were previously thought only to be found in the human gut.

CONCLUSIONS

The cattle slurry virome is complex, diverse and dominated by novel genera, many of which are not recovered using long or short-reads alone. Phages were found to encode a wide range of AMGs that are not constrained to particular groups or predicted hosts, including virulence determinants and putative ARGs. The application of agricultural slurry to land may therefore be a driver of bacterial virulence and antimicrobial resistance in the environment. Video abstract.

The communities and functional profiles of virioplankton along a salinity gradient in a subtropical estuary

Viruses are the major drivers shaping microorganismal communities, and impact marine biogeochemical cycling. They are affected by various environmental parameters, such as salinity. Although the spatiotemporal distribution and dynamics of virioplankton have been extensively studied in saline environments, few detailed studies of community structure and function of viruses along salinity gradients have been conducted. Here, we used the 16S and 18S rRNA gene amplicon and metagenomic sequencing from a subtropical estuary (Pearl River Estuary, PRE; located in Shenzhen, Guangdong Province, China) to explore how viral community composition and function vary along a salinity gradient. Results showed that the detected viruses were mainly bacteriophages. The double-stranded DNA viruses were the most abundant (especially Siphoviridae, Myoviridae, Mimiviridae, Phycodnaviridae, and Podoviridae), followed by a small number of single-stranded DNA (Circoviridae) and RNA (Retroviridae) viruses. Viral biodiversity significantly declined and community structure varied greatly along the salinity gradient. The salinity, ammonium and dissolved oxygen were dominated factors influencing the community composition of viruses. Association network analysis showed that viruses had a negative effect on multiple host taxa (prokaryotic and eukaryotic species). Metagenomic data revealed that the main viral functional potential was involved in organic matter metabolism by carbohydrate-active enzymes (CAZymes). Deeper comparative functional analyses showed that viruses in the low-salinity environment had more carbohydrate-binding module and glycosidase hydrolases activities than those under high-salinity conditions. However, an opposite pattern was observed for carbohydrate esterases. These results suggest that virus-encoded CAZyme genes may alter the bacterial metabolism in estuaries. Overall, our results demonstrate that there is a spatial heterogeneity in the composition and function of virioplankton along a salinity gradient. This study enhances our understanding of viral distribution and their contribution to regulating carbon degradation throughout environments with varying salinities in subtropical estuaries.

Metaviromics coupled with phage-host identification to open the viral 'black box'

Viruses are found in almost all biomes on Earth, with bacteriophages (phages) accounting for the majority of viral particles in most ecosystems. Phages have been isolated from natural environments using the plaque assay and liquid medium-based dilution culturing. However, phage cultivation is restricted by the current limitations in the number of culturable bacterial strains. Unlike prokaryotes, which possess universally conserved 16S rRNA genes, phages lack universal marker genes for viral taxonomy, thus restricting cultureindependent analyses of viral diversity. To circumvent these limitations, shotgun viral metagenome sequencing (i.e., metaviromics) has been developed to enable the extensive sequencing of a variety of viral particles present in the environment and is now widely used. Using metaviromics, numerous studies on viral communities have been conducted in oceans, lakes, rivers, and soils, resulting in many novel phage sequences. Furthermore, auxiliary metabolic genes such as ammonic monooxygenase C and β-lactamase have been discovered in viral contigs assembled from viral metagenomes. Current attempts to identify putative bacterial hosts of viral metagenome sequences based on sequence homology have been limited due to viral sequence variations. Therefore, culture-independent approaches have been developed to predict bacterial hosts using single-cell genomics and fluorescentlabeling. This review focuses on recent viral metagenome studies conducted in natural environments, especially in aquatic ecosystems, and their contributions to phage ecology. Here, we concluded that although metaviromics is a key tool for the study of viral ecology, this approach must be supplemented with phage-host identification, which in turn requires the cultivation of phage-bacteria systems.

Eukaryotic virus composition can predict the efficiency of carbon export in the global ocean

The biological carbon pump, in which carbon fixed by photosynthesis is exported to the deep ocean through sinking, is a major process in Earth's carbon cycle. The proportion of primary production that is exported is termed the carbon export efficiency (CEE). Based on in-lab or regional scale observations, viruses were previously suggested to affect the CEE (i.e., viral “shunt” and “shuttle”). In this study, we tested associations between viral community composition and CEE measured at a global scale. A regression model based on relative abundance of viral marker genes explained 67% of the variation in CEE. Viruses with high importance in the model were predicted to infect ecologically important hosts. These results are consistent with the view that the viral shunt and shuttle functions at a large scale and further imply that viruses likely act in this process in a way dependent on their hosts and ecosystem dynamics.

•

Eukaryotic virus community composition is shown to predict carbon export efficiency

•

Tens of viruses are highly important in the prediction of the efficiency

•

These viruses are inferred to infect ecologically important hosts

New viral biogeochemical roles revealed through metagenomic analysis of Lake Baikal

BACKGROUND

Lake Baikal is the largest body of liquid freshwater on Earth. Previous studies have described the microbial composition of this habitat, but the viral communities from this ecosystem have not been characterized in detail.

RESULTS

Here, we describe the viral diversity of this habitat across depth and seasonal gradients. We discovered 19,475 bona fide viral sequences, which are derived from viruses predicted to infect abundant and ecologically important taxa that reside in Lake Baikal, such as Nitrospirota, Methylophilaceae, and Crenarchaeota. Diversity analysis revealed significant changes in viral community composition between epipelagic and bathypelagic zones. Analysis of the gene content of individual viral populations allowed us to describe one of the first bacteriophages that infect Nitrospirota, and their extensive repertoire of auxiliary metabolic genes that might enhance carbon fixation through the reductive TCA cycle. We also described bacteriophages of methylotrophic bacteria with the potential to enhance methanol oxidation and the S-adenosyl-L-methionine cycle.

CONCLUSIONS

These findings unraveled new ways by which viruses influence the carbon cycle in freshwater ecosystems, namely, by using auxiliary metabolic genes that act upon metabolisms of dark carbon fixation and methylotrophy. Therefore, our results shed light on the processes through which viruses can impact biogeochemical cycles of major ecological relevance. Video Abstract.

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 elements and their potential influence on microbial processes along the permanently stratified Cariaco Basin redoxcline

Little is known about viruses in oxygen-deficient water columns (ODWCs). In surface ocean waters, viruses are known to act as gene vectors among susceptible hosts. Some of these genes may have metabolic functions and are thus termed auxiliary metabolic genes (AMGs). AMGs introduced to new hosts by viruses can enhance viral replication and/or potentially affect biogeochemical cycles by modulating key microbial pathways. Here we identify 748 viral populations that cluster into 94 genera along a vertical geochemical gradient in the Cariaco Basin, a permanently stratified and euxinic ocean basin. The viral communities in this ODWC appear to be relatively novel as 80 of these viral genera contained no reference viral sequences, likely due to the isolation and unique features of this system. We identify viral elements that encode AMGs implicated in distinctive processes, such as sulfur cycling, acetate fermentation, signal transduction, [Fe–S] formation, and N-glycosylation. These AMG-encoding viruses include two putative Mu-like viruses, and viral-like regions that may constitute degraded prophages that have been modified by transposable elements. Our results provide an insight into the ecological and biogeochemical impact of viruses oxygen-depleted and euxinic habitats.

Genomic Characteristics and Potential Metabolic Adaptations of Hadal Trench Roseobacter and Alteromonas Bacteria Based on Single-Cell Genomics Analyses

Heterotrophic bacteria such as those from the Roseobacter group and genus Alteromonas dominate the hadal zones of oceans; however, we know little about the genomic characteristics and potential metabolic adaptations of hadal trench-dwelling bacteria. Here, we report multiple single amplified genomes (SAGs) belonging to Roseobacter and Alteromonas, recovered from the hadal zone of the Mariana Trench. While phylogenetic analyses show that these hadal SAGs cluster with their surface relatives, an analysis of genomic recruitment indicates that they have higher relative abundances in the hadal zone of the Mariana Trench. Comparative genomic analyses between the hadal SAGs and reference genomes of closely related shallow-water relatives indicate that genes involved in the mobilome (prophages and transposons) are overrepresented among the unique genes of the hadal Roseobacter and Alteromonas SAGs; the functional proteins encoded by this category of genes also shows higher amino acid sequence variation than those encoded by other gene sets within the Roseobacter SAGs. We also found that genes involved in cell wall/membrane/envelope biogenesis, transcriptional regulation, and metal transport may be important for the adaptation of hadal Roseobacter and Alteromonas lineages. These results imply that the modification of cell surface-related proteins and transporters is the major direction of genomic evolution in Roseobacter and Alteromonas bacteria adapting to the hadal environment, and that prophages and transposons may be the key factors driving this process.

Assessing Viral Abundance and Community Composition in Four Contrasting Regions of the Southern Ocean

We explored how changes of viral abundance and community composition among four contrasting regions in the Southern Ocean relied on physicochemical and microbiological traits. During January-February 2015, we visited areas north and south of the South Orkney Islands (NSO and SSO) characterized by low temperature and salinity and high inorganic nutrient concentration, north of South Georgia Island (NSG) and west of Anvers Island (WA), which have relatively higher temperatures and lower inorganic nutrient concentrations. Surface viral abundance (VA) was highest in NSG (21.50 ± 10.70 × 10⁶ viruses mL^-1) and lowest in SSO (2.96 ± 1.48 × 10⁶ viruses mL^-1). VA was positively correlated with temperature, prokaryote abundance and prokaryotic heterotrophic production, chlorophyll a, diatoms, haptophytes, fluorescent organic matter, and isoprene concentration, and was negatively correlated with inorganic nutrients (NO_3-, SiO₄^2-, PO₄^3-), and dimethyl sulfide (DMS) concentrations. Viral communities determined by randomly amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR) were grouped according to the sampling location, being more similar within them than among regions. The first two axes of a canonical correspondence analysis, including physicochemical (temperature, salinity, inorganic nutrients-NO_3-, SiO₄^2-, and dimethyl sulfoniopropionate -DMSP- and isoprene concentrations) and microbiological (chlorophyll a, haptophytes and diatom, and prokaryote abundance and prokaryotic heterotrophic production) factors accounted for 62.9% of the variance. The first axis, temperature-related, accounted for 33.8%; the second one, salinity-related, accounted for 29.1%. Thus, different environmental situations likely select different hosts for viruses, leading to distinct viral communities.

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.