Metabolic Clusters of Early-Lactating Dairy Cows Based on Blood β-hydroxybutyrate Trajectories and Predicted from Milk Compounds

High-yielding dairy cows encounter metabolic challenges in early lactation. Typically, β-hydroxybutyrate (BHB), measured at a specific time point is employed to diagnose the metabolic status of cows based on a predetermined threshold. However, in early lactation, BHB is highly dynamic, and there is high interindividual variability in its time profile. This could limit the effectiveness of the single measurement and threshold-based diagnosis probably contributing to the disparities in reports linking metabolic status with productive and reproductive outcomes. This research delves into the examination of the trajectories of BHB to unveil inter-cow variations and identify latent metabolic groups. We compiled a data set from 2 observational studies involving a total of 195 lactations from multiparous Holstein Friesian cows. The data set encompasses measurements of BHB, NEFA, and insulin from blood samples collected at 3, 6, 9, and 21 d in milk (DIM), along with weekly determinations of milk composition and fatty acids (FA) proportions in milk fat. In both experiments, milk yield (MY) and feed intake were recorded daily during the first month of lactation. We explored interindividual and intraindividual variations in metabolic responses using the trajectories of blood BHB and evaluated the presence of distinct metabolic groups based on such variations. For this purpose, we employed the growth mixture model (GMM), a trajectory clustering technique. Our findings unveil novel insights into the diverse metabolic responses among cows, encompassing both trajectory patterns and the magnitude of blood BHB concentrations. Specifically, we identified 3 latent metabolic groups: the "QuiBHB" cluster (≈10%) exhibited a higher initial BHB concentration than other clusters, peaking on d 9 (average maximum BHB of 2.4 mM) and then declining by d 21; the "SloBHB" cluster (≈23%) started with a lower BHB concentration, gradually increasing until d 9, and at the highest BHB concentration at d 21 (1.6 mM serum BHB at the end of the experimental period); and the "LoBHB" cluster (≈67%) began with the lowest serum BHB concentration (serum BHB <0.75 mM), remaining relatively stable throughout the sampling period. Notably, the 3 metabolic groups exhibited significant physiological disparities, evident in blood NEFA and insulin concentrations. The QuiBHB and SloBHB cows exhibited higher NEFA and lower insulin concentrations as compared with the LoBHB cows. Interestingly, these metabolic differences extended to MY and DMI during the first month of lactation. The elevated BHB concentrations observed in QuiBHB cows were linked with lower DMI and MY as compared with SloBHB and LoBHB cows. Accordingly, these animals were considered metabolically impaired. Conversely, SloBHB cows displayed higher MY along with increased DMI, and thus the elevated BHB might be indicative of an adaptive response for these cows. The QuiBHB cows also displayed higher proportions of unsaturated FA (UFA), monounsaturated FA (MUFA), and total C18:1 FA in milk during the first week of lactation. Prediction of the QuiBHB cows using these FA and test day variables resulted in moderate predictive accuracy (ROCAUC > 0.7). Given the limited sample size for the development of prediction models, and the variation in DIM among samples in the same week, the result is indicative of the predictive potential of the model and room for model optimization. In summary, distinct metabolic groups of cows could be identified based on the trajectories of blood BHB in early lactation.