IMG/VR v.2.0: an integrated data management and analysis system for cultivated and environmental viral genomes

Genome analysis of vB_SupP_AX, a novel N4-like phage infecting Sulfitobacter

Sulfitobacter is a bacterium recognized for its production of AMP-independent sulfite oxidase, which is instrumental in the creation of sulfite biosensors. This capability underscores its ecological and economic relevance. In this study, we present a newly discovered phage, Sulfitobacter phage vB_SupP_AX, which was isolated from Maidao of Qingdao, China. The vB_SupP_AX genome is linear and double-stranded and measures 75,445 bp with a GC content of 49%. It encompasses four transfer RNA (tRNA) sequences and 79 open reading frames (ORFs), one of which is an auxiliary metabolic gene encoding thioredoxin. Consistent with other N4-like phages, vB_SupP_AX possesses three distinct RNA polymerases and is characterized by the presence of four tRNA molecules. Comparative genomic and phylogenetic analyses position vB_SupP_AX and three other viral genomes from the Integrated Microbial Genomes/Virus v4 database within the Rhodovirinae virus subfamily. The identification of vB_SupP_AX enhances our understanding of virus-host interactions within marine ecosystems.

Diversity and potential host-interactions of viruses inhabiting deep-sea seamount sediments

Seamounts are globally distributed across the oceans and form one of the major oceanic biomes. Here, we utilized combined analyses of bulk metagenome and virome to study viral communities in seamount sediments in the western Pacific Ocean. Phylogenetic analyses and the protein-sharing network demonstrate extensive diversity and previously unknown viral clades. Inference of virus-host linkages uncovers extensive interactions between viruses and dominant prokaryote lineages, and suggests that viruses play significant roles in carbon, sulfur, and nitrogen cycling by compensating or augmenting host metabolisms. Moreover, temperate viruses are predicted to be prevalent in seamount sediments, which tend to carry auxiliary metabolic genes for host survivability. Intriguingly, the geographical features of seamounts likely compromise the connectivity of viral communities and thus contribute to the high divergence of viral genetic spaces and populations across seamounts. Altogether, these findings provides knowledge essential for understanding the biogeography and ecological roles of viruses in globally widespread seamounts.

Phage-inclusive profiling of human gut microbiomes with Phanta

Due to technical limitations, most gut microbiome studies have focused on prokaryotes, overlooking viruses. Phanta, a virome-inclusive gut microbiome profiling tool, overcomes the limitations of assembly-based viral profiling methods by using customized k-mer-based classification tools and incorporating recently published catalogs of gut viral genomes. Phanta's optimizations consider the small genome size of viruses, sequence homology with prokaryotes and interactions with other gut microbes. Extensive testing of Phanta on simulated data demonstrates that it quickly and accurately quantifies prokaryotes and viruses. When applied to 245 fecal metagenomes from healthy adults, Phanta identifies ~200 viral species per sample, ~5× more than standard assembly-based methods. We observe a ~2:1 ratio between DNA viruses and bacteria, with higher interindividual variability of the gut virome compared to the gut bacteriome. In another cohort, we observe that Phanta performs equally well on bulk versus virus-enriched metagenomes, making it possible to study prokaryotes and viruses in a single experiment, with a single analysis.

Biogeographic patterns and drivers of soil viromes

Viruses are crucial in shaping soil microbial functions and ecosystems. However, studies on soil viromes have been limited in both spatial scale and biome coverage. Here we present a comprehensive synthesis of soil virome biogeographic patterns using the Global Soil Virome dataset (GSV) wherein we analysed 1,824 soil metagenomes worldwide, uncovering 80,750 partial genomes of DNA viruses, 96.7% of which are taxonomically unassigned. The biogeography of soil viral diversity and community structure varies across different biomes. Interestingly, the diversity of viruses does not align with microbial diversity and contrasts with it by showing low diversity in forest and shrubland soils. Soil texture and moisture conditions are further corroborated as key factors affecting diversity by our predicted soil viral diversity atlas, revealing higher diversity in humid and subhumid regions. In addition, the binomial degree distribution pattern suggests a random co-occurrence pattern of soil viruses. These findings are essential for elucidating soil viral ecology and for the comprehensive incorporation of viruses into soil ecosystem models.

COBRA improves the completeness and contiguity of viral genomes assembled from metagenomes

Viruses are often studied using metagenome-assembled sequences, but genome incompleteness hampers comprehensive and accurate analyses. Contig Overlap Based Re-Assembly (COBRA) resolves assembly breakpoints based on the de Bruijn graph and joins contigs. Here we benchmarked COBRA using ocean and soil viral datasets. COBRA accurately joined the assembled sequences and achieved notably higher genome accuracy than binning tools. From 231 published freshwater metagenomes, we obtained 7,334 bacteriophage clusters, ~83% of which represent new phage species. Notably, ~70% of these were circular, compared with 34% before COBRA analyses. We expanded sampling of huge phages (≥200 kbp), the largest of which was curated to completion (717 kbp). Improved phage genomes from Rotsee Lake provided context for metatranscriptomic data and indicated the in situ activity of huge phages, whiB-encoding phages and cysC- and cysH-encoding phages. COBRA improves viral genome assembly contiguity and completeness, thus the accuracy and reliability of analyses of gene content, diversity and evolution.

Discovering and exploring the hidden diversity of human gut viruses using highly enriched virome samples

Viruses are an abundant and crucial component of the human microbiome, but accurately discovering them via metagenomics is still challenging. Currently, the available viral reference genomes poorly represent the diversity in microbiome samples, and expanding such a set of viral references is difficult. As a result, many viruses are still undetectable through metagenomics even when considering the power of de novo metagenomic assembly and binning, as viruses lack universal markers. Here, we describe a novel approach to catalog new viral members of the human gut microbiome and show how the resulting resource improves metagenomic analyses. We retrieved >3,000 viral-like particles (VLP) enriched metagenomic samples (viromes), evaluated the efficiency of the enrichment in each sample to leverage the viromes of highest purity, and applied multiple analysis steps involving assembly and comparison with hundreds of thousands of metagenome-assembled genomes to discover new viral genomes. We reported over 162,000 viral sequences passing quality control from thousands of gut metagenomes and viromes. The great majority of the retrieved viral sequences (~94.4%) were of unknown origin, most had a CRISPR spacer matching host bacteria, and four of them could be detected in >50% of a set of 18,756 gut metagenomes we surveyed. We included the obtained collection of sequences in a new MetaPhlAn 4.1 release, which can quantify reads within a metagenome matching the known and newly uncovered viral diversity. Additionally, we released the viral database for further virome and metagenomic studies of the human microbiome.

Large language models improve annotation of prokaryotic viral proteins

Viral genomes are poorly annotated in metagenomic samples, representing an obstacle to understanding viral diversity and function. Current annotation approaches rely on alignment-based sequence homology methods, which are limited by the paucity of characterized viral proteins and divergence among viral sequences. Here we show that protein language models can capture prokaryotic viral protein function, enabling new portions of viral sequence space to be assigned biologically meaningful labels. When applied to global ocean virome data, our classifier expanded the annotated fraction of viral protein families by 29%. Among previously unannotated sequences, we highlight the identification of an integrase defining a mobile element in marine picocyanobacteria and a capsid protein that anchors globally widespread viral elements. Furthermore, improved high-level functional annotation provides a means to characterize similarities in genomic organization among diverse viral sequences. Protein language models thus enhance remote homology detection of viral proteins, serving as a useful complement to existing approaches.

Genomic insights into the adaptation of Synechococcus to the coastal environment on Xiamen

Synechococcus are widely distributed in the global ocean, from the pelagic zone to coastal waters. However, little is known about Synechococcus in coastal seawater due to limitations in isolation and culture conditions. In this study, a combination of metagenomic sequencing technology, flow cytometry sorting, and multiple displacement amplification was used to investigate Synechococcus in the coastal seawater of Xiamen Island. The results revealed 18 clades of Synechococcus and diverse metabolic genes that appear to contribute to their adaptation to the coastal environment. Intriguingly, some metabolic genes related to the metabolism of carbohydrates, energy, nucleotides, and amino acids were found in 89 prophage regions that were detected in 16,258 Synechococcus sequences. The detected metabolic genes can enable prophages to contribute to the adaptation of Synechococcus hosts to the coastal marine environment. The detection of prophages also suggests that the cyanophages have infected Synechococcus. On the other hand, the identification of 773 genes associated with antiviral defense systems indicates that Synechococcus in Xiamen have evolved genetic traits in response to cyanophage infection. Studying the community diversity and functional genes of Synechococcus provides insights into their role in environmental adaptation and marine ecosystems.

Characterization and genomic analysis of a novel lytic phage vB_PstM_ZRG1 infecting Stutzerimonas stutzeri, representing a new viral genus, Elithevirus

Stutzerimonas stutzeri is an opportunistic pathogen widely distributed in the environment and displays diverse metabolic capabilities. In this study, a novel lytic S. stutzeri phage, named vB_PstM_ZRG1, was isolated from the seawater in the East China Sea (29°09'N, 123°39'E). vB_PstM_ZRG1 was stable at temperatures ranging from -20°C to 65°C and across a wide range of pH values from 3 to 10. The genome of vB_PstM_ZRG1 was determined to be a double-stranded DNA with a genome size of 52,767 bp, containing 78 putative open reading frames (ORFs). Three auxiliary metabolic genes encoded by phage vB_PstM_ZRG1 were predicted, including Toll/interleukin-1 receptor (TIR) domain, proline-alanine-alanine-arginine (PAAR) protein and SGNH (Ser-Gly-Asn-His) family hydrolase, especially TIR domain is not common in isolated phages. Phylogenic and network analysis showed that vB_PstM_ZRG1 has low similarity to other phage genomes in the GenBank and IMG/VR database, and might represent a novel viral genus, named Elithevirus. Additionally, the distribution map results indicated that vB_PstM_ZRG1 could infect both extreme colds- and warm-type hosts in the marine environment. In summary, our finding provided basic information for further research on the relationship between S. stutzeri and their phages, and expanded our understanding of genomic characteristics, phylogenetic diversity and distribution of Elithevirus.

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.

Unexplored diversity and ecological functions of transposable phages

Phages are prevalent in diverse environments and play major ecological roles attributed to their tremendous diversity and abundance. Among these viruses, transposable phages (TBPs) are exceptional in terms of their unique lifestyle, especially their replicative transposition. Although several TBPs have been isolated and the life cycle of the representative phage Mu has been extensively studied, the diversity distribution and ecological functions of TBPs on the global scale remain unknown. Here, by mining TBPs from enormous microbial genomes and viromes, we established a TBP genome dataset (TBPGD), that expands the number of accessible TBP genomes 384-fold. TBPs are prevalent in diverse biomes and show great genetic diversity. Based on taxonomic evaluations, we propose the categorization of TBPs into four viral groups, including 11 candidate subfamilies. TBPs infect multiple bacterial phyla, and seem to infect a wider range of hosts than non-TBPs. Diverse auxiliary metabolic genes (AMGs) are identified in the TBP genomes, and genes related to glycoside hydrolases and pyrimidine deoxyribonucleotide biosynthesis are highly enriched. Finally, the influences of TBPs on their hosts are experimentally examined by using the marine bacterium Shewanella psychrophila WP2 and its infecting transposable phage SP2. Collectively, our findings greatly expand the genetic diversity of TBPs, and comprehensively reveal their potential influences in various ecosystems.

DETIRE: a hybrid deep learning model for identifying viral sequences from metagenomes

A metagenome contains all DNA sequences from an environmental sample, including viruses, bacteria, archaea, and eukaryotes. Since viruses are of huge abundance and have caused vast mortality and morbidity to human society in history as a type of major pathogens, detecting viruses from metagenomes plays a crucial role in analyzing the viral component of samples and is the very first step for clinical diagnosis. However, detecting viral fragments directly from the metagenomes is still a tough issue because of the existence of a huge number of short sequences. In this study a hybrid Deep lEarning model for idenTifying vIral sequences fRom mEtagenomes (DETIRE) is proposed to solve the problem. First, the graph-based nucleotide sequence embedding strategy is utilized to enrich the expression of DNA sequences by training an embedding matrix. Then, the spatial and sequential features are extracted by trained CNN and BiLSTM networks, respectively, to enrich the features of short sequences. Finally, the two sets of features are weighted combined for the final decision. Trained by 220,000 sequences of 500 bp subsampled from the Virus and Host RefSeq genomes, DETIRE identifies more short viral sequences (<1,000 bp) than the three latest methods, such as DeepVirFinder, PPR-Meta, and CHEER. DETIRE is freely available at Github (https://github.com/crazyinter/DETIRE).

A First Insight into the Microbial and Viral Communities of Comau Fjord—A Unique Human-Impacted Ecosystem in Patagonia (42∘ S)

While progress has been made in surveying the oceans to understand microbial and viral communities, the coastal ocean and, specifically, estuarine waters, where the effects of anthropogenic activity are greatest, remain partially understudied. The coastal waters of Northern Patagonia are of interest since this region experiences high-density salmon farming as well as other disturbances such as maritime transport of humans and cargo. Here, we hypothesized that viral and microbial communities from the Comau Fjord would be distinct from those collected in global surveys yet would have the distinctive features of microbes from coastal and temperate regions. We further hypothesized that microbial communities will be functionally enriched in antibiotic resistance genes (ARGs) in general and in those related to salmon farming in particular. Here, the analysis of metagenomes and viromes obtained for three surface water sites showed that the structure of the microbial communities was distinct in comparison to global surveys such as the Tara Ocean, though their composition converges with that of cosmopolitan marine microbes belonging to Proteobacteria, Bacteroidetes, and Actinobacteria. Similarly, viral communities were also divergent in structure and composition but matched known viral members from North America and the southern oceans. Microbial communities were functionally enriched in ARGs dominated by beta-lactams and tetracyclines, bacitracin, and the group macrolide–lincosamide–streptogramin (MLS) but were not different from other communities from the South Atlantic, South Pacific, and Southern Oceans. Similarly, viral communities were characterized by exhibiting protein clusters similar to those described globally (Tara Oceans Virome); however, Comau Fjord viromes displayed up to 50% uniqueness in their protein content. Altogether, our results indicate that microbial and viral communities from the Comau Fjord are a reservoir of untapped diversity and that, given the increasing anthropogenic impacts in the region, they warrant further study, specifically regarding resilience and resistance against antimicrobials and hydrocarbons.

Metagenomic Sequencing Identified Specific Bacteriophage Signature Discriminating between Healthy and Diarrheal Neonatal Piglets

Neonatal diarrhea is one of the most severe diseases in human beings and pigs, leading to high mortality and growth faltering. Gut microbiome-related studies mostly focus on the relationship between bacteria and neonatal diarrhea onset, and no research study has investigated the role of the gut virome in neonatal diarrhea. Here, using metagenomic sequencing, we characterized the fecal viral community of diarrheal and healthy neonatal piglets. We found that the viral community of diarrheal piglets showed higher individual heterogeneity and elevated abundance of Myoviridae. By predicting the bacterial host of the identified viral genomes, phages infecting Proteobacteria, especially E. coli, were the dominant taxa in neonatal diarrheal piglets. Consistent with this, the antibiotic resistance gene of E. coli origin was also enriched in neonatal diarrheal piglets. Finally, we established a random forest model to accurately discriminate between neonatal diarrheal piglets and healthy controls and identified genus E. coli- and genus listeria-infecting bacteriophages, including psa and C5 viruses, as key biomarkers. In conclusion, we provide the first glance of viral community and function characteristics in diarrheal and healthy neonatal piglets. These findings expand our understanding of the relationship among phages, bacteria and diarrhea, and may facilitate the development of therapeutics for the prevention and treatment of neonatal diarrhea.

Viruses Regulate Microbial Community Assembly Together with Environmental Factors in Acid Mine Drainage

Viruses are widespread in various ecosystems, and they play important roles in regulating the microbial community via host-virus interactions. Recently, metagenomic studies showed that there are extremely diverse viruses in different environments from the ocean to the human gut, but the influences of viral communities on microbial communities are poorly understood, especially in extreme environments. Here, we used metagenomics to characterize microbial communities and viral communities in acid mine drainage (AMD) and evaluated how viruses shape microbial community constrained by the harsh environments. Our results showed that AMD viral communities are significantly associated with the microbial communities, and viral diversity has positive correlations with microbial diversity. Viral community explained more variations of microbial community composition than environmental factors in AMD of a polymetallic mine. Moreover, we found that viruses harboring adaptive genes regulate a relative abundance of hosts under the modulation of environmental factors, such as pH. We also observed that viral diversity has significant correlations with the global properties of microbial cooccurrence networks, such as modularity. In addition, the results of null modeling analyses revealed that viruses significantly affect microbial community phylogeny and play important roles in regulating ecological processes of community assembly, such as dispersal limitation and homogenous dispersal. Together, these results revealed that AMD viruses are critical forces driving microbial network and community assembly via host-virus interactions. IMPORTANCE Viruses as mobile genetic elements play critical roles in the adaptive evolution of their hosts in extreme environments. However, how viruses further influence microbial community structure and assembly is still unclear. A recent metagenomic study observed diverse viruses unexplored in acid mine drainage, revealing the associations between the viral community and environmental factors. Here, we showed that viruses together with environmental factors can constrain the relative abundance of host and microbial community assembly in AMD of copper mines and polymetallic mines. Our results highlight the importance of viruses in shaping the microbial community from the individual host level to the community level.

Deep metagenomic mining reveals bacteriophage sequence motifs driving host specificity

Bacteriophages can adapt to new hosts by altering sequence motifs through recombination or convergent evolution. Where such motifs exist and what fitness advantage they confer remains largely unknown. We report a new method, Bacteriophage Library Informed Sequence Scoring (BLISS), to discover sequence motifs in metagenomic datasets governing phage activity. BLISS uses experimental deep mutational scanning data to create sequence profiles to enable deep mining of metagenomes for functional motifs which are otherwise invisible to searches. We experimentally tested 10,073 BLISS-derived sequence motifs for the receptor-binding protein of the T7 phage. The screen revealed hundreds of T7 variants with novel host specificity with functional motifs sourced from distant families besides other major phyla. Position, substitution and location preferences on T7 dictated different specificities. To demonstrate therapeutic utility, we engineered highly active T7 variants against urinary tract pathogens. BLISS is a powerful tool to unlock the functional potential encoded in phage metagenomes.

IMG/VR v4: an expanded database of uncultivated virus genomes within a framework of extensive functional, taxonomic, and ecological metadata

Viruses are widely recognized as critical members of all microbiomes. Metagenomics enables large-scale exploration of the global virosphere, progressively revealing the extensive genomic diversity of viruses on Earth and highlighting the myriad of ways by which viruses impact biological processes. IMG/VR provides access to the largest collection of viral sequences obtained from (meta)genomes, along with functional annotation and rich metadata. A web interface enables users to efficiently browse and search viruses based on genome features and/or sequence similarity. Here, we present the fourth version of IMG/VR, composed of >15 million virus genomes and genome fragments, a ≈6-fold increase in size compared to the previous version. These clustered into 8.7 million viral operational taxonomic units, including 231 408 with at least one high-quality representative. Viral sequences in IMG/VR are now systematically identified from genomes, metagenomes, and metatranscriptomes using a new detection approach (geNomad), and IMG standard annotation are complemented with genome quality estimation using CheckV, taxonomic classification reflecting the latest taxonomic standards, and microbial host taxonomy prediction. IMG/VR v4 is available at https://img.jgi.doe.gov/vr, and the underlying data are available to download at https://genome.jgi.doe.gov/portal/IMG_VR.

Marine DNA methylation patterns are associated with microbial community composition and inform virus-host dynamics

BACKGROUND

DNA methylation in prokaryotes is involved in many different cellular processes including cell cycle regulation and defense against viruses. To date, most prokaryotic methylation systems have been studied in culturable microorganisms, resulting in a limited understanding of DNA methylation from a microbial ecology perspective. Here, we analyze the distribution patterns of several microbial epigenetics marks in the ocean microbiome through genome-centric metagenomics across all domains of life.

RESULTS

We reconstructed 15,056 viral, 252 prokaryotic, 56 giant viral, and 6 eukaryotic metagenome-assembled genomes from northwest Pacific Ocean seawater samples using short- and long-read sequencing approaches. These metagenome-derived genomes mostly represented novel taxa, and recruited a majority of reads. Thanks to single-molecule real-time (SMRT) sequencing technology, base modification could also be detected for these genomes. This showed that DNA methylation can readily be detected across dominant oceanic bacterial, archaeal, and viral populations, and microbial epigenetic changes correlate with population differentiation. Furthermore, our genome-wide epigenetic analysis of Pelagibacter suggests that GANTC, a DNA methyltransferase target motif, is related to the cell cycle and is affected by environmental conditions. Yet, the presence of this motif also partitions the phylogeny of the Pelagibacter phages, possibly hinting at a competitive co-evolutionary history and multiple effects of a single methylation mark.

CONCLUSIONS

Overall, this study elucidates that DNA methylation patterns are associated with ecological changes and virus-host dynamics in the ocean microbiome. Video Abstract.

Insights into the global freshwater virome

Viruses are by far the most abundant life forms on this planet. Yet, the full viral diversity remains mostly unknown, especially in environments like freshwater. Therefore, we aimed to study freshwater viruses in a global context. To this end, we downloaded 380 publicly available viral metagenomes (>1 TB). More than 60% of these metagenomes were discarded based on their levels of cellular contamination assessed by ribosomal DNA content. For the remaining metagenomes, assembled contigs were decontaminated using two consecutive steps, eventually yielding 273,365 viral contigs longer than 1,000 bp. Long enough contigs (≥ 10 kb) were clustered to identify novel genomes/genome fragments. We could recover 549 complete circular and high-quality draft genomes, out of which 10 were recognized as being novel. Functional annotation of these genomes showed that most of the annotated coding sequences are DNA metabolic genes or phage structural genes. On the other hand, taxonomic analysis of viral contigs showed that most of the assigned contigs belonged to the order Caudovirales, particularly the families of Siphoviridae, Myoviridae, and Podoviridae. The recovered viral contigs contained several auxiliary metabolic genes belonging to several metabolic pathways, especially carbohydrate and amino acid metabolism in addition to photosynthesis as well as hydrocarbon degradation and antibiotic resistance. Overall, we present here a set of prudently chosen viral contigs, which should not only help better understanding of freshwater viruses but also be a valuable resource for future virome studies.

Extensive gut virome variation and its associations with host and environmental factors in a population-level cohort

Indigenous bacteriophage communities (virome) in the human gut have a huge impact on the structure and function of gut bacterial communities (bacteriome), but virome variation at a population scale is not fully investigated yet. Here, we analyse the gut dsDNA virome in the Japanese 4D cohort of 4198 deeply phenotyped individuals. By assembling metagenomic reads, we discover thousands of high-quality phage genomes including previously uncharacterised phage clades with different bacterial hosts than known major ones. The distribution of host bacteria is a strong determinant for the distribution of phages in the gut, and virome diversity is highly correlated with anti-viral defence mechanisms of the bacteriome, such as CRISPR-Cas and restriction-modification systems. We identify 97 various intrinsic/extrinsic factors that significantly affect the virome structure, including age, sex, lifestyle, and diet, most of which showed consistent associations with both phages and their predicted bacterial hosts. Among the metadata categories, disease and medication have the strongest effects on the virome structure. Overall, these results present a basis to understand the symbiotic communities of bacteria and their viruses in the human gut, which will facilitate the medical and industrial applications of indigenous viruses.

Here, Nishijima et al. perform a large-scale analysis of the human gut virome in the Japanese 4D cohort of 4198 deeply phenotyped individuals, revealing thousands of bacteriophage genomes, virus-bacteria interactions, and describing associations with various host and environmental factors.

Viral informatics: bioinformatics-based solution for managing viral infections

Several new viral infections have emerged in the human population and establishing as global pandemics. With advancements in translation research, the scientific community has developed potential therapeutics to eradicate or control certain viral infections, such as smallpox and polio, responsible for billions of disabilities and deaths in the past. Unfortunately, some viral infections, such as dengue virus (DENV) and human immunodeficiency virus-1 (HIV-1), are still prevailing due to a lack of specific therapeutics, while new pathogenic viral strains or variants are emerging because of high genetic recombination or cross-species transmission. Consequently, to combat the emerging viral infections, bioinformatics-based potential strategies have been developed for viral characterization and developing new effective therapeutics for their eradication or management. This review attempts to provide a single platform for the available wide range of bioinformatics-based approaches, including bioinformatics methods for the identification and management of emerging or evolved viral strains, genome analysis concerning the pathogenicity and epidemiological analysis, computational methods for designing the viral therapeutics, and consolidated information in the form of databases against the known pathogenic viruses. This enriched review of the generally applicable viral informatics approaches aims to provide an overview of available resources capable of carrying out the desired task and may be utilized to expand additional strategies to improve the quality of translation viral informatics research.

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

Insights into the human gut virome by sampling a population from the Indian subcontinent

Gut virome plays an important role in human physiology but remains poorly understood. This study reports an investigation of the human gut DNA-virome of a previously unexplored ethnic population through metagenomics of faecal samples collected from individuals residing in Northern India. Analysis shows that, similar to the populations investigated earlier, majority of the identified virome belongs to bacteriophages and a smaller fraction (<20 %) consists of viruses that infect animals, archaea, protists, multiple domains or plants. However, crAss-like phages, in this population, are dominated by the genera VI, VII and VIII. Interestingly, it also reveals the presence of a virus family, Sphaerolipoviridae, which has not been detected in the human gut earlier. Viral families, Siphoviridae, Myoviridae, Podoviridae, Microviridae, Herelleviridae and Phycodnaviridae are detected in all of the analysed individuals, which supports the existence of a core virome. Lysogeny-associated genes were found in less than 10 % of the assembled genomes and a negative correlation was observed in the richness of bacterial and free-viral species, suggesting that the dominant lifestyle of gut phage is not lysogenic. This is in contrast to some of the earlier studies. Further, several hundred high-quality viral genomes were recovered. Detailed characterization of these genomes would be useful for understanding the biology of these viruses and their significance in human physiology.

Revolutionized virome research using systems microbiology approaches

Currently, both pathogenic and commensal viruses are continuously being discovered and acknowledged as ubiquitous components of microbial communities. The advancements of systems microbiological approaches have changed the face of virome research. Here, we focus on viral metagenomic approach to study virus community and their interactions with other microbial members as well as their hosts. This review also summarizes challenges, limitations, and benefits of the current virome approaches. Potentially, the studies of virome can be further applied in various biological and clinical fields.

Metagenomic Analyses of Multiple Gut Datasets Revealed the Association of Phage Signatures in Colorectal Cancer

The association of colorectal cancer (CRC) and the human gut microbiome dysbiosis has been the focus of several studies in the past. Many bacterial taxa have been shown to have differential abundance among CRC patients compared to healthy controls. However, the relationship between CRC and non-bacterial gut microbiome such as the gut virome is under-studied and not well understood. In this study we conducted a comprehensive analysis of the association of viral abundances with CRC using metagenomic shotgun sequencing data of 462 CRC subjects and 449 healthy controls from 7 studies performed in 8 different countries. Despite the high heterogeneity, our results showed that the virome alpha diversity was consistently higher in CRC patients than in healthy controls (p-value &lt;0.001). This finding is in sharp contrast to previous reports of low alpha diversity of prokaryotes in CRC compared to healthy controls. In addition to the previously known association of Podoviridae, Siphoviridae and Myoviridae with CRC, we further demonstrate that Herelleviridae, a newly constructed viral family, is significantly depleted in CRC subjects. Our interkingdom association analysis reveals a less intertwined correlation between the gut virome and bacteriome in CRC compared to healthy controls. Furthermore, we show that the viral abundance profiles can be used to accurately predict CRC disease status (AUROC &gt;0.8) in both within-study and cross-study settings. The combination of training sets resulted in rather generalized and accurate prediction models. Our study clearly shows that subjects with colorectal cancer harbor a distinct human gut virome profile which may have an important role in this disease.

The pediatric virome in health and disease

Associations between the global microbiome and diseases of children have been studied extensively; however, research on the viral component of the microbiome, the "virome," is less advanced. The analysis of disease associations with the virome is often technically challenging, requiring a close examination of the "virome dark matter." The gut is a particularly rich source of viral particles, and now multiple studies have reported intriguing associations of the virome with childhood diseases. For example, virome studies have elucidated new lineages of gut viruses that appear to be tightly associated with childhood diarrhea, and consistent patterns are starting to emerge from virome studies in pediatric IBD. In this review, we summarize the methods for studying the virome and recent research on the nature of the virome during childhood, focusing on specific studies of the intestinal virome in pediatric diseases.

Host-Associated Phages Disperse across the Extraterrestrial Analogue Antarctica

Extreme Antarctic conditions provide one of the closest analogues of extraterrestrial environments. Since air and snow samples, especially from polar regions, yield DNA amounts in the lower picogram range, binning of prokaryotic genomes is challenging and renders studying the dispersal of biological entities across these environments difficult. Here, we hypothesized that dispersal of host-associated bacteriophages (adsorbed, replicating, or prophages) across the Antarctic continent can be tracked via their genetic signatures, aiding our understanding of virus and host dispersal across long distances. Phage genome fragments (PGFs) reconstructed from surface snow metagenomes of three Antarctic stations were assigned to four host genomes, mainly Betaproteobacteria, including Ralstonia spp. We reconstructed the complete genome of a temperate phage with nearly complete alignment to a prophage in the reference genome of Ralstonia pickettii 12D. PGFs from different stations were related to each other at the genus level and matched similar hosts. Metagenomic read mapping and nucleotide polymorphism analysis revealed a wide dispersal of highly identical PGFs, 13 of which were detected in seawater from the Western Antarctic Peninsula at a distance of 5,338 km from the snow sampling stations. Our results suggest that host-associated phages, especially of Ralstonia sp., disperse over long distances despite the harsh conditions of the Antarctic continent. Given that 14 phages associated with two R. pickettii draft genomes isolated from space equipment were identified, we conclude that Ralstonia phages are ideal mobile genetic elements to track dispersal and contamination in ecosystems relevant for astrobiology. IMPORTANCE Host-associated phages of the bacterium Ralstonia identified in snow samples can be used to track microbial dispersal over thousands of kilometers across the Antarctic continent, which functions as an extraterrestrial analogue because of its harsh environmental conditions. Due to the presence of these bacteria carrying genome-integrated prophages on space-related equipment and the potential for dispersal of host-associated phages demonstrated here, our work has implications for planetary protection, a discipline in astrobiology interested in preventing contamination of celestial bodies with alien biomolecules or forms of life.

Quantifying and Cataloguing Unknown Sequences within Human Microbiomes

Advances in genome sequencing technologies and lower costs have enabled the exploration of a multitude of known and novel environments and microbiomes. This has led to an exponential growth in the raw sequence data that are deposited in online repositories. Metagenomic and metatranscriptomic data sets are typically analysed with regard to a specific biological question. However, it is widely acknowledged that these data sets are comprised of a proportion of sequences that bear no similarity to any currently known biological sequence, and this so-called "dark matter" is often excluded from downstream analyses. In this study, a systematic framework was developed to assemble, identify, and measure the proportion of unknown sequences present in distinct human microbiomes. This framework was applied to 40 distinct studies, comprising 963 samples, and covering 10 different human microbiomes including fecal, oral, lung, skin, and circulatory system microbiomes. We found that while the human microbiome is one of the most extensively studied, on average 2% of assembled sequences have not yet been taxonomically defined. However, this proportion varied extensively among different microbiomes and was as high as 25% for skin and oral microbiomes that have more interactions with the environment. A rate of taxonomic characterization of 1.64% of unknown sequences being characterized per month was calculated from these taxonomically unknown sequences discovered in this study. A cross-study comparison led to the identification of similar unknown sequences in different samples and/or microbiomes. Both our computational framework and the novel unknown sequences produced are publicly available for future cross-referencing. Our approach led to the discovery of several novel viral genomes that bear no similarity to sequences in the public databases. Some of these are widespread as they have been found in different microbiomes and studies. Hence, our study illustrates how the systematic characterization of unknown sequences can help the discovery of novel microbes, and we call on the research community to systematically collate and share the unknown sequences from metagenomic studies to increase the rate at which the unknown sequence space can be classified.

Phage-encoded ribosomal protein S21 expression is linked to late-stage phage replication

The ribosomal protein S21 (bS21) gene has been detected in diverse viruses with a large range of genome sizes, yet its in situ expression and potential significance have not been investigated. Here, we report five closely related clades of bacteriophages (phages) represented by 47 genomes (8 curated to completion and up to 331 kbp in length) that encode a bS21 gene. The bS21 gene is on the reverse strand within a conserved region that encodes the large terminase, major capsid protein, prohead protease, portal vertex proteins, and some hypothetical proteins. Based on CRISPR spacer targeting, the predominance of bacterial taxonomic affiliations of phage genes with those from Bacteroidetes, and the high sequence similarity of the phage bS21 genes and those from Bacteroidetes classes of Flavobacteriia, Cytophagia and Saprospiria, these phages are predicted to infect diverse Bacteroidetes species that inhabit a range of depths in freshwater lakes. Thus, bS21 phages have the potential to impact microbial community composition and carbon turnover in lake ecosystems. The transcriptionally active bS21-encoding phages were likely in the late stage of replication when collected, as core structural genes and bS21 were highly expressed. Thus, our analyses suggest that the phage bS21, which is involved in translation initiation, substitutes into the Bacteroidetes ribosomes and selects preferentially for phage transcripts during the late-stage replication when large-scale phage protein production is required for assembly of phage particles.

Highly diverse flavobacterial phages isolated from North Sea spring blooms

It is generally recognized that phages are a mortality factor for their bacterial hosts. This could be particularly true in spring phytoplankton blooms, which are known to be closely followed by a highly specialized bacterial community. We hypothesized that phages modulate these dense heterotrophic bacteria successions following phytoplankton blooms. In this study, we focused on Flavobacteriia, because they are main responders during these blooms and have an important role in the degradation of polysaccharides. A cultivation-based approach was used, obtaining 44 lytic flavobacterial phages (flavophages), representing twelve new species from two viral realms. Taxonomic analysis allowed us to delineate ten new phage genera and ten new families, from which nine and four, respectively, had no previously cultivated representatives. Genomic analysis predicted various life styles and genomic replication strategies. A likely eukaryote-associated host habitat was reflected in the gene content of some of the flavophages. Detection in cellular metagenomes and by direct-plating showed that part of these phages were actively replicating in the environment during the 2018 spring bloom. Furthermore, CRISPR/Cas spacers and re-isolation during two consecutive years suggested that, at least part of the new flavophages are stable components of the microbial community in the North Sea. Together, our results indicate that these diverse flavophages have the potential to modulate their respective host populations.

Remarkably coherent population structure for a dominant Antarctic Chlorobium species

BACKGROUND

In Antarctica, summer sunlight enables phototrophic microorganisms to drive primary production, thereby "feeding" ecosystems to enable their persistence through the long, dark winter months. In Ace Lake, a stratified marine-derived system in the Vestfold Hills of East Antarctica, a Chlorobium species of green sulphur bacteria (GSB) is the dominant phototroph, although its seasonal abundance changes more than 100-fold. Here, we analysed 413 Gb of Antarctic metagenome data including 59 Chlorobium metagenome-assembled genomes (MAGs) from Ace Lake and nearby stratified marine basins to determine how genome variation and population structure across a 7-year period impacted ecosystem function.

RESULTS

A single species, Candidatus Chlorobium antarcticum (most similar to Chlorobium phaeovibrioides DSM265) prevails in all three aquatic systems and harbours very little genomic variation (≥ 99% average nucleotide identity). A notable feature of variation that did exist related to the genomic capacity to biosynthesize cobalamin. The abundance of phylotypes with this capacity changed seasonally ~ 2-fold, consistent with the population balancing the value of a bolstered photosynthetic capacity in summer against an energetic cost in winter. The very high GSB concentration (> 10⁸ cells ml^-1 in Ace Lake) and seasonal cycle of cell lysis likely make Ca. Chlorobium antarcticum a major provider of cobalamin to the food web. Analysis of Ca. Chlorobium antarcticum viruses revealed the species to be infected by generalist (rather than specialist) viruses with a broad host range (e.g., infecting Gammaproteobacteria) that were present in diverse Antarctic lakes. The marked seasonal decrease in Ca. Chlorobium antarcticum abundance may restrict specialist viruses from establishing effective lifecycles, whereas generalist viruses may augment their proliferation using other hosts.

CONCLUSION

The factors shaping Antarctic microbial communities are gradually being defined. In addition to the cold, the annual variation in sunlight hours dictates which phototrophic species can grow and the extent to which they contribute to ecosystem processes. The Chlorobium population studied was inferred to provide cobalamin, in addition to carbon, nitrogen, hydrogen, and sulphur cycling, as critical ecosystem services. The specific Antarctic environmental factors and major ecosystem benefits afforded by this GSB likely explain why such a coherent population structure has developed in this Chlorobium species. Video abstract.

Presence and Persistence of Putative Lytic and Temperate Bacteriophages in Vaginal Metagenomes from South African Adolescents

The interaction between gut bacterial and viral microbiota is thought to be important in human health. While fluctuations in female genital tract (FGT) bacterial microbiota similarly determine sexual health, little is known about the presence, persistence, and function of vaginal bacteriophages. We conducted shotgun metagenome sequencing of cervicovaginal samples from South African adolescents collected longitudinally, who received no antibiotics. We annotated viral reads and circular bacteriophages, identified CRISPR loci and putative prophages, and assessed their diversity, persistence, and associations with bacterial microbiota composition. Siphoviridae was the most prevalent bacteriophage family, followed by Myoviridae, Podoviridae, Herelleviridae, and Inoviridae. Full-length siphoviruses targeting bacterial vaginosis (BV)-associated bacteria were identified, suggesting their presence in vivo. CRISPR loci and prophage-like elements were common, and genomic analysis suggested higher diversity among Gardnerella than Lactobacillus prophages. We found that some prophages were highly persistent within participants, and identical prophages were present in cervicovaginal secretions of multiple participants, suggesting that prophages, and thus bacterial strains, are shared between adolescents. The number of CRISPR loci and prophages were associated with vaginal microbiota stability and absence of BV. Our analysis suggests that (pro)phages are common in the FGT and vaginal bacteria and (pro)phages may interact.

Unveiling Ecological and Genetic Novelty within Lytic and Lysogenic Viral Communities of Hot Spring Phototrophic Microbial Mats

Viruses exert diverse ecosystem impacts by controlling their host community through lytic predator-prey dynamics. However, the mechanisms by which lysogenic viruses influence their host-microbial community are less clear. In hot springs, lysogeny is considered an active lifestyle, yet it has not been systematically studied in all habitats, with phototrophic microbial mats (PMMs) being particularly not studied. We carried out viral metagenomics following in situ mitomycin C induction experiments in PMMs from Porcelana hot spring (Northern Patagonia, Chile). The compositional changes of viral communities at two different sites were analyzed at the genomic and gene levels. Furthermore, the presence of integrated prophage sequences in environmental metagenome-assembled genomes from published Porcelana PMM metagenomes was analyzed. Our results suggest that virus-specific replicative cycles (lytic and lysogenic) were associated with specific host taxa with different metabolic capacities. One of the most abundant lytic viral groups corresponded to cyanophages, which would infect the cyanobacteria Fischerella, the most active and dominant primary producer in thermophilic PMMs. Likewise, lysogenic viruses were related exclusively to chemoheterotrophic bacteria from the phyla Proteobacteria, Firmicutes, and Actinobacteria. These temperate viruses possess accessory genes to sense or control stress-related processes in their hosts, such as sporulation and biofilm formation. Taken together, these observations suggest a nexus between the ecological role of the host (metabolism) and the type of viral lifestyle in thermophilic PMMs. This has direct implications in viral ecology, where the lysogenic-lytic switch is determined by nutrient abundance and microbial density but also by the metabolism type that prevails in the host community.

IMPORTANCE Hot springs harbor microbial communities dominated by a limited variety of microorganisms and, as such, have become a model for studying community ecology and understanding how biotic and abiotic interactions shape their structure. Viruses in hot springs are shown to be ubiquitous, numerous, and active components of these communities. However, lytic and lysogenic viral communities of thermophilic phototrophic microbial mats (PMMs) remain largely unexplored. In this work, we use the power of viral metagenomics to reveal changes in the viral community following a mitomycin C induction experiment in PMMs. The importance of our research is that it will improve our understanding of viral lifestyles in PMMs via exploring the differences in the composition of natural and induced viral communities at the genome and gene levels. This novel information will contribute to deciphering which biotic and abiotic factors may control the transitions between lytic and lysogenic cycles in these extreme environments.

Perspective on taxonomic classification of uncultivated viruses

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

Unexpected myriad of co-occurring viral strains and species in one of the most abundant and microdiverse viruses on Earth

Viral genetic microdiversity drives adaptation, pathogenicity, and speciation and has critical consequences for the viral-host arms race occurring at the strain and species levels, which ultimately impact microbial community structure and biogeochemical cycles. Despite the fact that most efforts have focused on viral macrodiversity, little is known about the microdiversity of ecologically important viruses on Earth. Recently, single-virus genomics discovered the putatively most abundant ocean virus in temperate and tropical waters: the uncultured dsDNA virus vSAG 37-F6 infecting Pelagibacter, the most abundant marine bacteria. In this study, we report the cooccurrence of up to ≈1,500 different viral strains (>95% nucleotide identity) and ≈30 related species (80-95% nucleotide identity) in a single oceanic sample. Viral microdiversity was maintained over space and time, and most alleles were the result of synonymous mutations without any apparent adaptive benefits to cope with host translation codon bias and efficiency. Gene flow analysis used to delimitate species according to the biological species concept (BSC) revealed the impact of recombination in shaping vSAG 37-F6 virus and Pelagibacter speciation. Data demonstrated that this large viral microdiversity somehow mirrors the host species diversity since ≈50% of the 926 analyzed Pelagibacter genomes were found to belong to independent BSC species that do not significantly engage in gene flow with one another. The host range of this evolutionarily successful virus revealed that a single viral species can infect multiple Pelagibacter BSC species, indicating that this virus crosses not only formal BSC barriers but also biomes since viral ancestors are found in freshwater.

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.

Genome and Ecology of a Novel Alteromonas Podovirus, ZP6, Representing a New Viral Genus, Mareflavirus

Alteromonas is a ubiquitous, abundant, copiotrophic and phytoplankton-associated marine member of the Gammaproteobacteria with a range extending from tropical waters to polar regions and including hadal zones. Here, we describe a novel Alteromonas phage, ZP6, that was isolated from surface coastal waters of Qingdao, China. ZP6 contains a linear, double-stranded, 38,080-bp DNA molecule with 50.1% G+C content and 47 putative open reading frames (ORFs). Three auxiliary metabolic genes were identified, encoding metal-dependent phosphohydrolase, diaminopurine synthetase, and nucleotide pyrophosphohydrolase. The first two ORFs facilitate the replacement of adenine (A) by diaminopurine (Z) in phage genomes and help phages to evade attack from host restriction enzymes. The nucleotide pyrophosphohydrolase enables the host cells to stop programmed cell death and improves the survival rate of the host in a nutrient-depleted environment. Phylogenetic analysis based on the amino acid sequences of whole genomes and comparative genomic analysis revealed that ZP6 is most closely related to Enhodamvirus but with low similarity (shared genes, <30%, and average nucleotide sequence identity, <65%); it is distinct from other bacteriophages. Together, these results suggest that ZP6 could represent a novel viral genus, here named Mareflavirus. Combining its ability to infect Alteromonas, its harboring of a diaminopurine genome-biosynthetic system, and its representativeness of an understudied viral group, ZP6 could be an important and novel model system for marine virus research.

IMPORTANCEAlteromonas is an important symbiotic bacterium of phytoplankton, but research on its bacteriophages is still at an elementary level. Our isolation and genome characterization of a novel Alteromonas podovirus, ZP6, identified a new viral genus of podovirus, namely, Mareflavirus. The ZP6 genome, with a diaminopurine genome-biosynthetic system, is different from those of other isolated Alteromonas phages and will bring new impetus to the development of virus classification and provide important insights into novel viral sequences from metagenomic data sets.

MVIP: multi-omics portal of viral infection

Virus infections are huge threats to living organisms and cause many diseases, such as COVID-19 caused by SARS-CoV-2, which has led to millions of deaths. To develop effective strategies to control viral infection, we need to understand its molecular events in host cells. Virus related functional genomic datasets are growing rapidly, however, an integrative platform for systematically investigating host responses to viruses is missing. Here, we developed a user-friendly multi-omics portal of viral infection named as MVIP (https://mvip.whu.edu.cn/). We manually collected available high-throughput sequencing data under viral infection, and unified their detailed metadata including virus, host species, infection time, assay, and target, etc. We processed multi-layered omics data of more than 4900 viral infected samples from 77 viruses and 33 host species with standard pipelines, including RNA-seq, ChIP-seq, and CLIP-seq, etc. In addition, we integrated these genome-wide signals into customized genome browsers, and developed multiple dynamic charts to exhibit the information, such as time-course dynamic and differential gene expression profiles, alternative splicing changes and enriched GO/KEGG terms. Furthermore, we implemented several tools for efficiently mining the virus-host interactions by virus, host and genes. MVIP would help users to retrieve large-scale functional information and promote the understanding of virus-host interactions.

Viral Characteristics of the Warm Atlantic and Cold Arctic Water Masses in the Nordic Seas

Nordic Seas are the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean with complex water masses, experiencing an abrupt climate change. Though knowledge of the marine virosphere has expanded rapidly, the diversity of viruses and their relationships with host cells and water masses in the Nordic Seas remain to be fully revealed. Here, we establish the Nordic Sea DNA virome (NSV) data set of 55,315 viral contigs including 1,478 unique viral populations from seven stations influenced by both the warm Atlantic and cold Arctic water masses. Caudovirales dominated in the seven NSVs, especially in the warm Atlantic waters. The major giant nucleocytoplasmic large DNA viruses (NCLDVs) contributed a significant proportion of the classified viral contigs in the NSVs (32.2%), especially in the cold Arctic waters (44.9%). The distribution patterns of Caudovirales and NCLDVs were a reflection of the community structure of their hosts in the corresponding water masses and currents. Latitude, pH, and flow speed were found to be key factors influencing the microbial communities and coinfluencing the variation of viral communities. Network analysis illustrated the tight coupling between the variation of viral communities and microbial communities in the Nordic Seas. This study suggests a probable linkage between viromes, host cells, and surface water masses from both the cool Arctic and warm Atlantic Oceans. IMPORTANCE This is a systematic study of Nordic Sea viromes using metagenomic analysis. The viral diversity, community structure, and their relationships with host cells and the complex water masses from both the cool Arctic and the warm Atlantic oceans were illustrated. The NCLDVs and Caudovirales are proposed as the viral characteristics of the cold Arctic and warm Atlantic waters, respectively. This study provides an important background for the viromes in the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean and sheds light on their responses to abrupt climate change in the future.

A Previously Undescribed Highly Prevalent Phage Identified in a Danish Enteric Virome Catalog

Gut viruses are important, yet often neglected, players in the complex human gut microbial ecosystem. Recently, the number of human gut virome studies has been increasing; however, we are still only scratching the surface of the immense viral diversity. In this study, 254 virus-enriched fecal metagenomes from 204 Danish subjects were used to generate the Danish Enteric Virome Catalog (DEVoC) containing 12,986 nonredundant viral scaffolds, of which the majority was previously undescribed, encoding 190,029 viral genes. The DEVoC was used to compare 91 healthy DEVoC gut viromes from children, adolescents, and adults that were used to create the DEVoC. Gut viromes of healthy Danish subjects were dominated by phages. While most phage genomes (PGs) only occurred in a single subject, indicating large virome individuality, 39 PGs were present in more than 10 healthy subjects. Among these 39 PGs, the prevalences of three PGs were associated with age. To further study the prevalence of these 39 prevalent PGs, 1,880 gut virome data sets of 27 studies from across the world were screened, revealing several age-, geography-, and disease-related prevalence patterns. Two PGs also showed a remarkably high prevalence worldwide-a crAss-like phage (20.6% prevalence), belonging to the tentative AlphacrAssvirinae subfamily, and a previously undescribed circular temperate phage infecting Bacteroides dorei (14.4% prevalence), called LoVEphage because it encodes lots of viral elements. Due to the LoVEphage's high prevalence and novelty, public data sets in which the LoVEphage was detected were de novo assembled, resulting in an additional 18 circular LoVEphage-like genomes (67.9 to 72.4 kb). IMPORTANCE Through generation of the DEVoC, we added numerous previously uncharacterized viral genomes and genes to the ever-increasing worldwide pool of human gut viromes. The DEVoC, the largest human gut virome catalog generated from consistently processed fecal samples, facilitated the analysis of the 91 healthy Danish gut viromes. Characterizing the biggest cohort of healthy gut viromes from children, adolescents, and adults to date confirmed the previously established high interindividual variation in human gut viromes and demonstrated that the effect of age on the gut virome composition was limited to the prevalence of specific phage (groups). The identification of a previously undescribed prevalent phage illustrates the usefulness of developing virome catalogs, and we foresee that the DEVoC will benefit future analysis of the roles of gut viruses in human health and disease.

Comprehensive discovery of CRISPR-targeted terminally redundant sequences in the human gut metagenome: Viruses, plasmids, and more

Viruses are the most numerous biological entity, existing in all environments and infecting all cellular organisms. Compared with cellular life, the evolution and origin of viruses are poorly understood; viruses are enormously diverse, and most lack sequence similarity to cellular genes. To uncover viral sequences without relying on either reference viral sequences from databases or marker genes that characterize specific viral taxa, we developed an analysis pipeline for virus inference based on clustered regularly interspaced short palindromic repeats (CRISPR). CRISPR is a prokaryotic nucleic acid restriction system that stores the memory of previous exposure. Our protocol can infer CRISPR-targeted sequences, including viruses, plasmids, and previously uncharacterized elements, and predict their hosts using unassembled short-read metagenomic sequencing data. By analyzing human gut metagenomic data, we extracted 11,391 terminally redundant CRISPR-targeted sequences, which are likely complete circular genomes. The sequences included 2,154 tailed-phage genomes, together with 257 complete crAssphage genomes, 11 genomes larger than 200 kilobases, 766 genomes of Microviridae species, 56 genomes of Inoviridae species, and 95 previously uncharacterized circular small genomes that have no reliably predicted protein-coding gene. We predicted the host(s) of approximately 70% of the discovered genomes at the taxonomic level of phylum by linking protospacers to taxonomically assigned CRISPR direct repeats. These results demonstrate that our protocol is efficient for de novo inference of CRISPR-targeted sequences and their host prediction.

The evolution and origins of viruses are long-standing questions in the field of biology. Viral genomes provide fundamental information to infer the evolution and origin of viruses. However, viruses are extraordinarily diverse, and there are no single genes shared across entire species. Several methods were developed to collect viral genomes from metagenome. To infer viral genomes from metagenome, previous approaches relied on reference viral genomes. We thought that such reference-based methods may not be sufficient to uncover diverse viral genomes; therefore, we developed a pipeline that utilizes CRISPR, a prokaryotic adaptive immunological memory. Using this pipeline, we discovered more than 10,000 positively complete CRISPR-targeted genomes from human gut metagenome datasets. A substantial portion of the discovered genomes encoded various types of capsid proteins, supporting the contention that these sequences are viral. Although the majority of these capsid-protein-coding sequences were previously characterized, we notably discovered Inoviridae genomes that were previously difficult to infer as being viral. Furthermore, some of the remaining unclassified sequences without a detectable capsid-protein-encoding gene had a notably low protein-coding ratio. Overall, our pipeline successfully discovered viruses and previously uncharacterized presumably mobile genetic elements targeted by CRISPR.

Characterizing Enterotypes in Human Metagenomics: A Viral Perspective

The diversity and high genomic mutation rates of viral species hinder our understanding of viruses and their contributions to human health. Viral enterotypes as a description of the gut virome, its characteristics have not been thoroughly studied. Here we investigated the human gut virome composition using previously published sequencing data of 2,690 metagenomes from seven countries with various phenotypes. We found that the virome was dominated by double-stranded DNA viruses in our data, and young children and adults showed different stages in their fecal enterovirus composition. Beta diversity showed there were significantly less homogeneous in individuals with severe disorders of bile acid secretion, such as cirrhosis. In contrast, there were no significant differences in distances to centroids or viral components between patients with phenotypes unrelated to bile acid, such as hypertension. Enterotypes determined independently from various projects showed similar specific viruses and enrichment direction. Confounding factors, such as different sequencing platforms and library construction, did not confuse enterotyping. The gut virome composition pattern could be described by two viral enterotypes, which supported a discrete, rather than a gradient, distribution. Three main components, enterotype 1 and 2 specific viruses and the other, comprise the total viral variation in these sets. Compared with enterotype 2, enterotype 1 had a higher viral count, Shannon index, and similarity between samples. The relative abundance of enterotype-specific viruses is a crucial determinant of enterotype assignment. Samples not matching any of the defined enterotypes in the database did not necessarily correlate to sickness. Therefore, the background context must be carefully considered when using a viral enterotype as a feature for disease prediction. Our results highlight important insights into the human gut virome composition by exploring two-main viral enterotypes in population and providing an alternate covariate for early disease screening.

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.

The novel anti-CRISPR AcrIIA22 relieves DNA torsion in target plasmids and impairs SpyCas9 activity

To overcome CRISPR-Cas defense systems, many phages and mobile genetic elements (MGEs) encode CRISPR-Cas inhibitors called anti-CRISPRs (Acrs). Nearly all characterized Acrs directly bind Cas proteins to inactivate CRISPR immunity. Here, using functional metagenomic selection, we describe AcrIIA22, an unconventional Acr found in hypervariable genomic regions of clostridial bacteria and their prophages from human gut microbiomes. AcrIIA22 does not bind strongly to SpyCas9 but nonetheless potently inhibits its activity against plasmids. To gain insight into its mechanism, we obtained an X-ray crystal structure of AcrIIA22, which revealed homology to PC4-like nucleic acid-binding proteins. Based on mutational analyses and functional assays, we deduced that acrIIA22 encodes a DNA nickase that relieves torsional stress in supercoiled plasmids. This may render them less susceptible to SpyCas9, which uses free energy from negative supercoils to form stable R-loops. Modifying DNA topology may provide an additional route to CRISPR-Cas resistance in phages and MGEs.

Long-read metagenomics of multiple displacement amplified DNA of low-biomass human gut phageomes by SACRA pre-processing chimeric reads

The human gut bacteriophage community (phageome) plays an important role in the host's health and disease; however, the entire structure is poorly understood, partly owing to the generation of many incomplete genomes in conventional short-read metagenomics. Here, we show long-read metagenomics of amplified DNA of low-biomass phageomes with multiple displacement amplification (MDA), involving the development of a novel bioinformatics tool, split amplified chimeric read algorithm (SACRA), that efficiently pre-processed numerous chimeric reads generated through MDA. Using five samples, SACRA markedly reduced the average chimera ratio from 72% to 1.5% in PacBio reads with an average length of 1.8 kb. De novo assembly of chimera-less PacBio long reads reconstructed contigs of ≥5 kb with an average proportion of 27%, which was 1% in contigs from MiSeq short reads, thereby dramatically improving contig length and genome completeness. Comparison of PacBio and MiSeq contigs found MiSeq contig fragmentations frequently near local repeats and hypervariable regions in the phage genomes, and those caused by multiple homologous phage genomes coexisting in the community. We also developed a reference-independent method to assess the completeness of the linear phage genomes. Overall, we established a SACRA-coupled long-read metagenomics robust to highly diverse gut phageomes, identifying high-quality circular and linear phage genomes with adequate sequence quantity.

Taxonomy-aware, sequence similarity ranking reliably predicts phage-host relationships

Background

Characterizing phage–host interactions is critical to understanding the ecological role of both partners and effective isolation of phage therapeuticals. Unfortunately, experimental methods for studying these interactions are markedly slow, low-throughput, and unsuitable for phages or hosts difficult to maintain in laboratory conditions. Therefore, a number of in silico methods emerged to predict prokaryotic hosts based on viral sequences. One of the leading approaches is the application of the BLAST tool that searches for local similarities between viral and microbial genomes. However, this prediction method has three major limitations: (i) top-scoring sequences do not always point to the actual host; (ii) mosaic virus genomes may match to many, typically related, bacteria; and (iii) viral and host sequences may diverge beyond the point where their relationship can be detected by a BLAST alignment.

Results

We created an extension to BLAST, named Phirbo, that improves host prediction quality beyond what is obtainable from standard BLAST searches. The tool harnesses information concerning sequence similarity and bacteria relatedness to predict phage–host interactions. Phirbo was evaluated on three benchmark sets of known virus–host pairs, and it improved precision and recall by 11–40 percentage points over currently available, state-of-the-art, alignment-based, alignment-free, and machine-learning host prediction tools. Moreover, the discriminatory power of Phirbo for the recognition of virus–host relationships surpassed the results of other tools by at least 10 percentage points (area under the curve = 0.95), yielding a mean host prediction accuracy of 57% and 68% at the genus and family levels, respectively, and drops by 12 percentage points when using only a fraction of viral genome sequences (3 kb). Finally, we provide insights into a repertoire of protein and ncRNA genes that are shared between phages and hosts and may be prone to horizontal transfer during infection.

Conclusions

Our results suggest that Phirbo is a simple and effective tool for predicting phage–host relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01146-6.

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.

Mechanism for Cas4-assisted directional spacer acquisition in CRISPR-Cas

Prokaryotes adapt to challenges from mobile genetic elements by integrating spacers derived from foreign DNA in the CRISPR array¹. Spacer insertion is carried out by the Cas1-Cas2 integrase complex^2-4. A substantial fraction of CRISPR-Cas systems use a Fe-S cluster containing Cas4 nuclease to ensure that spacers are acquired from DNA flanked by a protospacer adjacent motif (PAM)^5,6 and inserted into the CRISPR array unidirectionally, so that the transcribed CRISPR RNA can guide target searching in a PAM-dependent manner. Here we provide a high-resolution mechanistic explanation for the Cas4-assisted PAM selection, spacer biogenesis and directional integration by type I-G CRISPR in Geobacter sulfurreducens, in which Cas4 is naturally fused with Cas1, forming Cas4/Cas1. During biogenesis, only DNA duplexes possessing a PAM-embedded 3'-overhang trigger Cas4/Cas1-Cas2 assembly. During this process, the PAM overhang is specifically recognized and sequestered, but is not cleaved by Cas4. This 'molecular constipation' prevents the PAM-side prespacer from participating in integration. Lacking such sequestration, the non-PAM overhang is trimmed by host nucleases and integrated to the leader-side CRISPR repeat. Half-integration subsequently triggers PAM cleavage and Cas4 dissociation, allowing spacer-side integration. Overall, the intricate molecular interaction between Cas4 and Cas1-Cas2 selects PAM-containing prespacers for integration and couples the timing of PAM processing with the stepwise integration to establish directionality.

PAM-repeat associations and spacer selection preferences in single and co-occurring CRISPR-Cas systems

BACKGROUND

The adaptive CRISPR-Cas immune system stores sequences from past invaders as spacers in CRISPR arrays and thereby provides direct evidence that links invaders to hosts. Mapping CRISPR spacers has revealed many aspects of CRISPR-Cas biology, including target requirements such as the protospacer adjacent motif (PAM). However, studies have so far been limited by a low number of mapped spacers in the database.

RESULTS

By using vast metagenomic sequence databases, we map approximately one-third of more than 200,000 unique CRISPR spacers from a variety of microbes and derive a catalog of more than two hundred unique PAM sequences associated with specific CRISPR-Cas subtypes. These PAMs are further used to correctly assign the orientation of CRISPR arrays, revealing conserved patterns between the last nucleotides of the CRISPR repeat and PAM. We could also deduce CRISPR-Cas subtype-specific preferences for targeting either template or coding strand of open reading frames. While some DNA-targeting systems (type I-E and type II systems) prefer the template strand and avoid mRNA, other DNA- and RNA-targeting systems (types I-A and I-B and type III systems) prefer the coding strand and mRNA. In addition, we find large-scale evidence that both CRISPR-Cas adaptation machinery and CRISPR arrays are shared between different CRISPR-Cas systems. This could lead to simultaneous DNA and RNA targeting of invaders, which may be effective at combating mobile genetic invaders.

CONCLUSIONS

This study has broad implications for our understanding of how CRISPR-Cas systems work in a wide range of organisms for which only the genome sequence is known.