|
|
||||||||
USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE 68933
1 Correspondence: P.O. Box 166 (phone: 402-762-4248; fax: 402-762-4209; E-mail:
williams{at}email.marc.usda.gov).
| Abstract |
|---|
|
|
|---|
Key Words: Beef Cattle Energy Metabolism Models Weight Gain
| Introduction |
|---|
|
|
|---|
| Materials and Methods |
|---|
|
|
|---|
Table 1
contains a list of acronyms used in this paper. The EBW of mature cattle that contains 25% ether-extractable lipid is referred to as the standard reference EBW (SREBW). Williams and Jenkins (1997) published estimates of SREBW for several breeds and breed crosses of beef cattle. Major components of EBW are ether-extractable lipid, referred to as empty body fat weight (FAT), and empty body fat-free weight (FFM). Only the protein fraction in FFM (PRO) contains energy; hence, ME is retained in FAT and PRO. Daily rates of change in these five components (FBW, EBW, FAT, FFM, PRO) are prefixed with the letter d (e.g., daily rate of change in EBW is dEBW). These daily rates of change are part of a system of differential equations that define the models in this paper and a solution for this system of equations is obtained using a fourth-order Runge-Kutta procedure (Shampine and Watts, 1970) for numerical integration.
|
Current systems that predict dEBW first predict recovered energy (RE), and then use the predicted values of RE to predict dEBW. The partitioning of MEI in these systems was discussed by Williams and Jenkins (2003) and is summarized with the following equation:
![]() | [1] |
where MEm is the amount of MEI used for maintenance, kg is the net efficiency of ME utilization for gain, and the term (MEI - MEm) is MEg. Two of the most widely used models to predict dEBW are from ARC (1980) and NRC (2000). Starting with a known MEI, these models predict dEBW by first estimating MEm, and then predicting kg from the dietary ME density. These estimates of MEm and kg are used as inputs to Eq. [1]
to estimate RE, and empirical equations are then used to predict dEBW directly from RE.
In the ARC (1980) system, RE is predicted with the following equation:
![]() |
where Q is the metabolizability of the diet, I is feed intake, FHP is fasting heat production, and km is the efficiency of ME utilization for maintenance. The term Q x I is MEI, and FHP/km is MEm, which makes this equation the same as Eq. [1]
. The estimate of RE is used to predict dEBW with the following equation:
![]() |
where EVG is the energy value of dEBW and is predicted in the ARC (1980) system with the following equation:
![]() |
The NRC (2000) system predicts RE with the following equation:
![]() |
where NERm is the NE requirement for maintenance, which is the same as FHP, and NEg and NEm are the dietary concentrations of NE for gain and maintenance, respectively. Using relationships (NRC, 2000) based on the concentration of ME in diet (MD), NEg and NEm are equivalent to kg x MD and km x MD, respectively, and we can rewrite the above equation as follows:
![]() |
This equation for RE is the same as Eq. [1]
, which makes the formulations in ARC (1980) and NRC (2000) the same. The systems differ in the assumptions used to estimate partial efficiencies of ME utilization for maintenance and gain and the requirements for maintenance.
To predict dEBW, NRC (2000) proposed the following allometric relationship between RE and dEBW:
![]() | [2] |
Estimates of a and b were obtained by linear regression using the log transformation of Eq. [2]
, thus:
![]() |
Estimates of ln(a) and b were -2.7566 and 1.097, respectively, for medium-frame steers. The reciprocal relationship in the above equation is used to predict dEBW as follows:
![]() |
Models developed by Notter (1977) and Oltjen et al. (1986) use a more mechanistic approach to predict dEBW. These models use the same method as in Eq. [1]
to obtain an estimate of RE, and the approach to predict dEBW is based on using an underlying protein growth curve to estimate the potential dPRO, and dFAT is calculated as a residual as shown in the following equations:
![]() |
![]() |
![]() |
![]() |
where FP is the fraction of dFFM that is protein, and 5.7 and 9.5 represent the amount of energy in Mcal/kg of DM of protein and fat, respectively (Brouwer, 1965).
Proposed Model to Predict dEBW Using MEg as the Input
Components of a model that uses an input value for MEg to predict dEBW are illustrated in Figure 1
. The first component is the prediction of RE from MEg, and this is accomplished through the prediction of kg, which is used to predict RE from MEg. The second component is the prediction of dEBW using RE as the input. This is accomplished by first partitioning RE into energy retained in protein and fat, and then estimating the weights of protein (dPRO) and fat (dFAT) that are retained. The fraction of protein (FP) in dFFM is then estimated and used to calculate dFFM from dPRO. Finally, dEBW is calculated as the sum of dFFM and dFAT. The approach used to predict dEBW based on these two component models is discussed below.
|
![]() | [3] |
The term HiEv is a component of heat production that is attributable to support metabolism. The system used by NRC (2000) assumes that kg is constant for a specific diet, and the system used by ARC (1980) assumes that kg is constant for a specific diet fed at a specific level of intake; hence, the relative contributions of protein and fat to RE have no impact on kg (i.e., composition of gain has no impact on kg). Conversely, experimental results with sheep and cattle show that kg varies with composition of gain.
Graham (1980) found that kg increased from 0.32 in weaned lambs at 2 mo to 0.55 for a similar diet given to the same sheep at age 10 mo, and that subsequently, it did not vary for sheep up to 6 yr old. Blaxter et al. (1966) found no significant difference in kg in cattle aged 15 to 81 wk. Average daily gain in these cattle contained 30.5 kcal of protein energy/100 kcal of retained energy at 15 wk of age, and 24.8 kcal of protein energy/100 kcal of retained energy at 81 wk of age. These results show that there was very little difference in composition of gain in the data, which may explain why differences in kg were not significant. Other evidence reviewed by Vermorel and Bickel (1980) and studies with Hereford and Holstein steers (Garrett, 1971) indicate that kg increased with the proportion of fat in gain. Metabolizable energy is retained in protein and fat, and as an animal grows from birth to maturity without any restriction in feed availability, the proportion of fat in dEBW increases. Hence, the above experimental results were interpreted to suggest that ME is retained in fat with a higher efficiency than in protein. This suggests that kg should be weighted by the relative contributions of protein and fat to RE, using separate partial net efficiencies for ME retention in protein (kp) and fat (kf).
Pullar and Webster (1977) measured energy balance during growth in lean and fatty zucker rats because the phenotypes differ markedly in the way they partition ME to energy retained in protein and fat. These authors obtained values of 0.45 and 0.75 for kp and kf, respectively, by multiple regression analysis. Accepting a kp value of 0.45 for ruminants limits kg to a minimal value of 0.45; however, this minimum value would most likely be greater than 0.45 since in growing ruminants, energy retained in protein is never 100% of RE. This minimal kg would be much greater than that in the results of Graham (1980) for 2-mo-old weaned lambs. Also, the equation published by NRC (2000) to calculate kg gives a value of 0.296 at a ME concentration of 2.0 Mcal/kg of DM. This suggests that in ruminants, energy is retained in protein much less efficiently than in nonruminants.
Ørskov and McDonald (1970), Bickel and Durrer (1974), and Rattray and Joyce (1976), in experiments with ruminants, obtained estimates for kp and kf that varied from 0.56 to 1.28 and from 0.18 to 0.36, respectively. Geay (1984) reviewed 52 sets of published data in ruminants on energy retention as protein and found that kg decreased curvilinearly as the proportion of protein energy in RE increased. Analysis of these data by Geay (1984) resulted in estimates of 0.2 and 0.75 for kp and kf, respectively. These estimates agree fairly well with the data collected in sheep by Ørskov and McDonald (1970), Rattray and Joyce (1976), and Notter et al. (1984). These results suggest that in ruminants, energy is retained in fat with about the same efficiency as that in nonruminants, but energy retention in protein in ruminants is less than half as efficient as that in nonruminants. Feed to gain ratio of cattle on high-energy diets in feedlots is about 6.5:1, and in chickens, it is less than 2:1. This supports the lower efficiency of energy retention as protein in ruminants. Partial net efficiencies of ME utilization of 0.2 for protein gain and 0.75 for fat gain will be used to derive an overall equation to predict kg as a function of the proportion of RE that is in protein.
Metabolizable energy that is available for gain is recovered in protein and fat, and using the partial net efficiencies for ME recovery as protein energy (REp) and fat energy (REf), RE can be expressed as follows (Geay, 1984):
![]() | [4] |
Solving for kg in terms of REp and RE (REf = RE - REp) we get:
![]() | [5] |
This estimate of kg is used in Eq. [3]
to calculate RE. With this equation, if REp = 0, then kg = 0.75, and if REp = RE, then kg = 0.2. The relationship in Eq. [5]
between kg and REp/RE is plotted in Figure 2
, which shows the curvilinear decrease in kg from 0.75 to 0.2, as REp/RE increases from 0 to 1. The unknowns in Eq. [5]
are REp and RE, and methods to estimate them are discussed in the next section.
|
The model of Williams and Jenkins (1998) was used to predict RE as a function of dEBW. This model uses dEBW as its input to predict the amount of dFFM in dEBW, and dFAT is obtained as a residual according to the following two equations:
![]() | [6] |
![]() | [7] |
The parameters b and c model the impact of breed and previous and present level of nutrition on the composition of dEBW. A complete description of the development and parameterization of the model can be found in Keele et al. (1992), Williams et al. (1995), and Williams and Jenkins (1997; 1998).
The components of dEBW in Eq. [6] and [7]![]()
can be used to estimate RE if the FP in dFFM can be calculated. Data from Buckley (1985) was used to develop an empirical equation to estimate FP. These data contained observed values on the chemical empty body composition of Hereford, Charolais, and Simmental heifers slaughtered at birth, 3, 7, 8, 10, and 14 mo of age, which showed that FP in FFM varied from 20.5% at birth to 24.5% at 14 mo of age. Values of FP in these data were regressed on linear and quadratic terms for UEBW (EBW/SREBW) to obtain the following equation:
![]() | [8] |
and dPRO was predicted as protein weight on day t minus the protein weight on day t - 1, according to the following equation:
![]() | [9] |
With this estimate of dPRO, RE was estimated as follows:
![]() | [10] |
where 5.7 x dPRO is REp and 9.5 x dFAT is REf.
This above system of equations (Eq. [6], [7], [8], [9], [10]![]()
![]()
![]()
![]()
) shows the prediction of REp, REf, and RE using dEBW as the input. The predicted values of REp and RE provide the input to predict kg with Eq. [5]
. The final step in the development of this component model is the prediction of dEBW using RE as the input, and the approach in this case will be based on reformulating Eq. [6]
through [10]
to derive reciprocal relationships that can be used to predict dEBW with RE as the input.
Recovered energy in dFFM is in the form of protein; hence, the energy density of dFFM (FFME) can be expressed as follows:
![]() | [11] |
and
![]() |
We can now rewrite Eq. [10]
as follows:
![]() | [12] |
Next, we substitute Eq. [6] and [7]![]()
for dFFM and dFAT in Eq. [12]
to derive the following:
![]() |
![]() |
Then, collecting terms for dEBW we obtain:
![]() |
Solving for dEBW, we arrive at the following:
![]() | [13] |
This equation expresses dEBW as a function of RE and other variables, and it completes the development of the component model to predict dEBW using RE as the input.
Integrating the Component Models.
The primary objective is to predict dEBW using MEg as the input. This is achieved by first predicting RE from MEg, and then predicting dEBW from RE. The processes contained in these two component models were integrated into a single model, which is illustrated in Figure 3
. The model requires initial values for EBW, PRO, FAT, and the previous level of gain. These initial values were used to calculate starting values for FFME, c, b, FP, and FFM. At the start of the simulation, an initial value for kg was calculated from the ME density of the diet according to NRC (2000). With a fixed input value for MEg, RE is determined according to Eq. [3]
. This predicted RE value, along with initial values for FFME, b, and c are used in Eq. [13]
to predict dEBW. The predicted dEBW is partitioned into dFFM and dFAT using Eq. [6] and [7]![]()
, respectively. Both dEBW and dFFM are used to update the integrals EBW and FFM, respectively. The value for EBW is used in Eq. [8]
to calculate a new value for FP, and this is used with the updated value for FFM in Eq. [9]
to calculate dPRO. Values for dPRO and dFFM are used in Eq. [11]
to calculate a new value for FFME. Values for dFAT and dPRO are multiplied by 9.5 and 5.7, respectively, to obtain estimates for REp and REf, and these are summed to obtain RE as shown in Eq. [10]
.
|
The input value of MEg would be negative when MEI is less than MEm. In this case, the model allows for the mobilization of body tissues to satisfy the negative MEg, and RE is predicted with the following equation:
![]() |
We assumed that when body tissues are mobilized, only 80% of the energetic value of the tissues would be available, and km is used instead of kg since the energy is being used to satisfy maintenance requirements. The value for RE in this case is negative, and when used in Eq. [13]
, the predicted value for dEBW would be negative, resulting in EBW loss.
On the beginning day of the simulation, initial conditions are set and usually step size has to be halved before the value of kg converges. On the next day, starting conditions are more accurate and the kg value usually converges without halving the step size. This part of the model is based on a very strong relationship between dEBW and composition of dEBW that responds to both previous and current levels of nutrition and breed composition. In other applications, this model, as illustrated in Figure 3
, can also be used in the reverse manner to predict MEg using dEBW as the input, thus providing a means to predict a part of the feed requirements to achieve a certain level of response.
| Results and Discussion |
|---|
|
|
|---|
These results show that kg decreased as the energy density of the diet decreased using constant values for kp and kf. The value of kg is calculated in the model as a function of composition of gain and not energy density of the diet; hence, the lower kg value obtained with the L diet compared with the H diet is a result of leaner EBW gains with the L diet. This is further illustrated with the lower kg value obtained with the L75 diet compared with the L diet. In this case, both diets had the same energy density, but steers on the L75 diet had leaner gains compared with steers on the L diet.
The model illustrated in Figure 3
was tested with experimental data (Smith et al., 1976; Cundiff et al., 1981, 1984) on MEI during the finishing period of Hereford x Angus steers in the first three cycles of the Germplasm Evaluation Program at the U.S. Meat Animal Research Center. Daily ME intake was partitioned to maintenance, support metabolism, and MEg according to the model developed by Williams and Jenkins (2003), and MEg was used to simulate growth and body composition of these steers with the model shown in Figure 3
. Energy retained in dFAT and dPRO was predicted on a daily basis and used to calculate kg according to Eq. [5]
and also to calculate dEBW using a predicted value for FP according to Eq. [8]
. Results of this test showed a predicted value of 233 kg for gain in EBW over the feeding period compared with an observed gain of 234 kg, and this is evidence that the model, as illustrated in Figure 3
, was formulated and programmed correctly.
Predicted values for kg in Figure 4
increased from 0.39 to 0.56 as UEBW increased from 0.34 to 0.7. The rate of increase in kg decreased up to a value of about 0.55 for UEBW and was almost constant between UEBW values of 0.5 and 0.7. This initial decrease in the rate of increase in kg is a result of the impact of previous nutrition on the composition of present gains in EBW. In these experiments, steers were conditioned at a low ADG for an average duration of 35 d before starting the finishing phase, and this resulted in high protein gains at the start of the finishing phase. A distributed lag function is used to model the transition from one nutritional level to another (Keele et al., 1992). At the start of the simulation, this results in a larger than normal decrease in the fraction of protein and a larger than normal increase in the fraction of fat in dEBW, and this causes a large increase in the value of kg. As days in the finishing phase increase, both the rate of decrease in the fraction of protein and rate of increase in the fraction of fat in dEBW decrease and gradually approach the values for the present plane of nutrition, and the impact of previous nutritional level decreases and eventually disappears.
|
| Implications |
|---|
|
|
|---|
Received for publication September 19, 2002. Accepted for publication February 20, 2003.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
F. C. P. Castro Bulle, P. V. Paulino, A. C. Sanches, and R. D. Sainz Growth, carcass quality, and protein and energy metabolism in beef cattle with different growth potentials and residual feed intakes J Anim Sci, April 1, 2007; 85(4): 928 - 936. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. H. M. R. Fernandes, K. T. Resende, L. O. Tedeschi, J. S. Fernandes Jr., H. M. Silva, G. E. Carstens, T. T. Berchielli, I. A. M. A. Teixeira, and L. Akinaga Energy and protein requirements for maintenance and growth of Boer crossbred kids J Anim Sci, April 1, 2007; 85(4): 1014 - 1023. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. B. Williams, G. L. Bennett, T. G. Jenkins, L. V. Cundiff, and C. L. Ferrell Using simulation models to predict feed intake: Phenotypic and genetic relationships between observed and predicted values in cattle J Anim Sci, June 1, 2006; 84(6): 1310 - 1316. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. B. Williams Technical Note: A dynamic model to predict the composition of fat-free matter gains in cattle J Anim Sci, June 1, 2005; 83(6): 1262 - 1266. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. B. Williams and T. G. Jenkins A dynamic model of metabolizable energy utilization in growing and mature cattle. I. Metabolizable energy utilization for maintenance and support metabolism J Anim Sci, June 1, 2003; 81(6): 1371 - 1381. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. B. Williams and T. G. Jenkins A dynamic model of metabolizable energy utilization in growing and mature cattle. III. Model evaluation J Anim Sci, June 1, 2003; 81(6): 1390 - 1398. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |