Next Article in Journal
Dietary Cold Pressed Watercress and Coconut Oil Mixture Enhances Growth Performance, Intestinal Microbiota, Antioxidant Status, and Immunity of Growing Rabbits
Previous Article in Journal
Using Radio-Frequency Identification Technology to Measure Synchronised Ranging of Free-Range Laying Hens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monensin Alters the Functional and Metabolomic Profile of Rumen Microbiota in Beef Cattle

College of Agriculture, Communities, and the Environment, Kentucky State University, Frankfort, KY 40601, USA
*
Author to whom correspondence should be addressed.
Submission received: 17 October 2018 / Revised: 13 November 2018 / Accepted: 14 November 2018 / Published: 17 November 2018
(This article belongs to the Section Cattle)

Abstract

:

Simple Summary

Monensin can enhance the efficiency of feed utilization by modulating rumen fermentation; however, its effects on rumen function has not been fully described. Thus, this study integrated metagenomics and metabolomics analysis to identify differences in functional attributes and metabolites of rumen microbiota in beef steers fed no or 200 mg/d of monensin. Our results showed differences in relative abundance of functional genes involved in lipid metabolism and amino acid metabolism as well as changes in rumen fluid metabolites and their metabolic pathways. This study revealed a better understanding of the effects of monensin, which may enable more effective use of this additive for beef cattle production.

Abstract

To identify differences in rumen function as a result of feeding monensin to beef cattle, rumen fluid metagenomics and metabolomics analyses were used to evaluate the functional attributes and metabolites of rumen microbiota in beef steers fed no or 200 mg/d of monensin. Eight rumen-fistulated steers were used in the study for a period of 53 days. Rumen fluid samples were collected on the last day of the experiment. Monensin increased the relative abundance of Selenomonas sp. ND2010, Prevotella dentalis, Hallella seregens, Parabacteroides distasonis, Propionispira raffinosivorans, and Prevotella brevis, but reduced the relative abundance of Robinsoniella sp. KNHs210, Butyrivibrio proteoclasticus, Clostridium botulinum, Clostridium symbiosum, Burkholderia sp. LMG29324, and Clostridium butyricum. Monensin increased the relative abundance of functional genes involved in amino acid metabolism and lipid metabolism. A total of 245 metabolites were identified. Thirty-one metabolites were found to be differentially expressed. Pathway analysis of the differentially expressed metabolites revealed upregulated metabolic pathways associated with metabolism of linoleic acid and some amino acids. These findings confirm that monensin affects rumen fermentation of forage-fed beef cattle by modulating the rumen microbiome, and by reducing amino acid degradation and biohydrogenation of linoleic acid in the rumen.

1. Introduction

The rumen microbiota plays a central role in the efficiency of digestion in ruminants [1], thus, the use of rumen fermentation modifiers, such as monensin, that can enhance the efficiency of feed utilization by increasing the amount of energy available has been the focus of research for several years [2]. Monensin is a carboxylic polyether ionophore fed to ruminants to modify rumen fermentation dynamics by selectively inhibiting growth of gram-positive bacteria, which produce most of the acetate, lactate, and hydrogen in the rumen [3,4]. This favors growth of gram-negative bacteria and production of propionate in the rumen [5]. Increased production of propionate in the rumen increases hepatic gluconeogenic flux [6], which improves the overall energy status of ruminants [7]. Monensin also decreases amino acid fermentation and ruminal ammonia concentration in the rumen [8], however, despite having a basic understanding of the effects of monensin and the significance of rumen microbiota, little is known about the effects of monensin on the genetic, functional, and metabolomic attributes of the rumen microbiota.
In recent years, the availability of high-throughput comparative metagenomics enabled by development of next-generation sequencing platforms has advanced our understanding of the composition and function of microbial populations in diverse environments [9,10]. Microbial characterization via 16S rRNA has been used to assess the effects of monensin on ruminal bacterial diversity in dairy and beef cattle [4,11]. However, this technique is limited because of the bias induced by PCR while amplifying the gene and it offers limited taxonomical and functional resolution; in general, OTUs are analyzed at the genera level, and can be less precise at the species level [12,13]. Shotgun metagenomic sequencing avoids PCR bias and provides a genetic and functional perspective of the microbiome [14]. This technique has already been used to reveal functional genes and microbiome associated with performance and disease conditions of cattle [15,16]. Application of metabolomics analysis provides an opportunity to measure large numbers of small molecule metabolites in cells, tissues and biofluids comprehensively [17]. Recent studies have applied metabolomics to predict feed efficiency and residual feed intake [18,19], examine disease conditions [20], evaluate dietary responses to different feeds [21], and assess milk quality of ruminants [22].
To date, the influence of dietary monensin on functional and metabolomic attributes of the rumen microbiota has not been studied in beef cattle. Therefore, the objective of this study is to determine the functional analysis of the rumen microbiota and its metabolomic attributes in forage-fed Holstein steers fed no or 200 mg d−1 of monensin. In this study, we applied ultra-performance liquid chromatography (UPLC)/time-of-flight mass spectrometry (MS) and multivariate/univariate statistical analysis to characterize rumen metabolites and whole-metagenomic shotgun sequencing to provide a snapshot of the rumen microbial population.

2. Materials and Methods

The research protocol (protocol number 18-001) was reviewed and approved by the Institutional Animal Care and Use Committee of Kentucky State University.

2.1. Animals, Housing, and Feeding

The trial consisted of a 14-d adaptation to diet and a 39-d treatment period. On day 0, eight rumen-fistulated Holstein steers (593 ± 24 kg body weight) maintained at the Kentucky State University Research Farm were blocked by body weight and assigned randomly to one of two treatments for 39 days: no supplementation (control) and monensin supplementation (200 mg hd−1 d−1). The steers were housed in individual pens and had ad libitum access to water and red clover/orchard grass hay (Table S1). Mineral mix (Hubbard feeds, Mankato, MN, USA) was fed free choice. Concentrate supplement containing corn gluten meal, soy hull, and cracked corn, with no or 200 mg of monensin (Elanco Animal Health, Greenfield, IN, USA) was supplemented at 4 kg/steer daily (Table S1). Steers were fed once daily at 0900 h.

2.2. Rumen Fluid Collection

Representative samples (300 mL) of the ruminal contents were collected via the cannula by spot sampling from the midpoint along the length of the ruminal contents and from the midpoint along the height of the ruminal contents (ventral rumen) at approximately 3, 6, and 9 h after feeding on the last day of the experiment. At the time of collection, ruminal contents were tightly hand-strained through 4 layers of sterile cheesecloth to separate the liquid and solid samples. Daily composited samples (solid and liquid samples) were mixed 1:1 (w/w) and stored at −80 °C for subsequent metagenomics and metabolomics analysis.

2.3. Whole-Genome Shotgun Sequencing

2.3.1. DNA Extraction and Sequencing

Frozen rumen fluid samples were thawed at room temperature and centrifuged at 15,000 × g followed by removing the supernatant. The DNA was then extracted and purified from the pellets using a PowerSoil DNA Isolation Kit (MO BIO Laboratories Inc., Carlsbad, CA, USA) according to the manufacturer’s instructions. The integrity of the DNA was verified by agarose (0.7%) gel electrophoresis, and the DNA was stored at −20 °C until further use. The sequencing library was constructed by DNA fragmentation following manufacturer’s instructions (Illumina) [23]. Following Library preparation, the sheared DNA fragments were then sequenced on an Illumina HiSeq 2500 platform by an Illumina HiSeq–PE150 bp strategy.

2.3.2. Bioinformatics and Statistical Analysis

After sequencing, the raw sequence data files were demultiplexed and stored as fastq format. The sequence data were assembled using IDBA-UD assembler [24] following quality control such as removal of reads containing barcode and adaptors, low-quality reads, and contaminative reads. MetaGene (http://metagene.cb.k.u-tokyo.ac.jp/) was employed to predict open reading frames from assembled contigs [25], and gene sequence with length more than 100 bp were kept and translated into amino acid sequences. All the predicted genes were then put together and clustered by CD-HIT software (http://www.bioinformatics.org/cd-hit/) with the parameters (95% identity, 90% coverage) [26]. The longest gene was used as a representative sequence to construct a non-redundant gene set. The non-redundant multi-source protein annotation database (M5NR [27]) was used for phylogenetic classification. Based on Kruskal–Wallis (KW) sum-rank test, the linear discriminant analysis effect size (LEfSe) method was used to identify the most differentially abundant taxonomic features at genus and species levels [28]. For this analysis, the significance threshold for the Kruskal–Wallis (KW) test was set to 0.05 and the logarithmic linear discriminant analysis (LDA) score cut-off was set to 2.0. The set of genes was aligned with KEGG (http://www.genome.jp/kegg/) gene database using BLAST to obtain functional annotation information [29]. Gene sequences were aligned with the CAZy database (http://www.cazy.org/) to get information on the functional classification of carbohydrate-active enzymes [30]. Differences in the functional features of the two metagenomes were analyzed using Mann Whitney test using p-value ≤ 0.05.

2.4. Non-Targeted Metabolomics Analysis

2.4.1. Sample Preparation and Analysis

The rumen fluid samples were thawed on ice at room temperature approximately 2 h before use. The samples were prepared using the procedure of Want et al. [31]. Briefly, 500 µL of the ruminal fluid samples were mixed with 2 mL of methanol-water (1:1, v/v) and then vortex-mixed for 2 min. Subsequently, the mixture was centrifuged at 15,000× g at 4 °C for 10 min. The supernatant was dried in a vacuum concentrator, and re-suspended in 200 µL methanol/water (1:1 vol/vol). The analysis was done using an UltiMateTM 3000 ultra-performance liquid-chromatography (UPLC) system (Thermo Fischer Scientific, Waltham, MA, USA) equipped with an autosampler and coupled with an Orbitrap-Velos mass spectrometer (MS). Chromatographic separation was carried out with an Agilent Extend C-18 column (3.0 × 150 mm, 3.5 µm; Agilent) with a temperature of 45 °C and sample manager temperature of 4 °C. The UPLC mobile phases were (A) 0.1% formic acid-water and (B) 0.1% formic acid-acetonitrile. The injection volume was 5 µL, and the flow rate was 0.5 mL/min. Mass spectrometry was performed in both positive and negative modes. The capillary voltage was 3.5 kV for both positive and negative modes. Capillary and source-heat temperature was set at 350 °C, sheath gas flow at 40 L h−1, and auxiliary gas flow at 10 L h−1. Quality control (QC) sample was prepared by mixing 10 µL of each rumen fluid sample. The QC was run in positive and negative modes every 4 samples to serve as technical replicates for validating stability and reproducibility of UPLC/MS system.
For analysis of volatile fatty acids (acetate, propionate, and butyrate), 12 μL of 50% H2SO4 were added to 12 mL of the liquid portion of ruminal content, and the mixture was centrifuged at 11,500× g for 20 min. The supernatant was stored at −20 °C until analyzed for volatile fatty acids using a Merck Hitachi Elite La-Chrome High-Performance Liquid Chromatograph system (Hitachi L2400, Tokyo, Japan) fitted with a Bio-Rad Aminex HPX-87H column (Bio-Rad Laboratories, Hercules, CA, USA) with a 0.015 M sulfuric acid mobile phase and a flow rate of 0.7 mL/min at 50 °C [32].

2.4.2. Data Processing and Statistical Analysis

The raw data were converted to Analysis Base File format by Reifycs ABF Converter (http://www.reifycs.com/AbfConverter/index.html). Processing of the raw data, including peak picking, deconvolution, compound identification, and peak alignment was done using the MS-DIAL version 2.84 [33] based on the mass to charge (m/z) value and the retention time of the ion signals. Identified metabolites from both positive and negative modes were merged and imported into the MetaboAnalyst 4.0 [34] for multivariate/univariate analysis. After normalization, Principal Components Analysis (PCA) was first used as an unsupervised method for data visualization and outlier identification. Supervised regression modeling was then performed by the use of orthogonal partial least squares discriminant analysis (OPLS-DA) to identify the significantly differential metabolites. The differential metabolites were filtered and confirmed by combining the results of t-test (p ≤ 0.10) and fold change (FC) of the peak intensities (mean value of peak intensity obtained from monensin group/mean value of peak intensity obtained from Control group). The chemical structures of differential metabolites were identified according to online databases such as the Human Metabolome Database (www.hmdb.ca), Metlin (www.metlin.scripps.edu) and the Mass Bank (www.massbank.jp) using the data of accurate masses and MS/MS fragments. Pathway analysis was conducted using MetaboAnalyst 4.0 software with a Bos taurus pathway library using hypergeometric test for over-representation analysis and relative-betweeness centrality for pathway topology analysis.

3. Results

3.1. Rumen Fluid Metagenomics Profiling

3.1.1. Sequencing Results

Metagenome sequencing of rumen fluid samples from the 8 Holstein beef steers yielded approximately 24.1 million reads per sample, which resulted in about 84,033 contigs per sample. Rarefaction analysis, which was used to assess the depth of sequencing, showed that the number of sequences used for all samples was sufficient to determine the total number of sequence types (Figure S1).

3.1.2. Phylogenetic Profile

Sequences from the rumen fluid samples predominantly contained phylotypes affiliated with Bacteroidetes (39.4 ± 5.3%), Firmicutes (15.9 ± 2.0%), Proteobacteria (1.9 ± 0.9%), Actinobacteria (0.6 ± 0.08%), and Euryarchaeota (1.0 ± 0.9%) (Table S2). Though no statistical difference was found, the relative abundance of Bacteroidetes was greater in steers fed monensin (43.5 ± 3.5% vs. 35.3 ± 3.1%) (Table S2). At the genus level, 1541 taxa were detected. Prevotella (24.7 ± 5.05%) was the most predominant, followed by Bacteroides (6.63 ± 1.38%) (Table S3). At the species level, 7105 taxa were detected (Table S4). The most predominant were Prevotella sp. FD3004 (6.41 ± 1.96%), followed by Prevotella ruminicola (2.93 ± 0.58%), Prevotella sp. MA2016 (1.86 ± 0.47%), and Prevotella brevis (1.70 ± 0.30%) (Table S4). To identify differentially abundant taxa that were mostly affected by treatment with monensin, we compared the rumen microbial population of both treatments using a metagenomic biomarker discovery approach, LEfSe, which performs a nonparametric Wilcoxon sum-rank test followed by linear discriminant analysis to assess the effect size of each differentially abundant taxon [34]. Using LEfSe, no treatment effects were found at the phylum level. At the genus level, the relative abundance of Mitsuokella, Hallella, and Propionispira was enriched in steers fed monensin, while Streptococcus, Sphaerochaeta, Burkholeria, Lachnoanaerobaculum, Terriglobus, Fusobacterium, and Methanobacterium were reduced (Figure 1A). At the species level, the relative abundance of 616 taxa were affected by treatment (data not shown). However, LEfSe analysis of the 500 most predominant species revealed that the relative abundance of Selenomonas spp. ND2010, Prevotella dentalis, Hallella seregens, Parabacteroides distasonis, Propionispira raffinosivorans, and Prevotella brevis were enriched in steers fed monensin, while the relative abundance of 12 taxa, such as Robinsoniella sp. KNHs210, Butyrivibrio proteoclasticus, Clostridium botulinum, Clostridium symbiosum, Burkholderia sp. LMG29324, and Clostridium butyricum, reduced (Figure 1B).

3.1.3. Functional Profile of the Ruminal Microbial Community

Using KEGG orthology level 2, the most predominant genes were those involved in metabolism, including amino acid metabolism, carbohydrate metabolism, glycan biosynthesis and metabolism, lipid metabolism, metabolism of cofactors and vitamins, xenobiotics biodegradation and metabolism, and nucleotide metabolism (Figure 2A). Other predominant genes were those involved in membrane transport, signal transduction, replication and repair, translation, and folding, sorting, and degradation (Figure 2A). The relative abundance of functional genes involved in amino acid metabolism, nucleotide metabolism, and lipid metabolism were enriched (p < 0.05) in the metagenome of steers fed dietary monensin (Figure 2B).
The genes aligned to the CAZy database were categorized into 6 types, auxiliary activities, carbohydrate-binding modules, carbohydrate esterases, glycoside hydrolases, glycoside transferases, and polysaccharide lyases. Carbohydrate hydrolases were the most enriched, followed by glycosyl transferases and carbohydrate esterases (Figure 3). No effects of dietary treatment on the relative abundance of carbohydrate-active enzymes were detected.

3.2. Rumen Fluid Metabolomics Profiling

Quality control samples were used to demonstrate the stability of the UPLC/MS system. The clustering of QC samples in the PCA scores plot (Figure S2) demonstrated satisfactory reproducibility and stability of the metabolic features and were subsequently used for statistical analysis.
In total, 245 metabolites were identified (Table S5). Principal component analysis revealed different distributions of the two treatments. PC1 explained 93.4% and PC2 explained 3.4% of the variance (Figure 4A). Score plot from OPLS-DA, a powerful tool used for dimension reduction and identification of spectral features that drive group separation [35], revealed a clear separation of the two treatments (Figure 4B), indicating that feeding dietary monensin altered the ruminal metabolome. Permutation analysis (p-value = 0.05, R2 = 0.99, Q2 = 0.42) of the OPLS-DA confirmed the validity of the model.
Thirty-one metabolites were differentially expressed (p ≤ 0.10; Table 1). Feeding dietary monensin increased (p ≤ 0.10; FC > 1.0) the relative concentrations of 24 metabolites, such as linoleic acid, l-phenylalanine, hypoxanthine, oxysporidinone, l-β-homomethionine, 2-methylgluric acid, isopongaflavone, catechol, pimelic acid, l-Histidine, tri-O-methylgenistein, oxysporidinone, formylindole, indole-3-carboxylic acid, and 8,15-DiHETE. The relative concentrations of acetate, 4-aminophenol, tyramine, 3-(4,4-dimethyl-4,5-dihydro-1,3-oxazol-2-yl)pyridine, 4-pyrodoxic acid, 1,3-dipropyl-7-methylxanthine, pentalenolactone, and vildagliptin, and apigenin were reduced (p ≤ 0.10; FC < 1.0) in steers fed monensin. Results of volatile fatty acid profile revealed that monensin reduced acetate (52.3 vs. 60 mM; p ≤ 0.05), but increased propionate (27.5 vs. 23.3 mM) compared with the control.
Pathway analysis of the differentially expressed metabolites revealed 14 associated metabolic pathways. Metabolic pathways associated with linoleic acid metabolism (impact value = 1.0), phenylalanine, tyrosine and tryptophan biosynthesis (impact value = 0.5), phenylalanine metabolism (impact value = 0.41), and histidine metabolism (impact value = 0.27) were upregulated in steers fed monensin (Table 2).

4. Discussion

The rumen microbiota plays a central role in the efficiency of digestion in ruminants [36]. Thus, determining the functional and metabolomic attributes of the rumen microbiome is essential for understanding their role on host metabolism and health [36,37]. Manipulation of the rumen microbiota via dietary intervention, such as the use of monensin with the aim of improving the host performance has been an active area of research for several years [38]. In this study, we used shotgun metagenomics sequencing and KEGG gene annotation database to evaluate the functional potential of the rumen microbiota, as well as UPLC/MS method and pattern recognition approaches, such as PCA and OPLS-DA to identify differences in rumen metabolites in beef cattle fed no or 200 mg d−1 of monensin. To the best of our knowledge, this is the first report evaluating the effects of dietary monensin on the functional and metabolomic attributes of the rumen microbiota in forage-fed beef cattle.
It is notable that the rumen microbial composition and the relative abundance of phyla is influenced by diet [16,39]. Similar to other studies [39,40], the predominant phyla identified in this study include Bacteroidetes, Firmicutes, Proteobacteria, and Euryarchaeota. The numerically higher relative abundance of gram-negative Bacteroidetes and lower relative abundance of gram-positive Firmicutes in steers fed monensin in this study is in accordance with the mode of action of monensin, which acts by reducing the population of gram-positive bacteria in the rumen. A recent study that characterized the gut microbiome of feedlot cattle supplemented with an antibiotic feed additive containing tylosin and monensin had similar results [16]. In addition, a study that investigated the ruminal bacterial diversity in dairy cattle fed monensin indicated the proportion of Firmicutes was reduced while that of Bacteroidetes was increased [11].
We used the LEfSe method to identify those taxa that are affected by the dietary treatment at the genus and species levels. LEfSe combines statistical significance, biological relevance, and effect size using linear discriminant analysis to identify genomic features such as genes, pathways, or taxa characterizing the differences between two or more biological conditions [28]. In this study, both gram-positive and gram-negative taxa were reduced in response to monensin. This agrees with the suggestion that the presence or absence of an outer membrane is not the only factor that determines susceptibility of a bacterium to monensin [11]. Future studies should evaluate other factors that may alter the effect of monensin on rumen bacteria such as chemicophysical features of solid feed particles and the ability to adapt to ionophores [11].
It should be noted that all of the taxa enriched in response to monensin treatment are gram-negative, confirming that the overall effects of monensin favor the growth and activities of gram-negative bacteria. Prevotella sp., Hallella seregens, and Parabacteroides distasonis produce succinate [41,42,43], Selenomonas species have the ability to decarboxylate succinate to form propionate [44], and Propionispira raffinosivorans is a major propionic acid-producer [45]. The shift in the microbial population in favor of the aforementioned bacteria supports the increased propionate synthesis observed in this study as a result of feeding monensin. Similar to the results of this study, several studies have reported the effectiveness of dietary monensin at altering the production and proportion of rumen fermentation acids. A recent study reported increased propionate concentrations in forage-fed beef steers consuming 200 mg d−1 of monensin [46]. Similarly, beef cattle consuming low-quality dry winter range grass fed 200 mg d−1 of monensin had a decrease in the molar percentage of acetate and an increase in the molar percentage of propionate in [47]. The relative proportion of fermentation acids produced in the rumen is significant due to the role of propionate in energy metabolism [7]. In ruminants, propionate is the most predominant substrate for gluconeogenesis, a major pathway for maintaining adequate glucose supply for ruminants [48]. Therefore, increased propionate production in the rumen as a result of feeding monensin is beneficial for cattle.
Another aspect of monensin inclusion worthy of attention is its effect on ruminal methane production. It is hypothesized that monensin can reduce methane production by inhibiting the growth of the bacteria responsible for supplying methanogens with the primary substrates, H2 and formate, for methanogenesis [49]. The reduced abundance of Methanobacterium by monensin in this study may be as a result of reduced availability of H2 and/or differences in growth kinetics of methanogen species because ionophores are not known to inhibit methanogens directly, but can change the population dynamics of methanogen species [50].
In a comprehensive study of the bovine ruminal fluid metabolome, Saleem et al. [21] combined nuclear magnetic resonance spectroscopy, inductively coupled plasma mass-spectroscopy, gas chromatography-mass spectrometry, direct flow injection mass spectrometry and lipidomics with computer-aided literature mining to identify and quantify about 246 ruminal fluid metabolites. In accordance with this study, the composition of bovine ruminal fluid is dominated primarily by microbial metabolites, including short-chain fatty acids, phospholipids, dicarboxylic acids, amino acids, and triglycerides [21]. The alteration in the rumen metabolome we observed as a result of feeding monensin was not surprising since alteration of the rumen microbial population generally influences the types of compounds produced by the rumen microbiota. In this study, the relative concentrations of 31 of the total 247 metabolites identified were altered by feeding monensin.
Pathway analysis of the metabolites expressed differentially in monensin-treated cattle, based on t-test and FC analysis, revealed fourteen pathways were affected: purine metabolism, tyrosine metabolism, vitamin B6 metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, nitrogen metabolism, linoleic acid metabolism, biosynthesis of unsaturated fatty acids, aminoacyl-tRNA biosynthesis, histidine metabolism, glycolysis or gluconeogenesis, propanoate metabolism, pyruvate metabolism, and caffeine metabolism. However, only four pathways (phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, histidine metabolism, and linoleic acid metabolism) had pathway impact values greater than 0.1, which represents the cut-off for relevance [51]. All of the relevant pathways were upregulated in steers fed monensin due to increased concentrations of the associated metabolites: phenylalanine, histidine, and linoleic acid. This is in agreement with the KEGG functional gene annotation of the metagenome, which revealed higher relative abundance of functional genes involved in amino acid metabolism and lipid metabolism in steers fed monensin.
The increased abundance of functional genes involved in lipid metabolism and the upregulation of the linoleic acid metabolism pathway in steers fed monensin confirm results from previous studies that suggest monensin can inhibit ruminal biohydrogenation of unsaturated fatty acids [52]. As shown in the current study, several lines of evidence indicate monensin can inhibit the growth of Butyrivibrio sp., a group of gram-positive bacteria that can biohydrogenate linoleic acid in the rumen [44,53]. Relative abundance of Butyrivibrio proteoclasticus, which is known to biohydrogenate conjugated linoleic acid (CLA) and vaccenic acid in the rumen [54], was reduced in the current study. Fellner et al. [52] observed higher concentrations of vaccenic acid and CLA in continuous culture of bacteria after treatment with monensin. In Korean native steers fed concentrated feed supplemented with monensin, concentrations of cis-9, trans-11 CLA were increased [55]. Consequently, ruminants receiving dietary monensin have higher concentrations of CLA in adipose tissue and milk [56,57]. Odongo et al. [58] and Silva-Kazama et al. [59] reported increased concentration of CLA in milk of dairy cows fed monensin. Likewise, the use of 230 mg/d of monensin increased the concentration of α-linolenic acids in subcutaneous fat of young bulls [60]. Due to the numerous health benefits of CLA for humans, nutritional strategies that increase the CLA content of animal products are considered desirable [58,61].
As shown in this study, monensin can reduce amino acid degradation and ammonia accumulation in the rumen [62,63] due to suppression of hyper-ammonia-producing bacteria, Peptostreptococcus and Clostridium, in the rumen [64,65,66]. A significant reduction of amino acid degradation was observed in mixed and pure cultures of rumen bacteria treated with monensin [67]. In another study, monensin decreased ammonia concentrations and increased the flow of amino acids from the rumen in cows fed timothy hay supplemented with casein, gelatin, or soy hydrolysate [66]. The reduced abundance of amino acid fermenting bacteria as well as upregulation of phenylalanine, tyrosine, and tryptophan biosynthesis, phenylalanine, and histidine metabolism pathways observed in this study suggests that monensin can reduce wasteful degradation of amino acids in the rumen.
Metagenomics approach requires sophisticated computational methods which are reliant on the integrity and presence of sequence data within established databases that may be subject to low taxonomic resolution [68,69]. Non-targeted metabolomics approach relies on comparing peak intensity to evaluate differences in relative abundance of metabolites which may often lack accuracy and precision [70]. Also, it is still a challenge to accurately identify metabolites due to the complexity and chemical diversity of the metabolome [71]. Another limitation of this study is the limited number of animals used due to high cost of shotgun metagenomic sequencing and LC-MS-based metabolomics. Despite these limitations, this study enhances our understanding of the effects of monensin and confirms the usefulness of metagenomics and metabolomics analyses in ruminant nutrition studies.

5. Conclusions

This study showed that feeding monensin altered functional and metabolomic attributes of the rumen microbiota in forage-fed beef cattle. Differences in relative abundance of functional genes involved in lipid metabolism and amino acid metabolism as well as changes in rumen fluid metabolites and their pathways revealed that dietary monensin affects rumen fermentation of forage-fed beef cattle by altering the rumen volatile fatty acid profiles in favor of propionate production, by reducing amino acid degradation, and by manipulating biohydrogenation of unsaturated fatty acid in the rumen.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2076-2615/8/11/211/s1, Table S1: Nutritional composition of the diet, Table S2: Average relative abundance of taxa at the genus level, Table S3: Average relative abundance of taxa at the genus level, Table S4: Average relative abundance of taxa at the genus level, Table S5: List of identified metabolites, Figure S1: Rarefaction curve of the 8 samples, Figure S2. The scores scatter plot of the PCA model showing the clustering of the quality control samples.

Author Contributions

Conceptualization, I.O. and K.A.; Formal analysis, J.A.; Investigation, I.O., H.S., J.L. and J.A.; Methodology, H.S.; Supervision, I.O. and J.L.; Writing—original draft, I.O.; Writing—review & editing, K.A., J.L. and J.A.

Funding

This study was supported by funds from the US Department of Agriculture’s National Institute of Food and Agriculture Evans-Allen project 1008985.

Acknowledgments

The authors are grateful to Christina Mcmanus for reviewing and editing this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Morgavi, D.P.; Kelly, W.J.; Janssen, P.H.; Attwood, G.T. Rumen microbial metagenomics and its application to ruminant production. Animal 2013, 7, 184–201. [Google Scholar] [CrossRef] [PubMed]
  2. Duffield, T.F.; Bagg, R.N. Use of ionophores in lactating dairy cattle: A review. Can. Vet. J. 2000, 41, 388–394. [Google Scholar] [PubMed]
  3. Haney, M.E.; Knox, N.G.; Hoehn, M.M. Monensin, a new biologically active compound. Antimicrob. Agents Chemother. 1967, 7, 353–358. [Google Scholar]
  4. Weimer, P.J.; Thomas, E.E.; Mertens, D.R.; Stevenson, D.M. Effect of monensin feeding and withdrawal on populations of individual bacterial species in the rumen of lactating dairy cows fed high-starch rations. Appl. Microbiol. Biotechnol. 2008, 80, 135–145. [Google Scholar] [CrossRef] [PubMed]
  5. Guan, H.; Wittenberg, K.M.; Ominski, K.H.; Krause, D.O. Efficacy of ionophores in cattle diets for mitigation of enteric methane. J. Anim. Sci. 2006, 84, 1896–1906. [Google Scholar] [CrossRef] [PubMed]
  6. Baird, D.G.; Lomax, M.A.; Symonds, H.W.; Shaw, S.R. Net hepatic and splanchnic metabolism of lactate, pyruvate and propionate in dairy cows in vivo in relation to lactation and nutrient supply. Biochem. J. 1980, 186, 47–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Aschenbach, J.R.; Penner, G.B.; Stump, F.; Gäbel, G. Role of fermentation acid absorption in the regulation of ruminal pH. J. Anim. Sci. 2011, 89, 1092–1107. [Google Scholar] [CrossRef] [PubMed]
  8. Russell, J.B.; Houlihan, A.J. Ionophore resistance of ruminal bacteria and its potential impact on human health. FEMS Microbiol. Rev. 2003, 27, 65–74. [Google Scholar] [CrossRef] [Green Version]
  9. Ercolini, D. High-throughput sequencing and metagenomics: Moving forward in the culture-independent analysis of food microbial ecology. Appl. Environ. Microbiol. 2013, 79, 3148–3155. [Google Scholar] [CrossRef] [PubMed]
  10. Jovel, J.; Patterson, J.; Wang, W.; Hotte, S.; O’Keefe, N.; Mitchel, T.; Perry, T.; Kao, D.; Mason, A.L.; Madsen, K.L.; et al. Characterization of the gut microbiome using 16S or shotgun metagenomics. Front. Microbiol. 2016, 7, 459. [Google Scholar] [CrossRef] [PubMed]
  11. Kim, M.; Eastridge, M.; Yu, Z. Investigation of ruminal bacterial diversity in dairy cattle fed supplementary monensin alone and in combination with fat, using pyrosequencing analysis. Can. J. Microbiol. 2013, 60, 65–71. [Google Scholar] [CrossRef] [PubMed]
  12. Aird, D.; Ross, M.G.; Chen, W.S.; Danielsson, M.; Fennell, T.; Russ, C.; Jaffe, D.; Nusbaum, C.; Gnirkel, A. Analyzing and minimizing PCR amplification bias in Illumina sequencing libraries. Genome Biol. 2011, 12, R18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ranjan, R.; Rani, A.; Metwally, A.; McGee, H.S.; Perkins, D.L. Analysis of the microbiome: Advantages of whole genome shotgun versus 16S amplicon sequencing. Biochem. Biophys. Res. Commun. 2016, 469, 967–977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Yang, X.; Noyes, N.R.; Doster, E.; Martin, J.N.; Linke, L.M.; Magnuson, R.J. Use of metagenomic shotgun sequencing technology to detect foodborne pathogens within the microbiome of the beef production chain. Appl. Environ. Microbiol. 2016, 82, 2433–2443. [Google Scholar] [CrossRef] [PubMed]
  15. Zinicola, M.; Higgins, H.; Lima, S.; Machado, V.; Guard, C.; Bicalho, R. Shotgun metagenomic sequencing reveals functional genes and microbiome associated with bovine digital dermatitis. PLoS ONE 2015, 10, e0133674. [Google Scholar] [CrossRef] [PubMed]
  16. Thomas, M.; Webb, M.; Ghimire, S.; Blair, A.; Olson, K.; Fenske, G.J.; Fonder, A.T.; Christopher-Hennings, J.; Brake, D.; Scaria, J. Metagenomic characterization of the effect of feed additives on the gut microbiome and antibiotic resistome of feedlot cattle. Sci. Rep. 2017, 7, 12257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Goldansaz, S.A.; Guo, A.C.; Sajed, T.; Steele, M.A.; Plastow, G.S.; Wishart, D.S. Livestock metabolomics and the livestock metabolome: A systematic review. PLoS ONE 2017, 12. [Google Scholar] [CrossRef] [PubMed]
  18. Karisa, B.K.; Thomson, J.; Wang, Z.; Li, C.; Montanholi, Y.R.; Miller, S.P.; Moore, S.S.; Plastow, G.S. Plasma metabolites associated with residual feed intake and other productivity performance traits in beef cattle. Livest. Sci. 2014, 165, 200–211. [Google Scholar] [CrossRef]
  19. Artegoitia, V.M.; Foote, A.P.; Lewis, R.M.; Freetly, H.C. Rumen fluid metabolomics analysis associated with feed efficiency on crossbred steers. Sci. Rep. 2017, 7, 2864. [Google Scholar] [CrossRef] [PubMed]
  20. Hailemariam, D.; Mandal, R.; Saleem, F.; Dunn, S.M.; Wishart, D.S.; Ametaj, B.N. Identification of predictive biomarkers of disease state in transition dairy cows. J. Dairy Sci. 2014, 97, 2680–2693. [Google Scholar] [CrossRef] [PubMed]
  21. Saleem, F.; Ametaj, B.N.; Bouatra, S.; Mandal, R.; Zebeli, Q.; Dunn, S.M.; Wishart, D.S. A metabolomics approach to uncover the effect of grain diets on rumen health in dairy cows. J. Dairy Sci. 2012, 95, 6606–6623. [Google Scholar] [CrossRef] [PubMed]
  22. Abarghuei, M.; Rouzbehan, J.Y.; Salem, A.Z.M.; Zamiri, M.J. Nitrogen balance, blood metabolites and milk fatty acid composition of dairy cows fed pomegranate-peel extract. Livest. Sci. 2014, 164, 72–80. [Google Scholar] [CrossRef]
  23. Qin, J.; Li, R.; Raes, J.; Arumugam, M.; Burgdorf, K.S.; Manichanh, C.; Nielsen, T.; Pons, N.; Levenez, F.; Yamada, T.; et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature 2010, 464, 59–65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Peng, Y.; Leung, H.C.M.; Yiu, S.M.; Chin, F.Y.L. IDBA-UD: A de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics 2012, 28, 1420–1428. [Google Scholar] [CrossRef] [PubMed]
  25. Noguchi, H.; Taniguchi, T.; Itoh, T. MetaGene Annotator: Detecting species-specific patterns of ribosomal binding site for precise gene prediction in anonymous prokaryotic and phage genomes. DNA Res. 2008, 15, 387–396. [Google Scholar] [CrossRef] [PubMed]
  26. Li, W.; Godzik, A. CD-HIT: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [PubMed]
  27. Wilke, A.; Harrison, T.; Wilkening, J.; Field, D.; Glass, E.M.; Kyrpides, N.; Mavrommatis, K.; Meyer, F. The M5nr: A novel non-redundant database containing protein sequences and annotations from multiple sources and associated tools. BMC Bioinform. 2012, 13, 141. [Google Scholar] [CrossRef] [PubMed]
  28. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, 60. [Google Scholar] [CrossRef] [PubMed]
  29. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2002, 8, 27–30. [Google Scholar] [CrossRef]
  30. Lombard, V.; Golaconda, R.H.; Drula, E.; Coutinho, P.M.; Henrissat, B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014, 42, 490–495. [Google Scholar] [CrossRef] [PubMed]
  31. Want, E.J.; Masson, P.; Michopoulos, F.; Wilson, I.D.; Theodoridis, G.; Plumb, R.S. Global metabolic profiling of animal and human tissues via UPLC-MS. Nat. Protoc. 2013, 8, 17–32. [Google Scholar] [CrossRef] [PubMed]
  32. Muck, R.E.; Dickerson, J.T. Storage temperature effects on proteolysis in alfalfa silage. Trans. ASAE 1988, 31, 1005–1009. [Google Scholar] [CrossRef]
  33. Tsugawa, H.; Cajka, T.; Kind, T.; Ma, Y.; Higgins, B.; Ikeda, K.; Mitsuhiro, K.; VanderGheynst, J.; Fiehn, O.; Arita, M. MS-DIAL: Data independent MS/MS deconvolution for comprehensive metabolome analysis. Nat. Methods 2015, 12, 523–526. [Google Scholar] [CrossRef] [PubMed]
  34. Xia, J.; Wishart, D.S. Using metaboAnalyst 3.0 for comprehensive metabolomics data analysis. Curr. Protoc. Bioinform. 2016, 55, 14. [Google Scholar] [CrossRef] [PubMed]
  35. Bylesjo, M.; Rantalainen, M.; Cloarec, O.; Nicholson, J.K.; Holmes, E.; Trygg, J. OPLS discriminant analysis: Combining the strengths of PLS-DA and SIMCA classification. J. Chemom. 2006, 20, 341–351. [Google Scholar] [CrossRef]
  36. Guarner, F. Enteric flora in health and disease. Digestion 2006, 73, 5–12. [Google Scholar] [CrossRef] [PubMed]
  37. Malmuthuge, N.; Guan, L.L. Understanding host-microbial interactions in rumen: Searching the best opportunity for microbiota manipulation. Microbiome 2017, 8, 8. [Google Scholar] [CrossRef] [PubMed]
  38. Weimer, P.J. Redundancy, resilience, and host specificity of the ruminal microbiota: Implications for engineering improved ruminal fermentations. Front. Microbiol. 2015, 6, 296. [Google Scholar] [CrossRef] [PubMed]
  39. Fernando, S.C.; Purvis, H.T.; Najar, F.Z.; Sukharnikov, L.O.; Krehbiel, C.R.; Nagaraja, T.G.; Roe, B.A.; DeSilva, U. Rumen microbial population dynamics during adaptation to a high-grain diet. Appl. Environ. Microbiol. 2010, 76, 7482–7490. [Google Scholar] [CrossRef] [PubMed]
  40. Iqbal, M.W.; Zhang, Q.; Yang, Y.; Li, L.; Zou, C.; Huang, C.; Lin, B. Comparative study of rumen fermentation and microbial community differences between water buffalo and Jersey cows under similar feeding conditions. J. Appl. Anim. Res. 2017, 46, 740–748. [Google Scholar] [CrossRef] [Green Version]
  41. Macy, J.M.; Ljungdahl, L.G.; Gottschalk, G. Pathway of succinate and propionate formation in Bacteroides fragilis. J. Bacteriol. 1978, 134, 84–91. [Google Scholar] [PubMed]
  42. Moore, L.V.; Moore, W.E.C. Oribaculum catoniae gen. nov., sp. nov., Catonella morbi gen. nov., sp. nov., Hallella seregens gen. nov., sp. nov., Johnsonella ignava gen. nov., sp. nov., and Dialister pneumosintes gen. nov., comb. nov., nom. rev., anaerobic gram-negative bacilli from the human gingival crevice. Int. J. Syst. Bacteriol. 1994, 44, 187–192. [Google Scholar] [PubMed]
  43. Sakamoto, M.; Suzuki, N.; Matsunaga, N.; Koshihara, K.; Seki, M.; Komiya, H.; Benno, Y. Parabacteroides gordonii sp. Nov., isolated from human blood cultures. Int. J. Syst. Evol. Microbiol. 2009, 59, 2843–2847. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, M.; Wolin, M.J. Effect of monensin and lasalocid-sodium on the growth of methanogenic and rumen saccharolytic bacteria. Appl. Environ. Microbiol. 1979, 38, 72–77. [Google Scholar] [PubMed]
  45. Ueki, A.; Watanabe, M.; Ohtaki, Y.; Kaku, N.; Ueki, K. Description of Propionispira arcuata sp. nov., isolated from a methanogenic reactor of cattle waste, reclassification of Zymophilus raffinosivorans and Zymophilus paucivorans as Propionispira raffinosivorans comb. nov. and Propionispira paucivorans comb. nov. and emended description of the genus Propionispira. Int. J. Syst. Evol. Microbiol. 2014, 64, 3571–3577. [Google Scholar] [PubMed]
  46. Bell, N.L.; Anderson, R.C.; Callaway, T.R.; Franco, M.O.; Sawyer, J.E.; Wickersham, T.E. Effect of monensin inclusion on intake, digestion, and ruminal fermentation parameters in Bos indicus and Bos taurus steers consuming bermudagrass hay. J. Anim. Sci. 2017, 95, 2736–2746. [Google Scholar] [CrossRef] [PubMed]
  47. Lemenager, R.P.; Owens, F.N.; Shockey, B.J.; Lusby, K.S.; Totusek, R. Monensin, forage intake and lactation of range beef cows. J. Anim. Sci. 1978, 47, 247–254. [Google Scholar] [CrossRef]
  48. Reynolds, C.K.; Aikman, P.C.; Lupoli, B.; Humphries, D.J.; Beever, D.E. Splanchnic metabolism of dairy cows during the transition from late gestation through early lactation. J. Dairy Sci. 2003, 86, 1201–1217. [Google Scholar] [CrossRef]
  49. Bergen, W.G.; Bates, D.B. Ionophores: Their effect on production efficiency and mode of action. J. Anim. Sci. 1984, 58, 1465–1483. [Google Scholar] [CrossRef] [PubMed]
  50. Patra, A.; Park, T.; Kim, M.; Yu, Z. Rumen methanogens and mitigation of methane emission by anti-methanogenic compounds and substances. J. Anim. Sci. Biotechnol. 2017, 8, 13. [Google Scholar] [CrossRef] [PubMed]
  51. Sun, H.Z.; Wang, D.M.; Wang, B.; Wang, J.K.; Liu, H.Y.; le Guan, L. Metabolomics of four biofluids from dairy cows: Potential biomarkers for milk production and quality. J. Proteome Res. 2015, 14, 1287–1298. [Google Scholar] [CrossRef] [PubMed]
  52. Fellner, V.; Sauer, F.D.; Kramer, J.K.G. Effect of nigericin, monensin, and tetronasin on biohydrogenation in continuous flow-through ruminal fermenters. J. Dairy Sci. 1997, 80, 921–928. [Google Scholar] [CrossRef]
  53. Van Nevel, C.J.; Demeyer, D.I. Lipolysis and biohydrogenation of soybean oil in the rumen in vitro: Inhibition by antimicrobials. J. Dairy Sci. 1995, 78, 2797–2806. [Google Scholar] [CrossRef]
  54. Jeyanathan, J.; Escobar, M.; Wallace, R.J.; Fievez, V.; Vlaeminck, B. Biohydrogenation of 22:6n-3 by Butyrivibrio proteoclasticus P18. BMC Microbiol. 2016, 16, 104. [Google Scholar] [CrossRef] [PubMed]
  55. Wang, J.H.; Choi, S.H.; Lim, K.W.; Kim, K.H.; Song, M.K. Effect of the mixed oil and monensin supplementation, and feeding duration of supplements on c9, t11-CLA contents in plasma and fat tissues of Korean native (Hanwoo) steers. Asian-Aust. J. Anim. Sci. 2006, 10, 1464–1469. [Google Scholar] [CrossRef]
  56. Sauer, F.D.; Fellner, V.; Kinsman, R.; Kramer, J.K.G.; Jackson, H.A.; Lee, A.J.; Chen, S. Methane output and lactation response in Holstein cattle with monensin or unsaturated fat added to the diet. J. Anim. Sci. 1998, 76, 906–914. [Google Scholar] [CrossRef] [PubMed]
  57. Da Silva-Kazama, D.C.; Santos, G.T.; Branco, A.F.; Damasceno, J.C.; Kazama, R.; Matsushita, M.; Horst, J.A.; dos Santos, W.B.R.; Petit, H.V. Production performance and milk composition of dairy cows fed whole or ground flaxseed with or without monensin. J. Dairy Sci. 2007, 90, 2928–2936. [Google Scholar] [CrossRef] [PubMed]
  58. Odongo, N.E.; Bagg, R.; Vessie, G.; Dick, P.; Or-Rashid, M.M.; Hook, S.E.; Gray, J.T.; Kebreab, E.; France, J.; MCbride, B.W. Long-term effects of feeding monensin on methane production in lactating dairy cows. J. Dairy Sci. 2007, 90, 1781–1788. [Google Scholar] [CrossRef] [PubMed]
  59. Da Silva-Kazama, D.C.; Côrtes, C.; Kazama, R.; Gagnon, N.; Benchaar, C.; Santos, G.T.; Zeoula, L.M.; Petit, H.V. Ruminal fermentation characteristics and fatty acid profile of ruminal fluid and milk of dairy cows fed flaxseed hulls supplemented with monensin. J. Dairy Res. 2011, 78, 56–62. [Google Scholar] [CrossRef] [PubMed]
  60. Ladeira, M.M.; Santarosa, L.C.; Chizzotti, M.L.; Ramos, E.M.; Machado Neto, O.R.; Oliveira, D.L. Fatty acid profile, color and lipid oxidation of meat from young bulls fed ground soybean or rumen protected fat with or without monensin. Meat Sci. 2014, 96, 597–605. [Google Scholar] [CrossRef] [PubMed]
  61. Bhattacharya, A.; Banu, J.; Rahman, M.; Causey, J.; Fernandes, G. Biological effects of conjugated linoleic acids in health and disease. J. Nutr. Biochem. 2006, 17, 789–810. [Google Scholar] [CrossRef] [PubMed]
  62. Eschenlauer, S.C.P.; McKain, N.; Walker, N.D.; McEwan, N.R.; Newbold, C.J.; Wallace, R.J. Ammonia production by rumen microorganisms and enumeration, isolation, and characterization of bacteria capable of growth on peptides and amino acids from the sheep rumen. Appl. Environ. Microbiol. 2002, 68, 4925–4931. [Google Scholar] [CrossRef] [PubMed]
  63. Benchaar, C.; Petit, H.V.; Berthiaume, R.; Whyte, T.D.; Chouinard, P.Y. Effects of addition of essential oils and monensin premix on digestion, ruminal fermentation, milk production, and milk composition in dairy cows. J. Dairy Sci. 2006, 89, 4352–4364. [Google Scholar] [CrossRef]
  64. Russell, J.; Bottje, W.; Cotta, M. Degradation of protein by mixed cultures of rumen bacteria: Identification of Streptococcus bovis as an actively proteolytic rumen bacterium. J. Anim. Sci. 1981, 53, 242–252. [Google Scholar] [CrossRef] [PubMed]
  65. Russell, J.B.; Onodera, R.; Hino, T. Ruminal protein fermentation: New perspectives on previous contradictions. In Physiological Aspects of Digestion and Metabolism in Ruminants; Tsuda, T., Sasaki, Y., Kawashima, R., Eds.; Academic Press: Tokyo, Japan, 1991; pp. 681–697. [Google Scholar]
  66. Yang, C.M.J.; Russell, J.B. Effect of monensin on the specific activity of ammonia production by ruminal bacteria and disappearance of amino nitrogen from the rumen. Appl. Environ. Microbiol. 1993, 59, 3250–3254. [Google Scholar] [PubMed]
  67. Callaway, T.R.; Carneiro De Melo, A.M.S.; Russell, J.B. The effects of monensin on ruminal fermentations in vitro. Curr. Microbiol. 1997, 35, 90–96. [Google Scholar] [CrossRef] [PubMed]
  68. Gevers, D.; Pop, M.; Schloss, P.D.; Huttenhower, C. Bioinformatics for the Human Microbiome Project. PLoS Comput. Biol. 2012, 8, e1002779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. McAllister, T.A.; Duniere, L.; Drouin, P.; Xu, S.; Wang, Y.; Munns, K.; Zaheer, R. Silage review: Using molecular approaches to define the microbial ecology of silage. J. Dairy Sci. 2018, 101, 4060–4074. [Google Scholar] [CrossRef] [PubMed]
  70. Veenstra, T.D. Metabolomics: The final frontier? Genome Med. 2012, 4, 40. [Google Scholar] [CrossRef] [PubMed]
  71. Aretz, I.; Meierhofer, D. Advantages and Pitfalls of mass spectrometry based metabolome profiling in systems biology. Int. J. Mol. Sci. 2016, 17, 632. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Linear discriminant analysis effect size (LEfSe) of rumen microbiota of beef steer fed no (control) or 200 mg/d of monensin. The linear discriminant analysis plot indicates the most differentially abundant taxa found by ranking according to their effect size (≥2.0) at the genus (A) and species level (B). The taxa enriched in steers fed the control diet are indicated with a positive score (green), and taxa enriched by the monensin treatment have a negative score (red). Only taxa meeting the significant threshold of 2.0 are shown.
Figure 1. Linear discriminant analysis effect size (LEfSe) of rumen microbiota of beef steer fed no (control) or 200 mg/d of monensin. The linear discriminant analysis plot indicates the most differentially abundant taxa found by ranking according to their effect size (≥2.0) at the genus (A) and species level (B). The taxa enriched in steers fed the control diet are indicated with a positive score (green), and taxa enriched by the monensin treatment have a negative score (red). Only taxa meeting the significant threshold of 2.0 are shown.
Animals 08 00211 g001
Figure 2. (A) Distribution of category of the predicted genes by KEGG. (B) Differential KEGG gene functions. Differences between control and monensin samples were tested for significance using a Mann Whitney test (p ≤ 0.05).
Figure 2. (A) Distribution of category of the predicted genes by KEGG. (B) Differential KEGG gene functions. Differences between control and monensin samples were tested for significance using a Mann Whitney test (p ≤ 0.05).
Animals 08 00211 g002
Figure 3. Relative abundance of the category of carbohydrate-active enzymes according to carbohydrate-active enzyme (CAZy) database. AA = Auxilliary Activities, CBM = Carbohydrate Binding Modules, CE = Carbohydrate Esterases, GH = Glycoside Hydrolases, GT = Glycosyl Transferases, PL = Polysaccharide Lyases.
Figure 3. Relative abundance of the category of carbohydrate-active enzymes according to carbohydrate-active enzyme (CAZy) database. AA = Auxilliary Activities, CBM = Carbohydrate Binding Modules, CE = Carbohydrate Esterases, GH = Glycoside Hydrolases, GT = Glycosyl Transferases, PL = Polysaccharide Lyases.
Animals 08 00211 g003
Figure 4. (A) The scores plot of PCA model showing the directions that best explain the variance between the two treatments. (B) OPLS-DA score plot of all metabolite features. Group 1 = steers fed Control diet, Group 2 = steers fed 200 mg d−1 of monensin. One data point represents one composite rumen fluid sample of each steer.
Figure 4. (A) The scores plot of PCA model showing the directions that best explain the variance between the two treatments. (B) OPLS-DA score plot of all metabolite features. Group 1 = steers fed Control diet, Group 2 = steers fed 200 mg d−1 of monensin. One data point represents one composite rumen fluid sample of each steer.
Animals 08 00211 g004
Table 1. Fold changes of differential rumen fluid metabolites in beef steers fed no or 200 mg d−1 of monensin.
Table 1. Fold changes of differential rumen fluid metabolites in beef steers fed no or 200 mg d−1 of monensin.
MetaboliteFCp-Value
Acetate (mM)0.870.05
4-aminophenol0.800.09
Trans-cinnamaldehyde1.310.07
Hypoxanthine1.200.08
Tyramine0.900.07
Formylindole1.370.09
Indole-3-carboxylic acid1.440.05
l-β-homomethionine1.750.03
3-(4,4-dimethyl-4,5-dihydro-1,3-oxazol-2-yl)pyridine0.810.05
4-pyridoxic acid0.800.07
Simazine1.230.03
1,3-dipropyl-7-methylxanthine0.650.02
2-(1,3-benzodioxol-5-yl)ethyl-4-methoxy-2H-pyranone1.450.09
Pentalenolactone0.650.02
l-histidine1.760.04
Vildagliptin0.700.06
Oxybuprocaine1.410.05
Tri-O-methylgenistein1.480.01
Isopongaflavone1.560.07
Daphnoretine acetate1.480.02
l-phenylalanine1.960.01
Oxysporidinone1.240.02
Catechol1.420.03
2-methylgluric acid1.130.08
Oxypurinol1.230.02
protocatechuic acid1.390.01
Pimelic acid1.380.01
Isobutylglutaric acid1.130.05
Apigenin0.800.06
Linoleic acid1.610.03
Eriodictyol 7,3′-dimethyl ether1.170.09
8,15-DiHETE1.480.02
FC = fold change, mean value of peak intensity obtained from monensin group ÷ mean value of peak intensity obtained from control group. FC values >1 means that metabolite is greater in steers fed monensin and FC values < 1 means that metabolite is lower in steers fed monensin.
Table 2. Metabolic pathway analysis of the differential metabolites.
Table 2. Metabolic pathway analysis of the differential metabolites.
AnnotationImpactp-ValueTCHit
Caffeine metabolism00.15121
Purine metabolism0.010.62681
Tyrosine metabolism0.030.44421
Aminoacyl-tRNA biosynthesis00.22642
Vitamin B6 metabolism00.1291
Phenylalanine, tyrosine and tryptophan biosynthesis0.50.0541
Phenylalanine metabolism0.410.1291
Linoleic acid metabolism1.000.0651
Biosynthesis of unsaturated fatty acids00.41421
Histidine metabolism0.270.17141
Nitrogen metabolism00.1191
Glycolysis or gluconeogenesis0.030.30261
Pyruvate metabolism0.060.26221
Propanoate metabolism00.24201
Impact = pathway impact value calculated from pathway topology analysis; TC = total number of compounds in the pathway; Hit = matched number from the uploaded metabolite data.

Share and Cite

MDPI and ACS Style

Ogunade, I.; Schweickart, H.; Andries, K.; Lay, J.; Adeyemi, J. Monensin Alters the Functional and Metabolomic Profile of Rumen Microbiota in Beef Cattle. Animals 2018, 8, 211. https://0-doi-org.brum.beds.ac.uk/10.3390/ani8110211

AMA Style

Ogunade I, Schweickart H, Andries K, Lay J, Adeyemi J. Monensin Alters the Functional and Metabolomic Profile of Rumen Microbiota in Beef Cattle. Animals. 2018; 8(11):211. https://0-doi-org.brum.beds.ac.uk/10.3390/ani8110211

Chicago/Turabian Style

Ogunade, Ibukun, Hank Schweickart, Kenneth Andries, Jerusha Lay, and James Adeyemi. 2018. "Monensin Alters the Functional and Metabolomic Profile of Rumen Microbiota in Beef Cattle" Animals 8, no. 11: 211. https://0-doi-org.brum.beds.ac.uk/10.3390/ani8110211

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop