Tracking microbial colonization in fecal microbiota transplantation experiments via genome-resolved metagenomics
Microbiome volume 5, Article number: 50 (2017)
Fecal microbiota transplantation (FMT) is an effective treatment for recurrent Clostridium difficile infection and shows promise for treating other medical conditions associated with intestinal dysbioses. However, we lack a sufficient understanding of which microbial populations successfully colonize the recipient gut, and the widely used approaches to study the microbial ecology of FMT experiments fail to provide enough resolution to identify populations that are likely responsible for FMT-derived benefits.
We used shotgun metagenomics together with assembly and binning strategies to reconstruct metagenome-assembled genomes (MAGs) from fecal samples of a single FMT donor. We then used metagenomic mapping to track the occurrence and distribution patterns of donor MAGs in two FMT recipients.
Our analyses revealed that 22% of the 92 highly complete bacterial MAGs that we identified from the donor successfully colonized and remained abundant in two recipients for at least 8 weeks. Most MAGs with a high colonization rate belonged to the order Bacteroidales. The vast majority of those that lacked evidence of colonization belonged to the order Clostridiales, and colonization success was negatively correlated with the number of genes related to sporulation. Our analysis of 151 publicly available gut metagenomes showed that the donor MAGs that colonized both recipients were prevalent, and the ones that colonized neither were rare across the participants of the Human Microbiome Project. Although our dataset showed a link between taxonomy and the colonization ability of a given MAG, we also identified MAGs that belong to the same taxon with different colonization properties, highlighting the importance of an appropriate level of resolution to explore the functional basis of colonization and to identify targets for cultivation, hypothesis generation, and testing in model systems.
The analytical strategy adopted in our study can provide genomic insights into bacterial populations that may be critical to the efficacy of FMT due to their success in gut colonization and metabolic properties, and guide cultivation efforts to investigate mechanistic underpinnings of this procedure beyond associations.
Fecal microbiota transplantation (FMT), the transference of fecal material from a healthy donor to a recipient, has gained recognition as an effective and relatively safe treatment for recurrent or refractory Clostridium difficile infection (CDI) [1,2,3,4,5,6,7,8]. Its success in treating CDI sparked interest in investigating FMT as a treatment for other medical conditions associated with intestinal dysbiosis, such as ulcerative colitis [9,10,11], Crohn’s disease (CD) [12,13,14], irritable bowel syndrome (IBS) [15, 16], and others, including metabolic syndrome , neurodevelopmental , and autoimmune disorders . Despite the excitement due to its therapeutic potential, FMT also presents challenges for researchers and clinicians with potential adverse outcomes, including the transfer of infectious organisms  or contaminants from the environment [21, 22]. A complete understanding of FMT from a basic science perspective is still lacking, as we have yet to determine the key microbial populations that are responsible for beneficial outcomes, as well as adverse effects.
Recent advances in high-throughput sequencing technologies, molecular approaches, and computation have dramatically increased our ability to investigate the ecology of microbial populations. Utilization of these advances at a proper level of resolution can lead to a better mechanistic understanding of FMT and identify new therapeutic opportunities or address potential risks. Most current studies on FMT use amplicons from marker genes, such as the 16S ribosomal RNA gene, to characterize the composition of microbial communities [23,24,25,26]. While providing valuable insights into the broad characteristics of FMTs, amplicons from the 16S ribosomal RNA gene do not offer the resolution to effectively identify populations that colonize recipients . Other studies use shotgun metagenomics to annotate short reads and map them to reference genomes in order to track changes in the functional potential or membership in the gut microbial communities of recipients [28,29,30]. In a recent study, Li et al.  demonstrated the coexistence of donors and recipients’ gut microbes 3 months after FMT by mapping short metagenomic reads to reference genomes. Although this approach provides more information than marker gene amplicons alone, it is subject to the limitations and biases of reference genomic databases, is unable to characterize populations that do not have closely related culture representatives, and does not provide direct access to the genomic context of relevant populations for more targeted follow-up studies.
Metagenomic assembly and binning [31, 32] is an alternative approach to characterizing microbial communities through marker gene amplicons or reference genomes. Here, we used the state-of-the-art metagenomic assembly and binning strategies to reconstruct microbial population genomes directly from a single FMT donor and tracked the occurrence of resulting metagenome-assembled genomes (MAGs) in two FMT recipients up to 8 weeks.
Sample collection, preparation, and sequencing
We collected a total of 10 fecal samples; four samples from a single donor “D” (a 30-year-old male) and three samples from each of the two recipients “R01” (a 23-year-old male) and “R02” (a 32-year-old female) before and after FMT. Recipient samples originated from time points pre-FMT, 4 weeks after FMT, and 8 weeks after FMT, while four samples from the donor were collected on four separate days 2 weeks prior to the transplantation. All fecal samples were handled under anaerobic conditions prior to transplantation, and the recipients had no genetic relationship to the donor. Through a single colonoscopy for each recipient, the donor sample DS 01 was transferred to R01, and the donor sample DS 04 was transferred to R02. All samples were stored at −80 °C until DNA extraction. We extracted the genomic DNA from frozen samples according to the centrifugation protocol outlined in MoBio PowerSoil kit with the following modifications: cell lysis was performed using a GenoGrinder to physically lyse the samples in the MoBio Bead Plates and Solution (5–10 min). After final precipitation, the DNA samples were resuspended in TE buffer and stored at −20 °C until further analysis. We prepared our shotgun metagenomic libraries with OVATION Ultralow protocol (NuGen) and used an Illumina NextSeq 500 platform to generate 2 × 150 nt paired-end sequencing reads.
Metagenomic assembly and binning
We removed the low-quality reads from the raw sequencing results using the program “iu-filter-quality-minoche” in illumina-utils  (available from https://github.com/merenlab/illumina-utils) according to Minoche et al. . We then co-assembled reads from the donor samples using MEGAHIT v1.0.6 , used Centrifuge v1.0.2-beta  to remove contigs that matched to human genome, and mapped short reads from each recipient and donor sample to the remaining contigs using Bowtie2 v2.0.5 . We then used anvi’o v2.3.1 (available from http://merenlab.org/software/anvio) to profile mapping results, finalize genomic bins, and visualize results following the workflow outlined in Eren et al. . Briefly, (1) the program “anvi-gen-contigs-database” profiled our contigs using Prodigal v2.6.3  with default settings to identify open reading frames and HMMER  to identify matching genes in our contigs to bacterial  and archaeal  single-copy core gene collections, (2) “anvi-init-bam” converted mapping results into BAM files, (3) “anvi-profile” processed each BAM file to estimate the coverage and detection statistics of each contig using samtools , and finally, (4) “anvi-merge” combined profiles from each sample to create a merged anvi’o profile for our dataset. We used “anvi-cluster-with-concoct” for the initial binning of contigs using CONCOCT  by constraining the number of clusters to 10 (“--num-clusters 10”) to minimize the “fragmentation error” (where multiple bins describe one population). We then interactively refined each CONCOCT bin that exhibit “conflation error” (where one bin describes multiple populations) using the program “anvi-refine” based on tetra-nucleotide frequency, taxonomy, mean coverage, and completion and redundancy estimates based on bacterial and archaeal single-copy genes. We classified a given genome bin as a “metagenome-assembled genome” (MAG) if it was more than 70% complete or larger than 2 Mbp, and its redundancy was estimated to be less than 10%. We used “anvi-interactive” to visualize the distribution of our bins across samples and “anvi-summarize” to generate static HTML output for binning results. Besides anvi’o, we also used CheckM v1.0.7  to assess the completion and contamination of our bins.
Taxonomic and functional annotation of MAGs
We employed multiple approaches to infer taxonomy. Besides the taxonomic annotations reported by CheckM, we searched amino acid sequences for the RecA gene for each MAG in the National Center for Biotechnology Information (NCBI) databases. We also used Phylosift v1.0.1  to determine the phylogenomic relationships between our MAGs and a collection of 1758 reference genomes (with no redundancy at the species level) that we acquired from the Ensembl database  (Additional file 1: Table S3). Briefly, Phylosift (1) identifies a set of 37 marker gene families in each genome, (2) concatenates the alignment of each marker gene family across genomes, and (3) computes a phylogenomic tree from the concatenated alignment using FastTree . Lastly, we used FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/) to finalize the phylogenomic tree for publication. We used RAST  to ascribe functions to our MAGs.
Criteria for detection and colonization of MAGs in the recipients
For each genome bin, anvi’o reports the percentage of nucleotide positions in all contigs that are covered by at least one short read based on mapping results. This statistic gives an estimate of “detection” regardless of the coverage of a given genome bin. We required the detection statistic of a genome bin to be at least 25% to consider it “detected” in a given sample. This prevented inflated detection rates due to non-specific mapping, which is not uncommon due to relatively well-conserved genes across gut populations. Finally, we conservatively decided that a MAG was transferred from the donor and colonized a given recipient successfully only if (1) it was detected in both samples that were collected from the recipient at 4 and 8 weeks after the FMT and (2) it was not detected in the pre-FMT sample from the same recipient.
The use of HMP metagenomes
We used 151 Human Microbiome Project (HMP) gut metagenomes  to estimate the detection of our MAGs and to compare the taxonomic profiles of our metagenomes in the context of the HMP participants. To estimate detection, we mapped HMP metagenomes to our MAGs using Bowtie2 with default parameters and considered a MAG to be detected in a given HMP metagenome when its level of detection surpassed 25%. We also annotated our metagenomes and HMP gut metagenomes using MetaPhlAn2 . Additionally, to estimate the contribution of donor MAGs to recipient taxonomic profiles, we used the 60 and 83 MAGs that were not detected pre-FMT in samples of R01 and R02 to recruit and remove reads from the post-FMT R01 and R02 metagenomes. MetaphlAn2 estimated the taxonomical profiles of the remaining reads. Additional file 1: Table S3 reports taxonomic annotations.
We performed cluster analyses on distribution profiles of MAGs using the R library vegan with Bray-Curtis distances of normalized values. We used the PERMANOVA (R adonis vegan)  test to measure the degree of similarity of the bacterial communities between the samples in the study. We further used similarity index (SIMPER) analysis to identify the taxa that contributed the highest dissimilarity between the samples. We classified the MAGs into four main groups based on their colonization characteristics in the recipients. We then performed a Wilcoxon signed-rank test (STAMP)  with Benjamini-Hochberg FDR (false discovery rate) correction for multiple tests on the total-sum normalized data to ascertain any significant differences in the functional potential between the groups and carried out canonical correspondence analysis based on functional potential and the MAGs’ colonization characteristics.
The shotgun sequencing of genomic DNA from 10 fecal samples resulted in a total of 269,144,211 quality-filtered 2 × 150 paired-end metagenomic reads (Additional file 2: Table S1). By co-assembling the donor samples, which corresponded to 115,037,928 of the quality-filtered reads, we recovered 51,063 contigs that were longer than 2.5 kbp and organized them into 444 genomic bins comprising a total of 442.64 Mbp at various levels of completion (Additional file 2: Table S1, Additional file 3: Figure S1). Using completion and size criteria, we designated 92 of our genomic bins as metagenome-assembled genomes (MAGs) (Fig. 1, Additional file 2: Table S1). Four major patterns emerged from the distribution of MAGs across individuals: MAGs that colonized both recipients R01 and R02 (group I, n = 20), MAGs that colonized only R01 (group II, n = 11), only R02 (group III, n = 8) and MAGs that did not colonize either of the recipients (group IV, n = 13) (Fig. 1). We found no correlation between the abundances of MAGs in donor samples and their success at colonizing recipients (ANOVA, F = 0.717, p = 0.543). Additional file 2: Table S1 reports the detection and mean coverage statistics for each MAG in each group. Both of our recipients had mild/moderate ulcerative colitis (R01 had proctitis and R02 had left-sided/distal colitis), and both of them showed a decrease in intestinal inflammation after FMT. R01 had an initial fecal calprotectin count of 26 prior to FMT, which decreased to <17 at week 4 (no data was collected at week 8), and R02’s fecal calprotectin count decreased from 947 (pre-FMT) to 68 and 34 in week 4 (W4) and week 8 (W8), respectively.
The taxonomy of 14 of the 20 MAGs that colonized both recipients resolved to the order Bacteroidales (Fig. 1). Besides Bacteroidales, group I also included five MAGs that were classified as order Clostridiales and one MAG as Coriobacteriales. CheckM partitioned the group I MAGs into Bacteroides (n = 5), Alistipes (n = 5), Odoribacter (n = 1), Paraprevotella (n = 1), and Barnesiella (n = 1). Seven MAGs were not assigned to a specific genus in this group. In contrast to the Bacteroidales-dominated group I, 10 of the 13 MAGs that did not colonize recipients (group IV) resolved to the order Clostridiales. The remaining three MAGs were not assigned any taxonomy at the order level. The only genus-level annotation for the MAGs in group IV was Ruminococcus (n = 2). Overall, CheckM did not assign any genus-level taxonomy to 18 of the 33 MAGs in group I and group IV. MAGs that colonized only one of the two recipients did not show a consistent taxonomic signal: while 9 of 11 MAGs that colonized only R01 (group II) were assigned to the order Clostridiales, only 4 of 8 MAGs that only colonized R02 (group III) were assigned to that order (Fig. 1, Additional file 2: Table S1). The remaining MAGs in group III were assigned to orders Bacteroidales (n = 2) and Burkholderiales (n = 1), and one of them did not resolve to any order-level taxon (n = 1). We also employed a phylogenomic approach to investigate the validity of the taxonomy of our MAGs with respect to reference genomes (Additional file 4: Figure S2). All order-level taxonomic annotations by CheckM were consistent with the phylogenomic placement of our MAGs. We resolved inconsistencies at the genus-level by removing annotations that differed between approaches.
We used the Human Microbiome Project gut metagenomes to investigate whether the differential colonization outcomes we observed for donor MAGs were representative of their occurrence in healthy individuals. Our analysis of 151 publicly available gut metagenomes showed that the donor MAGs in group I that colonized both recipients were more prevalent, and the ones in group IV that colonized neither of the recipients were more rare across the participants of the HMP (ANOVA, F = 10.04, p < 0.001; Fig. 1). In fact, while the group I MAGs occurred in 40.63% (±26.67%) of the individuals, the group IV MAGs showed more sporadic patterns as they were detected in only 10.39(±5.69%) of the HMP metagenomes (Fig. 1, Additional file 1: Table S3, Additional file 5: Figure S3). These results suggest that the relevance of the donor MAGs we recovered here are not necessarily limited to our study and represent a subset of populations with patterns of high- and low-colonization success.
We then investigated whether there was a link between the functional potential of MAGs and their success of colonization. The canonical correspondence analysis (CCA) of 500 functions and 107 sub-systems (Additional file 6: Table S2) revealed that the group I MAGs possessed a higher relative abundance of genes coding for quinone cofactors, along with functions involving lipoic acid synthesis and metabolism of aromatic compounds. In contrast, the group IV MAGs carried higher number of genes related to dormancy and sporulation, spore DNA protection, and motility and chemotaxis (pseudo-F = 2.156, p < 0.0001). Interestingly, MAGs that were transferred to only R01 (group II) also carried higher numbers of genes for dormancy and sporulation and spore DNA protection (Fig. 2, Additional file 6: Table S2). The functional similarity between group II and group IV MAGs was also mirrored with their detection rate across the HMP metagenomes (Additional file 1: Table S3, Additional file 5: Figure S3). This suggests that despite strong signal, the functional potential is not a robust predictor of colonization success and individuals may respond differently to the same donor.
We also investigated to what extent the recipients became donor-like after the FMT. As expected, the analysis based on the mean coverage of donor MAGs suggested an increased similarity between the donor and recipients after the FMT (Fig. 3a). Mapping metagenomic short reads to donor MAGs offers a limited insight into the recipient microbial populations. To minimize any potential biases, we also generated genus-level taxonomic profiles for donor and recipient samples by annotating all metagenomic short reads. This strategy also allowed us to compare our donor, as well as our recipients before and after FMT, to the larger context of the HMP cohort (Fig. 3b). This analysis suggested that the donor samples showed no discernible differences from the HMP cohort (PERMANOVA, pseudo-F = 1.489, p = 0.110). In contrast, the pre-FMT samples of both recipients differed significantly from the HMP cohort and from the donor samples (PERMANOVA, pseudo-F = 4.470, p = 0.001). However, the recipient samples after FMT no longer differed from the HMP cohort (PERMANOVA, pseudo-F = 1.395, p = 0.168) (Fig. 3b). Based on the taxonomic profiles from short metagenomic reads, both recipients were more than 60% similar to the donor microbiota after FMT. Similarity percentage analysis (SIMPER) of the taxonomic profiles suggested that the two recipients were 79.60% similar after FMT and that Bacteroides was responsible for the largest fraction (32.42%) of difference between recipient samples of pre-FMT and 4 weeks after FMT. There were no significant changes in the recipients’ taxonomic profiles between W4 and W8 after FMT (PERMANOVA, pseudo-F = 0.221, p = 0.631). These results show that our recipients, who were outliers to the HMP cohort, became more like the HMP cohort following FMT. Yet, one important question remains: to what extent this convergence is due to donor MAGs and not due to the changing abundance of the initially rare members of the recipient microbial populations? To investigate this, we first identified donor MAGs that were not detected in a given recipient’s pre-FMT metagenome and removed short reads from the recipient’s W4 and W8 metagenomes after FMT. We then characterized the taxonomic profiles with the remaining reads (Fig. 3c). Without reads matching to donor MAGs, the taxonomic profiles for the samples collected after FMT remained close to the pre-FMT status, and dissimilar to both the donor and HMP metagenomes, suggesting a more important contribution of the donor populations to the convergence rather than the emergence of initially rare recipient populations. This was particularly apparent for R01, for which the taxonomic profiles before and after FMT at W8 were highly dissimilar (R 2 of 0.122) yet became more alike after the subtraction of reads matching to donor MAGs (R 2 of 0.634). Taking into account the fact that 92 MAGs are only a subset of the 444 bins we initially recovered from the donor, we tested whether these trends remained similar when all 444 metagenomic bins were used for removal of short reads and confirmed that the outcome did not change (Additional file 7: Figure S4).
Here, we demonstrate that shotgun metagenomics can facilitate the tracking of bacterial populations in FMT experiments by linking their colonization trends to genomic contexts through metagenome-assembled genomes (MAGs). In our study, the relative abundance of MAGs in donor metagenomes did not predict whether they would colonize the FMT recipients or not. However, their success in colonizing the recipients mirrored their occurrence in the large cohort of the Human Microbiome Project (HMP): the MAGs that colonized both of our recipients were also found in many of the participants of the HMP, and the ones that colonized neither of our recipients were missing in many. Besides suggesting future directions to investigate the ecological and functional basis of gut colonization, this observation suggest that perhaps large-scale metagenomic surveys can be useful to predict the colonization properties of bacterial populations rapidly.
Previous studies reported an increase in the relative abundance of Alistipes [23,55,, 24, 54–56] and Bacteroides populations after FMT experiments [23,24,25,26, 30]. The success of the order Bacteroidales was also striking in our dataset: 14 of the 20 group I Bacteroidales MAGs we identified in the donor successfully colonized both recipient guts. Although the taxonomic signal was relatively strong, our results also showed that taxonomy is not the sole predictor of transfer, as MAGs that resolved to the same genera (i.e., Alistipes, Bacteroides, and Clostridium) showed different colonization properties. In addition, taxonomic annotation of a large fraction of MAGs in our study did not resolve to a genus name, which suggests that bacterial populations that have not yet been characterized in culture collections may be playing important roles in FMT treatments.
Although a substantial number of studies report successful medical outcomes of FMT experiments [3, 7, 57, 58], a complete understanding of this procedure from the perspective of microbial ecology is still lacking. Studying FMT as an ecological event and the identification of its key components that facilitate the procedure’s success as a treatment for intestinal disorders require a detailed characterization of the transferred microbial populations at an appropriate level of resolution. MAGs reconstructed directly from donor samples can provide enough resolution to guide cultivation efforts. A recent effort by Vineis et al.  demonstrated this principle by first identifying bacterial populations of interest using MAGs reconstructed from a gut metagenome and then using the genomic context of those MAGs to screen culture experiments from the same gut sample to bring the target population to the bench. A similar approach in the context of FMTs can provide opportunities to design experiments to explore the functional basis of colonization in controlled systems.
The complete transfer of fecal matter between individuals comes with various risks. For instance, a recent meta-analysis of 50 peer-reviewed FMT case reports reported 38 potentially transfer-related adverse effects in FMT patients in 35 studies, including fever, sore throat, vomiting, abdominal pain, bowel perforation, rhinorrhea, transient relapse of UC and CDI, and in one case, death, due to temporary systemic immune response to the applied bacteria . Besides bacteria, FMT can transfer viruses, archaea, and fungi, as well as other agents of the donor host such as colonocytes , which may affect the recipient’s biology in unexpected ways. A more complete understanding of the microbial ecology of FMTs could identify precisely what needs to be transferred so that recipients benefit from the positive outcomes of FMT without incurring medical risks from uncharacterized biological material.
A recent study by Khanna et al.  reported high rates of success with the treatment of patients with primary Clostridium difficile infection (CDI) using an investigational oral microbiome therapeutic, SER-109, which contains bacterial spores enriched and purified from healthy donors. However, Seres Therapeutics announced more recently that interim findings from the mid-stage clinical study of SER-109 failed to meet their primary goal of reducing the risk of recurrence for up to 8 weeks . In our study, the MAGs that failed to colonize any of the recipients were significantly enriched for spore-formation genes, and they also showed a very sporadic distribution across the HMP cohort. Interestingly, Nayfach et al.  recently made a similar observation regarding the transmission of bacteria and sporulation in a different system: vertical transmission of bacteria between mothers and their infants. Bacterial populations with high-vertical transmission rates had lower number of genes related to sporulation . These observations suggest that excluding non-spore forming bacteria may decrease the efficacy of FMT therapies due to limited colonization efficiency, and deeper insights into the functional basis of microbial colonization require further study.
Using bacterial populations associated with positive health outcomes and that harbor high colonization properties may result in more effective therapies compared to cleansing all but spore-forming bacteria to avoid the transfer of pathogens in FMT experiments. The analytical strategy adopted in our study can provide genomic insights into bacterial populations that may be critical to the efficacy of FMT due to their success in gut colonization and metabolic properties and guide cultivation efforts to investigate mechanistic underpinnings of this procedure beyond associations.
Fecal microbiota transplantation
Rupnik M, Wilcox MH, Gerding DN. Clostridium difficile infection: new developments in epidemiology and pathogenesis. Nat Rev Micro. 2009;7:526–36.
Rohlke F, Surawicz CM, Stollman N. Fecal flora reconstitution for recurrent Clostridium difficile infection: results and methodology. J Clin Gastroenterol. 2010;44:567–70.
Gough E, Shaikh H, Manges AR. Systematic review of intestinal microbiota transplantation (fecal bacteriotherapy) for recurrent Clostridium difficile infection. Clin Infect Dis. 2011;53:994–1002.
Landy J, Al-Hassi HO, McLaughlin SD, Walker AW, Ciclitira PJ, Nicholls RJ, et al. Review article: faecal transplantation therapy for gastrointestinal disease. Aliment Pharmacol Ther. 2011;34:409–15.
Rohlke F, Stollman N. Fecal microbiota transplantation in relapsing Clostridium difficile infection. Therap Adv Gastroenterol. 2012;5:403–20.
Kelly CP. Fecal microbiota transplantation—an old therapy comes of age. N Engl J Med. 2013;368:474–5.
Kassam Z, Lee CH, Yuan Y, Hunt RH. Fecal microbiota transplantation for Clostridium difficile infection: systematic review and meta-analysis. Am J Gastroenterol. 2013;108:500–8.
Lübbert C, John E, von Müller L. Clostridium difficile infection: guideline-based diagnosis and treatment. Dtsch Ärzteblatt Int. 2014;111:723–31.
Anderson JL, Edney RJ, Whelan K. Systematic review: faecal microbiota transplantation in the management of inflammatory bowel disease. Aliment Pharmacol Ther. 2012;36:503–16.
Kunde S, Pham A, Bonczyk S, Crumb T, Duba M, Conrad HJ, et al. Safety, tolerability, and clinical response after fecal transplantation in children and young adults with ulcerative colitis. J Pediatr Gastroenterol Nutr. 2013;56:597–601.
Moayyedi P, Surette MG, Kim PT, Libertucci J, Wolfe M, Onischi C, et al. Fecal microbiota transplantation induces remission in patients with active ulcerative colitis in a randomized controlled trial. Gastroenterology. 2015;149:102–9. e6.
Zhang F, Luo W, Shi Y, Fan Z, Ji G. Should we standardize the 1,700-year-Old fecal microbiota transplantation. Am J Gastroenterol. 2012;107:1755.
Kao D, Hotte N, Gillevet P, Madsen K. Fecal microbiota transplantation inducing remission in Crohn’s colitis and the associated changes in fecal microbial profile. J Clin Gastroenterol. 2014;48:625–8.
Gordon H, Harbord M. A patient with severe Crohn’s colitis responds to faecal microbiota transplantation. J Crohn's & Colitis. 2014;8:256. LP-257.
Borody TJ, George L, Andrews P, Brandi S, Noonan S, Cole P, et al. Bowel-flora alteration: a potential cure for inflammatory bowel disease and irritable bowel syndrome? Med J Aust. 1989;150(10):604.
Andrews P, Borody T, Shortis N, Thompson S. Bacteriotherapy for chronic constipation—a long term follow-up. Gastroenterology. 1995;108:A563.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4.
Finegold SM, Molitoris D, Song Y, Liu C, Vaisanen M-L, Bolte E, et al. Gastrointestinal microflora studies in late-onset autism. Clin Infect Dis. 2002;35:S6–16.
Çorapçioğlu D, Tonyukuk V, Kiyan M, Yilmaz AE, Emral R, Kamel N, et al. Relationship between thyroid autoimmunity and yersinia enterocolitica antibodies. Thyroid. 2002;12:613–7.
Chehoud C, Dryga A, Hwang Y, Nagy-Szakal D, Hollister EB, Luna RA, et al. Transfer of viral communities between human individuals during fecal microbiota transplantation. MBio. 2016;7:e00322.
Morduchowicz G, Pitlik SD, Huminer D, Alkan M, Drucker M, Rosenfeld JB, et al. Transfusion reactions due to bacterial contamination of blood and blood products. Rev Infect Dis. 1991;13:307–14.
Bihl F, Castelli D, Marincola F, Dodd RY, Brander C. Transfusion-transmitted infections. J Transl Med. 2007;5:1–11.
Angelberger S, Reinisch W, Makristathis A, Lichtenberger C, Dejaco C, Papay P, et al. Temporal bacterial community dynamics vary among ulcerative colitis patients after fecal microbiota transplantation. Am J Gastroenterol. 2013;108:1620–30.
Fuentes S, van Nood E, Tims S, Heikamp-de Jong I, Ter Braak CJ, Keller JJ, et al. Reset of a critically disturbed microbial ecosystem: faecal transplant in recurrent Clostridium difficile infection. ISME J. 2014;8:1621–33.
Weingarden A, González A, Vázquez-Baeza Y, Weiss S, Humphry G, Berg-Lyons D, et al. Dynamic changes in short- and long-term bacterial composition following fecal microbiota transplantation for recurrent Clostridium difficile infection. Microbiome. 2015;3:10.
Kumar R, Maynard CL, Eipers P, Goldsmith KT, Ptacek T, Grubbs JA, et al. Colonization potential to reconstitute a microbe community in patients detected early after fecal microbe transplant for recurrent C. difficile. BMC Microbiol. 2016;16:5. doi:10.1186/s12866-015-0622-2.
Janda JM, Abbott SL. 16S rRNA gene sequencing for bacterial identification in the diagnostic laboratory: pluses, perils, and pitfalls. J Clin Microbiol. 2007;45:2761–4.
Methé A, et al. A framework for human microbiome research. Nature. 2012;486:215–21.
Voigt AY, Costea PI, Kultima JR, Li SS, Zeller G, Sunagawa S, et al. Temporal and technical variability of human gut metagenomes. Genome Biol. 2015;16:1–12.
Li SS, Zhu A, Benes V, Costea PI, Hercog R, Hildebrand F, et al. Durable coexistence of donor and recipient strains after fecal microbiota transplantation. Science. 2016;352:586–9.
Tyson GW, Chapman J, Hugenholtz P, Allen EE, Ram RJ, Richardson PM, et al. Community structure and metabolism through reconstruction of microbial genomes from the environment. Nature. 2004;428:37–43.
Venter JC, Remington K, Heidelberg JF, Halpern AL, Rusch D, Eisen JA, et al. Environmental genome shotgun sequencing of the Sargasso Sea. Science. 2004;304:66. LP-74.
Eren AM, Vineis JH, Morrison HG, Sogin ML. A filtering method to generate high quality short reads using illumina paired-end technology. PLoS One. 2013;8:e66643.
Minoche AE, Dohm JC, Himmelbauer H. Evaluation of genomic high-throughput sequencing data generated on Illumina HiSeq and genome analyzer systems. Genome Biol. 2011;12:R112.
Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.
Kim D, Song L, Breitwieser FP, Salzberg SL. Centrifuge: rapid and sensitive classification of metagenomic sequences. bioRxiv. Genome Res. 2016;26(12):1721–9.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Meth. 2012;9:357–9.
Eren AM, Esen ÖC, Quince C, Vineis JH, Morrison HG, Sogin ML, et al. Anvi’o: an advanced analysis and visualization platform for ‘omics data. PeerJ. 2015;3:e1319.
Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:1–11.
Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39:W29–31.
Campbell JH, O’Donoghue P, Campbell AG, Schwientek P, Sczyrba A, Woyke T, et al. UGA is an additional glycine codon in uncultured SR1 bacteria from the human microbiota. Proc Natl Acad Sci. 2013;110:5540–5.
Rinke C, Schwientek P, Sczyrba A, Ivanova NN, Anderson IJ, Cheng J-F, et al. Insights into the phylogeny and coding potential of microbial dark matter. Nature. 2013;499:431–7.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Alneberg J, Bjarnason BS, de Bruijn I, Schirmer M, Quick J, Ijaz UZ, et al. Binning metagenomic contigs by coverage and composition. Nat Meth. 2014;11:1144–6.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.
Darling AE, Jospin G, Lowe E, Matsen IV FA, Bik HM, Eisen JA. PhyloSift: phylogenetic analysis of genomes and metagenomes. Moustafa A, editor. PeerJ. 2014;2:e243
Birney E, Andrews TD, Bevan P, Caccamo M, Chen Y, Clarke L, et al. An overview of Ensembl. Genome Res. 2004;14:925–8.
Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26:1641–50.
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:1–15.
Pasolli E, Truong DT, Malik F, Waldron L, Segata N. Machine learning meta-analysis of large metagenomic datasets: tools and biological insights. PLoS Comput Biol. 2016;12:1–26.
Truong DT, Franzosa EA, Tickle TL, Scholz M, Weingart G, Pasolli E, et al. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat Meth. 2015;12:902–3.
Oksanen JF, Blanchet G, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. R package version 2.4-0. https://CRAN.R-project.org/package=vegan. 2016. Accessed 1 May 2017.
Parks DH, Beiko RG. Identifying biologically relevant differences between metagenomic communities. Bioinformatics. 2010;26:715–21.
Shahinas D, Silverman M, Sittler T, Chiu C, Kim P, Allen-Vercoe E, et al. Toward an understanding of changes in diversity associated with fecal microbiome transplantation based on 16S rRNA gene deep sequencing. MBio. 2012;3:1–10.
Hamilton MJ, Weingarden AR, Unno T, Khoruts A, Sadowsky MJ. High-throughput DNA sequence analysis reveals stable engraftment of gut microbiota following transplantation of previously frozen fecal bacteria. Gut Microbes. 2013;4:125–35.
Seekatz AM, Aas J, Gessert CE, Rubin TA, Saman DM, Bakken JS, et al. Recovery of the gut microbiome following fecal microbiota transplantation. MBio. 2014;5(3):e00893–14. doi:10.1128/mBio.00893-14.
Brandt LJ, Aroniadis OC, Mellow M, Kanatzar A, Kelly C, Park T, et al. Long-term follow-up of colonoscopic fecal microbiota transplant for recurrent Clostridium difficile infection. Am J Gastroenterol. 2012;107:1079–87.
van Nood E, Vrieze A, Nieuwdorp M, Fuentes S, Zoetendal EG, de Vos WM, et al. Duodenal infusion of donor feces for recurrent Clostridium difficile. N Engl J Med. 2013;368:407–15.
Vineis JH, Ringus DL, Morrison HG, Delmont TO, Dalal S, Raffals LH. Patient-specific bacteroides genome variants in pouchitis. MBio. 2016;7:1–11.
Wang S, Xu M, Wang W, Cao X, Piao M, Khan S, et al. Systematic review: adverse events of fecal microbiota transplantation. PLoS One. 2016;11:e0161174.
Bojanova DP, Bordenstein SR. Fecal transplants: what is being transferred? PLoS Biol. 2016;14:e1002503.
Khanna S, Pardi DS, Kelly CR, Kraft CS, Dhere T, Henn MR, et al. A novel microbiome therapeutic increases gut microbial diversity and prevents recurrent Clostridium difficile infection. J Infect Dis. 2016;214:173–81.
Ratner M. Seres’s pioneering microbiome drug fails mid-stage trial. Nat Biotech. 2016;34:1004–5.
Nayfach S, Rodriguez-Mueller B, Garud N, Pollard KS. An integrated metagenomics pipeline for strain profiling reveals novel patterns of bacterial transmission and biogeography. Genome Res. 2016;26:1612–25.
We thank Rich Fox for his help with computer systems and gratefully acknowledge Sarah Goeppinger, Ruben Colman, Dylan Rodriquez, Yunwei Wang, and Joseph Vineis for their technical and administrative support, Sebastian Schmidt and Simone Li for their comments on our pre-print, and our anonymous reviewers for their efforts to improve our study.
AME was supported by the Frank R. Lillie Research Innovation Award and startup funds from the University of Chicago. This project was supported by the Mutchnik Family Charitable Fund and the University of Chicago Gastro-Intestinal Research Foundation.
Availability of data and materials
We made publicly available the anvi’o profiles to reproduce all findings and to interactively visualize and further investigate all MAGs in this study (doi:10.6084/m9.figshare.4792621), the FASTA files and distribution statistics for each MAG across samples (doi:10.6084/m9.figshare.4792627), and individual anvi’o profiles for the detection of the donor MAGs across HMP metagenomes (doi:10.6084/m9.figshare.4793761). See http://merenlab.org/data/ for visualization and analysis instructions. Raw metagenomic data are also stored at the NCBI Sequence Read Archive under the accession number SRP093449.
SAK, NAH, DAA, and DTR designed the study and collected and processed the patient samples. STML, TOD, HGM, and AME generated, processed, and analyzed the sequencing data. AS and OCE contributed to the bioinformatics analyses. STML, TOD, and AME wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study was reviewed and approved by the University of Chicago Ethics Committee and by the University of Chicago Institutional Review Board (IRB 132-0212). Written and informed consent was obtained for all participants.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Phylogenomic relationships between our MAGs and a collection of 1758 reference genomes. (XLSX 775 kb)
Additional information for samples, metagenomic bins, and metagenome-assembled genomes in our study. (XLSX 180 kb)
The distribution of all 444 donor bins across samples. Yellow markings represent the 92 MAGs. (PNG 7138 kb)
Phylogenomic analysis of 92 MAGs in the context of 1758 gold standard genomes. Squares and triangles next to MAGs indicate the agreements (green) and disagreements (red) between the phylogenomic analysis and CheckM results. The blank ones indicate that we did not assign any taxonomy for a given MAG at the specified level. (PNG 1729 kb)
Detection of the donor MAGs across 151 HMP metagenomes. (PNG 983 kb)
The canonical correspondence analysis (CCA) of 500 functions and 107 sub-systems. (XLSX 240 kb)
Non-metric multidimensional scaling based on microbial genus-level taxonomy of 4 donors, 6 recipients, and 151 HMP metagenomes after the removal of short reads matching to donor bins from recipient metagenomes. (PNG 2235 kb)
About this article
Cite this article
Lee, S.T.M., Kahn, S.A., Delmont, T.O. et al. Tracking microbial colonization in fecal microbiota transplantation experiments via genome-resolved metagenomics. Microbiome 5, 50 (2017). https://doi.org/10.1186/s40168-017-0270-x