|
|
||||||||
1 Danish Institute of Agricultural Sciences, Department of Animal Health and Welfare, Research Center Foulum, PO Box 50, DK-8830 Tjele, Denmark
2 Animal Biology Division, Scottish Agricultural College, West Mains Road, Edinburgh EH9 3JG, UK
| ABSTRACT |
|---|
|
|
|---|
Key Words: lactation adipose body condition score genotype
Abbreviation key: ME = metabolizable energy
| INTRODUCTION |
|---|
|
|
|---|
The view taken here is that both the mobilization of body reserves in early lactation, and the subsequent gain in body reserves during pregnancy, are to a large extent genetically driven. Seen in this way, body energy mobilization is not a response but rather a natural component of the reproductive cycle. It would be expected to be independent of food energy content. As shown in Figure 1
, ignoring this preprogrammed body mobilization has important consequences for the prediction of energy requirements and the intake necessary to meet these requirements. With very few exceptions, current methods to predict energy requirements are based on estimates of milk production and maintenance. They do not explicitly allow for any genetically driven body energy mobilization. Prediction of the cows energy requirements can be substantially improved, particularly in early lactation, by incorporating genetically driven body energy mobilization.
|
| RATIONALE |
|---|
|
|
|---|
The existence of environmentally driven body lipid mobilization, where mobilization is a response, is not in doubt. There are numerous examples showing that reduction in nutrient availability, caused by a decrease in the supply or quality of the food being offered, results in an increase in mobilization of body lipid (e.g., Friggens et al., 1998). The ability to buffer the consequences of a constraining nutritional environment is an important function of the body reserves. However, the existence of environmentally driven mobilization of the body reserves does not preclude the existence also of a genetically driven body mobilization. The existence of genetically driven body lipid mobilization is less well accepted, especially in relation to the postpartum period. It can be inferred from consideration of the evolutionary benefits of safeguarding reproductive success by strategic use of body reserves (Pond, 1984; Knight, 2001; Friggens, 2003).
The genetically driven accumulation of body reserves by the mother during pregnancy to support the energetic demands of the following lactation confers obvious advantages (Vernon and Pond, 1997; Oftedal, 2000). However, it is not so obvious that there should be a genetically driven loss of body lipid in early lactation. If the environmental pressure to mobilize body lipid in early lactation were removed, for instance by provision of super-abundant nutrition, why should the mother still mobilize? There are evolutionary arguments as to why it would be advantageous to decrease the size of the body reserves with increasing time from birth (Friggens, 2003) that are related to the cost of having body reserves. Carrying additional body mass increases the energy costs of movement. The reduction in mobility that results will, in a natural environment, increase the risk of predation in predated species and the risk of failing to predate successfully in predators.
The extent to which these arguments are of major relevance to domestic livestock such as dairy cattle is open to question. Experimental studies may offer more direct evidence, particularly those that have sought to modulate body reserves and rates of mobilization. The evolutionary arguments alluded to above imply that at any given time in pregnancy and lactation there is an optimal level of body fatness that the animal is genetically driven to attain. There is thus a trajectory or pattern of body fatness through the reproductive cycle. The notion of genetically determined trajectories in development of body tissues, including fat, has long been accepted in the study of growth (Hammond, 1940; Waddington, 1957; Emmans, 1988). An important feature of such trajectories is that they are defended, and that deviations from the trajectory are generally followed by compensatory or corrective adjustments to regain the trajectory once environmental conditions allow (Allden, 1970; Wright and Russel, 1991).
It is commonly found that dairy cows become increasingly fat as they progress through pregnancy (e.g., Koenen et al., 2001). However, at the same time milk yield is usually declining and it is usually only in the last 2 mo of pregnancy that the cow is not concurrently lactating. This makes it difficult to attribute the change in body fatness solely to the degree of pregnancy, and a definitive study in dairy cattle is lacking (Knight, 2001). There is, however, some evidence that the trajectory of body fatness through pregnancy, usually measured as body condition score, is defended. A striking feature of multiple-lactation studies that involve underfeeding is that, although underfeeding causes substantially greater loss of body reserves in early lactation, there is generally little difference in condition score between underfed and control cows by the end of the subsequent pregnancy (Broster et al., 1993; Chilliard et al., 2000). When given favorable feeding, cows compensate for prior nutritional insult by rapidly regaining body condition (Friggens et al., 1998; Morrison et al., 1999). The same was found to apply to cows that had undergone a period of draught work (Zerbini et al., 1996).
There is also some evidence of cows defending the trajectory relative to overfeeding of energy. We have previously fed lactating cows and nonlactating heifers, with the same starting condition score, either a high or a normal energy density diet during the last 24 wk of gestation (Ingvartsen et al., 1995). The cows responded to the overfeeding by diverting almost all the increased energy into milk and appeared to defend their genetically driven level of body fatness. Only when the option of diverting energy into milk was not available, i.e., in the heifers, did the energy overfeeding distort the genetically determined trajectory of body fatness.
It thus appears that the drive to attain the genetically determined trajectory of body fatness through pregnancy is accorded a high priority. A possible mechanism is suggested by the available literature on the changes that occur in the metabolism and endocrine sensitivity of adipose tissue as pregnancy progresses. For most of pregnancy, the rate of lipogenesis is enhanced and the rate of lipolysis is relatively low, favoring lipid accretion (Vernon and Flint, 1984). As parturition approaches, these relative rates are reversed such that lipid mobilization is favored (Vernon and Flint, 1984; Vernon et al., 2001). A number of hormones have been implicated in the control of adipose metabolism in pregnancy, in particular insulin, and changes in insulin sensitivity of adipose tissue appear to be important (Vernon et al., 2001). More recently, a marked increase in expression of leptin mRNA by adipose tissue with advancing pregnancy has been reported in sheep (Ehrhardt et al., 2001).
The concerted changes in the endocrine milieu that start at the end of pregnancy continue in early lactation (McNamara, 1997). In relation to adipose tissue it is well documented that lipogenesis, especially reesterification of fatty acids, is heavily down-regulated and lipolysis is substantially increased in early lactation (McNamara, 1997; Chilliard et al., 2000). With respect to genetically driven trajectories of body mobilization, an important finding is that the responsiveness of adipose tissue to lipolytic stimuli (usually ß-adrenergic agonists) is markedly enhanced in early lactation (Metz and van den Bergh, 1977; Theilgaard et al., 2002). Not only is this responsiveness affected by energy balance and body fatness, but there is a marked effect of physiological state (lactating vs. nonlactating) on lipolytic response (Chilliard et al., 1998; Theilgaard et al., 2002). Manipulation of key hormones in early lactation such as insulin, growth hormone, and leptin also cause marked changes in body energy mobilization (Ingvartsen and Boisclair, 2001).
There is strong evidence that the trajectories of body fatness through early lactation are strongly defended. It has been repeatedly shown that nutritional manipulation of cows to be fatter, or leaner, than normal at calving provokes a change in subsequent body lipid mobilization such that normal levels of body fatness are regained approximately 3 to 4 mo after calving (Garnsworthy and Topps, 1982; Broster and Broster, 1998). Further, in situations where very high levels of nutritional resources are provided, mobilization of body lipid reserves still occurs, even though there is no apparent need for the energy this provides (Friggens et al., 1993). In dairy cows, attempts to decrease lipid mobilization by supplementation with dietary lipid in early lactation have not been successful (Gagliostro and Chilliard, 1991a; Grummer et al., 1995). We have found no convincing evidence that cows, which have gained body lipid in pregnancy, retain these extra reserves in the subsequent lactation.
In summary, the available evidence from the literature is generally consistent with there being a genetically driven component to changes in body reserves through both pregnancy and lactation. It is therefore not surprising that there are strong genetic correlations between measures of body fatness at different time-points in lactation (Coffey et al., 2001; Pryce et al., 2002). Genetically driven processes are, by their nature, amenable to description and therefore should be able to be incorporated into predictions of energy requirements. However, it is clear from the above considerations of trajectories that it is necessary to recognize 2 types of genetically driven body mobilization. As shown in Figure 2
, a cow that is currently in a nonconstraining environment can be either following her genetically predetermined trajectory of body fatness or be in the process of reconverging on that trajectory, i.e., compensating for a previous nutritional insult.
|
The method presented below for predicting genetically driven body lipid mobilization aims to describe the general biological framework for this process. Such a model would provide the basis for biologically meaningful genetic evaluation of genetically driven body lipid mobilization, though this is beyond the scope of the current paper.
| MODEL DESCRIPTION |
|---|
|
|
|---|
A1 The cow is driven to have a specific amount of lipid, here called L' (kg), at a particular time in lactation, here called T' (expressed as days from calving).
A2 The rate of change of body lipid, dL/dt (kg/d), itself changes linearly with time between calving and T'.
A3 If pregnant, the cow is driven to have a specific amount of lipid at the next calving, here called Lnext (kg).
A4 The rate of change of body lipid (dL/dt) itself changes linearly with time between conception and the next calving.
A5 At times greater than T', and if the cow is not pregnant then dL/dt is assumed to be 0, i.e., the cow has no drive to increase body lipid.
These assumptions, when combined with the appropriate inputs (described below), are sufficient to generate the curves shown in Figures 3
and 4
. The need for additional assumptions to deal with 2 particular cases becomes apparent when the following questions are considered.
|
|
What happens if conception occurs before T'?
If no limit is put on the rate of lipid loss then, because T' is fixed, cows that are relatively fat at calving can supply all of their energy needs from body reserves for some time. A consequence is that no food intake would be needed to meet the requirements of the cow. As resources other than energy will always be needed from the food this is an unreasonable result. To deal with this without expressly considering resources other than energy a further assumption is made.
A6 There is a maximal rate of lipid loss, called maxLipLoss (kg/d). This is a negative value of dL/dt.
When the initial rate of change of body lipid, dLcalv, is calculated to be less than maxLipLoss, the cow cannot, given A6, meet both assumptions A1 and A2 unless provision is made for allowing either L' or T' to change. Relaxing the assumption that the change in the rate of lipid change was linear was considered but rejected. The solution adopted requires the following further assumption.
A7 T' can increase to allow L' to be achieved without exceeding maxLipLoss.
The effect of this assumption is shown in Figure 5
. The broken line shows the case in which no maximal lipid loss is imposed. The solid line shows what happens when assumptions A6 and A7 are applied. The exact value of T' needed for dLcalv = maxLipLoss is obtained by replacing dLcalv with maxLipLoss in equation (1)
.
|
A8 When T' occurs after conception, the rate of body lipid change at T' (dLT') is the rate that the cow would ordinarily have had given her degree of pregnancy.
A9 A cow still losing body lipid after conception is driven to achieve not L' but the level of body lipid that she would ordinarily have had given her stage of pregnancy, i.e., L' plus the pregnancy gain in lipid. This modified target amount of lipid is called Ltarget.
Equations and Inputs
The model requires the following inputs: the amount of body lipid the cow has at calving (Lcalv) and the time of conception (Con). Lcalv can be derived from live weight and condition score at calving (see below). The model parameters are listed in Table 1
.
|
The case in which T' occurs before conception.
Based on assumptions A1 and A2, the initial rate of lipid mobilization, dLcalv, is calculated on the basis that the total amount of lipid mobilized, the area of the triangle defined by coordinates (calving, 0), (T', 0), and (calving, dLcalv) in Figure 3
, must be equal to L' - Lcalv. Therefore:
![]() | (1) |
If dLcalv is more negative than maxLipLoss (see A6) then the following adjustments are made according to A7:
![]() | (2) |
![]() | (3) |
At any time between calving and T' the rate of lipid change (dL/dt) is given by:
![]() | (4) |
When Lcalv = L' there is no change in lipid between calving and T'. When Lcalv is less than L' body lipid increases between calving and T'. Similarly, using assumptions A3 and A4, the rate of change of body lipid between conception and next calving can be calculated on the basis that the area of the triangle under the line after conception (Figure 3
) must be equal to Lnext - L'. Therefore, between conception and next calving:
![]() | (5) |
where dLnext is the rate of lipid change at the next calving and Gest is the length of pregnancy. At any time between conception and the next calving the rate of lipid change is given by:
![]() | (6) |
It follows that:
![]() | (7) |
The case in which T' occurs after conception.
In this case, L' is no longer the target level of lipid that the cow is driven to achieve. Instead, the cow is "aiming" for Ltarget (A9), and the rate of body lipid change for that stage of pregnancy (A8). In other words, dL/dt at T', here called dLT', is not zero (Figure 5
). Thus some adjustments to equation 1
are necessary, resulting in:
![]() | (8) |
Ltarget and dLT' can be derived from adjustments of the equations relating to the period from conception to next calving as follows:
![]() | (9) |
![]() | (10) |
Lipid gain during DFCon is the area of the triangle between the line from conception to T' and dL/dt = 0 in Figure 5
. Thus Ltarget is:
![]() | (11) |
where the asterisk denotes the product of the 2 associated parameters.
The value of dLT' is calculated from equation 7
with DFCon substituting for (DFCalv - Con), therefore:
![]() | (12) |
Substituting equation 12
in equation 11
gives:
![]() | (13) |
it follows that:
![]() | (14) |
The equation for dLT' (12)
contains only one unknown, T'. The same is true for Ltarget (equation 14
). Therefore, equation 8
appears to contain 2 unknownsdLcalv and T'. However, given assumptions A1, A2, A6, and A7, only one of dLcalv and T' can be unknown in any given situation. If T' is not changed from the default value, then dLcalv is unknown. If T' is unknown, i.e., it is changed from the default value, this is because dLcalv has been set to maxLipLoss and is thus known.
Equation 8
can now be solved for dLcalv, assuming T' is not changed from the default value. If the resulting value of dLcalv is less than maxLipLoss, then dLcalv should be set to maxLipLoss and equation 8
solved for T'. The relevant rates of body lipid change before and after T' are then given by equations 4
and 7
.
Estimation of Body Lipid Mass from Weight and Condition Score
Body lipid mass, for example Lcalv, can be calculated as the product of empty BW and the proportion of lipid in the empty body. These quantities can be estimated from live weight and BCS.
To convert BCS (measured on a 0 to 5 scale (Lowman et al., 1976)) to body lipid proportion, the following equation derived from the serial slaughter experiments of Wright and Russel (1984) on nonpregnant mature cows was used:
![]() | (15) |
Empty BW can be calculated from live weight if estimates of gut fill are available. In a given cow, it is not expected that gut fill increases with increasing body fatness. However, as condition score increases so does live weight. Thus, if gut fill is estimated as a function of live weight then that live weight should be standardized for condition score (e.g., Zygoyiannis et al., 1997). Standardized live weight and empty BW were derived by assuming that the live weight associated with a unit change in condition score was 0.129 (live weight at condition score 5) and that the gut fill was 0.15 of live weight at condition score 3 (Zygoyiannis et al., 1997).
Parameter Values and Model Evaluation
Initial values for the parameters in Table 2
were estimated from available studies, where it was judged that the trajectories of condition score change through lactation were largely unaffected by feeding and thus reflected genetically driven mobilization. Because a number of different condition score scales exist and because lipid mass at a given condition score is a function of the size of animal (Zygoyiannis et al., 1997), the values for L', and Lnext are given as proportions in Table 2
. Using the 0 to 5 condition score scale of Lowman et al. (1976), the value of 0.26 g of lipid/g of empty body for L' equates to condition score 2.5 and the value of 0.32 g lipid/g of empty body for Lnext equates to condition score 3.0. For the purpose of model evaluation, the results are not sensitive to the assumptions used for converting condition score and live weight to body lipid mass.
|
For the first test, back fat thickness measured by ultrasound scanning and live weight data from 36 Danish Holstein heifers in the period from calving to 22 wk postcalving were used. In the 24 wk up to calving, the heifers were split into 3 groups receiving rations of high, medium, and low quality designed to result in different levels of body fatness at calving. After calving, all the heifers received the same adequate TMR (11.7 MJ of metabolizable energy/kg of DM) ad libitum. In a previous experiment under the same conditions, the ratio condition score:back fat thickness was found to be 0.25 (R2 = 98%; Ingvartsen, unpublished). This was used to convert the ultrasound measurements to condition scores, which together with live weight values were used to calculate Lcalv. The average live weights at calving were 555, 604, and 650 kg for low, medium, and high groups, respectively. Conception occurred on average at 92 DIM.
In the second data set, condition score and live weight data from 65 Jersey cows were used. The data were collected within an experiment in which the cows received either a high or a low quality feed ad libitum throughout lactation. Data from the high quality feeding treatment only were used. (Data from the low quality feeding treatment were excluded as this treatment was designed to be nutritionally limiting and thus cause environmentally driven mobilization.) The data comprised 26 first, 23 second, and 16 third lactations. The average live weight at calving for first, second, and third parity was 403, 443, and 476, respectively. The average DIM of conception of the first, second, and third parity cows was 106, 105, and 102, respectively. Further details of this experiment have been reported (Nielsen et al., 2003).
The ability of the model to fit the observed data was estimated using the error of prediction, which is the square root of the ratio between the mean square of the prediction error (MSPE) and the observed mean. The MSPE was defined as follows:
![]() |
where Oi is the ith observed value, Pi is the ith predicted value and n is the total number of observations.
The error in central tendency (ECT), which provides a measure of the bias of the prediction, was also calculated as the square of the difference between the average of the predicted values and the average of the observed values. The ECT was used to calculate the percentage bias:
![]() |
| RESULTS |
|---|
|
|
|---|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Over 20 yr ago Bauman and Currie (1980) applied the concept of homeorhesis (Waddington, 1957; Kennedy, 1967) to lactation. The notion of "orchestrated or coordinated changes in metabolism of body tissues necessary to support a physiological state" (Bauman and Currie, 1980) was an important recognition of the fact that homeostasis alone cannot explain the dynamic associated with changing physiological state through pregnancy and lactation. However, the challenge of incorporating homeorhesis into the pre-existing (homeostasis based) view of lactational physiology, especially in quantitative terms, has all too frequently been avoided.
One possible reason for avoiding the idea of homeorhesis is that it requires a driving force or goal. The rationale presented in this paper suggests that the driving force is the expression of genotype through time to meet the evolutionary goal of maximizing reproductive success (see Friggens, 2003). The recognition that an orchestrated process implies a driving force has prompted alternative names for this process such as teleophoresis (see Chilliard et al., 2000; Bauman, 2000). In this paper, an attempt to describe teleophoresis as it affects body reserves has been made. This description is a high level description, it does not describe the metabolic and endocrine processes involved. There are 2 main reasons for this. The hormones that are involved in lipogenesis and lipolysis are important mediators of lipid metabolism, and thus lipid mobilization, but they are not the driving force behind lipid mobilization. In the context of an improved understanding of (dairy cow) physiology, study of the mediators such as the hormones addresses an important question; "how does lipid mobilization occur?" but it does not adequately address the question "why does lipid mobilization occur?" The other reason for choosing a high level description is that at this level there is sufficient information available to parameterize such a model for the practical purpose of improving prediction of energy requirements.
The model presented here is a simple method for predicting the genetically driven changes in body lipid. The model describes both the predetermined genetically driven mobilization and genetically driven reconvergence when a cow is regaining her predetermined level of body fatness after a nutritional insult. The model was found to readily reproduce at least some observed patterns of change in body lipid reserves through lactation. The evaluation of the model used 2 data sets from 2 distinct breeds of dairy cattle, Holstein and Jersey. The average error of prediction was low, less than the variation normally associated with the recording of condition score (Broster and Broster, 1998), and was similar for the 2 data sets. Further, in one data set (Holstein) nutritional manipulation prior to calving had been used to create 3 different levels of body condition at calving (Figure 6
). There was no significant difference between the 3 levels of calving condition score in the error of prediction of the condition score change after calving. Thus the model appears to adequately describe compensatory trajectories in body lipid.
In the Jersey data set, where cows from first, second, and third parities were represented (Figure 7
), the model produced an adequate fit for all 3 parities. However, the error of prediction was greatest in parity 3, which could possibly reflect a parity effect on the model parameter values, although this was the parity that also contained the fewest cows. The error of prediction for the 3 Jersey parities was improved, although not significantly, by decreasing T' from the initial value of 112 to 70 DIM. No such adjustment was needed in the Holstein data set. Given that these results were derived from only 2 experiments with relatively few animals, any conclusions drawn from these data about breed, or genotype, differences in the values of the parameters of the model must be tentative. This is an issue warranting further study. If the model is to be used to produce improved predictions of energy requirements then it becomes important to quantify the extent to which the different parameters are affected by genotype and parity. The same applies to possible future genetic selection for body lipid usage (e.g., Coffey et al., 2001). In this context, the finding that the model readily reproduces observed patterns of change in body lipid reserves through lactation suggests that the model will provide the basis for a more biologically meaningful genetic evaluation of genetically driven body lipid mobilization.
It is relevant to discuss the limitations of the model presented here. The model is not designed to predict environmentally driven changes in body lipid. In this situation, by definition, the cow does not have sufficient nutritional resources available to allow the demands of both milk production and genetically driven changes in body lipid to be met. Predicting body mobilization when resources are limiting requires the development of rules to describe the partitioning of nutrients between "competing" functions. The large differences which can exist between individual animals in the way they partition nutrients has long been recognized (Kellner, 1926, p243). An adequate quantitative understanding of partition is still lacking (Bauman, 2000), and, therefore, this case is out with the scope of the present model. However, an important step towards generating such models is being able to describe the underlying genetically driven partition of nutrients between milk and body reserves (Ingvartsen et al., 1999). The model presented here provides a basis for achieving this.
Within the model presented, a central assumption is that the rate of change of body mobilization (dL/dt) changes in a linear way with days from calving. The assumption of linearity was made on the basis of choosing the simplest functional form that resulted in realistic patterns of changing body lipid reserves through lactation. We could find no strong evidence to justify adoption of a more complex function. The consequence of this assumption is that the patterns changing body lipid reserves are described by quadratic functions. These functions appear to be adequate for practical purposes and have some attractive statistical properties. Estimates of genetic heritabilities and breeding values for condition score curves have recently become available (e.g., Pryce et al., 2000; Coffey et al., 2001; Koenen et al., 2001) and are, for similar reasons, all based on polynomial models of condition score relative to days from calving. However, it would be worthwhile to test the assumption of linearity with a purpose-built data set where the frequency and accuracy of body lipid measurements were sufficiently high as to allow a meaningful test of the assumption.
The underlying biological processes are, to some extent, taken into account in the current model by allowing for the effect of prior nutritional insult (as measured by Lcalv) on the trajectory of body lipid change. The model also makes the distinction between the body lipid change postcalving (up to T') and body lipid change postconception. This is, however, a rather crude representation of physiological state. Although it is convenient to describe the changes in body lipid reserves as a function of time, this is unlikely to be the causal force behind changing body lipid reserves. The evolutionary arguments alluded to in the Rationale suggest that genetically driven changes in body lipid reserves should be described as a function of changing reproductive priority (Friggens, 2003). Clearly, this is an issue that should be addressed as it offers the opportunity to generalize the model to other mammals and to deal with different reproductive strategies. Further, we believe that the principles applied here to body reserves have the potential to improve our ability to quantify other teleophoretic processes.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication July 17, 2003. Accepted for publication August 18, 2003.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
E. Kristensen, S. Ostergaard, M. A. Krogh, and C. Enevoldsen Technical Indicators of Financial Performance in the Dairy Herd J Dairy Sci, February 1, 2008; 91(2): 620 - 631. [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. R. Roche, D. P. Berry, J. M. Lee, K. A. Macdonald, and R. C. Boston Describing the Body Condition Score Change Between Successive Calvings: A Novel Strategy Generalizable to Diverse Cohorts J Dairy Sci, September 1, 2007; 90(9): 4378 - 4396. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. E. Pryce and B. L. Harris Genetics of body condition score in New Zealand dairy cows. J Dairy Sci, November 1, 2006; 89(11): 4424 - 4432. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. J. Kononoff, S. K. Ivan, W. Matzke, R. J. Grant, R. A. Stock, and T. J. Klopfenstein Milk production of dairy cows fed wet corn gluten feed during the dry period and lactation. J Dairy Sci, July 1, 2006; 89(7): 2608 - 2617. [Abstract] [Full Text] [PDF] |
||||