|
|
||||||||

* Department of Population Medicine, and
Department of Animal and Poultry Science, University of Guelph, Guelph, Ontario, Canada, N1G 2W1
1 Corresponding author: rwalsh{at}uoguelph.ca
| ABSTRACT |
|---|
|
|
|---|
Key Words: anovulation time to pregnancy ketosis progesterone
| INTRODUCTION |
|---|
|
|
|---|
Following calving there is a normal period of anestrus during which uterine involution occurs. Thereafter, coordinated activity of the hypothalamicpituitary axis and ovaries reestablishes, culminating in ovulation and reinitiation of a regular cycle of estrus and ovulation (Wiltbank et al., 2002). Recent research has suggested that approximately 20% of dairy cows have neither displayed estrus nor ovulated before the start of the breeding period (Opsomer et al., 2000; Moreira et al., 2001). Ovulation synchronization protocols permit tighter control of time to first service; however, the probability of pregnancy in cows failing to return to a predictable ovulatory cycle is significantly reduced (Gumen et al., 2003).
There are many factors that influence time to first ovulation, and whether normal estrus and ovulatory cycles are maintained thereafter, including breed, parity, season, BCS, peripartum disease, and a variety of nutritional effects. Attention has focused on the association between postpartum negative energy balance and the resumption of cyclicity in dairy cows (Butler, 2000). Some studies were conducted in one herd at a time, in a single season, or have focused on characterization of the return to ovulatory cycles rather than identification of potential risk factors (Lamming and Darwash, 1998). The objectives of the current study were to estimate the cow-level prevalence of anovulatory condition in Ontario dairy herds, and identify potential herd- and cow-level risk factors for the condition.
| MATERIALS AND METHODS |
|---|
|
|
|---|
, agreed to enroll eligible cohorts of cows at the scheduled herd visit and take milk samples on the day of PGF2
injection, but before the scheduled PGF2
was administered.
Cows and heifers were assigned a BCS using a 5-point scale (Ferguson et al., 1994) in the 2 wk before calving and again between 52 and 66 DIM. Analysis was performed using BCS as an indicator variable, such that cows were classified as thin (
2.5), normal (
2.75 to
3.5) or fat (
3.75) at each of the 2 sampling periods. To account for differences in BCS classification between veterinarians, the change in BCS rather than the absolute BCS value was offered as a covariate associated with identification as anovulatory (Hady et al., 1994). Loss of greater than 1 BCS unit in early lactation was used to identify animals that suffer from severe negative energy balance (Loeffler et al., 1999). Parity and calving date were recorded. Calving date was subsequently classified into season of calving (fall = September through November; winter = December through February; spring = March through May; summer = June through August). The occurrence of dystocia (hard pull or veterinary-assisted calvings), twins, retained placenta (failure to pass fetal membranes by 24 h after parturition), and milk fever (cow down < 24 h after calving with no other illness) were recorded within 24 h of parturition. Milk BHBA concentration was determined using a previously validated milk test (KetoTest; Elanco Animal Health; Guelph, ON; Geishauser et al., 2000) between the second and eighth day postcalving and again 7 d later. Additional health records were maintained using a combination of on-farm paper and computer-based recording systems. Recorded disease events included metritis (severe infection of the uterus, cow off feed with a fever > 39.5°C, and foul-smelling uterine discharge), clinical ketosis (off feed and milk BHBA concentration > 500 µmol/L), displaced abomasum (diagnosed by the farm veterinarian), and clinical mastitis (visibly abnormal milk). Recent studies (Garbarino et al., 2004; Hernandez et al., 2005) indicated that moderate to severe lameness was associated with commencement of ovarian function and reproductive performance. Therefore, when cows were restrained for collection of milk samples, cows with an arched back while standing (corresponding to moderate to severe lameness; i.e., score 3 or 4 from Sprecher et al., 1997) were classified as lame, and cows without a back arch while standing were classified as nonlame. Skim milk progesterone concentrations were measured twice, 14 d apart (range 10 to 16), with the first sample taken between 37 and 51 DIM and the second sample between 52 and 66 DIM. Milk samples (approximately 40 mL) were collected by forestripping at the morning milking into vials containing a bacteriocidal tablet (6 mg of 2-bromo-2-nitropropane-1,3-diol and 0.3 mg of pimaricin). Samples were accumulated at the respective veterinary clinic and stored at 4°C. Milk samples were shipped to the University of Guelph between Monday and Wednesday each week to ensure samples were received within the same week. After receipt, whole milk was stored at 4°C until aspiration of the fat layer. A 2-mL sample of skim milk was stored at 20°C until progesterone concentrations were determined using a previously validated solid-phase RIA procedure (Gowan and Etches, 1979) and the Coat-a-Count system (Diagnostic Products Corporation, Los Angeles, CA). Animals were classified as anovulatory if the progesterone concentration was less than 1 ng/mL in both skim milk samples. The limit of detection of the progesterone assay was 0.12 ng/mL, and the inter- and intraassay coefficients of variation were 18.2 and 3.9%, respectively.
Reproductive records including insemination, pregnancy, and culling dates were retrieved from electronic files maintained by producers and Canwest DHIA every month. Herd file backups were obtained and stored permanently at 6-mo intervals to ensure access to all individual cow data. Pregnancy diagnosis was performed at 2-wk herd visits with the earliest time for insemination to pregnancy exam being 28 d, but the majority performed at 35 d or later. Pregnancy loss was identified at routine reconfirmation or if a cow was recorded as pregnant to an insemination performed after a previous diagnosis of pregnancy. A herd questionnaire was administered within 6 mo of herd enrollment to encourage accurate data recording, reconfirm the experimental protocol, and evaluate the impact of herd management on their respective performance. The questionnaire focused on cow management through the dry period up to and including the time of first insemination. Milking frequency (2 vs. 3 times/d), use of ovulation synchronization protocols for first insemination (> 90% of the herd vs. < 50%), hoof-trimming frequency, and footbath protocols were recorded.
General information including barn type (free-stall vs. tie-stall), stall dimensions (width and length), stall surface (pasture mat, concrete, or sand), and bedding material (sand, shavings, or straw), and the number of stall rows were recorded. Alley surface (rubber vs. smooth or ridged concrete) was recorded. Heat-abatement strategies were recorded for each herd without attempting to track variation in ambient temperatures experienced by individual herds.
Data were recorded concerning the feeding program used, the number of ration changes occurring between dry-off and 60 DIM, and any differences in management between primiparous and multiparous animals. The frequency of ration evaluation, forage DM determination, and number of times feed was delivered per day were recorded. Additionally, linear bunk space and water space per stall were measured. These variables were recorded for both the milking herd and dry cows.
Data Management and Statistical Analysis
All analyses were performed using Intercooled STATA 9.1 (STATA Corporation LP, 2005). Sample size estimates based on an estimated prevalence of anovulatory condition of 20%, an allowable error of 5% in the estimated prevalence, and cluster sampling suggested that a minimum of 750 cows were needed (Dohoo et al., 2003). Descriptive statistics were used to determine if cows that completed the sampling protocol differed from those cows calving over the study period that failed to complete the protocol. For animals from which 2 milk samples were collected, the time between the second milk sample and insemination was evaluated to estimate the potential for misclassification of anovulation (i.e., low progesterone associated with estrus resulting in false classification as anovulatory of animals that began to cycle between samples). Low progesterone (< 1 ng/mL) in both milk samples was the definition of an-ovulation; as such, cows inseminated within 5 d of the second milk sample were classified as anovulatory. In addition, milk collection dates were cross-referenced with reproductive hormone injections to avoid confounding of progesterone if prostaglandin was given (in deviation from the protocol) within 10 d of the milk collection date.
Distributions of specific variables were evaluated among animals that were classified as anestrus and those that were cycling. Accounting for correlation of cows within herds, the season of calving and parity (potential risk factors for anovulation) were evaluated one at a time. Continuous independent variables were offered to a multivariable logistic regression model unadjusted to assess their linear relationship. In instances where the linearity assumption was violated, these variables were evaluated using multiple linear relationships (splines) or were transformed into categorical variables. Risk factors that had univariate associations (P < 0.10) were subsequently evaluated through a backward stepwise multivariable logistic regression model using a generalized estimating equation to account for correlation of cows within herds. Variables were removed from the model if the associated P > 0.10, and potential 2-way interactions were evaluated if both main effects were significant.
Time from calving to first insemination and to pregnancy was estimated using a Kaplan-Meier life-table to account for enrolled cows that were not inseminated or not diagnosed pregnant before the end of the study. Kaplan-Meier (product limit) survival function estimates were used to calculate crude associations of an-ovulation with median time from calving to first insemination and median time to pregnancy, as well as to generate a graph of cumulative pregnancy risk over time. The restricted median time to event and 95% confidence intervals were determined as the first time the survival function was less than or equal to 0.50. The effect of anovulation on time to pregnancy was analyzed using multivariable survival analysis using Coxs proportional hazards regression. The Cox proportional hazard model measures the probability of pregnancy per unit time expressed as a hazard ratio. Robust standard errors were generated to account for clustering.
The main effect of interest was the impact of anovulation on time to first insemination and pregnancy. In addition to the main effect of cyclic status, the effects of parity, season of calving, periparturient disease, and milk BHBA concentration in each week postpartum were offered to the model as covariates. Estimation of the effect of cyclic status on time to first insemination was evaluated separately in herds systematically using an ovulation synchronization protocol and timed AI for first service. These models were built by manual backward stepwise elimination, removing the covariate with the largest probability value at each step, until all covariates were significant at P < 0.10. If 1 level of a categorical variable was significant, all levels were retained. Significant differences between class variables were assessed using a Wald postestimation test. Two-way interactions were assessed among significant variables by generating a cross product term. If the interaction term was significant at P
0.05, both fixed effects and the interaction term were retained in the model. The assumption of proportional hazard over time was assessed graphically using log-log plots. Violations of the proportional hazard assumption were corrected by allowing the covariate to vary with the log of time (Dohoo et al., 2003). This correction retains the standard hazard rate interpretation; however, it requires the calculation of the hazard rate at specific points in time.
Herd-level risk factors were investigated modeling the herd-specific seasonal prevalence of anestrus. Season was dichotomized into fall vs. all other seasons. This was done to account for different prevalence of anovulation in fall-calving cows. In addition, this division separated cows that had a dry period on pasture and captured a period with an increased number of calvings. Linear regression was used, accounting for the effect of clustering at the herd level. Herd size, barn design (free-stall vs. other), number of rows, stall design [large (> 127 cm wide and > 183 cm long) vs. small], linear feedbunk space per stall in the milking herd and dry groups, seasonal pasture use for the milking herd and dry cows, lameness prevention strategies, heat-abatement strategies, nutritional evaluation frequencies, feed delivery frequency, number of ration changes from dry-off to 60 DIM, milking frequency, and use of ovulation synchronization protocols were considered as covariates. Herd performance variables including proportion of cows culled before 60 DIM and herd prevalence of subclinical ketosis were evaluated in the first and second weeks.
| RESULTS |
|---|
|
|
|---|
Of the 18 herds included in the analysis, 5 reported no use of timed AI (TAI) protocols (31% of cows), 7 herds applied a TAI protocol to animals failing to display estrus by a herd-specific time postpartum (47% of cows), and the remaining 6 herds reported application of a TAI protocol for greater than 90% of animals at first AI (21% of cows). For analysis of time to first insemination, herds using observed estrus and occasional use of TAI protocols were combined and compared with herds enrolling > 90% of cows into an ovulation synchronization and TAI program.
A total of 266 animals were inseminated before the second milk sample being collected. Nine percent of these (24/266) were classified as anestrous based on milk progesterone concentrations. Of the apparently noncycling cows, 75% (18/24) were from herds that used estrus detection to identify cows for first insemination. Fifteen apparently noncycling cows (progesterone < 1 ng/mL in both milk samples) were excluded from the analysis because of hormonal therapy within 10 d before a milk sample.
Forty-three percent of cows had circulating progester-one concentrations < 1 ng/mL at the first milk sample. At the second milk sample 37% of cows had circulating progesterone concentrations less than 1 ng/mL. The overall prevalence of anovulation was 19.5% [95% confidence interval (CI): 17.4 to 21.7%]. A large amount of herd-level variation was evident [5% (95% CI: 1.0 to 13.7%) to 44.1% (95% CI: 33.2 to 55.3%); Figure 1
).
|
third: 33%; P = 0.47) and previous days dry (mean multiparous animals only = 74.1; 95% CI: 71.0 to 79.2; P = 0.58) of animals meeting the inclusion criteria relative to animals lost to follow up was not different. Peri-parturient disease information was more likely to be recorded in animals that met all inclusion criteria than in animals that did not. The culling pattern was weighted earlier in lactation for cows that did not have 2 milk samples collected relative to cows that met all inclusion criteria.
In the population of animals available for analysis, 8.4% of calvings were classified as a hard pull or required veterinary assistance. There were 49 (3.6%) twin calvings. Periparturient diseases included milk fever (2.4% among lactation
2), retained placenta (8.8%), metritis (7.9%), clinical ketosis (3.1%), and displaced abomasum (3.1%). Before 60 DIM, 8.9% of animals had a clinical case of mastitis and a lameness event was reported in 6.5% of cows.
Milk BHBA concentrations were determined in each of the first and second weeks postpartum in 1,041 animals. Twenty-six percent of animals (n = 1,041; 95% CI: 23.6 to 29.1) were subclinically ketotic (SCK; milk BHBA
100 µmol/L; Geishauser et al., 2000) in the first week postpartum. There was large variation in prevalence between herds (5.1 to 82.3%). In the second week postpartum 25.5% (95% CI: 22.8 to 28.2) of animals were SCK with equally large variation between herds (5.5 to 84.6%).
Cows experiencing dystocia or calving with twins were significantly more likely to be classified as anovulatory in a model accounting for clustering of cows within herd, season of calving, and parity, and in a multivariable logistic regression model including all significant factors. Similarly, retained placenta and displaced abomasum were significantly associated with anovulation, whereas lameness significantly increased the risk of anovulation when accounting for the effects of parity, season, and herd (Table 1
). Anovular cows were significantly more likely to be SCK in both the first and second week. But, after accounting for the effect of periparturient disease and milk production, animals with milk BHBA
100 µmol/L between 2 and 8 DIM were 1.5 times more likely classified as anovular. The probability of anovulation decreased by 2% for every 100-kg increase in the first 305-d, mature-equivalent milk projection between 5,000 and 9,500 kg. The risk of anovulation (17.2%, 95% CI: 14.5 to 19.8) did not change as first 305-d, mature-equivalent milk projection increased above 9,500 kg. The prevalence of an-ovulation changed with season. Animals calving in the fall were 46% less likely to be anovular than cows calving in the winter, whereas animals calving in the spring were significantly more likely to be anovulatory relative to animals calving in the fall. Cows calving in the summer were at an intermediate risk of anovulation between cows calving in the fall and those calving in the spring (Table 1
). Days dry, the fat to protein ratio at first DHIA milk test, BCS at calving, BCS at 60 DIM, and the change in BCS were not significantly associated with the risk of being classified as anovular.
|
|
|
|
|
The herd prevalence of anovulation increased by 2.1% for every 10% increase in the herd prevalence of SCK in the first week postpartum. The prevalence of SCK was 15.4, 38.6, and 29.8% in 2-row free-stalls, 3-row free-stalls, and tie-stalls, respectively. The proportion of cows culled before 50 DIM, the use of pasture during the summer, and all possible 2-way interactions were not significant. Residual analysis identified one leverage herd. Exclusion of this herd did not alter the coefficient estimates but it did alter the significance of the Wald test for ketosis in the first week.
| DISCUSSION |
|---|
|
|
|---|
The prevalences of peripartum diseases and SCK reported were similar to other observational studies in Ontario (McLaren et al., 2005). Dystocia significantly increased the risk of anovulation. This dichotomization of the 4-category calving score (unassisted, easy pull, hard pull, and surgery) is in agreement with other studies in which a composite variable including dystocia, stillbirth, and twin calvings was strongly associated with increased risk of anestrus (Opsomer et al., 2000; Shrestha et al., 2004a).
Peripartum disease events, including retained placenta, displaced abomasum, and lameness, occurring before the end of the voluntary waiting period significantly increased the probability of cows experiencing anovulation. Accounting for the time from calving to disease event, neither clinical ketosis nor mastitis was associated with the risk of failure to cycle. Another study found that after accounting for the impact of season of calving and parity, clinical ketosis increased the risk of anovulation more than 10-fold; however, this association was no longer significant after accounting for other risk factors such as dystocia and abnormal vaginal discharge (Opsomer et al., 2000).
Subclinical ketosis in the first week after calving significantly increased the risk of delayed ovulation. The association between the magnitude of negative energy balance and delayed ovulation is well established (Staples et al., 1990; Butler, 2000). A direct association between increased milk BHBA and delayed onset of luteal activity has not previously been reported; however, Cook et al. (2001) reported decreased probability of pregnancy at first insemination and increased time to pregnancy associated with hyperketonemia. The impact of SCK must be interpreted in light of when ketones were measured. In the present study, SCK in either wk 1 or 2 was associated with anovulation, but when other significant variables were accounted for, only SCK in wk 1 postpartum was associated with anovulation. The mechanism of this is not clear, but the peak incidence of SCK occurs in wk 1 of lactation (Duffield et al., 1998), and cows that are ketotic in wk 1 postpartum may be in more severe maladaptation to negative energy balance that cows that do not experience SCK until wk 2.
Body condition score before calving and at 60 ± 7 DIM, and the change between these 2 measurements were not significantly associated with anovulation in this study. Body condition score, and change in BCS had no association with reproductive performance measured by time to first observed estrus, first insemination, and time to pregnancy, although cyclicity was not measured (Ruegg and Milton, 1995). Conversely, a significant association between BCS and BCS loss between calving and first insemination was reported in several studies (Loeffler et al., 1999). The present results do not refute an association between BCS and risk of anovulation as reported by Lopez et al. (2005), but suggests that measurement of circulating BHBA concentrations in early lactation may provide additional information on the risk of anovulatory conditions.
We identified no association between duration of the dry period (74.1 d) and subsequent return to an ovulatory cycle. Others have suggested that cows with a dry period greater than 77 d were 2.9 times more likely to suffer prolonged postpartum anovulation than animals dry 64 to 70 d (Opsomer et al., 2000). The first mature-equivalent milk projection was included as a potential risk factor for anestrus because this value is consistently recorded for each animal before the end of the voluntary waiting period. The mean first projection was 9,500 kg. In animals below the mean (5,000 to 9,500 kg), the risk of anovulation decreased by 2% for each 100-kg increase in first projection. A first projection > 9,500 kg was not associated with a change in risk of anovulation, which is consistent with the results of a smaller study (Lopez et al., 2005). Animals with a low first projection failed to efficiently make the transition into lactation.
Animals calving in the fall were significantly less likely to be anovular relative to animals calving in the winter (P < 0.01) or spring (P = 0.05). Spring calving was associated with delayed commencement of luteal activity and is not unique to ruminants (Opsomer et al., 2000). The underlying mechanism cannot be determined in the current study; however, there were no differences based on use of pasture in the summer in the milking herd or during the dry period, suggesting that seasonal husbandry differences are not to blame. Furthermore, animals calving in late winter or spring will enter the breeding period during the hot months, which will negatively impact reproductive performance through increased probability of anovulation.
Anestrus was associated with an 8-d increase in the time to first service in herds in which most or all AI was based on detection of estrus. There were 42 anovular cows inseminated based on observed estrus before the second milk sample, suggesting that the estrus was not associated with an ovulation, cows were incorrectly identified as in estrus, or were correctly identified in estrus, but an irregular interovulatory interval followed this heat, consistent with abnormal ovarian activity. As anticipated, time to first AI was not different between cycling and anovular cows in herds using fixed-time insemination programs. Insemination in these herds was not associated with classification of cyclic status.
The probability of pregnancy at first insemination was 30.4, 20.3, 35.9, and 29.7% in cycling and anestrous cows bred at observed estrus, and cycling and anestrous cows bred at TAI, respectively. Animals bred at observed estrus were inseminated earlier than cows inseminated at TAI. The probability of pregnancy increases through the first 80 DIM independent of insemination protocol applied (Stevenson and Phatak, 2005).
The median time to pregnancy increased by 30 d in anovular cows relative to cycling cows (156 vs. 136 DIM). The impact of anestrus on time to pregnancy violated the proportional hazard assumption of semi-parametric survival analysis. To account for this, the impact of anestrus on time to pregnancy was permitted to vary with time. The instantaneous probability of pregnancy in anovulatory cows increased over time such that anovulatory cows returned to a daily probability of pregnancy similar to cycling animals around 165 DIM. The apparent risk of pregnancy loss was low (6%) and not different between cycling and anovular cows. Moreover, cows were typically diagnosed between 35 and 45 d after AI and were not systematically reexamined to measure pregnancy losses; losses were based on reexamination in some herds and pregnancy diagnosis to AI occurring subsequent to a previous pregnancy diagnosis. Therefore, these risks are likely underestimated and are not directly comparable to other studies in which pregnancy was first diagnosed earlier and was systematically followed for at least 28 d (Santos et al., 2004).
In this small sample of herds, free-stall barns had a lower prevalence of anovulation than did tie-stall barns. Within free-stall facilities, 3- or 6-row barns were more likely to have a higher prevalence of anovulation than 2-row barns. Independent of barn design, at the herd level, the prevalence of anestrus was linearly related to the prevalence of ketosis in the first week. The significance of the ketosis estimate was sensitive to the removal of a leverage herd; however, the direction and magnitude of the coefficient did not change. Given the sample size, this study was unable to identify the characteristics of free-stall barns that decreased the seasonal prevalence of anestrus or any potential interactions. Moreover, the finding that 3-row free-stall barns had a higher prevalence of anovulation than 2-row barns suggests that the potential reduction of linear bunk space per stall, and social dynamics surrounding feeding may account for some of the increase risk in anovulation (DeVries et al., 2004; Mentink and Cook, 2006). Further study of these herd-level effects in a large number of herds is warranted.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication June 6, 2006. Accepted for publication July 31, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. L. Bamber, G. E. Shook, M. C. Wiltbank, J. E. P. Santos, and P. M. Fricke Genetic parameters for anovulation and pregnancy loss in dairy cattle J Dairy Sci, November 1, 2009; 92(11): 5739 - 5753. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Bossaert, J. L. M. R. Leroy, S. De Vliegher, and G. Opsomer Interrelations Between Glucose-Induced Insulin Response, Metabolic Indicators, and Time of First Ovulation in High-Yielding Dairy Cows J Dairy Sci, September 1, 2008; 91(9): 3363 - 3371. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |