|
|
||||||||

* Centre for Genetic Improvement of Livestock Department of Animal and Poultry Science, University of Guelph, Ontario, Canada N1G 2W1,
Department of Animal Science Purdue University, West Lafayette, IN 47901
Corresponding author:
P. J. Boettcher; e-mail:
paulboettcher{at}anafi.it.
| ABSTRACT |
|---|
|
|
|---|
Key Words: genetics grazing sire evaluation
Abbreviation key: 305ME = 305-d mature equivalent, CI = calving interval, GxE = genotype x environmental interaction
| INTRODUCTION |
|---|
|
|
|---|
One factor that could cause GxE, however, is differences in feeding and management systems, in particular between herds that practice rotational grazing in the warm months versus herds that confine their cattle for the entire year. Despite the cold winters and short growing season, a small number of Canadian producers practice rotational or managed intensive grazing. Intensive rotational grazing and confinement feeding and management are different environments, and genetic evaluations of sires with daughters primarily in one environment may not be adequate predictors of progeny performance in another environment. In confinement systems, the ideal cows are those that produce the most milk with maximal energy inputs, while under pasture based systems the optimal genotypes are those that produce moderately high levels of milk given somewhat limited energy inputs. Records of cows from confinement systems dominate current genetic evaluations and selection decisions within the Canadian Holstein population.
Recent studies using data from different countries have shown evidence for GxE between grazing and conventional feeding systems. For example, Weigel et al. (2001) used lactation records and a sire model to estimate the genetic correlations of milk production across 17 countries. The correlation of milk yield in New Zealand (where grazing is nearly exclusively practiced) was only 0.89 and 0.88 with Canada and the US, respectively, versus 0.95 between the US and Canada. For comparison, the correlation was 0.96 between New Zealand and Australia (where grazing is practiced more commonly than in North America).
Some producers who use intensive grazing have become concerned that the current system does not produce and rank bulls adequately for their needs and have begun to use bulls proven in other countries (such as New Zealand), herd bulls, or crossbreeding. Recent studies (Weigel et al., 1999; Kearney et al., 2001a, 2001b) have examined records from grazing and conventionally managed herds in the United States and have concluded that effects of GxE on accuracy of sire rankings were not great enough to warrant a separate selection program directed toward graziers. The objective of this study was to evaluate effects of GxE on dairy traits among herds in Canada that practice intensive rotational grazing versus conventional methods that rely on stored feeds. Because of the importance of type traits in Canada (Canadian Dairy Network, 2001) in addition to production traits, effects of GxE for conformation were also examined.
| MATERIALS AND METHODS |
|---|
|
|
|---|
A control (conventional) group consisted of 34 large (>100 cows) free-stall operations in the same two provinces. The decision to use large free-stall herds for the comparison group was based on a number of factors. First, although these herds may not represent the present typical Canadian herd, the number of such herds in Canada (and particularly the number of cows in such herds) is continually increasing. Thus, these herds were expected to be more representative of the future typical herd in Canada. Second, these types of herds are now being targeted by AI companies for their progeny test programs (Welper, 2000), meaning that many of the future breeding decisions for the Canadian dairy population will be made based on phenotypes expressed in this type of environment. Finally, we expected that this type of comparison would have the greatest potential to reveal effects of GxE. The herds from Ontario in the conventional group had previously participated in a study on clinical lameness (Boettcher et al., 1998) and were chosen for this study because the authors were familiar with their management practices. In addition, the 10 largest free-stall herds in Nova Scotia were chosen to complete the conventional group. As a result of the requirement for pedigree data, all herds in both groups had registered Holsteins. Effects of GxE were studied for three types of traits: 1) production, 2) reproduction, and 3) conformation.
Production
For production, traits of primary interest were yields of milk, fat, and protein, but percentages of fat and protein were also considered in some of the analyses. Data were 305-d mature equivalent (305ME) DHI records from 1990 and later and included 6749 lactations from 2817 cows in the grazing herds and 29,371 lactations from 12,774 cows in the control herds. The number of records per cow was approximately the same in both groups, suggesting that longevity and culling rates were similar across environments.
To investigate GxE, two approaches were taken for production traits. The first approach tested for a scaling effect on yield traits, such as that observed by Veerkamp et al. (1995), in which the phenotypic differences in production of cows of differing genetic values varied according to the environment. To test for such a scaling effect, herds were separated into grazing and conventional groups and first lactation records of cows in each group (n = 2177) for grazing herds and 9502 for conventional) were regressed on the PTA of their respective sires. Significant differences between environments in the regression coefficients for PTA were assumed to indicate a scaling effect associated with the different environments. Weigel et al. (1999) applied a similar approach to production records from grazing and conventional herds in Wisconsin.
The model equation for this analysis was:
![]() | ([1]) |
where, yijk = 305ME record for milk, fat, or protein yields of cow k, calving at age j, in herd-year season i; HYSi = fixed effect of ith herd-year-season (4 seasons per year); Aj fixed effect of jth age group (j = 1,4); ß coefficient of regression of the production record on the PTA of the sire of cow k; and eijk = random residual. In theory, ß is expected to be approximately 1, under the assumption that the PTA represents the true transmitting ability of each bull. The PTA were one half times EBV from the May 2001 Canadian genetic evaluation.
A second analysis of the data was undertaken to examine whether yields of milk, fat, and protein, and percentages of fat and protein were actually different genetic traits in the grazing and conventional environments. This hypothesis was analyzed by estimating the genetic correlation between records for each trait expressed under the two different environmental conditions. Correlations significantly less than 1.0 were assumed to indicate GxE and statistically significant correlations <0.80 were assumed to indicate biologically important GxE (Robertson, 1959). Therefore, each cow had the information reported for a given trait in one environment and missing for the other. The genetic relationships from common sires was used to obtain the correlations. Among cows in the grazing herds, 83% had at least one paternal-half sib among the cows in the conventional group and 68% of the cows in the conventional herds had half-sibs in the grazing herds. The relationship matrix was constructed by tracing pedigrees for three generations.
A bivariate animal model with repeated records was assumed to describe the data. The model equation for this analysis was,
![]() | ([2]) |
where yijkt = 305ME milk, fat, or protein, or fat or protein percentage for cow k in environment t, herd-year-season-parity group i and parity-age group j; HYSPi = fixed effect of herd-year-season-parity (first or later) i; PAj = the fixed effect of the parity-age of calving group j (4 age groups in 3 parity groups, first, second, and later); Akt = the random additive genetic effect of cow k in environment t; PEk = the random permanent environmental effect of cow k; and eijk = random residual. Variance components were estimated by REML using the VCE 4.0 software of Groeneveld and Garcia-Cortes (1998).
Reproduction
For reproduction, effects of GxE were examined for calving interval (CI). The CI for a given lactation was simply recorded as the difference in days from the date of calving for that lactation and the following date of parturition. Because some cows did not have a subsequent calving (due to culling or censoring) the dataset included fewer records than had been available for production. In addition, CI was required to be >310 d, to avoid the inclusion of data from abortions or otherwise abnormally short pregnancies. This analysis included 3887 records (1784 cows) from the grazing herds and 16,461 records (8091 cows) from the conventional herds. the percentage of nonculled or noncensored records (obtained by dividing the numbers of usable CI records by the numbers of production records) was the same (63%) in both environments.
For CI, only the genetic correlation between environments was of interest and regression on sire PTA was not possible because no genetic evaluation for CI is performed in Canada. The model applied was the same as that used for production traits (Equation [2]
). In addition, the genetic parameters were reestimated using a model that considered phenotypic milk yield (305ME) as a covariate.
Conformation
Approximately 30 different conformation traits are recorded and evaluated in Canada, but many of these traits are correlated with each other. Therefore, we chose to concentrate our evaluation to the three traits that contribute directly to the Lifetime Profit Index, the official national selection index in Canada (Canadian Dairy Network, 2001). These traits are mammary system, feet and legs, and frame and capacity.
Although all of the herds participated in type classification, not all of the cows with production information had been classified. Only first-lactation records were used because classification in later lactations is not obligatory (i.e., data from later lactations are likely to be a biased sample) and only first-parity records contribute to the national genetic evaluation in Canada. For these reasons, the number of records for conformation were less than were available for production or reproduction. A total of 11,963 records were used, 2191 from the grazing herds and 9772 from the conventional herds.
Following the approach used for yield traits, regression was first used to test for scaling effects. Genetic correlations then were calculated across the different environments to test for other forms of GxE that could cause reranking of sires. The effects included in the models were similar to those included in the national evaluation for type.
For the regression of cow phenotype on sire PTA for type, the model equation was,
![]() | ([3]) |
where yijkl = conformation score of cow l in age at classification and stage of lactation group i, round-classifier group j (a different classifier visits the herd on each date of classification), and herd-round classifier group k; SAi = fixed effect of age and stage of lactation group i (4 age and 11 stage groups); RCj = fixed effect of round-classifier j; HRCk = the fixed effect of herd-round-classifier k; ß coefficient of regression of the type scores on the PTA of the sire of cow l; and eijkl = random residual. The PTA for type in Canada are expressed on a standardized scale (µ = 0 and
= 5), rather than on the phenotypic scale and, therefore, the regression coefficient was not expected to be 1.0. Students t-test was used to determine whether the regression coefficients differed in the two environments.
The model equation for the estimation of heritabilities and genetic correlations across the two groups was,
![]() | ([4]) |
where yijklt = conformation score of cow l in environment t, in age and stage of lactation group i, round-classifier group j, and herd-round-classifier group k; SAl, RCj, and HRCk are the same as in equation [3]
; Alt = the random additive genetic effect of cow l in environment t; and eijkl = random residual.
| RESULTS |
|---|
|
|
|---|
|
Means of EBV are in Table 2
for the traits for which EBV are calculated for cows in Canada. Means are significantly greater (P < 0.01) in the conventional herds for all traits except for frame and capacity. The differences in EBV explain approximately 15 to 20% of the differences in phenotypic means (Table 1
) for yield traits and about 25% for type traits.
|
|
Veerkamp et al. (1995) reported a similar scaling effect when comparing the performance of selected and control lines of dairy cattle when fed two different TMR (neither including pasture) with either high or low levels of concentrate. The regression of actual performance on pedigree index was slightly lower for the cows on the low concentrate diet and on the high concentrate diet.
No scaling effects were observed for the type traits (Table 3
). Regression coefficients for the type traits were not expected to be 1 because EBV for type in Canada are transformed to a standard scale.
Estimates of heritabilities within each environment and genetic correlations across environments are in Table 4
for all nine traits evaluated. Heritabilities were generally lower than the grazing herds, but these differences were nonsignificant for most of the traits. For the yield traits, heritabilities fell in the general range of 0.35 to 0.40 for the conventional herds versus only 0.30 to 0.35 for the grazing herds. The increased estimates of heritability in the conventional herds coincide with the trends observed in regression coefficients (Table 3
), which suggested that genetic differences among animals were expressed to a greater degree in the conventionally managed herds. Heritabilities of percentage traits were quite high but consistent with other estimates previously made on the Canadian Holstein population (Boettcher and Gibson, 1997; Z. Lui, 1996, unpublished technical report).
|
0.92. Even less evidence of GxE was observed for reproduction and conformation traits. Some indication for GxE was observed for CI, as estimated heritability in grazing herds (0.052) was nearly double that in conventional herds (0.027) and estimated genetic correlation across groups was only 0.64. However, these differences were not significant (P > 0.05). When milk yield was added to the model as a covariate, results changed very little except that the estimate of heritability in the grazing herds decreased to 0.036.
Among the type traits, heritability for mammary system was greater in the conventional herds (0.20 vs 0.13) but genetic correlations between groups was 1.00 (SE of zero was an artifact of the estimation procedure). Perhaps the greater heritability in conventional herds was associated with higher mean production, and thus more stress on udder attachments. For feet and legs and frame and capacity, heritabilities were not different across environments and genetic correlations were not different from unity. These two results were interesting considering that 1) with respect to feet and legs, cows in grazing herds would generally have been exposed to much different environments (less concrete) from cows in the free-stall herds, and 2) for frame and capacity, the diets of cows in the grazing herds were expected to be quite different (more forage) from those in the free-stall herds.
| DISCUSSION AND CONCLUSIONS |
|---|
|
|
|---|
A number of factors could explain this lack of observed GxE. First, due to the cool climate and short growing seasons in Canada, the herds in this study practiced grazing for only about half of the year. The rest of the time feeding and management would be expected to be similar to the rest of Canada. Therefore, the two systems of management were not as extremely different as could be observed in another country such as the United States or between countries such as New Zealand and Canada. Secondly, all of the herds in both groups were enrolled in DHI and had primarily registered animals. These factors may suggest that other management practices were in common across groups. Another possibility is that expression of genes is not particularly different in quite diverse environments.
The level of GxE is too low to justify a separate national sire proving program for graziers, particularly in Canada, where, due to the climate, the number of producers practicing is relatively small. Weigel et al. (1999) concluded that a separate sire proving program for graziers would require sampling of 600 to 700 bulls per year to achieve the same genetic progress as using sires based on a national program for which results were genetically correlated by 0.90. Weigel et al. (1999) also questioned the practice of using sires proven in typically grazing countries such as New Zealand or Ireland, because progeny testing programs in those countries are small and selection intensity would be lost.
One approach that could help resolve this problem is borderless international sire evaluations (Lohuis and Dekkers, 1998) in which herds are clustered (Weigel and Rekaya, 2000) by management practices, climate, and other factors instead of by country. Assuming that one of the clusters included primarily grazing herds and sires were progeny tested in multiple clusters, then perhaps graziers could select based on rankings for daughter performance in that cluster. However, adoption of such a system would require the navigation around various data recording, computational, and political obstacles.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Received for publication October 22, 2001. Accepted for publication December 11, 2001.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
W. J. Fulkerson, T. M. Davison, S. C. Garcia, G. Hough, M. E. Goddard, R. Dobos, and M. Blockey Holstein-Friesian Dairy Cows Under a Predominantly Grazing System: Interaction Between Genotype and Environment J Dairy Sci, February 1, 2008; 91(2): 826 - 839. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Haskell, S. Brotherstone, A. B. Lawrence, and I. M. S. White Characterization of the Dairy Farm Environment in Great Britain and the Effect of the Farm Environment on Cow Life Span J Dairy Sci, November 1, 2007; 90(11): 5316 - 5323. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. G. Fahey, M. M. Schutz, D. L. Lofgren, A. P. Schinckel, and T. S. Stewart Genotype by Environment Interaction for Production Traits While Accounting for Heteroscedasticity J Dairy Sci, August 1, 2007; 90(8): 3889 - 3899. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Bryant, N. Lopez-Villalobos, J. E. Pryce, C. W. Holmes, D. L. Johnson, and D. J. Garrick Environmental Sensitivity in New Zealand Dairy Cattle J Dairy Sci, March 1, 2007; 90(3): 1538 - 1547. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Rozzi, F. Miglior, and K. J. Hand A Total Merit Selection Index for Ontario Organic Dairy Farmers J Dairy Sci, March 1, 2007; 90(3): 1584 - 1593. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Konig, G. Dietl, I. Raeder, and H. H. Swalve Genetic Relationships for Dairy Performance Between Large-Scale and Small-Scale Farm Conditions J Dairy Sci, November 1, 2005; 88(11): 4087 - 4096. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Horan, P. Dillon, P. Faverdin, L. Delaby, F. Buckley, and M. Rath The Interaction of Strain of Holstein-Friesian Cows and Pasture-Based Feed Systems on Milk Yield, Body Weight, and Body Condition Score J Dairy Sci, March 1, 2005; 88(3): 1231 - 1243. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. R. Zwald, K. A. Weigel, and T. J. Lawlor Jr. Genetic Parameters for Conformation Traits in Herds that Differ in Mean Final Score and Completeness of Pedigree and Performance Data J Dairy Sci, January 1, 2005; 88(1): 386 - 391. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. F. Kearney, M. M. Schutz, P. J. Boettcher, and K. A. Weigel Genotype x Environment Interaction for Grazing Versus Confinement. I. Production Traits J Dairy Sci, February 1, 2004; 87(2): 501 - 509. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |