Abstract
BACKGROUND:
The present work stands as an endeavor to uncover the ideal harvesting time of leaves in which they exhibits the maximum contents of bioactive molecules such as essential fatty acids, tocopherols and carotenoids.
METHODS:
A large scale investigation was carried out for the leaves of Pistacia atlantica involving a large number of populations collected over a period of four months during the growing season. The antioxidant activity was evaluated using both DPPH and β-carotene assays. The chemical percentage variability of the fatty acids was investigated using statistical analysis methods (Agglomerative Hierarchical Clustering “AHC”, also cited as CAH).
OBJECTIVE:
During the growth period, the effects of harvesting date on the total contents of lipids, tocopherols, carotenoids, fatty acids composition and also the antioxidant activities of the lipids were investigated.
RESULTS:
The content of myristic acid and other saturated fatty acids increased during leaf development, while linoleic, linolenic acids and unsaturated fatty acids decreased. The highest percentages of both linolenic C18:3 (27.25±5.92%) and linoleic acids C18:2 (17.68±3.80%) were obtained for the month of May at the first stage of leaves development (young leaves), but higher percentage levels of C18:1, were obtained for both consecutive months of August & September (28.83±6.50%; 27.79±8.63%, respectively) at intermediate developing stage. The lipids, tocopherols, carotenoids contents and the saturated and unsaturated fatty acids were dependent on the harvest time. The antioxidant activity showed higher powers at the first developing stage (May). Two main clusters and two sub-clusters of the fatty acids were distinguished and were also depending clearly on the period of the collection.
CONCLUSIONS:
The results showed that the FA were dependent on the period of collection of the leaves. The main result of this study illustrate the nutritional potential (richness in MUFA such as C18:1, 2, 3) of the oil of P. atlantica leaves, which can provide opportunities for rational exploitation in the food industries or for medicinal purposes.
Introduction
One of the objectives of the food industry is to seek new resources of oil that preferably presents nutritional values. For human health, fats and particularly vegetable oils are considered as an important source of energy when glucose is not available. Lipids play an important role in food and nutrition serving as a source of essential fatty acids, vitamins, carotenoids and antioxidants compounds. As part of lipids, unsaturated fatty acids UFA are known for their multiple dietary, pharmaceutical and clinical benefits. The naturally occurring fatty acid: conjugated linoleic acid (C18:2), presents different properties such as anti-atherogenic and anti-carcinogenic [1]. In the same context and besides its role in reducing cardiovascular disease risk, studies evaluated the direct anti-atherogenic and anti-inflammatory effects of dietary α-linolenic acid (C18:3) polyunsaturated fatty acids (PUFA) in addition to their lipid-lowering effects [1, 2]. Vitamin E in the form of tocopherols or tocotrienols is an example of natural antioxidants found in lipids fractions [3]. Antioxidants serve as protective role by scavenging free radicals such as reactive oxygen species that can lead to the damage of cell membrane. Carotenoids are also of commercial significance due their application as food colorants.
In recent years, ethnobotanical studies performed in several countries indicated that Pistacia species play a vital role in the nutrition and agricultural economy of many communities, living in arid and semi-arid regions [4, 5]. The Pistacia nuts are edible, with abundance in unsaturated fatty and other bioactive compounds. For their nutritional profile, nuts have been reported to reduce the risk of coronary heart disease, blood pressure, cancer, inflammation, gallstones, and diabetes [4]. Pistachio acquired its economic importance from being one of the most favored nuts around the world. Moreover, pistachio has a high oil content consisting of about 50–62% of its weight [6]. Pistachio nut particularly is rich in protein, dietary fiber, potassium, magnesium, vitamin K, γ-tocopherol, and a number of phytochemicals. With renewed interest in relevance of ethnobotanical uses, several properties such as antioxidant, wound healing, antimicrobial, antiviral, anticholinesterase, anti-inflammatory, antinociceptive, antidiabetic, antitumor, antiancer, antihyperlipidemic, antiatherosclerotic, and hepatoprotective have been evaluated [4]. Ongoing studies on pistachio extracts has demonstrated lately the confirmed effects of its kernel on cellular viability, intracellular reactive oxygen species (ROS) production and cell death in MCF-7 breast cancer cells [7]. Although pistachio oils have different nutritional values according to the nut origin source, pistachio cold pressed oil has generally a great amount of unsaturated fatty acids and bioactive phytochemicals such as phenolics, stilbenes, carotenoids, phytosterols, and proanthocyanidins [6, 8]. Because of its special chemical composition, pistachio oil plays an important role in the industrial production of cosmetics, with positive gentle effects on skin and hair. Moreover, it can contribute significantly to the pharmaceuticals industry, as it is rich in bioactive ingredients, having both antioxidant and anti-inflammatory properties [6].
The Pistacia atlantica plant belongs to the Anacardiaceae family. Its origins are from the Mediterranean and the Middle East countries. The French botanist R. L. Desfontaines discovered it in 1798 [9]. This plant grows in the Algerian steppe and desert areas. It grows in clean forest. The Pistacia atlantica tree has significant importance as follows: it resists the high temperature and it does not need much water; the Pistacia atlantica tree protects the steppe region and humid the areas because it is drought resistant; it preserves the soil, stabilizes and cools it.
The Pistacia atalantica tree produces a huge quantity of fruits, leaves and galls utilized traditionally for culinary and medicinal purposes and not yet used in industrial applications. All these organs of the tree are rich in natural bioactive compounds. P. atlantica is considered a functional food, because it is rich on nutrients, including sugars, fatty acids, amino acids, minerals, vitamins, and antioxidants [5]. The fruits contain an important amount of oil with predominance of unsaturated fatty acids (73–83%) [10–12], and also the oil is rich in tocopherol compounds [13]. Frequent fruit consumption has shown a reduced risk of cardiovascular disease and lowering effects of total and LDL cholesterols as well as properties to eliminate stomach diseases [14]. Hispolone and luteolin are naturals antioxidants that has been isolated from the fungus Inonotus hispidus growing on Pistacia atlantica [15]. The different organs of Pistacia atlantica tree contain the essential oil with predominance of monoterpene hydrocarbon [16–23]. The essential oil of leaves of Pistacia atlantica is very important source of non-toxic natural antioxidants because of its significant antioxidant activity [24]. Also Pistacia atlantica leaves constitute a source of phenolic compounds with important antimicrobial, particularly antifungal, antioxidant, antidiabetic and antihypertensive [25]. The essential oil of the galls showed antimicrobial, antibacterial, antioxidant and low cytotoxic activities [26]. The amount of the phenolic content, the chemical composition of the essential oils of the Pistacia atlantica leaves and the antioxidants activities are dependent on seasons and vary the during stage of development of the leaves [22, 27].
Although the fatty acids composition, and the total contents of both tocopherols and carotenoids of Pistacia atlantica were recently reported by another team of our laboratory [28]. Nevertheless, this previous study was not conclusive since many parts were investigated at once (leaves, galls and fruits), using an extremely reduced number of samples (at number of four for each plant part). Actually, these samples were collected from the same location over three months of study period [28]. Upon these items, the variability factors were missing and hence the study results are considered preliminary (conclusions must be considered with precaution), and need to be further studied with more expanded regions and with larger number of samples over wider range of growth period.
Hence, this report aims to investigate at a large scale the change in the content in lipid fraction, the fatty acids composition, tocopherols, carotenoids content and the antioxidant activities of the lipids from leaves of Pistacia atlantica tree collected over an extended period of four months during the growing season and involving 22 different populations. Moreover, the antioxidant activitiy was determined using two different in vitro tests: DPPH and for the first time β-carotene bleaching assays (lipid peroxidation test). Data analysis using AHC was carried out in order to establish positive distinctive groups on the basis of the studied parameters. The present work stands as an endeavor to uncover the ideal harvesting time of leaves when rich, in bioactive molecules such as essential fatty acids, tocopherols and carotenoids.
Materials and methods
Chemicals
Absolute ethanol, methanol, n-hexane, chloroform, activated charcoal, sulfuric acid, acetic acid, ferric chloride, dichloromethane, Vitamin E and DPPH• (1,1-diphenyl, 2-picrylhydrazyl) were obtained from Sigma-Aldrich (Saint-Louis, Missouri, USA). Anhydrous sodium sulphate and sodium were purchased from Analar NORMAPUR Prolabo (France). Butylated hydroxyanisole (BHA) was purchased from Fluka chemie (Buchs, Switzerland). Acetic anhydride was obtained from Merck (Darmstadt, Germany). Ortho-phenanthroline was obtained from Biochem Chemopharma (Cosne-Cours-sur-Loire, France). All reagents and chemicals used were of analytical grade.
Plant material
The study was performed on Pistacia atlantica leaves obtained from six different trees belonging to three sites namely: Site 1 (33°48’7′′N and 2°53’19′′E with an altitude of 760 m), site 2 (33°45’2′′N and 2°48’06′′E with an altitude of 777 m) and site 3 (33°39’14′′N and 2°55’10′′E with an altitude of 810 m). The six different chosen populations were picked-up randomly. The samples were collected at the bottom of each tree, and at different locations all around the tree. All samples were collected regularly during growth (different phenological states) at the beginning of each month for a period of four months (May “04 samples”, August “06 samples”, September “06 samples” and November “06 samples”, in year 2014). All leaf samples were green, except those obtained in November (autumn) which were in yellow color. In May, the leaves were young and still growing bigger. In full developed stage (August/September) the leaves were fully expanded. Finally at last stage (November) the leaves were tending to yellow color before falling.
Although, extraction of lipids should be performed as soon as possible after the removal of tissues from the living organism, sometimes, nonpolar solvents, such as diethyl ether and n-hexane, do not easily penetrate the moist tissues (>8% moisture); therefore, effective lipid extraction does not occur. Therefore, reducing moisture content of the samples may facilitate lipid extraction [29]. In this context, and since the solvent employed was n-hexane, the collected leaf samples of the plants were air-dried in the shade at room temperature for two weeks, prior to lipid extraction.
Lipid extraction
The air dried leaves of Pistacia atlantica were milled into powder using a coffee grinder. Lipid extraction was performed with n-hexane using a Soxhlet extraction apparatus for 2–3 days. The lipid extract was dried by adding 2 g of anhydrous sodium sulphate, filtered and evaporated to dryness using a rotary evaporator (temperature set at 40°C), with a ratio of 1:5 (mass) activated carbon to initially air-dried leaves, chlorophyll removal was carried out by adsorption on activated charcoal by reflux using chloroform as a solvent for 2 hours. The obtained extract was cooled to room temperature then filtered. Chloroform was removed from the solution on a rotary vacuum evaporator at 40°C. The lipid extract was weighed to determine the lipid content. The lipid fractions were subjected to another operation which is waxes removal. The lipid extract were dissolved in absolute ethanol and stored in the refrigerator for 24 hours, then filtered and evaporated at 40°C. The lipid extract was stored at+4°C for further analysis.
Fatty acids composition
In order to determine the fatty acids (FA) composition, lipids were transformed into fatty acid methyl esters (FAME), the method consisted of refluxing the lipids (0.5 g) within 20 mL of 0.5% sodium methylate (NaOMe prepared by action of sodium on methanol) for 25 min. After cooling, 20 mL of distilled water was added. In the next step, methyl esters of FA were extracted with dichloromethane (DCM), and then washed with distilled water. Subsequently, the obtained solution was dried over anhydrous sodium sulfate Na2SO4 then evaporated under vacuum at 40°C. Prior to gas chromatography (GC) analysis and identification, the FAME were purified by chromatography over silica gel as stationary phase and DCM was used as a mobile phase.
Gas chromatographic (GC) analysis of FAME
FAME were analyzed on a gas chromatograph (Chrompack, CP 9002, Netherlands) equipped with a split-splitless mode injection system, flame ionization detector (FID) and a DB23 capillary column (30 m×0.32 mm i.d., 0.32μm film thickness). Oven temperature was in isotherm mode (250°C). Injector and detector temperatures were set at 250°C. Nitrogen was used as carrier gas at a flow rate of 1 mL/min. The injection volume was 1μL. The fatty acids were identified by comparison of their retention times with those of pure standards.
Gas chromatography/mass spectroscopy (GC/MS)
FAME were analyzed by GC/MS using a Shimadzu GC/MS-QP2010 ultra equipped with Auto sampler AOC-20i, Ion source: electronic impact High-performance Quadrupole Mass Filter. Separation of compounds was carried out in a DB-5 ms capillary (Length 30 m; diameter 0.25 mm; film thickness 0.25μm) using helium as the carrier gas (99.99% purity) with a flow rate of 1.2 mL/min. The oven temperature was programmed at 70°C for 5 min, then 4°C/min to 250°C, followed by temperature rate, 2°C/min up to 300°C which was maintained for 5 min.
Total tocopherol content
The total tocopherol content (TTC) of lipid fractions was determined using the spectrometric method of Emmerie, et al. [30] with some modifications. In order to draw a calibration curve relating the optical density to the concentration, a series of known concentration solutions of commercial vitamin E was prepared in ethanol. 1 mL of 1,10-phenanthroline reagent (0.4% in absolute ethanol) and 0.5 mL of ferric chloride reagent (0.12% in absolute ethanol) were added to 1 mL of each solution. The mixture was incubated for 5 minutes at room temperature. Absorbance was measured at 510 nm (Shimadzu UV/Vis 1601 apparatus) against a blank prepared using 1 mL of ethanol, 1 mL of 1,10-phenanthroline reagent and 0.5 mL of iron chloride reagent. The dewaxed lipid samples were treated with the same procedure described above. TTC in dewaxed lipids were calculated from the vitamin E calibration curve and results were expressed as mg α-tocopherol eq./g of lipid. All determinations of TTC were carried out in triplicate.
Total carotenoid content
Determination of total carotenoids (TCC) of lipid fractions was performed by the method described by Rutz, et al. [31], with modifications. A standard calibration curve was obtained from different concentrations of β-carotene solutions prepared in dichloromethane; the absorbance of each solution was determined at 460 nm. Quantification of carotenoids in dewaxed lipids was carried out by the same procedure. The results were expressed as mg β-carotene eq/g of lipids. All determinations of TCC were carried out in triplicate.
In vitro antioxidant activity
DPPH radical scavenging assay
DPPH assay was performed using the stable radical 2,2-diphenyl-1-picrylhydrazyl as described by Molyneux [32], with a few modifications. Different volumes of dewaxed lipid dilutions in absolute ethanol were added to 1 mL of 100μM solution of DPPH•. After 30 min of incubation at room temperature, the absorbance was read against a blank at 517 nm. The absorbance of the control was measured by replacing the sample with ethanol for each dilution volume. The percent scavenging activity I% was calculated as follows:
Where:
Ablank is the absorbance of the control reaction (DPPH with only ethanol) and Asample is the absorbance in the presence of the extract sample.
The antiradical activity was expressed as EC50 (mg/mL) which represents the concentration required to cause 50% initial DPPH• inhibition. Commercial vitamin E and BHA were used as reference standards. All experiments were performed in triplicate.
β-carotene bleaching assay
The β-carotene bleaching method was established for each fraction trough the method proposed by Sun, et al. [33] with some modifications. 1 mg of β-carotene was dissolved in 10 mL of dichloromethane. 5 mL of β-carotene solution was mixed with 20 mg of linoleic acid and 200 mg of Tween 80, then dichloromethane was evaporated under vacuum at 40°C and the resulting mixture was immediately diluted with 50 mL of distilled water, then the emulsion was vigorously shaken. 40μL of different solutions concentrations were added to 1 mL of β-carotene/linoleic acid emulsion. The absorbance of the preparations was measured at 490 nm before and after 120 min of incubation in a water bath at 50°C against negative control prepared by substituting the extract with an equal volume of methanol. The antioxidant activity (AA%) was performed using the following equation:
Where:
As120 and Ac120 are the absorbance measured for the test sample and control, respectively, after 120 min of incubation.
Ac0 is the absorbance for the control measured at zero time before the incubation.
The results were also expressed by AA50, the concentration providing 50% β-carotene bleaching inhibition). Vitamin E was used as positive control. All trials were performed in triplicate. All experiments were performed in triplicate.
Statistical analysis
Data analysis
Statistical analysis was carried out using XLSTAT Version 2016.02.27444. Experiments were conducted with six replicates (prior six randomly chosen tree populations). The results were expressed as mean±SD of three independent experiments (three leaves repetitions).
Raw data obtained were analyzed through analysis of variance (ANOVA) combined (as a posttest procedure) with Tukey’s post-hoc (HSD) multiple range tests. P-values of less than 5% (p < 0.05), were considered to be significant.
Cluster analysis
Cluster analysis was performed using Agglomerative Hierarchical Clustering AHC (Ward’s technique) with Euclidean distance measure. The calculus was performed using a set of individuals composed of 22 individuals (representing different collection sites and different months of collection) and the whole 11 fatty acids identified.
Results and discussion
Lipids extract
The lipid contents (range variation) obtained from air-dried Pistacia atlantica leaves harvested during several months are shown in Table 1 (the data are presented statistically). Dewaxed lipids content mean value increased gradually throughout leaf development, (Fig. 1); it seems that the maximum mean value was reached in November at the last stage of development (3.05% against 0.89% for full developing stage of leaves). In order to confirm the lipids yield evolution throughout this season, a statistical analysis involving ANOVA followed by Tukey’s tests were performed. In other words, the question was: is there significant difference between the mean yields for the four months? Indeed, results showed significant differences (95% of confidence at least) between the yields of May, July and November. The yield in September was not significantly different to both August and November. This close similarity could be explained by the fact that the September month is an intermediate state where the full developed phenological stage spread from August until the beginning of November. In this case, the lipids yields increased rapidly at the beginning of the phenological state May (young leaves), then increased very slowly to reach a steady state (a plateau) from the full until end of the developing phenological state (from August to November). Similar results were reported earlier, in which the total lipid content in Catharanthus roseus leaf, increased during expansion at full maturity. And with no marked changes occurred during aging [34].
Composition of the major fatty acids and lipid contents obtained from P. atlantica leaves
Composition of the major fatty acids and lipid contents obtained from P. atlantica leaves
SFA: saturated fatty acids; MUFA: monounsaturated fatty acids; PUFA: polyunsaturated fatty acids; UFA: unsaturated fatty acids; UFA/SFA: ratio “unsaturated/saturated” fatty acids; Mean: mean value relative to six different replicates (populations); sd: standard deviation. *UFA/SFA was not directly calculated by the ratio of values listed in the each column; rather that, this ratio was obtained by determining at a first step the ratio for each sample, in a second step, the tabulated ratio was determined using determined ratios in first step. Mean values followed by the same small letter within the same line are not significantly different (p < 0.05) according to Tukey’s (Multiple Range Test).

Variation in lipids content (mean value of each month) during development of Pistacia atlantica leaves. Vertical bars indicate standard deviation error (±). Mean values followed by the same small letter over experimental points are not significantly different (p > 0.05) according to Tukey’s (Multiple Range Test).
The low lipid content recorded at first stage of development can be attributed to the growing process of the young leaves, while at the last phenological state the biosynthetic activity of the lipids is at its highest rate. This final stage where leaf color is tending to yellow is considered a critical state since shortly afterwards leaf began to fall (aging and senescence). The gradual increase pattern may be explained by the fact that lipid content in Pistacia atlantica leaves might be regulated with the intervention of the enzymatic system of fatty acid synthetase (FA-Synthetase) which operates differently during this last stage. In other words, FA-Synthetase process was induced in some manner which resulted in higher production rate of lipids (accumulation). In the opposite way, at first stage, enzymes were qualified by a phenomenon called retro-inhibition in which they were more or less inactive. Hence, lipid production was at its lowest content. Accordingly, and in some recent works, it was also reported that the biosynthesis of lipids in fruits does not initiate at early stage of the development [35–37].
Changing weather conditions might also provide an explanation for the seasonal variation in the lipid content. Actually, lipids biosynthesis is known as sensitive to temperature. This parameter plays a crucial role in the composition of the oil FA by affecting the process regulating of FA-desaturase enzyme. As a consequence, in this study, lower temperatures at full and at last developing stages may be partially responsible of the higher lipids contents and especially in November where temperature reaches it lower value. This hypothesis was validated by Larkindale, et al. [38] who suggested that the lipid biosynthesis increases considerably at lower temperature.
The lipids accumulation may also be affected by other factors than temperature, such as light intensity exposure and daylight length, which are dependent on the season. It seems that decreasing day light and time exposure from May to November are both factors responsible for increasing lipid content through this growing period, in which lipid biosynthesis were more favorable.
The results of gas chromatography GC and gas chromatography coupled to the mass spectroscopy GC/MS analysis of the fatty acid methyl esters (FAME) of P. atlantica leaves obtained during different harvest months are shown in Table 1. In all samples, the range of main saturated fatty acids (SFA) detected were palmitic acid (C16:0 = 11.88–32.27%) followed by myristic acid (C14:0 = 2.94–18.55%), margaric acid (C17:0 = 0.0–10.43%) and stearic acid (C18:0 = 2.45–8.49%). On the other hand, the main unsaturated fatty acids (UFA) were oleic acid (C18:1 = 8.75–43.41%) followed by α-linolenic acid (C18:3 = 5.07–35.49%), linoleic acid (C18:2 = 6.38–21.48%), palmitoleic acid (C16:1 = 1.26–20.39%) and gondoic acid (C20:1 = 0.0–7.6%), although the others were found to be of low overall abundance.
For the overall samples, UFA fractions compositions showed the highest percentages ranging from 37.91 to 73.56% which were mostly represented by high amounts of poly unsaturated fatty acids (PUFA) in the range of (15.48–57.37%). Alternatively, the SFA were varying from 26.44 to 62.09%.
The results of ANOVA combined with Tukey’s tests showed significant difference between the contents of SFA in May and November. Same result was obtained with UFA. This means that there is a significant difference between the beginning and the end of the developing state for the leaves. This result is quite in agreement with the case of lipids yield, at the same extent. Ones again, the intermediate months were not significantly different, nevertheless, the mean variations values of SFA and UFA were evolving progressively: increasing for SFA and simultaneously decreasing for UFA (Fig. 2). Similar report finding was in accordance with obtained evolution paces for both SFA and UFA [34, 39]. Indeed, with the progress of leaf expansion, the saturated fatty acids (capric, lauric, myristic, palmitic and stearic) increased while unsaturated fatty acids (oleic, linoleic and stearic) were found to decrease when changes in lipid profile during growth and senescence of Catharanthus roseus leaf was investigated by Mishra, et al. [34].

Bars representation of the variation in unsaturated (UFA) and saturated (SFA) fatty acids compositions of the Pistacia atlantica leaves during development (mean values were considered). Vertical bars indicate standard deviation error (±). Mean values followed by the same letter over similar colored bars are not significantly different (p > 0.05) according to Tukey’s (Multiple Range Test).
For the first stage of leaf development (May), the average percentage of UFA was higher than SFA, and the difference was the highest one during the season, but after progressive evolution throughout the months until the last stage of development (November), the SFA and UFA percentages were converting to exactly same values (Table 1 & Fig. 2).
The identification of the FAME showed some clear variability of the compositions of the leaves during development. The mean percentage of May (first stage) was significantly different to the rest of the months (Fig. 3). Same results of previous analysis were confirmed for α-linolenic and linoleic acids. When mean values were considered, the total PUFA and its main components such as α-linolenic and linoleic acids showed the same accumulation gait ranging from 44.93±8.99% to 21.22±2.84% and from 27.25±5.92% to 12.53±3.75%, and from 17.68±3.80 to 8.69±1.67%, respectively, which decreased significantly from the first to the middle development stage (from May to August). Then it increased slightly towards the last development stage (November), or at least, remains practically steady until the end of the phenological state taking in consideration the non significant difference (p < 0.05) between the months of: August, September and November. It should be noted that PUFA are essential to human development and health. Moreover, they are highly important in cancer risk reduction, and also required to improve cardiovascular illnesses [40] and, due to their role in preventing inflammation and cerebro-vascular diseases, they are also studied as potential therapeutic agents for neuro-protective effects [41].

Bars representation of the variation in poly unsaturated fatty acids (PUFA), C18:3 and C18:2 compositions of the Pistacia atlantica leaves during development (mean values were considered). Vertical bars indicate standard deviation error (±). Mean values followed by the same letter over similar colored bars are not significantly different (p > 0.05) according to Tukey’s (Multiple Range Test).
It should be highlighted that α-linolenic acid was the dominant UFA at the beginning of leaf formation in May where it presented a maximum mean of 27.25±5.92%, then, its amounts approached closely those of linoleic acid for the remaining development stages. It is worth noting that α-linolenic acid shows many protective effects (cardiovascular-protective, anti-cancer, neuro-protective, anti-osteoporotic, anti-inflammatory, and antioxidative effects) against a myriad of diseases [42]. Oleic acid was the main MUFA in all samples. It presented contrasting accumulation pace compared to α-linolenic and linoleic acids (Fig. 4). It increases starting from May to reach its maximum value for a long middle development stage period (37.60±28.83%, 43.41±27.79% in August and September, respectively); after that its decrease was once again towards its practical initial value at last development stage when leaves were of autumn colour. This analysis was confirmed by the Tukey’s test which has confirmed the presence of two similar pairs of different groups: May/November and August/September.

Bars representation of the variation in saturated fatty acids (SFA) and poly unsaturated fatty acids (PUFA) compositions of the Pistacia atlantica leaves during development (mean values were considered). Vertical bars indicate standard deviation error (±).
According to previous studies, the oxidative stability of P. atlantica fruit oil was associated with high content of oleic acid [43]. Inversely, for this study, leaves lipids were associated with higher content of linolenic acid and especially in May. The hight contents of both linoleic and linolenic acids at the beginning of the season might be attributed to the increase in the activity of desaturase enzyme at this specific period. In a previous work which dealt with leaves part, it was demonstrated that the linolenic acid content of the leaves extracts decreased over heat stress [38]. This was in agreement with the results of this report in which the linolenic acid content decrease significantly in summer period (August) and then remained significantly steady in September and even in November which is not expected and where temperature decreased at the end of the phonological state. The alterations (deterioration of membranes) that occur within cells of senescing leaves at last phonological state might be responsible of the low contents of both linoleic and linolenic acids. Alterations in cell membrane composition and integrity reflect changes in membrane permeability. In addition, and similar to this last report, the loss of linolenic acid was correlated with significant increases of palmitic acid. In the present study, an increase in lipid saturation of leaves occurred during investigated months. This is primarily due to a decrease in linolenic acid and an increase in linoleic and palmitic acids. According to Larkindale, et al. [38] similar saturation increases have been observed in many plants. For a nutritional point a view comparison of main FA contents of investigated plant with a well known food rich nutrition plant date palm (Phoenix dactylifera L.) shows that investigated oil was more interesting. In fact, seed oil date palm was rich in oleic acid but its of content of linoleic and especially linolenic acid was very low [44].
Regarding the SFA, for the overall samples, palmitic acid showed its highest mean values at the last stage of leaf development (24.78±5.33%), followed by myristic acid detected with higher percentage (13.81±3.17%) at the same period. Indeed, for both components, it seemed that the mean accumulation increased gradually during leaf development. This observation was not confirmed by Tukey’s analysis for palmitic acid, since no significant differences was found between the means of the months. In this situation, the mean values of palmitic acid could be considered steady along the season.
At last, for the following components: C17:0, C18:0, C20:0 and C20:1, it seems that their mean values were practically steady over the season (Fig. 4).
We had noted with interest that the significant proportion of α-linolenic acid found in the leaves within first development stage (May) was very higher than those reported in previous published studies, thus, Samani, et al. [45], reported 1.09% for Atlas pistachio seeds from Iran. Trabelsi, et al. [46], registered 13.18% and 0.68% at the first and last stages of maturation for P. lentiscus fruit from Tunisia. Belyagoubi-Benhammou, et al. [47], obtained 0.44% for P. lentiscus fruit from Algeria. In fact, the obtained percentages of α-linolenic acid found in the leaves within first development stage (May) were similar to those of Chelghoum, et al. [28]. But, due to a reduced number of samples for this last report, no significant change of the whole reported FA and their different classes was noticed for the consecutive studied months of July, September and October. Similarly, in this study and except for the month of May, the remaining months showed no significance in the change of some specific FA components such as: C18:0, C18:2&3 and also for PUFA fraction.
The UFA/SFA ratio of unsaturated fatty acids to saturated fatty acids was not steady, it decreased from the maximum 2.19±0.44% observed in May to reach 1.04±0.32% in the last development stage. These values would supply to the oil beneficial properties in the first stage of development. As confirmed in this study, many reports which have analyzed either the entire lipid extract of photo-synthetic tissue or that of the chloroplast alone have shown a decline of more than 50% in the unsaturated/saturated esterified fatty acid ratio during senescence [34].
Finally, it should be noted that C12:0 was not reported existing in the leaf lipids by Chelghoum, et al. [28]. At the opposite side, both C20:0 and C20:1 are not identified in the current investigation, and which were reported by Chelghoum, et al. [28]. Except for C18:3, highest values for the following major components C18:0, 1, 2 were recorded in comparison with the previous study of Chelghoum, et al. [28].
In recent years it has become evident that the different forms of tocopherols (vitamin E isomers) have antioxidant properties. They play an important role in cancer disease prevention and health promotion Li, et al. [48]. Hence, the study of tocopherols in P. atlantica leaves is of a major interest. The results of total tocopherol content (TTC) range variations in the dewaxed lipid fractions of P. atlantica leaves are shown in Table 2. The range variation was satrting from 1.55 to 4.26 mg/g of lipids. This range was pretty similar to the one determined recently for leaves of P. lentiscus tree from Algeria (3.2–4.6 mg/g of lipids) [49], but notably higher than those recently published (0.21–1.53 mg/g of lipids) for the same part and the same tree [28]. This higher content might be explained by the fact that when the crude lipids were dewaxed, the tocopherol compounds were obtained with higher quantities. Nevertheless, the obtained range was higher than those obtained by Benalia, et al. [50] for the pumpkin Cucurbita pepo L. seeds oil (0.10–0.22 mg/g of lipids) and by Guenane, et al. [51] who also registered 0.34 and 0.51 mg/g of oil for mature and immature P. atlantica fruits, respectively. Furthermore it is very important to mention that the results were higher than those of fruit oil found in other plants of the anacardiaceae family such as P. terebinthus var. Chia growing wild in Turkey as reported by Matthäus, et al. [52] and which recorded only 0.39 and 0.51 mg/g of oil.
Total tocopherols (TTC) and total carotenoids (TCC) contents in P. atlantica leaves during development expressed as mg (α-tocopherol/β-carotene) equivalents/g lipids. Analysis of the differences between the categories performed with a confidence interval of 95%
Total tocopherols (TTC) and total carotenoids (TCC) contents in P. atlantica leaves during development expressed as mg (α-tocopherol/β-carotene) equivalents/g lipids. Analysis of the differences between the categories performed with a confidence interval of 95%
Mean: mean value relative to six different replicates (populations); sd: standard deviation. Mean values followed by the same small letter within the same line are not significantly different (p < 0.05) according to Tukey’s (Multiple Range Test).
As consequence, P. atlantica leaves seem to be rich in TTC. So, in order to evaluate the differences in TTC, the main previous factor was examined, which is harvest month (leaf development).
The obtained results showed significant differences in TTC according to the leaf development. It appears that the maximum of TTC is practically produced in September, towards last development stage, the TTC seemed to decrease. This observation was consolidated by Tukey’s test which showed significant difference between mean values for May and September (Table 2). Comparison of our results with other published data presented by Gourine, et al. [22] in which they reported that the seasonal variation showed that most of the main components of the essential oils reached their maximum values in September.
According to literature [53], tocopherols (Toc) play diversified roles in plant growth and physiological processes. Irrespective of form, Toc plays vital roles in development, signal transduction, interacting, or controlling phytohormonal regulation, senescence, and so on. The forms and levels of Toc are altered depending on plant growth stage and those perform different metabolic functions. It is also recognized that Toc plays a role in plants’ nutrient translocation, accumulation, and vascular systems and in carbohydrate translocation in various plant parts through phloem. Tocopherol is an antioxidant and thus has a clear role in plants’ tolerance to abiotic stress.
It has been recognized that a diet rich in carotenoids, brought by the regular consumption of fruits and vegetables suggests a protective effect against cancer and cardiovascular diseases [54, 55]. These compounds play an important role in human vision because of their antioxidant properties and can filter blue light from the eye macula [56, 57]. Tocopherols and carotenoids improve oil stability and thus oils naturally rich in these constituents are preferred [58, 59]. The results of total carotenoid contents (TCC) in the dewaxed lipid fractions obtained from the leaves of P. atlantica are summarized in Table 2. The whole range variation was from 0.1 to 0.6 mg/g of lipids, which were presenting important variations. Contrarily to TTC case, significant difference between first and last month (beginning and end of development stages) were determined. For middle development (August and September), TCC were not significantly different (Table 2). TCC mean value increased from the beginning stage (May) to reach a plateau (steady state) during full developing stage, and then decreased once again to reach a mean value very close but slightly higher than obtained at the beginning stage. Comparison with recent published work for both same tree part showed that obtained TCC range was clearly higher than reported range of 0.0186–0.0866 mg/g lipids [28]. Once again, this difference could be attributed to the pre-treatment of crude oil by the dewaxing step.
Antioxidant activity
To estimate the antioxidant capacity, two antioxidant assay methods were investigated which are DPPH and β-carotene bleaching methods. To our knowledge, no studies have been done on the impacts of the potential influences of the main factor (harvest month) on the total antioxidant activities of lipid fraction of P. atlantica leaves during development. Table 3 summarizes the results for the antioxidant properties (EC50 values) of P. atlantica leaves for both assays. The comparison of the different results showed that the antioxidant activities in all samples, at different regions over development stages had varied in the same way for both assays (DPPH and β-carotene bleaching methods).
Variation of antioxidant activity of lipid fractions of P. atlantica leaves using DPPH and β-carotene bleaching assays
Variation of antioxidant activity of lipid fractions of P. atlantica leaves using DPPH and β-carotene bleaching assays
sd: standard deviation. Mean values followed by the same small letter within the same line are not significantly different (p < 0.05) according to Tukey’s Multiple Range Test.
First, the Tukey’s test applied for both antioxidant assays (DPPH and β-carotene) showed significant difference between first month of May and the rest of the months. Explicitly, August, September and November were not significantly different (Table 3).
Regarding DPPH results, the values of EC50 found for the different samples during leaf development stages varied between 3.48 and 13.17 mg/mL. DPPH results were not far different from those obtained by Gourine, et al. [22] which reported EC50 values of leaf essential oils for P. atlantica ranging from 8.80–27.48 mg/mL. The antioxidant activities were significantly higher than those presented by Guenane, et al. [51] who obtained EC50 values for P. atlantica fruit oil, ranging from 17.1 to 51.3 mg/mL. In the same context, it useful to remind that, ethanolic extracts of Pistacia lentiscus berries growing in Algeria were investigated recently and showed strong in vitro antioxidant activities using DPPH, ABTS and reducing power assays [60]. Comparison with obtained DPPH activities with those of activities of a several varieties cultivars of date palm (Phoenix dactylifera L.) seed oils shows that investigated oil was by far more active than reported ones (46.42–77.58 mg/mL) [44].
According to different stages of leaf development, the decrease of EC50 mean values was significant from the beginning until the middle stage of leaf development (August); subsequently, the values of EC50 remained very close and steady (no significant difference). It is important to remind that the lowest antioxidant capacity is reflected by the highest EC50 value. In this case, the antioxidant activity in May was lower than the rest of the months. Lipid extracts were less active compared to selected synthetic antioxidants (vitamin E, BHA) which showed very strong anti DPPH activities (Table 3). Comparison of obtained DPPH results with the previous study which dealt with the same plant part is practically impossible since different concentrations of initial DPPH solutions were employed [28].
β-carotene bleaching assay was applied for assessing antioxidant power of P. atlantica lipids fraction. The bleaching of β-carotene is a reaction caused by free radicals resulting from the hydroperoxides formed from linoleic acid. The rate of β-carotene bleaching can be decelerate in the presence of antioxidants by neutralization of linoleate free radical and other free radicals resulted through the system [61]. The AA50 whole range variation values obtained with influence of the studied factor was ranged from 3.99 to 12.92 mg/mL (Table 3), it was practically the same range obtained with DPPH assay. The results of β-carotene bleaching test were very close to the range found for P. lentiscus leaves (Anacardiaceae) essential oil from Tunisia that was reported by Aissi, et al. [62] and who obtained a range variying from 5.57 to 11.20 mg/mL. The mean AA50 value of 11.11±1.07 mg/mL registered in May was in accordance with the range assessed by Bachrouch, et al. [63] who obtained a range values extending from 10 to 11 mg/mL for essential oils of P. lentiscus leaves collected at the flowering stage (May). The best capacity was recorded in August for the leaves (AA50 value 3.99 mg/mL) and the lowest capacity was observed for the same tree in fewer months earlier in May (AA50 value 12.92 mg/mL). According to AA50 mean values, the leaves exhibited their best antioxidant activity in August (AA50 mean value 6.98±1.62 mg/mL). But statistically, the lowest antioxidant activity was recorded for developed stages (from August to November) with no significant difference. Likewise to DPPH case study for different stages of leaf development, the decrease of EC50 mean values was significant from the beginning until the middle stage of leaf development (August); subsequently, the values of EC50 come very close and steady.
As for the lowest antioxidant capacity, the highest AA50 value was recorded for in the first development stage (May): 10.55±1.90 mg/mL. The obtained results of antioxidant capacity were substantially lower than those determined by the vitamin E, which presents a high capacity to prevent the bleaching of β-carotene.
For better understanding of similarities and/or differences between the studied samples in terms of FA compositions (Possible presence of region effect and/or month of collection on the FA composition), cluster analysis was performed using Agglomerative Hierarchical Clustering (AHC). The results of AHC schemed in the dendrogram of Fig. 5, showed two clearly distinguished clusters “or groups” of samples (cluster I: left cluster) and (cluster II: right one composed of two subgroups SG(1) & SG(2)). The presence of these groups and subgroups consolidated the initial theory of existing dissimilarity of the studied samples upon the variability of the chemical composition of the studied samples. A close examination of the first group (I) reveals that it enclosed whole samples collected for the consecutive months of August & September (middle developing stage: full expanded leaves). Alternatively, the second group (II) agglomerated mainly the samples of May (late stage of beginning of the formation of the young leaves) with those of November (last stage of losing leaves from the trees: yellow leaves). Moreover, the first subgroup SG(1) separates the samples of November from those of May which belongs this time to the SG(2). For SG(1), exception was made for one sample which belongs to September. So, for this sample it may be the presence of some over process for its phenological stage at this particular case. In order to better clarify and quantify the differences between the determined clusters, we have determined the range variations, the mean values and their standard deviation for all groups and subgroups. These results are summarized in Table 4.

Dendrogram obtained from a cluster analysis of 22 individuals (representing different regions and different months of collection) and 11 variables (fatty acids obtained from Algerian Pistacia atlantica leaves). Samples are clustered using Ward’s technique with a Euclidean distance measure. S1-4: samples collected in May, S5-10: samples collected in August, S11-16: samples collected in September, S17–22: are samples collected in November.
Range variation and mean values of the different clusters of the studied individuals obtained by CAH classification for the composition of the major fatty acids obtained from P. atlantica leaves
*UFA/SFA was not directly calculated by the ratio of values listed in the each column; rather that, this ratio was obtained by determining at a first step the ratio for each sample, in a second step, the tabulated ratio was determined using determined ratios in first step.
For the case study of main similarities between all clusters, the three clusters are practically presenting similar mean values in the composition of the following SFA: C12:0; C17:0; C18:0 and C20:0. Now, the differences or the dissimilarities between these clusters are discussed as follows. The first group was mainly differentiated with its high percentage of UFA C18:1ω9, with a mean and sd variation of 29.07±7.14% with comparison of SG(1) (15.97±2.98%) and SG(2) (12.36±4.47%). The SG(2), was characterized by relatively higher percentages of UFA C18:3ω3 (27.25±5.92% vs “15.98±4.50% and 11.22±4.18% ”) and C18:2ω6 (17.68±3.80% vs “12.38±3.44% and 9.83±2.30% ”) when compared with the rest of the subgroup 1 and group (I), respectively. This SG(2) was also characterized by lower mean percentage of SFA C14:0 (3.82% against 13.21% and 8.10%). Finally, the SG(1) is characterized by relatively high mean percentage of C14:0 (13.21%) and C16:0 (24.91%) (Table 4), but the lowest mean value for the minor UFA C20:1ω9 (0.75%).
Alternatively, if we consider the difference between the obtained clusters upon the different total chemical classes (SFA, MUFA, PUFA and USF/SFA) we can observe that the SG(2) is mainly characterized by the highest mean value of PUFA (44.93% vs “28.36% and 21.04% ”). For this SG, the UFA was again slightly higher than those of the rest of the clusters (68.16% vs “51.07% and 58.91% ”). At the opposite side, SG(1) is also characterized by lower mean percentage of SFA (31.82% vs 48.93% and 41.09%). Furthermore, the UFA/SFA for this SG was the highest one (2.19±0.44). For the first group (I), it was distinguished mainly by its highest content of MUFA (37.87%) and by its lowest content of PUFA (21.04%), when compared to the rest of the clusters. The SG(1) presented the highest obtained mean value of SFA (48.93%) and the lowest UFA/SFA value (1.08±0.31). Some similarities were detected for the two SG(1&2) referring to their practically similar mean contents of MUFA (22.71% and 23.25%).
The current study showed that the harvesting month has a significant effect on the present material, concerning the lipids yield, the composition of the fatty acids and total contents of tocopherols and carotenoids. The increasing in the lipid content during development of Pistacia atlantica leaves (from May to November) was clearly detected. The lipids were rich in different and very important UFA. The SFA increased during leaves development while UFA decreased rapidly from first month to reach a steady state. The highest antioxidant activities and the percentages of both C18:3 and C18:2 were obtained for the month of May at the first development stage of leaves (young leaves). This month is considered as an optimal period for harvesting the leaves, which provide the most interesting nutritional potential in the food industries or for medicinal purposes.
Footnotes
Acknowledgments
The authors also wish to acknowledge the DGRSDT (Direction Générale de la Recherche Scientifique et du Développement Technologique) of the MESRS (Ministère de l’Enseignement Supérieur et de la Recherche Scientifique) of Algeria for the financial support granted for our research laboratory.
Author contributions
M.I. Benguechoua collected the plant samples, performed experiments and optimization and drafted the manuscript. M. Benguechoua supervised experiments. N. Gourine corrected the draft manuscript, checked the results and discussion, and supervised the work. A.M.S. Silva provided analysis. M. Saidi and M. Yousfi designed the study, supervised the work.
Conflict of interest
All authors state that they have no conflict of interest to declare.
