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.

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.

Bayesian single-step genomic evaluations combining local and foreign information in Walloon Holsteins.

Most dairy cattle populations found in different countries around the world are small to medium sized and use many artificial insemination bulls imported from different foreign countries. The Walloon population in the southern part of Belgium is a good example for such a small-scale population. Wallonia has also a very active community of Holstein breeders requesting high level genetic evaluation services. Single-step Genomic BLUP (ssGBLUP) methods allow the simultaneous use of genomic, pedigree and phenotypic information and could reduce potential biases in the estimation of genomically enhanced breeding values (GEBV). Therefore, in the context of implementing a Walloon genomic evaluation system for Holsteins, it was considered as the best option. However, in contrast to multi-step genomic predictions, natively ssGBLUP will only use local phenotypic information and is unable to use directly important other sources of information coming from abroad, for example Multiple Across Country Evaluation (MACE) results as provided by the Interbull Center (Uppsala, Sweden). Therefore, we developed and implemented single-step Genomic Bayesian Prediction (ssGBayes), as an alternative method for the Walloon genomic evaluations. The ssGBayes method approximated the correct system of equations directly using estimated breeding values (EBV) and associated reliabilities (REL) without any explicit deregression step. In the Walloon genomic evaluation, local information refers to Walloon EBV and REL and foreign information refers to MACE EBV and associated REL. Combining simultaneously all available genotypes, pedigree, local and foreign information in an evaluation can be achieved but adding contributions to left-hand and right-hand sides subtracting double-counted contributions. Correct propagation of external information avoiding double counting of contributions due to relationships and due to records can be achieved. This ssGBayes method computed more accurate predictions for all types of animals. For example, for genotyped animals with low Walloon REL (<0.25) without MACE results but sired by genotyped bulls with MACE results, the average increase of REL for the studied traits was 0.38 points of which 0.08 points could be traced to the inclusion of MACE information. For other categories of genotyped animals, the contribution by MACE information was also high. The Walloon genomic evaluation system passed for the first time the Interbull GEBV tests for several traits in July 2013. Recent experiences reported here refer to its use in April 2016 for the routine genomic evaluations of milk production, udder health and type traits. Results showed that the proposed methodology should also be of interest for other, similar, populations.

Genetic evaluation for birth and conformation traits in dual-purpose Belgian Blue cattle using a mixed inheritance model.

The segregation of the causal mutation () in the muscular hypertrophy gene in dual-purpose Belgian Blue (dpBB) cattle is considered to result in greater calving difficulty (dystocia). Establishing adapted genetic evaluations might overcome this situation through efficient selection. However, the heterogeneity of dpBB populations at the locus implies separating the major gene and other polygenic effects in complex modeling. The use of mixed inheritance models may be an interesting option because they simultaneously assume both influences. A genetic evaluation in dpBB based on a mixed inheritance model was developed for birth and conformation traits: gestation length (GL), calving difficulty (CD), birth weight (BiW), and body conformation score (BC). A total of 27,362 animals having records were used for analyses. The total number of animals in the pedigree used to build the numerator relationship matrix was 62,617. Genotypes at the locus were available for 2,671 animals. Missing records at this locus were replaced with genotype probabilities. A total of 13,221 (48.3%) were registered as dpBB, 1,287 (4.7%) as beef Belgian Blue, and 12,854 (47.0%) were unknown. From those 13,221 dpBB animals, 650, 849, and 534 had double or single copies or no copy, respectively, of the causal mutation () in the muscular hypertrophy gene, whereas 11,188 had missing genotypes. This heterogeneity at the locus may be the reason for high variability in the studied traits, that is, high heritability estimates of 0.33, 0.30, 0.38, and 0.43 for GL, CD, BiW, and BC, respectively. In general, additive ( < 0.05) and dominance ( < 0.001) allele substitution for calves and dams had significant impact for all traits. The moderate coefficient of genetic variation (27.80%) and high direct heritability (0.28) for CD suggested genetic variability in dpBB and possible genetic improvement through selection. This variability has allowed dpBB breeders to successfully apply mass selection in the past. Genetic trend means from 1988 to 2016 showed that sire selection for CD within genotype was progressively applied by breeders. The selection intensity was more important for CD in double-muscled lines than in segregated lines. Our study illustrated the possible confusion caused by the use of major genes in selection and the importance of fitting appropriate models such as mixed inheritance models that combine polygenic and gene content information.

Assessing the effect of pregnancy stage on milk composition of dairy cows using mid-infrared spectra.

Changes in milk production traits (i.e., milk yield, fat, and protein contents) with the pregnancy stage are well documented. To our knowledge, the effect of pregnancy on the detailed milk composition has not been studied so far. The mid-infrared (MIR) spectrum reflects the detailed composition of a milk sample and is obtained by a nonexhaustive and widely used method for milk analysis. Therefore, this study aimed to investigate the effect of pregnancy on milk MIR spectrum in addition to milk production traits (milk yield, fat, and protein contents). A model including regression on the number of days pregnant was applied on milk production traits (milk yield, fat, and protein contents) and on 212 spectral points from the MIR spectra of 9,757 primiparous Holstein cows from Walloon herds. Effects of pregnancy stage were expressed on a relative scale (effect divided by the squared root of the phenotypic variance); this allowed comparisons between effects on milk traits and on 212 spectral points. Effect of pregnancy stage on production traits were in line with previous studies indicating that the model accounted well for the pregnancy effect. Trends of the relative effect of the pregnancy stage on the 212 spectral points were consistent with known and observed effect on milk traits. The highest effect of the pregnancy was observed in the MIR spectral region from 968 to 1,577 cm-1. For some specific wavenumbers, the effect was higher than for fat and protein contents in the beginning of the pregnancy (from 30 to 90 or 120 d pregnant). In conclusion, the effect of early pregnancy can be observed in the detailed milk composition through the analysis of the MIR spectrum of bovine milk. Further analyses are warranted to explore deeply the use of MIR spectra of bovine milk for breeding and management of dairy cow pregnancy.

Development of Fourier transform mid-infrared calibrations to predict acetone, ß-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network.

To manage negative energy balance and ketosis in dairy farms, rapid and cost-effective detection is needed. Among the milk biomarkers that could be useful for this purpose, acetone and ß-hydroxybutyrate (BHB) have been proved as molecules of interest regarding ketosis and citrate was recently identified as an early indicator of negative energy balance. Because Fourier transform mid-infrared spectrometry can provide rapid and cost-effective predictions of milk composition, the objective of this study was to evaluate the ability of this technology to predict these biomarkers in milk. Milk samples were collected in commercial and experimental farms in Luxembourg, France, and Germany. Acetone, BHB, and citrate contents were determined by flow injection analysis. Milk mid-infrared spectra were recorded and standardized for all samples. After edits, a total of 548 samples were used in the calibration and validation data sets for acetone, 558 for BHB, and 506 for citrate. Acetone content ranged from 0.020 to 3.355mmol/L with an average of 0.103mmol/L; BHB content ranged from 0.045 to 1.596mmol/L with an average of 0.215mmol/L; and citrate content ranged from 3.88 to 16.12mmol/L with an average of 9.04mmol/L. Acetone and BHB contents were log-transformed and a part of the samples with low values was randomly excluded to approach a normal distribution. The 3 edited data sets were then randomly divided into a calibration data set (3/4 of the samples) and a validation data set (1/4 of the samples). Prediction equations were developed using partial least square regression. The coefficient of determination (R(2)) of cross-validation was 0.73 for acetone, 0.71 for BHB, and 0.90 for citrate with root mean square error of 0.248, 0.109, and 0.70mmol/L, respectively. Finally, the external validation was performed and R(2) obtained were 0.67 for acetone, 0.63 for BHB, and 0.86 for citrate, with respective root mean square error of validation of 0.196, 0.083, and 0.76mmol/L. Although the practical usefulness of the equations developed should be further verified with other field data, results from this study demonstrated the potential of Fourier transform mid-infrared spectrometry to predict citrate content with good accuracy and to supply indicative contents of BHB and acetone in milk, thereby providing rapid and cost-effective tools to manage ketosis and negative energy balance in dairy farms.

Integration of external estimated breeding values and associated reliabilities using correlations among traits and effects.

Based on a Bayesian view of linear mixed models, several studies showed the possibilities to integrate estimated breeding values (EBV) and associated reliabilities (REL) provided by genetic evaluations performed outside a given evaluation system into this genetic evaluation. Hereafter, the term "internal" refers to this given genetic evaluation system, and the term "external" refers to all other genetic evaluations performed outside the internal evaluation system. Bayesian approaches integrate external information (i.e., external EBV and associated REL) by altering both the mean and (co)variance of the prior distributions of the additive genetic effects based on the knowledge of this external information. Extensions of the Bayesian approaches to multivariate settings are interesting because external information expressed on other scales, measurement units, or trait definitions, or associated with different heritabilities and genetic parameters than the internal traits, could be integrated into a multivariate genetic evaluation without the need to convert external information to the internal traits. Therefore, the aim of this study was to test the integration of external EBV and associated REL, expressed on a 305-d basis and genetically correlated with a trait of interest, into a multivariate genetic evaluation using a random regression test-day model for the trait of interest. The approach we used was a multivariate Bayesian approach. Results showed that the integration of external information led to a genetic evaluation for the trait of interest for, at least, animals associated with external information, as accurate as a bivariate evaluation including all available phenotypic information. In conclusion, the multivariate Bayesian approaches have the potential to integrate external information correlated with the internal phenotypic traits, and potentially to the different random regressions, into a multivariate genetic evaluation. This allows the use of different scales, heritabilities, variance components, measurement units, or trait definitions for external and internal traits. However, one possible issue for implementing multivariate Bayesian approaches could be the availability or estimation of genetic correlations between external and internal traits.

Mid-infrared prediction of lactoferrin content in bovine milk: potential indicator of mastitis.

Lactoferrin (LTF) is a milk glycoprotein favorably associated with the immune system of dairy cows. Somatic cell count is often used as an indicator of mastitis in dairy cows, but knowledge on the milk LTF content could aid in mastitis detection. An inexpensive, rapid and robust method to predict milk LTF is required. The aim of this study was to develop an equation to quantify the LTF content in bovine milk using mid-infrared (MIR) spectrometry. LTF was quantified by enzyme-linked immunosorbent assay (ELISA), and all milk samples were analyzed by MIR. After discarding samples with a coefficient of variation between 2 ELISA measurements of more than 5% and the spectral outliers, the calibration set consisted of 2499 samples from Belgium (n = 110), Ireland (n = 1658) and Scotland (n = 731). Six statistical methods were evaluated to develop the LTF equation. The best method yielded a cross-validation coefficient of determination for LTF of 0.71 and a cross-validation standard error of 50.55 mg/l of milk. An external validation was undertaken using an additional dataset containing 274 Walloon samples. The validation coefficient of determination was 0.60. To assess the usefulness of the MIR predicted LTF, four logistic regressions using somatic cell score (SCS) and MIR LTF were developed to predict the presence of mastitis. The dataset used to build the logistic regressions consisted of 275 mastitis records and 13 507 MIR data collected in 18 Walloon herds. The LTF and the interaction SCS × LTF effects were significant (P < 0.001 and P = 0.02, respectively). When only the predicted LTF was included in the model, the prediction of the presence of mastitis was not accurate despite a moderate correlation between SCS and LTF (r = 0.54). The specificity and the sensitivity of models were assessed using Walloon data (i.e. internal validation) and data collected from a research herd at the University of Wisconsin - Madison (i.e. 5886 Wisconsin MIR records related to 93 mastistis events - external validation). Model specificity was better when LTF was included in the regression along with SCS when compared with SCS alone. Correct classification of non-mastitis records was 95.44% and 92.05% from Wisconsin and Walloon data, respectively. The same conclusion was formulated from the Hosmer and Lemeshow test. In conclusion, this study confirms the possibility to quantify an LTF indicator from milk MIR spectra. It suggests the usefulness of this indicator associated to SCS to detect the presence of mastitis. Moreover, the knowledge of milk LTF could also improve the milk nutritional quality.

Genetic analysis of lactoferrin content in bovine milk.

Bovine lactoferrin (LF) is mainly present in milk and shows important physiological and biological functions. The aim of this study was to estimate the heritability and correlation values of LF content in bovine milk with different economic traits as milk yield (MY), fat and protein percentages, and somatic cell score (SCS). Variance components of the studied traits were estimated by REML using a multiple-trait mixed model. The obtained heritability (0.22) for LF content predicted using mid-infrared spectrometry (pLF) suggested the possibility of animal selection based on the increase of LF content in milk. The phenotypic and genetic correlation values calculated between pLF and SCS were moderate (0.31 and 0.24, respectively). Furthermore, a preliminary study of bovine LF gene polymorphism effects was performed on the same production traits. By PCR, all exons of the LF gene were amplified and then sequenced. Three new polymorphisms were detected in exon 2, exon 11, and intron 8. We examined the effects of LF gene polymorphisms of exons 2, 4, 9, 11, and 15, and intron 8 on pLF, MY, fat and protein percentages, and SCS. The different observed effects did not reach a significant level probably because of the characteristics of the studied population. However, the results were promising, and LF may be a potential indicator of mastitis. Further studies are necessary to evaluate the effect of genetic selection based on LF content on the improvement of mastitis resistance.

Genomic location of the bovine growth hormone secretagogue receptor (GHSR) gene and investigation of genetic polymorphism.

The growth hormone secretagogue receptor (GHSR) is involved in the regulation of energetic homeostasis and GH secretion. In this study, the bovine GHSR gene was mapped to BTA1 between BL26 and BMS4004. Two different bovine GHSR CDS (GHSR1a and GHSR1b) were sequenced. Six polymorphisms (five SNPs and one 3-bp indel) were also identified, three of them leading to amino acid variations L24V, D194N, and Del R242. These variations are located in the extracellular N-terminal end, the exoloop 2, and the cytoloop 3 of the receptor, respectively.

Genetic variability of lactoferrin content estimated by mid-infrared spectrometry in bovine milk.

The effects of lactoferrin (LF) on the immune system have already been shown by many studies. Unfortunately, the current methods used to measure LF levels in milk do not permit the study of the genetic variability of lactoferrin or the performance of routine genetic evaluations. The first aim of this research was to derive a calibration equation permitting the prediction of LF in milk by mid-infrared spectrometry (MIR). The calibration with partial least squares on 69 samples showed a ratio of standard error of cross-validation to standard deviation equal to 1.98. Based on this value, the calibration equation was used to establish an LF indicator trait (predicted LF; pLF) on a large number of milk samples (n = 7,690). A subsequent study of its variability was conducted, which confirmed that stage of lactation and lactation number influence the overall pLF level. Small differences in mean pLF among 7 dairy breeds were also observed. The pLF content of Jersey milk was significantly higher than that in Holstein milk. Therefore, the choice of breed could change the expected LF level. Heritability estimated for pLF was 19.7%. The genetic and phenotypic correlations between somatic cell score and pLF were 0.04 and 0.26, respectively. As somatic cell score increases in presence of mastitis, this observation seems to indicate that pLF, or a function of observed pLF, compared with expected LF might have potential as an indicator of mastitis. The negative genetic correlation (-0.36) between milk yield and pLF could indicate an undesirable effect of selection for high milk production on the overall LF level.