|
|
||||||||

* Statistical Animal Genetics Group, Institute of Animal Science, Swiss Federal Institute of Technology, ETH Zentrum, Zurich CH 8092, Switzerland and
Holstein Association of Switzerland, Grangeneuve, 1725 Posieux, Switzerland
Corresponding author: H. Kadarmideen; e-mail: haja.kadarmideen{at}inw.agrl.ethz.ch.
| ABSTRACT |
|---|
|
|
|---|
Key Words: body condition score body type genetic analysis milk production
| INTRODUCTION |
|---|
|
|
|---|
Dairy cow fertility is an economically important trait, but heritability is typically low, around 0.03, and many countries have data quality problems (Kadarmideen et al., 2000, 2003; Kadarmideen and Coffey 2001; Kadarmideen 2003). Body condition score is easily and cheaply recorded as a part of type classification, has a heritability similar to milk yield (range 0.25 to 0.35), and has a favorable genetic/phenotypic correlation with fertility traits (Veerkamp et al., 2001; Dechow et al., 2002; Lopez-Gatius et al., 2003). It may also be possible to use BCS as an indicator trait for (metabolic) diseases and disorders because BCS is strongly related to energy balance and longevity (Coffey et al., 2001).
In Switzerland, the Holstein Association started recording BCS in 2001, and these data have never been analyzed. Heritabilities for some type traits in Swiss Holsteins have been estimated, but their genetic/phenotypic correlations with BCS have never been investigated. Very few international studies have reported genetic correlation of BCS with as many conformation and production traits as are considered here. It is of interest to Holstein and other cattle-breeding associations of Switzerland (Brown Swiss and Simmental) to estimate genetic and phenotypic parameters of BCS data collected up to this point and to investigate the genetic and other relationships of BCS with production and type traits that are in the national economic selection index. Therefore, the main objectives of this study were to estimate genetic and phenotypic parameters of BCS and its genetic and phenotypic correlation with type and milk production traits.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
![]() |
where y1 was a vector of records of type or milk production traits and y2 was a vector of BCS records, respectively; X1 and X2 were design matrices relating fixed effects in b1 and b2 to y1 and y2, respectively. The fixed effects in the model for milk production (b1) and BCS (b2) were:
![]() |
and the fixed effects for the conformation traits were the same as those used for BCS (b2), where µ = the overall mean, HYC = herd-year of calving, MC = month of calving, ac = age (in days) at calving as a covariate, HYV = herd-year of visit by classifier, LS = stage of lactation (in months) at the time of classification, MT = month of type classification, at = age (in days) at condition scoring or type classification as a covariate, and hp = percentage of Holstein genes as a covariate.
The vector of random effects of sire was included in s1 and s2, where the design matrices, Z1 and Z2, relate records to these sire effects. The vector of random residuals were e1 and e2 for y1 and y2, respectively. The expectation (E) and variances (V) of model terms were as follows:
![]() |
where
![]() |
and where A is the additive relationship matrix and
is the Kronecker product. Also, R is residual variance-covariance matrix with
![]() |
where,
![]() |
and I is the identity matrix. Each trait was assumed to follow multivariate normal distributions with the above mean and variance-covariances. We assumed no covariance between residuals and genetic effects; cov (s,e) = 0. It should be noted that BCS data is missing on some animals; therefore, the residual elements corresponding to rows and columns of animals that has missing BCS observations in R0 are set to zero. Handling of missing observations for one or more traits in multi-trait mixed model equations are described by Henderson (1984).
Data were also analyzed by univariate models, which consisted of the same effects as those fitted for individual traits in bivariate settings (in b1 and s for production traits and b2 and s for BCS and type traits). Univariate results were used for obtaining starting values for variance components estimation by the bivariate models, for estimating least squares means for fixed effects of BCS and to examine the distribution of BCS EBV.
Implementation and software.
Variance-covariance parameters for all models were estimated using the software package, ASREML (Gilmour et al., 2001). Heritabilities (h2) and all correlations (genetic-rg, residual-re, and phenotypic-rp) were computed using estimated variance-covariance matrices. In ASREML, it was possible to get standard errors for all estimated parameters or its ratios. In total, there were 61 blocks of genetic analyses (31 for univariate and 30 for bivariate analyses). It was not our aim to compute correlations (covariances) among type traits or between type and milk production traits.
| RESULTS |
|---|
|
|
|---|
|
Lactation stage and month of classification.
Least squares means for the month of classification and stage of lactation, derived from univariate analysis of BCS, are given in Table 2
. Least squares means of BCS for month of classification effect (significant at P < 0.05 for all levels except the month of July) indicated that cows that were classified in January or February tended to have, on average, smaller scores than those condition scored in December. Raw means computed from the data showed May and August had lowest and highest scores, respectively.
|
The BCS loss or gain can be better visualized in weeks rather than months because some of the BCS changes during a crucial part of lactation seem to happen in a matter of weeks. Figure 2
shows that the mean BCS by stage of lactation (in weeks) from calving date. There were only 7 cows classified in wk 1, but a reasonable number of cows was classified in later weeks (98 in wk 2, 284 in wk 3, 386 in wk 4, etc.). As seen in Figure 2
, the drop in body condition was fast in the first 12 wk, moderate for the following 10 wk, and then cows seemed to put on fat reserves after 22 wk. Similar phenotypic trends (in months or weeks) have been observed in earlier studies (Jones et al., 1999; Pryce et al., 2001; Berry et al., 2002).
|
Sire EBV for BCS
The distribution of BCS EBV for 857 sires (in the pedigree) is plotted in Figure 3
. The minimum and maximum of sire EBV were -0.46 and +0.51, respectively (Figure 3
). The standard deviation of EBV was 0.09, and about 95% of the bulls EBV ranged between -2 to +2 standard deviation units. This can be seen in Figure 3
, where there is a higher peak at the center of the distribution (Figure 3
) compared with the normal distribution. When distribution of EBV was tested for deviation from normality, there was a positive kurtosis of 2.75, showing that most sires EBV were in fact very close to the mean. Although, there is not much ideal spread of BCS EBV among sires, the estimated heritability of 0.24 (Table 2
) would help in distinguishing worst/best sires for daughters BCS, and thus help in selection and breeding.
|
|
Body Condition Score and Milk Production Traits
Genetic correlations of BCS with milk production traits (with SE) are given in Table 4
. They were -0.12 for both milk and protein yield and -0.17 for fat yield. Estimates for milk and protein yield were not highly significant. However, these results indicate that increased genetic merit for milk, fat, and protein yield would result in decreased genetic merit for BCS in the long-term selection. Estimates of residual and phenotypic correlations are also given in Table 4
. Estimates of residual correlations were lower than estimates of genetic correlations (around 0.03) and negative. Estimate of phenotypic correlation between BCS and milk yield showed that increased milk yield, fat, and protein yield are associated with lower BCS. The standard error was 0.01 for all estimates of residual correlations and phenotypic correlations.
|
| DISCUSSION |
|---|
|
|
|---|
In Switzerland, the Holstein Association started recording BCS as part of the type classification scheme in 2001. Here, we report estimates of genetic and environmental effects that affect BCS and their association with those of body conformation and milk production traits, using the BCS data collected up to this point. Such genetic parameters (heritabilities and genetic correlations) are needed to construct new or modify existing Swiss selection indexes, which are aimed at selecting "genetically optimal" animals, giving due economic weight to production and type traits. In this study, genetic and phenotypic parameters were estimated for BCS, type, and production traits, for that purpose.
Regression of BCS on Holstein genes showed that an increase in Holstein genetics decreases body condition. This trend has been observed by others (e.g., Veerkamp et al., 2001; Berry et al., 2002). Regression estimates of age at classification showed that BCS increases with age, and many other type traits also changed as a function of age (in days). These results are to be expected, but here we quantified such changes in body composition or type as a function of age.
Least squares means for 9 stages of lactation (Table 2
) and BCS loss or gain in weeks (Figure 2
) indicated that, at the time of calving, cows start with good body condition, lose condition progressively up to approximately 22 wk (about 5 mo) of calving, and then tend to gain body condition to maintain later stages of pregnancy and calve again. Similar phenotypic trend (in months or weeks) has been observed in earlier studies (Jones et al., 1999; Pryce et al., 2001; Berry et al., 2002). All estimates of fixed effects for BCS reported here were derived from univariate analysis of BCS which when compared to estimates from series of bivariate analysis of BCS with type or production traits were similar.
Heritability of BCS obtained in this study was 0.24, which is similar to estimates reported by Pryce et al. (2001), Gallo et al. (2001), and Berry et al. (2002). There were some slight differences between our estimates and estimates of Dechow et al. (2002) and Veerkamp et al. (2001), but these differences may be due to differences in definition of BCS (some studies defined and estimated parameters for BCS as longitudinal measurements or used a different scale of measurements: 1 to 5 vs. 1 to 9), sire vs. animal models, standard deviations of the traits analyzed, and different data editing and estimation procedure used.
Estimates of heritability for many type traits are similar to most literature estimates (Rupp and Boichard 1999; Gallo et al., 2001). Heritabilities for type traits reported in this study were, in general, similar to those used in genetic evaluation by the Swiss Holstein Association (S. Wegmann, personal communication).
Direction and magnitude of genetic correlation of BCS with many type traits (e.g., heart girth, rear leg side view, fore udder attachment, udder depth, and teat position rear) found in this study were similar to earlier literature estimates (Veerkamp and Brotherstone 1997; Pryce et al., 2001). Genetic correlations between BCS and 305-d milk production traits were, in general, unfavorable, so a genetic decline in body condition is to be expected if selection for production alone continues in future. Some studies have reported these antagonistic genetic correlations between BCS and yield traits on a more severe scale (Pryce et al. [2001] for the U.K.; Berry et al. [2002] for Ireland; approximately -0.63, -0.41, and -0.54 for milk, fat, and protein yield, respectively) or on a moderate scale of approximately -0.20 to -0. 30 (Veerkamp et al. [2001] and Dechow et al. [2003] for Dutch and U.S. data, respectively). Our estimates (-0.12 to -0.17) are probably underestimated due to having missing BCS observations in bivariate analyses, using sire models instead of animal models or different estimation methods, or a combination of these reasons.
Estimates of genetic and environmental parameters for BCS reported in this study should be validated after collecting more years of BCS data and reanalysis. For example, very small differences in LS means for stage of lactation effect (Table 2
), narrow distribution of sire BCS EBVs (Figure 3
), nonsignificant genetic correlation of milk or protein yield with BCS (Table 4
) may change when analyses are repeated with few more years of additional data.
In this study, we considered a single observation of BCS per cow and correspondingly estimated genetic parameters for this cross-sectional data. It is also of interest to know how body condition/composition changes over the trajectory of lactation according to the demands for production and maintenance. In this case, we would have to record BCS at different stages of lactation at the cow level, and it would be almost impossible on a national scale. One approach to collecting these longitudinal data for BCS would be to collect data on daughters of sires, with condition score at different stages of lactation, such that sires genetic merit is estimated based on their daughters data sorted by stages of lactation. Random regression models for longitudinal BCS data was first applied by Jones et al. (1999) to U.K. data and later applied by Veerkamp et al. (2001) to Dutch data. As more data accumulate for BCS sufficient for sire by stage of lactation classes, this random regression model would be a preferred approach as this can overcome the problem of limited data at the cow level, as well as improve accuracy of predictions. In future, we plan to apply random regression models to our datasets.
Kadarmideen et al. (2003) showed that single-trait genetic evaluation for fertility will be biased and recommended that fertility breeding values be calculated from a multitrait analysis that includes milk yield as a correlated trait. A similar principle would also be applicable to genetic evaluation for BCS because farmers may select cows for breeding based on their genetic merit for milk yield (in which case, genetically thin animals are preferentially treated). This bivariate analysis of BCS with milk yield needs further investigation. Type traits and longevity have a moderate genetic correlation (Vollema et al., 1997). In our study, we showed that BCS is genetically related to stature, heart girth, capacity, and to overall scores for some conformation traits, such as udder, feet, and legs. This indicates that selection for optimal body condition together with these important conformation traits is likely to increase longevity. The use of BCS as a predictor of longevity, however, depends on how much information BCS provides over and above those found in these conformation traits.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication June 8, 2003. Accepted for publication August 4, 2003.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
R. Dal Zotto, M. De Marchi, C. Dalvit, M. Cassandro, L. Gallo, P. Carnier, and G. Bittante Heritabilities and Genetic Correlations of Body Condition Score and Calving Interval with Yield, Somatic Cell Score, and Linear Type Traits in Brown Swiss Cattle J Dairy Sci, December 1, 2007; 90(12): 5737 - 5743. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. N. Kadarmideen and L. L. G. Janss Population and systems genetics analyses of cortisol in pigs divergently selected for stress Physiol Genomics, March 14, 2007; 29(1): 57 - 65. [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] |
||||
![]() |
B. Karacaoren, F. Jaffrezic, and H. N. Kadarmideen Genetic Parameters for Functional Traits in Dairy Cattle from Daily Random Regression Models J Dairy Sci, February 1, 2006; 89(2): 791 - 798. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Neuenschwander, H. N. Kadarmideen, S. Wegmann, and Y. de Haas Genetics of Parity-Dependant Production Increase and its Relationship with Health, Fertility, Longevity, and Conformation in Swiss Holsteins J Dairy Sci, April 1, 2005; 88(4): 1540 - 1551. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sewalem, G. J. Kistemaker, and B. J. Van Doormaal Relationship Between Type Traits and Longevity in Canadian Jerseys and Ayrshires Using a Weibull Proportional Hazards Model J Dairy Sci, April 1, 2005; 88(4): 1552 - 1560. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sewalem, G. J. Kistemaker, F. Miglior, and B. J. Van Doormaal Analysis of the Relationship Between Type Traits and Functional Survival in Canadian Holsteins Using a Weibull Proportional Hazards Model J Dairy Sci, November 1, 2004; 87(11): 3938 - 3946. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |