Genetic parameters that considered tolerance for heat stress were estimated for production, udder health, and milk composition traits. Data included 202,733 test-day records for milk, fat, and protein yields, fat and protein percentages, somatic cell score (SCS), 10 individual milk fatty acids (FA) predicted by mid-infrared spectrometry, and 7 FA groups. Data were from 34,468 first-lactation Holstein cows in 862 herds in the Walloon region of Belgium and were collected between 2007 and 2010. Test-day records were merged with daily temperature-humidity index (THI) values based on meteorological records from public weather stations. The maximum distance between each farm and its corresponding weather station was 21km. Linear reaction norm models were used to estimate the intercept and slope responses of 23 traits to increasing THI values. Most yield and FA traits had phenotypic and genetic declines as THI increased, whereas SCS, C18:0, C18:1 cis-9, and 4 FA groups (unsaturated FA, monounsaturated FA, polyunsaturated FA, and long-chain FA) increased with THI. Moreover, the latter traits had the largest slope-to-intercept genetic variance ratios, which indicate that they are more affected by heat stress at high THI levels. Estimates of genetic correlations within trait between cold and hot environments were generally high (>0.80). However, lower estimates (<=0.67) were found for SCS, fat yield, and C18:1 cis-9, indicating that animals with the highest genetic merit for those traits in cold environments do not necessarily have the highest genetic merit for the same traits in hot environments. Among all traits, C18:1 cis-9 was the most sensitive to heat stress. As this trait is known to reflect body reserve mobilization, using its variations under hot conditions could be a very affordable milk biomarker of heat stress for dairy cattle expressing the equilibrium between intake and mobilization under warm conditions.

The main goal of this study was to develop, apply, and validate a new method to predict an indicator for CH4 eructed by dairy cows using milk mid-infrared (MIR) spectra. A novel feature of this model was the consideration of lactation stage to reflect changes in the metabolic status of the cow. A total of 446 daily CH4 measurements were obtained using the SF6 method on 142 Jersey, Holstein, and Holstein-Jersey cows. The corresponding milk samples were collected during these CH4 measurements and were analyzed using MIR spectroscopy. A first derivative was applied to the milk MIR spectra. To validate the novel calibration equation incorporating days in milk (DIM), 2 calibration processes were developed: the first was based only on CH4 measurements and milk MIR spectra (independent of lactation stage; ILS); the second included milk MIR spectra and DIM information (dependent on lactation stage; DLS) by using linear and quadratic modified Legendre polynomials. The coefficients of determination of ILS and DLS equations were 0.77 and 0.75, respectively, with standard error of calibration of 63g/d of CH4 for both calibration equations. These equations were applied to 1,674,763 milk MIR spectra from Holstein cows in the first 3 parities and between 5 and 365 DIM. The average CH4 indicators were 428, 444, and 448g/d by ILS and 444, 467, and 471g/d by DLS for cows in first, second, and third lactation, respectively. Behavior of the DLS indicator throughout the lactations was in agreement with the literature with values increasing between 0 and 100 DIM and decreasing thereafter. Conversely, the ILS indicator of CH4 emission decreased at the beginning of the lactation and increased until the end of the lactation, which differs from the literature. Therefore, the DLS indicator seems to better reflect biological processes that drive CH4 emissions than the ILS indicator. The ILS and DLS equations were applied to an independent data set, which included 59 respiration chamber measurements of CH4 obtained from animals of a different breed across a different production system. Results indicated that the DLS equation was much more robust than the ILS equation allowing development of indicators of CH4 emissions by dairy cows. Integration of DIM information into the prediction equation was found to be a good strategy to obtain biologically meaningful CH4 values from lactating cows by accounting for biological changes that occur throughout the lactation.

Renewed interest in heat stress effects on livestock productivity derives from climate change, which is expected to increase temperatures and the frequency of extreme weather events. This study aimed at evaluating the effect of temperature and humidity on milk production in highly selected dairy cattle populations across 3 European regions differing in climate and production systems to detect differences and similarities that can be used to optimize heat stress (HS) effect modeling. Milk, fat, and protein test day data from official milk recording for 1999 to 2010 in 4 Holstein populations located in the Walloon Region of Belgium (BEL), Luxembourg (LUX), Slovenia (SLO), and southern Spain (SPA) were merged with temperature and humidity data provided by the state meteorological agencies. After merging, the number of test day records/cows per trait ranged from 686,726/49,655 in SLO to 1,982,047/136,746 in BEL. Values for the daily average and maximum temperature-humidity index (THIavg and THImax) ranges for THIavg/THImax were largest in SLO (22-74/28-84) and shortest in SPA (39-76/46-83). Change point techniques were used to determine comfort thresholds, which differed across traits and climatic regions. Milk yield showed an inverted U-shaped pattern of response across the THI scale with a HS threshold around 73 THImax units. For fat and protein, thresholds were lower than for milk yield and were shifted around 6 THI units toward larger values in SPA compared with the other countries. Fat showed lower HS thresholds than protein traits in all countries. The traditional broken line model was compared with quadratic and cubic fits of the pattern of response in production to increasing heat loads. A cubic polynomial model allowing for individual variation in patterns of response and THIavg as heat load measure showed the best statistical features. Higher/lower producing animals showed less/more persistent production (quantity and quality) across the THI scale. The estimated correlations between comfort and THIavg values of 70 (which represents the upper end of the THIavg scale in BEL-LUX) were lower for BEL-LUX (0.70-0.80) than for SPA (0.83-0.85). Overall, animals producing in the more temperate climates and semi-extensive grazing systems of BEL and LUX showed HS at lower heat loads and more re-ranking across the THI scale than animals producing in the warmer climate and intensive indoor system of SPA.