Nuclear magnetic resonance-based serum metabolic profiling of dairy cows with footrot.

Footrot is a debilitating and contagious disease in dairy cows, caused by the Gram-negative anaerobe Dichelobacter nodosus. 1H-NMR (nuclear magnetic resonance)-based metabolomics has been previously used to understand the pathology and etiology of several diseases. The objective of this study was to characterize serum from dairy cows with footrot (n=10) using 1H-NMR-based metabolomics and chemometric analyses. 1H-NMR spectroscopy with multivariate pattern recognition (principal component analysis and orthogonal partial least-squares discriminant analysis) was performed to identify biomarkers in cows with footrot (F) and healthy controls (C). 1H-NMR analysis facilitated the identification of 21 metabolites. Among these metabolites, 4 metabolites were higher and 17 metabolites were lower in the F group than in the C group. The serum levels of 5 metabolites were significantly different (P<0.05) between the two groups. The results revealed that cows with footrot have altered carbohydrate, amino acid, lipid and energy metabolic pathways. Metabolomic approaches are a clinically useful diagnostic tool for understanding the biochemical alterations and mechanisms of several diseases.

Footrot is a debilitating and contagious disease in ruminants that is caused by Dichelobacter nodosus, a Gram-negative anaerobe [6, 23]. Footrot, characterized by odorous exudative inflammation and keratinous hoof separation, is associated with anorexia and reduced milk yield in cows [17]. This painful condition causes lameness due to swelling and purulence of the foot, cracking and sloughing of interdigital skin, and puffiness of the coronary band. Infectious, immune, hormonal, nutritional and environmental factors are involved in footrot pathogenesis [16]. Footrot, which has an average worldwide incidence rate of 15–25% in dairy cows, causes significant financial losses in the dairy and beef cattle industries owing to expensive preventive and therapeutic interventions, reduced animal welfare and decreased productivity [37].

Current footrot diagnostic procedures include clinical examination, in vitro agglutination tests and polymerase chain reaction (PCR) [8]. If footrot control and management were effective during the non-transmission period, the disease could be completely eradicated. Currently, the complex pathogenesis of the disease poses considerable challenges to diagnosis and monitoring.

Fields, such as systems biology, genomics, transcriptomics, proteomics and metabolomics, have facilitated the development of novel diagnostic technologies. Serum proteomic analyses of cows infected with footrot represent an effective method for disease diagnosis [27]. However, proteomic analyses fail to reveal the complete metabolic profile of an animal. On the other hand, metabolomics measures the metabolic responses to genetic modification or pathophysiological stimuli [24]. Additionally, metabolomics complements information obtained from genetic and proteomic analyses. Currently, metabolomic techniques include nuclear magnetic resonance (NMR) spectroscopy, mass spectrometry and metabolite chips. Among these techniques, 1H-NMR has been frequently used to assess the role of metabolites in the body and their association with various pathological conditions. Recently, 1H-NMR has been used to investigate the effects of fatty liver disease, ketosis, milk fever and displaced abomasum on cow metabolism [3, 28, 29, 35]. However, the biochemical changes that occur in cows with footrot are unknown.

The main objective of this study was to use 1H-NMR-based metabolomics to assess metabolic differences between dairy cows with footrot and healthy controls. This study will provide useful information on footrot metabolite biomarkers.

MATERIALS AND METHODS

Ethics statement: This study was conducted in accordance with the Guide for the Care and Use of Laboratory Animals. The study was approved by the Ethics Committee of Heilongjiang Bayi Agricultural University (Permit Number: 20120319-1).

Animals: A commercial Holstein dairy cow farm in China suffering from an outbreak of footrot was chosen for this study. The cows in this study were at the same age and parity (1–3 births and 3–5 years old), in the middle stage of lactation. Cows acutely affected by footrot (n=10) were assigned to a test group (F group), while healthy cows (n=10) were assigned to the control group (C group). Footrot diagnosis was based on the presence of lameness in one or more feet and swollen or necrotic skin in the digits. Cows with severe lameness, commonly attributed to other diseases (e.g., foot and toe abscesses, white line disease and foot-and-mouth disease) were excluded from the study. Similarly, cows with mastitis, respiratory system diseases and digestive system diseases were excluded. Definitive footrot diagnosis was achieved by PCR (results not shown). To examine healthy cows, control samples from unaffected cattle in the affected dairy herd were pooled using the same procedure. All cows in this study were fed a total mixed ration (TMR) at 05:00, 14:00 and 20:00 hr. The TMR consisted of 55.60% dry matter (DM), 16.00% crude protein, 34.30% neutral detergent fiber, 22.00% acid detergent fiber, 5.60% fat, 1.07% calcium, 0.49% phosphorus, 0.32% magnesium, 0.13% sodium, 1.40% potassium, 0.39% chloride and 0.22% sulfur [28]. Table 1 shows the age, parity, dry matter intake (DMI) and body condition score (BCS) of cows in this study. The day before blood sample collection, the food intake of cows used for this study was recorded. BCS consists of a five-point scale (1 represents thin, and 5 represents obese) commonly used to assess body fat stores in Holstein dairy cows [34].

Table 1.

Clinical data and biochemical parameters of healthy cows (C) and cows with footrot (F)

No.	Parametersa)	C group	F group	P-values
1	Number	10	10	-
2	Age (yr)	3.98 ± 0.26	3.82 ± 0.24	0.651
3	Parity	1.90 ± 0.23	2.10 ± 0.23	0.550
4	BCS	3.30 ± 0.10	2.72 ± 0.11	0.001
5	DMI (kg/day)	18.07 ± 0.28	17.23 ± 0.31	0.058
6	Glc (mmol×l-1)	3.01 ± 0.28	2.68 ± 0.08	0.026
7	TG (mmol×l-1)	0.70 ± 0.04	0.68 ± 0.03	0.665
8	ALT (U×l-1)	20.07 ± 0.53	20.92 ± 0.58	0.297
9	AST (U×l-1)	70.76 ± 1.62	71.49 ± 1.80	0.770
10	TP (g×l-1)	62.54 ± 1.55	61.49 ± 1.20	0.581
11	NO (μmol×l1)	1.86 ± 0.18	4.59 ± 0.78	<0.001
12	MDA (nmol×l-1)	5.39 ± 0.22	6.23 ± 0.37	0.077
13	T-AOC (U×ml-1)	48.17 ± 1.19	40.74 ± 0.65	<0.001
14	SOD (U×ml-1)	108.66 ± 2.90	124.38 ± 2.31	<0.001
15	GSH-Px (U×l-1)	292.23 ± 2.29	245.81 ± 2.62	<0.001

a) Data are expressed as mean ± SEM where appropriate. Abbreviations: BCS (body condition score), Glc (glucose), TG (triglycerides), ALT (alanine aminotransferase), AST (aspartate aminotransferase), TP (total protein), NO (nitric oxide), MDA (maleic dialdehyde), T-AOC (total antioxidant capacity), SOD (activities of superoxide dismutase) and GSH-Px (glutathione peroxidase).

Sample collection and preparation: Following an overnight fast, blood samples (10 ml) were collected from the vena caudalis between 06:00 and 07:00 hr and centrifuged (3,500×g, 20 min at 4°C) to obtain the serum. The resulting serum samples were transferred into 1.5 ml Eppendorf tubes and stored at −80°C.

Biochemical analyses: Serum levels of glucose (Glc), triglycerides (TG), total protein (TP), aspartate aminotransferase (AST) and alanine aminotransferase (ALT) were measured with a Beckman Synchron CX5CE analyzer (Beckman Coulter, Inc., Fullerton, CA, U.S.A.). Commercial kits (Diagnostic Product, Nanjing Jiancheng Bioengineering Institute, Nanjing, China) were used to analyze serum nitric oxide (NO) level and activities of glutathione peroxidase (GSH-Px) and superoxide dismutase (SOD). All kits were utilized according to the manufacturer’s recommendations. Serum maleic dialdehyde (MDA) concentration was determined by the double heating method [10]. Serum total antioxidant capacity (TAC) was determined by performing a thiobarbituric acid reactive substances radical cation colorization assay as described by Guzel et al. [13].

Serum metabolic profiling by : Serum was allowed to thaw at room temperature and was centrifuged at 10,000×g for 10 min at 4°C. An aliquot (200 µl) of the resulting supernatant was mixed with 150 µl buffer (pH 7.4, 0.2 mol/l Na2HPO4 and 0.2 mol/l NaH2PO4) and 150 µl sodium 3-trimethylsilyl-(2,2,3,3-D4) propionate (TSP, 1 mg/ml, Sigma-Aldrich, St. Louis, MO, U.S.A.). Subsequently, the samples were centrifuged again at 10,000 ×g for 10 min at 4°C. Finally, 550 µl of the resulting supernatant was transferred into 5 mm NMR tubes containing 100 µl D2O.

NMR was performed at 25°C in a Bruker Avance spectrometer with a 1H operating frequency of 500 MHz (Bruker, Ettlingen, Germany). One-dimensional (1D) spectra were recorded with the Carr-Purcell-Meiboom-Gill (CPMG) sequence of D–[–90°–(τ–180°–τ) n–acquisition], with a total spin-echo relaxation delay (2 nτ) of 10 msec. Free induction decays (FIDs) were obtained with 128 scans of 32 K data points of 2.1 sec and a spectral width of 7,500 Hz [31]. The FID signal was measured with an exponential line broadening of 0.3 Hz and was Fourier-transformed to obtain the 1H-NMR spectra. The broad NMR signals of macromolecules like lipoproteins and proteins were reduced with a fixed total spin-echo delay (2nτ) of 40 msec. 1H-NMR spectra were manually phased and baseline-corrected using TopSpin software version 3.0 (Bruker Biospin, Rheinstetten, Germany).

Pre-processing of NMR spectral data: Raw NMR data were processed as previously described [9, 30]. The 1H-NMR spectra were corrected by the TSP signal into δH 0.00 and automatically converted into ASCII files using MestReNova Version 8.0.1 (Mestrelab Research SL, Santiago de Compostela, Spain). The ASCII files were imported into “R” software (http://cran.r-project.org/) and aligned using an open-access software package. The spectra were binned into consecutive non-overlapping regions of δ 0.003 chemical shift bins between δ 0.4–9.0 to avoid the encompassment of multiple peaks. The chemical shift at δ 0.0 was referenced to the TSP signal. Regions of residual water and EDTA signals between δ 4.2–9.0 were removed. Based on least square minimization, the remaining binned data (δ 0.4–4.2) of each NMR spectrum were auto-scaled and mean-centered to reduce multicollinearity and to simplify the interpretation of the modeling coefficients. For each spectrum, we calculated the integral values for different sample dilutions to reduce the effects of any concentration variability among the samples.

Multivariate pattern recognition: Multivariate statistical analyses of NMR data were performed by SIMCA-P 15.0 software (Umetrics, UmeÅ, Sweden). Unsupervised principal component analysis (PCA) and supervised orthogonal partial least-squares discriminant analysis (OPLS-DA) were used to determine differences between groups.

Based on a Hotelling’s T2 range plot, PCA was performed on the mean-centered data to identify possible outliers and to assess any intrinsic clustering from the two groups [36]. A PCA model was performed with the first two principal components (PC1 and PC2) in the scores plot. Through an orthogonal filter, the OPLS-DA model was used to identify the most significant variations between the two groups and remove unrelated systematic variations. A 10-fold cross-validated method was performed in the OPLS-DA model. In the scores plot of OPLS-DA, t[1]P represents the first principal component, and t[2]O represents the second orthogonal component. R2 and Q2, which represent the predicted ability of the model and the explained variability, respectively, were used to describe the quality of the OPLS-DA model. A color-coded coefficient plot of the OPLS-DA model was generated to identify differential metabolites between the two groups. In color-coded coefficient plots, red signals in the spectra corresponded to metabolites with greater contribution to separation between the two groups than blue signals. Signals with a positive or negative direction of t[1]P corresponded to metabolites that were present at high or low concentrations, respectively.

Identification of metabolites and metabolic pathway analysis: Based on the chemical shifts, metabolites in the two groups were identified with the Chenomx NMR Suite 7.5 library (Chenomx Inc., Edmonton, Canada). The identified metabolites were defined based on data from a previous study and public databases, e.g., the Human Metabolome Database (HMDB, http://www.hmdb.ca) and the Biological Magnetic Resonance Data Bank (BMRB, www.bmrb.wisc.edu) [24]. Serum metabolites from the two groups are presented in Table 2. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.kegg.jp), the impact factor of the topology analyses of metabolic pathways was graphically presented using metabolomics data analysis software (Metaboanalyst 2.0, www.metaboanalyst.ca/MetaboAnalysts) [32]. The results were presented graphically as a bubble plot. Darker colored, larger areas of the bubbles represent more significant metabolite changes in the corresponding pathway.

Table 2.

Assignment of the identified metabolites in serum of healthy cows and cows with footrot

No.	Metabolites	Chemical shift δ 1Ha)	Expression change b)
1	3-Hydroxy-3-methylglutarate	1.32 (s), 2.39(d), 2.45(d)	↑
2	3-Hydroxybutyrate	1.20 (d), 2.30 (dd), 2.40 (dd), 4.13 (dd)	↑
3	3-Hydroxyisobutyrate	1.05 (d), 2.46 (m), 3.53 (q), 3.69 (q)	↓
4	Ethanol	1.19 (t), 3.67 (q)	↓
5	Acetoacetate	2.28 (s), 3.44 (s)	↓
6	Alanine	1.49 (d), 3.78 (dd)	↓
7	Acetone	2.24 (s)	↓
8	Carnitine	2.40 (q), 2.44 (dd), 3.21 (s), 3.38 (t), 3.43 (t)	↓
9	Pyruvate	2.38 (s)	↓*
10	Sarcosine	2.74 (s), 3.61 (s)	↓*
11	Succinate	2.39 (s)	↓**
12	Betaine	3.25 (s), 3.89 (s)	↓
13	Creatine	3.04 (s), 3.93 (s)	↓
14	1-Methylhistidine	3.66 (s)	↓
15	Glycine	3.55 (s)	↓
16	Glucose	3.26 (dd), 3.41 (t), 3.45 (m), 3.8 (t)	↓
17	Glycerate	3.7 (dd)	↓
18	Lactate	1.33 (d), 4.12 (q)	↑**
19	Serine	3.84 (dd), 3.96 (m)	↓
20	Valine	0.99 (d), 1.02 (d), 3.62 (d)	↓**
21	3-Methylhistidine	3.71 (s), 3.92 (dd)	↓

a) Multiplicity: s, singlet peak; d, doublet peaks; dd, double doublet peaks; t, triplet peaks; q, quartet peak; m, multiple peaks. b) “↑” and “↓” represent increased and decreased levels in the footrot cows compared with control cows, respectively. “*” and “**” represent significant differences (P<0.05) and highly significant differences (P<0.01), respectively.

Univariate analysis: Clinical data, serum biochemistry data and differential metabolites between the two groups were analyzed by one-way ANOVA using SPSS statistical software Version 11.0 (SPSS, Inc., Chicago, IL, U.S.A.). Data were presented as mean ± standard error of the mean (SEM). P<0.05 was considered statistically significant.

RESULTS

Clinical information and serum biochemistry analysis: The clinical data and serum biochemical differences between cows with footrot (F group) and the healthy controls (C group) are shown in Table 1. There were no differences in age, parity or DMI between the two groups (P>0.05). On the other hand, BCS was significantly different between F and C groups (P<0.05). Serum Glc, NO, T-AOC, SOD and GSH-Px were significantly different between the two groups (P<0.05). No significant differences were observed for TG, ALT, AST, TP or MDA levels (P>0.05).

Metabolic profiles of serum samples from dairy cows: The 1H-NMR (δ 0.4–4.2) spectra revealed a wide range of metabolites. Typical CPMG 1H-NMR profiles of sera obtained from the F and C groups are presented in Fig. 1. Twenty-one metabolites were identified and confirmed by public databases and two-dimensional NMR techniques (heteronuclear single quantum correlation and total correlation spectroscopy). The endogenous metabolites included glucose, ketone bodies (acetone, 3-hydroxybutyrate, 3-hydroxy-3-methylglutarate, 3-hydroxyisobutyrate and acetoacetate)and amino acids (alanine, serine, valine, glycine, sarcosine, 1-methylhistidine and 3-methylhistidine), among others. The main differential metabolites were identified through multivariate data analysis.

Fig. 1.

Typical 500 MHz 1H-NMR spectrum of serum samples from footrot (F, top spectrum) and control (C, bottom spectrum) and footrot (F, spectrum) groups. A total of 21 metabolites were assigned. Metabolites: 1, 3-hydroxy-3-methylglutarate; 2, 3-hydroxybutyrate; 3, 3-hydroxyisobutyrate; 4, ethanol; 5, acetoacetate; 6, alanine; 7, acetone; 8, carnitine; 9, pyruvate; 10, sarcosine; 11, succinate; 12, betaine; 13, creatine; 14, 1-methylhistidine; 15, glycine; 16, glucose; 17, glycerate; 18, lactate; 19, serine; 20, valine; 21, 3-methylhistidine.

Determining metabolic variations between groups using multivariate analysis: First, a PCA scores plot with pareto-scaled data was generated and analyzed (Fig. 2). PC1 and PC2 accounted for 43.91% and 14.49% of the total variation in the dataset, respectively. The OPLS-DA model was performed to discriminate samples according to their groups (footrot or control). The results are shown in Fig. 3. The OPLS-DA score plot (Fig. 3a) demonstrated a reasonable and considerable separation between the two groups (R2=0.79; Q2=0.75), suggesting a significant metabolic perturbation in the footrot group. To identify the 21 serum metabolites that contributed to group separation, OPLS-DA color-coded loadings plots (Fig. 3b and c) were generated. For each variable, information was provided by different colors (red representing high values; blue representing low values), based on the absolute correlation coefficient (|r|). The significant metabolites in the negative regions of the loading plot (Fig. 3b and 3c) were increased in cows with footrot. Conversely, the significant metabolites in the positive regions of the loading plot were decreased in cows with footrot. The results of OPLS-DA revealed that 21 different metabolites were present (Table 2). The OPLS-DA loading plot showed that cows with footrot had higher levels of 3-hydroxy-3-methylglutarate, 3-hydroxybutyrate and lactate, and lower levels of 3-hydroxyisobutyrate, ethanol, acetoacetate, alanine, acetone, carnitine, betaine, creatine, isoleucine, glycine, glucose, glycerate, serine, 3-methylhistidine, pyruvate, sarcosine, succinate and valine, compared to the control cows. The P-values of the 21 metabolites were obtained by t-tests (Table 2). Levels of pyruvate and sarcosine (P<0.05), and lactate, succinate and valine (P<0.01) were significantly different between cows from F and C groups.

Fig. 2.

PCA score plot of 1H-NMR spectra of serum samples from cows with footrot (box, ■) and healthy controls (dot, ●). Each point represents an individual sample on the score plot. The center and margin of the ellipse represent mean and standard deviation, respectively.

Fig. 3.

OPLS-DA comparison of serum spectra between the C and F groups. (a) OPLS-DA scores plot obtained from 1H-NMR CPMG spectral data for cows with footrot (F, box, ■) and controls (C, dot, ●), where one point represents one sample. (b and c) Corresponding loading plots (δ 0.4–2.5 and 2.6–4.2) color-coded according to correlation coefficient from blue to red. The upward and downward peaks represent decreased and increased levels in cows with footrot, respectively. The metabolites are presented in Table 2.

Identification of impacted metabolite pathways: The MetaboAnalyst 2.0 software was used to provide impact factors from the pathway topology analysis for all analyzed pathways. The 28 different metabolic pathways of the obtained serum metabolites were mapped (Fig. 4). The top 10 pathways for footrot included glycine, serine and threonine metabolism, ketone body synthesis and degradation, methane metabolism, valine, leucine and isoleucine biosynthesis, pyruvate metabolism, aminoacyl-tRNA biosynthesis, glycerolipid metabolism, butanoate metabolism, glycolysis or gluconeogenesis, and the TCA cycle (Table 3).

Fig. 4.

The bubble plot shows a metabolomic overview of altered metabolic pathways, mapping serum metabolites identified in footrot and control groups. Keys: (a) glycine, serine and threonine metabolism; (b) synthesis and degradation of ketone bodies; (c) methane metabolism; (d) valine, leucine and isoleucine biosynthesis; (e) pyruvate metabolism; (f) aminoacyl-tRNA biosynthesis; (g) glycerolipid metabolism; (h) butanoate metabolism; (i) glycolysis or gluconeogenesis; and (j) citrate cycle (TCA cycle).

Table 3.

Significantly altered metabolic pathways in cows with footrot

No.	Pathway	Metabolitesa)	Hitsb)	-log (g)	Impact
1	Glycine, serine and threonine metabolism	32	7	15.66	0.596
2	Synthesis and degradation of ketone bodies	5	2	6.16	0.600
3	Methane metabolism	9	2	4.92	0.400
4	Valine, leucine and isoleucine biosynthesis	11	2	4.51	0.338
5	Pyruvate metabolism	22	2	3.17	0.189
6	Aminoacyl-tRNA biosynthesis	64	3	2.68	0.140
7	Glycerolipid metabolism	18	1	1.42	0.105
8	Butanoate metabolism	20	4	8.74	0.101
9	Glycolysis or Gluconeogenesis	26	3	5.09	0.099
10	Citrate cycle (TCA cycle)	20	2	3.35	0.098

a) Total number of compounds in pathways based on KEGG databases. b) Matched number of identified metabolites in serum from F and C groups.

DISCUSSION

Dairy cows with acute footrot can cause considerable economic loss through a reduction in milk production [21]. If not treated in time, the infection can lead to blood poisoning and be life threatening [19]. Footrot is a multifactorial disease related to immune, hormonal, nutritional and environmental factors as well as infection [16]. However, a comprehensive study of the etiological effects of these causative factors in the pathogenesis of footrot has not yet been performed. Metabolomics provides a novel pathway for resolving this question. 1H-NMR-based metabolomics represents a powerful tool for identifying serum metabolites in ruminant animals [3, 22]. Further, preliminary studies by our group suggested that 1H-NMR-based metabolomics can assess ketosis and milk fever in dairy cows and identify metabolites [28, 29]. However, little is known about the metabolic impact of footrot in dairy cows. In this study, an 1H-NMR-based metabolomics approach combined with biochemical assays was used to comprehensively assess serum metabolites in dairy cows with and without footrot.

Antioxidants play an important physiological role in counteracting free radicals and preventing cellular damage. The TAC was measured to determine the cumulative action of all the antioxidants present in blood. A significant decrease (P<0.001) in TAC level was detected in cows with footrot compared to controls. The GSH-Px activity was significantly lower (P<0.001) in the F group compared to the controls, which suggests a reduction in antioxidative protection and superoxide radical scavenging action in cows with footrot. The high level of serum SOD seen in the F group may be a compensatory response to protect the body from oxidative damage. These observations may be because of alterations in the antioxidative and oxidative balance, which is in accordance with a previous study [1]. Activation of the immune system by footrot infection will release reactive oxygen and reactive nitrogen species, inducing oxidative stress and lipid peroxidation, which affects cow development [14]. MDA is produced by lipid peroxide, which is used as a marker for oxidative damage induction [20]. Higher levels of serum MDA indicated lipid peroxidation in cows with footrot. Large quantities of NO, synthesized by infectious disease-activated macrophages, play a major role in defense mechanisms [15]. NO also may cause tissue damage during inflammation, through the formation of peroxynitrite. Peroxynitrite is a powerful oxidant capable of initiating lipid peroxidation [4]. Significantly elevated NO levels in cows with footrot (P<0.001) suggest that lipid peroxidation is induced in this disease. Furthermore, liver damage is often induced under conditions of lipid peroxidation, as shown by a slight rise in serum levels of AST and ALT [32]. Higher AST levels were observed in cows with footrot, which may be attributed to an acute dystrophic muscle condition as described by a previous study [38].

Combined with statistical models, 1H-NMR identified alterations in specific metabolic pathways of cows with footrot. Footrot affected several metabolic pathways, including amino acid metabolism, carbohydrate metabolism, lipid metabolism and energy metabolism.

Alterations in amino acid metabolism: Cows with footrot had significant disturbances in amino acid metabolism. This was evident from the low levels of valine, betaine, glycine, alanine, creatine, serine, sarcosine, carnitine, ethanol, 1-methylhistidine and 3-hydroxyisobutyrate. Past studies by our group showed that the blood levels of 1-methylhistidine, 3-hydroxyisobutyrate, alanine, betaine, creatine and ethanol were also lower in dairy cows with ketosis, milk fever or fatty liver than in controls [28,29, 35]. Meanwhile, the blood levels of glycine and valine in cows with footrot and other diseases (ketosis and fatty liver) showed the opposite trend. Valine, glycine, 3-hydroxyisobutyrate (an intermediate of valine metabolism) and serine are gluconeogenic amino acids that produce oxaloacetate [12]. Alanine, a non-essential amino acid that can be synthesized from pyruvate, is one of the most important energy sources released by muscle tissues. Footrot often affects animal physiology in many ways by causing mild fever, lameness, severe pain and a resulting drop in feed intake. However, this study was performed on acute cases of footrot in cows. Compared with the control cows, the DMI showed an insignificant downward trend (P=0.058) in cows with footrot. Based on the results, the lack of gluconeogenic substrates may play an important role in the pathogenesis of acute footrot. Creatine increases the synthesis of adenosine triphosphate (ATP) in all body cells. Additionally, creatine plays a vital role in the urea cycle and in ATP homeostasis [33]. Decreased creatine levels suggest lack in energy of cows with footrot.

Sarcosine, an intermediate in the conversion of choline to glycine, is a protein constituent [2]. Additionally, sarcosine is a prime metabolic source of components like creatine and serine for living cells. Serum biochemistry analysis showed that the level of TP in F group serum was slightly decreased compared to the level in C group serum. Low levels of serum sarcosine in cows with footrot indicated high protein utilization.

1-Methylhistidine, which is converted into histidine, is closely linked to inflammation, oxidative stress and protein-energy wasting. 3-Methylhistidine, a product of peptide bond synthesis actin and myosin methylation, is an indicator of muscle protein catabolism. Previous studies have shown that 3-methylhistidine is a clinical marker of nitrogen loss, occurring in conditions, such as trauma and infection [26]. Changes of 1-methylhistidine and 3-methylhistidine levels in cows with footrot suggest insufficient dietary protein intake and high oxidative stress.

Alterations in carbohydrate metabolism: Compared with healthy controls, cows with footrot had reduced glucose and pyruvate levels. Glucose plays an important role in oxidation/reduction, and pyruvate is a key compound in the network of metabolic processes [5]. Glucose produces ATP via glycolysis, a pathway in which glucose is converted into pyruvate. Under aerobic conditions, pyruvate is subsequently converted into acetyl CoA, the compound that initiates the TCA cycle. Cows with footrot have reduced ATP synthesis from glycolysis, glycogenolysis and the TCA cycle. On the other hand, lactate is produced as a waste product during the process of converting glucose into pyruvate. Lactate is metabolized to form glucose, by the process known as gluconeogenesis. However, because of enzymes’ high affinity for glucose, a high rate of glycolysis was maintained even in the presence of oxygen, leading to increased lactate levels [7, 38]. In the present study, higher serum lactate levels in cows with footrot may be related to the inhibition of glycolysis and gluconeogenesis, indicating inefficient metabolism [37]. In contrast, decreased lactate levels were found in our study on cows with ketosis and fatty liver, which indicated an increase in gluconeogenesis [28, 35].

Phosphate derivatives of glycerate, such as 1,3-bisphosphoglyceric acid, 2-phosphoglyceric acid, 2,3-bisphosphoglyceric acid and 3-phosphoglyceric acid, are important biochemical intermediates of glycolysis. In humans, lack of glycerate kinase leads to D-glyceric aciduria and D-glycerate anemia, which contribute to failure to thrive, progressive neurological impairment, seizures, hypotonia and metabolic acidosis [25]. In this study, cows with footrot had low levels of serum glycerate, which revealed the presence of a carbohydrate metabolic disorder.

Succinate is an important biochemical intermediate, playing an important role in the TCA cycle. Furthermore, succinate is capable of donating electrons to the electron transport chain [32]. Low levels of serum succinate in cows with footrot demonstrated a reduced flux through the TCA cycle. However, we found that blood succinate levels were increased in a fatty-liver study by our group. This result suggests that the TCA cycle is upregulated in body glucose metabolism [28].

Alterations in lipid metabolism: Several studies have shown that carnitine plays an important role in fatty acid metabolism [11]. Carnitine increases fatty acid transport into mitochondria, thereby increasing pyruvate oxidation. Decreased carnitine levels were observed in cows with footrot, indicating enhanced fatty acid transport into mitochondria.

Ketone bodies (i.e., acetoacetate, 3-hydroxybutyrate and acetone) derived from acetyl-CoA are byproducts of fatty acid oxidation [28]. Cows with footrot had reduced levels of acetoacetate and acetone, and high levels of 3-hydroxybutyrate, indicative of increased lipolysis. As a hydroxy fatty acid, 3-hydroxy-3-methylglutarate inhibits cholesterol synthesis in non-insulin-dependent diabetics [18]. In this experiment, cows with footrot had elevated levels of 3-hydroxy-3-methylglutarate.

1H-NMR metabolomic profiling was performed on serum to explore metabolic alterations in cows with footrot. Footrot affected amino acid, carbohydrate and lipid metabolic processes. The serum levels of four metabolites were found to be significantly different (P<0.05) between cows with footrot and healthy controls. These metabolites may be useful for diagnosing acute or continuous impaired nutritional conditions in cows with footrot. Metabolomics provides a novel viewpoint and technological approach to link disease processes with biochemical pathways. 1H-NMR techniques can reveal the functions of different body systems. In summary, serum metabolite profiling demonstrated significant differences between cows with footrot and healthy controls.