Multielement Characterization and Antioxidant Activity of Italian Extra-Virgin Olive Oils

Food product safety and quality are closely related to the elemental composition of food. This study combined multielement analysis and chemometric tools to characterize 237 extra-virgin olive oil (EVOO) samples from 15 regions of Italy, and to verify the possibility of discriminating them according to different quality factors, such as varietal or geographical origin or whether they were organically or traditionally produced. Some elements have antioxidant properties, while others are toxic to humans or can promote oxidative degradation of EVOO samples. In particular, the antioxidant activity of oils’ hydrophilic fraction was estimated and the concentrations of 45 elements were determined by inductively coupled plasma mass spectrometry (ICP-MS). At first, univariate and multivariate analyses of variance were used to compare the element concentrations, and statistically significant differences were found among samples from different regions. Successively, discriminant classification approaches were used to build a model for EVOO authentication, considering, in turn, various possible categorizations. The results have indicated that chemometric methods coupled with ICP-MS have the potential to discriminate and characterize the different types of EVOO, and to provide “typical” elemental fingerprints of the various categories of samples.

1 Introduction

The elemental composition of foods is of toxicological and nutritional interest and can be considered an important quality parameter (Astolfi et al., 2021a; 2020a; 2020b). In particular, the concentrations of trace elements in extra-virgin olive oil (EVOO) are also one of the criteria for the assessment of the quality regarding storable period and freshness (Choe and Min, 2006). In fact, some elements, such as Ca, Co, Cu, Fe, Mg, Mn, Ni, and Sn, can promote the oxidative degradation of this important component of the Mediterranean diet appreciated among consumers for its nutritional properties and specific flavor (Choe and Min, 2006). Other elements (such as As, Cd, Cr, Cu, Hg, and Pb) present in EVOO are of great concern because they are toxic and potentially carcinogenic to humans even at low concentration (Tchounwou et al., 2012). The International Olive Council has established, as a quality criterion, a maximum residue level (MRL) for the content of As, Cu, Pb (0.1 mg kg⁻¹), and Fe (3 mg kg⁻¹) in olive oils and olive–pomace oils (International Olive Council, 2019), and the maximum levels of Cu and Fe in other vegetable oils have been also legislated (Codex Stan 33-1981, 2021), varying from 0.1 up to 5.0 mg kg⁻¹. Recently, element determination in EVOO samples has gained importance for oil geographical traceability and authentication (Cordella et al., 2002; Dugo et al., 2004; Benincasa et al., 2007; Cabrera Vique et al., 2012; Camin et al., 2010; Beltran et al., 2015; Bajoub et al., 2018; Aceto et al., 2019; Damak et al., 2019; Zaroual et al., 2021). In particular, elements are useful in the characterization of protected designations of origin (PDOs) or protected geographical indications (PGIs) (European Union (EU), 2012), and they can also contribute to determine EVOO geographical provenance of non-PDO oils (Beltran et al., 2015; Aceto et al., 2019). In fact, the presence of metals in EVOO varies according to their origin and can be due to natural contamination from the soil, environment, fertilizers, and genotype of the plant or to the production process and contact with storage materials (Zeiner et al., 2005; Chatzistathis et al., 2009; Kabata-Pendias, 2010; Lepri et al., 2011; Yaşar et al., 2012; Bakircioglu et al., 2013). A suitable statistical treatment of trace element data could allow a geographical characterization of different EVOO samples. Principal component analysis (PCA) and hierarchical cluster analysis (HCA) (Gumus et al., 2017; Luka and Akun, 2019; Russo et al., 2020; Savio et al., 2014), linear discriminant analysis (LDA) (Benincasa et al., 2007; Cabrera-Vique et al., 2012; Beltran et al., 2015; Damak et al., 2019), classification trees (CTs) (Gumus et al., 2017), and artificial neural networks (ANNs) (Farmaki et al., 2012; Gonzalez-Fernandez et al., 2019) have been used most.

Several beneficial implications of EVOO are derived from its antioxidant content (Dugo et al., 2020; Hannachi and Elfalleh, 2020). Intake of antioxidant compounds from oil, such as phenols, phenolic acids, and flavonoids (Capriotti et al., 2014), is usually related to health well-being. As well known, natural antioxidants play a key role in contrasting reactive species activity in living organisms, thus preventing oxidative stress-related diseases, such as cardiovascular and neurodegenerative illness and many other chronic disorders (Pérez-Jiménez et al., 2008; Cioffi et al., 2010; Šarolić et al., 2014). Moreover, antioxidants prevent lipid oxidations that cause quality degradation and unpleasant taste formation in edible oils (Christodouleas et al., 2015). Therefore, estimation of antioxidant capacity is crucial for evaluating oil’s healthy and organoleptic properties. One of the most widely used in vitro procedures to routinely and globally estimate oil antioxidant power is the 2,2-diphenyl-1-picrylhydrazyl spectrophotometric assay (DPPH) that has the possibility of being easily applied to a high number of samples, allowing a great level of reliability (Kedare and Singh, 2011; Frezzini et al., 2019). The assay is based on the quantitative measurement of the decrease of absorbance due to the scavenging capacity of antioxidants present in the sample toward DPPH free radicals (Christodouleas et al., 2015).

All the described aspects making trace element determination, as well as the antioxidant activity of EVOO samples, are very important for both economic and health contexts (Zaroual et al., 2021; Bajoub et al., 2018). In particular, the European Union is the first producer, consumer, and exporter of olive oil in the world (Eurostat, 2019; IOC, 2018a,b). Italy follows Spain, the first world producer with an average of 20% of the total European olive oil production. About two-thirds of total Italian production is represented by EVOO (Carbone et al., 2018). Therefore, the use of a rapid and accurate analytical method for trace element analysis in EVOO has a great importance in quality control and food analysis (Llorent-Martínez et al., 2011; Shah and Soylak, 2021). Unfortunately, the determination of trace elements in EVOO samples is particularly difficult to perform, as some of them are present at very low concentrations and due to high complexity of the matrix (Shah and Soylak, 2021; Trindade et al., 2015). Sample preparation of EVOO samples is a critical step, and the determination of trace elements in EVOO requires very sensitive instrumental techniques such as inductively coupled plasma–mass spectrometry (ICP-MS) (Astolfi et al., 2021b).

The main purpose of this study is to evaluate the most significant relationships between element levels in EVOO and different categorizations, mostly related to the geographical origin using chemometric tools coupled with the ICP-MS method. For this purpose, 45 elements from a total of 237 EVOO samples from 15 Italian regions were analyzed. Also, the antioxidant activity of oils’ hydrophilic fraction (HF) was estimated by the DPPH assay. The corresponding data set constituted the basis for building and validating classification models for the discrimination of the samples according to specific categorizations, which reflect possible quality attributes of the oils (and for which there could be a statistically significant number of individuals available). In particular, discriminant classification models were built using partial least square discriminant analysis (PLS-DA) to account for the possibility of dealing with correlated variables and low samples to variable ratios; moreover, to evaluate model stability and, at the same time, their reliability in an unbiased way, also in cases where the available number of samples per category was not too large, a repeated double cross-validation strategy (rDCV) was adopted.

2 Materials and methods

Sample collection

EVOO samples (N = 237) were collected between 2017 and 2018 from 15 production regions of Italy and different cultivars. In particular, a total of 64 EVOO samples were with PDOs and 21 with PGIs (European Union (EU), 2012). Table 1 summarizes the number of EVOO samples according to their geographical provenances in terms of the regions. All samples (∼100 mg) were kept in screw-capped glass vials in the dark at room temperature until analysis.

TABLE 1

TABLE 1. Number of samples of extra virgin olive oil for each considered category.

Chemicals

All the solutions were prepared with deionized water (18.3 MΩ cm resistivity) obtained from an Arioso (Human Corporation, Seoul, Korea) Power I RO-UP Scholar UV deionizer system. HNO₃ at 67% (suprapure; Carlo Erba Reagents, Milan, Italy), H₂O₂ at 30% (suprapure; Merck KgaA, Darmstadt, Germany), and Ar, He, and H₂ gases at 99.9995% (SOL Spa, Monza, Italy) were used.

For ICP-MS analysis, all calibration standard solutions were prepared from a 1,000-mg l⁻¹ multielement standard solution (VWR International, Milan, Italy) by dilution with 10% (v/v) HNO₃ and H₂O₂ (2:1 v/v). Single standard solutions of In, Rh, Sc, and Th (at 0.010 mg l⁻¹ from 1,000 ± 5 mg l⁻¹; Merck KGaA, Darmstadt, Germany) and Y (at 0.005 from 1,000 ± 2 mg l⁻¹; Panreac Química, Barcelona, Spain) were used as internal standards. A multielemental solution containing Ba, Be, Ce, Co, In, Pb, Mg, Tl, and Th (at 0.005 mg l⁻¹ from 10.00 ± 0.05 mg l⁻¹; Spectro Pure, Ricca Chemical Company, Arlington, TX, USA) was used to check the instrument performance.

For the estimation of the antioxidant activity of EVOO samples, DPPH was purchased from Sigma Aldrich Co. (St. Louis, MO, USA).

Sample preparation and analysis

2.1.1 Analysis of elements

Duplicate samples (∼0.5 g) of each EVOO variety were accurately weighed in 10-ml disposable graduated tubes (Artiglass, Due Carrare, PD, Italy). Then, 5 ml reagent mixture of 10% (v/v) HNO₃ and H₂O₂ (2:1 v/v) was added to each tube and heated in a water bath (WB12, Argo Lab, Modena, Italy) at 95°C for 40 min (Astolfi et al., 2021b). The lower aqueous phase was transferred into a clean tube and subjected to the ICP-MS (820-MS; Bruker, Bremen, Germany) analysis without further dilutions. The elements were monitored in standard and collision–reaction interface (CRI) modes to check and reduce possible polyatomic interference, and the following isotopes were used: ⁷Li, ⁹Be, ¹¹B, ²³Na, ²⁴Mg, ²⁷Al, ²⁸Si, ³¹P, ³⁹K, ⁴⁴Ca, ⁴⁹Ti, ⁵¹V, ⁵²Cr, ⁵⁵Mn, ⁵⁷Fe, ⁵⁹Co, ⁶⁰Ni, ⁶⁵Cu, ⁶⁶Zn, ⁷¹Ga, ⁷⁵As, ⁷⁸Se, ⁸⁵Rb, ⁸⁸Sr, ⁹⁰Zr, ⁹³Nb, ⁹⁸Mo, ¹⁰⁷Ag, ¹¹²Cd, ¹¹⁸Sn, ¹²¹Sb, ¹²⁵Te, ¹³³Cs, ¹³⁷Ba, ¹³⁹La, ¹⁴⁰Ce, ¹⁴¹Pr, ¹⁴⁶Nd, ¹⁵⁹Tb, ¹⁶³Dy, ¹⁸²W, ²⁰⁵Tl, ²⁰⁸Pb, ²⁰⁹Bi, and ²³⁸U. CRI was used with He (30 ml min⁻¹) and H₂ (70 ml min⁻¹) as cell gases. The ICP-MS operating conditions and parameters were as follows: radiofrequency power 1,400 W; plasma Ar flow rate 18 l min⁻¹; auxiliary Ar flow rate 1.8 l min⁻¹; nebulizer gas flow rate 0.9 l min⁻¹; peak hopping scanning mode; steady-state analysis mode; dwell time between 50 and 100 ms, pump rate 3 rpm; five scans/replicate; and three replicates/sample. For the quantitative analysis of EVOO samples, calibration curves were built on seven different concentrations between 0.00025 and 0.05 mg l⁻¹ and 0.0125 and 5 mg l⁻¹ for all trace and major elements, respectively.

2.1.2 Estimation of antioxidant activity

DPPH assay was performed according to the procedure described by Šarolić et al. (2014) with slight modifications. In detail, ∼0.5 g of each EVOO sample was mixed with 1 ml of 80:20 (v/v) CH₃OH:H₂O, and the mixture was blended in an ultrasonic bath (PROCLEAN 10.0 ultrasonic cleaner; Ulsonix, Berlin, Germany) for 15 min at 30°C. When the two phases appeared, the hydrophilic phase was collected, and the extraction was repeated another two times. Then, the hydrophilic extracts were combined to get a homogeneous sample. To perform DPPH assay, 50 µl of the HF sample was added to 2 ml of methanolic DPPH (0.04 mM), then the mixture was shaken for 30 min by rotating agitation (60 rpm; rotator; Glas-Col, Terre Haute, IN, USA) at room temperature in the dark and analyzed by UV-Vis spectrophotometry (Varian Cary 50 Bio UV-Vis; Varian Inc., Palo Alto, CA, USA) set at 517 nm, by measuring the sample absorbance decrease against the control (blank solution). Solutions were prepared daily and used fresh, and three replicates of each type of oil were performed. The DPPH radical scavenging activity was calculated in terms of percentage reduction of DPPH according to the following equation:

$\text{DPPH }\left[\text{\%}\right]=\frac{\left({\text{A}}_0-{\text{A}}_{\text{S}}\right)}{{\text{A}}_0}\times 100$

where A₀ represents the absorbance of the blank solution and A_S is the absorbance of the sample.

Quality assurance

The method accuracy for element determination was checked by recovery assays in the EVOO samples adding element at the low (0.005 and 0.02 mg l⁻¹) and high (0.2 and 1 mg l⁻¹) spike concentrations for all trace and major elements (B, Ca, K, Mg, Na, P, Si, and Sr) and always in the linear calibration range. In addition, accuracy was tested by a certified reference material (Conostan S-21; lot number: 21550100) obtained from SCP SCIENCE (Baie D’Urfé, Canada). The recoveries fell within 20% of the expected value and reproducibility lower than 20% (Astolfi et al., 2021b). The method detection and quantification limits (MDL and MQL, respectively) were in the range 0.004–510 and 2.5–5,000 μg kg⁻¹, respectively. Only the Ca, Cr, Mg, Mn, Ni, P, Rb, Ti, and Zn levels in the EVOO samples were 100% greater than the MDL. The possible instrumental drift for the ICP-MS analysis was checked and corrected using an internal standard solution of In, Rh, Sc, Th, and Y (Astolfi et al., 2021b; 2020c). Blank samples and control standards were tested every 20 samples in each run, and recalibration was performed every 100 samples.

Statistical analysis

The data were statistically evaluated according to the procedures of the software SPSS Statistics 25 (IBM Corp., Armonk, NY, USA) for univariate analysis. Analytical replicates were averaged prior to the successive elaboration. Non-parametric tests (Kruskal–Wallis and pairwise post-hoc) were applied because of the unequal numbers of samples per group and the not normal distribution (Soliani, 2003). The element concentrations measured below MDL were substituted by its half value (MDL/2) for the statistical elaboration (Farmaki et al., 2012). A p-value lower than 0.05 was considered statistically significant.

Partial least square discriminant analysis (PLS-DA; Ståhle and Wold, 1987; Barker and Rayens, 2003) implemented through in-house written functions running under the Matlab environment (R2015b, v.8.6, The MathWorks Inc., Natick, MA, USA) was used to build multivariate classification models. PLS-DA is a regression-based classification model which operates by coding class belonging by means of a dummy binary response matrix (or vector, when the problems involve only pairs of classes, as in the present study). In particular, if discrimination is sought between two categories, class belonging of the training samples is described by the vector y, having 1 in correspondence of all the individuals from the first class and 0 in all the remaining positions (i.e., those corresponding to the second group). A PLS model (Wold et al., 1983) is then built between the experimental data X and the dummy vector y, and the predicted value of the response ($\widehat{\mathit{y}})$ constitutes the basis for the classification of the samples: since the predicted responses are real-valued, an optimal threshold $y_{thres}$ has to be calculated so that, if the predicted response is greater than $y_{thres}$, the sample is predicted as class 1, otherwise as class 2. In the present study, the threshold was calculated by applying LDA on the predicted responses calculated on the training samples (Perez et al., 2009).

The reliability of the classification models was evaluated by means of a repeated double-cross-validation (rDCV) procedure (Filzmoser et al., 2009). Double cross-validation (DCV) is a validation strategy which involves two nested loops of cross-validation: an inner loop for model selection (i.e., for choosing the optimal number of latent variables) and an outer loop which mimics an external (i.e., not involved in any model building and/or optimization stage) test set, to be used for estimating the prediction and generalization ability. In order to avoid that the performances of the model depend on a particular sample splitting scheme, the procedure is repeated a sufficient number of times, changing the distribution of the individuals across the different cancelation groups, hence the name “repeated” DCV. Repeating the double-cross-validation procedure allows also having multiple predictions for the same samples, which translates to the possibility of estimating confidence intervals for all the classification figures of merit and model parameters.

3 Results and discussion

Levels of elements

Table 2 shows the concentration of the elements in EVOO from all over Italy. The content of Tl was below the respective MDL (0.06 µg kg⁻¹) in all the samples. Si, Te, and Nb were found above the MDL (270, 0.03, and 0.04 µg kg⁻¹, respectively) only in 1%, 3%, and 7% of all samples. Only the Ca, Cr, Mg, Mn, Ni, P, Ti, and Zn levels in the EVOO samples were 100% greater than the MDL. The maximum concentrations for As (4.0 μg kg⁻¹), Cu (41.6 μg kg⁻¹), Fe (582 μg kg⁻¹), and Pb (22.1 μg kg⁻¹) were lower than the MRLs established by the IOC for olive and pomace-olive oils, which are 100 μg kg⁻¹ for As, Cu, and Pb and 3,000 μg kg⁻¹ for Fe (International Olive Council, 2009). Calcium showed the highest concentration ranging from 1,230 to 35,700 μg kg⁻¹, whereas from 10- to 50-fold lower levels were found for Fe, Mg, Na, P, and Zn (median = 77, 91, 110, 272, and 111 μg kg⁻¹, respectively).

TABLE 2

TABLE 2. Method detection limits (MDL; μg kg⁻¹) and element levels [median, minimum (min) and maximum (max); μg kg⁻¹] in extra-virgin olive oils (EVOO; n = 237) from all over Italy.

Concentrations of elements obtained in this study were compared to levels measured in EVOO from several other Mediterranean countries (Supplementary Tables S1–S4). Levels of many elements showed wide variability even within the same country. The Ag, Ba, P, and Sn data were not considered because these elements are not completely extracted with the method used. As regards the content of B, Be, Dy, Nd, Pr, Si, Tb, and Te, we could not find other data for EVOO in the literature. Our results were similar to those reported by another study on Italian EVOO (Benincasa et al., 2007); on the contrary, they differed significantly from other data concerning most of the elements investigated in the EVOOs of Spain (Beltran et al., 2015; Llorent-Martínez et al., 2014), Croatia (Pošćić et al., 2019), Tunisia (Damak et al., 2019), and Turkey (Gumus et al., 2017). The concentrations of Ca (1,230–35,700 μg kg⁻¹), Cr (0.4–839 μg kg⁻¹), Mg (21–723 μg kg⁻¹), and Ni (2.1–49.7 μg kg⁻¹) found in this study were in the same range to that found in other Italian EVOO (Ca = 1850–26,900 μg kg⁻¹; Cr = 116–437 μg kg⁻¹; Mg = 56–1,030 μg kg⁻¹; and Ni = nd-46.9 μg kg⁻¹) as reported by Benincasa et al. (2007), but from 10 to 100 times higher than the levels reported in Croatian (Pošćić et al., 2019) and Turkish EVOO (Gumus et al., 2017). Fe concentrations (<12–582 μg kg⁻¹) varied from 100 times lower to 100 times higher than the level of Fe quantified in EVOO from Turkey (1–14,670 μg kg⁻¹) by Gumus et al. (2017) and Croatia (0.19–2.57 μg kg⁻¹) by Pošćić et al. (2019) or Spain (0.5–1.2 μg kg⁻¹) by Beltran et al. (2015), respectively. This variability in the concentrations of the elements present in EVOO samples may depend on various factors related to the geochemistry of the provenance soil but also to physiological aspects typical of the species from which a particular EVOO derives (Giaccio and Vicentini, 2008).

Grouping the data according to geographic origin as north (Emilia Romagna, Liguria, Lombardy, Trentino Alto Adige, and Veneto), center (Abruzzo, Lazio, Marche, Tuscany and Umbria), and south (Apulia, Calabria, Campania, Sardinia and Sicily) of Italy, it is possible to identify elements that differ significantly from one group to another (Table 3). In particular, the EVOO samples from northern Italy had significantly higher levels of Cs, Fe, Na, P, and Pr than those from central Italy and Fe, Pr, and U than those from southern Italy. Both Fe and Pr appear to provide a good tool for tracing the EVOO production chain in accord with other authors (Aceto et al., 2019; Damak et al., 2019). Iron is common in silicates and carbonates present in soil (Pohl, 2011); however, some authors reported that Fe may be present in edible oils as a result of storage and processing contaminations (Mendil et al., 2009; Zeiner et al., 2010). Praseodymium and the other lanthanides do not have a defined role in the metabolism of plants; therefore, their distribution remains almost unchanged in the passage from the soil to the fruits (Aceto et al., 2019). For this reason, these elements can be used as fingerprints to discriminate the geographic origin of the EVOO samples (Farmaki et al., 2012; Aceto et al., 2019). In addition, the analysis of some elements in EVOO, such as Cs and Rb, which can be easily mobilized in the soil, can be linked to a geogenic source rather than an anthropogenic origin (such as extraction process or cultivation practices) and can help in the geographical traceability of EVOO samples (Kelly, Heaton, & Hoogewerff, 2005).

TABLE 3

TABLE 3. Element levels [median, minimum (min) and maximum (max); μg kg⁻¹] in extra-virgin olive oils from north (n = 20), central (n = 132) and south (n = 85) Italy.

By comparing the concentrations of the elements in the EVOO samples from each region (Supplementary Tables S5–S7), the number of elements that differ significantly increases. Table 4 shows a summary of all the elements that differ significantly according to the region. Emilia Romagna was not considered for the comparison because there was only one EVOO sample to consider. EVOOs from Lombardy did not have levels of elements that are significantly different from those of oils from all other regions. Considering the other oils of northern Italy, the EVOOs from Trentino and Liguria differed significantly from the EVOOs from Marche only for the content of Na, which in the EVOOs from Marche (median = 38 μg kg⁻¹) was about four times lower, while the EVOOs from Veneto had a higher content of Fe (median = 218 μg kg⁻¹) than the oils from Abruzzo (median = 15 μg kg⁻¹) and a higher content of Fe and Na (median = 218 and 174 μg kg⁻¹, respectively) compared to the Marche. The EVOO samples from Lazio differed significantly for a large number of elements (Ba, Ca, Cd, Ce, Cs, Dy, Ga, La, Mg, Na, Nd, Pr, Pb, Rb, Sb, Sr, Tb, Ti, U) compared to Tuscany, Abruzzo, Campania, and Marche. In all cases, levels of Cd (median = 0.14 μg kg⁻¹), La (median = 0.20 μg kg⁻¹), and Rb (median = 0.48 μg kg⁻¹) were higher than those of oils from other regions mentioned above.

TABLE 4

TABLE 4. Summary of significant differences within medians of the 45 selected elements and antioxidant activity (DPPH%) among all samples from Italian regions by Kruskall–Wallis and pairwise post-hoc tests. A p-value lower than 0.05 was considered statistically significant.

Antioxidant activity

Following the extraction and storage of EVOO, it is inevitable that an oxidation process occurs, which leads to a deterioration of the oil (Bendini et al., 2007). Some factors such as temperature, light, oxygen and other chemical elements, unsaturated fatty acid composition, and the presence of antioxidants can affect the oxidation process differently (Frankel 1985). Phenolic compounds have antioxidant capacities in EVOO since they can eliminate peroxyl and alkoxy radicals and chelate transition metal ions present in traces (Visioli et al., 1998). Several elements are known for their antioxidant properties (Perna et al., 2012; Thiruvengadam et al., 2020). Indeed, in the present study (Supplementary Table S8), significant correlations were observed between the antioxidant activity and elements. A positive and significant moderate correlation (r = 0.500–0.768, p = 0.05) was observed between Al, Ca, Fe, V, and Zr in EVOOs from Abruzzo, Ba in EVOOs from Apulia, and B, Mn, Se, and V in EVOOs from Sardinia and the DPPH% data. Conversely, a low and positive correlation (r < 0.4) was recorded between Ba and Ni and the antioxidant activity of all samples. Other elements (Ag, Cr, Cu, Li, Sb, Si, and Tl) might not affect the antioxidant properties as non-significant correlations were observed between them.

Supplementary Tables S5–S7 show the antioxidant activity measured by the DPPH assay (DPPH%) in the EVOO samples from each region. EVOOs from central and southern Italy showed higher antioxidant activity than oils from northern Italy. In particular, Table 4 shows that EVOOs from Sicily had a significantly lower DPPH% (median = 18.2%) than oils from Abruzzo (median = 47%), Apulia (median = 37.6%), and Tuscany (median = 36.2%), while the EVOOs from Liguria had significantly lower DPPH% (median = 15.7%) compared to Tuscany. The highest data of DPPH% (67.3%) was found in the oils of Campania. Cioffi et al. (2010) demonstrated that oils from Campania have antioxidant properties, which are very likely due to the presence of high contents of phenolic compounds.

Classification of EVOOs according to geographical origin

At first, the possibility of discriminating the different EVOOs according to their geographical origin was considered. In particular, due to the unavailability of the information about the origin of all the samples and to the unbalancedness in the distribution of samples per class, when considering the oils of known origin, several two-class models (i.e., comparing two regions at a time) were built and validated. Here it must be further stressed that all the regions for which the available number of certified individuals was too low to be considered representative have not been included in the comparison.

In all cases, PLS-DA analysis was carried out on the matrix made up of the concentrations of the elements presenting at least 70% of the values above the limit of detection (so to avoid possible artifacts related to data imputation) and including also TEAC and DPPH. Models were built after autoscaling and validated by means of an rDCV procedure with 50 runs, 10 cancelation groups in the outer loop (the one mimicking the external test set) and 5 in the inner loop (the one used for model selection, i.e., definition of the optimal number of latent variables). The results obtained are summarized in Table 5, where the accuracy, the mean correct classification rate, and the sensitivities for the two compared classes are reported. Since two-class discriminant models were calculated, due to symmetry the sensitivity of a class (true positive rate) is the specificity (true negative rate) of the other category; this is why sensitivities only have been reported. Moreover, since the number of samples per class was, in some cases, highly unbalanced (Table 1), we have decided to report both classification accuracy (percentage of correctly classified samples irrespectively of the category over the total number of samples) and the mean correct classification rate, which is the average of the specificities of the two classes.

TABLE 5

TABLE 5. PLS-DA discrimination between pairs of geographical origin. Figures of merit estimated on the outer loop of the rDCV procedure (expressed as mean ± standard deviation).

As anticipated in the materials and methods section, the use of repeated double cross-validation allows obtaining not only a point estimate of the figures of merit on the validation (outer loop) samples but also their confidence intervals, so as to be able to evaluate the consistency of the results.

By looking at Table 5, it is evident how the different models result in different reliabilities, with some presenting rather low classification performances. On the other hand, there are some models which result in an overall accuracy higher than 75%, with a comparable mean correct classification rate (suggesting that the classification performances are not affected by the numerosity of the samples; Table 1). Additionally, the standard deviation of the figures of merit for these models is relatively low (corresponding roughly to one more sample being correctly or wrongly classified with respect to the reported averages), confirming the consistency of the obtained classification. Based on these considerations, only the best models will be discussed in detail in the remainder of this section, namely, Lazio vs. Abruzzo, Lazio vs. Sicily, Abruzzo vs. Calabria and Abruzzo vs. Sicily.

The first model to be examined is the one discriminating Lazio samples from the oils from Abruzzo, for which an overall 76.2 ± 3.9% classification accuracy on the outer loop samples was registered. By looking at the individual sensitivities together with their confidence intervals (73.1 ± 3.9% for Lazio and 81.0 ± 7.7% for Abruzzo), it can be stated that the two categories are predicted comparably well. These results can also be graphically appreciated in Figure 1, where the mean scores of the outer loop samples along the only canonical variate of the model together with their 95% confidence intervals are displayed. It is evident from Figure 1 how almost all the Abruzzo samples have positive scores, while the large majority of Lazio samples are characterized by negative coordinates on the component, indicating a good separation between the categories.

FIGURE 1

FIGURE 1. –PLS-DA model for the discrimination between Lazio and Abruzzo samples: mean scores of the rDCV outer loop samples along the only canonical variate of the model together with their 95% confidence intervals. Legend: red bars–Lazio; blue bars–Abruzzo.

For the sake of interpretation, another advantage of the rDCV procedure is that confidence intervals can also be calculated for model parameters, so as to be able to identify which are the variables that contribute significantly to the discrimination (e.g., by inspecting the values of the associated regression coefficients or of the VIP scores). Moreover, investigating the sign of the regression coefficients also allows postulating whether the associated predictor is more or less concentrated in a category with respect to the other. In particular, the variables found to significantly contribute to the discriminant model were V, Fe, Zn, Rb, antioxidant capacity (all higher in Lazio samples), and Ni and antioxidant activity in the DPPH assay (higher in the oils from Abruzzo).

As far as the Lazio vs. Sicily model is concerned, a slightly higher accuracy was obtained (79.4 ± 3.1%), the individual sensitivities being 81.9 ± 4.9% for Lazio and 76.1 ± 2.8% for Sicily. Analogously to that described above, the discrimination between the two classes can also be visually appreciated in Figure 2, where the mean scores of the outer loop samples along the only canonical variate of the model together with their 95% confidence intervals are displayed.

FIGURE 2

FIGURE 2. –PLS-DA model for the discrimination between Lazio and Sicily samples: mean scores of the rDCV outer loop samples along the only canonical variate of the model together with their 95% confidence intervals. Legend: red bars–Lazio; blue bars–Sicily.

In this case, based on the values of the PLS-DA regression coefficients, all the variables found to significantly contribute to the discriminant model (Na, Mg, P, Ti, Rb, antioxidant activity in the DPPH assay) should be, on average, higher in the oils from Lazio. When considering the model discriminating Abruzzo oils from the Calabrian ones, an overall 81.4 ± 6.1% accuracy on the outer loop samples was obtained, the mean correct classification rate (81.0 ± 7.1%) and the individual sensitivities for the two categories (82.9 ± 4.35% for Abruzzo and 79.1 ± 12.2% for Calabria) being almost equal. In particular, the higher standard deviation of the sensitivity for Calabria is due to the very limited number of samples in that class. When looking at the significant predictors, only five variables (P, V, Fe, Zn, and antioxidant activity) were identified, and the coefficients indicate that they all should be, on average, higher for the Abruzzo samples.

Lastly, the Abruzzo vs. Sicily model resulted in an accuracy of 75.0 ± 4.3%, the sensitivities being 81.2 ± 7.6% for Abruzzo and 69.7 ± 5.9% for Sicily. Inspection of the model parameters led to identifying as significant the contribution of Na, Ni, and antioxidant activity (in the DPPH assay with higher data in the oils from Abruzzo) and of V and Fe, more concentrated in the samples from Sicily.

Classification of EVOOs according to cultivar and to whether it was organically produced.

In a second stage of the study, the possibility of discriminating oil samples according to their cultivar was explored. In this case, given the available information about the samples and the fact that only a relatively small fraction of the analyzed oils was monovarietal, the investigation was restricted to the comparison of Coratina (21 samples) and Frantoio (12 samples) (Table 1). The PLS-DA classification approach was validated through an rDCV strategy as described in the previous section and resulted in an overall accuracy of 68.9 ± 6.2%, and 80.7 ± 8.8% and 61.9 ± 7.6% sensitivities for Frantoio and Coratina, respectively, corresponding to a mean correct classification rate of 71.3 ± 6.2%. Investigation of the model parameters suggested that five variables only, namely, P, Ti, Zn (higher in Coratina), Fe, and Ni (more concentrated in Frantoio), significantly contributed to the discriminant model.

Lastly, the possibility of discriminating whether the oil was organically produced or not was also attempted, but the classification model resulted in a very poor accuracy (close to 50%) suggesting that, at least for the investigated samples, organic cultivation has little impact on the elemental composition with respect to non-organic production.

4 Conclusion

This study showed that the As, Cu, Fe, and Pb levels in the analyzed samples were far below the MRLs, which certifies the high quality of Italian EVOO.

The element concentrations allow to distinguish well some geographical origins of the EVOO samples and also, although slightly less well, the two cultivars Coratina and Frantoio. On the other hand, given the high heterogeneity of the data set, it is not possible to distinguish organic oils from non-organic ones. This is probably due to the fact that within the two classes the variability related to geographical origin and cultivar is added.

This study can be used to create datasets for element levels in EVOOs for each production region to support geographic origin authentication. In the future, other information will have to be considered together with the elemental profile of EVOO such as climatic factors and bioavailable fraction of the total content of elements to further corroborate the use of the elements as a marker of provenance.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

ML-A: conceptualization, investigation, methodology, validation, formal analysis, visualization, data curation, writing—original draft preparation, writing—reviewing and editing and supervision. FM: methodology, validation, formal analysis, writing—original draft preparation and writing—reviewing and editing. MA-F: investigation. LM: formal analysis. ALC: resources. CMM: resources. SC: resources and supervision.

Funding

This work was partially supported by the “Agroalimentare e Ricerca” (AGER) program, project AGER2-Rif.20160169, “Valorization of Italian Olive products through INnovative analytical tools–VIOLIN”.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank Dr. Elisabetta Marconi, and Dr. Giulia Vitiello for their excellent support in the treatment and classification of the samples.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2021.769620/full#supplementary-material

References

Aceto, M., Calà, E., Musso, D., Regalli, N., and Oddone, M. (2019). A preliminary study on the authentication and traceability of extra virgin olive oil made from Taggiasca olives by means of trace and ultra-trace elements distribution. Food Chemistry 298, 125047. doi:10.1016/j.foodchem.2019.125047

PubMed Abstract | CrossRef Full Text | Google Scholar

Astolfi, M. L., Conti, M. E., Marconi, E., Massimi, L., and Canepari, S. (2020a). Effectiveness of Different Sample Treatments for the Elemental Characterization of Bees and Beehive Products. Molecules 25, 4263. doi:10.3390/molecules25184263

PubMed Abstract | CrossRef Full Text | Google Scholar

Astolfi, M. L., Marconi, E., Protano, C., and Canepari, S. (2020b). Comparative elemental analysis of dairy milk and plant-based milk alternatives. Food Control 116, 107327. doi:10.1016/j.foodcont.2020.107327

CrossRef Full Text | Google Scholar

Astolfi, M. L., Marconi, E., Vitiello, G., and Massimi, L. (2021b). An optimized method for sample preparation and elemental analysis of extra-virgin olive oil by inductively coupled plasma mass spectrometry. Food Chemistry 360, 130027. doi:10.1016/j.foodchem.2021.130027

PubMed Abstract | CrossRef Full Text | Google Scholar

Astolfi, M. L., Marotta, D., Cammalleri, V., Marconi, E., Antonucci, A., Avino, P., et al. (2021a). Determination of 40 Elements in Powdered Infant Formulas and Related Risk Assessment. Ijerph 18, 5073. doi:10.3390/ijerph18105073

PubMed Abstract | CrossRef Full Text | Google Scholar

Astolfi, M. L., Protano, C., Marconi, E., Massimi, L., Brunori, M., Piamonti, D., et al. (2020c). A new rapid treatment of human hair for elemental determination by inductively coupled mass spectrometry. Anal. Methods 12, 1906–1918. doi:10.1039/c9ay01871a

CrossRef Full Text | Google Scholar

Bajoub, A., Bendini, A., Fernández-Gutiérrez, A., and Carrasco-Pancorbo, A. (2018). Olive oil authentication: A comparative analysis of regulatory frameworks with especial emphasis on quality and authenticity indices, and recent analytical techniques developed for their assessment. A review. Crit. Rev. Food Sci. Nutr. 58, 832–857. doi:10.1080/10408398.2016.1225666

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakircioglu, D., Kurtulus, Y. B., and Yurtsever, S. (2013). Comparison of extraction induced by emulsion breaking, ultrasonic extraction and wet digestion procedures for determination of metals in edible oil samples in Turkey using ICP-OES. Food Chemistry 138, 770–775. doi:10.1016/j.foodchem.2012.10.089

PubMed Abstract | CrossRef Full Text | Google Scholar

Barker, M., and Rayens, W. (2003). Partial least squares for discrimination. J. Chemometrics 17, 166–173. doi:10.1002/cem.785

CrossRef Full Text | Google Scholar

Beltrán, M., Sánchez-Astudillo, M., Aparicio, R., and García-González, D. L. (2015). Geographical traceability of virgin olive oils from south-western Spain by their multi-elemental composition. Food Chemistry 169, 350–357. doi:10.1016/j.foodchem.2014.07.104

PubMed Abstract | CrossRef Full Text | Google Scholar

Bendini, A., Cerretani, L., Carrasco-Pancorbo, A., Gómez-Caravaca, A., Segura-Carretero, A., Fernández-Gutiérrez, A., et al. (2007). Phenolic Molecules in Virgin Olive Oils: a Survey of Their Sensory Properties, Health Effects, Antioxidant Activity and Analytical Methods. An Overview of the Last Decade Alessandra. Molecules 12 (8), 1679–1719. doi:10.3390/12081679

PubMed Abstract | CrossRef Full Text | Google Scholar

Benincasa, C., Lewis, J., Perri, E., Sindona, G., and Tagarelli, A. (2007). Determination of trace element in Italian virgin olive oils and their characterization according to geographical origin by statistical analysis. Analytica Chim. Acta 585, 366–370. doi:10.1016/j.aca.2006.12.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabrera-Vique, C., Bouzas, P. R., and Oliveras-López, M. J. (2012). Determination of trace elements in extra virgin olive oils: A pilot study on the geographical characterisation. Food Chemistry 134 (1), 434–439. doi:10.1016/j.foodchem.2012.02.088

CrossRef Full Text | Google Scholar

Camin, F., Larcher, R., Perini, M., Bontempo, L., Bertoldi, D., Gagliano, G., et al. (2010). Characterisation of authentic Italian extra-virgin olive oils by stable isotope ratios of C, O and H and mineral composition. Food Chemistry 118, 901–909. doi:10.1016/j.foodchem.2008.04.059

CrossRef Full Text | Google Scholar

Capriotti, A. L., Cavaliere, C., Crescenzi, C., Foglia, P., Nescatelli, R., Samperi, R., et al. (2014). Comparison of extraction methods for the identification and quantification of polyphenols in virgin olive oil by ultra-HPLC-QToF mass spectrometry. Food Chemistry 158, 392–400. doi:10.1016/j.foodchem.2014.02.130

PubMed Abstract | CrossRef Full Text | Google Scholar

Carbone, A., Cacchiarelli, L., and Sabbatini, V. (2018). Exploring quality and its value in the Italian olive oil market: a panel data analysis. Agric. Econ. 6, 6. doi:10.1186/s40100-018-0102-8

CrossRef Full Text | Google Scholar

Chatzistathis, T., Therios, I., and Alifragis, D. (2009). Differential uptake, distribution within tissues, and use efficiency of manganese, iron, and zinc by olive cultivars Kothreiki and Koroneiki. horts 44, 1994–1999. doi:10.21273/HORTSCI.44.7.1994

CrossRef Full Text | Google Scholar

Choe, E., and Min, D. B. (2006). Mechanisms and factors for edible oil oxidation. Comp. Rev. Food Sci. Food Saf. 5, 169–186. doi:10.1111/j.1541-4337.2006.00009.x

CrossRef Full Text | Google Scholar

Christodouleas, D. C., Fotakis, C., Nikokavoura, A., Papadopoulos, K., and Calokerinos, A. C. (2015). Modified DPPH and ABTS assays to assess the antioxidant profile of untreated oils. Food Anal. Methods 8 (5), 1294–1302. doi:10.1007/s12161-014-0005-6

CrossRef Full Text | Google Scholar

Cioffi, G., Pesca, M. S., De Caprariis, P., Braca, A., Severino, L., and De Tommasi, N. (2010). Phenolic compounds in olive oil and olive pomace from Cilento (Campania, Italy) and their antioxidant activity. Food Chemistry 121 (1), 105–111. doi:10.1016/j.foodchem.2009.12.013

CrossRef Full Text | Google Scholar

Codex Stan 33-1981 (2021). Standard for Olive Oils and Olive Pomace Oils. Adopted in 1981Revision: 1989, 2003, 2015. Amendment: 2009, 2013. URL: http://www.fao.org/input/download/standards/88/CXS_033e_2015.pdf (accessed on April, 2021).

Google Scholar

Cordella, C., Moussa, I., Martel, A.-C., Sbirrazzuoli, N., and Lizzani-Cuvelier, L. (2002). Recent developments in food characterization and adulteration detection: Technique-oriented perspectives. J. Agric. Food Chem. 50, 1751–1764. doi:10.1021/jf011096z

PubMed Abstract | CrossRef Full Text | Google Scholar

Damak, F., Asano, M., Baba, K., Suda, A., Araoka, D., Wali, A., et al. (2019). Interregional traceability of Tunisian olive oils to the provenance soil by multielemental fingerprinting and chemometrics. Food Chemistry 283, 656–664. doi:10.1016/j.foodchem.2019.01.082

PubMed Abstract | CrossRef Full Text | Google Scholar

Dugo, L., Russo, M., Cacciola, F., Mandolfino, F., Salafia, F., Vilmercati, A., et al. (2020). Determination of the phenol and tocopherol content in italian high-quality extra-virgin olive oils by using LC-MS and multivariate data analysis. Food Anal. Methods 13, 1027–1041. doi:10.1007/s12161-020-01721-7

CrossRef Full Text | Google Scholar

European Union (EU) (20122012). The Commission Regulation No. 1151/2012 of 21 November 2012 on quality schemes for agricultural products and foodstuffs OJ L 343/1, 14.12. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32012R1151&from=EN (accessed on August, 2021).

Google Scholar

Eurostat (2019). Olive trees - Area by age and density classes (area in ha). Available online: http://appsso.eurostat.ec.europa.eu/nui/show.do?dataset=orch_olives3&lang=en (accessed on August, 2021).

Google Scholar

Farmaki, E. G., Thomaidis, N. S., Minioti, K. S., Ioannou, E., Georgiou, C. A., and Efstathiou, C. E. (2012). Geographical characterization of Greek olive oils using rare earth elements content and supervised chemometric techniques. Anal. Lett. 45, 920–932. doi:10.1080/00032719.2012.655656

CrossRef Full Text | Google Scholar

Filzmoser, P., Hron, K., and Reimann, C. (2009). Principal Component Analysis for Compositional Data With Outliers. Environmetrics: The Official Journal of the International Environmetrics Society 20 (6), 621–632. doi:10.1002/env.966

CrossRef Full Text | Google Scholar

Frankel, E. N. (1985). “Chemistry of autoxidation: mechanism, products and flavor significance,” in Flavor Chemistry of Fats and Oils. Editors D. B. Min, and T. H. Smouse (Champaign, IL (USA): AOCS Press), 1–37.

Google Scholar

Frezzini, M. A., Castellani, F., De Francesco, N., Ristorini, M., and Canepari, S. (2019). Application of DPPH assay for assessment of particulate matter reducing properties. Atmosphere 10 (12), 816. doi:10.3390/atmos10120816

CrossRef Full Text | Google Scholar

Giaccio, M., and Vicentini, A. (2008). Determination of the geographical origin of wines by means of the mineral content and the stable isotope ratios: A review. J. Comm. Sci. Technol. Qual. 47 (I–IV), 267–284.

Google Scholar

Giacomo, D., Pera Lara, L., Daniele, G., Francesco, S., and Turco Vincenzo, L. (2004). Influence of the olive variety and the zone of provenience on selenium content determined by cathodic stripping potentiometry (CSP) in virgin olive oils. Food Chemistry 88, 135–140. doi:10.1016/j.foodchem.2003.12.036

CrossRef Full Text | Google Scholar

Gonzalez-Fernandez, I., Iglesias-Otero, M. A., Esteki, M., Moldes, O. A., Mejuto, J. C., and Simal-Gandara, J. (2019). A critical review on the use of artificial neural networks in olive oil production, characterization and authentication. Crit. Rev. Food Sci. Nutr. 59 (12), 1913–1926. doi:10.1080/10408398.2018.1433628

PubMed Abstract | CrossRef Full Text | Google Scholar

Gumus, Z. P., Celenk, V. U., Tekin, S., Yurdakul, O., and Ertas, H. (2017). Determination of trace elements and stable carbon isotope ratios in virgin olive oils from Western Turkey to authenticate geographical origin with a chemometric approach. Eur. Food Res. Technol. 243, 1719–1727. doi:10.1007/s00217-017-2876-4

CrossRef Full Text | Google Scholar

Hannachi, H., and Elfalleh, W. (2020). Enrichment of olive oil with polyphenols from oleaster leaves using central composite design for the experimental measurements. Anal. Lett. 54, 590–607. doi:10.1080/00032719.2020.1774599

CrossRef Full Text | Google Scholar

International Olive Council (IOC) (2018a). IOC data for the 2017/18 crop year show a year-on-year increase in the production of olive oil. Available online: http://www.internationaloliveoil.org/news/view/698-year-2018-news/1049-ioc-data-for-the-2017-18-crop-year-show-a-year-on-year-increase-in-the-production-of-olive-oil (accessed on August, 2021).

Google Scholar

International Olive Council (IOC) (2019). Trade standard applying to olive oils and olive-pomace oils. COI/T.15/NC No. 3/Rev. 15. Available online: https://www.internationaloliveoil.org/wp-content/uploads/2020/07/Trade-standard-T15-NC3-Rev15-EN.pdf (accessed on April, 2021).

Google Scholar

International Olive Council (IOC) (2018b). World Olive Oil Figures. Available online: http://www.internationaloliveoil.org/estaticos/view/131-world-olive-oilfigures (accessed on August, 2021).

Google Scholar

Kabata-Pendias, A. (2010). Trace elements in soils and plants. 4th ed. Boca Raton: CRC Press. doi:10.1201/b10158

CrossRef Full Text | Google Scholar

Kedare, S. B., and Singh, R. P. (2011). Genesis and development of DPPH method of antioxidant assay. J. Food Sci. Technol. 48 (4), 412–422. doi:10.1007/s13197-011-0251-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelly, S., Heaton, K., and Hoogewerff, J. (2005). Tracing the geographical origin of food: The application of multi-element and multi-isotope analysis. Trends Food Sci. Tech. 16, 555–567. doi:10.1016/j.tifs.2005.08.008

CrossRef Full Text | Google Scholar

Lepri, F. G., Chaves, E. S., Vieira, M. A., Ribeiro, A. S., Curtius, A. J., DeOliveira, L. C. C., et al. (2011). Determination of Trace Elements in Vegetable Oils and Biodiesel by Atomic Spectrometric Techniques-A Review. Appl. Spectrosc. Rev. 46, 175–206. doi:10.1080/05704928.2010.529628

CrossRef Full Text | Google Scholar

Llorent-Martínez, E. J., Fernández-de Córdova, M. L., Ortega-Barrales, P., and Ruiz-Medina, A. (2014). Quantitation of Metals During the Extraction of Virgin Olive Oil From Olives Using ICP-MS After Microwave-Assisted Acid Digestion. J. Am. Oil Chem. Soc. 91, 1823–1830. doi:10.1007/s11746-014-2511-5

CrossRef Full Text | Google Scholar

Llorent-Martínez, E. J., Ortega-Barrales, P., Fernández-de Córdova, M. L., and Ruiz-Medina, A. (2011). Analysis of the legislated metals in different categories of olive and olive-pomace oils. Food Control 22, 221–225. doi:10.1016/j.foodcont.2010.07.002

CrossRef Full Text | Google Scholar

Luka, M. F., and Akun, E. (2019). Investigation of trace metals in different varieties of olive oils from northern Cyprus and their variation in accumulation using ICP-MS and multivariate techniques. Environ. Earth Sci. 78, 578. doi:10.1007/s12665-019-8581-9

CrossRef Full Text | Google Scholar

Mendil, D., Uluözlü, Ö. D., Tüzen, M., and Soylak, M. (2009). Investigation of the levels of some element in edible oil samples produced in Turkey by atomic absorption spectrometry. J. Hazard. Mater. 165, 724–728. doi:10.1016/j.jhazmat.2008.10.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Pérez, N. F., Ferré, J., and Boqué, R. (2009). Calculation of the reliability of classification in discriminant partial least-squares binary classification. Chemometrics Intell. Lab. Syst. 95, 122–128. doi:10.1016/j.chemolab.2008.09.005

CrossRef Full Text | Google Scholar

Pérez-Jiménez, J., Arranz, S., Tabernero, M., Díaz- Rubio, M. E., Serrano, J., Goñi, I., et al. (2008). Updated methodology to determine antioxidant capacity in plant foods, oils and beverages: Extraction, measurement and expression of results. Food Res. Int. 41 (3), 274–285. doi:10.1016/j.foodres.2007.12.004

CrossRef Full Text | Google Scholar

Perna, A., Simonetti, A., Intaglietta, I., Sofo, A., and Gambacorta, E. (2012). Metal content of southern Italy honey of different botanical origins and its correlation with polyphenol content and antioxidant activity. Int. J. Food Sci. Technol. 47, 1909–1917. doi:10.1111/j.1365-2621.2012.03050.x

CrossRef Full Text | Google Scholar

Pohl, W. L. (2011). Economic geology: Principles and practice. New Jersey: Wiley-Blackwell. 978-1-444-33663-4. doi:10.1002/9781444394870

CrossRef Full Text | Google Scholar

Pošćić, F., Furdek Turk, M., Bačić, N., Mikac, N., Bertoldi, D., Camin, F., et al. (2019). Removal of pomace residues is critical in quantification of element concentrations in extra virgin olive oil. J. Food Compost. Anal. 77, 39–46. doi:10.1016/j.jfca.2019.01.002

CrossRef Full Text | Google Scholar

Russo, G., Beritognolo, I., Bufacchi, M., Stanzione, V., Pisanelli, A., Ciolfi, M., et al. (2020). Advances in biocultural geography of olive tree (Olea europaea L.) landscapes by merging biological and historical assays. Sci. Rep. 10, 7673. doi:10.1038/s41598-020-64063-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Šarolić, M., Gugić, M., Tuberoso, C., Jerković, I., Šuste, M., Marijanović, Z., et al. (2014). Volatile Profile, Phytochemicals and Antioxidant Activity of Virgin Olive Oils from Croatian Autochthonous Varieties Mašnjača and Krvavica in Comparison with Italian Variety Leccino. Molecules 19 (1), 881–895. doi:10.3390/molecules19010881

PubMed Abstract | CrossRef Full Text | Google Scholar

Shah, N. S., and Soylak, M. (2021). Advanced methodologies for trace elements in edible oil samples: a review. Crit. Rev. Anal. Chemistry 15, 1–20. doi:10.1080/10408347.2021.1895710

CrossRef Full Text | Google Scholar

Ståhle, L., and Wold, S. (1987). Partial least squares analysis with cross-validation for the two-class problem: A Monte Carlo study. J. Chemometrics 1, 185–196. doi:10.1002/cem.1180010306

CrossRef Full Text | Google Scholar

Tchounwou, P. B., Yedjou, C. G., Patlolla, A. K., and Sutton, D. J. (2012). Heavy Metal Toxicity and the Environment. Exp. Suppl. 101, 133–164. doi:10.1007/978-3-7643-8340-4_6

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiruvengadam, M., Ghimire, B. K., Kim, S.-H., Yu, C. Y., Oh, D.-H., Chelliah, R., et al. (2020). Assessment of mineral and phenolic profiles and their association with the antioxidant, cytotoxic effect, and antimicrobial potential of Lycium chinense Miller. Plants 9 (8), 1023. doi:10.3390/plants9081023

PubMed Abstract | CrossRef Full Text | Google Scholar

Trindade, A. S., Dantas, A. F., Lima, D. C., Ferreira, S .L., and Teixeira, L. S. (2015). Multivariate Optimization of Ultrasound-Assisted Extraction for Determination of Cu, Fe, Ni, and Zn in Vegetable Oils by High-Resolution Continuum Source Atomic Absorption Spectrometry. Food Chem. 15 (185), 145–150. doi:10.1016/j.foodchem.2015.03.118

CrossRef Full Text | Google Scholar

Visioli, F., Bellomo, G., and Galli, C. (1998). Free radical-scavenging properties of olive oil polyphenols. Biochem. Biophysical Res. Commun. 247 (1), 60–64. doi:10.1006/bbrc.1998.8735

PubMed Abstract | CrossRef Full Text | Google Scholar

Wold, S., Martens, H., and Wold, H. (1983). “The multivariate calibration problem in chemistry solved by the PLS method,” in Matrix Pencils. Lecture Notes in Mathematics. Editors B. Kågström, and A. Ruhe. 1st ed. (Germany: Springer, Berlin/Heidelberg), 973, 286–293. doi:10.1007/BFb0062108

CrossRef Full Text | Google Scholar

Yaşar, S. B., Baran, E. K., and Alkan, M. (2012). “Metal determinations in olive oil, Olive oil - Constituents, quality, health properties and bioconversions,” in InTech. Editor Dr. Dimitrios Boskou. Available online: http://www.intechopen.com/books/olive-oil-constituents-quality-healthproperties-and-bioconversions/metal-determinations-in-olive-oil (accessed on July, 2021).

Google Scholar

Zaroual, H., Chénè, C., El Hadrami, E. M., and Karoui, R. (2021). Application of new emerging techniques in combination with classical methods for the determination of the quality and authenticity of olive oil: a review. Crit. Rev. Food Sci. Nutr. 1, 1–24. doi:10.1080/10408398.2021.1876624

CrossRef Full Text | Google Scholar

Zeiner, M., Juranovic-Cindric, I., and Škevin, D. (2010). Characterization of extra virgin olive oils derived from the Croatian cultivar Oblica. Eur. J. Lipid Sci. Technol. 112, 1248–1252. doi:10.1002/ejlt.201000006

CrossRef Full Text | Google Scholar

Zeiner, M., Steffan, I., and Cindric, I. J. (2005). Determination of trace elements in olive oil by ICP-AES and ETA-AAS: A pilot study on the geographical characterization. Microchemical J. 81, 171–176. doi:10.1016/j.microc.2004.12.002

CrossRef Full Text | Google Scholar