Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Host-Microbe Biology

Divergent Relationships between Fecal Microbiota and Metabolome following Distinct Antibiotic-Induced Disruptions

Jocelyn M. Choo, Tokuwa Kanno, Nur Masirah Mohd Zain, Lex E. X. Leong, Guy C. J. Abell, Julie E. Keeble, Kenneth D. Bruce, A. James Mason, Geraint B. Rogers
Rosa Krajmalnik-Brown, Editor
Jocelyn M. Choo
aInfection and Immunity Theme, South Australia Health and Medical Research Institute, North Terrace, Adelaide, SA, Australia
cSchool of Medicine, Flinders University, Bedford Park, Adelaide, SA, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tokuwa Kanno
bKing’s College London, Institute of Pharmaceutical Science, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nur Masirah Mohd Zain
bKing’s College London, Institute of Pharmaceutical Science, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lex E. X. Leong
aInfection and Immunity Theme, South Australia Health and Medical Research Institute, North Terrace, Adelaide, SA, Australia
cSchool of Medicine, Flinders University, Bedford Park, Adelaide, SA, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Guy C. J. Abell
cSchool of Medicine, Flinders University, Bedford Park, Adelaide, SA, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julie E. Keeble
bKing’s College London, Institute of Pharmaceutical Science, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kenneth D. Bruce
bKing’s College London, Institute of Pharmaceutical Science, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
A. James Mason
bKing’s College London, Institute of Pharmaceutical Science, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Geraint B. Rogers
aInfection and Immunity Theme, South Australia Health and Medical Research Institute, North Terrace, Adelaide, SA, Australia
cSchool of Medicine, Flinders University, Bedford Park, Adelaide, SA, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rosa Krajmalnik-Brown
Arizona State University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphere.00005-17
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

The intestinal microbiome plays an essential role in regulating many aspects of host physiology, and its disruption through antibiotic exposure has been implicated in the development of a range of serious pathologies. The complex metabolic relationships that exist between members of the intestinal microbiota and the potential redundancy in functional pathways mean that an integrative analysis of changes in both structure and function are needed to understand the impact of antibiotic exposure. We used a combination of next-generation sequencing and nuclear magnetic resonance (NMR) metabolomics to characterize the effects of two clinically important antibiotic treatments, ciprofloxacin and vancomycin-imipenem, on the intestinal microbiomes of female C57BL/6 mice. This assessment was performed longitudinally and encompassed both antibiotic challenge and subsequent microbiome reestablishment. Both antibiotic treatments significantly altered the microbiota and metabolite compositions of fecal pellets during challenge and recovery. Spearman’s correlation analysis of microbiota and NMR data revealed that, while some metabolites could be correlated with individual operational taxonomic units (OTUs), frequently multiple OTUs were associated with a significant change in a given metabolite. Furthermore, one metabolite, arginine, can be associated with increases/decreases in different sets of OTUs under differing conditions. Taken together, these findings indicate that reliance on shifts in one data set alone will generate an incomplete picture of the functional effect of antibiotic intervention. A full mechanistic understanding will require knowledge of the baseline microbiota composition, combined with both a comparison and an integration of microbiota, metabolomics, and phenotypic data.

IMPORTANCE Despite the fundamental importance of antibiotic therapies to human health, their functional impact on the intestinal microbiome and its subsequent ability to recover are poorly understood. Much research in this area has focused on changes in microbiota composition, despite the interdependency and overlapping functions of many members of the microbial community. These relationships make prediction of the functional impact of microbiota-level changes difficult, while analyses based on the metabolome alone provide relatively little insight into the taxon-level changes that underpin changes in metabolite levels. Here, we used combined microbiota and metabolome profiling to characterize changes associated with clinically important antibiotic combinations with distinct effects on the gut. Correlation analysis of changes in the metabolome and microbiota indicate that a combined approach will be essential for a mechanistic understanding of the functional impact of distinct antibiotic classes.

INTRODUCTION

The intestinal microbiome plays an important role in regulating many aspects of host physiology, including glucose and fat metabolism (1–3), systemic immunity (4, 5), and central nervous system function (6, 7). Antibiotic-induced disruption of the intestinal microbiota has been known for some time to be associated with infection, such as by Clostridium difficile (8). For example, cefoperazone administration in mice results in the development of an intestinal metabolome that favors C. difficile germination and growth (9). In addition, antibiotic-induced perturbation of gut microbiota is being increasingly associated with both the loss of beneficial functions and the gain of microbiota functions that are proinflammatory (10).

The clear importance of the intestinal microbiota to host health has also led to growing concerns about the wider impact of antibiotic interventions on chronic conditions, particularly during human development (11–16). While establishing links between antibiotic exposure and outcome is challenging, large-cohort studies have concluded that the use of certain antibiotics, including vancomycin, cephalosporins, penicillins, and particularly macrolides, can be associated with weight gain (15–18). For example, weight gain in mice has been associated with an increase in Lactobacillus reuteri and a decrease in Escherichia coli cells as a result of vancomycin administration (19). Changes in short-chain fatty acids (SCFAs) resulting from exposure to vancomycin suggest a functional consequence of antibiotic exposure. A change in a given metabolite therefore may prove to be a more direct link to a particular pathophysiological change than loss or gain of individual taxa (20–22), especially in light of functional redundancy (19, 20, 23). SCFAs are immunomodulatory (21, 22), can interact with nerve cells to stimulate the sympathetic and autonomic nervous system (24, 25), regulate colonic energy metabolism and autophagy (26), and have a range of other important functions (1, 27, 28). Increasing the focus on the metabolome also reduces concern over the impact of variety in baseline gastrointestinal microbiota, which as we have recently shown, differs even in genetically identical mice within the same production facility (29). Studying antibiotic-induced changes by combining the microbiota and metabolome therefore provides a promising route for developing mechanistic insight into their impacts on the microbial community and, in the longer term, on host outcome.

Here, we investigate the relationship between the antibiotic-induced disruption and reestablishment of the intestinal microbiota and the corresponding changes in the fecal metabolome in genetically identical mice. We used two distinct antibiotic treatments, ciprofloxacin and a vancomycin-imipenem combination, that importantly have been shown to differ markedly in their therapeutic impacts on murine microbiota (30), with ciprofloxacin having minimal effects on the anaerobic gut microbiota and vancomycin-imipenem reducing intestinal anaerobes with a broad-spectrum activity (31). Our first aim therefore was to characterize how two substantially different antibiotic regimens have an impact on the same mouse gut microbiota and metabolome. By so doing, we also aimed to determine the additional value resulting from integrating these two approaches in terms of understanding the impact of antibiotic challenge on the murine gut microbial community.

RESULTS

Antibiotics alter microbial alpha diversity.The fecal microbiota and metabolome were assessed immediately prior to antibiotic treatment (baseline, or time 1 [T1]), after 14 days of antibiotic treatment (during the course of antibiotic treatment, T2), and 9 days following the termination of antibiotics (recovery, T3). Fecal microbiota alpha diversity was assessed across each time point for the control and antibiotic-treated groups. The taxon richness (Taxa_S), evenness (Simpson index), and diversity (Shannon index) of murine fecal microbiota are shown in Fig. 1A to C, respectively. Ciprofloxacin treatment resulted in a significant reduction in taxon richness (P = 0.008, Wilcoxon test) but no change in microbiota diversity and evenness. Vancomycin-imipenem treatment resulted in a significant decrease in taxon richness (P = 0.016), evenness (P = 0.016), and diversity (P = 0.016) (Fig. 1).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Alpha diversity analysis of fecal microbiota of the control, ciprofloxacin, and vancomycin-imipenem groups. The fecal microbial community of mice was analyzed immediately before treatment (baseline, T1) and after 14 days of antibiotic treatment (T2), as well as at 9 days after antibiotic treatment (T3), to analyze for recovery from antibiotics. Statistical significance between the T1, T2, and T3 time points were analyzed using the Wilcoxon test at a significance level of 0.05.

Antibiotic treatment resulted in shifts in the fecal microbial composition and structure.The fecal-microbiota compositions of the control and antibiotic groups were visualized by nonmetric multidimensional scaling (NMDS) based on the Bray-Curtis dissimilarity distances of square root-transformed relative abundances (Fig. 2A; see also the compositional data shown in Fig. S1 in the supplemental material). Permutational multivariate analysis of variance (PERMANOVA) indicated that the gut microbiota compositions between the groups did not differ significantly prior to the antibiotic treatment period (Table S1). Analysis of the distance to a group’s centroid (permutational analysis of the homogeneity of group dispersions [PERMDISP]) also revealed no significant difference in the levels of the homogeneity of dispersions of the microbial community between the control group and the ciprofloxacin {t = 1.245, permutation test P value [P(perm)] = 0.387} or vancomycin-imipenem [t = 1.451, P(perm) = 0.328] group at baseline (Fig. 2B).

FIG S1

Bacterial OTU plots for the control, ciprofloxacin, and vancomycin-imipenem groups at baseline (T1) (A), at the completion of 14 days of antibiotic treatment (T2) (B), and 9 days after cessation of antibiotic treatment (T3) (C). Download FIG S1, PDF file, 0.2 MB.
Copyright © 2017 Choo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

TABLE S1

PERMANOVA of the microbiota compositions among the control and antibiotic groups accompanied by pairwise comparisons at baseline (T1), at the end of 14 days of treatment (T2), and 9 days after treatment (T3) and the corresponding Q2 for pairwise OPLS-DA comparisons of metabolomes. Download TABLE S1, XLSX file, 0.01 MB.
Copyright © 2017 Choo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Shifts in the microbial community compositions of the control and antibiotic groups. (A) NMDS plot based on Bray-Curtis distances of the fecal microbiota. Each data point represents a mouse fecal sample from the control, ciprofloxacin, or vancomycin-imipenem group at each time point, with the corresponding labels indicated on the right. (B) The distance from the centroid within each group was analyzed using PERMDISP to determine the homogeneity of dispersion. Statistical significance between the control and antibiotic groups for each time point was analyzed based on the P value of the permutation tests at a significance level of 0.05. C, control; CIP, ciprofloxacin; VI, vancomycin-imipenem.

Within 14 days of antibiotic treatment, the fecal microbial compositions between the control and antibiotic groups were significantly different [P(perm) = 0.0001, square root estimated components of variation (ECV) = 41.73, 9,951 permutations] (Table S1). Pairwise comparison indicated that ciprofloxacin significantly altered the microbiota composition from that in the control group [t = 3.82, P(perm) = 0.0003, 5,069 permutations] (Table S1). Antibiotic-associated effects on the homogeneity of microbiota community dispersion were analyzed using PERMDISP (Fig. 2B), and linear discriminant analysis effect size (LEfSe) analysis was performed to determine bacterial taxa that significantly differed between the control and antibiotic groups (Fig. 3). The distance to the group’s centroid was significantly reduced in the ciprofloxacin group compared to that for controls [t = 3.466, P(perm) = 0.003] (Fig. 2B), suggesting that microbiota changes resulted from antibiotic selective pressure. LEfSe comparisons indicated that ciprofloxacin resulted in significant decreases in several taxa, including Odoribacter, Alistipes, Streptococcus, Lactobacillus, Clostridium, Turicibacter, and RC9 as well as the Prevotellaceae (uncultured) and RF9 families (Fig. 3A), which were completely depleted by ciprofloxacin administration compared to the taxon levels of the control group (Fig. S2). Additionally, increased relative abundances of Bacteroides, Marvinbryantia, and Coprococcus organisms were observed in the ciprofloxacin group.

FIG S2

Heatmap analysis of bacterial taxa that significantly differed between the control and ciprofloxacin or vancomycin-imipenem groups at the end of antibiotic treatment and 9 days after cessation of antibiotics, based on LEfSe analysis (LDA score, >3.0). The dendrogram generated represents hierarchical clustering of bacterial taxa based on Bray-Curtis distances and Ward’s method of clustering. The abundance of each bacterial taxon is depicted based on the square root OTU counts. Download FIG S2, PDF file, 1 MB.
Copyright © 2017 Choo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

LEfSe comparison analysis between the control and ciprofloxacin or vancomycin-imipenem groups at the end of antibiotic treatment (A or B, respectively) or 9 days after cessation of antibiotic treatment (C or D, respectively). The red or green shading depicts bacterial taxa that were significantly higher in either the control, ciprofloxacin, or vancomycin-imipenem group, as indicated. The yellow circles on the cladogram depict bacterial taxa that were not significantly changed. Selection of discriminative taxa between the control and antibiotic groups were based on an LDA score cutoff of 3.0, and differences in the relative abundances of taxa were statistically determined based on a Mann-Whitney test at a significance level of 0.05.

Vancomycin-imipenem also resulted in a significant difference in the fecal microbiota compositions [t = 10.43, P(perm) = 0.0003, 5,043 permutations] and significantly reduced the distance to the group’s centroid compared to that of the control group [t = 8.896, P(perm) = 0.0004] (Table S1 and Fig. 2B, respectively). Comparison of the t statistic values indicated that these changes were greater than those observed with ciprofloxacin. LEfSe analysis comparing the control and vancomycin-imipenem groups indicated significant depletion of large numbers of taxa, including in members of the Bacteroidetes phylum (Odoribacter, Alistipes, Bacteroides, Parabacteroides, RC9, S24-7), Firmicutes phylum (Roseburia, Candidatus arthromitus, Clostridium, and Turicibacter), and Enterorhabdus, as well as members of the Lachnospiraceae, Ruminococcaceae, Prevotellaceae, and RF9 families (Fig. 3B). These changes were accompanied by increases in the relative abundance of Proteobacteria, including Enterobacter, Escherichia-Shigella, Citrobacter, Achromobacter, and Salmonella, as well as the genera Streptococcus, Lactobacillus, and Planomicrobium.

Agent-specific recovery of baseline fecal microbiota following antibiotic treatment.The degree to which the fecal microbial community had recovered by 9 days following the completion of antibiotic treatment was assessed. Following the recovery period, levels of taxon richness in the ciprofloxacin group remained unchanged; however, microbiota evenness and diversity were significantly reduced compared to levels at antibiotic cessation (P = 0.016 in both cases, Wilcoxon test) and to baseline levels (P = 0.016 and P = 0.031, respectively) (Fig. 1). In addition, a PERMANOVA pairwise comparison indicated that the microbiota composition of the ciprofloxacin group remained significantly different from that of the controls [t = 3.04, P(perm) = 0.0002, 5,100 permutations] (Table S1). This difference from the microbiota composition at baseline was also significant [t = 4.29, P(perm) = 0.0003, 5,023 permutations] (Table S2). LEfSe analysis indicated a significantly increased relative abundance of Bacteroides organisms and a depletion of Alistipes spp., Odoribacter spp., and the RC9 gut group in the ciprofloxacin group (Fig. 3C), which failed to be reestablished in the gut microbiota (Fig. S2).

TABLE S2

Pairwise comparisons of fecal microbiota compositions and metabolomes at different time points, including at baseline (T1), at the end of 14 days of treatment (T2), and 9 days after treatment (T3) for the control, ciprofloxacin, and vancomycin-imipenem groups. Download TABLE S2, XLSX file, 0.01 MB.
Copyright © 2017 Choo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

In the vancomycin-imipenem group, the levels of microbial richness, evenness, and diversity significantly increased 9 days after antibiotic cessation compared to the levels at antibiotic cessation (P = 0.008 in all cases, Wilcoxon test) but did not reach baseline levels (P = 0.031 in all cases) (Fig. 1). The fecal microbiota composition of the vancomycin-imipenem group also remained significantly different from that of the control group [t = 3.45, P(perm) = 0.0001, 5,047 permutations] or from the baseline composition [t = 6.42, P(perm) = 0.0003, 5,045 permutations] (Table S1 and S2, respectively). The extent of the differences between the microbiota compositions in the vancomycin-imipenem group and the control (66.02% mean similarity) was greater than that observed with ciprofloxacin (70.88% mean similarity) based on Bray-Curtis distances, suggesting that recovery of the microbiota was slower with vancomycin-imipenem administration.

In keeping with the posttreatment effects observed with ciprofloxacin, LEfSe analysis revealed that the taxa Alistipes, Odoribacter, and the RC9 gut group, which were depleted during vancomycin-imipenem treatment, were not restored during the recovery period (Fig. 3D). The relative abundances of several additional taxa, including Family XII and the Ruminococcaceae, Prevotellaceae, and RF9 families, were also significantly reduced after antibiotic treatment in the vancomycin-imipenem group. The relative abundances of Escherichia-Shigella organisms remained elevated, while the relative abundances of Bacteroides, Enterorhabdus, and Enterococcus organisms increased posttreatment.

Antibiotic-induced dysbiosis results in antibiotic-specific alteration of the fecal metabolome.The influence of antibiotic-induced microbiota alterations on the fecal metabolome composition was assessed by NMR spectroscopy. Supervised multivariate analysis based on a cross-validated OPLS-DA (orthogonal projections to latent structures discriminant analysis) model (leave one out) demonstrated that the metabolic profile of the ciprofloxacin group was altered from that of the control group (predictive ability [Q2] = 0.476). A volcano plot comparing the ciprofloxacin group to the control group indicated that the levels of amino acids, such as valine, leucine, isoleucine, and phenylalanine, as well as α-aminobutyric acid, were significantly increased in the ciprofloxacin-treated group, while levels of the sugar glycerol decreased (Fig. 4A).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Volcano plot analysis of fecal metabolites altered in the ciprofloxacin or vancomycin-imipenem groups at the end of treatment (A) or 9 days after cessation of antibiotics (B). The x and y axes of the volcano plot represent the log2-fold changes (against the control group) and the corresponding log FDR-adjusted P value of all metabolites, respectively. The vertical and horizontal lines separate metabolites that had a 2-fold change and an FDR-adjusted P value at 0.05. AABA, α-aminobutyric acid; GABA, γ-aminobutyric acid.

The administration of vancomycin-imipenem resulted in a metabolome profile distinct from that the control group (Q2 = 0.981), and these differences were larger than those observed in the ciprofloxacin group (Fig. 4A). Several metabolites, including the amino acids alanine, methionine, tyrosine, glutamine, arginine, and asparagine, as well as the organic acids citrate, α-aminobutyric acid, and propionate, were found to be lower in the vancomycin-imipenem group than in the controls. The nucleoside uridine was also significantly lower in the treated group than in the controls. Increased levels of sucrose, sarcosine, arabinitol, and gluconate were observed.

Agent-specific reestablishment of the fecal metabolome following antibiotic dysbiosis.The OPLS-DA model for the recovery of the fecal metabolome indicated a moderate shift in the metabolite profile of the ciprofloxacin group relative to that of the controls (Q2 = 0.475), although the corresponding volcano plot analysis indicated that no individual metabolites differed significantly between the two groups (Fig. 4B).

In contrast, the OPLS-DA model describing the metabolic patterns of the vancomycin-imipenem group relative to those of the control group was very strong (Q2 = 0.872), although only four metabolites were significantly altered when considered individually (Fig. 4B). The vancomycin-imipenem group was again found to have lower levels of methionine but also higher levels of γ-aminobutyric acid, arginine, and lysine than the control group.

The relationship between altered microbiota composition and changes in the fecal metabolome.The extensive changes in the fecal microbiota and metabolome following exposure to vancomycin-imipenem suggested that these responses can be used to assess microbiota-metabolome associations. The functional correlation between alterations in the microbiota and metabolites was assessed using Spearman’s correlation, based on 11 bacterial taxa (up to a 70% contribution based on a similarity percentage [SIMPER] analysis) and 14 NMR metabolites (fold change significance at a P of <0.05) that contributed substantially to the differences between the control and vancomycin-imipenem groups at the end of the treatment (Fig. 5B). The correlation analysis indicated that all 14 metabolites could be correlated with changes in bacterial taxa. The increase in metabolites such as arabinitol and sucrose were positively correlated with increased relative abundances of Enterobacter, Escherichia-Shigella, and Lactobacillus organisms but negatively correlated with Ruminococcaceae, Lachnospiraceae, and Enterorhabdus. Increased levels of sarcosine also resulted in similar microbiota-metabolite relationships, although no correlations were observed with the taxa Lactobacillus, Enterobacter, and Roseburia. Tyrosine, methionine, citrate, asparagine, alanine, and γ-aminobutyric acid shared similar correlation patterns with each other (Fig. 5B). Fewer correlations could be found for the ciprofloxacin-treated group, with only four out of six metabolites correlated with OTU-level changes during treatment (Fig. 5A) and only arginine correlating with OTU during the recovery from vancomycin-imipenem challenge (Fig. 5D).

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Spearman’s correlation analysis of microbiota and NMR spectra of the control and ciprofloxacin groups (A, C) or the vancomycin-imipenem group (B, D) after 14 days of antibiotic treatment and 9 days after cessation of antibiotics, respectively. Bacterial taxa and metabolites in the correlation matrix were sorted based on the Euclidean distance and Ward’s method of hierarchical clustering. Significant microbiota-metabolite correlations were determined based on an r of less than −0.75 or more than 0.75 and an FDR adjusted P value of less than 0.01.

The potential for changes in the intestinal microbiota to explain associated changes in the metabolome were investigated using PICRUSt. This approach was used to predict shifts in the prevalences of genes that encode metabolic traits, based on changes in the relative abundances of detected taxa. Metabolic pathways predicted to alter significantly as a result of changes in microbiota composition are shown for antibiotic treatment and at recovery periods for ciprofloxacin and vancomycin-imipenem (Fig. S3A and B and S4A and B, respectively). While such projections do not take into consideration cross-metabolism between taxa or levels of gene transcription, these results were consistent in many cases with observed microbe-metabolite correlations. For example, treatment with vancomycin-imipenem resulted in a decrease in the abundance of metabolites involved in the alanine, aspartate, and glutamate metabolic pathways. In keeping with these observed changes, predicted shifts in the metagenome included a decrease in the carriage of genes that encode related enzymes. Decreases in aspartate aminotransferase (an enzyme involved in the conversion of l-aspartate to oxaloacetate and glutamine), aspartate 4-decarboxylase (which directly converts l-aspartate to l-alanine), and alanine dehydrogenase (which is involved in converting pyruvate to l-alanine) were predicted, while genes encoding alanine-synthesizing transaminase (involved in the interconversion between pyruvate and l-alanine) increased in prevalence. For enzymes that were involved in the conversion of succinate semialdehyde to either glutamate or succinate, significant reductions in 4-aminobutyrate aminotransferase/(S)-3-amino-2-methylpropionate transaminase (involved in the conversion to glutamate) were observed, while an increase in succinate-semialdehyde dehydrogenase (NADP+) (involved in the conversion to succinate) was observed. Decreases were seen in 1-pyrroline-5-carboxylate dehydrogenase (involved in the production of l-glutamate) and in glutamate synthase and glutaminase (involved in the direct conversion of glutamine to glutamate), and an increase in the prevalence of alanine-synthesizing transaminase was also predicted.

FIG S3

Predicted changes in bacterial contributions to the fecal metabolome as a result of ciprofloxacin treatment (A) and during recovery after ciprofloxacin treatment (B), estimated based on changes in the observed taxon relative abundances using PICRUSt. Mean proportions of pathways assigned to control and treatment groups (percentages) and Welch's t test with Benjamini-Hochberg false-discovery rate-corrected significance are shown. Download FIG S3, PDF file, 0.2 MB.
Copyright © 2017 Choo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

FIG S4

Predicted changes in bacterial contributions to the fecal metabolome as a result of vancomycin-imipenem treatment (A) and during recovery after vancomycin-imipenem treatment (B), estimated based on changes in the observed taxon relative abundances using PICRUSt. Mean proportions of pathways assigned to control and treatment groups (percentages) and Welch's t test with Benjamini-Hochberg false-discovery rate-corrected significance are shown. Download FIG S4, PDF file, 0.3 MB.
Copyright © 2017 Choo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

PICRUSt-derived predictions of changes in gene prevalence also included those involved in arginine biosynthesis and tyrosine metabolism, pathways which showed significant changes in NMR-derived metabolome profiles. For example, there was a decrease in the prevalence of acetylglutamate/acetylaminoadipate kinase, N-acetyl-gamma-glutamyl-phosphate reductase, and acetylornithine aminotransferase, each of which is involved in the conversion from N-acetylglutamate to N-acetylornithine (a substrate for arginine synthesis). Finally, an increase in genes encoding 4-hydroxyphenylpyruvate dioxygenase, an enzyme that converts 4-hydroxyphenylpyruvate (a substrate for tyrosine metabolism) to homogentisate, were consistent with the observed decrease in tyrosine levels. Furthermore, in keeping with the observed increase in the abundance of arabinitol in the vancomycin-imipenem-treated mice, there was a significant decrease in the predicted prevalence of genes encoding xylulokinases, which are involved in its catabolism.

Changes in microbiota and metabolome composition were less marked in mice treated with ciprofloxacin. However, a predicted decrease in the carriage of acyl coenzyme A (acyl-CoA) dehydrogenase was consistent with the observed increase in the abundance of valine, leucine, and isoleucine.

DISCUSSION

The gut microbiota is critical for human health (32–34), with a number of essential metabolites derived exclusively through the activities of these intestinal microbes (30, 35, 36). The gut microbiota also produces of a wide range of compounds that modulate host physiology (37), including immune regulation, metabolic control (32, 38), central nervous system function (39), and xenobiotic metabolism (40). Analysis of antibiotic-induced disturbances in the gut microbiota and its corresponding metabolome can therefore provide insight into both acute and chronic effects of antibiotics on the host and may yield a functional understanding of the development of associated conditions (41). A number of previous studies have used C57BL/6 mice to investigate links between antibiotic-induced changes in intestinal microbiota composition and host physiology (42–45), while assessment of metabolomic changes have been limited to other genetic backgrounds (10, 46).

The present study highlights a number of functional effects on the mouse gut microbiome and metabolome which warrant further investigation and are potential starting points for a mechanistic understanding of the role of various antibiotic regimens in, e.g., gut inflammation and/or obesity. More importantly, the present study reveals substantial variation in dynamic relationships between individual taxa, groups of taxa, and metabolites at different stages of antibiotic intervention. The following examples indicate that while some interesting information may be determined by examining either microbiota or the metabolome, the most informative mechanistic investigations will involve a combined assessment of the relationships between microbiota and the metabolome for a given phenotypic trait.

In a first example, vancomycin-imipenem treatment resulted in substantial changes in microbiota composition, including the depletion of members of the Ruminococcaceae and Lachnospiraceae families and increases in the relative abundances of Enterobacteriaceae and Lactobacillus spp., with concomitant changes in the levels of a large number of metabolites. The potential for microbiota-level changes to influence disease risk by altering the composition of the metabolome are well illustrated by such changes. For example, arabinitol is poorly absorbed in the intestine, and its removal is based on its conversion to pentose sugars by taxa, including Lachnospiraceae and Ruminococcaceae (47). Further, the prevalence of genes encoding arabinitol-catabolic enzymes was predicted to be significantly reduced based on changes in relative taxon abundances. We observed a negative correlation between the relative abundances of Lachnospiraceae and Ruminococcaceae and fecal arabinitol levels. Increased levels of arabinitol and sugars, such as sucrose, that were also increased during vancomycin-imipenem treatment have been associated with susceptibility to C. difficile infection in mice and might act as a growth substrate (9). Further, the increased relative abundance of Enterobacteriaceae that accompanied decreases in the abundances of Lachnospiraceae and Ruminococcaceae are strongly associated with gut inflammation (48). These changes correlated with decreased glutamine levels, which has anti-inflammatory effects potentially through its role in the maintenance of tissue permeability and its inhibitory action on NF-kB activation and p38 mitogen-activated protein kinase (MAPK) pathways based on animal and human intestinal studies (49, 50).

In contrast, the functional relationships between bacterial taxa and, e.g., arginine are more complex. The reestablishment of microbiota/metabolome composition 9 days after cessation of vancomycin-imipenem treatment was characterized by significant increases in arginine. These increases were correlated with increases in Enterobacter spp. and decreases in Alistipes spp. Again, the prevalence of genes encoding enzymes in the arginine biosynthesis pathway were predicted to be significantly increased, based on taxon relative abundances. Arginine plays an important role as a precursor for immune-modulatory compounds (51–53), and the functional significance of antibiotic-induced metabolite changes in arginine may therefore be explored in future studies. Changes in arginine were also observed, however, immediately at the end of the 14-day treatment with vancomycin-imipenem. Here a modest but significant reduction in arginine was instead correlated with depletion of Ruminococcaceae and Bacteroides spp. and increases in Escherichia-Shigella spp., and there was no correlation with either Enterobacter or Alistipes spp. Arginine levels in the gut may therefore respond to changes in the prevalence of multiple OTUs, and there may be considerable functional redundancy. As such, given that we have previously observed functional divergence in gastrointestinal microbiotas even in genetically identical mice that originate in the same production facility (29), it may be inaccurate to identify individual OTUs as responsible for a particular functional impact without considering the baseline microbiota composition or without supporting metabolomic data.

In a second example, the levels of severity of antibiotic impact can be substantially different when assessed by next-generation sequencing and NMR metabolomics. In keeping with studies of humans (54–56), the administration of ciprofloxacin also substantially reduced microbial richness and evenness. Here, normal microbiota composition was not restored in the mice within the 9 days after antibiotic treatment, a process that can take up to a month in humans (54, 55). Changes in microbiota compositions were accompanied by changes in levels of fecal metabolites. However, while levels of metabolites such as valine, leucine, isoleucine, and phenylalanine significantly increased, they returned to the levels observed for the control group in the 9 days following antibiotics. Notably, supplementation of these amino acids has been shown to promote insulin resistance in humans (57) and rats (58) and to increase the risk of type 2 diabetes (59).

We observed many instances of significant correlation between changes in bacterial taxa and metabolites. However, disparities in the dynamics of microbiota and metabolome disruption in response to antibiotics, and subsequent restoration of baseline levels, highlight the importance of assessing both systems. For example, several bacterial taxa, particularly members of the Bacteroidetes phylum, were substantially depleted by ciprofloxacin treatment and did not recover during the 9 days after cessation of the antibiotics. However, ciprofloxacin-driven disruption of the metabolome was less substantial, with its composition following the recovery period broadly in keeping with the baseline. Such discrepancies between the microbiota and metabolome-level antibiotic effects are a further indication of functional redundancy in the intestinal microbiome, the need to examine both composition and function in order to fully characterize the impacts of antibiotic therapy, and the potential influence of antibiotics on host physiology.

Our study had a number of limitations that must be taken into consideration. We examined the impact of antibiotic exposure on the fecal microbiota and correlated this with changes in the compositions of the fecal metabolome. However, establishing direct links between bacterial taxon relative abundances and levels of specific metabolites is extremely challenging. Cross-metabolism between species means that many different taxa can contribute to particular pathways, with functional redundancy between phylogenetically distant populations potentially offsetting reductions in other populations. Furthermore, antibiotics may influence aspects of host physiology, such as appetite, or the production of mucins and immunoglobulins, providing a host-mediated path to metabolome alterations.

Identifying links between antibiotic exposure, microbiota compositional change, an altered fecal metabolome, and disrupted host physiology is extremely challenging. However, the combined use of the sequencing-based bacterial community analysis and metabolomic profiling described here represents an important strategy in elucidating such relationships.

MATERIALS AND METHODS

Murine fecal samples.Feces were collected from female C57BL/6 mice at Charles River, Inc., United Kingdom Ltd. (Margate, United Kingdom) under a commercial license, and all mice were maintained and used in accordance with the Animal Scientific Procedures Act (60) and amendment regulations of 2012 (61). All mice were initially housed in one barrier room within the facility and then transferred to an isolator cage on commencement of antibiotic treatment, where they were divided into three groups of 8 mice each (4 mice per cage). The three groups represented a control group (no treatment) and two antibiotic treatment groups (either vancomycin-imipenem or ciprofloxacin). Six-week-old mice were placed on the same diet (a VRF1 diet, SDS), and antibiotics were dosed at 50 mg/kg of body weight/day delivered in drinking water. Fresh medicated solutions were prepared weekly and stored at 4°C. Water was changed twice weekly, and mice were observed closely for any changes in hydration or for adverse effects; mice were weighed once a week throughout the study. No significant differences in mass were observed between groups during the study, with mice in the control, ciprofloxacin, and vancomycin-imipenem groups attaining masses of 19.26 ± 1.38 g, 19.45 ± 1.11 g, and 19.03 ± 0.81 g, respectively (weights for individual animals are shown in Table S3 in the supplemental material). One fecal pellet was taken from each mouse at three time points: T1, prior to commencement of antibiotic treatment, T2, at termination of treatment after 14 days on antibiotics, and T3, 9 days after antibiotic treatment ceased. After collection, pellets were placed into separate, sterile collection tubes and frozen at −80°C prior to analysis.

TABLE S3

Weights of individual animals in the control and antibiotic treatment groups at each of the three assessment points. Download TABLE S3, XLSX file, 0.01 MB.
Copyright © 2017 Choo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

Sample processing.Known masses of mouse feces were immersed in 1 ml of cold (4°C) 1× sterile phosphate-buffered saline (PBS; pH 7.4) (Thermo Fisher Scientific, United Kingdom) and centrifuged at 13,000 × g for 10 min to form a pellet for DNA extraction. The supernatant was transferred to a fresh microcentrifuge tube for metabolomics study. DNA extraction was performed using a Mo Bio PowerLyzer PowerSoil DNA isolation kit (Mo Bio Laboratories, Carlsbad, CA, USA), as previously described (62).

16S rRNA gene amplicon sequencing.Amplicons of the v4 hypervariable region of the bacterial 16S rRNA gene were generated from DNA extracts, as described previously (62). Briefly, amplicons were generated from 25 PCR cycles, and indexes were attached to the amplicon with 8 PCR cycles. Sequencing was performed on an Illumina MiSeq platform at the David R. Gunn Genomics Facility, South Australian Health and Medical Research Institute. Full details are provided in the supplemental material. Paired-end 16S rRNA gene sequence reads were analyzed with the Quantitative Insights into Microbial Ecology (QIIME) software (v1.8.0) (63) using a bioinformatics pipeline described previously (64). Bar-coded forward and reverse sequencing reads were quality filtered and merged using Paired-End reAd mergeR (PEAR, v0.9.6) (65). Chimeras were detected and filtered from the paired-end reads using USEARCH (v6.1) (66) in a comparison with representative sequences from the Greengenes database (v13.8) that clustered with 97% similarity (67). Operational taxonomic units (OTUs) were assigned to the reads using an open-reference approach with the UCLUST algorithm (v1.2.22q) in a comparison with sequences in SILVA database release 111 (July 2012) (68) that clustered at 97% identity. During the OTU assignment, sequences preclustered at 80% similarity to the reference sequence prior to de novo clustering. All samples were subsampled to 6,250 reads based on the lowest read depth.

1H NMR metabolomics.As described in the paragraphs above, supernatants were obtained from the mouse feces resuspended in cold PBS. Supernatants were frozen by immersion in liquid nitrogen, lyophilized at −58°C overnight, and resuspended in 500 µl D2O. 1H NMR spectra were acquired under automation at 298 K and 700 MHz on a Bruker Avance II 700 NMR spectrometer (Bruker BioSpin, Coventry, United Kingdom) equipped with a 5-mm helium-cooled quadruple resonance cryoprobe and a cooled SampleJet sample changer. The temperature was allowed to stabilize for 3 min after insertion into the magnet. Tuning, matching, and shimming was performed for each sample, and the 1H pulse length was calibrated on each sample and was typically around 12 μs. One-dimensional (1D) Carr-Purcell-Meiboom-Gill presaturation (CPMG-presat) (cpmgpr1d) experiments were acquired with 128 transients, a spectral width of 20.5 ppm, 64,000 data points, a mixing time of 10 ms, a relaxation delay of 4 s, and a total echo delay of 78.7 ms. 1D nuclear Overhauser effect spectroscopy (NOESY)-presat (noesygppr1d) data were also acquired but with 32 transients. Free-induction decays were multiplied with an exponential function (line broadening of 0.3 Hz), Fourier transformed, and calibrated to a 2,2,3,3,-D4-3-(trimethylsilyl) propionic acid (TSP) reference signal at 0.0 ppm. Phase correction was performed manually, and automatic baseline correction was applied. To help in the assignment of metabolite resonances, 2D J-resolved (jresgpprqf) and correlation spectroscopy (COSY) (cosygpprqf) spectra were acquired for a subset of samples. In 2D J-resolved spectra, 64 transients of 8,000 data points were acquired for each of 40 increments, with a relaxation delay of 2 s, during which the solvent signal was suppressed by presaturation. Spectral widths of 16 ppm and 78 Hz were used. In 2D COSY spectra, 8 transients of 4,000 data points were acquired for each of 400 increments, with a relaxation delay of 2 s, during which the solvent signal was suppressed by presaturation. A spectral width of 16 ppm was used in each dimension.

Preprocessing and OPLS-DA (orthogonal projections to latent structures discriminant analysis) were carried out with both MVAPACK (69) and software that was developed in our laboratory for a previous study (70) using the Python programing language with NumPy and SciPy for calculations and matplotlib for visualization. The nonlinear iterative partial least-squares (NIPALS) algorithm (71) was used for OPLS-DA. For integrated microbiome and metabolomic analysis, consensus OPLS were created using the K-OPLS R package (72).

Regions above 8.5 ppm and below 0.5 ppm were excluded because of noise content. The water peak and TSP reference signal were also excluded. Spectra further aligned with the icoshift (73) algorithm and bucketed with the optimized bucketing algorithm (74), using a 0.005-ppm minimum bin size, leaving 1,652 data points per spectrum. These spectra were subjected to probabilistic quotient normalization (PQN) and Pareto scaled (75).

Data analysis and statistics.Microbial data were analyzed for alpha diversity measures (taxon richness, S; Shannon-Wiener index, H; Simpson diversity index, 1 to D) of the microbial community using PAST (v3.04) (76). Operational taxonomic unit (OTU) relative abundance was imported into the Primer-E software (v.6; Primer-E Ltd., Plymouth, United Kingdom) for beta diversity analysis. Bray-Curtis similarities were calculated based on the square root-transformed OTU relative abundances and were used in the nonmetric multidimensional scaling (NMDS) ordination plot. The permutational analysis of variance (PERMANOVA) model was used for testing the null hypothesis of no difference between the compared groups (54), based on the parameters’ permutation of residuals under a reduced model and a type III sum of squares. PERMDISP was used to assess the dispersion of the microbial community within the groups (54). Microbial differences between groups at the taxonomic level were tested for statistical significance based on the Kruskal-Wallis and Mann-Whitney tests using GraphPad Prism 6 (GraphPad Software, Inc., La Jolla, USA). The LEfSe algorithm was performed using the Galaxy application tool (http://huttenhower.sph.harvard.edu/galaxy/ ) with a linear discriminant analysis cutoff score of 3.0 and a P of <0.05 for statistical significance (55). Functional predictions of microbial community were performed using PICRUSt on closed-reference OTUs with 97% identity based on the Greengenes database (v13.5) (56). The OTUs were normalized on PICRUSt and used for the prediction of KEGG orthologs (KOs). The predicted metagenome was statistically analyzed on STAMP using Welch’s t test with the Benjamini-Hochberg correction for the false-discovery rate (FDR) and filtered to retain features with an effect size (ratio of proportions) of greater than two (77). A heatmap was generated using the ggplots2 package (v2.0.0), and a dendrogram of bacterial taxa was generated based on the Bray-Curtis distances and hierarchical clustering performed using Ward’s method (78, 79). Two and three mice, respectively, were excluded from the microbiota and metabolome analysis, as samples from them failed quality control thresholds.

For analysis of the metabolome data, cross-validation was performed in the same manner for OPLS-DA and consensus OPLS-DA. Seventy-five percent of the samples were used as a training set, and the remaining 25% were used as a test set, ensuring that the number of samples in the test set was proportional to the total number of samples from each class and that at least one sample from each class was present in the test set. To choose the number of components for the model, a leave-one-out cross-validation step was carried out on the samples in the training set, and the F1 was used to choose the number of components, with the additional constraint that a maximum of 8 components was used. A double cross-validation was repeated 2,000 times for each group with randomly chosen samples in the training and test sets to prevent bias due to the choice of training or test set.

This procedure was repeated with randomly generated class assignments to provide a reference value for Q2. The chosen number of components minus one was then used as an OPLS filter, and a PLS-DA with two components was carried out on the filtered data to yield one predictive and one orthogonal component. In the back-scaled loadings analysis, peaks that allow the models to distinguish between classes were assigned by comparing chemical-shift values and multiplicities from J-resolved NMR spectra to values from the BMRB (80), HMDB (81), and Chenomx software. Volcano plots were generated in R. Five samples were excluded from the metabolome analysis, as these samples failed quality control thresholds. Microbiota and metabolome associations based on Spearman’s correlation were performed in R using the corrplot package. Bacterial taxa and metabolites were sorted according to hierarchical clustering based on Euclidean distances, and Ward’s method was used to perform cluster analysis.

TEXT S1

16S rRNA gene amplicon sequencing. Download TEXT S1, PDF file, 0.02 MB.
Copyright © 2017 Choo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

Accession number(s).Microbial sequencing data have been deposited in the Sequence Read Archive (SRA) database under GenBank accession number SRP096906 . Metabolite data have been deposited in the MetaboLights database under accession number MTBLS422 .

ACKNOWLEDGMENTS

We thank R. A. Atkinson of the KCL Centre for Biomolecular Spectroscopy for his assistance with the NMR experiments.

The Wellcome Trust supported the KCL Centre with a Capital Award. This research was funded/supported by the National Institute for Health Research (NIHR) Biomedical Research Centre, based at Guy’s and St. Thomas' NHS Foundation Trust, and King’s College London (Ph.D. studentship to T.K.). The study was also supported by a BBSRC Impact Acceleration Award to King’s College London.

The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR, or the Department of Health.

FOOTNOTES

    • Received January 5, 2017.
    • Accepted January 13, 2017.
  • Copyright © 2017 Choo et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license .

REFERENCES

  1. 1.↵
    1. Yano JM,
    2. Yu K,
    3. Donaldson GP,
    4. Shastri GG,
    5. Ann P,
    6. Ma L,
    7. Nagler CR,
    8. Ismagilov RF,
    9. Mazmanian SK,
    10. Hsiao EY
    . 2015. Indigenous bacteria from the gut microbiota regulate host serotonin biosynthesis. Cell161:264–276. doi:10.1016/j.cell.2015.02.047.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Crane JD,
    2. Palanivel R,
    3. Mottillo EP,
    4. Bujak AL,
    5. Wang H,
    6. Ford RJ,
    7. Collins A,
    8. Blümer RM,
    9. Fullerton MD,
    10. Yabut JM,
    11. Kim JJ,
    12. Ghia JE,
    13. Hamza SM,
    14. Morrison KM,
    15. Schertzer JD,
    16. Dyck JRB,
    17. Khan WI,
    18. Steinberg GR
    . 2015. Inhibiting peripheral serotonin synthesis reduces obesity and metabolic dysfunction by promoting brown adipose tissue thermogenesis. Nat Med21:166–172. doi:10.1038/nm.3766.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Fu J,
    2. Bonder MJ,
    3. Cenit MC,
    4. Tigchelaar EF,
    5. Maatman A,
    6. Dekens JAM,
    7. Brandsma E,
    8. Marczynska J,
    9. Imhann F,
    10. Weersma RK,
    11. Franke L,
    12. Poon TW,
    13. Xavier RJ,
    14. Gevers D,
    15. Hofker MH,
    16. Wijmenga C,
    17. Zhernakova A
    . 2015. The gut microbiome contributes to a substantial proportion of the variation in blood lipids. Circ Res117:817–824. doi:10.1161/CIRCRESAHA.115.306807.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Valentini M,
    2. Piermattei A,
    3. Di Sante G,
    4. Migliara G,
    5. Delogu G,
    6. Ria F
    . 2014. Immunomodulation by gut microbiota: role of Toll-like receptor expressed by T cells. J Immunol Res2014:586939. doi:10.1155/2014/586939.
    OpenUrlCrossRef
  5. 5.↵
    1. Sorini C,
    2. Falcone M
    . 2013. Shaping the (auto)immune response in the gut: the role of intestinal immune regulation in the prevention of type 1 diabetes. Am J Clin Exp Immunol2:156–171.
    OpenUrlPubMed
  6. 6.↵
    1. Cryan JF,
    2. O’Mahony SM
    . 2011. The microbiome-gut-brain axis: from bowel to behavior. Neurogastroenterol Motil23:187–192. doi:10.1111/j.1365-2982.2010.01664.x.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Rogers GB,
    2. Keating DJ,
    3. Young RL,
    4. Wong ML,
    5. Licinio J,
    6. Wesselingh S
    . 2016. From gut dysbiosis to altered brain function and mental illness: mechanisms and pathways. Mol Psychiatry21:738–748. doi:10.1038/mp.2016.50.
    OpenUrlCrossRef
  8. 8.↵
    1. Britton RA,
    2. Young VB
    . 2014. Role of the intestinal microbiota in resistance to colonization by Clostridium difficile. Gastroenterology146:1547–1553. doi:10.1053/j.gastro.2014.01.059.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Theriot CM,
    2. Koenigsknecht MJ,
    3. Carlson PE Jr,
    4. Hatton GE,
    5. Nelson AM,
    6. Li B,
    7. Huffnagle GB,
    8. Li JZ,
    9. Young VB
    . 2014. Antibiotic-induced shifts in the mouse gut microbiome and metabolome increase susceptibility to Clostridium difficile infection. Nat Commun5:3114. doi:10.1038/ncomms4114.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Zhao Y,
    2. Wu J,
    3. Li JV,
    4. Zhou NY,
    5. Tang H,
    6. Wang Y
    . 2013. Gut microbiota composition modifies fecal metabolic profiles in mice. J Proteome Res12:2987–2999. doi:10.1021/pr400263n.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Turta O,
    2. Rautava S
    . 2016. Antibiotics, obesity and the link to microbes—what are we doing to our children?BMC Med14:57. doi:10.1186/s12916-016-0605-7.
    OpenUrlCrossRef
  12. 12.↵
    1. Cox LM,
    2. Blaser MJ
    . 2015. Antibiotics in early life and obesity. Nat Rev Endocrinol11:182–190. doi:10.1038/nrendo.2014.210.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Mueller NT,
    2. Whyatt R,
    3. Hoepner L,
    4. Oberfield S,
    5. Dominguez-Bello MG,
    6. Widen EM,
    7. Hassoun A,
    8. Perera F,
    9. Rundle A
    . 2015. Prenatal exposure to antibiotics, cesarean section and risk of childhood obesity. Int J Obes (Lond)39:665–670. doi:10.1038/ijo.2014.180.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Azad MB,
    2. Bridgman SL,
    3. Becker AB,
    4. Kozyrskyj AL
    . 2014. Infant antibiotic exposure and the development of childhood overweight and central adiposity. Int J Obes (Lond)38:1290–1298. doi:10.1038/ijo.2014.119.
    OpenUrlCrossRef
  15. 15.↵
    1. Saari A,
    2. Virta LJ,
    3. Sankilampi U,
    4. Dunkel L,
    5. Saxen H
    . 2015. Antibiotic exposure in infancy and risk of being overweight in the first 24 months of life. Pediatrics135:617–626. doi:10.1542/peds.2014-3407.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Schwartz BS,
    2. Pollak J,
    3. Bailey-Davis L,
    4. Hirsch AG,
    5. Cosgrove SE,
    6. Nau C,
    7. Kress AM,
    8. Glass TA,
    9. Bandeen-Roche K
    . 2016. Antibiotic use and childhood body mass index trajectory. Int J Obes (Lond)40:615–621. doi:10.1038/ijo.2015.218.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Million M,
    2. Thuny F,
    3. Angelakis E,
    4. Casalta JP,
    5. Giorgi R,
    6. Habib G,
    7. Raoult D
    . 2013. Lactobacillus reuteri and Escherichia coli in the human gut microbiota may predict weight gain associated with vancomycin treatment. Nutr Diabetes3:e87. doi:10.1038/nutd.2013.28.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Korpela K,
    2. Salonen A,
    3. Virta LJ,
    4. Kekkonen RA,
    5. Forslund K,
    6. Bork P,
    7. de Vos WM
    . 2016. Intestinal microbiome is related to lifetime antibiotic use in Finnish pre-school children. Nat Commun7:10410. doi:10.1038/ncomms10410.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Moya A,
    2. Ferrer M
    . 2016. Functional redundancy-induced stability of gut microbiota subjected to disturbance. Trends Microbiol24:402–413. doi:10.1016/j.tim.2016.02.002.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Ferrer M,
    2. Ruiz A,
    3. Lanza F,
    4. Haange SB,
    5. Oberbach A,
    6. Till H,
    7. Bargiela R,
    8. Campoy C,
    9. Segura MT,
    10. Richter M,
    11. von Bergen M,
    12. Seifert J,
    13. Suarez A
    . 2013. Microbiota from the distal guts of lean and obese adolescents exhibit partial functional redundancy besides clear differences in community structure. Environ Microbiol15:211–226. doi:10.1111/j.1462-2920.2012.02845.x.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Furusawa Y,
    2. Obata Y,
    3. Fukuda S,
    4. Endo TA,
    5. Nakato G,
    6. Takahashi D,
    7. Nakanishi Y,
    8. Uetake C,
    9. Kato K,
    10. Kato T,
    11. Takahashi M,
    12. Fukuda NN,
    13. Murakami S,
    14. Miyauchi E,
    15. Hino S,
    16. Atarashi K,
    17. Onawa S,
    18. Fujimura Y,
    19. Lockett T,
    20. Clarke JM,
    21. Topping DL,
    22. Tomita M,
    23. Hori S,
    24. Ohara O,
    25. Morita T,
    26. Koseki H,
    27. Kikuchi J,
    28. Honda K,
    29. Hase K,
    30. Ohno H
    . 2013. Commensal microbe-derived butyrate induces the differentiation of colonic regulatory T cells. Nature504:446–450. doi:10.1038/nature12721.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    1. Park BO,
    2. Kim SH,
    3. Kong GY,
    4. Kim DH,
    5. Kwon MS,
    6. Lee SU,
    7. Kim MO,
    8. Cho S,
    9. Lee S,
    10. Lee HJ,
    11. Han SB,
    12. Kwak YS,
    13. Lee SB,
    14. Kim S
    . 2016. Selective novel inverse agonists for human GPR43 augment GLP-1 secretion. Eur J Pharmacol771:1–9. doi:10.1016/j.ejphar.2015.12.010.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Ha CW,
    2. Lam YY,
    3. Holmes AJ
    . 2014. Mechanistic links between gut microbial community dynamics, microbial functions and metabolic health. World J Gastroenterol20:16498–16517. doi:10.3748/wjg.v20.i44.16498.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Kimura I,
    2. Inoue D,
    3. Maeda T,
    4. Hara T,
    5. Ichimura A,
    6. Miyauchi S,
    7. Kobayashi M,
    8. Hirasawa A,
    9. Tsujimoto G
    . 2011. Short-chain fatty acids and ketones directly regulate sympathetic nervous system via G protein-coupled receptor 41 (GPR41). Proc Natl Acad Sci U S A108:8030–8035. doi:10.1073/pnas.1016088108.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Nøhr MK,
    2. Pedersen MH,
    3. Gille A,
    4. Egerod KL,
    5. Engelstoft MS,
    6. Husted AS,
    7. Sichlau RM,
    8. Grunddal KV,
    9. Poulsen SS,
    10. Han S,
    11. Jones RM,
    12. Offermanns S,
    13. Schwartz TW
    . 2013. GPR41/FFAR3 and GPR43/FFAR2 as cosensors for short-chain fatty acids in enteroendocrine cells vs FFAR3 in enteric neurons and FFAR2 in enteric leukocytes. Endocrinology154:3552–3564. doi:10.1210/en.2013-1142.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    1. Donohoe DR,
    2. Garge N,
    3. Zhang X,
    4. Sun W,
    5. O’Connell TM,
    6. Bunger MK,
    7. Bultman SJ
    . 2011. The microbiome and butyrate regulate energy metabolism and autophagy in the mammalian colon. Cell Metab13:517–526. doi:10.1016/j.cmet.2011.02.018.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    1. den Besten G,
    2. van Eunen K,
    3. Groen AK,
    4. Venema K,
    5. Reijngoud D-J,
    6. Bakker BM
    . 2013. The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J Lipid Res54:2325–2340. doi:10.1194/jlr.R036012.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Macfarlane S,
    2. Macfarlane GT
    . 2003. Regulation of short-chain fatty acid production. Proc Nutr Soc62:67–72. doi:10.1079/PNS2002207.
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    1. Rogers GB,
    2. Kozlowska J,
    3. Keeble J,
    4. Metcalfe K,
    5. Fao M,
    6. Dowd SE,
    7. Mason AJ,
    8. McGuckin MA,
    9. Bruce KD
    . 2014. Functional divergence in gastrointestinal microbiota in physically separated genetically identical mice. Sci Rep4:5437. doi:10.1038/srep05437.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Berger M,
    2. Gray JA,
    3. Roth BL
    . 2009. The expanded biology of serotonin. Annu Rev Med60:355–366. doi:10.1146/annurev.med.60.042307.110802.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Hoentjen F,
    2. Harmsen HJ,
    3. Braat H,
    4. Torrice CD,
    5. Mann BA,
    6. Sartor RB,
    7. Dieleman LA
    . 2003. Antibiotics with a selective aerobic or anaerobic spectrum have different therapeutic activities in various regions of the colon in interleukin 10 gene deficient mice. Gut52:1721–1727. doi:10.1136/gut.52.12.1721.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Cani PD,
    2. Geurts L,
    3. Matamoros S,
    4. Plovier H,
    5. Duparc T
    . 2014. Glucose metabolism: focus on gut microbiota, the endocannabinoid system and beyond. Diabetes Metab40:246–257. doi:10.1016/j.diabet.2014.02.004.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Backhed F,
    2. Ding H,
    3. Wang T,
    4. Hooper LV,
    5. Koh GY,
    6. Nagy A,
    7. Semenkovich CF,
    8. Gordon JI
    . 2004. The gut microbiota as an environmental factor that regulates fat storage. Proc Natl Acad Sci U S A101:15718–15723. doi:10.1073/pnas.0407076101.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Turnbaugh PJ,
    2. Gordon JI
    . 2009. The core gut microbiome, energy balance and obesity. J Physiol587:4153–4158. doi:10.1113/jphysiol.2009.174136.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Ridlon JM,
    2. Kang DJ,
    3. Hylemon PB,
    4. Bajaj JS
    . 2014. Bile acids and the gut microbiome. Curr Opin Gastroenterol30:332–338. doi:10.1097/MOG.0000000000000057.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Marcobal A,
    2. Kashyap PC,
    3. Nelson TA,
    4. Aronov PA,
    5. Donia MS,
    6. Spormann A,
    7. Fischbach MA,
    8. Sonnenburg JL
    . 2013. A metabolomic view of how the human gut microbiota impacts the host metabolome using humanized and gnotobiotic mice. ISME J7:1933–1943. doi:10.1038/ismej.2013.89.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Sayin SI,
    2. Wahlström A,
    3. Felin J,
    4. Jäntti S,
    5. Marschall HU,
    6. Bamberg K,
    7. Angelin B,
    8. Hyötyläinen T,
    9. Orešič M,
    10. Bäckhed F
    . 2013. Gut microbiota regulates bile acid metabolism by reducing the levels of tauro-beta-muricholic acid, a naturally occurring FXR antagonist. Cell Metab17:225–235. doi:10.1016/j.cmet.2013.01.003.
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.↵
    1. Krajmalnik-Brown R,
    2. Ilhan ZE,
    3. Kang DW,
    4. DiBaise JK
    . 2012. Effects of gut microbes on nutrient absorption and energy regulation. Nutr Clin Pract27:201–214. doi:10.1177/0884533611436116.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Carabotti M,
    2. Scirocco A,
    3. Maselli MA,
    4. Severi C
    . 2015. The gut-brain axis: interactions between enteric microbiota, central and enteric nervous systems. Ann Gastroenterol28:203–209.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Maurice CF
    . 2013. Xenobiotics and the active gut microbiome: unknown effects unveiled. Med Sci (Paris)29:846–848. (In French.) doi:10.1051/medsci/20132910011.
    OpenUrlCrossRef
  41. 41.↵
    1. Scott FI,
    2. Horton DB,
    3. Mamtani R,
    4. Haynes K,
    5. Goldberg DS,
    6. Lee DY,
    7. Lewis JD
    . 2016. Administration of antibiotics to children before age 2 years increases risk for childhood obesity. Gastroenterology151:120.e5–129.e5. doi:10.1053/j.gastro.2016.03.006.
    OpenUrlCrossRef
  42. 42.↵
    1. Mahana D,
    2. Trent CM,
    3. Kurtz ZD,
    4. Bokulich NA,
    5. Battaglia T,
    6. Chung J,
    7. Müller CL,
    8. Li H,
    9. Bonneau RA,
    10. Blaser MJ
    . 2016. Antibiotic perturbation of the murine gut microbiome enhances the adiposity, insulin resistance, and liver disease associated with high-fat diet. Genome Med8:48. doi:10.1186/s13073-016-0297-9.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. van Opstal E,
    2. Kolling GL,
    3. Moore JH II,
    4. Coquery CM,
    5. Wade NS,
    6. Loo WM,
    7. Bolick DT,
    8. Shin JH,
    9. Erickson LD,
    10. Warren CA
    . 2016. Vancomycin treatment alters humoral immunity and intestinal microbiota in an aged mouse model of Clostridium difficile infection. J Infect Dis214:130–139. doi:10.1093/infdis/jiw071.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Bech-Nielsen GV,
    2. Hansen CHF,
    3. Hufeldt MR,
    4. Nielsen DS,
    5. Aasted B,
    6. Vogensen FK,
    7. Midtvedt T,
    8. Hansen AK
    . 2012. Manipulation of the gut microbiota in C57BL/6 mice changes glucose tolerance without affecting weight development and gut mucosal immunity. Res Vet Sci92:501–508. doi:10.1016/j.rvsc.2011.04.005.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Sekirov I,
    2. Tam NM,
    3. Jogova M,
    4. Robertson ML,
    5. Li Y,
    6. Lupp C,
    7. Finlay BB
    . 2008. Antibiotic-induced perturbations of the intestinal microbiota alter host susceptibility to enteric infection. Infect Immun76:4726–4736. doi:10.1128/IAI.00319-08.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Yap IKS,
    2. Li JV,
    3. Saric J,
    4. Martin FP,
    5. Davies H,
    6. Wang Y,
    7. Wilson ID,
    8. Nicholson JK,
    9. Utzinger J,
    10. Marchesi JR,
    11. Holmes E
    . 2008. Metabonomic and microbiological analysis of the dynamic effect of vancomycin-induced gut microbiota modification in the mouse. J Proteome Res7:3718–3728. doi:10.1021/pr700864x.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    1. Jump RLP,
    2. Polinkovsky A,
    3. Hurless K,
    4. Sitzlar B,
    5. Eckart K,
    6. Tomas M,
    7. Deshpande A,
    8. Nerandzic MM,
    9. Donskey CJ
    . 2014. Metabolomics analysis identifies intestinal microbiota-derived biomarkers of colonization resistance in clindamycin-treated mice. PLoS One9:e101267. doi:10.1371/journal.pone.0101267.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Arthur JC,
    2. Gharaibeh RZ,
    3. Mühlbauer M,
    4. Perez-Chanona E,
    5. Uronis JM,
    6. McCafferty J,
    7. Fodor AA,
    8. Jobin C
    . 2014. Microbial genomic analysis reveals the essential role of inflammation in bacteria-induced colorectal cancer. Nat Commun5:4724. doi:10.1038/ncomms5724.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Lecleire S,
    2. Hassan A,
    3. Marion-Letellier R,
    4. Antonietti M,
    5. Savoye G,
    6. Bole-Feysot C,
    7. Lerebours E,
    8. Ducrotte P,
    9. Dechelotte P,
    10. Coeffier M
    . 2008. Combined glutamine and arginine decrease proinflammatory cytokine production by biopsies from Crohn’s patients in association with changes in nuclear factor-kappaB and p38 mitogen-activated protein kinase pathways. J Nutr138:2481–2486. doi:10.3945/jn.108.099127.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    1. Vicario M,
    2. Amat C,
    3. Rivero M,
    4. Moreto M,
    5. Pelegri C
    . 2007. Dietary glutamine affects mucosal functions in rats with mild DSS-induced colitis. J Nutr137:1931–1937.
    OpenUrlAbstract/FREE Full Text
  51. 51.↵
    1. Wanasen N,
    2. Soong L
    . 2008. l-Arginine metabolism and its impact on host immunity against Leishmania infection. Immunol Res41:15–25. doi:10.1007/s12026-007-8012-y.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Ren W,
    2. Chen S,
    3. Yin J,
    4. Duan J,
    5. Li T,
    6. Liu G,
    7. Feng Z,
    8. Tan B,
    9. Yin Y,
    10. Wu G
    . 2014. Dietary arginine supplementation of mice alters the microbial population and activates intestinal innate immunity. J Nutr144:988–995. doi:10.3945/jn.114.192120.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    1. Zabaleta J,
    2. McGee DJ,
    3. Zea AH,
    4. Hernandez CP,
    5. Rodriguez PC,
    6. Sierra RA,
    7. Correa P,
    8. Ochoa AC
    . 2004. Helicobacter pylori arginase inhibits T cell proliferation and reduces the expression of the TCR zeta-chain (CD3zeta). J Immunol173:586–593. doi:10.4049/jimmunol.173.1.586.
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    1. Anderson MJ,
    2. Walsh DCI
    . 2013. PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: what null hypothesis are you testing?Ecol Monogr83:557–574. doi:10.1890/12-2010.1.
    OpenUrlCrossRef
  55. 55.↵
    1. Segata N,
    2. Izard J,
    3. Waldron L,
    4. Gevers D,
    5. Miropolsky L,
    6. Garrett WS,
    7. Huttenhower C
    . 2011. Metagenomic biomarker discovery and explanation. Genome Biol12:R60. doi:10.1186/gb-2011-12-6-r60.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Langille MGI,
    2. Zaneveld J,
    3. Caporaso JG,
    4. McDonald D,
    5. Knights D,
    6. Reyes JA,
    7. Clemente JC,
    8. Burkepile DE,
    9. Vega Thurber RL,
    10. Knight R,
    11. Beiko RG,
    12. Huttenhower C
    . 2013. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol31:814–821. doi:10.1038/nbt.2676.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Tremblay F,
    2. Krebs M,
    3. Dombrowski L,
    4. Brehm A,
    5. Bernroider E,
    6. Roth E,
    7. Nowotny P,
    8. Waldhausl W,
    9. Marette A,
    10. Roden M
    . 2005. Overactivation of S6 kinase 1 as a cause of human insulin resistance during increased amino acid availability. Diabetes54:2674–2684. doi:10.2337/diabetes.54.9.2674.
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    1. Xiao F,
    2. Huang Z,
    3. Li H,
    4. Yu J,
    5. Wang C,
    6. Chen S,
    7. Meng Q,
    8. Cheng Y,
    9. Gao X,
    10. Li J,
    11. Liu Y,
    12. Guo F
    . 2011. Leucine deprivation increases hepatic insulin sensitivity via GCN2/mTOR/S6K1 and AMPK pathways. Diabetes60:746–756. doi:10.2337/db10-1246.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    1. Wang TJ,
    2. Larson MG,
    3. Vasan RS,
    4. Cheng S,
    5. Rhee EP,
    6. McCabe E,
    7. Lewis GD,
    8. Fox CS,
    9. Jacques PF,
    10. Fernandez C,
    11. O’Donnell CJ,
    12. Carr SA,
    13. Mootha VK,
    14. Florez JC,
    15. Souza A,
    16. Melander O,
    17. Clish CB,
    18. Gerszten RE
    . 2011. Metabolite profiles and the risk of developing diabetes. Nat Med17:448–453. doi:10.1038/nm.2307.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Parliament, United Kingdom
    . 1986. Animals (Scientific Procedures) Act 1986. Parliament, London, United Kingdom.
  61. 61.↵
    1. Parliament, United Kingdom
    . 2012. Animals (Scientific Procedures) Act 1986 amendment regulations 2012. Parliament, London, United Kingdom.
  62. 62.↵
    1. Choo JM,
    2. Leong LEX,
    3. Rogers GB
    . 2015. Sample storage significantly influence faecal microbiome profiles. Sci Rep5:16350. doi:10.1038/srep16350.
    OpenUrlCrossRefPubMed
  63. 63.↵
    1. Caporaso JG,
    2. Kuczynski J,
    3. Stombaugh J,
    4. Bittinger K,
    5. Bushman FD,
    6. Costello EK,
    7. Fierer N,
    8. Gonzalez Pena AG,
    9. Goodrich JK,
    10. Gordon JI,
    11. Huttley GA,
    12. Kelley ST,
    13. Kinghts D,
    14. Koenig JE,
    15. Ley RE,
    16. Lozupone CA,
    17. McDonald D,
    18. Muegge BD,
    19. Pirrung M,
    20. Reeder J,
    21. JRS,
    22. Turnbaugh PJ,
    23. Walters WA,
    24. Widmann J,
    25. Yatsunenko T,
    26. Zaneveld J,
    27. Knight R
    . 2010. QIIME allows analysis of high-throughput community sequencing data. Nat Methods7:335–336. doi:10.1038/nmeth.f.303.
    OpenUrlCrossRefPubMedWeb of Science
  64. 64.↵
    1. Jervis-Bardy J,
    2. Leong LEX,
    3. Marri S,
    4. Smith RJ,
    5. Choo JM,
    6. Smith-Vaughan HC,
    7. Nosworthy E,
    8. Morris PS,
    9. O’Leary S,
    10. Rogers GB,
    11. Marsh RL
    . 2015. Deriving accurate microbiota profiles from human samples with low bacterial content through post-sequencing processing of Illumina MiSeq data. Microbiome3:19. doi:10.1186/s40168-015-0083-8.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Zhang J,
    2. Kobert K,
    3. Flouri T,
    4. Stamatakis A
    . 2014. PEAR: a fast and accurate Illumina paired-end reAd mergeR. Bioinformatics30:614–620. doi:10.1093/bioinformatics/btt593.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Edgar RC
    . 2010. Search and clustering orders of magnitude faster than BLAST. Bioinformatics26:2460–2461. doi:10.1093/bioinformatics/btq461.
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.↵
    1. McDonald D,
    2. Price MN,
    3. Goodrich J,
    4. Nawrocki EP,
    5. DeSantis TZ,
    6. Probst A,
    7. Andersen GL,
    8. Knight R,
    9. Hugenholtz P
    . 2012. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J6:610–618. doi:10.1038/ismej.2011.139.
    OpenUrlCrossRefPubMedWeb of Science
  68. 68.↵
    1. Quast C,
    2. Pruesse E,
    3. Yilmaz P,
    4. Gerken J,
    5. Schweer T,
    6. Yarza P,
    7. Peplies J,
    8. Glockner FO
    . 2013. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res41:D590–D596. doi:10.1093/nar/gks1219.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Worley B,
    2. Powers R
    . 2014. MVAPACK: a complete data handling package for NMR metabolomics. ACS Chem Biol9:1138–1144. doi:10.1021/cb4008937.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Vermeer LS,
    2. Fruhwirth GO,
    3. Pandya P,
    4. Ng T,
    5. Mason AJ
    . 2012. NMR metabolomics of MTLn3E breast cancer cells identifies a role for CXCR4 in lipid and choline regulation. J Proteome Res11:2996–3003. doi:10.1021/pr300111x.
    OpenUrlCrossRefPubMed
  71. 71.↵
    1. Andersson M
    . 2009. A comparison of nine PLS1 algorithms. J Chemometrics23:518–529. doi:10.1002/cem.1248.
    OpenUrlCrossRef
  72. 72.↵
    1. Boccard J,
    2. Rutledge DN
    . 2013. A consensus orthogonal partial least squares discriminant analysis (OPLS-DA) strategy for multiblock omics data fusion. Anal Chim Acta769:30–39. doi:10.1016/j.aca.2013.01.022.
    OpenUrlCrossRefPubMed
  73. 73.↵
    1. Savorani F,
    2. Tomasi G,
    3. Engelsen SB
    . 2010. icoshift: a versatile tool for the rapid alignment of 1D NMR spectra. J Magn Reson202:190–202. doi:10.1016/j.jmr.2009.11.012.
    OpenUrlCrossRefPubMedWeb of Science
  74. 74.↵
    1. Sousa FFO,
    2. Luzardo-Álvarez A,
    3. Blanco-Méndez J,
    4. Otero-Espinar FJ,
    5. Martín-Pastor M,
    6. Sández Macho I
    . 2013. Use of 1H NMR STD, waterLOGSY, and Langmuir monolayer techniques for characterization of drug-zein protein complexes. Eur J Pharm Biopharm85:790–798. doi:10.1016/j.ejpb.2013.07.008.
    OpenUrlCrossRef
  75. 75.↵
    1. Dieterle F,
    2. Ross A,
    3. Schlotterbeck G,
    4. Senn H
    . 2006. Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. Application in 1H NMR metabonomics. Anal Chem78:4281–4290. doi:10.1021/ac051632c.
    OpenUrlCrossRefPubMed
  76. 76.↵
    1. Hammer O,
    2. Harper DAT,
    3. Ryan PD
    . 2001. PAST: paleontological statistics software package for education and data analysis. Paleontol Electron4:9.
    OpenUrl
  77. 77.↵
    1. Parks DH,
    2. Tyson GW,
    3. Hugenholtz P,
    4. Beiko RG
    . 2014. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics30:3123–3124. doi:10.1093/bioinformatics/btu494.
    OpenUrlCrossRefPubMed
  78. 78.↵
    1. Wickham H
    . 2009. ggplot2: elegant graphics for data analysis. Springer-Verlag, New York, NY. http://ggplot2.org .
  79. 79.↵
    1. Team RC
    . 2015. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
  80. 80.↵
    1. Ulrich EL,
    2. Akutsu H,
    3. Doreleijers JF,
    4. Harano Y,
    5. Ioannidis YE,
    6. Lin J,
    7. Livny M,
    8. Mading S,
    9. Maziuk D,
    10. Miller Z,
    11. Nakatani E,
    12. Schulte CF,
    13. Tolmie DE,
    14. Kent Wenger R,
    15. Yao H,
    16. Markley JL
    . 2008. BioMagResBank. Nucleic Acids Res36:D402–D408. doi:10.1093/nar/gkm957.
    OpenUrlCrossRefPubMedWeb of Science
  81. 81.↵
    1. Wishart DS
    . 2010. Computational approaches to metabolomics. Methods Mol Biol593:283–313. doi:10.1007/978-1-60327-194-3_14.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Divergent Relationships between Fecal Microbiota and Metabolome following Distinct Antibiotic-Induced Disruptions
Jocelyn M. Choo, Tokuwa Kanno, Nur Masirah Mohd Zain, Lex E. X. Leong, Guy C. J. Abell, Julie E. Keeble, Kenneth D. Bruce, A. James Mason, Geraint B. Rogers
mSphere Feb 2017, 2 (1) e00005-17; DOI: 10.1128/mSphere.00005-17

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Divergent Relationships between Fecal Microbiota and Metabolome following Distinct Antibiotic-Induced Disruptions
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Divergent Relationships between Fecal Microbiota and Metabolome following Distinct Antibiotic-Induced Disruptions
Jocelyn M. Choo, Tokuwa Kanno, Nur Masirah Mohd Zain, Lex E. X. Leong, Guy C. J. Abell, Julie E. Keeble, Kenneth D. Bruce, A. James Mason, Geraint B. Rogers
mSphere Feb 2017, 2 (1) e00005-17; DOI: 10.1128/mSphere.00005-17
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

beta-lactam
fluoroquinolone
glycopeptide
metabolic activity
microbiota composition

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5042