Subscriber access provided by UCL Library Services
Omics Technologies Applied to Agriculture and Food
Metabolomics Profiling of Serum and Urine in Three Beef Cattle Breeds Revealed Different Levels of Tolerance to Heat Stress Yupeng Liao, Rui Hu, Zhisheng Wang, Quanhui Peng, Xianwen Dong, Xiangfei Zhang, Huawei Zou, Qijian Pu, Bai Xue, and Lizhi Wang J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b01794 • Publication Date (Web): 15 Jun 2018 Downloaded from http://pubs.acs.org on June 16, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 41
Journal of Agricultural and Food Chemistry
Metabolomics Profiling of Serum and Urine in Three Beef Cattle Breeds Revealed Different Levels of Tolerance to Heat Stress Yupeng Liao, Rui Hu, Zhisheng Wang*, Quanhui Peng, Xianwen Dong, Xiangfei Zhang, Huawei Zou, Qijian Pu, Bai Xue, and Lizhi Wang Institute of Animal Nutrition, Sichuan Agricultural University, Key Laboratory of Bovine Low-Carbon Farming and Safe Production, Chengdu, Sichuan, 611130, China. *Corresponding author, Tel: +86-0835-2885730; Fax: +86-0835-2885730; E-mail:
[email protected].
1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
1
ABSTRACT: This study was to determine differences in the global metabolic profiles
2
of serum and urine of Xuanhan yellow cattle, Simmental crossbred cattle (Simmental
3
× Xuanhan yellow cattle) and cattle-yaks (Jersey × Maiwa yak) under heat stress
4
(temperature-humidity index remained above 80 for 1 week). A total of 55 different
5
metabolites associated with the three breeds were identified in the serum and urine
6
samples by gas chromatography-mass spectrometry. The metabolic adaptations to heat
7
stress are heterogeneous. Cattle-yaks mobilize a greater amount of body protein to
8
release glucogenic amino acids to supply energy, whereas the tricarboxylic acid cycle
9
is inhibited. Simmental crossbred cattle mobilize a greater amount of body fat to use
10
free fatty acids as an energy source. In comparison with Simmental crossbred cattle and
11
cattle-yaks, Xuanhan yellow cattle have higher glycolytic activity, and possess a
12
stronger antioxidant defense system and are, in conclusion, more adapted to hot and
13
humid environments.
14
KEYWORDS: metabolomics, gas chromatography-mass spectrometry, heat stress,
15
beef cattle, metabolic pathway
16 17 18 19 20 21 2
ACS Paragon Plus Environment
Page 2 of 41
Page 3 of 41
Journal of Agricultural and Food Chemistry
22
INTRODUCTION
23
Heat stress is a worldwide problem for the livestock industry. In particular, cattle are
24
extremely sensitive to hot and humid environments.1 In the US, heat stress leads to
25
annual economic losses in the beef cattle industry of 369 million dollars.2 In southern
26
China, summer heat stress causes beef cattle farms to be unprofitable.3 Furthermore, as
27
global temperatures rise,4 the beef cattle industry will face more serious challenges from
28
heat stress in the future. Earlier studies suggested that a reduction in dry matter intake
29
(DMI) by cattle exposed to heat stress is responsible for decreases in weight gain and
30
milk yield.5, 6 However, a series of recent studies have proven that heat stress reduces
31
DMI only partially and accounts for about 50% of the decrease in productivity, of which
32
the remainder is due to health problems and metabolic disorders in cattle.7, 8
33
Breed is a natural factor that can affect many aspects of livestock production, such
34
as productivity, temperature adaptability, reproductive function, temperament, and
35
immune function.9,
36
associated with breed.1,
37
account for approximately 80% of the national herd.12 Yellow cattle have a long history
38
of providing draft power, leather, milk, and meat for people in China. Xuanhan yellow
39
cattle, which are a Chinese southern local yellow cattle breed lived in the hilly areas at
40
low altitude, are well adapted to coarse feedstuffs and adverse conditions. However,
41
since the 1980s, market changes have caused an increase in interest in beef breeds with
42
other characteristics, in particular, meat yield performance. In this context, China has
10
It has been reported that heat tolerance in beef cattle is closely 11
In China, yellow cattle are the predominant breeds and
3
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 41
43
introduced a large number of foreign cattle for crossbreeding, such as Simmentals, to
44
improve the productivity of local yellow cattle. The cattle-yak was bred to combine the
45
tolerance of high altitudes and cold of yaks with the excellent productivity of other
46
cattle. Each beef cattle breed has its unique nutrition metabolism characteristics to be
47
adapted to distinct environments in the long-term evolution. In comparison with yellow
48
cattle, Yaks have more efficient energy and nitrogen utilization.13 At present, the above
49
three beef breeds are widely farmed in southern China, where is hot and humid
50
environments in summer. Nutrition metabolism alteration is an important way to
51
response to heat stress in mammals.14 However, under the same environment and diet,
52
the different levels of tolerance and nutrition metabolism to response to heat stress of
53
three beef breeds (Xuanhan yellow cattle, Simmental crossbreeds, and cattle-yaks) are
54
still unclear.
55
Metabolomics is an emerging technological and analytical field in systems biology,
56
following genomics, transcriptomics, and proteomics. The global metabolic profile of
57
an individual can directly reflect the final result of the interaction of a variety of factors,
58
including genetic, physiological, and environmental factors.15 Metabolomics employs
59
high-throughput
60
chromatography-mass spectrometry, and gas chromatography-mass spectrometry (GC-
61
MS), to study metabolites in biological samples (biofluids or tissues) to further explain
62
the physiological status of an individual.16 Among metabolomics methods, GC-MS has
63
been appreciated by researchers owing to its advantages of high resolution and
approaches,
such
as
nuclear
magnetic
4
ACS Paragon Plus Environment
resonance,
liquid
Page 5 of 41
Journal of Agricultural and Food Chemistry
64
detection sensitivity.17 Previous studies have demonstrated the practicability and
65
advantages of using metabolomics methods to study tolerance of seasonal weight loss
66
in sheep breeds.18 However, to date global metabolic profiling of cattle has focused on
67
dairy cows,19, 20 and information on metabolomics research into beef cattle and, more
68
importantly, the metabolic adaptation of the above three beef breeds under heat stress
69
has not been reported.
70
In this study, we aimed to determine differences in the global metabolic profiles of
71
serum and urine in the local breed (Xuanhan yellow cattle), Simmental crossbred cattle
72
(Simmental ×Xuanhan yellow cattle), and cattle-yaks (Jersey ×Maiwa yak) under heat
73
stress. We utilized a GC-MS-based metabolomics platform in combination with
74
multivariate statistical analysis to study different metabolic adaptations of the above
75
three beef cattle breeds under heat stress. The results will be of great importance for
76
understanding the heat stress tolerance of the above three beef breeds and may be
77
beneficial in the selection of beef cattle breeds for hot and humid areas.
78
MATERIALS AND METHODS
79
Animal Experiment and Sample Collection
80
The animal experiment in this study was approved by the Animal Care and Use
81
Committee of the Animal Nutrition Institute at Sichuan Agricultural University and was
82
carried out according to the Guide for the Care and Use of Laboratory Animals of the
83
National Research Council. In brief, 24 twenty-six-month-old steers were used,
84
including eight Xuanhan yellow cattle (XHC), eight Simmental × Xuanhan yellow 5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
85
cattle crossbred cattle (SXC), and eight Jersey ×Maiwa yak crossbred cattle (JMY). All
86
animals were fed with the same diet, which was designed according to the 2004 Chinese
87
feeding standard for beef cattle (NY/T 815–2004), and the concentrate-to-roughage
88
ratio of the diet was 6:4. The ingredients and chemical composition of the diet are
89
shown in Table S-1. The cattle were fed twice daily at 9:00 AM and 5:00 PM, and water
90
was provided at all times. According to the US National Research Council (1971), the
91
THI in the beef cattle house was calculated using the formula: THI = 0.72 ×(Td + Tw)
92
+ 40.6, where Td is the dry-bulb temperature (°C) and Tw is the wet-bulb temperature
93
(°C). Traditionally, the temperature-humidity index (THI) has been employed to
94
estimate the degree of heat stress in dairy cows, which are in heat stress when the THI
95
exceeds 72, and severe heat stress occurs if the THI exceeds 80.21 Serum and urine
96
samples were obtained in August (summer season), and the THI in the beef cattle house
97
was approximately 76–82 during this month (Fig. S-1). As for the sampling time
98
according to the previous literature,22 after the THI in the beef cattle house remained
99
above 80 for 1 week, fasting blood samples were obtained before morning feeding from
100
the jugular vein of the beef cattle. After the blood was left to stand for 30 min, the serum
101
was extracted from the blood samples by centrifugation at 3000 rpm for 10 min at 4 °C
102
and immediately stored in the dark at −80 °C until analysis was carried out. Urine
103
samples were collected in 15 mL Falcon tubes with a homemade urine bag at the same
104
time. Then, the urine samples were centrifuged at 4000 rpm for 10 min at 4 °C and
105
immediately stored in the dark at −80 °C until analysis was carried out. The rectal 6
ACS Paragon Plus Environment
Page 6 of 41
Page 7 of 41
Journal of Agricultural and Food Chemistry
106
temperature, respiration rate, DMI, and average daily gain (ADG) of the three beef
107
cattle breeds are shown in Table S-2.
108
Sample Pretreatment and GC-MS Analysis
109
Procedures for the extraction and derivatization of the serum and urine samples were
110
employed as previously described,23, 24 and were described in detail in the Supporting
111
Information (Supplemental Methods).
112
For the quality control samples, 20 μL of each prepared sample extract was taken
113
and mixed as previously described25 to monitor the reproducibility and stability of the
114
instrument.
115
The derivatized samples were analyzed by GC/MS (Agilent 7890A/5975C, Agilent
116
Technologies, Santa Clara, CA, USA). The detailed parameters were described under
117
Supplemental Methods.
118
Data Processing
119
The original GC/MS data were converted into files in netCDF format26 and
120
subsequently preprocessed, including the identification, filtration, and alignment of
121
peaks by XCMS software (www.bioconductor.org). Finally, a data matrix of
122
information, including mass-to-charge ratios (m/z), retention times, and peak intensities,
123
was obtained. Identification of metabolites was performed using the AMIDS system
124
(Automated Mass Spectral Deconvolution and Identification System) by searching the
125
National Institute of Standards and Technology library (http://srdata.nist.gov/gateway/)
126
and the Wiley Chemical Structure Library.27 7
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
127
The normalized GC-MS datasets were imported into the SIMCA-P 13.0 software
128
package (Umetrics, Umea, Sweden) for multivariate analysis, including principal
129
component analysis (PCA) and partial least squares discriminant analysis (PLS-DA).
130
Prior to the multivariate analysis, mean centering and unit variance scaling were
131
performed on the GC-MS data to make the metabolites equal in importance. PCA was
132
utilized to visualize global clustering and display differences in metabolic profiles
133
between samples. To achieve the maximum separation between samples, PLS-DA was
134
used to identify differential metabolites that explained the separation between groups.
135
The model quality was evaluated using the R2X, R2Y, and Q2 parameters. The R2X and
136
R2Y parameters, which represent the fractions of explained X-variation and Y-variation,
137
respectively, can be used to evaluate the quality of the model. The Q2 parameter
138
represents the predictive ability of the model. In general, the model is acceptable when
139
the values of R2X, R2Y, and Q2 are greater than 0.5.24 In addition to cross-validation,
140
the PLS-DA models were verified via 100 repeated permutation tests.
141
Identification of Differential Metabolites and Metabolic Pathway
142
Analysis
143
As described above, PLS-DA was employed to identify differential metabolites
144
between groups. Furthermore, variable importance in the projection (VIP) values in the
145
PLS-DA model were used to rank the metabolites based on their importance in
146
discriminating between groups. Metabolites with the highest VIP values are the most
147
powerful group discriminators. Traditionally, VIP values >1 are significant. In this 8
ACS Paragon Plus Environment
Page 8 of 41
Page 9 of 41
Journal of Agricultural and Food Chemistry
148
study, metabolites with VIP values(VIP > 1.0)in the PLS-DA model and significance
149
(p < 0.05) in Student’s t-test (SPSS 16.0) were selected as differential metabolites. Then,
150
receiver operating characteristic (ROC) curves analysis was performed (SPSS 16.0),
151
and the areas under the curves (AUCs) were calculated to determine the diagnostic
152
value of these differential metabolites.28 The discriminatory power of the differential
153
metabolites was ranked and visualized using heat maps.28 Metabolic pathways were
154
obtained for these identified differential metabolites from the KEGG (Kyoto
155
Encyclopedia of Genes and Genomes) database (http://www.genome.jp/kegg).
156
Determination of Metabolism and Activity of Antioxidant Enzymes
157
The levels of lactate dehydrogenase (LDH), malondialdehyde (MDA), superoxide
158
dismutase (SOD), and glutathione peroxidase (GSH-Px) and the total antioxidant
159
capacity (T-AOC) in serum samples were measured using enzyme-linked
160
immunosorbent assay kits (Chengdu Anya Biotechnology Co., Ltd). Absorbance values
161
were determined at 450 nm with a microplate reader. The calibration standards are
162
assayed at the same time as the samples and allow the operator to produce a standard
163
curve of Optical Density. The levels of these enzymes in the samples is then determined
164
by comparing the O.D. of the samples to the standard curve.
165
Enzyme activity was analyzed by univariate analysis of variance with a post hoc
166
Tukey’s test using SPSS 16.0 software (Chicago, IL, USA). Statistical significance was
167
set at p < 0.05, and a high level of significance at p < 0.01.
168
RESULTS 9
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
169
Metabolomics Profiling of Serum and Urine Samples
170
Typical GC/MS chromatograms of the serum and urine samples from the three beef
171
cattle breeds are shown in Figures 1A and B, respectively. In total, 167 and 220 valid
172
peaks were identified in the chromatograms of serum and urine, respectively. A total of
173
78 and 81 metabolites were identified and quantified in the serum and urine samples,
174
respectively. These metabolites in serum and/or urine can be divided into eight major
175
groups, namely, carbohydrates, amines, amino acids, fatty acids, organic acids,
176
phosphoric acid, polyol, and nucleotides, on the basis of the characteristics of each
177
chemical. Furthermore, the five metabolites with the highest concentrations in serum
178
were glucose, urea, lactic acid, phosphoric acid, and cholesterol in all breeds. In contrast,
179
the five metabolites with the highest concentrations in urine were slightly different in
180
the three breeds: in the XHC group these were phosphoric acid, creatinine,
181
pseudouridine, uric acid, and benzoic acid, whereas in the SXC and JMY groups these
182
were phosphoric acid, creatinine, pseudouridine, allantoin, and benzoic acid.
183
Data for the serum and urine samples acquired by GC-MS were integrated and
184
analyzed using PCA (Fig. S-2). The PCA score plots for the serum samples showed a
185
distinct separation between the three breeds. For the urine samples, the PCA model was
186
not able to completely separate the three breeds, and we therefore focused mainly on
187
the results of supervised analysis. The PLS-DA score plots for the two sample types
188
(serum and urine) showed a clear separation between the XHC, SXC, and JMY groups,
189
which suggested that the serum and urine metabolic profiles of the three breeds were 10
ACS Paragon Plus Environment
Page 10 of 41
Page 11 of 41
Journal of Agricultural and Food Chemistry
190
significantly different under heat stress (Fig. 2A and E). Furthermore, to further
191
investigate the metabolic status and discover potential biomarkers of each breed, PCA
192
and PLS-DA analysis were also applied to each combination of two breeds (XHC
193
versus SXC, XHC versus JMY, and SXC versus JMY), all of which showed distinct
194
separations (Fig. S-2, Fig. 2B–D and F–H).
195
Identification of Different Metabolites
196
In the present study, different metabolites were selected according to their VIP values
197
in the PLS-DA model (VIP > 1) and p-values from Student’s t-test (p < 0.05). A total
198
of 42 and 23 different metabolites associated with serum and urine from the breeds are
199
listed in Tables 1 and 2, respectively. Specifically, 25, 30, and 23 metabolites from
200
serum were identified in the XHC group in comparison with the SXC group, in the XHC
201
group in comparison with the JMY group, and in the SXC group in comparison with
202
the JMY group, respectively. Moreover, 11, 16, and 12 metabolites from urine were
203
identified in the XHC group in comparison with the SXC group, in the XHC group in
204
comparison with the JMY group, and in the SXC group in comparison with the JMY
205
group, respectively.
206
In order to determine the diagnostic value of these different metabolites for
207
differentiating any two of the XHC, SXC, and JMY groups, ROC analysis was carried
208
out. A heat map shows any two groups of different metabolites ranked in order of their
209
AUC values (Fig. S-3). All the different metabolites displayed good diagnostic abilities,
210
with AUC values of 0.73–1.00. 11
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
211
Integration of Key Different Metabolic Pathways
212
The KEGG pathway database was utilized for analyzing related metabolic pathways
213
of 55 metabolites found in serum and urine. As summarized in Tables 1 and 2, we found
214
that these metabolites were involved in multiple biochemical pathways, such as
215
glycolysis, amino acid metabolism, fatty acid metabolism, inositol phosphate
216
metabolism, vitamin E metabolism, purine metabolism, and the tricarboxylic acid
217
(TCA) cycle. Finally, we combined these results to draw a metabolic network map (Fig.
218
3) to show a more intuitive correlation between these metabolites.
219
Different Activities of Metabolic and Antioxidant Enzymes
220
Figure 4 shows the activities of key enzymes related to glycolysis and the antioxidant
221
defense system. LDH activity was very significantly higher in the XHC group in
222
comparison with the SXC and JMY groups (Fig. 4A). Furthermore, MDA levels were
223
significantly lower in the XHC group in comparison with the SXC and JMY groups
224
(Fig. 4B), and SOD activity and T-AOC were significantly higher in the XHC group in
225
comparison with the SXC and JMY groups (Fig. 4B).
226
DISCUSSION
227
We found heat stress induced significantly reduction of dry matter intake per kg
228
BW (-18.27%, -24.75%, and -28.77%, respectively) and average daily gain (-17.33%,
229
-21.43%, and -24.68%, respectively) of Xuanhan yellow cattle, Simmental crossbred
230
cattle and cattle-yaks when compared with our previous data during non-heat stress
231
period. Accordingly, heat stress may had lowest impact on production performance of 12
ACS Paragon Plus Environment
Page 12 of 41
Page 13 of 41
Journal of Agricultural and Food Chemistry
232
Xuanhan yellow cattle among three breeds and cattle-yak suffered the greatest
233
negative interference. So we wanted to understand differences in the global metabolic
234
profiles of the three breeds under heat stress. In this study, we utilized GC-MS-based
235
metabolomics to determine differences in the metabolite profiles of serum and urine
236
in Xuanhan yellow cattle, Simmental crossbred cattle, and cattle-yaks under heat
237
stress. A total of 78 and 81 metabolites were identified and quantified in serum and
238
urine, respectively. In total, 55 different metabolites associated with the three breeds
239
were obtained under heat stress. These metabolites were involved in glycolysis, amino
240
acid metabolism, fatty acid metabolism, the TCA cycle, and purine metabolism. Our
241
results revealed the bases of different metabolic adaptations in the three breeds under
242
heat stress and could benefit breeding programs.
243
Glycolysis
244
A previous study showed that heat stress enhanced glycolysis, which led to
245
increases in the levels of pyruvic acid and lactic acid in dairy cows.22 In the present
246
study, we found that serum levels of glucose, pyruvic acid, and lactic acid were higher
247
in the XHC group in comparison with the SXC and JMY groups. Moreover, LDH
248
activity was significantly higher in the XHC group in comparison with the SXC and
249
JMY groups. All the above results indicate that glycolytic activity and anaerobic cell
250
respiration rates were higher in the XHC group in comparison with the SXC and JMY
251
groups. Pyruvic acid, which is formed by the breakdown of glucose, is an important
252
pathway junction in carbohydrate catabolism and can be further converted into acetyl13
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
253
CoA to supply energy via the TCA cycle under aerobic conditions.29 However, under
254
anaerobic conditions pyruvic acid is reduced to lactic acid by LDH. Lactic acid is thus
255
the end-product of glucose metabolism under anaerobic conditions,30 which prevent
256
the inhibition of glycolysis, but lactic acid can be oxidized to pyruvic acid in the
257
presence of oxygen. Therefore, the higher glycolytic activity in the XHC group could
258
be a metabolic adaptation to supply more energy under heat stress.
259
Amino Acid Metabolism
260
Heat stress increases the disintegration of proteins and mobilization of amino acids
261
to enable the production of additional energy. These processes produce higher levels
262
of glucogenic amino acids, creatinine, and urea.31, 32 In the present study, we found
263
that serum levels of glucogenic amino acids, such as methionine, glutamine,
264
phenylalanine, tyrosine, and asparagine, which could generate pyruvic acid, α-
265
ketoglutaric acid, fumaric acid, and oxaloacetic acid,33 were significantly higher in the
266
JMY group in comparison with the SXC and XHC groups. Moreover, levels of
267
ornithine and urea, which participate in the urea cycle,34 were also higher in the JMY
268
group. On the other hand, levels of glutamic acid, glutamine, asparagine, ornithine,
269
and urea were significantly higher in the SXC group in comparison with the XHC
270
group. In addition, serum creatinine levels were significantly higher in the JMY and
271
SXC groups in comparison with the XHC group. Creatinine, which is the degradation
272
product of creatine phosphate, plays an important role in energy balance when
273
animals have high energy requirements. The higher levels of creatinine in the JMY 14
ACS Paragon Plus Environment
Page 14 of 41
Page 15 of 41
Journal of Agricultural and Food Chemistry
274
and SXC groups could indicate an increase in the mobilization of creatine phosphate
275
in muscle tissue for energy supply.35 All the above results indicate that the catabolism
276
of proteins and mobilization of amino acids were highest in the JMY group, followed
277
by the SXC group, and lowest in the XHC group. Therefore, the mobilization of
278
muscle protein in the JMY group could be a metabolic adaptation to release
279
glucogenic amino acids that the body can subsequently convert into energy-rich
280
metabolites under heat stress.
281
Interestingly, serum levels of putrescine, which is a polyamine derived from
282
arginine, were significantly lower in the JMY group in comparison with the XHC and
283
SXC groups. Putrescine can improve the response of heat shock proteins36 and change
284
the topology of DNA to an extent that promotes survival.37 In combination with the
285
higher level of ornithine, we think that the JMY breed may require more putrescine to
286
resist heat stress. However, lower levels of putrescine may interfere with heat
287
adaptability in the JMY breed.
288
Lipid Metabolism
289
Heat stress leads to an increase in circulating free fatty acids, including linoleic
290
acid, oleic acid, and arachidonic acid, because of lipid catabolism.38 In the present
291
study, we found that serum and urine levels of glycerol and free fatty acids, including
292
palmitic acid, tetradecanoic acid, oleic acid, linoleic acid, and arachidonic acid, which
293
belong to the family of long-chain fatty acids, were significantly higher in the SXC
294
group in comparison with the XHC and JMY groups. Moreover, levels of glycerol and 15
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
295
free fatty acids were significantly higher in the JMY group in comparison with the
296
XHC group. In addition, urine levels of myo-inositol and nicotinic acid, which can
297
both play an important role in fat metabolism,39 were higher in the SXC and JMY
298
groups in comparison with the XHC group. All the above results suggest that fat
299
catabolism was highest in the SXC group, followed by the JMY group, and lowest in
300
the XHC group. In the liver, free fatty acids can be converted into acetyl-CoA by β-
301
oxidation to supply energy via the TCA cycle. Glycerol can be converted into
302
glycerol-3-phosphate and glyceric acid via glycolysis to supply energy.40 Therefore,
303
the mobilization of body fat in the SXC group could be a metabolic adaptation to use
304
free fatty acids and glycerol as energy sources under heat stress.
305
However, it is noteworthy that previous research has confirmed that palmitic acid is
306
associated with the activation of inflammatory responses in animal models,41 and
307
omega-6 linoleic acid and arachidonic acid also promote inflammation.42 Therefore,
308
the SXC and JMY breeds may have a greater risk of inflammation in comparison with
309
the XHC breed.
310
TCA Cycle
311
In previous studies, heat stress led to upregulation of the TCA cycle, which
312
increased levels of TCA cycle intermediates, such as citric acid, cis-aconitic acid, α-
313
ketoglutaric acid, succinic acid, fumaric acid, and malic acid.38, 43 In the present study,
314
we found that serum and urine levels of citric acid, aconitic acid, and fumaric acid
315
were significantly higher in the SXC group in comparison with the JMY and XHC 16
ACS Paragon Plus Environment
Page 16 of 41
Page 17 of 41
Journal of Agricultural and Food Chemistry
316
groups. Moreover, levels of citric acid, aconitic acid, and succinic acid were
317
significantly higher in the XHC group in comparison with the JMY group. These
318
results suggest that TCA cycle activity was highest in the SXC group, followed by the
319
XHC group, and lowest in the JMY group. This may result from the higher level of
320
acetyl-CoA due to the β-oxidation of fatty acids in the SXC group.
321
Interestingly, the JMY group had higher urine levels of methylmalonic acid and
322
methylcitric acid in comparison with the SXC and XHC groups. Medically, the
323
methylmalonic acid concentration is used to diagnose methylmalonic aciduria, which
324
is caused by an inherited deficiency of the mitochondrial enzyme methylmalonyl-CoA
325
mutase or a metabolic disorder involving its coenzyme cobalamin.44 A deficiency of
326
methylmalonyl-CoA mutase, which catalyzes the isomerization of methylmalonyl-
327
CoA to succinyl-CoA, results in an accumulation of methylmalonic acid and
328
methylcitric acid and thus inhibits the TCA cycle.45 It was first discovered that the
329
JMY breed has a greater risk of methylmalonic aciduria in comparison with the SXC
330
and XHC breeds. In brief, under heat stress the TCA cycle activity in the JMY breed
331
is the lowest in the three breeds, the reason for which needs further research.
332
Oxidative Stress
333
Heat stress increases the production of reactive oxygen species (ROS), including
334
superoxide ions (O2-), nitrogen monoxide (NO), and hydroxyl radicals (HO·), which
335
destroy the structure and function of the cell membrane and damage the mitochondria,
336
eventually leading to cell death.46 ROS are balanced by natural antioxidant 17
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
337
compounds such as glutathione, GSH-Px, SOD, vitamin C, and vitamin E.47 Previous
338
studies revealed that heat stress increased oxidative stress, which increased MDA
339
levels48 and decreased levels of SOD and GSH-Px in cattle exposed to heat stress.49
340
We found that serum MDA levels were significantly lower in the XHC group in
341
comparison with the SXC and JMY groups, and SOD activity and T-AOC were
342
significantly higher in the XHC group in comparison with the SXC and JMY groups.
343
Furthermore, in the present study we found that urine levels of uric acid were
344
significantly higher and those of allantoin were lower in the XHC group in
345
comparison with the SXC and JMY groups. As we know that cattle can express the
346
enzyme uricase, which converts uric acid into allantoin. Allantoin is thus the final
347
oxidation product of purine metabolism. In the presence of oxidative stress, uric acid
348
is oxidized by ROS to allantoin.50 Previous studies have shown that allantoin is a
349
potential marker for monitoring oxidative status in humans.51 The higher levels of uric
350
acid and lower levels of allantoin in the XHC group may indicate that this breed
351
suffers from less oxidative stress under heat stress. α-Tocopherol and γ-tocopherol are
352
forms of vitamin E and possess potent anti-inflammatory and antioxidant properties.52
353
In the present study, we found that serum levels of α-tocopherol and γ-tocopherol
354
were significantly higher in the XHC group in comparison with the SXC group.
355
Moreover, serum levels of γ-tocopherol were significantly higher in the XHC group in
356
comparison with the JMY group. The lower levels of tocopherols in the SXC and
357
JMY groups may result from more severe lipid inflammation and greater oxidative 18
ACS Paragon Plus Environment
Page 18 of 41
Page 19 of 41
Journal of Agricultural and Food Chemistry
358
stress under heat stress. All the above results suggest that the XHC breed suffers from
359
less oxidative stress and possesses a stronger antioxidant defense system in
360
comparison with the SXC and JMY breeds under heat stress.
361
In summary, under heat stress cattle-yaks (Jersey × Maiwa yak) mobilize a greater
362
amount of body protein to release glucogenic amino acids to supply energy, which
363
may represent their metabolic adaptation mechanism, whereas the TCA cycle is
364
inhibited; the above results indicate that the heat stress tolerance of cattle-yaks is the
365
lowest among the three breeds. In contrast, Simmental crossbred cattle (Simmental ×
366
Xuanhan yellow cattle) mobilize a greater amount of body fat in order to use free fatty
367
acids and glycerol as energy sources, and the TCA cycle activity is higher, which may
368
represent their metabolic adaptation mechanism. In comparison with Simmental
369
crossbred cattle and cattle-yaks, Xuanhan yellow cattle change their nutritional energy
370
source without significantly compromising the storage of body fat and protein and
371
have higher levels of pyruvic acid and lactic acid, which indicates that the glycolysis
372
process is enhanced. Furthermore, Xuanhan yellow cattle suffer from less oxidative
373
stress and possess a stronger antioxidant defense system under heat stress and are thus
374
more adapted to hot and humid environments.
375
ASSOCIATED CONTENT
376
Supporting Information
377
Figure S-1, THI in August; Figure S-2, PCA score plots for the serum and urine
378
samples; Figure S-3, heat maps showing the AUC values of different metabolites; 19
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
379
Table S-1, ingredients and chemical composition of the diet; Table S-2, physiological
380
index and production performance; supplemental methods.
381
AUTHOR INFORMATION
382
Corresponding Author
383
*Phone/fax: +86-0835-2885730. E-mail:
[email protected].
384
Author Contributions
385
Y.L. and R.H. contributed equally and should be considered co-first authors.
386
Notes
387
The authors declare no competing financial interest.
388
ACKNOWLEDGMENTS
389
The authors gratefully thank the members of the Institute of Animal Nutrition,
390
Sichuan Agricultural University for their assistance in the sampling and analysis of
391
the samples. We also thank Suzhou BioNovoGene Co., Ltd for their assistance in the
392
original data processing and related bioinformatics analysis of GC-MS.
393
FUNDING SOURCES
394
This work was funded by National key research and development program of China
395
(2017YFD0502005), Sichuan Science and Technology Program (2018NZ0002) and
396
China Agriculture (Beef Cattle/Yak) Research System (CARS-37)
397
(http://www.beefsys.com/).
398 399 20
ACS Paragon Plus Environment
Page 20 of 41
Page 21 of 41
400 401 402 403 404
Journal of Agricultural and Food Chemistry
REFERENCES (1) Silanikove, N. Effects of heat stress on the welfare of extensively managed domestic ruminants. Livestock Production Science. 2000, 67 (1), 1-18. (2) St-Pierre, N. R.; Cobanov, B.; Schnitkey, G. Economic losses from heat stress by US livestock industries. Journal of dairy science. 2003, 86 (5), E52–77.
405
(3) Song, X.; Luo, J.; Fu, D. Traditional Chinese medicine prescriptions enhance
406
growth performance of heat stressed beef cattle by relieving heat stress responses and
407
increasing apparent nutrient digestibility. Asian-Australia journal Animal Science.
408
2014, 27 (10), 1513.
409 410 411 412
(4) Luber, G.; McGeehin, M. Climate change and extreme heat events. American journal of preventive medicine. 2008, 35 (5), 429–435. (5) Fuquay, J. W. Heat stress as it affects animal production. Journal of Animal Science. 1981, 52 (1), 164–174.
413
(6) West, J. W.; Mullinix, B. G.; Bernard, J. K. Effects of hot, humid weather on milk
414
temperature, dry matter intake, and milk yield of lactating dairy cows. Journal of dairy
415
science. 2003, 86 (1), 232-242.
416
(7) Bernabucci, U.; Biffani, S.; Buggiotti, L.; Vitali, A.; Lacetera, N.; Nardone, A.
417
The effects of heat stress in Italian Holstein dairy cattle. Journal of dairy science. 2014,
418
97 (1), 471–486.
419
(8) Wheelock, J. B.; Rhoads, R. P.; Vanbaale, M. J.; Sanders, S. R.; Baumgard, L. H.
420
Effects of heat stress on energetic metabolism in lactating Holstein cows. Journal of 21
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
421
dairy science. 2010, 93 (2), 644–655.
422
(9) Gregory, K. E.; Cundiff, L.V.; Koch, R. M. Breed effects and heterosis in
423
advance generations of composite populations for reproduction and maternal traits of
424
beef cattle. Journal of Animal Science. 1991, 69 (7), 2795.
425
(10) Carroll, J. A.; Burdick, N. C.; Reuter, R. R.; Chase, J. C.; Spiers, D. E.;
426
Arthington, J. D.; Coleman, S. W. Differential acute phase immune responses by Angus
427
and Romosinuano steers following an endotoxin challenge. Domestic Animal
428
Endocrinology. 2011, 41 (4), 163–73.
429 430 431 432
(11) Gaughan, J. B.; Mader, T. L.; Holt, S. M. Heat tolerance of Boran and Tuli crossbred steers. Journal of Animal Science. 1999, 77 (9), 2398-405. (12) Zhou, G. H.; Liu, L.; Xiu, X. L. Productivity and carcass characteristics of pure and crossbred Chinese Yellow Cattle. Meat Science. 2001, 58 (4), 359-362.
433
(13) Qiu, Q.; Zhang, G.; Ma, T.; Qian, W.; Wang, J.; Ye, Z.; Cao, C.; Hu, Q.; Kim, J.;
434
Larkin, D. M.; Auvil, L.; Capitanu, B.; Ma, J.; Lewin, H. A.; Qian, X.; Lang, Y.; Zhou,
435
R.; Wang, L.; Wang, K.; Xia, J.; Liao, S.; Pan, S.; Lu, X.; Hou, H.; Wang, Y.; Zang, X.;
436
Yin, Y.; Ma, H.; Zhang, J.; Wang, Z.; Zhang, Y.; Zhang, D.; Yonezawa, T.; Hasegawa,
437
M.; Zhong, Y.; Liu, W.; Zhang, Y.; Huang, Z.; Zhang, S.; Long, R.; Yang, H.; Wang, J.;
438
Lenstra, J. A.; Cooper, D. N.; Wu, Y.; Wang, J.; Shi, P.; Wang, J.; Liu, J. The yak
439
genome and adaptation to life at high altitude. Nature genetics. 2012, 44, (8), 946-949.
440
(14) Ippolito, D. L.; Lewis, J. A.; Yu, C.; Leon, L. R.; Stallings, J. D. Alteration in
441
circulating metabolites during and after heat stress in the conscious rat: potential 22
ACS Paragon Plus Environment
Page 22 of 41
Page 23 of 41
442 443 444
Journal of Agricultural and Food Chemistry
biomarkers of exposure and organ-specific injury. BMC physiology. 2014, 14, 14. (15) Nicholson, J. K. Global systems biology, personalized medicine and molecular epidemiology. Molecular. Systems. Biology. 2006, 2 (1), 52.
445
(16) Zhao, S.; Zhao, J.; Bu, D.; Sun, P.; Wang, J.; Dong, Z. Metabolomics analysis
446
reveals large effect of roughage types on rumen microbial metabolic profile in dairy
447
cows. Microbiol. 2014, 59 (1), 79−85.
448
(17) Lindon, J. C.; Holmes, E.; Nicholson, J. K. Metabonomics techniques and
449
applications to pharmaceutical research & development. Pharmaceutical Research.
450
2006, 23 (6), 1075−1088.
451
(18) Palma, M.; Scanlon, T.; Kilminster, T.; Milton, J.; Oldham, C.; Greeff, J.;
452
Matzapetakis, M.; Almeida, A. M., The hepatic and skeletal muscle ovine metabolomes
453
as affected by weight loss: a study in three sheep breeds using NMR-metabolomics.
454
Scientific reports. 2016, 6, 39120.
455
(19) Sun, H. Z.; Wang, D. M.; Wang, B.; Wang, J. K.; Liu, H. Y.; Guan le, L.; Liu,
456
J. X., Metabolomics of four biofluids from dairy cows: potential biomarkers for milk
457
production and quality. Journal of proteome research. 2015, 14 (2), 1287-98.
458
(20) Dervishi, E.; Zhang, G.; Dunn, S. M.; Mandal, R.; Wishart, D. S.; Ametaj, B.
459
N., GC-MS Metabolomics Identifies Metabolite Alterations That Precede Subclinical
460
Mastitis in the Blood of Transition Dairy Cows. Journal of proteome research. 2017,
461
16 (2), 433-446.
462
(21) Bohmanova, J.; Misztal, I.; Cole, J. B. Temperature-Humidity Indices as 23
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
463
Indicators of Milk Production Losses due to Heat Stress. Journal of dairy science. 2007,
464
90 (4), 1947-1956.
465
(22) Tian, H.; Wang, W.; Zheng, N.; Cheng, J.; Li, S.; Zhang, Y.; Wang, J.,
466
Identification of diagnostic biomarkers and metabolic pathway shifts of heat-stressed
467
lactating dairy cows. Journal of proteomics. 2015, 125, 17-28.
468
(23) Dunn, Warwick B. Procedures for large-scale metabolic profiling of serum and
469
plasma using gas chromatography and liquid chromatography coupled to mass
470
spectrometry. Nature Protocols. 2011, 6 (7), 1060-1083.
471
(24) Chan, E. C.; Pasikanti, K. K.; Nicholson, J. K. Global urinary metabolic
472
profiling procedures using gas chromatography-mass spectrometry. Nature Protocols.
473
2011, 6 (10), 1483–1499.
474
(25) Sangster.; Timothy. A pragmatic and readily implemented quality control
475
strategy for HPLC-MS and GC-MS-based metabonomic analysis. Analyst. 2006, 131
476
(10), 1075-1078.
477
(26) Smith, C. A.; Want, E. J.; O’Maille, G.; Abagyan, R.; Siuzdak, G. XCMS:
478
processing mass spectrometry data for metabolite profiling using nonlinear peak
479
alignment, matching, and identification. Analytical Chemistry. 2006, 78 (3), 779–787.
480
(27) Oberacher, H; Whitley, G; Berger B. Evaluation of the sensitivity of the 'Wiley
481
registry of tandem mass spectral data, MSforID' with MS/MS data of the
482
'NIST/NIH/EPA mass spectral library'. Journal of Mass Spectrometry. 2013, 48 (4),
483
487-496. 24
ACS Paragon Plus Environment
Page 24 of 41
Page 25 of 41
Journal of Agricultural and Food Chemistry
484
(28) Dunn, W. B.; Broadhurst, D.; Begley, P.; Zelena, E.; Francis-McIntyre, S.;
485
Anderson, N. Human Serum Metabolome (HUSERMET) Consortium. Procedures for
486
large-scale metabolic profiling of serum and plasma using gas chromatography and
487
liquid chromatography coupled to mass spectrometry. Nature Protocols. 2011, 6 (7),
488
1060–1083.
489
(29) Birungi, G.; Chen, S. M.; Loy, B. P.; Ng, M. L.; Li, S. F. Metabolomics approach
490
for investigation of effects of dengue virus infection using the EA.hy926 cell line.
491
Journal of Proteome Research. 2010, 9 (12), 6523-6534.
492
(30) Xuan, J.; Pan, G.; Qiu, Y.; Yang, L.; Su, M.; Liu, Y.; Chen, J.; Feng, G.; Fang,
493
Y.; Jia, W.; Xing, Q.; He, L. Metabolomic Profiling to Identify Potential Serum
494
Biomarkers for Schizophrenia and Risperidone Action. Journal of proteome research.
495
2011, 10 (12), 5433-5443.
496
(31) Hou, Y.; Wang, X.; Lei, Z.; Ping, J.; Liu, J.; Ma, Z.; Zhang, Z.; Jia, C.; Jin, M.;
497
Li, X.; Li, X.; Chen, S.; Lv, Y.; Gao, Y.; Jia, W.; Su, J., Heat-stress-induced metabolic
498
changes and altered male reproductive function. Journal of proteome research. 2015,
499
14 (3), 1495-1503.
500 501
(32) Gao, S. T.; Guo, J.; Quan, S. Y. The effects of heat stress on protein metabolism in lactating Holstein cows. Journal of dairy science. 2017, 100, 1-10.
502
(33) Lapierre, H.; Lobley, G. E.; Doepel, L. Triennial Lactation Symposium:
503
Mammary metabolism of amino acids in dairy cows. Journal of Animal Science. 2012,
504
90 (5), 1708-1721. 25
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
505
(34) Motyl, T.; Barej, W. Plasma amino acid indices and urinary 3-methyl histidine
506
excretion in dairy cows in early lactation. Annales De Recherches Veterinaires Annals
507
of Veterinary Research. 1986, 17 (2), 153–157.
508
(35) Serrano-Contreras, J. I.; Garcfa-Perez, I.; Melendez-Camargo, M. E.; Zepeda,
509
L. G. NMR-based metabonomic analysis of physiological responses to starvation and
510
refeeding in the rat. Journal of proteome research. 2016, 15, 3241–3254.
511
(36) Robert Michaud, M.; Benoit, J. B.; Lopez-Martinez, G.; Elnitsky, M. A.; Lee,
512
R. E., Jr.; Denlinger, D. L., Metabolomics reveals unique and shared metabolic changes
513
in response to heat shock, freezing and desiccation in the Antarctic midge, Belgica
514
antarctica. Journal of insect physiology. 2008, 54 (4), 645-55.
515
(37) Tkachenko, A.G.; Pshenichnov, M. R.; Salakhetdnova, O. Y.; Nesterova, L. Y.
516
The role of putrescine and potassium transport in the regulation of DNA topology
517
during Escherichia coli adaptation to heat stress. Mikrobiologiya. 1998, 67 (5), 601–
518
606.
519
(38) Hou, Y.; Wang, X.; Ping, J.; Lei, Z.; Gao, Y.; Ma, Z.; Jia, C.; Zhang, Z.; Li, X.;
520
Jin, M.; Li, X.; Suo, C.; Zhang, Y.; Su, J., Metabonomics Approach to Assessing the
521
Modulatory Effects of Kisspeptin-10 on Liver Injury Induced by Heat Stress in Rats.
522
Scientific reports. 2017, 7 (1), 7020.
523
(39) Sun, L. W.; Zhang, H. Y.; Wu, L.; Shu, S.; Xia, C.; Xu, C.; Zheng, J. S., (1)H-
524
Nuclear magnetic resonance-based plasma metabolic profiling of dairy cows with
525
clinical and subclinical ketosis. Journal of dairy science. 2014, 97 (3), 1552-1562. 26
ACS Paragon Plus Environment
Page 26 of 41
Page 27 of 41
Journal of Agricultural and Food Chemistry
526
(40) Zhou, J.; Chen, J.; Hu, C.; Xie, Z.; Li, H.; Wei, S.; Wang, D.; Wen, C.; Xu, G.,
527
Exploration of the serum metabolite signature in patients with rheumatoid arthritis
528
using gas chromatography-mass spectrometry. Journal of pharmaceutical and
529
biomedical analysis. 2016, 127, 60-67.
530
(41) Sieber, J.; Lindenmeyer, M. T.; Kampe, K.; Campbell, K. N.; Cohen, C. D.;
531
Hopfer, H.; Mundel, P.; Jehle, A.W. Regulation of podocyte survival and
532
endoplasmicreticulum stress by fatty acids. American journal of physiology-renal
533
physiology. 2010, 299 (40, 821–829.
534
(42) Lazic, M.; Inzaugarat, M. E.; Povero, D. Reduced Dietary Omega-6 to Omega-
535
3 Fatty Acid Ratio and 12/15-Lipoxygenase Deficiency Are Protective against Chronic
536
High Fat Diet-Induced Steatohepatitis. Plos One. 2014, 9 (9), e107658.
537
(43) Shi, Y.; Wang, D., Implication of metabolomic profiles to wide thermoneutral
538
zone in Mongolian gerbils (Meriones unguiculatus). Integrative zoology. 2016, 11 (4),
539
282-94.
540
(44) Yang, X.; Sakamoto, O.; Matsubara, Y. Mutation analysis of the MMAA and
541
MMAB genes in Japanese patients with vitamin B12-responsive methylmalonic
542
acidemia:identification of a prevalent MMAA mutation. Molecular genetics and
543
metabolism. 2004, 82 (4), 329-333.
544
(45) Okun, J. G.; Horster, F.; Farkas, L. M. Neurodegeneration in methylmalonic
545
aciduria involves inhibition of complex Ⅱ and the tricarboxylic acid cycle, and
546
synergistically acting excitotoxicity. Journal of Biological Chemistry. 2002, 277 (17), 27
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
547 548 549 550 551
14674-14680. (46) Ray, P. D.; Huang, B. W.; Tsuji, Y. Reactive oxygen species (ROS) homeostasis and redox regulation in cellular signaling. Cell Signal. 2012, 24 (5), 981-990. (47) Miller, J. K.; Brzezinska-Slebodzinska, E.; Madsen, F. C. Oxidative stress, antioxidant and animals function. Journal of dairy science. 1993, 76 (9), 2812-2823.
552
(48) Wang, X.; Dong, W.; Yuan, B.; Yang, Y.; Yang, D.; Lin, X. I.; Chen, C.; Zhang,
553
W., Vitamin E confers cytoprotective effects on cardiomyocytes under conditions of
554
heat stress by increasing the expression of metallothionein. International Journal of
555
Molecular Medicine. 2016, 37 (5), 1429-1436.
556
(49) Sakatani, M.; Balboula, A. Z.; Yamanaka, K.; Takahashi, M., Effect of summer
557
heat environment on body temperature, estrous cycles and blood antioxidant levels in
558
Japanese Black cow. Animal science journal. 2012, 83 (5), 394-402.
559
(50) Kaur, H.; Halliwell, B. Action of biologically-relevant oxidizing species upon
560
uric acid. Identification of uric acid oxidation products. Chem. Biol. Interact. 1990, 73
561
(2-3), 235–247.
562
(51) Marklund, N.; Ostman, B.; Nalmo, L.; Persson, L.; Hillered, L. Hypoxanthine,
563
uric acid and allantoin as indicators of in vivo free radical reactions. Description of a
564
HPLC method and human brain microdialysis data. Acta Neurochirurgica. 2000, 142
565
(10), 1135-1142.
566
(52) Dietrich, M.; Traber, M. G.; Jacques, P. F.; Cross, C. E.; Hu, Y.; Block, G. Does
567
gamma-tocopherol play a role in the primary prevention of heart disease and cancer? A 28
ACS Paragon Plus Environment
Page 28 of 41
Page 29 of 41
568
Journal of Agricultural and Food Chemistry
review. Journal of the American College of Nutrition. 2006, 25 (4), 292–299.
29
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure captions Figure 1. Typical GC/MS chromatograms of the three beef cattle breeds from (A) serum and (B) urine. Figure 2. PLS-DA loading plots based on the serum and urine metabolic profilings. (A) The three breeds from serum (3 components, R2X= 0.387, R2Y=916, Q2 = 0.813). (B) XHC vs. SXC from serum (2 components, R2X= 0.352, R2Y=981, Q2 = 0.900). (C) XHC vs. JMY from serum (2 components, R2X= 0.344, R2Y=984, Q2 = 0.904). (D) SXC vs. JMY from serum (2 components, R2X= 0.347, R2Y=904, Q2 = 0.573). (E) The three breeds from urine (4 components, R2X= 0.579, R2Y=935, Q2 = 0.792). (F) XHC vs. SXC from urine (4 components, R2X= 0.658, R2Y=995, Q2 = 0.835). (G) XHC vs. JMY from urine (4 components, R2X= 0.658, R2Y=995, Q2 = 0.872). (H) SXC vs. JMY from urine (3 components, R2X= 0.562, R2Y=969, Q2 = 0.818). Figure 3. Different metabolic pathways from the three beef cattle breeds under heat stress. Red-colored symbols represent significant higher of metabolites in latter group compared with the former group, and black-colored symbols mean no difference, while blue-colored symbols indicate lower. Figure 4. Key enzymes related to the three beef cattle breeds in glycolysis, lipid metabolism and antioxidant defense system. (A) LDH activity. (B) MDA concentration and SOD, GSH-Px, T-AOC activity. * P < 0.05, ** P < 0.01.
30
ACS Paragon Plus Environment
Page 30 of 41
Page 31 of 41
Journal of Agricultural and Food Chemistry
Table 1. Identification of Different Metabolites from Any Two Groups of XHC, SXC and JMY in Serum XHC vs. SXC Metabolites
VIPa
Pb
XHC vs. JMY
FCc VIPa
Pb
FCc
SXC vs. JMY
VIPa
Pb
FCc
Metabolic pathway
Glucose
1.55 0.006 0.88 1.73 < 0.001 0.86
-
-
-
Glycolysis
Pyruvic acid
1.85 < 0.001 0.74 1.58 0.002
0.81
-
-
-
Glycolysis
Lactic acid
2.04 < 0.001 0.63 1.44 0.009
0.80
1.54 0.015 1.26 Glycolysis Pentose and glucuronate
Xylitol
-
-
-
-
-
-
2.00 < 0.001 0.71 interconversions Pentose and glucuronate
Glucuronic acid
-
-
-
-
-
-
1.58 0.007 0.79 interconversions Cysteine and methionine
Methionine
-
-
-
1.35 0.013
1.34
1.45 0.011 1.30 metabolism Cysteine and methionine
Cysteine
-
-
-
1.55 0.007
1.47
-
-
metabolism D-Glutamine and D-
Glutamic acid
1.38 0.025 1.50 1.50 0.011
1.64
-
-
glutamate metabolism D-Glutamine and D-
Glutamine
1.17 0.011 2.99 1.73 < 0.001 10.10 1.68 0.005 3.38 glutamate metabolism
2-amino-Butyric acid
1.33 0.028 1.95 1.49 0.009
2.25
-
-
-
Glutamate metabolism Phenylalanine and tyrosine
Phenylalanine
-
-
-
1.47 0.005
1.22
1.73 0.002 1.23 metabolism Phenylalanine and tyrosine
Tyrosine
-
-
-
1.27 0.025
2.27
1.29 0.035 2.19 metabolism Alanine, aspartate and
Asparagine
1.45 0.011 1.78 1.62 < 0.001 3.59
1.51 0.009 2.02 glutamate metabolism Arginine and proline
Ornithine
1.56 0.004 1.40 1.83 < 0.001 1.92
4-Hydroxyproline
1.40 0.006 1.41 1.24 0.030
1.24
1.71 0.003 1.37
-
-
-
metabolism Arginine and proline metabolism Arginine and proline
Creatinine
1.39 0.014 1.14 1.48 0.003
1.20
-
-
metabolism Arginine and proline
Putrescine
-
-
-
1.51 0.003
0.85
1.63 0.006 0.86 metabolism Arginine and proline
Urea
-
-
-
-
-
-
1.42 0.016 1.31 metabolism
Oleic acid
1.58 < 0.001 2.05
-
-
Linoleic acid
1.84 < 0.001 1.43 1.28 0.009
1.37
1.48 0.004 0.55 Fatty acid metabolism -
31
ACS Paragon Plus Environment
-
-
Fatty acid metabolism
Journal of Agricultural and Food Chemistry
Arachidonic acid
1.58 0.001 2.32 1.28 0.009
Palmitic acid
1.57 0.001 1.54
-
Tetradecanoic acid
1.24 0.012 1.55
-
Heptanoic acid
1.44 0.010 1.36 1.68 < 0.001 1.24
Glycerol
1.46 0.002 1.69 1.08 0.040
1.21
Glyceric acid
1.63 0.001 0.75 1.48 0.002
0.77
Cholesterol
-
-
-
myo-Inositol-1-phosphate
-
-
-
-
Page 32 of 41
1.32
1.42 0.014 0.57 Fatty acid metabolism
-
-
1.42 0.010 0.72 Fatty acid metabolism
-
-
1.77 < 0.001 0.78 Fatty acid metabolism
-
-
-
-
-
Fatty acid metabolism
1.27 0.021 0.71 Glycerolipid metabolism -
-
-
Glycerolipid metabolism
1.33 0.034 0.86 Steroid biosynthesis Inositol phosphate
1.00 0.050
0.67
1.53 0.005 0.51 metabolism Inositol phosphate
Myo-inositol
-
-
-
1.29 0.016
1.34
-
-
metabolism
Citric acid
1.23 0.032 1.23
-
-
-
Fumaric acid
1.92 < 0.001 1.85
-
-
-
1.99 < 0.001 0.54 TCA Cycle
-
-
-
1.30 0.031 0.86 TCA Cycle
α-Ketoglutaric acid
-
-
γ-Tocopherol
1.77 0.004 0.51 1.38 0.011
α-Tocopherol
1.71 0.001 0.71
-
-
TCA Cycle
0.68
-
-
-
Vitamin E metabolism
-
-
-
-
Vitamin E metabolism
1.27 0.013
3.71
-
-
-
Uric acid metabolism
Ribose
1.30 0.024 1.37 1.27 0.031
1.48
-
-
-
Purine metabolism
Oxalic acid
1.72 < 0.001 1.69 1.11 0.028
1.51
-
-
-
Others
1-Monooctadecanoylglycerol
1.21 0.025 1.30
Allantoin
aVIP
-
-
-
-
-
-
-
-
-
-
1.66 0.002 0.64 Others
Monomethylphosphate
-
-
-
1.67 < 0.001 0.66
1.52 0.008 0.73 Others
2-Oxoisocaproic acid
-
-
-
1.40 0.008
1.32
1.82 0.001 1.28 Others
Erythronic acid
-
-
-
1.48 0.004
1.23
-
-
-
Others
1,3-Di-tert-butylbenzene
-
-
-
1.31 0.013
0.90
-
-
-
Others
(variable importance in the projection) values were obtained from PLS-DA models with a threshold of 1. bP-
values were calculated from student's t-test with a threshold of 0.05. cFC (fold change) values were obtained from mean peak area of latter group/mean peak area of former group obtained. If the FC value is greater than 1, it means 32
ACS Paragon Plus Environment
Page 33 of 41
Journal of Agricultural and Food Chemistry
that metabolite level is higher in latter group compared with the former group.
33
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 34 of 41
Table 2. Identification of Different Metabolites from Any Two Groups of XHC, SXC and JMY in Urine XHC vs. SXC
XHC vs. JMY
SXC vs. JMY
Metabolites
Metabolic pathway VIPa
Pb
FCc VIPa
Pb
FCc VIPa
Lactic acid
1.49
0.037 0.70 1.56 0.014 0.63
Glycerol-3-phosphate
2.25
0.001 2.49
-
-
-
-
Pb
FCc
-
-
Glycolysis
1.76
0.012 0.57 Glycolysis 0.007 0.41 Glycolysis
Glyceric acid-3-phosphate
-
-
-
2.15 0.001 0.53 1.84
Fructose
-
-
-
1.45 0.039 0.49
-
-
-
Glycolysis
Galactose
-
-
-
1.37 0.044 0.65
-
-
-
Galactose metabolism
Xylitol
-
-
-
Pentose and glucuronate -
-
-
1.63
0.009 0.48 interconversions
2-Aminoadipic acid
-
-
-
-
-
-
1.48
0.022 0.51 Lysine degradation Arginine and proline
Urea
2.32 < 0.001 2.33 1.95 0.002 2.09
-
-
metabolism Arginine and proline
Putrescine
-
-
-
-
-
-
1.72
0.008 0.77 metabolism
Oleic acid
2.11
0.001 1.72 1.62 0.015 1.37
Linoleic acid
1.67
0.012 1.47
myo-Inositol
2.37 < 0.001 2.63 1.77 0.003 1.93 1.45
-
-
-
1.35
-
-
Fatty acid metabolism
0.049 0.78 Fatty acid metabolism Inositol phosphate 0.032 0.73 metabolism
Citric acid
-
-
-
1.55 0.041 0.70 1.80
0.004 0.53 TCA Cycle
Aconitic acid
2.04
0.002 1.54 1.76 0.008 0.66 2.24 < 0.001 0.43 TCA Cycle
Succinic acid
-
-
-
1.39 0.024 0.35
Methylmalonic acid
-
-
-
1.97 0.002 2.13 1.93
0.016 2.09 TCA Cycle
Methylcitric acid
-
-
-
1.57 0.023 1.82 1.70
0.005 1.95 TCA Cycle
-
-
-
TCA Cycle
Nicotinate and Nicotinic acid
1.87
0.004 1.62 1.61 0.006 1.97
-
-
nicotinamide metabolism
Uric acid
1.30
0.040 0.74 1.47 0.019 0.70
-
-
-
Purine metabolism
Allantoin
1.77
0.030 2.74 1.80 0.013 2.41
-
-
-
Uric acid metabolism
trans-Ferulic acid
2.07
0.004 0.59
-
-
-
Others
-
-
-
Others
1,3-Di-tert-butylbenzene
-
-
-
2,3-Dihydroxybutanedioic acid
-
-
-
aVIP
-
-
-
1.56 0.020 0.90 -
-
-
1.46
0.039 0.75 Others
(variable importance in the projection) values were obtained from PLS-DA models with a threshold of 1. bP-
values were calculated from student's t-test with a threshold of 0.05. cFC (fold change) values were obtained from 34
ACS Paragon Plus Environment
Page 35 of 41
Journal of Agricultural and Food Chemistry
mean peak area of latter group/mean peak area of former group obtained. If the FC value is greater than 1, it means that metabolite level is higher in latter group compared with the former group.
35
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure 1 A
B
36
ACS Paragon Plus Environment
Page 36 of 41
Page 37 of 41
Journal of Agricultural and Food Chemistry
Figure 2 A
B
C
D
E
F
37
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
G
H
38
ACS Paragon Plus Environment
Page 38 of 41
Page 39 of 41
Journal of Agricultural and Food Chemistry
Figure 3
39
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure 4 XHC
SXC
JMY
A
B
40
ACS Paragon Plus Environment
Page 40 of 41
Page 41 of 41
Journal of Agricultural and Food Chemistry
Table of Contents (TOC) Graphic
41
ACS Paragon Plus Environment