Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 22 July 2022
Sec. Analytical Chemistry

Characterization of Ikaria Heather Honey by Untargeted Ultrahigh-Performance Liquid Chromatography-High Resolution Mass Spectrometry Metabolomics and Melissopalynological Analysis

  • 1Laboratory of Pesticides’ Toxicology, Benaki Phytopathological Institute, Athens, Greece
  • 2General Chemical State Laboratory, Independent Public Revenue Authority (A.A.D.E), Athens, Greece

Honey represents a valuable food commodity, known since ancient times for its delicate taste and health benefits due to its specific compositional characteristics, mainly the phenolic compound content. “Anama” honey is a monofloral honey produced from the nectar of Erica manipuliflora plant, a heather bush of the Greek island of Ikaria, one of the Mediterranean’s longevity regions. “Anama” is characterized by a unique aroma and taste, with a growing demand for consumption and the potential to be included in the list of products with a protected designation of origin. The aim of this study was to determine the chemical and botanical profile of authentic Anama honey samples and find similarities and differences with honey samples of a different botanical origin from the same geographical area. Untargeted Ultrahigh-Performance Liquid Chromatography-Hybrid Quadrupole-Orbitrap High-Resolution Mass Spectrometry (UHPLC-HRMS) metabolomics study was conducted on authentic heather, pine, and thyme honey samples from Ikaria and neighboring islands. The Principal Component Analysis (PCA), Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA), and differential analysis were performed using the entire metabolic profile of the samples and allowed the identification of chemical markers for sample discrimination. Thirty-two characteristic secondary metabolites (cinnamic acids, phenolic acids, flavonoids, terpenes) and other bioactive phenolic compounds, some of them not previously reported in a heather honey (aucubin, catalpol, domesticoside, leonuriside A, picein among others), emerged as potential chemical indicators of Anama honey. Melissopalynological analysis was also carried out to decipher the botanical and geographical origin of Anama honey. The relative frequency of the pollen of dominant plants of the Ericaceae family and a multitude of nectariferous and nectarless plants contributing to the botanical profile of Anama was evaluated. The identification of the pollen sources enabled a potential correlation of differentially increased secondary metabolites and chemicals with their botanical origin. The physicochemical profile of Anama was also determined, including the parameters of pH, color, electrical conductivity, diastase, moisture, as well as sugars, supporting the high quality of this heather honey.

1 Introduction

Honey is a natural product acclaimed since ancient times both for its taste and its various beneficial effects on human health (Hills et al., 2019; Kamaruzzaman et al., 2019; Nguyen et al., 2019). Though not consumed in large quantities (the average consumption of honey per capita and per day in Europe is less than 5 g/capita/day) (SANTE/11956/2016, 2018), honey and other apiculture products (e.g., royal jelly) are still an important source of valuable nutrients for humans, as they are extremely rich in bioactive components. The main components of honey are fructose and glucose, and also many other disaccharides, such as maltose and sucrose. Honey also contains small amounts of vitamins and minerals, as well as organic acids, amino acids, proteins, enzymes, and polyphenols, mainly flavonoids, and phenolic acids. The chemical composition of honey and its content of bioactive compounds is expected to vary to a great extent depending on the floral nectar that bees forage on, which, in turn, is influenced by the plant species and geographical origin. Phenolic compounds responsible for antioxidant and anti-inflammatory capacity of honey are secondary metabolites synthesized by plants and could be used as botanical markers for classification and authentication, especially in the case of monofloral honeys (Ciulu et al., 2016). Monofloral (or unifloral) honeys are considered those that come wholly or mainly from the indicated plant source and possess the organoleptic, physicochemical, and microscopic characteristics of the source (Council-Directive, 2001; Oddo et al., 2004). They must also conform to existing legislative criteria (Thrasyvoulou et al., 2018). Some monofloral honeys are more appreciated than others due either to their flavor and aroma properties or to their pharmacological attributes. The most famous monofloral honey used for medical purposes is Manuka honey from New Zealand, but several studies have revealed that other monofloral honeys, such as heather honey, are also very effective (Salonen et al., 2017; Gourdomichali and Papakonstantinou, 2018). Heather honey, obtained from plants that belong to Ericaceae family, is relatively rare and very appreciated by consumers, due to its sensorial properties and potential health benefits (Cianciosi et al., 2018; Rodríguez-Flores et al., 2021). Greece is a Mediterranean country with a high production of heather honey. Several different species of Erica plants contribute to the production of this honey type, especially Erica arborea and E. manipuliflora, among others. In the last 5 years, the monofloral heather honey from Erica manipuliflora plant of Ikaria (and Fournoi) island, one of the Mediterranean’s longevity regions, called “Anama” (Greek, Άναμα), is at the forefront of the Greek apiculture and honey community due to its unique flavor and organoleptic features. It has an intense aroma and its taste is pungent, smoky, slightly bitter, and mildly sweet, which leaves a long aftertaste. It is on the list of three products from the North Aegean prefecture to be assessed for products of protected designation of origin and the relevant registration file will be prepared by the Directorate of Rural Development of the North Aegean Region (Anama-honey, 2021).

The current Regulation (EU) 625/2017 on official food control prioritizes the control of food authenticity (REG-625/2017, 2017), including protected designations of origin, geographical indications, and traditional specialties (Council-Directive, 2001). A wide range of analytical methods has been developed for food characterization, quality, and authenticity testing; yet, in the last few years, novel methods principally based on the advances in mass spectrometry (MS) techniques are applied (Cook and Nightingale, 2018; Braconi et al., 2021). These omics-based methods are moving toward the maximum collection of information (chemical context, markers) on sample composition, fostering substantially this domain with proven application in honey commodities (Koulis et al., 2021). Complementary to this, and for honey, in particular, one of the most important analyses is the melissopalynological (pollen) analysis, which comprises the study of the whole range of plants that bees visit, by identifying the pollen deposited in the honey and extracting information about the botanical and geographical origin of the honey. Pollen analysis can also provide evidence of the presence of exogenous elements in honey, such as those that can result from bee feed, sometimes of dubious quality (von der Ohe et al., 2004). In the same direction, metabolomic approaches in combination with robust statistical analysis are pivotal tools in distinguishing diverse types of honey or investigating their taste and maturation attributes, disclosing characteristic-unique chemicals for the discrete categories (Guo, 2021; Koulis et al., 2021; Montoro et al., 2021).

To our knowledge, the already published data regarding the Greek heather honey chemical profile is very limited. More specifically, a recently published work has incorporated Ultrahigh Performance Liquid Chromatography High-Resolution Mass Spectrometry (UHPLC-HRMS) targeted and untargeted metabolomics analysis to discriminate between honey samples of different botanical sources from Greece and Poland. However, the exact geographical origin of the Greek samples was not provided, while melissopalynological and physicochemical analyses were not conducted (Koulis et al., 2021). Another work performed by using solid-phase microextraction and gas chromatography–mass spectrometry has been focused on volatile constituents of Greek heather honey (Xagoraris et al., 2021).

In this study, heather honey samples from Ikaria Island were investigated through an untargeted UHPLC-HRMS metabolomics study to elucidate the profile of secondary metabolites and chemicals. Using advanced chemometrics tools, the results were compared to the metabolomic results of pine and thyme honey from the same geographical origin area collected at different times of the year, to unveil potential chemical markers of this unique type of Anama honey. Extensive melissopalynological analysis performed on all honey samples, revealed a diversity of plants contributing to Anama honey composition when compared to the other Aegean monofloral pine and thyme honeys. For the first time, the identification of the pollen sources enabled a potential correlation of differentially increased secondary metabolites and chemicals with their botanical sources. Physicochemical parameters including the determination of pH, color, electrical conductivity, diastase, moisture, as well as sugars, were also performed, supporting the high quality of heather Anama honey.

2 Materials and Methods

2.1 Honey Samples

Authentic honey samples (n = 50) were provided by the beekeeping and rural cooperatives of five Aegean islands for consecutive years 2018, 2019, and 2020. Samples were collected immediately after harvesting and kept at -18°C until analysis. Harvesting was conducted in three different periods of the year, according to the botanical origin of honey.

Anama honey samples (n = 10) were collected from the islands of Ikaria and Fournoi between November and January 2018–2019 and 2019–2020. These islands were considered as one group due to their vicinity, and transferring of hives from one island to the other by the beekeepers. Thyme honey samples (n = 10) were also collected from Fournoi (June to August of 2018, 2019, and 2020), accompanied by the same number of samples from Kea (n = 10) and Syros (n = 10) islands. Pine honey samples (n = 10) were collected from Samos Island (August 2018–2020), the closest island to Ikaria and Fournoi islands.

2.2 Chemicals and Reagents

Acetonitrile (ACN), methanol (MeOH), and formic acid of liquid chromatography–mass spectrometry (LC-MS) grade were obtained from Merck (Darmstadt, Germany). Ultra-pure water (H2O) was produced from SG Millipore apparatus. Abscisic acid was obtained from Sigma Aldrich (Seelze, Germany). Discovery® octadecylsilane (DSC-C18) solid-phase extraction (SPE) cartridges (500 mg, 6 ml) were purchased from Sigma Aldrich (Supelco, Bellefonte, PA, United States). Nylon filters (0.22 μm) were obtained from Macherey-Nagel (Duren, Germany). All reagents and chemicals were of analytical grade. (+)-Abscisic acid analytical standard (≥ 98%) was obtained from Sigma Aldrich (Seelze, Germany).

2.3 Sample Preparation for UHPLC-HRMS

Sample preparation was adapted from a published procedure (Di Marco et al., 2018). Briefly, honey (3 g) was homogenized with the aid of a glass rod in 7 ml of acidified ultrapure H2O (pH = 2.2). Then, the mixture was loaded on a pre-activated DSC-C18 SPE cartridge [activation was performed by succeeding elution of MeOH (3 ml) and H2O (3 ml)]. Consequently, the SPE cartridge was washed with acidified H2O (2 ml) and 5 ml of ultrapure H2O. Compounds of interest were eluted with 3.5 ml of a MeOH/ACN (2:1, v/v) solution. Then, the eluate was evaporated to dryness using a rotary evaporator (bath temperature not exceeding 30°C) and reconstituted in MeOH/H2O (7:3, v/v). After filtration with Nylon filters (0.22 μm), the extract was injected into UHPLC-HRMS system.

2.4 Ultrahigh Performance Liquid Chromatography Orbitrap High-Resolution Mass Spectrometry Analysis of Honey Extracts

Samples were analyzed with a Dionex Ultimate 3000 UHPLC system (Thermo Fisher Scientific, San Jose, CA, United States) linked to a Q-Exactive Orbitrap HRMS (Thermo Fisher Scientific, San Jose, CA, United States). Chromatographic separation was performed on a Hypersil Gold UPLC C18 (2.1 mm × 150 mm, 1.9 μm) reversed phased column (Thermo Fisher Scientific, San Jose, CA, United States). The ionization was performed using heated electrospray (ESI) in both positive and negative modes. Chromatographic and mass spectrometry conditions and parameters are described in the previous work of our group (Stavropoulou et al., 2021). To monitor the stability and repeatability of the instrumental procedure, a quality control (QC) sample was analyzed three times before and at the end of the sequence, and after every six experimental samples. The QC sample was prepared as a pool of all study samples.

2.5 Data Processing and Multivariate Statistical Analysis

Data were processed with Compound Discoverer 2.1 (Thermo Fisher Scientific, San Jose, CA, United States). The workflow was an untargeted metabolomics approach for detecting unknown compounds and searching in online databases. The software was used for peak detection, deconvolution, deisotoping, alignment, gap-filling, and composition prediction procedures. Peaks in which the QC sample coverage was less than 50% and the relative standard deviation of the areas under the peaks was more than 30% were excluded. The generated peak lists of accurate masses and retention times paired with corresponding intensities for all detected peaks were exported as a .csv file, imported to Microsoft Excel 2010, and manipulated appropriately using the “concatenate,” “round” and “transpose” commands. The R package MetaboDiff (Mock et al., 2018) was implemented for the differential analysis of metabolomic data and for the calculation of Log2 Fold Changes and adjusted p-values (Benjamini-Hochberg).

The MS data were then subjected to multivariate statistical analysis (MVA) i.e., Principal Component Analysis (PCA) and Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) using the ropls R Bioconductor package (Thévenot et al., 2015) to determine the optimal number of components, confirm the validity of the model by permutation testing, detect outliers and perform feature selection with Variable Importance in Projection (VIP) scoring from OPLS-DA models. Variables that exhibited a VIP score greater than 1.5, combined with an adjusted p-value ≤0.05 and Log2 Fold Change (Anama vs. Thyme or Anama vs. Pine) >1.5 from differential analysis, were considered as contributing the most to group separation. The common differentially increased variables on Anama honey compared to thyme and pine honey were forwarded to putatively annotation. Databases such as mzCloud (https://www.mzcloud.org/), MassBank (www.massbank.jp), and Metlin (http://metlin.scripps.edu) were used for putative annotation by applying m/z tolerance of 5 ppm and taking into consideration the isotope distribution similarity and MS/MS fragmentation pattern deriving from the Data Dependent Acquisition (DDA) methodology. For the compounds for which there was no record in the aforementioned databases in silico fragmentation, the approach was applied using the MetFrag online tool, using 5 ppm search tolerance and 0.001 mass deviations to match generated fragments against MS/MS peaks (Wolf et al., 2010). For the MetFrag workflow, the candidate structures were retrieved from the databases PubChem (PubChem, n.d.) and Kegg. For the selection of the compounds >0.8 final scores were applied. The score for every matching is presented in Supplementary Table S1.

2.6 Melissopalynological Analysis

For the Melissopalynological analysis, the methods of Louveaux & von der Ohe were applied (von der Ohe et al., 2004; Louveaux et al., 1978). In this respect, a solution of 10 g of honey in 20 ml of distilled water is centrifuged for 10 min at 1,000 g. After decanting the supernatant liquid, the sediment is suspended again in 20 ml of distilled water and centrifuged for 5 min at 1,000 g. The supernatant liquid is decanted, and the sediment is spread over a microscope slide on a heating plate (40°C). The sediment is left to dry and then is colored with some drops of a solution of fuchsine (∼0.05‰ w/v in ethanol/water 1/1). After drying (40°C), a coverslip with a drop of glycerin jelly (mounting medium; Kaiser’s Glycerol Gelatin TM Merck 1.09242.0100) is placed upon the sediment and the slide is left on the heating plate for another 5 min. The microscope slide is investigated under the microscope (400x). For the determination of relative frequencies of pollen types, 500 to 1000 pollen grains are counted (Behm et al., 1996). To detect all possible plant sources, especially those that exhibit very large pollen grains, which are not expected to be present in a substantial amount, the entire slide was scanned at 50x magnification. Quantitative melissopalynological analysis was achieved using Slide Grids, 20 mm × 20 mm. At least 10 mm2 were counted.

2.7 Physicochemical Characteristics

The determination of physicochemical parameters of Anama honey samples was achieved by applying the “Harmonised Methods of the International Honey Commission (IHC)” (IHC, 2009). Specifically, the following protocols of IHC methods were applied: The refractive index (RI) of honey was measured using a refractometer. The water content was determined from the refractive index of honey by reference to a standard table. The electrical conductivity of a solution of 20 g dry matter of honey in 100 ml distilled water was measured using an electrical conductivity cell, at 20°C. The pH of a solution of 10 g honey in 75 ml of carbon dioxide-free water was measured using a pH meter, accurate to 0.01 units. The determination of the hydroxymethylfurfural (HMF) content was based on the measurement of UV absorbance of HMF at 284 nm (after White). In order to avoid the interference of other components at this wavelength, the difference between the absorbances of a clear aqueous honey solution and the same solution after the addition of bisulfite was determined. The HMF content was calculated after the subtraction of the background absorbance at 336 nm. The determination of the diastase activity of honey was performed by a photometric method in which an insoluble blue-dyed cross-linked type of starch was used as the substrate (Phadebas tablets). This was hydrolyzed by the enzyme, yielding blue water-soluble fragments, determined photometrically at 620 nm. The absorbance of the solution is directly proportional to the diastase activity of the sample. The determination of sugars was performed in a filtrated solution of 10 g of honey in 100 ml of water by High-Pressure Liquid Chromatography (HPLC) using a RI-detector. Peaks were identified on the basis of their retention times. Quantitation was performed according to the external standard method on peak areas or peak heights.

3 Results

3.1 UHPLC-HRMS Data Analysis and Chemometrics

Application of the untargeted LC-MS screening in the three different honey types, as described in previous work (Stavropoulou et al., 2021), resulted in a peak list containing the accurate mass vs. intensity of the detected features. Then, an untargeted metabolomics workflow encompassing statistical analysis was performed to identify the most important markers that discriminate Anama from pine and thyme honey samples.

More specifically, aiming to identify potential chemical markers of Anama honey when paralleled with pine and thyme honey, multivariate statistical analysis (MVA) i.e., Principal Component Analysis (PCA), Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) and differential analysis were performed using the entire metabolic profile of the samples (all features obtained). QC clustering in PCA analysis in both modes confirmed the reliability of the system and the quality of the acquired data. Based on the results of PCA analysis, a tendency of separation of the three different honey types could be observed in both modes (Figure 1, PCA and OPLS-DA loading plots are presented in Supplementary Figures S1–S6). These findings indicate that the majority of their chemical content was different for the three honey types, whereas some metabolites were similar. Two representative chromatograms of the three types of honey are presented in Figure 2.

FIGURE 1
www.frontiersin.org

FIGURE 1. PCA score plots (A,B) of the honey Anama (orange), pine (green), thyme (blue) and QC (purple) samples after UHPLC-HRMS analyses in negative (R2X = 0.527) (A) and positive (R2X = 0.524) (B) ion mode, respectively.

FIGURE 2
www.frontiersin.org

FIGURE 2. UHPLC-HRMS representative chromatograms of Anama, thyme, and pine (from up to bottom respectively) honey extracts in negative [(A), ESI-] and positive [(B), ESI+] ionization mode.

To further verify the differences in metabolic profiles, OPLS-DA was carried out which showed a clear separation between Anama and other monofloral honeys (Figures 3, 4). Features with Variable Importance in Projection (VIP) scoring from OPLS-DA models greater than 1.5 combined with adjusted p-value ≤0.05 and Log2 Fold Change (Anama vs. Thyme and Anama vs. Pine) >1.5 from differential analysis, were considered as contributing the most to group separation.

FIGURE 3
www.frontiersin.org

FIGURE 3. PCA (Ai, Bi), OPLS-DA (Aii, Bii) score plots and permutation tests (Aiii, Biii) after the comparison of Anama (orange) and Thyme (blue) honeys employing UHPLC-HRMS analysis in negative (A) and positive (B) ion mode, respectively. The goodness of fit and prediction of the models were for PCA analysis (Ai) R2X = 0.504, (Bi) R2X = 0.586 and for OPLS-DA analysis (Aii) R2X = 0.591, R2Y = 0.963, Q2 = 0.862, (Bii) R2X = 0.52, R2Y = 0.909, Q2 = 0.831, x and y axes of panels (ii) correspond to T-score (t1) and Orthogonal T-score (to1) respectively.

FIGURE 4
www.frontiersin.org

FIGURE 4. PCA (Ai, Bi), OPLS-DA (Aii, Bii) score plots and permutation tests (Aiii, Biii) after the comparison of Anama (orange) and pine (green) honeys employing UHPLC-HRMS analysis in negative (A) and positive (B) ion mode, respectively. The goodness of fit and prediction of the models were for PCA analysis (Ai) R2X = 0.653, (Bi) R2X = 0.618 and for OPLS-DA analysis (Aii) R2X = 0.644, R2Y = 0.908, Q2 = 0.828, (Bii) R2X = 0.707, R2Y = 0.991, Q2 = 0.849, x and y axes of panels (ii) correspond to T-score (t1) and Orthogonal T-score (to1) respectively.

Based on the aforementioned statistical criteria, 32 characteristic compounds were found to be increased in all Anama samples analyzed, when compared to pine and to thyme honeys (using the obtained peak area for relative quantification). Similarly, according to the statistical analysis and the log2FC values, various compounds were found to be differentially increased in thyme and pine honey compared to Anama. Some key bioactive molecules putatively annotated in Anama honey that can serve as specific chemical markers of Anama honey, including their experimental monoisotopic mass, retention time (tR), MS/MS fragment ions, and molecular formula, are presented in Table 1. For the compounds that the putative annotation was based on in silico fragmentation, the final score is presented in Supplementary Table S1. Figure 5 presents some representative chemical structures of characteristics markers in Anama honey. (±)-Abscisic acid (Table 1) was identified after a comparison of MS/MS data with data produced from mzCloud database (Figure 6) and authentic analytical standard solution injection.

TABLE 1
www.frontiersin.org

TABLE 1. Monoisotopic mass and MS/MS fragment ions of selected secondary metabolites and compounds, potential chemical markers of Anama honey.

FIGURE 5
www.frontiersin.org

FIGURE 5. Chemical structures of selected differentially increased compounds in Anama honey.

FIGURE 6
www.frontiersin.org

FIGURE 6. Comparison of MS/MS data from Anama honey (top spectra) to MS/MS of (±)-abscisic acid (A) and isophorone (B) of mzCloud database (bottom spectra).

3.2 Melissopalynological Analysis

The results of melissopalynological analysis showed that all Anama honeys were compliant with the requirements of the National legislation of Greece (Author Anonymous, 2022) as well as with the regulations of other European countries about heather honey (Thrasyvoulou et al., 2018), having a percentage of Erica sp. pollen above 45% (Table 2, Supplementary Table S4). National regulations are important due to the lack of International or European standards for monofloral honeys. Indicative photos with pollen grains are presented in Figures 7A–D. Although Erica pollen dominates (av. 85%, range 63–99%, Supplementary Table S4), a large number of other species were detected (av. 24 taxa per sample, range 15–39, Table S4), characteristic of the rich flora of the island. The relative proportion of the rest of the taxa in honey is low, a fact that is natural due to the overabundance of pollen grains from Erica. However, their presence is crucial for the identification of the geographical origin of honey, as it is a good fingerprint of the environment where the honey comes from (von der Ohe et al., 2004; Louveaux et al., 1978). Alongside, melissopalynological analysis did not detect any exogenous elements resulting from bee feed, e.g., starch granules, demonstrating the authenticity of honey.

TABLE 2
www.frontiersin.org

TABLE 2. Pollen spectrum of honeys studied. Data presented as frequency classes.

FIGURE 7
www.frontiersin.org

FIGURE 7. Indicative pollen grains were identified in Anama honey. (A) Erica-Thymus, (B) Erica-Arbutus-Myrtus, (C) Erica-Crocus-Pistacia, (D) Erica-Ceratonia.

3.3 Physicochemical Analysis

The results of physicochemical analysis are shown in Table 3. All parameters were compliant with the requirements of European and International legislation (Council-Directive, 2001; CXS-12-1981, 1981), provided that heather honey “Αnama” is considered “honey with a low natural enzyme content”. Based on international regulations the diastase activity of honey is in general not less than 8 Schade units, and in the case of honeys with a low natural enzyme content not less than 3 Schade Units (CXS-12-1981, 1981). In the current study, Anama honey had a Diastase number of m.v. = 7.3 Schade units (SD = 3.2, range = 2.0–11.5 Schade units). All samples were freshly collected and the HMF value, checked for some of them, was always <15 mg/kg supporting the above categorization of Anama honey (Council-Directive, 2001). The low diastase content is not surprising, as it is also mentioned in the literature for Erica sp. honeys (Oddo et al., 2004; Thrasyvoulou et al., 2018).

TABLE 3
www.frontiersin.org

TABLE 3. Physicochemical parameters of Anama honeys.

As regards the electrical conductivity of heather honey, it is considered an exception by International and Greek legislation and therefore there is no legal limit (Council-Directive, 2001; Oddo et al., 2004; Thrasyvoulou et al., 2018; Author Anonymous, 2022; CXS-12-1981, 1981; Thrasyvoulou and Manikis, 1995). Nevertheless, the electrical conductivity determined for “Anama” honey in the current study, is consistent with the values given by the International literature (e.g. m.v. = 0.67, SD = 0.16 mS/cm, 20°C (Council-Directive, 2001; Oddo et al., 2004; Thrasyvoulou et al., 2018; CXS-12-1981, 1981; Thrasyvoulou and Manikis, 1995), Table 3: m.v. = 0.67, SD = 0.20 mS/cm, 20°C).

Among the rest of physicochemical parameters, the relatively high percentage of moisture is worthy of note: m.v. = 18.0%w/w, SD = 0,7 (Table 3), being in agreement with International Legislation, as the upper legal limit of heather honey is, by way of exception 23%, instead of 20% for the rest of honeys [48. 49]. Sugar analysis was also conducted. The values determined for Fructose, Glucose, Saccharose, and Maltose were within the legal limits for blossom honey (Council-Directive, 2001; CXS-12-1981, 1981). The color of Anama honey is initially reddish-bronze with dark green tones (m.v. = 105 mm Pfund), but it soon turns into light caramel color, after a fine crystallization occurs, 1–3 months after harvesting. Finally, the determined pH values ranged from 4.0 to 4.5 (Table 3).

4 Discussion

4.1 Anama Honey Secondary Metabolites—Chemical Markers

The UHPLC-HRMS profile of the phenolic compounds of the three different monofloral honeys exhibited some common compounds including p-coumaric acid, caffeic acid, syringic acid, and flavones such as apigenin and chrysin and the flavanones, hesperetin, pinocembrin, taxifolin and naringenin (Table 4). These phenolic compounds are also common in various types of honeys as reported in the literature. More specifically, syringic acid and chrysin were previously reported as thyme, pine, and fir Greek honey components (Karabagias et al., 2014), displaying substantial fluctuation in concentrations among these types of honeys. In the presented study, though quantitation was not performed, considering the relative abundances obtained, such a marked concentration increase for one particular honey was not observed. Chrysin and pinocembrin were also key components of Portuguese heather honey (Ferreres et al., 1994). Apigenin, p-coumaric acid, and taxifolin were additionally reported as Greek heather and thyme honey constituents, demonstrating differences in concentrations (Koulis et al., 2021). Caffeic acid and hesperetin identified in the presented work, are also constituents of several honey types, including heather honey (Sergiel et al., 2014). Generally, heather honeys from Erica are characterized by high amounts of abscisic acid (Ferreres et al., 1996; da Silva et al., 2021) a fact also confirmed in this study. Quantitation of abscisic acid due to its indisputable utility as a chemical marker and the availability of analytical standard was performed in all Anama honey samples by LC-MS, resulting in a mean concentration of 104.11 ± 7.15 μg/g honey, (see Supplementary Figure S7 and Supplementary Materials S7, S8). Overall, this compound was reported to be used for floral authentication, but the quantitative analysis is crucial since abscisic acid is a constituent of another type of honeys at lower concentrations reaching the ng/g scale (Wang et al., 1998).

TABLE 4
www.frontiersin.org

TABLE 4. Phenolic compounds common in Anama, Thyme and Pine honeys.

The results of untargeted HRMS metabolomic analysis and subsequent workflow unveiled 32 compounds as potential chemical markers of Anama honey, most of which have never been reported as Greek heather honey components (Table 1). One of the differentially increased molecules in “Anama” honey was annotated as the compound picein (Table 1). Picein is a phenolic compound, glucoside of piceol, reported as a constituent of several plants. Among them, Salvia (Lamiaceae) (Wang et al., 1998), Opuntia ficus-indica (Cactaceae) (Strid, 2016), and Helianthemum (Cistaceae) (Chemam et al., 2017), are plants of high beekeeping value, widespread in Ikaria and Fournoi islands (Table 2). For the mass feature with m/z 332.1101, tR = 8.69 min, the molecular formula of C14H20O9 was assigned to leonuriside A (Table 1). Leonuriside A is a phenolic glycoside reported to be encountered in plants of Paronychia genus (native to Ikaria Island), one of the genera of the family Caryophyllaceae, while pollen grains of these nectarless plants have been identified in Anama honey samples in this work (Table 2). Astragalin, a kaempferol O-glucoside, detected at tR = 10.5 min with m/z 448.1007 (Table 1) has been also reported as a Polish heather honey constituent (Goslinski et al., 2021), while it is its first report in Greek heather honey. Kaempferol-3-Ο-galactoside (monoisotopic mass at m/z 448.1007, molecular formula C21H20O11, Table 1) or trifolin, is a kaempferol derivative identified in Anama honey. A respective glucoside was detected in carob fruit (Ceratonia siliqua L.) (Kyriacou et al., 2021). Pollen of Ceratonia is identified in nine out of ten samples of Anama honey, with a relatively high pollen frequency (Table 2). Salicylic acid, a phenolic acid detected at tR10.80 with m/z 138.0306), has been also reported in a previous study using gas chromatography–mass spectrometry, as a component of Greek thyme honey at relatively low concentrations (Alissandrakis et al., 2009). 2-anisic acid (2-methoxy-benzoic acid) detected at tR = 10.81 with m/z 152.0474 (Table 1), has been previously reported as a volatile marker for Greek cotton honey (Alissandrakis et al., 2005).

In the frame of rarely reported compounds, domesticoside (tR = 8.88, with m/z 344.1109, Table 1), a phenolic glycoside, was putatively annotated for the first time in Anama honey. This molecule is a component of several plants and among them Artemisia and Prunus. To our knowledge, two species of the genus Artemisia are distributed in Ikaria: Artemisia arborescens and Artemisia vulgaris. With regard to Prunus, there are several Prunus species, mostly cultivated, and abundant in Ikaria Island. Catalpol (tR = 8.63, with m/z 362.1205, Table 1) is an iridoid glycoside, of the class of monoterpenoids, described as a nectar secondary compound (Michaud et al., 2019). Another iridoid glycoside aucubin (tR = 11.04, with m/z 346.1258, Table 1), along with catalpol, are pharmacologically active natural components of numerous plants (Sertic et al., 2015; Miehe-Steier et al., 2015). Ganolucidic acid B (tR = 17.73, with m/z 502.3296, Table 1), belongs also to the terpenoids family and it has been isolated from Ganoderma Lucidum fungus. Its distinctive increase in Anama honey is worth mentioning considering its first report in honey, and the pronounced bioactivity of Ganoderma Lucidum which is conceived as traditional medicine (Yang Y. et al., 2019). Methyl 3-aminobenzoate (M3AB) (tR = 6.67, with m/z 151.0634, Table 1) belongs to aminobenzoates, while its isomer methyl 2-aminobenzoate or methyl anthralinate is established as a chemical marker of citrus honey (Sesta et al., 2008). Consequently, M3AB can be reported as a potential chemical marker of Anama honey. Benzoate derivatives are also key synthons for the development of bioactive molecules (Kronenberger et al., 2020).

Regarding alkyl phenols, 4-propylphenol (tR = 11.02, with m/z 136.0889, Table 1) emerged as a putative chemical marker of Anama honey. The literature search showed that the methyl ether of 2-propylphenol has been a component of honeydew honey (Štefániková et al., 2021). 3-Propylphenol and 4-propylphenol have been identified as aroma compounds of fermented cocoa (Füllemann and Steinhaus, 2020), the latter being an inhibitor of cyclic AMP phosphodiesterase (Nikaido et al., 1981). Another alkyl phenol, 2-ethylphenol (tR = 7.98, with m/z 122.0732, Supplementary Table S1), was also increased in Anama samples. Ethylphenols are reported as aroma compounds in wine (Monje et al., 2002).

Amino acids, including tryptophan, are also described as components of honey (Leu et al., 2012). D-tryptophan (tR = 10.81, with m/z 204.0899, Table 1) was detected in the Anama honey samples, whereas it has already been proposed as a discriminating marker of Acacia honey (Pirini et al., 1992). Tryptophan is pivotal in protein-protein interactions and its metabolite, kynurenic acid, displays significant bioactivity, being assigned as a putative biomarker of chestnut honey (Turski et al., 2016).

Benzaldehyde and its derivatives are well acquainted with volatile components of honey (Karabagias et al., 2019). 4-ethylbenzaldehyde (tR = 7.73, with m/z 134.0731, Table 1) can be utilized as a discriminating chemical compound of Anama. Benzaldehyde has been identified as a Greek heather honey constituent using GC-MS (Xagoraris et al., 2021), while its ethylated derivative has been detected for the first time in this work in Anama heather honey samples by LC-HRMS. Cuminaldehyde (tR = 16.43, with m/z 148.0889, Table 1), another benzaldehyde derivative, was differentially increased in Anama honey. This aldehyde is reported as a constituent of Croatian lime tree honey (Lusic et al., 2007) and Jarrah honey, another important antimicrobial commercial honey, originating from remote parts of Western Australia (Anand et al., 2019). Ferulic acid (tR = 8.90, with m/z 194.0580, Table 1) and its derivative, 5-hydroxyferulic acid methyl ester (tR = 11.21, with m/z 210.0524, Table 1) were also differentially increased in Anama honey (Table 1).

6-Methoxyluteolin 7-rhamnoside (tR = 11.54, with m/z 462.1163, Table 1) is a flavonoid reported in Brazilian red propolis (Silva-Beltran et al., 2021; do Nascimento et al., 2019), known especially for its biological potential (dos Santos et al., 2021). Isophorone (tR = 7.74, with m/z 138.1045, Table 1), a norisoprenoid compound, was identified in Anama honey samples while its MS spectrum and MS/MS fragments are presented in Figure 6. Norisoprenoid compounds, usually identified by means of gas chromatography, have been reported as constituents-markers of heather honey from Greece (Xagoraris et al., 2021) and other countries (Castro-Vázquez et al., 2012); hence, this work comes to verify previous findings. Isophorone has also been reported as a constituent of essential oils of Eruca vesicaria, a plant found in Aegean Islands (Hichri et al., 2019) and in Crocus species from Greece (Lamari et al., 2018). Crocus is widespread in Ikaria island and blooms simultaneously with Erica manipuliflora, during autumn and winter (Table 2).

The presented metabolomics approach also disclosed (+)-7-isomethyl jasmonate (isomethyl ester of jasmonic acid) as a characteristic compound that was found to increase in Anama honeys. Another key ester of jasmonic acid, methyl jasmonate (MJ), is a volatile compound known for its utility in plant defense mechanisms, and as an elicitor to improve the phenol content of fruits. It is found in higher plants, so its detection is expected. Jasmonates have long been studied for their anticancer activity. MJ, the most potent of this family, exhibits a specific cell death-induction effect on numerous cancer cells (Raviv et al., 2013). Vitamins are also reported as honey components. Riboflavin (tR = 9.27, with m/z 376.1384, Table 1) and its degradation product, lumichrome (tR = 10.81, with m/z 242.0802, Table 1) were differentially increased in Anama honey compared to thyme and pine honey of the same region. Both compounds have been reported as characteristic molecules of Dalmatian sage honey (Tuberoso et al., 2012). But, in principle, this finding is in agreement with the results regarding Turkish heather honey (Kıvrak, 2016), where vitamins of the B-group were predominant, and Greek heather honey (Koulis et al., 2021), in which lumichrome was abundant. An additional compound worth mentioning is dehypoxanthine futalosine (tR = 7.08, with m/z 296.0897, Table 1). It is a biosynthetic intermediate in the menaquinone (Vitamin K) biosynthesis (Manion-Sommerhalter et al., 2021), first reported as a potential honey constituent.

Phenylacetic acid (tR = 8.48, monoisotopic mass at m/z 166.063, Table 1) is a hydroxy monocarboxylic acid, hitherto stated as a potential marker of thistle honey (Tuberoso et al., 2011), found also in Manuka and Kanuka nectars (Bong et al., 2018), and typical for cornflower honey (Kus et al., 2014). The results of this work come to verify its role as a “marker substance” due to its accumulation in Anama honey (Koulis et al., 2021). 3,4-Dimethoxycinnamic acid (tR = 11.02, monoisotopic mass at m/z 208.0734, molecular formula C11H12O4), which was increased in Anama, is a molecule described as a marker of black locust (Robinia pseudoacacia) honey (Margaoan et al., 2021), also found in Greek cotton honey (Alissandrakis et al., 2005). The current work enlarges its utility as a non-specific chemical marker. Another constituent tentatively annotated in Anama honey is trichocarpin (tR = 12.4, monoisotopic mass at m/z 406.1277, Table 1). This molecule is referenced as a volatile metabolite of balsam poplar (a type of poplar tree) (Mattes et al., 1987), but is not yet reported as a honey component (based on research on the accessible bibliography). It is noteworthy that certain poplar trees are a pivotal source of propolis for bees. In the same context, plastoquinol annotated in Anama honey (Table 1) is the reduced form of plastoquinone (present in photosynthetic pathways) and a dominant antioxidant (Borisova-Mubarakshina et al., 2019).

The results of this study indicate that the excess of nectariferous (n = 44) and nectarless (n = 23) plants detected in the Anama honey can “guarantee” to a large extent the natural foundation of several of the phytochemicals detected. In all cases, blooming-flowering period information is essential when trying to connect the findings with the botanical origin and bees’ foraging. With all these in view, it can be inferred that there is still room for additional phytochemical analysis on a plethora of species that are either limitedly, scarcely, or not investigated at all.

The latter is a limitation of the presented work in the effort to associate the chemicals identified in honey with specific plants. This will be overcome by conducting a systematic study of the flora of Ikaria and simultaneous chemical analysis of the dominant plants of these islands. With these in view, heather honey from other regions of Greece can also be collected and analyzed along with Anama honey to unveil any potential differences as regards chemical markers and melissopalynological analysis.

Overall, Anama honey possesses a distinctive chemical profile that discriminates it from other important honeys of the Aegean region, an area with important biodiversity of flora and fauna (Archipelagos, 2022). This is supported by the fact that an increased relative abundance of more than 32 bioactive compounds was found in Anama honey compared to pine and thyme honey from the same geographical origin. From a chromatography point of view, it is evident that UPLC-HRMS has the capacity to detect compounds that are usually identified by gas chromatography, such as benzaldehydes and norisoprenoids. It is also apparent that the following procedure managed to extract a diversity of compounds covering a wide range of polarities. This is the first holistic report on the phytochemical profile of this exceptional Greek honey, embracing cutting-edge high-resolution mass spectrometry and statistical analysis reinforced by fundamental melissopalynological analysis. Based on reports from the literature, the majority of the annotated compounds that belong to various chemical classes have exhibited antimicrobial, antioxidant, and antiviral properties (Margaoan et al., 2021). Specifically, leonuriside A has been reported to inhibit the growth of human breast cancer MCF-7 cell line (Thai et al., 2017), also acting as a potential anti-diabetic agent (Nguyen et al., 2018). Trifolin is a bioactive molecule that induces apoptosis in human lung cancer cells (Kim et al., 2016). Catalpol’s tentative characterization and subsequent calculation of the differential increase in Anama honey deserves particular attention. Catalpol demonstrates therapeutic potential for central nervous system diseases, exemplified by Alzheimer and Parkinson (Zhang et al., 2021). Its activity is attributed to the inhibition of inflammation, antioxidant stress, improvement of mitochondrial dysfunction, and synaptic damage. An additional bioactive iridoid annotated in this work is aucubin. Its beneficial character is confirmed by various studies, exerting anti-nephrotoxic (Potacnjak et al., 2020; Potocnjak et al., 2021), anti-inflammatory (Olajide et al., 2020), and anti-osteoporotic effects (Li et al., 2020). Finally, ferulic acid, a common component in several types of honeys, and its analogues possess antioxidant properties and are tested in therapeutic schemes against several diseases (Paiva et al., 2013). Consequently, the newly reported molecules and the rest of the potential chemical markers can set the basis for further investigation of the biological activity of Anama honey.

4.2 Pine and Thyme Honey Secondary Metabolites

The results of the increased secondary metabolites and compounds found in Thyme and Pine honeys employing UHPLC-HRMS analyses in comparison with Anama honey (Anama vs. Thyme, Anama vs. Pine) are presented in Supplementary Table S2.

4.2.1 Thyme

Citrusin C (Supplementary Table S2) is a phenolic glycoside putatively annotated for the first time in thyme honey. It is found in several plants, such as Ducrosia (aerial parts) (Morgan et al., 2014), and Ononis (aerial parts) (Mezrag et al., 2017). Other citrusins (A, B, D) are reported constituents of several plants such as Artemisia (Xiao et al., 2021), Cortex (Yang S. et al., 2019), Buxus (Saleem et al., 2019), and sweet oranges (Díaz et al., 2014). Citrusins exhibit anti-inflammatory (Yang S. et al., 2019) and antiviral activity (Tan et al., 2012).

Another bioactive molecule putatively annotated in thyme honey is syringin (Supplementary Table S2). It is a phenolic glycoside found in a multitude of foods and plants such as oats (Pretorius et al., 2021) and Daphne (Frezza et al., 2021; Ghanadian et al., 2020). It is a bioactive molecule displaying cardioprotective action in synergies with other compounds (Yao et al., 2021). (±)-Ribaline is a rare quinoline alkaloid (Supplementary Table S2), referenced as a constituent of Zanthoxylum mayu (Rutaceae family) (Barrera-Thomas, 2015). Gardenoside is a bioactive iridoid originating from Gardenia fruits. Norcoclaurine (Supplementary Table S2) is a 1-benzylisoquinoline that has displayed anti-HIV activity (Kashiwada et al., 2005), therefore its possible annotation increases the bioactive profile of any natural product; in this case, thyme honey. Its (S)-enantiomer displays partial neuroprotective activity (Peana et al., 2019).

4.2.2 Pine

The LC-MS analysis of pine honey sample extracts also displayed several bioactive molecules, some of them first reported in this type of honey. Piceid (Supplementary Table S2), amongst others, is acknowledged as the major resveratrol derivative in grape juices (Romero-Pérez et al., 1999). Piceid and resveratrol are recognized for their anti-aging properties (Shen et al., 2017), hence their annotation in natural products is an added value. 10-epi-Eupatoroxin (Supplementary Table S2) is a sesquiterpene lactone whose cytotoxic activity is promising, as demonstrated by a quantitative structure-activity relationship (QSAR) study (Scotti et al., 2007). 3-methylsuberic acid is its first putative annotated report in pine honey. Suberic acid and derivatives are found in plants such as Isatis tinctoria L. (Hartleb and Seifert, 1995), Viola odorata (Jasim et al., 2018), and Dalbergia odorifera (Sun et al., 2021). Interestingly, in pine honey, the biologically important 10-hyrodxydecenoic acid (royal jelly acid) was also differentially increased. A summary of the current findings for thyme and pine honey is being prepared to present the full results and possible correlation with plant sources.

4.3 Melissopalynological Analysis of Anama Honey

Despite the predominance of Erica pollen in all samples (av. 85%, range 63–99%, Supplementary Table S4), the study of the whole pollen spectra generates important information about the natural habitat of Ikaria and Fournoi and reflects the characteristics of the biome, which possesses a rich flora. As expected, the main contribution to nectar, apart from Erica manipuliflora, comes from plants that flower in the same season as heather, i.e., autumn and winter, especially madrones (Arbutus), carob trees (Ceratonia), asparagus, false yellowhead (Dittrichia), common ivy (Hedera), etc. (Table 2). In contrast, nectarless plants are dominated by those that flower in spring, indicating that they come from pollen stored in honeycombs. Thus, the families of Anacardiaceae, Cistaceae, Oleaceae, Scrophulariaceae, Fagaceae, and Fabaceae are dominant, although these families include spring and summer flowering plants. The only exception is the autumn-blooming Cyclamen sp. pollen, which is present in the 70% of the honey samples.

Nectariferous plants from previous blooming periods were also detected. A possible explanation is that they come from frames bearing honey and beebread from a spring or summer foraging, that has not been extracted (e.g., due to the presence of brood). This explains the high rates of plants from the Lamiaceae family that are observed in some cases, since these islands also produce excellent quality thyme honey. Pollen of Myrtus and Styrax is also distinguished, due to the wide distribution of these plants on the island. (Table 2, Supplementary Table S4). Pollen grains per 10 g of Anama honey, ranged from 50,000 to 210,000 (m.v. 103,314, SD = 45,310) and taking into account only nectariferous plants PG Nect./10 g ranged from 30,000 to 170,000 (m.v. 127,387, SD = 46,460) (Table 2, Supplementary Table S4). In addition to pollen grains, HDE were observed and counted. According to Louveaux et al. (Louveaux et al., 1978), the extent to which a given honey sample is derived from different plant sources can be deduced from the frequencies of the pollens and honeydew elements in it. This rule is valid only if the honey contains few honeydew elements (HDE/p < 1). In our study, the ratio of honeydew elements to pollen grains was determined and found to be HDE/p < 1 in all cases (Table 2, Supplementary Table S4).

Therefore, apart from proving the main botanical origin of honey, the melissosopalynological analysis of honey “Anama” can provide information about its geographical origin, as it is directly related to the flora of Ikaria and Fournoi. Both nectariferous and nectarless plants contribute to the Melissopalynological profile of the honey studied.

4.4 Melissopalynological Analysis of Thyme and Pine Honey

4.4.1 Thyme Honey

Thyme honey originating from Kea, Syros, and Fournoi islands, is dominated by pollen of Thymbra capitata (m.v. = 65%, range = 18%–97%), accompanied by pollen of the native flora of the Mediterranean maquis.

The most important families of nectariferous plants, observed in more than half of the samples, are reported here, along with some representative plant genera (in descending order of occurrence in honey): Lamiaceae (Thymbra/Thymus, Salvia, Lavandula, Teucrium, Ballota, Phlomis/Lamium, Satureja), Asteraceae (Taraxacum/Crepis, Centaurea spinosa, Dittrichia/Inula, Anthemis), Fabaceae (Melilotus, Anthyllis hermanniae, Trifolium, Ononis, Acacia, Lotus, Vicia), Brassicaceae (Sinapis, Eruca), Apiaceae (Daucus/Crithmum//Foeniculum, Smyrnium, Ferula, Tordyllium, Pimpinella), Boraginaceae (Echium, Cynoglossum, Heliotropium), Myrtaceae (Eucalyptus, Myrtus), Convolvulaceae (Convolvulus), Asphodelaceae (Asphodelus), Malvaceae (Malva), Aizoaceae, Oxalidaceae (Oxalis), Tamaricaceae [Tamarix (SYROS)], Rosaceae (KEA), Ericaceae [Erica (in specific cases as secondary pollen in FOURNOI thyme honey)], Aceraceae (Acer, in specific cases as secondary pollen, only in KEA thyme honey)].

As a result of quantitative analysis, the pollen grains per 10 g of thyme honey had a mean value of 29,000 (range = 3800–67,000). The number of honeydew elements was in all cases low or very low: HDE/p < 0.1 in 96% of samples.

4.4.2 Pine Honey

Although pine honey comes from the honeydew of Marchalina hellenica, an insect that lives in the Pinus brutia forests of Aegean islands, melissopalynological analysis revealed a plethora of plants represented by the pollen encountered in the honey sediment. The mean value of taxa per sample was 41, range = 26–55, reflecting the rich flora of Samos island. Pine honey originating from Samos island is dominated by pollen of forest plants, especially deciduous trees, shrubs, climbing plants, and wildflowers. The most important families and plant genera, observed in more than half of the samples, are reported here, in descending order of occurrence in honey.

Families of nectariferous plants: Asteraceae, Fabaceae, Lamiaceae, Rosaceae, Ericaceae, Fagaceae, Boraginaceae, Asparagaceae, Apiaceae, Caesalpiniaceae, Brassicaceae, Araliaceae, Anacardiaceae, Myrtaceae, Smilacaceae, Lauraceae, Liliaceae, Styracaceae, Oxalidaceae.

Important Genera of nectariferous plants: Castanea, Dittrichia,Inula, Rhus, Myrtus, Melilotus, Thymbra,Thymus, Erica, Asparagus, Hedera, Ceratonia, Taraxacum, Rubus, Anthemis, Sinapis, Trifolium, Smilax, Styrax, Laurus, Onopordum, Oxalis.

Nectarless plants also play an important role in the melissopalynological profile of honeydew honey, especially considering that the honeydew itself does not provide pollen to bees. Therefore, families and genera of nectarless plants are listed below in descending order: Fagaceae (Quercus coccifera/ilex, Quercus pubescens), Cistaceae, Anacardiaceae (Pistacia), Oleaceae (Olea), Scrophulariaceae (Verbascum), Pinaceae, Hypericaceae (Hypericum), Fabaceae (Calicotome-Genista), Chenopodiaceae/Amaranthaceae, Ephedraceae (Ephedra), Platanaceae (Platanus), Rosaceae (Poterium), Ranunculaceae (Anemone).

As expected, the number of honeydew elements was numerous (HDE/p = 3.0–4.5) or very numerous (HDE/p > 4.5) in the 67% of the samples, and of medium quantity (HDE/p = 1.5–2.9) in the 30% of the samples. As a result of quantitative analysis, the pollen grains per 10 g of pine honey had a mean value of 71,000 (range = 13,000–240,000) and the HDE per 10 g of honey had a mean value of 245,000 (range = 100,000–570,000).

It is obvious that the pollen spectrum in heather, thyme, and pine honey differs significantly due to the different botanical origins and the different periods of the year in which the honey is harvested. Nevertheless, it is the comparative metabolomic analysis that revealed important metabolites and bioactive compounds for each type of honey.

Conclusion

The application of an untargeted UHPLC-HRMS screening with a metabolomics workflow assisted by statistical analysis, provided a diversity of compounds as tentative, specific, and non-specific, chemical markers of Ikaria’s unique heather honey called “Anama”. 32 chemicals (among them, phenolic compounds and triterpenes) emerged as potential chemical markers of Anama honey, discriminating it from the other important Greek honeys of pine and thyme. Some of these chemicals, such as aucubin, catalpol, leonuriside A, picein, domesticoside, 4-ethylbenzaldehyde, 3,4-dimethoxycinnamic acid, dehypoxanthine futalosine, ganolucidic acid B, kaempferol-3-O-galactoside, 2-ethylphenol, and 4-propylphenol were not previously reported in Greek heather honey, and these compounds could be useful markers of the botanical origin of Anama which has the potential to be regarded as honey with protected designation of origin. Melissopalynological analysis also identified an excess of nectariferous (n = 44) and nectarless (n = 23) plants in Anama honey, and an attempt was made, in addition to the accurate HRMS spectra, to link the metabolites identified with the findings of the melissopalynological analysis. The latter is essential in the context of connecting secondary metabolites and compounds with the flora of this Mediterranean blue zone region, reinforcing the dependence of these compounds on the botanical origin of Anama honey. The identified metabolites exhibit significant antioxidant activity, offering Anama honey exceptional characteristics with advantageous effects on human health. The future perspective and steps of this work should include the targeted UHPLC-HRMS quantitative analysis of all commercially available chemical markers tentatively identified in this work. In the same context, nectar and pollen from the dominant plants of the region, which were identified in the presented study by extensive melissopalynological analysis, should be collected and subsequent phytochemical analysis pursued to fully elucidate the origin of the chemical substances found in Anama honey. Finally, heather honey from other regions of Greece can also be investigated in parallel with Anama honey to disclose any probable differences with respect to chemical markers and melissopalynological analysis.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Author Contributions

“Conceptualization, KMK, SI, DT, and KM; methodology, KMK, EB, SI, KB, IT, EM-K, TB; software, KMK, EB, SI, KB, IT; validation, KMK, SI; formal analysis, KMK, EB, SI, IT; investigation, KMK, SI, EB, IT, EM-K, TB; resources, KM, DT; data curation, KMK, EB, SI, KB, IT; writing—original draft preparation, KMK; writing—review and editing, KMK, EB, SI, DT, KM; visualization, EB, SI, IT, KMK; supervision, KMK, KM, DT; project administration, KMK, DT, KM; funding acquisition, KM, DT. All authors have read and agreed to the published version of the manuscript.”

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

Authors would like to thank the beekeeping and rural cooperatives of five Aegean islands for the provision of honey samples.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2022.924881/full#supplementary-material

References

Alissandrakis, E., Kibaris, A. C., Tarantilis, P. A., Harizanis, P. C., and Polissiou, M. (2005). Flavour Compounds of Greek Cotton Honey. J. Sci. Food Agric. 85, 1444–1452. doi:10.1002/jsfa.2124

CrossRef Full Text | Google Scholar

Alissandrakis, E., Tarantilis, P. A., Pappas, C., Harizanis, P. C., and Polissiou, M. (2009). Ultrasound-assisted Extraction Gas Chromatography-Mass Spectrometry Analysis of Volatile Compounds in Unifloral Thyme Honey from Greece. Eur. Food Res. Technol. 229, 365–373. doi:10.1007/s00217-009-1046-8

CrossRef Full Text | Google Scholar

Anand, S., Deighton, M., Livanos, G., Pang, E. C. K., and Mantri, N. (2019). Agastache Honey Has Superior Antifungal Activity in Comparison with Important Commercial Honeys. Sci. Rep. 9, 18197. doi:10.1038/s41598-019-54679-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Archipelagos (2022). Available at: https://archipelago.gr/en/our-work/terrestrial-conservation/forests-flora/(accessed January 12, 2022).

Author Anonymous, (2022). Greek Code for Foodstuffs and Beverages, Article 67a. Available at: https://www.aade.gr/sites/default/files/2020-03/67a-iss1.pdf (Accessed on March 30, 2022).

Google Scholar

Barrera-Thomas, M. (2015). Isolation and Chemical Modifications of Natural Products From Medicinal Plants. Quebec: Université du Québec à Montréal.

Google Scholar

Behm, F., von der Ohe, K., and Henrich, W. (1996). Zuverlässigkeit der Pollenanalyse von Honig. Bestimmung der Pollenhäufigkeit. Dtsch. Leb. 92, 183

Google Scholar

Bong, J., Loomes, K. M., Lin, B., and Stephens, J. M. (2018). New Approach: Chemical and Fluorescence Profiling of NZ Honeys. Food Chem. 267, 355–367. doi:10.1016/j.foodchem.2017.07.065

PubMed Abstract | CrossRef Full Text | Google Scholar

Borisova-Mubarakshina, M. M., Vetoshkina, D. V., and Ivanov, B. N. (2019). Antioxidant and Signaling Functions of the Plastoquinone Pool in Higher Plants. Physiol. Plant 166, 181–198. doi:10.1111/ppl.12936

PubMed Abstract | CrossRef Full Text | Google Scholar

Braconi, D., Millucci, L., Parisi, M. L., Spiga, O., and Santucci, A. (2021). “Omics-based Technologies for Food Authentication and Traceability,” in Food Authentication and Traceability. Editor CM Galanakis (Cambridge, Massachusetts: Academic Press), 215–245. doi:10.1016/b978-0-12-821104-5.00003-9

CrossRef Full Text | Google Scholar

Castro-Vázquez, L., Alañon, M. E., Gonzalez-Viñas, M. A., and Pérez-Coello, M. S. (2012). Changes in the Volatile Fractions and Sensory Properties of Heather Honey during Storage under Different Temperatures. Eur. Food Res. Technol. 235, 185–193. doi:10.1007/s00217-012-1756-1

CrossRef Full Text | Google Scholar

Chantarawong, W., Kuncharoen, N., Tanasupawat, S., and Chanvorachote, P. (2019). Lumichrome Inhibits Human Lung Cancer Cell Growth and Induces Apoptosis via a P53-dependent Mechanism. Nutr. Cancer 71, 1390–1402. doi:10.1080/01635581.2019.1610183

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaudhari, S. S., and Gokhale, D. V. (2016). Phenyllactic Acid: A Potential Antimicrobial Compound in Lactic Acid Bacteria. J. Bacteriol. Mycol. Open Access 2, 121–125. doi:10.15406/jbmoa.2016.02.00037

CrossRef Full Text | Google Scholar

Chemam, Y., Benayache, S., Marchioni, E., Zhao, M., Mosset, P., and Benayache, F. (2017). On-Line Screening, Isolation and Identification of Antioxidant Compounds of Helianthemum Ruficomum. Molecules 22, 239. doi:10.3390/molecules22020239

PubMed Abstract | CrossRef Full Text | Google Scholar

Cianciosi, D., Forbes-Hernández, T. Y., Afrin, S., Gasparrini, M., Reboredo-Rodriguez, P., Manna, P. P., et al. (2018). Phenolic Compounds in Honey and Their Associated Health Benefits: A Review. Molecules 23, 2322. doi:10.3390/molecules23092322

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciulu, M., Spano, N., Pilo, M., and Sanna, G. (2016). Recent Advances in the Analysis of Phenolic Compounds in Unifloral Honeys. Molecules 21, 451. doi:10.3390/molecules21040451

PubMed Abstract | CrossRef Full Text | Google Scholar

Cook, P. W., and Nightingale, K. K. (2018). Use of Omics Methods for the Advancement of Food Quality and Food Safety. Anim. Front. 8, 33–41. doi:10.1093/af/vfy024

PubMed Abstract | CrossRef Full Text | Google Scholar

Council-Directive (2001). COUNCIL DIRECTIVE 2001/110/EC of 20 December 2001 Relating to Honey.

Google Scholar

CXS-12-1981 (1981). Codex Alimentarius Standard for Honey.

Google Scholar

da Silva, L. R., Chiste, R. C., and Fernandes, E. (2021). Chemical and Antioxidant Characterization of the Portuguese Heather Honey from Calluna Vulgaris. Separations 8, 177. doi:10.3390/separations8100177

CrossRef Full Text | Google Scholar

Dawood, H. M., Shawky, E., Hammoda, H. M., Metwally, A. M., and Ibrahim, R. S. (2020). Chemical Constituents from Artemisia Annua and Vitex Agnus-Castus as New Aromatase Inhibitors: In-Vvitro and In-Silico Studies. J. Mex. Chem. Soc. 64, 316–326. doi:10.29356/jmcs.v64i4.1236

CrossRef Full Text | Google Scholar

Di Marco, G., Gismondi, A., Panzanella, L., Canuti, L., Impei, S., Leonardi, D., et al. (2018). Botanical Influence on Phenolic Profile and Antioxidant Level of Italian Honeys. J. Food Sci. Technol. 55, 4042–4050. doi:10.1007/s13197-018-3330-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Díaz, R., Pozo, O. J., Sancho, J. V., and Hernández, F. (2014). Metabolomic Approaches for Orange Origin Discrimination by Ultra-high Performance Liquid Chromatography Coupled to Quadrupole Time-Of-Flight Mass Spectrometry. Food Chem. 157, 84–93. doi:10.1016/j.foodchem.2014.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

do Nascimento, T. G., Arruda, R. E. D. S., Almeida, E. T. d. C., Oliveira, J. M. D. S., Basílio-Júnior, I. D., Porto, I. C. C. d. M., et al. (2019). Comprehensive Multivariate Correlations between Climatic Effect, Metabolite-Profile, Antioxidant Capacity and Antibacterial Activity of Brazilian Red Propolis Metabolites during Seasonal Study. Sci. Rep. 9, 18293. doi:10.1038/s41598-019-54591-3

PubMed Abstract | CrossRef Full Text | Google Scholar

dos Santos, C. M., Mesquita, L. M. d. S., Braga, A. R. C., and de Rosso, V. V. (2021). Red Propolis as a Source of Antimicrobial Phytochemicals: Extraction Using High-Performance Alternative Solvents. Front. Microbiol. 12, 659911. doi:10.3389/fmicb.2021.659911

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreres, F., Andrade, P., and Tomás-Barberán, F. A. (1996). Natural Occurrence of Abscisic Acid in Heather Honey and Floral Nectar. J. Agric. Food Chem. 44, 2053–2056. doi:10.1021/jf9507553

CrossRef Full Text | Google Scholar

Ferreres, F., Andrade, P., and Tomasbarberan, F. A. (1994). Flavonoids from Portuguese Heather Honey. Z Leb. Unters Forch 199, 32–37. doi:10.1007/bf01192949

CrossRef Full Text | Google Scholar

Frezza, C., Venditti, A., De Vita, D., Sciubba, F., Tomai, P., Franceschin, M., et al. (2021). Phytochemical Analysis and Biological Activities of the Ethanolic Extract of Daphne Sericea Vahl Flowering Aerial Parts Collected in Central Italy. Biomolecules 11, 379. doi:10.3390/biom11030379

PubMed Abstract | CrossRef Full Text | Google Scholar

Füllemann, D., and Steinhaus, M. (2020). Characterization of Odorants Causing Smoky Off-Flavors in Cocoa. J. Agric. Food Chem. 68, 10833–10841. doi:10.1021/acs.jafc.0c04633

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghanadian, M., Ali, Z., Khan, I. A., Balachandran, P., Nikahd, M., Aghaei, M., et al. (2020). A New Sesquiterpenoid from the Shoots of Iranian Daphne Mucronata Royle with Selective Inhibition of STAT3 and Smad3/4 Cancer-Related Signaling Pathways. DARU J. Pharm. Sci. 28, 253–262. doi:10.1007/s40199-020-00336-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Goslinski, M., Nowak, D., and Szwengiel, A. (2021). Multidimensional Comparative Analysis of Bioactive Phenolic Compounds of Honeys of Various Origin. Antioxidants-Basel. 10, 530. doi:10.3390/antiox10040530

PubMed Abstract | CrossRef Full Text | Google Scholar

Gourdomichali, T., and Papakonstantinou, E. (2018). Short-term Effects of Six Greek Honey Varieties on Glycemic Response: A Randomized Clinical Trial in Healthy Subjects. Eur. J. Clin. Nutr. 72, 1709–1716. Epub 2018/04/25. doi:10.1038/s41430-018-0160-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y. (2021). Metabolomics of Mature Honey Formation. Nat. Food 2, 223. doi:10.1038/s43016-021-00273-1

CrossRef Full Text | Google Scholar

Hartleb, I., and Seifert, K. (1995). Acid Constituents fromIsatis Tinctoria. Planta Med. 61, 95–96. doi:10.1055/s-2006-958019

PubMed Abstract | CrossRef Full Text | Google Scholar

Heitz, T., Smirnova, E., Widemann, E., Aubert, Y., Pinot, F., and Ménard, R. (2016). The Rise and Fall of Jasmonate Biological Activities. Subcell. Biochem. 86, 405–426. Epub 2016/03/31. doi:10.1007/978-3-319-25979-6_16

PubMed Abstract | CrossRef Full Text | Google Scholar

Hichri, F., Hichri, A. H., Maha, M., Hossan, A. S. M., Flamini, G., and Jannet, H. B. (2019). Chemical Composition, Antibacterial, Antioxidant and In Vitro Antidiabetic Activities of Essential Oils from Eruca Vesicaria. Chem. Biodivers. 16, e1900183. doi:10.1002/cbdv.201900183

PubMed Abstract | CrossRef Full Text | Google Scholar

Hills, S. P., Mitchell, P., Wells, C., and Russell, M. (2019). Honey Supplementation and Exercise: A Systematic Review. Nutrients 11, 1586. doi:10.3390/nu11071586

PubMed Abstract | CrossRef Full Text | Google Scholar

IHC (2009). Harmonised Methods of the International Honey Commission 2009. Available at: https://www.ihc-platform.net/ihcmethods2009.pdf (accessed on 20 April 2021)

Google Scholar

Jasim, S. F., Baqer, N. N., and Alraheem, E. A. (2018). Detection of Phytochemical Constituent in Flowers of Viola Odorata by Gas Chromatography-Mass Spectrometry. Asian J. Pharm. Clin. Res. 11, 262–269. doi:10.22159/ajpcr.2018.v11i5.24288

CrossRef Full Text | Google Scholar

Kamaruzzaman, M. A., Chin, K. Y., and Ramli, E. S. M (2019). A Review of Potential Beneficial Effects of Honey on Bone Health. Evid. Based Complement. Altern. Med. 2019, 8543618. doi:10.1155/2019/8543618

PubMed Abstract | CrossRef Full Text | Google Scholar

Karabagias, I. K., Nikolaou, C., and Karabagias, V. K. (2019). Volatile Fingerprints of Common and Rare Honeys Produced in Greece: In Search of PHVMs with Implementation of the Honey Code. Eur. Food Res. Technol. 245, 23–39. doi:10.1007/s00217-018-3137-x

CrossRef Full Text | Google Scholar

Karabagias, I. K., Vavoura, M. V., Nikolaou, C., Badeka, A. V., Kontakos, S., and Kontominas, M. G. (2014). Floral Authentication of Greek Unifloral Honeys Based on the Combination of Phenolic Compounds, Physicochemical Parameters and Chemometrics. Food Res. Int. 62, 753–760. doi:10.1016/j.foodres.2014.04.015

CrossRef Full Text | Google Scholar

Kashiwada, Y., Aoshima, A., Ikeshiro, Y., Chen, Y.-P., Furukawa, H., Itoigawa, M., et al. (2005). Anti-HIV Benzylisoquinoline Alkaloids and Flavonoids from the Leaves of Nelumbo nucifera, and Structure-Activity Correlations with Related Alkaloids. Bioorg. Med. Chem. 13, 443–448. doi:10.1016/j.bmc.2004.10.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Kesari, K. K., Dhasmana, A., Shandilya, S., Prabhakar, N., Shaukat, A., Dou, J., et al. (2020). Plant-Derived Natural Biomolecule Picein Attenuates Menadione Induced Oxidative Stress on Neuroblastoma Cell Mitochondria. Antioxidants (Basel) 9, 552. Epub 2020/07/08. doi:10.3390/antiox9060552

PubMed Abstract | CrossRef Full Text | Google Scholar

Kikuchi, T., Kanomi, S., Murai, Y., Kadota, S., Tsubono, K., and Ogita, Z.-I. (1986). Constituents of the Fungus Ganoderma Lucidum (Fr.) Karst. III Structures of Ganolucidic Acids A and B, New Lanostane-type Triterpenoids. Chem. Pharm. Bull. 34, 4030–4036. doi:10.1248/cpb.34.4030

CrossRef Full Text | Google Scholar

Kim, M.-J., Kwon, S.-B., Kim, M.-S., Jin, S. W., Ryu, H. W., Oh, S.-R., et al. (2016). Trifolin Induces Apoptosis via Extrinsic and Intrinsic Pathways in the NCI-H460 Human Non-small Cell Lung-Cancer Cell Line. Phytomedicine 23, 998–1004. Epub 2016/07/23. doi:10.1016/j.phymed.2016.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Kıvrak, Ş. (2016). Determination of B-Group Vitamins in Turkish Honey Using Ultra-performance Liquid Chromatography with Electrospray Ionization Coupled to Tandem Mass Spectrometry. J. Liq. Chromatogr. Relat. Technol. 39, 847–851. doi:10.1080/10826076.2016.1275001

CrossRef Full Text | Google Scholar

Koulis, G. A., Tsagkaris, A. S., Aalizadeh, R., Dasenaki, M. E., Panagopoulou, E. I., Drivelos, S., et al. (2021). Honey Phenolic Compound Profiling and Authenticity Assessment Using HRMS Targeted and Untargeted Metabolomics. Molecules 26, 2769. doi:10.3390/molecules26092769

PubMed Abstract | CrossRef Full Text | Google Scholar

Kowalski, S., Kopuncová, M., Ciesarová, Z., and Kukurová, K. (2017). Free Amino Acids Profile of Polish and Slovak Honeys Based on LC-MS/MS Method without the Prior Derivatisation. J. Food Sci. Technol. 54, 3716–3723. doi:10.1007/s13197-017-2838-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kronenberger, T., Ferreira, G. M., de Souza, A. D. F., Santos, S. d. S., Poso, A., Ribeiro, J. A., et al. (2020). Design, Synthesis and Biological Activity of Novel Substituted 3-benzoic Acid Derivatives as MtDHFR Inhibitors. Bioorg Med. Chem. 28, 115600. doi:10.1016/j.bmc.2020.115600

PubMed Abstract | CrossRef Full Text | Google Scholar

Kus, P. M., Jerkovic, I., Tuberoso, C. I., Marijanovic, Z., and Congiu, F. (2014). Cornflower (Centaurea Cyanus L.) Honey Quality Parameters: Chromatographic Fingerprints, Chemical Biomarkers, Antioxidant Capacity and Others. Food Chem. 142, 12–18. Epub 2013/09/05. doi:10.1016/j.foodchem.2013.07.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Kyriacou, M. C., Antoniou, C., Rouphael, Y., Graziani, G., and Kyratzis, A. (2021). Mapping the Primary and Secondary Metabolomes of Carob (Ceratonia Siliqua L.) Fruit and its Postharvest Antioxidant Potential at Critical Stages of Ripening. Antioxidants (Basel) 10, 57. doi:10.3390/antiox10010057

PubMed Abstract | CrossRef Full Text | Google Scholar

Lamari, F. N., Papasotiropoulos, V., Tsiris, D., Bariamis, S. E., Sotirakis, K., Pitsi, E., et al. (2018). Phytochemical and Genetic Characterization of Styles of Wild Crocus Species from the Island of Crete, Greece and Comparison to Those of Cultivated C. Sativus. Fitoterapia 130, 225–233. doi:10.1016/j.fitote.2018.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Leu, Y.-L., Hwang, T.-L., Kuo, P.-C., Liou, K.-P., Huang, B.-S., and Chen, G.-F. (2012). Constituents from Vigna Vexillata and Their Anti-inflammatory Activity. Ijms 13, 9754–9768. Epub 2012/09/06. doi:10.3390/ijms13089754

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Zhang, Y., Zhang, X., Lu, W., Liu, X., Hu, M., et al. (2020). Aucubin Exerts Anti-osteoporotic Effects by Promoting Osteoblast Differentiation. Aging 12, 2226–2245. doi:10.18632/aging.102742

PubMed Abstract | CrossRef Full Text | Google Scholar

Louveaux, J., Maurizio, A., and Vorwohl, G. (1978). Methods of Melissopalynology. Bee World, 59, 139–157. doi:10.1080/0005772x.1978.11097714

CrossRef Full Text | Google Scholar

Lusic, D., Koprivnjak, O., Curic, D., Sabatini, A. G., and Conte, L. S. (2007). Volatile Profile of Croatian Lime Tree (Tilia sp.), Fir Honeydew (Abies Alba) and Sage (Salvia Officinalis) Honey. Food Technol. Biotech. 45, 156–165.

Google Scholar

Manion-Sommerhalter, H. R., Fedoseyenko, D., Joshi, S., and Begley, T. P. (2021). Menaquinone Biosynthesis: The Mechanism of 5,8-Dihydroxy-2-Naphthoate Synthase (MqnD). Biochemistry 60, 1947–1951. doi:10.1021/acs.biochem.1c00257

PubMed Abstract | CrossRef Full Text | Google Scholar

Margaoan, R., Topal, E., Balkanska, R., Yucel, B., Oravecz, T., Cornea-Cipcigan, M., et al. (2021). Monofloral Honeys as a Potential Source of Natural Antioxidants, Minerals and Medicine. Antioxidants-Basel 10, 1023. doi:10.3390/antiox10071023

PubMed Abstract | CrossRef Full Text | Google Scholar

Mattes, B. R., Clausen, T. P., and Reichardt, P. B. (1987). Volatile Constituents of Balsam Poplar: The Phenol Glycoside Connection. Phytochemistry 26, 1361–1366. doi:10.1016/s0031-9422(00)81813-0

CrossRef Full Text | Google Scholar

Mezrag, A., Malafronte, N., Bouheroum, M., Travaglino, C., Russo, D., Milella, L., et al. (2017). Phytochemical and Antioxidant Activity Studies on Ononis Angustissima L. Aerial Parts: Isolation of Two New Flavonoids. Nat. Prod. Res. 31, 507–514. doi:10.1080/14786419.2016.1195381

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaud, K. M., Irwin, R. E., Barber, N. A., and Adler, L. S. (2019). Preinfection Effects of Nectar Secondary Compounds on a Bumble Bee Gut Pathogen. Environ. Entomol. 48, 685–690. doi:10.1093/ee/nvz018

PubMed Abstract | CrossRef Full Text | Google Scholar

Miehe-Steier, A., Roscher, C., Reichelt, M., Gershenzon, J., and Unsicker, S. B. (2015). Light and Nutrient Dependent Responses in Secondary Metabolites of Plantago Lanceolata Offspring Are Due to Phenotypic Plasticity in Experimental Grasslands. Plos One 10, e0136073. doi:10.1371/journal.pone.0136073

PubMed Abstract | CrossRef Full Text | Google Scholar

Mock, A., Warta, R., Dettling, S., Brors, B., Jäger, D., and Herold-Mende, C. (2018). MetaboDiff: An R Package for Differential Metabolomic Analysis. Bioinformatics 34, 3417–3418. doi:10.1093/bioinformatics/bty344

PubMed Abstract | CrossRef Full Text | Google Scholar

Monje, M.-C., Privat, C., Gastine, V., and Nepveu, F. (2002). Determination of Ethylphenol Compounds in Wine by Headspace Solid-phase Microextraction in Conjunction with Gas Chromatography and Flame Ionization Detection. Anal. Chim. Acta 458, 111–117. doi:10.1016/s0003-2670(01)01528-8

CrossRef Full Text | Google Scholar

Montoro, P., D'Urso, G., Kowalczyk, A., and Tuberoso, C. I. G. (2021). LC-ESI/LTQ-Orbitrap-MS Based Metabolomics in Evaluation of Bitter Taste of Arbutus Unedo Honey. Molecules 26, 2765. doi:10.3390/molecules26092765

PubMed Abstract | CrossRef Full Text | Google Scholar

Morgan, A. M. A., Lee, H. W., Lee, S.-H., Lim, C.-H., Jang, H.-D., and Kim, Y. H. (2014). Anti-osteoporotic and Antioxidant Activities of Chemical Constituents of the Aerial Parts of Ducrosia Ismaelis. Bioorg. Med. Chem. Lett. 24, 3434–3439. doi:10.1016/j.bmcl.2014.05.077

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, H. T. L., Panyoyai, N., Kasapis, S., Pang, E., and Mantri, N. (2019). Honey and its Role in Relieving Multiple Facets of Atherosclerosis. Nutrients 11, 167. doi:10.3390/nu11010167

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, V. B., Wang, S.-L., Nguyen, A. D., Vo, T. P. K., Zhang, L.-J., Nguyen, Q. V., et al. (2018). Isolation and Identification of Novel α-amylase Inhibitors from Euonymus Laxiflorus Champ. Res. Chem. Intermed. 44, 1411–1424. doi:10.1007/s11164-017-3175-1

CrossRef Full Text | Google Scholar

Nikaido, T., Ohmoto, T., Noguchi, H., Kinoshita, T., Saitoh, H., and Sankawa, U. (1981). Inhibitors of Cyclic AMP Phosphodiesterase in Medicinal Plants. Planta Med. 43, 18–23. doi:10.1055/s-2007-971466

PubMed Abstract | CrossRef Full Text | Google Scholar

Oddo, L. P., Piro, R., Bruneau, É., Guyot-Declerck, C., Ivanov, T., Piskulová, J., et al. (2004). Main European Unifloral Honeys: Descriptive Sheets. Apidologie 35, S38–S81. doi:10.1051/apido:2004049

CrossRef Full Text | Google Scholar

Olajide, O. A., and Sarker, S. D. (2020). “Anti-inflammatory Natural Products,” in Annual Reports in Medicinal Chemistry. Editors SD Sarker, and L Nahar (Cambridge, Massachusetts: Academic Press), 55, 153–177. doi:10.1016/bs.armc.2020.02.002

CrossRef Full Text | Google Scholar

Paiva, L. B. d., Goldbeck, R., Santos, W. D. d., and Squina, F. M. (2013). Ferulic Acid and Derivatives: Molecules with Potential Application in the Pharmaceutical Field. Braz. J. Pharm. Sci. 49, 395–411. doi:10.1590/s1984-82502013000300002

CrossRef Full Text | Google Scholar

Peana, A. T., Bassareo, V., and Acquas, E. (2019). Not Just from Ethanol. Tetrahydroisoquinolinic (TIQ) Derivatives: From Neurotoxicity to Neuroprotection. Neurotox. Res. 36, 653–668. doi:10.1007/s12640-019-00051-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pirini, A., Conte, L. S., Francioso, O., and Lercker, G. (1992). Capillary Gas Chromatographic Determination of Free Amino Acids in Honey as a Means of Discrimination between Different Botanical Sources. J. High. Resol. Chromatogr. 15, 165–170. doi:10.1002/jhrc.1240150306

CrossRef Full Text | Google Scholar

Potacnjak, I., Marinic, J., Baticic, L., Simic, L., Broznic, D., and Domitrovic, R. (2020). Aucubin Administered by Either Oral or Parenteral Route Protects against Cisplatin-Induced Acute Kidney Injury in Mice. Food Chem. Toxicol. 142, 111472. doi:10.1016/j.fct.2020.111472

PubMed Abstract | CrossRef Full Text | Google Scholar

Potocnjak, I., Marinic, J., Baticic, L., Simic, L., Broznic, D., and Domitrovic, R. (2021). Aucubin Administered by Either Oral or Parenteral Route Protects against Cisplatin-Induced Acute Kidney Injury in Mice. Febs Open Bio 11, 487.

Google Scholar

Pretorius, C. J., Tugizimana, F., Steenkamp, P. A., Piater, L. A., and Dubery, I. A. (2021). Metabolomics for Biomarker Discovery: Key Signatory Metabolic Profiles for the Identification and Discrimination of Oat Cultivars. Metabolites 11, 165. doi:10.3390/metabo11030165

PubMed Abstract | CrossRef Full Text | Google Scholar

PubChem (n.d.). National Center for Biotechnology Information (2016) PubChem Database. Available at: https://pubchem.ncbi.nlm.nih.gov/search/search.cgi (Accessed Jan 14, 2016).

Google Scholar

Randjelović, P., Veljković, S., Stojiljković, N., Sokolović, D., Ilić, I., Laketić, D., et al. (2015). The Beneficial Biological Properties of Salicylic Acid. Acta Fac. medicae Naissensis 32, 259–265. doi:10.1515/afmnai-2015-0026

CrossRef Full Text | Google Scholar

Raviv, Z., Cohen, S., and Reischer-Pelech, D. (2013). The Anti-cancer Activities of Jasmonates. Cancer Chemother. Pharmacol. 71, 275–285. doi:10.1007/s00280-012-2039-z

PubMed Abstract | CrossRef Full Text | Google Scholar

REG-625/2017 (2017). Regulation (EU) 625/2017, on Official Controls and Other Official Activities Performed to Ensure the Application of Food and Feed Law, Rules on Animal Health and Welfare, Plant Health and Plant Protection Products.

Google Scholar

Rodríguez-Flores, M. S., Falcão, S. I., Escuredo, O., Seijo, M. C., and Vilas-Boas, M. (2021). Description of the Volatile Fraction of Erica Honey from the Northwest of the Iberian Peninsula. Food Chem. 336, 127758. doi:10.1016/j.foodchem.2020.127758

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero-Pérez, A. I., Ibern-Gómez, M., Lamuela-Raventós, R. M., and de La Torre-Boronat, M. C. (1999). Piceid, the Major Resveratrol Derivative in Grape Juices. J. Agric. Food Chem. 47, 1533–1536. doi:10.1021/jf981024g

PubMed Abstract | CrossRef Full Text | Google Scholar

Saleem, H., Htar, T. T., Naidu, R., Ahmad, I., Zengin, G., Ahmad, M., et al. (2019). Investigations into the Therapeutic Effects of Aerial and Stem Parts of Buxus Papillosa C.K. Schneid.: In Vitro Chemical, Biological and Toxicological Perspectives. J. Pharm. Biomed. Analysis 166, 128–138. doi:10.1016/j.jpba.2019.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Salonen, A., Virjamo, V., Tammela, P., Fauch, L., and Julkunen-Tiitto, R. (2017). Screening Bioactivity and Bioactive Constituents of Nordic Unifloral Honeys. Food Chem. 237, 214–224. doi:10.1016/j.foodchem.2017.05.085

PubMed Abstract | CrossRef Full Text | Google Scholar

SANTE/11956/2016 (2018). Technical Guidelines for Determining the Magnitude of Pesticide Residues in Honey and Setting Maximum Residue Levels in Honey. rev 9

Google Scholar

Scotti, M. T., Fernandes, M. B., Ferreira, M. J. P., and Emerenciano, V. P. (2007). Quantitative Structure-Activity Relationship of Sesquiterpene Lactones with Cytotoxic Activity. Bioorg. Med. Chem. 15, 2927–2934. doi:10.1016/j.bmc.2007.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sergiel, I., Pohl, P., and Biesaga, M. (2014). Characterisation of Honeys According to Their Content of Phenolic Compounds Using High Performance Liquid Chromatography/tandem Mass Spectrometry. Food Chem. 145, 404–408. doi:10.1016/j.foodchem.2013.08.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Sertic, M., Crkvencic, M., Mornar, A., Pilepic, K. H., Nigovic, B., and Males, Z. (2015). Analysis of Aucubin and Catalpol Content in Different Plant Parts of Four Globularia Species. J. Appl. Bot. Food Qual. 88, 209–214. doi:10.5073/JABFQ.2015.088.030

CrossRef Full Text | Google Scholar

Sesta, G., Piana, M. L., Persano Oddo, L., Lusco, L., and Belligoli, P. (2008). Methyl Anthranilate inCitrushoney. Analytical Method and Suitability as a Chemical Marker. Apidologie 39, 334–342. doi:10.1051/apido:2008009

CrossRef Full Text | Google Scholar

Shen, C.-Y., Jiang, J.-G., Yang, L., Wang, D.-W., and Zhu, W. (2017). Anti-ageing Active Ingredients from Herbs and Nutraceuticals Used in Traditional Chinese Medicine: Pharmacological Mechanisms and Implications for Drug Discovery. Br. J. Pharmacol. 174, 1395–1425. doi:10.1111/bph.13631

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva-Beltran, N. P., Umsza-Guez, M. A., Rodrigues, D. M. R., Galvez-Ruiz, J. C., Castro, T. L. D., and Balderrama-Carmona, A. P. (2021). Comparison of the Biological Potential and Chemical Composition of Brazilian and Mexican Propolis. Appl. Sci-Basel. 11, 11417. doi:10.3390/app112311417

CrossRef Full Text | Google Scholar

Stavropoulou, M. I., Termentzi, A., Kasiotis, K. M., Cheilari, A., Stathopoulou, K., Machera, K., et al. (2021). Untargeted Ultrahigh-Performance Liquid Chromatography-Hybrid Quadrupole-Orbitrap Mass Spectrometry (UHPLC-HRMS) Metabolomics Reveals Propolis Markers of Greek and Chinese Origin. Molecules 26, 456. doi:10.3390/molecules26020456

PubMed Abstract | CrossRef Full Text | Google Scholar

Štefániková, J., Martišová, P., Šnirc, M., Šedík, P., and Vietoris, V. (2021). Screening of the Honey Aroma as a Potential Essence for the Aromachology. Appl. Sci. 11, 8177. doi:10.3390/app11178177

CrossRef Full Text | Google Scholar

Stefanikova, J., Martisova, P., Snirc, M., Sedik, P., and Vietoris, V. (2021). Screening of the Honey Aroma as a Potential Essence for the Aromachology. Appl. Sci-Basel. 11, 8177. doi:10.3390/app11178177

CrossRef Full Text | Google Scholar

Strid, A. (2016). Atlas of the Aegean Flora-Part 2: Maps: Botanic Garden and Botanical Museum Berlin.

Google Scholar

Sun, K., Su, C., Li, W., Gong, Z., Sha, C., and Liu, R. (2021). Quality Markers Based on Phytochemical Analysis and Anti-inflammatory Screening: An Integrated Strategy for the Quality Control of Dalbergia Odorifera by UHPLC-Q-Orbitrap HRMS. Phytomedicine 84, 153511. doi:10.1016/j.phymed.2021.153511

PubMed Abstract | CrossRef Full Text | Google Scholar

Suwannasom, N., Kao, I., Pruss, A., Georgieva, R., and Bäumler, H. (2020). Riboflavin: The Health Benefits of a Forgotten Natural Vitamin. Int. J. Mol. Sci. 21, 950. doi:10.3390/ijms21030950

PubMed Abstract | CrossRef Full Text | Google Scholar

Sytar, O., Hemmerich, I., Zivcak, M., Rauh, C., and Brestic, M. (2018). Comparative Analysis of Bioactive Phenolic Compounds Composition from 26 Medicinal Plants. Saudi J. Biol. Sci. 25, 631–641. doi:10.1016/j.sjbs.2016.01.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, Q.-W., Ouyang, M.-A., and Wu, Z.-J. (2012). A New Seco-Neolignan Glycoside from the Root Bark of Ailanthus Altissima. Nat. Prod. Res. 26, 1375–1380. doi:10.1080/14786419.2011.587187

PubMed Abstract | CrossRef Full Text | Google Scholar

Thai, T. H., Hai, N. T., Hien, N. T., Ha, C. T. T., Cuong, N. T., Binh, P. T., et al. (2017). Cytotoxic Constituents of Mallotus Microcarpus. Nat. Product. Commun. 12, 1934578X1701200325. doi:10.1177/1934578x1701200325

CrossRef Full Text | Google Scholar

Thévenot, E. A., Roux, A., Xu, Y., Ezan, E., and Junot, C. (2015). Analysis of the Human Adult Urinary Metabolome Variations with Age, Body Mass Index, and Gender by Implementing a Comprehensive Workflow for Univariate and OPLS Statistical Analyses. J. Proteome Res. 14, 3322–3335. doi:10.1021/acs.jproteome.5b00354

PubMed Abstract | CrossRef Full Text | Google Scholar

Thrasyvoulou, A., and Manikis, J. (1995). Some Physicochemical and Microscopic Characteristics of Greek Unifloral Honeys. Apidologie 26, 441–452. doi:10.1051/apido:19950601

PubMed Abstract | CrossRef Full Text | Google Scholar

Thrasyvoulou, A., Tananaki, C., Goras, G., Karazafiris, E., Dimou, M., Liolios, V., et al. (2018). Legislation of Honey Criteria and Standards. J. Apic. Res. 57, 88–96. doi:10.1080/00218839.2017.1411181

CrossRef Full Text | Google Scholar

Tsagkarakou, A. S., Chasapi, S. A., Koulas, S. M., Tsialtas, I., Kyriakis, E., Drakou, C. E., et al. (2021). Structure Activity Relationship of the Binding of P-Coumaroyl Glucose to Glycogen Phosphorylase and its Effect on Hepatic Cell Metabolic Pathways. Eur. J. Med. Chem. Rep. 3, 100011. doi:10.1016/j.ejmcr.2021.100011

CrossRef Full Text | Google Scholar

Tuberoso, C. I. G., Bifulco, E., Caboni, P., Sarais, G., Cottiglia, F., and Floris, I. (2011). Lumichrome and Phenyllactic Acid as Chemical Markers of Thistle (Galactites Tomentosa Moench) Honey. J. Agric. Food Chem. 59, 364–369. doi:10.1021/jf1039074

PubMed Abstract | CrossRef Full Text | Google Scholar

Tuberoso, C. I. G., Jerković, I., Bifulco, E., Marijanovic, Z., Congiu, F., and Bubalo, D. (2012). Riboflavin and Lumichrome in Dalmatian Sage Honey and Other Unifloral Honeys Determined by LC-DAD Technique. Food Chem. 135, 1985–1990. doi:10.1016/j.foodchem.2012.06.096

PubMed Abstract | CrossRef Full Text | Google Scholar

Tundis, R., Loizzo, M., Menichini, F., Statti, G., and Menichini, F. (2008). Biological and Pharmacological Activities of Iridoids: Recent Developments. Mrmc 8, 399–420. doi:10.2174/138955708783955926

CrossRef Full Text | Google Scholar

Turski, M. P., Chwil, S., Turska, M., Chwil, M., Kocki, T., Rajtar, G., et al. (2016). An Exceptionally High Content of Kynurenic Acid in Chestnut Honey and Flowers of Chestnut Tree. J. Food Compos. Analysis 48, 67–72. doi:10.1016/j.jfca.2016.02.003

CrossRef Full Text | Google Scholar

von der Ohe, W., Oddo, L. P., Piana, M. L., Morlot, M., and Martin, P. (2004). Harmonized Methods of Melissopalynology. Apidologie 35, S18–S25. doi:10.1051/apido:2004050

CrossRef Full Text | Google Scholar

Wang, M., Li, J., Rangarajan, M., Shao, Y., LaVoie, E. J., Huang, T.-C., et al. (1998). Antioxidative Phenolic Compounds from Sage (Salvia Officinalis). J. Agric. Food Chem. 46, 4869–4873. doi:10.1021/jf980614b

CrossRef Full Text | Google Scholar

Wolf, S., Schmidt, S., Müller-Hannemann, M., and Neumann, S. (2010). In Silico fragmentation for Computer Assisted Identification of Metabolite Mass Spectra. BMC Bioinforma. 11, 148. doi:10.1186/1471-2105-11-148

CrossRef Full Text | Google Scholar

Xagoraris, M., Chrysoulaki, F., Revelou, P. K., Alissandrakis, E., Tarantilis, P. A., and Pappas, C. S. (2021). Unifloral Autumn Heather Honey from Indigenous Greek Erica Manipuliflora Salisb.: SPME/GC-MS Characterization of the Volatile Fraction and Optimization of the Isolation Parameters. Foods 10, 2487. doi:10.3390/foods10102487

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, L.-H., Li, H.-B., Huang, Y.-X., Qin, D.-P., Zhang, C.-F., Wang, Z.-Z., et al. (2021). Research on Chemical Constituents from Artemisia Annua I. Zhongguo Zhongyao Zazhi 46, 1160–1167. doi:10.19540/j.cnki.cjcmm.20210107.601

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, S., Sun, F., Ruan, J., Yan, J., Huang, P., Wang, J., et al. (2019). Anti-inflammatory Constituents from Cortex Dictamni. Fitoterapia 134, 465–473. doi:10.1016/j.fitote.2019.03.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Battesti, M.-J., Paolini, J., Muselli, A., Tomi, P., and Costa, J. (2012). Melissopalynological Origin Determination and Volatile Composition Analysis of Corsican "Erica Arborea Spring Maquis" Honeys. Food Chem. 134, 37–47. doi:10.1016/j.foodchem.2012.02.026

CrossRef Full Text | Google Scholar

Yang, Y., Zhang, H., Zuo, J., Gong, X., Yi, F., Zhu, W., et al. (2019). Advances in Research on the Active Constituents and Physiological Effects of Ganoderma Lucidum. Biomed. Dermatol 3, 6. doi:10.1186/s41702-019-0044-0

CrossRef Full Text | Google Scholar

Yao, J., Li, Y., Jin, Y., Chen, Y., Tian, L., and He, W. (2021). Synergistic Cardioptotection by Tilianin and Syringin in Diabetic Cardiomyopathy Involves Interaction of TLR4/NF-Κb/nlrp3 and PGC1a/SIRT3 Pathways. Int. Immunopharmacol. 96, 107728. doi:10.1016/j.intimp.2021.107728

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Liu, K., Shi, M., Xie, L., Deng, M., Chen, H., et al. (2021). Therapeutic Potential of Catalpol and Geniposide in Alzheimer's and Parkinson's Diseases: A Snapshot of Their Underlying Mechanisms. Brain Res. Bull. 174, 281–295.doi:10.1016/j.brainresbull.2021.06.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: mediterranean honey, pollen, ikaria honey, metabolomics, UHPLC-HRMS, orbitrap, melissopalynology, chemometrics

Citation: Kasiotis KM, Baira E, Iosifidou S, Bergele K, Manea-Karga E, Theologidis I, Barmpouni T, Tsipi D and Machera K (2022) Characterization of Ikaria Heather Honey by Untargeted Ultrahigh-Performance Liquid Chromatography-High Resolution Mass Spectrometry Metabolomics and Melissopalynological Analysis. Front. Chem. 10:924881. doi: 10.3389/fchem.2022.924881

Received: 20 April 2022; Accepted: 18 May 2022;
Published: 22 July 2022.

Edited by:

Paola Donato, Università degli Studi di Messina, Italy

Reviewed by:

Siew Hua Gan, Monash University Malaysia, Malaysia
Francesca Aiello, University of Calabria, Italy

Copyright © 2022 Kasiotis, Baira, Iosifidou, Bergele, Manea-Karga, Theologidis, Barmpouni, Tsipi and Machera. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Konstantinos M. Kasiotis, Sy5LYXNpb3Rpc0BicGkuZ3I=; Despina Tsipi, ZC50c2lwaUBhYWRlLmdy

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.