Short communication: Genetic and environmental relationships between milk yield and kidding interval in dairy goats
Article Outline
Abstract
First-parity 305-d milk yield and intervals between first and second kiddings from 1975 through 2005 were analyzed to estimate genetic and environmental parameters for United States Alpine, LaMancha, Nubian, Saanen, and Toggenburg dairy goats. The data set included 43,612 does with first-parity milk yield and 25,863 does with first kidding interval. Parameters were estimated by REML using bivariate mixed models within and across breeds. Fixed effects were month and herd-year of kidding within breed and month and breed-herd-year of kidding across breeds. Random effects were animal and residual. The heritability estimate for first-parity milk yield was 0.36
±
0.01 across breeds and ranged from 0.35 to 0.38 within breed; heritability for first kidding interval was 0.05
±
0.01 across breeds and ranged from 0.00 to 0.15 within breed. The estimate of the genetic correlation between first-parity milk yield and first kidding interval was positive (unfavorable) across breeds (0.35
±
0.09). Estimates of environmental correlations within and across breeds were positive (0.16–0.25). The presence of unfavorable genetic relationships between milk yield and kidding interval indicates a need to include reproductive performance as a selection criterion.
Key words: genetic parameter, kidding interval, milk yield, dairy goat
Several studies (e.g., Pryce et al., 2004, Andersen-Ranberg et al., 2005) have indicated unfavorable genetic and environmental correlations between milk yield and fertility in dairy cattle. Because of the economic importance of fertility, reproduction traits have been incorporated as a selection criterion for dairy cattle in many countries (e.g., VanRaden, 2004). Knowledge of those relationships in dairy goats is limited. Estimates of heritabilities and correlations are needed (Montaldo and Manfredi, 2002). The objective of this study was to estimate heritabilities of milk yield and kidding interval for 5 goat breeds based on first kiddings as well as genetic, environmental, and phenotypic correlations between those traits.
Milk yield and kidding interval data from DHI lactation records in the USDA national dairy database (Beltsville, MD) and provided to the American Dairy Goat Association (Spindale, NC) were analyzed for kiddings between 1975 and 2005. Pedigree information was available for all does with records and their ancestors for animals born in 1971 or later. First-parity, mature-equivalent, 305-d milk yield (Wiggans and Hubbard, 2001), and interval between first and second kiddings were analyzed. Only data from lactations with ≥100 DIM and from herd-year-breed combinations and sires with ≥2 daughter records were retained for final analysis. Records with first kidding intervals of <170 or >730 d were considered to be abnormal and were excluded. Numbers of sires and does for first-parity milk yield and first kidding interval are shown in Table 1 along with means (unadjusted) and phenotypic SD (unadjusted and from the analysis).
Table 1. Numbers of sires and does, means (unadjusted), and phenotypic SD (unadjusted and from model analysis) for first-parity milk yield and interval between first and second kiddings of United States dairy goats by breed
| Item | Alpine | LaMancha | Nubian | Saanen | Toggenburg | All breeds |
|---|---|---|---|---|---|---|
| Sires, n | 2,346 | 994 | 2,544 | 1,123 | 1,036 | 8,043 |
| First-parity milk yield | ||||||
| 15,147 | 4,633 | 12,227 | 5,828 | 5,777 | 43,612 | |
| 1,090 | 978 | 928 | 1,187 | 1,050 | 1,040 | |
| 322 | 290 | 303 | 392 | 317 | 379 | |
| 245 | 216 | 210 | 278 | 235 | 237 | |
| First kidding interval | ||||||
| 9,222 | 2,623 | 7,075 | 3,451 | 3,492 | 25,863 | |
| 381 | 374 | 373 | 384 | 382 | 379 | |
| 89 | 90 | 87 | 93 | 89 | 89 | |
| 82 | 82 | 75 | 77 | 78 | 79 |
Data were analyzed within breed (Alpine, LaMancha, Nubian, Saanen, and Toggenburg) with bivariate mixed models that included fixed effects for month and herd-year of kidding. Random effects were animal (with all relationships considered through the inverse of the numerator relationship matrix) and residual. Covariance components were estimated with ASReml (Gilmour et al., 2006). All breeds also were analyzed together by including a breed-herd-year instead of herd-year effect in the model.
Heritability estimates for first-parity milk yield (Table 2) were large (0.35– 0.38) compared with previous estimates (Montaldo and Manfredi, 2002; Valencia et al., 2007; Torres-Vázquez et al., 2009), which shows the possibility for increased accuracy of selection and more rapid genetic improvement for productivity (Montaldo et al., 1994). The heritability estimate for first-parity milk yield from the analysis across breeds was 0.36
±
0.01.
Table 2. Estimates of heritability and genetic, environmental, and phenotypic correlations for first-parity milk yield and interval between first and second kiddings of United States dairy goats by breed
| Parameter | Alpine | LaMancha | Nubian | Saanen | Toggenburg | All breeds |
|---|---|---|---|---|---|---|
| Heritability | ||||||
| 0.36 | 0.38 | 0.38 | 0.35 | 0.36 | 0.36 | |
| 0.06 | 0.15 | 0.00 | 0.07 | 0.02 | 0.05 | |
| Genetic correlation | 0.31 | 0.36 | — | 0.26 | 0.38 | 0.35 |
| Environmental correlation | 0.22 | 0.21 | 0.21 | 0.16 | 0.25 | 0.20 |
| Phenotypic correlation | ||||||
| 0.21 | 0.22 | 0.21 | 0.17 | 0.23 | 0.20 | |
| 0.17 | 0.14 | 0.15 | 0.10 | 0.15 | 0.151 | |
1Weighted mean across breeds. |
Heritability estimates for first kidding interval (Table 2) ranged from 0.00 to 0.15. The heritability estimate of 0.00 for Nubians may be the result of sampling error; larger data sets may be needed to estimate genetic parameters accurately for first kidding interval. The across-breed heritability estimate for first kidding interval was 0.05
±
0.01, which was slightly larger than the estimates of 0.02 for Polish goats and 0.03 for Norwegian goats reported by Bagnicka et al. (2007). Breed variation in heritability estimates for first kidding interval may be caused by low heritability for all breeds and large sampling variances, as shown by the relatively large standard errors (Table 2). Previous heritability estimates for kidding interval from small samples (Odubote, 1996; Ribeiro et al., 2000; Lôbo and da Silva, 2005) probably had large standard errors; however, those estimates ranged from 0.01 to 0.09, which indicated that heritability is less than 0.10.
Estimates of genetic correlation between first-parity milk yield and first kidding interval (Table 2) were positive (unfavorable) for Alpines (0.31), LaManchas (0.36), Saanens (0.26), and Toggenburgs (0.38). A likelihood ratio test (Gilmour et al., 2006) indicated that the genetic correlation was different from zero for Alpines (P
<
0.01) and LaManchas (P
<
0.05) but not for Saanens or Toggenburgs (P
>
0.05). The genetic correlation was not estimable for Nubians because the heritability estimate for first kidding interval was near zero. The across-breed genetic correlation (0.35) based on data from all breeds also was different from zero (P
<
0.01). The low heritability estimates for first kidding interval compared with first-parity milk yield and positive genetic correlations between first kidding interval and first-parity milk yield agree with reports for reproductive traits of dairy cattle (Castillo-Juarez et al., 2000, Pryce et al., 2004, Andersen-Ranberg et al., 2005). VanRaden et al. (2004) estimated the genetic correlation between days open and milk yield to be 0.38 for United States Holsteins. Because of the positive genetic correlation between first-parity milk yield and first kidding interval, selection indexes should include first kidding interval to prevent a decline in fertility.
Estimates of the environmental correlations between first-parity milk yield and first kidding interval (Table 2) were positive (unfavorable) for all breeds (0.16–0.25), with small standard errors. The across-breed environmental correlation was 0.20. Because of the small heritability estimates for first kidding interval, phenotypic correlations (Table 2) were determined almost completely by environmental correlations. Similar estimates have been reported for dairy cattle (Pryce et al., 2004, Andersen-Ranberg et al., 2005, Mackey et al., 2007). Unadjusted phenotypic correlations between milk yield and kidding interval, which include the covariances across herd-years (Table 2), were also positive, but they were smaller than phenotypic correlations from model residuals. The environmental correlations indicate some environmental antagonism between milk yield and fertility within herd-year. However, the unadjusted phenotypic correlations, which include herd-year effects, were smaller than environmental correlations, indicating favorable management effects across herd-years for both traits.
Acknowledgments
The authors thank S. M. Hubbard (Animal Improvement Programs Laboratory, Agricultural Research Service, USDA, Beltsville, MD) for suggestions on manuscript improvement.
References
- . Heritabilities, genetic correlations, and genetic change for female fertility and protein yield in Norwegian dairy cattle. J. Dairy Sci. 2005;88:348–355
- . Heritability for reproduction traits in Polish and Norwegian populations of dairy goat. Small Rumin. Res. 2007;68:256–262
- . Effect of herd environment on the genetic and phenotypic relationships among milk yield, conception rate, and somatic cell score in Holstein cattle. J. Dairy Sci. 2000;83:807–814
- . ASReml User Guide. Hemel Hempstead, UK: Release 2.0. VSN International Ltd.; 2006;
- . Parâmetros genéticos para características de interesse econômico em cabras das raças Saanen e Anglo-nubiana. [Genetic parameters for economics traits in Saanen and Anglo-Nubian goats.]. Revista Ciência Agronômica. 2005;36:104–110
- . Associations between genetic merit for milk production and animal parameters and the fertility performance of dairy cows. Animal. 2007;1:29–43
- Montaldo, H. H., and E. Manfredi. 2002. Organisation of selection programmes for dairy goats. Commun. no. 01–35 in Proc. 7th World Congr. Genet. Appl. Livest. Prod., Montpellier, France. INRA, Paris, France.
- . Respuesta esperada a la selección para producción de leche en caprinos con pruebas de progenie, MOET y selección masal. [Expected genetic response for milk production in goats with progeny testing, MOET and mass selection.]. Arch. Zootec. 1994;43:335–343
- . Genetic parameters for litter size at birth and kidding interval in the West African Dwarf goats. Small Rumin. Res. 1996;20:261–265
- . Fertility in the high-producing dairy cow. Livest. Prod. Sci. 2004;86:125–135
- . Genetic and environmental effects on the age at first kidding and kidding interval in a Saanen goat herd. Ars Vet. 2000;16:192–197
- . Genetic and phenotypic parameters of milk yield, milk composition and age at first kidding in Saanen goats from Mexico. Livest. Sci. 2009;126:147–153
- . Genetic and phenotypic parameters for lactation traits in a flock of Saanen goats in Mexico. Small Rumin. Res. 2007;68:318–322
- . Invited review: Selection on net merit to improve lifetime profit. J. Dairy Sci. 2004;87:3125–3131
- . Development of a national genetic evaluation for cow fertility. J. Dairy Sci. 2004;87:2285–2292
- . Genetic evaluation of yield and type traits of dairy goats in the United States. J. Dairy Sci. 2001;84(E-Suppl):E69–E73
PII: S0022-0302(10)70296-4
doi:10.3168/jds.2009-2593
© 2010 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.
