INVITED REVIEW: Research on ruminal biohydrogenation: Achievements, gaps in knowledge, and future approaches from the perspective of dairy science

Scientific knowledge about ruminal biohydrogenation (BH) has improved greatly since this metabolic process was empirically confirmed in 1951. For years, BH had mostly been perceived as a process to be avoided to increase the post-ruminal flow of UFA from the diet. Two milestones changed this perception and stimulated great interest in BH intermediates themselves: In 1987, the in vitro anticarcinogenic properties of CLA were described, and in 2000, the inhibition of milk fat synthesis by trans-10 cis-12 CLA was confirmed. Since then, numerous BH metabolites have been described in small and large ruminants, and the major deviation from the common BH pathway (i.e., the trans-10 shift) has been reasonably well established. However, there are some less well-characterized alterations, and the comprehensive description of new BH intermediates (e.g., using isotopic tracers) has not been coupled with research on their biological effects. In this regard, the low quality of some published fatty acid profiles may also be limiting the advance of knowledge in BH. Furthermore, although BH seems to no longer be considered a metabolic niche inhabited by a few bacterial species with a highly specific metabolic capability, researchers have failed to elucidate which specific microbial groups are involved in the process and the basis for alterations in BH pathways (i.e., changes in microbial populations or their activity). Unraveling both issues may be beneficial for the description of new microbial enzymes involved in ruminal lipid metabolism that have industrial interest. From the perspective of diary science, other knowledge gaps that require additional research in the coming years are evaluation of the relationship between BH and feed efficiency and enteric methane emissions, as well as improving our understanding of how alterations in BH are involved in milk fat depression. Addressing these issues will have relevant practical implications in dairy science.

Review of equations to predict methane emissions in dairy cows from milk fatty acid profiles and their application to commercial dairy farms

Greenhouse gas emission from the activities of all productive sectors is currently a topic of foremost importance. The major contributors in the livestock sector are ruminants, especially dairy cows. This study aimed to evaluate and compare 21 equations for predicting enteric methane emissions (EME) developed on the basis of milk traits and fatty acid profiles, which were selected from 46 retrieved through a literature review. We compiled a reference database of the detailed fatty acid profiles, determined by GC, of 992 lactating cows from 85 herds under 4 different dairy management systems. The cows were classified according to DIM, parity order, and dairy system. This database was the basis on which we estimated EME using the selected equations. The EME traits estimated were methane yield (20.63 ± 2.26 g/kg DMI, 7 equations), methane intensity (16.05 ± 2.76 g/kg of corrected milk, 4 equations), and daily methane production (385.4 ± 68.2 g/d, 10 equations). Methane production was also indirectly calculated by multiplying the daily corrected milk yield by the methane intensity (416.6 ± 134.7 g/d, 4 equations). We also tested for the effects of DIM, parity, and dairy system (as a correction factor) on the estimates. In general, we observed little consistency among the EME estimates obtained from the different equations, with exception of those obtained from meta-analyses of a range of data from different research centers. We found all the EME predictions to be highly affected by the sources of variation included in the statistical model: DIM significantly affected the results of 19 of the 21 equations, and parity order influenced the results of 13. Different patterns were observed for different equations with only some of them in accordance with expectations based on the cow's physiology. Finally, the best predictions of daily methane production were obtained when a measure of milk yield was included in the equation or when the estimate was indirectly calculated from daily milk yield and methane intensity.

Predicting methane emissions of individual grazing dairy cows from spectral analyses of their milk samples

Data on the enteric methane emissions of individual cows are useful not just in assisting management decisions and calculating herd inventories but also as inputs for animal genetic evaluations. Data generation for many animal characteristics, including enteric methane emissions, can be expensive and time consuming, so being able to extract as much information as possible from available samples or data sources is worthy of investigation. The objective of the present study was to attempt to predict individual cow methane emissions from the information contained within milk samples, specifically the spectrum of light transmittance across different wavelengths of the mid-infrared (MIR) region of the electromagnetic spectrum. A total of 93,888 individual spot measures of methane (i.e., individual samples of an animal's breath when using the GreenFeed technology) from 384 lactations on 277 grazing dairy cows were collapsed into weekly averages expressed as grams per day; each weekly average coincided with a MIR spectral analysis of a morning or evening individual cow milk sample. Associations between the spectra and enteric methane measures were performed separately using partial least squares regression or neural networks with different tuning parameters evaluated. Several alternative definitions of the enteric methane phenotype (i.e., average enteric methane in the 6 d preceding or 6 d following taking the milk sample or the average of the 6 d before and after the milk sample, all of which also included the enteric methane emitted on the day of milk sampling), the candidate model features (e.g., milk yield, milk composition, and milk MIR) as well as validation strategy (i.e., cross-validation or leave-one-experimental treatment-out) were evaluated. Irrespective of the validation method, the prediction accuracy was best when the average of the milk MIR from the morning and evening milk sample was used and the prediction model was developed using neural networks; concurrently including milk yield and days in milk in the prediction model generated superior predictions relative to just the spectral information alone. Furthermore, prediction accuracy was best when the enteric methane phenotype was the average of at least 20 methane spot measures across a 6-d period flanking each side of the milk sample with associated spectral data. Based on the strategy that achieved the best accuracy of prediction, the correlation between the actual and predicted daily methane emissions when based on 4-fold cross-validation varied per validation stratum from 0.68 to 0.75; the corresponding range when validated on each of the 8 different experimental treatments focusing on alternative pasture grazing systems represented in the dataset varied from 0.55 to 0.71. The root mean square error of prediction across the 4-folds of cross-validation was 37.46 g/d, whereas the root mean square error averaged across all folds of leave-one-treatment-out was 37.50 g/d. Results suggest that even with the likely measurement errors contained within the MIR spectrum and gold standard enteric methane phenotype, enteric methane can be reasonably well predicted from the infrared spectrum of milk samples. What is yet to be established, however, is whether (a) genetic variation exists in this predicted enteric methane phenotype and (b) selection on estimates of genetic merit for this phenotype translate to actual phenotypic differences in enteric methane emissions.

Approaches for predicting dairy cattle methane emissions: from traditional methods to machine learning

Measuring dairy cattle methane (CH4) emissions using traditional recording technologies is complicated and expensive. Prediction models, which estimate CH4 emissions based on proxy information, provide an accessible alternative. This review covers the different modeling approaches taken in the prediction of dairy cattle CH4 emissions and highlights their individual strengths and limitations. Following the guidelines set out by the Preferred Reporting Items for Systematic Review and Meta-Analyses (PRISMA); Scopus, EBSCO, Web of Science, PubMed and PubAg were each queried for papers with titles that contained search terms related to a population of "Bovine," exposure of "Statistical Analysis or Machine Learning," and outcome of "Methane Emissions". The search was executed in December 2022 with no publication date range set. Eligible papers were those that investigated the prediction of CH4 emissions in dairy cattle via statistical or machine learning (ML) methods and were available in English. 299 papers were returned from the initial search, 55 of which, were eligible for inclusion in the discussion. Data from the 55 papers was synthesized by the CH4 emission prediction approach explored, including mechanistic modeling, empirical modeling, and machine learning. Mechanistic models were found to be highly accurate, yet they require difficult-to-obtain input data, which, if imprecise, can produce misleading results. Empirical models remain more versatile by comparison, yet suffer greatly when applied outside of their original developmental range. The prediction of CH4 emissions on commercial dairy farms can utilize any approach, however, the traits they use must be procurable in a commercial farm setting. Milk fatty acids (MFA) appear to be the most popular commercially accessible trait under investigation, however, MFA-based models have produced ambivalent results and should be consolidated before robust accuracies can be achieved. ML models provide a novel methodology for the prediction of dairy cattle CH4 emissions through a diverse range of advanced algorithms, and can facilitate the combination of heterogenous data types via hybridization or stacking techniques. In addition to this, they also offer the ability to improve dataset complexity through imputation strategies. These opportunities allow ML models to address the limitations faced by traditional prediction approaches, as well as enhance prediction on commercial farms.

Predicting enteric methane emission in lactating Holsteins based on reference methane data collected by the GreenFeed system

Methane emission is not included in the current breeding goals for dairy cattle mainly due to the expense and difficulty in obtaining sufficient data to generate accurate estimates of the relevant traits. While several models have been developed to predict methane emission from milk spectra using reference methane data obtained by the respiration chamber, SF6 and sniffer methods, the prediction of methane emission from milk mid-infrared (MIR) spectra using reference methane data collected by the GreenFeed system has not yet been explored. Methane emission was monitored for 151 cows using the GreenFeed system. Prediction models were developed for daily and average (for the trial period of 12 or 14 days) methane production (g/d), yield (g/kg DM intake (DMI)) and intensity (g/kg of fat- and protein-corrected milk) using partial least squares regression. The predictions were evaluated in 100 repeated validation cycles, where animals were randomly partitioned into training (80%) and testing (20%) populations for each cycle. The best performing model was observed for average methane intensity using MIR, parity and DMI with validation coefficient of determination (R²_val) and RMSE of prediction of 0.66 and 4.7 g/kg of fat- and protein-corrected milk, respectively. The accuracy of the best models for average methane production and average methane yield were poor (R²_val = 0.28 and 0.12, respectively). A lower accuracy of prediction was observed for methane intensity and production (R²_val = 0.42 and 0.17) when daily records were used while prediction for methane yield was comparable to that for average methane yield (R²_val = 0.16). Our results suggest the potential to predict methane intensity with moderate accuracy. In this case, prediction models for average methane values were generally better than for daily measures when using the GreenFeed system to obtain reference methane emission measurements.

Milk metabolome reveals variations on enteric methane emissions from dairy cows fed a specific inhibitor of the methanogenesis pathway

Metabolome profiling in biological fluids is an interesting approach for exploring markers of methane emissions in ruminants. In this study, a multiplatform metabolomics approach was used for investigating changes in milk metabolic profiles related to methanogenesis in dairy cows. For this purpose, 25 primiparous Holstein cows at similar lactation stage were fed the same diet supplemented with (treated, n = 12) or without (control, n = 13) a specific antimethanogenic additive that reduced enteric methane production by 23% with no changes in intake, milk production, and health status. The study lasted 6 wk, with sampling and measures performed in wk 5 and 6. Milk samples were analyzed using 4 complementary analytical methods, including 2 untargeted (nuclear magnetic resonance and liquid chromatography coupled to a quadrupole-time-of-flight mass spectrometer) and 2 targeted (liquid chromatography-tandem mass spectrometry and gas chromatography coupled to a flame ionization detector) approaches. After filtration, variable selection and normalization data from each analytical platform were then analyzed using multivariate orthogonal partial least square discriminant analysis. All 4 analytical methods were able to differentiate cows from treated and control groups. Overall, 38 discriminant metabolites were identified, which affected 10 metabolic pathways including methane metabolism. Some of these metabolites such as dimethylsulfoxide, dimethylsulfone, and citramalic acid, detected by nuclear magnetic resonance or liquid chromatography-mass spectrometry methods, originated from the rumen microbiota or had a microbial-host animal co-metabolism that could be associated with methanogenesis. Also, discriminant milk fatty acids detected by targeted gas chromatography were mostly of ruminal microbial origin. Other metabolites and metabolic pathways significantly affected were associated with AA metabolism. These findings provide new insight on the potential role of milk metabolites as indicators of enteric methane modifications in dairy cows.

Sequence-based genome-wide association study of individual milk mid-infrared wavenumbers in mixed-breed dairy cattle

BACKGROUND

Fourier-transform mid-infrared (FT-MIR) spectroscopy provides a high-throughput and inexpensive method for predicting milk composition and other novel traits from milk samples. While there have been many genome-wide association studies (GWAS) conducted on FT-MIR predicted traits, there have been few GWAS for individual FT-MIR wavenumbers. Using imputed whole-genome sequence for 38,085 mixed-breed New Zealand dairy cattle, we conducted GWAS on 895 individual FT-MIR wavenumber phenotypes, and assessed the value of these direct phenotypes for identifying candidate causal genes and variants, and improving our understanding of the physico-chemical properties of milk.

RESULTS

Separate GWAS conducted for each of 895 individual FT-MIR wavenumber phenotypes, identified 450 1-Mbp genomic regions with significant FT-MIR wavenumber QTL, compared to 246 1-Mbp genomic regions with QTL identified for FT-MIR predicted milk composition traits. Use of mammary RNA-seq data and gene annotation information identified 38 co-localized and co-segregating expression QTL (eQTL), and 31 protein-sequence mutations for FT-MIR wavenumber phenotypes, the latter including a null mutation in the ABO gene that has a potential role in changing milk oligosaccharide profiles. For the candidate causative genes implicated in these analyses, we examined the strength of association between relevant loci and each wavenumber across the mid-infrared spectrum. This revealed shared association patterns for groups of genomically-distant loci, highlighting clusters of loci linked through their biological roles in lactation and their presumed impacts on the chemical composition of milk.

CONCLUSIONS

This study demonstrates the utility of FT-MIR wavenumber phenotypes for improving our understanding of milk composition, presenting a larger number of QTL and putative causative genes and variants than found from FT-MIR predicted composition traits. Examining patterns of significance across the mid-infrared spectrum for loci of interest further highlighted commonalities of association, which likely reflects the physico-chemical properties of milk constituents.

Improving robustness and accuracy of predicted daily methane emissions of dairy cows using milk mid-infrared spectra

BACKGROUND

A robust proxy for estimating methane (CH₄ ) emissions of individual dairy cows would be valuable especially for selective breeding. This study aimed to improve the robustness and accuracy of prediction models that estimate daily CH₄ emissions from milk Fourier transform mid-infrared (FT-MIR) spectra by (i) increasing the reference dataset and (ii) adjusting for routinely recorded phenotypic information. Prediction equations for CH₄ were developed using a combined dataset including daily CH₄ measurements (n = 1089; g d^-1 ) collected using the SF₆ tracer technique (n = 513) and measurements using respiration chambers (RC, n = 576). Furthermore, in addition to the milk FT-MIR spectra, the variables of milk yield (MY) on the test day, parity (P) and breed (B) of cows were included in the regression analysis as explanatory variables.

RESULTS

Models developed based on a combined RC and SF₆ dataset predicted the expected pattern in CH₄ values (in g d^-1 ) during a lactation cycle, namely an increase during the first weeks after calving followed by a gradual decrease until the end of lactation. The model including MY, P and B information provided the best prediction results (cross-validation statistics: R² = 0.68 and standard error = 57 g CH₄ d^-1 ).

CONCLUSIONS

The models developed accounted for more of the observed variability in CH₄ emissions than previously developed models and thus were considered more robust. This approach is suitable for large-scale studies (e.g. animal genetic evaluation) where robustness is paramount for accurate predictions across a range of animal conditions. © 2020 Society of Chemical Industry.

Untargeted metabolomic investigate milk and ruminal fluid of Holstein cows supplemented with Perilla frutescens leaf

Milk compounds are important for human nutrient requirements and health. The ruminal metabolic profile is responsible for dietary nutrition and determines milk production. Perilla frutescens leaf (PFL) is a commonly used medicinal herb due to its bioactive metabolites. This study elucidated the effects of PFL on the metabolome of two biofluids (rumen fluid and milk) of 14 cows fed a basic total mixed ration diet (CON, n = 7) and supplemented with 300 g/d PFL per cow (PFL, n = 7) by ultra-performance liquid chromatography coupled to quadrupole time-of-flight mass spectrometry. Milk PE-NMe (18:1(9Z)/18:1(9Z)) and DG (18:0/20:4(5Z,8Z,11Z,14Z)/0:0), oleanolic acid, and nucleotides were upregulated, and milk medium-chain fatty acids (2-hydroxycaprylic acid) were down-regulated in response to PFL. The supplementation of PFL increased the abundance of pyrimidine nucleotides both in rumen fluid and milk. The pathways of pyrimidine metabolism and biosynthesis of unsaturated fatty acids were enriched both in the rumen fluid and milk. We also found the milk 2-hydroxycaprylic acid was positively correlated with ruminal uridine 5-monophosphate, and was negatively correlated with ruminal deoxycytidine, and the milk thymidine was positively correlated with ruminal icosenoic acid. This study found that the supplementation of PFL could alter the ruminal metabolic profiles and milk synthesis through regulation of the pathways of pyrimidine metabolism and biosynthesis of unsaturated fatty acids. Our new findings provide comprehensive insights into the metabolomics profile of rumen fluid and milk, supporting the potential production of Perilla frutescens milk in dairy cows.

Predicting the likelihood of conception to first insemination of dairy cows using milk mid-infrared spectroscopy

The objective of this study was to examine the ability of milk mid-infrared (MIR) spectroscopy and other on-farm data, such as milk yield, milk composition, stage of lactation, calving age, days in milk at insemination, and somatic cell count, to identify cows that were most or least likely to conceive to first insemination. A total of 16,628 spectral and milk production records of 7,040 cows from 29 commercial dairy herds across 3 Australian states were used. Three models, comprising different explanatory variables, were tested. Model 1 included features that are readily available on farms participating in milk recording, such as milk yield, milk composition, somatic cell count, days from calving to insemination, and calving season. Days in milk and age at calving were incorporated into model 1 to form model 2. In model 3, MIR was added to model 2, but to avoid double counting, milk composition traits of model 2 were removed. The models were first trained on extreme data [i.e., including cows that (1) conceived to first insemination and (2) cows with no conception event recorded and with only 1 insemination]. Then, the models were validated in a fresh data set with all cows regardless of conception outcomes present to test for their ability to identify cows that conceived or did not conceive to first insemination. To do this, we ranked the predicted probability of all cows in the validation set and then selected the top and bottom records in varying proportions from 5 to 40% (i.e., where the model predicted the highest versus lowest likelihood of conception to first insemination, respectively) and compared with the actual values. The model's performance was evaluated through herd-year by herd-year external validation and measured as the proportion of selected records being correct. The results show that when more cows are selected (i.e., descending confidence), the accuracy of the models was reduced, and selecting the 10% of cows with the highest confidence of predictions produces optimal accuracy. Irrespective of the proportions, none of the models could predict cows that conceived to first insemination, with an accuracy around 0.48. When attempting to predict the bottom 10% of cows, which had the least likelihood of conception to first insemination, model 1 had prediction accuracy around 0.64. Compared with model 1, the addition of days in milk and calving age (model 2) resulted in a negligible improvement in prediction accuracy (0.01 to 0.03). Model 3 had the highest prediction accuracy (0.76), which implies that in the models tested, MIR is of primary importance in the prediction of fertility of dairy cows. In conclusion, this study indicates that MIR and other milk recording data could be used to identify cows with potential difficulty in getting pregnant to first insemination with promising accuracy.

Genetic Parameters of Different FTIR-Enabled Phenotyping Tools Derived from Milk Fatty Acid Profile for Reducing Enteric Methane Emissions in Dairy Cattle

This study aimed to infer the genetic parameters of five enteric methane emissions (EME) predicted from milk infrared spectra (13 models). The reference values were estimated from milk fatty acid profiles (chromatography), individual model-cheese, and daily milk yield of 1158 Brown Swiss cows (85 farms). Genetic parameters were estimated, under a Bayesian framework, for EME reference traits and their infrared predictions. Heritability of predicted EME traits were similar to EME reference values for methane yield (CH₄/DM: 0.232-0.317) and methane intensity per kg of corrected milk (CH₄/CM: 0.177-0.279), smaller per kg cheese solids (CH₄/SO: 0.093-0.165), but greater per kg fresh cheese (CH₄/CU: 0.203-0.267) and for methane production (dCH₄: 0.195-0.232). We found good additive genetic correlations between infrared-predicted methane intensities and the reference values (0.73 to 0.93), less favorable values for CH₄/DM (0.45-0.60), and very variable for dCH₄ according to the prediction method (0.22 to 0.98). Easy-to-measure milk infrared-predicted EME traits, particularly CH₄/CM, CH₄/CU and dCH₄, could be considered in breeding programs aimed at the improvement of milk ecological footprint.

Infrared Spectrometry as a High-Throughput Phenotyping Technology to Predict Complex Traits in Livestock Systems

High-throughput phenotyping technologies are growing in importance in livestock systems due to their ability to generate real-time, non-invasive, and accurate animal-level information. Collecting such individual-level information can generate novel traits and potentially improve animal selection and management decisions in livestock operations. One of the most relevant tools used in the dairy and beef industry to predict complex traits is infrared spectrometry, which is based on the analysis of the interaction between electromagnetic radiation and matter. The infrared electromagnetic radiation spans an enormous range of wavelengths and frequencies known as the electromagnetic spectrum. The spectrum is divided into different regions, with near- and mid-infrared regions being the main spectral regions used in livestock applications. The advantage of using infrared spectrometry includes speed, non-destructive measurement, and great potential for on-line analysis. This paper aims to review the use of mid- and near-infrared spectrometry techniques as tools to predict complex dairy and beef phenotypes, such as milk composition, feed efficiency, methane emission, fertility, energy balance, health status, and meat quality traits. Although several research studies have used these technologies to predict a wide range of phenotypes, most of them are based on Partial Least Squares (PLS) and did not considered other machine learning (ML) techniques to improve prediction quality. Therefore, we will discuss the role of analytical methods employed on spectral data to improve the predictive ability for complex traits in livestock operations. Furthermore, we will discuss different approaches to reduce data dimensionality and the impact of validation strategies on predictive quality.

Large-scale phenotyping in dairy sector using milk MIR spectra: Key factors affecting the quality of predictions

Methods and technologies enabling the estimation at large scale of important traits for the dairy sector are of great interest. Those phenotypes are necessary to improve herd management, animal genetic evaluation, and milk quality control. In the recent years, the research was very active to predict new phenotypes from the mid-infrared (MIR) analysis of milk. Models were developed to predict phenotypes such as fine milk composition, milk technological properties or traits related to cow health, fertility and environmental impact. Most of models were developed within research contexts and often not designed for routine use. The implementation of models at a large scale to predict new traits of interest brings new challenges as the factors influencing the robustness of models are poorly documented. The first objective of this work is to highlight the impact on prediction accuracy of factors such as the variability of the spectral and reference data, the spectral regions used and the complexity of models. The second objective is to emphasize methods and indicators to evaluate the quality of models and the quality of predictions generated under routine conditions. The last objective is to outline the issues and the solutions linked with the use and transfer of models on large number of instruments. Based on partial least square regression and 10 datasets including milk MIR spectra and reference quantitative values for 57 traits of interest, the impact of the different factors is illustrated by evaluating the influence on the validation root mean square error of prediction (RMSEP). In the displayed examples, all factors, when well set up, increase the quality of predictions, with an improvement of the RMSEP ranging from 12% to 43%. This work also aims to underline the need for and the complementarity between different validation procedures, statistical parameters and quality assurance methods. Finally, when using and transferring models, the impact of the spectral standardization on the prediction reproducibility is highlighted with an improvement up to 86% with the tested models, and the monitoring of individual spectrometer stability over time appears essential. This list inspired from our experience is of course not exhaustive. The displayed results are only examples and not general rules and other aspects play a role in the quality of final predictions. However, this work highlights good practices, methods and indicators to increase and evaluate quality of phenotypes predicted at a large scale. The results obtained argue for the development of guidelines at international levels, as well as international collaborations in order to constitute large and robust datasets and enable the use of models in routine conditions.

The use of milk Fourier transform mid-infrared spectra and milk yield to estimate heat production as a measure of efficiency of dairy cows

Background

Transformation of feed energy ingested by ruminants into milk is accompanied by energy losses via fecal and urine excretions, fermentation gases and heat. Heat production may differ among dairy cows despite comparable milk yield and body weight. Therefore, heat production can be considered an indicator of metabolic efficiency and directly measured in respiration chambers. The latter is an accurate but time-consuming technique. In contrast, milk Fourier transform mid-infrared (FTIR) spectroscopy is an inexpensive high-throughput method and used to estimate different physiological traits in cows. Thus, this study aimed to develop a heat production prediction model using heat production measurements in respiration chambers, milk FTIR spectra and milk yield measurements from dairy cows.

Methods

Heat production was computed based on the animal’s consumed oxygen, and produced carbon dioxide and methane in respiration chambers. Heat production data included 168 24-h-observations from 64 German Holstein and 20 dual-purpose Simmental cows. Animals were milked twice daily at 07:00 and 16:30 h in the respiration chambers. Milk yield was determined to predict heat production using a linear regression. Milk samples were collected from each milking and FTIR spectra were obtained with MilkoScan FT 6000. The average or milk yield-weighted average of the absorption spectra from the morning and afternoon milking were calculated to obtain a computed spectrum. A total of 288 wavenumbers per spectrum and the corresponding milk yield were used to develop the heat production model using partial least squares (PLS) regression.

Results

Measured heat production of studied animals ranged between 712 and 1470 kJ/kg BW^0.75. The coefficient of determination for the linear regression between milk yield and heat production was 0.46, whereas it was 0.23 for the FTIR spectra-based PLS model. The PLS prediction model using weighted average spectra and milk yield resulted in a cross-validation variance of 57% and a root mean square error of prediction of 86.5 kJ/kg BW^0.75. The ratio of performance to deviation (RPD) was 1.56.

Conclusion

The PLS model using weighted average FTIR spectra and milk yield has higher potential to predict heat production of dairy cows than models applying FTIR spectra or milk yield only.

The evolving role of Fourier-transform mid-infrared spectroscopy in genetic improvement of dairy cattle

Over the last 100 years, significant advances have been made in the characterisation of milk composition for dairy cattle improvement programs. Technological progress has enabled a shift from labour intensive, on-farm collection and processing of samples that assess yield and fat levels in milk, to large-scale processing of samples through centralised laboratories, with the scope extended to include quantification of other traits. Fourier-transform mid-infrared (FT-MIR) spectroscopy has had a significant role in the transformation of milk composition phenotyping, with spectral-based predictions of major milk components already being widely used in milk payment and animal evaluation systems globally. Increasingly, there is interest in analysing the individual FT-MIR wavenumbers, and in utilising the FT-MIR data to predict other novel traits of importance to breeding programs. This includes traits related to the nutritional value of milk, the processability of milk into products such as cheese, and traits relevant to animal health and the environment. The ability to successfully incorporate these traits into breeding programs is dependent on the heritability of the FT-MIR predicted traits, and the genetic correlations between the FT-MIR predicted and actual trait values. Linking FT-MIR predicted traits to the underlying mutations responsible for their variation can be difficult because the phenotypic expression of these traits are a function of a diverse range of molecular and biological mechanisms that can obscure their genetic basis. The individual FT-MIR wavenumbers give insights into the chemical composition of milk and provide an additional layer of granularity that may assist with establishing causal links between the genome and observed phenotypes. Additionally, there are other molecular phenotypes such as those related to the metabolome, chromatin accessibility, and RNA editing that could improve our understanding of the underlying biological systems controlling traits of interest. Here we review topics of importance to phenotyping and genetic applications of FT-MIR spectra datasets, and discuss opportunities for consolidating FT-MIR datasets with other genomic and molecular data sources to improve future dairy cattle breeding programs.

Predicting milk fatty acids and energy balance of dairy cows in Australia using milk mid-infrared spectroscopy

Mid-infrared spectroscopy (MIRS) is traditionally used for analysing milk fat, protein and lactose concentrations in dairy production, but there is growing interest in using it to predict difficult, or expensive-to-measure, phenotypes on a large scale. The resulting prediction equations can be applied to MIRS data from commercial herd-testing, to facilitate management and feeding decisions, or for genomic selection purposes. We investigated the ability of MIRS of milk samples to predict milk fatty acids (FAs) and energy balance (EB) of dairy cows in Australia. Data from 240 Holstein lactating cows that were part of two 32-day experiments, were used. Milk FAs were measured twice during the experimental period. Prediction models were developed using partial least-square regression with a 10-fold cross-validation. Measures of prediction accuracy included the coefficient of determination (R2cv) and root mean-square error. Milk FAs with a chain length of ≤16 were accurately predicted (0.89 ≤ R2cv ≤ 0.95), while prediction accuracy for FAs with a chain length of ≥17 was slightly lower (0.72 ≤ R2cv ≤ 0.82). The accuracy of the model prediction was moderate for EB, with the value of R2cv of 0.48. In conclusion, the ability of MIRS to predict milk FAs was high, while EB was moderately predicted. A larger dataset is needed to improve the accuracy and the robustness of the prediction models.

Enteric Methane Emissions of Dairy Cows Predicted from Fatty Acid Profiles of Milk, Cream, Cheese, Ricotta, Whey, and Scotta

Enteric methane emissions (EME) of ruminants contribute to global climate change, but any attempt to reduce it will need an easy, inexpensive, and accurate method of quantification. We used a promising indirect method for estimating EMEs of lactating dairy cows based on the analysis of the fatty acid (FA) profile of their milk. The aim of this preliminary study was to assess milk from four single samplings (morning whole, evening whole, evening partially skimmed, and vat milks) as alternatives to reference whole milk samples from two milkings. Three fresh products (cream, cheese, and ricotta), two by-products (whey and scotta), and two long-ripened cheeses (6 and 12 months) were also assessed as alternative sources of information to reference milk. The 11 alternative matrices were obtained from seven experimental cheese- and ricotta-making sessions carried out every two weeks following the artisanal Malga cheese-making procedure using milk from 148 dairy cows kept on summer highland pastures. A total of 131 samples of milk, dairy products, and by-products were analyzed to determine the milk composition and to obtain detailed FA profiles using bi-dimensional gas-chromatography. Two equations taken from a published meta-analysis of methane emissions measured in the respiration chambers of cows on 30 different diets were applied to the proportions of butyric, iso-palmitic, iso-oleic, vaccenic, oleic, and linoleic acids out of total FAs to predict methane yield per kg of dry matter ingested and methane intensity per kg of fat and protein corrected milk produced by the cows. Methane yield and intensity could be predicted from single milk samples with good accuracy (trueness and precision) with respect to those predicted from reference milk. The fresh products (cream, cheese and ricotta) generally showed good levels of trueness but low precision for predicting both EME traits, which means that a greater number of samples needs to be analyzed. Among by-products, whey could be a viable alternative source of information for predicting both EME traits, whereas scotta overestimated both traits and showed low precision (due also to its very low fat content). Long-ripened cheeses were found to be less valuable sources of information, although six-month cheese could, with specific correction factors, be acceptable sources of information for predicting the methane yield of lactating cows. These preliminary results need to be confirmed by further study on different dairy systems and cheese-making technologies but offer new insight into a possible easy method for monitoring the EME at the field level along the dairy chain.

Classifying the fertility of dairy cows using milk mid-infrared spectroscopy

The objective of this study was to investigate the potential of milk mid-infrared (MIR) spectroscopy, MIR-derived traits including milk composition, milk fatty acids, and blood metabolic profiles (fatty acids, β-hydroxybutyrate, and urea), and other on-farm data for discriminating cows of good versus poor likelihood of conception to first insemination (i.e., pregnant vs. open). A total of 6,488 spectral and milk production records of 2,987 cows from 19 commercial dairy herds across 3 Australian states were used. Seven models, comprising different explanatory variables, were examined. Model 1 included milk production; concentrations of fat, protein, and lactose; somatic cell count; age at calving; days in milk at herd test; and days from calving to insemination. Model 2 included, in addition to the variables in model 1, milk fatty acids and blood metabolic profiles. The MIR spectrum collected before first insemination was added to model 2 to form model 3. Fat, protein, and lactose percentages, milk fatty acids, and blood metabolic profiles were removed from model 3 to create model 4. Model 5 and model 6 comprised model 4 and either fertility genomic estimated breeding value or principal components obtained from a genomic relationship matrix derived using animal genotypes, respectively. In model 7, all previously described sources of information, but not MIR-derived traits, were used. The models were developed using partial least squares discriminant analysis. The performance of each model was evaluated in 2 ways: 10-fold random cross-validation and herd-by-herd external validation. The accuracy measures were sensitivity (i.e., the proportion of pregnant cows that were correctly classified), specificity (i.e., the proportion of open cows that were correctly classified), and area under the curve (AUC) for the receiver operating curve. The results showed that in all models, prediction accuracy obtained through 10-fold random cross-validation was higher than that of herd-by-herd external validation, with the difference in AUC ranging between 0.01 and 0.09. In the herd-by-herd external validation, using basic on-farm information (model 1) was not sufficient to classify good- and poor-fertility cows; the sensitivity, specificity, and AUC were around 0.66. Compared with model 1, adding milk fatty acids and blood metabolic profiles (model 2) increased the sensitivity, specificity, and AUC by 0.01, 0.02, and 0.02 unit, respectively (i.e., 0.65, 0.63, and 0.678). Incorporating MIR spectra into model 2 resulted in sensitivity, specificity, and AUC values of 0.73, 0.63, and 0.72, respectively (model 3). The comparable prediction accuracies observed for models 3 and 4 mean that useful information from MIR-derived traits is already included in the spectra. Adding the fertility genomic estimated breeding value and animal genotypes (model 7) produced the highest prediction accuracy, with sensitivity, specificity, and AUC values of 0.75, 0.66, and 0.75, respectively. However, removing either the fertility estimated breeding value or animal genotype from model 7 resulted in a reduction of the prediction accuracy of only 0.01 and 0.02, respectively. In conclusion, this study indicates that MIR and other on-farm data could be used to classify cows of good and poor likelihood of conception with promising accuracy.

Validation strategy can result in an overoptimistic view of the ability of milk infrared spectra to predict methane emission of dairy cattle

Because of the environmental impact of methane (CH₄), it is of great interest to reduce CH₄ emission of dairy cattle and selective breeding might contribute to this. However, this approach requires a rapid and inexpensive measurement technique that can be used to quantify CH₄ emission for a large number of individual dairy cows. Milk infrared (IR) spectroscopy has been proposed as a predictor for CH₄ emission. In this study, we investigated the feasibility of milk IR spectra to predict breath sensor-measured CH₄ of 801 dairy cows on 10 commercial farms. To evaluate the prediction equation, we used random and block cross validation. Using random cross validation, we found a validation coefficient of determination (R²val) of 0.49, which suggests that milk IR spectra are informative in predicting CH₄ emission. However, based on block cross validation, with farms as blocks, a negligible R²val of 0.01 was obtained, indicating that milk IR spectra cannot be used to predict CH₄ emission. Random cross validation thus results in an overoptimistic view of the ability of milk IR spectra to predict CH₄ emission of dairy cows. The difference between the validation strategies could be due to the confounding of farm and date of milk IR analysis, which introduces a correlation between batch effects on the IR analyses and farm-average CH₄. Breath sensor-measured CH₄ is strongly influenced by farm-specific conditions, which magnifies the problem. Milk IR wavenumbers from water absorption regions, which are generally considered uninformative, showed moderate accuracy (R²val = 0.25) when based on random cross validation, but not when based on block cross validation (R²val = 0.03). These results indicate, therefore, that in the current study, random cross validation results in an overoptimistic view on the ability of milk IR spectra to predict CH₄ emission. We suggest prediction based on wavenumbers from water absorption regions as a negative control to identify potential dependence structures in the data.