http://orcid.org/0000-0003-2354-5903
Figures
Abstract
Antibiotic administration can be a cause of gastrointestinal disease in horses, creating a disruption in the normal population and function of bacteria found in the hindgut. The objective of this study was to describe the changes in the cecal and fecal microbiomes and metabolomes of clinically healthy horses before and after metronidazole administration. Metronidazole (15 mg/kg BID PO) was given to five horses with cecal cannulas. The study was suspended on Day 3 due to adverse gastrointestinal effects. Cecal and fecal samples were obtained before (Days minus52, m28, m14, and 0) and after (Days 7, 14, 28, and 52) metronidazole administration. DNA was extracted from the cecal and fecal samples, and 16S rRNA genes were sequenced. Richness and evenness indices were significantly decreased by metronidazole administration in both cecal and fecal samples, but the overall composition was only significantly changed in fecal samples on Day 3 (ANOSIM, p = 0.008). The most dominant phyla were Bacteroidetes and Firmicutes in all groups examined. In fecal samples, significant changes of the phyla Actinobacteria, Spirochaetes, Lentisphaerae, and Verrucomicrobia occurred on Day 3, which correlated with clinical signs of gastrointestinal disease. The metabolome was characterized by mass spectrometry-based methods and only named metabolites were included in the analysis. Fecal, but not cecal, metabolites were significantly affected by metronidazole. The fecal metabolites affected represent diverse metabolic pathways, such as the metabolism of amino acids, carbohydrates, lipids, nucleic acids and cofactors and vitamins. Metronidazole administration has potential to cause adverse effects in horses, alters the bacterial composition of the horse’s cecal and fecal content, and the metabolome of fecal samples.
Citation: Arnold CE, Isaiah A, Pilla R, Lidbury J, Coverdale JS, Callaway TR, et al. (2020) The cecal and fecal microbiomes and metabolomes of horses before and after metronidazole administration. PLoS ONE 15(5): e0232905. https://doi.org/10.1371/journal.pone.0232905
Editor: Ch Ratnasekhar, Queen's University Belfast, UNITED KINGDOM
Received: December 16, 2019; Accepted: April 23, 2020; Published: May 22, 2020
Copyright: © 2020 Arnold et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Data is available in: National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the accession number SRP119693 Metabolomicsworkbench.org under the submission ST001248
Funding: Funding for this project was provided by the Department of Large Animal Clinical Sciences. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript
Competing interests: The authors have declared that no competing interests exist.
Introduction
Amongst veterinary species, horses have a comparatively high degree of morbidity and mortality attributed to diseases of their gastrointestinal tract (GIT) [1]. An emerging body of evidence links GIT disease, such as colic and colitis, to alterations in the bacterial communities in the equine hindgut [2–5]. The term dysbiosis has been used to describe the changes in the microbial communities associated with disease states in other species [6, 7], a concept which appears relevant for the horse. While it is not always clear if dysbiosis is the cause or effect of disease, it is well established that the resident microbiota play a vital role in the health of the horse by providing energy and nutrient needs via the fermentation of plant products [2, 4, 5, 8]. For horses with colitis, the physiologic consequences of dysbiosis potentially include an overgrowth of pathogenic bacteria, bacterial toxin production, mucosal inflammation, alterations in the metabolism of carbohydrates, volatile fatty acids and bile acids [9].
Equine colitis is typically classified by etiology (e.g., infectious, antimicrobial-associated, non-steroidal anti-inflammatory-induced, inflammatory infiltrates, unknown, etc.) [10]. To date, there are few studies describing the microbial community structure associated with the various types of colitis. Antibiotic associated diarrhea (AAD) is a common sequelae of antimicrobial administration, and is associated with a higher risk of mortality than other forms of colitis [11]. Antimicrobial agents, which are commonly administered to horses to treat established infections and to provide prophylaxis for surgically created wounds, affect the microbiome of the hindgut of the horse. Currently, tetracyclines, macrolides, cephalosporins, fluoroquinolones, trimethoprim-sulphonamides, aminoglycosides, chloramphenicol, and β-lactams all have been reported to cause colitis in horses [3, 12].
Paradoxically, some antimicrobial agents, such as metronidazole, have been found useful in the treatment of colitis. Metronidazole is a bactericidal nitroimidazole with activity against anaerobic bacteria and protozoa. It is commonly used as a primary treatment for diarrhea in dogs [13], and has been used to treat foals with hemorrhagic enteritis caused by Clostridium difficile [14]. Metronidazole is also prescribed to Clostridial diarrhea in adult horses. Although historically metronidazole has been associated with a low incidence of diarrhea when used in adult horses [15], its effect on the bacterial communities of the GIT has not yet been reported. While culture-based methods have traditionally been used to assess the effects of antibiotics, the use of Next Generation Sequencing and metabolomics can potentially provide greater insight regarding the changes in the microbial populations and their functional effect on the metabolism of the horse’s hindgut [3, 16, 17]. The objective of this study was to characterize the changes in the cecal and fecal microbiome and metabolome of the horse before and after metronidazole administration.
Materials and methods
This study was approved by Texas A&M University Institutional Animal Care and Use Committee (IACUC; Protocol number 2014–0123).
Study population and sample collection
Five horses belonging to the Department of Animal Science at Texas A&M University with indwelling cecal cannulas were used as subjects. The horses had been cannulated approximately 10 years prior to this study, and had received no antibiotics or other medications (excluding routine vaccinations and anthelminthics) for the previous 12 months. The horses were housed on a dry lot, fed free choice coastal hay and a commercial pelleted concentrate (Nutrena Safe Choice, Cargill, Minnetonka, MN) at 0.5% of body weight for the duration of the study. Age, breed, sex, and weight of the study participants is summarized in Table 1.
Fecal samples were collected immediately after natural elimination. Cecal samples were collected by removing the plug of the cannula while horses were restrained in stocks and siphoning 5 mls of cecal ingesta, including both the particulate and liquid fractions. Samples were placed on ice for transport to the laboratory, where they were refrigerated at 4°C prior to DNA extraction. Samples were collected before (study days referred to as Dminus52, Dm28, Dm14, and 0) and after metronidazole administration (referred to as D7, D14, D28, and D52). Metronidazole (Unichem Pharmaceuticals, Hasbrouck Heights, NJ) was administered after sample collection on Day 0 at a dose of 15 mg/kg PO BID for 3 days instead of 7 days as planned in the study protocol due to development of complications (see below in result section for details). This resulted in an additional sampling point for feces on Day 3.
DNA extraction
One hundred mg of feces or cecal contents (particulate and liquid fractions) was aliquoted into a sterile 1.7 ml tube (Microtube, Sarstedt AG & Co, Numbrecht, Germany) containing 150 μl of 0.1 mm zirconia-silica beads and 100 μl of 0.5 mm zirconia-silica beads (BioSpec Products Inc., Bartlesville, OK, USA). Samples were then homogenized (FastPrep-24, MP Biomedicals, Irvine, CA, USA) for a duration of 1 minute at a speed of 4 m/s. DNA was extracted using the PowerSoil DNA Isolation Kit (MO BIO, Carlsbad, CA, USA) following the manufacturer’s instructions.
Sequencing of 16S rRNA genes
Sequencing of the V4 region of the 16S rRNA gene was performed at MR DNA (www.mrdnalab.com, Shallowater, TX, USA) on an Illumina MiSeq platform (Illumina Inc., San Diego, CA) following the manufacturer’s guidelines using 515F (5’-GTGCCAGCMGCCGCGGTAA-3’) and 806R (5’- GGACTACVSGGGTATCTAAT-3’). Briefly, the PCR reaction was performed in a single-step 30 cycle PCR using the HotStarTaq Plus Master Mix Kit (Qiagen, USA) under the following conditions: 94°C for 3 minutes, followed by 28 cycles (5 cycles used on PCR products) of 94°C for 30 seconds, 53°C for 40 seconds and 72°C for 1 minute, after which a final elongation step at 72°C for 5 minutes was performed. Using Illumina TruSeq DNA’s protocol, a DNA library was set up and Illumina MiSeq was used for sequencing according the manufacturer’s guidelines.
Analysis of sequences
QIIME v1.9 (Quantitative Insights into Microbial Ecology) was used for analysis of the sequences. After sequencing, barcodes and primers were removed and short (<150 bp), ambiguous, homopolymeric sequences were depleted from the dataset. USEARCH was used to identify and remove chimeric sequences. Operational taxonomic units (OTUs) were assigned based on at least 97% sequence similarity to the Greengenes database (v13.5) using an open reference approach. Sequences determined to be mitochondria, chloroplasts, unassigned, or those belonging to the phylum cyanobacteria were excluded from further analysis.
Enteric pathogen testing
Both cecal and fecal samples were tested for the presence of Salmonella and Clostridial toxins. For Salmonella PCR, a tetrathionate green broth enrichment method was used prior to extraction, and qPCR performed as previously described [18]. Commercially available ELISAs for Clostridium perfringens enterotoxin and Clostridium difficile toxins A and B were performed according to the manufacturer’s instructions (TechLab, Blacksburg, VA, USA).
Metabolomics analysis
Samples were stored at −80°C until shipped on dry ice to the West Coast Metabolomics Core (University of California, Davis, CA, USA) for untargeted analysis [19]. Samples were lyophilized for 24 hours and weighed in a 1.5 ml Eppendorf tube. One ml of a 3:3:2 extraction mixture of degassed acetonitrile (Fisher, Ottawa, Canada, A9554), isopropanol (Fisher, Ottawa, Canada, A461212) and water (Fisher, Ottawa, Canada, 7732-18-5) was used to re-suspend 4 mg of lyophilized sample. Three millimeter grinding beads (Next Advance, Troy, NY, USA) were added to disrupt the sample at 1500 rpm for 30 seconds followed by shaking at 4°C for 5 minutes and centrifugation for 2 minutes at 14,000 rcf. The supernatant was separated into two aliquots of 475 ul for dry down in a Centrivap cold trap vacuum overnight at room temperature (Labconco, Kansas City, MO, USA). Dried supernatant was resuspended in 500 ul of 50% aqueous acetonitrile (Fisher, Ottawa, Canada, A9554) and centrifuged at 14,000 rpm for 2 minutes. The supernatant was dried down overnight in a Centrivap overnight at room temperature. A volume of 10 ul of 40 mg/ml solution of methoxyamine hydrochloride (Sigma-Aldrich, St. Louis, MO, USA, 89803) in pyridine (99.99%) (Sigma-Aldrich, St. Louis, MO, USA 270970) was added to the samples and shaken for 90 minutes at 30°C. Next, 90 uL of N-methyl-N-trimethylsilyltrifluoroacetamide (MSTFA) with 1% trimethylchlorosilane (TMCS) (Sigma-Aldrich, St. Louis, MO, USA 69479) and Fatty Acid Methyl Esters (FAMEs) retention indexing markers (see next paragraph) were added with a ratio 100:0.001 was added to each sample and shaken for 37°C for 30 minutes for trimethylsilylation of acidic protons. The reaction mixture was placed in a 2 ml clear glass auto-sampler vial with micro-insert (Agilent, Santa Clara, CA, USA, 5185–5946) and closed with a 11mm T/S/T crimp cap (Thomas Scientific, Swedesboro, NJ, USA, 11-0038A).
Quality control was assessed by retention indexing with Fatty Acid Methyl Esters (FAMEs) and the use of blanks and quality control samples. A mixture of internal retention index markers using 13 FAME markers dissolved in chloroform (0.8mg/ml C8-C16 or 0.4 mg/ml C18-C30) (Acros, Morris Plains, New Jersey, USA, AC4235500100) was added to each sample (S1 File) as listed above during derivatization. During data acquisition, both blank samples and a quality control mix were injected every 10 samples. Blank samples were used to check for carry over and also in data processing for blank subtraction. The quality control mix included 28 compounds used in 6 concentrations (S1 File) that were dried down and derivatized as previously described.[19]
An Agilent 6890 GC equipped with a Gerstel automatic liner exchange system (ALEX) that includes a multipurpose sample (MPS2) dual rail, and a Gerstel CIS cold injection system (Gerstel, Muehlheim, Germany) with temperature program was used as follows: 50°C to 275°C final temperature at a rate of 12°C/s and hold for 3 minutes. Injection volume is 0.5 μl with 10 μl/s injection speed on a splitless injector with purge time of 25 seconds. Liner (Gerstel #011711-010-00) is changed after every 10 samples, (using the Maestro1 Gerstel software vs. 1.1.4.18). Before and after each injection, the 10 μl injection syringe is washed three times with 10 μl ethyl acetate. For gas chromatography, a 30 m long, 0.25 mm i.d. Rtx-5Sil MS column (0.25 μm 95% dimethyl 5% diphenyl polysiloxane film) with additional 10 m integrated guard column was used (Restek, Bellefonte, PA, USA). 99.9999% pure Helium with built-in purifier (Airgas, Radnor Pennsylvania, USA) was set at constant flow of 1 ml/min. The oven temperature was held constant at 50°C for 1 min and then ramped at 20°C/min to 330°C for 5 minutes. A Leco Pegasus IV time of flight mass spectrometer controlled by the Leco ChromaTOF software vs. 2.32 (St. Joseph, MI, USA) was used for mass spectrometry. The transfer line temperature between gas chromatograph and mass spectrometer was set to 280°C. Electron impact ionization at 70V was employed with an ion source temperature of 250°C. Acquisition rate was 17 spectra/second, with a scan mass range of 85–500 Da.
Raw data files were processed using ChromaTOF v.2.32. BinBase algorithm matched spectra to database compounds, and quantification was reported by peak height of an ion at the specific retention index characteristic of the compound across all samples. Peak heights were normalized by average total peak-sums for detected compounds across each sample group. Metabolites were cross referenced using their compound number to the KEGG database (https://www.genome.jp/kegg/pathway.html) in order to identify their metabolic pathways such as metabolism, cellular processes, and others. Metabolomic data has been submitted to metabolomicsworkbench.org under the submission ST001248.
Statistical analysis
Microbiome analysis.
Prior to analysis, data were tested for normality using the Shapiro-Wilk test (JMP Pro 14, SAS, Marlow, Buckinghamshire). As data followed a non-normal distribution, non-parametric measures were used throughout the study. Adjustments for multiple comparisons were made with either a Dunn’s post-test or Fischer’s least significant differences. P- and q-values <0.05 were considered statistically significant.
Alpha diversity was calculated using observed OTUs, Shannon, and Chao1 metrics to compare species richness and evenness. Statistical analysis of alpha diversity indices was performed using the software package PRISM (PRISM 7, GraphPad Software Inc., San Diego, CA). A Friedman’s test followed by a Dunn’s multiple comparison post-test were performed to assess differences in alpha diversity metrics between study days.
Beta diversity (bacterial community composition) was calculated using both weighted and unweighted UniFrac metrics to measure similarity between samples, and visualized for clustering with Principle Coordinate Analysis (PCoA) plots. An Analysis of Similarity test (ANOSIM) within the PRIMER 6 (PRIMER-E Ltd. Luton, UK) software package was performed on the beta diversity distance matrices to assess the significance of the differences in the bacterial community composition.
The abundance of bacterial taxa in the cecal and fecal samples was evaluated using a Friedman’s test (PRISM 7, GraphPad Software Inc., San Diego, CA) followed by a Dunn’s multiple comparison post-test. Only bacterial taxa present in at least 50% of the samples of at least one time point were included in the analysis. P values of <0.05 and q values of <0.1 were considered significant.
Linear discriminant analysis effect size (LEfSe) using the web-based program Calypso v8.62 (http://cgenome.net/wiki/index.php/Calypso) was performed to analyze the abundance of bacterial taxa and their associations with study day. A cut-off threshold of 3.5 was set for significance.
Metabolomics analysis.
MetaboAnalyst 4.0 (Xia Lab, McGill University, Canada) was used to analyze metabolomics data. The peak intensity data table contained peak heights normalized against the average total peak sums. Data was not filtered, normalized or transformed, but was subjected to Pareto scaling. Principle Component Analysis (PCA) plots were used to display metabolic composition of individual horses and a heat map was used to display hierarchical clustering of the significant metabolites. An ANOVA test followed by a Fisher’s least significant differences was used to determine, which time points in the cecal or fecal samples differed significantly over study days Dm28-D28.
Results
Study population and sample collection
Metronidazole administration began as planned on Day 0 after sample collection. By the evening treatment time on Day 3, all horses had become inappetent and had developed significant skin scalding associated with the canula site. Metronidazole was discontinued after a total of 5 doses had been administered (Day 3 of administration) due to the investigators’ concern of impending colitis. An extra fecal sample was collected on Day 3, while cecal samples were not obtained at this time point due to the horses’ level of discomfort with canula manipulation. The skin scalding was treated with local wound care including gentle washing of the skin and application of an emollient ointment around the canula site. All horses were normal on physical examination and their appetites had recovered within 48 hours of discontinuing metronidazole administration. Skin scalding resolved by Day 7. Horse 3 developed mild signs of colic (i.e., flank watching, pawing) on Day 13 and was moved from the research facility to the hospital for further diagnostics and monitoring. The horse’s vital parameters (i.e., temperature, heart rate, respiratory rate, mucous membrane color) remained within normal limits and it was producing normal feces. Palpation per rectum, passage of a nasogastric tube and preliminary bloodwork to assess hydration status (i.e., packed cell volume, total protein and lactate) were performed. No abnormalities were detected. Due to the mild nature of the colic, Horse 3 received 7 L of oral fluids with electrolyte supplementation and one dose of flunixin meglumine (1.1 mg/kg iv). The horse remained hospitalized for 24 hours during which time his vital parameters remained within normal limits, he was passing normal manure and the displayed no symptoms of colic. The investigator (CA) elected to return the horse to the research facility in the afternoon on Day 14. After scheduled sample collection, the horse was found deceased 4 hours later. The Institutional Animal Care and Use Committee was notified, and a necropsy was performed. The cause of death was determined to be due to typhylocolitis. Due to the death of this subject, data from Horse 3 is missing from the Day 28 and 52 time point data sets. The remaining 4 horses completed the study with no other adverse events. Due to these events, baseline samples include Days m52, m28, m14 and 0. Day3 represents the period when horses were treated with metronidazole and Days 7, 14, 28 and 52 represent the time period following discontinuation of metronidazole.
Sequence analysis
The total sequence analysis yielded 6,288,464 quality sequences for all analyzed samples (n = 81, mean ± SD = 77,635 ± 15, 531). For cecal content samples (n = 38), the mean quality sequences and standard deviation was 77,054 ± 18841. For fecal samples (n = 43), the mean quality sequences and standard deviation was 78,212 ± 12202. Samples were rarefied to an even depth of 40,005 reads per sample. The sequences were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the accession number SRP119693.
Beta diversity
Beta diversity (Fig 1A), as measured by Unifrac distances showed no significant differences for cecal samples (unweighted, R = 0.072, p = 0.108; weighted, R = 0.159, p = 0.001). Fecal samples (Fig 1B) showed visible clustering on Day 3 (ANOSIM, R = 0.152, p = 0.004), which was confirmed by pairwise comparisons. Beta diversity on Day 3 was significantly different from all other time points (Dm52, Dm28, Dm14, D0, D7, D14, D28, and D52) on both weighted and unweighted Unifrac distances in the fecal samples (Table 2).
(A) Cecal samples. (B) Fecal samples.
Alpha diversity
Alpha diversity indices (observed OTUs, Chao1, and Shannon) for cecal and fecal samples are reported in Fig 2A–2F. Due to the death of Horse 3 on Day 14, analysis by Friedman’s test on Days 28 and 52 could not be completed, but the data for the remaining 4 horses is provided. For cecal samples, the following significant differences were noted: Observed OTUs between Dm28 and D14 (p = 0.0007), Dm14 and D14 (p = 0.0007); Chao 1 between Dm28 and D14 (p = 0.001) and Dm28 and D7 (p = 0.001); and Shannon for Dm28 and D14 (p = 0.0021). For fecal samples, the following differences were noted: observed OTUs Dm52 and D3 (p = 0.015), Chao 1 between Dm14 and D3 (p = 0.008), and Shannon between D0 and D3 (p = 0.0209).
(A) Observed OTUs of cecal samples. (B) Chao1 of cecal samples. (C) Shannon of cecal samples. (D) Observed OTUs of fecal samples. (E) Chao1 of fecal samples. (F) Shannon of fecal samples.
Univariate analysis
In the cecal samples, only the phyla Firmicutes (increased in abundance from Dm52 to Dm28; p = 0026, q = 0.0208) and Tenericutes (decreased in abundance from Dm52 to D14; p = 0026, q = 0.0208) were significantly different between study days after adjustments for multiple comparisons had been made. Other phyla such as Fibrobacteres (p = 0.0154, q = 0.0778) and Planctomycetes (p = 0.0357, q = 0.1428) had significantly different p values, but did not achieve significance after adjustments for multiple comparisons were made (Fig 3A–3D). At the family level, there were 27 taxa that were significantly different on a Friedman’s test, but not after adjustment for multiple comparisons. Complete information regarding cecal content taxa at the phylum and family levels are included in S1 and S2 Tables.
(A) Fibrobacteres. (B) Firmicutes. (C) Planctomycetes. (D) Tenericutes.
In fecal samples, Lentisphaerae (p = 0.0035, q = 0.0438) and Spirochaetes (p = 0.0060, q = 0.05) were decreased on D3. Elusimicrobia (p = 0.0159, q = 0.0795) was significantly decreased on D14 compared to Dm28, but this difference was not evident once adjusted for multiple comparisons. Actinobacteria (p = 0.0117, q = 0.0731), Fibrobacteres (p = 0.0440, q = 0.1000), SR1 (p = 0.0212, q = 0.0850), Synergistetes (p = 0.0266, q = 0.0850), and Verrucomicrobia (p = 0.0387, q = 0.1000) were only significantly different on the Friedman’s but not the Dunn’s post-test (Fig 4A–4E). Complete information regarding fecal taxa at the phylum and family level is included in S3 and S4 Tables.
(A) Actinobacteria. (B) Elusimicrobia. (C) Lentisphearae. (D) Spirochaetes. (E) Verrucomicrobmia.
Linear discriminant analysis effect size
Linear discriminant analysis effect size (LEfSe) was used to elucidate taxa associated with study day. A linear discriminant analysis (LDA) score > 3.5 was considered significant. In the cecal content samples, only 2 phyla (i.e., Firmicutes and Fibrobacteres) met this threshold. At the family level in cecal samples, 5 taxa (i.e., Fibrobacteriaceae, Unclassified Y2, Ruminococcaceae, Unclassified Clostridiales, and Lachnospiraceae) were more abundantly expressed prior to metronidazole administration, whereas 3 taxa (i.e., Porphyromonadaceae, Veillonellaceae, and Succinivibrionaceae) were more abundantly expressed after metronidazole administration. Table 3 contains complete information regarding the taxa, study day, and LDA scores from cecal content samples.
Only LDA scores of >3.5 are shown.
In fecal samples, 2 taxa at the phylum level (i.e., Elusimicrobia and Euryarchea) and 3 taxa at the family level (i.e., Clostridicaceae, Methanocorpusculaceae, and Ruminococcaceae) were associated with an LDA score of greater than 3.5 at baseline. Fecal taxa at the phylum level significantly associated study days following metronidazole administration included: D3: Actinobacteria and Protobacteria; D14: Verrucomicrobia and Lentisphaerae; and D28: Planctommycetes and Synergistetes. At the family level, Alcaligenaceae, Corynebacteriaceae, Neisseriaceae, Actinomycetaceae, Porphyomonadaceae, Tisserellaceae, Enterobacteriaceae, Pastuerellaceae, Streptococcaceae, Aerococcaceae, Lacobacillaeceae, and Aeromondaceae were more abundant on D3, whereas Methanobacteriaceae and RFP12 were more abundant on D14. Bacillaceae and Planococcaceae were associated with D28. Table 4 displays the results of LEfSe analysis for fecal samples at both the phylum and family levels.
Only LDA scores of >3.5 are shown.
Testing for enteric pathogens
All horses were negative for Salmonella PCR in fecal samples across all time points. Four horses (i.e., horses 1, 2, 4, and 5) were positive for Salmonella by PCR testing in cecal content samples on Day 7. Horses 1, 2, and 4 remained positive on Day 14, but tested negative on Day 28. Horse 5 remained positive at D52, the final time point of the study. Horse 3 was positive on Day 14, and subsequently died later that day. Serotyping revealed Salmonella enterica serotype Newport in one horse (horse 1), Salmonella enterica serotype anatum in 3 horses (horses 2, 3, and 4) and multiple serotypes in one horse (horse 4). All horses were negative for Clostridium perfringens and C. difficile toxins in both cecal and fecal samples at all time points.
Metabolomic analysis of cecal and fecal samples
Using an untargeted approach, a total of 554 unique metabolites weredetected, of which 223 were named. Metabolomic data has been submitted to metabolomicsworkbench.org under the submission ST001248.
Only named metabolites were included in the analysis. Metabolites were examined in the cecum on Days m28, m14, 0, 7, 14 and 28. Metabolites were examined in the feces on Days m28, m14, 0, 3, 7, 14 and 28. Samples were not collected on Day m52 and 52 for metabolomic analysis, and on Day m28 Horse 5 did not have enough sample for analysis. Metabolites were analyzed using PCA score plots and heat maps. In the cecal samples, the PCA plot of all named metabolites (Fig 5A) indicated that samples on D14 were different from other study days. This is due to the presence of an outlier, Horse 3, whose metabolite profiles were presumably altered by the presence of gastrointestinal disease and subsequent death. No cecal metabolites were significantly different by study day after adjustments for multiple comparisons were made.
The shaded ellipses represent the 95% confidence interval of the metabolite profile. (A) Cecal metabolites. (B) Fecal metabolites.
In fecal samples, Day 3 samples clustered distinctly from all other study days in the PCA scores plot (Fig 5B). Mulitvariate analysis of all the named metabolites found 21 fecal metabolites that were significantly different across study days following adjustment for multiple comparisons. A heatmap of the significantly different named fecal metabolites (Fig 6) indicated a visible change in concentration relative to study day. Metabolites related to nucleic acid metabolism (thymine, uracil, xanthine), amino acid metabolism (putrescine, salicylic acid, serine, threonine, tryptophan, tyramine, valine), lipid metabolism (phosphoenolamine, ribonic acid) and carbohydrate metabolism (ribose) all had significantly elevated levels on Day 3. Other metabolites related to lipid metabolism (phytanic acid, stearic acid) and the metabolism of cofactors and vitamins (δ-tocopherol, γ-tocopherol) were decreased on Day 3 in relation to other study days. Information regarding the significant fecal metabolites and their respective KEGG pathways is provided in S5 Table. Information regarding the changes in significant fecal metabolites by study day is provided in S1–S5 Fig.
Each column represents a sample, sorted by study day, and each row represents one metabolite. Color intensity indicates intensity the increase (red) or the decrease (blue) in metabolite concentration.
Discussion
This study utilized a herd of horses with cecal cannulas to investigate the effects of metronidazole on the equine cecal content and fecal microbiome and metabolome. The horses had been cannulated approximately 10 years earlier, and had experienced no disruption in environment, feed, housing, or exposure to new horses for 6 months prior to and during the study period. While the use of cannulated animals is commonplace in nutritional studies, it is somewhat novel in microbiome and metabolome studies. The presence of a cannula in the cecum does allow ingresss of oxygen into the cecum during sampling periods. However, in other species it has been demonstrated that anaerobic conditions are quickly restored once the plug is replaced [20]. In ruminant species such as goats and cattle, the presence of the cannula has not significantly altered diversity indices or taxonomy of major phyla [21–23]. Because horses were sampled every 7–14 days, the short-term introduction of oxygen likely had a minimal effect. Although unlikely to significantly alter the results, the use of cannulated animals in this study should be interpreted accordingly.
The horses experienced adverse effects at two time points in the study, Days 3 and 14. The first adverse event occurred on Day 3 in all 5 horses and was manifested as inappetance and skin scalding at the canula site. Inappetance can be an early manifestation of gastrointestinal disease, and the skin scalding was presumably due to a change in the pH of the cecal fluid. Previously studies using these same 5 cannulated horses for grain overload studies indicated that decreasing the pH of the cecal fluid could result in fluid leaking from the canula and subsequent skin scalding [24]. Out of caution, the authors elected to suspend metronidazole administration after only 5 of the 14 doses had been given. All the horses returned to an apparently healthy state within 48 hours, and none of the horses showed symptoms of gastrointestinal disease until later in the study. The second adverse event occurred on Day 13 when Horse 3 experienced an episode of colic and died 24 hours later. Post-mortem examination indicated colitis as the cause of death. As the study had been suspended 11 days earlier and the remaining horses appeared healthy, the investigators chose to continue sampling from the remaining horses and use the data from the deceased horse. Due to the loss of Horse 3, samples from Days 28 and 52 only included horses 1, 2, 4, and 5. Friedman’s analysis for these two time points was not possible, although the data from these 4 horses is presented.
Similarly to reports in dogs, metronidazole decreased the diversity of the microbiome in the horses evaluated, but at different time points for cecal content and fecal samples [13, 25]. The OTU, Shannon, and Chao 1 metrics were consistently reduced in cecal samples on Day 14, whereas fecal samples had the lowest alpha diversity measures on Day 3. Recovery of alpha diversity indices occured by Day 28 in both sample types. Beta diversity was unaffected by metronidazole administration in the cecal samples. However fecal samples were distincly clustered on Day 3 compared to the other time points in the weighted and unweighted PCoA plots. Metronidazole had the greatest effects on alpha and beta diversity on Days 3 and 14, which is consistent with the clinical appearance of gastrointestinal disease in the study subjects.
At baseline, both cecal and fecal samples appeared to have a similar microbial community composition. Bacteroidetes and Firmicutes, important phyla for fiber degradation and affected by the presence of gastrointestinal disease in other studies [2, 3, 5, 17], constituted the majority of the phyla. Spirochaetes, Fibrobacteres, Proteobacteria, Tenericutes, and Verrucomicrobia comprised approximately 1% or more of the total bacteria. The cecum had a greater percentage of Bacteroidetes, while the feces contained more Fibrobacteres and Spirochaetes. This is consistent with other reports comparing the anatomic compartments of the GIT in normal horses [26, 27] and also in previous reports of cannulated horses [28].
Over time, the bacterial community composition was altered by metronidazole. In this study, the abundance of Bacteroidetes remained unchanged in either sample type. In the cecum, the phyla Firmicutes, Fibrobacteres, Planctomycetes, and Tenericutes were affected. There was a significant change in the abundance of these phyla from Dm52 to Dm28, which is likely related to normal variation in the microbiota over time, ambient temperature, and season [29–31]. The magnitude of this change likely overshadowed the downward trend in abundance of Fibrobacteres, Firmicutes, Planctomycetes, and Tenericutes on Days 7 and 14 due to metronidazole administration. The alterations in the Firmicutes occurred within the class Clostridia and the families Other, Unknown, Christensenellaceae, Dehalobacteriaceae, Lachnospiraceae, Ruminococcaceae, and Veillonellaceae. A similar effect was noted with the phylum Fibrobacteres and the family Fibrobacteraceae prior to metronidazole. These changes did not achieve significance after adjustments for multiple comparisons were made at the family level. This finding was reinforced by the LEfSe analysis. These changes may have reached significance if the number of horses used in the study would have been greater or if metronidazole administration had not been terminated early on Day 3.
In the feces, the phyla Actinobacteria, Elusimicrobia, Lentispherae, Spirochaetes, and Verrucomicrobia were affected by metronidazole administration. All of these phyla, except Actinobacteria, decreased in abundance on either Days 3, 7, or 14. The abundance of Actinobacteria was significantly elevated on Day 3 in the feces compared to other study days, with the family Bifidobacteriaceae accounting for majority of this increase. Actinobacteria is a gram-positive facultative anaerobe organism that accounts for a relatively small percentage of the total number of bacteria. However, a similar increase was reported after metronidazole use in dogs and rats [13, 32]. The significance of the decrease in the abundance Elusimicrobia, Lentisphereae, and Spirochaetes at various time points after metronidazole administration is unknown. Lentisphereae and Spirochaetes are anaerobes or facultative anaerobes, which were likely affected by metronidazole. Lentispherae is commonly found in the gut of mammals but comprises less than 1% of all taxa. Spirochaetes, however, represent 6.5% and 11.3% of the cecal content and fecal samples, respectively. Metronidazole caused a significant decrease in the Spirochaetes at both the phylum and family level on D3 but quickly recovered by the last sampling points. The functional role of these three phyla and their response to metronidazole are yet to be fully elucidated.
Univariate analysis and LEfSe indicated that Actinobacteria and Verrucomicrobia were affected on Day 3 in the feces. Verrucomicrobia showed a marked decline in both cecal content and fecal samples on Days 3 and 7, although this result did not reach significance with the Dunn’s post-test (q = 0.0717). Verrucomicrobia is a phylum of strict anaerobes that maintains the mucus layer between the gut lumen and the enterocytes. In some studies of healthy horses, Verrucomicrobia have been reported to account for up to 40% of the fecal microbiota [3]. In this study, Verrucomicrobia accounted for a much smaller percentage (i.e., 2.6–5.4%) of the total fecal microbiota at the four baseline sampling points, but even further decreased on Day 3 in fecal samples. This trend to decline after antibiotic use has been previously noted in both dogs and humans, and is thought to play a role in the loss of the intestinal barrier function in colitis.
All five of the horses in this study were PCR positive for Salmonella in their cecal content, but not in fecal samples, after metronidazole administration. Although not routinely tested for Salmonella before this study, the horses had never previously displayed symptoms of GIT disease. Identification of 3 serovars of Salmonella from the horses lends evidence against a herd outbreak of infectious Salmonella, which would typically include one serovar. All horses, except one, reverted to negative PCR status after the discontinuation of the metronidazole by the final study time point. Also, all horses were PCR negative at the 4 pre-treatment time points in both sample locations and remained negative in the fecal samples. Thus, the authors suggest that the metronidazole induced a degree of dysbiosis, which resulted in expansion of this enteric pathogen in the cecum, but not in the distal hindgut.
This study employed GC-MS methods to identify and quantitate changes in the end products of metabolism in the cecal and fecal samples before and after metronidazole administration. Of these 2 sample types, only the fecal metabolites were signficantly altered from baseline after adjustments for multiple comparisons were made. The fecal metabolites represented diverse metabolic pathways, such as nucleic acid, amino acid, carbohydrate, lipid and cofactor and vitamin metabolism. All significant metabolites belonging to the amino acid group were increased on Day 3. This could result from alterations in the commensal bacteria amd their role in the absorption or synthesis of these amino acids. Similar trends have been reported in dogs and humans with inflammatory bowel disease [33, 34]. Ribose, a metabolite of the pentose phosphate pathway in carbohydrate metabolism, was significantly elevated on Day 3, consistent with an oxidative stress response noted in the metabolomic profiles dogs with gastrointestintal disease [33]. Similarly, decreases in tocopherols, analogs of the anti-oxidant vitamin E, have been reported in horses suffereing from obesity [35].
Although metronidazole is reportedly safe to use in equine patients, the horses in this study appeared to develop early indications of colitis despite having normal feces. AAD is poorly defined in the veterinary literature in regards to the number or character of abnormal stools, but is generally regarded as a temporal association with the initiation or discontinuation of an antimicrobial agent and the development of diarrhea [9]. While diarrhea is characteristically considered pathognomonic for colitis, most horses exhibit prodromal symptoms associated with the gastrointestinal tract before diarrhea is clinically manifested. These symptoms often include inappetance, malaise, fever, and behavioral expressions of abdominal pain that precede the development of diarrhea. The clinical impression is that AAD is often acute in onset, occurring rapidly after the initiation of antimicrobial therapy. In horses, one study reported that the average time for development of diarrhea was 3.4 days after antibiotic administration (range: 1-11days) [36]. In humans, there is also evidence for the development of diarrhea associated with antibiotic discontinuation. This may have played a role in the death of Horse 3 on Day 14.
In this study, metronidazole decreased the diversity and altered the bacterial composition of cecal content and fecal samples. Subsequent functional alterations of the microbiome were reflected in the metabolite profile of the fecal samples. The timing of these changes coincides with the development of symptoms of GIT disease in these horses. Antibiotic administration, including metronidazole, is recognized as a risk factor for the development of diarrhea in species, such as humans, dogs, cats, and horses [3, 37–40].
Supporting information
S1 File. FAMEs, quality control samples used in GL-MS chromatography.
https://doi.org/10.1371/journal.pone.0232905.s001
(DOCX)
S1 Table. The abundance of cecal taxa at the phylum level by study day.
https://doi.org/10.1371/journal.pone.0232905.s002
(XLSX)
S2 Table. The abundance of cecal taxa at the family level by study day.
https://doi.org/10.1371/journal.pone.0232905.s003
(XLSX)
S3 Table. The abundance of fecal taxa at the phylum level by study day.
https://doi.org/10.1371/journal.pone.0232905.s004
(XLSX)
S4 Table. The abundance of fecal taxa at the family level by study day.
https://doi.org/10.1371/journal.pone.0232905.s005
(XLSX)
S5 Table. Significant fecal metabolites by study day and KEG pathway.
https://doi.org/10.1371/journal.pone.0232905.s006
(XLSX)
S1 Fig. Fecal metabolites in pathways of amino acid metabolism.
https://doi.org/10.1371/journal.pone.0232905.s007
(TIF)
S2 Fig. Fecal metabolites in pathways of nucleic acid metabolism.
https://doi.org/10.1371/journal.pone.0232905.s008
(TIF)
S3 Fig. Fecal metabolites in pathways of carbohydrate metabolism.
https://doi.org/10.1371/journal.pone.0232905.s009
(TIF)
S4 Fig. Fecal metabolites in pathways of lipid metabolism.
https://doi.org/10.1371/journal.pone.0232905.s010
(TIF)
S5 Fig. Fecal metabolites in pathways of the metabolism of cofactors and vitamins.
https://doi.org/10.1371/journal.pone.0232905.s011
(TIF)
References
- 1. Traub-Dargatz JL, Kopral CA, Seitzinger AH, Garber LP, Forde K, White NA. Estimate of the national incidence of and operation-level risk factors for colic among horses in the United States, spring 1998 to spring 1999. J Am Vet Med Assoc. 2001;219(1):67–71. Epub 2001/07/07. pmid:11439773.
- 2. Costa MC, Arroyo LG, Allen-Vercoe E, Stampfli HR, Kim PT, Sturgeon A, et al. Comparison of the fecal microbiota of healthy horses and horses with colitis by high throughput sequencing of the V3-V5 region of the 16S rRNA gene. PLoS One. 2012;7(7):e41484. Epub 2012/08/04. pmid:22859989; PubMed Central PMCID: PMC3409227.
- 3. Costa MC, Stampfli HR, Arroyo LG, Allen-Vercoe E, Gomes RG, Weese JS. Changes in the equine fecal microbiota associated with the use of systemic antimicrobial drugs. BMC Vet Res. 2015;11:19. Epub 2015/02/04. pmid:25644524; PubMed Central PMCID: PMC4323147.
- 4. Weese JS, Holcombe SJ, Embertson RM, Kurtz KA, Roessner HA, Jalali M, et al. Changes in the faecal microbiota of mares precede the development of post partum colic. Equine Vet J. 2015;47(6):641–9. Epub 2014/09/27. pmid:25257320.
- 5. Stewart HL, Southwood LL, Indugu N, Vecchiarelli B, Engiles JB, Pitta D. Differences in the equine faecal microbiota between horses presenting to a tertiary referral hospital for colic compared with an elective surgical procedure. Equine Vet J. 2019;51(3):336–42. Epub 2018/08/29. pmid:30153353.
- 6. Stewart AS, Pratt-Phillips S, Gonzalez LM. Alterations in intestinal permeability: the role of the "leaky gut" in health and disease. J Equine Vet Sci. 2017;52:10–22. Epub 2017/05/01. pmid:31000910; PubMed Central PMCID: PMC6467570.
- 7. Suchodolski JS. Diagnosis and interpretation of intestinal dysbiosis in dogs and cats. Vet J. 2016;215:30–7. Epub 2016/05/11. pmid:27160005.
- 8. Glinsky MJ SR, Spires HR, Davis CL. measurement of volatile fatty acid production rates in the cecum of the pony. J Anim Sci. 1976;42(6):1465–70. pmid:931822
- 9. McGorum BC, Pirie RS. Antimicrobial associated diarrhoea in the horse. Part 1: Overview, pathogenesis and risk factors. Equine Veterinary Education. 2009;21(11):610–6.
- 10. Feary DJ, Hassel DM. Enteritis and Colitis in Horses. Vet Clin North Am Equine Pract. 2006;22(2):437–79. https://doi.org/10.1016/j.cveq.2006.03.008. pmid:16882483
- 11. Cohen . Characteristics and risk factors for failure of horses with acute diarrhea to survive: 122 cases (1990–1996). J Am Vet Med Assoc. 1999;214(3):382–90. pmid:10023402
- 12. McGorum BC, Pirie RS. Antimicrobial associated diarrhoea in the horse. Part 2: Which antimicrobials are associated with AAD in the horse? Equine Vet Educ. 2010;22(1):43–50.
- 13. Olson E HJ, Waddle M, Steiner J, Suchodolski J, Gaschen F. Evaluation of the effects of a 2-week treatment with metronidazole on the fecal microbiome of healthy dogs. J Vet Intern Med. 2015;29(4):1184–5.
- 14. Weese JS, Staempfli HR, Prescott JF. A prospective study of the roles of Clostridium difficile and enterotoxigenic Clostridium perfringens in equine diarrhoea. Equine Vet J. 2001;33(4):403–9. Epub 2001/07/27. pmid:11469775.
- 15. Sweeney RW, Sweeney CR, Weiher J. Clinical use of metronidazole in horses: 200 cases (1984–1989). J Am Vet Med Assoc. 1991;198(6):1045–8. Epub 1991/03/15. pmid:2032913.
- 16. Gustafsson Båverud, Franklin Gunnarsson, Ögren Ingvast-Larsson. Repeated administration of trimethoprim/sulfadiazine in the horse–pharmacokinetics, plasma protein binding and influence on the intestinal microflora. J Vet Pharmacol Ther. 1999;22(1):20–6. pmid:10211713
- 17. Harlow BE, Lawrence LM, Flythe MD. Diarrhea-associated pathogens, lactobacilli and cellulolytic bacteria in equine feces: responses to antibiotic challenge. Vet Microbiol. 2013;166(1–2):225–32. Epub 2013/06/19. pmid:23769300.
- 18. Kurowski PB, Traub-Dargatz JL, Morley PS, Gentry-Weeks CR. Detection of Salmonella spp in fecal specimens by use of real-time polymerase chain reaction assay. Am J Vet Res. 2002;63(9):1265–8. pmid:12224858
- 19. Fiehn O, Wohlgemuth G, Scholz M, Kind T, Lee DY, Lu Y, et al. Quality control for plant metabolomics: reporting MSI-compliant studies. Plant J. 2008;53(4):691–704. pmid:18269577
- 20. Towne G, Nagaraja TG, Owensby C, Harmon D. Ruminal evacuation's effect on microbial activity and ruminal function. J Anim Sci. 1986;62(3):783–8. pmid:3084437
- 21. Lodge-Ivey SL, Browne-Silva J, Horvath MB. Technical note: Bacterial diversity and fermentation end products in rumen fluid samples collected via oral lavage or rumen cannula. J Anim Sci. 2009;87(7):2333–7. pmid:19329475
- 22. Terré M, Castells L, Fàbregas F, Bach A. Short communication: Comparison of pH, volatile fatty acids, and microbiome of rumen samples from preweaned calves obtained via cannula or stomach tube. J Dairy Sci. 2013;96. pmid:23706486
- 23. Wang L, Wu D, Yan T, Wang L. The impact of rumen cannulation on the microbial community of goat rumens as measured using 16S rRNA high-throughput sequencing. J Anim Physiol Anim Nutr (Berl). 2018;102(1):175–83. pmid:29057500
- 24. Warzecha CM, Coverdale JA, Janecka JE, Leatherwood JL, Pinchak WE, Wickersham TA, et al. Influence of short-term dietary starch inclusion on the equine cecal microbiome. J Anim Sci. 2017;95(11):5077–90. Epub 2018/01/03. pmid:29293739.
- 25. Igarashi H, Maeda S, Ohno K, Horigome A, Odamaki T, Tsujimoto H. Effect of oral administration of metronidazole or prednisolone on fecal microbiota in dogs. PLoS One. 2014;9(9):e107909. pmid:25229475
- 26. Costa MC, Silva G, Ramos RV, Staempfli HR, Arroyo LG, Kim P, et al. Characterization and comparison of the bacterial microbiota in different gastrointestinal tract compartments in horses. Vet J. 2015;205(1):74–80. Epub 2015/05/16. pmid:25975855.
- 27. Dougal K, Harris PA, Edwards A, Pachebat JA, Blackmore TM, Worgan HJ, et al. A comparison of the microbiome and the metabolome of different regions of the equine hindgut. FEMS Microbiol Ecol. 2012;82(3):642–52. Epub 2012/07/05. pmid:22757649.
- 28.
Jevit M. Microflora of the equine gut and its ramifications on the development of laminitis; A comparison of fecal and cecal diversity and Illumina and Roche 454 sequencers Pittsburgh, PA: Duquesne University; 2016.
- 29. Dougal K, Harris PA, Girdwood SE, Creevey CJ, Curtis GC, Barfoot CF, et al. Changes in the total fecal bacterial population in individual horses maintained on a restricted diet over 6 weeks. Front Microbiol. 2017;8:1502. Epub 2017/08/30. pmid:28848517; PubMed Central PMCID: PMC5554519.
- 30. Garcia-Mazcorro JF, Dowd SE, Poulsen J, Steiner JM, Suchodolski JS. Abundance and short-term temporal variability of fecal microbiota in healthy dogs. Microbiology Open. 2012;1(3):340–7. pmid:23170232
- 31. Martin EJ, Reed KJ, Kunz IGZ, Coleman RJ, Coleman SJ. Normal variation and changes over time in the equine intestinal microbiome. J Equine Vet Sci. 2017;52.
- 32. Pélissier M-A, Vasquez N, Balamurugan R, Pereira E, Dossou-Yovo F, Suau A, et al. Metronidazole effects on microbiota and mucus layer thickness in the rat gut. FEMS Microbiol Ecol. 2010;73(3):601–10. pmid:20579100
- 33. Minamoto Y, Otoni CC, Steelman SM, Büyükleblebici 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. 2015;6(1):33–47. pmid:25531678
- 34. Bjerrum JT, Wang Y, Hao F, Coskun M, Ludwig C, Günther U, et al. Metabonomics of human fecal extracts characterize ulcerative colitis, Crohn's disease and healthy individuals. Metabolomics. 2015;11:122–33. Epub 2014/06/01. pmid:25598765.
- 35. Coleman MC, Whitfield-Cargile CM, Madrigal RG, Cohen ND. Comparison of the microbiome, metabolome, and lipidome of obese and non-obese horses. PLoS One. 2019;14(4):e0215918. pmid:31013335
- 36. Wilson DA, MacFadden KE, Green EM, Crabill M, Frankeny RL, Thorne JG. Case control and historical cohort study of diarrhea associated with administration of trimethoprim-potentiated sulphonamides to horses and ponies. J Vet Intern Med. 1996;10(4):258–64. pmid:8819052
- 37. Costa MC, Weese JS. Understanding the intestinal microbiome in health and disease. Vet Clin North Am Equine Pract. 2018. Epub 2018/02/07. pmid:29402480.
- 38. AlShawaqfeh MK, Wajid B, Minamoto Y, Markel M, Lidbury JA, Steiner JM, et al. A dysbiosis index to assess microbial changes in fecal samples of dogs with chronic inflammatory enteropathy. FEMS Microbiol Ecol. 2017;93(11):fix136-fix. pmid:29040443
- 39. Honneffer JB, Minamoto Y, Suchodolski JS. Microbiota alterations in acute and chronic gastrointestinal inflammation of cats and dogs. World J Gastroenterol. 2014;20(44):16489–97. Epub 2014/11/28. pmid:25469017.
- 40. Bartlett JG. Antibiotic-associated diarrhea. N Engl J Med. 2002;346(5):334–9. pmid:11821511.