Next Article in Journal
Effect of Physical Exercise and Genetic Background on Glucose Homeostasis and Liver/Muscle Proteomes in Mice
Previous Article in Journal
In Vitro Metabolic Study of Four Synthetic Cathinones: 4-MPD, 2-NMC, 4F-PHP and bk-EPDP
Previous Article in Special Issue
Stability of Wheat Floret Metabolites during Untargeted Metabolomics Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimisation of the HS-SPME/GC-MS Approach by Design of Experiments Combined with Chemometrics for the Classification of Cretan Virgin Olive Oils

by
Artemis Lioupi
1,2,3,
Ioannis Sampsonidis
2,3,4,
Christina Virgiliou
1,2,3,
Vassiliki T. Papoti
5,
Kyriaki G. Zinoviadou
5,
Apostolos Spyros
6 and
Georgios Theodoridis
1,2,3,*
1
Laboratory of Analytical Chemistry, School of Chemistry, Aristotle University of Thessaloniki, GR-54124 Thessaloniki, Greece
2
Biomic AUTh, Center for Interdisciplinary Research and Innovation (CIRI-AUTH), Balkan Center B1.4, 10th km Thessaloniki-Thermi Rd, P.O. Box 8318, GR-57001 Thessaloniki, Greece
3
FoodOmicsGR Research Infrastructure, AUTh Node, Center for Interdisciplinary Research and Innovation (CIRI-AUTH), Balkan Center B1.4, 10th km Thessaloniki-Thermi Rd, P.O. Box 8318, GR-57001 Thessaloniki, Greece
4
Department of Nutritional Sciences and Dietetics, International Hellenic University, GR-57400 Thessaloniki, Greece
5
Department of Food Science and Technology, Perrotis College, American Farm School, GR-55102 Thessaloniki, Greece
6
NMR Laboratory, Department of Chemistry, University of Crete, Voutes Campus, P.O. Box 2208, GR-71003 Heraklion, Crete, Greece
*
Author to whom correspondence should be addressed.
Metabolites 2022, 12(2), 114; https://doi.org/10.3390/metabo12020114
Submission received: 5 January 2022 / Revised: 19 January 2022 / Accepted: 20 January 2022 / Published: 25 January 2022
(This article belongs to the Special Issue Metabolomics Methodologies and Applications II)

Abstract

:
A headspace-solid phase microextraction/gas chromatography-mass spectrometry (HS-SPME/GC-MS) method was developed herein for the analysis of virgin olive oil volatile metabolome. Optimisation of SPME conditions was performed by Design of Experiments (DoE) and Response Surface Methodology (RSM) approaches and factors, such as sample volume, sample stirring, extraction temperature and time, and desorption temperature and time, were examined to reach optimal microextraction conditions. The potential of the optimised method was then investigated for its use in the classification of Cretan virgin olive oil samples with the aid of multivariate statistical analysis. Certain markers were identified with significance in the geographical classification of Cretan extra-virgin olive oil (EVOO) samples. In total, 92 volatile organic compounds were tentatively identified and semi-quantified, and the data obtained confirm that the method is robust, reliable, and analytically powerful for olive oil classification.

1. Introduction

Extra virgin olive oil is one of the most significant ingredients of the Mediterranean diet, and its consumption has been significantly increased due to its high and special nutritional and sensory characteristics [1]. Virgin olive oil (VOO) is classified to quality categories as described by official standards (e.g., Codex Alimentarius, International Olive Council, and European Union) where sensory assessment is also included [2]. The sensory quality of VOO is determined by a panel test that evaluates positive and negative descriptors, according to European Union (EU) criteria (EU regulations) [3,4]. Olive oils possessing the fruity attribute and lacking any sensory defects can be categorized as “extra virgin”, whilst the existence and severity of sensory defects are used to verify whether the oil could be categorized as “virgin” or lampante [5]. Volatile compounds in olive oil can be linked to both positive and negative sensory characteristics and have a major effect on the oil’s quality and, as a result, on consumer preference. Hence, determining the VOO volatile fraction is important to ensure consumer protection, but also to avoid unfair competition, which can destabilize the market and disrupt regional/national economies.
In Greece, over 30 olive varieties are mainly produced [6]. Although Greece ranks third in the world in VOO production, the high-quality potential of Greek olive oils has been scarcely examined, and such information has not been communicated to the market. Crete is one of the main olive oil-producing areas in Greece and two oil cultivars, Koroneiki and Tsounati, mainly originate from there. Koroneiki, one of the most well-known and appraised olive cultivars in the world, is the predominant olive oil cultivar in Greece, accounting for about 60% of the country’s olive oil production (approximately 250,000 t) [6]. The Tsounati variety (also referred to as “Athinolia” or “Mastoeidis”) is less exploited and is cultivated in specific areas of Peloponnese, mainly in south Lakonia, Argolida, as well as in western Crete [7]. The quality potential of both examined cultivars can be partially revealed from the quality awards they gain from international competitions (e.g., Athena International Olive Oil Competition). It is necessary to verify the quality standards of Greek traditional high-added-value VOO to protect the final consumer from non-authentic food products and to help Greek food producers, suppliers, and service establishments compete in national and international markets [8,9,10].
Gas chromatography-mass spectrometry (GC-MS) has been commonly applied for the analysis of volatile organic compounds (VOCs) in EVOO, with an emphasis given in biomarker-detection, authenticity control, and issues associated with food quality and safety [11,12,13]. Headspace solid-phase microextraction (HS-SPME) has a number of advantages over other extraction methods, including automation and simplicity. Due to these advantages, it is usually the method of choice, with many applications in the field [14,15,16]. However, for this type of analysis, it is important to identify the role of all the possible factors that can influence the extraction of volatile compounds. In the present publication, sample volume, sample stirring, extraction temperature and time, and desorption time are the factors examined to reach optimal microextraction conditions. According to our knowledge, there are only two studies that focus on the optimisation of HS-SPME parameters by utilising advanced statistical techniques such as Design of Experiments (DOE) for VOO volatile compound analysis [17,18]. The presented approach is the first that optimises five HS-SPME factors simultaneously for EVOO analysis.
The overall goal was the development of an HS-SPME/GC-MS method for the analysis of a VOO volatile profile. For this purpose, an optimisation of SPME conditions was performed using DOE and Response Surface Methodology (RSM) approaches. The method was then applied for the analysis of VOOs originating from different regions of Crete, Southern Greece. This study was focused on the potential use of the optimised HS-SPME/GC-MS method combined with chemometrics for sample classification. The method proved to be robust, enabling the identification and semi-quantification of VOO VOCs.

2. Results and Discussion

2.1. Optimisation of SPME Extraction Conditions

SPME sorption is dependent among other parameters on the sample volume, and it has been observed that extraction efficiency is improved when the volume of HS in the vial is minimised [19]. However, the extraction efficiency of volatile compounds may show a non-linear relation toward the sample volume due to fiber overload. Hence, the sample mass of the sample was examined under four different quantities (1.000, 1.850, 3.000, and 4.000 g). As illustrated in Figure 1a, sample amounts of 1.850 and 3.000 g improved extraction efficiency. However, a further increase from 3.000 to 4.000 g of sample had a negative effect on the extraction efficiency, yielding an optimal sample amount of 1.850 g.
For the optimisation of stirring, four different stirring speeds (100, 150, 200, 250, and 300 rpm) were evaluated. As shown in Figure 1b, after sample stirring at 250 rpm, no further improvement in extraction efficiency was observed. Hence, a stirring speed of 250 rpm was used for further experiments.
Although included in the DOE-RSM optimisation as a studied factor, a preliminary screening experiment was performed to determine the behaviour of the most critical HS-SPME condition for EVOO analysis, which is extraction temperature. Selection of an appropriate extraction temperature is crucial to minimise the oxidation of the olive oil sample inside the SPME vial, preventing the release of unwanted oxidation components [20]. More specifically, in unheated or heated-at-a-lower-temperature (50 °C) olive oil samples, a higher content of alcohols, esters, and terpenes has been observed compared to oils heated at higher temperatures; in the latter, higher amounts of aldehydes and carboxylic acids were found [21]. During oxidation off-flavor compounds produced due to fatty acid deterioration can be detected [22]. Consequently, the initially pleasant sensory characteristics of the oil attributed to compounds derived from the lipoxygenase (LOX) pathway eventually give way to unpleasant sensory attributes due to volatile compound deterioration. The screening of this parameter was performed by monitoring the total peak area of the selected components that are indicative of olive oil oxidation, while also monitoring components that are generally desirable according to the literature [21,22,23,24]. The selected compounds that indicate increased oil oxidation and must be kept at a minimum during sample heating were heptanal, 2-heptenal, octanal, and 2,4-heptadienal. The compounds that were selected as desirable for the analysis and contribute to its green odour notes [25] were 3-hexenal, (2E)-hexenal, 2-hexene-1-ol, hexan-1-ol, and 1-penten-3-one.
The study of extraction temperature showed an increase for both the desirable and undesirable oxidation components (Figure 1c). As indicated, there is a steep increase for the desirable components maximizing at around 55 °C, while the increase for the oxidation components appears to take place also around 55 °C, showing continuous increasing trend for the whole tested temperature range. This is in accordance with other related studies [13,21], where extraction temperatures above 55 °C were not considered and were not studied as a means to avoid thermal volatile alterations and oxidative degradation of the olive oil.
Regarding DOE-RSM optimization, a quadratic model including all linear terms, two-way interaction terms and second-order terms was fitted to the response data. The Multiple R-squared, Adjusted R-squared, p-value, and output of the Analysis of Variance (ANOVA) table for the model can be seen in Table 1.
As indicated by the Multiple R-squared, Adjusted R-squared, and p-value, the quadratic model is a good fit to the data and is statistically significant. There is also no evidence of lack of fit in the model. By studying Table S1, one can see that the most significant terms are extraction temperature and extraction time. The interaction between desorption temperature and desorption time also appears to be of significance.
A contour plot for the response surface of extraction time and extraction temperature can be seen in Figure 2. The plot indicates that the optimal point for these two factors is close to their maximum values. It also indicates that extraction time matters less as extraction temperature increases and vice-versa. For example, at a temperature of 85 °C, the same effect is achieved with an extraction time between 40 and 65 min. However, the initial screening study indicated that at temperatures above 55 °C an increase in oxidation products is witnessed. Therefore, a safer factor level selection would be an extraction temperature of 55 °C and an extraction time of 50 min.
As the interaction between desorption time and desorption temperature appears significant, by studying their contour plot at Figure S1 at the selected extraction time and temperature, one can see a saddle-like plot, which suggests that optima lie with one factor maxed or the other. To keep desorption temperature as low as possible to avoid any oxidation or thermal degradation effects at the inlet, a selection of desorption temperature of 240 °C and desorption time of 12 min is made.
Regarding conditioning time, all the models’ terms for this factor appear not to be significant, so a choice of the shortest time is made, in this case 5 min. In summary, the SPME conditions chosen for analysis are 55 °C of extraction temperature, 5 min of conditioning time, 50 min of extraction time, and 240 °C and 12 min desorption temperature and time, respectively.
The aim of the present study was the development of a method for the comprehensive profiling of the volatile content of VOOs. The oven program seemed to be appropriate for the most volatile components, since at 17 min and 180 °C all components of interest were eluted. The DVB/CAR/PDMS fiber showed high extraction efficiency and was selected for further study. Extraction temperature and extraction time had a great effect on the detection of volatile components and their composition. With an extraction temperature of 55 °C and an extraction time of 50 min, an increased area of the total ion chromatogram was observed, preventing thermal oxidation or degradation of the compounds of interest. The optimum HS-SPME/GC-MS conditions were applied in the profiling method, facilitating the identification and semi-quantitation of 92 VOCs of different chemical classes, as shown in Table 2. A representative gas chromatogram of a Quality Control (QC) sample analysed with the optimal HS-SPME/GC-MS conditions is shown in Figure 3.

2.2. Identification of VOCs by HS-SPME/GC-MS Analysis

Olive oil volatile compounds are primarily derived from the enzymatic oxidation of the fatty acids linoleic and linolenic. Positive aroma perceptions in olive oil are produced by endogenous plant enzymes, through the LOX, while aroma defects are linked with chemical oxidation and exogenous enzymes [5]. The volatile fraction of VOOs mainly consists of aldehydes, ketones, aliphatic and aromatic hydrocarbons, ethers, esters, aliphatic and triterpenic alcohols, and furan derivatives. Aroma compounds responsible for green desirable notes, found in high-quality VOOs, are produced enzymatically from polyunsaturated fatty acids through the LOX pathway. It has been observed that in the aroma of these oils, C6 aldehydes and alcohols and their corresponding esters are the most abundant accumulation products [26]. Both major and minor VOCs are important to olive oil aroma, since minor compounds with low odor threshold values are frequently considered significant contributors. Even those whose levels are below their olfactory threshold can still provide useful quality marker information [27].
In all samples, C6 and C5 VOCs derived from the LOX pathway were the most abundant (Figure 4a,b). This is in accordance with the organoleptic results obtained by the tasting panel, where samples were identified as EVOOs with medium fruitiness attributes (3 < x < 6). C6 and C5 VOCs contribute to the delicate green aroma of VOO, which enhances the “green” fruitiness attribute and is highly valued by consumers. The formation of these compounds is mainly influenced by the cultivar, the degree of ripeness of the fruit, the storage period of fruits before oil extraction, and the processing method [28].
C6 aldehydes produced from the linolenic acid, E-2-hexanal, 3-hexenal, and Z-2-hexanal, were found in considerable amounts. In agreement with data reported in other studies examining European olive oil cultivars, (2E)-hexenal, responsible for the “green” attributes, was the most prevalent aroma compound in all VOO samples [20]. More specifically, concerning the Koroneiki cultivar, (2E)-hexenal ranged from 7.64 mg kg−1 in a conventional sample to 51.25 mg kg−1 in an organic sample from the Heraklion region unit. In a similar manner, Z-3-hexanal, which is positively correlated with fruity, leave/grass, and almond notes [29], was found in high levels in VOO samples, reaching 8.77 mg kg−1 in an organic sample from Chania. High amounts of hexanal, produced by breaking 13-hydroperoxide from linoleic acid, are frequently related to sensory notes sweet, apple, and green [30]. In the studied samples, hexanal reached 8.04 mg kg−1 in an organic sample. C6 alcohols such as Z-3-hexen-1-ol, E-2-hexen-1-ol, and hexan-1-ol, related to sweetness, were also identified [5]. Z-3-hexen-1-ol and hexan-1-ol were then esterified by the alcohol acetyltransferase to produce the corresponding esters (3Z)-hex-3-en-1-yl acetate and hexyl acetate, which also were found in considerable amounts in the studied samples ranging from 2.13 to 40.08 mg kg−1 and from 0.5 to 9.75 mg kg−1, respectively. The contribution of these volatiles is quite significant, as they complement aromas with sweet and pleasant notes [31].
Aside from C6 molecules, the headspace of VOO contains a reasonable amount of C5 alcohols and C5 carbonyl compounds [32]. The existence of these VOCs indicates the development of a new branch of the LOX pathway that produces C5 compounds. This extra branch is active during the biogeneration of olive oil aroma. C5 aldehydes and alcohols contribute to the positive flavor attributes of olive oil, providing pungent sensations, and are correlated to bitterness [5]. The molecule 1-penten-3-ol is related to “mouldy” and “rancid” defects associated with lawn, olive, leaf, and pungency and ranged from 0.90 mg kg−1 in conventional VOO of Heraklion to 2.94 mg kg−1 in an organic sample of Koroneiki cultivar in Chania. (Z)-2-penten-1-ol (0.60–1.87 mg kg−1) and 2-pentenal (E)- (0.1–0.51 mg kg−1) are minor compounds whose presence is correlated with banana, sweet, fruity and green, apple, and bitter almond attributes, respectively [33]. C5 ketones, pentene dimers, or monoterpenes affect olive oil aroma even in low concentrations. The molecule 1-penten-3-one reached high amounts at the Koroneiki samples of Heraklion (5.75 mg kg−1) whereas it was not detected in an organic VOO of the Koroneiki variety from Rethymnon. Its presence is linked to leaf, mustard, and pungent notes [34], whereas higher amounts have been also related to metallic notes. The isomer 3-ethyl-1,5-octadiene is produced over the oil extraction process through the same branch of the LOX pathway as the C5 compounds and was detected in all samples (0.12–4.91%) [35]. Other minor volatiles, terpene hydrocarbons (mono- and sesquiterpenes), were detected in Cretan VOOs, and despite the fact they do not contribute to VOO aroma due to their low concentration and high odor threshold, they constitute significant quality markers [29]. Among all terpenoid compounds, the most abundant in the studied VOO samples, were d-limonene (0.04–20.18 mg kg−1) with citrus, mint odor perceptions, α-farnesene (0.02–0.27 mg kg−1) with floral, herb, wood, and sweet odor qualities [36], and (E)-4,8-dimethyl-1,3,7-nonatriene (0.04–0.19 mg kg−1). Other terpenoid volatiles such as alpha- and beta-pinene, o-cymene, and copaene were also detected in lower amounts. Some rare terpenic compounds such as citral, (+)-cyclosativene, eremophilene, neral, 3-carene, and kessane, which so far have been reported only in a few studies [29,37,38], were also identified with the current method.
The detection of benzenes, toluene, xylenes, and styrene (BTEXS) in olive oil is of high importance. These VOCs are widely distributed in the environment and food products mainly due to emissions from vehicles, bonfires, and paints into the ambient air near orchards. In olive oil, BTEXS are frequently identified, and their presence is attributed to several factors, such as biological processes in the fruit, the production technology used, contamination by fuel vapors, etc. [39]. In this study, toluene, styrene, xylene, propylbenzene, 1-ethyl-2-methylbenzene, 4-ethyltoluene, 1,4-diethylbenzene, and naphthalene were identified. The presence of styrene in VOO samples has been attributed to natural causes, such as decarboxylation of the cinnamic acid that is naturally present in the olive pulp or by migration from the plastic package [39]. Naphthalene is usually hanged on trees by olive farmers, as it has a characteristic aroma that is theorised to repel the olive fly (Bactrocera oleae), and in this study was identified in low levels (0.01–0.92 mg kg−1) in 8 out of 63 samples [40].
According to our knowledge, VOCs such as 1-propanone-1-phenyl, benzene, 1,4-diethyl, 11-hexadecen-1-ol, (Z), and 2,4-di-tert-butylphenol are reported for the first time in VOO samples.

2.3. Classification of Cretan VOO Samples

EVOO flavor is considered a significant quality but may also be utilized as a classification criterion. VOO volatile profile plays a fundamental role in defining the geographical origin of oils. The majority (roughly 70%) of PDO/PGI olive oils in Greece are produced in the Peloponnese and Crete. Six of the ten Cretan PDOs come exclusively from Koroneiki, the main Greek cultivar and also the most prominent in Crete; three include blends of 90% Koroneiki VOO, and only one contains blends of 60% oil from Tsounati, the second most significant Cretan cultivar [41]. One PDO was registered for the entire Lasithi regional unit, whilst two PDOs from Chania come from two distinct terroirs but include comparable blends of Koroneiki and Tsounati oils. It is essential that these specific PDOs be subjected to thorough administrative control. These objective criteria, including organoleptic ones, should be accessible to support the authenticity and prevent attempts at fraud.
In the current study, the potential of the newly developed method was investigated for its use in the classification of Cretan VOO samples with the aid of multivariate statistical analysis. The obtained data were analysed using multivariate statistical analysis models to identify correlations between VOC profiles and VOO geographical origin; thus, chromatographic peak areas were used for modelling the differences between the studied groups. The overall performance of the analytical method was evaluated by assessing the variability of the detected compound in QC samples. All the compounds in QC samples presented % RSD < 30%, indicating satisfactory stability and reproducibility of the analytical system during the analytical batch. As shown in Figure S2, principal component analysis (PCA) was used to evaluate the analytical system’s suitability and reproducibility (QC cluster). Group differentiation and identification of potential markers related to geographical origin were performed by latent structures OPLS DA analysis.
Before the assessment of geographical origin to the volatile fingerprint of VOOs, the effect of farming type (organic vs. conventional) was examined. Multivariate statistical analysis was performed to samples from each area produced with different cultivation types. It was observed that farming type is not a statistically significant factor for the specific samples, and thus geographical origin assessment was further performed using both sample types (organic and conventional). Initially, VOO samples from nearby production regions such as Chania-Rethymnon, Rethymnon-Heraklion and Lasithi-Heraklion (Figure S3) did not present any statistically significant differences regarding their volatile content. Although the region of Chania is located on the western side of Crete and Heraklion roughly in the middle of the island, the content of VOOs was not affected by the different origins due probably to other stronger common factors that may affect the VOO profile, which may include climate, fruit maturity, and processing technology [5]. The impact of geographical origin was found significant when VOOs from Lasithi, the easternmost region of Crete, were compared with samples from Chania and Rethymnon. The respective orthogonal projection to latent structures discriminant analysis (OPLS-DA) models are illustrated in Figure 5a,b, where samples with the same origin are clustered together and show clear differentiation against samples from Lasithi. Clear separation was observed between VOOs from Chania and Lasithi and between VOOs from Lasithi and Rethymnon with cross-validated ANOVA analysis (CV Analysis), p = 0.0008 and p = 0.02, respectively. Permutation tests showed that the models were robust, with high predictability (R2Y (cum) and Q2 values of 0.827 and 0.656 and 0.881 and 0.453, respectively). Studying the multivariate loadings vectors in each case revealed various analytes related to the differentiation between the studied groups. As a final step, the unpaired Student’s t-test was performed to establish and confirm the significance of important metabolites. All statistically significant metabolites and their variations in each group are presented as box plots in Figure 6.
The compounds found significant (p < 0.05) for geographical differentiation consisted of four terpenic hydrocarbons (6-methyl-5-hepten-2-one, copaene, (E)-4,8-dimethyl-1,3,7-nonatriene, and (Z)-beta-ocimene), two esters ((3Z)-hex-3-en-1-yl acetate and hexyl acetate), two aldehydes ((E)-2-pentenal and pentanal), three hydrocarbons (dodecene, (E)-2-dodecene, undecene, and 1-ethyl-2-methylbenzene), one alcohol ((Z)-11-hexadecen-1-ol), and one ether (1-methoxy-2-propanol), as shown in Figure 6. Eight of them, 6-methyl-5-hepten-2-one, copaene, (E)-4,8-dimethyl-1,3,7-nonatriene, (Z)-beta-ocimene, (3Z)-hex-3-en-1-yl acetate, hexyl acetate, (E)-2-pentenal, and dodecene, have previously been reported as markers for the classification of Greek VOO according to geographical origin [42,43]. More specifically, efforts have been put on studies that focus on the differentiation of olive oils produced in different islands of north-western Greece [43] (Zakynthos, Kefalonia, Lefkada, Kerkyra, Preveza, and Messologi) and different regions [42] (Messinia and Lakonia in the Peloponnese, Irakleio on the island of Crete, and Etoloakarnania in central Greece).
In this study, when samples from Lasithi and Rethymnon were compared, seven volatiles were found to be responsible for the differentiation between the groups. Three of them ((3Z)-hex-3-en-1-yl acetate, hexyl acetate and (E)-2-pentenal) are produced through the LOX-pathway, contributing to olive oil aroma. The volatiles 1-methoxy-2-propanol and 1-dodecene have also been identified as potential markers in Croatian monovarietal/PDO EVOO [44], while pentanal has shown high discrimination ability for the classification of north Moroccan olive oils [45]. A ripening indicator with fruity-like odor, 6-Methyl-5-hepten-2-one, has been detected as a marker for the characterisation of Turkish olive oils according to origin (Mediterranean, Aegean, South-eastern Anatolia, Marmara, and Black Sea) [46]. Overall, (3Z)-hex-3-en-1-yl acetate, and hexyl acetate were found upregulated in samples from Lasithi while all the others were increased in samples from Rethymnon. Pentanal was found as statistically significant in both origin studies of VOOs. Samples from Lasithi showed lower levels of pentanal compared to Rethymnon and Chania.
In total, eight compounds were found to differ between VOOs from Chania and Lasithi. Three terpenic hydrocarbons, copaene, with a sweet-fruity description, (E)-4,8-dimethyl-1,3,7-nonatriene, and (Z)-beta-ocimene, with an herb aroma description, were among the most important discriminant volatiles. This is in accordance with several studies where terpene hydrocarbons were reported as suitable markers of the geographic origin and genotype of EVOOs [29,47,48]. (E)-2-Dodecene has been reported as a discriminant variable in a study investigating VOOs of Chemlali and Neb Jmel cultivars growing in central and south Tunisia [49]. In this study, 1-ethyl-2-methylbenzene, 1-undecene, and (Z)-11-hexadecen-1-ol are reported as discriminant volatiles for geographical origin classification for the first time. In general, pentanal and 1-undecene, were found in increased amounts in oils from Chania, while the other statistically significant compounds were found in decreased levels in oils from Lasithi.
The aroma profile of the monovarietal Koroneiki VOOs examined was mainly characterized by the presence of C6 LOX alcohols, aldehydes, and esters as well as the contribution of C5 LOX compounds similarly to literature reports for VOOs of the same cultivar [6,50]. VOCs concentrations in the monovarietal Koroneiki VOOs of the current study varied significantly either as total content, or as total group content, or even as content of most individual VOCs. Differences were considered significant even within the same farming type (organic or conventional) or production region, revealing the complexity of the VOO aroma profile. The latter is in accordance with reported findings on variable profiles expected due to the effect and synergism of various biotic and abiotic parameters, such as agricultural practices, altitude, harvest time, and conditions employed from harvest-to-oil production and storage [51]. Nonetheless, VOCs identified in the current study as positive contributors (e.g., (2E)-hexenal, hexanal, hexan-1-ol, 1-penten-3-one, hexyl acetate, and hexenyl acetate), in terms of amount and/or importance of contribution to overall profile (relevant to OTV), have been also recognized as such for monovarietal Koroneiki VOOs studied [6,50]. Moreover, terpenic compounds useful indicators for geographical and botanical VOO differentiation have been reported to be present in the aroma profile of Koroneiki VOOs [6,50], which is in line with current findings; limonene, copaene, and a farnesene were present in all monovarietal Koroneiki VOO samples studied.
The application potential of the method could be improved by analysing a larger dataset of EVOO samples and/or samples from oncoming harvest years from the same areas. We are planning to expand our data set in this direction, analysing samples produced in subsequent harvest years. In addition, we plan to use this technology to compare Cretan vs. non-Cretan VOOs, a comparison that could be of legislating and commercial significance.

3. Materials and Methods

3.1. Samples

A total of 63 bottled and branded VOO samples that were previously evaluated for quality parameters were collected for the 5th Cretan olive oil competition. The samples were produced from olives of both organic (n = 19) and conventional (n = 44) olive groves and different local olive cultivars (Koroneiki and Tsounati) harvested in 2018–2019 in different regions of Crete (Chania n = 15, Rethymnon n = 12, Lasithi n = 10, and Heraklion n = 26) (Table S2). The organoleptic assessment of the olive oil samples was performed according to the Official EU method EEC/2568/91 [4] and its amendments by the trained sensory panel of the accredited Sensory Evaluation Laboratory of Crete ACR in Rethymnon Crete, Greece, and stored at −20 °C in dark glass vials that were shielded from light and oxygen.

3.2. Chemicals

LC-MS grade methanol was purchased from Sigma-Aldrich (Darmstadt, Germany). The internal standard eucalyptol (1,8-cineole) was of analytical grade and was purchased from Sigma-Aldrich (Darmstadt, Germany).

3.3. Optimisation of SPME Extraction Conditions

Initially, the effect of 3 experimental conditions for the HS-SPME method was studied by triplicate analysis following the well-established one-factor-at-a-time approach (OFAT). Sample mass was studied in the range of (1.000–4.000 g) in order to investigate any fiber saturation effects. Additionally, a range of settings was investigated for stirring of the sample (100–300 rpm). The latter two were optimised using the total peak area. The effects of extraction temperature for VOCs that are both desirable and undesirable were screened in the range of 30 to 85 °C so that an initial assessment was made for this crucial experimental parameter. Regarding the optimisation of the rest SPME conditions, a modelling approach was performed using the DOE and RSM approaches for the total peak area of each chromatogram representing all the ionisable and detectable volatile components of the sample. The optimisation was performed using a base fractional factorial (FF) 25−1 design with 5 factors (resolution V) that was augmented to a circumscribed central composite (CCC) with an alpha value of 2. A total of 32 runs included 16 base runs for FF (cube portion), 6 centre points to measure the experimental error, and 10 runs for CCC (star portion). Runs were randomised for each portion and were run in a single analytical sequence. The SPME factors selected for optimisation were extraction temperature, conditioning time, extraction time, desorption temperature, and desorption time. Although a preliminary screened factor, extraction temperature was included in the model as a significant factor to have a more complete picture of the SPME condition effects and to increase the model’s applicability. The experimental area for the corresponding factors is given in Table 3, with the corresponding design provided in Table S3. Factor and level selection was made based upon instrumental limitations, fit-for-purpose needs, and by referring to relevant literature [18,52].
The response selected for optimisation was the total chromatogram area, as measured by the sum of the areas of each ion, after the removal of the ions that were attributed mainly to siloxane contamination, as described in Section 3.6.

3.4. Sample Preparation

The optimal SPME conditions were set following the approach described in Section 3.3. Briefly, 1.85 gr of olive oil were weighed into a 15 mL glass vial, and 4 μL of a 1000 μg/mL eucalyptol internal standard solution in methanol were added, reaching a final concentration of 2 μg/mL. The vials were closed with a PTFE/silicone septum and equilibration was performed with a PAL Shimadzu autosampler unit (AOC 6000, CTC Analytics, Zwingen, Switzerland) at 55 ± 0.1 °C for 15 min upon agitation at 250 rpm. The (DVB/CAR/PDMS) fiber (2 cm length, 50/30 thickness) (Sigma-Aldrich) was then introduced to the headspace for 50 min at 55 °C. The fiber was pre- and post-conditioned for 10 min at 260 °C, according to the manufacturer’s instructions.

3.5. Gas Chromatography-Mass Spectrometry

VOCs were analyzed using a Shimadzu GCMS-QP2020 instrument equipped with a PAL SHIMADZU autosampler unit (AOC 6000, CTC Analytics, Zwingen, Switzerland). Chromatographic separation was performed on a MEGA-5 MS capillary column (30 m × 0.25 mm, 0.25 μm) (MEGA, Legano, Italy). Injection was operated in split mode (1:2) at 250 °C; He as the carrier gas at a flow rate of 1.2 mL/min; the column was held for 2.5 min at 40 °C, then programmed at 10 °C min−1 to 230 °C and held for 5 min. The mass detector was operated in the electron impact mode at 70 eV. Temperatures of MS source and quadrupole were set to 240 and 200 °C, respectively. The mass spectra scanned at m/z 40–450 amu range.
A pooled sample (QC) was prepared as a representative sample by mixing equal volumes of each VOO sample. The QC sample was injected at the start, every seven samples, and the end of the run to evaluate instrument stability. Blank runs were performed during the study to reveal possible carryover. The samples were analysed in randomised order.

3.6. Data Processing and Chemometrics

The Shimadzu software LabSolutions, GCMS Solutions version 2.50 SU3LabSolution, was used for acquiring the GC-MS data. Assignment Validator and Integrator (GAVIN) script for MATlab was used for performing peak integration complementary to Automated Mass spectral Deconvolution and Identification System (AMDIS) for peak deconvolution and identification. Total peaks areas for RSM-related modeling were extracted from the raw chromatogram files using OpenChrom® Lablicate Edition version 1.4 (Lablicate GmbH, Hamburg, Germany) [53] including the area of all scanned ions minus those attributed to siloxane contamination from column: 207, 281, and 355 [54]. Compound identification was performed by comparing mass spectra of eluting compounds with those of commercial libraries (NIST17 and FFNSC3) and by calculating linear retention indices relative to a series of n-alkanes (C8-C24) (Sigma-Aldrich). The VOCs were semi-quantified by dividing the peak areas of the compounds of interest by the peak area of the internal standard (Eucalyptol) and multiplying this ratio by the initial concentration of the internal standard (expressed as mg kg−1). The peak areas were measured from the full scan chromatograph using total ion current (TIC).
DOE design and analysis were performed using R version 3.6.2 [55] with RStudio (RStudio, Inc., Boston, MA, USA) [56] and DOE/RSM-related functionality [57,58,59]. The SIMCA package (version 13.0.2.; Umetrics, Sweden) was used for multivariate statistical analysis and biomarker assessment via VIP plots (Variable Importance for the Projection), loading plots, S-plots, p (corr), and hotelling’s lines. Principal components analysis (PCA) and orthogonal projection to latent structures discriminant analysis (OPLS-DA) were performed to assess data in a multivariate setting. Models’ validation was evaluated using permutation plots and CV-ANOVA value. Two-tailed t-test, with unequal variance and a threshold of p < 0.05, and ANOVA were employed in Microsoft Excel Spreadsheets.

4. Conclusions

In this study, a HS-SPME/GC-MS method by Design of Experiments (DoE) and Response Surface Methodology (RSM) approaches was successfully optimized. The optimized parameters were sample amount, sample stirring, extraction temperature, conditioning time, extraction time, desorption temperature, and desorption time. The optimized HS-SPME/GC-MS method allowed the identification and semi-quantification of a high number of compounds. Of the 92 VOCs identified, 4 are found for the first time in EVOO while others, mainly terpenoid compounds, are rarely identified in EVOO, with the use of SPME/GC-MS. Statistical analysis allowed for identifying markers with significance in the geographical classification of Cretan EVOO samples. The method is deemed useful for the analysis of EVOO, providing easy, automated, efficient, and economic sample preparation and analysis along with information-rich GC-MS data.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/metabo12020114/s1, Table S1: Quadratic model parameters including term estimates, standard error (Std.Error), t values and p-values. Term significance is provided in the column Signif, in a scale given at the footnote of the table. Extraction temperature (x1), conditioning time (x2), extraction time (x3), desorption time (x4) and desorption temperature (x5). Example: x1 is a first order term, x1:x2 is a two-way interaction term and x1^2 is a quadratic term, Table S2: Quality characteristics of EVOO samples, Table S3: Complete Circumscribed Central Composite design indicating the run order along with the levels of each factor. Units for extraction temperature (Ex.temp) are °C, conditioning time (Co.time) mins, extraction time (Ex.time) mins, desorption temperature (Des.temp) °C and desorption time (Des.time) mins, Figure S1: Contour plot for factors desorption temperature (Des.temp) and desorption time (Des.time). Green colour indicates lower response values and pink colour higher response values. Contours include corresponding response levels. The contour plot is a slice at the selected points for extraction temperature (Ex.temp) and extraction time (Ex.time), and the centre point for conditioning time (Co.time), Figure S2: PCA score plot showing the examined groups together with QC samples. QC samples are clustering together, showing system’s stability, Figure S3: Geographical origin of Cretan EVOOs analysed in our study.

Author Contributions

A.L., I.S. and C.V., methodology, formal analysis, data curation, and writing—original draft preparation; V.T.P., A.S. and K.G.Z., Formal analysis, Writing—original draft; G.T., data curation, writing—original draft preparation, writing—review and editing, visualization, supervision and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the project FoodOmicsGR “Comprehensive Characterization of Foods” (MIS 5029057), which is implemented under the action “Reinforcement of the Research and Innovation Infrastructure”, funded by the Operational Programme “Competitiveness, Entrepreneurship and Innovation” (NSRF 2014–2020) and co-financed by Greece and the European Union (European Regional Development Fund). All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful to the Sensory Evaluation Laboratory of Crete ACR in Rethymnon Crete, Greece for providing the organoleptic data of the VOO samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arroyo-Manzanares, N.; Gabriel, F.; Carpio, A.; Arce, L. Use of Whole Electrophoretic Profile and Chemometric Tools for the Differentiation of Three Olive Oil Qualities. Talanta 2019, 197, 175–180. [Google Scholar] [CrossRef] [PubMed]
  2. Conte, L.; Bendini, A.; Valli, E.; Lucci, P.; Moret, S.; Maquet, A.; Lacoste, F.; Brereton, P.; García-González, D.L.; Moreda, W.; et al. Olive Oil Quality and Authenticity: A Review of Current EU Legislation, Standards, Relevant Methods of Analyses, Their Drawbacks and Recommendations for the Future. Trends Food Sci. Technol. 2020, 105, 483–493. [Google Scholar] [CrossRef]
  3. Da Silva, M.D.; Freitas, A.M.; Cabrita, M.J.; Garcia, R. Olive Oil Composition: Volatile Compounds; IntechOpen: London, UK, 2012; ISBN 978-953-307-921-9. [Google Scholar]
  4. Commission Regulation (EEC) No 2568/91 of 11 July 1991 on the Characteristics of Olive Oil and Olive-Residue Oil and on the Relevant Methods of Analysis; European Union: Luxembourg, 1991; Volume 248.
  5. Kalua, C.M.; Allen, M.S.; Bedgood, D.R.; Bishop, A.G.; Prenzler, P.D.; Robards, K. Olive Oil Volatile Compounds, Flavour Development and Quality: A Critical Review. Food Chem. 2007, 100, 273–286. [Google Scholar] [CrossRef]
  6. Theodosi, S.; Kosma, I.S.; Badeka, A.V. Quality Characteristics of Koroneiki Olive Oil from Zakynthos Island (Greece) and Differentiation Depending on the Altitude Level. Eur. Food Res. Technol. 2021, 247, 1235–1248. [Google Scholar] [CrossRef]
  7. Skiada, V.; Tsarouhas, P.; Varzakas, T. Comparison and Discrimination of Two Major Monocultivar Extra Virgin Olive Oils in the Southern Region of Peloponnese, According to Specific Compositional/Traceability Markers. Foods 2020, 9, 155. [Google Scholar] [CrossRef] [Green Version]
  8. Revelou, P.-K.; Pappa, C.; Kakouri, E.; Kanakis, C.D.; Papadopoulos, G.K.; Pappas, C.S.; Tarantilis, P.A. Discrimination of Botanical Origin of Olive Oil from Selected Greek Cultivars by SPME-GC-MS and ATR-FTIR Spectroscopy Combined with Chemometrics. J. Sci. Food Agric. 2021, 101, 2994–3002. [Google Scholar] [CrossRef] [PubMed]
  9. Pousinis, P.; Pechlivanis, A.; Lioupi, A.; Gika, H.; Virgiliou, C.; Kodra, D.; Marinaki, M.; Theodoridis, G. Innovation in Nutrition and Food Analysis: The FoodOmicsGR_Research Infrastructure. Column 2021, 17, 12–19. [Google Scholar]
  10. Theodoridis, G.; Pechlivanis, A.; Thomaidis, N.S.; Spyros, A.; Georgiou, C.A.; Albanis, T.; Skoufos, I.; Kalogiannis, S.; Tsangaris, G.T.; Stasinakis, A.S.; et al. FoodOmicsGR_RI: A Consortium for Comprehensive Molecular Characterisation of Food Products. Metabolites 2021, 11, 74. [Google Scholar] [CrossRef]
  11. Lioupi, A.; Nenadis, N.; Theodoridis, G. Virgin Olive Oil Metabolomics: A Review. J. Chromatogr. B 2020, 1150, 122161. [Google Scholar] [CrossRef]
  12. Castro-Puyana, M.; Pérez-Míguez, R.; Montero, L.; Herrero, M. Reprint of: Application of Mass Spectrometry-Based Metabolomics Approaches for Food Safety, Quality and Traceability. TrAC Trends Anal. Chem. 2017, 96, 62–78. [Google Scholar] [CrossRef]
  13. Reboredo-Rodríguez, P.; González-Barreiro, C.; Cancho-Grande, B.; Simal-Gándara, J. Dynamic Headspace/GC–MS to Control the Aroma Fingerprint of Extra-Virgin Olive Oil from the Same and Different Olive Varieties. Food Control. 2012, 25, 684–695. [Google Scholar] [CrossRef]
  14. Aprea, E.; Gasperi, F.; Betta, E.; Sani, G.; Cantini, C. Variability in Volatile Compounds from Lipoxygenase Pathway in Extra Virgin Olive Oils from Tuscan Olive Germoplasm by Quantitative SPME/GC-MS. J. Mass Spectrom. 2018, 53, 824–832. [Google Scholar] [CrossRef]
  15. Romero, I.; García-González, D.L.; Aparicio-Ruiz, R.; Morales, M.T. Validation of SPME-GCMS Method for the Analysis of Virgin Olive Oil Volatiles Responsible for Sensory Defects. Talanta 2015, 134, 394–401. [Google Scholar] [CrossRef]
  16. Aprea, E.; Gika, H.; Carlin, S.; Theodoridis, G.; Vrhovsek, U.; Mattivi, F. Metabolite Profiling on Apple Volatile Content Based on Solid Phase Microextraction and Gas-Chromatography Time of Flight Mass Spectrometry. J. Chromatogr. A 2011, 1218, 4517–4524. [Google Scholar] [CrossRef]
  17. Borges, T.H.; Ramalhosa, E.; Seiquer, I.; Pereira, J.A. Use of Response Surface Methodology (RSM) for the Identification of the Best Extraction Conditions for Headspace Solid-Phase Micro Extraction (HS-SPME) of the Volatile Profile of Cv. Arbequina Extra-Virgin Olive Oil. Eur. J. Lipid Sci. Technol. 2018, 120, 1700356. [Google Scholar] [CrossRef]
  18. De los Angeles Fernandez, M.; Assof, M.; Jofre, V.; Silva, M.F. Volatile Profile Characterization of Extra Virgin Olive Oils from Argentina by HS-SPME/GC-MS and Multivariate Pattern Recognition Tools. Food Anal. Methods 2014, 7, 2122–2136. [Google Scholar] [CrossRef]
  19. Wardencki, W.; Michulec, M.; Curyło, J. A Review of Theoretical and Practical Aspects of Solid-Phase Microextraction in Food Analysis. Int. J. Food Sci. Technol. 2004, 39, 703–717. [Google Scholar] [CrossRef]
  20. Kiritsakis, A.K. Flavor Components of Olive Oil—A Review. J. Am. Oil Chem. Soc. 1998, 75, 673–681. [Google Scholar] [CrossRef]
  21. Nunes, C.A.; Souza, V.R.; de Corrêa, S.C.; Silva, M.; de Cássia da Costa e Silva, M.; Bastos, S.C.; Pinheiro, A.C.M. Heating on the Volatile Composition and Sensory Aspects of Extra-Virgin Olive Oil. Ciênc. Agrotec. 2013, 37, 566–572. [Google Scholar] [CrossRef] [Green Version]
  22. Morales, M.T.; Rios, J.J.; Aparicio, R. Changes in the Volatile Composition of Virgin Olive Oil during Oxidation:  Flavors and Off-Flavors. J. Agric. Food Chem. 1997, 45, 2666–2673. [Google Scholar] [CrossRef]
  23. Morales, M.T.; Przybylski, R. Olive Oil Oxidation. In Handbook of Olive Oil: Analysis and Properties; Aparicio, R., Harwood, J., Eds.; Springer US: Boston, MA, USA, 2013; pp. 479–522. ISBN 978-1-4614-7777-8. [Google Scholar]
  24. Aparicio, R.; Morales, M.T. Characterization of Olive Ripeness by Green Aroma Compounds of Virgin Olive Oil. J. Agric. Food Chem. 1998, 46, 1116–1122. [Google Scholar] [CrossRef]
  25. Dhifi, W.; Angerosa, F.; Serraiocco, A.; Oumar, I.; Hamrouni, I.; Marzouk, B. Virgin Olive Oil Aroma: Characterization of Some Tunisian Cultivars. Food Chem. 2005, 93, 697–701. [Google Scholar] [CrossRef]
  26. Angerosa, F. Influence of Volatile Compounds on Virgin Olive Oil Quality Evaluated by Analytical Approaches and Sensor Panels. Eur. J. Lipid Sci. Technol. 2002, 104, 639–660. [Google Scholar] [CrossRef]
  27. Petronilho, S.; Lopez, R.; Ferreira, V.; Coimbra, M.A.; Rocha, S.M. Revealing the Usefulness of Aroma Networks to Explain Wine Aroma Properties: A Case Study of Portuguese Wines. Molecules 2020, 25, 272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Brkić Bubola, K.; Koprivnjak, O.; Sladonja, B.; Lukić, I. Volatile Compounds and Sensory Profiles of Monovarietal Virgin Olive Oils from Buža, Črna and Rosinjola Cultivars in Istria (Croatia). Food Technol. Biotechnol. 2012, 50, 192–198. [Google Scholar]
  29. Cecchi, T.; Alfei, B. Volatile Profiles of Italian Monovarietal Extra Virgin Olive Oils via HS-SPME–GC–MS: Newly Identified Compounds, Flavors Molecular Markers, and Terpenic Profile. Food Chem. 2013, 141, 2025–2035. [Google Scholar] [CrossRef]
  30. Mansouri, F.; Moumen, A.B.; Richard, G.; Fauconnier, M.-L.; Sindic, M.; Caid, H.S.; Elamrani, A. Flavor Profiles of Monovarietal Virgin Olive Oils Produced in the Oriental Region of Morocco. Oilseeds Fats Crops Lipids 2017, 24, A501. [Google Scholar] [CrossRef] [Green Version]
  31. Kosma, I.S.; Kontominas, M.G.; Badeka, A.V. The Application of Chemometrics to Volatile Compound Analysis for the Recognition of Specific Markers for Cultivar Differentiation of Greek Virgin Olive Oil Samples. Foods 2020, 9, 1672. [Google Scholar] [CrossRef]
  32. Angerosa, F.; Mostallino, R.; Basti, C.; Vito, R. Virgin Olive Oil Odour Notes: Their Relationships with Volatile Compounds from the Lipoxygenase Pathway and Secoiridoid Compounds. Food Chem. 2000, 68, 283–287. [Google Scholar] [CrossRef]
  33. Reboredo-Rodríguez, P.; González-Barreiro, C.; Cancho-Grande, B.; Simal-Gándara, J. Improvements in the Malaxation Process to Enhance the Aroma Quality of Extra Virgin Olive Oils. Food Chem. 2014, 158, 534–545. [Google Scholar] [CrossRef]
  34. Luna, G.; Morales, M.T.; Aparicio, R. Characterisation of 39 Varietal Virgin Olive Oils by Their Volatile Compositions. Food Chem. 2006, 98, 243–252. [Google Scholar] [CrossRef]
  35. Pérez, A.G.; de la Rosa, R.; Pascual, M.; Sánchez-Ortiz, A.; Romero-Segura, C.; León, L.; Sanz, C. Assessment of Volatile Compound Profiles and the Deduced Sensory Significance of Virgin Olive Oils from the Progeny of Picual × Arbequina Cultivars. J. Chromatogr. A 2016, 1428, 305–315. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Stilo, F.; Cordero, C.; Sgorbini, B.; Bicchi, C.; Liberto, E. Highly Informative Fingerprinting of Extra-Virgin Olive Oil Volatiles: The Role of High Concentration-Capacity Sampling in Combination with Comprehensive Two-Dimensional Gas Chromatography. Separations 2019, 6, 34. [Google Scholar] [CrossRef] [Green Version]
  37. Youssef, O.; Guido, F.; Manel, I.; Youssef, N.B.; Luigi, C.P.; Mohamed, H.; Daoud, D.; Mokhtar, Z. Volatile Compounds and Compositional Quality of Virgin Olive Oil from Oueslati Variety: Influence of Geographical Origin. Food Chem. 2011, 124, 1770–1776. [Google Scholar] [CrossRef]
  38. Cerretani, L.; Salvador, M.D.; Bendini, A.; Fregapane, G. Relationship Between Sensory Evaluation Performed by Italian and Spanish Official Panels and Volatile and Phenolic Profiles of Virgin Olive Oils. Chem. Percept. 2008, 1, 258. [Google Scholar] [CrossRef]
  39. Gilbert-López, B.; Robles-Molina, J.; García-Reyes, J.F.; Molina-Díaz, A. Rapid Determination of BTEXS in Olives and Olive Oil by Headspace-Gas Chromatography/Mass Spectrometry (HS-GC–MS). Talanta 2010, 83, 391–399. [Google Scholar] [CrossRef]
  40. Bempelou, E.; Anagnostopoulos, C.; Liapis, K. Investigation of Naphthalene Contamination in Olive Oil from Greece. Toxicol. Environ. Chem. 2019, 101, 45–58. [Google Scholar] [CrossRef]
  41. Nenadis, N.; Mastralexi, A.; Tsimidou, M.Z. Physicochemical Characteristics and Antioxidant Potential of the Greek PDO and PGI Virgin Olive Oils (VOOs). Eur. J. Lipid Sci. Technol. 2019, 121, 1800172. [Google Scholar] [CrossRef]
  42. Kosma, I.; Vatavali, K.; Kontakos, S.; Kontominas, M.; Kiritsakis, A.; Badeka, A. Geographical Differentiation of Greek Extra Virgin Olive Oil from Late-Harvested Koroneiki Cultivar Fruits. J. Am. Oil Chem Soc. 2017, 94, 1373–1384. [Google Scholar] [CrossRef]
  43. Pouliarekou, E.; Badeka, A.; Tasioula-Margari, M.; Kontakos, S.; Longobardi, F.; Kontominas, M.G. Characterization and Classification of Western Greek Olive Oils According to Cultivar and Geographical Origin Based on Volatile Compounds. J. Chromatogr. A 2011, 1218, 7534–7542. [Google Scholar] [CrossRef]
  44. Lukić, I.; Lukić, M.; Žanetić, M.; Krapac, M.; Godena, S.; Brkić Bubola, K. Inter-Varietal Diversity of Typical Volatile and Phenolic Profiles of Croatian Extra Virgin Olive Oils as Revealed by GC-IT-MS and UPLC-DAD Analysis. Foods 2019, 8, 565. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Bajoub, A.; Sánchez-Ortiz, A.; Ajal, E.A.; Ouazzani, N.; Fernández-Gutiérrez, A.; Beltrán, G.; Carrasco-Pancorbo, A. First Comprehensive Characterization of Volatile Profile of North Moroccan Olive Oils: A Geographic Discriminant Approach. Food Res. Int. 2015, 76, 410–417. [Google Scholar] [CrossRef] [PubMed]
  46. Şişik Oğraş, Ş.; Kaban, G.; Kaya, M. Volatile Compounds of Olive Oils from Different Geographic Regions in Turkey. Int. J. Food Prop. 2018, 21, 1833–1843. [Google Scholar] [CrossRef] [Green Version]
  47. Vichi, S.; Guadayol, J.M.; Caixach, J.; López-Tamames, E.; Buxaderas, S. Monoterpene and Sesquiterpene Hydrocarbons of Virgin Olive Oil by Headspace Solid-Phase Microextraction Coupled to Gas Chromatography/Mass Spectrometry. J. Chromatogr. A 2006, 1125, 117–123. [Google Scholar] [CrossRef] [PubMed]
  48. Zunin, P.; Boggia, R.; Salvadeo, P.; Evangelisti, F. Geographical Traceability of West Liguria Extravirgin Olive Oils by the Analysis of Volatile Terpenoid Hydrocarbons. J. Chromatogr. A 2005, 1089, 243–249. [Google Scholar] [CrossRef]
  49. Mansour, A.B.; Gargouri, B.; Flamini, G.; Bouaziz, M. Effect of Agricultural Sites on Differentiation between Chemlali and Neb Jmel Olive Oils. J. Oleo. Sci. 2015, 64, 381–392. [Google Scholar] [CrossRef] [Green Version]
  50. Kandylis, P.; Vekiari, A.S.; Kanellaki, M.; Grati Kamoun, N.; Msallem, M.; Kourkoutas, Y. Comparative Study of Extra Virgin Olive Oil Flavor Profile of Koroneiki Variety (Olea Europaea Var. Microcarpa Alba) Cultivated in Greece and Tunisia during One Period of Harvesting. LWT Food Sci. Technol. 2011, 44, 1333–1341. [Google Scholar] [CrossRef]
  51. Angerosa, F.; Servili, M.; Selvaggini, R.; Taticchi, A.; Esposto, S.; Montedoro, G. Volatile Compounds in Virgin Olive Oil: Occurrence and Their Relationship with the Quality. J. Chromatogr. A 2004, 1054, 17–31. [Google Scholar] [CrossRef]
  52. Ribeiro, L.H.; Costa Freitas, A.M.; Gomes da Silva, M.D.R. The Use of Headspace Solid Phase Microextraction for the Characterization of Volatile Compounds in Olive Oil Matrices. Talanta 2008, 77, 110–117. [Google Scholar] [CrossRef]
  53. Wenig, P.; Odermatt, J. OpenChrom: A Cross-Platform Open Source Software for the Mass Spectrometric Analysis of Chromatographic Data. BMC Bioinform. 2010, 11, 405. [Google Scholar] [CrossRef] [Green Version]
  54. Brocks, J.J.; Hope, J.M. Tailing of Chromatographic Peaks in GC–MS Caused by Interaction of Halogenated Solvents with the Ion Source. J. Chromatogr. Sci. 2014, 52, 471–475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ihaka, R.; Gentleman, R. R: A Language for Data Analysis and Graphics. J. Comput. Graph. Stat. 1996, 5, 299–314. [Google Scholar] [CrossRef]
  56. RStudio|Open Source & Professional Software for Data Science Teams. Available online: https://rstudio.com/ (accessed on 17 August 2021).
  57. Grömping, U. R Package FrF2 for Creating and Analyzing Fractional Factorial 2-Level Designs. J. Stat. Softw. 2014, 56, 1–56. [Google Scholar] [CrossRef] [Green Version]
  58. Grömping, U. R Package DoE.Base for Factorial Experiments. J. Stat. Softw. 2018, 85, 1–41. [Google Scholar] [CrossRef] [Green Version]
  59. Lenth, R.V. Response-Surface Methods in R., Using Rsm. J. Stat. Softw. 2009, 32, 1–17. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Influence of (a) sample amount, (b) stirring speed on total areas, and (c) extraction temperature on selected VOCs. (E + 08 = × 108).
Figure 1. Influence of (a) sample amount, (b) stirring speed on total areas, and (c) extraction temperature on selected VOCs. (E + 08 = × 108).
Metabolites 12 00114 g001aMetabolites 12 00114 g001b
Figure 2. Contour plot for factors extraction temperature (Ex.temp) and extraction time (Ex.time). Green colour indicates lower response values and pink colour higher response values. Contours include corresponding response levels. The contour plot is a slice at the centre points for the rest of the remaining factors (Conditioning time (Co.time), Desorption temperature (Des.temp), and Desorption time (Des.time)). (e + 08 = × 108).
Figure 2. Contour plot for factors extraction temperature (Ex.temp) and extraction time (Ex.time). Green colour indicates lower response values and pink colour higher response values. Contours include corresponding response levels. The contour plot is a slice at the centre points for the rest of the remaining factors (Conditioning time (Co.time), Desorption temperature (Des.temp), and Desorption time (Des.time)). (e + 08 = × 108).
Metabolites 12 00114 g002
Figure 3. Major VOCs identified in the TIC of a QC olive oil sample analysed with optimised HS-SPME/GC-MS conditions. Upper in the insets, zoomed areas are shown. The identified components are given in numbers. (1) 1-methoxy-2-propanol, (2) 1-penten-3-ol, (3) 1-penten-3-one (4) 3-pentanone, (5) pentanal, (6) 3-methylbutan-1-ol, (7) (E)-2-pentenal, (8) (E)-2-pentenal, (9) 1-pentanol, (10) 2-penten-1-ol, (Z)-, (11) 1-octene, (12) octane (13) 3-hexenal, (14) hexanal, (15) (Z)-hex-3-en-1-ol, (16) (2E)-hexenal, (17) (E)-2-hexen-1-ol, (18) 1-hexanol, (19) (prop-2-en-1-yl)cyclopentane, (20) o-xylene, (21) 3-ethyl-1,5-octadiene, (22) heptanal, (23) (Z)- 2-penten-1-ol, acetate, (24) 2,4-hexadienal, (25) 3-ethyl-1,5-octadiene, (26) 3-ethyl-1,5-octadiene, (27) 2,2-dimethyl-3-heptanone, (28) 5-ethyl-2(5H)-furanone, (29) hexanoic acid, (30) 6-methyl-5-hepten-2-one, (31) 3-ethyl-1,5-octadiene, (32) 3-ethyl-1,5-octadiene, (33) (3Z)-hex-3-en-1-yl acetate, (34) hexyl acetate, (35) eucalyptol (I.S), (36) benzyl alcohol, (37) o-cymene, (38) 2,2-dimethyl-3-heptanone, (39) 1-undecene, (40) nonanal, (41) (E)-4,8-dimethyl-1,3,7-nonatriene, (42) 2-ethylhexanoic acid, (43) octanoic acid, (44) 1-nonanol, (45) (E)-2-dodecene, (46) 11-hexadecen-1-ol, (Z)-, (47) nonanoic acid, (48) 1-dodecene, (49) copaene, (50) pentadecane, (51) eremophilene, (52) alpha-farnesene. (E + 07 = × 107).
Figure 3. Major VOCs identified in the TIC of a QC olive oil sample analysed with optimised HS-SPME/GC-MS conditions. Upper in the insets, zoomed areas are shown. The identified components are given in numbers. (1) 1-methoxy-2-propanol, (2) 1-penten-3-ol, (3) 1-penten-3-one (4) 3-pentanone, (5) pentanal, (6) 3-methylbutan-1-ol, (7) (E)-2-pentenal, (8) (E)-2-pentenal, (9) 1-pentanol, (10) 2-penten-1-ol, (Z)-, (11) 1-octene, (12) octane (13) 3-hexenal, (14) hexanal, (15) (Z)-hex-3-en-1-ol, (16) (2E)-hexenal, (17) (E)-2-hexen-1-ol, (18) 1-hexanol, (19) (prop-2-en-1-yl)cyclopentane, (20) o-xylene, (21) 3-ethyl-1,5-octadiene, (22) heptanal, (23) (Z)- 2-penten-1-ol, acetate, (24) 2,4-hexadienal, (25) 3-ethyl-1,5-octadiene, (26) 3-ethyl-1,5-octadiene, (27) 2,2-dimethyl-3-heptanone, (28) 5-ethyl-2(5H)-furanone, (29) hexanoic acid, (30) 6-methyl-5-hepten-2-one, (31) 3-ethyl-1,5-octadiene, (32) 3-ethyl-1,5-octadiene, (33) (3Z)-hex-3-en-1-yl acetate, (34) hexyl acetate, (35) eucalyptol (I.S), (36) benzyl alcohol, (37) o-cymene, (38) 2,2-dimethyl-3-heptanone, (39) 1-undecene, (40) nonanal, (41) (E)-4,8-dimethyl-1,3,7-nonatriene, (42) 2-ethylhexanoic acid, (43) octanoic acid, (44) 1-nonanol, (45) (E)-2-dodecene, (46) 11-hexadecen-1-ol, (Z)-, (47) nonanoic acid, (48) 1-dodecene, (49) copaene, (50) pentadecane, (51) eremophilene, (52) alpha-farnesene. (E + 07 = × 107).
Metabolites 12 00114 g003
Figure 4. Natural variation of VOCs derived from LOX pathway in a (a) Koroneiki and (b) Tsounati variety sample.
Figure 4. Natural variation of VOCs derived from LOX pathway in a (a) Koroneiki and (b) Tsounati variety sample.
Metabolites 12 00114 g004
Figure 5. Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) on EVOO samples volatile profile, using the ‘geographical origin’ as class membership criterion, (a) Lasithi vs. Rethymno, (b) Chania vs. Lasithi.
Figure 5. Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) on EVOO samples volatile profile, using the ‘geographical origin’ as class membership criterion, (a) Lasithi vs. Rethymno, (b) Chania vs. Lasithi.
Metabolites 12 00114 g005
Figure 6. Box plots of identified biomarkers in VOOs samples. Boxes are drawn from the 25th to 75th percentiles in the concentration distribution. (a) Lasithi vs. Rethymnon, (b) Chania vs. Lasithi. (e + 07 = × 107).
Figure 6. Box plots of identified biomarkers in VOOs samples. Boxes are drawn from the 25th to 75th percentiles in the concentration distribution. (a) Lasithi vs. Rethymnon, (b) Chania vs. Lasithi. (e + 07 = × 107).
Metabolites 12 00114 g006
Table 1. Model fit data along with corresponding output from the ANOVA analysis table.
Table 1. Model fit data along with corresponding output from the ANOVA analysis table.
Multiple R-Squared0.9028
Adjusted R-Squared0.7262
p-Value4.02 × 10−3
DfSum SqMean SqF ValuePr (>F)
FO (x1, x2, x3, x4, x5)51.2055 × 10182.4111 × 101716.37489.02 × 10−5
TWI (x1, x2, x3, x4, x5)101.9382 × 10171.9382 × 10161.31640.3283
PQ (x1, x2, x3, x4, x5)51.0584 × 10172.1168 × 10161.43760.2855
Residuals111.6197 × 10171.4724 × 1016
Lack of Fit64.5182 × 10167.5303 × 10150.32240.8996
Pure Error51.1679 × 10172.3357 × 1016
FO stands for First Order terms, TWI for Two Way Interaction terms and PQ for Pure Quadratic terms. Df are the Degrees of Freedom, Sum Sq the Sum of Squares and Mean Sq the Mean Square. Extraction temperature (x1), conditioning time (x2), extraction time (x3), desorption time (x4) and desorption temperature (x5).
Table 2. Chromatographic data (retention time, retention index) and mean values (mg/kg) relative to internal standard (Eucalyptol) and standard deviation for identified volatile compounds in Cretan EVOO samples.
Table 2. Chromatographic data (retention time, retention index) and mean values (mg/kg) relative to internal standard (Eucalyptol) and standard deviation for identified volatile compounds in Cretan EVOO samples.
Compound NameRtRIexpRI LitChania
Mean (n = 17)
±SDHeraklion
Mean (n = 26)
±SDLasithi Mean
(n = 10)
±SDRethymnon Mean
(n = 10)
±SD
Alcohols
Isobutanol26266220.060.050.060.050.10.040.070.06
1-Penten-3-ol, C52.546836861.960.381.690.491.450.491.050.37
3-Methylbutan-1-ol3.337387440.120.050.140.050.150.050.130.07
1-Pentanol3.897717660.090.080.080.060.040.050.170.05
(Z)-2-Penten-1-ol, C53.927747671.390.251.050.360.860.350.850.19
(Z)-Hex-3-en-1-ol, C6-LnA5.538578580.120.040.120.050.090.050.190.06
(E)-Hex-2-en-1-ol, C6-LnA5.828718661.991.752.232.490.722.367.963.45
Hexan-1-ol, C6-LA5.98758621.770.91.620.880.990.94.311.12
1-Octanol9.76107610740.080.020.080.030.070.020.070.04
1-Nonanol11.48117611720.060.030.10.050.10.040.050.04
(Z)-11-Hexadecen-1-ol12.311228-0.050.010.050.020.050.020.030.02
Aldehydes
3-Methylbutanal2.266526500.040.030.040.040.050.030.050.03
(E)-2-Pentenal, C53.527497540.080.030.060.030.060.020.040.02
(E)-2-Pentenal, C53.697607540.30.080.250.090.220.090.170.07
Pentanal2.757036950.990.440.870.40.750.390.370.22
3-Hexenal, C6-LnA4.488067952.852.781.391.850.831.661.840.87
Hexanal, C6-LA4.528087983.511.423.11.062.751.074.221.08
2-Hexenal, C6-LnA5.478548540.190.050.160.050.150.050.220.07
(2E)-Hexenal, C6-LnA5.6186285625.439.6222.169.5721.519.2332.6713.59
Heptanal6.69089030.190.070.210.070.170.060.120.06
2,4 Hexadienal6.799199161.950.711.430.611.220.571.660.34
Benzeneacetaldehyde9.39105510490.030.010.040.020.060.020.060.02
Nonanal10.38111110990.980.471.140.541.170.50.450.47
Acids
Acetic acid1.725965950.560.560.520.560.340.490.570.89
3-Methylbutanoic acid5.358528580.020.010.020.010.020.010.030.01
Hexanoic acid8.049819810.730.710.550.510.210.490.70.35
Heptanoic acid9.72107410760.10.030.090.030.060.030.090.03
2-Ethylhexanoic acid10.52111811230.580.550.370.4500.410.370.41
Octanoic acid11.36117011770.140.060.130.050.060.050.140.04
Nonanoic acid12.92126812700.250.10.250.10.10.090.260.09
Esters
Ethyl acetate1.896156140.220.120.270.130.290.120.220.16
3-Methylbutyl acetate6.048828830.130.080.130.070.180.0700.11
(Z)-Pent-2-en-1-yl acetate6.769189160.590.330.390.250.270.230.510.11
(3Z)-Hex-3-en-1-yl acetate, C6-LnA8.551009100518.499.0312.157.1913.467.227.974.09
Hexyl acetate, C6-LA8.68101610114.52.443.071.854.211.892.811.51
Ethers
1-Methoxy-2-propanol2.366546500.080.030.110.080.10.070.080.03
1-Methoxyhexane4.998308310000.0200.0600
(3Z)-1-Methoxy-3-hexene 5.05838832000000.0800
Ketones
1-Penten-3-one, C52.586876791.960.921.651.171.461.070.71.15
3-Pentanone2.717016941.650.921.260.721.040.751.510.77
2,2-Dimethyl-3-heptanone9.641070-5.075.212.073.540.853.143.221.13
Terpenic Compounds
.alpha.-Pinene79309320.010.040.020.0700.060.180.03
.beta.-Pinene8.059839800.070.260.070.3100.2700.15
6-Methyl-5-hepten-2-one8.29909850.610.410.560.330.510.310.780.34
Alpha-Terpinene8.79101810170000.0100.010.020
3-Carene8.791018101800.010.010.0100.010.020
(Z)-beta-Ocimene9.29105010480.050.020.070.030.070.030.090.02
o-Cymene8.95103110370.140.480.20.5900.541.330
d-Limonene9.02103510410.71.830.993.290.12.870.194.33
(+)-2-Carene9.0410351031000.060.2300.230.920
Neral12.611247124400.0100.0100.0100.01
Citral13.071276127600.0100.0200.0200.02
(+)-Cyclosativene14.6713841368000.010.0100.010.010.01
Copaene14.73138813920.030.010.050.030.060.030.10.02
Eremophilene16.34150515030.010.010.010.010.010.010.010.01
alpha-Farnesene16.39150915090.050.030.060.030.050.030.050.07
Kessane16.941551153000000000
Liguloxide17.0715611533000.010.0100.0100.01
Others
2-Methylpentane1.625885730.040.030.040.020.030.020.030.01
3-Methylpentane3.397417480.050.030.070.030.080.030.090.03
Toluene3.867707710.170.10.190.210.130.240.060.2
1-Octene4.257947920.170.10.150.080.140.070.210.05
Octane4.428038001.640.782.131.082.2813.860.79
Styrene6.368988990.040.040.040.0600.060.090.11
o-Xylene6.368988990.040.050.060.070.060.0700.1
(Prop-2-en-1-yl)cyclopentane6.318958980.310.060.270.070.230.080.220.06
Pentene dimer 1, C56.4899-0.220.050.190.060.160.060.150.05
Pentene dimer 2, C57.199409471.650.411.410.451.210.470.980.33
Pentene dimer 3, C57.339479491.90.461.540.531.350.531.220.34
4,8-Dimethyl-1,7-nonadiene7.389499980.010.0200.0100.0100
Propylbenzene7.5996096200000000.02
2,2,6-Trimethyloctane7.63961964000.010.0300.0200.15
2,2-Dimethyl-3-heptanone7.699659650.780.620.420.430.230.390.460.22
1-Ethyl-2-methylbenzene7.739719690.050.030.040.040.030.040.040.05
5-Ethyl-2(5H)-furanone7.8297196812.9112.725.398.72.227.766.53.27
2,2,4-Trimethylpentane8.2990-0.150.50.050.3400.3800.08
Pentene dimer 4, C58.27994-0.880.230.770.230.660.240.580.18
Pentene dimer 5, C58.34998-3.260.872.70.922.380.942.140.63
2-Propylfuran8.451004-0.120.050.10.040.080.040.130.04
4-Ethyltoluene8.891029-0.020.010.020.010.010.010.020.01
2,2,4,4,6-Pentamethylheptane8.931030-0000.0100.0400.07
5-Ethyl-2(5H)-furanone9.171044-0.640.660.270.440.120.390.310.14
1-Phenyl-1-propanone9.421059 0000.010000
1,4-Diethylbenzene9.53106310700000.0100.0100
1-Undecene10.07107810820.040.010.030.020.030.020.020.02
(E)-4,8-Dimethyl-1,3,7-nonatriene10.11109610980.060.020.090.040.090.040.050.03
Methyl benzoate10.27110411060.040.010.050.020.050.020.030.01
Naphthalene11.88120112060.060.2400.150.010.1300
(E)-2-Dodecene11.96120612061.780.322.140.642.020.671.050.58
1-Dodecene12.96127012650.040.020.040.020.030.020.070.02
Pentadecane16.28150115000.010.0100.0100.010.020.01
2,4-Di-tert-butylphenol16.44151315170.020.010.020.010.010.010.010.01
Total 108.1 81.42 92.67 85.77
C513.6 11.3 12.35 12.06
C654.37 42.93 47.86 42.74
C5 + C667.97 54.22 60.21 54.8
Rt: Retention Time, I.S.: Internal Standard, RIex: Retention Index experimental, RIlit: Retention Index literature, SD: Standard Deviation.
Table 3. Factors and their corresponding levels in each portion of the RSM design: cube portion for the FF base design and star portion for the CCC augmentation.
Table 3. Factors and their corresponding levels in each portion of the RSM design: cube portion for the FF base design and star portion for the CCC augmentation.
Abbr.Cube Level 1Cube Level 2Star Level 1Star Level 2Units
Ex.temp40702585°C
Co.time1020525mins
Ex.time2050565mins
Des.temp250270240280°C
Des.time5102.512.5mins
Ex.temp stands for extraction temperature, Co.time for conditioning time, Ex.time for extraction time, Des.temp for desorption temperature, and Des.time for desorption time. Units are standard SI units for each factor.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lioupi, A.; Sampsonidis, I.; Virgiliou, C.; Papoti, V.T.; Zinoviadou, K.G.; Spyros, A.; Theodoridis, G. Optimisation of the HS-SPME/GC-MS Approach by Design of Experiments Combined with Chemometrics for the Classification of Cretan Virgin Olive Oils. Metabolites 2022, 12, 114. https://doi.org/10.3390/metabo12020114

AMA Style

Lioupi A, Sampsonidis I, Virgiliou C, Papoti VT, Zinoviadou KG, Spyros A, Theodoridis G. Optimisation of the HS-SPME/GC-MS Approach by Design of Experiments Combined with Chemometrics for the Classification of Cretan Virgin Olive Oils. Metabolites. 2022; 12(2):114. https://doi.org/10.3390/metabo12020114

Chicago/Turabian Style

Lioupi, Artemis, Ioannis Sampsonidis, Christina Virgiliou, Vassiliki T. Papoti, Kyriaki G. Zinoviadou, Apostolos Spyros, and Georgios Theodoridis. 2022. "Optimisation of the HS-SPME/GC-MS Approach by Design of Experiments Combined with Chemometrics for the Classification of Cretan Virgin Olive Oils" Metabolites 12, no. 2: 114. https://doi.org/10.3390/metabo12020114

APA Style

Lioupi, A., Sampsonidis, I., Virgiliou, C., Papoti, V. T., Zinoviadou, K. G., Spyros, A., & Theodoridis, G. (2022). Optimisation of the HS-SPME/GC-MS Approach by Design of Experiments Combined with Chemometrics for the Classification of Cretan Virgin Olive Oils. Metabolites, 12(2), 114. https://doi.org/10.3390/metabo12020114

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop