Validation of mid-infrared spectrometry in milk for predicting body energy status in Holstein-Friesian cows

Energy balance of dairy cows predicted by mid-infrared spectra data of milk using Bayesian approaches

Information on dry matter intake (DMI) and energy balance (EB) at the animal and herd level is important for management and breeding decisions. However, routine recording of these traits at commercial farms can be challenging and costly. Fourier-transform mid-infrared (FT-MIR) spectroscopy is a noninvasive technique applicable to a large cohort of animals that is routinely used to analyze milk components and is convenient for predicting complex phenotypes that are typically difficult and expensive to obtain on a large scale. We aimed to develop prediction models for EB and use the predicted phenotypes for genetic analysis. First, we assessed prediction equations using 4,485 phenotypic records from 167 Holstein cows from an experimental station. The phenotypes available were body weight (BW), milk yield (MY) and milk components, weekly-averaged DMI, and FT-MIR data from all milk samples available. We implemented mixed models with Bayesian approaches and assessed them through 50 randomized replicates of a 5-fold cross-validation. Second, we used the best prediction models to obtain predicted phenotypes of EB (EBp) and DMI (DMIp) on 5 commercial farms with 2,365 phenotypic records of MY, milk components and FT-MIR data, and BW from 1,441 Holstein cows. Third, we performed a GWAS and estimated heritability and genetic correlations for energy content in milk (EnM), BW, DMIp, and EBp using the genomic information available on the cows from commercial farms. The highest correlation between the predicted and observed phenotype (ry,y^) was obtained with DMI (0.88) and EB (0.86), while predicting BW was, as anticipated, more challenging (0.69). In our study, models that included FT-MIR information performed better than models without spectra information in the 3 traits analyzed, with increments in prediction correlation ranging from 5% to 10%. For the predicted phenotypes calculated by the prediction equations and data from the commercial farms the heritability ranged between 0.11 and 0.16 for EnM, DMIp and EBp, and 0.42 for BW. The genetic correlation between EnM and BW was -0.17, with DMIp was 0.40 and with EBp was -0.39. From the GWAS, we detected one significant QTL region for EnM, and 3 for BW, but none for DMIp and EBp. The results obtained in our study support previous evidence that FT-MIR information from milk samples contribute to improve the prediction equations for DMI, BW, and EB, and these predicted phenotypes may be used for herd management and contribute to the breeding strategy for improving cow performance.

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.

Association between Days Open and Parity, Calving Season or Milk Spectral Data

Milk spectral data on 2118 cows from nine herds located in northern China were used to access the association of days open (DO). Meanwhile, the parity and calving season of dairy cows were also studied to characterize the difference in DO between groups of these two cow-level factors. The result of the linear mixed-effects model revealed that no significant differences were observed between the parity groups. However, a significant difference in DO exists between calving season groups. The interaction between parity and calving season presented that primiparous cows always exhibit lower DO among all calving season groups, and the variation in DO among parity groups was especially clearer in winter. Survival analysis revealed that the difference in DO between calving season groups might be caused by the different P/AI at the first TAI. In addition, the summer group had a higher chance of conception in the subsequent services than other groups, implying that the micro-environment featured by season played a critical role in P/AI. A weak linkage between DO and wavenumbers ranging in the mid-infrared region was detected. In summary, our study revealed that the calving season of dairy cows can be used to optimize the reproduction management. The potential application of mid-infrared spectroscopy in dairy cows needs to be further developed.

Comparing human milk macronutrients measured using analyzers based on mid-infrared spectroscopy and ultrasound and the application of machine learning in data fitting

OBJECTIVE

Fat, carbohydrates (mainly lactose) and protein in breast milk all provide indispensable benefits for the growth of newborns. The only source of nutrition in early infancy is breast milk, so the energy of breast milk is also crucial to the growth of infants. Some macronutrients composition in human breast milk varies greatly, which could affect its nutritional fulfillment to preterm infant needs. Therefore, rapid analysis of macronutrients (including lactose, fat and protein) and milk energy in breast milk is of clinical importance. This study compared the macronutrients results of a mid-infrared (MIR) analyzer and an ultrasound-based breast milk analyzer and unified the results by machine learning.

METHODS

This cross-sectional study included breastfeeding mothers aged 22-40 enrolled between November 2019 and February 2021. Breast milk samples (n = 546) were collected from 244 mothers (from Day 1 to Day 1086 postpartum). A MIR milk analyzer (BETTERREN Co., HMIR-05, SH, CHINA) and an ultrasonic milk analyzer (Honɡyanɡ Co,. HMA 3000, Hebei, CHINA) were used to determine the human milk macronutrient composition. A total of 465 samples completed the tests in both analyzers. The results of the ultrasonic method were mathematically converted using machine learning, while the Bland-Altman method was used to determine the limits of agreement (LOA) between the adjusted results of the ultrasonic method and MIR results.

RESULTS

The MIR and ultrasonic milk analyzer results were significantly different. The protein, fat, and energy determined using the MIR method were higher than those determined by the ultrasonic method, while lactose determined by the MIR method were lower (all p < 0.05). The consistency between the measured MIR and the adjusted ultrasound values was evaluated using the Bland-Altman analysis and the scatter diagram was generated to calculate the 95% LOA. After adjustments, 93.96% protein points (436 out of 465), 94.41% fat points (439 out of 465), 95.91% lactose points (446 out of 465) and 94.62% energy points (440 out of 465) were within the LOA range. The 95% LOA of protein, fat, lactose and energy were - 0.6 to 0.6 g/dl, -0.92 to 0.92 g/dl, -0.88 to 0.88 g/dl and - 40.2 to 40.4 kj/dl, respectively and clinically acceptable. The adjusted ultrasonic results were consistent with the MIR results, and LOA results were high (close to 95%).

CONCLUSIONS

While the results of the breast milk rapid analyzers using the two methods varied significantly, they could still be considered comparable after data adjustments using linear regression algorithm in machine learning. Machine learning methods can play a role in data fitting using different analyzers.

Comparison of Single-Breed and Multi-Breed Training Populations for Infrared Predictions of Novel Phenotypes in Holstein Cows

In general, Fourier-transform infrared (FTIR) predictions are developed using a single-breed population split into a training and a validation set. However, using populations formed of different breeds is an attractive way to design cross-validation scenarios aimed at increasing prediction for difficult-to-measure traits in the dairy industry. This study aimed to evaluate the potential of FTIR prediction using training set combining specialized and dual-purpose dairy breeds to predict different phenotypes divergent in terms of biological meaning, variability, and heritability, such as body condition score (BCS), serum β-hydroxybutyrate (BHB), and kappa casein (k-CN) in the major cattle breed, i.e., Holstein-Friesian. Data were obtained from specialized dairy breeds: Holstein (468 cows) and Brown Swiss (657 cows), and dual-purpose breeds: Simmental (157 cows), Alpine Grey (75 cows), and Rendena (104 cows), giving a total of 1461 cows from 41 multi-breed dairy herds. The FTIR prediction model was developed using a gradient boosting machine (GBM), and predictive ability for the target phenotype in Holstein cows was assessed using different cross-validation (CV) strategies: a within-breed scenario using 10-fold cross-validation, for which the Holstein population was randomly split into 10 folds, one for validation and the remaining nine for training (10-fold_HO); an across-breed scenario (BS_HO) where the Brown Swiss cows were used as the training set and the Holstein cows as the validation set; a specialized multi-breed scenario (BS+HO_10-fold), where the entire Brown Swiss and Holstein populations were combined then split into 10 folds, and a multi-breed scenario (Multi-breed), where the training set comprised specialized (Holstein and Brown Swiss) and dual-purpose (Simmental, Alpine Grey, and Rendena) dairy cows, combined with nine folds of the Holstein cows. Lastly a Multi-breed CV2 scenario was implemented, assuming the same number of records as the reference scenario and using the same proportions as the multi-breed. Within-Holstein, FTIR predictions had a predictive ability of 0.63 for BCS, 0.81 for BHB, and 0.80 for k-CN. Using a specific breed (Brown Swiss) as the training set for prediction in the Holstein population reduced the prediction accuracy by 10% for BCS, 7% for BHB, and 11% for k-CN. Notably, the combination of Holstein and Brown Swiss cows in the training set increased the predictive ability of the model by 6%, which was 0.66 for BCS, 0.85 for BHB, and 0.87 for k-CN. Using multiple specialized and dual-purpose animals in the training set outperforms the 10-fold_HO (standard) approach, with an increase in predictive ability of 8% for BCS, 7% for BHB, and 10% for k-CN. When the Multi-breed CV2 was implemented, no improvement was observed. Our findings suggest that FTIR prediction of different phenotypes in the Holstein breed can be improved by including different specialized and dual-purpose breeds in the training population. Our study also shows that predictive ability is enhanced when the size of the training population and the phenotypic variability are increased.

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.

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.

Fourier transform infrared spectroscopy of milk samples as a tool to estimate energy balance, energy- and dry matter intake in lactating dairy cows

The objective of the study was to evaluate the potential of Fourier transform infrared spectroscopy (FTIR) analysis of milk samples to predict body energy status and related traits (energy balance (EB), dry matter intake (DMI) and efficient energy intake (EEI)) in lactating dairy cows. The data included 2371 milk samples from 63 Norwegian Red dairy cows collected during the first 105 days in milk (DIM). To predict the body energy status traits, calibration models were developed using Partial Least Squares Regression (PLSR). Calibration models were established using split-sample (leave-one cow-out) cross-validation approach and validated using an external test set. The PLSR method was implemented using just the FTIR spectra or using the FTIR together with milk yield (MY) or concentrate intake (CONCTR) as predictors of traits. Analyses were conducted for the entire first 105 DIM and separately for the two lactation periods: 5 ≤ DIM ≤ 55 and 55 < DIM ≤ 105. To test the models, an external validation using an independent test set was performed. Predictions depending on the parity (1st, 2nd and 3rd-to 6th parities) in early lactation were also investigated. Accuracy of prediction (r) for both cross-validation and external test set was defined as the correlation between the predicted and observed values for body energy status traits. Analyzing FTIR in combination with MY by PLSR, resulted in relatively high r-values to estimate EB (r = 0.63), DMI (r = 0.83), EEI (r = 0.84) using an external validation. Only moderate correlations between FTIR spectra and traits like EB, EEI and dry matter intake (DMI) have so far been published. Our hypothesis was that improvements in the FTIR predictions of EB, EEI and DMI can be obtained by (1) stratification into different stages of lactations and different parities, or (2) by adding additional information on milking and feeding traits. Stratification of the lactation stages improved predictions compared with the analyses including all data 5 ≤ DIM ≤105. The accuracy was improved if additional data (MY or CONCTR) were included in the prediction model. Furthermore, stratification into parity groups, improved the predictions of body energy status. Our results show that FTIR spectral data combined with MY or CONCTR can be used to obtain improved estimation of body energy status compared to only using the FTIR spectra in Norwegian Red dairy cattle. The best prediction results were achieved using FTIR spectra together with MY for early lactation. The results obtained in the study suggest that the modeling approach used in this paper can be considered as a viable method for predicting an individual cow's energy status.

Genetic mechanisms underlying feed utilization and implementation of genomic selection for improved feed efficiency in dairy cattle

The economic importance of genetically improving feed efficiency has been recognized by cattle producers worldwide. It has the potential to considerably reduce costs, minimize environmental impact, optimize land and resource use efficiency, and improve the overall cattle industry’s profitability. Feed efficiency is a genetically complex trait that can be described as units of product output (e.g., milk yield) per unit of feed input. The main objective of this review paper is to present an overview of the main genetic and physiological mechanisms underlying feed utilization in ruminants and the process towards implementation of genomic selection for feed efficiency in dairy cattle. In summary, feed efficiency can be improved via numerous metabolic pathways and biological mechanisms through genetic selection. Various studies have indicated that feed efficiency is heritable, and genomic selection can be successfully implemented in dairy cattle with a large enough training population. In this context, some organizations have worked collaboratively to do research and develop training populations for successful implementation of joint international genomic evaluations. The integration of “-omics” technologies, further investments in high-throughput phenotyping, and identification of novel indicator traits will also be paramount in maximizing the rates of genetic progress for feed efficiency in dairy cattle worldwide.

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.

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

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

Genetic correlations between energy status indicator traits and female fertility in primiparous Nordic Red Dairy cattle

Inclusion of feed efficiency traits into the dairy cattle breeding programmes will require considering early lactation energy status to avoid deterioration in health and fertility of dairy cows. In this regard, energy status indicator (ESI) traits, for example, blood metabolites or milk fatty acids (FAs), are of interest. These indicators can be predicted from routine milk samples by mid-IR reflectance spectroscopy (MIR). In this study, we estimated genetic variation in ESI traits and their genetic correlation with female fertility in early lactation. The data consisted of 37 424 primiparous Nordic Red Dairy cows with milk test-day records between 8 and 91 days in milk (DIM). Routine test-day milk samples were analysed by MIR using previously developed calibration equations for blood plasma non-esterified FA (NEFA), milk FAs, milk beta-hydroxybutyrate (BHB) and milk acetone concentrations. Six ESI traits were considered and included: plasma NEFA concentration (mmol/l) either predicted by multiple linear regression including DIM, milk fat to protein ratio (FPR) and FAs C10:0, C14:0, C18:1 cis-9, C14:0 * C18:1 cis-9 (NEFAFA) or directly from milk MIR spectra (NEFAMIR), C18:1 cis-9 (g/100 ml milk), FPR, BHB (mmol/l milk) and acetone (mmol/l milk). The interval from calving to first insemination (ICF) was considered as the fertility trait. Data were analysed using linear mixed models. Heritability estimates varied during the first three lactation months from 0.13 to 0.19, 0.10 to 0.17, 0.09 to 0.14, 0.07 to 0.10, 0.13 to 0.17 and 0.13 to 0.18 for NEFAMIR, NEFAFA, C18:1 cis-9, FPR, milk BHB and acetone, respectively. Genetic correlations between all ESI traits and ICF were from 0.18 to 0.40 in the first lactation period (8 to 35 DIM), in general somewhat lower (0.03 to 0.43) in the second period (36 to 63 DIM) and decreased clearly (-0.02 to 0.19) in the third period (64 to 91 DIM). Our results indicate that genetic variation in energy status of cows in early lactation can be determined using MIR-predicted indicators. In addition, the markedly lower genetic correlation between ESI traits and fertility in the third lactation month indicated that energy status should be determined from the first test-day milk samples during the first 2 months of lactation.

Predicting milk fatty acids and energy balance of dairy cows in Australia using milk mid-infrared spectroscopy

Mid-infrared spectroscopy (MIRS) is traditionally used for analysing milk fat, protein and lactose concentrations in dairy production, but there is growing interest in using it to predict difficult, or expensive-to-measure, phenotypes on a large scale. The resulting prediction equations can be applied to MIRS data from commercial herd-testing, to facilitate management and feeding decisions, or for genomic selection purposes. We investigated the ability of MIRS of milk samples to predict milk fatty acids (FAs) and energy balance (EB) of dairy cows in Australia. Data from 240 Holstein lactating cows that were part of two 32-day experiments, were used. Milk FAs were measured twice during the experimental period. Prediction models were developed using partial least-square regression with a 10-fold cross-validation. Measures of prediction accuracy included the coefficient of determination (R2cv) and root mean-square error. Milk FAs with a chain length of ≤16 were accurately predicted (0.89 ≤ R2cv ≤ 0.95), while prediction accuracy for FAs with a chain length of ≥17 was slightly lower (0.72 ≤ R2cv ≤ 0.82). The accuracy of the model prediction was moderate for EB, with the value of R2cv of 0.48. In conclusion, the ability of MIRS to predict milk FAs was high, while EB was moderately predicted. A larger dataset is needed to improve the accuracy and the robustness of the prediction models.

Energy profiling of dairy cows from routine milk mid-infrared analysis

The balance of body energy within and across lactations can have health and fertility consequences for the dairy cow. This study aimed to create a large calibration data set of dairy cow body energy traits across the cow's productive life, with concurrent milk mid-infrared (MIR) spectral data, to generate a prediction tool for use in commercial dairy herds. Detailed phenotypic data from 1,101 Holstein Friesian cows from the Langhill research herd (SRUC, Scotland) were used to generate energy balance (EB) and effective energy intake (EI), both in megajoules per day. Pretreatment of spectral data involved standardization to account for drift over time and machine. Body energy estimates were aligned with their spectral data to generate a prediction of these traits based on milk MIR spectroscopy. After data edits, partial least squares analysis generated prediction equations with a coefficient of determination from split sample 10-fold cross validation of 0.77 and 0.75 for EB and EI, respectively. These prediction equations were applied to national milk MIR spectra on over 11 million animal test dates (January 2013 to December 2016) from 4,453 farms. The predictions generated from these were subject to phenotypic analyses with a fixed regression model highlighting differences between the main dairy breeds in terms of energy traits. Genetic analyses generated heritability estimates for EB and EI ranging from 0.12 to 0.17 and 0.13 to 0.15, respectively. This study shows that MIR-based predictions from routinely collected national data can be used to generate predictions of dairy cow energy turnover profiles for both animal management and genetic improvement of such difficult and expensive-to-record traits.

Classifying the fertility of dairy cows using milk mid-infrared spectroscopy

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

Body and milk traits as indicators of dairy cow energy status in early lactation

The inclusion of feed intake and efficiency traits in dairy cow breeding goals can lead to increased risk of metabolic stress. An easy and inexpensive way to monitor postpartum energy status (ES) of cows is therefore needed. Cows' ES can be estimated by calculating the energy balance from energy intake and output and predicted by indicator traits such as change in body weight (ΔBW), change in body condition score (ΔBCS), milk fat:protein ratio (FPR), or milk fatty acid (FA) composition. In this study, we used blood plasma nonesterified fatty acids (NEFA) concentration as a biomarker for ES. We determined associations between NEFA concentration and ES indicators and evaluated the usefulness of body and milk traits alone, or together, in predicting ES of the cow. Data were collected from 2 research herds during 2013 to 2016 and included 137 Nordic Red dairy cows, all of which had a first lactation and 59 of which also had a second lactation. The data included daily body weight, milk yield, and feed intake and monthly BCS. Plasma samples for NEFA were collected twice in lactation wk 2 and 3 and once in wk 20. Milk samples for analysis of fat, protein, lactose, and FA concentrations were taken on the blood sampling days. Plasma NEFA concentration was higher in lactation wk 2 and 3 than in wk 20 (0.56 ± 0.30, 0.43 ± 0.22, and 0.13 ± 0.06 mmol/L, respectively; all means ± standard deviation). Among individual indicators, C18:1 cis-9 and the sum of C18:1 in milk had the highest correlations (r = 0.73) with NEFA. Seven multiple linear regression models for NEFA prediction were developed using stepwise selection. Of the models that included milk traits (other than milk FA) as well as body traits, the best fit was achieved by a model with milk yield, FPR, ΔBW, ΔBCS, FPR × ΔBW, and days in milk. The model resulted in a cross-validation coefficient of determination (R²cv) of 0.51 and a root mean squared error (RMSE) of 0.196 mmol/L. When only milk FA concentrations were considered in the model, NEFA prediction was more accurate using measurements from evening milk than from morning milk (R²cv = 0.61 vs. 0.53). The best model with milk traits contained FPR, C10:0, C14:0, C18:1 cis-9, C18:1 cis-9 × C14:0, and days in milk (R²cv = 0.62; RMSE = 0.177 mmol/L). The most advanced model using both milk and body traits gave a slightly better fit than the model with only milk traits (R²cv = 0.63; RMSE = 0.176 mmol/L). Our findings indicate that ES of cows in early lactation can be monitored with moderately high accuracy by routine milk measurements.

The potential of Fourier transform infrared spectroscopy of milk samples to predict energy intake and efficiency in dairy cows

Knowledge of animal-level and herd-level energy intake, energy balance, and feed efficiency affect day-to-day herd management strategies; information on these traits at an individual animal level is also useful in animal breeding programs. A paucity of data (especially at the individual cow level), of feed intake in particular, hinders the inclusion of such attributes in herd management decision-support tools and breeding programs. Dairy producers have access to an individual cow milk sample at least once daily during lactation, and consequently any low-cost phenotyping strategy should consider exploiting measureable properties in this biological sample, reflecting the physiological status and performance of the cow. Infrared spectroscopy is the study of the interaction of an electromagnetic wave with matter and it is used globally to predict milk quality parameters on routinely acquired individual cow milk samples and bulk tank samples. Thus, exploiting infrared spectroscopy in next-generation phenotyping will ensure potentially rapid application globally with a negligible additional implementation cost as the infrastructure already exists. Fourier-transform infrared spectroscopy (FTIRS) analysis is already used to predict milk fat and protein concentrations, the ratio of which has been proposed as an indicator of energy balance. Milk FTIRS is also able to predict the concentration of various fatty acids in milk, the composition of which is known to change when body tissue is mobilized; that is, when the cow is in negative energy balance. Energy balance is mathematically very similar to residual energy intake (REI), a suggested measure of feed efficiency. Therefore, the prediction of energy intake, energy balance, and feed efficiency (i.e., REI) from milk FTIRS seems logical. In fact, the accuracy of predicting (i.e., correlation between predicted and actual values; root mean square error in parentheses) energy intake, energy balance, and REI from milk FTIRS in dairy cows was 0.88 (20.0MJ), 0.78 (18.6MJ), and 0.63 (22.0MJ), respectively, based on cross-validation. These studies, however, are limited to results from one research group based on data from 2 contrasting production systems in the United Kingdom and Ireland and would need to be replicated, especially in a range of production systems because the prediction equations are not accurate when the variability used in validation is not represented in the calibration data set. Heritable genetic variation exists for all predicted traits. Phenotypic differences in energy intake also exists among animals stratified based on genetic merit for energy intake predicted from milk FTIRS, substantiating the usefulness of such FTIR-predicted phenotypes not only for day-to-day herd management, but also as part of a breeding strategy to improve cow performance.