Abstract

Stingless bee honey is a good source of antioxidants, which is attributed to the phenolic compounds. The type and concentration of phenolic compounds in honey can be affected by botanical origin. Therefore, in this study, Heterotrigona itama honey from three botanical origins (gelam, acacia, and starfruit) was evaluated for its antioxidant activity and profile of phenolic compounds. Apis mellifera honey was used as a comparison. Antioxidant activity and profile of phenolic compounds in honey were determined using spectrophotometric and chromatographic methods, respectively. The total phenolic content (TPC), total flavonoids content (TFC), free radical scavenging activity (IC50), and ferric reducing antioxidant power (FRAP) of H. itama were ranged between 52.64 and 74.72 mg GAE/100 g honey, 10.70–25.71 mg QE/100 g honey, 11.27–24.09 mg/mL, and 77.88–164.88 µmol FeSO4.7H2O/100 g honey, respectively. The findings showed that the antioxidant activity and phenolic and flavonoid contents in H. itama honey were significantly higher than Apis honey. Benzoic acid and taxifolin were found as the predominant phenolic acid and flavonoid in all samples. However, chrysin was significantly highest in Apis honey than stingless bee honey. This result suggested that chrysin can be used as a chemical marker to distinguish Apis honey from stingless bee honey. Gallic acid and ellagic acid were found as the chemical marker for gelam honey, salicylic acid, benzoic acid, and 4-hydroxybenzoic acid for starfruit honey while ferulic acid for acacia honey.

1. Introduction

Honey is produced by stingless bees and honeybees mainly from the nectar of plants or plant sap [1]. It has been documented that honey is a complex substance and contains over 200 components, which may be either produced during the maturation process of honey, added by bees, or derived from plants [2]. The main components of honey are sugar and water. It also comprises important minor components such as organic acids, amino acids, vitamins, minerals, enzymes, and phenolic compounds. The composition of honey and phenolic compounds primarily depends on the botanical origin [3, 4] and external factors such as harvesting season, environment, storage, and processing method. Phenolic compounds are the most important antioxidant in honey that are responsible for the honey’s therapeutic properties. Thus, honey has been used in traditional and modern treatments to treat human illnesses. In modern treatment, honey has been used to treat diseases associated with oxidative stress such as diabetes mellitus, hypertension, atherosclerosis, cancer, and Alzheimer’s [5]. In addition, phenolic compounds also can affect the organoleptic properties (color, taste, or flavor) of honey [2]. In a recent study, phenolic compounds have been used as a chemical marker to differentiate honey from different botanical and geographical origins [6]. For instance, abscisic acid was reported as a potential floral marker for two Polish unifloral honey [7], and quercetin was suggested as a marker for sunflower honey [8].

Phenolic compounds (phenolic acids and flavonoids) are the secondary components of plants [9] that transfer to honey through nectar, pollen, or propolis by bees [10]. They are responsible for the most antioxidant activity in honey [11]. In recent years, there has been a growing interest in identifying phenolic compounds in stingless bee honey as it has been reported to have higher antioxidant activity and is considered as a potential source of antioxidants compared to Apis mellifera honey, which is good for human health and well-being [12]. Biluca et al. [13], da Silva et al. [14], Sousa et al. [15], and Olivera et al. [16] have identified numerous phenolic compounds in Melipona and Trigona honey. All the studies are from America Latin. However, in Malaysia, there is very limited scientific data on the profile of phenolic compounds in stingless bee honey. Ranneh et al. [17] investigated Malaysian stingless bee honey and Tualang honey. They found gallic, caffeic, syringic, catechin, cinnamic, p-coumaric, apigenin, and 4-hydroxybenzoic as major compounds. Recently, Majid et al. [18] reported six phenolic compounds (chlorogenic acid, p-coumeric acid, epicatechin, rutin, catechin, and protocatechuic acid in unifloral and multifloral honey produced by Heterotrigona itama from Johor. Moreover, there are no available studies on phenolic compounds in stingless bee honey from acacia, gelam, and starfruit. Therefore, more studies need to be conducted to explore and identify more phenolic acids and flavonoids present in stingless bee honey from different botanical origins. In this study, three types of stingless bee honey (acacia, gelam, and starfruit) that are commonly consumed in Malaysia were investigated. The main goal of this present study was to determine the phenolic profiles of stingless bee honey and Apis honey by identifying their phenolic acids and flavonoids, correlate antioxidant activity with phenolic compounds, and verify the botanical origin according to the specific phenolic compounds using multivariate statistical method.

2. Materials and Methods

2.1. Materials

Methanol (MeOH), hydrochloric acid (HCl), formic acid (CH2O2), acetic acid (CH3COOH), sodium acetate (CH3COONa), sodium Carbonate (NaCO3), sodium nitrite (NaNO2), aluminum chloride (AlCl3), sodium hydroxide (NaOH), ferric tripyridyltriazine (Fe3+-TPTZ), iron (III) chloride hexahydrate, 2,2-diphenyl-1-picrylhydrazyl (DPPH) radicals, Folin–Ciocalteu reagent, and amberlite XAD-2 resin were purchased from Fisher Scientific (Pittsburg, PA). Phenolic compound standards (gallic acid, 3, 4-dihydroxybenzoic acid, 4-hydroxybenzoic acid, chlorogenic acid, syringic acid, (-)-epicatechin, (+-)-catechin hydrate, p-coumaric acid, ferulic acid, ellagic acid, benzoic acid, salicylic acid, trans-cinnamic acid, vanillic acid, caffeic acid, taxifolin, myricetin, quercetin hydrate, (+-)-naringenin, and chrysin) were obtained from Sigma Chemical Co. (St. Louis, MO). For HPLC analysis, all chemicals and solvents were of HPLC grade.

2.2. Honey Samples

Natural stingless bee honey produced by H. itama from acacia (Acacia mangium), starfruit (Averrhoa carambola L.), and gelam (Meleleuca cajaputi Powell) were obtained from three different farms that cultivated acacia, gelam, and starfruit trees located in Johor, Terengganu, and Pahang, respectively (Figure 1). The samples were collected during the flowering season between August and November 2016. A. mellifera honey (acacia) was used as a comparison in this study. The botanical origin of the honey samples was determined based on the location of beehives and the availability of the floral source on the farm. The samples were collected using a vacuum pump and filtered to remove foreign materials. They were stored in air-tight containers and kept at 4°C until further analysis. The details of honey samples used are summarized in Table 1.

2.3. Total Phenolic Content

The total phenolic content (TPC) was determined using Folin–Ciocalteu assay according to the method described by Garjanovic et al. [19] with some modifications. A 200 µL of honey solution (0.05 g/mL methanol) was added to 1 mL of Folin–Ciocalteu reagent in a test tube. Then, the mixture was allowed to stand at room temperature. After 6 min, 800 µL of sodium carbonate solution (7.5% Na2CO3) was added to the mixture. Then, the mixture was left to stand for 2 h in a dark room at room temperature. The absorbance of the mixture was measured at 740 nm using a UV-visible spectrophotometer (GENESYS™ 10S, Thermo Fisher Scientific, MA). The calibration curve was prepared using gallic acid at concentrations between 10 and 50 µg/mL (r2 = 0.9999). Results were expressed as mg gallic acid equivalents (GAE) per 100 g. The analysis was performed in triplicate.

2.4. Total Flavonoids Content

The total flavonoid content (TFC) was measured according to the method suggested by Kamboj et al. [20] with minor modifications. One milliliter of honey solution (0.25 g/mL methanol) was mixed with 0.3 mL of 5% sodium nitrite (NaNO2). The mixture was left to stand for 5 min at room temperature. Afterward, 0.3 mL of 10% aluminum chloride (Al3Cl) was added, and the mixture was allowed to stand at room temperature for 6 min. Finally, 2 mL of 1 M of sodium hydroxide (NaOH) was added to the mixture. The absorbance of the mixture was measured at 510 nm using a UV-visible spectrophotometer (GENESYS™ 10S, Thermo Fisher Scientific, MA). Quercetin was used as a standard to plot the calibration curve at the concentrations of 10–50 µg/mL (r2 = 0.98). The results were expressed as mg quercetin equivalent (QE) per 100 g. The analysis was carried out in triplicate.

2.5. Free Radical Scavenging Activity

The scavenging activities of the honey samples were measured by DPPH assay as described by Meda et al. [21] with slight modifications. A 400 mg of honey sample was mixed with different amounts of methanol (3, 6, 9, 12, and 15 mL) separately. Then, 0.75 mL of honey in methanol was mixed with 1.5 mL of DPPH solution. The mixture was left in the dark area for 15 min at room temperature. The absorbance was measured against a blank solution containing water with a UV-visible spectrophotometer (A GENESYSTM 10S, Thermo Fisher Scientific, MA) at 517 nm. The antiradical activity of the sample was calculated using the following formula:where Ac is the absorbance of the control and As is the absorbance of the sample. Each honey sample was analyzed in five dilutions. The antiradical activity was expressed as IC50 (the concentration of honey solution needed to reduce the concentration of DPPH in the solution to 50% of its initial concentration). The concentration of honey sample required to scavenge 50% of DPPH (IC50) was determined from the plotted graph of scavenging activity against the honey dilutions.

2.6. Ferric Reducing Antioxidant Power

The reducing ability of the sample was carried out according to the procedure described by Khalil et al. [22] with slight modifications. The method was based on the reduction of Fe3+-TPTZ to a blue-colored Fe2+ FRAP reagent. It was prepared by mixing 0.3 M of acetate buffer (pH 3.6), 10 mM of TPTZ, and 20 mM FeCl3.6H2O in a ratio of 10:1:1 at 37°C. A 200 µL of the honey sample (0.05 g/mL) in methanol was added to 1.5 mL of FRAP reagent in a test tube. After 10 min of incubation at 37°C in a water bath, the absorbance was measured using the UV-visible spectrophotometer (GENESYS™ 10S, Thermo Fisher Scientific, MA) at 593 nm. The antioxidant potential of the sample was determined from a standard curve using ferrous sulfate heptahydrate FeSO4.7H2O at concentrations between 10 and 100 µg/mL (r2 = 0.98). Triplication of the test was performed, and the results were expressed as the mean average.

2.7. Extraction of Phenolic Compounds Using XAD-2 Resin

The extraction of phenolic compounds was prepared according to the method described by Kassim et al. [23] and Ferreira et al. [24] with some modifications. A 10 g of honey sample was mixed with 100 mL of acidified water (pH 2; 0.02 M HCl) and 150 g of Amberlite XAD-2 resin (Fluka Chemei; pore size = 9 nm; particle size = 0.3–1.2 mm). The mixture was then homogenized for 10 min. The mixture was poured into a column (35 × 2 cm) and washed with 100 mL of acid solution (0.02 M HCl) to remove all sugars and other polar constituents in honey. Subsequently, 300 mL of deionized water was used to rinse the mixture in the column. The phenolic fractions were eluted with 100 mL of methanol and evaporated using a rotary evaporator until dry at 40°C. Methanol was used based on its ability to extract more phenolic compounds as reported by Kassim et al. [23]. The residues were dissolved in 250 µL methanol for HPLC analysis.

2.8. HPLC Analysis

The detection of phenolic compounds was performed using high-performance liquid chromatography equipped with a diode array detector (HPLC-DAD). The column used was a reversed-phase C18 column, ACE (4.6 × 250 mm, particle size 5 µM) from Altmann Analytik GmbH & Co. KG (Munich, Germany). Approximately 10 µL of the phenolic extract was injected into the HPLC, and the separation was performed using 0.25% formic acid and 2% methanol in water (solvent A) and methanol (solvent B). The gradients used were: 10% methanol (B) for 15 min, 40% methanol (B) for 20 min, 45% methanol (B) for 30 min, 60% methanol (B) for 50 min, 80% methanol (B) for 52 min, 90% methanol (B) for 60 min, and followed by isocratic elution with 90% methanol (B) for 65 min. Then, the gradient was changed to 10% methanol for 68 min until 73 min. The flow rate used was 1 mL/min. Phenolic compounds were detected using UV absorption spectra monitored at 290 nm and 340 nm. The phenolic compounds were identified by comparing the chromatographic retention times with the standards. The calibration curve of the standards was used to determine the concentration of the phenolic compounds in the extracts [23].

2.9. Statistical Analysis

All analyses were conducted in triplicate, and the data were expressed as means ± standard deviation. Analysis of variance (ANOVA) and Turkey’s test at the 95% confidence level were performed using MINITAB software (State College, PA). Correlations between antioxidant activity and TPC, TFC, IC50, and FRAP were obtained by Pearson’s correlation analysis. Multivariate analysis was performed using MINITAB software (State College, PA).

3. Results and Discussion

3.1. Antioxidant Properties of Honey
3.1.1. Total Phenolic Content and Total Flavonoids Content

The antioxidant capacity of honey is usually measured by evaluating the contents of phenolic and flavonoids as they are the key compounds responsible for the antioxidant activity of honey [13]. TPC values in the honey samples ranged between 52.64 and 74.62 mg GAE/100 g honey (Table 2). A significant difference () was observed in TPC values between H. itama (gelam and starfruit honey) and A. mellifera honey. Gelam honey showed the significantly highest value of TPC, while acacia honey obtained the lowest value. However, no significant differences () were observed among H. itama honey from different botanical origins. In contrast, Abu Bakar et al. [25] demonstrated the TPC value of H. itama from Malaysia was higher than our results. The TPC values varied between 435.69 and 516.07 mg GAE/100 g honey. Recently, Imtiazah et al. [26] reported 368.11 mg GAE/100 g honey total phenolic content. However, Tufail Ahmad et al. [27] reported 5.86 mg GAE/100 g honey total phenolic content in Malaysian stingless bee honey produced by Geniotrigona thoracica, which was lower than our results. Ranneh et al. [17] also reported lower TPC in stingless bee honey collected from the forest in Kedah and Johor. The values varied between 228.09 and 235.28 mg GAE/kg honey. Ismail et al. [28] reported 33.2–60.2 mg GAE/100 g honey of total phenolic content, which agreed with our study. In comparison to other countries, the results found in this study were comparable to those found by Silva et al. [29] and da Silva et al. [14], where the TPC values were ranged between 1.30 and 66 mg GAE/100 g honey for stingless bee honey from Paraiba and Amazona, Brazil. Furthermore, Silva et al. [29] also reported the higher TPC values in nine Melipona subnitida honey from two semiarid regions, Paraiba, Brazil.

The TFC values of H. itama honey and Apis honey are shown in Table 2. The highest TFC value (25.71 mg QE/100 g honey) was found in starfruit honey, followed by gelam (20.67 mg QE/100 g honey), acacia (10.70 mg QE/100 g honey), and Apis (7.02 mg QE/100 g honey) honey. A significant difference () was observed in all honey samples investigated. In comparison with Malaysian stingless bee honey, Ranneh et al. [17] and Selvaraju et al. [30] reported lower TFC where the values ranged between 97.88 and 101.5 mg CE/kg and 10.18 and 12.68 mg CE/kg, respectively. In a recent study, Imtiazah et al. [26] also demonstrated a lower TFC in Malaysian stingless bee honey (64.25 mg CEQ/kg honey). However, higher TFC was found in multifloral stingless bee honey collected in Sabah and Kelantan with values between 43.2 and 65.9 mg QE/100 g honey [28]. The results from this study were consistent with those reported for stingless bee honey from Alagoas, Brazil, 11.69–49.50 mg QE/100 g honey [31]. Nevertheless, a higher TFC value was found in six species of Melipona honey from Brazil [32] and a lower TFC value in Melipona beecheii honey from Cuba [11].

Overall, stingless bee honey had significantly higher TPC and TFC as compared to Apis honey. A low concentration of phenolic acids and flavonoids may contribute to the low TPC and TFC values of Apis honey as shown in Table 3. In addition to that, the difference can be explained by the different bee species and the way honey is produced by two different bees [33]. While the variation observed between stingless bee honey is due to the different botanical origins of the honey as the botanical origin has a direct association with the phenolic compounds present in honey [34]. Comparison of TPC and TFC data with similar honey from Malaysia and other regions could not be done as there is no related publication available.

3.1.2. The Free Radical Scavenging Activity and Ferric Reducing Antioxidant Power

Antioxidant activities of honey samples were estimated using 2,2-diphenyl-1-picrylhydrazyl (DPPH) and ferric reducing antioxidant power (FRAP) assays. The free radical scavenging activity of honey samples was expressed as IC50, indicating the amount of antioxidant needed to reduce the initial concentration of DPPH solution by 50%. Honey with a low value of IC50 has greater antioxidant activity than honey with a high value of IC50 [35]. Results of the free radical scavenging activity of honey samples are presented in Table 2. The IC50 value of H. itama honey varied between 11.27 and 24.09 mg/mL with the gelam honey having the lowest value. Thus, gelam honey had greater antioxidant activity than starfruit and acacia honey, but no significant difference () was observed among them. Aljadi and Kamaruddin [36] also reported that gelam honey (honeybee honey) had higher antioxidant activity than coconut honey. Kek et al. [37] reported 26.63 mL/g in Malaysian stingless bee honey produced by H. itama. While Ismail et al. [28] reported a range of 10.6–19.7 mg/mL, which was consistent with our results. When compared to the stingless bee honey from Brazil, the IC50 value in this study was lower than the results reported (25.39–51.55 mg/mL). While the IC50 value of H. itama honey (11.27–24.09 mg/mL) was lower than Apis honey (53.65 mg/mL). A similar finding was reported by Duarte et al. [31]. However, there was no significant difference found among the honey samples except between gelam and Apis honey. In contrast, Alvarez-Suarez et al. [12] reported a higher IC50 in M. beecheii honey than Apis honey.

The ferric reducing antioxidant power (FRAP) analysis was performed to measure the ability of phenolic compounds as an antioxidant to reduce ferric ions (Fe3+) to ferrous ions (Fe2+). From the results in Table 2, the FRAP values of all honey samples varied between 72.78 and 164.88 µmol FeSO4.7H2O/100 g honey, where both gelam and starfruit honey displayed significantly () higher FRAP values than those from acacia and Apis honey. Among the analyzed honey samples, starfruit honey showed the highest FRAP value, which indicated that starfruit honey possesses a stronger reducing power than other honey samples. The values of FRAP in this study were higher than those reported by Alvarez-Suarez et al. [12] on Cuban honey (M. beecheii) of 38.54 µmol FeSO4.7H2O/100 g honey but lower than those reported by Chan et al. [38] on Malaysia honey (Trigona spp.) of 3,630.18–7,477.03 µmol FeSO4.7H2O/kg honey and by Nweze et al. [39] on Nigerian honey of 417.36–439.15 µmol FeSO4.7H2O/100 g honey. In recent studies, Majid et al. [23] and Tufail Ahmad et al. [27] reported higher FRAP values of Malaysian stingless bee honey ranging between 283.80 and 1401.80 μM Fe (II)/100 g honey and 27.18 µmol FeSO4.7H2O/g honey, respectively.

The variation in IC50 and FRAP values in the honey samples analyzed could attribute to the different amounts and types of phenolic compounds [40] presented in the honey samples as the samples used in this study belonged to different botanical origins. Therefore, honey from different botanical origins possesses distinct antioxidant activity [14].

3.2. Correlation Analysis

The Pearson’s correlation showed a significant correlation between FRAP and TFC (0.991). This indicates that flavonoids are the greatest contributor to the reducing power of honey. Nevertheless, no significant correlation was observed between TPC and TFC (0.922), between DPPH (IC50) and TPC (−0.813) and TFC (−0.658), and between FRAP and TPC (0.919) and DPPH (IC50) (−0.588). The findings from this study also suggest that the antioxidant activity of honey is not only influenced by the phenolic and flavonoids content but also other antioxidant compounds that partially contribute to the antioxidant activity of honey. As known, honey contains a wide variety of significant antioxidant compounds such as organic acids, amino acids, enzymes, carotenoids, Maillard reaction products, and proteins [34, 41, 42]. Duarte et al. [30] also reported a significant correlation between the FRAP and flavonoids content in Africanized and stingless bee honey from Alagoas, Brazil. Other studies also demonstrated a significant correlation between the FRAP and flavonoids content, suggesting flavonoids as one of the key components responsible for the antioxidant activity of honey [22, 40, 43, 44].

3.3. Limit of Detection and Limit of Quantitation of Phenolic Compounds

The limit of detection (LOD) and limit of quantitation (LOQ) values were calculated based on the 3standard deviation of blank response/slope and 10standard deviation of blank response/slope, respectively. The LOD of phenolic acids and flavonoids were ranged between 0.30 and 148.33 µg/100 g honey and 0.05–0.55 µg/100 g honey, respectively, while the LOQ of phenolic acids and flavonoids were ranged between 1.25 and 494.38 µg/100 g honey and 0.18–1.85 µg/100 g honey, respectively.

3.4. Phenolic Profile

Phenolic compounds identification was performed using HPLC analysis. The chromatograms of the standard of phenolic compounds and honey samples are shown in Figures 2(a) and 2(b), respectively. Table 4 shows the phenolic compounds present in all honey samples. Sixteen compounds were found in starfruit honey, while 15 compounds were identified in acacia and gelam honey. These results indicate that all honey samples possess almost similar phenolic compound profiles; however, different concentrations were quantified in the honey sample from different botanical origins. Ranneh et al. [17] found 13 phenolic compounds (gallic acid, caffeic acid, caffeic acid phenethyl ester, syringic acid, catechin, apigenin, chrysin, cinnamic acid, 2-hydroxybenzoic acid, kaempferol, p-coumaric acid, quercetin-3-O-rutinosid, and 4-hydroxybenzoic acid) in Malaysian stingless bee honey using LC-ESI-MS/MS. On the contrary, Majid et al. [18] identified 6 phenolic compounds (chlorogenic acid, p-coumaric acid, catechin, protocatechuic acid, epicatechin, and rutin) in H. itama honey collected from Johor, Malaysia. Biluca et al. [13] demonstrated the presence of mandelic acid, rosmarinic acid, caffeic acid, aromadendrin, vanillin, isoquercetin, umbelliferone, eriodictyol, sinapaldehyde, syringaldehyde, and carnosol in Brazilian stingless bee honey from different geographical origins. In a different study, Alvarez-Suarez et al. [12] found dihydrocaffeic acid, C-pentosyl-C-hexosyl-apigenin, quercetin deoxyhexosyl hexoside, apigenin trihexoside, kaempferol deoxyhexosyl hexoside, isohamnetin deoxyhexosyl hexoside, isohamnetin, bis-methylated quercetin, apigenin, kaempferol, methyl luteolin, and methyl quercetin in M. beecheii from Cuba. All these compounds were not detected in this present study, which could be explained by the difference in the origin of the honey, the bee species, and the method of analysis used.

Benzoic acid was the most abundant phenolic compound found in all honey samples with values ranging between 738.08 µg/100 g honey–12,626.00 µg/100 g honey. Starfruit honey had the highest benzoic acid (12,626.00 µg/100 g honey), followed by gelam (2,808.40 µg/100 g honey), acacia (2,434.80 µg/100 g honey), and Apis (738.08 µg/100 g honey) honey. Statistical differences were observed between all honey samples. In a recent study, Braghini et al. [45] reported that carnosol was identified as the major phenolic compound in Melipona bicolor honey. In previous studies, fraxin was found abundantly in Ecuadorian stingless bee honey [46] and taxifolin in Melipona (Michmelia) seminigra merrillae honey from Brazil [13]. Regarding the flavonoids content, taxifolin was identified as the major flavonoid (76.50–1,212.70 µg/100 g honey) in all honey samples investigated. However, stingless bee honey (248.69–1,212.70 µg/100 g honey) exhibited significantly higher taxifolin content than Apis honey (76.50 µg/100 g honey). Our result was in accordance with those previously found in Amazon honey from M. (Michmelia) seminigra merrillae [14]. On the other hand, luteolin was found as the predominant flavonoid in Melipona asilvai and Melipona quadrifascita honey; kaempferol in Melipona anthidioides, Melipona scutellaris, and M. subnitida honey; and apigenin in Melipona mandacaia honey from Sergipe state, Brazil [32]. Different types and concentrations of phenolic compounds in honey were reported from different studies, which might be attributed to the different botanical origins, bee species, and geographical locations of honey.

In this present study, ferulic acid, salicylic acid, ellagic acid, benzoic acid, 4-hydroxybenzoic acid, taxifolin, myricetin, quercetin, and naringenin were detected in stingless bee honey. To the best of our knowledge, this is the first report of the presence of these compounds in Malaysian stingless bee honey. The findings revealed that stingless bee honey contains a wide range of phenolic compounds. This finding could have a significant impact on the stingless bee industry sustainability as well as promote Malaysian stingless bee honey globally.

Overall, stingless bee honey exhibited more types and higher concentrations of phenolic compounds compared to Apis honey. This suggests that stingless bee honey is a good source of natural phenolic acids and flavonoids. These findings are consistent with those reported by Alvarez-Suarez et al. [12] and Guerrini et al. [46]. Phenolic acids and flavonoids have been proven to have antioxidant activity, which has a significant correlation with honey’s medicinal properties, and each of them has specific health benefits [47]. Furthermore, stingless bee honey has the potential in treating eye illnesses, wounds, cancer, diabetes mellitus, hypertension, microbial infection, fertility problems, and dysregulated lipid profiles [47]. This indicates that stingless bee has a great potential to be commercialized and consumed by a large population. Acacia, gelam, and starfruit honey can be a good choice as an antioxidant source for humans. Since they have high antioxidant content than Apis honey, we believe that stingless bee honey has a better effect on human health.

3.5. Multivariate Analysis

Multivariate analysis was performed to determine potential phenolic compounds, which can be used as a chemical marker to distinguish stingless bee honey from Apis honey. The principal component analysis (PCA) was used to identify honey based on the concentration of phenolic compounds. To establish the best parameters for classifying honey samples, two principal components with eigenvalues greater than one were extracted. The first principal component (PC1) accounted for 58.6% of the variance, while the second principal component (PC2) accounted for 26.3%. These two components explained 84.9% of the variation in the data. Based on the score plot (Figure 3), all honey samples were scattered away from each other, indicating that all honey samples contained different types and concentrations of phenolic compounds. From the score plot, Apis honey was located far away from the stingless bee honey, proving that Apis honey was significantly different from the stingless bee honey in terms of the phenolic compounds. This finding suggests that phenolic compounds can distinguish Apis honey from stingless bee honey. A biplot was further created to display the relationship between honey samples and phenolic compounds (Figure 4). Then, a loading plot was performed to assess which phenolic compounds contributed significantly to the separation (Figure 5). According to the PC2, chrysin was the only phenolic compound responsible to separate Apis honey from stingless bee honey. As a result, chrysin can be used as a chemical marker to distinguish Apis honey from stingless bee honey.

To study the potential to differentiate acacia, gelam, and starfruit honey, data of Apis honey were removed, and principal component analysis (PCA) was constructed. The first principal component (PC1) accounted for 71.2% of the variance, while the second principal component (PC2) accounted for 27.0%. These two components explained 98.2% of the variation in the data. Based on the score plot (Figure 6), all stingless bee honey were well separated, indicating that all honey samples were significantly different in phenolic compounds content.

The most influential phenolic compounds responsible for clustering to be identified as a potential marker for all honey samples were determined using biplot (Figure 7(a)) and loading plot (Figure 7(b)). Based on the loading plot, gallic acid and ellagic acid have a strong impact on the separation of gelam honey. Salicylic acid, benzoic acid, and 4-hydroxybenzoic acid are responsible for the separation of starfruit honey, while ferulic acid is the only one that separated acacia from other honey samples.

4. Conclusion

This study discovered that botanical origin significantly influences antioxidant activity and the composition of phenolic compounds in honey. Flavonoids were shown to have a significant relationship with the antioxidant activity (FRAP) of honey. In addition, HPLC analysis indicated and validated that stingless bee honey is rich in phenolic components as compared to Apis honey, showing that stingless bee honey is a valuable source of natural phenolic acids and flavonoids. Furthermore, the findings suggest that the botanical origin has a significant impact on the content of phenolic compounds in honey. A multivariate data analysis was used to find a potential marker that might be used to distinguish honey samples based on their phenolic compounds. The results showed that chrysin successfully distinguished the stingless bee honey from Apis honey. Gelam honey was discriminated by gallic acid and ellagic acid, starfruit honey by salicylic acid, benzoic acid, and 4-hydroxybenzoic acid, while acacia honey by ferulic acid. However, the number of samples used in this study was limited. Thus, a further investigation with a large number of samples is required in order to confirm the present results.

Data Availability

Data used to support the findings of this study are included in the article.

Conflicts of Interest

The authors declare that there are no potential conflicts of interest.

Acknowledgments

The authors like to thank the Ministry of Higher Education (MOHE) for the research support (FRGS/1/2020/STG04/UPM/01/1) and the HICoE awarded to the Institute of Tropical Agriculture and Food Security (ITAFoS), Universiti Putra Malaysia. The authors would also like to thank ITAFoS and the Faculty of Food Science and Technology for the research facilities rendered.