Amélie Vanlierde

Potential use of milk mid-infrared spectra to predict individual methane emission of dairy cows

This study investigates the feasibility to predict individual methane (CH(4)) emissions from dairy cows using milk mid-infrared (MIR) spectra. To have a large variability of milk composition, two experiments were conducted on 11 lactating Holstein cows (two primiparous and nine multiparous). The first experiment aimed to induce a large variation in CH(4) emission by feeding two different diets: the first one was mainly composed of fresh grass and sugar beet pulp and the second one of maize silage and hay. The second experiment consisted of grass and corn silage with cracked corn, soybean meal and dried pulp. For each milking period, the milk yields were recorded twice daily and a milk sample of 50 ml was collected from each cow and analyzed by MIR spectrometry. Individual CH(4) emissions were measured daily using the sulfur hexafluoride method during a 7-day period. CH(4) daily emissions ranged from 10.2 to 47.1 g CH(4)/kg of milk. The spectral data were transformed to represent an average daily milk spectrum (AMS), which was related to the recorded daily CH(4) data. By assuming a delay before the production of fermentation products in the rumen and their use to produce milk components, five different calculations were used: AMS at days 0, 0.5, 1, 1.5 and 2 compared with the CH(4) measurement. The equations were built using Partial Least Squares regression. From the calculated R(2)(cv), it appears that the accuracy of CH(4) prediction by MIR changed in function of the milking days. In our experimental conditions, the AMS at day 1.5 compared with the measure of CH(4) emissions gave the best results. The R(2) and s.e. of the cross-validation were equal to 0.79 and 5.14 g of CH(4)/kg of milk. The multiple correlation analysis performed in this study showed the existence of a close relationship between milk fatty acid (FA) profile and CH(4) emission at day 1.5. The lower R(2) (R(2) = 0.76) obtained between FA profile and CH(4) emission compared with the one corresponding to the obtained calibration (R(2)(c) = 0.87) shows the interest to apply directly the developed CH(4) equation instead of the use of correlations between FA and CH(4). In conclusion, our preliminary results suggest the feasibility of direct CH(4) prediction from milk MIR spectra. Additional research has the potential to improve the calibrations even further. This alternative method could be useful to predict the individual CH(4) emissions at farm level or at the regional scale and it also could be used to identify low-CH(4)-emitting cows.

Hot topic: Innovative lactation-stage-dependent prediction of methane emissions from milk mid-infrared spectra

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.

Capitalizing on fine milk composition for breeding and management of dairy cows.

The challenge of managing and breeding dairy cows is permanently adapting to changing production circumstances under socio-economic constraints. If managing and breeding address different timeframes of action, both need relevant phenotypes that allow for precise monitoring of the status of the cows, and their health, behavior, and well-being as well as their environmental impact and the quality of their products (i.e., milk and subsequently dairy products). Milk composition has been identified as an important source of information because it could reflect, at least partially, all these elements. Major conventional milk components such as fat, protein, urea, and lactose contents are routinely predicted by mid-infrared (MIR) spectrometry and have been widely used for these purposes. But, milk composition is much more complex and other nonconventional milk components, potentially predicted by MIR, might be informative. Such new milk-based phenotypes should be considered given that they are cheap, rapidly obtained, usable on a large scale, robust, and reliable. In a first approach, new phenotypes can be predicted from MIR spectra using techniques based on classical prediction equations. This method was used successfully for many novel traits (e.g., fatty acids, lactoferrin, minerals, milk technological properties, citrate) that can be then useful for management and breeding purposes. An innovation was to consider the longitudinal nature of the relationship between the trait of interest and the MIR spectra (e.g., to predict methane from MIR). By avoiding intermediate steps, prediction errors can be minimized when traits of interest (e.g., methane, energy balance, ketosis) are predicted directly from MIR spectra. In a second approach, research is ongoing to detect and exploit patterns in an innovative manner, by comparing observed with expected MIR spectra directly (e.g., pregnancy). All of these traits can then be used to define best practices, adjust feeding and health management, improve animal welfare, improve milk quality, and mitigate environmental impact. Under the condition that MIR data are available on a large scale, phenotypes for these traits will allow genetic and genomic evaluations. Introduction of novel traits into the breeding objectives will need additional research to clarify socio-economic weights and genetic correlations with other traits of interest.

Milk mid-infrared spectra enable prediction of lactation-stage-dependent methane emissions of dairy cattle within routine population-scale milk recording schemes

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.

Changes throughout lactation in phenotypic and genetic correlations between methane emissions and milk fatty acid contents predicted from milk mid-infrared spectra

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.

Genetic parameters of mid-infrared methane predictions and their relationships with milk production traits in Holstein cattle

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.

Short communication: Development of an equation for estimating methane emissions of dairy cows from milk Fourier transform mid-infrared spectra by using reference data obtained exclusively from respiration chambers

Evaluation and mitigation of enteric methane (CH4) emissions from ruminant livestock, in particular from dairy cows, have acquired global importance for sustainable, climate-smart cattle production. Based on CH4 reference measurements obtained with the SF6 tracer technique to determine ruminal CH4 production, a current equation permits evaluation of individual daily CH4 emissions of dairy cows based on milk Fourier transform mid-infrared (FT-MIR) spectra. However, the respiration chamber (RC) technique is considered to be more accurate than SF6 to measure CH4 production from cattle. This study aimed to develop an equation that allows estimating CH4 emissions of lactating cows recorded in an RC from corresponding milk FT-MIR spectra and to challenge its robustness and relevance through validation processes and its application on a milk spectral database. This would permit confirming the conclusions drawn with the existing equation based on SF6 reference measurements regarding the potential to estimate daily CH4 emissions of dairy cows from milk FT-MIR spectra. A total of 584 RC reference CH4 measurements (mean ± standard deviation of 400 ± 72 g of CH4/d) and corresponding standardized milk mid-infrared spectra were obtained from 148 individual lactating cows between 7 and 321 d in milk in 5 European countries (Germany, Switzerland, Denmark, France, and Northern Ireland). The developed equation based on RC measurements showed calibration and cross-validation coefficients of determination of 0.65 and 0.57, respectively, which is lower than those obtained earlier by the equation based on 532 SF6 measurements (0.74 and 0.70, respectively). This means that the RC-based model is unable to explain the variability observed in the corresponding reference data as well as the SF6-based model. The standard errors of calibration and cross-validation were lower for the RC model (43 and 47 g/d vs. 66 and 70 g/d for the SF6 version, respectively), indicating that the model based on RC data was closer to actual values. The root mean squared error (RMSE) of calibration of 42 g/d represents only 10% of the overall daily CH4 production, which is 23 g/d lower than the RMSE for the SF6-based equation. During the external validation step an RMSE of 62 g/d was observed. When the RC equation was applied to a standardized spectral database of milk recordings collected in the Walloon region of Belgium between January 2012 and December 2017 (1,515,137 spectra from 132,658 lactating cows in 1,176 different herds), an average ± standard deviation of 446 ± 51 g of CH4/d was estimated, which is consistent with the range of the values measured using both RC and SF6 techniques. This study confirmed that milk FT-MIR spectra could be used as a potential proxy to estimate daily CH4 emissions from dairy cows provided that the variability to predict is covered by the model.

Estimation of CH4 emissions from milk MIR spectra using respiration chamber as reference technique

AIM Identify signals of fat deposition and adaptation through genome-wide scan of the Barbaresca fat-tail sheep. ANIMALS Barbaresca in an ancient Sicilian fat-tail sheep, highly endangered at present. Of the 35 000 heads of 1980, abour 1 300 are left nowadays in 20 flocks. The breed originated from crosses between Barbary sheep from North Africa and the Pinzirita breed at times of the Arab settling in Sicily (9th century). The breed is reared in a very restricted area in central Sicily on smalland medium-sized farms under a semi-extensive farming system. It is a dual-purpose breed: milk for cheese and meat. Barbaresca is one of the only two fat-tail sheep of Italy. METHODS Genotypic data were obtained with the OvineSNP50K array. Fst values of differentiation for 43072 markers were calculated in pairwise comparisons of Barbaresca with each of 13 Italian thin tail breeds. Fat-tail sheep still represent twenty-five percent of the world sheep population; they are predominant in pastoral, transhumant and low input systems. In Western countries and in high input systems they are generally endangered. Fat-tail sheep preserved genetic variability for functional adaptation. The identification of the genes with a role in the fat-tail phenotype contributes to the understanding of the physiology of fat deposition as well as the mechanisms of adaptation and is essential for maintaining future breeding options. Heritability estimates for the 1st litter size, pregnancy rate and whelping success were low (0.05-0.14)  Grading size and quality had moderate heritability estimates 0.27 and 0.21, respectively  Genetic correlations between animal grading size and fertility traits were unfavourable (from -0.15 to -0.53)  Grading quality and guard hair coverage had antagonistic relationships with all the studied fertility traits (from -0.21 to -0.54) Genetic parameters of fertility and grading traits in Finnish blue foxTrabajo presentado al: 68th Annual Meeting of the European Federation of Animal Science (EAAP). (Tallin, Estonia. 28 agosto - 2 septiembre).Trabajo presentado al: 68th Annual Meeting of the European Federation of Animal Science (EAAP). (Tallin, Estonia. 28 agosto - 2 septiembre).

Relationships between milk mid-IR predicted gastro-enteric methane production and the technical and financial performance of commercial dairy herds

Considering economic and environmental issues is important in ensuring the sustainability of dairy farms. The objective of this study was to investigate univariate relationships between lactating dairy cow gastro-enteric methane (CH4) production predicted from milk mid-IR (MIR) spectra and technico-economic variables by the use of large scale and on-farm data. A total of 525 697 individual CH4 predictions from milk MIR spectra (MIR-CH4 (g/day)) of milk samples collected on 206 farms during the Walloon milk recording scheme were used to create a MIR-CH4 prediction for each herd and year (HYMIR-CH4). These predictions were merged with dairy herd accounting data. This allowed a simultaneous study of HYMIR-CH4 and 42 technical and economic variables for 1024 herd and year records from 2007 to 2014. Pearson correlation coefficients (r) were used to assess significant relationships (P<0.05). Low HYMIR-CH4 was significantly associated with, amongst others, lower fat and protein corrected milk (FPCM) yield (r=0.18), lower milk fat and protein content (r=0.38 and 0.33, respectively), lower quantity of milk produced from forages (r=0.12) and suboptimal reproduction and health performance (e.g. longer calving interval (r=-0.21) and higher culling rate (r=-0.15)). Concerning economic results, low HYMIR-CH4 was significantly associated with lower gross margin per cow (r=0.19) and per litre FPCM (r=0.09). To conclude, this study suggested that low lactating dairy cow gastro-enteric CH4 production tended to be associated with more extensive or suboptimal management practices, which could lead to lower profitability. The observed low correlations suggest complex interactions between variables due to the use of on-farm data with large variability in technical and management practices.

Variation of greenhouse gas emissions and identification of their drivers during the fattening of Belgian Blue White bulls based on a LCA approach

Greenhouse gas emission intensity (GHGI; kilograms carbon dioxide equivalents/kilograms liveweight gain) have to be reduced so as to limit the impact of human activities on global warming while furnishing food to human. In this respect, performances of 654 Belgian Blue double-muscled bulls (BBdm) during their fattening phase were recorded. On this basis, their greenhouse gas emissions were modelled to estimate variation in GHGI and investigate mitigation options at that level. The relevance of theses option is discussed, taking into account the whole life and production system scales. Large variations (mean (s.d.)) were observed (from 7.2 (0.4) to 10.0 (0.7) kg carbon dioxide equivalents/kg liveweight gain) for, respectively, the 1st- and 4th-quantile groups defined for GHGI. Early culling, low liveweight and age at start of the fattening phase of the bulls would lead to a reduction of GHGI. Nevertheless, more than 32% of the variation remained unexplained. However, decision leading to reduction of GHG intensity at this stage of the life may be compensated in the early stage of BBdm. Attention is drawn on the necessity to encompass the whole life of BBdm for investigating mitigation options and on the sensitivity of the results on models and methodological choices.