|
|
||||||||


* 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 |
|---|
|
|
|---|
Key Words: dairy bull energy balance genetics selection
Abbreviation key: conW = conceptus weight, HUK = Holstein UK, LWT = liveweight, PIN = profit index
| INTRODUCTION |
|---|
|
|
|---|
BCS is a technique of appraising the fat content of the body using visual and tactile methods (Lowman et al., 1976). There is a good relationship between BCS and total body fat content (Fox et al., 1999) such that BCS is a useful management aid for dairy farmers in monitoring the nutritional and metabolic status of cows. Furthermore, BCS levels and changes in BCS are associated with the health and fertility status of the cow (Veerkamp et al., 2000; De Vries et al., 2000; Collard et al., 2000; Pryce et al., 2001), so BCS is an obvious target for potential selection indices. Selection for yield alone has resulted in cows that have a lower BCS than cows of average genetic merit for production (Pryce et al., 1999). Using random regression techniques of analysis on field data with BCS measured at different times but only once on each animal, Jones et al. (1999) demonstrated that changes in BCS throughout lactation are under genetic control and that BCS has a heritability of around 0.3. De Vries et al. (1999) reported a similar analysis with data from the Netherlands. Coffey et al. (2001), using data from an experimental herd, showed that BCS and liveweight (LWT) could be combined into an overall genetic energy balance curve for sires based on measurements of their daughters. If energy balance curves can be calculated for sires using field data collected from their daughters, the opportunity arises to explore the genetic relationship between these curves and health, fertility, and survival traits. BCS is relatively easy and cheap to measure on large numbers of daughters via progeny testing and national conformation assessment schemes. Energy balance curves may provide data that could be included in a multitrait index aimed at improving health and fertility and thereby reducing wasteage from the dairy herd for both welfare and environmental reasons.
Liveweight is potentially a useful management aid (Maltz et al., 1997), yet it is not routinely recorded in UK commercial dairy herds. Liveweight can be predicted from linear type measurements and subsequent genetic evaluations for LWT can be obtained with little loss in accuracy (Koenen and Groen, 1998). Given that BCS is recorded as part of the UK linear type assessment scheme for Holsteins, there is the opportunity to calculate sire breeding values for LWT and BCS using large amounts of national data.
The objectives of this study were: 1) to model LWT predicted from type traits and condition score measured in the field to obtain daily breeding values for sires; 2) to combine those daily estimates into an overall genetic energy balance profile; and 3) to compare genetic energy balance curves among sires. To do this we used random regression techniques and restricted maximum likelihood (REML) on data from first lactation cows participating in the national conformation assessment scheme.
| MATERIALS AND METHODS |
|---|
|
|
|---|
![]() | ([1]) |
where yiklmt = liveweight recorded on cow i at day t of lactation. For animal i, yoik = the fixed effect of year of inspection, dietl = diet type (grazed grass or total mixed ration), aaim = age at inspection (grouped into 3-mo classes), and the linear covariate pch = percentage North American Holstein genes.
is the regression coefficient of LWT on percentage of North American Holstein genes. ßj are the regressions of LWT on the jth linear type trait, xj, and eiklmt is an error term. The cow was included in the model as a random effect since there were repeated observations on each cow.
National Linear Type Data and Condition Score
In the national conformation assessment scheme, information on 16 linear type traits and a number of other traits including BCS are scored on a 1 to 9 scale (Brotherstone et al., 1990). Participation in the overall scheme is voluntary although once a farm is committed, every heifer on the farm at the time of the visit must be inspected. Body condition score has been included in the scheme since August 1996, based on the scoring system of Lowman et al. (1976) but adapted to be scored on a 1 to 9 scale, where 1 is the lowest level of fatness and 9 is the highest level of fatness. Data for all heifers scored since August 1996 up to March 2001 were extracted from the HUK database. This produced 240,232 records. Scores for each trait were scaled by the ratio of the individual field officer standard deviation to the mean standard deviation for all field officers, in order to standardize the range of scores for all field officers (Brotherstone et al., 1990). After records were removed for animals inspected at day 0 or after day 305 of lactation, and from sire progeny groups of less than 10, there remained 58,784 records on daughters of 1240 young sires with first crop daughters and 10 proven sires. Proven sires were included in the analysis in order to improve genetic connections within the data but these sires did not contribute to the between-sire variance component. For each of the cows in the final dataset, LWT was predicted using equation [1]
.
Genetic and environmental variance components were estimated using a random regression sire model with the ASREML statistical package (Gilmour et al., 1998). The random regression model fitted in this study was:
![]() | ([2]) |
where yijkt is predicted LWT or BCS for daughters i of sire j in herd-visit k at day t of lactation. For daughter i of sire j, hvk = herd-date-of-visit interaction (a herd may be visited twice a year), pch = percentage North American Holstein genes (linear) and aac = age at calving in months (linear and quadratic).
1 and
2 are the regression coefficients for age at calving and proportion North American Holstein genes respectively. ßm are the fixed regression coefficients,
jm are the additive genetic random regression coefficients for sire j and
ijkt is the error associated with day t of lactation. Pm(t) is the mth legendre polynomial evaluated at day t of lactation and 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 (Draper and Smith, 1998). Different residual errors were associated with observations over time and, based on preliminary analyses, residual error classes were defined as days of lactation as shown in Table 1
. Within classes, residual errors were assumed to be homogeneous and between classes, residual covariances were assumed to be zero.
|
2 distribution with the number of degrees of freedom equal to the difference in the number of variance components to be estimated. Sire PTA obtained from the analysis were used to calculate daily values for each trait on the phenotypic scale for all sires in the dataset, for all traits for d of lactation 1 to 305. Energy balance was derived using predicted body protein and lipid changes 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. Details of the formulas used to convert to effective energy equivalents are given by Coffey et al. (2001).
The formula used to predict body lipid content from BCS is based on empty body weight (liveweight minus gut fill and weight of conceptus). The effects of pregnancy on the prediction of body lipid content were accounted for by modelling conceptus weight (ConW; fetus plus placenta plus fluid) at day p after conception using the formula from Bruce et al. (1984) which assumes a calf birth weight of 40 kg.
![]() |
The daily predicted weight of the conceptus was subtracted from LWT to ensure body lipid estimation was not biased upwards by the presence (weight) of conceptus. All sires were assumed to have average daughter fertility and the average daughter to be pregnant at 110 d of lactation. This expected day of conception corresponds to the average calving interval of around 390 d in the UK and ConW was modeled from day 110 of lactation. The effects of pregnancy and gut fill were assumed to be the same for all animals.
Dry matter intake (DMI) was predicted from the formula taken from NRC (1987)
![]() |
and was subsequently used to calculate gut fill (GF) for animal i on day t of lactation using the formula from Coffey et al. (2001)
![]() |
where MEC is the metabolizable energy content of the feed. A value of 11.82 was derived from the mean of all total mixed ration fed at Langhill for the last 10 years and was assumed to be indicative of the average MEC of dairy cattle feed available to a cow over one year.
| RESULTS |
|---|
|
|
|---|
![]() |
where aai was 0.0, 21.7, 34.5, 44.7, 69.2, or 68.7 for age at inspection in months of 23 to 25, 26 to 28, 29 to 31, 32 to 34, 35 to 37, and > 37, respectively. pch is percentage North American Holstein genes, Ang is angularity, CW is chest width, BD is body depth and Sta is stature. Data used to derive this formula were collected between 1990 and 2001. To generalize the formula so that it can be used with this period an average year of inspection solution has been included in the constant term. The variable diet had a value of 0 when the animals were eating grass and 34.79 when housed inside eating conserved forage. For national data it was assumed that animals were outside between mo 4 and 9 (April to September).
The correlation between actual and predicted LWT was 0.92, i.e., 85% of the variance in LWT was accounted for using the prediction formula. Predictions of LWT in the national population may be inferior in accuracy if the Langhill data used to derive the prediction formula are not representative of the general population. It would be useful to obtain additional data from other herds where LWT and linear type data were available to test the formula on cows not contributing data to the formula. Since these data were not available the correlation coefficient quoted is an upper limit.
For BCS, a significant improvement in fit resulted when the order of the polynomial for the random curve was increased from 3 to 4, indicating that a polynomial of order 4 (cubic) was better able to model the data than a quadratic polynomial. For LWT, there were convergence problems using order 4 for the random curve so order 3 was used.
Residual error classes (Table 1
) were based on earlier work (Coffey et al., 2002) and preliminary analysis of these data. Ten classes were defined ranging in size from 1053 records (d 1 to 15) to 8621 records (d 60 to 89). Values of error variance found in this analysis are larger than those obtained from a genetic analysis of experimental data (Coffey et al., 2001), but the data used here were taken from national records that are likely to be more variable than those from an experimental farm. Values of error variance for LWT are similar to those reported by Koenen and Groen (1998) and for BCS are similar to those reported by Koenen et al. (2001).
The daily average sire solutions for body energy loss, relative to the start of lactation (cumulative body energy state), for LWT and BCS are given in Figures 1
and 2
. For all sires, the overall mean maximum daughter body energy loss was 1499 MJ (SD = 144 MJ) and occurred at d 99 (SD = 12.8 d) of lactation, whereas the mean total daughter body energy loss at d 305 of lactation was 779 MJ (SD = 224 MJ).
|
|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Genetic variance for both traits used in this study rose abruptly towards the end of lactation. This had the effect of increasing the heritability in a similar fashion. At the minimum point of the trajectory, the heritability for LWT was about 0.30 and for BCS was about 0.25, which is similar to estimates reported elsewhere (Jones et al., 1999; Veerkamp et al., 2000). However, the rise in genetic variance at both ends of the trajectory increased the heritability to unusually high values of around 0.90 for LWT and 0.60 for CS. An increase in genetic variance at the right extreme of the trajectory could be because pregnancy effects had not been fully accounted for in this analysis and at both extremes because there were limited data at these points (Pool et al., 2000). Conversely, Fischer and van der Werf (2002) found, using simulated data, that the numbers of records at the extremities had little effect on the variance. In the study reported here, conceptus weight in pregnancy was modeled to remove the additional predicted weight of the conceptus. There may be other effects not accounted for, that alter either the shape of the animal, and hence the linear type scores used to predict liveweight, or the assessment of BCS, for example, fluid retention.
The relationship between body size (or LWT) and profitability is not clear. Sieber et al. (1988) reported that, in the US, taller cows produced more milk than shorter cows but lighter cows produced more fat corrected milk than did heavier cows in their first and later lactations. There is conflicting evidence regarding the benefit of using LWT in a selection index since its economic value is uncertain. In New Zealand, Spelman and Garrick (1997) show that inclusion of a negative economic weight for LWT optimizes profitability in a multitrait index and is expected to lead to a genetic reduction in LWT of 0.19 kg/year at a selection intensity of 1.95 SD. However, BCS is a measure that predicts the proportion of LWT (empty body weight) that is lipid, and reducing LWT at a fixed BCS will lead to a lower total amount of body lipid in the cow. Veerkamp (1998) concluded that the relationship between feed intake, LWT and efficiency is confounded by the relationship of yield and LWT with body tissue mobilization. The negative consequences of a reduction in body lipid content of cows at reduced LWT should be considered before incorporating size or LWT into a multitrait index.
Mean maximum body energy loss for progeny of all sires reported in this study (1499 MJ) occurred at d 99 of lactation and is similar to but, as expected, slightly higher than that reported by Tamminga et al. (1997) who found a total body energy loss of 1284 MJ up to d 56 of lactation. In study reported here, the profiles of body energy loss indicate that daughters of the majority of sires lost most body energy before d 150 and then recovered body energy, whereas a few continued to lose body energy through to the end of lactation. Current genetic evaluations for production do not account for energy contributions from body tissue mobilization, therefore some sires with high merit for production may have daughters with unacceptable body tissue mobilization profiles, leading to higher costs associated with poorer health and fertility. A method of accounting for this cost could be to correct yield for body tissue mobilization. In the energy system used in this study (Emmans, 1984), production of 1 kg milk (at 4.2% fat, 3.4% protein and 4% lactose) requires approximately 4.2 MJ effective energy. The cumulative body energy loss at d 305 of lactation was converted to the effective energy equivalent in kg milk for each sire. Using this method, the mean amount of milk accounted for by body energy loss was 189 kg milk, which is similar to the value of 324 kg milk quoted by Tamminga et al. (1997). Deducting the body energy equivalent of kg milk from the breeding value for milk for each sire resulted in a correlation of 0.98 between the ranking of sires before and after adjustment. However, some sires changed rank by large amounts, the largest being +355 positions. The four sires that moved up the ranking by the most places as a result of contributing body tissue to milk production when ranked on the adjusted milk PTA are given in Figure 10
. These four sires had 3211, 3380, 44, and 125 daughters, respectively. The four sires that decreased the most in the rankings are given in Figure 11
and these sires had 60, 20, 28, and 21 daughters, respectively. These groups of sires have clear differences in the profiles of body energy loss and gain of their daughters with those in Figure 10
actually gaining body energy at the end of lactation. Pryce et al. (2001) suggested that using BCS at some stages of lactation in an index could help alleviate the unfavorable effect of selection for yield on fertility. Further work is required to determine the most useful and informative parameters of the BCS curve over multiple lactations.
|
|
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication September 23, 2002. Accepted for publication January 13, 2003.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
J. K. Toshniwal, C. D. Dechow, B. G. Cassell, J. A. D. R. N. Appuhamy, and G. A. Varga Heritability of Electronically Recorded Daily Body Weight and Correlations with Yield, Dry Matter Intake, and Body Condition Score J Dairy Sci, August 1, 2008; 91(8): 3201 - 3210. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Oikonomou, G. E. Valergakis, G. Arsenos, N. Roubies, and G. Banos Genetic Profile of Body Energy and Blood Metabolic Traits Across Lactation in Primiparous Holstein Cows J Dairy Sci, July 1, 2008; 91(7): 2814 - 2822. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. Windig, B. Beerda, and R. F. Veerkamp Relationship Between Milk Progesterone Profiles and Genetic Merit for Milk Production, Milking Frequency, and Feeding Regimen in Dairy Cattle J Dairy Sci, July 1, 2008; 91(7): 2874 - 2884. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Wall, M. P. Coffey, and P. R. Amer A Theoretical Framework for Deriving Direct Economic Values for Body Tissue Mobilization Traits in Dairy Cattle J Dairy Sci, January 1, 2008; 91(1): 343 - 353. [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] |
||||
![]() |
S. McCarthy, D. P. Berry, P. Dillon, M. Rath, and B. Horan Influence of Holstein-Friesian Strain and Feed System on Body Weight and Body Condition Score Lactation Profiles J Dairy Sci, April 1, 2007; 90(4): 1859 - 1869. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Wall, M. P. Coffey, and S. Brotherstone The Relationship Between Body Energy Traits and Production and Fitness Traits in First-Lactation Dairy Cows J Dairy Sci, March 1, 2007; 90(3): 1527 - 1537. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Brotherstone, M. P. Coffey, and G. Banos Genetic Parameters of Growth in Dairy Cattle and Associations Between Growth and Health Traits J Dairy Sci, January 1, 2007; 90(1): 444 - 450. [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] |
||||
![]() |
D. Hojman, M. Gips, and E. Ezra Association Between Live Body Weight and Milk Urea Concentration in Holstein Cows J Dairy Sci, February 1, 2005; 88(2): 580 - 584. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Mark Applied Genetic Evaluations for Production and Functional Traits in Dairy Cattle J Dairy Sci, August 1, 2004; 87(8): 2641 - 2652. [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] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |