Marine DNA Viral Macro- and Microdiversity from Pole to Pole

Characteristics of a pseudolysogenic phage vB_YpM_HQ103 infecting Yersinia pestis

The plague, caused by Yersinia pestis, is a natural focal disease and the presence of Y. pestis in the environment is a critical ecological concern worldwide. The role of Y. pestis phages in the ecological life cycle of the plague is crucial. Previously, a temperature-sensitive phage named vB_YpM_HQ103 was isolated from plague foci in Yunnan province, China. Upon infecting the EV76 strain of Y. pestis, vB_YpM_HQ103 exhibits lysogenic behavior at 21 °C and lytic behavior at 37 °C. Various methods including continuous passage lysogenic tests, in vitro lysis tests, comparative genomic assays, fluorescence quantitative PCR and receptor identification tests were employed to demonstrate that the lysogenic life cycle of this phage is applicable to wild Y. pestis strains; its lysogeny is pseudolysogenic (carrying but not integrating), allowing it to replicate and proliferate within Y. pestis. Furthermore, we have identified the outer membrane protein OmpA of Y. pestis as the receptor for phage infection. In conclusion, our research provides insight into the characteristics and receptors of a novel Y. pestis phage infection with a pseudolysogenic cycle. The findings of this study enhance our understanding of Y. pestis phages and plague microecology, offering valuable insights for future studies on the conservation and genetic evolution of Y. pestis in nature.

The direct and indirect drivers shaping RNA viral communities in grassland soils

Recent studies have revealed diverse RNA viral communities in soils. Yet, how environmental factors influence soil RNA viruses remains largely unknown. Here, we recovered RNA viral communities from bulk metatranscriptomes sequenced from grassland soils managed for 5 years under multiple environmental conditions including water content, plant presence, cultivar type, and soil depth. More than half of the unique RNA viral contigs (64.6%) were assigned with putative hosts. About 74.7% of these classified RNA viral contigs are known as eukaryotic RNA viruses suggesting eukaryotic RNA viruses may outnumber prokaryotic RNA viruses by nearly three times in this grassland. Of the identified eukaryotic RNA viruses and the associated eukaryotic species, the most dominant taxa were Mitoviridae with an average relative abundance of 72.4%, and their natural hosts, Fungi with an average relative abundance of 56.6%. Network analysis and structural equation modeling support that soil water content, plant presence, and type of cultivar individually demonstrate a significant positive impact on eukaryotic RNA viral richness directly as well as indirectly on eukaryotic RNA viral abundance via influencing the co-existing eukaryotic members. A significant negative influence of soil depth on soil eukaryotic richness and abundance indirectly impacts soil eukaryotic RNA viral communities. These results provide new insights into the collective influence of multiple environmental and community factors that shape soil RNA viral communities and offer a structured perspective of how RNA virus diversity and ecology respond to environmental changes.

IMPORTANCE

Climate change has been reshaping the soil environment as well as the residing microbiome. This study provides field-relevant information on how environmental and community factors collectively shape soil RNA communities and contribute to ecological understanding of RNA viral survival under various environmental conditions and virus-host interactions in soil. This knowledge is critical for predicting the viral responses to climate change and the potential emergence of biothreats.

A novel flavobacterial phage abundant during green tide, representing a new viral family, Zblingviridae

Flavobacteriia are the dominant and active bacteria during algal blooms and play an important role in polysaccharide degradation. However, little is known about phages infecting Flavobacteriia, especially during green tide. In this study, a novel virus, vB_TgeS_JQ, infecting Flavobacteriia was isolated from the surface water of the Golden Beach of Qingdao, China. Transmission electron microscopy demonstrated that vB_TgeS_JQ had the morphology of siphovirus. The experiments showed that it was stable from -20°C to 45°C and pH 5 to pH 8, with latent and burst periods both lasting for 20 min. Genomic analysis showed that the phage vB_TgeS_JQ contained a 40,712-bp dsDNA genome with a GC content of 30.70%, encoding 74 open-reading frames. Four putative auxiliary metabolic genes were identified, encoding electron transfer-flavoprotein dehydrogenase, calcineurin-like phosphoesterase, phosphoribosyl-ATP pyrophosphohydrolase, and TOPRIM nucleotidyl hydrolase. The abundance of phage vB_TgeS_JQ was higher during Ulva prolifera (U. prolifera) blooms compared with other marine environments. The phylogenetic and comparative genomic analyses revealed that vB_TgeS_JQ exhibited significant differences from all other phage isolates in the databases and therefore was classified as an undiscovered viral family, named Zblingviridae. In summary, this study expands the knowledge about the genomic, phylogenetic diversity and distribution of flavobacterial phages (flavophages), especially their roles during U. prolifera blooms.

IMPORTANCE

The phage vB_TgeS_JQ was the first flavobacterial phage isolated during green tide, representing a new family in Caudoviricetes and named Zblingviridae. The abundance of phage vB_TgeS_JQ was higher during the Ulva prolifera blooms. This study provides insights into the genomic, phylogenetic diversity, and distribution of flavophages, especially their roles during U. prolifera blooms.

Multikingdom characterization of gut microbiota in patients with rheumatoid arthritis and rheumatoid arthritis-associated interstitial lung disease

Rheumatoid arthritis-associated interstitial lung disease (RA-ILD) is a serious and common extra-articular disease manifestation. Patients with RA-ILD experience reduced bacterial diversity and gut bacteriome alterations. However, the gut mycobiome and virome in these patients have been largely neglected. In this study, we performed whole-metagenome shotgun sequencing on fecal samples from 30 patients with RA-ILD, and 30 with RA-non-ILD, and 40 matched healthy controls. The gut bacteriome and mycobiome were explored using a reference-based approach, while the gut virome was profiled based on a nonredundant viral operational taxonomic unit (vOTU) catalog. The results revealed significant alterations in the gut microbiomes of both RA-ILD and RA-non-ILD groups compared with healthy controls. These alterations encompassed changes in the relative abundances of 351 bacterial species, 65 fungal species, and 4,367 vOTUs. Bacteria such as Bifidobacterium longum, Dorea formicigenerans, and Collinsella aerofaciens were enriched in both patient groups. Ruminococcus gnavus (RA-ILD), Gemmiger formicilis, and Ruminococcus bromii (RA-non-ILD) were uniquely enriched. Conversely, Faecalibacterium prausnitzii, Bacteroides spp., and Roseburia inulinivorans showed depletion in both patient groups. Mycobiome analysis revealed depletion of certain fungi, including Saccharomyces cerevisiae and Candida albicans, in patients with RA compared with healthy subjects. Notably, gut virome alterations were characterized by an increase in Siphoviridae and a decrease in Myoviridae, Microviridae, and Autographiviridae in both patient groups. Hence, multikingdom gut microbial signatures showed promise as diagnostic indicators for both RA-ILD and RA-non-ILD. Overall, this study provides comprehensive insights into the fecal virome, bacteriome, and mycobiome landscapes of RA-ILD and RA-non-ILD gut microbiota, thereby offering potential biomarkers for further mechanistic and clinical research.

Novel phage infecting the Roseobacter CHUG lineage reveals a diverse and globally distributed phage family

Bacteriophages play an essential role in shaping the diversity and metabolism of bacterial communities. Marine Roseobacter group is an abundant heterotrophic bacterial group that is involved in many major element cycles, especially carbon and sulfur. Members of the Roseobacter CHUG (Clade Hidden and Underappreciated Globally) lineage are globally distributed and are activated in pelagic marine environments. In this study, we isolated and characterized a phage, CRP-810, that infects the CHUG strain FZCC0198. The genome of CRP-810 was dissimilar to those of other known phages. Additionally, 251 uncultured viral genomes (UViGs) closely related to CRP-810 were obtained from the uncultivated marine viral contig databases. Comparative genomic and phylogenetic analyses revealed that CRP-810 and these related UViGs exhibited conserved genome synteny, representing a new phage family with at least eight subgroups. Most of the CRP-810-type phages contain an integrase gene, and CRP-810 can be integrated into the host genome. Further analysis revealed that three CRP-810-type members were prophages found in the genomes of marine SAR11, Poseidonocella, and Sphingomonadaceae. Finally, viromic read-mapping analysis showed that CRP-810-type phages were globally distributed and displayed distinct biogeographic patterns related to temperature and latitude. Many members with a lower G + C content were mainly distributed in the trade station, whereas members with a higher G + C content were mainly distributed in polar and westerlies station, indicating that the niche differentiation of phages was subject to host adaptation. Collectively, these findings identify a novel phage family and expand our understanding of phylogenetic diversity, evolution, and biogeography of marine phages.

IMPORTANCE

The Roseobacter CHUG lineage, affiliated with the Pelagic Roseobacter Cluster (PRC), is widely distributed in the global oceans and is active in oligotrophic seawater. However, knowledge of the bacteriophages that infect CHUG members is limited. In this study, a CHUG phage, CRP-810, that infects the CHUG strain FZCC0198, was isolated and shown to have a novel genomic architecture. In addition, 251 uncultured viral genomes closely related to CRP-810 were recovered and included in the analyses. Phylogenomic analyses revealed that the CRP-810-type phages represent a new phage family containing at least eight genus-level subgroups. Members of this family were predicted to infect various marine bacteria. We also demonstrated that the CRP-810-type phages are widely distributed in global oceans and display distinct biogeographic patterns related to latitude. Collectively, this study provides important insights into the genomic organization, diversity, and ecology of a novel phage family that infect ecologically important bacteria in the global ocean.

Sequencing-based analysis of microbiomes

Microbiomes occupy a range of niches and, in addition to having diverse compositions, they have varied functional roles that have an impact on agriculture, environmental sciences, and human health and disease. The study of microbiomes has been facilitated by recent technological and analytical advances, such as cheaper and higher-throughput DNA and RNA sequencing, improved long-read sequencing and innovative computational analysis methods. These advances are providing a deeper understanding of microbiomes at the genomic, transcriptional and translational level, generating insights into their function and composition at resolutions beyond the species level.

Temporal turnover of viral biodiversity and functional potential in intertidal wetlands

As the central members of the microbiome networks, viruses regulate the composition of microbial communities and drive the nutrient cycles of ecosystems by lysing host cells. Therefore, uncovering the dynamic patterns and the underlying ecological mechanisms mediating the tiniest viral communities across space and through time in natural ecosystems is of crucial importance for better understanding the complex microbial world. Here, the temporal dynamics of intertidal viral communities were investigated via a time-series sampling effort. A total of 1911 viral operational taxonomic units were recovered from 36 bimonthly collected shotgun metagenomes. Functionally important auxiliary metabolic genes involved in carbohydrate, sulfur, and phosphorus metabolism were detected, some of which (e.g., cysH gene) were stably present within viral genomes over time. Over the sampling period, strong and comparable temporal turnovers were observed for intertidal viromes and their host microbes. Winter was determined as the pivotal point for the shifts in viral diversity patterns. Notably, the viral micro-diversity covaried with the macro-diversity, following similar temporal patterns. The relative abundances of viral taxa also covaried with their host prokaryotes. Meanwhile, the virus-host relationships at the whole community level were relatively stable. Further statistical analyses demonstrated that the dynamic patterns of viral communities were highly deterministic, for which temperature was the major driver. This study provided valuable mechanistic insights into the temporal turnover of viral communities in complex ecosystems such as intertidal wetlands.

A genomic compendium of cultivated human gut fungi characterizes the gut mycobiome and its relevance to common diseases

The gut fungal community represents an essential element of human health, yet its functional and metabolic potential remains insufficiently elucidated, largely due to the limited availability of reference genomes. To address this gap, we presented the cultivated gut fungi (CGF) catalog, encompassing 760 fungal genomes derived from the feces of healthy individuals. This catalog comprises 206 species spanning 48 families, including 69 species previously unidentified. We explored the functional and metabolic attributes of the CGF species and utilized this catalog to construct a phylogenetic representation of the gut mycobiome by analyzing over 11,000 fecal metagenomes from Chinese and non-Chinese populations. Moreover, we identified significant common disease-related variations in gut mycobiome composition and corroborated the associations between fungal signatures and inflammatory bowel disease (IBD) through animal experimentation. These resources and findings substantially enrich our understanding of the biological diversity and disease relevance of the human gut mycobiome.

Exploring the roles of ribosomal peptides in prokaryote-phage interactions through deep learning-enabled metagenome mining

BACKGROUND

Microbial secondary metabolites play a crucial role in the intricate interactions within the natural environment. Among these metabolites, ribosomally synthesized and post-translationally modified peptides (RiPPs) are becoming a promising source of therapeutic agents due to their structural diversity and functional versatility. However, their biosynthetic capacity and ecological functions remain largely underexplored.

RESULTS

Here, we aim to explore the biosynthetic profile of RiPPs and their potential roles in the interactions between microbes and viruses in the ocean, which encompasses a vast diversity of unique biomes that are rich in interactions and remains chemically underexplored. We first developed TrRiPP to identify RiPPs from ocean metagenomes, a deep learning method that detects RiPP precursors in a hallmark gene-independent manner to overcome the limitations of classic methods in processing highly fragmented metagenomic data. Applying this method to metagenomes from the global ocean microbiome, we uncover a diverse array of previously uncharacterized putative RiPP families with great novelty and diversity. Through correlation analysis based on metatranscriptomic data, we observed a high prevalence of antiphage defense-related and phage-related protein families that were co-expressed with RiPP families. Based on this putative association between RiPPs and phage infection, we constructed an Ocean Virus Database (OVD) and established a RiPP-involving host-phage interaction network through host prediction and co-expression analysis, revealing complex connectivities linking RiPP-encoding prokaryotes, RiPP families, viral protein families, and phages. These findings highlight the potential of RiPP families involved in prokaryote-phage interactions and coevolution, providing insights into their ecological functions in the ocean microbiome.

CONCLUSIONS

This study provides a systematic investigation of the biosynthetic potential of RiPPs from the ocean microbiome at a global scale, shedding light on the essential insights into the ecological functions of RiPPs in prokaryote-phage interactions through the integration of deep learning approaches, metatranscriptomic data, and host-phage connectivity. This study serves as a valuable example of exploring the ecological functions of bacterial secondary metabolites, particularly their associations with unexplored microbial interactions. Video Abstract.

Profiles of phage in global hospital wastewater: Association with microbial hosts, antibiotic resistance genes, metal resistance genes, and mobile genetic elements

Hospital wastewater (HWW) is known to host taxonomically diverse microbial communities, yet limited information is available on the phages infecting these microorganisms. To fill this knowledge gap, we conducted an in-depth analysis using 377 publicly available HWW metagenomic datasets from 16 countries across 4 continents in the NCBI SRA database to elucidate phage-host dynamics and phage contributions to resistance gene transmission. We first assembled a metagenomic HWW phage catalog comprising 13,812 phage operational taxonomic units (pOTUs). The majority of these pOTUs belonged to the Caudoviricetes order, representing 75.29 % of this catalog. Based on the lifestyle of phages, we found that potentially virulent phages predominated in HWW. Specifically, 583 pOTUs have been predicted to have the capability to lyse 81 potentially pathogenic bacteria, suggesting the promising role of HWW phages as a viable alternative to antibiotics. Among all pOTUs, 1.56 % of pOTUs carry 108 subtypes of antibiotic resistance genes (ARGs), 0.96 % of pOTUs carry 76 subtypes of metal resistance genes (MRGs), and 0.96 % of pOTUs carry 22 subtypes of non-phage mobile genetic elements (MGEs). Predictions indicate that certain phages carrying ARGs, MRGs, and non-phage MGEs could infect bacteria hosts, even potential pathogens. This suggests that phages in HWW may contribute to the dissemination of resistance-associated genes in the environment. This meta-analysis provides the first global catalog of HWW phages, revealing their correlations with microbial hosts and pahge-associated ARGs, MRG, and non-phage MGEs. The insights gained from this research hold promise for advancing the applications of phages in medical and industrial contexts.

Long-read powered viral metagenomics in the oligotrophic Sargasso Sea

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

Phages in different habitats and their ability to carry antibiotic resistance genes

As the most abundant organisms on Earth, phages play a key role in the evolution of bacterial antibiotic resistance. Although previous studies have demonstrated the molecular mechanisms of horizontal gene transfer mediated by mobile genetic elements, our understanding of the intertwined relationships between antibiotic resistance genes (ARGs) and phages is limited. In this study, we analysed 2781 metagenomic samples to reveal the composition and species interactions of phage communities in different habitats as well as their capacity to carry ARGs with health risks. The composition of phage communities varies in different habitats and mainly depends on environmental conditions. Terrestrial habitats display more complex and robust interactions between phages than aquatic and human-associated habitats, resulting in the highest biodiversity of phages. Several types of phages in certain taxa (4.95-7.67%, mainly belonging to Caudoviricetes) have the capacity to carry specific ARGs and display a high potential risk to human health, especially in human-associated habitats. Overall, our results provide insights into the assembly mechanisms of phage communities and their effects on the dissemination of antibiotic resistance.

The bacteriome-coupled phage communities continuously contract and shift to orchestrate the traditional rice vinegar fermentation

Amounts of microbiome studies have uncovered the microbial communities of traditional food fermentations, while in which the phageome development with time is poorly understood. Here, we conducted a study to decipher both phageome and bacteriome of the traditional rice vinegar fermentation. The vinegar phageomes showed significant differences in the alpha diversity, network density and clustering coefficient over time. Peduoviridae had the highest relative abundance. Moreover, the phageome negatively correlated to the cognate bacteriome in alpha diversity, and undergone constantly contracting and shifting across the temporal scale. Nevertheless, 257 core virial clusters (VCs) persistently occurred with time whatever the significant impacts imposed by the varied physiochemical properties. Glycoside hydrolase (GH) and glycosyltransferase (GT) families genes displayed the higher abundances across all samples. Intriguingly, diversely structuring of toxin-antitoxin systems (TAs) and CRISPR-Cas arrays were frequently harbored by phage genomes. Their divergent organization and encoding attributes underlie the multiple biological roles in modulation of network and/or contest of phage community as well as bacterial host community. This phageome-wide mapping will fuel the current insights of phage community ecology in other traditional fermented ecosystems that are challenging to decipher.

Distribution patterns and functional diversity of DNA viruses determined by ecological niches in huge river ecosystems

While the vast number of DNA and RNA viruses participate in biogeochemical cycles in natural systems, little is known about virome in river ecosystems. Here, we analyzed the DNA viral composition and its metabolic potential in the Yangtze River, including freshwater (FW) and freshwater sediments (FWS). A total of 1237 river-derived virus contigs (RVCs) were obtained following de novo assembly from 62 metagenomics. We found that the viral diversity is significantly positively correlated longitudinally. Moreover, FW exhibited a greater viral variety and significantly different composition than FWS. The viral co-occurrence network suggested that positive correlations predominate between RVCs. Lastly, 1657 viral functions were predicted by gene ontology. Notably, 96 of 150 RVCs with higher weights identified by random-forest classier were more abundant in FW, which most engage organic cyclic compound metabolic processes and hydrolase activity. Together, this study highlights the previously unrecognized viruses and the importance of their distributions and functions in major river systems.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Community Structure, Drivers, and Potential Functions of Different Lifestyle Viruses in Chaohu Lake

Viruses, as the most prolific entities on Earth, constitute significant ecological groups within freshwater lakes, exerting pivotal ecological roles. In this study, we selected Chaohu Lake, a representative eutrophic freshwater lake in China, as our research site to explore the community distribution, driving mechanisms, and potential ecological functions of diverse viral communities, the intricate virus-host interaction systems, and the overarching influence of viruses on global biogeochemical cycling.

Temporal Dynamics and Contribution of Phage Community to the Prevalence of Antibiotic Resistance Genes in a Full-Scale Sludge Anaerobic Digestion Plant

Anaerobic digestion (AD) is an important biological resource recovery process, where microorganisms play key roles for material transformation. There has been some knowledge about the prokaryotic community and antibiotic resistance genes (ARGs) in AD, but there has been very limited knowledge of phages. In this study, samples from a full-scale AD plant were collected over 13 months, sequenced, and analyzed for viral and prokaryotic metagenomes. Totally, 3015 viral operational taxonomic units (vOTUs) were detected, mostly assigned to Caudoviricetes. The phage community had faster temporal variation than the prokaryotic community. Warm seasons harbored a higher abundance of both temperate phages and broad host-range phages. Seven ARGs of 6 subtypes were carried by 20 vOTUs, a representative ermT gene was synthesized and expressed, and the resistance activity in the host was examined, confirming the real activity of virus-carried ARGs in the AD process. Some of the ARGs were horizontally transferred between the phage and prokaryotic genomes. However, phage infection was not found to contribute to ARG transfer. This study provided an insight into the ecological patterns of the phage community, confirmed the antibiotic resistance activity of virus-carried ARGs, evaluated the contribution of phages on the ARG prevalence, and laid the foundation for the control strategies of the community and antibiotic resistance in the AD process.

Datascape: exploring heterogeneous dataspace

Data science is a powerful field for gaining insights, comparing, and predicting behaviors from datasets. However, the diversity of methods and hypotheses needed to abstract a dataset exhibits a lack of genericity. Moreover, the shape of a dataset, which structures its contained information and uncertainties, is rarely considered. Inspired by state-of-the-art manifold learning and hull estimations algorithms, we propose a novel framework, the datascape, that leverages topology and graph theory to abstract heterogeneous datasets. Built upon the combination of a nearest neighbor graph, a set of convex hulls, and a metric distance that respects the shape of the data, the datascape allows exploration of the dataset's underlying space. We show that the datascape can uncover underlying functions from simulated datasets, build predictive algorithms with performance close to state-of-the-art algorithms, and reveal insightful geodesic paths between points. It demonstrates versatility through ecological, medical, and simulated data use cases.

Discovery and description of novel phage genomes from urban microbiomes sampled by the MetaSUB consortium

Bacteriophages are recognized as the most abundant members of microbiomes and have therefore a profound impact on microbial communities through the interactions with their bacterial hosts. The International Metagenomics and Metadesign of Subways and Urban Biomes Consortium (MetaSUB) has sampled mass-transit systems in 60 cities over 3 years using metagenomics, throwing light into these hitherto largely unexplored urban environments. MetaSUB focused primarily on the bacterial community. In this work, we explored MetaSUB metagenomic data in order to recover and analyze bacteriophage genomes. We recovered and analyzed 1714 phage genomes with size at least 40 kbp, from the class Caudoviricetes, the vast majority of which (80%) are novel. The recovered genomes were predicted to belong to temperate (69%) and lytic (31%) phages. Thirty-three of these genomes have more than 200 kbp, and one of them reaches 572 kbp, placing it among the largest phage genomes ever found. In general, the phages tended to be site-specific or nearly so, but 194 genomes could be identified in every city from which phage genomes were retrieved. We predicted hosts for 48% of the phages and observed general agreement between phage abundance and the respective bacterial host abundance, which include the most common nosocomial multidrug-resistant pathogens. A small fraction of the phage genomes are carriers of antibiotic resistance genes, and such genomes tended to be particularly abundant in the sites where they were found. We also detected CRISPR-Cas systems in five phage genomes. This study expands the previously reported MetaSUB results and is a contribution to the knowledge about phage diversity, global distribution, and phage genome content.

Genome sequences of the first Autographiviridae phages infecting marine Roseobacter

The ubiquitous and abundant marine phages play critical roles in shaping the composition and function of bacterial communities, impacting biogeochemical cycling in marine ecosystems. Autographiviridae is among the most abundant and ubiquitous phage families in the ocean. However, studies on the diversity and ecology of Autographiviridae phages in marine environments are restricted to isolates that infect SAR11 bacteria and cyanobacteria. In this study, ten new roseophages that infect marine Roseobacter strains were isolated from coastal waters. These new roseophages have a genome size ranging from 38 917 to 42 634 bp and G+C content of 44.6-50 %. Comparative genomics showed that they are similar to known Autographiviridae phages regarding gene content and architecture, thus representing the first Autographiviridae roseophages. Phylogenomic analysis based on concatenated conserved genes showed that the ten roseophages form three distinct subgroups within the Autographiviridae, and sequence analysis revealed that they belong to eight new genera. Finally, viromic read-mapping showed that these new Autographiviridae phages are widely distributed in global oceans, mostly inhabiting polar and estuarine locations. This study has expanded the current understanding of the genomic diversity, evolution and ecology of Autographiviridae phages and roseophages. We suggest that Autographiviridae phages play important roles in the mortality and community structure of roseobacters, and have broad ecological applications.

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.

Benchmarking informatics approaches for virus discovery: caution is needed when combining in silico identification methods

Understanding the ecological impacts of viruses on natural and engineered ecosystems relies on the accurate identification of viral sequences from community sequencing data. To maximize viral recovery from metagenomes, researchers frequently combine viral identification tools. However, the effectiveness of this strategy is unknown. Here, we benchmarked combinations of six widely used informatics tools for viral identification and analysis (VirSorter, VirSorter2, VIBRANT, DeepVirFinder, CheckV, and Kaiju), called "rulesets." Rulesets were tested against mock metagenomes composed of taxonomically diverse sequence types and diverse aquatic metagenomes to assess the effects of the degree of viral enrichment and habitat on tool performance. We found that six rulesets achieved equivalent accuracy [Matthews Correlation Coefficient (MCC) = 0.77, Padj ≥ 0.05]. Each contained VirSorter2, and five used our "tuning removal" rule designed to remove non-viral contamination. While DeepVirFinder, VIBRANT, and VirSorter were each found once in these high-accuracy rulesets, they were not found in combination with each other: combining tools does not lead to optimal performance. Our validation suggests that the MCC plateau at 0.77 is partly caused by inaccurate labeling within reference sequence databases. In aquatic metagenomes, our highest MCC ruleset identified more viral sequences in virus-enriched (44%-46%) than in cellular metagenomes (7%-19%). While improved algorithms may lead to more accurate viral identification tools, this should be done in tandem with careful curation of sequence databases. We recommend using the VirSorter2 ruleset and our empirically derived tuning removal rule. Our analysis provides insight into methods for in silico viral identification and will enable more robust viral identification from metagenomic data sets.

IMPORTANCE

The identification of viruses from environmental metagenomes using informatics tools has offered critical insights in microbial ecology. However, it remains difficult for researchers to know which tools optimize viral recovery for their specific study. In an attempt to recover more viruses, studies are increasingly combining the outputs from multiple tools without validating this approach. After benchmarking combinations of six viral identification tools against mock metagenomes and environmental samples, we found that these tools should only be combined cautiously. Two to four tool combinations maximized viral recovery and minimized non-viral contamination compared with either the single-tool or the five- to six-tool ones. By providing a rigorous overview of the behavior of in silico viral identification strategies and a pipeline to replicate our process, our findings guide the use of existing viral identification tools and offer a blueprint for feature engineering of new tools that will lead to higher-confidence viral discovery in microbiome studies.

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

INTRODUCTION

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSION

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

Virome profiling of fig wasps (Ceratosolen spp.) reveals virus diversity spanning four realms

We investigated the virome of agaonid fig wasps (Ceratosolen spp.) inside syconia ("fruits") of various Ficus trees fed upon by frugivores such as pteropodid bats in Sub-Saharan Africa. This virome includes representatives of viral families spanning four realms and includes near-complete genome sequences of three novel viruses and fragments of five additional potentially novel viruses evolutionarily associated with insects, fungi, plants, and vertebrates. Our study provides evidence that frugivorous animals are exposed to a plethora of viruses by coincidental consumption of fig wasps, which are obligate pollinators of figs worldwide.

Viral potential to modulate microbial methane metabolism varies by habitat

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

Identification of HDV-like theta ribozymes involved in tRNA-based recoding of gut bacteriophages

Trillions of microorganisms, collectively known as the microbiome, inhabit our bodies with the gut microbiome being of particular interest in biomedical research. Bacteriophages, the dominant virome constituents, can utilize suppressor tRNAs to switch to alternative genetic codes (e.g., the UAG stop-codon is reassigned to glutamine) while infecting hosts with the standard bacterial code. However, what triggers this switch and how the bacteriophage manipulates its host is poorly understood. Here, we report the discovery of a subgroup of minimal hepatitis delta virus (HDV)-like ribozymes - theta ribozymes - potentially involved in the code switch leading to the expression of recoded lysis and structural phage genes. We demonstrate their HDV-like self-scission behavior in vitro and find them in an unreported context often located with their cleavage site adjacent to tRNAs, indicating a role in viral tRNA maturation and/or regulation. Every fifth associated tRNA is a suppressor tRNA, further strengthening our hypothesis. The vast abundance of tRNA-associated theta ribozymes - we provide 1753 unique examples - highlights the importance of small ribozymes as an alternative to large enzymes that usually process tRNA 3'-ends. Our discovery expands the short list of biological functions of small HDV-like ribozymes and introduces a previously unknown player likely involved in the code switch of certain recoded gut bacteriophages.

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

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

Shewanella phage encoding a putative anti-CRISPR-like gene represents a novel potential viral family

Shewanella is a prevalent bacterial genus in deep-sea environments including marine sediments, exhibiting diverse metabolic capabilities that indicate its significant contributions to the marine biogeochemical cycles. However, only a few Shewanella phages were isolated and deposited in the NCBI database. In this study, we report the isolation and characterization of a novel Shewanella phage, vB_SbaS_Y11, that infects Shewanella KR11 and was isolated from the sewage in Qingdao, China. Transmission electron microscopy revealed that vB_SbaS_Y11 has an icosahedral head and a long tail. The genome of vB_SbaS_Y11 is a linear, double-stranded DNA with a length of 62,799 bp and a G+C content of 46.9%, encoding 71 putative open reading frames. No tRNA genes or integrase-related feature genes were identified. An uncharacterized anti-CRISPR AcrVA2 gene was detected in its genome. Phylogenetic analysis based on the amino acid sequences of whole genomes and comparative genomic analyses indicate that vB_SbaS_Y11 has a novel genomic architecture and shares low similarity to Pseudomonas virus H66 and Pseudomonas phage F116. vB_SbaS_Y11 represents a potential new family-level virus cluster with eight metagenomic assembled viral genomes named Ranviridae.IMPORTANCEThe Gram-negative Shewanella bacterial genus currently includes about 80 species of mostly aquatic Gammaproteobacteria, which were isolated around the globe in a multitude of environments, such as freshwater, seawater, coastal sediments, and the deepest trenches. Here, we present a Shewanella phage vB_SbaS_Y11 that contains an uncharacterized anti-CRISPR AcrVA2 gene and belongs to a potential virus family, Ranviridae. This study will enhance the knowledge about the genome, diversity, taxonomic classification, and global distribution of Shewanella phage populations.

Prophage enhances the ability of deep-sea bacterium Shewanella psychrophila WP2 to utilize D-amino acid

Prophages are prevalent in the marine bacterial genomes and reshape the physiology and metabolism of their hosts. However, whether and how prophages influence the microbial degradation of D-amino acids (D-AAs), which is one of the widely distributed recalcitrant dissolved organic matters (RDOMs) in the ocean, remain to be explored. In this study, we addressed this issue in a representative marine bacterium, Shewanella psychrophila WP2 (WP2), and its integrated prophage SP1. Notably, compared to the WP2 wild-type strain, the SP1 deletion mutant of WP2 (WP2ΔSP1) exhibited a significantly lower D-glutamate (D-Glu) consumption rate and longer lag phase when D-Glu was used as the sole nitrogen source. The subsequent transcriptome analysis identified 1,523 differentially expressed genes involved in diverse cellular processes, especially that multiple genes related to inorganic nitrogen metabolism were highly upregulated. In addition, the dynamic profiles of ammonium, nitrate, and nitrite were distinct between the culture media of WP2 and WP2ΔSP1. Finally, we provide evidence that SP1 conferred a competitive advantage to WP2 when D-Glu was used as the sole nitrogen source and SP1-like phages may be widely distributed in the global ocean. Taken together, these findings offer novel insight into the influences of prophages on host metabolism and RDOM cycling in marine environments.IMPORTANCEThis work represents the first exploration of the impact of prophages on the D-amino acid (D-AA) metabolism of deep-sea bacteria. By using S. psychrophila WP2 and its integrated prophage SP1 as a representative system, we found that SP1 can significantly increase the catabolism rate of WP2 to D-glutamate and produce higher concentrations of ammonium, resulting in faster growth and competitive advantages. Our findings not only deepen our understanding of the interaction between deep-sea prophages and hosts but also provide new insights into the ecological role of prophages in refractory dissolved organic matter and the nitrogen cycle in deep oceans.

The structure and assembly mechanisms of T4-like cyanophages community in the South China Sea

Marine ecosystems contain an immense diversity of phages, many of which infect cyanobacteria (cyanophage) that are largely responsible for primary productivity. To characterize the genetic diversity and biogeographic distribution of the marine T4-like cyanophage community in the northern South China Sea, the T4-like cyanophage portal protein gene (g20) was amplified. Phylogenetic analysis revealed that marine T4-like cyanophages were highly diverse, with g20 operational taxonomic units being affiliated with five defined clades (Clusters I-V). Cluster II had a wide geographic distribution, Cluster IV was the most abundant in the open sea, and Cluster I was dominant in coastal shelf environments. Our results showed T4-like cyanophages (based on g20) community was generally shaped via heterogeneous selection. Highly variable environmental factors (such as salinity and temperature) can heterogeneously select different cyanophage communities. Nevertheless, the dominant drivers of the T4-like cyanophage community based on the g20 and g23 (T4-like phage major capsid protein gene) were different, probably due to different coverages by the primer sets. Furthermore, the community assembly processes of T4-like cyanophages were affected by host traits (abundance and distribution), viral traits (latent period, burst size, and host range), and environmental properties (temperature and salinity).IMPORTANCECyanophages are abundant and ubiquitous in the oceans, altering population structures and evolution of cyanobacteria, which account for a large portion of global carbon fixation, through host mortality, horizontal gene transfer, and the modulation of host metabolism. However, little is known about the biogeography and ecological drivers that shape the cyanophage community. Here, we use g20 and g23 genes to examine the biogeographic patterns and the assembly mechanisms of T4-like cyanophage community in the northern part of the South China Sea. The different coverages of primer sets might lead to the different dominant drivers of T4-like cyanophage community based on g20 and g23 genes. Our results showed that characteristics of viral traits (latent period, burst size, and host range) and host traits (abundance and distribution) were found to either limit or enhance the biogeographic distribution of T4-like cyanophages. Overall, both virus and host properties are critical to consider when determining rules of community assembly for viruses.

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.

Comparative evaluation of bioinformatic tools for virus-host prediction and their application to a highly diverse community in the Cuatro Ciénegas Basin, Mexico

Due to the enormous diversity of non-culturable viruses, new viruses must be characterized using culture-independent techniques. The associated host is an important phenotypic feature that can be inferred from metagenomic viral contigs thanks to the development of several bioinformatic tools. Here, we compare the performance of recently developed virus-host prediction tools on a dataset of 1,046 virus-host pairs and then apply the best-performing tools to a metagenomic dataset derived from a highly diverse transiently hypersaline site known as the Archaean Domes (AD) within the Cuatro Ciénegas Basin, Coahuila, Mexico. Among host-dependent methods, alignment-based approaches had a precision of 66.07% and a sensitivity of 24.76%, while alignment-free methods had an average precision of 75.7% and a sensitivity of 57.5%. RaFAH, a virus-dependent alignment-based tool, had the best overall performance (F1_score = 95.7%). However, when predicting the host of AD viruses, methods based on public reference databases (such as RaFAH) showed lower inter-method agreement than host-dependent methods run against custom databases constructed from prokaryotes inhabiting AD. Methods based on custom databases also showed the greatest agreement between the source environment and the predicted host taxonomy, habitat, lifestyle, or metabolism. This highlights the value of including custom data when predicting hosts on a highly diverse metagenomic dataset, and suggests that using a combination of methods and qualitative validations related to the source environment and predicted host biology can increase the number of correct predictions. Finally, these predictions suggest that AD viruses infect halophilic archaea as well as a variety of bacteria that may be halophilic, halotolerant, alkaliphilic, thermophilic, oligotrophic, sulfate-reducing, or marine, which is consistent with the specific environment and the known geological and biological evolution of the Cuatro Ciénegas Basin and its microorganisms.

Environmental viromes reveal global virosphere of deep-sea sediment RNA viruses

INTRODUCTION

Viruses are the most abundant and diverse life forms on the earth. Both DNA viruses and RNA viruses play important roles in marine ecosystems via regulating biogeochemical cycles.

OBJECTIVES

However, the virome of marine RNA viruses has been rarely explored so far. In this study, therefore, the environmental viromes of deep-sea sediment RNA viruses were characterized on a global scale to reveal the global virosphere of deep-sea RNA viruses.

METHODS

The viral particles were purified from each of 133 deep-sea sediment samples and then characterized based on metagenomes of RNA viruses.

RESULTS

In this study, we established the global virome dataset of deep-sea RNA viruses purified from 133 sediment samples that were collected from typical deep-sea ecosystems of three oceans. A total of 85,059 viral operational taxonomic units (vOTUs) were identified, of which only 1.72% were hitherto known, indicating that the deep-sea sediment is a repository of novel RNA viruses. These vOTUs were classified into 20 viral families, including prokaryotic (7.09%) and eukaryotic (65.81%) RNA viruses. Furthermore, 1,463 deep-sea RNA viruses with complete genomes were obtained. The differentiation of RNA viral communities was driven by the deep-sea ecosystems as opposed to geographical region. Specifically, the virus-encoded metabolic genes took great effects on the differentiation of RNA viral communities by mediating the energy metabolism in the deep-sea ecosystems.

CONCLUSIONS

Therefore, our findings indicate that the deep sea is a vast reservoir of novel RNA viruses for the first time, and the differentiation of RNA viral communities is driven by the deep-sea ecosystems through energy metabolism.

The long and short of it: benchmarking viromics using Illumina, Nanopore and PacBio sequencing technologies

Viral metagenomics has fuelled a rapid change in our understanding of global viral diversity and ecology. Long-read sequencing and hybrid assembly approaches that combine long- and short-read technologies are now being widely implemented in bacterial genomics and metagenomics. However, the use of long-read sequencing to investigate viral communities is still in its infancy. While Nanopore and PacBio technologies have been applied to viral metagenomics, it is not known to what extent different technologies will impact the reconstruction of the viral community. Thus, we constructed a mock bacteriophage community of previously sequenced phage genomes and sequenced them using Illumina, Nanopore and PacBio sequencing technologies and tested a number of different assembly approaches. When using a single sequencing technology, Illumina assemblies were the best at recovering phage genomes. Nanopore- and PacBio-only assemblies performed poorly in comparison to Illumina in both genome recovery and error rates, which both varied with the assembler used. The best Nanopore assembly had errors that manifested as SNPs and INDELs at frequencies 41 and 157 % higher than found in Illumina only assemblies, respectively. While the best PacBio assemblies had SNPs at frequencies 12 and 78 % higher than found in Illumina-only assemblies, respectively. Despite high-read coverage, long-read-only assemblies recovered a maximum of one complete genome from any assembly, unless reads were down-sampled prior to assembly. Overall the best approach was assembly by a combination of Illumina and Nanopore reads, which reduced error rates to levels comparable with short-read-only assemblies. When using a single technology, Illumina only was the best approach. The differences in genome recovery and error rates between technology and assembler had downstream impacts on gene prediction, viral prediction, and subsequent estimates of diversity within a sample. These findings will provide a starting point for others in the choice of reads and assembly algorithms for the analysis of viromes.

Structural changes in the gut virome of patients with atherosclerotic cardiovascular disease

IMPORTANCE

Existing studies have found that there is a close relationship between human virome and numerous diseases, and diseases may affect the diversity and composition of the virome; at the same time, changes in the virome will in turn affect the onset and progression of the disease. However, the composition and functional capabilities of the gut virome associated with atherosclerotic cardiovascular disease (ACVD) have not been systematically investigated. To our knowledge, this is the first study investigating the gut virome in patients with ACVD. We characterized the structural changes in the gut virome of ACVD patients, which may facilitate additional mechanistic, diagnostic, and interventional studies of ACVD and related diseases.

Diversities and interactions of phages and bacteria in deep-sea sediments as revealed by metagenomics

Phages are found virtually everywhere, even in extreme environments, and are extremely diverse both in their virion structures and in their genomic content. They are thought to shape the taxonomic and functional composition of microbial communities as well as their stability. A number of studies on laboratory culture and viral metagenomic research provide deeper insights into the abundance, diversity, distribution, and interaction with hosts of phages across a wide range of ecosystems. Although most of these studies focus on easily accessible samples, such as soils, lakes, and shallow oceans, little is known about bathypelagic phages. In this study, through analyzing the 16S rRNA sequencing and viral metagenomic sequencing data of 25 samples collected from five different bathypelagic ecosystems, we detected a high diversity of bacteria and phages, particularly in the cold seep and hydrothermal vent ecosystems, which have stable chemical energy. The relative abundance of phages in these ecosystems was higher than in other three abyssal ecosystems. The low phage/host ratios obtained from host prediction were different from shallow ecosystems and indicated the prevalence of prophages, suggesting the complexity of phage-bacteria interactions in abyssal ecosystems. In the correlation analysis, we revealed several phages-bacteria interaction networks of potential ecological relevance. Our study contributes to a better understanding of the interactions between bathypelagic bacteria and their phages.

High niche specificity and host genetic diversity of groundwater viruses

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

Megataxonomy and global ecology of the virosphere

Nearly all organisms are hosts to multiple viruses that collectively appear to be the most abundant biological entities in the biosphere. With recent advances in metagenomics and metatranscriptomics, the known diversity of viruses substantially expanded. Comparative analysis of these viruses using advanced computational methods culminated in the reconstruction of the evolution of major groups of viruses and enabled the construction of a virus megataxonomy, which has been formally adopted by the International Committee on Taxonomy of Viruses. This comprehensive taxonomy consists of six virus realms, which are aspired to be monophyletic and assembled based on the conservation of hallmark proteins involved in capsid structure formation or genome replication. The viruses in different major taxa substantially differ in host range and accordingly in ecological niches. In this review article, we outline the latest developments in virus megataxonomy and the recent discoveries that will likely lead to reassessment of some major taxa, in particular, split of three of the current six realms into two or more independent realms. We then discuss the correspondence between virus taxonomy and the distribution of viruses among hosts and ecological niches, as well as the abundance of viruses versus cells in different habitats. The distribution of viruses across environments appears to be primarily determined by the host ranges, i.e. the virome is shaped by the composition of the biome in a given habitat, which itself is affected by abiotic factors.

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.

PhageScope: a well-annotated bacteriophage database with automatic analyses and visualizations

Bacteriophages are viruses that infect bacteria or archaea. Understanding the diverse and intricate genomic architectures of phages is essential to study microbial ecosystems and develop phage therapy strategies. However, the existing phage databases are short of meticulous annotations. To this end, we propose PhageScope (https://phagescope.deepomics.org), an online phage database with comprehensive annotations. PhageScope harbors a collection of 873 718 phage sequences from various sources. Applying fifteen state-of-the-art tools to perform systematic annotations and analyses, PhageScope provides annotations on genome completeness, host range, lifestyle information, taxonomy classification, nine types of structural and functional genetic elements, and three types of comparative genomic studies for curated phages. Additionally, PhageScope incorporates automatic analyses and visualizations for curated and customized phages, serving as an efficient platform for phage study.

IPEV: identification of prokaryotic and eukaryotic virus-derived sequences in virome using deep learning

BACKGROUND

The virome obtained through virus-like particle enrichment contains a mixture of prokaryotic and eukaryotic virus-derived fragments. Accurate identification and classification of these elements are crucial to understanding their roles and functions in microbial communities. However, the rapid mutation rates of viral genomes pose challenges in developing high-performance tools for classification, potentially limiting downstream analyses.

FINDINGS

We present IPEV, a novel method to distinguish prokaryotic and eukaryotic viruses in viromes, with a 2-dimensional convolutional neural network combining trinucleotide pair relative distance and frequency. Cross-validation assessments of IPEV demonstrate its state-of-the-art precision, significantly improving the F1-score by approximately 22% on an independent test set compared to existing methods when query viruses share less than 30% sequence similarity with known viruses. Furthermore, IPEV outperforms other methods in accuracy on marine and gut virome samples based on annotations by sequence alignments. IPEV reduces runtime by at most 1,225 times compared to existing methods under the same computing configuration. We also utilized IPEV to analyze longitudinal samples and found that the gut virome exhibits a higher degree of temporal stability than previously observed in persistent personal viromes, providing novel insights into the resilience of the gut virome in individuals.

CONCLUSIONS

IPEV is a high-performance, user-friendly tool that assists biologists in identifying and classifying prokaryotic and eukaryotic viruses within viromes. The tool is available at https://github.com/basehc/IPEV.

Vibrio cyclitrophicus phage encoding gene transfer agent fragment, representing a novel viral family

Vibrio is a prevalent bacterial genus in aquatic environments and exhibits diverse metabolic capabilities, playing a vital role in marine biogeochemical cycles. This study isolated a novel virus infecting Vibrio cyclitrophicus, vB_VviC_ZQ26, from coastal waters near Qingdao, China. The vB_VviC_ZQ26 comprises a linear double-stranded DNA genome with a length of 42,982 bp and a G + C content of 43.21 %, encoding 72 putative open reading frames (ORFs). Transmission electron microscope characterization indicates a siphoviral-morphology of vB_VviC_ZQ26. Nucleic-acids-wide analysis indicates a tetranucleotide frequency deviation for genomic segments encoding putative gene transfer agent protein (GTA) and coil-containing protein, implying divergent origins occurred in different parts of viral genomes. Phylogenetic and genome-content-based analysis suggest that vB_VviC_ZQ26 represents a novel vibriophage-specific family designated as Coheviridae. From the result of biogeographic analysis, Coheviridae is mainly colonized in the temperate and tropical epipelagic zones. This study describes a novel vibriophage infecting V. cyclitrophicus, shedding light on the evolutionary divergence of different parts of the viral genome and its ecological footprint in marine environments.

Distinct and rich assemblages of giant viruses in Arctic and Antarctic lakes

Giant viruses (GVs) are key players in ecosystem functioning, biogeochemistry, and eukaryotic genome evolution. GV diversity and abundance in aquatic systems can exceed that of prokaryotes, but their diversity and ecology in lakes, especially polar ones, remain poorly understood. We conducted a comprehensive survey and meta-analysis of GV diversity across 20 lakes, spanning polar to temperate regions, combining our extensive lake metagenome database from the Canadian Arctic and subarctic with publicly available datasets. Leveraging a novel GV genome identification tool, we identified 3304 GV metagenome-assembled genomes, revealing lakes as untapped GV reservoirs. Phylogenomic analysis highlighted their dispersion across all Nucleocytoviricota orders. Strong GV population endemism emerged between lakes from similar regions and biomes (Antarctic and Arctic), but a polar/temperate barrier in lacustrine GV populations and differences in their gene content could be observed. Our study establishes a robust genomic reference for future investigations into lacustrine GV ecology in fast changing polar environments.

Breaking the Ice: A Review of Phages in Polar Ecosystems

Bacteriophages, or phages, are viruses that infect and replicate within bacterial hosts, playing a significant role in regulating microbial populations and ecosystem dynamics. However, phages from extreme environments such as polar regions remain relatively understudied due to challenges such as restricted ecosystem access and low biomass. Understanding the diversity, structure, and functions of polar phages is crucial for advancing our knowledge of the microbial ecology and biogeochemistry of these environments. In this review, we will explore the current state of knowledge on phages from the Arctic and Antarctic, focusing on insights gained from -omic studies, phage isolation, and virus-like particle abundance data. Metagenomic studies of polar environments have revealed a high diversity of phages with unique genetic characteristics, providing insights into their evolutionary and ecological roles. Phage isolation studies have identified novel phage-host interactions and contributed to the discovery of new phage species. Virus-like particle abundance and lysis rate data, on the other hand, have highlighted the importance of phages in regulating bacterial populations and nutrient cycling in polar environments. Overall, this review aims to provide a comprehensive overview of the current state of knowledge about polar phages, and by synthesizing these different sources of information, we can better understand the diversity, dynamics, and functions of polar phages in the context of ongoing climate change, which will help to predict how polar ecosystems and residing phages may respond to future environmental perturbations.

Bacteriophage Taxonomy: A Continually Evolving Discipline

While taxonomy is an often underappreciated branch of science, it serves very important roles. Bacteriophage taxonomy has evolved from a discipline based mainly on morphology, characterized by the work of David Bradley and Hans-Wolfgang Ackermann, to the sequence-based approach that is taken today. The Bacterial Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) takes a holistic approach to classifying prokaryote viruses by measuring overall DNA and protein similarity and phylogeny before making decisions about the taxonomic position of a new virus. The huge number of complete genomes being deposited with the National Center for Biotechnology Information (NCBI) and other public databases has resulted in a reassessment of the taxonomy of many viruses, and the future will see the introduction of new viral families and higher orders.

Driving through stop signs: predicting stop codon reassignment improves functional annotation of bacteriophages

The majority of bacteriophage diversity remains uncharacterized, and new intriguing mechanisms of their biology are being continually described. Members of some phage lineages, such as the Crassvirales, repurpose stop codons to encode an amino acid by using alternate genetic codes. Here, we investigated the prevalence of stop codon reassignment in phage genomes and its subsequent impacts on functional annotation. We predicted 76 genomes within INPHARED and 712 vOTUs from the Unified Human Gut Virome Catalogue (UHGV) that repurpose a stop codon to encode an amino acid. We re-annotated these sequences with modified versions of Pharokka and Prokka, called Pharokka-gv and Prokka-gv, to automatically predict stop codon reassignment prior to annotation. Both tools significantly improved the quality of annotations, with Pharokka-gv performing best. For sequences predicted to repurpose TAG to glutamine (translation table 15), Pharokka-gv increased the median gene length (median of per genome median) from 287 to 481 bp for UHGV sequences (67.8% increase) and from 318 to 550 bp for INPHARED sequences (72.9% increase). The re-annotation increased median coding capacity from 66.8% to 90.0% and from 69.0% to 89.8% for UHGV and INPHARED sequences predicted to use translation table 15. Furthermore, the proportion of genes that could be assigned functional annotation increased, including an increase in the number of major capsid proteins that could be identified. We propose that automatic prediction of stop codon reassignment before annotation is beneficial to downstream viral genomic and metagenomic analyses.

Host species and geography impact bee-associated RNA virus communities with evidence for isolation by distance in viral populations

Virus symbionts are important mediators of ecosystem function, yet we know little of their diversity and ecology in natural populations. The alarming decline of pollinating insects in many regions of the globe, especially the European honey bee, Apis mellifera, has been driven in part by worldwide transmission of virus pathogens. Previous work has examined the transmission of known honey bee virus pathogens to wild bee populations, but only a handful of studies have investigated the native viromes associated with wild bees, limiting epidemiological predictors associated with viral pathogenesis. Further, variation among different bee species might have important consequences in the acquisition and maintenance of bee-associated virome diversity. We utilized comparative metatranscriptomics to develop a baseline description of the RNA viromes associated with wild bee pollinators and to document viral diversity, community composition, and structure. Our sampling includes five wild-caught, native bee species that vary in social behavior as well as managed honey bees. We describe 26 putatively new RNA virus species based on RNA-dependent RNA polymerase phylogeny and show that each sampled bee species was associated with a specific virus community composition, even among sympatric populations of distinct host species. From 17 samples of a single host species, we recovered a single virus species despite over 600 km of distance between host populations and found strong evidence for isolation by distance in associated viral populations. Our work adds to the small number of studies examining viral prevalence and community composition in wild bees.

Quantification of Virion-Sense and Complementary-Sense DNA Strands of Circular Single-Stranded DNA Viruses

Circular ssDNA viruses are ubiquitous and can be found in both prokaryotes and eukaryotes. To understand the interaction of ssDNA viruses with their hosts, it is important to characterize the dynamics of viral sense (VS) and complementary-sense (CS) viral strands during the infection process. Here, we present a simple and rapid protocol that allows sensitive and accurate determination of the VS and CS strands generated during viral infection.The method consists of a two-step qPCR in which the first step uses a strand-specific (CS or VS) labeled primer and T4 DNA polymerase that lacks strand displacement activity and makes a single copy per VS or CS strand. Next, the T4 DNA polymerase and unincorporated oligonucleotides are removed by a silica membrane spin column. Finally, the purified VS or CS strands are quantified by qPCR in a second step in which amplification uses a tag primer and a specific primer. Absolute quantification of VS and CS strands is obtained by extrapolating the Cq data to a standard curve of ssDNA, which can be generated by phagemid expression. Quantification of VS and CS strands of two geminiviruses in infections of Solanum lycopersicum (tomato) and Nicotiana benthamiana plants using this method is shown.

Structural and Genomic Diversity of Bacteriophages

Bacteriophage diversity is a relatively unknown frontier that is rapidly being explored, leading to a wealth of new information. New bacteriophages are being discovered at an astounding rate via both phage isolation studies and metagenomic analyses. In addition, a nucleotide sequence-based viral taxonomic system has been developed to better handle this wealth of new information. As a result of these developments, phage scientists are transitioning from knowing that there must be huge numbers of diverse kinds of phage particles in natural environments to identifying the actual abundance and phage diversity that is present in specific environments. This review documents the beginning of this transition, offering a glimpse into the magnitude of change unfolding in the field. It stands as a testament to the expanding frontiers of phage research, illuminating the remarkable progress made in unraveling the intricate world of bacteriophage diversity and advancing our understanding of these enigmatic viral entities.