JDS
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Interpretive Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Beerda, B.
Right arrow Articles by Veerkamp, R. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Beerda, B.
Right arrow Articles by Veerkamp, R. F.
J. Dairy Sci. 90:219-228
© American Dairy Science Association, 2007.

Effects of Genotype by Environment Interactions on Milk Yield, Energy Balance, and Protein Balance

B. Beerda, W. Ouweltjes, L. B. J. Sebek, J. J. Windig and R. F. Veerkamp

Animal Sciences Group, Wageningen University and Research Centre, Animal Production Division, 8200 AB Lelystad, the Netherlands

1 Corresponding author: bonne.beerda{at}wur.nl


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Increases in genetic merit for milk yield are associated with increases in mobilization of body reserves. This study assessed the effects of genotype by environment (GxE) interactions on milk yield and energy and protein balances. Heifers (n = 100) with high or low genetic merit for milk yield were milked 2 or 3 times a day and received rations of low or high caloric density. The management factors were selected to induce substantial differences in milk production levels and model different management strategies. The 2 x 2 x 2 factorial arrangement enables the assessment of the effects of genotype, environment, and GxE interactions. Mean daily energy-corrected milk production in the first 100 d in milk varied between 21.8 and 35.2 kg among the groups. The experimental factors affected milk production in the presumed direction. Ration was the most determinant factor on milk production. Effects of milking frequency and genetic merit were significant only in the groups that were fed rations with high caloric density. Signs for severe negative energy balances, protein balances, and low body condition scores, all of which may be indicative of health risks, were not concentrated in the highest producing cows. Feed caloric density and milking frequency had stronger effects on energy balances and protein balances, with unfavorable effects of low caloric density feed and an extra milking. This emphasizes the possible effect of mismanagement on animal health risks. High genetic merit cows had significantly lower postpartum body condition scores. Genotype x environment interactions existed, but more information is needed to determine if cows of different genetic merit for milk yield are differently at risk for disease under specific conditions. High milk production levels per se will increase allostatic load, but need not compromise the health status of relatively young cows. Ongoing one-sided selection for high yield may be combined with good animal health, but because high genetic merit for milk yield seems intrinsically connected to the allocation of resources from maintenance toward milk, this puts increasing demands on farmers’ time and management skills.

Key Words: milk yield • energy balance • genotype by environment


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Milk production levels have increased around 2% per - year (since 1985) in the United States and the Netherlands (Dillon et al., 2006), in part because of selective breeding. Increases in genetic merit for milk yield go together with increases in feed intake, but the latter does not fully compensate the extra energy demands during early lactation, resulting in a more negative energy balance (NEB) and increased mobilization of - body reserves (van Arendonk et al., 1991). Severe NEB and the associated changes in metabolism play a central - role in the genesis of production diseases (Goff and Horst, 1997; Drackley, 1999). In the present study, the role of energy balances in yield–health relationships - is considered with emphasis on the roles of genotype, environment, and genotype x environment (GxE) interactions. The Danish Institute of Agricultural Sciences has conducted a long-term evaluation of the genetic factors breed and aptitude for milk yield, and TMR energy density on milk production, BW, BCS (Nielsen . et al., 2003), body fat, and lipolytic responses (Theilgaard et al., 2002). Results from this evaluation have demonstrated that genetic selection has resulted in . stronger lipolytic responses to energy deficiency in Danish Reds and Holsteins than in Jerseys. Feeding a TMR with higher energy density increased BCS, and breed x diet interactions exist for milk production. The present study compares Holstein-Friesian heifers with high or low genetic merit for milk, fat, and protein production . under 4 different conditions: combinations of milking frequency (2 or 3 times a day) and feed caloric density (high or low). The experimental factors were selected for their strong effects on milk production, thus causing substantial differences in milk production levels, and because they create environments that model different management strategies in the range from extensive to intensive. The experimental design allows discrimination between the roles of genotype, environment, and GxE interactions in relationships between milk production level and energy and protein balances. Here, the effects of GxE interactions on milk yield, energy balance, and protein balance are assessed and the results help to determine if potential health risks in dairy cows are especially linked to high yield per se, genetic merit for high yield, management, or a mismatch between the latter 2 factors.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Experimental Design
Treatments were applied in a 2 x 2 x 2 factorial arrangement. Heifers with a high genetic merit for fat and protein production (H) or a low genetic merit (L) were milked 2 times per day (2) or 3 times (3) and fed a relatively high caloric density ration (E, as derived from energy) or low caloric density ration (S, as derived from silage; Table 1Go). Prepartum and during the first 100 DIM, measurements were done to assess milk yield and energy and protein balances in H and L under different keeping conditions. Per combination of milking frequency and ration, 14 H cows were compared with 11 L cows. The slight imbalance was caused by difficulties in obtaining L heifers on farms that met our health requirements. Experimental groups were designated by frequency of milking (2 or 3), energy density of ration (E or S), and genetic merit group (H or L) to give 8 treatment designations: 2EH, 2EL, 2SH, 2SL, 3EH, 3EL, 3SH, and 3SL.


View this table:
[in this window]
[in a new window]

 
Table 1. The compositions of the 2 rations fed to cows1
 
Animals and Keeping Conditions
The heifers were purchased from 61 farms across the Netherlands. Measures were taken to reduce the risk of introducing novel pathogens into the research farm’s herd. Farms from which the heifers were purchased had demonstrably low risks for harboring paratuberculosis, leptospirosis, or infectious bovine rhinotracheitis, as evidenced by certificates from specific disease eradication and prevention programs. Blood samples were collected from candidate experimental cows and screened for possible infections of infectious bovine rhinotracheitis, bovine virus diarrhea, leptospirosis, and Salmonella. Heifers that showed signs of disease, including skin lesions and leg or claw disorders, were barred from the research farm. One hundred heifers were used in this study. They arrived at the research farm on average 10.6 ± 0.8 (SEM) wk before calving. The latter occurred at a mean age of 2.2 ± 0.02 (SEM) yr from May 2003 through August 2004.

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 ({euro}). 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 1Go). 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 1Go), 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:


Formula

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
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Predicted mean values during the first 14 wk of lactation for the different parameters are summarized in Table 2Go.


View this table:
[in this window]
[in a new window]

 
Table 2. Predicted means for all heifers and per experimental group1
 
Milk Production
Mean milk production varied between 21.8 ± 1.0 kg of ECM/d in group 2SL and 35.2 ± 0.9 kg of ECM/d in group 3EH (Table 2Go). Milk production levels were higher (P < 0.001) in cows in E groups (32.4 ± 0.5) than in those in S groups (22.8 ± 0.5). Milking frequency (P < 0.001) and genetic merit for milk yield (P < 0.01) affected milk production, but these effects depended on ration (P-values for 2-way interactions < 0.05). Only when cows were fed a high caloric density ration did 3 times milking (3E = 34.5 ± 0.7) increase milk production level compared with 2 times (2E = 30.4 ± 0.7), and did H cows (EH = 34.2 ± 0.7) have higher yields than L cows (EL = 30.7 ± 0.7). Interactions between milking frequency and genetic merit (P < 0.01) resulted in decreased milk production levels in the 2L groups compared with the 2H, 3L, and 3H groups. Milk production levels fluctuated between lactation weeks (P < 0.05), especially when the level of production was relatively high; for example, in the E groups (see Figure 1Go).


Figure 1
View larger version (16K):
[in this window]
[in a new window]

 
Figure 1. Mean ECM productions per group (n = 11 to 14) of cows of high (H) or low (L) genetic merit for milk yield that were milked 2 or 3 times daily and fed high (E) or low (S) caloric density rations. Milk yield was significantly affected by week of lactation, ration, milking frequency x ration, and genetic merit x ration.

 
Fat and Protein Concentrations
Milking 3 times compared with 2 times reduced (P < 0.05) fat content from 4.43 ± 0.08% to 4.20 ± 0.08%. Fat content was increased in H cows (4.46 ± 0.10%) compared with L cows (4.19 ± 0.11), but only for the cows that were fed a low caloric density ration (P-value for 2-way interaction < 0.05). Protein content in milk was reduced (P < 0.001) by milking 3 times a day from 3.20 ± 0.03% to 3.02 ± 0.03. Protein content was higher (P < 0.001) in the milk of cows in E groups (3.24 ± 0.03%) than in S groups (2.98 ± 0.03). Fat to protein ratios in milk were lower (P < 0.001) for cows in E groups (1.34 ± 0.02) than in S groups (1.45 ± 0.02). Ratios decreased gradually from 1.46 ± 0.02 in wk 1 to 1.38 ± 0.02 in wk 9 and thereafter stabilized around this value. Lactose content was higher (P < 0.001) in the milk of cows in E groups (4.78 ± 0.02%) than in S groups (4.63 ± 0.02). Fluctuations in milk content during lactation occurred for fat and protein (P-values < 0.001) and are illustrated in Figure 2Go.


Figure 2
View larger version (14K):
[in this window]
[in a new window]

 
Figure 2. Mean values (n = 100) of percentages of protein (•, left y-axis), fat ({square}, right y-axis), and lactose ({triangleup}, right y-axis) in milk during the first 14 wk of lactation. Changes in time were significant for protein and fat.

 
DMI
Dry matter intake was higher (P < 0.001) in cows in E groups (20.3 ± 0.3 kg/d) than in S groups (16.1 ± 0.3) and, during the first weeks of lactation, increased (P < 0.001) over time (Figure 3Go). Only within the groups that were milked 2 times daily did the H cows (18.7 ± 0.3 kg/d) have a higher (P-value for 2-way interaction < 0.05) DMI than the L cows (17.7 ± 0.4).


Figure 3
View larger version (12K):
[in this window]
[in a new window]

 
Figure 3. Mean (n = 50) DMI during the first 14 wk of lactation for heifers that were fed low caloric density rations (•) or high caloric density rations ({square}). Differences between the groups and changes in time were significant.

 
Energy and Protein Balance
Calculated EB were lower (P < 0.001) in cows that were milked 3 times (–13.3 ± 1.7 MJ/d) than in those that were milked 2 times (–5.3 ± 1.7 MJ/d). Cows in E groups had a mean EB of –6.0 ± 1.7 MJ/d, which was higher (P < 0.001) than the –12.6 ± 1.7 MJ/d in cows in S groups. In the E groups, but not the S groups, H cows tended (interaction genetic merit x ration: P = 0.07) to have a lower EB (–6.5 MJ/d) than the L cows. Like calculated EB, calculated PB were lower (P < 0.001) in cows that were milked 3 times (–85.4 ± 23.4 g/d) than in those that were milked 2 times (10.9 ± 23.5 g/d). Cows in E groups had higher (P < 0.01) PB (–7.8 ± 23.4 g/d) than cows in S groups (–66.6 ± 23.6 g/d). The effect of genetic merit (P < 0.01) on PB merely reflected that PB in the EL groups were on average 142 g/d higher (P-value for 2-way interaction < 0.01) than in the EH, SL, and SH groups. High caloric density rations improved PB, but only in L cows. Alternatively, the latter had a better PB than H cows, but only when a high caloric density ration was fed. Changes in calculated EB and PB over time (P-values < 0.001) are illustrated in Figure 4Go.


Figure 4
View larger version (12K):
[in this window]
[in a new window]

 
Figure 4. Mean values (n = 100) for calculated energy balances (EB, •, left y-axis) and calculated protein balances (PB, {square}, right y-axis) during the first 14 wk of lactation. Changes in time were significant for both EB and PB.

 
BCS
The mean postpartum BCS was 2.78 ± 0.04; with a lower (P < 0.01) mean score for H cows (2.71 ± 0.04) than for L cows (2.85 ± 0.05). The mean BCS prepartum, assessed on average 1.6 wk before calving, was 3.58 ± 0.03, with the same predicted scores for H and L cows. Postpartum EBW was on average 445.3 ± 4.2 kg and affected only by week of lactation (P < 0.001). Mean weights decreased from 476.8 ± 4.6 kg in wk 1 to 444.7 ± 4.3 in wk 7, after which weights increased gradually to 449.7 ± 4.5 kg in wk 14. Uncorrected postpartum BW was on average 518.1 ± 4.5 kg. Unlike EBW, uncorrected postpartum BW was not affected by week of lactation. Interactions between milking frequency and genetic merit, which were near significance for EBW (P = 0.07), affected BW (P < 0.05), with cows in 3L groups being heavier than those in 3H groups. Mean postpartum BW were on average 535.3 ± 8.5 and 508.3 ± 7.5 kg, respectively, compared with 511.7 ± 8.5 and 516.9 ± 7.7 kg for cows in the 2L and 2H groups.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Dairy cow health affects animal welfare, animal productivity, zoonotic risks, and international trade. This study investigated whether selection for milk yield affects the risk for negative energy and protein balances, which may be associated with health disorders. Heifers with high or low genetic merit for milk yield were exposed to experimental factors that, a priori, were assumed to create substantial differences in milk production levels and model different management strategies in the range from extensive to intensive. This allowed the assessment of the effects of genotype environment and GxE interactions in relationships between milk yield and animal health. The experimental conditions were successful in creating substantial contrasts among groups in terms of genotypes, environments and, subsequently, milk production levels. The difference in Inet EBV between the high and low genetic merit cows was {euro}195, on average. With an average genetic trend of about +{euro}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 ({alpha} = 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 1Go and 2Go 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 2Go). 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 tissue’s 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
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
The results support the view that high genetic merit for milk yield is intrinsically connected to increased risks for NEB-related disorders, although feed caloric density and milking frequency had the stronger impact on energy and protein balances, which underlines the impact that mismanagement can have on animal health risks. Ongoing one-sided selection for high yield can be combined with good animal health, but at the cost of increasing demands on farmers’ time and management skills. Genotype x environment interactions played a role in the present study, but specific health parameters need to be evaluated to determine if cows of different genetic merit for milk yield are differently at risk for disease under intensive and extensive management conditions. Present and earlier findings indicate that high milk production levels per se need not compromise the health status of cows. High yield will increase allostatic load, however, and it cannot be excluded that accumulated physiological load across multiple lactations compromises disease resistance in older cows. Prospective publications will deal with the specific health indictors that were measured during this study and link these to the basic findings reported here.


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
We are indebted to Hans Hopster, Gert van Duinkerken and Sikko Beukema for their parts in the preparation of the study. The support and assistance during the experiment by all staff of research farm Nij Bosma Zathe is highly appreciated and we want to especially thank Johannes Postma, Jan Zonderland, and Eltjo van Marum. Henk Sulkers and Yvonne van Hierden played a major role in performing this experiment and we are grateful for their significant contribution. Arie Klop assisted us in feed-related matters. Many students participated in the study and we are thankful for the help of Leonie Jelsma, Jeroen Nolles, Daan Jansen, Alet Wiersma, Nina Düffels, Gooitzen van der Veen, Jolanda Veldman, Anja de Wit, Sanne Bruining, Gerben Malda, Mirjam Wiersma, and Arno Bakker. This work was funded by the Dutch Ministry of Agriculture, Nature and Food Quality (LNV).

Received for publication April 17, 2006. Accepted for publication September 5, 2006.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 


Beam, S. W., and W. R. Butler. 1999. Effects of energy balance on follicular development and first ovulation in postpartum dairy cows. J. Reprod. Fertil. Suppl. 54:411–424.[Medline]

Beerda, B., J. E. Kornalijnslijper, J. T. van der Werf, E. N. Noordhuizen-Stassen, and H. Hopster. 2004. Effects of milk production capacity and metabolic status on HPA function in early postpartum dairy cows. J. Dairy Sci. 87:2094–2102.[Abstract/Free Full Text]

Berry, D. P., F. Buckley, P. Dillon, R. D. Evans, M. Rath, and R. F. Veerkamp. 2002. Genetic parameters for level and change of body condition score and body weight in dairy cows. J. Dairy Sci. 85:2030–2039.[Abstract/Free Full Text]

Brocard, V., and N. Bareille. 2004. Dairy cows health and production performances according to the concentrate level. Page 143 in Proc. 55th Annu. Mtg. Eur. Assoc. Anim. Prod., Bled, Slovenia.

Bruinenberg, M. H., Y. van der Honing, R. E. Agnew, T. Yan, A. M. van Vuuren, and H. Valk. 2002. Energy metabolism of dairy cows fed on grass. Livest. Prod. Sci. 75:117–128.

De Vries, M. J., S. van der Beek, L. Kaal-Lansbergen, W. Ouweltjes, and J. B. M. Wilmink. 1999. Modeling of energy balance in early lactation and the effect of energy deficits in early lactation on first detected estrus postpartum in dairy cows. J. Dairy Sci. 82:1927–1934.[Abstract]

Dechow, C. D., G. W. Rogers, U. Sander-Nielsen, L. Klei, T. J. Lawlor, J. S. Clay, A. E. Freeman, G. Abdel-Azim, A. Kuck, and S. Schnell. 2004. Correlations among body condition scores from various sources, dairy form, and cow health from the United States and Denmark. J. Dairy Sci. 87:3526–3533.[Abstract/Free Full Text]

Dillon, P., D. P. Berry, R. D. Evans, F. Buckley, and B. Horan. 2006. Consequences of genetic selection for increased milk production in European seasonal pasture-based systems of milk production. Livest. Sci. 99:141–158.

Drackley, J. K. 1999. Biology of dairy cows during the transition period: The final frontier? J. Dairy Sci. 82:2259–2273.[Abstract]

Falconer, D. S. 1952. The problem of environment and selection. Am. Nat. 86:293–298.

Flint, D. J., E. Tonner, C. H. Knight, C. B. A. Whitelaw, J. Webster, M. Barber, and G. Allan. 2001. Control of mammary involution by insulin-like growth factor binding proteins: Role of prolactin. Livest. Prod. Sci. 70:115–120.

Gillund, P., O. Reksen, Y. T. Grohn, and K. Karlberg. 2001. Body condition related to ketosis and reproductive performance in Norwegian dairy cows. J. Dairy Sci. 84:1390–1396.[Abstract]

Gilmour, A. R., B. R. Cullis, S. J. Welham, and R. Thompson. 2004. ASReml Reference Manual. 2nd edition (draft 1.62). NSW Department of Primary Industries, Locked Bag, Orange, NSW 2800, Australia.

Goff, J. P., and R. L. Horst. 1997. Physiological changes at parturition and their relationship to metabolic disorders. J. Dairy Sci. 80:1260–1268.[Abstract]

Hoogendoorn, C. J., S. N. McCutcheon, G. A. Lynch, B. W. Wickham, and A. K. H. Macgibbon. 1990. Production responses of New Zealand Friesian cows at pasture to exogenous recombinantly derived bovine somatotropin. Anim. Prod. 51:431–439.

Jarrige, R. 1989. Ruminant Nutrition: Recommended allowances and feed tables. Page 389. INRA Publications, Paris, John Libbey Eurotext, London, UK

Knight, C. H., M. A. Alamer, A. Sorensen, I. M. Nevison, D. J. Flint, and R. G. Vernon. 2004. Metabolic safety margins do not differ between cows of high and low genetic merit for milk production. J. Dairy Res. 71:141–153.[Medline]

Miglior, F., B. L. Muir, and B. J. Van Doormaal. 2005. Selection indices in Holstein cattle of various countries. J. Dairy Sci. 88:1255–1263.[Abstract/Free Full Text]

Nielsen, H. M., N. C. Friggens, P. Løvendahl, J. Jensen, and K. L. Ingvartsen. 2003. Influence of breed, parity, and stage of lactation on lactational performance and relationship between body fatness and live weight. Livest. Prod. Sci. 79:119–133.

Oltenacu, P. A., and B. Algers. 2005. Selection for increased production and the welfare of dairy cows: Are new breeding goals needed? Ambio 34:311–315.[Medline]

Pryce, J. E., R. J. Esslemont, R. Thompson, R. F. Veerkamp, M. A. Kossaibati, and G. Simm. 1998. Estimation of genetic parameters using health, fertility and production data from a management recording system for dairy cattle. Anim. Sci. 66:577–584.

Sehested, J., T. Kristensen, and K. Soegaard. 2003. Effect of concentrate supplementation level on production, health and efficiency in an organic dairy herd. Livest. Prod. Sci. 80:153–165.

Stefanon, B., M. Colitti, G. Gabai, C. H. Knight, and C. J. Wilde. 2002. Mammary apoptosis and lactation persistency in dairy animals. J. Dairy Res. 69:37–52.[Medline]

Tamminga, S., P. A. Luteijn, and R. G. M. Meijer. 1997. Changes in composition and energy content of liveweight loss in daily cows with time after parturition. Livest. Prod. Sci. 52:31–38.

Tamminga, S., W. M. Vanstraalen, A. P. J. Subnel, R. G. M. Meijer, A. Steg, C. J. G. Wever, and M. C. Blok. 1994. The Dutch protein evaluation system–The DVE/OEB system. Livest. Prod. Sci. 40:139–155.

Theilgaard, P., N. C. Friggens, K. H. Sloth, and K. L. Ingvartsen. 2002. The effect of breed, parity and body fatness on the lipolytic response of dairy cows. Anim. Sci. 75:209–219.

van Arendonk, J. A. M., G. J. Nieuwhof, H. Vos, and S. Korver. 1991. Genetic aspects of feed intake and efficiency in lactating dairy heifers. Livest. Prod. Sci. 29:263–275.

Van Es, A. J. H. 1978. Feed evaluation for ruminants: 1. The systems in use from May 1977 onwards in The Netherlands. Livest. Prod. Sci. 5:331–345.

Veerkamp, R. F., B. Beerda, and T. van der Lende. 2003. Effects of genetic selection for milk yield on energy balance, levels of hormones, and metabolites in lactating cattle, and possible links to reduced fertility. Livest. Prod. Sci. 83:257–275.

Veerkamp, R. F., G. Simm, and J. D. Oldham. 1994. Effects of interaction between genotype and feeding system on milk production, feed intake, efficiency and body tissue mobilization in dairy cows. Livest. Prod. Sci. 39:229–241.

Veerkamp, R. F., G. Simm, and J. D. Oldham. 1995. Genotype by environment interactions: Experience from Langhill. In Breeding and Feeding the High Genetic Merit Dairy Cow. Br. Soc. Anim. Sci. Occ. Publ. 19:59–66.

Walker, G. P., F. R. Dunshea, and P. T. Doyle. 2004. Effects of nutrition and management on the production and composition of milk fat and protein: A review. Aust. J. Agric. 55:1009–1028.

Windig, J. J., M. P. L. Calus, and R. F. Veerkamp. 2005. Influence of herd environment on health and fertility and their relationship with milk production. J. Dairy Sci. 88:335–347.[Abstract/Free Full Text]


This article has been cited by other articles:


Home page
J DAIRY SCIHome page
H. Hammami, B. Rekik, C. Bastin, H. Soyeurt, J. Bormann, J. Stoll, and N. Gengler
Environmental sensitivity for milk yield in Luxembourg and Tunisian Holsteins by herd management level
J Dairy Sci, September 1, 2009; 92(9): 4604 - 4612.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
G. Bobe, G. L. Lindberg, L. F. Reutzel, and M. D. Hanigan
Effects of lipid supplementation on the yield and composition of milk from cows with different {beta}-lactoglobulin phenotypes
J Dairy Sci, January 1, 2009; 92(1): 197 - 203.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
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]


Home page
J DAIRY SCIHome page
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]


Home page
J DAIRY SCIHome page
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]


Home page
J DAIRY SCIHome page
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]


Home page
J DAIRY SCIHome page
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]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Interpretive Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Beerda, B.
Right arrow Articles by Veerkamp, R. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Beerda, B.
Right arrow Articles by Veerkamp, R. F.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS