1. Introduction
Anestrous is one of the most common reproductive problems in cross-bred dairy cows in tropical countries. Prevalence estimates range from 20 to 50% [
1,
2] and result in subsequent financial losses [
3]. Anestrous is a multifactorial problem related to negative energy balance (NEB) during the transition period, heat stress, reproductive pathology, and low heat detection efficiency [
4]. Farm management interventions such as dietary management after calving [
5], reducing heat stress [
6], and improving estrus detection practices [
7] are commonly used to mitigate these factors.
The bovine estrous cycle is regulated by endogenous hormones produced by the hypothalamic-pituitary-gonadal axis. Exogenous reproductive hormones or their analogs (e.g. GnRH, PGF
2α, eCG, and progesterone) are commonly used to induce cyclicity in anestrous cows and help reproductive management in dairy herds. Various programs that use these hormones have been developed and can be divided into 1) estrus detection-based programs before artificial insemination (EAI), 2) fixed-time AI (TAI), and 3) combinations of EAI and TAI. The primary purposes of these hormone interventions among anestrous dairy cows are to increase the probabilities of estrous detection and insemination and to increase the probability of pregnancy within a desired period [
8]. The probability of pregnancy per cow eligible to be inseminated within the desired after treatment period is called pregnancy risk [
9]. The terms 'pregnancy risk' and 'pregnancy rate' were discussed in terms of appropriate terminology because they were not evaluated and expressed per unit of time [
9]. Therefore, the pregnancy rate could refer to pregnancy risk in some previous studies [
5,
10,
11]. Using hormone treatment can increase pregnancy risk among anestrous dairy cows [
5,
10]. Moreover, the TAI protocols allow all animals to be artificially inseminated without reliance on heat detection, and may be helpful in situations where heat detection efficiency is lacking or when estrous behavior is decreased [
5].
Pregnancy risk, after hormone treatment, has been shown to vary between different protocols [
11] and is also affected by several animal-related factors, such as the presence or absence of a corpus luteum (CL) at the initial time of hormone implementation, body condition score (BCS), parity, and days in milk (DIM) [
12]. Environmental factors such as heat stress influence pregnancy risk [
13]; pregnancy risk during the hot and humid season is typically lower than during the cool season [
14,
15]. Farm management factors such as estrus detection efficiency and proper AI practices also influence pregnancy risk [
16,
17].
Estrus induction and TAI have been shown to minimize estrus detection and insemination problems on dairy farms [
8]. Some farmers question the value of hormone treatment due to its higher cost than estrus detection followed by AI alone. The use of reproductive hormones, synchronization programs and TAI in anestrous bovine species (cows and buffaloes) in tropical countries has been studied and showed desired results in improving reproductive performance under tropical conditions [
11,
18,
19,
20,
21,
22]; however, information on the cost-efficiency of different hormonal treatment protocols is still limited.
Although several studies have been published evaluating the effect of hormonal treatments on anestrous bovines in recent years, it challenges the idea that most studies are experimental studies comparing difference protocols under the selected conditions or cohort studies on commercial farms in temperate regions. These may not reflect the real world situation in tropical countries when veterinarians applied treatment protocols on small-holder and semi-commercial farms where the performance of hormone treatment protocols might be affected by the factors mentioned above and veterinarian and farmer decisions for treatment.
Moreover, climate change is more likely to negatively affect the productivity, reproductive performance, and health of livestock in tropical countries and is becoming even more extreme than anyone predicted over the past 10 years. High environmental temperature has unfavorable effects on reproductive function among dairy cattle [
23]. Effective cooling systems, housing management, hormone treatment, and other reproductive management technologies have been applied to diminish climatic impacts on reproductive performance. However, anecdotal evidence showed that hormonal therapy appears to be less effective than expected under field conditions based on discussions with local Thai dairy veterinarians.
The present retrospective study is the first epidemiological study on hormone treatment protocols in anestrous dairy cattle in Thailand using clinical records from a veterinary teaching hospital. This study aimed to help farmers and clinicians in tropical countries make the most effective decision about hormone treatment in field conditions. The objectives of this study were: 1) to evaluate the performance of hormone treatment protocols on pregnancy risk; 2) to determine factors associated with pregnancy success after hormone treatment; and 3) to assess the cost-efficiency of hormone treatment protocols among tropical anestrous dairy cattle in Thailand.
2. Materials and Methods
The study is a retrospective observational study carried out by the Herd Health and Production Service (HHPS) Unit from the Large Animal Clinic of the Kasetsart University Veterinary Teaching Hospital, Nong Pho. This study was part of a project that aimed to improve cattle reproductive efficiency for small farmers in three provinces of the central-western region of Thailand (that is, Ratchaburi, Kanchanaburi, and Nakhon Pathom). Clinical records of anestrous dairy cows that received hormone treatments were obtained from the HHPS database from January to August 2017.
2.1. Dairy Farms and Management
Ratchaburi, Kanchanaburi, and Nakhon Pathom are in the tropical climate zone between 13 ° 50′ and 14 ° 10′ N and between 99 ° 50′ and 100 ° 00′ W. The average temperatures range from 21 ° C to 36 ° C and the average relative humidity range from 69% to 79% [
24]. It is a typical agricultural area containing 13% of the total population of dairy cows in Thailand [
25].
Thai crossbred Holstein dairy cattle, commonly implies roughly 87.5% Holstein-Friesian and 12.5% indigenous Thai crossbred dairy cattle [
26,
27], is the main population among participating dairy farms. The size of the herd was relatively small, ranging from five (small-holder herds) to 90 lactating cows (semi-commercial herds). Housing was tied-stall or loose-housing systems. Farmers fed their cows twice a day with a predominantly forage and concentrate-based diet, using roughage such as Napier grass (
Pennisetum purpureum), para grass (
Brachiaria mutica), and local crop by-products (eg, sweet and baby corn by-products, pineapple by-products, and rice straw hay). Forage and concentrate are commonly fed based on ad libitum forage supplemented with protein/energy concentrate using a 2:1 ratio of concentrate to kg of milk produced. The voluntary waiting period was generally set at 50-60 days. AI was the primary breeding method used. Estrus detection was performed on the basis of visual observation at least twice a day. The AI was usually carried out within 12-16 hours after standing heat was observed. The farms had monthly reproductive health checks performed by the Large Animal Unit, including diagnosis and treatment of reproductive disorders, reproductive cycle status and estrus monitoring, and pregnancy diagnosis. The pregnancy diagnosis was made by ultrasound or rectal palpation approximately 35 days after AI.
2.2. Data Management
Anestrous cases were defined as the absence of visible estrus detected by farmers for <30 days before the day of hormone treatment and diagnosed by clinicians. Furthermore, the anesthetized cases in this study had to meet all inclusion criteria: 1) lactating cows; 2) reproductive data were routinely recorded in the HHPS database; and 3) cows were treated with hormone treatment.
The following information was available for each cow on the day of hormone treatment: cow id, parity, BCS based on the scoring system of Furguson et al. [
28](< 3.00 or > 3.75 and 3.00-3.75), DIM, AI history (yes and no), history of hormone use (yes and no), season (Winter = 1 January to 2 March, Summer = 3 March to 25 May, and Rainy = 26 May to 31 August), herd id, herd size ( 20 lactating cows and > 20 lactating cows), housing system (tie-stall and loose housing), clinician experience (< 6 years and 6 years). Furthermore, an ovarian structure was collected at the initial time of hormone treatment and categorized as cows with CL present at the beginning of hormone treatment (cyclic cow) and cows without CL present (noncyclic cow). Data on hormone treatment protocols were also collected and explored for validity and reliability; incomplete and incorrect protocols were excluded.
We aimed to study the performance and cost-efficiency impact of single-treatment hormonal therapy in anestrous dairy cattle in Thailand. We did not aim to include analysis of repeated treatment protocols (that is, repeated treatments with hormonal therapy of the same animal in the same dry period) in this study. The hormone treatment protocols with less than 2% of the total records were excluded to assess common practices. The protocols were classified as 1) estrus detection before AI (EAI; Protocols 1-4, 9-10) and 2) fixed-time AI (TAI; Protocols 5-8) (
Figure 1)
Estrous expression, AI event, and pregnancy outcome after hormone treatment were collected to assess estrus detection risk, insemination risk, pregnancy risk, and pregnancy per AI. Cows were classified as nonpregnant and pregnant; They were considered as nonpregnant if they met at least one of these criteria: 1) cows were diagnosed as non-pregnant; 2) cows returned to estrus before pregnancy diagnosis; 3) cows were re-inseminated before pregnancy diagnosis; and 4) there was no AI event after treatment. The records were also excluded if pregnancy data was missing. Estrus detection risk (ER), the probability of estrus detection after treatment, was calculated by the number of cows detected in estrus after treatment, divided by the number of cows treated that were supposed to be detected in estrus after treatment. Insemination risk (IR), the probability of AI after treatment, was calculated by the number of cows inseminated after treatment, divided by the number of cows treated that were supposed to be inseminated after treatment. Pregnancy risk (PR), the probability of pregnant cows after treatment, was calculated as the number of cows diagnosed as pregnant with AI after treatment, divided by the number of cows treated that were supposed to be detected in estrus and inseminated after treatment. Pregnancy per AI was calculated as the number of cows diagnosed as pregnant with AI after treatment, divided by the number of inseminated cows after treatment.
2.3. Cost-Efficiency for Hormonal Treatment
Total treatment cost was calculated as the sum of hormone treatment cost, veterinary service fees, and frozen semen and AI costs. Costs that were not part of treatment costs, such as costs associated with routine daily management, including labour and estrus detection, and monthly costs of herd health service, which are standardized per cow in the study area, were not included in the analysis. All costs were estimated using the average price in the Thai market quoted here in US dollars. The exchange rate from May 8, 2020 was used (1 USD = 32.28 THB) for the whole analysis [
29]. The cost of hormone drugs per single dose used in the dairy population studied were
$6.20 for buserelin (Receptal
®) and
$2.48 for cloprostenol (Estrumate
®). The veterinary service fee was estimated to be
$3.10 per cow for the entire treatment program. The price for frozen semen and artificial insemination, calculated based on the regular price reported in the study area, was
$3.10 per service. Based on these data, the total cost of treatment per cow and the total cost per pregnancy were calculated for TAI and EAI among cyclic and noncyclic dairy cows. The total treatment cost per cow per protocol was calculated as the total treatment costs for a particular protocol (see
Figure 1) on all farms, divided by the number of cows treated in that protocol. The total cost per pregnancy was calculated similarly to the total cost of treatment per protocol divided by the number of pregnancies achieved in that protocol.
Figure 1.
Schematic diagram showing hormone protocols that were treated among 294 anestrous dairy cows from 65 herds under the Herd Health and Production Service Program, the Large Animals Unit, Veterinary Teaching Hospital, Nong Rho, Kasetsart University in 3 provinces of the central-western region of Thailand since January to August 2017; Protocol 1-4, 9-10 were recognized as hormone treatment with estrus detection before AI (EAI), and Protocol 5 to 8 were recognized as hormone treatment with fixed-time AI (TAI). TAI = fixed-time artificial insemination, AI = artificial insemination, GnRH= Gonadotropin-releasing hormone, PGF2α= Prostaglandin F2alpha.
Figure 1.
Schematic diagram showing hormone protocols that were treated among 294 anestrous dairy cows from 65 herds under the Herd Health and Production Service Program, the Large Animals Unit, Veterinary Teaching Hospital, Nong Rho, Kasetsart University in 3 provinces of the central-western region of Thailand since January to August 2017; Protocol 1-4, 9-10 were recognized as hormone treatment with estrus detection before AI (EAI), and Protocol 5 to 8 were recognized as hormone treatment with fixed-time AI (TAI). TAI = fixed-time artificial insemination, AI = artificial insemination, GnRH= Gonadotropin-releasing hormone, PGF2α= Prostaglandin F2alpha.
2.4. Statistical Analysis
Data analysis was carried out using the STATA version 16.0 statistical software package (Stata Corp., College Station, TX). The Chi-square test (P < 0.05) was used to compare ER, IR, and PR of two types of hormone protocols between cyclic and noncyclic dairy cows. A two-sample independent t-test (P < 0.05) was used to compare the total treatment cost, hormone cost, veterinary service fees, and frozen semen and AI costs of two types of hormone protocols between cyclic and noncyclic dairy cows.
The success of pregnancy was considered the outcome of interest for the multivariate logistic regression analysis. The analysis evaluated the performance of hormone treatment protocols on pregnancy success after treatment. All explanatory variables are shown in the causal diagram (
Figure 2). Descriptive statistics (mean, standard deviation, minimum, maximum, frequency, and percentage) were used to assess the characteristics of the outcome variable and all explanatory variables. Explanatory variables with missing values of more than 10% were excluded from the analysis. Continuous variables were assessed for linearity on a logit scale using the Lowess smoothing curve, Box-cox, and fractional polynomials [
30]. In addition, the correlation among the independent variables was assessed using Pearson's correlation, Kruskal-Wallis test, and the Chi-square test. Unconditional association was tested between explanatory variables and outcome, and variables that were unconditionally associated (P <0.020) with outcome were retained for use in the multivariate logistic regression model. A manual stepwise backward approach was applied to estimate the final multivariate model; The approach was initialized using all significant explanatory variables from the univariate analysis. The variables with the highest p-values (P > 0.05) were removed with the addition of previously removed variables. A variable was considered a confounder and retained in the final model when the removal variables altered the coefficient of the remaining variable by more than 20% and met the following criteria: 1) was associated with the outcome; 2) was not the consequence of the outcome; and 3) was not an intervening variable in the causal diagram. However, non-significant variables that were not confounders were permanently removed and the final multivariate mode was decided when all explanatory variables in the model showed a significance level (P <0.05). The biologically plausible interactions between the explanatory variables were tested and kept in the model if they were significant (P <0.05). The fit of the final model was tested using the Hosmer-Lemeshow Goodness-of-fit test. In addition, the sensitivity, specificity, and area under the ROC curve (AUC) of the model were calculated to estimate the predictive ability of the model. DeltaX
2 and Delta-beta were used to identify outliers. The outliers were removed from the data set and investigated for their effect on our result
3. Results
3.1. Study Population
A total of 489 hormone treatment records were obtained from 414 cows in 76 farms were obtained; reasons for treatment are presented in
Table 1. Thirty farms (39.5%) used hormone treatment to address reproductive problems before the study period. From this dataset, 369 records from 308 cows in 66 farms were identified as using hormone treatment therapy to address anestrous problems; 61 records were removed due to retreatment of the same cow. Subsequently, 14 records were removed because they did not meet the inclusion criteria and because of non-compliance protocols. 294 records were obtained from 294 cows on 65 farms treated with ten different protocols (
Figure 1). Additionally, three protocols (Protocols 8, 9, and 10) that amount to nine cows were excluded due to the few numbers of treatments (less than 2% of the total records). As a result, seven treatment protocols were considered in this study and classified as protocols with estrus detection before AI (EAI; Protocols 1-4) and protocols with fixed-time AI (TAI; Protocols 5-7). Six records were also deleted because pregnancy outcomes were missing. The final data set included 279 records of 279 cows on 64 farms. Thirty-eight farms (66% of cows), 22 farms (27% of cows), and four farms (7% of cows) were in Ratchaburi, Kanchanaburi, and Nakhon Pathom provinces, respectively. A total of 35 farms (62% of cows) were registered as employing a loose-housing system, and 29 farms (38% of cows) used tie-stalls.
3.2. Hormone Treatment Outcomes
Hormone treatment results are shown in
Table 2. The pregnancy risk for the TAI protocol was higher than for the EAI protocol, but pregnancy per AI did not differ significantly between the two protocols in cyclic and noncyclic cows.
3.3. Factors Associated with Pregnancy Success after Hormonal Treatment
The distribution of the level variables of interest and the univariate analysis can be found in
Table 3. According to the diagram in
Figure 2, the ovarian structure at the initial time of hormone treatment and season were considered confounding factors for the hormone treatment protocol in the final model. Although the ovarian structure at the initial time of hormone treatment was not related to pregnancy outcome in the univariate analysis, it was forced into the multivariate model to control for the known effect, explained by the previous study [
12]. As a result, the coefficient of the hormone treatment protocol was changed by more than 30%. Although the history of AI was related to pregnancy outcome in the univariate analysis, it was highly collated with the season. Moreover, only the coefficient of the season variable was affected and became nonsignificant when AI history and season remained in the final model. Therefore, the history of AI was considered the intervening variable for the season and was removed from the analysis. There was no interaction between explanatory variables found in this study. Farms and clinicians were not significant when accounting for random effects in the multilevel logistic regression model; therefore, a simple logistic regression model was used. The final results of the multivariate model for factors associated with pregnancy success are presented in
Table 4.
The Goodness of fit of the final model was evaluated using the Hosmer–Lemeshow test and revealed very weak evidence of lack of fit (P = 0.85). Sensitivity, specificity, and AUC were 0.70, 0.68, and 0.72, respectively. The cut-off point for the sensitivity and specificity measures was 0.18, resulting from a sensitivity-specificity plot. Evaluating outliers and influential observations did not reveal weaknesses in the model.
3.4. Cost-Efficiency for Hormonal Treatment
The total cost for the TAI protocol was higher than the EAI protocol in cyclic and noncyclic cows; however, the cost per pregnancy for the TAI protocol was lower than the EAI protocol in noncyclic cows (
Table 5).
4. Discussion
4.1. Hormone Treatment Outcomes
The PR for the TAI protocol was significantly higher than the EAI protocol in noncyclic cows; there was a trend towards significance in cyclic cows. However, when cows were inseminated after treatment, there were no differences in pregnancy according to AI between the EAI and TAI protocols for both noncyclic and cyclic cows, which was in accordance with previous studies [
31,
32]. The low PR for the EAI protocol in both cyclic and non-cyclic cows can be explained by the fact that the IR for the EAI protocol was significantly lower than that for the TAI protocol because the IR for the EAI protocol relies on ER that requires estrus detection by the farmer. In contrast, all cows that received the TAI protocol were inseminated after hormone treatment whether they showed estrous behavior or not. The results showed that the ER for the EAI protocol was low, especially in noncyclic cows, which can be explained by limited knowledge of estrus detection among Thai dairy farmers [
33]. Therefore, the majority of cows who received the EAI protocol lost the opportunity to be inseminated after hormone treatment. Moreover, GnRH analogs have a limited half-life which does not exceed 5-6 hours, which represents the main side effect of the GnRH analog [
21]. Therefore, its use in the EAI protocol for noncyclic cows may not be efficient for estrous behavior because the estrus induction efficiency varies and depends on the size and stage of ovarian follicles at the time of treatment [
21]. However, Amin and Said [
21] reported that adding chitosan to GnRH increases its efficacy, accelerates the growth of ovarian follicles, and improves hormonal insufficiency, resulting in an increased estrous induction rate and a conception rate. Therefore, GnRH-conjugated chitosan might be an appropriate choice to apply in the EAI protocol for noncyclic cows.
Estrus detection occurred more often in cyclic cows than noncyclic cows following the EAI protocol due to the higher ER. Previous research showed that cyclic cows display obvious expression of estrous behavior and a longer estrus duration than noncyclic cows [
32]. This phenomenon can be explained by the concept of the depth of anestrous; cows in deep anestrous conditions are more affected by many factors that have a deleterious effect on ovarian activity and are more difficult to treat [
34]. Noncyclic cows suffering from true anestrous or inactive ovaries had deeper anestrous than did cyclic cows and rarely showed estrous behavior after treatment. It was surprising that the ER for the TAI protocol was limited to about 0.3 in both cyclic and noncyclic cows, which is noticeably low. This might be because farmers may ignore estrus detection when TAI protocols are applied while paying more attention to cows treated with EAI protocols. Another reason could be the consequence of the last injection of GnRH before AI in the TAI protocol [
32] (
Figure 1). GnRH injection suppressed estrous behavior because GnRH directly inhibits natural estradiol synthesis from granulosa cells, resulting in an estradiol decline in blood concentrations [
35,
36,
37].
4.2. Factors Associated with Pregnancy Success after Hormonal Treatment
When multivariate logistic regression analysis was performed to evaluate factors associated with pregnancy success, cows that received the TAI protocol were about four times more likely to get pregnant than cows that received the EAI protocol. This is consistent with our study results mentioned above: the TAI protocol showed a higher proportion of pregnant cows than the EAI protocol in both noncyclic and cyclic cows. However, there is no interaction between the hormone treatment protocol and ovarian structure at the initial time of hormone treatment.
Cows with 3.00 ≤ BCS ≤ 3.75 were more likely to get pregnant than cows with BCS < 3.00 or BCS > 3.75. Consistent with the findings of the present study, previous research supports that cows with a BCS ≥ 3.00 are more likely to conceive and remain pregnant than cows with a BCS < 3.00 following hormone treatment [
38,
39]. Our study proposes that this could be the long-term adverse effect of NEB from the transition period [
40], especially in Thailand, where dairy farmers are not well trained in nutritional management [
41]. Unfavorable metabolic changes during NEB have many detrimental effects on the quality of oocytes and embryos, resulting in unsuccessful fertilization and pregnancy after insemination [
42]. Low BCS has an adverse impact on follicle development and reduces ovulation rate after hormone treatment [
43]. Furthermore, Madureira et al. [
44] reported that estrous expression was affected by BCS, since cows with BCS ≤ 2.50 expressed less intense estrous activity. This disturbed estrus detection efficacy when the EAI protocol was applied, resulting in estrus detection failure and loss of AI opportunities. On the other hand, the lower pregnancy success among cows with BCS > 3.75 compared to cows with 3.00 ≤ BCS ≤ 3.75 can be explained by the fact that cows with BCS > 3.75 can present insulin resistant complications; there is evidence that this effect results in reduced oocyte quality and impaired embryo development [
45].
The loose housing cows had better pregnancy success than the tie-stall cows after hormone treatment. This finding is consistent with previous research indicating that loose-housed cows showed better fertility than tethered cows; loose-housed cows have a better pregnancy rate after first A, the number of services per conception, calving to conception interval, and calving interval [
46,
47]. Additionally, loose-housed cows can freely show estrous behavior and provide greater opportunities to realize estrous behavior, while tie-stall housing systems greatly impair estrous expression because cows are restricted to individual stalls [
48]. Movement restriction in the tie-stall system might alter reproductive function due to physiological stress; tie-stall cows had a higher concentration of hair cortisol compared to loose-housed or free-stall cows [
49,
50].
The season negatively impacted pregnancy success after hormone treatment. Pregnancy probabilities increased from the rainy season to the summer season and reached the highest point in the winter. These findings agreed with the most recent Thai report that the probability of pregnancy in the winter season was higher than in the rainy season and the summer season after hormone implementation [
13]. A low proportion of pregnancies has also been reported in bovines after receiving hormone treatments during the rainy season in Thailand [
38]. Biological reasons for this can be explained by the effect of heat stress; cows usually face mild stress conditions in the winter season, after which they experience moderate heat stress during the summer and rainy seasons [
13]. Heat stress adversely affects fertility by altering reproductive physiology through hyperthermia, oxidative stress, hormone imbalance, and negative energy balance [
19,
51,
52].
The effects of heat stress result in reduced oocyte quality, decreased semen quality, and impaired embryo development, which contribute to poor estrous expression and infertility, especially in Thailand [
53]. Furthermore, our study result showed that pregnancy success reached the lowest point in the rainy season. Schüller et al. [
54] reported that the likelihood of dairy cows becoming pregnant is reduced by both short- and long-term heat stress. Long-term heat stress before the day of breeding leads to a reduced number of follicles, a reduced estradiol concentration in the follicles, and an earlier emergence of the dominant follicle [
55,
56]. In our study, we proposed that anestrous cows received hormone treatment in the rainy season; these cows experienced heat stress during the summer season, resulting in low pregnancy success after treatment. On the other hand, other environmental factors, such as rainfall, also affected pregnancy success. Heavy rain falling on the day before and after AI is related to low pregnancy success in dairy cows [
57]. Furthermore, heavy rain, strong wind, or high humidity indirectly affect pregnancy success due to suppressing estrous behaviors [
58], which will affect estrus detection efficiency and AI timing for cows receiving the EAI protocol.
The season and ovarian structure at the initial time of hormone treatment were also identified as confounding variables for the hormone treatment protocol. The reason for this is the clinician's decision, who were more likely to treat anestrous cows with the EAI protocol rather than the TAI protocol in winter, contrasting with the summer and rainy seasons due to the adverse effect of heat stress and economic concerns. Similarly, clinicians might also be aware of the low estrous expression ability in noncyclic cows [
31]. Therefore, far more cyclic cows were treated by the EAI protocol than noncyclic cows, and a majority of noncyclic cows were treated by the TAI protocol (
Table 2).
The findings of our study did have some limitations. First, enrollment in the study was selected from the hospital database using the anesthesizing definition and selection criteria. Therefore, there was the potential for selection bias to exist because the study population was not randomly selected from the target population. Second, the misclassification of the ovarian structure at the initial time of treatment may have occurred due to the misidentification between CL and follicles by clinicians and due to the recall bias by farmers when reporting estrous expression. All anesthetized cows were diagnosed and treated by six different clinicians in this study, and the cows came from different farms. However, due to the limited small sample in each category, the mixed model cannot account for farm and clinicians as a random effect. This study assumed that the sources of errors were equally distributed. Finally, the AUC value revealed that the overall diagnostic precision of the model was acceptable discrimination [
29]. The probability that the model correctly classified pregnancy success after hormone treatment was 72%.
4.3. Cost-Efficiency for Hormonal Treatment
The total average treatment cost per cow indicated that the TAI protocol was more expensive than the EAI protocol due to the higher cost of hormones and frozen semen and the AI service. The cost of hormones represented the majority of the treatment costs, the most expensive of which was GnRH. Ricci et al. [
59] reported that the GnRH expense had a greater impact on the net profit gain than prostaglandins. The average cost of frozen semen and AI service was higher than that of the EAI protocol. These findings also support the observation we initially reported that the cost of hormone therapy is a substantial component of total costs, probably contributing to the erroneous belief of non-adopter Thai dairy farmers that it may not be cost effective.
If only the total average treatment cost per cow is taken into account, fixed-time AI hormonal treatment appears more expensive than estrus detection alone. However, one must also take into account the resulting pregnancies that arise due to treatment. The cost per pregnancy for noncyclic cows receiving the EAI protocol was about twice as much compared to costs for the TAI protocol. The PR for the EAI protocol was lower than for the TAI protocol. However, the cost per pregnancy for cyclic cows receiving the EAI protocol was lower than the TAI protocol; the EAI protocol tended to have a lower PR than the TAI protocol. We propose that the TAI protocol is the most cost-efficient single-treatment hormonal therapy option for noncyclic cows, whereas the EAI protocol is the most cost-efficient option for cyclic cows.
We note that we were unable to collect a set of farm economic data as we would ideally have liked. Under field conditions, farmers rarely recorded milk production costs and income. Our calculations reflect hormone treatment costs without taking into wider consideration other production costs including labor costs, maintenance costs for non-pregnant cows (e.g., feed and feeding cost, and re-insemination cost), culling costs, and the cost of lost milk revenue due to extended periods of nonpregnancy.
5. Conclusions
Anestrous is a major reproduction problem in the central-western region of Thailand. This is the first study to integrate evaluation of hormone treatment outcomes, factors associated with pregnancy success after hormone treatment, and cost-efficiency analysis. These would help farmers and clinicians make the most effective decision to deal with the anestrous problem in dairy cows under tropical conditions, such as in Thailand. Pregnancy success was influenced by the hormone treatment protocol, BCS at the initial time of hormone treatment, housing conditions, and season. The TAI protocol seems to be a very useful method for solving anestrous problems in dairy cows due to high PR. But cost-efficiency analysis showed that the TAI protocol suits non-cyclic cows, whereas the EAI protocol suits cyclic cows due to the lower cost per pregnancy. Therefore, the hormone treatment decision to solve the anestrous problem should be based on the health of the animal, environmental conditions, farm and management conditions, and the cost effectiveness of treatment protocols.
Author Contributions
Conceptualization, T. Changtes, J. Sanchez, P. Arunvipas, T. Patanasatienkul, D. Hall, L. Heider and T. Rukkwamsuk; Methodology, T. Changtes, P. Arunvipas, P. Thammahakin, J. Jareonsawat and T. Rukkwamsuk; Software, T. Changtes, T. Patanasatienkul; Validation, T. Changtes, J. Sanchez, D. Hall and L. Heider; Formal analysis, T. Changtes, J. Sanchez, D. Hall and L. Heider and T. Patanasatienkul; Investigation, T. Changtes, T. Thammahakin, and J. Jareonsawat; Resources, J. Sanchez and T. Rukkwamsuk; Data curation, T. Changtes, T. Patanasatienkul, P. Thammahakin, J. Jareonsawat and T. Rukkwamsuk; Wring original draft, T. Changtes, and T. Patanasatienkul; Writing - review & editing, J. Sanchez, D. Hall, L. Heider, and T. Rukkwamsuk.; Visualization, T. Changtes; Supervision, J. Sanchez and T. Rukkwamsuk; Project administration, T. Rukkwamsuk; Funding acquisition, J. Sanchez and T. Rukkwamsuk.
Funding
This research was funded by the Faculty of Veterinary Medicine of Kasetsart University. The Canada-ASEAN Scholarships and Educational Exchanges for Development (SEED) provided a scholarship to Thitiwich Changtes to visit Canada for data analysis.
Institutional Review Board Statement
The study did not involve animals as an experimental setup. Data were from the routine services of the Herd Health and Production Service (HHPS) Unit from the Large Animal Clinic of the Kasetsart University Veterinary Teaching Hospital, Nong Pho. Required permission to conduct the study and to use the data has been received by the Faculty of Veterinary Medicine, Kasetsart University.
Informed Consent Statement
Not applicable.
Data Availability Statement
The data that support the findings of this study are available from the authors but restrictions apply to the availability of these data, which were used under permission from the Faculty of Veterinary Medicine, Kasetsart University. Data are, however, available from the authors upon reasonable request.
Acknowledgments
This study was mainly supported by a contribution from Faculty of Veterinary Medicine, Kasetsart University, and Canada-ASEAN Scholarships and Educational Exchanges for Development (SEED). We would like to express our gratitude to all bovine veterinarians and staff of the large animal unit of Kasetsart University Veterinary-Teaching Hospital Nong Pho for their responsibilities and hard work that provided us the valuable data reflecting the real outcomes of both the field-applied hormone program and the cost-efficiency for farmers.
Conflicts of Interest
The authors declare no competing interests.
References
- Khan, M.H.; Manoj, K.; Pramod, S. Reproductive disorders in dairy cattle under semi-intensive system of rearing in North-Eastern India. Vet. World 2016, 9, 512–8. [CrossRef]
- Regassa, T.; Ashebir, G. Major factors influencing the reproductive performance of dairy farms in Mekelle City, Tigray, Ethiopia. J. Dairy Vet. Anim. Res. 2016, 3, 145–9. [CrossRef]
- Kumar, P.R.; Singh S.K.; Kharche, S.D.; Govindaraju, C.S.; Behera, B.K.; Shukla, S.N.; Kumar, H.; Agarwal, S.K. Anestrus in cattle and buffalo: Indian perspective. Adv. Anim. Vet. Sci. 2014, 2,124–38. [CrossRef]
- Abraham, F. An overview on functional causes of infertility in cows. JFIV Reprod Med Genet. 2017, 5, 203. [CrossRef]
- Rhodes, F.M.; McDougall, S.; Burke, C.R.; Verkerk, G.A.; Macmillan K.L. Invited review: treatment of cows with an extended postpartum anestrous interval. J. Dairy Sci. 2003, 86, 1876–94. [CrossRef]
- Collier, R.J.; Dahl G.E.; VanBaale, M.J. Major advances associated with environmental effects on dairy cattle. J. Dairy Sci. 2006, 89,1244–53. [CrossRef]
- Cerri, R.L.A.; Burnett, T.A.; Madureira, A.M.L.; Silper, B.F.; Denis-Robichaud, J.; LeBlanc, S.; Cooke, R.F.; Vasconcelos, J.L.M. Symposium review: Linking activity-sensor data and physiology to improve dairy cow fertility. J. Dairy Sci. 2021, 104, 1220–31. [CrossRef]
- Lucy, M.C.; McDougall S.; Nation D.P. The use of hormonal treatments to improve the reproductive performance of lactating dairy cows in feedlot or pasture-based management systems. Anim. Reprod. Sci. 2004, 82, 495–512. [CrossRef]
- Melendez, P.; Duchens, M.; Perez, A.; Moraga, L.; Archbald, L. Characterization of estrus detection, conception and pregnancy risk of Holstein cattle from the central area of Chile. Theriogenology 2008, 70, 631–7. [CrossRef]
- Sajapitak, S.; Yawongsa, A.; Arunvipas, P. Estrus synchronization and fixed-time artificial insemination using combined progesterone and prostaglandin F2α together with hCG, or GnRH in cross-bred Indu-Brazil beef cows. J. Kasetsart Vet. 2015, 25, 115–22.
- Deesuk, P.; Wichaiwong, A.; Sirisatien, S.; Aiumlamai, S. The use of hormonal treatment for induction of estrus on conception rate and days open in postpartum dairy cows. KKU. Vet. J. 2020, 8, 23–32.
- Herlihy, M.M.; Crowe, M.A.; Berry, D.P.; Diskin, M.G.; Butler, S.T. Factors associated with fertility outcomes in cows treated with protocols to synchronize estrus and ovulation in seasonal-calving, pasture-based dairy production systems. J. Dairy Sci. 2013, 96, 1485–1498. [CrossRef]
- Kajaysri, J.; Wattananorasate, A. Influence of heat stress on pregnancy rate of beef cattle heifers and cows using ovsynch plus CIDR followed by timed artificial insemination in Thailand. Thai. J. Vet. Med. 2018, 48, 71–77.
- Alnimer, M.; De Rosa, G.; Grasso, F.; Napolitano, F.; Bordi, A. Effect of climate on the response to three oestrous synchronization techniques in lactating dairy cows. Anim. Reprod. Sci. 2002, 71, 157–68. [CrossRef]
- Morini, G.; López-Gatius, F.; Garcia-Ispierto, I.; De Rensis, F. Effect of utilization of single or double prostaglandin administration within an ovsynch fixed-time artificial insemination protocol during summer season in dairy cows. Ann. Anim. Sci. 2019, 19, 725–31. [CrossRef]
- Alnimer, M.A. Comparison of an oestrus synchronisation protocol with oestradiol benzoate and PGF2α and insemination at detected oestrus to a timed insemination protocol (Ovsynch) on reproductive performance of lactating dairy cows. Reprod. Nutr. Dev. 2005, 45, 699–708. [CrossRef]
- Baruselli, P.S.; Ferreira, R.M.; Sá Filho, M.F.; Bó G.A. Review: Using artificial insemination v. natural service in beef herds. Animals. 2018, 12, 45–52. [CrossRef]
- Pereira, M.H.C.; Wiltbank, M.C.; Barbosa, L.F.S.P.; Costa Jr, W.M.; Carvalho, M.A.P.; Vasconcelos, J.L.M. Effect of adding a gonadotropin-releasing-hormone treatment at the beginning and a second prostaglandin F2α treatment at the end of an estradiol-based protocol for timed artificial insemination in lactating dairy cows during cool or hot seasons of the year. J. Dairy Sci. 2015, 98, 947–59. [CrossRef]
- Amin, Y.A.; El-Naga, E.M.A.; Noseer, E.A.; Fouad, S.S.; Ali, R.A. Synchronization with controlled internal drug release (CIDR) and prostaglandin F2α (PGF2α) influences oxidant/antioxidant biomarkers and mineral profile in summer-stressed anoestrous buffalo (Bubalus bubalis). Theriogenology 2019, 134, 34–41. [CrossRef]
- Oliveira, F.A.; Almeida, Í.C.D.; Penitente-Filho, J.M.; Torres, C.A.A. Estradiol and GnRH on ovulation induction for estrus synchronized crossbred cows. Rev. Caatinga. 2020, 33, 815–23. [CrossRef]
- Amin, Y.A.; Said, A. The addition of chitosan to GnRH analog induces ovarian resumption and improves conception rates in buffaloes. Trop. Anim. Sci. J. 2021, 44, 1–9. [CrossRef]
- Bosolasco, D.; Nuñez-Olivera, R.; De Brun, V.; Meikle, A.; Menchaca, A. Estradiol cypionate administered at the end of a progesterone-based protocol for FTAI induces ovulation and improves postovulatory luteal function and uterine environment in anestrous beef cows. Theriogenology 2021, 162, 74–83. [CrossRef]
- Mittal, P.K.; Gottam, G.S.; Gupta, B.; Bilochi, D.R. The effect of climate change on productivity and reproductive and health performance of livestock: A review. J. Entomol. Zool. Stud. 2019, 7, 4–9.
- Sae-tiao, T.; Laodim, T.; Koonawootrittriron, S.; Suwanasopee, T.; Elzo, M.A. Tropical climate change and its effect on milk production of dairy cattle in Thailand. Livest. Res. Rural. Dev. 2019, 31.
- DLD. Information on the number of livestock in Thailand for the year 2020. Available online: https://ict.dld.go.th/webnew/images/stories/report/regislives/2020/T3-1-Milk.pdf (accessed on 17 March 2021).
- Endris, M.; Tumwasorn, S.; Sopannarath, P.; Prasanpanich, S. Genotype by region interaction on milk production traits of Holstein crossbred dairy cows in Thailand. Kasetsart J. (Nat. Sci.). 2013, 47, 228–37.
- Fathoni, A.; Boonkum, W.; Chankitisakul, V.; Duangjinda, M. An Appropriate Genetic Approach for Improving Reproductive Traits in Crossbred Thai–Holstein Cattle under Heat Stress Conditions. Vet. Sci. 2022, 9, 163. [CrossRef]
- Ferguson, J.D.; Galligan, D.T.; Thomsen, N. Principal descriptors of body condition score in Holstein cows. J. Dairy Sci. 1994, 77, 2695–703. [CrossRef]
- Bank of Thailand. Foreign Exchange Rate as of 8 May 2020. Available online: https://www.bot.or.th/english/statistics/_layouts/application/exchangerate/exchangerate.aspx.; 2020 (accessed on 8 May 2020).
- Dohoo, I.; Martin, W.; Stryhn, H. Veterinary Epidemiologic Research, 2nd ed.; VER Inc.: Charlottetown, PEI, Canada, 2010; pp. 395–426.
- Rabiee, A.R.; Lean, I.J.; Stevenson, M.A. Efficacy of Ovsynch program on reproductive performance in dairy cattle: A meta-analysis. J. Dairy Sci. 2005, 88, 2754–70. [CrossRef]
- Burnett, T.A.; Madureira, A.M.; Silper, B.F.; Fernandes, A.C.C.; Cerri, R.L. Integrating an automated activity monitor into an artificial insemination program and the associated risk factors affecting reproductive performance of dairy cows. J. Dairy Sci. 2017, 100, 5005–18. [CrossRef]
- Yeamkong, S.; Koonawootrittriron, S.; Elzo, M.A.; Suwanasopee, T. Effect of experience, education, record keeping, labor and decision making on monthly milk yield and revenue of dairy farms supported by a private organization in Central Thailand. Asian-Australas. J. Anim. Sci. 2010, 23, 814–24. [CrossRef]
- Opsomer, G.; Mijten, P.; Coryn, M.; de Kruif, A. Post-partum anoestrus in dairy cows: A review. Vet. Q. 1996, 18, 68–75. [CrossRef]
- Twagiramungu, H.; Guilbault, L.A.; Dufour, J.J. Synchronization of ovarian follicular waves with a gonadotropin-releasing hormone agonist to increase the precision of estrus in cattle: a review. J. Anim. Sci. 1995, 73, 3141–51. [CrossRef]
- Stevenson, J.S.; Tiffany, S.M.; Lucy, M.C. Use of estradiol cypionate as a substitute for GnRH in protocols for synchronizing ovulation in dairy cattle. J. Dairy Sci. 2004, 87, 3298–305. [CrossRef]
- Sellars, C.B.; Dalton, J.C.; Manzo, R.; Day, J.; Ahmadzadeh, A. Time and incidence of ovulation and conception rates after incorporating estradiol cypionate into a timed artificial insemination protocol. J. Dairy Sci. 2006, 89, 620–6. [CrossRef]
- Kaewlamun, W.; Suwimonteerabutr, J.; Chaimee, T.; Virakul, P.; Techakumphu, M. Low pregnancy rate in dairy cattle after fixed time artificial insemination using Norgestomet+ PGF2α+ eCG program during the hot and humid months in Thailand. Thai. J. Vet. Med. 2008, 38, 53–8.
- Esmaili-Tazangi, M.K.; Mirzaei, A. The effect of body condition loss and milk yield on the efficiency of Ovsynch in cyclic Holstein dairy cows. Rev. Med. Vet. 2015, 166, 345–9.
- Butler, W. Nutrition, negative energy balance and fertility in the postpartum dairy cow. Cattle Pract. 2005, 13, 13–8.
- Changtes, T.; Arunvipas, P.; Rukkwamsuk, T. A Review: Nutritional Factors Influencing Reproductive Performance of Crossbred Holstein Friesian Dairy Cows in Thailand. J. Kasetsart Vet. 2020, 30, 171–88.
- Leroy, J.L.M.R.; Bie, J.; Jordaens, L.; Desmet, K.; Smits, A.; Marei, W.F.A.; Bols, P.E.J.; Hoeck, V.V. Negative energy balance and metabolic stress in relation to oocyte and embryo quality: an update on possible pathways reducing fertility in dairy cows. Anim. Reprod. 2017, 14, 497–506. [CrossRef]
- Madureira, A.M.L.; Polsky, L.B.; Burnett, T.A.; Silper, B.F.; Soriano, S.; Sica, A.F.; Pohler, K. G.; Vasconcelos, J.L.M.; Cerri, R.L.A. Intensity of estrus following an estradiol-progester-one-based ovulation synchronization protocol influences fertility outcomes. J. Dairy Sci. 2019, 102, 3598–608. [CrossRef]
- Madureira, A.M.L.; Silper, B.F.; Burnett, T.A.; Polsky, L.; Cruppe, L.H.; Veira, D.M.; Vasconcelos, J.L.M.; Cerri, R.L.A. Factors affecting expression of estrus measured by activity monitors and conception risk of lactating dairy cows. J. Dairy Sci. 2015, 98, 7003–14. [CrossRef]
- Sinclair, K.D. Declining fertility, insulin resistance and fatty acid metabolism in dairy cows: Developmental consequences for the oocyte and pre-implantation embryo. Acta Sci. Vet. 2010, 38, 545–57.
- Petersson, K.J.; Strandberg, E.; Gustafsson, H.; Berglund, B. Environmental effects on progesterone profile measures of dairy cow fertility. Anim. Reprod. Sci. 2006, 91, 201–14. [CrossRef]
- Sawa, A.; Bogucki, M. Effect of housing system and milk yield on cow fertility. Arch. Anim. Breed. 2011, 54, 249–56. [CrossRef]
- Beaver, A.; Ritter, C.; von Keyserlingk, M.A. The dairy cattle housing dilemma: Natural behavior versus animal care. Vet. Clin. North Am Food Anim. Pract. 2019, 35, 11–27. [CrossRef]
- Endo, N.; Kitamura, T.; Okubo, M.; Tanaka, T. Hair cortisol concentration in pre-and postpartum dairy cows, and its association with body condition, hock health, and reproductive status. Anim. Sci. J. 2019, 90, 924–31. [CrossRef]
- Sharma, A.; Umapathy, G.; Kumar, V.; Phillips, C.J. Hair cortisol in sheltered cows and its association with other welfare indicators. Animals. 2019, 9, 248. [CrossRef]
- Sammad, A.; Umer, S.; Shi, R.; Zhu, H.; Zhao, X.; Wang, Y. Dairy cow reproduction under the influence of heat stress. J. Anim. Physiol. Anim. Nutr. 2019, 104, 978–86. [CrossRef]
- Amin, Y.A.; Youssef, N.A.M.; Mahmoud, A.E.Z.E.; Salah, M.; Khalil, A.M.; Shanab, O.; Hassaneen, A.S.A. Impact of polyherbal formulation oral administration on the estrus response, luteal activity, and oxidative stress in postpartum dairy cows with ovarian subfunction. Vet. World 2022, 1, 360–7. [CrossRef]
- Thammahakin, P.; Yawongsa, A.; Rukkwamsuk, T. Effect of heat stress on reproductive performance of dairy cows under tropical climate: a review. J. Kasetsart Vet. 2020, 30, 111–132.
- Schüller L.K.; Burfeind O.; Heuwieser, W. Impact of heat stress on conception rate of dairy cows in the moderate climate considering different temperature–humidity index thresholds, periods relative to breeding, and heat load indices. Theriogenology 2014, 81, 1050–7. [CrossRef]
- Wolfenson, D.; Lew, B.J.; Thatcher, W.W.; Graber, Y.; Meidan, R. Seasonal and acute heat stress effects on steroid production by dominant follicles in cows. Anim. Reprod. Sci. 1997, 47, 9–19. [CrossRef]
- Wolfenson, D.; Thatcher, W.W.; Badinga, L.; Savio, J.D.; Meidan, R.; Lew, B.J.; Braw-Tal, R.; Berman, A. Effect of heat stress on follicular development during the estrous cycle in lactating dairy cattle. Biol. Reprod. 1995, 52, 1106–13. [CrossRef]
- Villa-Mancera, A.; Méndez-Mendoza, M.; Huerta-Crispín, R.; Vázquez-Flores, F.; Córdova-Izquierdo, A. Effect of climate factors on conception rate of lactating dairy cows in Mexico. Trop. Anim. Health Prod. 2011, 43, 597–601. [CrossRef]
- Roelofs, J.; López-Gatius, F.; Hunter, R.H.F.; Van Eerdenburg, F.J.C.M.; Hanzen, C. When is a cow in estrus? Clinical and practical aspects. Theriogenology 2010, 74, 327–44. [CrossRef]
- Ricci, A.; Li, M.; Fricke, P.M.; Cabrera, V.E. Short communication: Economic impact among 7 reproductive programs for lactating dairy cows, including a sensitivity analysis of the cost of hormonal treatments. J. Dairy Sci. 2020, 103, 5654–61. [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).