Next Article in Journal
Raw Meat Contaminated with Cephalosporin-Resistant Enterobacterales as a Potential Source of Human Home Exposure to Multidrug-Resistant Bacteria
Previous Article in Journal
Insights into the Allosteric Effect of SENP1 Q597A Mutation on the Hydrolytic Reaction of SUMO1 via an Integrated Computational Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Profiling of Essential Oils from the Leaves of Pistacia lentiscus Collected in the Algerian Region of Tizi-Ouzou: Evidence of Chemical Variations Associated with Climatic Contrasts between Littoral and Mountain Samples

1
BIOPI-UPJV Laboratory UMRT BioEcoAgro INRAE1158, SFR Condorcet FR CNRS 3417, UFR de Sciences, Université de Picardie Jules Verne, 33 Rue Saint Leu, 80000 Amiens, France
2
Laboratory of Natural Resources, University Mouloud Mammeri of Tizi-Ouzou, Tizi Ouzou 15000, Algeria
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(13), 4148; https://doi.org/10.3390/molecules27134148
Submission received: 19 April 2022 / Revised: 22 June 2022 / Accepted: 23 June 2022 / Published: 28 June 2022

Abstract

:
Leaves of Pistacia lentiscus were collected from two Algerian sites in the mountains and the littoral of the Tizi-Ouzou region. The harvest was conducted in four consecutive seasons on the same selected set of trees. Essential oils (EOs) were extracted by hydrodistillation; then, they were analyzed by gas chromatography coupled mass spectrometry (GC-MS). Forty-seven constituents could be detected and quantified, including α-pinene (2–13%), β-caryophyllene (8–25%), β-myrcene (0.3–19%), bornyl acetate (0.8–7%), δ-cadinene (3–8%), bisabolol (1–9%), β-pinene (0.9–7%), caryophyllene oxide (4–9%), and α-cadinol (3–11%). Antioxidant (AOx) activities of the EOs were assessed by ferric reducing antioxidant power (FRAP), 2,2-diphenyl-1-picrylhydrazyl (DPPH), and 2,2′-azino-bis (3-ethylbenzothiazoline)-6-sulfonic acid (ABTS) assays. Significant differences in EO composition and AOx activities appeared dependent on the season and the site. Variations of AOx activities were significant for the FRAP and ABTS tests but not for DPPH. Characterization of the leaf fatty acyl (FA) profiles was performed by GC-MS. Variability appeared according to season and altitude. Polyunsaturated fatty acids levels were high (27–55%) at the coldest date and place. The levels of linolenic acyl in the leaves were significantly correlated with bisabolol levels in the EOs (Spearman’s correlation coefficient: 0.818). Such results will be useful for the sustainable local valorization of wild P. lentiscus. These data also open new routes for further studies on terpenoid biosynthesis using correlation networks and fluxomic approaches.

Graphical Abstract

1. Introduction

Belonging to the Anacardiaceae family, Pistacia lentiscus (P. lentiscus) is one of at least 11 plant species of the Pistacia genus. Commonly called “lentisk” or “mastic tree” (in reference to its resin or so-called mastic), P. lentiscus is an evergreen shrub, largely distributed in Mediterranean countries [1]. It is a branched perennial thermophile shrub, measuring one to five meters in height. Its leaves are pinnate and evergreen. It is a dioecious species with distinct male and female plants harboring unisexual flowers. These flowers are approximately three millimeters wide and appear in spring (from March to May). They are aromatic, either male or female, and occur in small clusters. The female flowers are yellowish green. The male flowers are dark red. On female trees, the fruit is a rounded mono-sperm drupe of four to six millimeters in diameter. The fruit is first green, then red, and becomes blackish when ripe in autumn [2]. A resin, also called mastic, is the best-known product of this plant. It is an aromatic and resinous substance that seeps from the trunk and the main branches [3]. The traditional use of different parts of the plant for therapeutic applications makes lentisk considered as a medicinal plant [4]. People in rural areas use the infusion of its leaves to treat gastric ulcers [5], respiratory disorders [6], and diabetes for its hypoglycemic effect [7]. They also use decoctions of lentisk aerial parts to treat tooth ache and throat infections [8]. The results of an ethnopharmacological study on lentisk, carried out in Algeria, have led to various therapeutic indications of this plant in the country [9]. The main components and pharmacological properties of the species were reviewed by Nahida (i.e., antimicrobial, hepatoprotective, antioxidant, and anticancer activities) [10] and Bozorgi (i.e., antimicrobial activity) [11]. The Pistacia species show variations in morphological, anatomical, chemical, and molecular traits across geographic origins and sexes [12].
Essential oils (EOs) have been produced from several parts of the Pistacia species [13]. Extracted from leaves, resin, mature or immature fruits, galls, buds, twigs, or flowers, lentisk EOs are generally a limpid liquid with a yellow color and an intense and piercing aromatic odor [13].
Different EOs from Pistacia lentiscus leaves have been produced by hydrodistillation [14]. Their chemical composition greatly depends on their origins (i.e., Algeria, Morocco, Tunisia, Italy, Spain, Egypt, Greece, or France). In Algeria, P. lentiscus EOs are mainly characterized by terpinen-4-ol (17%) and limonene (34%) [15,16]. Myrcene (39.2% [17]; 25.3% [18]) is the main compound reported in EOs from Morocco. In Tunisian and Italian EOs, α-pinene represented 9–18% of the total components [19,20]. Other studies reported the presence of monoterpene hydrocarbons [21,22,23] or sesquiterpene hydrocarbons [20,24] as the predominant fraction in lentisk EOs. The change in the chemical composition of the oils has been attributed (i) to the wide range of habitats in which lentisk grows [25]; (ii) to environmental factors; (iii) to seasons [26] and phenological stages [27]; (iv) to gender [2,28].
Lentisk EOs and derivatives have been used in perfumery [29] (e.g., Perfume “Accord Boisé” Minalys™ from the Olfactis Cie), in the food industry as flavoring in drinks [30,31] (e.g., liqueur Callicounis™ and Skinos™), and chewing gum [29] (e.g., Mastic Gum Elma™). P. lentiscus EOs have interesting biological activities. Indeed, antimicrobial effects against Helicobacter pylori, Staphylococcus aureus, Sarcina lutea, Escherichia coli, Torulopsis glabrata, and Listeria monocytogenes have been reported [21,29]. The use of P. lentiscus EOs as a natural fumigant against Lepidoptera, such as Ectomyelois ceratoniae and Ephesia kuehniella, had significant insecticidal effects [32]. The oil shows different levels of antioxidant properties [33].
These previous studies described various potential applications of P. lentiscus essential oils. Nevertheless, each use was based on the fine composition of essential oils and large variability occurred. It is also known that biological activity depends on the chemical composition of the EOs [34,35]. Most of the time, the main components are important, but sometimes synergic effects due to the different, eventually minor, compounds can lead to surprising enhanced (or decreased) biological activities.
Thus, in order to evaluate the potential of the natural P. lentiscus resources from the wide Algerian region of Tizi-Ouzou, it is of primary necessity to describe the profile of the essential oils from the leaves of wild plants. It is also needed to determine the variations of such profiles, and, if possible, to determine the main factors affecting such a profile.
This work aimed to explore the composition of wild P. lentiscus EOs extracted from leaves. Sample collection was conducted from two contrasting sites in the Tizi-Ouzou region (i.e., mountain and littoral); on the same selected tree from the bush, sampling was conducted four times a year to cover the four seasons (and yearly phenological stages of this perennial plant). We thus studied the variability of the terpenoid content in P. lentiscus and investigated the impact of this variability on its antioxidant activities as a function of both the harvesting period and altitude. In addition, the profiling of the fatty acyls from leaf lipids was performed. It allowed for proving that the main variations in the EOs’ composition were associated with the adaptation of the plants during the colder season in the Tizi-Ouzou Mountains.

2. Results and Discussion

2.1. Qualitative and Semi-Quantitative Composition of the Essential Oils

The different essential oils (EOs) were analyzed by GC-MS. Based on the previous literature and databases, 47 chromatographic peaks could be annotated (Table 1). These compounds corresponded to chromatographic peaks representing from 88% to 97% of the total area of all peaks. These oils were characterized by their richness in sesquiterpene hydrocarbons. This volatile fraction was the most abundant (29–46%) in almost all the different EOs (Table 1). Monoterpene hydrocarbons were present (5–42%) in all the samples. Oxygenated sesquiterpenes and oxygenated monoterpenes were in 13–37% and 3–9%, respectively. Within all samples, the major compounds detected were: α-pinene (2–13%), β-caryophyllene (8–25%), β-myrcene (0.1–19%), bornyl acetate (0.8–7%), δ-cadinene (3–8%), bisabolol (1–9%), β-pinene (0.9–7%), caryophyllene oxide (4–9%), and α-cadinol (3–11%).
Recent data from different Mediterranean countries have noted such a wide spectrum of chemical variability involving both the main compounds and the total amounts of the terpene classes [36,37]. For terpene classes, monoterpene hydrocarbons usually represented the main fraction: 75% in Egypt [37], 68% in Greece [3], 59% in Tunisia [24], and 56% in Italy [36]. However, in Tunisia, lentisk EOs were equally rich in monoterpene hydrocarbons (41%) and sesquiterpene hydrocarbons (40%) [19]. Lo-Presti et al. [20] and Llorens et al. [31] showed the prominence of sesquiterpenes in EOs from Italy (47%) and Spain (39%). Beyond terpene class variations, terpene levels varied as well. Thus, depending on the soil quality, different EO chemotypes have been described. In Spain, Llorens et al. [31], detected the presence of germacrene (10–17%), β-caryophyllene (11–12%), myrcene (1–11%), and terpinen-4-ol (5–6%) from calcicolous and siliceous soils. In Italy, Congiu et al. [38] detected β-caryophyllene (31%), germacrene (12%), and δ-cadinene (6%) as the main terpenes, while α-pinene and terpinen-4-ol were the chemotypes described in the (also Italian) EOs analyzed by Negro et al. [29]. EOs from Morocco varied in chemical composition; some were rich in myrcene (33–39%) [17,18,39], while others were mainly characterized by their 24% α-pinene level [40]. A number of other chemotypes rich in car-3-ene were also reported in Morocco [41]. The lentisk EOs in Tunisia were characterized by their high level of α-pinene (17%), terpinen-4-ol (12%), and δ-terpinene (9%) [24]. In addition, in Tunisia, germacrene D (11%), α-pinene (9%), limonene (8%), β-caryophyllene (8%), and δ-cadinene (8%) were found to be major compounds by Aissi et al. [19]. The (leaf) EOs from Egypt were described as rich in α-phellandrene [37]. It is thus clear that EOs from the leaves of Pistacia lentiscus presented variable compositions depending on the countries where the plants grew. Monoterpenes (i.e., α-pinene, terpineol, myrcene, limonene, and terpinen-4-ol) and sesquiterpenes (i.e., β-caryophyllene, δ-cadinene, and germacrene D) were the most common compounds. However, their levels greatly varied depending on the origin of the samples. The main factor contributing to this chemo-variability was attributed to the environmental conditions in the habitats of the Pistacia lentiscus species. Few or no data exist concerning genetic traits associated with the EO profiles. If P. lentiscus is an allogamous (dioic) species, we can assume that chemotypes were mostly (and naturally) selected essentially based on pedoclimatic traits. In our study, according to the harvest time (season) and place (altitude), various main chemotypes were detected. Except for the precise percentages, some of our main chemotypes correspond to other profiles previously described: caryophyllene [31,38], δ-cadinene [19], myrcene [18,39], and α-pinene [22]. Variations in EO profiles can be of interest. Nevertheless, it is important to determine the origins of variations before expecting a standardized production and a further economical use of such natural products, hence our interest in this littoral and mountain study.

2.2. Quantitative Variations of the Major Compounds and Terpenic Classes

The relative abundance of each compound in the different essential oils (EOs) was determined based on chromatographic peak areas measured in the total ion current (TIC). As previously described [2], each relative abundance was then expressed as a percentage of the peak total area.
Table 1 presents the mean contents (and standard deviations) of each EO based on their origin: season and geographical site. Variations among sample origins clearly occurred. For the main chromatographic peaks, major differences can be observed (Table 1). As such, the highest percentages of α-pinene (12.6%) and β-pinene (7.1%) were, respectively, found in spring and summer at the mountain site. The highest amount of β-myrcene (19%) was detected in the mountains in autumn. In all EOs, the levels of β-caryophyllene exceeded 8% at both sites and harvest periods from 8.8% (winter, mountain) to 25.4% (summer, coast). δ-Cadinene (8.5%) and bisabolol (9.7%) harbored their highest levels in the EOs collected in winter in the mountains. In addition to these compounds, other molecules (also formatted in bold in Table 1) presented, at least once, a level superior to 5%, e.g., bornyl acetate (0.8–7.0%), caryophyllene oxide (4.5–9.9%), and α-cadinol (3.7–11.9%).
A two-factor analysis of a variance ANOVA test was performed with all the values to evaluate the effects of site altitude (Factor 1), season (Factor 2), and the possible interaction between both factors and each molecule. The significance of the effects are presented in Table 1 (column: Altitude * Season Interaction) with the following symbols: *** p-value < 0.001; ** p-value < 0.01; * p-value < 0.05. Most of the time, altitude and/or season had a significant (p < 0.05) or a highly significant p < 0.01 to p < 0.001) effect on the variation of the EO component levels. The interaction between altitude and season (Table 1, column: Altitude * Season Interaction) had a highly significant effect (p < 0.001) for most of the molecules found in the different essential oils.
A Tukey test (post hoc) was applied for all the data. Significant differences (p < 0.05) between the mean levels of each compound are reported in each row of Table 1 as different letters. Homogeneous groups of the EO contents appeared among the eight altitudinal and seasonal conditions. These data confirm that differences and similarities occurred in the chemical profiles of our EOs. However, these data were not clear enough to determine the main rules on the implied variations in the individual compounds. The same analysis performed on the different terpene classes also showed variations over the four distinct classes of compounds (i.e., Table 1: monoterpene hydrocarbons (MHs); sesquiterpene hydrocarbons (SHs); oxygenated monoterpenes (OMs); oxygenated sesquiterpenes (OS)). MHs were higher in the mountains (28–42.8%) than in the littoral (6.4–27.1%), except for during winter (5.8%). For both sites, the highest rate of MHs was observed in autumn. MH levels in the mountains in autumn (42.8%) were one-and-a-half times the level of MHs from littoral samples in autumn (27.1%). For each site, the differences between levels of MHs were significant (p < 0.05) according to the season, except between spring and summer for both sites. Contrary to MHs, sesquiterpenes hydrocarbons (SHs) were higher (p < 0.05) in the littoral area (37.9–46.5%) than in the mountain (29.1–33.8%). SH levels were higher in the winter (mountain) and in the summer (littoral). Whatever the site or altitude, OS levels were higher in winter (mountain, 37.0%; littoral, 30.3%). OMs were present in smaller quantities than those of MH-SH-OS for all samples. Their variations were not significant according to the season and the altitude. For the four compound classes, the interaction between altitude and season (Table 1, Altitude * Season Interactions) had a significant (p < 0.05) or very significant (p < 0.001 for MH and SH) influence on the average levels. According to these results, we can conclude that the two parameters studied (i.e., altitude and harvesting season) contributed significantly to the variation in the chemical composition of the essential oils obtained from leaves of Pistacia lentiscus collected in the Algerian region of Tizi-Ouzou. Nevertheless, within all the variability detected for the composition of the different EOs, it was necessary (i) to determine the main traits of each oil and (ii) to compare them with other samples previously described in the literature. As such, the major EO constituents that were detected in lentisk leaves harvested at the vegetative overwintering stage (February in Turkey) were previously described as terpinen-4-ol (29%; [42]) and β-myrcene (39%; [17]). For leaves collected at the flowering stage, in different places (i.e., countries), various EO chemotypes have also been described. The major compounds described were: α-pinene (25%), terpinen-4-ol (21%), δ-cadinene (11%), and β-caryophyllene (19%) by Negro et al. [29]; germacrene D (22%), α-pinene (17%), and sabinene (15%) by El Idrissi et al. [43]; tricyclene (7%), terpinen-4-ol (7%), sabinene (6%), and caryophyllene (6%) by Haloui et al. [23]. At a late fruiting stage (i.e., ~October), Mecherara et al. [16] detected the predominance of limonene (43%), α-pinene (34%), and myrcene (33.1 %) in EOs from lentisk leaves collected in Algeria. However, Ismail et al., (2012) in Tunisia reported the presence of α–pinene (20%), limonene (15%), β-pinene (9%), and germacrene D (8%) [22]. At early fruiting (i.e., ~June), the main EO compounds from lentisk leaves growing in the west Algerian region of Oran were terpinen-4-ol (41%), α-terpineol (7%), α–pinene (9%), limonene (9%), β-myrcene (10%), and p-cymene (8%) [44]. When comparing the composition of our EOs with those previously described in the literature and according to the harvest period, some similarities and differences can be raised. We note that some compounds, such as α-pinene, were similar at the flowering and early fruiting stages as was the case of Negro et al. [29], El Idrissi et al. [43], Haloui et al. [23], and Hamiani et al. [44]; similarities also occurred with β–myrcene in the late fruiting stage as described by Mecherara et al. [16], and β-pinene in the early fruiting stage was a main compound as shown by Ismail et al. [22].
However, our results largely differ from those of Zaouali et al. [2], who studied variations in the EOs’ chemical composition according to the plant sex for different harvest periods in Tunisia. In female trees, α-limonene was mostly detected in leaves (26–29%) sampled at the early fruiting (August) and late fruiting (October) stages. Germacrene-D was dominant in leaves (20%) at the flowering and early fruiting stages, and δ-cadinene (15.6%) reached its highest content in leaves in the flowering (March) stage. In male trees, the highest amounts of germacrene-D (13%) were recorded in leaves at the flowering stage. Oppositely, for the variations in our results according to the (limited) literature regarding lentisk EOs, the influence of ecological factors may support an important part of the explanation. In other Mediterranean evergreen plant species, the influence of ecological factors on the chemical composition of EOs extracted from leaves has been clearly pointed out in Juniperus oxycedrus [45], Rosmarinus officinalis [46], Myrtus communis [47], and Pinus halepensis [48]. However, for P. lentiscus, no detailed information has until now been available to describe, and maybe explain, the precise relationship between environmental factors and the variations in EO composition. To our knowledge, only one unique study described the influence of the altitudinal gradient on the chemo-variation of P. lentiscus EO composition [26]. The samples were collected in Algeria, but with only one harvest period (September 2008). Data from our multisite/multiseason results are in accordance with Ait Said et al. [26]. The relationship between altitude and the chemotypes were commonly characterized by caryophyllene and by significant levels of monoterpene and sesquiterpene hydrocarbons.

2.3. Multivariate Analysis of All the Components Detected in the EOs

In order to better show the differentiation of the chemical compositions between the different EOs, a principal component analysis (PCA) was carried out on the 47 major compounds annotated above (Figure 1). The two axes of the PCA show 78.6% of the variation. The score plot (Figure 1A) on Axis 1 separates the samples of the two groups according to the sampling site. In each group, the samples for seasons are separated based on Axis 2 of the PCA. Winter and autumn appear mainly separated by Component 2. The mountain group occupies a larger space than the littoral one (i.e., distances between the four seasons were higher for the mountain EOs). This indicates that the variations in the compositions according to the harvesting season were higher at the mountain than at the coastal site. The loading plot (Figure 1B) shows that the most discriminating compounds were the “major” compounds, which appeared at least once at a level greater than 5%. Observations from the littoral samples, regardless of the season, were characterized by high quantities of β-caryophyllene (Group 1). The highest proportion of β-caryophyllene was found in the summer. The winter samples from the littoral group were characterized, in addition to β-caryophyllene, by significant contents α-cadinol and caryophyllene oxide (Group 2). EOs from the mountain site contained high amounts of δ-cadinene, bisabolol, and bornyl acetate in winter (Group 3). Samples collected in autumn had a high proportion of β-myrcene followed by α-pinene (Group 4), mainly from the mountain samples. The samples harvested in spring and summer at the mountain were characterized by high amounts of α-pinene and β-pinene (Group 5).
According to the score and loading plots, the compositions of the EOs analyzed could be subdivided into five sub-chemotypes: Groups 1, 2, 3, 4, and 5.
A difference in altitude is usually accompanied by changes in a range of environmental conditions such as temperature, water precipitation, wind exposure, sunlight intensity, UV radiation, and air humidity [49]. Climatic conditions at (relatively) high altitudes are mainly lower average temperatures, increased thermal amplitude between day and night, and higher light intensity. Such environmental conditions cause plants to change their morphology, physiology, and productivity in order to protect themselves and adapt to such stressful conditions [49,50].
Abreu and Mazzafera [51] reported that lower temperatures at high altitudes can induce an overproduction of phytochemicals (Hypericum brasiliense Choisy). Previous studies have suggested that increasing altitude and subsequent changes in solar radiation and temperature in plant habitats can strongly correlate with the content of secondary metabolites, especially phenolics, due to the fact of their protective function against oxidative damage [52]. In order to evaluate a potential link between our EOs’ profiles and their antioxidant properties, further studies were performed with a primary focus on antioxidant activity measurements.

2.4. Antioxidant Activity

The antioxidant activities (AOx) of the different essential oils, assessed by the ferric reducing power (FRAP), free radical scavenging activity (DPPH), and the ABTS radical cation reduction tests are presented in Table 2.
According to the origin of the essential oil (i.e., harvesting season and altitude), ANOVA analysis showed that the altitude and season interaction had a significant effect (p < 0.001) on variations in antioxidant activity. Significant variations (p < 0.05) in the different antioxidant activities can be observed using the post hoc test (Table 2). The FRAP test always led to higher values than the ABTS and the DPPH tests. As compared with the littoral samples, the highest antioxidant (AOx) effect was measured for the mountain samples among the three tests. The ferric reducing capacity (FRAP) of these EOs varied from 9.6 (spring) to 22.3 mgTE/gEO (summer) for the mountain site and from 6.4 (winter) to 15.7 mgTE/gEO (summer) for the littoral samples. Four conditions could be statistically separated (annotation as “a, b, c, and d” for the post hoc test in in Table 2, column: FRAP). Some of the highest FRAP values were recorded for samples collected during the summer period at both sites. However, they were higher at the mountain site (22.3 mgTE/gEO versus 15.7 for littoral). The antioxidant capacity determined by the two other tests, DPPH and ABTS, was quite comparable (but lower than FRAP). The extreme ABTS values were obtained from the mountain site. They were between 0.09 (summer) and 0.32 mgTE/gEO (spring). The post hoc test separated five significantly different homogeneous groups (annotated as “a, b, c, d, and e” in Table 2, column: ABTS). For the DPPH values, levels between 0.06 (winter, littoral) and 0.44 mgTE/gEO (spring, mountain) were found. Nevertheless, in the case of DPPH activity, the post hoc test distinguished significant differences (i.e., “a, b, c, and d”) between all the sample origins (altitude and season).
The antioxidant values of vitamin C (used as a positive control for validating the quantification method) showed that 1 g of Trolox reacted as ~1 g of vitamin C (i.e., 1.08, 0.89, and 1.34 gTE/g vitamin C for the FRAP, ABTS, and DPPH, respectively; Table 2). This confirmed that the measurements of antioxidant activity were correct.
According to the best knowledge of the authors, no report on this plant is available; therefore, comparisons were made with other species from the same genus. Studies on the antioxidant activity of essential oils extracted from the Pistacia species are not as numerous. Barra et al. [53] showed a low antiradical capacity of essential oils obtained from P. lentiscus leaves collected from different origins and harvesting periods. In addition, Bouyahya et al., showed fairly good antioxidant effects (ABTS and DPPH tests) by the EOs from whole aerial parts of P. lentiscus harvested in Morocco [39]. In addition, EOs from P. atlantica [13] presented a low capacity for free radical neutralization and an interesting reducing power. On the other hand, Bampouli et al. [54] showed that P. lentiscus var. chia EOs did not exhibit antioxidant activity with the DPPH test.
An exploration of the potential correlations between AOx activities and the levels of each component in the EOs was undertaken. Table 3 shows the correlation coefficients between the FRAP, ABTS, and DPPH values and the components previously identified as the main compounds in EOs, as far as the sum of the oxygenated monoterpenes and oxygenated sesquiterpenes. According to this correlation analysis (Table 3), several observations can be noted. The compounds contributing most to the variation in antioxidant activities among the major compounds of EOs were: (1) bisabolol, δ-cadinene, and ß-pinene for the FRAP test; (2) β-caryophyllene, bornyl acetate, caryophyllene oxide, α-cadinol, and oxygenated monoterpenes for the DPPH test; (3) caryophyllene oxide and α-cadinol for the ABTS test.
The change in antioxidant activity can be explained by the variation in the chemical composition of the EOs. However, it is nevertheless difficult to attribute these differences to a single component [53]. In general, strong correlations between individual metabolites and the antioxidant activity of plant extracts are rarely observed, even if multivariate models can show the predictive power of the metabolome for antioxidant activity and can indicate which groups of metabolites contribute most to this activity [55,56]. Sometimes, in addition to the major compounds, minor compounds can also play a significant role in the AOx activity of an EO. In fact, the reducing activity of EOs is obviously related to the oxygenated compounds and the terpenes with conjugated double bonds in these EOs [57].

2.5. Profiling of the Fatty Acids from the Leaf Acyl Lipids

Because we suspected a main link between climatic environment (i.e., sites and seasons) and EO composition, we decided to evaluate the fatty acyl (FA) profile of the leaves, with a more dedicated view toward the level of FA unsaturation, which is commonly associated with a physiological adaptation of plant tissues to lower temperatures. This was decided in order to correlate symptoms of low temperature exposure (fatty acid desaturases overexpression) to specific EO profiles, even if we do not suspect a direct metabolic link associating the terpene and fatty acid synthesis (i.e., no presumable fluxomic correlation). It was chosen to extract total lipids from leaves and to trans-esterify all the esterified acyls (from acyl lipids) using a TMAH reactant [58]. The resulting fatty acid methyl esters (FAMEs) were analyzed by GC-MS. The FA composition of the different leaves is presented in Table 4.
Ten FAMEs were identified (Table 4). They represented 94 to 97% (peak area % in TIC) of the extracts. The leaves were all characterized by important levels of polyunsaturated fatty acids (PUFAs): 27–55% linoleic and linolenic FA. The highest PUFA levels were detected in winter in the mountains (55%: linoleic 20% and linolenic 35%). Saturated fatty acids (SFAs), mainly palmitic and stearic FAs, reached 34 to 39%. The monounsaturated oleic acid varied between 3 and 31%. The relative amount of the FA classes were in the order PUFAs > SFAs > MUFAs. The main fatty acids detected in the leaf lipids were linolenic acid (17–35%), palmitic acid (24–31%), linoleic acid (9–28%), and oleic acid (2–30%).
In comparison with the few previous studies, our investigation is the first description of the variation in FA contents in leaves of Algerian lentisk according to the period and altitude at harvest. Our quantitative results are mostly comparable to those of Diamantoglou et al. [59], obtained in Greece, and those of Akdemir et al. [60], in Turkey. Our lentisk leaves from the Algerian Tizi-Ouzou region presented FA proportions comparable to those of the Greece and Turkish ones for C14:0 (myristic), C16:0 (palmitic), C16:1 (palmitoleic), C18:2 (linoleic), and C18:3 (linolenic). On the other hand, behenic (C22:0) was not detected in the samples from Greece and Turkey. Furthermore, the results of Harrat et al. [61] (P. lentiscus from Algeria) reported the presence of other minor fatty acids in the leaves such as C12:0, C13:0, C14:1, C15:0, and C20:0. We did not significantly observe the presence of these minor FAs in our samples. In our quantitative approach, the interaction between altitude and harvesting season was tested (Table 4, column: Altitude * Season Interaction). It had a significant effect on the variation of the levels of the main FAs. This was not true for C16:0, C18:0, and C22:0. This interaction was also significant on the variation of saturated fatty acid (SFA) and polyunsaturated fatty acid contents (p < 0.01). It was less significant for monounsaturated fatty acids (MUSFAs) (p < 0.05). The ANOVA results showed that the seasonal effect had more impact (p < 0.001) on the variation in the FA contents than the altitude effect. In order to develop a better view of our results, we performed principal component analysis (Figure 2).

2.6. Multivariate Analysis of the Fatty Acids Detected (as Methyl Ester: FAMEs) in the Leaf Acyl-Lipids

The two axes of the PCA analysis of the FAMEs show 94.7% of the variance of the data. The score plot (Figure 2A) separates the samples from mountain and littoral, principally, for the most extreme conditions (i.e., mountain in winter and littoral in summer) in the first dimension of the PCA. The second dimension of the PCA mainly separates the mountain samples in the summer. The loading plot (Figure 2B) shows that unsaturated fatty acids (i.e., C18:1, C18:2, and C18:3) were the main markers of seasonal and altitudinal variations within the analyzed samples. In the case of monounsaturated fatty acids, the oleic acid content increased significantly during the summer months and decreased again strongly in winter at both sites. However, the levels of this acid were higher at the lowest altitude (i.e., littoral). Polyunsaturated fatty acids (PUFAs) decreased during the warm season and increased in winter at the two growing areas. Linolenic acid (a triunsaturated PUFA, also qualified as an omega 6 fatty acid) reached its maximum (35%) of total FAs during the winter months at high altitude. In contrast, linoleic acid was characteristic of the winter samples at low altitude. During the winter months, PUFAs generally increased [62]. This can be attributed to a higher expression of delta12 and delta15 desaturases genes that drive the respective transformation of (i) oleic acid to linoleic acid and (ii) linoleic acid to linolenic acid [63]. This corresponds to an acclimatization of the plant to lower temperatures [59,62]. At the level of cell membranes, it allows a maintenance of phospholipid fluidity and, thus, of membrane activity. Such an increase in the level of polyunsaturated fatty acyls in plant cell membranes was associated with a progressive decrease in the daily temperatures at the end of autumn and the beginning of winter. It also includes the effects of decreasing temperature during the night. This phenomenon was amplified at higher altitudes, where night temperatures are usually colder (and earlier) compared to littoral areas.
The correlation coefficients between the main FAs in leaves and the main terpenoids from the EOs (Table 5) showed that β-caryophyllene and bisabolol were oppositely associated with linolenic (C18:3) and oleic (C18:1) acyls. When colder conditions occurred, C18:3 and bisabolol levels increased. When medium temperatures occurred, C18:1 and β-caryophyllene levels were higher. There was no direct link established between the biosynthetic pathway of fatty acids and terpenoids, but this correlation clearly demonstrated that the variation in temperatures from the mountain versus littoral site explained most of the variations in the EO compositions. Finally, two main EO chemotypes could be distinguished: a “cold winter chemotype”, characterized by relatively high levels of bisabolol (Group 3 in Figure 2); and a “temperate chemotype” characterized by high levels of β-caryophyllene. The chemotype of temperate samples can be subdivided into three sub-chemotypes: Groups 1, 4, and 5 as discussed above in Figure 2.

3. Materials and Methods

3.1. Sampling and Botanical Identification

Samples from Pistacia lentiscus were collected from plants in the Algerian region of Tizi-Ouzou. Locally called “Wilaya of Tizi-Ouzou”, this region covers approximately 3000 km2 (100 km EW; 50 km NS) and is limited at the south side by the Djurdjura Mountains and to the north by the Mediterranean Sea. Its relatively contrasted climatic regimes are representative of those of the larger zone of the north of Algeria in between the coast and the mountains. For this study, we selected two geographical sites in the region of Tizi-Ouzou: The first site was located in Ait-Irane, on a mountain in the Tizi-Ouzou region (longitude 36°29′58.3″ N, latitude 4°04′43.4″ E, altitude 876 m). Its Mediterranean climate is characterized by a dry summer and relatively low temperatures during the winter (Table S1; [64], see Supplementary Materials). The second sampling site was selected in Tigzirt, on the littoral of the Tizi-Ouzou region (longitude 36°53′43.0″ N, latitude 4°11′00.4″ E, altitude 13 m). Its summer and winter temperatures are less extreme (Table S2; [65]) and more humidity is present all year long (Table S2; [65]). The contrast between the night/day temperatures is also lower (Table S2; [65]). At each site, a set of three individual female trees were chosen for all the different samplings. The average height of each tree was between 3.1 and 3.7 m. The leaves were collected at human height (1.6–2 m, i.e., halfway up the crown of each shrub) in order to present a homogeneous set of leaves. Each leaf sample was taken from 4 branches located at four opposite sites all around each tree. Each sample was taken in a staggered manner so as not to resample areas that had already been taken in the previous–following season. Sampling was conducted in the morning at sunrise (6–8 a.m.). Samples were taken at both sites on the same plants over four seasons in 2019 at four different vegetative stages: vegetative overwintering (mid-February), full flowering (mid-April), early fruiting (mid-August), and late fruiting (mid-October).
The plant species, Pistacia lentiscus, was identified using the Flora of Quézel and Santa (1962–1963). Furthermore, samples of lentisk were deposited at the ENSA (National High School of Agronomy in Algeria) herbarium. This identification was confirmed by Benhouhou and this led to the delivery of a certificate (ref: 18112020/22:09/ENSA24112020) attesting that this plant was, indeed, Pistacia lentiscus.

3.2. Extraction of Essential Oils

Leaf samples were cleaned of debris and air dried in the dark for 7 days at a temperature of 25–30 °C. The essential oils (EOs) were extracted by a hydrodistillation process using a Clevenger-type apparatus. Then, the EOs were stored in opaque glass bottles at 4 °C until further use. According to the period and altitude of the harvest, eight groups of essential oils were obtained: winter mountain, spring mountain, summer mountain, autumn mountain, winter littoral, spring littoral, summer littoral, and autumn littoral.

3.3. GC-MS Analysis of Essential Oils

The chemical analysis of the EOs was carried out by gas chromatography mass spectrometry (GC-MS) using a Trace GC-Ultra instrument coupled to a DSQ-II (single-quadrupole) mass spectrometer supplied by Thermo Fisher Scientific (Asnières-sur-Seine, France). The chromatographic separation was conducted on a 30 m × 0.25 mm × 0.25 μm TR-5 MS capillary column (5% phenyl and 95% polysiloxane from Thermo-France). One microliter of sample was injected in the splitless mode at 200 °C. The oven temperature was set at 50 °C for three minutes; then, it was programmed to increase from 50 to 270 °C at a rate of 5/min and then to increase to 330 °C at a rate of 10/min. Helium was used as carrier gas (flow rate 1 mL/min). The transfer line was heated to 280 °C. Ionization in the electron impact mode was conducted at 220 °C with an energy of 70 eV. A full-scan detection for masses (m/z) between 35 and 350 UMA was performed at a rate of 3000 UMA/sec (4.18 scan per second). Each sample of essential oil was diluted 1000 times in hexane containing pentadecane (1/4000 v/v) as the internal standard.
The identification of oil components was assigned by comparison of their retention times and retention indices (RIs). Experimental RIs were calculated as the relative retention times of compounds compared with those of (C10–C40) n-alkanes. Comparison of experimental RIs (RIxp) was conducted with RIs from the literature [2,16,18,19,23,29,31,37,44,66] and with those of authentic compounds available in the authors’ laboratory. Moreover, identification was confirmed by comparing the mass spectra of each molecule with (i) those available in the NIST 05 library of the GC/MS data system and with(ii) the mass spectra published [67]. The relative quantification of each compound within a sample was based on the peak relative (%) surface measured with the total ion current (TIC) after normalization of the surface area to the internal standard.

3.4. Antioxidant Activity

The antioxidant activity of EOs was assessed using three techniques: the DPPH radical-scavenging test, the ferric reducing antioxidant power (FRAP) assay, and the ABTS radical-cation reduction test. All measurements were spectrophotometry determined (see below) and expressed in Trolox equivalent antioxidant capacity (TEAC). Trolox is the hydrophilic equivalent of vitamin E and it was used as the standard for all three tests. All measurements were performed in triplicate and the results are expressed as Trolox equivalent (i.e., mg of Trolox equivalent/g of essential oil; mg TE/g EO) for all tests (i.e., DPPH, FRAP, and ABTS). Vitamin C (ascorbic acid) was tested as the control for the three tests; the results are expressed as mg Trolox equivalent/g of vitamin C.

3.4.1. Radical Scavenging Test (DPPH)

An essential oil diluted in absolute ethanol (50 µL; 2 mg/mL) was added to 950 µL of freshly prepared ethanolic (EtOH 100%) solution of DPPH (100 mM) and diluted to 10 mM on the test day. The mixture was shaken and kept in the dark at room temperature for 20 min. The absorbance was measured at 517 nm. A blank solution was prepared in parallel: 50 µL of absolute ethanol and 950 µL of DPPH 10 mM. The Trolox standard range (i.e., 0, 5, 12.5, 20, and 25 µM) was prepared in absolute ethanol [19,68].

3.4.2. Ferric Reducing Power Test (FRAP)

The FRAP test was prepared by mixing a 10 mM solution of tripyridyltriazine (TPTZ) in hydrochloric acid (HCl 40 mM), a 20 mM ferric chloride (FeCl3) solution, and an acetate buffer (300 mM, pH: 3.6), the proportions of 1:1:10, respectively. The reaction was set up by mixing 900 µL of the above FRAP solution with 50 µL of the ethanol (100%)-diluted essential oil sample (2 mg/mL). The reaction mixture was incubated at 37 °C for 30 min. The absorbance was measured at 593 nm. A blank consisted of 50 µL of absolute ethanol and 950 µL of FRAP solution. The Trolox standard range (i.e., 0, 50, 100, 200 m, and 400 µM) was prepared in absolute ethanol [19,69].

3.4.3. Radical Cation Reduction Test (ABTS)

One hundred microliters of essential oil diluted in absolute ethanol (2 mg/mL) was added to 900 µL of the (100%) ethanolic solution of ABTS. ABTS+ was produced by the reaction of 7 mM of ABTS with 2.5 mM of sodium persulfate (K2S2O8) in distilled water for 16 h in the dark. The concentration of ABTS+ was then adjusted by dilution to an absorbance of 0.7 ± 0.02 (at 734 nm). A blank consisted of 100 µL of absolute ethanol and 900 µL of ABTS solution. The Trolox standard range (i.e., 0, 4, 10, 16, and 20 µM) was prepared in absolute ethanol. After incubation in the dark for 6 min, the absorbance of each sample was measured at 734 nm [70].

3.5. Preparation of Lipidic Extracts from Lentisk Leaves and GC-MS Analysis

The lipid extracts from leaves were prepared for four harvest periods and from both sites. The extraction of leaf lipids was performed mainly according to Akdemir et al. [60]. Powdered air-dried leaves (300 mg) were extracted in 1.5 mL of methanol–chloroform (1–2 V/V) for 30 min under stirring (tripentadecanoin was also added at 1 mg/mL as an internal standard (Sigma-Aldrich, Saint-Quentin-Fallavier, France) for further GC-MS analysis). After a 5 min centrifugation at 1100× g, the supernatant was treated by adding an equivalent volume of hexane. After stirring the mixture, 8 µL of tetramethylammonium hydroxide was added (TMAH; Sigma-Aldrich). The use of tetramethylammonium hydroxide (TMAH) allowed for the preconversion of fatty acids into fatty acid methyl esters (FAMEs) in order to increase their volatility for GC-MS analysis [58]. After resting, the supernatant was recovered and put in vials for GC-MS analysis. For each sample, the separation of the different compounds was performed on a TR-5 MS (5% phenyl and 95% polysiloxane) capillary GC column (30 m × 0.25 mm × 0.25 μm) from Thermo (Asnières-sur-Seine, France). One microliter of each GC-MS sample was injected in the splitless mode at 230 °C. The detector temperature was 250 °C. A gradient of oven temperatures from 230 °C to 290 °C was used (230 to 245 °C at 15 °C·min−1; 245 to 290 °C at 90 °C·min−1; 290 °C for 2 min). Helium was used as the carrier gas at a flow rate of 1 mL/min. The ionization in the electron impact mode was conducted with an energy of 70 eV and full-scan detection for masses (m/z) between 35 and 350 atomic mass unit (AMU). The retention indices were calculated according to a series of n-alkanes (C10–C40). The identification was performed by comparison of their mass spectra with those stored in the NIST 05 library (protocol mainly issuing from [71]).

3.6. Statistical Analysis

All determinations were performed in triplicate. The results were then expressed as the mean ± standard deviation of three measurements. Statistical analyses were performed in R (R core team, 2021) [72]. The data were analyzed by comparison of the means using a two-factor analysis of variance ANOVA model of Type III, followed by a Tukey test for the determination of significant differences in concentrations of all the detected compounds in the different samples: EOs (with monoterpenes, sesquiterpenes, and others), trans-esterified lipid extracts (fatty acid methyl esters). For each analyzed compound, the assumptions of the ANOVA model (normal distribution and homoscedasticity of the error terms) were checked using the diagnostic plots generated by the plot() function in R, and, if necessary, square root transformation was applied.
The same procedure was used to assess the variation in the antioxidant activities among samples. The correlations between components of the essential oils and antioxidant activities were estimated using Spearman’s correlation coefficient (rs). The same correlation analysis was conducted for major components of the essential oils and acyl composition. Differentiation between populations (based on their EO and acyl compositions) was assessed by principal component analysis (PCA) using R library FactoMiner R (R Core Team 2020) [72,73].

4. Conclusions

This paper showed a large variability in the terpenoid (i.e., essential oils) and fatty acid compositions of leaves from Pistacia lentiscus collected at two contrasting natural sites in the Algerian region of Tizi-Ouzou. Wild plants naturally growing at high altitudes (mountain site at 876 m) and at low altitudes (littoral site at 13 m) were selected, and some of their leaves were consecutively sampled four times during their different vegetative–reproductive stages over one year. The extraction of essential oils (EOs) from the leaves and the further gas chromatography analysis allowed for the identification and relative quantification of 47 compounds. Based on the main constituents of the EOs, different chemotypes could be described. EOs from the coastal site corresponded to four harvest periods: winter, spring, summer, and autumn. They were characterized by high amounts of β-caryophyllene (C26). EOs from the mountain site were seasonally variable. In winter, EOs contained high amounts of δ-cadinene (C36), bisabolol (C46), and bornyl acetate (C20). Samples collected in spring and summer were characterized by high amounts of α-pinene (C2) and β-pinene (C5). Samples collected in autumn had a high proportion of β-myrcene (C6) followed by α-pinene (C2). The chemical class of monoterpene hydrocarbons was higher in the mountains. In addition, the highest rate of this class was observed in autumn for both sites. Contrary to monoterpene hydrocarbons, sesquiterpene hydrocarbons were present at higher concentrations in the EOs issuing from the coastal area. The antioxidant (AOx) activity of the different EOs varied according to the harvesting period and the altitude. The highest antiradical effects were observed with the FRAP test. The variation in the AOx activities can partially be correlated to the composition of the essential oils. This study thus showed that the change in altitude and harvesting period affected the chemotypic nature of the lentisk EOs as well as their AOx potential, especially at the highest altitudes. The Pistacia lentiscus natural population of the Tizi-Ouzou littoral can thus be used for the production of a quite constant chemotyped EO. Oppositely, because of the seasonal variability, EOs from high-altitude populations may be better suited for obtaining dedicated original chemotypes and for the further cultivation and exploitation of their (potentially therapeutic) potential.
The fatty acyl (FA) profile of the leaf acyl lipids was analyzed by GC-MS after total lipid extraction and transesterification into FA methyl esters (FAMEs). The resulting composition was characterized by its richness in polyunsaturated fatty acids (PUFAs) especially in winter in the mountain. PUFA levels were followed by saturated fatty acids (SFAs) and monounsaturated fatty acids (MUFAs). The amounts of FAs decreased in the order PUFA > SFA > MUFA. The major FAs detected in the leaves were linolenic acid, palmitic acid, linoleic acid, and oleic acid. Unsaturated fatty acids (i.e., C18:1, C18:2, and C18:3) levels significantly varied with season and altitude. The colder the climatic conditions the higher the relative ratio of PUFAs. The variability in the EOs’ chemotypes can be associated with this physiological response of plants to temperatures. The leaves harvested at the coldest season were characterized by the highest levels of bisabolol. The EOs from other sites and seasons were characterized by three other chemotypes in which β-caryophyllene was always the main compound. From these results, a planification of leaf harvest can be conducted for producing EOs with specific chemotypes. It is the first step in the further sustainable and rational valorization of natural resources from the region of Tizi-Ouzou. It is also an opportunity to more deeply study the impact of climatic changes on the resulting terpenoid biosynthesis using correlation networks and fluxomic analysis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27134148/s1: Table S1: Weather data collected between 1999 and 2019 in Ait Irane, Algeria; Table S2: Weather data collected between 1999 and 2019 in Tigzirt, Algeria.

Author Contributions

The individual contribution of authors was as follows. C.S. was in charge of all phases (i.e., design, experimentation, extractions, analyses, statistical evaluations, interpretation of results, and writing of the article) of the work that is included in her PhD based on an Algeria–France collaboration. N.J. was in charge of part of the supervision of that work and contributed to the writing and revision of the manuscript. E.C. supervised the experiments devoted to antioxidant analysis and contributed to the specific interpretation and the revision of these aspects in the manuscript. R.D. supervised the statistical analysis and contributed to the interpretation and the revision of the manuscript. J.X.F. supervised the statistical analysis and contributed to the interpretation and the revision of the manuscript. R.M. supervised the statistical analysis and contributed to the interpretation and the revision of the manuscript. F.A. contributed to the initiation of the project with her process engineering expertise. F.F. contributed to the whole project through supervision of the dedicated sample collection, essential oil production, and participation in the interpretation of the results, the writing, and the improvement of this article. E.G. contributed to the whole project through supervision of and participation to the analytical work, the interpretation of results, the writing, and the improvement of this article. Conceptualization, C.S., F.F. and E.G.; methodology, C.S., N.J., E.C., R.D., J.X.F., R.M., F.A., F.F. and E.G.; software, C.S., R.D., J.X.F. and R.M.; validation, C.S., R.D., J.X.F., R.M., F.F. and E.G.; formal analysis, C.S., F.F. and E.G.; investigation, C.S., F.F. and E.G.; resources, F.F. and E.G.; data curation, C.S. and E.G.; writing—original draft preparation, C.S., F.F. and E.G.; writing—review and editing, C.S., F.F. and E.G.; visualization, C.S., N.J., E.C., R.D., J.X.F., R.M., F.A., F.F. and E.G.; supervision, N.J., F.A., F.F. and E.G.; project administration, F.F. and E.G.; funding acquisition, F.F. and E.G. Authorship here is limited to those who have contributed substantially to the work reported. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by (1) the French Ministry of Foreign Affairs through Campus France (file no. 973496A); (2) the “Conseil Régional de Picardie” through the MTV (Métabo-Typage Végétal) project; (3) the University of Picardie Jules Verne through the French Ministry of Higher Education and Research, region “Hauts de France”, the European FEDER funds within the so-called CPER programs devoted to equipment acquisitions. This work was also funded by Algerian authorities though the financial support of the Tizi-Ouzou UMMTO University and a social grant attributed to Chabha SEHAKI for her PhD.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Acknowledgments

The authors thank Salima Benhouhou (Ecole Nationale Supérieure Agronomique, ENSA, Algiers, Algeria) for the identification of the Pistacia lentiscus plants, the recording of samples, and the delivery of a certificate.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

References

  1. Benhammou, N.B.; Belyagoubi, L.; El Zerey-Belaskri, A.; Zitouni, A.; Ghembaza, N.; Benhassaini, H.; Atik-Bekkara, F.; Piras, A.; Falconieri, D.; Rosa, A. Fatty acid composition and antioxidant activity of Pistacia lentiscus L. fruit fatty oil from Algeria. J. Food Meas. Charact. 2018, 12, 1408–1412. [Google Scholar] [CrossRef]
  2. Zaouali, Y.; Bel HadjYahya, I.; Jaouadi, R.; Messaoud, C.; Boussaid, M. Sex-related differences in essential oil composition, phenol contents and antioxidant activity of aerial parts in Pistacia lentiscus L. during seasons. Ind. Crop. Prod. 2018, 121, 151–159. [Google Scholar]
  3. Gardeli, C.; Papageorgiou, V.; Mallouchos, A.; Theodosis, K.; Komaitis, M. Essential oilcomposition of Pistacia lentiscus L. and Myrtuscommunis L.: Evaluation of antioxidant capacity of methanolic extracts. Food Chem. 2008, 107, 1120–1130. [Google Scholar] [CrossRef]
  4. Landau, S.; Muklada, H.; Markovics, A.; Azaizeh, H. Traditional uses of Pistacia lentiscus in veterinary and human medicine. In Medicinal and Aromatic Plants of the Middle-East; Springer: Dordrecht, The Netherlands, 2014; Volume 2, pp. 163–180. [Google Scholar]
  5. Mehenni, C.; Atmani-Kilani, D.; Dumarçay, S.; Perrin, D.; Gérardin, P.; Atmani, D. Hepatoprotective and antidiabetic effects of Pistacia lentiscus leaf and fruit extracts. J. Food Drug Anal. 2016, 24, 653–669. [Google Scholar] [CrossRef] [PubMed]
  6. Said, O.; Khalil, K.; Fulder, S.; Azaizeh, H. Ethno-pharmacological survey of medicinal herbs in Israel, the Golan Heights and the West Bank region. J. Ethnopharmacol. 2002, 83, 251–265. [Google Scholar] [CrossRef]
  7. Jamila, F.; Mostafa, E. Ethnobotanical survey of medicinal plants used by people in Oriental Morocco to manage various ailments. J. Ethnopharmacol. 2014, 154, 76–87. [Google Scholar] [CrossRef]
  8. Guarrera, P.M.; Salerno, G.; Caneva, G. Folk phytotherapeutical plants from Maratea area (Basilicata, Italy). J. Ethnopharmacol. 2005, 99, 367–378. [Google Scholar] [CrossRef]
  9. Abdeldjelil, M.C.; Bensegueni, A.; Messaï, A.; Agabou, A.; Benazzouz, H. Medicinal use of Pistacia lentiscus fixed oil in Constantine province, north-east Algeria. J. Nat. Prod. Plant Resour. 2014, 4, 48–51. [Google Scholar]
  10. Nahida, A.S.; Siddiqui, A.N. Pistacia lentiscus: A review on phytochemistry and pharmacological properties. Int. J. Pharm. Pharm. Sci. 2012, 4, 16–20. [Google Scholar]
  11. Bozorgi, M.; Memariani, Z.; Mobli, M.; Surmaghi, S.M.H.; Reza Shams-Ardekan, M.; Rahimi, R. Five pistacia species (P. vera, P. atlantica, P. terebinthus, P. khinjuk and P. lentiscus): A review of their traditional uses, phytochemistry, and pharmacology. Sci. World J. 2013, 2013, 219815. [Google Scholar] [CrossRef] [Green Version]
  12. Gratani, L.; Catoni, R.; Varone, L. Morphological, anatomical and physiological leaf traits of Q. ilex, P. latifolia, P. Lentiscus. And M. communis and their response to Mediterranean climate stress factors. Bot. Stud. 2013, 54, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gourine, N.; Yousfi, M.; Bombarda, I.; Nadjemi, B.; Stocker, P.; Gaydou, E.M. Antioxidant activities and chemical composition of essential oil of Pistacia atlantica from Algeria. Ind. Crop. Prod. 2010, 31, 203–208. [Google Scholar] [CrossRef]
  14. Ahmed, Z.B.; Yousfi, M.; Viaene, J.; Dejaegher, B.; Demeyer, K.; Vander Heyden, Y. Four Pistacia atlantica subspecies (atlantica, cabulica, kurdica and mutica): A review of their botany, ethnobotany, phytochemistry and pharmacology. J. Ethnopharmacol. 2021, 265, 113329. [Google Scholar] [CrossRef] [PubMed]
  15. Ben Youssef, E.H.; Charchari, S.; Nacer-Bey, N. The essential oil of Pistacia lentiscus L. from Algeria. J. Essent. Oil Res. 2005, 17, 642–644. [Google Scholar] [CrossRef]
  16. Mecherara-Idjeri, S.; Hassani, A.; Castola, V.; Casanova, J. Composition and chemical variability of the essential oil from Pistacia lentiscus L. growing wild in Algeria Part 1: Leaf oil. J. Essent. Oil Res. 2008, 20, 32–38. [Google Scholar] [CrossRef]
  17. Amhamdi, H.; Aouinti, F.; Wathelet, J.P.; Elbachiri, A. Chemical composition of the essential oil of Pistacia lentiscus L. from eastern Morocco. Rec. Nat. Prod. 2009, 3, 90–95. [Google Scholar]
  18. Hafsé, M.; Benbrahim, K.F.; Saidi, A.; Farah, A. Volatile components and antibacterial profile of essential oils extracted from leaves and twigs of Pistacia lentiscus L. Br. Microbiol. Res. J. 2013, 3, 602–611. [Google Scholar] [CrossRef]
  19. Aissi, O.; Boussaid, M.; Messaoud, C. Essential oil composition in natural populations of Pistacia lentiscus L. from Tunisia: Effect of ecological factor and incidence on antioxidant and antiacetylcholinesterase activities. Ind. Crop. Prod. 2016, 91, 379–382. [Google Scholar] [CrossRef]
  20. Lo Presti, M.; Sciarrone, D.; Crupi, M.L.; Costa, R.; Ragusa, S.; Dugo, G.; Mondello, L. Evaluation of the volatile and chiral composition in Pistacia lentiscus L. essential oil. Flavour Fragr. J. 2008, 23, 249–257. [Google Scholar] [CrossRef]
  21. Djenane, D.; Yanguela, J.; Montanes, L.; Djerbal, M.; Roncales, P. Antimicrobial activity of Pistacia lentiscus and Satureja montana essential oils against Listeria monocytogenes CECT 935 using laboratory media: Efficacy and synergistic potential in minced beef. Food Control 2011, 22, 1046–1053. [Google Scholar] [CrossRef]
  22. Ismail, A.; Hamrouni, L.; Hanana, M.; Jamoussi, B. Chemical composition and herbicidal effects of Pistacia lentiscus L. essential oil against weeds. Int. J. Med. Aromat. Plants 2012, 2, 558–565. [Google Scholar]
  23. Haouli, A.; Seridi, R.; Djemli, S.; Bourdjiba, O.; Frih, H. Contribution to the analysis of Pistacia lentiscus extracted oil. Am.-Eur. J. Agric. Environ. Sci. 2015, 15, 1075–1081. [Google Scholar]
  24. Douissa, F.B.; Hayder, N.; Chekir-Ghedira, L.; Hammami, M.; Ghedira, K.; Mariotte, A.M.; Dijoux-Franca, M.G. New study of the essential oil from leaves of Pistacia lentiscus L. (Anacardiaceae) from Tunisia. Flavour Fragr. J. 2005, 20, 410–414. [Google Scholar] [CrossRef]
  25. Nahum, S.; Inbar, M.; Ne’eman, G.; Ben-Shlomo, R. Phenotypic plasticity and gene diversity in Pistacia lentiscus L. along environmental gradients in Israel. Tree Genet. Genomes 2008, 4, 777–785. [Google Scholar] [CrossRef]
  26. Ait Said, S.; Fernandez, C.; Greff, S.; Torre, F.; Derridj, A.; Gauquelin, T.; Mevy, J.P. Inter-population variability of terpenoid composition in leaves of Pistacia lentiscus L. from Algeria: A chemoecological approach. Molecules 2011, 16, 2646–2657. [Google Scholar] [CrossRef]
  27. Carvalho, S.; Macel, M.; Mulder, P.P.J.; Skidmore, A.; van der Putten, W.H. Chemical variation in Jacobaea vulgaris is influenced by the interaction of season and vegetation successional stage. Phytochemistry 2014, 99, 86–94. [Google Scholar] [CrossRef]
  28. Juvany, M.; Müller, M.; Pintó-Marijuan, M.; Munné-Bosch, S. Sex-related differences in lipid peroxidation and photoprotection in Pistacia lentiscus. J. Exp. Bot. 2014, 65, 1039–1049. [Google Scholar] [CrossRef] [Green Version]
  29. Negro, C.; De Bellis, L.; Miceli, A. Chemical composition and antioxidant activity of Pistacia lentiscus essential oil from Southern Italy (Apulia). J. Essent. Oil Res. 2015, 27, 23–29. [Google Scholar] [CrossRef]
  30. Foddai, M.; Kasabri, V.; Afif, F.U.; Azara, E.; Petretto, G.L.; Pintore, G. In vitro inhibitory effects of Sardinian Pistacia lentiscus L. and Pistacia terebinthus L. on metabolic enzymes Pancreatic lipase, -amylase, and -glucosidase. Starch-Stärke 2014, 67, 204–212. [Google Scholar] [CrossRef]
  31. Llorens-Molina, J.A.; Gonzalez, S.V.; Martinez, J.S. Essential oil composition of leaves of Pistacia lentiscus L. growing wild in Valencia (Spain). Nat. Volatiles Essent. Oils 2015, 2, 17–26. [Google Scholar]
  32. Bachrouch, O.; Jemâa, J.M.-B.; AidiWaness, W.; Talou, T.; Marzouk, B.; Abdrerraba, M. Composition and insecticidial activity of essential oil from Pistacia lentiscus L. against Ectomyelois ceratoniae Zeller and Ephestia kuehniella Zeller (Lepidoptera: Pyralidae). J. Stored Prod. Res. 2010, 46, 242–247. [Google Scholar] [CrossRef]
  33. Klibet, F.; Boumendjel, A.; Khiari, M.; El Feki, A.; Abdenour, C.; Messarah, M. Oxidative stress-related liver dysfunction by sodium arsente: Alleviatio by Pistacia lentiscus oil. Pharm. Biol. 2016, 54, 354–363. [Google Scholar] [CrossRef] [PubMed]
  34. Maes, C.; Meersmans, J.; Lins, L.; Bouquillon, S.; Fauconnier, M.L. Essential Oil-Based Bioherbicides: Human Health Risks Analysis. Int. J. Mol. Sci. 2021, 22, 9396. [Google Scholar] [CrossRef]
  35. Damiens, A.; Alebrahim, M.T.; Léonard, E.; Fayeulle, A.; Furman, C.; Hilbert, J.L.; Siah, A.; Billamboz, M. Sesamol–based terpenoids as promising bio-sourced crop protection compounds against the wheat pathogen Zymoseptoriatritici. Pest Manag. Sci. 2021, 77, 2403–2414. [Google Scholar] [CrossRef] [PubMed]
  36. Vidrich, V.; Fusi, P.; Graziano, A.; Silvestrini, E.; Michelozzi, M.; Marco, F. Chemical composition of the essential oil of Pistacia Lentiscus L. J. Essent. Oil Res. 2004, 16, 223–226. [Google Scholar] [CrossRef]
  37. El Bishbishy, M.H.; Gad, H.A.; Aborehab, N.M. Chemometric discrimination of three Pistacia species via their metabolic profiling and their possible in vitro effects on memory functions. J. Pharm. Biomed. Anal. 2020, 177, 112840. [Google Scholar] [CrossRef]
  38. Congiu, R.; Falconieri, D.; Marongiu, B.; Piras, A.; Porcedda, S. Extraction and isolation of Pistacia lentiscus L. essential oil by supercritical CO2. Flavour Fragr. J. 2002, 17, 239–244. [Google Scholar] [CrossRef]
  39. Bouyahya, A.; Assemian, I.C.C.; Mouzount, H.; Bourais, I.; Et-Touys, A.; Fellah, H.; Benjouad, A.; Dakka, N.; Bakri, Y. Could volatile compounds from leaves and fruits of Pistacia lentiscus constitute a novel source of anticancer, antioxidant, antiparasitic and antibacterial drugs? Ind. Crop. Prod. 2019, 128, 62–69. [Google Scholar] [CrossRef]
  40. Derwich, E.; Manar, A.; Benziane, Z.; Boukir, A. GC/MS analysis and In vitro antibacterial activity of the essential oil isolated from leaf of Pistacia lentiscus growing in morocoo. World Appl. Sci. J. 2010, 8, 1267–1276. [Google Scholar]
  41. Aouinti, F.; Zidane, H.; Tahri, M.; Wathelet, J.P.; El Bachiri, A. Chemical composition, mineral contents and antioxidant activity of fruits of Pistacia lentiscus L. from Eastern Morocco. J. Mater. Environ. Sci. 2014, 5, 199–206. [Google Scholar]
  42. Kıvçak, B.; Akay, S. Quantitative determination of α-tocopherol in Pistacia lentiscus, Pistacia lentiscus var. chia, and Pistacia terebinthus by TLC-densitometry and colorimetry. Fitoterapia 2005, 76, 62–66. [Google Scholar] [CrossRef]
  43. El Idrissi, M.; Barbouchi, M.; Choukrad, M.; Lhoussain, L. Chemical composition and antimicrobial activity of essential oils isolated from leaves and twigs of Pistacia lentiscus L. Growing wild in Morocco. World J. Pharm. Sci. 2016, 5, 516–524. [Google Scholar]
  44. Hamiani, A.; Kambouche, N.; Bellahouel, S. Chemical composition and antioxidant potential of Pistacia lentiscus L. essential oil from Oran (Algeria). Adv. Biosci. Biotechnol. 2016, 7, 539–544. [Google Scholar]
  45. Medini, H.; Elaissi, A.; Khouja, M.L.; Chraief, I.; Farhat, F.; Hammami, M.; Chemli, R.; Harzallah-Skhiri, F. Leaf essential oil of Juniperusoxycedrus L. (Cupressaceae) harvested in northern Tunisia: Composition and intra-Specific variability. Chem. Biodivers. 2010, 7, 1254–1266. [Google Scholar] [CrossRef]
  46. Lakušić, D.V.; Ristić, M.S.; Slavkovska, V.N.; Šinžar-Sekulić, J.B.; Lakušić, B.S. Environment-related variations of the composition of the essential oils of rosemary (Rosmarinus officinalis L.) in the balkan penninsula. Chem. Biodivers. 2012, 9, 1286–1302. [Google Scholar] [CrossRef]
  47. Rahimmalek, M.; Mirzakhani, M.; Pirbalouti, A.G. Essential oil variation among 21 wild myrtle (Myrtus communis L.) populations collected from different geographical regions in Iran. Ind. Crop. Prod. 2013, 51, 328–333. [Google Scholar] [CrossRef]
  48. Djerrad, Z.; Kadik, L.; Djouahri, A. Chemical variability and antioxidant activities among Pinushalepensis Mill. essential oils provenances, depending on geographic variation and environmental conditions. Ind. Crop. Prod. 2015, 74, 440–449. [Google Scholar] [CrossRef]
  49. Edreva, A.; Velikova, V.; Tsonev, T.; Dagnon, S.; Gurel, A.; Aktaş, L.; Gesheva, E. Stress-protective role of secondary metabolites: Diversity of functions and mechanisms. Gen. Appl. Plant Physiol. 2008, 34, 67–78. [Google Scholar]
  50. Kofidis, G.; Bosabalidis, A.M. Effects of altitude and season on glandular hairs and leaf structural traits of Nepetanuda L. Bot. Stud. 2008, 49, 363–372. [Google Scholar]
  51. Abreu, I.N.; Mazzafera, P. Effect of water and temperature stress on the content of active constituents of Hypericum brasiliense Choisy. Plant Physiol. Biochem. 2005, 43, 241–248. [Google Scholar] [CrossRef] [PubMed]
  52. Cirak, C.; Radusiene, J.; Jakstas, V.; Ivanauskas, L.; Seyis, F.; Yayla, F. Altitudinal changes in secondary metabolite contents of Hypericum androsaemum and Hypericum polyphyllum. Biochem. Syst. Ecol. 2017, 70, 108–115. [Google Scholar] [CrossRef]
  53. Barra, A.; Coroneo, V.; Dessi, S.; Cabras, P.; Angioni, A. Characterization of the volatile constituents in the essential oil of Pistacia lentiscus L. from different origins and its antifungal and antioxidant activity. J. Agric. Food Chem. 2007, 55, 7093–7098. [Google Scholar] [CrossRef] [PubMed]
  54. Bampouli, A.; Kyriakopoulou, K.; Papaefstathiou, G.; Louli, V.; Krokida, M.; Magoulas, K. Comparison of different extraction methods of Pistacia lentiscus var. chia leaves: Yield: Antioxidant activity and essential oil chemical composition. J. Appl. Res. Med. Aromat. Plants 2014, 1, 81–91. [Google Scholar] [CrossRef]
  55. Desmet, S.; Morreel, K.; Dauwe, R. Origin and Function of Structural Diversity in the Plant Specialized Metabolome. Plants 2021, 10, 2393. [Google Scholar] [CrossRef]
  56. Nguyen, T.K.O.; Nguyen, N.L.; Pham, H.N.; Sawada, Y.; Hirai, M.Y.; Dauwe, R.; Dijoux-Franca, M.G. Development of a Pteris vittata L. compound database by widely targeted metabolomics profiling. Biomed. Chromatogr. BMC 2021, 35, e5110. [Google Scholar] [CrossRef] [PubMed]
  57. Olszowy, M.; Dawidowicz, A.L. Essential oils as antioxidants: Their evaluation by DPPH, ABTS, FRAP, CUPRAC, and β-carotene bleaching methods. Monatsh. Chem. 2016, 147, 2083–2091. [Google Scholar] [CrossRef]
  58. Mannion, D.T.; Furey, A.; Kilcawley, K.N. Free fatty acids quantification in dairy products. Int. J. Dairy Technol. 2016, 69, 1–12. [Google Scholar] [CrossRef] [Green Version]
  59. Diamantoglou, S.; Meletiou-christou, M.S. Das jahres periodische Verhalten des Rohfetts und der Fettsäuren in Rinden und Bläuern von Pistacia lentiscus, Pistacia terebinthus und Pistacia vera: The Lipid Content and Fatty Acid Composition of Barks and Leaves of Pistacia lentiscus, Pistacia terebinthus and Pistacia vera During the Course of a Year. Z. Pflanzenphysiol. 1979, 93, 219–228. [Google Scholar]
  60. Akdemir, O.F.; Tilka, E.; Onay, A.; Keskin, C.; Bashan, M.; Kilinc, F.M.; Kizmaz, V.; Suzerer, V. Determination of the fatty acid composition of the fruits and different organs of lentisk (Pistacia lentiscus L.). J. Essent. Oil-Bear Plants 2015, 18, 1224–1233. [Google Scholar] [CrossRef]
  61. Harrat, M.; Benalia, M.; Gourine, N.; Yousfi, M. Variability of the chemical compositions of fatty acids, tocopherols and lipids antioxidant activities, obtained from the leaves of Pistacia lentiscus L. growing in Algeria. Mediterr. J. Nutr. Metab. 2018, 11, 199–215. [Google Scholar] [CrossRef]
  62. Martz, F.; KiviniemI, S.; Heino, P.; Palva, T.; Sutinen, M.L. Genetic control of fatty acid desaturation during cold acclimation and de-acclimation of birch. Acta Hortic. 2003, 618, 137–146. [Google Scholar] [CrossRef]
  63. Gontier, E.; Gougeon, S.; Guillot, X.; Thomasset, B.; Méjean, L.; Tran, T.L.M.; Bourgaud, F. Les plantes, sources d’acides gras essentiels oméga 3. Oléagineux Corp. Gras Lipides OCL-J. 2004, 11, 106–111. [Google Scholar] [CrossRef] [Green Version]
  64. Climate-data.org. Available online: https://fr.climate-data.org/afrique/algerie/tizi-ouzou/ait-irane-693594/ (accessed on 5 February 2022).
  65. Climate-data.org. Available online: https://fr.climate-data.org/afrique/algerie/tizi-ouzou/tigzirt-321206/ (accessed on 5 February 2022).
  66. Arabi, A.; Djibaoui, R.; Malihac, C.; Sisbane, I.; Lattab, A.; Bechelaghem, N.; Dahah, H.; Charef, R.; Ettalhi, M.; Taleb, F.; et al. Chemical composition and antibacterial activity of essential oil from leaves and twigs of Pistacia lentiscus growing in Mostaganem Province (Algeria). Int. J. Biosci. 2017, 10, 146–158. [Google Scholar]
  67. Adams, R.P. Identification of Essential Oil Components by Gas Chromatography/Mass Spectrometry, 4th ed.; Allured Publishing Corporation: Carol Stream, IL, USA, 2007; p. 811. [Google Scholar]
  68. Doshi, P.; Adsule, P.; Banerjee, K.; Oulkar, D. Phenolic compounds, antioxidant activity and insulinotropic effect of extracts prepared from grape (Vitis vinifera L.) by products. J. Food Sci. Technol. 2015, 52, 181–190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Messaoud, C.; Laabidi, A.; Boussaid, M. Myrtus communis L. infusions: The effect of infusion time on phytochemical composition antioxidant, and antimicrobial activities. J. Food Sci. 2012, 77, 941–947. [Google Scholar] [CrossRef]
  70. Choque, E.; Durrieu, V.; Alric, I.; Raynal, J.; Mathieu, F. Impact of spray-drying on biological properties of chitosan matrices supplemented with antioxidant fungal extracts for wine applications. Curr. Microbiol. 2020, 77, 210–219. [Google Scholar] [CrossRef] [PubMed]
  71. Jousse, C.; Schiltz, S.; Fournez, A.; Guillot, X.; Thomasset, B.; Gougeon, S.; Bourgaud, F.; Gontier, E. Rapid, cost-effective screening of flax genotypes to identify desirable fatty acid compositions. Electron. J. Plant Breed. 2010, 1, 1396–1404. [Google Scholar]
  72. 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 21 March 2022).
  73. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R package for multivariate analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef] [Green Version]
Figure 1. PCA analysis performed on all the compounds from the essential oils of Pistacia lentiscus leaf samples collected at the two altitude sites (△ for littoral versus ◯ for mountain) and in 4 consecutive seasons: (A) score plot of PC1 versus PC2 scores; (B) loading plot of PC1- and PC2-contributing EO compounds (the most discriminating compounds are represented in color).
Figure 1. PCA analysis performed on all the compounds from the essential oils of Pistacia lentiscus leaf samples collected at the two altitude sites (△ for littoral versus ◯ for mountain) and in 4 consecutive seasons: (A) score plot of PC1 versus PC2 scores; (B) loading plot of PC1- and PC2-contributing EO compounds (the most discriminating compounds are represented in color).
Molecules 27 04148 g001
Figure 2. PCA analysis of the FAMEs from the leaf acyl lipids of Pistacia lentiscus leaf samples collected at the two altitude sites (△ for littoral versus ◯ for mountain) and 4 consecutive seasons: (A) score plot of the PC1 versus PC2 scores; (B) loading plot of PC1- and PC2-contributing FAME compounds (the most discriminating compounds are represented in color).
Figure 2. PCA analysis of the FAMEs from the leaf acyl lipids of Pistacia lentiscus leaf samples collected at the two altitude sites (△ for littoral versus ◯ for mountain) and 4 consecutive seasons: (A) score plot of the PC1 versus PC2 scores; (B) loading plot of PC1- and PC2-contributing FAME compounds (the most discriminating compounds are represented in color).
Molecules 27 04148 g002
Table 1. Compounds detected in the essential oils of Pistacia lentiscus leaves collected at two altitude sites and in four consecutive seasons. The levels of each compound are expressed as the chromatographic relative peak surface measured as the total ion current (area % TIC). The extraction yield for each EO was similar: 0.125 ± 0.005 g/100 g of dry matter.
Table 1. Compounds detected in the essential oils of Pistacia lentiscus leaves collected at two altitude sites and in four consecutive seasons. The levels of each compound are expressed as the chromatographic relative peak surface measured as the total ion current (area % TIC). The extraction yield for each EO was similar: 0.125 ± 0.005 g/100 g of dry matter.
Means % Content ± Standard Deviation
MountainLittoralAltitude * Season Interaction
C CompoundsIRWinterSpringSummerAutumnWinterSpringSummerAutumn
C1MHtricyclene9230.2 ± 0.1 ab0.9 ± 0.1 b0.8 ± 0.1 ab0.5 ± 0.1 a0.1 ± 0.1 c0.3 ± 0.1 d0.2 ± 0.1 d0.2 ± 0.1 ab***
C2MHα-pinene9352.6 ± 0.2c13.1 ± 0.1b12.6 ± 0.3b9.4 ± 0.9a2.1 ± 0.4c5.4 ± 0.4d4.9 ± 0.5d8.3 ± 0.3a***
C3MHcamphene9520.9 ± 0.1 cd3.1 ± 0.1 b3.2 ± 0.0 d1.6 ± 0.0 a0.6 ± 0.1 e0.8 ± 0.1 f1.1 ± 0.1 f1.0 ± 0.0 bc***
C4MHsabinene9750.5 ± 0.0 e1.4 ± 0.0 c1.2 ± 0.1 c3.0 ± 0.4 a0.3 ± 0.1 f1.1 ± 0.1 d1.0 ± 0.1 cd2.0 ± 0.1 b*
C5MHβ-pinene9810.9 ± 0.0b4.1 ± 0.1a7.1 ± 0.0b3.6 ± 0.4a1.2 ± 0.2b1.2 ± 0.1b3.7 ± 0.3c3.9 ± 0.1a***
C6MHβ-myrcene9930.1 ± 0.0c0.5 ± 0.0c0.3 ± 0.0c19.0 ± 0.9b0.5 ± 0.1d5.0 ± 0.1b4.8 ± 0.9ab6.1 ± 0.4a***
C7MHα-phellandrene10110.0 ± 0.0 cd0.4 ± 0.0 bc0.0 ± 0.0 bc0.4 ± 0.0 b0.1 ± 0.1 e0.1 ± 0.0 de0.1 ± 0.0 a0.2 ± 0.0 a***
C8Otho-cymene10320.2 ± 0.1 c0.7 ± 0.1 ab0.9 ± 0.1 bc0.9 ± 0.1 ab0.3 ± 0.1 d0.3 ± 0.1 cd0.5 ± 0.1 d0.6 ± 0.1 a***
C9Othp-cymene10330.1 ± 0.1 ac0.2 ± 0.1 a0.4 ± 0.1 ab0.4 ± 0.1 a0.1 ± 0.1 bc0.1 ± 0.1 a0.2 ± 0.1 c0.3 ± 0.1 a
C10MHlimonene10360.5 ± 0.1 f3.8 ± 0.1 f2.7 ± 0.0 b4.2 ± 0.1 b1.4 ± 0.2 e4.7 ± 0.1 d1.1 ± 0.2 c4.9 ± 0.3 a***
C11MHɣ-terpinene10660.2 ± 0.0 d0.5 ± 0.0 c0.3 ± 0.0 c0.9 ± 0.0 a0.1 ± 0.0 e0.3 ± 0.0 d0.3 ± 0.0 bc0.4 ± 0.0 b***
C12MHα-terpinolene10940.0 ± 0.0 c0.2 ± 0.0 d0.0 ± 0.0 bc0.2 ± 0.0 b0.1 ± 0.0 d0.2 ± 0.0 d0.1 ± 0.0 a0.1± 0.0 a***
C13Oth2-nonanaone10990.1 ± 0.0 cd0.3 ± 0.0 b0.1 ± 0.0 bc0.6 ± 0.1 a0.1 ± 0.0 d0.3 ± 0.0 b0.4 ± 0.1 a0.5 ± 0.1 a***
C14Othnonanol11090.2 ± 0.1 ac0.4 ± 0.1 d0.3 ± 0.1 bcd0.1 ± 0.1 a0.1 ± 0.1 ab0.5 ± 0.1 cd0.4 ± 0.1 ad0.3 ± 0.1 ac
C15Othnonanal11140.0 ± 0.0 a0.2 ± 0.0 b0.0 ± 0.0 a0.0 ± 0.0 a0.0 ± 0.0 a0.2 ± 0.0 b0.0 ± 0.0 a0.0 ± 0.0 a
C16OMmenta-2-en-1-ol <E-p>11370.0 ± 0.0 ac0.1 ± 0.0 ac0.0 ± 0.0 bc0.0 ± 0.0 ab0.0 ± 0.0 c0.0 ± 0.0 c0.0 ± 0.0 c0.0 ± 0.0 a
C17OMborneol11870.1 ± 0.0 ab0.2 ± 0.0 ab0.1 ± 0.0 ab0.1 ± 0.0 a0.0 ± 0.0 ab0.1 ± 0.0 c0.1 ± 0.0 ab0.1 ± 0.0 b***
C18OMterpinen-4-ol12071.3 ± 0.1 ce1.4 ± 0.0 b1.4 ± 0.2 de1.7 ± 0.1 a0.5 ± 0.0 e0.9 ± 0.0 cd1.5 ± 0.2 cd1.5 ± 0.0 c***
C19OMα-terpineol12090.4 ± 0.1 d0.4 ± 0.1 cd0.9 ± 0.1 c0.5 ± 0.1 a0.3 ± 0.1 c0.5 ± 0.1 b0.5 ± 0.1 e0.6 ± 0.1 ab***
C20OMbornyl acetate12976.2 ± 0.2a7.0 ± 0.1b4.9 ± 0.3b1.5 ± 0.1c2.8 ± 0.1b1.9 ± 0.1f1.9 ± 0.0d0.8 ± 0.0e***
C21Oth2-undecanone13011.6 ± 0.0 d2.0 ± 0.0 bc0.9 ± 0.1 bc2.1 ± 0.1 b1.6 ± 0.3 cd2.0 ± 0.1 c1.9 ± 0.2 a2.5 ± 0.2 b***
C22OMα-terpinyl acetate13580.2 ± 0.0 d0.3 ± 0.0 f0.0 ± 0.0 c0.0 ± 0.0 e1.0 ± 0.1 a1.5 ± 0.1 c0.3 ± 0.0 a0.5 ± 0.1 b***
C23SHα-copaene13860.9 ± 0.0 b0.8 ± 0.0 c1.2 ± 0.0 c0.7 ± 0.0 a0.3 ± 0.0 d0.6 ± 0.0 e0.5 ± 0.0 g0.4 ± 0.1 f***
C24SHβ-bourbonene13950.0 ± 0.0 ab0.1 ± 0.0 de0.0 ± 0.0 ac0.0 ± 0.0 e0.1 ± 0.0 bcd0.1 ± 0.0 ce0.0 ± 0.0 a0.0 ± 0.0 a***
C25SHβ-elemene13981.1 ± 0.0 d1.2 ± 0.0 cd1.5 ± 0.0 f1.0 ± 0.0 e1.4 ± 0.1 a1.2 ± 0.0 c1.5 ± 0.0 ef1.3 ± 0.0 b***
C26SHβ-caryophyllene14358.8 ± 0.1e10.7 ± 0.2d10.8 ± 0.5f13.2 ± 0.9f24.9 ± 1.0c17.2 ± 0.6b25.4 ± 0.8b21.7 ± 0.1a***
C27SHβ-cubebene14440.8 ± 0.1 b0.8 ± 0.0 ab0.4 ± 0.1 ab0.8 ± 0.0 ab0.9 ± 0.1 ab0.9 ± 0.1 ab0.8 ± 0.1 a1.6 ± 0.9 ab
C28Othisoamyl benzoate14550.8 ± 0.1 ab0.7 ± 0.1 ab1.0 ± 0.1 a0.6 ± 0.1 b0.8 ± 0.1 ab0.6 ± 0.1 ab0.5 ± 0.1 c0.6 ± 0.1 b***
C29SHα-caryophyllene14722.1± 0.0 f2.0 ± 0.0 e2.3 ± 0.1 h1.7 ± 0.1 g3.3 ± 0.0 a2.5 ± 0.1 b3.5 ± 0.0 d2.8 ± 0.0 c***
C30SHallo-aromadendrene14770.8 ± 0.0 a0.7 ± 0.0 b0.7± 0.0 b0.5 ± 0.0 b0.6 ± 0.0 a0.7 ± 0.0 c0.6 ± 0.0 c0.4 ± 0.0 d***
C31SHɣ-muurolene14892.6 ± 0.0 a1.9 ± 0.0 d2.8 ± 0.1 c1.2 ± 0.1 c1.4 ± 0.1 b1.7 ± 0.1 e1.4 ± 0.0 g1.0 ± 0.0 f***
C32SHgermacrene D14982.5 ± 0.6 bc2.8 ± 0.1 d1.2 ± 0.3 d4.0 ± 0.1 d4.2 ± 0.3 cd4.2 ± 0.3 b4.2 ± 0.5 a3.2 ± 0.2 b***
C33SHα-muurolene15141.9 ± 0.2 b1.7 ± 0.0 bc2.5 ± 0.0 cd1.2 ± 0.0 d1.8 ± 0.1 a1.6 ± 0.0 cd1.7 ± 0.1 e1.4 ± 0.0 d***
C34SHβ-bisabolene15201.1 ± 0.0 a0.8 ± 0.0 a0.4 ± 0.0 a0.8 ± 0.4 a0.4 ± 0.4 a0.8 ± 0.2 a0.6 ± 0.0 a0.6 ± 0.0 a*
C35SHɣ-cadinene15300.6 ± 0.1 b0.8 ± 0.1 b0.9 ± 0.0 b0.4 ± 0.0 c1.1 ± 0.0 a0.8 ± 0.0 b0.7 ± 0.1 c0.6 ± 0.0 b***
C36SHδ-cadinene15348.5 ± 0.0a4.7 ± 0.1b3.6 ± 0.4b3.0 ± 0.3b4.8 ± 0.2a4.9 ± 0.1b4.6 ± 0.3a3.3 ± 0.1c***
C37SHcalamenene15402.3 ± 0.0 a1.6 ± 0.2 c1.5 ± 0.0 b0.7 ± 0.0 b0.9 ± 0.0 a1.0 ± 0.0 d0.9 ± 0.0 d0.6 ± 0.0 e***
C38OSelemol15650.8 ± 0.0 c0.2 ± 0.0 f0.5 ± 0.0 a0.5 ± 0.0 d0.4 ± 0.0 e0.6 ± 0.0 b0.2 ± 0.0 e0.3 ± 0.0 g***
C39OScaryophyllene oxide15738.7 ± 0.2c5.7 ± 1.2c7.9 ± 0.7bc4.5 ± 0.4d9.9 ± 0.5a7.8 ± 0.1ab7.5 ± 0.6c7.8 ± 1.0cd**
C40OShumulene epoxide II16330.9 ± 0.1 bc0.6 ± 0.1 df0.7 ± 0.1 b0.2 ± 0.1 ef0.9 ± 0.1 a0.8 ± 0.1 bd0.5 ± 0.1 cde0.5 ± 0.1 f***
C41OScubenol16483.0 ± 0.1 a1.2 ± 0.0 d2.0 ± 0.1 cd0.7 ± 0.1 b1.0 ± 0.1 a1.4 ± 0.0 bc1.2 ± 0.1 e0.8 ± 0.0 f***
C42OSepi-Cadinol16632.7 ± 0.0 c1.8 ± 0.0 cd2.8 ± 0.1 de1.3 ± 0.1 e2.8 ± 0.1 a2.5 ± 0.1 b2.7 ± 0.1 e2.3 ± 0.0 de***
C43OSδ-cadinol16682.3 ± 0.1 b1.1 ± 0.1 cd1.3 ± 0.2 cd0.6 ± 0.0 d1.5 ± 0.1 a1.4 ± 0.1 bc1.2 ± 0.1 cd0.9 ± 0.2 e***
C44OSα-cadinol16798.0 ± 0.1cd5.0 ± 0.9bc8.0 ± 0.2bc3.7 ± 0.3d11.9 ± 0.4a7.9 ± 0.1ab7.6 ± 0.8bc9.5 ± 3.2bc**
C45OSaromadendrene oxide16841.0 ± 0.1 a0.6± 0.0 a1.1 ± 0.4 a0.3 ± 0.0 a0.5 ± 0.0 a0.5 ± 0.0 a0.4 ± 0.0 b0.3 ± 0.0 b**
C46OSbisabolol17049.7 ± 0.3ab2.6 ± 0.0b1.2 ± 0.0a1.4 ± 0.1ab1.3 ± 0.1ab1.5 ± 0.0c1.1 ± 0.1ab1.2 ± 0.0d***
C47Othbenzyl benzoate17970.4 ± 0.0 ab0.2 ± 0.0 c0.5 ± 0.0 b0.3 ± 0.1 c0.8 ± 0.2 ab0.9 ± 0.0 a0.4 ± 0.0 b0.4 ± 0.0 ab***
Total (sum of C1-47)88.0 ± 0.691.1 ± 1.494.8 ± 0.694.0 ± 0.390.9 ± 0.890.9 ± 0.494.9 ± 0.396.6 ± 3.0
Terpene Classes
MHmonoterpene hydrocarbons5.8 ± 0.3 c28.0 ± 0.1 b28.2 ± 0.3 b42.8 ± 2.6 a6.4 ± 1.1 d19.1 ± 1.0 c17.4 ± 1.5 c27.1 ± 1.2 a***
SHsesquiterpene hydrocarbons33.8 ± 0.5 c30.5 ± 0.6 c29.8 ± 1.3 d29.1 ± 1.4 d45.9 ± 0.9 a37.9 ± 0.9 a46.5 ± 1.7 a38.3 ± 0.4 b***
OMoxygenated monoterpenes8 ± 0.3 a9.3 ± 0.1 a7.3 ± 0.5 a3.8 ± 0.2 a4.6 ± 0.4 a4.9 ± 0.2 a4.3 ± 0.3 a3.4 ± 0.1 a
OSoxygenated Sesquiterpenes37.0 ± 0.5 bc18.8 ± 1.1 c25.4 ± 0.6 bc13.2 ± 1.0 d30.3 ± 0.9 a24.3 ± 0.4 b22.4 ± 0.9 c23.6 ± 4.4 e***
OthOthers 3.3 ± 0.1 d4.5 ± 0.2 cd4.1 ± 0.2 bc5.1 ± 0.0 ab3.7 ± 0.5 d4.8 ± 0.1 cd4.3 ± 0.1 bc5.1 ± 0.1 a
The means followed by different letters (in italic) within the same row were significantly different (p < 0.05). The main compounds name and levels are written in bold. Significance effect of the interaction between the altitude and season as determined by ANOVA, p-values: *** p < 0.001; ** p < 0.01; * p < 0.05.
Table 2. Evaluation of the antioxidant activities of Pistacia lentiscus essential oils using FRAP, ABTS, and DPPH tests.
Table 2. Evaluation of the antioxidant activities of Pistacia lentiscus essential oils using FRAP, ABTS, and DPPH tests.
Antioxidant Assays
Essential Oil SamplesFRAP (mg TE/g EO)ABTS (mg TE/g EO)DPPH (mgTE/g EO)
MountainWinter11.3 ± 0.1 c0.28 ± 0.01 b0.12 ± 0.0 a
Spring9.6 ± 0.3 a0.32 ± 0.0 c0.44 ± 0.01 a
Summer22.3 ± 0.2 c0.09 ± 0.01 e0.08 ± 0.01 b
Autumn14.3 ± 0.1 a0.24 ± 0.01 a0.11 ± 0.0 a
LittoralWinter6.4 ± 0.9 c0.09 ± 0 c0.06 ± 0.0 c
Spring6.6 ± 1.1 b0.23 ± 0.01 e0.06 ± 0.0 d
Summer15.7 ± 0.5 d0.30 ± 0.01 a0.08 ± 0.0 b
Autumn14.8 ± 0.4 b0.13 ± 0.01 d0.06 ± 0.0 c
Altitude * Season Interaction*********
Vitamin C (mg TE/g vit C)1080 ± 144894 ± 881340 ± 192
Means followed by different letters (in italic) within the same row were significantly different (p < 0.05). Significance effect of the interaction between the altitude and season as determined by ANOVA, p-values: *** p < 0.001.
Table 3. Highest correlation coefficients between components of EOs and antioxidant activities.
Table 3. Highest correlation coefficients between components of EOs and antioxidant activities.
Major CompoundsFRAPABTSDPPH
α-Pinene0.370.090.39
ß-Pinene0.62 **−0.090.17
β-Myrcene0.070.02−0.33
Bornyl acetate−0.180.270.56 **
β-Caryophyllene−0.07−0.18−0.67 ***
δ-Cadinene−0.54 **0.330.15
Caryophyllene oxide−0.24−0.45 *−0.46 *
α-Cadinol−0.11−0.54 **−0.56 **
Bisabolol−0.63 ***0.360.43
Oxygenated monoterpenes−0.210.300.55 **
Oxygenated sesquiterpenes−0.17−0.40−0.27
Significant correlations: *** p < 0.001; ** p < 0.01; * p < 0.05.
Table 4. FAME analysis of the acyl lipids extracted from lentisk leaves. The levels of each compound are expressed as chromatographic relative peak surface measured as total ion current (area % TIC). The extraction yield for each lipid fraction was similar: 1.95 ± 0.05 g/100 g of dry matter.
Table 4. FAME analysis of the acyl lipids extracted from lentisk leaves. The levels of each compound are expressed as chromatographic relative peak surface measured as total ion current (area % TIC). The extraction yield for each lipid fraction was similar: 1.95 ± 0.05 g/100 g of dry matter.
Mountain (Site M)Littoral (Site L)Altitude * Season
Interaction
Fatty AcidCorresponding
acyl
WinterSpringSummerAutumnWinterSpringSummerAutumn
Methyl Ester (ME)
Myristic acid MEC14:06.6 ± 1.7 ac4.0 ± 0.2 bc1.8 ± 0.1 ac3.4 ± 0.1 ab2.5 ± 0.1 bc3.2 ± 0.1 c2.7 ± 0.3 a2.7 ± 0.5 d***
Hexadecenoic acid MEC16:11.2 ± 0.3 c0.9 ± 0.1 b1.2 ± 0.1 bc1.1 ± 0.1 a0.5 ± 0.1 b1.1 ± 0.1 ab1.3 ± 0.1 bc1.5 ± 0.2 bc***
Palmitic acid MEC16:024.8 ± 2.3 c28.3 ± 0.3 c27.4 ± 0.5 c30.3 ± 0.9 a26.3 ± 1.1 bc31.5 ± 0.4 ac31.4 ± 1.0 ab31.7 ± 2.0 a
Margaric acid MEC17:00.7 ± 0.1 c0.8 ± 0.1 b0.8 ± 0.1 bc0.9 ± 0.1 a0.6 ± 0.1 b0.9 ± 0.1 ab0.9 ± 0.1 ab1.1 ± 0.1 ab**
Linoleic acid MEC18:220.1 ± 1.7 ab20.0 ± 0.2 c28.2 ± 8.7 a13.9 ± 0.6 c26.2 ± 0.9 ab22.7 ± 0.5 bc9.7 ± 0.2 c14.1 ± 1.2 bc***
Oleic acid MEC18:12.8 ± 0.2 c5.7 ± 0.4 bc15.3 ± 8.3 d10.8 ± 0.3 c14.4 ± 0.3 bc10.4 ± 0.6 ab30.2 ± 1.2 c15.3 ± 1.5 a**
Linolenic acid MEC18:335.5 ± 1.8 d34.0 ±1.2 bc18.5 ± 0.2 a31.8 ± 1.1 c23.4 ± 0.2 e21.7 ± 0.6 e17.6 ± 0.7 ab27.2 ± 2.8 e***
Stearic acid MEC18:01.9 ± 0.2 bc2.3 ± 0.1 ab2.6 ± 0.1 c2.6 ± 0.1 a2.1 ± 0.1 bc2.2 ± 0.2 ab2.8 ± 0.1 bc2.6 ± 0.3 a
Docosanoic acid MEC22:00.9 ± 0.1 ac1.2 ± 0.1 c1.1 ± 0.1 bc1.2 ± 0.1 a0.9 ± 0.1 bc1.3 ± 0.1 c1.2 ± 0.2 ac1.2 ± 0.2 ab
Total identified94.3 ± 3.096.7 ± 0.796.8 ± 0.795.9 ± 0.696.9 ± 0.394.8 ± 0.997.8 ± 0.297.3 ± 0.3
C14; C16; C17; C18; C22Mean (ΣSFA) 34.8 ± 1.0 d36.6 ± 0.4 d33.6 ± 0.7 d38.3 ± 1.1 a32.4 ± 1.2 cd39.0 ± 0.7 bd39.0 ± 1.5 ab39.6 ± 2.9 abc**
C16:1; C18:1Mean (ΣMUFA)3.9 ± 0.5 c6.6 ± 0.4 ac16.5 ± 8.3 d11.9 ± 0.4 bc14.9 ± 0.2 ac11.5 ± 0.5 ab31.5 ± 1.2 c16.8 ± 1.7 a*
C18:2; C18:3Mean (ΣPUFA)55.5 ± 2.7 b54.0 ±1.1 bc46.7 ± 8.6 a45.7± 0.5 bd49.6 ± 0.9 bd44.3 ± 1.1 cd27.3 ± 0.5 bd41.2 ± 1.8 d**
Means followed by different letters (in italic) within the same row were significantly different (p < 0.05). Significance effect of the interaction between altitude and season as determined by ANOVA, p-values: *** p < 0.001; ** p < 0.01; * p < 0.05.
Table 5. Correlation coefficients between major components of EOs and FAMEs from the leaves of Pistacia lentiscus.
Table 5. Correlation coefficients between major components of EOs and FAMEs from the leaves of Pistacia lentiscus.
CompoundsC18:3C18:2C18:1C16:0C14:0
β-Caryophyllene−0.58 *−0.190.77 ***0.48 *0.31
β-Myrcene0.15−0.55 ***0.49 *0.76 ***0.08
α-Pinene0.10−0.08−0.180.210.009
Bornyl acetate0.290.44 *−0.57 *−0.68 **0.3
δ-Cadinene0.150.36−0.42 *−0.340.39
Bisabolol0.81 ***0.29−0.83 ***−0.370.80 ***
ß-Pinene−0.26 *−0.100.180.21−0.41 *
Caryophyllene oxide−0.140.400.09−0.46 *−0.32
α-Cadinol−0.30 0.340.21−0.33−0.48 *
Significant correlation (in bold): *** p < 0.001; ** p < 0.01; * p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sehaki, C.; Jullian, N.; Choque, E.; Dauwe, R.; Fontaine, J.X.; Molinie, R.; Ayati, F.; Fernane, F.; Gontier, E. Profiling of Essential Oils from the Leaves of Pistacia lentiscus Collected in the Algerian Region of Tizi-Ouzou: Evidence of Chemical Variations Associated with Climatic Contrasts between Littoral and Mountain Samples. Molecules 2022, 27, 4148. https://doi.org/10.3390/molecules27134148

AMA Style

Sehaki C, Jullian N, Choque E, Dauwe R, Fontaine JX, Molinie R, Ayati F, Fernane F, Gontier E. Profiling of Essential Oils from the Leaves of Pistacia lentiscus Collected in the Algerian Region of Tizi-Ouzou: Evidence of Chemical Variations Associated with Climatic Contrasts between Littoral and Mountain Samples. Molecules. 2022; 27(13):4148. https://doi.org/10.3390/molecules27134148

Chicago/Turabian Style

Sehaki, Chabha, Nathalie Jullian, Elodie Choque, Rebecca Dauwe, Jean Xavier Fontaine, Roland Molinie, Fadila Ayati, Farida Fernane, and Eric Gontier. 2022. "Profiling of Essential Oils from the Leaves of Pistacia lentiscus Collected in the Algerian Region of Tizi-Ouzou: Evidence of Chemical Variations Associated with Climatic Contrasts between Littoral and Mountain Samples" Molecules 27, no. 13: 4148. https://doi.org/10.3390/molecules27134148

APA Style

Sehaki, C., Jullian, N., Choque, E., Dauwe, R., Fontaine, J. X., Molinie, R., Ayati, F., Fernane, F., & Gontier, E. (2022). Profiling of Essential Oils from the Leaves of Pistacia lentiscus Collected in the Algerian Region of Tizi-Ouzou: Evidence of Chemical Variations Associated with Climatic Contrasts between Littoral and Mountain Samples. Molecules, 27(13), 4148. https://doi.org/10.3390/molecules27134148

Article Metrics

Back to TopTop