JDS
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J. Dairy Sci. 2008. 91:4424-4428. doi:10.3168/jds.2008-1116
© 2008 American Dairy Science Association ®

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Interpretive Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lassen, J.
Right arrow Articles by Mark, T.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Lassen, J.
Right arrow Articles by Mark, T.

Short Communication: Genotype by Housing Interaction for Conformation and Workability Traits in Danish Holsteins

J. Lassen*,1 and T. Mark{dagger}

* Department of Genetics and Biotechnology, Faculty of Agricultural Sciences, Aarhus University, PO Box 50, Blichers Alle, 8830 Tjele, Denmark
{dagger} 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
 TOP
 ABSTRACT
 ACKNOWLEDGEMENTS
 REFERENCES
 
A total of 30,190 first-parity Danish Holstein cows housed in free stalls or tie stalls were analyzed to quantify to what degree genotype by housing interaction existed for 21 conformation and 2 workability traits. Each trait measured in different housing systems was treated as 2 separate traits in a bivariate animal model. Genetic correlations between the 2 traits as well as differences in genetic and residual variance were used as measurements of whether or not genotype by housing interaction occurred. Genetic correlations were in general close to unity (>0.9), except for body width (0.87 ± 0.06). Genetic variance was greater (22% on average) for all traits except 2, whereas residual variance was greater for feet and leg traits and lower for body and udder traits when measured in tie-stall environments for nearly all traits. In total, this led to a greater heritability for traits measured in tie-stall barns. Conformation and workability traits measured in free stalls and tie stalls should not be considered as different traits when predicting breeding values, but heterogeneous variance could be accounted for.

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 1Go). 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 1Go). The pedigree was traced as far back as possible and the pedigree file contained 169,041 animals.


Figure 1
View larger version (12K):
[in this window]
[in a new window]

 
Figure 1. Trend in distribution of cows in tie and free stalls from 1997 to 2002.

 

View this table:
[in this window]
[in a new window]

 
Table 1. Number of records before and after data editing
 
A total of 23 traits were available in this study. Seven of the traits were related to the body of the cow, 5 were feet and leg traits, 9 were udder traits, and 2 were workability traits (i.e., milking speed and temperament). All traits except stature were measured on a scale from 1 to 9 with increments of 1 and with different optima (Table 2Go). Stature was measured in centimeters and had an optimum between 147 and 150 cm. Conformation traits were objectively measured by professional classifiers, whereas the 2 workability traits were reported by the farmer. The different optima were decided by the national breed associations. Some of the traits are measured today but are not a part of the present breeding goal. Breeding values for the traits are predicted and published.


View this table:
[in this window]
[in a new window]

 
Table 2. Trait optimum, mean, and SD from single-trait analysis1
 
Variance components were estimated using both a single-trait model (assuming performance in each housing system to be the same trait) and a bivariate model (assuming performance in each housing system to be different traits). The single-trait model was


Formula 1(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{sigma}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{sigma}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


Formula 1

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 3Go). 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.


View this table:
[in this window]
[in a new window]

 
Table 3. Estimates of heritability (h2), additive genetic variance (Va) and residual variance (Ve) for performance in free stall (free), tie stall (tie), and combined performance (all), and genetic correlation (rg) between performance in free and tie stalls for different traits with relative SE1
 
In general, heritabilities were lower (18% on average) in free-stall barns than in tie-stall barns, and estimates from univariate models were intermediate (Table 3Go). This was more pronounced for body and udder traits than for feet and legs traits. However, the heritability was greater in free-stall barns than in tie-stall barns for rear leg rear view, rear udder height, and body condition score. For rear leg rear view, one possible explanation could be that classifiers saw some cows walk in the free stalls during classification, and thereby had a better impression of the conformation of the rear legs. The lower heritability in free-stall barns comes both from lower additive genetic variance and greater residual variance in free-stall barns (Table 3Go). Standard errors for heritability and additive genetic variance were between 0.01 and 0.09 for all traits except for stature, for which the standard errors were between 0.24 and 0.41. Differences in additive genetic variance between free-stall and tie-stall barns may indicate genotype by environment interaction. Another explanation could be that the lower additive genetic variance was due to more cows having a wrong ID or incorrect pedigree information in free stalls compared with tie stalls. However, this has not been investigated further. The greater residual variance indicated that classifiers had more difficulty conducting accurate classification in free stalls than in tie stalls. Free-stall barns are typically larger and cows are not necessarily used to close contact with humans as occurs in a classification situation. There are no guidelines on whether cows should be fixed or not when classified in a free-stall environment, which could lead to inconsistencies in the classification and thereby greater residual variance.

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
 TOP
 ABSTRACT
 ACKNOWLEDGEMENTS
 REFERENCES
 
The authors would like to acknowledge Vibeke Duch-waider (LIFE, University of Copenhagen, Denmark) for preliminary work on the data and the Danish Agricultural Advisory Centre (Skejby, Denmark) for contributing data to this study.

Received for publication February 22, 2008. Accepted for publication July 1, 2008.


    REFERENCES
 TOP
 ABSTRACT
 ACKNOWLEDGEMENTS
 REFERENCES
 


Buenger, A., V. Ducrocq, and H. H. Swalwe. 2001. Analysis of survival in dairy cows with supplementary data on type scores and housing systems from a region of Northwest Germany. J. Dairy Sci. 84:1531–1541.[Abstract]

Burke, B. P., and D. A. Funk. 1993. Relationship of linear type traits and herd life under different management systems. J. Dairy Sci. 76:2772–2782.

Calus, M. P. L., A. F. Groen, and G. de Jong. 2002. Genotype x environment interaction for protein yield in Dutch dairy cattle as quantified by different models. J. Dairy Sci. 85:3115–3123.[Abstract/Free Full Text]

Carlén, E., K. Jansson, and E. Strandberg. 2005. Genotype by environment interaction for udder health traits in Swedish Holstein cows. Page 47 in Proc. 56th Annu. Mtg. Eur. Assoc. Anim. Prod., Uppsala, Sweden. Wageningen Acad. Publ., Wageningen, the Netherlands.

Fatehi, J., A. Stella, J. J. Shannon, and P. J. Boettcher. 2003. Genetic parameters for feet and leg traits evaluated in different environments. J. Dairy Sci. 86:661–666.[Abstract/Free Full Text]

Jensen, J., E. A. Mäntysaari, P. Madsen, and R. Thompson. 1997. Residual maximum likelihood estimation of (co)variance components in multivariate mixed linear models using average information. J. Ind. Soc. Agric. Stat. 49:215–236.

Kolmodin, R., E. Strandberg, P. Madsen, J. Jensen, and H. Jorjani.2002. Genotype by environment interaction in Nordic dairy cattle studied using reaction norms. Acta Agric. Scand. A 52:11–24.

Koots, K. R., K. M. Wade, B. W. Kennedy, J. C. M. Dekkers, G. C. Smith, and E. B. Burnside. 1994. Method and effect of adjustment for heterogeneous variance of Holstein conformation traits. J. Dairy Sci. 77:294–302.[Abstract]

Madsen, P., and J. Jensen. 2008. DMU: A user’s guide. A package for analysing multivariate mixed models. Version 6, release 4.7. DJF, Foulum, Denmark. http://www.dmu.agrsci.dk/

Mark, T. 2004. Applied genetic evaluations for production and functional traits in dairy cattle. J. Dairy Sci. 87:2641–2652.[Abstract/Free Full Text]

Mark, T., W. F. Fikse, P. G. Sullivan, and P. M. VanRaden. 2007. Prediction of genetic correlations and international breeding values for missing traits. J. Dairy Sci. 90:4805–4813.[Abstract/Free Full Text]

Mark, T., P. G. Sullivan, W. F. Fikse, and P. M. VanRaden. 2006. Prior genetic correlations and non-measured traits. Interbull Bull. 35:72–75.

Schaeffer, L. R. 1994. Multiple-country comparisons of dairy sires. J. Dairy Sci. 77:2671–2678.[Abstract]

Weigel, K. A., and T. J. Lawlor. 1994. Adjustment for heterogeneous variance in genetic evaluations for conformation of United States Holsteins. J. Dairy Sci. 77:1691–1701.[Abstract]

Windig, J. J., M. P. L. Calus, B. Beerda, and R. F. Veerkamp. 2006. Genetic correlations between production and health and fertility dependent on herd environment. J. Dairy Sci. 89:1765–1775.[Abstract/Free Full Text]


This article has been cited by other articles:


Home page
J DAIRY SCIHome page
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]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Interpretive Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lassen, J.
Right arrow Articles by Mark, T.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Lassen, J.
Right arrow Articles by Mark, T.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS