Skip to main content
  • Research article
  • Open access
  • Published:

Plasma metabolomic profiling of dairy cows affected with ketosis using gas chromatography/mass spectrometry



Ketosis is an important problem for dairy cows` production performance. However, it is still little known about plasma metabolomics details of dairy ketosis.


A gas chromatography/mass spectrometry (GC/MS) technique was used to investigate plasma metabolic differences in cows that had clinical ketosis (CK, n=22), subclinical ketosis (SK, n=32), or were clinically normal controls (NC, n=22). The endogenous plasma metabolome was measured by chemical derivatization followed by GC/MS, which led to the detection of 267 variables. A two-sample t-test of 30, 32, and 13 metabolites showed statistically significant differences between SK and NC, CK and NC, and CK and SK, respectively. Orthogonal signal correction-partial least-square discriminant analysis (OPLS-DA) revealed that the metabolic patterns of both CK and SK were mostly similar, with the exception of a few differences. The development of CK and SK involved disturbances in many metabolic pathways, mainly including fatty acid metabolism, amino acid metabolism, glycolysis, gluconeogenesis, and the pentose phosphate pathway. A diagnostic model arbitrary two groups was constructed using OPLS-DA and receiver–operator characteristic curves (ROC). Multivariate statistical diagnostics yielded the 19 potential biomarkers for SK and NC, 31 for CK and NC, and 8 for CK and SK with area under the curve (AUC) values. Our results showed the potential biomarkers from CK, SK, and NC, including carbohydrates, fatty acids, amino acids, even sitosterol and vitamin E isomers, etc. 2-piperidinecarboxylic acid and cis-9-hexadecenoic acid were closely associated with metabolic perturbations in ketosis as Glc, BHBA and NEFA for dealing with metabolic disturbances of ketosis in clinical practice. However, further research is needed to explain changes of 2,3,4-trihydroxybutyric acid, 3,4-dihydroxybutyric acid, α-aminobutyric acid, methylmalonic acid, sitosterol and α-tocopherol in CK and SK, and to reveal differences between CK and SK.


Our study shows that some new biomarkers of ketosis from plasma may find new metabolic changes to have clinically new utility and significance in diagnosis, prognosis, and prevention of ketosis in the future.


Ketosis is one of the most prevalent metabolic diseases of dairy cows during transition period [1]. Considering the serious consequences of ketosis such as fatty liver, abomasum displacement, infectious diseases, and reproductive diseases, this disease has become a health concern during the past decades [2, 3]. Therefore, elucidating the pathogenesis of ketosis is of great importance to the monitoring, prevention, and early treatment of ketosis. To this end, there have been many reports about clinical pathological changes in dairy ketosis. β-hydroxybutyrate (BHBA) is widely considered the golden standard for diagnosing ketosis in dairy cows. In addition, blood glucose levels (Glc), total triglycerides (TG), nonesterified fatty acids (NEFA), and aspartate aminotransferase (AST) can also be analyzed to monitor ketosis-related complications [35]. However, to the best of our knowledge, the plasma metabolomic profiling of ketosis has not yet been clarified.

Metabolomics technology has proven to be a powerful tool for biomarker screening, disease diagnosis, and characterization of biological pathways in humans, rats, and cattle [610]. In view of the high efficiency of chromatographic separation and the sensitive detection of separated components, gas chromatography/mass spectrometry (GC/MS) is widely employed for the plasma metabolomic profiling of diseases, in combination with multivariate statistical analysis including principal component analysis (PCA), partial least square-discriminant analysis (PLS-DA), and orthogonal partial least square-discriminant analysis (OPLS-DA), etc. [11, 12].

Although milk and plasma metabolome of ketosis has been reported using NMR technique [6, 7], it is little known about metabolome profiles of clinical and subclinical ketosis. In this study, a metabolomics approach by GC/MS in conjunction with multivariate statistical analysis was developed to show alterations in plasma metabolite profiles of cows associated with normal or enhanced ketone body formation. Plasma samples from clinically normal cows or those with subclinical or clinical ketosis were profiled by GC/MS. Receiver operator characteristic curves (ROC) and OPLS-DA, which finds new variables with optimal discriminatory ability and small redundancy, were used to explore perturbations in metabolic patterns and potential biomarkers for clinical and subclinical ketosis.

Results and discussions

Comparison of demographics, clinical information and blood chemical parameters

Clinical data from normal cows (NC) or those with subclinical (SK) or clinical ketosis (CK) were collected to uncover any biochemical differences among the three groups shown in Table 1. Postpartum days, ages, parity, daily milk yield (MY), and body condition score (BCS) were comparable among the three groups of cows. There were significant differences in MY and BCS plasma Glc and BHBA concentrations among the three groups of cows, but not in plasma AST, CHO, and TG. Our results indicated that MY and BCS decreased significantly in dairy cows with ketosis, which is consistent with some reports that the affected cows often loss more milk yield and body weight [13, 14].

Table 1 Main characteristics and blood biochemical parameters in normal control cows (NC), cows with subclinical ketosis (SK), and cows with clinical ketosis (CK)

All affected cows mainly presented low MY, loss of BCS, low blood Glc levels and high blood BHBA concentrations, which implies poor glucose regulation and lipid profile disorders. Although the above mentioned clinical biochemistry, such as Glc, BHBA, and AST levels, may to some extent reflect the metabolic state of ketosis, they only give limited information, as there are various metabolic abnormalities that can exist during ketosis progression, as well as other underlying metabolic changes associated with ketosis that can result in serious metabolic outcomes such as fatty liver and infectious diseases [13]. Thus, metabolomics may provide a more global and powerful approach for understanding the metabolic patterns and potential biomarkers of ketosis, and additionally, may shed light on the pathogenesis of ketosis in dairy cows.

GC/MS metabolomic profiling

Figure 1 shows the three representative chromatograms of derived plasma samples from the three groups of cows, which had almost all same peaks among the three chromatograms, but there was also difference in peaks’ area among them. After mean-centering, UV-scaling, and deconvolution with AMDIS software, 40 of 267 variables were only significantly different between groups by multivariate statistical analysis, and then were identified as endogenous metabolites by NIST mass spectra library, and were categorized as carbohydrates, amino acids, fatty acids, sitosterol and vitamin E isomers, etc. (Table 2). Each of these is involved in multiple biochemical processes, especially those processes related to energy balance and lipid metabolism.

Figure 1
figure 1

The representative GC/MS chromatograms of plasma samples from No. 8 of normal control (NC, purple peaks) cows, No. 35 of cows with subclinical ketosis (SK, brown peaks) and No. 68 of cows with clinical ketosis (CK, green peaks). X-axis is retention time (min); Y-axis is intensity of MS (%).

Table 2 Marker metabolites found in the GC/MS chromatograms of any two groups of CK, SK, and NC

Since the metabolic profiles of CK and SK were similar, multivariate statistical analyses were used to explore biomarker candidates and disturbances in metabolic patterns of ketosis. PCA was used to examine the clustering of samples for differences in metabolic patterns among the three groups. The PCA score plot could not differentiate between CK, SK, and NC because the samples from the three groups scattered into each other; however, they could be successfully discriminated using the PLS-DA model (data not shown). Subsequently, OPLS-DA, which reduces dimensionality of the original data, was applied to explore metabolic disturbances in, NC, SK, and CK. The plot of PC1 (first principle component, t[1]P) versus O (orthogonal component, t[2]O) is shown in Figure 2. As illustrated in Figure 2, samples from the CK group lay on the left side of t[1]P, while samples from SK and NC lay in the middle and on the right side of t[1]P. In addition, the parameters of OPLS-DA, R2Y and Q2, were 0.828 and 0.693, respectively. R2Y showed the explanative ability of the model, and Q2 was the result of a seven-fold cross-validation, and showed the predictive ability of the model for metabolic profiling of the data [15, 16]. These results indicate that the proposed metabolic process of ketosis was from NC to SK, and then to CK.

Figure 2
figure 2

Score plot of OPLS-DA among three groups after excluding three outlier samples from Group CK, SK and NC respectively.

To further investigate the metabolic status and potential biomarkers of two groups, samples from two groups were subjected to an OPLS-DA model. As can be seen in Figures 3 (a, b, c), samples from the corresponding two groups were distinctly separated on two sides of the PC1 baseline. Furthermore, the OPLS-DA parameters, R2Y and Q2, were 0.819 and 0.74 for CK vs. NC, 0.856 and 0.595 for SK vs. NC, and 0.931 and 0.301 for CK vs. SK, respectively. The results of the permutation test for the R2 and Q2 intercepts were 0.645 and −0.279 for CK vs. NC, 0.292 and −0.144 for SK vs. NC, and 0.323 and −0.105 for CK vs. SK, respectively. Thus, OPLS-DA was suitable for exploring differences between two groups in this study.

Figure 3
figure 3

The OPLS-DA scores plot and ROC plot of plasma obtained from arbitrary two groups. (a, b, c) Score plot of OPLS-DA of using Principal Component 1 (PC1) in Group CK vs. NC, Group SK vs. NC, and Group CK vs. SK respectively; (d, e, f) ROC plot of cis-9-hexadecenoic acid, 2-piperidinecarboxylic acid, and 3-hydroxybutyric acid in the diagnosis in Group CK vs. NC, Group SK vs. NC, and Group CK vs. SK respectively.

Metabolic disturbances of ketosis

Table 2 lists 21 decreased metabolites from the three groups of cows, the majority of which were amino acids, carbohydrates, and their metabolites; the minority of which were carboxyl acids and hydroxyl acids from fatty acids metabolism and the tricarboxylic acid (TCA) cycle. Meanwhile it also shows the 19 increased metabolites, the majority of which were non-esterified fatty acids (NEFA), amino acids, and their metabolites; the minority of which were carbohydrates, and even sitosterol and vitamin E isomers, etc. These results suggest that metabolic disturbances of ketosis involve in multi-biochemical pathways such as glycolysis, gluconeogenesis, amino acids metabolism, fatty acids metabolism, pentose phosphate pathway.

In the current study, potential metabolites from arbitrary two groups were found to be significantly changed based on the largest VIP and a students t-test (P < 0.05). Table 2 summarizes 40 potential biomarkers that were differentially found in the CK, SK, and NC groups, which 25 metabolites were the same in both CK and SK, including 7 metabolites in CK vs. SK, indicating that CK and SK have rather similar pathogenesis.

Figure 4a shows that nine of the 25 metabolites; namely, lactic acid (LA), glucuronic acid (GLCA), l-alanine (L-ala), glycolic acid (GA), ribitol, pyroglutamic acid (pGlu), galactose (Gal), 2,3,4-trihydroxybutyric acid (THBA), and glucose (Glc), decreased from low to high in both CK and SK compared to NC. With the exception of L-ala in CK vs. SK and Gal in SK vs. NC, fold changes of other decreased metabolites are in turn low from CK vs. NC to SK vs. NC, and CK vs. SK. LA, L-ala, and pGlu from glucogenic amino acids and THBA from glucogenic and ketogenic threonine, which may be converted into glucose in liver by gluconeogenesis [1720]. Low pGlu values may indicate glutathione deficiency due to oxidative stress [17]. Thus, gluconeogenesis decreased in both CK and SK because of a decrease in these glucogenic amino acids. In addition, GLCA, ribitol, Gal, and Glc can play a role in energy metabolism by entering glycolysis via the pentose phosphate pathway [2124]. Decreased ribitol levels may be related to riboflavin deficiency [22]. Vitamin C may also be involved in energy metabolism due to GLCA, which is derived from vitamin C [21]. Therefore, a low level of plasma carbohydrates may play an important role in the development of ketosis, and low vitamin C and B2 levels, as well as oxidative stress, may also be contributing factors. It is also worth noting that a decrease in GA may reduce milk components in both CK and SK, since GA is used to synthesize milk fat or lactose [25]. GLCA and LA may play important roles in maintaining energy balance because they were maximally down-regulated in both CK and SK. In brief, our study suggests that ketosis is closely associated with disturbances in carbohydrate metabolism due to hypoglycemia and a lack of glucogenic precursors such as LA and L-ala. Vitamin C and B2 deficiencies, oxidative stress, and low GLCA and LA, may also contribute to the development of this disease.

Figure 4
figure 4

Intensity of metabolites from CK compared to NC (blue column), SK compared to NC (red column), and CK compared to SK (yellow column). (a) Expression intensity of 25 identical metabolites from CK vs. NC and SK vs. NC, including 7 metabolites from CK vs. SK. 2PC, 2-piperidinecarboxylic acid; C-9-HA, cis-9-hexadecenoic acid; MA, myristic acid; T-9-OA, Trans-9-octadecenoic acid; BHBA, 3-hydroxybutyric acid; SA, stearic acid; HA, heptadecanoic acid; PA, palmitic acid; 3HV, 3-hydroxyvaleric acid; HMA, 3-hydroxy-3-methylglutaric acid; VE, vitamin E (α-tocopherol); AMA, aminomalonic acid; AABA, α-aminobutyric acid; L-ile, l-isoleucine; Gly, glycine; Glc, Glucose; THBA, 2,3,4-trihydroxybutyric acid; Gal, galactose; pGlu, pyroglutamic acid; GA, glycolic acid; L-ala, l-alanine; GLCA, glucuronic acid; LA, lactic acid. (b) Expression intensity of 15 different metabolites from three groups. GABA, 4-aminobutyric acid; L-ser, l-serine; 4HYP, 4-hydroxyproline; 3-HIV, 3-hydroxyisovaleric acid; 2Me3HB, 2-methyl-3-hydroxybutyric acid; L-orn, l-ornithine; MMA, methylmalonic acid; 2-KG, 2-ketoglutaric acid; L-leu, l-leucine; Cit, citrate; DHB, 3,4-dihydroxybutyric acid; L-pro, l-proline.

Figure 4a also shows that 16 of 25 metabolites increased in both CK and SK compared to NC. With the exception of aminomalonic acid (AMA), α-aminobutyric acid (AABA), l-isoleucine (L-ile), and glycine (Gly), fold changes of other metabolites increased in turn from CK vs. NC to SK vs. NC, and CK vs. SK. These metabolites were mainly increased 3-hydroxybutyric acid (BHBA) and nonesterified fatty acids (NEFAs), including palmitic acid (PA), heptadecanoic acid (HA), stearic acid (SA), trans-9-octadecenoic acid (T-9-OA), myristic acid (MA), cis-9-hexadecenoic acid (C-9-HA), which belong to the families of ketone bodies, long chain unsaturated fatty acids, and saturated acids [1, 26], confirming that a great amount of fat mobilization resulting from hypoglycemia may cause ketosis. In addition, some up-regulated amino acids and their catabolic products, such as L-ile, a glucogenic and ketogenic amino acid [25]; Gly, biosynthesized from serine [27]; AMA, a constituent of proteins before hydrolysis [28]; and 2-piperidinecarboxylic acid (2PC), a metabolite of the lysine metabolism [29]; suggest that proteolysis increases to meet body energy demand in both CK and SK. Since 2PC up-regulation was highest in both CK and SK, it may play an important role in ketone body synthesis as a catabolic product of ketogenic lysine [29]. In addition, the up-regulation of other metabolites, such as 3-hydroxyvaleric acid (3HV), 3-hydroxy-3-methylglutaric acid (HMG), and AABA, implicated abnormal metabolic changes or functional abnormalities in both CK and SK, since high 3HV concentrations can be found in methylmalonic acidemia [30], increased HMG values may be caused by decreased coenzyme Q10 synthesis [31], and elevated AABA levels can be indicative of liver disease [32]. However, it still remains unknown whether they play roles in the development of ketosis. Interestingly, we observed high sitosterol and vitamin E isomers in sick cows. In view of the anti-inflammatory property of sitosterol and the anti-oxidative role of α-tocopherol [33, 34], those cows can not effectively utilize them, rendering them susceptible to infectious diseases and oxidative stress. However, it might be that the host responds to sickness by increasing those metabolites. It might be also both ways. Therefore, Need more search to see why they are high.

Although CK and SK had many of the same metabolites in both CK and SK, some different metabolites were also found between the two groups, as is shown in Figure 4b. Specifically, seven different metabolites increased in CK vs. NC, including 2-methyl-3-hydroxybutyric acid (2Me3HB) is a metabolite from isoleucine catabolism, β-oxidation of fatty acids and ketogenesis [25], and xylitol, is a precursor of xylulose 5-phosphate, and an intermediate of the pentose phosphate and glycolytic pathways [35]. The decreased metabolites included 3-hydroxyisovaleric acid (3HIV) from catabolism of the ketogenic amino acid leucine [36], 4-aminobutyric acid (GABA) from L-glutamic acid catabolism [37], melibiose from galactose metabolism [25], erythritol, a precursor of fructose 6-phosphate [38], and l-serine (L-ser), a glucogenic and ketogenic amino acid [39]. The results suggest that these carbohydrates and amino acids may be important factors in the development of CK. In addition, five different metabolites were up-regulated in SK compared to NC, and included 4-hydroxyproline (4-HYP) from proline hydroxylation [40], and l-leucine (L-leu), which is a ketogenic amino acid [19]. Down-regulated metabolites between the two groups included 2-ketoglutaric acid (2-KG), which is involved in the TCA cycle [41], l-ornithine (L-ori), which participates in the urea cycle [42], and methylmalonic acid (MMA), which is from a C-methylated derivative of malonate [43]. These results indicate that SK may disturb the TCA cycle, urea cycle, and odd-numbered fatty acid metabolism. Therefore, a few different metabolites had been found in CK or SK in this study, suggesting that they may play certain roles in identifying development process or types of ketosis.

Finally, Figure 4b shows the six metabolites that decreased in CK compared to SK; namely, GABA; L-ser; 4HYP; citrate; 3,4-dihydroxybutyric acid (DHB); l-proline (L-pro). Citrate is an intermediate in the TCA cycle [44]. DHB may be a product of amino acid catabolism [45]. L-pro is an important component of collagen and is derived from the amino acid l-glutamate [46]. Therefore, a decrease of the above mentioned metabolites, especially citrate, can further aggravate ketosis due to interruption of TCA cycle and lack of amino acids. Our research may shed light on the clinical potential of some potential metabolites in the development of ketosis.

In Table 2, long or medium-chain fatty acids and hydroxyl acids such as C-9-HA, MA, BHBA, and 2PC, etc. increased significantly in SK and CK compared with NC while metabolites such as 3-HV and 3-HIV changed differently in SK and CK group. The decreased metabolites were mainly some compounds such as LA, L-ala, GLCA, GA, Gal, Glc, etc., which are relative to glycolysis and TCA cycle. Thus, those metabolites may become potential biomarkers for diagnosing ketosis including SK and CK in dairy cows.

However, there were less differential metabolites in CK vs SK groups. Long or medium-chain fatty acids and hydroxyl acids such as HA, SA, BHBA, 3-HV ect. increased significantly in CK compared with SK, but amino acids like L-ala, L-ser, L-pro, 4-HYP, and GABA & citrate and DHBA decreased significantly. Therefore, those metabolites may be potential biomarkers for distinguishing the serious degree of ketosis in dairy cows.

In general, many metabolic pathways are disturbed in the development of ketosis. The supply of amino acids for gluconeogenesis decreases because there is a negative energy balance (NEB) in both CK and SK. In the meantime, the huge energy requirement for cow milk production after calving causes increased availability of fatty acids for oxidation, which displaces glucose as the oxidative fuel. In NEB situations, hypoglycemia activates fat mobilization, leading to NEFAs accumulation in liver, production of advanced glycolic end products, and increased oxidative stress, which result in many metabolic dysfunctions. The role of certain intermediates like THBA, DHBA, AABA, and MMA in the development of ketosis is worthy of further exploration in order to understand metabolic disturbances in ketosis.

Diagnostic test and receiver–operator characteristic curves (ROC)

PCA, PLS-DA, and OPLS-DA, the most commonly used algorithms in metabolomics [33, 46], were employed to process the GC/MS data. The OPLS-DA score plot showed that the two groups were scattered into two different regions using PC1 (Figures 3a, 3b, 3c). ROC analysis using the cross-validated predicted Y (predicted class) values was performed to validate the robustness of the OPLS-DA model. Sensitivity and specificity tradeoffs were summarized for each variable in Table 3, such as C-9-HA for CK vs. NC, 2PC for SK vs. NC, and BHBA for CK vs. SK, using the area under the curve (AUC), and were calculated using the trapezoidal rule (Figure 3d, AUC 0.907 in CK vs. NC; Figure 3e, AUC 0.923 in SK vs. NC; Figure 3f, AUC 0.953 in CK vs. SK). It has been reported that ROC analysis is able to determine easily ability for identifying disease at any cutoff. AUC value from ROC is usually between 1.0 and 0.5. AUC is more close to 1, the higher the accuracy test is, and the bigger the diagnostic value is , then the less false positive or negative ratio is [47, 48]. In addition, the bigger positive likelihood ratio (+LR) is, the higher the true positive probability is [48, 49]. According to AUC values and +LR of any two groups in Table 3, 31 metabolites for CK vs NC, 19 for SK vs NC, and 8 for CK vs SK may be potential biomarkers for identifying ketosis. Of course, those metabolites with high AUC values and high fold changes will be better potential diagnostic biomarkers for detection and classification of ketosis in the future.

Table 3 The potential diagnostic biomarkers of any two groups of CK, SK, and NC

In this study, the detected metabolites that were differentially found between CK and SK were closely related to alterations in carbohydrates, fatty acids, amino acids, even sitosterol and vitamin E isomers, etc. Our study confirmed that BHBA is an excellent biomarker for ketosis including CK and SK, which has long been considered a gold indicator for diagnosing ketosis [2, 3]. Hypoglycemia was the main biochemical characteristics in CK and SK (Table 1), which is in accordance with the metabolomics results in Table 2. Furthermore, C-9-HA, which belongs to the NEFA family, was significantly increased in CK compared to NC, an important indicator of the NEB and fat mobilization. However, Glc and the NEFA family could not become potential biomarkers for distingushing CK and SK. In addition, the increased 2PC in SK compared to NC, which has not been previously reported as a potential biomarker, suggested that some new biomarkers may play important roles in the pathogenesis and development of ketosis in dairy cows.


The present study is an integrated analysis in dairy cattle ketosis based on plasma metabolomic profiling by GC/MS. We first discovered that 40 metabolites (i.e. fatty acids, amino acids, carbohydrates, and others) were differentially found among CK, SK and NC. This proved that metabolic patterns of ketosis could be reflected by metabolomics technology, which could, to some extent, reveal the development and progression of ketosis. To research ketosis, potential biomarkers of CK and SK could uncover the same or different modes of metabolites and metabolic pathways in the development and progression of ketosis. Furthermore, new potential metabolites could shed light on new strategies for the diagnosis, prognosis, and prevention of ketosis in the future.


Animals and sample collection

All experimental animals were treated according to the International Guiding Principles for Biomedical Research Involving Animals. Twenty-two cows with clinical ketosis (CK), 32 cows with subclinical ketosis (SK), and 22 normal control cows (NC) from an intensive 1000 dairy cattle farm in Mishan, Heilongjiang, China were used in this exploratory study. All cows were fed a total mixed ration (TMR) during the transition period. TMR consisted of 8.5 kg of concentrated feed, 18.5 kg of silage maize, 4 kg of alfalfa, and 0.35 kg of fat. TMR consisted of 55.60% dry matter (DM), 16% crude protein, 1.75 mcal/DM net energy for lactation (NEL), 5.60% fat, 39.10% neutral detergent fiber (NDF), 20.30% acid detergent fiber (ADF), 180 g Ca, and 116 g P. The cows were considered to have CK if they showed typical clinical symptoms including lack of appetite, apparent wasting, low milk yield, even nervous signs etc., and had high plasma BHBA concentrations ( > 1.60 mmol/L). If the cows had little clinical symptoms and high plasma BHBA concentrations (> 1.20 mmol/L), they were considered to have SK. If the cows had no clinical signs and normal plasma BHBA concentrations (< 1.00 mmol/L), they were considered normal controls (NC) [3, 4]. Within 12 h of calving, whole first-morning blood samples were collected and immediately centrifuged at 1400 × g for 10 min at room temperature. The supernatants were aliquoted into Eppendorf tubes (1 mL plasma/tube) and stored at −80°C until analysis. Both clinical parameters and plasma metabolomics were measured for all experimental animals. Main characteristics and blood biochemical parameters, which included age, parity, milk yield (MY, kg/d), body condition score (BCS), blood glucose (Glc), β-hydroxaybutyric acid (BHBA), triglycerides (TG), aspartate aminotransferase (AST), and cholesterol (CHO), are shown in Table 1.

Sample preparation

The sample preparation was modified according to literature [48, 49] and described as follows. To an Eppendorf tube were added 50 μL of thawed plasma at room temperature, 10 μL 13C6-Leucine (0.5 mg/mL in water, as internal standard), and 180 μL cold methanol. The mixture was vigorously vortexed for 30 s and centrifuged at 16 000×g for 15 min at 4°C, and then 192 μL of the supernatant was transferred to a glass vial for drying under gentle nitrogen stream. A volume of 80 μL of methoxylamine hydrochloride in pyridine (20 mg/mL) was added to the residue, and the solution was strongly vortexed for 30 s and incubated at 37°C for 90 min. The sample was derivatized using 80 μL N,O-bis (trimethylsilyl)-trifluoroacetamide (BSTFA) with 1% trimethylsilyl chloride (TMCS) at 70°C for 60 min. To avoid decomposition, all derivatized samples were analyzed within 48 h by GC/MS (Agilent 7890A GC /5975C MS).

Data preprocessing and multivariate statistical analysis

The extraction, alignment, deconvolution, and further processing of raw GC/MS data were referred to the protocols of literature [4749], except the mass window was set to 70–600 m/z. The preprocessed data table, including observations (sample name), variables (rt_mz), and summarized peak area, was imported into Simca-P 11.0 software (Umetrics AB, Umea, Sweden) for multivariate statistical analysis, where mean-centering and UV-scaling were employed prior to PCA, PLS-DA, and OPLS-DA. Principal component analysis (PCA) was initially performed for visualizing the trends of samples in this study. To remove the noise, partial least square-discriminant analysis (PLS-DA) and orthogonal partial least square-discriminant analysis (OPLS-DA) were utilized to explore differences among groups. PLS-DA was used to model the metabolic profiling data among CK, SK, and NC. OPLS-DA filtered unrelated factors to CK, SK, and NC in the metabolic profiling. Model parameters, R2Y and Q2, were applied to evaluate model stability and the ability to explain and predict the raw data. In addition, the possibility of over-fitting was tested using the permutation test [11, 12, 15].

Identification and qualification of differential metabolites

Differential metabolites (i.e., potential biomarkers) between arbitrary two groups were found by combining Variable Importance in Projection (VIP >1) values of the OPLS-DA model and students t-test (P <0.05). The structural qualification of potential biomarkers was performed according to the method of literature [57]. The automated mass spectral deconvolution and identification system (AMDIS, National Institute of Standards and Technology, Gaithersburg, MD) was first utilized to deconvolute mass peaks and search the self-constructed standard library by matching retention time and mass spectra pattern. The peaks which were not matched from standard library were introduced to the NIST MS 2.0 software for automatically searching from the NIST 08 library. The peaks with matching similarity larger than 80% were assigned as candidate compounds.

Statistical analysis

Variance was used to analyze the some information based on the diagnoses of CK, SK, and NC. All data processed by ANOVA (SPSS for Windows version 11.0) were presented as mean ± standard deviation (SD). P values of less than were considered to indicate a significant difference. In addition, the receiver operating characteristic (ROC) analysis was carried out using SPSS, and the AUC value, specificity, sensitivity were calculated to evaluate the diagnostic value of the potential biomarkers from the differential metabolites of the disease. P values of less than 0.05 were considered to indicate a significant difference [4749].



Aminomalonic acid


α-aminobutyric acid


3-hydroxybutyric acid


cis-9-hexadecenoic acid


3,4-dihydroxybutyric acid


Glycolic acid








Glucuronic acid


4-aminobutyric acid


Heptadecanoic acid


3-hydroxy-3-methylglutaric acid


Lactic acid














Myristic acid


Methylmalonic acid


Palmitic acid


pyroglutamic acid


Stearic acid


THBA: 2,3,4-trihydroxybutyric acid; Trans-9-octadecenoic acid


3-hydroxyvaleric acid




2-ketoglutaric acid


2-piperidinecarboxylic acid




3-hydroxyisovaleric acid


2-methyl-3-hydroxybutyric acid.


  1. Grummer RR: Etiology of lipid-related metabolic disorders in periparturient dairy cows. J Dairy Sci. 1993, 76: 3882-3896. 10.3168/jds.S0022-0302(93)77729-2.

    Article  PubMed  CAS  Google Scholar 

  2. Herdt TH: Ruminant adaptation to negative energy balance. Influences on the etiology of ketosis and fatty liver. Vet Clin North Am Food Anim Pract. 2000, 16: 215-230.

    PubMed  CAS  Google Scholar 

  3. Oetzel GR: Monitoring and testing dairy herds for metabolic diseases. Vet Clin North Amer Food Anim Pract. 2004, 20: 651-674. 10.1016/j.cvfa.2004.06.006.

    Article  Google Scholar 

  4. Xia C, Wang Z, Xu C, Zhang HY: Concentrations of plasma metabolites, hormones, and mRNA abundance of adipose leptin and hormone-sensitive lipase in ketotic and nonketotic dairy cows. J Vet Intern Med. 2012, 26: 415-417. 10.1111/j.1939-1676.2011.00863.x.

    Article  PubMed  CAS  Google Scholar 

  5. Ilves A, Harzia H, Ling K, Ots M, Soomets U, Kilk K: Alterations in milk and blood metabolomes during the first months of lactation in dairy cows. J Dairy Sci. 2012, 95: 5788-5797. 10.3168/jds.2012-5617.

    Article  PubMed  CAS  Google Scholar 

  6. Dunn WB, Bailey NJ, Johnson HE: Measuring the metabolome: current analytical technologies. Analyst. 2005, 130: 606-625. 10.1039/b418288j.

    Article  PubMed  CAS  Google Scholar 

  7. Klein MS, Buttchereit N, Miemczyk SP, Immervoll AK, Louis C, Wiedemann S, Junge W, Thaller G, Oefner PJ, Gronwald W: NMR metabolomic analysis of dairy cows reveals milk glycerophosphocholine to phosphocholine ratio as prognostic biomarker for risk of ketosis. J Proteome Res. 2012, 11: 1373-1381. 10.1021/pr201017n.

    Article  PubMed  CAS  Google Scholar 

  8. Burim NA, Qendrim Z, Fozia S, Psychogios N, Michael JL, Dunn SM, Jianguo X, Wishart DS: Metabolomics reveals unhealthy alterations in rumen metabolism with increased proportion of cereal grain in the diet of dairy cows. Metabolomics. 2010, 6 (4): 583-594. 10.1007/s11306-010-0227-6.

    Article  Google Scholar 

  9. Saleem F, Ametaj BN, Bouatra S, Mandal R, Zebeli Q, Dunn SM, Wishart DS: A metabolomics approach to uncover the effects of grain diets on rumen health in dairy cows. J Dairy Sci. 2012, 95 (11): 6606-6623. 10.3168/jds.2012-5403.

    Article  PubMed  CAS  Google Scholar 

  10. Saleem F, Bouatra S, Chi-Guo A, Psychogios N, Mandal R, Dunn SM, Burim NA, David SW: he bovine ruminal fluid metabolome. Metabolomics. 2013, 9 (2): 360-378. 10.1007/s11306-012-0458-9.

    Article  CAS  Google Scholar 

  11. Eriksson L, Johansson E, Kettaneh WN, Trygg J, Wikström C, Wold S: Multi- and megavariate data analysis—principles and applications. Umea: UmetricsAB: 2001.

    Google Scholar 

  12. Nicholson JK, Connelly J, Lindon JC, Holmes E: Metabolomics: a platform for studying drug toxicity and gene function. Nat Rev Drug Discov. 2002, 1: 153-10.1038/nrd728.

    Article  PubMed  CAS  Google Scholar 

  13. Rajala-Schultz PJ, Gröhn YT, McCulloch CE: Effects of milk fever, ketosis, and lameness on milk yield in dairy cows. J Dairy Sci. 1999, 82 (2): 288-294. 10.3168/jds.S0022-0302(99)75235-5.

    Article  PubMed  CAS  Google Scholar 

  14. Gillund P, Reksen O, Gröhn YT, Karlberg K: Body condition related to ketosis and reproductive performance in Norwegian dairy cows. J Dairy Sci. 2001, 84 (6): 1390-1396. 10.3168/jds.S0022-0302(01)70170-1.

    Article  PubMed  CAS  Google Scholar 

  15. Lindgren F, Hansen B, Karcher W, Eriksson L: Model validation by permutation tests: Applications to variable selection. J Chemom. 1996, 10: 521-532. 10.1002/(SICI)1099-128X(199609)10:5/6<521::AID-CEM448>3.0.CO;2-J.

    Article  CAS  Google Scholar 

  16. Wu H, Xue RY, Dong L, Liu TT, Deng CH, Zeng HZ, Shen XZ: Metabolomic profiling of human urine in hepatocellular carcinoma patients using gas chromatography/mass spectrometry. Anal Chim Acta. 2009, 648: 98-104. 10.1016/j.aca.2009.06.033.

    Article  PubMed  CAS  Google Scholar 

  17. Geenen S, Guallar-Hoyas C, Michopoulos F, Kenna JG, Kolaja KL, Westerhoff HV, Thomas P, Wilson ID: HPLC-MS/MS methods for the quantitative analysis of 5-oxoproline (pyroglutamate) in rat plasma and hepatic cell line culture medium. J Pharm Biomed Anal. 2011, 56: 655-663. 10.1016/j.jpba.2011.06.001.

    Article  PubMed  CAS  Google Scholar 

  18. Bogdanov M, Matson WR, Wang L, Matson T, Pullman RS, Bressman SS, Beal MF: Metabolomic profiling to develop blood biomarkers for Parkinson's disease. Brain. 2008, 131: 389-396. 10.1093/brain/awm304.

    Article  PubMed  Google Scholar 

  19. Nelson DL, Cox MM: Lehninger's Principles of Biochemistry. 4th edition. New York: WH. Freeman and Company; 2005.

    Google Scholar 

  20. Ametaj BN, Emmanuel DG, Zebeli Q, Dunn SM: Feeding high proportions of barley grain in a total mixed ration perturbs diurnal patterns of plasma metabolites in lactating dairy cows. J Dairy Sci. 2009, 92: 1084-1091. 10.3168/jds.2008-1465.

    Article  PubMed  CAS  Google Scholar 

  21. King C, Rios G, Green M, Tephly T: UDP-glucuronosyltransferases. Curr Drug Metab. 2000, 1: 143-161. 10.2174/1389200003339171.

    Article  PubMed  CAS  Google Scholar 

  22. Lankinen M, Schwab U, Seppänen-Laakso T, Mattila I, Juntunen K, Mykkänen H: Poutanen K, Gylling H, Oresic M: Metabolomic analysis of plasma metabolites that may mediate effects of rye bread on satiety and weight maintenance in postmenopausal women. J Nutr. 2011, 141: 31-36. 10.3945/jn.110.131656.

    Article  PubMed  CAS  Google Scholar 

  23. Bosch AM: Classical galactosaemia revisited. J Inherit Metab Dis. 2006, 29: 516-525. 10.1007/s10545-006-0382-0.

    Article  PubMed  CAS  Google Scholar 

  24. Huzzey JM, Nydam DV, Grant RJ, Overton TR: The effects of overstocking Holstein dairy cattle during the dry period on cortisol secretion and energy metabolism. J Dairy Sci. 2012, 95: 4421-4433. 10.3168/jds.2011-5037.

    Article  PubMed  CAS  Google Scholar 

  25. Salway JG: Metabolism at a glance. 3rd edition. 350 Main Street, Malden, Massachusetts, 02148-5020, USA: Blackwell Pub Ltd.; 2004.

    Google Scholar 

  26. Doepel L, Lapierre H, Kenneky JJ: Peripartum performance and metabolism of dairy cows in response to prepartum energy and protein intake. J Dairy Sci. 2002, 85: 2315-2334. 10.3168/jds.S0022-0302(02)74312-9.

    Article  PubMed  CAS  Google Scholar 

  27. Marsh DC, Vreugdenhil PK, Mack VE, Belzer FO, Southard JH: Glycine protects hepatocytes from injury caused by anoxia, cold ischemia and mitochondrial inhibitors, but not injury caused by calcium ionophores or oxidative stress. Hepatology. 1993, 17: 91-98.

    Article  PubMed  CAS  Google Scholar 

  28. Marvin DA, Francis B: A Metabolic Study of α-aminobutyric acid. J Biol Chem. 1949, 180: 1059-1063. Downloaded from

    Google Scholar 

  29. Hinko CN, Crider AM, Kliem MA, Steinmiller CL, Seo TH, Bin H, Venkatarangan P, EL-assadi AA, Hyejung C, Burns CM, Tietz EI, Andersen PH, Klitgaard H: Anticonvulsant activity of novel derivatives of 2- and 3-piperidinecarboxylic acid in mice and rats. Neuropharm. 1996, 35: 1721-1735. 10.1016/S0028-3908(96)00105-0.

    Article  CAS  Google Scholar 

  30. Renée PK, Takhar K, Kathryn AJ, Katherine RT, Adams J, Lisa NB, Gerd K, Brewer WY, Charles RR, Henri B: Parenteral and enteral metabolism of anaplerotic triheptanoin in normal rats. Am J Physiol Endocrinol Metab. 2006, 291: E860-E866. 10.1152/ajpendo.00366.2005.

    Article  Google Scholar 

  31. Mao YY, Bai JQ, Chen JH, Shou ZF, Qiang H, Wu JY, Chen Y, Cheng YY: A pilot study of GC/MS-based serum metabolic profiling of acute rejection in renal transplantation. Transplant Immunology. 2008, 19: 74-80. 10.1016/j.trim.2008.01.006.

    Article  PubMed  CAS  Google Scholar 

  32. Effros RM: Alpha aminobutyric acid, an alternative measure of hepatic injury in sepsis?. Transl Res. 2011, 158: 326-327. 10.1016/j.trsl.2011.07.003.

    Article  PubMed  CAS  Google Scholar 

  33. Ros MM, Sterk SS, Verhagen H, Stalenhoef AF, de Jong N: Phytosterol consumption and the anabolic steroid boldenone in humans: a hypothesis piloted. Food Addit Contam. 2007, 24: 679-684. 10.1080/02652030701216727.

    Article  PubMed  CAS  Google Scholar 

  34. Wefers HS: The protection of ascorbate and glutathione against microsomal lipid peroxidation is dependent on Vitamin E. Europ J Biochem. 1988, 174: 353-357. 10.1111/j.1432-1033.1988.tb14105.x.

    Article  PubMed  CAS  Google Scholar 

  35. Dupriez VJ, Rousseau GG: Glucose response elements in a gene that codes for 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase. DNA Cell Biol. 1997, 16: 1075-1085. 10.1089/dna.1997.16.1075.

    Article  PubMed  CAS  Google Scholar 

  36. Ferreira G, Weiss WP, Willett LB: Changes in measures of biotin status do not reflect milk yield responses when dairy cows are fed supplemental biotin. J Dairy Sci. 2007, 90: 1452-1459. 10.3168/jds.S0022-0302(07)71630-2.

    Article  PubMed  CAS  Google Scholar 

  37. Watanabe M, Maemura K, Kanbara K, Tamayama T, Hayasaki H: GABA and GABA receptors in the central nervous system and other organs. Int Rev Cytol. 2002, 213: 1-47.

    Article  PubMed  CAS  Google Scholar 

  38. Munro IC, Berndt WO, Borzelleca JF, Flamm G, Lynch BS, Kennepohl E, Bär EA, Modderman J: Erythritol: an interpretive summary of biochemical, metabolic, toxicological and clinical data. Food Chem Toxicol. 1998, 36: 1139-1174. 10.1016/S0278-6915(98)00091-X.

    Article  PubMed  CAS  Google Scholar 

  39. Hamana M, Ohtsuka H, Oikawa M, Kawamura S: Blood free amino acids in the postpartum dairy cattle with left displaced abomasums. J Vet Med Sci. 2010, 72: 1355-1358. 10.1292/jvms.09-0239.

    Article  PubMed  CAS  Google Scholar 

  40. Siddiqi NJ, Alhomida AS, Pandey VC: Hydroxyproline distribution in the plasma of various mammals. J Biochem Mol Biol Biophys. 2002, 6: 159-163. 10.1080/10258140290027324.

    Article  PubMed  CAS  Google Scholar 

  41. Postic C, Dentin R, Girard J: Role of the liver in the control of carbohydrate and lipid homeostasis. Diabetes Metab. 2004, 30: 398-408. 10.1016/S1262-3636(07)70133-7.

    Article  PubMed  CAS  Google Scholar 

  42. Motyl T, Barej W: Plasma amino acid indices and urinary 3-methyl histidine excretion in dairy cows in early lactation. Ann Rech Vet. 1986, 17: 153-157.

    PubMed  CAS  Google Scholar 

  43. James W, Croom JR, Bauman DE, Davis CL: Methylmalonic acid in low-fat milk syndrome. J Dairy Sci. 1981, 64: 649-654. 10.3168/jds.S0022-0302(81)82625-2.

    Article  Google Scholar 

  44. Yamashina S, Ikejima K, Enomoto N, Takei Y, Sato N: Glycine as a therapeutic immuno- nutrient for alcoholic liver disease alcohol. Clin Exp Res. 2005, 29: 162S-165S. 10.1097/01.alc.0000189281.82523.6c.

    Article  CAS  Google Scholar 

  45. Joanna KC: Noninvasive urinary organic acids test to assess biochemical and nutritional individuality in autistic children. Clinical Biochem. 2011, 44: 686-691. 10.1016/j.clinbiochem.2011.01.015.

    Article  Google Scholar 

  46. Kim D, Yamagishi N, Ueki A, Miura M, Saito F, Sato S, Furuhama K: Changes in plasma bone metabolic markers in periparturient dairy cows. J Vet Med Sci. 2010, 72: 773-776. 10.1292/jvms.09-0409.

    Article  PubMed  CAS  Google Scholar 

  47. Nishiumi S, Kobayashi T, Ikeda A, Yoshie T, Kibi M, Izumi Y, Okuno T, Hayashi N, Kawano S, Takenawa T, Azuma T, Yoshida M: A novel serum metabolomics-based diagnostic approach for colorectal cancer. PloS ONE. 2012, 7 (7): e40459-10.1371/journal.pone.0040459.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  48. Chan EC, Koh PK, Mal M, Cheah PY, Eu KW, Backshall A, Cavill R, Nicholson JK, Keun HC: Metabolic profiling of human colorectal cancer using high-resolution magic angle spinning nuclear magnetic resonance (HR-MAS NMR) spectroscopy and gas chromatography mass spectrometry (GC/MS). J Proteome Res. 2009, 8: 352-361. 10.1021/pr8006232.

    Article  PubMed  CAS  Google Scholar 

  49. Gao XF, Pujos-Guillot E, Sébédio JL: Development of a quantitative metabolomic approach to study clinical human fecal water metabolome based on trimethylsilylation derivatization and GC/MS analysis. Anal Chem. 2010, 82: 6447-6456. 10.1021/ac1006552.

    Article  PubMed  CAS  Google Scholar 

Download references


The authors thank the Animal Medicine Key Laboratory, at Heilongjiang Bayi Agricultural University for instrumental help with biochemical test and Shanghai Sensichip Infotech Co., Ltd (Shanghai, 200433, China) for the GS/MS work. This study was funded by National Science and Technology Program of China (2012BAD12B05-2) and National Nature Science Foundation of China (31072181) & the Foundation of Heilongjiang Province (No. GZ02B033). The authors declare that they have nothing to disclose regarding funding from industry or conflict of interest with respect to this manuscript.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Cheng Xia.

Additional information

Competing interests

There are no conflicts of interests.

Authors’ contributions

HYZ performed the metabolomic samples preparation, data generation, analyzed and wrote the draft of the manuscript. CX participated in design of this study and helped to interpret the metabolomic data and analysis. CX and LW conceived of the study, participated in its design and coordination, helped to interpret the statistical analysis of the metabolomic data and draft the manuscript. LWS and SS did sample collection, sample preparation, facilitated sample collection, contributed to design of the study, provided knowledge in dairy metabolic biology and helped to draft the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Zhang, H., Wu, L., Xu, C. et al. Plasma metabolomic profiling of dairy cows affected with ketosis using gas chromatography/mass spectrometry. BMC Vet Res 9, 186 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: