Next Article in Journal
Impacts of Hydrothermal Treatments on the Morphology, Structural Characteristics, and In Vitro Digestibility of Water Caltrop Starch
Next Article in Special Issue
Rapid and Accurate Approach for Honeybee Pollen Analysis Using ED-XRF and FTIR Spectroscopy
Previous Article in Journal
Microwave-Assisted Defibrillation of Microalgae
Previous Article in Special Issue
Assessment of the Botanical Origin of Polish Honeys Based on Physicochemical Properties and Bioactive Components with Chemometric Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Volatile Profile of Portuguese Monofloral Honeys: Significance in Botanical Origin Determination

by
Alexandra M. Machado
1,
Marília Antunes
2,
Maria Graça Miguel
3,
Miguel Vilas-Boas
4 and
Ana Cristina Figueiredo
1,*
1
Centro de Estudos do Ambiente e do Mar (CESAM Lisboa), Centro de Biotecnologia Vegetal (CBV), Faculdade de Ciências da Universidade de Lisboa, DBV, C2, Piso 1, Campo Grande, 1749-016 Lisboa, Portugal
2
Centro de Estatística e Aplicações (CEAUL), Departamento de Estatística e Investigação Operacional, Faculdade de Ciências da Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal
3
Faculdade de Ciências e Tecnologia, Mediterranean Institute for Agriculture, Environment and Development, Campus de Gambelas, Universidade do Algarve, 8005-139 Faro, Portugal
4
CIMO, Centro de Investigação de Montanha, Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(16), 4970; https://doi.org/10.3390/molecules26164970
Submission received: 29 July 2021 / Revised: 11 August 2021 / Accepted: 12 August 2021 / Published: 17 August 2021

Abstract

:
The volatile profiles of 51 samples from 12 monofloral-labelled Portuguese honey types were assessed. Honeys of bell heather, carob tree, chestnut, eucalyptus, incense, lavender, orange, rape, raspberry, rosemary, sunflower and strawberry tree were collected from several regions from mainland Portugal and from the Azores Islands. When available, the corresponding flower volatiles were comparatively evaluated. Honey volatiles were isolated using two different extraction methods, solid-phase microextraction (SPME) and hydrodistillation (HD), with HD proving to be more effective in the number of volatiles extracted. Agglomerative cluster analysis of honey HD volatiles evidenced two main clusters, one of which had nine sub-clusters. Components grouped by biosynthetic pathway defined alkanes and fatty acids as dominant, namely n-nonadecane, n-heneicosane, n-tricosane and n-pentacosane and palmitic, linoleic and oleic acids. Oxygen-containing monoterpenes, such as cis- and trans-linalool oxide (furanoid), hotrienol and the apocarotenoid α-isophorone, were also present in lower amounts. Aromatic amino acid derivatives were also identified, namely benzene acetaldehyde and 3,4,5-trimethylphenol. Fully grown classification tree analysis allowed the identification of the most relevant volatiles for discriminating the different honey types. Twelve volatile compounds were enough to fully discriminate eleven honey types (92%) according to the botanical origin.

1. Introduction

Honey is a food product that has been consumed by mankind for since at least 6000 years [1]. Carbohydrates are the main components, comprising about 95% of honey’s dry weight, including the dominant monosaccharides fructose and glucose, in addition to about 25 different oligosaccharides [2]. Other constituents, such as water (10–20% w/w), the second-largest component of honey, and other minor compounds are present, such as proteins, free amino acids, enzymes, organic acids, lipids, volatiles, vitamins, minerals and phenolic compounds [3,4]. Honey consumption has increased in the last few years mainly due to a growth in the world population as well as due to a preference for natural foods of a rising number of consumers [5].
According to both Codex Alimentarius and the European Honey Directive [6,7], honey must not have any added food ingredients, including food additives or other substances, apart from other forms of honey. The authenticity of this food product comprises two main aspects, the origin of the honey, which includes botanical source and geographical origin, and the production method related to the harvesting of honey and its processing [8]. Honey’s botanical source strongly affects its organoleptic properties [9,10], monofloral honey being more in demand compared to multifloral, due to its biological activities as well as the characteristic aroma and taste [11]. Honey is recognised as a high-quality food product, its falsification being a major problem around the world since it is one of the most adulterated products, after olive oil and milk [5,10,12,13,14]. The increase in the number of consumers concerned about honey’s authenticity and environmental sustainability [15] has increased the demand for reliable analytical methods to establish criteria to guarantee honey authenticity [10].
Honey’s aroma depends on the volatile compounds present, mostly in the nectar, which help to discriminate honey types from different botanical and geographical origins [16,17,18,19,20]. The mainland Portugal and Azores and Madeira Islands are characterised by a rich and varied honey flora, contributing to the production of a great diversity of local monofloral honeys. Most representative monofloral honeys include those obtained from the mainland from bell heather (Erica spp.), eucalyptus (Eucalyptus spp.) and lavender (Lavandula luisieri, L. pedunculata, L. stoechas). Other monofloral honeys are also important, such as those from chestnut (Castanea sativa), incense (Pittosporum undulatum) from Azores Islands, orange (Citrus spp.), pennyroyal (Mentha pulegium), rosemary (Rosmarinus officinalis), strawberry tree (Arbutus unedo), sunflower (Helianthus annuus) and viper’s bugloss (Echium spp.) [21]. In 2019, honey production in Portugal was about 10,104 t, including 127 t and 46 t from Azores and Madeira Islands, respectively [22]. Only a small number of studies concern the characterisation of volatiles of Portuguese honeys, namely bell heather, chestnut, eucalyptus and multifloral honeys [23,24,25,26].
Whereas gas chromatography–mass spectrometry (GC-MS) is the technique mostly used to identify volatiles, different volatile extraction techniques have been used, such as hydrodistillation (HD), solid-phase microextraction (SPME), liquid–liquid extraction, simultaneous distillation extraction, simultaneous steam distillation–dichloromethane extraction (using Likens–Nickerson apparatus), ultrasonic solvent extraction or the purge-and-trap technique [19,27,28].
With a view to integrate knowledge on the volatiles of Portuguese monofloral honeys and determine whether they can be used as a tool for honey authenticity determination, this work aimed at defining (a) the volatile composition of Portuguese monofloral honeys from different botanical origins by two different methods: SPME and HD; (b) whether there is a correspondence between the flower’s volatiles and the monofloral honey type’s volatiles; (c) the degree of correlation between the honey samples; and (d) the discriminatory potential of specific volatile compounds regarding the botanical origins of different monofloral honey types.

2. Results and Discussion

2.1. Honey and Flower Volatile Profiles

The volatiles from 12 monofloral-labelled Portuguese honey types were assessed in a total of 51 samples (Figure 1, Table 1, Supplementary Material Figure S1). Honey volatiles were isolated by two different methods, solid-phase microextraction (SPME) and hydrodistillation (HD) (Figure 1). To evaluate whether the presence of some honey volatiles is due to the volatiles from flowers visited by the bees, the main components obtained by HD from some selected flower species were also determined (Table 1, Supplementary Material Figure S2).

2.1.1. Volatiles Collected by Solid-Phase Microextraction (SPME)

The relative amounts of all the identified compounds obtained by SPME are listed in SM Table S1 according to the minimum and maximum percentage range of components found for honey samples. In total, 61 volatiles were identified in the 51 analysed honey samples. The identified compounds were dominated by benzoic acid derivatives and oxygen-containing monoterpenes, ranging from not detected (nd) to 59% and from nd to 47%, respectively. Alkanes and fatty acid esters were also identified in most of the samples. Apocarotenoids were found in great amounts in a specific group comprising strawberry tree honeys, ranging from 60 to 88%. Since with SPME, a lower number of volatiles was detected and their vast majority was also detected with HD isolation procedure, HD data were selected for further processing.

2.1.2. Volatiles Collected by Hydrodistillation (HD)

The relative amounts of all the identified volatile compounds isolated by HD from the 51 monofloral-labelled honey samples are listed in SM Table S2 in order of their elution on the DB-1 column. One hundred ninety-two (192) volatile compounds were identified in the analysed honey samples. The main identified compounds included hydrocarbons, mostly alkanes (11 to 77%). Alkenes (nd–24%), methyl-branched alkanes (nd–2%), long-chain fatty acids (1–45%), their esters and aliphatic fatty alcohols (1–16%) were also found in relatively high percentages. In lower amounts, oxygen-containing monoterpenes and sesquiterpenes were identified in all honey samples, ranging from <0.05 to 27% and from <0.05 to 4%, respectively. Aromatic amino acid derivatives were also found in all samples (0.1–25%), as well as green leaf volatiles (GLVs) in lower amounts, although GLVs were not identified in the strawberry tree honey samples. Apocarotenoids were the dominant group found in this honey type (5–44%).
The analysis of selected flower volatiles confirmed the presence of common compounds identified both in the flowers as well as in the corresponding monofloral honeys (Table 1).
The volatile profile obtained by the hydrodistillation of each of the honey types analysed is, in this section, compared with data on the honey’s volatile marker compounds from literature.
Carob tree honey samples showed predominance of cis- and trans-linalool oxides (furanoid and pyranoid), (0.1–15% and <0.05–5% correspond to the furanoid form, -cis and -trans, respectively, and (nd–0.3%, nd–0.6%) correspond to the pyranoid form, -cis and -trans, respectively, also identified in the volatile profile of carob tree flowers. These compounds were also identified in the flower volatiles from C. siliqua collected in the south of Portugal [29]. cis- and trans-Linalool oxides (furanoid) and hotrienol have also been identified as the main marker volatile compounds in carob tree honey samples from Morocco and Balearic Islands [30,31]. In the current work, cis- and trans-linalool oxides (pyranoid) and trans-β-damascenone (nd–0.5%) were also considered as possible marker compounds for this honey type.
2′-Aminoacetophenone (<0.05–0.7%) as well as acetophenone (0.1–0.4%) were the characteristic volatiles identified in chestnut honey samples, with the second being also detected in the corresponding flowers. Chestnut honey aroma has been characterised by the presence of 2′-aminoacetophenone, acetophenone, benzaldehyde, phenyl ethyl alcohol, linalool, n-nonanol and n-nonanal, as reported in Croatian, Italian and Spanish chestnut honeys [32,33,34]. A previous work with Portuguese chestnut honey [25] showed benzaldehyde, n-nonanol, n-nonanal and linalool in common with the current study. The results obtained herein indicate acetophenone and 2′-aminoacetophenone as volatile markers for this honey type, with only acetophenone identified in the flowers.
Eucalyptus honey samples evidenced 1-nonanol (<0.05–0.6%), α-eudesmol (0.1–0.5%), β-eudesmol (0.1–0.4%), n-nonanal (0.1–0.4%), aromadendrene (nd–0.4%), 2-undecanol (nd–0.3%) and p-cymen-8-ol (nd–0.2%) as characteristic compounds, with α- and β-eudesmol as well as aromadendrene also present in the respective flowers. 4-Keto-isophorone, 2-hydroxyisophorone, nerol oxide, p-cymene and p-cymen-8-ol contribute to the recognisable eucalyptus honey aroma [35,36,37], in agreement with the current results. The correspondence found with the sesquiterpene aromadendrene between eucalyptus flowers and the honey samples could indicate this compound as a possible volatile marker for Portuguese eucalyptus honey.
The distinguishing volatiles found for bell heather honey were benzene acetaldehyde (2–24%), cis-linalool oxide (furanoid, 2–12%), hotrienol (2–10%), trans-linalool oxide (furanoid, 1.0–5%), 1,2-dihydro-1,1,6-trimethyl-naphtalene (<0.05–3%), benzaldehyde (0.1–2%) and benzylalcohol (<0.05–0.9%). 1,2-Dihydro-1,1,6-trimethyl-naphtalene was likewise identified in bell heather flowers, with benzaldehyde and benzene acetaldehyde. Samples H1 and H3 were characterised by having benzene acetaldehyde (17–24%) and hotrienol (3–7%) as dominant compounds, while samples H4 and H6 exhibited cis-linalool oxide (furanoid, 9–12%), hotrienol (3–10%) and trans-linalool oxide (furanoid, 4–5%) as the main compounds. Samples H2 and H5 showed lower percentages of these compounds. Benzene acetaldehyde, cis-linalool oxide and hotrienol have also been mentioned as characteristic compounds in studies with Portuguese and Spanish bell heather honeys [24,25]. Although benzene acetaldehyde could be a Strecker degradation product of phenylalanine [28] or used as a bee repellent to facilitate honey collection [38], this compound was also identified in other studies with different honey types [39]. Phenylalanine is one main amino acid present in bell heather honey [40], and so its content in the correspondent plants should be determined to obtain a relationship between it and benzene acetaldehyde. Compounds shared with other honey types, such as isovaleric acid, 2-methylbutyric acid, benzyl alcohol and phenyl ethyl alcohol were also identified. Other studies have also reported the presence of apocarotenoids, such as α-isophorone [41] and decanoic acid and p-anisaldehyde, also identified in the current work, indicating a botanical origin within the Ericaceae family [42].
Benzyl salicylate was the distinguishing compound of incense honey (nd–0.3%), also detected in the corresponding flowers. This could be indicative of a marker compound for this honey type. However, other volatiles were also found, namely benzene acetaldehyde as well as cis-linalool oxide (furanoid) in the same amounts (1–3%), trans-linalool oxide (furanoid, 0.6–2%), benzaldehyde (<0.05–0.4%), α-eudesmol (nd–0.4%), α-terpineol (<0.05–0.3%) and limonene (nd–0.3%). Incense honey samples were characterised by the absence of hotrienol and by the presence of oleic acid as the main component (17–39%). There are few studies regarding this honey type, this being the first time that its volatile compounds were identified.
Lavender honey samples showed n-nonanal (nd–0.8%) and n-decanal (nd–0.4%) as differentiating compounds, although they were not found in the analysed flowers. Benzaldehyde, hexanoic acid, benzene acetaldehyde and phenyl ethyl alcohol were also detected, as mentioned in another work with Portuguese lavender honey specifically obtained from Lavandula stoechas [43]. The volatiles’ composition suggests that lavender honey is distinctive relative to the other analysed honey types; however, no specific markers were found.
Several distinctive volatile compounds were identified in orange honey, such as methyl anthranilate (nd–2%); lilac aldehydes A (nd–0.8%), B (nd–1%), C (nd–0.7%), D (nd–0.2%) and E (nd–0.05%); indole (nd–0.3%) and cis-myrcenol (nd–0.2%). With the exception of lilac aldehydes D and E, the remaining compounds were also identified in orange flowers, which could confirm their botanical origin. Lilac aldehydes A, B and C and methyl anthranilate have also been mentioned by other authors as markers for orange honeys from Greece and Spain [44,45,46]. Additionally, indole was identified in this study in both flowers and honey, which may indicate that this compound has a possible marker for Citrus spp. honey type, and has also been identified by other authors in the nectar of orange blossoms [39,47].
Rape honey samples only evidenced dimethyl trisulphide as a distinguishing compound in trace amounts (<0.05%), which is in agreement with the results reported by other works [48,49]. These samples were relatively poor in the number of volatiles, and they were common to several honeys, such as benzene acetaldehyde and hotrienol, also reported for rape honeys obtained from the Czech Republic [50].
Raspberry and rosemary honey samples did not evidence specific volatile compounds. There are only a few studies on raspberry honey volatiles and report compounds in common with the present work, for instance, cis- and trans-linalool oxide (furanoid); hotrienol; and lilac aldehydes A, B and C identified in Slovakian honey [51].
The alcohols 3-methyl-1-butanol and 3-methyl-3-buten-1-ol, the apocarotenoid 4-keto-isophorone and lilac aldehyde isomers have been previously identified in rosemary honeys [35,52]. Whereas these compounds were not detected in the present study, others were, which are considered ubiquitous honey constituents, such as benzaldehyde, benzene acetaldehyde and hotrienol. Probably because only one sample was studied in this work, no marker compounds were found for this honey type.
The main volatiles of sunflower honey’s were α-pinene (<0.05–0.2%), also detected in large amounts in the flowers, as well as in other honey types and flower volatiles. Benzene acetaldehyde (0.1–1%) and β-copaene (0.1–0.2%) were also detected in lower amounts. α-Pinene was the only compound identified in this honey, in common with the honeys of other countries, namely the Czech Republic, Romania and Spain [53]. However, β-copaene was only identified in this honey type, suggesting it as a possible marker for Portuguese sunflower honey.
Strawberry tree honey samples evidenced the presence of the following volatiles: α-isophorone (5–39%), 3,4,5-trimethylphenol (2–13%), 2,3,5-trimethylphenol (0.4–4%), 4-keto-isophorone (0.4–3%), 2-hydroxyisophorone (0.2–3%), veratrole (nd–2%), 1,2-dihydro-1,1,6-trimethyl-naphtalene (<0.05–1%), edulan (<0.05–0.7%), 5-methyl-3-hexen-2-one (<0.05–0.4%) and 2,6-dimethylpyrazine (nd–0.1%). In common between the flowers and the honey were α-isophorone, 2,3,5-trimethylphenol, 3,4,5-trimethylphenol and 4-keto-isophorone, suggesting them as possible markers for strawberry tree honey. The apocarotenoids identified, α-isophorone and 4-oxoisophorone, have also been reported by other authors as floral markers for this honey type produced in Greece and Sardinia [54,55,56]. Furthermore, in this study, other compounds were also identified, such as the aromatic amino acid derivatives veratrole and 2,3,5-trimethylphenol, as well as 3,4,5-trimethylphenol in higher amounts. Despite not being identified in the present work, benzene derivatives (dimethyl, trimethyl, tetramethyl) were identified by Karabagias et al. [54] in Greek strawberry tree honeys.
Notwithstanding their contribution to the honey volatile profile, some identified compounds cannot be considered as reliable botanical origin honey markers, not only due to their presence in the comb as well as in bees and cuticular wax, but also because of their general occurrence in different honey types. Among these are alkanes (n-tricosane, n-pentacosane, n-heptacosane), alkenes ((Z)-9-tricosene and (Z)-9-pentacosene)), methyl-branched hydrocarbons, long-chain fatty acids (oleic acid) and their esters (ethyl oleate), as well as aliphatic fatty alcohols (1-hexadecanol, oleyl alcohol, 1-octadecanol). The results reported herein are in line with those of different authors [57,58,59,60].
Among the compounds common to different honey types that cannot be considered as specific volatile markers are benzaldehyde, benzene acetaldehyde, benzyl alcohol, isovaleric acid, 2-methylbutyric acid and linalool derivatives, such as cis- and trans-linalool oxides (furan isomers). Compounds that could be generated in honey during heating were also identified, namely furan and pyran derivatives [28]. However, the amount of these compounds, which could be identified as artefacts, was quite small (≤2.3%), since the volatile extraction was performed for just 1 h to diminish the occurrence of hydrolysis reactions and the formation of sugar degradation products [61].
Despite several studies describing the identification of specific volatile compound markers in some honey types [30,35,44,62,63], only a few compounds seem really specific, with variable concentrations according to the different plant sources.

2.2. Statistical Analysis

Since HD was more effective in the number of isolated volatiles, statistical analysis was performed with the results obtained by this extraction method using all 51 honey samples.

2.2.1. Cluster Analysis

Hierarchical clustering was used to find whether the samples of each monofloral honey type were chemically correlated and whether there was correlation between honey types based on the composition of the volatiles (Figure 2, Table 1 and Table 2). Cluster analysis evidenced a dendrogram with two main clusters, one of which had nine sub-clusters. Despite the main clusters, I and II, showing very low correlation (Scorr < 0.12), in general, honeys from the same botanical origin were assigned to the same cluster and sub-cluster, highlighting their volatiles’ similarity.
Cluster I included nine moderately correlated sub-clusters, Ia–Ii (Scorr < 0.44), in which were included all honey types except three samples of strawberry tree honey (cluster II). Cluster I was characterised by the dominance of straight-chain hydrocarbons n-C23 (6–29%), n-C9 (1–26%), n-C21 (1–26%), n-C5 (0.3–16%) and fatty acids, as well as their derivatives oleic acid (nd–39%), ethyl oleate (nd–11%) and linoleic acid (nd–9%). Oxygen-containing monoterpenes such as cis- and trans-linalool oxide (furanoid) were also present, ranging from nd to 15% and nd to 5%, respectively.
Four carob tree honey samples evidenced high correlation (Scorr < 0.82) in sub-cluster Ia (Ct1, Ct2, Ct3, Ct5), while Ct4, positioned in sub-cluster Ig, showed moderate correlation with the former ones (Scorr < 0.58). Both chestnut (Scorr < 0.86) and eucalyptus honey samples were highly correlated in sub-cluster Ia (Scorr < 0.82).
Bell heather honey samples evidenced different main volatiles according to geographical region, highlighting its influence. The six bell heather honey samples formed three groups with two samples in each. Specifically, sub-cluster Ia showed very highly correlated (Scorr < 0.92) H2 and H5, both obtained from Trás-os-Montes e Alto Douro; sub-cluster Ig had highly correlated (Scorr < 0.80) H4 and H6 from the Minho region; and sub-cluster Ii had highly correlated (Scorr < 0.90) H1 and H3 from Beira Baixa. Benzene acetaldehyde was dominant in the Ii cluster.
Sub-cluster Ih comprised the four incense honey samples I1-I4, highly correlated (Scorr < 0.90), evidencing oleic acid (17–39%) and linoleic acid (nd–9%) as the main components.
Lavender honey samples were all highly correlated (Scorr < 0.72), with six of them in sub-cluster Ia (L1, L2, L4–L6, L8), while samples L3 and L7 were isolated in sub-clusters Ib and Ic, respectively.
Orange honey samples were positioned in three separate sub-clusters. Ia had the highly correlated (Scorr < 0.86) O2 to O6 and O9 samples; sub-cluster If had very highly correlated (Scorr < 0.98) O7 and O8; and O1 was placed alone in sub-cluster Ih.
Rape honey samples were positioned in different sub-clusters, R1 in Ie and R2 in Ia, showing moderate correlation between them (Scorr < 0.68).
Raspberry honey samples also exhibited moderate correlation (Scorr < 0.52), with Rb1 placed in sub-cluster Ih and Rb2 in Ia.
The only rosemary honey sample (Ro1) was found in sub-cluster Ia.
Sf1 and Sf3 sunflower honey samples were very highly correlated (Scorr < 0.94) in sub-cluster Ia, and sample Sf2 was placed in sub-cluster Id, moderately correlated with the other two (Scorr < 0.70). Sunflower honey samples produced in Alentejo remained together in the dendrogram, while the other, obtained from Algarve, was positioned in a different sub-cluster, indicating an influence of the geographical origin on the volatile composition. The compounds β-myrcene, benzene acetaldehyde, 2,3,5-trimethylphenol, palmitic acid and oleyl alcohol could be geographical region dependent for this honey type.
Cluster II, with three samples of strawberry tree honey, showed α-isophorone (36–39%) and 3,4,5-trimethylphenol (4–13%) as the main compounds. The three samples included (St1–St3) showed very high correlation (Scorr < 0.96), while the other strawberry tree honey sample, St4, placed in sub-cluster Ia, showed very low correlation with the former ones (Scorr < 0.12). The differences observed could be explained due to geographical reasons. Three of them were obtained from Alentejo and Algarve (St1–St3), while St4, with a lower amount of apocarotenoids and placed in a different cluster, was also produced in Algarve but in a region with peculiar geographical and climatic conditions, in Aljezur, surrounded by the sea and the hills, reflecting this dual influence also on honey volatile composition.

2.2.2. Classification Tree

A fully grown classification tree was built with the volatiles obtained by each honey’s hydrodistillation (Figure 3). It was found that 12 compounds were enough to fully discriminate 11 out of the 12 honey types, detailed in the bottom nodes, according to the botanical origin. The discriminant compounds were benzaldehyde, benzene acetaldehyde, β-copaene, cis-linalool oxide (furanoid), n-decane, ethyl hexadecanoate, α-eudesmol, 2-furfural, heptacosene, hotrienol, n-tricosane and UI 8. The square at each node shows the honey types present in the node, the number of honey samples included and the percentage of the data encompassed in that node. At the root (top node), the number of labelled honey samples for each honey type is shown in two rows. The first row (from left to right) exhibits two of the chestnut honeys, five of the carob tree honeys, five of the eucalyptus honeys, six of the bell heather honeys, four of the incense honeys and eight of the lavender honeys. From left to right, the second row shows the remaining samples according to honey types: nine of orange, two of rape, two of raspberry, one of rosemary, three of sunflower and four of strawberry tree. The root node contains 100% of the data. The first split (evaluation of the condition UI 8 < 0.25) puts 82% of the data in a group where the lavender honey is the most frequent honey type and 18% in a group where the orange honey is the most frequent. After this split, the classification tree continued the discrimination of the honey types, according to volatile compounds, along the internal nodes, and at the bottom of the tree, most of the samples were included in the node labelled with the correspondent honey type.
The classification tree suggests that the 12 identified volatile compounds are potentially good discriminators of honey types, since they were enough to fully discriminate 11 different monofloral honeys under analysis. These are important results because they have allowed pinpointing, from a large number of volatiles, some compounds able to discriminate honey types. Hotrienol, identified as one of the dominant compounds in heather honey, was used to discriminate this honey type, while cis-linalool oxide (furanoid) was used for carob tree and α-eudesmol for eucalyptus honey samples. β-Copaene, only identified in sunflower honey samples, was also useful to discriminate this honey type.
In the present global market, food product authenticity allied to quality assurance is a mandatory requirement. Honey, as a food product, linked with the growing use in the pharmaceutical industry, is widely consumed all over the world. The guarantee of honey authenticity can be achieved through a joint effort, with the integration of several analysis techniques. Analysis of honey’s volatiles is one additional technique that can be applied to the study of monofloral honeys to distinguish between those from different botanical sources as well as geographical origins. This work showed that is possible to differentiate monofloral honeys according to the identified volatile compounds and, in some cases, geographical locations.
The results of this study increased our knowledge about the dominant compounds found in different honey types produced in several regions in mainland Portugal and the Azores Islands. The detection of some specific compounds in the honeys, also identified in the respective flowers, allowed understanding the origin of marker compounds, characteristic of the analysed monofloral honeys. Further work on honey volatiles with a larger number of samples per honey type would help build a robust classification tool, allowing discrimination between honeys from different botanical sources. This knowledge could be useful for other studies on the nutritional and therapeutic potential of Portuguese honey.

3. Materials and Methods

3.1. Honey Sampling

Fifty-one honey samples labelled as monofloral according to botanical origin were collected from several regions of mainland Portugal, specifically Minho, Trás-os-Montes e Alto Douro, Beira Baixa, Estremadura, Alto Alentejo, Baixo Alentejo, Alentejo Litoral and Algarve, and from Azores Islands, São Miguel and Pico, as detailed in Table 1. Most of the samples were obtained from producers and the remaining from specialised shops, between 2015 and 2018. Honey samples were labelled with 12 different botanical origins when acquired. The honey samples were stored in a cool, dry place until further assessment. All the samples were subjected to pollen analysis to confirm their botanical origin, as previously described [64]. Flowers corresponding to the different honey types were collected from different geographical locations in Portugal, as detailed in Table 1.

3.2. Honey and Flower Volatile Sampling

The volatiles from 51 monofloral honeys were isolated by two different methodologies, (a) solid-phase microextraction (SPME) and (b) hydrodistillation, to compare the efficiency of volatile extraction. Flower volatiles were isolated by hydrodistillation.

3.2.1. Sampling of Honey Volatiles by SPME

Headspace volatiles were collected by SPME using a 100 μM polydimethylsiloxane (PDMS)-coated fibre (Supelco, Bellefonte, PA, USA) inserted into a manually operated SPME holder. Each SPME fibre was thermally conditioned for up to 20 min at 250 °C according to the manufacturer’s recommendations before use. Blank assays of the fibres were performed regularly.
For volatiles’ collection and after lid removal, each of the monofloral-labelled honey flask was inserted into a glass desiccator (ø 20 cm) (Figure 1, top row) and left to stand for 1 h at room temperature for atmosphere homogenisation. The conditioned fibres were exposed in the desiccator for 1 h. Three SPME fibres were used per desiccator, with two being analysed by gas chromatography (GC) with flame ionisation detection for quantification of volatiles and the other being analysed by gas chromatography–mass spectrometry (GC-MS) for chemical composition analysis. Preliminary experiments were run to ascertain the optimal period and temperature of fibre exposure.

3.2.2. Sampling of Honey and Flower Volatiles by Hydrodistillation (HD)

Volatiles were obtained by hydrodistillation using a Clevenger-type apparatus according to the European Pharmacopoeia (Council of Europe, 2010). In the case of honeys, hydrodistillation was run for 1 h (Figure 1, bottom row) and from 1 to 3 h for flowers, depending on the amount of the plant material available (Figure S1). Approximately 50 g of each honey, used in SPME analysis, was added to 100 mL of distilled water before hydrodistillation. The extracted organic compounds were recovered with in-lab distilled pentane from the supernatant layer on the distilled water (hydrolate) that remained after the hydrodistillation was complete. The pentane solution was then concentrated to a volume of ~100 µL at room temperature under a stream of nitrogen. Honey and flower volatiles were stored at −20 °C until analysis.

3.3. Analysis and Quantification of Compounds

Volatiles were analysed by GC for quantification and by GC-MS for component identification.

3.3.1. Gas Chromatography (GC)

SPME samples. Immediately after sampling, the SPME needle was introduced into the split/splitless injector of a PerkinElmer Clarus 400 gas chromatograph (PerkinElmer, Waltham, MA, USA) equipped with two flame ionisation detectors with a data handling system. Two columns of different polarities were inserted into the injector port: a DB-1 fused-silica column (100% dimethylpolysiloxane, 30 m × 0.25 mm i.d., film thickness 0.25 μm; J & W Scientific Inc., Folsom, CA, USA) and a DB-17HT fused-silica column ((50–phenyl)-methylpolysiloxane, 30 m × 0.25 mm i.d., film thickness 0.15 μm; J & W Scientific). The oven temperature was programmed to rise from 45 to 175 °C at 3 °C/min, then to 300 °C at 15 °C/min and finally to remain isothermal for 10 min, for a total run time of 61.67 min. The SPME fibre was desorbed in splitless mode for 1 min, and the injector and detector temperatures were 250 and 290 °C, respectively; the carrier gas was hydrogen, adjusted to a linear velocity of 30 cm/s. The percentage composition of the volatiles was computed by the normalisation method from the GC peak areas, without the use of correction factors, calculated as mean values of two SPME fibres from each sample.
Hydrodistillation samples. Volatiles obtained from honey and flowers were analysed using the same equipment and settings as described in this section for SPME samples, with the following exceptions: the split injector ratio was 1:40, and the injector temperature was 280 °C. The percentage composition of the volatiles was computed by the normalisation method from the GC peak areas, without the use of correction factors, calculated as mean values of two injections from each sample.

3.3.2. Gas Chromatography–Mass Spectrometry (GC-MS)

SPME samples. Directly after sampling, the SPME needle was introduced into the split/splitless injector of a PerkinElmer Clarus 600 gas chromatograph equipped with a DB-1 fused-silica column (100% dimethylpolysiloxane, 30 m × 0.25 mm i.d., film thickness 0.25 μm; J & W Scientific), interfaced with a PerkinElmer Clarus 600T mass spectrometer (software version 5.4.2.1617, PerkinElmer, Shelton, CT, USA). Injector and oven temperatures were as indicated in Section 3.3.1; the transfer line temperature was 280 °C; the ion source temperature was 220 °C; the carrier gas was helium, adjusted to a linear velocity of 30 cm/s; analyte desorption was achieved in splitless mode for 1 min; the ionisation energy was 70 eV; the scan range was 40–300 u; and the scan time was 1 s.
The identities of the components were assigned by a comparison of their retention indices (RIs) to C6–C31 n-alkane (Sigma) indices and GC-MS spectra from a laboratory-made library based upon the analyses of reference essential oils, laboratory-synthesised components and commercially available standards.
Hydrodistillation samples. The equipment and the analysis settings were as above, with the exceptions already mentioned in Section 3.3.2. In addition, the apocarotenoid edulan could be identified by comparing with one of the authentic samples present in the reference purple passion fruit [65].

3.4. Statistical Analysis

3.4.1. Cluster Analysis

The percentage composition of the isolated volatiles was used to determine the relationship between the different samples by cluster analysis using the Numerical Taxonomy Multivariate Analysis System (NTSYS PC software, version 2.2, Exeter Software, Exeter University, Exeter, UK). For cluster analysis, the correlation coefficient was selected as a measure of similarity among samples and the unweighted pair group method with arithmetical averages (UPGMA) was used for cluster definition. The degree of correlation was evaluated according to [66] as very high [0.90, 1.00], high [0.70, 0.90[, moderate [0.40, 0.70[, low [0.20, 0.40[ and very low (<0.20).

3.4.2. Classification Tree

Considering the botanical sources of honey, a classification tree was built using the volatile components obtained by hydrodistillation. The tree was fully grown with the purpose of investigating whether a small number of volatiles would be enough to nearly completely discriminate the honey types. This approach had a discovery purpose rather than the aim of building a classification tool. The analysis was performed in R (version 4.0.2) and RStudio (version 1.3.1093) using the package rpart [67].

Supplementary Materials

The following are available online: Figure S1: Geographical origin of the studied honey samples according to production locality. Figure S2: Detail of the flowers and extraction of flowers volatiles by hydrodistillation. A,B. Carob tree. C,D. Chestnut. E,F. Eucalyptus. G,H. Bell heather. I. Incense. J. Lavender. K,L. Orange. M,N. Sunflower. O,P. Strawberry tree. Table S1: Volatile profile of the honey samples obtained by solid-phase microextraction (SPME). The minimum and maximum values represent the lowest and highest percentages, respectively, of each identified component. Table S2: Volatile profiles of the honey samples obtained by hydrodistillation (HD). The minimum and maximum values represent the lowest and highest percentages, respectively, of each identified component.

Author Contributions

A.M.M. and A.C.F. conceived and designed the methodology, performed the formal analysis and investigated and wrote the manuscript. A.C.F. supervised the work. A.M.M., A.C.F., M.A., M.G.M. and M.V.-B. read and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Alexandra M. Machado was funded by Fundação para a Ciência e a Tecnologia (FCT) with PhD grant SFRH/BD/117013/2016. This study was partially funded by FCT/MCTES through national funds, under CESAM UIDB/50017/2020 + UIDP/50017/2020, CIMO UIDB/00690/2020, CEAUL UIDB/00006/2020, MED UIDB/05183/2020, FEDER, PT2020 PACompete 2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Alexandra M. Machado is grateful to the Fundação para a Ciência e a Tecnologia (FCT) for the PhD grant SFRH/BD/117013/2016. This study was partially funded by the FCT/MCTES through national funds under CESAM UIDB/50017/2020 + UIDP/50017/2020, CIMO UIDB/00690/2020, CEAUL UIDB/00006/2020, MED UIDB/05183/2020, FEDER and PT2020 PACompete 2020.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the honeys and/or honey extracts are available from the authors.

References

  1. Schievano, E.; Stocchero, M.; Morelato, E.; Facchin, C.; Mammi, S. An NMR-based metabolomic approach to identify the botanical origin of honey. Metabolomics 2012, 8, 679–690. [Google Scholar] [CrossRef]
  2. Bogdanov, S.; Jurendic, T.; Sieber, R.; Gallmann, P. Honey for nutrition and health: A review. J. Am. Coll. Nutr. 2008, 27, 677–689. [Google Scholar] [CrossRef]
  3. Naila, A.; Flint, S.H.; Sulaiman, A.Z.; Ajit, A.; Weeds, Z. Classical and novel approaches to the analysis of honey and detection of adulterants. Food Control 2018, 90, 152–165. [Google Scholar] [CrossRef]
  4. Escuredo, O.; Míguez, M.; Fernández-González, M.; Carmen Seijo, M. Nutritional value and antioxidant activity of honeys produced in a European Atlantic area. Food Chem. 2013, 138, 851–856. [Google Scholar] [CrossRef]
  5. García, N.L. The current situation on the international honey market. Bee World 2018, 95, 89–94. [Google Scholar] [CrossRef]
  6. Codex Alimentarius Commission. Revised Codex Standard for Honey, Codex Stan 12-1981, Rev. 1 (1987), Rev.2 (2001); Food and Agriculture Organization (FAO): Rome, Italy, 2001; pp. 1–7. [Google Scholar]
  7. EU. Council Directive 2001/110/CE concerning Honey. Off. J. Eur. Commun. 2002, 10, 47–52. [Google Scholar]
  8. Siddiqui, A.J.; Musharraf, S.G.; Choudhary, M.I.; Rahman, A. Application of analytical methods in authentication and adulteration of honey. Food Chem. 2017, 217, 687–698. [Google Scholar] [CrossRef]
  9. Zábrodská, B.; Vorlová, L. Adulteration of honey and available methods for detection—A review. Acta Vet. Brno 2015, 83, 85–102. [Google Scholar] [CrossRef] [Green Version]
  10. Ferreiro-González, M.; Espada-Bellido, E.; Guillén-Cueto, L.; Palma, M.; Barroso, C.G.; Barbero, G.F. Rapid quantification of honey adulteration by visible-near infrared spectroscopy combined with chemometrics. Talanta 2018, 188, 288–292. [Google Scholar] [CrossRef]
  11. Jaafar, M.B.; Othman, M.B.; Yaacob, M.; Talip, B.A.; Ilyas, M.A.; Ngajikin, N.H.; Fauzi, N.A.M. A review on honey adulteration and the available detection approaches. Int. J. Integr. Eng. 2020, 12, 125–131. [Google Scholar]
  12. Kilpinen, O.; Vejsnæs, F. New analytical methods against honey fraud are problematic for commercial beekeepers. Bee World 2021, 98, 14–16. [Google Scholar] [CrossRef]
  13. Wu, L.; Du, B.; Vander Heyden, Y.; Chen, L.; Zhao, L.; Wang, M.; Xue, X. Recent advancements in detecting sugar-based adulterants in honey—A challenge. Trends Analyt. Chem. 2017, 86, 25–38. [Google Scholar] [CrossRef]
  14. Fakhlaei, R.; Selamat, J.; Khatib, A.; Razis, A.F.A.; Sukor, R.; Ahmad, S.; Babadi, A.A. The toxic impact of honey adulteration: A review. Foods 2020, 9, 1538. [Google Scholar] [CrossRef] [PubMed]
  15. Phipps, R. International honey market. Am. Bee J. 2021, 161, 1–6. [Google Scholar]
  16. Trifković, J.; Andrić, F.; Ristivojević, P.; Guzelmeric, E.; Yesilada, E. Analytical methods in tracing honey authenticity. J. AOAC Int. 2017, 100, 827–839. [Google Scholar] [CrossRef]
  17. Manyi-Loh, C.E.; Ndip, R.N.; Clarke, A.M. Volatile compounds in honey: A review on their involvement in aroma, botanical origin determination and potential biomedical activities. Int. J. Mol. Sci. 2011, 12, 9514–9532. [Google Scholar] [CrossRef] [Green Version]
  18. Balkanska, R.; Stefanova, K.; Grigorova, R.S. Main honey botanical components and techniques for identification: A review. J. Apic. Res. 2020, 59, 852–861. [Google Scholar] [CrossRef]
  19. Kaškonienė, V.; Venskutonis, P.R. Floral markers in honey of various botanical and geographic origins: A review. Compr. Rev. Food Sci. Food Saf. 2010, 9, 620–634. [Google Scholar] [CrossRef]
  20. Stanimirova, I.; Üstün, B.; Cajka, T.; Riddelova, K.; Hajslova, J.; Buydens, L.M.C.; Walczak, B. Tracing the geographical origin of honeys based on volatile compounds profiles assessment using pattern recognition techniques. Food Chem. 2010, 118, 171–176. [Google Scholar] [CrossRef]
  21. Programa Apícola Nacional 2020–2022. Available online: https://www.gpp.pt/index.php/pan/programa-apicola-nacional-2020-2022 (accessed on 25 May 2021).
  22. Portal Oficial—Instituto Nacional de Estatística. Available online: https://www.ine.pt/xportal/xmain?xpid=INE&xpgid=ine_indicadores&indOcorrCod=0008957&contexto=bd&selTab=tab2 (accessed on 21 June 2021).
  23. Karabagias, I.K.; Karabagias, V.K.; Badeka, A.V. The honey volatile code: A collective study and extended version. Foods 2019, 8, 508. [Google Scholar] [CrossRef] [Green Version]
  24. Rodríguez-Flores, M.S.; Falcão, S.I.; Escuredo, O.; Seijo, M.C.; Vilas-Boas, M. Description of the volatile fraction of Erica honey from the northwest of the Iberian Peninsula. Food Chem. 2021, 336, 127758. [Google Scholar] [CrossRef] [PubMed]
  25. Karabagias, I.K.; Maia, M.; Karabagias, V.K.; Gatzias, I.; Badeka, A.V. Characterization of eucalyptus, chestnut and heather honeys from Portugal using multi-parameter analysis and chemo-calculus. Foods 2018, 7, 194. [Google Scholar] [CrossRef] [Green Version]
  26. Pontes, M.; Marques, J.C.; Câmara, J.S. Screening of volatile composition from Portuguese multifloral honeys using headspace solid-phase microextraction-gas chromatography–quadrupole mass spectrometry. Talanta 2007, 74, 91–103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Aman Dekebo, A.N.U.; Kwon, S.Y.; Kim, D.H.; Jung, C. Volatiles analysis of honey by gas chromatography-mass spectrometry (GC-MS): Comparison of SPME volatiles extraction methods. Korean J. Apic. 2018, 33, 117–128. [Google Scholar]
  28. Jerković, I.; Kuś, P.M.; Carbonell-Barrachina, Á.A. Volatile organic compounds as artefacts derived from natural phytochemicals sourced form plants and honey. Phytochem. Rev. 2019, 18, 871–891. [Google Scholar] [CrossRef]
  29. Custódio, L.; Serra, H.; Nogueira, J.M.F.; Gonçalves, S.; Romano, A. Analysis of the volatiles emitted by whole flowers and isolated flower organs of the carob tree using HS-SPME-GC/MS. J. Chem. Ecol. 2006, 32, 929–942. [Google Scholar] [CrossRef]
  30. Boi, M.; Llorens, J.A.; Cortés, L.; Lladó, G.; Llorens, L. Palynological and chemical volatile components of typically autumnal honeys of the western Mediterranean. Grana 2013, 52, 93–105. [Google Scholar] [CrossRef]
  31. Ouradi, H.; Hanine, H.; Fauconnier, M.-L.; Kenne, T.; Rizki, H.; Ennahli, S.; Hssaini, L. Determination of physico-biochemical proprieties and composition in volatile constituents by solid phase micro-extraction of honey samples from different botanical and geographical origins in Morocco. J. Apic. Res. 2021, 60, 84–98. [Google Scholar] [CrossRef]
  32. Castro-Vázquez, L.; Díaz-Maroto, M.C.; de Torres, C.; Pérez-Coello, M.S. Effect of geographical origin on the chemical and sensory characteristics of chestnut honeys. Food Res. Int. 2010, 43, 2335–2340. [Google Scholar] [CrossRef]
  33. Verzera, A.; Campisi, S.; Zappala, M.; Bonaccorsi, I. SPME-GC-MS analysis of honey volatile components for the characterization of different floral origin. Am. Lab. 2001, 33, 18–21. [Google Scholar]
  34. Siegmund, B.; Urdl, K.; Jurek, A.; Leitner, E. “More than honey”: Investigation on volatiles from monovarietal honeys using new analytical and sensory approaches. J. Agric. Food Chem. 2018, 66, 2432–2442. [Google Scholar] [CrossRef] [PubMed]
  35. Castro-Vázquez, L.; Díaz-Maroto, M.C.; González-Viñas, M.A.; Pérez-Coello, M.S. Differentiation of monofloral citrus, rosemary, eucalyptus, lavender, thyme and heather honeys based on volatile composition and sensory descriptive analysis. Food Chem. 2009, 112, 1022–1030. [Google Scholar] [CrossRef]
  36. Petretto, G.L.; Urgeghe, P.P.; Mascia, I.; Fadda, C.; Rourke, J.P.; Pintore, G. Stir bar sorptive extraction coupled with GC/MS applied to honey: Optimization of method and comparative study with headspace extraction techniques. Eur. Food Res. Technol. 2017, 243, 735–741. [Google Scholar] [CrossRef] [Green Version]
  37. Vázquez, L.C.; Díaz-Maroto, M.C.; Guchu, E.; Pérez-Coello, M.S. Analysis of volatile compounds of eucalyptus honey by solid phase extraction followed by gas chromatography coupled to mass spectrometry. Eur. Food Res. Technol. 2006, 224, 27–31. [Google Scholar] [CrossRef]
  38. Beckmann, K.; Beckh, G.; Lüllmann, C.; Speer, K. Phenylacetaldehyde in Honey—Residue or Natural Compound? Dtsch. Lebensm. Rundsch. 2007, 103, 154–158. [Google Scholar]
  39. Jerković, I.; Prđun, S.; Marijanović, Z.; Zekić, M.; Bubalo, D.; Svečnjak, L.; Tuberoso, C.I.G. Traceability of satsuma mandarin (Citrus unshiu Marc.) honey through nectar/honey-sac/honey pathways of the headspace, volatiles, and semi-volatiles: Chemical markers. Molecules 2016, 21, 1302. [Google Scholar] [CrossRef] [Green Version]
  40. Nozal, M.J.; Bernal, J.L.; Toribio, M.L.; Diego, J.C.; Ruiz, A. Rapid and sensitive method for determining free amino acids in honey by gas chromatography with flame ionization or mass spectrometric detection. J. Chromatogr. A 2004, 1047, 137–146. [Google Scholar] [CrossRef]
  41. Tan, S.T.; Wilkins, A.L.; Holland, P.T.; McGhie, T.K. Extractives from New Zealand unifloral honeys. 2. Degraded carotenoids and other substances from heather honey. J. Agric. Food Chem. 1989, 37, 1217–1221. [Google Scholar] [CrossRef]
  42. Guyot, C.; Scheirman, V.; Collin, S. Floral origin markers of heather honeys: Calluna vulgaris and Erica arborea. Food Chem. 1999, 64, 3–11. [Google Scholar] [CrossRef]
  43. Guyot-Declerck, C.; Renson, S.; Bouseta, A.; Collin, S. Floral quality and discrimination of Lavandula stoechas, Lavandula angustifolia, and Lavandula angustifolia×latifolia honeys. Food Chem. 2002, 79, 453–459. [Google Scholar] [CrossRef]
  44. Kadar, M.; Juan-Borrás, M.; Carot, J.M.; Domenech, E.; Escriche, I. Volatile fraction composition and physicochemical parameters as tools for the differentiation of lemon blossom honey and orange blossom honey. J. Sci. Food Agric. 2011, 91, 2768–2776. [Google Scholar] [CrossRef]
  45. Alissandrakis, E.; Tarantilis, P.A.; Harizanis, P.C.; Polissiou, M. Aroma investigation of unifloral Greek citrus honey using solid-phase microextraction coupled to gas chromatographic–mass spectrometric analysis. Food Chem. 2007, 100, 396–404. [Google Scholar] [CrossRef]
  46. Castro-Vázquez, L.; Díaz-Maroto, M.C.; Pérez-Coello, M.S. Aroma composition and new chemical markers of Spanish citrus honeys. Food Chem. 2007, 103, 601–606. [Google Scholar] [CrossRef]
  47. Jabalpurwala, F.A.; Smoot, J.M.; Rouseff, R.L. A comparison of citrus blossom bolatiles. Phytochemistry 2009, 70, 1428–1434. [Google Scholar] [CrossRef]
  48. Ruisinger, B.; Schieberle, P. Characterization of the key aroma compounds in rape honey by means of the molecular sensory science concept. J. Agric. Food Chem. 2012, 60, 4186–4194. [Google Scholar] [CrossRef] [PubMed]
  49. Seisonen, S.; Kivima, E.; Vene, K. Characterisation of the aroma profiles of different honeys and corresponding flowers using solid-phase microextraction and gas chromatography–mass spectrometry/olfactometry. Food Chem. 2015, 169, 34–40. [Google Scholar] [CrossRef]
  50. Kruží, V.; Grégrová, A.; Ziková, A.; Čížková, H. Rape honey: Determination of botanical origin based on volatile compound profiles. J. Food Nutr. Res. 2019, 58, 339–348. [Google Scholar]
  51. Pažitná, A.; Janáčová, A.; Špánik, I. The enantiomer distribution of major chiral volatile organic compounds in slovakian monofloral honeys. J. Food Nutr. Res. 2012, 51, 235–241. [Google Scholar]
  52. Soria, A.C.; Sanz, J.; Martínez-Castro, I. SPME Followed by GC–MS: A powerful technique for qualitative analysis of honey volatiles. Eur. Food Res. Technol. 2009, 228, 579–590. [Google Scholar] [CrossRef] [Green Version]
  53. Juan-Borrás, M.; Domenech, E.; Hellebrandova, M.; Escriche, I. Effect of country origin on physicochemical, sugar and volatile composition of acacia, sunflower and tilia honeys. Food Res. Int. 2014, 60, 86–94. [Google Scholar] [CrossRef]
  54. Graikou, K.; Andreou, A.; Chinou, I. Chemical profile of Greek Arbutus unedo honey: Biological properties. J. Apic. Res. 2021, 1–7. [Google Scholar] [CrossRef]
  55. Bianchi, F.; Careri, M.; Musci, M. Volatile norisoprenoids as markers of botanical origin of Sardinian strawberry-tree (Arbutus unedo L.) honey: Characterisation of aroma compounds by dynamic headspace extraction and gas chromatography–mass spectrometry. Food Chem. 2005, 89, 527–532. [Google Scholar] [CrossRef]
  56. Karabagias, I.K.; Nikolaou, C.; Karabagias, V.K. Volatile fingerprints of common and rare honeys produced in Greece: In search of PHVMs with implementation of the honey code. Eur. Food Res. Technol. 2019, 245, 23–39. [Google Scholar] [CrossRef]
  57. Fröhlich, B.; Tautz, J.; Riederer, M. Chemometric classification of comb and cuticular waxes of the honeybee Apis mellifera carnica. J. Chem. Ecol. 2000, 26, 123–137. [Google Scholar] [CrossRef]
  58. Teerawanichpan, P.; Robertson, A.J.; Qiu, X. A fatty acyl-CoA reductase highly expressed in the head of honey bee (Apis mellifera) involves biosynthesis of a wide range of aliphatic fatty alcohols. Insect Biochem. Mol. Biol. 2010, 40, 641–649. [Google Scholar] [CrossRef]
  59. Thom, C.; Gilley, D.C.; Hooper, J.; Esch, H.E. The scent of the waggle dance. PLoS Biol. 2007, 5, e228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Jerković, I.; Marijanović, Z.; Ljubičić, I.; Gugić, M. Contribution of the bees and combs to honey volatiles: Blank-trial probe for chemical profiling of honey biodiversity. Chem. Biodivers. 2010, 7, 1217–1230. [Google Scholar] [CrossRef] [PubMed]
  61. Castro-Vázquez, L.; Leon-Ruiz, V.; Alañon, M.E.; Pérez-Coello, M.S.; González-Porto, A.V. Floral origin markers for authenticating Lavandin honey (Lavandula angustifolia x latifolia). Discrimination from Lavender honey (Lavandula latifolia). Food Control 2014, 37, 362–370. [Google Scholar] [CrossRef]
  62. Bobis, O.; Moise, A.R.; Ballesteros, I.; Reyes, E.S.; Durán, S.S.; Sánchez-Sánchez, J.; Cruz-Quintana, S.; Giampieri, F.; Battino, M.; Alvarez-Suarez, J.M. Eucalyptus honey: Quality parameters, chemical composition and health-promoting properties. Food Chem. 2020, 325, 126870. [Google Scholar] [CrossRef]
  63. Jerkovic, I.; Kranjac, M.; Suste, M.; Kus, P.M.; Svecnjak, L. Rhamnus frangula honey: Screening of volatile organic compounds and their composition after short-term heating. Chem. Nat. Compd. 2015, 51, 1174–1177. [Google Scholar] [CrossRef]
  64. Machado, A.M.; Tomás, A.; Russo-Almeida, P.; Duarte, A.; Antunes, M.; Vilas-Boas, M.; Miguel, M.G.; Figueiredo, A.C. Portuguese monofloral honeys. Physicochemical parameters as tools in botanical source differentiation. Food Chem. 2021. under review. [Google Scholar]
  65. Whitfield, F.B.; Stanley, G. The structure and stereochemistry of Edulan I and II and the stereochemistry of the 2,5,5,8a-Tetramethyl-3,4,4a,5,6,7,8,8a-octahydro-2H-1-benzopyrans. Aust. J. Chem. 1977, 30, 1073–1091. [Google Scholar] [CrossRef]
  66. Pestana, M.H.; Gageiro, J.N. Análise De Dados Para Ciências Sociais: A Complementaridade do SPSS; Edições Sílabo: Lisboa, Portugal, 2000. [Google Scholar]
  67. Recursive Partitioning and Regression Trees [R Package Rpart Version 4.1-15]. Available online: https://CRAN.R-project.org/package=rpart (accessed on 4 June 2021).
Figure 1. Extraction of volatiles from different honey samples. (AC) SPME collection of volatiles from carob tree honey (A), eucalyptus honey (B) and lavender honey (C). (DF) Hydrodistillation of orange honey (D), sunflower honey (E) and strawberry tree honey (F).
Figure 1. Extraction of volatiles from different honey samples. (AC) SPME collection of volatiles from carob tree honey (A), eucalyptus honey (B) and lavender honey (C). (DF) Hydrodistillation of orange honey (D), sunflower honey (E) and strawberry tree honey (F).
Molecules 26 04970 g001
Figure 2. Dendrogram obtained by cluster analysis of the percentage composition of volatile compounds isolated by hydrodistillation from the 51 honey samples, based on correlation, and using the unweighted pair-group method with arithmetic average (UPGMA). Ct: carob tree; C: chestnut; E: eucalyptus; H: bell heather; I: incense; L: lavender; O: orange; R: rape; Rb: raspberry; Ro: rosemary; Sf: sunflower; St: strawberry tree.
Figure 2. Dendrogram obtained by cluster analysis of the percentage composition of volatile compounds isolated by hydrodistillation from the 51 honey samples, based on correlation, and using the unweighted pair-group method with arithmetic average (UPGMA). Ct: carob tree; C: chestnut; E: eucalyptus; H: bell heather; I: incense; L: lavender; O: orange; R: rape; Rb: raspberry; Ro: rosemary; Sf: sunflower; St: strawberry tree.
Molecules 26 04970 g002
Figure 3. Classification tree with the honey types separated in nodes, according to volatile compounds’ discrimination. * Identification based on mass spectra only. UI: unidentified compounds; C: chestnut; Ct: carob tree; E: eucalyptus; H: bell heather; I: incense. L: lavender; O: orange; R: rape; Rb: raspberry; Sf: sunflower; St: strawberry tree. Nodes should be read as follows: The letter at the top is the node label, corresponding to the most frequent honey type in the node; the numbers represent the number of samples of each of the 12 honey types in the node (honey types in alphabetic order); and in the fourth line is the percentage of honey samples in the node (out of the 51 samples under analysis).
Figure 3. Classification tree with the honey types separated in nodes, according to volatile compounds’ discrimination. * Identification based on mass spectra only. UI: unidentified compounds; C: chestnut; Ct: carob tree; E: eucalyptus; H: bell heather; I: incense. L: lavender; O: orange; R: rape; Rb: raspberry; Sf: sunflower; St: strawberry tree. Nodes should be read as follows: The letter at the top is the node label, corresponding to the most frequent honey type in the node; the numbers represent the number of samples of each of the 12 honey types in the node (honey types in alphabetic order); and in the fourth line is the percentage of honey samples in the node (out of the 51 samples under analysis).
Molecules 26 04970 g003
Table 1. Monofloral-labelled honeys studied, their production regions and analysed flower types, with corresponding geographical origins and volatile components found in each honey’s botanical source.
Table 1. Monofloral-labelled honeys studied, their production regions and analysed flower types, with corresponding geographical origins and volatile components found in each honey’s botanical source.
Honey SamplesLocality/Region of Honey ProductionFlower’s Geographical OriginVolatile Components Present in Selected Flowers and Simultaneously in Corresponding Honeys
Carob tree (Ct)Estremaduraα-Pinene, cis- and trans-linalool oxide (furanoid), cis- and trans-linalool oxide (pyranoid) *, UI 8, methyl anthranilate, 1,2-dihydro-1,1,6-trimethyl-naphtalene *, α-eudesmol
Ct1Lagos/Algarve
Ct2Tavira/Algarve
Ct3Olhão/Algarve
Ct4Tavira/Algarve
Ct5Tavira/Algarve
Chestnut (C)
C1Macedo de Cavaleiros/Trás-os-Montes e Alto DouroTrás-os-Montes e Alto DouroBenzaldehyde, α-pinene, benzene acetaldehyde, acetophenone, cis- and trans-linalool oxide (furanoid), UI 8, 1,2-dihydro-1,1,6-trimethyl-naphtalene *, α-eudesmol, β-eudesmol
C2Vila Pouca de Aguiar/Trás-os-Montes e Alto Douro
Eucalyptus (E)
E1Ponte de Lima/MinhoBeira Litoralα-Pinene, cis-linalool oxide (furanoid), aromadendrene, α-eudesmol, β-eudesmol
E2Penamacor/Beira Baixa
E3Marvão/Alto Alentejo
E4Vila Pouca de Aguiar/Trás-os-Montes e Alto Douro
E5Sesimbra/Estremadura
Bell heather (H)
H1Fundão/Beira BaixaTrás-os-Montes e Alto DouroBenzaldehyde, benzene acetaldehyde, UI 8, 1,2-dihydro-1,1,6-trimethyl-naphtalene *
H2Boticas/Trás-os-Montes e Alto Douro
H3Pampilhosa da Serra/Beira Baixa
H4Ponte da Barca/Minho
H5Vila Pouca de Aguiar/Trás-os-Montes e Alto Douro
H6Melgaço/Minho
Incense (I)
I1São Miguel/AçoresSão Miguel/AçoresBenzyl salicylate *
I2Pico/Açores
I3São Miguel/Açores
I4Pico/Açores
Lavender (L)
L1Castelo Branco/Beira BaixaBeira BaixaNo specific components found
L2Castelo Branco/Beira Baixa
L3Macedo de Cavaleiros/Trás-os-Montes e Alto Douro
L4Mértola/Baixo Alentejo
L5Penamacor/Beira Baixa
L6Évora/Alto Alentejo
L7Seixal/Estremadura
L8Sesimbra/Estremadura
Orange (O)
O1Faro/AlgarveEstremaduraα-Pinene; lilac aldehydes A *, B *, C *; UI 8; cis-myrcenol; indole; methyl anthranilate
O2Tavira/Algarve
O3Sesimbra/Estremadura
O4Aljezur/Algarve
O5Odemira/Alentejo Litoral
O6Odemira/Alentejo Litoral
O7Lagos/Algarve
O8Odemira/Alentejo Litoral
O9Portimão/Algarve
Rape (R)
R1Odemira/Alentejo Litoral
R2Odemira/Alentejo Litoral
Raspberry (Rb)
Rb1Sesimbra/Estremadura
Rb2Odemira/Alentejo Litoral
Rosemary (Ro)
Ro1Sesimbra/Estremadura
Sunflower (Sf)
Sf1Odemira/Alentejo LitoralBaixo Alentejoα-Pinene
Sf2Aljezur/Algarve
Sf3Ferreira do Alentejo/Baixo Alentejo
Strawberry tree (St)
St1Aljustrel/Baixo AlentejoEstremaduraα-Isophorone *, 4-keto-isophorone *, UI 8, 2,3,5-trimethylphenol, 3,4,5-trimethylphenol, 1,2-dihydro-1,1,6-trimethyl-naphtalene *
St2Ferreira do Alentejo/Baixo Alentejo
St3Lagos/Algarve
St4Aljezur/Algarve
Total51
* Identification based on mass spectra only; UI: unidentified compounds; † unavailable for extraction.
Table 2. Minimum and maximum percentage range of the main components (≥5%) isolated by hydrodistillation from the different monofloral honeys grouped according to cluster analysis. For samples grouped in each of the clusters and sub-clusters, see Figure 2. Fully detailed composition is provided in SM Table S2.
Table 2. Minimum and maximum percentage range of the main components (≥5%) isolated by hydrodistillation from the different monofloral honeys grouped according to cluster analysis. For samples grouped in each of the clusters and sub-clusters, see Figure 2. Fully detailed composition is provided in SM Table S2.
ComponentsRICluster ICluster II
IaIbIcIdIeIfIgIhIi
Min.Max. Min.Max.Min.Max.Min.Max.Min.Max.Min.Max.
Benzene acetaldehyde1002 9.98.83.60.70.10.10.41.13.80.23.117.324.40.10.3
cis-Linalool oxide (furanoid)1045 6.5t0.5t0.10.10.39.015.10.22.92.94.3 t
trans-Linalool oxide (furanoid)1059 2.4 t t0.10.13.95.1t1.71.62.0 t
Hotrienol *1074 5.10.10.50.10.1 2.99.8 0.33.26.9
α-Isophorone (2-cyclohexen-1-one, 3,5,5-trimethyl-)1074 4.4 0.1t 0.1 0.1 0.1 2.636.438.8
3,4,5-Trimethylphenol (3,4,5-hemimellitenol)1277 1.5 ttt0.1 t 0.2 0.64.413.3
Decanoic acid1356 2.2t t tt0.26.4 t0.53.4t0.4
n-Heptadecane17001.45.91.03.32.87.91.22.01.03.90.82.61.21.30.30.6
n-Nonadecane19002.414.71.84.56.025.62.13.81.46.12.05.12.93.10.61.8
Hexadecanoic acid (palmitic acid)1908 7.41.42.113.45.42.54.9 2.23.37.70.62.71.92.9
Oleyl alcohol * [(Z)-octadec-9-en-1-ol]2044 0.3 11.40.1
1-Octadecanol (stearyl alcohol)2095 4.50.9 0.3 3.1 1.11.43.4
Heneicosene *2088 5.4 t0.40.21.34.1 2.0 1.6 0.31.4
Methyl oleate (methyl cis-9-octadecenoate)2096 2.16.3 t t 1
n-Heneicosane21003.311.42.44.36.815.124.025.51.35.02.74.02.33.12.97.0
Oleic acid (cis-9-octadecenoic acid)2119 13.71.82.12.6 t1.01.33.917.338.80.42.60.88.7
Linoleic acid (cis-9, cis-12-octadecadienoic acid)2140 2.3 1.3 9.4
Ethyl oleate2151 10.81.11.72.00.6 0.3 9.3 2.1
(Z)-9-Tricosene2287 6.1 2.1
n-Tricosane230013.929.312.212.914.613.49.413.86.012.67.514.87.38.01.22.5
n-Pentacosane25000.316.16.06.26.77.08.711.32.76.04.47.83.05.01.52.3
Heptacosene2667 6.21.21.70.50.2tt0.81.40.15.30.52.0
n-Heptacosane2700 8.22.62.12.02.74.35.81.12.3t7.11.43.41.02.6
n-Octacosane2800 4.6 0.1 0.3 0.4
% Identification 69.696.470.068.887.095.674.190.273.289.581.590.486.589.580.288.7
Grouped components
Terpenes and derivatives
Hemiterpene hydrocarbons ttt t ttt
Monoterpene hydrocarbons 1.2 0.1 0.1t0.7 0.5t0.5
Oxygen-containing monoterpenes 12.10.31.00.10.21.92.219.927.41.0.5.410.511.60.10.6
Sesquiterpene hydrocarbons 0.70.31.80.2t 0.1 0.2t0.2 t
Oxygen-containing sesquiterpenes 1.80.5ttttt t 0.8t0.2t4.2
Oxygen-containing diterpenes 0.4 tt t
Apocarotenoids 5.2tt0.1t 0.10.40.7 0.10.63.741.644.3
Amino acid derivatives and phenylpropanoids
Benzoic acid derivatives 1.31.01.1tttt0.20.5 0.71.62.6t0.3
Phenylpropenes 1.01.1t tt 1.1 0.3t0.5tt
Aromatic amino acid derivatives0.19.98.83.61.10.11.42.01.13.80.93.117.325.27.418.6
Fatty acids and derivatives
Green leaf volatiles (GLVs) 0.1ttttttt0.1tt
Fatty acids1.018.63.54.218.05.43.55.19.011.420.645.33.79.62.812.1
Alkanes40.874.227.737.143.276.955.765.715.235.924.542.722.227.011.214.9
Methyl-branched hydrocarbons 1.20.80.80.50.20.30.8 0.9 1.50.30.4t0.1
Other fatty acid derivatives6.628.525.618.223.412.68.816.610.214.310.717.414.518.12.04.4
Carbohydrate derivatives 2.0t1.00.20.1t0.12.02.8t0.71.92.20.20.9
Nitrogen-containing compounds t t0.1
Sulphur-containing compounds t t
Others 2.30.4 t3.5 t0.41.2t1.4
RI: In-lab calculated retention index relative to C6-C31 n-alkanes on the DB-1 column. Min.: minimum; Max.: maximum; t: traces (<0.05%). * Identification based on mass spectra only.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Machado, A.M.; Antunes, M.; Miguel, M.G.; Vilas-Boas, M.; Figueiredo, A.C. Volatile Profile of Portuguese Monofloral Honeys: Significance in Botanical Origin Determination. Molecules 2021, 26, 4970. https://doi.org/10.3390/molecules26164970

AMA Style

Machado AM, Antunes M, Miguel MG, Vilas-Boas M, Figueiredo AC. Volatile Profile of Portuguese Monofloral Honeys: Significance in Botanical Origin Determination. Molecules. 2021; 26(16):4970. https://doi.org/10.3390/molecules26164970

Chicago/Turabian Style

Machado, Alexandra M., Marília Antunes, Maria Graça Miguel, Miguel Vilas-Boas, and Ana Cristina Figueiredo. 2021. "Volatile Profile of Portuguese Monofloral Honeys: Significance in Botanical Origin Determination" Molecules 26, no. 16: 4970. https://doi.org/10.3390/molecules26164970

APA Style

Machado, A. M., Antunes, M., Miguel, M. G., Vilas-Boas, M., & Figueiredo, A. C. (2021). Volatile Profile of Portuguese Monofloral Honeys: Significance in Botanical Origin Determination. Molecules, 26(16), 4970. https://doi.org/10.3390/molecules26164970

Article Metrics

Back to TopTop