- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
The rumen microbiome plays a central role in ruminant nutrition and health. To investigate the effects of different feeding regimens on it, this study employed multi-omics analysis to reveal how natural grazing versus intensive feeding alters the rumen microbiota and metabolites in Shorthorn cattle. A total of 18 male shorthorn cattle of about 17 months of age and similar body weight were selected and randomly divided into 3 groups: natural grazing bull group (DJCF), intensive feeding bull group (DJCY) and intensive feeding steer group (DJC). The experiment period was 361 days. After the fattening trial, rumen fluid was collected at slaughter. Microbiota and metabolites were analyzed by 16S rRNA sequencing and LC–MS, and correlations were assessed. The results indicate that different feeding regimens were strongly associated with shifts in rumen microbial diversity and community composition. The ACE and Shannon indices of DJCF group were significantly higher than those of DJCY and DJC group (p < 0.05). Bacteroidetes and Firmicutes were the dominant phyla, with relative abundances of 57.62% (DJCF), 54.11% (DJCY), 48.84% (DJC) and 34.07, 38.31, 43.08%, respectively, showing no significant differences (p > 0.05). At the genus level, Prevotella and Rikenellaceae_RC9_gut_group were dominant. The abundance of Prevotella was highest in DJCY (22.52%), significantly differing from DJC (12.43%; p < 0.05), while Rikenellaceae_RC9_gut_group abundances were 12.56% (DJCF), 9.92% (DJCY), and 11.89% (DJC). In the fungal community, Neocallimastigomycota and Ascomycota were the dominant phyla, and there were no significant differences among the three groups. At the genus level, Caecomyces, the highest in the DJC group, with a significant difference from the DJCF group (p < 0.05). Orpinomyces, the highest in the DJCF group, with significant differences from the DJCY and DJC groups (p < 0.05). There were significant differences in rumen metabolites between different groups, and a variety of different metabolites were identified, involving sucrose and starch metabolism, purine metabolism and other pathways (p < 0.05). In addition, there was a significant correlation between rumen microbes and metabolites (p < 0.05). Thus, an intensive feeding system altered the rumen microbiome, resulting in improvements of Shorthorn cattle growth. Nevertheless, the specific causal relationships and underlying regulatory mechanisms governing the interplay between rumen microbiota and metabolic processes remain to be further elucidated through in-depth investigations.
1 Introduction
The Shorthorn cattle, a prominent breed of beef cattle, originated in the historic counties of Northumberland, Durham, Yorkshire, and Lincolnshire in northeastern England. Empirical evidence suggests that through systematic selective breeding during the late Victorian agricultural revolution, Shorthorn cattle evolved into a globally recognized beef genotype by the early 20th century—a period characterized by proto-industrialization (Deerry, 2018). In comparison to indigenous Chinese cattle breeds, Shorthorn cattle demonstrate superior growth rates and meat quality. This makes them an essential genetic resource for enhancing beef production in China, particularly through hybrid breeding programs.
Traditionally, these cattle are raised on artificial grasslands. However, during the six-month forage withering period, their weight experiences a significant decline, which adversely affects subsequent fattening efficiency. Therefore, it is essential to optimize feeding strategies to mitigate these seasonal nutritional deficits in order to sustain Shorthorn cattle production. By aligning feed composition with the physiological needs of the cattle, intensive feeding regimes can reduce nutrient wastage, enhance metabolic efficiency, and minimize risks of environmental contamination (Moorby and Fraser, 2021).
The microbial ecosystem in the rumen plays a vital role in nutrient digestion and metabolism, with its composition influenced by factors such as host genetic background, developmental stage, dietary composition, and environmental conditions (Lin et al., 2019; Zhang et al., 2018). Changes in dietary nutrient composition can alter the abundance of specific microbial populations. For instance, studies have shown that the energy level of the ration significantly affects the rumen microbial composition in yaks (Ahmad et al., 2020). Additionally, feedlot fattening has been found to impact rumen fermentation, microbial diversity, and meat quality (Hu et al., 2021). Furthermore, Yang et al. (2024) observed that intensive feeding regimes led to an increase in metabolites associated with amino acid metabolism, such as L-glutamate (Yang et al., 2024). The proportion of dietary concentrate has also been reported to influence rumen microbial abundance, diversity, metabolite concentrations, and metabolic pathways (Pang et al., 2022). To gain deeper insights into these diet-driven microbial and metabolic dynamics, multi-omics approaches, integrating genomic, transcriptomic, and metabolomic data, have proven essential. These methods help link functional gene activity with microbial communities and metabolic processes within the rumen (Denman et al., 2018). By applying these approaches, researchers can better understand how dietary alterations affect microbial populations and metabolic functions, leading to a more comprehensive understanding of the rumen ecosystem in ruminants.
Despite significant advancements in ruminant nutrition, Shorthorn cattle, with their distinct genetic and metabolic characteristics, remain underexplored. Most existing studies on Shorthorn cattle have focused on genetic traits and growth performance, often overlooking the functional role of rumen microbiota and metabolites in nutrient utilization. To address this gap, our study integrates 16S/ITS sequencing and LC–MS metabolomics to investigate: (1) the impact of natural grazing versus intensive feeding on rumen microbial diversity and metabolite profiles in Shorthorn cattle, and (2) the relationships between microbial communities and key metabolites and metabolic pathways. This research represents the first multi-omics characterization of the Shorthorn cattle rumen ecosystem, providing valuable insights that can help refine feeding strategies and improve production efficiency.
2 Materials and methods
2.1 Animals, diets, and experimental design
The study was conducted at the Yunnan Provincial Breeding and Promotion Centre for Breeding Sheep in Xundian County, located in the eastern suburbs of Kunming, approximately 90 km from the city (longitude 103°11′ east, latitude 25°40′ north). The site is situated in a low hilly area at an altitude of 2,040 m, with an average annual temperature of 13.4 °C and average annual precipitation of 1,025 mm. A total of 18 male Shorthorn cattle, approximately 17 months old and with similar body weights, were selected and randomly assigned to three groups: the natural grazing bull group (DJCF, n = 6), the intensive feeding bull group (DJCY, n = 6), the intensive feeding castrated bull group (DJC, n = 6). Castration was performed in early March to allow full recovery before the commencement of the fattening trial. The DJCF group grazed on artificially improved grassland, while the DJCY and DJC groups were raised in enclosures and freely fed a total mixed ration (TMR). All three groups were provided with ad libitum access to water. The nutritional composition of the DJCY group and DJC group are shown in Table 1. The experimental period lasted 361 days, including a 10-day pretest phase and a 351-day main testing phase. During the pre-feeding period, bulls and castrated bulls in the intensive feeding groups were gradually acclimated to the TMR diet. Rumen fluid was collected at slaughter following the completion of the fattening trial.
2.2 Sample collection and measurements
On the day preceding the conclusion of the experiment, all experimental animals were subjected to a 24-h feed withdrawal, followed by an 8-h water restriction, in order to standardize the conditions prior to the collection of rumen fluid. Following the guidelines outlined in the Chinese National Standard GB/T 19477–2004 for bovine slaughtering protocols, qualified personnel performed dissection of the gastrointestinal tract. The rumen was quickly removed to collect the rumen fluid. After filtration through sterile gauze, rumen fluid was aliquoted into three 10 mL cryovials, snap-frozen in liquid nitrogen, and stored at −80 °C for microbiota and metabolite analysis.
2.3 DNA extraction and sequencing
DNA was extracted from the microbial community in rumen fluid, stored at −80 °C, using the TIANamp Stool DNA Kit (TianGen, Beijing, China, catalog: DP712), following the manufacturer’s protocol. The rumen fluid samples were subsequently transported under cryogenic conditions to Yunnan Pulis Biotechnology Co., Ltd. (Kunming, China) for high-throughput sequencing. The 16S rRNA and ITS genes from target regions were PCR-amplified with region-specific primers containing unique barcodes, using barcoded specific primers 515F (515F:5’-GTGYCAGCMGCCGCGGTAA-3′) and 806R (5′- GTGCCAGCMGCCGCGG-3′) to amplify the V3 – V4 hypervariable region of 16S rRNA gene. ITS1 region primers (ITS1-5F: GGAAGTAAAAGTCGTAACAAGG, GCTGCGTTCTTCATCGATGC): Identification of fungal diversity. PCR reaction system and conditions: Each PCR reaction system contains 15 μL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μM primers, and 10 ng of genomic DNA template. The first denaturation step was performed at 98 °C for 1 min, followed by 30 cycles at 98 °C (10 s), 50 °C (30 s), and 72 °C (30 s), with a final extension at 72 °C for 5 min. PCR products were detected by 2% agarose gel electrophoresis. Qualified PCR products were then purified using magnetic beads and quantified via enzyme-linked immunosorbent assay (ELISA). PCR products were mixed in equal volumes according to their concentrations, thoroughly mixed, and subsequently detected by 2% agarose gel electrophoresis. The target band was then recovered using a universal DNA purification and recovery kit (Tian Gen). Library construction was performed using the NEB Next® Ultra™ II FS DNA PCR-free Library Prep Kit (New England Biolabs). The constructed libraries were quantified using Qubit and qPCR. Qualified libraries were sequenced using the NovaSeq 6,000 sequencing platform with PE250 sequencing.
Individual sample data were extracted from the raw sequencing data based on barcode sequences and PCR amplification primers. After demultiplexing, barcode and primer sequences were trimmed to obtain preliminary sequence data. Using FLASH (Version 1.2.111), reads from each sample were merged to generate raw tag data (Raw Tags). Paired-end reads were merged and subjected to quality control with fastp (v0.23.1) to yield high-quality clean tags. Chimeric sequences were identified and removed using UCHIME with reference to the SILVA database for 16S rRNA and the UNITE database for ITS (Knight, 2011). Quality-filtered, non-chimeric sequences were clustered into operational taxonomic units (OTUs) at 97% sequence identity in QIIME 2 (via VSEARCH) (Callahan et al., 2019), and an OTU table and representative sequences were generated (Wang et al., 2021). Taxonomic annotation was conducted in QIIME 2 against the SILVA (v138.1) reference database for 16S rRNA and UNITE (v8.2) for ITS. A phylogenetic tree was inferred from multiple-sequence-aligned OTU representative sequences. Prior to downstream analyses, samples were rarefied to the minimum sequencing depth to standardize library sizes.
2.4 Determination and analysis of rumen metabolome
Rumen metabolomics was conducted using LC–MS/MS with a quadrupole-Orbitrap mass spectrometer (Q Exactive Orbitrap) coupled to UHPLC (Thermo Fisher Scientific, USA). Transfer 100 μL of rumen fluid to a 1.5 mL microcentrifuge tube. Add 400 μL of extraction solvent (acetonitrile:methanol, 1:1, v/v) containing an isotopically labeled internal standard mixture. Vortex for 30 s, sonicate in an ice–water bath for 5 min, and incubate at −40 °C for 1 h to precipitate proteins. Centrifuge at 12,000 rpm for 15 min at 4 °C. The resulting supernatant was transferred to a clean glass vial for downstream analyses. A pooled QC sample was prepared by mixing equal aliquots of supernatant from each sample. Polar metabolite profiling was conducted using a Vanquish UHPLC system coupled to an Orbitrap Exploris 120 mass spectrometer (Thermo Fisher Scientific). Chromatographic separation was performed on a Waters ACQUITY UPLC BEH Amide column (2.1 × 50 mm, 1.7 μm) with solvent A (25 mM ammonium acetate + 25 mM ammonium hydroxide in water) and solvent B (acetonitrile). The autosampler was held at 4 °C, and 2 μL were injected. Data acquisition was carried out on an Orbitrap Exploris 120 under Xcalibur v4.4 control (Thermo Fisher Scientific).
Raw data were converted to mzXML with ProteoWizard (msConvert). Metabolite identification used BiotreeDB (v3.0) in R, and results were visualized with R package (Smoot et al., 2011). Differential metabolites were defined by OPLS-DA VIP scores together with Student’s t-test (FC > 1, p < 0.05, VIP > 1), and KEGG enrichment was tested by a hypergeometric test.
2.5 Statistical analysis
Alpha diversity was estimated in QIIME 2 using observed OTUs, Chao1 richness, and Shannon diversity. Group differences were assessed with the non-parametric Kruskal–Wallis test. Boxplots of α-diversity were generated in R (v4.4.1) with ggplot2. Principal coordinate analysis (PCoA) was performed in R using the vegan package. Ordination scatterplots were created with ggplot2, with 95% confidence ellipses to delineate groups.
The phyloseq package was used to calculate the relative abundance of microorganisms at different taxonomic levels (e.g., phylum, genus), and stacked bar charts were plotted to demonstrate the compositional characteristics of the microbial communities by ggplot2. Processing of data and prediction of function by R v4.4.1 tidyverse, microeco, aplot, ggsci, reshape2 packages.
The raw 16S rRNA/ITS sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) under accession number PRJNA1277806. To examine associations between characteristic microbial taxa and key metabolites, Spearman’s rank correlations were computed in R using the psych package. A co-occurrence network was then constructed by retaining only significant associations (|r| > 0.7, p < 0.05) and visualized in Cytoscape (v3.10.1) (Smoot et al., 2011).
3 Results
3.1 Analysis of differences in the composition of bacterial microbial communities
3.1.1 Analysis of microbial diversity differences
Dietary composition exerted a pronounced effect on microbial diversity. Alpha-diversity analyses revealed differences among the three groups: the DJCF group showed significantly higher ACE, Pielou’s evenness, and Shannon indices than DJCY and DJC (p < 0.05), whereas DJCY and DJC did not differ (p > 0.05) (Figures 1A–D). Figure 1E shows that a total of 21,232 OTUs were identified across all dietary treatments: 6,041 in DJCF, 4,182 in DJCY, and 3,727 in DJC. Pairwise overlaps included 967 OTUs shared by DJCF-DJCY, 1,330 by DJCY-DJC, and 964 by DJCF-DJC, with 4,021 OTUs common to all three groups. Principal coordinate analysis (PCoA) further separated DJCF from the other two groups based on bacterial community composition, with PC1 and PC2 explaining 24.76 and 19.76% of the variance, respectively (Figure 1F).
 
  Figure 1. Differential analysis of bacterial microbial communities and diversity. (A) Ace index; (B) good_coverage index; (C) pielou index; (D) Shannon index, significant is marked by *p < 0.05, **p < 0.01; (E) Venn diagram analysis among groups; (F) PCA analysis; The relative abundance of fungal at the phylum level (G) and genus level (H).
3.1.2 Analysis of bacterial communities in shorthorn cattle
Across dietary treatments, Bacteroidetes and Firmicutes were the dominant phyla, as shown in Figure 1G and Table 2. The relative abundances of Bacteroidetes in the DJCF, DJCY, and DJC groups were 57.62, 54.11, and 48.84%, respectively. Firmicutes was highest in the DJC group (43.08%), followed by DJCY (38.31%) and DJCF (34.07%). However, the differences in the relative abundances of Bacteroidetes and Firmicutes among the three groups were not statistically significant (p > 0.05). Among the other phyla, the relative abundance of Patescibacteria was significantly higher in DJCY than in DJCF (p < 0.05), whereas no significant difference was detected between DJCY and DJC (p > 0.05). By contrast, Fibrobacterota was significantly more abundant in DJCF than in both DJCY and DJC (p < 0.05). By contrast, Desulfobacterota was significantly less abundant in DJCF than in DJCY and DJC (p < 0.05). Dietary variation influenced the genus-level composition of the microbiota in the DJCF, DJCY, and DJC groups. Across the groups, Prevotella, Rikenellaceae_RC9_gut_group, and F082 were the most abundant taxa. Prevotella was most abundant in the DJCY group (22.52%) and was significantly more abundant than in the DJC group (p < 0.05). Rikenellaceae_RC9_gut_group was most abundant in the DJCF group (12.56%), the highest among the three groups (Figure 1H and Table 3).
 
  Table 2. Analysis of community composition at the rumen bacterial phylum level of different feeding practices in Shorthorn cattle.
 
  Table 3. Analysis of community composition at the rumen bacterial genus level of different feeding practices in Shorthorn cattle.
3.1.3 Microbial interactions between rumen bacterial communities in shorthorn cattle
Correlation analysis (Figures 2A–C) depicts co-occurrence networks of the rumen microbiota in Shorthorn cattle under different diets (DJCF, DJCY, DJC). Network analysis revealed significantly more positive than negative correlations among congeneric taxa (p < 0.05). Notably, taxa within the same genus showed strong positive associations (p < 0.05).
 
  Figure 2. Network diagram. Different nodes represent different genera, node size represents the average relative abundance of the genus, nodes of the same phylum have the same color (as indicated in the figure legend), the absolute value of the correlation coefficient between line thickness and species interactions between nodes is positive, and the correlation coefficient between line color and correlation is positive (red is positively, blue is negatively). (A) DJCF, (B) DJCY, (C) DJC.
3.2 Analysis of differences in the composition of fungal microbial communities
3.2.1 Analysis of fungal diversity differences
α-Diversity analysis (Figures 3A–D) revealed no significant differences in α-diversity indices among the three groups (p > 0.05). Venn-diagram analysis of shared and unique features (Figure 3E) identified 23,478 features across all groups, of which 939 were shared. The DJCF, DJCY, and DJC groups contained 7,483, 6,939, and 5,208 unique OTUs, respectively.
 
  Figure 3. Differential analysis of fungal microbial communities and diversity. (A) Ace index; (B) good_coverage index; (C) pielou index; (D) Shannon index, significant is marked by*p < 0.05, **p < 0.01; (E) Venn diagram analysis among groups; (F) PCA analysis; The relative abundance of fungal at the phylum level (G) and genus level (H).
β-Diversity analysis using PCoA (Figure 3F) showed that the fungal community composition of the DJCF group was distinct from the other two groups, while DJCY and DJC did not differ. The first two principal coordinates (PC1 and PC2) accounted for 46.74 and 11.02% of the total variation, respectively. Samples from the DJCF and DJCY groups formed tight clusters, indicating lower within-group variability. Conversely, samples from the DJC group were more dispersed, suggesting higher within-group variability.
3.2.2 Analysis of fungal microbial communities
A comparison of the dominant phyla among the three groups is shown in Figure 3G and Table 4. We analyzed the 10 identified fungal phyla. In the DJCF, DJCY, and DJC groups, the dominant phyla were Neocallimastigomycota and Ascomycota. Neocallimastigomycota accounted for 86.74, 75.68, and 84.88% in DJCF, DJCY, and DJC, respectively; Ascomycota accounted for 2.54, 14.25, and 2.09%, respectively. No significant differences in the relative abundances of these dominant phyla were detected among the three groups (p > 0.05). As shown in Figure 3H and Table 5, Caecomyces accounted for 14.22, 24.73, and 27.85% in DJCF, DJCY, and DJC, respectively, and was highest in DJC; its abundance in DJC was significantly higher than in DJCF (p < 0.05). Orpinomyces was most abundant in DJCF (29.83%) and was significantly higher than in DJCY (7.63%) and DJC (15.98%) (p < 0.05). Cyllamyces accounted for 10.13, 10.24, and 7.24% in DJCF, DJCY, and DJC, respectively, and Piromyces accounted for 14.26, 6.62, and 6.30%, respectively, there were no significant differences among the three groups (p > 0.05).
 
  Table 4. Rumen fungi community composition at phylum in the different feeding systems in Shorthorn cattle.
 
  Table 5. Rumen fungi community composition at genus in the different feeding systems in Shorthorn cattle.
3.2.3 Microbial interactions between rumen fungal communities in Shorthorn cattle
To elucidate the relationships among the DJCF, DJCY, and DJC microbial communities in Shorthorn cattle under different diets, we constructed co-occurrence networks (Figures 4A–C). Network analysis showed that positive correlations between species were more frequent than negative ones, particularly among congeneric taxa.
 
  Figure 4. Network diagram. Different nodes represent different genera, node size represents the average relative abundance of the genus, nodes of the same phylum have the same color (as indicated in the figure legend), the absolute value of the correlation coefficient between line thickness and species interactions between nodes is positive, and the correlation coefficient between line color and correlation is positive (red is positively, blue is negatively). (A) DJCF, (B) DJCY, (C) DJC.
3.3 Rumen metabolome associated with different feeding patterns
3.3.1 Metabolomic profiles in rumen
We profiled rumen metabolites using untargeted metabolomics. Orthogonal partial least squares–discriminant analysis (OPLS-DA) showed marked differences among the DJCF, DJCY, and DJC groups (Figures 5A–C). Pairwise comparisons were well separated in the score plots, and all samples fell within the 95% confidence ellipse. As shown in Figures 5D–F, the R2Y/Q2 values for the comparisons between the DJCF and DJCY group, the DJCF and DJC group, and the DJCY and DJC group were 0.995/0.723, 0.992/0.673, and 0.973/0.620, respectively. These results indicate good model fit (R2Y) and predictive ability (Q2), supporting the reliability of subsequent statistical analyses.
![Six-part image with scatter plots A to C and validation plots D to F. Plots A, B, and C display different group separations using t[1]P and t[1]O components. Groups DJCF, DJCY, and DJC are shown with various color codes. Plots D, E, and F display R2Y and Q2 validation with correlation coefficients. Each includes intercept values and dotted lines highlighting trends.](https://www.frontiersin.org/files/Articles/1657402/fmicb-16-1657402-HTML/image_m/fmicb-16-1657402-g005.jpg) 
  Figure 5. Corresponding validation plots and orthogonal projections to latent structures-discriminant analysis (OPLS-DA) score plots derived from gas chromatography time-of-flight mass spectrometer (GC-TOF/MS) metabolite profiles of Shorthorn cattle rumen samples with Different feeding methods. Corresponding validation plots and OPLS-DA score plots for (A,D) group DJCF vs. group DJCY, (B,E) group DJCF vs. group DJC, and (C,F) group DJCY vs. DJC.
3.3.2 Analysis of differential metabolites under different feeding practices
As shown in Figures 6A–C, 116 metabolites were downregulated and 1741 metabolites were upregulated in DJCF vs. DJCY group. Seven hundred and twenty-eight metabolites were downregulated and 2040 metabolites were upregulated in DJCF vs. DJC group. Seven hundred and seventy-three metabolites were downregulated and 1,468 metabolites were upregulated in DJCY vs. DJC group. DJCY group contained 45 metabolites with significantly elevated concentrations compared to DJCF group (p < 0.05, VIP > 1), such as Daltogen, Cytidine 5′-monophosphate (CMP), 1-Methyladenosine, Turanose, Isomaltose, Lactose, Cellobiose, Zidovudine. There were 28 metabolites with significantly different properties in the DJCF group that were significantly higher than in the DJCY group (p < 0.05), such as 2-methylbenzoic acid, 4-methylbenzoic acid, phenylacetic acid, trimethylamine N-oxide, and caprolactam. The comparison analysis in the DJCF vs. DJC group revealed that the relative concentrations of 47 metabolites [e.g., Octanoylglucuronide, 2’-Deoxyguanosine 5′-monophosphate (dGMP), Adenosine monophosphate (AMP), 2’-Deoxyadenosine 5′-monophosphate (dAMP) and Adenosine] were significantly higher in the rumen of DJC group, and the relative concentrations of 18 metabolites [e.g., 2-Oxoadipic acid, Luteolin, Prostaglandin A2, and Prostaglandin B2 (PGB2) were significantly higher in the DJCF group] (p < 0.05, VIP > 1). A total of 27 metabolites, such as Octanoylglucuronide, Eupatilin, Astragalin, 5-Hydroxyvalproic acid and Kynurenic acid were abundant in the DJC group, and 2-Hydroxybutyric acid, Prostaglandin A2, Prostaglandin B2 (PGB2), 13(S)-HODE, and PC (16:0/P-16:0) were abundant in the DJCY group by the comparison analysis in the DJCY vs. DJC group (Figures 6D–F).
 
  Figure 6. Ruminal differential metabolites identified by volcano plots (A–C) and hierarchical clustering analysis (E,F) in Shorthorn cattle under different feeding practices. Comparisons: (A,D) DJCF vs. DJCY, (B,E) DJCF vs. DJC, (C,F) DJCY vs. DJC.
3.3.3 Metabolic pathways of differential metabolites
Using MetPA, we identified 14 metabolic pathways in the DJCF vs. DJCY comparison; among these, the starch and sucrose metabolism pathway and the purine metabolism pathway were significant. In the DJCY vs. DJC comparison, six pathways were identified, but none reached significance (p < 0.01, pathway impact > 0.1) (Figures 7A–C). According to KEGG mapping, cellobiose, trehalose, maltose, and isomaltose mapped to the starch and sucrose metabolism pathway, whereas adenosine 5′-diphosphate (ADP), adenosine monophosphate (AMP), adenosine, 2′-deoxyadenosine 5′-monophosphate (dAMP), deoxyinosine, deoxyguanosine, and adenine mapped to purine metabolism.
 
  Figure 7. Pathway enrichment analysis using significant differences in rumen metabolites of Shorthorn cattle under different feeding practices. (A) DJCF vs. DJCY, (B) DJCF vs. DJC, (C) DJCY vs. DJC.
3.3.4 Correlations between the ruminal metabolomes and microbiomes
To investigate links between the rumen microbiota and host metabolism, we assessed correlations between the 10 most abundant bacterial genera and the 10 most abundant fungal genera and metabolites detected in rumen fluid using Spearman’s rank correlation analysis (Figures 8A,B). Among bacterial genera, Prevotella showed significantly positive correlations with eight metabolites, including daltogen, cytidine 5′-monophosphate (CMP), N6-acetyllysine, betonicine, and 2′-deoxyguanosine 5′-monophosphate (dGMP), and significantly negative correlations with 20 metabolites, including yohimbine, PC (14:0/P-16:0), sebacic acid, 5-hydroxyindole-3-acetic acid, trimethoprim, Zeran, and adenine (p < 0.05). gF082 showed significant positive correlations with phytanic acid, eupatilin, 2-oxoadipic acid, luteolin, ethyl 4-hydroxybenzoate, and caprolactam (p < 0.05). g_NK4A214_group showed significant positive correlations with six differentially abundant metabolites, and g_[Eubacterium]_coprostanoligenes_group with 27 differentially abundant metabolites (p < 0.05) (Figure 8A).
 
  Figure 8. Correlation analysis between genera and metabolite concentrations affected by the different feeding methods. (A) Correlation between differential bacteria and differential metabolites. (B) Correlation between differential fungi and differential metabolites.
Among fungal genera (Figure 8B), Caecomyces, Neocallimastigaceae gen. Incertae sedis, and Neocallimastix showed positive correlations with 14 metabolites, including eupatilin, PC (22:0/20:2(11Z,14Z)), PC (16:0/P-16:0), astragalin, suberic acid, 5-aminopentanoic acid, S-adenosylmethioninamine, and ethyl 4-hydroxybenzoate. Conversely, Orpinomyces was significantly negatively associated with these differentially abundant metabolites (p < 0.05). Saccharomycetales_gen_Incertae_sedis was positively associated with 24 metabolites (p < 0.05), such as cytidine 5′-monophosphate (CMP), zidovudine, 8-hydroxycarbostyril, trehalose, PC (16:0/16:0), diphenylamine, and turanose.
4 Discussion
The microbial community within the rumen is essential for sustaining the physiological health and normal digestive function of ruminants (Ley et al., 2008). In this study, our multi-omics analysis revealed distinct shifts in rumen microbiota and metabolomes in Shorthorn cattle under different feeding regimens, thereby providing novel insights into the interplay between dietary factors, microbial dynamics, and metabolic outcomes. The taxonomically rich rumen bacterial community drives anaerobic digestion to generate bioavailable substrates for host metabolism (Saleem et al., 2012; Simon et al., 2017). Alpha diversity serves as a key indicator of microbial community robustness, with higher values indicating greater structural stability and functional resilience against environmental perturbations (Jiang F. et al., 2021). Our results showed that Chao1 and Shannon indices were significantly higher in the naturally grazed group (DJCF) than in the intensively fed groups (DJCY, DJC), consistent with previous findings that high-fiber diets enhance community stability. Liu et al. observed that increasing dietary concentrates reduced rumen microbial diversity in yaks (Liu et al., 2019); similarly, He et al. reported higher microbial diversity in naturally grazed cattle (He et al., 2024a). This effect may stem from a decrease in rumen pH due to high-concentrate diets, which in turn inhibit the growth of fiber-degrading bacteria (Sonnenburg and Bäckhed, 2016). The PCoA results further confirmed the distinct bacterial and fungal community structures between DJCF and the intensive-feeding groups, with the first two principal components explaining 24.76 and 19.76% of bacterial variance, respectively, which indicates a strong dietary impact on microbial assembly.
At the phylum level, we found that Bacteroidota and Firmicutes were dominant; together, these phyla have been reported to comprise more than 80% of the bacteria in the gastrointestinal tract (Eckburg et al., 2005; Jiang C. X. et al., 2021; Liu et al., 2021). These observations are consistent with previous research (Chen et al., 2015; Ma et al., 2023; Zhou et al., 2017). Bacteroidota utilize carbohydrate-active enzymes (CAZymes) to ferment dietary polysaccharides into short-chain fatty acids (SCFAs) while mediating protein, bile acid, and oligosaccharide metabolism (He et al., 2024b; Hess et al., 2011; Thomas et al., 2011). Synergistically, Firmicutes employ specialized hydrolases to bioconvert complex substrates (e.g., cellulose, proteins) into bioavailable nutrients (He et al., 2024b; Huo et al., 2014). This cross-phylum partnership optimizes host digestive efficiency through complementary metabolic pathways (He et al., 2024a; Nagaraja and Titgemeyer, 2007). The DJCF group had the highest Bacteroidota abundance (57.62%) and the lowest Firmicutes abundance (34.07%), while the DJC group showed the opposite trend (Bacteroidota 48.84% and Firmicutes 43.08%). This shift coincided with changes in metabolic profiles: higher Bacteroidota in DJCF was associated with elevated fibre-degradation metabolites (e.g., cellobiose), whereas increased Firmicutes in the intensive-feeding groups correlated with higher purine metabolites (AMP, dAMP) and carbohydrate intermediates (maltose, glucose). These findings suggest that the Bacteroidota–Firmicutes balance acts as a key regulator of nutrient-utilization pathways in response to dietary change.
At the genus level, Prevotella and Rikenellaceae_RC9_gut_group were the dominant taxa, both belonging to Bacteroidetes, in agreement with the results of the present study (He et al., 2024a; Yang et al., 2024). Functionally, Prevotella is a key rumen bacterium mediating proteolysis and amylolysis via secreted enzymes (proteases, amylases, hemicellulases), thereby enabling the degradation of dietary proteins, starch, and hemicellulose; its abundance typically increases with higher dietary protein (Bekele et al., 2010; Cholewińska et al., 2020; He et al., 2024a). Rikenellaceae_RC9_gut_group, in turn, has been linked to lipid metabolism and the production of short-chain fatty acids, notably propionate, and may also contribute to protein degradation (Conte et al., 2022; Fan et al., 2017). Diet appears to modulate the relative dominance of these genera: Fernando et al. reported that concentrate-fed beef cows had a higher proportion of Bacteroidetes than those fed green hay, with Rikenellaceae_RC9_gut_group as the dominant genus (Fernando et al., 2010; He et al., 2024a), whereas in our study Prevotella predominated in the DJCY and DJC groups. This discrepancy likely reflects the higher crude-protein and starch contents of the DJCY/DJC total mixed-ration diets, which favor protein- and starch-degrading microbes such as Prevotella. Correlation analyses further supported these functional inferences: Prevotella, most abundant in DJCY (22.52%), showed significant positive correlations with cytidine 5′-monophosphate (CMP) and 2′-deoxyguanosine 5′-monophosphate (dGMP), and negative correlations with phytanic acid and caprolactam, suggesting a context-specific role in nucleotide metabolism under intensive feeding. By contrast, Rikenellaceae_RC9_gut_group, most abundant in DJCF (12.56%), associated with higher levels of propionate precursors, consistent with its contribution to lipid metabolism and short-chain fatty acid (SCFA) synthesis (Conte et al., 2022). Taken together, these patterns align with evidence that altering feed type and concentrate proportion (e.g., TMR with higher crude protein and starch) reshapes the Bacteroidetes/Firmicutes balance and modulates carbohydrate and purine-related pathways, thereby favoring Prevotella-driven metabolism under intensive diets.
Anaerobic fungi are pivotal decomposers in the rumen, and Neocallimastigomycota typically predominate; in our dataset they accounted for 75.68–86.74% of the fungal community, in line with reports that this phylum is the dominant fungal group in ruminants (Boots et al., 2013; Denman et al., 2018; Wallace et al., 2019). Diet remodels fungal diversity and composition: intensive feeding reshapes the rumen microbiota, including the mycobiome, whereas seasonal and forage quality shifts also modulate fungal taxa (Cengiz and John, 2002; Faichney et al., 1997; He et al., 2024b; Wang et al., 2023). In our study, the natural-grazing group (DJCF) displayed higher Neocallimastigomycota abundance than the intensive-feeding groups, consistent with literature showing context-dependent fungal responses to dietary inputs in high-altitude ruminants. Functionally, Neocallimastigomycota enhance nutrient bioavailability by secreting potent cellulases and cellulosomes that depolymerize recalcitrant plant cell walls (e.g., Orpinomyces systems) (Abhaya et al., 2017; Boots et al., 2013). Beyond classic fiber degradation, our integrative analyses revealed fungus–metabolite couplings that highlight ecological specialization: Orpinomyces (enriched in DJCF) correlated positively with 2-oxoadipic acid and luteolin—metabolites linked to plant-matrix breakdown—whereas Caecomyces (elevated under intensive feeding) correlated with phosphatidylcholine (PC) species, suggesting adjustment of lipid-related metabolism under high-concentrate diets. These findings extend current evidence that feed type and concentrate proportion restructure the rumen community and its metabolic outputs, and they pinpoint fungal lineages that track specific metabolite pathways, thereby underscoring fungi’s contribution to host energy harvest beyond lignocellulose turnover (Denman et al., 2018; He et al., 2024a; He et al., 2024b; Wallace et al., 2019).
Intensive or high-concentrate feeding reshapes rumen carbohydrate metabolism and fermentation; in yaks, multi-omics shows that intensive systems elevate dietary starch/protein and reconfigure the microbiome–metabolome interface (Yang et al., 2024). Consistent with this paradigm, our intensive-feeding groups (DJCY/DJC) accumulated carbohydrate intermediates (e.g., maltose, glucose), and pathway analysis indicated significant shifts in starch and sucrose metabolism together with purine metabolism, a pattern widely observed when feed type or forage-to-concentrate ratios change (Liu et al., 2019; Penner et al., 2009). Similar to the results of this experiment, we observed a significant increase of carbohydrate metabolites in the DJCY group, such as maltose, glucose, which are important intermediates of carbohydrate metabolism. The metabolic pathway is primarily sucrose and starch metabolism. In the rumen, microbial enzymes break down starch and cellulose into glucose, which then undergoes glycolysis to produce pyruvate. This pyruvate, derived from carbohydrate metabolites, acts as the key substrate for VFA production (Xue et al., 2018). Metabolic pathway analysis revealed that sucrose and starch metabolism and purine metabolism were significantly affected by feeding regimens. Four metabolites (Cellobiose, Trehalose, Maltose, Isomaltose) in the sucrose and starch pathway were more abundant in DJCF, aligning with the higher fiber-degrading microbiota in this group. Meanwhile, purine metabolites (AMP, dAMP, Adenosine) were enriched in DJCY and DJC, which corresponded to the elevated Prevotella abundance and were consistent with the role of these metabolites as biomarkers for microbial protein synthesis (Ametaj et al., 2010). The significant correlation between Prevotella and purine metabolites further confirms the direct involvement of specific taxa in regulating these pathways. Notably, our study identified several low-abundance microbial taxa with potential functional significance. LEfSe analysis revealed nine differential species in DJCF, including Fibrobacter succinogenes, which showed positive correlations with phytanic acid and eupatilin despite its low relative abundance. Similarly, the DJC group had specific low-abundance oscillospiraceae taxa associated with octanoylglucuronide, a metabolite involved in lipid conjugation. These findings highlight that even low-abundance microbes may play critical roles in metabolic regulation, but their actual metabolic functions remain unclear due to limitations in current multi-omics approaches.
Of particular interest is the discovery of changes in several metabolites associated with the purine metabolic pathway. Ruminal microbes degrade nitrogenous compounds from feed and channel them into nucleic-acid metabolism through de novo and salvage pathways; nucleosides are then catabolized to bases (e.g., xanthine, hypoxanthine, guanine, adenine) and further to degradation products (Fujihara and Shem, 2011; Stentoft et al., 2015). We found higher levels of purine-related metabolites in the concentrate-fed groups, consistent with previous reports of elevated xanthine and hypoxanthine under high-concentrate diets (Ametaj et al., 2010; Zhang et al., 2017). In addition, adenine and xanthine—and likely pseudouridine—identified among the feed-intake groups likely originate from bacterial nucleic-acid degradation (Noel et al., 2017; Sutton et al., 1975). In the above results, it was shown that concentrate feed contributes to rumen purine metabolism. Arachidonic acid functions as a sensitive and specific plasma biomarker for the average daily gain (ADG) in steers (Artegoitia et al., 2016). Furthermore, Naturally occurring in all mammalian cells, arachidonic acid (an ω-6 tetra-unsaturated fatty acid) maintains membrane fluidity while serving as the primary biosynthetic precursor for lipid-derived signaling molecules, including prostaglandins and leukotrienes (Martin et al., 2016). In this study, we found a higher concentration of prostaglandin in the DJCF, DJCY group compared with that in the DJC group. It may be related to the effect of castration of Shorthorn cattle on androgen. Phosphatidylcholine (PC), a glycerophospholipid, serves as a critical structural element in cellular lipid bilayers. Choline constitutes a primary component of PC, and adequate PC levels are essential for supporting lipid transport mechanisms (Shiau and Cho, 2002; Yeh et al., 2015). PC functions as a choline reservoir, with the stored choline being further processed into organic osmolytes through metabolic pathways involving glycine, serine, and threonine (Athamena et al., 2011; Jiang et al., 2018). Ren et al. showed that PC concentration increased linearly with increasing concentrate levels and that this metabolite, which is related to lipid metabolism, plays an important role in yak growth and rumen health. The results of the present experimental study showed that the PC concentration in the two fattening groups were significantly higher than that in the grazing group, indicating that the PC concentration increased with the nutrient level of the diets. To further explore the potential connection between rumen dominant microflora, and metabolites, Pearson correlation analysis was carried out. Dai et al. indicated that differential metabolites are significantly correlated, either positively or negatively, with dominant microorganisms (Dai et al., 2023). Emerging evidence from integrated multi-omics investigations has revealed utilization or production associations between bacteria and metabolomes in the rumen (Liu et al., 2019; Ma et al., 2021). Saccharomycetales_gen_Incertae_sedis in the DJCY group showed positive correlations with 24 metabolites (including CMP and turanose), implying a potential role in nucleotide and carbohydrate metabolism. However, the regulatory mechanisms driving these associations—whether through direct metabolic conversion or indirect ecological interactions—remain to be elucidated.
Overall, these patterns highlight the effects of different feeding strategies in Shorthorn cattle and provide insights into the mechanisms by which feeding practices shape rumen microbial processes. While our study integrates 16S rRNA gene and ITS amplicon sequencing with metabolomics, directly inferring in situ metabolic functions and the specific regulatory roles of low-abundance bacteria and fungi remains challenging. In future work, we will combine metagenomic and metatranscriptomic approaches to delineate links between microbial gene expression and metabolic pathways, and we will elucidate organism-level functions in the rumen microecosystem through isolation and in-vitro metabolic validation of key strains, thereby providing a theoretical basis for precision nutrition in Shorthorn cattle.
5 Conclusion
This study used multi-omics analysis to reveal the adaptive mechanisms of the rumen microbiome and metabolites of Shorthorn cattle in intensive feeding systems. High-concentrate diets enriched proteolytic bacteria Prevotella and suppressed fiber-degrading fungi Neocallimastigomycota, activating purine metabolism pathways, indicating enhanced microbial protein synthesis efficiency. Concurrently, Shorthorn cattle maintained Bacteroidota functionality under high-starch diets through host genetic adaptation. Thus, an intensive feeding system altered the rumen microbiome and rumen metabolism, resulting in improvements of Shorthorn cattle growth. These findings provide novel strategies for precision feeding based on metabolic biomarkers, though further validation is required to address limitations such as limited sample size and potential confounding effects of castration.
Data availability statement
The datasets presented in this study can be found in online repositories. This data can be found here: https://www.ncbi.nlm.nih.gov/, PRJNA1277806.
Ethics statement
The animal study was approved by Animal Protection and Utilization Committee of Yunnan Agricultural University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
SD: Writing – original draft. SH: Writing – original draft, Data curation. SZ: Writing – original draft, Visualization. QL: Writing – review & editing, Investigation. HM: Formal analysis, Methodology, Writing – review & editing. DW: Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
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.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.
Footnotes
References
Abhaya, R., Smith, O. P., Youssef, N. H., Struchtemeyer, C. G., Atiyeh, H. K., and Elshahed, M. S. (2017). Utilizing anaerobic fungi for two-stage sugar extraction and biofuel production from lignocellulosic biomass. Front. Microbiol. 8:635. doi: 10.3389/fmicb.2017.00635
Ahmad, A. A., Yang, C., Zhang, J., Kalwar, Q., and Ding, X. (2020). Effects of dietary energy levels on rumen fermentation, microbial diversity, and feed efficiency of yaks (Bos grunniens). Front. Microbiol. 11:625. doi: 10.3389/fmicb.2020.00625
Ametaj, B. N., Zebeli, Q., Saleem, F., Psychogios, N., and Wishart, D. S. (2010). Metabolomics reveals unhealthy alterations in rumen metabolism with increased proportion of cereal grain in the diet of dairy cows. Metabolomics 6, 583–594. doi: 10.1007/s11306-010-0227-6
Artegoitia, V. M., Foote, A. P., Lewis, R. M., and Freetly, H. C. (2016). 1605 rumen fluid metabolomics analysis associated with feed efficiency on crossbred steers. J. Anim. Sci. 94, 780–781. doi: 10.2527/jam2016-1605
Athamena, A., Brichon, G., Trajkovic-Bodennec, S., Péqueux, A., Chapelle, S., Bodennec, J., et al. (2011). Salinity regulates N-methylation of phosphatidylethanolamine in euryhaline crustaceans hepatopancreas and exchange of newly-formed phosphatidylcholine with hemolymph. J. Comp. Physiol. B 181, 731–740. doi: 10.1007/s00360-011-0562-6
Bekele, A. Z., Koike, S., and Kobayashi, Y. (2010). Genetic diversity and diet specificity of ruminal Prevotella revealed by 16S rrna gene-based analysis. FEMS Microbiol. Lett. 305, 49–57. doi: 10.1111/j.1574-6968.2010.01911.x
Boots, B., Lillis, L., Clipson, N., Petrie, K., Kenny, D. A., Boland, T. M., et al. (2013). Responses of anaerobic rumen fungal diversity (phylum Neocallimastigomycota) to changes in bovine diet. J. Appl. Microbiol. 114, 626–635. doi: 10.1111/jam.12067
Callahan, B. J., Wong, J., Heiner, C., Oh, S., Theriot, C. M., Gulati, A. S., et al. (2019). High-throughput amplicon sequencing of the full-length 16S rrna gene with single-nucleotide resolution. Nucleic Acids Res. 47:e103. doi: 10.1093/nar/gkz569
Cengiz, A., and John, W. R. (2002). De novo synthesis of amino acids by the ruminal anaerobic fungi, Piromyces communis and Neocallimastix frontalis. FEMS Microbiol. Lett. 2, 243–247. doi: 10.1016/S0378-1097(02)00752-8
Chen, Y. B., Lan, D. L., Tang, C., Yang, X. N., and Li, J. (2015). Effect of Dna extraction methods on the apparent structure of yak rumen microbial communities as revealed by 16S rdna sequencing. Pol. J. Microbiol. 64, 29–36. doi: 10.33073/pjm-2015-004
Cholewińska, P., CzyŻ, K., Nowakowski, P., and Wyrostek, A. (2020). The microbiome of the digestive system of ruminants–a review. Anim. Health Res. Rev. 21, 3–14. doi: 10.1017/S1466252319000069
Conte, G., Dimauro, C., Daghio, M., Serra, A., Mannelli, F., Mcammond, B., et al. (2022). Exploring the relationship between bacterial genera and lipid metabolism in bovine rumen. Animal 16:100520. doi: 10.1016/j.animal.2022.100520
Dai, R., Ma, X., Dingkao, R., Huang, C., La, Y., Li, X., et al. (2023). Effects of dietary crude protein levels in the concentrate supplement after grazing on rumen microbiota and metabolites by using metagenomics and metabolomics in Jersey-yak. Front. Microbiol. 14:1124917. doi: 10.3389/fmicb.2023.1124917
Deerry, M. (2018). Genetics, biotechnology, and breeding: north American shorthorn production in the twenty-first century. Agric. Hist. 92:54. doi: 10.3098/ah.2018.092.1.054
Denman, S. E., Morgavi, D. P., and Mcsweeney, C. S. (2018). Review: the application of omics to rumen microbiota function. Animal 12, s233–s245. doi: 10.1017/S175173111800229X
Eckburg, P. B., Bik, E. M., Bernstein, C. N., Purdom, E., Dethlefsen, L., Sargent, M., et al. (2005). Diversity of the human intestinal microbial Flora. Science 308, 1635–1638. doi: 10.1126/science.1110591
Faichney, G., Poncet, C., Lassalas, B., Jouany, J., Millet, L., Doré, J., et al. (1997). Effect of concentrates in a hay diet on the contribution of anaerobic fungi, protozoa and bacteria to nitrogen in rumen and duodenal digesta in sheep. Anim. Feed Sci. Technol. 64, 193–213. doi: 10.1016/s0377-8401(96)01059-0
Fan, P., Liu, P., Song, P., Chen, X., and Ma, X. (2017). Moderate dietary protein restriction alters the composition of gut microbiota and improves ileal barrier function in adult pig model. Sci. Rep. 7:43412. doi: 10.1038/srep43412
Fernando, S. C., Purvis, H. T., Najar, F. Z., Sukharnikov, L. O., Krehbiel, C. R., Nagaraja, T. G., et al. (2010). Rumen microbial population dynamics during adaptation to a high-grain diet. Appl. Environ. Microbiol. 76, 7482–7490. doi: 10.1128/Aem.00388-10
Fujihara, T., and Shem, M. N. (2011). Metabolism of microbial nitrogen in ruminants with special reference to nucleic acids. Anim. Sci. J. 82, 198–208. doi: 10.1111/j.1740-0929.2010.00871.x
He, S. C., Yuan, Z. M., Dai, S. F., Wang, Z. B., Zhao, S. S., Wang, R. J., et al. (2024a). Intensive feeding alters the rumen microbiota and its fermentation parameters in natural grazing yaks. Front. Vet. Sci. 11:1365300. doi: 10.3389/fvets.2024.1365300
He, S. C., Zhao, S. S., Wang, Z. B., Dai, S. F., Mao, H. M., and Wu, D. W. (2024b). Impact of seasonal variation in pasture on rumen microbial community and volatile fatty acids in grazing yaks: insights from high-altitude environments. Microorganisms 12:1701. doi: 10.3390/microorganisms12081701
Hess, M., Sczyrba, A., Egan, R., Kim, T. W., Chokhawala, H., Schroth, G., et al. (2011). Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science 331, 463–467. doi: 10.1126/science.1200387
Hu, C. S., Ding, L. M., Jiang, C. X., Ma, C. F., Liu, B. T., Li, D. L., et al. (2021). Effects of management, dietary intake, and genotype on rumen morphology, fermentation, and microbiota, and on meat quality in yaks and cattle. Front. Nutr. 8:755255. doi: 10.3389/fnut.2021.755255
Huo, W. J., Zhu, W. Y., and Mao, S. Y. (2014). Impact of subacute ruminal acidosis on the diversity of liquid and solid-associated bacteria in the rumen of goats. World J. Microbiol. Biotechnol. 30, 669–680. doi: 10.1007/s11274-013-1489-8
Jiang, C. X., Ding, L. M., Dong, Q. M., Wang, X. J., Wei, H. Y., Hu, C. S., et al. (2021). Effects of root extracts of three traditional Chinese herbs as dietary supplements on dry matter intake, average daily gain, rumen fermentation and ruminal microbiota in early weaned yak calves. Anim. Feed Sci. Technol. 278:115002. doi: 10.1016/j.anifeedsci.2021.115002
Jiang, F., Gao, H. M., Qin, W., Song, P. F., Wang, H. J., Zhang, J. J., et al. (2021). Marked seasonal variation in structure and function of gut microbiota in forest and alpine musk deer. Front. Microbiol. 12:2460. doi: 10.3389/fmicb.2021.699797
Jiang, W. W., Tian, X. L., Fang, Z. H., Li, L., Dong, S. L., Li, H. D., et al. (2018). Metabolic responses in the gills of tongue sole (Cynoglossus semilaevis) exposed to salinity stress using NMR-based metabolomics. Sci. Total Environ. 653, 465–474. doi: 10.1016/j.scitotenv.2018.10.404
Knight, R. (2011). Uchime improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Ley, R. E., Lozupone, C. A., Hamady, M., Knight, R., and Gordon, J. I. (2008). Worlds within worlds: evolution of the vertebrate gut microbiota. Nat. Rev. Microbiol. 6, 776–788. doi: 10.1038/nrmicro1978
Lin, L. M., Xie, F., Sun, D. M., Liu, J. H., Zhu, W. Y., and Mao, S. Y. (2019). Ruminal microbiome-host crosstalk stimulates the development of the ruminal epithelium in a lamb model. Microbiome 7:83. doi: 10.1186/s40168-019-0701-y
Liu, H., Jiang, H., Hao, L. Z., Cao, X. L., Degen, A., and Zhou, J. W. (2021). Rumen bacterial community of grazing lactating yaks (Poephagus grunniens) supplemented with concentrate feed and/or rumen-protected lysine and methionine. Animals 11:2425. doi: 10.3390/ani11082425
Liu, C., Wu, H., Liu, S. J., Chai, S. T., Meng, Q. X., and Zhou, Z. M. (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
Ma, X. Y., La, Y. F., Yang, G. W., Dai, R. F., Zhang, J. X., and Zhang, Y. H. (2023). Multi-omics revealed the effects of dietary energy levels on the rumen microbiota and metabolites in yaks under house-feeding conditions. Front. Microbiol. 14:1309535. doi: 10.3389/fmicb.2023.1309535
Ma, Y., Wang, C., Zhang, H., Yu, L. H., Dong, L., and Gong, D. Q. (2021). Illumina sequencing and metabolomics analysis reveal thiamine modulation of ruminal microbiota and metabolome characteristics in goats fed a high-concentrate diet. Front. Microbiol. 12:653283. doi: 10.3389/fmicb.2021.653283
Martin, S. A., Brash, A. R., and Murphy, R. C. (2016). The discovery and early structural studies of arachidonic acid. J. Lipid Res. 57:jlr.R068072. doi: 10.1194/jlr.R068072
Moorby, J. M., and Fraser, M. D. (2021). Review: New feeds and new feeding systems in intensive and semi-intensive forage-fed ruminant livestock systems. Animal. 15:100297. doi: 10.1016/j.animal.2021.100297
Nagaraja, T. G., and Titgemeyer, E. C. (2007). Ruminal acidosis in beef cattle: the current microbiological and nutritional outlook. J. Dairy Sci. 90, E17–E38. doi: 10.3168/jds.2006-478
Noel, S. J., Attwood, G. T., Rakonjac, J., Moon, C. D., and Janssen, P. H. (2017). Seasonal changes in the digesta-adherent rumen bacterial communities of dairy cattle grazing pasture. PLoS One 12:e0173819. doi: 10.1371/journal.pone.0173817
Pang, K. Y., Chai, S. T., Yang, Y. K., Wang, X., Liu, S. J., and Wang, S. X. (2022). Dietary forage to concentrate ratios impact on yak ruminal microbiota and metabolites. Front. Microbiol. 13:964564. doi: 10.3389/fmicb.2022.964564
Penner, G. B., Taniguchi, M., Guan, L. L., Beauchemin, K. A., and Oba, M. (2009). Effect of dietary forage to concentrate ratio on volatile fatty acid absorption and the expression of genes related to volatile fatty acid absorption and metabolism in ruminal tissue. J. Dairy Sci. 92, 2767–2781. doi: 10.3168/jds.2008-1716
Saleem, F., Bouatra, S., Guo, A. C., Psychogios, N., and Wishart, D. S. (2012). The bovine ruminal fluid metabolome. Metabolomics 9, 360–378. doi: 10.1007/s11306-012-0458-9
Shiau, S. Y., and Cho, W. H. (2002). Choline requirements of grass shrimp (Penaeus monodon) as affected by dietary lipid level. Anim. Sci. 75, 97–102. doi: 10.1017/S1357729800052875
Simon, D., Amélia, C. S., Jürgen, C., Uwe, B., Markus, R., Jana, S., et al. (2017). A structural and functional elucidation of the rumen microbiome influenced by various diets and microenvironments. Front. Microbiol. 8:1605. doi: 10.3389/fmicb.2017.01605
Smoot, M. E., Ono, K., Ruscheinski,, Wang, P. L., and Ideker, T. (2011). Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics 27, 431–432. doi: 10.1093/bioinformatics/btq67
Sonnenburg, J. L., and Bäckhed, F. (2016). Diet–microbiota interactions as moderators of human metabolism. Nature 535, 56–64. doi: 10.1038/nature18846
Stentoft, C., Røjen, B. A., Jensen, S. K., Kristensen, N. B., and Vestergaard, M. (2015). Absorption and intermediary metabolism of purines and pyrimidines in lactating dairy cows. Briti. J. Nutr. 113, 560–573. doi: 10.1017/S0007114514004000
Sutton, J. D., Smith, R. H., Mcallan, A. B., Storry, J. E., and Corse, D. A. (1975). Effect of variations in dietary protein and of supplements of cod-liver oil on energy digestion and microbial synthesis in the rumen of sheep fed hay and concentrates. J. Agric. Sci. 84, 317–326. doi: 10.1017/S0021859600052461
Thomas, F., Hehemann, J. H., Rebuffet, E., Czjzek, M., and Michel, G. (2011). Environmental and gut Bacteroidetes: the food connection. Front. Microbiol. 2:93. doi: 10.3389/fmicb.2011.00093
Wallace, R. J., Sasson, G., Garnsworthy, P. C., Tapio, I., and Mizrahi, I. (2019). A heritable subset of the core rumen microbiome dictates dairy cow productivity and emissions. Sci. Adv. 5:eaav8391. doi: 10.1126/sciadv.aav8391
Wang, Y. Y., Guo, H., Gao, X. G., and Wang, J. H. (2021). The intratumor microbiota signatures associate with subtype, tumor stage, and survival status of esophageal carcinoma. Front. Oncol. 11:754788. doi: 10.3389/fonc.2021.754788
Wang, R. J., He, S. C., Huang, D., He, S. Y., Guo, T. Q., Chen, T., et al. (2023). Differences in composition and diversity of rumen fungi in buffalo fed different diets. Anim. Biotechnol. 34, 5075–5086. doi: 10.1080/10495398.2023.2276974
Xue, F. G., Pan, X. H., Jiang, L. S., Guo, Y. M., and Xiong, B. H. (2018). GC–MS analysis of the ruminal metabolome response to thiamine supplementation during high grain feeding in dairy cows. Metabolomics 14, 1–12. doi: 10.1007/s11306-018-1362-8
Yang, S. L., Zheng, J. Y., Mao, H. M., Vinitchaikul, P., Wu, D. W., and Chai, J. M. (2024). Multiomics of yaks reveals significant contribution of microbiome into host metabolism. NPJ Biofilms Microbiomes 10:133. doi: 10.1038/s41522-024-00609-2
Yeh, S. P., Shiu, P. J., Guei, W. C., Lin, Y. H., and Liu, C. H. (2015). Improvement in lipid metabolism and stress tolerance of juvenile giant grouper, Epinephelus lanceolatus (Bloch), fed supplemental choline. Aqua C. Res. 46, 1810–1821. doi: 10.1111/are.12334
Zhang, T., Mu, Y., Zhang, D., Lin, X., Wang, Z. H., Hou, Q. L., et al. (2018). Determination of microbiological characteristics in the digestive tract of different ruminant species. Microbiology 8:e00769. doi: 10.1002/mbo3.769
Zhang, R. Y., Zhu, W. Y., Jiang, L. S., and Mao, S. Y. (2017). Comparative metabolome analysis of ruminal changes in Holstein dairy cows fed low- or high-concentrate diets. Metabolomics 13:74. doi: 10.1007/s11306-017-1204-0
Keywords: Shorthorn cattle, rumen microbial community, metabolomics, different feeding methods, 16S rRNA/ITS, multi-omics
Citation: Dai S, He S, Zhao S, Li Q, Mao H and Wu D (2025) Feeding regimens reshape rumen microbiota and metabolome in Shorthorn cattle: a multi-omic insight into microbial diversity and metabolic pathway dynamics. Front. Microbiol. 16:1657402. doi: 10.3389/fmicb.2025.1657402
Edited by:
P. K. Malik, National Institute of Animal Nutrition and Physiology (ICAR), IndiaReviewed by:
Shraddha Trivedi, National Institute of Animal Nutrition and Physiology (ICAR), IndiaBrijesh Yadav, U.P. Pandit Deen Dayal Upadhyaya Veterinary Science University, India
Copyright © 2025 Dai, He, Zhao, Li, Mao and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: DongWang Wu, ZGFud2V5QDE2My5jb20=; HuaMing Mao, MTk4NTAzOUB5bmF1LmVkdS5jbg==
 ShuSheng Zhao
ShuSheng Zhao