Phenotypic and Genetic Analyses of Mastitis, Endometritis, and Ketosis on Milk Production and Reproduction Traits in Chinese Holstein Cattle

Mastitis (MAS), endometritis (MET), and ketosis (KET) are prevalent diseases in dairy cows that result in substantial economic losses for the dairy farming industry. This study gathered 26,014 records of the health and sickness of dairy cows and 99,102 data of reproduction from 13 Holstein dairy farms in Central China; the milk protein and milk fat content from 56,640 milk samples, as well as the pedigree data of 37,836 dairy cows were obtained. The logistic regression method was used to analyze the variations in the prevalence rates of MAS, MET, and KET among various parities; the mixed linear model was used to examine the effects of the three diseases on milk production, milk quality, and reproductive traits. DMU software (version 5.2) utilized the DMUAI module in conjunction with the single-trait and two-trait animal model, as well as best linear unbiased prediction (BLUP), to estimate the genetic parameters for the three diseases, milk production, milk quality, and reproductive traits in dairy cows. The primary findings of the investigation comprised the following: (1) The prevalence rates of MAS, MET, and KET in dairy farms were 20.04%, 10.68%, and 7.33%, respectively. (2) MAS and MET had a substantial impact (p < 0.01) on milk production, resulting in significant decreases of 112 kg and 372 kg in 305-d Milk Yield (305-d MY), 4 kg and 12 kg in 305-d Protein Yield (305-d PY), and 6 kg and 16 kg in 305-d Fat Yield (305-d FY). As a result of their excessive 305-d MY, some cows were diagnosed with KET due to glucose metabolism disorder. The 305-d MY of cows with KET was significantly higher than that of healthy cows (205 kg, p < 0.01). (3) All three diseases resulted in an increase in the Interval from Calving to First Service (CTFS, 0.60-1.50 d), Interval from First Service to Conception (FSTC, 0.20-16.20 d), Calving Interval (CI, 4.00-7.00 d), and Number of Services (NUMS, 0.07-0.35). (4) The heritabilities of cows with MAS, MET, and KET were found to be low, with values of 0.09, 0.01, and 0.02, respectively. The genetic correlation between these traits ranged from 0.14 to 0.44. This study offers valuable insights on the prevention and control of the three diseases, as well as feeding management and genetic breeding.

Serum Metabolomics and Ionomics Analysis of Hoof-Deformed Cows Based on LC-MS/MS and ICP-OES/MS

In order to explore the metabolic and ionic changes of hoof-deformed cows, the serum samples of 10 healthy cows (group C) and 10 hoof-deformed cows (group T) were analyzed by LC-MS/MS and ICP-OES/MS. The pathway enrichment of differential metabolites was analyzed by screening and identifying differential metabolites and ions and using a bioinformatics method. The integration of metabolomics and ionics was analyzed with ggplot2 software in R language, and verified by MRM target metabolomics. The results showed that 127 metabolites were screened by metabolomics, of which 81 were up-regulated (p < 0.05) and 46 were down-regulated (p < 0.05). The results of ICP-OES/MS showed that 13 kinds of ions such as K, Li, and Pb in serum of dairy cows were up-regulated, while 18 kinds of ions such as Al, Cu and Sb were down-regulated. The integrated analysis of metabolomics and ionics found that potassium ions were positively correlated with L-tyrosine, L-proline, thiamine and L-valine. Sodium ions were positively correlated with L-valine and negatively correlated with α-D-glucose. The results of high-throughput target metabolomics showed that the contents of L-proline, L-phenylalanine and L-tryptophan in serum of dairy cows increased significantly, which was consistent with the results of non-target metabolomics. In a word, the metabolism and ion changes in dairy cows with hoof deformation were revealed by metabolomics and ionics.

Understanding the phenotypic and genetic background of the lactose content in Sarda dairy sheep

Lactose, the principal carbohydrate found in milk, plays an important role in the physiological processes of milk production because it is related to milk volume, and it is responsible for the osmotic equilibrium between blood and milk in the mammary gland. In this study, factors affecting lactose content (LC) in sheep milk are investigated. For this purpose, 2,358 test-day records were sampled from 509 ewes (3-7 records per animal). The LC and other main milk traits were analyzed using a mixed linear model that included days in milk (DIM) class, parity, lambing month, and type of lambing as fixed effects and animal, permanent environment, and flock test day as random effects. The pedigree-based approach was used to estimate the heritability and repeatability of LC. Moreover, the genomic background of LC was investigated through a GWAS. The LC was affected by all tested factors (i.e., DIM class, parity, lambing month, and type of lambing). Low heritability (0.10 ± 0.05) and moderate repeatability (0.42 ± 0.02) were estimated for LC. High negative genetic correlations were estimated between LC and NaCl (-0.99 ± 0.01) and between LC and somatic cell count (-0.94 ± 0.05). Only 2 markers passed the chromosome-wide Bonferroni threshold. Results of the present study, although obtained on a relatively small sample, suggest the possibility to include LC in the breeding programs, particularly because of its strong relationship with NaCl and somatic cell count.

Association of herd hyperketolactia prevalence with transition management practices and herd productivity on Canadian dairy farms-A retrospective cross-sectional study

The objective of this observational study was to assess the relationship between herd-level prevalence of hyperketolactia (HPH) with management practices of the transition period and herd milk production. Dairy herds (n = 71) were selected based on their inclusion in a herd management risk assessment study (August 2014-March 2018) using a Vital 90 (Elanco) Risk Assessment tool (one assessment per farm). Data from multiple milk recording test-days (Dairy Herd Improvement, DHI; Lactanet) were included in the analysis. Tests performed within ±6 mo relative to each farm's risk assessment date were included (10 ± 2 SD tests per farm). The majority of the farms were located in Ontario (83%). For each farm DHI test, the data set included herd average milk yield (kg/cow per day), average milk fat and protein (%), average somatic cell count (cells/mL), average days in milk (DIM), number of cows tested for ketosis, number of ketosis-positive tests (milk β-hydroxybutyrate ≥0.15 mmol/L), and proportion of cows by parity groups. Overall HPH (5-21 DIM) was calculated based on data available per farm (sum of all positive tests within 5-21 DIM/sum of all cows tested within 5-21 DIM). Each farm average was obtained by considering all test-days. A logit-transformation was applied to hyperketolactia prevalence. Linear regression models (PROC GLM and MIXED of SAS, Version 9.4) were used to predict herd HPH (milk β-hydroxybutyrate ≥0.15 mmol/L within 5 to 21 DIM; the outcome of interest). Four initial models (far-off, close-up, and fresh periods, and DHI) were separately built to assess associations between their variables and HPH; a final model considered variables selected in the initial models. Univariable (liberal P < 0.25) followed by multivariable models were used to build specific models for each period of the risk assessment. Herd prevalence of hyperketolactia was 27 ± 14%, with an average herd size of 141 ± 110 cows. The final HPH model (R² = 24.8%) included weighted milk yield, the proportion of primiparous cows, water access in the close-up period, and access to rest areas or stall access in the fresh period. Herd prevalence of hyperketolactia was negatively associated with milk yield [odds ratio, OR = 0.96 (95% confidence interval 0.92-0.99)] and proportion of primiparous cows [OR = 0.98 (0.96-0.99)]. The odds of hyperketolactia were greater with poor water access and quality (<5 cm of linear access per cow; dirty water; only 1 water location in pen) than with ≥10.2 cm of linear access per cow; clean water; >2 water locations in pen [1.23 (1.11-2.39)] in the close-up period. The odds of hyperketolactia were greater in farms providing limited access to rest areas in the fresh period than in farms providing constant access to rest areas, without dead-ends [1.64 (1.03-2.80)]. In Canadian dairy herds, HPH in early lactation was associated with certain transition-period management practices and was negatively associated with herd productivity.

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.

Characteristics of Holstein cows predisposed to ketosis during the post-partum transition period

BACKGROUND

Ketosis is a common metabolic disorder during the post-partum transition period of dairy cattle. How the method of reproduction, parturition time, and calf birth weight affect the occurrence of ketosis on dairy herds remains elusive.

OBJECTIVES

This study investigated factors associated with the severity of ketosis.

METHODS

We divided 186 Holstein cows into three classifications based on the highest β-hydroxybutyrate (BHBA) concentration during the post-partum transition period, namely non-ketosis (<1.2 mmol/L, n = 94), subclinical ketosis (1.2-2.9 mmol/L, n = 58), and clinical ketosis (≥3.0 mmol/L, n = 34). We evaluated characteristics of cows associated with the severity of ketosis.

RESULTS

Ketosis was not associated with the method of reproduction, parturition time, pregnancy wastage, premature delivery, retained placenta, and type of calf. Cows calving in spring and especially summer were at higher risk of severe ketosis (p < 0.01). Cows with increased body condition score (BCS) at parturition, age, lactation number, and calving interval were more likely to develop severe ketosis (p < 0.05). Cows with clinical ketosis produced most milk (29.9 ± 1.0 kg) from days four to six, whereas cows without ketosis produced the least (21.3 ± 0.8 kg) (p < 0.001). Heavier calf birth weight resulted in high risk of severe ketosis (p < 0.01), due to increased milk yield during the early lactation.

CONCLUSIONS

The severity of ketosis is associated with the calving season, BCS at parturition, age, lactation number, calving interval, milk yield in the early lactation period, and calf birth weight. Nonetheless, it was not associated with the method of reproduction, parturition time, pregnancy wastage, premature delivery, retained placenta, and type of calf. This study is the first to investigate the associations between ketosis and calf birth weight. Our findings could help predict cows at risk of ketosis and take precautions.

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.

Investigation of β-hydroxybutyrate in early lactation of Simmental cows: Genetic parameters and genomic predictions

Genomic information allows for a more accurate calculation of relationships among animals than the pedigree information, leading to an increase in accuracy of breeding values. Here, we used pedigree-based and single-step genomic approaches to estimate variance components and breeding values for β-hydroxybutyrate milk content (BHB). Additionally, we performed a genome-wide association study (GWAS) to depict its genetic architecture. BHB concentrations within the first 90 days of lactation, estimated from milk medium infrared spectra, were available for 30,461 cows (70,984 records). Genotypes at 42,152 loci were available for 9,123 animals. Low heritabilities were found for BHB using pedigree-based (0.09 ± 0.01) and genomic (0.10 ± 0.01) approaches. Genetic correlation between BHB and milk traits ranged from -0.27 ± 0.06 (BHB and protein percentage) to 0.13 ± 0.07 (BHB and fat-to-protein ratio) using pedigree and from -0.26 ± 0.05 (BHB and protein percentage) to 0.13 ± 0.06 (BHB and fat-to-protein ratio) using genomics. Breeding values were validated for 344 genotyped cows using linear regression method. The genomic EBV (GEBV) had greater accuracy (0.51 vs. 0.45) and regression coefficient (0.98 vs. 0.95) compared to EBV. The correlation between two subsequent evaluations, without and with phenotypes for validation cows, was 0.85 for GEBV and 0.82 for EBV. Predictive ability (correlation between (G)EBV and adjusted phenotypes) was greater when genomic information was used (0.38) than in the pedigree-based approach (0.31). Validation statistics in the pairwise two-trait models (milk yield, fat and protein percentage, urea, fat/protein ratio, lactose and logarithmic transformation of somatic cells count) were very similar to the ones highlighted for the single-trait model. The GWAS allowed discovering four significant markers located on BTA20 (57.5-58.2 Mb), where the ANKH gene is mapped. This gene has been associated with lactose, alpha-lactalbumin and BHB. Results of this study confirmed the usefulness of genomic information to provide more accurate variance components and breeding values, and important insights about the genomic determination of BHB milk content.

GC/MS and LC/MS Based Serum Metabolomic Analysis of Dairy Cows With Ovarian Inactivity

Metabolic disorders may lead to the inactive ovaries of dairy cows during early lactation. However, the detailed metabolic profile of dairy cows with inactive ovaries around 55 days postpartum has not been clearly elucidated. The objective of this study was to investigate the metabolic difference in cows with inactive ovaries and estrus from the perspective of serum metabolites. According to clinical manifestations, B-ultrasound scan, rectal examination, 15 cows were assigned to the estrus group (E; follicular diameter 15–20 mm) and 15 to the inactive ovary group (IO; follicular diameter <8 mm and increased <2 mm within 5 days over two examinations). The blood was collected from the tail vein of the cow to separate serum 55–60 days postpartum, and then milked and fasted in the morning. Serum samples were analyzed using gas chromatography time-of-flight mass spectrometry technology (GC-TOF-MS) and ultra-high-pressure liquid chromatography-quadrupole-time-of-flight mass spectrometry (UHPLC-QTOF-MS). Differences in serum metabolites were identified using multivariate statistical analysis and univariate analysis. Thirty differentially abundant metabolites were identified between the two groups. In cows with inactive ovaries compared with cows in estrus, 20 serum metabolites were significantly higher (beta-cryptoxanthin (p = 0.0012), 9-cis-retinal (p = 0.0030), oxamic acid (p = 0.0321), etc.) while 10 metabolites were significantly lower (monostearin (p = 0.0001), 3-hydroxypropionic acid (p = 0.0005), D-talose (p = 0.0018), etc.). Pathway analysis indicated that the serum differential metabolites of multiparous cows in estrus obtained by the two metabolomics techniques were mainly involved in β-alanine metabolism and steroid biosynthesis metabolism, while other involved metabolic pathways were related to metabolism of glyoxylate; dicarboxylate metabolism; fructose, mannose, glutathione, glycerolipid, glycine, serine, threonine, propanoate, retinol, and pyrimidine metabolism. This indicates that the abnormalities in glucose metabolism, lipid metabolism, amino acid metabolism, and glutathione metabolism of postpartum dairy cows obstructed follicular development.

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.

Invited review: The influence of immune activation on transition cow health and performance-A critical evaluation of traditional dogmas

The progression from gestation into lactation represents the transition period, and it is accompanied by marked physiological, metabolic, and inflammatory adjustments. The entire lactation and a cow's opportunity to have an additional lactation are heavily dependent on how successfully she adapts during the periparturient period. Additionally, a disproportionate amount of health care and culling occurs early following parturition. Thus, lactation maladaptation has been a heavily researched area of dairy science for more than 50 yr. It was traditionally thought that excessive adipose tissue mobilization in large part dictated transition period success. Further, the magnitude of hypocalcemia has also been assumed to partly control whether a cow effectively navigates the first few months of lactation. The canon became that adipose tissue released nonesterified fatty acids (NEFA) and the resulting hepatic-derived ketones coupled with hypocalcemia lead to immune suppression, which is responsible for transition disorders (e.g., mastitis, metritis, retained placenta, poor fertility). In other words, the dogma evolved that these metabolites and hypocalcemia were causal to transition cow problems and that large efforts should be enlisted to prevent increased NEFA, hyperketonemia, and subclinical hypocalcemia. However, despite intensive academic and industry focus, the periparturient period remains a large hurdle to animal welfare, farm profitability, and dairy sustainability. Thus, it stands to reason that there are alternative explanations to periparturient failures. Recently, it has become firmly established that immune activation and the ipso facto inflammatory response are a normal component of transition cow biology. The origin of immune activation likely stems from the mammary gland, tissue trauma during parturition, and the gastrointestinal tract. If inflammation becomes pathological, it reduces feed intake and causes hypocalcemia. Our tenet is that immune system utilization of glucose and its induction of hypophagia are responsible for the extensive increase in NEFA and ketones, and this explains why they (and the severity of hypocalcemia) are correlated with poor health, production, and reproduction outcomes. In this review, we argue that changes in circulating NEFA, ketones, and calcium are simply reflective of either (1) normal homeorhetic adjustments that healthy, high-producing cows use to prioritize milk synthesis or (2) the consequence of immune activation and its sequelae.

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.

Integration of Wet-Lab Measures, Milk Infrared Spectra, and Genomics to Improve Difficult-to-Measure Traits in Dairy Cattle Populations

The objective of this study was to evaluate the contribution of Fourier-transformed infrared spectroscopy (FTIR) data for dairy cattle breeding through two different approaches: (i) estimating the genetic parameters for 30 measured milk traits and their FTIR predictions and investigating the additive genetic correlation between them and (ii) evaluating the effectiveness of FTIR-derived phenotyping to replicate a candidate bull’s progeny testing or breeding value prediction at birth. Records were available from 1,123 cows phenotyped using gold standard laboratory methodologies (LAB data). This included phenotypes related to fine milk composition and milk technological characteristics, milk acidity, and milk protein fractions. The dataset used to generate FTIR predictions comprised 729,202 test-day records from 51,059 Brown Swiss cows (FIELD data). A first approach consisted of estimating genetic parameters for phenotypes available from LAB and FIELD datasets. To do so, a set of bivariate animal models were run, and genetic correlations between LAB and FIELD phenotypes were estimated using FIELD information obtained at the population level. Heritability estimates were generally higher for FIELD predictions than for the corresponding LAB measures. The additive genetic correlations (r_a) between LAB and FIELD phenotypes had different magnitudes across traits but were generally strong. Overall, these results demonstrated the potential of using FIELD information as indicator traits for the indirect genetic improvement of LAB measures. In the second approach, we included genotype information for 1,011 cows from the LAB dataset, 1,493 cows from the FIELD dataset, 181 sires with daughters in both LAB and FIELD datasets, and 540 sires with daughters in the FIELD dataset only. Predictions were obtained using the single-step GBLUP method. A four fold cross-validation was used to assess the predictive ability of the different models, assessed as the ability to predict masked LAB records from daughters of progeny testing bulls. The correlation between observed and predicted LAB measures in validation was averaged over the four training-validation sets. Different sets of phenotypic information were used sequentially in cross-validation schemes: (i) LAB cows from the training set; (ii) FIELD cows from the training set; and (iii) FIELD cows from the validation set. Models that included FIELD records showed an improvement for the majority of traits. This study suggests that breeding programs for difficult-to-measure traits could be implemented using FTIR information. While these programs should use progeny testing, acceptable values of accuracy can be achieved also for bulls without phenotyped progeny. Robust calibration equations are, deemed as essential.

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.

Large-scale phenotyping in dairy sector using milk MIR spectra: Key factors affecting the quality of predictions

Methods and technologies enabling the estimation at large scale of important traits for the dairy sector are of great interest. Those phenotypes are necessary to improve herd management, animal genetic evaluation, and milk quality control. In the recent years, the research was very active to predict new phenotypes from the mid-infrared (MIR) analysis of milk. Models were developed to predict phenotypes such as fine milk composition, milk technological properties or traits related to cow health, fertility and environmental impact. Most of models were developed within research contexts and often not designed for routine use. The implementation of models at a large scale to predict new traits of interest brings new challenges as the factors influencing the robustness of models are poorly documented. The first objective of this work is to highlight the impact on prediction accuracy of factors such as the variability of the spectral and reference data, the spectral regions used and the complexity of models. The second objective is to emphasize methods and indicators to evaluate the quality of models and the quality of predictions generated under routine conditions. The last objective is to outline the issues and the solutions linked with the use and transfer of models on large number of instruments. Based on partial least square regression and 10 datasets including milk MIR spectra and reference quantitative values for 57 traits of interest, the impact of the different factors is illustrated by evaluating the influence on the validation root mean square error of prediction (RMSEP). In the displayed examples, all factors, when well set up, increase the quality of predictions, with an improvement of the RMSEP ranging from 12% to 43%. This work also aims to underline the need for and the complementarity between different validation procedures, statistical parameters and quality assurance methods. Finally, when using and transferring models, the impact of the spectral standardization on the prediction reproducibility is highlighted with an improvement up to 86% with the tested models, and the monitoring of individual spectrometer stability over time appears essential. This list inspired from our experience is of course not exhaustive. The displayed results are only examples and not general rules and other aspects play a role in the quality of final predictions. However, this work highlights good practices, methods and indicators to increase and evaluate quality of phenotypes predicted at a large scale. The results obtained argue for the development of guidelines at international levels, as well as international collaborations in order to constitute large and robust datasets and enable the use of models in routine conditions.

Heritability estimates of predicted blood β-hydroxybutyrate and nonesterified fatty acids and relationships with milk traits in early-lactation Holstein cows

At the beginning of lactation, high-producing cows commonly experience an unbalanced energy status that is often responsible for the onset of metabolic disorders and impaired health and performance. Blood β-hydroxybutyrate (BHB) and nonesterified fatty acids (NEFA) are indicators of excessive fat mobilization and circulating ketone bodies. Recently, prediction models based on mid-infrared (MIR) spectroscopy have been developed to assess blood BHB and NEFA from routinely collected individual milk samples. This study aimed to estimate genetic parameters of blood BHB and NEFA predicted from milk MIR spectra and to assess their phenotypic and genetic correlations with milk production and composition traits in early-lactation Holstein cows. The data set comprised the first test-day record within lactation and spectra of individual milk samples (n = 22,718) of 13,106 Holstein cows collected from 5 to 35 d in milk (DIM). Blood BHB and NEFA were predicted from milk MIR spectra using previously developed prediction models. Genetic parameters of blood metabolites and milk traits were estimated for the whole observational period (5-35 DIM) and within 6 classes of DIM. Blood BHB and NEFA showed similar genetic variation across DIM, with the highest heritability in the first 10 d after calving (0.31 ± 0.06 and 0.19 ± 0.05 for BHB and NEFA, respectively). The genetic correlation between BHB and NEFA was moderate (0.51 ± 0.05). Genetic correlations of BHB with milk yield, SCS, protein percentage, lactose percentage, and urea nitrogen content were similar to, or at least in the same direction as, the correlations of NEFA with the same traits, whereas opposite correlations were observed with fat percentage and fat-to-protein ratio. Results of the current study suggest that blood BHB and NEFA predicted from milk MIR spectra have genetic variation that is potentially exploitable for breeding purposes. Therefore, they could be used as indicator traits of hyperketonemia in a selection index aimed to reduce the susceptibility of dairy cows to metabolic disorders in early lactation.

The evolving role of Fourier-transform mid-infrared spectroscopy in genetic improvement of dairy cattle

Over the last 100 years, significant advances have been made in the characterisation of milk composition for dairy cattle improvement programs. Technological progress has enabled a shift from labour intensive, on-farm collection and processing of samples that assess yield and fat levels in milk, to large-scale processing of samples through centralised laboratories, with the scope extended to include quantification of other traits. Fourier-transform mid-infrared (FT-MIR) spectroscopy has had a significant role in the transformation of milk composition phenotyping, with spectral-based predictions of major milk components already being widely used in milk payment and animal evaluation systems globally. Increasingly, there is interest in analysing the individual FT-MIR wavenumbers, and in utilising the FT-MIR data to predict other novel traits of importance to breeding programs. This includes traits related to the nutritional value of milk, the processability of milk into products such as cheese, and traits relevant to animal health and the environment. The ability to successfully incorporate these traits into breeding programs is dependent on the heritability of the FT-MIR predicted traits, and the genetic correlations between the FT-MIR predicted and actual trait values. Linking FT-MIR predicted traits to the underlying mutations responsible for their variation can be difficult because the phenotypic expression of these traits are a function of a diverse range of molecular and biological mechanisms that can obscure their genetic basis. The individual FT-MIR wavenumbers give insights into the chemical composition of milk and provide an additional layer of granularity that may assist with establishing causal links between the genome and observed phenotypes. Additionally, there are other molecular phenotypes such as those related to the metabolome, chromatin accessibility, and RNA editing that could improve our understanding of the underlying biological systems controlling traits of interest. Here we review topics of importance to phenotyping and genetic applications of FT-MIR spectra datasets, and discuss opportunities for consolidating FT-MIR datasets with other genomic and molecular data sources to improve future dairy cattle breeding programs.

Fasciola hepatica seroprevalence in Northern German dairy herds and associations with milk production parameters and milk ketone bodies

Infections with the liver fluke Fasciola hepatica remain a serious problem in dairy herds causing significant production losses. In sheep, a strong relationship between F. hepatica infections and an increase in serum ketone bodies due to reduced feed intake and liver damage was demonstrated. We hypothesized that F. hepatica infections might contribute to an increase in milk ketone bodies in dairy herds. Thus, the objective of the study was to estimate the association between F. hepatica bulk tank milk (BTM) antibodies and milk production parameters (milk yield, milk protein, fat yield), somatic cell count (SCC) and the milk ketone bodies ß-hydroxybutyrate (BHB) and acetone, inferred from Fourier transform infrared (FTIR) spectrometry, via linear mixed model analysis. A further aim was to follow up the F. hepatica seroprevalence in dairy herds in the northern German region East Frisia. We collected BTM samples between October and December from 1022 herds in 2017 and 1318 herds in 2018. Overall, 33.1 % of the herds tested positive in 2017 and 37.0 % in 2018, showing decreased F. hepatica seroprevalences compared to prior seroprevalence studies in the same region in 2010, 2008 and 2006 (> 45 % positive herds). We estimated a significant negative association (P < 0.001) between herd F. hepatica infection category and average milk yield with a loss of -1.62 kg per cow per day in strongly infected herds compared to BTM ELISA negative herds. Moreover, F. hepatica infection category had a significant effect on herd average milk protein and fat yield (P < 0.001), showing a decrease of 0.06 kg for both parameters from BTM ELISA negative herds to strongly infected herds. No significant association with milk SCC was found (P = 0.664). Regarding ketone bodies, we estimated significant higher average BHB values in strongly infected herds compared to the other three infection categories in the model analysis (P = 0.002). The association between F. hepatica infection category and acetone values was not significant (P = 0.079). Besides primary ketosis, fasciolosis should be considered as differential diagnosis in dairy herds with increased BHB values.

Blood β-hydroxybutyrate concentrations and early lactation management strategies on pasture-based dairy farms in Colombia

The increasing global demand for food requires sustainable solutions to close the gap in agricultural yield between industrialized and non-industrialized countries. Our objectives in this cross-sectional study were to: 1) characterize farm populations, milk yield, and early lactation management strategies of dairy cows in three different regions of Colombia, and 2) determine the association of these management strategies with blood β-hydroxybutyrate (BHB) concentrations in the first 42 days in milk (DIM). Dairy herds (n = 56) in the Antioquia, Caldas, and Cundinamarca regions of Colombia were visited once from May through July 2018. A survey was administered to farm owners to collect demographic, management, and herd nutrition information. Blood samples from dairy cows (n = 880) between calving and 42 DIM were used to measure blood BHB concentration. Associations between management and nutritional strategies and blood BHB concentration were examined using mixed models. Prevalence of hyperketonemia was calculated as the number of samples with BHB concentration ≥1.2 mmol/L divided by the total number of samples. The estimated diet composition for early lactation dairy cows was 65.5% pasture and 31.8% commercial concentrates. The farm median milk yield, protein concentration, and fat concentration were 21.0 kg (range = 13.1-36 kg), 3.2% (range = 2.7-4.1%), and 3.5% (range = 3.0-4.1%), respectively. Milk yield least squares means (95% confidence interval; CI) differed by region: 21.7 (20.3, 23.2), 18.5 (17.0, 20.2), and 20.3 (18.5, 22.4) kg in Antioquia, Caldas, and Cundinamarca, respectively. Median blood BHB concentration was 0.5 and ranged from 0.1-4.4 mmol/L; blood BHB concentration was not different among the three regions. Pasture fertilization, increased parity, and BCS were associated with changes in blood BHB concentration. The overall prevalence of hyperketonemia was 4.5%. Geographical region affected the prevalence of hyperketonemia at 2.5%, 4.0%, and 10.2% in Antioquia, Caldas, and Cundinamarca, respectively. Mean stocking density (95% CI) was greater in Cundinamarca than Antioquia or Caldas at 3.3 (2.2, 5.0), 2.8 (2.1, 3.9) and 1.7 (1.2, 2.6) animals per ha, respectively, and was associated with hyperketonemia prevalence. Farms that abruptly stop milking cows at dry-off had 80% of the hyperketonemia events in the study. Pasture-based dairies in Colombia had lower blood BHB concentrations and estimated milk yield compared with confined production systems in temperate zones. However, geographical region, stocking density, and abrupt cessation of milking at dry-off were associated with prevalence of hyperketonemia in pasture-based dairies.

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.

Invited review: Milk lactose-Current status and future challenges in dairy cattle

Lactose is the main carbohydrate in mammals' milk, and it is responsible for the osmotic equilibrium between blood and alveolar lumen in the mammary gland. It is the major bovine milk solid, and its synthesis and concentration in milk are affected mainly by udder health and the cow's energy balance and metabolism. Because this milk compound is related to several biological and physiological factors, information on milk lactose in the literature varies from chemical properties to heritability and genetic associations with health traits that may be exploited for breeding purposes. Moreover, lactose contributes to the energy value of milk and is an important ingredient for the food and pharmaceutical industries. Despite this, lactose has seldom been included in milk payment systems, and it has never been used as an indicator trait in selection indices. The interest in lactose has increased in recent years, and a summary of existing information about lactose in the dairy sector would be beneficial for the scientific community and the dairy industry. The present review collects and summarizes knowledge about lactose by covering and linking several aspects of this trait in bovine milk. Finally, perspectives on the use of milk lactose in dairy cattle, especially for selection purposes, are outlined.

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.

Invited review: β-hydroxybutyrate concentration in blood and milk and its associations with cow performance

Hyperketonemia (HYK) is one of the most frequent and costly metabolic disorders in high-producing dairy cows and its diagnosis is based on β-hydroxybutyrate (BHB) concentration in blood. In the last 10 years, the number of papers that have dealt with the impact of elevated BHB levels in dairy cattle has increased. Therefore, this paper reviewed the recent literature on BHB concentration in blood and milk, and its relationships with dairy cow health and performance, and farm profitability. Most studies applied the threshold of 1.2 mmol/l of BHB concentration in blood to indicate HYK; several authors considered BHB concentrations between 1.2 and 2.9 mmol/l as subclinical ketosis, and values ⩾3.0 mmol/l as clinical ketosis. Results on HYK frequency (prevalence and incidence) and cow performance varied according to parity and days in milk, being greater in multiparous than in primiparous cows, and in the first 2 weeks of lactation than in later stages. Hyperketonemia has been associated with greater milk fat content, fat-to-protein ratio and energy-corrected milk, and lower protein and urea nitrogen in milk. The relationships with milk yield and somatic cell count are still controversial. In general, HYK impairs health of dairy cows by increasing the risk of the onset of other early lactation diseases, and it negatively affects reproductive performance. The economic cost of HYK is mainly due to impaired reproductive performance and milk loss. From a genetic point of view, results from the literature suggested the feasibility of selecting cows with low susceptibility to HYK. The present review highlights that milk is the most promising matrix to identify HYK, because it is easy to sample and allows a complete screening of the herd through BHB concentration predicted using mid-IR spectroscopy during routine milk recording. Further research is needed to validate accurate and convenient methods to discriminate between cows in risk of HYK and healthy animals in field conditions and to support farmers to achieve an early detection and minimise the economic losses.

A landscape of the heritability of Fourier-transform infrared spectral wavelengths of milk samples by parity and lactation stage in Holstein cows

Fourier-transform near- and mid-infrared (FTIR) milk spectral data are routinely collected in many countries worldwide. Establishing an optimal strategy to use spectral data in genetic evaluations requires knowledge of the heritabilities of individual FTIR wavelength absorbances. Previous FTIR heritability estimates have been based on relatively small sample sizes and have not considered the possibility that heritability may vary across parities and stages of the lactation. We used data from ∼370,000 test-day records of Canadian Holstein cows to produce a landscape of the heritability of FTIR spectra, 1,060 wavelengths in the near- and mid-infrared spectrum (5,011-925 cm^-1), by parity and month of the lactation (mo 1 to 3 and mo 1 to 6, respectively). The 2 regions of the spectrum associated with absorption of electromagnetic energy by water molecules were estimated to have very high phenotypic variances, very low heritabilities, and very low proportion of variance explained by herd-year-season (HYS) subclasses. The near- or short-wavelength infrared (SWIR: 5,066-3,672 cm^-1) region was also characterized by low heritability estimates, whereas the estimated proportion of the variance explained by HYS was high. The mid-wavelength infrared region (MWIR: 3,000-2,500 cm^-1) and the transition between mid and long-wavelength infrared region (MWIR-LWIR: 1,500-925 cm^-1) harbor several waves characterized by moderately high (≥0.4) heritabilities. Most of the high-heritability regions contained wavelengths that are reported to be associated with important milk metabolites and components. Interestingly, these 2 same regions tended to show more variability in heritabilities between parity and lactation stage. Second parity showed heritability patterns that were distinctly different from those of the first and third parities, whereas the first 2 mo of the lactation had clearly distinct heritability patterns compared with mo 3 to 6.

Genetic variation in serum protein pattern and blood β-hydroxybutyrate and their relationships with udder health traits, protein profile, and cheese-making properties in Holstein cows

The aim of this study was to investigate in Holstein cows the genetic basis of blood serum metabolites [i.e., total protein, albumin, globulin, albumin:globulin ratio (A:G), and blood β-hydroxybutyrate (BHB)], a set of milk phenotypes related to udder health, milk quality technological characteristics, and genetic relationships among them. Samples of milk were collected from 498 Holstein cows belonging to 28 herds. All animal welfare and milk phenotypes were assessed using standard analytical methodology. A set of Bayesian univariate and bivariate animal models was implemented via Gibbs sampling, and statistical inference was based on the marginal posterior distributions of parameters of concern. We observed a small additive genetic influence for serum albumin concentrations, moderate heritability (≥0.20) for total proteins, globulins, and A:G, and high heritability (0.37) for blood BHB. Udder health traits (somatic cell score, milk lactose, and milk pH) showed low or moderate heritabilities (0.15-0.20), whereas variations in milk protein fraction concentrations were confirmed as mostly under genetic control (heritability: 0.21-0.71). The moderate and high heritabilities observed for milk coagulation properties and curd firming modeling parameters provided confirmation that genetic background exerts a strong influence on the cheese-making ability of milk, largely due to genetic polymorphisms in the major milk protein genes. Blood BHB showed strong negative genetic correlations with globulins (-0.619) but positive correlations with serum albumin (0.629) and A:G (0.717), which suggests that alterations in the serum protein pattern and BHB blood levels are likely to be genetically related. Strong relationships were found between albumin and fat percentages (-0.894), between globulin and α_S2-CN (-0.610), and, to a lesser extent, between serum protein pattern and milk technological characteristics. Genetic relationships between blood BHB and traits related to udder health and milk quality and technological characteristics were mostly weak. This study provides evidence that there is exploitable additive genetic variation for traits related to animal health and welfare and throws light on the shared genetic basis of these traits and the phenotypes related to the quality and cheese-making ability of milk.