Metabolomic profiling of Fiore Sardo cheese: Investigation of the influence of thermal treatment and ripening time using univariate and multivariate classification techniques.

The effect of different sub-pasteurization heat treatments and different ripening times was investigated in this work. The metabolite profiles of 95 cheese samples were analyzed using GC-MS in order to determine the effects of thermal treatment (raw milk, 57 °C and 68 °C milk thermization) and ripening time (105 and 180 days). ANOVA test on GC-MS peaks complemented with false discovery rate correction was employed to identify the compounds whose levels significantly varied over different ripening times and thermal treatments. The univariate t-test classifier and Partial Least Square Discriminant Analysis (PLS-DA) provided acceptable classification results, with an overall accuracy in cross-validation of 76% for the univariate model and 72% from the PLS-DA. The metabolites that mostly changed with ripening time were amino acids and one endocannabinoid (i.e., arachidonoyl amide), while compounds belonging to the classes of biogenic amines and saccharides resulted in being strongly affected by the thermization process.

Investigation of Seasonal Variation in Fatty Acid and Mineral Concentrations of Pecorino Romano PDO Cheese: Imputation of Missing Values for Enhanced Classification and Metabolic Profile Reconstruction.

Seasonal variation in fatty acids and minerals concentrations was investigated through the analysis of Pecorino Romano cheese samples collected in January, April, and June. A fraction of samples contained missing values in their fatty acid profiles. Probabilistic principal component analysis, coupled with Linear Discriminant Analysis, was employed to classify cheese samples on a production season basis while accounting for missing data and quantifying the missing fatty acid concentrations for the samples in which they were absent. The levels of rumenic acid, vaccenic acid, and omega-3 compounds were positively correlated with the spring season, while the length of the saturated fatty acids increased throughout the production seasons. Concerning the classification performances, the optimal number of principal components (i.e., 5) achieved an accuracy in cross-validation equal to 98%. Then, when the model was tasked with imputing the lacking fatty acid concentration values, the optimal number of principal components resulted in an R2 value in cross-validation of 99.53%.