|
|
||||||||
ebekAnimal Sciences Group, Wageningen University and Research Centre, Animal Production Division, 8200 AB Lelystad, the Netherlands
1 Corresponding author: bonne.beerda{at}wur.nl
| ABSTRACT |
|---|
|
|
|---|
Key Words: milk yield energy balance genotype by environment
| INTRODUCTION |
|---|
|
|
|---|
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
The heifers were sired by 53 different bulls and the size of paternal half-sibling groups was controlled to 6 individuals or less. The bloodlines of the heifers was 100% North American Holstein or mixed breed with 87.5% North American Holstein and 12.5% European Friesian-Holstein (n = 8). Cows were selected on the basis of the Dutch "Inet" EBV. Inet is the Dutch production index for milk, fat, and protein and is calculated as: 0.06 x EBV kg of milk + 0.7 x EBV kg of fat + 4.2 x EBV kg of protein; and is expressed in euros (
). The H cows had an average Inet value of 171, ranging from 112 to 241, and L cows had a mean Inet value of 24, ranging from 117 to 34. In a similar way, the mean EBV for kilograms of milk production were 771 (ranging from 6 to 1,704) and 68 (from 729 to 561) for H and L cows, respectively. The mean reliability of the EBV of the heifers sires was 95%.
The heifers were housed indoors in 4 adjacent sections of a free-stall barn at research farm Nij Bosma Zathe, the Netherlands. Each of the 4 sections, representing a combination of milking frequency and ration, consisted of 16 cubicles and 16 feeding troughs and housed both H and L heifers. Nonexperimental cows were added to each group up to a group total of 16 cows. The cows in the 2-times-daily groups were milked at 0600 and 1700 h, whereas those in the 3-times-daily groups were milked at 0600, 1400, and 2200 h. The free-stall floor was slatted concrete and cleaned by a manure scraper. Cubicles were equipped with rubber mats and a thin layer of sawdust. Water and partially mixed ration (PMR) were available ad libitum. The PMR was provided via a roughage intake control system, which allowed the monitoring of the daily roughage intake per individual. The roughage intake control system (Insentec, Marknesse, the Netherlands) comprises a feeding gate (which recognizes individual cows by their neck collars), and weighing troughs, which register the amount of PMR that cows eat during a meal. Concentrates were provided individually via feeding boxes and in the milking parlor. Cows in the E groups had access to a high caloric density PMR and 8 kg of concentrate per day (Table 1
). Those in the S groups had access to a low caloric density PMR and received 3 kg of concentrate per day. The composition of the ration for the E groups, which was well balanced and attuned to expected mean milk production, was adapted for the cows in the S groups in such a way that energy content would be reduced with minimal effects on other aspects of the ration. Obviously, the 2 rations differed for multiple characteristics (see Table 1
), but here these differences are referred to as variation in "caloric density."
Measurements
On an individual level, milk production was recorded every milking. Fat, protein, and lactose contents in milk were measured weekly in composite milk samplings from successive milkings over 2 d. Weighing platforms were used for automatically monitoring BW daily. Body condition was scored every 2 wk starting in lactation wk 1 by a single experienced farm assistant. Intake of roughage and concentrates was automatically registered per individual cow and used to calculate daily intakes. Additionally, a set of parameters was measured to assess blood metabolite and hormone profiles, locomotory function, udder health, and immune function, but these will be presented elsewhere.
Calculations
Milk production was expressed as ECM: (0.337 + 0.116 x fat content + 0.06 x protein content) x kg of milk. Energy balance (EB) was expressed as NEL in megajoules and calculated as: DMI per day x NEL in total ration {[(42.4 x BW0.75 + 442 x ECM) x (1 + (ECM 15) x 0.00165)] + 660} x 0.0069, where +660 is the correction for growth in heifers and the factor x0.0069 converts the Dutch energy unit VEM into megajoules. For information on the Dutch energy (VEM) system, see Van Es (1978). Protein balance (PB, in g/d) was calculated as intake of true protein digested in the small intestine (DVE) minus requirement DVE: {54 + [0.1 x BW]} + {1.396 x [(milk production in kg/d) x (milk protein in %) x 10]} + {0.000195 x [(milk production in kg/ d) x (milk protein in %) x 10]2}. For information on the Dutch protein (DVE/OEB) system, see Tamminga et al. (1994). Postpartum empty body weight (EBW) was expressed as BW (DMI x 4) to correct for gut fill (Jarrige, 1989).
Statistical Analyses
The data were analyzed with REML-type analyses by use of ASREML (Gilmour et al., 2004). Fixed effects in the model were milking frequency (2 or 3), caloric density of the ration (S or E), genetic merit (L or H), and interactions among these, week of lactation (1 to 14), and the time since the start of the experiment (wk 1 to 81). The latter functioned to minimize bias originating from time-related variations (e.g., seasonal fluctuations) and will not be discussed further. Time effects were fitted as splines and effects of lactation week were fitted per combination of milking frequency, caloric density of the ration and genetic merit. Thus, the model for statistical analyses was:
![]() |
where Yijklmn is the variable for milking frequency i, caloric density j, genetic merit k in lactation week m, in week l since the start of the experiment and animal n; MFi is the milking frequency (2 or 3 times a day); QRj is the caloric density of the feed ration (relatively high or low); GMk is the genetic merit (high or low genetic merit for fat and protein production); inijk represents all interactions between the 3 fixed effects; f(WSl) is a spline describing the effect of the weeks passed since the start of the experiment; f(WLijkm) is a spline describing the effect of the weeks passed since the start of the lactation; and ANn is the individual animal (included as a random effect).
Results presented here are the means (±SEM) predicted by the described statistical model. In the results section, effects and differences between groups are reported only if significant (P < 0.05), unless explicitly stated otherwise.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
195, on average. With an average genetic trend of about +
20/yr in the Netherlands, this represents a period of 10 yr of ongoing selection. The 2 rations differed considerably in energy content although they were within the range of normal rations used on Dutch farms, and caused a significant difference in DMI (4.2 kg/d). The differences in rations caused substantial differences in EB, PB, and milk production levels (9.6 kg of ECM/d). Milking 3 times daily instead of 2 times increased yield (4.1 kg of ECM/d), although like the stimulating effects of high genetic merit (3.5 kg of ECM/ d), this required the provision of high caloric density feed. The experimental factors caused mean daily ECM productions in the first 100 DIM to vary between 21.8 and 35.2 kg, which represents a substantial contrast as mean daily ECM production levels of Dutch Black and White dairy cows (across 305 d) were 19.5 and 28.8 kg in 1985 and 2004, respectively. The experimental factors affected milk production in the presumed direction. Milk productions were lowest in the L cows that were fed the S ration and milked 2 times daily and highest in the H cows that were fed the E ration and milked 3 times a day. The discriminative power of the study may be illustrated by the minimal difference that lead to significance for both main effects and GxE interactions. Power calculations (
= 0.05, ß = 0.10) resulted in minimal differences of about 10 and 20% of the mean ECM, respectively, with similar findings for EB. For ECM this means that contrasts between levels of main factors of about 3 kg/d or more were significant. For differences between combinations of genotype, milking frequency, and feed caloric density, this was about 6 kg/d or more and this relatively low power may have attributed to the absence of significant 3-way interactions.
Effects of Genotype
Significant effects of genetic merit for milk yield, independent of other factors, occurred for postpartum BCS with significantly lower scores in H cows than in L cows. A positive genetic correlation exists between milk production level and loss of body condition in early lactation (Berry et al., 2002). Changes in BCS reflect the energy status of cows and several studies associated high genetic merit with strong NEB (Veerkamp et al., 1995). Increased growth hormone concentrations in high genetic merit cows are likely to play a role (Veerkamp et al., 2003), for example, as growth hormone inhibits insulin-induced lipogenesis and thus influences the partitioning of resources in the body. Strikingly, the lower BCS in H cows than L cows was not accompanied by a stronger calculated NEB in H cows. Inaccuracies in calculating EB will have played a role (see next section), although a biological cause cannot be excluded. The H cows may have relatively high energy expenditures for maintenance and differ from L cows in the way food resources are utilized and allocated to body functions.
Effects of Environment
Ration.
Ration was the most determining factor affecting milk production, with decreased milk production when the caloric density was low (9.6 kg of ECM/ d). In part, this decrease in ECM was caused by significantly reduced protein content (0.26%). The low caloric density ration significantly reduced milk lactose (0.15%), but not fat content, resulting in an increased fat to protein ratio (+0.11). The results are not surprising as the provision of less concentrate typically decreases DMI, increases fat content, and decreases protein content (Brocard and Bareille, 2004). Higher energy intakes from supplementation with cereal grains, as for the cows in the E groups in this study, promotes ruminal propionate and microbial protein more than acetate and this increases the mammary synthesis of lactose and protein to a higher degree than fat (Walker et al., 2004). In addition to the energy supply, the direct protein supply will have affected protein content. Dry matter intake and calculated PB were significantly higher in cows in E groups than in S groups (+58.8 g/d), and increased intake of CP, at least relative to conditions of low protein intake, will increase protein content (Walker et al., 2004). The temporal patterns in milk production and content as shown in Figures 1
and 2
are in line with earlier studies reviewed by Walker et al. (2004). Because of dilution, fat and protein concentrations tend to be lowest around 40 to 60 d postpartum when lactose and milk productions are relatively high.
In line with the fat to protein ratios and calculated PB, the calculated EB were significantly lower in cows in the S groups than E groups (6.6 MJ/d). The impact of the feed treatment overruled the typical inverse relationship between milk yield and energy balance within a more or less homogeneous setting (Veerkamp et al., 1995). Others have reported that EB was more strongly controlled by feed intake than by milk production (Beam and Butler, 1999); this could apply especially to early lactation when cows have difficulty matching feed intake to energy expenditure. Calculated EB is typically most negative within the first 12 d postpartum and increases above zero after approximately 70 d postpartum (De Vries et al., 1999). The latter corresponds with EB becoming positive in lactation wk 11 in our cows in the E groups. Cows in the S groups typically remained in NEB throughout the 100-d experimental period and this marks these animals as having strong NEB compared with the normal situation.
The mean cumulative balances for energy and protein were 930 MJ/100 d and 3,720 g/100 d, respectively, and during the first 100 DIM the heifers lost, on average, 27.1 kg of EBW (as derived from Table 2
). This amounts to energy and protein values per kilogram of mobilized body reserves of 34.3 MJ and 137 g, respectively. Such values are within the wide range of expected values (Tamminga et al., 1997), although the calculated energy deficit is high relative to loss of EBW. Compare, for example, the 34.3 MJ/kg of mobilized body reserves (present study) with the 30.9 MJ/kg during the first 8 wk of lactation as reported by Tamminga et al. (1997). One unit of fat from body tissue contains more energy than one unit of protein from body tissue and the relatively little loss of EBW compared with calculated NEB indicates that the heifers mainly mobilized fat. The variation in the aforementioned values will in part result from inaccuracies. The systems used for calculating energy balances, like the Dutch VEM system (Van Es, 1978), are based on averages found for cows used in given feeding trials. Variation in feed composition and individual traits of cows, like their rumen flora and papillae, will cause inaccuracies in the estimation of the cows energy balance (Bruinenberg et al., 2002). Also, cows differ with regard to the partitioning of energy between different physiological systems; cows with similar energy intakes and expenditures via the milk may actually experience a different EB.
Milking Frequency.
Milking frequency affects the rate of cell death by apoptosis, like the degree of oxidative stress imposed by feed and the tissues ability to handle damage induced by reactive oxygen species (Stefanon et al., 2002). Milking more frequently may shift the balance between cell proliferation and cell removal in a way that increases the number of mammary cells and favors milk yield. For example, by local enhancement of prolactin receptors and by increased prolactin-mediated inhibition of local IGF binding protein-5 production, which would otherwise abolish the cell-saving actions of growth hormone-induced IGF-I (Flint et al., 2001). Thus, increased milking frequency promotes milk yield and persistency of lactation; for example, by counteracting the negative effects of oxidative stress. Here, we found that increased milking frequency increased milk yield, although only significantly when the cows received high caloric density rations. This indicates that in the S groups, nutrient intake was limiting for milk production. Milking the cows 3 times a day instead of 2 times significantly decreased EB. More-frequent milking lowered protein content (0.18%), fat content (0.23%), calculated EB (8 MJ/d), and calculated PB (96.3 g/d). Environmental factors like milking frequency and feeding management, more than milk production level per se, strongly predicted the cows energy balances.
Effects of GxE Interactions
Different genotypes are not always equally affected by different environments (Falconer, 1952), and across different environments such GxE interactions can result in different magnitudes of variance (scaling) as well as differences in ranking of the genotypes (reranking). Genotype x environment interactions in dairy cattle have been demonstrated at Langhill Dairy Cattle Research Centre in Scotland (Veerkamp et al., 1994) and Moorepark Research Centre in Ireland (Dillon et al., 2006). Especially under conditions that promote high-yield, high genetic merit cows have increased tissue mobilization, DMI, and, as a result, milk production (Veerkamp et al., 1995). This corresponds with our findings. High genetic merit compared with low genetic merit (mean difference of 3.5 kg of ECM/d), as well as 3 times milking compared with 2 times milking (mean difference of 4.1 kg of ECM/d), significantly increased milk production, but only in the groups that received high caloric density rations. The GxE interactions indicate that in the cows that received low caloric density rations, the low intake of resources limited the effects of factors that normally promote milk production. Similarly, Hoogendoorn et al. (1990) found that grazed cows showed milk responses to growth hormone only when the quality of the grass was relatively good. Earlier findings predict that our H cows should have had a stronger NEB than L cows, at least in the E groups (Veerkamp et al., 1995). This effect (6.5 MJ/d) proved to be near significance (P = 0.07). The expected higher DMI in H cows compared with L cows under yield-promoting conditions (i.e., in the E groups) was not found. The H cows did have higher DMI than the L cows, but only in the groups that were milked 2 times daily. This influence of milking frequency is difficult to explain. Interactions between milking frequency and genetic merit were detected also for postpartum BW, with L cows being heavier than H cows in the groups that were milked 3 times daily. This was near significance for EBW and seems to correspond with the relatively low BCS of H cows. The higher fat content in H cows than L cows, but only in the S groups, mirrors relatively high fat output in the H cows under severe conditions. Differences in fat content between cows in the HS and LS groups existed from lactation wk 1 and increased somewhat over time, indicating a tendency of high genetic merit animals to mobilize fat reserves for fat production, although diet-related variation in ruminal acetate production may have played a role. The relatively favorable calculated PB for cows in the EL groups means that under favorable conditions of sufficient nutrient intake, the balance between allocating protein to maintenance or milk production tips more toward the former in low than in high genetic merit cows. High genetic merit cows seem more prone to allocate resources toward high fat output in times of nutrient-deficient diets and toward high protein output with high nutrient diets. The relatively low milk productions by cows in the 2L groups compared with the 2H, 3L, and 3H groups indicates that a first yield-promoting measure, be it high genetic merit or an extra milking, renders the most marginal benefit compared with additional ones.
High Yield and Animal Health
High genetic merit for milk yield has been associated with poor energy balance, increased disease susceptibility (Pryce et al., 1998), and decreased fertility (Veerkamp et al., 2003), although not in all studies (Knight et al., 2004). The high genetic merit cows in this study were more prone to allocate resources toward milk production than low genetic cows and had lowered postpartum BCS. Low BCS and high dairy form (angularity) are genetically correlated with reduced cow health (Dechow et al., 2004) and strong decreases in postpartum BCS decrease fertility (Gillund et al., 2001). High genetic merit for milk yield seems intrinsically connected to increased risks for disorders related to severe NEB, although in this study management factors such as feed caloric density and milking frequency had the stronger impact on calculated EB and PB. Low nutrient input reduced energy balances, as did an extra milking, but high nutrient input may have negative side effects as well. Sehested et al. (2003) compared groups of cows fed a diet of clover grass supplemented with different amounts of concentrates (0, 19, or 38% of DMI), and the lowest incidences of health disorders were associated with the group that did not receive concentrates. Mismanagement can be a more potent risk factor for animal health than genetic ability for milk yield, although one-sided selection for high milk yield is likely to render cows that are increasingly more vulnerable to disease and poor fertility (Oltenacu and Algers, 2005). Keeping such animals by farmers that have increasingly less time available per cow because of scaling-up of farms, at least in the Netherlands, demands increasingly more management skills. In recent years, the tendency for one-sided selection on high yield has been weakened (Miglior et al., 2005) and this change is supported by the present findings.
Clear GxE interactions existed for milk production level, which increases the likelihood that such interactions play a role in the relationship between milk yield and animal health. Others reported on such interactions. Windig et al. (2005) found that cows in high-producing herds had better udder health than those in low-producing herds, but within the former, the highest producing cows had the poorest udder health. In situations of inferior management, unfavorable genetic correlations between milk yield and somatic cells scores or conception rate are relatively strong. From different studies, Oltenacu and Algers (2005) concluded that negative genetic correlations between production and fitness traits increase in less favorable environments and that selection for increased yield has decreased the adaptability of modern dairy cows. The present findings indicate a relatively strong tendency of high genetic merit cows to allocate resources toward milk production at the cost of maintenance, but there was little evidence that these animals are especially at risk for severe negative EB and PB during low nutrient supply. High genetic merit cows may have been selected to reach a certain limit in terms of BCS loss and NEB, but like low genetic merit cows, reduce yield in response to low nutrient supply beyond that point. The assumption that high genetic merit cows especially have difficulty in coping with extensive keeping conditions is not confirmed by the results presented here because, for example, they lower milk production levels in response to low nutrient intake in the same way as low genetic merit cows.
Signs for increased health risks, like severe negative EB and PB or low BCS, were not concentrated in the highest producing cows that were in the 3EH group or E groups. The production levels were apparently within the cows metabolic capacities, which more than their mammary glands capacities limits milk yield. High milk production levels per se need not compromise the health status of cows, although a negative correlation with EB exists and this is linked to minor changes in the hypothalamic-pituitary axis (Beerda et al., 2004). The latter indicates that high yield does represent allostatic load and triggers allostatic responses, but does not compromise the functioning of biological systems such as the immune system. In older cows, the situation may be different because of accumulated load.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication April 17, 2006. Accepted for publication September 5, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. Hammami, B. Rekik, H. Soyeurt, C. Bastin, J. Stoll, and N. Gengler Genotype x Environment Interaction for Milk Yield in Holsteins Using Luxembourg and Tunisian Populations J Dairy Sci, September 1, 2008; 91(9): 3661 - 3671. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. Windig, B. Beerda, and R. F. Veerkamp Relationship Between Milk Progesterone Profiles and Genetic Merit for Milk Production, Milking Frequency, and Feeding Regimen in Dairy Cattle J Dairy Sci, July 1, 2008; 91(7): 2874 - 2884. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. C. Friggens, P. Berg, P. Theilgaard, I. R. Korsgaard, K. L. Ingvartsen, P. Lovendahl, and J. Jensen Breed and Parity Effects on Energy Balance Profiles Through Lactation: Evidence of Genetically Driven Body Energy Change J Dairy Sci, November 1, 2007; 90(11): 5291 - 5305. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Ruiz-Sanchez, R. W. Blake, H. M. A. Castro-Gamez, F. Sanchez, H. H. Montaldo, and H. Castillo-Juarez Short Communication: Changes in the Association Between Milk Yield and Age at First Calving in Holstein Cows with Herd Environment Level for Milk Yield J Dairy Sci, October 1, 2007; 90(10): 4830 - 4834. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Ouweltjes, B. Beerda, J. J. Windig, M. P. L. Calus, and R. F. Veerkamp Effects of Management and Genetics on Udder Health and Milk Composition in Dairy Cows J Dairy Sci, January 1, 2007; 90(1): 229 - 238. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |