RESUMEN
Heat stress profoundly affects the reproductive success of buffaloes, which are vital for the dairy industry due to their unique anatomical and physiological characteristics, necessitating careful evaluation under such conditions. Hence, this guided our search for quantifying heat stress' impact on Mehsana buffaloes using the best THI model and evaluating sires' performance. Fertility records (days open converted to daughter pregnancy rate) were collected in the span of over 24 years, w.e.f. 1989 to 2012. Finally, 3070 records of first lactation cows, daughters of 117 sires from DURDA, Gujarat, India, were used in the analysis. Meteorological data were retrieved from IMD, Pune, to understand the relationship between daughter pregnancy rate (DPR) and heat stress indicators. Several heat stress models were compared based on R2, adjusted R2, AIC, and BIC values, and the impact of heat stress was quantified. The year was classified into different heat stress zones, viz., Non heat stress zone (NHSZ), Heat stress zone (HSZ), and critical heat stress zone (CHSZ), drawing from the findings of DPR and THI. The THI 4th model based on dry and wet bulb temperature was identified as the best-fit model, and DPR significantly changed (P < 0.01) by 1.14% per unit change in THI value based on the month of calving. The average EBVs of the sires for DPR were found to be 20.78% (NHSZ), 38.09% (HSZ), and 39.08% (CHSZ) using BLUP-SM and 20.78% (NHSZ), 37.30% (HSZ), and 38.87% (HSZ) using BLUP-AM. Subsequently, the optimum sire for each of the zones was prioritized. It is noteworthy that bulls that performed better in NHSZ did not perform as well in HSZ and CHSZ, and vice versa. This supports the possibility of evaluating bulls independently in each heat stress zone.
