Metagenomic compendium of 189,680 DNA viruses from the human gut microbiome

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

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

Viral activation and ecological restructuring characterize a microbiome axis of spaceflight-associated immune activation

Maintenance of astronaut health during spaceflight will require monitoring and potentially modulating their microbiomes, which play a role in some space-derived health disorders. However, documenting the response of microbiota to spaceflight has been difficult thus far due to mission constraints that lead to limited sampling. Here, we executed a six-month longitudinal study centered on a three-day flight to quantify the high-resolution microbiome response to spaceflight. Via paired metagenomics and metatranscriptomics alongside single immune profiling, we resolved a microbiome "architecture" of spaceflight characterized by time-dependent and taxonomically divergent microbiome alterations across 750 samples and ten body sites. We observed pan-phyletic viral activation and signs of persistent changes that, in the oral microbiome, yielded plaque-associated pathobionts with strong associations to immune cell gene expression. Further, we found enrichments of microbial genes associated with antibiotic production, toxin-antitoxin systems, and stress response enriched universally across the body sites. We also used strain-level tracking to measure the potential propagation of microbial species from the crew members to each other and the environment, identifying microbes that were prone to seed the capsule surface and move between the crew. Finally, we identified associations between microbiome and host immune cell shifts, proposing both a microbiome axis of immune changes during flight as well as the sources of some of those changes. In summary, these datasets and methods reveal connections between crew immunology, the microbiome, and their likely drivers and lay the groundwork for future microbiome studies of spaceflight.

Choice of Ultrafilter Affects Recovery Rate of Bacteriophages

Studies into the viral fraction of complex microbial communities, like in the mammalian gut, have recently garnered much interest. Yet there is still no standardized protocol for extracting viruses from such samples, and the protocols that exist employ procedures that skew the viral community of the sample one way or another. The first step of the extraction pipeline often consists of the basic filtering of macromolecules and bacteria, yet even this affects the viruses in a strain-specific manner. In this study, we investigate a protocol for viral extraction based on ultrafiltration and how the choice of ultrafilter might influence the extracted viral community. Clinical samples (feces, vaginal swabs, and tracheal suction samples) were spiked with a mock community of known phages (T4, c2, Φ6, Φ29, Φx174, and Φ2972), filtered, and quantified using spot and plaque assays to estimate the loss in recovery. The enveloped Φ6 phage is especially severely affected by the choice of filter, but also tailed phages such as T4 and c2 have a reduced infectivity after ultrafiltration. We conclude that the pore size of ultrafilters may affect the recovery of phages in a strain- and sample-dependent manner, suggesting the need for greater thought when selecting filters for virus extraction.

Phables: from fragmented assemblies to high-quality bacteriophage genomes

MOTIVATION

Microbial communities have a profound impact on both human health and various environments. Viruses infecting bacteria, known as bacteriophages or phages, play a key role in modulating bacterial communities within environments. High-quality phage genome sequences are essential for advancing our understanding of phage biology, enabling comparative genomics studies and developing phage-based diagnostic tools. Most available viral identification tools consider individual sequences to determine whether they are of viral origin. As a result of challenges in viral assembly, fragmentation of genomes can occur, and existing tools may recover incomplete genome fragments. Therefore, the identification and characterization of novel phage genomes remain a challenge, leading to the need of improved approaches for phage genome recovery.

RESULTS

We introduce Phables, a new computational method to resolve phage genomes from fragmented viral metagenome assemblies. Phables identifies phage-like components in the assembly graph, models each component as a flow network, and uses graph algorithms and flow decomposition techniques to identify genomic paths. Experimental results of viral metagenomic samples obtained from different environments show that Phables recovers on average over 49% more high-quality phage genomes compared to existing viral identification tools. Furthermore, Phables can resolve variant phage genomes with over 99% average nucleotide identity, a distinction that existing tools are unable to make.

AVAILABILITY AND IMPLEMENTATION

Phables is available on GitHub at https://github.com/Vini2/phables.

The multi-kingdom microbiome of the goat gastrointestinal tract

BACKGROUND

Goat is an important livestock worldwide, which plays an indispensable role in human life by providing meat, milk, fiber, and pelts. Despite recent significant advances in microbiome studies, a comprehensive survey on the goat microbiomes covering gastrointestinal tract (GIT) sites, developmental stages, feeding styles, and geographical factors is still unavailable. Here, we surveyed its multi-kingdom microbial communities using 497 samples from ten sites along the goat GIT.

RESULTS

We reconstructed a goat multi-kingdom microbiome catalog (GMMC) including 4004 bacterial, 71 archaeal, and 7204 viral genomes and annotated over 4,817,256 non-redundant protein-coding genes. We revealed patterns of feeding-driven microbial community dynamics along the goat GIT sites which were likely associated with gastrointestinal food digestion and absorption capabilities and disease risks, and identified an abundance of large intestine-enriched genera involved in plant fiber digestion. We quantified the effects of various factors affecting the distribution and abundance of methane-producing microbes including the GIT site, age, feeding style, and geography, and identified 68 virulent viruses targeting the methane producers via a comprehensive virus-bacterium/archaea interaction network.

CONCLUSIONS

Together, our GMMC catalog provides functional insights of the goat GIT microbiota through microbiome-host interactions and paves the way to microbial interventions for better goat and eco-environmental qualities. Video Abstract.

Viruses in deep-sea cold seep sediments harbor diverse survival mechanisms and remain genetically conserved within species

Deep sea cold seep sediments have been discovered to harbor novel, abundant, and diverse bacterial and archaeal viruses. However, little is known about viral genetic features and evolutionary patterns in these environments. Here, we examined the evolutionary ecology of viruses across active and extinct seep stages in the area of Haima cold seeps in the South China Sea. A total of 338 viral operational taxonomic units are identified and linked to 36 bacterial and archaeal phyla. The dynamics of host-virus interactions are informed by diverse antiviral defense systems across 43 families found in 487 microbial genomes. Cold seep viruses are predicted to harbor diverse adaptive strategies to persist in this environment, including counter-defense systems, auxiliary metabolic genes, reverse transcriptases, and alternative genetic code assignments. Extremely low nucleotide diversity is observed in cold seep viral populations, being influenced by factors including microbial host, sediment depth, and cold seep stage. Most cold seep viral genes are under strong purifying selection with trajectories that differ depending on whether cold seeps are active or extinct. This work sheds light on the understanding of environmental adaptation mechanisms and evolutionary patterns of viruses in the sub-seafloor biosphere.

Bacteriophage communities are a reservoir of unexplored microbial diversity in neonatal health and disease

Acquisition and development of the gut microbiome are vital for immune education in neonates, especially those born preterm. As such, microbial communities have been extensively studied in the context of postnatal health and disease. Bacterial communities have been the focus of research in this area due to the relative ease of targeted bacterial sequencing and the availability of databases to align and validate sequencing data. Recent increases in high-throughput metagenomic sequencing accessibility have facilitated research to investigate bacteriophages within the context of neonatal gut microbial communities. Focusing on unexplored viral diversity, has identified novel bacteriophage species and previously uncharacterised viral diversity. In doing so, studies have highlighted links between bacteriophages and bacterial community structure in the context of health and disease. However, much remains unknown about the complex relationships between bacteriophages, the bacteria they infect and their human host. With a particular focus on preterm infants, this review highlights opportunities to explore the influence of bacteriophages on developing microbial communities and the tripartite relationships between bacteriophages, bacteria and the neonatal human host. We suggest a focus on expanding collections of isolated bacteriophages that will further our understanding of the growing numbers of bacteriophages identified in metagenomes.

Phage-bacteria dynamics during the first years of life revealed by trans-kingdom marker gene analysis

Humans are colonized with commensal bacteria soon after birth, and, while this colonization is affected by lifestyle and other factors, bacterial colonization proceeds through well-studied phases. However, less is known about phage communities in early human development due to small study sizes, inability to leverage large databases, and lack of appropriate bioinformatics tools. In this study, whole genome shotgun sequencing data from the TEDDY study, composed of 12,262 longitudinal samples from 887 children in 4 countries, is reanalyzed to assess phage and bacterial dynamics simultaneously. Reads from these samples were mapped to marker genes from both bacteria and a new database of tens of thousands of phage taxa from human microbiomes. We uncover that each child is colonized by hundreds of different phages during the early years, and phages are more transitory than bacteria. Participants' samples continually harbor new phage species over time whereas the diversification of bacterial species begins to saturate. Phage data improves the ability for machine learning models to discriminate samples by country. Finally, while phage populations were individual-specific, striking patterns arose from the larger dataset, showing clear trends of ecological succession amongst phages, which correlated well with putative host bacteria. Improved understanding of phage-bacterial relationships may reveal new means by which to shape and modulate the microbiome and its constituents to improve health and reduce disease, particularly in vulnerable populations where antibiotic use and/or other more drastic measures may not be advised.

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.

Probiotics alleviate constipation and inflammation in late gestating and lactating sows

Constipation and systemic inflammation are common in late pregnant and lactating sows, which cause health problems like uteritis, mastitis, dystocia, or even stillbirth, further influencing piglets' survival and growth. Probiotic supplementation can improve such issues, but the beneficial mechanism of relieving constipation and enhancing gut motility remains underexplored. This study aimed to investigate the effects and mechanism of probiotic supplementation in drinking water to late pregnant sows on constipation, inflammation, and piglets' growth performance. Seventy-four sows were randomly allocated to probiotic (n = 36) and control (n = 38) groups. Probiotic treatment significantly relieved sow constipation, enhanced serum IL-4 and IL-10 levels while reducing serum IL-1β, IL-12p40, and TNF-α levels, and increased piglet daily gain and weaning weight. Furthermore, probiotic administration reshaped the sow gut bacteriome and phageome structure/diversity, accompanied by increases in some potentially beneficial bacteria. At 113 days of gestation, the probiotic group was enriched in several gut microbial bioactive metabolites, multiple carbohydrate-active enzymes that degrade pectin and starch, fecal butyrate and acetate, and some serum metabolites involved in vitamin and amino acid metabolism. Our integrated correlation network analysis revealed that the alleviation of constipation and inflammation was associated with changes in the sow gut bacteriome, phageome, bioactive metabolic potential, and metabolism.

HiFun: homology independent protein function prediction by a novel protein-language self-attention model

Protein function prediction based on amino acid sequence alone is an extremely challenging but important task, especially in metagenomics/metatranscriptomics field, in which novel proteins have been uncovered exponentially from new microorganisms. Many of them are extremely low homology to known proteins and cannot be annotated with homology-based or information integrative methods. To overcome this problem, we proposed a Homology Independent protein Function annotation method (HiFun) based on a unified deep-learning model by reassembling the sequence as protein language. The robustness of HiFun was evaluated using the benchmark datasets and metrics in the CAFA3 challenge. To navigate the utility of HiFun, we annotated 2 212 663 unknown proteins and discovered novel motifs in the UHGP-50 catalog. We proved that HiFun can extract latent function related structure features which empowers it ability to achieve function annotation for non-homology proteins. HiFun can substantially improve newly proteins annotation and expand our understanding of microorganisms' adaptation in various ecological niches. Moreover, we provided a free and accessible webservice at http://www.unimd.org/HiFun, requiring only protein sequences as input, offering researchers an efficient and practical platform for predicting protein functions.

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.

Phables: from fragmented assemblies to high-quality bacteriophage genomes

Microbial communities influence both human health and different environments. Viruses infecting bacteria, known as bacteriophages or phages, play a key role in modulating bacterial communities within environments. High-quality phage genome sequences are essential for advancing our understanding of phage biology, enabling comparative genomics studies, and developing phage-based diagnostic tools. Most available viral identification tools consider individual sequences to determine whether they are of viral origin. As a result of the challenges in viral assembly, fragmentation of genomes can occur, leading to the need for new approaches in viral identification. Therefore, the identification and characterisation of novel phages remain a challenge. We introduce Phables, a new computational method to resolve phage genomes from fragmented viral metagenome assemblies. Phables identifies phage-like components in the assembly graph, models each component as a flow network, and uses graph algorithms and flow decomposition techniques to identify genomic paths. Experimental results of viral metagenomic samples obtained from different environments show that Phables recovers on average over 49% more high-quality phage genomes compared to existing viral identification tools. Furthermore, Phables can resolve variant phage genomes with over 99% average nucleotide identity, a distinction that existing tools are unable to make. Phables is available on GitHub at https://github.com/Vini2/phables.

An extended catalog of integrated prophages in the infant and adult fecal microbiome shows high prevalence of lysogeny

Background and aims

The acquisition and gradual maturation of gut microbial communities during early childhood is central to an individual's healthy development. Bacteriophages have the potential to shape the gut bacterial communities. However, the complex ecological interactions between phages and their bacterial host are still poorly characterized. In this study, we investigated the abundance and diversity of integrated prophages in infant and adult gut bacteria by detecting integrated prophages in metagenome assembled genomes (MAGs) of commensal bacteria.

Methods

Our study included 88 infants sampled at 3 weeks, 3 months, 6 months, and 12 months (n = 323 total samples), and their parents around delivery time (n = 138 total samples). Fecal DNA was extracted and characterized by using shotgun metagenomic sequencing, and a collection of prokaryotic MAGs was generated. The MAG collection was screened for the presence of integrated bacteriophage sequences, allowing their taxonomic and functional characterization.

Results

A large collection of 6,186 MAGs from infant and adult gut microbiota was obtained and screened for integrated prophages, allowing the identification of 7,165 prophage sequences longer than 10 kb. Strikingly, more than 70% of the near-complete MAGs were identified as lysogens. The prevalence of prophages in MAGs varied across bacterial families, with a lower prevalence observed among Coriobacteriaceae, Eggerthellaceae, Veillonellaceae and Burkholderiaceae, while a very high prevalence of lysogen MAGs were observed in Oscillospiraceae, Enterococcaceae, and Enterobacteriaceae. Interestingly for several bacterial families such as Bifidobacteriaceae and Bacteroidaceae, the prevalence of prophages in MAGs was higher in early infant time point (3 weeks and 3 months) than in later sampling points (6 and 12 months) and in adults. The prophage sequences were clustered into 5,616 species-like vOTUs, 77% of which were novel. Finally, we explored the functional repertoire of the potential auxiliary metabolic genes carried by these prophages, encoding functions involved in carbohydrate metabolism and degradation, amino acid metabolism and carbon metabolism.

Conclusion

Our study provides an enhanced understanding of the diversity and prevalence of lysogens in infant and adult gut microbiota and suggests a complex interplay between prophages and their bacterial hosts.

Long-Read Sequencing Reveals Extensive DNA Methylations in Human Gut Phagenome Contributed by Prevalently Phage-Encoded Methyltransferases

DNA methylation plays a crucial role in the survival of bacteriophages (phages), yet the understanding of their genome methylation remains limited. In this study, DNA methylation patterns are analyzed in 8848 metagenome-assembled high-quality phages from 104 fecal samples using single-molecule real-time sequencing. The results demonstrate that 97.60% of gut phages exhibit methylation, with certain factors correlating with methylation densities. Phages with higher methylation densities appear to have potential viability advantages. Strikingly, more than one-third of the phages possess their own DNA methyltransferases (MTases). Increased MTase copies are associated with higher genome methylation densities, specific methylation motifs, and elevated prevalence of certain phage groups. Notably, the majority of these MTases share close homology with those encoded by gut bacteria, suggesting their exchange during phage-bacterium interactions. Furthermore, these MTases can be employed to accurately predict phage-host relationships. Overall, the findings indicate the widespread utilization of DNA methylation by gut DNA phages as an evasion mechanism against host defense systems, with a substantial contribution from phage-encoded MTases.

Interrogating the viral dark matter of the rumen ecosystem with a global virome database

The diverse rumen virome can modulate the rumen microbiome, but it remains largely unexplored. Here, we mine 975 published rumen metagenomes for viral sequences, create a global rumen virome database (RVD), and analyze the rumen virome for diversity, virus-host linkages, and potential roles in affecting rumen functions. Containing 397,180 species-level viral operational taxonomic units (vOTUs), RVD substantially increases the detection rate of rumen viruses from metagenomes compared with IMG/VR V3. Most of the classified vOTUs belong to Caudovirales, differing from those found in the human gut. The rumen virome is predicted to infect the core rumen microbiome, including fiber degraders and methanogens, carries diverse auxiliary metabolic genes, and thus likely impacts the rumen ecosystem in both a top-down and a bottom-up manner. RVD and the findings provide useful resources and a baseline framework for future research to investigate how viruses may impact the rumen ecosystem and digestive physiology.

Terabase-Scale Coassembly of a Tropical Soil Microbiome

Petabases of environmental metagenomic data are publicly available, presenting an opportunity to characterize complex environments and discover novel lineages of life. Metagenome coassembly, in which many metagenomic samples from an environment are simultaneously analyzed to infer the underlying genomes' sequences, is an essential tool for achieving this goal. We applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 terabases (Tbp) of metagenome data from a tropical soil in the Luquillo Experimental Forest (LEF), Puerto Rico. The resulting coassembly yielded 39 high-quality (>90% complete, <5% contaminated, with predicted 23S, 16S, and 5S rRNA genes and ≥18 tRNAs) metagenome-assembled genomes (MAGs), including two from the candidate phylum Eremiobacterota. Another 268 medium-quality (≥50% complete, <10% contaminated) MAGs were extracted, including the candidate phyla Dependentiae, Dormibacterota, and Methylomirabilota. In total, 307 medium- or higher-quality MAGs were assigned to 23 phyla, compared to 294 MAGs assigned to nine phyla in the same samples individually assembled. The low-quality (<50% complete, <10% contaminated) MAGs from the coassembly revealed a 49% complete rare biosphere microbe from the candidate phylum FCPU426 among other low-abundance microbes, an 81% complete fungal genome from the phylum Ascomycota, and 30 partial eukaryotic MAGs with ≥10% completeness, possibly representing protist lineages. A total of 22,254 viruses, many of them low abundance, were identified. Estimation of metagenome coverage and diversity indicates that we may have characterized ≥87.5% of the sequence diversity in this humid tropical soil and indicates the value of future terabase-scale sequencing and coassembly of complex environments. IMPORTANCE Petabases of reads are being produced by environmental metagenome sequencing. An essential step in analyzing these data is metagenome assembly, the computational reconstruction of genome sequences from microbial communities. "Coassembly" of metagenomic sequence data, in which multiple samples are assembled together, enables more complete detection of microbial genomes in an environment than "multiassembly," in which samples are assembled individually. To demonstrate the potential for coassembling terabases of metagenome data to drive biological discovery, we applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 Tbp of reads from a humid tropical soil environment. The resulting coassembly, its functional annotation, and analysis are presented here. The coassembly yielded more, and phylogenetically more diverse, microbial, eukaryotic, and viral genomes than the multiassembly of the same data. Our resource may facilitate the discovery of novel microbial biology in tropical soils and demonstrates the value of terabase-scale metagenome sequencing.

Maast: genotyping thousands of microbial strains efficiently

Existing single nucleotide polymorphism (SNP) genotyping algorithms do not scale for species with thousands of sequenced strains, nor do they account for conspecific redundancy. Here we present a bioinformatics tool, Maast, which empowers population genetic meta-analysis of microbes at an unrivaled scale. Maast implements a novel algorithm to heuristically identify a minimal set of diverse conspecific genomes, then constructs a reliable SNP panel for each species, and enables rapid and accurate genotyping using a hybrid of whole-genome alignment and k-mer exact matching. We demonstrate Maast's utility by genotyping thousands of Helicobacter pylori strains and tracking SARS-CoV-2 diversification.

Presence and role of viruses in anaerobic digestion of food waste under environmental variability

BACKGROUND

The interaction among microorganisms in the anaerobic digestion of food waste (ADFW) reactors lead to the degradation of organics and the recycling of energy. Viruses are an important component of the microorganisms involved in ADFW, but are rarely investigated. Furthermore, little is known about how viruses affect methanogenesis.

RESULTS

Thousands of viral sequences were recovered from five full-scale ADFW reactors. Gene-sharing networks indicated that the ADFW samples contained substantial numbers of unexplored anaerobic-specific viruses. Moreover, the viral communities in five full-scale reactors exhibited both commonalities and heterogeneities. The lab-scale dynamic analysis of typical ADFW scenarios suggested that the viruses had similar kinetic characteristics to their prokaryotic hosts. By associating with putative hosts, a majority of the bacteria and archaea phyla were found to be infected by viruses. Viruses may influence prokaryotic ecological niches, and thus methanogenesis, by infecting key functional microorganisms, such as sulfate-reducing bacteria (SRB), syntrophic acetate-oxidizing bacteria (SAOB), and methanogens. Metabolic predictions for the viruses suggested that they may collaborate with hosts at key steps of sulfur and long-chain fatty acid (LCFA) metabolism and could be involved in typical methanogenesis pathways to participate in methane production.

CONCLUSIONS

Our results expanded the diversity of viruses in ADFW systems and suggested two ways that viral manipulated ADFW biochemical processes. Video Abstract.

Phage diversity in cell-free DNA identifies bacterial pathogens in human sepsis cases

Bacteriophages, viruses that infect bacteria, have great specificity for their bacterial hosts at the strain and species level. However, the relationship between the phageome and associated bacterial population dynamics is unclear. Here we generated a computational pipeline to identify sequences associated with bacteriophages and their bacterial hosts in cell-free DNA from plasma samples. Analysis of two independent cohorts, including a Stanford Cohort of 61 septic patients and 10 controls and the SeqStudy cohort of 224 septic patients and 167 controls, reveals a circulating phageome in the plasma of all sampled individuals. Moreover, infection is associated with overrepresentation of pathogen-specific phages, allowing for identification of bacterial pathogens. We find that information on phage diversity enables identification of the bacteria that produced these phages, including pathovariant strains of Escherichia coli. Phage sequences can likewise be used to distinguish between closely related bacterial species such as Staphylococcus aureus, a frequent pathogen, and coagulase-negative Staphylococcus, a frequent contaminant. Phage cell-free DNA may have utility in studying bacterial infections.

Small protein modules dictate prophage fates during polylysogeny

Most bacteria in the biosphere are predicted to be polylysogens harbouring multiple prophages1-5. In studied systems, prophage induction from lysogeny to lysis is near-universally driven by DNA-damaging agents6. Thus, how co-residing prophages compete for cell resources if they respond to an identical trigger is unknown. Here we discover regulatory modules that control prophage induction independently of the DNA-damage cue. The modules bear little resemblance at the sequence level but share a regulatory logic by having a transcription factor that activates the expression of a neighbouring gene that encodes a small protein. The small protein inactivates the master repressor of lysis, which leads to induction. Polylysogens that harbour two prophages exposed to DNA damage release mixed populations of phages. Single-cell analyses reveal that this blend is a consequence of discrete subsets of cells producing one, the other or both phages. By contrast, induction through the DNA-damage-independent module results in cells producing only the phage sensitive to that specific cue. Thus, in the polylysogens tested, the stimulus used to induce lysis determines phage productivity. Considering the lack of potent DNA-damaging agents in natural habitats, additional phage-encoded sensory pathways to lysis likely have fundamental roles in phage-host biology and inter-prophage competition.

Petascale Homology Search for Structure Prediction

The recent CASP15 competition highlighted the critical role of multiple sequence alignments (MSAs) in protein structure prediction, as demonstrated by the success of the top AlphaFold2-based prediction methods. To push the boundaries of MSA utilization, we conducted a petabase-scale search of the Sequence Read Archive (SRA), resulting in gigabytes of aligned homologs for CASP15 targets. These were merged with default MSAs produced by ColabFold-search and provided to ColabFold-predict. By using SRA data, we achieved highly accurate predictions (GDT_TS > 70) for 66% of the non-easy targets, whereas using ColabFold-search default MSAs scored highly in only 52%. Next, we tested the effect of deep homology search and ColabFold's advanced features, such as more recycles, on prediction accuracy. While SRA homologs were most significant for improving ColabFold's CASP15 ranking from 11th to 3rd place, other strategies contributed too. We analyze these in the context of existing strategies to improve prediction.

Ultra-deep sequencing of Hadza hunter-gatherers recovers vanishing gut microbes

The gut microbiome modulates immune and metabolic health. Human microbiome data are biased toward industrialized populations, limiting our understanding of non-industrialized microbiomes. Here, we performed ultra-deep metagenomic sequencing on 351 fecal samples from the Hadza hunter-gatherers of Tanzania and comparative populations in Nepal and California. We recovered 91,662 genomes of bacteria, archaea, bacteriophages, and eukaryotes, 44% of which are absent from existing unified datasets. We identified 124 gut-resident species vanishing in industrialized populations and highlighted distinct aspects of the Hadza gut microbiome related to in situ replication rates, signatures of selection, and strain sharing. Industrialized gut microbes were found to be enriched in genes associated with oxidative stress, possibly a result of microbiome adaptation to inflammatory processes. This unparalleled view of the Hadza gut microbiome provides a valuable resource, expands our understanding of microbes capable of colonizing the human gut, and clarifies the extensive perturbation induced by the industrialized lifestyle.

Virus impacted community adaptation in oligotrophic groundwater environment revealed by Hi-C coupled metagenomic and viromic study

Viruses play a crucial role in microbial mortality, diversity and biogeochemical cycles. Groundwater is the largest global freshwater and one of the most oligotrophic aquatic systems on Earth, but how microbial and viral communities are shaped in this special habitat is largely unexplored. In this study, we collected groundwater samples from 23 to 60 m aquifers at Yinchuan Plain, China. In total, 1920 non-reductant viral contigs were retrieved from metagenomes and viromes constructed by Illumina and Nanopore hybrid sequencing. Only 3% of them could be clustered with known viruses, most of which were Caudoviricetes. Coupling 1.2 Tb Hi-C sequencing with CRISPR matching and homology search, we connected 469 viruses with their hosts while some viral clusters presented a broad-host-range trait. Meanwhile, a large proportion of biosynthesis related auxiliary metabolism genes were identified. Those characteristics might benefit viruses for a better survival in this special oligotrophic environment. Additionally, the groundwater virome showed genomic features distinct from those of the open ocean and wastewater treatment facilities in GC distribution and unannotated gene compositions. This paper expands the current knowledge of the global viromic records and serves as a foundation for a more thorough understanding of viruses in groundwater.

Enrichment Culture but Not Metagenomic Sequencing Identified a Highly Prevalent Phage Infecting Lactiplantibacillus plantarum in Human Feces

Lactiplantibacillus plantarum (previously known as Lactobacillus plantarum) is increasingly used as a probiotic to treat human diseases, but its phages in the human gut remain unexplored. Here, we report its first gut phage, Gut-P1, which we systematically screened using metagenomic sequencing, virus-like particle (VLP) sequencing, and enrichment culture from 35 fecal samples. Gut-P1 is virulent, belongs to the Douglaswolinvirus genus, and is highly prevalent in the gut (~11% prevalence); it has a genome of 79,928 bp consisting of 125 protein coding genes and displaying low sequence similarities to public L. plantarum phages. Physiochemical characterization shows that it has a short latent period and adapts to broad ranges of temperatures and pHs. Furthermore, Gut-P1 strongly inhibits the growth of L. plantarum strains at a multiplicity of infection (MOI) of 1e-6. Together, these results indicate that Gut-P1 can greatly impede the application of L. plantarum in humans. Strikingly, Gut-P1 was identified only in the enrichment culture, not in our metagenomic or VLP sequencing data nor in any public human phage databases, indicating the inefficiency of bulk sequencing in recovering low-abundance but highly prevalent phages and pointing to the unexplored hidden diversity of the human gut virome despite recent large-scale sequencing and bioinformatics efforts. IMPORTANCE As Lactiplantibacillus plantarum (previously known as Lactobacillus plantarum) is increasingly used as a probiotic to treat human gut-related diseases, its bacteriophages may pose a certain threat to their further application and should be identified and characterized more often from the human intestine. Here, we isolated and identified the first gut L. plantarum phage that is prevalent in a Chinese population. This phage, Gut-P1, is virulent and can strongly inhibit the growth of multiple L. plantarum strains at low MOIs. Our results also show that bulk sequencing is inefficient at recovering low-abundance but highly prevalent phages such as Gut-P1, suggesting that the hidden diversity of human enteroviruses has not yet been explored. Our results call for innovative approaches to isolate and identify intestinal phages from the human gut and to rethink our current understanding of the enterovirus, particularly its underestimated diversity and overestimated individual specificity.

Longitudinal Investigation of Enteric Virome Signatures from Parental-Generation to Offspring Pigs

To date, studies on the swine gut microbiome have focused almost exclusively on bacteria. Despite recent advances in the understanding of the swine gut bacteriome at different growth stages, a comprehensive longitudinal study of the lifetime dynamics of the swine gut virome is lacking. Here, we used metagenomic sequencing combined with bioinformatic analysis techniques to characterize the gut viromes of parental-generation and offspring pigs at different biological classification levels. We collected 54 fecal samples from 36 parental-generation pigs (18 breeding boars [Duroc] and 18 pregnant/lactating sows [Landrace]) and 108 fecal samples from 18 offspring pigs during the lactation (day 3), nursery (days 26, 35, and 49), growing (day 120), and finishing (day 180) stages. Alpha diversity, including community richness (richness index) and diversity (Shannon index), showed an overall increasing trend in offspring pigs. Distinct shifts (beta diversity) in the microbiome structure along different growth stages were observed. The linear discriminant analysis effect size (LEfSe) algorithm revealed 53 viral genus that are stage specific. Host prediction results showed that enteric viruses are probably correlated with carbohydrate decomposition. We identified abundant auxiliary carbohydrate-active enzyme (CAZyme) genes from enteric viruses, most of which are glycoside hydrolase genes and participate in the biolysis of complex polysaccharides. IMPORTANCE This study shows that distinct stage-associated swine gut viromes may be determined by age and/or gut physiology at different growth stages, and enteric viruses probably manipulate carbohydrate decomposition by abundant glycoside hydrolases. These findings fill a gap in the longitudinal pattern of the swine gut virome and lay the foundation for research on the function of swine enteric viruses.

Centenarians have a diverse gut virome with the potential to modulate metabolism and promote healthy lifespan

Distinct gut microbiome ecology may be implicated in the prevention of aging-related diseases as it influences systemic immune function and resistance to infections. Yet, the viral component of the microbiome throughout different stages in life remains unexplored. Here we present a characterization of the centenarian gut virome using previously published metagenomes from 195 individuals from Japan and Sardinia. Compared with gut viromes of younger adults (>18 yr) and older individuals (>60 yr), centenarians had a more diverse virome including previously undescribed viral genera, such as viruses associated with Clostridia. A population shift towards higher lytic activity was also observed. Finally, we investigated phage-encoded auxiliary functions that influence bacterial physiology, which revealed an enrichment of genes supporting key steps in sulfate metabolic pathways. Phage and bacterial members of the centenarian microbiome displayed an increased potential for converting methionine to homocysteine, sulfate to sulfide and taurine to sulfide. A greater metabolic output of microbial hydrogen sulfide in centenarians may in turn support mucosal integrity and resistance to pathobionts.

Geospatially-resolved public-health surveillance via wastewater sequencing

Wastewater, which contains everything from pathogens to pollutants, is a geospatially-and temporally-linked microbial fingerprint of a given population. As a result, it can be leveraged for monitoring multiple dimensions of public health across locales and time. Here, we integrate targeted and bulk RNA sequencing (n=1,419 samples) to track the viral, bacterial, and functional content over geospatially distinct areas within Miami Dade County from 2020-2022. First, we used targeted amplicon sequencing (n=966) to track diverse SARS-CoV-2 variants across space and time, and we found a tight correspondence with clinical caseloads from University students (N = 1,503) and Miami-Dade County hospital patients (N = 3,939 patients), as well as an 8-day earlier detection of the Delta variant in wastewater vs. in patients. Additionally, in 453 metatranscriptomic samples, we demonstrate that different wastewater sampling locations have clinically and public-health-relevant microbiota that vary as a function of the size of the human population they represent. Through assembly, alignment-based, and phylogenetic approaches, we also detect multiple clinically important viruses (e.g., norovirus ) and describe geospatial and temporal variation in microbial functional genes that indicate the presence of pollutants. Moreover, we found distinct profiles of antimicrobial resistance (AMR) genes and virulence factors across campus buildings, dorms, and hospitals, with hospital wastewater containing a significant increase in AMR abundance. Overall, this effort lays the groundwork for systematic characterization of wastewater to improve public health decision making and a broad platform to detect emerging pathogens.

Metagenomics: An Effective Approach for Exploring Microbial Diversity and Functions

Various fields have been identified in the "omics" era, such as genomics, proteomics, transcriptomics, metabolomics, phenomics, and metagenomics. Among these, metagenomics has enabled a significant increase in discoveries related to the microbial world. Newly discovered microbiomes in different ecologies provide meaningful information on the diversity and functions of microorganisms on the Earth. Therefore, the results of metagenomic studies have enabled new microbe-based applications in human health, agriculture, and the food industry, among others. This review summarizes the fundamental procedures on recent advances in bioinformatic tools. It also explores up-to-date applications of metagenomics in human health, food study, plant research, environmental sciences, and other fields. Finally, metagenomics is a powerful tool for studying the microbial world, and it still has numerous applications that are currently hidden and awaiting discovery. Therefore, this review also discusses the future perspectives of metagenomics.

The human gut virome: composition, colonization, interactions, and impacts on human health

The gut virome is an incredibly complex part of the gut ecosystem. Gut viruses play a role in many disease states, but it is unknown to what extent the gut virome impacts everyday human health. New experimental and bioinformatic approaches are required to address this knowledge gap. Gut virome colonization begins at birth and is considered unique and stable in adulthood. The stable virome is highly specific to each individual and is modulated by varying factors such as age, diet, disease state, and use of antibiotics. The gut virome primarily comprises bacteriophages, predominantly order Crassvirales, also referred to as crAss-like phages, in industrialized populations and other Caudoviricetes (formerly Caudovirales). The stability of the virome's regular constituents is disrupted by disease. Transferring the fecal microbiome, including its viruses, from a healthy individual can restore the functionality of the gut. It can alleviate symptoms of chronic illnesses such as colitis caused by Clostridiodes difficile. Investigation of the virome is a relatively novel field, with new genetic sequences being published at an increasing rate. A large percentage of unknown sequences, termed 'viral dark matter', is one of the significant challenges facing virologists and bioinformaticians. To address this challenge, strategies include mining publicly available viral datasets, untargeted metagenomic approaches, and utilizing cutting-edge bioinformatic tools to quantify and classify viral species. Here, we review the literature surrounding the gut virome, its establishment, its impact on human health, the methods used to investigate it, and the viral dark matter veiling our understanding of the gut virome.

Exposing new taxonomic variation with inflammation - a murine model-specific genome database for gut microbiome researchers

BACKGROUND

The murine CBA/J mouse model widely supports immunology and enteric pathogen research. This model has illuminated Salmonella interactions with the gut microbiome since pathogen proliferation does not require disruptive pretreatment of the native microbiota, nor does it become systemic, thereby representing an analog to gastroenteritis disease progression in humans. Despite the value to broad research communities, microbiota in CBA/J mice are not represented in current murine microbiome genome catalogs.

RESULTS

Here we present the first microbial and viral genomic catalog of the CBA/J murine gut microbiome. Using fecal microbial communities from untreated and Salmonella-infected, highly inflamed mice, we performed genomic reconstruction to determine the impacts on gut microbiome membership and functional potential. From high depth whole community sequencing (~ 42.4 Gbps/sample), we reconstructed 2281 bacterial and 4516 viral draft genomes. Salmonella challenge significantly altered gut membership in CBA/J mice, revealing 30 genera and 98 species that were conditionally rare and unsampled in non-inflamed mice. Additionally, inflamed communities were depleted in microbial genes that modulate host anti-inflammatory pathways and enriched in genes for respiratory energy generation. Our findings suggest decreases in butyrate concentrations during Salmonella infection corresponded to reductions in the relative abundance in members of the Alistipes. Strain-level comparison of CBA/J microbial genomes to prominent murine gut microbiome databases identified newly sampled lineages in this resource, while comparisons to human gut microbiomes extended the host relevance of dominant CBA/J inflammation-resistant strains.

CONCLUSIONS

This CBA/J microbiome database provides the first genomic sampling of relevant, uncultivated microorganisms within the gut from this widely used laboratory model. Using this resource, we curated a functional, strain-resolved view on how Salmonella remodels intact murine gut communities, advancing pathobiome understanding beyond inferences from prior amplicon-based approaches. Salmonella-induced inflammation suppressed Alistipes and other dominant members, while rarer commensals like Lactobacillus and Enterococcus endure. The rare and novel species sampled across this inflammation gradient advance the utility of this microbiome resource to benefit the broad research needs of the CBA/J scientific community, and those using murine models for understanding the impact of inflammation on the gut microbiome more generally. Video Abstract.

Candida expansion in the gut of lung cancer patients associates with an ecological signature that supports growth under dysbiotic conditions

Candida species overgrowth in the human gut is considered a prerequisite for invasive candidiasis, but our understanding of gut bacteria promoting or restricting this overgrowth is still limited. By integrating cross-sectional mycobiome and shotgun metagenomics data from the stool of 75 male and female cancer patients at risk but without systemic candidiasis, bacterial communities in high Candida samples display higher metabolic flexibility yet lower contributional diversity than those in low Candida samples. We develop machine learning models that use only bacterial taxa or functional relative abundances to predict the levels of Candida genus and species in an external validation cohort with an AUC of 78.6-81.1%. We propose a mechanism for intestinal Candida overgrowth based on an increase in lactate-producing bacteria, which coincides with a decrease in bacteria that regulate short chain fatty acid and oxygen levels. Under these conditions, the ability of Candida to harness lactate as a nutrient source may enable Candida to outcompete other fungi in the gut.

Highly host-linked viromes in the built environment possess habitat-dependent diversity and functions for potential virus-host coevolution

Viruses in built environments (BEs) raise public health concerns, yet they are generally less studied than bacteria. To better understand viral dynamics in BEs, this study assesses viromes from 11 habitats across four types of BEs with low to high occupancy. The diversity, composition, metabolic functions, and lifestyles of the viromes are found to be habitat dependent. Caudoviricetes species are ubiquitous on surface habitats in the BEs, and some of them are distinct from those present in other environments. Antimicrobial resistance genes are identified in viruses inhabiting surfaces frequently touched by occupants and in viruses inhabiting occupants' skin. Diverse CRISPR/Cas immunity systems and anti-CRISPR proteins are found in bacterial hosts and viruses, respectively, consistent with the strongly coupled virus-host links. Evidence of viruses potentially aiding host adaptation in a specific-habitat manner is identified through a unique gene insertion. This work illustrates that virus-host interactions occur frequently in BEs and that viruses are integral members of BE microbiomes.

Large language models improve annotation of viral proteins

Viral sequences are poorly annotated in environmental samples, a major roadblock to understanding how viruses influence microbial community structure. Current annotation approaches rely on alignment-based sequence ho-mology methods, which are limited by available viral sequences and sequence divergence in viral proteins. Here, we show that protein language model representations capture viral protein function beyond the limits of remote sequence homology by targeting two axes of viral sequence annotation: systematic labeling of protein families and function identification for biologic discovery. Protein language model representations capture protein functional properties specific to viruses and expand the annotated fraction of ocean virome viral protein sequences by 37%. Among unannotated viral protein families, we identify a novel DNA editing protein family that defines a new mobile element in marine picocyanobacteria. Protein language models thus significantly enhance remote homology detection of viral proteins and can be utilized to enable new biological discovery across diverse functional categories.

Structural atlas of a human gut crassvirus

CrAssphage and related viruses of the order Crassvirales (hereafter referred to as crassviruses) were originally discovered by cross-assembly of metagenomic sequences. They are the most abundant viruses in the human gut, are found in the majority of individual gut viromes, and account for up to 95% of the viral sequences in some individuals1-4. Crassviruses are likely to have major roles in shaping the composition and functionality of the human microbiome, but the structures and roles of most of the virally encoded proteins are unknown, with only generic predictions resulting from bioinformatic analyses4,5. Here we present a cryo-electron microscopy reconstruction of Bacteroides intestinalis virus ΦcrAss0016, providing the structural basis for the functional assignment of most of its virion proteins. The muzzle protein forms an assembly about 1 MDa in size at the end of the tail and exhibits a previously unknown fold that we designate the 'crass fold', that is likely to serve as a gatekeeper that controls the ejection of cargos. In addition to packing the approximately 103 kb of virus DNA, the ΦcrAss001 virion has extensive storage space for virally encoded cargo proteins in the capsid and, unusually, within the tail. One of the cargo proteins is present in both the capsid and the tail, suggesting a general mechanism for protein ejection, which involves partial unfolding of proteins during their extrusion through the tail. These findings provide a structural basis for understanding the mechanisms of assembly and infection of these highly abundant crassviruses.

iPHoP: An integrated machine learning framework to maximize host prediction for metagenome-derived viruses of archaea and bacteria

The extraordinary diversity of viruses infecting bacteria and archaea is now primarily studied through metagenomics. While metagenomes enable high-throughput exploration of the viral sequence space, metagenome-derived sequences lack key information compared to isolated viruses, in particular host association. Different computational approaches are available to predict the host(s) of uncultivated viruses based on their genome sequences, but thus far individual approaches are limited either in precision or in recall, i.e., for a number of viruses they yield erroneous predictions or no prediction at all. Here, we describe iPHoP, a two-step framework that integrates multiple methods to reliably predict host taxonomy at the genus rank for a broad range of viruses infecting bacteria and archaea, while retaining a low false discovery rate. Based on a large dataset of metagenome-derived virus genomes from the IMG/VR database, we illustrate how iPHoP can provide extensive host prediction and guide further characterization of uncultivated viruses.

Large-scale phage cultivation for commensal human gut bacteria

Phages are highly abundant in the human gut, yet most of them remain uncultured. Here, we present a gut phage isolate collection (GPIC) containing 209 phages for 42 commensal human gut bacterial species. Genome analysis of the phages identified 34 undescribed genera. We discovered 22 phages from the Salasmaviridae family that have small genomes (∼10-20 kbp) and infect Gram-positive bacteria. Two phages from a candidate family, Paboviridae, with high prevalence in the human gut were also identified. Infection assays showed that Bacteroides and Parabacteroides phages are specific to a bacterial species, and strains of the same species also exhibit substantial variations in phage susceptibility. A cocktail of 8 phages with a broad host range for Bacteroides fragilis strains effectively reduced their abundance in complex host-derived communities in vitro. Our study expands the diversity of cultured human gut bacterial phages and provides a valuable resource for human microbiome engineering.

Expanding known viral diversity in the healthy infant gut

The gut microbiome is shaped through infancy and impacts the maturation of the immune system, thus protecting against chronic disease later in life. Phages, or viruses that infect bacteria, modulate bacterial growth by lysis and lysogeny, with the latter being especially prominent in the infant gut. Viral metagenomes (viromes) are difficult to analyse because they span uncharted viral diversity, lacking marker genes and standardized detection methods. Here we systematically resolved the viral diversity in faecal viromes from 647 1-year-olds belonging to Copenhagen Prospective Studies on Asthma in Childhood 2010, an unselected Danish cohort of healthy mother-child pairs. By assembly and curation we uncovered 10,000 viral species from 248 virus family-level clades (VFCs). Most (232 VFCs) were previously unknown, belonging to the Caudoviricetes viral class. Hosts were determined for 79% of phage using clustered regularly interspaced short palindromic repeat spacers within bacterial metagenomes from the same children. Typical Bacteroides-infecting crAssphages were outnumbered by undescribed phage families infecting Clostridiales and Bifidobacterium. Phage lifestyles were conserved at the viral family level, with 33 virulent and 118 temperate phage families. Virulent phages were more abundant, while temperate ones were more prevalent and diverse. Together, the viral families found in this study expand existing phage taxonomy and provide a resource aiding future infant gut virome research.

Effect of the probiotic strain, Lactiplantibacillus plantarum P9, on chronic constipation: a randomized, double-blind, placebo-controlled study

Chronic constipation (CC) is a common gastrointestinal condition associated with intestinal inflammation, and the condition considerably impairs patients' quality of life. We conducted a large-scale 42-day randomized, double-blind, placebo-controlled trial to investigate the effect of probiotics in alleviating CC. 163 patients diagnosed with CC (following Rome IV criteria) were randomly divided into probiotic (n = 78; received Lactiplantibacillus plantarum P9 [P9]; 1×10¹¹ CFU/day) and placebo (n = 85; received placebo material) groups. Ingesting P9 significantly improved the weekly mean frequency of complete spontaneous bowel movements (CSBMs) and spontaneous bowel movements (SBMs), while significantly reducing the level of worries and concerns (WO; P < 0.05). Comparing with the placebo group, P9 group was significantly enriched in potentially beneficial bacteria (Lactiplantibacillus plantarum and Ruminococcus_B gnavus), while depriving of several bacterial and phage taxa (Oscillospiraceae sp., Lachnospiraceae sp., and Herelleviridae; P < 0.05). Interesting significant correlations were also observed between some clinical parameters and subjects' gut microbiome, including: negative correlation between Oscillospiraceae sp. and SBMs; positive correlation between WO and Oscillospiraceae sp., Lachnospiraceae sp. Additionally, P9 group had significantly (P < 0.05) more predicted gut microbial bioactive potential involved in the metabolism of amino acids (L-asparagine, L-pipecolinic), short-/medium-chain fatty acids (valeric acid and caprylic acid). Furthermore, several metabolites (p-cresol, methylamine, trimethylamine) related to the intestinal barrier and transit decreased significantly after P9 administration (P < 0.05). In short, the constipation relief effect of P9 intervention was accompanied by desirable changes in the fecal metagenome and metabolome. Our findings support the notion of applying probiotics in managing CC.

The genomic landscape of reference genomes of cultivated human gut bacteria

Culture-independent metagenomic studies have revolutionized our understanding of the gut microbiota. However, the lack of full genomes from cultured species is still a limitation for in-depth studies of the gut microbiota. Here we present a substantially expanded version of our Cultivated Genome Reference (CGR), termed CGR2, providing 3324 high-quality draft genomes from isolates selected from a large-scale cultivation of bacterial isolates from fecal samples of healthy Chinese individuals. The CGR2 classifies 527 species (179 previously unidentified species) from 8 phyla, and uncovers a genomic and functional diversity of Collinsella aerofaciens. The CGR2 genomes match 126 metagenome-assembled genomes without cultured representatives in the Unified Human Gastrointestinal Genome (UHGG) collection and harbor 3767 unidentified secondary metabolite biosynthetic gene clusters, providing a source of natural compounds with pharmaceutical potentials. We uncover accurate phage-bacterium linkages providing information on the evolutionary characteristics of interaction between bacteriophages and bacteria at the strain level.

Diversity and potential function of pig gut DNA viruses

Viruses are ubiquitous in the gut of animals and play an important role in the ecology of the gut microbiome. The potential effects of these substances on the growth and development of the body are not fully known. Little is known about the effects of breeding environment on pig gut virome. Here, there are 3584 viral operational taxonomic units (vOTUs) longer than 5 kb identified by virus-enriched metagenome sequencing from 25 pig fecal samples. Only a small minority of vOTUs (11.16%) can be classified at the family level, and ∼50% of the genes could be annotated, supporting the concept of pig gut as reservoirs of substantial undescribed viral genetic diversity. The composition of pig gut virome in the six regions may be related to geography. There are only 20 viral clusters (VCs) shared among pig gut virome in six regions of Shanxi Province. These viruses rarely carry antibiotic resistance genes (ARGs). At the same time, they possess abundant auxiliary metabolic genes (AMGs) potentially involved in carbon, sulfur metabolism and cofactor biosynthesis, etc. This study has revealed the unique characteristics and potential function of pig gut DNA virome and established a foundation for the recognition of the viral roles in gut environment.

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

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

Interpersonal variability of the human gut virome confounds disease signal detection in IBD

Viruses are increasingly recognised as important components of the human microbiome, fulfilling numerous ecological roles including bacterial predation, immune stimulation, genetic diversification, horizontal gene transfer, microbial interactions, and augmentation of metabolic functions. However, our current view of the human gut virome is tainted by previous sequencing requirements that necessitated the amplification of starting nucleic acids. In this study, we performed an original longitudinal analysis of 40 healthy control, 19 Crohn's disease, and 20 ulcerative colitis viromes over three time points without an amplification bias, which revealed and highlighted the interpersonal individuality of the human gut virome. In contrast to a 16 S rRNA gene analysis of matched samples, we show that α- and β-diversity metrics of unamplified viromes are not as efficient at discerning controls from patients with inflammatory bowel disease. Additionally, we explored the intrinsic properties of unamplified gut viromes and show there is considerable interpersonal variability in viral taxa, infrequent longitudinal persistence of intrapersonal viruses, and vast fluctuations in the abundance of temporal viruses. Together, these properties of unamplified faecal viromes confound the ability to discern disease associations but significantly advance toward an unbiased and accurate representation of the human gut virome.

What Lies Beneath? Taking the Plunge into the Murky Waters of Phage Biology

The sequence revolution revealed that bacteria-infecting viruses, known as phages, are Earth's most abundant biological entities. Phages have far-reaching impacts on the form and function of microbial communities and play a fundamental role in ecological processes. However, even well into the sequencing revolution, we have only just begun to explore the murky waters around the phage biology iceberg. Many viral reads cannot be assigned to a culturable isolate, and reference databases are biased toward more easily collectible samples, which likely distorts our conclusions. This minireview points out alternatives to mapping reads to reference databases and highlights innovative bioinformatic and experimental approaches that can help us overcome some of the challenges in phage research and better decipher the impact of phages on microbial communities. Moving beyond the identification of novel phages, we highlight phage metabolomics as an important influencer of bacterial host cell physiology and hope to inspire the reader to consider the effects of phages on host metabolism and ecosystems at large. We encourage researchers to report unassigned/unknown sequencing reads and contigs and to continue developing alternative methods to investigate phages within sequence data.

Tundra Soil Viruses Mediate Responses of Microbial Communities to Climate Warming

The rise of global temperature causes the degradation of the substantial reserves of carbon (C) stored in tundra soils, in which microbial processes play critical roles. Viruses are known to influence the soil C cycle by encoding auxiliary metabolic genes and infecting key microorganisms, but their regulation of microbial communities under climate warming remains unexplored. In this study, we evaluated the responses of viral communities for about 5 years of experimental warming at two depths (15 to 25 cm and 45 to 55 cm) in the Alaskan permafrost region. Our results showed that the viral community and functional gene composition and abundances (including viral functional genes related to replication, structure, infection, and lysis) were significantly influenced by environmental conditions such as total nitrogen (N), total C, and soil thawing duration. Although long-term warming did not impact the viral community composition at the two depths, some glycoside hydrolases encoded by viruses were more abundant at both depths of the warmed plots. With the continuous reduction of total C, viruses may alleviate methane release by altering infection strategies on methanogens. Importantly, viruses can adopt lysogenic and lytic lifestyles to manipulate microbial communities at different soil depths, respectively, which could be one of the major factors causing the differences in microbial responses to warming. This study provides a new ecological perspective on how viruses regulate the responses of microbes to warming at community and functional scales. IMPORTANCE Permafrost thawing causes microbial release of greenhouse gases, exacerbating climate warming. Some previous studies examined the responses of the microbial communities and functions to warming in permafrost region, but the roles of viruses in mediating the responses of microbial communities to warming are poorly understood. This study revealed that warming induced changes in some viral functional classes and in the virus/microbe ratios for specific lineages, which might influence the entire microbial community. Furthermore, differences in viral communities and functions, along with soil depths, are important factors influencing microbial responses to warming. Collectively, our study revealed the regulation of microbial communities by viruses and demonstrated the importance of viruses in the microbial ecology research.

Characterizations of the Gut Bacteriome, Mycobiome, and Virome in Patients with Osteoarthritis

The gut microbiota plays an essential role in the regulation of the immune system and the etiology of human autoimmune diseases. However, a holistic understanding of the gut bacteriome, mycobiome, and virome in patients with osteoarthritis (OA) remains lacking. Here, we explored the gut microbiotas of 44 OA patients and 46 healthy volunteers via deep whole-metagenome shotgun sequencing of their fecal samples. The gut bacteriome and mycobiome were analyzed using a reference-based strategy. Gut viruses were identified from the metagenomic assembled contigs, and the gut virome was profiled based on 6,567 nonredundant viral operational taxonomic units (vOTUs). We revealed that the gut microbiome (including bacteriome, mycobiome, and virome) of OA patients is fundamentally altered, characterized by a panel of 279 differentially abundant bacterial species, 10 fungal species, and 627 vOTUs. The representative OA-enriched bacteria included Anaerostipes hadrus (GENOME147149), Prevotella sp900313215 (GENOME08259), Eubacterium_E hallii (GENOME000299), and Blautia A (GENOME001004), while Bacteroides plebeius A (GENOME239725), Roseburia inulinivorans (GENOME 001770), Dialister sp900343095 (GENOME075103), Phascolarctobacterium faecium (GENOME233517), and several members of Faecalibacterium and Prevotella were depleted in OA patients. Fungi such as Debaryomyces fabryi (GenBank accession no. GCA_003708665), Candida parapsilosis (GCA_000182765), and Apophysomyces trapeziformis (GCA_000696975) were enriched in the OA gut microbiota, and Malassezia restricta (GCA_003290485), Aspergillus fumigatus (GCA_003069565), and Mucor circinelloides (GCA_010203745) were depleted. The OA-depleted viruses spanned Siphoviridae (95 vOTUs), Myoviridae (70 vOTUs), and Microviridae (5 vOTUs), while 30 Siphoviridae vOTUs were enriched in OA patients. Functional analysis of the gut bacteriome and virome also uncovered their functional signatures in relation to OA. Moreover, we demonstrated that the OA-associated gut bacterial and viral signatures are tightly interconnected, suggesting that they may impact disease together. Finally, we showed that the multikingdom signatures are effective in discriminating the OA patients from healthy controls, suggesting the potential of gut microbiota for the prediction of OA and related diseases. Our results delineated the fecal bacteriome, mycobiome, and virome landscapes of the OA microbiota and provided biomarkers that will aid in future mechanistic and clinical intervention studies. IMPORTANCE The gut microbiome of OA patients was completely altered compared to that in healthy individuals, including 279 differentially abundant bacterial species, 10 fungal species and 627 viral operational taxonomic units (vOTUs). Functional analysis of the gut bacteriome and virome also revealed their functional signatures in relation to OA. We found that OA-associated gut bacterial and viral signatures were tightly interconnected, indicating that they may affect the disease together. The OA patients can be discriminated effectively from healthy controls using the multikingdom signatures, suggesting the potential of gut microbiota for the prediction of OA and related diseases.

Evaluation of computational phage detection tools for metagenomic datasets

Introduction

As new computational tools for detecting phage in metagenomes are being rapidly developed, a critical need has emerged to develop systematic benchmarks.

Methods

In this study, we surveyed 19 metagenomic phage detection tools, 9 of which could be installed and run at scale. Those 9 tools were assessed on several benchmark challenges. Fragmented reference genomes are used to assess the effects of fragment length, low viral content, phage taxonomy, robustness to eukaryotic contamination, and computational resource usage. Simulated metagenomes are used to assess the effects of sequencing and assembly quality on the tool performances. Finally, real human gut metagenomes and viromes are used to assess the differences and similarities in the phage communities predicted by the tools.

Results

We find that the various tools yield strikingly different results. Generally, tools that use a homology approach (VirSorter, MARVEL, viralVerify, VIBRANT, and VirSorter2) demonstrate low false positive rates and robustness to eukaryotic contamination. Conversely, tools that use a sequence composition approach (VirFinder, DeepVirFinder, Seeker), and MetaPhinder, have higher sensitivity, including to phages with less representation in reference databases. These differences led to widely differing predicted phage communities in human gut metagenomes, with nearly 80% of contigs being marked as phage by at least one tool and a maximum overlap of 38.8% between any two tools. While the results were more consistent among the tools on viromes, the differences in results were still significant, with a maximum overlap of 60.65%. Discussion: Importantly, the benchmark datasets developed in this study are publicly available and reusable to enable the future comparability of new tools developed.

PhageTailFinder: A tool for phage tail module detection and annotation

Decades of overconsumption of antimicrobials in the treatment and prevention of bacterial infections have resulted in the increasing emergence of drug-resistant bacteria, which poses a significant challenge to public health, driving the urgent need to find alternatives to conventional antibiotics. Bacteriophages are viruses infecting specific bacterial hosts, often destroying the infected bacterial hosts. Phages attach to and enter their potential hosts using their tail proteins, with the composition of the tail determining the range of potentially infected bacteria. To aid the exploitation of bacteriophages for therapeutic purposes, we developed the PhageTailFinder algorithm to predict tail-related proteins and identify the putative tail module in previously uncharacterized phages. The PhageTailFinder relies on a two-state hidden Markov model (HMM) to predict the probability of a given protein being tail-related. The process takes into account the natural modularity of phage tail-related proteins, rather than simply considering amino acid properties or secondary structures for each protein in isolation. The PhageTailFinder exhibited robust predictive power for phage tail proteins in novel phages due to this sequence-independent operation. The performance of the prediction model was evaluated in 13 extensively studied phages and a sample of 992 complete phages from the NCBI database. The algorithm achieved a high true-positive prediction rate (>80%) in over half (571) of the studied phages, and the ROC value was 0.877 using general models and 0.968 using corresponding morphologic models. It is notable that the median ROC value of 992 complete phages is more than 0.75 even for novel phages, indicating the high accuracy and specificity of the PhageTailFinder. When applied to a dataset containing 189,680 viral genomes derived from 11,810 bulk metagenomic human stool samples, the ROC value was 0.895. In addition, tail protein clusters could be identified for further studies by density-based spatial clustering of applications with the noise algorithm (DBSCAN). The developed PhageTailFinder tool can be accessed either as a web server (http://www.microbiome-bigdata.com/PHISDetector/index/tools/PhageTailFinder) or as a stand-alone program on a standard desktop computer (https://github.com/HIT-ImmunologyLab/PhageTailFinder).