Greek Graviera Cheese Assessment through Elemental Metabolomics-Implications for Authentication, Safety and Nutrition.

This study presents the comprehensive elemental profile of Greek Graviera (Gruyère) cheeses. In total, 105 samples from nine different geographic regions produced from sheep, goat and cow milk and their mixtures were assessed. Elemental signatures of 61 elements were investigated for determination of geographic origin and milk type. Regional and milk type classification through Linear Discriminant Analysis was successful for almost all cases, while a less optimistic cross validation exercise presented lower classification rates. That points to further research using a much larger sample set, increasing confidence for cheese authentication utilizing also bioinformatics tools under development. This is the first study reporting signatures of 61 elements in dairy products including all sixteen rare earth elements and all seven precious metals. Safety and quality were assessed regarding toxic and nutritive elements. According to both EU and USA regulations and directives, Graviera is a nutritional source for trace and macro elements with low levels of toxic elements.

1. Introduction

Food pan class="Chemical">authentication’s impn>ortance is inpan class="Chemical">creasing due to the consumers’ interest in accurate food labeling, forcing producers and retailers to follow. Regulatory authorities are interested in analytical methods for food authenticity to support law enforcement [1,2]. Dairy products play a central role in a nutritious and balanced diet and their consumption has been associated with several health benefits due to their high contents of protein, essential fatty acids and minerals [3].

“Graviera” (Gruyère) is a hard type cheepan class="Chemical">se, holding the n>an class="Chemical">second place, after “Feta”, in the Greek dairy production volume. It is produced mainly from a mixture of sheep and goat milk although Gravieras can be found made of solely sheep, or goat or cow or a mixture of all three kinds of milk. Most of the 700 Greek cheese-making dairies are small-medium size, collecting milk over a radius of about 30 km, although very few can collect milk over substantially longer distances. The majority of dairy sheep and goat farms are also of small to medium size, grazing their animals near the farm [4]. Most of the Gravieras are commercialized with a geographical denomination, but only three of them are registered under the Protected Designation of Origin (PDO) EU scheme, namely: “Graviera Agrafon”, “Graviera Kritis” and “Graviera Naxou”. However, the variety of local climates, the highly rich flora of endemic plants, the predominant microflora, and the processing of the milk along with the traditional cheese-making practices are crucial for Graviera’s quality characteristics [5].

pan class="Chemical">Authenticity of dairy products can n>an class="Chemical">be assessed by various analytical techniques such as molecular [6], chromatographic [7], vibrational and fluorescence spectroscopy [8], elemental fingerprinting [9,10,11], isotopic [12], non-chromatographic mass spectrometry [13] and Nuclear Magnetic Resonance (NMR) [14]. Cheese authentication through elemental fingerprinting has been highlighted in several articles [9,10,11,15,16,17,18]. However, these publications do not take into account the comprehensive elemental profile as proposed by elemental metabolomics, including Rare Earth Elements (REE) and precious metals [19].

Results from earlier studies revepan class="Chemical">al a relation n>an class="Chemical">between the elemental profile of soil and the derived food products [15,20]. The elemental content of animal products depends, beyond feed-vegetation intake, on various factors such as animal species (e.g., cow, sheep or goat), mineral supplementation, drinking water and production practices. Moreover, the elemental profile of soil and mineral pollution are associated with each specific geographical area and characterize the origin [9,12,21,22].

Elementpan class="Chemical">al transfer along the food chain is a complex procedure that is also element specific. Elemental fingerprint in the first part of the food chain, i.e., the “plant”, might be substantially different from the elemental fingerprint in the last part, i.e., milk, meat or cheese. That is due to the impact of the aforementioned factors, all metabolic processes and different absorption rates in different organisms. In this respect, the REE fingerprint is directly linked to the geology of the area and could be minimally affected by other factors. REEs have been proven to be very reliable and authentic markers in various products, like split peas [23], wines [24] and game meat [25] and are minimally affected by harvesting variations [26].

Elementpan class="Chemical">al metabolomin>an class="Chemical">cs is emerging as a new technology with applications in various fields such as nutrition, agriculture and links food science to health [19]. Inductively coupled plasma mass spectroscopy (ICP-MS) is the main choice due to its capabilities for rapid ultra-trace level multi-element determinations. Basic principles of elemental metabolomics include but are not limited to proper sample preparation for elemental analysis, use of standard reference materials during ICP-MS analysis and proper data handling, statistical analysis and reporting. Elemental metabolomics Standard Operation Principles (SOPs), as well as the minimal information reporting standards (MIRSs), are reported in literature [19].

This study determined, through ICP-MS anpan class="Chemical">alyn>an class="Chemical">sis, the elemental signatures of 61 elements of Greek Graviera cheese with the aim to investigate the use of the elemental metabolome for the assessment of geographical origin, safety and nutritional quality of this type of cheese.

2. Results & Discussion

2.1. Elemental Concentrations in Greek Graviera and Implications for Authentication

2.1.1. Geographical Origin

The databpan class="Chemical">ase of elementn>an class="Chemical">al signatures, comprised of 61 elements, is shown Table 1, Table 2, Table 3 and Table 4. Table 1 presents rare earth elements & actinides (16 & 2 = 18 elements), Table 2 precious metals & ultra-trace elements (7 & 8 = 15 elements), Table 3 trace elements (15 elements) and Table 4 trace elements of high abundance & macro elements (9 & 4 = 13 elements). Ten rare earth elements (Dy, Er, Eu, Nd, Pr, Sc, Sm, Y, Yb) and one actinide (U) showed statistically significant different values between Greek regions (administrative regions). The same was found for 4 precious metals (Au, Pd, Re, Ru) and 5 ultra-trace elements (Nb, Ta, Tl, W, Zr). All 28 trace and macro elements showed statistically significantly different values between Greek regions in addition to Ag, Al, Bi, Cd, Cu, Mo, Ni and Pb.

Table 1

Rare earth elements and actinides in Greek Graviera, μg kg−1, mean ± Standard Error of the Mean, SEM (samples).

Region		Ce	Dy	Er	Eu	Gd	Ho	La	Lu	Nd	Pr	Sc	Sm	Tb	Tm	Y	Yb	Th	U
Central Greece	Mean	2.5	0.17	0.19	3.1	0.5	0.06	1.7	0.07	1.2	0.33	26	2.9	0.16	0.07	0.8	0.22	1.4	3.7
	SEM (8)	0.5	0.04	0.04	1.0	0.1	0.01	0.4	0.02	0.2	0.07	7	0.7	0.04	0.02	0.2	0.04	0.2	0.4
Crete	Mean	4.1	0.24	0.26	1.0	0.9	0.08	2.7	0.08	2.0	0.5	27	1.0	0.21	0.11	1.1	0.29	1.7	2.6
	SEM (22)	0.9	0.04	0.05	0.1	0.2	0.01	0.6	0.01	0.6	0.1	4	0.1	0.03	0.02	0.4	0.06	0.1	0.3
Epirus	Mean	2.8	0.24	0.31	1.1	0.69	0.086	1.9	0.092	1.6	0.43	34	0.89	0.18	0.097	0.9	0.	1.56	2.2
	SEM (9)	0.3	0.02	0.02	0.09	0.05	0.006	0.2	0.005	0.1	0.04	2	0.06	0.01	0.009	0.1	0.02	0.07	0.1
Macedonia	Mean	2.2	0.18	0.25	0.7	0.56	0.071	1.3	0.085	1.0	0.24	2	0.8	0.25	0.085	0.51	0.25	1.6	2.9
	SEM (13)	0.4	0.02	0.02	0.1	0.06	0.005	0.3	0.006	0.1	0.03	1	0.1	0.04	0.007	0.08	0.02	0.1	0.2
North Aegean	Mean	2.0	0.147	0.217	0.8	0.48	0.065	1.3	0.080	0.9	0.23	26	0.69	0.23	0.078	0.35	0.21	1.60	2.7
	SEM (8)	0.4	0.007	0.008	0.1	0.05	0.004	0.3	0.003	0.1	0.02	2.0	0.09	0.04	0.005	0.06	0.01	0.07	0.2
Peloponnese	Mean	3.6	0.26	0.28	1.2	0.7	0.09	1.8	0.09	1.5	0.45	28	1.3	0.18	0.08	1.2	0.28	1.7	3.3
	SEM (11)	0.7	0.04	0.04	0.2	0.1	0.01	0.3	0.01	0.2	0.08	5	0.1	0.03	0.01	0.2	0.03	0.2	0.4
South Aegean	Mean	3.7	0.22	0.26	1.1	1.0	0.09	2.3	0.086	1.7	0.41	28	1.1	0.19	0.09	0.9	0.28	1.52	2.4
	SEM (21)	0.8	0.02	0.04	0.1	0.3	0.01	0.4	0.007	0.3	0.07	3	0.2	0.02	0.01	0.2	0.03	0.09	0.3
Thessaly	Mean	2.2	0.12	0.13	0.9	0.7	0.06	3	0.07	0.7	0.18	19	1.4	0.12	0.05	0.58	0.16	1.3	4.2
	SEM (7)	0.4	0.02	0.02	0.2	0.3	0.01	2	0.01	0.1	0.02	5	0.3	0.03	0.02	0.08	0.03	0.2	0.6
Thrace	One sample **	1.7	0.13	0.07	0.6	0.3	0.03	0.9	0.04	0.5	0.16	14	1.0	0.06	0.02	0.8	0.12	1.2	4.3
p Value *		0.098	>0.001	0.001	>0.001	0.190	0.075	0.209	0.413	0.002	>0.001	0.028	>0.001	0.236	0.075	0.004	0.007	0.337	>0.001

* p Values > 0.05 mean that the element concentration is not statistically different between regions. ** One digestion, measured in triplicate.

Table 2

Precious metals and ultra-trace elements in Greek Graviera, μg kg−1, mean ± Standard Error of the Mean, SEM (samples).

Region		Au	Ir	Pd	Pt	Re	Rh	Ru	Hf	Nb	Sb	Sn	Ta	Tl	W	Zr
Central Greece	Mean	6	0.9	2.1	1.4	1.6	5	10	0.5	1.2	6	12	0.9	2.0	6	7
	SEM (8)	1	0.1	0.4	0.3	0.4	2	2	0.1	0.3	1	4	0.3	0.6	2	1
Crete	Mean	4.0	0.65	1.7	1.6	0.7	5.6	6	0.52	1.2	6	10	0.7	1.32	7	5
	SEM (22)	0.8	0.07	0.3	0.2	0.2	0.8	1	0.06	0.2	1	2	0.2	0.08	1	1
Epirus	Mean	3.8	0.68	1.7	1.73	2.2	4.3	5.3	0.56	1.30	6.2	13	0.50	1.5	7.6	5.4
	SEM (9)	0.4	0.04	0.1	0.07	0.6	0.4	0.2	0.04	0.09	0.4	2	0.03	0.1	0.5	0.4
Macedonia	Mean	3.8	0.64	1.34	1.6	0.6	5	3.8	0.49	0.6	4.4	7	0.44	1.4	4.4	3.6
	SEM (13)	0.5	0.04	0.09	0.1	0.2	2	0.7	0.03	0.1	0.2	1	0.08	0.1	0.6	0.3
North Aegean	Mean	3.1	0.57	1.37	1.5	0.5	6	4.0	0.47	0.89	5.0	9	0.50	1.31	5.4	3.8
	SEM (8)	0.3	0.04	0.09	0.1	0.1	1	0.5	0.02	0.09	0.2	2	0.07	0.05	0.5	0.3
Peloponnese	Mean	6	0.9	2.9	1.5	0.8	4.7	7.1	0.6	1.3	6.2	13	0.6	1.5	6.0	6.0
	SEM (11)	2	0.1	0.2	0.2	0.2	0.9	0.8	0.1	0.2	0.7	4	0.1	0.1	1	0.5
South Aegean	Mean	3.7	0.64	1.8	1.8	0.8	5.1	5.8	0.57	1.3	6.1	11	0.51	2.1	6.5	5.6
	SEM (21)	0.6	0.04	0.2	0.1	0.1	0.7	0.9	0.05	0.1	0.6	2	0.07	0.3	0.9	0.8
Thessaly	Mean	7	0.7	2.2	1.4	1.2	4	11	1.1	1.8	5.8	13	3	2.1	4	7
	SEM (7)	1	0.1	0.3	0.3	0.4	1	2	0.6	0.7	0.8	3	2	0.2	1	1
Thrace	One sample **	5	0.4	1.4	0.5	0.2	1	11	0.1	0.8	4.6	9	1	1.3	0.9	6
p Value *		0.005	0.075	0.015	0.389	0.026	0.845	>0.001	0.111	0.029	0.217	0.419	0.002	0.013	0.022	0.003

* p Values > 0.05 mean that the element concentration is not statistically different between regions. ** One digestion, measured in triplicate.

Table 3

Trace elements in Greek Graviera, μg kg−1, mean ± Standard Error of the Mean, SEM (samples).

Region		Ag	As	Ba	Bi	Cd	Co	Cr	Cs	Cu	Ga	Mo	Ni	Pb	Se	V
Central Greece	Mean	3.4	25	3000	17	5.2	80	650	14	800	35	100	430	38	110	500
	SEM (8)	0.6	20	500	6	0.7	20	40	3	100	4	30	40	17	20	70
Crete	Mean	4.4	212	1000	16	5.8	39	590	4.5	780	13	110	370	28	65	390
	SEM (22)	0.7	8	100	2	0.8	6	30	0.7	50	1	30	20	7	9	20
Epirus	Mean	4.2	202	920	15.3	6.2	36.0	580	5.3	820	12.4	130	380	27	75	339
	SEM (9)	0.3	6	70	0.6	0.3	0.7	10	0.4	40	0.5	10	20	3	3	9
Macedonia	Mean	2.9	220	1000	14.8	4.3	28	500	3.0	530	12	60	290	21	51	470
	SEM (13)	0.2	20	200	0.8	0.2	5	30	0.4	80	2	10	40	4	9	40
North Aegean	Mean	3.8	220	800	14.5	5.4	31	540	2.8	670	11.6	100	320	20	60	440
	SEM (8)	0.4	9	60	0.7	0.3	3	20	0.3	90	0.8	10	30	2	4	30
Peloponnese	Mean	3.9	240	1000	49	5.8	70	640	5	640	17	110	367	24	100	470
	SEM (11)	0.6	10	200	35	0.6	20	30	1	60	2	20	9	4	20	50
South Aegean	Mean	4.0	220	1000	14	5.8	38	570	6	780	14	110	350	32	63	400
	SEM (21)	0.4	10	100	1	0.4	5	20	2	40	2	20	20	6	8	30
Thessaly	Mean	5	170	1155	14	5.3	70	700	5	800	19	80	410	25	110	340
	SEM (7)	2	30	207	2	0.9	10	100	1	100	4	20	30	7	20	60
Thrace	One sample **	2	240	697	11	3.8	85	730	3	740	17	90	420	20	210	560
p Value *		0.184	0.044	>0.001	0.416	0.140	>0.001	>0.001	>0.001	0.130	>0.001	0.137	0.058	0.432	>0.001	0.007

* p Values > 0.05 mean that the element concentration is not statistically different between regions. ** One digestion, measured in triplicate.

Table 4

Trace elements of high abudance mg kg-1 and macro elements (Ca, K, Mg & P) g kg-1, mean ± Standard Error of the Mean, SEM (samples).

Region		Al	B	Fe	Mn	Rb	Si	Sr	Ti	Zn	Ca	K	Mg	P
Central Greece	Mean	1.6	62	32	0.42	1.5	157	7.6	62	34	8.7	0.56	0.47	7.3
	SEM (8)	0.2	8	3	0.04	0.1	14	0.9	2	2	0.6	0.04	0.02	0.3
Crete	Mean	1.8	48	29	0.39	0.98	134	4.7	62	30	9.9	0.45	0.52	7.0
	SEM (22)	0.3	5	1	0.02	0.07	8	0.5	1	1	0.4	0.02	0.01	0.2
Epirus	Mean	1.3	42.9	31.9	0.44	1.03	126	3.8	68	34.9	10.9	0.56	0.56	7.5
	SEM (9)	0.2	0.6	0.6	0.02	0.06	3	0.3	1	0.8	0.2	0.02	0.01	0.1
Macedonia	Mean	1.3	44	23	0.32	0.8	128	2.4	50	28	7.4	0.7	0.38	5.8
	SEM (13)	0.2	4	2	0.03	0.1	10	0.3	5	4	0.8	0.1	0.03	0.5
North Aegean	Mean	0.9	44	26	0.27	0.63	122	2.5	56	35	8.7	0.68	0.39	6.3
	SEM (8)	0.1	2	1	0.02	0.04	4	0.2	2	1	0.5	0.03	0.02	0.3
Peloponnese	Mean	2	58	29	0.47	0.90	151	3.6	62	33	9.3	0.64	0.50	7.2
	SEM (11)	1	6	1	0.05	0.08	11	0.3	1	1	0.6	0.06	0.02	0.2
South Aegean	Mean	1.8	45	31	0.36	0.88	124	4.2	61	32	10.0	0.51	0.49	6.9
	SEM (21)	0.4	3	1	0.03	0.06	5	0.4	2	1	0.4	0.02	0.03	0.2
Thessaly	Mean	1.6	79	35	0.41	1.02	166	4.3	68	36	9.8	0.66	0.55	8.1
	SEM (7)	0.3	19	3	0.04	0.08	30	0.5	2	1	0.8	0.06	0.03	0.4
Thrace	One sample **	3.0	85	40	0.29	0.65	177	3.7	66	42	7.4	0.70	0.41	8.3
p Value *		0.162	>0.001	>0.001	>0.001	>0.001	0.009	>0.001	>0.001	0.016	>0.001	0.020	>0.001	>0.001

* p Values > 0.05 mean that the element concentration is not statistically different between regions. ** One digestion, measured in triplicate.

pan class="Chemical">Sixty-one elements were upan class="Chemical">sed as predictor variables to develop a method for assigning Graviera origin to nine geographic regions. The classification table (Table 5) shows that one sample from the Epirus region (sample 15, Table S1, from Arta) and another from Thessaly (sample 98, Table S1, from Larisa) were misclassified as coming from Central Greece. However, both Arta and Larisa are adjacent to the Central Greece region. That could provide an adequate explanation for the misclassification as goats and sheep movement between adjacent regions could modify their elemental content and subsequently the milk they produce. However, after further investigation we found that the two dairies in Arta and Larisa, to fulfill their needs, purchased milk from Central Greece, Amfiloxia and Lamia, respectively. The third misclassified sample was from Macedonia region (sample 45, Table S1, from Grevena, prepared from cow milk) classified as South Aegean. Cows in the Macedonia region obtain a large portion of their feed from grass that grows locally that is rich in REEs and aluminum [27] (aluminum ores are usually accompanied by REEs [28]). This is reflected in the content of REEs and aluminum found in Gravieras produced in that region. Sample 45 shows less than 50% REEs and aluminum content compared to all other samples from this region. This points to a different feeding scheme using mostly imported feeds as in the South Aegean. Through a cross validation exercise using the leave-one-out approach, the above results are quite optimistic and the classification rate was just 32.7%. This warrants further research enhancing the sample bank with much more samples resulting from different production periods.

Table 5

Geographical region classification table based on the 61 element signature.

Actual	Group Size	Predicted
		Central Greece	Crete	Epirus	Macedonia	North Aegean	Peloponnese	South Aegean	Thessaly	Thrace
Central Greece	13	13	0	0	0	0	0	0	0	0
		(100.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)
Crete	22	0	22	0	0	0	0	0	0	0
		(0.0%)	(100.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)
Epirus	9	1	0	8	0	0	0	0	0	0
		(11.1%)	(0.0%)	(88.9%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)
Macedonia	13	0	0	0	12	0	0	1	0	0
		(0.0%)	(0.0%)	(0.0%)	(92.3%)	(0.0%)	(0.0%)	(7.7%)	(0.0%)	(0.0%)
North Aegean	8	0	0	0	0	8	0	0	0	0
		(0.0%)	(0.0%)	(0.0%)	(0.0%)	(100.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)
Peloponnese	11	0	0	0	0	0	11	0	0	0
		(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(100.0%)	(0.0%)	(0.0%)	(0.0%)
South Aegean	21	0	0	0	0	0	0	21	0	0
		(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(100.0%)	(0.0%)	(0.0%)
Thessaly	7	1	0	0	0	0	0	0	6	0
		(14.3%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(0.0%)	(85.7%)	(0.0%)

Percent of cases correctly classified: 95.9%.

Elementpan class="Chemical">al metabolomin>an class="Chemical">cs has potential for detecting production method (feed with pasture vs imported/dried feeds). This needs further research with feeding experiments. A useful aspect of elemental metabolomics applied to dairy products could be a bioinformatics tool to detect the feeding scheme utilizing soil composition analysis. Details on the tool/algorithm can be extracted from the discussion above on the Macedonia region’s cows and soil.

The most pan class="Chemical">significant predictor variables are the rare earths Ce, Er, Eu, Ho, La, Sm, Tm, Yb, the n>an class="Chemical">actinide Th, the precious metals Pt, Re, Ru, the ultra trace elements Hf, Nb, Sb, W, the trace elements Ag, Al, As, B, Ba, Cd, Co, Cr, Cs, Fe, Ga, Hf, Mo, Nb, Ni, Pb, Sb, Se, Sr, V, W, Zn and the macro elements Ca, Mg and P The classification was more successful using the comprehensive elemental signature as proposed by elemental metabolomics [19]. This result on cheeses is in contrast to previous studies on authentication of game meat [25], wines [24] and split-peas [23], where specific groups of elements such as REEs were sufficient.

Levels of pan class="Chemical">REEs were higher in Crete cheeses, probably reflecting the vegetation and soil composition [21]. This is most pronounced for the light REEs (LREEs) Pr and Nd (Table 1). This is in accordance with previous findings that Crete is enriched in LREEs, due to monazite and allanite ores [29]. It is interesting to note that another couple of LREEs, Eu & Sm, were enriched, by three times, in cheese from Central Greece, pointing to further authentication markers. These findings about Eu and Sm need further research such as a check of soil composition differences and the influence of different flora grown there.

Usually n>an class="Chemical">all rare earths in different materials follow the same pattern, i.e., they are all enriched or depleted: fava Santorini’s [26]; Italian milk [30]; mushrooms substrates [31]. However, this pattern is differentiated by genetic factors as seen in two different mushroom species [31]. The production method could also differentiate the pattern as seen in game and farmedrabbits [25].

North Aegean cheepan class="Chemical">ses showed much higher levels of n>an class="Chemical">Rb, Cs (alkali metals) and Sr, Ba (alkaline earth metals) in agreement with previous studies [9,10,11,12,15,18,32,33,34], where alkali and alkaline earth metals, especially Rb, Cs, Sr and Ba, were proven reliable cheese authenticity markers. Another interesting result shown in Table 2 is that the precious metals Au, Pd & Ru were found in higher amounts in Thessaly’s cheeses. Geological data [35] explain the increased content of precious metals transported from the Pindus mountain range by the Pinios River to Thessaly. This is in accordance with the view that precious metals are potential authenticity markers [22].

Our data are in line with previous findings. Camin et pan class="Chemical">al. n>an class="Chemical">also stated [12] that levels of Cu, Mo, Ni, Fe, Mn, Ga and Se showed significant differences between grated hard cheeses. Osorio et al. [15] found that Ag, Ba, Ca, K, Mg, Mn, P and Sr presented different profiles for different Halloumi cheese production locations and also highlighted the potential of Sr for traceability information from soil as it cannot be added from cheese making equipment. All elements commended by Osorio in addition to Ag presented significant differences between regions in our work. Korenovska and Suhaj [11], working with Slovakian, Polish, and Romanian Bryndza cheeses, found also that Cr, Hg, Mn and V along alkali and alkaline earth metals were the best elemental indicators.

In agreement with our results, Pillonel et pan class="Chemical">al. [18] working with Emmentn>an class="Chemical">al cheese support the view that the elemental profile allows the discrimination of close regions of production, where the distances are in the order of a few tenths of kilometers up to 150 km. It should be noted that the distance range for adjacent regions in our study is between 10 to 150 km. This is the first study reporting precious metals, rare earth and ultratrace element assessment for dairy product authentication.

2.1.2. Milk Type

The Gravierpan class="Chemical">as upan class="Chemical">sed in this study were manufactured from mixed sheep and goat milk (78), sheep (10), goat (8) and cow (8) and one from sheep, goat and cow milk. The comprehensive elemental fingerprint (Table S2) of 61 elements was used for the classification according to the milk type. Best markers of milk type were: Bi, Cr, Fe, Mn, Ni, Se, Sr, Zn, Mg and P (Table S2). The classification table (Table 6), shows that all 78 cheeses from sheep and goat milk were correctly categorized. Our data did not contain information on the % percentage of goat milk used. This explains the misclassification of two out of the ten sheep milk Graviera samples into the sheep and goat class. This group is not so well defined and homogenous as even the same producer uses different percentages of goat milk, according to its availability. Only one cow Graviera sample was misclassified into the sheep + goat group. This Graviera was from the Macedonia region (sample 41, Table S1, from Grevena). Through a cross validation exercise using the leave-one-out approach, the above results are quite optimistic and the classification rate was only 50.5%.

Table 6

Classification table of milk type based on the total elemental fingerprint.

Actual Type of Milk	Group Size	Predicted Type of Milk
		Sheep + Goat	Sheep	Goat	Cow
Sheep + goat	78	78	0	0	0
		(100.0%)	(0.0%)	(0.0%)	(0.0%)
Sheep	10	2	8	0	0
		(20.0%)	(80.0%)	(0.0%)	(0.0%)
Goat	8	0	0	8	0
		(0.0%)	(0.0%)	(100.0%)	(0.0%)
Cow	8	1	0	0	7
		(12.5%)	(0.0%)	(0.0%)	(87.5%)

Percent of cases correctly classified: 91.9%.

Previous attempan class="Chemical">pts to cln>an class="Chemical">assify cheeses according to the milk type, in comparison to our study, were restricted to small elemental fingerprints. Fresno et al. [34] were the first to report differences in P, K, Mg, Zn, Fe and Mn, concerning various ripened and unripe Spanish cheeses. Necemer et al. [32] found that, for Slovenian cheeses, the best milk-type indicators were Ca, Br, Zn and Sr. Our study is the first concerning the determination of cheese milk type using elemental metabolomics.

2.2. Contribution to Total Diet—Safety Aspects and Nutritional Value

Repan class="Chemical">garding toxic elements such n>an class="Chemical">as, Cd, Pb, Sn and Sb, the examined samples presented low values for most of them. In more detail, Pb levels were determined to range from 19.8 μg kg−1 (North Aegean) to 38.3 μg kg−1 (Central Greece), i.e., similar to the Khozam et al. study [36] on Lebanese cheese (32.4 μg kg−1 wet weight), but noticeably lower than that reported by Vural et al. study [37] for south-eastern Anatolia-Turkey cheese (4600–7700 μg kg−1 wet weight), the Lante et al. study [38] for Crescenza and Squacquerone cheeses (600 μg kg−1 fresh weight) and the Mendil et al. study [39] for Turkish cheeses (110–960 μg kg−1 wet weight).

pan class="Chemical">Concerning n>an class="Chemical">As levels, they were determined to range from 171 μg kg−1 (Thessaly) to 247 μg kg−1 (North Aegean). Khozam et al. [36] found lower values in Lebanese cheese (2.2 μg kg−1 wet weight). Regarding Cd concentration, it ranged from 3.8 μg kg−1 (Thrace) to 6.2 μg kg−1 (Kriti) in the present study, i.e., lower levels compared to Vural et al. [37] (100 μg kg−1 to 300 μg kg−1 wet weight). Further, Khozam et al. [36] determined Cd in Lebanese cheese to be 0.14 μg kg−1 (wet weight). As regards Sn levels, in the present study they were determined to range from 7.0 μg kg−1 (Makedonia) to 13.4 μg kg−1 (Peloponnisos), while in the Khozam et al. study [36], they were much lower, i.e., 0.037 μg kg−1 (wet weight). Finally, concerning Sb levels, they were measured to range from 4.4 μg kg−1 (Makedonia) to 6.5 μg kg−1 (North Aegean), whereas Khozam et al. [36] determined them to be 0.44 μg kg−1 (wet weight). Regarding other elements, Greek Graviera samples presented lower levels of Thallium, Vanadium, Silver, Antimony, Aluminium and higher levels of Titanium and Barium in comparison to Turkish milk and yogurt [40]. As mentioned above, Hg was not determined.

The Food and Agripan class="Chemical">culture Organization/World Health Organization Joint Expert Committee on Food Additives (JECFA) has established a Provisional Tolerable Weekly Intake (PTWI) for several toxic elements and especially heavy metals. According to the Hellenic Statistical Authority, the daily intake of cheese in the Greek population is 94.7 g/person. For an adult person (e.g., a man of 75 kg body weight), the percentage intake of each toxic element is provided in Table 7. The highest intake is observed for As (12.8%) while the lowest is for Sn (0.0007%), reflecting the high allowed limit for Sn. Column 4 shows the % intake according to the Reasonable Daily Intake of cheese based on the Canadian Food Inspection Agency, 57 g cheese consumption per day. Here, it must be mentioned that reasonable intake has been estimated considering the food habits of Canadians. As regards for these calculations, a mean value of all analyzed Graviera cheese samples was taken into account for each element.

Table 7

Provisional Tolerable Weekly Intake (PTWI) of the studied toxic elements along with the % intake in Greek & Canadian populations.

Element	PTWI (μg kg−1)	% Intake (Greece)	% Intake (Canada)
Al	1000	1.34%	0.81%
As	15	12.8%	7.7%
Cd	7	0.70%	0.42%
Pb	25	0.91%	0.55%
Sn	14,000	0.0007%	0.0004%

pan class="Chemical">Concerning nutrition, trace amounts of n>an class="Chemical">Fe, Mn, Mo, Zn, Co, Ni, Cr, Se, Cu, Si, I, and F are necessary for proper human health, apart from H, C, N, O, Na, K, S, Cl, Mg, Ca, and P which are required in relatively large quantities in a diet. There is also a group of elements called ultra-trace minerals, including V, Sn, Ni, As, and B, that are being investigated for possible biological function but currently do not have clearly defined biochemical roles [19]. Thus, in order to prevent nutrient deficiencies, but also to reduce the risk of chronic diseases such as osteoporosis, cancer and cardiovascular disease, scientific food committees around the world have established specific limits for each element intake with values adapted to different population groups (children, adolescents, pregnant women or older people). The European Commission has established Nutritive Reference Values for adults, according to Regulation (EU) No 1169/2011 (25 October 2011) that are presented in Table 8. WHO/FAO and USDA (United States Department of Agriculture) have also established Recommended Dietary Allowances (RDA) indicating the amount of an individual nutrient that people need for good health depending on their age and gender.

Table 8

Nutritive Reference Values (NRV) for adults, Recommended Daily Allowance (RDA) and Adequate Intake (AI) * for male 31–50 years old, female in pregnancy 31–50 years old and lactation 31–50 years old along with % intake according to Greek and Canadian consumptions.

Element	NRV	% Intake		RDA/AI * Males 31–50 years of age	% Intake		RDA/AI * Pregnancy 31–50 years of age	% Intake		RDA & AI * Lactation 31–50 years of age	% Intake
		Greece	Canada		Greece	Canada		Greece	Canada		Greece	Canada
Ca	800 mg	111	67%	1000 mg	89%	54%	1000 mg	89%	54%	1000 mg	89%	54%
P	700 mg	94%	56%	700 mg	94%	56%	700 mg	94%	56%	700 mg	94%	56%
Mg	375 mg	12%	7.3%	420 mg	11%	6.5%	360 mg	13%	7.6%	320 mg	14%	8.5%
Fe	14 mg	20%	12%	8 mg	35%	21%	27 mg	10%	6.2%	9 mg	31%	19%
Zn	10 mg	31%	19%	11 mg	29%	17%	11 mg	29%	17%	12 mg	26%	16%
Cu	1 mg	6.8%	4.1%	900 μg	7.6%	4.6%	1000 μg	6.8%	4.1%	1300 μg	5.2%	3.2%
Mn	2 mg	1.8%	1.1%	2.3 * mg	1.5%	0.9%	2.0 * mg	1.8%	1.1%	2.6 * mg	1.4%	0.8%
Se	55 μg	13%	7.8%	55 μg	13%	7.8%	60 μg	12%	7.2%	70 μg	10%	6.1%
Cr	40 μg	139%	83%	35 * μg	158%	95%	30 * μg	185%	111%	45 * μg	123%	74%
Mo	50 μg	20%	12%	45 μg	22%	13%	50 μg	20%	12%	50 μg	20%	12%

pan class="Chemical">As shown in Table 8, % Ca intake from Graviera wn>an class="Chemical">as sufficient and ranged from 54% to 111%, while P ranged from 56% to 94%. Moreover, the Ca–P ratio was 1.4:1, so consumption of Graviera cheese is one the most convenient ways for proper intake of both minerals through the diet. High dietary Ca–P ratios play important role in bone health [41]. The % zinc intake ranged from 16% to 31%. Iron, Cr and Mo % intakes ranged from 6.2% to 35%, 74% to 185% and 12% to 22%, respectively. Regarding Mg, Cu, Mn and Se the % intake is less significant. These results highlight Graviera cheese as good source of trace and macro elements, especially for Ca, P, Zn, Cr, Fe and Mo.

pan class="Chemical">Compn>ared with other studies like the Moreno-Rojn>an class="Chemical">as et al. study of different cheese types [10], nutritive elements were found at similar levels. Camin et al. [12] found lower Fe and higher Se and Mo in various European hard cheeses such as PDO Parmigiano Reggiano. Suhaj et al. [33] determined Cr in lower levels, Mo and Ca in slightly lower levels and Mn in slightly higher levels in some European Emmental and Edam hard cheeses than the Graviera samples in our study.

3. Materials and Methods

3.1. Instrumentation and Reagents

Chemicpan class="Chemical">als used were nitric acid (Suprapur®, 65% w/v, Merck, Darmstadt, Germany), hydrogen peroxide (Suprapur®, 30% w/v, Merck, Darmstadt, Germany), ICP internal standards of Ge and In and ICP-MS certified multi-element standards (all from Inorganic Ventures, NJ, USA). Ultrapure water with a resistance of 18.2 MΩ cm−1 obtained from a MilliQ plus system (Millipore, Saint Quentin Yvelines, France) was used in all procedures.

Elementpan class="Chemical">al n>an class="Chemical">content was determined using a Perkin Elmer (SCIEX, Toronto, ON, Canada) 9000 Series ICP-MS. Inductively coupled plasma mass spectroscopy is predominantly used in authentication studies due to its capability for rapid ultra-trace level multi-element determinations [42].

3.2. Sample Collection, Preparation and Digestion

One hundred and five Graviera cheepan class="Chemical">se sampn>les were upan class="Chemical">sed for the purposes of this study. The geographical origin of the samples is reported in Table S1 and depicted in the map shown in Figure 1. Most of the samples were collected from small-medium dairies and the rest from the respective local markets. The sampling strategy excluded large dairies that are able to collect milk from different Greek regions and bulk it in their premises for the production of their own trade mark. Samples were taken from a lot and after grinding, they were preserved in a freezer (−32 °C) before analysis.

Figure 1

Map depicting origin of samples.

Sample digestion wpan class="Chemical">as performed with a mipan class="Chemical">crowave-assisted digestion system (CEM, Mars X-Press, Matthews, NC, USA). Approximately 0.50 g of cheese was weighted in an analytical balance in a polypropylene tube. Then, 4.0 mL of HNO3 was added to pre-digest samples for 30 min. The resulting cheese suspension was transferred quantitatively, with the use of 4.0 mL HNO3 and 2.0 mL H2O2 to the microwave digestion PTFE vessel. The samples were heated in the microwave accelerated digestion system according the following program: the power was ramped during 20 min from 100 to 1200 W and held for 15 min. The temperature reached a maximum of 200 °C and followed by a cool-down cycle for 15 min. PTFE vessels were sealed throughout the aforementioned cycle to avoid volatilization losses. Although all samples were completely brought to solution, to disregard any small particle passing optical inspection entering the ICP-MS, solutions were filtered with polyester disposable syringe filters 0.20 μm/ 15 mm (Chromafil, Macherey-Nagel, Düren, Germany). Before injection in the ICP-MS, sample solutions were diluted, as required, with ultrapure water.

3.3. ICP-MS Analysis

The studied elements pan class="Chemical">asn>an class="Chemical">sessed were:

pan class="Chemical">REEs: Ce, Dy, Er, Eu, Gd, Ho, La, Lu, Nd, Pr, Sc, Sm, Tb, Tm, Y, Yb

pan class="Chemical">Actinides: Th, U

Precious pan class="Chemical">metals: n>an class="Chemical">Au, Ir, Pd, Pt, Re, Rh, Ru

Ultra-trace elements: Hf, Nb, pan class="Chemical">Sb, n>an class="Chemical">Sn, Ta, Tl, W, Zr

Trace elements: Ag, pan class="Chemical">Al, n>an class="Chemical">As, B, Ba, Bi, Cd, Co, Cr, Cs, Cu, Fe, Ga, Mn, Mo, Ni, Pb, Rb, Se, Si, Sr, Ti, V, Zn

Mapan class="Chemical">cro elements: Ca, K, n>an class="Chemical">Mg, P

Limits of quantification for pan class="Chemical">all were lower than thon>an class="Chemical">se determined in the samples (Table S3). Operating conditions of the ICP-MS were as follows: nebulizer gas flow of 0.75 L min−1, ICP RF power of 950 W, lens voltage of 7 V, pulse stage voltage of 950 V and sample uptake rate of 26 rpm. Calibration curves ranges were from 1 ng kg−1 to 1000 μg kg−1 for rare earths, precious metals and ultra-trace elements, while from 0.01 μg kg−1 to 10 mg kg−1 for trace and macro elements. Indium was used as internal standard for rare earths, precious and ultra-trace elements, while germanium was used for trace and macro elements. In detail, the daily analytical procedure is:

Start

Daily perforpan class="Species">mance check with n>an class="Chemical">set up solution from Perkin Elmer that contains 10 μg L−1 of the following elements: Be, Mg, Co, In, Pb

Standard solutions for pan class="Chemical">REEs, n>an class="Chemical">actinides, precious metals and ultra-trace elements of 0.001, 0.01, 0.1, 0.5, 1, 5, 10, 50, 100 and 1000 μg L−1

Standard solutions for mapan class="Chemical">cro and trace elements of 0.01, 0.1, 0.5, 1, 5, 10, 50, 100, 1000 and 10000 μg L−1

Steps to pan class="Chemical">be repn>eated after 4 h:

Standard repan class="Chemical">ference materin>an class="Chemical">als

Blank

Samples

End

Standard repan class="Chemical">ference materin>an class="Chemical">als

3.4. Calibration and Quality Assurance

To pan class="Chemical">asn>an class="Chemical">sess the accuracy of the process the following standard reference materials were obtained from the European Commission, Joint Research Center, institute for reference materials and measurements IRMM, Belgium and the National Institute of Standards & Technology (NIST), USA:

Trace and mapan class="Chemical">cro elements, ERM-BD151 skimmed n>an class="Gene">milk powder (IRMM), RM 8414 bovine muscle powder (NIST) and RM 1573a tomato leaves (NIST).

Rare earth elements and pan class="Chemical">actinides n>an class="Chemical">CRM-668 mussel tissue (IRMM).

The standard repan class="Chemical">ference materin>an class="Chemical">als were subjected the same analytical process: Digested three different times, each digestate measured in triplicate (Table 9). Recoveries were in the range 67–121% for all elements other than Se. In order to overcome Ar2+ interferences we measured Se 82.

Table 9

Standard reference material results, % recoveries and % RSD, n = 9 (three different digestions, measured in triplicate).

	BCR 668					ERM-BD151
REEs	Certified μg kg−1	Found μg kg−1	% Recoveries	% RSD	TREs	Certified mg kg−1	Found mg kg−1	% Recoveries	%RSD
Ce	89 ± 7	82 ± 3	92	4	Cd	0.106 ± 0.013	0.121 ± 0.041	114	34
Dy	8.9 ± 0.6	7.9 ± 0.9	89	11	Cu	5.00 ± 0.23	5.77 ± 0.97	115	17
Er	4.5 ± 0.5	3.7 ± 0.7	82	19	Fe	53 ± 4	50 ± 13	94	26
Eu	2.79 ± 0.16	2.4 ± 0.10	86	4	Mn	0.29 ± 0.03	0.34 ± 0.09	117	26
Gd	13.0 ± 0.6	12.1 ± 0.7	93	6	Pb	0.207 ± 0.014	0.201 ± 0.049	97	24
Ho	1.8 ± 0.6 a	1.2 ± 0.2	67	17	Se	0.19 ± 0.04	0.27 ± 0.09	142	33
La	80 ± 6	71 ± 5	89	7	Zn	44.9 ± 2.3	54.0 ± 7.9	120	15
Lu	0.389 ± 0.024	0.379 ± 0.010	97	3	Macro elements
Nd	54 ± 4	49 ± 1	91	2	Ca	13.9 ± 0.7	9.9 ± 1.4	71	14
Pr	12.3 ± 1.1	11.1 ± 0.2	90	2	K	17.0 ± 0.8	14.0 ± 1.7	82	12
Sc	8.5 ± 1.8	10.3 ± 1.5	121	14	Mg	1.26 ± 0.07	1.20 ± 0.13	95	11
Sm	11.2 ± 0.8	10.4 ± 0.4	93	4	Na	4.19 ± 0.23	3.89 ± 0.7	93	18
Tb	1.62 ± 0.12	1.60 ± 0.07	99	4	P	11.0 ± 0.6	12.3 ± 2.4	112	20
Tm	0.48 ± 0.08	0.571 ± 0.03	119	5
Y	59 ± 5	53 ± 3	90	6
Yb	2.8 ± 0.5 a	2.1 ± 0.1	75	5
Actinides
Th	10.7 ± 1.2	9.7 ± 0.3	91	3
U	56 ± 5	51 ± 4	91	8

a Indicative values.

3.5. Statistical Analysis

Statisticpan class="Chemical">al anpan class="Chemical">alysis was performed using SPSS software (IBM, Armonk, NY, USA) for the descriptive statistics and cross validation and Statgraphics Centurion XV software (Statpoint technologies, Warrenton, VA, USA) in order to analyze the data using statistical models and predictive analyses (Linear Discriminant Analysis).

4. Conclusions

We prepan class="Chemical">sent results from 61 elements in cheen>an class="Chemical">se for the first time with implications in food authentication, safety and nutrition. Further work is in progress for data selection, increasing confidence for food authentication using bioinformatics tools under development. This is the first study reporting signatures of 61 elements including rare earth elements and all the precious metals in cheese. We highlight the application of elemental metabolomics to human nutrition assessing both nutritive and toxic elements. The results demonstrate that elemental metabolomics could be potentially used for discrimination of cheeses produced in different geographical zones and milk type. The method needs further improvement by bioinformatics tools to automate data cleaning done manually.

In pan class="Chemical">compn>arison to molepan class="Chemical">cular analysis, elemental metabolomics is simple and convenient. The first step is accurately weighing samples in capped polypropylene tubes to analyze when convenient, when adequate samples are collected and when instrumentation is available. There are no requirements concerning temperature, time, or any other storage condition. The only requirement is creation of comprehensive elemental metabolome databases for food authentication, quality and safety. Elemental metabolomics are becoming more affordable by lowering the ICP-MS purchasing cost and increasing capabilities concerning interferences [19]. We envisage open access elemental databases for improvement of human nutrition and health.