Evidence of additive genetic variation for major milk proteins in dairy cows: A meta-analysis

In the past, there have been reports of genetic parameters for milk proteins in various dairy cattle populations. The high variability among genetic parameter estimates has been caused by this. This study aimed to use a random-effects meta-analysis model to compile published estimates of genetic parameter for major milk proteins of α-lactalbumin, β-lactoglobulin, sum of whey proteins, casein, αs1-casein, αs2-casein, β-casein, and κ-casein in dairy cows. The study used a total of 140 heritability and 256 genetic correlation estimates from 23 papers published between 2004 and 2022. The estimated range of milk protein heritability is from 0.284 (for α-lactalbumin in milk) to 0.596 (for sum of whey proteins). The genetic correlation estimates between casein and milk yield, milk fat and protein percentages were -0.461, 0.693, and 0.976, respectively (p < 0.05). The genetic correlation estimates between milk proteins expressed as a percentage of milk were significant and varied from 0.177 (between β-lactoglobulin and κ-casein) to 0.892 (between αs1-casein and αs2-casein). Moderate-to-high heritability estimates for milk proteins and their low genetic associations with milk yield and composition indicated the possibility for improving milk proteins in a genetic selection plan with negligible correlated effects on production traits in dairy cows.

Prediction of key milk biomarkers in dairy cows through milk mid-infrared spectra and international collaborations

At the individual cow level, suboptimum fertility, mastitis, negative energy balance, and ketosis are major issues in dairy farming. These problems are widespread on dairy farms and have an important economic impact. The objectives of this study were (1) to assess the potential of milk mid-infrared (MIR) spectra to predict key biomarkers of energy deficit (citrate, isocitrate, glucose-6 phosphate [glucose-6P], free glucose), ketosis (β-hydroxybutyrate [BHB] and acetone), mastitis (N-acetyl-β-d-glucosaminidase activity [NAGase] and lactate dehydrogenase), and fertility (progesterone); (2) to test alternative methodologies to partial least squares (PLS) regression to better account for the specific asymmetric distribution of the biomarkers; and (3) to create robust models by merging large datasets from 5 international or national projects. Benefiting from this international collaboration, the dataset comprised a total of 9,143 milk samples from 3,758 cows located in 589 herds across 10 countries and represented 7 breeds. The samples were analyzed by reference chemistry for biomarker contents, whereas the MIR analyses were performed on 30 instruments from different models and brands, with spectra harmonized into a common format. Four quantitative methodologies were evaluated to address the strongly skewed distribution of some biomarkers. Partial least squares regression was used as the reference basis, and compared with a random modification of distribution associated with PLS (random-downsampling-PLS), an optimized modification of distribution associated with PLS (KennardStone-downsampling-PLS), and support vector machine (SVM). When the ability of MIR to predict biomarkers was too low for quantification, different qualitative methodologies were tested to discriminate low versus high values of biomarkers. For each biomarker, 20% of the herds were randomly removed within all countries to be used as the validation dataset. The remaining 80% of herds were used as the calibration dataset. In calibration, the 3 alternative methodologies outperform the PLS performances for the majority of biomarkers. However, in the external herd validation, PLS provided the best results for isocitrate, glucose-6P, free glucose, and lactate dehydrogenase (coefficient of determination in external herd validation [R2v] = 0.48, 0.58, 0.28, and 0.24, respectively). For other molecules, PLS-random-downsampling and PLS-KennardStone-downsampling outperformed PLS in the majority of cases, but the best results were provided by SVM for citrate, BHB, acetone, NAGase, and progesterone (R2v = 0.94, 0.58, 0.76, 0.68, and 0.15, respectively). Hence, PLS and SVM based on the entire dataset provided the best results for normal and skewed distributions, respectively. Complementary to the quantitative methods, the qualitative discriminant models enabled the discrimination of high and low values for BHB, acetone, and NAGase with a global accuracy around 90%, and glucose-6P with an accuracy of 83%. In conclusion, MIR spectra of milk can enable quantitative screening of citrate as a biomarker of energy deficit and discrimination of low and high values of BHB, acetone, and NAGase, as biomarkers of ketosis and mastitis. Finally, progesterone could not be predicted with sufficient accuracy from milk MIR spectra to be further considered. Consequently, MIR spectrometry can bring valuable information regarding the occurrence of energy deficit, ketosis, and mastitis in dairy cows, which in turn have major influences on their fertility and survival.

Predicting milk protein fractions using infrared spectroscopy and a gradient boosting machine for breeding purposes in Holstein cattle

In recent years, increasing attention has been focused on the genetic evaluation of protein fractions in cow milk with the aim of improving milk quality and technological characteristics. In this context, advances in high-throughput phenotyping by Fourier transform infrared (FTIR) spectroscopy offer the opportunity for large-scale, efficient measurement of novel traits that can be exploited in breeding programs as indicator traits. We took milk samples from 2,558 Holstein cows belonging to 38 herds in northern Italy, operating under different production systems. Fourier transform infrared spectra were collected on the same day as milk sampling and stored for subsequent analysis. Two sets of data (i.e., phenotypes and FTIR spectra) collected in 2 different years (2013 and 2019-2020) were compiled. The following traits were assessed using HPLC: true protein, major casein fractions [α_S1-casein (CN), α_S2-CN, β-CN, κ-CN, and glycosylated-κ-CN], and major whey proteins (β-lactoglobulin and α-lactalbumin), all of which were measured both in grams per liter (g/L) and proportion of total nitrogen (% N). The FTIR predictions were calculated using the gradient boosting machine technique and tested by 3 different cross-validation (CRV) methods. We used the following CRV scenarios: (1) random 10-fold, which randomly split the whole into 10-folds of equal size (9-folds for training and 1-fold for validation); (2) herd/date-out CRV, which assigned 80% of herd/date as the training set with independence of 20% of herd/date assigned as the validation set; (3) forward/backward CRV, which split the data set in training and validation set according with the year of milk sampling (FTIR and gold standard data assessed in 2013 or 2019-2020) using the "old" and "new" databases for training and validation, and vice-versa with independence among them; (4) the CRV for genetic parameters (CRV-gen), where animals without pedigree as assigned as a fixed training population and animals with pedigree information was split in 5-folds, in which 1-fold was assigned to the fixed training population, and 4-folds were assigned to the validation set (independent from the training set). The results (i.e., measures and predictions) of CRV-gen were used to infer the genetic parameters for gold standard laboratory measurements (i.e., proteins assessed with HPLC) and FTIR-based predictions considering the CRV-gen scenario from a bi-trait animal model using single-step genomic BLUP. We found that the prediction accuracies of the gradient boosting machine equations differed according to the way in which the proteins were expressed, achieving higher accuracy when expressed in g/L than when expressed as % N in all CRV scenarios. Concerning the reproducibility of the equations over the different years, the results showed no relevant differences in predictive ability between using "old" data as the training set and "new" data as the validation set and vice-versa. Comparing the additive genetic variance estimates for milk protein fractions between the FTIR predicted and HPLC measures, we found reductions of -19.7% for milk protein fractions expressed in g/L, and -21.19% expressed as % N. Although we found reductions in the heritability estimates, they were small, with values ranging from -1.9 to -7.25% for g/L, and -1.6 to -7.9% for % N. The posterior distributions of the additive genetic correlations (r_a) between the FTIR predictions and the laboratory measurements were generally high (>0.8), even when the milk protein fractions were expressed as % N. Our results show the potential of using FTIR predictions in breeding programs as indicator traits for the selection of animals to enhance milk protein fraction contents. We expect acceptable responses to selection due to the high genetic correlations between HPLC measurements and FTIR predictions.

Breeding for improved protein fractions and free amino acids concentration in bovine milk

Considerable resources are required to routinely measure detailed milk compositional traits. Hence, an insufficient volume of phenotypic data can hinder genetic progress in these traits within dairy cow breeding programmes. The objective of the present study was to quantify the opportunities for breeding for improved milk protein and free amino acid (FAA) composition by exploiting mid‐infrared spectroscopy (MIRS) predictions routinely recorded from milk samples. Genetic parameters for protein fractions and FAA composition were estimated using 134,546 test‐day records from 16,166 lactations on 9,572 cows using linear mixed models. Heritability of MIRS‐predicted protein fractions ranged from 0.19 (α‐lactalbumin) to 0.55 (β‐lactoglobulin A), while heritability of MIRS‐predicted FAA ranged from 0.08 for glycine to 0.29 for glutamic acid. Genetic correlations among the MIRS‐predicted FAA were moderate to strong ranging from −0.44 (aspartic acid and lysine) to 0.97 (glutamic acid and total FAA). Adjustment of the genetic correlations for the genetic merit of 24‐h milk yield did not greatly affect the correlations. Results from the current study highlight the presence of exploitable genetic variation for both protein fractions and FAA in dairy cow milk. Besides, the direction of genetic correlations reveals that breeding programmes directly selecting for greater milk protein concentration carry with them favourable improvement in casein and whey fractions.

Infrared Predictions Are a Valuable Alternative to Actual Measures of Dry-Cured Ham Weight Loss in the Training of Genome-Enabled Prediction Models

Selection to reduce ham weight losses during dry-curing (WL) requires individual traceability of hams throughout dry-curing, with high phenotyping costs and long generation intervals. Infrared spectroscopy enables cost-effective, high-throughput phenotyping for WL 24 h after slaughter. Direct genomic values (DGV) of crossbred pigs and their purebred sires were estimated, for observed (OB) and infrared-predicted WL (IR), through models developed from 640 and 956 crossbred pigs, respectively. Five Bayesian models and two pseudo-phenotypes (estimated breeding value, EBV, and adjusted phenotype) were tested in random cross-validation and leave-one-family-out validation. The use of EBV as pseudo-phenotypes resulted in the highest accuracies. Accuracies in leave-one-family-out validation were much lower than those obtained in random cross-validation but still satisfactory and very similar for both traits. For sires in the leave-one-family-out validation scenario, the correlation between the DGV for IR and EBV for OB was slightly lower (0.32) than the correlation between the DGV for OB and EBV for OB (0.38). While genomic prediction of OB and IR can be equally suggested to be incorporated in future selection programs aiming at reducing WL, the use of IR enables an early, cost-effective phenotyping, favoring the construction of larger reference populations, with accuracies comparable to those achievable using OB phenotype.

A meta-analysis of genetic parameter estimates for milk and serum minerals in dairy cows

This study aimed to conduct a meta-analysis based on a random-effects model to combine different published heritability estimates and genetic correlations for milk and serum minerals in dairy cows. In total, 59 heritability and 25 genetic correlation estimates from 12 articles published between 2009 and 2021 were used. The heritability estimates for milk macro-minerals were moderate to high and ranged from 0.311 (for Na) to 0.420 (for Ca). On the other hand, milk micro-minerals had lower heritabilities with a range from 0.013 (for Fe) to 0.373 (for Zn). The heritability estimates for serum macro-minerals were generally low and varied from 0.126 (for K) to 0.206 (for Mg). The estimates of genetic correlation between milk macro-minerals varied from -0.024 (between Na and K) to 0.625 (between Mg and P). The genetic correlations of milk Ca and P with milk yield were -0.171 and -0.211, respectively. The estimates of genetic parameters reported in this meta-analysis study are appropriate to utilize in breeding plans when valid estimates are not available for milk minerals in dairy cow populations.

Application of Optical Quality Control Technologies in the Dairy Industry: An Overview

Sustainable development of the agricultural industry, in particular, the production of milk and feed for farm animals, requires accurate, fast, and non-invasive diagnostic tools. Currently, there is a rapid development of a number of analytical methods and approaches that meet these requirements. Infrared spectrometry in the near and mid-IR range is especially widespread. Progress has been made not only in the physical methods of carrying out measurements, but significant advances have also been achieved in the development of mathematical processing of the received signals. This review is devoted to the comparison of modern methods and devices used to control the quality of milk and feed for farm animals.

Effects of breed, farm intensiveness, and cow productivity level on cheese-making ability predicted using infrared spectral data at the population level

Fourier-transform infrared (FTIR) spectra collected during milk recording schemes at population level can be used for predicting novel traits of interest for farm management, cows' genetic improvement, and milk payment systems. The aims of this study were as follows. (1) To predict cheese yield traits using FTIR spectra from routine milk recordings exploiting previously developed calibration equations. (2) To compare the predicted cheese-making abilities of different dairy and dual-purpose breeds. (3) To analyze the effects of herds' level of intensiveness (HL) and of the cow's level of productivity (CL). (4) To compare the patterns of predicted cheese yields with the patterns of milk composition in different breeds to discern the drivers of cheese-making efficiency. The major sources of variation of FTIR predictions of cheese yield ability (fresh cheese or cheese solids produced per unit milk) of individual milk samples were studied on 115,819 cows of 4 breeds (2 specialized dairy breeds, Holstein and Brown Swiss, and 2 dual-purpose breeds, Simmental and Alpine Grey) from 6,430 herds and exploiting 1,759,706 FTIR test-day spectra collected over 7 yr of milk sampling. Calibration equations used were previously developed on 1,264 individual laboratory model cheese procedures (cross-validation R² 0.85 and 0.95 for fresh and solids cheese yields, respectively). The linear model used for statistical analysis included the effects of parity, lactation stage, year of calving, month of sampling, HL, CL, breed of cow, and the interactions breed × HL and breed × CL. The HL and CL stratifications (5 classes each) were based on average daily secretion of milk net energy per cow. All effects were highly significant (P < 0.001). The major conclusions were as follows. (1) The FTIR-based prediction of cheese yield of milk goes beyond the knowledge of fat and protein content, partially explaining differences in cheese-making ability in different cows, breeds and herds. (2) Differences in cheese yields of different breeds are only partially explained by milk fat and protein composition, and less productive breeds are characterized by a higher milk nutrient content as well as a higher recovery of nutrients in the cheese. (3) High-intensive herds not only produce much more milk, but the milk has a higher nutrient content and a higher cheese yield, whereas within herds, compared with less productive cows, the more productive cows have a much greater milk yield, milk with a greater content of fat but not of protein, and a moderate improvement in cheese yield, differing little from expectations based on milk composition. Finally, (4) the effects of HL and CL on milk quality and cheese-making ability are similar but not identical in different breeds, the less productive ones having some advantage in terms of cheese-making ability. We can obtain FTIR-based prediction of cheese yield from individual milk samples retrospectively at population level, which seems to go beyond the simple knowledge of milk composition, incorporating information on nutrient retention ability in cheese, with possible advantages for management of farms, genetic improvement of dairy cows, and milk payment systems.

Meta-analysis-based estimates of efficiency of calcium utilisation by ruminants

The most abundant mineral in the body of animals is Ca, which has crucial importance for the regulation of various processes. The maintenance of Ca balance has become more challenging, especially in lactating ruminants, owing to the increased milk yields and thus Ca requirement. To determine the Ca requirement, factors such as Ca secretion via milk or Ca deposition in body tissues and conception products are summed up to the net Ca requirement. Nevertheless, dietary Ca cannot be completely utilised by the animal to cover the net Ca requirement, therefore a value for the efficiency of Ca utilisation is applied, which is the maximum proportion of Ca from the feed that the animal can use for covering the net requirement. However, current estimates for the efficiency of Ca utilisation are inconsistent. Therefore, the objective of the present meta-analysis was to estimate the efficiency of Ca utilisation for ruminants, considering the Ca supply of the animal. A data set of 223 observations was compiled from 37 studies, including data on cattle and small ruminants. Standardised Ca digestibility was calculated from data on Ca intake and faecal Ca excretion, corrected for faecal endogenous losses. Furthermore, a data subset on only lactating ruminants was created. For this subset, Ca excretion via faeces and urine and standardised Ca digestibility were related to the Ca supply of the animal. An exponential function was fitted to standardised digestibility data in response to Ca concentration in the diet and Ca supply, revealing that standardised Ca digestibility decreased with increasing dietary Ca concentration and Ca supply. The median for standardised Ca digestibility was 40%, with a remarkable variation between 9% and 88%. In response to Ca supply, faecal Ca excretion increased in a strong linear manner (slope = 0.76, R² = 0.96). Excretion of Ca via urine was very low even when Ca supply was very high. To conclude, Ca digestibility is a suitable indicator for the efficiency of Ca utilisation, since excessive Ca is almost completely excreted in faeces; however, Ca digestibility has to be determined at a Ca supply level below the requirement of the animal. To date, only very limited number of data have been reported for such supply conditions. Comparative studies using various Ca sources are suggested for future studies but should be conducted using a marginal Ca supply level.

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.

Genetic analysis on infrared-predicted milk minerals for Danish dairy cattle

A group of milk components that has shown potential to be predicted with milk spectra is milk minerals. Milk minerals are important for human health and cow health. Having an inexpensive and fast way to measure milk mineral concentrations would open doors for research, herd management, and selective breeding. The first aim of this study was to predict milk minerals with infrared milk spectra. Additionally, milk minerals were predicted with infrared-predicted fat, protein, and lactose content. The second aim was to perform a genetic analysis on infrared-predicted milk minerals, to identify QTL, and estimate variance components. For training and validating a multibreed prediction model for individual milk minerals, 264 Danish Jersey cows and 254 Danish Holstein cows were used. Partial least square regression prediction models were built for Ca, Cu, Fe, K, Mg, Mn, Na, P, Se, and Zn based on 80% of the cows, selected randomly. Prediction models were externally validated with 8 herds based on the remaining 20% of the cows. The prediction models were applied on a population of approximately 1,400 Danish Holstein cows with 5,600 infrared spectral records and 1,700 Danish Jersey cows with 7,200 infrared spectral records. Cows from this population had 50k imputed genotypes. Prediction accuracy was good for P and Ca, with external R² ≥ 0.80 and a relative prediction error of 5.4% for P and 6.3% for Ca. Prediction was moderately good for Na with an external R² of 0.63, and a relative error of 18.8%. Prediction accuracies of milk minerals based on infrared-predicted fat, protein, and lactose content were considerably lower than those based on the infrared milk spectra. This shows that the milk infrared spectrum contains valuable information on milk minerals, which is currently not used. Heritability for infrared-predicted Ca, Na, and P varied from low (0.13) to moderate (0.36). Several QTL for infrared-predicted milk minerals were observed that have been associated with gold standard milk minerals previously. In conclusion, this study has shown infrared milk spectra were good at predicting Ca, Na, and P in milk. Infrared-predicted Ca, Na, and P had low to moderate heritability estimates.

Sequence-based GWAS and post-GWAS analyses reveal a key role of SLC37A1, ANKH, and regulatory regions on bovine milk mineral content

The mineral composition of bovine milk plays an important role in determining its nutritional and cheese-making value. Concentrations of the main minerals predicted from mid-infrared spectra produced during milk recording, combined with cow genotypes, provide a unique opportunity to decipher the genetic determinism of these traits. The present study included 1 million test-day predictions of Ca, Mg, P, K, Na, and citrate content from 126,876 Montbéliarde cows, of which 19,586 had genotype data available. All investigated traits were highly heritable (0.50-0.58), with the exception of Na (0.32). A sequence-based genome-wide association study (GWAS) detected 50 QTL (18 affecting two to five traits) and positional candidate genes and variants, mostly located in non-coding sequences. In silico post-GWAS analyses highlighted 877 variants that could be regulatory SNPs altering transcription factor (TF) binding sites or located in non-coding RNA (mainly lncRNA). Furthermore, we found 47 positional candidate genes and 45 TFs highly expressed in mammary gland compared to 90 other bovine tissues. Among the mammary-specific genes, SLC37A1 and ANKH, encoding proteins involved in ion transport were located in the most significant QTL. This study therefore highlights a comprehensive set of functional candidate genes and variants that affect milk mineral content.

Prediction of meat quality traits in the abattoir using portable near-infrared spectrometers: heritability of predicted traits and genetic correlations with laboratory-measured traits

Background

The possibility of assessing meat quality traits over the meat chain is strongly limited, especially in the context of selective breeding which requires a large number of phenotypes. The main objective of this study was to investigate the suitability of portable infrared spectrometers for phenotyping beef cattle aiming to genetically improving the quality of their meat. Meat quality traits (pH, color, water holding capacity, tenderness) were appraised on rib eye muscle samples of 1,327 Piemontese young bulls using traditional (i.e., reference/gold standard) laboratory analyses; the same traits were also predicted from spectra acquired at the abattoir on the intact muscle surface of the same animals 1 d after slaughtering. Genetic parameters were estimated for both laboratory measures of meat quality traits and their spectra-based predictions.

Results

The prediction performances of the calibration equations, assessed through external validation, were satisfactory for color traits (R² from 0.52 to 0.80), low for pH and purge losses (R² around 0.30), and very poor for cooking losses and tenderness (R² below 0.20). Except for lightness and purge losses, the heritability estimates of most of the predicted traits were lower than those of the measured traits while the genetic correlations between measured and predicted traits were high (average value 0.81).

Conclusions

Results showed that NIRS predictions of color traits, pH, and purge losses could be used as indicator traits for the indirect genetic selection of the reference quality phenotypes. Results for cooking losses were less effective, while the NIR predictions of tenderness were affected by a relatively high uncertainty of estimate. Overall, genetic selection of some meat quality traits, whose direct phenotyping is difficult, can benefit of the application of infrared spectrometers technology.

On-site visible-near IR prediction of iodine number and fatty acid composition of subcutaneous fat of raw hams as phenotypes for a heavy pig breeding program

The quality of subcutaneous fat of raw hams is a trait of interest in selective breeding programs for pig lines used in dry-cured ham production, and rapid, non-invasive methods for its assessment are available. However, the efficacy of such methods to provide indicator traits for breeding programs needs to be proven. The study investigated the accuracy of on-site visible-near IR spectroscopy predictions of iodine number and fatty acid (FA) composition of raw ham subcutaneous fat, and it evaluated their effectiveness as indicator traits of ham fat quality in a pig breeding program. Prediction equations were developed using visible-near IR spectra acquired at the slaughterhouse from five sites in subcutaneous fat of raw hams of 1025 crossbred pigs. Pigs were raised, under standardized rearing and feeding conditions, in the sib-testing program of the Goland C21 boar line and slaughtered at nine months of age and average body weight of 166 ± 15 kg. Accuracy was generally relatively poor, but R² in external validation was >0.7 for iodine number and concentration of C18:2n-6, polyunsaturated FAs and omega-6 FAs. To assess the effectiveness of the on-site predictions as indicator traits in a breeding program, (co)variance components of the measured traits (OBS) and of their predictions using in-lab (in-lab-PR) or on-site (on-site-PR) spectrometers were estimated. Available records for OBS were 6814 and 2048, for iodine number and FA composition, respectively. Predictions using in-lab were available for pigs slaughtered between 2006 and 2014, for a total of 10 153 records. Predictions using on-site were obtained from spectra collected since 2011, for a total of 10 296 records. The estimated heritabilities for the investigated traits ranged from 0.34 to 0.50 and were greater for on-site-PR than for OBS. Genetic correlations between OBS and in-lab-PR were very close to 1.00 for all the investigated traits, whereas those between OBS and on-site-PRED ranged from 0.86 to 0.94. On-site visible-IR predictions are accurate enough to support the use of this technique for large-scale phenotyping of raw ham fat quality, even when dealing with animals of a single genetic line raised in standardized conditions, and may be implemented as indicator traits in breeding programs.

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.

Prediction of Milk Coagulation Properties and Individual Cheese Yield in Sheep Using Partial Least Squares Regression

Simple Summary

Considered that all sheep milk in Italy is destined for cheese processing, traits describing rennet coagulation aptitude should be among the most important selection goals for dairy breeds. To reduce the costs and logistics related to the large-scale recording of these traits, mid-infrared (MIR) spectroscopy could be conveniently used to generate reliable predictions without any additional cost. The aims of this research were to predict the milk coagulation properties (MCP) and individual cheese yield (ILCY) in sheep by MIR spectrometry using partial least squares regression (PLS), and to compare different data pre-treatment procedures. The prediction results observed in the present study, although moderate, suggest the possibility of adding novel phenotypes (e.g., MCP and ILCY) in breeding schemes for dairy sheep breeds. Mid-infrared spectroscopy coupled with PLS regression could allow the prediction of phenotypes at the population level without additional costs.

Abstract

The objectives of this study were (i) the prediction of sheep milk coagulation properties (MCP) and individual laboratory cheese yield (ILCY) from mid-infrared (MIR) spectra by using partial least squares (PLS) regression, and (ii) the comparison of different data pre-treatments on prediction accuracy. Individual milk samples of 970 Sarda breed ewes were analyzed for rennet coagulation time (RCT), curd-firming time (k20), and curd firmness (a30) using the Formagraph instrument; ILCY was measured by micro-manufacturing assays. An Furier-transform Infrared (FTIR) milk-analyzer was used for the estimation of the milk gross composition and the recording of MIR spectrum. The dataset (n = 859, after the exclusion of 111 noncoagulating samples) was divided into two sub-datasets: the data of 700 ewes were used to estimate prediction model parameters, and the data of 159 ewes were used to validate the model. Four prediction scenarios were compared in the validation, differing for the use of whole or reduced MIR spectrum and the use of raw or corrected data (locally weighted scatterplot smoothing). PLS prediction statistics were moderate. The use of the reduced MIR spectrum yielded the best results for the considered traits, whereas the data correction improved the prediction ability only when the whole MIR spectrum was used. In conclusion, PLS achieves good accuracy of prediction, in particular for ILCY and RCT, and it may enable increasing the number of traits to be included in breeding programs for dairy sheep without additional costs and logistics.

Prediction of blood β-hydroxybutyrate content and occurrence of hyperketonemia in early-lactation, pasture-grazed dairy cows using milk infrared spectra

The objective of this study was to evaluate the ability of milk infrared spectra to predict blood β-hydroxybutyrate (BHB) concentration for use as a management tool for cow metabolic health on pasture-grazed dairy farms and for large-scale phenotyping for genetic evaluation purposes. The study involved 542 cows (Holstein-Friesian and Holstein-Friesian × Jersey crossbreds), from 2 farms located in the Waikato and Taranaki regions of New Zealand that operated under a seasonal-calving, pasture-based dairy system. Milk infrared spectra were collected once a week during the first 5 wk of lactation. A blood "prick" sample was taken from the ventral labial vein of each cow 3 times a week for the first 5 wk of lactation. The content of BHB in blood was measured immediately using a handheld device. After outlier elimination, 1,910 spectra records and corresponding BHB measures were used for prediction model development. Partial least square regression and partial least squares discriminant analysis were used to develop prediction models for quantitative determination of blood BHB content and for identifying cows with hyperketonemia (HYK). Both quantitative and discriminant predictions were developed using the phenotypes and infrared spectra from two-thirds of the cows (randomly assigned to the calibration set) and tested using the remaining one-third (validation set). A moderate accuracy was obtained for prediction of blood BHB. The coefficient of determination (R²) of the prediction model in calibration was 0.56, with a root mean squared error of prediction of 0.28 mmol/L and a ratio of performance to deviation, calculated as the ratio of the standard deviation of the partial least squares model calibration set to the standard error of prediction, of 1.50. In the validation set, the R² was 0.50, with root mean squared error of prediction values of 0.32 mmol/L, which resulted in a ratio of performance to deviation of 1.39. When the reference test for HYK was defined as blood concentration of BHB ≥1.2 mmol/L, discriminant models indicated that milk infrared spectra correctly classified 76% of the HYK-positive cows and 82% of the HYK-negative cows. The quantitative models were not able to provide accurate estimates, but they could differentiate between high and low BHB concentrations. Furthermore, the discriminant models allowed the classification of cows with reasonable accuracy. This study indicates that the prediction of blood BHB content or occurrence of HYK from milk spectra is possible with moderate accuracy in pasture-grazed cows and could be used during routine milk testing. Applicability of infrared spectroscopy is not likely suited for obtaining accurate BHB measurements at an individual cow level, but discriminant models might be used in the future as herd-level management tools for classification of cows that are at risk of HYK, whereas quantitative models might provide large-scale phenotypes to be used as an indicator trait for breeding cows with improved metabolic health.

Strategies for noise reduction and standardization of milk mid-infrared spectra from dairy cattle

The use of Fourier-transform mid-infrared (FTIR) spectroscopy is of interest to the dairy industry worldwide for predicting milk composition and other novel traits that are difficult or expensive to measure directly. Although there are many valuable applications for FTIR spectra, noise from differences in spectral responses between instruments is problematic because it reduces prediction accuracy if ignored. The purpose of this study was to develop strategies to reduce the impact of noise and to compare methods for standardizing FTIR spectra in order to reduce between-instrument variability in multiple-instrument networks. Noise levels in bands of the infrared spectrum caused by the water content of milk were characterized, and a method for identifying and removing outliers was developed. Two standardization methods were assessed and compared: piecewise direct standardization (PDS), which related spectra on a primary instrument to spectra on 5 other (secondary) instruments using identical milk-based reference samples (n = 918) analyzed across the 6 instruments; and retroactive percentile standardization (RPS), whereby percentiles of observed spectra from routine milk test samples (n = 2,044,094) were used to map and exploit primary- and secondary-instrument relationships. Different applications of each method were studied to determine the optimal way to implement each method across time. Industry-standard predictions of milk components from 2,044,094 spectra records were regressed against predictions from spectra before and after standardization using PDS or RPS. The PDS approach resulted in an overall decrease in root mean square error between industry-standard predictions and predictions from spectra from 0.190 to 0.071 g/100 mL for fat, from 0.129 to 0.055 g/100 mL for protein, and from 0.143 to 0.088 g/100 mL for lactose. Reductions in prediction error for RPS were similar but less consistent than those for PDS across time, but similar reductions were achieved when PDS coefficients were updated monthly and separate primary instruments were assigned for the North and South Islands of New Zealand. We demonstrated that the PDS approach is the most consistent method to reduce prediction errors across time. We also showed that the RPS approach is sensitive to shifts in milk composition but can be used to reduce prediction errors, provided that secondary-instrument spectra are standardized to a primary instrument with samples of broadly equivalent milk composition. Appropriate implementation of either of these approaches will improve the quality of predictions based on FTIR spectra for various downstream applications.

Variation of Detailed Protein Composition of Cow Milk Predicted from a Large Database of Mid-Infrared Spectra

Simple Summary

Milk proteins are one of the most valuable milk components. The objective of the present study was to assess sources of variation of detailed protein composition predicted from infrared spectra in milk of dairy and dual-purpose cattle breeds. Results showed that protein fractions were primarily influenced by days in milk, and the relative proportion of each fraction through lactation was not constant. Protein fractions correlated with crude protein, total casein, fat and milk urea nitrogen. In perspective, mid-infrared predictions of milk fractions could be useful for the dairy sector to improve nutritional and technological properties of milk.

Abstract

This study aimed to investigate factors affecting protein fractions, namely α-casein (α-CN), β-casein (β-CN), κ-casein (κ-CN), β-lactoglobulin (β-LG) and α-lactalbumin (α-LA) predicted from milk infrared spectra in milk of dairy and dual-purpose cattle breeds. The dataset comprised 735,328 observations from 49,049 cows in 1782 herds. Results highlighted significant differences of protein fractions in milk of the studied breeds. Significant variations of protein fractions were found also through parities and lactation, with the latter thoroughly influencing protein fractions percentage. Interesting correlations (r) were estimated between β-CN, κ-CN and β-LG, expressed as percentage of crude protein, and milk urea nitrogen (r = 0.31, −0.20 and −0.26, respectively) and between α-LA and fat percentage (r = 0.41). The present study paves the way for future studies on the associations between protein fractions and milk technological properties, and for the estimation of genetic parameters of predicted protein composition.

Standardization of milk mid-infrared spectrometers for the transfer and use of multiple models

An increasing number of models are being developed to provide information from milk Fourier transform mid-infrared (FT-MIR) spectra on fine milk composition, technological properties of milk, or even cows' physiological status. In this context, and to take advantage of these existing models, the purpose of this work was to evaluate whether a spectral standardization method can enable the use of multiple equations within a network of different FT-MIR spectrometers. The piecewise direct standardization method was used, matching "slave" instruments to a common reference, the "master." The effect of standardization on network reproducibility was assessed on 66 instruments from 3 different brands by comparing the spectral variability of the slaves and the master with and without standardization. With standardization, the global Mahalanobis distance from the slave spectra to the master spectra was reduced on average from 2,655.9 to 14.3, representing a significant reduction of noninformative spectral variability. The transfer of models from instrument to instrument was tested using 3 FT-MIR models predicting (1) the quantity of daily methane emitted by dairy cows, (2) the concentration of polyunsaturated fatty acids in milk, and (3) the fresh cheese yield. The differences, in terms of root mean squared error, between master predictions and slave predictions were reduced after standardization on average from 103 to 17 g/d, from 0.0315 to 0.0045 g/100 mL of milk, and from 2.55 to 0.49 g of curd/100 g of milk, respectively. For all the models, standard deviations of predictions among all the instruments were also reduced by 5.11 times for methane, 5.01 times for polyunsaturated fatty acids, and 7.05 times for fresh cheese yield, showing an improvement of prediction reproducibility within the network. Regarding the results obtained, spectral standardization allows the transfer and use of multiple models on all instruments as well as the improvement of spectral and prediction reproducibility within the network. The method makes the models universal, thereby offering opportunities for data exchange and the creation and use of common robust models at an international level to provide more information to the dairy sector from direct analysis of milk.

Comparison of Bayesian regression models and partial least squares regression for the development of infrared prediction equations

The objective of this study was to compare the prediction accuracy of 92 infrared prediction equations obtained by different statistical approaches. The predicted traits included fatty acid composition (n = 1,040); detailed protein composition (n = 1,137); lactoferrin (n = 558); pH and coagulation properties (n = 1,296); curd yield and composition obtained by a micro-cheese making procedure (n = 1,177); and Ca, P, Mg, and K contents (n = 689). The statistical methods used to develop the prediction equations were partial least squares regression (PLSR), Bayesian ridge regression, Bayes A, Bayes B, Bayes C, and Bayesian least absolute shrinkage and selection operator. Model performances were assessed, for each trait and model, in training and validation sets over 10 replicates. In validation sets, Bayesian regression models performed significantly better than PLSR for the prediction of 33 out of 92 traits, especially fatty acids, whereas they yielded a significantly lower prediction accuracy than PLSR in the prediction of 8 traits: the percentage of C18:1n-7 trans-9 in fat; the content of unglycosylated κ-casein and its percentage in protein; the content of α-lactalbumin; the percentage of α_S2-casein in protein; and the contents of Ca, P, and Mg. Even though Bayesian methods produced a significant enhancement of model accuracy in many traits compared with PLSR, most variations in the coefficient of determination in validation sets were smaller than 1 percentage point. Over traits, the highest predictive ability was obtained by Bayes C even though most of the significant differences in accuracy between Bayesian regression models were negligible.