Effect of Dietary Protein Levels on Dynamic Changes and Interactions of Ruminal Microbiota and Metabolites in Yaks on the Qinghai-Tibetan Plateau

To improve performance and optimize rumen function in yaks (Bos grunniens), further knowledge on the appropriate dietary protein levels for ruminal microbiota and the metabolite profiles of yaks in feedlot feeding is necessary. Current understanding of dietary protein requirements, ruminal microbiota, and metabolites is limited. In this study, yaks were fed a low-protein diet (L; 9.64%), middle low-protein diet (ML; 11.25%), middle high-protein diet (MH; 12.48%), or a high-protein diet (H; 13.87%), and the effects of those diets on changes and interactions in ruminal microbiota and metabolites were investigated. Twenty-four female yaks were selected, and the effects on ruminal microbiota and metabolites were investigated using 16s rRNA gene sequencing and gas chromatography time-of-flight/mass spectrometry (GC-TOF/MS). Diets containing different protein levels changed the composition of the rumen bacterial community, the H group significantly reduced the diversity of ruminal microbiota (p < 0.05), and the number of shared amplicon sequence variants (ASVs) between the H group and the other three groups was lower, suggesting that the ruminal microbiota community fluctuated more with a high-protein diet. In rumen, Bacteroidetes, Firmicutes, and Proteobacteria were the most abundant bacteria at the phylum level, and Bacteroidetes was significantly less abundant in the MH group than in the L and ML groups (p < 0.05). Prevotella_1, Rikenellaceae_RC9_gut_group, and Christensenellaceae_R-7_group had the highest abundance at the genus level. Prevotellaceae was enriched in the low-protein groups, whereas Bacteroidales_BS11_gut_group was enriched in the high-protein groups. Rumen metabolite concentrations and metabolic patterns were altered by dietary protein levels: organic acid metabolites, antioxidant-related metabolites, and some plant-derived metabolites showed variation between the groups. Enrichment analysis revealed that significant changes were concentrated in six pathways, including the citrate cycle (TCA cycle), glyoxylate and dicarboxylate metabolism, and butanoate metabolism. Network analysis showed promotion or restraint relationships between different rumen microbiota and metabolites. Overall, the rumen function was higher in the MH group. This study provides a reference for appropriate dietary protein levels and improves understanding of rumen microbes and metabolites.

Introduction

The yak (Bos grunniens) is a unique species that lives on the Qinghai-Tibetan Plateau (QTP). It is adapted to high-altitude and low-oxygen environments and provides resources such as meat, milk, and fur for local pastoralists (Ding et al., 2014). Female yaks are important for the entire population; they provide enough milk for three meals for the pastoralists and are crucial to the reproduction of the herd. The health and nutritional status of females has a direct impact on young animals. As yaks rely on natural pasture, seasonal nutritional imbalances have a significant impact, especially in winter when pastures are insufficient for their maintenance requirements, causing weight loss and even death (Long et al., 2008). Pastoralists remain heavily impacted by traditional values that affect yak productivity and the ability to resist disasters (such as snowstorms). Therefore, some pastoralists are gradually changing their traditional grazing to indoor feeding or supplements in winter. This has had a positive effect on livestock growth performance, grazing stress, and pastoralist income (Liu P. et al., 2019). However, the lack of nutritional requirements of yak is another major knowledge gap that hinders development of the animal husbandry. Resources for diet and the nutrition (especially the protein and energy content) are extremely important. Ministry of Agriculture of the PRC (2004) was launched due to market demand, allowing the centralized management of beef cattle, but the standard of yaks still needs to be explored.

Rumen function plays a critical role in ruminant production, animal health, and metabolism. The rumen microbiome consists of bacteria, fungi, archaea, and protozoa, which are required to degrade and transform feedstuffs, fermented plant proteins, and polysaccharides. These provide metabolites such as volatile fatty acids (VFAs), amino acids, and saccharides used by organisms to promote growth and enable microbial reproduction (Shabat et al., 2016). This process is equally important for humans because it harnesses the solar energy stored in plant fibers and converts it into milk and meat (Jami et al., 2013). A number of studies have shown that diet is a key factor influencing rumen microbes (Kim et al., 2011; Yeoman et al., 2018). Microbiome studies are becoming increasingly easier, and big data-driven omics research has led to a steep rise in investigations involving two of the most functional omes, the microbiome and metabolome (Misra, 2020). The former includes both the microbiota and its genome, while the latter is a useful method to identify the health status of an organism in nutrition regulation studies, and is widely used in nutritional studies and more frequently used in ruminant studies (Zeisel et al., 2005; Drackley et al., 2006). Therefore, a combination of the microbiota and metabolites could offer important insights. A study of both rumen microbes and metabolites in yak found that the microbiota composition was significantly different in the forge group compared to the concentrate group. This affected the concentration of metabolites and metabolic way, as well as the concentrations of ruminal metabolites participating in protein digestion and absorption (Liu C. et al., 2019). Another study on dairy cow metabolomics found that a high-concentrate diet (70%) increased the concentration of bacterial degradation products, some toxic compounds and 15 amino acids, which changed the rumen metabolic pattern, with more complicated metabolites compared to the low concentrate diet (40%) (Zhang et al., 2017). Although microbiome and metabolome studies have been widely applied in animal nutrition, there is a lack of information on yaks. Thus, improving research on yaks will help us to further understand the rumen environment.

16S rRNA sequencing is an effective method analyzing microbiota with techniques improving steadily with the development of microbiology. Second-generation technology is mature and available to researchers to achieve the research objectives efficiently. Metabolomic analysis methods include nuclear magnetic resonance (NMR), liquid chromatography–mass spectrometry (LC-MS), and gas chromatography–mass spectrometry (GC-MS). Techniques can be selected based on the physical and chemical properties of the sample (Krone et al., 2010; Drexler et al., 2011; Oakman et al., 2011). In recent years, LC-MS and GC-MS techniques have been used extensively. LC-MS is most suited to investigating higher-molecular-weight compounds, but GC-MS is widely used because of its high-throughput capability and high sensitivity (Styczynski et al., 2007; Liu et al., 2010).

In this study, 16S rRNA sequencing and gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF/MS) were used to explore the effect of dietary protein levels on ruminal microbiota and metabolite profiles in yaks. An association analysis was conducted to identify the relationship between ruminal microbiota and metabolites. It is hypothesized that different protein levels impact the ruminal microbiota and metabolites.

Methods

Animals, Diets, and Experimental Design

This study was conducted in January 2019 at the Haibei Demonstration Zone of Plateau Modern Ecological Animal Husbandry Science and Technology, in Qinghai Province, China.

Twenty-four two-year-old healthy female yaks (body weight: 107.54 ± 4.72 kg) were selected from the grazing grassland. Each yak was ear-marked and randomly assigned into four groups. Yaks were fed a low-protein diet (L; 9.64%), middle low-protein diet (ML; 11.25%), middle high-protein diet (MH; 12.48%), or a high-protein diet (H; 13.87%). All yaks were initially fed 1% of their mean body weight. The amount was increased daily until day 15 when the volume of food reached 1.5% of the yak’s body weight. After day 15, the amount of feed was adjusted every 2 weeks based on body weight. The diet was formulated according to the nutrient requirements for dairy heifers (Council, 2007), with varying protein levels and the same metabolizable energy (ME) level. The nutrient composition is provided in Table 1. Yaks were allowed free access to water and fed twice daily at 8:00 and 17:00.

TABLE 1

Table 1. Ingredients and nutrient composition of each diet.

The crude protein (CP), neutral detergent fiber (NDF), acid detergent fiber (ADF), Ca, and P in each sample were measured in the laboratory, and the ME was calculated. Mixed feed (100 g) was collected and dried in a forced-air oven at 60°C for 48 h, then ground through a 1-mm sieve before analysis. DM and N contents were determined according to AOAC (2007). NDF and ADF contents were analyzed according to Vansoest et al. (1991).

Rumen Sample Collection and Measurements

At the end of the trial (day 135), rumen fluid samples were collected before the first morning feeding using an oral stomach tube. Rumen samples were strained through four layers of cheesecloth to obtain rumen liquids. Disposable sterile gloves were worn throughout the sample collection process to avoid contamination. Twenty-four samples were collected and placed in 2-ml frozen tubes to avoid cross-contamination. The samples were immediately frozen in liquid nitrogen and stored at −80°C.

Microbial DNA Extraction, PCR Amplification, Sequencing, and Sequencing Data Processing

Total genome DNA was extracted using the CTAB/SDS method. DNA concentration and purity was monitored on 1% agarose gels. According to the concentration, DNA was diluted to 1 ng/μl using sterile water. The 16S rRNA genes of distinct regions (16S V3–V4) were amplified using specific primers 341F (5′-CCTAYGGGRBGCASCAG-3′) and 806R (5′-GGACTACNNGGGTATCTAAT-3′) with barcode. All PCR reactions were carried out using 15 μl of Phusion^® High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μM forward and reverse primers, and approximately 10 ng of template DNA. Thermal cycling consisted of initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, elongation at 72°C for 30 s, followed by 72°C for 5 min. The same volume of 1 × loading buffer (containing SYB green) was mixed with the PCR products and electrophoresis was operated on 2% agarose gel for detection. PCR products were mixed in equidensity ratios and then purified with Qiagen Gel Extraction Kit (Qiagen, Germany). Sequencing libraries were generated using TruSeq^® DNA PCR-Free Sample Preparation Kit (Illumina, United States) following the manufacturer’s recommendations, and index codes were added. The library quality was assessed on the Qubit@2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. The library was sequenced on an Illumina NovaSeq platform and 250-bp paired-end reads were generated by Novogene Bioinformatics Technology Co., Ltd. (Beijing, China).

The 16S rRNA gene sequences were processed using USEARCH v10.0 (Edgar, 2010) and scripts were written by Liu et al. (2020). Metadata are given in Supplementary Table 1. The quality of the paired-end Illumina reads was checked by FastQC v.0.11.5 (Andrews, 2013) and processed in the following steps by USEARCH: merging paired reads and relabeling of sequencing names, removal of barcodes and primers, filtering of low-quality reads, and finding non-redundancy reads. Unique reads were denoised into amplicon sequence variants (ASVs) by unoise3 in USEARCH (Edgar and Flyvbjerg, 2015). The feature table (ASV table) was generated by VSEARCH (Rognes et al., 2016). The SILVA v123 (Quast et al., 2013) database was used to classified the taxonomy of the representative sequences, and then plastid and non-bacteria were removed.

The alpha diversity and Bray--Curtis distance-based constrained principal coordinate analysis (CPCoA) were accessed using USEARCH v10.0. The alpha diversity boxplot, CPCoA plot, taxonomy barplot, and UpSet plot visualization were performed with ImageGP website¹. The EdgeR package was used for bacteria gene differential analysis, and group-to-group differential analysis data are displayed as a Manhattan plot.

All the raw sequences after assembling and filtering were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive², under accession number PRJNA706154.

Preparation of Rumen Sample for GC-TOF/MS and Identification of Compounds

Rumen fluid samples were dried completely in a vacuum concentrator without heating, and then 60 μl of methoxyamination hydrochloride (20 mg/ml in pyridine) was added and incubated for 30 min at 80°C. BSTFA reagent (80 μl, 1% TMCS, v/v) was added to the sample aliquots and incubated for 1.5 h at 70°C. Samples were analyzed by an Agilent 7890 gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF/MS). The system utilized a DB-5MS capillary column coated with 5% diphenyl cross-linked with 95% dimethylpolysiloxane (30 m × 250 μm inner diameter, 0.25 μm film thickness; J&W Scientific, Folsom, CA, United States). A 1-μl aliquot of analyte was injected in splitless mode. Helium was used as the carrier gas, the front inlet purge flow was 3 ml min^–1, and the gas flow rate through the column was 1 ml min^–1. The initial temperature was kept at 50°C for 1 min, then raised to 310°C at a rate of 20°C min^–1, and then kept at 310°C for 6 min. The injection, transfer line, and ion source temperatures were 280°C, 280°C, and 250°C, respectively. The energy was −70 eV in electron impact mode. The mass spectrometry data were acquired in full-scan mode with m/z of 50–500, at a rate of 12.5 spectra per second after a solvent delay of 4.78 min.

Chroma TOF 4.3X software of LECO Corporation and the LECO-Fiehn Rtx5 database (Kind et al., 2009) were used for raw peaks exacting, data baselines filtering, and calibration of the baseline, peak alignment, deconvolution analysis, peak identification, and integration of the peak area. Both mass spectrum match and retention index match were considered in the metabolite identification.

The resulting data, with the compound name, sample label, and normalized peak area were imported into SIMCA-P software (V14.1, Umetrics AB, Umea, Sweden) for orthogonal projections to latent structures-discriminant analysis (OPLS-DA). Differentially expressed metabolites (DEMs) were identified, combing Variable Importance in Projection (VIP) obtained from OPLS-DA analysis and t-test (VIP > 1 and p < 0.01). DEMs were further identified and validated by the Bovine Metabolome Database (BMDB³) and the Kyoto Encyclopedia of Genes and Genomes (KEGG⁴). To identify the patterns of different metabolites with increasing dietary protein levels, the linear, quadratic, and cubic effects among treatments were evaluated among the different metabolites (Zhang et al., 2020). DEMs were imported into the MetaboAnalyst web server⁵ to view their metabolic pathway distribution and enrichment analysis (Xia et al., 2009).

Correlations Between Microbial Communities and Rumen Metabolites

The DEMs selected by the previous step and the microbes screened by Edge R were used for interactive analysis in R. Spearman correlations were calculated using the Psych packages. The metabolic correlation network was visualized using the Fruchterman Reingold Algorithm in Gephi 0.9.2 software⁶ (Bastian et al., 2009).

Statistical Analysis

Statistical analysis was performed using R (v3.6.1) and SPSS 17.0. The linear, quadratic, and cubic effects of treatments were evaluated by the lm function of Estimability package in R, according to Zhang et al. (2020).

Ethical Consideration

All procedures involving animal care adhered to the guidelines provided by the Institution of Animal Care and the Ethics Committee of the Northwest Institute of Plateau Biology, Chinese Academy of Sciences (NWIPB20160302).

Results

Variation in Bacterial Diversity in Different Treatment Groups

In total, 2,105,160 raw reads were obtained from the 24 rumen samples, after screening, and 1,500,085 high-quality sequences were obtained (average, 62,503; range, 40,833–71,120 reads per sample). A total of 5600 ASVs were produced based on denoising. The rarefaction curves and richness barplot for the ASVs showed that the quality of the observed species increased with sequencing depth, the H group had the lowest richness of all groups (Supplementary Figure 1). The ACE index was similar in the L, ML, and MH groups, but significantly lower in the H group (p < 0.05). The Shannon index had a similar tendency as ACE but was not significant (Figures 1A,B). Analysis of beta diversity with CPCoA (Bray–Curtis distance) showed that rumen microbiota formed four distinct clusters: the L and ML groups separated along the second coordinate axis, while the MH and H groups separated along both coordinate axis with other groups (p = 0.001) (Figure 1C). The CPCoA analysis accounted for 18.5% of variance.

FIGURE 1

Figure 1. Bacterial diversity analysis of rumen samples between yak fed diets with different protein levels. The alpha diversity between different groups (A) ACE index; (B) Shannon diversity; beta diversity (C) CPCoA plot using Bray–Curtis dissimilarity based on ASVs in different groups.

Variation in Bacterial Composition in Different Treatment Groups

Across all groups, 22 phyla, 44 classes, 73 orders, 102 families, and 183 genera were identified. The relative abundance of dominant taxa at the phylum and genus levels are shown in Figures 2A–C. At the phylum level, the most abundant phyla were Bacteroidetes (49.73%), which were significantly less abundance in the MH group than in the L and ML groups (p < 0.05), but the H group was not significantly different from the other groups. Firmicutes (41.02%) and Proteobacteria (2.63%) were enriched in the MH group, but were not significantly differentiated between the other groups (p > 0.05). The Firmicutes/Bacteroidetes (F/B) ratio in the MH group was higher than other groups and significantly higher than the L group (p < 0.05). Less abundant phyla included Spirochaetae (1.38%), Fibrobacteres (0.96%), Tenericutes (0.93%), and Cyanobacteria (0.79%), whereas Fibrobacteres and Tenericutes showed no significant differences between the groups. Spirochaetae were higher in the H group and significantly higher than in the ML group. Cyanobacteria was lower in the H group compared to other groups, and significantly lower than in the L group. At the genus level, Prevotella_1 (16.09%), Rikenellaceae_RC9_gut_group (8.95%), and Christensenellaceae_R-7_group (7.73%) were the predominant genera, Prevotella_1 was significantly enriched in the L and ML groups compared to the MH and H groups. Rikenellaceae_RC9_gut_group was significantly more abundant in the H group than in the ML group. Christensenellaceae_R-7_group was highest in the MH group, but not significantly. Other less abundant genera included Prevotellaceae_UCG-003 (4.43%), Ruminococcaceae_UCG-014 (2.82%), and Ruminococcaceae_UCG-010 (2.48%). Prevotellaceae_UCG-003 was significantly more abundant in the L and ML groups than in the MH group. However, Ruminococcaceae_UCG-014 showed no significant differences between the groups. Ruminococcaceae_UCG-010 was significantly higher in the H group compared to the L group (Figure 2D).

FIGURE 2

Figure 2. Bacterial composition of rumen samples collected from yaks fed diets with different protein levels. Bacterial composition at the (A) phylum and (B) genus level, (C) further details for the top phylum and genus; (D) bacterial genera with significant differences among groups, different letters indicate a significant difference between the groups.

Variation in Bacterial Genus in Different Treatment Groups

In terms of specific bacterial microbiota, 5,600 ASVs were identified; 47 ASVs were shared by all groups, 27 ASVs existed only in the L and ML groups, 15 ASVs were colonized in the MH and H groups, 5 ASVs occurred in just the L and MH groups, and 1 ASV existed only in the ML and H groups. Different groups had different bacterial microbiota. The specific ASVs in each group from large to small are as follows: H (51) > L (37) > MH (31) > ML (23) (Figure 3A). To explore the key microbiota between the groups, the bacterial genera at different protein levels were compared (Figures 3B–G). As protein increased, microbial fluctuations became larger, with changes mainly occurring in Bacteroidetes, Firmicutes, and Proteobacteria. The ASVs belonging to Bacteroidetes varied more widely than the other microbes. In low-protein groups, a striking enrichment occurred in most Prevotellaceae (family level), whereas there were less Bacteroidales_BS11_gut_group and Rikenellaceae. Moreover, the Bacteroidales_UCG-001 and Christensenellaceae showed a decreasing trend in the low-protein group (Supplementary Table 2).

FIGURE 3

Figure 3. Bacterial variation comparison analysis of rumen samples between yaks fed diets with different protein levels. (A) UpSet plot of the ASV count in each group. Comparative analysis between two groups: (B) L vs. ML; (C) L vs. MH; (D) L vs. H; (E) ML vs. MH; (F) ML vs. H; (G) MH vs. H.

Identification, Quantification, and Statistical Comparison of GC-TOF/MS Metabolites in the Rumen

Supplementary Figure 2 shows the Pearson correlation coefficient between the QC samples. The results confirm the reliable repeatability and precision of the data obtained in this study.

Through GC-MS detection and Chroma TOF software (used for quality control and identification), 202 valid peaks were obtained from the four groups (Supplementary Table 3). These valid peaks were mainly organic acids, lipids and lipid-like molecules, organoheterocyclic compounds, organic oxygen compounds, phenylpropanoids and polyketides, and benzenoids.

Score plots were created using the OPLS-DA model (Figures 4A–E) to verify the differentiated metabolites between the four groups, and the validation plots (Supplementary Figure 3) are also provided. The corresponding R²Y-values of the OPLS-DA models in L vs. ML, L vs. MH, L vs. H, ML vs. MH, and ML vs. H were 0.999, 0.987, 0.849, 0.999, and 0.996, respectively. The MH vs. H the model was not available (Q²Y < 0.5). The R²-values, except MH and H, were all > 0.80, indicating a satisfactory effectiveness of the model. All the samples in the score plots were within the 95% Hotelling T2 ellipse. The score plots revealed clear separation and discrimination between the different protein levels. This indicates that the OPLS-DA model can be used to identify differences between the groups (except MH and H groups).

FIGURE 4

Figure 4. OPLS-DA score plots derived from the GC-TOF/MS metabolite profiles of rumen samples between yak fed diets with different protein levels. (A) The L group vs. ML group; (B) the L group vs. MH group; (C) the L group vs. H group; (D) the ML group vs. MH group; (E) the ML group vs. H group.

Metabolomic Profiles in the Rumen

The combined statistical analysis and VIP value identified 44 differential metabolites (p < 0.01 and VIP > 1). Metabolites were classified according to the properties of the compounds, the super class, the mean value, the standard error of the mean (SEM), and the q values represented in Table 2. In general, our results show that the main metabolite differences between the four groups were alterations of organic acids and derivatives (12 metabolites), lipids and lipid-like molecules (6 metabolites), organic oxygen compounds (3 metabolites), phenylpropanoids and polyketides (7 metabolites), benzenoids (4 metabolites), organoheterocyclic compounds (2 metabolites), and others (12 metabolites). With increased protein, 3 metabolites (thymol, piceatannol, and hesperitin) cubically increased (q < 0.05), biuret significantly decreased, and 12 metabolites (hydroxyurea, 3-hydroxynorvaline, 2-ketobutyric acid, glutaric acid, n, n-dimethylarginine, 4-hydroxybutyrate, dihydrocoumarin, benzoic acid, 4-hydroxybenzoic acid, atrazine-2-hydroxy, l-threose, and n, n-dimethyl-p-phenylenediamine) linearly decreased; the ML group had the highest value and the MH group had the lowest value (q < 0.05). Six metabolites [capric acid, 3-phenyllactic acid, hydrocinnamic, 3-(4-hydroxyphenyl) propionic acid, norleucine, and allylmalonic] linearly changed (q < 0.05), three metabolites (citramalic acid, maleimide, and 9-fluorenone) quadratically changed, and 4-aminobutyric acid cubically changed (q < 0.05) (Table 2). Other responses included a mix of linear, quadratic, and cubic increases and decreases, as shown in Table 2.

TABLE 2

Table 2. Comparisons of ruminal metabolites that significantly changed with yak fed diets containing different protein level ratios.

Metabolic Pathways of Differential Metabolites

Differential metabolites in rumen samples from the four groups were analyzed using MetaboAnalyst 4.0 see text footnote 5to reveal their association with metabolic pathways. Six pathways, namely, the citrate cycle (TCA cycle), glyoxylate and dicarboxylate metabolism, butanoate metabolism, synthesis and degradation of ketone bodies, tyrosine metabolism, and alanine, aspartate, and glutamate metabolism, were identified as significantly different pathways (p < 0.05) (Figure 5). According to the KEGG pathway identification, five metabolites (3-hydroxybutyric acid, 4-aminobutyric acid, isocitric acid, L-DOPA, and noradrenaline) were mapped into these six significant pathways and were therefore identified as key metabolites among the total 44 DEMs.

FIGURE 5

Figure 5. Metabolome view map of the differentially expressed metabolites identified. The deep color shows higher pathway impact values; the larger size demonstrates higher pathway enrichment.

Correlations Between the Ruminal Metabolome and Microbiome

The correlation analysis of ruminal microbiota and metabolites is provided in Figure 6. Seventeen metabolites (from 44 DEMs) and 48 ASVs (from difference analysis) were screened. The isocitric acid at the center was associated with 25 ASVs (14 positive and 11 negative) and exhibited a positive correlation with most Prevotellaceae, but a negative correlation with Bacteroidales_BS11_gut_group, Bacteroidales_S24-7_group, Christensenellaceae, and Lachnospiraceae. 3-Phenyllactic acid and 4-Hydroxybenzoic acid showed a positive correlation with Prevotellaceae. 3-Phenyllactic acid was negatively associated with Bacteroidales_BS11_gut_group but 4-Hydroxybenzoic acid was positively associated. In addition, galactonic acid was positively associated with Bacteroidales_BS11_gut_group and Christensenellaceae, whereas it was negatively associated with most Prevotellaceae. The correlations between other metabolites are shown in Figure 6.

FIGURE 6

Figure 6. Correlation analysis between the Ruminal Metabolome and Microbiome. The correlations were based on Spearman correlation coefficients (p < 0.05). The brown lines to green lines provide the corresponding change from positive to negative correlations.

Discussion

Crude protein is a vital factor in the growth and reproduction of ruminants. Improving the dietary protein level not only promotes yak growth, and the activity of several enzymes, but improves meat composition in Tibetan sheep (Zhang et al., 2014; Wang et al., 2020). In our study, we found that the H group had a negative effect on bacterial richness and alpha diversity, suggesting that a high-quality diet may support lower bacterial diversity, but that lower diversity is not necessarily conducive to a healthy rumen environment (Latham et al., 2018), which is consistent with previous research (Petri et al., 2013). Beta diversity showed a more pronounced specificity cluster, but there was a partial overlap between the L and ML groups, indicating that the two low-protein groups had similar rumen microbiota. Conversely, the MH group and the H group formed distinct rumen flora.

Similar to previous studies (de Menezes et al., 2011; Li et al., 2018; Qiu et al., 2019), Bacteroidetes, Firmicutes, and Proteobacteria were found to be the most abundant bacteria; their higher proportion allowed their members an advantage in the use of nutrients. Among them, Bacteroidetes was significantly lower in the MH group than in the L and ML groups, which concurs with previous results where a high nutrient diet decreased the abundance of Bacteroidetes (Ley et al., 2006). However, Firmicutes and Proteobacteria did not vary significantly with protein levels in our study. The F/B ratio is considered a useful indicator of the overall gut bacteria balance and may be involved in maintaining energy homeostasis and suppressing opportunistic pathogens (Guo et al., 2019). For ruminants, the F/B ratio can reflect the host’s ability to absorb and store energy (Turnbaugh et al., 2006; Guo et al., 2008). In this study, the F/B ratio was highest in the MH group; the protein and other nutrients in the diet may have improved the energy absorption capacity of the ruminal bacteria.

At the genus level, Prevotella_1, Rikenellaceae_RC9_gut_group, and Christensenellaceae_R-7_group were the top three species. Prevotellaceae_UCG-003, Ruminococcaceae_UCG-014, and Ruminococcaceae_UCG-010 were also detected. Prevotella is known as the most abundant bacterial genus in the rumen (Stevenson and Weimer, 2007; Dan et al., 2016) and can produce large amounts of SCFAs by metabolizing dietary fiber in plant cell walls and promote host absorption of monosaccharides by interacting with complex dietary polysaccharides (Ellekilde et al., 2014; Ramayo-Caldas et al., 2016). In the present study, Prevotella_1 was significantly enriched in the low-protein groups, which could be the result of the higher fiber in the diets, promoting Prevotella to efficiently use hemicelluloses to provide more energy for the yaks (Osborne and Dehority, 1989). The Rikenellaceae_RC9_gut_group belongs to Rikenellaceae, whose members are producers of acetic and propionic acid short-chain fatty acids and ensure that the rumen contains a high percentage of propionic acid (Su et al., 2014). An increase in Rikenellaceae_RC9_gut_group in the MH and H groups confirms its important role in protein fermentation (Yang et al., 2020). The Christensenellaceae_R-7_group was highest in the MH group. Previous research showed that this genus is more abundant in an alfalfa group than a rice straw group, which indicates a broad link between Christensenellaceae and health with improved digestive function (Liu et al., 2016). Ruminococcaceae_UCG-014 and Ruminococcaceae_UCG-010 are known to degrade cellulose and protein (FitzGerald et al., 2015), which may explain their enrichment in the high-protein groups.

To further investigate the bacterial characterization in each group, we conducted a differential analysis and discovered that there were more bacterial microbiota present in the H group than in any other group. With this high level of protein, the rumen bacteria were different and less consistent than in the other groups. We also found that in the low-protein groups, most of the Prevotellaceae were enriched, while some Bacteroidales_BS11_gut_group were depleted. Prevotellaceae are responsible for metabolizing dietary fiber from plant cell walls, producing succinic acids for host absorption, and enhancing utilization of the diet (Purushe et al., 2010). Prevotellaceae are more abundant in a forage diet (Liu C. et al., 2019). Our study supports this finding. Bacteroidales_BS11_gut_group plays an important role in the fermentation of active hemicellulose monomeric sugars and the production of short-chain fatty acids (Solden et al., 2017). We speculated that high-protein diets have an effect on this process, which is consistent with previous research (Liu C. et al., 2019). Additionally, Rikenellaceae have the ability to ferment proteins, explaining their higher abundance in the high-protein groups.

The OPLS-DA analysis showed that dietary protein significantly altered the ruminal metabolite composition, suggesting that rumen metabolism is closely related to nutrition. The trends or regression lines between the different groups showed that organic acids were the most significantly affected metabolites in the present study. Lipids, phenolic acids, and aromatic compounds were also detected. In the intestine, proteins are degraded into amino acids, which are then degraded by the microbial flora into a series of downstream products that enter the liver. Amino acids that are not absorbed by the intestine can be fermented by the bacteria to produce organic acids, which are thought to be the intermediate products in the fermentation of sugars by bacteria (Zhao, 2013). The large variation in organic acids in the present study indicates that organic acid metabolism is highly sensitive to nutrient changes. Phenolic acids, such as 3-phenyllactic acid and 3-(4-hydroxyphenyl) propionic acid, are important antioxidants in the host and changes in their metabolic levels imply the changes in the bacteria and the host’s antioxidant capacity (Zhao, 2013). Saleem et al. (2012) reported a decreasing trend in hydrocinnamic acid in high-concentrate grain diets, which is similar to our results. Studies have verified that phytanic acid is derived from chlorophyll and ruminal microorganisms can form phytanic acid by hydrogenation and oxidation of the intermediate product, dihydrophytol. This allows for the slow accumulation of phytanic acid in the body of ruminants, which, as detected in our study, should be attributed to degrade cellulose bacteria (Hansen, 1966; Dawson et al., 1974).

Another metabolite, L-DOPA, is produced by plants and affects the secretion of growth hormones in steers (Kasuya et al., 2017). In the rumen, tyrosine hydroxylase may catalyze the conversion of tyrosine to L-DOPA, allowing the bacteria to use the nitrogen in L-DOPA as a nutrient source in the rumen. L-DOPA increases when the concentration of easily degradable components in the feed is reduced, which may explain the increase on the H diet (Chikagwa-Malunga et al., 2009). However, our knowledge of ruminant metabolites is far less than that of other creatures (Foroutan et al., 2020), and some DEMs detected in our study have not been reported and will require further research.

In terms of further pathway enrichment, the TCA cycle had the largest pathway impact. Isocitric acid and citramalic acid have been established as key players in this pathway, and the two metabolites were significantly decreased on the high-protein diet, appearing that they may be rapidly produced and utilized in this group, promoting the TCA cycle as well as energy conversion (Zhang et al., 2020). Furthermore, isocitric acid also participates in glyoxylate and dicarboxylate metabolism. This pathway is closely associated with the TCA cycle (Lin et al., 2018). A high-protein diet can therefore promote the conversion and absorption of energy in yaks. In addition, the synthesis and degradation of ketone bodies is one of the metabolic pathways associated with the degradation of fatty acids (Zhou et al., 2020). Other pathways related to amino acid metabolism were significantly enriched in this metabolic pathway, indicating that varying dietary protein has a considerable impact on the metabolism of amino acids.

In ruminants, the TCA cycle is central to cellular energy metabolism and chemical synthesis and synthesizes and converts glucose, amino acids, and fatty acids, inter alia (Zhu, 2015). In the present study, we observed that isocitric acid (an intermediate product of the TCA cycle) was positively correlated with Prevotellaceae. Some members of Prevotellaceae can improve feed utilization and produce VFAs to provide more energy, which has a positive effect on the TCA cycle. Therefore, it makes sense that Prevotellaceae were enriched in the L group and were correlated with isocitric acid. Bacteroidales_BS11_gut_group belongs to Bacteroidales, an important specialized anaerobe component of the rumen microbiota that converts pyruvate to acetyl-CoA (Hwang et al., 2017), a key intermediate metabolite of the TCA cycle. This family was negatively correlated with isocitric acid, probably because there is a depletion relationship between the TCA cycle and acetyl-CoA. Galactonic acid (the main component of pectin) was positively associated with Christensenellaceae in this study, suggesting that Christensenellaceae may be useful in degrading dietary polysaccharides. Moreover, 3-phenyllactic acid has antimicrobial properties and is related to phenylalanine metabolism, and phenylalanine can be catalyzed by phenylalanine hydroxylase to produce tyrosine (Yang et al., 2018; Trost et al., 2020). Our results suggest that Prevotellaceae and Bacteroidales exert positive and negative effects on the tyrosine metabolic pathway, respectively. A bioactive ingredient in some plants, 4-hydroxybenzoic acid was found to have anti-inflammatory effects in animal trials and human intestinal research (Abbas et al., 2010; Khan et al., 2016). In our study, we also detected these bioactive ingredients, which may be contained in some components of the diet and metabolized in the yak. Unfortunately, their relationships with bacteria require further investigation.

Throughout this study, there were intricate relationships between yak rumen microbes and their metabolites, and both were influenced by dietary protein levels. A previous study indicated that approximately 55–60% of rumen fluid metabolites are associated with the microbiota (Saleem et al., 2013). Other studies have also demonstrated a utilization or productive link between the composition of rumen bacteria and the rumen metabolome (Liu C. et al., 2019; Wang et al., 2019). The varying changes and relationships within the yak rumen in this study revealed important features under the impact of dietary protein level.

In summary, this study combined the microbiome and metabolomics to analyze the correlation between differential bacteria and differential metabolites in the rumen of yaks, and these associations and variations were directly related to dietary protein concentrations. This information enhances our understanding of yak ruminal bacteria and metabolites and provides more information for the development of this field and the protein level requirements within yak diets. In addition, the causes and mechanisms driving the interactions between rumen bacteria and rumen metabolism are worth further investigation.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA706154, accession: PRJNA706154.

Ethics Statement

The animal study was reviewed and approved by Institution of Animal Care and the Ethics Committee of the Northwest Institute of Plateau Biology, Chinese Academy of Sciences (NWIPB20160302).

Author Contributions

SX, TX, and XZ: Conception and experiment design. TX, XZ, XW, and YG: experiment conduction. XZ and XW: statistical analysis. SK: resources. XZ: writing – original draft preparation. SX, TX, XW, NZ, LH, and LH: writing – review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant/Award Numbers: XDA2005010406 and XDA23060603), the National Key Research and Development Program of China (Grant/Award Number: 2018YFD0502301), the Joint research project of Sanjiangyuan National Park (Grant/Award Number: YHZX-2020-7), and the Science and Technology Support Program of Qinghai Province (Grant/Award Number: 2019-SF-149).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.684340/full#supplementary-material

Footnotes

1. ^ http://www.ehbio.com/ImageGP/
2. ^ https://www.ncbi.nlm.nih.gov/
3. ^ https://bovinedb.ca/
4. ^ https://www.kegg.jp/
5. ^ https://www.metaboanalyst.ca/
6. ^ https://gephi.org/

References

Abbas, S., Greige-Gerges, H., Karam, N., Piet, M. H., Netter, P., and Magdalou, J. (2010). Metabolism of parabens (4-hydroxybenzoic acid esters) by hepatic esterases and UDP-glucuronosyltransferases in man. Drug Metab. Pharmacokinet 25, 568–577. doi: 10.2133/dmpk.dmpk-10-rg-013

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrews, S. (2013). FastQC A Quality Control tool for High Throughput Sequence Data. Babraham: Babraham Bioinformatics.

Google Scholar

AOAC (2007). Official Methods of Analysis. Gaithersburg, MD: Association of Official Analytical Chemists.

Google Scholar

Bastian, M., Heymann, S., and Jacomy, M. (2009). “Gephi: An Open Source Software for Exploring and Manipulating Networks,” in Proceedings of the Third International Conference on Weblogs and Social Media, (San Jose: ICWSM).

Google Scholar

Chikagwa-Malunga, S. K., Adesogan, A. T., Salawu, M. B., Szabo, N. J., Littell, R. C., Kim, S. C., et al. (2009). Nutritional characterization of Mucuna pruiriens. Anim. Feed Sci. Technol. 148, 51–67. doi: 10.1016/j.anifeedsci.2008.03.005

CrossRef Full Text | Google Scholar

Council, N. R. (2007). Nutrient Requirements of Small Ruminants: Sheep, Goats, Cervids, and New World Camelids. Washington, DC: The National Academies Press.

Google Scholar

Dan, X., Chen, H., Chen, F., He, Y., Zhao, C., Zhu, D., et al. (2016). Analysis of the rumen bacteria and methanogenic archaea of yak (Bos grunniens) steers grazing on the Qinghai-Tibetan Plateau. Livestock Sci. 188, 61–71. doi: 10.1016/j.livsci.2016.04.009

CrossRef Full Text | Google Scholar

Dawson, R. M. C., Hemington, N., Grime, D., Lander, D., and Kemp, P. (1974). LIPOLYSIS AND HYDROGENATION OF GALACTOLIPIDS AND ACCUMULATION OF PHYTANIC ACID IN RUMEN. Biochem. J. 144, 169–171. doi: 10.1042/bj1440169

PubMed Abstract | CrossRef Full Text | Google Scholar

de Menezes, A. B., Lewis, E., O’Donovan, M., O’Neill, B. F., Clipson, N., and Doyle, E. M. (2011). Microbiome analysis of dairy cows fed pasture or total mixed ration diets. FEMS Microbiol. Ecol. 78, 256–265. doi: 10.1111/j.1574-6941.2011.01151.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, L. M., Wang, Y. P., Brosh, A., Chen, J. Q., Gibb, M. J., Shang, Z. H., et al. (2014). Seasonal heat production and energy balance of grazing yaks on the Qinghai-Tibetan plateau. Anim. Feed Sci. Technol. 198, 83–93. doi: 10.1016/j.anifeedsci.2014.09.022

CrossRef Full Text | Google Scholar

Drackley, J. K., Donkin, S. S., and Reynolds, C. K. (2006). Major advances in fundamental dairy cattle nutrition. J. Dairy Sci. 89, 1324–1336. doi: 10.3168/jds.S0022-0302(06)72200-7

CrossRef Full Text | Google Scholar

Drexler, D. M., Reily, M. D., and Shipkova, P. A. (2011). Advances in mass spectrometry applied to pharmaceutical metabolomics. Analyt. Bioanalyt. Chem. 399, 2645–2653. doi: 10.1007/s00216-010-4370-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C., and Flyvbjerg, H. (2015). Error filtering, pair assembly and error correction for next-generation sequencing reads. Bioinformatics 31, 3476–3482. doi: 10.1093/bioinformatics/btv401

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellekilde, M., Krych, L., Hansen, C. H., Hufeldt, M. R., Dahl, K., Hansen, L. H., et al. (2014). Characterization of the gut microbiota in leptin deficient obese mice - Correlation to inflammatory and diabetic parameters. Res. Vet. Sci. 96, 241–250. doi: 10.1016/j.rvsc.2014.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

FitzGerald, J. A., Allen, E., Wall, D. M., Jackson, S. A., Murphy, J. D., and Dobson, A. D. (2015). Methanosarcina Play an Important Role in Anaerobic Co-Digestion of the Seaweed Ulva lactuca: Taxonomy and Predicted Metabolism of Functional Microbial Communities. PLoS One 10:e0142603. doi: 10.1371/journal.pone.0142603

PubMed Abstract | CrossRef Full Text | Google Scholar

Foroutan, A., Fitzsimmons, C., Mandal, R., Piri-Moghadam, H., Zheng, J., Guo, A., et al. (2020). The Bovine Metabolome. Metabolites 10, metabo10060233. doi: 10.3390/metabo10060233

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, B., Li, D., Zhou, B., Jiang, Y., Bai, H., Zhang, Y., et al. (2019). Comparative characterization of bacterial communities in geese consuming of different proportions of ryegrass. PLoS One 14:e0223445. doi: 10.1371/journal.pone.0223445

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Xia, X., Tang, R., and Wang, K. (2008). Real-time PCR quantification of the predominant bacterial divisions in the distal gut of Meishan and Landrace pigs. Anaerobe 14, 224–228. doi: 10.1016/j.anaerobe.2008.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, R. P. (1966). PHYTANIC ACID AND OTHER BRANCHED-CHAIN FATTY ACID CONSTITUENTS OF BOVINE RUMEN BACTERIA. J. Dairy Res. 33:333. doi: 10.1017/s0022029900012024

CrossRef Full Text | Google Scholar

Hwang, N., Eom, T., Gupta, S. K., Jeong, S. Y., Jeong, D. Y., Kim, Y. S., et al. (2017). Genes and Gut Bacteria Involved in Luminal Butyrate Reduction Caused by Diet and Loperamide. Genes 8:genes8120350. doi: 10.3390/genes8120350

PubMed Abstract | CrossRef Full Text | Google Scholar

Jami, E., Israel, A., Kotser, A., and Mizrahi, I. (2013). Exploring the bovine rumen bacterial community from birth to adulthood. ISME J. 7, 1069–1079. doi: 10.1038/ismej.2013.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasuya, E., Sutoh, M., and Yayou, K. I. (2017). The effects of l-DOPA and sulpiride on growth hormone secretion at different injection times in Holstein steers. Anim. Sci. J. 88, 1842–1848. doi: 10.1111/asj.12850

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, S. A., Chatterjee, S. S., and Kumar, V. (2016). Low dose aspirin like analgesic and anti-inflammatory activities of mono-hydroxybenzoic acids in stressed rodents. Life Sci. 148, 53–62. doi: 10.1016/j.lfs.2016.02.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, M., Morrison, M., and Yu, Z. (2011). Status of the phylogenetic diversity census of ruminal microbiomes. FEMS Microbiol. Ecol. 76, 49–63. doi: 10.1111/j.1574-6941.2010.01029.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kind, T., Wohlgemuth, G., Lee, D. Y., Lu, Y., Palazoglu, M., Shahbaz, S., et al. (2009). FiehnLib: Mass Spectral and Retention Index Libraries for Metabolomics Based on Quadrupole and Time-of-Flight Gas Chromatography/Mass Spectrometry. Analyt. Chem. 81, 10038–10048. doi: 10.1021/ac9019522

PubMed Abstract | CrossRef Full Text | Google Scholar

Krone, N., Hughes, B. A., Lavery, G. G., Stewart, P. M., Arlt, W., and Shackleton, C. H. L. (2010). Gas chromatography/mass spectrometry (GC/MS) remains a pre-eminent discovery tool in clinical steroid investigations even in the era of fast liquid chromatography tandem mass spectrometry (LC/MS/MS). J. Steroid Biochem. Mol. Biol. 121, 496–504. doi: 10.1016/j.jsbmb.2010.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Latham, E. A., Weldon, K. K., Wickersham, T. A., Coverdale, J. A., and Pinchak, W. E. (2018). Responses in the rumen microbiome of Bos taurus and indicus steers fed a low-quality rice straw diet and supplemented protein. J. Anim. Sci. 96, 1032–1044. doi: 10.1093/jas/sky023

PubMed Abstract | CrossRef Full Text | Google Scholar

Ley, R. E., Turnbaugh, P. J., Klein, S., and Gordon, J. I. (2006). Microbial ecology - Human gut microbes associated with obesity. Nature 444, 1022–1023. doi: 10.1038/4441022a

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Mehmood, K., Zhang, H., Jiang, X., Shahzad, M., Dong, X., et al. (2018). Characterization of fungus microbial diversity in healthy and diarrheal yaks in Gannan region of Tibet Autonomous Prefecture. Acta Trop. 182, 14–26. doi: 10.1016/j.actatropica.2018.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Z., Ye, W., Zu, X., Xie, H., Li, H., Li, Y., et al. (2018). Integrative metabolic and microbial profiling on patients with Spleen-yang-deficiency syndrome. Sci. Rep. 8:6619. doi: 10.1038/s41598-018-24130-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Wu, H., Liu, S., Chai, S., Meng, Q., and Zhou, Z. (2019). Dynamic Alterations in Yak Rumen Bacteria Community and Metabolome Characteristics in Response to Feed Type. Front. Microbiol. 10:1116. doi: 10.3389/fmicb.2019.01116

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Zhang, M., Xue, C., Zhu, W., and Mao, S. (2016). Characterization and comparison of the temporal dynamics of ruminal bacterial microbiota colonizing rice straw and alfalfa hay within ruminants. J. Dairy Sci. 99, 9668–9681. doi: 10.3168/jds.2016-11398

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L. S., Aa, J. Y., Wang, G. J., Yan, B., Zhang, Y., Wang, X. W., et al. (2010). Differences in metabolite profile between blood plasma and serum. Analyt. Biochem. 406, 105–112. doi: 10.1016/j.ab.2010.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Ding, L., Zhou, Y., Jing, X., and Degen, A. A. (2019). Behavioural characteristics of yaks grazing summer and winter pastures on the Qinghai-Tibetan Plateau. Appl. Anim. Behav. Sci. 218:007. doi: 10.1016/j.applanim.2019.06.007

CrossRef Full Text | Google Scholar

Liu, Y. X., Qin, Y., Chen, T., Lu, M., Qian, X., Guo, X., et al. (2020). A practical guide to amplicon and metagenomic analysis of microbiome data. Protein Cell 12, 315–330. doi: 10.1007/s13238-020-00724-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, R. J., Ding, L. M., Shang, Z. H., and Guo, X. H. (2008). The yak grazing system on the Qinghai-Tibetan plateau and its status. Rangeland J. 30, 241–246. doi: 10.1071/rj08012

CrossRef Full Text | Google Scholar

Ministry of Agriculture of the PRC (2004). Feeding Standard for Beef Cattle. Beijing: Ministry of Agriculture of the PRC.

Google Scholar

Misra, B. B. (2020). The Connection and Disconnection Between Microbiome and Metabolome: A Critical Appraisal in Clinical Research. Biol. Res. Nurs. 22, 561–576. doi: 10.1177/1099800420903083

PubMed Abstract | CrossRef Full Text | Google Scholar

Oakman, C., Tenori, L., Biganzoli, L., Santarpia, L., Cappadona, S., Luchinat, C., et al. (2011). Uncovering the metabolomic fingerprint of breast cancer. Int. J. Biochem. Cell Biol. 43, 1010–1020. doi: 10.1016/j.biocel.2010.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Osborne, J. M., and Dehority, B. A. (1989). Synergism in degradation and utilization of intact forage cellulose, hemicellulose, and pectin by 3 pure cultures of ruminal bacteria. Appl. Environ. Microbiol. 55, 2247–2250. doi: 10.1128/aem.55.9.2247-2250.1989

PubMed Abstract | CrossRef Full Text | Google Scholar

Petri, R. M., Schwaiger, T., Penner, G. B., Beauchemin, K. A., Forster, R. J., McKinnon, J. J., et al. (2013). Characterization of the Core Rumen Microbiome in Cattle during Transition from Forage to Concentrate as Well as during and after an Acidotic Challenge. PLoS One 8:0083424. doi: 10.1371/journal.pone.0083424

PubMed Abstract | CrossRef Full Text | Google Scholar

Purushe, J., Fouts, D. E., Morrison, M., White, B. A., Mackie, R. I., et al. (2010). Comparative genome analysis of Prevotella ruminicola and Prevotella bryantii: insights into their environmental niche. Microb. Ecol. 60, 721–729. doi: 10.1007/s00248-010-9692-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, Q., Gao, C., Gao, Z., Rahman, M. A. U., He, Y., Cao, B., et al. (2019). Temporal Dynamics in Rumen Bacterial Community Composition of Finishing Steers during an Adaptation Period of Three Months. Microorganisms 7:microorganisms7100410. doi: 10.3390/microorganisms7100410

PubMed Abstract | CrossRef Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramayo-Caldas, Y., Mach, N., Lepage, P., Levenez, F., Denis, C., Lemonnier, G., et al. (2016). Phylogenetic network analysis applied to pig gut microbiota identifies an ecosystem structure linked with growth traits. ISME J. 10, 2973–2977. doi: 10.1038/ismej.2016.77

PubMed Abstract | CrossRef Full Text | Google Scholar

Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahe, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:22. doi: 10.7717/peerj.2584

PubMed Abstract | CrossRef Full Text | Google Scholar

Saleem, F., Ametaj, B. N., Bouatra, S., Mandal, R., Zebeli, Q., Dunn, S. M., et al. (2012). A metabolomics approach to uncover the effects of grain diets on rumen health in dairy cows. J. Dairy Sci. 95, 6606–6623. doi: 10.3168/jds.2012-5403

PubMed Abstract | CrossRef Full Text | Google Scholar

Saleem, F., Bouatra, S., Guo, A. C., Psychogios, N., Mandal, R., Dunn, S. M., et al. (2013). The Bovine Ruminal Fluid Metabolome. Metabolomics 9, 360–378. doi: 10.1007/s11306-012-0458-9

CrossRef Full Text | Google Scholar

Shabat, S. K., Sasson, G., Doron-Faigenboim, A., Durman, T., Yaacoby, S., Berg Miller, M. E., et al. (2016). Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 10, 2958–2972. doi: 10.1038/ismej.2016.62

PubMed Abstract | CrossRef Full Text | Google Scholar

Solden, L. M., Hoyt, D. W., Collins, W. B., Plank, J. E., Daly, R. A., Hildebrand, E., et al. (2017). New roles in hemicellulosic sugar fermentation for the uncultivated Bacteroidetes family BS11. ISME J. 11, 691–703. doi: 10.1038/ismej.2016.150

PubMed Abstract | CrossRef Full Text | Google Scholar

Stevenson, D. M., and Weimer, P. J. (2007). Dominance of Prevotella and low abundance of classical ruminal bacterial species in the bovine rumen revealed by relative quantification real-time PCR. Appl. Microbiol. Biotechnol. 75, 165–174. doi: 10.1007/s00253-006-0802-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Styczynski, M. P., Moxley, J. F., Tong, L. V., Walther, J. L., Jensen, K. L., and Stephanopoulos, G. N. (2007). Systematic identification of conserved metabolites in GC/MS data for metabolomics and biomarker discovery. Analyt. Chem. 79, 966–973. doi: 10.1021/ac0614846

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, X. L., Tian, Q., Zhang, J., Yuan, X. Z., Shi, X. S., Guo, R. B., et al. (2014). Acetobacteroides hydrogenigenes gen. nov., sp. nov., an anaerobic hydrogen-producing bacterium in the family Rikenellaceae isolated from a reed swamp. Int. J. Syst. Evol. Microbiol. 64(Pt 9), 2986–2991. doi: 10.1099/ijs.0.063917-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Trost, K., Ahonen, L., Suvitaival, T., Christiansen, N., Nielsen, T., Thiele, M., et al. (2020). Describing the fecal metabolome in cryogenically collected samples from healthy participants. Sci. Rep. 10:885. doi: 10.1038/s41598-020-57888-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Turnbaugh, P. J., Ley, R. E., Mahowald, M. A., Magrini, V., Mardis, E. R., and Gordon, J. I. (2006). An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444, 1027–1031. doi: 10.1038/nature05414

PubMed Abstract | CrossRef Full Text | Google Scholar

Vansoest, P. J., Robertson, J. B., and Lewis, B. A. (1991). METHODS FOR DIETARY FIBER, NEUTRAL DETERGENT FIBER, AND NONSTARCH POLYSACCHARIDES IN RELATION TO ANIMAL NUTRITION. J. Dairy Sci. 74, 3583–3597. doi: 10.3168/jds.S0022-0302(91)78551-2

CrossRef Full Text | Google Scholar

Wang, B., Ma, M. P., Diao, Q. Y., and Tu, Y. (2019). Saponin-Induced Shifts in the Rumen Microbiome and Metabolome of Young Cattle. Front. Microbiol. 10:356. doi: 10.3389/fmicb.2019.00356

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Xu, T., Zhang, X., Geng, Y., Kang, S., and Xu, S. (2020). Effects of Dietary Protein Levels on Growth Performance, Carcass Traits, Serum Metabolites, and Meat Composition of Tibetan Sheep during the Cold Season on the Qinghai-Tibetan Plateau. Animals 10:ani10050801. doi: 10.3390/ani10050801

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, J., Psychogios, N., Young, N., and Wishart, D. S. (2009). MetaboAnalyst: a web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 37, W652–W660. doi: 10.1093/nar/gkp356

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C., Tsedan, G., Liu, Y., and Hou, F. (2020). Shrub coverage alters the rumen bacterial community of yaks (Bos grunniens) grazing in alpine meadows. J. Anim. Sci. Technol. 62, 504–520. doi: 10.5187/jast.2020.62.4.504

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Dong, G., Wang, Z., Wang, J., Zhang, Z., and Liu, J. (2018). Rumen and plasma metabolomics profiling by UHPLC-QTOF/MS revealed metabolic alterations associated with a high-corn diet in beef steers. PLoS One 13:e0208031. doi: 10.1371/journal.pone.0208031

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeoman, C. J., Ishaq, S. L., Bichi, E., Olivo, S. K., Lowe, J., and Aldridge, B. M. (2018). Biogeographical Differences in the Influence of Maternal Microbial Sources on the Early Successional Development of the Bovine Neonatal Gastrointestinal tract. Sci. Rep. 8:3197. doi: 10.1038/s41598-018-21440-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeisel, S. H., Freake, H. C., Bauman, D. E., Bier, D. M., Burrin, D. G., German, J. B., et al. (2005). The nutritional phenotype in the age of metabolomics. J. Nutrit. 135, 1613–1616.

Google Scholar

Zhang, H. B., Wang, Z. S., Peng, Q. H., Tan, C., and Zou, H. W. (2014). Effects of different levels of protein supplementary diet on gene expressions related to intramuscular deposition in early-weaned yaks. Anim. Sci. J. 85, 411–419. doi: 10.1111/asj.12161

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Shi, H. T., Wang, Y. C., Li, S. L., Cao, Z. J., Yang, H. J., et al. (2020). Carbohydrate and amino acid metabolism and oxidative status in Holstein heifers precision-fed diets with different forage to concentrate ratios. Animal 14, 2315–2325. doi: 10.1017/S1751731120001287

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Zhu, W., Jiang, L., and Mao, S. (2017). Comparative metabolome analysis of ruminal changes in Holstein dairy cows fed low- or high-concentrate diets. Metabolomics 13, 1204–1200. doi: 10.1007/s11306-017-1204-0

CrossRef Full Text | Google Scholar

Zhao, Y. (2013). Antibiotics-induced mice microbial alteration and their relationship with host metabolic profiles. Huazhong: Huazhong University of Science and Technology.

Google Scholar

Zhou, Y., Li, M., Song, J., Shi, Y., Qin, X., Gao, Z., et al. (2020). The cardioprotective effects of the new crystal form of puerarin in isoproterenol-induced myocardial ischemia rats based on metabolomics. Sci. Rep. 10:17787. doi: 10.1038/s41598-020-74246-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, S. Q. (2015). TCA cycle for forward and reverse rotation. Bull. Biol. 50, 16–19.

Google Scholar