Discovery of different metabotypes in overconditioned dairy cows by means of machine learning

Abstract

Using data from targeted metabolomics in serum in combination with machine learning (ML) approaches, we aimed at (1) identifying divergent metabotypes in overconditioned cows and at (2) exploring how metabotypes are associated with lactation performance, blood metabolites, and hormones. In a previously established animal model, 38 pregnant multiparous Holstein cows were assigned to 2 groups that were fed differently to reach either high (HBCS) or normal (NBCS) body condition score (BCS) and backfat thickness (BFT) until dryoff at -49 d before calving [NBCS: BCS < 3.5 (3.02 ± 0.24) and BFT < 1.2 cm (0.92 ± 0.21), mean ± SD; HBCS: BCS > 3.75 (3.82 ± 0.33) and BFT > 1.4 cm (2.36 ± 0.35)]. Cows were then fed the same diets during the dry period and the subsequent lactation, and maintained the differences in BFT and BCS throughout the study. Blood samples were collected weekly from 7 wk antepartum (ap) to 12 wk postpartum (pp) to assess serum concentrations of metabolites (by targeted metabolomics and by classical analyses) and metabolic hormones. Metabolic clustering by applying 4 supervised ML-based classifiers [sequential minimal optimization (SMO), random forest (RF), alternating decision tree (ADTree), and naïve Bayes-updatable (NB)] on the changes (d 21 pp minus d 49 ap) in concentrations of 170 serum metabolites resulted in 4 distinct metabolic clusters: HBCS predicted HBCS (HBCS-PH, n = 13), HBCS predicted NBCS (HBCS-PN, n = 6), NBCS predicted NBCS (NBCS-PN, n = 15), and NBCS predicted HBCS (NBCS-PH, n = 4). The accuracies of SMO, RF, ADTree, and NB classifiers were >70%. Because the number of NBCS-PH cows was low, we did not consider this group for further comparisons. Dry matter intake (kg/d and percentage of body weight) and energy intake were greater in HBCS-PN than in HBCS-PH in early lactation, and HBCS-PN also reached a positive energy balance earlier than did HBCS-PH. Milk yield was not different between groups, but milk protein percentage was greater in HBCS-PN than in HBCS-PH cows. The circulating concentrations of fatty acids (FA) increased during early lactation in both groups, but HBCS-PN cows had lower concentrations of β-hydroxybutyrate, indicating lower ketogenesis compared with HBCS-PH cows. The concentrations of insulin, insulin-like growth factor 1, leptin, adiponectin, haptoglobin, glucose, and revised quantitative insulin sensitivity check index did not differ between the groups, whereas serum concentrations of glycerophospholipids were lower before calving in HBCS-PH than in HBCS-PN cows. Glycine was the only amino acid that had higher concentration after calving in HBCS-PH than in HBCS-PN cows. The circulating concentrations of some short- (C2, C3, and C4) and long-chain (C12, C16:0, C18:0, and C18:1) acylcarnitines on d 21 pp were greater in HBCS-PH than in HBCS-PN cows, indicating incomplete FA oxidation. In conclusion, the use of ML approaches involving data from targeted metabolomics in serum is a promising method for differentiating divergent metabotypes from apparently similar BCS phenotypes. Further investigations, using larger numbers of cows and farms, are warranted for confirmation of this finding.