Metabolic disorders in early lactation have negative effects on dairy cow health and farm profitability. One method for monitoring the metabolic status of cows is metabolic profiling, which uses associations between the concentrations of several metabolites in serum and the presence of metabolic disorders. In this cross-sectional study, we investigated the use of midinfrared (MIR) spectroscopy of milk for predicting the concentrations of these metabolites in serum. Between July and October 2017, serum samples were taken from 773 early-lactation Holstein Friesian cows located on 4 farms in the Gippsland region of southeastern Victoria, Australia, on the same day as milk recording. The concentrations in sera of β-hydroxybutyrate (BHB), fatty acids, urea, Ca, Mg, albumin, and globulins were measured by a commercial diagnostic laboratory. Optimal concentration ranges for each of the 7 metabolites were obtained from the literature. Animals were classified as being either affected or unaffected with metabolic disturbances based on these ranges. Milk samples were analyzed by MIR spectroscopy. The relationships between serum metabolite concentrations and MIR spectra were investigated using partial least squares regression. Partial least squares discriminant analyses (PLS-DA) were used to classify animals as being affected or not affected with metabolic disorders. Calibration equations were constructed using data from a randomly selected subset of cows (n = 579). Data from the remaining cows (n = 194) were used for validation. The coefficient of determination (R 2) of serum BHB, fatty acids, and urea predictions were 0.48, 0.61, and 0.90, respectively. Predictions of Ca, Mg, albumin, and globulin concentrations were poor (0.06 ≤ R 2 ≤ 0.17). The PLS-DA models could predict elevated fatty acid and urea concentrations with an accuracy of approximately 77 and 94%, respectively. A second independent validation data set was assembled in March 2018, comprising blood and milk samples taken from 105 autumn-calving cows of various breeds. The accuracies of BHB and fatty acid predictions were similar to those obtained using the first validation data set. The PLS-DA results were difficult to interpret due to the low prevalence of metabolic disorders in the data set. Our results demonstrate that MIR spectroscopy of milk shows promise for predicting the concentration of BHB, fatty acids, and urea in serum; however, more data are needed to improve prediction accuracies.

The objective of this study was to investigate the potential of milk mid-infrared (MIR) spectroscopy, MIR-derived traits including milk composition, milk fatty acids, and blood metabolic profiles (fatty acids, β-hydroxybutyrate, and urea), and other on-farm data for discriminating cows of good versus poor likelihood of conception to first insemination (i.e., pregnant vs. open). A total of 6,488 spectral and milk production records of 2,987 cows from 19 commercial dairy herds across 3 Australian states were used. Seven models, comprising different explanatory variables, were examined. Model 1 included milk production; concentrations of fat, protein, and lactose; somatic cell count; age at calving; days in milk at herd test; and days from calving to insemination. Model 2 included, in addition to the variables in model 1, milk fatty acids and blood metabolic profiles. The MIR spectrum collected before first insemination was added to model 2 to form model 3. Fat, protein, and lactose percentages, milk fatty acids, and blood metabolic profiles were removed from model 3 to create model 4. Model 5 and model 6 comprised model 4 and either fertility genomic estimated breeding value or principal components obtained from a genomic relationship matrix derived using animal genotypes, respectively. In model 7, all previously described sources of information, but not MIR-derived traits, were used. The models were developed using partial least squares discriminant analysis. The performance of each model was evaluated in 2 ways: 10-fold random cross-validation and herd-by-herd external validation. The accuracy measures were sensitivity (i.e., the proportion of pregnant cows that were correctly classified), specificity (i.e., the proportion of open cows that were correctly classified), and area under the curve (AUC) for the receiver operating curve. The results showed that in all models, prediction accuracy obtained through 10fold random cross-validation was higher than that of herd-by-herd external validation, with the difference in AUC ranging between 0.01 and 0.09. In the herd-byherd external validation, using basic on-farm information (model 1) was not sufficient to classify good-and poor-fertility cows; the sensitivity, specificity, and AUC were around 0.66. Compared with model 1, adding milk fatty acids and blood metabolic profiles (model 2) increased the sensitivity, specificity, and AUC by 0.01, 0.02, and 0.02 unit, respectively (i.e., 0.65, 0.63, and 0.678). Incorporating MIR spectra into model 2 resulted in sensitivity, specificity, and AUC values of 0.73, 0.63, and 0.72, respectively (model 3). The comparable prediction accuracies observed for models 3 and 4 mean that useful information from MIR-derived traits is already included in the spectra. Adding the fertility genomic estimated breeding value and animal genotypes (model 7) produced the highest prediction accuracy, with sensitivity, specificity, and AUC values of 0.75, 0.66, and 0.75, respectively. However, removing either the fertility estimated breeding value or ...

Increased concentrations of some serum biomarkers are known to be associated with impaired health of dairy cows. Therefore, being able to predict these biomarkers, especially in the early stage of lactation, would enable preventive management decision. Some health biomarkers may also be used as phenotypes for genetic improvement for improved animal health. In this study, we validated the accuracy and robustness of models for predicting serum concentrations of β-hydroxybutyrate (BHB), fatty acids, and urea nitrogen, using milk midinfrared (MIR) spectroscopy. The data included 3,262 blood samples of 3,027 lactating Holstein-Friesian cows from 19 dairy herds in Southeastern Australia, collected in the period from July 2017 to April 2020. The models were developed using partial least squares regression and were validated using 10-fold random cross-validation, herd-year by herd-year external validation, and year by year validation. The coefficients of determination (R 2 ) for prediction of serum BHB, fatty acids, and urea obtained through random cross-validation were 0.60, 0.42, and 0.87, respectively. For the herd-year by herd-year external validation, the prediction accuracies held up comparatively well, with R 2 values of 0.49, 0.33, and 0.67 for of serum BHB, fatty acids, and urea, respectively. When the models were developed using data from a single year to predict data collected in future years, the R 2 remained comparable, however, the root mean squared errors increased substantially (4-10 times larger than compared with that of herd-year by herd-year external validation) which could be due to machine differences in spectral response, the change in spectral response of individual machines over time, or other differences associated with farm management between seasons. In conclusion, the mid-infrared equations for predicting serum BHB, fatty acids, and urea have been validated. The prediction equations could be used to help farmers detect cows with metabolic disorders in early lactation in addition to generating novel phenotypes for genetic improvement purposes.

respectively), and moderate for albumin, globulins, and albumin-to-globulin ratio (0.27 ± 0.06, 0.46 ± 0.06, and 0.41 ± 0.06, respectively). The heritability of haptoglobin concentration was close to 0. The magnitude of genetic correlations between traits (estimated using bivariate models) varied considerably (0.01 to 0.96), and standard errors of these correlations were high (0.02 to 0.44). Interestingly, the direction of most genetic correlations was favorable, suggesting that selecting for more optimal concentrations of one biomarker may result in more optimal concentrations of other biomarkers. Correlations between biomarker GEBV and existing breeding values for survival, somatic cell count, and daughter fertility were small to moderate (0.07 to 0.45) and favorable, whereas correlations with breeding values for milk production traits were small (≤0.15). Accuracies of GEBV were evaluated by using 5-fold cross validation, and by calculating accuracies from prediction error variances associated with the GEBV. Accuracies of GEBV predicted using 5-fold cross validation were low (0.05 to 0.27), whereas the means of individual accuracies were greater, ranging from 0.31 to 0.51. Although increasing the size of the reference population should theoretically improve accuracies, our results suggest that genomic prediction of health biomarkers may allow identification of cows that are less susceptible to diseases in early lactation.

Breeding objectives in the dairy industry have shifted from being solely focused on production to including fertility, animal health, and environmental impact. Increased serum concentrations of candidate biomarkers of health and fertility, such as β-hydroxybutyric acid (BHB), fatty acids, and urea are difficult and costly to measure, and thus limit the number of records. Accurate genomic prediction requires a large reference population. The inclusion of milk mid-infrared (MIR) spectroscopic predictions of biomarkers may increase genomic prediction accuracy of these traits. Our objectives were to (1) estimate the heritability of, and genetic correlations between, selected serum biomarkers and their respective MIR predictions, and (2) evaluate genomic prediction accuracies of either only measured serum traits, or serum traits plus MIR-predicted traits. The MIR-predicted traits were either fitted in a single trait model, assuming the measured trait and predicted trait were the same trait, or in a multitrait model, where measured and predicted trait were assumed to be correlated traits. We performed all analyses using relationship matrices constructed from pedigree (A matrix), genotypes (G matrix), or both pedigree and genotypes (H matrix). Our data set comprised up to 2,198 and 9,657 Holstein cows with records for serum biomarkers and MIR-predicted traits, respectively. Heritabilities of measured serum traits ranged from 0.04 to 0.07 for BHB, from 0.13 to 0.21 for fatty acids, and from 0.10 to 0.12 for urea. Heritabilities for MIR-predicted traits were not significantly different from those for the measured traits. Genetic correlations between measured traits and MIR-predicted traits were close to 1 for urea. For BHB and fatty acids, genetic correlations were lower and had large standard errors. The inclusion of MIR predicted urea substantially increased prediction accuracy for urea. For BHB, including MIR-predicted BHB reduced the genomic prediction accuracy, whereas for fatty acids, prediction accuracies were similar with either measured fatty acids, MIR-predicted fatty acids, or both. The high genetic correlation between urea and MIR-predicted urea, in combination with the increased prediction accuracy, demonstrated the potential of using MIR-predicted urea for genomic prediction of urea. For BHB and fatty acids, further studies with larger data sets are required to obtain more accurate estimates of genetic correlations.

Most livestock metabolomic studies involve relatively small, homogenous populations of animals. However, livestock farming systems are non-homogenous, and large and more diverse datasets are required to ensure that biomarkers are robust. The aims of this study were therefore to (1) investigate the feasibility of using a large and diverse dataset for untargeted proton nuclear magnetic resonance (1H NMR) serum metabolomic profiling, and (2) investigate the impact of fixed effects (farm of origin, parity and stage of lactation) on the serum metabolome of early-lactation dairy cows. First, we used multiple linear regression to correct a large spectral dataset (707 cows from 13 farms) for fixed effects prior to multivariate statistical analysis with principal component analysis (PCA). Results showed that farm of origin accounted for up to 57% of overall spectral variation, and nearly 80% of variation for some individual metabolite concentrations. Parity and week of lactation had much smaller effects on both the spectra as a whole and individual metabolites (<3% and <20%, respectively). In order to assess the effect of fixed effects on prediction accuracy and biomarker discovery, we used orthogonal partial least squares (OPLS) regression to quantify the relationship between NMR spectra and concentrations of the current gold standard serum biomarker of energy balance, β-hydroxybutyrate (BHBA). Models constructed using data from multiple farms provided reasonably robust predictions of serum BHBA concentration (0.05 ≤ RMSE ≤ 0.18). Fixed effects influenced the results biomarker discovery; however, these impacts could be controlled using the proposed method of linear regression spectral correction.

The objective of this study was to evaluate the ability of milk infrared spectra to predict blood β-hydroxybutyrate (BHB) concentration for use as a management tool for cow metabolic health on pasturegrazed dairy farms and for large-scale phenotyping for genetic evaluation purposes. The study involved 542 cows (Holstein-Friesian and Holstein-Friesian × Jersey crossbreds), from 2 farms located in the Waikato and Taranaki regions of New Zealand that operated under a seasonal-calving, pasture-based dairy system. Milk infrared spectra were collected once a week during the first 5 wk of lactation. A blood "prick" sample was taken from the ventral labial vein of each cow 3 times a week for the first 5 wk of lactation. The content of BHB in blood was measured immediately using a handheld device. After outlier elimination, 1,910 spectra records and corresponding BHB measures were used for prediction model development. Partial least square regression and partial least squares discriminant analysis were used to develop prediction models for quantitative determination of blood BHB content and for identifying cows with hyperketonemia (HYK). Both quantitative and discriminant predictions were developed using the phenotypes and infrared spectra from two-thirds of the cows (randomly assigned to the calibration set) and tested using the remaining one-third (validation set). A moderate accuracy was obtained for prediction of blood BHB. The coefficient of determination (R 2 ) of the prediction model in calibration was 0.56, with a root mean squared error of prediction of 0.28 mmol/L and a ratio of performance to deviation, calculated as the ratio of the standard deviation of the partial least squares model calibration set to the standard error of predic-tion, of 1.50. In the validation set, the R 2 was 0.50, with root mean squared error of prediction values of 0.32 mmol/L, which resulted in a ratio of performance to deviation of 1.39. When the reference test for HYK was defined as blood concentration of BHB ≥1.2 mmol/L, discriminant models indicated that milk infrared spectra correctly classified 76% of the HYK-positive cows and 82% of the HYK-negative cows. The quantitative models were not able to provide accurate estimates, but they could differentiate between high and low BHB concentrations. Furthermore, the discriminant models allowed the classification of cows with reasonable accuracy. This study indicates that the prediction of blood BHB content or occurrence of HYK from milk spectra is possible with moderate accuracy in pasture-grazed cows and could be used during routine milk testing. Applicability of infrared spectroscopy is not likely suited for obtaining accurate BHB measurements at an individual cow level, but discriminant models might be used in the future as herd-level management tools for classification of cows that are at risk of HYK, whereas quantitative models might provide large-scale phenotypes to be used as an indicator trait for breeding cows with improved metabolic health.

Disorders of energy metabolism, which can result from a failure to adapt to the period of negative energy balance immediately after calving, have significant negative effects on the health, welfare and profitability of dairy cows. The most common biomarkers of energy balance in dairy cows are β-hydroxybutyrate (BHBA) and non-esterified fatty acids (NEFA). While elevated concentrations of these biomarkers are associated with similar negative health and production outcomes, the phenotypic and genetic correlations between them are weak. In this study, we used an untargeted 1H NMR metabolomics approach to investigate the serum metabolomic fingerprints of BHBA and NEFA. Serum samples were collected from 298 cows in early lactation (calibration dataset N = 248, validation N = 50). Metabolomic fingerprinting was done by regressing 1H NMR spectra against BHBA and NEFA concentrations (determined using colorimetric assays) using orthogonal partial least squares regression. Prediction accuracies were high for BHBA models, and moderately high for NEFA models (R2 of external validation of 0.88 and 0.75, respectively). We identified 16 metabolites that were significantly (variable importance of projection score > 1) correlated with the concentration of one or both biomarkers. These metabolites were primarily intermediates of energy, phospholipid, and/or methyl donor metabolism. Of the significant metabolites identified; (1) two (acetate and creatine) were positively correlated with BHBA but negatively correlated with NEFA, (2) nine had similar associations with both BHBA and NEFA, (3) two were correlated with only BHBA concentration, and (4) three were only correlated with NEFA concentration. Overall, our results suggest that BHBA and NEFA are indicative of similar metabolic states in clinically healthy animals, but that several significant metabolic differences exist that help to explain the weak correlations between them. We also identified several metabolites that may be useful intermediate phenotypes in genomic selection for improved metabolic health.