|
|
||||||||
1 Canadian Dairy Network, 150 Research Lane, Guelph, ON, Canada N1G 4T2
2 Centre for the Genetic Improvement of Livestock, Department of Animal and Poultry Science, University of Guelph, ON, Canada N1G 2W1
Corresponding author: B. L. Muir; e-mail: muir{at}cdn.ca.
| ABSTRACT |
|---|
|
|
|---|
Key Words: genetic parameter persistency reproduction
Abbreviation key: AFS = age at first insemination as a heifer (in days), CD = calving difficulty, CI = calving interval, MILK = estimated 305-d milk yield in first lactation, MTP = multiple-trait prediction, NRRC = nonreturn rate at 56 d after first insemination in first lactation, NRRH = nonreturn rate at 56 d after first insemination as a heifer, PeakD = day in milk of peak milk yield
| INTRODUCTION |
|---|
|
|
|---|
Increasingly, research indicates an undesirable relationship between milk yield and fertility (Pryce et al., 1997; Dematawewa and Berger, 1998; Abdallah and McDaniel, 2000). How total yield is spread throughout the lactation may influence this antagonistic relationship between yield and fertility. Persistency can be defined as the ability to maintain the level of milk production after peak milk yield. A possible way to increase total yields without increasing the occurrence of disease or reproductive failure is to select for increased lactation persistency in addition to total production. This selection strategy could be attempted by decreasing stressful peak yields and maintaining a high level of production after peak yield, thereby flattening and extending the lactation curve. The Wilmink function (Wilmink, 1987) has been used in the past to model lactation curves and define persistency (Jamrozik et al., 1997; Schaeffer et al., 2000; van der Linde et al., 2000).
In general, producers breed cows as early as possible in lactation, believing that this strategy will enhance profitability. This management strategy forces higher producing cows to be inseminated when they are producing 45 to 50 kg daily (Arbel et al., 2001). Poor conception rates and reproductive inefficiency may contribute to the increasingly antagonistic relationship between yield and fertility that has been found (Pryce and Veerkamp, 2001).
Estimates of genetic parameters for persistency (Jamrozik et al., 1998; Rekaya et al., 2000; Van der Linde et al., 2000; Jakobsen et al., 2002) and fertility (Hoekstra et al., 1994; Pryce et al., 1997; Kadarmideen et al., 2000; Veerkamp et al., 2001) are abundant in the literature, as are relationships between fertility and milk production (Van Arendonk et al., 1989; Dematawewa and Berger, 1998; Veerkamp et al., 2001; Haile-Mariam et al., 2003). Investigation into relationships between persistency and reproductive performance is lacking. Therefore, the objective of this study was to estimate genetic relationships between lactation persistency and reproductive performance in first-lactation Holsteins. The timing and height of the lactation curve at peak yield have been suggested as alternative measures to persistency to genetically manipulate the shape of the lactation curve (Ferris et al., 1985). The genetic correlation between peak milk yield and 305-d milk yield was estimated to be high, and both traits had similar estimated genetic correlations with persistency and day of peak milk yield (Muir, 2004). In addition, 305-d milk yield is a more widely used production trait. Therefore, day in milk of peak milk yield (PeakD) and estimated 305-d milk yield in first lacation (MILK) were selected as additional traits for analysis to relate with persistency and reproductive performance. The hypothesis was that cows with persistent lactations have better reproductive performance.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Lactation data obtained from the Canadian Dairy Network were first-parity test day yields from lactations made between 1990 and 2003 by Canadian Holsteins. Days in milk ranged from 5 to 305 d. Daily milk yield was required, but records with missing fat and protein yield were allowed. Lactations having the date of first test >50 DIM and/or average interval between tests >50 d were excluded. Daily milk yields >90 kg or 1.5 kg were also excluded to avoid possible erroneous data.
Fertility and production data were merged, and a sample of data was created by randomly selecting herds and then selecting all records within that herd for genetic parameter estimation. In doing so, approximately 50% of herds on milk recording were selected. To be included in the analysis, each heifer was required to have all reproductive performance traits of interest (AFS, NRRH, CD, NRRC, and CI) and, therefore, must survive until second calving. Routine collection of complete reproductive information has only recently been adopted in Canada; therefore, the oldest animals in the data set were born in 1995. The final data set contained all 8 traits of interest on 33,312 first-lactation Holsteins.
Models
A multiple-trait prediction (MTP) procedure was used to estimate 305-d lactation yields from the individual milk, fat, and protein test day yields (Schaeffer and Jamrozik, 1996). The MTP method is based on Wilminks function (Wilmink, 1987) and incorporates information about standard lactation curves and covariances between yields for milk, fat, and protein. Wilminks model for one trait is
![]() |
where y is yield on day t of lactation, a is a regression coefficient related to the maximum daily level of production, b is a regression coefficient related to production decrease after peak yield (persistency), c is a regression coefficient related to production increase toward peak, and the factor 0.05 is related to the approximate DIM when peak milk yield occurs. Test day yields were weighted by their relative variances, and standard lactation curves of cows from similar breed, region, lactation, age at calving, and season of calving were used for the estimation of lactation curve parameters for each cow. For this study, only first lactations were considered. Using the 3 estimated Wilmink curve parameters for each animal, persistency, MILK, and PeakD for first lactation were also estimated. The definition of persistency selected for use in this study was the Wilmink b parameter, which was the estimated slope of the lactation curve after peak milk yield. This trait was chosen as a phenotypic measure of persistency because of its ease of calculation and moderate heritability (Muir, 2004) and its ability to be related with reproductive performance and other economically important traits on a large number of animals. In addition, this measure was easily computed with a software package currently used by milk-recording agencies to provide on-farm management information. The MILK was calculated by summing all estimated individual test day yields from d 5 to 305 DIM.
An 8-trait animal model was applied to calculate genetic parameters. Three production traits (persistency, PeakD, and MILK) and 5 reproductive performance traits (AFS, NRRH, CD, NRRC, and CI) constituted the 8-trait analysis. In matrix notation, the model for the 8 traits can be written as
![]() |
where b includes all fixed effects, a is a vector of additive genetic animal effects, p is a vector of herd within year-season of birth effects, q is a vector of service sire or sire of calf effects, e is the vector of residual effects, and X, Z, W, and M are the incidence matrices. Assume that
![]() |
and
![]() |
with
![]() |
where G is a genetic covariance matrix (8 x 8) among traits, assumed to be the same for all heifers; P is a herd within year-season of birth covariance matrix (8 x 8) among traits, assumed to be the same for all heifers; A is the additive genetic relationship matrix among the animals;
is Kronecker product function (Searle, 1982); R is a residual covariance matrix (8 x 8) between traits; Q1, Q2, and Q3 are matrices of zeros, with a 1 in the diagonal element corresponding to NRRH, CD, and NRRC, respectively; and
q12,
q22, and
q32 are the service sire effect for NRRH, sire of calf effect for CD, and service sire effect for NRRC variances, respectively.
Fixed and random effects included in the genetic models for each trait are shown in Table 1
. Year of birth ranged from 1995 to 2000. Four seasons of calving were defined, December to February, March to June, July to September, and October to November. Age at first calving ranged from 19 to 36 mo. Random effect of herd within year-season of birth had 21,265 levels (5725 herds) and was included in all models to account for contemporary group effect properly. Year and season of birth was included as a fixed effect in every model to account for time trends. Random effect of service sire was included in the models for NRRH and NRRC, and these effects were assumed to be uncorrelated. Effect of sire of calf was included in the model for CD. Service sires and sires of calves were all assumed to be unrelated. Residuals were assumed to be uncorrelated. Although NRRH and NRRC were recorded as binary traits, and CD as a categorical trait, records for these traits were assumed to be normally distributed. It was assumed that there was no interaction between herd within year-season of birth and animal genetic effects. Covariances between effect of service sire and effects of sire of calf were assumed to be zero. All environmental factors not included the model were assumed to be the same for all records or contributed to the residual.
|
Flat priors were assumed for all fixed effects. Inverted Wishart distributions were assumed for genetic, herd within year-season of birth, and residual covariances. Starting values for all variances and covariances were found by compiling prior knowledge of the variances of all traits and assuming a correlation of 0.05 between all traits to derive covariances. The analysis was performed on 33,312 animals with records with a total of 96,218 pedigree animals.
Heritability was calculated as additive genetic variance (
2a) divided by phenotypic variance (
2P), where
2P =
2a +
2YS:H +
2e for AFS, CI, persistency, PeakD, and Milk, and
2P =
2a +
2YS:H +
2SS1,
2SC, or
2SS2 +
2e for NRRH, NRRC, and CD, respectively, and where YS:H, SS1, SC, and SS2 were random effects of herd within year-season of birth, effect of service sire on NRRF, effect of sire of calf on CD, and effect of service sire on NRRC, respectively. Additive genetic correlations were calculated as
i,j ÷ (
i x
j), where
i,j was the covariance between trait i and trait j, respectively, and
i and
j were the standard deviations of trait i and trait j, respectively.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
Genetic Correlations Between Reproductive Performance Traits
Estimated genetic correlations between reproductive performance traits were variable and lower (range = 0.50 to 0.78) than other recently reported estimates (Pryce et al., 1997; Dematawewa and Berger, 1998; Kadarmideen et al., 2003). Estimated genetic correlation between the 2 heifer reproductive performance traits (AFS and NRRH) was small and positive (0.08 ± 0.12) with a large standard error, indicating that the age of the animal when it is first inseminated had little or nothing to do with conception success at that insemination. Estimated genetic correlation of AFS with NRRC was negative and stronger than that with NRRH (0.20 ± 0.08), indicating that heifers younger than average when first inseminated (that came into estrus for the first time earlier in life) tended to have a better chance of conceiving at first insemination after first calving. Estimated genetic correlation between NRRH and NRRC was only moderately positive, indicating that conception success as a virgin heifer could not accurately predetermine conception success as a first lactation cow and that different genetic and physiological determination of fertility between parities may exist. Batra et al. (1987b) showed that genetic correlations between heifer fertility and cow fertility were close to zero, indicating that heifer and cow fertility seem to be genetically independent. Nonreturn rate as a heifer and as a first lactation cow (NRRH and NRRC) were genetically unrelated to CI.
Moderate relationships existed between CD and other reproductive traits. Heifers younger than average when first inseminated tended to have greater dystocia at first calving (rg = 0.35 ± 0.06). First-lactation cows that had greater CD also tended to conceive more successfully at first insemination after calving (rg = 0.41 ± 0.08). A longer CI was genetically associated with a more difficult first calving. Perhaps breeding of difficult calving heifers was delayed until after stressful high peak milk yields, when a successful insemination would be more probable and lactations were subsequently prolonged longer than average. Kadarmideen et al. (2000) reported that genetic correlations between dystocia and other reproductive traits were close to zero.
Genetic Correlations Between Production Traits
Estimated genetic correlation between PeakD and persistency was 0.54 (± 0.07), indicating that as the interval from initiation of lactation to peak yield increased and persistency improved. Selection on one trait should improve the other, but only moderately so. Other researchers have found the genetic correlation between persistency and timing of peak yield to be positive and high (Batra et al., 1987b; Rekaya et al., 2000) and some to be positive and moderate (Ferris et al., 1985). The PeakD was also positively genetically correlated with MILK (0.63 ± 0.06), indicating that as total lactation yield increased, so did the interval from initiation of lactation to peak yield.
Estimated genetic correlation between MILK and persistency in first lactation was 0.21 (± 0.06) and was higher than expected. Several researchers have found correlations (positive and negative) of persistency with 305-d yield ranging from 0 to over 0.60, varying with the way persistency was defined (Ferris et al. 1985; Jamrozik et al., 1998; Rekaya et al., 2000; van der Linde et al., 2000; Jakobsen et al., 2002). Gengler (1996) recommended the use of a measure of persistency that is genetically independent of 305-d milk yield; however, in this study, a small positive genetic correlation was estimated between yield and persistency, indicating that selection for increased MILK would slightly improve persistency. Clearly, relationships exist between the shape of the lactation curve and total production. Minimizing this relationship will allow more efficient selection for total lactation yield and persistency simultaneously.
Genetic Correlations Between Reproduction and Production
Estimated genetic correlations between reproduction traits and persistency ranged from 0.17 to 0.43 (Table 3
). Age at first insemination was negatively correlated with persistency (0.17 ± 0.07), indicating that heifers that were younger when first inseminated tended to have better persistency in first lactation. The correlation between persistency in first lactation and NRRH was positive (0.16 ± 0.13) but had a large standard error, so caution should be used when making statements about this possible relationship. Reproductive performance as a virgin heifer could indicate future shapes of lactation curves, and these relationships should be further investigated. Lean et al. (1989) reported that cows with higher than average peak milk yields and poorer than average persistency were less likely to conceive in 1 or 2 services. Bar-Anan et al. (1985) defined persistency as average daily yield divided by peak yield and found that cows with higher lactation persistency had better reproductive performance. Haile-Mariam et al. (2003) reported the absence of genetic antagonism between fertility and persistency. Nonreturn rate at 56 d after first insemination in first lactation was also positively correlated with persistency (0.32 ± 0.12), indicating that first lactation cows that tended to conceive on first insemination also tended to have better persistency.
The genetic correlation between CI and persistency was estimated to be 0.17 (± 0.09). As expected, the greater persistency was in first lactation, the longer the interval was from first to second calving. Longer CI are generally viewed as undesirable; therefore, genetically an antagonistic relationship existed between persistency and CI in first lactation. However, Arbel et al. (2001) reported that there was an economic advantage to extending lactations by 60 d in high yielding dairy heifers in first lactation. Stott et al. (1999) reported that the marginal value for CI at the optimum estrous detection rate was $11.75/d rising to $14.06/d if breeding was delayed.
The estimated genetic correlation between CD and persistency was strong (0.43 ± 0.08). Heifers that had a difficult first calving tended to have more persistent first lactations. Heifers that had difficulty calving would also tend to have lower peak yields, possibly causing the antagonistic relationship between CD and persistency. Previously, speculation was made that heifers with greater CD may have longer intervals from calving to first insemination, and, therefore, fertility would be better later in lactation after peak milk yield. If this were true, one would expect the correlation between CI and persistency to be > 0.17. To test this hypothesis, DIM at insemination could be included as a covariate in the model for CI and correlated with persistency in a bivariate animal model or number of days from calving to first insemination could be included as a correlated trait in a multiple trait animal model. Alternatively, one could speculate that persistency is genetically correlated with less severe early lactation negative energy balance. Higher persistency appears to be positively correlated with later DIM at peak milk yield and lower peak yield (Ferris et al., 1985), both of which are indicators of less severe negative energy balance. Veerkamp et al. (2000) reported that the interval from calving to first service decreased when negative energy balance was less severe.
Estimated genetic correlations between reproductive performance and PeakD ranged from 0.11 to 0.31 (Table 3
). The strongest genetic correlation was between CI and PeakD and suggested that the longer into lactation an animal experienced peak daily milk yield, the longer the interval was between first and second calving. A negative genetic relationship between AFS and PeakD existed and indicated that heifers that were younger when first inseminated tended to have a later peak yield in first lactation. Heifers that were younger when first inseminated also tended to have better persistency. Later peak yields were genetically associated with better persistency, so these two estimated genetic correlations were expected. Cows that had difficulty calving tended to have a later PeakD (rg = 0.15 ± 0.08).
Estimated genetic correlations between reproductive performance and MILK were similar in sign, but different in magnitude, to those with persistency. The strongest genetic correlation existed between MILK and CI and indicated that a longer CI was associated with greater MILK in first lactation (0.51 ± 0.11). Similar genetic correlations have been reported between CI and MILK in the literature, ranging from 0.23 to 0.67 (Hoekstra et al., 1994; Veerkamp et al., 2001; Haile-Mariam et al., 2003; Kadarmideen et al., 2003).
The estimated genetic correlation between AFS and MILK was 0.37 (± 0.05) and indicated that heifers that were younger when first inseminated tended to have greater MILK in first lactation. Heifers that were younger when first inseminated probably came into heat for the first time earlier in life and therefore showed a greater propensity to mature. Genetic ability to mature faster and reach a greater proportion of their mature body weight at an earlier age may have been the biological driving force causing AFS to be genetically correlated to greater total yield expected in first lactation. Early maturing heifers with a lower AFS also had better persistency and later peak yields in first lactation. Cow fertility (NRRC) was uncorrelated with MILK. Other estimates of genetic correlations of MILK with conception success in the literature have been much stronger (Hoekstra et al., 1994; Veerkamp et al., 2001). Correlations between persistency and conception success were stronger than those found with MILK. Therefore, conception success as a first-lactation cow tended to be more related to the dispersion of yield in first lactation rather than total yield.
Age at first insemination was moderately heritable and negatively genetically correlated with persistency and with MILK. Therefore, selecting heifers that come into heat for the first time at a younger age could improve persistency and MILK in first lactation. Heifers inseminated at a younger age also tended to have better conception success (NRRC) in first lactation. However, an unfavorable correlation between AFS and CD existed. Selecting for earlier AFS may result in improved persistency, 305-d yield, and NRRC, but early maturing heifers may need to be mated to calving ease sires to decrease dystocia at first calving.
A positive estimated genetic correlation was shown between CI and persistency in first lactation. Possibly, first-lactation heifers that were open longer would seem to be more persistent because they would be directing less energy toward the unborn calf. Additionally, low persistency could have been detrimental to conception after first calving. However, the estimated genetic correlation between persistency and NRRC was 0.32 and indicated that conception success at first insemination after calving tended to be related to better persistency in first lactation. More persistent lactations were somewhat associated with higher 305-d yields. Perhaps insemination of first-lactation cows that experienced dystocia was delayed longer than average, possibly long after stressful peak yields, and, therefore, better conception success resulted. Having data available on the exact number and timing of all inseminations and conception rates for each insemination, in addition to having the reasons for and length of voluntary waiting periods, would help assign cause and effect to some of these relationships.
Several researchers have shown a consistently antagonistic relationship between reproductive performance and milk yield. In addition, studies show that genetically, persistency can vary with total yield; minimizing this relationship will allow for more effective selection of total yield and shape of lactation curves. Prior to conducting this study, research lacked investigation into how these relationships all fit together, i.e., the relationships between reproductive performance and persistency, under conditions of continuous selection for milk yield.
In this study, persistency, MILK, and PeakD were calculated by fitting lactation curves to phenotypic test day yields for individual cows (using MTP) and were not directly observed, which might have resulted in less accurate estimates of covariance components for these traits. Using a hierarchical model (Rekaya et al., 2000) or a random regression test day model to estimate lactation curve parameters directly (Jamrozik et al., 1997) would properly account for all covariances among observations (test days), which should improve the accuracy of the estimates of covariance components. Once genetic parameters of lactation curves are estimated, relationships among PeakD, persistency, and MILK could be derived analytically.
Some differences in estimated relationships between reproductive performance and persistency exist between these results and those that have been previously reported (Bar-Anan et al., 1985; Lean et al., 1989; Haile-Mariam et al., 2003). In this study, heifers were required to have all 5 reproductive performance traits recorded (and therefore survive to a second calving). Selection of data in this way might have been biased toward healthy, more fertile heifers, as the poorest heifers (in terms of health and fertility) were probably culled before being able to calve for a second time. Although means for the selected data set in this study were similar to those of the entire Canadian population (Table 2
), selection bias might have affected genetic parameter estimates. The magnitude of this bias was not quantified but was expected to be minimal. Estimates of heritability of reproductive traits found in this multiple-trait analysis using a selected population were nearly identical to single trait estimates found using an unselected data set and similar models (Fatehi, personal communication, 2003). Selection bias could be avoided by creating a data set allowing for the possibility of missing observations for fertility traits. Computational packages that allow for missing observations and complicated covariance structures on a large number of animals should be investigated for calculation of genetic parameters for several fertility traits simultaneously with production traits. Persistency was defined differently here, than in the previously mentioned studies and might have contributed to some differences found. Only first lactation was investigated in this study. Relationships between persistency and fertility in later lactations may be different than those found here in first lactation and should be investigated further.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication March 26, 2004. Accepted for publication June 3, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. Togashi and C. Y. Lin Genetic Modification of the Lactation Curve by Bending the Eigenvectors of the Additive Genetic Random Regression Coefficient Matrix J Dairy Sci, December 1, 2007; 90(12): 5753 - 5758. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. D. R. N. Appuhamy, B. G. Cassell, C. D. Dechow, and J. B. Cole Phenotypic Relationships of Common Health Disorders in Dairy Cows to Lactation Persistency Estimated from Daily Milk Weights J Dairy Sci, September 1, 2007; 90(9): 4424 - 4434. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. P. P. Macciotta, D. Vicario, and A. Cappio-Borlino Use of multivariate analysis to extract latent variables related to level of production and lactation persistency in dairy cattle. J Dairy Sci, August 1, 2006; 89(8): 3188 - 3194. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Harder, J. Bennewitz, D. Hinrichs, and E. Kalm Genetic parameters for health traits and their relationship to different persistency traits in German Holstein dairy cattle. J Dairy Sci, August 1, 2006; 89(8): 3202 - 3212. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. I. Weller, E. Ezra, and G. Leitner Genetic analysis of persistency in the israeli holstein population by the multitrait animal model. J Dairy Sci, July 1, 2006; 89(7): 2738 - 2746. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Jamrozik, J. Fatehi, G. J. Kistemaker, and L. R. Schaeffer Estimates of Genetic Parameters for Canadian Holstein Female Reproduction Traits J Dairy Sci, June 1, 2005; 88(6): 2199 - 2208. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Wall, I. M. S. White, M. P. Coffey, and S. Brotherstone The Relationship Between Fertility, Rump Angle, and Selected Type Information in Holstein-Friesian Cows J Dairy Sci, April 1, 2005; 88(4): 1521 - 1528. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |