|
|
||||||||


* Division of Animal Physiology, School of Biosciences, University of Nottingham, Sutton Bonington Campus, Leics, LE12 5RD, United Kingdom
Animal Biology Division, Scottish Agricultural College, Edinburgh, EH9 3JT, United Kingdom
Roslin Institute (Edinburgh), Roslin, Midlothian, EH25 9PS, United Kingdom
Corresponding author:
M. D. Royal; e-mail:
melissa.royal{at}nottingham.ac.uk.
| ABSTRACT |
|---|
|
|
|---|
Key Words: milk progesterone calving interval body condition score linear type traits
Abbreviation key: CI = calving interval, CLA = interval to commencement of luteal activity postpartum, HUKI = Holstein UK and Ireland, PCH = percentage Holstein, PIN95 = 1995 profit index, ra = additive genetic correlation
| INTRODUCTION |
|---|
|
|
|---|
The analyses of milk progesterone profiles have shown that a number of atypical ovarian patterns are phenotypically associated with reduced fertility in the UK dairy cow (Darwash et al., 1997a; Royal et al., 2000a). Three such patterns—interval to commencement of luteal activity (CLA) postpartum, persistent corpus luteum type I, and length of the first luteal phase postpartum—have moderate levels of heritability in todays Holstein Friesian population of dairy cattle (Darwash et al., 1997b; Royal, 1999; Veerkamp et al., 2000; Royal et al., 2000b, 2002). Genetic regressions of lnCLA on sires PTA for milk, fat, and protein yield and one of the UKs national dairy selection indices related to production, the 1995 profit index (PIN95), are highly unfavorable (Royal et al., 2000b, 2002). Furthermore, the additive genetic correlation (ra) between predicted peak milk yield (d 56) and lnCLA is large (0.36) and unfavorable (Royal et al., 2002). Other studies (Pryce et al., 2000) have reported a moderate and negative ra (–0.40) between a traditional measure of fertility, calving interval, (CI), and average BCS (an indicator of energy balance), in addition to moderate ra with a number of linear type traits (e.g., 0.33, 0.47 with stature and angularity, respectively), indicating that taller, more angular cows have longer CI.
To assess the usefulness of CLA in future breeding programs to improve fertility, it is necessary to determine ra with traits in the breeding goal or used in selection indices, such as CI, production, and linear type traits. In the UK, linear type traits and BCS are routinely recorded by Holstein UK and Ireland (HUKI). For herds that are part of the linear type classification scheme, it is compulsory to have heifers classified to provide information for sire type PTA estimation (Brotherstone and Hill, 1991). Trained classifiers record type traits on a linear scale from 1 to 9 according to biological extremes (e.g., for stature 1 is very short and 9 is very tall). BCS is used to assess the amount of body fat. The assessment is made visually and/or by touch, on the tail head and/or the loin area. It is used by many farmers as a management aid in deducing how well the current feeding needs are being met relative to stage of lactation in dairy cows. BCS and angularity are genetically correlated with CI (Pryce et al., 2000) and may also be good indicators of CLA. In addition, ra between BCS and angularity is high. Where ra between type and breeding goal traits are high, type information can be used to provide additional information. For example, other type traits (udder depth, foot angle, fore udder attachment, and teat length) are currently used to predict lifespan PTA (lifespan is the measure of longevity used in the UK) where actual lifespan information is not available or is sparse (Brotherstone et al., 1997). Therefore, type information may also be valuable in predicting fertility.
Little information is available on the genetic relationships between endocrine and traditional reproductive traits, and linear type traits commonly used in selection programs. As an initial step towards filling this gap in knowledge, a database of over 1200 records on cows with CLA and pedigree information was related via common sires to estimated PTA for CI, production, linear type traits, and average BCS. The objectives of the current analyses were to estimate genetic regressions of lnCLA on the described sire PTA values and to obtain some preliminary evidence on the sign and magnitude of genetic correlations.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
As there were relatively few sires in common between the milk progesterone and the sire PTA database, no further restrictions were made. Sires represented in the final data were generally bulls that were used extensively, and thus the accuracies were >0.7. A summary of statistics for the sire PTA used in these analyses are presented in Table 2
.
|
Pedigree and performance records.
Lactation and reproductive performance data from the milk progesterone database, in addition to information from two commercial databases (National Milk Records Plc, Chippenham, UK and HUKI Ricksmanworth, UK) were combined. Information from the milk progesterone database, related to 1212 lactations, was collected between October 1996 and March 1999.
Three-generation pedigrees were obtained for all cows in the study (i.e., up to and including great-grandparents). Wherever possible this was extracted from the HUKI database. Milk records were obtained from the NMR database for all cows, except those of Roslin Institute, which were directly available.
The Holstein percentages (PCH) of all cows in the study were either obtained from the HUKI database or calculated from the known pedigree using sire information and the origin of maternal ancestors. The distribution of PCH (for the cows monitored and their sires) in the database is illustrated in Royal et al. (2002b). Approximately 85% of the sires and 16% of the cows monitored in the current database were 100% North American Holstein. Approximately 7% of the sires and 5% of the dams were 100% British Friesian. The average percentage of North American Holsteins in the current database was 70.5%.
Data Integration
The milk progesterone database formed 169 paternal half-sib groups (group sizes between one and 40 daughters). The distribution of paternal half-sib family sizes in the database is illustrated in Royal et al. (2002b). A total of 923 maternal half-sib groups were present (including singleton groups), with the largest maternal half-sib group size of three. It was not possible to identify two of the 169 sires and 31 of the dams. Those unidentified were assumed to be unique.
Restricting analyses to those with sire PTA for production reduced the number of lactations to 1207, restricting analyses to those with sire PTA for type information reduced the number to 1094, and restricting analyses to those with sire PTA for CI and average BCS reduced the number to 1023.
Statistical Analysis
CLA has a long-tailed distribution (see Figure 1
); therefore to improve the properties of the distribution, it is common to log-transform the data. Thus, CLA was analyzed on the natural log-transformed scale as in previous publications (Darwash et al., 1997b; Royal et al., 2002).
|
|
| RESULTS |
|---|
|
|
|---|
Assumed estimates for heritability (h2) and phenotypic and genetic standard deviations (
P and
a, respectively) are presented in Table 3
. Genetic regression coefficients (b) of lnCLA on sire PTA, their standard errors, and levels of significance in addition to estimated ra are presented in Table 4
.
|
|
The regression of lnCLA on milk, fat, and protein yield were all positive (appropriate to ra range of 0.33 to 0.69). The only coefficient that was significantly different from zero was the regression of lnCLA on sire PTA for fat yield (P < 0.005). The magnitude was such that for every 10 kg increase in fat yield, CLA increased by 6.01% (1.6 d). The regression on sire PTA for milk yield approached significance (P < 0.1).
| DISCUSSION |
|---|
|
|
|---|
The genetic relationship between fertility and BCS has been investigated from several different angles over the past few years. Pryce et al. (2000) investigated the relationship between CI and average BCS, BCS in mo 1 to 10 of lactation and BCS change between mo 1 and 2, 1 and 3, and 1 and 4 of lactation. Veerkamp et al. (2000) investigated the genetic relationship of CLA to live weight and energy balance and, in a more recent analysis (Veerkamp et al., 2001), BCS to CI and interval to first service. Finally, the current analysis estimated the genetic correlation between lnCLA and average BCS. In previously published work described above (Pryce et al., 2000; Veerkamp et al., 2000, 2001), the genetic correlations estimated were negative (ranging from –0.06 to –0.88). Although it has often been suggested that these relationships reflect a yield effect upon fertility, attempts were made to account for yield by both Veerkamp et al. (2000) and Pryce et al. (2000), albeit by phenotypic adjustment. In both cases, correlations remained negative and of substantial magnitude, supporting the hypothesis that body tissue mobilization and partitioning for the purpose of production may occur as a result of some of the same processes that affect reproduction. The genetic correlations between BCS (energy balance), milk yield, and fertility are likely to reflect pleiotrophic gene action (i.e., a gene affecting two or more characters [Falconer and Mackay, 1996]). Thus, the genetic correlation could lie through either 1) hormones such as insulin, growth hormone, and insulin-like growth factors controlling intermediary metabolism (metabolic pathways by which the basic molecular building blocks in a cell are interconverted and incorporated into larger molecules) having direct effects on ovarian function or 2) reproductive hormones regulating ovarian function having direct effects on intermediary metabolism.
The physiological basis for these relationships may be as follows. Negative energy balance, responsible for the loss of BCS as occurs during peak lactation, is associated with changes in the secretion and circulating levels of a wide range of hormones. In particular, hormones controlling intermediary metabolism lead to the onset of gluconeogenesis and the mobilization of energy stored in the form of carbohydrates, fat, and protein. Selective breeding for high milk yield affects the same hormone systems: in high-yielding cows, circulating levels of growth hormone and corticosteroids are raised, while those of insulin and IGF-I are lowered (Breier et al., 1986; Rutter et al., 1989; Richards et al., 1991). The decrease in circulating insulin may cause the decrease in IGF-I, despite the increase in growth hormone, through a reduction in hepatic growth hormone receptor expression. All of these hormone changes may cause a delay in CLA, as all have effects on follicle development.
Selection for high yield, in addition to reducing overall energy balance, liveweight, and BCS (Veerkamp and Koenen, 1999), leads to alterations in circulating levels of reproductive hormones. For instance, in severe cases of negative energy balance, gonadotrophin secretion is blocked and ruminants become anoestrus (Jolly et al., 1995). It appears that gonadotrophin secretion is affected by selection for milk yield, since Gong et al. (2000) demonstrated a decreased LH pulse frequency and a decreased response to GnRH in cows selected for high yield. Furthermore, Royal (1999) and Royal et al. (2000c) provided evidence to suggest that a substantial proportion of phenotypic variation in the LH response to a GnRH challenge in prepubertal Holstein-Friesian heifers is due to additive genetic variation (h2 = 0.51). Since LH pulse frequency is one of the critical factors inducing onset of ovarian cyclicity postpartum (Jolly et al., 1995), this may be significant in the aetiology of subfertility. Hence, this trait may be of potential interest to dairy cattle selection programs to improve fertility. In addition to LH, other reproductive hormones (FSH and prolactin) may also be involved. Leptin may have a role in mediating these effects because it can act directly to induce the release of LH and FSH from bovine pituitary glands (Liou et al., 1997; Yu et al., 1997). Recent work by Kadokawa et al. (2000) has reported a phenotypic correlation (0.83) between the interval from parturition to the leptin nadir and the interval to first postpartum ovulation, suggesting that a delay in leptin recovery postpartum is associated with an extended CLA. There are no reports to our knowledge suggesting differences in leptin levels between high and low genetic merit dairy cattle, and a recent analysis by Royal et al. (2002a) found that additive genetic variation was not responsible for any variation observed in plasma leptin concentrations in prepubertal Holstein-Friesian heifers. However, they concluded that since high heritabilities were obtained in pigs at later stages of maturity (Cameron et al., 2000), genetic variation in dairy heifers may increase following expression of different genes or following metabolic challenge later in life.
The current work to investigate the relationships between linear type traits and CLA is relevant because linear type traits feature in most breeding programs and are good predictors of live weight and condition score (Veerkamp and Brotherstone, 1997). A number of the regression coefficients between fertility (as measured by CLA) and linear type traits reported in this study are of similar magnitude and direction to those published for CI and linear type traits by Pryce et al. (1998, 2000) and Dadati et al. (1986). The current results suggest that "frailer" cows (e.g., those possessing one or more of the following characteristics: angular, thin, narrow-chested, narrow-rumped, or high-pinned) on average will have reduced fertility, as measured by extended CLA. Cows with the opposite attributes are often referred to as "stronger," and such cows appear to have better fertility, which is in agreement with Dadati et al. (1986).
This work thus supports and furthers that of Pryce et al. (2000) and Veerkamp et al. (2000), showing that average BCS and some type traits may also be used to predict fertility, measured by CLA in this study, an important prerequisite for high fertility.
The genetic regression of lnCLA on PTA CI was moderate and favorable, appropriate to a genetic correlation of 0.36. Although CLA is free from management decisions, CI proofs may be biased by many managerial influences. Therefore, it is possible that the regression may be either inflated or deflated. For example, interval to first service in high genetic merit cows may be extended for management reasons (over and above the ra with yield) to prolong lactation in comparison to their lower-yielding herd mates. If this were the case, sires of high genetic merit cows may be given a lower ranking for CI than their true breeding value, that is, their PTA would be biased upwards. The regression of lnCLA on sire PTA for CI was affected by a small number of cows in the sire groups of PTA
6.5 d, which had a relatively low CLA. This may reflect the inclusion in the database of some cows with exceptionally long CI for management reasons, in view of the relatively long (600-d) restriction applied to CI in estimating (co)variance compounds for this parameter. Following exclusion of the 46 animals in this range (PTA
6.5), the regression of lnCLA on sire PTA for CI was 0.021 (SE = 0.010; P < 0.05) and appropriate to a genetic correlation of 0.69.
Alternatively, cows that are truly highly infertile will fail to have a subsequent calving due to culling for infertility. Under these circumstances, only the best of the sires daughters will have a CI and the bull PTA for CI may be biased downwards. This has been partially corrected for when estimating CI PTA used in this analysis by using bivariate models for the estimation of PTA values that include traits such as yield to help correct for missing data (Pryce et al., 2000). Therefore, we would predict that in this particular analysis, it is more likely that the regression of lnCLA on PTA CI is biased downwards.
Genetic regressions of lnCLA on updated PTA for milk, fat, and protein yield were positive and agreed with those reported previously by Royal et al. (2002b) and with estimates by Veerkamp et al. (2000). Cows with the genetic potential to produce higher yields tend to have an extended CLA. Therefore, the drive towards milk yield has resulted, in part, in reduced fertility. However, although the genetic relationships between fertility and milk yield are unfavorable, they are not 1, so selection for fertility using a designed index does not have to reduce milk yield or compromise selection for other traits of economic importance (Pryce et al., 2000b also demonstrates this).
Important implications of these findings are that, because CLA has a moderate heritability and is measurable in all animals rather than only those that complete a lactation and calve again, it may be possible to produce more accurate sire PTA for fertility (compared to PTA for traditional fertility traits such as CI) while using smaller progeny test groups. Therefore, CLA should be considered for inclusion in a fertility index.
This research has added to our understanding of CLA. A substantial component of the observed variation in CLA is an attribute of the individual and is genetic in origin (i.e., not simply a result of chance events experienced by the cow and errors in measurement). The observed variation is not only phenotypically correlated to traditional measures of fertility (Darwash et al., 1997a; Royal et al., 2002), but we now have the first evidence of a potential genetic correlation between CLA and CI, a traditional measure of fertility. The genetic relationship of CLA to other components of a proposed fertility index, such as production and type (as investigated here), appears to be substantial and significant, which has positive implications for the evaluation of CLA as a candidate trait for fertility in selection indexes. Furthermore, with the availability of on-line milk progesterone monitoring systems (likely within the next 5 yr) providing rapid progesterone profiles, the measurement of endocrine fertility parameters, and hence, estimation of breeding values using these measurements, will become more cost effective and more incisive. Since CLA appears to be a good predictor of average BCS, endocrine measurements of fertility may reduce the need for routine condition score assessment. In the meantime, type, production, average BCS, and other traditional measures of fertility such as CI should be used to halt or at least reduce the genetic decline in fertility.
| CONCLUSIONS |
|---|
|
|
|---|
| APPENDIX |
|---|
|
|
|---|
![]() |
where fixed effects are:
| µ | = | overall mean
|
| Lj | = | lactation number (i = 1 to 9)
|
| Hj | = | herd (j = 1 to 9)
|
| Yk | = | year of calving (k = 1995 to 1998)
|
| Sl | = | season of cavling (l = 1 to 4)
|
| Dm | = | diet (m = 1 to 23)
|
| Un | = | uterine infection (n = 0 or 1)
|
| Ro | = | retained placenta (o = 0 or 1)
|
| Cp | = | dystocia (p = 0 or 1)
|
| BqXq | = | regression variable (q = 1 or 2) where:
|
| X1 | = | percentage of Holstein genes, and
|
| X2 | = | sire PTA for the trait to be investigated
|
and where random effects are:
| HYkj | = | herd-year interaction (N(0, HY2)
|
| HSjl | = | herd-season interaction (N(0, HS2)
|
| YSkl | = | year-season interaction (N(0, YS2)
|
| HYSkl | = | herd-year-season interaction (N(0, HYS2)
|
| Ar | = | breeding value (N(0, A2A) where A is the numerator relationship matrix of cows available in the data
|
| Pr | = | individual (N(0, µ2)
|
| Eijklmnopqrs | = | error (N(0, E2)
|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication January 11, 2002. Accepted for publication March 27, 2002.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. M. Bewley and M. M. Schutz Review: An Interdisciplinary Review of Body Condition Scoring for Dairy Cattle Professional Animal Scientist, December 1, 2008; 24(6): 507 - 529. [Abstract] [PDF] |
||||
![]() |
J. J. Windig, B. Beerda, and R. F. Veerkamp Relationship Between Milk Progesterone Profiles and Genetic Merit for Milk Production, Milking Frequency, and Feeding Regimen in Dairy Cattle J Dairy Sci, July 1, 2008; 91(7): 2874 - 2884. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Dal Zotto, M. De Marchi, C. Dalvit, M. Cassandro, L. Gallo, P. Carnier, and G. Bittante Heritabilities and Genetic Correlations of Body Condition Score and Calving Interval with Yield, Somatic Cell Score, and Linear Type Traits in Brown Swiss Cattle J Dairy Sci, December 1, 2007; 90(12): 5737 - 5743. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. M. Chagas, J. J. Bass, D. Blache, C. R. Burke, J. K. Kay, D. R. Lindsay, M. C. Lucy, G. B. Martin, S. Meier, F. M. Rhodes, et al. Invited Review: New Perspectives on the Roles of Nutrition and Metabolic Priorities in the Subfertility of High-Producing Dairy Cows J Dairy Sci, September 1, 2007; 90(9): 4022 - 4032. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Hayhurst, M. K. Sorensen, M. D. Royal, and P. Lovendahl Metabolic Regulation in Danish Bull Calves and the Relationship to the Fertility of Their Female Offspring J Dairy Sci, August 1, 2007; 90(8): 3909 - 3916. [Abstract] [Full Text] [PDF] |
||||
![]() |
K.-J. Petersson, B. Berglund, E. Strandberg, H. Gustafsson, A. P. F. Flint, J. A. Woolliams, and M. D. Royal Genetic Analysis of Postpartum Measures of Luteal Activity in Dairy Cows J Dairy Sci, January 1, 2007; 90(1): 427 - 434. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Z. Caraviello, K. A. Weigel, M. Craven, D. Gianola, N. B. Cook, K. V. Nordlund, P. M. Fricke, and M. C. Wiltbank Analysis of Reproductive Performance of Lactating Cows on Large Dairy Farms Using Machine Learning Algorithms J Dairy Sci, December 1, 2006; 89(12): 4703 - 4722. [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] |
||||
![]() |
T. van der Lende, L. M. T. E. Kaal, R. M. G. Roelofs, R. F. Veerkamp, C. Schrooten, and H. Bovenhuis Infrequent Milk Progesterone Measurements in Daughters Enable Bull Selection for Cow Fertility J Dairy Sci, November 1, 2004; 87(11): 3953 - 3957. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. D. Dechow, G. W. Rogers, L. Klei, T. J. Lawlor, and P. M. VanRaden Body Condition Scores and Dairy Form Evaluations as Indicators of Days Open in US Holsteins J Dairy Sci, October 1, 2004; 87(10): 3534 - 3541. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. H. van der Waaij A resource allocation model describing consequences of artificial selection under metabolic stress J Anim Sci, April 1, 2004; 82(4): 973 - 981. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. P. Berry, F. Buckley, P. Dillon, R. D. Evans, M. Rath, and R. F. Veerkamp Genetic Parameters for Body Condition Score, Body Weight, Milk Yield, and Fertility Estimated Using Random Regression Models J Dairy Sci, November 1, 2003; 86(11): 3704 - 3717. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. P. Berry, F. Buckley, P. Dillon, R. D. Evans, M. Rath, and R. F. Veerkamp Genetic Relationships among Body Condition Score, Body Weight, Milk Yield, and Fertility in Dairy Cows J Dairy Sci, June 1, 2003; 86(6): 2193 - 2204. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |