Next Article in Journal
Alpha-Glucosidase Inhibitory Assay-Screened Isolation and Molecular Docking Model from Bauhinia pulla Active Compounds
Next Article in Special Issue
Chemical Variability and In Vitro Anti-Inflammatory Activity of Leaf Essential Oil from Ivorian Isolona dewevrei (De Wild. & T. Durand) Engl. & Diels
Previous Article in Journal
Design, Synthesis, and Antibacterial Evaluation of Novel Ocotillol Derivatives and Their Synergistic Effects with Conventional Antibiotics
Previous Article in Special Issue
Chemical Composition, Antioxidant and Anti-Inflammatory Activities of Clary Sage and Coriander Essential Oils Produced on Polluted and Amended Soils-Phytomanagement Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Volatile Organic Compounds of Sideritis romana L. and Sideritis montana L. from Croatia

1
Department of Pharmaceutical Botany, University of Zagreb Faculty of Pharmacy and Biochemistry, Schrottova 39, 10 000 Zagreb, Croatia
2
Department of Food Technology and Biotechnology, Faculty of Chemistry and Technology, University of Split, Ruđera Boškovića 35, 21 000 Split, Croatia
3
Department of Organic Chemistry, Faculty of Chemistry and Technology, University of Split, Ruđera Boškovića 35, 21 000 Split, Croatia
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(19), 5968; https://doi.org/10.3390/molecules26195968
Submission received: 31 August 2021 / Revised: 22 September 2021 / Accepted: 28 September 2021 / Published: 1 October 2021

Abstract

:
A study on the headspace volatile organic compounds (VOCs) profile of native populations of Sideritis romana L. and Sidertis montana L., Lamiaceae, from Croatia is reported herein, to elucidate the phytochemical composition of taxa from this plant genus, well-known for traditional use in countries of the Mediterranean and the Balkan region. Headspace solid-phase microextraction (HS-SPME), using divinylbenzene/carboxene/polydimethylsiloxane (DVB/CAR/PDMS) or polydimethylsiloxane/divinylbenzene (PDMS/DVB) fiber, coupled with gas chromatography-mass spectrometry (GC-MS) was applied to analyze the dried aerial parts of six native populations in total. Furthermore, principal component analysis (PCA) was conducted on the volatile constituents with an average relative percentage ≥1.0% in at least one of the samples. Clear separation between the two species was obtained using both fiber types. The VOCs profile for all investigated populations was characterized by sesquiterpene hydrocarbons, followed by monoterpene hydrocarbons, except for one population of S. romana, in which monoterpene hydrocarbons predominated. To our knowledge, this is the first report on the VOCs composition of natural populations of S. romana and S. montana from Croatia as well as the first reported HS-SPME/GC-MS analysis of S. romana and S. montana worldwide.

Graphical Abstract

1. Introduction

Essential oils (EOs) are natural, volatile and complex mixtures of lipophilic compounds, often including terpenes, phenol-derived aromatic compounds and aliphatic compounds, that are biosynthesized by aromatic plants as their secondary metabolites and represented by a strong odor [1]. Ecologically, they are important for the plant defense system by providing protection against grazers, taking part in fire and drought tolerance and attracting pollinators and other animals for seed dispersal. The most often reported bioactivities of EOs and their components are antimicrobial, antiviral, antinociceptive, anticancer, anti-inflammatory, digestive, semiochemical and free radical scavenging [2], while new potential applications in human health, nutrition, agriculture and the environment are continuously emerging [1].
The volatile components of EOs are biosynthesized in specialized secretory structures such as glandular trichomes located on the surface, or secretory cells and ducts situated within the tissues of various plant organs of aromatic plants. The conventional techniques for volatile organic compounds (VOCs) extraction include, e.g., hydrodistillation and steam distillation, expression (cold pressing), solvent extraction [3] and enfleurage (a classical method for extracting volatile oil from flowers using a layer of animal fat) [4]. Along with conventional methods used for the extraction of volatile compounds from plants, solid-phase microextraction (SPME) is considered a convenient alternative technique for a variety of sample matrices (liquid/solid samples or headspace (HS) vapors). This technique bypasses the distillation or solvent extraction procedure [2]. It uses a fused silica fiber that is externally coated with an appropriate stationary phase such as polydimethylsiloxane (PDMS), polyacrylate (PA), PDMS/divinylbenzene (DVB), carbowax (CW)/DVB, and carboxene (CAR)/PDMS on which the headspace VOCs are adsorbed [5]. After the extraction, the fiber is usually directly inserted into the injector for gas chromatography analysis usually coupled with mass spectrometry (GC-MS), which is the most often used technique for the chemical characterization of EOs [2]. SPME has been widely used for chemical analysis of volatile components of different types of food substances, flavors, and medicinal plant materials because it is a simple, rapid, inexpensive and solvent-free technique that demands only a small amount of sample and reduces the degradation of VOCs [6,7]. Furthermore, headspace solid-phase microextraction (HS-SPME) requires no sample pre-treatment, allows high selectivity for small compounds and reusability of various commercially available SPME fibers [8], which makes this technique generally well accepted for VOCs extraction.
The genus Sideritis L. (family Lamiaceae) includes more than 150 species of annual or perennial herbs and shrubs, distributed in temperate and tropical regions of the Northern Hemisphere, with the majority of species found in the Mediterranean area, North Africa, the Iberian Peninsula, the Middle East and in the Macaronesian region [9,10,11]. The use of aerial parts of many Sideritis species has been well reported, especially in Albania, Bulgaria, Greece, Macedonia and Turkey, where they are consumed as herbal teas and applied in traditional medicine for the treatment of various disorders of the gastrointestinal and respiratory system, as well as for burns and wounds [10,12].
In recent decades, many researchers have been trying to justify the traditional uses of Sideritis species by elucidating their phytochemical composition and investigating potentially useful pharmacological activities of these species [10,13,14]. Although species from Lamiaceae family are well known for their aromatic properties and are recognized as the most important source of EOs with economical interest among Angiosperms [3], the species of the genus Sideritis have been reported as poor in EOs [10,15,16]. Nevertheless, for the pharmacological activity of Sideritis species, terpenes and flavonoids seem to be the constituents that are most important. Other chemical components which have been identified within the genus are iridoids, coumarins, lignans, phenylpropanoid glycosides and sterols [15].
Only three taxa of this genus are reported to be native to Croatia: Sideritis romana L., S. romana L. subsp. purpurea (Talbot ex Benth.) Heywood and S. montana L. [17]. While several reports on EOs of S. romana L. [18,19], S. romana subsp. purpurea [20,21] and S. montana L. [18,22,23,24,25,26] isolated by hydrodistillation have already been published, as well as that of S. romana subsp. purpurea EO constituents obtained by steam distillation [27], the composition of the headspace volatile components extracted using HS-SPME has never been evaluated for these species. Furthermore, to the best of our knowledge, the only report on using the HS-SPME technique for the detection of volatile components within the genus Sideritis has been published for S. ozturkii Aytac and Aksoy [28]. Additionally, the volatile HS profiles of S. scardica Griseb. and S. raeseri Boiss. and Heldr. were analyzed using GC-FID/MS [29], while EO composition of S. albiflora Hub.-Mor. was analyzed by both HS GC-MS and thermal desorption GC-MS [30].
The aim of the present study was to compare the VOCs composition of the aerial parts of S. romana L. and S. montana L. populations from different parts of Croatia and to determine characteristic volatile compounds for each species. Having in mind that the composition of volatile compounds may vary depending on environmental conditions [31,32], in comparison to similar studies conducted on hydrodistilled EOs of these species, which were based on the analysis of single populations [18,19,22,23,24,25,26], the advantage of our study is a greater number of investigated S. romana (four) and S. montana (two) populations. Considering that, on one hand, Sideritis species are generally poor in essential oil, and the limited availability of plant material in natural habitats from which the samples had been collected, and, on the other hand, the multiple advantages of HS-SPME, which have been mentioned above, instead of the conventionally used hydrodistillation, the latter technique was chosen for the isolation of volatile components, which were further analyzed by using GC-MS, with the aim to obtain more reliable results. Previous studies have shown that hydrodistillation and HS-SPME can detect similar major volatile constituents of individual plant species, although their relative percentages may vary depending on the technique used [33]. Having in mind that HS-SPME may, among other things, depend on the type of fiber used for extraction [34], the volatile components were extracted by HS-SPME using two different fibers suitable for untargeted analysis [35], divinylbenzene/carboxene/polydimethylsiloxane (DVB/CAR/PDMS) and polydimethylsiloxane/divinylbenzene (PDMS/DVB) fiber. To account for possible genetically and/or ecologically dependent variabilities in VOCs composition, several populations of S. romana (four populations) and S. montana (two populations) have been included in the analysis and principal component analysis (PCA) was conducted on the major volatile constituents (average relative percentage ≥ 1.0% in at least one of the samples) obtained by each fiber to examine the interrelationships among the investigated populations. Furthermore, the obtained results have been compared to previous findings for the same two species that were mostly based on hydrodistillation. Finally, the VOCs profile was compared to related species and discussed considering the reported biological activities of major constituents. To our knowledge, this study is the first report on the VOCs profile of natural populations of S. romana and S. montana from Croatia.

2. Results

2.1. HS-SPME/GC-MS Analysis

2.1.1. DVB/CAR/PDMS Fiber

VOCs of four populations of Sideritis romana and two populations of S. montana were extracted by HS-SPME, using DVB/CAR/PDMS fiber and their composition was analyzed by GC-MS. In total, 69 VOCs were identified, accounting for 86.18–90.29% of the total VOCs content (Table 1). The VOCs compositions of investigated populations of both investigated species were characterized by sesquiterpene hydrocarbons, followed by monoterpene hydrocarbons, except for that of S. romana from Morinje Bay (S. romana M), in which monoterpene hydrocarbons and other compounds were more abundant than sesquiterpene hydrocarbons. Additionally, for the latter, oxygenated sesquiterpenes were more represented than in other investigated samples (Figure 1). The main individual compounds that were identified for S. romana were bicyclogermacrene (15.16–23.04%) (except S. romana M), trans-caryophyllene (3.60–23.88%) (except S. romana from Kamenjak (S. romana K) and Blato (S. romana B)), trans-β-farnesene (1.00–12.79%), limonene (1.68–7.58%), alloaromadendrene (2.35–7.32%), β-pinene (2.61–8.01%), and isocaryophyllene (2.39–4.70%) (except S. romana M). Exceptionally, for S. romana M, the main identified compounds were benzyl alcohol (13.28%), viridiflorol (10.29%), ledene (8.20%), limonene (7.10%), β-pinene (6.67%), and p-cymene (5.49%). For S. montana, the main identified individual compounds were germacrene D (17.04–23.23%), trans-caryophyllene (6.56–11.89%), δ-cadinene (7.57–8.85%), bicyclogermacrene (4.19–11.76%), limonene (6.31–8.14%), and trans-β-farnesene (3.56–4.83%) (Table 1).

2.1.2. PDMS/DVB Fiber

VOCs of the same four populations of Sideritis romana and two populations of S. montana were extracted by HS-SPME, using also PDMS/DVB fiber and their composition was analyzed by GC-MS. In total, 56 VOCs were identified, accounting for 89.35–97.09% of the total VOCs content (Table 2). As with DVB/CAR/PDMS fiber, the VOCs compositions of investigated populations of both species were characterized by sesquiterpene hydrocarbons, followed by monoterpene hydrocarbons, except for that of S. romana M, in which monoterpene hydrocarbons and oxygenated sesquiterpenes were more abundant than sesquiterpene hydrocarbons. Moreover, other compounds were more abundant in this population than in other investigated samples (Figure 2). By applying PDMS/DVB fiber, more sesquiterpene hydrocarbons, and less oxygenated monoterpenes and other compounds were extracted than with the application of DVB/CAR/PDMS fiber. The main individual compounds that were identified for S. romana were bicyclogermacrene (10.21–57.50%), β-pinene (6.78–13.26%), isocaryophyllene (5.17–7.10%) (except S. romana M), trans-β-farnesene (3.55–7.71%), germacrene D (1.22–12.18%), γ-terpinene (0.11-8.96%), limonene (1.97–3.86%) and α-pinene (2.09–3.91%). For S. romana M, the main individual compounds were viridiflorol (24.19%), β-pinene (13.26%), bicyclogermacrene (10.21%), γ-terpinene (8.96%) and trans-caryophyllene (6.35%). In S. montana, the main individual compounds were germacrene D (51.08–53.63%), bicyclogermacrene (9.54–18.03%), trans-caryophyllene (6.65–15.49%), limonene (2.91–3.74%) and trans-β-farnesene (2.17–4.02%) (Table 2).

2.2. Principal Component Analysis

In order to analyze the differences between VOCs compositions of investigated S. romana and S. montana populations, principal component analysis of the major identified components of both taxa was performed, including only those compounds showing at least 1.0% of the total VOCs content in at least one of the samples (38 compounds in total), i.e., 37 compounds identified using DVB/CAR/PDMS fiber and 17 compounds identified using PDMS/DVB fiber (Table 3).

2.2.1. DVB/CAR/PDMS Fiber

The biplot constructed by the first two principal components that are showing the distribution of investigated S. romana and S. montana populations and VOCs identified with DVB/CAR/PDMS fiber is presented in Figure 3. Principal component 1 (PC1) accounted for 49.37% and principal component 2 (PC2) for 30.62% of total variance in the data. Clear separation between S. romana and S. montana was obtained. For all populations of S. romana, except for S. romana M, distinctive components for discrimination observed by using DVB/CAR/PDMS fiber were bicyclogermacrene (C29), alloaromadendrene (C21), trans-β-farnesene (C23), isocaryophyllene (C24), selina-5,11-diene (C22), α-gurjunene (C17) and spathulenol (C36), while S. romana M, was characterized by benzyl alcohol (C11), viridiflorol (C37), ledene (C28), p-cymene (C9), acetic acid (C1), benzaldehyde (C4) and β-myrcene (C6). Distinctive components for discrimination of S. montana populations were germacrene D (C26), δ-cadinene (C32), γ-cadinene (C31), α-amorphene (C25), α-muurolene (C30), α-copaene (C16), α-cadinene (C34), cadina-1,4-diene (C33) and α-calacorene (C35).

2.2.2. PDMS/DVB Fiber

The biplot constructed by the first two principal components that are showing the distribution of investigated S. romana and S. montana populations and VOCs identified with PDMS/DVB fiber is presented in Figure 4. Principal component 1 (PC1) accounted for 55.88% and principal component 2 (PC2) for 33.21% of total variance in the data. Clear separation between S. romana and S. montana was obtained. Distinctive components observed by using PDMS/DVB fiber for all S. romana populations, except for S. romana M, were bicyclogermacrene (C15′), isocaryophyllene (C13′), trans-β-farnesene (C12′) and bicycloelemene (C9′). S. romana M population was characterized by viridiflorol (C17′), γ-terpinene (C7′), α-pinene (C1′), p-cymene (C4′), benzyl alcohol (C6′) and α-terpinene (C3′). Distinctive components for discrimination of S. montana were germacrene D (C14′), α-copaene (C10′) and δ-cadinene (C16′).

3. Discussion

Depending on the polarity of VOCs of an analyte, several types of fibers may be used for extracting different groups of compounds such as PDMS (a non-polar fiber used for volatile compounds), PA (a polar fiber used for polar semi-volatile compounds), CW/DVB (a polar fiber used for alcohols and volatiles), PDMS/DVB (a bipolar fiber used for volatile compounds of medium polarity, amines and nitroaromatics), CAR/PDMS (a bipolar fiber used for low molecular weight volatile compounds) and DVB/CAR/PDMS (a bipolar fiber used for polar and non-polar volatile compounds) [34,36]. In this study, DVB/CAR/PDMS and PDMS/DVB fiber, which were previously shown to be most suitable for untargeted HS-SPME analysis of volatiles [35], were used for the analysis of VOCs from aerial parts of S. romana and S. montana from Croatia. In our study, the number of extracted VOCs was higher using DVB/CAR/PDMS fiber (69 compounds) than it was by using PDMS/DVB fiber (56 compounds), which is in accordance with the work of Wang et al. [37], in which DVB/CAR/PDMS fiber showed better efficiency to extract volatile compounds (50 compounds) from the samples of Aquilegia japonica Nakai and H.Hara, compared to CAR/PDMS (47 compounds) and PDMS/DVB fibers (45 compounds), as well as with the work of Sukkaew et al. [38], where more volatile components were extracted using DVB/CAR/PDMS fiber (51 compounds) than with PDMS/DVB (38 compounds), PDMS (38 compounds) and CAR/PDMS (37 compounds) fibers in Murraya koenigii (L.) Sprengel fresh leaves. In our study, the dominant chemical classes of volatile compounds in all samples of the two investigated Sideritis species extracted using DVB/CAR/PDMS and PDMS/DVB fibers were sesquiterpene hydrocarbons followed by monoterpene hydrocarbons, except for S. romana M in whose extract monoterpene hydrocarbons were more abundantly present. The obtained results for S. montana are in accordance with the results reported by Venditti et al. [23] for EO of S. montana subsp. montana from central Italy obtained by hydrodistillation, where sesquiterpene hydrocarbons led by germacrene D and bicyclogermacrene were observed as the most abundant EO components. According to the latter, oil-poor species from Lamioideae subfamily produce EOs rich in sesquiterpene hydrocarbons. Furthermore, our results are in accordance with the classification of Sideritis species from Turkey given by Kirimer et al. [18], based on the main components of hydrodistilled EOs, where S. montana ssp. montana was classified into the sesquiterpene hydrocarbons-rich group, and S. romana ssp. romana into the oxygenated monoterpenes-rich group. Considering the relative percentages of extracted chemical classes, in comparison to PDMS/DVB fiber, DVB/CAR/PDMS fiber extracted more oxygenated monoterpenes and other compounds (mostly non-terpenes), which were more abundantly present in investigated S. romana populations. Conversely, in accordance with previously reported findings [36], PDMS/DVB fiber was observed to be more efficient in extracting sesquiterpene hydrocarbons, which were relatively higher in investigated S. montana populations when comparing to the results obtained for the S. romana populations included in this study. In fact, for both species, the most abundant VOCs, namely the sesquiterpene hydrocarbons germacrene D (S. montana) and bicyclogermacrene (S. romana, except for S. romana M) were extracted more efficiently by PDMS/DVB fiber. The same was observed for viridiflorol, an oxygenated sesquiterpene that was found as the major characteristic compound of the chemically distinct population of S. romana from Morinje Bay (S. romana M).
The dominant compounds of S. romana (except S. romana M) extracted with DVB/CAR/PDMS fiber were bicyclogermacrene, trans-caryophyllene, trans-β-farnesene, limonene and alloaromadendrene, while bicyclogermacrene, β-pinene, isocaryophyllene, trans-β-farnesene, germacrene D and γ-terpinene were dominant after the extraction with PDMS/DVB fiber. A different composition of EO was reported for S. romana from Italy, where carvacrol, limonene and 1,8-cineole were reported as the dominant compounds isolated by hydrodistillation [19], while neither carvacrol and 1,8-cineole nor thymol, oct-1-en-3-ol and borneol, reported as the dominant compounds isolated by hydrodistillation from S. romana subsp. romana from Turkey [18], were detected in our samples of S. romana from Croatia. For S. romana M, the dominant compounds extracted with DVB/CAR/PDMS were benzyl alcohol, viridiflorol, ledene and limonene, while with PDMS/DVB these were viridiflorol, β-pinene, bicyclogermacrene and γ-terpinene. These results are not in accordance with the previously mentioned studies of S. romana from Italy and Turkey. For other subspecies of S. romana, such as S. romana subsp. purpurea from Greece, bicyclogermacrene, β-caryophyllene (=trans-caryophyllene or isocaryophyllene), γ-muurolene, β-pinene, (E)-β-farnesene (=trans-β-farnesene) and spathulenol were found as the major components obtained by hydrodistillation [21], while for the same subspecies from Montenegro, the major EOs constituents obtained by steam distillation were γ-elemene and spathulenol [27]. Moreover, the most abundant EO components obtained by hydrodistillation from the same subspecies, which was also collected in Montenegro, included bicyclogermacrene, germacrene D, (E)-caryophyllene and spathulenol [20]. According to the PCA analysis, even though it was present with more than 1% only in one of the analyzed samples, spathulenol appeared to be one of the compounds distinctive for S. romana. However, it was not found in S. romana M, which was observed to contain much lower amounts of bicyclogermacrene, which were detected only using PDMS/DVB fiber. However, spathulenol may be an artifact, having in mind that bicyclogermacrene is easily converted to spathulenol at room temperature [39]. The compound may also be formed during the process of hydrodistillation, as seen in the comparison of EO composition and VOCs composition of S. scardica samples prepared by other techniques of isolation such as solvent extraction and supercritical carbon dioxide extraction [40]. On the other hand, distinctive VOCs of S. romana M such as viridiflorol (detected in higher amount using PDMS/DVB fiber) and ledene (detected only using DVB/CAR/PDMS fiber) may be the result of bicyclogermacrene hydration and/or rearrangement [41,42]. According to Carvalho et al. [43], who studied the volatile fractions of in natura, fresh, and dried Casearia sylvestris var. sylvestris Sw. and var. lingua (Cambess.) Eichler leaves, viridiflorol is an artifact most likely formed from bicyclogermacrene, whose content increases during the drying process. As it can be seen from Table 4, the distinct sample of S. romana from Morinje Bay (S. romana M) was collected prior to other samples of the same species and, therefore, may have been more susceptible to transformations during processing and storage of plant material, having in mind that all collected samples were analyzed simultaneously. The sample from Morinje Bay also contained more other compounds (non-terpenes), which were mostly oxygenated compounds (alcohols, aldehydes and a carboxylic acid), indicating its possible oxidative degradation. Other compounds that were present and/or were more abundant in S. romana populations other than S. romana M, such as alloaromadendrene and aromadendrene, may also have been formed by isomerization from bicyclogermacrene [44] and could have produced spathulenol by further oxidation [43].
Comparable to our results, variations in bicyclogermacrene content were observed in aroma compounds isolated by hydrodistillation of dried aerial parts from different populations of S. scardica and S. raeseri from Macedonia [29]. Moreover, it was reported by Kirimer et al. [46] in their extensive study conducted on 50 Sideritis taxa from the section Empedoclia that viridiflorol was one of the major components isolated by hydrodistillation from S. perfoliata L. from Turkey, while the same compound was not observed for this species in a subsequent study done by the same authors conducted on only two species [47]. Additionally, viridiflorol was found as one of the major constituents of EO isolated by hydrodistillation from S. montana subsp. montana from Turkey together with germacrene D, bicyclogermacrene and (E)-β-farnesene [18]. In the present study, the predominant compounds for S. montana extracted with DVB/CAR/PDMS fiber were germacrene D, trans-caryophyllene, δ-cadinene, bicyclogermacrene, limonene and trans-β-farnesene, while with PDMS/DVB fiber germacrene D was followed by bicyclogermacrene, trans-caryophyllene, limonene and trans-β-farnesene. In investigated samples of S. montana, as well as in the samples of S. romana other than S. romana M, viridiflorol was found only in trace amounts (<1%), while ledene was observed only in S. montana from Ježević (S. montana J), the latter being the previously collected and thus older of the two investigated samples of S. montana (stored for a slightly longer period before GC-MS analysis). Considering that a significantly lower amount of bicyclogermacrene was detected in the same sample compared to S. montana from Mosor (S. montana M), ledene may have been formed during the drying and storage period from the aforementioned bicyclogermacrene. Our results considering the major VOCs of investigated populations of S. montana are in accordance with the results of previous studies from Bulgaria [22], Turkey [18], Serbia [25] and Italy [23], in which germacrene D predominated in the chemical composition of EOs isolated from aerial parts by hydrodistillation, usually being followed by bicyclogermacrene. Different major compounds were reported for S. montana from Iran, in which geraniol particularly predominated in the EO obtained by hydrodistillation from the flowering spikes of the plant [26], while it was not detected in our samples of S. montana aerial parts, as well as for S. montana subsp. montana from Turkey, in which β-caryophyllene, α-pinene and β-pinene were found to be the major compounds of the hydrodistilled EO obtained from aerial parts of the plant [24]. The present study also indicated that S. montana, unlike S. romana, is characterized by several sesquiterpene hydrocarbons having a cadalane skeleton (δ-cadinene, γ-cadinene, α-amorphene, α-muurolene, α-cadinene, and cadina-1,4-diene). Clear separation between S. montana and S. romana samples was observed in the PCA with both fiber types used, with S. romana M showing a distinct VOCs profile.
Observed differences in the chemical composition of investigated Croatian Sideritis species from previously published data on VOCs composition of aerial parts from samples of different geographical origins may be explained by different ecological conditions (Table 4). According to Köppen climate classification [45], the collection areas from which the populations included in the present study were sampled belong to different climates including the Mediterranean climate with hot summer (S. romana B, S. romana M), the Mediterranean climate with warm summer (S. romana from Pražnica (S. romana P)), the temperate humid climate with hot summer (S. romana K, S. montana M) and the temperate humid climate with warm summer (S. montana J). This might have party affected their VOCs production. As it is known, VOCs composition can change depending on the type of stimulus received from the environment [32]. For example, variation of VOCs content may have been caused by environmental stress such as drought stress, as it was detected in the study of different cultivars of Thymus vulgaris L., in which α-phellandrene, o-cymene, γ-terpinene and β-caryophyllene were recognized as the compounds included in drought stress adaptation [48]. These compounds were observed as one of the major identified VOCs in the present study. For example, γ-terpinene was found to be especially dominant in S. romana M. Additionally, individual VOCs production might have been affected by the plant species that grew in the vicinity of investigated specimens in each habitat (place of collection) [49,50]. Although this has not been recorded in the present study, insight into the diversity of plant species in the harvesting locations can be partially gained through floristic studies that were conducted in the areas that are overlapping with or are near to the harvesting locations [51,52,53,54,55,56,57,58,59].
Besides ecological conditions, the observed differences may also be explained by the applied methods of extraction, having in mind that none of the previously mentioned studies used the HS-SPME method for extraction of VOCs from S. romana and S. montana. Instead, as already mentioned above, most authors used hydrodistillation, with the exception of the study done by Garzolli et al. [27], in which steam distillation was used. A few comparative analyses of plant volatile compounds isolated by hydrodistillation and HS-SMPE were published. A study on Myrtus communis L. showed that the applied microextraction techniques (HS-SPME and HS-SDME) extracted compounds that were more volatile, such as α-pinene and limonene, in comparison to hydrodistillation, which, in turn, obtained higher peak areas for low volatile compounds [60]. Moreover, analysis of volatile compounds from fruits of Seseli libanotis (L.) W.D.J.Koch also showed that higher contents of low-boiling compounds such as sabinene, β-phellandrene, α-pinene, β-pinene, β-myrcene, γ-terpinene and α-phellandrene may be extracted by using HS-SPME than by hydrodistillation, while higher amounts of high-boiling compounds (higher molecular mass and lower volatility) were observed using hydrodistillation at high temperature [61]. Qualitative and quantitative differences of volatile compounds of Melissa officinalis L. were observed after applying hydrodistillation and HS-SPME, which may be attributed to the formation of artifacts [6]. Furthermore, variations in the percentages and nature of compounds adsorbed on a SPME fiber compared to those found in hydrodistilled EOs were reported for Petroselinum crispum Mill. [62], while the comparison of the same techniques on Bupleurum plantagineum Desf. showed differences mainly in minor components [63]. Therefore, despite some differences, the results obtained by HS-SPME considering the major VOCs may be comparable to the results obtained by hydrodistillation [61,62,63].
In the studied samples, the most abundant VOCs, depending on the fiber used for solid-phase microextraction, were bicyclogermacrene (S. romana K, S. romana B and S. romana P) and/or trans-caryophyllene (S. romana P), benzyl alcohol or viridiflorol (S. romana M), and germacrene D (S. montana J and S. montana M). These compounds were previously shown to possess certain biological activities that may be of interest to the pharmaceutical, food and/or cosmetics industry. For example, bicyclogermacrene displayed cytotoxic activity against a range of cancer cell lines (IC50 = 3.1–21 μg/mL [64] and IC50 = 1.5–4.4 μg/mL). The same was observed for germacrene D (IC50 = 2.7–8.0 μg/mL) [65]. Additionally, germacrene D and trans-caryophyllene were recognized as the major EO constituents of selected species from the genera Siparuna Aublet and Piper L. that possess antiradical activity [66]. According to Dahham et al. [67], β-caryophyllene possesses antibacterial activity against S. aureus (MIC = 3 ± 1.0 µM) as well as anti-fungal (MIC = 4–7 µM), antioxidant (IC50 = 1.25–3.23 µM), and anti-proliferative activity against colorectal cancer cells (IC50 = 19 µM). Benzyl alcohol is frequently used as a bacteriostatic agent in liquid pharmaceutical products [68] as well as a preservative in cosmetic products [69]. On the other hand, antiradical (IC50 = 57.55–74.7 μg/mL), anti-mycobacterial (MIC = 190.0 μg/mL) and anti-inflammatory activity (60% reduction of carrageenan-induced mice paw edema after 3–30 mg/kg oral administration) were reported for viridiflorol [70]. Moreover, viridiflorol was recently observed to induce intracellular Ca2+ mobilization in human neutrophils and C20 microglial cells [71].
Related species of the genus Sideritis, namely S. scardica [22,29,72], S. clandestina (Bory and Chaub.) Hayek [21], S. raeseri [72,73,74], and S. syriaca L. [22], whose aerial parts have been traditionally used for the treatment of inflammation, gastrointestinal disorders and cough associated with cold and whose usage has been approved by the European Medicines Agency (EMA) Committee on Herbal Medicinal Products (HMPC) [75], were observed to possess, as one of their major volatile constituents, the same VOCs that were found in S. romana and S. montana in the present study, thus, indicating a promising potential for future investigations of these species.

4. Materials and Methods

4.1. Plant Material

Aerial parts of S. romana (four samples, consisting of 5–15 cm long shoots) and S. montana (two samples, consisting of 7–26 cm long shoots) were collected during the flowering period in June 2020 from natural populations on several locations in Croatia (Figure 5). In the field, the plant material was stored in paper bags and later spread out into one layer to dry out, in a place with room temperature. The mass of dried samples ranged between 2 and 22 g for S. romana and between 3 and 11 g for S. montana. Taxonomic identification was performed by T. M. Voucher specimens have been deposited within the Herbarium of the Department of Pharmaceutical Botany, University of Zagreb Faculty of Pharmacy and Biochemistry, Croatia. Sample codes, voucher numbers, harvesting locations and geographic coordinates are given in Table 4.

4.2. SPME Fibers and Extraction Procedure

HS-SPME was achieved with a manual SPME holder by using two fibers, (divinylbenzene/carboxene/polydimethylsiloxane (DVB/CAR/PDMS) and polydimethylsiloxane/divinylbenzene (PDMS/DVB)), that were conditioned according to Supelco Co. instructions before extraction. Cut samples (1 g) were placed separately in glass vials (5 mL) and sealed hermetically using PTFE/silicone septa. The vials were maintained in a water bath (60 °C) during equilibration (15 min) and extraction by HS-SPME (45 min). After extraction, the SPME fiber was withdrawn and inserted into the GC-MS injector (250 °C) for thermal desorption (6 min). The treatment was similar as previously reported [76]. HS-SPME was performed in duplicate and average values are presented in Table 1, Table 2 and Table 3.

4.3. GC-MS Analysis

GC-MS analyses were performed on a gas chromatograph model 7820A (Agilent Technologies, Palo Alto, CA, USA) containing a HP-5MS capillary column (5% phenyl-methylpolysiloxane, Agilent J and W; 30 m × 0.25 mm i.d., coating thickness 0.25 µm) and a mass selective detector (MSD) model 5977E (Agilent Technologies, Palo Alto, CA, USA). The GC conditions were described previously [76,77]. The carrier gas was helium (He 1.0 mL/min). The oven temperature was set at 70 °C for 2 min, then it was increased from 70 °C to 200 °C at a rate of 3 °C/min, and held at 200 °C for 15 min. The MSD (EI mode) was used at 70 eV, and 30–300 amu mass range was applied.
The compounds identification was based on the retention indices (RIs) determined relative to retention times of n-alkanes (C9-C25) and their comparison with literature data (National Institute of Standards and Technology) as well as by their mass spectra compared with the spectra from Wiley 9 (Wiley, New York, NY, USA) and NIST 17 (Gaithersburg, MD, USA) mass spectral libraries. The percentage composition was calculated using the normalization method (without correction factors). The average component percentages in Table 1–3 were calculated from duplicate GC-MS analyses [76,77].

4.4. Principal Component Analysis

Principal component analysis (PCA) was conducted on the volatile constituents having an average relative percentage ≥ 1.0% in at least one of the samples (in total, six observation of 37 variables for DVB/CAR/PDMS fiber and six observations of 17 variables for PDMS/DVB fiber), in order to examine the interrelationships among the investigated populations of S. romana (four populations) and S. montana (two populations), in RStudio version 1.4.1717 [78] using the function prcomp in R version 4.1.1. [79]. Scaling was set to “TRUE” to perform the analysis on normalized data. Plotting was performed using the function autoplot from the package ggfortify [80,81].

5. Conclusions

Headspace volatile organic compounds (VOCs) profile of aerial parts of altogether six native populations of Sideritis romana and S. montana from Croatia were determined by HS-SPME/GC-MS, using DVB/CAR/PDMS and PDMS/DVB fibers for compound extraction, which resulted in identification of 69 and 56 compounds, respectively. The most abundant VOCs in both species were those belonging to the class of sesquiterpene hydrocarbons, which were extracted more efficiently using PDMS/DVB fiber. The performed PCA analyses highlighted the major volatile constituents for each population and revealed clear separation between the investigated species using both fiber types. The most abundant VOCs found in the analyzed samples of S. romana were bicyclogermacrene and/or trans-caryophyllene or viridiflorol and benzyl alcohol, while germacrene D was the most abundant VOC found in the analyzed samples of S. montana, all of which were previously reported to possess certain biological activities that are potentially interesting for the pharmaceutical, food and/or cosmetics industry. The observed differences in the presence and amount of detected volatile compounds among the researched populations of the same species may be a result of different environmental and ecological conditions existing on the collection sites, while specific differences may also be attributed to changes that might have occurred during drying and storage of sampled plant material. The obtained results, according to which S. romana and S. montana show a promising potential for future utilization, suggest that the origin of plant material and/or growing as well as processing conditions should be considered as possible factors affecting VOCs production if these species were to be exploited.

Author Contributions

Conceptualization, T.M. and M.F.; collection of plant material, T.M.; HS-SPME and GC-MS analysis, Z.M. and I.J.; PCA analysis, M.F; data analysis, T.M. and M.F.; supervision, Ž.M. and I.J.; visualization, M.F.; writing—original draft, T.M. and M.F.; writing—review and editing, T.M., M.F., I.J. and Ž.M.; funding acquisition, M.F. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by the University of Zagreb (Z-209).

Data Availability Statement

The data presented in this study are available in article.

Acknowledgments

The authors would like to thank the University of Zagreb for financial support.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are available from the authors.

References

  1. Bakkali, F.; Averbeck, S.; Averbeck, D.; Idaomar, M. Biological effects of essential oils—A review. Food Chem. Toxicol. 2008, 46, 446–475. [Google Scholar] [CrossRef]
  2. Sadgrove, N.; Jones, G. A contemporary introduction to essential oils: Chemistry, bioactivity and prospects for Australian agriculture. Agriculture 2015, 5, 48–102. [Google Scholar] [CrossRef] [Green Version]
  3. de Mesquita, L.S.S.; Luz, T.R.S.A.; de Mesquita, J.W.C.; Coutinho, D.F.; do Amaral, F.M.M.; de Sousa Ribeiro, M.N.; Malik, S. Exploring the anticancer properties of essential oils from family Lamiaceae. Food Rev. Int. 2019, 35, 105–131. [Google Scholar] [CrossRef]
  4. Paibon, W.; Yimnoi, C.-A.; Tembab, N.; Boonlue, W.; Jampachaisri, K.; Nuengchamnong, N.; Waranuch, N.; Ingkaninan, K. Comparison and evaluation of volatile oils from three different extraction methods for some Thai fragrant flowers. Int. J. Cosmet. Sci. 2011, 33, 150–156. [Google Scholar] [CrossRef] [PubMed]
  5. Sagratini, G.; Maggi, F.; Bílek, T.; Papa, F.; Vittori, S. Analysis of the volatile compounds of Teucrium flavum L. subsp. flavum (Lamiaceae) by headspace solid-phase microextraction coupled to gas chromatography with flame ionisation and mass spectrometric detection. Nat. Prod. Res. 2012, 26, 1339–1347. [Google Scholar] [CrossRef] [PubMed]
  6. Rehman, S.-U.; Latief, R.; Bhat, K.A.; Khuroo, M.A.; Shawl, A.S.; Chandra, S. Comparative analysis of the aroma chemicals of Melissa officinalis using hydrodistillation and HS-SPME techniques. Arab. J. Chem. 2017, 10, S2485–S2490. [Google Scholar] [CrossRef] [Green Version]
  7. El-Sakhawy, F.S.; Kassem, H.A.; El-Gayed, S.H.; Mostafa, M.M. Headspace solid phase microextraction analysis of volatile compounds of the aerial parts and flowers of Plectranthus neochilus Schltr. and Salvia farinacea Benth. J. Essent. Oil Bear. Plants 2018, 21, 674–686. [Google Scholar] [CrossRef]
  8. Jerković, I.; Marijanović, Z.; Radonić, A.; Zekić, M.; Kranjac, M. The application of headspace solid-phase microextraction as a preparation approach for gas chromatography with mass spectrometry. Kem. Ind. 2020, 69, 515–520. [Google Scholar] [CrossRef]
  9. Fraga, B.M. Phytochemistry and chemotaxonomy of Sideritis species from the Mediterranean region. Phytochemistry 2012, 76, 7–24. [Google Scholar] [CrossRef]
  10. González-Burgos, E.; Carretero, M.E.; Gómez-Serranillos, M.P. Sideritis spp.: Uses, chemical composition and pharmacological activities—A review. J. Ethnopharmacol. 2011, 135, 209–225. [Google Scholar] [CrossRef]
  11. Sideritis L. Available online: http://www.worldfloraonline.org/taxon/wfo-4000035314 (accessed on 15 July 2021).
  12. Aneva, I.; Zhelev, P.; Kozuharova, E.; Danova, K.; Nabavi, S.F.; Behzad, S. Genus Sideritis, section Empedoclia in southeastern Europe and Turkey—studies in ethnopharmacology and recent progress of biological activities. DARU J. Pharm. Sci. 2019, 27, 407–421. [Google Scholar] [CrossRef]
  13. Cavalcanti, M.R.M.; Passos, F.R.S.; Monteiro, B.S.; Gandhi, S.R.; Heimfarth, L.; Lima, B.S.; Nascimento, Y.M.; Duarte, M.C.; Araujo, A.A.S.; Menezes, I.R.A.; et al. HPLC-DAD-UV analysis, anti-inflammatory and anti-neuropathic effects of methanolic extract of Sideritis bilgeriana (lamiaceae) by NF-κB, TNF-α, IL-1β and IL-6 involvement. J. Ethnopharmacol. 2021, 265, 113338. [Google Scholar] [CrossRef] [PubMed]
  14. Deveci, E.; Tel-Çayan, G.; Duru, M.E.; Öztürk, M. Phytochemical contents, antioxidant effects, and inhibitory activities of key enzymes associated with Alzheimer’s disease, ulcer, and skin disorders of Sideritis albiflora and Sideritis leptoclada. J. Food Biochem. 2019, 43, e13078. [Google Scholar] [CrossRef] [PubMed]
  15. Bojović, D.; Janković, S.; Potpara, Z.; Tadić, V. Summary of the phytochemical research performed to date on Sideritis species. Ser. J. Exp. Clin. Res. 2011, 12, 109–122. [Google Scholar] [CrossRef] [Green Version]
  16. Solomou, A.D.; Skoufogianni, E.; Mylonas, C.; Germani, R.; Danalatos, N.G. Cultivation and utilization of “Greek mountain tea” (Sideritis spp.): Current knowledge and future challenges. Asian J. Agric. Biol. 2019, 7, 289–299. [Google Scholar]
  17. Domac, R. Flora Hrvatske, 2nd ed.; Školska Knjga: Zagreb, Croatia, 2002; pp. 291–292. [Google Scholar]
  18. Kirimer, N.; Tabanca, N.; Özek, T.; Tümen, G.; Başer, K.H.C. Essential oils of annual Sideritis species growing in Turkey. Pharm. Biol. 2000, 38, 106–111. [Google Scholar] [CrossRef]
  19. Flamini, G.; Cioni, P.L.; Morelli, I. Characterization of the volatile fraction of a Sideritis romana population of Montemarcello (Eastern Liguria). J. Essent. Oil Res. 1994, 6, 239–242. [Google Scholar] [CrossRef]
  20. Tadić, V.; Oliva, A.; Božović, M.; Cipolla, A.; De Angelis, M.; Vullo, V.; Garzoli, S.; Ragno, R. Chemical and antimicrobial analyses of Sideritis romana L. subsp. purpurea (Tal. ex Benth.) Heywood, an endemic of the Western Balkan. Molecules 2017, 22, 1395. [Google Scholar] [CrossRef] [Green Version]
  21. Koutsaviti, A.; Bazos, I.; Milenković, M.; Pavlović-Drobac, M.; Tzakou, O. Antimicrobial activity and essential oil composition of five Sideritis taxa of Empedoclia and Hesiodia sect. from Greece. Rec. Nat. Prod. 2013, 7, 6–14. [Google Scholar]
  22. Todorova, M.N.; Christov, R.C. Essential oil composition of three Sideritis species from Bulgaria. J. Essent. Oil Res. 2000, 12, 418–420. [Google Scholar] [CrossRef]
  23. Venditti, A.; Bianco, A.; Frezza, C.; Serafini, M.; Giacomello, G.; Giuliani, C.; Bramucci, M.; Quassinti, L.; Lupidi, G.; Lucarini, D.; et al. Secondary metabolites, glandular trichomes and biological activity of Sideritis montana L. subsp. montana from Central Italy. Chem. Biodivers. 2016, 13, 1380–1390. [Google Scholar] [CrossRef] [PubMed]
  24. Kilic, O. Essential oil composition of two Sideritis L. taxa from Turkey: A chemotaxonomic approach. Asian J. Chem. 2014, 26, 2466–2470. [Google Scholar] [CrossRef]
  25. Miladinović, D.L.; Ilić, B.S.; Mihajilov-Krstev, T.M.; Nikolić, N.D.; Milosavljević, V.N.; Nikolić, D.M. Antibakterijski potencijal etarskog ulja Sideritis montana L. (Lamiaceae). Hem. Ind. 2012, 66, 541–545. [Google Scholar] [CrossRef] [Green Version]
  26. Meshkatalsadat, M.H.; Sarabi, R.S.; Akbari, N.; Pireai, M. Composition of essential oils from Sideritis montana of Iran. Asian J. Chem. 2007, 19, 5769–5771. [Google Scholar]
  27. Garzoli, S.; Božović, M.; Baldisserotto, A.; Andreotti, E.; Pepi, F.; Tadić, V.; Manfredini, S.; Ragno, R. Sideritis romana L. subsp. purpurea (Tal. ex Benth.) Heywood, a new chemotype from Montenegro. Nat. Prod. Res. 2017, 32, 1056–1061. [Google Scholar] [CrossRef]
  28. Uysal, T.; Sezer, E.N.Ş.; Bozkurt, M. Headspace solid phase microextraction (HS-SPME) analysis of Sideritis ozturkii Aytac & Aksoy. In Proceedings of the 5th International Mediterranean Symposium on Medicinal and Aromatic Plants Book, Cappadocia, Turkey, 24–26 April 2019. [Google Scholar]
  29. Qazimi, B.; Stefkov, G.; Karapandzova, M.; Cvetkovikj, I.; Kulevanova, S. Aroma compounds of mountain tea (Sideritis scardica and S. raeseri) from Western Balkan. Nat. Prod. Commun. 2014, 9, 1369–1372. [Google Scholar] [CrossRef] [Green Version]
  30. Topçu, G.; Barla, A.; Gören, A.C.; Bilsel, G.; Bilsel, M.; Tümen, G. Analysis of the essential oil composition of Sideritis albiflora using direct thermal desorption and headspace GC-MS techniques. Turk. J. Chem. 2005, 29, 525–529. [Google Scholar]
  31. Aboukhalid, K.; Al Faiz, C.; Douaik, A.; Bakha, M.; Kursa, K.; Agacka-Mołdoch, M.; Machon, N.; Tomi, F.; Lamiri, A. Influence of environmental factors on essential oil variability in Origanum compactum Benth. growing wild in Morocco. Chem. Biodivers. 2017, 14, e1700158. [Google Scholar] [CrossRef]
  32. Kumar, A.; Memo, M.; Mastinu, A. Plant behaviour: An evolutionary response to the environment? Plant Biol. 2020, 22, 961–970. [Google Scholar] [CrossRef]
  33. Wajs, A.; Pranovich, A.; Reunanen, M.; Willför, S.; Holmbom, B. Characterisation of volatile organic compounds in stemwood using solid-phase microextraction. Phytochem. Anal. 2006, 17, 91–101. [Google Scholar] [CrossRef] [PubMed]
  34. Dadalı, C.; Elmacı, Y. Optimization of headspace-solid phase microextraction (HS-SPME) technique for the analysis of volatile compounds of margarine. J. Food Sci. Technol. 2019, 56, 4834–4843. [Google Scholar] [CrossRef]
  35. Pinheiro, G.P.; Galbiatti, M.I.; Carneiro, M.J.; Sawaya, A.C.H.F. Comparison of four different solid-phase microextraction fibers for analysis of Plectranthus amboinicus (Lour.) Spreng. leaf volatiles. Adv. Med. Plant. Res. 2019, 7, 38–43. [Google Scholar] [CrossRef]
  36. Jerković, I.; Marijanović, Z. A short review of headspace extraction and ultrasonic solvent extraction for honey volatiles fingerprinting. Croat. J. Food. Sci. Technol. 2009, 1, 28–34. [Google Scholar]
  37. Wang, H.-Y.; Zhang, W.; Dong, J.-H.; Wu, H.; Wang, Y.-H.; Xiao, H.-X. Optimization of SPME–GC–MS and characterization of floral scents from Aquilegia japonica and A. amurensis flowers. BMC Chem. 2021, 15, 26. [Google Scholar] [CrossRef] [PubMed]
  38. Sukkaew, S.; Pripdeevech, P.; Thongpoon, C.; Machan, T.; Wongchuphan, R. Volatile constituents of Murraya koenigii fresh leaves using headspace solid phase microextraction—gas chromatography—mass spectrometry. Nat. Prod. Commun. 2014, 9, 1783–1786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Toyota, M.; Koyama, H.; Mizutani, M.; Asakawa, Y. (−)-ent-spathulenol isolated from liverworths is an artefact. Phytochemistry 1996, 41, 1347–1350. [Google Scholar] [CrossRef]
  40. Tadić, V.; Bojović, D.; Arsić, I.; Đorđević, S.; Aksentijevic, K.; Stamenić, M.; Janković, S. Chemical and antimicrobial evaluation of supercritical and conventional Sideritis scardica Griseb., Lamiaceae extracts. Molecules 2012, 17, 2683–2703. [Google Scholar] [CrossRef] [Green Version]
  41. Jakab, E.; Sebestyén, Z.; Babinszki, B.; Barta-Rajnai, E.; Czégény, Z.; Nicol, J.; Clayton, P.; Liu, C. Thermo-oxidative decomposition of lovage (Levisticum officinale) and davana (Artemisia pallens) essential oils under simulated tobacco heating product conditions. Beitr. Tabakforsch. Int. 2020, 29, 27–43. [Google Scholar] [CrossRef]
  42. Tran, D.N.; Cramer, N. Biomimetic synthesis of (+)-ledene, (+)-viridiflorol, (−)-palustrol, (+)-spathulenol, and psiguadial A, C, and D via the platform terpene (+)-bicyclogermacrene. Chem. Eur. J. 2014, 20, 10654–10660. [Google Scholar] [CrossRef]
  43. Carvalho, F.A.; Oda, F.B.; Crevelin, E.J.; Crotti, A.E.M.; dos Santos, A.G. Casearia sylvestris essential oil degradation products generated by leaf processing. Chem. Biodivers. 2021, 18, e2000880. [Google Scholar] [CrossRef]
  44. Njoroge, S.M.; Ukeda, H.; Sawamura, M. Changes in the volatile composition of yuzu (Citrus junos Tanaka) cold-pressed oil during storage. J. Agric. Food Chem. 1996, 44, 550–556. [Google Scholar] [CrossRef]
  45. Šegota, T.; Filipčić, A. Köppenova podjela klima i hrvatsko nazivlje. Geoadria 2003, 8, 17–37. [Google Scholar] [CrossRef] [Green Version]
  46. Kirimer, N.; Baser, K.H.C.; Demirci, B.; Duman, H. Essential oils of Sideritis species of Turkey belonging to the section Empedoclia. Chem. Nat. Comp. 2004, 40, 19–23. [Google Scholar] [CrossRef]
  47. Kirimer, N.; Demirci, B.; Iscan, G.; Baser, K.H.C.; Duman, H. Composition of the essential oils of two Sideritis species from Turkey and antimicrobial activity. Chem. Nat. Comp. 2008, 44, 121–123. [Google Scholar] [CrossRef]
  48. Mahdavi, A.; Moradi, P.; Mastinu, A. Variation in terpene profiles of Thymus vulgaris in water deficit stress response. Molecules 2020, 25, 1091. [Google Scholar] [CrossRef] [Green Version]
  49. Kigathi, R.N.; Weisser, W.W.; Reichelt, M.; Gershenzon, J.; Unsicker, S.B. Plant volatile emission depends on the species composition of the neighboring plant community. BMC Plant Biol. 2019, 19, 58. [Google Scholar] [CrossRef] [PubMed]
  50. Perrino, E.V.; Valerio, F.; Gannouchi, A.; Trani, A.; Mezzapesa, G. Ecological and plant community implication on essential oils composition in useful wild officinal species: A pilot case study in Apulia (Italy). Plants 2021, 10, 574. [Google Scholar] [CrossRef]
  51. Topić, J.; Šegulja, N. Floristic and ecological characteristics of the southernmost part of Istria (Croatia). Acta Bot. Croat. 2000, 59, 179–200. [Google Scholar]
  52. Bogdanović, S.; Iveša, N.; Temunović, M.; Ljubičić, I. Plant diversity of Gornji Kamenjak (Istria, Croatia). Agric. Conspec. Sci. 2018, 83, 195–204. [Google Scholar]
  53. Milović, M.; Pandža, M.; Jasprica, N.; Tafra, D.; Krpina, V. The vascular flora of Mt Svilaja (Outer Dinarides, South Croatia). Nat. Croat. 2021, 30, 85–144. [Google Scholar] [CrossRef]
  54. Milović, M. The flora of Šibenik and its surroundings. Nat. Croat. 2002, 11, 171–223. [Google Scholar]
  55. Tafra, D.; Pandža, M.; Milović, M. Vacular flora of the town of Omiš. Nat. Croat. 2012, 21, 301–334. [Google Scholar]
  56. Vladović, D.; Ilijanić, L.J. Treći prilog flori planine Mosor (Hrvatska). Acta Bot. Croat. 1995, 54, 41–46. [Google Scholar]
  57. Štamol, V.; Marković, L.J. Prilog flori otoka Brača. Acta Bot. Croat. 1985, 44, 99–106. [Google Scholar]
  58. Trinajstić, I. Vaskularna flora otoka Hvara. Acta Bot. Croat. 1993, 52, 113–143. [Google Scholar]
  59. Bedalov, M. Prilog flori otoka Šolte. Acta Bot. Croat. 1989, 48, 121–127. [Google Scholar]
  60. Moradi, M.; Kaykhaii, M.; Ghiasvand, A.R.; Shadabi, S.; Salehinia, A. Comparison of headspace solid-phase microextraction, headspace single-drop microextraction and hydrodistillation for chemical screening of volatiles in Myrtus communis L. Phytochem. Anal. 2012, 23, 379–386. [Google Scholar] [CrossRef]
  61. Skalicka-Wozniak, K.; Los, R.; Glowniak, K.; Malm, A. Comparison of hydrodistillation and headspace solid-phase microextraction techniques for antibacterial volatile compounds from the fruits of Seseli libanotis. Nat. Prod. Commun. 2010, 5, 1427–1430. [Google Scholar] [CrossRef] [Green Version]
  62. Farouk, A.; Ali, H.; Al-Khalifa, A.R.; Mohsen, M.; Fikry, R. Aroma volatile compounds of parsley cultivated in the Kingdom of Saudi Arabia and Egypt extracted by hydrodistillation and headspace solid-phase microextraction. Int. J. Food Prop. 2017, 20, S2868–S2877. [Google Scholar] [CrossRef]
  63. Mékaoui, R.; Benkaci-Ali, F.; Alsafra, Z.; Eppe, G. Volatiles profile of the different parts of Algerian Bupleurum plantagineum Desf. by headspace solid-phase microextraction and hydrodistillation. Nat. Prod. Res. 2020, 34, 3134–3138. [Google Scholar] [CrossRef]
  64. Grecco, S.S.; Martins, E.G.A.; Girola, N.; de Figueiredo, C.R.; Matsuo, A.L.; Soares, M.G.; Bertoldo, B.C.; Sartorelli, P.; Lago, J.H.G. Chemical composition and in vitro cytotoxic effects of the essential oil from Nectandra leucantha leaves. Pharm. Biol. 2015, 53, 133–137. [Google Scholar] [CrossRef]
  65. da Silva, E.B.P.; Matsuo, A.L.; Figueiredo, C.R.; Chaves, M.H.; Sartorelli, P.; Lago, J.H.G. Chemical constituents and cytotoxic evaluation of essential oils from leaves of Porcelia macrocarpa (Annonaceae). Nat. Prod. Commun. 2013, 8, 277–279. [Google Scholar] [CrossRef] [Green Version]
  66. Noriega, P.; Guerrini, A.; Sacchetti, G.; Grandini, A.; Ankuash, E.; Manfredini, S. Chemical composition and biological activity of five essential oils from the Ecuadorian Amazon rain forest. Molecules 2019, 24, 1637. [Google Scholar] [CrossRef] [Green Version]
  67. Dahham, S.S.; Tabana, Y.M.; Iqbal, M.A.; Ahamed, M.B.K.; Ezzat, M.O.; Majid, A.S.A.; Majid, A.M.S.A. The anticancer, antioxidant and antimicrobial properties of the sesquiterpene β-caryophyllene from the essential oil of Aquilaria crassna. Molecules 2015, 20, 11808–11829. [Google Scholar] [CrossRef]
  68. Mashayekhi, H.A.; Rezaee, M.; Garmaroudi, S.S.; Montazeri, N.; Ahmadi, S.J. Rapid and sensitive determination of benzaldehyde arising from benzyl alcohol used as preservative in an injectable formulation solution using dispersive liquid–liquid microextraction followed by gas chromatography. Anal. Sci. 2011, 27, 865–868. [Google Scholar] [CrossRef] [Green Version]
  69. Almoughrabie, S.; Ngari, C.; Guillier, L.; Briandet, R.; Poulet, V.; Dubois-Brissonnet, F. Rapid assessment and prediction of the efficiency of two preservatives against S. aureus in cosmetic products using High Content Screening—Confocal Laser Scanning Microscopy. PLoS ONE 2020, 15, e0236059. [Google Scholar] [CrossRef]
  70. Trevizan, L.N.T.; do Nascimento, K.F.; Santos, J.A.; Kassuya, C.A.L.; Cardoso, C.A.L.; Vieira, M.C.; Moreira, F.M.F.; Croda, J.; Formagio, A.S.N. Anti-inflammatory, antioxidant and anti-Mycobacterium tuberculosis activity of viridiflorol: The major constituent of Allophylus edulis (A. St.-Hil., A. Juss. & Cambess.) Radlk. J. Ethnopharmacol. 2016, 192, 510–515. [Google Scholar] [PubMed]
  71. Schepetkin, I.A.; Özek, G.; Özek, T.; Kirpotina, L.N.; Khlebnikov, A.I.; Quinn, M.T. Chemical composition and immunomodulatory activity of essential oils from Rhododendron albiflorum. Molecules 2021, 26, 3652. [Google Scholar] [CrossRef] [PubMed]
  72. Trikka, F.; Michailidou, S.; Makris, A.M.; Argiriou, A. Biochemical fingerprint of Greek Sideritis spp.: Implications for potential drug discovery and advanced breeding strategies. Med. Aromat. Plants 2019, 8, 335. [Google Scholar] [CrossRef] [Green Version]
  73. Koedam, A. Volatile oil composition of Greek mountain tea (Sideritis spp.). J. Sci. Food Agric. 1986, 36, 681–684. [Google Scholar] [CrossRef]
  74. Pljevljakušić, D.; Šavikin, K.; Janković, T.; Zdunić, G.; Ristić, M.; Godjevac, D.; Konić-Ristić, A. Chemical properties of the cultivated Sideritis raeseri Boiss. & Heldr. subsp. raeseri. Food Chem. 2011, 124, 226–233. [Google Scholar]
  75. Assessment report on Sideritis scardica Griseb.; Sideritis clandestina (Bory & Chaub.) Hayek; Sideritis raeseri Boiss. & Heldr.; Sideritis syriaca L., herba. Available online: https://www.ema.europa.eu/en/documents/herbal-report/final-assessment-report-sideritis-scardica-griseb-sideritis-clandestina-bory-chaub-hayek-sideritis_en.pdf (accessed on 24 August 2021).
  76. Jerković, I.; Kranjac, M.; Marijanović, Z.; Roje, M.; Jokić, S. Chemical diversity of headspace and volatile oil composition of two brown algae (Taonia atomaria and Padina pavonica) from the Adriatic Sea. Molecules 2019, 24, 495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Jerković, I.; Cikoš, A.-M.; Babić, S.; Čižmek, L.; Bojanić, K.; Aladić, K.; Ul’yanovskii, N.V.; Kosyakov, D.S.; Lebedev, A.T.; Čož-Rakovac, R.; et al. Bioprospecting of less-polar constituents from endemic brown macroalga Fucus virsoides J. Agardh from the Adriatic Sea and targeted antioxidant effects in vitro and in vivo (zebrafish model). Mar. Drugs 2021, 19, 235. [Google Scholar] [CrossRef] [PubMed]
  78. RStudio Team. RStudio: Integrated Development for R. RStudio; PBC: Boston, MA, USA, 2021; Available online: https://www.rstudio.com/ (accessed on 24 August 2021).
  79. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.r-project.org/ (accessed on 24 August 2021).
  80. Tang, Y.; Horikoshi, M.; Li, W. ggfortify: Unified Interface to Visualize Statistical Results of Popular R Packages. R J. 2016, 8, 474–485. [Google Scholar] [CrossRef] [Green Version]
  81. Horikoshi, M.; Tang, Y. ggfortify: Data Visualization Tools for Statistical Analysis Results. 2016. Available online: https://CRAN.R-project.org/package=ggfortify (accessed on 24 August 2021).
Figure 1. Average relative percentages of different classes of headspace volatile organic compounds obtained from various populations of S. romana and S. montana using DVB/CAR/PDMS fiber.
Figure 1. Average relative percentages of different classes of headspace volatile organic compounds obtained from various populations of S. romana and S. montana using DVB/CAR/PDMS fiber.
Molecules 26 05968 g001
Figure 2. Average relative percentages of different classes of headspace volatile organic compounds obtained from various populations of S. romana and S. montana using PDMS/DVB fiber.
Figure 2. Average relative percentages of different classes of headspace volatile organic compounds obtained from various populations of S. romana and S. montana using PDMS/DVB fiber.
Molecules 26 05968 g002
Figure 3. Biplot obtained by principal component analysis of VOCs composition of investigated populations of S. romana and S. montana, based on major components with average percentages ≥ 1% in at least one of the samples, detected with DVB/CAR/PDMS fiber.
Figure 3. Biplot obtained by principal component analysis of VOCs composition of investigated populations of S. romana and S. montana, based on major components with average percentages ≥ 1% in at least one of the samples, detected with DVB/CAR/PDMS fiber.
Molecules 26 05968 g003
Figure 4. Biplot obtained by principal component analysis of VOCs composition of investigated populations of S. romana and S. montana, based on major components with average percentages ≥ 1% in at least one of the samples, detected with PDMS/DVB fiber.
Figure 4. Biplot obtained by principal component analysis of VOCs composition of investigated populations of S. romana and S. montana, based on major components with average percentages ≥ 1% in at least one of the samples, detected with PDMS/DVB fiber.
Molecules 26 05968 g004
Figure 5. Harvesting locations for investigated Sideritis species.
Figure 5. Harvesting locations for investigated Sideritis species.
Molecules 26 05968 g005
Table 1. Volatile organic compounds (VOCs) composition (%) of S. romana and S. montana obtained by HS-SPME/GC-MS with DVB/CAR/PDMS fiber.
Table 1. Volatile organic compounds (VOCs) composition (%) of S. romana and S. montana obtained by HS-SPME/GC-MS with DVB/CAR/PDMS fiber.
No.CompoundCD 1RI 2S. romana
K
S. romana
B
S. romana
M
S. romana
P
S. montana
J
S. montana
M
1Acetic acidC1<900-0.853.111.330.320.55
2Pentanal <9000.150.40----
3Hexanal <9000.070.080.520.240.020.04
4(E)-Hex-2-enal <9000.240.130.900.420.140.47
5α-ThujeneC29360.471.031.080.210.160.05
6α-PineneC39441.320.942.512.193.421.23
7Camphene 9600.070.070.010.190.07-
8BenzaldehydeC49690.440.322.580.710.370.42
9Hexanoic acid 9750.12----0.09
10Sabinene 9810.270.440.210.090.130.09
11Oct-1-en-3-ol 9830.60---0.130.46
12β-PineneC59853.692.616.678.012.811.52
13β-MyrceneC69950.800.841.270.880.910.64
14(E,E)-Hepta-2,4-dienal 1001-0.070.63-0.05-
15α-PhellandreneC710110.323.372.540.200.100.09
16α-TerpineneC810230.620.611.320.380.070.07
17p-CymeneC910320.991.815.491.320.040.05
18LimoneneC1010366.347.587.101.688.146.31
19Benzyl alcoholC1110402.963.8213.283.840.842.14
20Phenylacetaldehyde 10510.12----0.09
21γ-TerpineneC1210661.750.983.671.400.080.08
22Octan-1-olC1310751.372.223.552.500.350.29
23α-Terpinolene 10930.330.520.340.050.08-
24Linalool 1103-0.26----
25Nonanal 11070.110.070.420.19-0.09
26Benzeneethanol 11170.170.13--0.04-
27α-Campholenal 11320.07-----
28trans-Pinocarveol 3 1145-0.360.770.40--
29Pinocarvone 11690.200.230.710.31--
304-Terpineol 11820.600.950.890.150.050.07
31Cryptone 1191-0.20----
32α-TerpineolC1411950.503.851.250.020.04-
33MyrtenalC1511990.390.541.030.47--
34Decanal 12080.160.10-0.23--
35Verbenone 12120.110.15----
362-Phenoxyethanol 12230.150.06----
37Carvone 12490.12-----
38Bicycloelemene 13420.750.53-0.400.220.41
39α-Cubebene 13530.12---0.560.56
40α-CopaeneC1613750.18---1.672.02
41Isoledene 13770.370.29-0.230.34-
42Longifolene 13950.560.34-0.280.37-
43α-GurjuneneC1714121.170.950.610.850.120.41
44AristoleneC1814221.080.91----
45trans-CaryophylleneC191423--3.6023.8811.896.56
46Calarene 14360.330.31-0.34-0.19
47AromadendreneC2014381.111.05-0.960.232.87
48AlloaromadendreneC2114437.326.692.356.270.760.90
49Selina-5,11-diene 3C2214472.522.28-1.91--
50γ-Muurolene 14550.990.71-0.15--
51trans-β-FarneseneC2314615.1212.791.005.584.833.56
52IsocaryophylleneC2414714.702.39-3.22--
53γ-Gurjunene 14760.930.62-0.45--
54α-AmorpheneC2514800.58---2.843.31
55Germacrene DC2614841.500.06--23.2317.04
56trans-β-Ionone 14890.310.180.370.22--
57β-CadineneC2714942.191.68-1.240.981.47
58LedeneC281497--8.20-2.27-
59BicyclogermacreneC29149823.0420.48-15.164.1911.76
60α-MuuroleneC3015020.610.18--1.782.93
61γ-CadineneC3115170.58---3.333.37
62δ-CadineneC3215271.840.15--7.578.85
63Dihydroactinidiolide 15320.410.410.600.330.110.35
64Cadina-1,4-dieneC3315360.20---0.941.24
65α-CadineneC3415410.24---1.411.70
66α-CalacoreneC3515470.470.04--1.051.36
67SpathulenolC3615800.581.30-0.520.150.47
68ViridiflorolC3715950.420.3610.290.200.080.21
69α-Cadinol 16580.34-----
Total identified [%] 86.1890.2988.8789.6089.2886.38
1 CD: Compound designation of the major components (average percentages ≥ 1.0% in at least one of the samples), which were included in the principal component analysis. 2 RI: Retention index determined relative to a homologous series of n-alkanes (C9–C25) on a HP-5MS column; 3 Tentatively identified.
Table 2. Volatile organic compounds (VOCs) composition (%) of S. romana and S. montana obtained by HS-SPME/GC-MS with PDMS/DVB fiber.
Table 2. Volatile organic compounds (VOCs) composition (%) of S. romana and S. montana obtained by HS-SPME/GC-MS with PDMS/DVB fiber.
No.CompoundCD 1RI 2S. romana
K
S. romana
B
S. romana
M
S. romana
P
S. montana
J
S. montana
M
1Pentanal <9000.02--0.03--
2Hexanal <9000.020.030.120.100.010.02
3(E)-Hex-2-enal <9000.020.030.100.060.030.05
4α-Thujene 9360.550.410.790.580.070.04
5α-PineneC1′9442.092.313.912.681.691.14
6Camphene 9600.050.040.100.060.02-
7Benzaldehyde 9690.060.080.400.150.040.10
8Sabinene 9810.560.350.460.430.120.10
9β-PineneC2′9856.788.0113.2610.351.801.46
10β-Myrcene 9950.550.400.990.800.280.20
11(E,E)-Hepta-2,4-dienal 1001-0.050.120.070.020.03
12α-Phellandrene 10110.150.830.270.220.02-
13α-TerpineneC3′10230.53-1.151.01--
14p-CymeneC4′10320.810.163.101.57--
15LimoneneC5′10363.861.973.312.273.742.91
16Benzyl alcoholC6′10400.630.962.551.120.160.55
17(Z)-β-ocymene 10430.03-----
18Phenylacetaldehyde 10510.03----0.02
19(E)-β-ocymene 10540.02-0.11---
20γ-TerpineneC7′10662.530.118.964.49--
21Octan-1-olC8′10751.171.332.462.240.120.16
22α-Terpinolene 10930.110.070.190.140.01-
23Linalool 1103-0.860.250.22--
24Nonanal 11070.03-----
25Benzeneethanol 11170.02-----
26trans-Pinocarveol 3 1145-0.050.35---
27Pinocarvone 11690.100.040.26---
284-Terpineol 11820.060.030.350.15--
29α-Terpineol 11950.220.580.41---
30Myrtenal 11990.150.070.700.110.02-
31Decanal 1208--0.190.11--
32Verbenone 12120.04-0.10---
332-Phenoxyethanol 12230.01-----
34BicycloelemeneC9′13422.182.400.431.730.340.68
35α-Cubebene 1353----0.040.05
36α-CopaeneC10′13750.07---1.331.55
37α-Gurjunene 1412--0.79---
38trans-CaryophylleneC11′14230.540.916.352.3015.496.65
39Aromadendrene 14380.25-----
40Alloaromadendrene 1443---0.250.060.14
41Selina-5,11-diene 3 1447-0.26----
42γ-Muurolene 1455----0.090.32
43trans-β-FarneseneC12′14614.887.713.554.904.022.17
44IsocaryophylleneC13′14717.106.85-5.17--
45α-Amorphene 1480----0.19-
46Germacrene DC14′148412.182.161.771.2253.6351.08
47trans-β-Ionone 1489-0.270.390.16--
48β-Cadinene 1494----0.010.01
49BicyclogermacreneC15′149847.3257.5010.2148.089.5418.03
50α-Muurolene 1502----0.07-
51γ-Cadinene 1517----0.250.32
52δ-CadineneC16′15270.10---0.911.32
53Dihydroactinidiolide 15320.170.260.420.220.120.20
54Cadina-1,4-diene 1536-----0.05
55Spathulenol 1580----0.34-
56ViridiflorolC17′1595--24.19---
Total identified [%] 95.9997.0993.0692.9994.5889.35
1 CD: Compound designation of the major components (average percentages ≥ 1.0% in at least one of the samples), which were included in the principal component analysis. 2 RI: Retention index determined relative to a homologous series of n-alkanes (C9–C25) on a HP-5MS column; 3 Tentatively identified.
Table 3. Composition (%) of major identified VOCs of S. romana and S. montana obtained by HS-SPME/GC-MS with DVB/CAR/PDMS and PDMS/DVB fiber.
Table 3. Composition (%) of major identified VOCs of S. romana and S. montana obtained by HS-SPME/GC-MS with DVB/CAR/PDMS and PDMS/DVB fiber.
No.CompoundCD 1(A 2/B 3)S. romana K
(A 2/B 3)
S. romana B
(A 2/B 3)
S. romana M
(A 2/B 3)
S. romana P
(A 2/B 3)
S. montana J
(A 2/B 3)
S. montana M
(A 2/B 3)
1Acetic acidC1/--/-0.85/-3.11/-1.33/-0.32/-0.55/-
2α-ThujeneC2/-0.47/0.551.03/0.411.08/0.790.21/0.580.16/0.070.05/0.04
3α-PineneC3/C1′1.32/2.090.94/2.312.51/3.912.19/2.683.42/1.691.23/1.14
4BenzaldehydeC4/-0.44/0.060.32/0.082.58/0.400.71/0.150.37/0.040.42/0.10
5β-PineneC5/C2′3.69/6.782.61/8.016.67/13.268.01/10.352.81/1.801.52/1.46
6β-MyrceneC6/-0.80/0.550.84/0.401.27/0.990.88/0.800.91/0.280.64/0.20
7α-PhellandreneC7/-0.32/0.153.37/0.832.54/0.270.20/0.220.10/0.020.09/-
8α-TerpineneC8/C3′0.62/0.530.61/-1.32/1.150.38/1.010.07/-0.07/-
9p-CymeneC9/C4′0.99/0.811.81/0.165.49/3.101.32/1.570.04/-0.05/-
10LimoneneC10/C5′6.34/3.867.58/1.977.10/3.311.68/2.278.14/3.746.31/2.91
11Benzyl alcoholC11/C6′2.96/0.633.82/0.9613.28/2.553.84/1.120.84/0.162.14/0.55
12γ-TerpineneC12/C7′1.75/2.530.98/0.113.67/8.961.40/4.490.08/-0.08/-
13Octan-1-olC13/C8′1.37/1.172.22/1.333.55/2.462.50/2.240.35/0.120.29/0.16
14α-TerpineolC14/-0.50/0.223.85/0.581.25/0.410.02/-0.04/--/-
15MyrtenalC15/-0.39/0.150.54/0.071.03/0.700.47/0.11-/0.02-/-
16Bicycloelemene-/C9′0.75/2.180.53/2.40-/0.430.40/1.730.22/0.340.41/0.68
17α-CopaeneC16/C10′0.18/0.07-/--/--/-1.67/1.332.02/1.55
18α-GurjuneneC17/-1.17/-0.95/-0.61/0.790.85/-0.12/-0.41/-
19AristoleneC18/-1.08/-0.91/--/--/--/--/-
20trans-CaryophylleneC19/C11′-/0.54-/0.913.60/6.3523.88/2.3011.89/15.496.56/6.65
21AromadendreneC20/-1.11/0.251.05/--/-0.96/-0.23/-2.87/-
22AlloaromadendreneC21/-7.32/-6.69/-2.35/-6.27/0.250.76/0.060.90/0.14
23Selina-5,11-diene 4C22/-2.52/-2.28/0.26-/-1.91/--/--/-
24trans-β-FarneseneC23/C12′5.12/4.8812.79/7.711.00/3.555.58/4.904.83/4.023.56/2.17
25IsocaryophylleneC24/C13′4.70/7.102.39/6.85-/-3.22/5.17-/--/-
26α-AmorpheneC25/-0.58/--/--/--/-2.84/0.193.31/-
27Germacrene DC26/C14′1.50/12.180.06/2.16-/1.77-/1.2223.23/53.6317.04/51.08
28β-CadineneC27/-2.19/-1.68/--/-1.24/-0.98/0.011.47/0.01
29LedeneC28-/-/--/-8.20/--/-2.27/--/-
30BicyclogermacreneC29/C15′23.04/47.3220.48/57.50-/10.2115.16/48.084.19/9.5411.76/18.03
31α-MuuroleneC30/-0.61/-0.18/--/--/-1.78/0.072.93/-
32γ-CadineneC31/-0.58/--/--/--/-3.33/0.253.37/0.32
33δ-CadineneC32/C16′1.84/0.100.15/--/--/-7.57/0.918.85/1.32
34Cadina-1,4-dieneC33/-0.20/--/--/--/-0.94/-1.24/0.05
35α-CadineneC34/-0.24/--/--/--/-1.41/-1.70/-
36α-CalacoreneC35/-0.47/-0.04/--/--/-1.05/-1.36/-
37SpathulenolC36/-0.58/-1.30/--/-0.52/-0.15/0.340.47/-
38ViridiflorolC37/C17′0.42/-0.36/-10.29/24.190.20/-0.08/-0.21/-
1 CD: Compound designation of the major components (average percentages ≥ 1.0% in at least one of the samples), which were included in the principal component analysis; 2 A: Compound extracted using DVB/CAR/PDMS fiber; 3 B: Compound extracted using PDMS/DVB fiber; 4 Tentatively identified.
Table 4. Sample codes, voucher numbers, harvesting locations, geographic coordinates and climatic affiliations according to Köppen climate classification for investigated Sideritis species.
Table 4. Sample codes, voucher numbers, harvesting locations, geographic coordinates and climatic affiliations according to Köppen climate classification for investigated Sideritis species.
SpeciesSample CodeVoucher No.LocationDate of CollectionLatitudeLongitudeKöppen Climate
Classification 1
S. romanaS. romana K18 021Istria,
Kamenjak
13 June 202044°47′46.15″ N13°54′25.96″ ECfa 2
S. romana B18 023Dalmatia,
Blato
21 June 202043°10′12.71″ N17°11′45.75″ ECsa 3
S. romana M18 020Dalmatia, Morinje Bay7 June 202043°40′44.79″ N15°57′41.25″ ECsa 3
S. romana P18 022Dalmatia,
Brač, Pražnica
20 June 202043°19′1.3″ N16°40′24.45″ ECsb 4
S. montanaS. montana J18 010Dalmatia,
Ježević
6 June 202043°55′1.51″ N16°28′8.41″ ECfb 5
S. montana M18 011Dalmatia, Mosor, Gornje Sitno22 June 202043°31′13.56″ N16°36′10.44″ ECfa 2
1 Data obtained from Šegota and Filipčić [45]; 2 Cfa: temperate humid climate with hot summer; 3 Csa: Mediterranean climate with hot summer; 4 Csb: Mediterranean climate with warm summer; 5 Cfb: temperate humid climate with warm summer.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marić, T.; Friščić, M.; Marijanović, Z.; Maleš, Ž.; Jerković, I. Comparison of Volatile Organic Compounds of Sideritis romana L. and Sideritis montana L. from Croatia. Molecules 2021, 26, 5968. https://doi.org/10.3390/molecules26195968

AMA Style

Marić T, Friščić M, Marijanović Z, Maleš Ž, Jerković I. Comparison of Volatile Organic Compounds of Sideritis romana L. and Sideritis montana L. from Croatia. Molecules. 2021; 26(19):5968. https://doi.org/10.3390/molecules26195968

Chicago/Turabian Style

Marić, Tihana, Maja Friščić, Zvonimir Marijanović, Željan Maleš, and Igor Jerković. 2021. "Comparison of Volatile Organic Compounds of Sideritis romana L. and Sideritis montana L. from Croatia" Molecules 26, no. 19: 5968. https://doi.org/10.3390/molecules26195968

Article Metrics

Back to TopTop