Genetic covariance components for measures of nitrogen utilization in grazing dairy cows

Improved nitrogen utilization of dairy production systems should improve not only the economic output of the systems but also the environmental metrics. One strategy to improve efficiency is through breeding programs. Improving a trait through breeding is conditional on the presence of exploitable genetic variability. Using a database of 1,291 deeply phenotyped grazing dairy cows, the genetic variability for 2 definitions of nitrogen utilization was studied: nitrogen use efficiency (i.e., nitrogen output in milk and meat divided by nitrogen available) and nitrogen balance (i.e., nitrogen available less nitrogen output in milk and meat). Variance components for both variables were estimated using animal repeatability linear mixed models. Genetic variability was detected for both nitrogen utilization metrics, even though their heritability estimates were low (<0.10). Validation of genetic evaluations revealed that animals divergent for nitrogen use efficiency or nitrogen balance indeed differed phenotypically, further demonstrating that breeding for improved nitrogen efficiency should result in a shift in the population mean toward better efficiency. Nitrogen use efficiency and nitrogen balance were not genetically correlated with each other (<|0.28|), and neither metric was correlated with milk urea nitrogen (<|0.12|). Nitrogen balance was unfavorably correlated with milk yield, showing the importance of including the nitrogen utilization metrics in a breeding index to improve nitrogen utilization without negatively impacting milk yield. In conclusion, improvement of nitrogen utilization through breeding is possible, even if more nitrogen utilization phenotypic data need to be collected to improve the selection accuracy considering the low heritability estimates.

Cow-level factors associated with nitrogen utilization in grazing dairy cows using a cross-sectional analysis of a large database

Reducing nitrogen pollution while maintaining milk production is a major challenge of dairy production. One of the keys to delivering on this challenge is to improve the efficiency of how dairy cows use nitrogen. Thus, estimating the nitrogen utilization of lactating grazing dairy cows and exploring the association between animal factors and productivity with nitrogen utilization are the first steps to understanding the nitrogen utilization complex in dairy cows. Nitrogen utilization metrics were derived from milk and body weight records from 1,291 grazing dairy cows of multiple breeds and crossbreeds; all cows had sporadic information on nitrogen intake concurrent with information on nitrogen sinks (and other nitrogen sources, such as body tissue mobilization). Several nitrogen utilization metrics were investigated, including nitrogen use efficiency (nitrogen output as products such as milk and meat divided by nitrogen intake) and nitrogen excreted (nitrogen intake less the nitrogen output as products such as milk and meat). In the present study, a primiparous Holstein-Friesian used, on average, 20.6% of the nitrogen it ate, excreting the surplus as feces and urine, representing 402 g of nitrogen per day. Intercow variability existed, with a between-cow standard deviation of 0.0094 for nitrogen use efficiency and 24 g of nitrogen per day for nitrogen excretion. As lactation progressed, nitrogen use efficiency declined and nitrogen excretion increased. Nevertheless, nitrogen use efficiency improved (i.e., decreased) from first to second parity, even though it did not improve from second to third parity or greater. Furthermore, nitrogen excretion continued to increase from first to third parity or greater. Nitrogen use efficiency and nitrogen excretion were negatively correlated (-0.56 to -0.40), signifying that dairy cows who partition more of the ingested nitrogen into products such as milk and meat, on average, also excrete less nitrogen. Milk urea nitrogen was, at best, weakly correlated with nitrogen use efficiency and nitrogen excretion; the correlations were between -0.01 and 0.06. In conclusion, several cow-level factors such as parity, stage of lactation, and breed were associated with the range of different nitrogen efficiency metrics investigated; moreover, even after accounting for such effects, 4.8% to 6.3% of the remaining variation in the nitrogen use efficiency and nitrogen balance metrics were attributable to intercow differences.

Usefulness of mid-infrared spectroscopy as a tool to estimate body condition score change from milk samples in intensively fed dairy cows

Directly measuring individual cow energy balance is not trivial. Other traits such as body condition score (BCS) and BCS change (ΔBCS) can, however, be used as an indicator of cow energy status. Body condition score is a metric used worldwide to estimate cow body reserves, but the estimation of ΔBCS was, until now, conditional on the availability of multiple BCS assessments. The aim of the present study was to estimate ΔBCS from milk mid-infrared (MIR) spectra and days in milk (DIM) in intensively fed dairy cows using statistical prediction methods. Daily BCS was interpolated from cubic splines fitted through the BCS records and daily ΔBCS was calculated from these splines. The ΔBCS records were merged with milk MIR spectra recorded on the same week. The dataset comprised 37,077 ΔBCS phenotypes across 9,403 lactations from 6,988 cows in 151 herds based in Quebec, Canada. Partial least squares regression (PLSR) and a neural network (NN) were then used to estimate ΔBCS from (1) MIR spectra only, (2) DIM only, or (3) MIR spectra and DIM together. The ΔBCS data in both the first 120 and 305 DIM of lactation were used to develop the estimates. Daily ΔBCS had a standard deviation of 4.40 × 10-3 BCS units in the 120-d dataset and of 3.63 × 10-3 BCS units in the 305-d dataset. A 4-fold cross-validation was used to calibrate and test the prediction equations. External validation was also conducted using more recent years of data. Irrespective of whether based on the first 120 or 305 DIM, or when MIR spectra only, DIM only or MIR spectra and DIM were jointly used as prediction variables, NN produced the lowest root mean square error (RMSE) of cross-validation (1.81 × 10-3 BCS units and 1.51 × 10-3 BCS units, respectively, using the 120-d and 305-d dataset). Relative to predictions for the entire 305 DIM, the RMSE of cross-validation was 15.4% and 1.5% lower in the first 120 DIM when using PLSR and NN, respectively. Predictions from DIM only were more accurate than those using just MIR spectra data but, irrespective of the dataset and of the prediction model used, combining DIM information with MIR spectral data as prediction variables reduced the RMSE compared with the inclusion of DIM alone, albeit the benefit was small (the RMSE from cross-validation reduced by up to 5.5% when DIM and spectral data were jointly used as model features instead of DIM only). However, when predicting extreme ΔBCS records, the MIR spectral data were more informative than DIM. Model performance when predicting ΔBCS records in future years was similar to that from cross-validation demonstrating the ability of MIR spectra of milk and DIM combined to estimate ΔBCS, particularly in early lactation. This can be used to routinely generate estimates of ΔBCS to aid in day-to-day individual cow management.

Estimation of body condition score change in dairy cows in a seasonal calving pasture-based system using routinely available milk mid-infrared spectra and machine learning techniques

Body condition score (BCS) is a subjective estimate of body reserves in cows. Body condition score and its change in early lactation have been associated with cow fertility and health. The aim of the present study was to estimate change in BCS (ΔBCS) using mid-infrared spectra of the milk, with a particular focus on estimating ΔBCS in cows losing BCS at the fastest rate (i.e., the cows most of interest to the producer). A total of 73,193 BCS records (scale 1 to 5) from 6,572 cows were recorded. Daily BCS was interpolated from cubic splines fitted through the BCS records, and subsequently used to calculate daily ΔBCS. Body condition score change records were merged with milk mid-infrared spectra recorded on the same week. Both morning (a.m.) and evening (p.m.) spectra were available. Two different statistical methods were used to estimate ΔBCS: partial least squares regression and a neural network (NN). Several combinations of variables were included as model features, such as days in milk (DIM) only, a.m. spectra only and DIM, p.m. spectra only and DIM, and a.m. and p.m. spectra as well as DIM. The data used to estimate ΔBCS were either based on the first 120 DIM or all 305 DIM. Daily ΔBCS had a standard deviation of 1.65 × 10^-3 BCS units in the 305 DIM data set and of 1.98 × 10^-3 BCS units in the 120 DIM data set. Each data set was divided into 4 sub-data sets, 3 of which were used for training the prediction model and the fourth to test it. This process was repeated until all the sub-data sets were considered as the test data set once. Using all 305 DIM, the lowest root mean square error of validation (RMSEV; 0.96 × 10^-3 BCS units) and the strongest correlation between actual and estimated ΔBCS (0.82) was achieved with NN using a.m. and p.m. spectra and DIM. Using the 120 DIM data, the lowest RMSEV (0.98 × 10^-3 BCS units) and the strongest correlation between actual and estimated ΔBCS (0.87) was achieved with NN using DIM and either a.m. spectra only or a.m. and p.m. spectra together. The RMSEV for records in the lowest 2.5% ΔBCS percentile per DIM in early lactation was reduced up to a maximum of 13% when spectra and DIM were both considered in the model compared with a model that considered just DIM. The performance of the NN using DIM and a.m. spectra only with the 120 DIM data was robust across different strata of farm, parity, year of sampling, and breed. Results from the present study demonstrate the ability of mid-infrared spectra of milk coupled with machine learning techniques to estimate ΔBCS; specifically, the inclusion of spectral data reduced the RMSEV over and above using DIM alone, particularly for cows losing BCS at the fastest rate. This approach can be used to routinely generate estimates of ΔBCS that can subsequently be used for farm decisions.

Predicting cow milk quality traits from routinely available milk spectra using statistical machine learning methods

Numerous statistical machine learning methods suitable for application to highly correlated features, as those that exist for spectral data, could potentially improve prediction performance over the commonly used partial least squares approach. Milk samples from 622 individual cows with known detailed protein composition and technological trait data accompanied by mid-infrared spectra were available to assess the predictive ability of different regression and classification algorithms. The regression-based approaches were partial least squares regression (PLSR), ridge regression (RR), least absolute shrinkage and selection operator (LASSO), elastic net, principal component regression, projection pursuit regression, spike and slab regression, random forests, boosting decision trees, neural networks (NN), and a post-hoc approach of model averaging (MA). Several classification methods (i.e., partial least squares discriminant analysis (PLSDA), random forests, boosting decision trees, and support vector machines (SVM)) were also used after stratifying the traits of interest into categories. In the regression analyses, MA was the best prediction method for 6 of the 14 traits investigated [curd firmness at 60 min, α_S1-casein (CN), α_S2-CN, κ-CN, α-lactalbumin, and β-lactoglobulin B], whereas NN and RR were the best algorithms for 3 traits each (rennet coagulation time, curd-firming time, and heat stability, and curd firmness at 30 min, β-CN, and β-lactoglobulin A, respectively), PLSR was best for pH, and LASSO was best for CN micelle size. When traits were divided into 2 classes, SVM had the greatest accuracy for the majority of the traits investigated. Although the well-established PLSR-based method performed competitively, the application of statistical machine learning methods for regression analyses reduced the root mean square error compared with PLSR from between 0.18% (κ-CN) to 3.67% (heat stability). The use of modern statistical machine learning methods for trait prediction from mid-infrared spectroscopy may improve the prediction accuracy for some traits.

Clustering high-dimensional mixed data to uncover sub-phenotypes: joint analysis of phenotypic and genotypic data

The LIPGENE-SU.VI.MAX study, like many others, recorded high-dimensional continuous phenotypic data and categorical genotypic data. LIPGENE-SU.VI.MAX focuses on the need to account for both phenotypic and genetic factors when studying the metabolic syndrome (MetS), a complex disorder that can lead to higher risk of type 2 diabetes and cardiovascular disease. Interest lies in clustering the LIPGENE-SU.VI.MAX participants into homogeneous groups or sub-phenotypes, by jointly considering their phenotypic and genotypic data, and in determining which variables are discriminatory. A novel latent variable model that elegantly accommodates high dimensional, mixed data is developed to cluster LIPGENE-SU.VI.MAX participants using a Bayesian finite mixture model. A computationally efficient variable selection algorithm is incorporated, estimation is via a Gibbs sampling algorithm and an approximate BIC-MCMC criterion is developed to select the optimal model. Two clusters or sub-phenotypes ('healthy' and 'at risk') are uncovered. A small subset of variables is deemed discriminatory, which notably includes phenotypic and genotypic variables, highlighting the need to jointly consider both factors. Further, 7 years after the LIPGENE-SU.VI.MAX data were collected, participants underwent further analysis to diagnose presence or absence of the MetS. The two uncovered sub-phenotypes strongly correspond to the 7-year follow-up disease classification, highlighting the role of phenotypic and genotypic factors in the MetS and emphasising the potential utility of the clustering approach in early screening. Additionally, the ability of the proposed approach to define the uncertainty in sub-phenotype membership at the participant level is synonymous with the concepts of precision medicine and nutrition. Copyright © 2017 John Wiley & Sons, Ltd.