Next Article in Journal
A Survey of Archaeal Restriction–Modification Systems
Previous Article in Journal
SARS-CoV-2 Rapid Antigen Test Based on a New Anti-Nucleocapsid Protein Monoclonal Antibody: Development and Real-Time Validation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 11
Issue 10
10.3390/microorganisms11102423
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Preliminary Study of the Potential Molecular Mechanisms of Individual Growth and Rumen Development in Calves with Different Feeding Patterns

by
Jie Wang
1,†,
Kaisen Zhao
2,†,
Mianying Li
2,
Huimei Fan
2,
Meigui Wang
2,
Siqi Xia
2,
Yang Chen
2,
Xue Bai
2,
Zheliang Liu
2,
Jiale Ni
2,
Wenqiang Sun
1,
Xianbo Jia
1 and
Songjia Lai
1,*
1
Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China
2
College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Microorganisms 2023, 11(10), 2423; https://doi.org/10.3390/microorganisms11102423
Submission received: 29 August 2023 / Revised: 19 September 2023 / Accepted: 22 September 2023 / Published: 28 September 2023
(This article belongs to the Section Veterinary Microbiology)
Browse Figures
Review Reports Versions Notes

Abstract

:
At present, it is common to feed calves with “Concentrate”, “Concentrate + hay” and TMR “Total Mixed Rations” feeding patterns in China, which achieved well feeding efficiency, but the three feeding patterns molecular regulation mechanism in actual production is still unclear. The study aimed to explore the most suitable feeding pattern for Chinese Holstein calves to improve the rumen fermentation function and growth performance of calves. In this regard, the interactions between rumen microorganisms and host metabolism were investigated. The rumen volume and weight of calves in the GF group were significantly higher than those in the GFF and TMR groups (p < 0.05), and the rumen pH of calves in the GF group was 6.47~6.79. Metagenomics analysis revealed that the rumen microbiome of GF and GFF calves had higher relative abundances of Methanobrevibacter, Methanosphaera, and Methanolacinia (p < 0.05). Prevotella multisaccharivorax was significantly more abundant in the rumen of GF calves (p < 0.05), indicating that GF group calves had a stronger ability to ferment sugars. Notably, in the pyruvate metabolic pathway, phosphoenolpyruvate carboxylase was significantly up-regulated in GF calves compared with the TMR group, and pyruvate-phosphate dikinase was significantly down-regulated. Metabolomic results showed that Ursodeoxycholic acid was significantly up-regulated in GF calves, and most of the differential metabolites were enriched in Bile secretion pathways. The association analysis study found that the microorganisms of Prevotella and Ruminococcaceae might cooperate with the host, which was helpful for the digestion and absorption of lipids and made the calves have better growth. The three feeding modes had similar effects, but the ‘GF’ feeding pattern was more beneficial to the individual growth and ruminal development regarding ruminal morphology, contents physiology and microorganisms. Furthermore, the synergistic effect of rumen microorganisms and the host could more effectively hydrolyze lipid substances and promote the absorption of lipids, which was of great significance to the growth of calves.

1. Introduction

The health and biological functioning of livestock are often prioritized. In recent years, nutritional strategies have emerged. It has been proposed as a critical factor in improving the health status and welfare of animals as well as enhancing productivity and performance in livestock [1,2], and it has been well demonstrated that nutritional strategies could increase the level of production and improve the health of animals and products obtained from them. Feed type, feeding method and nutrition level are related to the rumen development of calves [3], and the growth and development of calves during lactation determine the performance of final milk production to a certain extent [4]. At present, the feeding patterns of calves during lactation are different in China, mainly including “concentrate”, “TMR” and “concentrate plus hay” feeding patterns. Feeding only concentrate might occasionally lead to a lower pH value of ruminal liquid, aggregation of rumen papilla, incomplete keratinization, and other symptoms in calves [5]. A few papers reported that due to the high moisture content of TMR, it is not easy to achieve the ideal growth effect, and it is easy to cause calves diarrhea [6]. Premature hay intake by calves might lead to insufficient energy intake of animals and diseases (such as hay belly) due to their rumen not fully developed and their ability to digest roughage was limited. Too late intake of hay might affect the digestion ability of calves to roughage in adulthood, which was not conducive to rumen development and optimal production performance [7,8,9]. However, the molecular mechanism of individual growth and rumen development of calves with three feeding patterns remains unclear.
Microorganisms in the rumen interact with the host, regulate the nutrient metabolism of the organism, participate in the material transport process, and play an important role in the organism’s immune media, which could feedback affect the feed utilization rate and the growth performance of the host [10,11]. The research showed that diet was used as an exogenous dietary signal to mediate the interaction between host and rumen microbiota [12]. It is believed that diet is the main factor affecting the composition of rumen microbial structure [13], and the change in feed nutrient level might cause an adaptive change in the rumen microbial community. Feeding beef cattle with a high concentrate-to-feed ratio significantly increases the number of Streptococcus bovis, Megasphaera elsdenii and Prevotella bryantii and reduces the relative abundance of cellulolytic bacteria [14]. Microbes in the rumen of ruminants mainly include Ruminococcus, Fibrobacteres, Prevotellariae, etc., and they are closely related to the host organism’s production performance and health status of the host organism [15].
Additionally, short-chain fatty acids (SCFAs), metabolites of rumen microorganisms, are the main energy source, have anti-inflammatory and immunomodulatory effects, and could strengthen the intestinal barrier [16]. In recent years, more and more studies have explored the metabolic differences of calves based on metabolomics [17,18,19]. A few papers identified serum metabolites associated with calf nutrition and healthy development [20].
The study aimed to explore the effects of three different feeding patterns of “GF”, “GFF”, and “TMR” on the growth, rumen development, rumen fluid short-chain fatty acids, rumen microbes and organism metabolism of calves, and finally, select the most suitable feeding pattern for calves, reveal their mechanisms at the molecular level, which to lay a theoretical foundation for further explaining the structure and function of rumen microorganism and its degradation mechanism to feed, and provide a theoretical basis for calf feeding.

2. Materials and Methods

2.1. Ethics Statement

The study was approved and conducted in accordance with the ethical standards of the Institutional Animal Care and Use Committee of the College of Animal Science and Technology, Sichuan Agricultural University, Sichuan, China.

2.2. Animals

The study selected nine 7-day-old Chinese Holstein bull calves with healthy, and the average body weight was 41.58 ± 0.79 kg. The nine calves were randomly divided into three groups of GFF, GF and TMR, which, following the completely random grouping principles, feed and water were provided ad libitum, and each individual was raised in a single pen and fed milk. Based on feeding milk, the GFF group was fed with a “ concentrate” pattern. In the GF group, 10% hay was added based on fine feed, in which the mixture ratio of hay was alfalfa grass: oat grass = 3:2, cut into 1–2 cm long and mixed evenly. The TMR group was fed with a “TMR” feed pattern. Other feeding and management methods are carried out according to the existing methods of Sichuan Xuebao Dairy Group Co., Ltd., and the body weight of the calves at 80 days old and slaughter them. The formula and nutritional composition of the diet are shown in Table S1.

2.3. Sample Collection

The sampling time for calves was from 10:00 to 11:00 a.m. Calves fasted for 24 h, jugular vein blood collection and timely centrifugation (4 °C, 3000 r/min, 5 min) to obtain serum. After slaughter, the rumen contents of calves were quickly collected and transported in a 2 mL cryopreservation tube and placed in liquid nitrogen. Finally, they were stored in a refrigerator at −80 °C for subsequent experiments, and all preliminary sample processing was completed within seven days after sample collection. At the same time, the pH value of rumen contents was measured with a PHS-25C PH tester (Ostell Industrial Co., Ltd., Shanghai, China), and rumen weight and volume were determined.

2.4. DNA Extraction, Metagenomic Sequencing

According to the kit instructions, microbial DNA was extracted from rumen content samples (the samples were liquid) using the Tiangen Magnetic Bead Kit (Tiangen Biotech, Beijing, China). The purity and integrity of DNA were detected by 1% agarose gel electrophoresis. DNA was quantified using Qubit® dsDNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies, Foster City, CA, USA), and appropriate samples were taken into centrifuge tubes, dilute the sample with sterile water to an OD value between 1.8~2.0. Then, taking 1 μg of genomic DNA from the sample, a metagenomic library was constructed using NEBNext® Ultra DNA Library Prep Kit for Illumina (NEB, Illumina Inc., San Diego, CA, USA) and sequenced using the Illumina HiSeq Xten platform. The library construction and sequencing work were performed by Beijing Novogene Technology Co., Ltd. Company (Beijing, China). Preprocessing the Raw Data obtained from the Illumina HiSeq sequencing platform using Readfq (V8, https://github.com/cjfields/readfq, accessed on March 2022) to remove reads with low-quality bases (the default quality threshold < 38) exceeding 40 bp, N bases reaching 10 bp, and the overlap between adapters exceeding 15 bp. Then, using bowtie2 software (version 2.2.4, http://bowtie-bio.sourceforge.net/bowtie2/index.shtml, accessed on March 2022) to filter out DNA interference that may originate from calves [21,22]. Clean Data for each sample was assembled and analyzed using MEGAHIT software [23,24] (v1.0.4-beta, https://github.com/voutcn/megahit, accessed on March 2022). MetaGeneMark software [21] (V2.10, http://topaz.gatech.edu/GeneMark/, accessed on March 2022) was used to predict the open reading region (ORF) of scaftigs ≥ 500 bp in each sample, and the information with length less than 100 nt in the prediction results were filtered [25,26]. Finally, the CD-HIT software [27,28] (V4.5.8, http://www.bioinformatics.org/cd-hit, accessed on March 2022) was used to de-redundant the ORF prediction results to obtain a non-redundant initial gene catalogue, using Bowtie2 (Bowtie2.2.4) Compare the Clean Data of each sample to the initial gene catalogue, calculate the number of reads of the genes aligned in each sample [29,30], and obtain the Unigenes for subsequent analysis [31].

2.5. Classification and Functional Annotation of Rumen Metagenomes

According to NCBI’s NR database (Version 2018-01-02, https://www.ncbi.nlm.nih.gov/, accessed on April 2022), using DIAMOND software [32] (V0.9.9.110, https://github.com/bbuchfink/diamond/, accessed on April 2022)for taxonomic assessment of rumen microbes in calves. The LCA algorithm (Applied in Megan software for systematic classification, version 4.0, https://en.wikipedia.org/wiki/Lowest_common_ancestor, accessed on April 2022) was used to determine the species annotation information of the aligned sequences and to obtain the abundance information of each sample at the Kingdom, phylum, genus, and species. Starting from the abundances at each taxonomic level, PCA [33] (RADE4 package, version 2.15.3) dimension reduction analysis was performed, followed by lefse analysis using lefse software [34] (v1.0, LDA score default to 3) to look for differential species between groups.
Unigenes were compared with the KEGG [35] database (version 2018-01-01, http://www.kegg.jp/kegg/, accessed on April 2022) and CAZy [36] databases (version 201801, http://www.cazy.org/, accessed on April 2022) by using DIAMOND software (v0.9.9.110, https://github.com/bbuchfink/diamond/, accessed on April 2022), and alignments were performed with parameters set to blastp, -e 1e-5 [29]. Then, the abundances at different functional levels were counted, starting from the abundances at each taxonomic level, annotated gene number statistics, relative abundance profile display, metabolic pathway comparison analysis, and component difference analysis using lefse.

2.6. Rumen Fluid SCFA and Serum Metabolomic Analysis

The extraction of short-chain fatty acids (SCFA) in rumen fluid and serum metabolites, data preprocessing, and quantitative and qualitative analysis of metabolites were completed by Nuohe Zhiyuan Technology Co., Ltd. (Beijing, China). The concentration of SCFA in the rumen was detected and calculated by the standard curve of acetic acid, propionic acid, butyric acid, isobutyric acid, valeric acid, isovaleric acid and caproic acid. Then, metaX [37] was used to convert the data, principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were performed, and the fold change (FC value) of metabolites between the two groups was calculated. The default criteria for screening differential metabolites were VIP > 1, p < 0.05 and FC ≥ 2 or FC ≤ 0.5. Next, the functions of different metabolites and metabolic pathway analysis by the KEGG database (https://www.kegg.jp/kegg/pathway accessed on May 2022).

2.7. Statistical Analysis

Statistical analysis was performed using SPSS 27.0 software (IBM, Chicago, IL, USA), and the significant difference analysis among the three groups were performed by one-way ANOVA. It was statistically significant when p < 0.05.
To determine the relationship between rumen microorganisms and serum metabolites, spearman correlation analysis was used to analyze the correlation between 30 microbial species with the highest abundance and serum metabolites by R software(v4.1.3, and R package “pheatmap” package in R was used to draw a clustering heat map.

3. Results

3.1. Growth Traits Characteristics of Calves

The results of growth performance and rumen development indexes among the three patterns of calves were showed in Table 1. The slaughter weight of calves in the GF and GFF groups was significantly greater than that in the TMR group (p < 0.05), and the weight gain of calves in the GF and GFF groups also increased but not significantly. The rumen weight and volume of calves in the GF group were significantly higher than those in the other two groups (p < 0.05), while the rumen pH of calves in the GFF group decreased significantly compared to GF and TMR groups (p < 0.05). It was worth noting that t the ratio of rumen to live weight and the ratio of rumen to total stomach weight of calves in the GFF group were higher than those in the other two groups, indicated that the rumen development status of this group of calves was better.
The dominant short-chain fatty acids in rumen fluid were acetic, propionic and butyric acids in the three groups, and the contents of propionic acid and valeric acid in the rumen of calves in the GF group were significantly higher than those in GFF and TMR groups (p < 0.05, Figure 1A). Therefore, the three different feeding patterns achieved different growth promotion effects, but the ‘GF’ feeding pattern was more helpful to the rumen development of calves. Calves in the TMR group were accompanied by slight diarrhea in the late test experiment (Figure 1B), and ‘hairy bulbs’ were found in the rumen of the GFF group (Figure 1C).

3.2. Analysis of the Rumen Metagenome

A total of 247,127,428.4 Mbp of raw data and the clean data of 110,839.92 Mbp were obtained. The GF, GFF and TMR groups had 12,215.79 Mbp, 12,345.66 Mbp and 12,385.19 Mbp average of clean data, respectively. Q20 and Q30 of each sample were above 97.50% and 93.50%. A total of 1,536,918 contiguous clonal lines (N50 length of 30,779 ± 139 bp) were generated, yielding 170,769 ± 18,589 per sample (mean ± standard error) after Metagename assembly (Table S2). The rumen microbes consisted of 84.374% bacteria (320,295,748 sequences), 0.898% archaea (3,408,383 sequences), 0.014% eukaryotes (51,112 sequences), and 0.146% viruses (552,601 sequences) (Table 2).
The results of the domains of calf rumen microbiota under three different feeding patterns are shown in Figure S1. There was a significant difference in Archaea between GF and TMR (p < 0.05). Bacteria and Eukaryota were significantly different between GFF and TMR groups, while Eukaryota was not different (p > 0.05) between GF and TMR groups. Furthermore, there was no significant difference in microbial composition between GF and GFF at the domain level.

3.3. Compositional Characteristics and Taxonomic Differences of the Microbiome

The study found that the dominant bacterial phylum was independent of the three feeding patterns. Bacteroidetes (40.46 ± 3.02%) was dominant, followed by Firmicutes (32.99 ± 3.23%), Proteobacteria (2.20 ± 0.47%) and Actinobacteria (1.91 ± 0.53%) (Table 3). However, its relative abundance was different. Compared with the TMR group, the relative abundance of Firmicutes of GF calves was significantly increased, but Bacteroidetes was significantly decreased (p < 0.05). Similar changes were also observed in the GFF group, but not significantly (Figure 2A). The dominant genera were Prevotella (27.76 ± 2.04%), Clostridium (7.31 ± 2.94%), Olsenella (1.20 ± 0.36%) and Butyrivibrio (1.09 ± 0.17%) (Table S3). At the species level, the analysis showed the relative abundances of Prevotella multisaccharivorax and Prevotella sp. AGR2160 were significantly increased in the GF and GFF groups compared with the TMR group (p < 0.05). The relative abundance of Ruminococcaceae flavefaciens, Prevotella sp. Tc 2-28, Prevotella brevis and Alistipes sp. CAG: 435 was significantly decreased in the GF group (p < 0.05) but not significantly decreased in GFF, and Prevotella ruminicola was significantly decreased in the GF and GFF groups (p < 0.05, Figure 2B).
For the Archaea, compared with the TMR group, the abundance of Euryarchaeota (99.61 ± 0.11%) was extremely significantly higher in both the GF and GFF groups (p < 0.01, Figure S2A). The abundance of Candidatus Woesearchaeota and Candidatus Bathyarchaeota decreased significantly in GF and GFF (p < 0.05). At the genus level, except for Caldivirga, the relative abundance of other genera decreased in GF and GFF groups (p < 0.05, Figure S2B). At the species level, the relative abundances of Methanobrevibacter wolinii, Methanobrevibacter boviskoreani and Methanobrevibacter sp. AbM4 were significantly increased in GF calves, and the relative abundances of Methanobrevibacter ruminantium, Methanobrevibacter olleyae and Methanobrevibacter sp. YE315 were significantly decreased compared with the TMR group (p < 0.05). The calves in the GFF group also had similar but not significant changes (Figure S2C).

3.4. Functional Characteristics and Functional Differences of Microbiome

The results of the KEGG function of the rumen microbial community showed a total of 266 pathways at the third level in GF, GFF and TMR groups (Table S4). They were respectively for ‘Metabolism’ (52.75 ± 1.00%), ‘Genetic Information Processing’ (22.06 ± 1.12%), ‘Environmental Information Processing’ (9.14 ± 0.22%), ‘Cellular Processes’ (7.54 ± 0.13%), ‘Human Diseases’ (5.77 ± 0.07%) and ‘Organismal Systems’ (2.75 ± 0.04%). At the second level, a total of 42 categories were observed, including ‘Carbohydrate metabolism’ (3.31 ± 0.06%), ‘Energy metabolism’ (1.95 ± 0.04%), ‘neurodegenerative diseases’ (0.04 ± 0.01%), ‘Amino acid metabolism’ (2.8 ± 0.09%), ‘Translation’ (2.31 ± 0.15%) and ‘Membrane transport’ (1.31 ± 0.05%) which were the most abundant. At the third level, compared with the TMR group, there were 18 extremely significant enriched third-level pathways (including 5 “metabolic” pathways, 2 “environmental information processing” pathways, 4 “cellular process” pathways, 3 “human disease” pathways, 1 “genetic information processing” pathway and 3 “Organismal Systems” pathways) in the GF group, and only five third-level pathways were extremely significantly enriched in the GFF group (p < 0.01), including “Type II diabetes mellitus”, “Human papillomavirus infection”, “Mineral absorption”, “Nitrogen metabolism” and “Viral carcinogenesis”. Only one Caprolactam degradation was significantly enriched in the GF group compared with the GFF group (p < 0.01, Figure 3).
For the differential KEGG modules involved, 16 GF-enriched and 25 GFF-enriched modules were identified compared to the TMR group (p < 0.01, Figure 4A,B). Nine extremely significantly enriched modules were found between the GF and GFF groups, including 3 GF-enriched modules and 6 GFF-enriched modules (p < 0.01, Figure 4C). Regarding energy metabolism, only and M00172 (C4-dicarboxylic acid cycle) was enriched in the GF calves compared with GFF calves. However, about carbohydrate metabolism, only M00010 (oxaloacetate ≥ 2-oxoglutarate) among the above extremely significantly enriched in the GFF group, six modules enriched in the TMR group (LDA > 3 and p < 0.05, Figure 4D).
For the CAZymes map, a total of 278 genes encoding CAZymes were identified (Table S5), including 124 glycoside hydrolases (GHs), 56 glycosyltransferases (GTs), 17 polysaccharide lyases (PLs), five auxiliary oxidoreductases (AAs) and 58 carbohydrate-binding modules (CBMs). Among these, GH2 (8.14 ± 0.41%) was the dominant gene, followed by GH3 (4.26 ± 0.09%), GH2 (4.10 ± 0.17%), GH43 (3.84 ± 0.32%), GH13 (3.71 ± 0.17%) and GT4 (3.02 ± 0.13%) genes. And in the GF vs. TMR comparison, 26 were enriched in the rumen of GF calves (9GH, 10GT, 2CE and 5CBM) and 26 were also enriched in the rumen of TMR group (16GH, 6CE, 3CBM and 1PL) (Figure 5A); in the GFF vs. TMR comparison, 19 were enriched in the GFF group (9GT, 8GH and 2CE), while 13 were enriched in the TMR group (6GH, 5CE, 1PL and 1CBM) (Figure 5B). Rumen microbial degraded amylum, mainly including EC 3.2.1.1 (α-amylase), EC 2.4.1.1 (starch phosphorylase), EC 3.2.1.3 (glucoamylase), EC 3.2.1.68 (isoamylase) and EC 3.2.1.41 (pullulanase) key enzymes, and the relative abundance of the EC 3.2.1.4 (cellulase), EC 3.2.1.91 (cellobiohydrolase) and EC 3.2.1.21 (β-glucosidase) for cellulose degradation decreased significantly in GF and GFF calves compared to TMR group (p < 0.05, Figure 5C).

3.5. The Relationship between the Types of Rumen Microorganisms and Their Functions

As the fermentation and digestion ability of feed nutrients is an important criterion for rumen development, the study further focused on the function of lipid metabolism and carbohydrate metabolism of rumen microorganisms. It was found that pyruvate metabolism was significantly enriched in GF (p < 0.05), C5 branched dibasic acid metabolism and propanoate metabolism were significantly enriched in GFF (p < 0.05), while pentose and gluconate interconversion and amino sugar and nucleoside sugar metabolism were significantly enriched in the rumen of TMR group (p < 0.05, Figure S3A–C). In the secondary function of Lipid metabolism, Glycolipid metabolism, Glycophospholipid metabolism, Primary bile acid biosynthesis and Secondary bile acid biosynthesis were significantly enriched in GF and GFF (p < 0.05), and Arabic acid metabolism was significantly enriched in TMR (p < 0.05). Interestingly, the two pathways involved in fatty acid metabolism, which were “Fatty acid biosynthesis” and “Fatty acid degradation”, were enriched in the rumen of the TMR group, while Ether lipid metabolism, Synthesis and degradation of ketone bodies were enriched in the rumen of GF and GFF groups. Furthermore, compared with the GF group, Fatty acid biosynthesis was significantly enriched in the rumen of calves in the GFF group (p < 0.05, Figure S3D–F). Notably, in the pyruvate metabolic pathway, the abundance of genes encoding enzymes also changed significantly. The 2.7.1.40 (phosphoenolpyruvate carboxylase) was significantly up-regulated in GF and GFF calves compared with the TMR group, while 2.7.9.1 (pyruvate-phosphate dikinase) was significantly down-regulated (p < 0.05, Figure 6).
By constructing the correlation analysis between bacterial species with lipids and carbohydrate metabolism. Forty-one (41) species were significantly correlated with amino sugar and nucleoside sugar metabolism and pentose and gluconate interconversion (R > 0.5 and p < 0.05, Figure S4). In the secondary functional pathway of lipid metabolism, 59 bacteria were significantly correlated with bile acid biosynthesis (R > 0.5 and p < 0.05), of which 38 were negatively correlated with bile acid biosynthesis. Moreover, a total of 43 species were significantly associated with two fatty acid pathways (R > 0.5 and p < 0.05), of which 22 bacteria were positively correlated with fatty acid biosynthesis (Figure S5). Among the 22 bacterial species positively associated with fatty acid biosynthesis, the five bacteria with the strongest correlation were Prevotellaceae bacterium HUN156, Ruminococcus sp. NK3A76, Treponema berlinense, Ruminococcus flavefaciens and Succinimonas amylolytica.

3.6. Analysis of Serum Metabolomics

The PLS-DA pattern recognition method was used to identify the overall metabolic differences between each pair of groups. The results showed significant differences in metabolites between each pair of groups. The model’s quality parameters were R2 > Q2 and Q2 < 0 (Figure 7A–C), indicating that the PLS-DA model was stable and reliable without overfitting. In addition, PCA analysis showed that the serum metabolites of the three groups of GF, GFF and TMR were well separated (Figure 7D), and the QC samples were tightly bound, indicating that the reproducibility of this experiment was good.
A total of 121 significantly differential metabolites were identified in the GF and TMR groups with VIP score > 1, Fold Change > 1.5 (p < 0.05), and 22 metabolites were classified as lipids and lipid-like molecules (Table S6). Compared with the TMR group, the level of 9 lipids and lipid-like molecules was significantly decreased in the GFF group (Table S7). Furthermore, regarding the key significant metabolites in the GF group, the level of 7 lipids and lipid-like molecules (including Androsterone, Palmitoylcarnitine and 4-Methylvaleric Acid, etc.) were significantly decreased in the GF group compared to those in the GFF group (Table S8). Most of the differential metabolites in the GF and TMR groups were enriched in the vitamin digestion and absorption and bile secretion pathways by KEGG enrichment analysis (Figure 8A,B). Notably, As shown in Figure 8C, significant changes in Steroid hormone biosynthesis occurred in the GF group compared to the GFF calves (p < 0.05).

3.7. Correlation Analysis of the Metabolites and Microorganisms

Pearson correlation analysis was performed on the top 30 microbial taxa in relative abundance and 24 lipids and lipid-like modules metabolites (Figure 9). The results showed a significant correlation between 212 pairs (29.44%) out of 720 pairs (p < 0.05). Among these pairs, 95 pairs (44.81%) were negatively correlated, and 117 (55.19%) were positively correlated. Most of the microbes that were significantly associated with lipids and lipid-like molecule metabolites belong to the genus Prevotella and Ruminococcaceae.

4. Discussion

Through comprehensive analysis of phenotypic indexes, rumen metagenome and serum metabolism of Chinese Holstein bull calves, the study preliminarily explored the microbial and host metabolism molecular mechanism of interaction possibly affecting the animal growth and rumen development of bull calves. It evaluated the structural composition of rumen microorganisms under different feeding patterns and the effects of host metabolism on bull calves’ individual growth and rumen development, which selected the most suitable feeding pattern for calves from a molecular perspective.
The stability of rumen pH is the premise of rumen microbial function and the guarantee for the normal function of the rumen [38]. The rumen pH range of 6.0~7.0 was suitable for rumen nutrient digestion and absorption [39], when pH < 6.5, is not conducive to the digestion of cellulose, and the buffering capacity of the rumen is significantly reduced. When pH < 6.0, the activity of cellulose decomposition bacteria is inhibited or even killed, and the buffering capacity of the rumen is significantly reduced [40]. The study’s results on rumen pH value indicated that the digestibility of GF group calves was the strongest, followed by TMR group calves, and GFF group calves were the weakest. Although the rumen fluid pH has no direct impact on the development of the rumen, Krehbiel et al. [41] found that it could change the proportion of butyric acid in the rumen fluid to cause changes in the proportion of SCFA, which was the best stimulus for the development of the rumen. The absorption of butyric acid also increased with the decrease in pH [42]. Simultaneously, VFA in the rumen was an essential component in maintaining the growth and reproduction of ruminants, which could meet about 70% of the energy demand [43,44,45], so the concentration of VFA was an important reference index for rumen fermentation. The results of the study showed that the concentrations of short-chain fatty acids such as acetic acid, propionic acid, and butyric acid in the rumen of three groups of calves were significantly different, indicating that the abundance of VFA-related microorganisms in the rumen of calves was changed under different feeding patterns, further affecting the content of VFA in the rumen. Particularly, the A/P (acetic acid to propionic acid ratios) in the rumen of the GF group calves was significantly reduced compared with the TMR group, indicating that the calves with GF feeding pattern were more conducive to rumen fermentation, increasing the yield of propionic acid and reducing the ethylene propylene ratio, which was beneficial to improving feed utilization, similar to the research results of Kang et al. [46].
Bacteria are the most abundant microbial community in the rumen of dairy cows, which is consistent with the results found in previous studies [47]. In the study, the difference in rumen microorganisms among GF, GFF and TMR calves was mainly shown in bacteria. The dominant phylum in the rumen of calves were Bacteroidetes, Firmicutes and Proteobacteria, and the dominant genus was Prevotella, which is consistent with the results of Jami et al. [48] analysis of dominant microorganisms in the rumen of Holstein cows. It was found that bacteria were the main participants in the degradation and fermentation of most macromolecular substances in feed, indicating that bacteria played an important role in promoting the growth and development of the host rumen than archaea, viruses, and eukaryotes. The results showed that the relative abundance of Prevotella multisaccharivorax and Prevotella sp. AGR 2160 were significantly increased in the rumen of GF and GFF calves compared to the TMR group. Prevotella might improve the glucose metabolism of the host by promoting increased glycogen storage [49], Prevotella multisaccharivorax fermented sugars (including glucose, glycerol, lactose, maltose, etc.). It also variably fermented arabinose and salicin, and the products of fermentation are acids; two of the significant products formed are acetic acid and succinic acid [50], which indicated that the rumen fermentation ability of GF and GFF groups of calves for carbohydrates in feed was stronger. Notably, at the species level, most of the species that showed higher abundance in the rumens of the three groups of calves belonged to the Prevotella genus. This genus could produce succinic acid and acetate from starch and proteins in feed and was one of the core bacteria in the rumen of dairy cows [51].
Regarding archaea, Metanobrevibacter, Metanosphaera and Metanolacia had high relative abundance in GF and GFF groups, which could use hydrogen, carbon dioxide and formate as the substrate for methane production [52], and use hydrogen in the rumen anaerobic environment as the final way of carbohydrate fermentation, so that microorganisms involved in fermentation could play the best function and support the complete oxidation of the substrate. If hydrogen for carbohydrate fermentation is not removed, it might inhibit the metabolism of rumen microorganisms [53]. This meant that rumen fermentation of calves with GF and GFF feeding patterns produced more hydrogen, which increased the consumption of carbohydrates and provided more energy [54].
Previous research has shown that GTs were responsible for catalyzing the glycosidic bond cleavage and were the second abundant CAZy family (accounting for 19–24% of the total CAZymes) [55]. CBMs didn’t have enzymatic activity, but they could help Cazymes (such as GHs, CEs, and AAs) combine with polysaccharides to enhance their activity [56]. In the study, GF and GFF calves had higher GT abundance than TMR calves, indicating that the rumen of calves with these two feeding patterns was more conducive to the phosphorylation of carbohydrate nutrients, which was helpful to the digestion of carbohydrates. In combination with the abundance changes of key enzymes that degrade starch and cellulose, the study speculated that the GF group increased the amount of the binding activity reduces the content of starch and cellulose in the rumen, which reduced the gene abundance of glucoamylase (EC 3.2.1.3), cellulase (EC 3.2.1.4), cellobiohydrolase (EC 3.2.1.91) and β-The gene abundance of glucosidase (EC 3.2.1.21) in the rumen.
An important function identified in the study was pyruvate metabolism. For ruminants, pyruvate was an intermediate product of carbohydrate fermentation by rumen microorganisms [57], and the ATP released by the conversion of pyruvate into volatile fatty acids (VFAs) in the rumen could promote the growth of rumen microorganisms. The results showed that pyruvate metabolism and glycolysis/gluconeogenesis pathway were enriched in the rumen of GF calves. Phosphoenolpyruvate (PEP) and pyruvate could be transformed mutually, and this step was not only an important part of the pyruvate metabolic pathway but also the last step of the glycolysis pathway [58]. In previous studies, dietary supplementation of calcium pyruvate increased rumen VFAs content, reduced dietary protein degradation, and reduced rumen NH3-H, suggesting that calcium pyruvate was beneficial to rumen metabolism [59,60]. In pyruvate metabolism, the down-regulation of pyruvate-phosphate dikinase (EC: 2.7.9.1) and the up-regulation of phosphoenol-pyruvate carboxylase (EC: 4.1.1.31) promoted the last step of the glycolytic pathway (conversion of PEP to pyruvate), and a large amount of ATP was produced, which could promote rumen microorganisms Valine, leucine, and isoleucine biosynthesis. The results indicated that calves could promote rumen microbial protein synthesis under GF feeding mode [61], improve the utilization rate of rumen nitrogen and promote the rumen fermentation efficiency of calves. Additionally, the results showed that the species with a significant positive correlation with the pyruvate metabolic pathway belong to the Clostridiales and the Methanobrevibacter genus, which indicated that Clostridiales and Methanobrevibacter play an important role in pyruvate metabolism.
Besides carbohydrate metabolism, the functions of rumen microorganisms related to lipid metabolism also have important physiological roles in rumen development [62,63,64,65]. Therefore, the study further focused on the lipid metabolism function of rumen microorganisms. It was found that Arachidonic acid metabolism was significantly enriched in the rumen microbiota of TMR group calves. Arachidonic acid metabolism was an important part of the inflammatory cascade reaction, which could lead to the production of prostaglandins and leukotrienes, most of which were strong proinflammatory factors [66,67,68]. And it might cause abnormal secretion of gastric acid and gastrointestinal hormones, leading to insufficient gastrointestinal motility, and extremely prone to feeding intolerance [69], which indicates that feeding the TMR diet to calves might cause and aggravate rumen inflammation reaction, thereby inhibiting the development of rumen, reducing the fermentation capacity of rumen, and leading to significant reduction in weight gain of calves in TMR group. Bile acid was a major component of bile, which could be secreted into the gastrointestinal tract after eating to promote the digestion and absorption of triacylglycerol, cholesterol and fat-soluble vitamins [70], and it was essential to maintain a healthy gastrointestinal microbiota, balance lipid and carbohydrate metabolism, insulin sensitivity and innate immunity [71]. BAs had an emulsifying effect on dietary fat, increasing the surface area of fat to make it easier to be digested by trypsin [72]. Some essential vitamins (A, D, E and K) were non-polar lipids and could be absorbed by the body only with the participation of BAs [73,74].
Additionally, the bile acid cycle plays a key role in physiological cholesterol homeostasis [75]. In the study, the primary and secondary bill acid biosynthesis were significantly enriched in GF and GFF groups, indicating that GF and GFF calves might have a stronger emulsifying effect on the fat in the feed, which was more conducive to the digestion and absorption of lipids, and provide nutrients needed for the host organism metabolism. Interestingly, the two pathways related to fatty acid metabolism were enriched in the rumen of TMR calves. Previous studies have found that fatty acids were precursors of the synthesis of inflammatory mediators, which played an important role in regulating innate immunity and might contribute to the development of some autoimmune diseases [76,77,78]. Indicated that feeding TMR to calves might cause rumen inflammation, which was also consistent with the previous analysis results.
Untargeted metabolomics analysis of serum showed significant differences in metabolites among GF, GFF and TMR groups. Compared with the TMR group, Hexanoylcarnitine was significantly decreased in GF and GFF calves, and Propionylcarnitine was significantly increased in GF calves. Carnitine is a type B fatty acid-produced process of fatty acid esterification, which plays an important role in regulating muscle oxidative metabolism [79]. Caprylcarnitine and propionyl carnitine are the long-chain acylcarnitines (LCACs) fatty acid derivatives of carnitine [80]. The accumulation of most LCACs during fatty acid oxidation might lead to insulin resistance [81]. In the study, the level of Hexanoylcarnitine in TMR calves was significantly higher than that in GF and GFF calves, which might lead to the disturbance of fatty acids and glucose circulation in the body, and further inhibit the growth rate of calves. Concomitantly, the level of Propionylcarnitine in the organism might be a sign of cobalamin deficiency and had beneficial effects on heart metabolism and blood circulation [82]. Propionylcarnitine had a higher level in GF calves, indicating that the metabolism of GF calves was more vigorous. KEGG enrichment of differential metabolites showed that most of the differential metabolites of GF and GFF were enriched in Bile secret pathways, and the related metabolite Ursodeoxic acid was significantly up-regulated in the GF calf group. Ursodeoxycholic acid was formed by the transformation of deoxycholic acid, which helped regulate lipid metabolism and the integrity of the intestinal barrier [83]. And bile acids could be used as a signal molecule to participate in energy metabolism and other regulatory processes by activating bile acid receptors [84]. Thus, the growthiness of calves might be related to the activation of bile acid receptors in the organism under the GF feeding pattern.
Importantly, it was found that there was a very close relationship between rumen microorganisms and lipid metabolites. Most of the microorganisms that have a significant correlation with most lipid metabolites belong to the Prevotella and Ruminococcaceae genus. Prevotella was beneficial to microbial metabolism to produce short-chain fatty acids [85,86], which could protect the gastrointestinal mucosal barrier, regulate immune function [87,88], and alleviate the effects of rumen structure and metabolic disorders and immune disorders caused by the feeding process. Ruminococcaceae was the main cellulolytic bacteria in the rumen, which was related to the release of inflammatory and cytotoxic factors and could participate in immune regulation and maintain a healthy homeostasis environment [89]. It was a potentially beneficial bacterium. Some studies showed that microorganisms could act as regulators of lipid digestion, absorption, storage and secretion in the gastrointestinal tract [90,91,92], which indicated that their synergistic action with the host could hydrolyze lipid substances more efficiently, which was conducive to the absorption of fat by the organism, and was of great significance to the growth and development of calves.
The study further demonstrated that feeding only concentrated feed to calves could cause a series of inflammatory reactions, lead to immune system disorders, and affect the normal growth and development of the rumen, but there were certain differences from previous research results at the genetic level [93]. Although the study found reliable results in rumen microorganisms body metabolism, due to the small number of test samples and the lack of relevant data on rumen microbial metabolism, further research will be needed in the individual growth, rumen development and disease of calves.

5. Conclusions

Calves achieved similar feeding results with “concentrate plus hay”, “concentrate” and “TMR” feeding patterns, but the rumen pH of calves was more suitable for the function of microorganisms with “concentrate plus hay” feeding patterns. The composition of microorganisms in “concentrate plus hay” feeding pattern calves was more helpful to the digestion of carbohydrates and alleviated inflammation. Simultaneously, the synergistic effect with the host could hydrolyze lipid substances more efficiently and promote the absorption of fat by the organism, which was of great significance to the growth of calves. Due to the small sample size of this study, further research should increase the sample size of each group and focus on the interaction between rumen microorganisms, microbial metabolism, and host metabolism.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/microorganisms11102423/s1, Figure S1: Comparison of microbial domains of rumen in three groups; Figure S2: Comparison of microbial in rumen archaea among the three groups. (A) Comparison of archaea at phylum level. (B) Comparison of archaea at the genus level. (C) Comparison of archaea at the species level; Figure S3: Differential KEGG functions among the three groups. (A, B and C) Comparison of Carbohydrate metabolism of rumen microbial among the three groups. (D, E and F) Comparison of lipid metabolism of rumen microbial among the three groups. * p < 0.05, significant difference, ** p < 0.01, extremely significant difference; Figure S4: Correlation analysis between top 100 microorganisms and carbohydrate metabolism pathways at species level; Figure S5: Correlation analysis between top 100 microorganisms and lipid metabolism pathways at species level; Table S1: The ingredients in the diets and the nutritional composition; Table S2: Sequence data summary of rumen contents samples generated; Table S3: Relative abundances of the top 30 components of the microbiota at the genus level; Table S4: Composition of metabolic pathways based on the first-level and second-level in the KEGG; Table S5: Composition of CAZymes based on the class-level and family-level enzymes; Table S6: Significant metabolites between the GF and TMR groups; Table S7: Significant metabolites between the GFF and TMR groups; Table S8: Significant metabolites between the GF and GFF groups.

Author Contributions

S.L., J.W., W.S. and X.J. contributed to the conception and design of the study. M.L., H.F. and X.B. collected samples. K.Z., S.X. and Y.C. conducted the research and performed the statistical analysis. K.Z. wrote the first draft of the manuscript. J.W., Z.L., M.W. and J.N. wrote sections and modified the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by high-quality cattle breeding materials and methods innovation and new variety breeding, a key R & D project of Sichuan Province (2021YFYZ0001); National modern agricultural industry technology system Sichuan innovation team technical R&D and integrated application post of beef cattle industry in agricultural area (SCCXTD-2022-13).

Institutional Review Board Statement

The authors confirm that this study was performed in accordance with the Guidelines of Good Experimental Practices adopted by the Institute of Animal Science of the Sichuan Agricultural University, Chengdu, China. All experimental protocols involving animals were approved by the Animal Care and Use Committee for Biological Studies, Sichuan Province, China (DKY-B2019302083).

Informed Consent Statement

Not applicable.

Data Availability Statement

All the figures and tables used to support the results of this study are included.

Acknowledgments

We thank other authors and employees of Hongfeng Dairy Farm of Sichuan Xuebao Dairy Group for help with the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Monteverde, V.; Congiu, F.; Vazzana, I.; Dara, S.; Di Pietro, S.; Piccione, G. Serum lipid profile modification related to polyunsaturated fatty acid supplementation in thoroughbred horses. J. Appl. Anim. Res. 2017, 45, 615–618. [Google Scholar] [CrossRef]
  2. Liu, F.; Hogg, J.; Kracht, S.; Brewster, C.J.; Henman, D.J.; Athorn, R.Z.; Morrison, R.S.; Smits, R.J.; Campbell, R.G. Supplementing 2 g per day bovine lactoferrin from late gestation until weaning did not improve lactation performance of mixed parity sows. Anim. Prod. Sci. 2019, 59, 615–618. [Google Scholar] [CrossRef]
  3. Lin, L.; Xie, F.; Sun, D.; Liu, J.; Zhu, W.; Mao, S. Ruminal microbiome-host crosstalk stimulates the development of the ruminal epithelium in a lamb model. Microbiome 2019, 7, 83. [Google Scholar] [CrossRef] [PubMed]
  4. Soberon, F.; Raffrenato, E.; Everett, R.W.; Van Amburgh, M.E. Preweaning milk replacer intake and effects on long-term productivity of dairy calves. J. Dairy Sci. 2012, 95, 783–793. [Google Scholar] [CrossRef]
  5. Zitnan, R.; Voigt, J.; Schönhusen, U.; Wegner, J.; Kokardová, M.; Hagemeister, H.; Levkut, M.; Kuhla, S.; Sommer, A. Influence of dietary concentrate to forage ratio on the development of rumen mucosa in calves. Arch. Tierernahr. 1998, 51, 279–291. [Google Scholar] [CrossRef]
  6. Overvest, M.A.; Bergeron, R.; Haley, D.B.; DeVries, T.J. Effect of feed type and method of presentation on feeding behavior, intake, and growth of dairy calves fed a high level of milk. J. Dairy Sci. 2016, 99, 317–327. [Google Scholar] [CrossRef]
  7. Baldwin, R.L.; McLeod, K.R.; Klotz, J.L.; Heitmann, R.N. Rumen Development, Intestinal Growth and Hepatic Metabolism in the Pre- and Postweaning Ruminant. J. Dairy Sci. 2004, 87, E55–E65. [Google Scholar] [CrossRef]
  8. Mojahedi, S.; Khorvash, M.; Ghorbani, G.R.; Ghasemi, E.; Mirzaei, M.; Hashemzadeh-Cigari, F. Performance, nutritional behavior, and metabolic responses of calves supplemented with forage depend on starch fermentability. J. Dairy Sci. 2018, 101, 7061–7072. [Google Scholar] [CrossRef]
  9. Hosseini, S.M.; Ghorbani, G.R.; Rezamand, P.; Khorvash, M. Determining optimum age of Holstein dairy calves when adding chopped alfalfa hay to meal starter diets based on measures of growth and performance. Animal 2016, 10, 607–615. [Google Scholar] [CrossRef]
  10. Jami, E.; White, B.A.; Mizrahi, I. Potential role of the bovine rumen microbiome in modulating milk composition and feed efficiency. PLoS ONE 2014, 9, e85423. [Google Scholar] [CrossRef]
  11. Wallace, R.J.; Sasson, G.; Garnsworthy, P.C.; Tapio, I.; Gregson, E.; Bani, P.; Huhtanen, P.; Bayat, A.R.; Strozzi, F.; Biscarini, F.; et al. A heritable subset of the core rumen microbiome dictates dairy cow productivity and emissions. Sci. Adv. 2019, 5, eaav8391. [Google Scholar] [CrossRef] [PubMed]
  12. Ye, Z.; Xu, Y.-J.; Liu, Y. Influences of dietary oils and fats, and the accompanied minor content of components on the gut microbiota and gut inflammation: A review. Trends Food Sci. Technol. 2021, 113, 255–276. [Google Scholar] [CrossRef]
  13. Lazuka, A.; Auer, L.; O’Donohue, M.; Hernandez-Raquet, G. Anaerobic lignocellulolytic microbial consortium derived from termite gut: Enrichment, lignocellulose degradation and community dynamics. Biotechnol. Biofuels 2018, 11, 284. [Google Scholar] [CrossRef]
  14. Fernando, S.C.; Purvis, H.T., 2nd; Najar, F.Z.; Sukharnikov, L.O.; Krehbiel, C.R.; Nagaraja, T.G.; Roe, B.A.; Desilva, U. Rumen microbial population dynamics during adaptation to a high-grain diet. Appl. Environ. Microbiol. 2010, 76, 7482–7490. [Google Scholar] [CrossRef]
  15. Xie, F.; Anderson, C.L.; Timme, K.R.; Kurz, S.G.; Fernando, S.C.; Wood, J.R. Obesity-Dependent Increases in Oocyte mRNAs Are Associated With Increases in Proinflammatory Signaling and Gut Microbial Abundance of Lachnospiraceae in Female Mice. Endocrinology 2016, 157, 1630–1643. [Google Scholar] [CrossRef] [PubMed]
  16. Poudel, P.; Froehlich, K.; Casper, D.P.; St-Pierre, B. Feeding Essential Oils to Neonatal Holstein Dairy Calves Results in Increased Ruminal Prevotellaceae Abundance and Propionate Concentrations. Microorganisms 2019, 7, 120. [Google Scholar] [CrossRef] [PubMed]
  17. Leal, L.N.; Doelman, J.; Keppler, B.R.; Steele, M.A.; Martín-Tereso, J. Preweaning nutrient supply alters serum metabolomics profiles related to protein and energy metabolism and hepatic function in Holstein heifer calves. J. Dairy Sci. 2021, 104, 7711–7724. [Google Scholar] [CrossRef]
  18. Qi, Y.; Zhao, X.; Huang, D.; Pan, X.; Yang, Y.; Zhao, H.; Hu, H.; Cheng, G. Exploration of the Relationship between Intestinal Colostrum or Milk, and Serum Metabolites in Neonatal Calves by Metabolomics Analysis. J. Agric. Food Chem. 2018, 66, 7200–7208. [Google Scholar] [CrossRef]
  19. Terler, G.; Poier, G.; Klevenhusen, F.; Zebeli, Q. Replacing concentrates with a high-quality hay in the starter feed in dairy calves: I. Effects on nutrient intake, growth performance, and blood metabolic profile. J. Dairy Sci. 2022, 105, 2326–2342. [Google Scholar] [CrossRef]
  20. Taylor, E.N.; Han, J.; Fan, C.; Beckmann, M.; Hewinson, G.; Rooke, D.; Koets, A.P.; Mur, L.A.J. Defining Fatty Acid Changes Linked to Rumen Development, Weaning and Growth in Holstein-Friesian Heifers. Metabolites 2022, 12, 374. [Google Scholar] [CrossRef]
  21. Karlsson, F.H.; Fåk, F.; Nookaew, I.; Tremaroli, V.; Fagerberg, B.; Petranovic, D.; Bäckhed, F.; Nielsen, J. Symptomatic atherosclerosis is associated with an altered gut metagenome. Nat. Commun. 2012, 3, 1245. [Google Scholar] [CrossRef]
  22. Fosso, B.; Santamaria, M.; D’Antonio, M.; Lovero, D.; Corrado, G.; Vizza, E.; Passaro, N.; Garbuglia, A.R.; Capobianchi, M.R.; Crescenzi, M.; et al. MetaShot: An accurate workflow for taxon classification of host-associated microbiome from shotgun metagenomic data. Bioinformatics 2017, 33, 1730–1732. [Google Scholar] [CrossRef] [PubMed]
  23. Nielsen, H.B.; Almeida, M.; Juncker, A.S.; Rasmussen, S.; Li, J.; Sunagawa, S.; Plichta, D.R.; Gautier, L.; Pedersen, A.G.; Le Chatelier, E.; et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat. Biotechnol. 2014, 32, 822–828. [Google Scholar] [CrossRef] [PubMed]
  24. Karlsson, F.H.; Tremaroli, V.; Nookaew, I.; Bergström, G.; Behre, C.J.; Fagerberg, B.; Nielsen, J.; Bäckhed, F. Gut metagenome in European women with normal, impaired and diabetic glucose control. Nature 2013, 498, 99–103. [Google Scholar] [CrossRef]
  25. Zeller, G.; Tap, J.; Voigt, A.Y.; Sunagawa, S.; Kultima, J.R.; Costea, P.I.; Amiot, A.; Böhm, J.; Brunetti, F.; Habermann, N.; et al. Potential of fecal microbiota for early-stage detection of colorectal cancer. Mol. Syst. Biol. 2014, 10, 766. [Google Scholar] [CrossRef] [PubMed]
  26. Sunagawa, S.; Coelho, L.P.; Chaffron, S.; Kultima, J.R.; Labadie, K.; Salazar, G.; Djahanschiri, B.; Zeller, G.; Mende, D.R.; Alberti, A.; et al. Ocean plankton. Structure and function of the global ocean microbiome. Science 2015, 348, 1261359. [Google Scholar] [CrossRef]
  27. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 2012, 28, 3150–3152. [Google Scholar] [CrossRef]
  28. Qin, J.; Li, Y.; Cai, Z.; Li, S.; Zhu, J.; Zhang, F.; Liang, S.; Zhang, W.; Guan, Y.; Shen, D.; et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 2012, 490, 55–60. [Google Scholar] [CrossRef]
  29. Oh, J.; Byrd, A.L.; Deming, C.; Conlan, S.; Kong, H.H.; Segre, J.A. Biogeography and individuality shape function in the human skin metagenome. Nature 2014, 514, 59–64. [Google Scholar] [CrossRef]
  30. Feng, Q.; Liang, S.; Jia, H.; Stadlmayr, A.; Tang, L.; Lan, Z.; Zhang, D.; Xia, H.; Xu, X.; Jie, Z.; et al. Gut microbiome development along the colorectal adenoma-carcinoma sequence. Nat. Commun. 2015, 6, 6528. [Google Scholar] [CrossRef]
  31. Villar, E.; Farrant, G.K.; Follows, M.; Garczarek, L.; Speich, S.; Audic, S.; Bittner, L.; Blanke, B.; Brum, J.R.; Brunet, C.; et al. Ocean plankton. Environmental characteristics of Agulhas rings affect interocean plankton transport. Science 2015, 348, 1261447. [Google Scholar] [CrossRef] [PubMed]
  32. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef] [PubMed]
  33. Avershina, E.; Frisli, T.; Rudi, K. De novo semi-alignment of 16S rRNA gene sequences for deep phylogenetic characterization of next generation sequencing data. Microbes Environ. 2013, 28, 211–216. [Google Scholar] [CrossRef]
  34. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [PubMed]
  35. Kanehisa, M.; Goto, S.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. Data, information, knowledge and principle: Back to metabolism in KEGG. Nucleic Acids Res. 2014, 42, D199–D205. [Google Scholar] [CrossRef]
  36. Scher, J.U.; Sczesnak, A.; Longman, R.S.; Segata, N.; Ubeda, C.; Bielski, C.; Rostron, T.; Cerundolo, V.; Pamer, E.G.; Abramson, S.B.; et al. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis. eLife 2013, 2, e01202. [Google Scholar] [CrossRef]
  37. Wen, B.; Mei, Z.; Zeng, C.; Liu, S. metaX: A flexible and comprehensive software for processing metabolomics data. BMC Bioinform. 2017, 18, 183. [Google Scholar] [CrossRef]
  38. Liu, S.J.; Bu, D.P.; Wang, J.Q.; Liu, L.; Liang, S.; Wei, H.Y.; Zhou, L.Y.; Li, D.; Loor, J.J. Effect of incremental levels of fish oil supplementation on specific bacterial populations in bovine ruminal fluid. J. Anim. Physiol. Anim. Nutr. 2012, 96, 9–16. [Google Scholar] [CrossRef]
  39. Wales, W.J.; Kolver, E.S.; Thorne, P.L.; Egan, A.R. Diurnal variation in ruminal pH on the digestibility of highly digestible perennial ryegrass during continuous culture fermentation. J. Dairy Sci. 2004, 87, 1864–1871. [Google Scholar] [CrossRef]
  40. DePeters, E.J.; Bath, D.L. Canola Meal Versus Cottonseed Meal as the Protein Supplement in Dairy Diets. J. Dairy Sci. 1986, 69, 148–154. [Google Scholar] [CrossRef]
  41. Krehbiel, C.R.; Harmon, D.L.; Schneider, J.E. Effect of increasing ruminal butyrate on portal and hepatic nutrient flux in steers. J. Anim. Sci. 1992, 70, 904–914. [Google Scholar] [CrossRef] [PubMed]
  42. McLeod, K.R.; Baldwin, R.L.T. Effects of diet forage:concentrate ratio and metabolizable energy intake on visceral organ growth and in vitro oxidative capacity of gut tissues in sheep. J. Anim. Sci. 2000, 78, 760–770. [Google Scholar] [CrossRef] [PubMed]
  43. Ma, J.; Wang, C.; Wang, Z.; Cao, G.; Hu, R.; Wang, X.; Zou, H.; Kang, K.; Peng, Q.; Xue, B.; et al. Active dry yeast supplementation improves the growth performance, rumen fermentation, and immune response of weaned beef calves. Anim. Nutr. 2021, 7, 1352–1359. [Google Scholar] [CrossRef] [PubMed]
  44. Spears, J.W.; Schlegel, P.; Seal, M.C.; Lloyd, K.E. Bioavailability of zinc from zinc sulfate and different organic zinc sources and their effects on ruminal volatile fatty acid proportions. Livest. Prod. Sci. 2004, 90, 211–217. [Google Scholar] [CrossRef]
  45. Matthews, C.; Crispie, F.; Lewis, E.; Reid, M.; O’Toole, P.W.; Cotter, P.D. The rumen microbiome: A crucial consideration when optimising milk and meat production and nitrogen utilisation efficiency. Gut Microbes 2019, 10, 115–132. [Google Scholar] [CrossRef]
  46. Kang, J.; Zeng, B.; Tang, S.; Wang, M.; Han, X.; Zhou, C.; Yan, Q.; Liu, J.; Tan, Z. Effects of Momordica charantia polysaccharide on in vitro ruminal fermentation and cellulolytic bacteria. Ital. J. Anim. Sci. 2017, 16, 226–233. [Google Scholar] [CrossRef]
  47. Shabat, S.K.; Sasson, G.; Doron-Faigenboim, A.; Durman, T.; Yaacoby, S.; Berg Miller, M.E.; White, B.A.; Shterzer, N.; Mizrahi, I. Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 2016, 10, 2958–2972. [Google Scholar] [CrossRef]
  48. Jami, E.; Israel, A.; Kotser, A.; Mizrahi, I. Exploring the bovine rumen bacterial community from birth to adulthood. ISME J. 2013, 7, 1069–1079. [Google Scholar] [CrossRef]
  49. Kovatcheva-Datchary, P.; Nilsson, A.; Akrami, R.; Lee, Y.S.; De Vadder, F.; Arora, T.; Hallen, A.; Martens, E.; Björck, I.; Bäckhed, F. Dietary Fiber-Induced Improvement in Glucose Metabolism Is Associated with Increased Abundance of Prevotella. Cell Metab. 2015, 22, 971–982. [Google Scholar] [CrossRef]
  50. Sakamoto, M.; Umeda, M.; Ishikawa, I.; Benno, Y. Prevotella multisaccharivorax sp. nov., isolated from human subgingival plaque. Int. J. Syst. Evol. Microbiol. 2005, 55, 1839–1843. [Google Scholar] [CrossRef]
  51. Stevenson, D.M.; Weimer, P.J. 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. 2007, 75, 165–174. [Google Scholar] [CrossRef] [PubMed]
  52. Balch, W.E.; Fox, G.E.; Magrum, L.J.; Woese, C.R.; Wolfe, R.S. Methanogens: Reevaluation of a unique biological group. Microbiol. Rev. 1979, 43, 260–296. [Google Scholar] [CrossRef] [PubMed]
  53. Sharp, R.; Ziemer, C.J.; Stern, M.D.; Stahl, D.A. Taxon-specific associations between protozoal and methanogen populations in the rumen and a model rumen system. FEMS Microbiol. Ecol. 1998, 26, 71–78. [Google Scholar] [CrossRef]
  54. Hook, S.E.; Wright, A.D.; McBride, B.W. Methanogens: Methane producers of the rumen and mitigation strategies. Archaea 2010, 2010, 945785. [Google Scholar] [CrossRef]
  55. Mota, T.R.; Oliveira, D.M.D.; Marchiosi, R.; Ferrarese-Filho, O.; Santos, W.D.D. Plant cell wall composition and enzymatic deconstruction. AIMS Bioeng. 2018, 5, 63–77. [Google Scholar] [CrossRef]
  56. Bernardes, A.; Pellegrini, V.O.A.; Curtolo, F.; Camilo, C.M.; Mello, B.L.; Johns, M.A.; Scott, J.L.; Guimaraes, F.E.C.; Polikarpov, I. Carbohydrate binding modules enhance cellulose enzymatic hydrolysis by increasing access of cellulases to the substrate. Carbohydr. Polym. 2019, 211, 57–68. [Google Scholar] [CrossRef]
  57. Petzel, E.A.; Titgemeyer, E.C.; Smart, A.J.; Hales, K.E.; Foote, A.P.; Acharya, S.; Bailey, E.A.; Held, J.E.; Brake, D.W. What is the digestibility and caloric value of different botanical parts in corn residue to cattle?1. J. Anim. Sci. 2019, 97, 3056–3070. [Google Scholar] [CrossRef]
  58. Chen, J.; Mitra, R.; Zhang, S.; Zuo, Z.; Lin, L.; Zhao, D.; Xiang, H.; Han, J. Unusual Phosphoenolpyruvate (PEP) Synthetase-Like Protein Crucial to Enhancement of Polyhydroxyalkanoate Accumulation in Haloferax mediterranei Revealed by Dissection of PEP-Pyruvate Interconversion Mechanism. Appl. Environ. Microbiol. 2019, 85, e00984-19. [Google Scholar] [CrossRef]
  59. Chen, X.; Wang, X.; Wang, X.; Yang, J.; Shen, X.; Jie, C. Effect of Calcium Pyruvate on Rumen Metabolism of Goats. Chin. Agric. Sci. Bull. 2006, 12, 17–20. [Google Scholar]
  60. Ran, L.W.; Shen, X.Z.; Jie, T.U.; Lin, P.F.; Chen, J. Study on the Effect of Calcium Pyruvate on the Dynamics of Rumen VFA Concentration in Goats. Acta Ecol. Anim. Domastici 2008, 2, 41–44. [Google Scholar]
  61. Allison, M.J.; Bucklin, J.A.; Robinson, I.M. Importance of the isovalerate carboxylation pathway of leucine biosynthesis in the rumen. Appl. Microbiol. 1966, 14, 807–814. [Google Scholar] [CrossRef] [PubMed]
  62. Mackie, R.I.; White, B.A.; Bryant, M.P. Lipid metabolism in anaerobic ecosystems. Crit. Rev. Microbiol. 1991, 17, 449–479. [Google Scholar] [CrossRef] [PubMed]
  63. Dewanckele, L.; Toral, P.G.; Vlaeminck, B.; Fievez, V. Invited review: Role of rumen biohydrogenation intermediates and rumen microbes in diet-induced milk fat depression: An update. J. Dairy Sci. 2020, 103, 7655–7681. [Google Scholar] [CrossRef] [PubMed]
  64. Maia, M.R.; Chaudhary, L.C.; Bestwick, C.S.; Richardson, A.J.; McKain, N.; Larson, T.R.; Graham, I.A.; Wallace, R.J. Toxicity of unsaturated fatty acids to the biohydrogenating ruminal bacterium, Butyrivibrio fibrisolvens. BMC Microbiol. 2010, 10, 52. [Google Scholar] [CrossRef]
  65. Bauchart, D.; Legay-Carmier, F.; Doreau, M.; Gaillard, B. Lipid metabolism of liquid-associated and solid-adherent bacteria in rumen contents of dairy cows offered lipid-supplemented diets. Br. J. Nutr. 1990, 63, 563–578. [Google Scholar] [CrossRef] [PubMed]
  66. Stanke-Labesque, F.; Pofelski, J.; Moreau-Gaudry, A.; Bessard, G.; Bonaz, B. Urinary leukotriene E4 excretion: A biomarker of inflammatory bowel disease activity. Inflamm. Bowel Dis. 2008, 14, 769–774. [Google Scholar] [CrossRef]
  67. Milton-Thompson, G.J.; Cummings, J.H.; Newman, A.; Billings, J.A.; Misiewicz, J.J. Colonic and small intestinal response to intravenous prostaglandin F2 alpha and E2 in man. Gut 1975, 16, 42–46. [Google Scholar] [CrossRef]
  68. Sartor, R.B. Mechanisms of disease: Pathogenesis of Crohn’s disease and ulcerative colitis. Nat. Clin. Pract. Gastroenterol. Hepatol. 2006, 3, 390–407. [Google Scholar] [CrossRef]
  69. Hamilton, E.; Massey, C.; Ross, J.; Taylor, S. Early enteral feeding in very low birth weight infants. Early Hum. Dev. 2014, 90, 227–230. [Google Scholar] [CrossRef]
  70. Joyce, S.A.; Gahan, C.G. Bile Acid Modifications at the Microbe-Host Interface: Potential for Nutraceutical and Pharmaceutical Interventions in Host Health. Annu. Rev. Food Sci. Technol. 2016, 7, 313–333. [Google Scholar] [CrossRef]
  71. Jia, W.; Xie, G.; Jia, W. Bile acid-microbiota crosstalk in gastrointestinal inflammation and carcinogenesis. Nat. Rev. Gastroenterol. Hepatol. 2018, 15, 111–128. [Google Scholar] [CrossRef] [PubMed]
  72. Lombardo, D. Bile salt-dependent lipase: Its pathophysiological implications. Biochim. Biophys. Acta 2001, 1533, 1–28. [Google Scholar] [CrossRef] [PubMed]
  73. Thompson, G.R. Absorption of fat-soluble vitamins and sterols. J. Clin. Pathol. Suppl. 1971, 5, 85. [Google Scholar] [CrossRef]
  74. Saeed, A.; Hoekstra, M.; Hoeke, M.O.; Heegsma, J.; Faber, K.N. The interrelationship between bile acid and vitamin A homeostasis. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 2017, 1862, 496–512. [Google Scholar] [CrossRef]
  75. Li, T.; Apte, U. Bile Acid Metabolism and Signaling in Cholestasis, Inflammation, and Cancer. Adv. Pharmacol. 2015, 74, 263–302. [Google Scholar]
  76. Yaqoob, P. Fatty acids as gatekeepers of immune cell regulation. Trends Immunol. 2003, 24, 639–645. [Google Scholar] [CrossRef]
  77. Batista-Gonzalez, A.; Vidal, R.; Criollo, A.; Carreño, L.J. New Insights on the Role of Lipid Metabolism in the Metabolic Reprogramming of Macrophages. Front. Immunol. 2019, 10, 2993. [Google Scholar] [CrossRef]
  78. Frommer, K.W.; Schäffler, A.; Rehart, S.; Lehr, A.; Müller-Ladner, U.; Neumann, E. Free fatty acids: Potential proinflammatory mediators in rheumatic diseases. Ann. Rheum. Dis. 2015, 74, 303–310. [Google Scholar] [CrossRef]
  79. Fan, H.; Li, Y.; Wang, J.; Shao, J.; Tang, T.; Elzo, M.A.; Wang, L.; Lai, T.; Ma, Y.; Gan, M.; et al. Untargeted Metabolomic Characteristics of Skeletal Muscle Dysfunction in Rabbits Induced by a High Fat Diet. Animals 2021, 11, 1722. [Google Scholar] [CrossRef]
  80. Stephens, F.B.; Constantin-Teodosiu, D.; Greenhaff, P.L. New insights concerning the role of carnitine in the regulation of fuel metabolism in skeletal muscle. J. Physiol. 2007, 581, 431–444. [Google Scholar] [CrossRef]
  81. Kim, T.; He, L.; Johnson, M.S.; Li, Y.; Zeng, L.; Ding, Y.; Long, Q.; Moore, J.F.; Sharer, J.D.; Nagy, T.R.; et al. Carnitine Palmitoyltransferase 1b Deficiency Protects Mice from Diet-Induced Insulin Resistance. J. Diabetes Metab. 2014, 5, 361. [Google Scholar] [CrossRef] [PubMed]
  82. Campbell, C.D.; Ganesh, J.; Ficicioglu, C. Two newborns with nutritional vitamin B12 deficiency: Challenges in newborn screening for vitamin B12 deficiency. Haematologica 2005, 90, Ecr45. [Google Scholar] [PubMed]
  83. Goossens, J.F.; Bailly, C. Ursodeoxycholic acid and cancer: From chemoprevention to chemotherapy. Pharmacol. Ther. 2019, 203, 107396. [Google Scholar] [CrossRef] [PubMed]
  84. Molinaro, A.; Wahlström, A.; Marschall, H.U. Role of Bile Acids in Metabolic Control. Trends Endocrinol. Metab. 2018, 29, 31–41. [Google Scholar] [CrossRef]
  85. Kong, C.; Gao, R.; Yan, X.; Huang, L.; Qin, H. Probiotics improve gut microbiota dysbiosis in obese mice fed a high-fat or high-sucrose diet. Nutrition 2019, 60, 175–184. [Google Scholar] [CrossRef]
  86. Song, X.; Zhong, L.; Lyu, N.; Liu, F.; Li, B.; Hao, Y.; Xue, Y.; Li, J.; Feng, Y.; Ma, Y.; et al. Inulin Can Alleviate Metabolism Disorders in ob/ob Mice by Partially Restoring Leptin-related Pathways Mediated by Gut Microbiota. Genom. Proteom. Bioinform. 2019, 17, 64–75. [Google Scholar] [CrossRef]
  87. Burger-van Paassen, N.; Vincent, A.; Puiman, P.J.; van der Sluis, M.; Bouma, J.; Boehm, G.; van Goudoever, J.B.; van Seuningen, I.; Renes, I.B. The regulation of intestinal mucin MUC2 expression by short-chain fatty acids: Implications for epithelial protection. Biochem. J. 2009, 420, 211–219. [Google Scholar] [CrossRef]
  88. Inagaki, A.; Sakata, T. Dose-dependent stimulatory and inhibitory effects of luminal and serosal n-butyric acid on epithelial cell proliferation of pig distal colonic mucosa. J. Nutr. Sci. Vitaminol. 2005, 51, 156–160. [Google Scholar] [CrossRef]
  89. Wang, Y.; Nan, X.; Zhao, Y.; Wang, Y.; Jiang, L.; Xiong, B. Ruminal Degradation of Rumen-Protected Glucose Influences the Ruminal Microbiota and Metabolites in Early-Lactation Dairy Cows. Appl. Environ. Microbiol. 2021, 87, e01908-20. [Google Scholar] [CrossRef]
  90. Rabot, S.; Membrez, M.; Bruneau, A.; Gérard, P.; Harach, T.; Moser, M.; Raymond, F.; Mansourian, R.; Chou, C.J. Germ-free C57BL/6J mice are resistant to high-fat-diet-induced insulin resistance and have altered cholesterol metabolism. FASEB J. 2010, 24, 4948–4959. [Google Scholar]
  91. Sommer, F.; Bäckhed, F. Know your neighbor: Microbiota and host epithelial cells interact locally to control intestinal function and physiology. Bioessays 2016, 38, 455–464. [Google Scholar] [CrossRef] [PubMed]
  92. Martinez-Guryn, K.; Hubert, N.; Frazier, K.; Urlass, S.; Musch, M.W.; Ojeda, P.; Pierre, J.F.; Miyoshi, J.; Sontag, T.J.; Cham, C.M.; et al. Small Intestine Microbiota Regulate Host Digestive and Absorptive Adaptive Responses to Dietary Lipids. Cell Host Microbe 2018, 23, 458–469.e455. [Google Scholar] [CrossRef] [PubMed]
  93. Wang, J.; Fan, H.; Li, M.; Zhao, K.; Xia, S.; Chen, Y.; Shao, J.; Tang, T.; Bai, X.; Liu, Z.; et al. Integration of Non-Coding RNA and mRNA Profiles Reveals the Mechanisms of Rumen Development Induced by Different Types of Diet in Calves. Genes 2023, 14, 1093. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 11 02423 g001
Figure 1. Overview of rumen contents in calves. (A) The concentration of sCFA in rumen fluid and different lowercase letters indicated significant differences between groups (p < 0.05). (B) “Hair bulbs” in the rumen of GFF calves. (C) Calves in the TMR group had slight diarrhea.
Figure 1. Overview of rumen contents in calves. (A) The concentration of sCFA in rumen fluid and different lowercase letters indicated significant differences between groups (p < 0.05). (B) “Hair bulbs” in the rumen of GFF calves. (C) Calves in the TMR group had slight diarrhea.
Microorganisms 11 02423 g001
Microorganisms 11 02423 g002
Figure 2. Comparison of microbial in rumen bacteria among the three groups. (A) Comparison of bacteria at the phylum level. (B) Comparison of bacteria at the species level.
Figure 2. Comparison of microbial in rumen bacteria among the three groups. (A) Comparison of bacteria at the phylum level. (B) Comparison of bacteria at the species level.
Microorganisms 11 02423 g002
Microorganisms 11 02423 g003
Figure 3. The log2FC of extremely significantly enriched metabolic pathways among the three groups, p < 0.01.
Figure 3. The log2FC of extremely significantly enriched metabolic pathways among the three groups, p < 0.01.
Microorganisms 11 02423 g003
Microorganisms 11 02423 g004
Figure 4. Comparison of rumen microbial extremely significantly different KEGG modules among the three groups. (A) Significantly different modules between GF and TMR calves. (B) Significantly different modules between GFF and TMR calves. (C) Significantly different modules between GF and GFF calves. (D) Carbohydrate metabolism. Significant different pathways were tested by Linear discriminant analysis effect size (LEfSe) analysis with linear discriminant analysis (LDA) score of >3 and p < 0.05.
Figure 4. Comparison of rumen microbial extremely significantly different KEGG modules among the three groups. (A) Significantly different modules between GF and TMR calves. (B) Significantly different modules between GFF and TMR calves. (C) Significantly different modules between GF and GFF calves. (D) Carbohydrate metabolism. Significant different pathways were tested by Linear discriminant analysis effect size (LEfSe) analysis with linear discriminant analysis (LDA) score of >3 and p < 0.05.
Microorganisms 11 02423 g004
Microorganisms 11 02423 g005
Figure 5. Differential CAZyme functions among three groups. (A,B) Significantly different CAZyme functions in GF and GFF groups compared with the TMR group, significantly different CAZymes were tested by Linear discriminant analysis effect size (LEfSe) analysis with linear discriminant analysis (LDA) score of >3 and p value of <0.05. (C) Comparisons of the gene abundance of KO enzymes for amylase and cellulase, * significantly different among groups, * p < 0.05, ** p < 0.01.
Figure 5. Differential CAZyme functions among three groups. (A,B) Significantly different CAZyme functions in GF and GFF groups compared with the TMR group, significantly different CAZymes were tested by Linear discriminant analysis effect size (LEfSe) analysis with linear discriminant analysis (LDA) score of >3 and p value of <0.05. (C) Comparisons of the gene abundance of KO enzymes for amylase and cellulase, * significantly different among groups, * p < 0.05, ** p < 0.01.
Microorganisms 11 02423 g005
Microorganisms 11 02423 g006
Figure 6. Differential ECs and ECs involved in pyruvate metabolism (EC indicated the number of Enzyme in the pathway). (A) The ECs of calves in the GF and GFF groups were significantly different from the TMR group. (B) Microbial functions and species involved in Pyruvate Metabolism in the rumen of calves.
Figure 6. Differential ECs and ECs involved in pyruvate metabolism (EC indicated the number of Enzyme in the pathway). (A) The ECs of calves in the GF and GFF groups were significantly different from the TMR group. (B) Microbial functions and species involved in Pyruvate Metabolism in the rumen of calves.
Microorganisms 11 02423 g006
Microorganisms 11 02423 g007
Figure 7. Multivariate statistical analysis of metabolomic. (A) PLS−DA score plots for GF and TMR groups. (B) PLS−DA score plots for GFF and TMR groups. (C) PLS−DA score plots for GF and GFF groups. (D) PCA score plots for serum samples in the three groups. Each point represents an individual sample.
Figure 7. Multivariate statistical analysis of metabolomic. (A) PLS−DA score plots for GF and TMR groups. (B) PLS−DA score plots for GFF and TMR groups. (C) PLS−DA score plots for GF and GFF groups. (D) PCA score plots for serum samples in the three groups. Each point represents an individual sample.
Microorganisms 11 02423 g007
Microorganisms 11 02423 g008
Figure 8. KEGG enrichment analysis of differential metabolites between the three groups. (A) KEGG enrichment analysis of differential metabolites between GF and TMR groups. (B) KEGG enrichment analysis of differential metabolites between GFF and TMR groups. (C) KEGG enrichment analysis of differential metabolites between GF and GFF groups.
Figure 8. KEGG enrichment analysis of differential metabolites between the three groups. (A) KEGG enrichment analysis of differential metabolites between GF and TMR groups. (B) KEGG enrichment analysis of differential metabolites between GFF and TMR groups. (C) KEGG enrichment analysis of differential metabolites between GF and GFF groups.
Microorganisms 11 02423 g008
Microorganisms 11 02423 g009
Figure 9. Correlation analysis between the significantly differential metabolites and rumen bacterial taxa identified * p < 0.05, significant correlation, ** p < 0.01, extremely significant correlation.
Figure 9. Correlation analysis between the significantly differential metabolites and rumen bacterial taxa identified * p < 0.05, significant correlation, ** p < 0.01, extremely significant correlation.
Microorganisms 11 02423 g009
Table 1. Growth performance and rumen development parameters of calves in GF, GFF and TMR Groups.
Table 1. Growth performance and rumen development parameters of calves in GF, GFF and TMR Groups.
ParametersGFGFFTMR
birth weight (kg)41.38 ± 1.3740.95 ± 1.5842.45 ± 1.63
slaughtering weight (kg)100.62 ± 1.94 a100.07 ± 1.92 a90.77 ± 2.42 b
liveweight gain (kg)59.23 ± 7.7759.12 ± 8.4448.32 ± 11.82
weight of rumen (kg)3.35 ± 0.11 a2.35 ± 0.69 b2.30 ± 0.0.68 b
volume of rumen (L)32.62 ± 0.41 a20.48 ± 3.22 b21.00 ± 6.07 b
The ratio of rumen to live weight (%)2.14 ± 0.03 ab2.35 ± 0.43 a1.57 ± 0.17 b
The ratio of rumen to total stomach weight (%)60.95 ± 2.40 a61.13 ± 3.34 a52.38 ± 0.59 b
the pH of the rumen6.63 ± 0.16 a5.33 ± 0.26 b6.30 ± 0.13 a
Note: There was no same lowercase letter in the same row of data shoulder mark to indicate that the difference was significant (p < 0.05), and the same letter or no letter indicated that the difference was not significant (p > 0.05), GF: “Concentrate” feeding pattern, GFF: “Concentrate + hay” feeding pattern, TMR: “Total Mixed Rations” feeding pattern.
Table 2. Composition of rumen microorganisms of calves.
Table 2. Composition of rumen microorganisms of calves.
SampleArchaeaBacteriaEukaryotaVirusesOthers
GF1545,05635,433,043371025,2376,172,506
1.292%84.005%0.009%0.060%14.634%
GF2460,36535,799,595289143,3925,873,309
1.091%84.874%0.007%0.103%13.925%
GF3363,25534,506,6761810172,2467,135,566
0.861%81.809%0.004%0.408%16.917%
GFF1221,64135,363,072249060,7956,531,555
0.525%83.839%0.006%0.144%15.485%
GFF2998,99834,991,380271543,1466,143,314
2.368%82.958%0.006%0.102%14.565%
GFF3368,13735,470,400635055,2136,279,452
0.873%84.094%0.015%0.131%14.887%
TMR176,60436,590,887565231,7715,474,639
0.182%86.750%0.013%0.075%12.979%
TMR2158,47736,270,6491112140,2525,699,055
0.376%85.991%0.026%0.095%13.511%
TMR3215,85035,870,04714,37380,5485,998,736
0.512%85.041%0.034%0.191%14.222%
mean378,70935,588,416567961,4006,145,348
0.898%84.374%0.013%0.146%14.569%
SEM92,075.71211,1401449.2714,898.98162,144.56
0.218%0.501%0.003%0.035%0.384%
Table 3. Relative abundances of the top 10 components of the microbiota at the phylum level.
Table 3. Relative abundances of the top 10 components of the microbiota at the phylum level.
PhylumRelative Abundance (%)MeanSEM
GFGFFTMR
Firmicutes43.6829.6925.6132.990.0323
Bacteroidetes31.1141.3848.8940.460.0302
Proteobacteria10.703.082.4421.960.0047
Actinobacteria2.722.560.4419.100.0053
Chlamydiae0.320.170.190.230.0003
Tenericutes0.560.620.440.540.0008
Euryarchaeota1.081.250.350.890.0022
Spirochaetes0.110.371.130.540.0016
Fibrobacteres0.080.230.680.330.0010
Synergistetes0.070.140.030.080.0002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, J.; Zhao, K.; Li, M.; Fan, H.; Wang, M.; Xia, S.; Chen, Y.; Bai, X.; Liu, Z.; Ni, J.; et al. A Preliminary Study of the Potential Molecular Mechanisms of Individual Growth and Rumen Development in Calves with Different Feeding Patterns. Microorganisms 2023, 11, 2423. https://doi.org/10.3390/microorganisms11102423

AMA Style

Wang J, Zhao K, Li M, Fan H, Wang M, Xia S, Chen Y, Bai X, Liu Z, Ni J, et al. A Preliminary Study of the Potential Molecular Mechanisms of Individual Growth and Rumen Development in Calves with Different Feeding Patterns. Microorganisms. 2023; 11(10):2423. https://doi.org/10.3390/microorganisms11102423

Chicago/Turabian Style

Wang, Jie, Kaisen Zhao, Mianying Li, Huimei Fan, Meigui Wang, Siqi Xia, Yang Chen, Xue Bai, Zheliang Liu, Jiale Ni, and et al. 2023. "A Preliminary Study of the Potential Molecular Mechanisms of Individual Growth and Rumen Development in Calves with Different Feeding Patterns" Microorganisms 11, no. 10: 2423. https://doi.org/10.3390/microorganisms11102423

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

Article Metrics

Back to TopTop