Next Article in Journal
Phenotypic Characteristics Associated with Shelter Dog Adoption in the United States
Next Article in Special Issue
Predicting Daily Dry Matter Intake Using Feed Intake of First Two Hours after Feeding in Mid and Late Lactation Dairy Cows with Fed Ration Three Times Per Day
Previous Article in Journal
Preserving Porphyra umbilicalis and Saccharina latissima as Silages for Ruminant Feeding
Previous Article in Special Issue
Effect of Feeding Improved Grass Hays and Eragrostis tef Straw Silage on Milk Yield, Nitrogen Utilization, and Methane Emission of Lactating Fogera Dairy Cows in Ethiopia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Franziska Dettmann
1,2,
Daniel Warner
1,
Bart Buitenhuis
3,
Morten Kargo
3,4,
Anne Mette Hostrup Kjeldsen
4,
Niels Henning Nielsen
5,
Daniel M. Lefebvre
1 and
Debora E. Santschi
1,*
1
Lactanet, Sainte-Anne-de-Bellevue, QC H9X3R4, Canada
2
LKV Niedersachsen e.V., 26789 Leer, Germany
3
Center for Quantitative Genetics and Genomics, Aarhus University, 8830 Tjele, Denmark
4
SEGES, 8200 Aarhus N, Denmark
5
RYK, 8200 Aarhus N, Denmark
*
Author to whom correspondence should be addressed.
Submission received: 15 September 2020 / Revised: 15 October 2020 / Accepted: 20 October 2020 / Published: 23 October 2020
(This article belongs to the Collection Smart Farming in Dairy Production)

Abstract

:

Simple Summary

Dairy cows mobilize their body reserves to maintain milk production after parturition. If excessive, this may be detrimental for their health and reproduction performance. Detecting cows with an important bodyweight loss early on is essential to treat or manage these cows accordingly. Routine milk recording samples were used to develop a prediction model for early lactation body reserve mobilization in commercial dairy farms using a machine-learning approach. Bodyweight loss was mainly explained by decreased short-chain fatty acids and increased C18:0 fatty acids in milk. An early warning system for cows with an important bodyweight loss may be implemented routinely.

Abstract

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

1. Introduction

In early lactation, most dairy cows face energy deficits, due to an imbalance between feed intake and milk production requirements. Negative energy balance varies between cows in duration and extent [1], but postpartum cows typically mobilize body reserves [2,3] to maintain milk fat production [4]. Due to mobilization of body fat reserves, a certain loss in bodyweight (BW) is to be expected during early lactation, and a post-calving deterioration of 0.5 to 1.0 body condition scores (BCS) units on the 5-point scale are generally accepted in practice [5]. However, excessive mobilization of body fat might increase metabolic disorders and reproductive issues [3].
Mobilization of body fat is also reflected in a modified milk fatty acid (FA) profile. De novo FA synthesis in the mammary gland is decreased in early lactation, whereas greater proportions of preformed FA from body fat mobilization are incorporated into milk [6,7]. While BW is not typically routinely recorded in dairy farms, milk composition, and milk FA profiles can be routinely determined by Fourier-transform infrared (FTIR) spectroscopy. Therefore, this rapid analysis offers an opportunity to develop a regular, early indicator for bodyweight change (BWC) based on the FA profile of milk recording samples. The objective of this study was to evaluate whether milk FA by FTIR analysis can be used to predict BWC in early lactating cows from milk recording samples collected in commercial dairy farms.

2. Materials and Methods

Cow information, test day milk yields, and milk components, and BW records of individual cows were obtained by the Danish Cattle database (SEGES, Skejby, Denmark). All herds were equipped with an automatic milking system (AMS; Lely Nordic A/S, Fredricia, Denmark). No animal handling was required, and all producers gave their informed consent for the use of their data.

2.1. Test Day Records for Milk Components and Fatty Acids

The initial dataset included 197,058 test day records of 34,870 Danish Holstein cows from 169 herds in Denmark, collected from March 2015 to March 2017. Cows were between 5 and 305 days in milk (DIM) with parity 1 to 3. Test day records comprised 17 variables: Test day milk yield; fat; protein; somatic cell count; beta-hydroxybutyrate (BHB); the individual milk FA myristic (C14:0), palmitic (C16:0), stearic (C18:0) and oleic (C18:1) acid; groups of milk FA according to their degree of saturation (saturated FA, SFA; mono-unsaturated FA, MUFA; and poly-unsaturated FA, PUFA), and their chain length (short-chain FA, SCFA; medium-chain FA, MCFA; and long-chain FA, LCFA). Milk samples were analyzed as regular milk recording samples within a dairy herd improvement (DHI) service, using FTIR spectroscopy on Foss MilkoScan FT+ and FT 6000 (Foss Electric A/S, Hillerød, Denmark), equipped with special software (Foss Application Note 64; Foss, Hillerød, Denmark) for prediction of aforementioned FA.
Following procedures previously described [8], observations with one or more missing FA fractions and abnormal observations were removed, defined as having a PUFA concentration greater or equal to the MUFA concentration; a ratio of the sum of SFA, MUFA, and PUFA contents to total fat content of less than 0.825 or greater than 1.075 (values were chosen such that 5% of remaining observations were removed); and a ratio of the sum of SCFA, MCFA, and LCFA contents to the sum of the SFA, MUFA, and PUFA contents of less than 0.84 or greater than 1.04 (values were chosen such that 1% of remaining observations were removed).
Test day records between 5 and 35 DIM were retained for the model to predict BWC in early lactation. As test day records were generally collected monthly, for most cows (87%), only the first test day during early lactation was retained. The remainder was tested twice during early lactation.

2.2. Bodyweight Records

Bodyweight data were recorded by floor scales in Lely AMS (Lely Astronaut A3 and A4; Lely Nordic A/S, Fredricia, Denmark) at each milking. Floor scales were installed at the bottom of the AMS and were composed of four load cells in each corner. The expected tolerance for each load cell was ±10% based on the recommendations by Lely (Lely Nordic A/S, Fredericia, Denmark). The participating producers were expected to calibrate the floor scales regularly in line with recommendations by Lely. However, due to the nature of the study using large scale farm records, it was not possible to ascertain whether calibrations were done.
The initial dataset included 28,581,762 BW records of 35,787 Danish Holstein cows from 168 herds in Denmark, collected from January 2015 to September 2017. Cows were between 5 to 305 DIM with parity 1 to 3. Individual BW was corrected for potential discrepancies between AMS within a herd for herds equipped with more than one AMS, and smoothed using locally weighted smoothing procedures to correct for daily differences in BW (Figure 1), following procedures described elsewhere [9]. Bodyweight records were removed if the BW was less than 300 kg or greater than 1100 kg, and if the BW was more than 3 × the standard error of the smoothed BW. Smoothed BW records were then transformed to a relative daily BW change per DIM for each cow (expressed as g BWC/kg of BW), with the aim to generate a standardized metric across cows with different BW, as follows:
BWCDIM(x) = (BWDIM(x) − BWDIM(x − 1))/BWDIM(x − 1)
where, BWCDIM(x) is the daily BW change (g/kg of BW) at DIM x; BWDIM(x) is the smoothed daily BW (kg) at DIM x; and BWDIM(x − 1) is the smoothed BW (kg) of the previous day.

2.3. Data Editing

Records of BWC were merged with FA profiles of the respective test date. As milk fatty acid profiles were shown to change mostly within the first two days up to four days following a dietary transition [10,11], it was assumed that the FA profiles at test date were the results of the BWC observed during the previous three days. An average BWC was therefore calculated over the three days preceding the test date record (i.e., BWC from day 3 throughout test day) and merged with the respective test day record. As the average daily BCW for DIM 5 and 6 could not be computed, the earliest test day records retained were those for DIM 7. The final data set included 19,138 test day records between DIM 7 and 35 with the respective average BWC. In total, 16,847 early lactation Danish Holstein cows across 165 herds were retained for analysis (Table 1). Cows had an average milk yield of 35.9 kg at test day (standard deviation (SD) of 10.90; range from 2.1 to 76.8 kg).

2.4. Prediction Model

Feature selection was conducted through a principal component analysis with the aim of keeping only the most important information in the dataset. All variables included in the first three principal components (explaining 80% of the total variance) after plotting eigenvalues according to their size were retained [12]. Somatic cell count was omitted from the final dataset, due to its low contribution (<5%) and low correlation (r < 0.50) with any of the three principal components. The final dataset comprised 16 variables, namely, the milk components and milk FA presented in Table 1, milk yield, parity (1, 2, and 3), DIM, test date year, and month. None of the predictors had a near-zero variance. Variables were normalized to a mean of 0 and a standard deviation of 1.
Prediction of BWC at test day was based on decision tree induction using a random forest algorithm [13]. A random forest is an ensemble of multiple decision trees, and thus, essentially the sum of piecewise linear functions. The dataset is split into smaller sets that become more manageable by repeatedly selecting a random sample with replacement from the training set and fitting multiple decision trees to these samples (i.e., each decision tree is shown a different training set). In addition, predictors are randomly selected to prevent trees from becoming correlated. Thus, a random forest model may deal with a multitude of linear and nonlinear relationships among predictors, and the complexity inherent to a high-dimensional dataset [14].

2.5. Model Development

The model was implemented in R (version 3.5.0; R Foundation for Statistical Computing, Vienna, Austria) using the modeling package workflow caret [15]. The dataset was randomly split into a training set (80%) used for model preparation, and an evaluation set (20%) containing output values withheld from the algorithm with the aim to evaluate predictions from the trained model with previously unseen BWC records. Data splitting was done based on two different stratification approaches to ensure that data of a specific cow (grouped by animal) or specific herd (grouped by herd) were used exclusively either to train or to validate the model. The aim of the latter approach was to ensure that herd-level information was not inadvertently passed to both the training and evaluation set when cows from the same herd were assigned.
The model was prepared using 10-fold cross-validation with three iterations. For each iteration, a model was trained on nine splits of the data set and cross-validated on the remaining part of the data set (i.e., one split). Like the data splitting approach, described above, either cows or herds were randomly assigned to the cross-validation folds—such that data of a specific cow or herd are used exclusively either to train or to cross-validate the model. The optimal hyper-parameter configuration for the random forest algorithm was evaluated based on the repeated cross-validation through a random parameter search [16] of the number of classification trees (500, 1000, 1500, and 2000 trees) and the number of randomly selected predictors (1 to 10). The final model configuration was 500 trees with three randomly selected predictors. Model accuracy was estimated on the average of the 10-fold repeated cross-validation and compared to the validation set based on Pearson’s r correlation coefficient, Lin’s concordance correlation coefficient (CCC), the coefficient of determination (R2), the root mean squared error of prediction (RMSEP), the mean absolute error (MAE), the ratio of performance to deviation (RPD; defined as the SD of observed values divided by the RMSEP [17]), and the ratio of performance to the interquartile range (RPIQ; defined as the interquartile range of the observed values divided by the RMSEP [18]). The predictive performance of the model is best with small values for RMSEP and MAE, and large values for r, CCC, R2, RDP, and RPIQ. The variables contributing most to the model development (i.e., the global variable importance scores) were extracted, and overall variable importance was determined for each tree based on the mean decrease in node impurity as determined by the Gini index and averaged across trees. The model response was visualized with an alluvial diagram for the variables with the highest feature importance with variables normalized to a mean of 0 and a standard deviation of 1. In addition, localized variable importance scores were extracted for randomly selected cows to interpret the model response for individual cows by using the break-down approach implemented in the DALEX package in R 3.5.0 [19]. This procedure sequentially attributes variables to the model output, thus, decomposing predictions and interpreting the link between input variables and model output [20].

3. Results and Discussion

3.1. Milk Fatty Acid Profiles

Milk FA profiles differed with parity and lactation stage, as shown in Figure 2. In early lactation (defined as 7 to 35 DIM in the present study), the concentrations of SCFA and MCFA increased, whereas the initially high levels of LCFA decreased rapidly as DIM progressed. These changes in the milk FA profiles can be explained by increased body fat mobilization [2,3] to maintain milk production [4]. With increasing parity, the concentration of FA increased, due to a greater milk fat production and a stronger body fat mobilization, due to a stronger negative energy balance in cows at later parity [21]. Our findings are in line with those reported in the literature. Molar proportions of the SCFA C10:0 to C14:0 were significantly lower in early lactation than in mid-lactation [7]. De novo synthesis of SCFA in the mammary gland might be thereby inhibited by LCFA from body fat [6] and linked to the adaptation of the rumen microflora to the lactation diet.
In adipose tissue, the LCFA C18:1c9 and C18:0 and the MCFA C16:0 account for nearly 90% of the FA at approximately equal molar proportions [22]. With the mobilization of body reserves during early lactation, these FA would be more likely incorporated into milk fat. Nonetheless, while daily yields of C18:0 and C18:1 were indeed greater in early lactation, daily yields of C16:0 were smaller in early lactation in the present study. These results confirmed previous findings [23,24].

3.2. Relative BW Change

Bodyweight curves were distinctive across parities, with a lower average BW for cows in first parity (603 ± 66.1 kg; mean ± SD) than cows in second (656 ± 67.7 kg), or third parity (689 ± 70.2 kg). The relative daily BWC in early lactation was −0.52 ± 2.65 g/kg of BW, −0.64 ± 2.82 g/kg of BW, and −0.82 ± 5.53 g/kg of BW for cows in first, second and third parity, respectively (Figure 3). These results are comparable to the range of BW losses reported for different parities [21] with average BW losses of 6.5%, 8.5%, and 8.4% from calving to nadir BW for cows in first, second, and third or greater parity, respectively. Overall, these results suggest that the level and duration of negative energy balance increased with later parity in line with findings in the literature [21]. However, these may not be valid for all cows as cows with a high BW loss in the previous lactation were shown to have increased BW losses at later parity (5.87%, 6.49%, and 7.89% BW loss for cows in first, second and third parity, respectively, from calving until five weeks of lactation), but cows with a low BW loss in the previous lactation had decreased BW losses at later parity (7.27%, 4.83%, and 5.45% BW loss) [25]. Variation in BWC may be also be explained by differences in gastrointestinal size and content [26]. In the present study, BWC reflects variation in BW as mobilization or deposition and changes in gastrointestinal fill [3].
In practice, the extent of BW loss might be more easily interpreted through a deterioration of the body condition. While a BCS loss of no more than 0.5 to 1.0 units on a 5-point scale should be targeted in practice [5], an average daily BW loss of 3 g/kg of BW over a 28-day period from DIM 7 to 35 can be translated into a theoretical cumulative loss of 1.2 to 1.5 BCS on the 5-point scale whilst considering the 10% range of error of the floor scales (an average increase of 41 kg BW per unit of increase in BCS for Holstein and mixed herds was assumed [27]).

3.3. Model Performance and Evaluation

Daily BWC during early lactation was predicted from test-day milk components and FA using decision-tree induction through a random forest algorithm. Decision-tree induction is a flexible approach using all available variables. However, the variables that most influenced the predictions with an overall importance of over 70% across all decision trees are SCFA, C18:0, and MUFA (Figure 4). Due to the presence of collinearity among the individual FA, care should be taken as to the ranking of the individual variable as their importance might not be well defined. Milk BHB was not considered an important feature in the present study despite its known association with hyperketonemia, a metabolic disorder commonly affecting cows in negative energy balance after parturition [28,29]. However, the prevalence of 8.4% for elevated milk BHB (>0.15 mmol/L) in the present study was considerably lower than that observed in dairy cows in Canada (21% and 22.6%) [28,29], and across Europe (21.8%) [30]. In addition to possible genetic differences, the low BHB prevalence might explain its low overall importance for BWC in our model.
Using a model-agnostic approach, the model predictions can be investigated in more detail (Figure 5). For cows with a large BW loss in early lactation, the predicted BWC were explained by generally lower concentrations of SCFA and higher concentrations of C18:0 and MUFA. In contrast, for cows with a BW gain in early lactation, predicted BWC was explained by generally higher concentrations of SCFA and lower concentrations of C18:0 and MUFA. Decision rules may thereby vary for each prediction instance using decision-tree induction. For instance, the predicted BWC was influenced by different predictors for two randomly selected cows from the same herd with a relatively large BW loss in early lactation (Figure 6). For a cow at DIM 25 with a predicted BWC of −3.69 g/kg of BW (observed BWC of −4.05), BWC was mainly positively influenced by LCFA and negatively influenced by SCFA and C18:0 (Figure 6A). In contrast, for a cow at DIM 10 with a predicted BWC of −1.98 g/kg of BW (observed BWC of −2.22), BWC was mainly negatively influenced by DIM, C14:0, but positively influenced by LCFA and milk BHB (Figure 6B). This example further shows that decision rules can be interpreted for single prediction instances despite the complexity of a random forest algorithm relying on the prediction of a multitude of decision trees. By integrating additional information, such as diet information, to test day milk fatty acids, it might be possible to investigate the specific pathways shown below in more detail.
Predicted BWC after cross-validation (median ± absolute deviation; −0.44 ± 2.25 g/kg of BW per day) did not significantly differ (p > 0.05) from observed BWC (−0.41 ± 1.64 g/kg of BW per day), suggesting that a random forest algorithm was capable of predicting BWC. Nonetheless, predictions of BWC were off the range of the observed BWC (interquartile range of 3.04 for predicted values as compared to 2.24 g/kg of BW for observed values). The prediction model tended to underestimate BWC for cows with a large BW loss (Figure 7). However, prediction of an absolute value of BWC might be less relevant for dairy producers if cows with a high BW loss can be correctly classified as such. Therefore, future applications should focus on classifying cows with a high versus a low BWC rather than predicting an absolute value for BWC.
Evaluation with previously unseen BWC records diminished the prediction performance model (Table 2). It is assumed that the model performance after cross-validation might be over-optimistic, due to the expected collinearity among FA. This limits the applicability of this model in the current state. Moreover, although data were randomly assigned to either the training or the evaluation set, the observed range in the training set was somewhat smaller than that in the evaluation set (see Table A1 in the Appendix).
Other limitations to the model likely originate from underlying milk FA and BW records used for the prediction model. The milk FA was estimated from FTIR measurements with partial least squares algorithms that have a different degree of accuracy [31]. The coefficient of variation ranged from 3.6 to 12.8% for the main milk FA, with SFA and SCFA algorithms performing best with a coefficient of variation below 5%. The algorithms for the minor milk FA were weaker, with a coefficient of variation of 14.4% for PUFA and 27.2% for TFA. Milk FA measured by a reference method, such as gas chromatography, is expected to be more reliable for BWC predictions, but is not routinely used for milk recording samples, due to the high analytical costs and low analytical throughput.
The prediction performance may also suffer from the daily fluctuations of BW measurements from floor scales, as previously observed [26], in particular in commercial dairy farms. These fluctuations commonly occur, due to daily variations in milk yield, feed intake, water intake, and time of milking in relation to the time of BW measurements, and might be thereby confounded with changes of physiological significance. In addition, due to the nature of the study using a large set of farm records, it was not possible to ascertain whether the floor scales were regularly calibrated.
The overall prediction performance might be further improved by a more frequent milk sampling scheme than the conventional once-a-month sampling by feeding the model with more relevant cow-side information during early lactation. This would better intercept the dynamics of milk composition change, due to changes in the metabolic status of the cow occurring after the onset of lactation.
Further improvements to the prediction performance can be expected by including additional biomarkers in milk, such as acetone and citrate, that can be measured through FTIR and were previously associated with ketosis and negative energy balance in early lactation cows, respectively [32]. Future work should likewise focus on the use of the FTIR milk spectrum that may include additional model features to the milk composition and milk FA profile currently available. Despite the lack of a direct biochemical association between the FTIR spectrum and bodyweight change, as opposed to milk FA that peak in various regions of the FTIR spectrum, the potential of using FTIR spectra for physiological processes was previously shown, among others, for predicting the energy status of cows in early lactation [33,34].

4. Conclusions

The results of this study suggest that BWC might be estimated by FTIR milk FA profiles in milk recording samples. Nonetheless, before this prediction model can be used in commercial dairy farms, the model needs to be validated for the different country- or region-specific conditions, breeds, and herd management and feeding strategies. Further work analyses are needed to assess the impact of the level of BWC on milk production, reproductive performance, and health. For use in a real-time decision-support tool, predictions of BWC using early lactation test day may not be practical as any suggested modification in management strategies will likely not improve reproductive performance in the current lactation. However, the concept described here might help identify deficiencies on the herd level and improve overall herd management as a preventive strategy to improve herd performance in future lactations.

Author Contributions

Conceptualization, D.W., B.B., M.K., A.M.H.K., N.H.N., D.M.L. and D.E.S.; Data curation, D.W.; Formal analysis, F.D. and D.W.; Funding acquisition, D.M.L. and D.E.S.; Investigation, F.D. and D.W.; Methodology, F.D., D.W. and D.E.S.; Project administration, D.E.S.; Resources, B.B., M.K., A.M.H.K., N.H.N., D.M.L. and D.E.S.; Software, D.W.; Supervision, D.E.S.; Validation, D.W., B.B., M.K., A.M.H.K., N.H.N. and D.E.S.; Visualization, D.W.; Writing—Original draft, F.D. and D.W.; Writing—Eview & editing, D.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We thank Lely Nordic A/S (Fredericia, Denmark) for supplying the bodyweight records, and H. Martinussen (SEGES, Denmark), I. B. Skaarup Hansen (RYK, Denmark), and R. Lacroix (Lactanet, QC, Canada) for support with data analysis. Parts of this study have been presented at the ICAR (International Committee for Animal Recording) conference in Prague, Czech Republic, 17–21 June 2019, and at the ADSA Annual Meeting in Cincinnati, OH, 23–26 June 2019.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Distribution of predictors among the training and evaluation set with herds randomly assigned to either of the two datasets.
Table A1. Distribution of predictors among the training and evaluation set with herds randomly assigned to either of the two datasets.
VariableTraining SetSDP1P99CountEvaluation SetSDP1P99Count
MeanMean
Parity 1 42.8% 40.9%
Parity 2 34.9% 35.4%
Parity 3 22.3% 23.7%
DIM21.78.10735 21.78.19735
Milk yield (kg)36.110.9514.858.9 35.510.7113.260.6
Fat (%)4.370.9592.496.97 4.380.9672.467.26
Protein (%)3.380.3242.714.21 3.380.3202.754.26
SCC22570173878 24075563359
BHB0.060.0610.000.32 0.060.0720.000.27
BWC (g/kg BW)−0.623.714−8.995.70 −0.662.889−8.795.47
Fatty acids (g/100 g milk)
SFA2.620.5831.434.09 2.620.5761.404.29
MUFA1.280.3990.652.54 1.290.4150.632.58
PUFA0.170.0460.070.28 0.170.0460.070.30
SCFA0.430.1040.220.69 0.430.1010.220.73
LCFA1.830.5760.853.65 1.850.6000.803.64
MCFA1.550.3810.792.52 1.550.3780.792.62
C16:01.080.2560.601.76 1.080.2530.571.83
C14:00.360.0900.190.59 0.360.0870.190.62
C18:00.560.1620.281.03 0.570.1660.261.06
C18:11.150.3820.542.38 1.160.3980.512.37
Abbreviations used: SD = standard deviation, P1 = 1st percentile, P99 = 99th percentile, SFA = saturated fatty acids, MUFA = mono-unsaturated fatty acid, PUFA = poly-unsaturated fatty acids, SCFA = short-chain fatty acid, MCFA = medium-chain fatty acid, LCFA = long-chain fatty acid, SCC = somatic cell count, BHB = beta-hydroxybutyrate.

References

  1. Jorritsma, R.; Wensing, T.; Kruip, T.A.M.; Vos, P.L.A.M.; Noordhuizen, J.P.T.M. Metabolic changes in early lactation and impaired reproductive performance in dairy cows. Vet. Res. 2003, 34, 11–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Andrew, S.M.; Waldo, D.R.; Erdman, R.A. Direct analysis of body composition of dairy cows at three physiological stages. J. Dairy Sci. 1994, 77, 3022–3033. [Google Scholar] [CrossRef]
  3. Grummer, R.R.; Rastani, R.R. Review: When should lactating dairy cows reach positive energy balance? Prof. Anim. Sci. 2003, 19, 197–203. [Google Scholar] [CrossRef]
  4. Bar-Peled, U.; Lehrer, A.R.; Folman, Y.; Bruckental, I.; Kali, J.; Gacitua, H.; Maltz, E.; Tagari, H.; Robinzon, B. Milk production, and some related nutritional and reproductive variables, of dairy cows under two different milking regimes in early lactation. In Prospects for Automatic Milking; EAAP Publ. No. 65; Ipema, A.H., Lipus, A.C., Metz, J.H.M., Rossing, W., Eds.; PU-DOC Scientific Publishers: Wageningen, The Netherlands, 1992; pp. 219–226. [Google Scholar]
  5. Roche, J.R.; Friggens, N.C.; Kay, J.K.; Fisher, M.W.; Stafford, K.J.; Berry, D.P. Invited review: Body condition score and its association with dairy cow productivity, health, and welfare. J. Dairy Sci. 2009, 92, 5769–5801. [Google Scholar] [CrossRef] [Green Version]
  6. Palmquist, D.L.; Beaulieu, D.; Barbano, D.M. Feed and animal factors influencing milk fat composition. J. Dairy Sci. 1993, 76, 753–1771. [Google Scholar] [CrossRef]
  7. Garnsworthy, P.C.; Masson, L.L.; Lock, A.L.; Mottram, T.T. Variation of milk citrate with stage of lactation and de novo fatty acid synthesis in dairy cows. J. Dairy Sci. 2006, 89, 1604–1612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Hein, L.; Sørensen, L.P.; Kargo, M.; Buitenhuis, A.J. Genetic analysis of predicted fatty acid profiles of milk from Danish Holstein and Danish Jersey cattle populations. J. Dairy Sci. 2018, 101, 2148–2157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Bossen, D.; Hlidasdottir, T.; Thorup, V.M. Diurnal variation in live weight for evaluation of feed ration allowance and intake. In Proceedings of the Precision Dairy Conference and Expo, Rochester, MA, USA, 26–27 June 2013; pp. 179–180. [Google Scholar]
  10. Elgersma, A.; Ellen, G.; van der Horst, H.; Boer, H.; Dekker, P.R.; Tamminga, S. Quick changes in milk fat composition from cows after transition from fresh grass to a silage diet. Anim. Feed Sci. Technol. 2004, 117, 13–27. [Google Scholar] [CrossRef]
  11. Gervais, R.; Université Laval, Québec, QC, Canada. Personal communication, 2020.
  12. Peres-Neto, P.R.; Jackson, D.A.; Somers, K.M. How many principal components? Stopping rules for determining the number of non-trivial axes revisited. Br. J. Stat. Psych. 2005, 49, 974–997. [Google Scholar] [CrossRef]
  13. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar]
  14. Breiman, L. Statistical modeling: The two cultures. Stat. Sci. 2001, 16, 199–231. [Google Scholar]
  15. Kuhn, M. Building predictive models in R using the caret package. J. Stat. Softw. 2008, 28, 1–26. [Google Scholar]
  16. Bergstra, J.; Bengio, Y. Random search for hyper-parameter optimization. J. Mach. Learn. Res. 2012, 13, 281–305. [Google Scholar]
  17. Williams, P.C.; Sobering, D.C. Comparison of commercial near infrared transmittance and reflectance instruments for analysis of whole grains and seeds. J. Near Infrared Spectrosc. 1993, 1, 25–32. [Google Scholar]
  18. Bellon-Maurel, V.; Fernandez-Ahumada, E.; Palagos, B.; Roger, J.-M.; McBratney, A. Critical review of chemometric indicators commonly used for assessing the quality of the prediction of soil attributes by NIR spectroscopy. Trends Anal. Chem. 2010, 29, 1073–1081. [Google Scholar] [CrossRef]
  19. Biecek, P. DALEX: Explainers for Complex Predictive Models in R. J. Mach. Learn. Res. 2018, 19, 1–5. [Google Scholar]
  20. Gosiewska, A.; Biecek, P. Do not trust additive explanations. arvXi 2019, arXiv:1903.11420. [cs.LG]. Available online: https://arxiv.org/abs/1903.11420 (accessed on 22 October 2020).
  21. Van Straten, M.; Shpigel, N.Y.; Friger, M. Analysis of daily body weight of high-producing dairy cows in the first one hundred twenty days of lactation and associations with ovarian inactivity. J. Dairy Sci. 2008, 91, 3353–3362. [Google Scholar] [CrossRef] [Green Version]
  22. Christie, W.W. The composition, structure and function of lipids in the tissues of ruminant animals. In Lipid Metabolism in Ruminant Animals; Christie, W.W., Ed.; Pergamon Press: Oxford, UK, 1981; pp. 95–191. [Google Scholar]
  23. Artegoitia, V.; Meikle, A.; Olazabal, L.; Damián, J.P.; Adrien, M.L.; Mattiauda, D.A.; Bermudez, J.; Torre, A.; Carriquiry, M. Milk casein and fatty acid fractions in early lactation are affected by nutritional regulation of body condition score at the beginning of the transition period in primiparous and multiparous cows under grazing conditions. J. Anim. Physiol. Anim. (Berl.) 2014, 97, 919–932. [Google Scholar] [CrossRef]
  24. Mäntysaari, P.; Mäntysaari, E.A.; Kokkonen, T.; Mehtiö, T.; Kajava, S.; Grelet, C.; Lidauer, P.; Lidauer, M.H. Body and milk traits as indicators of dairy cow energy status in early lactation. J. Dairy Sci. 2019, 102, 7904–7916. [Google Scholar] [CrossRef] [Green Version]
  25. Zachut, M.; Moallem, U. Consistent magnitude of postpartum body weight loss within cows across lactations and the relation to reproductive performance. J. Dairy Sci. 2017, 100, 3143–3154. [Google Scholar] [CrossRef] [PubMed]
  26. Maltz, E.; Devir, S.; Metz, J.H.M.; Hegeveen, H. The body weight of the dairy cow: I. Introductory study into body weight changes in dairy cows as a management aid. Livest. Prod. Sci. 1997, 48, 175–186. [Google Scholar] [CrossRef]
  27. Enevoldsen, C.; Kristensen, T. Estimation of body weight from body size measurements and body conditions scores in dairy cows. J. Dairy Sci. 1997, 80, 1988–1995. [Google Scholar] [CrossRef]
  28. Santschi, D.E.; Lacroix, R.; Durocher, J.; Duplessis, M.; Moore, R.K.; Lefebvre, D.M. Prevalence of elevated milk β-hydroxybutyrate concentrations in Holstein cows measured by Fourier-transform infrared analysis in Dairy Herd Improvement milk samples and association with milk yield and components. J. Dairy Sci. 2016, 99, 9263–9270. [Google Scholar] [CrossRef]
  29. Tatone, E.H.; Duffield, T.F.; LeBlanc, S.J.; DeVries, T.J.; Gordon, J.L. Investigating the within-herd prevalence and risk factors for ketosis in dairy cattle in Ontario as diagnosed by the test-day concentration of β-hydroxybutyrate in milk. J. Dairy Sci. 2017, 100, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Suthar, V.S.; Canelas-Raposo, J.; Deniz, A.; Heuwieser, W. Prevalence of subclinical ketosis and relationships with postpartum diseases in European dairy cows. J. Dairy Sci. 2013, 96, 2925–2938. [Google Scholar] [CrossRef] [Green Version]
  31. FOSS. Application Note 64; FOSS Analytical A/S: Hillerød, Denmark, 2009; pp. 11–12. [Google Scholar]
  32. Grelet, C.; Bastin, C.; Gelé, M.; Davière, J.-B.; Johan, M.; Werner, A.; Reding, R.; Fernandez Pierna, J.A.; Colinet, F.G.; Dardenne, P.; et al. Development of Fourier transform mid-infrared calibrations to predict acetone, β-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network. J. Dairy Sci. 2016, 99, 4816–4825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Smith, S.L.; Denholm, S.J.; Coffey, M.P.; Wall, E. Energy profiling of dairy cows from routine milk mid-infrared analysis. J. Dairy Sci. 2019, 102, 11169–11179. [Google Scholar] [CrossRef]
  34. McParland, S.; Banos, G.; Wall, E.; Coffey, M.P.; Soyeurt, H.; Veerkamp, R.F.; Berry, D.P. The use of mid-infrared spectrometry to predict body energy status of Holstein cows. J. Dairy Sci. 2011, 94, 3651–3661. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Daily smoothed bodyweight of dairy cows across 165 Danish Holstein herds from 7 to 305 days in milk.
Figure 1. Daily smoothed bodyweight of dairy cows across 165 Danish Holstein herds from 7 to 305 days in milk.
Animals 10 01958 g001
Figure 2. Distribution of short-chain (SCFA), medium-chain (MCFA), and long-chain (LCFA) fatty acids in g/day of Danish Holstein cows across 165 herds. The shaded vertical bars represent early lactation from DIM 7 to 35.
Figure 2. Distribution of short-chain (SCFA), medium-chain (MCFA), and long-chain (LCFA) fatty acids in g/day of Danish Holstein cows across 165 herds. The shaded vertical bars represent early lactation from DIM 7 to 35.
Animals 10 01958 g002
Figure 3. Mean daily bodyweight (BW) change of dairy cows across 165 Danish Holstein herds.
Figure 3. Mean daily bodyweight (BW) change of dairy cows across 165 Danish Holstein herds.
Animals 10 01958 g003
Figure 4. Relative overall importance after cross-validation of the top 10 predictors with the highest feature importance used by the random forest algorithm.
Figure 4. Relative overall importance after cross-validation of the top 10 predictors with the highest feature importance used by the random forest algorithm.
Animals 10 01958 g004
Figure 5. Model response, as shown by an alluvial diagram with predicted bodyweight change (BWC), ranked from negative to positive BWC as affected by the predictors with the highest feature importance. Nodes on the left denote categories of predicted BWC with the respective range for groups of cows in the bottom 10% (bottom node), intermediate (10–90%), and upper 10% (upper node) predicted BWC; nodes for predictors with respective means were fixed to five nodes based on the frequency distribution.
Figure 5. Model response, as shown by an alluvial diagram with predicted bodyweight change (BWC), ranked from negative to positive BWC as affected by the predictors with the highest feature importance. Nodes on the left denote categories of predicted BWC with the respective range for groups of cows in the bottom 10% (bottom node), intermediate (10–90%), and upper 10% (upper node) predicted BWC; nodes for predictors with respective means were fixed to five nodes based on the frequency distribution.
Animals 10 01958 g005
Figure 6. Contribution of variables for predicted bodyweight change (g/kg of BW per day; x-axis) of two randomly selected cows within the same herd experiencing a large BW loss during early lactation. Cow 36362 (DIM 25): −4.05 g/kg of BW per day (observed) versus −3.69 g/kg of BW per day (predicted) (A); cow 36357 (DIM 10): −2.22 g kg of BW per day (observed) versus −1.98 g/kg of BW per day (predicted) (B).
Figure 6. Contribution of variables for predicted bodyweight change (g/kg of BW per day; x-axis) of two randomly selected cows within the same herd experiencing a large BW loss during early lactation. Cow 36362 (DIM 25): −4.05 g/kg of BW per day (observed) versus −3.69 g/kg of BW per day (predicted) (A); cow 36357 (DIM 10): −2.22 g kg of BW per day (observed) versus −1.98 g/kg of BW per day (predicted) (B).
Animals 10 01958 g006
Figure 7. Predicted versus observed bodyweight change (BWC) with Pearson’s r and Lin’s concordance correlation coefficient (CCC) after cross-validation.
Figure 7. Predicted versus observed bodyweight change (BWC) with Pearson’s r and Lin’s concordance correlation coefficient (CCC) after cross-validation.
Animals 10 01958 g007
Table 1. Test day records and bodyweights (BW) of 16,847 Danish Holstein cows in early lactation (7‒35 DIM) across 165 herds.
Table 1. Test day records and bodyweights (BW) of 16,847 Danish Holstein cows in early lactation (7‒35 DIM) across 165 herds.
TraitParity 1 (n = 8113) 1Parity 2 (n = 6698) 1Parity 3 (n = 4327) 1
MeanSDP1P99MeanSDP1P99MeanSDP1P99
Milk yield 1 (kg)27.46.511.642.941.18.817.160.643.39.616.763.3
Fat 1 (%)4.500.972.567.444.250.912.436.864.330.982.417.25
Protein 1 (%)3.410.312.784.223.370.332.724.283.340.342.714.26
Ratio fat-to-protein 11.320.280.772.191.260.270.732.101.300.290.752.28
BHB 1 (mmol/L)0.0480.0510.0000.2300.0670.0660.0000.2930.0770.0750.0000.337
Fatty acid 1 (g/100 g of milk)
C14:00.370.090.200.630.360.090.190.590.360.090.180.61
C16:01.110.250.621.871.050.250.561.741.050.260.561.81
C18:00.590.170.291.090.540.150.250.980.560.170.261.08
C18:11.200.400.542.461.110.350.512.221.150.400.512.48
SFA2.680.581.484.392.570.571.334.032.600.601.404.32
MUFA1.330.410.662.671.230.370.622.421.280.420.622.65
PUFA0.170.050.080.310.160.040.010.280.160.040.070.30
SCFA0.430.110.230.740.430.100.220.690.440.100.230.72
MCFA1.590.370.842.681.520.380.752.511.520.390.762.64
LCFA1.900.590.853.741.760.540.783.421.830.600.793.77
Initial BW 2 (kg)57260.343271764164.849180468170.7523854
Average BW 3 (kg)60366.145877265667.750583468970.2531870
BWC 4 (g/kg of BW per day)−0.522.65−7.775.17−0.642.82−9.065.59−0.825.53−10.508.57
Abbreviations used: SD = standard deviation, P1 = 1st percentile, P99 = 99th percentile. 1 Test-day records (abbreviations used: BHB = beta-hydroxybutyrate, SFA = saturated fatty acids, MUFA = mono-unsaturated fatty acid, PUFA = poly-unsaturated fatty acids, SCFA = short-chain fatty acid, MCFA = medium-chain fatty acid, LCFA = long-chain fatty acid); 2 At day in milk 7; 3 Between day in milk 7 and 35; 4 Bodyweight change.
Table 2. Performance of a random forest algorithm to predict bodyweight change.
Table 2. Performance of a random forest algorithm to predict bodyweight change.
Model Metric 1Grouped by Animal 2EvaluationGrouped by Herd 3Evaluation
Cross-ValidationCross-Validation
Pearson’s r0.950.490.940.55
Lin’s CCC0.850.390.830.46
R20.910.240.890.31
RMSEP1.692.571.792.33
MAE0.761.710.731.71
RPD2.181.152.081.20
RPIQ1.811.191.701.35
1 CCC = concordance correlation coefficient; RMSEP = root mean square error of the prediction; MAE = mean absolute error; RPD = ratio of performance to deviation; RPIQ = ratio of performance to the interquartile range. 2 Data splitting for model training (80%) and model validation (20%) stratified by cow. 3 Data splitting for model training (80%) and model validation (20%) stratified by the herd.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dettmann, F.; Warner, D.; Buitenhuis, B.; Kargo, M.; Kjeldsen, A.M.H.; Nielsen, N.H.; Lefebvre, D.M.; Santschi, D.E. Fatty Acid Profiles from Routine Milk Recording as a Decision Tool for Body Weight Change of Dairy Cows after Calving. Animals 2020, 10, 1958. https://0-doi-org.brum.beds.ac.uk/10.3390/ani10111958

AMA Style

Dettmann F, Warner D, Buitenhuis B, Kargo M, Kjeldsen AMH, Nielsen NH, Lefebvre DM, Santschi DE. Fatty Acid Profiles from Routine Milk Recording as a Decision Tool for Body Weight Change of Dairy Cows after Calving. Animals. 2020; 10(11):1958. https://0-doi-org.brum.beds.ac.uk/10.3390/ani10111958

Chicago/Turabian Style

Dettmann, Franziska, Daniel Warner, Bart Buitenhuis, Morten Kargo, Anne Mette Hostrup Kjeldsen, Niels Henning Nielsen, Daniel M. Lefebvre, and Debora E. Santschi. 2020. "Fatty Acid Profiles from Routine Milk Recording as a Decision Tool for Body Weight Change of Dairy Cows after Calving" Animals 10, no. 11: 1958. https://0-doi-org.brum.beds.ac.uk/10.3390/ani10111958

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop