Milk Mid-infrared Spectroscopy Research Articles

Transferability across countries of equations developed using milk mid-infrared spectroscopy to estimate daily body condition score change in dairy cows

Routine milk samples are commonly subjected to spectroscopic analysis within the mid-infrared (MIR) region of the electromagnetic spectrum to estimate macro-constituents of milk like fat, protein, lactose, and urea content. These spectra, however, can also be used to predict other traits, such as daily body condition score (BCS) change. The objective of the present study was to assess the transferability across countries of equations to predict daily body condition score change (ΔBCS) developed using milk MIR data collected in Ireland and in Canada. Body condition was scored on a scale from 1 (emaciated) to 5 (obese) in both countries. A total of 347,254 BCS records from 80,400 Canadian cows were available along with 73,193 BCS records from 6,572 Irish cows. Partial least squares regression (PLSR) and neural networks (NN) were separately used to predict daily ΔBCS. Two scenarios were studied 1) using Canadian and Irish data combined as the calibration data set to predict daily ΔBCS in Canada and in Ireland separately, and 2) Canadian and Irish data used separately to predict daily ΔBCS in each country separately. These prediction methods were applied to data with and without pretreatment (i.e., first derivative of the spectrum) as well as with and without standardizing daily ΔBCS across countries. For all the scenarios investigated, the correlation between actual and predicted daily ΔBCS when calibrated and validated (using cross-validation) in the same country ranged from 0.92 to 0.94, and from 0.85 to 0.87 for the Canadian and Irish data sets, respectively. When the data from Canada and Ireland were combined in the calibration process to predict daily ΔBCS, the correlations between actual and predicted ΔBCS were ≥ 0.90 and ≥ 0.80 for Canadian and Irish daily ΔBCS, respectively indicating no improvement in predictive ability. Predictive performance when calibrated using just Canadian data and validated using just Irish data was poor, and vice versa. Nonetheless, when developing equations for a country for which a limited database (i.e., 100 records) of gold standard and MIR data were available, predictive performance improved when the limited database was supplemented with the large data set from the other country. In general, for some of the investigated scenarios, standardizing the daily ΔBCS data within country before undertaking the calibration improved prediction accuracy. In conclusion, the benefit of merging data from different countries, at least based on the trait (i.e., daily ΔBCS) and countries (i.e., Ireland and Canada) considered in the present study were limited and, in cases, counter-productive.

Identification and quantification of goat milk adulteration using mid-infrared spectroscopy and chemometrics

The fraudulent adulteration of goat milk with cheaper and more available milk of other species such as cow milk is occurrence. The aims of the present study were to investigate the effect of goat milk adulteration with cow milk on the mid-infrared (MIR) spectrum and further evaluate the potential of MIR spectroscopy to identify and quantify the goat milk adulterated. Goat milk was adulterated with cow milk at 5 different levels including 10%, 20%, 30%, 40%, and 50%. Statistical analysis showed that the adulteration had significant effect on the majority of the spectral wavenumbers. Then, the spectrum was preprocessed with standard normal variate (SNV), multiplicative scattering correction (MSC), Savitzky-Golay smoothing (SG), SG plus SNV, and SG plus MSC, and partial least squares discriminant analysis (PLS-DA) and partial least squares regression (PLSR) were used to establish classification and regression models, respectively. PLS-DA models obtained good results with all the sensitivity and specificity over 0.96 in the cross-validation set. Regression models using raw spectrum obtained the best result, with coefficient of determination (R2), root mean square error (RMSE), and the ratio of performance to deviation (RPD) of cross-validation set were 0.98, 2.01, and 8.49, respectively. The results preliminarily indicate that the MIR spectroscopy is an effective technique to detect the goat milk adulteration with cow milk. In future, milk samples from different origins and different breeds of goats and cows should be collected, and more sophisticated adulteration at low levels should be further studied to explore the potential and effectiveness of milk mid-infrared spectroscopy and chemometrics.

Using milk mid-infrared spectroscopy to estimate cow-level nitrogen efficiency metrics

Minimizing pollution from the dairy sector is paramount; one potential cause of such pollution is excess nitrogen. Nitrogen pollution contributes to a deterioration in water quality as well as an increase in both eutrophication and greenhouse gases. It is therefore essential to minimize the loss of nitrogen from the sector, including excretion from the cow. Breeding programs are one potential strategy to improve the efficiency with which nitrogen is used by dairy cows but relies on routine access to individual cow information on how efficiently each cows uses the nitrogen it ingests. A total of 3,497 test-day records for individual cow nitrogen efficiency metrics along with milk yield and the associated milk spectra were used to investigate the ability of milk infrared spectral data to predict these nitrogen traits; both traditional partial least squares regression and neural networks were used in the prediction process. The data originated from 4 farms across 11 years. The nitrogen traits investigated were nitrogen intake, nitrogen use efficiency, and nitrogen balance. Both nitrogen use efficiency and nitrogen balance were calculated considering nitrogen intake, nitrogen in milk, nitrogen in the conceptus, nitrogen used for the growth, nitrogen stored in the reserves, and nitrogen mobilized from the reserves. Irrespective of the nitrogen-related trait being investigated, the best prediction from 4-fold cross-validation were achieved using neural networks that considered both the morning and evening milk spectra along with milk yield, parity, and days in milk in the prediction process. The coefficient of determination in the cross-validation was 0.61, 0.74, and 0.58 for nitrogen intake, nitrogen use efficiency, and nitrogen balance, respectively. In a separate series of validation approaches, the calibration and validation was stratified by herd (n = 4) and separately by year. For these scenarios, partial least squares regression generated more accurate predictions compared with neural networks; the coefficient of determination was always lower than 0.29 and 0.60 when validation was stratified by herd and year, respectively. Therefore, if the variability of the data being predicted in the validation data sets is similar to that in the data used to develop the predictions, then nitrogen-related traits can be predicted with reasonable accuracy. In contrast, where the variability of the data that exists in the validation data set is poorly represented in the calibration data set, then poor predictions will ensue.

Is it possible to predict the methane emission intensity of Swedish dairy cows from milk spectra?

Emissions of greenhouse gases (GHG), especially methane (CH4), from dairy production have received much research attention in the past 15 years, with the main focus being to identify factors affecting CH4 production and measures to reduce CH4 emissions from dairy cows. However, measurement of CH4 production by dairy cows in commercial herds is time-consuming and requires expensive equipment, so there is a need to find alternative ways to estimate individual and herd CH4 emissions. Regular milk analyses are performed for many cows, so data from mid-infrared spectroscopy (MIRS) on individual milk samples could perhaps be utilised to predict CH4 emissions intensity (MI, CH4 g/kg milk production). This study investigated the potential and limitations of predicting individual MI by integrating data from CH4 measurements made by an infrared sniffer (IRS) and milk MIRS data taken from fortnightly morning milkings during the full lactation records of 37 multiparous cows. Partial least square regression was used to create prediction models for six-week lactation sub-periods and for full lactations, which were validated using leave-one-cow-out cross-validation. Coefficient of determination in predicting MI was low, indicating that the method is not suitable for predicting variations in individual MI, although it should further be evaluated at herd level.

Using machine learning methods to predict dry matter intake from milk mid-infrared spectroscopy data on Swedish dairy cattle.

In this research communication we compare three different approaches for developing dry matter intake (DMI) prediction models based on milk mid-infrared spectra (MIRS), using data collected from a research herd over five years. In dairy production, knowledge of individual DMI could be important and useful, but DMI can be difficult and expensive to measure on most commercial farms as cows are commonly group-fed. Instead, this parameter is often estimated based on the age, body weight, stage of lactation and body condition score of the cow. Recently, milk MIRS have also been used as a tool to estimate DMI. There are different methods available to create prediction models from large datasets. The main data used were total DMI calculated as a 3-d average, coupled with milk MIRS data available fortnightly. Data on milk yield and lactation stage parameters were also available for each animal. We compared the performance of three prediction approaches: partial least-squares regression, support vector machine regression and random forest regression. The full milk MIRS alone gave low to moderate prediction accuracy (R2 = 0.07-0.40), regardless of prediction modelling approach. Adding more variables to the model improved R2 and decreased the prediction error. Overall, partial least-squares regression proved to be the best method for predicting DMI from milk MIRS data, while MIRS data together with milk yield and concentrate DMI at 3-30 d in milk provided good prediction accuracy (R2 = 0.52-0.65) regardless of the prediction tool used.

Predicting methane emission in Canadian Holstein dairy cattle using milk mid-infrared reflectance spectroscopy and other commonly available predictors via artificial neural networks.

Interest in reducing eructed CH4 is growing, but measuring CH4 emissions is expensive and difficult in large populations. In this study, we investigated the effectiveness of milk mid-infrared spectroscopy (MIRS) data to predict CH4 emission in lactating Canadian Holstein cows. A total of 181 weekly average CH4 records from 158 Canadian cows and 217 records from 44 Danish cows were used. For each milk spectra record, the corresponding weekly average CH4 emission (g/d), test-day milk yield (MY, kg/d), fat yield (FY, g/d), and protein yield (PY, g/d) were available. The weekly average CH4 emission was predicted using various artificial neural networks (ANN), partial least squares regression, and different sets of predictors. The ANN architectures consisted of 3 training algorithms, 1 to 10 neurons with hyperbolic tangent activation function in the hidden layer, and 1 neuron with linear (purine) activation function in the hidden layer. Random cross-validation was used to compared the predictor sets: MY (set 1); FY (set 2); PY (set 3); MY and FY (set 4); MY and PY (set 5); MY, FY, and PY (set 6); MIRS (set 7); and MY, FY, PY, and MIRS (set 8). All predictor sets also included age at calving and days in milk, in addition to country, season of calving, and lactation number as categorical effects. Using only MY (set 1), the predictive accuracy (r) ranged from 0.245 to 0.457 and the root mean square error (RMSE) ranged from 87.28 to 99.39 across all prediction models and validation sets. Replacing MY with FY (set 2; r = 0.288-0.491; RMSE = 85.94-98.04) improved the predictive accuracy, but using PY (set 3; r = 0.260-0.468; RMSE = 86.95-98.47) did not. Adding FY to MY (set 4; r = 0.272-0.469; RMSE = 87.21-100.76) led to a negligible improvement compared with sets 1 and 3, but it slightly decreased accuracy compared with set 2. Adding PY to MY (set 5; r = 0.250-0.451; RMSE = 87.66-100.94) did not improve prediction ability. Combining MY, FY, and PY (set 6; r = 0.252-0.455; RMSE = 87.74-101.93) yielded accuracy slightly lower than sets 2 and 3. Using only MIRS data (set 7; r = 0.586-0.717; RMSE = 69.09-96.20) resulted in superior accuracy compared with all previous sets. Finally, the combination of MIRS data with MY, FY, and PY (set 8; r = 0.590-0.727; RMSE = 68.02-87.78) yielded similar accuracy to set 7. Overall, sets including the MIRS data yielded significantly better predictions than the other sets. To assess the predictive ability in a new unseen herd, a limited block cross-validation was performed using 20 cows in the same Canadian herd, which yielded r = 0.229 and RMSE = 154.44, which were clearly much worse than the average r = 0.704 and RMSE = 70.83 when predictions were made by random cross-validation. These results warrant further investigation when more data become available to allow for a more comprehensive block cross-validation before applying the calibrated models for large-scale prediction of CH4 emissions.

Phenotypic and genetic characterization of the occurrence of noncoagulating milk in dairy sheep

Milk coagulation ability is of central importance for the sheep dairy industry because almost all sheep milk is destined for cheese processing. The occurrence of milk with impaired coagulation properties is an obstacle to cheese processing and, in turn, to the profitability of the dairy companies. In this work, we investigated the causes of noncoagulation of sheep milk; specifically, we studied the effect of milk physicochemical properties on milk coagulation status [coagulating and noncoagulating (NC) milk samples, which do or do not coagulate within 30 min, respectively], and whether mid-infrared spectroscopy (MIR) could be used to assess variability in coagulation status. We also investigated the genetic background of milk coagulation ability. Individual milk samples were collected from 996 Sarda ewes farmed in 47 flocks located in Sardinia (Italy). Considered traits were daily milk yield, milk composition traits, and milk coagulation properties (rennet coagulation time, curd firming time, and curd firmness), and MIR spectra were acquired. About 9% of samples did not coagulate within 30 min. A logistic regression approach was used to test the effect of milk-related traits on milk coagulation status. A principal component (PC) analysis was carried out on the milk MIR spectra, and PC scores were then used as covariates in a logistic regression model to assess their relationship with milk coagulation status. Results of the present work demonstrated that the probability of having NC samples increases as milk contents of proteins and chlorides and somatic cell score increase. The analysis of PC extracted from milk spectra that influenced coagulation status highlighted key regions associated with lactose and protein concentrations, and others not associated with routinely collected milk composition traits. These results suggest that the occurrence of NC is mostly related to damage of the epithelium secretory mammary cells, which occurs with the advancement of a lactation or due to unhealthy mammary gland status. Genetic analysis of milk coagulation status and of the extracted PC confirmed the genetic background of the milk coagulability of sheep milk.

Identifying Health Status in Grazing Dairy Cows from Milk Mid-Infrared Spectroscopy by Using Machine Learning Methods.

Simple SummaryDiseases in dairy livestock farming can lead to important economic losses. Several studies have been conducted to identify illness such as lameness by using MIR spectrometry data and relying on partial least squares discriminant analysis. However, this method suffers some limitations. In this study, random forest, support vector machine, neural network, convolutional neural network and ensemble models were used to test the feasibility of identifying cow sickness among 1909 milk sample MIR spectra from Holstein-Friesian, Jersey and Holstein-Friesian × Jersey crossbreed cows. The results obtained show that it is possible to identify a health problem with a reasonable level of accuracy using a neural network.The early detection of health problems in dairy cattle is crucial to reduce economic losses. Mid-infrared (MIR) spectrometry has been used for identifying the composition of cow milk in routine tests. As such, it is a potential tool to detect diseases at an early stage. Partial least squares discriminant analysis (PLS-DA) has been widely applied to identify illness such as lameness by using MIR spectrometry data. However, this method suffers some limitations. In this study, a series of machine learning techniques—random forest, support vector machine, neural network (NN), convolutional neural network and ensemble models—were used to test the feasibility of identifying cow sickness from 1909 milk sample MIR spectra from Holstein-Friesian, Jersey and crossbreed cows under grazing conditions. PLS-DA was also performed to compare the results. The sick cow records had a time window of 21 days before and 7 days after the milk sample was analysed. NN showed a sensitivity of 61.74%, specificity of 97% and positive predicted value (PPV) of nearly 60%. Although the sensitivity of the PLS-DA was slightly higher than NN (65.6%), the specificity and PPV were lower (79.59% and 15.25%, respectively). This indicates that by using NN, it is possible to identify a health problem with a reasonable level of accuracy.

Mining data from milk mid-infrared spectroscopy and animal characteristics to improve the prediction of dairy cow's liveweight using feature selection algorithms based on partial least squares and Elastic Net regressions

Body weight (BW) of dairy cows is relevant for breeding programs and farm management to assess the maintenance requirements, reproduction performance, or health status of cow. Currently, it is still difficult to follow BW changes of individual cows routinely in large herds. Combined with animal characteristics, milk mid-infrared (MIR) spectrum was proposed as an additional source of information to predict BW under the framework of dairy herd improvement (DHI) programs. However, the presence of less informative variables in the prediction equation could impact negatively its robustness. This research aims to improve the robustness of BW regression models by applying a feature selection before modeling. A total of 5,920 BW records composed of animal characteristics and milk MIR spectrum were collected from Holstein cows. Three feature selection algorithms were applied to select the most informative variables: partial least squares regression (PLS) combined with sum of ranking difference (PLS-SRD), PLS combined with uninformative variables elimination (PLS-UVE), and the output of Elastic Net regression (EN). Four herd independent validation sets and the corresponding remained calibration datasets having on average 163 and 1,708 records, respectively, were used to develop models using PLS or EN approaches. Ten-fold cross-validation was conducted to parametrize each model. Parity, days in milk (DIM), milk yield (MY), and two MIR spectral points were selected as relevant variables to predict BW. PLS (root mean square error of validation, RMSEp = 60 kg) and EN (RMSEp = 60 kg) regressions employing these 5 predictors were more robust than the models developed without MIR or using MIR without feature selection. The EN models had a cross-validation root mean square error of around 53 kg. The 2 MIR points explained up to 4.20% variation in predicting BW. The RMSE of validation sets using another brand of spectrometer were around 64 kg. This study confirms the possibility to predict an indicator of BW from animal characteristics and MIR variables. The variable selection procedures improved the model's robustness and transferability. The accuracy of BW prediction seems to be sufficient to provide useful information for breeding program and farm management decisions under a DHI framework.

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 (R2) 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 R2 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 R2 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.

The use of milk mid-infrared spectroscopy to improve genomic prediction accuracy of serum biomarkers

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

Prediction of pregnancy state from milk mid-infrared (MIR) spectroscopy in dairy cows

Prediction of pregnancy state from milk mid-infrared (MIR) spectroscopy in dairy cows

Predicting the likelihood of conception to first insemination of dairy cows using milk mid-infrared spectroscopy

The objective of this study was to examine the ability of milk mid-infrared (MIR) spectroscopy and other on-farm data, such as milk yield, milk composition, stage of lactation, calving age, days in milk at insemination, and somatic cell count, to identify cows that were most or least likely to conceive to first insemination. A total of 16,628 spectral and milk production records of 7,040 cows from 29 commercial dairy herds across 3 Australian states were used. Three models, comprising different explanatory variables, were tested. Model 1 included features that are readily available on farms participating in milk recording, such as milk yield, milk composition, somatic cell count, days from calving to insemination, and calving season. Days in milk and age at calving were incorporated into model 1 to form model 2. In model 3, MIR was added to model 2, but to avoid double counting, milk composition traits of model 2 were removed. The models were first trained on extreme data [i.e., including cows that (1) conceived to first insemination and (2) cows with no conception event recorded and with only 1 insemination]. Then, the models were validated in a fresh data set with all cows regardless of conception outcomes present to test for their ability to identify cows that conceived or did not conceive to first insemination. To do this, we ranked the predicted probability of all cows in the validation set and then selected the top and bottom records in varying proportions from 5 to 40% (i.e., where the model predicted the highest versus lowest likelihood of conception to first insemination, respectively) and compared with the actual values. The model's performance was evaluated through herd-year by herd-year external validation and measured as the proportion of selected records being correct. The results show that when more cows are selected (i.e., descending confidence), the accuracy of the models was reduced, and selecting the 10% of cows with the highest confidence of predictions produces optimal accuracy. Irrespective of the proportions, none of the models could predict cows that conceived to first insemination, with an accuracy around 0.48. When attempting to predict the bottom 10% of cows, which had the least likelihood of conception to first insemination, model 1 had prediction accuracy around 0.64. Compared with model 1, the addition of days in milk and calving age (model 2) resulted in a negligible improvement in prediction accuracy (0.01 to 0.03). Model 3 had the highest prediction accuracy (0.76), which implies that in the models tested, MIR is of primary importance in the prediction of fertility of dairy cows. In conclusion, this study indicates that MIR and other milk recording data could be used to identify cows with potential difficulty in getting pregnant to first insemination with promising accuracy.

Short communication: Time-dependent genetic parameters and single-step genome-wide association analyses for predicted milk fatty acid composition in Ayrshire and Jersey dairy cattle

Milk fat content and fatty acid (FA) composition have great economic value to the dairy industry as they are directly associated with taste and chemical-physical characteristics of milk and dairy products. In addition, consumers' choices are not only based on the nutritional aspects of food, but also on products known to promote better health. Milk FA composition is also related to the metabolic status and physiological stages of cows and thus can also be used as indicator for other novel traits of interest (e.g., metabolic diseases and methane yield). Genetic selection is a promising alternative to manipulate milk FA composition. In this study, we aimed to (1) estimate time-dependent genetic parameters for 5 milk FA groups (i.e., short-chain, medium-chain, long-chain, saturated, and unsaturated) predicted based on milk mid-infrared spectroscopy, for Canadian Ayrshire and Jersey breeds, and (2) conduct a time-dependent, single-step genome-wide association study to identify genomic regions, candidate genes, and metabolic pathways associated with milk FA. We analyzed 31,709 test-day records of 9,648 Ayrshire cows from 268 herds, and 34,341 records of 11,479 Jersey cows from 883 herds. The genomic database contained a total of 2,330 Ayrshire and 1,019 Jersey animals. The average daily heritability ranged from 0.18 (long-chain FA) to 0.34 (medium-chain FA) in Ayrshire, and from 0.25 (long-chain and unsaturated FA) to 0.52 (medium-chain and saturated FA) in Jersey. Important genomic regions were identified in Bos taurus autosomes BTA3, BTA5, BTA12, BTA13, BTA14, BTA16, BTA18, BTA20, and BTA21. The proportion of the variance explained by 20 adjacent SNP ranged from 0.71% (saturated FA) to 1.11% (long-chain FA) in Ayrshire, and from 0.70% (unsaturated FA) to 3.09% (medium-chain FA) in Jersey cattle. Important candidate genes and pathways were also identified, such as the PTK2 and TRAPPC9 genes, associated with milk fat percentage, and HMGCS, FGF10, and C6 genes, associated with fertility traits and immune response. Our findings on the genetic parameters and candidate genes contribute to a better understanding of the genetic architecture of milk FA composition in Ayrshire and Jersey dairy cattle.

Diagnosing the pregnancy status of dairy cows: How useful is milk mid-infrared spectroscopy?

Pregnancy diagnosis is an essential part of successful breeding programs on dairy farms. Milk composition alters with pregnancy, and this is well documented. Fourier-transform mid-infrared (MIR) spectroscopy is a rapid and cost-effective method for providing milk spectra that reflect the detailed composition of milk samples. Therefore, the aim of this study was to assess the ability of MIR spectroscopy to predict the pregnancy status of dairy cows. The MIR spectra and insemination records were available from 8,064 Holstein cows of 19 commercial dairy farms in Australia. Three strategies were studied to classify cows as open or pregnant using partial least squares discriminant analysis models with random cow-independent 10-fold cross-validation and external validation on a cow-independent test set. The first strategy considered 6,754 MIR spectra after insemination used as independent variables in the model. The results showed little ability to detect the pregnancy status as the area under the receiver operating characteristic curve was 0.63 and 0.65 for cross-validation and testing, respectively. The second strategy, involving 1,664 records, aimed to reduce noise in the MIR spectra used as predictors by subtracting a spectrum before insemination (i.e., open spectrum) from the spectrum after insemination. The accuracy was comparable with the first approach, showing no superiority of the method. Given the limited results for these models when using combined data from all stages after insemination, the third strategy explored separate models at 7 stages after insemination comprising 348 to 1,566 records each (i.e., progressively greater gestation) with single MIR spectra after insemination as predictors. The models developed using data recorded after 150 d of pregnancy showed promising prediction accuracy with the average value of area under the receiver operating characteristic curve of 0.78 and 0.76 obtained through cross-validation and testing, respectively. If this can be confirmed on a larger data set and extended to somewhat earlier stages after insemination, the model could be used as a complementary tool to detect fetal abortion.

Usefulness of milk mid-infrared spectroscopy for predicting lameness score in dairy cows

The objective of this study was to evaluate the ability of milk infrared spectra to predict cow lameness score (LMS) for use as an indicator of cow health on Australian dairy farms, or as an indicator trait for genetic evaluation purposes. The study involved 3,771 cows from 10 farms in Australia. Milk infrared spectra collected during the monthly herd testing were available in all the farms involved in the study. Lameness score was measured once in each herd, within 72 h from a test day, and merged to the closest spectra records. Lameness score was expressed on a scale from 0 to 3, where 0 is assigned to sound cows and scores 1 to 3 are assigned to cows with increased lameness severity. Partial least squares discriminant analysis was used to develop prediction models for classifying sound (score 0) and not-sound cows (i.e., cows walking unevenly, score greater than 0). Discriminant models were tested in a 10-fold random cross-validation process. Milk infrared spectra correctly classified only 57% of the cows walking unevenly and only 59% of the sound cows. When additional predictors (parity, age at calving, days in milk, and milk yield) were included in the prediction model, the model correctly classified 57% of the cows walking unevenly and 62% of the sound cows. The same model applied only to the cows in the first third of lactation correctly classified 66% of the cows walking unevenly and 57% of the sound cows. When the prediction model was used to identify lame cows (scores 2 and 3), only 49% of them were classified as such. These results are considered to be too poor to envisage a practical application of these models in the near future as on-farm tools to provide an indication of LMS. To investigate whether, at this stage, predictions of the LMS could be useful as large-scale phenotypes for animal breeding purposes, we estimated (co)variance components for actual and predicted LMS using 2,670 and 24,560 records, respectively. As the genetic correlation between actual and predicted LMS was not significantly different from zero, predictions of lameness from milk spectra and additional on-farm variables cannot be used, at this stage, as an indicator trait for actual LMS. More research is needed to find better strategies to predict lameness.

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.

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.

Invited review: Large-scale indirect measurements for enteric methane emissions in dairy cattle: A review of proxies and their potential for use in management and breeding decisions

Efforts to reduce the carbon footprint of milk production through selection and management of low-emitting cows require accurate and large-scale measurements of methane (CH4) emissions from individual cows. Several techniques have been developed to measure CH4 in a research setting but most are not suitable for large-scale recording on farm. Several groups have explored proxies (i.e., indicators or indirect traits) for CH4; ideally these should be accurate, inexpensive, and amenable to being recorded individually on a large scale. This review (1) systematically describes the biological basis of current potential CH4 proxies for dairy cattle; (2) assesses the accuracy and predictive power of single proxies and determines the added value of combining proxies; (3) provides a critical evaluation of the relative merit of the main proxies in terms of their simplicity, cost, accuracy, invasiveness, and throughput; and (4) discusses their suitability as selection traits. The proxies range from simple and low-cost measurements such as body weight and high-throughput milk mid-infrared spectroscopy (MIR) to more challenging measures such as rumen morphology, rumen metabolites, or microbiome profiling. Proxies based on rumen samples are generally poor to moderately accurate predictors of CH4, and are costly and difficult to measure routinely on-farm. Proxies related to body weight or milk yield and composition, on the other hand, are relatively simple, inexpensive, and high throughput, and are easier to implement in practice. In particular, milk MIR, along with covariates such as lactation stage, are a promising option for prediction of CH4 emission in dairy cows. No single proxy was found to accurately predict CH4, and combinations of 2 or more proxies are likely to be a better solution. Combining proxies can increase the accuracy of predictions by 15 to 35%, mainly because different proxies describe independent sources of variation in CH4 and one proxy can correct for shortcomings in the other(s). The most important applications of CH4 proxies are in dairy cattle management and breeding for lower environmental impact. When breeding for traits of lower environmental impact, single or multiple proxies can be used as indirect criteria for the breeding objective, but care should be taken to avoid unfavorable correlated responses. Finally, although combinations of proxies appear to provide the most accurate estimates of CH4, the greatest limitation today is the lack of robustness in their general applicability. Future efforts should therefore be directed toward developing combinations of proxies that are robust and applicable across diverse production systems and environments.