|
|
||||||||
Centre for Nutrition Modelling, Department of Animal and Poultry Science, University of Guelph, Ontario, Canada, N1G 2W1
1 Corresponding author: jellis{at}uoguelph.ca
| ABSTRACT |
|---|
|
|
|---|
Key Words: maintenance energy expenditure energy balance dairy cow
| INTRODUCTION |
|---|
|
|
|---|
In the process of developing approaches to predict daily DMI in a dynamic simulation of dairy cow performance, it is important to first evaluate and correct, if necessary, the energy accumulation in body stores implied from observed DMI and milk production rates. McNamara and Baldwin (2000) found, using a model of the metabolic transactions in the lactating cow (Baldwin, 1995), that in response to dietary changes, body fat accumulation was either over- or underpredicted in simulations of more than a few weeks. They concluded that the models precision in predicting long-term dynamic changes in energy-utilizing reactions was inadequate. Heuer et al. (2001) also observed, in a model designed to predict herd mean NEL energy balance for the first 12 wk of lactation (WOL), that the predicted NEL balance was generally higher than calculated NEL balance in all weeks of lactation. However, at the individual cow level, they observed that the standard deviation was often larger than the mean difference.
A companion paper (Ellis et al., 2006) describes the performance of selected DMI prediction equations in a dynamic model of energy flows in the lactating dairy cow. The test variable used to compare model predictions with observations, and thereby evaluate the DMI prediction equations, is the integrated instantaneous NEL balance, expressed as weekly BW throughout lactation. Published data on DMI, BW, and milk production over a full lactation were collected for this goal. To use a predicted NEL balance as the test variable, it was necessary to first evaluate NEL balances obtained from the observations alone. The purpose of the work reported here was to determine the degree of deviation of predicted from observed BW according to NEL balance calculated from observations in the collated data sets and to determine an appropriate correction of the bias for subsequent DMI prediction evaluations (Ellis et al., 2006). The bias was placed in maintenance expenditures so that a time-independent, as well as a time-dependent, equation relating maintenance expenditures to WOL could be developed.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
![]() | [1] |
The NEL equivalents of body mass, according to NRC (1988), were assumed so that BW change (
BWi) in kilograms per day was 0.203 times NEL balance (Mcal/d; equation 1) if the NEL balance was negative and 0.195 times NEL balance if positive. The older NRC (1988) factors were used because information on BCS, required for calculating the newer NRC (2001) factors, was not available from the data sets. Assuming a BCS of 3.5, BW change would be 0.196 and 0.171 times NEL balance for negative and positive balances, respectively, according to NRC (2001). Analysis showed that the model was relatively insensitive to a 15% change in NEL equivalents of body mass, and so the single NRC (1988) factors were used.
Predicted BW of cows starting at wk 3 of lactation was obtained from the initial observed BW at wk 2 of lactation adjusted for the predicted BW changes for that week (
BWn). Similarly, the BW of cows at progressive WOL (BWn+1) was obtained from the previous predicted BW (BWn) adjusted for that weeks predicted BW changes (
BWn):
![]() | [2] |
Maintenance NEL expenditures were calculated as 0.08 Mcal of NEL/kg of BW0.75 (NRC, 2001) where BW was the predicted, not the observed value, to simulate a dynamic system. The NEL intake was calculated as observed DMI (kg/d) multiplied by the energy content of the feed, estimated for each WOL by the NRC (2001) computer formulation program, from information on diet composition provided in each publication. Assumptions were that pregnancy commenced at 90 DIM, age at first calving was 24 mo, calving interval was 12 mo, BCS was 3.5, calf birth weight was 45 kg, and no animals were grazing. If lactation number was greater than 1, mature BW was set to the observed value at the onset of lactation. For first-lactation cows, mature BW was calculated as 1.25 times observed BW at the start of lactation (NRC, 2001). Milk production was entered as kg of FCM, 4% milk fat, 3% true protein, and 4.78% lactose.
The NRC (2001) feed library was used, with modifications, when chemical composition of the feed was provided. When diet NEL values were provided in the publications, the NRC (2001) program was still used to calculate NEL. This ensured that all diet NEL values were obtained from the same set of assumptions.
Net energy (NEL) for milk production was the observed FCM production (kg/d) multiplied by 0.749 Mcal of NEL/kg (NRC, 2001).
Extra energy expenditures by first-lactation heifers was calculated according to the NRC (2001) for a BCS of 3.5 as 0.0103 Mcal of NEL/kg of SBW0.75, where SBW is shrunk BW = 0.96 x BW.
Statistical Analyses
Mean square prediction error (MSPE) for each of the 21 data sets was calculated as:
![]() | [3] |
where n is the number of observations, Oi is the observed value, and Pi is the predicted value. Square root of the MSPE, expressed as a proportion of the observed mean, gave an estimate of the overall prediction error. The MSPE was decomposed into random error, error due to deviation of the predicted vs. observed regression slope from unity, and error due to overall bias (Bibby and Toutenburg, 1977).
Predictions were also evaluated by examining the slope and intercept of predicted vs. observed BW for the 21 sets of lactation data. Using PROC MEANS in SAS (SAS Institute, 2000), the slopes and intercepts were tested for significant difference from the line of unity (slope of one, intercept of zero). To identify patterns of bias in the predictions, residual BW (predicted observed) was regressed against WOL and predicted BW, and the average linear and quadratic coefficients of the plots were tested against zero using PROC MEANS in SAS. Mean predicted BW over the full lactation was compared against the observed means by t-test, where n = 21.
To account for the bias detected in predicted energy balance, maintenance energy expenditures, in Mcal/kg of BW 0.75, were estimated for each of the 21 data sets as the value that minimized residual sum of squares (RSS) between predicted and observed BW from wk 4 through 37 of lactation. Weekly maintenance energy expenditures were also calculated from the weekly DMI, FCM yield, and BW observations as the value that yielded a NEL balance of zero; that is, predicted BWn = observed BWn.
| RESULTS |
|---|
|
|
|---|
|
New Maintenance Descriptions
In an attempt to correct the bias in BW prediction, an optimum estimate of maintenance energy expenditures, as a function of BW0.75, was obtained for each of the 21 data sets. The average expenditure was 0.097 ± 0.004 Mcal of NEL/kg of BW0.75 (Table 2
). An adjusted mean, eliminating 2 values that were likely sources of errors (data sets 10 and 19), was 0.096 ± 0.003 Mcal of NEL/kg of BW0.75, a 1% change that reduced SE by 25%. In addition, second-order polynomial equations relating weekly maintenance estimates to WOL were parameterized for each of the 21 data sets (Table 2
), and for the average of weekly estimates obtained (Table 2
, Figure 2
). The average maintenance energy expenditure at the onset of lactation was approximately 0.08 Mcal of NEL/kg of BW0.75. This value increased until wk 15 of lactation where it plateaued at approximately 0.098 Mcal of NEL/kg of BW0.75 (Figure 2
).
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Weighting time-independent and time-dependent maintenance equations by the number of cows per treatment result in a time-independent maintenance equation of 0.94 Mcal of NEL/kg of BW0.75, and a time-dependent equation of: Maintenance NEL = [0.00002014 x WOL2 + 0.00122316 x WOL + 0.07643808] Mcal of NEL/kg of BW0.75, in which the same 2 data sets were removed as in the original scenario. Thus, weighting had little effect on the parameter estimates. Although weighting can remove bias in parameter estimates introduced by data sets with a low n, the goal of our subsequent work (Ellis et al., 2006) was to attempt to correct bias through a biological feedback of energy stores onto DMI predictions instead of a statistical weighting procedure. Using weighted maintenance parameters would result in maintenance expenditures further away from what some individual data sets need to correct the energy balance problem, and this additional error would have to be absorbed by the feedback mechanism. For this reason, although we present the weighted factors, the unweighted factors were focused on for model development.
Potential Sources of Error
Components of the energy balance model other than the maintenance calculation could have contributed to the error in BW prediction. Milk production, expressed as kilograms of 4% FCM (Gaines, 1928), allows a constant energy content of 0.749 Mcal of NEL/kg to be used (NRC, 2001), except when milk fat is less than 3% (NRC, 2001, such as during milk fat depression due to dietary manipulation (Weiss, 2002). In the performance data used for model evaluation here, none of the average milk fat percentages fell below 3%, and milk fat depression was not evident or reported. It is unlikely, then, that using 4% FCM to describe milk production substantially contributed to the error in energy balance seen here.
A second potential source of error is calculation of the energy content of the feed. The NRC (2001) computer model program is commonly used to calculate dietary NEL in North America, and is an improvement over the NRC (1988) method of determination. In 2 separate evaluations of accuracy, NEL values for feeds according to NRC (1989) were 5% (Weiss, 1998) and 5 to 7% (Vermorel and Coulon, 1998) too high, whereas with the NRC (2001) method, overprediction of dietary NEL was decreased to 1.2% (Weiss, 1998). Depending on the feed used, prediction of dietary energy content with this methodology may partially contribute to poor prediction of BW within an energy balance model. However, the 1.2% overprediction reported by Weiss (1998) is not large enough to account for the degree of BW overprediction evident in this study. Adjustment downwards of dietary NEL within the model by 1.2% decreased the average predicted BW from 602.7 to 593.9 kg and decreased the total RSS by approximately 15%. Adjusting dietary NEL accounted for only 8.8 kg (20%) of the average 44.5-kg BW overprediction, leaving 80% of the BW overprediction unaccounted for. Although calculation of dietary NEL using the NRC (2001) methodology might have contributed to error in predicted BW within the energy model, it was not the major contributing factor. However, until systems that can more accurately represent the true dietary NEL are available, these limitations must be acknowledged and accepted as likely having a role in contributing to errors in energy balance calculations.
In this study, the energetic cost of pregnancy was not included in the estimation of energy balance. There is usually little or no increase in energy requirements due to pregnancy during lactation, and thus any increases in energy expenditures due to pregnancy would not come into play until after the cessation of lactation (Bauman and Currie, 1980). Because none of the data sets used here extended beyond 37 WOL, excluding pregnancy did not contribute to overprediction of BW. In addition, calf weight should not significantly affect average daily gain of the cow until days pregnant is greater than 190, or approximately 40 WOL (NRC, 2001). Thus, the actual weight of the fetus would not affect observed BW against which the predicted values, which do not consider fetus weight, were tested.
Another source of error to be considered is that changes in BW may not always reflect changes in retained tissue energy (NRC, 1988). From wk 5 to 12, energy reserves can differ by 40% with no change in BW (NRC, 1988), due to rapid changes in DMI and gut fill. Chilliard et al. (1991) found that a 4-kg increase in gut fill accompanies a 1-kg increase in DMI, whereas more recent studies (Komaragiri and Erdman, 1997; Komaragiri et al., 1998) suggest a 2.5-kg increase per 1-kg increase in DMI. Because tissue mobilization occurs at the same time that DMI is rapidly increasing in early lactation, decreases in BW due to depleting energy stores are masked by a simultaneous increase in gut fill. After peak lactation, DMI and gut fill decrease at the same time as body energy stores begin to increase. Although BW may not always accurately represent changes in body tissue weight, the model predicts similar BW at the end of lactation whether the effect of gut fill is included or not and including this effect does not correct the substantial overprediction of BW. If observed BW is adjusted downwards to account for a 2.5-kg increase in gut fill per 1-kg increase in DMI, the overall mean BW changes by 0.03%.
The NEL content of BW gain/loss is another potential source of error. The NRC (1988) values were used here because information on BCS was not available, and is required for the NRC (2001) system. The energy content of gain/loss definitely changes across BCS, and the error in using average equations to describe all cows is acknowledged. However, analysis here showed that this model is fairly insensitive to these values. The energy content of gain would have to be increased by over 200% to correct the discrepancy between predicted and observed BW. Although these factors contribute to overall error, they are likely not the major contributing factor.
Additional maintenance energy expenditures for growing animals were added to the energy balance model for first-lactation data sets only. Unfortunately, not all data sets reported age of the animals used, and all these studies grouped animals in any lactation greater than one together. It is acknowledged that there will be some growth in cows in their second lactation but, determination of which, if any, data sets should be given additional energy expenditures was impossible here. This error is acknowledged, and although it may contribute to some of the error in energy balance calculation, it will remain unaddressed until more descriptive data come available.
The above review suggests that, although errors may exist within components of energy balance calculation other than the maintenance cost, none by itself is substantial enough to account for the degree of BW overprediction reported here. Furthermore, even when adjustments were applied simultaneously to the energy content of feed, gut fill, growth, and the NEL content of BW gain and loss, the maintenance energy expenditure equation still has to be increased to approximately 0.090 Mcal/kg of BW0.75 to correct the remaining BW overprediction. The body of literature suggesting that maintenance energy expenditures during lactation are currently underestimated (for examples, see Yan et al., 1997; Kirkland and Gordon, 2001; Kebreab et al., 2003) supports the idea that the BW overprediction in this model is caused at least partially by underestimation of maintenance energy expenditures.
Correction of the Bias
Maintenance energy expenditure by the lactating dairy cow, defined here as the NEL equivalent used to maintain the animal, as it is, at zero BW content change, is described by the NRC (2001) as 0.08 Mcal/kg of BW0.75, based on measured fasting heat production estimates of 0.073 Mcal/kg of BW0.75 for dry nonpregnant dairy cows housed in tie stalls or metabolic chambers (Flatt et al., 1965). A 10% allowance for normal voluntary activity of cows that would be housed in dry lot or freestall systems increased the factor to 0.08 Mcal/kg of BW0.75 (NRC, 2001). The dairy estimates agree with beef cattle estimates (NRC, 1996), in which the base maintenance requirement, 0.065 Mcal/kg of BW0.75, with a Holstein/Jersey breed adjustment factor of 1.2, results in a maintenance energy requirement of 0.079 Mcal/kg of BW0.75.
Several studies have shown that energy expenditures for maintenance are substantially higher than this during lactation. By nonlinear and linear regression analysis of milk energy output on energy intake adjusted to zero energy balance, Kebreab et al. (2003) estimated maintenance energy expenditures to be 16 to 22% higher during lactation than predicted by NRC (2001). Similarly, Moe et al. (1970) reported a 22% increase in maintenance energy expenditures based on multiple regression analysis of energy balance data. Although such estimates were based on extrapolation outside the range of observations to zero energy retention, indirect calorimetry has also yielded estimates of the lactation effect on maintenance energy expenditures as an increase of 49% (Ritzman and Benedict, 1938), 23% (Flatt et al., 1969), 27% (Kirkland and Gordon, 2001), and 29% (Yan, et al., 1997). According to observed changes in BW, maintenance energy expenditures have been estimated to be increased more than 30% (Neville and McCullough, 1969) or 10 to 20% (Ferrell and Jenkins, 1987) during lactation.
Smith and Baldwin (1974) showed that organ weights were higher in lactating compared with nonlactating cows for liver (25%), heart (22%), mammary gland (73%), lungs (22%), rumen (20%), abomasums (35%), intestines (31%), spleen (20%), and adrenals (19%). Assuming no change in energy expended/weight of tissue, a 10% increase in energy expenditures during lactation could be attributed to differences in organ tissue weight alone. According to our estimates, maintenance energy expenditures rose until approximately 15 WOL, where it plateaued at approximately 0.098 Mcal of NEL/kg of BW0.75 (Figure 2
). The 0.096 time-independent equation represents the weighted average of this curve. Interestingly, the 0.08 Mcal of NEL/kg of BW0.75 maintenance equation (NRC, 2001) seems to agree with values obtained here in early lactation (Figure 2
). It is only as the demands of lactation increase and lactation persists that maintenance energy expenditures rise to a higher value. It appears that the 0.08 maintenance equation still describes well the energy expenditures of a nonlactating dairy cow. Maintenance energy expenditures rise in early lactation because milk production and DMI both peak within the first 15 WOL, placing increased demands on the internal organs to increase in both size and activity. In comparative slaughter measurements, rumen, small intestine, and liver weights of dairy cattle increased as a percentage of empty BW from 2 to 17 WOL, and then decreased by wk 34 (Baldwin et al., 2004). Similarly, weights of ewe liver and rumen reached a maximum at 6 WOL, and small intestines at 4 WOL, and weights declined thereafter (Fell et al., 1972). The increase, plateau, and eventual decline in metabolically active organ size and activity are likely associated with a similar pattern in maintenance energy expenditures. Although our data showed no decline in maintenance energy expenditures, it is possible, had the data extended further into late lactation, that a decline would be observed. Some individual data sets did show a decline in maintenance energy expenditures after peak lactation, but averaging curves smoothes out any rise to peak, peak, or eventual decline observed. More extensive data would be required to determine the exact pattern of maintenance energy expenditure change over the course of lactation.
Increased metabolic activity of organs to support lactation could be considered either as a maintenance cost or an inefficiency of conversion of ME to NEL. The 2 are actually indistinguishable mathematically and the selection of one over the other remains a matter primarily of philosophical preference. In the NRC (2001) system, an efficiency adjustment to correct the bias we encountered would be accommodated by changing NEL values of feedstuffs, which is rather more difficult to do than the alternative maintenance adjustment undertaken here.
The variation in energy balance bias represents a problem in the use of one equation to describe all lactating dairy cows. The between-data set variation in maintenance energy expenditures around the updated 0.096 maintenance equation was 25%. The highest maintenance value was 0.135 and the lowest was 0.068 Mcal of NEL/kg of BW0.75. Although this variation likely encompasses multiple sources of error, others have reported similar, if not higher, degrees of variation. Ritzman and Benedict (1938) found between-animal variation in maintenance energy requirements under uniform conditions of up to 50% in Holstein and up to 71% in Jersey cows. Carstens et al. (1989) and Taylor and Young (1968) found that maintenance requirements varied by 20 to 30% due to genetic differences alone, which appeared moderately to highly heritable.
Genetic potential for production, previous plane of nutrition, estimation of dietary energy content, environment, BCS, and physiological status likely contribute cumulatively to a high degree of variation between animals. It was evident from a review of the literature that maintenance is not a constant function of BW (Webster et al., 1982; Ferrell and Jenkins, 1985). The possible roles of several factors contributing to the variation in energy balance bias reported here were examined. The relationship between maintenance energy expenditures and average BW, average DMI, peak FCM yield, total FCM yield, and dietary NEL content were not significant. In contrast, Ferrell and Jenkins (1987) found that maintenance energy expenditures increased about 9.6 kcal/kg of BW0.75 per d for each 1-kg increase in milk yield at peak lactation. Interestingly, they reported that 40% of the variation in maintenance energy expenditures could be explained by variation in milk production potential alone.
It is important to reiterate that the maintenance description developed here is essentially an error term, containing errors associated with the other aspects of energy balance calculations discussed above. It was not the goal of this paper to assign individual correction factors to all of these components, but to provide a feasible solution to the model, agreeable with data in the literature, such that the model could be used for further evaluations of DMI prediction equations.
| CONCLUSIONS |
|---|
|
|
|---|
The 0.096 Mcal of NEL/kg of BW0.75 equation is 20% higher than the original 0.08 Mcal of NEL/kg of BW0.75. This increase is in agreement with estimates of maintenance energy expenditures presented in the reviewed literature, which ranged from 10 to 49%, but averaged from 20 to 30% (for examples see Moe et al., 1970; Yan et al., 1997; Kirkland and Gordon, 2001). This new time-independent equation, as well as a time-dependent quadratic equation relating maintenance expenditures to WOL, are useful in better describing energy balance within a dynamic model, but do little to account for the substantial variation that is evident. It is the goal of the companion paper (Ellis et al., 2006) to address between-data set variation and correction of errors that persist in accumulating over time.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Received for publication October 2, 2005. Accepted for publication January 12, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. P. Weiss, D. J. Wyatt, and T. R. McKelvey Effect of Feeding Propionibacteria on Milk Production by Early Lactation Dairy Cows J Dairy Sci, February 1, 2008; 91(2): 646 - 652. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. C. Friggens, C. Ridder, and P. Lovendahl On the Use of Milk Composition Measures to Predict the Energy Balance of Dairy Cows J Dairy Sci, December 1, 2007; 90(12): 5453 - 5467. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. C. Friggens, P. Berg, P. Theilgaard, I. R. Korsgaard, K. L. Ingvartsen, P. Lovendahl, and J. Jensen Breed and Parity Effects on Energy Balance Profiles Through Lactation: Evidence of Genetically Driven Body Energy Change J Dairy Sci, November 1, 2007; 90(11): 5291 - 5305. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. Ellis, F. Qiao, and J. P. Cant Prediction of Dry Matter Intake Throughout Lactation in a Dynamic Model of Dairy Cow Performance J Dairy Sci, May 1, 2006; 89(5): 1558 - 1570. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |