Next Article in Journal
HPLC–(Q)-TOF-MS-Based Study of Plasma Metabolic Profile Differences Associated with Age in Pediatric Population Using an Animal Model
Previous Article in Journal
Correlation of Cerebral Microdialysis with Non-Invasive Diffuse Optical Cerebral Hemodynamic Monitoring during Deep Hypothermic Cardiopulmonary Bypass
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 12
Issue 8
10.3390/metabo12080738
metabolites-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

Changes in Blood Metabolic Profiles Reveal the Dietary Deficiencies of Specific Nutrients and Physiological Status of Grazing Yaks during the Cold Season in Qinghai Province of China

by
Jian Gao
1,†,
Deyu Yang
2,†,
Zhanying Sun
1,
Jianzhang Niu
2,
Yuhong Bao
3,
Suozhu Liu
4,
Zhankun Tan
4,
Lizhuang Hao
2,*,
Yanfen Cheng
1,* and
Shujie Liu
2
1
Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing 210095, China
2
Key Laboratory of Plateau Grazing Animal Nutrition and Feed Science of Qinghai Province, Qinghai Plateau Yak Research Center, Qinghai Academy of Animal Science and Veterinary Medicine of Qinghai University, Xining 810016, China
3
Institute of Grassland Science, Tibet Academy of Agricultural and Animal Husbandry Sciences, Lhasa 850000, China
4
College of Animal Science, Tibet Agricultural and Animal Husbandry University, Nyingchi 860000, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Metabolites 2022, 12(8), 738; https://doi.org/10.3390/metabo12080738
Submission received: 29 June 2022 / Revised: 22 July 2022 / Accepted: 9 August 2022 / Published: 11 August 2022
(This article belongs to the Section Animal Metabolism)
Browse Figures
Versions Notes

Abstract

:
This study aimed to investigate the changes in the blood metabolic profiles of grazing yaks during the cold season to reveal their physiological status and seek the nutrients needed to be supplemented. Six castrated yaks (3 years old) with 166.8 kg (standard deviation = 5.3) of liveweight grazed in the Qinghai-Tibetan Plateau were used as experimental animals without supplementary feeding. Blood samples of each animal were collected in October and December 2015, and March 2016 for the analysis of serum biochemicals and metabolome. Results showed serum indices involved in protein metabolism in grazing yaks showed greater differences during the cold season than the metabolisms of energy or minerals. Cold stress in December had minor effects on the serum metabolic profiles of yaks compared with those in October. Yaks in October and December shared seven differential serum metabolites and enrichments of the “arachidonic acid metabolism” and “glycine, serine, and threonine metabolism” pathways compared with those in March caused by the shortage of feeds. Summarily, the nutrient deficiency would be influential on the physiological status of grazing yaks during the cold season, especially on the protein metabolism, which could be improved by supplementary feeds.

Graphical Abstract

1. Introduction

Yaks (Bos grunniens) are special livestock inhabited at the Qinghai-Tibet Plateau in China. These animals are key economic incomes for the local herders for providing fur, meat, and milk. Although yaks can tolerate the harsh conditions of this ecosystem such as hypoxia, intense ultraviolet radiation, and feed deficiency [1], the breeding levels of yaks are still low. Most of them are fully grazed in the grassland of Qinghai-Tibet Plateau [2], because the harsh environment in this plateau could not meet the cultivation of common feed crops for the local livestock. Due to the feeding pattern, yaks suffer nutrient deficiencies and cold stress during the cold season, which directly affects their blood metabolic profiles and performances. Previous research reported that both above-ground biomasses and nutritive values (especially energy and nitrogen) of forages in the Qinghai-Tibet Plateau decreased hugely during the cold season [2,3]. Besides, yaks also suffer cold stress caused by the low temperature during the cold season in the highlands. Previous research reported that extreme cold conditions have profound effects on serum metabolic profiles [4] and productivity [5] in ruminants. Thus, it is not surprising that the shortage of feeds and the cold stress cause the liveweight loss of grazing yaks, even up to 25% of total liveweight without supplementary feeds during the cold season [6].
During long-term evolution, yaks could partly adapt to the harsh environment in the highlands. The characteristics and functions of the rumen microbiome in yaks changed by altitude [7,8] and seasonal dynamics [2,9] to adapt to the forage composition in the highlands. This adaptation contributed to no differences in the ruminal efficiency of microbial protein synthesis under the shortage of feed sources [10] and greater ruminal utilization of lignocellulosic plants compared with other domestic animals [11,12]. Furthermore, yaks have greater tolerance to cold stress than other ruminants. The longer and denser fur covered the body surface, and the smaller and fewer sweat glands in the skin of yaks would abate the harmful effects of cold stress during the cold season [13]. Unfortunately, the adaptation of yaks to the highlands mentioned above could not help them to maintain their liveweight during the cold season [6]. Researchers reported that supplementary feeds can attenuate the harmful effects of nutrient deficiencies on the grazing yaks [14]. However, of the nutrient ingredients, which one is the most important for grazing yaks is still unclear. Thus, it is difficult to avoid the liveweight losses and maintain the health status of grazing yaks by supplementary feeding precisely.
Blood parameters are important indicators of animal health. The study of hematological and biochemical parameters helps in understanding the relationship of blood characteristics to the habitat and the adaptability of the animals to their environment. Changes in blood parameters depend on the species, sex, age, and nutritional and physical condition of animals [15,16,17,18,19]. The investigation of blood parameters in animals represents a crucial tool in the advancement of management conditions, especially concerning the detection of healthy versus infected or stressed animals [20,21,22,23,24]. Blood metabolic profiles including the biochemicals, low-molecular metabolites, cations, and anions are the indicators of the metabolic status in ruminants [25]. Puppel and Kuczyńska [26] reviewed the changes in the blood metabolic profiles of perinatal cows and concluded that the increased concentrations of non-esterified fatty acid and β-hydroxybutyric acid (a ketone body) in the blood are the most commonly used indicators for the negative energy balance of hosts. Yaks showed similar results during starvation with lower concentrations of glucose and triglyceride and greater concentrations of ketone body and non-esterified fatty acid in their serum [27,28]. For the protein status of hosts, blood protein and urea nitrogen are important indices to indicate the protein intake and protein absorption of digestive tracts [29]. Cold stress also affects the carbohydrate metabolism and protein status in ruminants, which could reflect using the biochemical parameters of their blood [4]. The above biochemicals are useful for predicting the requirement of major nutrients of ruminants. Metabolomic approaches are powerful tools to investigate the low-molecule metabolites in biological fluids which could reflect the metabolic changes of animals precisely. Recent studies have investigated the energy and protein metabolism of yaks in response to supplementary feeds during the warm season [30] or the cold season [31] using metabolomics technologies. However, the dynamic changes in the blood metabolome of grazing yaks in response to the cold season are still unclear, which causes difficulties for supplementary nutrients precisely for these yaks.
Based on the above results, we hypothesize that cold stress and dietary deficiencies during the cold season affect the blood metabolic profiles of grazing yaks, which reflect the shortage of specific nutrients for these animals. The purpose of the current study was to investigate the changes in serum biochemical indices including minerals as well as the low-molecule metabolites identified by the untargeted metabolome in grazing yaks during the cold season. Results obtained in this study would help to improve the supplementary feed strategies and the precise nutrition for grazing yaks in the highlands with harsh conditions.

2. Materials and Methods

2.1. Experimental Design and Sampling

Six castrated yaks (3 years old) with an average liveweight of 166.8 kg (standard deviation = 5.3) were used as the experimental animals. This experiment contained 15 d for adaption and 150 d for the trial period from 11 October 2015 to 25 March 2016. Before the beginning of this experiment, yaks were kept free of ecto- and endo-parasites using albendazole. All animals were grazed in the same meadow by local farmers (101°4′46″ E, 37°4′46″ N, Haibei, Qinghai) from 07:00 to 18:00 h daily without any supplementary feeding and were given free access to fresh drinking water. The experimental site is 3160 m above sea level with an average annual temperature of 7.6 °C. During the experimental period, the average temperature was −3 °C.
Blood samples were collected from the jugular vein using the vacutainers for serum collection with coagulants (BD Medical Instrument Co., Ltd., Shanghai, China) and ruminal contents were sampled simultaneously through the esophagus using a tube from each yak on 25 October 2015 (Oct, −1 to 11 °C), 25 December 2015 (Dec, −24 to −4 °C), and 25 March 2016 (Mar, −12 to 4 °C) before grazing in the morning. Serum samples were achieved by centrifuging the blood samples at 2000× g for 10 min using the Pico-21 tabletop centrifuge of ThermoFisher Scientific (China) Co., Ltd. (Shanghai, China). Liveweight of the yaks was obtained using the weighing scale platform (0.5 kg accuracy) purchased from Zhengfeng Electronic Technology Ltd. (Suzhou, Jiangsu province, China). Meanwhile, the plant biomass of the experimental meadow was recorded and forage samples were collected in several areas of the meadow on the same sampling days. Forage samples were stored at −20 °C, while the samples of serum and ruminal content were stored at −80 °C for later analysis.

2.2. Analysis of Forage, Rumen Content, and Serum Indices

The forage samples were air-dried at 65 °C and ground for the analysis of nutrient composition. Dry matter, ether extract, and crude ash in forages were determined using the methods of AOAC [32]. The content of organic matter in forages was the difference between the contents of dry matter minus crude ash. The nitrogen content of forages was analyzed by the Kjeldahl method as described in no. 988.05 of AOAC [32] and crude protein was calculated by nitrogen% × 6.25. Fiber components of forages including neutral detergent fiber and acid detergent fiber were analyzed according to Van Soest et al. [33]. Contents of non-fiber carbohydrates in forages were calculated by the contents of dry matter minus the percentages of crude protein, ash, ether extract, and neutral detergent fiber. The contents of calcium and total phosphorus in forages were analyzed using the methods of AOAC [32]. Ruminal ammonia-N and microbial crude protein were analyzed using the methods of Broderick and Kang [34] and Makkar et al. [35], respectively.
The concentrations of serum biochemical indices were all analyzed using the Roche Cobas® 8000 modular analyzer in Qinghai Provincial People’s Hospital and related kits purchased from Roche Diagnostics (Basel, Switzerland), including indices involved with energy metabolism (glucose, total cholesterol, triglyceride), protein metabolism (total protein, albumin, globulin, urea nitrogen), and metabolic enzymes (alanine transaminase, aspartate aminotransferase, alkaline phosphatase, lactate dehydrogenase). For the determination of mineral elements, serum samples were digested with concentrated HNO3 for a night and then heated at 80 °C for 5 h on a hot plate. Concentrations of mineral elements in the digested serum including K, Na, Ca, Mg, and Fe were determined using the inductively coupled plasma optical emission spectroscopy (iCAP, ThermoFisher Scientific (China) Co., Ltd., Shanghai, China) as per the method of Harrington et al. [36].

2.3. Pretreatment and Identification of Serum Metabolites by GC-TOF/MS

A 100 μL serum sample of each yak was fully mixed with 0.35 mL methanol and 20 μL L-2-chlorophenylalanin (purity ≥ 98%, Shanghai Hengbai Biotech Co., Ltd., Shanghai, China) aqueous solution (1 mg/mL) as the internal standard. After fully mixing, the mixture was then centrifuged at 13,000× g under 4 °C for 15 min using the 5425R microcentrifuge (Eppendorf, Hamburg, Germany). The supernatant was concentrated and incubated with methoxy amination hydrochloride (20 mg/mL in pyridine) at 80 °C for 30 min. Then, the solution was subsequently incubated with bis(trimethylsilyl)trifluoroacetamide containing 1% trimethylchlorosilane (v/v) (Regis Technologies, Inc., Morton Grove, IL, USA) at 70 °C for 2 h to conduct the pre-column derivatization of the serum metabolites. Untargeted serum metabolome was conducted using the Agilent 7890 gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF/MS, LECO, St. Joseph, MI, USA). A DB-5MS capillary column (30 m × 250 μm inner diameter, 0.25 μm film thickness) coated with 5% diphenyl cross-linked with 95% dimethylpolysiloxane (J&W Scientific, Folsom, CA, USA) was used to separate the serum metabolites. Test conditions of the gas chromatograph and mass spectrometer were described in detail by Sun et al. [37].

2.4. Bioinformatic Analysis of Untargeted Serum Metabolome

Raw peaks were processed using Chroma TOF software (version 4.3X) and identified with the LECO-Fiehn Rtx5 database (LECO, St. Joseph, MI, USA) as the methods described by Kind et al. [38]. Peak areas of metabolites in each sample were normalized by the peak area of internal standard (L-2-chlorophenylalanin) and subsequently analyzed using MetaboAnalyst (version 5.0, https://www.metaboanalyst.ca, accessed on 13 April 2022) [39]. Metabolites were removed when they were not annotated or more than two-thirds of each sample were missing values. The method of k-nearest neighbor was used to estimate the missing values and data were normalized by the log transformation for later analysis. Chemometrics analysis of serum metabolites for all groups including the principal component analysis (PCA) and sparse partial least squares-discriminant analysis (PLS-DA) were conducted using MetaboAnalyst (version 5.0). The differential serum metabolites between every two groups were identified based on the variable importance in projection (VIP > 1) from orthogonal PLS-DA and the false discovery rate (FDR) from t-tests (FDR < 0.05). Fold changes were calculated by the area of each metabolite between groups (former month/later month), while pathway impact analysis was conducted with the peak area of differential serum metabolites using the MetaboAnalyst (version 5.0) and the Kyoto encyclopedia of genes and genomes (updated October 2019, https://www.genome.jp/kegg/, accessed on 13 April 2022) [40]. The enriched pathways with pathway impact > 0.1 and FDR < 0.05 were deemed as the key differential pathways between groups. The random forest analysis of differential serum indices or metabolites was conducted using the Wekemo Bioincloud (https://www.bioincloud.tech, accessed on 17 April 2022) and only the results with FDR < 0.05 were shown.

2.5. Statistical Analysis

Nutritional values of forage samples among different sampling days were analyzed by one-way analysis of variance, while serum biochemical indices were analyzed by the one-way analysis of covariance with the ID number of yaks as the covariant using SPSS software (version 26, IBM Corporation, Chicago, IL, USA). Multiple comparisons of Bonferroni correction were used to identify the difference between groups under the limitation of the family-wise error rate. All data were tested for normality before statistical analysis. The results were presented as mean ± standard deviation or mean ± pooled SEM. p-values or FDR < 0.05 were identified as the significant differences among groups.

3. Results

3.1. Nutritional Value of Forages, Yaks’ Liveweight, and Ruminal N Metabolism

During this experiment, the average temperatures in October, December, and March were −4 to 11 °C, −20 to 0 °C, and −9 to 6 °C, respectively. In October and December, the plant biomasses in the experimental meadow were 120.63 g/m2 and 124.53 g/m2 on an air-dried basis, while it was much lower in March (3.32 g/m2). The contents of organic matter and fibrous materials showed no difference in the forages collected in different months of the cold season (p > 0.10, Table 1). However, forages in October had greater content of crude protein than those collected in December (p = 0.010), while forages collected in March had lower content of ether extract than the other two months (p = 0.002). During the cold season, the liveweights of grazing yaks decreased significantly (p = 0.027) and ruminal fluid collected in October had greater concentrations of ammonia-N (p = 0.014) and microbial crude protein (p < 0.001) than those collected in the other two months (Figure 1).

3.2. Serum Biochemical Indices and Mineral Elements of Yaks

Table 2 shows that grazing yaks in October had the greatest serum concentrations of total protein (p < 0.001), globulin (p = 0.001), and activity of aspartate aminotransferase (p = 0.007), while they had the lowest total cholesterol in serum (p < 0.001) compared with those in the other two months. Significant differences were found in the serum albumin of yaks between December and March (p = 0.044), albumin to globulin ratio between October and December (p = 0.012), and activity of lactate dehydrogenase between October and March (p = 0.001). The grazing yaks in March had the lowest serum urea nitrogen compared with other months (p = 0.002). Furthermore, no differences were found in the serum concentrations of glucose, triglyceride, and the activities of alanine transaminase and alkaline phosphatase (p > 0.05, Table 2). For the serum minerals of grazing yaks during the cold season, yaks in October had greater serum sodium than in December (p = 0.038), greater magnesium than in December and March (p = 0.005), and tended to increase serum phosphorus compared with other sampling times (p = 0.053). Yaks in March had the lowest serum potassium than in other months (p < 0.001), while no differences were found in serum calcium and iron (p > 0.05, Table 2).

3.3. Serum Metabolome of Yaks

Raw results of serum samples based on GC-TOF/MS obtained 436 valid peaks. After filtering, 393 peaks were kept and 163 metabolites were annotated effectively. Results of PCA based on GC-TOF/MS metabolic profiles showed no clear separation (Figure 2A), while results of sparse PLS-DA showed a clear separation with a 0.111 error rate among all groups (Figure 2B). Models of orthogonal PLS-DA were used for the analysis of differential between groups. Results showed that the corresponding R2Y and Q2 in the models between each two groups were R2Y = 0.913, Q2 = 0.699 (Oct vs. Dec), R2Y = 0.932, Q2 = 0.697 (Oct vs. Mar), R2Y = 0.962, Q2 = 0.646 (Dec vs. Mar), respectively, suggesting orthogonal PLS-DA models were effective for the identification of differential serum metabolites between each of the two groups.
A total of 17 differential metabolites were obtained in the serum samples based on the VIP obtained from orthogonal PLS-DA (VIP > 1) and the FDR from the student’s t-test (FDR < 0.05, Table 3). Of them, 5 (Oct vs. Dec), 9 (Oct vs. Mar), and 12 (Dec vs. Mar) differential metabolites were identified in the serum of yaks during the cold season, respectively. Compared with December, yaks in October had greater serum concentrations of lyxose (FDR = 0.019), threitol (FDR = 0.027), tyrosine (FDR = 0.027), leucrose (FDR = 0.035), and lower threonic acid (FDR = 0.019). Compared with March, yaks in October had greater serum threitol (FDR = 0.003) and lactic acid (FDR = 0.048), while they had lower serum monostearin (FDR < 0.001), 1-monopalmitin (FDR < 0.001), threonic acid (FDR < 0.001), D-(glycerol 1-phosphate) (FDR = 0.001), arachidonic acid (FDR = 0.004), glycine (FDR = 0.027), and gluconic acid (FDR = 0.033). All 12 differential metabolites of yak serum were lower in December compared with March, including monostearin (FDR < 0.001), 1-monopalmitin (FDR < 0.001), D-(glycerol 1-phosphate) (FDR < 0.001), arachidonic acid (FDR < 0.001), α-ketoglutaric acid (FDR = 0.011), gluconic acid (FDR = 0.012), serine (FDR = 0.019), glycine (FDR = 0.019), phosphate (FDR = 0.023), 3-hydroxybenzoic acid (FDR = 0.033), glucuronic acid (FDR = 0.034), and threonic acid (FDR = 0.047).
The differential metabolites of yak serum were further analyzed by the pathway impact. Figure 3B showed that only the pathway of “phenylalanine, tyrosine, and tryptophan biosynthesis” was enriched between October and December (pathway impact = 0.50, FDR < 0.001). The differential serum metabolites were similar between October to March, and December to March, which shared seven differential metabolites (Figure 3A). The pathways of “arachidonic acid metabolism” (pathway impact = 0.32, FDR < 0.001) and “glycine, serine, and threonine metabolism” (pathway impact = 0.27, FDR = 0.002) were enriched between October and March (Figure 3C), while three pathways were enriched between December and March including “arachidonic acid metabolism” (pathway impact = 0.32, FDR < 0.001), “glycine, serine, and threonine metabolism” (pathway impact = 0.27, FDR = 0.001), and “ascorbate and aldarate metabolism” (pathway impact = 0.25, FDR = 0.002) (Figure 3D). Random forest analysis was conducted to clarify the contribution of serum indices with significant differences and differential serum metabolites (Figure 4). In total, three key serum indices (i.e., potassium, total cholesterol, and urea nitrogen) and six key differential serum metabolites (i.e., threonic acid, threitol, lyxose, monostearin, D-(glycerol 1-phosphate), and glycine) among the grazing yaks during the cold season under the conditions of FDR < 0.01 and mean decrease accuracy ≥ 0.04.

4. Discussion

4.1. Effects of Dietary Deficiencies on the Blood Metabolic Profiles of Grazing Yaks between October and March during the Cold Season

Due to the lack of grain sources and outdated animal husbandry, grazing is the major feeding pattern and the forages in the meadow are the direct feeds for yaks on the Qinghai-Tibet Plateau [6]. However, the harsh environment on this plateau decreases the production and nutritive levels of forages during the long cold season (October to the following May) [2,41]. In the present experiment, the content of ether extract in the forages collected in March was significantly lower than that of the forages collected in October and December, which was in agreement with Guo et al. [2]. However, the contents of organic matter and fiber materials had no difference among sampling days. The reason could be that the forages in the experiment pasture were highly lignified and low-quality, since more than 65% of their composition was neutral detergent fiber. Thus, the seasonal changes caused minor effects on the nutrient compositions of these forages.
In the present experiment, yaks suffered starvation in March compared with October and December because of the shortage in plant biomasses, which could cause the changes in blood biochemicals and metabolites in response to the deficiency of specific nutrients [25]. Serum biochemicals involved in the protein metabolism showed a larger variation than the carbohydrate metabolism between October and March. Both limited intakes of nitrogen and energy had serious impacts on the rumen microbes of yaks, especially the key role of energy on the synthesis of microbial crude protein [42], which could be supported by the significantly lower concentrations of ruminal ammonia-N and microbial crude protein in March than those in October. The blood urea derives from the ruminal ammonia that is transformed in the liver [43], while the blood protein levels of ruminants are affected by the dietary nitrogen intake [29]. The decrease of ruminal ammonia-N and microbial crude protein in March caused by starvation in yaks would affect their protein metabolism. Zhou et al. reported that blood concentrations of urea, total protein, and albumin in grazing yaks decreased without supplementary feeds [31]. Thus, the limitation of both dietary energy and nitrogen would further result in lower serum concentrations of total protein, globulin, and urea-nitrogen of yaks in March in the present experiment, which was in agreement with Zhou et al. [31].
Dietary energy density is positive in the blood concentration of glucose in cattle [44]. In the present experiment, serum glucose was not altered among different months. The reason could be that specific gut bacteria genera altered with seasons would help yaks to digest efficiently the lignocellulosic materials in winter forages compared with other cattle [2,9], which could provide more precursors for gluconeogenesis in the liver. Besides, the serum concentration of triglyceride in yaks showed no differences among months during the cold season in the present results, which was in agreement with Ding et al. [45]. Previous research reported that serum concentration of total cholesterol in yaks during the cold season was much lower than that during the warm season, suggesting lower serum cholesterol might indicate the shortage of intakes of dry matter and energy in yaks [46,47]. However, in the present experiment, yaks had lower serum concentrations of total cholesterol in October than those in December and March, which was not in agreement with the above results. The reason could be that the nutrient composition of low-quality forages in the pasture changed slightly in the present experiment. In the present experiment, dietary deficiencies resulted in the significant decrease of serum potassium and magnesium, and the decreased tendency of serum calcium and phosphorus in grazing yaks between October and March. Zhou et al. reported that the content of potassium in pasture forages of the Qinghai showed no difference between the withered (December to April of next year) and withering periods (September to November), while magnesium levels were significantly higher in the forages during the withering period [48]. In addition, the increased calcium content and similar phosphorus content were observed in the forages collected in March compared with those collected in October. Thus, shortage of feeds or low absorption of these minerals would be the reason for the decrease of them in grazing yaks during the cold season.
For serum differential metabolites, gluconic acid (gluconate) is a metabolite in the pentose-phosphate pathway. Previous research reported that the plasma gluconate was increased by long-term caloric restriction [49]. In the present experiment, the increased serum gluconic acid in March could be an indicator to identify yaks suffering from starvation or dietary deficiencies of energy and protein in March. Furthermore, monostearin and 1-monopalmitin are monoglycerides and their concentrations are negatively correlated to blood glucose [50]. Thus, the greater concentrations of these two monoglycerides in the serum of yaks collected in March than in October and December could be attributed to the serum glucose among different months. Results also showed that the greater serum arachidonic acid and glycine of yaks in March affected the pathways of “arachidonic acid metabolism” and “glycine, serine, and threonine metabolism” according to the analysis of pathway impact. Arachidonic acid is an essential fatty acid obtained from foods and its concentration in blood is affected by the intake of arachidonic acid [51]. Thus, the different content of arachidonic acid in forages among months could be the reason for the greater serum arachidonic acid of yaks in March. Besides, glycine is a proteinogenic amino acid that is involved in the body protein turnover [52] and conditionally affects the regeneration efficiency of skeletal muscle regeneration in mammals [53]. Serum glycine of grazing yaks increased in March compared with the other two months in the present experiment. The reason could be that dietary energy or protein deficiency caused the muscle breakdown in these yaks and subsequently increased glycine turnover in the body, which was in agreement with the results in cows suffering malnutrition during early lactation [54]. Results of the random forest showed that the concentrations of total protein, urea nitrogen, and glycine were key metabolic profiles in yaks between October and March, suggesting that dietary deficiency of nutrients caused by starvation would be the major detrimental factor for the health status of yaks during the cold season.

4.2. Effects of Cold Stress on the Blood Metabolic Profiles of Grazing Yaks between October and December during the Cold Season

Grazing yaks also suffer cold stress under extremely low temperatures in winter. For example, the temperatures in December were much lower than those in October in the present experiment (−4 to 11 °C vs. −20 to 0 °C) with no difference in the plant biomass of pasture for the yaks. Thus, the differences in blood metabolic profiles of yaks between October and December would be attributed to the cold stress mainly. A previous study reported that cold exposure increased the dry matter intake and decreased the apparent digestibility of dry matter, organic matter and nitrogen, as well as the concentration of ruminal ammonia in sheep [55]. These results were in agreement with our results that the concentrations of ruminal ammonia and microbial crude protein in yaks were both lower in December than in October, suggesting a decreased nitrogen digestion of them. Nazifi et al. [4] reported that the cold stress (4 °C) caused the decrease of urea nitrogen, triglyceride, activities of alanine aminotransferase and alkaline phosphatase, and the increase of total protein, glucose, cholesterol, and total lipid in sheep blood compared to those under the optimum temperature (21 °C). In the present experiment, yaks that suffered cold stress in December had lower serum concentrations of total protein, globulin, and activity of aspartate aminotransferase, but greater serum cholesterol than those in October. These results were contradictory to Nazifi et al. [4]. Serum activity of aspartate aminotransferase is positive with the intensity of metabolic changes in the hosts, especially protein metabolism [56]. The acute decrease of feeds for grazing yaks from September to November [2] could cause negative energy balance and muscle damage, and subsequently increase the activities of lactate dehydrogenase and aspartate aminotransferase on 25 October rather than on 25 December or 25 March. However, the actual mechanism still needs to be investigated. Cold stress would also affect the absorption efficiency of minerals in mammals, such as the depressed absorption of calcium in young turkeys [57], which could be the reason for the decreased serum magnesium level and the similar calcium level of grazing yaks in December compared with October of the present experiment, even though the calcium content in forages significantly increased in December.
For the small metabolites, the results identified only five differential metabolites in the serum of yaks between October and December, suggesting yaks would have a stronger adaption to the cold environment than other ruminants due to their long-term evolution. Based on these differential metabolites, the pathway “phenylalanine, tyrosine, and tryptophan biosynthesis” of yaks was greatly inhibited in December compared with that in October, due to the decreased serum concentration of tyrosine in December. Previous research reported that cold stress would raise the activity of tyrosine hydroxylase, which is the rate-limiting enzyme for the synthesis of catecholamines including epinephrine, norepinephrine, and dopamine in different tissues of mammals [58]. Thus, the decreased serum tyrosine and inhibition of the pathway “phenylalanine, tyrosine, and tryptophan biosynthesis” could be attributed to the increased activity of tyrosine hydroxylase in yaks’ tissues caused by cold stress in the present experiment. Tyrosine supplementation would improve the ability of mammals to maintain their core temperature [59]. Though serum metabolites of grazing yaks were slightly changed by the cold stress, supplementary tyrosine could be beneficial for them to resist the harsh environment during the cold season. According to the results of random forest, the key metabolic profiles of yaks’ serum (FDR < 0.05, mean decrease accuracy > 0.04) did not have obvious differences between October and December. This suggests that yaks would tolerate the cold stress from −4~11 °C in October to −20~0 °C in December, which could be attributed to their physiological structure and adaptation to the harsh environment in the highlands. However, the dietary deficiencies of nutrients especially the energy and protein had greater effects on the metabolic profiles in the serum of yaks as well as their metabolic status during the cold season. It is worth noting that the percentage of crude protein in the loss of empty body weight of grazing yaks during the cold season was much higher than that of fat (44.0% vs. 24.3%) [6], suggesting the protein degradation in the body would be more hugely than fat degradation to compensate for the lack of energy during this period. Supplementing energy feeds (barley) would be more useful to prevent the weight loss of grazing yaks during the cold season than supplementing the protein source (rapeseed meal) [31].

5. Conclusions

Summarily, dietary nutrient deficiency would be more influential on the physiological status of grazing yaks than the cold stress during the cold season, since the metabolic profiles related to the protein and amino acid metabolism showed great changes in the serum of grazing yaks. The metabolism of protein and amino acids of grazing yaks was affected more greatly by nutrient deficiency than the metabolism of energy and minerals during the cold season based on the analysis of the random forest. Supplementary feeding would be beneficial for improving the metabolic status of grazing yaks during the cold season, but it still needs to be further investigated.

Author Contributions

Conceptualization, L.H., S.L. (Shujie Liu) and Y.C.; methodology, L.H. and Y.B.; software, J.G.; validation, D.Y. and J.N.; formal analysis, D.Y. and L.H.; investigation, D.Y. and J.N.; resources, L.H., S.L. (Suozhu Liu), and Z.T.; data curation, D.Y., J.G., L.H., and Y.C.; writing—original draft preparation, J.G. and Y.C.; writing—review and editing, J.G. and L.H.; visualization, J.G.; supervision, S.L. (Shujie Liu); project administration, L.H. and Z.S.; funding acquisition, L.H. and S.L. (Shujie Liu). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Special Topics of the Second Comprehensive Scientific Expedition of the Qinghai-Tibet Plateau, grant number 2019QZKK0606; National Natural Science Foundation of China, grant number 32060766; Qinghai Province Key Laboratory of Animal Nutrition and Feed Science for Plateau Grazing Livestock, grant number 2022-ZJ-Y17; Qinghai Province “Kunlun Talents · High-end Innovation and Entrepreneurial Talents” Top Talent Project (2020).

Institutional Review Board Statement

The study was approved and supervised by Committee of Animal Use, Academy of Animal Science and Veterinary Medicine, Qinghai University (approval number QHUA-201508).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data that are reported in the present experiment were not deposited in an official repository, but are available from the authors upon request.

Acknowledgments

The authors thank Shanghai Biotree Biotech Co., Ltd. for the technical help with the analysis of serum untargeted metabolome.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PCA: Principal component analysis; PLS-DA, partial least squares-discriminant analysis; VIP, Variable important in projection; FDR, false discovery rate; GC-TOF/MS, gas chromatograph coupled with time-of-flight mass spectrometer.

References

  1. Miao, F.; Guo, Z.; Xue, R.; Wang, X.; Shen, Y. 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 2015, 10, e0127275. [Google Scholar] [CrossRef]
  2. Guo, N.; Wu, Q.; Shi, F.; Niu, J.; Zhang, T.; Degen, A.A.; Fang, Q.; Ding, L.; Shang, Z.; Zhang, Z. Seasonal dynamics of diet–gut microbiota interaction in adaptation of yaks to life at high altitude. npj Biofilms Microbiomes 2021, 7, 38. [Google Scholar] [CrossRef] [PubMed]
  3. Long, R.; Apori, S.; Castro, F.; Ørskov, E. Feed value of native forages of the Tibetan Plateau of China. Anim. Feed Sci. Technol. 1999, 80, 101–113. [Google Scholar] [CrossRef]
  4. Nazifi, S.; Saeb, M.; Rowghani, E.; Kaveh, K. The influences of thermal stress on serum biochemical parameters of Iranian fat-tailed sheep and their correlation with triiodothyronine (T3), thyroxine (T4) and cortisol concentrations. Comp. Clin. Pathol. 2003, 12, 135–139. [Google Scholar] [CrossRef]
  5. Young, B. Ruminant cold stress: Effect on production. J. Anim. Sci. 1983, 57, 1601–1607. [Google Scholar] [CrossRef] [PubMed]
  6. Xue, B.; Zhao, X.Q.; Zhang, Y.S. Seasonal changes in weight and body composition of yak grazing on alpine-meadow grassland in the Qinghai-Tibetan plateau of China. J. Anim. Sci. 2005, 83, 1908–1913. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, Z.; Xu, D.; Wang, L.; Hao, J.; Wang, J.; Zhou, X.; Wang, W.; Qiu, Q.; Huang, X.; Zhou, J. Convergent evolution of rumen microbiomes in high-altitude mammals. Curr. Biol. 2016, 26, 1873–1879. [Google Scholar] [CrossRef]
  8. Wu, D.; Vinitchaikul, P.; Deng, M.; Zhang, G.; Sun, L.; Wang, H.; Gou, X.; Mao, H.; Yang, S. Exploration of the effects of altitude change on bacteria and fungi in the rumen of yak (Bos grunniens). Arch. Microbiol. 2021, 203, 835–846. [Google Scholar] [CrossRef]
  9. Ma, L.; Xu, S.; Liu, H.; Xu, T.; Hu, L.; Zhao, N.; Han, X.; Zhang, X. Yak rumen microbial diversity at different forage growth stages of an alpine meadow on the Qinghai-Tibet Plateau. PeerJ 2019, 7, e7645. [Google Scholar] [CrossRef]
  10. Shi, F.; Guo, N.; Degen, A.; Niu, J.; Wei, H.; Jing, X.; Ding, L.; Shang, Z.; Long, R. Effects of level of feed intake and season on digestibility of dietary components, efficiency of microbial protein synthesis, rumen fermentation and ruminal microbiota in yaks. Anim. Feed Sci. Technol. 2020, 259, 114359. [Google Scholar] [CrossRef]
  11. Dai, X.; Zhu, Y.; Luo, Y.; Song, L.; Liu, D.; Liu, L.; Chen, F.; Wang, M.; Li, J.; Zeng, X. Metagenomic insights into the fibrolytic microbiome in yak rumen. PLoS ONE 2012, 7, e40430. [Google Scholar] [CrossRef]
  12. Jing, X.; Ding, L.; Zhou, J.; Huang, X.; Degen, A.; Long, R. The adaptive strategies of yaks to live in the Asian highlands. Anim. Nutr. 2022, 9, 249–258. [Google Scholar] [CrossRef]
  13. Das, P.; Ranjan, R.; Paul, S. A comparative histological study on the sweat gland of cattle (B. indicus) and yak (P. poephagus). Explor. Anim. Med. Res. 2014, 4, 183–187. [Google Scholar]
  14. Long, R.; Dong, S.; Wei, X.; Pu, X. The effect of supplementary feeds on the bodyweight of yaks in cold season. Livest. Prod. Sci. 2005, 93, 197–204. [Google Scholar] [CrossRef]
  15. Fiore, E.; Giambelluca, S.; Morgante, M.; Piccione, G.; Vazzana, I.; Contiero, B.; Orefice, T.; Arfuso, F.; Gianesella, M. Changes in thyroid hormones levels and metabolism in dairy cows around calving. Acta Vet. 2017, 67, 318–330. [Google Scholar] [CrossRef]
  16. Arfuso, F.; Fazio, F.; Levanti, M.; Rizzo, M.; Di Pietro, S.; Giudice, E.; Piccione, G. Lipid and lipoprotein profile changes in dairy cows in response to late pregnancy and the early postpartum period. Arch. Anim. Breed. 2016, 59, 429–434. [Google Scholar] [CrossRef]
  17. Arfuso, F.; Rizzo, M.; Giannetto, C.; Giudice, E.; Fazio, F.; Piccione, G. Age-related changes of serum mitochondrial uncoupling 1, rumen and rectal temperature in goats. J. Therm. Biol. 2016, 59, 47–51. [Google Scholar] [CrossRef]
  18. Arfuso, F.; Fazio, F.; Rizzo, M.; Marafioti, S.; Zanghì, E.; Piccione, G. Factors affecting the hematological parameters in different goat breeds from Italy. Ann. Anim. Sci. 2016, 16, 743–757. [Google Scholar] [CrossRef]
  19. Fiore, E.; Arfuso, F.; Gianesella, M.; Vecchio, D.; Morgante, M.; Mazzotta, E.; Badon, T.; Rossi, P.; Bedin, S.; Piccione, G. Metabolic and hormonal adaptation in Bubalus bubalis around calving and early lactation. PLoS ONE 2018, 13, e0193803. [Google Scholar] [CrossRef]
  20. Carcangiu, V.; Arfuso, F.; Luridiana, S.; Giannetto, C.; Rizzo, M.; Bini, P.P.; Piccione, G. Relationship between different livestock managements and stress response in dairy ewes. Arch. Anim. Breed. 2018, 61, 37–41. [Google Scholar] [CrossRef]
  21. Arfuso, F.; Zumbo, A.; Castronovo, C.; Giudice, E.; Piccione, G.; Monteverde, V.; Giannetto, C. The housing system influences daily total locomotor activity (TLA) in dairy cows. Biol. Rhythm Res. 2022, in press. [Google Scholar] [CrossRef]
  22. Perillo, L.; Arfuso, F.; Piccione, G.; Dara, S.; Tropia, E.; Cascone, G.; Licitra, F.; Monteverde, V. Quantification of some heavy metals in hair of dairy cows housed in different areas from Sicily as a bioindicator of environmental exposure—A preliminary study. Animals 2021, 11, 2268. [Google Scholar] [CrossRef]
  23. Arfuso, F.; Fazio, F.; Chikhi, L.; Aymond, G.; Piccione, G.; Giannetto, C. Acute stress response of sheep to shearing procedures: Dynamic change of cortisol concentration and protein electrophoretic pattern. Animals 2022, 12, 862. [Google Scholar] [CrossRef]
  24. Arfuso, F.; Acri, G.; Piccione, G.; Sansotta, C.; Fazio, F.; Giudice, E.; Giannetto, C. Eye surface infrared thermography usefulness as a noninvasive method of measuring stress response in sheep during shearing: Correlations with serum cortisol and rectal temperature values. Physiol. Behav. 2022, 250, 113781. [Google Scholar] [CrossRef]
  25. Zaitsev, S.Y.; Bogolyubova, N.V.; Zhang, X.; Brenig, B. Biochemical parameters, dynamic tensiometry and circulating nucleic acids for cattle blood analysis: A review. PeerJ 2020, 8, e8997. [Google Scholar] [CrossRef]
  26. Puppel, K.; Kuczyńska, B. Metabolic profiles of cow’s blood; a review. J. Sci. Food Agric. 2016, 96, 4321–4328. [Google Scholar] [CrossRef]
  27. Zou, H.; Hu, R.; Wang, Z.; Shah, A.M.; Zeng, S.; Peng, Q.; Xue, B.; Wang, L.; Zhang, X.; Wang, X. Effects of nutritional deprivation and re-alimentation on the feed efficiency, blood biochemistry, and rumen microflora in yaks (Bos grunniens). Animals 2019, 9, 807. [Google Scholar] [CrossRef]
  28. Hu, R.; Zou, H.; Wang, H.; Wang, Z.; Wang, X.; Ma, J.; Shah, A.M.; Peng, Q.; Xue, B.; Wang, L. Dietary energy levels affect rumen bacterial populations that influence the intramuscular fat fatty acids of fattening yaks (Bos grunniens). Animals 2020, 10, 1474. [Google Scholar] [CrossRef]
  29. Law, R.; Young, F.; Patterson, D.; Kilpatrick, D.; Wylie, A.; Mayne, C. Effect of dietary protein content on animal production and blood metabolites of dairy cows during lactation. J. Dairy Sci. 2009, 92, 1001–1012. [Google Scholar] [CrossRef]
  30. Xue, B.; Zhang, J.; Wang, Z.; Wang, L.; Peng, Q.; Da, L.; Bao, S.; Kong, X.; Xue, B. Metabolism response of grazing yak to dietary concentrate supplementation in warm season. Animal 2021, 15, 100175. [Google Scholar] [CrossRef]
  31. Zhou, J.; Yue, S.; Peng, Q.; Wang, L.; Wang, Z.; Xue, B. Metabonomic responses of grazing yak to different concentrate supplementations in cold season. Animals 2020, 10, 1595. [Google Scholar] [CrossRef] [PubMed]
  32. AOAC. Official Methods of Analysis, 19th ed.; Association of Official Analytical Chemists: Washington, DC, USA, 2012. [Google Scholar]
  33. Van Soest, P.J.; Robertson, J.B.; Lewis, B.A. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 1991, 74, 3583–3597. [Google Scholar] [CrossRef]
  34. Broderick, G.A.; Kang, J.H. Automated simultaneous determination of ammonia and total amino acids in ruminal fluid and in vitro media. J. Dairy Sci. 1980, 63, 64–75. [Google Scholar] [CrossRef]
  35. Makkar, H.; Sharma, O.; Dawra, R.; Negi, S. Simple determination of microbial protein in rumen liquor. J. Dairy Sci. 1982, 65, 2170–2173. [Google Scholar] [CrossRef]
  36. Harrington, J.M.; Young, D.J.; Essader, A.S.; Sumner, S.J.; Levine, K.E. Analysis of human serum and whole blood for mineral content by ICP-MS and ICP-OES: Development of a mineralomics method. Biol. Trace Elem. Res. 2014, 160, 132–142. [Google Scholar] [CrossRef]
  37. Sun, H.Z.; Wang, D.M.; Wang, B.; Wang, J.K.; Liu, H.Y.; Guan, L.L.; Liu, J.X. Metabolomics of four biofluids from dairy cows: Potential biomarkers for milk production and quality. J. Proteome Res. 2015, 14, 1287–1298. [Google Scholar] [CrossRef]
  38. Kind, T.; Wohlgemuth, G.; Lee, D.Y.; Lu, Y.; Palazoglu, M.; Shahbaz, S.; Fiehn, O. FiehnLib: Mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal. Chem. 2009, 81, 10038–10048. [Google Scholar] [CrossRef]
  39. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef]
  40. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  41. Long, R.; Ding, L.; Shang, Z.; Guo, X. The yak grazing system on the Qinghai-Tibetan plateau and its status. Rangel. J. 2008, 30, 241–246. [Google Scholar] [CrossRef]
  42. Lu, Z.; Xu, Z.; Shen, Z.; Tian, Y.; Shen, H. Dietary energy level promotes rumen microbial protein synthesis by improving the energy productivity of the ruminal microbiome. Front. Microbiol. 2019, 10, 847. [Google Scholar] [CrossRef]
  43. Bach, A.; Calsamiglia, S.; Stern, M.D. Nitrogen metabolism in the rumen. J. Dairy Sci. 2005, 88, E9–E21. [Google Scholar] [CrossRef]
  44. Zhang, Q.; Su, H.; Wang, F.; Cao, Z.; Li, S. Effects of energy density in close-up diets and postpartum supplementation of extruded full-fat soybean on lactation performance and metabolic and hormonal status of dairy cows. J. Dairy Sci. 2015, 98, 7115–7130. [Google Scholar] [CrossRef]
  45. Ding, X.; Guo, X.; Yan, P.; Liang, C.; Bao, P.; Chu, M. Seasonal and nutrients intake regulation of lipoprotein lipase (LPL) activity in grazing yak (Bos grunniens) in the Alpine Regions around Qinghai Lake. Livest. Sci. 2012, 143, 29–34. [Google Scholar] [CrossRef]
  46. Kumagai, H.; Nakajima, M.; Anzai, H.; Sakai, T.; Oishi, K.; Hirooka, H.; Shah, M.K. Health and mineral nutrition status of yaks in southern Mustang, Nepal. Anim. Sci. J. 2017, 88, 1156–1161. [Google Scholar] [CrossRef]
  47. Zheng, J.; Du, M.; Zhang, J.B.; Yi, L.Z.; Ahmad, A.A.; Shen, J.; Hosseini Salekdeh, G.; Ding, X. Transcriptomic and metabolomic analyses reveals inhibition of hepatic adipogenesis and fat catabolism in yak for adaptation to forage shortage during cold season. Front. Cell Dev. Biol. 2022, 9, 759521. [Google Scholar] [CrossRef]
  48. Zhou, Y.; Hao, L.; Liu, S. Milk yield and milk mineral element contents of grazing yaks in different phenological periods in alpine pastures in the source area of three rivers. Chin. J. Anim. Nutr. 2020, 32, 3742–3749. (In Chinese) [Google Scholar] [CrossRef]
  49. Rezzi, S.; Martin, F.-P.J.; Shanmuganayagam, D.; Colman, R.J.; Nicholson, J.K.; Weindruch, R. Metabolic shifts due to long-term caloric restriction revealed in nonhuman primates. Exp. Gerontol. 2009, 44, 356–362. [Google Scholar] [CrossRef]
  50. Zeng, M.; Xiao, Y.; Liang, Y.; Wang, B.; Chen, X.; Cao, D.; Li, H.; Wang, M.; Zhou, Z. Metabolic alterations of impaired fasting glucose by GC/MS based plasma metabolic profiling combined with chemometrics. Metabolomics 2010, 6, 303–311. [Google Scholar] [CrossRef]
  51. Nelson, G.J.; Schmidt, P.C.; Bartolini, G.; Kelley, D.S.; Phinney, S.D.; Kyle, D.; Silbermann, S.; Schaefer, E.J. The effect of dietary arachidonic acid on plasma lipoprotein distributions, apoproteins, blood lipid levels, and tissue fatty acid composition in humans. Lipids 1997, 32, 427–433. [Google Scholar] [CrossRef]
  52. Matthews, D.; Conway, J.; Young, V.; Bier, D. Glycine nitrogen metabolism in man. Metabolism 1981, 30, 886–893. [Google Scholar] [CrossRef]
  53. Gheller, B.J.; Blum, J.E.; Lim, E.W.; Handzlik, M.K.; Fong, E.H.H.; Ko, A.C.; Khanna, S.; Gheller, M.E.; Bender, E.L.; Alexander, M.S. Extracellular serine and glycine are required for mouse and human skeletal muscle stem and progenitor cell function. Mol. Metab. 2021, 43, 101106. [Google Scholar] [CrossRef]
  54. Xu, W.; Vervoort, J.; Saccenti, E.; Kemp, B.; van Hoeij, R.J.; van Knegsel, A.T. Relationship between energy balance and metabolic profiles in plasma and milk of dairy cows in early lactation. J. Dairy Sci. 2020, 103, 4795–4805. [Google Scholar] [CrossRef]
  55. Kennedy, P.; Milligan, L. Effects of cold exposure on digestion, microbial synthesis and nitrogen transformations in sheep. Br. J. Nutr. 1978, 39, 105–117. [Google Scholar] [CrossRef]
  56. Sakowski, T.; Kuczyńska, B.; Puppel, K.; Metera, E.; Słoniewski, K.; Barszczewski, J. Relationships between physiological indicators in blood, and their yield, as well as chemical composition of milk obtained from organic dairy cows. J. Sci. Food Agric. 2012, 92, 2905–2912. [Google Scholar] [CrossRef]
  57. Wolfenson, D.; Sklan, D.; Graber, Y.; Kedar, O.; Bengal, I.; Hurwitz, S. Absorption of protein, fatty acids and minerals in young turkeys under heat and cold stress. Br. Poult. Sci. 1987, 28, 739–742. [Google Scholar] [CrossRef]
  58. Talas, Z.S.; Yurekli, M. The effects of enalapril maleate and cold stress exposure on tyrosine hydroxylase activity in some rat tissues. Cell Biochem. Funct. 2006, 24, 537–540. [Google Scholar] [CrossRef]
  59. Lang, J.A.; Krajek, A.C.; Schwartz, K.S.; Rand, J.E. Oral L-tyrosine supplementation improves core temperature maintenance in older adults. Med. Sci. Sports Exerc. 2020, 52, 928–934. [Google Scholar] [CrossRef]
Metabolites 12 00738 g001 550
Figure 1. Dynamic changes in (A) liveweights, concentrations of (B) ammonia-N and (C) microbial crude protein in the ruminal fluid of grazing yaks during the cold season. Values with different letters represent significant differences among groups (p < 0.05). Results are presented as mean ± standard deviation.
Figure 1. Dynamic changes in (A) liveweights, concentrations of (B) ammonia-N and (C) microbial crude protein in the ruminal fluid of grazing yaks during the cold season. Values with different letters represent significant differences among groups (p < 0.05). Results are presented as mean ± standard deviation.
Metabolites 12 00738 g001
Metabolites 12 00738 g002 550
Figure 2. (A) Principal component analysis (PCA) and (B) sparse partial least squares-discriminant analysis (sPLS-DA) of serum metabolites in grazing yaks during the cold season.
Figure 2. (A) Principal component analysis (PCA) and (B) sparse partial least squares-discriminant analysis (sPLS-DA) of serum metabolites in grazing yaks during the cold season.
Metabolites 12 00738 g002
Metabolites 12 00738 g003 550
Figure 3. (A) summary of the differential serum metabolites in October or December compared with March; pathway impact analysis based on the differential serum metabolites between (B) October and December; (C) October and March; (D) December and March. Pathway impact > 0.20 and false discovery rate (FDR) < 0.05 were deemed as the key enriched pathways.
Figure 3. (A) summary of the differential serum metabolites in October or December compared with March; pathway impact analysis based on the differential serum metabolites between (B) October and December; (C) October and March; (D) December and March. Pathway impact > 0.20 and false discovery rate (FDR) < 0.05 were deemed as the key enriched pathways.
Metabolites 12 00738 g003
Metabolites 12 00738 g004 550
Figure 4. Random forest models of (A) serum indices with significant differences and (B) differential serum metabolites in grazing yaks during the cold season. Only the indices with a false discovery rate (FDR) < 0.05 are presented in the figures.
Figure 4. Random forest models of (A) serum indices with significant differences and (B) differential serum metabolites in grazing yaks during the cold season. Only the indices with a false discovery rate (FDR) < 0.05 are presented in the figures.
Metabolites 12 00738 g004
Table
Table 1. Nutritional value of forages in different months of cold season (dry matter basis).
Table 1. Nutritional value of forages in different months of cold season (dry matter basis).
Item (% of Dry Matter)OctDecMarp-Value
Organic matter97.75 ± 0.2497.73 ± 0.1997.93 ± 0.040.554
Non-fiber carbohydrate *25.59 ± 0.78 ab23.65 ± 1.34 b26.33 ± 0.73 a0.011
Crude protein4.48 ± 0.40 a3.80 ± 0.26 b4.17 ± 0.28 ab0.010
Ether extract2.24 ± 0.19 a2.42 ± 0.11 a1.32 ± 0.21 b0.002
Neutral detergent fiber65.44 ± 1.6667.86 ± 3.9666.11 ± 2.740.609
Acid detergent fiber36.63 ± 0.5837.18 ± 1.1337.01 ± 0.610.707
Calcium0.376 ± 0.047 a0.817 ± 0.022 b0.862 ± 0.038 b<0.001
Phosphorus0.026 ± 0.001 a0.017 ± 0.001 b0.024 ± 0.006 a0.035
Values with different superscripts represented significant differences among groups (p < 0.05). Results were presented as mean ± standard deviation. * Contents of non-fiber carbohydrates in forages were calculated by the contents of dry matter minus the percentages of crude protein, ash, ether extract, and neutral detergent fiber in forages.
Table
Table 2. Dynamic changes in serum biochemical indices of grazing yaks during cold season.
Table 2. Dynamic changes in serum biochemical indices of grazing yaks during cold season.
ItemsOctDecMarSEMp-Value
Protein metabolism
   Total Protein (g/L)78.4 a69.8 b64.8 b1.82<0.001
   Albumin (g/L)32.9 ab34.2 a30.5 b0.960.044
   Globulin (g/L)45.5 a35.6 b34.3 b1.860.001
   Albumin/Globulin0.735 b0.972 a0.895 ab0.0490.012
   Urea nitrogen (mmol/L)7.06 a6.04 a3.56 b0.5530.002
Energy metabolism
   Glucose (mmol/L)4.584.533.960.2230.130
   Total cholesterol (mmol/L)1.94 b2.53 a2.79 a0.111<0.001
   Triglyceride (mmol/L)0.2470.3500.3630.0550.291
Metabolic enzymes
   Alanine transaminase (μkat/L)0.4880.6530.5280.0530.104
   Aspartate aminotransferase (μkat/L)1.40 a1.07 b1.10 b0.0680.007
   Alkaline phosphatase (μkat/L)2.341.841.380.2560.057
   Lactate dehydrogenase (μkat/L)16.8 a14.4 ab12.4 b0.650.001
Minerals utilization
   Potassium (mmol/L)5.06 a5.24 a4.50 b0.098<0.001
   Sodium (mmol/L)142.3 a138.7 b140.0 ab0.910.038
   Calcium (mmol/L)2.482.412.260.0650.088
   Magnesium (mmol/L)0.987 a0.770 b0.745 b0.0470.005
   Iron (μmol/L)26.220.821.62.090.182
   Phosphorus (mmol/L)2.882.022.090.2270.053
Values with different superscripts represent significant differences among groups (p < 0.05).
Table
Table 3. Differential serum metabolites in yaks grazing in different months of the cold season.
Table 3. Differential serum metabolites in yaks grazing in different months of the cold season.
Metabolite NameRetention Time (min)MassSimilarityVIP aFold Change bFDR c
Oct vs. Dec
   Lyxose10.76217631.81.941.3580.019
   Threonic acid10.19292888.61.970.6520.019
   Threitol9.87217873.61.911.5270.027
   Tyrosine12.39218950.91.971.5230.027
   Leucrose16.8073506.81.912.1970.035
Oct vs. Mar
   Monostearin16.5757581.82.210.265<0.001
   1-Monopalmitin15.4057751.82.130.363<0.001
   Threonic acid10.19292888.62.110.522<0.001
   D-(glycerol 1-phosphate)11.39299793.42.060.3850.001
   Threitol9.87217873.62.021.5960.003
   Arachidonic acid14.4291737.31.980.3320.004
   Glycine7.23102819.31.860.2610.027
   Gluconic acid12.66292602.41.830.4060.033
   Lactic acid6.7759806.51.791.4730.048
Dec vs. Mar
   Monostearin16.5757581.82.370.278<0.001
   1-Monopalmitin15.4057751.82.300.487<0.001
   D-(glycerol 1-phosphate)11.39299793.42.340.339<0.001
   Arachidonic acid14.4291737.32.310.421<0.001
   α-Ketoglutaric acid10.33198845.42.080.5140.011
   Gluconic acid12.66292602.42.080.3460.012
   Serine8.93204934.12.100.5250.015
   Glycine7.23102819.32.020.4880.019
   Phosphate8.33158818.81.980.5910.023
   3-Hydroxybenzoic acid10.35120546.41.950.8430.033
   Glucuronic acid12.43333479.61.920.8650.034
   Threonic acid10.19292888.61.900.8000.047
a Variable important in projection (VIP) was obtained from the orthogonal partial least squares discriminant analysis (PLS-DA) of serum metabolites in yaks. b Fold changes were calculated by the area of each metabolite between groups (former month/later month). c FDR represents the false discovery rate.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gao, J.; Yang, D.; Sun, Z.; Niu, J.; Bao, Y.; Liu, S.; Tan, Z.; Hao, L.; Cheng, Y.; Liu, S. Changes in Blood Metabolic Profiles Reveal the Dietary Deficiencies of Specific Nutrients and Physiological Status of Grazing Yaks during the Cold Season in Qinghai Province of China. Metabolites 2022, 12, 738. https://doi.org/10.3390/metabo12080738

AMA Style

Gao J, Yang D, Sun Z, Niu J, Bao Y, Liu S, Tan Z, Hao L, Cheng Y, Liu S. Changes in Blood Metabolic Profiles Reveal the Dietary Deficiencies of Specific Nutrients and Physiological Status of Grazing Yaks during the Cold Season in Qinghai Province of China. Metabolites. 2022; 12(8):738. https://doi.org/10.3390/metabo12080738

Chicago/Turabian Style

Gao, Jian, Deyu Yang, Zhanying Sun, Jianzhang Niu, Yuhong Bao, Suozhu Liu, Zhankun Tan, Lizhuang Hao, Yanfen Cheng, and Shujie Liu. 2022. "Changes in Blood Metabolic Profiles Reveal the Dietary Deficiencies of Specific Nutrients and Physiological Status of Grazing Yaks during the Cold Season in Qinghai Province of China" Metabolites 12, no. 8: 738. https://doi.org/10.3390/metabo12080738

APA Style

Gao, J., Yang, D., Sun, Z., Niu, J., Bao, Y., Liu, S., Tan, Z., Hao, L., Cheng, Y., & Liu, S. (2022). Changes in Blood Metabolic Profiles Reveal the Dietary Deficiencies of Specific Nutrients and Physiological Status of Grazing Yaks during the Cold Season in Qinghai Province of China. Metabolites, 12(8), 738. https://doi.org/10.3390/metabo12080738

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