|
|
||||||||

* Animal Biology Division, Scottish Agricultural College, West Mains Road, Edinburgh EH9 3JG, UK
Institute of Cell, Animal and Population Biology, University of Edinburgh, West Mains Road, Edinburgh, EH9 3JT, UK
Corresponding author:
M. P. Coffey; e-mail:
m.coffey{at}ed.sac.ac.uk.
| ABSTRACT |
|---|
|
|
|---|
Energy balance over three lactations was modelled using sinusoidal functions which were associated with individual cows and allowed to vary between cows. The parameters of these curves are potentially useful since they have a biological interpretation. The phase relates to the period from calving to return to positive energy balance, and the amplitude relates to the degree of body energy loss (and recovery). The sinusoidal functions fitted to the curve removed a significant proportion of the variation, but accounted for only 45% and 40% of the variation in EB1 and EB2, respectively.
The relationship between energy balance in the first three lactations is likely to be more complex than a simple linear function, but the profile of energy balance over the first three lactations may be a useful selection criteria in a multi-trait index. Energy balance profile over lactations one to three can be modelled with moderate accuracy using sinusoidal functions, and this warrants further research.
Abbreviation key: CS = condition score, EB1 = energy balance calculated from feed intake and milk production records, EB2 = energy balance calculated from body energy state changes, FI = feed intake, LWT = live weight, MY = Milk yield
| INTRODUCTION |
|---|
|
|
|---|
The replacement of body tissue in modern dairy cows usually takes place later in the same lactation, once milk output begins to decline. If the cow is pregnant, yield is further depressed (Olori et al., 1997) as the fetus develops. Replenishment of body tissue in preparation for the next lactation creates competition for nutrients among the demands for current yield and the developing fetus. The use of body lipid as a nutritional buffer is a normal mammalian physiological process (Pond et al., 1999), but the biological rules determining priorities of use, or partitioning, of available nutrients in cows of different genetic merit for production and at different stages of lactation and gestation are unknown.
The term "energy balance" is often used to describe the body energy state of dairy cows, which is the outcome of daily energy flux; negative energy balance is associated with body energy loss and positive energy balance with body energy gain. Cows which lose body tissue, and hence energy, in early lactation usually return to positive energy balance at around 40 to 80 days post partum (Sutter et al., 2000; Coffey et al., 2001; Veerkamp et al., 2000). Kendrick et al. (1999) found that return to positive energy balance was at 21 days or 49 days with high or low energy diets, respectively. However, cumulative body energy loss in the first lactation is, on average, only fully recovered at around day 200 (Coffey et al., 2001). Continued growth is required through the dry period and into the second lactation in order to reach physiological maturity. Negative energy balance is related to some health traits (Collard et al., 2000), to resumption of reproductive activity (De Vries et al., 1999; Veerkamp et al., 2000), and to oocyte size and quality (Beam and Butler, 1999; Kendrick et al., 1999). Individual cows may not regain all lost body energy in the first lactation leading to a greater deficit to be replenished in the second. This may result in carryover effects from one lactation to the next, and on both yield and nonyield traits such as health and fertility.
It is important to view productivity over the cows entire lifetime and not simply on the basis of a single lactation. Current and previous body energy state changes may create a legacy that affects current and future productivity and health and fertility. Parameters of the lactation curve are different between at least the first three lactations, although curve parameters from one lactation can be used to predict those of subsequent lactations. Friggens et al. (2000) found that parameters associated with peak lactation and the rate of decline after peak in the third lactation can be described as simple ratios of those parameters in the first and second lactations. This is important when considering energy balance across lactations for the same animals when measures are available only in the first lactation, or where subsequent lactation measures may be biased by selection. The use of random regressions and Fourier series allows modelling of cyclical changes in a trait over seasons (Meyer, 2000) and could be used to model lifetime (or at least multiple lactation) energy balance changes in dairy cows. Parameters of these cyclical functions might then be analyzed for relationships with traits of economic importance such as health, fertility, and survival.
The objectives of this study were: 1) to model phenotypic daily milk yield, fresh feed intake, live weight and condition score using random regression techniques and restricted maximum likelihood (REML) for traits measured at different times and with different frequencies on dairy cows from the Langhill Dairy Cattle Research Centre; 2) to combine daily estimates obtained in the first objective into an overall energy balance for each day of the first three successive lactations; 3) to compare energy balance curves over three lactations; and 4) to investigate the feasibility of modelling energy balance in the first three lactations using harmonic analysis.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
![]() |
where Fit represents fixed effects of genetic line (2 groups), feed group (2 groups), time of measurement (year and week of measurement), the covariates percentage North American Holstein genes (linear) and age at calving in months (linear and quadratic) for animal i. ßm are the fixed regression coefficients;
im are the random regression coefficients associated with the animal plus its permanent environment; and
it is the residual error associated with days since calving t. Pm(t) is the mth Legendre polynomial evaluated at t. The parameters f and k are the order of the fixed and random polynomials, respectively. Legendre polynomials were used because they are easy to manipulate, have good convergence properties and, being orthogonal polynomials, correlations between coefficients are lower than between the coefficients of ordinary polynomials. The four traits analyzed separately by lactation were: milk yield, condition score, feed intake, and live weight. Residual, or measurement, errors were expected to have heterogeneous variances over each lactation, with larger variances at the beginning of lactation and around peak yield. Different residual errors were, therefore, associated with observations over time. Based on preliminary analyses, residual error classes were defined as days 1 to 6, 7 to 9, 10 to 12, 13 to 15, 16 to 29, 30 to 99, and 100 to 305. Within classes, residual errors were assumed to be homogeneous. Ideally, we would have preferred to model the residual variances using a continuous function, but an attempt to do so failed due to problems with convergence.
Based on previous genetic analyses of these data (Coffey et al., 2001), fixed regressions, which model the general shape of the curve and are common to all animals, were fitted for all traits as polynomials of order 5. Third and fourth order polynomials were used to model the animal effect, and the goodness of fit was compared using a likelihood ratio test based on the
2 distribution, and a significance level of 0.05 with the number of degrees of freedom equal to the difference in the number of variance components to be estimated.
Animal solutions obtained from the analysis were used to calculate daily phenotypic values for all animals in the data set, for all traits for days of lactation 1 to 305. Energy balance was derived in two different ways for the same cows after converting all measures to energy equivalents using the effective energy system of Emmans (1994) with additional terms for the major organic components of milk. The first method of deriving energy balance (EB1) was based on predicted energy available from recorded feed consumption minus energy required for recorded daily milk production, and maintenance predicted from LWT and gut fill. The second method (EB2) used body protein and lipid changes predicted from LWT and CS. Details of the formulas used to convert traits to effective energy equivalents are given by Coffey et al. (2001). Using these formulas allows a comparison of methods of calculation and sources of data. Applying these formulas to data from successive lactations also allows an assessment of the usefulness of first lactation energy balance in predicting energy balance in the second and third lactation for each of the methods of calculation.
The component values used in the calculation of energy content of feed and milk were determined in two ways: as the average over the whole data set used for each trait, and from results of feed and milk analyses concurrent with the time period in which the food was eaten or milk produced. This enables a comparison of energy balance calculated using different data sources for the subsequent purpose of using national data where feed intake and feed energy analysis are not available. Feed samples were taken daily and pooled for weekly analysis of components. Milk samples were obtained weekly from both daily milkings, and analyzed separately by an official recording agency (Cattle Information Service, Scotland) using a Foss Milkoscan 605 (Foss, Denmark). Milk energy output for each day was calculated from daily yields, and weekly milk sample analysis and feed energy intake was calculated from feed intake and the latest feed analysis.
The effects of pregnancy on the prediction of body weight were accounted for in part by modelling conceptus total weight (fetus plus placenta plus fluid) using an exponential growth curve from day of conception. Jakobsen et al. (1957) showed that energy requirements for the conceptus of heifers rose exponentially through gestation. For lactations 1 to 3, this curve resulted in a weight of conceptus at 280 days of gestation of 80, 90, and 100 kg respectively to account for assumed increases in weight of this component in larger cows. The daily predicted weight of conceptus was subtracted from empty body weight to reduce any upwards bias on body lipid estimation by the presence of conceptus.
A visual appraisal of a sample of energy balance curves over three lactations for individual cows suggested that a possible method of describing the fluctuations in energy balance across lactations might be by using sinusoidal functions. All cows calved approximately every 12 months in order to remain on the trial, so energy balance data were collected over the three lactation cycles approximately every 365 days. As cows regain body condition after calving at different rates, (i.e., the phase varies between cows), and cows have different total amounts of body energy loss, (i.e., the amplitude varies between cows), sinusoidal functions were associated with individual cows and allowed to vary between cows.
The model for the harmonic analysis of the data was:
![]() |
where yi = energy balance for animal i on day T (either EB1 or EB2); T = days since first calving; a and b are fixed regression coefficients which cater for a trend over time; ci and d iare the random regression coefficients; and ei = error term for animali.
| RESULTS |
|---|
|
|
|---|
Residual error variances are given in Table 2
for all error classes (days) for each trait in each lactation. For MY, the residual error variance increased for each class by lactation number in line with increased yield from later parity cows. For MY and FI, residual error variance declined as lactation progressed; whereas, for LWT it declined and then increased in the last error class. In lactation 1, the residual error variance for CS declined in early lactation and increased up to day 99. In the last measurement error class (100–305 days), it declined dramatically. In lactations 2 and 3, residual variation for CS followed a similar pattern to the other traits, although it did not decline in late lactation. Values of residual error variance found in this analysis are similar in magnitude to those obtained from a genetic analysis (Coffey et al., 2001) for the last two classes. These latter classes are the only ones that approximated to similar time periods in lactation in both analyses.
|
|
|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The condition scoring method of Lowman et al. (1976) uses a scale of 0 to 5 for cows that may vary in weight. Therefore, CS predicts body lipid in a proportionate manner; therefore, body lipid estimation is sensitive to prediction or measurement of body weight. In the formula used in this study, body lipid is predicted from empty body weight which is dependent on unbiased estimates of gut fill. Towards the end of gestation, the conceptus accounts for an increasingly larger amount of body weight and may affect gut fill. Therefore, the conceptus must be properly accounted for in order to minimize error in prediction of body lipid using CS and empty body weight. In this study, the predicted weight of conceptus was subtracted from body weight before EB2 calculation. Energetic consumption by the conceptus was ignored, but this would only have a small effect on EB1 calculation.
The effect of pregnancy on the individual traits was ignored in this study since the data were obtained from block-calving cows. All cows used in this study were within 4 months of calving of each other and were at approximately the same stage of lactation when pregnant. Therefore, the interaction between stage of lactation and gestation was assumed to be the same for all cows, leading to no bias. Olori et al. (1997) showed that yield was significantly reduced after the fifth month of gestation, but that yield was only reduced by more than 1 kg milk/day after 6 months of gestation. However, their study reported only the effects of pregnancy on milk yield and its components. It would be useful to know if there are significant effects of pregnancy on other traits, such as feed intake and particularly on body energy state. Of relevance to this study would be any three-way interaction of day of lactation, current or previous body energy state, and day of pregnancy. This would provide information on the genetic control of nutrient partitioning in animals of different energy status when those animals depart from their preferred energy status as defined by their genetic merit for body energy at that stage of lactation (Jones et al., 1999).
It has been reported that post partum reproductive activity may only resume once the nadir of negative energy balance has been reached (Butler and Smith, 1989; De Vries et al., 1999; Veerkamp et al., 2000), indicating that the rate of return to positive energy balance may be a useful indicator of resumption of reproductive activity and, by implication, possibly a useful selection objective. In this study, using EB1, cows returned to positive energy balance at days 72, 75, and 95 in lactations 1, 2, and 3, respectively. Using body state changes as an indicator, EB2 return to positive energy balance was at days 77, 83, and 73 for lactations 1, 2, and 3. This discrepancy between the methods for calculating energy balance might be due to the inadequacy of the formulas to predict body lipid content from CS observations, the impact of estimation of gut fill on the prediction of body lipid from CS observations, bias in the observation of CS in different parity animals, inaccuracies in the estimation of dietary energy content, inaccuracies in the estimation of milk energy or maintenance energy output, or a combination of these. It would be useful for future research in this area to have revised values for the relationship between body lipid content and CS observations in modern dairy cows of different weight, stage of lactation, yield level, and conformation. This would enable more precise use of CS observations for farm management purposes and for research where energy balance is needed.
The rate of body energy loss and gain relative to the start of each lactation appears to change as lactation number increases (Figure 7
). In the first lactation there will be a component of lipid-free growth that will require nutrients and increase the weight of the animal. This may slow down the apparent accretion of body fat since ingested energy will be partitioned to growth rather than deposited in fat depots that are measured by CS.
The correlation between energy balance on the same day of lactation in different parities is an indication of how accurately energy balance calculated from data gathered in one lactation predicts the subsequent lactation energy balance. The correlations between days in successive lactations change in their profile across lactation (Table 3
). For EB1, correlations were lower at the beginning of lactation, rose to peak at around day 150, then dropped slightly as lactation progressed. This peak in correlation corresponds approximately to the point of minimum condition score and may indicate that minimum condition score is less dependent on management or environment than is early or late lactation condition score when yields are lower. For EB2, the highest correlation occurred around 90 days, but followed a similar pattern as EB1.
Changes in energy balance over the lifetime of a cow may be a useful future selection objective since there are genetic differences in energy balance profiles between bull daughter groups in their first lactation (Coffey et al., 2001). A preliminary investigation indicated that lifetime energy balance may be modelled with low accuracy using sinusoidal functions. The parameters of the curves are useful since they have a biological interpretation. The phase relates to the period from calving to return to positive energy balance, and the amplitude relates to the degree of body energy loss (and recovery). Clearly, a declining amplitude over successive lactations means the cow is losing body energy over time and eventually may be compromised due to insufficient body energy reserves.
Here we used first order sinusoidal functions to model the periodic fluctuations over three lactations in an attempt to model energy balance over three lactations as a continuous trait even though the measurements do not span the dry period between lactations. We fitted a single sine and cosine function that did not produce a perfect fit. The lack of fit was associated with the fact that while the curve of the underlying biological process was cyclic, its shape was not symmetrical nor exactly sinusoidal, as there are differences in slope for the ascending and descending phases of the curve. Fitting a sine/cosine function creates a symmetrical curve and smooth continuum between lactation and dry period when the trait may, in fact, be discontinuous. Adding further sine and cosine terms with higher frequencies did not successfully eliminate the lack of fit. On the contrary, although with the addition of higher order sine and cosine terms the fit improved, the fitted curve became ripply in appearance. A further complication is the lack of data throughout the dry period. Additional data on traits affecting energy balance collected during the dry period should enable us to identify a more suitable function and improve the fit.
A more detailed analysis is warranted in order to refine the fit of the curve, and a possible method of achieving this is to use circular splines. Smoothing splines are constrained by boundary conditions so that in the interval (a,b) the value at a is equal to the value atb . Analysis using circular splines is the same as harmonic analysis except that the higher frequency terms are given a lower weighting, resulting in a smooth curve. Although splines allow flexibility in the shape of the curve (White et al., 1999), the inclusion of these higher frequency terms makes it difficult to attribute any biological meaning to the coefficients of the curve. A biological interpretation of the coefficients is useful in that it allows us to make meaningful associations between them and other characteristics of importance, such as longevity, health, and reproductive success.
In this study, weekly milk sample analysis and feed analysis was available, and so, EB1 and EB2 were calculated using average values (Figures 5
and 6
) and actual values for each cow for each day (Figures 8
and 9
). This comparison enabled an assessment of the potential use of national measures of energy balance derived from linear type traits, condition score, and average values of milk and feed energy. The similarity in shape of these two sets of graphs demonstrates that on average, a single lactation value for feed energy and milk energy can be used when estimating average bull daughters energy balance. This is important if EB2 is used to calculate energy balance from national data, since the estimation of empty body weight depends, in part, on an estimation of gut fill, which itself depends on estimation of feed analysis of digestibility.
A notable difference between EB1 and EB2 in both sets of graphs (Figures 5
and 6
and Figures 8
and 9
) is that the rate of return to positive energy balance is lower for EB2. This would imply a lag between the animals true energy state derived from estimates of energy flux using feed energy intake and milk energy output and its energy state predicted from body fat changes. An alternative suggestion is that body fat is being mobilized from fat stores assessed using CS, and then retained within the body. This repartitioning of nutrients may serve to support internal organs associated with increased metabolic activity and lactation.
Figures 8
and 9
also demonstrate more clearly the harmonic nature of the energy balance curve across lactations. At the end of each lactation, the curve tends towards zero energy balance as the animal approaches its desired energetic state. However, these data do not contain measurements taken during the dry period and, therefore, curve fitting is more difficult. Future experimentation to record appropriate measures through the dry period on individual cows would enable verification of the usefulness of fitting sinusoidal curves for energy balance throughout the animals productive life.
| CONCLUSIONS |
|---|
|
|
|---|
|
|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication June 26, 2001. Accepted for publication November 26, 2001.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
E. Kramer, E. Stamer, J. Spilke, G. Thaller, and J. Krieter Analysis of water intake and dry matter intake using different lactation curve models J Dairy Sci, August 1, 2009; 92(8): 4072 - 4081. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Bewley and M. M. Schutz Review: An Interdisciplinary Review of Body Condition Scoring for Dairy Cattle Professional Animal Scientist, December 1, 2008; 24(6): 507 - 529. [Abstract] [PDF] |
||||
![]() |
D. P. Berry, B. Horan, M. O'Donovan, F. Buckley, E. Kennedy, M. McEvoy, and P. Dillon Genetics of Grass Dry Matter Intake, Energy Balance, and Digestibility in Grazing Irish Dairy Cows J Dairy Sci, October 1, 2007; 90(10): 4835 - 4845. [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] |
||||
![]() |
G. Banos, M. P. Coffey, E. Wall, and S. Brotherstone Genetic relationship between first-lactation body energy and later-life udder health in dairy cattle. J Dairy Sci, June 1, 2006; 89(6): 2222 - 2232. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. P. Coffey, J. Hickey, and S. Brotherstone Genetic Aspects of Growth of Holstein-Friesian Dairy Cows from Birth to Maturity J Dairy Sci, January 1, 2006; 89(1): 322 - 329. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Wall, M. P. Coffey, and S. Brotherstone Body Trait Profiles in Holstein-Friesians Modeled Using Random Regression J Dairy Sci, October 1, 2005; 88(10): 3663 - 3671. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Banos, S. Brotherstone, and M. P. Coffey Genetic Profile of Total Body Energy Content of Holstein Cows in the First Three Lactations J Dairy Sci, July 1, 2005; 88(7): 2616 - 2623. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Banos, M. P. Coffey, and S. Brotherstone Modeling Daily Energy Balance of Dairy Cows in the First Three Lactations J Dairy Sci, June 1, 2005; 88(6): 2226 - 2237. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. P. Coffey, G. Simm, J. D. Oldham, W. G. Hill, and S. Brotherstone Genotype and Diet Effects on Energy Balance in the First Three Lactations of Dairy Cows J Dairy Sci, December 1, 2004; 87(12): 4318 - 4326. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Banos, S. Brotherstone, and M. P. Coffey Evaluation of Body Condition Score Measured Throughout Lactation as an Indicator of Fertility in Dairy Cattle J Dairy Sci, August 1, 2004; 87(8): 2669 - 2676. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Wall, S. Brotherstone, J. A. Woolliams, G. Banos, and M. P. Coffey Genetic Evaluation of Fertility Using Direct and Correlated Traits J Dairy Sci, December 1, 2003; 86(12): 4093 - 4102. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. P. Berry, F. Buckley, P. Dillon, R. D. Evans, M. Rath, and R. F. Veerkamp Genetic Parameters for Body Condition Score, Body Weight, Milk Yield, and Fertility Estimated Using Random Regression Models J Dairy Sci, November 1, 2003; 86(11): 3704 - 3717. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. P. Coffey, G. Simm, W. G. Hill, and S. Brotherstone Genetic Evaluations of Dairy Bulls for Daughter Energy Balance Profiles Using Linear Type Scores and Body Condition Score Analyzed Using Random Regression J Dairy Sci, June 1, 2003; 86(6): 2205 - 2212. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |