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: All mechanistic rumen models cover the main drivers of variation in rumen function, which are feed intake, the differences between feedstuffs and feeds in their intrinsic rumen degradation characteristics, and fractional outflow rate of fluid and particulate matter. Dynamic modeling approaches are best suited to the prediction of more nuanced responses in rumen metabolism, and represent the dynamics of the interactions between substrates and micro-organisms and inter-microbial interactions. The concepts of dynamics are discussed for the case of rumen starch digestion as influenced by starch intake rate and frequency of feed intake, and for the case of fermentation of fiber in the large intestine. Adding representations of new functional classes of micro-organisms (i.e., with new characteristics from the perspective of whole rumen function) in rumen models only delivers new insights if complemented by the dynamics of their interactions with other functional classes. Rumen fermentation conditions have to be represented due to their profound impact on the dynamics of substrate degradation and microbial metabolism. Although the importance of rumen pH is generally acknowledged, more emphasis is needed on predicting its variation as well as variation in the processes that underlie rumen fluid dynamics. The rumen wall has an important role in adapting to rapid changes in the rumen environment, clearing of volatile fatty acids (VFA), and maintaining rumen pH within limits. Dynamics of rumen wall epithelia and their role in VFA absorption needs to be better represented in models that aim to predict rumen responses across nutritional or physiological states. For a detailed prediction of rumen N balance there is merit in a dynamic modeling approach compared to the static approaches adopted in current protein evaluation systems. Improvement is needed on previous attempts to predict rumen VFA profiles, and this should be pursued by introducing factors that relate more to microbial metabolism. For rumen model construction, data on rumen microbiomes are preferably coupled with knowledge consolidated in rumen models instead of relying on correlations with rather general aspects of treatment or animal. This helps to prevent the disregard of basic principles and underlying mechanisms of whole rumen function.