Long-Term Ripening Evaluation of Ewes ’ Cheeses by Fourier-Transformed Infrared Spectroscopy under Real Industrial Conditions

The ripening changes over time of special cheeses (Pecorino, ewes’ ripe, and Gouda) made with ewes’ milk were evaluated using FTIR/ATR spectroscopy during approximately one year. The midinfrared FTIR/ATR analyses were carried out in different ripening times between the cheese varieties and processed by means of multivariate statistical approaches. Overall, during the maturation, we observed a downward trend of the absorbance intensity of the amide group peaks (1700 to 1500 cm), which is linked to the breakdown of peptide bonds. Similar behavior was obtained for the lipidic region (3000 to 2800 cm and 1765 to 1730 cm). Hierarchical cluster analysis and principal component analysis allowed the evaluation of the physicochemical changes of the cheeses. The proteolysis occurs in a fast pace during the first trimester of the ripening process, and the lipids are converted to smaller species as the times goes by. Our results indicate that infrared spectroscopy can be a useful tool in determining optimal temporal parameters in stages involving the development, production, and even a possible estimation of shelf life of cheeses.


Introduction
Cheese is the most abundant processed dairy product in the world.The global production is around 23 million tons, representing an increase of 791,000 tons between 2013 and 2014 [1].The United States, for example, was responsible for approximately a quarter of the world production in 2014 and expected to reach 286,000 tons in export trades by the end of 2017 [1,2].Thus, the quality control of cheeses is of great importance, especially when a long period of ripening is required before marketing the product.About 70% of the total milk produced in Brazil is processed into cheeses, where approximately 1100 tons are consumed within the domestic market per year.The special cheese market is growing throughout the years, as well as the research and development of products adapted for the Brazilian market [3].
The physical-chemical and sensorial characteristics of dairy products can be changed by several factors including herd genetics, health, feeding, and environmental conditions, among others [4].Cheese ripening involves several biochemical reactions, leading to the development of the unique sensorial attributes for each cheese variety [5][6][7][8][9][10].The formation of smaller peptides and free amino acids has a direct impact on flavor and texture, due to the breakdown of the protein network [7].In addition to that, lipolysis generates free fatty acids and organic acids that are further metabolized to volatile and nonvolatile flavor components [11].
Usual methods for analyzing cheese composition and ripening index are time-consuming and expensive.There is a need for new methodologies and instrumentation to observe the main changes during the ripening process.In fact, nowadays, several noninvasive spectroscopic techniques have gained more and more space due to several reasons: (i) they are environmentally friendly due to the absence of reagents; (ii) they are fast due to the absence of sample preparation; and (iii) they are susceptible to miniaturization [12][13][14][15].In particular, instrumental technique such as midinfrared (MIR) spectroscopy has been widely used in the analysis of food products [16].Encompassing the spectral range of 4000 to 400 cm −1 , MIR spectroscopy allows access to the vibrational modes of the samples which may result in bands that are assigned to specific chemical groups.In 2005, AOAC international has approved the Fourier transform infrared (FTIR) method for determination of fat, protein, lactose, and total solid contents in some dairy products [17].Since then, studies on cheese ripening using spectroscopic methods have been published [18][19][20][21][22][23].However, they were all carried out during a relatively short period of time, that is, less than 3 months.A better exploration of the maturation time requires a longer period of analysis, once the development of sensory attributes such as texture and flavor is directly related to the rate of breakdown of the casein matrix over time, especially in hard and semihard-type cheeses [11].
This paper deals with three ewes' cheeses under real industrial conditions and named "Pecorino," "Gouda," and "ewes' ripe."Although present in the European and North American markets, the "Pecorino" and "Sheep Gouda," the information about ewes' cheeses are in general scarce [24].Concerning Brazil, the three cheeses are promising novelties in the market, especially the "ewes' ripe," as it is a variation of "Minas" cheese (made originally with cow's milk).To evaluate them, attenuated total reflection (ATR) unit was coupled to FTIR spectrometer to follow the changes in the cheeses over a long period (approximately one year).As a result, we show that infrared spectroscopy provides secure hints to determine optimal temporal parameters in ripening stages involving the development of the aforementioned cheeses.

Cheese Samples.
A dairy industry with Brazilian federal inspection seal (SIF), located in the municipality of Soledade de Minas, Minas Gerais, Brazil, supplied three different varieties of cheese (Pecorino, Gouda, and ewes' ripe) made with ewes' milk.The cheeses had different production dates and at aging times (1 to 12 months), so that represents real industrial process conditions.

Cheese Composition.
The compositions of the supplied cheeses were carried out in terms of moisture, fat, and protein content to classify them according to the Brazilian legislation.The measurements were carried out following Lutz [25] and Brasil [26], which establish the physical and chemical official analytical methods for control of milk and dairy products.All samples were tested in duplicate.
2.3.FTIR/ATR Spectroscopy Measurements.Infrared spectroscopy was performed using a FT-MIR spectrophotometer (Bruker, model Vertex 70), with ATR accessory and OPUS 6.5 software.Infrared spectra were recorded in absorbance mode, between 4000 and 400 cm −1 , with a resolution of 4 cm −1 and 32 scans accumulated by sample.The data were analyzed in the software OriginPro 8.0, Excel, and Minitab.
Cheese slices were prepared with approximately 4 cm 2 of area and thickness of 1 cm.The samples were taken from the center of the cheeses to avoid any contamination during the storage time.The measurements were performed in quadruplicate, using two pieces of each cheese.
2.4.Statistical Analysis.The differences in group mean values were submitted to analysis of variance (one-way ANOVA).Tukey HSD test was carried out when significant differences among groups were obtained.The statistical tests were performed with a significance level equals to 5% using Statistica 7.1 software.The spectra obtained between 3000 and 2800 cm −1 , 1766 and 1730 cm −1 , and 1700 and 1000 cm −1 were submitted to multivariate analysis using Minitab 16.1.1software.Hierarchical cluster analysis (HCA) was performed with different linkage methods (data not shown), where the clustering using the McQuitty linkage provided the best overview of the results.Finally, the principal component analysis (PCA) was applied to the spectral sets to evaluate the evolution during ripening.

Results and Discussion
3.1.Physical-Chemical Characterization of the Cheeses.The earliest stage of maturation of each cheese received, that is Pecorino (2nd month), Gouda (2nd month), and ewes' ripe (1st month), was submitted to physical-chemical analysis to estimate values of moisture, total protein, and fat in dry matter (FDM) as presented in Table 1.All cheese samples were classified as "fatty" since their FDM values were between 45.0% and 59.9%.In relation to moisture, the Pecorino cheese revealed being the variety with the lowest value, "hard," while ewes' ripe and Gouda were classified as "average moisture" or "semihard" [27].Moisture values were within the allowed range according to the local regulation [26].Analysis of variance revealed that the cheeses presented significant difference between them, regarding the above parameters.2 Journal of Spectroscopy cheeses over some months (M).In general, for all samples, the spectra presented bands around the same spectral positions which, nevertheless, differ in intensity a priori as well as along the time.The difference in intensity per se is related

FTIR/ATR-MIR Spectroscopy
to the nature of the cheese and is weighted to different functional groups [20,28,29].The differences along the time have their origins linked to the biochemical processes occurring during ripening.The biochemistry of the cheese ripening is usually evaluated under three general subjects: proteolysis, lipolysis, and glycolysis [11,30].Among these processes, the first two are considered the most important for the great majority of the matured cheeses [5][6][7][8][9][10].Therefore, the spectra were subdivided in narrow regions, to ease the interpretation of the data over time.These regions were conveniently handled by multivariate statistical analysis.
3.2.1.Amide I and II Region.The area between 1700 and 1500 cm −1 are related to peptide bonds (CO-NH), representing two major peaks: amide I in 1640 cm −1 (νC=O, νC-N) and amide II in 1550 cm −1 (δN-H, νC-N) [29,31].The behavior of these peaks during the ripening processes of the three types of cheese can be observed in Figure 2.
In general, all varieties studied have shown a decrease in the absorbance in both bands during the development of the ripening process.This indicates the occurrence of the hydrolysis of proteins by breakdown of the peptide bonds.This conclusion is supported by Fox et al. [6].
Martín-del-Campo et al. [29] have used midinfrared spectroscopy to study the ripening phenomena of two cheese zones (core and under-rind) of Camembert cheeses.Significant changes were recorded for amide I and II bands in the under-rind and only for amide II in the core.The maturation time analyzed was of 27 days where the absorbance values decreased until the tenth day and then tended to increase up to the end of the maturation.Mazerolles et al. [32] studied sixteen experimental semihard cheeses in the range of amide I and II at 1, 21, 51, and 81 days of maturation.Their results have shown a continuous decrease in the amide I band and an increase of amide II.
In Figure 2, one can observe that even though the absorbance intensities have decreased systematically over-ripening time, they do not occur in the same way among the cheese varieties analyzed and not with the same rate in the months of ripening.Regarding the Pecorino cheese (Figure 2(a)), the highest variation in absorbance levels (linked to proteolysis) appeared in the months between 4 and 5 of the maturation process.In ewes' ripe cheeses (Figure 2(c)), this fact could be verified in the months between 1 and 3, while in Gouda cheeses (Figure 2(b)), the proteolysis seemed to happen in a more constant pace, meaning that the variations among the months were similar up to the 5th month.The difference in the intensity decay of the peaks could be explained by the fact that the casein hydrolysis (proteolysis), lipid hydrolysis (lipolysis), and lactose fermentation processes are dependent on the milk composition, the fabrication process, pH, salt content, water activity, storage temperature, and specially, the type of yeast present in the cheese [8,11,30].
Studies evaluating the amide peaks in cheeses have stated that changes in the intensity of the amide I peak might be associated with modification in the secondary structure of casein, protein aggregation, and protein/water interaction      Journal of Spectroscopy [29,[32][33][34].Therefore, since the cheeses studied here had different maturation times and were produced with different batches of milk, the factors mentioned above may have influenced the proteolysis process in different ways.Farkye and Fox [35] stated that for the most of the semihard and hard cheeses, the proteolysis process is commonly used as a ripening index.Thus, through the behavior and evolution of the peaks obtained in the FTIR/ATR spectra, the technique can be extremely useful to verify the occurrence of the proteolysis phenomenon during the cheese ripening.

Lipid and Carbohydrate
Regions. Figure 3 shows three peaks associated with the lipid fraction.According to Dufour et al. [28], the determination of fat is commonly used by the association of the C = O (~1750 cm −1 ) group and acyl chain C-H (3000-2800 cm −1 ) stretch wavenumbers.In Figure 3, in general, the signal of methylene peaks (bands around 2922 and 2852 cm −1 ) and those associated with ester and organic acid carbonyl (νC=O in 1743 cm −1 ) have decreased as the cheeses matured.Similar behavior was found by Martín-del-campo et al. [29].However, the cheeses with ripening time of 5 months (Pecorino), 5 months (Gouda), and 3 months (ewes' ripe) have shown the highest absorbance bands.According to Ticiani et al. [36], at the highest lactation stage, the ewes lose the ability of keeping the fat content stable.A study developed by Brito et al. [37] showed that the ewes raised in Brazil have the peak of lactation occurring in spring (Sept-Nov), with increased fat content of the milk in approximately 2.5%.Once the above-mentioned cheeses were made during the springtime, the deviation observed in the absorbance peaks may rely in the fluctuations of fat content during the peak of the lactation.
Collins et al. [38] stated that lipolysis occurs in most or all cheeses, but its intensity is more extensive in some hard Italian varieties and in Blue cheese.In addition, regarding the ripening of Cheddar, Gouda, and Swiss cheeses, the contribution of lipolysis is relatively low.Thus, among the cheese varieties studied here, the Gouda made with ewes' milk was the one that has had the least influence of lipolysis and smaller decreases in absorption during its maturation process (Figure 3(b)).
The peaks located around 1377 cm −1 (δ s CH 3 ) are assigned to glucose and galactose.Surrounding 1161 cm −1 are located peaks related to the sum of lactose (tertiary alcohol νC-OH) and monosaccharide (endocyclic oxygen νC-O) signals (Figure 1) [29].The evolution of these peaks during the ripening process has shown a similar behavior as the lipid bands, that is, a decrease in the absorbance values.

Ripening Time Discrimination by Multivariate Analysis.
The different cheeses and the correspondent ripening times were submitted to the multivariate approaches HCA and PCA to observe similarity patterns and possible correlations amongst the samples.The wavenumber range used for multivariate analysis was chosen to include a wide variety of chemical groups.The spectra were subdivided in two parts: one containing information regarding amide groups, alkanes, carbon-oxygen ester bonds (between 1700 and 1000 cm −1 ), and another with wavenumbers characterized with absorption peaks typical of lipid sources (3000-2800 cm −1 and 1765-1730 cm −1 ) [22,39].
Figure 4 shows the resulting dendrogram of the cluster analysis that compared the different cheese types and ripening times to obtain the similarities regarding their spectra.It was possible to observe a clear separation of the early stage samples among the 3 cheese types.A small cluster containing the earliest Gouda samples was the one that revealed to be less similar.Regarding the ewes' ripe cheese, the 1st month of maturation and the 12th month were related despite their differences.The other samples of ewes' ripe cheese, in between the 1st and 12th month, were closely related with Pecorino samples and the 11th maturation month of Gouda.These samples were all grouped together in a large cluster.The similarities between the oldest Pecorino (11 months of maturation) and Gouda (11 months of maturation), with the intermediate months of the ewes' ripe variety, suggest that the different ways that the cheeses were processed did Figure 4: Hierarchical cluster analysis of the 15 samples of Pecorino (PEC), Gouda (GOU), and ewes' ripe (R) varieties of cheese.The ripening times are represented by its number followed by M (months).
not affect the physicochemical reactions driving the maturation of the cheeses.PCA analysis resulted in PC1 and PC2 describing above 75% of total variance of the spectra (Figure 5(a)).A separation between the different cheeses and their ripening times agreed with the pattern obtained with cluster analysis.In relation to PC1 (54.1% of total variance), the main difference was regarding the Gouda and ewes' ripe cheeses.PC2 accounted for 21.9% of the variance of the data.In the loading plot (Figure 5(b)), the spectral region linked to lipids (3000-2800 cm −1 ) and other chemical groups related to carbohydrates and to the fingerprint of the samples (1280-1000 cm −1 ) weighted positively to PC1, while the region associated to ester groups (1330-1280 cm −1 ) and proteins (1700-1500 cm −1 ) was present in opposition.The PC2 on the other hand seemed to be more related to bonds linked to proteins, where more intense loadings were obtained.
Regarding Pecorino cheese, the samples of the 2nd and 3rd months of maturation were correlated to both PCs.Next to them was also located at the sample of the 4th month.The amide I and II regions were the main responsible for that proximity according to the loading plot (Figure 5(b)).The higher levels of proteolysis during the ripening of this cheese led to a decrease in number of peptide bonds (Figure 2(a)).The samples of the 5th and 11th months of Pecorino cheese were not well correlated in relation to PC1, and there was a significant difference between the earliest samples in relation to PC2.The 5th-month sample revealed high absorbance The altered profile observed in this sample might be related, as already mentioned, to the changes in the milk fat level during the lactation period.The Gouda cheese samples were all well correlated in relation to the PC2, while on PC1, the most variability between them could be explained.PCA validated the spectral observations, where the protein region was not significantly different as in the other cheeses (Figures 2(b) and 5(a)).In case of the peaks linked to lipids (Figure 3(b)), the Gouda samples showed higher absorbance values compared to the other cheeses (Figures 3(a) and 3(c)).This fact corroborates with the positioning of these samples on the positive scale of PC1.The positioning of the 5th month of Gouda might be related with the same issue observed on Pecorino's.
In the ewes' ripe cheeses (Figure 5(a)), the lipolysis was the main reason for the variability of the samples.The lipid region, which was linked to PC1, showed a reduction pattern, except for the sample of the 3rd month that might have also been affected by the fat content deviation during the lactation stage (Figure 3(c)).Regarding PC2, the sample of the 1st and 3rd months, on the other hand, revealed the highest peaks in the protein region (Figure 2(c)).These wavenumbers, once again, weighted significantly to the earliest samples, just as in Pecorino cheese.
Different authors also have demonstrated discrimination between different cheeses by means of PCA [29,32,40].They have found a link between ripening and absorbance changes of amide I and II bands studying the earliest stages of ripening; however, the changes observed were not as evident as the ones found in this study.
Our results agreed to the fact that proteolysis and lipolysis are the main factors driving the maturation of cheeses.A timeline could be observed in relation to both principal components, and a trendline of the ripening process was drawn in Figure 5(a).It was concluded that the proteolysis was more intense during the earliest stages up to the end of the first trimester.Thus, different cheese types and maturation times resulted in different spectral information which must be analyzed with caution during the association with ripening index of cheeses.

Conclusions
In this study, we have interpreted and differentiated cheeses during long maturation times by means of midinfrared spectroscopy technique (FTIR-MIR-ATR).The main factors influencing the ripening process (lipolysis and proteolysis) were stressed by multivariate analysis.The infrared spectroscopy proved to be efficient on following the changes in time occurring in the cheeses regarding proteolysis and lipolysis during the maturation period.Therefore, FTIR shows potential to assist in the quality control of dairy products.While evaluating the maturation degree, the method can be envisaged with the robustness of the chemometric tools.Implemented in dairy plants, the FTIR/ATR will allow automation of the cheese ripening process without the need of sample treatment as well as subsidize a better understanding about their shelf life.

Table 1 :
Physicochemical characterization of Pecorino (2nd month), Gouda (2nd month), and ewes' ripe (1st month) cheeses.Values obtained from duplicates and expressed as percentages ± standard deviation.FDM: fat in dry matter.Superscript letters were used to indicate the group classification.