J. Dairy Sci. 87:1457-1468
© American Dairy Science Association, 2004.
Genetic Correlations Among Production, Body Size, Udder, and Productive Life Traits Over Time in Holsteins
S. Tsuruta1,
I. Misztal1 and
T. J. Lawlor2
1 Animal and Dairy Science Department, University of Georgia, Athens 30602
2 Holstein Association USA Inc., Brattleboro, VT 05301
Corresponding author: S. Tsuruta; E-mail: shogo{at}uga.edu.
 |
ABSTRACT
|
|---|
Genetic correlations among milk, fat, and protein yields; body size composite (BSC); udder composite (UDC); and productive life (PL) in Holsteins were investigated over time. The data set contained 25,280 records of cows born in Wisconsin between 1979 and 1993. The multiple trait random regression (MT-RR) animal model included registration status, herd-year, age group, and stage of lactation as fixed effects; additive genetic effects with random regressions (RR) on year of birth using the first-order Legendre polynomial; and residual effects. Heterogeneous residual variances were considered in the model. Estimates of variance components and genetic correlations among traits from MT-RR were compared with those estimated with a multiple trait interval (MT-I) model, which assumed that every 3-yr interval was a separate trait and included the same effects as in the MT-RR model except for the RR. Genetic correlations estimated with MT-RR and MT-I models over time among all traits were compared with correlations among breeding values predicted with the single trait (ST) model without RR. Correlations among breeding values predicted with MT-RR, ST, and MT models were also calculated.
Additive genetic and residual variances for all traits except PL increased over time; those for PL were constant. As a result, heritability estimates had no significant changes during the 15 yr. Genetic correlations of PL with milk, fat, protein, and BSC declined to zero or negative; those with UDC remained positive. Correlations among breeding values predicted with ST, MT, and MT-RR models were relatively high for all traits except PL.
Genetic correlations between PL and other traits varied over time, with some correlations changing sign. For accurate indirect prediction of PL from other traits, the genetic correlations among the traits need to be re-estimated periodically.
Key Words: genetic correlation random regression Holstein
Abbreviation key: BSC = body size composite, HPD = highest posterior density interval, MT = multiple trait, MT-I = MT interval, PL = productive life, RR = random regression, ST = single trait, UDC = udder composite
 |
INTRODUCTION
|
|---|
The genetic evaluation of type traits in Holsteins has been conducted with a multiple trait (MT) model by Holstein Association USA, Inc. (Brattleboro, VT). The MT model assumes that additive genetic and environmental variances and genetic and environmental correlations among traits are constant over time. Indirect prediction of productive life (PL) from type traits (Weigel et al., 1998; VanRaden, 2001) also uses the same assumption. However, this assumption may not be valid. Additive genetic variance could be changed not only by selection but also by other factors such as migration, segregation, mutation, inbreeding, and re-definition of the recording system. Dairy producers decisions on culling cows are based on various reasons, such as milk, reproduction problems, mastitis, and labor costs. All of these factors can change over time. Environmental or residual variances could change because of economic and management reasons. Selection on one trait may impede the genetic improvement of other traits because of linkage disequilibrium. Weigel et al. (1998) showed weights and maximum reliability for direct and indirect PL. If genetic parameters change over time, those weights for PL may also vary over time. Some linear type traits may have an intermediate optimum. Selection on a linear index assumes that the extremes are best.
Lawlor et al. (2002) presented changes in genetic parameters of productions traits, linear type traits, and PL over time using a random regression (RR) model. Tsuruta et al. (2003) also investigated changes in genetic parameters of final score by RR models with regression on year of birth. As verified by simulation, such estimation can be successful; however, changes in genetic parameters must be gradual over time, and residual variances need to be modeled as heterogeneous.
The objective of this study was to estimate changes in genetic parameters for PL, production, and selected type traits over time.
 |
MATERIALS AND METHODS
|
|---|
Data
Production data (305-d milk, fat, and protein yields) and PL were provided by USDA-AIPL, and linear type data were obtained from Holstein Association USA, Inc. The final data set included 25,280 records of Wisconsin Holstein cows born over 15 yr between 1979 and 1993, with a total of 40,838 in the pedigree that included 3 generations preceding the recorded animals. The numbers of records; means; standard deviations for milk, fat, and protein yields; PL; body size composite (BSC); and udder composite (UDC) are shown in Table 1
. Milk, fat, and protein yields for 305 d were pre-adjusted by age at calving to mature equivalent records. Productive life was defined as months in milk at age 84 mo as described by VanRaden and Klaaskate (1993) but was not adjusted for milk yield. According to formulas in Sire Summaries (2003), BSC and UDC were defined with weighted linear type traits as follows:
View this table:
[in this window]
[in a new window]
|
Table 1. Data description for milk, fat, and protein yields; productive life (PL); body size composite (BSC); and udder composite (UDC).
|
|
Multiple-Trait Animal Model with RR
Multiple trait RR (MT-RR) models were defined as follows.
For milk, fat, and protein yields:
for BSC and UDC,
for PL,
where yijklmn = observations for registration status i, management group j, age group k, month of calving l, stage of lactation m, and cow n; regi = effects of registration status i (registered or grade); hyj = effects of herd-year j; agk = effects of age group k; mocl = month of calving l, slm = effects of stage of lactation x year at classification m; anp = RR coefficient p on year of birth for additive genetic effect on animal n; znp = Legendre polynomials; and eijklmn = residual effects. The first-order Legendre polynomial standardized to the range of 1 to +1 for the years from 1979 to 1993 was applied to the MT-RR model. In matrix notation, the model is
 |
where y = a vector of final scores, ß = a vector of fixed effects, a = a vector of intercepts and linear RR coefficients for additive genetic effects, e = a vector of residual effects, X = an incidence matrix for fixed effects, and Z = a matrix of covariates of Legendre polynomials for additive genetic effects. The variances are defined as
where G = an additive genetic covariance matrix with size 12 [6 traits x (intercept + linear RR) for linear RR],
= the direct matrix product, Ap = the additive genetic relationship matrix for p animals, R = a residual variance matrix for 6 traits with
for interval I, and Ii = an identity matrix for interval i. Heterogeneous residual variances were estimated in 5 intervals of 3 yr.
MT Interval Animal Model
A multiple trait interval (MT-I) model was used to compare the genetic parameters with those from the MT-RR model. The MT-I model treated 3-yr intervals as separate traits for a total of 30 traits (5 intervals x 6 original traits). Effects included in this model were the same as those in the MT-RR model except additive genetic effects did not include RR. Therefore, the additive genetic covariance matrix consisted of 30 x 30 parameters including covariances among the intervals, and the residual variance was the same size without the covariances.
Single-Trait Animal Model Without RR
A single trait (ST) model without RR was used to predict breeding values. Effects in the model were the same as those in the MT-RR model without RR. In the ST model, breeding values were predicted ignoring correlations among 6 traits. Afterward, correlations among all traits were calculated using predicted breeding values on all animals to compare with genetic correlations estimated with the MT-RR model.
MT Animal Model Without RR
A MT model without RR was also used to predict breeding values. The MT model assumes constant correlations among traits over time. Effects in the model were the same as those in the MT-RR model without RR. Correlations of these predicted breeding values with those from the ST model and from the MT-RR model were calculated for bulls with
10 daughters.
Computer Programs
Variance components were estimated with GIBBS3F90 (Misztal, 2003), which is the Gibbs sampling program with heterogeneous residual variances. After 20,000 Gibbs samples were discarded as burn-in, 80,000 samples were used to calculate posterior means and standard deviations for variance components, heritability, and genetic correlations. Breeding values were calculated with BLUP90IOD (Tsuruta et al., 2001), which uses the pre-conditioned conjugate gradient algorithm with iteration on data.
 |
RESULTS AND DISCUSSION
|
|---|
Data description
Means and standard deviations for 305-d milk, fat, and protein yields; PL; BSC; and UDC are presented in Table 1
. Means and standard deviations for milk, fat, and protein increased gradually over the years. Means for PL decreased by 4.6 mo in 15 yr, but the standard deviations were relatively stable; means and standard deviations for BSC and UDC fluctuated over the years.
Variance Components and Heritability
Minimum and maximum additive genetic variances and residual variances in 15 yr for milk, fat, and protein yields; PL; BSC; and UDC estimated with the MT-RR model are presented in Table 2
. That table also shows minimum and maximum 1.96 x posterior standard deviations (at 95%) in 15 yr for all traits. The ranges for those variances were large, indicating heterogeneity over the years. However, the ranges for heritability estimates were relatively small except for PL. Figure 1
(af) demonstrates posterior means and highest posterior density intervals (HPD) at 95% for additive genetic variances for milk (a), fat (b), protein (c), PL (d), BSC (e), and UDC (f) with the MT-RR model and with the MT-I model. Those means for all traits from the MT-I model were approximately analogous to those from the MT-RR model, although some of the additive genetic variances were out of the HPD.
View this table:
[in this window]
[in a new window]
|
Table 2. Minimum and maximum additive genetic and residual variances and heritability (1.96 x SD) for milk, fat, and protein yields; productive life (PL); body size composite (BSC); and udder composite (UDC) in 15 yr estimated with a multiple trait random regression model.
|
|

View larger version (33K):
[in this window]
[in a new window]
|
Figure 1. Posterior means and highest posterior density intervals at 95% (upper and lower) for additive genetic variances using multiple trait random regression and posterior means using multiple trait interval (MT-I) models: A) milk, B) fat, C) protein, D) productive life, E) body size composite, and F) udder composite. Legend: mean (), lower limit (), upper limit (- - -), and MT-I results ( ).
|
|
Additive genetic variances for all traits except for PL tended to increase over time; those for PL were nearly constant with larger HPD, indicating that increases in additive genetic variances could be due to scaling effects. As means of production traits increased, those variances increased, but mean and variances of PL did not change. Posterior means and HPD for residual variances for all traits are presented in Figure 2
(af). Residual variances increased dramatically for milk and fat and increased slightly for protein and BSC. Residual and additive genetic variances for PL were constant. The rise in UDC from 1983 to 1986 might have been due to definition changes for udder traits.

View larger version (30K):
[in this window]
[in a new window]
|
Figure 2. Posterior means and highest posterior density intervals at 95% (upper and lower) for residual variances using multiple trait random regression and posterior means using multiple trait interval (MT-I) models: A) milk, B) fat, C) protein, D) productive life, E) body size composite, and F) udder composite. Legend: mean (), lower limit (), upper limit (- - -), and MT-I results ( ).
|
|
Increases in additive genetic and residual variance for milk, fat, and protein yields were consistent with increases in phenotypic standard deviations shown in Table 1
. The MT-I model clearly showed heterogeneity of the residual variance, suggesting that it was essential to consider heterogeneous residual variances in the MT-RR model. Posterior means for heritability for all traits from the MT-I model and the MT-RR model are presented in Figure 3
(af). Because changes in additive genetic and residual variances were similar in direction, heritability estimates from the MT-I model for all traits were relatively constant over time and within HPD from the MT-RR model. However, those from the MT-RR model for type traits were underestimated in comparison with those from the MT-I model.

View larger version (30K):
[in this window]
[in a new window]
|
Figure 3. Posterior means and highest posterior density intervals at 95% (upper and lower) for heritability using multiple trait random regression and posterior means using multiple trait interval (MT-I) models: A) milk, B) fat, C) protein, D) productive life, E) body size composite, and F) udder composite. Legend: mean (), lower limit (), upper limit (- - -), and MT-I results ( ).
|
|
Genetic Correlations
Minimum and maximum genetic correlations with 1.96 x standard deviations in 15 yr among all traits are presented in Table 3
. Genetic correlations among milk, fat, and protein yields were highly or moderately positive. Genetic correlations of production traits with BSC were low and positive over the years; those with UDC were low and negative. Genetic correlations of PL with production traits and BSC were low with large standard deviations, ranging from negative to positive. However, those with UDC were moderately positive.
View this table:
[in this window]
[in a new window]
|
Table 3. Minimum and maximum genetic correlations (1.96 x SD) among milk, fat and protein yields, productive life (PL), body size composite (BSC) and udder composite (UDC) in 15 years estimated with multiple trait random regression model.
|
|
Genetic correlations and their HPD at 95% among milk, fat, protein, PL, BSC, and UDC for 15 yr using MT-RR and MT-I models are presented in Figure 4 (ao)
. Genetic correlations between milk and fat were positive and slightly decreased from 0.7 to 0.5 over time (Figure 4a
). Genetic correlations of milk with protein were around 0.8 and constant (Figure 4b
). Miller et al. (1967) estimated high genetic correlations (from 0.54 to 0.77) between herd life and milk production. Genetic correlations of milk with PL decreased from positive to zero over the years with large HPD, demonstrating dramatic changes over time and implying that higher producing cows used to live longer (Figure 4c
).

View larger version (31K):
[in this window]
[in a new window]
|
Figure 4. (af)Posterior means and highest posterior density intervals at 95% (upper and lower) for genetic correlations between traits using multiple trait random regression and posterior means using multiple trait interval (MT-I) models and correlations between predicted breeding values for those traits with a single-trait (ST) model: A) milk:fat, B) milk:protein, C) milk:productive life, D) milk:body size composite, E) milk:udder composite, F) fat:protein, G) fat:productive life, H) fat:body size composite, I) fat:udder composite, J) protein:productive life, K) protein:body size composite, L) protein:udder composite, M) productive life:body size composite, N) productive life:udder composite, and O) body size composite:udder composite. Legend: mean (), lower limit (), upper limit (- - -), MT-I results ( ), and ST results ( ).
|
|
Genetic correlations of milk with BSC were positive but low as shown in Table 3
, indicating that larger cows might produce more milk, although the trend was not strong (Figure 4d
). Low negative genetic correlations of milk with UDC imply that higher UDC does not reflect more milk production (Figure 4e
). Genetic correlations of fat with protein were high and decreased over the years (Figure 4f
). Genetic correlations of fat with PL also declined and approached zero or turned negative, but those with type traits were low and increased (Figure 4gi
). Genetic correlations of protein with PL decreased, and those with type traits were low and slightly increased or remained constant (Figure 4jl
). Genetic correlations of PL with BSC were low and decreased from positive to negative as did those between fat and PL (Figure 4m
). Hansen et al. (1999) reported that larger cows had shorter PL in selected lines for large vs. small body size. Our results in the latter years showed negative genetic correlations.

View larger version (32K):
[in this window]
[in a new window]
|
Figure 4. (gl)Posterior means and highest posterior density intervals at 95% (upper and lower) for genetic correlations between traits using multiple trait random regression and posterior means using multiple trait interval (MT-I) models and correlations between predicted breeding values for those traits with a single-trait (ST) model: A) milk:fat, B) milk:protein, C) milk:productive life, D) milk:body size composite, E) milk:udder composite, F) fat:protein, G) fat:productive life, H) fat:body size composite, I) fat:udder composite, J) protein:productive life, K) protein:body size composite, L) protein:udder composite, M) productive life:body size composite, N) productive life:udder composite, and O) body size composite:udder composite. Legend: mean (), lower limit (), upper limit (- - -), MT-I results ( ), and ST results ( ).
|
|

View larger version (17K):
[in this window]
[in a new window]
|
Figure 4. (mo)Posterior means and highest posterior density intervals at 95% (upper and lower) for genetic correlations between traits using multiple trait random regression and posterior means using multiple trait interval (MT-I) models and correlations between predicted breeding values for those traits with a single-trait (ST) model: A) milk:fat, B) milk:protein, C) milk:productive life, D) milk:body size composite, E) milk:udder composite, F) fat:protein, G) fat:productive life, H) fat:body size composite, I) fat:udder composite, J) protein:productive life, K) protein:body size composite, L) protein:udder composite, M) productive life:body size composite, N) productive life:udder composite, and O) body size composite:udder composite. Legend: mean (), lower limit (), upper limit (- - -), MT-I results ( ), and ST results ( ).
|
|
Genetic correlations between PL and UDC were low but positive (Figure 4n
), suggesting that better udder quality may be required to keep producing milk longer. However, because of low negative correlations between milk and UDC, cows may not be able to produce more milk in their lifetime. A producers main reasons to cull cows are changing over time. In the past, milk yield was influential, but now reproduction and mastitis are more important. These changes in trait definition for PL must have occurred gradually over the years. Vollema and Groen (1997) concluded that genetic correlations between conformation and longevity traits differed among years of birth and that genetic correlation between herd life and udder depth was positive. As they indicated, genetic relationship between 2 traits could be nonlinear. In addition, Larroque and Ducrocq (2001) and Pérez-Cabal and Alenda (2002) reported nonlinearity of relationships between type traits and PL or longevity in Holsteins; Fuerst-Waltl et al. (1998) found nonlinear genetic relationships between milk and type traits. Productive life is a trait that strongly deviates from the normal distribution, and some of the nonlinearity could be due to the use of linear methodology, which assumes a normal distribution. Pasman and Reinhardt (1999) also presented similar genetic correlations when functional (milk-adjusted) PL was obtained by survival analysis. Their correlations of PL with body size and udder traits were 0.17 and 0.33, respectively.
Correlation between breeding values for 2 traits may be a good indicator of the genetic correlation estimated directly from additive genetic (co)variances and, thus, may be used for validation. Figure 4
also shows correlations between breeding values predicted with the ST model for the 6 traits over the years. These correlations were in agreement with those from the MT-RR model.
Correlation Among Predicted Breeding Values
Correlations among predicted breeding values of milk, fat, protein, PL, BSC, and UDC for 435 bulls with
10 daughters using ST and MT models and the MT-RR model in 1992, 1986, and 1980 are presented in Table 4
. Breeding values predicted with the MT-RR model were compared with those with the ST model, which ignores genetic correlations among traits, and the MT model, which assumes that genetic correlations were constant over time. Overall, correlations among breeding values for all traits with different models were high (>0.9) except correlations for PL (0.4 to 0.9). Those low correlations might have been due to large standard deviations of those genetic parameters estimated with the MT-RR model and changes of genetic correlations with other traits. Table 4
also presents correlations between additive genetic effects in 1992, 1986, and 1980 derived from estimates of variance components. Those genetic correlations for PL were relatively low (0.22 between in 1992 and 1980 and 0.74 between in 1992 and 1986, respectively), but those for other traits were from 0.62 to 0.72 between 1992 and 1980 and >0.9 between 1992 and 1986, respectively, which suggests that although other traits changed moderately between 1980 and 1992, PL became a different trait. Estimates from RR models may be inflated particularly at the extremes, and it is possible that the actual correlations are somewhat higher.
View this table:
[in this window]
[in a new window]
|
Table 4. Correlation among predicted breeding values for 435 bulls with 10 daughters using single-trait (ST), multiple-trait (MT), and multiple trait random-regression (MT-RR) models for milk, fat, and protein yields; productive life (PL); body size composite (BSC); and udder composite (UDC) and estimated genetic correlation between additive genetic effects in 1992 and 1986 and in 1992 and 1980 using the MT-RR model.
|
|
Comparison of Models
One of the most important questions in this study is whether the MT-RR can provide reliable estimates. Shapes of correlations in this study, where RR were limited to linear, were constrained to approximately linear changes in correlations and quadratic changes in additive variances. Constraints in the MT-I model were smaller, but at higher computing costs. The estimates by MT-RR and MT-I models were mostly similar, although some differences were outside the 95% HPD interval. Also, the comparison of breeding values under MT-RR, MI, and ST models were in agreement. Thus, the MT-RR model with linear regression can be a reasonable tool for looking at genetic changes over time if such changes are gradual. Modeling of changes over time with higher flexibility can be accomplished by the use of nonlinear RR. With higher order polynomials, computing becomes more expensive, and modeling artifacts may be stronger (Tsuruta et al., 2003).
 |
CONCLUSIONS
|
|---|
Genetic correlations among traits can increase or decrease over time and even change sign. A RR model can estimate such changes if they are gradual. A producers culling decisions may be the primary reason for the change of genetic correlations between PL and some linear type traits. Periodical re-estimation of genetic parameters would be required for more accurate indirect prediction of PL.
 |
ACKNOWLEDGEMENTS
|
|---|
Financial support from Holstein Association USA Inc. is greatly appreciated. We are grateful to George Wiggans of USDA for providing production and PL data. We are thankful to Jan-Thijs van Kaam for his valuable comments.
Received for publication December 3, 2003.
Accepted for publication January 12, 2004.
 |
REFERENCES
|
|---|
Fuerst-Waltl, B., J. Sölkner, A. Essl, I. Hoeschele, and C. Fuerst. 1998. Non-linearity in the genetic relationship between milk yield and type traits in Holstein cattle. Livest. Prod. Sci. 57:4147.
Hansen, L. B., J. B. Cole, G. D. Marx, and A. J. Seykora. 1999. Productive life and reasons for disposal of Holstein cows selected for large versus small body size. J. Dairy Sci. 82:795801.[Abstract]
Larroque, H., and V. Ducrocq. 2001. Relationships between type and longevity in the Holstein breed. Genet. Sel. Evol. 33:3959.[Medline]
Lawlor, T. J., S. Tsuruta, L. Klei, and I. Misztal. 2002. Use of a random regression model to investigate changes in genetic parameters over time. Proc. 7th WCGALP, Montpellier, France. CD-ROM Communication 17:06.
Miller, P., L. D. Van Vleck, and C. R. Henderson. 1967. Relationships among herd life, milk production, and calving interval. J. Dairy Sci. 50:12831287.[Abstract/Free Full Text]
Misztal, I. 2003. BLUPF90 Manual. [Online]. Available at http://nce.ads.uga.edu/~ignacy/newprograms.html. Accessed Nov. 17, 2003.
Pasman, E., and F. Reinhardt. 1999. Genetic relationship between type composites and length of productive life of Black-and White Holstein cattle in Germany. Interbull Bull. 21:117121.
Pérez-Cabal, M. A., and R. Alenda. 2002. Genetic relationships between lifetime profit and type traits in Spanish Holstein cows. J. Dairy Sci. 85:34803491.[Abstract/Free Full Text]
Sire Summaries. 2003. Linear Type Evaluations. Holstein Association USA. Brattleboro, VT.
Tsuruta, S., I. Misztal, and I. Stranden. 2001. Use of the preconditioned conjugate gradient algorithm as a generic solver for mixed model-equations in animal breeding applications. J. Anim. Sci. 79:11661172.[Abstract/Free Full Text]
Tsuruta, S., I. Misztal, T. J. Lawlor, and L. Klei. Modeling final scores in US Holsteins as a function of year of classification using a random regression model. Livest. Prod. Sci. (accepted)
VanRaden, P. M. 2001. Methods to combine estimated breeding values obtained from separate sources. J. Dairy Sci. 34(E. Suppl.):E47E55.
VanRaden, P. M., and E. J. H. Klaaskate. 1993. Genetic evaluation of length of productive life including predicted longevity of life cows. J. Dairy Sci. 76:27582764.[Abstract]
Vollema, A. R., and A. F. Groen. 1997. Genetic correlations between longevity and conformation traits in an upgrading dairy cattle population. J. Dairy Sci. 80:30063014.[Abstract]
Weigel, K. A., T. J. Lawlor, Jr., P. M. VanRaden, and G. R. Wiggans. 1998. Use of linear type and production data to supplement early predicted transmitting abilities for productive life. J. Dairy Sci. 81:20402044.[Abstract]
This article has been cited by other articles:

|
 |

|
 |
 
H. D. Norman, J. L. Hutchison, J. R. Wright, M. T. Kuhn, and T. J. Lawlor
Selection on Yield and Fitness Traits When Culling Holsteins During the First Three Lactations
J Dairy Sci,
February 1, 2007;
90(2):
1008 - 1020.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
E. Hare, H. D. Norman, and J. R. Wright
Survival rates and productive herd life of dairy cattle in the United States.
J Dairy Sci,
September 1, 2006;
89(9):
3713 - 3720.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. M. VanRaden, C. M. B. Dematawewa, R. E. Pearson, and M. E. Tooker
Productive life including all lactations and longer lactations with diminishing credits.
J Dairy Sci,
August 1, 2006;
89(8):
3213 - 3220.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
S. Tsuruta, I. Misztal, and T. J. Lawlor
Changing Definition of Productive Life in US Holsteins: Effect on Genetic Correlations
J Dairy Sci,
March 1, 2005;
88(3):
1156 - 1165.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. M. VanRaden
Invited Review: Selection on Net Merit to Improve Lifetime Profit
J Dairy Sci,
October 1, 2004;
87(10):
3125 - 3131.
[Abstract]
[Full Text]
[PDF]
|
 |
|