Short communication: ketone body concentration in milk determined by Fourier transform infrared spectroscopy: value for the detection of hyperketonemia in dairy cows

Genome-Wide Association Study of Milk Composition in Karachai Goats

This study is first to perform a genome-wide association study (GWAS) to investigate the milk quality traits in Karachai goats. The objective of the study was to identify candidate genes associated with milk composition traits based on the identification and subsequent analysis of all possible SNPs, both genome-wide (high-confidence) and suggestive (subthreshold significance). To estimate the milk components, 22 traits were determined, including several types of fatty acids. DNA was extracted from ear tissue or blood samples. A total of 167 Karachai goats were genotyped using an Illumina GoatSNP53K BeadChip panel (Illumina Inc., San Diego, CA, USA). Overall, we identified 167 highly significant and subthreshold SNPs associated with the milk components of Karachai goats. A total of 10 SNPs were located within protein-coding genes and 33 SNPs in close proximity to them (±0.2 Mb). The largest number of genome-wide significant SNPs was found on chromosomes 2 and 8 and some of them were associated with several traits. The greatest number of genome-wide significant SNPs was identified for crude protein and lactose (6), and the smallest number-only 1 SNP-for freezing point depression. No SNPs were identified for monounsaturated and polyunsaturated fatty acids. Functional annotation of all 43 SNPs allowed us to identify 66 significant candidate genes on chromosomes 1, 2, 3, 4, 5, 8, 10, 13, 16, 18, 21, 23, 25, 26, and 27. We considered these genes potential DNA markers of the fatty acid composition of Karachai goat milk. Also, we found 12 genes that had a polygenic effect: most of them were simultaneously associated with the dry matter content and fatty acids (METTL, SLC1A 8, PHACTR1, FMO2, ECI1, PGP, ABCA3, AMDHD2). Our results suggest that the genes identified in our study affecting the milk components in Karachai goats differed from those identified in other breeds of dairy goats.

Prediction of metabolic status of dairy cows in early lactation using milk fatty acids and test-day variables

Early lactation metabolic imbalance is an important physiological change affecting the health, production, and reproduction of dairy cows. The aims of this study were (1) to evaluate the potential of test-day (TD) variables with or without milk fatty acids (FA) content to classify metabolically imbalanced cows and (2) to evaluate the robustness of the metabolic classification with external data. A data set was compiled from 3 experiments containing plasma β-hydroxybutyrate, nonesterified FA, glucose, insulin-like growth factor-I, FA proportions in milk fat, and TD variables collected from 244 lactations in wk 2 after calving. Based on the plasma metabolites, 3 metabolic clusters were identified using fuzzy c-means clustering and the probabilistic membership value of each cow to the 3 clusters was determined. Comparing the mean concentration of the plasma metabolites, the clusters were differentiated into metabolically imbalanced, moderately impacted, and balanced. Following this, the 2 metabolic status groups identified were imbalanced cows (n = 42), which were separated from what we refer to as "others" (n = 202) based on the membership value of each cow for the imbalanced cluster using a threshold of 0.5. The following 2 FA data sets were composed: (1) FA (groups) having high prediction accuracy by Fourier-transform infrared spectroscopy and, thus, have practical significance, and (2) FA (groups) formerly identified as associated with metabolic changes in early lactation. Metabolic status prediction models were built using FA alone or combined with TD variables as predictors of metabolic groups. Comparison was made among models and external evaluations were performed using an independent data set of 115 lactations. The area under the receiver operating characteristics curve of the models was between 75 and 91%, indicating their moderate to high accuracy as a diagnostic test for metabolic imbalance. The addition of FA groups to the TD models enhanced the accuracy of the models. Models with FA and TD variables showed high sensitivities (80-88%). Specificities of these models (73-79%) were also moderate and acceptable. The accuracy of the FA models on the external data set was high (area under the receiver operating characteristics curve between 76 and 84). The persistently good performance of models with Fourier-transform infrared spectroscopy-quantifiable FA on the external data set showed their robustness and potential for routine screening of metabolically imbalanced cows in early lactation.

Importance of Mid-Infrared Spectra Regions for the Prediction of Mastitis and Ketosis in Dairy Cows

Mid-infrared (MIR) spectroscopy is routinely applied to determine major milk components, such as fat and protein. Moreover, it is used to predict fine milk composition and various traits pertinent to animal health. MIR spectra indicate an absorbance value of infrared light at 1060 specific wavenumbers from 926 to 5010 cm−1. According to research, certain parts of the spectrum do not contain sufficient information on traits of dairy cows. Hence, the objective of the present study was to identify specific regions of the MIR spectra of particular importance for the prediction of mastitis and ketosis, performing variable selection analysis. Partial least squares discriminant analysis (PLS-DA) along with three other statistical methods, support vector machine (SVM), least absolute shrinkage and selection operator (LASSO), and random forest (RF), were compared. Data originated from the Austrian milk recording and associated health monitoring system (GMON). Test-day data and corresponding MIR spectra were linked to respective clinical mastitis and ketosis diagnoses. Certain wavenumbers were identified as particularly relevant for the prediction models of clinical mastitis (23) and ketosis (61). Wavenumbers varied across four distinct statistical methods as well as concerning different traits. The results indicate that variable selection analysis could potentially be beneficial in the process of modeling.

A meta-analysis of the genetic contribution estimates to major indicators for ketosis in dairy cows

The present study aimed to perform a meta-analysis using the random-effects model to merge published genetic parameter estimates for major indicators of ketosis [milk concentrations of acetone (ACETm) and β-hydroxybutyrate (BHBAm), and blood concentration of β-hydroxybutyrate (BHBAb)] in dairy cows. Overall, 51 heritability estimates and 130 genetic correlations from 19 papers published between 2012 and 2022 were used in this study. The average heritability estimates for ACETm, BHBAm, and BHBAb were 0.164, 0.123, and 0.141, respectively. The genetic correlation estimates between BHBAm and milk yield (MY), milk protein percentage (PP), and body condition score (BCS) were negative and moderate (-0.252, -0.200, and - 0.314, respectively). Genetic correlation estimates between BHBAm and milk fat percentage (FP), milk fat to protein ratio (FPR), and ketosis (KET) were moderate to high (0.411, 0.512, and 0.614, respectively). The genetic correlation estimates between BHBAb and MY and FP were low and equal to 0.128 and 0.035, respectively. The genetic correlation estimates between ACETm-MY and ACETm-PP were negative and moderate (-0.374 and - 0.398, respectively). Estimates of genetic correlation between ACETm and FP, FPR, and KET were moderate to high (0.455, 0.626, and 0.876, respectively). The results of this meta-analysis indicated the existence of additive genetic variation for ketosis indicator metabolites which could be exploited in genetic selection programs to reduce ketosis in dairy cows. Moreover, the results propose that selection for lower concentrations of indicator traits could be an effective plan for indirect improvement of production and reproduction performance, and health in dairy cows.

Prediction of Acute and Chronic Mastitis in Dairy Cows Based on Somatic Cell Score and Mid-Infrared Spectroscopy of Milk

Simple Summary

Mid-infrared (MIR) spectroscopy is the method of choice to determine milk components like fat, protein and urea. We examined the potential of MIR spectra analyses for the prediction of clinical mastitis events of dairy cows additionally, or alternatively, to somatic cell count, which is routinely used as an indicator for mastitis monitoring. Prediction models based on MIR spectra and a somatic cell count-derived score (SCS) were developed and compared. A model based on MIR spectra and SCS proved more accurate at predicting mastitis than models based on either indicator alone. Consequently, MIR spectra analyses add extra value in the prediction of clinical mastitis, making them potentially useful for dairy farm management and as an auxiliary trait for the genetic evaluation of udder health.

Abstract

Monitoring for mastitis on dairy farms is of particular importance, as it is one of the most prevalent bovine diseases. A commonly used indicator for mastitis monitoring is somatic cell count. A supplementary tool to predict mastitis risk may be mid-infrared (MIR) spectroscopy of milk. Because bovine health status can affect milk composition, this technique is already routinely used to determine standard milk components. The aim of the present study was to compare the performance of models to predict clinical mastitis based on MIR spectral data and/or somatic cell count score (SCS), and to explore differences of prediction accuracies for acute and chronic clinical mastitis diagnoses. Test-day data of the routine Austrian milk recording system and diagnosis data of its health monitoring, from 59,002 cows of the breeds Fleckvieh (dual purpose Simmental), Holstein Friesian and Brown Swiss, were used. Test-day records within 21 days before and 21 days after a mastitis diagnosis were defined as mastitis cases. Three different models (MIR, SCS, MIR + SCS) were compared, applying Partial Least Squares Discriminant Analysis. Results of external validation in the overall time window (−/+21 days) showed area under receiver operating characteristic curves (AUC) of 0.70 when based only on MIR, 0.72 when based only on SCS, and 0.76 when based on both. Considering as mastitis cases only the test-day records within 7 days after mastitis diagnosis, the corresponding areas under the curve were 0.77, 0.83 and 0.85. Hence, the model combining MIR spectral data and SCS was performing best. Mastitis probabilities derived from the prediction models are potentially valuable for routine mastitis monitoring for farmers, as well as for the genetic evaluation of the trait udder health.

The Use of Multilayer Perceptron Artificial Neural Networks to Detect Dairy Cows at Risk of Ketosis

Simple Summary

Ketosis is a serious metabolic disease in high-yield dairy cows, that affects productive herds throughout the world. Subclinical ketosis is one of the most dominant metabolic disorders in dairy herds during early lactation, so early detection and prevention are important for both economic and animal welfare reasons. Neural networks, which offer a high degree of accuracy in predicting various phenomena and processes where there is no clear causal correlation or there are no rules that allow the establishment of a logical cause-and-effect relationship, can be used to address problems related to prediction, classification, or control. A Multi-Layer perceptron (MLP) is a feedforward artificial neural network model that takes input data for a set of proper output. This study investigated the performance of four algorithms used to train MLP networks. The experimental results demonstrate that the MLP network model improved the accuracy of process recognition of subclinical ketosis in dairy cows. The received artificial model’s results were saved in the predictive model markup language (PMML) and can be used to describe the learning set, the algorithm used in the data mining application and related information.

Abstract

Subclinical ketosis is one of the most dominant metabolic disorders in dairy herds during lactation. Cows suffering from ketosis experience elevated ketone body levels in blood and milk, including β-hydroxybutyric acid (BHB), acetone (ACE) and acetoacetic acid. Ketosis causes serious financial losses to dairy cattle breeders and milk producers due to the costs of diagnosis and management as well as animal welfare reasons. Recent years have seen a growing interest in the use of artificial neural networks (ANNs) in various fields of science. ANNs offer a modeling method that enables the mapping of highly complex functional relationships. The purpose of this study was to determine the relationship between milk composition and blood BHB levels associated with subclinical ketosis in dairy cows, using feedforward multilayer perceptron (MLP) artificial neural networks. The results were verified based on the estimated sensitivity and specificity of selected network models, an optimum cut-off point was identified for the receiver operating characteristic (ROC) curve and the area under the ROC curve (AUC). The study demonstrated that BHB, ACE and lactose (LAC) levels, as well as the fat-to-protein ratio in milk, were important input variables in the network training process. For the identification of cows at risk of subclinical ketosis, variables such as BHB and ACE levels in milk were of particular relevance, with a sensitivity and specificity of 0.84 and 0.61, respectively. It was found that the back propagation algorithm offers opportunities to integrate artificial intelligence and dairy cattle welfare within a computerized decision support tool.

Genetic analyses of blood β-hydroxybutyrate predicted from milk infrared spectra and its association with longevity and female reproductive traits in Holstein cattle

Ketosis is one of the most prevalent and complex metabolic disorders in high-producing dairy cows and usually detected through analyses of β-hydroxybutyrate (BHB) concentration in blood. Our main objectives were to evaluate genetic parameters for blood BHB predicted based on Fourier-transform mid-infrared spectra from 5 to 305 d in milk, and estimate the genetic relationships of blood BHB with 7 reproduction traits and 6 longevity traits in Holstein cattle. Predicted blood BHB records of 11,609 Holstein cows (after quality control) were collected from 2016 to 2019 and used to derive 4 traits based on parity number, including predicted blood BHB in all parities (BHBp), parity 1 (BHB1), parity 2 (BHB2), and parity 3+ (BHB3). Single- and multitrait repeatability models were used for estimating genetic parameters for the 4 BHB traits. Random regression test-day models implemented via Bayesian inference were used to evaluate the daily genetic feature of BHB variability. In addition, genetic correlations were calculated for the 4 BHB traits with reproduction and longevity traits. The heritability estimates of BHBp, BHB1, BHB2, and BHB3 ranged from 0.100 ± 0.026 (± standard error) to 0.131 ± 0.023. The BHB in parities 1 to 3+ were highly genetically correlated and ranged from 0.788 (BHB1 and BHB2) to 0.911 (BHB1 and BHB3). The daily heritability of BHBp ranged from 0.069 to 0.195, higher for the early and lower for the later lactation periods. A similar trend was observed for BHB1, BHB2, and BHB3. There are low direct genetic correlations between BHBp and selected reproductive performance and longevity traits, which ranged from -0.168 ± 0.019 (BHBp and production life) to 0.157 ± 0.019 (BHBp and age at first calving) for the early lactation stage (5 to 65 d). These direct genetic correlations indicate that cows with higher BHBp (greater likelihood of having ketosis) in blood usually have shorter production life (-0.168 ± 0.019). Cows with higher fertility and postpartum recovery, such as younger age at first calving (0.157 ± 0.019) and shorter interval from calving to first insemination in heifer (0.111 ± 0.006), usually have lower BHB concentration in the blood. Furthermore, the direct genetic correlations change across parity and lactation stage. In general, our results suggest that selection for lower predicted BHB in early lactation could be an efficient strategy for reducing the incidence of ketosis as well as indirectly improving reproductive and longevity performance in Holstein cattle.

Developing a predictive model for beta-hydroxybutyrate and non-esterified fatty acids using milk fourier-transform infrared spectroscopy in dairy cows

Negative energy balance following parturition predisposes dairy cattle to numerous metabolic disorders. Current diagnostics are limited by their labor requirements and inability to measure multiple metabolic markers simultaneously. Fourier-transform Infrared spectroscopy (FTIR) data, measured from milk samples, could improve the detection of metabolic disorders using routine milk samples from dairy farms. The objective of this study was to develop a predictive model for numeric values of blood beta-hydroxybutyrate (BHB) and blood non-esterified fatty acids (NEFA). The study utilized a dataset comprised of 622 observations with known blood BHB and blood NEFA values measured concurrently with the milk FTIR data. Using ElasticNet regression on milk FTIR data and production information, we built regression models for numeric blood BHB and blood NEFA prediction and classification models for blood BHB values greater than 1.2 mmol/L and blood NEFA values greater than 0.7 mmol/L. The R² of the best fitting model was 0.56 and 0.51 for log-transformed BHB and log-transformed NEFA, respectively. The BHB classification model had a 90 % sensitivity and 83 % specificity and the NEFA classification model achieved a sensitivity of 73 % and specificity of 74 %. We applied our numeric prediction models to an external dataset (n = 9660) with known blood metabolites to validate the prediction accuracy of log-transformed blood BHB and log-transformed blood NEFA models. Log-transformed BHB root mean square error (RMSE) was 0.4018 and log-transformed NEFA RMSE was 0.4043. The second objective of this study was to develop a categorization for cows as either metabolically disordered or healthy. Responses to negative energy balance in transition cows are related to blood levels of BHB and NEFA. Cows suffering from metabolic disorders without elevated blood BHB values remain unidentified when detection is focused on blood BHB alone. To account for this differentiated metabolic response, we categorized cows as either 'metabolically healthy' or suffering a 'metabolic disorder' by using a combination of blood BHB, blood NEFA, and milk fat to protein quotient. We obtained a balanced accuracy of 94 % for the prediction of cow metabolic status. Direct prediction of metabolic status can be used to identify hyperketonemic cows in addition to cows exhibiting metabolic response patterns not detected by elevated blood BHB alone.

Characterization of ketolactia in dairy cows during early lactation

Fourier transform infrared spectroscopy (FTIR) allows for the determination of milk acetone (mACE) and β-hydroxybutyrate (mBHB) concentrations, providing a potential herd monitoring tool for hyperketolactia, defined as elevated milk ketone bodies. The study aim was to characterize mACE and mBHB concentration dynamics during early lactation in Polish Holstein-Friesian cows. Milk samples (n = 3,867,390) were collected within 6 to 60 days in milk (DIM) over a 4-yr period (April 1, 2013 to March 31, 2017) from approximately 21,300 dairy herds (average 38.7 cows/herd). Fixed effects of parity, DIM, and their interaction on mACE and mBHB concentrations were determined using a mixed model with a herd-year-season fixed effect and random cow effect. Published hyperketolactic mACE (≥0.15 mmol/L) and mBHB (≥0.10 mmol/L) threshold concentrations were used to classify study milk samples into ketolactia groups of normal (mACE <0.15 mmol/L and mBHB <0.10 mmol/L) and hyperketolactic (HYKL; either mACE ≥0.15 mmol/L or mBHB ≥0.10 mmol/L). Additionally, HYKL samples were categorized into subpopulations as having elevated mBHB and mACE (HYKL_ACEBHB, mACE ≥0.15 mmol/L and mBHB ≥0.10 mmol/L), only elevated mBHB (HYKL_BHB; mACE <0.15 mmol/L and mBHB ≥0.10 mmol/L), or only elevated mACE (HYKL_ACE; mACE ≥0.15 mmol/L and mBHB <0.10 mmol/L). Effects of parity, DIM, ketolactia group or subpopulation, and their interactions on mACE and mBHB concentrations were also determined using the mixed model that included ketolactia group or subpopulation as an independent variable. Across all data, mACE and mBHB concentrations were influenced by effects of parity, DIM, and their interaction as well as parity, DIM, ketolactia group or subpopulation, and their interactions. For all samples, mACE and mBHB concentrations decreased with increasing DIM, with mACE concentration declining more rapidly compared with mBHB. In the data set, 68% and 32% of all samples were defined as normal or HYKL, respectively. Among HYKL samples, mACE was elevated soon after calving and declined over time. In contrast, mBHB started lower after calving and increased reaching peak concentrations around 30 DIM, and then decreased. Within HYKL samples, 50.8, 41.3, and 7.9% were categorized as HYKL_ACEBHB, HYKL_BHB, and HYKL_ACE respectively. Between 6 and 21 DIM, 11.3% of HYKL were classified as HYKL_ACE. Primiparous cows had greater (14.8%) HYKL_ACE samples in this time period. In conclusion, this study has characterized mACE and mBHB concentrations during early lactation and determined effects of parity, DIM, and their interaction. Using published criteria interpreting mACE and mBHB concentrations, it was intriguing to identify a unique population of samples having elevated mACE without mBHB in early lactation, especially in primiparous cows. Further research is needed to determine if this sample population represents an unhealthy metabolic status that adversely affects cow health and performance.

Raw Cow Milk Protein Stability under Natural and Technological Conditions of Environment by Analysis of Variance

Heat stability (HS) is substantial technology property of raw milk. Analysis of sources of HS variation and its regular monitoring can contribute to creating higher added value in the dairy industry. The goal of this analysis was to assess the practice sources of raw cow milk HS variability on the results of an extensive data set of bulk tank milk samples. There was implemented neither a compositional technology modification nor acidity adjustment of milk, just original raw milk was used for the analysis. A total 2634 HS analyses were performed, including other milk indicators, during three years of an experimental period. The log HS mean and standard deviation were 1.273654 ± 0.144189, equal to the HS geometric mean of 18.8 min. Explanation of the HS variability through the linear model used was 41.1% (p < 0.0001). According to the results of the variance analysis, the milk HS was influenced (p = 0.0033 and mostly <0.0001) by all the farm factors such as year; season; calendar month; altitude; total annual rainfall; herd size by the number of cows; milk yield; cow breed; type of milking; litter type in the stable; summer grazing application; farm effect. During the calendar months (p < 0.0001), milk HS values suggest similar seasonal dynamics with the somatic cell count, total count of mesophilic microorganisms, coli bacteria count and urea and lactose concentration and opposite configuration pattern to fat, crude protein, solids-not-fat and total solids content and milk freezing point depression. Here performed quantification of these effects by analyzing the variance may allow efficient raw milk selection to be processed into specific dairy products.

Predicting Subclinical Ketosis in Dairy Cows Using Machine Learning Techniques

Simple Summary

The maintenance of cows in good health and physical condition is an important component of dairy cattle management. One of the major metabolic disorders in dairy cows is subclinical ketosis. Due to financial and organizational reasons it is often impossible to test all cows in a herd for ketosis using standard blood examination method. Using milk data from test-day records, obtained without additional costs for breeders, we found diagnostic models identifying cows-at-risk of subclinical ketosis. In addition, to select the best models, we present a general scoring approach for various machine learning models. With our models, breeders can identify dairy cows-at-risk of subclinical ketosis and implement appropriate management strategies and prevent losses in milk production.

Abstract

The diagnosis of subclinical ketosis in dairy cows based on blood ketone bodies is a challenging and costly procedure. Scientists are searching for tools based on results of milk performance assessment that would allow monitoring the risk of subclinical ketosis. The objective of the study was (1) to design a scoring system that would allow choosing the best machine learning models for the identification of cows-at-risk of subclinical ketosis, (2) to select the best performing models, and (3) to validate them using a testing dataset containing unseen data. The scoring system was developed using two machine learning modeling pipelines, one for regression and one for classification. As part of the system, different feature selections, outlier detection, data scaling and oversampling methods were used. Various linear and non-linear models were fit using training datasets and evaluated on holdout, testing the datasets. For the assessment of suitability of individual models for predicting subclinical ketosis, three β-hydroxybutyrate concentration in blood (bBHB) thresholds were defined: 1.0, 1.2 and 1.4 mmol/L. Considering the thresholds of 1.2 and 1.4, the logistic regression model was found to be the best fitted model, which included independent variables such as fat-to-protein ratio, acetone and β-hydroxybutyrate concentrations in milk, lactose percentage, lactation number and days in milk. In the cross-validation, this model showed an average sensitivity of 0.74 or 0.75 and specificity of 0.76 or 0.78, at the pre-defined bBHB threshold 1.2 or 1.4 mmol/L, respectively. The values of these metrics were also similar in the external validation on the testing dataset (0.72 or 0.74 for sensitivity and 0.80 or 0.81 for specificity). For the bBHB threshold at 1.0 mmol/L, the best classification model was the model based on the SVC (Support Vector Classification) machine learning method, for which the sensitivity in the cross-validation was 0.74 and the specificity was 0.73. These metrics had lower values for the testing dataset (0.57 and 0.72 respectively). Regression models were characterized by poor fitness to data (R² < 0.4). The study results suggest that the prediction of subclinical ketosis based on data from test-day records using classification methods and machine learning algorithms can be a useful tool for monitoring the incidence of this metabolic disorder in dairy cattle herds.

Milk Beta-Hydroxybutyrate and Fat to Protein Ratio Patterns during the First Five Months of Lactation in Holstein Dairy Cows Presenting Treated Left Displaced Abomasum and Other Post-Partum Diseases

Simple Summary

This study aimed to evaluate the 5-month pattern (averaged days in milk; DIM1 to 5) of milk beta-hydroxybutyrate (BHB) concentration and fat to protein content (F:P) ratio patterns from Holstein cows presenting postpartum diseases which have been treated. Cows presenting left displaced abomasum (LDA) and concomitant diseases within the first three months had higher concentrations of BHB than the control group (cows without diseases) in the first, but not in the second month postpartum. The F:P ratio had a similar evolution pattern also for DIM2. Animals with LDA were four to six times more likely to have a F:P ratio ≥ 1.29 than the control group during DIM1 and DIM2, respectively. Moderate and high correlations were also observed between the F:P ratio and BHB in DIM1 and DIM2, respectively. We concluded that animals suffering from LDA within the first three months postpartum have a significantly higher concentration of BHB and F:P ratio in milk than cows without postpartum diseases during the first two months. The treated cows with LDA quickly recovered normal levels, up to DIM3. The F:P ratio is a viable and economic indicator, mainly between the first two months postpartum, to estimate BHB concentration and energy balance in cows presenting LDA and in recovery.

Abstract

The main objective of the present study was to evaluate the beta-hydroxybutyrate (BHB) and fat to protein content (F:P) ratio patterns in the milk of Holstein cows with postpartum diseases throughout the first five months of lactation. This prospective study was performed at Vestjyske Dyrlaeger ApS (Nørre Nebel, Denmark). The milk fat, protein, and BHB were evaluated in the Danish Eurofins laboratory according to the monthly averaged days in milk (DIM1 to 5). According to clinical records, five groups were formed: A (control group; cows without diseases; n = 32), B (cows with left displaced abomasum -LDA- and concomitant diseases; n = 25); C (cows with other diseases up to DIM3; n = 13); D (cows with foot disorders up to DIM3; n = 26); and E (cows with disease manifestations in DIM4 and DIM5; n = 26). All the sick cows were treated after diagnosis, and laparoscopy was performed on cows with LDA. In group B, a higher concentration of BHB (0.18 ± 0.02 mmol/L; p < 0.001) was observed than in the control group (0.07 ± 0.02 mmol/L; p < 0.001) in DIM1, presenting an odds ratio (OR) = 8.9. In all groups, BHB decreased to 0.03–0.05 mmol/L (p < 0.05) since DIM3. The F:P ratio was higher in group B (1.77 ± 0.07) than in group A (1.32 ± 0.06; p < 0.05) in DIM1. A similar profile is observed in DIM2. It was observed that animals in group B were four to six times more likely to have a F:P ratio ≥1.29 during DIM1 (OR = 4.0; 95% CI:1.3–14.4; p = 0.01) and DIM2 (OR = 5.9; 95% CI %:1.9–21.9; p < 0.01), than cows in group A. There were also moderate and high correlations between the F:P ratio and the BHB for DIM1 (r = 0.57; r² = 0.33; RSD = 0.09; p < 0.001) and DIM2 (r = 0.78; r² = 0.60; RSD = 0.07; p < 0.001), respectively. We concluded that animals affected by LDA in the postpartum period have a higher concentration of BHB in milk in DIM1 and all treated animals quickly recover BHB levels up to DIM3. The F:P ratio is a viable and economic indicator, mainly in DIM1 and DIM2, to estimate BHB concentration and energy balance in cows with LDA and other postpartum diseases.

Genetic and nongenetic profiling of milk β-hydroxybutyrate and acetone and their associations with ketosis in Holstein cows

Ketosis is a metabolic disorder of increasing importance in high-yielding dairy cows, but accurate population-wide binary health trait recording is difficult to implement. Against this background, proper Gaussian indicator traits, which can be routinely measured in milk, are needed. Consequently, we focused on the ketone bodies acetone and β-hydroxybutyrate (BHB), measured via Fourier-transform infrared spectroscopy (FTIR) in milk. In the present study, 62,568 Holstein cows from large-scale German co-operator herds were phenotyped for clinical ketosis (KET) according to a veterinarian diagnosis key. A sub-sample of 16,861 cows additionally had first test-day observations for FTIR acetone and BHB. Associations between FTIR acetone and BHB with KET and with test-day traits were studied phenotypically and quantitative genetically. Furthermore, we estimated SNP marker effects for acetone and BHB (application of genome-wide association studies) based on 40,828 SNP markers from 4,384 genotyped cows, and studied potential candidate genes influencing body fat mobilization. Generalized linear mixed models were applied to infer the influence of binary KET on Gaussian-distributed acetone and BHB (definition of an identity link function), and vice versa, such as the influence of acetone and BHB on KET (definition of a logit link function). Additionally, linear models were applied to study associations between BHB, acetone and test-day traits (milk yield, fat percentage, protein percentage, fat-to-protein ratio and somatic cell score) from the first test-day after calving. An increasing KET incidence was statistically significant associated with increasing FTIR acetone and BHB milk concentrations. Acetone and BHB concentrations were positively associated with fat percentage, fat-to-protein ratio and somatic cell score. Bivariate linear animal models were applied to estimate genetic (co)variance components for KET, acetone, BHB and test-day traits within parities 1 to 3, and considering all parities simultaneously in repeatability models. Pedigree-based heritabilities were quite small (i.e., in the range from 0.01 in parity 3 to 0.07 in parity 1 for acetone, and from 0.03-0.04 for BHB). Heritabilites from repeatability models were 0.05 for acetone, and 0.03 for BHB. Genetic correlations between acetone and BHB were moderate to large within parities and considering all parities simultaneously (0.69-0.98). Genetic correlations between acetone and BHB with KET from different parities ranged from 0.71 to 0.99. Genetic correlations between acetone across parities, and between BHB across parities, ranged from 0.55 to 0.66. Genetic correlations between KET, acetone, and BHB with fat-to-protein ratio and with fat percentage were large and positive, but negative with milk yield. In genome-wide association studies, we identified SNP on BTA 4, 10, 11, and 29 significantly influencing acetone, and on BTA 1 and 16 significantly influencing BHB. The identified potential candidate genes NRXN3, ACOXL, BCL2L11, HIBADH, KCNJ1, and PRG4 are involved in lipid and glucose metabolism pathways.

Comparison of Fatty Acid Proportions Determined by Mid-Infrared Spectroscopy and Gas Chromatography in Bulk and Individual Milk Samples

Simple Summary

Information about fatty acid proportions in milk fat is important for many purposes, such as animal breeding, animal health control, as well as human nutrition. The routine methods for determining fatty acid proportions (e.g., mid-infrared spectroscopy) are rapid and relatively cheap, but there is a need to compare them with the reference analytical method (gas chromatography) to ensure their validity and suitability for various milk samples. The aim of this study is to compare the proportions of single fatty acids and their sums determined by utilizing both of these analytical methods and the resulting correlation coefficients. Our results show that the mid-infrared spectroscopy method is more appropriate (both for bulk and individual milk samples) for fatty acids present in high proportions of the total fat and for the sum of fatty acids (such as saturated and unsaturated) than for fatty acids with low proportions.

Abstract

Rapid analytical methods can contribute to the expansion of milk fatty acid determination for various important practical purposes. The reliability of data resulting from these routine methods plays a crucial role. Bulk and individual milk samples (60 and 345, respectively) were obtained from Czech Fleckvieh and Holstein dairy cows in the Czech Republic. The correlation between milk fatty acid (FA) proportions determined by the routine method (infrared spectroscopy in the mid-region in connection with Fourier transformation; FT-MIR) and the reference method (gas chromatography; GC) was evaluated. To validate the calibration of the FT-MIR method, a linear regression model was used. For bulk milk samples, the correlation coefficients between these methods were higher for the saturated (SFAs) and unsaturated FAs (UFAs) (r = 0.7169 and 0.9232; p < 0.001) than for the trans isomers of UFAs (TFAs) and polyunsaturated FAs (PUFAs) (r = 0.5706 and 0.6278; p < 0.001). Similar results were found for individual milk samples: r = 0.8592 and 0.8666 (p < 0.001) for SFAs and UFAs, 0.1690 (p < 0.01) for TFAs, and 0.3314 (p < 0.001) for PUFAs. The correlation coefficients for TFAs and PUFAs were statistically significant but too low for practical analytical application. The results indicate that the FT-MIR method can be used for routine determination mainly for SFAs and UFAs.

Short communication: Prediction of hyperketonemia in dairy cows in early lactation using on-farm cow data and net energy intake by partial least square discriminant analysis

The objectives of this study were (1) to evaluate if hyperketonemia in dairy cows (defined as plasma β-hydroxybutyrate ≥1.0 mmol/L) can be predicted using on-farm cow data either in current or previous lactation week, and (2) to study if adding individual net energy intake (NEI) can improve the predictive ability of the model. Plasma β-hydroxybutyrate concentration, on-farm cow data (milk yield, percentage of fat, protein and lactose, fat- and protein-corrected milk yield, body weight, body weight change, dry period length, parity, and somatic cell count), and NEI of 424 individual cows were available weekly through lactation wk 1 to 5 postpartum. To predict hyperketonemia in dairy cows, models were first trained by partial least square discriminant analysis, using on-farm cow data in the same or previous lactation week. Second, NEI was included in models to evaluate the improvement of the predictability of the models. Through leave-one trial-out cross-validation, models were evaluated by accuracy (the ratio of the sum of true positive and true negative), sensitivity (68.2% to 84.9%), specificity (61.5% to 98.7%), positive predictive value (57.7% to 98.7%), and negative predictive value (66.2% to 86.1%) to predict hyperketonemia of dairy cows. Through lactation wk 1 to 5, the accuracy to predict hyperketonemia using data in the same week was 64.4% to 85.5% (on-farm cow data only), 66.1% to 87.0% (model including NEI), and using data in the previous week was 58.5% to 82.0% (on-farm cow data only), 59.7% to 85.1% (model including NEI). An improvement of the accuracy of the model due to including NEI ranged among lactation weeks from 1.0% to 4.4% when using data in the same lactation week and 0.2% to 6.6% when using data in the previous lactation week. In conclusion, trained models via partial least square discriminant analysis have potential to predict hyperketonemia in dairy cows not only using data in the current lactation week, but also using data in the previous lactation week. Net energy intake can improve the accuracy of the model, but only to a limited extent. Besides NEI, body weight, body weight change, milk fat, and protein content were important variables to predict hyperketonemia, but their rank of importance differed across lactation weeks.

The use of milk Fourier transform mid-infrared spectra and milk yield to estimate heat production as a measure of efficiency of dairy cows

Background

Transformation of feed energy ingested by ruminants into milk is accompanied by energy losses via fecal and urine excretions, fermentation gases and heat. Heat production may differ among dairy cows despite comparable milk yield and body weight. Therefore, heat production can be considered an indicator of metabolic efficiency and directly measured in respiration chambers. The latter is an accurate but time-consuming technique. In contrast, milk Fourier transform mid-infrared (FTIR) spectroscopy is an inexpensive high-throughput method and used to estimate different physiological traits in cows. Thus, this study aimed to develop a heat production prediction model using heat production measurements in respiration chambers, milk FTIR spectra and milk yield measurements from dairy cows.

Methods

Heat production was computed based on the animal’s consumed oxygen, and produced carbon dioxide and methane in respiration chambers. Heat production data included 168 24-h-observations from 64 German Holstein and 20 dual-purpose Simmental cows. Animals were milked twice daily at 07:00 and 16:30 h in the respiration chambers. Milk yield was determined to predict heat production using a linear regression. Milk samples were collected from each milking and FTIR spectra were obtained with MilkoScan FT 6000. The average or milk yield-weighted average of the absorption spectra from the morning and afternoon milking were calculated to obtain a computed spectrum. A total of 288 wavenumbers per spectrum and the corresponding milk yield were used to develop the heat production model using partial least squares (PLS) regression.

Results

Measured heat production of studied animals ranged between 712 and 1470 kJ/kg BW^0.75. The coefficient of determination for the linear regression between milk yield and heat production was 0.46, whereas it was 0.23 for the FTIR spectra-based PLS model. The PLS prediction model using weighted average spectra and milk yield resulted in a cross-validation variance of 57% and a root mean square error of prediction of 86.5 kJ/kg BW^0.75. The ratio of performance to deviation (RPD) was 1.56.

Conclusion

The PLS model using weighted average FTIR spectra and milk yield has higher potential to predict heat production of dairy cows than models applying FTIR spectra or milk yield only.

Ketosis an Old Story Under a New Approach

Ketosis, characterized by high concentrations of ketone bodies in the blood, urine, and milk, affects a considerable number of cows immediately after calving. Although much is known about ketosis, dairy cows continue to be affected in every herd world-wide. Cows affected by ketosis are treated with palliative treatments after the disease is diagnosed. This is a very expensive approach and costs the dairy industry extra expenses, contributing to lower profitability of dairy herds. In this review article, we summarize the mainstream view on ketosis, classification of ketosis into three types, current diagnostic approaches to ketosis, and the economic impact of ketosis on dairy farms. Additionally, we discuss the most recent applications of the new ‘omics’ science of metabolomics in studying the etiopathology of ketosis as well as its contribution in identification of novel screening or diagnostic biomarkers of ketosis.

Opportunities and limitations of milk mid-infrared spectra-based estimation of acetone and β-hydroxybutyrate for the prediction of metabolic stress and ketosis in dairy cows

Subclinical (SCK) and clinical (CK) ketosis are metabolic disorders responsible for big losses in dairy production. Although Fourier-transform mid-infrared spectrometry (FTIR) to predict ketosis in cows exposed to great metabolic stress was studied extensively, little is known about its suitability in predicting hyperketonemia using individual samples, e.g. in small dairy herds or when only few animals are at risk of ketosis. The objective of the present research was to determine the applicability of milk metabolites predicted by FTIR spectrometry in the individual screening for ketosis. In experiment 1, blood and milk samples were taken every two weeks after calving from Holstein (n = 80), Brown Swiss (n = 72) and Swiss Fleckvieh (n = 58) cows. In experiment 2, cows diagnosed with CK (n = 474) and 420 samples with blood β-hydroxybutyrate [BHB] <1.0 mmol/l were used to investigate if CK could be detected by FTIR-predicted BHB and acetone from a preceding milk control. In experiment 3, correlations between data from an in farm automatic milk analyser and FTIR-predicted BHB and acetone from the monthly milk controls were evaluated. Hyperketonemia occurred in majority during the first eight weeks of lactation. Correlations between blood BHB and FTIR-predicted BHB and acetone were low (r = 0.37 and 0.12, respectively, P < 0.0001), as well as the percentage of true positive values (11.9 and 16.6%, respectively). No association of FTIR predicted ketone bodies with the interval of milk sampling relative to CK diagnosis was found. Data obtained from the automatic milk analyser were moderately correlated with the same day FTIR-predicted BHB analysis (r = 0.61). In conclusion, the low correlations with blood BHB and the small number of true positive samples discourage the use of milk mid-infrared spectrometry analyses as the only method to predict hyperketonemia at the individual cow level.

Short communication: Investigation of the temporal relationships between milk mid-infrared predicted biomarkers and lameness events in later lactation

This study reports on the exploration of temporal relationships between milk mid-infrared predicted biomarkers and lameness events. Lameness in dairy cows is an issue that can vary greatly in severity and is of concern for both producers and consumers. Metabolic disorders are often associated with lameness. However, lameness can arise weeks or even months after the metabolic disorder, making the detection of causality difficult. We already use mid-infrared technology to predict major milk components, such as fat or protein, during routine milk recording and for milk payment. It was recently shown that this technology can also be used to predict novel biomarkers linked to metabolic disorders in cows, such as oleic acid (18:1 cis-9), β-hydroxybutyrate, acetone, and citrate in milk. We used these novel biomarkers as proxies for metabolic issues. Other studies have explored the possibility of using mid-infrared spectra to predict metabolic diseases and found it (potentially) usable for indicating classes of metabolic problems. We wanted to explore the possible relationship between mid-infrared-based metabolites and lameness over the course of lactation. In total, data were recorded from 6,292 cows on 161 farms in Austria. Lameness data were recorded between March 2014 and March 2015 and consisted of 37,555 records. Mid-infrared data were recorded between July and December 2014 and consisted of 9,152 records. Our approach consisted of fitting preadjustments to the data using fixed effects, computing pair-wise correlations, and finally applying polynomial smoothing of the correlations for a given biomarker at a certain month in lactation and the lameness events scored on severity scale from sound or non-lame (lameness score of 1) to severely lame (lameness score of 5) throughout the lactation. The final correlations between biomarkers and lameness scores were significant, but not high. However, for the results of the present study, we should not look at the correlations in terms of absolute values, but rather as indicators of a relationship through time. When doing so, we can see that metabolic problems occurring in mo 1 and 3 seem more linked to long-term effects on hoof and leg health than those in mo 2. However, the quantity (only 1 pair-wise correlation exceeded 1,000 observations) and the quality (due to limited data, no separation according to more metabolic-related diseases could be done) of the data should be improved.

Symposium review: Big data, big predictions: Utilizing milk Fourier-transform infrared and genomics to improve hyperketonemia management

Negative animal health and performance outcomes are associated with disease incidences that can be labor-intensive, costly, and cumbersome for many farms. Amelioration of unfavorable outcomes through early detection and treatment of disease has emphasized the value of improving health monitoring. Although the value is recognized, detecting hyperketonemia (HYK) is still difficult for many farms to do practically and efficiently. Increasing data streams available to farms presents opportunities to use data to better monitor cow and herd health; however, challenges remain with regard to validating, integrating, and interpreting data. During the transition to lactation period, useful data are presented in the form of milk production and composition, milk Fourier-transform infrared (FTIR) wavelength absorbance, cow management records, and genomics, which have been employed to monitor postpartum onset of HYK. Attempts to predict postpartum HYK from test-day milk and performance variables incorporated into multiple linear regression models have demonstrated sufficient accuracy to monitor monthly herd prevalence; however, they lacked the sensitivity and specificity for individual cow diagnostics. Subsequent artificial neural network prediction models employing FTIR data and milk composition variables achieved 83 and 81% sensitivity and specificity for individual cow diagnostics. Although these results fail to reach the diagnostic goals of 90%, they are achieved without individual cow blood samples, which may justify acceptance of lower performance. The caveat is that these models depend on milk analysis, which is traditionally done every 4 weeks. This infrequent sampling allows for a single diagnostic sample for about half of the fresh cows. Benefits to farms are greatly improved if postpartum cows can be milk-tested weekly. Additionally, this allows for close monitoring of somatic cell count and may open the door for use of other herd health monitoring tools. Future improvements in these models may be achievable by maximizing sensitivity at the expense of specificity and may be most economical in disorders for which the cost of treatment is less than that of mistreating (e.g., HYK). Genomic predictions for HYK may be improved by incorporating genome-wide associated SNP and further utilized for precision management of HYK risk groups. Development and validation of HYK prediction models may provide producers with individual cow and herd-level management tools.

Disease structure of milk cows and the effect of the mass fractions ratio of fat and protein in milk on the level of the metabolites

The analysis of diseases occurring is given in animals in the first 100 DIM. Determined milk content: MF, MP, pH, Ur, BHBA, Ac, calculated FPR. INCD in milk cows in the first 100 DIM is the second most widespread – 30.72 %. Among INCD, DSD and RSD predominate in the former – 38.8 % each, DMEO accounts for 20.4 % of disease cases. In 17.05 % of the examined animals, the FPR corresponded to optimal values, and in 82.95 % it was 1.10 or less, which may indicate the spread of subacute subclinical rumen acidosis in the animals of the studied population. Exceeding the upper limit of FPR, indicating the presence of ketosis in animals, has not been established. Studies revealed a double excess of Ur content in milk, and in animals with normal FPR values, the Ur content was 11.15 % (p < 0.001) higher than in animals with reduced FPR. The pH of milk generally corresponded to the values of the physiological norm. The level of BHBA in milk was below the threshold values, but in the animal’s group with FPR 1.10 or less, the BHBA content in milk significantly (by 80.0 %, p < 0.01) exceeded the BHBA content in animals with normal FPR values. Cows with a normal FPR value, the Ac level in milk was found to exceed the threshold value by 28.57 %, and in animals with low FPR values, the established excess was 141.43 % (p <0.05).

Combined use of milk infrared spectra and genotypes can improve prediction of milk fat composition

It has been shown that milk infrared (IR) spectroscopy can be used to predict detailed milk fat composition. In addition, polymorphisms with substantial effects on milk fat composition have been identified. In this study, we investigated the combined use of milk IR spectroscopy and genotypes of dairy cows on the accuracy of predicting milk fat composition. Milk fat composition data based on gas chromatography and milk IR spectra were available for 1,456 Dutch Holstein Friesian cows. In addition, genotypes for the diacylglycerol acyltransferase 1 (DGAT1) K232A and stearoyl-CoA desaturase 1 (SCD1) A293V polymorphisms and a SNP located in an intron of the fatty acid synthase (FASN) gene were available. Adding SCD1 genotypes to the milk IR spectra resulted in a considerable improvement of the prediction accuracy for the unsaturated fatty acids C10:1, C12:1, C14:1 cis-9, and C16:1 cis-9 and their corresponding unsaturation indices. Adding DGAT1 genotypes to the milk IR spectra resulted in an improvement of the prediction accuracy for C16:1 cis-9 and C16 index. Adding genotypes of the FASN SNP to the IR spectra did not improve prediction of milk fat composition. This study demonstrated the potential of combining milk IR spectra with genotypic information from 3 polymorphisms to predict milk fat composition. We hypothesize that prediction accuracy of milk fat composition can be further improved by combining milk IR spectra with genomic breeding values.

Prediction of blood metabolites from milk mid-infrared spectra in early-lactation cows

Dairy cows commonly experience an unbalanced energy status in early lactation, and this condition can lead to the onset of several metabolic disorders. Blood metabolic profile testing is a valid tool to monitor and detect the most common early lactation disorders, but blood sampling and analysis are time-consuming and expensive, and the procedure is invasive and stressful for the cows. Mid-infrared (MIR) spectroscopy is routinely used to analyze milk composition, being a cost-effective and nondestructive method. The present study aimed to assess the feasibility of using routine milk MIR spectra for the prediction of main blood metabolites in dairy cows, and to investigate associations between measured blood metabolites and milk traits. Twenty herds of Holstein Friesian, Brown Swiss, or Simmental cows located in Northeast Italy were visited 1 to 4 times between December 2017 and June 2018, and blood and milk samples were collected from all lactating cows within 35 d in milk. Concentrations of main blood metabolites and milk MIR spectra were recorded from 295 blood and milk samples and used to develop prediction models for blood metabolic traits through backward interval partial least squares analysis. Blood β-hydroxybutyrate (BHB), urea, and nonesterified fatty acids were the most predictable traits, with coefficients of determination of 0.63, 0.58, and 0.52, respectively. On the contrary, predictive performance for blood glucose, triglycerides, cholesterol, glutamic oxaloacetic transaminase, and glutamic pyruvic transaminase were not accurate. Associations of blood BHB and urea with their respective contents in milk were moderate to strong, whereas all other correlations were weak. Predicted blood BHB showed an improved performance in detecting cows with hyperketonemia (blood BHB ≥ 1.2 mmol/L), compared with commercial calibration equation for milk BHB. Results highlighted the opportunity of using milk MIR spectra to predict blood metabolites and thus to collect routine information on the metabolic status of early-lactation cows at a population level.

Prediction of metabolic status of dairy cows in early lactation with on-farm cow data and machine learning algorithms

Metabolic status of dairy cows in early lactation can be evaluated using the concentrations of plasma β-hydroxybutyrate (BHB), free fatty acids (FFA), glucose, insulin, and insulin-like growth factor 1 (IGF-1). These plasma metabolites and metabolic hormones, however, are difficult to measure on farm. Instead, easily obtained on-farm cow data, such as milk production traits, have the potential to predict metabolic status. Here we aimed (1) to investigate whether metabolic status of individual cows in early lactation could be clustered based on their plasma values and (2) to evaluate machine learning algorithms to predict metabolic status using on-farm cow data. Through lactation wk 1 to 7, plasma metabolites and metabolic hormones of 334 cows were measured weekly and used to cluster each cow into 1 of 3 clusters per week. The cluster with the greatest plasma BHB and FFA and the lowest plasma glucose, insulin, and IGF-1 was defined as poor metabolic status; the cluster with the lowest plasma BHB and FFA and the greatest plasma glucose, insulin, and IGF-1 was defined as good metabolic status; and the intermediate cluster was defined as average metabolic status. Most dairy cows were classified as having average or good metabolic status, and a limited number of cows had poor metabolic status (10-50 cows per lactation week). On-farm cow data, including dry period length, parity, milk production traits, and body weight, were used to predict good or average metabolic status with 8 machine learning algorithms. Random Forest (error rate ranging from 12.4 to 22.6%) and Support Vector Machine (SVM; error rate ranging from 12.4 to 20.9%) were the top 2 best-performing algorithms to predict metabolic status using on-farm cow data. Random Forest had a higher sensitivity (range: 67.8-82.9% during wk 1 to 7) and negative predictive value (range: 89.5-93.8%) but lower specificity (range: 76.7-88.5%) and positive predictive value (range: 58.1-78.4%) than SVM. In Random Forest, milk yield, fat yield, protein percentage, and lactose yield had important roles in prediction, but their rank of importance differed across lactation weeks. In conclusion, dairy cows could be clustered for metabolic status based on plasma metabolites and metabolic hormones. Moreover, on-farm cow data can predict cows in good or average metabolic status, with Random Forest and SVM performing best of all algorithms.

Hyperketolactia occurrence before or after artificial insemination is associated with a decreased pregnancy per artificial insemination in dairy cows

The reproductive parameters of dairy cattle have continuously declined worldwide over the last 50 years. Nutritional imbalances are identified as risk factors for this decrease of reproductive performance. The present paper aims to quantify the decrease in the pregnancy per artificial insemination (P/AI) in the case of high milk ketones before and after AI. A total of 388,731 test-day from the Brittany Milk Recording Program in France from 226,429 cow-lactations were provided for this trial. For each test-day, information about lactation characteristics, date of AI, date of the following calving, and acetone and β-hydroxybutyrate (BHB) values were included. Ketones were predicted by Fourier transform mid-infrared spectroscopy using MilkoScan Foss analyzers (Foss, Hillerød, Denmark). Many thresholds were evaluated to define cows with hyperketolactia. Hyperketolactia statuses were then categorized into 1 of 4 possible classes according to the milk ketone dynamics for each AI and each threshold of acetone or BHB values (low-low, high-low, low-high, and high-high) within 20 d before and after AI. Similarly, the dynamics of udder health were characterized by changes in somatic cell counts measured at the same test day as ketone bodies. A logistic regression with a Poisson correction was performed to explain the relationship of P/AI with milk ketones and somatic cell count dynamics. Predicted acetone and BHB ranged from -0.51 to 4.92 mM (mean = 0.08 mM, SD = 0.10 mM) and -0.62 to 5.85 mM (mean = 0.07 mM, SD = 0.1 mM), respectively. Hyperketolactia defined by high acetone levels before AI was not associated with decreased P/AI, but high acetone levels after AI were associated with a >10% reduction in P/AI for all thresholds >0.10 mM. Hyperketolactia, defined by high BHB values before, after, or before and after AI, was associated with a 6 to 14% reduction in P/AI compared with cows with low BHB values. These associations are lower than those reported in previous trials in which blood ketones were used. High ketones in advanced lactation are likely to be the result of various primary disorders (secondary ketosis). Because the present work demonstrated that this situation is considered a risk factor for deteriorated reproductive performance, we suggest that high ketones in early and advanced lactation should be of interest to farm advisors.

Prediction of blood β-hydroxybutyrate content and occurrence of hyperketonemia in early-lactation, pasture-grazed dairy cows using milk infrared spectra

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 pasture-grazed 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²) 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 prediction, of 1.50. In the validation set, the R² 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.

Phenotypic relationships, genetic parameters, genome-wide associations, and identification of potential candidate genes for ketosis and fat-to-protein ratio in German Holstein cows

Energy demand for milk production in early lactation exceeds energy intake, especially in high-yielding Holstein cows. Energy deficiency causes increasing susceptibility to metabolic disorders. In addition to several blood parameters, the fat-to-protein ratio (FPR) is suggested as an indicator for ketosis, because a FPR >1.5 refers to high lipolysis. The aim of this study was to analyze phenotypic, quantitative genetic, and genomic associations between FPR and ketosis. In this regard, 8,912 first-lactation Holstein cows were phenotyped for ketosis according to a veterinarian diagnosis key. Ketosis was diagnosed if the cow showed an abnormal carbohydrate metabolism with increased content of ketone bodies in the blood or urine. At least one entry for ketosis in the first 6 wk after calving implied a score = 1 (diseased); otherwise, a score = 0 (healthy) was assigned. The FPR from the first test-day was defined as a Gaussian distributed trait (FPRgauss), and also as a binary response trait (FPRbin), considering a threshold of FPR = 1.5. After imputation and quality controls, 45,613 SNP markers from the 8,912 genotyped cows were used for genomic studies. Phenotypically, an increasing ketosis incidence was associated with significantly higher FPR, and vice versa. Hence, from a practical trait recording perspective, first test-day FPR is suggested as an indicator for ketosis. The ketosis heritability was slightly larger when modeling the pedigree-based relationship matrix (pedigree-based: 0.17; SNP-based: 0.11). For FPRbin, heritabilities were larger when modeling the genomic relationship matrix (pedigree-based: 0.09; SNP-based: 0.15). For FPRgauss, heritabilities were almost identical for both pedigree and genomic relationship matrices (pedigree-based: 0.14; SNP-based: 0.15). Genetic correlations between ketosis with FPRbin and FPRgauss using either pedigree- or genomic-based relationship matrices were in a moderate range from 0.39 to 0.71. Applying genome-wide association studies, we identified the specific SNP rs109896020 (BTA 5, position: 115,456,438 bp) significantly contributing to ketosis. The identified potential candidate gene PARVB in close chromosomal distance is associated with nonalcoholic fatty liver disease in humans. The most important SNP contributing to FPRbin was located within the DGAT1 gene. Different SNP significantly contributed to ketosis and FPRbin, indicating different mechanisms for both traits genomically.

Symposium review: Building a better cow-The Australian experience and future perspectives

Genomic selection has led to opportunities for developing new breeding values that rely on phenotypes in dedicated reference populations of genotyped cows. In Australia, it has been applied to 2 novel traits: feed efficiency, which was released in 2015 as feed saved breeding values, and heat tolerance genomic breeding values, released for the first time in 2017. Feed saved is already included in the national breeding objective, which is focused on profitability and designed to be in line with farmer preferences. Our future focus is on traits associated with animal health, either directly or in combination with predictor traits, such as mid-infrared spectral data and, into the future, automated data capture. Although it is common for many evaluated traits to have genomic reliabilities ranging between 60 and 75%, many new, genomic information-only traits are likely to have reliabilities of less than 50%. Pooling of phenotype data internationally and investing in maintenance of reference populations is one option to increase the reliability of these traits; the other is to apply improved genomic prediction methods. For example, advances in the use of sequence data, in addition to gene expression studies, can lead to improved persistence of genomic breeding values across breeds and generations and potentially lead to greater reliabilities. Lower genomic reliabilities of novel traits could reduce the overall index reliability. However, provided these traits contribute to the overall breeding objective (e.g., profit), they are worth including. Bull selection tools and personalized genetic trends are already available, but increased access to economic and automatic capture farm data may see even better use of data to improve farm management and selection decisions.

Predicting hyperketonemia by logistic and linear regression using test-day milk and performance variables in early-lactation Holstein and Jersey cows

Although cowside testing strategies for diagnosing hyperketonemia (HYK) are available, many are labor intensive and costly, and some lack sufficient accuracy. Predicting milk ketone bodies by Fourier transform infrared spectrometry during routine milk sampling may offer a more practical monitoring strategy. The objectives of this study were to (1) develop linear and logistic regression models using all available test-day milk and performance variables for predicting HYK and (2) compare prediction methods (Fourier transform infrared milk ketone bodies, linear regression models, and logistic regression models) to determine which is the most predictive of HYK. Given the data available, a secondary objective was to evaluate differences in test-day milk and performance variables (continuous measurements) between Holsteins and Jerseys and between cows with or without HYK within breed. Blood samples were collected on the same day as milk sampling from 658 Holstein and 468 Jersey cows between 5 and 20 d in milk (DIM). Diagnosis of HYK was at a serum β-hydroxybutyrate (BHB) concentration ≥1.2 mmol/L. Concentrations of milk BHB and acetone were predicted by Fourier transform infrared spectrometry (Foss Analytical, Hillerød, Denmark). Thresholds of milk BHB and acetone were tested for diagnostic accuracy, and logistic models were built from continuous variables to predict HYK in primiparous and multiparous cows within breed. Linear models were constructed from continuous variables for primiparous and multiparous cows within breed that were 5 to 11 DIM or 12 to 20 DIM. Milk ketone body thresholds diagnosed HYK with 64.0 to 92.9% accuracy in Holsteins and 59.1 to 86.6% accuracy in Jerseys. Logistic models predicted HYK with 82.6 to 97.3% accuracy. Internally cross-validated multiple linear regression models diagnosed HYK of Holstein cows with 97.8% accuracy for primiparous and 83.3% accuracy for multiparous cows. Accuracy of Jersey models was 81.3% in primiparous and 83.4% in multiparous cows. These results suggest that predicting serum BHB from continuous test-day milk and performance variables could serve as a valuable diagnostic tool for monitoring HYK in Holstein and Jersey herds.

Standardization of milk mid-infrared spectrometers for the transfer and use of multiple models

An increasing number of models are being developed to provide information from milk Fourier transform mid-infrared (FT-MIR) spectra on fine milk composition, technological properties of milk, or even cows' physiological status. In this context, and to take advantage of these existing models, the purpose of this work was to evaluate whether a spectral standardization method can enable the use of multiple equations within a network of different FT-MIR spectrometers. The piecewise direct standardization method was used, matching "slave" instruments to a common reference, the "master." The effect of standardization on network reproducibility was assessed on 66 instruments from 3 different brands by comparing the spectral variability of the slaves and the master with and without standardization. With standardization, the global Mahalanobis distance from the slave spectra to the master spectra was reduced on average from 2,655.9 to 14.3, representing a significant reduction of noninformative spectral variability. The transfer of models from instrument to instrument was tested using 3 FT-MIR models predicting (1) the quantity of daily methane emitted by dairy cows, (2) the concentration of polyunsaturated fatty acids in milk, and (3) the fresh cheese yield. The differences, in terms of root mean squared error, between master predictions and slave predictions were reduced after standardization on average from 103 to 17 g/d, from 0.0315 to 0.0045 g/100 mL of milk, and from 2.55 to 0.49 g of curd/100 g of milk, respectively. For all the models, standard deviations of predictions among all the instruments were also reduced by 5.11 times for methane, 5.01 times for polyunsaturated fatty acids, and 7.05 times for fresh cheese yield, showing an improvement of prediction reproducibility within the network. Regarding the results obtained, spectral standardization allows the transfer and use of multiple models on all instruments as well as the improvement of spectral and prediction reproducibility within the network. The method makes the models universal, thereby offering opportunities for data exchange and the creation and use of common robust models at an international level to provide more information to the dairy sector from direct analysis of milk.

Invited review: Opportunities for genetic improvement of metabolic diseases

Metabolic disorders are disturbances to one or more of the metabolic processes in dairy cattle. Dysfunction of any of these processes is associated with the manifestation of metabolic diseases or disorders. In this review, data recording, incidences, genetic parameters, predictors, and status of genetic evaluations were examined for (1) ketosis, (2) displaced abomasum, (3) milk fever, and (4) tetany, as these are the most prevalent metabolic diseases where published genetic parameters are available. The reported incidences of clinical cases of metabolic disorders are generally low (less than 10% of cows are recorded as having a metabolic disease per herd per year or parity/lactation). Heritability estimates are also low and are typically less than 5%. Genetic correlations between metabolic traits are mainly positive, indicating that selection to improve one of these diseases is likely to have a positive effect on the others. Furthermore, there may also be opportunities to select for general disease resistance in terms of metabolic stability. Although there is inconsistency in published genetic correlation estimates between milk yield and metabolic traits, selection for milk yield may be expected to lead to a deterioration in metabolic disorders. Under-recording and difficulty in diagnosing subclinical cases are among the reasons why interest is growing in using easily measurable predictors of metabolic diseases, either recorded on-farm by using sensors and milk tests or off-farm using data collected from routine milk recording. Some countries have already initiated genetic evaluations of metabolic disease traits and currently most of these use clinical observations of disease. However, there are opportunities to use clinical diseases in addition to predictor traits and genomic information to strengthen genetic evaluations for metabolic health in the future.