Next Article in Journal
Carpachromene Ameliorates Insulin Resistance in HepG2 Cells via Modulating IR/IRS1/PI3k/Akt/GSK3/FoxO1 Pathway
Next Article in Special Issue
Effect of Ultrasound-Assisted Extraction of Carotenoids from Papaya (Carica papaya L. cv. Sweet Mary) Using Vegetable Oils
Previous Article in Journal
Metabolite Fingerprinting Using 1H-NMR Spectroscopy and Chemometrics for Classification of Three Curcuma Species from Different Origins
Previous Article in Special Issue
Schisandrin B Attenuates Hepatic Stellate Cell Activation and Promotes Apoptosis to Protect against Liver Fibrosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Botanical Origin Differentiation of Malaysian Stingless Bee Honey Produced by Heterotrigona itama and Geniotrigona thoracica Using Chemometrics

Faculty of Science, Universiti Tunku Abdul Rahman (UTAR), Kampar 31900, Malaysia
*
Author to whom correspondence should be addressed.
Submission received: 19 November 2021 / Revised: 8 December 2021 / Accepted: 12 December 2021 / Published: 16 December 2021
(This article belongs to the Special Issue Plant Foods Ingredients as Functional Foods and Nutraceuticals)

Abstract

:
Stingless bee honey, specifically honeydew honey, is generally valued for its better health benefits than those of most blossom types. However, scientific studies about the differentiation of stingless bee honey based on honeydew and blossom origins are very limited. In this study, 13C NMR spectroscopy was employed to quantify the seven major sugar tautomers in stingless bee honey samples, and the major sugar compositions of both honeydew and blossom types were found not significantly different. However, several physicochemical properties of honeydew honey including moisture content, free acidity, electrical conductivity, ash content, acetic acid, diastase, hydrogen peroxide, and mineral elements levels were significantly higher; while total soluble solid, proline, and hydroxymethylfurfural were significantly lower than blossom honey. Greater antioxidant capacity in honeydew honey was proven with higher total phenolic compounds, ABTS, DPPH, superoxide radical scavenging activities, peroxyl radical inhibition, iron chelation, and ferric reducing power. Using principal component analysis (PCA), two clusters of stingless bee honey from the honeydew and blossom origin were observed. PCA also revealed that the differentiation between honeydew and blossom origin of stingless bee honey is possible with certain physicochemical and antioxidant parameters. The combination of NMR spectroscopy and chemometrics are suggested to be useful to determine the authenticity and botanical origin of stingless bee honey.

1. Introduction

Stingless bee is the largest group of eusocial bees on Earth with more than 500 described species. It is distributed in tropical and subtropical regions including tropical America, Africa, Australia, and Southeast Asia, where it is 50 times greater in number than honey bees [1]. Less attention and study have been focused on these bees and related products presumably because they are not found in North America and European regions [2].
The production of honey begins with nectar collection by bees from plants, whereby the bees transform and process the nectar with specific substances of their own. Instead of hexagonal-shaped combs, stingless bees store their nectar in small resin pots. To date, 32 stingless bee species have been documented in Malaysia [3], and Heterotrigona itama and Geniotrigona thoracica are the most common domesticated stingless bees for meliponiculture purpose. Recently, stingless bee honey, which is locally known as kelulut honey, is receiving more attention and demand from the public due to its distinctive sensory characteristics and health-promoting properties. Other than floral nectar, bees also harvest honeydew honey from plant secretions. Honeydew honey is highly valued because it is considered to possess more health benefits than floral or blossom honey [4]. Several studies also demonstrated that the antibacterial and antioxidant activities of honeydew honey are superior to those of blossom honey [5].
Although the quality parameters for honey have been defined by Codex Alimentarius Commission [6] and Council Directive of the European Union [7], they are not applicable for honey originating from tropical and subtropical regions, especially stingless bee honey [8]. Despite a specification standard for stingless bee honey formulated by the Department of Standards Malaysia [9] to regulate the quality of locally produced stingless bee honey, the quality requirements do not include differentiation of honeydew origin from blossom origin honeys. Furthermore, scientific studies about the differentiation of stingless bee honey based on honeydew and blossom origins are very limited. Hence, it is crucial to differentiate stingless bee honey to ensure authenticity and avoid adulteration. A systematic and comprehensive research that can profile the physicochemical properties and bioactivities of stingless bee honey is needed.
Thus, the main interest of this study was to differentiate the botanical origins of stingless bee honey using chemometric analysis. Such analysis is not only crucial for the characterization of stingless bee honey using physicochemical characteristics and antioxidant capacity but also important to ensure the quality of local stingless bee honey products for both domestic consumption and exports.

2. Results

2.1. Physicochemical Properties

2.1.1. Major Sugar Composition

13C NMR chemical shifts of the seven tautomers relative to 1,4-dioxane are reported in Supplementary Material Table S1. Tables S2 and S3 show the 13C NMR chemical shifts of the isoglucose and artificial honey, respectively. The comparison of actual concentration and 13C NMR determined concentration of sugars and tautomers is also shown in the same tables. The 13C NMR spectra of model compounds (Figure S1a–c), isoglucose (Figure S2), and artificial honey (Figure S3) are also provided as Supplementary Material. No significant changes in chemical shifts for each tautomer in isoglucose and artificial honey relative to the individual model compounds, indicating that the 13C chemical shifts of a sugar tautomer are not influenced by the presence of other sugars in a solution including honey [10,11]. The assignment of all carbon signals for honeydew and blossom honey (Figure S4) is based on the assigned 13C NMR spectra of the model compounds, isoglucose, and artificial honey. Table 1 summarizes the average concentration of each tautomer in honeydew and blossom honey. The content of glucose and fructose was not significantly different between honeydew honey and blossom honey. The content of sucrose in honeydew honey was too low to be integrated by 13C NMR spectroscopy.

2.1.2. Maturity of Honey

As shown in Table 2, although the color intensity of honeydew honey with an average value of 173.75 mAU was darker than blossom honey (146.67 mAU), the difference was not significant. The moisture content in stingless bee honey samples ranged from 22.50% to 28.00%, with significant higher moisture content in honeydew honey (27.03%) than blossom honey (25.67%). The water activity in both honeydew honey (0.63) and blossom honey (0.62) were rather similar and slightly higher than 0.6. Contrary to moisture content, the total soluble solid in honeydew honey (72.97 °Brix) was detected to be significantly lower than blossom honey (74.32 °Brix). Proline content in honeydew honey (522.38 mg/kg) was found to be significantly lower than blossom honey (563.11 mg/kg).

2.1.3. Purity of Honey

In Table 3, the values of both electrical conductivity and ash content showed great variability among botanical origins whereby values found in honeydew honey (0.45 mS/cm; 0.14%, g/100 g) were significantly higher than blossom honey (0.36 mS/cm; 0.04%, g/100 g). The diastase level in stingless bee honey samples reported in this study ranged from 1.83 to 3.04 Schade units. The hydrogen peroxide level of stingless bee honey samples in this study ranged from 111.07 to 192.82 µmol/L with a significantly higher value detected in honeydew honey (183.51 µmol/L) than blossom honey (137.61 µmol/L).

2.1.4. Deterioration State of Honey

As shown in Table 4, although the pH values of honeydew honey and blossom honey were not to be significantly different, the free acidity value of honeydew honey (90.96 meq/kg) was found to the significantly higher than blossom honey (76.64 meq/kg). Although higher levels of organic acids were determined in honeydew honey, only the acetic acid level of honeydew honey (0.10 g/kg) was significantly higher than blossom honey (0.07 g/kg). The HMF content of both honeydew and blossom honeys were below 20 mg/kg.

2.1.5. Mineral Profile

As displayed in Table 5, significantly higher levels of sodium (Na), potassium (K), magnesium (Mg), calcium (Ca), iron (Fe), and zinc (Zn) were measured in honeydew honey. However, the quantity of copper (Cu) and aluminum (Al) was too little to be quantified in both honey types. A significantly higher total content of mineral elements was determined in honeydew honey (732.91 mg/kg) than blossom honey (573.86 mg/kg).

2.2. Antioxidant Properties

Table 6 shows honeydew honey possessed greater antioxidant capacities than blossom honey. Other than significantly higher total phenolic content (104.09 mg GAE/kg), scavenging activities against ABTS, DPPH and superoxide radicals (63.15%, 35.65%, and 76.10%), peroxyl radical inhibition (5.71 µmol TE/g), iron chelation (18.94%), and ferric reducing power (3.18 mmol Fe(II)/kg) were all significantly higher than blossom honey.

2.3. Chemometric Analysis

As shown in Table 7, both first principal component (PC1) and second principal component (PC2) represented 56.40% of the variance. According to correlation coefficient, the parameters that most associated with PC1 were ash content (−0.822), hydrogen peroxide (−0.886), free acidity (−0.842), total mineral elements (−0.916), K (−0.817), Mg (−0.876), Ca (−0.863), total phenolic compounds (−0.817), and ferric reducing power (−0.907). Considering only these parameters, another PCA was generated with 84.60% of the total data variance. The correlation coefficients of these parameters were ash content (−0.891), hydrogen peroxide (−0.878), free acidity (−0.811), total mineral elements (−0.945), K (−0.892), Mg (−0.876), Ca (−0.810), total phenolic compounds (−0.851), and ferric reducing power (−0.945). This analysis was able to differentiate all the analyzed honey samples into two clusters based on honeydew and blossom origins (Figure 1).

3. Discussion

1H NMR spectroscopy is commonly used in food analysis; however, highly complex 1H NMR spectra can complicate the interpretation and profiling of food metabolites including honey [10]. Such an issue can be overcome by 13C NMR spectroscopy which has been largely employed for structural studies. The 13C NMR spectra obtained under proton decoupling spreads over a larger chemical shift range, and this makes it easier to identify the sugar resonances in honey [10]. To our best knowledge, the relevant scientific reports about using 13C NMR for the analysis of sugars in honey are still limited. Hence, the major sugar molecules of stingless bee honey, including mono- and di-saccharides, and their respective tautomers were identified and quantified in this study using 13C NMR spectroscopy.
Reducing sugars undergo mutarotation and form two or more species known as tautomers. Mutarotation happens when the hemiketal ring opens and closes and with α or β configuration [12,13,14]. For instance, d-glucose in aqueous solution appears in two forms, which are α-d-glucopyranose and β-d-glucopyranose; fructose appears in four forms in aqueous solution, namely α-d-fructopyranose, β-d-fructopyranose, α-d-fructofuranose, and β-d-fructofuranose. However, the composition of tautomers varies for each respective reducing sugar in different solvents [12,13,14]. Hence, accurate determination of the composition of tautomers of raw honey in aqueous solution is crucial. Table 1 shows the composition of the tautomers in stingless bee honey, determined by 13C NMR spectroscopy, which is in agreement with the literature [10,11,12,13]. Both glucose and fructose are the major monosaccharides found in honey formed by the hydrolysis of sucrose by enzyme invertase [15]. The reducing sugar content in stingless bee honey samples ranging from 73% to 73.61% also met the requirement set by the Department of Standards Malaysia, which is not more than 85.0% [9]. Honeydew honey was found to possess very minimal sucrose to be integrated by 13C NMR spectroscopy and sucrose content obtained from blossom honey in this study (0.40%) met the established requirement, which is not more than 9.5%. The values obtained by other authors in stingless bee honeys also met the requirements [16,17]. However, the presence of trehalulose in tested stingless bee honey samples was not detected by 13C NMR spectroscopy. The chemical shift for C1 of α-d-glucopyranose in trehalulose (99.3–101.3 ppm) [18,19,20] was not observed in stingless bee honey samples in this study. Moreover, the obtained chemical shift for C1 of α-d-glucopyranose tautomer was consistent with the value from a standard glucose compound as reported in another study (92.7 ppm) [10]. Hence, contrary to an earlier study [18], such disaccharide could be absent or present only in very low quantity in the tested stingless bee honey samples produced by H. itama and G. thoracica.
It has been reported that honeydew honey presented a lower value of the sum of glucose plus fructose than blossom honey because blossom honey is usually tested with higher glucose content [21,22]. However, this study and other authors showed different results in which there were no significant differences for both types of honey [23]. Glucose tends to precipitate and crystalize because it is less water soluble than fructose. The tendency of crystallization in honey can be evaluated based on the fructose-to-glucose ratio (F/G) and the glucose-to-moisture ratio (G/M). According to a previous study, F/G not more than 1.14 is associated with honey crystallization, whereas values 1.33 and above indicate a slower crystallization [24]. Other authors suggested that the G/M ratio could be a better indicator for honey crystallization with values not more than 1.7 indicate slow or no crystallizations, while a higher value up to 2 suggests rapid crystallization [21,22,25]. In this study, the F/G values of honeydew and blossom honeys were 1.06 and 1.08, respectively, whereas the G/M values of honeydew and blossom honeys were 1.31 and 1.38, respectively. Hence, these values suggest a slower crystallization of both honeydew and blossom honeys. Although the average value of total sugars in honeydew honey (73%) was lower than blossom honey (74.01%), there was no significant different in the total content of glucose and fructose between honeydew honey and blossom honey. Hence, the major sugar composition is suggested to be consistent among stingless bee honey.
A previous study also stated that honeydew honey is generally darker than the blossom origin honey [22]. The color intensity of raw honey is mainly related to the botanical origin and composition including mineral, pollen, and phenolic compounds [26]. In this study, honeydew honey with darker color (Table 2) was tested with higher mineral content (Table 5) and phenolic compounds (Table 6). The moisture content in stingless bee honey samples ranged from 22.50% to 28.00%, clearly exceeding the limits established by Codex Alimentarius Commission and Council Directive of the European Union, which stipulated the moisture content must be below 20% [6,7]. However, the values still comply with the Malaysian Standard for stingless bee honey established by the Department of Standards Malaysia that the moisture content of stingless bee honey should not exceed 35.0% [9]. The moisture content in stingless bee honey is one of the unique parameters demonstrating the importance of an international legislation geared toward stingless bee honeys. The significant difference in the moisture contents of honey bee honey and stingless bee honey is mainly due to the different strategies used to produce honey. Honey bees remove moisture by using their wings and store honey in sealed honeycombs, while stingless bees do not vaporize honey using their wings, and the honey has contact with the air until the honey pots are fully enclosed [17,27]. Higher moisture content of stingless bee honey than honey bee honey was also shown in other studies [16,28,29]. In addition to that, the difference in moisture content between honeydew honey and blossom honey could be due to a different nectar source [30,31].
Despite several articles reporting no significant difference for moisture content and water activity between honeydew and blossom honeys [21,22], stingless bee honeydew honey in this study was tested to have significantly higher moisture content. Moisture content is highly associated with water activity; water activity is a better indicator of microbial growth because it elucidates the amount of free water available for microorganisms [17]. The values of water activity obtained in this study were lower than previously reported values of 0.76–0.87 from stingless bee honey [17]. Water activity lower than 0.60 is able to prevent microbial growth including osmophilic yeasts; hence, stingless bee honey is generally more susceptible to microbial fermentation [17]. This parameter can be considered to be included in legislation to ensure the stability of stingless bee honey toward microbial fermentation. However, there was no significant difference in water activity between honeydew honey and blossom honey [21]. The total soluble solid in stingless bee honey (72.00–77.50 °Brix) was found to be in agreement with the values reported by previous studies from 60.85 to 72.25 °Brix [17] and from 55.20 to 76.10 °Brix [32]. Although different °Brix values in honey samples could be due to different amounts of sugars, the °Brix value is directly related to the amount of sugar in the honey [24,25]. The content of glucose and fructose is not significantly different between honeydew honey and blossom honey in this study as described above. Hence, it was suggested that lower total soluble solid could be linked to a significantly lower amount of proline and hydroxymethylfurfural found in honeydew honey.
Proline is the major component with 50–85% of the total amino acids in honey [33]. Proline content in honey has to be at least 180 mg/kg to be accepted internationally [16]. All stingless bee honey samples in this study met this requirement with the range of 454.80–670.27 mg/kg. The proline content of stingless bee honey in this study was found to be higher than the reported values by another study, ranging from 276.11 to 498.52 mg/kg [16]. Since the presence of proline in honey is related to its botanical source, it can be suggested that such a variation is due to the botanical preference of bees during nectar collection [34]. However, the high variability of the proline content in honey is unable to differentiate between the honeydew and blossom origins [22,32].
Both electrical conductivity and ash content are reliable parameters to differentiate between honeydew and blossom honeys [22]. Although there is no reference range recommended by Codex Alimentarius Commission and Council Directive of the European Union for the electrical conductivity of stingless bee honey [6,7], the electrical conductivity of stingless bee honey as shown in Table 3 (0.28–0.52 mS/cm) was found to be similar to the previously reported values, ranging from 0.32 to 1.10 mS/cm [28], ranging from 0.35 to 0.76 mS/cm [35], and ranging from 0.30 to 0.67 mS/cm [30]. The ash content of stingless bee honey samples fulfilled the requirement set by the Department of Standards Malaysia, which is not more than 1.00% [9]. Variations in electrical conductivity and ash content among different honey types were suggested due to different botanical sources foraged by the bees [28]. In this study, the values of both electrical conductivity and ash content showed great variability among botanical origins, whereby values found in honeydew honey were significantly higher than blossom honey, which indicate that honeydew honey is richer in both organic and inorganic substances including minerals (Table 5) and organic acids (Table 4).
Diastase is the major enzyme in honey, followed by invertase and glucose oxidase. Enzyme diastase catalyzes the breakdown of starch into maltose which originated from nectar, exudate, and the bees [36]. It is commonly used as an indicator of honey purity and freshness [37]. This enzyme is sensitive to heat and prolonged poor storage conditions; hence, it can be also used to monitor the deterioration of honey [8]. Diastase activity is commonly used in Europe as the indication of freshness of honey. To date, there is no standard for the diastase level of stingless bee honey. The diastase level in stingless bee honey samples reported in this study was rather low, ranging from 1.83 to 3.04 Schade units as compared with the diastase levels reported by a previous study ranging from 0.9 to 23.0 Schade units [38]. Low diastase level in honey samples could be due to acidic pH and hot climate. Diastase activity could be inactivated at pH values lower than 5.3 [39]. Therefore, low diastase level could be explained due to high acidity, as the pH values of stingless bee honey samples in this study range from 3.16 to 3.54. The higher temperature in tropical regions may also diminish the diastase number in honey [8]. Although a previous study showed there were no significant differences in the mean values of diastase activity between honeydew and blossom honeys [21], the honeydew honey was found to have significant higher enzymatic activity than blossom honey in this study.
Hydrogen peroxide is known to be a strong predictor for the antimicrobial activity in honey. To date, there is no legislation available for this parameter. The hydrogen peroxide level of stingless bee honey samples in this study ranged from 111.07 to 192.82 µmol/L which was slightly higher than the hydrogen peroxide level of stingless bee honey reported by another study ranging from 91.50 to 155.80 µmol/L [40]. A significantly higher level of hydrogen peroxide in honeydew honey may be associated with greater antimicrobial activity [41]. Although hydrogen peroxide is an important inhibitory factor of microbial growth, it is not the sole parameter that determines the antimicrobial potency of honey [41]. Although honey with no hydrogen peroxide had little or no antimicrobial activity, some honey samples also failed to exhibit antimicrobial activity despite having a relatively high level of hydrogen peroxide [42]. Studies showed that phenolic compounds and their interaction with hydrogen peroxide are the key factors responsible for the antimicrobial activity of honey [42,43,44]. The information about hydrogen peroxide level in stingless bee honey is still very limited. Hence, more studies are required to investigate the correlation of hydrogen peroxide with the antimicrobial effects of stingless bee honey.
As shown in Table 4, the pH values (from 3.16 to 3.54) of stingless bee honey were in compliance with the Malaysian Standard, in which they should fall between 2.5 and 3.8 [9]. The pH values of stingless bee honey also were similar to the values reported by other studies, in the range of 3.29–3.71 and 3.22–4.03, respectively [45,46]. The high free acidity found in stingless bee honey (60.67–95.33 meq/kg) appears to be a specific characteristic, but stingless bee honey is not included in the international standards for honey. This is due to the presence of organic acids, particularly gluconic acid and acetic acid [22]. Another article stated honeydew honey produced by honey bee showed a significantly higher pH value and acidity than blossom honey [21]. In this study, the free acidity of honeydew honey produced by stingless bee was also significantly higher than the blossom honey, but the pH value was not significantly higher.
Organic acids are minor constituents in honey, accounting for less than 0.5% of the fresh weight of honey. Still, the acidity of honey is strongly related with the organic acids in honey, which are derived from glucose oxidase enzymatic pathway, such as gluconic acid or from microbial fermentation, such as acetic acid [30,47]. Similarly, another study also stated that gluconic acid was the major organic acid in stingless bee honey with the range from 0.07 to 1.48 g/kg [17]. Such findings indicate the importance of gluconic acid in the acidity characteristic of honey. Studies showed that honey with a higher water content can promote the enzymatic activity of glucose oxidase to produce gluconic acid [48,49]. A wide variation of gluconic acid content is possible due to different amount of glucose and glucose oxidase enzymatic activity in honey [17]. The values of acetic acid evaluated in this study were also within the range of values reported by a previous study, ranging from 0.01 to 0.39 g/kg [17]. The higher acidity of honeydew honey in this study could be explained with significant higher level of acetic acid than blossom honey. Excessive acetic acid content is a good indicator of deterioration due to microbial fermentation, although acetic acid is generally found in most honeys [50]. To date, there is no study of normal and fermentation levels of acetic acid in honey. Other than gluconic and acetic acids, the acidity of honey is also contributed to by other organic acids such as succinic, formic, and malic acids [17,51]. Despite the contributions of organic acids to organoleptic properties, especially flavor and to physicochemical properties such as pH and free acidity, information about organic acids of honey is still very limited.
Hydroxymethylfurfural (HMF) content is one of the parameters used to assess the deterioration of honey, because HMF is produced from the degradation of sugars present in honey [8]. HMF is formed slowly during storage or at a faster rate if the honey is heated [15]. The HMF content of stingless bee honey samples were in accordance with the reference value established by the Department of Standards Malaysia for stingless bee honey at maximum level of 30 mg/kg [9]. Such expected values indicating the honey samples were harvested and stored properly, without undergoing any heating process or exposure to high ambient temperature during storage [32].
Table 5 shows that Na is the most abundant element in all stingless bee honey samples, ranging from 223.87 to 326.75 mg/kg, followed by K ranging from 165.27 to 302.40 mg/kg. To our knowledge, this is the first report that Na but not K is the most abundant element detected in honey. Other studies stated K was the most abundant amount of mineral element detected ranging from less than 105.50 to 761.22 mg/kg [32,52]. Na content recorded in the study was found to be within the range values of other studies as well, ranging from 73.00 to 589.46 mg/kg in stingless bee honey [32,52]. Other major mineral elements in stingless bee honey are Ca and Mg, ranging from 37.72 to 70.73 mg/kg and from 32.77 to 59.95 mg/kg, respectively. Other studies also reported similar values for Ca and Mg, ranging from 11.20 to 352.00 mg/kg and from 4.10 to 173.00 mg/kg, respectively [32,52]. Other mineral elements including Fe, Zn, Mn, and Cr were present in low quantities in stingless bee honey [32,52]. Instead of bee types, botanical and geographical origins are the primary factors of mineral content [53]. A previous study also found that the total mineral elements of honeydew honey including K, Mg, and Ca were significantly higher than those in blossom honey [54]. Significantly higher total mineral elements in honeydew honey were also reflected in electrical conductivity and ash content analyses (Table 3) in this study.
Analyzing the data, it can be noted that the data of certain parameters including moisture content, free acidity and diastase activity are considerably different between stingless bee honey and international legislations for honey bee honey specifically Apis mellifera (Codex Alimentarius and Council Directive of the European Union). Such data were not compliant with the established standard limits. Therefore, the current legislations should review the standard limits to fit the majority of studied stingless bee honeys from tropical and subtropical regions.
The total phenolic content of stingless bee honey (56.78–120.06 mg GAE/kg) as shown in Table 6 was generally higher than the previous reported values, from 27.33 to 55.86 mg GAE/kg; hence, higher ferric reducing power was observed in this study ranging from 1.33 to 3.58 mmol Fe(II)/kg compared with 0.54–1.64 mmol Fe(II)/kg [24]. Interestingly, another study reported the antioxidant properties of honey bee honey were found to be higher (total phenolic content: 409.73–475.68 mg GAE/kg; ferric reducing power: 5.00–5.26 mmol Fe(II)/kg) than a stingless bee honey that was produced by Melipona sp. (total phenolic content: 354.03–386.09 mg GAE/kg; ferric reducing power: 4.17–4.39 mmol Fe(II)/kg), but lower than another stingless bee honey which was produced by Hypotrigona sp. (total phenolic content: 375.82–522.91 mg GAE/kg; ferric reducing power: 6.65–6.69 mmol Fe(II)/kg) [22]. Hence, it can be said that there is no definite comparison outcome regarding the antioxidant properties in honey produced by different bee types. Since phenolic compounds in honey originate from plants, the phenolic content is greatly affected by the nectar source harvested by the bees [55]. In this study, honeydew honey was found not only to have significantly higher total phenolic content, and also significantly higher antioxidant capacities including ABTS, DPPH and superoxide radical scavenging activities, peroxyl radical inhibition, iron chelation, and ferric reducing power than blossom honey. Several authors also stated honeydew honey, which is usually darker in color possesses higher antioxidant activities [22,51]. A previous study also found that the mean values of total phenolic content (105.42 mg GAE/100 g), DPPH radical scavenging activity (41.54 mg AAE/100 g), and ferric reducing power (861.06 µmol Fe(II)/100 g) of honeydew honey produced by honey bees were higher than blossom honey, 60.50 mg GAE/100 g, 18.60 mg AAE/100 g, and 354 µmol Fe(II)/100 g, respectively [8]. There are different methods available to assess the antioxidant activity of honey; hence, it is necessary to use several tests to obtain good reliability [22].
PCA was performed to analyze the physicochemical and antioxidant properties of stingless bee honey samples. As shown in Table 7, this analysis managed to highlight the most suitable parameters to be used in the differentiation of honey samples based on botanical origin. Moreover, this statistical analysis was able to differentiate stingless bee honey samples into two distinctive clusters, which are honeydew and blossom origins (Figure 1). Hence, based on the correlation coefficient, it was suggested that parameters including ash content, hydrogen peroxide, free acidity, total mineral elements, K, Mg, Ca, total phenolic compounds, and ferric reducing power are able to differentiate stingless bee honey between honeydew and blossom origins. PCA was also employed in a previous study to differentiate the honeybee honey into honeydew and blossom types [8]. In their study, free acidity, total phenolic compounds, and ferric reducing power plus electrical conductivity, glucose level, and DPPH scavenging activity were suggested to be the suitable parameters.

4. Materials and Methods

4.1. Honey Samples

A total of 23 stingless bee honey samples (Table 8) produced by Heterotrigona itama (Figure S5a) and Geniotrigona thoracica (Figure S5b) were harvested from different areas in the jungle and secondary forest of Southern Negeri Sembilan, Northern Johor, and South Western Pahang in peninsular Malaysia (Figure S6). Raw honey samples were manually filtered and bottled without processing and heat treatment.

4.2. Physicochemical Properties

For 13C NMR analysis, 0.2 moles of each model compound (glucose, fructose, and sucrose) were prepared in deuterated water (D2O). Isoglucose (45% glucose and 55% fructose) and artificial honey (41.27% glucose, 50.79% fructose, and 7.94% sucrose) were prepared with D2O. Each honey sample was prepared by dissolving 200 μL (~260 mg) of honey in 300 μL of D2O. For internal reference, 0.2 moles of 1,4-dioxane (δ 67.19) was used for quantification. All samples were left overnight to fully equilibrate prior to analysis. 13C NMR analysis was conducted on JEOL JNM-ECX400 spectrometer operating at 100 MHz for carbon-13 nuclei. 13C NMR spectra of model compounds were obtained with 90° pulse width 7.25 μs: relaxation delay 2 s, 1000–2000 scans, and four pre-scans. For the isoglucose, artificial honey, and honey samples, the number of scans were increased to 10,000, to achieve a better resolution and sensitivity.
The assignment of 13C NMR chemical shifts of all sugar model compounds including glucose, fructose, and sucrose and their respective tautomers was completed in the chemical shifts reported in the literature [10,11]. The integration of 13C signals allows the direct quantitative determination of tautomers of glucose, fructose, and sucrose. The concentration of each sugar tautomer was calculated based on the signals that were unique for this particular tautomer, and overlapped signals were not used for quantification [10,11]. The quantification of the sugar molecules was achieved upon integration of nonoverlapping signals with the known concentration of 1,4-dioxane, for the internal standard. The applicability of 13C NMR method to quantify sugar molecules was validated by correlating the amount of glucose and fructose in isoglucose, and the amount of glucose, fructose, and sucrose in the artificial honey determined by 13C NMR with the actual weighed amount. The masses of all sugar compounds were calculated in g/100 g, taking into account the moisture content in stingless bee honey.

4.2.1. Color Intensity

The absorbance of filtered honey solution (50%, w/v) was measured at 450 and 720 nm, and the difference in the absorbance readings is expressed as mAU [56].

4.2.2. Total Soluble Solid and Moisture Content

The refractive index of honey was measured against distilled water using a refractometer (Atago, Japan). The reading was recorded in percent °Brix. Moisture content (%, g/100 g) of honey was calculated [57].

4.2.3. Water Activity

The water activity (aw) of honey was measured using a water activity meter (Novasina, Switzerland).

4.2.4. Proline

Honey solution (0.05 g/mL) was mixed with 50% (v/v) formic acid and 3% (w/v) ninhydrin solution, incubated in a boiling water bath for 15 min then at 70 °C for another 10 min, and added with 50% (v/v) 2-propanol. The mixture was left at room temperature before the absorbance was measured at 510 nm. A standard curve was constructed using proline solution (100–500 µg/mL), and the proline level was calculated based on the equation obtained from the standard curve. The final value of proline content (mg/kg) was determined [58].

4.2.5. Electrical Conductivity

According to harmonized methods of the International Honey Commission [56], electrical conductivity of honey solution (0.2 g/mL) was measured using a multiparameter tester (Oakton Instruments, Vernon Hills, IL, USA), and the result was expressed as millisiemens per centimeter (mS/cm).

4.2.6. Ash Content

Two grams of honey was put in a porcelain crucible and dried in an oven at 110 °C for 4 h. The honey was ashed in an electrical furnace (Nabertherm, Germany) at 600 °C for 6 h and weighed. The ash content % (g/100 g) was calculated [58].

4.2.7. Diastase

Honey solution (0.05 g/mL in 100 mM sodium maleate buffer) was incubated at 40 °C water bath for 5 min. An Amylazyme tablet (Megazyme, Ireland) was added into the sample, incubated for another 10 min and added with Trizma base solution, and left at room temperature for 5 min. The mixture was filtered, and the absorbance of sample solution was measured at 590 nm. The diastase number (DN) of honey was calculated.

4.2.8. Hydrogen Peroxide

Each honey sample and hydrogen peroxide standard (7.81–1000 µmol/L) was mixed individually with a working reagent that was composed of ammonium ferrous (II) sulfate, sulfuric acid, sorbitol, and xylenol orange. The mixture was incubated at room temperature for 20 min, and the absorbance was determined at 595 nm. The concentration of hydrogen peroxide (µmol/L) in each honey was calculated with reference to the equation obtained from the standard curve (Thermo Fisher Scientific Inc., Waltham, MA, USA).

4.2.9. pH

A solution containing 10 g of honey dissolved in 75 mL of distilled water was analyzed using a calibrated pH meter (Sartorius, Germany).

4.2.10. Free Acidity

Ten grams of honey was dissolved in 75 mL of distilled water. Titration to pH 8.30 was completed using a standardized 0.1 M NaOH solution. Free acidity (meq/kg) was calculated [58].

4.2.11. d-gluconic Acid

Each honey and standard solution was individually diluted with distilled water, added with buffer (pH 7.6) plus sodium azide (0.02% w/v) and NADP+/ATP solution. Next, 6-phosphogluconate dehydrogenase suspension was added. After 4 min, the absorbance was read at 340 nm. Lastly, gluconate kinase suspension was added. The absorbance was read again at 340 nm after 6 min. The D-gluconic acid level (g/L) was calculated (Megazyme, Ireland). The final value was expressed as g/kg of honey.

4.2.12. Acetic Acid

Each honey and standard solution was individually diluted with distilled water, added with buffer (pH 7.6), sodium azide (0.02% w/v) and NADP/ATP/phosphoenolpyruvate (PEP)/polyvinylpyrrolidone (PVP) solution. Next, coenzyme A solution and a mixture composed of d-lactate dehydrogenase, phosphotransacetylase, and pyruvate kinase were added. After 2 min, the absorbance was read at 340 nm. Lastly, acetate kinase suspension was added. The absorbance was read again at 340 nm after 4 min. The acetic acid level (g/L) was calculated (Megazyme, Ireland). The final value was expressed as g/kg of honey.

4.2.13. Hydroxymethylfurfural (HMF)

Carrez solution I (150 mg/mL potassium ferrocyanide) and Carrez solution II (300 mg/mL zinc acetate) were added into 0.2 g/mL of honey solution. The absorbance was measured at 284 and 336 nm. HMF content of honey (mg/kg) was calculated [59].

4.2.14. Mineral Content

Together with standards, the quantitative determination of minerals, including Na, K, Mg, Ca, Fe, Zn, Mn, and Cr, in each honey was analyzed with flame atomic absorption spectroscopy (AAS) (Agilent Technologies, Santa Clara, CA, USA). Honey was digested with 70% nitric acid and 30% hydrogen peroxide overnight and then heated at 100 °C for 10 min before analysis. The concentration of each mineral element (mg/kg) was obtained [60].

4.3. Antioxidant Properties

4.3.1. Total Phenolic Compounds

Each honey solution (0.2 g/mL) and gallic acid standard solution (200–1000 µg/mL) was mixed with Folin and Ciocalteu’s phenol reagent. After 3 min, 10% (w/v) sodium carbonate solution was added. The reaction mixture was incubated in the dark at room temperature for 90 min. The absorbance of each reaction mixture was read at 725 nm. The final value was calculated and expressed in milligram gallic acid equivalents per kg of sample (mg GAE/kg) [61].

4.3.2. ABTS (2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid)) Radical Scavenging Activity

The ABTS radical cation (ABTS+) solution (absorbance of 0.70 ± 0.02 at 734 nm) was incubated in the dark at room temperature for 12–16 h before use. The ABTS+ solution was added into honey solution (0.2 g/mL), and the reduction of absorbance was determined after 6 min. The free radical scavenging activity (% RSA) was attained [35].

4.3.3. DPPH (2,2-diphenyl-1-picrylhydrazyl) Radical Scavenging Activity

Honey solution (0.2 g/mL) was added with methanolic solution containing DPPH radicals (0.024 mg/mL), and the absorbance of the mixture was read at 517 nm after 15 min in the dark at room temperature. The free radical scavenging activity (% RSA) was attained [62].

4.3.4. Superoxide Anion Radical Scavenging Activity

The reaction started with peroxymonosulfate solution (60 μM) added with superoxide radical solution (nitroblue tetrazolium, 150 μM and NADH, 468 μM) and honey solution (0.2 g/mL). After incubation at room temperature for 5 min, absorbance was read at 560 nm. The free radical scavenging activity (% RSA) was attained [63].

4.3.5. Peroxyl Radical Inhibitory Activity

Honey solution (0.2 g/mL) and Trolox standard solution were added with fluorescein solution, followed with incubation for 30 min at 37 °C. Then, peroxyl radical solution was added. The absorbance of each honey and standard was read with an excitation wavelength of 480 nm and an emission wavelength of 520 nm, every 5 min for a total of 60 min. By using the calculated area under the curve (AUC), the antioxidant activity of honey was calculated and expressed in micromoles Trolox equivalent per liter sample (µmol TE/L) (Cell Biolabs, San Diego, CA, USA).

4.3.6. Iron Chelating Activity

A reaction mixture was prepared after the addition of honey solution (0.2 g/mL) into 0.10 mM ferrous sulfate and 0.25 mM ferrozine. After standing for 10 min at room temperature, the absorbance was read at 562 nm. The chelating activity (%) was attained [63].

4.3.7. Ferric Reducing Activity

FRAP reagent was added into honey solution (0.2 g/mL) and ferrous sulfate standard solution (0.2–1.0 mmol/L), individually. After the reaction mixture was incubated at 37 °C for 4 min, the absorbance was read at 593 nm. The FRAP value (mmol Fe [II]/kg) was attained [64].

4.4. Chemometric Analysis

Each honey was analyzed in triplicates and conducted at 25 °C unless stated otherwise. The data were expressed as means ± standard deviation. Independent t-test was performed to determine the mean value differences at level of significance of 0.05 between honeydew and blossom origins. Principal component analysis (PCA) was employed to interpret interdependence and visualize relatedness between data. The software Microsoft Excel Analyse-it Standard Edition v5.50 was used to perform the statistical analyses.

5. Conclusions

This study demonstrated the potential of 13C NMR spectroscopy to identify and quantify the major sugar molecules in stingless bee honey samples. Consistent major sugar composition observed in honeydew and blossom honey samples could be useful for the identification of stingless bee honey. Furthermore, several parameters including ash content, hydrogen peroxide, free acidity, total mineral elements, K, Mg, Ca, total phenolic compounds, and ferric reducing power were identified with principal component analysis that can differentiate stingless bee honey samples based on botanical origin. The chemometric analysis demonstrated the potentials of using these parameters to evaluate the authenticity of stingless bee honey in Malaysia. This is the first report characterizing stingless bee honey that collected over the recent three years for the differentiation of botanical origins using chemometrics. Obtained data could be useful for legislations for consideration to review the existing standard limits to fit the stingless bee honeys from tropical and subtropical regions.

Supplementary Materials

The following are available. Table S1: 13C chemical shifts of sugars in each model compound. Table S2: Assignment of the carbon resonances and comparison of the measured amount (%, g/100 g) to that of the actual weight amount of each tautomer in the 13C NMR spectra of isoglucose. Table S3: Assignment of the carbon resonances and comparison of the measured amount (%, g/100 g) to that of the actual weight amount of each tautomer in the 13C NMR spectra of artificial honey. Figure S1a: Glucose. Figure S1b: Fructose. Figure S1c: Sucrose. Figure S2: Isoglucose. Figure S3: Artificial honey mixture. Figure S4: 13C NMR spectra of stingless bee honey samples originating from (A–B) honeydew and (C–D) blossom. S5a: Heterotrigona itama. Figure S5b: Geniotrigona thoracica. Figure S6. The location of honey sample collection in peninsular Malaysia.

Author Contributions

Conceptualization, methodology, validation, software, formal analysis, investigation, data curation, writing—original draft preparation, and funding acquisition, W.-J.N.; writing—review and editing, visualization, supervision, project administration, resources, and funding acquisition, N.-W.S., P.A.-C.O. and K.-Y.E.; writing—review and editing, visualization, supervision, project administration, resources, funding acquisition, software, formal analysis, investigation, and data curation, T.-M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by Universiti Tunku Abdul Rahman (UTAR) and partially funded by Eco Bee Shop Sdn. Bhd., Malaysia, grant number 8053/000.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Supplementary Materials here.

Acknowledgments

All technical supports and materials provided by the Faculty of Science, UTAR.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Stingless bee honey samples are available from the authors.

References

  1. Michener, C.D. The meliponini. In Pot-honey; Vit, P., Pedro, S.R.M., Roubik, D., Eds.; Springer: New York, NY, USA, 2013; pp. 3–17. [Google Scholar]
  2. Hrncir, M.; Jarau, S.; Barth, F.G. Stingless bees (Meliponini): Senses and behavior. J. Comp. Physiol. A 2016, 202, 597–601. [Google Scholar] [CrossRef] [PubMed]
  3. Norowi, M.H.; Mohd, F.; Sajap, A.S.; Rosliza, J.; Suri, R. Conservation and sustainable utilization of stingless bees for pollination services in agricultural ecosystems in Malaysia. In Proceedings of the International Seminar on Enhancement of Functional Biodiversity Relevant to Sustainable Food Production in ASPAC, Tsukuba, Japan, 9–11 November 2010. [Google Scholar]
  4. González-Paramás, A.M.; García-Villanova, R.J.; Bárez, J.A.G.; Sánchez, J.S.; Albajar, R.A. Botanical origin of monovarietal dark honeys (from heather, holm oak, pyrenean oak and sweet chestnut) based on their chromatic characters and amino acid profiles. Eur. Food Res. Technol. 2007, 226, 87–92. [Google Scholar] [CrossRef]
  5. Seraglio, S.K.T.; Silva, B.; Bergamo, G.; Brugnerotto, P.; Gonzaga, L.V.; Fett, R.; Costa, A.C.O. An overview of physicochemical characteristics and health-promoting properties of honeydew honey. Int. Food Res. J. 2019, 119, 44–66. [Google Scholar] [CrossRef] [PubMed]
  6. Codex Alimentarius Commission. Standard for Honey. In International Food Standards; Codex Alimentarius Commission: Rome, Italy, 2019; pp. 1912–1981. [Google Scholar]
  7. Council Directive of the European Union. Council Directive 2001/110/EC of 20 December 2001 relating to honey. Off. J. Eur. Communities 2002, 10, 47–52. [Google Scholar]
  8. Bergamo, G.; Seraglio, S.K.T.; Gonzaga, L.V.; Fett, R.; Costa, A.C.O. Physicochemical characteristics of bracatinga honeydew honey and blossom honey produced in the state of Santa Catarina: An approach to honey differentiation. Food Res. Int. 2019, 116, 745–754. [Google Scholar] [CrossRef]
  9. Department of Standards Malaysia. Kelulut (Stingless Bee) Honey-Specification: MS 2683: 2017. In Malaysian Standard; Department of Standards Malaysia: Shah Alam, Malaysia, 2017; pp. 1–33. [Google Scholar]
  10. Kazalaki, A.; Misiak, M.; Spyros, A.; Dais, P. Identification and quantitative determination of carbohydrate molecules in Greek honey by employing 13C NMR spectroscopy. Anal. Methods 2015, 7, 5962–5972. [Google Scholar] [CrossRef]
  11. Mazzoni, V.; Bradesi, P.; Tomi, F.; Casanova, J. Direct qualitative and quantitative analysis of carbohydrate mixtures using 13C NMR spectroscopy: Application to honey. Magn. Reson. Chem. 1997, 35, S81–S90. [Google Scholar] [CrossRef]
  12. Angyal, S.J. The composition of reducing sugars in solution. Adv. Carbohydr. Chem. Biochem. 1984, 42, 5–68. [Google Scholar]
  13. Angyal, S.J. The composition of reducing sugars in solution: Current aspects. Adv. Carbohydr. Chem. Biochem. 1991, 49, 19–35. [Google Scholar]
  14. Angyal, S.J. The composition of reducing sugars in dimethyl sulfoxide solution. Carbohydr. Res. 1994, 1, 1–11. [Google Scholar] [CrossRef]
  15. da Silva, P.M.; Gauche, C.; Gonzaga, L.V.; Costa, A.C.O.; Fett, R. Honey: Chemical composition, stability and authenticity. Food Chem. 2016, 196, 309–323. [Google Scholar] [CrossRef] [PubMed]
  16. Nweze, J.A.; Okafor, J.I.; Nweze, E.I.; Nweze, J.E. Evaluation of physicochemical and antioxidant properties of two stingless bee honeys: A comparison with Apis mellifera honey from Nsukka, Nigeria. BMC Res. Notes 2017, 10, 566. [Google Scholar] [CrossRef] [Green Version]
  17. Shamsudin, S.; Selamat, J.; Sanny, M.; Abd. Razak, S.B.; Jambari, N.N.; Mian, Z.; Khatib, A. Influence of origins and bee species on physicochemical, antioxidant properties and botanical discrimination of stingless bee honey. Int. J. Food Prop. 2019, 22, 239–264. [Google Scholar] [CrossRef] [Green Version]
  18. Fletcher, M.T.; Hungerford, N.L.; Webber, D.; de Jesus, M.C.; Zhang, J.; Stone, I.S.; Blanchfield, J.T.; Zawawi, N. Stingless bee honey, a novel source of trehalulose: A biologically active disaccharide with health benefits. Sci. Rep. 2020, 10, 12128. [Google Scholar] [CrossRef] [PubMed]
  19. Bates, R.B.; Byrne, D.N.; Kane, V.V.; Miller, W.B.; Taylor, S.R. NMR characterization of trehalulose from the excrement of the sweet potato whitefly, Bemisia tabaci. Carbohydr. Res. 1990, 201, 342–345. [Google Scholar] [CrossRef]
  20. Thompson, J.; Robrish, S.A.; Pikis, A.; Brust, A.; Lichtenthaler, F.W. Phosphorylation and metabolism of sucrose and its five linkage-isomeric α-D-glucosyl-D-fructoses by Klebsiella pneumoniae. Carbohydr. Res. 2001, 331, 149–161. [Google Scholar] [CrossRef]
  21. Manzanares, A.B.; García, Z.H.; Galdón, B.R.; Rodríguez, E.R.; Romero, C.D. Differentiation of blossom and honeydew honeys using multivariate analysis on the physicochemical parameters and sugar composition. Food Chem. 2011, 126, 664–672. [Google Scholar] [CrossRef]
  22. Pita-Calvo, C.; Vázquez, M. Differences between honeydew and blossom honeys: A review. Trends Food Sci. Technol. 2017, 59, 79–87. [Google Scholar] [CrossRef]
  23. Vela, L.; de Lorenzo, C.; Perez, R.A. Antioxidant capacity of Spanish honeys and its correlation with polyphenol content and other physicochemical properties. J. Sci. Food Agric. 2007, 87, 1069–1075. [Google Scholar] [CrossRef]
  24. Smanalieva, J.; Senge, B. Analytical and rheological investigations into selected unifloral German honey. Eur. Food Res. Technol. 2009, 229, 107–113. [Google Scholar] [CrossRef]
  25. Manikis, I.; Thrasivoulou, A. The relation of physicochemical characteristics of honey and the crystallization sensitive parameters. Apiacta 2001, 36, 106–112. [Google Scholar]
  26. Ramalhosa, E.; Gomes, T.; Pereira, A.P.; Dias, T.; Estevinho, L.M. Mead production: Tradition versus modernity. Adv. Food Nutr. Res. 2011, 63, 101–118. [Google Scholar] [PubMed]
  27. Souza, E.C.A.; Menezes, C.; Flach, A. Stingless bee honey (Hymenoptera, Apidae, Meliponini): A review of quality control, chemical profile, and biological potential. Apidologie 2021, 52, 113–132. [Google Scholar] [CrossRef]
  28. Suntiparapop, K.; Prapaipong, P.; Chantawannakul, P. Chemical and biological properties of honey from Thai stingless bee (Tetragonula leaviceps). J. Apic. Res. 2012, 51, 45–52. [Google Scholar] [CrossRef]
  29. Kek, S.P.; Chin, N.L.; Yusof, Y.A.; Tan, S.W.; Chua, L.S. Classification of entomological origin of honey based on its physicochemical and antioxidant properties. Int. J. Food Prop. 2017, 20 (Suppl. S3), S2723–S2738. [Google Scholar] [CrossRef]
  30. de Sousa, J.M.B.; de Souza, E.L.; Marques, G.; de Toledo Benassi, M.; Gullón, B.; Pintado, M.M.; Magnani, M. Sugar profile, physicochemical and sensory aspects of monofloral honeys produced by different stingless bee species in Brazilian semi-arid region. LWT—Food Sci. Technol. 2016, 65, 645–651. [Google Scholar] [CrossRef] [Green Version]
  31. Silvano, M.F.; Varela, M.S.; Palacio, M.A.; Ruffinengo, S.; Yamul, D.K. Physicochemical parameters and sensory properties of honeys from Buenos Aires region. Food Chem. 2014, 152, 500–507. [Google Scholar] [CrossRef]
  32. Biluca, F.C.; Braghini, F.; Gonzaga, L.V.; Costa, A.C.O.; Fett, R. Physicochemical profiles, minerals and bioactive compounds of stingless bee honey (Meliponinae). J. Food Compost. Anal. 2016, 50, 61–69. [Google Scholar] [CrossRef]
  33. Czipa, N.; Borbely, M.; Gyori, Z. Proline content of different honey types. Acta Aliment. 2012, 41, 26–32. [Google Scholar] [CrossRef]
  34. Wen, Y.Q.; Zhang, J.; Li, Y.; Chen, L.; Zhao, W.; Zhou, J.; Jin, Y. Characterization of Chinese unifloral honeys based on proline and phenolic content as markers of botanical origin, using multivariate analysis. Molecules 2017, 22, 735. [Google Scholar] [CrossRef] [Green Version]
  35. Moniruzzaman, M.; Khalil, M.I.; Sulaiman, S.A.; Gan, S.H. Physicochemical and antioxidant properties of Malaysian honeys produced by Apis cerana, Apis dorsata and Apis mellifera. BMC Complement. Altern. Med. 2013, 13, 43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Ahmed, M.; Djebli, N.; Aissat, S.; Bacha, S.; Meslem, A.; Khiati, B. Synergistic inhibition of natural honey and potato starch and their correlation with diastase number and sugar content against Klebsiella pneumoniae ATCC 27736. Nat. Prod. Chem. Res. 2012, 1, 102. [Google Scholar] [CrossRef] [Green Version]
  37. Machado De-Melo, A.A.; Almeida-Muradian, L.B.D.; Sancho, M.T.; Pascual-Maté, A. Composition and properties of Apis mellifera honey: A review. J. Apic. Res. 2018, 57, 5–37. [Google Scholar] [CrossRef]
  38. Souza, B.A.; Roubik, D.W.; Barth, O.M.; Heard, T.A.; Enríquez, E.; Carvalho, C.; Villas-Bôas, J.; Marchini, L.; Locatelli, J.; Persano-Oddo, L.; et al. Composition of stingless bee honey: Setting quality standards. Interciencia 2006, 31, 867–875. [Google Scholar]
  39. Babacan, S.; Rand, A.G. Characterization of honey amylase. J. Food Sci. 2007, 72, C050–C055. [Google Scholar] [CrossRef] [PubMed]
  40. Massaro, C.F.; Shelley, D.; Heard, T.A.; Brooks, P. In vitro antibacterial phenolic extracts from “sugarbag” pot-honeys of Australian stingless bees (Tetragonula carbonaria). J. Agric. Food Chem. 2014, 62, 12209–12217. [Google Scholar] [CrossRef] [PubMed]
  41. Ng, W.J.; Sit, N.W.; Ooi, P.A.C.; Ee, K.Y.; Lim, T.M. The antibacterial potential of honeydew honey produced by stingless bee (Heterotrigona itama) against antibiotic resistant bacteria. Antibiotics 2020, 9, 871. [Google Scholar] [CrossRef]
  42. Chen, C.; Campbell, L.; Blair, S.E.; Carter, D.A. The effect of standard heat and filtration processing procedures on antimicrobial activity and hydrogen peroxide levels in honey. Front. Microbiol. 2012, 3, 265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Brudzynski, K.; Abubaker, K.; Castle, A. Re-examining the role of hydrogen peroxide in bacteriostatic and bactericidal activities of honey. Front. Microbiol. 2011, 2, 213. [Google Scholar] [CrossRef] [Green Version]
  44. Bucekova, M.; Sojka, M.; Valachova, I.; Martinotti, S.; Ranzato, E.; Szep, Z.; Majtan, V.; Klaudiny, J.; Majtan, J. Bee-derived antibacterial peptide, defensin-1, promotes wound re-epithelialisation in vitro and in vivo. Sci. Rep. 2017, 7, 7340. [Google Scholar] [CrossRef] [Green Version]
  45. Chan, B.K.; Haron, H.; Talib, R.A.; Subramaniam, P. Physical properties, antioxidant content and anti-oxidative activities of Malaysian stingless kelulut (Trigona spp.) honey. J. Agric. Sci. 2017, 9, 32–40. [Google Scholar]
  46. Selvaraju, K.; Vikram, P.; Soon, J.M.; Krishnan, K.T.; Mohammed, A. Melissopalynological, physicochemical and antioxidant properties of honey from West Coast of Malaysia. J. Food Sci. Technol. 2019, 56, 2508–2521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Dardon, M.J.; Maldonado-Aguilera, C.; Enríquez, E. The pot-honey of Guatemalan bees. In Pot-Honey; Vit, P., Pedro, S.R.M., Roubik, D., Eds.; Springer: New York, NY, USA, 2013; pp. 395–408. [Google Scholar]
  48. Olaitan, P.B.; Adeleke, O.E.; Iyabo, O.O. Honey: A reservoir for microorganisms and an inhibitory agent for microbes. Afr. Health Sci. 2007, 7, 159–165. [Google Scholar]
  49. Kretavičius, J.; Kurtinaitienė, B.; Račys, J.; Čeksterytė, V. Inactivation of glucose oxidase during heat-treatment de-crystallization of honey. Zemdirb. Agric. 2010, 97, 115–122. [Google Scholar]
  50. Mato, I.; Huidobro, J.F.; Simal-Lozano, J.; Sancho, M.T. Significance of nonaromatic organic acids in honey. J. Food Prot. 2003, 66, 2371–2376. [Google Scholar] [CrossRef]
  51. Tezcan, F.; Kolayli, S.; Sahin, H.; Ulusoy, E.; Erim, F.B. Evaluation of organic acid, saccharide composition and antioxidant properties of some authentic Turkish honeys. J. Food Nutr. Res. 2011, 50, 33–40. [Google Scholar]
  52. Cheng, M.Z.S.Z.; Ismail, M.; Chan, K.W.; Ooi, D.J.; Ismail, N.; Zawawi, N.; Mohd Esa, N. Comparison of sugar content, mineral elements and antioxidant properties of Heterotrigona itama honey from suburban and forest in Malaysia. Malays. J. Med. Health Sci. 2019, 15, 104–112. [Google Scholar]
  53. Kek, S.P.; Chin, N.L.; Tan, S.W.; Yusof, Y.A.; Chua, L.S. Classification of honey from its bee origin via chemical profiles and mineral content. Food Anal. Methods 2017, 10, 19–30. [Google Scholar] [CrossRef]
  54. Olga, E.; María, F.G.; Carmen, S.M. Differentiation of blossom honey and honeydew honey from Northwest Spain. Agriculture 2012, 2, 25–37. [Google Scholar] [CrossRef] [Green Version]
  55. Truchado, P.; Vit, P.; Heard, T.A.; Tomás-Barberán, F.A.; Ferreres, F. Determination of interglycosidic linkages in O-glycosyl flavones by high-performance liquid chromatography/photodiode-array detection coupled to electrospray ionization ion trap mass spectrometry. Its application to Tetragonula carbonaria honey from Australia. Rapid Commun. Mass Spectrum. 2015, 29, 948–954. [Google Scholar]
  56. Beretta, G.; Granata, P.; Ferrero, M.; Orioli, M.; Facino, R.M. Standardization of antioxidant properties of honey by a combination of spectrophotometric/fluorimetric assays and chemometrics. Anal. Chim. Acta 2005, 533, 185–191. [Google Scholar] [CrossRef]
  57. Association of Official Analytical Chemists (AOAC). Official Methods of Analysis of the Association of Official Analytical Chemists. J. AOAC Int. 1990, 934, 1–2. [Google Scholar]
  58. Bogdanov, S.; Martin, P.; Lullmann, C. Harmonised Methods of the International Honey Commission; Swiss Bee Research Centre, FAM: Liebefeld, Switzerland, 2002; pp. 1–63. [Google Scholar]
  59. Khalil, M.I.; Sulaiman, S.A.; Gan, S.H. High 5-hydroxymethylfurfural concentrations are found in Malaysian honey samples stored for more than one year. Food Chem. Toxicol. 2010, 48, 2388–2392. [Google Scholar] [CrossRef]
  60. Vanhanen, L.P.; Emmertz, A.; Savage, G.P. Mineral analysis of mono-floral New Zealand honey. Food Chem. 2011, 128, 236–240. [Google Scholar] [CrossRef]
  61. Khalil, M.; Moniruzzaman, M.; Boukraa, L.; Benhanifia, M.; Islam, M.; Sulaiman, S.A.; Gan, S.H. Physicochemical and antioxidant properties of Algerian honey. Molecules 2012, 17, 11199–11215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Ferreira, I.C.; Aires, E.; Barreira, J.C.; Estevinho, L.M. Antioxidant activity of Portuguese honey samples: Different contributions of the entire honey and phenolic extract. Food Chem. 2009, 114, 1438–1443. [Google Scholar] [CrossRef]
  63. Chai, T.; Mohan, M.; Ong, H.; Wong, F. Antioxidant, iron-chelating and anti-glucosidase activities of Typha domingensis Pers (Typhaceae). Trop. J. Pharm. Res. 2014, 13, 67–72. [Google Scholar] [CrossRef] [Green Version]
  64. Benzie, I.F.; Strain, J.J. Ferric reducing/antioxidant power assay: Direct measure of total antioxidant activity of biological fluids and modified version for simultaneous measurement of total antioxidant power and ascorbic acid concentration. Meth. Enzymol. 1999, 299, 15–27. [Google Scholar]
Figure 1. (a) Plot of principal component loadings stingless bee blossom and honeydew honey samples and the descriptors including physicochemical and antioxidant properties. (b) Plot of principal component loadings stingless bee blossom and honeydew honey samples and the descriptors including physicochemical and antioxidant properties.
Figure 1. (a) Plot of principal component loadings stingless bee blossom and honeydew honey samples and the descriptors including physicochemical and antioxidant properties. (b) Plot of principal component loadings stingless bee blossom and honeydew honey samples and the descriptors including physicochemical and antioxidant properties.
Molecules 26 07628 g001
Table 1. Quantification of sugar tautomers present in stingless bee honey samples using 13C NMR spectroscopy.
Table 1. Quantification of sugar tautomers present in stingless bee honey samples using 13C NMR spectroscopy.
SugarHoneydew HoneyBlossom Honey
Tautomer %Average Integration Valueg/100 gTautomer %Average Integration Valueg/100 g
α-d-glucopyranose (α-GP)41.700.282514.41 ± 0.4240.500.247514.45 ± 0.31
β-d-glucopyranose (β-GP)58.300.395020.99 ± 0.6259.500.368820.96 ± 0.52
Total glucose 35.39 ± 0.20 35.36 ± 0.76
α-d-fructopyranose (α-FP)2.590.01881.02 ± 0.054.170.02781.60 ± 0.01
β-d-fructopyranose (β-FP)69.590.505025.99 ± 0.4068.990.460026.19 ± 0.14
α-d-fructofuranose (α-FF)6.010.04362.29 ± 0.026.200.04142.39 ± 0.01
β-d-fructofuranose (β-FF)21.810.15838.32 ± 0.1220.650.13788.08 ± 0.04
Total fructose 37.61 ± 0.20 38.25 ± 0.21
Total sucrose----0.00380.40 ± 0.57
Fructose-to-glucose ratio (F/G) 1.06 ± 0.01 1.08 ± 0.03
Glucose-to-moisture ratio (G/M) 1.31 ± 0.01 1.38 ± 0.03
Each total content of glucose and fructose between honeydew honey and blossom honey was not significantly different at p < 0.05.
Table 2. The maturity parameters of honey samples.
Table 2. The maturity parameters of honey samples.
Honey SampleColor Intensity (mAU)Moisture
(%, g/100 g)
Water ActivityTotal Soluble Solid
(°Brix)
Proline (mg/kg)
S1180 ± 027.30 ± 00.62 ± 072.70 ± 0537.13 ± 6.86
S2190 ± 027.13 ± 0.060.60 ± 072.87 ± 0.06513.13 ± 5.00
S3200 ± 026.40 ± 00.62 ± 073.60 ± 0522.93 ± 20.01
S4190 ± 026.37 ± 0.060.63 ± 073.63 ± 0.06542.37 ± 37.50
S5190 ± 026.33 ± 0.060.63 ± 073.67 ± 0.06547.43 ± 16.01
S6150 ± 027.87 ± 0.060.62 ± 072.13 ± 0.06514.50 ± 3.90
S7140 ± 027.40 ± 00.65 ± 072.60 ± 0454.80 ± 11.39
S8150 ± 027.43 ± 0.060.65 ± 072.57 ± 0.06546.73 ± 10.83
S960 ± 023.50 ± 0.100.54 ± 076.50 ± 0.10577.00 ± 2.00
S1080 ± 022.57 ± 0.060.54 ± 077.43 ± 0.06546.00 ± 2.69
S1160 ± 022.50 ± 00.55 ± 077.50 ± 0578.37 ± 4.47
S1270 ± 026.60 ± 00.54 ± 073.40 ± 0588.80 ± 35.76
S1360 ± 026.73 ± 0.060.57 ± 073.27 ± 0.06514.23 ± 14.70
S1450 ± 024.57 ± 0.060.64 ± 075.43 ± 0.06525.27 ± 7.28
S1550 ± 024.30 ± 00.63 ± 075.70 ± 0525.70 ± 4.62
S1670 ± 024.20 ± 00.63 ± 075.80 ± 0525.10 ± 7.82
S17250 ± 026.20 ± 00.64 ± 073.80 ± 0670.27 ± 6.78
S18250 ± 027.20 ± 00.67 ± 072.80 ± 0598.50 ± 3.50
S19240 ± 027.10 ± 00.67 ± 072.87 ± 0.06612.40 ± 4.25
S20280 ± 027.30 ± 00.67 ± 072.70 ± 0627.57 ± 7.10
S21220 ± 027.37 ± 0.060.67 ± 072.63 ± 0.06510.67 ± 5.86
S22250 ± 028.00 ± 00.69 ± 072.00 ± 0510.47 ± 1.50
S23210 ± 026.97 ± 0.060.67 ± 073.03 ± 0.06535.70 ± 4.96
Average156.09 ± 76.5926.14 ± 1.640.62 ± 0.0573.85 ± 1.64548.94 ± 47.74
Honeydew173.75 ± 22.2327.03 ± 0.56 *0.63 ± 0.0272.97 ± 0.56 *522.38 ± 32.61 *
Blossom146.67 ± 92.4425.67 ± 1.830.62 ± 0.0674.32 ± 1.83563.11 ± 48.74
*—Significant difference between honeydew honey and blossom honey at p < 0.05.
Table 3. The purity parameters of honey samples.
Table 3. The purity parameters of honey samples.
Honey SampleElectrical Conductivity (mS/cm)Ash Content
(%, g/100 g)
Diastase
(Schade Unit/g)
Hydrogen Peroxide
(µmol/L)
S10.36 ± 00.18 ± 0.013.04 ± 0.08177.58 ± 19.80
S20.39 ± 00.16 ± 0.012.98 ± 0.05185.43 ± 7.29
S30.42 ± 00.15 ± 02.45 ± 0.03187.05 ± 20.39
S40.45 ± 00.13 ± 0.012.51 ± 0.09177.19 ± 18.63
S50.41 ± 00.13 ± 0.012.71 ± 0.04189.77 ± 18.09
S60.52 ± 00.13 ± 02.22 ± 0.08184.07 ± 13.14
S70.51 ± 00.13 ± 0.012.05 ± 0.02192.82 ± 20.58
S80.52 ± 00.14 ± 0.012.34 ± 0.07174.18 ± 6.22
S90.38 ± 00.04 ± 02.63 ± 0.59111.40 ± 1.77
S100.35 ± 00.04 ± 02.09 ± 0.02111.97 ± 1.82
S110.29 ± 00.04 ± 0.012.25 ± 0.03118.20 ± 10.78
S120.29 ± 00.04 ± 0.012.15 ± 0.04112.40 ± 2.50
S130.41 ± 00.04 ± 0.012.30 ± 0.04112.95 ± 3.49
S140.42 ± 00.04 ± 0.012.17 ± 0.18119.53 ± 8.88
S150.45 ± 00.04 ± 0.012.39 ± 0.02114.15 ± 9.31
S160.45 ± 00.04 ± 02.38 ± 0.03111.07 ± 11.28
S170.38 ± 00.05 ± 0.011.83 ± 0.15168.37 ± 21.51
S180.35 ± 00.05 ± 0.012.23 ± 0.35146.17 ± 30.41
S190.30 ± 00.05 ± 0.012.18 ± 0.14166.73 ± 12.74
S200.28 ± 00.05 ± 0.012.16 ± 0.14160.17 ± 17.36
S210.36 ± 00.05 ± 0.012.43 ± 0.02169.63 ± 4.57
S220.38 ± 00.05 ± 0.012.36 ± 0.03178.61 ± 2.01
S230.38 ± 00.05 ± 02.29 ± 0.08162.77 ± 21.68
Average0.39 ± 0.070.08 ± 0.052.35 ± 0.31153.57 ± 33.12
Honeydew0.45 ± 0.06 *0.14 ± 0.02*2.54 ± 0.34 *183.51 ± 15.11 *
Blossom0.36 ± 0.050.04 ± 0.012.26 ± 0.24137.61 ± 28.74
*—Significant difference between honeydew honey and blossom honey at p < 0.05.
Table 4. The deterioration parameters of honey samples.
Table 4. The deterioration parameters of honey samples.
Honey SamplepHFree Acidity (meq/kg)Gluconic Acid (g/kg)Acetic Acid (g/kg)Hydroxymethylfurfural
(mg/kg)
S13.26 ± 0.0289.33 ± 0.580.34 ± 0.010.09 ± 0.0119.69 ± 0.04
S23.26 ± 0.0187.33 ± 0.580.31 ± 0.010.11 ± 0.0219.41 ± 0.56
S33.18 ± 0.0285.00 ± 1.000.38 ± 0.060.13 ± 0.0319.48 ± 0.36
S43.16 ± 0.0293.00 ± 2.000.34 ± 0.010.12 ± 0.0417.31 ± 0.81
S53.21 ± 0.0290.33 ± 0.580.37 ± 0.030.09 ± 0.0316.62 ± 0.22
S63.54 ± 0.0194.67 ± 0.580.71 ± 0.060.09 ± 0.0311.45 ± 0.37
S73.54 ± 0.0192.67 ± 0.580.62 ± 0.030.08 ± 0.0214.85 ± 0.32
S83.51 ± 0.0195.33 ± 0.580.62 ± 0.040.09 ± 0.0315.64 ± 0.11
S93.46 ± 0.0172.00 ± 1.000.57 ± 0.010.09 ± 0.0319.87 ± 0.46
S103.51 ± 0.0170.67 ± 0.580.54 ± 0.050.16 ± 0.0519.17 ± 0.38
S113.28 ± 0.0172.00 ± 00.58 ± 0.030.12 ± 0.0217.22 ± 0.72
S123.27 ± 0.0178.00 ± 00.45 ± 00.09 ± 0.0227.10 ± 0.92
S133.22 ± 0.0175.00 ± 00.35 ± 0.010.08 ± 0.0113.13 ± 1.51
S143.25 ± 0.0262.00 ± 2.000.56 ± 0.050.10 ± 0.0213.54 ± 0.44
S153.17 ± 0.0260.67 ± 0.580.67 ± 0.020.10 ± 0.0214.22 ± 0.41
S163.22 ± 0.0361.67 ± 0.580.72 ± 0.060.13 ± 0.0212.71 ± 0.55
S173.21 ± 080.67 ± 0.580.34 ± 0.010.02 ± 0.0120.24 ± 0.46
S183.25 ± 0.0181.67 ± 0.580.35 ± 0.030.05 ± 0.0220.39 ± 0.46
S193.21 ± 0.0181.67 ± 0.580.43 ± 0.040.03 ± 0.0119.53 ± 0.27
S203.28 ± 0.0185.00 ± 00.25 ± 0.010.03 ± 0.0221.10 ± 1.28
S213.51 ± 086.67 ± 0.580.34 ± 0.020.04 ± 0.0124.54 ± 0.04
S223.51 ± 091.00 ± 1.000.30 ± 0.030.03 ± 0.0219.69 ± 0.19
S233.50 ± 091.00 ± 00.65 ± 0.050.04 ± 0.0118.88 ± 1.04
Average3.33 ± 0.1481.62 ± 10.700.47 ± 0.150.08 ± 0.0418.08 ± 3.74
Honeydew3.33 ± 0.1690.96 ± 3.54 *0.46 ± 0.150.10 ± 0.03 *16.81 ± 2.74 *
Blossom3.32 ± 0.1376.64 ± 9.880.47 ± 0.150.07 ± 0.0518.76 ± 4.05
*—Significant different between honeydew honey and blossom honey at p < 0.05.
Table 5. (1) The mineral elements profile of honey samples. (2) The mineral elements profile of honey samples.
Table 5. (1) The mineral elements profile of honey samples. (2) The mineral elements profile of honey samples.
(1)
Honey SampleNa
(mg/kg)
K
(mg/kg)
Mg
(mg/kg)
Ca
(mg/kg)
Fe
(mg/kg)
S1283.40 ± 6.01298.27 ± 11.0650.51 ± 0.6467.27 ± 2.8512.17 ± 1.13
S2295.42 ± 1.95285.92 ± 7.1755.39 ± 1.2970.24 ± 2.8911.82 ± 1.51
S3326.75 ± 22.74269.35 ± 40.9158.19 ± 0.4269.76 ± 2.9612.86 ± 0.35
S4305.63 ± 14.03296.93 ± 6.3959.95 ± 0.3865.54 ± 3.0514.29 ± 0.89
S5312.14 ± 25.09274.47 ± 28.4654.54 ± 0.4970.73 ± 0.9612.77 ± 1.22
S6300.13 ± 9.81302.40 ± 4.2752.64 ± 0.8666.43 ± 3.8313.77 ± 1.97
S7308.57 ± 10.66285.60 ± 14.2052.99 ± 0.4469.03 ± 2.4312.26 ± 0.86
S8316.87 ± 32.21273.73 ± 59.2654.29 ± 1.4669.77 ± 2.3612.50 ± 0.56
S9256.60 ± 27.52190.69 ± 11.0042.02 ± 1.8947.23 ± 5.7212.34 ± 0.80
S10242.68 ± 12.61178.61 ± 14.9135.95 ± 3.0337.72 ± 0.7910.41 ± 0.74
S11223.87 ± 16.54167.20 ± 11.0232.77 ± 1.0438.24 ± 2.7010.73 ± 1.50
S12263.65 ± 11.12166.67 ± 9.6234.57 ± 3.5947.45 ± 2.1210.68 ± 1.21
S13243.53 ± 48.50165.27 ± 14.5542.09 ± 5.2145.45 ± 4.3910.90 ± 1.57
S14286.20 ± 6.91191.37 ± 6.2542.53 ± 5.1253.90 ± 7.3712.73 ± 0.81
S15277.03 ± 23.49204.00 ± 16.5247.41 ± 4.0458.13 ± 2.7512.50 ± 0.66
S16283.70 ± 5.79179.83 ± 17.8244.78 ± 2.6963.27 ± 5.1611.87 ± 0.61
S17278.83 ± 12.71223.65 ± 14.9844.58 ± 1.5161.75 ± 2.0811.26 ± 0.92
S18278.59 ± 28.74232.70 ± 16.2246.36 ± 1.3563.73 ± 0.5011.38 ± 0.78
S19260.30 ± 43.61209.05 ± 18.2044.82 ± 2.0663.39 ± 1.1710.71 ± 0.93
S20263.88 ± 28.77224.23 ± 13.8940.25 ± 1.8563.99 ± 2.9312.09 ± 0.70
S21275.60 ± 25.84189.27 ± 12.6646.13 ± 5.1670.30 ± 2.0912.33 ± 1.46
S22274.60 ± 7.45184.07 ± 11.6050.83 ± 4.4969.33 ± 0.9010.57 ± 1.10
S23292.27 ± 6.55179.93 ± 1.1547.02 ± 2.3463.77 ± 6.7711.57 ± 0.80
Average280.45 ± 31.09224.92 ± 51.1746.98 ± 7.5660.71 ± 10.7011.94 ± 1.34
Honeydew306.11 ± 19.79 *285.83 ± 26.51 *54.81 ± 3.00 *68.60 ± 2.96 *12.80 ± 1.26 *
Blossom266.76 ± 27.16192.44 ± 24.0042.81 ± 5.6856.51 ± 10.9811.47 ± 1.14
(2)
Honey SampleZn
(mg/kg)
Mn
(mg/kg)
Cr
(mg/kg)
Cu + Al (mg/kg)Total Mineral Elements
(mg/kg)
S13.40 ± 0.260.44 ± 0.140.42 ± 0.08<LOQ715.88 ± 7.39
S23.74 ± 0.320.29 ± 0.210.66 ± 0.38<LOQ723.47 ± 5.73
S33.78 ± 0.380.59 ± 0.250.27 ± 0.20<LOQ741.56 ± 38.62
S43.63 ± 1.100.75 ± 0.160.44 ± 0.10<LOQ747.15 ± 16.63
S52.66 ± 0.580.63 ± 0.040.27 ± 0.20<LOQ728.21 ± 54.06
S63.88 ± 0.300.63 ± 0.060.63 ± 0.32<LOQ740.52 ± 3.59
S74.25 ± 0.280.72 ± 0.280.53 ± 0.40<LOQ733.95 ± 2.96
S84.14 ± 0.460.67 ± 0.150.60 ± 0.26<LOQ732.56 ± 26.36
S93.47 ± 0.650.62 ± 0.130.61 ± 0.11<LOQ553.58 ± 13.19
S103.41 ± 0.190.54 ± 0.070.56 ± 0.12<LOQ509.88 ± 9.36
S112.27 ± 0.850.42 ± 0.140.45 ± 0.06<LOQ475.94 ± 24.77
S122.98 ± 0.110.40 ± 0.130.51 ± 0.07<LOQ526.91 ± 15.85
S132.60 ± 0.380.52 ± 0.060.67 ± 0.11<LOQ511.02 ± 44.73
S143.03 ± 0.390.68 ± 0.100.82 ± 0.07<LOQ591.27 ± 13.90
S152.71 ± 0.330.69 ± 0.120.42 ± 0.12<LOQ602.90 ± 42.12
S163.31 ± 0.440.67 ± 0.120.61 ± 0.36<LOQ588.03 ± 20.94
S172.39 ± 0.780.45 ± 0.190.53 ± 0.25<LOQ623.44 ± 17.96
S182.00 ± 0.920.37 ± 0.050.63 ± 0.12<LOQ635.77 ± 31.56
S191.81 ± 0.190.39 ± 0.180.48 ± 0.09<LOQ590.96 ± 54.17
S201.76 ± 0.240.42 ± 0.270.54 ± 0.04<LOQ607.17 ± 19.59
S212.94 ± 0.240.60 ± 0.170.73 ± 0.16<LOQ597.89 ± 46.93
S223.32 ± 0.490.65 ± 0.050.59 ± 0.08<LOQ593.96 ± 9.66
S233.09 ± 0.170.80 ± 00.79 ± 0.17<LOQ599.24 ± 8.09
Average3.07 ± 0.810.56 ± 0.190.55 ± 0.21-629.19 ± 88.03
Honeydew3.68 ± 0.65 *0.59 ± 0.210.48 ± 0.27-732.91 ± 23.96 *
Blossom2.74 ± 0.700.55 ± 0.170.60 ± 0.17-573.86 ± 51.74
*—Significantly different between honeydew honey and blossom honey at p < 0.05. LOQ—Limits of quantitation.
Table 6. The antioxidant properties of honey samples.
Table 6. The antioxidant properties of honey samples.
Honey SampleTotal Phenolic Compounds
(mg GAE/kg)
ABTS Radical Scavenging Activity
(%)
DPPH Radical Scavenging Activity
(%)
Superoxide Radical Scavenging Activity
(%)
Peroxyl Radical Inhibition (µmol TE/g)Iron Chelation (%)Ferric Reducing Power
(mmol Fe(II)/kg)
S1105.60 ± 1.6359.13 ± 1.8734.83 ± 0.5676.60 ± 0.365.47 ± 0.0417.38 ± 1.533.20 ± 0.02
S2105.60 ± 4.3763.28 ± 0.4136.03 ± 0.3377.51 ± 0.345.61 ± 0.3014.29 ± 0.443.18 ± 0.03
S3120.06 ± 1.2970.90 ± 0.3636.26 ± 0.2576.00 ± 0.105.81 ± 0.1616.02 ± 0.813.31 ± 0.02
S4109.74 ± 1.0368.17 ± 0.8234.52 ± 1.2775.30 ± 0.265.61 ± 0.1216.08 ± 0.473.29 ± 0.02
S5113.59 ± 1.6968.69 ± 0.6637.37 ± 1.5076.50 ± 0.406.12 ± 0.3316.18 ± 0.703.58 ± 0.31
S695.39 ± 2.6355.68 ± 0.7035.44 ± 0.4073.98 ± 0.035.88 ± 0.3222.51 ± 0.933.06 ± 0.04
S783.28 ± 16.3658.00 ± 0.7035.25 ± 0.6976.02 ± 0.795.48 ± 0.2421.74 ± 1.872.96 ± 0.05
S899.51 ± 1.6261.33 ± 0.5835.49 ± 0.9376.92 ± 0.245.70 ± 0.3327.29 ± 1.522.87 ± 0.05
S965.58 ± 2.0353.89 ± 0.4536.14 ± 0.4270.80 ± 0.804.58 ± 0.1816.43 ± 1.561.33 ± 0.01
S1073.54 ± 2.1162.07 ± 0.5332.62 ± 2.1771.07 ± 0.674.90 ± 0.1713.08 ± 0.911.71 ± 0.01
S1162.25 ± 0.9453.44 ± 0.6835.48 ± 2.2069.73 ± 0.705.19 ± 0.2710.83 ± 0.261.64 ± 0.03
S1275.51 ± 2.2554.33 ± 0.3134.33 ± 0.3174.13 ± 2.805.18 ± 0.208.26 ± 1.231.63 ± 0.04
S1361.01 ± 2.2244.77 ± 0.6634.67 ± 0.5173.87 ± 2.345.97 ± 0.2010.54 ± 0.561.56 ± 0.04
S1457.35 ± 0.9851.28 ± 0.2433.38 ± 0.7274.06 ± 1.435.60 ± 0.1712.91 ± 1.671.57 ± 0.05
S1556.78 ± 2.2456.38 ± 0.7032.84 ± 1.7475.03 ± 0.454.85 ± 0.2710.18 ± 1.021.83 ± 0.01
S1659.32 ± 0.8551.82 ± 0.3533.09 ± 0.5172.45 ± 1.334.44 ± 0.1011.62 ± 0.691.56 ± 0.04
S1797.04 ± 1.4558.37 ± 0.4731.59 ± 0.8071.13 ± 0.815.61 ± 0.2412.37 ± 0.681.62 ± 0.08
S1899.94 ± 2.6960.44 ± 0.8733.44 ± 1.5575.07 ± 0.505.60 ± 0.1814.79 ± 0.262.25 ± 0.05
S1994.67 ± 2.0556.01 ± 0.2231.56 ± 0.6979.36 ± 0.555.65 ± 0.1119.15 ± 0.481.90 ± 0.07
S20101.13 ± 1.2162.00 ± 0.2632.83 ± 1.8073.73 ± 1.425.66 ± 0.0611.78 ± 0.382.20 ± 0.02
S2178.22 ± 0.5966.28 ± 0.4833.16 ± 0.6977.22 ± 0.705.64 ± 0.2012.89 ± 0.252.00 ± 0.01
S2288.14 ± 0.8566.74 ± 0.8232.17 ± 0.6978.89 ± 1.405.51 ± 0.1115.07 ± 0.612.29 ± 0.08
S2377.11 ± 2.6856.92 ± 0.5831.90 ± 1.6378.90 ± 1.505.47 ± 0.2214.85 ± 1.572.25 ± 0.34
Average86.10 ± 19.6959.13 ± 6.3334.10 ± 1.8974.97 ± 2.805.46 ± 0.4515.05 ± 4.432.29 ± 0.71
Honeydew104.09 ± 12.06 *63.15 ± 5.38 *35.65 ± 1.12 *76.10 ± 1.08 *5.71 ± 0.30 *18.94 ± 4.36 *3.18 ± 0.23 *
Blossom76.51 ± 15.8956.98 ± 5.7733.28 ± 1.6974.36 ± 3.235.32 ± 0.4712.98 ± 2.791.82 ± 0.30
*—Significant different between honeydew honey and blossom honey at p < 0.05.
Table 7. Factor loadings for parameters of stingless bee honey samples.
Table 7. Factor loadings for parameters of stingless bee honey samples.
Data Related to Figure 1aData Related to Figure 1b
Principal Component (PC) Number1212
Eigenvalue11.6935.2206.9520.666
% variance39.00%17.40%77.20%7.40%
Component score correlation
Color intensity−0.634−0.698
Moisture content−0.764−0.407
Water activity−0.560−0.390
Total soluble solid0.7630.408
Proline0.285−0.608
Electrical conductivity−0.4650.720
Ash content−0.8220.334−0.8910.310
Diastase−0.3300.206
Hydrogen peroxide−0.886−0.147−0.878−0.324
pH−0.1070.179
Free acidity−0.842−0.201−0.811−0.438
Gluconic acid0.2880.686
Acetic acid0.1640.715
HMF−0.040−0.666
Na−0.7100.268
K−0.8170.235−0.8920.307
Mg−0.8760.245−0.8760.158
Ca−0.863−0.127−0.810−0.229
Fe−0.4520.413
Zn−0.3750.654
Mn−0.1370.450
Cr0.1590.005
Total mineral elements−0.9160.250−0.9450.219
Total phenolic compounds−0.817−0.285−0.851−0.204
ABTS radical scavenging activity−0.653−0.196
DPPH radical scavenging activity−0.3000.478
Superoxide radical scavenging activity−0.632−0.259
Peroxyl radical inhibition−0.598−0.302
Iron chelation−0.6580.290
Ferric reducing power−0.9070.178−0.9450.110
Table 8. Details of stingless bee honey samples.
Table 8. Details of stingless bee honey samples.
SampleBee SpeciesNectar SourceBotanical OriginHarvest Time
S1Heterotrigona itamaAcacia tree (Acacia mangium)HoneydewAugust 2016
S2Heterotrigona itamaAcacia tree (Acacia mangium)HoneydewNovember 2016
S3 Heterotrigona itamaAcacia tree (Acacia mangium)HoneydewApril 2017
S4Heterotrigona itamaAcacia tree (Acacia mangium)HoneydewJuly 2017
S5Heterotrigona itamaAcacia tree (Acacia mangium)HoneydewSeptember 2017
S6Heterotrigona itamaAcacia tree (Acacia mangium)HoneydewApril 2018
S7Heterotrigona itamaAcacia tree (Acacia mangium)HoneydewJuly 2018
S8Heterotrigona itamaAcacia tree (Acacia mangium)HoneydewSeptember 2018
S9Heterotrigona itamaMultifloralBlossomAugust 2016
S10Heterotrigona itamaMultifloralBlossomNovember 2016
S11Heterotrigona itamaMultifloralBlossomMay 2017
S12Heterotrigona itamaMultifloralBlossomJuly 2017
S13Heterotrigona itamaMultifloralBlossomSeptember 2017
S14Heterotrigona itamaMultifloralBlossomApril 2018
S15Heterotrigona itamaMultifloralBlossomMay 2018
S16Heterotrigona itamaMultifloralBlossomJuly 2018
S17Geniotrigona thoracicaMultifloralBlossomOctober 2016
S18Geniotrigona thoracicaMultifloralBlossomDecember 2016
S19Geniotrigona thoracicaMultifloralBlossomApril 2017
S20Geniotrigona thoracicaMultifloralBlossomJuly 2017
S21Geniotrigona thoracicaMultifloralBlossomMarch 2018
S22Geniotrigona thoracicaMultifloralBlossomJune 2018
S23Geniotrigona thoracicaMultifloralBlossomOctober 2018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ng, W.-J.; Sit, N.-W.; Ooi, P.A.-C.; Ee, K.-Y.; Lim, T.-M. Botanical Origin Differentiation of Malaysian Stingless Bee Honey Produced by Heterotrigona itama and Geniotrigona thoracica Using Chemometrics. Molecules 2021, 26, 7628. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26247628

AMA Style

Ng W-J, Sit N-W, Ooi PA-C, Ee K-Y, Lim T-M. Botanical Origin Differentiation of Malaysian Stingless Bee Honey Produced by Heterotrigona itama and Geniotrigona thoracica Using Chemometrics. Molecules. 2021; 26(24):7628. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26247628

Chicago/Turabian Style

Ng, Wen-Jie, Nam-Weng Sit, Peter Aun-Chuan Ooi, Kah-Yaw Ee, and Tuck-Meng Lim. 2021. "Botanical Origin Differentiation of Malaysian Stingless Bee Honey Produced by Heterotrigona itama and Geniotrigona thoracica Using Chemometrics" Molecules 26, no. 24: 7628. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26247628

Article Metrics

Back to TopTop