Characterization of Microbial Dysbiosis and Metabolomic Changes in Dogs with Acute Diarrhea

Abstract

Limited information is available regarding the metabolic consequences of intestinal dysbiosis in dogs with acute onset of diarrhea. The aim of this study was to evaluate the fecal microbiome, fecal concentrations of short-chain fatty acids (SCFAs), as well as serum and urine metabolites in healthy dogs (n=13) and dogs with acute diarrhea (n=13). The fecal microbiome, SCFAs, and serum/urine metabolite profiles were characterized by 454-pyrosequencing of the 16S rRNA genes, GC/MS, and untargeted and targeted metabolomics approach using UPLC/MS and HPLC/MS, respectively. Significantly lower bacterial diversity was observed in dogs with acute diarrhea in regards to species richness, chao1, and Shannon index (p=0.0218, 0.0176, and 0.0033; respectively). Dogs with acute diarrhea had significantly different microbial communities compared to healthy dogs (unweighted Unifrac distances, ANOSIM p=0.0040). While Bacteroidetes, Faecalibacterium, and an unclassified genus within Ruminococcaceae were underrepresented, the genus Clostridium was overrepresented in dogs with acute diarrhea. Concentrations of fecal propionic acid were significantly decreased in acute diarrhea (p=0.0033), and were correlated to a decrease in Faecalibacterium (ρ=0.6725, p=0.0332). The predicted functional gene content of the microbiome (PICRUSt) revealed overrepresentations of genes for transposase enzymes as well as methyl accepting chemotaxis proteins in acute diarrhea. Serum concentrations of kynurenic acid and urine concentrations of 2-methyl-1H-indole and 5-Methoxy-1H-indole-3-carbaldehyde were significantly decreased in acute diarrhea (p=0.0048, 0.0185, and 0.0330, respectively). These results demonstrate that the fecal dysbiosis present in acute diarrhea is associated with altered systemic metabolic states.

About the Authors:

Blake C. Guard

Affiliation: Gastrointestinal Laboratory, Department of Small Animal Clinical Sciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, Texas, United States of America

James W. Barr

Affiliation: Emergency and Critical Care, Department of Small Animal Clinical Sciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, Texas, United States of America

Lavanya Reddivari

Affiliation: Department of Plant Science, Penn State University, University Park, Pennsylvania, United States of America

Cory Klemashevich

Affiliation: Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas, United States of America

Arul Jayaraman

Affiliations Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas, United States of America, Department of Microbial Pathogenesis and Immunology, Texas A&M Health Science Center, College Station, Texas, United States of America

Jörg M. Steiner

Affiliation: Gastrointestinal Laboratory, Department of Small Animal Clinical Sciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, Texas, United States of America

Jairam Vanamala

Affiliations Department of Food Science, Penn State University, University Park, Pennsylvania, United States of America, The Penn State Hershey Cancer Institute, Hershey, Pennsylvania, United States of America

Jan S. Suchodolski

* E-mail: [email protected]

Affiliation: Gastrointestinal Laboratory, Department of Small Animal Clinical Sciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, Texas, United States of America

Introduction

Recent molecular studies have greatly increased our knowledge about the microbiota in the gastrointestinal tract (GIT) of dogs, mice, and humans [1–4]. The GIT microbiota plays an important role in host health by stimulating the immune system, influencing gut structure, aiding in the defense against pathogens, and providing nutritional benefits to the host (e.g., production of SCFAs) [5–10]. Despite recent advances in characterizing microbial communities using sequencing technology, there remains a rudimentary understanding of the complex interactions that occur between the host and intestinal microbes and their metabolic end-products. The analysis of host and bacterial metabolites may give additional insights into the pathophysiology of gastrointestinal diseases, including acute diarrhea. Metabolomics is a relatively new field that aims to characterize qualitatively and/or quantitatively the presence of small molecules in biological samples [11]. This approach may enhance our understanding of the host-microbe interactions, as well as the metabolic pathways that are involved in health and disease. Initial metagenomic and metabolomic studies performed in humans and animal models suggest that the metabolites derived from diverse microbial communities may have a direct role in health and disease [12]. Therefore, more in-depth studies are needed to understand the relationships between the microbiome and the host [8]. Phylogenetic changes in intestinal microbiota have been previously described in dogs with acute and chronic GI disease [13–15]. A previous study by our group has compared the fecal microbiome between healthy dogs, dogs with IBD, and dogs with acute diarrhea, but reported only phylogenetic data based on 16S rRNA gene sequences and it is currently unknown whether these microbiome changes are associated with metabolic changes in the host [14]. Therefore, the aim of the present study was to evaluate a new cohort of dogs with acute diarrhea to confirm the previous findings, and also to elucidate whether microbiome changes based on 16S rRNA genes are associated with metabolic and functional changes. We profiled the fecal microbiome using 16S rRNA sequencing, measured and correlated fecal metabolic end products (i.e., SCFAs) with bacterial groups, and inferred metagenomics using PICRUSt [16]. Furthermore, this study used Ultra Performance Liquid Chromatography-Mass Spectrometry (UPLC-MS) in an untargeted metabolomics approach coupled with High Performance Liquid Chromatography-Mass Spectrometry (HPLC-MS) in a targeted metabolomics approach to characterize and confirm metabolomic alterations in the serum and urine of dogs with acute diarrhea.

Materials and Methods

Animal enrollment and sample collection

Naturally passed feces, serum samples (collected by venipuncture), and urine samples (collected by cystocentesis) were obtained from healthy dogs as well as dogs with acute diarrhea (AD). Dogs with AD were further classified as having non-hemorrhagic diarrhea (NHD) or hemorrhagic diarrhea (AHD) (Table 1). Feces were refrigerated immediately after collection, transferred within a few hours to a -80°C freezer, and stored frozen until processing for DNA extraction. Serum and urine samples were aliquoted and stored frozen at -80°C until processing. Owners provided written consent for their dogs to be used in this study. The collection of feces, serum, and urine was approved by the Texas A&M University Institutional Animal Care and Use Committee (IACUC): Protocol Number; 2012–101. None of the healthy dogs or dogs with acute diarrhea had been used in a previous study [14].

[Figure omitted. See PDF.]

Table 1. Summary of basic characteristics and alpha diversity measures.

https://doi.org/10.1371/journal.pone.0127259.t001

The control group consisted of 13 healthy pet dogs (Table 1, S1 Table). All dogs were privately owned, lived in various home environments, and were fed a variety of commercial diets. None of the dogs had a history of gastrointestinal signs or administration of antibiotics for at least a month prior to collection of fecal samples. All healthy dogs lived in Texas, USA.

The diseased group consisted of 13 dogs in total with acute diarrhea (6 dogs with NHD and 7 dogs with AHD) that presented to the Veterinary Medical Teaching Hospital at Texas A&M University with acute, non-hemorrhagic or hemorrhagic diarrhea (defined as duration of diarrhea <3 days) (Table 1, S1 Table). None of the dogs had a previous history of GI signs or had received antibiotics within the previous three months. Diagnostic evaluation included a complete blood count (CBC), serum chemistry profiles (SIRRUS Clinical Chemistry Analyzer), serum concentrations of canine trypsin-like immunoreactivity (cTLI), serum concentrations of canine pancreatic lipase immunocreactivity (cPLI), serum concentrations of cobalamin and folate (Immulite 2000 Vitamin B12, Folic Acid, Siemens Medical Solutions Diagnostics), and serum concentrations of C-reactive protein (CRP; Phase Range Canine C-reactive Protein Assay, Tridelta Development Ltd), and urine analysis. None of the dogs with acute diarrhea had been readmitted to the clinic at Texas A&M University after one month of initial presentation and for the purpose of this study their clinical signs were considered resolved.

DNA isolation

100 mg of feces were aliquoted into a sterile 1.7 ml tube (Microtube, Sarstedt AG & Co, Nümbrecht, Germany) containing 150 μl of 0.1 mm zirconia-silica beads and 100 μl of 0.5 mm zirconia-silica beads (BioSpec Products Inc., Barlesville, OK, USA). Samples were then homogenized (FastPrep-24, MP Biomedicals, USA) for a duration of 1 minute at a speed of 4 m/s. DNA was then extracted with the ZR fecal DNA Mini Prep kit following the manufacturer’s instructions (Zymo Research, Irvine CA, USA).

454-Pyrosequencing

Bacterial tag-encoded FLX-titanium amplicon pyrosequencing was performed targeting the V4–V6 region of the 16S rRNA gene using forward and reverse primers: 530F (5’-GTGCCAGCMGCNGCGG-3’) and 1100R (5’-GGGTTNCGNTCGTTG-3’), respectively [17]. Raw sequence data were screened, trimmed, de-noised, filtered, and depleted of chimeras using the QIIME v1.7 open-source pipeline [18]. Operational taxonomic units (OTUs) were assigned based on at least 97% sequence similarity using QIIME. The sequences were deposited in the Sequence Read Archive under the following accession number: SRP040310.

Quantitative PCR (qPCR)

To evaluate specific bacterial groups of interest to intestinal health (i.e., Lactobacillus and Bifidobacterium) and bacterial species with potential pathogenic roles (i.e., Escherichia coli, and Clostridium perfringens) qPCR was used as described previously [19–21] and primers/probe can be viewed in the S4 Table. Quantitative PCR reactions were performed using two reaction chemistries. For a subset of assays SYBR-green based reaction mixtures were used, with a total reaction volume of 10 μl. The final mix contained 5 μl SsoFast EvaGreen supermix (Bio-Rad Laboratories, CA, USA), 0.4 μl each of a forward and reverse primer (final concentration: 400 nM), 2.6 μl of high quality PCR water, and 2 μl of normalized DNA (final concentration: 5 ng/μl). Conditions for PCR were as follows: initial denaturation at 98°C for 2 min, then 40 cycles with denaturation at 98°C for 3 sec and annealing (S4 Table) for 3 sec. Post-amplification, a melt curve analysis was performed using these conditions: 95°C for 1 min, 55°C for 1 min, and increasing incremental steps of 0.5°C for 80 cycles for 5 sec each. All samples were run in duplicate fashion. TaqMan based reaction mixtures were used in a total reaction volume of 10 μl. The final mix contained 5 μl TaqMan Fast Universal PCR master mix (Life Technologies, NY, USA), 0.4 μl of a forward and reverse primer (final concentration: 400 nM), 2 μl of high quality PCR water, and 2 μl of normalized DNA (final concentration: 5 ng/μl). Conditions for PCR were as follows: initial denaturation at 95°C for 20 sec then 40 cycles with denaturation at 95°C and annealing (S4 Table) for 3 sec. Post-amplification, a melt curve analysis was performed using these conditions: 95°C for 1 min, 55°C for 1 min, and increasing incremental steps of 0.5°C for 80 cycles for 5 sec each. All samples were run in duplicate fashion.

Measurement of short-chain fatty acids

Analysis of short-chain fatty acids (SCFAs; i.e., acetate, propionate, and butyrate), and branched chain fatty acids (BCFAs; i.e., isobutyrate, isovalerate, and valerate) in feces was performed on a subgroup of dogs, as not all dogs had sufficient quantity of feces available. The SCFA and BCFA were measured using a dilution gas chromatography-mass spectrometry (GC-MS) assay as previously described [22] with some modifications. Briefly, the fecal samples were weighed and diluted 1:5 in extraction solution (2N hydrochloric acid). After homogenization using a multi-tube vortexer for 30 min at room temperature, fecal suspensions were centrifuged for 20 min at 2,100 x g and 4°C. Supernatants were then collected using serum filters (Fisherbrand serum filter system, Fisher Scientific Inc, Pittsburgh, Pa). From each sample, 500 μl of supernatant were mixed with 10 μl of internal standard (200 mM heptadeuterated butyric acid) and extracted using a C18 solid phase extraction column (Sep-Pak C18 1 cc Vac Cartridge, Waters Corporation, Milford, MA). Samples were derivatized using N-tert-butyldimethylsilyl-N-methyltrifluoroacetamide (MTBSTFA) at room temperature for 60 min. A gas chromatographer (Agilent 6890N, Agilent Technologies Inc, Santa Clara, CA) coupled with a mass spectrometer (Agilent 5975C, Agilent Technologies Inc, Santa Clara, CA) was used for chromatographic separation and quantification of the derivatized samples. Separation was achieved using a DB-1ms capillary column (30 m x 0.25 mm, 0.25 μm film thickness, Agilent Technologies Inc, Santa Clara, CA). The GC temperature program was as follows: 40°C held for 0.1 min, increased to 70°C by 5°C/min, 70°C held for 3.5 min, increased to 160°C by 20°C/min, and finally increased to 280°C for 3 min by 35°C/min. The total run time was 20.53 min. Solvent delay was 5 min. The mass spectrometer was operated in electron impact positive-ion mode with selective ion monitoring at mass-to-charge ratios (M/Z) of 117 (acetate), 131 (propionate), 145 (butyrate and isobutyrate), 152 (deuterated butyrate; internal standard), and 159 (valerate and isovalerate). The estimate of the relative concentrations was based on the ratio of the area under the curve of the internal standard for each fatty acid.

Measurement of dry weight

For determination of fecal dry weight, a 100 mg of feces from each sample was aliquoted into a sterile 1.7 ml serum tube (Microtube, Sarstedt AG & Co, Nümbrecht, Germany). This aliquot was solely used in the determination of dry weight and thus was weighed and dried at 105°C in an oven (Symphony Gravity Convection Oven, VWR) overnight. The percent of dry weight was calculated and SCFA/BCFA concentrations were adjusted (measured in μmol/g of dry feces for all SCFAs/BCFAs) to normalize by dry weight.

Untargeted metabolite detection in serum

Serum samples were thawed on ice and diluted 1:5 in cold methanol. The mixture was vortexed and incubated at -20°C for two hours to precipitate protein. The solution was centrifuged at 14,000 x g for 20 minutes at 4°C, and then the supernatant was collected and transferred to a Waters 96 well autosampler plate. 1 μl injections were performed on a Waters Acquity UPLC system. Separation was performed using a Waters Acquity UPLC C8 column (1.8 μm, 1.0 x 100 mm), using a gradient from solvent A (95% water, 5% methanol, 0.1% formic acid) to solvent B (95% methanol, 5% water, 0.1% formic acid). Injections were made in 100% A, which was held for 0.1 min, ramped to 40% B in 0.9 min, to 70% B over two min, and to 100% B over 8 min. Mobile phase was held at 100% B for 6 min, returned to starting conditions over 0.1 min, and allowed to re-equilibrate for 5.9 min. Flow rate was constant at 140 μl/min for the duration of the run. The column was held at 50°C and samples were held at 10°C. Column Eluent was infused into a Waters Xevo G2 Q-TOF MS fitted with an electrospray source. Data was collected in positive ion mode, scanning from 50–1200 at a rate of 0.2 sec per scan. Collision energy was set to 6 Volts for MS mode, and ramped from 15 to 30 volts for MSE mode. Calibration was performed prior to sample analysis via infusion of sodium formate solution, with mass accuracy within 1 ppm. The capillary voltage was held at 2200 volts, the source temp at 150°C, and the desolvation temperature at 350°C at a nitrogen desolvation gas flow rate of 800 L/hr.

Untargeted metabolite detection in urine

Analysis of metabolites was performed on a subgroup of dogs, as urine was unable to be collected from all dogs. Urine samples were diluted 1:2 in methanol, centrifuged, and transferred to a Waters 96 well autosampler plate. The peak areas were first normalized based on total extracted ion signal and then normalized based on urine creatinine concentrations measured by a SIRRUS Clinical Chemistry Analyzer. A volume of 5 μl was injected into a Waters Acquity UPLC system. Separation was performed using a Waters Acquity UPLC T3 column (1.8 μm, 1.0 x 100 mm), using a gradient from solvent A (water, 0.1% formic acid) to solvent B (acetonitrile, 0.1% formic acid). Injections were made in 100% A, which was held for 1 min, a 12 min linear gradient to 95% B was applied, and held at 95% B for 3 min, returned to starting conditions over 0.05 min, and allowed to re-equilibrate for 3.95 min. Flow rate was constant at 200 μl/min for the duration of the run. The column was held at 50°C and the samples were held at 5°C. Column eluent was infused into a Waters Xevo G2 Q-TOF MS fitted with an electrospray source. Data was collected in positive ion mode, scanning from 50–1200 at a rate of 0.2 sec per scan, alternating between MS and MSE mode. Collision energy was set to 6 Volts for MS mode, and ramped from 15 to 30 volts for MSE mode. Calibration was performed prior to sample analysis via infusion of sodium formate solution, with mass accuracy within 1 ppm. The capillary voltage was held at 2200V, the source temp at 150°C, and the desolvation temperature at 350°C at a nitrogen desolvation gas flow rate of 800 L/hr.

Targeted metabolite detection in serum and urine

Standards for Kynurenic Acid (IUPAC: 4-hydroxyquinooline-2-carboxyclic acid) and 2-methylindole were purchased from Sigma-Aldrich and were ≥98% pure. The standard for 5-Methoxy-1H-indole-3-carbaldehyde was purchased from Matrix Scientific and was ≥98% pure. The target compounds in samples were detected and quantified on a triple quadrupole linear ion trap mass spectrometer (3200 QTRAP, AB SCIEX, Foster City, CA) coupled to a binary pump HPLC (Prominence LC-20, Shimazu, Concord, Ontario, Canada). Peak identification and integration were performed using Analyst software (version 5, Agilent, Foster City, CA). Samples were maintained at 4°C on an autosampler prior to injection. Chromatographic separation was achieved on a C-18 column (Synergi Fusion 4 μm 80 Å 150 mm × 2 mm, Phenomenex, Torrance, CA) using a solvent gradient method [23]. Solvent A was formic acid solution in water (0.1%). Solvent B was acetonitrile with formic acid (0.1%). Injection volume was 10 μL. For urine, concentrations of compound within each urine sample were normalized to urine creatinine concentration and expressed as ratios. Serum was normalized by volume.

Statistical analysis

A subset of 6,900 sequences per sample was randomly selected to account for unequal sequencing depth across samples. Differences in microbial communities between healthy dogs, dogs with NHD, and dogs with AHD were analyzed using the phylogeny-based unweighted UniFrac distance metric, PCoA plots, and rarefaction curves generated by QIIME [18]. Rarefaction curves and PCoA plots show alpha (i.e., Chao 1, Shannon Index, and Observed Species) and beta (i.e., microbial community distance matrix) diversity, respectively. ANOSIM (Analysis of Similarity) within the software package PRIMER 6 (PRIMER-E Ltd., Luton, UK) was used to determine significant differences in microbial communities between healthy dogs and diseased dogs. To visualize the relative abundance of bacterial families for individual fecal samples, heat maps were generated in NCSS 2007 (NCSS, Kaysville, Utah).

All datasets were tested for normality using the Shapiro-Wilk test (JMP 10, SAS software Inc.). Because most datasets did not meet the assumptions of normal distribution, comparisons between healthy and disease groups were determined using non-parametric Kruskal-Wallis tests (healthy dogs vs. dogs with NHD vs. dogs with AHD) or a Mann-Whitney U test (healthy dogs vs. dogs with acute diarrhea AD [dogs with NHD and dogs with AHD combined]). Taxa that were present in at least 70% of dogs (either healthy or diseased) were included in 454-pyrosequencing data analysis. The resulting p-values of the Kruskal-Wallis tests or Mann-Whitney U test were adjusted for multiple comparisons using the Benjamini & Hochberg’s False Discovery Rate (FDR), and an adjusted p<0.05 was considered statistically significant [24]. A Dunn’s post-test was used to determine which disease types were significantly different if applicable.

PICRUSt

The software PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) was used to make functional gene content predictions based on 16S rRNA gene data present in the Greengenes database [25]. PICRUSt is freely available online in the Galaxy workflow framework and can also be used through the QIIME open-source pipeline [26, 27].

LEfSe

Linear discriminant analysis effect size (LEfSe) was used to elucidate taxa and genes associated with healthy or diseased states. For bacterial groups, the LDA score threshold was set to > 3.5; for functional genes and their specific KEGG orthologs the LDA score threshold was set to > 2.5. LEfSe is freely available online in the Galaxy workflow framework [26, 27].

Short-chain fatty acids

Fecal SCFA concentrations were tested for normality and compared between each group using a Mann-Whitney U test. SCFAs that were found to be significantly different between groups were then correlated to differentially abundant bacterial groups (i.e., Spearman’s rank correlation; JMP 10, SAS software Inc.).

Serum and urine untargeted metabolomics

Waters raw files were converted to netCDF format, and XCMS [28] peak detection, retention time alignment, and feature grouping were performed on both the low and high collision energy channels (MS and MSe). The datasets were separated following alignment, normalized to total extracted ion signal, averaged by injection replicate, and the idMS/MS workflow described previously was applied to all significant features for generation of indiscriminant MS/MS (idMS/MS) spectra for library searching and compound identification [29]. For identification of spectra, recreated idMS/MS spectra were searched against the NIST ‘12 MS/MS database, Metlin, and Massbank. idMS/MS spectra which matched database MS/MS spectra were assigned a level 2 identification confidence [30]. For urine, concentrations of compound within each urine sample were normalized to urine creatinine concentration and expressed as ratios. Serum was normalized by volume. After data pre-processing, principle component analysis and univariate analysis was conducted using the web-based metabolomic data processing tool MetaboAnalyst 2.0 [31, 32].

Results

Sequencing analysis

The analysis pipeline yielded 297,315 quality sequences for the 19 samples analyzed (mean ± standard deviation [of all samples] = 9,013 ± 1,203). Fig 1 illustrates the rarefaction curve for observed species (a count of all unique operational taxonomic units (OTUs)). All alpha diversity measures are summarized in Table 1. Observed species were significantly decreased in dogs with acute diarrhea (p = 0.0218). The Shannon Index was significantly decreased in dogs with AD and NHD compared to healthy dogs (p = 0.0033 and p < 0.0500, respectively). The Chao1 predictive diversity measure was also significantly decreased in dogs with acute diarrhea (p = 0.0176).

[Figure omitted. See PDF.]

Fig 1. Rarefaction analysis of 16S rRNA gene sequences obtained from canine fecal samples.

Lines represent the mean and error bars represent standard deviations. The analysis was performed on a randomly selected subset of 6,900 sequences per sample. NHD = acute non-hemorrhagic diarrhea; AHD = acute hemorrhagic diarrhea.

https://doi.org/10.1371/journal.pone.0127259.g001

Microbial communities

PCoA plots (Fig 2) based on unweighted Unifrac distances showed significant differences between healthy dogs and dogs with AD (ANOSIM; p = 0.0040). Furthermore dogs with NHD and dogs with AHD differed significantly from healthy dogs (ANOSIM; p = 0.0020 for both). There was no difference in microbial communities, however, between dogs with hemorrhagic and non-hemorrhagic diarrhea. Based on LDA effect size (LEfSe), Clostridium spp. was significantly associated with AD, while Prevotella spp., Blautia spp., Faecalibacterium spp., Eubacterium spp., and unclassified genera within the following families: Ruminococcaceae, Lachnospiraceae, Clostridia, Ruminococcaceae, and Coprobacillaceae were significantly associated with healthy dogs (Fig 3). A heatmap in Fig 4 was employed to illustrate these changes.

[Figure omitted. See PDF.]

Fig 2. Principal Coordinate Analysis (PCoA) of unweighted UniFrac distances of 16S rRNA genes.

Blue squares = healthy dogs, green circles = acute non-hemorrhagic diarrhea (NHD), and red triangles = acute hemorrhagic diarrhea (AHD); ANOSIM for healthy dogs vs. dogs with AD (NHD and AHD combined), p = 0.0040; and ANOSIM for NHD or AHD vs. healthy dogs, p = 0.0020 for both.

https://doi.org/10.1371/journal.pone.0127259.g002

[Figure omitted. See PDF.]

Fig 3. Differentially abundant bacterial groups.

Groups differentially abundant between healthy dogs and dogs with acute diarrhea. Red bars represent bacterial groups associated with dogs with acute diarrhea, while green bars represent bacterial groups associated with healthy dogs.

https://doi.org/10.1371/journal.pone.0127259.g003

[Figure omitted. See PDF.]

Fig 4. Heatmap illustrating the relative abundance of predominant bacterial genera in fecal samples.

Healthy = healthy dogs; NHD = acute non-hemorrhagic diarrhea; AHD = acute hemorrhagic diarrhea. “Unclass.” denotes an unclassified genus within the respective taxa.

https://doi.org/10.1371/journal.pone.0127259.g004

Table 2 and Fig 5A summarize differences in bacterial groups between disease groups. Sequences belonging to the phylum Bacteroidetes were significantly decreased in dogs with AD compared to healthy dogs (p = 0.0280). Sequences belonging to the genus Faecalibacterium and an unclassified genus within Ruminococcaceae were both significantly decreased in dogs with AD compared to healthy dogs (p = 0.0319 and 0.0368, respectively). Sequences belonging to the genus Clostridium were significantly increased in dogs with AD compared to healthy dogs (p = 0.0476). No differences were identified between dogs with hemorrhagic and non-hemorrhagic acute diarrhea when evaluating bacterial groups.

[Figure omitted. See PDF.]

Table 2. Percentages of the most abundant bacterial groups.

https://doi.org/10.1371/journal.pone.0127259.t002

[Figure omitted. See PDF.]

Fig 5. Groups of interest based on A) Sequencing (% of sequences) and B) quantitative PCR results. The qPCR data was expressed as log amount of DNA for each particular bacterial group per 10 ng of isolated total DNA.

Acute Diarrhea = both groups combined (NHD and AHD). Bars represent the median value for each group. P-values adjusted based on the Benjamini and Hochberg false discovery rate.

https://doi.org/10.1371/journal.pone.0127259.g005

Quantitative PCR

The abundance of Clostridium perfringens was significantly increased in dogs with AD (log DNA median [range]: 7 [5.8–7.5]) compared to healthy dogs (log DNA median [range]: 4.6 [3–6.1]; p = 0.0088). Clostridium perfringens was also significantly increased in the subgroup dogs with AHD (log DNA median [range]: 7.1 [6.9–7.4]) compared to healthy dogs (p<0.0500). No significant difference in Bifidobacterium, Lactobacillus, or E. coli was identified between dogs with acute diarrhea and healthy dogs (Fig 5B).

Functional genes

Univariate statistics revealed no significant differences in the percentage of KEGG orthologs belonging to functional gene families at all levels (e.g., 1, 2, and 3) among all groups of dogs after correcting for multiple comparisons (S2 Table). However, PICRUSt provided a snapshot of the distribution of genes across functional categories. At level 1, approximately 50% of genes belonged to metabolism, 19% belonged to genetic information processing, and 15% belonged to environmental information processing. Next, functional gene categories were analyzed using LEfSe. The cladogram in Fig 6A shows functional genes associated with healthy dogs and dogs with AD at all levels in the KEGG database hierarchy (LDA score threshold > 2.5). The following pathways were associated with acute diarrhea: benzoate degradation within xenobiotics biodegradation/metabolism, lipid metabolism, two component system within signal transduction, flagellar assembly/bacterial motility proteins/bacterial chemotaxis within cell motility, and metabolism. The following pathways were associated with healthy dogs: methane metabolism/oxidative phosphorylation within energy metabolism and chaperones/folding catalysts within folding sorting/degradation.

[Figure omitted. See PDF.]

Fig 6. A) Differentially abundant gene families identified within healthy dogs and dogs with acute diarrhea. B) Taxa contributing to the MAC protein gene and Tranposase enzyme gene (K03406 and K07483, respectively) (Table 3) in all dogs.

H.__ = Healthy, A.__ = acute non-hemorrhagic diarrhea, R.__ acute hemorrhagic diarrhea.

https://doi.org/10.1371/journal.pone.0127259.g006

In a separate analysis, KEEG orthologs (KOs) were analyzed by LEfSe without being categorized by functional hierarchy (this strategy allows for more finite discovery of specific genes associated with healthy dogs or dogs with AD) (Table 3). KOs found to be associated with healthy dogs were K07133, K07720, K02027, K02026, K0205, K10439, and K07718. These genes are responsible for the multiple sugar transport system permease protein, multiple sugar transport system substrate-binding protein, ribose transport system substrate-binding protein, two-component system/sensor histidine kinase YesM, ATPase enzyme, and two-component system/response regulator YesN, respectively. KOs found to be associated with dogs with acute diarrhea were K03406 and K07483. These genes are responsible for the methyl-accepting chemotaxis protein and transposase enzyme, respectively. PICRUSt allows the user to determine which taxa contribute to functional gene predictions (Fig 6B). To gain an overall understanding of taxa contributing to identification of genes associated with acute diarrhea, only taxa contributing 500 or more gene counts were included in the bar charts. K03406 and K07483 were both identified as genes that were significantly associated with acute diarrhea. The charts in Fig 6B show that the family Clostridiaceae contributed mainly to genes coding for the methyl accepting chemotaxis protein as well as the transposase enzyme (K03406 and K07483, respectively). In addition, more gene counts appeared to be contributed by the family Clostridiaceae in dogs with acute diarrhea compared to healthy dogs.

[Figure omitted. See PDF.]

Table 3. KEGG orthologs and functions that are highly abundant.

https://doi.org/10.1371/journal.pone.0127259.t003

Short-chain fatty acids

Butyric acid, acetic acid, and propionic acid were each expressed as a percentage of the total amount of fecal SCFA concentrations (Fig 7). The proportion of propionic acid was significantly decreased in dogs with AD (median [range]:12% [0–25%]) compared to healthy dogs (median [range]: 30% [20–44%]; p = 0.0033). In contrast, the proportion of butyric acid was significantly increased in dogs with AD (median [range]: 12% [8–26%]) compared to healthy dogs (median [range]: 6% [4–8%]; p = 0.0048). There were no significant differences in total fecal SCFA or BCFA concentrations (μmol/g of dry feces) observed between groups of dogs (S3 Table).

[Figure omitted. See PDF.]

Fig 7. Proportion of total fecal SCFA concentrations.

Circles = healthy dogs; Squares = acute diarrhea both groups combined (NHD and AHD). Bars represent the median value for each group. P-values adjusted based on the Benjamini and Hochberg false discovery rate.

https://doi.org/10.1371/journal.pone.0127259.g007

Correlations between SCFA and bacterial groups

A positive correlation was observed between propionic acid and an unclassified genus within Ruminococcaceae as well as Faecalibacterium (ρ = 0.8377 [p = 0.0042] and ρ = 0.6725 [p = 0.0332], respectively) (Fig 8). Butyric acid was found to have a negative correlation with an unclassified genus within Ruminococcaceae (ρ = -0.8265 [p = 0.0027]. Furthermore, propionic acid was negatively correlated with Lactobacillus (ρ = -0.8333 [p = 0.0424]).

[Figure omitted. See PDF.]

Fig 8. Correlations between SCFAs and bacterial groups.

A) Correlations between butyric acid and select bacterial groups. B) Correlations between propionic acid and select bacterial groups.

https://doi.org/10.1371/journal.pone.0127259.g008

Serum and urine metabolites

PCA plots showed no clustering in serum or urine metabolite profiles (Figs 9 and 10, respectively) between healthy dogs and dogs with acute diarrhea. A total of 82 unique compounds were identified in serum, while 362 unique compounds were identified in urine. Of these compounds, none were significantly changed after adjusting for multiple comparisons in either of the analyzed sample types. However, in serum, unadjusted p-values showed that the concentration of kynurenic acid in serum was significantly decreased in dogs with acute diarrhea compared to healthy dogs (p = 0.0048). In addition to this comparison, the ratio between kynurenic acid (K) and tryptophan (T) concentration in the serum of healthy dogs and dogs with acute diarrhea was compared. The ratio of K/T was significantly decreased in dogs with acute diarrhea (p = 0.0036).

[Figure omitted. See PDF.]

Fig 9. Analysis of serum metabolites.

A) Principle component analysis plot where red dots represent dogs with acute diarrhea while green dots represent healthy dogs. Ellipses show the 95% confidence distribution for each group. B) Relative concentration of kynurenic acid, tryptophan, and their ratio to one another in serum (p = 0.0048, 0.8282, and 0.0036; respectively).

https://doi.org/10.1371/journal.pone.0127259.g009

[Figure omitted. See PDF.]

Fig 10. Analysis of urine metabolites.

A) Principle component analysis plot where red dots represent dogs with acute diarrhea while green dots represent healthy dogs. Ellipses show the 95% confidence distribution for each group. B) Relative concentration of 2-methylindole and 5-methoxy-1H-indole-3-carbaldehyde in serum (p = 0.0185 and 0.0330, respectively).

https://doi.org/10.1371/journal.pone.0127259.g010

In urine, unadjusted p-values showed that the concentrations of 2-methyl-1H-indole and 5-Methoxy-1H-indole-3-carbaldehyde were significantly decreased in dogs with acute diarrhea compared to healthy dogs (p = 0.0185 and 0.0330, respectively).

A targeted approach was used to confirm the identity and validity of compounds found to be significantly altered from the untargeted metabolomics approach. The targeted approach for serum kynurenic acid yielded a strong correlation to the untargeted approach (ρ = 0.7213 and p = 0.0001) (S2 Fig). The targeted approach for urine 2-methyl-1H-indole and 5-Methoxy-1H-indole-3-carbaldehyde yielded a moderate to strong significant correlation to the untargeted approach (ρ = 0.6569 and p = 0.0042, ρ = 0.8162 and p<0.0001; respectively) (S3 Fig)

Discussion

In this study we compared the fecal microbiome as well as serum and urine metabolomic profiles between healthy dogs and dogs with acute diarrhea. Our results revealed significant changes in microbial communities in dogs with AD. Consistent with our previous findings in dogs with acute diarrhea and similar to data observed in humans and mice, alpha diversity measures were significantly decreased [14, 33, 34]. Principal coordinates analysis (PCoA) plots revealed significant differences in microbiome composition between healthy dogs and dogs with AD. We identified specific bacterial groups that were altered in disease. A significant increase in Clostridium perfringens in dogs with AD was observed as reported previously [5, 14]. In contrast, bacterial groups decreased in AD were Bacteroidetes, Faecalibacterium, and an unclassified genus within Ruminococcaceae. Some of these bacterial groups are believed to be important producers of various metabolites including SCFAs [35], and consequently, we observed significant differences in SCFAs patterns between healthy and diseased dogs, with propionic acid being significantly decreased in dogs with AD. However, unexpectedly butyric acid was significantly increased in dogs with AD. Studies in humans have demonstrated the importance of SCFAs such as butyrate and propionate. For example, propionic acid has recently garnered attention for its protective effects against carcinogenesis and colorectal cancer in humans [36]. Butyrate protects against colitis by inducing apoptosis in cells with DNA damage, and also increases the expression of tight junction proteins thereby reinforcing colonic defense barriers [37–39]. Based on a decrease in Faecalibacterium, which has been identified as a butyrate producing bacteria in humans [40, 41], a decrease in fecal butyrate concentrations was expected in this study. This, however, was not observed. It is possible that dogs with acute diarrhea experienced decreased utilization of butyric acid by the colonic epithelium and consequently butyric acid was excreted and found in higher concentrations in the feces. Faecalibacterium and an unclassified genus within Ruminococcaceae shared a positive correlation with propionic acid. Alterations in fecal propionic acid could possibly be due to decreased production and/or increased absorption into the gut epithelium during stages of acute diarrhea. Therefore, it appears that the dysbiosis in AD had direct impact on concentration of SCFA and this warrants further research into potential therapeutic applications. Despite these findings, it is important to note that measuring fecal SCFAs may not directly represent colonic conditions and further studies are needed to investigate SCFA concentrations in vivo.

PICRUSt was used to infer putative metagenomes from16S rRNA gene profiles [16]. In general, there was a homogenous distribution of genes across broad functional categories, which has been previously noted in literature [42]. For instance, the abundance of genes belonging to metabolism (50%), Genetic Information Processing (20%), Environmental Information Processing (15%), and Cellular Processes (1–2%) were similar across healthy and diseased dogs. This similarity between both groups could be due to convergence in gene content among bacteria that inhabit the GIT. LEfSe was used to identify differentially abundant functional genes (i.e., KEGG orthologs) that had not yet been grouped into families and some notable differences were identified. The methyl-accepting chemotaxis (MAC) protein gene (K03406) was significantly increased in AD. The MAC protein is an integral component of ribose and galactose chemoreceptors in Escherichia coli and has been suggested to play a role in the pathogenesis of E. coli [43]. It may be hypothesized that an increase in the number of genes encoding the MAC protein may contribute to the survival of pathogenic microorganisms that ultimately cause acute diarrhea [44, 45]. A gene encoding transposases (K07483) was also found to be significantly more abundant in AD. Transposases can induce various types of genome rearrangements and are a major cause of mutations [46]. In part, they can be detrimental by inactivating housekeeping genes or impairing the chromosomes integrity [47]. Adversely, transposases are essential to horizontal gene transfer (HGT), which in turn is largely responsible for the diversity, strong selection, and co-evolution of microbiota in humans [48]. Transposases were recently identified as enzymes largely effected by HGT in a large-scale comparative metagenomic analysis of human fecal samples. Evaluating transposase and HGT activity may help determine functional factors that that contribute to GI disease. Alterations in functional genes greatly increase our understanding of the disease processes during acute diarrhea. However, given the limitations of PICRUSt as solely a predictor of metagenomic function, a true metagenomics approach is warranted and may yield more in-depth resolution.

Serum and urine metabolites were evaluated using an untargeted metabolomics approach. While we observed major changes in microbial communities, only a few serum or urine metabolites were significantly altered before adjusting for multiple comparisons. However, after adjustment this significance was lost. A possibility is that acute onset of diarrhea leads to only minor systemic changes detectable in serum or urine due to the short time of disease. This may be in contrast to chronic GI inflammation, where alterations of serum metabolites associated with oxidative stress have been observed in dogs with IBD [20]. Another explanation is possibly due to a limited sample size in conjunction with a large set of comparisons which were used for statistical adjustment for multiple comparisons. This adjustment is inherently conservative given the breadth of metabolites discovered by an untargeted metabolomics approach. Therefore, targeted assays with standards were developed to validate the serum and urine metabolites that were significantly altered based on unadjusted p-values. The results of the targeted assays correlated significantly with the untargeted metabolomics data verifying the results.

Although it has long been recognized that the enteric community of bacteria broadly impacts host health, the biochemical products they produce are only recently becoming better understood. In this study, derivatives and metabolic end products of tryptophan (i.e., 2-methylindole, 5-methoxy-1H-indole-3-carbaldehyde) were significantly decreased in the urine of dogs with acute diarrhea. These alterations may reflect a difference in uptake or utilization. While the relationship between indole and indole-like compounds is unclear in the context of dogs with acute diarrhea, it is interesting to note that indole has been shown to increase expression of genes involved in strengthening of the mucosal barrier and mucin production [49]. Also, the same study showed that indole decreased TNF-α mediated activation of NF-κB, as well as increased expression of the anti-inflammatory cytokine IL-10 [49]. Kynurenic acid, a catabolic derivative of tryptophan, was significantly decreased in the serum of dogs with AD [49]. The importance of kynurenic acid in GI disease is not yet fully understood but warrants further attention. The intestinal mucosa can monitor microbial ligands via pattern recognition receptors and microbial metabolites via G-protein coupled receptors (GPRs), and GPR35 is thought to be activated by kynurenic acid [50–52]. Recently, a study reported decreased mucosal but increased systemic levels of kynurenic acid in patients with irritable bowel syndrome (IBS) [53]. Indoleamine 2,3-dioxygenase (IDO) and tryptophan 2,3,-dioxygenase (TDO) are two enzymes implicated in catabolism of tryptophan and production of kynurenine. These enzymes play a central role in the physiological regulation of tryptophan alterations in the human body and catalyze the first and rate limiting step of tryptophan degradation along the kynurenine pathway [54]. Studies in humans with other forms of GI disease such as Crohn’s disease (CD) have suggested that aberrant handling of luminal bacteria by the innate immune system results in an activated inflammatory cascade and secretion of cytokines (i.e., TNF-α and IFN-γ) [55, 56]. These cytokines are thought to be responsible for triggering enzymatic activation of IDO and subsequent catabolism of tryptophan to form kynurenine. It may be hypothesized that the enzymatic activation of IDO could negatively impact the immune system; however this conjecture warrants further research.

The kynurenine/tryptophan ratio is used as a surrogate marker for IDO1 (i.e., the gene responsible for IDO) activity in patients with CD [57]. In a well-characterized cohort of CD patients and controls, one study found that serum tryptophan was decreased while the kynurenine/tryptophan ratio was increased in active CD [58]. While kynurenine was not measured in this study, the ratio of serum kynurenic acid to tryptophan was investigated and found to be significantly decreased in dogs with acute diarrhea. It’s not surprising that these results do not align as the kynurenine/tryptophan ratio is often used to measure up- or down-regulation of the kynurenine pathway and unfortunately in this case we were unable to readily measure kynurenine, only one of its derivatives. It may be useful to target and monitor changes in kynurenine and tryptophan in a longitudinal study which may provide a clearer picture of the role these metabolites play in disease processes.

This study had limitations. Only a small number of animals were enrolled in the disease and healthy group. In this respect, we acknowledge that results should be viewed as descriptive and are useful in justifying larger scale studies to further confirm our findings. Evaluating pet dogs living in various home environments is inherently difficult due to differences in size, weight, obesity status and sex (intact vs. not-intact). In this study, the potential influence of gender on microbial communities (S1 Fig) was investigated using the unweighted Unifrac distance matrix on all dogs and there were no significant differences between males or females (ANOSIM: p = 0.6400) regardless of sexual status (ANOSIM: p = 0.7900). Studies in humans and mice have indicated that obesity plays a role in microbiome composition [59, 60]. The current data in dogs regarding the intestinal microbiome in dogs is conflicting, and a previous study has not shown major differences between lean and obese pet dogs [61]. The body condition score (BCS, on a scale of 1–9 with 4–5 being ideal) used to differentiate between lean and obese dogs (S1 Table) was fairly heterogeneous amongst all dogs and clustering analysis of unweighted Unifrac distances revealed no significant clustering by BCS. Additionally, in some instances, only a limited amount of sample was available; therefore some samples could not be included for all analyses. Therefore, it is possible that some alterations may have been overlooked. All dogs were housed in different home environments before presentation, which may have a minor effect on GI microbiota. Diet has also been shown to alter the intestinal microbiota [62]. However, in this study, dogs in both groups were on a maintenance diet, so this fact likely had a minimal effect on microbial groups.

Conclusions

Results of this study revealed a bacterial dysbiosis in fecal samples of dogs with acute diarrhea. Microbial diversity was significantly decreased and microbial communities differed significantly in dogs with AD. Predicted gene families were conserved across both groups of dogs. Furthermore, propionic acid was significantly correlated with a loss in SCFA producing bacteria (i.e., Faecalibacterium). The fecal dysbiosis was associated with significant changes in profiles of fecal SCFA and serum and urine metabolites. This suggests that acute episodes of diarrhea, even when uncomplicated, have an impact on the overall metabolic profile of the host.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. Principal Coordinate Analysis (PCoA) of unweighted UniFrac distances of 16S rRNA genes representing clustering of samples categorized by sex (A) and sexual status (B).

https://doi.org/10.1371/journal.pone.0127259.s001

(PDF)

S2 Fig. Comparison between serum untargeted and targeted metabolomics approach for kynurenic acid.

A) Untargeted kynurenic acid. B) Targeted kynurenic acid. C) Correlation between both approaches.

https://doi.org/10.1371/journal.pone.0127259.s002

(TIF)

S3 Fig. Comparison between urine untargeted and targeted metabolomics approach for 2-methyl-1H-indole and 5-methoxy-1H-indole-3-carbaldehyde.

A-C) Represent untargeted, targeted, and correlation results for 2-methyl-1H-indole, respectively. D-F) Represent untargeted, targeted, and correlation results for 5-methoxy-1H-indole-3-carbaldehyde, respectively.

https://doi.org/10.1371/journal.pone.0127259.s003

(TIF)

S1 Table. Signalment of dogs enrolled into this study.

https://doi.org/10.1371/journal.pone.0127259.s004

(PDF)

S2 Table. Percentages of KEGG orthologs that belong to gene families at levels 1, 2, and 3.

https://doi.org/10.1371/journal.pone.0127259.s005

(PDF)

S3 Table. Total concentrations of SCFAs and BCFAs.

https://doi.org/10.1371/journal.pone.0127259.s006

(PDF)

S4 Table. Oligonucleotide primers/probe used for this study.

https://doi.org/10.1371/journal.pone.0127259.s007

(PDF)

Author Contributions

Conceived and designed the experiments: BCG JWB JSS. Performed the experiments: BCG LR JV. Analyzed the data: BCG JSS. Contributed reagents/materials/analysis tools: JWB LR JMS JV JSS AJ CK. Wrote the paper: BCG LR JV JSS.

Citation: Guard BC, Barr JW, Reddivari L, Klemashevich C, Jayaraman A, Steiner JM, et al. (2015) Characterization of Microbial Dysbiosis and Metabolomic Changes in Dogs with Acute Diarrhea. PLoS ONE 10(5): e0127259. https://doi.org/10.1371/journal.pone.0127259

References

1. Frank DN, Zhu W, Sartor RB, Li E. Investigating the biological and clinical significance of human dysbioses. Trends in microbiology. 2011;19(9):427–34. pmid:21775143; PubMed Central PMCID: PMC3164499.

2. Handl S, Dowd SE, Garcia-Mazcorro JF, Steiner JM, Suchodolski JS. Massive parallel 16S rRNA gene pyrosequencing reveals highly diverse fecal bacterial and fungal communities in healthy dogs and cats. FEMS microbiology ecology. 2011;76(2):301–10. pmid:21261668.

3. Suchodolski JS, Dowd SE, Westermarck E, Steiner JM, Wolcott RD, Spillmann T, et al. The effect of the macrolide antibiotic tylosin on microbial diversity in the canine small intestine as demonstrated by massive parallel 16S rRNA gene sequencing. BMC microbiology. 2009;9:210. pmid:19799792; PubMed Central PMCID: PMC2759960.

4. Zoetendal EG, Collier CT, Koike S, Mackie RI, Gaskins HR. Molecular ecological analysis of the gastrointestinal microbiota: a review. The Journal of nutrition. 2004;134(2):465–72. pmid:14747690.

5. Bell JA, Kopper JJ, Turnbull JA, Barbu NI, Murphy AJ, Mansfield LS. Ecological characterization of the colonic microbiota of normal and diarrheic dogs. Interdisciplinary perspectives on infectious diseases. 2008;2008:149694. pmid:19282974; PubMed Central PMCID: PMC2648299.

6. Inness VL, McCartney AL, Khoo C, Gross KL, Gibson GR. Molecular characterisation of the gut microflora of healthy and inflammatory bowel disease cats using fluorescence in situ hybridisation with special reference to Desulfovibrio spp. Journal of animal physiology and animal nutrition. 2007;91(1–2):48–53. pmid:17217390.

7. Janeczko S, Atwater D, Bogel E, Greiter-Wilke A, Gerold A, Baumgart M, et al. The relationship of mucosal bacteria to duodenal histopathology, cytokine mRNA, and clinical disease activity in cats with inflammatory bowel disease. Veterinary microbiology. 2008;128(1–2):178–93. pmid:18054447.

8. Suchodolski JS. Companion animals symposium: microbes and gastrointestinal health of dogs and cats. Journal of animal science. 2011;89(5):1520–30. pmid:21075970.

9. Suchodolski JS, Xenoulis PG, Paddock CG, Steiner JM, Jergens AE. Molecular analysis of the bacterial microbiota in duodenal biopsies from dogs with idiopathic inflammatory bowel disease. Veterinary microbiology. 2010;142(3–4):394–400. pmid:19959301.

10. Xenoulis PG, Palculict B, Allenspach K, Steiner JM, Van House AM, Suchodolski JS. Molecular-phylogenetic characterization of microbial communities imbalances in the small intestine of dogs with inflammatory bowel disease. FEMS microbiology ecology. 2008;66(3):579–89. pmid:18647355.

11. Nicholson JK, Lindon JC. Systems biology: Metabonomics. Nature. 2008;455(7216):1054–6. pmid:18948945.

12. Wikoff WR, Anfora AT, Liu J, Schultz PG, Lesley SA, Peters EC, et al. Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. Proceedings of the National Academy of Sciences of the United States of America. 2009;106(10):3698–703. pmid:19234110; PubMed Central PMCID: PMC2656143.

13. Rossi G, Pengo G, Caldin M, Palumbo Piccionello A, Steiner JM, Cohen ND, et al. Comparison of microbiological, histological, and immunomodulatory parameters in response to treatment with either combination therapy with prednisone and metronidazole or probiotic VSL#3 strains in dogs with idiopathic inflammatory bowel disease. PloS one. 2014;9(4):e94699. pmid:24722235; PubMed Central PMCID: PMC3983225.

14. Suchodolski JS, Markel ME, Garcia-Mazcorro JF, Unterer S, Heilmann RM, Dowd SE, et al. The fecal microbiome in dogs with acute diarrhea and idiopathic inflammatory bowel disease. PloS one. 2012;7(12):e51907. pmid:23300577; PubMed Central PMCID: PMC3530590.

15. Foster ML, Dowd SE, Stephenson C, Steiner JM, Suchodolski JS. Characterization of the fungal microbiome (mycobiome) in fecal samples from dogs. Veterinary medicine international. 2013;2013:658373. pmid:23738233; PubMed Central PMCID: PMC3655676.

16. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nature biotechnology. 2013;31(9):814–21. pmid:23975157; PubMed Central PMCID: PMC3819121.

17. Dowd SE, Sun Y, Secor PR, Rhoads DD, Wolcott BM, James GA, et al. Survey of bacterial diversity in chronic wounds using pyrosequencing, DGGE, and full ribosome shotgun sequencing. BMC microbiology. 2008;8:43. pmid:18325110; PubMed Central PMCID: PMC2289825.

18. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nature methods. 2010;7(5):335–6. pmid:20383131; PubMed Central PMCID: PMC3156573.

19. Garcia-Mazcorro JF, Suchodolski JS, Jones KR, Clark-Price SC, Dowd SE, Minamoto Y, et al. Effect of the proton pump inhibitor omeprazole on the gastrointestinal bacterial microbiota of healthy dogs. FEMS microbiology ecology. 2012;80(3):624–36. pmid:22324305.

20. Minamoto Y, Otoni CC, Steelman SM, Buyukleblebici O, Steiner JM, Jergens AE, et al. Alteration of the fecal microbiota and serum metabolite profiles in dogs with idiopathic inflammatory bowel disease. Gut microbes. 2014:1–15. Epub 2014/12/23. pmid:25531678.

21. Minamoto Y, Dhanani N, Markel ME, Steiner JM, Suchodolski JS. Prevalence of Clostridium perfringens, Clostridium perfringens enterotoxin and dysbiosis in fecal samples of dogs with diarrhea. Veterinary microbiology. 2014;174(3–4):463–73. Epub 2014/12/03. pmid:25458422.

22. Moreau NM, Goupry SM, Antignac JP, Monteau FJ, Le Bizec BJ, Champ MM, et al. Simultaneous measurement of plasma concentrations and 13C-enrichment of short-chain fatty acids, lactic acid and ketone bodies by gas chromatography coupled to mass spectrometry. Journal of chromatography B, Analytical technologies in the biomedical and life sciences. 2003;784(2):395–403. pmid:12505787.

23. Zhu W, Stevens A, Dettmer K, Gottfried E, Hoves S, Kreutz M, et al. Quantitative profiling of tryptophan metabolites in serum, urine, and cell culture supernatants by liquid chromatography—tandem mass spectrometry. Anal Bioanal Chem. 2011;401(10):3249–61. pmid:21983980

24. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological). 1995:289–300.

25. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nature biotechnology. 2013. pmid:23975157.

26. Blankenberg D, Von Kuster G, Coraor N, Ananda G, Lazarus R, Mangan M, et al. Galaxy: a web-based genome analysis tool for experimentalists. Current protocols in molecular biology / edited by Frederick M Ausubel [et al]. 2010;Chapter 19:Unit 19 0 1-21 pmid:20069535.

27. Goecks J, Nekrutenko A, Taylor J, Galaxy T. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome biology. 2010;11(8):R86. pmid:20738864; PubMed Central PMCID: PMC2945788.

28. Smith CA, Want EJ, O'Maille G, Abagyan R, Siuzdak G. XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Analytical chemistry. 2006;78(3):779–87. pmid:16448051.

29. Broeckling CD, Heuberger AL, Prince JA, Ingelsson E, Prenni JE. Assigning precursor-product ion relationships in indiscriminant MS/MS data from non-targeted metabolite profiling studies. Metabolomics. 2013;9(1):33–43. WOS:000313736700005.

30. Sumner LW, Amberg A, Barrett D, Beale MH, Beger R, Daykin CA, et al. Proposed minimum reporting standards for chemical analysis. Metabolomics. 2007;3(3):211–21. WOS:000249945500006. pmid:24039616

31. Xia J, Mandal R, Sinelnikov IV, Broadhurst D, Wishart DS. MetaboAnalyst 2.0—a comprehensive server for metabolomic data analysis. Nucleic acids research. 2012;40(Web Server issue):W127–33. pmid:22553367; PubMed Central PMCID: PMC3394314.

32. Xia JG, Broadhurst DI, Wilson M, Wishart DS. Translational biomarker discovery in clinical metabolomics: an introductory tutorial. Metabolomics. 2013;9(2):280–99. WOS:000316913200003. pmid:23543913

33. Frank DN, St Amand AL, Feldman RA, Boedeker EC, Harpaz N, Pace NR. Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proceedings of the National Academy of Sciences of the United States of America. 2007;104(34):13780–5. pmid:17699621; PubMed Central PMCID: PMC1959459.

34. Lynch SV, Goldfarb KC, Wild YK, Kong W, De Lisle RC, Brodie EL. Cystic fibrosis transmembrane conductance regulator knockout mice exhibit aberrant gastrointestinal microbiota. Gut microbes. 2013;4(1):41–7. pmid:23060053; PubMed Central PMCID: PMC3555885.

35. Topping DL, Clifton PM. Short-chain fatty acids and human colonic function: roles of resistant starch and nonstarch polysaccharides. Physiological reviews. 2001;81(3):1031–64. pmid:11427691.

36. Hinnebusch BF, Meng S, Wu JT, Archer SY, Hodin RA. The effects of short-chain fatty acids on human colon cancer cell phenotype are associated with histone hyperacetylation. The Journal of nutrition. 2002;132(5):1012–7. pmid:11983830.

37. Cook SI, Sellin JH. Review article: short chain fatty acids in health and disease. Alimentary pharmacology & therapeutics. 1998;12(6):499–507. pmid:9678808.

38. Wong JM, de Souza R, Kendall CW, Emam A, Jenkins DJ. Colonic health: fermentation and short chain fatty acids. Journal of clinical gastroenterology. 2006;40(3):235–43. pmid:16633129.

39. Antharam VC, Li EC, Ishmael A, Sharma A, Mai V, Rand KH, et al. Intestinal Dysbiosis and Depletion of Butyrogenic Bacteria in Clostridium difficile Infection and Nosocomial Diarrhea. Journal of clinical microbiology. 2013;51(9):2884–92. pmid:23804381; PubMed Central PMCID: PMC3754663.

40. Pryde SE, Duncan SH, Hold GL, Stewart CS, Flint HJ. The microbiology of butyrate formation in the human colon. FEMS microbiology letters. 2002;217(2):133–9. pmid:12480096.

41. Duncan SH, Hold GL, Harmsen HJ, Stewart CS, Flint HJ. Growth requirements and fermentation products of Fusobacterium prausnitzii, and a proposal to reclassify it as Faecalibacterium prausnitzii gen. nov., comb. nov. International journal of systematic and evolutionary microbiology. 2002;52(Pt 6):2141–6. pmid:12508881.

42. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480–4. pmid:19043404; PubMed Central PMCID: PMC2677729.

43. Kondoh H, Ball CB, Adler J. Identification of a Methyl-Accepting Chemotaxis Protein for the Ribose and Galactose Chemoreceptors of Escherichia-Coli. Proceedings of the National Academy of Sciences of the United States of America. 1979;76(1):260–4. WOS:A1979GF88500057. pmid:370826

44. Schott T, Kondadi PK, Hanninen ML, Rossi M. Comparative genomics of Helicobacter pylori and the human-derived Helicobacter bizzozeronii CIII-1 strain reveal the molecular basis of the zoonotic nature of non-pylori gastric Helicobacter infections in humans. BMC genomics. 2011;12:534. pmid:22039924; PubMed Central PMCID: PMC3234257.

45. Benov L, Fridovich I. Escherichia coli exhibits negative chemotaxis in gradients of hydrogen peroxide, hypochlorite, and N-chlorotaurine: products of the respiratory burst of phagocytic cells. Proceedings of the National Academy of Sciences of the United States of America. 1996;93(10):4999–5002. pmid:8643518; PubMed Central PMCID: PMC39395.

46. Berg DE, Howe MM, Ajioka JW. Mobile DNA: American Society for Microbiology Washington, DC; 1989.

47. Aziz RK, Breitbart M, Edwards RA. Transposases are the most abundant, most ubiquitous genes in nature. Nucleic acids research. 2010;38(13):4207–17. WOS:000280538600009. pmid:20215432

48. Ley RE, Peterson DA, Gordon JI. Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell. 2006;124(4):837–48. pmid:16497592.

49. Bansal T, Alaniz RC, Wood TK, Jayaraman A. The bacterial signal indole increases epithelial-cell tight-junction resistance and attenuates indicators of inflammation. Proceedings of the National Academy of Sciences of the United States of America. 2010;107(1):228–33. pmid:19966295; PubMed Central PMCID: PMC2806735.

50. Mogi C, Tobo M, Tomura H, Murata N, He XD, Sato K, et al. Involvement of proton-sensing TDAG8 in extracellular acidification-induced inhibition of proinflammatory cytokine production in peritoneal macrophages. Journal of immunology. 2009;182(5):3243–51. pmid:19234222.

51. Khor B, Gardet A, Xavier RJ. Genetics and pathogenesis of inflammatory bowel disease. Nature. 2011;474(7351):307–17. pmid:21677747; PubMed Central PMCID: PMC3204665.

52. Maslowski KM, Vieira AT, Ng A, Kranich J, Sierro F, Yu D, et al. Regulation of inflammatory responses by gut microbiota and chemoattractant receptor GPR43. Nature. 2009;461(7268):1282–6. pmid:19865172; PubMed Central PMCID: PMC3256734.

53. Keszthelyi D, Troost FJ, Jonkers DM, Kruimel JW, Leue C, Masclee AA. Decreased levels of kynurenic acid in the intestinal mucosa of IBS patients: relation to serotonin and psychological state. Journal of psychosomatic research. 2013;74(6):501–4. pmid:23731747.

54. Pilotte L, Larrieu P, Stroobant V, Colau D, Dolusic E, Frederick R, et al. Reversal of tumoral immune resistance by inhibition of tryptophan 2,3-dioxygenase. Proceedings of the National Academy of Sciences of the United States of America. 2012;109(7):2497–502. pmid:22308364; PubMed Central PMCID: PMC3289319.

55. Xavier RJ, Podolsky DK. Unravelling the pathogenesis of inflammatory bowel disease. Nature. 2007;448(7152):427–34. pmid:17653185.

56. Abraham C, Cho JH. Inflammatory bowel disease. The New England journal of medicine. 2009;361(21):2066–78. pmid:19923578; PubMed Central PMCID: PMC3491806.

57. Gupta NK, Thaker AI, Kanuri N, Riehl TE, Rowley CW, Stenson WF, et al. Serum analysis of tryptophan catabolism pathway: correlation with Crohn's disease activity. Inflammatory bowel diseases. 2012;18(7):1214–20. pmid:21823214; PubMed Central PMCID: PMC3235239.

58. Schrocksnadel K, Wirleitner B, Winkler C, Fuchs D. Monitoring tryptophan metabolism in chronic immune activation. Clinica chimica acta; international journal of clinical chemistry. 2006;364(1–2):82–90. pmid:16139256.

59. Schwiertz A, Taras D, Schafer K, Beijer S, Bos NA, Donus C, et al. Microbiota and SCFA in lean and overweight healthy subjects. Obesity. 2010;18(1):190–5. pmid:19498350.

60. Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006;444(7122):1027–131. http://www.nature.com/nature/journal/v444/n7122/suppinfo/nature05414_S1.html. pmid:17183312

61. Handl S, German AJ, Holden SL, Dowd SE, Steiner JM, Heilmann RM, et al. Faecal microbiota in lean and obese dogs2013 2013-05-01 00:00:00. 332–43 p. pmid:23301868

62. Turnbaugh PJ, Ridaura VK, Faith JJ, Rey FE, Knight R, Gordon JI. The effect of diet on the human gut microbiome: a metagenomic analysis in humanized gnotobiotic mice. Science translational medicine. 2009;1(6):6ra14. pmid:20368178; PubMed Central PMCID: PMC2894525.