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

Background Ketosis is an important problem for dairy cows` production performance. However, it is still little known about plasma metabolomics details of dairy ketosis. Results 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. Conclusion 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.


Background
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 [3][4][5]. 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 [6][7][8][9][10]. 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].
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 [1][2][3]. 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. 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.

GC/MS metabolomic profiling
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  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, R 2 Y and Q 2 , were 0.828 and 0.693, respectively. R 2 Y showed the explanative ability of the model, and Q 2 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.
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, R 2 Y and Q 2 , 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 R 2 and Q 2 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. 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   Significant changes in the levels of corresponding metabolites between arbitrary two groups by the ANOVA and t-test (P < 0.05). The p-value calculated from two-tailed student's t test.

Metabolic disturbances of ketosis
f Fold change (FC) was calculated as binary logarithm of average mass response ratio between two groups, where the positive value means that the average mass response of the metabolite in the former is larger than that in the latter and vice versa. g Information from KEGG database or NuGOwiki database.
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 [17][18][19][20]. 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 [21][22][23][24]. 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 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-9octadecenoic acid (T-9-OA), myristic acid (MA), cis-9hexadecenoic 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 Lglutamic acid catabolism [37], melibiose from galactose metabolism [25], erythritol, a precursor of fructose 6phosphate [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.
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.

Conclusion
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 13 C 6 -Leucine (0.5 mg/mL in water, as internal standard), and 180 μL cold methanol.

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 [47][48][49], 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 squarediscriminant 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, R 2 Y and Q 2 , 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 [5][6][7]. 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 [47][48][49].