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.

Selection indices for residual feed intake derived from milk spectra

Improving production efficiency and minimizing the environmental impact of dairy farming are 2 important goals of the dairy industry. Achieving these objectives requires improving the feed-to-milk conversion efficiency. One way to achieve this goal is through genetic selection. However, measuring feed efficiency in commercial herds is currently not feasible. As such, we conducted a study to evaluate the genetic accuracy of various selection indices derived from Fourier transform mid-infrared (FTIR)-spectra or milk composition. We use 7,793 weekly records on 537 genotyped cows (78,964 SNPs), with information on residual feed intake (RFI), and FTIR-spectra. We fitted various types of selection indexes using the complete FTIR-spectra of milk samples. The estimated heritability of RFI was 0.12 ± 0.02. The accuracy of indirect selection using the FTIR-spectra was maximized using a principal components selection index (0.16 ± 0.07), followed by a Lasso-type penalized selection index (0.14 ± 0.06). We determined that an index based on milk spectral data recorded on ~25 daughters produced a progeny average with an accuracy comparable to direct phenotypic selection for RFI. We conclude that indirect selection for RFI using FTIR-spectra data can be effective for sires with progeny; however, future studies with a larger sample size are needed to validate these results.

Infrared spectroscopy coupled with machine learning algorithms for predicting the detailed milk mineral profile in dairy cattle

Milk minerals are not only essential components for human health, but they can be informative for milk quality and cow's health. Herein, we investigated the feasibility of Fourier Transformed mid Infrared (FTIR) spectroscopy for the prediction of a detailed panel of 17 macro, trace, and environmental elements in bovine milk, using partial least squares regression (PLS) and machine learning approaches. The automatic machine learning significantly outperformed the PLS regression in terms of prediction performances of the mineral elements. For macrominerals, the R2 ranged from 0.59 to 0.78. Promising predictability was achieved for Cu and B (R2 = 0.66 and 0.74, respectively) and more moderate ones for Fe, Mn, Zn, and Al (R2 from 0.48 to 0.58). These results provide a reliable basis for a rapid and cost-effective quantification of these traits, serving as a resource for dairy farmers seeking to enhance the quality of milk production and optimize cheese properties.

Rapid detection and quantification of melamine, urea, sucrose, water, and milk powder adulteration in pasteurized milk using Fourier transform infrared (FTIR) spectroscopy coupled with modern statistical machine learning algorithms

There is an evident requirement for a rapid, efficient, and simple method to screen the authenticity of milk products in the market. Fourier transform infrared (FTIR) spectroscopy stands out as a promising solution. This work employed FTIR spectroscopy and modern statistical machine learning algorithms for the identification and quantification of pasteurized milk adulteration. Comparative results demonstrate modern statistical machine learning algorithms will improve the ability of FTIR spectroscopy to predict milk adulteration compared to partial least square (PLS). To discern the types of substances utilized in milk adulteration, a top-performing multiclassification model was established using multi-layer perceptron (MLP) algorithm, delivering an impressive prediction accuracy of 97.4 %. For quantification purposes, bayesian regularized neural networks (BRNN) provided the best results for the determination of both melamine, urea and milk powder adulteration, while extreme gradient boosting (XGB) and projection pursuit regression (PPR) gave better results in predicting sucrose and water adulteration levels, respectively. The regression models provided suitable predictive accuracy with the ratio of performance to deviation (RPD) values higher than 3. The proposed methodology proved to be a cost-effective and fast tool for screening the authenticity of pasteurized milk in the market.

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.

Lactation modeling and the effects of rotational crossbreeding on milk production traits and milk-spectra-predicted enteric methane emissions

Rotational crossbreeding has not been widely studied in relation to the enteric methane emissions of dairy cows, nor has the variation in emissions during lactation been modeled. Milk infrared spectra could be used to predict proxies of methane emissions in dairy cows. Therefore, the objective of this work was to study the effects of crossbreeding on the predicted infrared proxies of methane emissions and the variation in the latter during lactation. Milk samples were taken once from 1,059 cows reared in 2 herds, and infrared spectra of the milk were used to predict milk fat (mean ± SD; 3.79 ± 0.81%) and protein (3.68 ± 0.36%) concentrations, yield (21.4 ± 1.5 g/kg dry matter intake), methane intensity (14.2 ± 2.0 g/kg corrected milk), and daily methane production (358 ± 108 g/d). Of these cows, 620 were obtained from a 3-breed (Holstein, Montbéliarde, and Viking Red) rotational mating system, and the rest were purebred Holsteins. Milk production data and methane traits were analyzed using a nonlinear model that included the fixed effects of herd, genetic group, and parity, and the 4 parameters (a, b, c, and k) of a lactation curve modeled using the Wilmink function. Milk infrared spectra were found to be useful for direct prediction of qualitative proxies, such as methane yield and intensity, but not quantitative traits, such as daily methane production, which appears to be better estimated (450 ± 125 g/d) by multiplying a measured daily milk yield by infrared-predicted methane intensity. Lactation modeling of methane traits showed daily methane production to have a zenith curve, similar to that of milk yield but with a delayed peak (53 vs. 37 d in milk), whereas methane intensity is characterized by an upward curve that increases rapidly during the first third of lactation and then slowly till the end of lactation (10.5 g/kg at 1 d in milk to 15.2 g/kg at 300 d in milk). However, lactation modeling was not useful in explaining methane yield, which is almost constant during lactation. Lastly, the methane yield and intensity of cows from 3-breed rotational crossbreeding are not greater, and their methane production is lower than that of purebred Holsteins (452 vs. 477 g/d). Given the greater longevity of crossbred cows, and their lower replacement rate, rotational crossbreeding could be a way of mitigating the environmental impact of milk production.

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.

Cheese yield and nutrients recovery in the curd predicted by Fourier-transform spectra from individual sheep milk samples

The objectives of this study were to explore the use of Fourier-transform infrared (FTIR) spectroscopy on individual sheep milk samples for predicting cheese-making traits, and to test the effect of the farm variability on their prediction accuracy. For each of 121 ewes from 4 farms, a laboratory model cheese was produced, and 3 actual cheese yield traits (fresh cheese, cheese solids, and cheese water) and 4 milk nutrient recovery traits (fat, protein, total solids, and energy) in the curd were measured. Calibration equations were developed using a Bayesian approach with 2 different scenarios: (1) a random cross-validation (80% calibration; 20% validation set), and (2) a leave-one-out validation (3 farms used as calibration, and the remaining one as validation set) to assess the accuracy of prediction of samples from external farms, not included in calibration set. The best performance was obtained for predicting the yield and recovery of total solids, justifying for the practical application of the method at sheep population and dairy industry levels. Performances for the remaining traits were lower, but still useful for the monitoring of the milk processing in the case of fresh curd and recovery of energy. Insufficient accuracies were found for the recovery of protein and fat, highlighting the complex nature of the relationships among the milk nutrients and their recovery in the curd. The leave-one-out validation procedure, as expected, showed lower prediction accuracies, as a result of the characteristics of the farming systems, which were different between calibration and validation sets. In this regard, the inclusion of information related to the farm could help to improve the prediction accuracy of these traits. Overall, a large contribution to the prediction of the cheese-making traits came from the areas known as "water" and "fingerprint" regions. These findings suggest that, according to the traits studied, the inclusion of water regions for the development of the prediction equation models is fundamental to maintain a high prediction accuracy. However, further studies are necessary to better understand the role of specific absorbance peaks and their contribution to the prediction of cheese-making traits, to offer reliable tools applicable along the dairy ovine chain.

Recent Advances in the Determination of Milk Adulterants and Contaminants by Mid-Infrared Spectroscopy

The presence of chemical contaminants, toxins, or veterinary drugs in milk, as well as the adulteration of milk from different species, has driven the development of new tools to ensure safety and quality. Several analytical procedures have been proposed for the rapid screening of hazardous substances or the selective confirmation of the authenticity of milk. Mid-infrared spectroscopy and Fourier-transform infrared have been two of the most relevant technologies conventionally employed in the dairy industry. These fingerprint methodologies can be very powerful in determining the trait of raw material without knowing the identity of each constituent, and several aspects suggest their potential as a screening method to detect adulteration. This paper reviews the latest advances in applying mid-infrared spectroscopy for the detection and quantification of adulterants, milk dilution, the presence of pathogenic bacteria, veterinary drugs, and hazardous substances in milk.

Prediction of detailed blood metabolic profile using milk infrared spectra and machine learning methods in dairy cattle

The adoption of preventive management decisions is crucial to dealing with metabolic impairments in dairy cattle. Various serum metabolites are known to be useful indicators of the health status of cows. In this study, we used milk Fourier-transform mid-infrared (FTIR) spectra and various machine learning (ML) algorithms to develop prediction equations for a panel of 29 blood metabolites, including those related to energy metabolism, liver function/hepatic damage, oxidative stress, inflammation/innate immunity, and minerals. For most traits, the data set comprised observations from 1,204 Holstein-Friesian dairy cows belonging to 5 herds. An exception was represented by β-hydroxybutyrate prediction, which contained observations from 2,701 multibreed cows pertaining to 33 herds. The best predictive model was developed using an automatic ML algorithm that tested various methods, including elastic net, distributed random forest, gradient boosting machine, artificial neural network, and stacking ensemble. These ML predictions were compared with partial least squares regression, the most commonly used method for FTIR prediction of blood traits. Performance of each model was evaluated using 2 cross-validation (CV) scenarios: 5-fold random (CV_r) and herd-out (CV_h). We also tested the best model's ability to classify values precisely in the 2 extreme tails, namely, the 25th (Q25) and 75th (Q75) percentiles (true-positive prediction scenario). Compared with partial least squares regression, ML algorithms achieved more accurate performance. Specifically, elastic net increased the R² value from 5% to 75% for CV_r and 2% to 139% for CV_h, whereas the stacking ensemble increased the R² value from 4% to 70% for CV_r and 4% to 150% for CV_h. Considering the best model, with the CV_r scenario, good prediction accuracies were obtained for glucose (R² = 0.81), urea (R² = 0.73), albumin (R² = 0.75), total reactive oxygen metabolites (R² = 0.79), total thiol groups (R² = 0.76), ceruloplasmin (R² = 0.74), total proteins (R² = 0.81), globulins (R² = 0.87), and Na (R² = 0.72). Good prediction accuracy in classifying extreme values was achieved for glucose (Q25 = 70.8%, Q75 = 69.9%), albumin (Q25 = 72.3%), total reactive oxygen metabolites (Q25 = 75.1%, Q75 = 74%), thiol groups (Q75 = 70.4%), total proteins (Q25 = 72.4%, Q75 = 77.2.%), globulins (Q25 = 74.8%, Q75 = 81.5%), and haptoglobin (Q75 = 74.4%). In conclusion, our study shows that FTIR spectra can be used to predict blood metabolites with relatively good accuracy, depending on trait, and are a promising tool for large-scale monitoring.

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.

Phenotypic and genetic characterization of the occurrence of noncoagulating milk in dairy sheep

Milk coagulation ability is of central importance for the sheep dairy industry because almost all sheep milk is destined for cheese processing. The occurrence of milk with impaired coagulation properties is an obstacle to cheese processing and, in turn, to the profitability of the dairy companies. In this work, we investigated the causes of noncoagulation of sheep milk; specifically, we studied the effect of milk physicochemical properties on milk coagulation status [coagulating and noncoagulating (NC) milk samples, which do or do not coagulate within 30 min, respectively], and whether mid-infrared spectroscopy (MIR) could be used to assess variability in coagulation status. We also investigated the genetic background of milk coagulation ability. Individual milk samples were collected from 996 Sarda ewes farmed in 47 flocks located in Sardinia (Italy). Considered traits were daily milk yield, milk composition traits, and milk coagulation properties (rennet coagulation time, curd firming time, and curd firmness), and MIR spectra were acquired. About 9% of samples did not coagulate within 30 min. A logistic regression approach was used to test the effect of milk-related traits on milk coagulation status. A principal component (PC) analysis was carried out on the milk MIR spectra, and PC scores were then used as covariates in a logistic regression model to assess their relationship with milk coagulation status. Results of the present work demonstrated that the probability of having NC samples increases as milk contents of proteins and chlorides and somatic cell score increase. The analysis of PC extracted from milk spectra that influenced coagulation status highlighted key regions associated with lactose and protein concentrations, and others not associated with routinely collected milk composition traits. These results suggest that the occurrence of NC is mostly related to damage of the epithelium secretory mammary cells, which occurs with the advancement of a lactation or due to unhealthy mammary gland status. Genetic analysis of milk coagulation status and of the extracted PC confirmed the genetic background of the milk coagulability of sheep milk.

Multivariate and Genome-Wide Analysis of Mid-Infrared Spectra of Non-Coagulating Milk of Sarda Sheep Breed

Milk coagulation ability is crucial for the dairy sheep industry since the whole amount of milk is processed into cheese. Non-coagulating milk (NCM) is defined as milk not forming a curd within the testing time. In sheep milk, it has been reported in literature that up to 20% of milk is NCM. Although the clotting properties of individual milk have been widely studied, little attention has been given to NCM and genomic dissection of this trait. Mid-infrared (MIR) spectra can be exploited both to predict cheese-making aptitude and to discriminate between coagulating milk and NCM. The main goals of this work were (i) to assess the predictivity of MIR spectra for NCM classification and (ii) to conduct a genome-wide association study on coagulation ability. Milk samples from 949 Sarda ewes genotyped and phenotyped for milk coagulation properties (MCPs) served as the training dataset. The validation dataset included 662 ewes. Three classical MCPs were measured: rennet coagulation time (RCT), curd firmness (a30), and curd firming time (k20). Moreover, MIR spectra were acquired and stored in the region between 925.92 and 5,011.54 cm−1. The probability of a sample to be NCM was modeled by step-wise logistic regression on milk spectral information (LR-W), logistic regression on principal component (LR-PC), and canonical discriminant analysis of spectral wave number (DA-W). About 9% of the samples did not coagulate at 30 min. The use of LR-W gave a poorer classification of NCM. The use of LR-PC improved the percentage of correct assignment (45 ± 9%). The DA-W method allows us to reach 75.1 ± 10.3 and 76.5 ± 18.4% of correct assignments of the inner and external validation datasets, respectively. As far as GWA of NCM, 458 SNP associations and 45 candidate genes were detected. The genes retrieved from public databases were mostly linked to mammary gland metabolism, udder health status, and a milk compound also known to affect the ability of milk to coagulate. In particular, the potential involvement of CAPNs deserves further investigation.

Challenging Sustainable and Innovative Technologies in Cheese Production: A Review

It is well known that cheese yield and quality are affected by animal genetics, milk quality (chemical, physical, and microbiological), production technology, and the type of rennet and dairy cultures used in production. Major differences in the same type of cheese (i.e., hard cheese) are caused by the rennet and dairy cultures, which affect the ripening process. This review aims to explore current technological advancements in animal genetics, methods for the isolation and production of rennet and dairy cultures, along with possible applications of microencapsulation in rennet and dairy culture production, as well as the challenge posed to current dairy technologies by the preservation of biodiversity. Based on the reviewed scientific literature, it can be concluded that innovative approaches and the described techniques can significantly improve cheese production.

Invited review: A comprehensive review of visible and near-infrared spectroscopy for predicting the chemical composition of cheese

Substantial research has been carried out on rapid, nondestructive, and inexpensive techniques for predicting cheese composition using spectroscopy in the visible and near-infrared radiation range. Moreover, in recent years, new portable and handheld spectrometers have been used to predict chemical composition from spectra captured directly on the cheese surface in dairies, storage facilities, and food plants, removing the need to collect, transport, and process cheese samples. For this review, we selected 71 papers (mainly dealing with prediction of the chemical composition of cheese) and summarized their results, focusing our attention on the major sources of variation in prediction accuracy related to cheese variability, spectrometer and spectra characteristics, and chemometrics techniques. The average coefficient of determination obtained from the validation samples ranged from 86 to 90% for predicting the moisture, fat, and protein contents of cheese, but was lower for predicting NaCl content and cheese pH (79 and 56%, respectively). There was wide variability with respect to all traits in the results of the various studies (standard deviation: 9-30%). This review draws attention to the need for more robust equations for predicting cheese composition in different situations; the calibration data set should consist of representative cheese samples to avoid bias due to an overly specific field of application and ensure the results are not biased for a particular category of cheese. Different spectrometers have different accuracies, which do not seem to depend on the spectrum extension. Furthermore, specific areas of the spectrum-the visible, infrared-A, or infrared-B range-may yield similar results to broad-range spectra; this is because several signals related to cheese composition are distributed along the spectrum. Small, portable instruments have been shown to be viable alternatives to large bench-top instruments. Last, chemometrics (spectra pre-treatment and prediction models) play an important role, especially with regard to difficult-to-predict traits. A proper, fully independent, validation strategy is essential to avoid overoptimistic results.

Experimental Treatment of Hazardous Ash Waste by Microbial Consortium Aspergillus niger and Chlorella sp.: Decrease of the Ni Content and Identification of Adsorption Sites by Fourier-Transform Infrared Spectroscopy

Despite the negative impact on the environment, incineration is one of the most commonly used methods for dealing with waste. Besides emissions, the production of ash, which usually shows several negative properties, such as a higher content of hazardous elements or strongly alkaline pH, is problematic from an environmental viewpoint as well. The subject of our paper was the assessment of biosorption of Ni from ash material by a microbial consortium of Chlorella sp. and Aspergillus niger. The solid substrate represented a fraction of particles of size <0.63 mm with a Ni content of 417 mg kg-1. We used a biomass consisting of two different organisms as the sorbent: a non-living algae culture of Chlorella sp. (an autotrophic organism) and the microscopic filamentous fungus A. niger (a heterotrophic organism) in the form of pellets. The experiments were conducted under static conditions as well as with the use of shaker (170 rpm) with different modifications: solid substrate, Chlorella sp. and pellets of A. niger; solid substrate and pellets of A. niger. The humidity-temperature conditions were also changed. Sorption took place under dry and also wet conditions (with distilled water in a volume of 30-50 ml), partially under laboratory conditions at a temperature of 25°C as well as in the exterior. The determination of the Ni content was done using inductively coupled plasma optical emission spectrometry (ICP-OES). The removal of Ni ranged from 13.61% efficiency (Chlorella sp., A. niger with the addition of 30 ml of distilled water, outdoors under static conditions after 48 h of the experiment) to 46.28% (Chlorella sp., A. niger with the addition of 30 ml of distilled water, on a shaker under laboratory conditions after 48 h of the experiment). For the purpose of analyzing the representation of functional groups in the microbial biomass and studying their interaction with the ash material, we used Fourier-transform infrared (FTIR) spectroscopy. We observed that the amount of Ni adsorbed positively correlates with absorbance in the spectral bands where we detect the vibrations of several organic functional groups. These groups include hydroxyl, aliphatic, carbonyl, carboxyl and amide structural units. The observed correlations indicate that, aside from polar and negatively charged groups, aliphatic or aromatic structures may also be involved in sorption processes due to electrostatic attraction. The correlation between absorbance and the Ni content reached a maximum in amide II band (r = 0.9; P < 0.001), where vibrations of the C=O, C-N, and N-H groups are detected. The presented results suggest that the simultaneous use of both microorganisms in biosorption represents an effective method for reducing Ni content in a solid substrate, which may be useful as a partial process for waste disposal.

Application of machine-learning methods to milk mid-infrared spectra for discrimination of cow milk from pasture or total mixed ration diets

The prevalence of "grass-fed" labeled food products on the market has increased in recent years, often commanding a premium price. To date, the majority of methods used for the authentication of grass-fed source products are driven by auditing and inspection of farm records. As such, the ability to verify grass-fed source claims to ensure consumer confidence will be important in the future. Mid-infrared (MIR) spectroscopy is widely used in the dairy industry as a rapid method for the routine monitoring of individual herd milk composition and quality. Further harnessing the data from individual spectra offers a promising and readily implementable strategy to authenticate the milk source at both farm and processor levels. Herein, a comprehensive comparison of the robustness, specificity, and accuracy of 11 machine-learning statistical analysis methods were tested for the discrimination of grass-fed versus non-grass-fed milks based on the MIR spectra of 4,320 milk samples collected from cows on pasture or indoor total mixed ration-based feeding systems over a 3-yr period. Linear discriminant analysis and partial least squares discriminant analysis (PLS-DA) were demonstrated to offer the greatest level of accuracy for the prediction of cow diet from MIR spectra. Parsimonious strategies for the selection of the most discriminating wavelengths within the spectra are also highlighted.

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.

Evaluating the performance of machine learning methods and variable selection methods for predicting difficult-to-measure traits in Holstein dairy cattle using milk infrared spectral data

Fourier-transform infrared (FTIR) spectroscopy is a powerful high-throughput phenotyping tool for predicting traits that are expensive and difficult to measure in dairy cattle. Calibration equations are often developed using standard methods, such as partial least squares (PLS) regression. Methods that employ penalization, rank-reduction, and variable selection, as well as being able to model the nonlinear relations between phenotype and FTIR, might offer improvements in predictive ability and model robustness. This study aimed to compare the predictive ability of 2 machine learning methods, namely random forest (RF) and gradient boosting machine (GBM), and penalized regression against PLS regression for predicting 3 phenotypes differing in terms of biological meaning and relationships with milk composition (i.e., phenotypes measurable directly and not directly in milk, reflecting different biological processes which can be captured using milk spectra) in Holstein-Friesian cattle under 2 cross-validation scenarios. The data set comprised phenotypic information from 471 Holstein-Friesian cows, and 3 target phenotypes were evaluated: (1) body condition score (BCS), (2) blood β-hydroxybutyrate (BHB, mmol/L), and (3) κ-casein expressed as a percentage of nitrogen (κ-CN, % N). The data set was split considering 2 cross-validation scenarios: samples-out random in which the population was randomly split into 10-folds (8-folds for training and 1-fold for validation and testing); and herd/date-out in which the population was randomly assigned to training (70% herd), validation (10%), and testing (20% herd) based on the herd and date in which the samples were collected. The random grid search was performed using the training subset for the hyperparameter optimization and the validation set was used for the generalization of prediction error. The trained model was then used to assess the final prediction in the testing subset. The grid search for penalized regression evidenced that the elastic net (EN) was the best regularization with increase in predictive ability of 5%. The performance of PLS (standard model) was compared against 2 machine learning techniques and penalized regression using 2 cross-validation scenarios. Machine learning methods showed a greater predictive ability for BCS (0.63 for GBM and 0.61 for RF), BHB (0.80 for GBM and 0.79 for RF), and κ-CN (0.81 for GBM and 0.80 for RF) in samples-out cross-validation. Considering a herd/date-out cross-validation these values were 0.58 (GBM and RF) for BCS, 0.73 (GBM and RF) for BHB, and 0.77 (GBM and RF) for κ-CN. The GBM model tended to outperform other methods in predictive ability around 4%, 1%, and 7% for EN, RF, and PLS, respectively. The prediction accuracies of the GBM and RF models were similar, and differed statistically from the PLS model in samples-out random cross-validation. Although, machine learning techniques outperformed PLS in herd/date-out cross-validation, no significant differences were observed in terms of predictive ability due to the large standard deviation observed for predictions. Overall, GBM achieved the highest accuracy of FTIR-based prediction of the different phenotypic traits across the cross-validation scenarios. These results indicate that GBM is a promising method for obtaining more accurate FTIR-based predictions for different phenotypes in dairy cattle.

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.

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.

Association between days open and milk spectral data in dairy cows

Data on 19,489 Brown Swiss cows reared in northeastern Italy were used to associate absorbances of individual wavenumbers within the mid-infrared range with days open (DO). Different postcalving days in milk (DIM) intervals were studied to determine the most informative milk sampling periods for predicting DO. Milk samples were analyzed using a MilkoScan (Foss Electric, Hillerød, Denmark) Fourier-transform infrared (FTIR) spectrometer for 1,060 wavenumbers (wn) ranging from 5,011 to 925 cm^-1. To determine DO, we considered an insemination to lead to conception when there was no return of heat (i.e., no successive insemination) and the cow had a subsequent calving date whereby gestation length was required to be within ±30 d of 290 d. Only milk records within the first 90 DIM were considered. Associations were inferred by (1) fitting linear regression models between the DO and each individual wavenumber or milk component, and (2) fitting a Bayesian regression model that included the complete FTIR spectral data. The effects of including systematic effects (parity number, year-season, herd) in the model on these associations were also studied. These analyses were performed for the complete data (5-90 DIM) and for data stratified by DIM period (5 to 30, 31 to 60, and 61 to 90 DIM). Overall, regions of wavenumbers of the milk FTIR spectra that were associated with DO included wn 2,973 to 2,830 cm^-1 [related to fat-B (C-H stretch)], wn 2,217 to 1,769 cm^-1 [related to fat-A (C = O stretch)], wn 1,546 cm^-1 (related to protein), wn 1,465 cm^-1 (related to urea and fat), wn 1,399 to 1,245 cm^-1 (related to acetone), and wn 1,110 cm^-1 (related to lactose). Estimated effects depended on the DIM period, with milk samples drawn during DIM intervals 31 to 60 d and 61 to 90 d being most strongly associated with DO. These DIM intervals are also typically most associated with negative energy balance and peak lactation.

Predicting the Content of 20 Minerals in Beef by Different Portable Near-Infrared (NIR) Spectrometers

The aim of this study was to test the predictability of a detailed mineral profile of beef using different portable near-infrared spectrometers (NIRS). These devices are rapid, chemical waste-free, cheap, nondestructive tools that can be used directly on the meat surface in the work environment without the need to take samples. We compared a transportable Visible-NIRS (weight 5.6 kg; wavelength 350-1830 nm), a portable NIRS (2.0 kg; 950-1650 nm), and a hand-held Micro-NIRS (0.06 kg; 905-1649 nm) to predict the contents of 20 minerals (measured by ICP-OES) in 178 beef samples (Longissimus thoracis muscle) using different mathematical pretreatments of the spectra and partial least square regressions. The externally validated results show that Fe, P, Mg, S, Na, and Pb have some potential for prediction with all instruments (R2VAL: 0.40-0.83). Overall, the prediction performances of the three instruments were similar, although the smallest (Micro-NIRS) exhibited certain advantages.

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

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

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.

Ability of near-infrared spectroscopy and chemometrics to predict the age of mosquitoes reared under different conditions

Background

Practical, field-ready age-grading tools for mosquito vectors of disease are urgently needed because of the impact that daily survival has on vectorial capacity. Previous studies have shown that near-infrared spectroscopy (NIRS), in combination with chemometrics and predictive modeling, can forecast the age of laboratory-reared mosquitoes with moderate to high accuracy. It remains unclear whether the technique has utility for identifying shifts in the age structure of wild-caught mosquitoes. Here we investigate whether models derived from the laboratory strain of mosquitoes can be used to predict the age of mosquitoes grown from pupae collected in the field.

Methods

NIRS data from adult female Aedes albopictus mosquitoes reared in the laboratory (2, 5, 8, 12 and 15 days-old) were analysed against spectra from mosquitoes emerging from wild-caught pupae (1, 7 and 14 days-old). Different partial least squares (PLS) regression methods trained on spectra from laboratory mosquitoes were evaluated on their ability to predict the age of mosquitoes from more natural environments.

Results

Models trained on spectra from laboratory-reared material were able to predict the age of other laboratory-reared mosquitoes with moderate accuracy and successfully differentiated all day 2 and 15 mosquitoes. Models derived with laboratory mosquitoes could not differentiate between field-derived age groups, with age predictions relatively indistinguishable for day 1–14. Pre-processing of spectral data and improving the PLS regression framework to avoid overfitting can increase accuracy, but predictions of mosquitoes reared in different environments remained poor. Principal components analysis confirms substantial spectral variations between laboratory and field-derived mosquitoes despite both originating from the same island population.

Conclusions

Models trained on laboratory mosquitoes were able to predict ages of laboratory mosquitoes with good sensitivity and specificity though they were unable to predict age of field-derived mosquitoes. This study suggests that laboratory-reared mosquitoes do not capture enough environmental variation to accurately predict the age of the same species reared under different conditions. Further research is needed to explore alternative pre-processing methods and machine learning techniques, and to understand factors that affect absorbance in mosquitoes before field application using NIRS.

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.

Chemometric authentication of farming systems of origin of food (milk and ripened cheese) using infrared spectra, fatty acid profiles, flavor fingerprints, and sensory descriptions

Milk samples from 1264 cows in 85 farms were authenticated for different farming-systems using a 10-fold cross-validated linear-discriminant-analysis using Fourier-transform infrared spectra (FTIRS) and gas-chromatographic fatty-acid (FA) profiles. FTIRS gave correct classification greater than FAs (97.4% vs. 81.1%) during calibration, but slightly worse in validation (73.5% vs 77.3%) and their combination improved the results. All milk samples were processed into ripened model-cheeses, and analyzed by near-infrared-spectrometry (NIRS), by proton-transfer-reaction time-of-flight mass-spectrometry for their volatile organic compound (VOCs) fingerprint and by panel sensory profiling (SENS). Farming-system authentication on cheese samples was less efficient than on milk, but still possible. The instrumental methods yielded similar validation results, better than SENS, and their combination improved the correct classification rate. The efficiency of the different technics was affected by specific farming systems. In conclusion, dairy products could be discriminated for farming-systems with acceptable accuracy, but the methods tested differ in sampling procedure, rapidity and costs.

Genetic variance of metabolomic features and their relationship with body weight and body weight gain in Holstein cattle1

In recent years, metabolomics has been used to clarify the biology underlying biological samples. In the field of animal breeding, investigating the magnitude of genetic control on the metabolomic profiles of animals and their relationships with quantitative traits adds valuable information to animal improvement schemes. In this study, we analyzed metabolomic features (MFs) extracted from the metabolomic profiles of 843 male Holstein calves. The metabolomic profiles were obtained using nuclear magnetic resonance (NMR) spectroscopy. We investigated 2 alternative methods to control for peak shifts in the NMR spectra, binning and aligning, to determine which approach was the most efficient for assessing genetic variance. Series of univariate analyses were implemented to elucidate the heritability of each MF. Furthermore, records on BW and ADG from 154 to 294 d of age (ADG154-294), 294 to 336 d of age (ADG294-336), and 154 to 336 d of age (ADG154-336) were used in a series of bivariate analyses to establish the genetic and phenotypic correlations with MFs. Bivariate analyses were only performed for MFs that had a heritability significantly different from zero. The heritabilities obtained in the univariate analyses for the MFs in the binned data set were low (<0.2). In contrast, in the aligned data set, we obtained moderate heritability (0.2 to 0.5) for 3.5% of MFs and high heritability (more than 0.5) for 1% of MFs. The bivariate analyses showed that ~12%, ~3%, ~9%, and ~9% of MFs had significant additive genetic correlations with BW, ADG154-294, ADG294-336, and ADG154-336, respectively. In all of the bivariate analyses, the percentage of significant additive genetic correlations was higher than the percentage of significant phenotypic correlations of the corresponding trait. Our results provided insights into the influence of the underlying genetic mechanisms on MFs. Further investigations in this field are needed for better understanding of the genetic relationship among the MFs and quantitative traits.

Accuracy and biases in predicting the chemical and physical traits of many types of cheeses using different visible and near-infrared spectroscopic techniques and spectrum intervals

Near-infrared spectroscopy (NIRS) has been widely used to determine various composition traits of many dairy products in the industry. In the last few years, near-infrared (NIR) instruments have become more and more accessible, and now, portable devices can be easily used in the field, allowing the direct measurement of important quality traits. However, the comparison of the predictive performances of different NIR instruments is not simple, and the literature is lacking. These instruments may use different wavelength intervals and calibration procedures, making it difficult to establish whether differences are due to the spectral interval, the chemometric approach, or the instrument's technology. Hence, the aims of this study were (1) to evaluate the prediction accuracy of chemical contents (5 traits), pH, texture (2 traits), and color (5 traits) of 37 categories of cheese; (2) to compare 3 instruments [2 benchtop, working in reflectance (R) and transmittance (T) mode (NIRS-R and NIRS-T, respectively) and 1 portable device (VisNIRS-R)], using their entire spectral ranges (1100-2498, 850-1048, and 350-1830 nm, respectively, for NIRS-R, NIRS-T and VisNIRS-R); (3) to examine different wavelength intervals of the spectrum within instrument, comparing also the common intervals among the 3 instruments; and (4) to determine the presence of bias in predicted traits for specific cheese categories. A Bayesian approach was used to develop 8 calibration models for each of 13 traits. This study confirmed that NIR spectroscopy can be used to predict the chemical composition of a large number of different cheeses, whereas pH and texture traits were poorly predicted. Color showed variable predictability, according to the trait considered, the instrument used, and, within instrument, according to the wavelength intervals. The predictive performance of the VisNIRS-R portable device was generally better than the 2 laboratory NIRS instruments, whether with the entire spectrum or selected intervals. The VisNIRS-R was found suitable for analyzing chemical composition in real time, without the need for sample uptake and processing. Our results also indicated that instrument technology is much more important than the NIR spectral range for accurate prediction equations, but the visible range is useful when predicting color traits, other than lightness. Specifically for certain categories (i.e., caprine, moldy, and fresh cheeses), dedicated calibrations seem to be needed to obtain unbiased and more accurate results.

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.

RESEARCH OF ORGANOLEPTIC PARAMETERS OF DUTCH CHEESE, PRODUCED FROM MILK OF COWS OF DIFFERENT BREEDS

The aim of this research is to study an influence of types of different breeds of milk cows on organoleptic properties of Dutch cheese at different rations of their food. It allows to receive cheese with prognosticated quality characteristics and to manage marketing strategies at cheese realization. The work presents the results of studying organoleptic properties of Dutch cheese, produced of milk of cows of different breeds, traditionally bred in the Northern-Eastern region of Ukraine. For studying an influence of cows feeding on exclusive properties of hard cheese, traditional food rations of cows were added with forages of a silage-hay type with Lucerne as a main source of proteins (44 % of the daily norm). Rations were practically identical with traditional ones by the content of energy and main food value factors and in general corresponded to norms of cows feeding. It has been demonstrated, that the high mark of cheese, up to 45 points for the taste and smell and total point mark 99 of 100 possible was obtained as a result of changing food rations at the expanse of introducing Lucerne silage in them. The data on optimization of food rations of certain cow breeds allow managing quality characteristics of milk and products of and are expedient for making cheese with unique regional characteristics. At conducting the comparative organoleptic assessment, there has been revealed a distinct dependence between a point mark of quality parameters of cheese on a breed of milk cows, and also on food rations of them. At the silage-hay food ration of animals, the received Dutch cheese had higher quality characteristics after 60 days of storage than cheese samples, obtained of cow milk at the traditional food ration.

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.