|
|
||||||||
,1




* Agriculture and Agri-Food Canada, Dairy and Swine Research and Development Centre, Sherbrooke, Quebec, Canada, J1M 1Z3
Canadian Dairy Network, Guelph, Ontario, Canada, N1G 4T2
1 Corresponding author: sewalem{at}cdn.ca
| ABSTRACT |
|---|
|
|
|---|
Key Words: herd life longevity genetic evaluation multiple-trait animal model
Longevity or herd life is a highly desirable trait that considerably affects overall profitability (Congleton and King, 1984; Allaire and Gibson, 1992). With increased longevity, the mean production of the herd increases. Furthermore, the economic importance of herd life compared with milk production is considered higher than other nonproduction traits (van Arendonk, 1991; Allaire and Gibson, 1992; Dekkers, 1993).
Longevity can be measured in several ways (VanRaden and Klaaskate, 1993; Brotherstone et al., 1998; Vollema et al., 2000; Ducrocq, 2002; VanRaden et al., 2006) and genetic evaluation systems are not standardized across countries, which makes the comparison of sire rankings difficult. In Canada, the survival of cows in each of the first 3 lactations [survival through first lactation (S1), survival through second lactation (S2), and survival through third lactation (S3)] are recorded as a binary traits and evaluated with a multiple-trait linear animal model in which survival in each lactation is considered as a separate trait (Jairath et al., 1998). For the traits S1, S2, and S3 the cow must initiate the second, third, and fourth lactations, respectively. Hereafter, we refer to this as the old system of evaluation of herd life in Canada.
An international genetic evaluation for longevity using direct herd life (DHL) carried out by the Interbull Centre in February 2005 (International Bull Evaluation Service, 2005) showed that, for Canada, the correlations with other countries were higher than correlations among most other countries. The report also showed that, compared with other countries, the Canadian system results are in favorable EBV correlation with other countries. However, with the old system of genetic evaluation for herd life, a 2-yr opportunity window was used to determine if a subsequent calving occurred, thus delaying when Canadian bulls might get their first genetic evaluations for herd life. The outcome of this extended period, before newly proven sires have direct survival data in their EBV, was that newly proven Canadian sires did not receive an international evaluation for DHL as quickly as newly proven bulls in other countries. The objective of the present study was to improve the genetic evaluation of herd life by breaking the first lactation into 3 periods; from first calving to 120 DIM (S1a), from 120 to 240 DIM (S1b), and from 240 DIM to second calving (S1c), thus including direct survival data into a bulls genetic evaluation earlier. Hereafter, we refer to this as the new system of evaluation of herd life in Canada.
Data for genetic parameter estimations were obtained from lactation records extracted for the official Canadian February 2005 genetic evaluation of herd life in the Holstein breed. The extract was modified to contain all production traits such as breed-corrected average (BCA), mature-equivalent yields, 305-d yields, and total cumulative yields for milk, fat, and protein (in kg). The primary production information used was BCA of the cow; if BCA yields were missing, then mature-equivalent yields, 305-d yields, or total cumulative yields were used. In all cases, the yields were then converted to BCA units. First-lactation records with calving age from 18 to 40 mo were used.
Censored records included cows that were sold for dairy purposes, exported, or leased, as well as cows still present in the herd. A lifetime record was considered completed if the cow received a termination code, indicating that she was culled or had died for any reason. For cases in which cows completed a given lactation but failed to initiate the next lactation, cows were considered culled if the latest recorded test-day record was more than 6 mo before the end of the study. Herds with at least 5 calvings in a given quota year (July to June) were assumed to be on test during that year.
The methodology that was developed for computing genetic evaluations for herd life consisted of 3 major steps: 1) computation of genetic evaluations for DHL based on survival of daughters; 2) computation of genetic evaluations for indirect herd life (IHL) based on an index of EBV for conformation, fertility, and udder health traits; and 3) combining genetic evaluations for DHL and IHL into an overall genetic evaluation for herd life (HL). The present study describes changes made in the genetic evaluation model for DHL, in which the first lactation is divided into 3 periods of survival. The revised model for IHL includes other traits such as fertility and health traits. The method to combine DHL and IHL in the overall index (HL) was as previously described by Jairath et al. (1998).
| Trait Definition and Model for DHL |
|---|
|
|
|---|
The survival data were analyzed using a multiple-trait animal model in which the 5 survival traits were considered as different traits. The following model was used for genetic parameter estimation:
![]() |
where yijkm = observation for survival (0 or 1) in trait m (m = 1 to 5) on cow k that calved in herd-year i; hyim = fixed effect of herd-year of calving i for trait m; rhsjm = fixed effect of subclass j for registry status (2 classes) x herd size change (3 classes) x season of calving (4 classes) for trait m; (agek)m = linear and quadratic regressions of survival in trait m on age at first calving k; trait m on protkm by rhsjm subclasses, where protk is the normal rank of the cow for protein yield in first lactation within herd-year of first calving; (protk)m and (fatk)m = linear, quadratic, and cubic regressions of survival in trait m on the cows normal ranks for protein and fat yield in lactation 1 within herd-year of first calving; animalkm = random additive genetic effect of animal k for survival in trait m; and eijkm = random residual. Contemporary groups of herd-year of calving were formed in relation to the dairy quota year (July to June) rather than a calendar year because of the impact of the quota system in Canada on culling. Genetic parameters for Holsteins were estimated using the software package VCE-5 (Groeneveld, 1998).
Estimates of genetic parameters for S1a, S1b, S1c, S2, and S3 are presented in Table 1
. The genetic variances for the first 2 traits (S1a and S1b) were small compared with the other traits. This was expected because the phenotypic variation within the first 2 traits was very small. Generally, similar low heritability values of herd life were reported (Jairath et al., 1998; Sewalem et al., 2005; VanRaden, et al., 2006). Genetic correlations among the 5 traits were moderate to high ranging from 0.73 to 0.97 (Table 1
).
|
| Genetic Evaluation for DHL |
|---|
|
|
|---|
Let Si = the mean survival rate for trait i + sire ETA for trait i, and the cumulative survival rate for trait i is
Si and C6 = 0.
![]() |
where DHL is the direct herd life, Ni = mean calving interval in days for trait i (n = 0 for trait 1, 2, and 3), Di are the means days of production for cows culled in trait i, and K is a constant; that is, the expected the mean calving intervals and average days of production for cows culled after the fourth calving.
The correlations among the 5 traits and the correlations of each trait with combined DHL are presented in Table 2
. Correlations among EBV of the 5 traits were very high, ranging from 0.90 to 0.99. Correlations of EBV of the 5 traits to the overall DHL were all > 0.95. Parallel genetic evaluation runs of the new (5-trait) and old (3-trait) systems for genetic evaluation of herd life were performed using the August 2006 data set. The correlations of DHL between the old and new systems of genetic evaluations ranged from 0.89 to 0.98 depending on the number of daughters per sire. As the number of daughters per sire increased, the correlation increased.
|
| Genetic Evaluations of IHL |
|---|
|
|
|---|
| Overall Index of HL |
|---|
|
|
|---|
In conclusion, the new system of herd life evaluation with 5 traits has been used for national evaluation of herd life for all Canadian dairy breeds (Ayrshire, Brown Swiss, Canadienne, Guernsey, Holstein, Jersey, and Milking Shorthorn) since February 2006. A DHL EBV is obtained for each bull by 7 mo after first calving of their daughters instead of 24 mo in the previous old system of evaluation herd life. The correlation of EBV between the new system and old system of evaluation was high, ranging from 0.89 to 0.98.
Received for publication November 1, 2006. Accepted for publication December 20, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Holtsmark, B. Heringstad, and J. Odegard Predictive abilities of different statistical models for analysis of survival data in dairy cattle J Dairy Sci, November 1, 2009; 92(11): 5730 - 5738. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |