|
|
||||||||

* Department of Genetics and Biotechnology, Faculty of Agricultural Sciences, Aarhus University, PO Box 50, Blichers Alle, 8830 Tjele, Denmark
Department of Large Animal Sciences, Faculty of Life Sciences, University of Copenhagen, Grønnegårdsvej 8, 1870 Frederiksberg, Denmark
1 Corresponding author: jan.lassen{at}agrsci.dk
| ABSTRACT |
|---|
|
|
|---|
Key Words: heterogeneous variance genetic interaction free stall tie stall
Over the last decades, a considerable structural development has occurred in Danish dairy cattle production, which has resulted in a change in the way cows are housed. A large number of free-stall barns have been built and today most cows are housed this way (an increase from 48 to 75%). The majority of cows entering these new barns are, however, sired by bulls that are selected based on daughter performance in tie stalls. Therefore, the EBV ranking of these bulls could change based on additional offspring in a different environment.
In this study, focus was on conformation and workability traits and their relation to 2 housing environments. It is of interest to investigate these traits because cows might develop differently in the 2 environments to function and work efficiently. This may be true for all traits but probably more pronounced for feet and leg traits because cows need to be active to function in a free stall.
Dairy cattle are kept in a range of different environments and many examples of genotype by environment interaction have been documented in the scientific literature. Genotype by environment interaction can come either from heterogeneous variances or from heterogeneous genetic correlations between traits across environments (Calus et al., 2002; Kolmodin et al., 2002; Carlén et al., 2005). The genotypes can be measured at different levels; for example, as certain breeds, groups of animals, individual animals, or genes. Likewise, the environment can be defined in many ways such as country, region, herd, cage, climatic descriptor, or housing system. In this study, the additive genetic value of the animal is the considered genotype, and the environment in question is the housing system.
A situation in which genotype by environment interaction occurs but is ignored could lead to biased breeding values and suboptimal selection decisions. One way to account for genotype by environment interaction in terms of heterogeneous variance is to preadjust phenotypically the data (Koots et al., 1994; Weigel and Lawlor, 1994). A way to account for genotype by environment interaction in terms of reranking is to conduct multivariate analyses of data in which the covariance between a trait measured in different environments is taken into account. In international genetic evaluations for bulls, genotype by environment interaction is accounted for by considering data from each country as different but correlated traits (Schaeffer, 1994). Bulls rank differently in each country because of across-country genetic correlations of around 0.6 to 0.9 depending on trait (Mark, 2004). Across-country genetic correlations less than unity are due, in part, to genotype by environment interactions (Mark et al., 2006, 2007).
In an analysis of Canadian Holstein data for 7 feet and leg traits measured in free stalls and tie stalls, genetic correlations significantly different from unity were reported on 6 of the traits (Fatehi et al., 2003). The study also reported that heritabilities tended to be greater in tie stalls than in free stalls. Burke and Funk (1993) reported significantly different genetic regression coefficients in American Holsteins for herd life on different conformation traits depending on housing category (tie-stall or free-stall systems). This indicates that genotype by environment interaction exists for conformation traits measured in different housing systems. In a study of longevity in German Holsteins the effect of housing system was included in the model (Buenger et al., 2001). The results showed that risk of culling was lower in free stalls than in tie stalls. The effect of straw bedding, however, was more pronounced.
Genotype by environment interaction in the form of reranking was found for SCC and mastitis across different environments of SCC in Swedish Holstein cows (Carlén et al., 2005). Based on principal component scores for fertility, substantial differences in genetic correlation between 305-d milk production and fertility were found in Dutch Holstein cows (Windig et al., 2006). Genotype by environment interaction for fertility traits was found in 4 Nordic Red breeds using reaction norm models (Kolmodin et al., 2002).
To date no studies have been conducted in terms of genotype by environment interaction in relation to housing systems in Danish Holsteins. Given the latest extreme structural development in dairy cattle production, this study is important to ensure that the genetic parameters used in the population corresponds to the actual production circumstances. The main aim of this study was to investigate the extent of genotype by housing type interaction for conformation and workability traits in Danish Holsteins and to discuss consequences for genetic evaluation and selection.
Information about housing systems for Danish dairy cattle herds from 1997 to 2002 was made available by the Danish Cattle Federation. Several free-stall barns was built during the last 10 yr, and the proportion of cows classified in free-stall barns increased from 48% in 1997 to 75% in 2002 (Figure 1
). These numbers reflect the proportion of cows in free-stall barns in the whole population. Conformation and workability records on first-parity Holstein cows obtained in the same period were extracted from the Danish Cattle database. Records of cows with missing sire information were omitted. Records in herds with <25 classified cows in the period were omitted as well as records made by classifiers with <100 registrations in the period. After this editing of the data, a total of 30,190 cows were available. Of these, 36% were registered in tie stalls (herds = 391) and 64% were registered in free stall barns (herds = 529). Not all cows had observations for all traits. Rear udder height and rear teat placement were only measured since 2000 and body condition since 2001 (Table 1
). The pedigree was traced as far back as possible and the pedigree file contained 169,041 animals.
|
|
|
![]() | (1) |
where h was a fixed effect of the ith herd (920 herds); ys was a fixed effect of the jth year and season of classification (4 seasons per year); dim was the days after calving and b1 a fixed regression coefficient (22
dim
658 d); age was the age of the cow when the record was taken and b2 a fixed regression coefficient (579
age
1,379 d); cy was a fixed effect of year within each classifier (5 classifiers in 6 yr = 30 cy effects; this effect was not included for the 2 workability traits because they were scored by the farmer); and a was a random additive genetic effect. The vector of additive genetic effects (a) was assumed to be N(0, A
2a), where A is the additive genetic relationship matrix; e was a random residual effect. The vector of residual effects (e) was assumed to be N(0, I
2e), and covar(a, e) = 0.
The bivariate model considered performance in tie stalls and free stalls as different but correlated traits and included the same explanatory variables as model [1]. The (co)variance structure among vectors of random variables in the bivariate model was
![]() |
Genetic parameters were estimated using average information (AI)-REML (Jensen et al., 1997) as implemented in the DMU software (Madsen and Jensen, 2008) for both the univariate and bivariate animal models.
Little evidence for genotype by housing interaction was found in this study in terms of reranking caused by genetic correlations different from unity (Table 3
). All genetic correlations were greater than 0.90, except for body width (0.87 ± 0.06). Genetic correlations for the 2 workability traits were 0.94 and 0.95 and were a little lower than for most of the conformation traits.
|
When editing the data, restrictions on number of observations per classifier (at least 100) and number of classifications per herd in the period (at least 25) were made. The restriction on number of observations per classifier removed very little of the data, so the main reason for data deletion was number of classifications per herd. Tie stalls also tend to be smaller than free stalls, so this restriction meant that the distribution of animals in each housing system was different than in the received data. The main reason for having this restriction was to obtain better estimates of the herd effect by avoiding small groups and thereby obtain lower standard errors on variance components. Using less stringent restrictions on number of classifications per herd to see if genetic parameters were biased was not performed in this study. Also, no restrictions were made on herds that changed from tie stalls to free stalls in the studied period. These herds still contributed with information to the data set according to their housing system code in the specific year.
Burke and Funk (1993) showed that cows with intermediate curvature of rear legs had longer herd life than cows with extreme curvature in both free-stall and tie-stall environments, although the absolute difference between optimum and extreme was greater in tie stalls than in free stalls. Burke and Funk (1993) also reported greater genetic regressions of herd life on nonlocomotion traits in tie-stall barns than in free-stall barns. The reason for this might be personnel (farmer) preference for different type traits in different environments.
The results in this study will have limited implications for routine evaluations of conformation and workability traits in Danish Holsteins. The traits measured in the 2 different housing systems can be considered the same trait in genetic evaluations because the genetic correlations between traits measured in different housing systems were very high and most were close to unity. Moreover, if the current trend continues, most cows in the future will be in free-stall barns; therefore, splitting the breeding values in 2 will be of decreased interest. However, heterogeneous genetic and residual variances in free and tie stalls were found for conformation and workability traits in Danish Holsteins. Heritabilities were on average 18% greater for traits measured in tie stalls than in free stalls. Therefore, heterogeneous variances could be accounted for in genetic evaluation of these traits. This was also the recommendation from a similar study made in Canadian Holsteins with focus on feet and leg traits (Fatehi et al., 2003). For 6 of 7 feet and legs traits, however, they found genetic correlations between registrations in free-stall and tie-stall environments that were significantly different from 1. A bivariate model to predict breeding values would improve the breeding value estimation if evidence of genetic correlations different from unity between traits measured in 2 environments were found. In this study, difference in terms of heterogeneous variance was found but it was insufficient to justify the implementation of a bivariate model for prediction of breeding values. Moreover, 2 breeding values for the same trait might lead to more confusion for the farmer when choosing bulls for mating.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication February 22, 2008. Accepted for publication July 1, 2008.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. V. Laursen, D. Boelling, and T. Mark Genetic parameters for claw and leg health, foot and leg conformation, and locomotion in Danish Holsteins J Dairy Sci, April 1, 2009; 92(4): 1770 - 1777. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |