Probing the "Dark Matter" of the Human Gut Phageome: Culture Assisted Metagenomics Enables Rapid Discovery and Host-Linking for Novel Bacteriophages

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.

Characterization of the gut bacterial and viral microbiota in latent autoimmune diabetes in adults

Latent autoimmune diabetes in adults (LADA) is a heterogeneous disease characterized by autoantibodies against insulin producing pancreatic beta cells and initial lack of need for insulin treatment. The aim of the present study was to investigate if individuals with LADA have an altered gut microbiota relative to non-diabetic control subjects, individuals with type 1 diabetes (T1D), and individuals with type 2 diabetes (T2D). Bacterial community profiling was performed with primers targeting the variable region 4 of the 16S rRNA gene and sequenced. Amplicon sequence variants (ASVs) were generated with DADA2 and annotated to the SILVA database. The gut virome was sequenced, using a viral particle enrichment and metagenomics approach, assembled, and quantified to describe the composition of the viral community. Comparison of the bacterial alpha- and beta-diversity measures revealed that the gut bacteriome of individuals with LADA resembled that of individuals with T2D. Yet, specific genera were found to differ in abundance in individuals with LADA compared with T1D and T2D, indicating that LADA has unique taxonomical features. The virome composition reflected the stability of the most dominant order Caudovirales and the families Siphoviridae, Podoviridae, and Inoviridae, and the dominant family Microviridae. Further studies are needed to confirm these findings.

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.

The gut virome is associated with stress-induced changes in behaviour and immune responses in mice

The microbiota-gut-brain axis has been shown to play an important role in the stress response, but previous work has focused primarily on the role of the bacteriome. The gut virome constitutes a major portion of the microbiome, with bacteriophages having the potential to remodel bacteriome structure and activity. Here we use a mouse model of chronic social stress, and employ 16S rRNA and whole metagenomic sequencing on faecal pellets to determine how the virome is modulated by and contributes to the effects of stress. We found that chronic stress led to behavioural, immune and bacteriome alterations in mice that were associated with changes in the bacteriophage class Caudoviricetes and unassigned viral taxa. To determine whether these changes were causally related to stress-associated behavioural or physiological outcomes, we conducted a faecal virome transplant from mice before stress and autochthonously transferred it to mice undergoing chronic social stress. The transfer of the faecal virome protected against stress-associated behaviour sequelae and restored stress-induced changes in select circulating immune cell populations, cytokine release, bacteriome alterations and gene expression in the amygdala. These data provide evidence that the virome plays a role in the modulation of the microbiota-gut-brain axis during stress, indicating that these viral populations should be considered when designing future microbiome-directed therapies.

Bacteriophage-Host Interactions and Coevolution

Bacteriophages are the most abundant entity on the planet and play very relevant roles in the diversity and abundance of their bacterial hosts. These interactions are subject to several factors, such as the first encounter of the phage with its host bacterium, in which molecular interactions are fundamental. Along with this, these interactions depend on the environment and other communities present. This chapter focuses on these phage-bacteria interactions, reviewing the knowledge of the early stage (receptor-binding proteins), host responses (resistance and counter-resistance), and ecological and evolutionary models described to date. In general, knowledge has focused on a few phage-bacteria models and has been deepened by sequencing and metagenomics. The study of phage-bacteria interactions is an essential step for the development of therapies and other applications of phages in the clinical and productive environment.

Does the Human Gut Virome Contribute to Host Health or Disease?

The human gastrointestinal (GI) tract harbors eukaryotic and prokaryotic viruses and their genomes, metabolites, and proteins, collectively known as the "gut virome". This complex community of viruses colonizing the enteric mucosa is pivotal in regulating host immunity. The mechanisms involved in cross communication between mucosal immunity and the gut virome, as well as their relationship in health and disease, remain largely unknown. Herein, we review the literature on the human gut virome's composition and evolution and the interplay between the gut virome and enteric mucosal immunity and their molecular mechanisms. Our review suggests that future research efforts should focus on unraveling the mechanisms of gut viruses in human homeostasis and pathophysiology and on developing virus-prompted precision therapies.

Genomes of a Novel Group of Phages That Use Alternative Genetic Code Found in Human Gut Viromes

Metagenomics provides detection of phage genome sequences in various microbial communities. However, the use of alternative genetic codes by some phages precludes the correct analysis of their genomes. In this study, the unusual phage genome (phAss-1, 135,976 bp) was found after the de novo assembly of the human gut virome. Genome analysis revealed the presence of the TAG stop codons in 41 ORFs, including characteristic phage ORFs, and three genes of suppressor tRNA. Comparative analysis indicated that no phages with similar genomes were described. However, two phage genomes (BK046881_ctckW2 and BK025033_ct6IQ4) with substantial similarity to phAss-1 were extracted from the human gut metagenome data. These two complete genomes demonstrated 82.7% and 86.4% of nucleotide identity, respectively, similar genome synteny to phAss-1, the presence of suppressor tRNA genes and suppressor TAG stop codons in many characteristic phage ORFs. These data indicated that phAss-1, BK046881_ctckW2, and BK025033_ct6IQ4 are distinct species within the proposed Phassvirus genus. Moreover, a monophyletic group of divergent phage genomes containing the proposed Phassvirus genus was found among metagenome data. Several phage genomes from the group also contain ORFs with suppressor TAG stop codons, indicating the need to use various translation tables when depositing phage genomes in GenBank.

Gut Phageome-An Insight into the Role and Impact of Gut Microbiome and Their Correlation with Mammal Health and Diseases

The gut microbiota, including bacteria, archaea, fungi, and viruses, compose a diverse mammalian gut environment and are highly associated with host health. Bacteriophages, the viruses that infect bacteria, are the primary members of the gastrointestinal virome, known as the phageome. However, our knowledge regarding the gut phageome remains poorly understood. In this review, the critical role of the gut phageome and its correlation with mammalian health were summarized. First, an overall profile of phages across the gastrointestinal tract and their dynamic roles in shaping the surrounding microorganisms was elucidated. Further, the impacts of the gut phageome on gastrointestinal fitness and the bacterial community were highlighted, together with the influence of diets on the gut phageome composition. Additionally, new reports on the role of the gut phageome in the association of mammalian health and diseases were reviewed. Finally, a comprehensive update regarding the advanced phage benchwork and contributions of phage-based therapy to prevent/treat mammalian diseases was provided. This study provides insights into the role and impact of the gut phagenome in gut environments closely related to mammal health and diseases. The findings provoke the potential applications of phage-based diagnosis and therapy in clinical and agricultural fields. Future research is needed to uncover the underlying mechanism of phage-bacterial interactions in gut environments and explore the maintenance of mammalian health via phage-regulated gut microbiota.

Microviruses: A World Beyond phiX174

Two decades of metagenomic analyses have revealed that in many environments, small (∼5 kb), single-stranded DNA phages of the family Microviridae dominate the virome. Although the emblematic microvirus phiX174 is ubiquitous in the laboratory, most other microviruses, particularly those of the gokushovirus and amoyvirus lineages, have proven to be much more elusive. This puzzling lack of representative isolates has hindered insights into microviral biology. Furthermore, the idiosyncratic size and nature of their genomes have resulted in considerable misjudgments of their actual abundance in nature. Fortunately, recent successes in microvirus isolation and improved metagenomic methodologies can now provide us with more accurate appraisals of their abundance, their hosts, and their interactions. The emerging picture is that phiX174 and its relatives are rather rare and atypical microviruses, and that a tremendous diversity of other microviruses is ready for exploration.

The gut virome and the relevance of temperate phages in human health

Alterations in the gut virome impact human health. Bacteriophages, viruses that infect bacteria, dominate the gut virome and are mainly composed by virulent and temperate phages. While virulent phages exclusively replicate within and lyse their bacterial host's cell, temperate phages switch from an integrated state residing within their bacterial host's chromosome to an induced free virion state via an induction event. How often do these induction events occur and what are their implications on gut homeostasis? Here, we summarize the current knowledge of the gut virome based on metagenomics and present how the proportion of induced temperate phages varies amongst individuals, age, and disease states. Finally, we highlight the importance of building upon classical culture-dependent techniques and sequencing approaches to improve our understanding of temperate phages to enable their potential therapeutic use.

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.

Statistical optimization of a podoviral anti-MRSA phage CCASU-L10 generated from an under sampled repository: Chicken rinse

Introduction

The insurgence of antimicrobial resistance is an imminent health danger globally. A wide range of challenging diseases are attributed to methicillin-resistant Staphylococcus aureus (MRSA) as it is weaponized with a unique array of virulence factors, and most importantly, the resistance it develops to most of the antibiotics used clinically. On that account, the present study targeted the optimization of the production of a bacteriophage active against MRSA, and evaluating some of its characters.

Methods and results

The bacteriophage originated from a quite peculiar environmental source, raw chicken rinse and was suggested to belong to Podoviridae, order Caudovirales. It withstood a variety of extreme conditions and yield optimization was accomplished via the D-optimal design by response surface methodology (RSM). A reduced quadratic model was generated, and the ideal production conditions recommended were pH 8, glycerol 0.9% v/v, peptone 0.08% w/v, and 107 CFU/ml as the host inoculum size. These conditions led to a two-log fold increase in the phage titer (1.17x10¹² PFU/ml), as compared to the regular conditions.

Discussion

To conclude, statistical optimization successfully enhanced the output of the podoviral phage titer by two-log fold and therefore, can be regarded as a potential scale-up strategy. The produced phage was able to tolerate extreme environmental condition making it suitable for topical pharmaceutical preparations. Further preclinical and clinical studies are required to ensure its suitability for use in human.

Bacteriophages of the Order Crassvirales: What Do We Currently Know about This Keystone Component of the Human Gut Virome?

The order Crassvirales comprises dsDNA bacteriophages infecting bacteria in the phylum Bacteroidetes that are found in a variety of environments but are especially prevalent in the mammalian gut. This review summarises available information on the genomics, diversity, taxonomy, and ecology of this largely uncultured viral taxon. With experimental data available from a handful of cultured representatives, the review highlights key properties of virion morphology, infection, gene expression and replication processes, and phage-host dynamics.

Antibiotic resistance genes, mobile elements, virulence genes, and phages in cultivated ESBL-producing Escherichia coli of poultry origin in Kwara State, North Central Nigeria

The paucity of information on the genomic diversity of drug-resistant bacteria in most food-producing animals, including poultry in Nigeria, has led to poor hazard characterization and the lack of critical control points to safeguard public health. Hence, this study used whole genome sequencing (WGS) to assess the presence and the diversity of antibiotic resistance genes, mobile genetic elements, virulence genes, and phages in Extended Spectrum Beta Lactamase producing Escherichia coli (ESBL - E. coli) isolates obtained from poultry via the EURL guideline of 2017 in Ilorin, Nigeria. The prevalence of ESBL - E. coli in poultry was 10.5 % (n = 37/354). The phenotypic antibiotic susceptibility testing showed that all the ESBL- E. coli isolates were multi-drug resistant (MDR). The in-silico analysis of the WGS raw-read data from 11 purposively selected isolates showed that the isolates had a wide array of ARGs that conferred resistance to beta-lactam antibiotics, and 8 other classes of antibiotics (fluoroquinolones, foliate pathway antagonists, aminoglycoside, phenicol, tetracycline, epoxide, macrolides, and rifamycin). All the ARGs were in the bacterial chromosome except in two isolates where plasmid-mediated quinolone resistance (PMQR) was detected. Two isolates carried the gyrAp.S83L mutation which confers resistance to certain fluoroquinolones. The mobilome consisted of several Col-plasmids and the predominant IncF plasmids belonged to the IncF64:A-:B27 sequence type. The virulome consisted of genes that function as adhesins, iron acquisition genes, toxins, and protectins. Intact phages were found in 8 of the 11 isolates and the phageome consisted of representatives of four families of viruses: Myoviridae (62.5 %, n = 5/8), Siphoviridae (37.5 %, n = 3/8), Inoviridae (12.5 %, n = 1), and Podoviridae (12.5 %, n = 1/8). ESBL - E. coli isolates harboured 1-5 intact phages and no ARGs were identified on any of the phages. Although five of the isolates belonged to phylogroup A, the isolates were diverse as they belonged to different serotype and sequence types. Our findings demonstrate the high genomic diversity of ESBL - E. coli of poultry origin in Ilorin, Nigeria. These diverse isolates harbor clinically relevant ARGs, mobile elements, virulence genes, and phages that may have detrimental zoonotic potentials on human health.

Identification of conserved genomic signatures specific to Bifidobacterium species colonising the human gut

Bifidobacterium species are known for their ability to inhabit various habitats and are often regarded as the first colonisers of the human gut. In the present work, we have used comparative genomics to identify conserved genomic signatures specific to Bifidobacterium species associated with the human gut. Our approach discovered five genomic signatures with varying lengths and confidence. Among the predicted five signatures, a 1790 bp multi-drug resistance (MDR) signature was found to be remarkably specific to only those species that can colonise the human gut. The signature codes for a membrane transport protein belonging to the major facilitator superfamily (MFS) generally involved in MDR. Phylogenetic analyses of the MDR signature suggest a lineage-specific evolution of the MDR signature in bifidobacteria colonising the human gut. Functional annotation led to the discovery of two conserved domains in the protein; a catalytic MFS domain involved in the efflux of drugs and toxins, and a regulatory cystathionine-β-synthase (CBS) domain that can interact with adenosyl-carriers. Molecular docking simulation performed with the modelled tertiary structure of the MDR signature revealed the putative functional role of the covalently linked domains. The MFS domain displayed a high affinity towards various protein synthesis inhibitor antibiotics and human bile acids, whereas the C-terminally linked CBS domain exhibited favourable binding with molecular structures of ATP and AMP. Therefore, we believe that the predicted signature represents a niche-specific survival trait involved in bile and antibiotic resistance, imparting an adaptive advantage to the Bifidobacterium species colonising the human gut.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03492-4.

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.

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).

A Genome of Temperate Enterococcus Bacteriophage Placed in a Space of Pooled Viral Dark Matter Sequences

In the human gut, temperate bacteriophages interact with bacteria through predation and horizontal gene transfer. Relying on taxonomic data, metagenomic studies have associated shifts in phage abundance with a number of human diseases. The temperate bacteriophage VEsP-1 with siphovirus morphology was isolated from a sample of river water using Enterococcus faecalis as a host. Starting from the whole genome sequence of VEsP-1, we retrieved related phage genomes in blastp searches of the tail protein and large terminase sequences, and blastn searches of the whole genome sequences, with matches compiled from several different databases, and visualized a part of viral dark matter sequence space. The genome network and phylogenomic analyses resulted in the proposal of a novel genus "Vespunovirus", consisting of temperate, mainly metagenomic phages infecting Enterococcus spp.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

Mutualistic interplay between bacteriophages and bacteria in the human gut

Bacteriophages (phages) are often described as obligate predators of their bacterial hosts, and phage predation is one of the leading forces controlling the density and distribution of bacterial populations. Every 48 h half of all bacteria on Earth are killed by phages. Efficient killing also forms the basis of phage therapy in humans and animals and the use of phages as food preservatives. In turn, bacteria have a plethora of resistance systems against phage attack, but very few bacterial species, if any, have entirely escaped phage predation. However, in complex communities and environments such as the human gut, this antagonistic model of attack and counter-defence does not fully describe the scope of phage-bacterium interactions. In this Review, we explore some of the more mutualistic aspects of phage-bacterium interactions in the human gut, and we suggest that the relationship between phages and their bacterial hosts in the gut is best characterized not as a fight to the death between enemies but rather as a mutualistic relationship between partners.