Short communication: Selecting the most informative mid-infrared spectra wavenumbers to improve the accuracy of prediction models for detailed milk protein content

SVR Chemometrics to Quantify β-Lactoglobulin and α-Lactalbumin in Milk Using MIR

Protein content variation in milk can impact the quality and consistency of dairy products, necessitating access to in-line real time monitoring. Here, we present a chemometric approach for the qualitative and quantitative monitoring of β-lactoglobulin and α-lactalbumin, using mid-infrared spectroscopy (MIR). In this study, we employed Hotelling T2 and Q-residual for outlier detection, automated preprocessing using nippy, conducted wavenumber selection with genetic algorithms, and evaluated four chemometric models, including partial least squares, support vector regression (SVR), ridge, and logistic regression to accurately predict the concentrations of β-lactoglobulin and α-lactalbumin in milk. For the quantitative analysis of these two whey proteins, SVR performed the best to interpret protein concentration from 197 MIR spectra originating from 42 Cornell University samples of preserved pasteurized modified milk. The R2 values obtained for β-lactoglobulin and α-lactalbumin using leave one out cross-validation (LOOCV) are 92.8% and 92.7%, respectively, which is the highest correlation reported to date. Our approach introduced a combination of preprocessing automation, genetic algorithm-based wavenumber selection, and used Optuna to optimize the framework for tuning hyperparameters of the chemometric models, resulting in the best chemometric analysis of MIR data to quantitate β-lactoglobulin and α-lactalbumin to date.

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

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

Dairy Cows' Health during Alpine Summer Grazing as Assessed by Milk Traits, Including Differential Somatic Cell Count: A Case Study from Italy

Simple Summary

Dairy herds in alpine areas usually adopt summer grazing, mainly to reduce feeding costs. This practice is related to the maintenance of local traditions and to the manufacturing of niche dairy products. However, it is important to assess the impact of this practice on cattle health. This case study investigated how milk-related health traits vary across extensive grazing during the summer period, using data collected in a dairy herd whose cows were repeatedly controlled for individual milk samples. Although the transition from barn farming to pasture led to a reduction in milk production, proper grazing management can make dairy cows more resilient in terms of udder health and metabolic efficiency. Findings of the present research report suggested that pasture can be adopted to maintain dairy herd sustainability without impairing animal health.

Abstract

Extensive summer grazing is a dairy herd management practice frequently adopted in mountainous areas. Nowadays, this activity is threatened by its high labour demand, but it is fundamental for environmental, touristic and economic implications, as well as for the preservation of social and cultural traditions. Scarce information on the effects of such low-input farming systems on cattle health is available. Therefore, the present case study aimed at investigating how grazing may affect the health status of dairy cows by using milk traits routinely available from the national milk recording scheme. The research involved a dairy herd of 52 Simmental and 19 Holstein × Simmental crossbred cows. The herd had access to the pasture according to a rotational grazing scheme from late spring up to the end of summer. A total of 616 test day records collected immediately before and during the grazing season were used. Individual milk yield was registered during the milking procedure. Milk samples were analysed for composition (fat, protein, casein and lactose contents) and health-related milk indicators (electrical conductivity, urea and β-hydroxybutyrate) using mid-infrared spectroscopy. Somatic cell count (SCC) and differential SCC were also determined. Data were analysed with a linear mixed model, which included the fixed effects of the period of sampling, cow breed, stage of lactation and parity, and the random effects of cow nested within breed and the residual. The transition from barn farming to pasture had a negative effect on milk yield, together with a small deterioration of fat and protein percentages. Health-related milk indicators showed a minor deterioration of the fat to protein ratio, differential SCC and electrical conductivity, particularly towards the end of the grazing season, whereas the somatic cell score and β-hydroxybutyrate were relatively constant. Overall, the study showed that, when properly managed, pasture grazing does not have detrimental effects on dairy cows in terms of udder health and efficiency. Therefore, the proper management of cows on pasture can be a valuable solution to preserve the economic, social and environmental sustainability of small dairy farms in the alpine regions, without impairing cows’ health.

Genetic and Non-Genetic Variation of Milk Total Antioxidant Activity Predicted from Mid-Infrared Spectra in Holstein Cows

Simple Summary

The total antioxidant activity (TAA) of food is important for human health and results from the contribution of different nutraceutical compounds. Direct determination of TAA in food is time-consuming and expensive. Infrared technologies allow the prediction of difficult-to-measure traits with certain accuracy in several organic matrices, including TAA of bovine milk. In order to understand the background of TAA and identify potential strategies to improve this feature in bovine milk, we explored its non-genetic sources of variation and estimated heritability and correlations with traits of economic interest in a large database of Holstein cows.

Abstract

Food antioxidants enhance products shelf life and stability during technological treatments through the maintenance of their physical and chemical properties. Moreover, they are endowed with several positive effects on human health, including cell membranes preservation, enzyme functionality, and DNA integrity. Milk has been described in relation to a wide array of fat soluble and water-soluble antioxidant compounds, in particular vitamin A, C, and E, lactoferrin and peptides derived from casein and whey proteins. The total antioxidant activity (TAA) of milk is a novel and scarcely explored trait, defined as the sum of antioxidant contributions of the aforementioned compounds. On this background, the aims of the present study were to investigate the variability of milk TAA on a large scale exploiting predictions obtained through mid-infrared (MIR) spectroscopy and to estimate genetic parameters of this trait in Holstein cows. Individual milk samples were collected between January 2011 and December 2018 during the routine milk recording procedure. Samples were analysed for gross composition through MIR spectroscopy and MIR spectra were stored. Milk TAA was then predicted (pTAA) from the stored milk MIR spectra (111,653 test-day records of 9519 cows in 344 herds) using the previously developed prediction model; considering the prediction accuracy, pTAA might be considered a proxy of the TAA determined through the reference method. Overall, pTAA averaged 7.16 mmoL/L of Trolox equivalents, showed a nadir around 40 days after calving and increased thereafter, following a linear trend up to the end of lactation. The lowest pTAA was observed in milk sampled from June to September. Milk pTAA was heritable (0.401 ± 0.015) and genetically associated to fat yield (0.366 ± 0.049), crude protein (CP) yield (0.238 ± 0.052), fat percentage (0.616 ± 0.022) and CP percentage (0.754 ± 0.015). The official selection index of Italian Holstein put the 49% of the emphasis on fat and protein yield and percentage; therefore, it derives that an indirect favourable selection for milk pTAA should be already in progress in Italian Holstein population.

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

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

Optical Methods Based on Ultraviolet, Visible, and Near-Infrared Spectra to Estimate Fat and Protein in Raw Milk: A Review

The present manuscript focuses on reviewing the optical techniques proposed to monitor milk quality in dairy farms to increase productivity and reduce costs. As is well known, the quality is linked to the fat and protein concentration; in addition, this issue is crucial to maintaining a healthy herd and preventing illnesses such as mastitis and ketosis. Usually, the quality of the milk is carried out with invasive methods employing chemical reagents that increase the time analysis. As a solution, several spectroscopy optical methods have been proposed, here, the benefits such as non-invasive measurement, online implementation, rapid estimation, and cost-effective execution. The most attractive optical methods to estimate fat and protein in cow's milk are compared and discussed considering their performance. The analysis is divided considering the wavelength operation (ultraviolet, visible, and infrared). Moreover, the weaknesses and strengths of the methods are fully analyzed. Finally, we provide the trends and a recent technique based on spectroscopy in the visible wavelength.

Effects of somatic cell score on milk yield and mid-infrared predicted composition and technological traits of Brown Swiss, Holstein Friesian, and Simmental cattle breeds

High milk somatic cell count (SCC) influences milk production and quality; however, very little is known about the effect of low SCC on milk quality, especially in terms of mineral content and coagulation properties. Thus, the present study aimed to investigate the effects of somatic cell score (SCS), calculated as log₂(SCC/100) + 3, on milk yield, composition (fat, crude protein, casein, lactose, milk urea nitrogen, protein fractions, and mineral contents), and coagulation properties of Brown Swiss, Holstein Friesian, and Simmental cows from multibreed herds. Milk composition and coagulation traits were predicted using mid-infrared spectroscopy. The data set comprised 95,591 observations of 6,940 cows in 313 multibreed herds, collected from January 2011 to December 2017. Observations were divided into 8 classes based on SCS. Statistical analysis was performed using a linear mixed model, which included breed, parity, stage of lactation, SCS class, and their interactions as fixed effects, and herd test day, cow, and residual as random effects. The probability that cows experienced SCS > 4.00 at 30 ± 5, 60 ± 5, or 90 ± 5 d after the observation test day was calculated for each SCS class, and odds ratios to the reference class (-1.00 < SCS ≤ 0.00) were reported. Results showed that the relationship between SCS and milk traits followed a third-order polynomial regression. The average loss of milk, fat, and crude protein yields were 0.43, 0.01, and 0.01 kg/d, respectively, for each SCS unit higher than 1.00. Very low SCS (<-1.00) had detrimental effects on milk yield and quality traits similar to or even stronger than high SCS (>4.00). Moreover, cows with SCS lower than -1.00 on a test day were about 7 times more likely to present high SCS within the following 90 ± 5 d than cows with SCS between -1.00 and 0.00. Breeds responded similarly to the increase of SCS, but the overall loss of fat and crude protein yields, and several minerals among Holstein Friesian were lower with increasing SCS. The best milk yield and quality were observed between SCS 0.00 and 1.00, but milk quality of Holstein Friesians started to decrease at lower SCS compared with milk quality of Brown Swiss and Simmental cows. Results suggest a breed-dependent optimum of SCS, and highlighted that very low SCS can be an indicator of udder health problems and, thus, may be used for early detection of mastitis.

A New Method for Total Fat Detection in Raw Milk Based on Dual Low-Coherence Interferometer

The present work experimentally demonstrates a multimode fiber optic sensing setup for total fat detection in raw milk samples. The optical fiber arrangement incorporates a low-coherence Fabry-Perot cavity operating in dual response. The system provides a phase modulation for a total fat range from 0.97% to 4.36%. Here, the protein remains constant at 3%. The data indicate that maximum sensitivity close to 616 pm/%fat could be achieved at optimal wavelength operation (500 nm). In addition, the system presented a minimal repeatability error measurement of 0.08%, cross-sensitivity between protein and fat of 0.134, and a regression coefficient of r2=0.9763. A thermal analysis was also performed, which indicate the temperature immunity of the system. The proposed method represents a low-cost alternative to detect minimal fat variations in raw cow milk.

Comparison of Bayesian and partial least squares regression methods for mid-infrared prediction of cheese-making properties in Montbéliarde cows

Assessing the cheese-making properties (CMP) of milks with a rapid and cost-effective method is of particular interest for the Protected Designation of Origin cheese sector. The aims of this study were to evaluate the potential of mid-infrared (MIR) spectra to estimate coagulation and acidification properties, as well as curd yield (CY) traits of Montbéliarde cow milk. Samples from 250 cows were collected in 216 commercial herds in Franche-Comté with the objectives to maximize the genetic diversity as well as the variation in milk composition. All coagulation and CY traits showed high variability (10 to 43%). Reference analyses performed for soft (SC) and pressed cooked (PCC) cheese technology were matched with MIR spectra. Prediction models were built on 446 informative wavelengths not tainted by the water absorbance, using different approaches such as partial least squares (PLS), uninformative variable elimination PLS, random forest PLS, Bayes A, Bayes B, Bayes C, and Bayes RR. We assessed equation performances for a set of 20 CMP traits (coagulation: 5 for SC and 4 for PCC; acidification: 5 for SC and 3 for PCC; laboratory CY: 3) by comparing prediction accuracies based on cross-validation. Overall, variable selection before PLS did not significantly improve the performances of the PLS regression, the prediction differences between Bayesian methods were negligible, and PLS models always outperformed Bayesian models. This was likely a result of the prior use of informative wavelengths of the MIR spectra. The best accuracies were obtained for curd yields expressed in dry matter (CY_DM) or fresh (CY_FRESH) and for coagulation traits (curd firmness for PCC and SC) using the PLS regression. Prediction models of other CMP traits were moderately to poorly accurate. Whatever the prediction methodology, the best results were always obtained for CY traits, probably because these traits are closely related to milk composition. The CY_DM predictions showed coefficient of determination (R²) values up to 0.92 and 0.87, and RS_y,x values of 3 and 4% for PLS and Bayes regressions, respectively. Finally, we divided the data set into calibration (2/3) and validation (1/3) sets and developed prediction models in external validation using PLS regression only. In conclusion, we confirmed, in the validation set, an excellent prediction for CY_DM [R² = 0.91, ratio of performance to deviation (RPD) = 3.39] and a very good prediction for CY_FRESH (R² = 0.84, RPD = 2.49), adequate for analytical purposes. We also obtained good results for both PCC and SC curd firmness traits (R² ≥ 0.70, RPD ≥1.8), which enable quantitative prediction.

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

Simple Summary

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

Abstract

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

The Strong Antioxidant Sheep/Goat Whey Protein Protects Against mTOR Overactivation in Rats: A Mode of Action Mimicking Fasting

Whey protein, a by-product of the cheese industry, can be putatively used as a functional food due to its beneficial health properties. The main objective of the present study was to assess in vivo the effect of a sheep/goat whey protein on the plasma amino acid profile and mammalian target of rapamycin (mTOR), a regulator of skeletal myogenesis. A control group was fed with a standard commercial diet while the experimental group received a standard commercial diet plus sheep/goat whey protein for 28 days. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) was conducted to determine plasma amino acid levels while the expression of p70-S6 Kinase 1 (p70-S6K1) in liver and quadriceps muscles was quantified and used as a biomarker of mTOR activity. The results obtained showed a decrease in the levels of essential and branched-chain amino acids (BCAAs) in the experimental group. Furthermore, p70-S6K1 expression was decreased in the liver of rats consumed whey protein. In conclusion, the reduction of amino acid levels and the concomitant inactivation of mTOR imply that whey could potentially act protectively against disorders induced by mTOR overactivation. Intriguingly, this mode of action mimics fasting, an approach with established advantageous health effects.

Predictive ability of mid-infrared spectroscopy for major mineral composition and coagulation traits of bovine milk by using the uninformative variable selection algorithm

Milk minerals and coagulation properties are important for both consumers and processors, and they can aid in increasing milk added value. However, large-scale monitoring of these traits is hampered by expensive and time-consuming reference analyses. The objective of the present study was to develop prediction models for major mineral contents (Ca, K, Mg, Na, and P) and milk coagulation properties (MCP: rennet coagulation time, curd-firming time, and curd firmness) using mid-infrared spectroscopy. Individual milk samples (n=923) of Holstein-Friesian, Brown Swiss, Alpine Grey, and Simmental cows were collected from single-breed herds between January and December 2014. Reference analysis for the determination of both mineral contents and MCP was undertaken with standardized methods. For each milk sample, the mid-infrared spectrum in the range from 900 to 5,000cm(-1) was stored. Prediction models were calibrated using partial least squares regression coupled with a wavenumber selection technique called uninformative variable elimination, to improve model accuracy, and validated both internally and externally. The average reduction of wavenumbers used in partial least squares regression was 80%, which was accompanied by an average increment of 20% of the explained variance in external validation. The proportion of explained variance in external validation was about 70% for P, K, Ca, and Mg, and it was lower (40%) for Na. Milk coagulation properties prediction models explained between 54% (rennet coagulation time) and 56% (curd-firming time) of the total variance in external validation. The ratio of standard deviation of each trait to the respective root mean square error of prediction, which is an indicator of the predictive ability of an equation, suggested that the developed models might be effective for screening and collection of milk minerals and coagulation properties at the population level. Although prediction equations were not accurate enough to be proposed for analytic purposes, mid-infrared spectroscopy predictions could be evaluated as phenotypic information to genetically improve milk minerals and MCP on a large scale.