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

Deciphering microbial metabolic interactions and their implications for community dynamics in acid mine drainage sediments

The microbially-mediated reduction processes have potential for the bioremediation of acid mine drainage (AMD), which represents a worldwide environment problem. However, we know little about the microbial interactions in anaerobic AMD sediments. Here we utilized genome-resolved metagenomics to uncover the nature of cooperative and competitive metabolic interactions in 90 AMD sediments across Southern China. Our analyses recovered well-represented prokaryotic communities through the reconstruction of 2625 population genomes. Functional analyses of these genomes revealed extensive metabolic handoffs which occurred more frequently in nitrogen metabolism than in sulfur metabolism, as well as stable functional redundancy across sediments resulting from populations with low genomic relatedness. Genome-scale metabolic modeling showed that metabolic competition promoted microbial co-occurrence relationships, suggesting that community assembly was dominated by habitat filtering in sediments. Notably, communities colonizing more extreme conditions tended to be highly competitive, which was typically accompanied with increased network complexity but decreased stability of the microbiome. Finally, our results demonstrated that heterotrophic Thermoplasmatota associated with ferric iron and sulfate reduction contributed most to the elevated levels of competition. Our study shed light on the cooperative and competitive metabolisms of microbiome in the hazardous AMD sediments, which may provide preliminary clues for the AMD bioremediation in the future.

Stable coexistence between an archaeal virus and the dominant methanogen of the human gut

The human gut virome, which is mainly composed of bacteriophages, also includes viruses infecting archaea, yet their role remains poorly understood due to lack of isolates. Here, we characterize a temperate archaeal virus (MSTV1) infecting Methanobrevibacter smithii, the dominant methanogenic archaeon of the human gut. The MSTV1 genome is integrated in the host chromosome as a provirus which is sporadically induced, resulting in virion release. Using cryo-electron tomography, we capture several intracellular virion assembly intermediates and confirm that only a small fraction of the host population actively produces virions in vitro. Similar low frequency of induction is observed in a mouse colonization model, using mice harboring a stable consortium of 12 bacterial species (OMM12). Transcriptomic analysis suggests a regulatory lysogeny-lysis switch involving an interplay between viral proteins to maintain virus-host equilibrium, ensuring host survival and viral persistence. Thus, our study sheds light on archaeal virus-host interactions and highlights similarities with bacteriophages in establishing stable coexistence with their hosts in the gut.

Improving viral annotation with artificial intelligence

Viruses of bacteria, "phages," are fundamental, poorly understood components of microbial community structure and function. Additionally, their dependence on hosts for replication positions phages as unique sensors of ecosystem features and environmental pressures. High-throughput sequencing approaches have begun to give us access to the diversity and range of phage populations in complex microbial community samples, and metagenomics is currently the primary tool with which we study phage populations. The study of phages by metagenomic sequencing, however, is fundamentally limited by viral diversity, which results in the vast majority of viral genomes and metagenome-annotated genomes lacking annotation. To harness bacteriophages for applications in human and environmental health and disease, we need new methods to organize and annotate viral sequence diversity. We recently demonstrated that methods that leverage self-supervised representation learning can supplement statistical sequence representations for remote viral protein homology detection in the ocean virome and propose that consideration of the functional content of viral sequences allows for the identification of similarity in otherwise sequence-diverse viruses and viral-like elements for biological discovery. In this review, we describe the potential and pitfalls of large language models for viral annotation. We describe the need for new approaches to annotate viral sequences in metagenomes, the fundamentals of what protein language models are and how one can use them for sequence annotation, the strengths and weaknesses of these models, and future directions toward developing better models for viral annotation more broadly.

Massive expansion of the pig gut virome based on global metagenomic mining

The pig gut virome plays a vital role in the gut microbial ecosystem of pigs. However, a comprehensive understanding of their diversity and a reference database for the virome are currently lacking. To address this gap, we established a Pig Virome Database (PVD) that comprised of 5,566,804 viral contig sequences from 4650 publicly available gut metagenomic samples using a pipeline designated "metav". By clustering sequences, we identified 48,299 viral operational taxonomic units (vOTUs) genomes of at least medium quality, of which 92.83% of which were not found in existing major databases. The majority of vOTUs were identified as Caudoviricetes (72.21%). The PVD database contained a total of 2,362,631 protein-coding genes across the above medium-quality vOTUs genomes that can be used to explore the functional potential of the pig gut virome. These findings highlight the extensive diversity of viruses in the pig gut and provide a pivotal reference dataset for forthcoming research concerning the pig gut virome.

Metagenomic discovery of microbial eukaryotes in stool microbiomes

Host-associated microbiota form complex microbial communities that are increasingly associated with host behavior and disease. While these microbes include bacterial, archaeal, viral, and eukaryotic constituents, most studies have focused on bacteria due to their dominance in the human host and available tools for investigation. Accumulating evidence suggests microbial eukaryotes in the microbiome play pivotal roles in host health, but our understandings of these interactions are limited to a few readily identifiable taxa because of technical limitations in unbiased eukaryote exploration. Here, we combined cell sorting, optimized eukaryotic cell lysis, and shotgun sequencing to accelerate metagenomic discovery and analysis of host-associated microbial eukaryotes. Using synthetic communities with a 1% microbial eukaryote representation, the eukaryote-optimized cell lysis and DNA recovery method alone yielded a 38-fold increase in eukaryotic DNA. Automated sorting of eukaryotic cells from stool samples of healthy adults increased the number of microbial eukaryote reads in metagenomic pools by up to 28-fold compared to commercial kits. Read frequencies for identified fungi increased by 10,000× on average compared to the Human Microbiome Project and allowed for the identification of novel taxa, de novo assembly of contigs from previously unknown microbial eukaryotes, and gene prediction from recovered genomic segments. These advances pave the way for the unbiased inclusion of microbial eukaryotes in deciphering determinants of health and disease in the host-associated microbiome.IMPORTANCEMicrobial eukaryotes are common constituents of the human gut where they can contribute to local ecology and host health, but they are often overlooked in microbiome studies. The lack of attention is due to current technical limitations that are heavily biased or poorly recovered DNA from microbial eukaryotes. We developed a method to increase the representation of these eukaryotes in metagenomic sequencing of microbiome samples that allows to improve their detection compared to prior methods and allows for the identification of new species. Application of the technique to gut microbiome samples improved detection of fungi, protists, and helminths. New eukaryotic taxa and their encoded genes could be identified by sequencing a small number of samples. This approach can improve the inclusion of eukaryotes into microbiome research.

Lactiplantibacillus plantarum P9 for chronic diarrhea in young adults: a large double-blind, randomized, placebo-controlled trial

Current treatments for chronic diarrhea have limited efficacy and several side effects. Probiotics have the potential to alleviate symptoms of diarrhea. This randomized, double-blind, placebo-controlled trial evaluates the effects of administering the probiotic Lactiplantibacillus plantarum P9 (P9) strain in young adults with chronic diarrhea (Clinical Trial Registration Number: ChiCTR2000038410). The intervention period lasts for 28 days, followed by a 14-day post-intervention period. Participants are randomized into the P9 (n = 93) and placebo (n = 96) groups, with 170 individuals completing the double-blind intervention phase (n = 85 per group). The primary endpoint is the diarrhea symptom severity score. Both intention-to-treat (n = 189) and per-protocol (n = 170) analyses reveal a modest yet statistically significant reduction in diarrhea severity compared to the placebo group (20.0%, P = 0.050; 21.4%, P = 0.048, respectively). In conclusion, the results of this study support the use of probiotics in managing chronic diarrhea in young adults. However, the lack of blood parameter assessment and the short intervention period represent limitations of this study.

Dietary Effects on the Gut Phageome

As knowledge of the gut microbiome has expanded our understanding of the symbiotic and dysbiotic relationships between the human host and its microbial constituents, the influence of gastrointestinal (GI) microbes both locally and beyond the intestine has become evident. Shifts in bacterial populations have now been associated with several conditions including Crohn's disease (CD), Ulcerative Colitis (UC), irritable bowel syndrome (IBS), Alzheimer's disease, Parkinson's Disease, liver diseases, obesity, metabolic syndrome, anxiety, depression, and cancers. As the bacteria in our gut thrive on the food we eat, diet plays a critical role in the functional aspects of our gut microbiome, influencing not only health but also the development of disease. While the bacterial microbiome in the context of disease is well studied, the associated gut phageome-bacteriophages living amongst and within our bacterial microbiome-is less well understood. With growing evidence that fluctuations in the phageome also correlate with dysbiosis, how diet influences this population needs to be better understood. This review surveys the current understanding of the effects of diet on the gut phageome.

Unveiling the unknown viral world in groundwater

Viruses as the prevailing biological entities are poorly understood in underground realms. Here, we establish the first metagenomic Groundwater Virome Catalogue (GWVC) comprising 280,420 viral species ( ≥ 5 kb) detected from 607 monitored wells in seven geo-environmental zones throughout China. In expanding ~10-fold the global portfolio of known groundwater viruses, we uncover over 99% novel viruses and about 95% novel viral clusters. By linking viruses to hosts from 119 prokaryotic phyla, we double the number of microbial phyla known to be virus-infected in groundwater. As keystone ultrasmall symbionts in aquifers, CPR bacteria and DPANN archaea are susceptible to virulent viruses. Certain complete CPR viruses even likely infect non-CPR bacteria, while partial CPR/DPANN viruses harbor cell-surface modification genes that assist symbiont cell adhesion to free-living microbes. This study reveals the unknown viral world and auxiliary metabolism associated with methane, nitrogen, sulfur, and phosphorus cycling in groundwater, and highlights the importance of subsurface virosphere in viral ecology.

A multi-kingdom collection of 33,804 reference genomes for the human vaginal microbiome

The human vagina harbours diverse microorganisms-bacteria, viruses and fungi-with profound implications for women's health. Genome-level analysis of the vaginal microbiome across multiple kingdoms remains limited. Here we utilize metagenomic sequencing data and fungal cultivation to establish the Vaginal Microbial Genome Collection (VMGC), comprising 33,804 microbial genomes spanning 786 prokaryotic species, 11 fungal species and 4,263 viral operational taxonomic units. Notably, over 25% of prokaryotic species and 85% of viral operational taxonomic units remain uncultured. This collection significantly enriches genomic diversity, especially for prevalent vaginal pathogens such as BVAB1 (an uncultured bacterial vaginosis-associated bacterium) and Amygdalobacter spp. (BVAB2 and related species). Leveraging VMGC, we characterize functional traits of prokaryotes, notably Saccharofermentanales (an underexplored yet prevalent order), along with prokaryotic and eukaryotic viruses, offering insights into their niche adaptation and potential roles in the vagina. VMGC serves as a valuable resource for studying vaginal microbiota and its impact on vaginal health.

Next-Generation Sequencing for Characterizing Respiratory Tract Virome and Improving Detection of Viral Pathogens in Children With Pneumonia

BACKGROUND

Pneumonia is typically caused by a variety of pathogenic microorganisms. Traditional research often focuses on the infection of a few microorganisms, whereas metagenomic studies focus on the impact of the bacteriome and mycobiome on respiratory diseases. Reports on the virome characteristics of pediatric pneumonia remain relatively scarce.

METHODS

We employed de novo assembly and combined homology- and feature-based methods to characterize the respiratory virome in whole-genome DNA sequencing samples from oropharynx (OP) swabs, nasopharynx (NP) swabs, and bronchoalveolar lavage fluids (BALF) of children with pneumonia.

RESULTS

Significant differences were observed in the alpha and beta diversity indexes, as well as in the composition of the oropharyngeal virome, between pneumonia cases and controls. We identified 1137 viral operational taxonomic units (vOTUs) with significant differences, indicating a preference of pneumonia-reduced vOTUs for infecting Prevotella, Neisseria, and Veillonella, whereas pneumonia-enriched vOTUs included polyomavirus, human adenovirus, and phages targeting Staphylococcus, Streptococcus, Granulicatella, and Actinomyces. Comparative analysis revealed higher relative abundances and prevalence rates of pneumonia-enriched OP vOTUs in NP and BALF samples compared to pneumonia-reduced vOTUs. Additionally, virome analysis identified six pediatric patients with severe human adenovirus or polyomavirus infections, five of whom might have been undetected by targeted polymerase chain reaction (PCR)-based testing.

CONCLUSIONS

This study offers insights into pediatric pneumonia respiratory viromes, highlighting frequent transmission of potentially pathogenic viruses and demonstrating virome analysis as a valuable adjunct for pathogen detection.

Protein Set Transformer: A protein-based genome language model to power high diversity viromics

Exponential increases in microbial and viral genomic data demand transformational advances in scalable, generalizable frameworks for their interpretation. Standard homology-based functional analyses are hindered by the rapid divergence of microbial and especially viral genomes and proteins that significantly decreases the volume of usable data. Here, we present Protein Set Transformer (PST), a protein-based genome language model that models genomes as sets of proteins without considering sparsely available functional labels. Trained on >100k viruses, PST outperformed other homology- and language model-based approaches for relating viral genomes based on shared protein content. Further, PST demonstrated protein structural and functional awareness by clustering capsid-fold-containing proteins with known capsid proteins and uniquely clustering late gene proteins within related viruses. Our data establish PST as a valuable method for diverse viral genomics, ecology, and evolutionary applications. We posit that the PST framework can be a foundation model for microbial genomics when trained on suitable data.

Phylogeny and disease associations of a widespread and ancient intestinal bacteriophage lineage

Viruses are core components of the human microbiome, impacting health through interactions with gut bacteria and the immune system. Most human microbiome viruses are bacteriophages, which exclusively infect bacteria. Until recently, most gut virome studies focused on low taxonomic resolution (e.g., viral operational taxonomic units), hampering population-level analyses. We previously identified an expansive and widespread bacteriophage lineage in inhabitants of Amsterdam, the Netherlands. Here, we study their biodiversity and evolution in various human populations. Based on a phylogeny using sequences from six viral genome databases, we propose the Candidatus order Heliusvirales. We identify heliusviruses in 82% of 5441 individuals across 39 studies, and in nine metagenomes from humans that lived in Europe and North America between 1000 and 5000 years ago. We show that a large lineage started to diversify when Homo sapiens first appeared some 300,000 years ago. Ancient peoples and modern hunter-gatherers have distinct Ca. Heliusvirales populations with lower richness than modern urbanized people. Urbanized people suffering from type 1 and type 2 diabetes, as well as inflammatory bowel disease, have higher Ca. Heliusvirales richness than healthy controls. We thus conclude that these ancient core members of the human gut virome have thrived with increasingly westernized lifestyles.

Large-scale investigation of species-specific orphan genes in the human gut microbiome elucidates their evolutionary origins

Species-specific genes, also known as orphans, are ubiquitous across life's domains. In prokaryotes, species-specific orphan genes (SSOGs) are mostly thought to originate in external elements such as viruses followed by horizontal gene transfer, whereas the scenario of native origination, through rapid divergence or de novo, is mostly dismissed. However, quantitative evidence supporting either scenario is lacking. Here, we systematically analyzed genomes from 4644 human gut microbiome species and identified more than 600,000 unique SSOGs, representing an average of 2.6% of a given species' pangenome. These sequences are mostly rare within each species yet show signs of purifying selection. Overall, SSOGs use optimal codons less frequently, and their proteins are more disordered than those of conserved genes (i.e., non-SSOGs). Importantly, across species, the GC content of SSOGs closely matches that of conserved ones. In contrast, the ∼5% of SSOGs that share similarity to known viral sequences have distinct characteristics, including lower GC content. Thus, SSOGs with similarity to viruses differ from the remaining SSOGs, contrasting an external origination scenario for most of them. By examining the orthologous genomic region in closely related species, we show that a small subset of SSOGs likely evolved natively de novo and find that these genes also differ in their properties from the remaining SSOGs. Our results challenge the notion that external elements are the dominant source of prokaryotic genetic novelty and will enable future studies into the biological role and relevance of species-specific genes in the human gut.

Gut virome-wide association analysis identifies cross-population viral signatures for inflammatory bowel disease

BACKGROUND

The gut virome has been implicated in inflammatory bowel disease (IBD), yet a full understanding of the gut virome in IBD patients, especially across diverse geographic populations, is lacking.

RESULTS

In this study, we conducted a comprehensive gut virome-wide association study in a Chinese cohort of 71 IBD patients (15 with Crohn's disease and 56 with ulcerative colitis) and 77 healthy controls via viral-like particle (VLP) and bulk virome sequencing of their feces. By utilizing an integrated gut virus catalog tailored to the IBD virome, we revealed fundamental alterations in the gut virome in IBD patients. These characterized 139 differentially abundant viral signatures, including elevated phages predicted to infect Escherichia, Klebsiella, Enterococcus_B, Streptococcus, and Veillonella species, as well as IBD-depleted phages targeting Prevotella, Ruminococcus_E, Bifidobacterium, and Blautia species. Remarkably, these viral signatures demonstrated high consistency across diverse populations such as those in Europe and the USA, emphasizing their significance and broad relevance in the disease context. Furthermore, fecal virome transplantation experiments verified that the colonization of these IBD-characterized viruses can modulate experimental colitis in mouse models.

CONCLUSIONS

Building upon these insights into the IBD gut virome, we identified potential biomarkers for prognosis and therapy in IBD patients, laying the foundation for further exploration of viromes in related conditions. Video Abstract.

VirRep: a hybrid language representation learning framework for identifying viruses from human gut metagenomes

Identifying viruses from metagenomes is a common step to explore the virus composition in the human gut. Here, we introduce VirRep, a hybrid language representation learning framework, for identifying viruses from human gut metagenomes. VirRep combines a context-aware encoder and an evolution-aware encoder to improve sequence representation by incorporating k-mer patterns and sequence homologies. Benchmarking on both simulated and real datasets with varying viral proportions demonstrates that VirRep outperforms state-of-the-art methods. When applied to fecal metagenomes from a colorectal cancer cohort, VirRep identifies 39 high-quality viral species associated with the disease, many of which cannot be detected by existing methods.

Integrating functional metagenomics to decipher microbiome-immune interactions

Microbial metabolites can be viewed as the cytokines of the microbiome, transmitting information about the microbial and metabolic environment of the gut to orchestrate and modulate local and systemic immune responses. Still, many immunology studies focus solely on the taxonomy and community structure of the gut microbiota rather than its functions. Early sequencing-based microbiota profiling approaches relied on PCR amplification of small regions of bacterial and fungal genomes to facilitate identification of the microbes present. However, recent microbiome analysis methods, particularly shotgun metagenomic sequencing, now enable culture-independent profiling of microbiome functions and metabolites in addition to taxonomic characterization. In this review, we showcase recent advances in functional metagenomics methods and applications and discuss the current limitations and potential avenues for future development. Importantly, we highlight a few examples of key areas of opportunity in immunology research where integrating functional metagenomic analyses of the microbiome can substantially enhance a mechanistic understanding of microbiome-immune interactions and their contributions to health and disease states.

Virus ecology and 7-year temporal dynamics across a permafrost thaw gradient

Soil microorganisms are pivotal in the global carbon cycle, but the viruses that affect them and their impact on ecosystems are less understood. In this study, we explored the diversity, dynamics, and ecology of soil viruses through 379 metagenomes collected annually from 2010 to 2017. These samples spanned the seasonally thawed active layer of a permafrost thaw gradient, which included palsa, bog, and fen habitats. We identified 5051 virus operational taxonomic units (vOTUs), doubling the known viruses for this site. These vOTUs were largely ephemeral within habitats, suggesting a turnover at the vOTU level from year to year. While the diversity varied by thaw stage and depth-related patterns were specific to each habitat, the virus communities did not significantly change over time. The abundance ratios of virus to host at the phylum level did not show consistent trends across the thaw gradient, depth, or time. To assess potential ecosystem impacts, we predicted hosts in silico and found viruses linked to microbial lineages involved in the carbon cycle, such as methanotrophy and methanogenesis. This included the identification of viruses of Candidatus Methanoflorens, a significant global methane contributor. We also detected a variety of potential auxiliary metabolic genes, including 24 carbon-degrading glycoside hydrolases, six of which are uniquely terrestrial. In conclusion, these long-term observations enhance our understanding of soil viruses in the context of climate-relevant processes and provide opportunities to explore their role in terrestrial carbon cycling.

A global atlas of soil viruses reveals unexplored biodiversity and potential biogeochemical impacts

Historically neglected by microbial ecologists, soil viruses are now thought to be critical to global biogeochemical cycles. However, our understanding of their global distribution, activities and interactions with the soil microbiome remains limited. Here we present the Global Soil Virus Atlas, a comprehensive dataset compiled from 2,953 previously sequenced soil metagenomes and composed of 616,935 uncultivated viral genomes and 38,508 unique viral operational taxonomic units. Rarefaction curves from the Global Soil Virus Atlas indicate that most soil viral diversity remains unexplored, further underscored by high spatial turnover and low rates of shared viral operational taxonomic units across samples. By examining genes associated with biogeochemical functions, we also demonstrate the viral potential to impact soil carbon and nutrient cycling. This study represents an extensive characterization of soil viral diversity and provides a foundation for developing testable hypotheses regarding the role of the virosphere in the soil microbiome and global biogeochemistry.

Longitudinal multi-omics analysis of host microbiome architecture and immune responses during short-term spaceflight

Maintenance of astronaut health during spaceflight will require monitoring and potentially modulating their microbiomes. However, documenting microbial shifts during spaceflight has been difficult due to mission constraints that lead to limited sampling and profiling. Here we executed a six-month longitudinal study to quantify the high-resolution human microbiome response to three days in orbit for four individuals. Using paired metagenomics and metatranscriptomics alongside single-nuclei immune cell profiling, we characterized time-dependent, multikingdom microbiome changes across 750 samples and 10 body sites before, during and after spaceflight at eight timepoints. We found that most alterations were transient across body sites; for example, viruses increased in skin sites mostly during flight. However, longer-term shifts were observed in the oral microbiome, including increased plaque-associated bacteria (for example, Fusobacteriota), which correlated with immune cell gene expression. Further, microbial genes associated with phage activity, toxin-antitoxin systems and stress response were enriched across multiple body sites. In total, this study reveals in-depth characterization of microbiome and immune response shifts experienced by astronauts during short-term spaceflight and the associated changes to the living environment, which can help guide future missions, spacecraft design and space habitat planning.

Distinct characteristics of the gut virome in patients with osteoarthritis and gouty arthritis

BACKGROUND/PURPOSE(S)

The gut microbiota and its metabolites play crucial roles in pathogenesis of arthritis, highlighting gut microbiota as a promising avenue for modulating autoimmunity. However, the characterization of the gut virome in arthritis patients, including osteoarthritis (OA) and gouty arthritis (GA), requires further investigation.

METHODS

We employed virus-like particle (VLP)-based metagenomic sequencing to analyze gut viral community in 20 OA patients, 26 GA patients, and 31 healthy controls, encompassing a total of 77 fecal samples.

RESULTS

Our analysis generated 6819 vOTUs, with a considerable proportion of viral genomes differing from existing catalogs. The gut virome in OA and GA patients differed significantly from healthy controls, showing variations in diversity and viral family abundances. We identified 157 OA-associated and 94 GA-associated vOTUs, achieving high accuracy in patient-control discrimination with random forest models. OA-associated viruses were predicted to infect pro-inflammatory bacteria or bacteria associated with immunoglobulin A production, while GA-associated viruses were linked to Bacteroidaceae or Lachnospiraceae phages. Furthermore, several viral functional orthologs displayed significant differences in frequency between OA-enriched and GA-enriched vOTUs, suggesting potential functional roles of these viruses. Additionally, we trained classification models based on gut viral signatures to effectively discriminate OA or GA patients from healthy controls, yielding AUC values up to 0.97, indicating the clinical utility of the gut virome in diagnosing OA or GA.

CONCLUSION

Our study highlights distinctive alterations in viral diversity and taxonomy within gut virome of OA and GA patients, offering insights into arthritis etiology and potential treatment and prevention strategies.

Structural Requirements for Reverse Transcription by a Diversity-generating Retroelement

Diversity-generating retroelements (DGRs) create massive protein sequence variation in ecologically diverse microbes. Variation occurs during reverse transcription of a protein-encoding RNA template coupled to misincorporation at adenosines. In the prototypical Bordetella bacteriophage DGR, the template must be surrounded by upstream and downstream RNA segments for cDNA synthesis by the reverse transcriptase bRT and associated protein Avd. The function of the surrounding RNA was unknown. Cryo-EM revealed that this RNA enveloped bRT and lay over barrel-shaped Avd, forming an intimate ribonucleoprotein (RNP). An abundance of essential interactions between RNA structural elements and bRT-Avd precisely positioned an RNA homoduplex for initiation of cDNA synthesis by cis -priming. Our results explain how the surrounding RNA primes cDNA synthesis, promotes processivity, terminates polymerization, and strictly limits mutagenesis to select proteins through mechanisms that are likely conserved in DGRs from distant taxa.

Characterizing the gut phageome and phage-borne antimicrobial resistance genes in pigs

BACKGROUND

Mammalian intestine harbors a mass of phages that play important roles in maintaining gut microbial ecosystem and host health. Pig has become a common model for biomedical research and provides a large amount of meat for human consumption. However, the knowledge of gut phages in pigs is still limited.

RESULTS

Here, we investigated the gut phageome in 112 pigs from seven pig breeds using PhaBOX strategy based on the metagenomic data. A total of 174,897 non-redundant gut phage genomes were assembled from 112 metagenomes. A total of 33,487 gut phage genomes were classified and these phages mainly belonged to phage families such as Ackermannviridae, Straboviridae, Peduoviridae, Zierdtviridae, Drexlerviridae, and Herelleviridae. The gut phages in seven pig breeds exhibited distinct communities and the gut phage communities changed with the age of pig. These gut phages were predicted to infect a broad range of 212 genera of prokaryotes, such as Candidatus Hamiltonella, Mycoplasma, Colwellia, and Lactobacillus. The data indicated that broad KEGG and CAZy functions were also enriched in gut phages of pigs. The gut phages also carried the antimicrobial resistance genes (ARGs) and the most abundant antimicrobial resistance genotype was diaminopyrimidine resistance.

CONCLUSIONS

Our research delineates a landscape for gut phages in seven pig breeds and reveals that gut phages serve as a key reservoir of ARGs in pigs. Video Abstract.

A panoramic view of the virosphere in three wastewater treatment plants by integrating viral-like particle-concentrated and traditional non-concentrated metagenomic approaches

Wastewater biotreatment systems harbor a rich diversity of microorganisms, and the effectiveness of biotreatment systems largely depends on the activity of these microorganisms. Specifically, viruses play a crucial role in altering microbial behavior and metabolic processes throughout their infection phases, an aspect that has recently attracted considerable interest. Two metagenomic approaches, viral-like particle-concentrated (VPC, representing free viral-like particles) and non-concentrated (NC, representing the cellular fraction), were employed to assess their efficacy in revealing virome characteristics, including taxonomy, diversity, host interactions, lifestyle, dynamics, and functional genes across processing units of three wastewater treatment plants (WWTPs). Our findings indicate that each approach offers unique insights into the viral community and functional composition. Their combined use proved effective in elucidating WWTP viromes. We identified nearly 50,000 viral contigs, with Cressdnaviricota and Uroviricota being the predominant phyla in the VPC and NC fractions, respectively. Notably, two pathogenic viral families, Asfarviridae and Adenoviridae, were commonly found in these WWTPs. We also observed significant differences in the viromes of WWTPs processing different types of wastewater. Additionally, various phage-derived auxiliary metabolic genes (AMGs) were active at the RNA level, contributing to the metabolism of the microbial community, particularly in carbon, sulfur, and phosphorus cycling. Moreover, we identified 29 virus-carried antibiotic resistance genes (ARGs) with potential for host transfer, highlighting the role of viruses in spreading ARGs in the environment. Overall, this study provides a detailed and integrated view of the virosphere in three WWTPs through the application of VPC and NC metagenomic approaches. Our findings enhance the understanding of viral communities, offering valuable insights for optimizing the operation and regulation of wastewater treatment systems.

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.

Characterizations of gut bacteriome, mycobiome, and virome of healthy individuals living in sea-level and high-altitude areas

BACKGROUND

The contribution of gut microbiota to human high-altitude adaptation remains inadequately understood.

METHODS

Here a comparative analysis of gut microbiota was conducted between healthy individuals living at sea level and high altitude using deep whole-metagenome shotgun sequencing, to investigate the adaptive mechanisms of gut microbiota in plateau inhabitants.

RESULTS

The results showed the gut bacteriomes in high-altitude individuals exhibited greater within-sample diversity and significant alterations in both bacterial compositional and functional profiles when compared to those of sea-level individuals, indicating the potential selection of unique bacteria associated with high-altitude environments. The strain-level investigation revealed enrichment of Collinsella aerofaciens and Akkermansia muciniphila in high-altitude populations. The characteristics of gut virome and gut mycobiome were also investigated. Compared to sea-level subjects, high-altitude subjects exhibited a greater diversity in their gut virome, with an increased number of viral operational taxonomic units (vOTUs) and unique annotated genes. Finally, correlation analyses revealed 819 significant correlations between 42 bacterial species and 375 vOTUs, while no significant correlations were observed between bacteria and fungi or between fungi and viruses.

CONCLUSION

The findings have significantly contributed to an enhanced comprehension of the mechanisms underlying the high-altitude geographic adaptation of the human gut microbiota.

Assessing phage-host population dynamics by reintroducing virulent viruses to synthetic microbiomes

Microbiomes feature complex interactions between diverse bacteria and bacteriophages. Synthetic microbiomes offer a powerful way to study these interactions; however, a major challenge is obtaining a representative bacteriophage population during the bacterial isolation process. We demonstrate that colony isolation reliably excludes virulent viruses from sample sources with low virion-to-bacteria ratios such as feces, creating "virulent virus-free" controls. When the virulent dsDNA virome is reintroduced to a 73-strain synthetic gut microbiome in a bioreactor model of the human colon, virulent viruses target susceptible strains without significantly altering community structure or metabolism. In addition, we detected signals of prophage induction that associate with virulent predation. Overall, our findings indicate that dilution-based isolation methods generate synthetic gut microbiomes that are heavily depleted, if not devoid, of virulent viruses and that such viruses, if reintroduced, have a targeted effect on community assembly, metabolism, and prophage replication.

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.

Petabase-Scale 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.

Genome engineering of the human gut microbiome

The human gut microbiome, a complex ecosystem, significantly influences host health, impacting crucial aspects such as metabolism and immunity. To enhance our comprehension and control of the molecular mechanisms orchestrating the intricate interplay between gut commensal bacteria and human health, the exploration of genome engineering for gut microbes is a promising frontier. Nevertheless, the complexities and diversities inherent in the gut microbiome pose substantial challenges to the development of effective genome engineering tools for human gut microbes. In this comprehensive review, we provide an overview of the current progress and challenges in genome engineering of human gut commensal bacteria, whether executed in vitro or in situ. A specific focus is directed towards the advancements and prospects in cargo DNA delivery and high-throughput techniques. Additionally, we elucidate the immense potential of genome engineering methods to enhance our understanding of the human gut microbiome and engineer the microorganisms to enhance human health.

Improved sampling and DNA extraction procedures for microbiome analysis in food-processing environments

Deep investigation of the microbiome of food-production and food-processing environments through whole-metagenome sequencing (WMS) can provide detailed information on the taxonomic composition and functional potential of the microbial communities that inhabit them, with huge potential benefits for environmental monitoring programs. However, certain technical challenges jeopardize the application of WMS technologies with this aim, with the most relevant one being the recovery of a sufficient amount of DNA from the frequently low-biomass samples collected from the equipment, tools and surfaces of food-processing plants. Here, we present the first complete workflow, with optimized DNA-purification methodology, to obtain high-quality WMS sequencing results from samples taken from food-production and food-processing environments and reconstruct metagenome assembled genomes (MAGs). The protocol can yield DNA loads >10 ng in >98% of samples and >500 ng in 57.1% of samples and allows the collection of, on average, 12.2 MAGs per sample (with up to 62 MAGs in a single sample) in ~1 week, including both laboratory and computational work. This markedly improves on results previously obtained in studies performing WMS of processing environments and using other protocols not specifically developed to sequence these types of sample, in which <2 MAGs per sample were obtained. The full protocol has been developed and applied in the framework of the European Union project MASTER (Microbiome applications for sustainable food systems through technologies and enterprise) in 114 food-processing facilities from different production sectors.

Ten computational challenges in human virome studies

In recent years, substantial advancements have been achieved in understanding the diversity of the human virome and its intricate roles in human health and diseases. Despite this progress, the field of human virome research remains nascent, primarily hindered by the absence of effective methods, particularly in the domain of computational tools. This perspective systematically outlines ten computational challenges spanning various types of virome studies. These challenges arise due to the vast diversity of viromes, the absence of a universal marker gene in viral genomes, the low abundance of virus populations, the remote or minimal homology of viral proteins to known proteins, and the highly dynamic and heterogeneous nature of viromes. For each computational challenge, we discuss the underlying reasons, current research progress, and potential solutions. The resolution of these challenges necessitates ongoing collaboration among computational scientists, virologists, and multidisciplinary experts. In essence, this perspective serves as a comprehensive guide for directing computational efforts in human virome studies.

Molecular characterization of plasma virome of hepatocellular carcinoma (HCC) patients

Hepatocellular carcinoma (HCC) stands as the most common cancer type, arising from various causes, and responsible for a substantial number of cancer-related fatalities. Recent advancements in viral metagenomics have empowered scientists to delve into the intricate diversity of the virosphere, viral evolution, interactions between viruses and their hosts, and the identification of viral causes behind disease outbreaks, the development of specific symptoms, and their potential role in altering the host's physiology. The present study had the objective of "Molecular Characterization of HBV, HCV, anelloviruses, CMV, SENV-D, SENV-H, HEV, and HPV viruses among individuals suffering from HCC." A total of 381 HCC patients contributed 10 cc of blood each for this study. The research encompassed the assessment of tumor markers, followed by molecular characterization of HBV, HCV, Anelloviruses (TTV, TTMV, and TTMDV), SENV-H and SENV-D viruses, HEV, CMV, and HPV, as well as histopathological examinations. The outcomes of this study revealed that majority of the HCC patients 72.4% (276/381) were male as compared to females. HCV infection, at 76.4% (291 out of 381), exhibited a significant association (p < 0.05) with HCC. Most patients displayed singular lesions in the liver, with Child Pugh Score Type B being the predominant finding in 45.2% of cases. Plasma virome analysis indicated the prevalence of TTMDV (75%), followed by TTMV (70%) and TTV (42.1%) among anelloviruses in HCC patients. Similarly, SENV-H (52%) was followed by SENV-D (20%), with co-infections at 15%. The presence of CMV and HEV among the HCC patients was recorded 5% each however 3.5% of the patients showed the presence of HPV. In conclusion, this study underscores that HCC patients serve as reservoirs for various pathogenic and non-pathogenic viruses, potentially contributing to the development, progression, and severity of the disease.

Deep metagenomic characterization of the gut virome in pregnant women with preeclampsia

Preeclampsia (PE), a pregnancy-specific syndrome, has been associated with the gut bacteriome. Here, to investigate the impact of the gut virome on the development of PE, we identified over 8,000 nonredundant viruses from the fecal metagenomes of 40 early-onset PE and 37 healthy pregnant women and profiled their abundances. Comparison and correlation analysis showed that PE-enriched viruses frequently connected to Blautia species enriched in PE. By contrast, bacteria linked to PE-depleted viruses were often the Bacteroidaceae members such as Bacteroides spp., Phocaeicola spp., Parabacteroides spp., and Alistipes shahii. In terms of viral function, PE-depleted viruses had auxiliary metabolic genes that participated in the metabolism of simple and complex polysaccharides, sulfur metabolism, lipopolysaccharide biosynthesis, and peptidoglycan biosynthesis, while PE-enriched viruses had a gene encoding cyclic pyranopterin monophosphate synthase, which seemed to be special, that participates in the biosynthesis of the molybdenum cofactor. Furthermore, the classification model based on gut viral signatures was developed to discriminate PE patients from healthy controls and showed an area under the receiver operating characteristic curve of 0.922 that was better than that of the bacterium-based model. This study opens up new avenues for further research, providing valuable insights into the PE gut virome and offering potential directions for future mechanistic and therapeutic investigations, with the ultimate goal of improving the diagnosis and management of PE.IMPORTANCEThe importance of this study lies in its exploration of the previously overlooked but potentially critical role of the gut virome in preeclampsia (PE). While the association between PE and the gut bacteriome has been recognized, this research takes a pioneering step into understanding how the gut virome, represented by over 8,000 nonredundant viruses, contributes to this condition. The findings reveal intriguing connections between PE-enriched viruses and specific gut bacteria, such as the prevalence of Blautia species in individuals with PE, contrasting with bacteria linked to PE-depleted viruses, including members of the Bacteroidaceae family. These viral interactions and associations provide a deeper understanding of the complex dynamics at play in PE.

Carbon starvation raises capacities in bacterial antibiotic resistance and viral auxiliary carbon metabolism in soils

Organic carbon availability in soil is crucial for shaping microbial communities, yet, uncertainties persist concerning microbial adaptations to carbon levels and the ensuing ecological and evolutionary consequences. We investigated organic carbon metabolism, antibiotic resistance, and virus-host interactions in soils subjected to 40 y of chemical and organic fertilization that led to contrasting carbon availability: carbon-poor and carbon-rich soils, respectively. Carbon-poor soils drove the enrichment of putative genes involved in organic matter decomposition and exhibited specialization in utilizing complex organic compounds, reflecting scramble competition. This specialization confers a competitive advantage of microbial communities in carbon-poor soils but reduces their buffering capacity in terms of organic carbon metabolisms, making them more vulnerable to environmental fluctuations. Additionally, in carbon-poor soils, viral auxiliary metabolic genes linked to organic carbon metabolism increased host competitiveness and environmental adaptability through a strategy akin to "piggyback the winner." Furthermore, putative antibiotic resistance genes, particularly in low-abundance drug categories, were enriched in carbon-poor soils as an evolutionary consequence of chemical warfare (i.e., interference competition). This raises concerns about the potential dissemination of antibiotic resistance from conventional agriculture that relies on chemical-only fertilization. Consequently, carbon starvation resulting from long-term chemical-only fertilization increases microbial adaptations to competition, underscoring the importance of implementing sustainable agricultural practices to mitigate the emergence and spread of antimicrobial resistance and to increase soil carbon storage.

A compendium of ruminant gastrointestinal phage genomes revealed a higher proportion of lytic phages than in any other environments

BACKGROUND

Ruminants are important livestock animals that have a unique digestive system comprising multiple stomach compartments. Despite significant progress in the study of microbiome in the gastrointestinal tract (GIT) sites of ruminants, we still lack an understanding of the viral community of ruminants. Here, we surveyed its viral ecology using 2333 samples from 10 sites along the GIT of 8 ruminant species.

RESULTS

We present the Unified Ruminant Phage Catalogue (URPC), a comprehensive survey of phages in the GITs of ruminants including 64,922 non-redundant phage genomes. We characterized the distributions of the phage genomes in different ruminants and GIT sites and found that most phages were organism-specific. We revealed that ~ 60% of the ruminant phages were lytic, which was the highest as compared with those in all other environments and certainly will facilitate their applications in microbial interventions. To further facilitate the future applications of the phages, we also constructed a comprehensive virus-bacteria/archaea interaction network and identified dozens of phages that may have lytic effects on methanogenic archaea.

CONCLUSIONS

The URPC dataset represents a useful resource for future microbial interventions to improve ruminant production and ecological environmental qualities. Phages have great potential for controlling pathogenic bacterial/archaeal species and reducing methane emissions. Our findings provide insights into the virome ecology research of the ruminant GIT and offer a starting point for future research on phage therapy in ruminants. Video Abstract.

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

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

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.

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.

Biogeographic patterns and drivers of soil viromes

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

Correlation between the gut microbiome and neurodegenerative diseases: a review of metagenomics evidence

A growing body of evidence suggests that the gut microbiota contributes to the development of neurodegenerative diseases via the microbiota-gut-brain axis. As a contributing factor, microbiota dysbiosis always occurs in pathological changes of neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. High-throughput sequencing technology has helped to reveal that the bidirectional communication between the central nervous system and the enteric nervous system is facilitated by the microbiota's diverse microorganisms, and for both neuroimmune and neuroendocrine systems. Here, we summarize the bioinformatics analysis and wet-biology validation for the gut metagenomics in neurodegenerative diseases, with an emphasis on multi-omics studies and the gut virome. The pathogen-associated signaling biomarkers for identifying brain disorders and potential therapeutic targets are also elucidated. Finally, we discuss the role of diet, prebiotics, probiotics, postbiotics and exercise interventions in remodeling the microbiome and reducing the symptoms of neurodegenerative diseases.

VirusPredictor: XGBoost-based software to predict virus-related sequences in human data

MOTIVATION

Discovering disease causative pathogens, particularly viruses without reference genomes, poses a technical challenge as they are often unidentifiable through sequence alignment. Machine learning prediction of patient high-throughput sequences unmappable to human and pathogen genomes may reveal sequences originating from uncharacterized viruses. Currently, there is a lack of software specifically designed for accurately predicting such viral sequences in human data.

RESULTS

We developed a fast XGBoost method and software VirusPredictor leveraging an in-house viral genome database. Our two-step XGBoost models first classify each query sequence into one of three groups: infectious virus, endogenous retrovirus (ERV) or non-ERV human. The prediction accuracies increased as the sequences became longer, i.e. 0.76, 0.93, and 0.98 for 150-350 (Illumina short reads), 850-950 (Sanger sequencing data), and 2000-5000 bp sequences, respectively. Then, sequences predicted to be from infectious viruses are further classified into one of six virus taxonomic subgroups, and the accuracies increased from 0.92 to >0.98 when query sequences increased from 150-350 to >850 bp. The results suggest that Illumina short reads should be de novo assembled into contigs (e.g. ∼1000 bp or longer) before prediction whenever possible. We applied VirusPredictor to multiple real genomic and metagenomic datasets and obtained high accuracies. VirusPredictor, a user-friendly open-source Python software, is useful for predicting the origins of patients' unmappable sequences. This study is the first to classify ERVs in infectious viral sequence prediction. This is also the first study combining virus sub-group predictions.

AVAILABILITY AND IMPLEMENTATION

www.dllab.org/software/VirusPredictor.html.

Analysis and culturing of the prototypic crAssphage reveals a phage-plasmid lifestyle

The prototypic crAssphage (Carjivirus communis) is one of the most abundant, prevalent, and persistent gut bacteriophages, yet it remains uncultured and its lifestyle uncharacterized. For the last decade, crAssphage has escaped plaque-dependent culturing efforts, leading us to investigate alternative lifestyles that might explain its widespread success. Through genomic analyses and culturing, we find that crAssphage uses a phage-plasmid lifestyle to persist extrachromosomally. Plasmid-related genes are more highly expressed than those implicated in phage maintenance. Leveraging this finding, we use a plaque-free culturing approach to measure crAssphage replication in culture with Phocaeicola vulgatus, Phocaeicola dorei, and Bacteroides stercoris, revealing a broad host range. We demonstrate that crAssphage persists with its hosts in culture without causing major cell lysis events or integrating into host chromosomes. The ability to switch between phage and plasmid lifestyles within a wide range of hosts contributes to the prolific nature of crAssphage in the human gut microbiome.

Expanding the human gut microbiome atlas of Africa

Population studies are crucial in understanding the complex interplay between the gut microbiome and geographical, lifestyle, genetic, and environmental factors. However, populations from low- and middle-income countries, which represent ~84% of the world population, have been excluded from large-scale gut microbiome research. Here, we present the AWI-Gen 2 Microbiome Project, a cross-sectional gut microbiome study sampling 1,803 women from Burkina Faso, Ghana, Kenya, and South Africa. By intensively engaging with communities that range from rural and horticultural to urban informal settlements and post-industrial, we capture population diversity that represents a far greater breadth of the world's population. Using shotgun metagenomic sequencing, we find that study site explains substantially more microbial variation than disease status. We identify taxa with strong geographic and lifestyle associations, including loss of Treponema and Cryptobacteroides species and gain of Bifidobacterium species in urban populations. We uncover a wealth of prokaryotic and viral novelty, including 1,005 new bacterial metagenome-assembled genomes, and identify phylogeography signatures in Treponema succinifaciens. Finally, we find a microbiome signature of HIV infection that is defined by several taxa not previously associated with HIV, including Dysosmobacter welbionis and Enterocloster sp. This study represents the largest population-representative survey of gut metagenomes of African individuals to date, and paired with extensive clinical biomarkers, demographic data, and lifestyle information, provides extensive opportunity for microbiome-related discovery and research.

Comparative analysis of the vaginal bacteriome and virome in healthy women living in high-altitude and sea-level areas

The vaginal microbiota plays an important role in the health of the female reproductive tract and is closely associated with various pregnancy outcomes and sexually transmitted diseases. Plenty of internal and external factors have strong influence on the changes in a woman's vaginal microbiome. However, the effect of a high-altitude on female vaginal microbiota has not been described. In this study, we characterized the vaginal bacteriome and virome of 13 and 34 healthy women living in high-altitude and sea-level areas, using whole-metagenome shotgun sequencing of their vaginal mucus samples. The results revealed that the vaginal bacteriomes of high-altitude individuals are featured by a significant increase of species diversity, depletion of Lactobacillus crispatus, and more abundant of some anaerobic bacteria, such as Chlamydia trachomatis, Mageeibacillus indolicus, Dialister micraerophilus, and Sneathia amnii). In addition, the vagina samples of sea-level subjects harbor more Lactobacillus strains, whereas the anaerobic bacteroidetes strains mostly appeared in high-altitude subjects. Identified and assembled 191 virus operational taxonomic units (vOTUs), there were significant differences in the abundance of 107 vOTUs between the two groups. Together, the results of this study raised the understanding of bacteriome and virome in the vagina of women at different elevations, and demonstrated that the vaginal microbiome is related to the high-altitude geographic adaptation.

A metagenomic catalog of the early-life human gut virome

Early-life human gut microbiome is a pivotal driver of gut homeostasis and infant health. However, the viral component (known as "virome") remains mostly unexplored. Here, we establish the Early-Life Gut Virome (ELGV), a catalog of 160,478 non-redundant DNA and RNA viral sequences from 8130 gut virus-like particles (VLPs) enriched or bulk metagenomes in the first three years of life. By clustering, 82,141 viral species are identified, 68.3% of which are absent in existing databases built mainly from adults, and 64 and 8 viral species based on VLPs-enriched and bulk metagenomes, respectively, exhibit potentials as biomarkers to distinguish infants from adults. With the largest longitudinal population of infants profiled by either VLPs-enriched or bulk metagenomic sequencing, we track the inherent instability and temporal development of the early-life human gut virome, and identify differential viruses associated with multiple clinical factors. The mother-infant shared virome and interactions between gut virome and bacteriome early in life are further expanded. Together, the ELGV catalog provides the most comprehensive and complete metagenomic blueprint of the early-life human gut virome, facilitating the discovery of pediatric disease-virome associations in future.

Exploring the gut DNA virome in fecal immunochemical test stool samples reveals associations with lifestyle in a large population-based study

Stool samples for fecal immunochemical tests (FIT) are collected in large numbers worldwide as part of colorectal cancer screening programs. Employing FIT samples from 1034 CRCbiome participants, recruited from a Norwegian colorectal cancer screening study, we identify, annotate and characterize more than 18000 DNA viruses, using shotgun metagenome sequencing. Only six percent of them are assigned to a known taxonomic family, with Microviridae being the most prevalent viral family. Linking individual profiles to comprehensive lifestyle and demographic data shows 17/25 of the variables to be associated with the gut virome. Physical activity, smoking, and dietary fiber consumption exhibit strong and consistent associations with both diversity and relative abundance of individual viruses, as well as with enrichment for auxiliary metabolic genes. We demonstrate the suitability of FIT samples for virome analysis, opening an opportunity for large-scale studies of this enigmatic part of the gut microbiome. The diverse viral populations and their connections to the individual lifestyle uncovered herein paves the way for further exploration of the role of the gut virome in health and disease.

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

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

Insights into gut microbiomes in stem cell transplantation by comprehensive shotgun long-read sequencing

The gut microbiome is a diverse ecosystem, dominated by bacteria; however, fungi, phages/viruses, archaea, and protozoa are also important members of the gut microbiota. Exploration of taxonomic compositions beyond bacteria as well as an understanding of the interaction between the bacteriome with the other members is limited using 16S rDNA sequencing. Here, we developed a pipeline enabling the simultaneous interrogation of the gut microbiome (bacteriome, mycobiome, archaeome, eukaryome, DNA virome) and of antibiotic resistance genes based on optimized long-read shotgun metagenomics protocols and custom bioinformatics. Using our pipeline we investigated the longitudinal composition of the gut microbiome in an exploratory clinical study in patients undergoing allogeneic hematopoietic stem cell transplantation (alloHSCT; n = 31). Pre-transplantation microbiomes exhibited a 3-cluster structure, characterized by Bacteroides spp. /Phocaeicola spp., mixed composition and Enterococcus abundances. We revealed substantial inter-individual and temporal variabilities of microbial domain compositions, human DNA, and antibiotic resistance genes during the course of alloHSCT. Interestingly, viruses and fungi accounted for substantial proportions of microbiome content in individual samples. In the course of HSCT, bacterial strains were stable or newly acquired. Our results demonstrate the disruptive potential of alloHSCTon the gut microbiome and pave the way for future comprehensive microbiome studies based on long-read metagenomics.

Dynamic nature of viral and bacterial communities in human faeces

Bacteriophages are a major component of the gut microbiome and are believed to play a role in establishment and stabilization of microbial communities by influencing taxonomic and functional diversity. We show that the activity of lytic and temperate phages can also significantly affect bacterial community structure in a model of extended colonic retention. Intact fresh human feces were incubated anaerobically at 37°C without homogenization and subjected to metagenomic sequencing. We observed subject-specific blooms and collapses of selected bacteriophage and bacterial populations within some individuals. Most notable were striking collapses of Prevotella populations accompanied by increases in specific bacteriophages. In a number of cases, we even observed a shift from one bacterial "enterotype" to another within 48 h. These results confirm that intact feces represents a highly dynamic ecological system and suggests that colonic retention time could have a profound effect on microbiome composition, including a significant impact by bacteriophages.