Insights on early response to acute heat shock of bovine mammary epithelial cells through a multimethod approach

Heat stress is a significant challenge in dairy cattle herds, affecting milk production and quality, and generating important changes at the cellular level. Most in vitro research on heat shock (HS) effects on dairy cow mammary cells was focused on medium-long-term effects. In recent years, Fourier transform-infrared (FT-IR) micro-spectroscopy has been increasingly used to study the effects of several external stresses on different cell lines, down to the level of single cellular components, such as DNA/RNA, lipids, and proteins. In this study, the possible changes at the biochemical and molecular level induced by acute (30 min-2 h) HS in bovine mammary epithelial (BME-UV1) cells were investigated. The cells were exposed to different temperatures, thermoneutral (TN, 37 °C) and HS (42 °C), and FT-IR spectra were acquired to analyse the effects of HS on biochemical characteristics of BME-UV1 cellular components (proteins, lipids, and DNA/RNA). Moreover, cell viability assay, reactive oxygen species production, and mRNA expression of heat shock proteins (HSPA1A, HSP90AA1, GRP78, GRP94) and antioxidant genes (SOD1, SOD2) by RT-qPCR were also analysed. The FT-IR results showed a change already at 30 min of HS exposure, in the content of long-chain fatty acids, which probably acted as a response to a modification of membrane fluidity in HS cells compared with TN cells. After 2 h of HS exposure, modification of DNA/RNA activity and accumulation of aggregated proteins was highlighted in HS cells. The gene expression analyses showed the overexpression of HSPA1A and HSP90AA1 starting from 30 min up to 2 h in HS cells compared with TN cells. At 2 h of HS exposure, also the overexpression of GRP94 was observed in HS cells. Acute HS did not affect cell viability, reactive oxygen species level, and SOD1 and SOD2 gene expression of BME-UV1 cells. According to the results obtained, cells initiate early defence mechanisms in case of acute HS and probably this efficient response capacity may be decisive for tolerance to heat stress of dairy cattle.

Combining genetic markers, on-farm information and infrared data for the in-line prediction of blood biomarkers of metabolic disorders in Holstein cattle

BACKGROUND

Various blood metabolites are known to be useful indicators of health status in dairy cattle, but their routine assessment is time-consuming, expensive, and stressful for the cows at the herd level. Thus, we evaluated the effectiveness of combining in-line near infrared (NIR) milk spectra with on-farm (days in milk [DIM] and parity) and genetic markers for predicting blood metabolites in Holstein cattle. Data were obtained from 388 Holstein cows from a farm with an AfiLab system. NIR spectra, on-farm information, and single nucleotide polymorphisms (SNP) markers were blended to develop calibration equations for blood metabolites using the elastic net (ENet) approach, considering 3 models: (1) Model 1 (M1) including only NIR information, (2) Model 2 (M2) with both NIR and on-farm information, and (3) Model 3 (M3) combining NIR, on-farm and genomic information. Dimension reduction was considered for M3 by preselecting SNP markers from genome-wide association study (GWAS) results.

RESULTS

Results indicate that M2 improved the predictive ability by an average of 19% for energy-related metabolites (glucose, cholesterol, NEFA, BHB, urea, and creatinine), 20% for liver function/hepatic damage, 7% for inflammation/innate immunity, 24% for oxidative stress metabolites, and 23% for minerals compared to M1. Meanwhile, M3 further enhanced the predictive ability by 34% for energy-related metabolites, 32% for liver function/hepatic damage, 22% for inflammation/innate immunity, 42.1% for oxidative stress metabolites, and 41% for minerals, compared to M1. We found improved predictive ability of M3 using selected SNP markers from GWAS results using a threshold of > 2.0 by 5% for energy-related metabolites, 9% for liver function/hepatic damage, 8% for inflammation/innate immunity, 22% for oxidative stress metabolites, and 9% for minerals. Slight reductions were observed for phosphorus (2%), ferric-reducing antioxidant power (1%), and glucose (3%). Furthermore, it was found that prediction accuracies are influenced by using more restrictive thresholds (-log10(P-value) > 2.5 and 3.0), with a lower increase in the predictive ability.

CONCLUSION

Our results highlighted the potential of combining several sources of information, such as genetic markers, on-farm information, and in-line NIR infrared data improves the predictive ability of blood metabolites in dairy cattle, representing an effective strategy for large-scale in-line health monitoring in commercial herds.

Prediction of likelihood of conception in dairy cows using milk mid-infrared spectra collected before the first insemination and machine learning algorithms

Accurate and ex-ante prediction of cows' likelihood of conception (LC) based on milk composition information could improve reproduction management on dairy farms. Milk composition is already routinely measured by mid-infrared (MIR) spectra, which are known to change with advancing stages of pregnancy. For lactating cows, MIR spectra may also be used for predicting the LC. Our objectives were to classify the LC at first insemination using milk MIR spectra data collected from calving to first insemination and to identify the spectral regions that contribute the most to the prediction of LC at first insemination. After quality control, 4,866 MIR spectra, milk production, and reproduction records from 3,451 Holstein cows were used. The classification accuracy and area under the curve (AUC) of 6 models comprising different predictors and 3 machine learning methods were estimated and compared. The results showed that partial least square discriminant analysis (PLS-DA) and random forest had higher prediction accuracies than logistic regression. The classification accuracy of good and poor LC cows and AUC in herd-by-herd validation of the best model were 76.35 ± 10.60% and 0.77 ± 0.11, respectively. All wavenumbers with values of variable importance in the projection higher than 1.00 in PLS-DA belonged to 3 spectral regions, namely from 1,003 to 1,189, 1,794 to 2,260, and 2,300 to 2,660 cm-1. In conclusion, the model can predict LC in dairy cows from a high productive TMR system before insemination with a relatively good accuracy, allowing farmers to intervene in advance or adjust the insemination schedule for cows with a poor predicted LC.

Genetic analysis of milk citrate predicted by milk mid-infrared spectra of Holstein cows in early lactation

Milk citrate is regarded as an early biomarker of negative energy balance in dairy cows during early lactation and serves as a suitable candidate phenotype for genomic selection due to its wide availability across a large number of cows through milk mid-infrared spectra prediction. However, its genetic background is not well known. Therefore, the objectives of this study were to (1) analyze the genetic parameters of milk citrate; (2) identify genomic regions associated with milk citrate; and (3) analyze the functional annotation of candidate genes and quantitative trait loci (QTL) related to milk citrate in Walloon Holstein cows. In total, 134,517 test-day milk-citrate phenotypes (mmol/L) collected within the first 50 d in milk on 52,198 Holstein cows were used. These milk-citrate phenotypes, predicted by milk mid-infrared spectra, were divided into 3 traits according to the first (citrate1), second (citrate2), and third to fifth parity (citrate3+). Genomic information for 566,170 SNPs was available for 4,479 animals. A multiple-trait repeatability model was used to estimate genetic parameters. A single-step GWAS was used to identify candidate genes for citrate and post-GWAS analysis was done to investigate the relationship and function of the identified candidate genes. The heritabilities estimated for citrate1, citrate2, and citrate3+ were 0.40, 0.37, and 0.35, respectively. The genetic correlations among the 3 traits ranged from 0.98 to 0.99. The genomic correlations among the 3 traits were also close to 1.00 across the genomic regions (1 Mb) in the whole genome, which means that citrate can be considered as a single trait in the first 5 parities. In total, 603 significant SNPs located on 3 genomic regions (chromosome 7, 68.569-68.575 Mb; chromosome 14, 0.15-1.90 Mb; and chromosome 20, 54.00-64.28 Mb), were identified to be associated with milk citrate. We identified 89 candidate genes including GPT, ANKH, PPP1R16A, and 32 QTL reported in the literature related to the identified significant SNPs. These identified QTL were mainly reported associated with milk fatty acids and metabolic diseases in dairy cows. This study suggests that milk citrate in Holstein cows is highly heritable and has the potential to be used as an early proxy for the negative energy balance of Holstein cows in a breeding objective.

Prediction of key milk biomarkers in dairy cows through milk mid-infrared spectra and international collaborations

At the individual cow level, suboptimum fertility, mastitis, negative energy balance, and ketosis are major issues in dairy farming. These problems are widespread on dairy farms and have an important economic impact. The objectives of this study were (1) to assess the potential of milk mid-infrared (MIR) spectra to predict key biomarkers of energy deficit (citrate, isocitrate, glucose-6 phosphate [glucose-6P], free glucose), ketosis (β-hydroxybutyrate [BHB] and acetone), mastitis (N-acetyl-β-d-glucosaminidase activity [NAGase] and lactate dehydrogenase), and fertility (progesterone); (2) to test alternative methodologies to partial least squares (PLS) regression to better account for the specific asymmetric distribution of the biomarkers; and (3) to create robust models by merging large datasets from 5 international or national projects. Benefiting from this international collaboration, the dataset comprised a total of 9,143 milk samples from 3,758 cows located in 589 herds across 10 countries and represented 7 breeds. The samples were analyzed by reference chemistry for biomarker contents, whereas the MIR analyses were performed on 30 instruments from different models and brands, with spectra harmonized into a common format. Four quantitative methodologies were evaluated to address the strongly skewed distribution of some biomarkers. Partial least squares regression was used as the reference basis, and compared with a random modification of distribution associated with PLS (random-downsampling-PLS), an optimized modification of distribution associated with PLS (KennardStone-downsampling-PLS), and support vector machine (SVM). When the ability of MIR to predict biomarkers was too low for quantification, different qualitative methodologies were tested to discriminate low versus high values of biomarkers. For each biomarker, 20% of the herds were randomly removed within all countries to be used as the validation dataset. The remaining 80% of herds were used as the calibration dataset. In calibration, the 3 alternative methodologies outperform the PLS performances for the majority of biomarkers. However, in the external herd validation, PLS provided the best results for isocitrate, glucose-6P, free glucose, and lactate dehydrogenase (coefficient of determination in external herd validation [R2v] = 0.48, 0.58, 0.28, and 0.24, respectively). For other molecules, PLS-random-downsampling and PLS-KennardStone-downsampling outperformed PLS in the majority of cases, but the best results were provided by SVM for citrate, BHB, acetone, NAGase, and progesterone (R2v = 0.94, 0.58, 0.76, 0.68, and 0.15, respectively). Hence, PLS and SVM based on the entire dataset provided the best results for normal and skewed distributions, respectively. Complementary to the quantitative methods, the qualitative discriminant models enabled the discrimination of high and low values for BHB, acetone, and NAGase with a global accuracy around 90%, and glucose-6P with an accuracy of 83%. In conclusion, MIR spectra of milk can enable quantitative screening of citrate as a biomarker of energy deficit and discrimination of low and high values of BHB, acetone, and NAGase, as biomarkers of ketosis and mastitis. Finally, progesterone could not be predicted with sufficient accuracy from milk MIR spectra to be further considered. Consequently, MIR spectrometry can bring valuable information regarding the occurrence of energy deficit, ketosis, and mastitis in dairy cows, which in turn have major influences on their fertility and survival.

Imputation of missing milk Fourier transform mid-infrared spectra using existing milk spectral databases: A strategy to improve the reliability of breeding values and predictive models

The use of milk Fourier transform mid-infrared (FT-MIR) spectrometry to develop management and breeding tools for dairy farmers and industry is growing and supported by the availability of numerous new predicted phenotypes to assess the nutritional quality of milk and its technological properties, but also the animal health and welfare status and its environmental fingerprint. For genetic evaluations, having a long-term and representative spectral dairy herd improvement (DHI) database improves the reliabilities of estimated breeding values (EBV) from these phenotypes. Unfortunately, most of the time, the raw spectral data used to generate these estimations are not stored. Moreover, many reference measurements of those phenotypes, needed during the FT-MIR calibration step, are available from past research activities but lack spectra records. So, it is impossible to use them to improve the FT-MIR models. Consequently, there is a strong interest in imputing those missing spectra. The innovative objective of this study was to use the existing large spectral DHI database to estimate missing spectra by selecting probable spectra using, as the match criteria, common dairy traits recorded for a long time by DHI organizations. We tested 4 match criteria combinations. Combination 1 required to have equal fat and protein contents between the sample for which a spectrum was to be estimated and the reference samples in the DHI database. Combination 2 also required an equal urea content. Combination 3 requested equal fat, protein, and lactose contents. Finally, combination 4 included all criteria. When more than one spectrum was found during the search, their average was the estimated spectrum for the query sample. Concretely, this study estimated missing spectra for 1,700 samples using 2,000,000 spectral DHI records. For assessing the effect of this spectral estimation on the prediction quality, FT-MIR equations were used to predict 11 phenotypes, selected as their quantification used different FT-MIR regions. They were related to the milk fat and mineral composition, lactoferrin content, quantity of eructed methane, body weight (BW), and dry matter intake. The accuracy between predictions obtained from actual and estimated spectra was evaluated by calculating the mean absolute error (MAE). The criteria in the fourth and second combinations were too strict to estimate a spectrum for most samples. Indeed, for many samples, no spectra with the same values for those matching criteria was found. The third match criteria combination had a poorer prediction performance for all studied traits and spectral absorptions than the first combination due to fewer matched samples available to compute the missing spectrum. By allowing a range for matching lactose content (±0.1 g/dL milk), we showed that this new combination increased the number of selected samples to compute missing spectra and predict better the infrared absorption at different wavenumbers, especially those related to the lactose quantification. The prediction performance was further improved by performing queries on the entire Walloon DHI spectral database (6,625,570 spectra), and it varied among the studied phenotypes. Without considering the traits used for the matching, the best predictions were obtained for the content of saturated fatty acids (MAE = 0.15 g/dL milk) and BW (MAE = 12.80 kg). Yet, the predictions for the unsaturated fatty acids were less accurate (MAE = 0.13 and 0.018 g/dL milk for monounsaturated and polyunsaturated fatty acids), likely because of the poorer predictions of spectral regions related to long-chain fatty acids. Similarly, poorer predictions were observed for the amount of methane eructed by dairy cows (MAE = 47.02 g/d), likely because it is not directly related to fat content or composition. Prediction accuracies for the remaining traits were also low. In conclusion, we observed that increasing the number of relevant matching criteria helps improve the quality of FT-MIR predicted phenotypes and the number of spectra used during the search. So, it would be of great interest to test in the future the suitability of the developed methodology with large-scale international spectral databases to improve the reliability of EBV from these FT-MIR-based phenotypes and the robustness of FT-MIR predictive models.

Mid-Level Data Fusion Combined with the Fingerprint Region for Classification DON Levels Defect of Fusarium Head Blight Wheat

In this study, a method of mid-level data fusion with the fingerprint region was proposed, which was combined with the characteristic wavelengths that contain fingerprint information in NIR and FT-MIR spectra to detect the DON level in FHB wheat during wheat processing. NIR and FT-MIR raw spectroscopy data on normal wheat and FHB wheat were obtained in the experiment. MSC was used for pretreatment, and characteristic wavelengths were extracted by CARS, MGS and XLW. The variables that can effectively reflect fingerprint information were retained to build the mid-level data fusion matrix. LS-SVM and PLS-DA were applied to investigate the performance of the single spectroscopic model, mid-level data fusion model and mid-level data fusion with fingerprint information model, respectively. The experimental results show that mid-level data fusion with a fingerprint information strategy based on fused NIR and FT-MIR spectra represents an effective method for the classification of DON levels in FHB wheat samples.

Effectiveness of mid-infrared spectroscopy for the prediction of cow milk metabolites

Proton nuclear magnetic resonance (¹H NMR) spectroscopy is acknowledged as one of the most powerful analytical methods with cross-cutting applications in dairy foods. To date, the use of ¹H NMR spectroscopy for the collection of milk metabolic profile is hindered by costly and time-consuming sample preparation and analysis. The present study aimed at evaluating the accuracy of mid-infrared spectroscopy (MIRS) as a rapid method for the prediction of cow milk metabolites determined through ¹H NMR spectroscopy. Bulk milk (n = 72) and individual milk samples (n = 482) were analyzed through one-dimensional ¹H NMR spectroscopy and MIRS. Nuclear magnetic resonance spectroscopy identified 35 milk metabolites, which were quantified in terms of relative abundance, and MIRS prediction models were developed on the same 35 milk metabolites, using partial least squares regression analysis. The best MIRS prediction models were developed for galactose-1-phosphate, glycerophosphocholine, orotate, choline, galactose, lecithin, glutamate, and lactose, with coefficient of determination in external validation from 0.58 to 0.85, and ratio of performance to deviation in external validation from 1.50 to 2.64. The remaining 27 metabolites were poorly predicted. This study represents a first attempt to predict milk metabolome. Further research is needed to specifically address whether developed prediction models may find practical application in the dairy sector, with particular regard to the screening of dairy cows' metabolic status, the quality control of dairy foods, and the identification of processed milk or incorrectly stored milk.

Seasonal Changes in the Antioxidant Activity and Biochemical Parameters of Goat Milk

Goats are ubiquitous, including in hot and dry regions, while also being very sensitive to climate fluctuations, expressed in temperature differences. This affects their productivity and milk quality. Adaptation to heat requires high energy costs, affects "neurohumoral" regulation and is accompanied by oxidative stress with the increased production of free radicals. The aim was to study the main biochemical parameters of goat milk and its antioxidant activity depending on the season of the year. Sampling was carried out in April, June, August and October. Analysis of the biochemical components and antioxidant activity of goat milk was performed using modern analytical systems. From spring to autumn, the mass fraction of true or crude proteins in goat milk increased by 14.6-63.7% or by 12.3-52.1%, and the mass fraction of caseins also increased by 13.6-60.6%. For vitamin C level and the total amount of water-soluble antioxidants, a pronounced gradual decrease from spring to autumn was observed. In the summer period, a small increase in the carotene level in milk (by 3.0-6.1% compared to April) was established. Vitamin A content increased by 86.5% (June) or by 70.3% (October) compared to April. Thus, the numerous significant changes in the major parameters of goat's milk depending on the season were revealed.

Prediction of detailed blood metabolic profile using milk infrared spectra and machine learning methods in dairy cattle

The adoption of preventive management decisions is crucial to dealing with metabolic impairments in dairy cattle. Various serum metabolites are known to be useful indicators of the health status of cows. In this study, we used milk Fourier-transform mid-infrared (FTIR) spectra and various machine learning (ML) algorithms to develop prediction equations for a panel of 29 blood metabolites, including those related to energy metabolism, liver function/hepatic damage, oxidative stress, inflammation/innate immunity, and minerals. For most traits, the data set comprised observations from 1,204 Holstein-Friesian dairy cows belonging to 5 herds. An exception was represented by β-hydroxybutyrate prediction, which contained observations from 2,701 multibreed cows pertaining to 33 herds. The best predictive model was developed using an automatic ML algorithm that tested various methods, including elastic net, distributed random forest, gradient boosting machine, artificial neural network, and stacking ensemble. These ML predictions were compared with partial least squares regression, the most commonly used method for FTIR prediction of blood traits. Performance of each model was evaluated using 2 cross-validation (CV) scenarios: 5-fold random (CV_r) and herd-out (CV_h). We also tested the best model's ability to classify values precisely in the 2 extreme tails, namely, the 25th (Q25) and 75th (Q75) percentiles (true-positive prediction scenario). Compared with partial least squares regression, ML algorithms achieved more accurate performance. Specifically, elastic net increased the R² value from 5% to 75% for CV_r and 2% to 139% for CV_h, whereas the stacking ensemble increased the R² value from 4% to 70% for CV_r and 4% to 150% for CV_h. Considering the best model, with the CV_r scenario, good prediction accuracies were obtained for glucose (R² = 0.81), urea (R² = 0.73), albumin (R² = 0.75), total reactive oxygen metabolites (R² = 0.79), total thiol groups (R² = 0.76), ceruloplasmin (R² = 0.74), total proteins (R² = 0.81), globulins (R² = 0.87), and Na (R² = 0.72). Good prediction accuracy in classifying extreme values was achieved for glucose (Q25 = 70.8%, Q75 = 69.9%), albumin (Q25 = 72.3%), total reactive oxygen metabolites (Q25 = 75.1%, Q75 = 74%), thiol groups (Q75 = 70.4%), total proteins (Q25 = 72.4%, Q75 = 77.2.%), globulins (Q25 = 74.8%, Q75 = 81.5%), and haptoglobin (Q75 = 74.4%). In conclusion, our study shows that FTIR spectra can be used to predict blood metabolites with relatively good accuracy, depending on trait, and are a promising tool for large-scale monitoring.

Novel prediction models for hyperketonemia using bovine milk Fourier-transform infrared spectroscopy

Metabolic diseases driven by negative energy balance in dairy cattle contribute to reduced milk production, increased disease incidence, culling, and death. Cow side tests for negative energy balance markers are available but are labor-intensive. Milk sample analysis using Fourier transform infrared spectroscopy (FTIR) allows for sampling numerous cows simultaneously. FTIR prediction models have moderate accuracy for hyperketonemia diagnosis (beta-hydroxybutyrate (BHB) ≥ 1.2 mmol/L). Most research using FTIR has focused on homogenous datasets and conventional prediction models, including partial least squares, linear discriminant analysis, and ElasticNet. Our objective was to evaluate more diverse modeling options, such as deep learning, gradient boosting machine models, and model ensembles for hyperketonemia classification. We compiled a sizable, heterogeneous dataset including milk FTIR and concurrent blood samples. Blood samples were tested for blood BHB, and wavenumber data was obtained from milk FTIR analysis. Using this dataset, we trained conventional prediction models and other options listed above. We demonstrate prediction model performance is similar for convolutional neural networks and ensemble models to simpler algorithm options. Results obtained from this study indicate that deep learning and model ensembles are potential algorithm options for predicting hyperketonemia in dairy cattle. Additionally, our results indicate hyperketonemia prediction models can be developed using heterogeneous datasets.

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.

Establishment and Validation of Fourier Transform Infrared Spectroscopy (FT–MIR) Methodology for the Detection of Linoleic Acid in Buffalo Milk

Buffalo milk is a dairy product that is considered to have a higher nutritional value compared to cow’s milk. Linoleic acid (LA) is an essential fatty acid that is important for human health. This study aimed to investigate and validate the use of Fourier transform mid-infrared spectroscopy (FT-MIR) for the quantification of the linoleic acid in buffalo milk. Three machine learning models were used to predict linoleic acid content, and random forest was employed to select the most important subset of spectra for improved model performance. The validity of the FT-MIR methods was evaluated in accordance with ICH Q2 (R1) guidelines using the accuracy profile method, and the precision, the accuracy, and the limit of quantification were determined. The results showed that Fourier transform infrared spectroscopy is a suitable technique for the analysis of linoleic acid, with a lower limit of quantification of 0.15 mg/mL milk. Our results showed that FT-MIR spectroscopy is a viable method for LA concentration analysis.

Predicting milk protein fractions using infrared spectroscopy and a gradient boosting machine for breeding purposes in Holstein cattle

In recent years, increasing attention has been focused on the genetic evaluation of protein fractions in cow milk with the aim of improving milk quality and technological characteristics. In this context, advances in high-throughput phenotyping by Fourier transform infrared (FTIR) spectroscopy offer the opportunity for large-scale, efficient measurement of novel traits that can be exploited in breeding programs as indicator traits. We took milk samples from 2,558 Holstein cows belonging to 38 herds in northern Italy, operating under different production systems. Fourier transform infrared spectra were collected on the same day as milk sampling and stored for subsequent analysis. Two sets of data (i.e., phenotypes and FTIR spectra) collected in 2 different years (2013 and 2019-2020) were compiled. The following traits were assessed using HPLC: true protein, major casein fractions [α_S1-casein (CN), α_S2-CN, β-CN, κ-CN, and glycosylated-κ-CN], and major whey proteins (β-lactoglobulin and α-lactalbumin), all of which were measured both in grams per liter (g/L) and proportion of total nitrogen (% N). The FTIR predictions were calculated using the gradient boosting machine technique and tested by 3 different cross-validation (CRV) methods. We used the following CRV scenarios: (1) random 10-fold, which randomly split the whole into 10-folds of equal size (9-folds for training and 1-fold for validation); (2) herd/date-out CRV, which assigned 80% of herd/date as the training set with independence of 20% of herd/date assigned as the validation set; (3) forward/backward CRV, which split the data set in training and validation set according with the year of milk sampling (FTIR and gold standard data assessed in 2013 or 2019-2020) using the "old" and "new" databases for training and validation, and vice-versa with independence among them; (4) the CRV for genetic parameters (CRV-gen), where animals without pedigree as assigned as a fixed training population and animals with pedigree information was split in 5-folds, in which 1-fold was assigned to the fixed training population, and 4-folds were assigned to the validation set (independent from the training set). The results (i.e., measures and predictions) of CRV-gen were used to infer the genetic parameters for gold standard laboratory measurements (i.e., proteins assessed with HPLC) and FTIR-based predictions considering the CRV-gen scenario from a bi-trait animal model using single-step genomic BLUP. We found that the prediction accuracies of the gradient boosting machine equations differed according to the way in which the proteins were expressed, achieving higher accuracy when expressed in g/L than when expressed as % N in all CRV scenarios. Concerning the reproducibility of the equations over the different years, the results showed no relevant differences in predictive ability between using "old" data as the training set and "new" data as the validation set and vice-versa. Comparing the additive genetic variance estimates for milk protein fractions between the FTIR predicted and HPLC measures, we found reductions of -19.7% for milk protein fractions expressed in g/L, and -21.19% expressed as % N. Although we found reductions in the heritability estimates, they were small, with values ranging from -1.9 to -7.25% for g/L, and -1.6 to -7.9% for % N. The posterior distributions of the additive genetic correlations (r_a) between the FTIR predictions and the laboratory measurements were generally high (>0.8), even when the milk protein fractions were expressed as % N. Our results show the potential of using FTIR predictions in breeding programs as indicator traits for the selection of animals to enhance milk protein fraction contents. We expect acceptable responses to selection due to the high genetic correlations between HPLC measurements and FTIR predictions.

Spectroscopic technologies and data fusion: Applications for the dairy industry

Increasing consumer awareness, scale of manufacture, and demand to ensure safety, quality and sustainability have accelerated the need for rapid, reliable, and accurate analytical techniques for food products. Spectroscopy, coupled with Artificial Intelligence-enabled sensors and chemometric techniques, has led to the fusion of data sources for dairy analytical applications. This article provides an overview of the current spectroscopic technologies used in the dairy industry, with an introduction to data fusion and the associated methodologies used in spectroscopy-based data fusion. The relevance of data fusion in the dairy industry is considered, focusing on its potential to improve predictions for processing traits by chemometric techniques, such as principal component analysis (PCA), partial least squares regression (PLS), and other machine learning algorithms.

Predicting nitrogen use efficiency, nitrogen loss and dry matter intake of individual dairy cows in late lactation by including mid-infrared spectra of milk samples

BACKGROUND

Nitrate leaching to groundwater and surface water and ammonia volatilization from dairy farms have negative impacts on the environment. Meanwhile, the increasing demand for dairy products will result in more pollution if N losses are not controlled. Therefore, a more efficient, and environmentally friendly production system is needed, in which nitrogen use efficiency (NUE) of dairy cows plays a key role. To genetically improve NUE, extensively recorded and cost-effective proxies are essential, which can be obtained by including mid-infrared (MIR) spectra of milk in prediction models for NUE. This study aimed to develop and validate the best prediction model of NUE, nitrogen loss (NL) and dry matter intake (DMI) for individual dairy cows in China.

RESULTS

A total of 86 lactating Chinese Holstein cows were used in this study. After data editing, 704 records were obtained for calibration and validation. Six prediction models with three different machine learning algorithms and three kinds of pre-processed MIR spectra were developed for each trait. Results showed that the coefficient of determination (R²) of the best model in within-herd validation was 0.66 for NUE, 0.58 for NL and 0.63 for DMI. For external validation, reasonable prediction results were only observed for NUE, with R² ranging from 0.58 to 0.63, while the R² of the other two traits was below 0.50. The infrared waves from 973.54 to 988.46 cm^-1 and daily milk yield were the most important variables for prediction.

CONCLUSION

The results showed that individual NUE can be predicted with a moderate accuracy in both within-herd and external validations. The model of NUE could be used for the datasets that are similar to the calibration dataset. The prediction models for NL and 3-day moving average of DMI (DMI_a) generated lower accuracies in within-herd validation. Results also indicated that information of MIR spectra variables increased the predictive ability of models. Additionally, pre-processed MIR spectra do not result in higher accuracy than original MIR spectra in the external validation. These models will be applied to large-scale data to further investigate the genetic architecture of N efficiency and further reduce the adverse impacts on the environment after more data is collected.

Fourier transform mid-infrared milk screening to improve milk production and processing

Milk mid-infrared spectrometry has been used for many years to quantify major milk compounds. Recently, much research has been conducted to extend the use of this technology to predict new, relevant phenotypes to assess the animals' welfare and the nutritional quality of milk, as well as its technological quality and environmental footprint. The transition from the research stage to field implementation is not easy, due to intrinsic and extrinsic constraints, but some developments can be considered to address these issues.

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.

In-line near-infrared analysis of milk coupled with machine learning methods for the daily prediction of blood metabolic profile in dairy cattle

Precision livestock farming technologies are used to monitor animal health and welfare parameters continuously and in real time in order to optimize nutrition and productivity and to detect health issues at an early stage. The possibility of predicting blood metabolites from milk samples obtained during routine milking by means of infrared spectroscopy has become increasingly attractive. We developed, for the first time, prediction equations for a set of blood metabolites using diverse machine learning methods and milk near-infrared spectra collected by the AfiLab instrument. Our dataset was obtained from 385 Holstein Friesian dairy cows. Stacking ensemble and multi-layer feedforward artificial neural network outperformed the other machine learning methods tested, with a reduction in the root mean square error of between 3 and 6% in most blood parameters. We obtained moderate correlations (r) between the observed and predicted phenotypes for γ-glutamyl transferase (r = 0.58), alkaline phosphatase (0.54), haptoglobin (0.66), globulins (0.61), total reactive oxygen metabolites (0.60) and thiol groups (0.57). The AfiLab instrument has strong potential but may not yet be ready to predict the metabolic stress of dairy cows in practice. Further research is needed to find out methods that allow an improvement in accuracy of prediction equations.

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.

Opportunities to Harness High-Throughput and Novel Sensing Phenotypes to Improve Feed Efficiency in Dairy Cattle

Simple Summary

Sensors, routinely collected on-farm tests, and other repeatable, high-throughput measurements can provide novel phenotype information on a frequent basis. Information from these sensors and high-throughput measurements could be harnessed to monitor or predict individual dairy cow feed intake. Predictive algorithms would allow for genetic selection of animals that consume less feed while producing the same amount of milk. Improved monitoring of feed intake could reduce the cost of milk production, improve animal health, and reduce the environmental impact of the dairy industry. Moreover, data from these information sources could aid in animal management (e.g., precision feeding and health detection). In order to implement tools, the relationship of measurements with feed intake needs to be established and prediction equations developed. Lastly, consideration should be given to the frequency of data collection, the need for standardization of data and other potential limitations of tools in the prediction of feed intake. This review summarizes measurements of feed efficiency, factors that may impact the efficiency and feed consumption of an animal, tools that have been researched and new traits that could be utilized for the prediction of feed intake and efficiency, and prediction equations for feed intake and efficiency presented in the literature to date.

Abstract

Feed for dairy cattle has a major impact on profitability and the environmental impact of farms. Sustainable dairy production relies on continued improvement in feed efficiency as a way to reduce costs and nutrient loss from feed. Advances in breeding, feeding and management have led to the dilution of maintenance energy and thus more efficient dairy cattle. Still, many additional opportunities are available to improve individual animal feed efficiency. Sensing technologies such as wearable sensors, image-based and high-throughput phenotyping technologies (e.g., milk testing) are becoming more available on commercial farm. The application of these technologies as indicator traits for feed intake and efficiency related traits would be advantageous to provide additional information to predict and manage feed efficiency. This review focuses on precision livestock technologies and high-throughput phenotyping in use today as well as those that could be developed in the future as possible indicators of feed intake. Several technologies such as milk spectral data, activity, rumen measures, and image-based phenotypes have been associated with feed intake. Future applications will depend on the ability to repeatably measure and calibrate these data across locations, so that they can be integrated for use in predicting and managing feed intake and efficiency on farm.

Milk metabolites and fatty acids as noninvasive biomarkers of metabolic status and energy balance in early-lactation cows

The objective was to study the effects of week of lactation (WOL) and experimental nutrient restriction on concentrations of selected milk metabolites and fatty acids (FA), and assess their potential as biomarkers of energy status in early-lactation cows. To study WOL effects, 17 multiparous Holstein cows were phenotyped from calving until 7 WOL while allowed ad libitum intake of a lactation diet. Further, to study the effects of nutrient restriction, 8 of these cows received a diet containing 48% straw (high-straw) for 4 d starting at 24 ± 3 days in milk (mean ± SD), and 8 cows maintained on the lactation diet were sampled to serve as controls. Blood and milk samples were collected weekly for the WOL data set, and daily from d -1 to 3 of nutrient restriction (or control) for the nutritional challenge data set. Milk β-hydroxybutyrate (BHB), isocitrate, glucose, glucose-6-phosphate (glucose-6P), galactose, glutamate, creatinine, uric acid, and N-acetyl-β-d-glucosaminidase activity (NAGase) were analyzed in p.m. and a.m. samples, and milk FA were analyzed in pooled p.m. and a.m. samples. Average energy balance (EB) per day ranged from -27 MJ/d to neutral when cows received the lactation total mixed ration, and from -109 to -87 ± 7 MJ/d for high-straw (least squares means ± standard error of the mean). Plasma nonesterified FA concentration was 1.67 ± 0.13 mM and BHB was 2.96 ± 0.39 mM on the d 3 of high-straw (least squares means ± standard error of the mean). Milk concentrations of BHB, glucose, glucose-6P, glutamate, and uric acid differed significantly between p.m. and a.m. milkings. Milk isocitrate, glucose-6P, creatinine, and NAGase decreased, whereas milk glucose and galactose increased with WOL. Changes in milk BHB, isocitrate, glucose, glucose-6P, and creatinine were concordant during early lactation and in response to nutrient restriction. Milk galactose and NAGase were modulated by WOL only, whereas glutamate and uric acid concentrations responded to nutrient restriction only. The high-straw increased milk concentrations of FA potentially mobilized from adipose tissue (e.g., C18:0 and cis-9 C18:1 and sum of odd- and branched-chain FA (OBCFA) with carbon chain greater than 16; ∑ OBCFA >C16), and decreased concentrations of FA synthesized de novo by the mammary gland (e.g., sum of FA with 6 to 15 carbons; ∑ C6:0 to C15:0). Similar observations were made during early lactation. Plasma nonesterified FA concentrations had the best single linear regression with EB (R² = 0.62). Milk isocitrate, Σ C6:0 to C15:0. and cis-9 C18:1 had the best single linear regressions with EB (R² ≥ 0.44). Milk BHB, isocitrate, galactose, glutamate, and creatinine explained up to 64% of the EB variation observed in the current study using multiple linear regression. Milk concentrations of ∑ C6:0 to C15:0, C18:0, cis-9 C18:1, and ∑ OBCFA >C16 presented some of the best correlations and regressions with other indicators of metabolic status, lipomobilization, and EB, and their responses were concordant during early lactation and during experimental nutrient restriction. Metabolites and FA secreted in milk may serve as noninvasive indicators of metabolic status and EB of early-lactation cows.

Multiple Breeds and Countries' Predictions of Mineral Contents in Milk from Milk Mid-Infrared Spectrometry

Measuring the mineral composition of milk is of major interest in the dairy sector. This study aims to develop and validate robust multi-breed and multi-country models predicting the major minerals through milk mid-infrared spectrometry using partial least square regressions. A total of 1281 samples coming from five countries were analyzed to obtain spectra and in ICP-AES to measure the mineral reference contents. Models were built from records coming from four countries (n = 1181) and validated using records from the fifth country, Austria (n = 100). The importance of including local samples was tested by integrating 30 Austrian samples in the model while validating with the remaining 70 samples. The best performances were achieved using this second set of models, confirming the need to cover the spectral variability of a country before making a prediction. Validation root mean square errors were 54.56, 63.60, 7.30, 59.87, and 152.89 mg/kg for Na, Ca, Mg, P, and K, respectively. The built models were applied on the Walloon milk recording large-scale spectral database, including 3,510,077. The large-scale predictions on this dairy herd improvement database provide new insight regarding the minerals' variability in the population, as well as the effect of parity, stage of lactation, breeds, and seasons.

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.

Full model selection using regression trees for numeric predictions of biomarkers for metabolic challenges in dairy cows

Dairy cows suffer poor metabolic adaptation syndrome (PMAS)¹ during early post-calving periods caused by negative energy balance. Measurement of blood beta-hydroxy butyric acid (BHBA)² and blood non-esterified fatty acids (NEFA)³ allow early and accurate detection of negative energy balance. Machine learning prediction of blood BHBA and blood NEFA using milk testing samples represents an opportunity to identify at-risk animals, using less labor than direct blood testing methods. Routine milk testing on modern dairies and computer record keeping provide an immense amount of data which can then be used in machine learning models. Previous research for predicting blood metabolites using Fourier-transform infrared spectroscopy (FTIR)⁴ milk data has focused mainly on individual models rather than a comparison among the models. Full model selection is the process of comparing different combinations of pre-processing methods, variable selection, and statistical learning algorithms to determine which model results in the lowest prediction error for a given dataset. For this project we used a full model selection approach with regression trees (rtFMS)⁵ . rtFMS uses the cross-validated performance of different model configurations to feed a regression tree for selecting a final model. A total of 384 possible model configurations (algorithms, predictors and data preprocessing options) for each outcome (blood BHBA and blood NEFA) were considered in the rtFMS technique. rtFMS allows direct comparison of multiple modeling approaches reducing bias due to empirical knowledge, modeling habits, or preferences, identifying the model with minimal root mean squared prediction error (RMSE)⁶ . An elastic net regression model was selected as the best performing model for both biomarkers. The input data for blood BHBA predictions were FTIR milk spectra, with a second derivative pre-processing, and a filter with 212 wave numbers, obtaining RMSE = 0.354 (0.328-0.392). The best performing model for blood NEFA had input data of FTIR milk spectra, with a second derivative pre-processing, and a filter with 212 wave numbers filter along with the time of milking, obtaining RMSE = 0.601 (0.564-0.654). The comparison of multiple modeling strategies, conducted by rtFMS, present an option for improved FTIR prediction models of blood BHBA and blood NEFA by reducing error due to human bias. The implementation of rtFMS to design future prediction models can guide model inputs and features. Our prediction models have the potential to increase early detection of metabolic disorders in dairy cows during the transition period.

Invited Review: Ketosis Diagnosis and Monitoring in High-Producing Dairy Cows

This work reviews the current impact and manifestation of ketosis (hyperketonemia) in dairy cattle, emphasizing the practical use of laboratory methods, field tests, and milk data to monitoring this disease. Ketosis is a major issue in high-producing cows, easily reaching a prevalence of 20% during early postpartum when the negative energy balance is well established. Its economic losses, mainly related to decreasing milk yield, fertility, and treatment costs, have been estimated up to €250 per case of ketosis/year, which can double if associated diseases are considered. A deep relationship between subclinical or clinical ketosis and negative energy balance and related production diseases can be observed mainly in the first two months postpartum. Fourier transform infrared spectrometry methods gradually take place in laboratory routine to evaluates body ketones (e.g., beta-hydroxybutyrate) and probably will accurately substitute cowside blood and milk tests at a farm in avenir. Fat to protein ratio and urea in milk are largely evaluated each month in dairy farms indicating animals at risk of hyperketonemia. At preventive levels, other than periodical evaluation of body condition score and controlling modifiable or identifying non-modifiable risk factors, the ruminatory activity assessment during the peripartum seems to be a valuable tool at farms. We conclude that a technological advance progressively takes place to mitigate the effects of these metabolic diseases, which challenge the high-yielding cows.

Influence of Estrus on the Milk Characteristics and Mid-Infrared Spectra of Dairy Cows

Simple Summary

Some studies have confirmed the variation in milk profiles when dairy cows show estrus. However, only a few milk components, such as fat, protein, and lactose, have been investigated so far, and thus any changes in the many other parts of milk’s composition due to estrus are unknown. Milk mid-infrared (MIR) spectra consist of wavenumbers, which provide insight into the chemical composition of milk. The MIR spectrum reflects the global composition of milk, but this information is currently underused. In this study, we considered MIR wavenumbers as traits, and directly studied the spectral information as a way to study the estrus of dairy cows linked to milk composition. This research provides a deeper understanding of the milk MIR spectrum and may lead to new approaches for estrus detection in dairy cows from routine milk analysis, thereby guiding an opportune insemination time.

Abstract

Milk produced by dairy cows is a complex combination of many components. However, at present, changes in only a few milk components (e.g., fat, protein, and lactose) during the estrus cycle in dairy cows have been documented. Mid-infrared (MIR) spectroscopy is a worldwide method routinely used for milk analysis, as MIR spectra reflect the global composition of milk. Therefore, this study aimed to investigate the changes in milk MIR spectra and milk production traits (fat, protein, lactose, urea, total solids (TS), and solid not fat (SnF)) due to estrus. Cows that were successfully inseminated, leading to conception, were included. Cows confirmed to be pregnant were considered to be in estrus at the day of insemination (day 0). A general linear mixed model, which included the random effect of cows, the fixed classification effects of parity number, days in relation to estrus, as well as the interaction between parity number and days in relation to estrus, was applied to investigate the changes in milk production traits and 1060 milk infrared wavenumbers, ranging from 925 to 5011 cm⁻¹, of 371 records from 162 Holstein cows on the days before (day −3, day −2, and day −1) and on the day of estrus (day 0). The days in relation to estrus had a significant effect on fat, protein, urea, TS, and SnF, whose contents increased from day −3 to day 0. Lactose did not seem to be significantly influenced by the occurrence of estrus. The days in relation to estrus had significant effects on the majority of the wavenumbers. Besides, we found that some of the wavenumbers in the water absorption regions were significantly changed on the days before and on the day of estrus. This suggests that these wavenumbers may contain useful information. In conclusion, the changes in the milk composition due to estrus can be observed through the analysis of the milk MIR spectrum. Further analyses are warranted to more deeply explore the potential use of milk MIR spectra in the detection of estrus.

Evaluating the performance of machine learning methods and variable selection methods for predicting difficult-to-measure traits in Holstein dairy cattle using milk infrared spectral data

Fourier-transform infrared (FTIR) spectroscopy is a powerful high-throughput phenotyping tool for predicting traits that are expensive and difficult to measure in dairy cattle. Calibration equations are often developed using standard methods, such as partial least squares (PLS) regression. Methods that employ penalization, rank-reduction, and variable selection, as well as being able to model the nonlinear relations between phenotype and FTIR, might offer improvements in predictive ability and model robustness. This study aimed to compare the predictive ability of 2 machine learning methods, namely random forest (RF) and gradient boosting machine (GBM), and penalized regression against PLS regression for predicting 3 phenotypes differing in terms of biological meaning and relationships with milk composition (i.e., phenotypes measurable directly and not directly in milk, reflecting different biological processes which can be captured using milk spectra) in Holstein-Friesian cattle under 2 cross-validation scenarios. The data set comprised phenotypic information from 471 Holstein-Friesian cows, and 3 target phenotypes were evaluated: (1) body condition score (BCS), (2) blood β-hydroxybutyrate (BHB, mmol/L), and (3) κ-casein expressed as a percentage of nitrogen (κ-CN, % N). The data set was split considering 2 cross-validation scenarios: samples-out random in which the population was randomly split into 10-folds (8-folds for training and 1-fold for validation and testing); and herd/date-out in which the population was randomly assigned to training (70% herd), validation (10%), and testing (20% herd) based on the herd and date in which the samples were collected. The random grid search was performed using the training subset for the hyperparameter optimization and the validation set was used for the generalization of prediction error. The trained model was then used to assess the final prediction in the testing subset. The grid search for penalized regression evidenced that the elastic net (EN) was the best regularization with increase in predictive ability of 5%. The performance of PLS (standard model) was compared against 2 machine learning techniques and penalized regression using 2 cross-validation scenarios. Machine learning methods showed a greater predictive ability for BCS (0.63 for GBM and 0.61 for RF), BHB (0.80 for GBM and 0.79 for RF), and κ-CN (0.81 for GBM and 0.80 for RF) in samples-out cross-validation. Considering a herd/date-out cross-validation these values were 0.58 (GBM and RF) for BCS, 0.73 (GBM and RF) for BHB, and 0.77 (GBM and RF) for κ-CN. The GBM model tended to outperform other methods in predictive ability around 4%, 1%, and 7% for EN, RF, and PLS, respectively. The prediction accuracies of the GBM and RF models were similar, and differed statistically from the PLS model in samples-out random cross-validation. Although, machine learning techniques outperformed PLS in herd/date-out cross-validation, no significant differences were observed in terms of predictive ability due to the large standard deviation observed for predictions. Overall, GBM achieved the highest accuracy of FTIR-based prediction of the different phenotypic traits across the cross-validation scenarios. These results indicate that GBM is a promising method for obtaining more accurate FTIR-based predictions for different phenotypes in dairy cattle.

Dairy Cows' Health during Alpine Summer Grazing as Assessed by Milk Traits, Including Differential Somatic Cell Count: A Case Study from Italy

Simple Summary

Dairy herds in alpine areas usually adopt summer grazing, mainly to reduce feeding costs. This practice is related to the maintenance of local traditions and to the manufacturing of niche dairy products. However, it is important to assess the impact of this practice on cattle health. This case study investigated how milk-related health traits vary across extensive grazing during the summer period, using data collected in a dairy herd whose cows were repeatedly controlled for individual milk samples. Although the transition from barn farming to pasture led to a reduction in milk production, proper grazing management can make dairy cows more resilient in terms of udder health and metabolic efficiency. Findings of the present research report suggested that pasture can be adopted to maintain dairy herd sustainability without impairing animal health.

Abstract

Extensive summer grazing is a dairy herd management practice frequently adopted in mountainous areas. Nowadays, this activity is threatened by its high labour demand, but it is fundamental for environmental, touristic and economic implications, as well as for the preservation of social and cultural traditions. Scarce information on the effects of such low-input farming systems on cattle health is available. Therefore, the present case study aimed at investigating how grazing may affect the health status of dairy cows by using milk traits routinely available from the national milk recording scheme. The research involved a dairy herd of 52 Simmental and 19 Holstein × Simmental crossbred cows. The herd had access to the pasture according to a rotational grazing scheme from late spring up to the end of summer. A total of 616 test day records collected immediately before and during the grazing season were used. Individual milk yield was registered during the milking procedure. Milk samples were analysed for composition (fat, protein, casein and lactose contents) and health-related milk indicators (electrical conductivity, urea and β-hydroxybutyrate) using mid-infrared spectroscopy. Somatic cell count (SCC) and differential SCC were also determined. Data were analysed with a linear mixed model, which included the fixed effects of the period of sampling, cow breed, stage of lactation and parity, and the random effects of cow nested within breed and the residual. The transition from barn farming to pasture had a negative effect on milk yield, together with a small deterioration of fat and protein percentages. Health-related milk indicators showed a minor deterioration of the fat to protein ratio, differential SCC and electrical conductivity, particularly towards the end of the grazing season, whereas the somatic cell score and β-hydroxybutyrate were relatively constant. Overall, the study showed that, when properly managed, pasture grazing does not have detrimental effects on dairy cows in terms of udder health and efficiency. Therefore, the proper management of cows on pasture can be a valuable solution to preserve the economic, social and environmental sustainability of small dairy farms in the alpine regions, without impairing cows’ health.

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.

Fourier Transform Infrared Spectroscopy as a Tool to Study Milk Composition Changes in Dairy Cows Attributed to Housing Modifications to Improve Animal Welfare

Animal welfare status is assessed today through visual evaluations requiring an on-farm visit. A convenient alternative would be to detect cow welfare status directly in milk samples, already routinely collected for milk recording. The objective of this study was to propose a novel approach to demonstrate that Fourier transform infrared (FTIR) spectroscopy can detect changes in milk composition related to cows subjected to movement restriction at the tie stall with four tie-rail configurations varying in height and position (TR1, TR2, TR3 and TR4). Milk mid-infrared spectra were collected on weekly basis. Long-term average spectra were calculated for each cow using spectra collected in weeks 8-10 of treatment. Principal component analysis was applied to spectral averages and the scores of principal components (PCs) were tested for treatment effect by mixed modelling. PC7 revealed a significant treatment effect (p = 0.01), particularly for TR3 (configuration with restricted movement) vs. TR1 (recommended configuration) (p = 0.03). The loading spectrum of PC7 revealed high loadings at wavenumbers that could be assigned to biomarkers related to negative energy balance, such as β-hydroxybutyrate, citrate and acetone. This observation suggests that TR3 might have been restrictive for cows to access feed. Milk FTIR spectroscopy showed promising results in detecting welfare status and housing conditions in dairy cows.

Development of a subacute ruminal acidosis risk score and its prediction using milk mid-infrared spectra in early-lactation cows

A routine monitoring for subacute ruminal acidosis (SARA) on the individual level could support the minimization of economic losses and the ensuring of animal welfare in dairy cows. The objectives of this study were (1) to develop a SARA risk score (SRS) by combining information from different data acquisition systems to generate an integrative indicator trait, (2) the investigation of associations of the SRS with feed analysis data, blood characteristics, performance data, and milk composition, including the fatty acid (FA) profile, (3) the development of a milk mid-infrared (MIR) spectra-based prediction equation for this novel reference trait SRS, and (4) its application to an external data set consisting of MIR data of test day records to investigate the association between the MIR-based predictions of the SRS and the milk FA profile. The primary data set, which was used for the objectives (1) to (3), consisted of data collected from 10 commercial farms with a total of 100 Holstein cows in early lactation. The data comprised barn climate parameters, pH and temperature logging from intrareticular measurement boluses, as well as jaw movement and locomotion behavior recordings of noseband-sensor halters and pedometers. Further sampling and data collection included feed samples, blood samples, milk performance, and milk samples, whereof the latter were used to get the milk MIR spectra and to estimate the main milk components, the milk FA profile, and the lactoferrin content. Because all measurements were characterized by different temporal resolutions, the data preparation consisted of an aggregation into values on a daily basis and merging it into one data set. For the development of the SRS, a total of 7 traits were selected, which were derived from measurements of pH and temperature in the reticulum, chewing behavior, and milk yield. After adjustment for fixed effects and standardization, these 7 traits were combined into the SRS using a linear combination and directional weights based on current knowledge derived from literature studies. The secondary data set was used for objective (4) and consisted of test day records of the entire herds, including performance data, milk MIR spectra and MIR-predicted FA. At farm level, it could be shown that diets with higher proportions of concentrated feed resulted in both lower daily mean pH and higher SRS values. On the individual level, an increased SRS could be associated with a modified FA profile (e.g., lower levels of short- and medium-chain FA, higher levels of C17:0, odd- and branched-chain FA). Furthermore, a milk MIR-based partial least squares regression model with a moderate predictability was established for the SRS. This work provides the basis for the development of routine SARA monitoring and demonstrates the high potential of milk composition-based assessment of the health status of lactating cows.

Validation of milk mid-infrared spectroscopy for predicting the metabolic status of lactating dairy cows in Australia

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

A comparison of 4 different machine learning algorithms to predict lactoferrin content in bovine milk from mid-infrared spectra

Lactoferrin (LF) is a glycoprotein naturally present in milk. Its content varies throughout lactation, but also with mastitis; therefore it is a potential additional indicator of udder health beyond somatic cell count. Condequently, there is an interest in quantifying this biomolecule routinely. First prediction equations proposed in the literature to predict the content in milk using milk mid-infrared spectrometry were built using partial least square regression (PLSR) due to the limited size of the data set. Thanks to a large data set, the current study aimed to test 4 different machine learning algorithms using a large data set comprising 6,619 records collected across different herds, breeds, and countries. The first algorithm was a PLSR, as used in past investigations. The second and third algorithms used partial least square (PLS) factors combined with a linear and polynomial support vector regression (PLS + SVR). The fourth algorithm also used PLS factors, but included in an artificial neural network with 1 hidden layer (PLS + ANN). The training and validation sets comprised 5,541 and 836 records, respectively. Even if the calibration prediction performances were the best for PLS + polynomial SVR, their validation prediction performances were the worst. The 3 other algorithms had similar validation performances. Indeed, the validation root mean squared error (RMSE) ranged between 162.17 and 166.75 mg/L of milk. However, the lower standard deviation of cross-validation RMSE and the better normality of the residual distribution observed for PLS + ANN suggest that this modeling was more suitable to predict the LF content in milk from milk mid-infrared spectra (R²v = 0.60 and validation RMSE = 162.17 mg/L of milk). This PLS +ANN model was then applied to almost 6 million spectral records. The predicted LF showed the expected relationships with milk yield, somatic cell score, somatic cell count, and stage of lactation. The model tended to underestimate high LF values (higher than 600 mg/L of milk). However, if the prediction threshold was set to 500 mg/L, 82% of samples from the validation having a content of LF higher than 600 mg/L were detected. Future research should aim to increase the number of those extremely high LF records in the calibration set.

Fatty Acid Profiles from Routine Milk Recording as a Decision Tool for Body Weight Change of Dairy Cows after Calving

Cows mobilize body reserves during early lactation, which is reflected in the milk fatty acid (FA) profile. Milk FA can be routinely predicted by Fourier-transform infrared (FTIR) spectroscopy, and be, thus, used to develop an early indicator for bodyweight change (BWC) in early lactating cows in commercial dairy farms. Cow records from 165 herds in Denmark between 2015 and 2017 were used with bodyweight (BW) records at each milking from floor scales in automatic milking systems. Milk FA in monthly test-day samples was predicted by FTIR. Predictions of BWC were based on a random forest model and included parity, stage of lactation, and test day milk production and components (fat, protein, and FA). Bodyweight loss was mainly explained by decreased short-chain FA (C4:0-C10:0) and increased C18:0 FA. The root mean square error (RMSE) of prediction after cross-validation was 1.79 g/kg of BW (R2 of 0.94). Model evaluation with previously unseen BWC records resulted in reduced prediction performance (RMSE of 2.33 g/kg of BW; R2 of 0.31). An early warning system may be implemented for cows with a large BW loss during early lactation based on milk FA profiles, but model performance should be improved, ideally by using the full FTIR milk spectra.

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.

Predicting bovine tuberculosis status of dairy cows from mid-infrared spectral data of milk using deep learning

Bovine tuberculosis (bTB) is a zoonotic disease in cattle that is transmissible to humans, distributed worldwide, and considered endemic throughout much of England and Wales. Mid-infrared (MIR) analysis of milk is used routinely to predict fat and protein concentration, and is also a robust predictor of several other economically important traits including individual fatty acids and body energy. This study predicted bTB status of UK dairy cows using their MIR spectral profiles collected as part of routine milk recording. Bovine tuberculosis data were collected as part of the national bTB testing program for Scotland, England, and Wales; these data provided information from over 40,500 bTB herd breakdowns. Corresponding individual cow life-history data were also available and provided information on births, movements, and deaths of all cows in the study. Data relating to single intradermal comparative cervical tuberculin (SICCT) skin-test results, culture, slaughter status, and presence of lesions were combined to create a binary bTB phenotype labeled 0 to represent nonresponders (i.e., healthy cows) and 1 to represent responders (i.e., bTB-affected cows). Contemporaneous individual milk MIR spectral data were collected as part of monthly routine milk recording and matched to bTB status of individual animals on the single intradermal comparative cervical tuberculin test date (±15 d). Deep learning, a sub-branch of machine learning, was used to train artificial neural networks and develop a prediction pipeline for subsequent use in national herds as part of routine milk recording. Spectra were first converted to 53 × 20-pixel PNG images, then used to train a deep convolutional neural network. Deep convolutional neural networks resulted in a bTB prediction accuracy (i.e., the number of correct predictions divided by the total number of predictions) of 71% after training for 278 epochs. This was accompanied by both a low validation loss (0.71) and moderate sensitivity and specificity (0.79 and 0.65, respectively). To balance data in each class, additional training data were synthesized using the synthetic minority over sampling technique. Accuracy was further increased to 95% (after 295 epochs), with corresponding validation loss minimized (0.26), when synthesized data were included during training of the network. Sensitivity and specificity also saw a 1.22- and 1.45-fold increase to 0.96 and 0.94, respectively, when synthesized data were included during training. We believe this study to be the first of its kind to predict bTB status from milk MIR spectral data. We also believe it to be the first study to use milk MIR spectral data to predict a disease phenotype, and posit that the automated prediction of bTB status at routine milk recording could provide farmers with a robust tool that enables them to make early management decisions on potential reactor cows, and thus help slow the spread of bTB.

Infrared Spectrometry as a High-Throughput Phenotyping Technology to Predict Complex Traits in Livestock Systems

High-throughput phenotyping technologies are growing in importance in livestock systems due to their ability to generate real-time, non-invasive, and accurate animal-level information. Collecting such individual-level information can generate novel traits and potentially improve animal selection and management decisions in livestock operations. One of the most relevant tools used in the dairy and beef industry to predict complex traits is infrared spectrometry, which is based on the analysis of the interaction between electromagnetic radiation and matter. The infrared electromagnetic radiation spans an enormous range of wavelengths and frequencies known as the electromagnetic spectrum. The spectrum is divided into different regions, with near- and mid-infrared regions being the main spectral regions used in livestock applications. The advantage of using infrared spectrometry includes speed, non-destructive measurement, and great potential for on-line analysis. This paper aims to review the use of mid- and near-infrared spectrometry techniques as tools to predict complex dairy and beef phenotypes, such as milk composition, feed efficiency, methane emission, fertility, energy balance, health status, and meat quality traits. Although several research studies have used these technologies to predict a wide range of phenotypes, most of them are based on Partial Least Squares (PLS) and did not considered other machine learning (ML) techniques to improve prediction quality. Therefore, we will discuss the role of analytical methods employed on spectral data to improve the predictive ability for complex traits in livestock operations. Furthermore, we will discuss different approaches to reduce data dimensionality and the impact of validation strategies on predictive quality.

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.

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: Potential prediction of vitamin B12 concentration based on mid-infrared spectral data using Holstein Dairy Herd Improvement milk samples

The purpose of this study was (1) to predict the quantitative concentration of vitamin B₁₂ in milk using mid-infrared (MIR) spectrometry, and (2) to evaluate the potential of MIR spectra to discriminate different clusters of records based on their B₁₂ concentration. Milk samples were collected from 4,340 Holstein cows between 3 and 592 d in milk and located in 100 herds. Samples were taken using in-line milk meters and divided into 2 aliquots: one for MIR spectrometry and the other for B₁₂ concentration reference analyses by radioassay. Analyses were performed on 311 selected spectral wavelengths. A partial least squares regression model was built to quantify B₁₂ concentration. Discriminant analysis was used to isolate B₁₂ concentration clusters. A B₁₂ concentration threshold was set at 442 ng/dL, because this represents the cutoff value for a 250-mL glass of milk to fulfill 46% of the daily vitamin B₁₂ recommended dietary allowance for individuals 14 yr or older. For each analysis, records coming from two-thirds of herds were used to calibrate prediction equations, and the remaining records (one-third of herds for validation) were used to assess the prediction performance. In the case of discriminant analysis, validation sets were divided into evaluation sets (one-third of herds) to obtain alternate probability cutoffs and in test sets (two-thirds of herds) to validate equations. Spectral and B₁₂ concentration outliers were identified by calculating standardized Mahalanobis distance and with a residual analysis, respectively (n = 3,154). Regarding quantitative B₁₂ concentration, cross-validation and validation coefficients of determination averaged 0.51 and 0.46, respectively, which are relatively low, which would limit the potential use of the developed quantitative equations. In addition, root mean square errors of prediction of cross validation and validation sets averaged 88.9 and 94.7 ng/dL, respectively. Area under the receiver operating characteristic curve of test sets averaged 0.81 based on the 442 ng/dL threshold, which could be considered to represent good accuracy of classification. However, the false discovery rate averaged 36%. In summary, models predicting quantitative B₁₂ concentration had low cross-validation and validation coefficients of determination, limiting their use, but the proposed discriminant models could be used to identify milk samples with naturally high B12.

Economic value of information from an alert system on physiological imbalance in fresh cows

Physiological imbalance is an abnormal physiological condition that cannot be directly observed but is assumed to precede subclinical and clinical diseases in the beginning of lactation. Alert systems to detect the physiological imbalance in a cow using Fourier transform mid-infrared spectroscopy in milk have been developed. The objective of this study was to estimate the value of information provided from such system with different indicator accuracies, herd prevalence and prices. A decision tree was created to model the probabilities of detection and associated costs of test outcome, intervention and occurrence of disease. We assumed that the negative effect of physiological imbalance was the development of subclinical ketosis and that this negative effect was prevented by drenching the cows with propylene glycol for 5 days. We simulated the economic impact of subclinical ketosis mediated through physiological imbalance to be $194 per case. The results showed that if the alert system was highly accurate (Se = 0.99/Sp = 0.99), and the prevalence of physiological imbalance was 30 %, the value of information provided from the system is $19 per cow-year. In case the prevalence is 5 % or 50 %, the value of information is $3 and $13, respectively. These estimates for the value do not cover the capital costs and operational costs of the alert system. This study furthermore clearly demonstrated that in order to estimate the value of information correctly, it is important to consider that drenching all cows and not drenching any of the cows are the two relevant alternative options in the absence of the alert system. In conclusion, the decision tree and sensitivity analysis developed in this study show that final economic results are highly variable to the prevalence of physiological imbalance and highest at an intermediate prevalence. Other relevant factors are the costs associated with drenching and the cost associated with treating false positives and not treating false negatives. In addition, this study highlights the benefits of simulation to pinpoint where additional information is needed to further quantify the economic value and required accuracy of an indication-based intervention system.