Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • 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
    • Special Issues
    • COVID-19 Special Collection
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • 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
  • About the Journal
    • About mSystems
    • 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
    • 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
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • 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
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Observation | Novel Systems Biology Techniques

Evaluating Organism-Wide Changes in the Metabolome and Microbiome following a Single Dose of Antibiotic

Alison Vrbanac, Kathryn A. Patras, Alan K. Jarmusch, Robert H. Mills, Samuel R. Shing, Robert A. Quinn, Fernando Vargas, David J. Gonzalez, Pieter C. Dorrestein, Rob Knight, Victor Nizet
Mariana X. Byndloss, Editor
Alison Vrbanac
aDepartment of Pediatrics, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kathryn A. Patras
bDepartment of Molecular Virology and Microbiology, Baylor College of Medicine, Houston, Texas, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kathryn A. Patras
Alan K. Jarmusch
cSkaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Alan K. Jarmusch
Robert H. Mills
aDepartment of Pediatrics, UC San Diego, La Jolla, California, USA
cSkaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, La Jolla, California, USA
dDepartment of Pharmacology, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Samuel R. Shing
aDepartment of Pediatrics, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert A. Quinn
eDepartment of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fernando Vargas
cSkaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David J. Gonzalez
cSkaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, La Jolla, California, USA
dDepartment of Pharmacology, UC San Diego, La Jolla, California, USA
fCenter for Microbiome Innovation, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for David J. Gonzalez
Pieter C. Dorrestein
aDepartment of Pediatrics, UC San Diego, La Jolla, California, USA
cSkaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, La Jolla, California, USA
fCenter for Microbiome Innovation, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rob Knight
aDepartment of Pediatrics, UC San Diego, La Jolla, California, USA
fCenter for Microbiome Innovation, UC San Diego, La Jolla, California, USA
gDepartment of Computer Science and Engineering, UC San Diego, La Jolla, California, USA
hDepartment of Bioengineering, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Rob Knight
Victor Nizet
aDepartment of Pediatrics, UC San Diego, La Jolla, California, USA
cSkaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, La Jolla, California, USA
fCenter for Microbiome Innovation, UC San Diego, La Jolla, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Victor Nizet
Mariana X. Byndloss
Vanderbilt University Medical Center
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00340-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Antibiotics are a mainstay of modern medicine, but as they kill their target pathogen(s), they often affect the commensal microbiota. Antibiotic-induced microbiome dysbiosis is a growing research focus and health concern, often assessed via analysis of fecal samples. However, such analysis does not inform how antibiotics influence the microbiome across the whole host or how such changes subsequently alter host chemistry. In this study, we investigated the acute (1 day postadministration) and delayed (6 days postadministration) effects of a single parenteral dose of two common antibiotics, ampicillin or vancomycin, on the global metabolome and microbiome of mice across 77 different body sites from 25 different organs. The broader-spectrum agent ampicillin had the greatest impact on the microbiota in the lower gastrointestinal tract (cecum and colon), where microbial diversity is highest. In the metabolome, the greatest effects were seen 1 day posttreatment, and changes in metabolite abundances were not confined to the gut. The local abundance of ampicillin and its metabolites correlated with increased metabolome effect size and a loss of alpha diversity versus control mice. Additionally, small peptides were elevated in the lower gastrointestinal tract of mice 1 day after antibiotic treatment. While a single parenteral dose of antibiotic did not drastically alter the microbiome, nevertheless, changes in the metabolome were observed both within and outside the gut. This study provides a framework for how whole-organism -omics approaches can be employed to understand the impact of antibiotics on the entire host.

IMPORTANCE We are just beginning to understand the unintended effects of antibiotics on our microbiomes and health. In this study, we aimed to define an approach by which one could obtain a comprehensive picture of (i) how antibiotics spatiotemporally impact commensal microbes throughout the gut and (ii) how these changes influence host chemistry throughout the body. We found that just a single dose of antibiotic altered host chemistry in a variety of organs and that microbiome alterations were not uniform throughout the gut. As technological advances increase the feasibility of whole-organism studies, we argue that using these approaches can provide further insight on both the wide-ranging effects of antibiotics on health and how to restore microbial communities to mitigate these effects.

OBSERVATION

Antibiotics, a cornerstone of modern medicine, have saved countless lives from infectious disease since their widespread introduction in the 1940s (1). However, in recent years, the collateral damage that antibiotics inflict on the commensal microbiome has become an increasing health concern. By ablating commensal microbes, antibiotics can allow secondary infections, such as Clostridioides difficile colitis (2) and candidiasis (3). Early life antibiotic exposure is particularly damaging (4) and has been associated with greater risk of asthma (5), eczema (6), inflammatory bowel disease (7), and obesity (8).

For antibiotic-induced microbiome perturbation to be transduced into altered disease risk, changes in the metabolism and/or molecular composition of the host must occur. However, little is known about how antibiotic-induced microbiome changes alter the overall chemical makeup of the body. Numerous studies have used fecal samples to assess microbiome composition (9–11) and stool metabolites (12–14) after antibiotic treatment. Here, we deployed novel mass spectrometry informatics and three-dimensional whole-organism data visualization approaches (15) to understand how a single dose of parenteral antibiotics affects the gastrointestinal (GI) tract microbiome and the local metabolome of every organ. We treated 10 mice with commonly prescribed parenteral antibiotics, ampicillin (AMP), a broad-spectrum beta-lactam, and vancomycin (VAN), a narrower-spectrum glycopeptide with activity against Gram-positive bacteria (16), and 10 mice were treated with phosphate-buffered saline (PBS; control) by intraperitoneal (i.p.) injection. One day and 6 days later, we dissected mice to sample 77 different body sites from 25 different organs to evaluate the whole-organism impact of antibiotics acutely and after microbiome recovery. We analyzed all samples by liquid chromatography-tandem mass spectrometry (LC-MS/MS) and the gastrointestinal tract samples (n = 38 per mouse) by 16S rRNA sequencing (see Fig. S1 in the supplemental material).

FIG S1

Overview of study design. Thirty-seven-week-old mice were ordered from Jackson Laboratories and randomly housed 2 or 3 mice per cage. After a 6-day rest period, mice were administered either 100 μl PBS, 100 mg/kg AMP in 100 μl PBS, or 100 mg/kg VAN in 100 μl PBS via intraperitoneal injection. Twenty-four hours later, 5 mice per group were dissected by following the protocol detailed in the methods. Six days postantibiotic administration, the remaining 5 mice per group were dissected by following the protocol detailed in the methods. Download FIG S1, TIF file, 1.4 MB.
Copyright © 2020 Vrbanac et al.

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

Antibiotic impacts on the microbiome.The greatest effects on the microbiome occurred in mice that had received antibiotics 24 h prior (AMP.d1 and VAN.d1). Compared to control mice, the centered log ratio (clr)-transformed proportion of Gram-positive bacteria was significantly reduced in the fecal pellet, certain sections of the stomach and colon of AMP.d1 and VAN.d1 mice, and certain sections of the cecum in AMP.d1 mice (Fig. S2).

FIG S2

Antibiotics reduce levels of Gram-positive bacteria. Asterisks indicate a significant difference in the centered log ratio (clr)-transformed proportion of Gram-positive bacteria compared to the control (Mann-Whitney U test, P < 0.05 after Benjamini-Hochberg correction for multiple comparisons). Error bars represent the 95% confidence interval (CI). Download FIG S2, TIF file, 1.1 MB.
Copyright © 2020 Vrbanac et al.

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

Compositionally, an increase in proteobacteria at the phylum level was observed in both AMP.d1 (cecum and fecal) and VAN.d1 mice (duodenum and stomach), differentially abundant sub-operational taxonomic units (sOTUs), also known as amplicon sequence variants or exact sequence variants, were seen in both AMP.d1 and VAN.d1 mice, and the level of one Clostridium sOTU was significantly lower in VAN.d6 mice (Fig. 1a and b). Antibiotic treatment also reduced microbiome alpha diversity (Fig. 1d). Specifically, Shannon diversity was reduced in the lower GI tract (cecum, colon, and fecal pellet) of AMP.d1 mice and in the upper GI tract (jejunum) of VAN.d1 and AMP.d1 mice. We next calculated unweighted UniFrac distance (17) for each group and compared effect size to those of controls down the GI tract (Fig. 1e). AMP.d1 mice had significantly altered beta diversity with high effect sizes, particularly in the lower GI tract, while beta diversity changes for VAN.d1 mice were only significant in a few GI sections. Effect size and diminished Shannon diversity were significantly correlated (Fig. 1c). AMP had a greater effect on the gut microbiota overall than VAN, consistent with its broader spectrum, recognizing that the biodistribution of parenteral VAN into the gut is poorly studied. Both AMP and VAN are renally excreted, and it is possible that some antibiotic was consumed orally via coprophagy through urine contacting fecal pellets.

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

Antibiotics affect the GI microbiome. (a) Stacked bar plots of the average relative abundance of bacterial phyla by organ in control mice and in mice 1 day after antibiotic treatment. (b) Count plot of significantly different sOTUs compared to the control; significance testing was performed with ANCOM (44). (c) Spearman correlation of changes in Shannon diversity from the control with pairwise effect size (e) for every noncontrol sample; the shaded area represents the 95% confidence interval (CI). (d) Shannon diversity down the GI tract. Asterisks indicate Shannon diversity significantly different from the control (Mann-Whitney U test with a Benjamini-Hochberg FDR control level of 0.1); error bars represent the 95% CI. (e) Pairwise effect size for each group compared to control samples, progressing down the GI tract. Significance was determined by PERMANOVA with a Benjamini-Hochberg FDR control level of 0.1 (48). The effect size is the PERMANOVA R2 value.

The global metabolome.To evaluate the global impact of antibiotics on metabolome composition, we calculated effect size using Bray-Curtis distance comparing each treatment group to controls and mapped the effect-size data onto a three-dimensional (3D) mouse model to provide a spatial visualization (Fig. 2a). For each body site/organ, averaging those with multiple samples (e.g., colon with 6 sections), only AMP.d1 and VAN.d1 showed significant differences in beta diversity (Fig. 2b). Both AMP.d1 and VAN.d1 had significantly different Bray-Curtis distance versus control mice in fecal samples. In peripheral organs, AMP.d1 mice differed significantly from the control in Bray-Curtis distance in the gallbladder, while VAN.d1 mice differed significantly in adrenal glands. We observed differentially abundant features in various organs compared to control mice (Fig. 2a, inset barplots), some putatively annotated with Global Natural Product Social Molecular Networking (GNPS) (18). Gallbladders of AMP.d1 mice had many differentially abundant features, including elevated bile acids (cholic, taurocholic, taurodeoxycholic, and cholenic acids). Uterine tissue of VAN.d1 and AMP.d1 mice had significantly higher levels of arachidonic acid, its metabolites, including 8-HETE and prostaglandin E2, and linoleic acid metabolites, such as 9-OxoODE. Significantly different features in adrenal glands of all antibiotic treatment groups compared to the control included a variety of elevated acylcarnitines and oxidized glutathione (GSSG). Glutathione is a key antioxidant in the adrenal cortex, where steroidogenesis generates reactive oxygen species (19), and GSSG is a biomarker of oxidative stress (20). VAN is nephrotoxic (21), and only kidneys of VAN.d1 mice had differentially abundant features, including elevated palmitoylcarnitine, a potential marker of renal toxicity in rats (22).

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

Organism-wide impact of antibiotics on the metabolome. (a) Pairwise effect size for each group compared to control samples, mapped onto body site using a 3D mouse model and a 2D illustration of the GI tract and percentage of significantly different metabolites compared to the control; significance was determined by dsFDR (47) (FDR control level, 0.1). (b) Pairwise effect size for each group compared to control samples. For body sites with multiple samples (i.e., colon cut into 6 sections), the metabolites were averaged for each mouse prior to calculating Bray-Curtis distance. Significance was determined by PERMANOVA with a Benjamini-Hochberg FDR control level of 0.1 (48), and effect size is the PERMANOVA R2 value.

Molecular networking analysis (18) of the raw metabolomics data revealed a cluster of features (molecules or metabolites) associated with AMP (Fig. S3). Many of these network features were present in our MS1 feature table and only in AMP.d1 mice, so these were summed to examine the collective distribution of AMP and its metabolites. Concentrated in the lower GI tract of AMP.d1 mice (Fig. S4A), the abundance of these AMP network features correlated with metabolome effect size (Bray-Curtis distance) and reduced Shannon diversity versus control mice (Fig. S3C and S4B), indicating that metabolome and microbiome effects were greatest in areas of high local AMP concentration. VAN (molecular weight, 1,449.2 g/mol) was not detected in our samples, limited by the LC-MS/MS MS1 scan range of m/z 100 to 1,500.

FIG S3

Molecular network of ampicillin and ampicillin network features. (a) GNPS molecular network displaying the molecular family in which ampicillin was annotated based on spectral library matching (blue) with unannotated analogs (green). Circles outlined in orange indicate a putative ampicillin analog. Neutral-charge chemical structures proposed are indicated. (b) MS/MS spectra for AMP. Putative fragments ions resulting from fragmentation of the proposed neutral-charge chemical structures are indicated using color. Download FIG S3, TIF file, 2.6 MB.
Copyright © 2020 Vrbanac et al.

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

FIG S4

Ampicillin network features correlate with changes in the metabolome and microbiome. (a) The sum of the peak area of putative ampicillin network features moving down the GI tract; error bars represent the 95% CI. (b) Pearson correlation of ampicillin network features with metabolome pairwise effect size to control (Bray-Curtis) for gut samples (stomach to fecal, down the GI) only from the ampicillin 1-day group. (c) Pearson correlation of ampicillin network features with the change in Shannon diversity (alpha diversity) relative to control samples for gut samples (stomach to fecal, down the GI) only from the AMP.d1 group. The shaded area represents 95% CI. Download FIG S4, TIF file, 1.7 MB.
Copyright © 2020 Vrbanac et al.

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

In the lower GI tract, differential abundance testing revealed an increase in small di- and tripeptides in AMP.d1 and VAN.d1 mice. To investigate this further, the metabolomics data were processed through PEAKs (23) to yield peptidomics data. For most sections of the lower GI tract and fecal pellet, the total number of peptides was significantly higher in day 1 postantibiotic mice than day 6 postantibiotic mice and untreated controls (Fig. S5A). These peptides aligned to several host proteins, including histones, which were particularly enriched in the AMP.d1 and VAN.d1 mice. Peptide fragments from histones H2A and H2B were significantly elevated in multiple sections of the lower GI (Fig. S5B). These elevated histones could simply be a marker of increased cell death, although histones also play an understudied role in innate immunity (reviewed in reference 24). H2A and H2B can contribute to host defense by direct antimicrobial activity (25, 26) or as components of neutrophil extracellular traps (27). H2A is also elevated in the ileum of chickens with experimental gut damage (28).

FIG S5

Small peptides are elevated in the lower GI tract day 1 postantibiotic exposure. (a) Sum of peptide spectral abundances along the lower GI tract comparing AMP.d1- and VAN.d1 (acute antibiotic)-treated mice to control mice and mice at day 6 posttreatment (nonacute). Significance indicates a P value of <0.05 by Mann-Whitney U test. (b) Abundances of peptides from histones, both H2A and H2B, were elevated in the lower GI tract of acute antibiotic-treated mice (AMP.d1 and VAN.d1) compared to other groups; significance testing was by dsFDR. Error bars represent the 95% CI. Download FIG S5, TIF file, 1.5 MB.
Copyright © 2020 Vrbanac et al.

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

To look at microbial and metabolite associations with histones, we ran a random forest regression on colon and cecum samples for both the microbiome and metabolome data, with each body site averaged for each individual mouse, for both histone H2A and histone H2B. Colon and cecum samples were evenly split between the train and test data sets. Regression performance for the microbiome data were modest for H2A (R2 = 0.59, P = 0.0033) and H2B (R2 = 0.41, P = 0.023), and the most important features in the regression model for H2A and H2B were mainly from the Clostridiales order and positively associated with histone abundance, except for one Anaeroplasma sOTU, the 6th most important feature for H2B, which was negatively correlated with both H2A and H2B abundance. For the metabolomics data, the regression performance was improved for H2A (R2 = 0.92, P = 7.3e−07) and H2B (R2 = 0.75, P = 0.00029), although annotated features were not among the top 20 most important features. The regression was repeated with the metabolome data filtered to only contain putatively annotated features for H2A (R2 = 0.81, P = 0.000068) and H2B (R2 = 0.737723, P = 0.000346). As expected, di- and tripeptides were among the top 10 annotated features for both H2A and H2B. Another top feature for H2A was phytomonic acid, a saturated fatty acid present in bacterial plasma membranes, which was positively associated with histone abundance; accumulation of phytomonic acid is associated with ethanol stress in Oenococcus oeni (29). Furthermore, 9,10-epoxy-12-octadecenoate, a metabolite that neutrophils synthesize during oxidative burst (30), was also among the top features for H2A and positively associated with histone abundance, hinting at neutrophil involvement.

Metabolome-microbiome interactions.For gut samples with paired 16S sequencing and metabolomics data, mmvec (31) was used to find cooccurrence probabilities between metabolites and microbes. Metabolites of the same general class, defined by putative GNPS annotations, had similar patterns of cooccurrence with microbes (Fig. S6A). To investigate these relationships in more detail, a multinomial regression model (Songbird [32]) was used to find sOTUs that were highly associated with AMP.d1 and VAN.d1 mice versus controls. Coprococcus and Sutterella were associated with AMP.d1 mice, and two family S24-7 bacteria were associated with VAN.d1 mice. Looking just at putatively annotated metabolites, all of these sOTUs had high cooccurrence probabilities with tri- or dipeptides and lowest cooccurrence probabilities with 13-docosenamide, a compound equally abundant in blanks (Fig. S6B). Coprococcus and Sutterella, the sOTUs associated with AMP.d1 mice, had high cooccurrence probabilities with several AMP network features and low cooccurrence probabilities with certain bile acids. Both the S24-7 sOTUs and Coprococcus had low cooccurrence probabilities with certain oligosaccharide- or sugar-related molecules, such as maltose. Marker gene sequencing data limits the exploration of the metabolic capacity of these bacteria to produce or consume these metabolites, although PICRUSt2 (33) predictions of functional potential for these two S24-7 and Coprococcus sOTUs included maltose O-acetyltransferase (KEGG entry K00661), an enzyme that acetylates maltose and other sugars, and a variety of other sugar transport systems for the Coprococcus sOTU (KEGG entries K10112, K10117, K10118, K10119, K03435, and K03436), suggesting the potential for the consumption of sugars. Differential abundance testing recapitulates some of these findings (elevated small peptides and AMP network features), although, in general, we did not observe alterations in the abundance of bile acids or sugar-related molecules. Studies with longer antibiotic courses report altered differential abundance of bile acids and sugars/carbohydrates (34, 35), and our cooccurrence data hint at a relationship between specific microbes enriched with antibiotic treatment and these metabolites.

FIG S6

Cooccurrence probabilities for microbes and metabolites cluster by metabolite class. (a) Log conditional cooccurrence probabilities for all microbes and putatively annotated metabolites. The microbial phylum is indicated on the left by row, and general metabolite class is indicated on the top by column (nonbiological metabolites are those that were highly abundant in blanks). (b) A multinomial regression model was used to find sOTUs that were most associated with the AMP.d1 mice (Coprococcus and Sutterella) and VAN.d1 mice (both family S24-7) compared to the control. Annotated metabolites that had the highest and lowest cooccurrence probabilities with these sOTUs are displayed. Download FIG S6, TIF file, 2.5 MB.
Copyright © 2020 Vrbanac et al.

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

Administering a single parenteral dose of either of the two common antibiotics did not drastically alter the microbiome or metabolome beyond recovery; in fact, microbiome alterations were primarily resolved by day 6. Nevertheless, a single dose of AMP or VAN still altered the metabolome at day 1 at a variety of body sites, both within and outside the GI tract. Presumably, extended courses of antibiotics, as typically prescribed in multiple daily doses, would exert even greater effects on the global metabolome, as previously seen in studies focused on the GI tract (36). Notably, we observed differential microbiome impacts with antibiotic type and by location in the GI tract. Successful restoration of microbial communities after antibiotic administration may require tailoring for antibiotic type and an understanding of which local GI communities are most impacted.

Organism-wide analyses of the scale we have undertaken remain ambitious and still encumber significant cost and manpower. However, as continual advances in technology have markedly reduced the costs of -omics platform analyses, this work provides a roadmap for understanding the real-time whole-organism-scale biological and chemical effects of antibiotic administration.

Methods.(i) Experimental design. Thirty-seven-week-old female C57BL/6 mice (Jackson Laboratory) were randomized to 3 or 2 mice per cage upon receipt and allowed to acclimatize for 6 days before antibiotic administration. AMP and VAN were administered by i.p. injection in doses of 100 mg/kg of body weight, dissolved in 100 μl PBS. Control mice were given 100 μl PBS via i.p. injection. Cohoused mice were administered the same antibiotic or PBS alone to avoid coprophagy across treatment groups. All mice, including controls, were dissected at both day 1 and day 6 time points. See Fig. S1 for a graphical representation of experimental design.

(ii) Sample collection and processing. Mice were dissected at both the day 1 and day 6 postantibiotic administration time points, alternating between treatment groups. The dissection protocol was as follows. Mice were sacrificed by CO2 asphyxiation, followed by cardiac puncture. Dissection was carried out under an open flame with forceps and scissors cleaned with 70% ethanol after the removal of each organ. Fresh fecal samples were collected prior to CO2 asphyxiation and organs removed in the following order: spleen, gallbladder, liver, pancreas, stomach, duodenum, jejunum, ileum, cecum, colon, kidneys, adrenal glands, urinary bladder, ovaries, uterine horns, uterine body, cervix, vagina, thymus, heart, lungs, esophagus, trachea, and brain. Organ subsections (e.g., lobes of the liver) were collected in separate tubes for a total of 77 samples per mouse.

Samples were flash frozen in 2-ml Eppendorf tubes in an isopropanol dry ice bath (−77°C) and stored at −80°C. All samples were weighed, and 10 μl of sterile water per milligram of tissue was added to each tube (i.e., 20 mg tissue, 200 μl water). For samples under 15 mg, 150 μl of water was added. Samples were homogenized in a Qiagen TissueLyser II using stainless steel beads.

(iii) Sample plating scheme. As the large number of samples required an extensive amount of LC-MS/MS instrument time, which was not available in a continuous block, we separated the samples into three organ groups to avoid issues with comparing treatment groups across metabolomics runs. These groups included “gut” samples (fecal, esophagus, stomach, duodenum, jejunum, ileum, cecum, colon, pancreas, liver, and gallbladder), “reproductive” samples (kidneys, adrenal glands, ovaries, vagina, cervix, uterus, uterine horns, and bladder), and “circulatory” samples (blood, spleen, heart, lung, trachea, thymus, and brain). To avoid plate effects biasing the results, samples were plated by cycling through treatment groups, resulting in every 96-well plate containing samples from multiple treatment groups. Both the metabolome and sequencing plates were prepped at the same time for the gut samples to avoid extra freeze-thaw cycles.

(iv) Metabolomics preparation and data acquisition. Five milligrams of tissue homogenate (50 μl) was added to 96-well plates containing 150 μl of 70% high-performance liquid chromatography (HPLC)-grade methanol (chilled) spiked with 5 μM sulfamethoxine (final concentration, ∼50% methanol). Plates were sealed, vortexed, and stored at 4°C overnight (12 h) for extraction. To remove insoluble material, plates were centrifuged at 2,000 rpm at 4°C for 10 min, and 130 μl from each well was aliquoted to a new 96-well plate. Plates were stored at −80°C, and, immediately prior to running on the LC-MS/MS instrument, were centrifuged at 2,000 rpm at 4°C for 10 min; 100 μl was aliquoted to a new 96-well plate. LC-MS/MS data acquisition was performed on a Vanquish ultrahigh-performance liquid chromatography (UPLC) system using a core-shell silica C18 column (50 by 2 mm, 1.7-μm particle size, 100-Å pore size; Kinetex, Phenomenex) coupled to a Q Exactive Orbitrap mass spectrometer (Thermo Fisher Scientific, Bremen, Germany). Five microliters of sample was injected and run at 0.5 ml/min on a gradient of solvent A (HPLC-grade water with 0.1% formic acid) and solvent B (HPLC-grade acetonitrile with 0.1% formic acid). The column was maintained at 40°C. The UPLC elution gradient ran for 12.5 min per sample: 5% B from 0 min to 1 min, a linear gradient of 5 to 100% B over 8 min, a hold at 100% B for 2 min, a return to 5% B over 0.5 min, and a hold at 5% B for 2 min to equilibrate the column for the next sample. The flow was directed into a heated electrospray ionization source operated in positive ionization mode with the following parameters: an auxiliary gas flow rate of 14 arbitrary units (a.u.), sweep gas flow rate of 3 a.u., sheath gas flow rate of 52 a.u., spray voltage of +3.5 kV, capillary temperature of 270°C, auxiliary gas heater temperature of 435°C, and S-Lens RF level of 50. The data-dependent acquisition mode was used to acquire the data in which MS1 scans from m/z 100 to 1,500 (scan rate, 7 Hz) were followed by an MS2 scan, specifically a product ion scan produced using stepped normalized collision energy higher-energy collisional dissociation, of the five most abundant ions from the prior MS1 scan.

(v) LC-MS/MS raw data processing. Raw data were uploaded to MassIVE (https://massive.ucsd.edu/), converted to .mzML files, imported into MZmine2 (37), and truncated at m/z 1500 and a 9.5-min retention time (RT). Due to some retention time shift while running the larger gut sample set, parameters are slightly modified versus non-gut sample sets, as indicated in parentheses. The parameters for the gut samples used in MZmine2 are as follows. Mass detection was performed with a noise threshold of 2.0e5 for MS1 and 2.0e3 for MS2 (1.0e3 for non-gut samples) in centroid mode, and chromatograms were built with a 0.05-min time span, 1.0e6 minimum height, and 10-ppm m/z tolerance. Chromatograms were deconvoluted with the baseline cutoff algorithm and a minimum peak height of 1.0e6, peak duration of 0.05 to 1.0 min, and baseline of 1.0e4. Isotope peak removal was performed with 15-ppm m/z tolerance, 0.3 retention time tolerance (0.05 RT for non-gut samples), and maximum charge of 4, and peaks were aligned with 10-ppm m/z tolerance, 75 weight m/z tolerance, and 0.4 RT tolerance (0.05 for non-gut samples). Gap filling was also performed with 10% intensity tolerance, 15-ppm m/z tolerance, and 0.3 RT tolerance (0.1 for non-gut samples). Peaks were also filtered to remove singletons found in only one sample. Both MS1 and MS2 feature tables were exported, and the “export for GNPS” feature was used to generate a .mgf file for GNPS (18). The signal intensities of the MS1 features were normalized (probabilistic quotient normalization) (38) to the sulfamethoxine internal standard.

(vi) Metabolite annotation. GNPS (18) was used to obtain putative annotations for the MS1 feature table via spectral library matching. The following parameters were used: no MS-Cluster, network was filtered to have edges with a cosine score above 0.7 and at least 4 matched peaks, and matches between network spectra and library spectra had a minimum of 6 matched peaks and a minimum cosine score of 0.7. For molecular networking with the raw files, the following parameters were used: MS-Cluster with a parent mass tolerance of 0.02 Da and an MS/MS fragment ion tolerance of 0.02 Da; consensus spectra with fewer than 5 spectra were discarded. The network filtered to have edges with a cosine score above 0.7 and at least 6 matched peaks, and matches between network spectra and library spectra had a minimum cosine score of 0.7 and at least 5 matched peaks.

(vii) Peptidomics data processing. LC-MS/MS .mzXML formatted files were loaded into PEAKS Studio 8.51 (23) for the identification of peptidic spectra matching the UniProt mouse protein database (www.uniprot.org; accessed 28 April 2018). Data were imported and refined according to PEAKS settings for Orbitrap instruments under the collision-induced dissociation fragmentation mode. Error tolerance parameters were set to 15-ppm parent mass error tolerance and 0.02-Da fragment mass error tolerance. The search settings included no added restriction enzymes as well as variable modifications for dehydration, acetylation (N terminal), oxidation (M), and formylation (N terminal). The maximum number of variable posttranslational modifications per peptide was set to 3. The label-free quantification mode was used, and the final false discovery rate (FDR) for peptide spectral matches was reported to be 2.9%. Quantification was normalized to the total ion chromatograph.

(viii) Microbiome sample processing. DNA was extracted from 100 μl of each gut sample using a MoBio PowerMag soil DNA isolation kit (Qiagen, Carlsbad, CA) by following the Earth Microbiome Project protocol (39). Samples were sequenced on an Illumina MiSeq.

(ix) Microbiome and metabolome data processing and analysis. Raw sequencing data were transferred to Qiita (40), where it was demultiplexed, trimmed to 150-bp reads, and denoised to sOTUs using Deblur (41). QIIME2 v2019.4 (42) was used for rarefaction (5,000 sequences per sample), to calculate beta diversity (unweighted UniFrac) and alpha diversity (Shannon) on microbiome data, and to calculate beta diversity (Bray-Curtis) on the metabolome data. Proportions of Gram-negative and Gram-positive bacteria were determined using Bugbase (https://doi.org/10.1101/133462). Bugbase requires an OTU table picked against the GreenGenes database, and closed reference OTU picking was performed with the GreenGenes database (v13_8) (43) at a 97% identity threshold. Using QIIME2 v2019.4 (42), the closed reference OTU table was filtered to exclude blanks and samples with fewer than 5,000 reads. The resulting .biom table was uploaded to BugBase for analysis. For ANCOMv (44) analysis for differentially abundant sOTUs, the nonrarefied feature table was filtered to remove samples with fewer than 5,000 reads/sample and all singleton sOTUs. GI body sites with multiple samples (i.e., colon, cut into 6 sections) were averaged for each mouse prior to performing ANCOM. PICRUSt2 (33) was used to predict functional abundances for sOTUs.

(x) Statistics and data visualization. SciPy (45) was used to run the Mann-Whitney U test for comparing proportions of Gram-positive bacteria, delta Shannon distance to control, and peptide sum down the GI tract. SciPy was also used to run Pearson and Spearman correlations. Statsmodels (46) was used to perform P value correction for multiple comparisons. As we did not flush the luminal contents, we observed some metabolome variation in GI sections depending on the volume of luminal contents in the section. For differential abundance and beta diversity analyses, these GI sections were averaged per body site per mouse (i.e., colon metabolome was the average of 6 colon sections per mouse). dsFDR (47) was used for differential abundance testing of metabolome data. Beta diversity significance was calculated with the QIIME2 (42) beta-group significance function using PERMANOVA and a Benjamini-Hochberg correction for multiple comparisons. Effect size was calculated by computing the PERMANOVA R2 values in python, analogous to the r-function of adonis (48). Figures were generated with matplotlib (49) and seaborn (50). The 3D mouse model was generated as described in Quinn et al. (15) from a magnetic resonance image of a mouse. Data visualization of metabolome effect size mapped onto the 3D mouse model was performed with ‘ili (https://ili.embl.de/) (51). Random forest regressions were performed with QIIME2 (42) using 100 trees, 20% of samples were used to test and 80% used to train, and 5 k-fold cross validations were performed.

Data availability.Raw data are available on MassIVE under the data sets MSV000082049, MSV000082157, and MSV000082048. Sequencing data are available under EBI accession number ERP121045.

ACKNOWLEDGMENTS

We thank Anupriya Tripathi and Yoshiki Vazquez-Baeza for their invaluable data analysis advice.

We acknowledge funding from the National Institutes of Health (NIH), grants DP1 AT010885, U01 AI124316, and P30 DK120515. R.H.M. is supported through a UCSD training grant from the NIH/NIDDK Gastroenterology Training Program (T32 DK007202). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

FOOTNOTES

    • Received April 16, 2020.
    • Accepted September 10, 2020.
  • Copyright © 2020 Vrbanac 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. Gaynes R
    . 2017. The discovery of penicillin—new insights after more than 75 years of clinical use. Emerg Infect Dis 23:849–853. doi:10.3201/eid2305.161556.
    OpenUrlCrossRef
  2. 2.↵
    1. Brown KA,
    2. Khanafer N,
    3. Daneman N,
    4. Fisman DN
    . 2013. Meta-analysis of antibiotics and the risk of community-associated Clostridium difficile Infection. Antimicrob Agents Chemother 57:2326–2332. doi:10.1128/AAC.02176-12.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. Xu J,
    2. Schwartz K,
    3. Bartoces M,
    4. Monsur J,
    5. Severson RK,
    6. Sobel JD
    . 2008. Effect of antibiotics on vulvovaginal candidiasis: a MetroNet study. J Am Board Fam Med 21:261–268. doi:10.3122/jabfm.2008.04.070169.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Bokulich NA,
    2. Chung J,
    3. Battaglia T,
    4. Henderson N,
    5. Jay M,
    6. Li H,
    7. Lieber AD,
    8. Wu F,
    9. Perez-Perez GI,
    10. Chen Y,
    11. Schweizer W,
    12. Zheng X,
    13. Contreras M,
    14. Dominguez-Bello MG,
    15. Blaser MJ
    . 2016. Antibiotics, birth mode, and diet shape microbiome maturation during early life. Sci Transl Med 8:343ra82. doi:10.1126/scitranslmed.aad7121.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Kozyrskyj AL,
    2. Ernst P,
    3. Becker AB
    . 2007. Increased risk of childhood asthma from antibiotic use in early life. Chest 131:1753–1759. doi:10.1378/chest.06-3008.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Hoskin-Parr L,
    2. Teyhan A,
    3. Blocker A,
    4. Henderson AJW
    . 2013. Antibiotic exposure in the first two years of life and development of asthma and other allergic diseases by 7.5 yr: a dose-dependent relationship. Pediatr Allergy Immunol 24:762–771. doi:10.1111/pai.12153.
    OpenUrlCrossRef
  7. 7.↵
    1. Kronman MP,
    2. Zaoutis TE,
    3. Haynes K,
    4. Feng R,
    5. Coffin SE
    . 2012. Antibiotic exposure and IBD development among children: a population-based cohort study. Pediatrics 130:e794–e803. doi:10.1542/peds.2011-3886.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Trasande L,
    2. Blustein J,
    3. Liu M,
    4. Corwin E,
    5. Cox LM,
    6. Blaser MJ
    . 2013. Infant antibiotic exposures and early-life body mass. Int J Obes 37:16–23. doi:10.1038/ijo.2012.132.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Dethlefsen L,
    2. Huse S,
    3. Sogin ML,
    4. Relman DA
    . 2008. The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16s rRNA sequencing. PLoS Biol 6:e280. doi:10.1371/journal.pbio.0060280.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Pérez-Cobas AE,
    2. Artacho A,
    3. Knecht H,
    4. Ferrús ML,
    5. Friedrichs A,
    6. Ott SJ,
    7. Moya A,
    8. Latorre A,
    9. Gosalbes MJ
    . 2013. Differential effects of antibiotic therapy on the structure and function of human gut microbiota. PLoS One 8:e80201. doi:10.1371/journal.pone.0080201.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Fouhy F,
    2. Guinane CM,
    3. Hussey S,
    4. Wall R,
    5. Ryan CA,
    6. Dempsey EM,
    7. Murphy B,
    8. Ross RP,
    9. Fitzgerald GF,
    10. Stanton C,
    11. Cotter PD
    . 2012. High-throughput sequencing reveals the incomplete, short-term recovery of infant gut microbiota following parenteral antibiotic treatment with ampicillin and gentamicin. Antimicrob Agents Chemother 56:5811–5820. doi:10.1128/AAC.00789-12.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Antunes LCM,
    2. Han J,
    3. Ferreira RBR,
    4. Lolić P,
    5. Borchers CH,
    6. Finlay BB
    . 2011. Effect of antibiotic treatment on the intestinal metabolome. Antimicrob Agents Chemother 55:1494–1503. doi:10.1128/AAC.01664-10.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    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 One 9:e101267. doi:10.1371/journal.pone.0101267.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Vrieze A,
    2. Out C,
    3. Fuentes S,
    4. Jonker L,
    5. Reuling I,
    6. Kootte RS,
    7. van Nood E,
    8. Holleman F,
    9. Knaapen M,
    10. Romijn JA,
    11. Soeters MR,
    12. Blaak EE,
    13. Dallinga-Thie GM,
    14. Reijnders D,
    15. Ackermans MT,
    16. Serlie MJ,
    17. Knop FK,
    18. Holst JJ,
    19. van der Ley C,
    20. Kema IP,
    21. Zoetendal EG,
    22. de Vos WM,
    23. Hoekstra JBL,
    24. Stroes ES,
    25. Groen AK,
    26. Nieuwdorp M
    . 2014. Impact of oral vancomycin on gut microbiota, bile acid metabolism, and insulin sensitivity. J Hepatol 60:824–831. doi:10.1016/j.jhep.2013.11.034.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Quinn RA,
    2. Melnik AV,
    3. Vrbanac A,
    4. Fu T,
    5. Patras KA,
    6. Christy MP,
    7. Bodai Z,
    8. Belda-Ferre P,
    9. Tripathi A,
    10. Chung LK,
    11. Downes M,
    12. Welch RD,
    13. Quinn M,
    14. Humphrey G,
    15. Panitchpakdi M,
    16. Weldon KC,
    17. Aksenov A,
    18. da Silva R,
    19. Avila-Pacheco J,
    20. Clish C,
    21. Bae S,
    22. Mallick H,
    23. Franzosa EA,
    24. Lloyd-Price J,
    25. Bussell R,
    26. Thron T,
    27. Nelson AT,
    28. Wang M,
    29. Leszczynski E,
    30. Vargas F,
    31. Gauglitz JM,
    32. Meehan MJ,
    33. Gentry E,
    34. Arthur TD,
    35. Komor AC,
    36. Poulsen O,
    37. Boland BS,
    38. Chang JT,
    39. Sandborn WJ,
    40. Lim M,
    41. Garg N,
    42. Lumeng JC,
    43. Xavier RJ,
    44. Kazmierczak BI,
    45. Jain R,
    46. Egan M,
    47. Rhee KE,
    48. Ferguson D,
    49. Raffatellu M,
    50. Vlamakis H,
    51. Haddad GG,
    52. Siegel D,
    53. Huttenhower C, et al
    . 2020. Global chemical effects of the microbiome include new bile-acid conjugations. Nature 579:123–129. doi:10.1038/s41586-020-2047-9.
    OpenUrlCrossRef
  16. 16.↵
    1. Wilhelm MP,
    2. Estes L
    . 1999. Symposium on antimicrobial agents. Part XII. Vancomycin. Mayo Clin Proc 74:928–935. doi:10.4065/74.9.928.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Lozupone C,
    2. Knight R
    . 2005. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 71:8228–8235. doi:10.1128/AEM.71.12.8228-8235.2005.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Wang M,
    2. Carver JJ,
    3. Phelan VV,
    4. Sanchez LM,
    5. Garg N,
    6. Peng Y,
    7. Nguyen DD,
    8. Watrous J,
    9. Kapono CA,
    10. Luzzatto-Knaan T,
    11. Porto C,
    12. Bouslimani A,
    13. Melnik AV,
    14. Meehan MJ,
    15. Liu W-T,
    16. Crüsemann M,
    17. Boudreau PD,
    18. Esquenazi E,
    19. Sandoval-Calderón M,
    20. Kersten RD,
    21. Pace LA,
    22. Quinn RA,
    23. Duncan KR,
    24. Hsu C-C,
    25. Floros DJ,
    26. Gavilan RG,
    27. Kleigrewe K,
    28. Northen T,
    29. Dutton RJ,
    30. Parrot D,
    31. Carlson EE,
    32. Aigle B,
    33. Michelsen CF,
    34. Jelsbak L,
    35. Sohlenkamp C,
    36. Pevzner P,
    37. Edlund A,
    38. McLean J,
    39. Piel J,
    40. Murphy BT,
    41. Gerwick L,
    42. Liaw C-C,
    43. Yang Y-L,
    44. Humpf H-U,
    45. Maansson M,
    46. Keyzers RA,
    47. Sims AC,
    48. Johnson AR,
    49. Sidebottom AM,
    50. Sedio BE, et al
    . 2016. Sharing and community curation of mass spectrometry data with Global Natural Products Social Molecular Networking. Nat Biotechnol 34:828–837. doi:10.1038/nbt.3597.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Prasad R,
    2. Kowalczyk JC,
    3. Meimaridou E,
    4. Storr HL,
    5. Metherell LA
    . 2014. Oxidative stress and adrenocortical insufficiency. J Endocrinol 221:R63–R73. doi:10.1530/JOE-13-0346.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Rossi R,
    2. Dalle-Donne I,
    3. Milzani A,
    4. Giustarini D
    . 2006. Oxidized forms of glutathione in peripheral blood as biomarkers of oxidative stress. Clin Chem 52:1406–1414. doi:10.1373/clinchem.2006.067793.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Bamgbola O
    . 2016. Review of vancomycin-induced renal toxicity: an update. Ther Adv Endocrinol Metab 7:136–147. doi:10.1177/2042018816638223.
    OpenUrlCrossRef
  22. 22.↵
    1. Zhao YY,
    2. Wang HL,
    3. Cheng XL,
    4. Wei F,
    5. Bai X,
    6. Lin RC,
    7. Vaziri ND
    . 2015. Metabolomics analysis reveals the association between lipid abnormalities and oxidative stress, inflammation, fibrosis, and Nrf2 dysfunction in aristolochic acid-induced nephropathy. Sci Rep 5:12936. doi:10.1038/srep12936.
    OpenUrlCrossRef
  23. 23.↵
    1. Zhang J,
    2. Xin L,
    3. Shan B,
    4. Chen W,
    5. Xie M,
    6. Yuen D,
    7. Zhang W,
    8. Zhang Z,
    9. Lajoie GA,
    10. Ma B
    . 2012. PEAKS DB: de novo sequencing assisted database search for sensitive and accurate peptide identification. Mol Cell Proteomics 11:M111.010587. doi:10.1074/mcp.M111.010587.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Parseghian MH,
    2. Luhrs KA
    . 2006. Beyond the walls of the nucleus: the role of histones in cellular signaling and innate immunity. Biochem Cell Biol 84:589–604. doi:10.1139/o06-082.
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.↵
    1. Patat SA,
    2. Carnegie RB,
    3. Kingsbury C,
    4. Gross PS,
    5. Chapman R,
    6. Schey KL
    . 2004. Antimicrobial activity of histones from hemocytes of the Pacific white shrimp. Eur J Biochem 271:4825–4833. doi:10.1111/j.1432-1033.2004.04448.x.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    1. Hoeksema M,
    2. Van Eijk M,
    3. Haagsman HP,
    4. Hartshorn KL
    . 2016. Histones as mediators of host defense, inflammation and thrombosis. Future Microbiol 11:441–453. doi:10.2217/fmb.15.151.
    OpenUrlCrossRef
  27. 27.↵
    1. Brinkmann V,
    2. Reichard U,
    3. Goosmann C,
    4. Fauler B,
    5. Uhlemann Y,
    6. Weiss DS,
    7. Weinrauch Y,
    8. Zychlinsky A
    . 2004. Neutrophil extracellular traps kill bacteria. Science 303:1532–1535. doi:10.1126/science.1092385.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. De Meyer F,
    2. Eeckhaut V,
    3. Ducatelle R,
    4. Dhaenens M,
    5. Daled S,
    6. Dedeurwaerder A,
    7. De Gussem M,
    8. Haesebrouck F,
    9. Deforce D,
    10. Van Immerseel F
    . 2019. Host intestinal biomarker identification in a gut leakage model in broilers. Vet Res 50:46. doi:10.1186/s13567-019-0663-x.
    OpenUrlCrossRef
  29. 29.↵
    1. Teixeira H,
    2. Gonçalves MG,
    3. Rozès N,
    4. Ramos A,
    5. San Romão MV
    . 2002. Lactobacillic acid accumulation in the plasma membrane of Oenococcus oeni: a response to ethanol stress? Microb Ecol 43:146–153. doi:10.1007/s00248-001-0036-6.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Hayakawa M,
    2. Sugiyama S,
    3. Takamura T,
    4. Yokoo K,
    5. Iwata M,
    6. Suzuki K,
    7. Taki F,
    8. Takahashi S,
    9. Ozawa T
    . 1986. Neutrophils biosynthesize leukotoxin, 9, 10-epoxy-12-octadecenoate. Biochem Biophys Res Commun 137:424–430. doi:10.1016/0006-291x(86)91227-1.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Morton JT,
    2. Aksenov AA,
    3. Nothias LF,
    4. Foulds JR,
    5. Quinn RA,
    6. Badri MH,
    7. Swenson TL,
    8. Van Goethem MW,
    9. Northen TR,
    10. Vazquez-Baeza Y,
    11. Wang M,
    12. Bokulich NA,
    13. Watters A,
    14. Song SJ,
    15. Bonneau R,
    16. Dorrestein PC,
    17. Knight R
    . 2019. Learning representations of microbe–metabolite interactions. Nat Methods 16:1306–1314. doi:10.1038/s41592-019-0616-3.
    OpenUrlCrossRef
  32. 32.↵
    1. Morton JT,
    2. Marotz C,
    3. Washburne A,
    4. Silverman J,
    5. Zaramela LS,
    6. Edlund A,
    7. Zengler K,
    8. Knight R
    . 2019. Establishing microbial composition measurement standards with reference frames. Nat Commun 10:2719. doi:10.1038/s41467-019-10656-5.
    OpenUrlCrossRef
  33. 33.↵
    1. Douglas GM,
    2. Maffei VJ,
    3. Zaneveld JR,
    4. Yurgel SN,
    5. Brown JR,
    6. Taylor CM,
    7. Huttenhower C,
    8. Langille MGI
    . 2020. PICRUSt2 for prediction of metagenome functions. Nat Biotechnol 38:685–688. doi:10.1038/s41587-020-0548-6.
    OpenUrlCrossRef
  34. 34.↵
    1. Theriot CM,
    2. Koenigsknecht MJ,
    3. Carlson PE,
    4. Hatton GE,
    5. Nelson AM,
    6. Li B,
    7. Huffnagle GB,
    8. Z Li J,
    9. Young VB,
    10. Young VB
    . 2014. Antibiotic-induced shifts in the mouse gut microbiome and metabolome increase susceptibility to Clostridium difficile infection. Nat Commun 5:3114. doi:10.1038/ncomms4114.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Gutierrez D,
    2. Weinstock A,
    3. Antharam VC,
    4. Gu H,
    5. Jasbi P,
    6. Shi X,
    7. Dirks B,
    8. Krajmalnik-Brown R,
    9. Maldonado J,
    10. Guinan J,
    11. Thangamani S
    . 2020. Antibiotic-induced gut metabolome and microbiome alterations increase the susceptibility to Candida albicans colonization in the gastrointestinal tract. FEMS Microbiol Ecol 96:fiz187. doi:10.1093/femsec/fiz187.
    OpenUrlCrossRef
  36. 36.↵
    1. Theriot CM,
    2. Bowman AA,
    3. Young VB
    . 2016. Antibiotic-induced alterations of the gut microbiota alter secondary bile acid production and allow for Clostridium difficile spore germination and outgrowth in the large intestine. mSphere 1:e00045-15. doi:10.1128/mSphere.00045-15.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Pluskal T,
    2. Castillo S,
    3. Villar-Briones A,
    4. Orešič M
    . 2010. MZmine 2: modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics 11:395. doi:10.1186/1471-2105-11-395.
    OpenUrlCrossRefPubMed
  38. 38.↵
    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 in1H NMR metabonomics. Anal Chem 78:4281–4290. doi:10.1021/ac051632c.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Marotz C,
    2. Amir A,
    3. Humphrey G,
    4. Gaffney J,
    5. Gogul G,
    6. Knight R
    . 2017. DNA extraction for streamlined metagenomics of diverse environmental samples. Biotechniques 62:290–293. doi:10.2144/000114559.
    OpenUrlCrossRef
  40. 40.↵
    1. Gonzalez A,
    2. Navas-Molina JA,
    3. Kosciolek T,
    4. McDonald D,
    5. Vázquez-Baeza Y,
    6. Ackermann G,
    7. DeReus J,
    8. Janssen S,
    9. Swafford AD,
    10. Orchanian SB,
    11. Sanders JG,
    12. Shorenstein J,
    13. Holste H,
    14. Petrus S,
    15. Robbins-Pianka A,
    16. Brislawn CJ,
    17. Wang M,
    18. Rideout JR,
    19. Bolyen E,
    20. Dillon M,
    21. Caporaso JG,
    22. Dorrestein PC,
    23. Knight R
    . 2018. Qiita: rapid, web-enabled microbiome meta-analysis. Nat Methods 15:796–798. doi:10.1038/s41592-018-0141-9.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Amir A,
    2. McDonald D,
    3. Navas-Molina JA,
    4. Kopylova E,
    5. Morton JT,
    6. Zech Xu Z,
    7. Kightley EP,
    8. Thompson LR,
    9. Hyde ER,
    10. Gonzalez A,
    11. Knight R
    . 2017. Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems 2:e00191-16. doi:10.1128/mSystems.00191-16.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Bolyen E,
    2. Rideout JR,
    3. Dillon MR,
    4. Bokulich NA,
    5. Abnet CC,
    6. Al-Ghalith GA,
    7. Alexander H,
    8. Alm EJ,
    9. Arumugam M,
    10. Asnicar F,
    11. Bai Y,
    12. Bisanz JE,
    13. Bittinger K,
    14. Brejnrod A,
    15. Brislawn CJ,
    16. Brown CT,
    17. Callahan BJ,
    18. Caraballo-Rodríguez AM,
    19. Chase J,
    20. Cope EK,
    21. Da Silva R,
    22. Diener C,
    23. Dorrestein PC,
    24. Douglas GM,
    25. Durall DM,
    26. Duvallet C,
    27. Edwardson CF,
    28. Ernst M,
    29. Estaki M,
    30. Fouquier J,
    31. Gauglitz JM,
    32. Gibbons SM,
    33. Gibson DL,
    34. Gonzalez A,
    35. Gorlick K,
    36. Guo J,
    37. Hillmann B,
    38. Holmes S,
    39. Holste H,
    40. Huttenhower C,
    41. Huttley GA,
    42. Janssen S,
    43. Jarmusch AK,
    44. Jiang L,
    45. Kaehler BD,
    46. Bin Kang K,
    47. Keefe CR,
    48. Keim P,
    49. Kelley ST,
    50. Knights D, et al
    . 2019. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol 37:852–857. doi:10.1038/s41587-019-0209-9.
    OpenUrlCrossRefPubMed
  43. 43.↵
    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 J 6:610–618. doi:10.1038/ismej.2011.139.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Mandal S,
    2. Van Treuren W,
    3. White RA,
    4. Eggesbø M,
    5. Knight R,
    6. Peddada SD
    . 2015. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis 26:27663. doi:10.3402/mehd.v26.27663.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Virtanen P,
    2. Gommers R,
    3. Oliphant TE,
    4. Haberland M,
    5. Reddy T,
    6. Cournapeau D,
    7. Burovski E,
    8. Peterson P,
    9. Weckesser W,
    10. Bright J,
    11. van der Walt SJ,
    12. Brett M,
    13. Wilson J,
    14. Millman KJ,
    15. Mayorov N,
    16. Nelson ARJ,
    17. Jones E,
    18. Kern R,
    19. Larson E,
    20. Carey CJ,
    21. Polat İ,
    22. Feng Y,
    23. Moore EW,
    24. VanderPlas J,
    25. Laxalde D,
    26. Perktold J,
    27. Cimrman R,
    28. Henriksen I,
    29. Quintero EA,
    30. Harris CR,
    31. Archibald AM,
    32. Ribeiro AH,
    33. Pedregosa F,
    34. van Mulbregt P,
    35. Vijaykumar A,
    36. Bardelli A,
    37. Pietro Rothberg A,
    38. Hilboll A,
    39. Kloeckner A,
    40. Scopatz A,
    41. Lee A,
    42. Rokem A,
    43. Woods CN,
    44. Fulton C,
    45. Masson C,
    46. Häggström C,
    47. Fitzgerald C,
    48. Nicholson DA,
    49. Hagen DR,
    50. Pasechnik DV
    , SciPy 1.0 Contributors. 2020. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Methods 17:261–272. doi:10.1038/s41592-019-0686-2.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Seabold S,
    2. Perktold J
    . 2010. Statsmodels: econometric and statistical modeling with python. Proceedings of the 9th Python in Science Conference. Scientific Computing with Python, Austin, TX. doi:10.25080/Majora-92bf1922-011.
    OpenUrlCrossRef
  47. 47.↵
    1. Jiang L,
    2. Amir A,
    3. Morton JT,
    4. Heller R,
    5. Arias-Castro E,
    6. Knight R
    . 2017. Discrete false-discovery rate improves identification of differentially abundant microbes. mSystems 2:e00092-17. doi:10.1128/mSystems.00092-17.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Anderson MJ
    . 2001. A new method for non-parametric multivariate analysis of variance. Austral Ecol 26:32–46. doi:10.1046/j.1442-9993.2001.01070.x.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Hunter JD
    . 2007. Matplotlib: a 2D graphics environment. Comput Sci Eng 9:90–95. doi:10.1109/MCSE.2007.55.
    OpenUrlCrossRef
  50. 50.↵
    1. Waskom M,
    2. Botvinnik O,
    3. Ostblom J,
    4. Lukauskas S,
    5. Hobson P,
    6. MaozGelbart Gemperline DC,
    7. Augspurger T,
    8. Halchenko Y,
    9. Cole JB,
    10. Warmenhoven J,
    11. de Ruiter J,
    12. Pye C,
    13. Hoyer S,
    14. Vanderplas J,
    15. Villalba S,
    16. Kunter G,
    17. Quintero E,
    18. Bachant P,
    19. Martin M,
    20. Meyer K,
    21. Swain C,
    22. Miles A,
    23. Brunner T,
    24. O’Kane D,
    25. Yarkoni T,
    26. Williams ML,
    27. Evans C
    . 24 January 2020. mwaskom/seaborn: v0.10.0. https://zenodo.org/record/3629446#.X2Oca4t7mUk.
  51. 51.↵
    1. Protsyuk I,
    2. Melnik AV,
    3. Nothias LF,
    4. Rappez L,
    5. Phapale P,
    6. Aksenov AA,
    7. Bouslimani A,
    8. Ryazanov S,
    9. Dorrestein PC,
    10. Alexandrov T
    . 2018. 3D molecular cartography using LC–MS facilitated by Optimus and 'ili software. Nat Protoc 13:134–154. doi:10.1038/nprot.2017.122.
    OpenUrlCrossRef
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Evaluating Organism-Wide Changes in the Metabolome and Microbiome following a Single Dose of Antibiotic
Alison Vrbanac, Kathryn A. Patras, Alan K. Jarmusch, Robert H. Mills, Samuel R. Shing, Robert A. Quinn, Fernando Vargas, David J. Gonzalez, Pieter C. Dorrestein, Rob Knight, Victor Nizet
mSystems Oct 2020, 5 (5) e00340-20; DOI: 10.1128/mSystems.00340-20

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 mSystems 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.
Evaluating Organism-Wide Changes in the Metabolome and Microbiome following a Single Dose of Antibiotic
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Evaluating Organism-Wide Changes in the Metabolome and Microbiome following a Single Dose of Antibiotic
Alison Vrbanac, Kathryn A. Patras, Alan K. Jarmusch, Robert H. Mills, Samuel R. Shing, Robert A. Quinn, Fernando Vargas, David J. Gonzalez, Pieter C. Dorrestein, Rob Knight, Victor Nizet
mSystems Oct 2020, 5 (5) e00340-20; DOI: 10.1128/mSystems.00340-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • OBSERVATION
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

3D data visualization
antibiotics
mass spectrometry
metabolome
microbiome

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • 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 #mSystemsJ

@ASMicrobiology

       

 

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

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

Online ISSN: 2379-5077