Next Article in Journal
CO2-Driven N-Formylation/N-Methylation of Amines Using C-Scorpionate Metal Complexes
Previous Article in Journal
Matrix Selection Strategies for MALDI-TOF MS/MS Characterization of Cyclic Tetrapyrroles in Blood and Food Samples
Previous Article in Special Issue
Antimony in Polyethylene Terephthalate-Bottled Beverages: The Migration Puzzle
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 29
Issue 4
10.3390/molecules29040869
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elemental Fingerprinting of Pecorino Romano and Pecorino Sardo PDO: Characterization, Authentication and Nutritional Value

by
Andrea Mara
1,*,
Marco Caredda
2,
Margherita Addis
2,
Francesco Sanna
3,
Mario Deroma
4,
Constantinos A. Georgiou
5,6,
Ilaria Langasco
1,
Maria I. Pilo
1,
Nadia Spano
1 and
Gavino Sanna
1,*
1
Department of Chemical, Physical, Mathematical and Natural Sciences, University of Sassari, Via Vienna 2, I-07100 Sassari, Italy
2
Department of Animal Science, Agris Sardegna, S.S. 291 Sassari-Fertilia, Km. 18,600, I-07040 Sassari, Italy
3
Department of Environmental Studies, Crop Protection and Production Quality Agris Sardegna, Viale Trieste 111, I-09123 Cagliari, Italy
4
Department of Agriculture, University of Sassari, Viale Italia, 39A, I-07100 Sassari, Italy
5
Chemistry Laboratory, Department of Food Science and Human Nutrition, Agricultural University of Athens, 75 Iera Odos, 118 55 Athens, Greece
6
FoodOmics.GR Research Infrastructure, Agricultural University of Athens, 118 55 Athens, Greece
*
Authors to whom correspondence should be addressed.
Molecules 2024, 29(4), 869; https://doi.org/10.3390/molecules29040869
Submission received: 15 January 2024 / Revised: 6 February 2024 / Accepted: 8 February 2024 / Published: 16 February 2024
(This article belongs to the Special Issue Food Analysis in the 21st Century: Challenges and Possibilities)
Browse Figures
Versions Notes

Abstract

:
Sardinia, located in Italy, is a significant producer of Protected Designation of Origin (PDO) sheep cheeses. In response to the growing demand for high-quality, safe, and traceable food products, the elemental fingerprints of Pecorino Romano PDO and Pecorino Sardo PDO were determined on 200 samples of cheese using validated, inductively coupled plasma methods. The aim of this study was to collect data for food authentication studies, evaluate nutritional and safety aspects, and verify the influence of cheesemaking technology and seasonality on elemental fingerprints. According to European regulations, one 100 g serving of both cheeses provides over 30% of the recommended dietary allowance for calcium, sodium, zinc, selenium, and phosphorus, and over 15% of the recommended dietary intake for copper and magnesium. Toxic elements, such as Cd, As, Hg, and Pb, were frequently not quantified or measured at concentrations of toxicological interest. Linear discriminant analysis was used to discriminate between the two types of pecorino cheese with an accuracy of over 95%. The cheese-making process affects the elemental fingerprint, which can be used for authentication purposes. Seasonal variations in several elements have been observed and discussed.

Graphical Abstract

1. Introduction

Dairy products and milk are among the most valuable foods due to their high nutritional value. Milk is the primary food of mammals at birth and contains, on average 4.9% of lactose, 6.2% of proteins, 7.9% of fats, 1.93 g kg−1 of calcium, 1.58 g kg−1 of phosporus and vitamins (vitamin A, 146 IU; vitamin D, 0.18 IU) [1]. Dairy products are derived from the milk of major and minor ruminant species, such as cows, buffaloes, sheep, and goats. In addition to their nutritional properties, dairy products are important for the economy and traditions of many countries. Among them, Italy is one of the most recognized in the world to produce protected designation of origin (PDO) products. The most widely renowned dairy products include Parmigiano Reggiano, Grana Padano, and Pecorino Romano.
Pecorino Romano (PR) is a sheep cheese primarily produced in Sardinia, an Italian region where the sheep dairy industry is economically relevant [2]. In fact, two other PDO sheep’s milk cheeses are produced here: Pecorino Sardo PDO (PS) and Fiore Sardo PDO. The production of PDO cheeses from sheep’s milk is the main source of income for the livestock industry on the island. Semi-extensive farming is the primary method for rearing milk sheep in Sardinia. This results in cheeses that exhibit unique sensorial properties due to the distinguishing features of pastures and climate.
Despite its relevance to the industry, the economic model is fragile due to its dependence on price fluctuations of PR. This often results in farms failing to cover production costs during times of price drops [2]. For these reasons, recent studies have aimed to enhance economic performance and sustainability of the supply chain [3], develop new marketing strategies, improve farm technologies [4], and diversify the production [5,6]. This latter strategy appears to be the most promising in reducing the reliance of milk costs on PR. Another key method for enhancing product value is to capitalize on consumer awareness of dairy quality and nutritional properties [7,8]. For instance, vitamins and minerals are of significant interest because of their association with various health benefits [7,9].
Elements such as Na, Mg, K, Ca, Fe, Cu, Zn and Se are essential for supporting the immune response, cellular processes, and antioxidant defenses [10,11]. These elements can be obtained naturally through a balanced and tailored diet, which prevents health complications caused by deficiency or excess. On the other hand, toxic elements such as As, Cd, Hg and Pb pose health risks to humans at any concentration [12]. Anthropogenic activities often lead to pollution by toxic elements, which contaminate food through water, soil, and air [13]. To ensure food safety, regulations and safety measures have been implemented to monitor and limit the presence of toxic elements in food. For instance, the European community has established maximum levels of toxic elements in food [14,15]. However, most of the toxic elements in milk and milk products are not currently regulated, such as As, Hg and Cd, while the limit for Pb in such matrices is 0.020 mg kg−1 [14]. Previous studies have investigated the dietary intake of toxic metals from milk and its derivatives [16,17,18].
In addition to nutritional and safety aspect, the elemental composition of foods [19] can provide valuable information for authentication [20], traceability [21], and origin assessment of dairy products [22]. The concentration of elements in foods can be affected by various factors, including climate and translocation from soil, water, and air [19]. Research has shown that elemental fingerprints of dairy products can be used to discriminate their geographical origin [23,24,25], verify their PDO authenticity [20,26,27,28,29], identify breeding methods [30], assess production processes [31], and trace the production chain [32]. Accurate sampling is necessary to achieve these goals and encompass all variables that can influence the elemental fingerprint. This includes seasonality, product processing steps, soil characteristics, pollution sources, and class variability [33].
For several decades, this research group has concentrated on the valorization [34,35,36,37,38], quality protection [39,40,41,42], classification [43], and food safety [44,45,46] of dairy products from Sardinia. Therefore, the potential of elemental analysis in food valorization and authentication was evaluated in this study by measuring the elemental fingerprint of two Sardinian PDO sheep cheeses. The concentrations of 31 elements in 200 samples of PR and PS were determined using inductively coupled plasma-optical emission spectroscopy (ICP-OES) and inductively coupled plasma-mass spectrometry (ICP-MS). The main objective was to evaluate nutritional properties and food safety to enhance products and record data for food authentication studies, e.g., for geographical discrimination. In addition, the effects of cheesemaking and seasonality on elemental fingerprints were evaluated.

2. Results and Discussion

2.1. Elemental Composition of Pecorino Romano PDO and Pecorino Sardo PDO

The elemental analysis of PR and PS was conducted using ICP-OES to determine the concentrations of macroelements (Ca, K, Mg, Na, P, and S), and ICP-MS to measure the amounts of trace elements (Zn, Fe, Mn, Cu, Se, Rb, Sr, Al, B, Co, Ni, Cr, V, Li, and Ag) and toxic elements (As, Cd, Hg, Pb, Sn, Sb, Tl, Te, Bi, and U). The results are presented in Table 1 and are expressed in kg of dry matter in the cheese.
Both cheeses had comparable levels of Ca, Mg, K and P in terms of macroelements, with an order of abundance that reflected their concentration in the original milk: Ca > P > K > S ≥ Mg [1,9]. However, due to the distinct salting process employed, the Na concentration in PR was higher than that in PS. Typically, the NaCl concentration in PR ranges from 3% to 7%, whereas in PS it rarely exceeds 2%. Even in terms of trace elements, both cheeses have a similar elemental content. Consistent with the initial milk composition, the trace elements found in highest abundance were Zn, Fe, and Cu. Both PR and PS contained similar amounts of Se, Rb, Sr, and Al. Other trace elements, such as Co, Ni, Cr, V, Li, and Ag, were present in both cheeses at levels near or below the limit of quantification.
Regarding toxic elements, both cheeses contained low levels of As, Cd and Pb. Hg was never detected. Other toxic elements, such as Sn, Sb, Te, Tl, Bi, and U, were generally either not quantified or present at very low levels. It is worth noting that the level of Te in PS was significantly higher than that in PR. However, the European Food Safety Agency (EFSA) is currently investigating the potential toxicity of Te [47].
Finally, a semi-quantitative analysis of rare earth elements (REEs) was preliminarily performed. The REEs were seldom detected above the instrumental detection limit, with a few exceptions for the LREEs. Further investigations will be conducted to optimize the limits of quantification of the analytical method and enable the determination of markers for traceability of the production chain [32].
To the best of our knowledge, the determination of trace elements in PR and/or PS has rarely been accomplished. Previous literature has mainly focused on quantifying macroelements [28,48,49], with only occasional attention given to trace elements such as Zn, Fe, Se, Cu [48], Ba [28], and Al, Ba, Cd, Co, Cr, Cu, Fe, Mn, Ni, Pb, Pt, Sr, and Zn [49]. The data obtained in this study are in good agreement with those of previous studies. Table S1 enables a comparison of the elemental compositions of PS and PR as measured in this study and in the literature.

2.2. Differentiation Due to Cheese-Making Process Technology

Elemental fingerprinting has been reported in the literature as a method for authenticating cheeses [23,24,26,27,28,50]. This technique is commonly used to discriminate cheeses made from milk of different animal origins [26,27,50], from different geographical areas [23,28], or from significantly different cheese-making processes, such as final moisture content and salting [24,26]. Samples were collected from various dairies in Sardinia (Italy) and varied in cheese-making technologies and production period (seasonality).
Principal Component Analysis (PCA) was used for data visualization. The data was cleaned by removing any elements that were not quantified in at least 90% of the samples or were not significant for the analysis. Additionally, Na was excluded as a variable to eliminate the influence of the salting process. Outliers were identified and removed using T2 and Q statistics after performing a preliminary PCA with p > 0.05. The results of the PCA are presented in Figure 1.
PC1 and PC2 accounted for 24.5% and 19.7% of the total variance, respectively. Figure 1A shows that the most abundant trace elements, such as Zn and Cu, are characterized by positive values of PC1, while macroelements, including Ca, Mg, P, and K, are characterized by positive PC2 values. Looking at the score plot (Figure 1B), positive PC1 values tended to occur in the PS cluster (red samples), which was associated with a higher concentration of trace elements, while the PR cluster (black samples) tended to have negative PC1 values. The differentiation between the two clusters was accentuated upon observing PC3, which explained 12.1% of the variance (3D score plot, Figure 1C). This evidence suggests that the elemental fingerprint may be used to discriminate between the two types of cheese. Linear Discriminant Analysis (LDA) was used for classification. The MANOVA test showed a significant difference between the two groups (F (15, 180), Wilks = 0.175, approx. F = 56.78, p < 0.001). Prior to LDA, the dataset was randomized and split into a training set (n = 140) and a test set (n = 55). The results obtained from cross-validation and prediction are reported in Table 2.
The levels of discrimination achieved in cross-validation (97.1%) and prediction (95.7%) were highly accurate. The elemental fingerprint can discriminate between PR and PS using macro-elements (i.e., Ca, K, Mg, P, and S) and trace elements (i.e., Zn, Fe, Mn, Cu, Se, Rb, Sr, Al, Co, and V). These findings were consistent when LDA was used to analyze data from samples produced by three farms that yielded both PR and PS during the same period. The statistical significance of the data was reduced due to the smaller sample size (n = 62). However, the Principal Component Analysis (PCA) in Figure S1 showed that the samples were distinguishable, and the LDA successfully classified them with an accuracy of 98.1% in cross-validation.
Among the literature reviewed, the study by Di Donato et al. [28] was the most comparable to the present study as it adopted a similar approach for authenticating Italian Pecorino cheese. However, it is important to note that the samples in their study were geographically diverse, collected from three different regions of Italy. In contrast, the differences found in our study were solely attributed to the cheese production method used.
The study confirmed the impact of cheese-making technology on the elemental fingerprint. The differences in the chemical equilibrium of milk’s various elements, as well as the chemical form (ionic soluble or colloidal bound to milk proteins and fats), and the variations in cheese-making techniques (coagulation, wheying, and salting) could explain the distinct elemental fingerprints of the two cheeses analyzed in this research. During cheese production, the main stages that lead to a significant variance in the concentration of elements found in the cheese are salting and coagulation.
Salting triggers an osmotic phenomena, resulting in fluctuations in the levels of unbound minerals, which can cause a loss of water and cationic elements such as Na, K, Al, Cd, Co, and Rb. On the other hand, the elements bound to caseins and fats, such as Ca, P, Fe, Mg, Mn, Ni, Pt, and Zn, become concentrated.
During coagulation, there is a non-uniform distribution of elements between the curd and whey. Na and K are soluble elements and tend to be distributed in the aqueous phase (whey). On the other hand, Ca, Mg and P are associated in different proportions with the colloidal suspension of casein micelles and are more concentrated in the curd during cheese production [51]. Several studies have shown that the soluble phase of sheep’s milk may contain different percentages of Ca, Mg and P, with deviations from the total content ranging from 20–25%, 35–64% and 35–40%, respectively [52,53,54]. Currently, there is limited information available on the distribution of trace elements during sheep milk cheese production. However, data indicate that Zn and Mn are primarily distributed in the curd, accounting for about 90%, while Fe and Cu account for about 70% [52]. Fluctuations in pH, temperature, and milk storage conditions affect the equilibrium between soluble and colloidal forms. Generally, a decrease in pH and temperature shifts the balance towards soluble ionic forms, while an increase in pH and temperature favors solubilization. This results in an increase in the retention of certain elements in the curd during cheese production.
Based on these considerations, the various conditions of acidification (which are more intense in PR compared to PS), curd breaking (which is more extensive in PR compared to PS), cooking (with curd cooking at 45 °C in PR and 43 °C in PS), whey removal, heating, curd cooling, and salting may have induced differences in the balance between curd and whey. As a result, the different cheesemaking technologies favor the retention of certain elements, especially Zn, Fe, Cu, and Mn in PS compared to PR (see Figure 1). Milder acidification and breaking of the curd may have resulted in less demineralization during the production of PS. This could have led to greater retention of trace elements that are partially bound to casein micelles, such as Zn, Mg, Fe, and Cu, as previously observed. Additionally, rapid lactic fermentation, followed by effective whey removal, promotes curd demineralization [55].

2.3. Effect of Seasonality

In Sardinian sheep farming births are synchronized and the lactation period starts in November and ends in June–July. The chemical composition of sheep’s milk changes during this period, depending on the diet, lactation stage, and climate [4]. The concentration of minerals and major classes of compounds can be affected by the lactation stage [56]. Therefore, we evaluated the effect of seasonality on the elemental composition of PR and PS cheeses using PCA and ANOVA.
The PCA analysis in Figure S2 shows a distinct trend for Pecorino Romano PDO cheese. The loading plot indicates that PC1 describes the correlation between trace element concentrations and seasonality. Negative values indicate elements that are more abundant in PR produced in summer (V, Al, Rb, and Fe), while positive values indicate elements that are more abundant in PR produced in winter (Zn and Cu). The ANOVA results confirmed that seasonality had an impact on 12 out of 15 elements (see Figure S3). Winter-produced cheeses had the highest concentrations of Zn, Cu, K, and Mn, while spring-produced cheeses had the highest concentrations of Ca, Mg, P, and S. Summer-produced cheeses had the highest concentrations of Rb, Fe, Al, and V. The concentrations of Se, Na, and Sr were not affected by seasonal variations.
In contrast, the impact of seasonality on Pecorino Sardo PDO is relatively minor. Although the PCA did not reveal any clear trends (see Figure S4), the ANOVA indicates a significant effect of seasonality on 8 out of 15 elements (Zn, Ca, K, P, S, Cu, Rb, and V), as shown in Figure S5. Notably, the trends for these elements were like those observed for PR, with higher concentrations of Ca, P, and Mg found in PS during the spring season. Additionally, the concentration of Cu was highest in winter cheese.
As expected, both types of pecorinos exhibit similar seasonal variations in their composition, reflecting the composition of sheep’s milk. Ca is closely linked with P in casein micelles, which provide the structure and stability of the micelles. Colloidal calcium phosphate links the casein submicelles together, occupying 6% of the micellar structure. Therefore, there is a positive correlation between Ca, P, and casein content in ruminant milk [57]. Both cheeses (PS and PR) exhibit the highest concentrations of Ca, Mg, P, and S during the spring season when the amount of casein in the ewe’s milk reaches its daily maximum level [4]. S is not directly involved in the stabilization of micelles, but it is present in proteins, specifically whey proteins (cysteine and cystine amino acids). Therefore, the higher concentration of sulfur in spring cheeses could be linked to the protein concentration found in the milk of spring sheep. Despite this, the cheese still contains only small amounts of whey proteins. When examining the alkaline elements, no evident trends were found. Na cannot be evaluated due to the salting process. K remained constant in winter and spring but decreased significantly in summer. Rb in PR increased from winter to summer, but this trend was not observed in PS. This difference was likely due to different production methods. Regarding trace elements, there were two opposite trends observed. At the start of the lactation period, the concentrations of Zn, Mn, and Cu were highest. However, towards the end of the lactation period, the concentrations of Fe and V tended to increase.
The elemental composition of sheep’s milk during lactation is subject to seasonal variations, which are influenced by various factors such as the lactation stage, nutritional status of the animal, as well as environmental and genetic factors [51,58]. The mineral content in milk is weakly affected by ruminant feeding because the maternal skeleton tends to demineralize during periods when dietary mineral intake does not meet the mineral requirements of the newborn, thus, compensating for the deficit [59]. Skeletal demineralization typically occurs during periods of high mineral demand, such as early lactation and colostrum production [60]. The influence of the lactation stage on the mineral composition of milk is not well-documented. In bovine milk, Ca, P, Mg, and Na levels tend to increase towards the end of the lactation period [61]. This is likely due to increased permeability of the mammary epithelium as lactation progresses [62]. Finally, the mineral content of milk can also be influenced by the animal’s health status and genetic type. The concentration of most minerals in milk decreases when mastitis is present in the mammary gland, except for sodium and chloride ions, which increase instead [63].

2.4. Nutritional and Safety Aspects

Milk and dairy products are considered highly nutritious. Mineral content is an important factor in determining food value, according to consumer preferences. Cheese is a well-known source of minerals, especially Ca, P, and Mg. Casein peptides in milk or cheese prevent the precipitation of calcium in the intestine, making it easily bioavailable [64]. Although the etiology of osteoporosis is complex, adequate calcium intake during childhood and adolescence is important for developing high peak bone mass. Maximizing bone mass early in life is considered a crucial preventive factor against osteoporosis [65]. This study’s results indicate that Pecorino Romano PDO and Pecorino Sardo PDO are sources of several nutritional elements. Figure 2 displays the daily mineral intake for both cheeses across genders. The data were calculated based on the Dietary Recommended Intakes (DRI) published by the US Department of Health and Human Services, National Institutes of Health [66].
Figure 2 shows that the daily consumption of both cheeses meets the Recommended Dietary Allowances (RDAs) and Adequate Intakes (AIs) for many elements (daily portion of 100 g). According to European guidelines [67,68], both cheeses are rich in Ca, P, Zn, and Se (DRI > 30%). Additionally, PR and PS cheeses are potential sources of Mg (females, DRI > 15%) and Cu (DRI > 15%). The Cu content in cheeses may vary seasonally, as shown in Figures S3–S5. Therefore, it may be possible to produce cheeses with high Cu content during the winter season. These findings are significant as they allow dairies and protection Consortia to implement nutritional labeling in accordance with European regulations [69]. As for toxic elements such as As, Cd, Hg, Pb, and Sn, both PR and PS cheeses have a high level of food safety. Additionally, the concentrations of Tl, Bi, and U were frequently below their limits of quantification. Moreover, the cheeses analyzed in this study were obtained from various locations (see Figure 3), indicating the high level of food safety in Sardinian sheep’s milk production.

3. Materials and Methods

3.1. Samples

A total of 200 samples of Pecorino cheese produced in Sardinia in 2021 were obtained from 16 dairy farms that collected milk from livestock farming in surrounding areas. Two PDO sheep cheeses were considered: PR (n = 103) and PS (n = 97). The production technologies adhere closely to the specifications outlined in their respective consortium regulations [70,71]. Pasteurized whole milk is curdled using cultures of milk enzymes from the milk’s place of origin. To produce PR, the milk is coagulated at 38–40 °C and the curd is cooked at 45–48 °C. The resulting 20–35 kg wheels mature for 5–18 months. On the other hand, PS is produced by coagulating the milk at 35–39 °C, cooking the curd at 43 °C and ripening the 1.7–4.0 kg wheels for 2–6 months. The salting process can occur through both dry and wet methods. For PR dry salting, the side and plate of the cheese receive about 2800 g of NaCl distributed thrice over a maximum span of 50 days. Alternatively, the wheels are immersed in a dynamic brine (at 12 °C and a NaCl concentration of 20–22%) for 20 days. Similarly, PS can also be salted dry or wet, but generally, it is kept in brine (24% of NaCl) for 10 h kg−1 of cheese at a temperature of 10–12 °C.
In this study, cheeses were produced by each dairy in three periods of the year: winter, 37%; spring, 35%; and summer, 28%. Thus, samples differed in seasoning and cheesemaking. Additionally, the cheeses varied in maturation (PR: 5–18 months, PS: 2–6 months). The wheels were first divided according to what reported in previous studies ISO 707:2008 [72]. Then, aliquots obtained were aggregated, homogenized, and stored at a temperature between −18 °C and −24 °C until analysis. Figure 3 shows the selected information on Pecorino cheese samples.

3.2. Instrumentation and Reagents

Elemental analysis was performed on a NexION 350X spectrometer equipped with an S10 autosampler, a glass concentric nebulizer, a glass cyclonic spray chamber, and a kinetic energy discrimination (KED) collision cell, all from Perkin Elmer (Milan, Italy). The most abundant elements were determined using an OPTIMA 7300 DV spectrometer (Perkin Elmer, Waltham, MA, USA) equipped with a GemTip Cross-Flow II nebulizer (Perkin Elmer, Waltham, MA, USA) and an autosampler (SC-2 DX, Elemental Scientific Inc., Omaha, NE, USA). To determine macroelements (Ca, K, Mg, Na, P, and S) cheese samples were previously dried in a drying oven (Memmert, Schwabach, Germany) and then calcined in a muffle furnace (Gelman Instrument, Opera, Italy). To detect trace elements (Zn, Fe, Mn, Cu, Se, Rb, Sr, Al, B, Co, Ni, Cr, V, Li, and Ag) and toxic elements (As, Cd, Hg, Pb, Sn, Sb, Tl, Te, Bi, and U) samples were digested using an ultraWAVE™ microwave single reaction chamber (SCR) system (Milestone, Sorisole, Italy) equipped with a rotor at 15 positions and 15 cm3 polytetrafluoroethylene (PTFE) vessels. In all procedures, type I water (resistivity > 18 MΩ cm−1) was produced using a MilliQ Plus System (Millipore, Milan, Italy). Hydrochloric acid (37% w/w), nitric acid (67–69% w/w, NORMATON® for ultra-trace analysis), hydrogen peroxide (30% w/w, NORMATON® for ultra-trace metal analysis), syringes (HENKE-JECT®, 20 cm3), syringe filters (25 mm diameter, 0.22 μm pores, nylon), paper filters (ashless, Grade 40, Whatman®) and metal-free falcon tubes (15 cm3 and 50 cm3) were purchased from VWR (Milan, Italy). Periodic table mix 1 (TraceCERT®, 33 elements, 10 mg dm−3 in 10% HNO3), periodic table mix 2 (TraceCERT®, 17 elements, 10 mg dm−3 in 5% HCl), periodic table mix 3 (TraceCERT®, 16 elements, 10 mg dm−3 in 5% HNO3), and certified skimmed milk powder ERM-BD151 were obtained from Sigma-Aldrich (St. Louis, MO, USA). Single standard solutions of Sc, Y, Ge, Rh, Ir, Mo, Sb, Sn, Hg, and U (100–1000 mg dm−3 in 2–5% HNO3) were obtained from LabKings (Hilversum, The Netherlands).

3.3. Sample Preparation

Sample preparation for macroelements analysis was performed as previously described [6]. The sample digestion for trace and toxic elements was made using an ultraWAVE™ microwave single reaction chamber (SRC) system (Milestone, Sorisole). With respect to conventional microwave instruments, SRC technology reaches higher temperatures and pressures and can manage higher amounts of samples using lower amounts of reagents [73]. In accordance with previous studies [73,74], nitric acid and hydrogen peroxide were used as oxidizing agents in this study. Approximately 0.450 g of the sample (exactly weighted on the analytical balance) was treated with 1 cm3 of HNO3 (67–69%), 2 cm3 of H2O2 (30%), and 4 cm3 of ultrapure H2O.
The digestion program is listed in Table 3. After cooling, the samples were collected, diluted to 15 cm3 using ultrapure H2O, and filtered using a syringe filter. The final residual acidity determined by titration with 0.1 mol dm−3 sodium hydroxide was 2.5 ± 0.2%. To ensure the quality of the analytical data, each digestion batch included a blank and a sample of certified reference material (CRM), ERM-BD151. The same CRM was used to assess the efficiency of microwave acid digestion in terms of matrix effect and trueness.

3.4. Elemental Analysis, Validation, Quality Control and Assurance

Macroelements (Ca, K, Mg, Na, P, and S) were determined using inductively coupled plasma-optical emission spectroscopy (ICP-OES), whereas trace and toxic elements were analyzed using inductively coupled plasma-mass spectrometry (ICP-MS). The instrumental parameters used for the analysis are reported in Table S2 (ICP-OES) and Table S3 (ICP-MS). Further details regarding the ICP-OES method have been previously reported [6], whereas the ICP-MS method was fully developed and validated in this study. For each PDO cheese, three samples were randomly selected and analyzed using the semi-quantitative TotalQuant® method (Syngistix software v 2.3). This preliminary analysis allowed the assessment of the elements that were always within the instrumental detection limit. Excluding these, the elements of interest for a possible ICP-MS quantification were Ag, Al, As, B, Bi, Cd, Co, Cr, Cu, Fe, Hg, Li, Mn, Ni, Pb, Rb, Sb, Se, Sn, Sr, Te, Tl, U, V, and Zn. Subsequently, the possible presence of polyatomic interferences in the real matrix was ascertained, and for each element, the most suitable analysis mode (STD mode or KED mode) was determined. Validation was accomplished in terms of limits of detection and quantification, precision, and trueness. The validation parameters are listed in Table S4. The limits of detection (LoD) and quantification (LoQ) were calculated according to Currie [75]. Method repeatability (CV%r) was assessed by analyzing samples in triplicate within the same analytical session, whereas intermediate precision (CV%IP) was calculated using data obtained from different analytical sessions. Finally, trueness was evaluated by analyzing certified milk CRM ERM BD-151 and spiking tests. In the last case, for each analyte, samples were spiked three times at increasing concentration levels. The trueness measured by analyzing the CRM (Table S5) was between 89 ± 5% (P) and 114 ± 5% (Na) for the macroelements (ICP-OES method), whereas that measured for the trace and toxic elements (ICP-MS method) ranged from 92 ± 5% (Cd) to 110 ± 5% (Se). Furthermore, the recoveries measured by spiking tests ranged from 86 ± 1% (Ag) to 149± 7% (As). The recovery results (Table S4) show that the determination of more than 80% of the elements was bias-free (criteria: t-test, p = 95%). Moderate underestimation was observed for Ag and Bi, and slight overestimation for Sr and V. However, the observed bias are acceptable according to the AOAC guidelines [76]. Finally, for Na and As, a meaningful overestimation bias was observed.
Quantification was performed by external calibration using single- and multi-element standard solutions in 2.5% HNO3. The calibration was performed according to the expected analyte concentrations. Additionally, measurements were performed in triplicate and the data were blank-corrected. To account for signal fluctuations and matrix effects, Rh (50 µg dm−3) and Ir (1 µg dm−3) were used as internal standards. A 60-s wash with a 2% aqueous solution of HNO3 was introduced between consecutive samples to prevent memory effects.

3.5. Statistical Analysis

Data analysis was performed using R-Studio (v. 4.3.1) and Chemometric Agile Tool (CAT) [77]. The Shapiro–Wilk test was used to confirm the normal distribution of the data. ANOVA and MANOVA tests were used to compare groups, and Tukey’s HSD test was used as a post hoc test. The data set was cleaned by removing all elements that were not quantified in at least 90% of the samples or were not significant for the analysis. Elements retained for chemometric analysis were Mg, Ca, P, K, Zn, Cu, Fe, Mg, Al, Rb, Sr, S, Se, and V. Principal Component Analysis (PCA) was performed for data visualization. Linear Discriminant Analysis (LDA) was used to discriminate samples from different categories. Statistical significance was set at α < 0.05.

4. Conclusions

Pecorino Romano PDO and Pecorino Sardo PDO are two popular and highly regarded sheep’s cheeses. This study aimed to address the lack of knowledge regarding their elemental composition and trace element content. To achieve this, 31 elements (macro, trace, and toxic) were analyzed in a comprehensive sample of cheeses using ICP-based methods. The results showed that these cheeses are rich in essential minerals such as calcium, phosphorus, zinc, and selenium, and could potentially be a source of copper and magnesium.
Additionally, the data allowed for an assessment of the impact of cheesemaking and seasonality on the elemental composition. Linear discriminant analysis showed that the elemental fingerprint can effectively distinguish dairy products based on their production method. This result could be a useful methodology for detecting food fraud involving other cheaper pecorinos.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29040869/s1, Figure S1: Principal component analysis performed on Pecorino Sardo and Pecorino Romano produced by 3 farms in the same period: (a) loading plot; (b) score plot.; Figure S2. Principal component analysis performed on Pecorino Romano samples and 14 elements: (a) loading plot; (b) score plot. Object colored according to seasonality; Figure S3: ANOVA analysis of macro and trace elements in Pecorino Romano PDO as a function of the seasonality; Figure S4: Principal component analysis performed on Pecorino Sardo samples and 14 elements: (a) loading plot; (b) score plot. Object colored according to seasonality; Figure S5: ANOVA analysis of macro and trace elements in Pecorino Sardo as a function of the seasonality; Table S1: Average elemental composition of Pecorino cheeses measured in this study and from literature data; Table S2: Instrumental conditions of the ICP-OES OPTIMA 7300 DV, Perkin Elmer; Table S3: Instrumental conditions of the ICP-MS NexION 350X, Perkin Elmer; Table S4: Validation parameters of the ICP-MS method for the elemental analysis of Pecorino cheeses; Table S5: Analysis of the CRM ERM BD-151 (skimmed milk powder).

Author Contributions

Conceptualization, M.A. and G.S.; methodology, A.M., M.A., C.A.G. and G.S.; software, A.M. and M.C.; validation, A.M., M.C. and G.S.; formal analysis, A.M., M.C., M.A., F.S. and G.S.; investigation, A.M., M.D. and I.L.; resources, M.A. and G.S.; data curation, A.M. and M.C.; writing—original draft preparation, A.M., M.C. and G.S.; writing—review and editing, A.M., M.C., M.A., C.A.G., I.L., M.I.P., N.S. and G.S.; visualization, A.M. and M.C.; supervision, M.A., C.A.G. and G.S.; project administration, M.A. and G.S.; funding acquisition, M.A. and G.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Sardegna Ricerche, “POR FESR—Azione 1.1.4 Sostegno alle attività collaborative di R&S per lo sviluppo di nuove tecnologie sostenibili, di nuovi prodotti e servizi, Progetti Cluster—La Diversificazione di Prodotto nell’ambito del Pecorino Romano DOP—CUP G77D13000010002.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in article and Supplementary Materials.

Acknowledgments

The authors would like to thank the Consorzio per la Tutela del formaggio Pecorino Romano DOP and the Consorzio per la Tutela del formaggio Pecorino Sardo DOP for their availability and collaboration in the involvement of the dairy companies and in the procurement of cheese samples. The authors greatly thank Demetrios Sotirchos and Georgios Danezis for their support.

Conflicts of Interest

The authors declare no conflicts 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. Pereira, P.C. Milk Nutritional Composition and Its Role in Human Health. Nutrition 2014, 30, 619–627. [Google Scholar] [CrossRef] [PubMed]
  2. Pulina, G.; Milán, M.J.; Lavín, M.P.; Theodoridis, A.; Morin, E.; Capote, J.; Thomas, D.L.; Francesconi, A.H.D.; Caja, G. Invited Review: Current Production Trends, Farm Structures, and Economics of the Dairy Sheep and Goat Sectors. J. Dairy Sci. 2018, 101, 6715–6729. [Google Scholar] [CrossRef] [PubMed]
  3. Vagnoni, E.; Franca, A.; Porqueddu, C.; Duce, P. Environmental Profile of Sardinian Sheep Milk Cheese Supply Chain: A Comparison between Two Contrasting Dairy Systems. J. Clean. Prod. 2017, 165, 1078–1089. [Google Scholar] [CrossRef]
  4. Pulina, G.; Atzori, A.S.; Dimauro, C.; Ibba, I.; Gaias, G.F.; Correddu, F.; Nudda, A. The Milk Fingerprint of Sardinian Dairy Sheep: Quality and Yield of Milk Used for Pecorino Romano P.D.O. Cheese Production on Population-Based 5-Year Survey. Ital. J. Anim. Sci. 2021, 20, 171–180. [Google Scholar] [CrossRef]
  5. Camanzi, L.; Arba, E.; Rota, C.; Zanasi, C.; Malorgio, G. A Structural Equation Modeling Analysis of Relational Governance and Economic Performance in Agri-Food Supply Chains: Evidence from the Dairy Sheep Industry in Sardinia (Italy). Agric. Food Econ. 2018, 6, 4. [Google Scholar] [CrossRef]
  6. Lai, G.; Caboni, P.; Piras, C.; Pes, M.; Sitzia, M.; Addis, M.; Pirisi, A.; Scano, P. Development and Chemico-Physical Characterization of Ovine Milk-Based Ingredients for Infant Formulae. Appl. Sci. 2023, 13, 653. [Google Scholar] [CrossRef]
  7. Kraus, A. Development of Functional Food with the Participation of the Consumer. Motivators for Consumption of Functional Products. Int. J. Consum. Stud. 2015, 39, 2–11. [Google Scholar] [CrossRef]
  8. Sajdakowska, M.; Gębski, J.; Guzek, D.; Gutkowska, K.; Żakowska-Biemans, S. Dairy Products Quality from a Consumer Point of View: Study among Polish Adults. Nutrients 2020, 12, 1503. [Google Scholar] [CrossRef]
  9. de la Fuente, M.A.; Juárez, M. Milk and Dairy Products. In Handbook of Mineral Elements in Food; Wiley: Hoboken, NJ, USA, 2015; pp. 645–668. ISBN 9781118654316. [Google Scholar]
  10. Cannas, D.; Loi, E.; Serra, M.; Firinu, D.; Valera, P.; Zavattari, P. Relevance of Essential Trace Elements in Nutrition and Drinking Water for Human Health and Autoimmune Disease Risk. Nutrients 2020, 12, 2074. [Google Scholar] [CrossRef]
  11. Islam, M.R.; Akash, S.; Jony, M.H.; Alam, M.N.; Nowrin, F.T.; Rahman, M.M.; Rauf, A.; Thiruvengadam, M. Exploring the Potential Function of Trace Elements in Human Health: A Therapeutic Perspective. Mol. Cell. Biochem. 2023, 478, 2141–2171. [Google Scholar] [CrossRef]
  12. Jaishankar, M.; Tseten, T.; Anbalagan, N.; Mathew, B.B.; Beeregowda, K.N. Toxicity, Mechanism and Health Effects of Some Heavy Metals. Interdiscip. Toxicol. 2014, 7, 60–72. [Google Scholar] [CrossRef]
  13. de Almeida Ribeiro Carvalho, M.; Botero, W.G.; de Oliveira, L.C. Natural and Anthropogenic Sources of Potentially Toxic Elements to Aquatic Environment: A Systematic Literature Review. Environ. Sci. Pollut. Res. 2022, 29, 51318–51338. [Google Scholar] [CrossRef]
  14. EC No 1881/2006; Setting Maximum Levels for Certain Contaminants in Foodstuffs. European Commission: Luxembourg, 2006; pp. 5–24.
  15. EU No 15/2011; Amending Regulation (EC) No 2074/2005 as Regards Recognised Testing Methods for Detecting Marine Biotoxins in Live Bivalve Molluscs. European Commission: Luxembourg, 2011; pp. 3–6.
  16. Koch, W.; Czop, M.; Iłowiecka, K.; Nawrocka, A.; Wiącek, D. Dietary Intake of Toxic Heavy Metals with Major Groups of Food Products—Results of Analytical Determinations. Nutrients 2022, 14, 1626. [Google Scholar] [CrossRef]
  17. Crupi, R.; Lo Turco, V.; Gugliandolo, E.; Nava, V.; Potortì, A.G.; Cuzzocrea, S.; Di Bella, G.; Licata, P. Mineral Composition in Delactosed Dairy Products: Quality and Safety Status. Foods 2022, 11, 139. [Google Scholar] [CrossRef]
  18. Giri, A.; Bharti, V.K.; Kalia, S.; Kumar, B.; Chaurasia, O.P. Health Risk Assessment of Heavy Metals Through Cow Milk Consumption in Trans-Himalayan High-Altitude Region. Biol. Trace Elem. Res. 2021, 199, 4572–4581. [Google Scholar] [CrossRef]
  19. Zhang, P.; Georgiou, C.A.; Brusic, V. Elemental Metabolomics. Brief. Bioinform. 2018, 19, 524–536. [Google Scholar] [CrossRef]
  20. Danezis, G.P.; Tsagkaris, A.S.; Camin, F.; Brusic, V.; Georgiou, C.A. Food Authentication: Techniques, Trends & Emerging Approaches. TrAC—Trends Anal. Chem. 2016, 85, 123–132. [Google Scholar] [CrossRef]
  21. Aceto, M. The Use of ICP-MS in Food Traceability. In Advances in Food Traceability Techniques and Technologies; Elsevier: Amsterdam, The Netherlands, 2016; pp. 137–164. ISBN 978-0-08-100310-7. [Google Scholar]
  22. Drivelos, S.A.; Georgiou, C.A. Multi-Element and Multi-Isotope-Ratio Analysis to Determine the Geographical Origin of Foods in the European Union. TrAC—Trends Anal. Chem. 2012, 40, 38–51. [Google Scholar] [CrossRef]
  23. de Andrade, B.M.; Margalho, L.P.; Batista, D.B.; Lucena, I.O.; Kamimura, B.A.; Balthazar, C.F.; Brexó, R.P.; Pia, A.K.R.; Costa, R.A.S.; Cruz, A.G.; et al. Chemometric Classification of Brazilian Artisanal Cheeses from Different Regions According to Major and Trace Elements by ICP-OES. J. Food Compos. Anal. 2022, 109, 104519. [Google Scholar] [CrossRef]
  24. Magdas, D.A.; Feher, I.; Cristea, G.; Voica, C.; Tabaran, A.; Mihaiu, M.; Cordea, D.V.; Bâlteanu, V.A.; Dan, S.D. Geographical Origin and Species Differentiation of Transylvanian Cheese. Comparative Study of Isotopic and Elemental Profiling vs. DNA Results. Food Chem. 2019, 277, 307–313. [Google Scholar] [CrossRef] [PubMed]
  25. Tedesco, R.; Villoslada Hidalgo, M.D.C.; Vardè, M.; Kehrwald, N.M.; Barbante, C.; Cozzi, G. Trace and Rare Earth Elements Determination in Milk Whey from the Veneto Region, Italy. Food Control 2021, 121, 107595. [Google Scholar] [CrossRef]
  26. Danezis, G.P.; Pappas, A.C.; Tsiplakou, E.; Pappa, E.C.; Zacharioudaki, M.; Tsagkaris, A.S.; Papachristidis, C.A.; Sotirakoglou, K.; Zervas, G.; Georgiou, C.A. Authentication of Greek Protected Designation of Origin Cheeses through Elemental Metabolomics. Int. Dairy J. 2020, 104, 104599. [Google Scholar] [CrossRef]
  27. Nečemer, M.; Potočnik, D.; Ogrinc, N. Discrimination between Slovenian Cow, Goat and Sheep Milk and Cheese According to Geographical Origin Using a Combination of Elemental Content and Stable Isotope Data. J. Food Compos. Anal. 2016, 52, 16–23. [Google Scholar] [CrossRef]
  28. Di Donato, F.; Foschi, M.; Vlad, N.; Biancolillo, A.; Rossi, L.; D’archivio, A.A. Multi-Elemental Composition Data Handled by Chemometrics for the Discrimination of High-Value Italian Pecorino Cheeses. Molecules 2021, 26, 6875. [Google Scholar] [CrossRef] [PubMed]
  29. Santarcangelo, C.; Baldi, A.; Ciampaglia, R.; Dacrema, M.; Di Minno, A.; Pizzamiglio, V.; Tenore, G.C.; Daglia, M. Long-Aged Parmigiano Reggiano PDO: Trace Element Determination Targeted to Health. Foods 2022, 11, 172. [Google Scholar] [CrossRef] [PubMed]
  30. Rodríguez-Bermúdez, R.; López-Alonso, M.; Miranda, M.; Fouz, R.; Orjales, I.; Herrero-Latorre, C. Chemometric Authentication of the Organic Status of Milk on the Basis of Trace Element Content. Food Chem. 2018, 240, 686–693. [Google Scholar] [CrossRef] [PubMed]
  31. de Oliveira Filho, E.F.; Miranda, M.; Ferreiro, T.; Herrero-Latorre, C.; Castro Soares, P.; López-Alonso, M. Concentrations of Essential Trace and Toxic Elements Associated with Production and Manufacturing Processes in Galician Cheese. Molecules 2022, 27, 4938. [Google Scholar] [CrossRef]
  32. Aceto, M.; Musso, D.; Calà, E.; Arieri, F.; Oddone, M. Role of Lanthanides in the Traceability of the Milk Production Chain. J. Agric. Food Chem. 2017, 65, 4200–4208. [Google Scholar] [CrossRef]
  33. Danezis, G.P.; Georgiou, C.A. Elemental Metabolomics: Food Elemental Assessment Could Reveal Geographical Origin. Curr. Opin. Food Sci. 2022, 44, 100812. [Google Scholar] [CrossRef]
  34. Urgeghe, P.P.; Piga, C.; Addis, M.; Di Salvo, R.; Piredda, G.; Scintu, M.F.; Wolf, I.V.; Sanna, G. SPME/GC-MS Characterization of the Volatile Fraction of an Italian PDO Sheep Cheese to Prevalent Lypolitic Ripening: The Case of Fiore Sardo. Food Anal. Methods 2012, 5, 723–730. [Google Scholar] [CrossRef]
  35. Caredda, M.; Addis, M.; Ibba, I.; Leardi, R.; Scintu, M.F.; Piredda, G.; Sanna, G. Prediction of Fatty Acid Content in Sheep Milk by Mid-Infrared Spectrometry with a Selection of Wavelengths by Genetic Algorithms. LWT 2016, 65, 503–510. [Google Scholar] [CrossRef]
  36. Idda, I.; Spano, N.; Ciulu, M.; Nurchi, V.M.; Panzanelli, A.; Pilo, M.I.; Sanna, G. Gas Chromatography Analysis of Major Free Mono- and Disaccharides in Milk: Method Assessment, Validation, and Application to Real Samples. J. Sep. Sci. 2016, 39, 4577–4584. [Google Scholar] [CrossRef] [PubMed]
  37. Pulinas, L.; Spanu, C.; Idda, I.; Ibba, I.; Nieddu, G.; Virdis, S.; Scarano, C.; Piras, F.; Spano, N.; Sanna, G.; et al. Production of Farmstead Lactose-Free Pecorino Di Osilo and Ricotta Cheeses from Sheep’s Milk. Ital. J. Food Saf. 2017, 6, 1. [Google Scholar] [CrossRef]
  38. Idda, I.; Spano, N.; Addis, M.; Galistu, G.; Ibba, I.; Nurchi, V.M.; Pilo, M.I.; Scintu, M.F.; Piredda, G.; Sanna, G. Optimization of a Newly Established Gas-Chromatographic Method for Determining Lactose and Galactose Traces: Application to Pecorino Romano Cheese. J. Food Compos. Anal. 2018, 74, 89–94. [Google Scholar] [CrossRef]
  39. Piga, C.; Urgeghe, P.P.; Piredda, G.; Scintu, M.F.; Sanna, G. Assessment and Validation of Methods for the Determination of γ-Glutamyltransferase Activity in Sheep Milk. Food Chem. 2009, 115, 1519–1523. [Google Scholar] [CrossRef]
  40. Piga, C.; Urgeghe, P.P.; Piredda, G.; Scintu, M.F.; Sanna, G. Analytical Methods for the Evaluation of α-l-Fucosidase Activity in Sheep Milk. Food Anal. Methods 2010, 3, 17–21. [Google Scholar] [CrossRef]
  41. 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]
  42. 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]
  43. Caredda, M.; Addis, M.; Ibba, I.; Leardi, R.; Scintu, M.F.; Piredda, G.; Sanna, G. Building of Prediction Models by Using Mid-Infrared Spectroscopy and Fatty Acid Profile to Discriminate the Geographical Origin of Sheep Milk. LWT 2017, 75, 131–136. [Google Scholar] [CrossRef]
  44. Zazzu, C.; Addis, M.; Caredda, M.; Scintu, M.F.; Piredda, G.; Sanna, G. Biogenic Amines in Traditional Fiore Sardo Pdo Sheep Cheese: Assessment, Validation and Application of an Rp-Hplc-Dad-Uv Method. Separations 2019, 6, 11. [Google Scholar] [CrossRef]
  45. Guiso, M.F.; Battacone, G.; Canu, L.; Deroma, M.; Langasco, I.; Sanna, G.; Tsiplakou, E.; Pulina, G.; Nudda, A. Essential and Toxic Mineral Content and Fatty Acid Profile of Colostrum in Dairy Sheep. Animals 2022, 12, 2730. [Google Scholar] [CrossRef] [PubMed]
  46. Spano, N.; Bortolu, S.; Addis, M.; Langasco, I.; Mara, A.; Pilo, M.I.; Sanna, G.; Urgeghe, P.P. An Analytical Protocol for the Differentiation and the Potentiometric Determination of Fluorine-Containing Fractions in Bovine Milk. Molecules 2023, 28, 1349. [Google Scholar] [CrossRef] [PubMed]
  47. Kowalczyk, E.; Givelet, L.; Amlund, H.; Sloth, J.J.; Hansen, M. Risk Assessment of Rare Earth Elements, Antimony, Barium, Boron, Lithium, Tellurium, Thallium and Vanadium in Teas. EFSA J. 2022, 20, e200410. [Google Scholar] [CrossRef] [PubMed]
  48. Manuelian, C.L.; Currò, S.; Penasa, M.; Cassandro, M.; De Marchi, M. Characterization of Major and Trace Minerals, Fatty Acid Composition, and Cholesterol Content of Protected Designation of Origin Cheeses. J. Dairy Sci. 2017, 100, 3384–3395. [Google Scholar] [CrossRef]
  49. Coni, E.; Bocca, B.; Caroli, S. Minor and Trace Element Content of Two Typical Italian Sheep Dairy Products. J. Dairy Res. 1999, 66, 589–598. [Google Scholar] [CrossRef]
  50. Danezis, G.; Theodorou, C.; Massouras, T.; Zoidis, E.; Hadjigeorgiou, I.; Georgiou, C.A. Greek Graviera Cheese Assessment through Elemental Metabolomics—Implications for Authentication, Safety and Nutrition. Molecules 2019, 24, 670. [Google Scholar] [CrossRef]
  51. Recio, I.; de la Fuente, M.A.; Juárez, M.; Ramos, M. Bioactive Components in Sheep Milk. In Bioactive Components in Milk and Dairy Products; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2009; pp. 83–104. ISBN 978-0-8138-2150-4. [Google Scholar]
  52. de la Fuente, M.A.; Olano, A.; Juârez, M. Distribution of Calcium, Magnesium, Phosphorus, Zinc, Manganese, Copper and Iron between the Soluble and Colloïdal Phases of Ewe’s and Goat’s Milk. Lait 1997, 77, 515–520. [Google Scholar] [CrossRef]
  53. Polychroniadou, A.; Vafopoulou, A. Variations of Major Mineral Constituents of Ewe Milk during Lactation. J. Dairy Sci. 1985, 68, 147–150. [Google Scholar] [CrossRef]
  54. Pellegrini, O.; Remeuf, F.; Rivemal, M. Evolution of Physico-Chemical Characteristics and Renneting Properties of Ewe’s Milk Collected in the “Roquefort Area”. Lait 1994, 74, 425–442. [Google Scholar] [CrossRef]
  55. Martín-Hernández, M.C.; Amigo, L.; Martín-Álvarez, P.J.; Juárez, M. Differentiation of Milks and Cheeses According to Species Based on the Mineral Content. Eur. Food Res. Technol. 1992, 194, 541–544. [Google Scholar] [CrossRef]
  56. Timlin, M.; Tobin, J.T.; Brodkorb, A.; Murphy, E.G.; Dillon, P.; Hennessy, D.; O’donovan, M.; Pierce, K.M.; O’callaghan, T.F. The Impact of Seasonality in Pasture-Based Production Systems on Milk Composition and Functionality. Foods 2021, 10, 607. [Google Scholar] [CrossRef]
  57. Holt, C. Milk Salts|Interaction with Caseins. In Encyclopedia of Dairy Sciences, 2nd ed.; Fuquay, J.W., Ed.; Academic Press: Cambridge, MA, USA, 2011; pp. 917–924. ISBN 978-0-12-374407-4. [Google Scholar]
  58. Fox, P.; McSweeney, P. Advanced Dairy Chemistry; Springer: Berlin/Heidelberg, Germany, 2009; pp. 351–389. ISBN 978-0-387-84864-8. [Google Scholar]
  59. Woodrow, J.; Sharpe, C.; Fudge, N.; Hoff, A.; Gagel, R.; Kovacs, C.S. Calcitonin Plays a Critical Role in Regulating Skeletal Mineral Metabolism during Lactation. Endocrinology 2006, 147, 4010–4021. [Google Scholar] [CrossRef] [PubMed]
  60. Liesegang, A.; Eicher, R.; Sassi, M.-L.; Risteli, J.; Kraenzlin, M.; Riond, J.-L.; Wanner, M. Biochemical Markers of Bone Formation and Resorption Around Parturition and during Lactation in Dairy Cows with High and Low Standard Milk Yields. J. Dairy Sci. 2000, 83, 1773–1781. [Google Scholar] [CrossRef]
  61. Gulati, A.; Galvin, N.; Lewis, E.; Hennessy, D.; O’Donovan, M.; McManus, J.J.; Fenelon, M.A.; Guinee, T.P. Outdoor Grazing of Dairy Cows on Pasture versus Indoor Feeding on Total Mixed Ration: Effects on Gross Composition and Mineral Content of Milk during Lactation. J. Dairy Sci. 2018, 101, 2710–2723. [Google Scholar] [CrossRef] [PubMed]
  62. Hettinga, K. Lactose in the Dairy Production Chain. In Lactose: Evolutionary Role, Health Effects, and Applications; Elsevier: Amsterdam, The Netherlands; Academic Press: Cambridge, MA, USA, 2019; pp. 231–266. ISBN 978-0-12-811720-0. [Google Scholar]
  63. Frédéric, G. The Minerals of Milk. Reprod. Nutr. Dev. 2005, 45, 473–483. [Google Scholar] [CrossRef]
  64. Ebringer, L.; Ferencik, M.; Krajcovic, J. Beneficial Health Effects of Milk and Fermented Dairy Products. Folia Microbiol. 2008, 53, 378–394. [Google Scholar] [CrossRef] [PubMed]
  65. Theobald, H. Dietary Calcium and Health. Nutr. Bull. 2005, 30, 237–277. [Google Scholar] [CrossRef]
  66. U.S. Department of Health & Human Services. National Institutes of Health Nutrient Recommendations and Databases. 2019. Available online: https://ods.od.nih.gov/HealthInformation/nutrientrecommendations.aspx (accessed on 15 January 2024).
  67. European Union Council Directive 90/496/EEC on Nutrition Labelling Rules of Foodstuffs. 1990. Available online: https://eur-lex.europa.eu/legal-content/EN/ALL/?uri=celex%3A31990L0496 (accessed on 15 January 2024).
  68. European Union Regulation (EC) No. 1924/2006 on Nutrition and Health Claims Made on Foods. 2006. Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2006:404:0009:0025:En:PDF (accessed on 15 January 2024).
  69. European Union Regulation (EC) No. 432/2012 on the Compilation of a List of Permitted Health Claims on Food, Other Than Those Referring to the Reduction of Disease Risks and to the Development and Health of Children. 2012. Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2012:136:0001:0040:en:PDF (accessed on 15 January 2024).
  70. Consorzio per la Tutela del Formaggio Pecorino Romano, Disciplinare Pecorino Romano DOP. 2009. Available online: https://www.pecorinoromano.com/pecorino-romano/disciplinare-e-norme (accessed on 15 January 2024).
  71. Consorzio per la Tutela del Formaggio Pecorino Sardo, Disciplinare Pecorino Sardo DOP. 1996. Available online: https://www.politicheagricole.it/flex/files/9/b/b/D.2baf7e91401d2f673d3e/Disciplinare_pecorino_sardo.pdf (accessed on 15 January 2024).
  72. Subramanian, A.; Rodriguez-Saona, L. Chemical and instrumental approaches to cheese analysis. Adv. Food Nutr. Res. 2010, 59, 167–213. [Google Scholar] [CrossRef]
  73. Mara, A.; Deidda, S.; Caredda, M.; Ciulu, M.; Deroma, M.; Farinini, E.; Floris, I.; Langasco, I.; Leardi, R.; Pilo, M.I.; et al. Multi-Elemental Analysis as a Tool to Ascertain the Safety and the Origin of Beehive Products: Development, Validation, and Application of an ICP-MS Method on Four Unifloral Honeys Produced in Sardinia, Italy. Molecules 2022, 27, 2009. [Google Scholar] [CrossRef]
  74. Astolfi, M.L.; Conti, M.E.; Marconi, E.; Massimi, L.; Canepari, S. Effectiveness of Different Sample Treatments for the Elemental Characterization of Bees and Beehive Products. Molecules 2020, 25, 4263. [Google Scholar] [CrossRef]
  75. Currie, L.A. Nomenclature in Evaluation of Analytical Methods Including Detection and Quantification Capabilities (IUPAC Recommendations 1995). Anal. Chim. Acta 1999, 391, 105–126. [Google Scholar] [CrossRef]
  76. Appendix F, Guidelines for Standard Method Performance Requirements, AOAC Official Methods of Analysis. 2016. Available online: https://www.aoac.org/wp-content/uploads/2019/08/app_f.pdf (accessed on 15 January 2024).
  77. Leardi, R.; Melzi, C.; Polotti, G. CAT (Chemometric Agile Tool Software). Available online: http://gruppochemiometria.it/index.php/software (accessed on 15 January 2024).
Molecules 29 00869 g001
Figure 1. Principal component analysis performed on data obtained from the elemental determination of 14 elements on 196 pecorino samples: (A) loading plot; (B) score plot; (C) 3D score plot. Object colored according to cheese type.
Figure 1. Principal component analysis performed on data obtained from the elemental determination of 14 elements on 196 pecorino samples: (A) loading plot; (B) score plot; (C) 3D score plot. Object colored according to cheese type.
Molecules 29 00869 g001
Molecules 29 00869 g002
Figure 2. Mineral daily intakes for adult males and females (17–70 years old) for the consumption of 100 g of Pecorino Romano PDO and Pecorino Sardo PDO.
Figure 2. Mineral daily intakes for adult males and females (17–70 years old) for the consumption of 100 g of Pecorino Romano PDO and Pecorino Sardo PDO.
Molecules 29 00869 g002
Molecules 29 00869 g003
Figure 3. Description of pecorino sampling in terms of dairies, samples, and period of production.
Figure 3. Description of pecorino sampling in terms of dairies, samples, and period of production.
Molecules 29 00869 g003
Table 1. Elemental analysis of Pecorino Sardo PDO (n = 97) and Pecorino Romano PDO (n = 103).
Table 1. Elemental analysis of Pecorino Sardo PDO (n = 97) and Pecorino Romano PDO (n = 103).
ElementPecorino Sardo PDOPecorino Romano PDO
MinMean ± st.dev.MaxMinMean ± st.dev.Max
Macro
(mg kg−1)		Ca10,00014,000 ± 100016,00010,00014,000 ± 100016,000
K10001300 ± 20016007001000 ± 1001300
Mg600700 ± 50800500600 ± 40800
Na50008000 ± 100011,00017,00025,000 ± 500038,000
P80009000 ± 50010,00070009000 ± 70010,000
S500700 ± 10010005001000 ± 2001400
Trace elements
(μg kg−1)		Zn37,00056,000 ± 900078,00018,00047,000 ± 750060,000
Fe26007000 ± 300014,30024006000 ± 9508500
Mn560850 ± 1001250310800 ± 1001150
Cu6001200 ± 50021005001000 ± 3501800
Se210340 ± 90530220400 ± 100580
Rb10001700 ± 50026006001600 ± 5002600
Sr780013,400 ± 250017,900440014,300 ± 250020,500
Al2006000 ± 300013,90025006000 ± 200011,300
B<542000 ± 20007100<548000 ± 800030,000
Co0.94 ± 191.34 ± 16
Ni<1030 ± 1060<1027 ± 540
Cr<3.140 ± 2095<3.120 ± 1050
V510 ± 521715 ± 524
Li<55<55<55<55<55<55
Ag<1.65 ± 58<1.65 ± 510
Toxic elements
(μg kg−1)		As<3.36 ± 18.45.78 ± 110.8
Cd0.51 ± 0.51.70.51.2 ± 0.51.9
Hg<30<30<30<30<30<30
Pb<3.420 ± 1040<3.420 ± 1045
Sn<2.420 ± 1054<2.410 ± 1032
Sb<3.612 ± 516<3.610 ± 518
Tl<0.51.9 ± 0.52.3<0.5<0.5<0.5
Te<1.2130 ± 50220<1.29 ± 515
Bi<0.5<0.5<0.5<0.52 ± 12.8
U<0.191 ± 15.2<0.192 ± 17.1
Table 2. Results of the LDA performed for the discrimination based on cheese type. Confusion matrix and accuracy in cross-validation (training) and prediction (testing).
Table 2. Results of the LDA performed for the discrimination based on cheese type. Confusion matrix and accuracy in cross-validation (training) and prediction (testing).
Confusion Matrix
TrainingRomanoSardoTestingRomanoSardo
Romano673Romano320
Sardo169Sardo221
Accuracy
RomanoSardoTotalRomanoSardoTotal
95.7%98.6%97.1%100%91.3%95.7%
Table 3. Cheese digestion conditions using an ultraWAVE™ SRC system (Milestone).
Table 3. Cheese digestion conditions using an ultraWAVE™ SRC system (Milestone).
StepTime (min)Temperature (°C)
1Heating25240
2Holding10240
3Coolingca. 30<40
Initial pressure: 4 MPa; Release pressure rate: 0.8 MPa min−1; Rotor: 15 positions; Vessels: 15 cm3 (PTFE), Method: 0.450 g cheese + 1 cm3 HNO3 (67–69%) + 2 cm3 H2O2 (30%) + 4 cm3 H2O.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mara, A.; Caredda, M.; Addis, M.; Sanna, F.; Deroma, M.; Georgiou, C.A.; Langasco, I.; Pilo, M.I.; Spano, N.; Sanna, G. Elemental Fingerprinting of Pecorino Romano and Pecorino Sardo PDO: Characterization, Authentication and Nutritional Value. Molecules 2024, 29, 869. https://doi.org/10.3390/molecules29040869

AMA Style

Mara A, Caredda M, Addis M, Sanna F, Deroma M, Georgiou CA, Langasco I, Pilo MI, Spano N, Sanna G. Elemental Fingerprinting of Pecorino Romano and Pecorino Sardo PDO: Characterization, Authentication and Nutritional Value. Molecules. 2024; 29(4):869. https://doi.org/10.3390/molecules29040869

Chicago/Turabian Style

Mara, Andrea, Marco Caredda, Margherita Addis, Francesco Sanna, Mario Deroma, Constantinos A. Georgiou, Ilaria Langasco, Maria I. Pilo, Nadia Spano, and Gavino Sanna. 2024. "Elemental Fingerprinting of Pecorino Romano and Pecorino Sardo PDO: Characterization, Authentication and Nutritional Value" Molecules 29, no. 4: 869. https://doi.org/10.3390/molecules29040869

APA Style

Mara, A., Caredda, M., Addis, M., Sanna, F., Deroma, M., Georgiou, C. A., Langasco, I., Pilo, M. I., Spano, N., & Sanna, G. (2024). Elemental Fingerprinting of Pecorino Romano and Pecorino Sardo PDO: Characterization, Authentication and Nutritional Value. Molecules, 29(4), 869. https://doi.org/10.3390/molecules29040869

Article Metrics

Back to TopTop