Abstract: The data set for this study comprised 1,488,474 test-day records for milk, fat, and protein yields and fat and protein percentages from 191,012 first-, second-, and third-parity Holstein cows from 484 farms. Data were collected from 2001 through 2007 and merged with meteorological data from 35 weather stations. A linear model (M1) was used to estimate the effects of the temperature-humidity index (THI) on production traits. Least squares means from M1 were used to detect the THI thresholds for milk production in all parities by using a 2-phase linear regression procedure (M2). A multiple-trait repeatability test-model (M3) was used to estimate variance components for all traits and a dummy regression variable (t) was defined to estimate the production decline caused by heat stress. Additionally, the estimated variance components and M3 were used to estimate traditional and heat-tolerance breeding values (estimated breeding values, EBV) for milk yield and protein percentages at parity 1. An analysis of data (M2) indicated that the daily THI at which milk production started to decline for the 3 parities and traits ranged from 65 to 76. These THI values can be achieved with different temperature/humidity combinations with a range of temperatures from 21 to 36°C and relative humidity values from 5 to 95%. The highest negative effect of THI was observed 4 d before test day over the 3 parities for all traits. The negative effect of THI on production traits indicates that first-parity cows are less sensitive to heat stress than multiparous cows. Over the parities, the general additive genetic variance decreased for protein content and increased for milk yield and fat and protein yield. Additive genetic variance for heat tolerance showed an increase from the first to third parity for milk, protein, and fat yield, and for protein percentage. Genetic correlations between general and heat stress effects were all unfavorable (from -0.24 to -0.56). Three EBV per trait were calculated for each cow and bull (traditional EBV, traditional EBV estimated with the inclusion of THI covariate effect, and heat tolerance EBV) and the rankings of EBV for 283 bulls born after 1985 with at least 50 daughters were compared. When THI was included in the model, the ranking for 17 and 32 bulls changed for milk yield and protein percentage, respectively. The heat tolerance genetic component is not negligible, suggesting that heat tolerance selection should be included in the selection objectives.

Abstract: A retrospective study on seasonal variations in the characteristics of cow’s milk and temperature-humidity index (THI) relationship was conducted on bulk milk data collected from 2003 to 2009. The THI relationship study was carried out on 508 613 bulk milk data items recorded in 3328 dairy farms form the Lombardy region, Italy. Temperature and relative humidity data from 40 weather stations were used to calculate THI. Milk characteristics data referred to somatic cell count (SCC), total bacterial count (TBC), fat percentage (FA%) and protein percentage (PR%). Annual, seasonal and monthly variations in milk composition were evaluated on 656 064 data items recorded in 3727 dairy farms. The model highlighted a significant association between the year, season and month, and the parameters analysed (SCC, TBC, FA%, PR%). The summer season emerged as the most critical season. Of the summer months, July presented the most critical conditions for TBC, FA% and PR%, (52 054 ± 183 655, 3.73% ± 0.35% and 3.30% ± 0.15%, respectively), and August presented higher values of SCC (369 503 ± 228 377). Each milk record was linked to THI data calculated at the nearest weather station. The analysis demonstrated a positive correlation between THI and SCC and TBC, and indicated a significant change in the slope at 57.3 and 72.8 maximum THI, respectively. The model demonstrated a negative correlation between THI and FA% and PR% and provided breakpoints in the pattern at 50.2 and 65.2 maximum THI, respectively. The results of this study indicate the presence of critical climatic thresholds for bulk tank milk composition in dairy cows. Such indications could facilitate the adoption of heat management strategies, which may ensure the health and production of dairy cows and limit related economic losses.