Oil Content, Fatty Acid and Phenolic Profiles of Some Olive Varieties Growing in Lebanon

Olive growing in Lebanon plays an important role at both a social and economic level. Nevertheless, the quality of olive oil produced in the country is rarely addressed. In this study, oil content, fatty acid, and phenolic profiles were studied along four different ripening stages for 11 varieties of olives, including two clones of the local variety “Baladi,” in addition to nine foreign varieties (“Ascolana Tenera,” “Bella di Cerignola,” “Itrana,” “Jabaa,” “Kalamata,” “Nabali,” “Salonenque,” “Sigoise,” and “Tanche”). Oil content was determined using the Soxhlet method and Abencor system. Fatty acid composition was determined using a GC-FID, total phenols using spectrophotometry, and the phenolic profile using HPLC-DAD. Results showed that variety, fruit ripening and their interaction have a significant effect on the overall studied oil parameters. Among the studied varieties, “Kalamata” presented the higher oil content on dry matter (OCDM = 48.24%), “Baladi 1” the highest oil content on humid matter (OCHM = 27.86%), and “Tanche” the highest oil industrial yield (OIY = 19.44%). While “Tanche” recorded the highest C18:1 (71.75%), “Ascolana Tenera” showed the highest total phenols (TP = 539 mg GAE/Kg of oil), “Salonenque” the highest oleacein (121.57 mg/Kg), and “Itrana” the highest oleocanthal contents (317.68 mg/Kg). On the other hand, oil content together with C18:2 and C18:0 increased along ripening while C18:1, total phenols and the main individual phenols decreased. Although preliminary, this study highlights the good quality of olive oil produced from both local and foreign varieties growing in Lebanon and encourages further investigations on the characterization and authentication of Lebanese olive oil.

The interest in olive varieties with higher oil content (OC), improved fatty acid composition, mainly high monounsaturated fatty acids (MUFAs) and a high content of phenolic compounds, has increased due to its stability and health benefits (9,10). While OC is associated with oil quantity and olive growing profitability; the proportions of the different fatty acids and phenolic compounds are associated with oil quality. For example, a high percentage of MUFA, mainly oleic acid, is a primordial factor in determining the nutritional value of the oil as it reduces the risk of atherosclerosis (11) and protects against different kinds of cancers (12). In addition, fatty acid composition influences the stability of the oil through the contribution of polyunsaturated fatty acids (PUFAs) to oil rancidity (13). On the other hand, the amount of olive oil phenolic compounds, such as oleuropein derivatives is of primary importance when evaluating its quality, as these natural antioxidants improve oil resistance to oxidation and are responsible for its sharp bitter taste (14). The pharmacological interest of olive phenolic compounds is also well-known (15,16).
Several agro-industrial parameters may modify the OC, the fatty acid composition and phenolic content of virgin olive oil (VOO). Indeed, previous studies showed that OC, fatty acid and phenolic profiles are a built-in genetic factor. Values between 10 and 30% of oil were reported when evaluating the different accessions of the World Olive Germplasm Bank in Cordoba-Spain (17). Moreover, in the Germplasm Banks of Catalonia and Cordoba, Tous et al. (18) and Uceda et al. (19), respectively, showed that more than 70% of the variation in the fatty acids (except for linolenic acid) and several minor components, such as phenolic compounds, bitter index (K225), and oil stability, was due to genetic effects. In addition, evaluation of new varieties obtained by breeding programs showed that genotypic variance was the main contributor to the total variance of fatty acids (20,21) and of phenolic compounds (21,22). In addition, fruit ripening affects oil quantity and composition as early harvested olives, especially green olives, give lower oil quantity but higher content in MUFAs and antioxidants than late harvested ones (23)(24)(25). Additionally, OC and composition may vary according to the climatic conditions of the year, the irrigation regimes and the processing systems among others (23,(25)(26)(27).
In the past few years, authenticity and characterization of olive oils have been the object of numerous studies due to the importance of the protection of consumers. Various physicochemical determinations in association with chemometric analyses have been applied: fatty acids (28)(29)(30)(31), fatty acids and triacylglycerols (32,33), sterols (34), phenolic compounds (35), and aromas (36). Yet little attention has been given to the authentication and characterization of Lebanese olive oils. Chehade et al. (6,37)  in Lebanon are still not comprehensively assessed. In this manuscript, we report the characterization of monovarietal olive oils for 11 varieties cultivated in northern Lebanon for their OC, fatty acid and phenolic composition along fruit ripening with the perspective of evaluating and valorizing Lebanese olive oil. Expected findings will allow the development of a set of practical recommendations providing farmers with knowledge of the best varieties to be grown and the best time for olive harvesting, to further improve Lebanese production of olive oil in terms of quantity and quality.

Plant Material and Extraction of Olive Oil
Olive fruits were collected from 11 varieties growing in Abdeh station of the Lebanese Agricultural Research Institute (LARI) during the 2015 olive harvest season. This station, located at 18 m a.s.l., 34 • 31'0" N and 35 • 58'0" E, has been hosting an international olive collection of 72 olive varieties, initially collected from Lebanon, Arabic, and European countries, since the early seventies. The Abdeh region is characterized by a typical Mediterranean climate with a dry summer from June to September. The average annual precipitation sums 870 mm, the soil is a clay type soil with 25% sand, 15% silt, and 60% clay with 2% organic matter content. The olive trees in the Abdeh collection are rainfed.
The 11 genotypes considered included two clones of the local traditional variety, "Baladi 1" and "Baladi 2, " as previously differentiated by Chalak et al. (38), in addition to nine foreign varieties i.e., "Nabali" and "Jabaa" (from Palestine), "Kalamata" (derived from Greece), "Salonenque" and "Tanche" (from France), "Ascolana Tenera, " "Bella di Cerignola" and "Itrana" (originating from Italy), and "Sigoise" (derived from Algeria) (2,(39)(40)(41). For each variety, three fruit samples were collected from three trees, at different harvesting dates, based on fruit skin color (0= deep or dark green; 1 = yellowish-green; 2 = yellowish with reddish spots; 3 = reddish or light violet; and 4 = black), and the ripening index (RI). RI was determined on samples of 100 fruits according to the method described by Uceda and Hermoso (42) and modified by El Riachy et al. (43), following the formula: where, n 0 to n 4 , is the total number of olive fruits in each category. Then, RI was categorized to the following five categories: RI ≤ 0.5 (RI0), 0.5 ≥ RI ≤ 1.5 (RI1), 1.5 ≥ RI ≤ 2 (RI2), 2.5 ≥ RI ≤ 3.5 (RI3), and finally RI ≥ 3.5 (RI4). The olive oil was obtained using the Abencor system (MC2 Ingenierías y Sistemas, Sevilla, Spain). For each sample, ∼500 g of fruits were ground to a paste using a hammer mill, stirred in the thermobeater for 30 min at 28 ± 1 • C, and then centrifuged for 2 min to separate the oil, which was collected and left to decant in graduated cylinders. OIY was calculated using the following formula (44): where, V is the volume of oil obtained; and W is the weight of the processed olive paste. Finally, the oil extracted was collected and stored in glass vials in darkness at −20 • C until analysis.

Chemical Reagents
The chemical reagents used for VOO characterization were of GC or HPLC grade. Methanol, acetonitrile, o-phosphoric acid and the Folin-Ciocateu (F-C) reagent were provided from Sigma-Aldrich (Steinheim, Germany), n-hexane from Hipersolv Chromanorm (Pennsylvania, USA), Na 2 CO 3 anhydrous from Acros organics (New Jersey, USA) and KOH from Fisher Scientific (Loughborough, UK

Determination of the Oil Content (OC) in the Olive Paste Using Soxhlet Method
OC in the olive paste was determined using a Gerhardt Soxhlet instrument (Gerhardt Soxtherm SE−416, Germany) (45). From each sample, an aliquot of 50-60 g of the olive paste was dried at 105 • C to weight stability using a rapid moisture analyzer (Metler Toledo HS 153, Switzerland). One gram of the dry paste was placed on a filter paper which was folded and closed tightly using a cotton wire. Each prepared sample was placed in previously weighed soxhlet beakers containing three boiling chips. Petroleum benzene was used as a solvent. The total program length was 2 h and 20 min. When the process ended, the samples were thrown out, and the remaining solvent was eliminated by placing the beakers under a fume hood overnight, then, the beakers were dried in an oven at 105 • C for 1.5 h in a desiccator for 45 min before being weighed. OC was calculated from both dry matter (OCDM) and humid matter basis (OCHM) as follows: where, m 1 = beaker weight before extraction; m 2 = beaker weight after extraction and MC = moisture content.

Determination of Fatty Acid Profile
Fatty acid methyl esters (FAMEs) were prepared based on a cold transmethylation reference method (46). A sample of 0.1 g oil was shaken manually with 2 mL n-hexane for 2 s then 0.2 mL of methanolic solution (2 N) of potassium hydroxide was added. The sample was mixed with a vortex for 1 min (1,400 rpm), before resting for 5 min. A volume of 975 µL of the upper phase that contains the FAME was transferred to 1.5 mL vials with 25 µL of external standard (nonadecanoate methyl ester 1,000 ppm). The separation of FAMEs was carried out using a Shimadzu gas chromatograph (GC-2010 Plus) equipped with a flame ionization detector (FID). A fused silica capillary column (DB-wax; Agilent Technologies, Wilmington, DE; 30 m length × 0.25 mm i.d. and 0.25 µm of film thickness) was used. Nitrogen was used as a carrier gas with a flow rate of 1.69 mL/min. As for the chromatographic gradient, the initial oven temperature was kept at 165 • C for 15 min and then programmed to rise at 5 • C/min up to 200 • C, maintained for 2 min, and followed by a second gradient of 5 • C/min to a final temperature of 240 • C, which was held for 5 min. The injector and detector temperatures were 250 and 280 • C, respectively. Hydrogen and compressed air were used for the flame detector. Finally, the injection volume was 1 µL with a split ratio of 50. Identification of the different fatty acids was achieved by a comparison of their retention times with those of authentic standards. Results were expressed as percentages of total fatty acids. Finally, the calculated sums and ratios were; saturated fatty acid (SFA), MUFA, PUFA, MUFA/PUFA, MUFA/SFA, PUFA/SFA, C16:0/C18:2, C18:1/C18:2, and C18:2/C18:3.

Extraction of Phenolic Compounds
Prior analysis samples were left to thaw at room temperature. Meanwhile, an internal standard solution was prepared by dissolving 15 mg of syringic acid in 10 mL of 60:40 v/v methanolwater. Then, 1 mL from this solution was diluted in a 25 mL volumetric flask with 60:40 v/v methanol-water. The phenolic compounds in the olive oil were extracted using a modification of the procedure described by Montedoro et al. (47). An amount of 3 g of olive oil was shaken manually with 2 mL of n-hexane for 15 s. Then, a volume of 1.75 mL of 60:40 (v/v) methanol-water mixture was added together with 0.25 mL of the internal standard solution and shaken for 2 min to undergo the first extraction. For the second extraction, 2 mL of methanol/water (60/40) was added and shaken for 2 min. The extracts from both extractions were combined and placed in the dark at −20 • C for further determinations.

Spectrophotometric Estimation of Total Phenols
The total phenols (TP) content of the oil extracts was determined according to the Folin-Ciocalteu spectrophotometric method (48). Briefly, 20 µL of the sample (with prior 1:50 dilution with water) was, in this order, mixed with 1.58 mL of water, 0.3 mL of 20% (w/v) Na 2 CO 3 aqueous solution, and 0.1 mL of F-C reagent, and heated in an oven for 5 min at 50 • C. Then, the resulting solution was allowed to stand for 30 min. The reaction product was spectrophotometrically monitored at 765 nm by a jenway UV/Vis spectrophotometer (Staffordshire, ST15 OSA, UK). A nine level calibration curve was prepared using gallic acid as a commercial standard (r 2 = 0.998), and the results were expressed as mg gallic acid equivalent (GAE)/Kg of oil.

Chromatographic Analysis of Phenols by HPLC-DAD
The produced extract was also used to determine the following eight individual phenolic compounds: hydroxytyrosol, tyrosol, vanillic acid, p-coumaric acid, oleacein, oleocanthal, luteolin, and apigenin. The extracted phenolic fraction was injected in triplicate in a Shimadzu High Performance Liquid Chromatograph (HPLC) equipped with an automatic injector, a column oven and a diode array UV detector (DAD). Separation of individual phenols was achieved on a Microsorb-MV 100 C18 column (250 × 4.6 id mm, 5 µ particle size), maintained at 40 • C. The injection volume was 20 µL and the flow rate 1.0 mL/min. Mobile phases were 0.2% o-phosphoric acid in water (mobile phase A) and a mixture methanol-acetonitrile (50:50, v/v) (B). The initial concentrations were 96% of A and 4% of B and the gradient was changed as follows: the concentration of B was increased to 50% in 40 min, increased to 60% in 5 min, and to 100% in 15 min, and maintained for 10 min. Initial conditions were reached in 7 min. The identification of individual olive oil phenols was performed at 280 nm, on the basis of their maximum absorption and retention times compared to those of commercial standard compounds. Phenolic compounds quantification was achieved using syringic acid as an internal standard and 9 point calibration curves of authentic standards. Results were expressed as mg of the target analyte per Kg of oil.

Statistical Analysis
The total number of samples was 129 samples. Each sample was loaded in triplicate for individual phenol analysis, in duplicate for fat content and for fatty acid determination, and just once for the TP estimation. Multivariate analysis of variance (MANOVA) was used to assess the combined effect of variety, fruit ripening, and their interaction on the different sets of variables; and one-way analysis of variance (ANOVA) was used to test the significance of the effect of the studied factors and their interaction on each parameter studied. For all traits studied, trait means, coefficient of variation, and standard error were calculated; while means were compared using Tukey's HSD test at (p < 0.05). A principal component analysis (PCA) was also performed in order to determine the degree of contribution of each of the characters to the total variation and to highlight the effect of the studied factors on those traits (49). Data processing was performed using the Statistics (Analytical Software, Tallahasse, FL, USA) and Unscrambler (CAMO A/S, Trodheim, Norway) statistical packages.

RESULTS
The combined effects of the variety, the fruit ripening and their interaction were assessed for the sets of variables related to OC, fatty acids and phenolic composition. The results of the MANOVA ( Table 1) showed a highly significant effect of the interaction variety × fruit ripening (p < 0.001) on the three sets of variables. Similarly, each factor alone also showed a highly significant effect (p < 0.001) on these sets of variables. In this case, the variety and its associated error exhibited a more remarkable effect, expressed as partial η², while 64% of the change in the OC set can be accounted for by the variety and its associated error; 83% in the fatty acid composition set; and 69% in the phenolic composition set. Moreover, the two studied factors and their interaction showed a very high power (=1) sufficient to detect such effects.
OC traits and the majority of fatty acids and phenolic compounds appeared to be significantly affected by the fruit ripening ( Table 4). OC assessed as OCDM, OCHM, and OIY, showed a general increase along the ripening process; although, OCHM and OIY decreased with the increase of moisture Frontiers in Nutrition | www.frontiersin.org The rest of the phenols increased until RI = 2 in case of hydroxytyrosol, p-coumaric acid, luteolin, and apigenin and until RI = 3 in case of tyrosol and vanillic acid, before decreasing at RI = 4; although the difference was only significant (p < 0.005) in the case of vanillic acid and apigenin. The general evolution pattern described above along ripening was in some studied traits different between varieties (Figure 1). For instance, OIY increased along ripening in "Baladi 2, " "Itrana, " and "Nabali"; however, it increased until RI = 2 in case of "Ascolana, " "Baladi 1, " "Jabaa, " and Salonenque" and until RI = 3 in "Bella di Cerignola, " "Tanche, " "Kalamata, " and "Sigoise" before decreasing later on. Among fatty acids, C16:0 decreased along ripening in "Ascolana, " "Baladi 1, " "Baladi 2, " "Kalamata, " "Sigoise, " and "Tanche"; decreased until RI = 3 and then increased in "Bella di Cerignola" and "Nabali"; but increased until RI = 3 and then decreased slightly in "Itrana, " and "Salonenque." While TP generally decreased along ripening in almost all varieties; it increased in "Bella di Cerignola" and Jabaa. The evolution of TP was noticeable in "Baladi 1, " "Baladi 2, " "Itrana, " and "Kalamata" as it increased from RI = 1 to RI = 2 and then decreased to values lower than those recorded at RI = 1. As per the most concentrated individual phenols, oleacein decreased in "Sigoise"; decreased in "Bella di Cerignola" till RI = 2 and then increased; decreased in "Baladi 1" and "Nabali" till RI = 3 and then increased; and decreased till RI = 2, increased till RI = 3 and then decreased in "Kalamata." In contrast, oleacein content increased in "Baladi 2"; increased till RI=2 in "Ascolana Tenera, " "Itrana, " and "Jabaa" and then decreased; and increased till RI = 3 and then decreased in "Salonenque" and "Tanche." Oleocanthal showed the same evolution pattern as oleacein in the case of "Ascolana Tenera, " "Baladi 2, " "Bella di cerigonala, " "Itrana, " "Nabali, " and "Sigoise." However, it increased till RI = 2, decreased at RI = 3 and increased beyond in "Baladi 1" and "Jabaa"; increased till RI = 2 and then decreased in "Salonenque"; increased till RI = 3 and then decreased in "Kalamata"; and decreased till RI = 2, increased at RI = 3 and then decreased beyond in "Tanche." A Principal component analysis (PCA) was performed for fatty acid traits which are recognized as the most stable indicators for characterizing a given VOO (39). The results showed that the first two PCs explained 78.25% of the total variance (Figure 2). PC1 accounted for 58.49% of the total variance with a high positive correlation with PUFA and C18:2; and a high negative correlation with MUFA and C18:1. However, PC2 accounted for 19.76%, with a high positive correlation with SFA and C16:0. In addition, in the correlation circle, a negative correlation was observed between C18:1 and MUFA from one side and C18:2 and PUFA from the other side (data not shown).
The PCA was able to differentiate more between harvesting times rather than between varieties, as it can be seen that in the individual score plot (Figure 2) the samples from harvest 1, with higher content of C16:0 and C18:1 in comparison with other harvests, are all grouped in the negative side of the PC1 and in the positive side of PC2 except the samples of "Jabaa" and "Salonenque" as these varieties are characterized by the lowest content of C18:1 (57.14 and 57.59%, respectively) and the highest content of C18:2 (19.47 and 17.94%, respectively) among all varieties studied.

DISCUSSION
In this study, the characterization of 11 olive varieties cultivated in Lebanon was reported at different ripening stages using 28 traits relevant to OC, fatty acid, and phenolic profile. Results showed the strong effects of variety, fruit ripening and their interaction on the majority of the studied traits. This strong genetic variability revealed here is in parallel with the results previously reported for advanced selections of the olive breeding program in Cordoba (Spain) (43,50). Moreover, these results confirm the physical and chemical modifications occurring during fruit ripening, affecting OC, fatty acid and phenolic composition and consequently oil quality and oxidative stability (24,51). The OC traits, mainly OCDM (31.42-48.24%), OCHM (14.18-27.86%), and OIY (7.14-19.44%), showed high variability among varieties. Similar wide ranges of variation were observed by Chehade et al. (37) while characterizing the oil content and composition of five main Lebanese olive cultivars ("Aayrouni, " "Abou Chawkeh, " "Baladi, " "Del, " and "Soury"). These authors reported an OCHM between 17.8 and 36.4% and an OCDM between 42.5 and 56.2% in olives at spotted stage of ripening, with comparable values recorded for "Baladi" (28 ± 1% and 42.5 ± 1.8%, respectively). OCDM values increased along ripening reaching a significantly higher value at RI = 4; however, OCHM and OIY reached the maximum values at RI = 3 before decreasing later on probably due to rainfall in this period. This pattern of evolution of OC traits along ripening was also similar to the one described in other olive varieties in Spain (23,24) and Tunisia (25,52).
Regarding fatty acid profile, all percentages of fatty acids obtained in the present study fit, more or less, with the requirements established by the IOC for VOO, except for C18:3 which slightly exceeded the limit of IOC (C18 : 3 < 1%) in case of "Ascolana Tenera" (1.02%), "Jabaa" (1.12%), "Salonenque" (1.10%), and "Sigoise" (1.08%). Similarly, previous work found that the amounts of some fatty acids could be outside the ranges listed by IOC such as in French (53) (37). Oleic acid is well-recognized to be the most important fatty acid in olive oil, associated with its high nutritional value and oxidative stability (61,62). An olive variety is considered TABLE 3 | OC, fatty acids and phenolic compounds traits variability among the studied varieties cultivated in Lebanon.
Concerning the evolution of the different fatty acids along ripening, the results of this study showed a general decrease of C16:0 and C18:1 together with a significant increase of C18:2; in concordance with several previous studies (24,64,68). According to Gutiérrez et al. (23), the C16:0 level fell during the ripening process, possibly as a result of a dilution effect. However, the increase in C18:2 content could be linked to the activity of the enzyme oleate desaturase, transforming oleic acid into linoleic acid (23,69). These changes will produce a decrease of the ratio C18:1/C18:2 together with an increase of the ratio PUFA/SFA which will have a detrimental effect on the oil oxidative stability reducing the shelf-life of olive oils obtained from late harvest.
Results of total and individual phenols also expressed a large variability among the varieties studied. These results are in agreement with previous works indicating that genetic variability is the main factor affecting the phenolic composition (22,70). The TP values obtained for "Baladi 1" and "Baladi 2" are higher than those reported by Chehade et al. (37) and El Riachy et al. (63) for the same variety (167-249 mg GAE/kg of oil and 194-236 mg GAE/kg of oil, respectively). It is worth pointing out here the low content of secoiridoids in the local varieties in comparison to the foreign ones, which highlights the need for breeding to improve the traditional local varieties as these phenolic compounds have a detrimental effect on olive oil quality and nutraceutical properties and are very useful as selection criteria in breeding programs. Additionally, several works indicate a significant decrease of TP content during ripening (43,54,71) similar to the pattern described in the present paper. This decrease was most likely correlated with the increased activity of hydrolytic enzymes observed during ripening (72). In parallel, oleacein and oleocanthal also decreased from 61.87 to 46.57 mg/Kg and from 136.28 to 81.13 mg/Kg, respectively. It is worth noting that, the other studied individual phenols recorded the FIGURE 1 | Evolution pattern along ripening of the different studied traits. MC, moisture content; OIY, oil industrial yield; TP, total phenols; RI, ripening index. , "Ascolana Tenera;" ×, "Baladi 1;" , "Baladi 2;" •, "Bella di Cerignola;" *, "Itrana;" •, "Jabaa;" +, "Kalamata;" △, "Nabali;" ♦", "Salonenque;" , "Sigoise;" , "Tanche". FIGURE 2 | Biplot of principal components 1 and 2 based on fatty acid profile components recorded for each variety at different ripening indexes. Asc, "Ascolana Tenera;" Bal1, "Baladi 1;" Bal2, "Baladi 2;" Bel, "Bella di Cerignola;" Itr, "Itrana;" Jab, "Jabaa;" Kal, "Kalamata;" Nab, "Nabali;" Sal, "Salonenque;" Sig, "Sigoise;" Tan, "Tanche." 1, 2, 3, and 4 are the different harvesting times. maximum values at RI between 2 and 3 before decreasing later on. Therefore, early harvesting (before black pigmentation, RI = 4) should be recommended to obtain higher content of phenolic compounds in addition to higher content of C18:1. These results are in partial agreement with other studies that reported an increase in TP content to a maximum level at the "spotted" and "purple" pigmentation, the content decreasing drastically as ripening progressed (24,52).
Finally, PCA analysis showed that PC1 and PC2 explained the highest total variance (78%) when using the set of 15 traits related to fatty acid composition alone. This PCA allowed the classification of the samples according to the ripening stage with samples at RI = 1, having higher percentages of C16:0, C18:1, and MUFA, clustering together in the negative side of PC1 and positive side of PC2 (Figure 2). These results are in line with those obtained by León et al. (73) who separated 18 varieties growing in the World Olive Germplasm Bank (WOGB) of IFAPA (Cordoba, Spain) into four groups according to their fatty acid composition and to the country of origin, where the percentages of C18:1, C18:2, and saturated fatty acids were the main contributors of variation.
In conclusion, the present study reports the characterization of 11 local and foreign varieties growing under Lebanon climatic conditions. In comparison to foreign varieties, the local "Baladi 1" showed outstanding OCHM; however, both "Baladi 1" and "Baladi 2" recorded similar fatty acid composition but low phenolic composition. Indeed, many of the studied foreign varieties presented several outstanding characteristics such as "Kalamata" in terms of higher OCHM; "Tanche" in terms of higher OIY; "Tanche" and "Sigoise" in terms of higher C18:1, higher MUFA, higher MUFA/PUFA ratio, higher C18:1/C18:2, higher MUFA/SFA, and higher C18:1/C18:2; "Ascolana Tenera" and "Itrana" in terms of higher total phenols content; "Itrana" in terms of oleocanthal; and "Salonenque" in terms of oleacein. Taking into consideration that the appropriate time for harvesting is when a good balance between oil quality and oil quantity is achieved, and based on the fact that OCHM and OIY increased significantly until RI = 2 while C18:1, MUFA, MUFA/PUFA, C16:0/C18:2, C18:1/C18:2, TP, oleacein, and oleocanthal decreased along ripening, we can conclude that the best time for harvesting is at RI = 2 (when the color is yellowish green with reddish spots). This is true for all studied varieties, except "Sigoise" that must be harvested later at RI = 3 for better oil extractability. These findings on monovarietal olive oils should be confirmed over years in order to evaluate and valorize the local olive germplasm compared to the foreign one, and for the characterization and authentication of Lebanese olive oil. This is also critical in identifying the main default of the local varieties in order to start their genetic improvement through breeding programs.

DATA AVAILABILITY
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.