Volatile Compounds Analysis of Tunisian Propolis and Its Antifungal Activity

In this study, four samples of Tunisian propolis were analyzed for their volatile compounds. Static Headspace technique coupled with gas chromatography-mass spectrometry (GC/MS) was used for the volatile profile characterization. Statistical investigation of the samples was made applying hierarchical clustering method, K-Means cluster and principal component analysis. Results showed that 47, 36, 30 and 24 different compounds were identified, respectively, in the samples from Zouarine, Zelligua, Bizerte and Beni Khalled. Propolis volatiles were dominated by monoterpene hydrocarbons. α-pinene was the major compound representing 81.14%, 82.67% and 90.74%, respectively, of the total propolis volatiles collected from Zelligua, Beni Khalled and Zouarine and only 45.22% of the sample from Bizerte which had a very different composition. The in vitro antifungal activity of the volatiles from all samples against Candida albicans was also assayed and reported.


Introduction
Propolis is a resinous adhesive natural substance, collected by honeybees like (Apis mellifera L.) from buds, leaves of trees and plants then mixed with pollen and enzymes secreted by bees [1].
Propolis is used as a purpose sealer to smooth out the hive internal walls and as a barrier against intruders [2].It is a valuable product from ancient times to nowadays.Propolis is mostly used, in folk medicine, in pharmacy and in food technology, as a remedy with high antioxidant, antimicrobial, antifungal, anti-inflammatory, antiviral and antitumor activities [3] [4].
Chemical composition of propolis is dependent on its botanical and geographical origin [5] [6], races of bees [7] [8] and collection season where it can be collected all year [9] [10].
Despite the fact that volatiles are found in low contents in propolis, their pleasant aroma and significant biological activities make them very valuable for the propolis characterization [12].Several studies showed that spread of volatile compounds identified in propolis is very wide [13] [14] [15] [16] [17]; the chemical composition depends on the specificity of the local flora at the collection region and thus on its geographic and climatic characteristics [2].
Indeed, in contrast to propolis of Europe, Algeria and Egypt, Tunisian propolis has a different botanical origin due to unique Tunisian flora which presents a generally known biodiversity with a high percentage of endemic plants [13] [18] [19] [20].
Studies regarding Tunisian propolis volatiles are very scare and limited to that of Martos et al. [21] investigation of flavonoids composition in Tunisian propolis and that of Kouidhi et al. [22] who showed that propolis ethanol extract possessed anti-cariogenic and anti-biofilm activities, and had potential protective effect against cancer cells proliferation.To the best of our knowledge, no study has been conducted concerning Tunisian propolis volatile compounds.
Candida albicans is the most well-known fungal species associated with the clinical manifestation of this pathology.It is the most prevalent and pathogenic species among all Candida infections which have become a serious health problem.As fungi grow increasingly resistant to the available drugs; new drugs must be researched and evaluated for their effectiveness in antifungal treatment [23] [24].This strain has been selected for the basis of its application purpose of further formulation study.
So, the present work focused on the determination of Tunisian propolis volatiles and a comparative study of the volatile composition of propolis collected from different regions was performed.Samples were analyzed using static Headspace coupled with GC/MS, applying statistical data mining techniques: principal component, hierarchical clustering and K-Means cluster analysis for classification, afterwards.The in vitro antifungal activity of the volatiles from all samples against human pathogen Candida albicans was also assayed.

Propolis Samples
In this study, four samples of Tunisian propolis were collected during 2014's spring and analyzed for their volatiles composition: two samples were collected from Kef region (200 g from Zouarine, 100 g from Zelligua), 100 g from Nabeul (Beni Khalled) and 50 g from Bizerte.
Sampling was performed based on diversity in all four regions; this difference is noticed especially with the botanical origins.
Eucalyptus, Rosemary, Thym, Citrus and Almond trees were found in both Zouarine and Zelligua sites.The Kef region is characterized by the dominance of the forests and the mountains which promoted trees diversity such as: Conifers, Birch, Ash, Yardstick, Elm, Chestnut, Poplar, Beech, Pine, Spruce, Fir, Plums and Cactus [25].
In Bizerte, bees are collecting propolis only from Orange trees and the dominant climate was steppe.
In Beni Khalled region, Rosemary, Thym, Cactus and small plants are found but there was no Eucalyptus.

Volatile Compounds Analysis by Headspace Coupled with GC-MS
Volatile compounds were analyzed using a Headspace (TELEDYNE TEKEMAR HT3TM) coupled with an Agilent GC-MS system (GC with 7890A, mass detector 5975C with Triple-Axis, insert XL MSD).
1 g of propolis samples were introduced in a 30 ml headspace vials incubated for 30 min at 80˚C in headspace oven then transferred in heated line at 115˚C to avoid condensation of volatiles which were injected in the GC inlet during 1 min with a static mode.A HP-5 ms column (5% phenylmethylsiloxane) was used (30 mm × 250 µm × 0.25 µm).The carrier gas was Helium (N60 = 99.99%) with the flow rate of 1 mL/min.
Each run was performed during 24 min.The temperature was programmed at 40˚C for 1 min, raised to 100˚C with a rate of 10˚C /min then raised to 200˚C with a rate of 20˚C/min and kept constant for 5 min.It was finally raised to 300˚C with a rate of 20˚C and kept constant for 2 min.Injection was realized within 250˚C inlet with a splitless mode.The auxiliary temperature was 250˚C and the mass spectrometer was operating in EI mode (70 eV).
Quadripole and source temperature were respectively fixed at 150˚C and 250˚C with a full scan mode from 40 m/z to 550 m/z.The identification of the propolis components was based on a long work with NIST02 mass spectra search library.
Quantification was based on the areas of total ion current (TIC) peaks.Relative area values (percentage of total volatiles) were used for characterization purposes.

Antifungal Activity: Propolis Effect on Candida albicans Growth
In order to study the effect of Propolis volatiles on Candida albicans growth, an initial suspension of the yeast was done with an optical density of 0.5 on Mac-Farland scale (corresponding to 10 8 CFU•ml −1 ) and was decimally diluted with 0.9% (w/v) NaCl to obtain an inoculum concentration of 10 5 CFU•ml −1 .The antifungal assay was performed using the double-dish chamber method.100 µl of overnight C. albicans culture were spread on one half of the Petri dishes containing Sabouraud dextrose agar medium.The second half containing a scale range of Propolis concentration 1 mg, 10 mg, 20 mg and 50 mg.
Plates were parafilmed and incubated for 24 h at 30˚C.The percentage of growth inhibition (GI) was calculated according to the following formula: ( )

Statistical Data Analysis
For metrological repeatability of the measurements, each propolis sample was analyzed 3 times using Headspace coupled with GC-MS.According to GC-MS results, data matrix of [12 × 87] was formed for further analysis.Measured data were processed using XLSTAT2015 software.After data preparation, the following statistical methods were applied: principal component analysis (PCA), for reduction of the space of the variables; hierarchical clustering analysis (HCA) using average linkage, in order to estimate the similarity between any pair of the clusters applying average distance between all pairs of objects in any two clusters calculation; K-Means cluster analysis (KMCA), which is an iterative technique for classification that minimizes the sum, over all clusters, of the within-cluster sums of point-to-cluster-centroid Euclidean distances [26].

Volatile Compounds Composition of Tunisian Propolis
Visual analysis of the propolis samples showed that their color varied from golden yellow (Zouarine), opaque green (Bizerte), glossy dark brown (Zelligua) to dark (Beni Khalled) but all samples have a strong pleasant aroma.Indeed, the propolis color and odor is dependent on its botanical source [27] [28].Bankova et al. [12] concluded that the main botanical source of European, non-tropical Asian and North American propolis is Populus spp.
Besides the main source, bees also use pine, linden, willow, cherry, apple and other trees [27].In tropical regions there are no poplars; therefore bees use other plant sources for bee glue.Propolis volatiles give to the bee glue its specific pleasant aroma.It is well known that bees (Apis mellifera) respond to odors in several behavioral contexts [29] [30].
In Tunisia, different botanical sources can be found such as: Eucalyptus, Rosemary, Thym, Citrus and Orange trees.It is also characterized by the dominance of the forests around, which promotes diversity of trees where bees can collect propolis such as: conifers, birch, ash, yardstick, elm, chestnut, poplar, beech, apple, pine, spruce, fir, plums and cactus [20] [25].Headspace-GC-MS analysis (Table 1) revealed the presence of 83 different compounds in all samples.Only 6 volatile compounds present with low proportions (0.53% -0.082%) were not identified.Propolis sample collected from Zouarine was characterized by the highest number of compounds identified (47), followed by Zelligua (36) and Bizerte (30), while the lowest one was observed in the sample collected from Beni Khalled (24).The low quantity of volatile compounds identified in the Beni Khalled propolis could be due to the presence of waxes or may be related to flora.
compounds which were: camphene, verbenene, α-terpinene, ∆.3 carene, transpinocarveol and limonene.These similarities may be due to the close floral origin of both sites especially the distance between them that not exceed 20 km.
Our results on the propolis from Bizerte are in accordance with those mentioned by Melliou et al. [14] who used hydrodistillation and GC-MS for volatile extraction and analysis of Greek propolis; they showed that α-pinene was the major constituent (45.8%).These results were lower than our findings for the three other samples.In addition, our findings are totally different from those found on tropical propolis in which α-pinene was present in low proportion up to 1.6% [32] [33].However, Nunes and Guerreiro [34] performed Headspace analysis and attributed higher amounts of α-pinene, β-pinene, sulcatone, carene, limonene, eucalyptol, α-ocimene, β-ocimene, acetophenone and nonanal to the green Brazilian propolis.Moreover, α-pinene composed more than a half of all volatiles in Brazilian and Uruguayan propolis and over 40% in Iranian, Kerman and Greek ones which is lower than our results in Tunisian propolis.In the other hand, β-pinene was found in Uruguayan propolis with a higher amount (27.44%) than α-pinene (22.96%), in Brazilian propolis (20.85%), in Estonian propolis (8.86%) and even in Greek propolis (2.2%); these results were higher than ours [35].Our results showed that Cedrol was identified in the Bizerte sample with an amount of 8.23%.This was in accordance with Greek propolis [14], higher than Ethiopian propolis (2.48%) and much lower than Brazilian one (33%).

4-Terpineol was found in all Tunisian regions in low amounts which is different
from Ethiopian propolis where it composed one of the predominant compounds (8.57%) [15].
In Tunisian propolis, dl-limonene was found in small amounts of 0.54%, 1.42% and 1.20%, respectively, in the samples of Zouarine, Zelligua and Beni Khalled.These results were different comparing to Croatian [13] and Brazilian propolis from Jaguari region [36], where high percentages of dl-limonene were detected, using hydrodistillation method.This compound is known with its anti-anxious and antidepressant benefits [35] [36] [37].β-methyl crotonaldehyde was identified as an important compound in Brazilian (10.1%) and also in Chinese propolis [28], using Headspace analysis, but was not found in our samples.
Moreover, Eucalyptol presented a high amount of 25.95% in Estonian propolis but wasn't found in our samples.Benzaldehyde was only identified in the Tunisian propolis from Zelligua with a low amount (0.360%).Greek propolis presented a close amount of this compound (0.1% -0.3%), whereas Estonian sample contained high proportion about 10.85%.Chinese propolis was totally different than samples from Europe and South America but characterized by high amounts of 3-methyl-3-buten-1-ol and 3-methyl-2-buten-1-ol, (40.33% and 11.57%, respectively), 4-penten-1-yl acetate (9.04%) and α-longipinene (9.41%) [28].That last compound was found only in Chinese and Greek propolis and was also reported in Brazilian green propolis in high amount (24.89%) [38].All these compounds were not identified in our study and this might be due to the botanical, geographical and floral differences.Comparing to different European samples, a lot of similarities were found with the volatile composition of propolis from Brazil, Uruguayan, China and Greece but the nearest one to that of Tunisian propolis is the volatile composition of the Greek propolis [5] [37].In fact, from 83 compounds, 37 are similar.

Chemometric Analysis
For the statistical analysis standardized, data matrix [12 × 87] was prepared as described above.
First, the PCA was applied to reduce variables in space and therefore the correlation between them.The number of principal components that will represent the data, by explaining the corresponding part of the total variance of the initial variables and will be used in further statistical analysis was chosen according to Kaiser's eigenvalue criterion [29].In order to cluster the samples, two techniques were used: Hierarchical clustering Analysis (HCA) and K-Means clustering.The calculated distances and the results of the HCA are presented as a dendrogram in Figure 1 and Figure 2.
In fact, Figure 1    Khalled and Zelligua have close similarities in their composition.Propolis sample from Zouarine is different but still have some similarities with both previous samples.Finally, this analysis proved that propolis sample from Bizerte form a unique cluster, far from the others which mean, that its volatile composition is so different.
Before classification using KMCA, it is required to specify the clusters number and the iterations number to the algorithm of this analysis.The first parameter was calculated using a simple rule of thumb ( ) where operator ⌈ ⌉ means rounding the result to the nearest larger integer number and n is the number of the initial variables, n = 87.The iterations number was not constant as it depends on more critical parameter (minimum Euclidean distances sum of particular classification result).
In order to find the global minimum, KMCA was programmed to run 500 times.The minimum value of the calculated minimums of the Euclidean distances sum of each run was chosen as the best classification result.The results of Journal of Biosciences and Medicines KMCA, when k = 7, are presented as data scatter plots in three-dimensional space of first three PCA components in Figure 3.
Dendrogram in Figure 2, scatter plot in Figure 3 and diagram in Figure 4 presented clustering results.It is clearly seen that Bizerte, Zouarine, Zelligua and Beni Khalled propolis samples fall into 3 separate clusters, meaning that their volatile composition was very different.Beni Khalled and Zelligua samples formed one cluster since their volatiles composition was very similar.The sample from Zouarine formed the second cluster which is not very far from this group.
Bizerte sample had a very different composition from all other samples and constitute the third cluster.This classification is confirmed by looking to data scatter plots of K-Means cluster analysis Figure 3 which showed that Zelligua, Beni Khalled and Zouarine samples were near to each other and only the sample of Bizerte was removed from the group.For better certitude, the 9 common compounds cited above were selected for the statistical analysis and another data matrix [12 × 9] was made.PCA for common compounds was applied and results are presented in Figure 5.
In this case, all non-common compounds with small contributions take values close to zero, and the highest common compounds having the biggest impact on the classification results are analyzed.The samples of Bizerte, Zouarine, Zelligua and Beni Khalled propolis again fall into three separate clusters.Zelligua and Beni Khalled samples are close to each other despite the difference of geographical positions of the collection sites.However, Zelligua and Zouarine samples regardless of their belonging to the same region, presented some differences in volatiles composition.

Antifungal Activity: Propolis Effect on Candida albicans Growth
The volatiles of all samples were also studied for their antifungal activity against a human pathogen fungus (Candida albicans).
Results in Table 2 and Figure 5 and Figure 6 showed that only volatiles from Beni Khalled and Zelligua were active against C. albicans.The chemometric analysis revealed that volatile composition of propolis from Beni Khalled and Zelligua were similar and have 18 common compounds, as proved in Figure 4 from the PCA analysis, this similarity may explain the antifungal activity of both propolis samples compared to the propolis samples from other regions in which no activity were observed even in higher concentrations.The observed antifungal activity is probably due to the high amounts of α-pinene in propolis from Zelligua and Beni Khalled, respectively (81.140% -82.673%) which is well known to possess similar antifungal and antimicrobial activities [39].This hypothesis cannot be confirmed since volatiles in propolis from the region of Zouarine contained higher amount of α-pinene reaching 90.739% and although it was not active against Candida.Those conclusions may give us the confirmation that only a synergy between compounds found in propolis from both regions can lead to the antifungal activity observed.
Further tests on the anti candida albicans activity of common compounds in active propolis samples may prove to us weather this antifungal activity is due only to one of those compounds, a synergy of all three or a synergy between all volatile compounds in propolis.

Conclusions
The present work provides first data about the Tunisian propolis volatiles and reveals its interesting character.

N 1 :
Number of Candida colonies/box.N C : Number of Candida colonies in the control box.

A
. Jihene et al.DOI: 10.4236/jbm.2018.66009123 Journal of Biosciences and Medicines shows hierarchical cluster analysis HCA of all propolis volatile compositions.Dendograms revealed that propolis samples from Beni DOI: 10.4236/jbm.2018.66009124 Journal of Biosciences and Medicines

Figure 2 .
Figure 2. Plots of K-Means cluster analysis.

Figure 3 .
Figure 3. PCA presentation of compounds in propolis.

Figure 4 .
Figure 4. PCA for common compounds in all propolis samples.

Figure 5 .Figure 6 .
Figure 5.Effect of VOCs of Beni Khalled propolis against C. albicans.(a) Effect of 10 mg of propolis against C. albicans growth; (b) Effect of 50 mg of propolis against C. albicans growth.

Table 1 .
Percentages of Tunisian propolis volatile compounds.
The two propolis samples (Zouarine and Zelligua) collected from the same region of Kef presented 6 other common DOI: 10.4236/jbm.2018.66009122 Journal of Biosciences and Medicines

Table 2 .
Growth inhibition of C. albicans by propolis VOCs.Values are the means of inhibitory rate SD, means based on 3 replicates.