Viral dark matter and virus-host interactions resolved from publicly available microbial genomes

Healthy and moribund Zhikong scallops (Chlamys farreri) developed different viral communities during a mass mortality event

Viral assemblages of scallops are still relatively unknown. Here, metagenomic analysis was used to study virus communities in the gut of scallops to establish the first scallop virome data set (SVD); this contains 7,447 viral operational taxonomic units. Protein-sharing networks and phylogenetic analyses demonstrated the high diversity and novelty of the SVD, which is very different from viromes from other habitats. Potentially pathogenic viruses are prevalent in the gut of scallops. In particular, the novel smacoviruses were identified, indicating that scallops may be a potential hotspot for this viral group. Inference of virus-host associations found extensive interactions between viruses and major prokaryotic lineages. Intriguingly, moribund scallops showed a higher diversity of auxiliary metabolic genes (AMGs) related to amino acid metabolism and cofactor and vitamin genes, while healthy scallops had fewer AMGs, with those present focusing on secondary metabolite biosynthesis and carbohydrate metabolism. These findings provide the first landscape of scallop gut viruses based on metagenomes and highlight the potential roles of diverse and unique gut viruses for the health of filter-feeding bivalves.

IMPORTANCE

This study uses metagenome sequencing to establish the first scallop virome database. The study reveals previously unknown diversity of scallop-associated viruses and provides insights into links between disease status and viral diversity and genome content. The study will interest many aquatic virologists and could have important implications in managing marine resources.

Unveiling the dynamic viral landscape across developmental stages of cold seep ecosystems: Implications for global marine biogeochemistry

Cold seeps are methane-rich ecosystems playing pivotal roles in global biogeochemical cycling, yet their viral communities remain underexplored. We present the first comprehensive viral metagenomic analysis across developmental stages of the Haima Cold Seep. Characterizing viral assemblages from chemoautotrophic, mature, and extinct seep sediments revealed 4272 viral operational taxonomic units, with 77 % representing novel lineages, highlighting cold seeps' unique viral diversity. Viral community structure and diversity varied significantly by seep stage, with highest diversity in the chemoautotrophic stage. While Siphoviridae and Microviridae dominated, their relative abundances shifted with maturity. Gammaproteobacteria emerged as predominant viral hosts, exhibiting distinct interaction patterns across stages. Notably, the chemoautotrophic stage harbored the highest abundance and diversity of virus-encoded auxiliary metabolic genes (AMGs; ∼450 AMGs), with significantly enriched carbohydrate metabolism and central carbon pathway genes (2.2-fold and 1.8-fold higher respectively, p < 0.005), amino acid metabolism (1.9-fold, p = 0.003), and sulfur relay system genes (2.0-fold, p = 0.002). In contrast, the mature stage exhibited distinct enrichment in energy metabolism genes (up to 3.9-fold difference between sites, p < 0.001) and xenobiotics degradation pathways, suggesting stage-specific viral impacts on biogeochemical cycling. Lytic lifestyles prevailed across stages, indicating dynamic virus-host interactions during seep development. These findings unveil complex viral ecology in cold seeps, with potential influences on microbial community structure and biogeochemical processes. Providing a foundation for understanding viral roles in cold seep ecosystem functioning and biogeochemical cycles, this study has implications for marine microbial ecology and environmental biotechnology.

Exploring Mobile Genetic Elements in Vibrio cholerae

Members of the bacterial species Vibrio cholerae are known both as prominent constituents of marine environments and as the causative agents of cholera, a severe diarrheal disease. While strains responsible for cholera have been extensively studied over the past century, less is known about their environmental counterparts, despite their contributions to the species' pangenome. This study analyzed the genome compositions of 46 V. cholerae strains, including pandemic and nonpandemic, toxigenic, and environmental variants, to investigate the diversity of mobile genetic elements (MGEs), embedded bacterial defense systems, and phage-associated signatures. Our findings include both conserved and novel MGEs across strains, pointing to shared evolutionary pathways and ecological niches. The defensome analysis revealed a wide array of antiphage/antiplasmid mechanisms, extending well beyond the traditional CRISPR-Cas and restriction-modification systems. This underscores the dynamic arms race between V. cholerae and MGEs and suggests that nonpandemic strains may act as reservoirs for emerging defense strategies. Moreover, the study showed that MGEs are integrated into genomic hotspots, which may serve as critical platforms for the exchange of defense systems, thereby enhancing V. cholerae's adaptive capabilities against phage attacks and other invading MGEs. Overall, this research offers new insights into V. cholerae's genetic complexity and potential adaptive strategies, offering a better understanding of the differences between environmental strains and their pandemic counterparts, as well as the possible evolutionary pathways that led to the emergence of pandemic strains.

Molecular piracy in deep-sea hydrothermal vent: phage-plasmid interactions revealed by phage-FISH in Marinitoga piezophila

Prokaryotes and mobile genetic elements (MGEs, such as viruses and plasmids) interact extensively, leading to horizontal gene transfer (HGT) and consequent microbial evolution and diversity. However, our knowledge of the interactions between MGEs in deep-sea hydrothermal ecosystems is limited. In this study, we adapted a phage-fluorescence in situ hybridization (phage-FISH) approach to visualize and quantify the dynamics of phage-plasmid interactions in an anaerobic, thermophilic deep-sea bacterium, Marinitoga piezophila. Notably, our results revealed that plasmid signals were detected in viral particles released from lysed cells, indicating that mitomycin C not only induced plasmid replication but also its packaging into phage particles. Further analysis of the DNA content in purified virions showed that the phage capsids incorporated plasmid DNA even without induction, and the majority of capsids (up to 70%) preferentially packaged plasmid DNA rather than viral DNA after induction. Therefore, this study provided direct evidence of molecular piracy in the deep-sea hydrothermal ecosystem, highlighting the important roles of selfish MGEs in virus-host interactions and HGT in extreme marine environments.

IMPORTANCE

Deep-sea hydrothermal vents are hotspots for microbes. Several studies revealed that virus-mediated horizontal gene transfer (HGT) in deep-sea hydrothermal vent ecosystems may be crucial to the survival and stability of prokaryotes in these extreme environments. However, little is known about the interaction between viruses and other mobile genetic elements (MGEs, such as plasmids), and how their interactions influence virus-mediated HGT in these ecosystems. In this study, we adapted a phage-fluorescence in situ hybridization approach to directly monitor the dynamics of phage-plasmid-host interactions at the single-cell level in the Marinitoga piezophila model. Interestingly, our results indicate that plasmid DNA could not only be induced by mitomycin C to a great extent but also hijacked viral assembly machinery to facilitate its propagation and spread. Therefore, the data presented here imply that the interaction between the viruses and other MGEs could play profound roles in virus-host interaction and virus-mediated HGT in the deep-sea hydrothermal ecosystem.

Examining the healthy human microbiome concept

Human microbiomes are essential to health throughout the lifespan and are increasingly recognized and studied for their roles in metabolic, immunological and neurological processes. Although the full complexity of these microbial communities is not fully understood, their clinical and industrial exploitation is well advanced and expanding, needing greater oversight guided by a consensus from the research community. One of the most controversial issues in microbiome research is the definition of a 'healthy' human microbiome. This concept is complicated by the microbial variability over different spatial and temporal scales along with the challenge of applying a unified definition to the spectrum of healthy microbiome configurations. In this Perspective, we examine the progress made and the key gaps that remain to be addressed to fully harness the benefits of the human microbiome. We propose a road map to expand our knowledge of the microbiome-health relationship, incorporating epidemiological approaches informed by the unique ecological characteristics of these communities.

A novel roseosiphovirus infecting dinoroseobacter shibae DFL12T represents a new genus

Bacteria belonging to the Roseobacter clade are key players in marine ecosystems, contributing significantly to carbon and sulfur cycles. Marine viruses, particularly those targeting Roseobacter, play crucial roles in regulating microbial communities and biogeochemical processes. Despite their importance, phages infecting organisms of the Roseobacter clade remain poorly understood. In this study, a novel roseophage, vB_DshS-R26L (R26L), infecting Dinoroseobacter shibae DFL12T, was isolated and characterized in terms of physiological and genomic properties. R26L has siphovirus morphology with an elongated head and a long, non-flexible tail. The phage has a narrow host range and demonstrates a long infection cycle with a latent period of 3.5 h and a burst size of 22 plaque-forming units (PFU cell- 1). R26L possesses a circular, double-stranded DNA genome of 79,534 bp with a G + C content of 62.6%, encoding a total of 116 open reading frames. Notably, seven auxiliary metabolic genes (AMGs), including those related to phosphate metabolism and queuosine biosynthesis, were identified. Phylogenetic and taxonomic analyses revealed that R26L represents a new genus, with its highest intergenomic similarities being 54.7% to another roseophage (R5C). By elucidating the unique characteristics of R26L, this study highlights the complexity of phage infections and the genomic diversity of roseophages, offering valuable insights into the ecological significance of Roseobacter-phage interactions in marine environments.

Are You My Host? An Overview of Methods Used to Link Bacteriophages with Hosts

Until recently, the only methods for finding out if a particular strain or species of bacteria could be a host for a particular bacteriophage was to see if the bacteriophage could infect that bacterium and kill it, releasing progeny phages. Establishing the host range of a bacteriophage thus meant infecting many different bacteria and seeing if the phage could kill each one. Detection of bacterial killing can be achieved on solid media (plaques, spots) or broth (culture clearing). More recently, additional methods to link phages and hosts have been developed. These include methods to show phage genome entry into host cells (e.g., PhageFISH); proximity of phage and host genomes (e.g., proximity ligation, polonies, viral tagging); and analysis of genomes and metagenomes (e.g., CRISPR spacer analysis, metagenomic co-occurrence). These methods have advantages and disadvantages. They also are not measuring the same interactions. Host range can be divided into multiple host ranges, each defined by how far the phage can progress in the infection cycle. For example, the ability to effect genome entry (penetrative host range) is different than the ability to produce progeny (productive host range). These different host ranges reflect bacterial defense mechanisms that block phage growth and development at various stages in the infection cycle. Here, I present a comparison of the various methods used to identify bacteriophage-host relationships with a focus on what type of host range is being measured or predicted.

Identification of a novel papillomavirus in oral swabs from giant pandas (Ailuropoda melanoleuca)

To fully characterize papillomavirus diversity in giant pandas (Ailuropoda melanoleuca), we identified a novel papillomavirus (named AmPV5, GenBank accession number MZ357114) in oral swabs from giant pandas with the help of viral metagenomics technology in this study. The complete circular genome of AmPV5 is 7,935 bp in length, with a GC content of 39.1%. It encodes five early genes (E1, E2, E4, E6, and E7), two late genes (L1 and L2), and features conserved zinc-binding domains (CXXC- (X)28/29-CXXC) in E6 and E7 genes. E7 protein has an LxCxE domain (pRB binding) in its N-terminal region. The nucleotide sequence of AmPV5 L1 gene shares < 70% identity with other related sequences available in the GenBank database. Phylogenetic analysis indicated that AmPV5 fell within the Lambdapapillomavirus genus but formed a monophyletic branch away from other papillomaviruses found in Ailuropoda melanoleuca, Canis, Felis catus, Panthera uncia, Enhydra lutris, and Procyon lotor. According to the International Committee on Taxonomy of Viruses (ICTV) classification guidelines, AmPV5 is classified as a new species within the Lambdapapillomavirus genus. The discovery provides valuable insights into the viral diversity in giant pandas and highlights the need for continued surveillance of wildlife pathogens. Future studies should explore the potential role of AmPV5 in the health and disease ecology of this endangered species.

Combined metabolomic and genomic analyses reveal phage-specific and infection stage-specific alterations to marine Roseobacter metabolism

Phages can reshape the metabolic network of hosts to support specific requirements for replication during infection. However, metabolomic profiling of phage-elicited host global metabolic alterations and the linkage of phage-encoded auxiliary metabolic genes to these alterations are understudied. In this study, the dynamics of intracellular metabolites of Dinoroseobacter shibae DFL12, a member of marine environmentally and biogeochemically relevant Roseobacter clade, in response to four distinct lytic roseophage infections were investigated. Metabolomic profiling indicated that roseophage infections significantly altered host metabolism in a phage-specific manner. Pathway enrichment analyses showed that the central carbon pathway and DNA, amino acid, and coenzyme metabolism were commonly altered by roseophages, revealing a central role of these pathways in phage replication. Furthermore, clear infection stage-specific host responses were observed, corresponding to different metabolic demands of phage replication in the early and late infection stages. Interestingly, the content of host vitamin B1, which is the essential nutrient provided by D. shibae to its symbiotic microalgae, increased in the early infection stage for most roseophages, implying that phage infection may impact the symbiosis of D. shibae with microalgae. Finally, combined metabolomic and phage genomics analyses showed that roseophages adopt different strategies to expand the host pyrimidine pool (recycling or de novo synthesis of pyrimidine nucleotides), and this difference was likely related to variation in the GC content between phage and host genomes. Collectively, these results highlight the potential importance of phage-specific and infection stage-specific host metabolic reprogramming in marine phage-host interactions, bacteria-microalgae symbiosis, and biogeochemical cycles.

Cave Pools in Carlsbad Caverns National Park Contain Diverse Bacteriophage Communities and Novel Viral Sequences

Viruses are the most abundant biological entities on Earth, and they play a critical role in the environment and biosphere where they regulate microbial populations and contribute to nutrient cycling. Environmental viruses have been the most studied in the ocean, but viral investigations have now spread to other environments. Here, viral communities were characterized in four cave pools in Carlsbad Caverns National Park to test the hypotheses that (i) viral abundance is ten-fold higher than prokaryotic cell abundance in cavern pools, (ii) cavern pools contain novel viral sequences, and (iii) viral communities in pools from developed portions of the cave are distinct from those of pools in undeveloped parts of the same cave. The relationship between viral and microbial abundance was determined through direct epifluorescence microscopy counts. Viral metagenomes were constructed to examine viral diversity among pools, identify novel viruses, and characterize auxiliary metabolic genes (AMGs). Bacterial communities were characterized by 16S rRNA gene amplicon sequencing. Epifluorescence microscopy showed that the ratio of viral-like particles (VLPs) to microorganisms was approximately 22:1 across all sites. Viral communities from pools with higher tourist traffic were more similar to each other than to those from less visited pools, although surprisingly, viruses did not follow the same pattern as bacterial communities, which reflected pool geochemistry. Bacterial hosts predicted from viral sequences using iPHoP showed overlap with both rare and abundant genera and families in the 16S rRNA gene dataset. Gene-sharing network analysis revealed high viral diversity compared to a reference viral database as well as to other aquatic environments. AMG presence showed variation in metabolic potential among the four pools. Overall, Carlsbad Cavern harbors novel viruses with substantial diversity among pools within the same system, indicating that caves are likely an important repository for unexplored viromes.

Method evaluation for viruses in activated sludge: Concentration, sequencing, and identification

Activated sludge (AS) in wastewater treatment plants is one of the largest artificial microbial ecosystems on earth and it makes enormous contributions to human societies. Viruses are an important component in AS with a high abundance. However, their communities and functionalities have not been as widely explored as those of other microorganisms, such as bacteria. This gap is mainly due to technical challenges in effective viral concentration, extraction, and sequencing. In this study, we compared four kinds of concentration methods, two sequencing approaches, and four identification bioinformatic tools to evaluate the whole analysis workflow for viruses in AS. Results showed flocculation, filtration, and resuspension (FFR) could get the longest DNA lengths and ultracentrifugation obtained the highest DNA yields for viruses in AS. Based on the results of present study, FFR and tangential flow filtration with the membrane pore size of 100 kDa were most recommended to concentrate viruses in AS samples with huge volumes. Besides, different concentration methods could get different viral catalogs and thus multiple methods should be combined to get the whole picture of viruses in the system. In addition, geNomad was the most recommended identification tool for viruses in the present study and the long-read sequencing could improve the assembly statistics of viruses when compared with the short-read sequencing. For the 8192 viral operational taxonomic units in this study, 95.1 % of them were phages and belonged to the same lineage at the order level of Caudovirales. Virulent phages dominated the AS system and Pseudomonadota were the main host. Taken together, this study provides new insights into methods selection for virus research of AS.

Metagenomic characterization of viruses and mobile genetic elements associated with the DPANN archaeal superphylum

The archaeal superphylum DPANN (an acronym formed from the initials of the first five phyla discovered: Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanohaloarchaeota and Nanoarchaeota) is a group of ultrasmall symbionts able to survive in extreme ecosystems. The diversity and dynamics between DPANN archaea and their virome remain largely unknown. Here we use a metagenomic clustered regularly interspaced short palindromic repeats (CRISPR) screening approach to identify 97 globally distributed, non-redundant viruses and unclassified mobile genetic elements predicted to infect hosts across 8 DPANN phyla, including 7 viral groups not previously characterized. Genomic analysis suggests a diversity of viral morphologies including head-tailed, tailless icosahedral and spindle-shaped viruses with the potential to establish lytic, chronic or lysogenic infections. We also find evidence of a virally encoded Cas12f1 protein (probably originating from uncultured DPANN archaea) and a mini-CRISPR array, which could play a role in modulating host metabolism. Many metagenomes have virus-to-host ratios >10, indicating that DPANN viruses play an important role in controlling host populations. Overall, our study illuminates the underexplored diversity, functional repertoires and host interactions of the DPANN virome.

The GC% landscape of the Nucleocytoviricota

Genomic studies on sequence composition employ various approaches, such as calculating the proportion of guanine and cytosine within a given sequence (GC% content), which can shed light on various aspects of the organism's biology. In this context, GC% can provide insights into virus-host relationships and evolution. Here, we present a comprehensive gene-by-gene analysis of 61 representatives belonging to the phylum Nucleocytoviricota, which comprises viruses with the largest genomes known in the virosphere. Parameters were evaluated not only based on the average GC% of a given viral species compared to the entire phylum but also considering gene position and phylogenetic history. Our results reveal that while some families exhibit similar GC% among their representatives (e.g., Marseilleviridae), others such as Poxviridae, Phycodnaviridae, and Mimiviridae have members with discrepant GC% values, likely reflecting adaptation to specific biological cycles and hosts. Interestingly, certain genes located at terminal regions or within specific genomic clusters show GC% values distinct from the average, suggesting recent acquisition or unique evolutionary pressures. Horizontal gene transfer and the presence of potential paralogs were also assessed in genes with the most discrepant GC% values, indicating multiple evolutionary histories. Taken together, to the best of our knowledge, this study represents the first global and gene-by-gene analysis of GC% distribution and profiles within genomes of Nucleocytoviricota members, highlighting their diversity and identifying potential new targets for future studies.

Comparative metagenomics of microbial communities and resistome in southern farming systems: implications for antimicrobial stewardship and public health

Agricultural practices significantly influence microbial diversity and the distribution of virulence and antimicrobial resistance (AMR) genes, with implications for ecosystem health and food safety. This study used metagenomic sequencing to analyze 60 samples (30 per state) including water, soil, and manure (10 each) from Alabama (a mix of cattle and poultry sources) and Tennessee (primarily from cattle). The results highlighted a rich microbial diversity, predominantly comprising Bacteria (67%) and Viruses (33%), with a total of over 1,950 microbial species identified. The dominant bacterial phyla were Proteobacteria, Cyanobacteria, Actinobacteria, Firmicutes, and Bacteroidetes, with the viral communities primarily represented by Phixviricota and Uroviricota. Distinct state-specific microbial profiles were evident, with Alabama demonstrating a higher prevalence of viral populations and unique bacterial phyla compared to Tennessee. The influence of environmental and agricultural practices was reflected in the microbial compositions: soil samples were notably rich in Actinobacteria, water samples were dominated by Proteobacteria and Cyanobacteria, and manure samples from Alabama showed a predominance of Actinobacteria. Further analyses, including diversity assessment and enterotype clustering, revealed complex microbial structures. Tennessee showed higher microbial diversity and phylogenetic complexity across most sample types compared to Alabama, with poultry-related samples displaying distinct diversity trends. Principal Coordinate Analysis (PCoA) highlighted notable state-specific variations, particularly in manure samples. Differential abundance analysis demonstrated elevated levels of Deinococcus and Ligilactobacillus in Alabama, indicating regional effects on microbial distributions. The virulome analysis revealed a significant presence of virulence genes in samples from Alabama. The community resistome was extensive, encompassing 109 AMR genes across 18 antibiotic classes, with manure samples displaying considerable diversity. Ecological analysis of the interactions between AMR gene subtypes and microbial taxa revealed a sophisticated network, often facilitated by bacteriophages. These findings underscore the critical role of agricultural practices in shaping microbial diversity and resistance patterns, highlighting the need for targeted AMR mitigation strategies in agricultural ecosystems to protect both public health and environmental integrity.

Discovery and biological confirmation of a highly divergent Tacaribe virus in metatranscriptomic data from neotropical bats

First isolated from neotropical fruit bats in Trinidad in 1956, Tacaribe virus (TCRV) has rarely been detected since. We searched for New World arenavirus reads in roughly 5.7 million sequencing runs available on public databases using Serratus. We recovered a complete genome of a divergent TCRV in metatranscriptomic data derived from heart and eye tissue of an adult male Jamaican fruit-eating bat sampled in the Dominican Republic, 2014. In total, 2,733 reads were mapped resulting in mean coverages of 7.4-fold for the L and 10.2-fold for the S segment. Re-testing original bat specimens showed the highest viral loads in liver tissue (245 copies/mg). Sanger sequencing of PCR amplicons from liver confirmed correctness of and completed the genome recovered from metatranscriptomic data, revealing conserved arenavirus genomic organization, length, intergenic regions, and genome termini. The newly found TCRV strain tentatively named DOM2014 clustered in a basal sister relationship to all other known TCRV strains with which it shared between 83.3%-86.0% genomic and 91.8%-93.7% translated amino acid sequence identity across protein-coding regions. DOM2014 showed a conserved glycine, proline, proline, threonine (GPPT) nucleoprotein motif, which is essential for TCRV interferon β antagonism. Our data confirm the association of TCRV with the bat genus Artibeus put into question by lethal experimental infections and scarce bat-derived TCRV genomic data. Broad genetic diversity and geographic spread require assessments of TCRV strain-associated pathogenicity, particularly for DOM2014 as a highly divergent TCRV strain. Confirmation of genomic database findings by testing original specimens provides robustness to our findings and supports the usefulness of metatranscriptomic studies.

IMPORTANCE

Clade B New World arenaviruses (NWA) include rodent-borne lethal hemorrhagic fever viruses, whereas Tacaribe virus (TCRV) stands out because of its detection in bats and its presumably low zoonotic potential. However, the bat association of TCRV was put into question by lethal experimental neotropical fruit bat infections and rare TCRV detection in bats. Scarce genomic data include near-identical viruses from Caribbean bats and ticks from the US sampled 50 years later. The prototype TCRV isolate used for experimental risk assessments has an extensive passage history in suckling mouse brains. Exploring the true genetic diversity, geographic distribution, and host range of bat-borne NWA is pivotal to assess their zoonotic potential and transmission cycles. We analyzed metatranscriptomic data for evidence of NWA identifying a highly divergent TCRV in bats and confirmed virus detection in original biological materials, supporting the association of TCRV with neotropical bats and warranting investigation of strain-associated TCRV pathogenicity.

Unveiling the Virome of Wild Birds: Exploring CRESS-DNA Viral Dark Matter

Amid global health concerns and the constant threat of zoonotic diseases, this study delves into the diversity of circular replicase-encoding single-stranded DNA (CRESS-DNA) viruses within Chinese wild bird populations. Employing viral metagenomics to tackle the challenge of "viral dark matter," the research collected and analyzed 3,404 cloacal swab specimens across 26 bird families. Metagenomic analysis uncovered a rich viral landscape, with 67.48% of reads classified as viral dark matter, spanning multiple taxonomic levels. Notably, certain viral families exhibited host-specific abundance patterns, with Galliformes displaying the highest diversity. Diversity analysis categorized samples into distinct groups, revealing significant differences in viral community structure, particularly noting higher diversity in terrestrial birds compared to songbirds and unique diversity in migratory birds versus perching birds. The identification of ten novel Circoviridae viruses, seven Smacoviridae viruses, and 167 Genomoviridae viruses, along with 100 unclassified CRESS-DNA viruses, underscores the expansion of knowledge on avian-associated circular DNA viruses. Phylogenetic and structural analyses of Rep proteins offered insights into evolutionary relationships and potential functional variations among CRESS-DNA viruses. In conclusion, this study significantly enhances our understanding of the avian virome, shedding light on the intricate relationships between viral communities and host characteristics in Chinese wild bird populations. The diverse array of CRESS-DNA viruses discovered opens avenues for future research into viral evolution, spread factors, and potential ecosystem impacts.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

BACKGROUND

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

RESULTS

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes had greater species richness and total viral genome abundances than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

CONCLUSIONS

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits. Video Abstract.

Biodiversity of mudflat intertidal viromes along the Chinese coasts

Viruses constitute the most diverse and abundant biological entities on Earth. However, our understanding of this tiniest life form in complex ecosystems remains limited. Here, we recover 20,102 viral OTUs from twelve intertidal zones along the Chinese coasts. Our analysis demonstrates high viral diversity and functional potential in intertidal zones, encoding important functional genes that can be potentially transferred to microbial hosts and mediate elemental biogeochemical cycles, especially carbon, phosphate and sulfur. Virus-host abundance dynamics vary among different microbial lineages. Viral community composition is closely associated with environmental conditions, including dissolved organic matter. Concordant biogeographic patterns are observed for viruses and microbes. Viral communities are generally habitat specific with low overlaps between intertidal and other habitats. Environmental factors and geographic distance dominate the compositional variation of intertidal viromes. Overall, these findings expand our understanding of intertidal viromes within an ecological framework, providing insights into the virus-host coevolutionary arms race.

Characterization and comparative genomic analysis of a marine Bacillus phage reveal a novel viral genus

Bacillus pumilus exhibits substantial economic significance, with its metabolism, adaptability, and ecological functions regulated by its bacteriophages. Here, we isolated and characterized a novel temperate phage vB_BpuM-ZY1 from B. pumilus derived from mangrove sediments by mitomycin C induction. Phage vB_BpuM-ZY1 is a typical myophage, which has an icosahedral head with a diameter of 43.34 ± 2.14 nm and a long contractible tail with a length of 238.58 ± 5.18 nm. Genomic analysis indicated that vB_BpuM-ZY1 encodes genes for lysogeny control, and its life cycle may be intricately regulated by multiple mechanisms. vB_BpuM-ZY1 was predicted to employ P2-like 5'-extended-cos packaging strategy. In addition, genome-wide phylogenetic tree and proteome tree analyses indicated that vB_BpuM-ZY1 belongs to the Peduoviridae family but forms a separate branch at a deeper taxonomic level. Particularly, the comparative genomic analysis showed that vB_BpuM-ZY1 has less than 70% intergenomic similarities with its most similar phages. Thus, we propose that vB_BpuM-ZY1 is a novel Bacillus phage belonging to a new genus under the Peduoviridae family. The protein-sharing network analysis identified 44 vB_BpuM-ZY1-related phages. Interestingly, these evolutionarily related myophages infect a broad range of hosts across different phyla, which may be explained by the high structural variations of the host recognition domain in their central spike proteins. Collectively, our study will contribute to our understanding of Bacillus phage diversity and Bacillus-phage interactions, as well as provide essential knowledge for the industrial application of B. pumilus.

IMPORTANCE

Although recent metagenomics research has obtained a wealth of phage genetic information, much of it is considered "dark matter" because of the lack of similarity with known sequences in the database. Therefore, the isolation and characterization of novel phages will help to interpret the vast unknown viral metagenome data and improve our understanding of phage diversity and phage-host interactions. Bacillus pumilus shows high economic relevance due to its wide applications in biotechnology, industry, biopharma, and environmental sectors. Since phages influence the abundance, metabolism, evolution, fitness, and ecological functions of bacteria through complex interactions, the significance of isolation and characterization of novel phages infecting B. pumilus is apparent. In this study, we isolated and characterized a B. pumilus phage belonging to a novel viral genus, which provides essential knowledge for phage biology as well as the industrial application of B. pumilus.

Another piece of puzzle for the human microbiome: the gut virome under dietary modulation

The virome is the most abundant and highly variable microbial consortium in the gut. Because of difficulties in isolating and culturing gut viruses and the lack of reference genomes, the virome has remained a relatively elusive aspect of the human microbiome. In recent years, studies on the virome have accumulated growing evidence showing that the virome is diet-modulated and widely involved in regulating health. Here, we review the responses of the gut virome to dietary intake and the potential health implications, presenting changes in the gut viral community and preferences of viral members to particular diets. We further discuss how viral-bacterial interactions and phage lifestyle shifts shape the gut microbiota. We also discuss the specific functions conferred by diet on the gut virome and bacterial community in the context of horizontal gene transfer, as well as the import of new viral members along with the diet. Collating these studies will expand our understanding of the dietary regulation of the gut virome and inspire dietary interventions and health maintenance strategies targeting the gut microbiota.

Metagenomic analysis reveals high diversity of gut viromes in yaks (Bos grunniens) from the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP), renowned for its exceptional biological diversity, is home to numerous endemic species. However, research on the virology of vulnerable vertebrates like yaks remains limited. In this study, our objective was to use metagenomics to provide a comprehensive understanding of the diversity and evolution of the gut virome in yak populations across different regions of the QTP. Our findings revealed a remarkably diverse array of viruses in the gut of yaks, including those associated with vertebrates and bacteriophages. Notably, some vertebrate-associated viruses, such as astrovirus and picornavirus, showed significant sequence identity across diverse yak populations. Additionally, we observed differences in the functional profiles of genes carried by the yak gut virome across different regions. Moreover, the virus-bacterium symbiotic network that we discovered holds potential significance in maintaining the health of yaks. Overall, this research expands our understanding of the viral communities in the gut of yaks and highlights the importance of further investigating the interactions between viruses and their hosts. These data will be beneficial for revealing the crucial role that viruses play in the yak gut ecology in future studies.

Control of lysogeny and antiphage defense by a prophage-encoded kinase-phosphatase module

The filamentous 'Pf' bacteriophages of Pseudomonas aeruginosa play roles in biofilm formation and virulence, but mechanisms governing Pf prophage activation in biofilms are unclear. Here, we identify a prophage regulatory module, KKP (kinase-kinase-phosphatase), that controls virion production of co-resident Pf prophages and mediates host defense against diverse lytic phages. KKP consists of Ser/Thr kinases PfkA and PfkB, and phosphatase PfpC. The kinases have multiple host targets, one of which is MvaU, a host nucleoid-binding protein and known prophage-silencing factor. Characterization of KKP deletion and overexpression strains with transcriptional, protein-level and prophage-based approaches indicates that shifts in the balance between kinase and phosphatase activities regulate phage production by controlling MvaU phosphorylation. In addition, KKP acts as a tripartite toxin-antitoxin system that provides defense against some lytic phages. A conserved lytic phage replication protein inhibits the KKP phosphatase PfpC, stimulating toxic kinase activity and blocking lytic phage production. Thus, KKP represents a phosphorylation-based mechanism for prophage regulation and antiphage defense. The conservation of KKP gene clusters in >1000 diverse temperate prophages suggests that integrated control of temperate and lytic phage infection by KKP-like regulatory modules may play a widespread role in shaping host cell physiology.

Long-read sequencing reveals extensive gut phageome structural variations driven by genetic exchange with bacterial hosts

Genetic variations are instrumental for unraveling phage evolution and deciphering their functional implications. Here, we explore the underlying fine-scale genetic variations in the gut phageome, especially structural variations (SVs). By using virome-enriched long-read metagenomic sequencing across 91 individuals, we identified a total of 14,438 nonredundant phage SVs and revealed their prevalence within the human gut phageome. These SVs are mainly enriched in genes involved in recombination, DNA methylation, and antibiotic resistance. Notably, a substantial fraction of phage SV sequences share close homology with bacterial fragments, with most SVs enriched for horizontal gene transfer (HGT) mechanism. Further investigations showed that these SV sequences were genetic exchanged between specific phage-bacteria pairs, particularly between phages and their respective bacterial hosts. Temperate phages exhibit a higher frequency of genetic exchange with bacterial chromosomes and then virulent phages. Collectively, our findings provide insights into the genetic landscape of the human gut phageome.

A natural ANI gap that can define intra-species units of bacteriophages and other viruses

Despite the importance of intra-species variants of viruses for causing disease and/or disrupting ecosystem functioning, there is no universally applicable standard to define these. A (natural) gap in whole-genome average nucleotide identity (ANI) values around 95% is commonly used to define species, especially for bacteriophages, but whether a similar gap exists within species that can be used to define intra-species units has not been evaluated yet. Whole-genome comparisons among members of 1,016 bacteriophage (Caudoviricetes) species revealed a region of low frequency of ANI values around 99.2%-99.8%, showing threefold or fewer pairs than expected for an even distribution. This second gap is prevalent in viruses infecting various cultured or uncultured hosts from a variety of environments, although a few exceptions to this pattern were also observed (3.7% of total species) and are likely attributed to cultivation biases or other factors. Similar results were observed for a limited set of eukaryotic viruses that are adequately sampled, including SARS-CoV-2, whose ANI-based clusters matched well with the WHO-defined variants of concern, indicating that our findings from bacteriophages might be more broadly applicable and the ANI-based clusters may represent functionally and/or ecologically distinct units. These units appear to be predominantly driven by (high) ecological cohesiveness coupled to either frequent recombination for bacteriophages or selection and clonal evolution for other viruses such as SARS-CoV-2, indicating that fundamentally different underlying mechanisms could lead to similar diversity patterns. Accordingly, we propose the ANI gap approach outlined above for defining viral intra-species units, for which we propose the term genomovars.

IMPORTANCE

Viral species are composed of an ensemble of intra-species variants whose individual dynamics may have major implications for human and animal health and/or ecosystem functioning. However, the lack of universally accepted standards to define these intra-species variants has led researchers to use different approaches for this task, creating inconsistent intra-species units across different viral families and confusion in communication. By comparing hundreds of mostly bacteriophage genomes, we show that there is a widely distributed natural gap in whole-genome average nucleotide identity values in most, but not all, of these species that can be used to define intra-species units. Therefore, these results advance the molecular toolbox for tracking viral intra-species units and should facilitate future epidemiological and environmental studies.

An enterococcal phage-derived enzyme suppresses graft-versus-host disease

Changes in the gut microbiome have pivotal roles in the pathogenesis of acute graft-versus-host disease (aGVHD) after allogenic haematopoietic cell transplantation (allo-HCT)1-6. However, effective methods for safely resolving gut dysbiosis have not yet been established. An expansion of the pathogen Enterococcus faecalis in the intestine, associated with dysbiosis, has been shown to be a risk factor for aGVHD7-10. Here we analyse the intestinal microbiome of patients with allo-HCT, and find that E. faecalis escapes elimination and proliferates in the intestine by forming biofilms, rather than by acquiring drug-resistance genes. We isolated cytolysin-positive highly pathogenic E. faecalis from faecal samples and identified an anti-E. faecalis enzyme derived from E. faecalis-specific bacteriophages by analysing bacterial whole-genome sequencing data. The antibacterial enzyme had lytic activity against the biofilm of E. faecalis in vitro and in vivo. Furthermore, in aGVHD-induced gnotobiotic mice that were colonized with E. faecalis or with patient faecal samples characterized by the domination of Enterococcus, levels of intestinal cytolysin-positive E. faecalis were decreased and survival was significantly increased in the group that was treated with the E. faecalis-specific enzyme, compared with controls. Thus, administration of a phage-derived antibacterial enzyme that is specific to biofilm-forming pathogenic E. faecalis-which is difficult to eliminate with existing antibiotics-might provide an approach to protect against aGVHD.

Asgard archaea defense systems and their roles in the origin of eukaryotic immunity

Dozens of new antiviral systems have been recently characterized in bacteria. Some of these systems are present in eukaryotes and appear to have originated in prokaryotes, but little is known about these defense mechanisms in archaea. Here, we explore the diversity and distribution of defense systems in archaea and identify 2610 complete systems in Asgardarchaeota, a group of archaea related to eukaryotes. The Asgard defense systems comprise 89 unique systems, including argonaute, NLR, Mokosh, viperin, Lassamu, and CBASS. Asgard viperin and argonaute proteins have structural homology to eukaryotic proteins, and phylogenetic analyses suggest that eukaryotic viperin proteins were derived from Asgard viperins. We show that Asgard viperins display anti-phage activity when heterologously expressed in bacteria. Eukaryotic and bacterial argonaute proteins appear to have originated in Asgardarchaeota, and Asgard argonaute proteins have argonaute-PIWI domains, key components of eukaryotic RNA interference systems. Our results support that Asgardarchaeota played important roles in the origin of antiviral defense systems in eukaryotes.

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.

Genome Characterization of Two Novel Lactococcus lactis Phages vL_296 and vL_20A

Fermented dairy products are produced using starter cultures. They ferment milk to create products with a certain texture, aroma, and taste. However, the lactic acid bacteria used in this production are prone to bacteriophage infection. We examined the genomes of two newly discovered bacteriophage species that were isolated from cheese whey during the cheesemaking process. We have determined the species and the lytic spectrum of these bacteriophages. Phages vL_20A and vL_296 were isolated using lactococcal indicator cultures. They have unique lytic spectra: of the 21 possible identified host bacteria, only four are shared amongst them. The vL_20A and vL_296 genomes comprise linear double-stranded DNA lengths with 21,909 and 22,667 nucleotide pairs, respectively. Lactococcus phage bIL67 (ANI 93.3 and 92.6, respectively) is the closest to the phages vL_20A and vL_296. The analysis of the CRISPR spacers in the genomes of starter cultures did not reveal any phage-specific vL_20A or vL_296 among them. This study highlights the biodiversity of L. lactis phages, their widespread presence in dairy products, and their virulence. However, the virulence of phages is balanced by the presence of a significant number of bacterial strains with different sensitivities to phages in the starter cultures due to the bacterial immune system.

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.

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.

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.

Deep mining of the Sequence Read Archive reveals major genetic innovations in coronaviruses and other nidoviruses of aquatic vertebrates

Virus discovery by genomics and metagenomics empowered studies of viromes, facilitated characterization of pathogen epidemiology, and redefined our understanding of the natural genetic diversity of viruses with profound functional and structural implications. Here we employed a data-driven virus discovery approach that directly queries unprocessed sequencing data in a highly parallelized way and involves a targeted viral genome assembly strategy in a wide range of sequence similarity. By screening more than 269,000 datasets of numerous authors from the Sequence Read Archive and using two metrics that quantitatively assess assembly quality, we discovered 40 nidoviruses from six virus families whose members infect vertebrate hosts. They form 13 and 32 putative viral subfamilies and genera, respectively, and include 11 coronaviruses with bisegmented genomes from fishes and amphibians, a giant 36.1 kilobase coronavirus genome with a duplicated spike glycoprotein (S) gene, 11 tobaniviruses and 17 additional corona-, arteri-, cremega-, nanhypo- and nangoshaviruses. Genome segmentation emerged in a single evolutionary event in the monophyletic lineage encompassing the subfamily Pitovirinae. We recovered the bisegmented genome sequences of two coronaviruses from RNA samples of 69 infected fishes and validated the presence of poly(A) tails at both segments using 3'RACE PCR and subsequent Sanger sequencing. We report a genetic linkage between accessory and structural proteins whose phylogenetic relationships and evolutionary distances are incongruent with the phylogeny of replicase proteins. We rationalize these observations in a model of inter-family S recombination involving at least five ancestral corona- and tobaniviruses of aquatic hosts. In support of this model, we describe an individual fish co-infected with members from the families Coronaviridae and Tobaniviridae. Our results expand the scale of the known extraordinary evolutionary plasticity in nidoviral genome architecture and call for revisiting fundamentals of genome expression, virus particle biology, host range and ecology of vertebrate nidoviruses.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Forty-nine metagenomic-assembled genomes from an aquatic virome expand Caudoviricetes by 45 potential new families and the newly uncovered Gossevirus of Bamfordvirae

Twenty complete genomes (29-63 kb) and 29 genomes with an estimated completeness of over 90 % (30-90 kb) were identified for novel dsDNA viruses in the Yangshan Harbor metavirome. These newly discovered viruses contribute to the expansion of viral taxonomy by introducing 46 potential new families. Except for one virus, all others belong to the class Caudoviricetes. The exception is a novel member of the recently characterized viral group known as Gossevirus. Fifteen viruses were predicted to be temperate. The predicted hosts for the viruses appear to be involved in various aspects of the nitrogen cycle, including nitrogen fixation, oxidation and denitrification. Two viruses were identified to have a host of Flavobacterium and Tepidimonas fonticaldi, respectively, by matching CRISPR spacers with viral protospacers. Our findings provide an overview for characterizing and identifying specific viruses from Yangshan Harbor. The Gossevirus-like virus uncovered emphasizes the need for further comprehensive isolation and investigation of polinton-like viruses.

Tips and tools to obtain and assess mosquito viromes

Due to their vectorial capacity, mosquitoes (Diptera: Culicidae) receive special attention from health authorities and entomologists. These cosmopolitan insects are responsible for the transmission of many viral diseases, such as dengue and yellow fever, causing huge impacts on human health and justifying the intensification of research focused on mosquito-borne diseases. In this context, the study of the virome of mosquitoes can contribute to anticipate the emergence and/or the reemergence of infectious diseases. The assessment of mosquito viromes also contributes to the surveillance of a wide variety of viruses found in these insects, allowing the early detection of pathogens with public health importance. However, the study of mosquito viromes can be challenging due to the number and complexities of steps involved in this type of research. Therefore, this article aims to describe, in a straightforward and simplified way, the steps necessary for obtention and assessment of mosquito viromes. In brief, this article explores: the capture and preservation of specimens; sampling strategies; treatment of samples before DNA/RNA extraction; extraction methodologies; enrichment and purification processes; sequencing choices; and bioinformatics analysis.

Large-scale single-virus genomics uncovers hidden diversity of river water viruses and diversified gene profiles

Environmental viruses (primarily bacteriophages) are widely recognized as playing an important role in ecosystem homeostasis through the infection of host cells. However, the majority of environmental viruses are still unknown as their mosaic structure and frequent mutations in their sequences hinder genome construction in current metagenomics. To enable the large-scale acquisition of environmental viral genomes, we developed a new single-viral genome sequencing platform with microfluidic-generated gel beads. Amplification of individual DNA viral genomes in mass-produced gel beads allows high-throughput genome sequencing compared to conventional single-virus genomics. The sequencing analysis of river water samples yielded 1431 diverse viral single-amplified genomes, whereas viral metagenomics recovered 100 viral metagenome-assembled genomes at the comparable sequence depth. The 99.5% of viral single-amplified genomes were determined novel at the species level, most of which could not be recovered by a metagenomic assembly. The large-scale acquisition of diverse viral genomes identified protein clusters commonly detected in different viral strains, allowing the gene transfer to be tracked. Moreover, comparative genomics within the same viral species revealed that the profiles of various methyltransferase subtypes were diverse, suggesting an enhanced escape from host bacterial internal defense mechanisms. Our use of gel bead-based single-virus genomics will contribute to exploring the nature of viruses by accelerating the accumulation of draft genomes of environmental DNA viruses.

Temperature modulates dominance of a superinfecting Arctic virus in its unicellular algal host

Complex virus-virus interactions can arise when multiple viruses coinfect the same host, impacting infection outcomes with broader ecological and evolutionary significance for viruses and host. Yet, our knowledge regarding virus competition is still limited, especially for single-celled eukaryotic host-virus systems. Here, we report on mutual interference of two dsDNA viruses, MpoV-45T and MpoV-46T, competing for their Arctic algal host Micromonas polaris. Both viruses affected each other's gene expression and displayed reduced genome replication during coinfection. MpoV-45T was the dominant virus, likely due to interference in the DNA replication of is competitor. Even when its coinfection was delayed, the dominant virus still prevailed while genome production of the other virus was strongly suppressed. This contrasts with typical superinfection exclusion, where the primary infection prevents secondary infection by other viruses. Higher temperature made the suppressed virus a stronger competitor, signifying that global warming is likely to alter virus-virus interactions in Arctic waters.

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.

The Diversity of Bacteriophages in the Human Gut

Bacteriophages, commonly referred to as phages, are viruses that infect bacteria and are among the most numerous microorganisms on the planet. They occur throughout nature occupying every habitat where their bacterial hosts can be found. Within these communities, phages are responsible for shaping the bacterial community structure and function through their interactions. Phages shape the community structure and function within the human gut but are also able to influence the human host. As such, there is increased interest in understanding the composition and activity of the gastrointestinal phages, although these studies have been hindered by the difficulties accompanying the study of the human gut. Here, we summarize the methods and findings pertaining to the diversity of the human gastrointestinal phages.

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.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

Background

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

Results

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes were more abundant and species rich than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

Conclusions

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits.

The virome in allergy and asthma: A nascent, ineffable player

Allergic diseases can be affected by virus-host interactions and are increasingly linked with the tissue-specific microbiome. High-throughput metagenomic sequencing has offered the opportunity to study the presence of viruses as an ecologic system, namely, the virome. Even though virome studies are technically challenging conceptually and analytically, they are already producing novel data expanding our understanding of the pathophysiologic mechanisms related to chronic inflammation and allergy. The importance of interspecies and intraspecies interactions is becoming apparent, as they can significantly, directly or indirectly, affect the host's response and antigenic state. Here, we emphasize the challenges and potential insights related to study of the virome in the context of allergy and asthma. We review the limited number of studies that have investigated the virome in these conditions, underlining the need for prospective, repeated sampling designs to unravel the virome's impact on disease development and its interplay with microbiota and immunity. The potential therapeutic use of bacteriophages, which are highly complex components of the virome, is discussed. There is clearly a need for further in-depth investigation of the virome as a system in allergic diseases.

VIGA: a one-stop tool for eukaryotic virus identification and genome assembly from next-generation-sequencing data

Identification of viruses and further assembly of viral genomes from the next-generation-sequencing data are essential steps in virome studies. This study presented a one-stop tool named VIGA (available at https://github.com/viralInformatics/VIGA) for eukaryotic virus identification and genome assembly from NGS data. It was composed of four modules, namely, identification, taxonomic annotation, assembly and novel virus discovery, which integrated several third-party tools such as BLAST, Trinity, MetaCompass and RagTag. Evaluation on multiple simulated and real virome datasets showed that VIGA assembled more complete virus genomes than its competitors on both the metatranscriptomic and metagenomic data and performed well in assembling virus genomes at the strain level. Finally, VIGA was used to investigate the virome in metatranscriptomic data from the Human Microbiome Project and revealed different composition and positive rate of viromes in diseases of prediabetes, Crohn's disease and ulcerative colitis. Overall, VIGA would help much in identification and characterization of viromes, especially the known viruses, in future studies.

Characterization of phage vB_EcoS-EE09 infecting E. coli DSM613 Isolated from Wastewater Treatment Plant Effluent and Comparative Proteomics of the Infected and Non-Infected Host

Phages influence microbial communities, can be applied in phage therapy, or may serve as bioindicators, e.g., in (waste)water management. We here characterized the Escherichia phage vB_EcoS-EE09 isolated from an urban wastewater treatment plant effluent. Phage vB_EcoS-EE09 belongs to the genus Dhillonvirus, class Caudoviricetes. It has an icosahedral capsid with a long non-contractile tail and a dsDNA genome with an approximate size of 44 kb and a 54.6% GC content. Phage vB_EcoS-EE09 infected 12 out of the 17 E. coli strains tested. We identified 16 structural phage proteins, including the major capsid protein, in cell-free lysates by protein mass spectrometry. Comparative proteomics of protein extracts of infected E. coli cells revealed that proteins involved in amino acid and protein metabolism were more abundant in infected compared to non-infected cells. Among the proteins involved in the stress response, 74% were less abundant in the infected cultures compared to the non-infected controls, with six proteins showing significant less abundance. Repressing the expression of these proteins may be a phage strategy to evade host defense mechanisms. Our results contribute to diversifying phage collections, identifying structural proteins to enable better reliability in annotating taxonomically related phage genomes, and understanding phage-host interactions at the protein level.

Extracellular vesicles are the main contributor to the non-viral protected extracellular sequence space

Environmental virus metagenomes, commonly referred to as "viromes", are typically generated by physically separating virus-like particles (VLPs) from the microbial fraction based on their size and mass. However, most methods used to purify VLPs, enrich extracellular vesicles (EVs) and gene transfer agents (GTAs) simultaneously. Consequently, the sequence space traditionally referred to as a "virome" contains host-associated sequences, transported via EVs or GTAs. We therefore propose to call the genetic material isolated from size-fractionated (0.22 µm) and DNase-treated samples protected environmental DNA (peDNA). This sequence space contains viral genomes, DNA transduced by viruses and DNA transported in EVs and GTAs. Since there is no genetic signature for peDNA transported in EVs, GTAs and virus particles, we rely on the successful removal of contaminating remaining cellular and free DNA when analyzing peDNA. Using marine samples collected from the North Sea, we generated a thoroughly purified peDNA dataset and developed a bioinformatic pipeline to determine the potential origin of the purified DNA. This pipeline was applied to our dataset as well as existing global marine "viromes". Through this pipeline, we identified known GTA and EV producers, as well as organisms with actively transducing proviruses as the source of the peDNA, thus confirming the reliability of our approach. Additionally, we identified novel and widespread EV producers, and found quantitative evidence suggesting that EV-mediated gene transfer plays a significant role in driving horizontal gene transfer (HGT) in the world's oceans.

A systematic analysis of marine lysogens and proviruses

Viruses are ubiquitous in the oceans, exhibiting high abundance and diversity. Here, we systematically analyze existing genomic sequences of marine prokaryotes to compile a Marine Prokaryotic Genome Dataset (MPGD, consisting of over 12,000 bacterial and archaeal genomes) and a Marine Temperate Viral Genome Dataset (MTVGD). At least 40% of the MPGD genomes contain one or more proviral sequences, indicating that they are lysogens. The MTVGD includes over 12,900 viral contigs or putative proviruses, clustered into 10,897 viral genera. We show that lysogens and proviruses are abundant in marine ecosystems, particularly in the deep sea, and marine lysogens differ from non-lysogens in multiple genomic features and growth properties. We reveal several virus-host interaction networks of potential ecological relevance, and identify proviruses that appear to be able to infect (or to be transferred between) different bacterial classes and phyla. Auxiliary metabolic genes in the MTVGD are enriched in functions related to carbohydrate metabolism. Finally, we experimentally demonstrate the impact of a prophage on the transcriptome of a representative marine Shewanella bacterium. Our work contributes to a better understanding of the ecology of marine prokaryotes and their viruses.

Factors Affecting Variation of the Human Gut Phageome

The gut microbiome is a dense and metabolically active consortium of microorganisms and viruses located in the lower gastrointestinal tract of the human body. Bacteria and their viruses (phages) are the most abundant members of the gut microbiome. Investigating their biology and the interplay between the two is important if we are to understand their roles in human health and disease. In this review, we summarize recent advances in resolving the taxonomic structure and ecological functions of the complex community of phages in the human gut-the gut phageome. We discuss how age, diet, and geography can all have a significant impact on phageome composition. We note that alterations to the gut phageome have been observed in several diseases such as inflammatory bowel disease, irritable bowel syndrome, and colorectal cancer, and we evaluate whether these phageome changes can directly or indirectly contribute to disease etiology and pathogenesis. We also highlight how lack of standardization in studying the gut phageome has contributed to variation in reported results.

Characterization and Comparative Genomic Analysis of a Deep-Sea Bacillus Phage Reveal a Novel Genus

As the most abundant biological entities, viruses are the major players in marine ecosystems. However, our knowledge on virus diversity and virus-host interactions in the deep sea remains very limited. In this study, vB_BteM-A9Y, a novel bacteriophage infecting Bacillus tequilensis, was isolated from deep-sea sediments in the South China Sea. vB_BteM-A9Y has a hexametric head and a long, complex contractile tail, which are typical features of myophages. vB_BteM-A9Y initiated host lysis at 60 min post infection with a burst size of 75 PFU/cell. The phage genome comprises 38,634 base pairs and encodes 54 predicted open reading frames (ORFs), of which 27 ORFs can be functionally annotated by homology analysis. Interestingly, abundant ORFs involved in DNA damage repair were identified in the phage genome, suggesting that vB_BteM-A9Y encodes multiple pathways for DNA damage repair, which may help to maintain the stability of the host/phage genome. A BLASTn search of the whole genome sequence of vB_BteM-A9Y against the GenBank revealed no existing homolog. Consistently, a phylogenomic tree and proteome-based phylogenetic tree analysis showed that vB_BteM-A9Y formed a unique branch. Further comparative analysis of genomic nucleotide similarity and ORF homology of vB_BteM-A9Y with its mostly related phages showed that the intergenomic similarity between vB_BteM-A9Y and these phages was 0-33.2%. Collectively, based on the comprehensive morphological, phylogenetic, and comparative genomic analysis, we propose that vB_BteM-A9Y belongs to a novel genus under Caudoviricetes. Therefore, our study will increase our knowledge on deep-sea virus diversity and virus-host interactions, as well as expanding our knowledge on phage taxonomy.