|
|
||||||||
1 Dairy Production Department, Teagasc, Dairy Production Research Centre Moorepark, Fermoy, Co. Cork, Ireland
2 Department of Animal Science, Faculty of Agriculture, University College Dublin, Belfield, Ireland
3 3INRA, UMR Production du Lait, 35590 St Gilles, France
Corresponding author: B. Horan; e-mail: bhoran{at}moorepark.teagasc.ie.
| ABSTRACT |
|---|
|
|
|---|
Key Words: Holstein-Friesian strain grass-based feed system
Abbreviation key: FS = feed system, HCFS = high concentrate feed system, HD = high durability, HF = Holstein-Friesian, HP = high production, HSFS = high stocking rate feed system, MPFS = high grass feed system, NZ = New Zealand, PGSSH = postgrazing sward surface height.
| INTRODUCTION |
|---|
|
|
|---|
Up to the mid1980s, the predominant breed of dairy cow in Ireland was the British Friesian. In Ireland over the last 20 yr, the use of North American Holstein-Friesian (HF) genetics has dominated, increasing from 9% in 1990 to 65% in 2001 (Evans et al., 2004). These animals have been aggressively selected for high milk production in a predominantly confinement environment (Rauw et al., 1998). The popularity of the North American HF was most likely due to its increased productivity over dairy breeds in a market that rewarded for milk volume with lesser incentive for milk solids content. However, aggressive genetic improvement within the North American HF for increased milk yield has resulted in a now well-documented negative effect on cow fertility (Pryce and Veerkamp, 2001). There is a known antagonistic genetic relationship between milk yield and reproductive traits (Evans et al., 2002; Berry et al., 2003). Pryce et al. (1998) estimated an increase of the calving interval by 5 to 10 d per 1000 kg of additional milk yield when selecting for milk yield only. These values are supported by the results from controlled experiments comparing genetic groups across feeding systems (Buckley et al., 2000; Kennedy et al., 2002).
In contrast, the New Zealand HF has been selected for milk production, feed efficiency, and survivability on a predominantly grass-based diet with very limited concentrate supplementation (Harris and Kolver, 2001). In New Zealand, the North American HF cows are of higher BW, produce more milk volume and protein yield, and have lower concentrations of milk fat and protein than the New Zealand HF (Harris and Kolver, 2001).
It has until recently been widely assumed that if a strain is the best for one environment (e.g., confinement), it will also be the best in all other environments. Although previous studies, across a relatively narrow range of genetic merit and/or feeding systems, have failed to show a genotype x feed system (FS) interaction, there is now evidence that high merit North American HF cows may be unable to express their full genetic potential when in a grass-based system of milk production (Kennedy et al., 2002). The importance of interactions between genotype and environment are substantial considering the diverse range of environments worldwide in which animals are now being managed.
The objective of the current study was therefore to investigate the effects of 3 strains of HF and 3 grass-based feed systems and their interactions on the milk production, BW, and BCS within Irish pasture-based systems of milk production.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Animals
Three strains of HF cow were compared: high production North American (HP), high durability North American (HD), and New Zealand (NZ). The mean pedigree index for each strain is displayed in Table 1
. The pedigree index for each cow was calculated as 0.50 x sire predicted difference + 0.25 x maternal grandsire predicted difference + 0.125 x maternal great-grandsire predicted difference. The predicted difference of the sires and maternal grandsires were from the February 2004 international evaluations of the INTERBULL Animal Center (Uppsala, Sweden) using the technique known as MACE (multiple-trait across-country evaluation).
|
A total of 99, 117, and 117 animals were used in yr 1, 2, and 3, respectively, divided equally between strains of HF and FS. The number of animals of each strain of HF and FS over the 3 yr of the study is displayed in Table 2
. In 2001, all 99 cows were 3 yr old; in 2002, 45 animals were 3 yr old and 72 were parity 2; and in 2003, 9 animals were 3 yr old, 45 were parity 2, and 63 were parity 3. Fifty-nine animals remained in the same FS for the duration of the study (3 yr). All primiparous cows were on a similar feeding regimen for the first 4 wk of lactation. Animals were selected within strain into groups of 3, on the basis of calving date, milk production in the first 4 wk of lactation, and BW, and then randomly assigned to 1 of 3 systems. Once allocated, animals were retained on the same FS in subsequent lactations.
|
|
The subpaddock residency time of the cows in both the HSFS and MPFS was the same, i.e., similar grazing rotation; therefore, the cows in the HSFS were forced to graze to a lower PGSSH throughout the grazing season. The grazing rotation in the HCFS moved independently to that in the MPFS and HSFS, however the HCFS group grazed to a similar PGSSH as the animals on MPFS. Pasture quality was maintained in this system by removing surplus grass throughout the experiment as silage (ODonovan, 2000).
Animal Measurements
During the 3 yr of the study, individual milk yields were recorded on 5 consecutive days/week. The milk fat, protein, and lactose concentrations were determined using a Milkoscan 203 (Foss Electric DK-3400, Hillerød, Denmark) from successive morning and evening samples collected once weekly. The BW of each animal was recorded weekly. Each BW was recorded electronically, using portable weighing scales and Winweigh software package (Tru-Test Ltd., Auckland, New Zealand). The scales were calibrated weekly against known weights. The BCS was recorded every 3 wk during the lactation on a 0-to-5 scale (0 = emaciated, 5 = extremely fat) with increments of 0.25 as outlined by Lowman et al. (1976). The first BW and BCS data were recorded at 24 h postpartum.
Feed Measurements
Pregrazing herbage yield (above 4 cm horizon) was determined on each grazing paddock based on 4 strips (0.80 m wide; 4.5 to 5.5 m long) of grass cut with an Agria mower (Agria-Werke GmbH, Mockmuhl/Wurtt, Germany). The grass from each strip was weighed and sampled and a subsample was dried overnight at 90°C for DM determination. The remaining herbage samples from each paddock were bulked and a further subsample taken (~100 g) and freeze-dried and used for chemical analysis.
During each rotation, 30 pregrazing sward surface heights were recorded for each strain within each paddock immediately before grazing with a further 30 post-grazing sward surface heights taken immediately after grazing (Hutchings, 1991).
Chemical analysis.
The composite herbage samples for each week were analyzed for modified acid detergent fiber (Clancy and Wilson, 1966), OM digestibility (Morgan et al., 1989), and Kjeldahl N. Similarly, in periods of silage supplementation, a composite grass silage sample for each week was analyzed for residual moisture, DM digestibility (Tilley and Terry, 1963), Kjeldahl N, modified ADF, and NDF (Van Soest et al., 1966). Concentrates were sampled weekly, bulked over each month, and analyzed for DM, total N, crude fiber, NDF, oil, and ash.
Data Handling
The data set was divided into 2 periods for the purpose of statistical analysis. Firstly, the production, BW, and BCS data were analyzed over the complete lactation (wk 3 of year to wk 53, inclusive) in all 3 yr. A second data set was created to examine the performance of each group exclusively during the period in lactation when the 3 feeding systems were in place (wk 10 of year to wk 45, inclusive) in all 3 yr. In early lactation (February and March) and late lactation (mid-November/December), when cows were housed by night or by night and day, the various feed systems were not in place.
Statistical Analysis
Data were analyzed using a repeated measures model (PROC MIXED) described below using the statistical procedures of SAS (SAS Institute, 2002). Cow was included as a random effect, and year, parity, strain of HF, and feed system were included as fixed effects. Each year, the newly introduced animals were selected within strain and parity on calving date and pre-experimental milk yield. Animals present from the previous year were maintained in the same feed system. To improve the accuracy of the model, pre-experimental milk yield, BW, and BCS were used as covariates specific to the traits being analyzed. The following model was used:
![]() |
where Yijkl = the performance of the animal in yr i of parity j and strain k on feed l, Ri = the effect of ith yr (i = 1, 2, 3), Lj = parity (j = 1, 2, 3), Sk = strain of HF (k = HP, HD, or NZ); Fl = the effect of the lth feed system (l = MPFS, HSFS, HCFS), and eijkl = the residual error term.
For the milk production variables, the covariates used were pre-experimental milk production, individual animal predicted difference for milk production, and calving date for each animal on the study. Models for BW and BCS included covariates for pre-experimental BW and condition score, gestation length, and calving date. Owing to the differences between strains in terms of the pre-experimental values, these covariates were centered within strain before inclusion. That is, the deviations from the strain mean were used as the covariates. The incorporation of individual animal covariates within the model reduced the residual error term, therefore explaining more of the variation within strains.
Orthogonal contrasts (Snedecor and Cochran, 1980) were used to make a number of comparisons. The PRODL contrast refers to the comparison of animals selected purely for high milk production potential (HP strain) with animals that are not selected aggressively for milk production (HD and NZ strains)[PRODL; HP vs. (HD +NZ)/2]. The LOWP contrast refers to the comparison of exclusively lower production groups, therefore comparing the lower production North American HF (HD strain) to the NZ strain (LOWP; HD vs. NZ). The effect of FS was determined by making the following comparisons with the MPFS (control): MPFS vs. HSFS and MPFS vs. HCFS. The effect of strain of HF on response to stocking rate and concentrate supplementation was also tested using the Students t-test to compare the differential in performance of each strain in the MPFS with the performance in the HSFS and HCFS during the feed system period.
| RESULTS |
|---|
|
|
|---|
|
|
|
Parity had a significant effect on total lactational milk production (significance levels are included in Table 6
). Third-parity cows had the highest (P < 0.001) yield of milk (7436 kg), SCM (7449 kg), peak daily milk (37.1 kg), fat (307 kg), protein (259 kg), and lactose (340 kg), the first-parity cows had the lowest (P < 0.001) (5533, 5322, 26.8, 233, 199, and 266 kg, respectively), and the second-parity cows were intermediate for all variables (6850, 6698, 33.7, 281, 242, and 317 kg, respectively). First-parity animals also had lower (P < 0.01) protein concentration (35.2 g/kg) compared with animals of second and third parity (35.9 and 36.0 g/kg, respectively).
Feed system period analysis.
During the feed system period, the interaction of strain of HF and FS as well as effects of strain of HF and FS were observed (Table 7
). Feed system had a significant effect on all variables measured. The cows in the HCFS produced the highest (P < 0.001) daily yield of milk (23.4 kg/d), SCM (22.6 kg/d), fat (933 g/d), protein (851 g/d), and lactose (1111 g/d), the cows in the HSFS the lowest (P < 0.05) (19.5 kg, 18.9 kg, 806 g, 694 g, and 908 g, per d, respectively), and that of the cows in MPFS was intermediate (20.2 kg, 19.7 kg, 831 g, 726 g, and 939g, per day, respectively). The HCFS group had lower (P < 0.01) fat composition (40.5 g/kg) compared with the HSFS group (41.6 g/kg), with MPFS intermediate (41.1 g/kg). The cows in the HCFS produced the highest protein and lactose content (36.8 and 47.5 g/kg, respectively), with no significant difference between the MPFS (36.1 and 46.4 g/kg, respectively) and HSFS (36.0 and 46.7 g/kg, respectively).
|
|
|
|
Animals of third parity were heaviest (P < 0.001) immediately postpartum, at nadir BW, and at drying-off (580, 522, and 605 kg, respectively), first-lactation animals the lightest (501, 438, and 531 kg, respectively), and second-lactation animals were intermediate (545, 495, and 573 kg, respectively). Primiparous animals had the highest BCS immediately postpartum (3.26), and at drying-off (2.93), whereas the second-parity animals had the lowest BCS immediately postpartum (3.18), and the third-parity animals had the lowest BCS at drying-off (2.87). Second-parity animals had lower (P < 0.05) BW and BCS loss between calving and nadir (50.6 kg and 0.56, respectively) compared with the primiparous (63.5 kg and 0.63, respectively) and third-parity animals (57.8 kg and 0.69, respectively).
Feed system period analysis.
Strain of HF, FS, and parity had significant effects on BW and BCS during the feed system period (Table 10
). The cows in the HCFS were heaviest at the beginning and end of the feed system period (524 and 583 kg, respectively), and had the highest BW gain over the 27-wk FS period (60 kg). The HSFS were lighter (P < 0.10) than the MPFS at wk 1 (508 and 514 kg, respectively), and at wk 26 (563 and 566 kg, respectively). The BW change did not differ between the HSFS and MPFS over the period (54.5 and 52.0 kg, respectively). The cows in the HCFS had higher BCS at the start and end of the FS period (2.90 and 3.03, respectively), compared with HSFS (2.82 and 2.87, respectively) and MPFS (2.82 and 2.83, respectively). The HCFS had the highest (P < 0.001) BCS gain over the period (0.15), with no difference between the MPFS (0.02) and HSFS (0.04) systems.
|
| DISCUSSION |
|---|
|
|
|---|
A genotype x environment interaction occurs when animals differ in their ability to perform in different environments (Falconer, 1990). A number of studies have shown a genotype x environment interaction in dairy production systems (Veerkamp et al., 1994; Boettcher et al., 2003). The existence of a genotype x environment interaction can be in the form of a "scaling effect," where the difference in the level of performance between 2 genotypes is smaller at a low level than a high feeding level, or a "re-ranking" where 2 genotypes rank differently in contrasting feeding systems. It is only possible to investigate the "scaling effect" with the present data set.
Kennedy et al. (2003) showed that with a high proportion of forage in the diet, high genetic merit dairy cows are not capable of eating much more than low genetic merit dairy cows, whereas on high concentrate diets, high genetic merit cows have the advantage of higher DM intakes. Hence, the differences in milk yield between genotypes are smaller in a high grass diet because intake is limited by constraining factors in the diet such as physical bulk, or by the relatively slow rate of intake that can be achieved by grazing cows compared with a system that entails a high level of concentrate supplementation.
Effect of Strain
The overall level of milk production achieved in the present study is similar to that obtained in previous studies at this center (Kennedy et al., 2002). Kennedy et al. (2002) showed total lactation milk yields of 6421 and 7389 kg with medium and high genetic merit HF cows, respectively, on low concentrate, and 7196 and 8461 kg, respectively, on a high concentrate system. These levels of production are much higher than that achieved in other pasture-based systems with high genetic merit cows. Kolver (2001) obtained 4678 kg of milk per cow with high genetic merit North American HF cows in a New Zealand grazing system with a low stocking rate. The milk protein composition achieved in the present experiment may be considered high for cows on a high forage system. Similar milk composition has been reported previously at this center (Buckley et al., 2000; Kennedy et al., 2002) when animals have been given access to high quality perennial ryegrass swards over a long grazing season.
The difference in observed milk production between the HP strain and the 2 low production strains (HD and NZ) is in agreement with the differences in predicted difference for milk, and similar to that observed in previous studies (Buckley et al., 2000; Boettcher et al., 2003). This higher milk production has been achieved successfully because milk production is a highly heritable selection criterion (Evans et al., 2002). The lower overall milk production and higher composition of the NZ strain compared with North American strain has also been observed previously (Kolver et al., 2000). Kolver (2001) showed that in New Zealand, North American-derived Holstein-Friesian cows were heavier, produced more milk volume and protein yield, and had lower concentration of fat and protein than New Zealand-derived Holstein-Friesian cows.
Until recently, the breeding schemes of the North American HF cow population within Europe have been based solely on increased milk production. Such breeding programs increase the gap between energy intake and output in early lactation. For this reason a greater mobilization of BW and BCS is expected among the HP strain. The higher BW of the HD strain relative to the HP animals is understandable given their lower level of milk production, coupled with their higher muscularity.
The reduced BCS displayed by the HP cows is also expected because many studies have shown lower BCS in animals selected for higher milk yield (Berry et al., 2003) and because the correlated response in feed intake in early lactation, due to selection for milk yield, can cover only 40 to 48% of the extra requirement (Veerkamp et al., 1994). Selection on milk yield alone is expected to increase the mobilization of body reserves in early lactation, the magnitude and duration of which is related to reduced health and fertility (Buckley et al., 2003). Many studies have shown that cows genetically superior for milk production tend to have genetically lower BCS throughout lactation and greater BCS change in early lactation than those of lower genetic merit (Buckley et al., 2000). Dechow et al. (2002) suggests that as cows are selected for higher milk production and become genetically thinner, the amount of body condition lost postpartum is likely to increase at any given level of BCS at parturition.
The lower BW and higher BCS of the NZ strain is a consequence of the negative weighting on milk volume and BW within the New Zealand Breeding Worth index (Harris, 1998). Previous studies in New Zealand showed that offspring of North American HF sires were heavier and taller than the New Zealand HF. Kolver et al. (2000) showed that North American HF heifers had a higher BW and they were unable to maintain acceptable body condition and BW on a pasture-based diet when compared with the New Zealand Holstein HF.
Effect of FS
The higher total lactation milk, SCM, fat, protein, and lactose yield achieved with the HCFS group is expected, given the large increase in energy supply with this FS. This is consistent with previous studies, which showed that total DMI increased with increasing proportion of concentrate in the diet (Bargo et al., 2002). The decline in milk fat content observed with HCFS was in agreement with previous findings (Bargo et al., 2002), and can be largely attributed to the reduced fiber content in the diet. The general increase in milk protein content on HCFS is likely attributable to increased energy intake and has been reported previously (Dillon and Crosse, 1997; Bargo et al., 2002).
The significant genotype x environment interaction observed for milk production over the complete lactation and during the feed system period agrees with similar farmlet studies (Kennedy et al., 2002), as well as large database studies (Boettcher et al., 2003). These interactions are explained by the large difference in response to additional energy intake between strains. The higher milk production response by the HP animals to a greater grass allowance in the MPFS demonstrates their greater nutritional stress, relative to the lower production strains. Differences in the response to additional energy input have been attributed to such factors as physiological stage, body condition, and milk production potential. As an animal approaches its maximum milk production potential, an increasing proportion of the energy supplied is retained as body lipid (Broster et al., 1985).
In the present study, the low response to concentrate supplementation from the NZ strain coupled with their higher BCS suggests that these animals achieved a greater proportion of their potential milk production in the MPFS than the HP strain. Coulon and Remond (1991) found that when energy requirements are not being met, concentrate supplementation appreciably influences the energy balance of the animal resulting in increased animal performance. Although body reserves usually buffer shortfalls in energy supply, prolonged mobilization is not possible due to innate metabolic adaptations within the animal (homeostasis). Consequently, high yielding cows are of lower body condition and display a more pronounced response to concentrate supplementation, as their energy demands for milk production are not fully met on a predominantly pasture-based diet. Many other studies using high genetic merit dairy cows for milk production have shown large responses to concentrate on pasture-based systems (Delaby et al., 2001; Kennedy et al., 2002).
The interaction of strain of HF and FS approached significance for BW change. This is mainly due to the lower BW gain of the HP animals in both the MPFS and HCFS. The lower BW gain among the HP strain over the feed system period may have arisen as a consequence of the higher milk production relative to the HD group during the feed system period. The energy intake resulting from increased grass allowance (in the MPFS) and additional concentrate supplementation (in the HCFS) resulted in a larger milk production response in these systems rather than in BW gain. The lack of agreement between BW and BCS change has been observed previously (Sutter and Beever, 2000). These authors concluded that BW change was an inadequate index of body tissue mobilization in early lactation as it might be biased by increased gut content or water repletion of body tissue.
Historically, genetic selection in dairy cattle has been within country. However, global selection can increase the rate of genetic progress in dairy cattle by up to 17% compared with within-country selection (Lohuis and Dekkers, 1998). However, such success has in turn raised new doubts relating to the problem of genotype x environment interaction, as exemplified in this study. In a seasonal pasture-based system, characteristics other than milk production, such as reproductive performance and animal health, are important (Veerkamp et al., 2002). Ultimately, the optimum cow for pasture-based systems can only be identified by combining all traits (production and health) of economic significance in a weighted index of economic merit and choosing sires at the top of this index. The results of this study show that accurate prediction of milk production performance must be based on knowledge of both dairy cow genotype and feed environment.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication August 20, 2004. Accepted for publication November 19, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. A. Macdonald, G. A. Verkerk, B. S. Thorrold, J. E. Pryce, J. W. Penno, L. R. McNaughton, L. J. Burton, J. A. S. Lancaster, J. H. Williamson, and C. W. Holmes A Comparison of Three Strains of Holstein-Friesian Grazed on Pasture and Managed Under Different Feed Allowances J Dairy Sci, April 1, 2008; 91(4): 1693 - 1707. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. McEvoy, E. Kennedy, J. P. Murphy, T. M. Boland, L. Delaby, and M. O'Donovan The Effect of Herbage Allowance and Concentrate Supplementation on Milk Production Performance and Dry Matter Intake of Spring-Calving Dairy Cows in Early Lactation J Dairy Sci, March 1, 2008; 91(3): 1258 - 1269. [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] |
||||
![]() |
E. S. Kolver, J. R. Roche, C. R. Burke, J. K. Kay, and P. W. Aspin Extending Lactation in Pasture-Based Dairy Cows: I. Genotype and Diet Effect on Milk and Reproduction J Dairy Sci, December 1, 2007; 90(12): 5518 - 5530. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Bytyqi, J. Odegard, H. Mehmeti, M. Vegara, and G. Klemetsdal Environmental Sensitivity of Milk Production in Extensive Environments: A Comparison of Simmental, Brown Swiss, and Tyrol Grey Using Random Regression Models J Dairy Sci, August 1, 2007; 90(8): 3883 - 3888. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Kennedy, M. O'Donovan, F. P. O'Mara, J. P. Murphy, and L. Delaby The Effect of Early-Lactation Feeding Strategy on the Lactation Performance of Spring-Calving Dairy Cows J Dairy Sci, June 1, 2007; 90(6): 3060 - 3070. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. McCarthy, D. P. Berry, P. Dillon, M. Rath, and B. Horan Influence of Holstein-Friesian Strain and Feed System on Body Weight and Body Condition Score Lactation Profiles J Dairy Sci, April 1, 2007; 90(4): 1859 - 1869. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. McCarthy, B. Horan, P. Dillon, P. O'Connor, M. Rath, and L. Shalloo Economic Comparison of Divergent Strains of Holstein-Friesian Cows in Various Pasture-Based Production Systems J Dairy Sci, March 1, 2007; 90(3): 1493 - 1505. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Roche, D. P. Berry, and E. S. Kolver Holstein-friesian strain and feed effects on milk production, body weight, and body condition score profiles in grazing dairy cows. J Dairy Sci, September 1, 2006; 89(9): 3532 - 3543. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||