Abstract
Phages, short for bacteriophages, are viruses that specifically infect bacteria and are the most abundant biological entities on earth found in every explored environment, from the deep sea to the Sahara Desert. Phages are abundant within the human biome and are gaining increasing recognition as potential modulators of the gut ecosystem. For example, they have been connected to gastrointestinal diseases and the treatment efficacy of Fecal Microbiota Transplant. The ability of phages to modulate the human gut microbiome has been attributed to the predation of bacteria or the promotion of bacterial survival by the transfer of genes that enhance bacterial fitness upon infection. In addition, phages have been shown to interact with the human immune system with variable outcomes. Despite the increasing evidence supporting the importance of phages in the gut ecosystem, the extent of their influence on the shape of the gut ecosystem is yet to be fully understood. Here, we discuss evidence for phage modulation of the gut microbiome, postulating that phages are pivotal contributors to the gut ecosystem dynamics. We therefore propose novel research questions to further elucidate the role(s) that they have within the human ecosystem and its impact on our health and well-being.
Introduction
The gastrointestinal tract (GIT) of humans and many other animals hosts a complex ecosystem inhabited by a plethora of different microorganisms, that include bacteria, fungi, archaea, protozoa, and viruses (). Multiple factors affect gut microbial communities and contribute to the complexity of this ecosystem. These factors include (but are not limited to) the anatomy of the GIT, peristaltic movements, the mucus layer and its shedding, host-produced compounds (e.g. bile acids or gastric juice), the constant influx of new microorganisms and nutrients through diet, and the host immune system (; ).
The microbial communities (microbiomes) in the gut are involved in regulating many aspects of the host’s well-being through the mediation of nutrient absorption, synthesis of vitamins and neurotransmitters, the development and modulation of the immune system, and modifying resistance against pathogens, among others (). Therefore, it is unsurprising that disturbances in this ecosystem, reflected by alterations in the microbial communities that inhabit it, have been connected to multiple diseases, from gut inflammation to neurological disorders (; ).
To date, the majority of the reported microbial impacts on the gut ecosystem and host health have been connected to the bacterial component of the microbiome (). Recently, viruses that infect bacteria, namely bacteriophages or, phages for short, have been gaining attention as potential modulators of the gut ecosystem due to their ability to affect bacterial communities. Most peculiarly, the influence of phages on the gut ecosystem seems to extend beyond their direct impacts on bacterial populations, extending to modulation of the host immune system (). Furthermore, alterations in the gut phage population have been connected to gastrointestinal diseases on multiple occasions (; ; ; Zuo et al., 2018; ; ), highlighting their contribution to gastrointestinal health. In this review, we postulate that phages play a major role in the gut ecosystem dynamics through an intricate network of interactions with both the gut bacterial community and the host immune system.
Phages: An Overview
Phages are obligate parasites that require a bacterial host for reproduction and are the most abundant and diverse biological entities on earth (). Structurally, the majority of phages are composed of a nucleic-acid genome packaged inside a protein shell (i.e., “capsid”; Figures 1A, B). Phage capsids are highly variable, both in size and morphology (i.e., polyhedral, filamentous, or pleomorphic). Some phages present an outer lipid membrane in addition to their protein capsid, while others only have the lipid membrane ().
Figure 1
Phage genomes are variable in size (ranging between ~3.5 kb and ∼540 kb) and are composed of either single or double-stranded DNA (ssDNA, dsDNA), or RNA (ssRNA, dsRNA) (; ; ; ; ; ; Yang et al., 2021; Zhai et al., 2021). Consequently, the genomes of phages infecting different hosts seldom share sequence similarities (). Additionally, phage genomes present a mosaic structure as a result of recombination events with bacteria () and other phage (). This genome mosaicism is characterized by highly similar sequences juxtaposed with sequences with which they share no apparent similarities (; ; ). The extreme variability and mosaicism of phage genomes significantly complicate taxonomic classification, which was originally performed using shared, sequence agnostic features (e.g., capsid morphology and type of nucleic acid composing the genome) (). For these reasons, virus taxonomy is under review due to new classification methods that include genomic sequences, genome organization, and host range (). (Up-to-date viral taxonomy information can be found at https://ictv.global/vmr/). Nonetheless, the original nomenclature is still commonly used, and phages are often distinguished based on the type of nucleic acids they carry (e.g., dsDNA) and structural conformation (e.g., Caudovirales).
Phages are generally classified as virulent and temperate phages based on the life cycle they follow. After the recognition and subsequent attachment to a specific receptor on the bacterial cell surface, the phage delivers its genome into the bacterial cell. Here, the phage genome is replicated and expressed using host cellular resources before new complete viral particles (virions) are assembled and released from the bacterial cell. Newly assembled virions can be released by phage-mediated lysis of the bacterial cell in the lytic cycle, which is common to most known phages, or through a bacterial secretion apparatus in the chronic cycle, of the filamentous Inoviridae phages [Figure 1B; reviewed by ]. By contrast, a temperate phage can undertake a lysogenic cycle in which, following the insertion of its genome into the bacterial cell, the phage enters a quiescent state. During this state, the phage genome, or prophage, is replicated with the host genome, either as a result of being integrated into the bacterial chromosome or as an extrachromosomal molecule. The lysogenic cycle typically ends when a specific stimulus [e.g., the bacterial SOS response ()] initiates either the lytic or chronic cycle and promotes the production of new virions and their release from the bacterial cell. However, prophages may also become defective and lose the ability to be induced or to excise from the bacterial chromosome (; ).
Phages in the Human Gastrointestinal Tract
The density of phages increases through the gastrointestinal tract from the small to the large intestine (). The density of phages in the large intestine ranges between 108 and 8×1010 phage virions per gram of feces, measured by phage particles count (; ) and estimation of viral genome numbers in feces of healthy adults (). In a mixed cohort of healthy and unhealthy individuals, phages appeared to make up the vast majority (97.7%) of gut viral genomes, with eukaryotic (2.1%), and archaeal viruses (0.1%) accounting for the remainder (). Notably, approximately 90% of the phage component was unclassified, while the remainder were non-enveloped DNA phages, belonging to the dsDNA order Caudovirales or the ssDNA families of Microviridae and Inoviridae (Figure 1) ().
The genomic diversity of gut phages remains largely unknown (; ; ) due to different impediments such as 1) the lack of a universal marker gene (i.e. an equivalent to the bacterial 16S rRNA gene) for targeted phage amplicon analyses and taxonomic assignments, 2) the high variability of phage genomes, 3) and the difficulty of cultivating gut phages (). The impact of these limitations is illustrated by the fact that the two most abundant fecal phage clades, crAssphage and Gubaphage, were only identified in 2014 () and 2021 (), respectively. Both crAssphage (Cross Assembly phage) and Gubaphage (Gut Bacteroidales phage) are dsDNA phages characterized by long genomes (~ 97 and ~80 kb, respectively) which infect bacteria of the genus Bacteroides (; ). This association has been predicted in silico for Gubaphage (), but it has been demonstrated in vitro for crAssphage (). In vitro cultivation has also revealed that crAssphage has a Podovirus-like conformation(Figure 1B) and a temperate life cycle (). The discovery of crAssphage has led to the identification of numerous crAss-like phages in the human gut microbiome ().
RNA phages are rare, if not absent, in the gut. Instead, the majority of viral RNA genomes within the gastrointestinal tract originate from plant viruses acquired in diet (Zhang et al., 2006; ). However, it remains possible that the low abundance of RNA gut phage genomes observed is due to the limited numbers of studies that have analysed the RNA viral component, and the scarcity of reference phage genomes for contig identification (Zhang et al., 2019; ). Moreover, the abundance of dsDNA Caudovirales phages and ssDNA phages(Figure 1) found in the human GIT, could be due to methodological biases in virion extraction and metagenomic sequencing that favour their identification at the expense of other phages [e.g. chloroform extraction for Caudovirales () or multiple displacement amplification for ssDNA phages ()].
The phage composition of the gut has been reported to remain stable for up to 1 year period in healthy adults (), and alterations have been associated with gastrointestinal diseases, such as Clostridioides (formerly Clostridium) difficile infections (CDI) and inflammatory bowel diseases (IBD) (; ; ; Zuo et al., 2018; ; ). Notably, different diseases (i.e., CDI and norovirus-associated diarrhoea, or ulcerative colitis and Crohn’s disease) were associated with specific gut phage compositions (; Zuo et al., 2018). More specifically, compared to healthy subjects, CDI patients showed an increased abundance of Caudovirales phages and a reduction in their diversity, richness, and eveness, while norovirus-associated diarrhoea patients showed, alongside reduction in Caudovirales richness and diversity, also a reduction in their abundance (Zuo et al., 2018). Furthermore, Crohn’s disease patients presented an increased Caudovirales richness when compared to healthy controls, while ulcerative colitis patients did not (), overall suggesting a connection between phage occupancy and the health state of the gut.
Phage Interactions in the Human GIT and Population Dynamics
The alterations of the gut phage composition in association with different diseases suggest a potential ecological influence of phages on the gut ecosystem. This influence has been attributed to their ability to interact with and modulate the gut bacterial community and the host immune system (), as both affect the homeostasis of the GIT (; ).
The strongest contribution that phages exert on the shape of the gut ecosystem is arguably through the modulation of the gut bacterial community. This is dependent on the phage life cycle. During the lytic cycle, phage predation follows predator-prey-like dynamics (). As such, there are strong selective pressures for bacteria to evolve resistance mechanisms against lytic phages (; ; ; ), and for the phages to develop strategies to evade those mechanisms (; ). This creates an arms race in which bacteria and phages are constantly evolving mechanisms to prevent and promote infection (), respectively.
However, interactions between phages and bacteria extend beyond predator-prey dynamics, as lysogenic phage infection has been suggested to potentially have beneficial effects on the bacterial host (). While a prophage remains dormant in the bacterial cell, its survival is directly linked to that of the host. Therefore, it is evolutionarily advantageous for the prophage to contribute to the host survival (), in what could be arguably referred to as a mutualistic interaction (). Such mutualism is exemplified by virions carrying genes that have no direct impact on the phage life cycle but can enhance the fitness of the bacterial host, termed “morons” () in a phenomenon called lysogenic conversion. Common morons include, for example, bacterial virulence () or metabolic genes (; Zuppi et al., 2020). Alternatively, morons may provide resistance from other infecting virions, in the superinfection exclusion phenomenon (). Interestingly, genes that enhance bacterial fitness have also been identified in virulent phages (i.e. phages that follow only the lytic cycle) (; ; ), suggesting that they may also promote limited bacterial survival or ‘cultivation’ to favour their reproduction.
In addition to modulating bacterial communities, phages influence the gut ecosystem by interacting directly with the immune cells and thereby modulating host immune activity [reviewed in (; ; )]. Phage particles can cross the epithelial barrier through a process known as transcytosis () and interact directly with the mammalian immune cells. T4 phage reduced production of reactive oxygen species from peripheral blood polymorphonuclear leukocytes and reduced NF-kB activity in mouse model (; ) demonstrating that phage particles can exert a dampening effect on the mammalian immune system. By contrast, phage particles were also shown to stimulate an immune response in mice via recognition from Toll-like Receptor 9 (TLR9) (), suggesting an ambivalent effect of phages on the mammalian immune response. Importantly, studies investigating the potential uses of phage therapy in treating bacterial infections demonstrated that phage particles can stimulate the production of specific neutralizing antibodies that could dampen phage activity (; ), showing a reciprocal influence between phages and the mammalian immune system in the GIT. Interestingly, specific phage-encoded proteins have been shown to modulate the interaction between bacteria and the immune system. The tail adhesin Gp12 was shown to bind bacterial lipopolysaccharide [or LPS, a bacterial endotoxin ()] and impede its recognition from the human immune system (). Additionally, the presence of Immunoglobulin-like domains on phage capsid was shown to mediate the binding to the intestinal mucin layer, impeding its colonization by bacterial communities and providing a non-host-derived immunity ().
Phages’ contributions to gut ecosystem dynamics are further modulated by their biological interactions with other phages upon infecting the bacterial host. These interactions can either be antagonistic or cooperative (), affecting the efficacy of their reproduction and, therefore, the phage-mediated modulation of bacterial communities and host immune system. Phage antagonism occurs when two phages with a common host, either compete for receptor adsorption () or impede each other’s infection and life cycle through superinfection exclusion, described above. Also, cooperative interactions between phages have been observed in overcoming bacterial defense mechanisms, such as CRISPR-Cas systems (). Moreover, Bacillus phages were shown to present a communication mechanism used to determine which life cycle to follow (; ), highlighting the surprisingly social nature of some phages.
As a result of these multiple interactions, the human GIT is likely to be home to different population dynamics between phages and bacteria. Population dynamics is the study of the changes of a population in size and structure over time, and the factors behind them. These changes are described through mathematical models (Box 1).
Box 1
Models proposed to describe the phage-bacteria population dynamics can be divided into two groups:
Group 1) models that are characterized by a low variability and diversity of phage and bacterial populations [e. Arms-race Dynamic () and Piggyback-the-Winner dynamic ()];
The Arms-race dynamic (ARD) is characterized by competition for survival driven by predator-prey interactions, with the bacterial population developing counter-infection defenses, while the phage population develops methods to evade these defenses [Reviewed in ()]. This leads to an arms race that causes the extinction of not “up to date” phage and bacterial populations, while the competent ones flourish (), resulting in few dominant phage and bacterial populations.
The Piggyback-the-winner dynamic (PtW) is driven by mutualistic interactions that occur in the lysogenic life cycle. More specifically, describing prophages that contribute to the survival of hosts through lysogenic conversion and superinfection exclusion. This dynamic is characterized by low variability and diversity of the phage and bacterial populations (; ) and leads to a few populations out-competing the others thanks to the advantages provided by this mutually beneficial interaction. PtW is the only proposed model of phage-bacteria population dynamics that is connected to the lysogenic cycle.
Group 2) models that are characterized by high phage and bacterial population variability and diversity [i.e. Fluctuating selection dynamic () and Kill-the-winner dynamic ()]. These models are driven by “negative frequency-dependent selection”, a form of natural selection in which the fitness of a genotype is inversely proportional to its frequency ().
The Fluctuating-selection dynamic (FSD) results from the bacterial fitness costs for developing phage-resistance mechanisms being disadvantageous in an environment in which multiple bacterial species are competing for resources. Specifically, effective phage-resistance mechanisms lead to a decrease in the number of infecting virions. As predation decreases, phage-resistant bacterial populations will be outcompeted by bacterial populations which did not invest in defense mechanisms. Ultimately, this leads to a switch in the bacterial communities from phage-resistant to non-phage-resistant bacteria. In this newly permissive environment, phage predation increases, making the niche favorable again for the growth of phage-resistant bacteria populations, and the cycle begins anew ().
The Kill-the-winner dynamic (KtW) occurs when the abundance of the “winning” bacterial species (i.e., the most competitive) is controlled by phage predation. This allows the coexistence of multiple bacterial and phage populations by limiting the expansion of the most competitive populations ().
Alternatively, these population dynamics models can be grouped based on the prevalent phage life cycle and modalities of interaction between phages and bacteria populations in the environment, namely lytic or predator-prey interactions (ARD, KtW, and FSD) and lysogenic or mutualistic interactions (PtW). The type of life cycle followed by phages in an environment is inferred based on the Virus-to-Microbe ratio (VMR), as VMR lower than one suggests the prevalence of the lysogenic cycle, while higher VMR suggests the prevalence of the lytic cycle.
Population dynamics are strictly dependent on the environment in which they take place. In GIT, factors influencing the population dynamics include pH and bile acid levels, structural conformations, mucin layer, oxygen levels, and nutrient availability; all varying across the length of the GIT. This variability influences gut bacteria () and, in certain cases, phage populations (; ; ). Mirroring these physiological gradients, phage-bacterial population dynamics have also been reported to vary throughout the gut (). For example, the prevalent population dynamic in the lumen of the colon is thought to be the PtW dynamic, as suggested by the low VMR reported in feces (≥ 1:1) (; ). In the mucin layer of the colon, the high VMR (~20:1) () suggests lytic cycle-driven dynamics (e.g. ARD, KtW, or FSD) (Box 1). This has been connected to a reduction of bacterial cell densities in the mucin layer compared to the lumen of the colon (). This is consistent with what was reported by Knowles et al., who showed that, in marine environments, the prevalent life cycle in a phage population was connected to bacterial densities (). Similarly, it was suggested that the prevalent population dynamic changes toward the distal part of the colon, as an increasing number of stressors reduce bacterial densities and promote induction of prophages, resulting in a shift toward lytic cycle-driven dynamics (ARD, KtW, or FSD) ().
Importantly, the population dynamics in the human GIT are likely to differ from the models used to describe them, as these were originally proposed for other environments [i.e. in vitro () and aquatic environments (; )], which are significantly different from the human GIT (). This has been exemplified by Park et al. who reported that, in patients with recurrent CDI, an increase in gut phageome (i.e., the phage component of the microbiome) diversity was accompanied by a reduction of bacterial diversity, contrary to predictions from the proposed models (Box 1) (). Despite being unclear whether these changes in microbial diversity might be the cause or the results of CDI, the increase of phage diversity not being coupled by an increase in bacterial diversity suggests that other factors in the human GIT contribute to these dynamics. These modifiers of the phage-bacteria population dynamics may include the host immune system, as it has been shown to both strongly influence and be influenced by both phages and bacteria. This hypothesis is supported by observations from Clooney and colleagues who identified an increase in induced temperate phages in IBD patients and connected the phage induction to an increase in inflammation biomarkers, such as reactive oxygen species ().
The Ecological Role of Phages in the Gut Ecosystem
The organisms in macroscopic ecosystems are connected by physical interactions and by participation in the flux of materials and energy (). Similar to macroscopic ecosystems, Hsu and colleagues illustrated the connection between different organisms in the gut by showing that, in gnotobiotic mice, phage-mediated lysis of bacterial populations had cascading effects on bacterial populations that were not directly targeted by phage predation (). Fluctuations in the abundance of individual species have distinct consequences on the ecosystem depending on their ecological role (). In the human GIT, as resident microbial communities exert a strong influence on the host immune system, phage-mediated alterations of these communities could promote cascading effects potentially leading to profound ecosystem alterations and disease () (Figure 2).
Figure 2
Phages alter the gut commensal bacterial composition through their parasitic activity. This may lead to the development of gastrointestinal diseases through altered colonization resistance and proliferation of opportunistic or obligate pathogens in vacated niches (). Furthermore, phage-mediated depletion of commensal bacterial populations may interfere with the production of immunomodulatory compounds. IBD patients were shown to have reduced levels of Faecalibacterium prausnitzii () and increased levels of its phage (). Notably, F. prausnitzii is a major contributor of colonic butyrate (), which is an anti-inflammatory short-chain fatty acid (SCFA) (). Therefore, increased phage predation may indirectly contribute to IBD inflammation by reducing the amounts of butyrate produced by F. prausnitzii.
Similarly, Parkinson’s disease patients were observed to have increased levels of lytic Lactococcus phages and depleted levels of Lactococcus bacteria, the latter of which regulate gut permeability and produce dopamine. These factors are implicated in Parkinson’s disease pathogenesis, suggesting that phage predation within the GIT could contribute to the development of this disease (). More broadly, these data demonstrate that phage-mediated depletion of commensal bacterial populations may lead to decreased production of bacterially-derived immunomodulatory compounds, such as SCFAs, thereby modulating gut inflammation. By contrast, phage predatory activity may contribute to GIT health by keeping bacterial abundances within tolerable levels. Consistent with this, phages were shown to provide a non-host derived immunity by binding to the mucin layer surrounding the intestinal epithelium via an Ig-like protein in the capsid. Through their lytic activity and predation in the mucin layer, phages prevent its colonization from bacterial cells () and activation of the immune system ().
Temperate phages are capable of strongly contributing to the bacterial host’s virulence and fitness by providing virulence genes, such as the phage-encoded toxins or immune evasion genes (). In this way, temperate phages can indirectly exert a pro-inflammatory effect in the GIT and cause alterations to the gut ecosystem. Prophages may also increase the fitness of commensal bacteria through superinfection exclusion or lysogenic conversion, contributing to the maintenance of a healthy gut environment. Importantly, the ecological impact of temperate phages goes beyond increasing the fitness of the bacterial cells they infect. Prophages can be induced, leading to the initiation of the lytic cycle and the death of their bacterial host, potentially changing the population dynamics and the impact on the GIT. Therefore, stimuli that cause prophage induction can affect their influence on the gut ecosystem. Following this reasoning, the induction of prophages has been proposed to participate in the establishment of positive feedback loops of GIT inflammation in humans ().
Lin et al. theorized that prophage induction promotes gut inflammation by spreading integrated virulence factors. Gut inflammation augments intestinal permeability, increasing the luminal oxygen level. This in turn promotes prophage induction through a mechanism that involves oxidative stress establishing a positive feedback loop (). Similarly, suggested that phage-mediated lysis of bacterial cells contributes to IBD associated gut inflammation by promoting the release of Pathogen Associated Molecular Patterns (PAMP), such as bacterial DNA, lipopolysaccharide, and peptidoglycan (), and their exposure to Pattern Recognition Receptors (PRRs) on gastrointestinal epithelial cells (e.g. Toll-like Receptors or TLR) (). This triggers the production of Reactive Oxygen and Reactive Nitrogen Species (ROS and RNS, respectively) which promote a host inflammatory response () that can indirectly increase prophage induction, as suggested previously. Moreover, ROS and RNS can promote prophage induction directly (; ), suggesting yet another mechanism through which phages may promote a positive inflammatory feedback loop (Figure 3).
Figure 3
The intense predatory activity of phages can have profound influences on the microbial communities and the gut ecosystem at large. Therefore, phages have been compared to apex predators in macroscopic environments (), namely predators whose activity significantly shapes the predated ecosystem. Nonetheless, the ambivalence of the interactions of temperate phages with bacterial communities, which can occur both in a predator-prey and mutualistic manner, suggests that their ecological role is more complex. Intriguingly, the modality of interaction with bacterial communities, and therefore the interaction with the host immune system and the population dynamics, appear to shift in response to environmental changes. This suggests that the ecological role of phages changes depending on the surrounding environment.
The known complexity of the direct and indirect impacts of phages on the gut ecosystem suggests that their ecological role has no direct comparison in the macroscopic environments. However, it is clear that phages strongly contribute to the temporal patterns, directionality, frequency, and magnitude of population changes within the gut microbial community, with profound repercussions on the gut ecosystem and health.
Future Research
Despite the increasing interest in gut phages and their roles in the gut microbiome, research into this field is still in its infancy. Arguably, a key question is the role of the phageome in the gut microbial ecosystem at large. This can be further divided into more specific questions regarding 1) consistency or variability of gut phageomes across individuals and populations; 2) host-specificity and influence of gut phages on their bacterial host; 3) direct influence of gut phageome on the human immune system; 4) shifts in the gut population dynamics either driven by or contributed by individual phage or phages(Figure 4).
Figure 4
To begin addressing these questions, it is important to survey the diversity and stability of gut phage composition across populations (rural, urban and international; healthy and unwell). This would allow the determination of a “core phageome” (i.e., a set of phages consistently identified in the gut microbiomes of individuals) and those that are positively associated with gut health. Identifying such a set of phages would allow a more targeted analysis, facilitating the understanding of the role of phages in the gut microbiome. Current attempts to address this question have reported highly subject-specific phageomes with only a minimal proportion of phage genomes (less than 1%) shared across more than half of the studied population (; ). As phage genomes are highly mosaic and variable, a “core phageome” could be better identified by focusing on shared phage functions between individuals, rather than attempting to identify common phage taxa or genomes. While most phageome genes are involved in the phage life cycle or encode structural proteins (), morons are involved in bacterial phenotypes and ecosystem dynamics in a variety of ways. Targeted identification and annotation of morons could therefore provide an additional view to how conserved or individualized the phageomes are across individuals and populations (i.e., core functional phageome). Since homology-based methods can only annotate a minority of phage gene clusters (), functional characterization of phageomes will benefit from recent neural network-based protein structure and function prediction methodologies (; ; ; ). Specific morons to evaluate are genes that increase the fitness of bacterial hosts. These include, but are not limited to, antibiotic resistance genes [from databases such as CARD () or ARDB ()], virulence genes [VFDB (), MvirDB (Zhou et al., 2007)], genes implicated in bacterial motility and other auxiliary metabolic genes (AMGs); specific tools to assess metagenomic AMGs include DRAM-v () and VIBRANT (). Additionally, other databases and tools developed for cataloging or analyzing prokaryotic genomes, such as DRAM (prokaryotic version) and METABOLIC (Zhou et al., 2020), could enable the identification of additional morons that were overlooked by the phage-specific tools or databases above. Genes coding for structural proteins can also be included in the set of phage-encoded genes that impact the gut ecosystem. An example of this is provided by the tail protein Gp12 () or capsid Ig-like domain (), structural proteins with potential anti-inflammatory outcomes.
Identifying such genes in the phageome also allows the determination of the impact that phages have on the gut ecosystem. The presence of lysogenic conversion morons in a phage genome describes its potential influence on its bacterial hosts, while the presence of genes such as the ones coding for tail protein Gp12 or the capsid Ig-like domains describe its potential influence on the human immune system. Nonetheless, the role of phages in the gut ecosystem is not entirely dependent on the morons or the immune-modulation genes they encode, but it is also dependent on the different life cycles they can follow and the bacterial hosts they infect. The influence that phages have on the bacterial host is in fact strictly dependent on their lifestyle, as temperate phages can establish a mutualistic relationship with their host, while virulent phages establish a predator-prey relationship. In addition to this, the ecosystemic outcomes of the phage influence on the bacterial host depend on the ecological role of that host. For example, a temperate phage increasing the fitness of the bacterial host by carrying AMGs, or a virulent phage infecting a bacterial population, will impact the gut ecosystem differently depending on whether the bacterial host is a commensal or a pathogenic bacterium. Moreover, phages able to infect multiple hosts can favour the spread of virulence genes with negative outcomes on gastrointestinal health. To investigate the lifestyle and host-specificity of phages, multiple tools have been developed. To determine the lifestyle of the sequenced phages, tools such as BACPHLIP (), Deephage (Wu et al., 2020), or PHACTS () can be used. BACPHLIP determines the phage lifestyle based on the presence or absence of genes associated with the lysogenic cycle in its genome. Deephage and PHACTS predict the phages’ lifestyle by identifying, in their genome, nucleotide features shared with a set of phages with known life cycle. To determine the host specificity of the sequenced phages, different tools have been proposed, and have been collectively described in ref (). These approaches can be divided into 1) alignment-dependent approaches, which determine the host specificity of a phage by identifying the sequences it shares with either bacteria [e.g., SpacePHARER (Zhang et al., 2021)] or phages with known hosts [e.g., RaFAH ()] in databases; 2) alignment-free approaches, which determine the host specificity based on genomic convergence between a phage and its host through machine learning models [e.g., WIsH () or PHP ()]; 3) integrative approaches that combine alignment-dependent and alignment-free approaches [e.g., VirHostMatcher Net () or PHISDetector (Zhang F et al., 2020)].
These approaches present significant impediments, as determining the phage lifestyle based on the presence or absence of specific genes is strictly dependent on the completeness of the assembled contig. Differently, tools that, to assess the phage lifestyle, rely on nucleotide similarities between the query and a set of phage genomes might exclude phages that significantly differ from that set. Furthermore, the tools to determine the host specificity present impediments based on the prediction approach. Alignment-based approaches strongly depend on the reference databases, while alignment-free approaches have a high occurrence of false positives, and the results often need to be confirmed through statistical analysis. Integrative approaches appear to overcome these issues but are still relatively new and their reliability still needs to be assessed (). Moreover, none of the different approaches that predict the host specificity describe the efficacy of the infection, i.e., the duration of the lytic cycle and the number of virions produced. The efficacy of phage infection has been reported to vary depending on the bacterial host infected (), with potentially different outcomes on the ecosystem. While numerous tools are being developed to determine the phage host specificity, a similar effort is not being made in the prediction of the phage lifestyle. The currently available tools either present significant impediments (i.e., BACPHLIP requires the query contigs to represent complete genomes) or are not user-friendly (i.e., Deephage). PHACTS is a web-based tool, which is not suitable for metagenomic analysis of large data. Therefore, the development of bioinformatic tools to predict the phage lifestyle that are more apt at metagenomic analysis would strongly contribute to optimizing the study of gut phageome.
Alternatively, in addition to in silico analyses of metagenomes, targeted in vitro work could contribute to the validation of computationally derived predictions of phage-host interactions, phage lifestyle, phage host-specificity, and efficacy of phage infection. A significant impediment to in vitro cultivation of gut phages is posed by the lifestyle of most gut bacteria, which are facultative or obligate anaerobes (). Nonetheless, different approaches for in vitro cultivation have been provided to overcome this issue (; ; ). Moreover, in vitro and in vivo description of phage interactions with human immune cells holds promising potential to study their impact on the host immune system in the GIT, as suggested by previous results (; ; ).
Changes in phageomes and phage-induced changes in bacteria may, in extreme cases, reorganize the gut ecosystem through shifts in population dynamics. Population dynamics in an environment have been usually described by determining the alpha-diversity of phage and bacterial populations and the VMR, an indication of the prevalent phage life cycle. However, analyzing gut ecosystems using sequencing-based technologies overlooks both changes in absolute virion and bacterial cell abundances as well as different microenvironments within the gut. To overcome the former limitation, measuring absolute viral abundance using a separate assay, such as quantitative microscopy () or phage spike-ins (), would help assess changes in population dynamics. To further address this issue, researchers should consider conducting, alongside purification and metagenomic sequencing of phages, also metagenomic sequencing of the bacterial component of the same sample. This can be performed either through direct analysis of the isolated bacterial component following viral enrichment. Such analysis will capture and distinguish between phages undergoing the lytic cycle and dormant prophages integrated into bacterial genomes. Together with assessing the alpha diversity of phage and bacterial communities in a gut microbiome, identifying the prevalent life cycle would help identify the population dynamics between the two (). Finally, integration of phageome and bacteriome analysis also enables the construction of bacterial CRISPR arrays (bacterial immune memory against phages) which can be compared to the phage genome repertoire. Evaluating the composition and evolution of bacterial CRISPR spacers could provide additional phage-host pair prediction. Together with other tools specifically designed for evaluating bacterial anti-phage defense systems, such as PADS Arsenal (Zhang Y et al., 2020), such analyses will provide insights into the phage-bacteria arms race, other interactions, and gut ecosystem dynamics at large.
Overall, improving computational phageome analyses, integration of phage virions and bacterial metagenomes with in vitro validations and follow-up studies will lead to a comprehensive understanding of the gut phageome, ultimately revolutionizing the way we think about our inner ecosystems.
Conclusions
Phages are pivotal components of the human gastrointestinal tract and changes in the phage composition and abundance have been associated with multiple gastrointestinal diseases. They are major drivers of bacterial evolution and important modulators of the host immune system, thereby contributing to the gut ecosystem dynamics. Despite much is yet to be understood about the gut phageome, focusing on the interactions between phages, gut bacteria and the host immune system holds promising potential for a more complete understanding of the gut microbiome and its connection with human health.
Funding
This work was funded by The Royal Society Te Apārangi Marsden Fund (MFP-UOA1901).
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Author contributions
MZ wrote the manuscript with comments from all other authors. JO’S and TV supervised the work. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors thank Prof. Forest Rohwer (San Diego State University, San Diego, California, USA) for his comments and suggestions on the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AlcockB. P.RaphenyaA. R.LauT. T. Y.TsangK. K.BouchardM.EdalatmandA.et al (2020). CARD 2020: Antibiotic Resistome Surveillance With the Comprehensive Antibiotic Resistance Database. Nucleic Acids Res.48, D517–D525. doi: 10.1093/nar/gkz935
2
AvraniS.SchwartzD. A.LindellD. (2012). Virus-Host Swinging Party in the Oceans. Mob. Genet. Elements2, 88–95. doi: 10.4161/mge.20031
3
BaekM.DiMaioF.AnishchenkoI.DauparasJ.OvchinnikovS.LeeG. R.et al (2021). Accurate Prediction of Protein Structures and Interactions Using a Three-Track Neural Network. Science373, 871–876. doi: 10.1126/science.abj8754
4
BarrJ. J.AuroR.FurlanM.WhitesonK. L.ErbM. L.PoglianoJ.et al (2013). Bacteriophage Adhering to Mucus Provide a Non-Host-Derived Immunity. Proc. Natl. Acad. Sci. U.S.A.110, 10771–10776. doi: 10.1073/pnas.1305923110
5
BernardC.LiY.LopezP.BaptesteE. (2021). Beyond Arbitrium: Identification of a Second Communication System in Bacillus Phage Phi3t That May Regulate Host Defense Mechanisms. ISME J.15, 545–549. doi: 10.1038/s41396-020-00795-9
6
Bondy-DenomyJ.QianJ.WestraE. R.BucklingA.GuttmanD. S.DavidsonA. R.et al (2016). Prophages Mediate Defense Against Phage Infection Through Diverse Mechanisms. ISME J.10, 2854–2866. doi: 10.1038/ismej.2016.79
7
BreitbartM.BonnainC.MalkiK.SawayaN. A. (2018). Phage Puppet Masters of the Marine Microbial Realm. Nat. Microbiol.3, 754–766. doi: 10.1038/s41564-018-0166-y
8
BronsteinJ. L. (2015). Mutualism. Oxford, UK: Oxford University Press.
9
Camarillo-GuerreroL. F.AlmeidaA.Rangel-PinerosG.FinnR. D.LawleyT. D. (2021). Massive Expansion of Human Gut Bacteriophage Diversity. Cell184, 1098–1109.e9. doi: 10.1016/j.cell.2021.01.029
10
Carroll-PortilloA.LinH. C. (2019). Bacteriophage and the Innate Immune System: Access and Signaling. Microorganisms7, 1–12. doi: 10.3390/microorganisms7120625
11
ChenL.YangJ.YuJ.YaoZ.SunL.ShenY.et al (2005). VFDB: A Reference Database for Bacterial Virulence Factors. Nucleic Acids Res.33, D325–D328. doi: 10.1093/nar/gki008
12
ChevallereauA.MeadenS.FradetO.LandsbergerM.MaestriA.BiswasA.et al (2020). Exploitation of the Cooperative Behaviors of Anti-CRISPR Phages. Cell Host Microbe27, 189–198.e6. doi: 10.1016/j.chom.2019.12.004
13
ChowJ.LeeS. M.ShenY.KhosraviA.MazmanianS. K. (2010). Chapter 8 - Host–Bacterial Symbiosis in Health and Disease. Adv. Immunol.107, 243–274. Academic Press. doi: 10.1016/B978-0-12-381300-8.00008-3
14
ClarkeB. (1962). “Balanced Polymorphism and the Diversity of Sympatric Species,” in Taxonomy and Geography. Oxford, UK: Systematic Association.
15
ClooneyA. G.SuttonT. D. S.ShkoporovA. N.HolohanR. K.DalyK. M.O'ReganO.et al (2019). Whole-Virome Analysis Sheds Light on Viral Dark Matter in Inflammatory Bowel Disease. Cell Host Microbe26, 764–778.e5. doi: 10.1016/j.chom.2019.10.009
16
CocletC.RouxS. (2021). Global Overview and Major Challenges of Host Prediction Methods for Uncultivated Phages. Curr. Opin. Virol.49, 117–126. doi: 10.1016/j.coviro.2021.05.003
17
CornuaultJ. K.PetitM.MariadassouM.BenevidesL.MoncautE.LangellaP.et al (2018). Phages Infecting Faecalibacterium Prausnitzii Belong to Novel Viral Genera That Help to Decipher Intestinal Viromes. Microbiome6, 65. doi: 10.1186/s40168-018-0452-1
18
CoutinhoF. H.Zaragoza-SolasA.López-PérezM. (2021). RaFAH: A Superior Method for Virus-Host Prediction. Patterns2, 100274. doi: 10.1016/j.patter.2021.100274
19
CumbyN.DavidsonA. R.MaxwellK. L. (2012). The Moron Comes of Age. Bacteriophage2, 225–228. doi: 10.4161/bact.23146
20
De SordiL.LourençoM.DebarbieuxL. (2019). The Battle Within: Interactions of Bacteriophages and Bacteria in the Gastrointestinal Tract. Cell Host Microbe25, 210–218. doi: 10.1016/j.chom.2019.01.018
21
DiardM.BakkerenE.CornuaultJ. K.MoorK.HausmannA.SellinM. E.et al (2017). Inflammation Boosts Bacteriophage Transfer Between Salmonella Spp. Science355, 1211–1215. doi: 10.1126/science.aaf8451
22
DionM. B.OechslinF.MoineauS. (2020). Phage Diversity, Genomics and Phylogeny. Nat. Rev. Microbiol.18, 125–138. doi: 10.1038/s41579-019-0311-5
23
Domingo-CalapP.Mora-QuilisL.SanjuánR. (2020). Social Bacteriophages. Microorganisms8, 533. doi: 10.3390/microorganisms8040533
24
DonaldsonG. P.LeeS. M.MazmanianS. K. (2016). Gut Biogeography of the Bacterial Microbiota. Nat. Rev. Microbiol.14, 20–32. doi: 10.1038/nrmicro3552
25
DraperL. A.RyanF. J.SmithM. K.JalankaJ.MattilaE.ArkkilaP. A.et al (2018). Long-Term Colonisation With Donor Bacteriophages Following Successful Faecal Microbial Transplantation. Microbiome6, 1–9. doi: 10.1186/s40168-018-0598-x
26
DutilhB. E.CassmanN.McNairK.SanchezS. E.SilvaG. G. Z.BolingL.et al (2014). A Highly Abundant Bacteriophage Discovered in the Unknown Sequences of Human Faecal Metagenomes. Nat. Commun.5, 4498. doi: 10.1038/ncomms5498
27
ErezZ.Steinberger-LevyI.ShamirM.DoronS.Stokar-AvihailA.PelegY.et al (2017). Communication Between Viruses Guides Lysis–Lysogeny Decisions. Nature541, 488–493. doi: 10.1038/nature21049
28
FeinerR.ArgovT.RabinovichL.SigalN.BorovokI.HerskovitsA. A. (2015). A New Perspective on Lysogeny: Prophages as Active Regulatory Switches of Bacteria. Nat. Rev. Microbiol.13, 641–650. doi: 10.1038/nrmicro3527
29
ForsterS. C.KumarN.AnonyeB. O.AlmeidaA.VicianiE.StaresM. D.et al (2019). A Human Gut Bacterial Genome and Culture Collection for Improved Metagenomic Analyses. Nat. Biotechnol.37, 186–192. doi: 10.1038/s41587-018-0009-7
30
GaliezC.SiebertM.EnaultF.VincentJ.SödingJ. (2017). WIsH: Who is the Host? Predicting Prokaryotic Hosts From Metagenomic Phage Contigs. Bioinformatics33, 3113–3114. doi: 10.1093/bioinformatics/btx383
31
GligorijevicV.RenfrewP. D.KosciolekT.LemanJ. K. (2021). Structure-Based Function Prediction Using Graph Convolutional Networks. Nat. Commun.12, 3168. doi: 10.1038/s41467-021-23303-9
32
GogokhiaL.BuhrkeK.BellR.HoffmanB.BrownD. G.Hanke-GogokhiaC.et al (2019). Expansion of Bacteriophages Is Linked to Aggravated Intestinal Inflammation and Colitis Article Expansion of Bacteriophages Is Linked to Aggravated Intestinal Inflammation and Colitis. Cell Host Microbe25, 285–299.e8. doi: 10.1016/j.chom.2019.01.008
33
GórskiA.KniotekM.Perkowska-PtasińskaA.MrózA.PrzerwaA.GorczycaW.et al (2006). Bacteriophages and Transplantation Tolerance. Transplant. Proc.38, 331–333. doi: 10.1016/j.transproceed.2005.12.073
34
GregoryA. C.ZablockiO.ZayedA. A.HowellA.BolducB.SullivanM. B. (2020). The Gut Virome Database Reveals Age-Dependent Patterns of Virome Diversity in the Human Gut. Cell Host Microbe28, 1–17. doi: 10.1016/j.chom.2020.08.003
35
GroseJ. H.CasjensS. R. (2014). Understanding the Enormous Diversity of Bacteriophages: The Tailed Phages That Infect the Bacterial Family Enterobacteriaceae. Virology468-470, 421–443. doi: 10.1016/j.virol.2014.08.024
36
HallA. R.ScanlanP. D.MorganA. D.BucklingA. (2011). Host-Parasite Coevolutionary Arms Races Give Way to Fluctuating Selection. Ecol. Lett.14, 635–642. doi: 10.1111/j.1461-0248.2011.01624.x
37
HamptonH. G.WatsonB. N. J.FineranP. C. (2020). The Arms Race Between Bacteria and Their Phage Foes. Nature577, 327–336. doi: 10.1038/s41586-019-1894-8
38
HarrisonE.BrockhurstM. A. (2017). Ecological and Evolutionary Benefits of Temperate Phage: What Does or Doesn’t Kill You Makes You Stronger. Bioessays39, 12. doi: 10.1002/bies.201700112
39
Heintz-BuschartA.WilmesP. (2018). Human Gut Microbiome: Function Matters. Trends Microbiol.26, 563–574. doi: 10.1016/j.tim.2017.11.002
40
HendrixR. W.SmithM. C. M.BurnsR. N.FordM. E.HatfullG. F. (1999). Evolutionary Relationships Among Diverse Bacteriophages and Prophages: All the World’s a Phage. Proc. Natl. Acad. Sci. U.S.A.96, 2192–2197. doi: 10.1073/pnas.96.5.2192
41
HobbsZ.AbedonS. T. (2016). Diversity of Phage Infection Types and Associated Terminology: The Problem With ‘Lytic or Lysogenic’. FEMS Microbiol. Lett.363, 1–8. doi: 10.1093/femsle/fnw047
42
HockenberryA. J.WilkeC. O. (2021). BACPHLIP: Predicting Bacteriophage Lifestyle From Conserved Protein Domains. PeerJ9, e11396. doi: 10.7717/peerj.11396
43
Hodyra-StefaniakK.MiernikiewiczP.DrapałaJ.DrabM.Jończyk-MatysiakE.LecionD.et al (2015). Mammalian Host-Versus-Phage Immune Response Determines Phage Fate In Vivo. Sci. Rep.5, 1–13. doi: 10.1038/srep14802
44
HolmfeldtK.SolonenkoN.ShahM.CorrierK.RiemannL.VerBerkmoesN. C.et al (2013). Twelve Previously Unknown Phage Genera are Ubiquitous in Global Oceans. Proc. Natl. Acad. Sci. U.S.A.110, 12798–12803. doi: 10.1073/pnas.1305956110
45
HoylesL.McCartneyA. L.NeveH.GibsonG. R.SandersonJ. D.HellerK. J.et al (2014). Characterization of Virus-Like Particles Associated With the Human Faecal and Caecal Microbiota. Res. Microbiol.165, 803–812. doi: 10.1016/j.resmic.2014.10.006
46
HsuB. B.GibsonT. E.YeliseyevV.LiuQ.LyonL.BryL.et al (2019). Dynamic Modulation of the Gut Microbiota and Metabolome by Bacteriophages in a Mouse Model. Cell Host Microbe25, 803–814.e5. doi: 10.1016/j.chom.2019.05.001
47
HuloC.de CastroE.MassonP.BougueleretL.BairochA.XenariosI.et al (2011). ViralZone: A Knowledge Resource to Understand Virus Diversity. Nucleic Acids Res.39, D576–D582. doi: 10.1093/nar/gkq901
48
JończykE.KłakM.MiędzybrodzkiR.GórskiA. (2011). The Influence of External Factors on Bacteriophages–Review. Folia Microbiol.56, 191–200. doi: 10.1007/s12223-011-0039-8
49
JumperJ.EvansR.PritzelA.GreenT.FigurnovM.RonnebergerO.et al (2021). Highly Accurate Protein Structure Prediction With AlphaFold. Nature596, 590–596. doi: 10.1038/s41586-021-03819-2
50
KauffmanK. M.HussainF. A.YangJ.ArevaloP.BrownJ. M.ChangW. K.et al (2018). A Major Lineage of Non-Tailed dsDNA Viruses as Unrecognized Killers of Marine Bacteria. Nature554, 118–122. doi: 10.1038/nature25474
51
KieftK.ZhouZ.AnantharamanK. (2020). VIBRANT: Automated Recovery, Annotation and Curation of Microbial Viruses, and Evaluation of Viral Community Function From Genomic Sequences. Microbiome8, 90. doi: 10.1186/s40168-020-00867-0
52
KimM. S.ParkE. J.RohS. W.BaeJ. W. (2011). Diversity and Abundance of Single-Stranded DNA Viruses in Human Feces. Appl. Environ. Microbiol.77, 8062–8070. doi: 10.1128/AEM.06331-11
53
KnowlesB.SilveiraC. B.BaileyB. A.BarottK.CantuV. A.Cobián-GüemesA. G.et al (2016). Lytic to Temperate Switching of Viral Communities. Nature531, 466–470. doi: 10.1038/nature17193
54
LabrieS. J.SamsonJ. E.MoineauS. (2010). Bacteriophage Resistance Mechanisms. Nat. Rev. Microbiol.8, 317–327. doi: 10.1038/nrmicro2315
55
LiangG.ZhaoC.ZhangH.MatteiL.Sherrill-MixS.BittingerK.et al (2020). The Stepwise Assembly of the Neonatal Virome Is Modulated by Breastfeeding. Nature581, 470–474. doi: 10.1038/s41586-020-2192-1
56
LimE. S.ZhouY.ZhaoG.BauerI. K.DroitL.NdaoI.M.et al (2015). Early Life Dynamics of the Human Gut Virome and Bacterial Microbiome in Infants. Nat. Med.21, 1228–1234. doi: 10.1038/nm.3950
57
LindellD.SullivanM. B.JohnsonZ. I.TolonenA. C.RohwerF.ChisholmS. W. (2004). Transfer of Photosynthesis Genes to and From Prochlorococcus Viruses. Proc. Natl. Acad. Sci. U.S.A.101, 11013–11018. doi: 10.1073/pnas.0401526101
58
LinD. M.LinH. C. (2019). A Theoretical Model of Temperate Phages as Mediators of Gut Microbiome Dysbiosis. F1000Res8, 997. doi: 10.12688/f1000research.18480.1
59
LiuH.WangJ.HeT.BeckerS.ZhangG.LiD.et al (2018). Butyrate: A Double-Edged Sword for Health? Adv. Nutr.9, 21–29. doi: 10.1093/advances/nmx009
60
LiuB.PopM. (2009). ARDB—Antibiotic Resistance Genes Database. Nucleic Acids Res.37, D443–D447. doi: 10.1093/nar/gkn656
61
ŁośJ. M.ŁośM.WęgrzynA. (2010). Hydrogen Peroxide-Mediated Induction of the Shiga Toxinconverting Lambdoid Prophage ST2-8624 in Escherichia Coli O157:H7. FEMS Immunol. Med. Microbiol.58, 322–329. doi: 10.1111/j.1574-695X.2009.00644.x
62
LuC.ZhangZ.CaiZ.ZhuZ.QiuY.WuA.et al (2021). Prokaryotic Virus Host Predictor: A Gaussian Model for Host Prediction of Prokaryotic Viruses in Metagenomics. BMC Biol.19, 5. doi: 10.1186/s12915-020-00938-6
63
MaY.PacanJ. C.WangQ.XuY.HuangX.KorenevskyA.et al (2008). Microencapsulation of Bacteriophage Felix O1 Into Chitosan-Alginate Microspheres for Oral Delivery. Appl. Environ. Microbiol.74, 4799–4805. doi: 10.1128/AEM.00246-08
64
MajewskaJ.BetaW.LecionD.Hodyra-StefaniakK.KłopotA.KaźmierczakZ.et al (2015). Oral Application of T4 Phage Induces Weak Antibody Production in the Gut and in the Blood. Viruses7, 4783–4799. doi: 10.3390/v7082845
65
ManriqueP.BolducB.WalkS. T.van der OostJ.de VosW. M.YoungM. J. (2016). Healthy Human Gut Phageome. Proc. Natl. Acad. Sci.113, 10400–10405. doi: 10.1073/pnas.1601060113
66
MartínR.MiquelS.UlmerJ.LangellaP.Bermúdez-HumaránL. G. (2014). Gut Ecosystem: How Microbes Help Us. Benef. Microbes5, 219–233. doi: 10.3920/BM2013.0057
67
McNairK.BaileyB. A.EdwardsR. A. (2012). PHACTS, A Computational Approach to Classifying the Lifestyle of Phages. Bioinformatics28, 614–618. doi: 10.1093/bioinformatics/bts014
68
MiedzybrodzkiR.Switala-JelenK.FortunaW.Weber-DabrowskaB.PrzerwaA.Lusiak-SzelachowskaM.et al (2008). Bacteriophage Preparation Inhibition of Reactive Oxygen Species Generation by Endotoxin-Stimulated Polymorphonuclear Leukocytes. Virus Res.131, 233–242. doi: 10.1016/j.virusres.2007.09.013
69
MiernikiewiczP.KłopotA.SoluchR.SzkutaP.KęskaW.Hodyra-StefaniakK.et al (2016). T4 Phage Tail Adhesin Gp12 Counteracts LPS-Induced Inflammation In Vivo. Front. Microbiol.7, 1–8. doi: 10.3389/fmicb.2016.01112
70
MirzaeiM. K.MauriceC. F. (2017). Ménage À Trois in the Human Gut: Interactions Between Host, Bacteria and Phages. Nat. Rev. Microbiol.15, 397–408. doi: 10.1038/nrmicro.2017.30
71
MulakA.BonazB. (2015). Brain-Gut-Microbiota Axis in Parkinson’s Disease. World J. Gastroenterol.21, 10609. doi: 10.3748/wjg.v21.i37.10609
72
NayfachS.Páez-EspinoD.CallL.LowS. J.SberroH.IvanovaN. N.et al (2021). Metagenomic Compendium of 189,680 DNA Viruses From the Human Gut Microbiome. Nat. Microbiol.6, 960–970. doi: 10.1038/s41564-021-00928-6
73
NguyenS.BakerK.PadmanB. S.PatwaR.DunstanR. A.WestonT. A.et al (2017). Bacteriophage Transcytosis Provides a Mechanism To Cross Epithelial Cell Layers. MBio8, 6. doi: 10.1128/mBio.01874-17
74
NormanJ. M.HandleyS. A.BaldridgeM. T.DroitL.LiuC. Y.KellerB. C.et al (2015). Disease-Specific Alterations in the Enteric Virome in Inflammatory Bowel Disease. Cell160, 447–460. doi: 10.1016/j.cell.2015.01.002
75
OfirG.SorekR. (2018). Contemporary Phage Biology: From Classic Models to New Insights. Cell172, 1260–1270. doi: 10.1016/j.cell.2017.10.045
76
OliveiraH.SampaioM.MeloL. D. R.DiasO.PopeW. H.HatfullG. F.et al (2019). Staphylococci Phages Display Vast Genomic Diversity and Evolutionary Relationships. BMC Genomics20, 357. doi: 10.1186/s12864-019-5647-8
77
OppenheimA. B.KobilerO.StavansJ.CourtD. L.AdhyaS. (2005). Switches in Bacteriophage Lambda Development. Annu. Rev. Genet.39, 409–429. doi: 10.1146/annurev.genet.39.073003.113656
78
ParkH.LaffinM. R.JovelJ.MillanB.HyunJ. E.HotteN.et al (2019). The Success of Fecal Microbial Transplantation in Clostridium Difficile Infection Correlates With Bacteriophage Relative Abundance in the Donor: A Retrospective Cohort Study. Gut Microbes10, 676–687. doi: 10.1080/19490976.2019.1586037
79
PenadésJ. R.ChenJ.Quiles-PuchaltN.CarpenaN.NovickR. P. (2015). Bacteriophage-Mediated Spread of Bacterial Virulence Genes. Curr. Opin. Microbiol.23, 171–178. doi: 10.1016/j.mib.2014.11.019
80
PopeW. H.BowmanC. A.RussellD. A.Jacobs-SeraD.AsaiD. J.CresawnS. G.et al (2015). Whole Genome Comparison of a Large Collection of Mycobacteriophages Reveals a Continuum of Phage Genetic Diversity. Elife4, e06416. doi : 10.7554/eLife.06416
81
RamisettyB. C. M.SudhakariP. A. (2019). Bacterial ‘Grounded’ Prophages: Hotspots for Genetic Renovation and Innovation. Front. Genet.10, 65. doi: 10.3389/fgene.2019.00065
82
RastallR. A. (2004). Bacteria in the Gut: Friends and Foes and How to Alter the Balance. J. Nutr.134, 2022S–2026S. doi: 10.1093/jn/134.8.2022S
83
ReyesA.HaynesM.HansonN.AnglyF. E.HeathA. C.RohwerF.et al (2010). Viruses in the Fecal Microbiota of Monozygotic Twins and Their Mothers. Nature466, 334–338. doi: 10.1038/nature09199
84
RohwerF.SegallA.StewardG.SeguritanV.BreitbartM.WolvenF.et al (2000). The Complete Genomic Sequence of the Marine Phage Roseophage SIO1 Shares Homology With Nonmarine Phages. Limnol. Oceanogr.45, 408–418. doi: 10.4319/lo.2000.45.2.0408
85
RokutanK.KawaharaT.KuwanoY.TominagaK.NishidaK.Teshima-KondoS. (2008). Nox Enzymes and Oxidative Stress in the Immunopathology of the Gastrointestinal Tract. Semin. Immunopathol.30, 315–327. doi: 10.1007/s00281-008-0124-5
86
RouxS.SolonenkoN. E.DangV. T.PoulosB. T.SchwenkS. M.GoldsmithD. B.et al (2016). Towards Quantitative Viromics for Both Double-Stranded and Single-Stranded DNA Viruses. PeerJ2016, 1–17. doi: 10.7717/peerj.2777
87
SamsonJ. E.MagadánA. H.SabriM.MoineauS. (2013). Revenge of the Phages: Defeating Bacterial Defences. Nat. Rev. Microbiol.11, 675–687. doi: 10.1038/nrmicro3096
88
SardelliL.PerottoniS.TunesiM.BoeriL.FuscoF.PetriniP.et al (2021). Technological Tools and Strategies for Culturing Human Gut Microbiota in Engineered In Vitro Models. Biotechnol. Bioeng.118, 2886–2905. doi: 10.1002/bit.27816
89
SaussetR.PetitM. A.Gaboriau-RouthiauV.De PaepeM. (2020). New Insights Into Intestinal Phages. Mucosal Immunol.13, 205–215. doi: 10.1038/s41385-019-0250-5
90
SchmererM.MolineuxI. J.BullJ. J. (2014). Synergy as a Rationale for Phage Therapy Using Phage Cocktails. PeerJ2, e590. doi: 10.7717/peerj.590
91
SegainJ.-P.Raingeard de la BlétièreD.BourreilleA.LerayV.GervoisN.RosalesC.et al (2000). Butyrate Inhibits Inflammatory Responses Through Nfκb Inhibition: Implications for Crohn’s Disease. Gut47, 397–403. doi: 10.1136/gut.47.3.397
92
ShafferM.BortonM. A.McGivernB. B.ZayedA. A.Leanti La RosaS.SoldenL. M.et al (2020). DRAM for Distilling Microbial Metabolism to Automate the Curation of Microbiome Function. Nucleic Acids Res.48, 8883–8900. doi: 10.1093/nar/gkaa621
93
SherwoodL.WilleyJ.WoolvertonC. J. (2013). Prescott’s Microbiology. London, UK: McGraw-Hill Education.
94
ShiN.LiN.DuanX.NiuH. (2017). Interaction Between the Gut Microbiome and Mucosal Immune System. Mil. Med. Res.4, 14. doi: 10.1186/s40779-017-0122-9
95
ShkoporovA. N.KhokhlovaE. V.FitzgeraldC. B.StockdaleS. R.DraperL. A.RossP.et al (2018a). Φcrass001 Represents the Most Abundant Bacteriophage Family in the Human Gut and Infects Bacteroides Intestinalis. Nat. Commun.9, 4781. doi: 10.1038/s41467-018-07225-7
96
ShkoporovA. N.RyanF. J.DraperL. A.FordeA.StockdaleS. R.DalyK. M.et al (2018b). Reproducible Protocols for Metagenomic Analysis of Human Faecal Phageomes. Microbiome6, 68. doi: 10.1186/s40168-018-0446-z
97
ShkoporovA. N.ClooneyA. G.SuttonT. D. S.RyanF. J.DalyK. M.NolanJ. A.et al (2019). The Human Gut Virome Is Highly Diverse, Stable, and Individual Specific. Cell Host Microbe26, 527–541.e5. doi: 10.1016/j.chom.2019.09.009
98
ShkoporovA. N.HillC. (2019). Bacteriophages of the Human Gut: The ‘Known Unknown’ of the Microbiome. Cell Host Microbe25, 195–209. doi: 10.1016/j.chom.2019.01.017
99
SimmondsP.AdamsM. J.BenkóM.BreitbartM.BristerJ. R.CarstensE. B.et al (2017). Consensus Statement: Virus Taxonomy in the Age of Metagenomics. Nat. Rev. Microbiol.15, 161–168. doi: 10.1038/nrmicro.2016.177
100
SinhaA.MauriceC. F. (2019). Bacteriophages: Uncharacterized and Dynamic Regulators of the Immune System. Mediators Inflamm.2019, 3730519. doi: 10.1155/2019/3730519
101
SokolH.SeksikP.FuretJ. P.FirmesseO.Nion-LarmurierI.BeaugerieL.et al (2009). Low Counts of Faecalibacterium Prausnitzii in Colitis Microbiota. Inflamm. Bowel Dis.15, 1183–1189. doi: 10.1002/ibd.20903
102
StecherB.HardtW. D. (2008). The Role of Microbiota in Infectious Disease. Trends Microbiol.16, 107–114. doi: 10.1016/j.tim.2007.12.008
103
Stuart ChapinF.IIIMatsonP. A.VitousekP. (2011). Principles of Terrestrial Ecosystem Ecology. New York, USA: Springer Science & Business Media.
104
SuttleC. A. (2005). Viruses in the Sea. Nature437, 356–361. doi: 10.1038/nature04160
105
TakeuchiO.AkiraS. (2010). Pattern Recognition Receptors and Inflammation. Cell140, 805–820. doi: 10.1016/j.cell.2010.01.022
106
TetzG.BrownS. M.HaoY.TetzV. (2018). Parkinson’s Disease and Bacteriophages as Its Overlooked Contributors. Sci. Rep.8, 10812. doi: 10.1038/s41598-018-29173-4
107
TetzG.TetzV. (2018). Bacteriophages as New Human Viral Pathogens. Microorganisms6, 54. doi: 10.3390/microorganisms6020054
108
ThingstadT. F. (2000). Elements of a Theory for the Mechanisms Controlling Abundance, Diversity, and Biogeochemical Role of Lytic Bacterial Viruses in Aquatic Systems. Limnol. Oceanogr.45, 1320–1328. doi: 10.4319/lo.2000.45.6.1320
109
ThurberR. V.HaynesM.BreitbartM.WegleyL.RohwerF. (2009). Laboratory Procedures to Generate Viral Metagenomes. Nat. Protoc.4, 470–483. doi: 10.1038/nprot.2009.10
110
ThursbyE.JugeN. (2017). Introduction to the Human Gut Microbiota. Biochem. J.474, 1823–1836. doi: 10.1042/BCJ20160510
111
TunyasuvunakoolK.AdlerJ.WuZ.GreenT.ZielinskiM.ŽídekA.et al (2021). Highly Accurate Protein Structure Prediction for the Human Proteome. Nature596, 590–596. doi: 10.1038/s41586-021-03828-1
112
Van BelleghemJ. D.DąbrowskaK.VaneechoutteM.BarrJ. J.BollykyP. L. (2019). Interactions Between Bacteriophage, Bacteria, and the Mammalian Immune System. Viruses11, 10. doi: 10.3390/v11010010
113
van HouteS.BucklingA.WestraE. R. (2016). Evolutionary Ecology of Prokaryotic Immune Mechanisms. Microbiol. Mol. Biol. Rev.80, 745–763. doi: 10.1128/MMBR.00011-16
114
VerthéK.PossemiersS.BoonN.VaneechoutteM.VerstraeteW. (2004). Stability and Activity of an Enterobacter Aerogenes-Specific Bacteriophage Under Simulated Gastro-Intestinal Conditions. Appl. Microbiol. Biotechnol.65, 465–472. doi: 10.1007/s00253-004-1585-7
115
WagnerP. L.WaldorM. K. (2002). Bacteriophage Control of Bacterial Virulence. Infect. Immun.70, 3985–3993. doi: 10.1128/IAI.70.8.3985-3993.2002
116
WangX.KimY.MaQ.HongS. H.PokusaevaK.SturinoJ. M.et al (2010). Cryptic Prophages Help Bacteria Cope With Adverse Environments. Nat. Commun.1, 147. doi: 10.1038/ncomms1146
117
WangW.RenJ.TangK.DartE.Ignacio-EspinozaJ. C.FuhrmanJ. A.et al (2020). A Network-Based Integrated Framework for Predicting Virus-Prokaryote Interactions. NAR Genom. Bioinform.2, lqaa044. doi: 10.1093/nargab/lqaa044
118
WangX.QuinnP. J. (2010). Lipopolysaccharide: Biosynthetic Pathway and Structure Modification. Prog. Lipid Res.49, 97–107. doi: 10.1016/j.plipres.2009.06.002
119
WellsJ. M.BrummerR. J.DerrienM.MacDonaldT. T.TroostF.CaniP. D.et al (2017). Homeostasis of the Gut Barrier and Potential Biomarkers. Am. J. Physiol. Gastrointest. Liver Physiol.312, G171–G193. doi: 10.1152/ajpgi.00048.2015
120
WuS.FangZ.TanJ.LiM.WangC.GuoQ.et al (2020). Distinguish Virulent and Temperate Phage-Derived Sequences in Metavirome Data With a Deep Learning Approach. GigaScience10, 9. doi: 10.1101/2020.12.25.424404
121
YangM.XiaQ.SenD.ZhangZ.QinF.ZhaoY. (2021). Genomic Characterization and Distribution Pattern of a Novel Marine OM43 Phage. Front. Microbiol.12, 651326. doi: 10.3389/fmicb.2021.651326
122
ZhaiZ.ZhangZ.ZhaoG.LiuX.QinF.ZhaoY. (2021). Genomic Characterization of Two Novel RCA Phages Reveals New Insights Into the Diversity and Evolution of Marine Viruses. Microbiol. Spectr.9, e0123921. doi: 10.1128/Spectrum.01239-21
123
ZhangY.-Z.ChenY.-M.WangW.QinX.-C.HolmesE. C. (2019). Expanding the RNA Virosphere by Unbiased Metagenomics. Annu. Rev. Virol.6, 119–139. doi: 10.1146/annurev-virology-092818-015851
124
ZhangT.BreitbartM.LeeW. H.RunJ. Q.WeiC. L.SohS. W. L.et al (2006). RNA Viral Community in Human Feces: Prevalence of Plant Pathogenic Viruses. PloS Biol.4, e3. doi: 10.1371/journal.pbio.0040003
125
ZhangF.ZhouF.GanR.RenC.JiaY.YuL.et al (2020). PHISDetector: A Tool to Detect Diverse in Silico Phage-Host Interaction Signals for Virome Studies. bioRxiv. Preprint. doi: 10.1101/661074
126
ZhangY.ZhangZ.ZhangH.ZhaoY.ZhangZ.XiaoJ. (2020). PADS Arsenal: A Database of Prokaryotic Defense Systems Related Genes. Nucleic Acids Res.48, D590–D598. doi: 10.1093/nar/gkz916
127
ZhangR.MirditaM.KarinE. L.NorroyC.GaliezC.SödingJ. (2021). SpacePHARER: Sensitive Identification of Phages From CRISPR Spacers in Prokaryotic Hosts. Bioinformatics37, 3364–3366. doi: 10.1093/bioinformatics/btab222
128
ZhouZ.TranP. Q.BreisterA. M.LiuY.KieftK.CowleyE. S.et al (2020). METABOLIC: High-Throughput Profiling of Microbial Genomes for Functional Traits, Biogeochemistry, and Community-Scale Metabolic Networks. Res. Sq. Preprint. doi: 10.21203/rs.3.rs-113327/v1
129
ZhouC. E.SmithJ.LamM.ZemlaA. (2007). MvirDB—A Microbial Database of Protein Toxins, Virulence Factors and Antibiotic Resistance Genes for Bio-Defence Applications. Nucleic Acids35, D391–D394. doi: 10.1093/nar/gkl791
130
ZuoT.WongS. H.LamK.LuiR.CheungK.TangW.et al (2018). Bacteriophage Transfer During Faecal Microbiota Transplantation in Clostridium Difficile Infection is Associated With Treatment Outcome. Gut67, 634–643. doi: 10.1136/gutjnl-2017-313952
131
ZuppiM.TozzoliR.ChianiP.QuirosP.Martinez-VelazquezA.MichelacciV.et al (2020). Investigation on the Evolution of Shiga Toxin-Converting Phages Based on Whole Genome Sequencing. Front. Microbiol.11, 1–14. doi: 10.3389/fmicb.2020.01472
Summary
Keywords
gut microbiome, bacteriophage, prophage, microbial ecology, metagenomics
Citation
Zuppi M, Hendrickson HL, O’Sullivan JM and Vatanen T (2022) Phages in the Gut Ecosystem. Front. Cell. Infect. Microbiol. 11:822562. doi: 10.3389/fcimb.2021.822562
Received
25 November 2021
Accepted
10 December 2021
Published
04 January 2022
Volume
11 - 2021
Edited by
Yuan Gao, Beijing Institute of Genomics (CAS), China
Reviewed by
Felix Broecker, Idorsia Pharmaceuticals Ltd, Switzerland; Yunxue Guo, South China Sea Institute of Oceanology (CAS), China
Updates
Copyright
© 2022 Zuppi, Hendrickson, O’Sullivan and Vatanen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tommi Vatanen, t.vatanen@auckland.ac.nz
This article was submitted to Microbiome in Health and Disease, a section of the journal Frontiers in Cellular and Infection Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.