Xiaojing Liu1†
Jie Li2,3†
Lizhuang Hao2*
Allan Degen4Dongyang Wang1Yonggui Ma5Jianzhang Niu2
Yanfen Cheng1,6*
Shujie Liu2- 1Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing, China
- 2Key Laboratory of Plateau Grazing Animal Nutrition and Feed Science of Qinghai Province, Qinghai Academy of Animal Science and Veterinary Medicine of Qinghai University, Xining, China
- 3Gansu Polytechnic College of Animal Husbandry & Engineering, Wuwei, China
- 4Desert Animal Adaptations and Husbandry, Wyler Department of Dryland Agriculture, Blaustein Institutes for Desert Research, Ben-Gurion University of the Negev, Beer Sheva, Israel
- 5Academy of Plateau Science and Sustainability, People’s Government of Qinghai Province and Beijing Normal University, Key Laboratory of Medicinal Animal and Plant Resources of Qinghai-Tibetan Plateau in Qinghai Province, College of Life Science, Qinghai Normal University, Xining, China
- 6State Key Laboratory of Grassland Agro-Ecosystems, Center for Grassland Microbiome, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, China
This study examined the effect of the ratio of dietary metabolizable energy (MJ) to nitrogen (g) content (ME:N) on average daily gain (ADG), blood biochemical indices, rumen fermentation parameters, and rumen bacterial community in yaks. Thirty-six male yaks, aged 2–3 years, were divided into three groups and received a ME:N ratio of 0.42 (HY), 0.36 (MY,) or 0.32 (LY) MJ/g. Dry matter intake ranged between 3.16 and 3.63 kg/d and was lesser (p < 0.001) in the LY group than the other two groups. ME intake increased (p < 0.001) with an increase in the ME:N ratio, while N intake did not differ among groups. The ADG was 660 g/day for the MY group, which was higher (p < 0.005) than the 430 g/day in the LY group, while the HY group gained 560 g/day and did not differ from the other two groups. Feed intake to ADG ratio ranged between 5.95 and 7.95, and numerically was highest in the LY group and lowest in the MY group. In general, the concentration of ruminal total volatile fatty acids (p < 0.03) and molar proportions of propionate (p < 0.04), increased, while the molar proportion of acetate (p < 0.005) and the acetate:propionate ratio decreased (p < 0.001) with a decrease in the ME:N ratio. The molar proportion of butyrate did not differ among groups (p = 0.112). Group MY had higher ruminal NH3-N content than group HY and had a higher serum glucose content but lower urea content, lactate dehydrogenase, and creatine kinase content than group LY. In ruminal bacteria at the phylum level, the relative abundance of Firmicutes (F) was greater and of Bacteroidetes (B) was lesser, while the F:B ratio was greater in group MY than in groups HY an LY. We concluded that the yaks consuming the diet containing a ME:N ratio of 0.36 MJ/g had the best performance of the three groups.
Introduction
Yaks are multifunctional and dominant livestock on the Qinghai-Tibetan plateau (QTP), providing meat, milk, fiber, dung, and transportation (Long et al., 2008). They are vital to the livelihood and play important roles in the cultural and religious traditions of local residents (Joshi et al., 2020). Yaks have adapted well to the harsh climate of the QTP, which is characterized by severe environmental conditions, namely low air temperature and oxygen concentration, high ultraviolet light, and a short growing season of forage (Miao et al., 2015). Traditionally, yaks graze in natural pasture all year round without supplements. Much of the grassland is degraded and forage is sparse and of poor quality in the winter, resulting in considerable body weight loss in the yaks. With the drive for ecological protection and restoring the degraded grassland, yak production is gradually introducing feedlot fattening.
Dietary energy and protein are important factors affecting the digestion and metabolism of nutrients in livestock production (Schroeder and Titgemeyer, 2008). Protein and energy levels are closely related to nutrient degradation, rumen fermentation, and intestinal microbial structure. Appropriate energy-to-nitrogen (ME:N) ratio can improve energy and N utilization efficiency and reduce N losses in urine and feces (Marini and Amburgh, 2001; Yan et al., 2006; Hossein-Zadeh and Ardalan, 2011). Furthermore, it has been reported that the ratio can affect average daily gain (ADG), rumen fermentation parameters, and bacteria diversity in ruminants (Gabler and Heinrichs, 2003).
With the increase in feedlot fattening of yaks, more information is needed on their energy and protein requirements. Often, requirements are adapted from cattle, but these two bovine species are different in their requirements. Moreover, the research on yaks has focused mainly on the effects of different energy levels on their metabolism and rumen microbial diversity (Ahmad et al., 2020; Hu et al., 2020), and more information is needed on optimal dietary energy and protein levels. Therefore, in this study, we offered yaks different ratios of ME:N and determined their dry matter intake (DMI), average daily gain (ADG), blood biochemistry, rumen fermentation parameters, and the diversity and structure of rumen bacteria. The results could provide a basis for optimal ME:N ratios for feedlot-raised yaks.
Materials and methods
The study was conducted at the Animal Nutrition Laboratory, Academy of Animal Husbandry and Veterinary Sciences, Qinghai University, Xining, Qinghai Province, China (101.52°E, 37.4°N, 2998 m above sea level). The site has a plateau continental climate. The protocol and all experimental procedures on the animals were approved by the Committee of Animal Use, Academy of Science and Veterinary Medicine, Qinghai University (Protocol QMKY202002).
Animal feeding and sampling
Thirty-six 2–3 year old male yaks were divided randomly into three groups based on the dietary metabolizable energy:nitrogen content (ME:N), namely 0.42 MJ/g (HY), 0.36 MJ/g (MY), and 0.32 MJ/g (LY). The dietary rations were prepared based on the NRC beef cattle feeding standard and relevant research on the nutritional needs of yaks (Han et al., 1993; Xue et al., 1994). Dietary composition and nutritional level of each diet are presented in Table 1.
Table 1. Dietary components and chemical composition of diets with different ME:N ratios (DM basis).
The yaks were kept in an open-sided roofed structure. They were penned (1.5 m × 3 m) individually, which allowed measurement of individual feed intake, and were fed ad libitum at 08:00 and 17:00 daily, with water freely available. The diet was offered for 28 days, and on day 29, before morning feed, the yaks were weighed, and jugular vein blood and rumen fluid were collected. Blood samples were collected by vacuum tubes, centrifuged at 3000 × g for 20 min and serum samples were stored at −20°C. Approximately 300 ml of rumen fluid was collected using a flexible oral stomach tube with a metal strainer (Anscitech Co. Ltd., Wuhan, China). The first 200 ml of fluid was discarded to eliminate saliva contamination, and the equipment was washed thoroughly with clean warm water between collections. The rumen fluid was filtered through four layers of gauze, pH was measured with a pH meter (EcoScan pH 5, Singapore), and then, the rumen fluid was stored at liquid nitrogen tank for determination of fermentation parameters and microbial extraction.
Determination of the chemical composition of the diets, blood biochemical indices, and rumen fermentation parameters
Feed samples of the three treatment diets were collected in two consecutive days each week (n = 8 for each diet), and stored at −20°C. The samples were dried in an oven at 65°C for 48 h, passed through a 2 mm sieve, and then dried at 105°C for 8 h to determine dry matter content. Gross energy of the diets was measured by an oxygen bomb calorimeter (German IKA C6000 gs, Staufen, Germany), and EE by Soxhlet extraction using petroleum ether (AOAC, 2010). Total N was determined by the Kjeldahl method and crude protein was calculated as N × 6.25. Neutral detergent fiber (NDF) concentration was assayed with heat stable α-amylase and sodium sulfite, and acid detergent fiber (ADF) concentration was determined using fiber bags and a fiber analyzer (ANKOM 200ifiber analyzer; ANKOM Technologies, Inc., Fairport, NY, United States) following Van Soest et al. (1991).
Serum biochemical determination kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China) were used to determine the concentrations of serum glucose (Glucose Assay Kit, F006-1-1), triglycerides (Triglyceride Assay Kit, A110-1-1), alkaline phosphatase (Alkaline Phosphatase Assay Kit, A059-2-2), total protein (Total Protein Quantitative Assay Kit, A045-2-2), albumin (Albumin Assay Kit, A028-2-1), urea (Urea Assay Kit, C013-2-1), creatine kinase (Creatine Kinase Assay Kit, A032-1-1), cholesterol (Total Cholesterol Assay Kit, A111-1-1), aspartate aminotransferase (Aspartate Aminotransferase Assay Kit, C010-2-1), alanine aminotransferase (Alanine Aminotransferase Assay Kit, C009-2-1), and lactate dehydrogenase (Lactate Dehydrogenase Assay Kit, A020-2-2). The colorimetric method was used for the detection of indicators. The concentration of globulin was taken as the difference between total protein and albumin (Jiang et al., 2015).
The concentration of ruminal NH3-N was determined by the colorimetric method (Shen et al., 2017), and of microbial crude protein by Coomassie brilliant blue method (Makkar et al., 1982). The concentrations of volatile fatty acids, acetate, propionate, isobutyrate, butyrate, isovalerate, and valerate were determined by gas chromatography (Daojin GC2014AFsc instrument, Shimadzu, Japan) at 40°C column temperature, 220°C injection temperature, and 230°C TCD temperature as detailed by Wang et al. (2011). The pressures of air, N2 carrier gas, and H2 were 0.05 Mpa. The standard curves were established and crotonic acid was used as the internal standard to calculate the concentration of the measured SCFA.
Microbiome sequencing and bioinformatic analysis
Total DNA from rumen fluid was extracted using Power DNA Isolation Kit (Mo Bio Laboratories, Carlsbad, CA). The DNA quality and concentration were checked by NanoDrop (Thermo Fisher Scientific, Waltham, MA, United States), and the V3-V4 region of 16S rRNA was amplified by PCR. High-throughput sequencing primers were 338-F (5″- ACTCCTACGGGAGGCAGCAG-3″) and 806-R (5”-GGACTACCVGGGTATCTAAT-3″; Klindworth et al., 2013). Agencourt AMPure XP magnetic beads (Beckman Coulter, Milan, Italy) amplified the product into purification, library quality inspection was carried out, and qualified libraries were sequenced using Illumina Miseq platform. After obtaining the original results, cutadapt v2.6 software removed the primer sequence and barcode (Bolger et al., 2014), and QIIME2 was used for quality filtering, double-end stitching and chimera removal (Caporaso et al., 2010). Operational taxonomic units (OTUs) were clustered with a 99% similarity cutoff standard into different characteristic sequences using the Vsearch plug-in unit in QIIME2. The relative abundances of bacteria at the phylum and genus levels were obtained by species annotation using Silva 16S rRNA database SILVA_128. Alpha diversity and Beta diversity were calculated with QIIME 2 software. The principal coordinate analysis (PCoA) was used to compare the bacteria communities among groups based on the unweighted UniFrac distance (Lozupone et al., 2011). The common and special OTUs among groups were analyzed using the Venn diagram. The linear discriminant analysis (LDA) effect size (LEfSe) was used to detect the dominant bacterial community differences among groups (Segata et al., 2011). To predict taxa functions, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway hierarchy level 3 was applied using PICRUSt software (Langille et al., 2013). A Spearman’s correlation heatmap analyzed the correlations between the different phenotypic indices and bacteria. Afterward, the raw sequence data were subjected to the Sequence Read Archive (SRA) of the NCBI under accession number PRJNA835682.
Data processing and analysis
The ME of the three diets was calculated as: DE = GE × [91.6694–91.3359 × (ADF-OM)]; ME = DE*0.82 (Chinese feeding standard NY/T815-2004). One-way ANOVA and LSD for multiple comparisons tested for differences among dietary treatments in dry matter intake (DMI), ADG, blood biochemical indices, and rumen fermentation parameters (SPSS 26.0, Chicago, IL, United States). A level of p < 0.05 was accepted as significant, and results are expressed as mean ± standard error of mean (SEM). Stacked histograms were plotted by Graph Prism v6.1 Bar and box graphs, PCoA, the Venn diagram, LEfSe, and heatmap graphs were plotted using the online website,2,3 in which the LDA threshold of comparative analysis of rumen bacterial abundance was greater than 3, and of bacterial differential function was greater than 2.
Results
Dry matter intake and average daily gain
The DMI ranged between 3.16 and 3.63 kg/d and was lesser in the LY group than the other two groups (p < 0.01). The ME intake increased (p < 0.001) with an increase in the ME:N ratio, while there was no difference among groups in N intake (Table 2). The ADG was 660 g/day for group MY, which was greater (p < 0.005) than the 430 g/day in group LY, while group HY gained 560 g/day and did not differ from the other two groups. Feed intake to ADG ratio ranged between 5.95 and 7.95, and numerically was the highest in group LY and the lowest in group MY.
Table 2. Effect of dietary ME:N ratio on dry matter intake, average daily gain, and feed:gain ratio in yaks.
Rumen fermentation parameters
Rumen pH in HY yaks was higher than in the MY and LY yaks (p < 0.001). Concentration of ruminal NH3-N decreased (p < 0.001) with an increase in the ME:N ratio, while microbial crude protein followed the same trend numerically, albeit the differences were not significant. In general, the concentrations of ruminal total volatile fatty acids (TVFAs; p < 0.03) and molar proportions of propionate (p < 0.04) increased, while the molar proportion of acetate (p < 0.005) and the acetate:propionate ratio decreased (p < 0.001) with a decrease in the ME:N ratio. The molar proportion of butyrate did not differ among groups (p = 0.112; Table 3).
Table 3. Effect of dietary ME:N ratio on rumen fermentation parameters in yaks.
Blood biochemical indices
The concentration of serum glucose in group MY was greater (p < 0.02) than in groups LY and HY; total protein in group MY and group LY was greater (p < 0.04) than in group HY; and, urea, lactate dehydrogenase, and creatine kinase in group LY were greater (p < 0.05) than in group HY and MY (Table 4). Concentrations of serum triglycerides, cholesterol, alkaline phosphatase, aspartate transaminase, alanine aminotransferase, albumin, and globulin did not differ among groups.
Table 4. Effect of dietary ME:N ratio on blood biochemical indices in yaks.
Ruminal bacteria
For the α diversity indices (Table 5; Figure 1), observed-species, Chao1 and Shannon were the lowest (p < 0.001) in group HY, but the Simpson index did not differ among groups.
Table 5. Changes of different ME:N ratios on alpha diversity indices of rumen bacteria in yaks.
Figure 1. Differences in ruminal bacterial diversity and richness among groups. Bacterial richness was estimated by the observed-species and the Chao1 indices (A,B). Bacterial diversity was estimated by Shannon and Simpson indices (C,D). HY, high ME:N (0.42 MJ/g); MY, medium ME:N (0.36 MJ/g); LY, low ME:N (0.32 MJ/g).
Based on Venn diagrams, there were 243 common taxa of bacteria, and 27, 14, and 26 unique taxa for the HY, MY, and LY groups, respectively (Figure 2A). There were also significant differences in rumen bacterial community structure and diversity among groups. For beta diversity, PCoA analysis revealed that the bacterial community structure of rumen fluid had clear structural separation among groups based on the Bray-Curtis index (Figure 2B). Differential bacteria were identified from the kingdom to genus levels using the LEfSe analysis. By setting the LDA value to 3 and analyzing the difference of rumen bacterial communities among groups, the rumen bacterial community responded significantly to the ME:N ratio (Figures 2C,D).
Figure 2. Venn diagrams (A) and PCoA plots (B) of the rumen bacterial community for the different ME:N ratios. LEfSe analysis of the bacterial community for the different ME:N. Cladogram (C) showing significantly enriched bacterial taxa from phylum to genus level. In the histogram (D), the red color illustrates that the abundance of taxa was highest in the high ME:N (0.42 MJ/g); the green color illustrates that the abundance of taxa was highest in the medium ME:N (0.36 MJ/g); and the blue color illustrates that the abundance of taxa was highest in the low ME:N (0.32 MJ/g). HY, high ME:N (0.42 MJ/g); MY, medium ME:N (0.36 MJ/g); LY, low ME:N (0.32 MJ/g).
Rumen bacteria were dominated by the phyla Firmicutes, Bacteroidetes, Proteobacteria, and Tenericutes (Figures 3A–D). The relative abundance of Firmicutes was greater (p < 0.01), and of Bacteroidetes was lesser (p < 0.05), while the ratio of Firmicutes to Bacteroidetes was greater (p < 0.01) in the MY group than the LY and HY groups. The dominant genera were Prevotella 1, which was least in the MY group, Ruminococcaceae_NK4A214_group, which was least in the LY group, Christensenellaceae_R-7_group and Rikenellaceae_RC9_gut_group, which were greatest in the MY group, and Prevotellaceae_UCG-003, which was greatest in the LY group (p < 0.001). In addition, differences among groups emerged for Unidentified, Ruminococcaceae UCG-014, Eubacterium coprostanoligenes group, Ruminococcus 2, Acetitomaculum, Lachnospiraceae XPB1014 group, and Veillonellaceae UCG-001 (Figure 3E; Table 6). Based on Spearman’s correlation analysis, the 10 different genera of rumen bacteria and the 10 phenotypic indicators correlation significantly (| r | > 0.6, p < 0.05; Figure 4). The KEGG function prediction of rumen bacterial community indicated that there were significant differences in rumen bacterial functions among groups (LDA > 2, p < 0.05; Figure 5).
Table 6. Effect of the ME:N ratio on the top 20 bacterial genera in yak rumen fluid.
Figure 3. Changes in the relative abundance of the dominant phylum and genera in the rumen bacteria among of yak groups. These are (A) Firmicutes, (B) Bacteroidetes, and (C) the ratio of Firmicutes to Bacteroidetes. Relative abundances of bacteria at phylum level (D) and genus level (E) for different ME:N ratios. HY, high ME:N (0.42 MJ/g); MY, medium ME:N (0.36 MJ/g); LY, low ME:N(0.32 MJ/g). *p < 0.05 and **p < 0.01.
Figure 4. The correlation heat maps between the production performance index, rumen fermentation, blood biochemical indices, and significant different bacteria.
Figure 5. Analysis of LEfSe difference in KEGG metabolic pathway of rumen bacteria. HY, high ME:N (0.42 MJ/g); MY, medium ME:N (0.36 MJ/g); LY, low ME:N (0.32 MJ/g).
Discussion
Diet is the main factor shaping the rumen bacterial community and determining the utilization rate of available energy and N by rumen microorganisms (Henderson et al., 2015). The appropriate ratio of rumen substrate and N availability is not only conducive to understanding the relationship between rumen fermentation and host supply and demand, but is also vital in optimizing energy and N use efficiency and reducing environmental pollution (Ahmad et al., 2020). The results in the present study demonstrated that the ADG, DMI, blood biochemical indices, rumen fermentation parameters, and rumen bacterial structure of yaks differed when consuming diets of different ME: N ratios.
Average daily weight gain and dry matter intake
In the present study, DMI did not differ between the ME:N ratios of 0.42 and 0.36 MJ/g, but then decreased at 0.32 MJ/g. Reports on the effect of ME:N ratio on DMI have been inconsistent. In a study by Gabler and Hendricks (2003) on Holstein heifers offered diets with an ME:N ratio between 0.54 and 0.34 MJ/g, DMI increased initially and then decreased, with a decrease in the ME:N ratio. In Chinese Holstein heifers consuming diets of ME:N ratio ranging between 0.60 and 0.44 MJ/g, DMI was not affected by the ratio (Dong et al., 2017). Some of the differences among studies could be due to differences in the composition of the diets and in the range of ME:N ratios offered the animals.
The ADG in group MY was greater than in group LY, which was consistent with the findings of Gabler and Hendricks (2003), that is, with a decrease in the ME:N (increase in proportion of N) ratio, ADG first increased and then decreased. In the study of Gabler and Hendricks (2003), Holstein heifers were offered 4 diets ranging in ME:N ratio from 0.58 to 0.34 MJ/g, and the highest ADG was achieved in the group consuming a ratio of 0.41 MJ/g. However, the effect of ME:N ratio on ADG has not been equivocal. In another study by Gabler and Hendricks (2003), Holstein heifers were offered diets differing in ME:N ratio between 0.54 and 0.34 MJ/g, and there was no difference in ADG among groups. Similarly, Dong et al. (2017) reported no difference in ADG among the three groups of Chinese Holstein heifers consuming diets of different ME:N ratios. Wittayakun et al. (2019) examined the effect of the ME:N ratio in Holstein heifers fed a pineapple waste silage-based diet with ratios between 0.65 and 0.46 MJ/g, and found no difference among groups in ADG.
Rumen fermentation parameters
As a large anaerobic fermentation tank, the rumen provides a stable environment for the growth and proliferation of microorganisms (Flint, 1997). The microorganisms help the host digest crude fiber and crude protein and produce ammonia nitrogen (NH3-N), amino acids, volatile fatty acids (VFAs), and other substances, and provide energy for yaks (Shabat et al., 2016). Consequently, the pH and concentrations of NH3-N and VFAs are the main indicators of the status of rumen fermentation (Liang et al., 2019). The VFAs lower the rumen pH, while ammonia increases the rumen pH (Zebeli and Metzler-Zebeli, 2012). In the present study, the TVFAs and the NH3-N increased with a decrease in the ratio of dietary ME:N, that is, they had opposing effects on the rumen pH. Concomitantly, the rumen pH decreased, which indicated that the effect of TVAs was greater than that of NH3-N on the rumen pH. The TVFAs were greater in the LY than in the HY group, which was consistent with the report of Dong et al. (2017), in which the Chinese Holstein heifers increased TVFAs with an increase in the proportion of dietary N. However, Gabler and Hendricks (2003) reported no difference in TVFAs in Holstein heifers consuming diets differing in ME:N ratios.
The concentrations of ruminal NH3-N in the three groups in the present study were all within the optimal concentrations of 5 to 25 mg/dl (Preston and Leng, 1987). However, there was an increase in ruminal NH3-N concentration with a decrease in ME:N ratio, that is, with an increase in dietary N, which is consistent with results in cattle reported by Gabler and Hendricks (2003) and Dong et al. (2017). It is possible that the rate at which microorganisms degraded nitrogen-containing substances and released NH3-N in the rumen exceeded the rate at which microorganisms used NH3-N to synthesize proteins (Mahmoudi-Abyane et al., 2020). This resulted in an accumulation of NH3-N in rumen fluid, and more so in yaks with greater dietary N intake.
In the present study the molar proportion of ruminal acetate decreased, whereas of propionate increased with a decrease in the dietary ME:N ratio. As the ratio of dietary ME:N decreased, the dietary protein content and non-structural carbohydrates increased and of fiber decreased. Production of acetate is linked closely to fiber degradation (Wang et al., 2017), while propionate, a major substrate for gluconeogenesis, increases with a low-fiber diet (Zou et al., 2019). These relationships could explain the decrease in the molar proportions of acetate and increase in propionate with a decrease in the ME:N ratio. However, the molar concentrations of VFAs have been inconsistent when related to the ME:N ratio in cattle, with reports of no effect of the ratio on acetate, propionate, and butyrate in one study (Gabler and Hendricks 2003), but a decrease in these VFAs with an increase in the ME:N ratio in another study (Dong et al., 2017). In the present study, the ratio of acetate to propionate ranged between 3.36 and 2.60. A ratio greater than 2.2 was reported to have a positive effect on the rumen and productive performance of cattle (Russell, 1998), which would indicate that all the yaks had a well-functioning rumen, regardless of the dietary ME:N ratio (Polyorach et al., 2014).
Blood biochemical indices
Urea is the final product of protein and amino acid metabolism, and it was reported that the concentration of serum urea was negatively correlated with nitrogen use efficiency in ruminants (Lewis, 1957). In the present study, the concentration of serum urea was highest at the highest proportion of dietary crude protein to ME, which was also reported by Gabler and Hendricks (2003) and Dong et al. (2017). This yak group also had the highest concentration of ruminal NH3-N, which can be explained by the strong correlation between the concentrations of ruminal NH3-N and serum urea (DePeters and Ferguson, 1992). The lowest serum urea concentration was in the MY group, which indicated an increased nitrogen deposition and protein synthesis (Van Bibber-Krueger et al., 2015), and was consistent with the greatest ADG in this group.
The liver is the central tissue of metabolism in animals and is involved in digestion, absorption, excretion, detoxification, material transportation, energy metabolism, and immunity (Zhang et al., 2019). Aspartate transaminase, which occurs mainly in hepatocytes, reflects the function of liver metabolism, and the serum concentration is an important indicator of the health status of the liver. When protein content in the diet exceeds requirements, it increases the liver load, and cells can be damaged (Shi et al., 2010). The permeability of cell membrane increases, aspartate transaminase in cytoplasm is released into the blood, and the concentration in serum increases. In the present study, the serum aspartate transaminase concentrations in group HY (86 U/l), group MY (85 U/l) and group LY (107 U/l) were within the normal range of 78 to 132 U/l for healthy cattle (Ingvartsen, 2006; Sun et al., 2015). However, there was a trend that the concentration in group LY was higher than in group MY, which suggests that the higher proportion of N (lower ME:N ratio) could cause liver inflammation. Concentration of serum creatine kinase is used as an indicator of gastrointestinal tract disorder (Payne and Payne, 1987). In the present study, the serum concentration of creatine kinase in group LY (363 U/l) was higher than in group MY (284 U/l) and group HY (283 U/l), and exceeded the normal range (14.7–309.8 U/l) reported for healthy water buffalo heifers (Abd Ellah et al., 2014). It is possible that group LY may have gastrointestinal irritations, but this requires further research in the expression of related genes and the structural morphology of the gastrointestinal tract. In addition, it was reported that an increase in the concentration of serum lactate dehydrogenase indicated hepatic lesions (Klein et al., 2020). In the present study, when the dietary ME:N ratio decreased from 0.42 MJ/g to 0.36 MJ/g, the concentration of serum lactate dehydrogenase was not affected, but when the ratio was further decreased to 0.32 MJ/g, the concentration of serum lactate dehydrogenase increased. Therefore, the low ratio of dietary metabolic energy to nitrogen content can potentially cause abnormal concentrations of serum biochemical-related indicators such as aspartate transaminase, creatine kinase, and lactate dehydrogenase, which could indicate liver damage.
Blood glucose is an energy source for tissues and cells in the body. In the present study, the 0.36 ME/N group had the highest serum glucose concentration and the highest ADG. It can be inferred that the increase of serum glucose concentration promoted glucose metabolism, thus improving ADG. High concentration of total protein in serum can improve the body’s metabolic level and immunity, but when the protein intake is insufficient, the serum total protein content decreases (Liu et al., 2015), which is consistent with the results of this study. With the highest ME:N ratio, and, therefore, lowest proportion of dietary N, serum concentration of total proteins decreased. However, the serum concentrations of albumin (A) and globulin (G) proteins did not differ among dietary groups in the present study and were similar to values reported for healthy cows (Alberghina et al., 2011). The A:G ratio range of 0.81 to 0.84 in the present study was within the range of 0.6 to 0.9 (Sun et al., 2015) reported for healthy cows.
Rumen bacterial community
The rumen microbiota in ruminants is affected by dietary intake and the external environment (Weimer, 2015). In the present study, when the dietary ME:N ratio decreased from 0.42 MJ/g to 0.36 MJ/g, the α diversity index in rumen bacteria, including the observed-species, Chao1 and Shannon, increased. The difference in the results of this study may be due to the high content of corn, soybean meal, and bran in the diet. These fermentable carbohydrates were degraded rapidly, resulting in a large amount of organic substrates, and, ultimately increasing rumen bacterial abundance. However, the competition for niches among bacteria may limit the bacteria abundance, which may be the reason why the α diversity from group MY to group LY did not increase. Bacteroidetes and Firmicutes were the dominant ruminal bacterial phyla, as has been reported in buffaloes (Singh et al., 1978), steers (Myer et al., 2015), Tibetan sheep (Fan et al., 2021), and other yaks (Fan et al., 2020). Based on the bacterial functions, group MY was enriched in the TCA cycle pathway compared with groups HY and LY, which is consistent with the concentrations of serum glucose content in group MY. This may be related to the relative abundance of Firmicutes and ratio of Bacteroidetes to Firmicutes, which were greater in group MY than groups HY and LY, and they correlated positively with serum glucose concentration. Bacteroidetes can degrade carbohydrate and protein, utilize polysaccharide, and enhance fattening (Tremaroli and Bäckhed, 2012). Firmicutes are important fibrous degrading bacteria that correlates positively with an increase in body weight (Ley et al., 2006), which is consistent with the greater ADG in group MY than in groups LY and HY. Firmicutes produce large amounts of energy-rich short-chain fatty acids that are associated with digestion efficiency, and are crucial for the energy balance of ruminants (Turnbaugh et al., 2006). In addition, studies have reported that a high ratio of Firmicutes to Bacteroidetes is linked with efficient feed utilization in ruminants (Myer et al., 2015) and this could also explain, at least in part, the greater ADG in the MY than in the LY group. The general trend of increasing the relative abundance of Firmicutes and decreasing the relative abundance of Bacteroidetes suggests that the MY group was promoting fiber digestion.
At the genus level, Prevotella has enzymes and gene clusters necessary for fermentation and utilization of complex polysaccharides, which help to degrade protein and carbohydrate (Xue et al., 2020). This genus requires fibrous tissue or substrates after fiber degradation for optimal growth (McLoughlin et al., 2020). Therefore, the higher relative abundance of Prevotella-1 in the HY group than the MY and LY groups, which may be due to group HY consuming the highest dietary NDF content. However, the reason for the higher relative abundance of Prevotella-1 in group LY than group MY, although it has a lower NDF content, is uncertain and warrants further research. In the current study, there was a positive correlation between Lachnospiraceae XPB1014 group and average daily gain. Lachnospiraceae is important for host metabolites and maintaining a healthy ruminant intestinal environment, with some members exhibiting strong hydrolyzing activities (Louis and Flint, 2009; Sagheddu et al., 2016), which may be part of the reason for the highest ADG in the 0.36 MJ/g group. In addition, the changes of genera between groups are also reflected in Christensenellaceae R-7 group and Rikenellaceae RC9 gut group. The bacteria in the Christensenellaceae family secrete α-arabinoside enzyme, β- galactosidase, and β-glucosidase, which are related to feed efficiency (Perea et al., 2017). Rikenellaceae RC9 plays a role in the degradation of plant derived polysaccharides (He et al., 2015; Peng et al., 2015), which is consistent with the significant positive correlations between Christensenellaceae R-7 and the concentration of serum glucose and between Rikenellaceae RC9 gut group and the concentration of ruminal TVFAs. Ruminococcus is the major cellulolytic bacteria and plays an important role in fiber degradation, being rich in genes encoding cellulase and hemicellulase (Qin et al., 2022), and more abundant in ruminants fed forage-based diets (Henderson et al., 2015). With the decrease of the ratio of ME:N in the present study, the content of dietary NDF decreased, and this may be the reason for the decrease of the abundance of Ruminococcus 2, Ruminococcaceae NK4A214 group, and Ruminococcaceae UCG-014 in the LY group. Ruminococcaceae UCG-014 was correlated positively with the molar concentration of acetate in rumen, which is consistent with the study of Fan et al. (2021). It was reported that the relative abundance of Acetitomaculum was high in ruminants fed a high concentrate diet (Wang et al., 2021), and that the relative abundance increased with increasing dietary energy (Wang et al., 2019a). In the present study, the relative abundance of this genus increased with an increase in dietary energy, as expected, but then decreased with a further increase in dietary energy. The reason for this inconsistent pattern in the relative abundance of Acetitomaculum is uncertain and requires further research. Studies have demonstrated that there are substantial differences in rumen microbial structure and diversity between yaks and cattle (An et al., 2005; Huang et al., 2012; Wang et al., 2019b). Therefore, there are differences in the utilization of nutrient substrates between bovine species and different relationships may be formed between the host and microorganisms. Studies are needed using multi-omics technologies, including metagenomics and proteomics to determine the rumen microbial function of yaks.
Conclusion
Average daily gain (ADG), blood biochemical indices, rumen fermentation parameters, and rumen bacterial structure differed among yaks when the actual dietary ME:N ratios were 0.42 (H), 0.36 (M), and 0.32 (L) MJ/g. The highest ADG was achieved with a ratio of 0.36 MJ/g. The relative abundance of the rumen bacteria phylum Firmicutes (F) was greater and of Bacteroidetes (B) was lesser, while the F:B ratio was greater in group M than in groups H an L. The concentration of serum glucose and the relative abundances of the bacteria genera Christensenellaceae R-7 group and Rikenellaceae RC9 gut group were greatest in this group. The 0.32 group had the highest concentration of ruminal NH3-N, but there were indications of liver and gastrointestinal inflammations. It was concluded that a ME:N ratio of approximately 0.36 MJ/g gave optimal results for growing yaks. However, more studies are needed to refine the optimal ME:N ratio with different feeds and on yaks of different ages, productive stages, and physiological conditions.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: NCBI-PRJNA835682.
Ethics statement
The protocol and all experimental procedures on the animals were approved by the Committee of Animal Use, Academy of Science and Veterinary Medicine, Qinghai University (Protocol QMKY202002). Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
JL, JN, and YM completed the experiment. XL, DW, and AD completed the data analysis. YC, SL, and LH conceived and designed the paper. XL, YC, AD, and LH wrote and revised the article. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by The Second Tibetan Plateau Scientific Expedition and Research Program (2019QZKK0606), Key Laboratory of Plateau Grazing Animal Nutrition and Feed Science of Qinghai Province (2022-ZJ-Y17), and Top Talent project of "Kunlun Talents–High-level Innovation and Entrepreneurship Talents" in Qinghai Province.
Acknowledgments
We thank two reviewers for very helpful suggestions on the paper.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
Abd Ellah, M. R., Hamed, M. I., Ibrahim, D. R., and Rateb, H. Z. (2014). Serum biochemical and haematological reference intervals for water buffalo Bubalus bubalis heifers. J. S. Afr. Vet. Assoc. 85, e1–e7. doi: 10.4102/jsava.v85i1.962
Ahmad, A. A., Yang, C., Zhang, J., Kalwar, Q., Liang, Z., Li, C., et al. (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
Alberghina, D., Giannetto, C., Vazzana, I., Ferrantelli, V., and Piccione, G. (2011). Reference intervals for total protein concentration, serum protein fractions, and albumin/globulin ratios in clinically healthy dairy cows. J. Vet. Diagn. Investig. 23, 111–114. doi: 10.1177/104063871102300119.
An, D., Dong, X., and Dong, Z. (2005). Prokaryote diversity in the rumen of yak (Bos grunniens) and Jinnan cattle (Bos taurus) estimated by 16S rDNA homology analyses. Anaerobe 11, 207–215. doi: 10.1016/j.anaerobe.2005.02.001
AOAC. (2010). Official methods of analysis. 19th Edn. Arlington: Association of Official Analytical Chemists.
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
DePeters, E. J., and Ferguson, J. D. (1992). Nonprotein nitrogen and protein distribution in the milk of cows. J. Dairy Sci. 75, 3192–3209. doi: 10.3168/jds.S0022-0302(92)78085-0
Dong, L. F., Zhang, W. B., Zhang, N. F., Tu, Y., and Diao, Q. Y. (2017). Feeding different dietary protein to energy ratios to Holstein heifers: effects on growth performance, blood metabolites and rumen fermentation parameters. J Anim Physiol Anim Nutr (Berl) 101, 30–37. doi: 10.1111/jpn.12493
Fan, Q., Cui, X., Wang, Z., Chang, S., Wanapat, M., Yan, T., et al. (2021). Rumen microbiota of Tibetan sheep (Ovis aries) adaptation to extremely cold season on the Qinghai-Tibetan plateau. Front Vet Sci 8:673822. doi: 10.3389/fvets.2021.673822
Fan, Q., Wanapat, M., Yan, T., and Hou, F. (2020). Altitude influences microbial diversity and herbage fermentation in the rumen of yaks. BMC Microbiol. 20:370. doi: 10.1186/s12866-020-02054-5
Flint, H. J. (1997). The rumen microbial ecosystem—some recent developments. Trends Microbiol. 5, 483–488. doi: 10.1016/S0966-842X(97)01159-1
Gabler, M. T., and Heinrichs, A. J. (2003). Dietary protein to metabolizable energy ratios on feed efficiency and structural growth of prepubertal Holstein heifers. J. Dairy Sci. 86, 268–274. doi: 10.3168/jds.S0022-0302(03)73605-4
Han, X. T., Hu, L. H., Xie, A. Y., Liu, S. J., and Bi, X. H. (1993). Energy metabolism of growing yaks. Chinese Qinghai J. Animal Veterinary Sci. 1:4.
He, B., Nohara, K., Ajami, N. J., Michalek, R. D., Tian, X., Wong, M., et al. (2015). Transmissible microbial and metabolomic remodeling by soluble dietary fiber improves metabolic homeostasis. Sci. Rep. 5:10604. doi: 10.1038/srep10604
Henderson, G., Cox, F., Ganesh, S., Jonker, A., Young, W., and Janssen, P. H. (2015). Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 5:14567. doi: 10.1038/srep14567
Hossein-Zadeh, N. G., and Ardalan, M. (2011). Estimation of genetic parameters for milk urea nitrogen and its relationship with milk constituents in Iranian Holsteins. Livest. Sci. 135, 274–281. doi: 10.1016/j.livsci.2010.07.020
Hu, R., Zou, H., Wang, H., Wang, Z., Wang, X., Ma, J., et al. (2020). Dietary energy levels affect rumen bacterial populations that influence the intramuscular fat fatty acids of fattening yaks (Bos grunniens). Animals (Basel) 10:1474. doi: 10.3390/ani10091474
Huang, X. D., Tan, H. Y., Long, R., Liang, J. B., and Wright, A. D. (2012). Comparison of methanogen diversity of yak (Bos grunniens) and cattle (Bos taurus) from the Qinghai-Tibetan plateau. China. BMC Microbiol 12:237. doi: 10.1186/1471-2180-12-237
Ingvartsen, K. L. (2006). Feeding- and management-related diseases in the transition cow: physiological adaptations around calving and strategies to reduce feeding-related diseases. Anim. Feed Sci. Technol. 126, 175–213. doi: 10.1016/j.anifeedsci.2005.08.003
Jiang, S. Z., Yang, Z. B., Yang, W. R., Li, Z., Zhang, C. Y., Liu, X. M., et al. (2015). Diets of differentially processed wheat alter ruminal fermentation parameters and microbial populations in beef cattle. J. Anim. Sci. 93, 5378–5385. doi: 10.2527/jas.2015-9547
Joshi, S., Shrestha, L., Bisht, N., Wu, N., Ismail, M., Dorji, T., et al. (2020). Ethnic and cultural diversity amongst yak herding communities in the Asian highlands. Sustain. For. 12:957. doi: 10.3390/su12030957
Klein, R., Nagy, O., Tóthová, C., and Chovanová, F. (2020). Clinical and diagnostic significance of lactate dehydrogenase and its isoenzymes in animals. Vet. Med. Int. 2020:5346483. doi: 10.1155/2020/5346483
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41:e1. doi: 10.1093/nar/gks808
Langille, M. G., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676
Lewis, D. (1957). Blood-urea concentration in relation to protein utilization in the ruminant. J. Agric. Sci. 48, 438–446. doi: 10.1017/S0021859600032962
Ley, R. E., Turnbaugh, P. J., Klein, S., and Gordon, J. I. (2006). Microbial ecology: human gut microbes associated with obesity. Nature 444, 1022–1023. doi: 10.1038/4441022a
Liang, P., Liu, S., Xu, F., Jiang, S., Yan, J., He, Q., et al. (2019). Corrigendum: powdery mildews are characterized by contracted carbohydrate metabolism and diverse effectors to adapt to obligate biotrophic lifestyle. Front. Microbiol. 10:1. doi: 10.3389/fmicb.2019.00001
Liu, Y., Kong, X., Jiang, G., Tan, B., Deng, J., Yang, X., et al. (2015). Effects of dietary protein/energy ratio on growth performance, carcass trait, meat quality, and plasma metabolites in pigs of different genotypes. J Anim Sci Biotechnol 6:36. doi: 10.1186/s40104-015-0036-x
Long, R. J., Ding, L. M., Shang, Z. H., and Guo, X. H. (2008). The yak grazing system on the Qinghai-Tibetan plateau and its status. Rangelands 30, 241–246. doi: 10.1071/RJ08012
Louis, P., and Flint, H. J. (2009). Diversity, metabolism and microbial ecology of butyrate-producing bacteria from the human large intestine. FEMS Microbiol. Lett. 294, 1–8. doi: 10.1111/j.1574-6968.2009.01514.x
Lozupone, C., Lladser, M. E., Knights, D., Stombaugh, J., and Knight, R. (2011). UniFrac: an effective distance metric for microbial community comparison. ISME J. 5, 169–172. doi: 10.1038/ismej.2010.133
Mahmoudi-Abyane, M., Alipour, D., and Moghimi, H. R. (2020). Effects of different sources of nitrogen on performance, relative population of rumen microorganisms, ruminal fermentation and blood parameters in male feedlotting lambs. Animal 14, 1438–1446. doi: 10.1017/s175173111900291x
Makkar, H. P., Sharma, O. P., Dawra, R. K., and Negi, S. S. (1982). Simple determination of microbial protein in rumen liquor. J. Dairy Sci. 65, 2170–2173. doi: 10.3168/jds.S0022-0302(82)82477-6
Marini, J. C., and Amburgh, M. S. (2001). Effect of different levels of protein on nitrogen metabolism of heifers. J. Dairy Sci. 84, 160–161.
McLoughlin, S., Spillane, C., Claffey, N., Smith, P. E., O'Rourke, T., Diskin, M. G., et al. (2020). Rumen microbiome composition is altered in sheep divergent in feed efficiency. Front. Microbiol. 11:1981. doi: 10.3389/fmicb.2020.01981
Miao, F., Guo, Z., Xue, R., Wang, X., and Shen, Y. (2015). Effects of grazing and precipitation on herbage biomass, herbage nutritive value, and yak performance in an alpine meadow on the Qinghai-Tibetan plateau. PLoS One 10:e0127275. doi: 10.1371/journal.pone.0127275
Myer, P. R., Smith, T. P., Wells, J. E., Kuehn, L. A., and Freetly, H. C. (2015). Rumen microbiome from steers differing in feed efficiency. PLoS One 10:e0129174. doi: 10.1371/journal.pone.0129174
Payne, J.M., and Payne, S. (1987). The metabolic profile test. United States: Oxford University Press.
Peng, B., Huang, S., Liu, T., and Geng, A. (2015). Bacterial xylose isomerases from the mammal gut Bacteroidetes cluster function in Saccharomyces cerevisiae for effective xylose fermentation. Microb. Cell Factories 14:70. doi: 10.1186/s12934-015-0253-1
Perea, K., Perz, K., Olivo, S. K., Williams, A., Lachman, M., Ishaq, S. L., et al. (2017). Feed efficiency phenotypes in lambs involve changes in ruminal, colonic, and small-intestine-located microbiota. J. Anim. Sci. 95, 2585–2592. doi: 10.2527/jas.2016.1222
Polyorach, S., Wanapat, M., and Cherdthong, A. (2014). Influence of yeast fermented cassava chip protein (YEFECAP) and roughage to concentrate ratio on ruminal fermentation and microorganisms using in vitro gas production technique. Asian-Australas J Anim Sci 27, 36–45. doi: 10.5713/ajas.2013.13298
Preston, T. R., and Leng, R. A. (1987). Matching ruminant production systems with available resources in the tropics and sub-tropics. Livest. Proc enec 19, 3–4. doi: 10.1016/0301-6226(88)90019-X
Qin, X., Zhang, D., Qiu, X., Zhao, K., Zhang, S., Liu, C., et al. (2022). 2-Hydroxy-4-(Methylthio) butanoic acid isopropyl ester supplementation altered ruminal and cecal bacterial composition and improved growth performance of finishing beef cattle. Front. Nutr. 9:833881. doi: 10.3389/fnut.2022.833881
Russell, J. B. (1998). The importance of pH in the regulation of ruminal acetate to propionate ratio and methane production in vitro. J. Dairy Sci. 81, 3222–3230. doi: 10.3168/jds.S0022-0302(98)75886-2
Sagheddu, V., Patrone, V., Miragoli, F., Puglisi, E., and Morelli, L. (2016). Infant early gut colonization by lachnospiraceae: high frequency of ruminococcus gnavus. Front. Pediatr. 4:57. doi: 10.3389/fped.2016.00057
Schroeder, G. F., and Titgemeyer, E. C. (2008). Interaction between protein and energy supply on protein utilization in growing cattle: a review. Livest. Sci. 114, 1–10. doi: 10.1016/j.livsci.2007.12.008
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Shabat, S. K., Sasson, G., Doron-Faigenboim, A., Durman, T., Yaacoby, S., Berg Miller, M. E., et al. (2016). Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 10, 2958–2972. doi: 10.1038/ismej.2016.62
Shen, J., Liu, Z., Yu, Z., and Zhu, W. (2017). Monensin and Nisin affect rumen fermentation and microbiota differently in vitro. Front. Microbiol. 8:1111. doi: 10.3389/fmicb.2017.01111
Shi, L., Yue, W., Zhang, C., Ren, Y., Zhu, X., Wang, Q., et al. (2010). Effects of maternal and dietary selenium (se-enriched yeast) on oxidative status in testis and apoptosis of germ cells during spermatogenesis of their offspring in goats. Anim. Reprod. Sci. 119, 212–218. doi: 10.1016/j.anireprosci.2010.02.012
Singh, J., Bhati, N. N., Kuila, R. K., and Ranganathan, B. (1978). Activity of streptococcus diacetilactis and its UV- induced mutant in dahi prepared from different types of milk 1. J. Food Prot. 41, 275–276. doi: 10.4315/0362-028X-41.4.275
Sun, Y., Wang, B., Shu, S., Zhang, H., Xu, C., Wu, L., et al. (2015). Critical thresholds of liver function parameters for ketosis prediction in dairy cows using receiver operating characteristic (ROC) analysis. Vet. Q. 35, 159–164. doi: 10.1080/01652176.2015.1028657
Tremaroli, V., and Bäckhed, F. (2012). Functional interactions between the gut microbiota and host metabolism. Nature 489, 242–249. doi: 10.1038/nature11552
Turnbaugh, P. J., Ley, R. E., Mahowald, M. A., Magrini, V., Mardis, E. R., and Gordon, J. I. (2006). An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444, 1027–1031. doi: 10.1038/nature05414
Van Bibber-Krueger, C. L., Miller, K. A., Parsons, G. L., Thompson, L. K., and Drouillard, J. S. (2015). Effects of zilpaterol hydrochloride on growth performance, blood metabolites, and fatty acid profiles of plasma and adipose tissue in finishing steers. J. Anim. Sci. 93, 2419–2427. doi: 10.2527/jas.2014-8771
Van Soest, P. J., Robertson, J. B., and Lewis, B. A. (1991). Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 74, 3583–3597. doi: 10.3168/jds.S0022-0302(91)78551-2
Wang, Y., Cao, P., Wang, L., Zhao, Z., Chen, Y., and Yang, Y. (2017). Bacterial community diversity associated with different levels of dietary nutrition in the rumen of sheep. Appl. Microbiol. Biotechnol. 101, 3717–3728. doi: 10.1007/s00253-017-8144-5
Wang, H., He, Y., Li, H., Wu, F., Qiu, Q., Niu, W., et al. (2019a). Rumen fermentation, intramuscular fat fatty acid profiles and related rumen bacterial populations of Holstein bulls fed diets with different energy levels. Appl. Microbiol. Biotechnol. 103, 4931–4942. doi: 10.1007/s00253-019-09839-3
Wang, H., Li, P., Liu, X., Zhang, C., Lu, Q., Xi, D., et al. (2019b). The composition of fungal communities in the rumen of Gayals (Bos frontalis), yaks (Bos grunniens), and Yunnan and Tibetan yellow cattle (Bos taurs). Pol. J. Microbiol. 68, 505–514. doi: 10.33073/pjm-2019-050
Wang, X. F., Mao, S. Y., Liu, J. H., Zhang, L. L., Cheng, Y. F., Jin, W., et al. (2011). Effect of the gynosaponin on methane production and microbe numbers in a fungus-methanogen co-culture. J. Anim. Feed Sci. 20, 272–284. doi: 10.22358/jafs/66184/2011
Wang, Y., Nan, X., Zhao, Y., Jiang, L., Wang, H., Hua, D., et al. (2021). Dietary supplementation with inulin improves lactation performance and serum lipids by regulating the rumen microbiome and metabolome in dairy cows. Anim Nutr 7, 1189–1204. doi: 10.1016/j.aninu.2021.09.007
Weimer, P. J. (2015). Redundancy, resilience, and host specificity of the ruminal microbiota: implications for engineering improved ruminal fermentations. Front. Microbiol. 6:296. doi: 10.3389/fmicb.2015.00296
Wittayakun, S., Innaree, W., Innaree, S., Chainetr, W., and Kongngoen, N. (2019). Effect of protein to metabolizable energy ratio in pineapple waste silage-based diets on performance of Holstein heifers. J. Anim. Health Prod. 7, 158–165. doi: 10.17582/journal.jahp/2019/7.4.158.165
Xue, B., Chai, S. T., Liu, S. J., Wang, W. B., Xie, A. Y., Hu, L. H., et al. (1994). Study on the protein requirement of growing yak. Chinese Qinghai J. Animal Veterinary Sci. 24:5.
Xue, M. Y., Sun, H. Z., Wu, X. H., Liu, J. X., and Guan, L. L. (2020). Multi-omics reveals that the rumen microbiome and its metabolome together with the host metabolome contribute to individualized dairy cow performance. Microbiome 8:64. doi: 10.1186/s40168-020-00819-8
Yan, T., Frost, J. P., Agnew, R. E., Binnie, R. C., and Mayne, C. S. (2006). Relationships among manure nitrogen output and dietary and animal factors in lactating dairy cows. J. Dairy Sci. 89, 3981–3991. doi: 10.3168/jds.S0022-0302(06)72441-9
Zebeli, Q., and Metzler-Zebeli, B. U. (2012). Interplay between rumen digestive disorders and diet-induced inflammation in dairy cattle. Res. Vet. Sci. 93, 1099–1108. doi: 10.1016/j.rvsc.2012.02.004
Zhang, R. Y., Liu, Y. J., Yin, Y. Y., Jin, W., Mao, S. Y., and Liu, J. H. (2019). Response of rumen microbiota, and metabolic profiles of rumen fluid, liver and serum of goats to high-grain diets. Animal 13, 1855–1864. doi: 10.1017/s1751731118003671
Keywords: yak, ratio of metabolizable energy to nitrogen, average daily gain, blood biochemistry, rumen fermentation, rumen bacteria
Citation: Liu X, Li J, Hao L, Degen A, Wang D, Ma Y, Niu J, Cheng Y and Liu S (2022) Effect of the ratio of dietary metabolizable energy to nitrogen content on production performance, serum metabolites, rumen fermentation parameters, and bacterial diversity in yaks. Front. Microbiol. 13:1013980. doi: 10.3389/fmicb.2022.1013980
Edited by:
Maria Siwek, University of Science and Technology (UTP), PolandCopyright © 2022 Liu, Li, Hao, Degen, Wang, Ma, Niu, Cheng and Liu. 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: Yanfen Cheng, yanfencheng@njau.edu.cn; Lizhuang Hao, lizhuanghao1122@foxmail.com
†These authors share first authorship