High-throughput, single-microbe genomics with strain resolution, applied to a human gut microbiome.

Characterizing complex microbial communities with single-cell resolution has been a long-standing goal of microbiology. We present Microbe-seq, a high-throughput method that yields the genomes of individual microbes from complex microbial communities. We encapsulate individual microbes in droplets with microfluidics and liberate their DNA, which we then amplify, tag with droplet-specific barcodes, and sequence. We explore the human gut microbiome, sequencing more than 20,000 microbial single-amplified genomes (SAGs) from a single human donor and coassembling genomes of almost 100 bacterial species, including several with multiple subspecies strains. We use these genomes to probe microbial interactions, reconstructing the horizontal gene transfer (HGT) network and observing HGT between 92 species pairs; we also identify a significant in vivo host-phage association between crAssphage and one strain of Bacteroides vulgatus. Microbe-seq contributes high-throughput culture-free capabilities to investigate genomic blueprints of complex microbial communities with single-microbe resolution.

Long-term persistence of crAss-like phage crAss001 is associated with phase variation in Bacteroides intestinalis.

BACKGROUND: The crAss-like phages are ubiquitous and highly abundant members of the human gut virome that infect commensal bacteria of the order Bacteroidales. Although incapable of lysogeny, these viruses demonstrate long-term persistence in the human gut microbiome, dominating the virome in some individuals. RESULTS: Here we show that rapid phase variation of alternate capsular polysaccharides in Bacteroides intestinalis cultures plays an important role in a dynamic equilibrium between phage sensitivity and resistance, allowing phage and bacteria to multiply in parallel. The data also suggests the role of a concomitant phage persistence mechanism associated with delayed lysis of infected cells, similar to carrier state infection. From an ecological and evolutionary standpoint, this type of phage-host interaction is consistent with the Piggyback-the-Winner model, which suggests a preference towards lysogenic or other "benign" forms of phage infection when the host is stably present at high abundance. CONCLUSION: Long-term persistence of bacteriophage and host could result from mutually beneficial mechanisms driving bacterial strain-level diversity and phage survival in complex environments.

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.

Phase-variable capsular polysaccharides and lipoproteins modify bacteriophage susceptibility in Bacteroides thetaiotaomicron.

A variety of cell surface structures dictate interactions between bacteria and their environment, including their viruses (bacteriophages). Members of the human gut Bacteroidetes characteristically produce several phase-variable capsular polysaccharides (CPSs), but their contributions to bacteriophage interactions are unknown. To begin to understand how CPSs have an impact on Bacteroides-phage interactions, we isolated 71 Bacteroides thetaiotaomicron-infecting bacteriophages from two locations in the United States. Using B. thetaiotaomicron strains that express defined subsets of CPSs, we show that CPSs dictate host tropism for these phages and that expression of non-permissive CPS variants is selected under phage predation, enabling survival. In the absence of CPSs, B. thetaiotaomicron escapes bacteriophage predation by altering expression of eight distinct phase-variable lipoproteins. When constitutively expressed, one of these lipoproteins promotes resistance to multiple bacteriophages. Our results reveal important roles for Bacteroides CPSs and other cell surface structures that allow these bacteria to persist under bacteriophage predation, and hold important implications for using bacteriophages therapeutically to target gut symbionts.

The crAss-like Phage Group: How Metagenomics Reshaped the Human Virome.

Metagenomics is currently the primary means for identifying new viruses. One of the most impactful metagenomic discoveries is that of crAssphage, the most abundant human-associated virus that is found in about 50% of human gut viromes where it can comprise up to 90% of the virus sequences. Although initial genome analysis of crAssphage failed to detect related phages, or functionally annotate most of the genes, subsequent reanalysis with powerful computational methods and larger databases led to the identification of an expansive group of crAss-like phages. The functions of most crAssphage proteins were predicted, including unusual ones such as giant RNA polymerase polyproteins. The host range of the crAss-like phages consists of various members of the bacterial phylum Bacteroidetes as demonstrated by CRISPR spacer analysis and by analysis of genes acquired by phages from the hosts. New metagenomic studies vastly expanded the crAss-like phage group and demonstrated its global spread and ancient association with primates. The first members of the crAss-like group was recently isolated and shown to infect the bacterium Bacteroides intestinales. Characterization of this phage validated the predicted podovirus-like virion structure and the identity of the major capsid protein and other predicted virion proteins, including three RNA polymerase subunits.

Acquisition, transmission and strain diversity of human gut-colonizing crAss-like phages.

CrAss-like phages are double-stranded DNA viruses that are prevalent in human gut microbiomes. Here, we analyze gut metagenomic data from mother-infant pairs and patients undergoing fecal microbiota transplantation to evaluate the patterns of acquisition, transmission and strain diversity of crAss-like phages. We find that crAss-like phages are rarely detected at birth but are increasingly prevalent in the infant microbiome after one month of life. We observe nearly identical genomes in 50% of cases where the same crAss-like clade is detected in both the mother and the infant, suggesting vertical transmission. In cases of putative transmission of prototypical crAssphage (p-crAssphage), we find that a subset of strains present in the mother are detected in the infant, and that strain diversity in infants increases with time. Putative tail fiber proteins are enriched for nonsynonymous strain variation compared to other genes, suggesting a potential evolutionary benefit to maintaining strain diversity in specific genes. Finally, we show that p-crAssphage can be acquired through fecal microbiota transplantation.

Biogeographic study of human gut-associated crAssphage suggests impacts from industrialization and recent expansion.

CrAssphage (cross-assembly phage) is a bacteriophage that was first discovered in human gut metagenomic data. CrAssphage belongs to a diverse family of crAss-like bacteriophages thought to infect gut commensal bacteria belonging to Bacteroides species. However, not much is known about the biogeography of crAssphage and whether certain strains are associated with specific human populations. In this study, we screened publicly available human gut metagenomic data from 3,341 samples for the presence of crAssphage sensu stricto (NC_024711.1). We found that crAssphage prevalence is low in traditional, hunter-gatherer populations, such as the Hadza from Tanzania and Matses from Peru, as compared to industrialized, urban populations. Statistical comparisons showed no association of crAssphage prevalence with variables such as age, sex, body mass index, and health status of individuals. Phylogenetic analyses show that crAssphage strains reconstructed from the same individual over multiple time-points, cluster together. CrAssphage strains from individuals from the same study population do not always cluster together. Some evidence of clustering is seen at the level of broadly defined geographic regions, however, the relative positions of these clusters within the crAssphage phylogeny are not well-supported. We hypothesize that this lack of strong biogeographic structuring is suggestive of an expansion event within crAssphage. Using a Bayesian dating approach, we estimate that this expansion has occurred fairly recently. Overall, we determine that crAssphage presence is associated with an industrialized lifestyle and the absence of strong biogeographic structuring within global crAssphage strains is likely due to a recent population expansion within this bacteriophage.

Evolution of BACON Domain Tandem Repeats in crAssphage and Novel Gut Bacteriophage Lineages.

The human gut contains an expanse of largely unstudied bacteriophages. Among the most common are crAss-like phages, which were predicted to infect Bacteriodetes hosts. CrAssphage, the first crAss-like phage to be discovered, contains a protein encoding a Bacteroides-associated carbohydrate-binding often N-terminal (BACON) domain tandem repeat. Because protein domain tandem repeats are often hotspots of evolution, BACON domains may provide insight into the evolution of crAss-like phages. Here, we studied the biodiversity and evolution of BACON domains in bacteriophages by analysing over 2 million viral contigs. We found a high biodiversity of BACON in seven gut phage lineages, including five known crAss-like phage lineages and two novel gut phage lineages that are distantly related to crAss-like phages. In three BACON-containing phage lineages, we found that BACON domain tandem repeats were associated with phage tail proteins, suggestive of a possible role of these repeats in host binding. In contrast, individual BACON domains that did not occur in tandem were not found in the proximity of tail proteins. In two lineages, tail-associated BACON domain tandem repeats evolved largely through horizontal transfer of separate domains. In the third lineage that includes the prototypical crAssphage, the tandem repeats arose from several sequential domain duplications, resulting in a characteristic tandem array that is distinct from bacterial BACON domains. We conclude that phage tail-associated BACON domain tandem repeats have evolved in at least two independent cases in gut bacteriophages, including in the widespread gut phage crAssphage.

The Human Gut Virome Is Highly Diverse, Stable, and Individual Specific.

The human gut contains a vast array of viruses, mostly bacteriophages. The majority remain uncharacterized, and their roles in shaping the gut microbiome and in impacting on human health remain poorly understood. We performed longitudinal metagenomic analysis of fecal viruses in healthy adults that reveal high temporal stability, individual specificity, and correlation with the bacterial microbiome. Using a database-independent approach that uses most of the sequencing data, we uncovered the existence of a stable, numerically predominant individual-specific persistent personal virome. Clustering of viral genomes and de novo taxonomic annotation identified several groups of crAss-like and Microviridae bacteriophages as the most stable colonizers of the human gut. CRISPR-based host prediction highlighted connections between these stable viral communities and highly predominant gut bacterial taxa such as Bacteroides, Prevotella, and Faecalibacterium. This study provides insights into the structure of the human gut virome and serves as an important baseline for hypothesis-driven research.

Genomic variation and strain-specific functional adaptation in the human gut microbiome during early life.

The human gut microbiome matures towards the adult composition during the first years of life and is implicated in early immune development. Here, we investigate the effects of microbial genomic diversity on gut microbiome development using integrated early childhood data sets collected in the DIABIMMUNE study in Finland, Estonia and Russian Karelia. We show that gut microbial diversity is associated with household location and linear growth of children. Single nucleotide polymorphism- and metagenomic assembly-based strain tracking revealed large and highly dynamic microbial pangenomes, especially in the genus Bacteroides, in which we identified evidence of variability deriving from Bacteroides-targeting bacteriophages. Our analyses revealed functional consequences of strain diversity; only 10% of Finnish infants harboured Bifidobacterium longum subsp. infantis, a subspecies specialized in human milk metabolism, whereas Russian infants commonly maintained a probiotic Bifidobacterium bifidum strain in infancy. Groups of bacteria contributing to diverse, characterized metabolic pathways converged to highly subject-specific configurations over the first two years of life. This longitudinal study extends the current view of early gut microbial community assembly based on strain-level genomic variation.

ΦCrAss001 represents the most abundant bacteriophage family in the human gut and infects Bacteroides intestinalis.

CrAssphages are an extensive and ubiquitous family of tailed bacteriophages, predicted to infect bacteria of the order Bacteroidales. Despite being found in ~50% of individuals and representing up to 90% of human gut viromes, members of this viral family have never been isolated in culture and remain understudied. Here, we report the isolation of a CrAssphage (ΦCrAss001) from human faecal material. This bacteriophage infects the human gut symbiont Bacteroides intestinalis, confirming previous in silico predictions of the likely host. DNA sequencing demonstrates that the bacteriophage genome is circular, 102 kb in size, and has unusual structural traits. In addition, electron microscopy confirms that ΦcrAss001 has a podovirus-like morphology. Despite the absence of obvious lysogeny genes, ΦcrAss001 replicates in a way that does not disrupt proliferation of the host bacterium, and is able to maintain itself in continuous host culture during several weeks.

A diversity-generating retroelement encoded by a globally ubiquitous Bacteroides phage.

BACKGROUND: Diversity-generating retroelements (DGRs) are genetic cassettes that selectively mutate target genes to produce hypervariable proteins. First characterized in Bordetella bacteriophage BPP-1, the DGR creates a hypervariable phage tail fiber that enables host tropism switching. Subsequent surveys for DGRs conclude that the majority identified to date are bacterial or archaeal in origin. This work examines bacteriophage and bacterial genomes for novel phage-encoded DGRs. RESULTS: This survey discovered 92 DGRs that were only found in phages exhibiting a temperate lifestyle. The majority of phage-encoded DGRs were identified as prophages in bacterial hosts from the phyla Bacteroidetes, Proteobacteria, and Firmicutes. Sequence reads from these previously unidentified prophages were present in viral metagenomes (viromes), indicating these prophages can produce functional viruses. Five phages possessed hypervariable proteins with structural similarity to the tail fiber of BPP-1, whereas the functions of the remaining DGR target proteins were unknown. A novel temperate phage that harbors a DGR cassette targeting a protein of unknown function was induced from Bacteroides dorei. This phage, here named Bacteroides dorei Hankyphage, lysogenizes 13 different Bacteroides species and was present in 34% and 21% of whole-community metagenomes and human-associated viromes, respectively. CONCLUSIONS: Here, the number of known DGR-containing phages is increased from four to 92. All of these phages exhibit a temperate lifestyle, including a cosmopolitan human-associated phage. Targeted hypervariation by temperate phages may be a ubiquitous mechanism underlying phage-bacteria interaction in the human microbiome.

Isolation of Bacteriophages of the Anaerobic Bacteria Bacteroides.

Here we describe the detection, enumeration, and isolation of bacteriophages infecting Bacteroides. The method is based on the infection of Bacteroides host strains and the production of visible plaques in a confluent lawn of the host strain using the double-layer agar method. This is a straightforward methodology that can be applied for the detection, enumeration and isolation of bacteriophages for other anaerobic bacteria, using an appropriate host strain and culture conditions. In the case of bacteriophages of Bacteroides the results can be obtained in less than 24 h, although the time could vary depending on the growth rate of the host strain.

Bacteriophages as indicators of faecal pollution and enteric virus removal.

Bacteriophages are an attractive alternative to faecal indicator bacteria (FIB), particularly as surrogates of enteric virus fate and transport, due to their closer morphological and biological properties. Based on a review of published data, we summarize densities of coliphages (F+ and somatic), Bacteroides spp. and enterococci bacteriophages (phages) in individual human waste, raw wastewater, ambient fresh and marine waters and removal through wastewater treatment processes utilizing traditional treatments. We also provide comparisons with FIB and enteric viruses whenever possible. Lastly, we examine fate and transport characteristics in the aquatic environment and provide an overview of the environmental factors affecting their survival. In summary, concentrations of bacteriophages in various sources were consistently lower than FIB, but more reflective of infectious enteric virus levels. Overall, our investigation indicates that bacteriophages may be adequate viral surrogates, especially in built systems, such as wastewater treatment plants. SIGNIFICANCE AND IMPACT OF THE STUDY: Bacteriophage are alternative fecal indicators that may be better surrogates for viral pathogens than fecal indicator bacteria (FIB). This report offers a summary of the existing literature concerning the utility of bacteriophage as indicators of viral presence (fecal sources and surface waters) and persistence (in built infrastructure and aquatic environments). Our findings indicate that bacteriophage levels in all matrices examined are consistently lower than FIB, but similar to viral pathogens. Furthermore, in built infrastructure (e.g. wastewater treatment systems) bacteriophage closely mimic viral pathogen persistence suggesting they may be adequate sentinels of enteric virus removal.

Imbalance of bacteriome profiles within the Finnish Diabetes Prediction and Prevention study: Parallel use of 16S profiling and virome sequencing in stool samples from children with islet autoimmunity and matched controls.

BACKGROUND: We set out to explore associations between the stool bacteriome profiles and early-onset islet autoimmunity, taking into account the interactions with the virus component of the microbiome. METHODS: Serial stool samples were longitudinally collected from 18 infants and toddlers with early-onset islet autoimmunity (median age 17.4 months) followed by type 1 diabetes, and 18 tightly matched controls from the Finnish Diabetes Prediction and Prevention (DIPP) cohort. Three stool samples were analyzed, taken 3, 6, and 9 months before the first detection of serum autoantibodies in the case child. The risk of islet autoimmunity was evaluated in relation to the composition of the bacteriome 16S rDNA profiles assessed by mass sequencing, and to the composition of DNA and RNA viromes. RESULTS: Four operational taxonomic units were significantly less abundant in children who later on developed islet autoimmunity as compared to controls-most markedly the species of Bacteroides vulgatus and Bifidobacterium bifidum. The alpha or beta diversity, or the taxonomic levels of bacterial phyla, classes or genera, showed no differences between cases and controls. A correlation analysis suggested a possible relation between CrAssphage signals and quantities of Bacteroides dorei. No apparent associations were seen between development of islet autoimmunity and sequences of yet unknown origin. CONCLUSIONS: The results confirm previous findings that an imbalance within the prevalent Bacteroides genus is associated with islet autoimmunity. The detected quantitative relation of the novel "orphan" bacteriophage CrAssphage with a prevalent species of the Bacteroides genus may exemplify possible modifiers of the bacteriome.

Occurrence of bacteriophages infecting Bacteroides host strains (ARABA 84 and GB-124) in fecal samples of human and animal origin.

Bacteriophage-based microbial source-tracking studies are an economical and simple way of identifying fecal sources in polluted water systems. Recently isolated Bacteroides spp. strains ARABA 84, and GB-124 have been shown to detect bacteriophages exclusively in aquatic systems impacted by human fecal material. To date, limited examination of the occurrence or concentration of phages capable of infecting Bacteroides fragilis strain GB-124 or B. thetaiotaomicron strain ARABA 84 in human and animal feces has been carried out. This study reports the prevalence rates and concentrations of phages infecting ARABA 84 and GB-124 host strains in human and a range of animal feces. Discrete human fecal samples (n=55) and pooled animal samples (n=46, representing the feces of over 230 animals) were examined for phages infecting the host strains ARABA 84, GB-124, and E. coli strain WG5. Both human Bacteroides host strains were highly specific (95% and 100% for ARABA 84 and GB-124, respectively), challenging results from previous studies. This study supports the use of Bacteroides strains GB-124 and ARABA 84 in fecal source tracking studies for the detection of human fecal contamination.

A highly abundant bacteriophage discovered in the unknown sequences of human faecal metagenomes.

Metagenomics, or sequencing of the genetic material from a complete microbial community, is a promising tool to discover novel microbes and viruses. Viral metagenomes typically contain many unknown sequences. Here we describe the discovery of a previously unidentified bacteriophage present in the majority of published human faecal metagenomes, which we refer to as crAssphage. Its ~97 kbp genome is six times more abundant in publicly available metagenomes than all other known phages together; it comprises up to 90% and 22% of all reads in virus-like particle (VLP)-derived metagenomes and total community metagenomes, respectively; and it totals 1.68% of all human faecal metagenomic sequencing reads in the public databases. The majority of crAssphage-encoded proteins match no known sequences in the database, which is why it was not detected before. Using a new co-occurrence profiling approach, we predict a Bacteroides host for this phage, consistent with Bacteroides-related protein homologues and a unique carbohydrate-binding domain encoded in the phage genome.

Bacteriophages infecting Bacteroides as a marker for microbial source tracking.

Bacteriophages infecting certain strains of Bacteroides are amid the numerous procedures proposed for tracking the source of faecal pollution. These bacteriophages fulfil reasonably well most of the requirements identified as appropriate for a suitable marker of faecal sources. Thus, different host strains are available that detect bacteriophages preferably in water contaminated with faecal wastes corresponding to different animal species. For phages found preferably in human faecal wastes, which are the ones that have been more extensively studied, the amounts of phages found in waters contaminated with human fecal samples is reasonably high; these amounts are invariable through the time; their resistance to natural and anthropogenic stressors is comparable to that of other relatively resistant indicator of faecal pollution such us coliphages; the abundance ratios of somatic coliphages and bacteriophages infecting Bacteroides thetaiotaomicron GA17 are unvarying in recent and aged contamination; and standardised detection methods exist. These methods are easy, cost effective and provide data susceptible of numerical analysis. In contrast, there are some uncertainties regarding their geographical stability, and consequently suitable hosts need to be isolated for different geographical areas. However, a feasible method has been described to isolate suitable hosts in a given geographical area. In summary, phages infecting Bacteroides are a marker of faecal sources that in our opinion merits being included in the "toolbox" for microbial source tracking. However, further research is still needed in order to make clear some uncertainties regarding some of their characteristics and behaviour, to compare their suitability to the one of emerging methods such us targeting Bacteroidetes by qPCR assays; or settling molecular methods for their determination.

Gnotobiotic mouse model of phage-bacterial host dynamics in the human gut.

Bacterial viruses (phages) are the most abundant biological group on Earth and are more genetically diverse than their bacterial prey/hosts. To characterize their role as agents shaping gut microbial community structure, adult germ-free mice were colonized with a consortium of 15 sequenced human bacterial symbionts, 13 of which harbored one or more predicted prophages. One member, Bacteroides cellulosilyticus WH2, was represented by a library of isogenic transposon mutants that covered 90% of its genes. Once assembled, the community was subjected to a staged phage attack with a pool of live or heat-killed virus-like particles (VLPs) purified from the fecal microbiota of five healthy humans. Shotgun sequencing of DNA from the input pooled VLP preparation plus shotgun sequencing of gut microbiota samples and purified fecal VLPs from the gnotobiotic mice revealed a reproducible nonsimultaneous pattern of attack extending over a 25-d period that involved five phages, none described previously. This system allowed us to (i) correlate increases in specific phages present in the pooled VLPs with reductions in the representation of particular bacterial taxa, (ii) provide evidence that phage resistance occurred because of ecological or epigenetic factors, (iii) track the origin of each of the five phages among the five human donors plus the extent of their genome variation between and within recipient mice, and (iv) establish the dramatic in vivo fitness advantage that a locus within a B. cellulosilyticus prophage confers upon its host. Together, these results provide a defined community-wide view of phage-bacterial host dynamics in the gut.

Use of abundance ratios of somatic coliphages and bacteriophages of Bacteroides thetaiotaomicron GA17 for microbial source identification.

Water contaminated with human faeces is a risk to human health and management of water bodies can be improved by determining the sources of faecal pollution. Field studies show that existing methods are insufficient and that different markers are required. This study proposes the combined use of two microbial indicators, the concentrations of which are presented as ratios. This provides a more reliable approach to identifying faecal sources as it avoids variation due to treatment or ageing of the contamination. Among other indicators, bacteriophages have been proposed as rapid and cheap indicators of faecal pollution. Samples analysed in this study were derived from wastewater treatment plants (raw sewage, secondary and tertiary effluents and raw sewage sludge) river water, seawater and animal related wastewater. The abundance ratios of faecal coliforms and Bacteroides phages, either strain RYC2056 (non-specific for faecal origin) or strain GA17 (specific for human pollution), and among somatic coliphages and phages infecting both Bacteroides strains, were evaluated. The results indicate that the ratio of somatic coliphages and phages infecting Bacteroides strain GA17, which is specific to human faecal sources, provides a robust method for discriminating samples, even those presenting different levels and ages of pollution, and allows samples polluted with human faeces to be distinguished from those containing animal faecal pollution. This method allows the generation of numerical data that can be further applied to numerical methods for faecal pollution discrimination.