Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses

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.

Genomic analysis and characterization of phages infecting the marine Roseobacter CHAB-I-5 lineage reveal a globally distributed and abundant phage genus

Marine phages play an important role in marine biogeochemical cycles by regulating the death, physiological metabolism, and evolutionary trajectory of bacteria. The Roseobacter group is an abundant and important heterotrophic bacterial group in the ocean, and plays an important role in carbon, nitrogen, sulfur and phosphorus cycling. The CHAB-I-5 lineage is one of the most dominant Roseobacter lineages, but remains largely uncultured. Phages infecting CHAB-I-5 bacteria have not yet been investigated due to the lack of culturable CHAB-I-5 strains. In this study, we isolated and sequenced two new phages (CRP-901 and CRP-902) infecting the CHAB-I-5 strain FZCC0083. We applied metagenomic data mining, comparative genomics, phylogenetic analysis, and metagenomic read-mapping to investigate the diversity, evolution, taxonomy, and biogeography of the phage group represented by the two phages. The two phages are highly similar, with an average nucleotide identity of 89.17%, and sharing 77% of their open reading frames. We identified several genes involved in DNA replication and metabolism, virion structure, DNA packing, and host lysis from their genomes. Metagenomic mining identified 24 metagenomic viral genomes closely related to CRP-901 and CRP-902. Genomic comparison and phylogenetic analysis demonstrated that these phages are distinct from other known viruses, representing a novel genus-level phage group (CRP-901-type). The CRP-901-type phages do not contain DNA primase and DNA polymerase genes, but possess a novel bifunctional DNA primase-polymerase gene with both primase and polymerase activities. Read-mapping analysis showed that the CRP-901-type phages are widespread across the world's oceans and are most abundant in estuarine and polar waters. Their abundance is generally higher than other known roseophages and even higher than most pelagiphages in the polar region. In summary, this study has greatly expanded our understanding of the genetic diversity, evolution, and distribution of roseophages. Our analysis suggests that the CRP-901-type phage is an important and novel marine phage group that plays important roles in the physiology and ecology of roseobacters.

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.

Significant Differences in Planktonic Virus Communities Between "Cellular Fraction" (0.22 ~ 3.0 µm) and "Viral Fraction" (< 0.22 μm) in the Ocean

Compared to free-living viruses (< 0.22 m) in the ocean, planktonic viruses in the "cellular fraction" (0.22 ~ 3.0 μm) are now far less well understood, and the differences between them remain largely unexplored. Here, we revealed that even in the same seawater samples, the "cellular fraction" comprised significantly distinct virus communities from the free virioplankton, with only 13.87% overlap in viral contigs at the species level. Compared to the viral genomes deposited in NCBI RefSeq database, 99% of the assembled viral genomes in the "cellular fraction" represented novel genera. Notably, the assembled (near-) complete viral genomes within the "cellular fraction" were significantly larger than that in the "viral fraction," and the "cellular fraction" contained three times more species of giant viruses or jumbo phages with genomes > 200 kb than the "viral fraction." The longest complete genomes of jumbo phage (~ 252 kb) and giant virus (~ 716 kb) were both detected only in the "cellular fraction." Moreover, a relatively higher proportion of proviruses were predicted within the "cellular fraction" than "viral fraction." Besides the substantial divergence in viral community structure, the different fractions also contained their unique viral auxiliary metabolic genes; e.g., those potentially participating in inorganic carbon fixation in deep sea were detected only in the "cellular-fraction" viromes. In addition, there was a considerable divergence in the community structure of both "cellular fraction" and "viral fraction" viromes between the surface and deep-sea habitats, suggesting that they might have similar environmental adaptation properties. The findings deepen our understanding of the complexity of viral community structure and function in the ocean.

Identification of Novel Viruses and Their Microbial Hosts from Soils with Long-Term Nitrogen Fertilization and Cover Cropping Management

Soils are the largest organic carbon reservoir and are key to global biogeochemical cycling, and microbes are the major drivers of carbon and nitrogen transformations in the soil systems. Thus, virus infection-induced microbial mortality could impact soil microbial structure and functions. In this study, we recovered 260 viral operational taxonomic units (vOTUs) in samples collected from soil taken from four nitrogen fertilization (N-fertilization) and cover-cropping practices at an experimental site under continuous cotton production evaluating conservation agricultural management systems for more than 40 years. Only ~6% of the vOTUs identified were clustered with known viruses in the RefSeq database using a gene-sharing network. We found that 14% of 260 vOTUs could be linked to microbial hosts that cover key carbon and nitrogen cycling taxa, including Acidobacteriota, Proteobacteria, Verrucomicrobiota, Firmicutes, and ammonia-oxidizing archaea, i.e., Nitrososphaeria (phylum Thermoproteota). Viral diversity, community structure, and the positive correlation between abundance of a virus and its host indicate that viruses and microbes are more sensitive to N-fertilization than cover-cropping treatment. Viruses may influence key carbon and nitrogen cycling through control of microbial function and host populations (e.g., Chthoniobacterales and Nitrososphaerales). These findings provide an initial view of soil viral ecology and how it is influenced by long-term conservation agricultural management. IMPORTANCE Bacterial viruses are extremely small and abundant particles that can control the microbial abundance and community composition through infection, which gradually showed their vital roles in the ecological process to influence the nutrient flow. Compared to the substrate control, less is known about the influence of soil viruses on microbial community function, and even less is known about microbial and viral diversity in the soil system. To obtain a more complete knowledge of microbial function dynamics, the interaction between microbes and viruses cannot be ignored. To fully understand this process, it is fundamental to get insight into the correlation between the diversity of viral communities and bacteria which could induce these changes.

Viral community in landfill leachate: Occurrence, bacterial hosts, mediation antibiotic resistance gene dissemination, and function in municipal solid waste decomposition

A municipal solid waste (MSW) landfill is a significant source of antibiotic resistance, pathogens and viruses and also a habitat for microbial consortia that perform MSW decomposition. Viruses are of great significance in ecological interactions such as MSW decomposition and antibiotic resistance gene (ARG) transmission. In this study, the viral community structure and the associated driver, the linkage of viruses and their bacterial hosts, the virus-associated ARG dissemination and virtual community function on MSW decomposition were investigated in landfill leachate from seven cities, China. The seven cities include four megacities, two large-scale cities and one small-scale city, representing the leachate characters of China. The results showed that the leachates were dominated by the phage families Siphoviridae, Myoviridae and Podoviridae (91.7 ± 3.6) %. Their putative hosts were the important MSW decomposers Lactobacillus, Pseudomonas, Clostridium, Proteiniphilum, and Bacteroides. The structure of the viral community was significantly affected by pH (P = 0.007, analyzed by RDA) and the bacterial community (R = 0.83, P < 0.001, analyzed by Mantel test). The relative abundance of ARGs showed a strong correlation (R > 0.8, P < 0.01) with viral family, suggesting that viruses play an important role in ARGs dissemination. Phage regulate bacterial population abundance through top-down effects, thus participating in MSW decomposition. These results demonstrate that viral community are involve in ARGs transmission and dissemination and mediate MSW decomposition in landfill.

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.

Discovering marine biodiversity in the 21st century

We review the current knowledge of the biodiversity of the ocean as well as the levels of decline and threat for species and habitats. The lack of understanding of the distribution of life in the ocean is identified as a significant barrier to restoring its biodiversity and health. We explore why the science of taxonomy has failed to deliver knowledge of what species are present in the ocean, how they are distributed and how they are responding to global and regional to local anthropogenic pressures. This failure prevents nations from meeting their international commitments to conserve marine biodiversity with the results that investment in taxonomy has declined in many countries. We explore a range of new technologies and approaches for discovery of marine species and their detection and monitoring. These include: imaging methods, molecular approaches, active and passive acoustics, the use of interconnected databases and citizen science. Whilst no one method is suitable for discovering or detecting all groups of organisms many are complementary and have been combined to give a more complete picture of biodiversity in marine ecosystems. We conclude that integrated approaches represent the best way forwards for accelerating species discovery, description and biodiversity assessment. Examples of integrated taxonomic approaches are identified from terrestrial ecosystems. Such integrated taxonomic approaches require the adoption of cybertaxonomy approaches and will be boosted by new autonomous sampling platforms and development of machine-speed exchange of digital information between databases.

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.

Characterization of the human gut virome in metabolic and autoimmune diseases

The intestinal microbiome is dominated by bacteria and plays a pivotal role in the occurrence and development of disease, including several metabolic and autoimmune disorders. While intestinal viral communities, primarily made up of bacteriophages, are also thought to play a role in disease pathogenesis in the gastrointestinal tract, they have received much less attention than intestinal bacteria. Thus, there is limited information about the relationship between bacteriophages and disease. This review explores a potential role for the intestinal viral microbiome in various metabolic and autoimmune diseases.

Viruses in astrobiology

Viruses are the most abundant biological entities on Earth, and yet, they have not received enough consideration in astrobiology. Viruses are also extraordinarily diverse, which is evident in the types of relationships they establish with their host, their strategies to store and replicate their genetic information and the enormous diversity of genes they contain. A viral population, especially if it corresponds to a virus with an RNA genome, can contain an array of sequence variants that greatly exceeds what is present in most cell populations. The fact that viruses always need cellular resources to multiply means that they establish very close interactions with cells. Although in the short term these relationships may appear to be negative for life, it is evident that they can be beneficial in the long term. Viruses are one of the most powerful selective pressures that exist, accelerating the evolution of defense mechanisms in the cellular world. They can also exchange genetic material with the host during the infection process, providing organisms with capacities that favor the colonization of new ecological niches or confer an advantage over competitors, just to cite a few examples. In addition, viruses have a relevant participation in the biogeochemical cycles of our planet, contributing to the recycling of the matter necessary for the maintenance of life. Therefore, although viruses have traditionally been excluded from the tree of life, the structure of this tree is largely the result of the interactions that have been established throughout the intertwined history of the cellular and the viral worlds. We do not know how other possible biospheres outside our planet could be, but it is clear that viruses play an essential role in the terrestrial one. Therefore, they must be taken into account both to improve our understanding of life that we know, and to understand other possible lives that might exist in the cosmos.

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.

Crop management shapes the diversity and activity of DNA and RNA viruses in the rhizosphere

BACKGROUND

The rhizosphere is a hotspot for microbial activity and contributes to ecosystem services including plant health and biogeochemical cycling. The activity of microbial viruses, and their influence on plant-microbe interactions in the rhizosphere, remains undetermined. Given the impact of viruses on the ecology and evolution of their host communities, determining how soil viruses influence microbiome dynamics is crucial to build a holistic understanding of rhizosphere functions.

RESULTS

Here, we aimed to investigate the influence of crop management on the composition and activity of bulk soil, rhizosphere soil, and root viral communities. We combined viromics, metagenomics, and metatranscriptomics on soil samples collected from a 3-year crop rotation field trial of oilseed rape (Brassica napus L.). By recovering 1059 dsDNA viral populations and 16,541 ssRNA bacteriophage populations, we expanded the number of underexplored Leviviricetes genomes by > 5 times. Through detection of viral activity in metatranscriptomes, we uncovered evidence of "Kill-the-Winner" dynamics, implicating soil bacteriophages in driving bacterial community succession. Moreover, we found the activity of viruses increased with proximity to crop roots, and identified that soil viruses may influence plant-microbe interactions through the reprogramming of bacterial host metabolism. We have provided the first evidence of crop rotation-driven impacts on soil microbial communities extending to viruses. To this aim, we present the novel principal of "viral priming," which describes how the consecutive growth of the same crop species primes viral activity in the rhizosphere through local adaptation.

CONCLUSIONS

Overall, we reveal unprecedented spatial and temporal diversity in viral community composition and activity across root, rhizosphere soil, and bulk soil compartments. Our work demonstrates that the roles of soil viruses need greater consideration to exploit the rhizosphere microbiome for food security, food safety, and environmental sustainability. Video Abstract.

Expansion of the global RNA virome reveals diverse clades of bacteriophages

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

Viruses direct carbon cycling in lake sediments under global change

Global change is altering the vast amount of carbon cycled by microbes between land and freshwater, but how viruses mediate this process is poorly understood. Here, we show that viruses direct carbon cycling in lake sediments, and these impacts intensify with future changes in water clarity and terrestrial organic matter (tOM) inputs. Using experimental tOM gradients within sediments of a clear and a dark boreal lake, we identified 156 viral operational taxonomic units (vOTUs), of which 21% strongly increased with abundances of key bacteria and archaea, identified via metagenome-assembled genomes (MAGs). MAGs included the most abundant prokaryotes, which were themselves associated with dissolved organic matter (DOM) composition and greenhouse gas (GHG) concentrations. Increased abundances of virus-like particles were separately associated with reduced bacterial metabolism and with shifts in DOM toward amino sugars, likely released by cell lysis rather than higher molecular mass compounds accumulating from reduced tOM degradation. An additional 9.6% of vOTUs harbored auxiliary metabolic genes associated with DOM and GHGs. Taken together, these different effects on host dynamics and metabolism can explain why abundances of vOTUs rather than MAGs were better overall predictors of carbon cycling. Future increases in tOM quantity, but not quality, will change viral composition and function with consequences for DOM pools. Given their importance, viruses must now be explicitly considered in efforts to understand and predict the freshwater carbon cycle and its future under global environmental change.

Ecogenomics reveals viral communities across the Challenger Deep oceanic trench

Despite the environmental challenges and nutrient scarcity, the geographically isolated Challenger Deep in Mariana trench, is considered a dynamic hotspot of microbial activity. Hadal viruses are the least explored microorganisms in Challenger Deep, while their taxonomic and functional diversity and ecological impact on deep-sea biogeochemistry are poorly described. Here, we collect 13 sediment cores from slope and bottom-axis sites across the Challenger Deep (down to ~11 kilometers depth), and identify 1,628 previously undescribed viral operational taxonomic units at species level. Community-wide analyses reveals 1,299 viral genera and distinct viral diversity across the trench, which is significantly higher at the bottom-axis vs. slope sites of the trench. 77% of these viral genera have not been previously identified in soils, deep-sea sediments and other oceanic settings. Key prokaryotes involved in hadal carbon and nitrogen cycling are predicted to be potential hosts infected by these viruses. The detected putative auxiliary metabolic genes suggest that viruses at Challenger Deep could modulate the carbohydrate and sulfur metabolisms of their potential hosts, and stabilize host’s cell membranes under extreme hydrostatic pressures. Our results shed light on hadal viral metabolic capabilities, contribute to understanding deep sea ecology and on functional adaptions of hadal viruses for future research.

Analysis of 13 sediment cores from the Challenger Deep of Marian Trench (down to 11 kilometers depth) identified distinct operational taxonomic units and relevant auxiliary metabolic genes, providing further insight into deep-sea viral metabolic capabilities and ecology.

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.

Nutrient levels and prokaryotes affect viral communities in plateau lakes

Viruses are the most abundant organisms in aquatic environments. Recent advances of viral metagenomic have greatly expanded our understanding of aquatic viral communities. However, little is known about the difference of viral communities and driving factors in freshwater lake. This study seeks to understand the spatio-temporal variation, differences, and driving factors of viral communities in two plateau lakes (Dianchi and Fuxian Lakes) with significant nutritional differences. The viral communities exhibited apparent seasonal variation in Dianchi Lake, while seasonal influences on the viral communities were greater than location-based influences. Two-thirds of all detected viral taxa were shared in two lakes, but there was variation in the composition of viral communities. Correlations between prokaryotic communities, environmental factors and viral communities were analyzed. The nutrients, chlorophyll a were primarily environmental parameters affecting viral communities, and the prokaryotic community was significantly correlated with the viral community. In addition, several viruses infecting humans were identified in two lakes, with the most abundant being Herpesviridae and Poxviridae. Overall, these findings provide information on the dynamics, composition, and differences of viral and prokaryotic communities in plateau lakes with different nutrient levels. These results suggest that nutritional levels and prokaryotic communities could play an important role in shaping viral communities in freshwater lakes.

Novel lineages of single-stranded DNA phages that coevolved with the symbiotic bacteria Rhizobium

This study describes novel single-stranded DNA phages isolated from common bean agriculture soils by infection of the nitrogen-fixing symbiotic bacteria Rhizobium etli and R. phaseoli. A total of 29 phages analyzed have 4.3–6 kb genomes in size and GC 59–60%. They belong to different clades unrelated to other Microviridae subfamilies. Three-dimensional models of the major capsid protein (MCP) showed a conserved β-barrel structural “jelly-roll” fold. A variable-length loop in the MCPs distinguished three Rhizobium microvirus groups. Microviridae subfamilies were consistent with viral clusters determined by the protein-sharing network. All viral clusters, except for Bullavirinae, included mostly microviruses identified in metagenomes from distinct ecosystems. Two Rhizobium microvirus clusters, chaparroviruses, and chicoviruses, were included within large viral unknown clusters with microvirus genomes identified in diverse metagenomes. A third Rhizobium microvirus cluster belonged to the subfamily Amoyvirinae. Phylogenetic analysis of the MCP confirms the divergence of the Rhizobium microviruses into separate clades. The phylogeny of the bacterial hosts matches the microvirus MCP phylogeny, suggesting a coevolutionary history between the phages and their bacterial host. This study provided essential biological information on cultivated microvirus for understanding the evolution and ecological diversification of the Microviridae family in diverse microbial ecosystems.

Periphytic biofilms-mediated microbial interactions and their impact on the nitrogen cycle in rice paddies

Rice paddies are unique waterlogged wetlands artificially constructed for agricultural production. Periphytic biofilms (PBs) at the soil-water interface play an important role in rice paddies characterized by high nutrient input but low utilization efficiency. PBs are composed of microbial aggregates, including a wide variety of microorganisms (algae, bacteria, fungi, protozoa, and metazoa), extracellular polymeric substances and minerals (iron, aluminum, and calcium), which form an integrated food web and energy flux within a relatively stable micro-ecosystem. PBs are crucial to regulate and streamline the nitrogen cycle by neutralizing nitrogen losses and improving rice production since PBs can serve as both a sink by capturing surplus nitrogen and a source by slowly re-releasing this nitrogen for reutilization. Here the ecological advantages of PBs in regulating the nitrogen cycle in rice paddies are illustrated. We summarize the key functional importance of PBs, including the intricate and delicate community structure, microbial interactions among individual phylotypes, a wide diversity of self-produced organics, the active adaptation of PBs to constantly changing environments, and the intricate mechanisms by which PBs regulate the nitrogen cycle. We also identify the future challenges of microbial interspecific cooperation in PBs and their quantitative contributions to agricultural sustainability, optimizing nitrogen utilization and crop yields in rice paddies.

Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems

Biogeographical studies have traditionally focused on readily visible organisms, but recent technological advances are enabling analyses of the large-scale distribution of microscopic organisms, whose biogeographical patterns have long been debated. Here we assessed the global structure of plankton geography and its relation to the biological, chemical, and physical context of the ocean (the ‘seascape’) by analyzing metagenomes of plankton communities sampled across oceans during the Tara Oceans expedition, in light of environmental data and ocean current transport. Using a consistent approach across organismal sizes that provides unprecedented resolution to measure changes in genomic composition between communities, we report a pan-ocean, size-dependent plankton biogeography overlying regional heterogeneity. We found robust evidence for a basin-scale impact of transport by ocean currents on plankton biogeography, and on a characteristic timescale of community dynamics going beyond simple seasonality or life history transitions of plankton.

Oceans are brimming with life invisible to our eyes, a myriad of species of bacteria, viruses and other microscopic organisms essential for the health of the planet. These ‘marine plankton’ are unable to swim against currents and should therefore be constantly on the move, yet previous studies have suggested that distinct species of plankton may in fact inhabit different oceanic regions. However, proving this theory has been challenging; collecting plankton is logistically difficult, and it is often impossible to distinguish between species simply by examining them under a microscope. However, within the last decade, a research schooner called Tara has travelled the globe to gather thousands of plankton samples. At the same time, advances in genomics have made it possible to identify species based only on fragments of their DNA sequence.

To understand the hidden geography of plankton communities in Earth’s oceans, Richter et al. pored over DNA from the Tara Oceans expedition. This revealed that, despite being unable to resist the flow of water, various planktonic species which live close to the surface manage to occupy distinct, stable provinces shaped by currents. Different sizes of plankton are distributed in different sized provinces, with the smallest organisms tending to inhabit the smallest areas. Comparing DNA similarities and speeds of currents at the ocean surface revealed how these might stretch and mix plankton communities.

Plankton play a critical role in the health of the ocean and the chemical cycles of planet Earth. These results could allow deeper investigation by marine modellers, ecologists, and evolutionary biologists. Meanwhile, work is already underway to investigate how climate change might impact this hidden geography.

Virioplankton assemblages from challenger deep, the deepest place in the oceans

Hadal ocean biosphere, that is, the deepest part of the world’s oceans, harbors a unique microbial community, suggesting a potential uncovered co-occurring virioplankton assemblage. Herein, we reveal the unique virioplankton assemblages of the Challenger Deep, comprising 95,813 non-redundant viral contigs from the surface to the hadal zone. Almost all of the dominant viral contigs in the hadal zone were unclassified, potentially related to Alteromonadales and Oceanospirillales. 2,586 viral auxiliary metabolic genes from 132 different KEGG orthologous groups were mainly related to the carbon, nitrogen, sulfur, and arsenic metabolism. Lysogenic viral production and integrase genes were augmented in the hadal zone, suggesting the prevalence of viral lysogenic life strategy. Abundant rve genes in the hadal zone, which function as transposase in the caudoviruses, further suggest the prevalence of viral-mediated horizontal gene transfer. This study provides fundamental insights into the virioplankton assemblages of the hadal zone, reinforcing the necessity of incorporating virioplankton into the hadal biogeochemical cycles.

•

The unique virioplankton assemblages of the Challenger Deep were revealed

•

Virus encoded auxiliary metabolic genes relating to the biogeochemical cycling

•

Viruses in deep and hadal zone tend to be lysogenic, and potentially mediate the horizontal gene transfer

Capturing the environment of the Clostridioides difficile infection cycle

Clostridioides difficile (formerly Clostridium difficile) infection is a substantial health and economic burden worldwide. Great strides have been made over the past several years in characterizing the physiology of C. difficile infection, particularly regarding how gut microorganisms and their host work together to provide colonization resistance. As mammalian hosts and their indigenous gut microbiota have co-evolved, they have formed a complex yet stable relationship that prevents invading microorganisms from establishing themselves. In this Review, we discuss the latest advances in our understanding of C. difficile physiology that have contributed to its success as a pathogen, including its versatile survival factors and ability to adapt to unique niches. Using discoveries regarding microorganism-host and microorganism-microorganism interactions that constitute colonization resistance, we place C. difficile within the fiercely competitive gut environment. A comprehensive understanding of these relationships is required to continue the development of precision medicine-based treatments for C. difficile infection.

New Microviridae isolated from Sulfitobacter reveals two cosmopolitan subfamilies of single-stranded DNA phages infecting marine and terrestrial Alphaproteobacteria

The Microviridae family represents one of the major clades of single-stranded DNA (ssDNA) phages. Their cultivated members are lytic and infect Proteobacteria, Bacteroidetes, and Chlamydiae. Prophages have been predicted in the genomes from Bacteroidales, Hyphomicrobiales, and Enterobacteriaceae and cluster within the 'Alpavirinae', 'Amoyvirinae', and Gokushovirinae. We have isolated 'Ascunsovirus oldenburgi' ICBM5, a novel phage distantly related to known Microviridae. It infects Sulfitobacter dubius SH24-1b and uses both a lytic and a carrier-state life strategy. Using ICBM5 proteins as a query, we uncovered in publicly available resources sixty-five new Microviridae prophages and episomes in bacterial genomes and retrieved forty-seven environmental viral genomes (EVGs) from various viromes. Genome clustering based on protein content and phylogenetic analysis showed that ICBM5, together with Rhizobium phages, new prophages, episomes, and EVGs cluster within two new phylogenetic clades, here tentatively assigned the rank of subfamily and named 'Tainavirinae' and 'Occultatumvirinae'. They both infect Rhodobacterales. Occultatumviruses also infect Hyphomicrobiales, including nitrogen-fixing endosymbionts from cosmopolitan legumes. A biogeographical assessment showed that tainaviruses and occultatumviruses are spread worldwide, in terrestrial and marine environments. The new phage isolated here sheds light onto new and diverse branches of the Microviridae tree, suggesting that much of the ssDNA phage diversity remains in the dark.

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.

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.

A snapshot of the global drinking water virome: Diversity and metabolic potential vary with residual disinfectant use

Viruses are important drivers of microbial community ecology and evolution, influencing microbial mortality, metabolism, and horizontal gene transfer. However, the effects of viruses remain largely unknown in many environments, including in drinking water systems. Drinking water metagenomic studies have offered a whole community perspective of bacterial impacts on water quality, but have not yet considered the influences of viruses. In this study, we address this gap by mining viral DNA sequences from publicly available drinking water metagenomes from distribution systems in six countries around the world. These datasets provide a snapshot of the taxonomic diversity and metabolic potential of the global drinking water virome; and provide an opportunity to investigate the effects of geography, climate, and drinking water treatment practices on viral diversity. Both environmental conditions and differences in sample processing were found to influence the viral composition. Using free chlorine as the residual disinfectant was associated with clear differences in viral taxonomic diversity and metabolic potential, with significantly fewer viral populations and less even viral community structures than observed in distribution systems without residual disinfectant. Additionally, drinking water viruses carry antibiotic resistance genes (ARGs), as well as genes to survive oxidative stress and nitrogen limitation. Through this study, we have demonstrated that viral communities are diverse across drinking water systems and vary with the use of residual disinfectant. Our findings offer directions for future research to develop a more robust understanding of how virus-bacteria interactions in drinking water distribution systems affect water quality.

Potential metabolic and genetic interaction among viruses, methanogen and methanotrophic archaea, and their syntrophic partners

The metabolism of methane in anoxic ecosystems is mainly mediated by methanogens and methane-oxidizing archaea (MMA), key players in global carbon cycling. Viruses are vital in regulating their host fate and ecological function. However, our knowledge about the distribution and diversity of MMA viruses and their interactions with hosts is rather limited. Here, by searching metagenomes containing mcrA (the gene coding for the α-subunit of methyl-coenzyme M reductase) from a wide variety of environments, 140 viral operational taxonomic units (vOTUs) that potentially infect methanogens or methane-oxidizing archaea were retrieved. Four MMA vOTUs (three infecting the order Methanobacteriales and one infecting the order Methanococcales) were predicted to cross-domain infect sulfate-reducing bacteria. By facilitating assimilatory sulfur reduction, MMA viruses may increase the fitness of their hosts in sulfate-depleted anoxic ecosystems and benefit from synthesis of the sulfur-containing amino acid cysteine. Moreover, cell-cell aggregation promoted by MMA viruses may be beneficial for both the viruses and their hosts by improving infectivity and environmental stress resistance, respectively. Our results suggest a potential role of viruses in the ecological and environmental adaptation of methanogens and methane-oxidizing archaea.

Computational Tools for the Analysis of Uncultivated Phage Genomes

Over a century of bacteriophage research has uncovered a plethora of fundamental aspects of their biology, ecology, and evolution. Furthermore, the introduction of community-level studies through metagenomics has revealed unprecedented insights on the impact that phages have on a range of ecological and physiological processes. It was not until the introduction of viral metagenomics that we began to grasp the astonishing breadth of genetic diversity encompassed by phage genomes. Novel phage genomes have been reported from a diverse range of biomes at an increasing rate, which has prompted the development of computational tools that support the multilevel characterization of these novel phages based solely on their genome sequences. The impact of these technologies has been so large that, together with MAGs (Metagenomic Assembled Genomes), we now have UViGs (Uncultivated Viral Genomes), which are now officially recognized by the International Committee for the Taxonomy of Viruses (ICTV), and new taxonomic groups can now be created based exclusively on genomic sequence information. Even though the available tools have immensely contributed to our knowledge of phage diversity and ecology, the ongoing surge in software programs makes it challenging to keep up with them and the purpose each one is designed for. Therefore, in this review, we describe a comprehensive set of currently available computational tools designed for the characterization of phage genome sequences, focusing on five specific analyses: (i) assembly and identification of phage and prophage sequences, (ii) phage genome annotation, (iii) phage taxonomic classification, (iv) phage-host interaction analysis, and (v) phage microdiversity.

Investigation of the Phageome and Prophages in French Cider, a Fermented Beverage

Phageomes are known to play a key role in the functioning of their associated microbial communities. The phageomes of fermented foods have not been studied thoroughly in fermented foods yet, and even less in fermented beverages. Two approaches were employed to investigate the presence of phages in cider, a fermented beverage made from apple, during a fermentation process of two cider tanks, one from an industrial producer and one from a hand-crafted producer. The phageome (free lytic phages) was explored in cider samples with several methodological developments for total phage DNA extraction, along with single phage isolation. Concentration methods, such as tangential flow filtration, flocculation and classical phage concentration methods, were employed and tested to extract free phage particles from cider. This part of the work revealed a very low occurrence of free lytic phage particles in cider. In parallel, a prophage investigation during the fermentation process was also performed using a metagenomic approach on the total bacterial genomic DNA. Prophages in bacterial metagenomes in the two cider tanks seemed also to occur in low abundance, as a total of 1174 putative prophages were identified in the two tanks overtime, and only two complete prophages were revealed. Prophage occurrence was greater at the industrial producer than at the hand-crafted producer, and different dynamics of prophage trends were also observed during fermentation. This is the first report dealing with the investigation of the phageome and of prophages throughout a fermentation process of a fermented beverage.

Diversity and ecological footprint of Global Ocean RNA viruses

DNA viruses are increasingly recognized as influencing marine microbes and microbe-mediated biogeochemical cycling. However, little is known about global marine RNA virus diversity, ecology, and ecosystem roles. In this study, we uncover patterns and predictors of marine RNA virus community- and "species"-level diversity and contextualize their ecological impacts from pole to pole. Our analyses revealed four ecological zones, latitudinal and depth diversity patterns, and environmental correlates for RNA viruses. Our findings only partially parallel those of cosampled plankton and show unexpectedly high polar ecological interactions. The influence of RNA viruses on ecosystems appears to be large, as predicted hosts are ecologically important. Moreover, the occurrence of auxiliary metabolic genes indicates that RNA viruses cause reprogramming of diverse host metabolisms, including photosynthesis and carbon cycling, and that RNA virus abundances predict ocean carbon export.

Diversity and origins of bacterial and archaeal viruses on sinking particles reaching the abyssal ocean

Sinking particles and particle-associated microbes influence global biogeochemistry through particulate matter export from the surface to the deep ocean. Despite ongoing studies of particle-associated microbes, viruses in these habitats remain largely unexplored. Whether, where, and which viruses might contribute to particle production and export remain open to investigation. In this study, we analyzed 857 virus population genomes associated with sinking particles collected over three years in sediment traps moored at 4000 m in the North Pacific Subtropical Gyre. Particle-associated viruses here were linked to cellular hosts through matches to bacterial and archaeal metagenome-assembled genome (MAG)-encoded prophages or CRISPR spacers, identifying novel viruses infecting presumptive deep-sea bacteria such as Colwellia, Moritella, and Shewanella. We also identified lytic viruses whose abundances correlated with particulate carbon flux and/or were exported from the photic to abyssal ocean, including cyanophages. Our data are consistent with some of the predicted outcomes of the viral shuttle hypothesis, and further suggest that viral lysis of both autotrophic and heterotrophic prokaryotes may play a role in carbon export. Our analyses revealed the diversity and origins of prevalent viruses found on deep-sea sinking particles and identified prospective viral groups for future investigation into processes that govern particle export in the open ocean.

PHISDetector: A Tool to Detect Diverse In Silico Phage-host Interaction Signals for Virome Studies

Phage-microbe interactions are appealing systems to study coevolution, and have also been increasingly emphasized due to their roles in human health, disease, and the development of novel therapeutics. Phage-microbe interactions leave diverse signals in bacterial and phage genomic sequences, defined as phage-host interaction signals (PHISs), which include clustered regularly interspaced short palindromic repeats (CRISPR) targeting, prophage, and protein-protein interaction signals. In the present study, we developed a novel tool phage-host interaction signal detector (PHISDetector) to predict phage-host interactions by detecting and integrating diverse in silico PHISs, and scoring the probability of phage-host interactions using machine learning models based on PHIS features. We evaluated the performance of PHISDetector on multiple benchmark datasets and application cases. When tested on a dataset of 758 annotated phage-host pairs, PHISDetector yields the prediction accuracies of 0.51 and 0.73 at the species and genus levels, respectively, outperforming other phage-host prediction tools. When applied to on 125,842 metagenomic viral contigs (mVCs) derived from 3042 geographically diverse samples, a detection rate of 54.54% could be achieved. Furthermore, PHISDetector could predict infecting phages for 85.6% of 368 multidrug-resistant (MDR) bacteria and 30% of 454 human gut bacteria obtained from the National Institutes of Health (NIH) Human Microbiome Project (HMP). The PHISDetector can be run either as a web server (http://www.microbiome-bigdata.com/PHISDetector/) for general users to study individual inputs or as a stand-alone version (https://github.com/HIT-ImmunologyLab/PHISDetector) to process massive phage contigs from virome studies.

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.

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.

CHERRY: a Computational metHod for accuratE pRediction of virus-pRokarYotic interactions using a graph encoder-decoder model

Prokaryotic viruses, which infect bacteria and archaea, are key players in microbial communities. Predicting the hosts of prokaryotic viruses helps decipher the dynamic relationship between microbes. Experimental methods for host prediction cannot keep pace with the fast accumulation of sequenced phages. Thus, there is a need for computational host prediction. Despite some promising results, computational host prediction remains a challenge because of the limited known interactions and the sheer amount of sequenced phages by high-throughput sequencing technologies. The state-of-the-art methods can only achieve 43% accuracy at the species level. In this work, we formulate host prediction as link prediction in a knowledge graph that integrates multiple protein and DNA-based sequence features. Our implementation named CHERRY can be applied to predict hosts for newly discovered viruses and to identify viruses infecting targeted bacteria. We demonstrated the utility of CHERRY for both applications and compared its performance with 11 popular host prediction methods. To our best knowledge, CHERRY has the highest accuracy in identifying virus–prokaryote interactions. It outperforms all the existing methods at the species level with an accuracy increase of 37%. In addition, CHERRY’s performance on short contigs is more stable than other tools.

A mixed blessing of viruses in wastewater treatment plants

Activated sludge of wastewater treatment plants harbors a very high diversity of both microorganisms and viruses, wherein the latter control microbial dynamics and metabolisms by infection and lysis of cells. However, it remains poorly understood how viruses impact the biochemical processes of activated sludge, for example in terms of treatment efficiency and pollutant removal. Using metagenomic and metatranscriptomic deep sequencing, the present study recovered thousands of viral sequences from activated sludge samples of three conventional wastewater treatment plants. Gene-sharing network indicated that most of viruses could not be assigned to known viral genera, implying activated sludge as an underexplored reservoir for new viruses and viral diversity. In silico predictions of virus-host linkages demonstrated that infected microbial hosts, mostly belonging to bacteria, were transcriptionally active and able to hydrolyze polymers including starches, celluloses, and proteins. Some viruses encode auxiliary metabolic genes (AMGs) involved in carbon, nitrogen, and sulfur cycling, and antibiotic resistance genes (ARGs) for resistance to multiple drugs. The virus-encoded AMGs may enhance the biodegradation of contaminants like starches and celluloses, suggesting a positive role for viruses in strengthening the performance of activated sludge. However, ARGs would be disseminated to different microorganisms using viruses as gene shuttles, demonstrating the possibility for viruses to facilitate the spread of antibiotic resistance in the environment. Collectively, this study highlights the mixed blessing of viruses in wastewater treatment plants, and deciphers how they manipulate the biochemical processes in the activated sludge, with implications for both environmental protection and ecosystem security.

A Synthesis of Viral Contribution to Marine Nitrogen Cycling

Nitrogen is an essential component of major cellular macromolecules, such as DNA and proteins. Its bioavailability has a fundamental influence on the primary production of both terrestrial and oceanic ecosystems. Diverse marine microbes consume nitrogen, while only a limited taxon could replenish it, leaving nitrogen one of the most deficient nutrients in the ocean. A variety of microbes are involved in complex biogeochemical transformations of nitrogen compounds, and their ecological functions might be regulated by viruses in different manners. First and foremost, viruses drive marine nitrogen flow via host cell lysis, releasing abundant organic nitrogen into the surrounding environment. Secondly, viruses can also participate in the marine nitrogen cycle by expressing auxiliary metabolic genes (AMGs) to modulate host nitrogen metabolic pathways, such as nitrification, denitrification, anammox, and nitrogen transmembrane transport. Additionally, viruses also serve as a considerable reservoir of nitrogen element. The efficient turnover of viruses fundamentally promotes nitrogen flow in the oceans. In this review, we summarize viral contributions in the marine nitrogen cycling in different aspects and discuss challenges and issues based on recent discoveries of novel viruses involved in different processes of nitrogen biotransformation.

Expanding the Database of Signal-Anchor-Release Domain Endolysins Through Metagenomics

Endolysins are bacteriophage-derived lytic enzymes with antimicrobial activity. The action of endolysins against Gram-negative bacteria remains a challenge due to the physical protection of the outer membrane. However, recent research has demonstrated that signal-anchor-release (SAR) endolysins permeate the outer membrane of Gram-negative bacteria. This study investigates 2628 putative endolysin genes identified in 183,298 bacteriophage genomes. Previously, bioinformatic approaches resulted in a database of 66 SAR endolysins. This manuscript almost doubles the list with 53 additional SAR endolysin candidates. Forty-eight of the putative SAR endolysins described in this study contained one muramidase catalytic domain, and five included additional cell wall-binding domains at the C-terminus. For the moment, SAR domains are found in four protein families: glycoside hydrolase family 19 (GH19), glycoside hydrolase family 24 (GH24), glycoside hydrolase family 25 (GH25), and glycoside hydrolase family 108 (GH108). These SAR lysis are clustered in eight groups based on biochemical properties and domain presence/absence. Therefore, in this study, we expand the arsenal of endolysin candidates that might act against Gram-negative bacteria and develop a consult database for antimicrobial proteins derived from bacteriophages.

Patterns and ecological drivers of viral communities in acid mine drainage sediments across Southern China

Recent advances in environmental genomics have provided unprecedented opportunities for the investigation of viruses in natural settings. Yet, our knowledge of viral biogeographic patterns and the corresponding drivers is still limited. Here, we perform metagenomic deep sequencing on 90 acid mine drainage (AMD) sediments sampled across Southern China and examine the biogeography of viruses in this extreme environment. The results demonstrate that prokaryotic communities dictate viral taxonomic and functional diversity, abundance and structure, whereas other factors especially latitude and mean annual temperature also impact viral populations and functions. In silico predictions highlight lineage-specific virus-host abundance ratios and richness-dependent virus-host interaction structure. Further functional analyses reveal important roles of environmental conditions and horizontal gene transfers in shaping viral auxiliary metabolic genes potentially involved in phosphorus assimilation. Our findings underscore the importance of both abiotic and biotic factors in predicting the taxonomic and functional biogeographic dynamics of viruses in the AMD sediments.

The biogeography of viral communities in extreme environments remains understudied. Here, the authors use metagenomic sequencing on 90 acid mine drainage sediments sampled across Southern China, showing the predominant effects of prokaryotic communities and the influence of environmental variables on viral taxonomy and function.

Newly identified HMO-2011-type phages reveal genomic diversity and biogeographic distributions of this marine viral group

Viruses play critical roles in influencing biogeochemical cycles and adjusting host mortality, population structure, physiology, and evolution in the ocean. Marine viral communities are composed of numerous genetically distinct subfamily/genus-level viral groups. Among currently identified viral groups, the HMO-2011-type group is known to be dominant and broadly distributed. However, only four HMO-2011-type cultivated representatives that infect marine SAR116 and Roseobacter strains have been reported to date, and the genetic diversity, potential hosts, and ecology of this group remain poorly elucidated. Here, we present the genomes of seven HMO-2011-type phages that were isolated using four Roseobacter strains and one SAR11 strain, as well as additional 207 HMO-2011-type metagenomic viral genomes (MVGs) identified from various marine viromes. Phylogenomic and shared-gene analyses revealed that the HMO-2011-type group is a subfamily-level group comprising at least 10 discernible genus-level subgroups. Moreover, >2000 HMO-2011-type DNA polymerase sequences were identified, and the DNA polymerase phylogeny also revealed that the HMO-2011-type group contains diverse subgroups and is globally distributed. Metagenomic read-mapping results further showed that most HMO-2011-type phages are prevalent in global oceans and display distinct geographic distributions, with the distribution of most HMO-2011-type phages being associated with temperature. Lastly, we found that members in subgroup IX, represented by pelagiphage HTVC033P, were among the most abundant HMO-2011-type phages, which implies that SAR11 bacteria are crucial hosts for this viral group. In summary, our findings substantially expand current knowledge regarding the phylogenetic diversity, evolution, and distribution of HMO-2011-type phages, highlighting HMO-2011-type phages as major ecological agents that can infect certain key bacterial groups.

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.

Distribution Characteristics of Soil Viruses Under Different Precipitation Gradients on the Qinghai-Tibet Plateau

Viruses are extremely abundant in the soil environment and have potential roles in impacting on microbial population, evolution, and nutrient biogeochemical cycles. However, how environment and climate changes affect soil viruses is still poorly understood. Here, a metagenomic approach was used to investigate the distribution, diversity, and potential biogeochemical impacts of DNA viruses in 12 grassland soils under three precipitation gradients on the Qinghai-Tibet Plateau, which is one of the most sensitive areas to climate change. A total of 557 viral operational taxonomic units were obtained, spanning 152 viral families from the 30 metagenomes. Both virus-like particles (VLPs) and microbial abundance increased with average annual precipitation. A significant positive correlation of VLP counts was observed with soil water content, total carbon, total nitrogen, soil organic matter, and total phosphorus. Among these biological and abiotic factors, SWC mainly contributed to the variability in VLP abundance. The order Caudovirales (70.1% of the identified viral order) was the predominant viral type in soils from the Qinghai-Tibet Plateau, with the Siphoviridae family being the most abundant. Remarkably, abundant auxiliary carbohydrate-active enzyme (CAZyme) genes represented by glycoside hydrolases were identified, indicating that soil viruses may play a potential role in the carbon cycle on the Qinghai-Tibet Plateau. There were more diverse hosts and abundant CAZyme genes in soil with moderate precipitation. Our study provides a strong evidence that changes in precipitation impact not only viral abundance and virus–host interactions in soil but also the viral functional potential, especially carbon cycling.

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.

Picorna-Like Viruses of the Havel River, Germany

To improve the understanding of the virome diversity of riverine ecosystems in metropolitan areas, a metagenome analysis was performed with water collected in June 2018 from the river Havel in Berlin, Germany. After enrichment of virus particles and RNA extraction, paired-end Illumina sequencing was conducted and assignment to virus groups and families was performed. This paper focuses on picorna-like viruses, the most diverse and abundant group of viruses with impact on human, animal, and environmental health. Here, we describe altogether 166 viral sequences ranging in size from 1 to 11.5 kb. The 71 almost complete genomes are comprised of one candidate iflavirus, one picornavirus, two polycipiviruses, 27 marnaviruses, 27 dicistro-like viruses, and 13 untypeable viruses. Many partial picorna-like virus sequences up to 10.2 kb were also investigated. The sequences of the Havel picorna-like viruses represent genomes of seven of eight so far known Picornavirales families. Detection of numerous distantly related dicistroviruses suggests the existence of additional, yet unexplored virus groups with dicistronic genomes, including few viruses with unusual genome layout. Of special interest is a clade of dicistronic viruses with capsid protein-encoding sequences at the 5′-end of the genome. Also, monocistronic viruses with similarity of their polymerase and capsid proteins to those of dicistroviruses are interesting. A second protein with NTP-binding site present in the polyprotein of solinviviruses and related viruses needs further attention. The results underline the importance to study the viromes of fluvial ecosystems. So far acknowledged marnaviruses have been isolated from marine organisms. However, the present study and available sequence data suggest that rivers and limnic habitats are relevant ecosystems with circulation of marnaviruses as well as a plethora of unknown picorna-like viruses.

Novel Viral DNA Polymerases From Metagenomes Suggest Genomic Sources of Strand-Displacing Biochemical Phenotypes

Viruses are the most abundant and diverse biological entities on the planet and constitute a significant proportion of Earth's genetic diversity. Most of this diversity is not represented by isolated viral-host systems and has only been observed through sequencing of viral metagenomes (viromes) from environmental samples. Viromes provide snapshots of viral genetic potential, and a wealth of information on viral community ecology. These data also provide opportunities for exploring the biochemistry of novel viral enzymes. The in vitro biochemical characteristics of novel viral DNA polymerases were explored, testing hypothesized differences in polymerase biochemistry according to protein sequence phylogeny. Forty-eight viral DNA Polymerase I (PolA) proteins from estuarine viromes, hot spring metagenomes, and reference viruses, encompassing a broad representation of currently known diversity, were synthesized, expressed, and purified. Novel functionality was shown in multiple PolAs. Intriguingly, some of the estuarine viral polymerases demonstrated moderate to strong innate DNA strand displacement activity at high enzyme concentration. Strand-displacing polymerases have important technological applications where isothermal reactions are desirable. Bioinformatic investigation of genes neighboring these strand displacing polymerases found associations with SNF2 helicase-associated proteins. The specific function of SNF2 family enzymes is unknown for prokaryotes and viruses. In eukaryotes, SNF2 enzymes have chromatin remodeling functions but do not separate nucleic acid strands. This suggests the strand separation function may be fulfilled by the DNA polymerase for viruses carrying SNF2 helicase-associated proteins. Biochemical data elucidated from this study expands understanding of the biology and ecological behavior of unknown viruses. Moreover, given the numerous biotechnological applications of viral DNA polymerases, novel viral polymerases discovered within viromes may be a rich source of biological material for further in vitro DNA amplification advancements.

The "Regulator" Function of Viruses on Ecosystem Carbon Cycling in the Anthropocene

Viruses act as "regulators" of the global carbon cycle because they impact the material cycles and energy flows of food webs and the microbial loop. The average contribution of viruses to the Earth ecosystem carbon cycle is 8.6‰, of which its contribution to marine ecosystems (1.4‰) is less than its contribution to terrestrial (6.7‰) and freshwater (17.8‰) ecosystems. Over the past 2,000 years, anthropogenic activities and climate change have gradually altered the regulatory role of viruses in ecosystem carbon cycling processes. This has been particularly conspicuous over the past 200 years due to rapid industrialization and attendant population growth. The progressive acceleration of the spread and reproduction of viruses may subsequently accelerate the global C cycle.

The hidden world within plants: metatranscriptomics unveils the complexity of wood microbiomes

The importance of plants as complex entities influenced by genomes of the associated microorganisms is now seen as a new source of variability for a more sustainable agriculture, also in the light of ongoing climate change. For this reason, we investigated through metatranscriptomics whether the taxa profile and behaviour of microbial communities associated with the wood of 20-year-old grapevine plants are influenced by the health status of the host. We report for the first time a metatranscriptome from a complex tissue in a real environment, highlighting that this approach is able to define the microbial community better than referenced transcriptomic approaches. In parallel, the use of total RNA enabled the identification of bacterial taxa in healthy samples that, once isolated from the original wood tissue, displayed potential biocontrol activities against a wood-degrading fungal taxon. Furthermore, we revealed an unprecedented high number of new viral entities (~120 new viral species among 180 identified) associated with a single and limited environment and with potential impact on the whole holobiont. Taken together, our results suggest a complex multitrophic interaction in which the viral community also plays a crucial role in raising new ecological questions for the exploitation of microbial-assisted sustainable agriculture.

Heterogeneity of soil bacterial and bacteriophage communities in three rice agroecosystems and potential impacts of bacteriophage on nutrient cycling

BACKGROUND

As genetic entities infecting and replicating only in bacteria, bacteriophages can regulate the community structure and functions of their host bacteria. The ecological roles of bacteriophages in aquatic and forest environments have been widely explored, but those in agroecosystems remains limited. Here, we used metagenomic sequencing to analyze the diversity and interactions of bacteriophages and their host bacteria in soils from three typical rice agroecosystems in China: double cropping in Guangzhou, southern China, rice-wheat rotation cropping in Nanjing, eastern China and early maturing single cropping in Jiamusi, northeastern China. Enterobacter phage-NJ was isolated and its functions on soil nitrogen cycling and effect on soil bacterial community structure were verified in pot inoculation experiments and 16S rRNA gene sequencing.

RESULTS

Soil bacterial and viral diversity and predicted functions varied among the three agroecosystems. Genes detected in communities from the three agroecosystems were associated with typical functions: soil bacteria in Jiamusi were significantly enriched in genes related to carbohydrate metabolism, in Nanjing with xenobiotics biodegradation and metabolism, and in Guangzhou with virulence factors and scarce in secondary metabolite biosynthesis, which might lead to a significant occurrence of rice bacterial diseases. The virus community structure varies significantly among the three ecosystems, only 13.39% of the total viral species were shared by the three rice agroecosystems, 59.56% of the viral species were specific to one agroecosystem. Notably, over-represented auxiliary carbohydrate-active enzyme (CAZyme) genes were identified in the viruses, which might assist host bacteria in metabolizing carbon, and 67.43% of these genes were present in Jiamusi. In bacteriophage isolation and inoculation experiments, Enterobacter bacteriophage-NJ reduced the nitrogen fixation capacity of soil by lysing N-fixing host bacteria and changed the soil bacterial diversity and community structure.

CONCLUSION

Our results showed that diversity and function predicted of paddy soil bacteria and viruses varied in the three agroecosystems. Soil bacteriophages can affect nutrient cycling by boosting host metabolism through the carried auxiliary metabolic genes (AMGs) and lysing the host bacteria that are involved in biogeochemical cycles. These findings form a basis for better understanding bacterial and bacteriophage diversity in different rice agroecosystems, laying a solid foundation for further studies of soil microbial communities that support ecofriendly production of healthy rice.