If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
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.
The process of lactogenesis is energetically demanding for dairy cows. If the cow does not consume enough energy to cover her energy requirements, she has to mobilize her own energy reserves to meet the deficit and can enter a state called negative energy balance (NEB;
). This is because, in addition to the energy requirements for lactation, maintenance, and mobility, energy intake is also a metabolic and hormonal function that underlies reproductive and immune performance (
demonstrated a significant negative genetic correlation between days to first observed estrous and energy balance (−0.16 to −0.28). Therefore, the use of particular cow energy traits for inclusion in selection criteria to improve resumption of cycling in subsequent lactations could be advocated. However, calculating energy traits frequently enough to monitor NEB levels and loss of condition, at both the herd and the individual cow level across time, for cows in commercial milking herds has been an inherently difficult thing to do in practice. This is due to the fact that recording feed intake for individual cows is impractical in terms of both time and expense (
). Using the traditional approach of subjective BCS is reactive in the sense that, by the time the change in phenotype can be recorded, the damage may already have been done. Prediction equations for energy balance have previously been developed based directly on mid-infrared (MIR) spectral data, a by-product generated from routine milk recording, and using these experimental data to validate calculations (
). Monitoring NEB and feed intake on a regular basis using MIR analyses of the milk would have distinct advantages: (1) spectral data are already gathered for predicting milk components, so additional predicted traits have only a small marginal cost; (2) any problems can be identified quickly and mitigation strategies deployed immediately.
Mid-infrared spectroscopy uses infrared light from the electromagnetic spectrum and works because the particular frequency absorbed by a molecule is determined by its shape and chemical bonds. The technique is used routinely for assessing the fat and protein content of milk as well as lactose, urea, and casein content (
). It is becoming increasingly apparent that different physiological processes can leave molecular signatures in the milk of a dairy cow that could potentially be detected by this technique.
used a reference set of bovine milk fatty acids (FA) and their milk MIR spectra to generate multivariate prediction equations for over 30 FA, which have been routinely updated as part of an EU-funded project, Optimir (www.optimir.eu), with accuracies ranging from 68% to 100% (
Creation of universal MIR calibrations by standardization of milk spectra: Example of fatty acids.
in: Book of abstracts of the 65th annual meeting of the European Federation of Animal Science. Wageningen Academic Publishers,
Copenhagen, Denmark2014: 108
Development of Fourier transform mid-infrared calibrations to predict acetone, β-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network.
The objective of this study was to use reference phenotypic data obtained from the dairy research herd at Scotland's Rural College (SRUC) and combine it with concurrent MIR spectral data to generate prediction equations for energy balance (EB) and energy intake (EI) such that predictions of these traits can be generated for cows in the national herd as part of routine milk recording and genetic analyses performed on MIR-derived predictions.
MATERIALS AND METHODS
Sample Data and Generation of Energy Phenotypes
Phenotypic data were collected from the Langhill lines of Holstein Friesian (HOL) dairy cows (n = 1,101) housed at SRUC Dairy Research Centre at Crichton Royal Farm in Dumfries, Scotland, and recorded between 2005 and 2018.
Dairy cattle in this herd have been subject to a long-term 2 × 2 factorial experiment for divergent feeding systems (high- and low-concentrate) and genetic lines (“Control” cows selected to be in line with the UK average genetic merit for kg of fat and protein, and “Select” cows bred to maximize these traits;
). At various stages over the course of the experiment, the components of the diets may change to account for changes in management system, yield levels, and, in the instance of one of these diet changes (home-grown vs. by-products), main ingredient types. However, these experiments have always run for 7 or more years, so are long-term in nature as well. Further, even if ingredient types have changed, we still aim to differentiate the diets by energy content. A summary of the diets is presented in the Supplemental Information (https://doi.org/10.3168/jds.2018-16112). Langhill cows are extensively monitored for production, health, and fertility, and data relevant to the present study were extracted for analysis. Traits of interest included daily milk yield (MY, kg; the sum of morning, midday, and evening milkings), daily BW (kg; the average of 3 daily measurements), weekly fat (FC) and protein (PROT) contents (%), weekly BCS (scored 0 to 5, inclusive, in increments of 0.25, with higher scores indicating fatter animals;
), and daily DMI (kg) for 3 consecutive days a week. Feed composition data included CP (g/kg DM), OM (g/kg DM), organic matter digestibility (OMD, g/kg DM), and ME (MJ/kg DM). Records with DIM between 5 and 305 [i.e., with week in milk (WIM) 1 to 44] were retained for analysis.
Not all observations were collected or available at a daily level; thus, daily estimates for MY, FC, PROT, DMI, BW, and BCS were generated using the random regression model in Equation [1]. This ensured a complete longitudinal data set for each cow, lactation, and DIM (
where yitk is the trait observation of cow i on DIM t (t = 5,…, 305) of lactation k (k = 1,…, 5); µ is the overall mean; Fitk are the fixed effects of diet group (diet rations, both high- and low-concentrate, were modified 2 times between 2005 and 2018 and thus fitted as 6 separate diets, or 6 classes); genetic group (2 classes, i.e., Control or Select line cows); calving age (3 classes: early-, mid-, and late-calving cows); year-by-season of calving interaction (2 calving seasons per year, i.e., 28 classes); year-by-month of record interaction (168 classes); ΣnβnPn(t) represents the fixed regression of the overall lactation curve (βn are fixed regression coefficients) described by fourth-order orthogonal polynomials (Pn) on DIM (t); ΣnλinPn(t) represents the random regression of the individual animal effect (λin are random regression coefficients) described by fourth-order orthogonal polynomials (Pn) on DIM (t) and describes the deviation from the overall lactation curve of the ith individual cow (i = 1, …, 1,101); and εitk is the random residual effect. All effects were fitted within each individual lactation k. All statistical analyses were conducted using the software package ASReml, version 3 (
). This generated a total complete longitudinal data set of 691,232 daily animal test-date records.
The body energy traits of interest in this study, as highlighted above, were effective EI (MJ/d), constituting the daily energy acquired through feeding and calculated using OM intake, CP, and ME; and EB (MJ/d), representing the difference between the amount of energy taken in minus that spent or contributed as a function of MY, FC, PROT, DMI, BW, and BCS. Formulae for the calculation of these energy traits are summarized in Equations [2] and [3]:
EI = 0.001 × (DMI × OM) × EEC
[2]
EB = EI − EEout,
[3]
where DMI × OM is a cow's daily intake of OM (kg); EEC is effective energy content of feed per kilogram of OM; and EEout is the energy used by a cow for milk production, maintenance, and general activity. Full details of these formulae and their underlying components can be found in
. These body energy traits constitute the reference phenotypic values used to calibrate concurrent spectral data described subsequently.
Mid-Infrared Data
Milk Sampling and Analysis
Mid-infrared analysis of milk samples obtained from dairy cows at SRUC's dairy research herd was conducted by National Milk Records and tested on their suite of Foss MilkoScan machines (Foss Analytics, Hillerød, Denmark), based at National Milk Laboratories (Hillington and Wolverhampton), for routine milk recording services. For the purposes of this study, only the MIR analysis of evening samples from Langhill recording were used, following previous analyses carried out on the same herd that highlighted higher accuracies of prediction when using evening compared with morning and midday milk samples (
). The Langhill herd at Crichton Royal Farm is milk recorded weekly; in contrast, the UK milk recording system takes place on a monthly basis.
Pretreatment and Standardization of Mid-Infrared Data
Following MIR analysis of the milk sample, a spectrum of 1,060 data points is generated; each point represents the absorption of mid-infrared light by the milk sample at a particular wavelength within the range of 900 cm−1 to 5,000 cm−1 (
). Before use in generating prediction equations, the spectral data were subject to several pretreatments. First, the transmittance data were converted to a linear absorbance scale using the transformation log10−0.5. Second, spectral data were standardized against a master spectrometer (
Creation of universal MIR calibrations by standardization of milk spectra: Example of fatty acids.
in: Book of abstracts of the 65th annual meeting of the European Federation of Animal Science. Wageningen Academic Publishers,
Copenhagen, Denmark2014: 108
), ensuring that resultant prediction tools developed using such spectral data can be applied to other machines within the same standardization network over time. Finally, a first-derivative calculation was applied to the data by subtracting from each wavepoint value the immediately preceding wavepoint value, similar to the process described in
Alignment of Spectra with Trait to Create Reference Data Set
Mid-infrared spectral data from the Langhill research herd were collected from early 2012 onward. Daily modeled energy traits were then aligned with pretreated spectral data based on animal ID and date of sampling (within ±1 d of milk sampling), such that only those spectra with a corresponding phenotypic record for both energy traits were retained. The final reference (calibration) data set consisted of over 30,594 animal test-date records of body energy trait measures, combined with concurrent pretreated MIR spectra, for lactations 1 to 5.
Data Preprocessing
Biological edits were applied to the base phenotypes used in the calculation of energy traits (MY, FC, PROT, DMI, BW, and BCS), retaining only those records with values within 3 standard deviations (SD) of their respective annual means. Regarding the spectra, and as a form of quality assurance, the Mahalanobis distance was calculated for each composite of spectral points, and only those with a distance less than 20 were retained (
Development of Fourier transform mid-infrared calibrations to predict acetone, β-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network.
Further preprocessing of the calibration data (via data thinning) was required to ensure a balanced representation between low and high reference energy values (
Development of Fourier transform mid-infrared calibrations to predict acetone, β-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network.
Development of Fourier transform mid-infrared calibrations to predict acetone, β-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network.
, such that data were divided across the densest parts of their distribution into equal-sized windows with a proportion (%) of total records selected randomly from each window. Defining the densest area on the bell-shaped distribution within which to apply thinning was a balance between representing as much of the range as possible and not restricting each interval to too few values. Thinning was found to be optimal when selecting 50% of records from each of 10 (15-MJ/d) intervals across the distributions of EB and EI within the ranges −75 ≤ x ≤ 75 and 75 ≤ x ≤ 225, respectively. This resulted in a reduced data set that spanned the same range of energy values as the full data set but with an increased balance between low and high values. Window size, frequency, and proportion of records retained were all varied to optimize the accuracy of the calibration data set.
Finally, records were considered outliers and removed from each of the prediction data sets if their predicted residual sum of squares (PRESS), obtained from the analysis of all reference data, exceeded 1 SD of the mean PRESS value. This process highlighted data from 2 sample dates in particular, where the majority of the predictions (79 and 80%, respectively) generated from spectra on these dates were greater than 3 SD of the mean PRESS; thus, these data were removed. This generated final calibration data sets for EB (n = 11,526) and EI (n = 11,941). A summary of the data is presented in Table 1.
Table 1Summary statistics of the energy balance (EB) and energy intake (EI) calibration data sets
Partial least squares (PLS) regression was used to calibrate the spectral data, using modeled energy phenotypes to generate prediction equations for EB and EI, respectively, based on a subset of the most informative wavepoints; namely 212 specific wavepoints in the range 968 to 2,966 cm−1 (
) and Microsoft SQL (Microsoft Corp., Redmond, WA), using the NIPALS method. The Python script was created to assess the optimal number of latent factors to use in the PLS model, up to a maximum of 20; that is, the number of factors (or variables) that maximize the variation explained while minimizing the prediction error.
Accuracy was measured using the coefficient of determination (R2), generated from the regression of predicted values on actual values by internal split-sample 10-fold cross-validation and external 4-way cross-validation; whereby the data set was split randomly into quarters and the PLS run 4 times, each time one of the quarters serving as the validation file and the remaining quarters as the calibration file. To investigate the effects of having repeated observations per cow, a second external 4-way cross-validation was carried out masking cows (and their corresponding data) instead of records. Accuracy of prediction was assessed using the R2 from validation methods, in addition to the root mean square error and the ratio of SD of the reference values to the standard error of prediction (RPD; a model is deemed to not predict the reference values well if the RPD value is less than or equal to 1.5;
) from the split-sample validation only. The optimal calibration data set generated for each of the 2 energy traits with their concurrent spectra generated linear prediction equations, which were then applied to a data set of national spectra.
Application of Tool to National Data: Phenotypic and Genetic Analysis
Since 2013, National Milk Records has routinely collected spectral data from over 360 commercial farms in the UK as part of an Innovate UK–funded project. From 2016 National Milk Records has retained spectral data generated by all relevant Foss machines; these data encompass MIR-analyzed milk samples from 4,822 herds, 1,484,377 animals, and over 20 million test dates. Since the start of 2018, it has been routine for SRUC to receive almost 15 million spectral data points per day. These data points are fed through a bespoke software platform, which ensures these data are processed (transformed, standardized, first derivative calculations, according to protocols) and subject to prediction equations from the calibration data sets discussed above, which generate a series of MIR-derived phenotypes, such as the energy traits described here. Subsequently, predictions for each of the 2 traits were scaled according to the mean and SD of the reference data set from which they are derived (Langhill), to account for the fact that data from this research herd represents such a broad range of phenotypes due to the nature of the 2 × 2 long-term factorial experiment performed at SRUC on diet and genetic line. To meet the objective of demonstrating the current tool's application and performing genetic analyses, MIR-derived predictions were obtained for 11,717,032 animal test-dates from almost a 4-year period (2013 to 2016).
Phenotypic Analyses
A fixed regression model was used to generate fitted curves for each predicted energy trait for each of the 5 main dairy breeds represented in the national data. These analyses focused on records from HOL, Shorthorn (SHO), Ayrshire (AYR), Jersey (JER), and Guernsey (GUE) cows. To ensure robustness of phenotypes (no attempt was made here to partition the animal genetic component from the permanent environment at this stage), the following edits were applied to the MIR-derived predictions before analysis: records were retained from DIM between 5 and 305, from cows whose age at first calving was between 16 and 48 mo, and at least 5 records per animal per lactation from the first 5 lactations were used. Further, energy predictions within ±3 SD of the overall mean were retained. The final data set totaled 6,888,665 test-date records on 626,237 animals from 4,198 herds. Fitted curves using predicted phenotypes were estimated using the following fixed regression model for each breed individually:
[4]
where yiot is the predicted trait of cow i belonging to breed o at WIM t (t = 1,…, 44); µ is the overall mean; Fiot are fixed effects of month and year of phenotypic record, lactation number, age at first calving (month), herd-year-season of calving, and breed (HOL, SHO, AYR, JER, or GUE); ΣnβnPn(t) represents the fixed regression of the overall lactation curve (βn are fixed regression coefficients) described by third-order Legendre polynomials (Pn) on WIM (t); ΣnβnPn(t) × Eo represents the interaction between the fixed regression of the overall lactation curve and breed and describes the deviation from the lactation curve of the oth individual breed (o = 1,…, 5); ai is the random additive genetic effect of the ith individual cow; and εiot is the random residual effect. All analyses were carried out using ASReml version 3 (
In addition to the data edits highlighted above, HOL records obtained from the national data were further preprocessed by applying the following constraints, such that all cows used for the purposes of the present study (1) were from sires who had sired at least 10 daughters, (2) were the first 250 daughters from each sire, and (3) had at least 8 records per animal per lactation from the first 4 lactations only (defined as lactations 1, 2, and 3 and above). These additional constraints were applied to help represent the population, ensuring that well-proven bulls did not dominate the data set and removing any potential selection bias between young and proven bulls. We did not have as much latitude, in terms of record numbers or overdominance of bulls in the other breeds, to require such additional edits above those applied to all breeds (i.e., number of test records per cow, minimum number of records in a contemporary group, minimum and maximum trait values). These data were linked to a pedigree consisting of over 422,000 individuals. The final data set totaled 1,079,160 test-date records on 93,735 animals from 2,977 herds. Variance components were estimated for each of the 5 breeds via the univariate repeated measures model in Equation [5], with respective pedigree relationship matrices fitted to account for the genetic relationships between individuals within breed:
[5]
where yit is the predicted trait at time t WIM (t = 1,…, 44); µ is the overall mean; Fit are the fixed effects of month and year of phenotypic record, lactation number, age at first calving (month), and herd-year-season of calving; ΣnβnPn(t) represents the fixed regression of the overall lactation curve (βn are fixed regression coefficients) described by fifth-order Legendre polynomials (Pn) on WIM (t); ai and pi are the random additive genetic and permanent environmental effects of the ith individual cow, respectively; and εit is the random residual effect. All analyses were carried out using ASReml version 3 (
Regarding the final data set (reduced, with outliers removed) the R2 obtained from internal split-sample 10-fold cross-validation were 0.77 and 0.75 for EB and EI, respectively. Further, from the external 4-way cross-validation, the R2 were 0.60 and 0.53 for EB and EI, respectively (with external 4-way cross-validation = 0.60 and 0.52 for EB and EI, respectively). The optimal numbers of latent factors selected in these analyses were 20 for each trait (the maximum set). These values, along with additional information regarding the calibration and validation, are presented in Table 2.
Table 2Calibration statistics for energy balance (EB) and energy intake (EI) from 3 independent data sets (full, reduced, reduced with outliers removed)
Nr = number of records in the data set; Nf = number of factors used in prediction equation; RMSE = root mean square error; RPD = ratio of SD of the reference values to the SE of prediction.
Outliers were values with a predicted residual sum of squares (PRESS) exceeding ± 1 SD of the mean PRESS value.
EB
11,525
20
0.77
0.60
0.60
17.96
2.1
EI
11,940
20
0.75
0.53
0.52
22.35
2.01
4Obtained using external 4-way cross-validation (on cows).
1 Nr = number of records in the data set; Nf = number of factors used in prediction equation; RMSE = root mean square error; RPD = ratio of SD of the reference values to the SE of prediction.
2 Obtained using internal split-sample 10-fold cross-validation.
3 Obtained using external 4-way cross-validation (on records).
5 Covering the same range of values as the full data set with a balanced representation between low and high reference values.
6 Outliers were values with a predicted residual sum of squares (PRESS) exceeding ± 1 SD of the mean PRESS value.
The RPD from split-sample validation for EB and EI was greater than 2, suggesting value in these prediction models at discriminating low from high values (
). This is perhaps unsurprising, given common use of some of the Langhill research herd data between both studies.
In this study we have applied a MIR-based prediction for energy traits from HOL cows. In reality, the prediction data set uses divergent lines of HOL for milk production traits interacting with 2 contrasting planes of nutrition (2 × 2), with the nutritional planes altering during the data window (6 diet rations in total). Although this does not reflect the other breeds, it is, by design, meant to cover a wide range of production potentials and feeding systems, and, as such, we expect divergent ranges of energy profiles within and between animals and experimental units. In the development of the model, the number of factors was tested until the minimum PRESS was attained, to choose the optimal number of PLS factors to diminish the possibility of overfitting (
Partial least-squares methods for spectral analyses. 1. Relation to other quantitative calibration methods and the extraction of qualitative information.
). Although a relatively high number of factors were used, the cross-validation error was not excessively high compared with a lower number of factors. Furthermore, it should be noted that the R2 could be, in part, influenced by the fact that EI and EB traits were generated from random regression modeling on the component trait elements (i.e., MY, DMI, BW, BCS) from within the research herd. In the development of a MIR spectral predictive tool, we explored the difference in R2 for MIR predictions of repeated records of BCS, a component trait of EB. When comparing MIR predictions obtained using actual BCS data during calibration with predictions obtained using BCS solutions from random regression modeling, the influence was minimal. This suggested that R2 generated from the compiled EB and EI traits may not necessarily be overinflated. Ideally the cross-validation would be performed on an independent data set, but unfortunately we did not have access to an external set of data to test MIR predictions with our 4-fold validations. The wavepoints studied here are those more informative for the major FA in milk, and changes in the fat profile have shown to be an indicator of EB profiles (
As we continue to assimilate spectra (now over 20 million animal test-dates worth of spectra and MIR-derived phenotypes) and collect and improve the accuracy of energy phenotypes (e.g., the potential use of calorimetric chambers), we should be able to improve system-specific equations and, by incorporating data from across systems, countries, and climates, improve the development of multicountry, multisystem tools (
Application of Tool to Utilize National Data: Phenotypic and Genetic Analysis
Application of the calibration data set to the national spectra generated 11,529,162 and 11,501,992 animal test-dates worth of predictions for EB and EI, respectively. Records were lost (1.6% and 1.84% for EB and EI, respectively) due to quality assurance, applied to ensure that predictions considered outside the ability of these tools to predict usefully were not reported.
Phenotypic Analysis
Figure 1 depicts the trend of EB, EI, and MY across WIM. The profiles of EB and EI generated from fixed phenotypic regression curves are an attempt to model the balance between energy acquired from both feed and body fat mobilization with the energy demands of lactation, growth, and maintenance, with all other fixed effects considered.
Figure 1Predicted milk yield and energy traits from a fixed regression model, plotted by week in milk (WIM) for (A) milk yield, (B) energy balance, and (C) energy intake for first-lactation Holstein Friesian (HOL), Shorthorn (SHO), Ayrshire (AYR), Jersey (JER), and Guernsey (GUE) cows.
Across all breeds, the trend in EB showed an initial recovery from a loss of energy in early lactation, reaching a peak toward the end of mid-lactation and subsequent decline before second calving. The difference in EI profiles between breeds was greatest between JER and HOL cows, suggesting that, at 1 WIM, JER were consuming 12.5 MJ/wk more than HOL. Moreover, at the greatest disparity (22 WIM), JER were consuming 17.4 MJ/wk more than HOL. Overall, JER had the most positive EB profile and the highest EI profile, despite producing the second-to-lowest average MY (34.4 kg/wk). This is in contrast to studies that have suggested JER cows to be more energy efficient than HOL (
). It has been shown, however, that JER cows eat more per unit of BW than HOL and, hence, have a higher intake capacity, likely due to their relatively larger digestive tract (
). Moreover, JER cows tend to produce milk with higher FC, which is energy dense, highlighting the fact that methods of measuring efficiency can be different. In contrast to JER, GUE cows showed the most negative EB profile across first lactation and the second-lowest EI profile in early lactation, consistent with their delivery of the lowest MY among all breeds (34.2 kg/wk).
At the onset of lactation, EB levels in JER cows were 6.7 MJ/wk higher than in GUE cows. Jersey and SHO animals tended to experience their nadir in EB at the start of lactation (1 WIM), much earlier than AYR (7 WIM), GUE (9 WIM), and HOL cows (10 WIM). The initial rate of EB loss in early lactation was fastest for HOL (−0.19 MJ/wk) and slowest for SHO cows (−0.004 MJ/wk). Energy intake was observed increasing from the start of lactation to a maximum at wk 13 (SHO), 15 (AYR), 22 (HOL and JER), and 25 (GUE). Despite the same duration to maximum EI of 22 wk observed in JER and HOL, JER showed the fastest rate of gain from start of lactation (0.22 MJ/wk) and HOL the lowest (0.01 MJ/wk). Despite the fact that HOL cows displayed the fastest decline of EB to nadir in early lactation and had the lowest rate of gain in EI to its maximum, they produced on average the greatest amount of milk of all 5 breeds (36.7 kg/wk). Moreover, HOL cows reached their peak EI in the same amount of time as JER cows reached their peak MY (22 WIM). The duration between peak MY and peak EB was second-lowest for HOL at 24 wk; strong selection for yield may have predisposed this breed to more readily utilize body energy through mobilization of their fat reserves to support lactation (
Genetic evaluations of dairy bulls for daughter energy balance profiles using linear type scores and body condition score analyzed using random regression.
), supporting the lower EI profile for HOL shown here despite highest MY.
Recovery from the initial nadir of EB was fastest in JER animals (0.13 MJ/wk), followed by HOL (0.07 MJ/wk), and lowest for SHO (0.04 MJ/wk). When looking at both energy trends across lactation, it is apparent that the peak in EI occurred 1 wk before a peak in EB for JER animals, 2 wk before for HOL, 3 wk before for GUE, 7 wk before for SHO, and 8 wk before for AYR. Furthermore, peak MY occurred 1 wk before the nadir of EB in early lactation for SHO and JER, 7 wk before in AYR, 9 wk before in GUE, and 10 wk before in HOL. The most positive EB across all breeds was achieved in mid-lactation. Maximum EB levels were observed at wk 20 for SHO, 23 for AYR and JER, 24 for HOL, and 28 for GUE. Shorthorn and AYR cows showed similar energy profiles and had similar average MY (35.2 kg/wk and 35.3 kg/wk for SHO and AYR, respectively). Cows belonging to these 2 breeds tended to be larger and heavier than their counterparts from the other breeds, and as such it may be that a greater proportion of their available energy goes into processes such as growth, at the expense of production. Such is consistent with the negative relationship between body size and efficiency (
Following the peaks observed in mid-lactation, the rate of EI decline leading to subsequent calving was fastest for GUE (−0.62 MJ/wk), followed by HOL (−0.6 MJ/wk), and slowest for SHO (−0.29 MJ/wk). This tail-off toward the end of the milk and energy distributions for all breeds may reflect the influence of fewer records available toward the end of lactation in the random regression rather than biological differences in the rate of decline between breeds. Alternatively, it may reinforce the idea that EB represents a different biological trait at different stages of lactation with differing efficiencies (
), such that prediction tools may need to be lactation-stage specific. With the reference data and concurrent spectra available from the Langhill herd, developing and testing such tools is highly plausible. It is important to bear in mind, however, that the tools established in the present study, and any future tools developed using the Langhill research herd, are based on one breed only (HOL), and, as such, validation data from other breeds would be necessary to assess the accuracy of tools applied across breeds before any deployment in the field.
Genetic Analysis
Estimates of additive genetic and permanent environmental variance, heritability, and repeatability of EB and EI for the different breeds are reported in Table 3. The heritability of EB for all breeds in the present study is within the range of other estimates from the literature (
, h2 = 0.10 to 0.34). However, unlike previous studies, the estimates here are based on MIR-derived predictions over all lactation, which have tended to be lower (
, who estimated the h2 of EB and EI as 0.1 and 0.2, respectively. Further, our h2 estimates of EB and EI are still greater than the estimates of some functional traits incorporated in national and international selection criteria and, therefore, may be more effective for use in selection criteria due to their significant genetic correlation with such traits (
Estimated additive genetic variation for these traits, across all breeds, highlights that energy partitioning is under at least some genetic control and, therefore, can facilitate the selection of animals with appropriate energy profiles that could contribute to the genetic improvement of fertility (
Development of Fourier transform mid-infrared calibrations to predict acetone, β-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network.
); as such, going forward the data should be modeled using test-day or random regression models and accounting for fewer records and error variance in early and late lactation when cows are drying off.
Lastly, optimization via genetic selection can only go so far; management of appropriate nutrition and levels of concentrate (see
, for greater effects of low-concentrate diet on body energy loss) and their cost will also determine the effectiveness and profitability of cows with sustainable energy profiles.
The prediction of energy traits from milk in this study relies on several assumptions: that the research herd from which reference data were sourced is representative of the herds in which it is intended to predict (
)—is accurate; that changes in the energy profiles of dairy cows leave some sort of molecular signal in their milk; and that we are able to detect this signal and calibrate it. We have demonstrated that by exploiting the relationship between the MIR absorption spectra (ultimately, relative milk fat and protein molecular concentrations) and concurrent phenotypes, a robust tool for energy traits can be applied to herds in the UK with cost and labor advantages over more manual (e.g., measurement of DMI) or subjective methods (e.g., BCS) or those that do not accurately reflect underlying energy stores (e.g., BW).
CONCLUSIONS
The automated prediction of energy traits deployed via routine milk recording and as part of DHI schemes could provide farmers with a low-cost and, more importantly, noninvasive tool to guide nutritional and transition cow management within herds. This would enable farmers to maximize number of lactations and reduce early culling, as well as optimize feed efficiency (due to the strong genetic correlation of EB with residual feed intake;
) while ameliorating the effects of health and fertility issues associated with NEB. This could be particularly effective if contemporary selection objectives (primarily for yield) favor cows that mobilize their body tissue and show compromised fertility (
The authors gratefully acknowledge collaboration with National Milk Records (Martin Busfield, Victoria Hicks; Harrogate, UK) and Marks and Spencer (London, UK); funding from Innovate UK (previously TSB; Swindon, UK; grant no: TP12483-81287) and Biotechnology and Biological Sciences Research Council (BBSRC, Swindon, UK; grant no: BB/K002260/1); and cooperation with the participating farms. Ian Archibald (SRUC, Edinburgh, UK) is acknowledged for curating the Langhill data, Jenni Flockhart (SRUC) for details of the Langhill diets, and all staff at Crichton farm for collecting samples and managing animals. The Langhill experiment at Crichton Dairy Research Centre and authors MC, SJD, and EW are supported by the Scottish Government Rural Affairs, Food and the Environment (RAFE, Edinburgh, UK) Strategic Research Portfolio 2016-2021.
Genetic evaluations of dairy bulls for daughter energy balance profiles using linear type scores and body condition score analyzed using random regression.
Development of Fourier transform mid-infrared calibrations to predict acetone, β-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network.
Creation of universal MIR calibrations by standardization of milk spectra: Example of fatty acids.
in: Book of abstracts of the 65th annual meeting of the European Federation of Animal Science. Wageningen Academic Publishers,
Copenhagen, Denmark2014: 108
Partial least-squares methods for spectral analyses. 1. Relation to other quantitative calibration methods and the extraction of qualitative information.