|
|
||||||||

* Department of Animal Production, and
Department of Clinics, Faculty of Veterinary Medicine, Aristotle University of Thessaloniki, Greece
1 Corresponding author: banos{at}vet.auth.gr
| ABSTRACT |
|---|
|
|
|---|
Key Words: body energy dairy cow energy balance genetic profile
| INTRODUCTION |
|---|
|
|
|---|
The direct calculation of cow energy balance from estimates of feed intake and milk yield is quite difficult under field conditions, usually because feed intake records at the individual cow level are mostly unavailable. Thus, the use of various energy balance indicator traits has been proposed for the study of negative energy balance and its effect on health and fertility.
Evaluation of body condition is a subjective yet reliable and widely accepted method of the assessment of a cows fat tissue and body mass reserves (Edmonson et al., 1989; Fox et al., 1999) and, therefore, BCS can be viewed as an indicator trait to energy balance. In addition, total body energy content (EC) and cumulative effective energy balance (CEEB), calculated from BCS and live weight records, have been proposed for the study of energy balance (Banos et al., 2006). Furthermore, various metabolic traits have been associated with energy balance. Blood serum concentrations of glucose, BHBA, and NEFA are examples of such traits (Reist et al., 2002; Clark et al., 2005).
Before genetic selection for any of the above mentioned energy balance-related traits becomes feasible, their genetic variation and heritability must be studied. Body condition score has already been genetically analyzed in several studies (e.g., Coffey et al., 2001; Berry et al., 2003; Banos et al., 2005). Banos et al. (2006) also reported genetic parameters for EC and CEEB derived from research farm data. However, to the best of our knowledge, the genetic profile of blood serum glucose, BHBA, and NEFA throughout a lactation period, or even for a part of this period, has not been reported.
The objective of this study was to obtain a deeper understanding of the genetic variation in energy balance and metabolic regulation during lactation by a) characterizing the systematic change of BCS, EC, CEEB, and blood serum glucose, BHBA, and NEFA concentrations just before first calving and throughout the first lactation of cows kept on a commercial farm setting and b) estimating genetic parameters for these traits.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Cows calved at an average age of 27.6 ± 2.6 mo. These animals were either born on the farm or had been imported as pregnant heifers from 3 other European countries. The latter is a common practice among many commercial dairy farms in Greece.
Body condition score of these cows was assessed after the morning milking, by the same trained veterinarian, weekly from calving to wk 13 of lactation and thereafter monthly until the end of a 305-d lactation. A 5-point scale (1 = emaciated, 5 = obese, scored in 0.25-point intervals) and the method described by Ferguson et al. (1994) were used. At the same time, the cows live weight was estimated using a heart girth tape.
Blood was drawn from the coccygeal vein or artery in a subset of 365 randomly selected cows. Blood samples were left to clot at room temperature for approximately 30 min and then centrifuged at 2,000 x g. The obtained serum samples were stored at –20°C until analyzed for glucose, BHBA, and NEFA. Serum glucose and NEFA were assayed colorimetrically using commercial kits (Glucose GOD-PAP method, P. Zafiropoulos S.A., Attiki, Greece, and Wako NEFA C kit, Wako Chemicals GmbH, Neuss, Germany, respectively). The serum concentration of BHBA was assayed with the use of an enzymatic kinetic method based on the oxidation of β-hydroxybutyrate to acetoacetate by β-hydroxybutyrate dehydrogenase (Bruss, 1997):
![]() |
The intraassay coefficient of variation was 1 to 3% for glucose and 1.5 to 3% for NEFA whereas the interassay coefficient was 3 to 5% for both glucose and NEFA. For BHBA, the intra- and interassay coefficients were 2 to 4% and 4 to 8%, respectively.
The final data set consisted of 8,094 BCS, 8,087 estimated live weight, and 6,015 blood serum glucose, BHBA, and NEFA concentration records. Some cows did not complete the 305-d lactation because they were involuntarily culled. However, all animals in the data set had at least one BCS and estimated live weight record. Body condition score and estimated live weight records were used for the calculation of body EC and CEEB according to the procedure described in detail by Banos et al. (2006). Briefly, empty body weight and total lipid and protein weights were predicted using the method of the National Research Council (NRC, 2001) and were then combined to calculate EC. The latter was an estimate of the total energy in a cows body at any given time of lactation. Changes in predicted lipid and protein weights from one week of lactation to the next were converted to effective energy as described by Banos et al. (2006), yielding estimates of CEEB. The latter represented body energy changes as accumulated throughout lactation.
Single measurements of BCS, estimated live weight, EC, and serum glucose, BHBA, and NEFA concentrations taken approximately 2 mo before first calving were also available for a subgroup of the studied animals.
Statistical Analysis
Repeated cow BCS, EC, CEEB, glucose, BHBA, and NEFA records taken throughout first lactation were analyzed with the following random regression model; each trait was analyzed separately:
![]() |
where Yijkmn = record of cow k in week of lactation m, YSi = fixed effect of year-season of calving i (4 levels), Cj = fixed effect of country of origin j (4 levels), a1 = linear regression coefficients on age at calving (age), Wm = week of lactation m (44 levels), bn = fixed regression coefficient on week of lactation, ckn = random regression coefficient on week of lactation associated with the genetic effect of cow k including the full animal pedigree relationship matrix, dkn = random regression coefficient on week of lactation associated with the permanent environment effect of cow k, Pn = nth orthogonal polynomial of week m (n + 1 = order of polynomial), and eijkmn = random residual term (4 classes as described below).
In the model the fixed regression coefficient was associated with an overall, average lactation curve for each trait, whereas the random regression was associated with each individual cows deviation from the overall curve. The full model was applied to the analysis of BCS, EC, and CEEB. For the analysis of blood serum glucose, BHBA, and NEFA concentration, the effects of country of origin and age at calving were not included because preliminary analyses showed that they were not significantly different from zero (P > 0.05).
Fourth-order orthogonal polynomial was chosen to model the genetic effect following preliminary analyses of different (lower) orders and comparisons between models using the log-likelihood, until the latter stopped changing significantly (P > 0.05) with increasing polynomial order. This was true for all traits except BHBA, where third-order was the maximum significant.
Permanent environment was modeled with second-order orthogonal polynomial. Efforts to increase the order were unsuccessful as they led to convergence problems. This may be a function of the data size, although the number of observations is not deemed too small for a controlled study. Another explanation can be that higher order regression might have actually over fitted first-lactation permanent environment effects. For the analysis of CEEB, glucose, and NEFA, a first-order permanent environment polynomial was the highest significant.
Depending on lactation stage, 4 measurement error classes were defined as follows: wk 1 to 4, 5 to 9, 10 to 21, and >21. Different residual variances were estimated for each measurement error class, while covariances between classes were zero.
The REML estimates of (co)variance components from the model were used to calculate heritabilities for each trait and week of lactation as well as genetic correlations between different weeks of lactation. All analyses were performed with the use of the ASREML software package (Gilmour et al., 2002). Single records, taken on heifers before calving, were analyzed with a similar model that included the effect of days to calving but excluded the fixed and random regressions on week of lactation. Product-moment correlations among estimated animal breeding values for all studied traits were calculated to derive estimates of the relationship between these traits.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
Estimated genetic variance of serum glucose levels was greater at the beginning and the end of lactation (Figure 4
). All estimates were significantly different from zero (P < 0.05) while heritability ranged from 0.12 (±0.04) on wk 38 to 0.39 (±0.11) at the end of lactation (Figure 5
). In this particular trait, there was a sharp increase of the genetic variance estimate toward the very end of the trajectory, which may be due to fitting a high-order polynomial. Estimates at the end of lactation were associated with relatively greater standard errors (0.003–0.004) compared with other lactation stages (0.001–0.002) but were still significantly different from zero (P < 0.05). Despite this observation, fitting a fourth-order polynomial was justified by the significant (P < 0.05) increase in the log-likelihood.
Genetic variance estimates for BHBA serum concentration were significantly greater than zero (P < 0.05) during the first 7 wk of lactation (Figure 4
). Heritability estimates for this period ranged from 0.08 (±0.04) to 0.40 (±0.06) (Figure 5
).
Estimated genetic variance of serum concentration of NEFA was significantly different from zero (P < 0.05) until wk 24 of lactation (Figure 4
). Heritability estimates for this period ranged from 0.08 (±0.04) to 0.35 (±0.05) in the beginning of lactation (Figure 5
).
Genetic correlations of the above traits in different weeks of lactation were near unity for adjacent weeks and decreased for weeks further apart. Compared with the onset of lactation, genetic correlations decreased most rapidly for BHBA, followed by glucose, NEFA, and the body energy traits. This is illustrated in Figure 6
presenting genetic correlations of each studied trait measured on wk 4 of lactation with subsequent week measurements. Average genetic correlation of body energy traits measured in the first 4 wk and last 4 wk of lactation was 0.25 and 0.10 for BCS and EC, respectively. Also, the average genetic correlation between the first 4 wk and all remaining weeks for the same 2 traits were 0.42 and 0.34, respectively. For CEEB, where genetic variation in the first 6 wk was not statistically significant (P > 0.05), genetic correlations of records in wk 7 to 10 with all remaining and the last 4 wk of lactation were 0.71 and 0.56, respectively. Genetic correlations of the blood serum concentration traits in different weeks of lactation were, in general, lower than those of the other traits. These correlations became practically zero for measurements taken more than 3 to 4 mo apart. For example, glucose measurements in the first 4 wk of lactation had significant (P < 0.05) positive genetic correlations with measurements up to wk 16 with an average estimate of 0.54. Thereafter, genetic correlations were not significantly different from zero. For BHBA, serum levels in the first 4 wk were significantly correlated with measurements until wk 9 (average estimate 0.72) and for NEFA the corresponding values were 15 wk and an average genetic correlation of 0.58.
|
Finally, average product-moment correlations among animal breeding values for the studied traits are presented in Table 3
. These are averages of correlations pertaining to the first 3 mo (13 wk) of lactation.
|
| DISCUSSION |
|---|
|
|
|---|
The genetic profile of BCS has been investigated in several previous studies. Our results are quite similar to those from studies conducted on single farms using observations on a few hundreds cows. For example, Coffey et al. (2001) reported heritability estimates for BCS across the first lactation that ranged from 0.38 to 0.81, using data from an experimental research station. In another study with similar data, Veerkamp and Brotherstone (1997) found that the heritability estimate was 0.43 for the first 26 wk of the first lactation. In the present study, data were from a commercial rather than a research farm, and the fact that results were consistent with the literature suggests that the type of farm may not be the crucial factor in a genetic analysis, as long as all other sources of systematic variation are properly recorded and accounted for.
Lower BCS heritability estimates have, in general, been reported in studies conducted using large numbers of cows raised in many different herds. This is expected because, in single farm studies, more controlled environmental conditions can lead to greater heritability estimates. Body condition score heritability estimates from large-scale studies have been reported to range from 0.23 to 0.28 (Jones et al., 1999), 0.39 to 0.51 (Berry et al., 2003), 0.22 to 0.38 (Banos et al., 2004), and 0.23 to 0.37 (Koenen et al., 2001).
Genetic variance and heritability estimates for EC and CEEB were very similar to those reported by Banos et al. (2006) using data from an experimental farm in Scotland. Among body energy traits measured in pregnant heifers, only BCS had a significant heritability estimate, suggesting possible utility in a genetic evaluation and selection program. Its genetic correlation with postpartum energy balance levels still needs to be established.
As has already been mentioned, genetic parameters for glucose, BHBA, and NEFA serum concentrations throughout lactation were not found in the literature. Therefore, no direct comparisons can be made with the results of this study. Ingvartsen and Friggens (2005) reported significant between-cow variation for glucose, BHBA, and NEFA concentrations. Although not a genetic analysis, this was an indication of genetic variation in these traits. Furthermore, genetic selection for milk yield has been found to affect these traits (Veerkamp et al., 2003), which is another indication of the existence of genetic variation. Such indications were confirmed and quantified in the present study. More comparable to our results were heritability estimates of blood serum glucose and NEFA concentration that were reported by Hayhurst et al. (2007) using data on Holstein male calves that were approximately 9 mo old. These estimates were 0.15 for both traits and were close to the lower limit of our results. The greater heritability values that were observed in our study during the first week of lactation are probably related to the metabolic challenges a dairy cow faces immediately after calving that may be connected to a different expression of her genetic component for the specific traits. This could also explain the greater heritability estimates for BHBA levels that were observed in our study during the first week of lactation.
Heritability estimates of blood metabolic traits in pregnant heifers were not statistically significant. From these estimates, only glucose approached significance (P = 0.08) and was slightly greater than the estimate reported for male calves in the study by Hayhurst et al. (2007).
Genetic correlations of the body energy traits studied here (BCS, EC, CEEB) in different weeks of lactation were very high for adjacent weeks and decreased for weeks further apart. This pattern was expected and consistent with that reported by Banos et al. (2005) for a trait associated with body energy balance accumulated throughout lactation. These results indicate that the way cows start accumulating lipid and protein changes at the first stages of lactation sets the pace for energy changes in the remainder of lactation. The utility of this result lies in the ability to predict future energy balance indicators, especially CEEB, from values obtained early in lactation, potentially assisting farm management practices. Furthermore, it could enable the estimation of CEEB at the first stages of lactation from values obtained in mid to late lactation.
Genetic correlations of the blood serum concentration traits in different weeks of lactation were, in general, lower than those of the other traits. These correlations became practically zero for measurements taken more than 3 to 4 mo apart. Results suggest that although early-lactation body energy measurements are decent genetic predictors of late lactation, this is not necessarily true for blood serum concentration traits, whose predictive capacity ends 11 to 16 wk later. For BHBA and NEFA, this is consistent with near-zero genetic variance levels observed at mid to late lactation. However, it is important to note that a genetic evaluation for the early-lactation profile of these traits is possible with a single measurement obtained at any time during the first 2 to 3 mo of lactation.
Product-moment correlations among animal breeding values for the studied traits were calculated as the minimum values of the genetic correlation between these traits. Average correlation estimates of the first 3 mo of lactation ranged from weak to strong. This early stage of lactation is the period when weekly phenotypic measurements were available. This period is also of interest to farm management because it precedes the onset of inseminations. Body condition score had a considerable correlation with the other body energy traits (EC and CEEB) but its correlation with metabolic traits was very low or close to zero. This result may have implications in the use of such traits for selection aiming at improving cow health and fertility. In the presence of significant genetic correlations with the latter, measuring only BCS will not provide as much information as measurement of more than one of the studied traits.
Previous studies on the genetic profile of BCS, EC, and CEEB, and their genetic correlation with economically important traits such as fertility or somatic cell count, have already alluded to the opportunity for genetic selection of dairy cows for these traits (Pryce et al., 2000; Banos et al., 2004, 2006). Results presented here support this argument. Furthermore, new selection opportunities arise with the genetic characterization of metabolic traits related to dairy cow energy balance. Genetic selection of cows that have more desirable metabolic profiles, especially at the critical early postpartum period, seems to be possible because the heritabilities of glucose, NEFA, and BHBA concentration, at least during the first weeks of lactation, were found to be significantly greater than zero. Heritability estimates of BCS before first calving suggest that the genetic evaluation of heifers could be possible, giving useful information early in the animals life. Validation of these data with similar studies that could also involve observations on multiparous cows would be desirable, while estimates of genetic correlations between these metabolic traits and health and fertility traits are largely missing from the literature. On the other hand, phenotypic relationships of metabolic energy balance indicators with health (LeBlanc et al., 2005; Hammon et al., 2006) and fertility (Reist et al., 2000; Westwood et al., 2002; Reist et al., 2003; Walsh et al., 2007; Wathes et al., 2007b) have already been reported. The substantial genetic variation presented here for these metabolic traits enhances the importance of an attempt to quantify the proportion of these phenotypic relationships that is due to genetic factors. If significant genetic correlations of blood metabolic traits with health and fertility are established, then measurements of the former can be used to facilitate genetic improvement of health and fertility traits, which are known to have low heritability.
Arguably, recording blood metabolic traits in a commercial dairy cattle population can be difficult. It requires a properly equipped laboratory, and the cost for consumables is not negligible. In our case, the consumables cost for the measurement of glucose, BHBA, and NEFA was approximately
3.5 (US$5.20) and analysis took 7 to 9 min per sample. An additional cost is required for the acquisition of blood samples. In this regard, possible first applications of such trait recording and monitoring may be restricted to selected herds, where elite breeding animals are principally identified and evaluated.
At the same time, metabolic profiles during the transition period are considered to be very useful tools for the refinement of nutrition and health management procedures, especially in high production herds (Herdt, 2000; Oetzel, 2004). This fact, together with the development of more automated techniques for frequent sampling and biochemical analyses and probably the reduction of the cost of consumables, could present a possibility for the introduction of metabolic profiles in dairy herd improvement schemes and the use of these measurements for management and selection purposes.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication December 20, 2007. Accepted for publication April 21, 2008.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |