Search

The aim of this study was to estimate phenotypic and genetic correlations between methane production (Mp) and milk fatty acid contents of first-parity Walloon Holstein cows throughout lactation. Calibration equations predicting daily Mp (g/d) and milk fatty acid contents (g/100 dL of milk) were applied on milk mid-infrared spectra related to Walloon milk recording. A total of 241,236 predictions of Mp and milk fatty acids were used. These data were collected between 5 and 305 d in milk in 33,555 first-parity Holstein cows from 626 herds. Pedigree data included 109,975 animals. Bivariate (i.e., Mp and a fatty acid trait) random regression test-day models were developed to estimate phenotypic and genetic parameters of Mp and milk fatty acids. Individual short-chain fatty acids (SCFA) and groups of saturated fatty acids, SCFA, and medium-chain fatty acids showed positive phenotypic and genetic correlations with Mp (from 0.10 to 0.16 and from 0.23 to 0.30 for phenotypic and genetic correlations, respectively), whereas individual long-chain fatty acids (LCFA), and groups of LCFA, monounsaturated fatty acids, and unsaturated fatty acids showed null to positive phenotypic and genetic correlations with Mp (from −0.03 to 0.13 and from −0.02 to 0.32 for phenotypic and genetic correlations, respectively). However, these correlations changed throughout lactation. First, de novo individual and group fatty acids (i.e., C4:0, C6:0, C8:0, C10:0, C12:0, C14:0, SCFA group) showed low phenotypic or genetic correlations (or both) in early lactation and higher at the end of lactation. In contrast, phenotypic and genetic correlations between Mp and C16:0, which could be de novo synthetized or derived from blood lipids, were more stable during lactation. This fatty acid is the most abundant fatty acid of the saturated fatty acid and medium-chain fatty acid groups of which correlations with Mp showed the same pattern across lactation. Phenotypic and genetic correlations between Mp and C17:0 and C18:0 were low in early lactation and increased afterward. Phenotypic and genetic correlations between Mp and C18:1 cis-9 originating from the blood lipids were negative in early lactation and increased afterward to become null from 18 wk until the end of lactation. Correlations between Mp and groups of LCFA, monounsaturated fatty acids, and unsaturated fatty acids showed a similar or intermediate pattern across lactation compared with fatty acids that compose them. Finally, these results indicate that correlations between Mp and milk fatty acids vary following lactation stage of the cow, a fact still often ignored when trying to predict Mp from milk fatty acid profile.

Body weight (BW) of dairy cows can be estimated using linear conformation traits (calculated BW; CBW), which are generally recorded only once during a lactation. However, predicted BW (PBW) throughout the lactation would be useful, e.g., at milk-recording dates allowing feed-intake prediction for advisory purposes. Therefore, a 2-step approach was developed to obtain PBW for each milk-recording date. In the first step, a random-regression test-day model was used with CBW as observations to predict PBW. The second step consisted in changing means and (co)variances of prior distributions for the additive genetic random effects of the test-day model by using priors derived from results of the first step to predict again PBW. A total of 25,061 CBW from 24,919 primiparous Holstein cows were computed using equations from literature. Using CBW as observations, PBW was then predicted over the whole lactation for 232,436 dates corresponding to 207,375milk-recording dates and 25,061 classification dates. Results showed that using both steps (the 2-step approach) provided more accurate predictions than using only the first step (the one-step approach). Based on the results of this preliminary study, BW of dairy cows could be predicted throughout the lactation using this procedure. These predictions could be useful in milk-recording systems to compute traits of interest (e.g., feed-intake prediction). The developed novel method is also flexible because actual direct measurements of BW can also be used together with CBW, the prediction model being able to accommodate different levels of accuracies of used BW phenotypes.

Many countries have pledged to reduce greenhouse gases. In this context, the dairy sector is one of the identified sectors to adapt production circumstances to address socio-environmental constraints due to its large carbon footprint related to CH4 emission. This study aimed mainly to estimate (1) the genetic parameters of 2 milk mid-infrared-based CH4 proxies [predicted daily CH4 emission (PME, g/d), and log-transformed predicted CH4 intensity (LMI)] and (2) their genetic correlations with milk production traits [milk (MY), fat (FY), and protein (PY) yields] from first- and second-parity Holstein cows. A total of 336,126 and 231,400 mid-infrared CH4 phenotypes were collected from 56,957 and 34,992 first- and second-parity cows, respectively. The PME increased from the first to the second lactation (433 vs. 453 g/d) and the LMI decreased (2.93 vs. 2.86). We used 20 bivariate random regression test-day models to estimate the variance components. Moderate heritability values were observed for both CH4 traits, and those values decreased slightly from the first to the second lactation (0.25 ± 0.01 and 0.22 ± 0.01 for PME; 0.18 ± 0.01 and 0.17 ± 0.02 for LMI). Lactation phenotypic and genetic correlations were negative between PME and MY in both first and second lactations (-0.07 vs. -0.07 and -0.19 vs. -0.24, respectively). More close scrutiny revealed that relative increase of PME was lower with high MY levels even reverting to decrease, and therefore explaining the negative correlations, indicating that higher producing cows could be a mitigation option for CH4 emission. The PME phenotypic correlations were almost equal to 0 with FY and PY for both lactations. However, the genetic correlations between PME and FY were slightly positive (0.11 and 0.12), whereas with PY the correlations were slightly negative (-0.05 and -0.04). Both phenotypic and genetic correlations between LMI and MY or PY or FY were always relatively highly negative (from -0.21 to -0.88). As the genetic correlations between PME and LMI were strong (0.71 and 0.72 in first and second lactation), the selection of one trait would also strongly influence the other trait. However, in animal breeding context, PME, as a direct quantity CH4 proxy, would be preferred to LMI, which is a ratio trait of PME with a trait already in the index. The range of PME sire estimated breeding values were 22.1 and 29.41 kg per lactation in first and second parity, respectively. Further studies must be conducted to evaluate the effect of the introduction of PME in a selection index on the other traits already included in this index, such as, for instance, fertility or longevity.

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. © 2016 American Dairy Science Association.

Methane (CH4) emission is an important environmental trait in dairy cows. Breeding aiming to mitigate CH4 emissions require the estimation of genetic correlations with other economically important traits and the prediction of their selection response. In this study, test-day CH4 emissions were predicted from milk mid-infrared spectra of Holstein cows. Predicted CH4 emissions (PME) and log-transformed CH4 intensity (LMI) computed as the natural logarithm of PME divided by milk yield (MY). Genetic correlations of PME and LMI with traits used currently were approximated from correlations between estimated breeding values of sires. Values were for PME with MY 0.06, fat yield (FY) 0.09, protein yield (PY) 0.13, fertility 0.17; body condition score (BCS) –0.02; udder health (UDH) 0.22; and longevity 0.22. As expected by its definition, values were negative for LMI with production traits (MY –0.61; FY –0.15 and PY –0.40) and positive with fertility (0.36); BCS (0.20); UDH (0.08) and longevity (0.06). The genetic correlations of 33 type traits with PME ranged from –0.12 to 0.25 and for LMI ranged from –0.22 to 0.18. Without selecting PME and LMI (status quo) the relative genetic change through correlated responses of other traits were in PME by 2% and in LMI by –15%, but only due to the correlated response to MY. Results showed for PME that direct selection of this environmental trait would reduce milk carbon foot print but would also affect negatively fertility. Therefore, more profound changes in current indexes will be required than simply adding environmental traits as these traits also affect the expected progress of other traits.

Methane (CH4) emission is an important environmental trait in dairy cows. Breeding aiming to mitigate CH4 emissions require the estimation of genetic correlations with other economically important traits and the prediction of their selection response. In this study, test-day CH4 emissions were predicted from milk mid-infrared spectra of Holstein cows. Predicted CH4 emissions (PME) and log-transformed CH4 intensity (LMI) computed as the natural logarithm of PME divided by milk yield (MY). Genetic correlations of PME and LMI with traits used currently were approximated from correlations between estimated breeding values of sires. Values were for PME with MY 0.06, fat yield (FY) 0.09, protein yield (PY) 0.13, fertility 0.17; body condition score (BCS) –0.02; udder health (UDH) 0.22; and longevity 0.22. As expected by its definition, values were negative for LMI with production traits (MY –0.61; FY –0.15 and PY –0.40) and positive with fertility (0.36); BCS (0.20); UDH (0.08) and longevity (0.06). The genetic correlations of 33 type traits with PME ranged from –0.12 to 0.25 and for LMI ranged from –0.22 to 0.18. Without selecting PME and LMI (status quo) the relative genetic change through correlated responses of other traits were in PME by 2% and in LMI by –15%, but only due to the correlated response to MY. Results showed for PME that direct selection of this environmental trait would reduce milk carbon foot print but would also affect negatively fertility. Therefore, more profound changes in current indexes will be required than simply adding environmental traits as these traits also affect the expected progress of other traits.

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 (

Mitigating the proportion of energy intake lost as methane could improve the sustainability and profitability of dairy production. As widespread measurement of methane emissions is precluded by current in vivo methods, the development of an easily measured proxy is desirable. An equation has been developed to predict methane from the mid-infrared (MIR) spectra of milk within routine milk-recording programs. The main goals of this study were to improve the prediction equation for methane emissions from milk MIR spectra and to illustrate its already available usefulness as a high throughput phenotypic screening tool. A total of 532 methane measurements considered as reference data (430 ± 129 g of methane/day) linked with milk MIR spectra were obtained from 165 cows using the SF6 technique. A first derivative was applied to the MIR spectra. Constant (P0), linear (P1) and quadratic (P2) modified Legendre polynomials were computed from each cows stage of lactation (days in milk), at the day of SF6 methane measurement. The calibration model was developed using a modified partial least-squares regression on first derivative MIR data points × P0, first derivative MIR data points × P1, and first derivative MIR data points × P2 as variables. The MIR-predicted methane emissions (g/day) showed a calibration coefficient of determination of 0.74, a cross-validation coefficient of determination of 0.70 and a standard error of calibration of 66 g/day. When applied to milk MIR spectra recorded in the Walloon Region of Belgium (≈2 000 000 records), this equation was useful to study lactational, annual, seasonal, and regional methane emissions. We conclude that milk MIR spectra has potential to be used to conduct high throughput screening of lactating dairy cattle for methane emissions. The data generated enable monitoring of methane emissions and production characteristics across and within herds. Milk MIR spectra could now be used for widespread screening of dairy herds in order to develop management and genetic selection tools to reduce methane emissions.