Identification and quantification of goat milk adulteration using mid-infrared spectroscopy and chemometrics

The fraudulent adulteration of goat milk with cheaper and more available milk of other species such as cow milk is occurrence. The aims of the present study were to investigate the effect of goat milk adulteration with cow milk on the mid-infrared (MIR) spectrum and further evaluate the potential of MIR spectroscopy to identify and quantify the goat milk adulterated. Goat milk was adulterated with cow milk at 5 different levels including 10%, 20%, 30%, 40%, and 50%. Statistical analysis showed that the adulteration had significant effect on the majority of the spectral wavenumbers. Then, the spectrum was preprocessed with standard normal variate (SNV), multiplicative scattering correction (MSC), Savitzky-Golay smoothing (SG), SG plus SNV, and SG plus MSC, and partial least squares discriminant analysis (PLS-DA) and partial least squares regression (PLSR) were used to establish classification and regression models, respectively. PLS-DA models obtained good results with all the sensitivity and specificity over 0.96 in the cross-validation set. Regression models using raw spectrum obtained the best result, with coefficient of determination (R2), root mean square error (RMSE), and the ratio of performance to deviation (RPD) of cross-validation set were 0.98, 2.01, and 8.49, respectively. The results preliminarily indicate that the MIR spectroscopy is an effective technique to detect the goat milk adulteration with cow milk. In future, milk samples from different origins and different breeds of goats and cows should be collected, and more sophisticated adulteration at low levels should be further studied to explore the potential and effectiveness of milk mid-infrared spectroscopy and chemometrics.

Genetic parameters of milk mid-infrared spectra and their genetic relationships with milk production and feed efficiency traits in French Lacaune dairy sheep

In French dairy sheep, Fourier transform infrared (FTIR) milk spectral data routinely predict the major milk components used in national genetic evaluations. The direct influence of genetic and environmental factors on milk FTIR spectra has been widely studied in dairy cattle, and relatively little in dairy ewes. In this study, 36,873 milk test-day records were available for 4,712 French Lacaune ewes farmed on 8 commercial farms. Our main goals were to provide the first description of spectral data and estimate the genetic parameters of French Lacaune dairy sheep during lactation. Principal component analysis (PCA) results demonstrated the impact of the lactation period on specific wavenumbers, allowing the identification of FTIR spectra collected at early (mo 2-4) and late (mo 5-7) lactation stages. The average estimated heritability (±mean SE) of the FTIR milk spectra from 2,971 to 926 cm-1 (446 wavenumbers) was 0.29 ± 0.02, ranging from 0.13 ± 0.01 to 0.42 ± 0.02. Furthermore, the heritabilities of spectra collected at the beginning or end of lactation changed at each point of the spectrum. However, at each wavenumber, the genomic correlation of transmittance values between these 2 lactation periods was high (>0.77), indicating the absence of a genotype-environment interaction. The genomic correlations between spectral regions and milk production traits (i.e., daily milk yield, fat and protein content, somatic cell score) varied from moderate to high. The results suggested that the most heritable areas of the spectrum were also genetically associated with dairy traits. Finally, the genomic correlations observed between the ewes' feed efficiency traits and the FTIR spectrum were moderate to high, while the genomic correlations between the change in body condition score and spectral data were rather low to moderate. This study confirmed that spectral data from Lacaune ewe milk were heritable, evolved phenotypically and genetically during lactation and were genetically correlated with traits included in breeding goals or traits of interest to the dairy industry.

Spectral Profiling (Fourier Transform Infrared Spectroscopy) and Machine Learning for the Recognition of Milk from Different Bovine Breeds

The Podolica cattle breed is widespread in southern Italy, and its productivity is characterized by low yields and an extraordinary quality of milk and meats. Most of the milk produced is transformed into "Caciocavallo Podolico" cheese, which is made with 100% Podolica milk. Fourier Transform Infrared Spectroscopy (FTIR) is the technique that, in this research work, was applied together with machine learning to discriminate 100% Podolica milk from contamination of other Calabrian cattle breeds. The analysis on the test set produced a misclassification percentage of 6.7%. Among the 15 non-Podolica samples in the test set, 2 were misclassified and recognized as Podolica milk even though the milk was from other species. The correct classification rate improved to 100% when the same method was applied to the recognition of Podolica and Pezzata Rossa milk produced by the same farm. Furthermore, this technique was tested for the recognition of Podolica milk mixed with milk from other bovine species. The multivariate model and the respective confusion matrices obtained showed that all the 14 Podolica samples (test set) mixed with 40% non-Podolica milk were correctly classified. In addition, Pezzata Rossa milk produced by the same farm was detected as a contaminant in Podolica milk from the same farm down to concentrations as little as 5% with a 100% correct classification rate in the test set. The method described yielded higher accuracy values when applied to the discrimination of milks from different breeds belonging to the same farm. One of the reasons for this phenomenon could be linked to the elimination of the environmental variable. However, the results obtained in this work demonstrate the possibility of using FTIR to discriminate between milks from different breeds.

The use of milk Fourier-transform infrared spectra for predicting cheesemaking traits in Grana Padano Protected Designation of Origin cheese

The prediction of the cheese yield (%CY) traits for curd, solids, and retained water and the amount of fat, protein, solids, and energy recovered from the milk into the curd (%REC) by Bayesian models, using Fourier-transform infrared spectroscopy (FTIR), can be of significant economic interest to the dairy industry and can contribute to the improvement of the cheese process efficiency. The yields give a quantitative measure of the ratio between weights of the input and output of the process, whereas the nutrient recovery allows to assess the quantitative transfer of a component from milk to cheese (expressed in % of the initial weight). The aims of this study were: (1) to investigate the feasibility of using bulk milk spectra to predict %CY and %REC traits, and (2) to quantify the effect of the dairy industry and the contribution of single-spectrum wavelengths on the prediction accuracy of these traits using vat milk samples destined to the production of Grana Padano Protected Designation of Origin cheese. Information from 72 cheesemaking days (in total, 216 vats) from 3 dairy industries were collected. For each vat, the milk was weighed and analyzed for composition (total solids [TS], lactose, protein, and fat). After 48 h from cheesemaking, each cheese was weighed, and the resulting whey was sampled for composition as well (TS, lactose, protein, and fat). Two spectra from each milk sample were collected in the range between 5,011 and 925 cm-1 and averaged before the data analysis. The calibration models were developed via a Bayesian approach by using the BGLR (Bayesian Generalized Linear Regression) package of R software. The performance of the models was assessed by the coefficient of determination (R2VAL) and the root mean squared error (RMSEVAL) of validation. Random cross-validation (CVL) was applied [80% calibration and 20% validation set] with 10 replicates. Then, a stratified cross-validation (SCV) was performed to assess the effect of the dairy industry on prediction accuracy. The study was repeated using a selection of informative wavelengths to assess the necessity of using whole spectra to optimize prediction accuracy. Results showed the feasibility of using FTIR spectra and Bayesian models to predict cheesemaking traits. The R2VAL values obtained with the CVL procedure were promising in particular for the %CY and %REC for protein, ranging from 0.44 to 0.66 with very low RMSEVAL (from 0.16 to 0.53). Prediction accuracy obtained with the SCV was strongly influenced by the dairy factory industry. The general low values gained with the SCV do not permit a practical application of this approach, but they highlight the importance of building calibration models with a dataset covering the largest possible sample variability. This study also demonstrated that the use of the full FTIR spectra may be redundant for the prediction of the cheesemaking traits and that a specific selection of the most informative wavelengths led to improved prediction accuracy. This could lead to the development of dedicated spectrometers using selected wavelengths with built-in calibrations for the online prediction of these innovative traits.

Assessing the relationship between somatic cell count and the milk mid-infrared spectrum in Chinese Holstein cows

BACKGROUND

Milk produced by dairy cows is a complex combination of many components, but the effect of mastitis has only been investigated for a few of these components. Milk mid-infrared (MIR) spectra can reflect the global composition of milk, and this study aimed to detect the relationships between milk MIR spectral wavenumbers and milk somatic cell count (SCC)-a sensitive biomarker for mastitis.

METHODS

Pearson correlation analysis was used to calculate the correlation coefficient between somatic count score (SCS) and spectral wavenumbers. A general linear mixed model was applied to investigate the effect of three different classes of SCC (low, middle and high) on spectral wavenumbers.

RESULTS

The mean correlation coefficient between the 'fingerprint region' (wavenumbers 925-1582 cm-1 ) and the SCS was higher than that for other regions of the MIR spectrum, and the specific wavenumber with the strongest correlation with the SCS was within the 'fingerprint region'. SCC class had a significant (p < 0.05) effect on 639 spectral wavenumbers. In particular, some spectral wavenumbers within the 'fingerprint region' were highly affected by the SCC class.

LIMITATION

The data were collected from only one province in China, so the generalisability of the findings may be limited.

CONCLUSION

SCC had close relationships with milk spectral wavenumbers related to important milk components or chemical bonds.

Integrating on-farm and genomic information improves the predictive ability of milk infrared prediction of blood indicators of metabolic disorders in dairy cows

BACKGROUND

Blood metabolic profiles can be used to assess metabolic disorders and to evaluate the health status of dairy cows. Given that these analyses are time-consuming, expensive, and stressful for the cows, there has been increased interest in Fourier transform infrared (FTIR) spectroscopy of milk samples as a rapid, cost-effective alternative for predicting metabolic disturbances. The integration of FTIR data with other layers of information such as genomic and on-farm data (days in milk (DIM) and parity) has been proposed to further enhance the predictive ability of statistical methods. Here, we developed a phenotype prediction approach for a panel of blood metabolites based on a combination of milk FTIR data, on-farm data, and genomic information recorded on 1150 Holstein cows, using BayesB and gradient boosting machine (GBM) models, with tenfold, batch-out and herd-out cross-validation (CV) scenarios.

RESULTS

The predictive ability of these approaches was measured by the coefficient of determination (R²). The results show that, compared to the model that includes only FTIR data, integration of both on-farm (DIM and parity) and genomic information with FTIR data improves the R² for blood metabolites across the three CV scenarios, especially with the herd-out CV: R² values ranged from 5.9 to 17.8% for BayesB, from 8.2 to 16.9% for GBM with the tenfold random CV, from 3.8 to 13.5% for BayesB and from 8.6 to 17.5% for GBM with the batch-out CV, and from 8.4 to 23.0% for BayesB and from 8.1 to 23.8% for GBM with the herd-out CV. Overall, with the model that includes the three sources of data, GBM was more accurate than BayesB with accuracies across the CV scenarios increasing by 7.1% for energy-related metabolites, 10.7% for liver function/hepatic damage, 9.6% for oxidative stress, 6.1% for inflammation/innate immunity, and 11.4% for mineral indicators.

CONCLUSIONS

Our results show that, compared to using only milk FTIR data, a model integrating milk FTIR spectra with on-farm and genomic information improves the prediction of blood metabolic traits in Holstein cattle and that GBM is more accurate in predicting blood metabolites than BayesB, especially for the batch-out CV and herd-out CV scenarios.

Importance of Mid-Infrared Spectra Regions for the Prediction of Mastitis and Ketosis in Dairy Cows

Mid-infrared (MIR) spectroscopy is routinely applied to determine major milk components, such as fat and protein. Moreover, it is used to predict fine milk composition and various traits pertinent to animal health. MIR spectra indicate an absorbance value of infrared light at 1060 specific wavenumbers from 926 to 5010 cm−1. According to research, certain parts of the spectrum do not contain sufficient information on traits of dairy cows. Hence, the objective of the present study was to identify specific regions of the MIR spectra of particular importance for the prediction of mastitis and ketosis, performing variable selection analysis. Partial least squares discriminant analysis (PLS-DA) along with three other statistical methods, support vector machine (SVM), least absolute shrinkage and selection operator (LASSO), and random forest (RF), were compared. Data originated from the Austrian milk recording and associated health monitoring system (GMON). Test-day data and corresponding MIR spectra were linked to respective clinical mastitis and ketosis diagnoses. Certain wavenumbers were identified as particularly relevant for the prediction models of clinical mastitis (23) and ketosis (61). Wavenumbers varied across four distinct statistical methods as well as concerning different traits. The results indicate that variable selection analysis could potentially be beneficial in the process of modeling.

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.

Phenotypic Analysis of Fourier-Transform Infrared Milk Spectra in Dairy Goats

The infrared spectrum of bovine milk is used to predict many interesting traits, whereas there have been few studies on goat milk in this regard. The objective of this study was to characterize the major sources of variation in the absorbance of the infrared spectrum in caprine milk samples. A total of 657 goats belonging to 6 breeds and reared on 20 farms under traditional and modern dairy systems were milk-sampled once. Fourier-transform infrared (FTIR) spectra were taken (2 replicates per sample, 1314 spectra), and each spectrum contained absorbance values at 1060 different wavenumbers (5000 to 930 × cm^-1), which were treated as a response variable and analyzed one at a time (i.e., 1060 runs). A mixed model, including the random effects of sample/goat, breed, flock, parity, stage of lactation, and the residual, was used. The pattern and variability of the FTIR spectrum of caprine milk was similar to those of bovine milk. The major sources of variation in the entire spectrum were as follows: sample/goat (33% of the total variance); flock (21%); breed (15%); lactation stage (11%); parity (9%); and the residual unexplained variation (10%). The entire spectrum was segmented into five relatively homogeneous regions. Two of them exhibited very large variations, especially the residual variation. These regions are known to be affected by the absorbance of water, although they also exhibited wide variations in the other sources of variation. The average repeatability of these two regions were 45% and 75%, whereas for the other three regions it was about 99%. The FTIR spectrum of caprine milk could probably be used to predict several traits and to authenticate the origin of goat milk.

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.

Prediction of fatty acid composition using milk spectral data and its associations with various mid-infrared spectral regions in Michigan Holsteins

Fatty acid composition in milk is not only reflective of nutritional quality but also potentially predictive of other attributes (e. g. including the cow's energy balance and its relative output of methane emissions). Furthermore, a higher ratio of long-chain to short-chain fatty acids or mean carbon number has been associated with negative energy balance in dairy cows, whereas enhanced nutritional properties have been generally associated with higher levels of unsaturation. We set out to directly compare Bayesian regression strategies with partial least squares for the prediction of various milk fatty acids using Fourier-transform infrared spectrum data on 777 milk samples taken from 579 cows on 4 Michigan dairy herds between 5 and 90 d in milk. We also set out to identify those spectral regions that might be associated with fatty acids and whether carbon number or level of unsaturation might contribute to the strength of these associations. These associations were based on adaptively clustered windows of wavenumbers to mitigate the distorting effects of severe multicollinearity on marginal associations involving individual wavenumbers. In general, Bayesian regression methods, particularly the variable selection method BayesB, outperformed partial least squares regression for cross-validation prediction accuracy for both individual fatty acids and fatty acid groups. Strong signals for wavenumber associations using BayesB were well distributed throughout the mid-infrared spectrum, particularly between 910 and 3,998 cm^-1. Carbon number appeared to be linearly related to strength of wavenumber associations for 38 moderately to highly predicted fatty acids within the spectral regions of 2,286 to 2,376 and 2,984 to 3,100 cm^-1, whereas nonlinear associations were determined within 1,141 to 1,205; 1,570 to 1,630; and 1,727 to 1,768 cm^-1. However, no such associations were detected with level of unsaturation. Spectral regions where there were significant relationships between strength of association and carbon number may be useful targets for inferring the relative proportion of long-chain to short-chain fatty acids, and hence energy balance.

Phenotypic and genetic variation of ultraviolet-visible-infrared spectral wavelengths of bovine meat

Spectroscopic predictions can be used for the genetic improvement of meat quality traits in cattle. No information is however available on the genetics of meat absorbance spectra. This research investigated the phenotypic variation and the heritability of meat absorbance spectra at individual wavelengths in the ultraviolet-visible and near-infrared region (UV-Vis-NIR) obtained with portable spectrometers. Five spectra per instrument were taken on the ribeye surface of 1185 Piemontese young bulls from 93 farms (13,182 Herd-Book pedigree relatives). Linear animal model analyses of 1481 single-wavelengths from UV-Vis-NIRS and 125 from Micro-NIRS were carried out separately. In the overlapping regions, the proportions of phenotypic variance explained by batch/date of slaughter (14 ± 6% and 17 ± 7%,), rearing farm (6 ± 2% and 5 ± 3%), and the residual variances (72 ± 10% and 72 ± 5%) were similar for the UV-Vis-NIRS and Micro-NIRS, but additive genetics (7 ± 2% and 4 ± 2%) and heritability (8.3 ± 2.3% vs 5.1 ± 0.6%) were greater with the Micro-NIRS. Heritability was much greater for the visible fraction (25.2 ± 11.4%), especially the violet, blue and green colors, than for the NIR fraction (5.0 ± 8.0%). These results allow a better understanding of the possibility of using the absorbance of visible and infrared wavelengths correlated with meat quality traits for the genetic improvement in beef cattle.

Influence of Estrus on the Milk Characteristics and Mid-Infrared Spectra of Dairy Cows

Simple Summary

Some studies have confirmed the variation in milk profiles when dairy cows show estrus. However, only a few milk components, such as fat, protein, and lactose, have been investigated so far, and thus any changes in the many other parts of milk’s composition due to estrus are unknown. Milk mid-infrared (MIR) spectra consist of wavenumbers, which provide insight into the chemical composition of milk. The MIR spectrum reflects the global composition of milk, but this information is currently underused. In this study, we considered MIR wavenumbers as traits, and directly studied the spectral information as a way to study the estrus of dairy cows linked to milk composition. This research provides a deeper understanding of the milk MIR spectrum and may lead to new approaches for estrus detection in dairy cows from routine milk analysis, thereby guiding an opportune insemination time.

Abstract

Milk produced by dairy cows is a complex combination of many components. However, at present, changes in only a few milk components (e.g., fat, protein, and lactose) during the estrus cycle in dairy cows have been documented. Mid-infrared (MIR) spectroscopy is a worldwide method routinely used for milk analysis, as MIR spectra reflect the global composition of milk. Therefore, this study aimed to investigate the changes in milk MIR spectra and milk production traits (fat, protein, lactose, urea, total solids (TS), and solid not fat (SnF)) due to estrus. Cows that were successfully inseminated, leading to conception, were included. Cows confirmed to be pregnant were considered to be in estrus at the day of insemination (day 0). A general linear mixed model, which included the random effect of cows, the fixed classification effects of parity number, days in relation to estrus, as well as the interaction between parity number and days in relation to estrus, was applied to investigate the changes in milk production traits and 1060 milk infrared wavenumbers, ranging from 925 to 5011 cm⁻¹, of 371 records from 162 Holstein cows on the days before (day −3, day −2, and day −1) and on the day of estrus (day 0). The days in relation to estrus had a significant effect on fat, protein, urea, TS, and SnF, whose contents increased from day −3 to day 0. Lactose did not seem to be significantly influenced by the occurrence of estrus. The days in relation to estrus had significant effects on the majority of the wavenumbers. Besides, we found that some of the wavenumbers in the water absorption regions were significantly changed on the days before and on the day of estrus. This suggests that these wavenumbers may contain useful information. In conclusion, the changes in the milk composition due to estrus can be observed through the analysis of the milk MIR spectrum. Further analyses are warranted to more deeply explore the potential use of milk MIR spectra in the detection of estrus.

Breed of goat affects the prediction accuracy of milk coagulation properties using Fourier-transform infrared spectroscopy

The prediction of traditional goat milk coagulation properties (MCP) and curd firmness over time (CF_t) parameters via Fourier-transform infrared (FTIR) spectroscopy can be of significant economic interest to the dairy industry and can contribute to the breeding objectives for the genetic improvement of dairy goat breeds. Therefore, the aims of this study were to (1) explore the variability of milk FTIR spectra from 4 goat breeds (Camosciata delle Alpi, Murciano-Granadina, Maltese, and Sarda), and to assess the possible discriminant power of milk FTIR spectra among breeds, (2) assess the viability to predict coagulation traits by using milk FTIR spectra, and (3) quantify the effect of the breed on the prediction accuracy of MCP and CF_t parameters. In total, 611 individual goat milk samples were used. Analysis of variance of measured MCP and CF_t parameters was carried out using a mixed model including the farm and pendulum as random factors, and breed, parity, and days in milk as fixed factors. Milk spectra for each goat were collected over the spectral range from wavenumber 5,011 to 925 × cm^-1. Discriminant analysis of principal components was used to assess the ability of FTIR spectra to identify breed of origin. A Bayesian model was used to calibrate equations for each coagulation trait. The accuracy of the model and the prediction equation was assessed by cross-validation (CRV; 80% training and 20% testing set) and stratified CRV (SCV; 3 breeds in the training set, one breed in the testing set) procedures. Prediction accuracy was assessed by using coefficient of determination of validation (R²_VAL), the root mean square error of validation (RMSE_VAL), and the ratio performance deviation. Moreover, measured and FTIR predicted traits were compared in the SCV procedure by assessing their least squares means for the breed effect, Pearson correlations, and variance heteroscedasticity. Results showed the feasibility of using FTIR spectra and multivariate analyses to correctly assign milk samples to their breeds of origin. The R²_VAL values obtained with the CRV procedure were moderate to high for the majority of coagulation traits, with RMSE_VAL and ratio performance deviation values increasing as the coagulation process progresses from rennet addition. Prediction accuracy obtained with the SCV were strongly influenced by the breed, presenting general low values restricting a practical application. In addition, the low Pearson correlation coefficients of Sarda breed for all the traits analyzed, and the heteroscedastic variances of Camosciata delle Alpi, Murciano-Granadina, and Maltese breeds, further indicated that it is fundamental to consider the differences existing among breeds for the prediction of milk coagulation traits.

The heritable landscape of near-infrared and Raman spectroscopic measurements to improve lipid content in Atlantic salmon fillets

Background

Product quality and production efficiency of Atlantic salmon are, to a large extent, influenced by the deposition and depletion of lipid reserves. Fillet lipid content is a heritable trait and is unfavourably correlated with growth, thus genetic management of fillet lipid content is needed for sustained genetic progress in these two traits. The laboratory-based reference method for recording fillet lipid content is highly accurate and precise but, at the same time, expensive, time-consuming, and destructive. Here, we test the use of rapid and cheaper vibrational spectroscopy methods, namely near-infrared (NIR) and Raman spectroscopy both as individual phenotypes and phenotypic predictors of lipid content in Atlantic salmon.

Results

Remarkably, 827 of the 1500 individual Raman variables (i.e. Raman shifts) of the Raman spectrum were significantly heritable (heritability (h²) ranging from 0.15 to 0.65). Similarly, 407 of the 2696 NIR spectral landscape variables (i.e. wavelengths) were significantly heritable (h² = 0.27–0.40). Both Raman and NIR spectral landscapes had significantly heritable regions, which are also informative in spectroscopic predictions of lipid content. Partial least square predicted lipid content using Raman and NIR spectra were highly concordant and highly genetically correlated with the lipid content values (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${r}_{\text{g}}$$\end{document}rg = 0.91–0.98) obtained with the reference method using Lin’s concordance correlation coefficient (CCC = 0.63–0.90), and were significantly heritable (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h}^{2}$$\end{document}h2 = 0.52–0.67).

Conclusions

Both NIR and Raman spectral landscapes show substantial additive genetic variation and are highly genetically correlated with the reference method. These findings lay down the foundation for rapid spectroscopic measurement of lipid content in salmonid breeding programmes.

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.

Predicting milk mid-infrared spectra from first-parity Holstein cows using a test-day mixed model with the perspective of herd management

The use of test-day models to model milk mid-infrared (MIR) spectra for genetic purposes has already been explored; however, little attention has been given to their use to predict milk MIR spectra for management purposes. The aim of this paper was to study the ability of a test-day mixed model to predict milk MIR spectra for management purposes. A data set containing 467,496 test-day observations from 53,781 Holstein dairy cows in first lactation was used for model building. Principal component analysis was implemented on the selected 311 MIR spectral wavenumbers to reduce the number of traits for modeling; 12 principal components (PC) were retained, explaining approximately 96% of the total spectral variation. Each of the retained PC was modeled using a single trait test-day mixed model. The model solutions were used to compute the predicted scores of each PC, followed by a back-transformation to obtain the 311 predicted MIR spectral wavenumbers. Four new data sets, containing altogether 122,032 records, were used to test the ability of the model to predict milk MIR spectra in 4 distinct scenarios with different levels of information about the cows. The average correlation between observed and predicted values of each spectral wavenumber was 0.85 for the modeling data set and ranged from 0.36 to 0.62 for the scenarios. Correlations between milk fat, protein, and lactose contents predicted from the observed spectra and from the modeled spectra ranged from 0.83 to 0.89 for the modeling set and from 0.32 to 0.73 for the scenarios. Our results demonstrated a moderate but promising ability to predict milk MIR spectra using a test-day mixed model. Current and future MIR traits prediction equations could be applied on the modeled spectra to predict all MIR traits in different situations instead of developing one test-day model separately for each trait. Modeling MIR spectra would benefit farmers for cow and herd management, for instance through prediction of future records or comparison between observed and expected wavenumbers or MIR traits for the detection of health and management problems. Potential resulting tools could be incorporated into milk recording systems.

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.

Genome-wide association study on Fourier transform infrared milk spectra for two Danish dairy cattle breeds

Background

Infrared spectral analysis of milk is cheap, fast, and accurate. Infrared light interacts with chemical bonds present inside the milk, which means that Fourier transform infrared milk spectra are a reflection of the chemical composition of milk. Heritability of Fourier transform infrared milk spectra has been analysed previously. Further genetic analysis of Fourier transform infrared milk spectra could give us a better insight in the genes underlying milk composition. Breed influences milk composition, yet not much is known about the effect of breed on Fourier transform infrared milk spectra. Improved understanding of the effect of breed on Fourier transform infrared milk spectra could enhance efficient application of Fourier transform infrared milk spectra. The aim of this study is to perform a genome wide association study on a selection of wavenumbers for Danish Holstein and Danish Jersey. This will improve our understanding of the genetics underlying milk composition in these two dairy cattle breeds.

Results

For each breed separately, fifteen wavenumbers were analysed. Overall, more quantitative trait loci were observed for Danish Jersey compared to Danish Holstein. For both breeds, the majority of the wavenumbers was most strongly associated to a genomic region on BTA 14 harbouring DGAT1. Furthermore, for both breeds most quantitative trait loci were observed for wavenumbers that interact with the chemical bond C-O. For Danish Jersey, wavenumbers that interact with C-H were associated to genes that are involved in fatty acid synthesis, such as AGPAT3, AGPAT6, PPARGC1A, SREBF1, and FADS1. For wavenumbers which interact with –OH, associations were observed to genomic regions that have been linked to alpha-lactalbumin.

Conclusions

The current study identified many quantitative trait loci that underlie Fourier transform infrared milk spectra, and thus milk composition. Differences were observed between groups of wavenumbers that interact with different chemical bonds. Both overlapping and different QTL were observed for Danish Holstein and Danish Jersey.

Genetic Analysis of Milk Production Traits and Mid-Infrared Spectra in Chinese Holstein Population

Simple Summary

Usually, spectral data are used as predictors to predict milk components, animal characteristics, and even reproductive status. Another innovative way to use spectral data involves considering spectral wavenumbers as traits and then analyzing from the genetic perspective. In this study, we considered milk spectral data directly as traits, then detected the influence of some non-genetic factors on spectral wavenumbers and estimated the genetic parameters of spectral points. The result of the present study could be used as a management tool for dairy farm and also provides a further understanding of genetic background of milk mid-infrared (MIR) spectra. In future, milk spectral data could be applied more effective. For example, some sub-clinical diseases might be detected based on the difference between the expected and observed values of the spectral traits. In addition, we could also use genetic correlation between wavenumbers and a trait of interest, which are difficult and expensive to measure, to apply for the genetic improvement of dairy species.

Abstract

Milk composition always serves as an indicator for the cow’s health status and body condition. Some non-genetic factors such as parity, days in milk (DIM), and calving season, which obviously affect milk performance, therefore, need to be considered in dairy farm management. However, only a few milk compositions are used in the current animal selection programs. The mid-infrared (MIR) spectroscopy can reflect the global composition of milk, but this information is currently underused. The objectives of this study were to detect the effect of some non-genetic factors on milk production traits as well as 1060 individual spectral points covering from 925.92 cm⁻¹ to 5011.54 cm⁻¹, estimate heritabilities of milk production traits and MIR spectral wavenumbers, and explore the genetic correlations between milk production traits and 1060 individual spectral points in a Chinese Holstein population. The mixed models procedure of SAS software was used to test the non-genetic factors. Single-trait animal models were used to estimate heritabilities and bivariate animal models were used to estimate genetic correlations using the package of ASReml in R software. The results showed that herd, parity, calving season, and lactation stage had significant effects on the percentages of protein and lactose, whereas herd and lactation stage had significant effects on fat percentage. Moreover, the herd showed a significant effect on all of the 1060 individual wavenumbers, whereas lactation stage, parity, and calving season had significant effect on most of the wavenumbers of the lactose-region (925 cm⁻¹ to 1200 cm⁻¹), protein-region (1240 cm⁻¹ to 1600 cm⁻¹), and fat-regions (1680 cm⁻¹ to 1770 cm⁻¹ and 2800 cm⁻¹ to 3015 cm⁻¹). The estimated heritabilities for protein percentage (PP), fat percentage (FP), and lactose percentage (LP) were 0.08, 0.05, and 0.09, respectively. Further, the milk spectrum was heritable but low for most individual points. Heritabilities of 1060 individual spectral points were 0.04 on average, ranging from 0 to 0.11. In particular, heritabilities for wavenumbers of spectral regions related to water absorption were very low and even null, and heritabilities for wavenumbers of specific MIR regions associated with fat-I, fat-II, protein, and lactose were 0.04, 0.06, 0.05, and 0.06 on average, respectively. The genetic correlations between PP and FP, PP and LP, FP, and LP were 0.78, −0.29, and −0.14, respectively. In addition, PP, FP, and LP shared the similar patterns of genetic correlations with the spectral wavenumbers. The genetic correlations between milk production traits and spectral regions related to important milk components varied from weak to very strong (0.01 to 0.94, and −0.01 to −0.96). The current study could be used as a management tool for dairy farms and also provides a further understanding of the genetic background of milk MIR spectra.

Combined use of milk infrared spectra and genotypes can improve prediction of milk fat composition

It has been shown that milk infrared (IR) spectroscopy can be used to predict detailed milk fat composition. In addition, polymorphisms with substantial effects on milk fat composition have been identified. In this study, we investigated the combined use of milk IR spectroscopy and genotypes of dairy cows on the accuracy of predicting milk fat composition. Milk fat composition data based on gas chromatography and milk IR spectra were available for 1,456 Dutch Holstein Friesian cows. In addition, genotypes for the diacylglycerol acyltransferase 1 (DGAT1) K232A and stearoyl-CoA desaturase 1 (SCD1) A293V polymorphisms and a SNP located in an intron of the fatty acid synthase (FASN) gene were available. Adding SCD1 genotypes to the milk IR spectra resulted in a considerable improvement of the prediction accuracy for the unsaturated fatty acids C10:1, C12:1, C14:1 cis-9, and C16:1 cis-9 and their corresponding unsaturation indices. Adding DGAT1 genotypes to the milk IR spectra resulted in an improvement of the prediction accuracy for C16:1 cis-9 and C16 index. Adding genotypes of the FASN SNP to the IR spectra did not improve prediction of milk fat composition. This study demonstrated the potential of combining milk IR spectra with genotypic information from 3 polymorphisms to predict milk fat composition. We hypothesize that prediction accuracy of milk fat composition can be further improved by combining milk IR spectra with genomic breeding values.

The use of mid-infrared spectra to map genes affecting milk composition

The aim of this study was to investigate the feasibility of using mid-infrared (MIR) spectroscopy analysis of milk samples to increase the power and precision of genome-wide association studies (GWAS) for milk composition and to better distinguish linked quantitative trait loci (QTL). To achieve this goal, we analyzed phenotypic data of milk composition traits, related MIR spectra, and genotypic data comprising 626,777 SNP on 5,202 Holstein, Jersey, and crossbred cows. We performed a conventional GWAS on protein, lactose, fat, and fatty acid concentrations in milk, a GWAS on individual MIR wavenumbers, and a partial least squares regression (PLS), which is equivalent to a multi-trait GWAS, exploiting MIR data simultaneously to predict SNP genotypes. The PLS detected most of the QTL identified using single-trait GWAS, usually with a higher significance value, as well as previously undetected QTL for milk composition. Each QTL tends to have a different pattern of effects across the MIR spectrum and this explains the increased power. Because SNP tracking different QTL tend to have different patterns of effect, it was possible to distinguish closely linked QTL. Overall, the results of this study suggest that using MIR data through either GWAS or PLS analysis applied to genomic data can provide a powerful tool to distinguish milk composition QTL.

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.

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.

Multi-population GWAS and enrichment analyses reveal novel genomic regions and promising candidate genes underlying bovine milk fatty acid composition

Background

The power of genome-wide association studies (GWAS) is often limited by the sample size available for the analysis. Milk fatty acid (FA) traits are scarcely recorded due to expensive and time-consuming analytical techniques. Combining multi-population datasets can enhance the power of GWAS enabling detection of genomic region explaining medium to low proportions of the genetic variation. GWAS often detect broader genomic regions containing several positional candidate genes making it difficult to untangle the causative candidates. Post-GWAS analyses with data on pathways, ontology and tissue-specific gene expression status might allow prioritization among positional candidate genes.

Results

Multi-population GWAS for 16 FA traits quantified using gas chromatography (GC) in sample populations of the Chinese, Danish and Dutch Holstein with high-density (HD) genotypes detects 56 genomic regions significantly associated to at least one of the studied FAs; some of which have not been previously reported. Pathways and gene ontology (GO) analyses suggest promising candidate genes on the novel regions including OSBPL6 and AGPS on Bos taurus autosome (BTA) 2, PRLH on BTA 3, SLC51B on BTA 10, ABCG5/8 on BTA 11 and ALG5 on BTA 12. Novel genes in previously known regions, such as FABP4 on BTA 14, APOA1/5/7 on BTA 15 and MGST2 on BTA 17, are also linked to important FA metabolic processes.

Conclusion

Integration of multi-population GWAS and enrichment analyses enabled detection of several novel genomic regions, explaining relatively smaller fractions of the genetic variation, and revealed highly likely candidate genes underlying the effects. Detection of such regions and candidate genes will be crucial in understanding the complex genetic control of FA metabolism. The findings can also be used to augment genomic prediction models with regions collectively capturing most of the genetic variation in the milk FA traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5573-9) contains supplementary material, which is available to authorized users.

A landscape of the heritability of Fourier-transform infrared spectral wavelengths of milk samples by parity and lactation stage in Holstein cows

Fourier-transform near- and mid-infrared (FTIR) milk spectral data are routinely collected in many countries worldwide. Establishing an optimal strategy to use spectral data in genetic evaluations requires knowledge of the heritabilities of individual FTIR wavelength absorbances. Previous FTIR heritability estimates have been based on relatively small sample sizes and have not considered the possibility that heritability may vary across parities and stages of the lactation. We used data from ∼370,000 test-day records of Canadian Holstein cows to produce a landscape of the heritability of FTIR spectra, 1,060 wavelengths in the near- and mid-infrared spectrum (5,011-925 cm^-1), by parity and month of the lactation (mo 1 to 3 and mo 1 to 6, respectively). The 2 regions of the spectrum associated with absorption of electromagnetic energy by water molecules were estimated to have very high phenotypic variances, very low heritabilities, and very low proportion of variance explained by herd-year-season (HYS) subclasses. The near- or short-wavelength infrared (SWIR: 5,066-3,672 cm^-1) region was also characterized by low heritability estimates, whereas the estimated proportion of the variance explained by HYS was high. The mid-wavelength infrared region (MWIR: 3,000-2,500 cm^-1) and the transition between mid and long-wavelength infrared region (MWIR-LWIR: 1,500-925 cm^-1) harbor several waves characterized by moderately high (≥0.4) heritabilities. Most of the high-heritability regions contained wavelengths that are reported to be associated with important milk metabolites and components. Interestingly, these 2 same regions tended to show more variability in heritabilities between parity and lactation stage. Second parity showed heritability patterns that were distinctly different from those of the first and third parities, whereas the first 2 mo of the lactation had clearly distinct heritability patterns compared with mo 3 to 6.

Prediction accuracy of direct and indirect approaches, and their relationships with prediction ability of calibration models

Milk infrared spectra are routinely used for phenotyping traits of interest through links developed between the traits and spectra. Predicted individual traits are then used in genetic analyses for estimated breeding value (EBV) or for phenotypic predictions using a single-trait mixed model; this approach is referred to as indirect prediction (IP). An alternative approach [direct prediction (DP)] is a direct genetic analysis of (a reduced dimension of) the spectra using a multitrait model to predict multivariate EBV of the spectral components and, ultimately, also to predict the univariate EBV or phenotype for the traits of interest. We simulated 3 traits under different genetic (low: 0.10 to high: 0.90) and residual (zero to high: ±0.90) correlation scenarios between the 3 traits and assumed the first trait is a linear combination of the other 2 traits. The aim was to compare the IP and DP approaches for predictions of EBV and phenotypes under the different correlation scenarios. We also evaluated relationships between performances of the 2 approaches and the accuracy of calibration equations. Moreover, the effect of using different regression coefficients estimated from simulated phenotypes (β_p), true breeding values (β_g), and residuals (β_r) on performance of the 2 approaches were evaluated. The simulated data contained 2,100 parents (100 sires and 2,000 cows) and 8,000 offspring (4 offspring per cow). Of the 8,000 observations, 2,000 were randomly selected and used to develop links between the first and the other 2 traits using partial least square (PLS) regression analysis. The different PLS regression coefficients, such as β_p, β_g, and β_r, were used in subsequent predictions following the IP and DP approaches. We used BLUP analyses for the remaining 6,000 observations using the true (co)variance components that had been used for the simulation. Accuracy of prediction (of EBV and phenotype) was calculated as a correlation between predicted and true values from the simulations. The results showed that accuracies of EBV prediction were higher in the DP than in the IP approach. The reverse was true for accuracy of phenotypic prediction when using β_p but not when using β_g and β_r, where accuracy of phenotypic prediction in the DP was slightly higher than in the IP approach. Within the DP approach, accuracies of EBV when using β_g were higher than when using β_p only at the low genetic correlation scenario. However, we found no differences in EBV prediction accuracy between the β_p and β_g in the IP approach. Accuracy of the calibration models increased with an increase in genetic and residual correlations between the traits. Performance of both approaches increased with an increase in accuracy of the calibration models. In conclusion, the DP approach is a good strategy for EBV prediction but not for phenotypic prediction, where the classical PLS regression-based equations or the IP approach provided better results.

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.