|
|
||||||||



* Department of Animal Science, Texas A&M University, College Station 77843-2471
Department of Animal Science, Cornell University, Ithaca, NY 14853
Elanco Animal Health, Guadalajara, Jalisco 44620, Mexico
1 Corresponding author: luis.tedeschi{at}tamu.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: body condition score fat mobilization fat repletion modeling
| INTRODUCTION |
|---|
|
|
|---|
Body weight changes reflect the use of energy reserves either to supplement ration deficiencies during early lactation or to store energy consumed above the requirements (Moe et al., 1972; NRC, 2001). The BW gain and loss after maturity is similar in composition to changes during growth (NRC, 2001) and can be used to predict changes in energy balance over the reproductive cycle (Reynoso-Campos et al., 2004). However, most dairy and beef producers monitor BCS changes in cows to manage energy reserves because frequent measurements of BW are not feasible under practical conditions. The Cornell Net Carbohydrate and Protein System (CNCPS; Fox et al., 2004) and the NRC (2000 NRC (2001) models use the body reserves model as devised by Fox et al. (1999), which was developed from data on the chemical body composition and BCS of 106 mature beef cows of diverse breed types and BW. As applied to dairy cattle, the model was evaluated with the data of Otto et al. (1991) and accounted for 95% of the variation in body fat, with only a 1.6% bias (Fox et al., 1999). The model predicted 80 kg of BW change per BCS compared with 84.6 kg observed in Holstein cows slaughtered over the range of dairy BCS.
For lactating dairy cows, the CNCPS and NRC models estimate energy and protein requirements for maintenance and pregnancy, and the amount remaining above intake is used to estimate ME- and MP-allowable milk production, respectively (Fox et al., 2004). The changes in BCS are not accounted for in predicting ME and MP balances. The objective of this study was to develop and compare 2 empirical models to eliminate biases in energy retention by adjusting the predicted ME- and MP-allowable milk production after consecutive changes in observed BCS have been accounted for.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Milk Energy and Body Content of Energy and Protein.
For both empirical models, the energy contained in milk production is computed using milk fat and milk true protein contents, as described by the NRC (2001). This energy in milk (MkE), as shown in Equation [1], is assumed to be the NEL:
![]() | [1] |
where MkEi is the energy content of milk (Mcal/kg), MkFi is the milk fat content (g/100 g), MkTPi is the milk true protein content (g/100 g), and the subscript i is the ith time period.
The BCS is a 5- (Wildman et al., 1982) or 9-point (Cantrell et al., 1981; Herd and Sprott, 1986) scale system that is highly related to body fat in cows (Houghton et al., 1990; Buskirk et al., 1992; NRC, 2000, 2001). Other scale systems that are used around the world (CSIRO, 1990) can be interconverted. Because the body reserves model used by the NRC (2001) is based on that developed by the NRC (2000), a BCS scale of 1 to 9 is used. Equation [2] is used to convert a BCS scale of 1 to 8 (BCS[18]), as used by the Commonwealth Scientific and Industrial Research Organisation (CSIRO, 1990), to a BCS scale of 1 to 5 (BCS[15]) and Equation [3] converts BCS[15] to a BCS scale of 1 to 9 (BCS[19]), and vice versa. As adopted by the NRC (2000 NRC (2001), shrunk BW (SBW) is computed from BW as shown in Equation [4] and empty BW (EBW) is estimated from SBW as shown in Equation [5], which is used to predict body reserves:
![]() | [2] |
![]() | [3] |
where BCS[15],i is the BCS on a scale of 1 to 5, BCS[18],i is the BCS on a scale of 1 to 8, and BCS[19],i is the BCS on a scale of 1 to 9,
![]() | [4] |
![]() | [5] |
where SBWi is shrunk BW (kg) and EBWi is empty BW (kg).
Empirical Reserves Model 1.
This model is based on the equations published by Fox et al. (1999) in which BCS[19] and EBW are used to compute the amount of body fat (TF; Equation [6]) and protein (TP; Equation [7]):
![]() | [6] |
![]() | [7] |
where EBWi is empty BW (kg), TFi is the amount of body fat (kg), BCS[19],i is the BCS (on a scale of 1 to 9), and TPi is the amount of body protein (kg).
For mature lactating cows, a change in BW does not necessarily indicate changes in tissue reserves, and vice versa. Andrew et al. (1994) and Gibb et al. (1992) analyzed slaughter data of dairy cows and reported as much as 40% variation in energy with no change in BW. This is likely because the gut fill varies from 2.5 (Komaragiri and Erdman, 1997) to 4 kg/kg of increase in DMI (Chilliard et al., 1991), which may offset the weight loss attributable to tissue mobilization by an increase in DMI during early lactation. Because of this disconnection between actual BW and energy reserves, ERM1 uses BCS changes to estimate EBW and energy reserves.
As discussed by Fox et al. (1999), the database used at the NRC (2000) to develop the body reserves model indicated that the mean BW change associated with a BCS change was equivalent to 6.85% of the mean BW. The Commonwealth Scientific and Industrial Research Organisation (1990) uses 8% of the standard reference weight per change in BCS[18], which is equivalent to 7% per change in BCS[19]. Therefore, a weight adjustment factor (WAF; Equation [8]) is computed from the BCS. Adjusted EBW values (aEBW; Equation [9]) are then computed for all other periods (i
2) based on their respective WAF values (Equation [8]), which are a function of the measured BCS for each period:
![]() | [8] |
![]() | [9] |
where BCS[19],i is the BCS (on a scale of 1 to 9).
The (EBWi=1/WAFi=1) factor in Equation [9] computes the expected BW at BCS 5. The aEBW for each period (aEBWi) is then used to assess the variation in tissue energy, which is added or subtracted from the tissue energy computed using the previous aEBW (aEBWi1) and Equations [6] and [7].
Empirical Reserves Model 2.
This model is based on the equations derived by Otto et al. (1991) to compute the proportion of fat (Equation [10]) and protein (Equation [11]) in the 9th to 11th rib section of Holstein cows:
![]() | [10] |
![]() | [11] |
where EE is ether extract (%) and BCS[15] is the BCS on a scale of 1 to 5.
Conversion of the fat and protein contents from the 9th to the 11th rib to EBW is performed using the relationship between 9th and 11th rib and carcass composition as developed by Hankins and Howe (1946; Equations [12] and [13], respectively) and from the carcass to EBW as derived by Garrett and Hinman (1969; Equations [14] and [15], respectively):
![]() | [12] |
![]() | [13] |
![]() | [14] |
![]() | [15] |
where EE is ether extract (%) and EBW is empty BW.
Equation [16] was derived by combining Equations [3], [10], [12], and [14], and was solved for fat in the EBW. Similarly, Equation [17] was derived by combining Equations [3], [11], [13], and [15], and was solved for protein in the EBW:
![]() | [16] |
![]() | [17] |
Like the ERM1, the ERM2 assumes that changes in BW may reflect changes in gut fill and may not represent a true change in tissue. Therefore, BCS changes are used to compute changes in the tissue as a function of EBW. However, unlike the ERM1, the ERM2 relies on a fixed variation of EBW per change in BCS based on the analyses of Otto et al. (1991), who found that each unit of BCS change was associated with a 56-kg change in live BW. When converted to EBW, assuming a gut fill of 15%, this change was 47.7 kg of EBW. Therefore, WAF values are compute for each interval based on the initial EBW and changes in the BCS using Equation [18] and aEBW are computed based on EBW and WAF values, as shown in Equation [19]:
![]() | [18] |
![]() | [19] |
where BCS[15],i is the body condition score (on a scale of 1 to 5).
For both models (ERM1 and ERM2), total energy (TE; Equation [20]) is computed from TF and TP multiplied by their respective heat of combustion. For growing animals, the heat of combustion of fat has been assumed to be 9.367 Mcal/kg (Blaxter and Rook, 1953) and protein has varied from 5.554 to 5.686 Mcal/kg (Lofgreen, 1965; Garrett, 1987). For dairy cows, Andrew et al. (1994) derived values of 9.2 and 5.57 Mcal/kg for fat and protein, respectively. The NRC (2000 The NRC (2001) has adopted the growing animal values (9.367 and 5.554 Mcal/kg). Because both estimates are nearly identical, the NRC (2000, 2001) values are used in these models:
![]() | [20] |
where TFi is the amount of body fat (kg), TPi is the amount of body protein (kg), TEi is the total energy (Mcal), and the subscript i is the ith period.
Assessing Changes in Body Energy.
The TE of the first period, which uses the current EBW of the cow, remains unchanged; however, the TE of subsequent periods is computed with the aEBW and Equation [20]. The variation in TE is computed using Equation [21]:
![]() | [21] |
where
TEi is the change in total energy (Mcal), and subscripts i and i1 represent actual and previous TE values, respectively.
When the
TE value is negative, there is a mobilization of reserve energy for milk production. The amount of milk production supported from mobilized reserves is added to the diet-allowable milk production. When the
TE value is positive, the intake of energy is greater than the energy required for milk production. In this case, part of the available energy is used for reserve deposition rather than milk production; therefore, the amount of energy deposited has to be used to reduce the diet-allowable milk production. The variation in body protein is computed using Equation [22]:
![]() | [22] |
where
TPi is total protein variation (kg), and the subscripts i and i1 represent actual and previous TE values, respectively.
The mobilized (
TE < 0 or
TP < 0) or deposited (
TE > 0 or
TP > 0) energy and protein are converted to milk equivalents using efficiencies of energy and protein conversion factors, as described in the next section.
Energy and Protein Efficiencies.
The coefficients of energy interconversion used in our model were derived by Moe et al. (1970) using a multiple regression analysis of respiration chamber data from 126 and 224 lactating dairy cows in negative and positive energy balances, respectively. The confidence interval of the coefficients reported by Moe et al. (1970) indicated a statistical difference (P < 0.05) between the efficiency of ME to net energy of reserves (NER; 72.6%), ME to NEL (63.5%), and NER to NEL (84%), suggesting there are significant differences in the metabolism of lactating cows at a positive and negative energy balance that affect the efficienct use of energy.
When
TE < 0, it indicates a negative energy balance, and energy reserves (NER) were used for milk production. An efficiency for NER to NEL of 82% is generally used (Moe, 1981; Fox et al., 1999; NRC, 2001). The Commonwealth Scientific and Industrial Research Organisation (1990) and the Agricultural Research Council (ARC, 1980) assume an efficiency of 84%, which is supported by the study of Vermorel and Bickel (1980). Similarly, Moe et al. (1970) observed an efficiency of 84% for lactating cows in negative energy balance, which was used in these models. Analysis of the variation in the Moe et al. (1970) data indicated that the true efficiency value was between 81.7 and 86%, assuming
= 5% and a less rigid combination of coefficients; we used the average efficiency of 84% in our model. The milk from mobilized reserves is added to the predicted diet ME-allowable milk using Equation [23]:
![]() | [23] |
where
TE is tissue energy variation (Mcal NEL/d),
Milk is milk variation (kg/d), and MkE is energy content of the milk (Mcal of NEL/kg).
A
TE > 0 indicates a positive energy balance in which diet energy was used for reserves rather than milk production. Therefore, the first step is to convert the NER to ME, the second step is to convert this amount of ME to NEL, and finally this NEL is divided by milk energy to compute the amount of milk that was not produced (Equation [24]). Commonly, an efficiency of ME to NER of 75% and ME to NEL of 64.4% are assumed (Moe, 1981; Fox et al., 1999; NRC, 2001). Moe et al. (1970) reported that lactating cows in positive energy balance had an efficiency of 63.5% for ME to NEL. An analysis of the possible combinations of the confidence intervals of the coefficients reported by Moe et al. (1970) suggested that the true efficiency value was between 61.2 and 65.9%, assuming
= 5%. Moe et al. (1970) reported an efficiency of 72.6% for ME to NER for lactating cows in positive energy balance; the confidence interval was 67.3 to 78.7% (
= 5%). Therefore, for
TE > 0 (positive energy balance), we assumed 63.5% for ME to NEL and 72.6% for ME to NER (Moe et al., 1970) for our models. We then calculated the amount of milk from this amount of TE and subtracted it from the predicted diet ME-allowable milk using Equation [24]:
![]() | [24] |
where
TE is tissue energy variation (Mcal of NEL/d),
Milk is milk variation (kg/d), and MkE is the energy content of the milk (Mcal of NEL/kg).
The mobilization of body tissue will also release AA that can be used directly and indirectly for milk production. The indirect form is through the recycled N into the gastrointestinal tract; this form is accounted for by the CNCPS model (Fox et al., 2004). The direct form, which is discussed here, is related to the incorporation of the AA produced from mobilized reserves into milk protein. In lactating sows, daily mobilization of protein was previously shown to be approximately 11 to 63 g and 45 to 195 g for first and third lactations, respectively, with an efficiency of use of 68.9 to 73.4% for milk production. These values are much greater than the values for efficiency of use of feed protein for milk production (42.7 and 46% for the first and third lactations, respectively; Lahrssen, 1988).
Dairy cows can mobilize between 7 and 13 kg of body protein within the first 2 to 4 wk of lactation (Journet et al., 1983). Increasing the quantity of RUP fed prepartum had a positive effect by decreasing the loss of BW and increasing the milk protein content (Van Saun et al., 1993).
Information regarding the efficiency with which these AA are utilized for milk production in cattle is scarce. Not all mobilized AA appear in milk protein, and the AA profile of muscle does not match the milk protein profile; therefore, adequate accounting for this phenomenon is needed (McNamara, 2000). The efficiency of use of mobilized AA for milk production is likely to vary depending on the nutritional status of the animal and the N balance. With N-deficient diets, the efficiency of use of N can be as great as 75% (NRC, 1985). Ruiz et al. (2002) indicated that using an efficiency of MP to net protein for milk production (NPL) of 75% resulted in no bias in the CNCPS-predicted MP-allowable milk. The Commonwealth Scientific and Industrial Research Organisation (1990) suggests that the mobilized net protein (NPR) has an efficiency of 80% for NPL. Therefore, when
TP < 0 (mobilization of body protein), we assume an efficiency of NPR to NPL of 80% and the amount of milk attributable to protein mobilization (Equation [25]) is added to the MP-allowable milk production:
![]() | [25] |
where
TP is tissue protein variation (g/d), MkTPi is the milk true protein content (g/100 g), and
Milk is milk variation (kg/d).
This efficiency depends on the profile of the AA of the mobilized protein and on the profile of the milk protein (Newbold, 1994). Ionophores may have an impact on sparing AA from gluconeogenesis (Tedeschi et al., 2003), leading to a greater efficiency of use.
Like energy, a
TP > 0 indicates a positive protein balance, in which protein was used for reserves rather than milk production. Therefore, to adjust the diet MP-allowable milk, we convert NPR to MP, which is converted to NPL. There is little information on the efficiency of conversion of MP to NPR and vice versa, likely because of the difficulty in measuring it. We assume an efficiency of 70% for our model. The efficiency of MP to NPL is assumed to be 64% in the Institut National de la Recherche Agronomique (1989) system, 65% in the beef NRC (2000) system, 67% in the dairy NRC (2001) system, and 68% in the Agricultural and Food Research Council (AFRC, 1992) system. The 65% efficiency was used in this model.
The amount of milk change attributable to deposition of tissue protein is computed as shown in Equation [26], which is used to adjust the diet MP-allowable milk:
![]() | [26] |
where
TP is tissue protein variation (g/d), MkTPi is milk true protein content (g/100 g), and
Milk is milk variation (kg/d).
Adjusting Predicted ME- and MP-Allowable Milk.
Figure 1
depicts a flowchart of the calculation for the 2 methods used to adjust the predicted ME- and MP-allowable milk production. In method 1 (mean), the adjusted milk production is computed for each BCS measured for the same lactating cow within a given time period, whereas method 2 (period) computes the adjusted milk production based on the first and last BCS within a given time period.
|
Evaluation of the Model Predictions
Assessment of the adequacy of mathematical models is only possible through the use of a combination of several statistical and empirical analyses and proper investigation regarding the purposes of the model initially conceptualized (Tedeschi, 2006). Therefore, several techniques were used to evaluate and compare the models in this study. Briefly, linear regression analysis observed on model-predicted values was conducted to identify precision and accuracy in the prediction of the rate of milk production. The coefficient of determination (r2; Neter et al., 1996), confidence intervals for the parameters (Mitchell, 1997), and a simultaneous test for the intercept and slope (Dent and Blackie, 1979; Mayer et al., 1994) were used. The deviation plot (model-predicted minus observed values against observed values) were used to study the behavior of model prediction compared with observed values (Mitchell and Sheehy, 1997). All residual analysis for outliers (extreme and influential points), homoscedasticity, and normal distribution assumptions were performed on the residual (regression-predicted minus observed values) against regression-predicted values (Neter et al., 1996).
Additional techniques were also used, as discussed by Tedeschi (2006), including accuracy from the concordance correlation coefficient (CCC by Lin, 1989; Cb by Nickerson, 1997; and A
by Liao, 2003), mean bias (Cochran and Cox, 1957), and mean square error of prediction (MSEP; Bibby and Toutenburg, 1977). The MSEP values were expanded into 3 fractions to represent errors in central tendency, errors attributable to regression, and errors attributable to disturbances (or random errors), that is, unexplained variance that could not be accounted for by the linear regression (Theil, 1961). Equation [27] has the equation modified by Tedeschi (2006):
![]() | [27] |
where MSEP is the mean square error of prediction, Xi is the ith model-predicted value, Yi is the ith observed value, s2 is the variance associated with observed and model-predicted values, and r2 is the coefficient of determination.
The 3 terms in Equation [27] represent errors in central tendency (mean bias), errors attributable to regression (systematic or slope bias), and errors attributable to disturbances (or random errors), that is, unexplained variance that cannot be accounted for by the linear regression. Each error term is commonly evaluated as a proportion of the total MSEP, thus indicating which terms have a greater influence in the MSEP. Detailed information about all these assessment techniques and the software used to perform the model comparison analyses are discussed in Tedeschi (2006).
The MSEP was used to compare the model accuracy among different combinations of nutritional models (CNCPS vs. NRC), reserves models (ERM1 vs. ERM2), and methods of calculation (mean vs. period) compared with observed values. Equation [28] was used to compute the difference between MSEP (
MSEP) of any 2 sets of model predictions for each data point, and Equation [29] was used to compute the mean and variance of the
MSEP. A t-test was conducted to statistically verify the difference of
MSEP from zero. In a practical application, if the standard deviation of
MSEP were greater than the
MSEP value, it would indicate that the
MSEP was not different from zero:
![]() | [28] |
![]() | [29] |
where
MSEP is the difference between the 2 sets of model predictions for each data point, f(X)i are model-predicted values using one set of predicted values and g(X)i are the model-predicted values using another set of predicted values, Y is the observed milk production,
is the mean of
MSEP, and
2MSEP is the variance of
MSEP.
A meta-analysis was conducted according to Equation [30] (Littell et al., 1999) and indicated no random influence of the trials on the intercept (P = 0.20) and slope (P = 0.36); therefore, the trials were pooled in the analysis:
![]() | [30] |
where Yij is the observed milk production (kg/d); Xij is the model-predicted milk production (kg/d); eij is the random error, independently, identically, and normally distributed with mean zero and variance
2; and a and b are fixed variables, with variancecovariance represented by
. A variance component structure was used in this analysis based on the 2 REML log likelihood values.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
Table 2
lists the energy content of repletion and depletion of body reserves for changes in the BCS of cows for 5 different SBW. These calculations were made with Equations [6] to [9] for ERM1, Equations [16] to [19] for ERM2, and Equation [20] for both models. The energy content per change in kilogram of SBW (
SBW) varied from 4.3 to 8.1 Mcal/
SBW for ERM2 and 2.8 to 6.3 Mcal/
SBW for ERM1. These values are in agreement with several studies (Houghton et al., 1990; Buskirk et al., 1992). Buskirk et al. (1992) reported values ranging from 2.8 to 7.8 Mcal/
SBW. However, for ERM2, within a BCS the megacalories per
SBW did not change with SBW, whereas for ERM1 this value changed with SBW. The mean of ERM1 is greater than the value of 3.6 Mcal/
SBW reported by Schwager-Suter et al. (2001b).
|
SBW obtained for ERM1 (Table 2
Figure 3
depicts the relationship between observed milk yield and milk yield predicted by the NRC (2001) and CNCPS (Fox et al., 2004) models without adjustment for changes in BCS. The NRC (2001) model had greater precision (based on the coefficient of determination of the regression between observed and model-predicted values) than the CNCPS model (0.90 vs. 0.85, respectively; Table 3
) but had a greater mean bias (12.3 vs. 5.34%, respectively; Table 3
). The root of MSEP values of the unadjusted (original) predictions by the CNCPS and NRC models differed (4.77 vs. 5.42, respectively; P = 0.03). Within the CNCPS and NRC predictions, the mean method for both ERM1 and ERM2 had similar MSEP (P > 0.20) and were lower (P < 0.01) than the unadjusted (original) and the period method MSEP values (Table 3
). Analysis of the MSEP decomposition indicated that the mean and systematic biases in the NRC (2001) predictions were a greater source of disparity than in the CNCPS prediction and that the random errors were greater in the CNCPS than in the NRC (2001) predictions. These findings suggest that the NRC (2001) predictions tended to be more homogeneous (precise) than those from the CNCPS model but were less accurate, as also indicated by the concordance correlation coefficient analysis (Table 3
). We concluded that both models adequately predicted the first-limiting ME-or MP-allowable milk.
|
|
BCS; Garnsworthy, 1988; Komaragiri and Erdman, 1997) and depends on several factors (e.g., age, breed, plane of nutrition, parity, etc.). This might explain why a simple adjustment (method 2 of calculation) can be more accurate than sequential adjustments (method 1 of calculation).
Figure 4
shows the relationship of observed milk production and first-limiting ME- or MP-allowable milk adjusted for changes in BCS. In agreement with Table 3
, the adjustment shifted the prediction to the right, decreasing the overall underprediction of both models. However, the improvements in model precision and accuracy were small.
|
BCS) across BCS categories (Garnsworthy, 1988; Komaragiri and Erdman, 1997). We expect that ERM1 is a more robust model than ERM2 for situations in which the animals are different from the ones used in the study by Otto et al. (1991). Schwager-Suter et al. (2001a) found that changes in BW were a better indicator of body tissue change than changes in BCS. The ERM2 can be implemented with a variable change in BW rather than a fixed value of 47.7 kg of EBW, as reported by Tennant et al. (2002). However, as we have stated, variation in BW is difficult to assess at the farm level because of changes in gut fill. As our model evaluation indicated, BCS can be used to account for energy mobilization and repletion.
In our model, we assumed fixed values for the fat and protein content of mobilized and replenished tissue. Williams et al. (1989) suggested that the contribution of energy from fat and protein (Equation [20]) be adjusted by stage of lactation. The authors proposed multiplicative factors for fat and protein (1.4 and 0.6, respectively) to account for differences in the amount of fat and protein mobilized during early lactation, as opposed to the composition of gain during late lactation.
More complex models of adipose tissue biochemistry have been discussed by McNamara (2000) and Baldwin (1995, chapters 12 and 16). These models use substrate saturation and enzyme kinetics or mass action to modulate the mechanisms of lipogenesis and lipolysis biochemically. Simpler models that work within the NRC and CNCPS frameworks are needed to account for BCS changes of lactating cows and predict production responses with different feeding systems given an animals potential performance, energy and protein supply, environmental conditions, and body reserve management strategies. Our results indicate that the model presented can be used to account for changes in BCS in formulating diets on farms and in evaluating differences in milk production with different experimental diets.
Received for publication August 26, 2005. Accepted for publication July 11, 2006.
| REFERENCES |
|---|
|
|
|---|
) of lactating beef cows. J. Anim. Sci. 70:38673876.[Abstract]This article has been cited by other articles:
![]() |
J. R. Roche, N. C. Friggens, J. K. Kay, M. W. Fisher, K. J. Stafford, and D. P. Berry Invited review: Body condition score and its association with dairy cow productivity, health, and welfare J Dairy Sci, December 1, 2009; 92(12): 5769 - 5801. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |