Genome signature-based dissection of human gut metagenomes to extract subliminal viral sequences

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

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

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

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

Inference of the Life Cycle of Environmental Phages from Genomic Signature Distances to Their Hosts

The environmental impact of uncultured phages is shaped by their preferred life cycle (lytic or lysogenic). However, our ability to predict it is very limited. We aimed to discriminate between lytic and lysogenic phages by comparing the similarity of their genomic signatures to those of their hosts, reflecting their co-evolution. We tested two approaches: (1) similarities of tetramer relative frequencies, (2) alignment-free comparisons based on exact k = 14 oligonucleotide matches. First, we explored 5126 reference bacterial host strains and 284 associated phages and found an approximate threshold for distinguishing lysogenic and lytic phages using both oligonucleotide-based methods. The analysis of 6482 plasmids revealed the potential for horizontal gene transfer between different host genera and, in some cases, distant bacterial taxa. Subsequently, we experimentally analyzed combinations of 138 Klebsiella pneumoniae strains and their 41 phages and found that the phages with the largest number of interactions with these strains in the laboratory had the shortest genomic distances to K. pneumoniae. We then applied our methods to 24 single-cells from a hot spring biofilm containing 41 uncultured phage-host pairs, and the results were compatible with the lysogenic life cycle of phages detected in this environment. In conclusion, oligonucleotide-based genome analysis methods can be used for predictions of (1) life cycles of environmental phages, (2) phages with the broadest host range in culture collections, and (3) potential horizontal gene transfer by plasmids.

ViralCC retrieves complete viral genomes and virus-host pairs from metagenomic Hi-C data

The introduction of high-throughput chromosome conformation capture (Hi-C) into metagenomics enables reconstructing high-quality metagenome-assembled genomes (MAGs) from microbial communities. Despite recent advances in recovering eukaryotic, bacterial, and archaeal genomes using Hi-C contact maps, few of Hi-C-based methods are designed to retrieve viral genomes. Here we introduce ViralCC, a publicly available tool to recover complete viral genomes and detect virus-host pairs using Hi-C data. Compared to other Hi-C-based methods, ViralCC leverages the virus-host proximity structure as a complementary information source for the Hi-C interactions. Using mock and real metagenomic Hi-C datasets from several different microbial ecosystems, including the human gut, cow fecal, and wastewater, we demonstrate that ViralCC outperforms existing Hi-C-based binning methods as well as state-of-the-art tools specifically dedicated to metagenomic viral binning. ViralCC can also reveal the taxonomic structure of viruses and virus-host pairs in microbial communities. When applied to a real wastewater metagenomic Hi-C dataset, ViralCC constructs a phage-host network, which is further validated using CRISPR spacer analyses. ViralCC is an open-source pipeline available at https://github.com/dyxstat/ViralCC .

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

Virsearcher: Identifying Bacteriophages from Metagenomes by Combining Convolutional Neural Network and Gene Information

Metagenome sequencing provides an unprecedented opportunity for the discovery of unknown microbes and viruses. A large number of phages and prokaryotes are mixed together in metagenomes. To study the influence of phages on human bodies and environments, it is of great significance to isolate phages from metagenomes. However, it is difficult to identify novel phages because of the diversity of their sequences and the frequent presence of short contigs in metagenomes. Here, virSearcher is developed to identify phages from metagenomes by combining the convolutional neural network (CNN) and the gene information of input sequences. Firstly, an input sequence is encoded in accordance with the different functions of its coding and the non-coding regions and then is converted into word embedding code through a word embedding layer before a convolutional layer. Meanwhile, the hit ratio of the virus genes is combined with the output of the CNN to further improve the performance of the network. The genes used by virSearcher consist of complete and incomplete genes. Experiments on several metagenomes have showed that, compared with others, virSearcher can significantly improve the performance for the identification of short sequences, while maintaining the performance for long ones. The source code of virSearcher is freely available from http://github.com/DrJackson18/virSearcher.

DeephageTP: a convolutional neural network framework for identifying phage-specific proteins from metagenomic sequencing data

Bacteriophages (phages) are the most abundant and diverse biological entity on Earth. Due to the lack of universal gene markers and database representatives, there about 50-90% of genes of phages are unable to assign functions. This makes it a challenge to identify phage genomes and annotate functions of phage genes efficiently by homology search on a large scale, especially for newly phages. Portal (portal protein), TerL (large terminase subunit protein), and TerS (small terminase subunit protein) are three specific proteins of Caudovirales phage. Here, we developed a CNN (convolutional neural network)-based framework, DeephageTP, to identify the three specific proteins from metagenomic data. The framework takes one-hot encoding data of original protein sequences as the input and automatically extracts predictive features in the process of modeling. To overcome the false positive problem, a cutoff-loss-value strategy is introduced based on the distributions of the loss values of protein sequences within the same category. The proposed model with a set of cutoff-loss-values demonstrates high performance in terms of Precision in identifying TerL and Portal sequences (94% and 90%, respectively) from the mimic metagenomic dataset. Finally, we tested the efficacy of the framework using three real metagenomic datasets, and the results shown that compared to the conventional alignment-based methods, our proposed framework had a particular advantage in identifying the novel phage-specific protein sequences of portal and TerL with remote homology to their counterparts in the training datasets. In summary, our study for the first time develops a CNN-based framework for identifying the phage-specific protein sequences with high complexity and low conservation, and this framework will help us find novel phages in metagenomic sequencing data. The DeephageTP is available at https://github.com/chuym726/DeephageTP.

DeepHost: phage host prediction with convolutional neural network

Next-generation sequencing expands the known phage genomes rapidly. Unlike culture-based methods, the hosts of phages discovered from next-generation sequencing data remain uncharacterized. The high diversity of the phage genomes makes the host assignment task challenging. To solve the issue, we proposed a phage host prediction tool-DeepHost. To encode the phage genomes into matrices, we design a genome encoding method that applied various spaced $k$-mer pairs to tolerate sequence variations, including insertion, deletions, and mutations. DeepHost applies a convolutional neural network to predict host taxonomies. DeepHost achieves the prediction accuracy of 96.05% at the genus level (72 taxonomies) and 90.78% at the species level (118 taxonomies), which outperforms the existing phage host prediction tools by 10.16-30.48% and achieves comparable results to BLAST. For the genomes without hits in BLAST, DeepHost obtains the accuracy of 38.00% at the genus level and 26.47% at the species level, making it suitable for genomes of less homologous sequences with the existing datasets. DeepHost is alignment-free, and it is faster than BLAST, especially for large datasets. DeepHost is available at https://github.com/deepomicslab/DeepHost.

Ecogenomics of Groundwater Phages Suggests Niche Differentiation Linked to Specific Environmental Tolerance

Viruses are ubiquitous microbiome components, shaping ecosystems via strain-specific predation, horizontal gene transfer and redistribution of nutrients through host lysis. Viral impacts are important in groundwater ecosystems, where microbes drive many nutrient fluxes and metabolic processes; however, little is known about the diversity of viruses in these environments. We analyzed four groundwater plasmidomes (the entire plasmid content of an environment) and identified 200 viral sequences, which clustered into 41 genus-level viral clusters (approximately equivalent to viral genera) including 9 known and 32 putative new genera. We used publicly available bacterial whole-genome sequences (WGS) and WGS from 261 bacterial isolates from this groundwater environment to identify potential viral hosts. We linked 76 of the 200 viral sequences to a range of bacterial phyla, the majority associated with Proteobacteria, followed by Firmicutes, Bacteroidetes, and Actinobacteria. The publicly available WGS enabled mapping bacterial hosts to several viral sequences. The WGS of groundwater isolates increased the depth of host prediction by allowing host identification at the strain level. The latter included 4 viruses that were almost entirely (>99% query coverage, >99% identity) identified as integrated in the genomes of Pseudomonas, Acidovorax, and Castellaniella strains, resulting in high-confidence host assignments. Lastly, 21 of these viruses carried putative auxiliary metabolite genes for metal and antibiotic resistance, which might drive their infection cycles and/or provide selective advantage to infected hosts. Exploring the groundwater virome provides a necessary foundation for integration of viruses into ecosystem models where they are key players in microbial adaption to environmental stress. IMPORTANCE To our knowledge, this is the first study to identify the bacteriophage distribution in a groundwater ecosystem shedding light on their prevalence and distribution across metal-contaminated and background sites. Our study is uniquely based on selective sequencing of solely the extrachromosomal elements of a microbiome followed by analysis for viral signatures, thus establishing a more focused approach for phage identifications. Using this method, we detected several novel phage genera along with those previously established. Our approach of using the whole-genome sequences of hundreds of bacterial isolates from the same site enabled us to make host assignments with high confidence, several at strain levels. Certain phage genes suggest that they provide an environment-specific selective advantage to their bacterial hosts. Our study lays the foundation for future research on directed phage isolations using specific bacterial host strains to further characterize groundwater phages, their life cycles, and their effects on groundwater microbiome and biogeochemistry.

The Human Gut Phageome: Origins and Roles in the Human Gut Microbiome

The investigation of the microbial populations of the human body, known as the microbiome, has led to a revolutionary field of science, and understanding of its impacts on human development and health. The majority of microbiome research to date has focussed on bacteria and other kingdoms of life, such as fungi. Trailing behind these is the interrogation of the gut viruses, specifically the phageome. Bacteriophages, viruses that infect bacterial hosts, are known to dictate the dynamics and diversity of bacterial populations in a number of ecosystems. However, the phageome of the human gut, while of apparent importance, remains an area of many unknowns. In this paper we discuss the role of bacteriophages within the human gut microbiome. We examine the methods used to study bacteriophage populations, how this evolved over time and what we now understand about the phageome. We review the phageome development in infancy, and factors that may influence phage populations in adult life. The role and action of the phageome is then discussed at both a biological-level, and in the broader context of human health and disease.

Infection with Bacteroides Phage BV01 Alters the Host Transcriptome and Bile Acid Metabolism in a Common Human Gut Microbe

Gut-associated phages are hypothesized to alter the abundance and activity of their bacterial hosts, contributing to human health and disease. Although temperate phages constitute a significant fraction of the gut virome, the effects of lysogenic infection are underexplored. We report that the temperate phage, Bacteroides phage BV01, broadly alters its host's transcriptome, the prominent human gut symbiont Bacteroides vulgatus. This alteration occurs through phage-induced repression of a tryptophan-rich sensory protein (TspO) and represses bile acid deconjugation. Because microbially modified bile acids are important signals for the mammalian host, this is a mechanism by which a phage may influence mammalian phenotypes. Furthermore, BV01 and its relatives in the proposed phage family Salyersviridae are ubiquitous in human gut metagenomes, infecting a broad range of Bacteroides hosts. These results demonstrate the complexity of phage-bacteria-mammal relationships and emphasize a need to better understand the role of temperate phages in the gut microbiome.

Swine virome on rural backyard farms in Mexico: communities with different abundances of animal viruses and phages

Domestic swine have been introduced by humans into a wide diversity of environments and have been bred in different production systems. This has resulted in an increased risk for the occurrence and spread of diseases. Although viromes of swine in intensive farms have been described, little is known about the virus communities in backyard production systems around the world. The aim of this study was to describe the viral diversity of 23 healthy domestic swine maintained in rural backyards in Morelos, Mexico, through collection and analysis of nasal and rectal samples. Next-generation sequencing was used to identify viruses that are present in swine. Through homology search and bioinformatic analysis of reads and their assemblies, we found that rural backyard swine have a high degree of viral diversity, different from those reported in intensive production systems or under experimental conditions. There was a higher frequency of bacteriophages and lower diversity of animal viruses than reported previously. In addition, sapoviruses, bocaparvoviruses, and mamastroviruses that had not been reported previously in our country were identified. These findings were correlated with the health status of animals, their social interactions, and the breeding/rearing environment (which differed from intensive systems), providing baseline information about viral communities in backyard swine.

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

Bacteriophages have important roles in the ecology of the human gut microbiome but are under-represented in reference databases. To address this problem, we assembled the Metagenomic Gut Virus catalogue that comprises 189,680 viral genomes from 11,810 publicly available human stool metagenomes. Over 75% of genomes represent double-stranded DNA phages that infect members of the Bacteroidia and Clostridia classes. Based on sequence clustering we identified 54,118 candidate viral species, 92% of which were not found in existing databases. The Metagenomic Gut Virus catalogue improves detection of viruses in stool metagenomes and accounts for nearly 40% of CRISPR spacers found in human gut Bacteria and Archaea. We also produced a catalogue of 459,375 viral protein clusters to explore the functional potential of the gut virome. This revealed tens of thousands of diversity-generating retroelements, which use error-prone reverse transcription to mutate target genes and may be involved in the molecular arms race between phages and their bacterial hosts.

Characterization of the gut DNA and RNA Viromes in a Cohort of Chinese Residents and Visiting Pakistanis

Trillions of viruses inhabit the gastrointestinal tract. Some of them have been well-studied on their roles in infection and human health, but the majority remains unsurveyed. It has been established that the composition of the gut virome is highly variable based on the changes of diet, physical state, and environmental factors. However, the effect of host genetic factors, for example ethnic origin, on the gut virome is rarely investigated. Here, we characterized and compared the gut virome in a cohort of local Chinese residents and visiting Pakistani individuals, each group containing twenty-four healthy adults and six children. Using metagenomic shotgun sequencing and assembly of fecal samples, a huge number of viral operational taxonomic units (vOTUs) were identified for profiling the DNA and RNA viromes. National background contributed a primary variation to individuals' gut virome. Compared with the Chinese adults, the Pakistan adults showed higher macrodiversity and different compositional and functional structures in their DNA virome and lower diversity and altered composition in their RNA virome. The virome variations of Pakistan children were not only inherited from that of the adults but also tended to share similar characteristics with the Chinese cohort. We also analyzed and compared the bacterial microbiome between two cohorts and further revealed numerous connections between viruses and bacterial host. Statistically, the gut DNA and RNA viromes were covariant to some extent (P < 0.001), and they both correlated the holistic bacterial composition and vice versa. This study provides an overview of the gut viral community in Chinese and visiting Pakistanis and proposes a considerable role of ethnic origin in shaping the virome.

Genome Characterization of a Novel Wastewater Bacteroides fragilis Bacteriophage (vB_BfrS_23) and its Host GB124

Bacteroides spp. are part of the human intestinal microbiota but can under some circumstances become clinical pathogens. Phages are a potentially valuable therapeutic treatment option for many pathogens, but phage therapy for pathogenic Bacteroides spp. including Bacteroides fragilis is currently limited to three genome-sequenced phages. Here we describe the isolation from sewage wastewater and genome of a lytic phage, vB_BfrS_23, that infects and kills B. fragilis strain GB124. Transmission electron microscopy identified this phage as a member of the Siphoviridae family. The phage is stable when held at temperatures of 4 and 60°C for 1 h. It has a very narrow host range, only infecting one host from a panel of B. fragilis strains (n = 8). Whole-genome sequence analyses of vB_BfrS_23 determined it is double-stranded DNA phage and is circularly permuted, with a genome of 48,011 bp. The genome encodes 73 putative open reading frames. We also sequenced the host bacterium, B. fragilis GB124 (5.1 Mb), which has two plasmids of 43,923 and 4,138 bp. Although this phage is host specific, its isolation together with the detailed characterization of the host B. fragilis GB124 featured in this study represent a useful starting point from which to facilitate the future development of highly specific therapeutic agents. Furthermore, the phage could be a novel tool in determining water (and water reuse) treatment efficacy, and for identifying human fecal transmission pathways within contaminated environmental waters and foodstuffs.

Probing the Mobilome: Discoveries in the Dynamic Microbiome

There has been an explosion of metagenomic data representing human, animal, and environmental microbiomes. This provides an unprecedented opportunity for comparative and longitudinal studies of many functional aspects of the microbiome that go beyond taxonomic classification, such as profiling genetic determinants of antimicrobial resistance, interactions with the host, potentially clinically relevant functions, and the role of mobile genetic elements (MGEs). One of the most important but least studied of these aspects are the MGEs, collectively referred to as the 'mobilome'. Here we elaborate on the benefits and limitations of using different metagenomic protocols, discuss the relative merits of various sequencing technologies, and highlight relevant bioinformatics tools and pipelines to predict the presence of MGEs and their microbial hosts.

Rumen Virus Populations: Technological Advances Enhancing Current Understanding

The rumen contains a multi-kingdom, commensal microbiome, including protozoa, bacteria, archaea, fungi and viruses, which enables ruminant herbivores to ferment and utilize plant feedstuffs that would be otherwise indigestible. Within the rumen, virus populations are diverse and highly abundant, often out-numbering the microbial populations that they both predate on and co-exist with. To date the research effort devoted to understanding rumen-associated viral populations has been considerably less than that given to the other microbial populations, yet their contribution to maintaining microbial population balance, intra-ruminal microbial lysis, fiber breakdown, nutrient cycling and genetic transfer may be highly significant. This review follows the technological advances which have contributed to our current understanding of rumen viruses and drawing on knowledge from other environmental and animal-associated microbiomes, describes the known and potential roles and impacts viruses have on rumen function and speculates on the future directions of rumen viral research.

Bacteriophage Isolation and Characterization: Phages of Escherichia coli

Here we introduce methods for the detection, enumeration, and isolation of bacteriophages from Escherichia coli. In bacteria, horizontal gene transfer may be mediated by virulent and temperate phages. Strict virulent phages, able to propagate in a suitable strain following the lytic pathway, can be isolated directly from different natural environments. In temperate phages, the lytic cycle must be activated, and phages are detected after their induction. In both cases, detection is based on the production of visible plaques in a confluent lawn of the host strain using a double agar layer method. Further purification and characterization are achieved by density gradients, electron microscopy studies, and genomic analysis. This straightforward methodology can be applied to the detection, enumeration, and isolation of bacteriophages from any bacterial species, using the appropriate host strain, media, and culture conditions.

Genomic and Ecogenomic Characterization of Proteus mirabilis Bacteriophages

Proteus mirabilis often complicates the care of catheterized patients through the formation of crystalline biofilms which block urine flow. Bacteriophage therapy has been highlighted as a promising approach to control this problem, but relatively few phages infecting P. mirabilis have been characterized. Here we characterize five phages capable of infecting P. mirabilis, including those shown to reduce biofilm formation, and provide insights regarding the wider ecological and evolutionary relationships of these phages. Transmission electron microscopy (TEM) imaging of phages vB_PmiP_RS1pmA, vB_PmiP_RS1pmB, vB_PmiP_RS3pmA, and vB_PmiP_RS8pmA showed that all share morphologies characteristic of the Podoviridae family. The genome sequences of vB_PmiP_RS1pmA, vB_PmiP_RS1pmB, and vB_PmiP_RS3pmA showed these are species of the same phage differing only by point mutations, and are closely related to vB_PmiP_RS8pmA. Podophages characterized in this study were also found to share similarity in genome architecture and composition to other previously described P. mirabilis podophages (PM16 and PM75). In contrast, vB_PimP_RS51pmB showed morphology characteristic of the Myoviridae family, with no notable similarity to other phage genomes examined. Ecogenomic profiling of all phages revealed no association with human urinary tract viromes, but sequences similar to vB_PimP_RS51pmB were found within human gut, and human oral microbiomes. Investigation of wider host-phage evolutionary relationships through tetranucleotide profiling of phage genomes and bacterial chromosomes, indicated vB_PimP_RS51pmB has a relatively recent association with Morganella morganii and other non-Proteus members of the Morganellaceae family. Subsequent host range assays confirmed vB_PimP_RS51pmB can infect M. morganii.

Beyond cells - The virome in the human holobiont

Viromics, or viral metagenomics, is a relatively new and burgeoning field of research that studies the complete collection of viruses forming part of the microbiota in any given niche. It has strong foundations rooted in over a century of discoveries in the field of virology and recent advances in molecular biology and sequencing technologies. Historically, most studies have deconstructed the concept of viruses into a simplified perception of viral agents as mere pathogens, which demerits the scope of large-scale viromic analyses. Viruses are, in fact, much more than regular parasites. They are by far the most dynamic and abundant entity and the greatest killers on the planet, as well as the most effective geo-transforming genetic engineers and resource recyclers, acting on all life strata in any habitat. Yet, most of this uncanny viral world remains vastly unexplored to date, greatly hindered by the bewildering complexity inherent to such studies and the methodological and conceptual limitations. Viromic studies are just starting to address some of these issues but they still lag behind microbial metagenomics. In recent years, however, higher-throughput analysis and resequencing have rekindled interest in a field that is just starting to show its true potential. In this review, we take a look at the scientific and technological developments that led to the advent of viral and bacterial metagenomics with a particular, but not exclusive, focus on human viromics from an ecological perspective. We also address some of the most relevant challenges that current viral studies face and ponder on the future directions of the field.

Mining, analyzing, and integrating viral signals from metagenomic data

BACKGROUND

Viruses are important components of microbial communities modulating community structure and function; however, only a couple of tools are currently available for phage identification and analysis from metagenomic sequencing data. Here we employed the random forest algorithm to develop VirMiner, a web-based phage contig prediction tool especially sensitive for high-abundances phage contigs, trained and validated by paired metagenomic and phagenomic sequencing data from the human gut flora.

RESULTS

VirMiner achieved 41.06% ± 17.51% sensitivity and 81.91% ± 4.04% specificity in the prediction of phage contigs. In particular, for the high-abundance phage contigs, VirMiner outperformed other tools (VirFinder and VirSorter) with much higher sensitivity (65.23% ± 16.94%) than VirFinder (34.63% ± 17.96%) and VirSorter (18.75% ± 15.23%) at almost the same specificity. Moreover, VirMiner provides the most comprehensive phage analysis pipeline which is comprised of metagenomic raw reads processing, functional annotation, phage contig identification, and phage-host relationship prediction (CRISPR-spacer recognition) and supports two-group comparison when the input (metagenomic sequence data) includes different conditions (e.g., case and control). Application of VirMiner to an independent cohort of human gut metagenomes obtained from individuals treated with antibiotics revealed that 122 KEGG orthology and 118 Pfam groups had significantly differential abundance in the pre-treatment samples compared to samples at the end of antibiotic administration, including clustered regularly interspaced short palindromic repeats (CRISPR), multidrug resistance, and protein transport. The VirMiner webserver is available at http://sbb.hku.hk/VirMiner/ .

CONCLUSIONS

We developed a comprehensive tool for phage prediction and analysis for metagenomic samples. Compared to VirSorter and VirFinder-the most widely used tools-VirMiner is able to capture more high-abundance phage contigs which could play key roles in infecting bacteria and modulating microbial community dynamics.

TRIAL REGISTRATION

The European Union Clinical Trials Register, EudraCT Number: 2013-003378-28 . Registered on 9 April 2014.

Tracking the Rules of Transmission and Introgression with Networks

Understanding how an animal organism and its gut microbes form an integrated biological organization, known as a holobiont, is becoming a central issue in biological studies. Such an organization inevitably involves a complex web of transmission processes that occur on different scales in time and space, across microbes and hosts. Network-based models are introduced in this chapter to tackle aspects of this complexity and to better take into account vertical and horizontal dimensions of transmission. Two types of network-based models are presented, sequence similarity networks and bipartite graphs. One interest of these networks is that they can consider a rich diversity of important players in microbial evolution that are usually excluded from evolutionary studies, like plasmids and viruses. These methods bring forward the notion of "gene externalization," which is defined as the presence of redundant copies of prokaryotic genes on mobile genetic elements (MGEs), and therefore emphasizes a related although distinct process from lateral gene transfer between microbial cells. This chapter introduces guidelines to the construction of these networks, reviews their analysis, and illustrates their possible biological interpretations and uses. The application to human gut microbiomes shows that sequences present in a higher diversity of MGEs have both biased functions and a broader microbial and human host range. These results suggest that an "externalized gut metagenome" is partly common to humans and benefits the gut microbial community. We conclude that testing relationships between microbial genes, microbes, and their animal hosts, using network-based methods, could help to unravel additional mechanisms of transmission in holobionts.

Phage particles harboring antibiotic resistance genes in fresh-cut vegetables and agricultural soil

Bacteriophages are ubiquitously distributed prokaryotic viruses that are more abundant than bacteria. As a consequence of their life cycle, phages can kidnap part of their host's genetic material, including antibiotic resistance genes (ARGs), which released phage particles transfer in a process called transduction. The spread of ARGs among pathogenic bacteria currently constitutes a serious global health problem. In this study, fresh vegetables (lettuce, spinach and cucumber), and cropland soil were screened by qPCR for ten ARGs (bla_TEM, bla_CTX-M-1 group, bla_CTX-M-9 group, bla_OXA-48, bla_VIM, mecA, sul1, qnrA, qnrS and armA) in their viral DNA fraction. The presence of ARGs in the phage DNA was analyzed before and after propagation experiments in an Escherichia coli host strain to evaluate the ability of the phage particles to infect a host. ARGs were found in the phage DNA fraction of all matrices, although with heterogeneous values. ARG prevalence was significantly higher in lettuce and soil, and the most common overall were β-lactamases. After propagation experiments, an increase in ARG densities in phage particles was observed in samples of all four matrices, confirming that part of the isolated phage particles were infectious. This study reveals the abundance of free, replicative ARG-containing phage particles in vegetable matrices and cropland soil. The particles are proposed as vehicles for resistance transfer in these environments, where they can persist for a long time, with the possibility of generating new resistant bacterial strains. Ingestion of these mobile genetic elements may also favor the emergence of new resistances, a risk not previously considered.

Global Distribution of Human-Associated Fecal Genetic Markers in Reference Samples from Six Continents

Numerous bacterial genetic markers are available for the molecular detection of human sources of fecal pollution in environmental waters. However, widespread application is hindered by a lack of knowledge regarding geographical stability, limiting implementation to a small number of well-characterized regions. This study investigates the geographic distribution of five human-associated genetic markers (HF183/BFDrev, HF183/BacR287, BacHum-UCD, BacH, and Lachno2) in municipal wastewaters (raw and treated) from 29 urban and rural wastewater treatment plants (750-4 400 000 population equivalents) from 13 countries spanning six continents. In addition, genetic markers were tested against 280 human and nonhuman fecal samples from domesticated, agricultural and wild animal sources. Findings revealed that all genetic markers are present in consistently high concentrations in raw (median log₁₀ 7.2-8.0 marker equivalents (ME) 100 mL^-1) and biologically treated wastewater samples (median log₁₀ 4.6-6.0 ME 100 mL^-1) regardless of location and population. The false positive rates of the various markers in nonhuman fecal samples ranged from 5% to 47%. Results suggest that several genetic markers have considerable potential for measuring human-associated contamination in polluted environmental waters. This will be helpful in water quality monitoring, pollution modeling and health risk assessment (as demonstrated by QMRAcatch) to guide target-oriented water safety management across the globe.

A human gut phage catalog correlates the gut phageome with type 2 diabetes

BACKGROUND

Substantial efforts have been made to link the gut bacterial community to many complex human diseases. Nevertheless, the gut phages are often neglected.

RESULTS

In this study, we used multiple bioinformatic methods to catalog gut phages from whole-community metagenomic sequencing data of fecal samples collected from both type II diabetes (T2D) patients (n = 71) and normal Chinese adults (n = 74). The definition of phage operational taxonomic units (pOTUs) and identification of large phage scaffolds (n = 2567, ≥ 10 k) revealed a comprehensive human gut phageome with a substantial number of novel sequences encoding genes that were unrelated to those in known phages. Interestingly, we observed a significant increase in the number of gut phages in the T2D group and, in particular, identified 7 pOTUs specific to T2D. This finding was further validated in an independent dataset of 116 T2D and 109 control samples. Co-occurrence/exclusion analysis of the bacterial genera and pOTUs identified a complex core interaction between bacteria and phages in the human gut ecosystem, suggesting that the significant alterations of the gut phageome cannot be explained simply by co-variation with the altered bacterial hosts.

CONCLUSIONS

Alterations in the gut bacterial community have been linked to the chronic disease T2D, but the role of gut phages therein is not well understood. This is the first study to identify a T2D-specific gut phageome, indicating the existence of other mechanisms that might govern the gut phageome in T2D patients. These findings suggest the importance of the phageome in T2D risk, which warrants further investigation.

Resolution of habitat-associated ecogenomic signatures in bacteriophage genomes and application to microbial source tracking

Just as the expansion in genome sequencing has revealed and permitted the exploitation of phylogenetic signals embedded in bacterial genomes, the application of metagenomics has begun to provide similar insights at the ecosystem level for microbial communities. However, little is known regarding this aspect of bacteriophage associated with microbial ecosystems, and if phage encode discernible habitat-associated signals diagnostic of underlying microbiomes. Here we demonstrate that individual phage can encode clear habitat-related 'ecogenomic signatures', based on relative representation of phage-encoded gene homologues in metagenomic data sets. Furthermore, we show the ecogenomic signature encoded by the gut-associated ɸB124-14 can be used to segregate metagenomes according to environmental origin, and distinguish 'contaminated' environmental metagenomes (subject to simulated in silico human faecal pollution) from uncontaminated data sets. This indicates phage-encoded ecological signals likely possess sufficient discriminatory power for use in biotechnological applications, such as development of microbial source tracking tools for monitoring water quality.

The human gut virome: form and function

Advances in next-generation sequencing technologies and the application of metagenomic approaches have fuelled an exponential increase in our understanding of the human gut microbiome. These approaches are now also illuminating features of the diverse and abundant collection of viruses (termed the virome) subsisting with the microbial ecosystems residing within the human holobiont. Here, we focus on the current and emerging knowledge of the human gut virome, in particular on viruses infecting bacteria (bacteriophage or phage), which are a dominant component of this viral community. We summarise current insights regarding the form and function of this 'human gut phageome' and highlight promising avenues for future research. In doing so, we discuss the potential for phage to drive ecological functioning and evolutionary change within this important microbial ecosystem, their contribution to modulation of host-microbiome interactions and stability of the community as a whole, as well as the potential role of the phageome in human health and disease. We also consider the emerging concepts of a 'core healthy gut phageome' and the putative existence of 'viral enterotypes' and 'viral dysbiosis'.

Review article: the human intestinal virome in health and disease

BACKGROUND

The human virome consists of animal-cell viruses causing transient infections, bacteriophage (phage) predators of bacteria and archaea, endogenous retroviruses and viruses causing persistent and latent infections. High-throughput, inexpensive, sensitive sequencing methods and metagenomics now make it possible to study the contribution dsDNA, ssDNA and RNA virus-like particles make to the human virome, and in particular the intestinal virome.

AIM

To review and evaluate the pioneering studies that have attempted to characterise the human virome and generated an increased interest in understanding how the intestinal virome might contribute to maintaining health, and the pathogenesis of chronic diseases.

METHODS

Relevant virome-related articles were selected for review following extensive language- and date-unrestricted, electronic searches of the literature.

RESULTS

The human intestinal virome is personalised and stable, and dominated by phages. It develops soon after birth in parallel with prokaryotic communities of the microbiota, becoming established during the first few years of life. By infecting specific populations of bacteria, phages can alter microbiota structure by killing host cells or altering their phenotype, enabling phages to contribute to maintaining intestinal homeostasis or microbial imbalance (dysbiosis), and the development of chronic infectious and autoimmune diseases including HIV infection and Crohn's disease, respectively.

CONCLUSIONS

Our understanding of the intestinal virome is fragmented and requires standardised methods for virus isolation and sequencing to provide a more complete picture of the virome, which is key to explaining the basis of virome-disease associations, and how enteric viruses can contribute to disease aetiologies and be rationalised as targets for interventions.

The Human Gut Phage Community and Its Implications for Health and Disease

In this review, we assess our current understanding of the role of bacteriophages infecting the human gut bacterial community in health and disease. In general, bacteriophages contribute to the structure of their microbial communities by driving host and viral diversification, bacterial evolution, and by expanding the functional diversity of ecosystems. Gut bacteriophages are an ensemble of unique and shared phages in individuals, which encompass temperate phages found predominately as prophage in gut bacteria (prophage reservoir) and lytic phages. In healthy individuals, only a small fraction of the prophage reservoir is activated and found as extracellular phages. Phage community dysbiosis is characterized by a shift in the activated prophage community or an increase of lytic phages, and has been correlated with disease, suggesting that a proper balance between lysis and lysogeny is needed to maintain health. Consequently, the concept of microbial dysbiosis might be extended to the phage component of the microbiome as well. Understanding the dynamics and mechanisms to restore balance after dysbiosis is an active area of research. The use of phage transplants to re-establish health suggests that phages can be used as disease treatment. Such advances represent milestones in our understanding of gut phages in human health and should fuel research on their role in health and disease.

Lysogeny in nature: mechanisms, impact and ecology of temperate phages

Viruses that infect bacteria (phages) can influence bacterial community dynamics, bacterial genome evolution and ecosystem biogeochemistry. These influences differ depending on whether phages establish lytic, chronic or lysogenic infections. Although the first two produce virion progeny, with lytic infections resulting in cell destruction, phages undergoing lysogenic infections replicate with cells without producing virions. The impacts of lysogeny are numerous and well-studied at the cellular level, but ecosystem-level consequences remain underexplored compared to those of lytic infections. Here, we review lysogeny from molecular mechanisms to ecological patterns to emerging approaches of investigation. Our goal is to highlight both its diversity and importance in complex communities. Altogether, using a combined viral ecology toolkit that is applied across broad model systems and environments will help us understand more of the diverse lifestyles and ecological impacts of lysogens in nature.

Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses

Ocean microbes drive biogeochemical cycling on a global scale. However, this cycling is constrained by viruses that affect community composition, metabolic activity, and evolutionary trajectories. Owing to challenges with the sampling and cultivation of viruses, genome-level viral diversity remains poorly described and grossly understudied, with less than 1% of observed surface-ocean viruses known. Here we assemble complete genomes and large genomic fragments from both surface- and deep-ocean viruses sampled during the Tara Oceans and Malaspina research expeditions, and analyse the resulting 'global ocean virome' dataset to present a global map of abundant, double-stranded DNA viruses complete with genomic and ecological contexts. A total of 15,222 epipelagic and mesopelagic viral populations were identified, comprising 867 viral clusters (defined as approximately genus-level groups). This roughly triples the number of known ocean viral populations and doubles the number of candidate bacterial and archaeal virus genera, providing a near-complete sampling of epipelagic communities at both the population and viral-cluster level. We found that 38 of the 867 viral clusters were locally or globally abundant, together accounting for nearly half of the viral populations in any global ocean virome sample. While two-thirds of these clusters represent newly described viruses lacking any cultivated representative, most could be computationally linked to dominant, ecologically relevant microbial hosts. Moreover, we identified 243 viral-encoded auxiliary metabolic genes, of which only 95 were previously known. Deeper analyses of four of these auxiliary metabolic genes (dsrC, soxYZ, P-II (also known as glnB) and amoC) revealed that abundant viruses may directly manipulate sulfur and nitrogen cycling throughout the epipelagic ocean. This viral catalog and functional analyses provide a necessary foundation for the meaningful integration of viruses into ecosystem models where they act as key players in nutrient cycling and trophic networks.

Diversity in a Polymicrobial Community Revealed by Analysis of Viromes, Endolysins and CRISPR Spacers

The polymicrobial biofilm communities in Mushroom and Octopus Spring in Yellowstone National Park (YNP) are well characterized, yet little is known about the phage populations. Dominant species, Synechococcus sp. JA-2-3B'a(2–13), Synechococcus sp. JA-3-3Ab, Chloroflexus sp. Y-400-fl, and Roseiflexus sp. RS-1, contain multiple CRISPR-Cas arrays, suggesting complex interactions with phage predators. To analyze phage populations from Octopus Spring biofilms, we sequenced a viral enriched fraction. To assemble and analyze phage metagenomic data, we developed a custom module, VIRITAS, implemented within the MetAMOS framework. This module bins contigs into groups based on tetranucleotide frequencies and CRISPR spacer-protospacer matching and ORF calling. Using this pipeline we were able to assemble phage sequences into contigs and bin them into three clusters that corroborated with their potential host range. The virome contained 52,348 predicted ORFs; some were clearly phage-like; 9319 ORFs had a recognizable Pfam domain while the rest were hypothetical. Of the recognized domains with CRISPR spacer matches, was the phage endolysin used by lytic phage to disrupt cells. Analysis of the endolysins present in the thermophilic cyanophage contigs revealed a subset of characterized endolysins as well as a Glyco_hydro_108 (PF05838) domain not previously associated with sequenced cyanophages. A search for CRISPR spacer matches to all identified phage endolysins demonstrated that a majority of endolysin domains were targets. This strategy provides a general way to link host and phage as endolysins are known to be widely distributed in bacteriophage. Endolysins can also provide information about host cell wall composition and have the additional potential to be used as targets for novel therapeutics.

Phages rarely encode antibiotic resistance genes: a cautionary tale for virome analyses

Antibiotic resistance genes (ARGs) are pervasive in gut microbiota, but it remains unclear how often ARGs are transferred, particularly to pathogens. Traditionally, ARG spread is attributed to horizontal transfer mediated either by DNA transformation, bacterial conjugation or generalized transduction. However, recent viral metagenome (virome) analyses suggest that ARGs are frequently carried by phages, which is inconsistent with the traditional view that phage genomes rarely encode ARGs. Here we used exploratory and conservative bioinformatic strategies found in the literature to detect ARGs in phage genomes, and experimentally assessed a subset of ARG predicted using exploratory thresholds. ARG abundances in 1181 phage genomes were vastly overestimated using exploratory thresholds (421 predicted vs 2 known), due to low similarities and matches to protein unrelated to antibiotic resistance. Consistent with this, four ARGs predicted using exploratory thresholds were experimentally evaluated and failed to confer antibiotic resistance in Escherichia coli. Reanalysis of available human- or mouse-associated viromes for ARGs and their genomic context suggested that bona fide ARG attributed to phages in viromes were previously overestimated. These findings provide guidance for documentation of ARG in viromes, and reassert that ARGs are rarely encoded in phages.

Freshwater Metaviromics and Bacteriophages: A Current Assessment of the State of the Art in Relation to Bioinformatic Challenges

Advances in bioinformatics and sequencing technologies have allowed for the analysis of complex microbial communities at an unprecedented rate. While much focus is often placed on the cellular members of these communities, viruses play a pivotal role, particularly bacteria-infecting viruses (bacteriophages); phages mediate global biogeochemical processes and drive microbial evolution through bacterial grazing and horizontal gene transfer. Despite their importance and ubiquity in nature, very little is known about the diversity and structure of viral communities. Though the need for culture-based methods for viral identification has been somewhat circumvented through metagenomic techniques, the analysis of metaviromic data is marred with many unique issues. In this review, we examine the current bioinformatic approaches for metavirome analyses and the inherent challenges facing the field as illustrated by the ongoing efforts in the exploration of freshwater phage populations.

Repertoire of human gut microbes

In 1675, Antoni Van Leeuwenhoeck was the first to observe several forms using an optical microscope that he named "animalcules", realizing later that these were microorganisms. The first classification of living organisms proposed by Ehrenberg in 1833 was based on what we could visualize. The failure of this kind of classification arises from viral culture, which preceded direct observations that were finally achieved during the 20th century by electron microscopy. The number of prokaryotic species is estimated at approximately 10 million, although only 1800 were known in 1980, and 14,000 to date, thanks to the advent of 16S rRNA amplification and sequencing. This highlights our inability to access the entire diversity. Indeed, a large number of bacteria are only, known as Operational Taxonomic Units (OTUs) and detected as a result of metagenomics studies, revealing an unexplored world known as the "dark matter". Recently, the rebirth of bacterial culture through the example of culturomics has dramatically increased the human gut repertoire as well as the 18SrRNA sequencing allowed to largely extend the repertoire of Eukaryotes. Finally, filtration and co-culture on free-living protists associated with high-throughput culture elucidated a part of the megavirome. While the majority of studies currently performed on the human gut microbiota focus on bacterial diversity, it appears that several other prokaryotes (including archaea) and eukaryotic populations also inhabit this ecosystem; their detection depending exclusively on the tools used. Rational and comprehensive establishment of this ecosystem will allow the understanding of human health associated with gut microbiota and the potential to change this.

Long-term microbiota and virome in a Zürich patient after fecal transplantation against Clostridium difficile infection

Fecal microbiota transplantation (FMT) is an emerging therapeutic option for Clostridium difficile infections that are refractory to conventional treatment. FMT introduces fecal microbes into the patient's intestine that prevent the recurrence of C. difficile, leading to rapid expansion of bacteria characteristic of healthy microbiota. However, the long-term effects of FMT remain largely unknown. The C. difficile patient described in this paper revealed protracted microbiota adaptation processes from 6 to 42 months post-FMT. Ultimately, bacterial communities were donor similar, suggesting sustainable stool engraftment. Since little is known about the consequences of transmitted viruses during C. difficile infection, we also interrogated virome changes. Our approach allowed identification of about 10 phage types per sample that represented larger viral communities, and phages were found to be equally abundant in the cured patient and donor. The healthy microbiota appears to be characterized by low phage abundance. Although viruses were likely transferred, the patient established a virome distinct from the donor. Surprisingly, the patient had sequences of algal giant viruses (chloroviruses) that have not previously been reported for the human gut. Chloroviruses have not been associated with intestinal disease, but their presence in the oropharynx may influence cognitive abilities. The findings suggest that the virome is an important indicator of health or disease. A better understanding of the role of viruses in the gut ecosystem may uncover novel microbiota-modulating therapeutic strategies.

The information science of microbial ecology

A revolution is unfolding in microbial ecology where petabytes of 'multi-omics' data are produced using next generation sequencing and mass spectrometry platforms. This cornucopia of biological information has enormous potential to reveal the hidden metabolic powers of microbial communities in natural and engineered ecosystems. However, to realize this potential, the development of new technologies and interpretative frameworks grounded in ecological design principles are needed to overcome computational and analytical bottlenecks. Here we explore the relationship between microbial ecology and information science in the era of cloud-based computation. We consider microorganisms as individual information processing units implementing a distributed metabolic algorithm and describe developments in ecoinformatics and ubiquitous computing with the potential to eliminate bottlenecks and empower knowledge creation and translation.

Spatial disturbances in altered mucosal and luminal gut viromes of diet-induced obese mice

Gut microbial biogeography is a key feature of host-microbe relationships. In gut viral ecology, biogeography and responses to dietary intervention remain poorly understood. Here, we conducted a metagenomic study to determine the composition of the mucosal and luminal viromes of the gut and to evaluate the impact of a Western diet on gut viral ecology. We found that mucosal and luminal viral assemblages comprised predominantly temperate phages. The mucosal virome significantly differed from the luminal virome in low-fat diet-fed lean mice, where spatial variation correlated with bacterial microbiota from the mucosa and lumen. The mucosal and luminal viromes of high-fat, high-sucrose 'Western' diet-fed obese mice were significantly enriched with temperate phages of the Caudovirales order. Interestingly, this community alteration occurred to a greater extent in the mucosa than lumen, leading to loss of spatial differences; however, these changes recovered after switching to a low-fat diet. Temperate phages enriched in the Western diet-induced obese mice were associated with the Bacilli, Negativicutes and Bacteroidia classes and temperate phages from the Bacteroidia class particularly encoded stress and niche-specific functions advantageous to bacterial host adaptation. This study illustrates a biogeographic view of the gut virome and phage-bacterial host connections under the diet-induced microbial dysbiosis.

Computational approaches to predict bacteriophage-host relationships

Metagenomics has changed the face of virus discovery by enabling the accurate identification of viral genome sequences without requiring isolation of the viruses. As a result, metagenomic virus discovery leaves the first and most fundamental question about any novel virus unanswered: What host does the virus infect? The diversity of the global virosphere and the volumes of data obtained in metagenomic sequencing projects demand computational tools for virus–host prediction. We focus on bacteriophages (phages, viruses that infect bacteria), the most abundant and diverse group of viruses found in environmental metagenomes. By analyzing 820 phages with annotated hosts, we review and assess the predictive power of in silico phage–host signals. Sequence homology approaches are the most effective at identifying known phage–host pairs. Compositional and abundance-based methods contain significant signal for phage–host classification, providing opportunities for analyzing the unknowns in viral metagenomes. Together, these computational approaches further our knowledge of the interactions between phages and their hosts. Importantly, we find that all reviewed signals significantly link phages to their hosts, illustrating how current knowledge and insights about the interaction mechanisms and ecology of coevolving phages and bacteria can be exploited to predict phage–host relationships, with potential relevance for medical and industrial applications.

New viruses infecting bacteria are increasingly being discovered in many environments through sequence-based explorations. To understand their role in microbial ecosystems, computational tools are indispensable to prioritize and guide experimental efforts. This review assesses and discusses a range of bioinformatic approaches to predict bacteriophage–host relationships when all that is known is their genome sequence.

Metagenome Analysis Exploiting High-Throughput Chromosome Conformation Capture (3C) Data

Microbial communities are complex and constitute important parts of our environment. Genomic analysis of these populations is a dynamic research area but remains limited by the difficulty in assembling full genomes of individual species. Recently, a new method for metagenome assembly/analysis based on chromosome conformation capture has emerged (meta3C). This approach quantifies the collisions experienced by DNA molecules to identify those sharing the same cellular compartments, allowing the characterization of genomes present within complex mixes of species. The exploitation of these chromosome 3D signatures holds promising perspectives for genome sequencing of discrete species in complex populations. It also has the potential to assign correctly extra-chromosomal elements, such as plasmids, mobile elements and phages, to their host cells.

The human gut virome: a multifaceted majority

Here, we outline our current understanding of the human gut virome, in particular the phage component of this ecosystem, highlighting progress, and challenges in viral discovery in this arena. We reveal how developments in high-throughput sequencing technologies and associated data analysis methodologies are helping to illuminate this abundant 'biological dark matter.' Current evidence suggests that the human gut virome is a highly individual but temporally stable collective, dominated by phages exhibiting a temperate lifestyle. This viral community also appears to encode a surprisingly rich functional repertoire that confers a range of attributes to their bacterial hosts, ranging from bacterial virulence and pathogenesis to maintaining host-microbiome stability and community resilience. Despite the significant advances in our understanding of the gut virome in recent years, it is clear that we remain in a period of discovery and revelation, as new methods and technologies begin to provide deeper understanding of the inherent ecological characteristics of this viral ecosystem. As our understanding increases, the nature of the multi-partite interactions occurring between host and microbiome will become clearer, helping us to more rationally define the concepts and principles that will underpin approaches to using human gut virome components for medical or biotechnological applications.

Metagenomic Sequencing Indicates That the Oropharyngeal Phageome of Individuals With Schizophrenia Differs From That of Controls

Mucosal sites such as the oropharynx contain a wide range of microorganisms, collectively designated as the microbiome. The microbiome can affect behavior through a number of neurobiological and immunological mechanisms. Most previous studies have focused on the bacterial components of the microbiome. However, the microbiome also includes viruses such as bacteriophages, which are viruses that infect bacteria and alter their metabolism and replication. We employed metagenomic analysis to characterize bacteriophage genomes in the oral pharynx of 41 individuals with schizophrenia and 33 control individuals without a psychiatric disorder. This analysis was performed by the generation of more than 100,000,000 sequence reads from each sample and the mapping of these reads to databases. We identified 79 distinct bacteriophage sequences in the oropharyngeal samples. Of these, one bacteriophage genome, Lactobacillus phage phiadh, was found to be significantly different in individuals with schizophrenia (P < .00037, q < 0.03 adjusted for multiple comparisons). The differential levels of Lactobacillus phage phiadh remained significant when controlling for age, gender, race, socioeconomic status, or cigarette smoking (P < .006). Within the group of individuals with schizophrenia, the level of Lactobacillus phage phiadh correlated with the prevalence of immunological disorders as well as with the administration of valproate, which has been shown in animal models to alter the microbiome. The bacteriophage composition of the oropharynx in individuals with schizophrenia differs from that of controls. The biological consequences of this difference and the potential effects of altering bacteriophage levels through therapeutic interventions are worthy of further investigation.

Methods for virus classification and the challenge of incorporating metagenomic sequence data

The division of viruses into orders, families, genera and species provides a classification framework that seeks to organize and make sense of the diversity of viruses infecting animals, plants and bacteria. Classifications are based on similarities in genome structure and organization, the presence of homologous genes and sequence motifs and at lower levels such as species, host range, nucleotide and antigenic relatedness and epidemiology. Classification below the level of family must also be consistent with phylogeny and virus evolutionary histories. Recently developed methods such as PASC, DEMaRC and NVR offer alternative strategies for genus and species assignments that are based purely on degrees of divergence between genome sequences. They offer the possibility of automating classification of the vast number of novel virus sequences being generated by next-generation metagenomic sequencing. However, distance-based methods struggle to deal with the complex evolutionary history of virus genomes that are shuffled by recombination and reassortment, and where taxonomic lineages evolve at different rates. In biological terms, classifications based on sequence distances alone are also arbitrary whereas the current system of virus taxonomy is of utility precisely because it is primarily based upon phenotypic characteristics. However, a separate system is clearly needed by which virus variants that lack biological information might be incorporated into the ICTV classification even if based solely on sequence relationships to existing taxa. For these, simplified taxonomic proposals and naming conventions represent a practical way to expand the existing virus classification and catalogue our rapidly increasing knowledge of virus diversity.

Identifying and analyzing bacteriophages in human fecal samples: what could we discover?

The human gut is a complex ecosystem, densely populated with microbes including enormous amounts of phages. Metagenomic studies indicate a great diversity of bacteriophages, and because of the variety of gut bacterial species, the human or animal gut is probably a perfect ecological niche for phages that can infect and propagate in their bacterial communities. In addition, some phages have the capacity to mobilize genes, as demonstrated by the enormous fraction of phage particles in feces that contain bacterial DNA. All these facts indicate that, through predation and horizontal gene transfer, bacteriophages play a key role in shaping the size, structure and function of intestinal microbiomes, although our understanding of their effects on gut bacterial populations is only just beginning.

Replicating phages in the epidermal mucosa of the eel (Anguilla anguilla)

In this work, we used the eel (Anguilla anguilla) as an animal model to test the hypothesis of Barr et al. (2013a,b) about the putative role of the epidermal mucosa as a phage enrichment layer. To this end, we analyzed the microbial content of the skin mucus of wild and farmed eels by using a metagenomic approach. We found a great abundance of replicating phage genomes (concatemers) in all the samples. They were assembled in four complete genomes of three Myovirus and one Podovirus. We also found evidences that ΦKZ and Podovirus phages could be part of the resident microbiota associated to the eel mucosal surface and persist on them over the time. Moreover, the viral abundance estimated by epiflorescent counts and by metagenomic recruitment from eel mucosa was higher than that of the surrounding water. Taken together, our results support the hypothesis that claims a possible role of phages in the animal mucus as agents controlling bacterial populations, including pathogenic species, providing a kind of innate immunity.

Tetranucleotide usage highlights genomic heterogeneity among mycobacteriophages

Background

The genomic sequences of mycobacteriophages, phages infecting mycobacterial hosts, are diverse and mosaic. Mycobacteriophages often share little nucleotide similarity, but most of them have been grouped into lettered clusters and further into subclusters. Traditionally, mycobacteriophage genomes are analyzed based on sequence alignment or knowledge of gene content. However, these approaches are computationally expensive and can be ineffective for significantly diverged sequences. As an alternative to alignment-based genome analysis, we evaluated tetranucleotide usage in mycobacteriophage genomes. These methods make it easier to characterize features of the mycobacteriophage population at many scales.

Description

We computed tetranucleotide usage deviation (TUD), the ratio of observed counts of 4-mers in a genome to the expected count under a null model. TUD values are comparable between members of a phage subcluster and distinct between subclusters. With few exceptions, neighbor joining phylogenetic trees and hierarchical clustering dendrograms constructed using TUD values place phages in a monophyletic clade with members of the same subcluster. Regions in a genome with exceptional TUD values can point to interesting features of genomic architecture. Finally, we found that subcluster B3 mycobacteriophages contain significantly overrepresented 4-mers and 6-mers that are atypical of phage genomes.

Conclusions

Statistics based on tetranucleotide usage support established clustering of mycobacteriophages and can uncover interesting relationships within and between sequenced phage genomes. These methods are efficient to compute and do not require sequence alignment or knowledge of gene content. The code to download mycobacteriophage genome sequences and reproduce our analysis is freely available at https://github.com/bsiranosian/tango_final.

Direct sequencing of human gut virome fractions obtained by flow cytometry

The sequence assembly of the human gut virome encounters several difficulties. A high proportion of human and bacterial matches is detected in purified viral samples. Viral DNA extraction results in a low DNA concentration, which does not reach the minimal limit required for sequencing library preparation. Therefore, the viromes are usually enriched by whole genome amplification (WGA), which is, however, prone to the development of chimeras and amplification bias. In addition, as there is a very wide diversity of gut viral species, very extensive sequencing efforts must be made for the assembling of whole viral genomes. We present an approach to improve human gut virome assembly by employing a more precise preparation of a viral sample before sequencing. Particles present in a virome previously filtered through 0.2 μm pores were further divided into groups in accordance with their size and DNA content by fluorescence activated cell sorting (FACS). One selected viral fraction was sequenced excluding the WGA step, so that unbiased sequences with high reliability were obtained. The DNA extracted from the 314 viral particles of the selected fraction was assembled into 34 contigs longer than 1,000 bp. This represents an increase to the number of assembled long contigs per sequenced Gb in comparison with other studies where non-fractioned viromes are sequenced. Seven of these contigs contained open reading frames (ORFs) with explicit matches to proteins related to bacteriophages. The remaining contigs also possessed uncharacterized ORFs with bacteriophage-related domains. When the particles that are present in the filtered viromes are sorted into smaller groups by FACS, large pieces of viral genomes can be recovered easily. This approach has several advantages over the conventional sequencing of non-fractioned viromes: non-viral contamination is reduced and the sequencing efforts required for viral assembly are minimized.