|
|
||||||||
1 Statistical Animal Genetics Group, Institute of Animal Sciences, Swiss Federal Institute of Technology, ETH Zentrum CH-8092 Zurich, Switzerland
2 Holstein Association of Switzerland, Grangeneuve CH-1725 Posieux, Switzerland
Corresponding author: Y. de Haas; e-mail: Yvette.deHaas{at}inw.agrl.ethz.ch.
| ABSTRACT |
|---|
|
|
|---|
Key Words: production increase maturity rate functional trait conformation
Abbreviation key: DFS = days to first service, ECM = energy-corrected milk, ECMI = production increase for energy-corrected milk, ECMIij = production increase for energy-corrected milk from lactation i to j, ID = identification number, MilkI = production increase for milk, MilkIij = production increase for milk from lactation i to j, NRR = nonreturn rate, ProdI = production increase.
| INTRODUCTION |
|---|
|
|
|---|
Based on results of the conformation evaluation, adulthood of a dairy cow can be assumed to be reached at 60 mo of age (S. Wegmann, personal communication, 2004), but cows first calve at a much younger age. Depending on the housing conditions (varying from intensive rearing to extensive alpine pasture conditions in Switzerland), age at first calving varies between 18 and 36 mo. Therefore, most of the milking heifers are still growing during their first lactation and part of the energy and protein intake will be assigned to this purpose (Veerkamp, 2000). Milk production of most dairy cows will increase from first to second and from second to third lactation as the cow grows close to mature size and is able to use more nutrients for milk production. However, the amount of increase is different among cows. There is genetic variance for maturity rate (Krogmeier et al., 2003), i.e., some animals have a genetic predisposition to produce moderately in first lactation and to increase production levels in second and third lactations.
Some breeders try to breed cows with moderate first lactation production and high yields in subsequent lactations so as to reach the full production potential at 60 mo of age. In breeders circles, this gradual production increase (ProdI), a measure of the maturity rate, is believed to be found in trouble-free cows. Geneticists became interested in this new trait and tried to define it in a way that would enable selection. Two definitions have been proposed: 1) the definition used in Switzer-land, based on the EBV for each of the first 3 lactations (Schleppi and Bigler, 2002) and 2) the use of real lactation yields of the cows instead of EBV (Krogmeier et al., 2003). However, no study has estimated genetic parameters for ProdI computed with actual lactation records. The hypothesis of possible links of ProdI with functional traits at a genetic or phenotypic level has never been investigated. Because health and fertility are important components of productive life (Pryce and Brotherstone, 1999), it is worthwhile to estimate genetic and phenotypic correlations of ProdI with these functional traits.
The objectives of this study were 1) to investigate the formula used to estimate ProdI, using the real lactation yields, and to find possibilities of improving it; 2) to estimate the variance components of ProdI; and 3) to quantify genetic and phenotypic correlations between ProdI and udder health, longevity, fertility, and conformation.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Thirteen years of production records with calving dates from 1991 through 2004 were in the production data set. This data set consisted of all lactations with at least one test-day, resulting in 221,850 cows with 551,784 lactations. The data set contained the cow identification number (ID); parity; herd ID; date of calving; milk, fat, and protein yields in 305-d and in total lactation; DIM; and status of the lactation (i.e., completed or not). Records of SCC were given as test-day results, recorded since 1994. In total 1,172,244 insemination records were available from 1994 to 2004. Conformation data, recorded since 1991 in Switzerland, included 21 linear traits, 5 composites, and final score for 140,325 cows. Linear traits were measured on a scale of 1 to 9; for first lactation heifers, composites were on a scale from 50 to 87; for dairy character, the scale went up to 90.
The pedigree file contained ID of the animal and of its sire and dam. Number of bulls and cows in the pedigree file were 153,011 and 589,225, respectively. Pedigrees were traced as far back as 1940 for the bulls and 1942 for the cows.
Trait Definitions and Data Editing
A data set containing ProdI, SCS, days to first service (DFS), NRR, longevity, and conformation was constructed with the available data.
ProdI was calculated as follows:
![]() | ([1]) |
where Mx is the milk yield or energy-corrected milk (ECM) yield in lactation x, x = 1, 2 or 3. The ProdI for milk yield was defined as MilkI and, for ECM yield, was defined as ECMI. This formula is adapted from Schleppi and Bigler (2002). For milk yield, standard 305-d lactation yields were used, which consisted of cows with at least 270 DIM with lactation recorded as "ended." The standard 305-d milk yields meant that lactations with 270 to 304 DIM were multiplied by a correction factor used by the Holstein Association of Switzerland (Wegmann, 2004, personal communication).
The ECM was computed with the following equation:
![]() | ([2]) |
where c is the correction factor for lactations with <305 d, M is the milk yield, Fd is the fat deviation, and Pd is the protein deviation (Kirchgessner, 1997). The following restrictions were applied to exclude outliers: milk and ECM yield had to be between 2000 and 16,000 kg, age at first calving had to be at least 500 d and < 2500 d for third calving.
In a later analysis, the ProdI formula was reconstructed to try to improve its accuracy in describing the increase of production. The difference between the first and second lactation, respectively, was defined as MilkI12 and ECMI12; the difference between the first and third lactation was defined as MilkI13 and ECMI13, respectively; and, finally, the difference between the second and third lactation, respectively, was defined as MilkI23 and ECMI23. These ProdI formulas aim to assess the increase of production in a more detailed way.
Lactation SCS was computed using test-day SCC data and was used for the variance component estimates. Test-day SCS was calculated as
![]() | ([3]) |
Lactation SCS was computed as the average of test-day SCS within lactations (i.e., made up to the 305th d of the lactation). The use of lactation SCS instead of SCC is based on the linearity of SCS (Ali and Shook, 1980) and has been used to analyze test-day SCC in many studies (Mrode and Swanson, 1996; de Haas et al., 2002). For this data set, plausibility restrictions were set at 0 and 8 for lactation SCS.
The first trait used for fertility was DFS and was computed as the interval between calving and first recorded insemination (in days). Plausibility values were set at 20 d as the lower limit and 180 d as the upper limit to avoid possible biases from cows that were physiologically abnormal or incorrect recording of calving or insemination dates. Days to first service after each of the first 3 calvings were computed. The second fertility trait used, NRR, is defined in Switzerland as a binary trait recording the absence or presence of insemination during the 56 d following the first insemination. The EBV of the sires were used for the evaluations with NRR to avoid the problem related to the analysis of binary traits. The NRR EBV are estimated by the Holstein Association with an animal model (Schnyder and Stricker, 2002).
The EBV of sires for their daughters longevity were used in calculating relationship between ProdI and longevity. In Switzerland, longevity is defined as the productive life, which is the life span between first calving and the last test date on official milk recording. Longevity EBV are estimated with a sire-maternal grandsire survival model (Vukasinovic et al., 2001).
Conformation analyses were made with classification results of first lactation heifers. Age at first calving had to be <1250 d, and the classification had to be made within 1 yr after calving.
The data were edited in a stepwise manner. First, only lactations 1 to 3 of the production data set were kept. Only cows with records in all 3 lactations and with known sires were taken for further analysis. Records on other traits (SCS, DFS, and conformation) were added when available (i.e., not all cows had records on all other traits). Final editing was done by excluding daughters of bulls with <6 daughters and cows in herds with <6 cows. This resulted in a data set of 42,807 cows with an observation for ProdI. Number of cows with records for other traits is shown in Table 1
. Calving dates of these cows were between January 1991 and July 2003.
|
Statistical Models and Analyses
Variance components were estimated with sire models. Considering the number of traits for which the correlation had to be computed, sire model was preferred over animal model for the genetic analyses.
The univariate sire model used was as follows:
![]() | ([4]) |
where y is the observation, µ is the overall mean, s is the random effect of the sire, and e is the residual. Fixed effects differed from trait to trait. The use of each one of them is summarized in Table 2
. The program ASREML (Gilmour et al., 2001) was used for estimation of all variance components and regression coefficients. This software estimates variance components by REML.
|
Because ProdI has information from all 3 lactations, all fixed effects related to production were used. The other traits had only those fixed effects related to the corresponding lactation. Moreover, classification effects (classifier and time to classification) were included as fixed effects for analyses of conformation traits. Bivariate sire model analyses were conducted to estimate the (co)variances and correlations between ProdI and the other traits. Genetic and phenotypic parameters were computed based on the estimated variance components.
Regressions (b) of daughter ProdI records on their sire EBV for longevity and NRR were estimated by a univariate analysis. The model was
![]() | ([5]) |
where EBVx is a sire EBV of trait x (longevity or NRR), fitted as a covariate, and all other terms are as described previously under equation [4].
| RESULTS |
|---|
|
|
|---|
|
|
p) ranged from 0.73 to 1.67. The heritability estimates were between 0.10 (heel depth) and 0.49 (stature). Composites had heritabilities between 0.14 (feet and legs) and 0.41 (frame/capacity). Final score had a mean of 78.35 and a heritability of 0.36. Standard errors for all functional traits were <0.05.
Correlation of ProdI with Functional Traits
Additive genetic and phenotypic correlations between MilkI and SCS, DFS, and conformation traits are presented in Table 4
. Each of the 3 lactations shows favorable additive genetic correlations between SCS and MilkI (0.30, 0.30, and 0.29). Some rump traits [loin (0.19) and pin width (0.20)] have significant negative genetic correlations with MilkI. Chest width (0.30) and heart girth (0.24) have moderate positive genetic correlations with MilkI. Phenotypic correlations were generally small but in the same direction as genetic correlations and had high standard errors.
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
ProdI Definitions
Evaluation of the standard ProdI formula was made with 3 changes on the formula of Schleppi and Bigler (2002). These were 1) the use of actual production yields instead of EBV to avoid the correlations used in the multivariate estimations of test-day EBV, 2) the absence of correction factors used to account for the differences in standard deviations in each lactation (because it was not used by the Holstein Association of Switzerland) (Wegmann, 2004, personal communication), and 3) the replacement of fat and protein yields by milk (MilkI) or ECM (ECMI) yields. The latter change was made to evaluate maturity rate as a functional trait: milk is a measure of the amount of fluid put in the udder of the heifer and, therefore, describes the constraints put on the udder during its early development. Energy-corrected milk is a measure of the energy content of milk and is, therefore, an indicator of the energy needed for milk production that cannot be used for growth; for the same reason, it is also a measure of the metabolic stress of the cow. Fat and protein production, which is currently used in Swiss breeding schemes, would have described these effects less accurately.
Later on, a second group of ProdI formulas was made to define more precisely the increase between each of the lactations. The formula of Schleppi and Bigler (2002) showed the global increase of the production from the first to the third lactations, but the increase between any 2 of them was not described. The standard ProdI does not reveal whether the increase was made at the second lactation, at the third lactation, or at both. The second group of formulas allowed the evaluation of each increase, i.e., between first and second lactation, first and third lactation, and second and third lactation. The similarities of variance components between the 3 ProdI in which first lactation was the base of comparison (MilkI12, MilkI13, and standard MilkI) show the importance of this lactation for maturity rate. Increases from first lactation to second, third, or both have the same patterns. For maturity rate, the emphasis must be put on what was produced in first lactation compared with any of the subsequent lactations. The MilkI23 and ECMI23 seemed to be different in both means and variances to standard MilkI and ECMI. The large coefficient of variation found for MilkI23 and ECMI23 shows that the increase in production for this trait has a high error, and, therefore, the increase is not as important as in ProdI, including first lactation. Moreover, MilkI23 and ECMI23 had heritabilities that were not different from zero, which resulted in high standard errors for all correlations estimated between these traits and other functional traits. For the latter reason, MilkI23 and ECMI23 will not be discussed any further.
Heritability estimates of the other ProdI were in the range of functional traits, not unlike udder health, fertility, and resistance to disease (Uribe et al., 1995; Kadarmideen et al., 2000; VanRaden et al., 2004). All ECMI showed higher means than MilkI. Therefore, a limiting factor for production in first lactation is energy needs (ECM) rather than udder capacity (milk). A cow with enough energy resources will produce the amount of milk related to it. Udder capacity will only play a minor role for the limitation of milk production. Therefore, it is good to breed for tight udders (i.e., high udder depth), as milk production will not be hindered much. Moreover, tight udders are positively correlated with longevity (Pasman and Reinhardt, 1999).
Correlations between MilkI and first lactation yield were computed to see whether there was a bias in the data set caused by the use of cows with at least 3 lactations ended. Cows culled in the first lactation because of poor milk production or on the contrary because of problems caused by a high first lactation milk yield could have biased the results. There is a slightly negative genetic correlation (0.23), but the standard error is relatively high (0.08). Therefore, we may assume that the bias caused by daughters culled during their first lactation is not big enough to have an impact on the analyses.
Correlations Between ProdI and Some Functional Traits
Genetic and phenotypic correlations were estimated between ProdI and SCS, DFS, NRR, longevity, and conformation. All ProdI have similar genetic correlation estimates with these traits. Therefore, to avoid redundancy, only the results of one ProdI (ECMI13) will be given, and when differences occur, they will be mentioned.
Udder health.
Somatic cell score was used as an indicator of udder health. Genetic correlations between lactation SCS and occurrence of clinical mastitis during the same lactation were medium to high (0.7) as reviewed by Mrode and Swanson (1996). Lund et al. (1994) reported a high estimate of 0.97 between SCS and clinical mastitis. The mean SCS of the current study (2.48) is in the same range as those estimated by de Haas et al. (2002). Weller et al. (1992) found a lower increase in SCS from first to third parity (2.00 to 2.28) than in the present study. Heritability of SCS (0.17 to 0.19) is high for this trait: Mrode and Swanson (1996) reported average heritability estimates of 0.11 for each of the first 3 lactations based on a review of studies. Weller et al. (1992) observed a heritability estimate of 0.19 for first lactation SCS, similar to the present study. The genetic correlation of ECMI with SCS in third parity is larger than with SCS in first parity. The amount of milk produced during the first lactation has an unfavorable effect on udder health in the same lactation (Kadarmideen et al., 2000), but, based on the correlations computed between ProdI and the different SCS, it has a stronger negative effect on SCS in later lactations. High milk production in the first lactation reduces the action of the immune system in the udder and leads to damage that will last for the subsequent lactations (Rupp et al., 2000). A moderate first lactation production, compared with adult production, is part of having a healthy udder; it allows enough protection of the udder and, therefore, allows time for it to develop into a full producing mammary system. Another understanding of the correlation between SCS and maturity rate would be that the damage to the udder by a case of clinical mastitis in first lactation hinders a high production increase in the following lactations.
Average additive genetic correlations between ProdI and SCS of the first 3 lactations (0.29) are slightly larger than those computed by Krogmeier et al. (2003), who worked with Fleckvieh (dual-purpose German Simmental) and Brown Swiss breeds (0.20 and 0.25, respectively). One of the reasons for the smaller values found by Krogmeier et al. (2003) could be a consequence of the breeds considered. Fleckvieh and Brown Swiss have a lower genetic potential for milk production and produce, on average, less milk in first lactation; therefore, the effects on udder health are less important.
Fertility traits.
The 2 traits used in the fertility analyses (DFS and NRR) were those recommended by Groen et al. (1998b) to describe fertility. Fertility was defined by Darwash et al. (1997) as "the ability of the animal to conceive and maintain pregnancy if served at the appropriate time in relation to ovulation." The first trait (DFS) is an indicator of observable cyclicity; the second one (NRR) indicates the probability of becoming pregnant at insemination.
Mean DFS was similar to the value reported by Kadarmideen et al. (2000). Heritability of DFS was in the same range as Pryce et al. (1998), who also computed a DFS for each of the first 3 lactations. However, it was lower than the heritability of 0.10 found by Schnyder and Stricker (2002) in Swiss Holsteins. Additive genetic correlations between ProdI and DFS in first lactation were negative, except with MilkI12. Correlations estimated with ECM had a larger negative value than those computed with milk. Once again, the importance of the energy needs is seen at the genetic level: cows with genetic potential for low first lactations, when compared with following lactations, start their first ovarian cycle of that lactation sooner than genetically high-producing cows. This result stresses the importance of having good maturity rate values for good fertility. Genetic correlations between ProdI and DFS in second and third lactations were highly variable and had higher standard errors.
Regression of ProdI on NRR was negative. This means that cows that are not re-inseminated in the 56 d following the first insemination have a lower maturity rate than cows that are; in other words, a good ability of becoming pregnant at first insemination is related to a slower maturity rate. This phenomenon may be explained by the positive genetic correlation between NRR and DFS in Swiss Holsteins (Schnyder and Stricker, 2002). This positive correlation means that cows bred soon after calving are more subject to not becoming pregnant at that particular insemination. However, Pryce et al. (1998) and Kadarmideen et al. (2000) found a negative correlation between these traits in UK dairy cows. It must be kept in mind that the NRR correlation in the present study is based on the EBV of the sires, giving high standard errors to the estimates because the EBV are based on binary observations. Moreover, these EBV are computed with insemination records from all lactations and are not values for each of the first 3 lactations. Correlation estimated between ProdI and NRR should be analyzed based on actual cow NRR records of the lactation concerned to get more accurate results. Krogmeier et al. (2003) found no correlation between ProdI and 90-d NRR (0.05 and 0.03 for Fleckvieh and Brown Swiss, respectively) using EBV values for both traits.
Longevity.
Regression coefficients of ProdI on sire EBV for longevity were positive and significant. Cows whose sires transmit good longevity had a slightly higher maturity rate than those from sires transmitting short productive life. Moreover, traits correlated to longevity were also correlated to ProdI. One of these traits is SCS, which was reported to have a correlation of 0.27 to life span (i.e., number of lactations a cow completes) by Pryce and Brotherstone (1999). Another trait positively correlated to both ProdI and longevity is udder depth (Weigel et al., 1997; Pasman and Reinhardt, 1999; Vukasinovic et al., 2002). The structure of the data set reduced the regression of ProdI on longevity because cows culled before the end of the third lactation were not included. Further research should find a way to eliminate the bias in the data set to find a more accurate estimate of the regression or even estimate the correlation between ProdI and longevity.
Conformation traits.
Conformation indicates how far a cow "conforms" to the True Type model (Holstein Canada, 2004). This model describes the cow that has the ability to produce large amounts of milk and last in the herd. For this reason, conformation is used as an indicator of longevity. Heritability estimates computed in the present study for conformation traits were in the same range as earlier estimates on Swiss Holsteins (Kadarmideen and Wegmann, 2003) except for stature and heart girth, which were defined differently in the current study (measures in scores instead of centimeters).
Of all conformation traits, chest width has the highest positive correlation to ProdI. This trait is also negatively correlated to milk yield in first lactation (Short and Lawlor, 1992). As first lactation milk yield has a negative weight in the ProdI formula, these cows might often have high ProdI. Moreover, first lactation cows showing width through the front end have a better energy balance throughout lactation and, therefore, less metabolic stress (Veerkamp, 1998). The latter is a hindrance to the increase in production in later lactations.
Of all conformation traits, angularity had the strongest negative correlation to ProdI. Angularity is the trait that often characterizes "milky" animals, as is corroborated by the high correlation between angularity and milk yield in first lactation (0.54; Short and Lawlor, 1992). Conformation of high-producing animals is de-fined by having more angularity than animals with low production. However, cows showing lots of angularity, because they produce more milk in first lactation, have difficulties in keeping enough energy to grow and, therefore, to be prepared for the production of more milk in the following lactations. Hence, production increases less than for cows showing less angularity.
Selection on ProdI
It would be interesting to know how to include maturity rate in breeding goals. The trait, ProdI, has some potentially useful correlations for selection. Moreover, there is variance in the PTA of bulls, which enables selection. Correlations to other traits make ProdI a useful functional trait. However, ProdI should not be selected as a single trait. This would lead to cows with low first lactation yield, therefore less profitable cows. Selection on actual production yield should be supplemented by selection on ProdI. The main goal of further research should be the estimation of a precise correlation between ProdI and longevity. This correlation and those already found in this work would enable estimation of the economic value of this trait. Time needed to compute ProdI is also a restriction to the use of that trait. With the current formulas, 2, or eventually 3, lactations have to be finished before a value can be calculated. Selection choices have already been made at that time.
It is important to note that the trait ProdI is defined as the ability of the cow to increase its production gradually from early to later lactations, regardless of actual amount of milk produced within lactations. Persistency, another functional trait, is defined as the ability of cows to sustain a constant level of production at all stages within lactation, whereas ProdI is related to an across-parity gradual increase in yield. Persistency within lactation is a trait that was not considered in the present study. This indication of the "flatness" of the lactation curve is important to select cows with fewer problems at the start of lactation because of a high metabolic stress. Krogmeier et al. (2003) reported a high positive correlation between ProdI and lactation persistency (0.61 for Brown Swiss).
It seems that cows with balanced milk yield throughout lactation show an increase in milk yield over lactations. Flat lactation curves and increases after each lactation should be the ultimate breeding goal from a metabolic stress point of view.
| CONCLUSIONS |
|---|
|
|
|---|
Received for publication August 26, 2004. Accepted for publication November 11, 2004.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |