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

Diversity of Bathyarchaeia viruses in metagenomes and virus-encoded CRISPR system components

Bathyarchaeia represent a class of archaea common and abundant in sedimentary ecosystems. Here we report 56 metagenome-assembled genomes of Bathyarchaeia viruses identified in metagenomes from different environments. Gene sharing network and phylogenomic analyses led to the proposal of four virus families, including viruses of the realms Duplodnaviria and Adnaviria, and archaea-specific spindle-shaped viruses. Genomic analyses uncovered diverse CRISPR elements in these viruses. Viruses of the proposed family "Fuxiviridae" harbor an atypical Type IV-B CRISPR-Cas system and a Cas4 protein that might interfere with host immunity. Viruses of the family "Chiyouviridae" encode a Cas2-like endonuclease and two mini-CRISPR arrays, one with a repeat identical to that in the host CRISPR array, potentially allowing the virus to recruit the host CRISPR adaptation machinery to acquire spacers that could contribute to competition with other mobile genetic elements or to inhibit host defenses. These findings present an outline of the Bathyarchaeia virome and offer a glimpse into their counter-defense mechanisms.

Exploring the Archaeal Virosphere by Metagenomics

During the past decade, environmental research has demonstrated that archaea are abundant and widespread in nature and play important ecological roles at a global scale. Currently, however, the majority of archaeal lineages cannot be cultivated under laboratory conditions and are known exclusively or nearly exclusively through metagenomics. A similar trend extends to the archaeal virosphere, where isolated representatives are available for a handful of model archaeal virus-host systems. Viral metagenomics provides an alternative way to circumvent the limitations of culture-based virus discovery and offers insight into the diversity, distribution, and environmental impact of uncultured archaeal viruses. Presently, metagenomics approaches have been successfully applied to explore the viromes associated with various lineages of extremophilic and mesophilic archaea, including Asgard archaea (Asgardarchaeota), ANME-1 archaea (Methanophagales), thaumarchaea (Nitrososphaeria), altiarchaea (Altiarchaeota), and marine group II archaea (Poseidoniales). Here, we provide an overview of methods widely used in archaeal virus metagenomics, covering metavirome preparation, genome annotation, phylogenetic and phylogenomic analyses, and archaeal host assignment. We hope that this summary will contribute to further exploration and characterization of the enigmatic archaeal virome lurking in diverse environments.

Phylogenomic analysis expands the known repertoire of single-stranded DNA viruses in benthic zones of the South Indian Ocean

Single-stranded (ss) DNA viruses are ubiquitous and constitute some of the most diverse entities on Earth. Most studies have focused on ssDNA viruses from terrestrial environments resulting in a significant deficit in benthic ecosystems including aphotic zones of the South Indian Ocean (SIO). Here, we assess the diversity and phylogeny of ssDNA in deep waters of the SIO using a combination of established viral taxonomy tools and a Hidden Markov Model based approach. Replication initiator protein-associated (Rep) phylogenetic reconstruction and sequence similarity networks were used to show that the SIO hosts divergent and as yet unknown circular Rep-encoding ssDNA viruses. Several sequences appear to represent entirely novel families, expanding the repertoire of known ssDNA viruses. Results suggest that a small proportion of these viruses may be circular genetic elements, which may strongly influence the diversity of both eukaryotes and prokaryotes in the SIO. Taken together, our data show that the SIO harbours a diverse assortment of previously unknown ssDNA viruses. Due to their potential to infect a variety of hosts, these viruses may be crucial for marine nutrient recycling through their influence of the biological carbon pump.

Host translation machinery is not a barrier to phages that interact with both CPR and non-CPR bacteria

IMPORTANCE

Here, we profiled putative phages of Saccharibacteria, which are of particular importance as Saccharibacteria influence some human oral diseases. We additionally profiled putative phages of Gracilibacteria and Absconditabacteria, two Candidate Phyla Radiation (CPR) lineages of interest given their use of an alternative genetic code. Among the phages identified in this study, some are targeted by spacers from both CPR and non-CPR bacteria and others by both bacteria that use the standard genetic code as well as bacteria that use an alternative genetic code. These findings represent new insights into possible phage replication strategies and have relevance for phage therapies that seek to manipulate microbiomes containing CPR bacteria.

Zea mays genotype influences microbial and viral rhizobiome community structure

Plant genotype is recognized to contribute to variations in microbial community structure in the rhizosphere, soil adherent to roots. However, the extent to which the viral community varies has remained poorly understood and has the potential to contribute to variation in soil microbial communities. Here we cultivated replicates of two Zea mays genotypes, parviglumis and B73, in a greenhouse and harvested the rhizobiome (rhizoplane and rhizosphere) to identify the abundance of cells and viruses as well as rhizobiome microbial and viral community using 16S rRNA gene amplicon sequencing and genome resolved metagenomics. Our results demonstrated that viruses exceeded microbial abundance in the rhizobiome of parviglumis and B73 with a significant variation in both the microbial and viral community between the two genotypes. Of the viral contigs identified only 4.5% (n = 7) of total viral contigs were shared between the two genotypes, demonstrating that plants even at the level of genotype can significantly alter the surrounding soil viral community. An auxiliary metabolic gene associated with glycoside hydrolase (GH5) degradation was identified in one viral metagenome-assembled genome (vOTU) identified in the B73 rhizobiome infecting Propionibacteriaceae (Actinobacteriota) further demonstrating the viral contribution in metabolic potential for carbohydrate degradation and carbon cycling in the rhizosphere. This variation demonstrates the potential of plant genotype to contribute to microbial and viral heterogeneity in soil systems and harbors genes capable of contributing to carbon cycling in the rhizosphere.

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

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

Enhancing plant growth in biofertilizer-amended soil through nitrogen-transforming microbial communities

Biofertilizers have immense potential for enhancing agricultural productivity. However, there is still a need for clarification regarding the specific mechanisms through which these biofertilizers improve soil properties and stimulate plant growth. In this research, a bacterial agent was utilized to enhance plant growth and investigate the microbial modulation mechanism of soil nutrient turnover using metagenomic technology. The results demonstrated a significant increase in soil fast-acting nitrogen (by 46.7%) and fast-acting phosphorus (by 88.6%) upon application of the bacterial agent. This finding suggests that stimulated soil microbes contribute to enhanced nutrient transformation, ultimately leading to improved plant growth. Furthermore, the application of the bacterial agent had a notable impact on the accumulation of key genes involved in nitrogen cycling. Notably, it enhanced nitrification genes (amo, hao, and nar), while denitrification genes (nir and nor) showed a slight decrease. This indicates that ammonium oxidation may be the primary pathway for increasing fast-acting nitrogen in soils. Additionally, the bacterial agent influenced the composition and functional structure of the soil microbial community. Moreover, the metagenome-assembled genomes (MAGs) obtained from the soil microbial communities exhibited complementary metabolic processes, suggesting mutual nutrient exchange. These MAGs contained widely distributed and highly abundant genes encoding plant growth promotion (PGP) traits. These findings emphasize how soil microbial communities can enhance vegetation growth by increasing nutrient availability and regulating plant hormone production. This effect can be further enhanced by introducing inoculated microbial agents. In conclusion, this study provides novel insights into the mechanisms underlying the beneficial effects of biofertilizers on soil properties and plant growth. The significant increase in nutrient availability, modulation of key genes involved in nitrogen cycling, and the presence of MAGs encoding PGP traits highlight the potential of biofertilizers to improve agricultural practices. These findings have important implications for enhancing agricultural sustainability and productivity, with positive societal and environmental impacts.

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

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

Metagenomic data from surface seawater of the east coast of South Korea

The East Sea, also known as the Sea of Japan, is a marginal sea located in the western Pacific Ocean, displaying comparable characteristics to Earth's oceans, thereby meriting its recognition as a "miniature ocean". The East Sea exhibits a range of annually-recurring biogeochemical features in accordance with seasonal fluctuations, such as phytoplankton blooms during the spring and autumn seasons. Despite ongoing monitoring efforts focused on water quality and physicochemical parameters, the investigation of prokaryotic assemblages in the East Sea, encompassing seasonal variations, has been infrequently pursued. Here, we present a monthly time-series metagenomic dataset spanning a one-year period in 2009, obtained from surface (10 m) seawater samples collected off the coast of the East Sea. The dataset encompasses 12 metagenomes, amounting 195 Gbp, with 14.73-22.52 Gbp per sample. This dataset is accompanied by concurrently measured physicochemical parameters. Our anticipation is that these metagenomes will facilitate extensive investigations aimed at elucidating various aspects of the marine microbial ecosystems in the East Sea.

Diversity of Bathyarchaeia viruses in metagenomes and virus-encoded CRISPR system components

Bathyarchaeia represent a class of archaea common and abundant in sedimentary ecosystems. The virome of Bathyarchaeia so far has not been characterized. Here we report 56 metagenome-assembled genomes of Bathyarchaeia viruses identified in metagenomes from different environments. Gene sharing network and phylogenomic analyses led to the proposal of four virus families, including viruses of the realms Duplodnaviria and Adnaviria, and archaea-specific spindle-shaped viruses. Genomic analyses uncovered diverse CRISPR elements in these viruses. Viruses of the proposed family 'Fuxiviridae' harbor an atypical type IV-B CRISPR-Cas system and a Cas4 protein that might interfere with host immunity. Viruses of the family 'Chiyouviridae' encode a Cas2-like endonuclease and two mini-CRISPR arrays, one with a repeat identical to that in the host CRISPR array, potentially allowing the virus to recruit the host CRISPR adaptation machinery to acquire spacers that could contribute to competition with other mobile genetic elements or to inhibition of host defenses. These findings present an outline of the Bathyarchaeia virome and offer a glimpse into their counter-defense mechanisms.

MArVD2: a machine learning enhanced tool to discriminate between archaeal and bacterial viruses in viral datasets

Our knowledge of viral sequence space has exploded with advancing sequencing technologies and large-scale sampling and analytical efforts. Though archaea are important and abundant prokaryotes in many systems, our knowledge of archaeal viruses outside of extreme environments is limited. This largely stems from the lack of a robust, high-throughput, and systematic way to distinguish between bacterial and archaeal viruses in datasets of curated viruses. Here we upgrade our prior text-based tool (MArVD) via training and testing a random forest machine learning algorithm against a newly curated dataset of archaeal viruses. After optimization, MArVD2 presented a significant improvement over its predecessor in terms of scalability, usability, and flexibility, and will allow user-defined custom training datasets as archaeal virus discovery progresses. Benchmarking showed that a model trained with viral sequences from the hypersaline, marine, and hot spring environments correctly classified 85% of the archaeal viruses with a false detection rate below 2% using a random forest prediction threshold of 80% in a separate benchmarking dataset from the same habitats.

ViWrap: A modular pipeline to identify, bin, classify, and predict viral-host relationships for viruses from metagenomes

Viruses are increasingly being recognized as important components of human and environmental microbiomes. However, viruses in microbiomes remain difficult to study because of the difficulty in culturing them and the lack of sufficient model systems. As a result, computational methods for identifying and analyzing uncultivated viral genomes from metagenomes have attracted significant attention. Such bioinformatics approaches facilitate the screening of viruses from enormous sequencing datasets originating from various environments. Though many tools and databases have been developed for advancing the study of viruses from metagenomes, there is a lack of integrated tools enabling a comprehensive workflow and analyses platform encompassing all the diverse segments of virus studies. Here, we developed ViWrap, a modular pipeline written in Python. ViWrap combines the power of multiple tools into a single platform to enable various steps of virus analysis, including identification, annotation, genome binning, species- and genus-level clustering, assignment of taxonomy, prediction of hosts, characterization of genome quality, comprehensive summaries, and intuitive visualization of results. Overall, ViWrap enables a standardized and reproducible pipeline for both extensive and stringent characterization of viruses from metagenomes, viromes, and microbial genomes. Our approach has flexibility in using various options for diverse applications and scenarios, and its modular structure can be easily amended with additional functions as necessary. ViWrap is designed to be easily and widely used to study viruses in human and environmental systems. ViWrap is publicly available via GitHub (https://github.com/AnantharamanLab/ViWrap). A detailed description of the software, its usage, and interpretation of results can be found on the website.

Beyond pathogenesis: Detecting the full spectrum of ecological interactions in the virosphere

The public perception of viruses has historically been negative. We are now at a stage where the development of tools to study viruses is at an all-time high, but society's perception of viruses is at an all-time low. The literature regarding viral interactions has been skewed towards negative (i.e., pathogenic) symbioses, whereas viral mutualisms remain relatively underexplored. Viral interactions with their hosts are complex and some non-pathogenic viruses could have potential benefits to society. However, viral research is seldom designed to identify viral mutualists, a gap that merits considering new experimental designs. Determining whether antagonisms, mutualisms, and commensalisms are equally common ecological strategies requires more balanced research efforts that characterize the full spectrum of viral interactions.

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.

The smallest in the deepest: the enigmatic role of viruses in the deep biosphere

It is commonly recognized that viruses control the composition, metabolism, and evolutionary trajectories of prokaryotic communities, with resulting vital feedback on ecosystem functioning and nutrient cycling in a wide range of ecosystems. Although the deep biosphere has been estimated to be the largest reservoir for viruses and their prokaryotic hosts, the biology and ecology of viruses therein remain poorly understood. The deep virosphere is an enigmatic field of study in which many critical questions are still to be answered. Is the deep virosphere simply a repository for deeply preserved, non-functioning virus particles? Or are deep viruses infectious agents that can readily infect suitable hosts and subsequently shape microbial populations and nutrient cycling? Can the cellular content released by viral lysis, and even the organic structures of virions themselves, serve as the source of bioavailable nutrients for microbial activity in the deep biosphere as in other ecosystems? In this review, we synthesize our current knowledge of viruses in the deep biosphere and seek to identify topics with the potential for substantial discoveries in the future.

Unexpected complexity of the ammonia monooxygenase in archaea

Ammonia oxidation, as the first step of nitrification, constitutes a critical process in the global nitrogen cycle. However, fundamental knowledge of its key enzyme, the copper-dependent ammonia monooxygenase, is lacking, in particular for the environmentally abundant ammonia-oxidizing archaea (AOA). Here the structure of the enzyme is investigated by blue-native gel electrophoresis and proteomics from native membrane complexes of two AOA. Besides the known AmoABC subunits and the earlier predicted AmoX, two new protein subunits, AmoY and AmoZ, were identified. They are unique to AOA, highly conserved and co-regulated, and their genes are linked to other AMO subunit genes in streamlined AOA genomes. Modeling and in-gel cross-link approaches support an overall protomer structure similar to the distantly related bacterial particulate methane monooxygenase but also reveals clear differences in extracellular domains of the enzyme. These data open avenues for further structure-function studies of this ecologically important nitrification complex.

Abundance, classification and genetic potential of Thaumarchaeota in metagenomes of European agricultural soils: a meta-analysis

BACKGROUND

For a sustainable production of food, research on agricultural soil microbial communities is inevitable. Due to its immense complexity, soil is still some kind of black box. Soil study designs for identifying microbiome members of relevance have various scopes and focus on particular environmental factors. To identify common features of soil microbiomes, data from multiple studies should be compiled and processed. Taxonomic compositions and functional capabilities of microbial communities associated with soils and plants have been identified and characterized in the past few decades. From a fertile Loess-Chernozem-type soil located in Germany, metagenomically assembled genomes (MAGs) classified as members of the phylum Thaumarchaeota/Thermoproteota were obtained. These possibly represent keystone agricultural soil community members encoding functions of relevance for soil fertility and plant health. Their importance for the analyzed microbiomes is corroborated by the fact that they were predicted to contribute to the cycling of nitrogen, feature the genetic potential to fix carbon dioxide and possess genes with predicted functions in plant-growth-promotion (PGP). To expand the knowledge on soil community members belonging to the phylum Thaumarchaeota, we conducted a meta-analysis integrating primary studies on European agricultural soil microbiomes.

RESULTS

Taxonomic classification of the selected soil metagenomes revealed the shared agricultural soil core microbiome of European soils from 19 locations. Metadata reporting was heterogeneous between the different studies. According to the available metadata, we separated the data into 68 treatments. The phylum Thaumarchaeota is part of the core microbiome and represents a major constituent of the archaeal subcommunities in all European agricultural soils. At a higher taxonomic resolution, 2074 genera constituted the core microbiome. We observed that viral genera strongly contribute to variation in taxonomic profiles. By binning of metagenomically assembled contigs, Thaumarchaeota MAGs could be recovered from several European soil metagenomes. Notably, many of them were classified as members of the family Nitrososphaeraceae, highlighting the importance of this family for agricultural soils. The specific Loess-Chernozem Thaumarchaeota MAGs were most abundant in their original soil, but also seem to be of importance in other agricultural soil microbial communities. Metabolic reconstruction of Switzerland_1_MAG_2 revealed its genetic potential i.a. regarding carbon dioxide (CO[Formula: see text]) fixation, ammonia oxidation, exopolysaccharide production and a beneficial effect on plant growth. Similar genetic features were also present in other reconstructed MAGs. Three Nitrososphaeraceae MAGs are all most likely members of a so far unknown genus.

CONCLUSIONS

On a broad view, European agricultural soil microbiomes are similarly structured. Differences in community structure were observable, although analysis was complicated by heterogeneity in metadata recording. Our study highlights the need for standardized metadata reporting and the benefits of networking open data. Future soil sequencing studies should also consider high sequencing depths in order to enable reconstruction of genome bins. Intriguingly, the family Nitrososphaeraceae commonly seems to be of importance in agricultural microbiomes.

Awaking the dormant virome in the rhizosphere

The rhizosphere is a vital soil compartment providing key plant-beneficial functions. However, little is known about the mechanisms driving viral diversity in the rhizosphere. Viruses can establish lytic or lysogenic interactions with their bacterial hosts. In the latter, they assume a dormant state integrated in the host genome and can be awakened by different perturbations that impact host cell physiology, triggering a viral bloom, which is potentially a fundamental mechanism driving soil viral diversity, as 22%-68% of soil bacteria are predicted to harbour dormant viruses. Here we assessed the viral bloom response in rhizospheric viromes by exposing them to three contrasting soil perturbation agents: earthworms, herbicide and antibiotic pollutant. The viromes were next screened for rhizosphere-relevant genes and also used as inoculant on microcosms incubations to test their impacts on pristine microbiomes. Our results show that while post-perturbation viromes diverged from control conditions, viral communities exposed to both herbicide and antibiotic pollutant were more similar to each other than those influenced by earthworms. The latter also favoured an increase in viral populations harbouring genes involved in plant-beneficial functions. Post-perturbation viromes inoculated on soil microcosms changed the diversity of pristine microbiomes, suggesting that viromes are important components of the soil ecological memory driving eco-evolutionary processes that determine future microbiome trajectories according to past events. Our findings demonstrate that viromes are active players in the rhizosphere and need to be considered in efforts to understand and control the microbial processes towards sustainable crop production.

ViWrap: A modular pipeline to identify, bin, classify, and predict viral-host relationships for viruses from metagenomes

Viruses are increasingly being recognized as important components of human and environmental microbiomes. However, viruses in microbiomes remain difficult to study because of difficulty in culturing them and the lack of sufficient model systems. As a result, computational methods for identifying and analyzing uncultivated viral genomes from metagenomes have attracted significant attention. Such bioinformatics approaches facilitate screening of viruses from enormous sequencing datasets originating from various environments. Though many tools and databases have been developed for advancing the study of viruses from metagenomes, there is a lack of integrated tools enabling a comprehensive workflow and analyses platform encompassing all the diverse segments of virus studies. Here, we developed ViWrap, a modular pipeline written in Python. ViWrap combines the power of multiple tools into a single platform to enable various steps of virus analysis including identification, annotation, genome binning, species- and genus-level clustering, assignment of taxonomy, prediction of hosts, characterization of genome quality, comprehensive summaries, and intuitive visualization of results. Overall, ViWrap enables a standardized and reproducible pipeline for both extensive and stringent characterization of viruses from metagenomes, viromes, and microbial genomes. Our approach has flexibility in using various options for diverse applications and scenarios, and its modular structure can be easily amended with additional functions as necessary. ViWrap is designed to be easily and widely used to study viruses in human and environmental systems. ViWrap is publicly available via GitHub (https://github.com/AnantharamanLab/ViWrap). A detailed description of the software, its usage, and interpretation of results can be found on the website.

Diverse viruses of marine archaea discovered using metagenomics

During the past decade, metagenomics became a method of choice for the discovery of novel viruses. However, host assignment for uncultured viruses remains challenging, especially for archaeal viruses, which are grossly undersampled compared to viruses of bacteria and eukaryotes. Here, we assessed the utility of CRISPR spacer targeting, tRNA gene matching and homology searches for viral signature proteins, such as major capsid proteins, for the assignment of archaeal hosts and validated these approaches on metaviromes from Yangshan Harbor (YSH). We report 35 new genomes of viruses which could be confidently assigned to hosts representing diverse lineages of marine archaea. We show that the archaeal YSH virome is highly diverse, with some viruses enriching the previously described virus groups, such as magroviruses of Marine Group II Archaea (Poseidoniales), and others representing novel groups of marine archaeal viruses. Metagenomic recruitment of Tara Oceans datasets on the YSH viral genomes demonstrated the presence of YSH Poseidoniales and Nitrososphaeria viruses in the global oceans, but also revealed the endemic YSH-specific viral lineages. Furthermore, our results highlight the relationship between the soil and marine thaumarchaeal viruses. We propose three new families within the class Caudoviricetes for the classification of the five complete viral genomes predicted to replicate in marine Poseidoniales and Nitrososphaeria, two ecologically important and widespread archaeal groups. This study illustrates the utility of viral metagenomics in exploring the archaeal virome and provides new insights into the diversity, distribution and evolution of marine archaeal viruses.

Propagation of viral genomes by replicating ammonia-oxidising archaea during soil nitrification

Ammonia-oxidising archaea (AOA) are a ubiquitous component of microbial communities and dominate the first stage of nitrification in some soils. While we are beginning to understand soil virus dynamics, we have no knowledge of the composition or activity of those infecting nitrifiers or their potential to influence processes. This study aimed to characterise viruses having infected autotrophic AOA in two nitrifying soils of contrasting pH by following transfer of assimilated CO₂-derived ¹³C from host to virus via DNA stable-isotope probing and metagenomic analysis. Incorporation of ¹³C into low GC mol% AOA and virus genomes increased DNA buoyant density in CsCl gradients but resulted in co-migration with dominant non-enriched high GC mol% genomes, reducing sequencing depth and contig assembly. We therefore developed a hybrid approach where AOA and virus genomes were assembled from low buoyant density DNA with subsequent mapping of ¹³C isotopically enriched high buoyant density DNA reads to identify activity of AOA. Metagenome-assembled genomes were different between the two soils and represented a broad diversity of active populations. Sixty-four AOA-infecting viral operational taxonomic units (vOTUs) were identified with no clear relatedness to previously characterised prokaryote viruses. These vOTUs were also distinct between soils, with 42% enriched in ¹³C derived from hosts. The majority were predicted as capable of lysogeny and auxiliary metabolic genes included an AOA-specific multicopper oxidase suggesting infection may augment copper uptake essential for central metabolic functioning. These findings indicate virus infection of AOA may be a frequent process during nitrification with potential to influence host physiology and activity.

Biogeochemical sulfur cycling of virus auxiliary metabolic genes involved in Napahai plateau wetland

Virus plays important roles in regulating microbial community structure, horizontal gene transfer, and promoting biological evolution, also augmenting host metabolism during infection via the expression of auxiliary metabolic genes (AMGs), and thus affect biogeochemical cycling in the oceans. As the "kidney of the earth," wetlands have rich biodiversity and abundant resources. Based on metagenomic data, 10 AMGs associated with sulfur cycling, i.e., tusA, moaD, dsrE, soxA, soxB, soxC, soxD, soxX, soxY, and soxZ, were analyzed in Napahai plateau wetland. The phylogenetic trees of AMGs involved in sulfur metabolism from different habitats and host origins were constructed. Combined with principal coordinate analysis, it revealed that most AMGs associated with sulfur metabolism clustered separately, indicating the abundance and uniqueness in this region. The sulfur metabolism pathways involved by AMGs were mainly SOX systems, among which sulfur oxidation was associated with moaD and dsrE genes, while sulfur transport was related to tusA genes. It provides an insight into the biogeochemical sulfur cycling in plateau wetlands and lays the foundation for further study on the co-evolution of virus and host.

Diversity, distribution, and functional potentials of magroviruses from marine and brackish waters

Marine group II (MGII) archaea (Ca. Poseidoniales) are among the most abundant microbes in global oceanic surface waters and play an important role in driving marine biogeochemical cycles. Magroviruses - the viruses of MGII archaea have been recently found to occur ubiquitously in surface ocean. However, their diversity, distribution, and potential ecological functions in coastal zones especially brackish waters are unknown. Here we obtained 234 non-redundant magroviral genomes from brackish surface waters by using homology searches for viral signature proteins highlighting the uncovered vast diversity of this novel viral group. Phylogenetic analysis based on these brackish magroviruses along with previously reported marine ones identified six taxonomic groups with close evolutionary connection to both haloviruses and the viruses of Marine Group I archaea. Magroviruses were present abundantly both in brackish and open ocean samples with some showing habitat specification and others having broad spectrums of distribution between different habitats. Genome annotation suggests they may be involved in regulating multiple metabolic pathways of MGII archaea. Our results uncover the previously overlooked diversity and ecological potentials of a major archaeal virial group in global ocean and brackish waters and shed light on the cryptic evolutionary history of archaeal viruses.

Virus diversity and interactions with hosts in deep-sea hydrothermal vents

BACKGROUND

The deep sea harbors many viruses, yet their diversity and interactions with hosts in hydrothermal ecosystems are largely unknown. Here, we analyzed the viral composition, distribution, host preference, and metabolic potential in different habitats of global hydrothermal vents, including vent plumes, background seawater, diffuse fluids, and sediments.

RESULTS

From 34 samples collected at eight vent sites, a total of 4662 viral populations (vOTUs) were recovered from the metagenome assemblies, encompassing diverse phylogenetic groups and defining many novel lineages. Apart from the abundant unclassified viruses, tailed phages are most predominant across the global hydrothermal vents, while single-stranded DNA viruses, including Microviridae and small eukaryotic viruses, also constitute a significant part of the viromes. As revealed by protein-sharing network analysis, hydrothermal vent viruses formed many novel genus-level viral clusters and are highly endemic to specific vent sites and habitat types. Only 11% of the vOTUs can be linked to hosts, which are the key microbial taxa of hydrothermal habitats, such as Gammaproteobacteria and Campylobacterota. Intriguingly, vent viromes share some common metabolic features in that they encode auxiliary genes that are extensively involved in the metabolism of carbohydrates, amino acids, cofactors, and vitamins. Specifically, in plume viruses, various auxiliary genes related to methane, nitrogen, and sulfur metabolism were observed, indicating their contribution to host energy conservation. Moreover, the prevalence of sulfur-relay pathway genes indicated the significant role of vent viruses in stabilizing the tRNA structure, which promotes host adaptation to steep environmental gradients.

CONCLUSIONS

The deep-sea hydrothermal systems hold untapped viral diversity with novelty. They may affect both vent prokaryotic and eukaryotic communities and modulate host metabolism related to vent adaptability. More explorations are needed to depict global vent virus diversity and its roles in this unique ecosystem. Video Abstract.

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.

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.

Ecosystem-specific microbiota and microbiome databases in the era of big data

The rapid development of sequencing methods over the past decades has accelerated both the potential scope and depth of microbiota and microbiome studies. Recent developments in the field have been marked by an expansion away from purely categorical studies towards a greater investigation of community functionality. As in-depth genomic and environmental coverage is often distributed unequally across major taxa and ecosystems, it can be difficult to identify or substantiate relationships within microbial communities. Generic databases containing datasets from diverse ecosystems have opened a new era of data accessibility despite costs in terms of data quality and heterogeneity. This challenge is readily embodied in the integration of meta-omics data alongside habitat-specific standards which help contextualise datasets both in terms of sample processing and background within the ecosystem. A special case of large genomic repositories, ecosystem-specific databases (ES-DB's), have emerged to consolidate and better standardise sample processing and analysis protocols around individual ecosystems under study, allowing independent studies to produce comparable datasets. Here, we provide a comprehensive review of this emerging tool for microbial community analysis in relation to current trends in the field. We focus on the factors leading to the formation of ES-DB's, their comparison to traditional microbial databases, the potential for ES-DB integration with meta-omics platforms, as well as inherent limitations in the applicability of ES-DB's.

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.

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.

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.

Phage Genome Diversity in a Biogas-Producing Microbiome Analyzed by Illumina and Nanopore GridION Sequencing

The microbial biogas network is complex and intertwined, and therefore relatively stable in its overall functionality. However, if key functional groups of microorganisms are affected by biotic or abiotic factors, the entire efficacy may be impaired. Bacteriophages are hypothesized to alter the steering process of the microbial network. In this study, an enriched fraction of virus-like particles was extracted from a mesophilic biogas reactor and sequenced on the Illumina MiSeq and Nanopore GridION sequencing platforms. Metagenome data analysis resulted in identifying 375 metagenome-assembled viral genomes (MAVGs). Two-thirds of the classified sequences were only assigned to the superkingdom Viruses and the remaining third to the family Siphoviridae, followed by Myoviridae, Podoviridae, Tectiviridae, and Inoviridae. The metavirome showed a close relationship to the phage genomes that infect members of the classes Clostridia and Bacilli. Using publicly available biogas metagenomic data, a fragment recruitment approach showed the widespread distribution of the MAVGs studied in other biogas microbiomes. In particular, phage sequences from mesophilic microbiomes were highly similar to the phage sequences of this study. Accordingly, the virus particle enrichment approach and metavirome sequencing provided additional genome sequence information for novel virome members, thus expanding the current knowledge of viral genetic diversity in biogas reactors.

DNA Viral Diversity, Abundance, and Functional Potential Vary across Grassland Soils with a Range of Historical Moisture Regimes

Soil viruses are abundant, but the influence of the environment and climate on soil viruses remains poorly understood. Here, we addressed this gap by comparing the diversity, abundance, lifestyle, and metabolic potential of DNA viruses in three grassland soils with historical differences in average annual precipitation, low in eastern Washington (WA), high in Iowa (IA), and intermediate in Kansas (KS). Bioinformatics analyses were applied to identify a total of 2,631 viral contigs, including 14 complete viral genomes from three deep metagenomes (1 terabase [Tb] each) that were sequenced from bulk soil DNA. An additional three replicate metagenomes (∼0.5 Tb each) were obtained from each location for statistical comparisons. Identified viruses were primarily bacteriophages targeting dominant bacterial taxa. Both viral and host diversity were higher in soil with lower precipitation. Viral abundance was also significantly higher in the arid WA location than in IA and KS. More lysogenic markers and fewer clustered regularly interspaced short palindromic repeats (CRISPR) spacer hits were found in WA, reflecting more lysogeny in historically drier soil. More putative auxiliary metabolic genes (AMGs) were also detected in WA than in the historically wetter locations. The AMGs occurring in 18 pathways could potentially contribute to carbon metabolism and energy acquisition in their hosts. Structural equation modeling (SEM) suggested that historical precipitation influenced viral life cycle and selection of AMGs. The observed and predicted relationships between soil viruses and various biotic and abiotic variables have value for predicting viral responses to environmental change.

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

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

Microbial Ecology of Oxygen Minimum Zones Amidst Ocean Deoxygenation

Oxygen minimum zones (OMZs) have substantial effects on the global ecology and biogeochemical processes of marine microbes. However, the diversity and activity of OMZ microbes and their trophic interactions are only starting to be documented, especially in regard to the potential roles of viruses and protists. OMZs have expanded over the past 60 years and are predicted to expand due to anthropogenic climate change, furthering the need to understand these regions. This review summarizes the current knowledge of OMZ formation, the biotic and abiotic factors involved in OMZ expansion, and the microbial ecology of OMZs, emphasizing the importance of bacteria, archaea, viruses, and protists. We describe the recognized roles of OMZ microbes in carbon, nitrogen, and sulfur cycling, the potential of viruses in altering host metabolisms involved in these cycles, and the control of microbial populations by grazers and viruses. Further, we highlight the microbial community composition and roles of these organisms in oxic and anoxic depths within the water column and how these differences potentially inform how microbial communities will respond to deoxygenation. Additionally, the current literature on the alteration of microbial communities by other key climate change parameters such as temperature and pH are considered regarding how OMZ microbes might respond to these pressures. Finally, we discuss what knowledge gaps are present in understanding OMZ microbial communities and propose directions that will begin to close these gaps.

Imaging Techniques for Detecting Prokaryotic Viruses in Environmental Samples

Viruses are the most abundant biological entities on Earth with an estimate of 1031 viral particles across all ecosystems. Prokaryotic viruses-bacteriophages and archaeal viruses-influence global biogeochemical cycles by shaping microbial communities through predation, through the effect of horizontal gene transfer on the host genome evolution, and through manipulating the host cellular metabolism. Imaging techniques have played an important role in understanding the biology and lifestyle of prokaryotic viruses. Specifically, structure-resolving microscopy methods, for example, transmission electron microscopy, are commonly used for understanding viral morphology, ultrastructure, and host interaction. These methods have been applied mostly to cultivated phage-host pairs. However, recent advances in environmental genomics have demonstrated that the majority of viruses remain uncultivated, and thus microscopically uncharacterized. Although light- and structure-resolving microscopy of viruses from environmental samples is possible, quite often the link between the visualization and the genomic information of uncultivated prokaryotic viruses is missing. In this minireview, we summarize the current state of the art of imaging techniques available for characterizing viruses in environmental samples and discuss potential links between viral imaging and environmental genomics for shedding light on the morphology of uncultivated viruses and their lifestyles in Earth's ecosystems.

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

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

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.

Expanding Archaeal Diversity and Phylogeny: Past, Present, and Future

The discovery of the Archaea is a major scientific hallmark of the twentieth century. Since then, important features of their cell biology, physiology, ecology, and diversity have been revealed. Over the course of some 40 years, the diversity of known archaea has expanded from 2 to about 30 phyla comprising over 20,000 species. Most of this archaeal diversity has been revealed by environmental 16S rRNA amplicon sequencing surveys using a broad range of universal and targeted primers. Of the few primers that target a large fraction of known archaeal diversity, all display a bias against recently discovered lineages, which limits studies aiming to survey overall archaeal diversity. Induced by genomic exploration of archaeal diversity, and improved phylogenomics approaches, archaeal taxonomic classification has been frequently revised. Due to computational limitations and continued discovery of new lineages, a stable archaeal phylogeny is not yet within reach. Obtaining phylogenetic and taxonomic consensus of archaea should be a high priority for the archaeal research community. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Methane-derived carbon flows into host-virus networks at different trophic levels in soil

The concentration of atmospheric methane (CH₄) continues to increase with microbial communities controlling soil-atmosphere fluxes. While there is substantial knowledge of the diversity and function of prokaryotes regulating CH₄ production and consumption, their active interactions with viruses in soil have not been identified. Metagenomic sequencing of soil microbial communities enables identification of linkages between viruses and hosts. However, this does not determine if these represent current or historical interactions nor whether a virus or host are active. In this study, we identified active interactions between individual host and virus populations in situ by following the transfer of assimilated carbon. Using DNA stable-isotope probing combined with metagenomic analyses, we characterized CH₄-fueled microbial networks in acidic and neutral pH soils, specifically primary and secondary utilizers, together with the recent transfer of CH₄-derived carbon to viruses. A total of 63% of viral contigs from replicated soil incubations contained homologs of genes present in known methylotrophic bacteria. Genomic sequences of ¹³C-enriched viruses were represented in over one-third of spacers in CRISPR arrays of multiple closely related Methylocystis populations and revealed differences in their history of viral interaction. Viruses infecting nonmethanotrophic methylotrophs and heterotrophic predatory bacteria were also identified through the analysis of shared homologous genes, demonstrating that carbon is transferred to a diverse range of viruses associated with CH₄-fueled microbial food networks.

Hostile Takeover: How Viruses Reprogram Prokaryotic Metabolism

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

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.

Lytic archaeal viruses infect abundant primary producers in Earth's crust

The continental subsurface houses a major portion of life's abundance and diversity, yet little is known about viruses infecting microbes that reside there. Here, we use a combination of metagenomics and virus-targeted direct-geneFISH (virusFISH) to show that highly abundant carbon-fixing organisms of the uncultivated genus Candidatus Altiarchaeum are frequent targets of previously unrecognized viruses in the deep subsurface. Analysis of CRISPR spacer matches display resistances of Ca. Altiarchaea against eight predicted viral clades, which show genomic relatedness across continents but little similarity to previously identified viruses. Based on metagenomic information, we tag and image a putatively viral genome rich in protospacers using fluorescence microscopy. VirusFISH reveals a lytic lifestyle of the respective virus and challenges previous predictions that lysogeny prevails as the dominant viral lifestyle in the subsurface. CRISPR development over time and imaging of 18 samples from one subsurface ecosystem suggest a sophisticated interplay of viral diversification and adapting CRISPR-mediated resistances of Ca. Altiarchaeum. We conclude that infections of primary producers with lytic viruses followed by cell lysis potentially jump-start heterotrophic carbon cycling in these subsurface ecosystems.

Ecogenomics of Groundwater Phages Suggests Niche Differentiation Linked to Specific Environmental Tolerance

Viruses are ubiquitous microbiome components, shaping ecosystems via strain-specific predation, horizontal gene transfer and redistribution of nutrients through host lysis. Viral impacts are important in groundwater ecosystems, where microbes drive many nutrient fluxes and metabolic processes; however, little is known about the diversity of viruses in these environments. We analyzed four groundwater plasmidomes (the entire plasmid content of an environment) and identified 200 viral sequences, which clustered into 41 genus-level viral clusters (approximately equivalent to viral genera) including 9 known and 32 putative new genera. We used publicly available bacterial whole-genome sequences (WGS) and WGS from 261 bacterial isolates from this groundwater environment to identify potential viral hosts. We linked 76 of the 200 viral sequences to a range of bacterial phyla, the majority associated with Proteobacteria, followed by Firmicutes, Bacteroidetes, and Actinobacteria. The publicly available WGS enabled mapping bacterial hosts to several viral sequences. The WGS of groundwater isolates increased the depth of host prediction by allowing host identification at the strain level. The latter included 4 viruses that were almost entirely (>99% query coverage, >99% identity) identified as integrated in the genomes of Pseudomonas, Acidovorax, and Castellaniella strains, resulting in high-confidence host assignments. Lastly, 21 of these viruses carried putative auxiliary metabolite genes for metal and antibiotic resistance, which might drive their infection cycles and/or provide selective advantage to infected hosts. Exploring the groundwater virome provides a necessary foundation for integration of viruses into ecosystem models where they are key players in microbial adaption to environmental stress. IMPORTANCE To our knowledge, this is the first study to identify the bacteriophage distribution in a groundwater ecosystem shedding light on their prevalence and distribution across metal-contaminated and background sites. Our study is uniquely based on selective sequencing of solely the extrachromosomal elements of a microbiome followed by analysis for viral signatures, thus establishing a more focused approach for phage identifications. Using this method, we detected several novel phage genera along with those previously established. Our approach of using the whole-genome sequences of hundreds of bacterial isolates from the same site enabled us to make host assignments with high confidence, several at strain levels. Certain phage genes suggest that they provide an environment-specific selective advantage to their bacterial hosts. Our study lays the foundation for future research on directed phage isolations using specific bacterial host strains to further characterize groundwater phages, their life cycles, and their effects on groundwater microbiome and biogeochemistry.

RaFAH: Host prediction for viruses of Bacteria and Archaea based on protein content

Culture-independent approaches have recently shed light on the genomic diversity of viruses of prokaryotes. One fundamental question when trying to understand their ecological roles is: which host do they infect? To tackle this issue we developed a machine-learning approach named Random Forest Assignment of Hosts (RaFAH), that uses scores to 43,644 protein clusters to assign hosts to complete or fragmented genomes of viruses of Archaea and Bacteria. RaFAH displayed performance comparable with that of other methods for virus-host prediction in three different benchmarks encompassing viruses from RefSeq, single amplified genomes, and metagenomes. RaFAH was applied to assembled metagenomic datasets of uncultured viruses from eight different biomes of medical, biotechnological, and environmental relevance. Our analyses led to the identification of 537 sequences of archaeal viruses representing unknown lineages, whose genomes encode novel auxiliary metabolic genes, shedding light on how these viruses interfere with the host molecular machinery. RaFAH is available at https://sourceforge.net/projects/rafah/.

•

RaFAH was developed to predict the hosts of viruses of Bacteria and Archaea

•

RaFAH displayed comparable or superior performance to other host-prediction tools

•

RaFAH performed well across viromes from eight different ecosystems

•

RaFAH identified hundreds of genomic sequences as derived from viruses of Archaea

Viruses that infect Bacteria and Archaea are ubiquitous and extremely abundant. Recent advances have led to the discovery of many thousands of complete and partial genomes of these biological entities. Understanding the biology of these viruses and how they influence their ecosystems depends on knowing which hosts they infect. We developed a tool that uses data from complete or fragmented genomes to predict the hosts of viruses using a machine-learning approach. Our tool, RaFAH, displayed performance comparable with or superior to that of other host-prediction tools. In addition, it identified hundreds of sequences as derived from the genomes of viruses of Archaea, which are one of the least characterized fractions of the global virosphere.

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.

Ecology of inorganic sulfur auxiliary metabolism in widespread bacteriophages

Microbial sulfur metabolism contributes to biogeochemical cycling on global scales. Sulfur metabolizing microbes are infected by phages that can encode auxiliary metabolic genes (AMGs) to alter sulfur metabolism within host cells but remain poorly characterized. Here we identified 191 phages derived from twelve environments that encoded 227 AMGs for oxidation of sulfur and thiosulfate (dsrA, dsrC/tusE, soxC, soxD and soxYZ). Evidence for retention of AMGs during niche-differentiation of diverse phage populations provided evidence that auxiliary metabolism imparts measurable fitness benefits to phages with ramifications for ecosystem biogeochemistry. Gene abundance and expression profiles of AMGs suggested significant contributions by phages to sulfur and thiosulfate oxidation in freshwater lakes and oceans, and a sensitive response to changing sulfur concentrations in hydrothermal environments. Overall, our study provides fundamental insights on the distribution, diversity, and ecology of phage auxiliary metabolism associated with sulfur and reinforces the necessity of incorporating viral contributions into biogeochemical configurations.

Benthic prokaryotic microbial community assembly and biogeochemical potentials in E. coli - Stressed aquatic ecosystems during plant decomposition

Benthic microbes play a crucial role in maintaining the biogeochemical balance of aquatic ecosystems especially the material cycling during plant decomposition. However, those systems in agricultural area are always threatened by agricultural run-off containing a mass of typical pathogenic invader- Escherichia coli. It is therefore vital to clarify the turnover, assembly, and geochemical functions of the E. coli invaded benthic prokaryotic microbial community during plant decomposition. During the decaying process, the key filtering factors of benthic community assembly were NH₄⁺-N (P < 0.001), NO₂^--N (P < 0.01), and Organic-N (P < 0.05). The E. coli colonized significantly in sediments (P < 0.001) and drove the turnover of the bacterial community (P = 0.001), which enhanced archaeal dominance in the benthic microbial network. E. coli also triggered niche structural variations. The biomass (%) of benthic nutrient cycling genera including Dechloromonas, Pseudomonas, Bacteroides, Candidatus_Methanofastidiosum, and Desulfomicrobium (P < 0.05) was altered by E. coli stress. The structural equation model illustrated that E. coli critically affected the benthic microbial geochemical functions in multiple pathways (P < 0.05). Our results provide new insights into benthic prokaryotic microbial community assembly and nutrient cycling and management under pollution stress.