|
|
||||||||

* Department of Dairy and Animal Science, The Pennsylvania State University, University Park 16802
Animal Improvement Programs Laboratory, Agricultural Research Service, USDA, Beltsville, MD 20705-2350
1 Corresponding author: cdechow{at}psu.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: herd heritability daughter-dam regression daughter-sire regression
Dechow and Norman (2007) estimated individual herd heritabilities for milk yield, fat yield, protein yield, and SCS for 1,939 herds using regression techniques. Correlations of regression heritabilities with individual herd heritability estimates generated with REML ranged from 0.45 to 0.68 for a subset of the 45 largest herds. However, attempts to adjust records to a constant genetic variance using the individual herd heritability estimates did not appreciably improve the accuracy of genetic evaluations. Lack of improvements when data were adjusted for individual herd heritabilities could indicate that the true genetic variance was not different between herds with high and low heritability.
Data quality issues, such as sire misidentification rate, would also be expected to impact herd heritability estimates regardless of the herds genetic variance. Individual herd heritabilities were generated for 20,902 herds and were merged with sire-misidentification rate as determined by DNA maker analysis for a subset of 230 herds (Dechow et al., 2008). The correlation between sire misidentification rate and a principal component for all measures of herd heritability was –0.50.
The first objective of this study was to compare (co)-variance parameter estimates for milk yield, fat yield, protein yield, and SCS among subsets of data that were pooled from herds with high, medium, or low individual herd heritability estimates. The second objective was to assess how individual herd heritability estimates compare with REML heritability estimates when data were pooled for herds with similar individual heritability estimates.
Individual herd heritabilities for 20,902 herds were estimated (Dechow et al., 2008). Briefly, the model used to generate individual herd heritabilities in ASREML (Gilmour et al., 2006) was
![]() | [1] |
where yijklmno = mature-equivalent milk yield, mature-equivalent fat yield, mature-equivalent protein yield, or SCS for the ith record of cow j of breed k, parity l, calving in herd-year-season m, in state n and herd o; BD = the fixed effect of breed, bl = coefficient for fixed regression on age (A) nested within parity, HYS = fixed effect of herd-year-season, bd = coefficient for fixed regression on dam record nested within state (D), bs = coefficient for fixed regression on sire PTA (S) nested within state, bsd = coefficient for fixed regression on the interaction between sire PTA and herd standard deviation (SD), bdo = coefficient for random regression on dam record nested within herd (F), bso = coefficient for random regression on sire PTA nested within herd (G), and eijklmno = effect of random residual. The random regression coefficients (bdo and bso) were assumed to be correlated.
Individual herd heritability from daughter-dam regression was 2(bd + bdo). Individual herd heritability for daughter-sire regression was:
![]() |
where bs, bsd, SDo, and bso are as defined in model 1; SDUS = genetic standard deviation assumed for USDA-DHIA Holstein genetic evaluations (655 kg for mature-equivalent milk); and R = average sire PTA reliability for all cows in a herd.
A principal component for heritability was generated for each herd using 8 herd heritability measures (daughter-dam and daughter-sire estimates for each of the 4 traits), and herds were stratified according to the principal component for herd heritability. Mature-equivalent milk, fat, and protein yields, and SCS from Holstein cows calving between August 2000 and August 2005 were extracted from the national dairy database if the cow calved in a randomly chosen herd with a principal component of heritability that fell into one of the following percentiles: lowest 5% (P5), 47 to 53% (P50), or above 95% (P95). The target number of cows in each subset was 20,000, and only cows with a first-lactation record available, an age at calving ranging from 20 to 120 mo, and that were from a sire with at least 10 daughters in the subset were retained. The total number of herds, cows, records, and sires represented in each subset in addition to average milk yield, fat yield, protein yield, and SCS are reported in Table 1
.
|
![]() |
where Y = a vector of records from P5, P50, and P95 for either milk yield, fat yield, protein yield, or SCS; β = a vector of fixed herd-year-season of calving effects with 6 bimonthly calving seasons, and linear and quadratic effects for age at calving nested with parity groups 1, 2, and
3; a = a vector of random animal effects; p = a vector of random permanent environmental effects; X, Z, and W are the corresponding incidence matrices; and
= a vector of random errors. All traits were also analyzed with single-trait models that included an additional random effect for sire by herd interaction.
The average daughter-dam and daughter-sire heritability estimates for P5, P50, and P95 are reported in Table 2
. The principal component for herd heritability was successful in stratifying herds according to herd heritability for all traits. Daughter-dam heritability estimates were higher than daughter-sire heritability estimates for all traits and in all subsets. Average daughter-dam heritability estimates ranged from 0.21 for protein yield in P5 to 0.53 for fat yield in P95. Daughter-sire herd heritability estimates ranged from an average of 0.09 for fat yield, protein yield, and SCS for P5 to 0.40 for protein yield in P95.
|
|
The highest phenotypic variance for all traits was observed for P95. Several studies have described a tendency toward higher heritability estimates as phenotypic variance increases (Lofgren et al., 1985; Vinson, 1987; Van Tassell et al., 1999). However, phenotypic variance was not generally higher for P50 than for P5.
Dechow and Norman (2007) speculated that daughter-sire PTA regression could be depressed in the presence of genotype by environment interaction. Sire by herd interactions might be expected to be higher in low daughter-sire heritability herds if genotype by environment interactions were responsible for lower individual herd heritability estimates. However, estimates of sire by herd interaction ranged from 0.84% (SCS) to 2.15% (milk yield) in P5, 1.16% (milk yield) to 1.98% (protein yield) in P50, and from 1.78% (milk yield) to 2.72% (fat yield) in P95. Genetic correlations among P5, P50, and P95 were also high, ranging from 0.86 to 1. There was little evidence that genotype by environment interaction depressed daughter-sire herd heritability estimates.
Genetic parameter estimates for pooled samples of records from herds that were stratified into high, medium, and low herd heritability groups indicated that the regression methods used to estimate individual herd heritabilities were effective. Daughter-dam heritability estimates were generally much higher than REML heritability estimates, whereas the magnitude of daughter-sire and REML heritability estimates were of a similar magnitude. Additive genetic variance increased progressively as herd heritability increased, whereas permanent environmental variance decreased as herd heritability increased.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication August 17, 2007. Accepted for publication January 1, 2008.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |