Next Article in Journal
Nutritional Properties of Camelids and Equids Fresh and Fermented Milk
Next Article in Special Issue
GC-MS Metabolomics and Antifungal Characteristics of Autochthonous Lactobacillus Strains
Previous Article in Journal
New Treatment Option for Clinical and Subclinical Mastitis in Dairy Cows Using Acoustic Pulse Technology (APT)
Previous Article in Special Issue
LC-QTOF/MS Untargeted Metabolomics of Sheep Milk under Cocoa Husks Enriched Diet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quality Control in Fiore Sardo PDO Cheese: Detection of Heat Treatment Application and Production Chain by MRI Relaxometry and Image Analysis

Porto Conte Ricerche s.r.l., S.P. 55 Porto Conte-Capo Caccia, Km 8.400 Loc. Tramariglio, 07041 Alghero, SS, Italy
*
Author to whom correspondence should be addressed.
Submission received: 4 March 2021 / Revised: 14 April 2021 / Accepted: 13 May 2021 / Published: 26 May 2021
(This article belongs to the Special Issue Innovation Meets Tradition in the Sheep and Goat Dairy Industry)

Abstract

:
Fiore Sardo (FS), a traditional Italian cheese, is present in the market as a heterogeneous variety of products. The use of heat-treated (HT) milk is forbidden by the official production protocol, but no official analytical method able to detect heat application is yet available. Here, a combined magnetic resonance imaging (MRI) relaxometry and image analysis approach to recognize FS made from raw milk is presented. Artisanal FS cheeses were produced from raw milk (RC) by five shepherds in accordance with the official protocol. They were compared to HT-milk counterparts (HTC). Additionally, industrially manufactured commercial FS cheeses (I) were also purchased and compared to RC and HTC. Relaxometry data of FS indicated the presence of two water populations; the ratio of characteristic relaxation time constant T2 and area fraction (Score, Ṩ) of the fastest relaxing population was used to compare RC, HTC and I samples. RC from HTC were successfully discriminated, the latter exhibiting lower Ṩ (enhanced protein hydration). I cheeses exhibited the lowest Ṩ values, sometimes comparable to HTC. Since visual appearance of RC and HTC is appreciably different, an image analysis deep learning approach using MRI and photographic pictures was adopted to discriminate the two productions, with promising percentages (>93%).

Graphical Abstract

1. Introduction

Protected Designation of Origin (PDO) describes a group of agricultural/food products originating from a specific geographical area, for which production, processing and preparation must take place according to a recognised know-how [1].
In the European Union, 185 cheeses are registered as PDO products; this number rises to 231 if cheeses registered as Protected Geographical Indications (PGI) are considered [2]. A first classification of these products can be made based on milk processing. A total of 91 cheeses out of 231 (PDO + PGI) compulsorily require the use of raw milk; in other cases, the use of both raw and pasteurized/heat-treated milk is allowed. The lack of a specific requirement allows producers to freely decide how to process milk, leading to cheeses characterized by very different attributes [3].
The effects of milk heat treatment on cheese features have been extensively investigated and explained. From a microbiological and sensory point of view, the heat treatment depletes pathogens along with the characteristic microbiota, resulting in a more uniform final product that exhibits a less characteristic flavour profile compared to the raw milk counterpart [4]. From a technological point of view, though the application of a heat treatment guarantees improved hygienic conditions during processing and a better safety of the final product, it can result in an impaired coagulation and syneresis and in an alteration of cheese rheology and texture [5,6].
Fiore Sardo PDO (FS) cheese is the oldest cheese of Sardinia (Italy), being historically produced by shepherds in very small artisanal cheese factories. It must be obtained exclusively using raw whole milk from the Sarda Sheep breed. The official cheesemaking protocol still follows the ancient traditional process [7]. The peculiarities of the manufacturing process confer to Fiore Sardo a firm, crumbly and grainy texture, which easily breaks into flakes, and a savoury, piquant and smoky flavour [8,9,10,11]. Market offer is characterized by a quite heterogeneous variety of Fiore Sardo PDO, with very different quality attributes, with almost 90% of the product supplied by the dairy industry (Agris Sardegna, 2015).
In particular, the use of heat-treated milk for FS production has been associated with common industrial processing practices [12,13,14,15]. Previous studies concluded that this practice cannot be excluded [16] and that it is not unrealistic that some FS are manufactured using heat-treated milk, against the specifications indicated in the EC Regulation. From a scientific point of view, several chemical and physical characteristics of some productions were found to be reasonably explained by assuming the effects of thermal treatments on milk [16]. This is also supported by the fact that no official analytical methods able to discriminate between ewe milk cheese obtained from raw and thermised milk are available [17].
Some analytical techniques able to discriminate products made from raw milk and heat-treated milk have been proposed for product quality safeguards [17,18,19,20,21].
Analytical approaches based on nuclear magnetic resonance (NMR) relaxometry are particularly interesting. This technique was previously applied to describe the thermal denaturation of whey proteins and other milk components. Lambelet et al. [22], for example, demonstrated that irreversible thermal denaturation of whey proteins can be detected by low-field NMR. They observed that when NMR transverse relaxation rates were plotted versus heating temperature, a sigmoid curve was obtained, with relaxation rate values showing a steep increase from about 60 °C to about 80 °C, and a flex corresponding to milk thermisation/pasteurization temperatures usually adopted in dairy industrial processing (64 °C to 72 °C range). Later work by the same research group extended the results to other milk components and dyalized skim milk [23], suggesting the suitability of NMR relaxometry for measuring the thermal denaturation of milk proteins. Recent works, based on a high sample numerosity and exploiting an experimental setup more closely resembling real industrial processing conditions, demonstrated that the effect of heat treatment on milk is transferred to NMR relaxometry characteristics of both rennet curd and aged cheese. Raw and heat-treated milk [24] and derived curds [25,26] could be successfully discriminated by benchtop low-field NMR relaxometry. A similar discriminative ability of NMR relaxometry was also observed in ewe milk mature cheese by high-field magnetic resonance imaging (MRI) [19].
The aforementioned investigations demonstrated the sensitivity of NMR relaxation parameters to key cheesemaking processes (mainly thermal treatments of milk) in real dairy systems (whole milk, curds, mature cheese), suggesting the suitability of NMR relaxometry as a valid and useful tool in quality control activities in the dairy industry and for the safeguard of typical dairy productions.
In this study, heat treatments of milk were intentionally applied to produce FS cheese, and NMR relaxometry of raw and heat-treated FS were compared. Moreover, a comparison between artisanal and some industrial FS, the latter purchased from the market and visibly different from the typical FS PDO, was carried out for a more comprehensive evaluation of the effect of production chain on relaxometry features and aiming to optimize the exploitation of NMR relaxometry to assess the adherence to the official production protocol. MRI was used to compare heat-treated samples to their raw counterparts. A supplemental investigation was also carried out by exploiting image analysis protocols on a selected subset of MRI images and digital pictures in order to correlate process-induced changes with MRI image features and the visual appearance of samples.

2. Materials and Methods

Two sets of experiments were carried out, as summarized in Table S1. The first dataset (Dataset 1) was acquired by analysing samples made by different artisanal producers in different seasons. Additionally, a selection of industrial cheeses was purchased from the local market and subjected to the same analytical protocol used to characterize samples from Dataset 1.

2.1. Production of Fiore Sardo (FS) Cheese

2.1.1. Dataset 1

Five shepherds (S1, S2, S3, S4, S5) who were artisanal manufacturers of Fiore Sardo were selected for cheese sample production. Each shepherd produced a total of 20 cheese wheels in two different seasons (March–April 2019, Season 1; January–February 2020, Season 2; 100 cheese wheels in total for the two seasons). Ten wheels out of 20 were produced from raw milk, while the remaining ten were produced from the same milk after thermisation. Milk thermisation was carried out in different industrial plants, using plate heat exchanger modules. Cheesemaking was performed according to the official cheesemaking protocol [7,9], but milk thermisation was used for selected samples. Aged cheese samples were collected after 105 and 180 days of ripening in each shepherd’s cellar. One quarter portion of each cheese wheel was sampled at each of the abovementioned time points, for a total of 10 wheels per producer per ripening time, 5 of which were made from raw milk (RC) and 5 from heat-treated milk (HTC).

2.1.2. Dataset 2

Industrial FS cheeses with similar ripening (6–8 months) were purchased from the local market in different seasons (January–November). Both cheeses manufactured by industrial dairy factories (11 samples) and from a maturer industry (2 samples) were purchased. Maturer industries buy Fiore Sardo from artisan producers and ripen these cheeses in their industrial cellars for several months, until they can be released to the market, i.e., after at least 105 days of ripening, according to the official cheesemaking protocol. Fiore Sardo PDO cheeses from the maturer industry are sold with the industrial brands, but milk is processed by shepherds and small artisans according to traditional processes [9]. A schematic representation of sample production and experimental plans (relative to both Dataset 1 and Dataset 2) is reported in Figure S1 (Supplementary Material).

2.2. Magnetic Resonance Imaging (MRI) Analysis

Each analytical sample of cheese subjected to MRI analysis represented the most central part of the cheese wheel, as schematically depicted in Figure S2.
In brief, each quarter was carefully cut in the portion corresponding to the wheel centre to obtain one cylindrical sample (diameter 2.3 cm, height 3–4 cm) to easily fit into a 50 mL Falcon tube. All samples were equilibrated to 22 °C for 1 h before MRI analysis. MRI experiments were performed using a Bruker Avance 300 MHz equipped with a microimaging Micro2.5 probe (Bruker Biospin, Karlsruhe, Germany) at room temperature (approximately 22–24 °C). A conventional Carr Purcell Meiboom Gill (CPMG) spin echo sequence was used with the following parameters: single 1 mm slice; echo time (TE) = 7.907 ms; repetition time (TR) = 3 s; echoes = 64; matrix = 128 × 128; number of excitations = 1–3; acquisition time = 25 min. To selectively analyse signals from water, the fat signal was suppressed, taking advantage of the different chemical shifts of fat and water at 7.05 T, with a 90° selective pulse applied before the conventional CPMG spin echo sequence at a frequency offset of 3.5 ppm with respect to water (preparatory fat saturation scheme provided by Bruker Topspin). All MR images were characterized by a signal to noise ratio higher than 100.
The AnalyzeNNLS software package was used to deconvolute the multi-exponential transversal relaxation time (T2) decay of water proton signals [27]. Bruker MRI image data were converted into multiecho image data (MEID) files (Matlab, R2020b, The Mathworks, Sherborn, MA). Regions of interest (ROIs) were manually drawn on the MEID file (five ROIs for each image) and automatically processed to obtain the geometric mean T2 (T2 relaxation time constant, corresponding to the maximum intensity of the peak, on a logarithmic scale) and the area fraction (obtained by dividing the sum of the T2 distribution within a desired region between T2 min and T2 max by the sum of all T2 distributions).

2.3. Moisture Content

Moisture content was measured using the rapid determination method for cheese (130 °C, 90 min; [28]).

2.4. Statistical Analysis

The effect of thermisation was evaluated by comparing MRI parameters (geometric mean T2 and area fractions) between RC and HTC by using univariate Student’s t-test. An additional statistical comparison (one-way ANOVA, followed by Tukey post-hoc test) was performed to compare each shepherd’s (i.e., artisanal) production to the industrial and maturer industry cheeses. Industrial cheeses (11 samples) were considered as a unique group of data, as were maturer cheeses (2 samples). Both tests were carried out with the MetaboAnalyst tool (https://www.metaboanalyst.ca, v.5, accessed 20 November 2020) [29], considering a significance threshold limit of p < 0. 05. Box plots were obtained using MetaboAnalyst.

2.5. MRI Data Conversion

T2 MRI Bruker Paravision raw 2dseq data from full Dataset 1 were first converted to Neuroimaging Informatics Technology Initiative (NIfTI, .nii) format using Bru2nii software (Bruker2NIfTI v1.0.20180303: by Matthew Brett, Andrew Janke, Mikaël Naveau, Chris Rorden, Windows 64-bit) with a resizing factor of 10. Then, the 64 slices from each NIfTI file were converted in JPEG format using a Matlab script first developed by Alex Laurence [30]. Only MRI images derived from the first MSME slice (TE = 7.907 ms) were used for deep transfer learning classification purposes.

2.6. Photographic Acquisition and Processing of Cheese Paste Surfaces

Snapshots of the surface of Fiore Sardo cheese from Dataset 1 produced in Season 1 at 105 days of ripening were taken using a Machine Vision System (MVS) assembled in our laboratory. Briefly, our MVS consisted of a portable plastified (PVC) white light-box (24 × 23 × 22 cm) equipped with a USB LED strip of white-light-emitting diode lamp (input 5 V, lumen 550, colour temperature 6500 K), assembled on a metallic support and directed upwards at a 90° angle with a cardboard box wall (Figure S3). Photos were taken without flash using a 25-megapixel-wide camera integrated in a Galaxy A50 smartphone. Cheese quarters were placed under a dark background and manually adjusted to be below the camera with a fixed distance of 25 cm. All operations were carried out in a dark room to minimize the effect of outside light. The acquired images (Joint Photographic Expert Group—JPEG; 24-bit colour depth, spatial resolution 3072 × 2048 pixels) were stored on a removable memory card. Images were finally processed by removing black background and manually cropped by using in-house-modified scripts available in the Matlab Image Processing Toolbox (Matlab, R2020b, The Mathworks, Sherborn, MA, USA).
Both digital pictures and MRI images were reordered in a concatenated folder structure according to the scheme illustrated in Figure S4. Each subfolder, related to the two ripening stages for each season, was separately compressed in zip format as input data for the deep learning classification described in the following section.

2.7. Deep Transfer Learning Based Classification

The proposed classification method relies on the deep transfer learning approach based on deep neural network (DNN) architectures. DNNs consist of a series of functions that take an input and return a predicted label. For this purpose, labeled input data were initially split into a training dataset (i.e., the set of data used to fit the model) and a validation dataset (i.e., the set of data used to provide an unbiased evaluation of a model fit on the training dataset while optimizing, or tuning, the model parameters). In brief, labeled input data pass through the DNN network using a loss function to evaluate how well the model performs at correctly identifying the true classes. The model optimizes for this loss function by computing the gradient of the loss function with respect to the model parameters. In parallel, the validation process optimizes the model parameters iteratively to minimize the loss and perform validation of the trained model [31]. Six pre-trained DNN frameworks were tested: ResNet [32], AlexNet [33], VggNet [34], SqueezeNet [35], DenseNet [36] and Inception-Net [37]. All models were implemented using the learning framework Pytorch, written in Python v. 3.7 and integrated into the open source GPU web server Google Colaboratory [38]. For the training step, all six DNN frameworks were trained under the following conditions: number of epoch = 1000; batch size = 8; learning rate = 0.001; momentum = 0.9. When the loss function converged and stabilized, the training step was automatically stopped and the training model saved. For the validation set, 20% of the training dataset data was used, and a 5-fold cross validation was performed to evaluate the performance of each model. Each model was trained 5 times on each dataset, and related performance parameters (training and validation accuracy and loss as well as training time) were computed. Finally, for all pre-trained models, the average values of the five replicates were calculated and listed in tables in order to compare classification performance.

3. Results and Discussion

3.1. MRI Analysis of Dataset 1

NMR relaxation of protons can occur only if stimulated by a fluctuating field of proper frequency that would induce the spin transition, in a non-spontaneous relaxation process that finally leads to equilibrium magnetization. Longitudinal relaxation, defined by the constant T1, is triggered by field fluctuations at the Larmor frequency. Contrarily, in solid or semi-solid systems such as cheese, the major contribution of transverse relaxation, described by T2, comes from the slow molecular dynamics of the studied liquid at zero Larmor frequency. For this reason, observed T1 values usually present a very broad distribution, while populations are sharper and generally better defined for T2. We chose the latter NMR relaxometry constant to describe the analysed cheeses.
Informative relaxometry features of the samples can be derived by analysing the multiexponential decay of the transverse magnetization from CPMG experiments at variable interpulse spacing [19,27,39]. Fitting of CPMG decays of MRI images of FS cheeses indicated a multiexponential behaviour (i.e., two T2 populations), deriving from the presence of two water proton pools (Figure 1).
Multiexponential decay of the transverse component of NMR magnetization arises from the combined effect of diffusive and chemical exchange phenomena. In this condition, water molecules in cheese could be described, at a first approximation, as either free water or hydration water molecules (i.e., water embedded in the protein matrix and strongly interacting with it). Of course, the term “free” could be misleading, since water dynamics in this fraction are in fact much slower than in bulk water. NMR relaxometry is able to clearly detect the presence of different water pools, since water molecules sample different domains within the time scale of the CPMG experiment (i.e., the interpulse spacing). When diffusion is slow (and/or when relaxation domains are large), water molecules in the free water domain and hydration water retain their characteristic relaxation times, and multiexponential distribution is observed (e.g., in aged cheese); when diffusion is fast (and/or heterogeneity of domains is small) (e.g., in protein solutions, milk and curd), relaxation of water in different compartments is averaged, and proton transverse relaxation is mainly monoexponential [24,39,40].
At high field, diffusive exchange of water and fast chemical exchange between water protons, labile protein protons, and exchangeable protons in small molecules (such as sugars, vitamins, etc.) play a significant role in characterizing the observed relaxation distributions (Figure 1).
It was shown that in skim milk chemical exchange is mainly modulated by whey proteins [40], but in whole milk, the role of fat globules should be taken into account, since fat influences NMR relaxometry results [24]. In addition, it is clear that in cheese fat, protons might reasonably have a role in defining the observed NMR relaxometry profiles [41]. However, in the following discussion, the effect of fat protons on the observed T2 distributions can be reasonably neglected, since fat signal has been saturated by a solvent suppression pulse in the MRI sequence, as previously explained [19].
For each pool, the mean area fractions (AF1 and AF2) and the corresponding T2 values (i.e., the maximum intensity of each peak, T21 and T22) were obtained (Table 1).
These data are in agreement with previous studies on aged ewe milk cheese and FS PDO [19], which also reported a first T2 population relaxing in the range of 7–12 ms and a second one at 45–53 ms. The main novelty with respect to the mentioned data consists of the fact that the cheeses analysed in the present work were all made by following the official cheesemaking protocol of FS PDO, except for the thermal treatment of milk. Furthermore, a larger amount of data were used here to preliminarily validate observed data trends; the present results were indeed acquired during a much larger lag time (two production seasons), and two ripening times (105 and 180 days) were compared. Moreover, products from five different artisanal producers plus industrial Fiore Sardo PDO cheeses were analysed.
Statistical analysis of MRI relaxometry values (Table 1) indicated significant differences in T21 and T22 values between RC and HTC cheeses. Evaluation of AF suggested the same consistent trend, with a significantly larger amount of protons in the faster relaxing population (AF1) in HTC cheeses than in RC. This is in quite good agreement with previous results on MRI relaxometry of heat-treated milk cheeses [19] and on the sensorial evaluation of Fiore Sardo cheese, suggesting that sensorial attributes and relaxometric features may correlate to each other and that they are both more appreciably affected by milk properties and cheesemaking processes than ripening [11].
Since for each production reported in Table 1 (S1–S5) the same milk was used, the manufacturing process was carried out by the same dairyman, and ripening time was exactly the same for both treatments (RC and HTC), the only factor that accounts for the observed relaxometry differences is milk thermisation.
From a previous investigation [19] and in accordance with pertinent literature [42,43,44], the fast relaxing proton pool (described by the pair T21; AF1) can be associated with a fraction of water molecules strongly interacting with protons and entrapped in the para-casein network (protein hydration water); the slower population (T22; AF2) can be associated with more mobile water molecules, which less strongly interact with the cheese protein network (which is sometimes referred to as bulk or “free” water).
Several different mechanisms may reasonably affect observed transverse relaxation rates. The first mechanism is related to an altered proportion between more mobile (free) water pools and water in the vicinity of the protein surface (hydration water), which is reasonably well described by AF1 [16,19,45]. This mechanism is in turn modulated by the diffusive exchange rate between water entrapped within the protein network of cheese, confined in cheese holes or pores and hydration water, and chemical exchange [39,46] between water protons and labile protein protons (mainly whey proteins [40]). The second mechanism involves internal water molecules in casein micelles, sometimes referred to as “bound” water molecules, which represent structural water molecules in the internal cavities, engaged in slowly modulated intermolecular dipole couplings with protein protons [46]. A significant effect of these “bound” water molecules on the observed proton transverse relaxation rate has been already suggested in dairy systems [40]. A third mechanism is explained by considering changes in the correlation time of macromolecules as a function of the ripening and processing methods [16,47]. Finally, relevant contribution could be expected from the mineral equilibrium of milk and its effect on micelles hydration. Equilibrium between colloidal (solid) and soluble calcium phosphate is in fact affected by temperature [48,49]; this has also been observed in Fiore Sardo cheeses manufactured from raw or heat-treated milk [6]. In particular, higher total calcium was found in Fiore Sardo produced from raw milk than from thermised milk, and it was hypothesized that mineral equilibrium in milk was shifted from soluble calcium to insoluble colloidal calcium as an effect of a high treatment temperature, leading to a higher retention of calcium in Fiore Sardo from raw milk [6]. Micellar calcium phosphate is a highly hydrated colloid, which greatly influences proton transverse relaxation. Ca and P association to casein and micelle hydration are strictly related phenomena, having a marked effect of NMR transverse relaxation of protons [40,48,49]. Heat treatment of milk was in fact found to cause precipitation of calcium phosphate and correspondingly affect proton T2 [48].

3.2. MRI Analysis of Dataset 2

Industrial Fiore Sardo samples were characterized by the presence of two water proton pools, in agreement with results on artisanal samples RC and HTC discussed above. Mean area fractions (AF1 and AF2) and T2 values (T21 and T22) are shown in Table 2.
Both relaxation time constants T2 and area values AF were consistent with the HTC artisanal cheeses analysed in Dataset 1, with the exception of samples from the maturer industry (I12, I13).
Both the increase of AF1 and the decrease of T21 are the result of a strong interaction between water protons and the para-casein protein network of cheese, i.e., of cheese protein hydration. From a certain point of view, AF1 more directly describes the level of protein hydration, being determined by the amount of protons in this pool. However, a decrease of the T21 value, being the result of the weighted average between the (long) T2 of protons in free water pools and that (very short) of labile protein protons, characterized by a different chemical shift and affected by the typically long reorientational times of biological macromolecules, can be indirectly influenced by the amount of protein hydrating water.
Aiming to find suitable and objective descriptors of the relaxometric changes induced in cheese by milk thermisation, we define a new parameter, namely the Score (Ṩ) factor, as a representative indicator of the state of cheese microstructure. The Score (Ṩ) is defined as follows:
= T 21 / AF 1
where T21 and AF1 are the relaxation time and the area fraction of population 1, respectively, of the MRI T2 distribution profile (Figure 1). Based on what was discussed above, the lower the Score, the more affected cheese microstructure is by milk thermisation. Correspondingly, the higher the Score, the less influenced are cheese proteins by any thermal effect on milk, i.e., milk can be considered as raw or untreated.
The calculated Ṩ values (Season 1 and Season 2) are presented in Figure 2.
Even at first glance, two considerations can be made when observing Figure 2: first, thermisation of milk in HTC cheeses (red symbols) leads to a reduction of Ṩ with respect to their raw RC counterparts (blue symbols); second, RC samples show much higher Ṩ variability than their HTC counterparts.
Cheeses from one producer (S5) showed subtle, though significant, NMR relaxometry changes upon milk thermisation (blue and red circles are superimposed in Figure 2). In fact, for producer S5, both AF1 and T21 values always showed only slightly significant changes when RC and HTC cheese samples were compared at 105 days of ripening and at 180 days in Season 1 (Table 1). It should be considered that, at the shortest ripening time allowed for Fiore Sardo commercialization (105 days of ripening), the effect of water redistribution and protein hydration might still be not sufficient in order to clearly discriminate RC from HTC. In these cases, NMR relaxometry analysis should be replicated after a longer ripening time to obtain clearer results (Figure 2b,d). If a discrimination between RC and HTC is still not present at six months of ageing, then very likely milk has been overheated to some extent, and cheese features are not comparable to raw milk cheeses anymore. This result can be explained by assuming that, though thermisation was not intentionally performed in S5 production, raw milk was likely subjected to excessive heat before rennet addition. In fact, it is worth recalling here that, even at an artisanal level, milk is always treated to some extent during cheesemaking, since rennet is usually added when milk reaches, upon heating, a temperature in the range of approximately 35 and 38 °C, depending on the season, site, climate and usual practices of production. Among Fiore Sardo artisan producers, milk heating, according to traditional cheesemaking processes, is mostly carried out by directly warming up (often by means of an industrial stove gas cooker) a copper boiler containing raw milk. In some other cases, but still at artisanal level, more modern stainless steel multipurpose cheesemaking tanks are used to heat milk by means of water vapour in somewhat larger dairy plants. Only bigger industrial dairy plants adopt flow pasteurizers equipped with sophisticated plate heat exchanger modules, in order to perform more controlled high temperature–short time heat treatments. This consideration could explain the very subtle difference between the NMR relaxometry parameters of RC and HTC cheeses at 105 days and 180 days of ripening in S5 samples.
Notably, S5 and other shepherd’s productions show superimposed error bars in Figure 2. This is mainly due to the large variability of RC and to similar technological considerations as already discussed for S5.
As explained above, the score value Ṩ represents and describes the status of water molecules more tightly interacting with the paracasein network of cheese. Low score values (i.e., higher AF1 and/or lower T21) indicate a high degree of protein hydration, arising from either ripening [16,45] or from the effect of heat treatments [19]. For cheeses having the same ripening time, the effect of the thermal treatment of milk on the following processing steps [5] dominates relaxation.
In fact, previous studies carried out at variable fields by fast field cycling NMR relaxometry (FFC-NMR) have demonstrated that the hydration of cheese proteins is affected by cheese ageing [45,47]. Further FFC-NMR studies have demonstrated that, for a given ageing time, cheese protein hydration can also be significantly affected by the different processing methods preceding cheese maturation [16]. As discussed above, two similar Fiore Sardo cheeses, which only differ in their heat treatment of milk, also differ in their NMR relaxometry characteristics. The reasons for the observed relaxometry differences can be found in cheese microstructure (such as fractal dimension, protein proton to water proton ratio, water mobility and protein hydration) [16].
The variability of Score values (Ṩ) within each group and among groups of cheeses is presented in the form of box plots in Figure 3. Standard deviations of RC of shepherds (blue boxplot) were generally larger than the other groups. Median values of RC are higher than all other groups, followed by maturer cheeses, HTC and industrial samples, with the latter showing the lowest median values. RC samples show considerably higher interquartile range than HTC and much larger whiskers. This could be reasonably ascribed to both the effect of the artisanal manufacturing process on cheese microstructure and to the microbiological complexity of the systems arising from the use of raw milk. In fact, while artisanal cheesemaking entails mainly manual practices, industrial production implies more standardized protocols and a higher degree of process automation. In addition, previous reports on raw and heat-treated milk cheeses demonstrated a larger variability of raw-milk cheeses than pasteurized-milk cheese counterparts, and it was suggested that the effect of heat on milk proteins, indigenous enzymes and microbial flora of milk may reduce variability, resulting in a more homogeneous peptide profile, sensory properties and NMR relaxation parameters for pasteurized-milk cheese [19,50,51].
Differences in the manufacturing process are reflected in different physico-chemical and sensory features, resulting in the distinctiveness of each wheel from the same shepherd. While heating milk to 35–38 °C, necessary for a proper action of rennet, is usually not strictly standardized in terms of temperature control, industrial processing is carefully controlled. Mild thermal treatments of raw milk, typical of artisanal practices, allow for the preservation of the microorganisms that represent a sort of fingerprint of each shepherd’s production. The development of indigenous milk microflora and the related molecular changes (e.g., changes in pH, proteolysis and lipolysis during ripening, etc. [5]) occurring during the following processing steps, together with the exclusively manual manufacturing operations on curd, are reflected in the microstructural, molecular and sensory properties of the final cheese, making every wheel a unique product. That is why artisanal Fiore Sardo PDO cheeses are very different from each other, being strongly linked to the territory and to the specific production practices. It is interesting to note, in this regard, that milk and curd handling practices for making Fiore Sardo are passed down from generation to generation and often differ from each other in subtle, but likely relevant, details. On the other hand, when heat is applied to the same milk, though (almost) the same manufacturing practices are carried out by the same cheesemongers, part of the characteristic microflora is depleted and milk functionality is altered [51,52]. The Ṩ data distribution in HTC cheeses indeed indicated a reduced heterogeneity in comparison to RC. The box plot of industrial cheeses more clearly shows that the spreading of Ṩ values is considerably reduced and very narrow in comparison to both RC and HTC. This is mostly ascribable to industrial practices, where cheesemaking is strictly controlled and standardized in terms of both microbiological conditions and technological operations. It is worth noting that industrial Fiore Sardo samples were purchased from different industrial producers, so the observed low interquartile range should be considered quite representative of the standardization of the industrial production in the market and cannot be ascribed to a specific producer.
Mean and standard deviations of Score values Ṩ for all groups of samples (RC, HTC, I and maturer), together with their statistical significances, are presented in Table 3.
In general, ANOVA results suggested a consistent differentiation among the four groups of cheeses. Industrial Fiore Sardo showed the lowest Ṩ values in all comparisons, except for one case (i.e., S3 180 days, Season 1), which had Ṩ values higher than HTC only. RC often had significantly higher Ṩ values than other groups of cheeses, except for two cases (i.e., S4 and S5, 105 days, Season 1), when maturer cheeses showed higher Ṩ. This result can be explained by considering that maturer cheeses were purchased after approximately 6–8 months of ripening. In some other cases, Ṩ values were comparable for RC, HTC and maturer cheeses. Statistical comparison of HTC and industrial cheeses usually resulted in a significantly lower Ṩ for the latter, but they were comparable in five cases. Maturer samples were never comparable to industrial samples.
Based on the aforementioned statistical evaluations, some conclusions can be drawn on the manufacturing process of each group of samples and on the quality features of the resulting cheeses. RC and HTC were entirely produced following artisanal procedures. Maturer cheeses share the same artisanal processing but were ripened in industrial-like conditions, in refrigerated cellars with controlled humidity and temperature. Despite different ripening conditions, RC and maturer cheeses were described by very similar Ṩ values, clearly suggesting that technological transformations that precede ripening confer to cheeses a common relaxometric behaviour described by comparable Ṩ. RC and maturer cheeses were generally described by higher Ṩ than industrial cheeses, although ANOVA suggested that RC, HTC and maturer cheeses can still differ to some extent among each other in terms of the new quality-related parameter, Ṩ. Since the RC and HTC cheeses of each shepherd were technologically processed in the same way, the decrease in Ṩ observed between RC and HTC (Figure 3), supported by ANOVA comparisons (Table 3), can be reasonably and solely related to the application of a heat treatment of the milk. On the other hand, industrial samples, besides having a much narrower Ṩ variability (Figure 3) clearly exhibited considerably and significantly lower Ṩ values with respect to the other three groups of cheeses. These relaxometric behaviours could be reasonably associated with the features of industrial productions, for which milk and curd processing and handling are more mechanized and standardized, ensuring more uniform features in the final products, as previously reported [16]. Moreover, a contribution to the low Ṩ may arise from the application of a heat treatment to the milk that is consistent with common industrial processing practices, as realistically suggested [12,13,14,15]. Interestingly, the discrimination ability of Ṩ and its potential in identifying relevant quality features of Fiore Sardo cheeses does not seem to be affected by different seasons of lactation (i.e., to milk compositional quality), being instead more strictly related to milk processing. This consideration derives from the significant statistical differences among the four groups of samples, regardless the season (Figure 3, Table 3). As far as ripening conditions are concerned, the significant role of ripening cannot be excluded. However, the effect of ripening conditions seems to be less influential than the milk and curd processing practices preceding cheese ageing. This makes NMR relaxometry a suitable analytical tool for characterizing the effect of heat treatments of milk on quality features of Fiore Sardo cheese.

3.3. Moisture Content

The moisture content of all cheese samples are reported in Table S2.
Moisture content was not influenced by heat treatment, indeed showing no consistent trend in RC and HTC samples. This is in contrast to MRI parameters, which showed consistent changes according to heat treatments. Since MRI parameters T21 and AF1 could be considered representative of the portion of water in the samples, such changes could not be related to total moisture content [19]. In addition, for industrial samples, no noticeable correlation was found between MRI parameters and moisture content.
Water status can be described at different scales of investigation—at a macroscopic level by moisture content and at a molecular level by MRI—and can lead to different information about the analyzed matrix. Moisture content is representative of the extractable water molecules by drying in certain conditions and is among the most common parameters used to describe food quality. However, it should not be considered as an exhaustive parameter as it does not give any information on water proton pools in the system (if any) and their interactions with the protein matrix in cheese [53].

3.4. Image Analysis of Fiore Sardo Pictures and MRI Images

Representative pictures of the visual appearance of RC and HTC at the minimum allowed ripening time for Fiore Sardo PDO (i.e., 105 days) are shown in Figure 4.
Even at first sight, RC, HTC and industrial cheeses appear very different from each other. The HTC samples surface appears more homogeneous, smoother and with fewer holes, while RC cheeses are more grainy, friable and show more and bigger holes. When cut with a knife, RC cheeses are prone to easily break and crumble, while HTC appear more creamy and can be cut very easily and without breaking into pieces.
The different visual appearances of HTC samples with respect to RC cheeses is reasonably ascribable to the modifications resulting from the heat treatment of milk, which is known to have, in Fiore Sardo PDO, multiple consequences on cheese composition, microstructure and texture, which in turn affect appearance and sensory properties in general [6,10,52].
In general, MRI images confirmed, by a non-invasive approach, what already appeared at a first visual inspection of samples. RC cheeses exhibited a coarser paste, often with multiple cracks and holes, while HTC samples exhibited a smoother and more homogenous paste (Figure S5).
Two kinds of scientific considerations can be made based on the MRI images acquired. On one hand, as discussed above, MRI images give relaxometric information that is able to provide useful insights on the state of biopolymers and on the mobility of the system at molecular level. On the other hand, important information can be derived from the texture of MRI images by performing a detailed image analysis.
In recent years, image analysis techniques rapidly emerged not only at the academic level but also as rapid, economic, consistent and objective inspection techniques in the agricultural and food industry [54,55]. In particular, for the dairy industry sector, modern computer vision and image analysis methods have been used to evaluate cheese external or internal quality attributes such as colour, shape, texture, amount and distribution of added ingredients or production defects [56,57,58].
In this study, we pave the way for a deep-learning-based method for the classification of Fiore Sardo cheese images. The method consists of using data collected both from MRI and from a computer vision system based on photographic images taken with a smartphone camera. Deep learning is a subfield of artificial intelligence (AI) concerned with specific algorithms for thinking simulations and data processing inspired by the neurons of the brain [31]. Deep-learning-based image classification has proven to be highly effective in image classifying, object detection and other computer vision related problems [59,60].
All pre-trained DNN models run well in the open web GPU_Google colab service, with compatible and relatively short training times. For both MRI and photo datasets, the best classification result was achieved by the SqueezeNet model (Tables S3 and S4). In particular, for the T2 weighted MRI images, the SqueezeNet model exhibits a training accuracy from 93% to 98%, with the highest validation score (100%) (Table S3). As far as the photo dataset is concerned, SqueezeNet shows a training accuracy of 96%, with the same validation of 100% (Table S4).
The other pre-trained models (ResNet, AlexNet, VggNet, Densenet and Inception-net) also exhibit satisfactory performance results, ranging from a minimum of 72% to a maximum of 92%, with related validation ranging from 77% to 100% (Tables S3 and S4).
We propose here an IA approach to classify photos and T2 weighted MRI images of Fiore Sardo cheese made with raw or thermised milk, by applying a deep neural network (DNN), which includes several chained layers of processing between the input (photos or MRI images) and the output (classification/label of the input image).
As with most state-of-the-art IA methods, deep learning models usually require fitting thousands or millions of input data; thus, their application to those studies where the sample size is limited appears prohibitive. To address this challenge, the deep transfer learning approach offers a more suitable alternative. Transfer learning is the method by which the model uses the knowledge gained during the training of a relatively large dataset in a different but related problem (i.e., image classification and recognition). During the transfer learning process, only the selected classifier(s) is trained in the new network, while the features learned from the large dataset are transferred. Therefore, with the deep transfer learning approach, we do not need to retrain the entire network for a new dataset, thus allowing the training of deep learning models with relatively few data, reduced computational power and less training time [61].
With the main disadvantage of the present study being the relatively low number of available cheese samples, compared to the large data size usually needed for training deep learning frameworks, we chose to transfer the learning technique by testing six pre-trained DNN frameworks (ResNet, AlexNet, VggNet, SqueezeNet, DenseNet and Inception-Net), a widely used approach for addressing classification and image recognition challenges of MRI data.
In the end, all tested pre-trained DNNs yielded promising performances, giving an average classification accuracy of over 72% and over 93% in particular for the best performing SqueezeNet model (Tables S3 and S4). We did not perform further accuracy tests of the DNN pre-trained model in this study; however, we are proposing to implement and optimize this model in a dedicated study when more data are available.

4. Conclusions

Our data demonstrate that NMR relaxometry is able to assess molecular mobility changes induced in cheese by the most common industrial thermal treatments carried out on milk. In the present work, the experimental plan was optimized to exclude possible confounding factors. The same milk was used for producing both raw-milk or thermised-milk cheese. Moreover, analyses were repeated in different seasons, and the products from different producers were compared. Results were compared with a selection of industrial Fiore Sardo cheeses purchased in the local market.
NMR relaxometry data of Fiore Sardo cheese indicated the presence of two water populations, described by characteristic T2 relaxation times and area fractions. In particular, the first and fastest relaxing population showed a good predictive value. In order to obtain a relaxometric descriptor of FS cheese microstructure and molecular mobility, a new quality-related parameter was introduced, namely the Score factor (Ṩ). The Ṩ describes cheese in relation to the hydration status of the casein network and the state of the water-fat-protein network in cheese. We were able to discriminate FS cheeses produced from raw (RC) and heat-treated milk (HTC), with the latter exhibiting a lower Ṩ, indicating an enhanced hydration of the casein network. Industrially manufactured FS cheeses had generally lower Ṩ values than RC and HTC. Samples from the maturer industry, which buys cheeses before ripening and ages cheeses in industrial-like refrigerated cellars, showed similar values to artisanal cheeses, suggesting that NMR relaxometry is more sensitive to changes occurring in the cheesemaking steps preceding ripening. A further discriminative approach was finally carried out by processing acquired MRI images and photos with an IA method based on pre-trained deep learning algorithms. Preliminary classification results suggested that this innovative approach is promising in providing suitable discrimination between RC and HTC Fiore Sardo cheese.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/dairy2020023/s1, Figure S1: sample production and sampling scheme for MRI analysis, Figure S2: MRI sample preparation scheme, Figure S3: Computer Vision System (CVS) for photographic analysis, Figure S4: Image Dataset Directory and File Structure used, Figure S5: Representative T2 weighted MRI images of Fiore Sardo cheeses, Table S1: Production of FS cheeses and datasets, Table S2: Moisture content data of cheeses, Table S3: Performance and classification results of pre-trained DL frameworks on MRI images, Table S4: Performance and classification results of pre-trained DL frameworks on photographic images.

Author Contributions

Conceptualization, R.A.; Methodology, R.A., E.C. and R.M.; Software, E.C. and R.M; Formal Analysis, R.A., E.C. and R.M; Investigation, R.A., E.C. and R.M.; Data Curation, E.C. and R.M.; Writing—Original Draft Preparation, R.A., E.C. and R.M.; Writing—Review and Editing, R.A., E.C. and R.M.; Visualization, E.C. and R.M.; Supervision, R.A.; Funding Acquisition, R.A.; Project Scientific Coordination, R.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Sardinia Regional Government by means of Sardegna Ricerche (art 9, L.R. n. 15, 17.11.2010).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. European Commission. Quality Schemes Explained. Available online: https://ec.europa.eu/info/food-farming-fisheries/food-safety-and-quality/certification/quality-labels/quality-schemes-explained (accessed on 16 February 2021).
  2. European Commission eAmbrosia—The EU Geographical Indications Register. Available online: https://ec.europa.eu/info/food-farming-fisheries/food-safety-and-quality/certification/quality-labels/geographical-indications-register/# (accessed on 11 April 2021).
  3. Sardo, P.; Barletta, M. European Designations between Identity Values and Market. 2019. Available online: https://www.slowfood.com/wp-content/uploads/2019/09/European-designations-between-identity-values-and-market_ENG.pdf (accessed on 26 February 2021).
  4. Fox, P.; McSweeney, P.; Cogan, T.; Guinee, T. Cheese—Chemistry, Physics and Microbiology, 3rd ed.; Guinee, P., Fox, P., McSweeney, T., Cogan, T., Eds.; Academic Press: Cambridge, MA, USA, 2004; ISBN 978-1-4613-6138-1. [Google Scholar]
  5. Fox, P.; Uniacke-Lowe, T.; McSweeney, P.; O’Mahony, J. (Eds.) Heat-Induced Changes in Milk. In Dairy Chemistry and Biochemistry; Springer International Publishing: Basel, Switzerland, 2015; pp. 345–375. ISBN 978-3-319-14891-5. [Google Scholar]
  6. Pirisi, A.; Pinna, G.; Papoff, C.M. Effect of milk thermisation on Fiore Sardo PDO cheese: 1. Physicochemical characteristics. Sci. Tecn. Latt. Cas. 1999, 50, 353–366. [Google Scholar]
  7. EUR-Lex—31996R1236—EN—EUR-Lex. Available online: https://eur-lex.europa.eu/legal-content/IT/TXT/?uri=CELEX%3A31996R1236 (accessed on 26 February 2021).
  8. Anedda, R. Magnetic Resonance Analysis of Dairy Processing Suitable Tools for the Dairy Industry. In Magnetic Resonance in Food Science: Defining Food by Magnetic Resonance; The Royal Society of Chemistry: London, UK, 2015; pp. 49–64. ISBN 978-1-78262-031-0. [Google Scholar]
  9. Anedda, R.; Pala, M.; Curti, E. Cheeses That Made History in Italian Dairy Tradition. In Cheeses around the World: Types, Production, Properties and Cultural and Nutritional Relevance; Ferrão, A.C., Guiné, R.P.F., Correia, P.R., Eds.; Nova Science Publishers: New York City, NY, USA, 2019; pp. 203–228. ISBN 978-1-53615-419-1. [Google Scholar]
  10. Mendia, C.; Ibañez, F.C.; Torre, P.; Barcina, Y. Effect of pasteurization on the sensory characteristics of a ewe’s-milk cheese. J. Sens. Stud. 1999, 14, 415–424. [Google Scholar] [CrossRef]
  11. Scintu, M.F.; Del Caro, A.; Urgeghe, P.P.; Piga, C.; Di Salvo, R. Sensory profile development for an Italian PDO ewe’s milk cheese at two different ripening times. J. Sens. Stud. 2010, 25, 577–590. [Google Scholar] [CrossRef]
  12. Siniscalchi, V. Environment, regulation and the moral economy of food in the Slow Food movement. J. Polit. Ecol. 2013, 20, 295–305. [Google Scholar] [CrossRef] [Green Version]
  13. Pittalis, P.; Cappellacci, U.; Contu, M.I.; Fasolino, G.; Peru, A.; Tedde, M.; Tocco, E.; Tunis, S.; Zedda, A. XV Legislature—Regional Question N.1380. Available online: http://www3.consregsardegna.it/XVlegislatura/interrogazioni/1380 (accessed on 26 February 2021).
  14. Ledda A XV Legislature—Regional Question N.1706/A. Available online: http://www3.consregsardegna.it/XVlegislatura/interrogazioni/1706 (accessed on 26 February 2021).
  15. Siniscalchi, V. Food, Slow Food and Middle Class Activism. In Food Activism: Agency, Democracy and Economy; Bloomsbury Academic: London, UK, 2014; pp. 9–24. [Google Scholar]
  16. Anedda, R.; Pardu, A.; Korb, J.P.; Curti, E. Effect of the manufacturing process on Fiore Sardo PDO cheese microstructure by multi-frequency NMR relaxometry. Food Res. Int. 2021, 140, 110079. [Google Scholar] [CrossRef]
  17. Caboni, P.; Maxia, D.; Scano, P.; Addis, M.; Dedola, A.; Pes, M.; Murgia, A.; Casula, M.; Profumo, A.; Pirisi, A. A gas chromatography-mass spectrometry untargeted metabolomics approach to discriminate Fiore Sardo cheese produced from raw or thermized ovine milk. J. Dairy Sci. 2019, 102, 5005–5018. [Google Scholar] [CrossRef]
  18. Dedola, A.S.; Piras, L.; Addis, M.; Pirisi, A.; Piredda, G.; Mara, A.; Sanna, G. New analytical tools for unmasking frauds in raw milk-based dairy products: Assessment, validation and application to fiore sardo PDO cheese of a RP-HPLC method for the evaluation of the α-l-fucosidase activity. Separations 2020, 7, 40. [Google Scholar] [CrossRef]
  19. Mulas, G.; Roggio, T.; Uzzau, S.; Anedda, R. A new magnetic resonance imaging approach for discriminating Sardinian sheep milk cheese made from heat-treated or raw milk. J. Dairy Sci. 2013, 96, 7393–7403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Mazza, M.; Guglielmetti, C.; Brusadore, S.; Sciuto, S.; Esposito, G.; Caramelli, M.; Peletto, S.; Acutis, P.L.; Marengo, E.; Manfredi, M.; et al. A proteomic approach to the safeguard of a typical agri-food product: Fiore sardo PDO. Adv. Dairy Res. 2019, 7. [Google Scholar] [CrossRef]
  21. Piga, C.; Urgeghe, P.P.; Piredda, G.; Scintu, M.F.; Di Salvo, R.; Sanna, G. Thermal inactivation and variability of γ-glutamyltransferase and α-l-fucosidase enzymatic activity in sheep milk. LWT Food Sci. Technol. 2013, 54, 152–156. [Google Scholar] [CrossRef]
  22. Lambelet, P.; Berrocal, R.; Ducret, F. Low resolution NMR spectroscopy: A tool to study protein denaturation: I. Application to diamagnetic whey proteins. J. Dairy Res. 1989, 56, 211–222. [Google Scholar] [CrossRef] [Green Version]
  23. Lambelet, P.; Berrocal, R.; Renevey, F. Low-field nuclear magnetic resonance relaxation study of thermal effects on milk proteins. J. Dairy Res. 1992, 59, 517–526. [Google Scholar] [CrossRef]
  24. Curti, E.; Pardu, A.; Melis, R.; Addis, M.; Pes, M.; Pirisi, A.; Anedda, R. Molecular mobility changes after high-temperature, short-time pasteurization: An extended time-domain nuclear magnetic resonance screening of ewe milk. J. Dairy Sci. 2020, 103, 9881–9892. [Google Scholar] [CrossRef]
  25. Curti, E.; Pardu, A.; Del Vigo, S.; Sanna, R.; Anedda, R. Non-invasive monitoring of curd syneresis upon renneting of raw and heat-treated cow’s and goat’s milk. Int. Dairy J. 2019, 90. [Google Scholar] [CrossRef]
  26. Curti, E.; Pardu, A.; Del Vigo, S.; Sanna, R.; Anedda, R. A low-field Nuclear Magnetic Resonance dataset of whole milk during coagulation and syneresis. Data Br. 2019, 26. [Google Scholar] [CrossRef] [PubMed]
  27. Bjarnason, T.A.; Mitchell, J.R. AnalyzeNNLS: Magnetic resonance multiexponential decay image analysis. J. Magn. Reson. 2010, 206, 200–204. [Google Scholar] [CrossRef] [PubMed]
  28. AOAC International. AOAC 948.12-2002, Moisture in Cheese. Method II (Rapid Screening Method): AOAC Official Method. Available online: http://www.aoacofficialmethod.org/index.php?main_page=product_info&products_id=1145 (accessed on 26 February 2021).
  29. Pang, Z.; Chong, J.; Li, S.; Xia, J. MetaboAnalystR 3.0: Toward an optimized workflow for global metabolomics. Metabolites 2020, 10, 186. [Google Scholar] [CrossRef]
  30. Laurence, A. NIfTI Image Converter (nii2png) for Python and Matlab|NIfTI-Image-Converter. Available online: https://alexlaurence.github.io/NIfTI-Image-Converter/ (accessed on 26 February 2021).
  31. Shrestha, A.; Mahmood, A. Review of deep learning algorithms and architectures. IEEE Access 2019, 7, 53040–53065. [Google Scholar] [CrossRef]
  32. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep residual learning for image recognition. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2016; pp. 770–778. [Google Scholar]
  33. Krizhevsky, A.; Sutskever, I.; Hinton, G.E. ImageNet classification with deep convolutional neural networks. In Proceedings of the 25th International Conference on Neural Information Processing Systems, Siem Reap, Cambodia, 13–16 December 2018; Curran Associates Inc.: Red Hook, NY, USA, 2012; Volume 1, pp. 1097–1105. [Google Scholar]
  34. Liu, S.; Deng, W. Very deep convolutional neural network based image classification using small training sample size. In Proceedings of the 2015 3rd IAPR Asian Conference on Pattern Recognition (ACPR), Kuala Lumpur, Malaysia, 3–6 November 2015; pp. 730–734. [Google Scholar]
  35. Iandola, F.N.; Moskewicz, M.W.; Ashraf, K.; Han, S.; Dally, W.J.; Keutzer, K. SqueezeNet: AlexNet-level accuracy with 50x fewer parameters and <1MB model size. arXiv 2016, arXiv:1602.07360. [Google Scholar]
  36. Huang, G.; Liu, Z.; Van Der Maaten, L.; Weinberger, K.Q. Densely connected convolutional networks. In Proceedings of the 2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Honolulu, HI, USA, 21–26 July 2017; pp. 2261–2269. [Google Scholar]
  37. Szegedy, C.; Vanhoucke, V.; Ioffe, S.; Shlens, J.; Wojna, Z. Rethinking the inception architecture for computer vision. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2016; pp. 2818–2826. [Google Scholar]
  38. Google Google Colaboratory. Available online: https://colab.research.google.com/notebooks/intro.ipynb (accessed on 26 February 2021).
  39. Hills, B.P.; Takacs, S.F.; Belton, P.S. A new interpretation of proton NMR relaxation time measurements of water in food. Food Chem. 1990, 37, 95–111. [Google Scholar] [CrossRef]
  40. Mariette, F.; Tellier, C.; Brule, G.; Marchal, P. Multinuclear Nmr-Study of the Ph dependent water state in skim milk and caseinate solutions. J. Dairy Res. 1993, 60, 175–188. [Google Scholar] [CrossRef]
  41. Mariette, F. NMR Relaxation of Dairy Products. In Modern Magnetic Resonance; Webb, G.A., Ed.; Springer: Dordrecht, The Netherlands; Berlin, Germany, 2006; pp. 1697–1701. ISBN 978-1-4020-3894-5. [Google Scholar]
  42. Gianferri, R.; D’Aiuto, V.; Curini, R.; Delfini, M.; Brosio, E. Proton NMR transverse relaxation measurements to study water dynamic states and age-related changes in Mozzarella di Bufala Campana cheese. Food Chem. 2007, 105, 720–726. [Google Scholar] [CrossRef]
  43. Gianferri, R.; Maioli, M.; Delfini, M.; Brosio, E. A low-resolution and high-resolution nuclear magnetic resonance integrated approach to investigate the physical structure and metabolic profile of Mozzarella di Bufala Campana cheese. Int. Dairy J. 2007, 17, 167–176. [Google Scholar] [CrossRef]
  44. Chaland, B.; Mariette, F.; Marchal, P.; De Certaines, J. 1H nuclear magnetic resonance relaxometric characterization of fat and water states in soft and hard cheese. J. Dairy Res. 2000, 64, 609–618. [Google Scholar] [CrossRef]
  45. Godefroy, S.; Korb, J.-P.; Creamer, L.K.; Watkinson, P.J.; Callaghan, P.T. Probing protein hydration and aging of food materials by the magnetic field dependence of proton spin-lattice relaxation times. J. Colloid Interface Sci. 2003, 267, 337–342. [Google Scholar] [CrossRef]
  46. Venu, K.; Denisov, V.P.; Halle, B. Water 1H magnetic relaxation dispersion in protein solutions. A quantitative assessment of internal hydration, proton exchange, and cross-relaxation. J. Am. Chem. Soc. 1997, 119, 3122–3134. [Google Scholar] [CrossRef]
  47. Conte, P.; Cinquanta, L.; Lo Meo, P.; Mazza, F.; Micalizzi, A.; Corona, O. Fast field cycling NMR relaxometry as a tool to monitor Parmigiano Reggiano cheese ripening. Food Res. Int. 2021, 139, 109845. [Google Scholar] [CrossRef]
  48. Boiani, M.; Fenelon, M.; FitzGerald, R.J.; Kelly, P.M. Use of 31P NMR and FTIR to investigate key milk mineral equilibria and their interactions with micellar casein during heat treatment. Int. Dairy J. 2018, 81, 12–18. [Google Scholar] [CrossRef] [Green Version]
  49. Wahlgren, N.M.; Dejmek, P.; Drakenberg, T. A 43Ca and 31P NMR study of the calcium and phosphate equilibria in heated milk solutions. J. Dairy Res. 1990, 57, 355–364. [Google Scholar] [CrossRef]
  50. Pisanu, S.; Pagnozzi, D.; Pes, M.; Pirisi, A.; Roggio, T.; Uzzau, S.; Addis, M.F. Differences in the peptide profile of raw and pasteurised ovine milk cheese and implications for its bioactive potential. Int. Dairy J. 2015, 42, 26–33. [Google Scholar] [CrossRef]
  51. Grappin, R.; Beuvier, E. Possible implications of milk pasteurization on the manufacture and sensory quality of ripened cheese. Int. Dairy J. 1997, 7, 751–761. [Google Scholar] [CrossRef]
  52. Singh, H.; Waungana, A. Influence of heat treatment of milk on cheesemaking properties. Int. Dairy J. 2001, 11, 543–551. [Google Scholar] [CrossRef]
  53. Slade, L.; Levine, H. Beyond water activity: Recent advances based on an alternative approach to the assessment of food quality and safety. Crit. Rev. Food Sci. Nutr. 1991, 30, 115–360. [Google Scholar] [CrossRef]
  54. Brosnan, T.; Sun, D.-W. Inspection and grading of agricultural and food products by computer vision systems—A review. Comput. Electron. Agric. 2002, 36, 193–213. [Google Scholar] [CrossRef]
  55. Narendra, V.G.; Hareesha, K.S. Quality inspection and grading of agricultural and food products by computer vision—A review. Int. J. Comput. Appl. 2010, 2. [Google Scholar] [CrossRef]
  56. Wang, H.-H.; Sun, D.-W. Melting characteristics of cheese: Analysis of effect of cheese dimensions using computer vision techniques. J. Food Eng. 2002, 52, 279–284. [Google Scholar] [CrossRef]
  57. Caccamo, M.; Melilli, C.; Barbano, D.M.; Portelli, G.; Marino, G.; Licitra, G. Measurement of gas holes and mechanical openness in cheese by image analysis. J. Dairy Sci. 2004, 87, 739–748. [Google Scholar] [CrossRef] [Green Version]
  58. Ni, H.; Guansekaran, S. Image processing algorithm for cheese shred evaluation. J. Food Eng. 2004, 61, 37–45. [Google Scholar] [CrossRef]
  59. Zhou, L.; Zhang, C.; Liu, F.; Qiu, Z.; He, Y. Application of deep learning in food: A review. Compr. Rev. Food Sci. Food Saf. 2019, 18, 1793–1811. [Google Scholar] [CrossRef] [Green Version]
  60. Lundervold, A.S.; Lundervold, A. An overview of deep learning in medical imaging focusing on MRI. Z. Med. Phys. 2019, 29, 102–127. [Google Scholar] [CrossRef]
  61. Weiss, K.; Khoshgoftaar, T.M.; Wang, D. A survey of transfer learning. J. Big Data 2016, 3, 9. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Representative T2 distribution profile of an aged FS cheese at 300 MHz. Two proton populations are observed, characterized by geometric mean T2 values (T21 and T22) and area fractions of each population (AF1 and AF2).
Figure 1. Representative T2 distribution profile of an aged FS cheese at 300 MHz. Two proton populations are observed, characterized by geometric mean T2 values (T21 and T22) and area fractions of each population (AF1 and AF2).
Dairy 02 00023 g001
Figure 2. Score factors (Ṩ) of Fiore Sardo cheeses from Dataset 1: 105 days (a) and 180 days (b) of ripening, Season 1 (March–April 2019); 105 days (c) and 180 days (d) of ripening, Season 2 (January–February 2020). Blue circles represent samples made from raw milk (RC) and red circles represent cheeses made starting from heat-treated milk (HTC).
Figure 2. Score factors (Ṩ) of Fiore Sardo cheeses from Dataset 1: 105 days (a) and 180 days (b) of ripening, Season 1 (March–April 2019); 105 days (c) and 180 days (d) of ripening, Season 2 (January–February 2020). Blue circles represent samples made from raw milk (RC) and red circles represent cheeses made starting from heat-treated milk (HTC).
Dairy 02 00023 g002
Figure 3. Box plots of the Score values (Ṩ) of each group of cheeses (RC, HTC, industrial and maturer cheeses).
Figure 3. Box plots of the Score values (Ṩ) of each group of cheeses (RC, HTC, industrial and maturer cheeses).
Dairy 02 00023 g003
Figure 4. Representative images of the visual appearance of RC, HTC and industrial samples: Producer S1, Season 2, 105 days of ripening, Fiore Sardo from raw milk (upper row) and from heat-treated milk (lower row) (a); details of RC (b) and HTC (c) samples ripened 105 days upon cutting; industrial Fiore Sardo purchased from a local market, with at least 105 days of ripening (d).
Figure 4. Representative images of the visual appearance of RC, HTC and industrial samples: Producer S1, Season 2, 105 days of ripening, Fiore Sardo from raw milk (upper row) and from heat-treated milk (lower row) (a); details of RC (b) and HTC (c) samples ripened 105 days upon cutting; industrial Fiore Sardo purchased from a local market, with at least 105 days of ripening (d).
Dairy 02 00023 g004
Table 1. T21 and T22 values (ms) and corresponding mean area fractions AF1 and AF2 (%) of artisanal FS cheeses produced in Season 1 and Season 2; asterisks (*) in each row indicate significant differences (t-test, p < 0.05) between RC and HTC parameters (T21 and T22, AF1 and AF2) for each sample (M = Mean; SD = Standard Deviation).
Table 1. T21 and T22 values (ms) and corresponding mean area fractions AF1 and AF2 (%) of artisanal FS cheeses produced in Season 1 and Season 2; asterisks (*) in each row indicate significant differences (t-test, p < 0.05) between RC and HTC parameters (T21 and T22, AF1 and AF2) for each sample (M = Mean; SD = Standard Deviation).
T21 T22 AF1 AF2
RCHTC RCHTC RCHTC RCHTC
MSDMSD MSDMSD MSDMSD MSDMSD
Season 1105 daysS111.160.479.480.40*52.583.7056.011.56*0.750.070.850.02*0.240.070.140.02*
S211.280.949.240.42*49.053.8452.151.85*0.700.030.780.03*0.320.050.210.03*
S310.740.789.130.58*52.373.5447.510.77*0.820.020.820.02 0.170.020.180.02
S48.460.318.640.51 48.334.1552.673.12*0.710.050.770.03*0.280.050.210.03*
S58.860.319.330.52*51.731.4849.691.46*0.770.020.810.02*0.230.020.190.02*
180 daysS111.941.178.760.24*47.833.9148.312.23 0.520.090.620.03*0.440.060.370.03*
S28.830.148.180.23*44.921.3853.341.05*0.740.020.770.02*0.250.020.220.02*
S310.520.917.790.47*51.981.9244.311.81*0.790.010.790.03 0.200.010.200.02
S49.741.647.980.28*40.673.2947.272.17*0.560.030.700.03*0.420.080.290.03*
S58.500.318.800.74*49.421.4946.661.93*0.670.030.690.05*0.320.030.300.05*
Season 2105 daysS112.701.6012.700.69 52.472.8557.722.80*0.680.050.780.04*0.310.050.220.04*
S214.111.4512.230.93*54.401.5953.042.98*0.630.060.780.04*0.360.060.210.04*
S310.260.2610.020.39*45.931.0552.221.46*0.790.010.840.03*0.200.020.160.03*
S410.201.4810.490.31 46.953.2850.162.18*0.670.070.790.01*0.320.070.200.01*
S59.291.0810.740.42*43.553.0853.181.05*0.680.020.820.01*0.310.020.170.01*
180 daysS19.080.188.820.50*49.262.6752.374.38*0.620.040.780.06*0.370.040.210.06*
S29.940.669.210.44*53.102.3552.331.93 0.720.050.800.03*0.270.050.190.03*
S310.390.508.780.35*50.971.0152.901.61*0.720.030.760.03*0.270.030.230.03*
S49.090.748.310.51*50.873.8748.033.28*0.680.050.730.03*0.310.050.260.03*
S514.552.2410.070.38*44.943.3453.001.14*0.570.030.720.02*0.420.030.270.02*
Table 2. T21 and T22 values and correspondent mean area fractions AF1 and AF2 of industrial FS cheeses (M = Mean; SD = Standard Deviation).
Table 2. T21 and T22 values and correspondent mean area fractions AF1 and AF2 of industrial FS cheeses (M = Mean; SD = Standard Deviation).
T21T22Area 1Area 2
MSDMSDMSDMSD
I18.3010.05540.3370.1000.8350.0020.1570.002
I210.0320.05652.0460.1410.8910.0010.1020.001
I39.2140.03845.0780.0350.8140.0010.1780.001
I48.5230.03029.6460.0430.8710.0020.1200.002
I58.0180.02451.9360.0470.6970.0020.2940.002
I68.7420.02044.3190.0420.8360.0010.1560.001
I79.3180.02351.3520.1330.8580.0120.1300.000
I89.4620.02251.0540.0580.8130.0010.1800.000
I99.4630.01851.0580.0370.8130.0010.1800.001
I109.4660.11448.4700.2920.8760.0020.1170.002
I119.1880.04746.3140.1310.8430.0010.1500.001
I127.7020.09834.8870.0840.5750.0040.4150.004
I137.3010.02934.7720.0130.5450.0020.4420.002
Table 3. Statistical comparison of Score (Ṩ) values (ANOVA, p < 0.05) of the four groups of Fiore Sardo samples, namely RC, HTC (produced in Season 1 and Season 2), maturer and industrial cheeses. Each shepherd production (RC and HTC) was individually compared to maturer and industrial cheeses; data are expressed as Mean ± Standard Deviation; means sharing the same superscript letters are not significantly different from each other; the letter “a” is assigned to the higher mean.
Table 3. Statistical comparison of Score (Ṩ) values (ANOVA, p < 0.05) of the four groups of Fiore Sardo samples, namely RC, HTC (produced in Season 1 and Season 2), maturer and industrial cheeses. Each shepherd production (RC and HTC) was individually compared to maturer and industrial cheeses; data are expressed as Mean ± Standard Deviation; means sharing the same superscript letters are not significantly different from each other; the letter “a” is assigned to the higher mean.
RCHTCMaturerIndustrial
S1Season 1105 days14.99 ± 1.70 a11.17 ± 0.62 c13.39 ± 0.18 b10.92 ± 0.62 c
180 days23.83 ± 4.83 a14.22 ± 1.22 b13.39 ± 0.18 b10.92 ± 0.62 c
Season 2105 days20.25 ± 2.38 a16.43 ± 1.60 b13.39 ± 0.18 c10.92 ± 0.62 d
180 days14.78 ± 1.00 a11.38 ± 1.45 c13.39 ± 0.18 b10.92 ± 0.62 c
S2Season 1105 days16.18 ± 1.85 a11.87 ± 0.90 c13.39 ± 0.18 b10.92 ± 0.62 d
180 days11.90 ± 0.81 b10.57 ± 0.48 c13.39 ± 0.18 a10.92 ± 0.62 c
Season 2105 days24.10 ± 3.44 a15.74 ± 1.59 b13.39 ± 0.18 c10.92 ± 0.62 d
180 days13.86 ± 1.74 a11.58 ± 0.86 b13.39 ± 0.18 a10.92 ± 0.62 c
S3Season 1105 days13.12 ± 1.19 a11.19 ± 0.90 b13.39 ± 0.18 a10.92 ± 0.62 b
180 days13.33 ± 1.29 a9.86 ± 0.73 c13.39 ± 0.18 a10.92 ± 0.62 b
Season 2105 days13.02 ± 0.51 a11.98 ± 0.82 b13.39 ± 0.18 a10.92 ± 0.62 c
180 days14.47 ± 0.94 a11.53 ± 0.70 c13.39 ± 0.18 b10.92 ± 0.62 d
S4Season 1105 days12.00 ± 0.97 b11.18 ± 0.75 c13.39 ± 0.18 a10.92 ± 0.62 c
180 days17.74 ± 4.27 a11.48 ± 0.75 bc13.39 ± 0.18 b10.92 ± 0.62 c
Season 2105 days15.65 ± 3.68 a13.23 ± 0.40 b13.39 ± 0.18 b10.92 ± 0.62 c
180 days13.41 ± 1.77 a11.42 ± 0.96 b13.39 ± 0.18 a10.92 ± 0.62 b
S5Season 1105 days11.58 ± 0.56 b11.57 ± 0.81 b13.39 ± 0.18 a10.92 ± 0.62 c
180 days12.74 ± 0.81 a12.76 ± 1.39 a13.39 ± 0.18 a10.92 ± 0.62 b
Season 2105 days12.93 ± 1.10 a13.08 ± 0.53 a13.39 ± 0.18 a10.92 ± 0.62 b
180 days25.61 ± 2.79 a14.05 ± 0.69 b13.39 ± 0.18 b10.92 ± 0.62 c
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Anedda, R.; Melis, R.; Curti, E. Quality Control in Fiore Sardo PDO Cheese: Detection of Heat Treatment Application and Production Chain by MRI Relaxometry and Image Analysis. Dairy 2021, 2, 270-287. https://0-doi-org.brum.beds.ac.uk/10.3390/dairy2020023

AMA Style

Anedda R, Melis R, Curti E. Quality Control in Fiore Sardo PDO Cheese: Detection of Heat Treatment Application and Production Chain by MRI Relaxometry and Image Analysis. Dairy. 2021; 2(2):270-287. https://0-doi-org.brum.beds.ac.uk/10.3390/dairy2020023

Chicago/Turabian Style

Anedda, Roberto, Riccardo Melis, and Elena Curti. 2021. "Quality Control in Fiore Sardo PDO Cheese: Detection of Heat Treatment Application and Production Chain by MRI Relaxometry and Image Analysis" Dairy 2, no. 2: 270-287. https://0-doi-org.brum.beds.ac.uk/10.3390/dairy2020023

Article Metrics

Back to TopTop