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

Towards on-line cheese monitoring: Exploration of semi-hard cheeses using NIR and 1H NMR spectroscopy

This study aims to investigate the potential of using advanced spectroscopies for cheese quality monitoring. For this purpose, six semi-hard cheeses manufactured using lactic acid bacteria (LAB) and/or propionic acid bacteria (PAB) were explored using near-infrared spectroscopy (NIRS) and Proton Nuclear Magnetic Resonance (1H NMR) spectroscopy. The spectral data were analyzed using principal component analysis for extraction of possible discriminative patterns in quality parameters. The results show that the green analytical, but primarily bulk-sensitive, NIRS method was able to discriminate the cheese varieties primarily due to differences in the first overtone CH stretching region between 1650 and 1720 nm, in particular by the lactate methylene absorption at 1674 nm. A total of 25 metabolites were identified in the 1H NMR spectra of the cheese extracts, several of which were associated with the LAB and PAB metabolic pathways. PAB-associated metabolites include propionate, acetate, and glutamate, while LAB-associated metabolites include lactate and acetoin among others.

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.

Interrelationships among physical and chemical traits of cheese: Explanatory latent factors and clustering of 37 categories of cheeses

Cheese presents extensive variability in physical, chemical, and sensory characteristics according to the variety of processing methods and conditions used to create it. Relationships between the many characteristics of cheeses are known for single cheese types or by comparing a few of them, but not for a large number of cheese types. This case study used the properties recorded on 1,050 different cheeses from 107 producers grouped into 37 categories to analyze and quantify the interrelationships among the chemical and physical properties of many cheese types. The 15 cheese traits considered were ripening length, weight, firmness, adhesiveness, 6 different chemical characteristics, and 5 different color traits. As the 105 correlations between the 15 cheese traits were highly variable, a multivariate analysis was carried out. Four latent explanatory factors were extracted, representing 86% of the covariance matrix: the first factor (38% of covariance) was named Solids because it is mainly linked positively to fat, protein, water-soluble nitrogen, ash, firmness, adhesiveness, and ripening length, and negatively to moisture and lightness; the second factor (24%) was named Hue because it is linked positively to redness/blueness, yellowness/greenness, and chroma, and negatively to hue; the third factor (17%) was named Size because it is linked positively to weight, ripening length, firmness, and protein; and the fourth factor (7%) was named Basicity because it is linked positively to pH. The 37 cheese categories were grouped into 8 clusters and described using the latent factors: the Grana Padano cluster (characterized mainly by high Size scores); hard mountain cheeses (mainly high Solids scores); very soft cheeses (low Solids scores); blue cheeses (high Basicity scores), yellowish cheeses (high Hue scores), and 3 other clusters (soft cheeses, pasta filata and treated rind, and firm mountain cheeses) according to specific combinations of intermediate latent factors and cheese traits. In this case study, the high variability and interdependence of 15 major cheese traits can be substantially explained by only 4 latent factors, allowing us to identify and characterize 8 cheese type clusters.

Exploring the use of NIR and Raman spectroscopy for the prediction of quality traits in PDO cheeses

The aims of this proof of principle study were to compare two different chemometric approaches using a Bayesian method, Partial Least Square (PLS) and PLS-discriminant analysis (DA), for the prediction of the chemical composition and texture properties of the Grana Padano (GP) and Parmigiano Reggiano (PR) PDO cheeses by using NIR and Raman spectra and quantify their ability to distinguish between the two PDO and among their ripening periods. For each dairy chain consortium, 9 cheese samples from 3 dairy industries were collected for a total of 18 cheese samples. Three seasoning times were chosen for each dairy industry: 12, 20, and 36 months for GP and 12, 24, and 36 months for PR. A portable NIR instrument (spectral range: 950-1,650 nm) was used on 3 selected spots on the paste of each cheese sample, for a total of 54 spectra collected. An Alpha300 R confocal Raman microscope was used to collect 10 individual spectra for each cheese sample in each spot for a total of 540 Raman spectra collected. After the detection of eventual outliers, the spectra were also concatenated together (NIR + Raman). All the cheese samples were assessed in terms of chemical composition and texture properties following the official reference methods. A Bayesian approach and PLS-DA were applied to the NIR, Raman, and fused spectra to predict the PDO type and seasoning time. The PLS-DA reached the best performances, with 100% correctly identified PDO type using Raman only. The fusion of the data improved the results in 60% of the cases with the Bayesian and of 40% with the PLS-DA approach. A Bayesian approach and a PLS procedure were applied to the NIR, Raman, and fused spectra to predict the chemical composition of the cheese samples and their texture properties. In this case, the best performance in validation was reached with the Bayesian method on Raman spectra for fat (R2VAL = 0.74). The fusion of the data was not always helpful in improving the prediction accuracy. Given the limitations associated with our sample set, future studies will expand the sample size and incorporate diverse PDO cheeses.

Portable NIR spectroscopy: the route to green analytical chemistry

There is a growing interest for cost-effective and nondestructive analytical techniques in both research and application fields. The growing approach by near-infrared spectroscopy (NIRs) pushes to develop handheld devices devoted to be easily applied for in situ determinations. Consequently, portable NIR spectrometers actually result definitively recognized as powerful instruments, able to perform nondestructive, online, or in situ analyses, and useful tools characterized by increasingly smaller size, lower cost, higher robustness, easy-to-use by operator, portable and with ergonomic profile. Chemometrics play a fundamental role to obtain useful and meaningful results from NIR spectra. In this review, portable NIRs applications, published in the period 2019-2022, have been selected to indicate starting references. These publications have been chosen among the many examples of the most recent applications to demonstrate the potential of this analytical approach which, not having the need for extraction processes or any other pre-treatment of the sample under examination, can be considered the "true green analytical chemistry" which allows the analysis where the sample to be characterized is located. In the case of industrial processes or plant or animal samples, it is even possible to follow the variation or evolution of fundamental parameters over time. Publications of specific applications in this field continuously appear in the literature, often in unfamiliar journal or in dedicated special issues. This review aims to give starting references, sometimes not easy to be found.

Hyperspectral imaging through vacuum packaging for monitoring cheese biochemical transformation caused by Clostridium metabolism

This study developed a novel method for monitoring cheese contamination with Clostridium spores non-invasively using hyperspectral imaging (HSI). The ability of HSI to quantify Clostridium metabolites was investigated with control cheese and cheese manufactured with milk contaminated with Clostridium tyrobutyricum, Clostridium butyricum and Clostridium sporogenes. Microbial count, HSI and SPME-GC-MS data were obtained over 10 weeks of storage. The developed method using HSI successfully quantified butyric acid (R² = 0.91, RPD = 3.38) a major compound of Clostridium metabolism in cheese. This study creates a new venue to monitor the spatial and temporal development of late blowing defect (LBD) in cheese using fast and non-invasive measurement.

Potential of Near-Infrared Spectroscopy (NIRS) for Efficient Classification Based on Postharvest Storage Time, Cultivar and Maturity in Coconut Water

Coconut water (CW) is a popular and healthful beverage, and ensuring its quality is crucial for consumer satisfaction. This study aimed to explore the potential of near-infrared spectroscopy (NIRS) and chemometric methods for analyzing CW quality and distinguishing samples based on postharvest storage time, cultivar, and maturity. CW from nuts of Wenye No. 2 and Wenye No. 4 cultivars in China, with varying postharvest storage time and maturities, were subjected to NIRS analysis. Partial least squares regression (PLSR) models were developed to predict reducing sugar and soluble sugar contents, revealing moderate applicability but lacking accuracy, with the residual prediction deviation (RPD) values ranging from 1.54 to 1.83. Models for TSS, pH, and TSS/pH exhibited poor performance with RPD values below 1.4, indicating limited predictability. However, the study achieved a total correct classification rate exceeding 95% through orthogonal partial least squares discriminant analysis (OPLS-DA) models, effectively discriminating CW samples based on postharvest storage time, cultivar, and maturity. These findings highlight the potential of NIRS combined with appropriate chemometric methods as a valuable tool for analyzing CW quality and efficiently distinguishing samples. NIRS and chemometric techniques enhance quality control in coconut water, ensuring consumer satisfaction and product integrity.

How Chemometrics Can Fight Milk Adulteration

Adulteration and fraud are amongst the wrong practices followed nowadays due to the attitude of some people to gain more money or their tendency to mislead consumers. Obviously, the industry follows stringent controls and methodologies in order to protect consumers as well as the origin of the food products, and investment in these technologies is highly critical. In this context, chemometric techniques proved to be very efficient in detecting and even quantifying the number of substances used as adulterants. The extraction of relevant information from different kinds of data is a crucial feature to achieve this aim. However, these techniques are not always used properly. In fact, training is important along with investment in these technologies in order to cope effectively and not only reduce fraud but also advertise the geographical origin of the various food and drink products. The aim of this paper is to present an overview of the different chemometric techniques (from clustering to classification and regression applied to several analytical data) along with spectroscopy, chromatography, electrochemical sensors, and other on-site detection devices in the battle against milk adulteration. Moreover, the steps which should be followed to develop a chemometric model to face adulteration issues are carefully presented with the required critical discussion.

Major Causes of Variation of External Appearance, Chemical Composition, Texture, and Color Traits of 37 Categories of Cheeses

Cheeses are produced by many different procedures, giving rise to many types differing in ripening time, size, shape, chemical composition, color, texture, and sensory properties. As the first step in a large project, our aim was to characterize and quantify the major sources of variation in cheese characteristics by sampling 1050 different cheeses manufactured by over 100 producers and grouped into 37 categories (16 with protected designation of origin, 4 traditional cheese categories, 3 pasta filata cheese categories, 5 flavored cheese categories, 2 goat milk categories, and 7 other categories ranging from very fresh to very hard cheeses). We obtained 17 traits from each cheese (shape, height, diameter, weight, moisture, fat, protein, water soluble nitrogen, ash, pH, 5 color traits, firmness, and adhesiveness). The main groups of cheese categories were characterized and are discussed in terms of the effects of the prevalent area of origin/feeding system, species of lactating females, main cheese-making technologies, and additives used. The results will allow us to proceed with the further steps, which will address the interrelationships among the different traits characterizing cheeses, detailed analyses of the nutrients affecting human health and sensorial fingerprinting.

Genomic evaluation methods to include intermediate correlated features such as high-throughput or omics phenotypes

Gene expression is supposed to be an intermediate between DNA and the phenotype, and it can be measured. Thus, for a trait, we may have intermediate measures, which are in fact a series of genetically controlled traits. Similarly, several traits may be measured or predicted using infrared spectra, accelerometers, and similar high-throughput measures that we will call "omics." Although these measurements have errors, many of them are heritable, and they may be more accurate or easier to record than the trait of interest. It is therefore important to develop methods to use intermediate measurements in selection. Here, we present methods and perspectives for selection based on massively recorded intermediate traits (omics). Recent developments allow a hierarchical integrated framework for prediction, in which a trait is partially controlled by omics. In addition, the omics measures are themselves partly controlled by genetics ("mediated breeding values") and partly by environment or residual factors. Thus, a part of the genetic determinism of a trait is mediated by omics, whereas the remaining part is not mediated, which results in "residual breeding values." In such a framework, genetic evaluations consist of 2 nested genomic BLUP-based models. In the first, the effect of omics on the trait (which can be seen as an improved estimate of the phenotype) and the residual breeding values are estimated. The second model extracts the mediated breeding values from the improved estimate of the phenotype, considering that omics themselves are heritable. The whole procedure is called GOBLUP (genomics omics BLUP) and it allows measures in only some individuals; that is, it is a "single-step"-like method. In this model, heritability is split into "mediated" and "not mediated" parts. This decomposition allows us to predict how accurate the omics measure of the trait would be compared with the direct measure. The ideal omics measure is heritable and explains a large part of the phenotypic variation of the trait. Ideally, this could be the case for some traits with low heritability. However, even if the omics measure explains only a small part of the phenotypic variation, when omics measurement themselves are heritable, the use of such a model would lead to more accurate selection. Expressions for upper bounds of reliability given omics measurements are also presented. More studies are needed to confirm the usefulness of omics or high-throughput prediction. Usefulness of the technology likely needs to be checked on a case-by-case basis.

Authenticity and Typicity of Traditional Cheeses: A Review on Geographical Origin Authentication Methods

Food fraud, corresponding to any intentional action to deceive purchasers and gain an undue economical advantage, is estimated to result in a 10 to 65 billion US dollars/year economical cost worldwide. Dairy products, such as cheese, in particular cheeses with protected land- and tradition-related labels, have been listed as among the most impacted as consumers are ready to pay a premium price for traditional and typical products. In this context, efficient food authentication methods are needed to counteract current and emerging frauds. This review reports the available authentication methods, either chemical, physical, or DNA-based methods, currently used for origin authentication, highlighting their principle, reported application to cheese geographical origin authentication, performance, and respective advantages and limits. Isotope and elemental fingerprinting showed consistent accuracy in origin authentication. Other chemical and physical methods, such as near-infrared spectroscopy and nuclear magnetic resonance, require more studies and larger sampling to assess their discriminative power. Emerging DNA-based methods, such as metabarcoding, showed good potential for origin authentication. However, metagenomics, providing a more in-depth view of the cheese microbiota (up to the strain level), but also the combination of methods relying on different targets, can be of interest for this field.

Non-Invasive Method to Predict the Composition of Requeijão Cremoso Directly in Commercial Packages Using Time Domain NMR Relaxometry and Chemometrics

Low Field Time-Domain Nuclear Magnetic Resonance (TD-NMR) relaxometry was used to determine moisture, fat, and defatted dry matter contents in “requeijão cremoso” (RC) processed cheese directly in commercial packaged (plastic cups or tubes with approximately 200 g). Forty-five samples of commercial RC types (traditional, light, lactose-free, vegan, and fiber) were analyzed using longitudinal (T₁) and transverse (T₂) relaxation measurements in a wide bore Halbach magnet (0.23 T) with a 100 mm probe. The T₁ and T₂ analyses were performed using CWFP-T₁ (Continuous Wave Free Precession) and CPMG (Carr-Purcell-Meiboom-Gill) single shot pulses. The scores of the principal component analysis (PCA) of CWFP-T₁ and CPMG signals did not show clustering related to the RC types. Optimization by variable selection was carried out with ordered predictors selection (OPS), providing simpler and predictive partial least squares (PLS) calibration models. The best results were obtained with CWFP-T₁ data, with root-mean-square errors of prediction (RMSEP) of 1.38, 4.71, 3.28, and 3.00% for defatted dry mass, fat in the dry and wet matter, and moisture, respectively. Therefore, CWFP-T₁ data modeled with chemometrics can be a fast method to monitor the quality of RC directly in commercial packages.