|
|
||||||||
Animal Improvement Programs Laboratory, Agricultural Research Service, USDA, Beltsville, MD 20705-2350
1 Corresponding author: paul{at}aipl.arsusda.gov
| ABSTRACT |
|---|
|
|
|---|
Key Words: genetic evaluation multibreed cross-breeding
| INTRODUCTION |
|---|
|
|
|---|
Crossbreeding is of increasing interest to dairy producers and dairy geneticists (McAllister, 2002; Weigel and Barlass, 2003; Heins et al., 2006). Cole et al. (2005) included crossbred and purebred Brown Swiss and Holsteins in US evaluations for calving ease. The number of first-generation (F1) crossbred dairy cows with usable yield records was about 10,000 in 2001, the latest birth year with complete data. This exceeds the numbers of purebred Brown Swiss, Guernsey, Ayrshire, or Milking Shorthorn cows. Holsteins became popular in many countries because of superior milk production, but some crossbreds have economic merit that is comparable with purebred Holsteins and may exceed Holstein merit if calving ease, calf livability, cow fertility, and cheese yield pricing are considered.
Inclusion of data from crossbred animals can lead to more reliable evaluations of purebred relatives, more accurate comparisons of genetic merit among all potential mates, and improved breeding programs that identify the best gene combinations. Goals of this research were to compare methods for evaluating mixed-breed populations and then to apply the best methods for routine evaluation of US dairy cattle.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Pedigrees for 46,603,162 dairy cattle were traced to the earliest ancestors recorded electronically, with a lower birth year limit of 1950 because earlier ancestors were not stored. Most animals (99%) had ancestors of only 1 breed, but 431,000 had ancestors of >1 breed. Of those, >350,000 had breed compositions with <94% of 1 breed and >6% of another breed because the cross-breeding occurred within the most recent 4 generations of the pedigree. Beginning in November 2005, the percentage of primary breed was reported for bulls and cows with pedigrees that contain more than 1 breed.
Breed composition of the cows with first lactations in 2004 included 90.9% Holsteins, 6.2% Jerseys, 0.8% Brown Swiss, 0.4% Guernseys, 0.3% Ayrshires, <0.1% Milking Shorthorns, 1.2% F1 crossbreds (coefficients of heterosis >50%), and 0.3% backcross cows (coefficients of heterosis >25%). Nearly all F1 cows had Holstein as one parent breed, and contributions from the other breeds were proportional to population size as reported by VanRaden and Sanders (2003). The number of F1 crossbreds doubled in the last 3 yr. For bulls born since 1997, only 4 Jerseys and 1 Brown Swiss had >25 cross-bred daughters, and each of these bulls had >200 pure-bred daughters. More recently, semen from Scandinavian Red and French breeds was imported and the resulting daughters are nearly all F1 crossbreds. Since 1987, over 5,000 herds had at least 1 crossbred cow, and currently 1,377 herds were coded as mixed-breed herds containing >25% crossbreds or cows of a different breed.
Unknown-parent groups in the animal model were separated by breed, pedigree path (dams of cows, sires of cows, and parents of sires), national origin (US or foreign), and birth year. Groups were formed when they included at least 500 animals within a time period and at least 2,000 animals across all years. The grouping pattern was similar to that for Dutch evaluations (NRS, 2005), except that they required only 40 animals per group. Larger numbers are needed for traits with lower heritability. Crossbred ancestors with no records and only 1 progeny were kept in the relationship matrix and treated as known so that the system of equations could link to animals with records back to purebred ancestor groups.
Heterogeneous variance adjustments (Wiggans and VanRaden, 1991) were modified for all-breed analysis of production traits and DPR. For mixed-breed herds, variance within herd would be overestimated if no account were taken of breed differences. Variance adjustments for milk, fat, and protein were previously based on ratios of milk variances, but variances of fat yield were used in the all-breed analysis. Variance adjustments were not used for all-breed PL and SCS evaluations because they had not been used previously in official within-breed evaluations.
Means, sums of squared deviations, and degrees of freedom were accumulated separately within herd, year, and breed. Those variance estimates within herd, year, and breed were combined with estimates within region, year, and breed to produce final variance adjustment factors. The regional estimate acted as a prior and received credit equal to 20 degrees of freedom. Data for other breeds were adjusted to make genetic variance equal to Holstein base cows. Variance adjustment factors are reported in Table 1
.
|
Additional age-parity-region-time factors were included in the animal models to account for gradual changes that might occur after the multiplicative preadjustments for age-parity were developed in 1995. These were estimated separately in the within-breed animal model from the data for each breed, but were estimated as uniform effects across breeds in the all-breed model. Recent age effects indicated that cows of all breeds are more productive at early ages relative to mature ages when compared with estimates from past decades, continuing the trend reported by Norman et al. (1995) for Holsteins. Multiplicative preadjustments may need to be updated in the future if relative maturity rates of breeds change.
Research to simultaneously estimate breed-age-parity-region-time effects in the all-breed model was abandoned because of very poor convergence or divergence and after learning that the Netherlands had attempted similar research, also with poor results (G. de Jong, NRS, Arnhem, the Netherlands; personal communication, 2006). Differences among breeds in recent residual age-parity effects were large only when comparing SCS of Holsteins to other breeds. These differences were corrected by applying multiplicative factors of 1.00, 1.00, 0.98, 0.96, and 0.94 for parities 1 to 5 of Holstein SCS data after the previous official age-parity factors were already applied. Simultaneous estimation might be possible for the largest breeds but not for all breeds.
Management groups in the within-breed evaluation were separate for registered and grade Holsteins if at least 5 cows of each type were present, whereas cows within the other breeds were grouped together regardless of registry status. In the all-breed evaluation, crossbreds were grouped together with registered or grade cows to allow estimation of breed differences. Crossbred cows sired by Holstein bulls were treated as grades, but all cows sired by bulls of other breeds were treated as purebreds and grouped with purebred cows. Management groups for herds that maintain separate herd codes for cows of different breeds in theory could be combined if the owner name and ZIP code match (Garcia-Peniche et al., 2005), but groups in the present study were combined only if herd codes matched. For within-breed evaluations, heritability of yield traits for Jerseys and Brown Swiss was higher (0.35) than for other breeds (0.30). For the all-breed model, the higher heritability for daughters of Jersey and Brown Swiss sires was accounted for by adjusting their lactation-length weights.
Convergence of solutions was tested by comparing results after 300 rounds of iteration to results after additional rounds. Priors for unknown-parent groups were set to 0 initially, and group solutions after 300 rounds were used as priors for remaining rounds. Final PTA from the all-breed system was compared with the official USDA within-breed PTA from August 2005. The comparisons were not exact because the all-breed analysis included about 2 mo more data for yield traits and SCS, and about 1 mo more data for PL, than did August 2005 evaluations.
Separate evaluations that included information from crossbred cows based on sire breed were also tested and compared with the within-breed evaluations. In the within-breed system, information from crossbred cows was used only if it was recorded in a grading-up program of a breed association. In the sire-breed evaluation, information from all crossbred daughters was included, but their pedigrees were truncated at the nearest purebred ancestor of another breed; more distant ancestors of other breeds were treated as unknown because the data file based on sire breed included no information for them. The minimum number of unknown dams per birth-year group was reduced to 150; separate groups for less numerous breeds were not estimated but instead were assumed to equal the primary breed.
Evaluations from an all-breed model can be reported with different genetic base options and including or excluding heterosis. An all-breed base was calculated using the mean of all cows born in 2000. Within-breed bases were calculated from the PTA of cows with coefficients of heterosis of <50% (i.e., F1 and backcross cows were not included). The PTA for each breed was adjusted to the within-breed base, as is done for goat evaluations and for current dairy cattle evaluations. Evaluations for crossbred animals with breed code XX will be reported on the base of the sire breed, which may cause some confusion because evaluations of animals from reciprocal crosses will be on different bases.
Conversions between all-breed and within-breed bases involved both a mean and the standard deviation ratio from Table 1
for traits with variance adjustment that differed by breed:
![]() |
![]() |
The within-breed, sire-breed, and all-breed models compared can be described with notation similar to Wiggans (1989):
![]() |
where y is a record preadjusted for multiplicative factors to account for heterogeneous variance and for age, parity, season, lactation length, and number of times milked daily that can vary by trait, time period, region, and breed; m is the management group mean; c is the interaction of herd with sire; p is the permanent environmental effect; a is the additive genetic effect; and e is the random residual. The animal model for bovines includes 3 new terms added since 1989: v is an age-parity-region-time effect (implemented in 1995); bFF is a regression of y on the cows inbreeding coefficient (implemented in 2005); and bHH is a regression of y on the cows coefficient of heterosis (used only in the all-breed model). The 3 models do not differ much in the terms they include, but rather in the data in y, the cows grouped in m, and the definitions of unknown-parent groups for missing ancestors in the relationship matrix.
| RESULTS |
|---|
|
|
|---|
|
99%), correlations exceeded 0.999 for Holsteins and 0.97 within other breeds. The PTA changed more for breeds with smaller populations, for bulls with fewer daughters, and for cows. Among other breeds, correlations tended to be lowest for Milking Shorthorns and highest for Jerseys. Correlations for PTA fat and PTA protein are not shown but were similar to those for PTA milk. Recent bulls were defined as those born since 1995 with daughters in
10 herds and within-breed reliability of
70% for yield or
40% for SCS, PL, or DPR. Recent cows were defined as those born since 1998 and reliability of
40% for yield or
30% for SCS, PL, or DPR.
|
10 crossbred daughters. Because many pure-bred animals have no crossbred progeny, changes in their PTA might also be explained by changes in the grouping of unknown dams and the addition of other breeds and crossbred cows to the management groups in mixed-breed herds. Those additional herdmates should increase accuracy but might also cause some bias if management of different breeds is not the same within herd. Crossbred animal PTA should be much more accurate with an all-breed than with a within-breed analysis because the relationship matrix can link to reliable sire PTA for breeds in both the maternal and paternal ancestry. For example, crossbred cows with the highest PTA PL each lived for several lactations and had a sire with high PTA PL and a grandsire with high PTA PL of another breed. If crossbred data are included only within sire breed for the data set, management group mates and genetic evaluations for maternal ancestors of other breeds are ignored.
The sire-breed model produced PTA very similar to official PTA from the within-breed model. Correlations were generally >0.999 for bulls with high reliability and also for recent bulls. Correlations for recent cows were about 0.995, and all of the largest changes were for F1 crossbred cows. A few crossbred cows were officially evaluated if enrolled in breed association grade-up programs, but most (93%) were not. The sire-breed model is an intermediate step between the current model and the all-breed model because records of crossbred cows are used, but many of their known ancestors and relatives of other breeds are treated as missing. Both the sire-breed and all-breed models add information from crossbred relatives but introduce some possibility of bias as compared with strictly purebred models.
Genetic Trends
Genetic trends for each breed and trait in the all-breed system are presented in Figures 1
to 6![]()
![]()
![]()
![]()
. Three trend validation tests (Interbull, 2004) were performed for each of 5 breeds (excluding Milking Shorthorn) and 5 traits (except that test 1 does not apply to PL). Interbull requires trend tests to be within 2 standard errors of 0.01 genetic standard deviations per year. Few biases were detected. For 64 of the 70 tests, 95% confidence intervals included the range of 0.01 to +0.01 genetic standard deviations per year.
|
|
|
|
|
|
|
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication October 25, 2006. Accepted for publication December 18, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. D. Norman, J. R. Wright, S. M. Hubbard, R. H. Miller, and J. L. Hutchison Reproductive status of Holstein and Jersey cows in the United States J Dairy Sci, July 1, 2009; 92(7): 3517 - 3528. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. B. Cole and D. J. Null Genetic evaluation of lactation persistency for five breeds of dairy cattle J Dairy Sci, May 1, 2009; 92(5): 2248 - 2258. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. D. Dechow, H. D. Norman, N. R. Zwald, C. M. Cowan, and O. M. Meland Relationship Between Individual Herd-Heritability Estimates and Sire Misidentification Rate J Dairy Sci, April 1, 2008; 91(4): 1640 - 1647. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. J. Heins, L. B. Hansen, A. J. Seykora, D. G. Johnson, J. G. Linn, J. E. Romano, and A. R. Hazel Crossbreds of Jersey x Holstein Compared with Pure Holsteins for Production, Fertility, and Body and Udder Measurements During First Lactation J Dairy Sci, March 1, 2008; 91(3): 1270 - 1278. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |