# Characterization of Various Honey Samples from Different Regions of Morocco Using Physicochemical Parameters, Minerals Content, Antioxidant Properties, and Honey-Specific Protein Pattern

Authors: Badiaa Lyoussi; Meryem Bakour; Redouan El-Haskoury; Hamada Imtara; Christophe Hano; Katarína Bíliková
Journal: Journal of Food Quality
 (2022)

---

## Abstract

Honey is a bee product relatively expensive; therefore, it has been a target of adulteration by many sweeteners. In this work, we evaluated the good quality, authenticity, and content in bioactive molecules of twenty-two Moroccan honey from different botanical origins and geographical areas. For that, the following analyses were determined: the content in total protein and especially the major royal jelly protein (apalbumin 1), the analysis of total acidity, free acidity, lactonic acidity, pH, ash, Pfund, electrical conductivity, and moisture. In addition, the content of sodium, potassium, calcium, and magnesium, the dosage of polyphenols, flavones, and flavonols, and the antioxidant activities were assessed. All analyzed samples had good antioxidant activities and present a source of antioxidant compounds, the predominant mineral in all honey samples was potassium, and the physicochemical parameters are in line with the standards’ recommended limits. The content of honey samples in total protein and apalbumin 1 ranged between 212μg/g and 4121.2 μg/g and between 27.4 μg/g and 790.82 μg/g, respectively. Overall, the detection of apalbumin 1 in all honey samples and the results of physicochemical parameters, minerals, bioactive compounds, and antioxidant activities confirm the authenticity and no adulteration of Moroccan honey.


---

## Body

## 1. Introduction

Honey is a sweet product produced by bees from the nectar of plants. It is nutritious and has traditionally been consumed by humans since the oldest times [1]. Honey is a complex mixture that comprises carbohydrates (60–85%) mainly glucose and fructose, water (12–23%), and other minor constituents such as proteins, enzymes, free amino acids, lipids, vitamins, phenolic acids, flavonoids, and mineral salts [2]. The biochemical composition of honey mostly depends on its floral source, the honey bee species, weather conditions, and geographical origin [3].There is a large volume of published studies describing the role of honey as functional food, it has been reported that honey has several pharmacological effects such as antioxidant, anticancer, antimicrobial, and inflammatory effects [4, 5]. It is widely used in wound healing and can counteract inflammation [6].Honey is relatively expensive; thus, it has been a target of adulteration by many adulterants such as corn syrup, sugar, and cane [7]. To verify the authenticity of honey, many tools are recommended such as physicochemical parameters, sensory analysis, microscopical examination, and the analysis of chemical composition [1, 8]. In addition to that, apalbumin 1 is among the main protein that exists in royal jelly and also in honey; it has been suggested to use it as a marker of authenticity and quality of honey since this protein is specific to bees and cannot be replaced by other ingredients [9]. Therefore, the objective of this work was to determine the authenticity and quality of honey samples of different botanical and geographical origins in Morocco by the analysis of physicochemical parameters (pH, free acidity, lactone acidity, total acidity, moisture, electrical conductivity, ash, and Pfund), minerals content (sodium, potassium, calcium, magnesium), antioxidant content and activities (polyphenols, flavones/flavonols, DPPH, ABTS, RP), and the analysis of total protein and more specifically apalbumin 1.


## 2. Material and Methods

### 2.1. Honey Samples

Twenty-two honey samples were obtained from beekeepers, who installed their hives in seven different regions of Morocco. Three samples were multifloral, and nineteen were monofloral (the pollen grains of the predominant plant are higher than 45%). The honey samples S1, S4, S8, S11, S12, S18, S20, and S21 were obtained from the Fez-Meknes region. The honey samples S2, S6, S13, S14, S17, and S22 were obtained from the Rabat-Salé-Kénitra region. The honey sample S3 was obtained from the Souss-Massa region. The honey samples S5 and S7 were presented from the Oriental region. The sample S9 was presented from the Drâa-Tafilalet region. The honey samples S10, S15, and S16 were obtained from the Beni Mellal-Khenefra region, and honey sample S19 was obtained from the Tangier-Tetouan-Al Hoceïma. The honey samples were produced by three different bee species:Apis mellifera intermissa, Apis mellifera sahariensis, and Apis mellifera major (Table 1 and Figure 1).Table 1
Information on the harvested year, the botanical, bee species, and geographical origin of the honey samples.
SampleLocal nameBotanical origin (Predominant pollen grains)Bee speciesLocationRegionYear of harvestS1AzirSalvia rosmarinus Spenn.Apis mellifera intermissaTanditFez-Meknes2015S2BouchnikhaAmmi visnaga L.Apis mellifera intermissaHedkourtRabat-Salé-Kénitra2015S3DaghmoussEuphorbia resinifera Berg.Apis mellifera sahariensisTiznitSouss-Massa2015S4SadraZiziphus lotus L.Apis mellifera intermissaMarmouchaFez-Meknes2015S5LatchinCitrus sinensis L.Apis mellifera intermissaNadorOriental2015S6LatchinCitrus sinensis L.Apis mellifera intermissaKhnichatRabat-Salé-Kénitra2015S7MultifloralMultifloralApis mellifera intermissaNadorOriental2015S8MultifloralMultifloralApis mellifera intermissaElmersFez-Meknes2015S9MultifloralMultifloralApis mellifera sahariensisSkouraDrâa-Tafilalet2015S10KharobCeratonia siliqua L.Apis mellifera intermissaKhenifraBéniMellal-Khénifra2015S11KharobCeratonia siliqua L.Apis mellifera intermissaTaounateFez-Meknes2015S12KharobCeratonia siliqua L.Apis mellifera intermissaTaounateFez-Meknes2015S13KhzamaLavandula angustifolia L.Apis mellifera intermissaOulmèsRabat-Salé-Kénitra2017S14HamdCitrus limon (L.) Burm. F.Apis mellifera intermissaKhnichatRabat-Salé-Kénitra2017S15SadraZiziphus lotus L.Apis mellifera intermissaKhenifraBéniMellal-Khénifra2017S16ChoukSilybum marianum (L.) Gaertn.Apis mellifera intermissaKhenifraBéniMellal-Khénifra2017S17BouchnikhaAmmi visnaga L.Apis mellifera intermissaSidi KacemRabat-Salé-Kénitra2017S18KebbarCapparis spinosa L.Apis mellifera intermissaMy yacoubFez-Meknes2017S19BakhenouArbutus unedo L.Apis mellifera majorTetouanTanger-Tétouan-AlHoceïma2017S20ZandazBupleurum spinosum GouanApis mellifera intermissaBoulemaneFez-Meknes2017S21Z’îtraThymus vulgaris L.Apis mellifera intermissaTimhditFez-Meknes2017S22ZaatarOriganum vulgare L.Apis mellifera intermissaErrachidiaRabat-Salé-​Kénitra2017Figure 1
Map of Morocco showing honey sampling regions.


### 2.2. Melissopalynological Analysis

The botanical origin of honey samples was determined using the method described by Louveaux et al. [10]. A minimum of 1000 pollen grains were counted for each honey sample under a microscope. If the percentage of any type of pollen grains found in honey exceeds 45% of the total pollen grains content, it is classified as the predominant and the honey is classified as monofloral.


### 2.3. Physicochemical Analysis

Total acidity, free acidity, lactonic acidity, pH, ash, electrical conductivity, and moisture, were analyzed using the methods recommended by the International Honey Commission [11]. Pfund was determined as described previously by Laaroussi et al. [12].


### 2.4. Minerals Content

The analysis of minerals elements (Na, K, Ca, and Mg) of honey samples was carried out using ICP- AES after the calcination method as described by Laaroussi et al. [12]. Briefly, the honey ashes were mixed with 5 ml of nitric acid 0.1 M and stirred on a heating plate until the total evaporation of nitric acid. Then, 10 ml of nitric acid was added and the mixture was made up to 25 ml with ultrapure water. All samples were analyzed in triplicate.


### 2.5. Polyphenols Content

The polyphenols content was assessed using the Folin-Ciocalteu method [13]. Briefly, 100 μL of aqueous extract of honey was mixed with 500 μl of Folin-Ciocalteau reagent solution (10 g of sodium tungstate and 2.5 g of sodium molybdate (2.5 g) were dissolved in 70 ml of distilled water; then, 5 ml of phosphoric acid (85%) and 10 ml of concentrated hydrochloric acid were added. After 10 hours, 15 g of lithium sulfate and 5 ml of distilled water were added and brought to 100 ml with distilled water) for 6 min, and then 400 μl of sodium carbonate (75 g/l) was added to the mixture. The absorbance was measured at 760 nm after 15 min of incubation. Gallic acid calibration was used as a standard for calibration. The results were expressed as milligrams of Gallic acid equivalents per gram (mg GAE/g).


### 2.6. Flavones and Flavonols Content

The content of flavones and flavonols was quantified by a colorimetric assay described previously by Bakour et al. [14]. Briefly, 500 μl of honey sample or standard was added to 500 μl of 20% AlCl3. After 1 h at room temperature, the absorbance was measured at 420 nm. A quercetin calibration curve was prepared, and total content was expressed as mg of quercetin equivalents per 100 g of honey (mg QE/100 g).


### 2.7. Radical Scavenging Activity (DPPH Assay)

The radical scavenging activity of the honey solution against DPPH free radical was measured using the method described by Bakour et al. [15]. Briefly, 100 μl of the aqueous honey extract was mixed with 900 μl of a 100 μM solution of DPPH radical prepared in ethanol (96%). The absorbance of the solution was measured at 540 nm after 30 min of incubation in the dark. Several concentrations of samples were made, and the IC50 (concentration of sample able to scavenge 50% of DPPH free radical) was determined graphically using the curve plotted by the percentage of DPPH inhibition as a function of the sample concentration:(1)%inhibition=Abs control−Abs sampleAbs control×100.


### 2.8. Azino-Bis (3-Ethylbenzothiazoline-6-Sulphonic Acid (ABTS)

Azino-bis (3-ethylbenzothiazoline-6-sulphonic acid (ABTS) free radical scavenging activity was analyzed using the method described by Bakour et al. [15]. 75 μl of aqueous extract of honey or standard control (BHT) was mixed with 825 μl of ABTS solution, and the absorbance of the mixture was measured after 6 min at 734 nm. The tests were carried out in triplicate, and the IC50 (concentration of sample able to scavenge 50% of ABTS free radical) was determined graphically using the curve plotted by the percentage of ABTS inhibition as a function of the sample concentration (equation (1)).


### 2.9. Ferric Reducing Power

The reducing power of the aqueous honey extract was determined using the method described by Bakour et al. [15]. 50 μl of aqueous honey extract (50% W/V) was mixed with 200 μl of 0.2 M sodium phosphate buffer (pH 6.6) and 200 μl of 1% potassium ferricyanide. The mixture was incubated at 50°C for 20 min. Then, 200 μl of 10% trichloroacetic acid (TCA) was added, and the mixture was centrifuged at 3000 rpm for 10 min. 500 μl of the above solution from each reaction was diluted with 500 μl of distilled water, and 100 μl of 0.1% ferric chloride (FeCl3) was added. The absorbance was measured at 700 nm, and ascorbic acid was used as a reference standard. The results were represented in EC50 values, corresponding to the concentration providing 50% of the antioxidant activity or 0.5 of absorbance in the reducing power assay measured at 700 nm.


### 2.10. Determination of Protein Concentration

The total protein content of the honey samples was determined by microplate assay according to Bradford [16]. To 100 μl of the sample or its dilution in physiological solution was added 100 μl of QuickStart Bradford reagent (BioRad, Laboratories, Inc., USA). The absorbance was measured at 595 nm. Bovine serum albumin (BSA, Sigma, USA) was used as standard. Each sample was analyzed at three dilutions and each dilution in three parallel analyses.


### 2.11. Determination of Apalbumin 1 in Honey by ELISA

Honey samples were analyzed using enzyme-linked immunosorbent assay (ELISA) for apalbumin 1 quantification as described previously in detail [17]. The 96 well/flat-bottom microtiter plates (Brand, Germany) were coated with antigen-diluted honey samples at dilution of 0.05% and/or 0.001% in Milli-Q water and/or standard solution of apalbumin 1 and incubated overnight at 4°C. After washing with TBS buffer (100 mmol/L Tris and 150 mmol/L NaCl, pH 7.5), the plates were incubated with polyclonal rabbit anti-apa1 antibody in milk buffer (2% nonfat milk in TBS) and then with peroxidase-conjugated anti-rabbit IgG in milk buffer for 1 h. Detection was done by adding 3% ABTS (2,2′-azino-bis-(3-benzthiazoline-6-sulfonic acid, Southern Biotech, USA), in 50 mmol/L citrate buffer pH 4.3; supplemented by hydrogen peroxide. The absorbance at 405 nm was read in a Microplate Spectrophotometer Power Wave TM XS (BioTek Instruments, INC, Winooski, Vermont, USA).


### 2.12. Statistical Analysis

All data are presented as mean ± SD (standard deviation). Graphpad prism (version 5.0; GraphPad Software, Inc., San Diego, USA) was used to compare honey samples using a one-way analysis of variance (ANOVA) followed by Tukey’s test, andp<0.05 was considered significant. Correlations between the parameters studied were achieved by the Pearson correlation coefficient (r). The principal component analysis (PCA) was accomplished using Past: paleontological statistics software package for education and data analysis, version 3.20.


## 2.1. Honey Samples

Twenty-two honey samples were obtained from beekeepers, who installed their hives in seven different regions of Morocco. Three samples were multifloral, and nineteen were monofloral (the pollen grains of the predominant plant are higher than 45%). The honey samples S1, S4, S8, S11, S12, S18, S20, and S21 were obtained from the Fez-Meknes region. The honey samples S2, S6, S13, S14, S17, and S22 were obtained from the Rabat-Salé-Kénitra region. The honey sample S3 was obtained from the Souss-Massa region. The honey samples S5 and S7 were presented from the Oriental region. The sample S9 was presented from the Drâa-Tafilalet region. The honey samples S10, S15, and S16 were obtained from the Beni Mellal-Khenefra region, and honey sample S19 was obtained from the Tangier-Tetouan-Al Hoceïma. The honey samples were produced by three different bee species:Apis mellifera intermissa, Apis mellifera sahariensis, and Apis mellifera major (Table 1 and Figure 1).Table 1
Information on the harvested year, the botanical, bee species, and geographical origin of the honey samples.
SampleLocal nameBotanical origin (Predominant pollen grains)Bee speciesLocationRegionYear of harvestS1AzirSalvia rosmarinus Spenn.Apis mellifera intermissaTanditFez-Meknes2015S2BouchnikhaAmmi visnaga L.Apis mellifera intermissaHedkourtRabat-Salé-Kénitra2015S3DaghmoussEuphorbia resinifera Berg.Apis mellifera sahariensisTiznitSouss-Massa2015S4SadraZiziphus lotus L.Apis mellifera intermissaMarmouchaFez-Meknes2015S5LatchinCitrus sinensis L.Apis mellifera intermissaNadorOriental2015S6LatchinCitrus sinensis L.Apis mellifera intermissaKhnichatRabat-Salé-Kénitra2015S7MultifloralMultifloralApis mellifera intermissaNadorOriental2015S8MultifloralMultifloralApis mellifera intermissaElmersFez-Meknes2015S9MultifloralMultifloralApis mellifera sahariensisSkouraDrâa-Tafilalet2015S10KharobCeratonia siliqua L.Apis mellifera intermissaKhenifraBéniMellal-Khénifra2015S11KharobCeratonia siliqua L.Apis mellifera intermissaTaounateFez-Meknes2015S12KharobCeratonia siliqua L.Apis mellifera intermissaTaounateFez-Meknes2015S13KhzamaLavandula angustifolia L.Apis mellifera intermissaOulmèsRabat-Salé-Kénitra2017S14HamdCitrus limon (L.) Burm. F.Apis mellifera intermissaKhnichatRabat-Salé-Kénitra2017S15SadraZiziphus lotus L.Apis mellifera intermissaKhenifraBéniMellal-Khénifra2017S16ChoukSilybum marianum (L.) Gaertn.Apis mellifera intermissaKhenifraBéniMellal-Khénifra2017S17BouchnikhaAmmi visnaga L.Apis mellifera intermissaSidi KacemRabat-Salé-Kénitra2017S18KebbarCapparis spinosa L.Apis mellifera intermissaMy yacoubFez-Meknes2017S19BakhenouArbutus unedo L.Apis mellifera majorTetouanTanger-Tétouan-AlHoceïma2017S20ZandazBupleurum spinosum GouanApis mellifera intermissaBoulemaneFez-Meknes2017S21Z’îtraThymus vulgaris L.Apis mellifera intermissaTimhditFez-Meknes2017S22ZaatarOriganum vulgare L.Apis mellifera intermissaErrachidiaRabat-Salé-​Kénitra2017Figure 1
Map of Morocco showing honey sampling regions.


## 2.2. Melissopalynological Analysis

The botanical origin of honey samples was determined using the method described by Louveaux et al. [10]. A minimum of 1000 pollen grains were counted for each honey sample under a microscope. If the percentage of any type of pollen grains found in honey exceeds 45% of the total pollen grains content, it is classified as the predominant and the honey is classified as monofloral.


## 2.3. Physicochemical Analysis

Total acidity, free acidity, lactonic acidity, pH, ash, electrical conductivity, and moisture, were analyzed using the methods recommended by the International Honey Commission [11]. Pfund was determined as described previously by Laaroussi et al. [12].


## 2.4. Minerals Content

The analysis of minerals elements (Na, K, Ca, and Mg) of honey samples was carried out using ICP- AES after the calcination method as described by Laaroussi et al. [12]. Briefly, the honey ashes were mixed with 5 ml of nitric acid 0.1 M and stirred on a heating plate until the total evaporation of nitric acid. Then, 10 ml of nitric acid was added and the mixture was made up to 25 ml with ultrapure water. All samples were analyzed in triplicate.


## 2.5. Polyphenols Content

The polyphenols content was assessed using the Folin-Ciocalteu method [13]. Briefly, 100 μL of aqueous extract of honey was mixed with 500 μl of Folin-Ciocalteau reagent solution (10 g of sodium tungstate and 2.5 g of sodium molybdate (2.5 g) were dissolved in 70 ml of distilled water; then, 5 ml of phosphoric acid (85%) and 10 ml of concentrated hydrochloric acid were added. After 10 hours, 15 g of lithium sulfate and 5 ml of distilled water were added and brought to 100 ml with distilled water) for 6 min, and then 400 μl of sodium carbonate (75 g/l) was added to the mixture. The absorbance was measured at 760 nm after 15 min of incubation. Gallic acid calibration was used as a standard for calibration. The results were expressed as milligrams of Gallic acid equivalents per gram (mg GAE/g).


## 2.6. Flavones and Flavonols Content

The content of flavones and flavonols was quantified by a colorimetric assay described previously by Bakour et al. [14]. Briefly, 500 μl of honey sample or standard was added to 500 μl of 20% AlCl3. After 1 h at room temperature, the absorbance was measured at 420 nm. A quercetin calibration curve was prepared, and total content was expressed as mg of quercetin equivalents per 100 g of honey (mg QE/100 g).


## 2.7. Radical Scavenging Activity (DPPH Assay)

The radical scavenging activity of the honey solution against DPPH free radical was measured using the method described by Bakour et al. [15]. Briefly, 100 μl of the aqueous honey extract was mixed with 900 μl of a 100 μM solution of DPPH radical prepared in ethanol (96%). The absorbance of the solution was measured at 540 nm after 30 min of incubation in the dark. Several concentrations of samples were made, and the IC50 (concentration of sample able to scavenge 50% of DPPH free radical) was determined graphically using the curve plotted by the percentage of DPPH inhibition as a function of the sample concentration:(1)%inhibition=Abs control−Abs sampleAbs control×100.


## 2.8. Azino-Bis (3-Ethylbenzothiazoline-6-Sulphonic Acid (ABTS)

Azino-bis (3-ethylbenzothiazoline-6-sulphonic acid (ABTS) free radical scavenging activity was analyzed using the method described by Bakour et al. [15]. 75 μl of aqueous extract of honey or standard control (BHT) was mixed with 825 μl of ABTS solution, and the absorbance of the mixture was measured after 6 min at 734 nm. The tests were carried out in triplicate, and the IC50 (concentration of sample able to scavenge 50% of ABTS free radical) was determined graphically using the curve plotted by the percentage of ABTS inhibition as a function of the sample concentration (equation (1)).


## 2.9. Ferric Reducing Power

The reducing power of the aqueous honey extract was determined using the method described by Bakour et al. [15]. 50 μl of aqueous honey extract (50% W/V) was mixed with 200 μl of 0.2 M sodium phosphate buffer (pH 6.6) and 200 μl of 1% potassium ferricyanide. The mixture was incubated at 50°C for 20 min. Then, 200 μl of 10% trichloroacetic acid (TCA) was added, and the mixture was centrifuged at 3000 rpm for 10 min. 500 μl of the above solution from each reaction was diluted with 500 μl of distilled water, and 100 μl of 0.1% ferric chloride (FeCl3) was added. The absorbance was measured at 700 nm, and ascorbic acid was used as a reference standard. The results were represented in EC50 values, corresponding to the concentration providing 50% of the antioxidant activity or 0.5 of absorbance in the reducing power assay measured at 700 nm.


## 2.10. Determination of Protein Concentration

The total protein content of the honey samples was determined by microplate assay according to Bradford [16]. To 100 μl of the sample or its dilution in physiological solution was added 100 μl of QuickStart Bradford reagent (BioRad, Laboratories, Inc., USA). The absorbance was measured at 595 nm. Bovine serum albumin (BSA, Sigma, USA) was used as standard. Each sample was analyzed at three dilutions and each dilution in three parallel analyses.


## 2.11. Determination of Apalbumin 1 in Honey by ELISA

Honey samples were analyzed using enzyme-linked immunosorbent assay (ELISA) for apalbumin 1 quantification as described previously in detail [17]. The 96 well/flat-bottom microtiter plates (Brand, Germany) were coated with antigen-diluted honey samples at dilution of 0.05% and/or 0.001% in Milli-Q water and/or standard solution of apalbumin 1 and incubated overnight at 4°C. After washing with TBS buffer (100 mmol/L Tris and 150 mmol/L NaCl, pH 7.5), the plates were incubated with polyclonal rabbit anti-apa1 antibody in milk buffer (2% nonfat milk in TBS) and then with peroxidase-conjugated anti-rabbit IgG in milk buffer for 1 h. Detection was done by adding 3% ABTS (2,2′-azino-bis-(3-benzthiazoline-6-sulfonic acid, Southern Biotech, USA), in 50 mmol/L citrate buffer pH 4.3; supplemented by hydrogen peroxide. The absorbance at 405 nm was read in a Microplate Spectrophotometer Power Wave TM XS (BioTek Instruments, INC, Winooski, Vermont, USA).


## 2.12. Statistical Analysis

All data are presented as mean ± SD (standard deviation). Graphpad prism (version 5.0; GraphPad Software, Inc., San Diego, USA) was used to compare honey samples using a one-way analysis of variance (ANOVA) followed by Tukey’s test, andp<0.05 was considered significant. Correlations between the parameters studied were achieved by the Pearson correlation coefficient (r). The principal component analysis (PCA) was accomplished using Past: paleontological statistics software package for education and data analysis, version 3.20.


## 3. Results and Discussion

### 3.1. Physicochemical Parameters

Twenty-two honey samples of different botanical and geographical origins of Morocco were analyzed for physicochemical parameters. The results presented in Table2 showed that the pH values ranged between 3.59 ± 0.10 in Euphorbia resinifera honey from Tiznit and 4.30 ± 0.20 in Ziziphus lotus honey from Marmoucha; these values are within the range recommended by Council Directive 2001/110/EC [18]. Free acidity ranged between 11.45 ± 0.12 mEq/kg in Salvia rosmarinus honey and 30.81 ± 1.20 mEq/kg in Thymus vulgaris from Timhdit, lactone acidity ranged between 5.09 ± 0.96 mEq/kg in Ziziphus lotus honey and 17.50 ± 0.22 mEq/kg in Bupleurum spinosum honey, and total acidity ranged between 18.30 ± 0.30 mEq/kg in Citrus sinensis honey and 42.28 ± 0.30 mEq/kg in Ceratonia siliqua honey. These values are in line with those recommended by Codex Alimentarius Commission [19]. For moisture values, all honey samples are below the maximum limit (20%) except the Arbutus unedo honey sample from Tetouan which was slightly higher (20.90 ± 0.15%). The recommended value of moisture in honey indicates its maturity and reflects the desirable density [20]. For electrical conductivity, all honey samples are bellowing the maximum limit allowed (800 μS/cm) except Ceratonia siliqua honey from Taounate (948.33 ± 1.63 μS/cm) and Arbutus unedo honey from Tetouan (913.51 ± 6.21 μS/cm). Similarly, the ash contents in all honey samples are below the limits (0.6%) except for the multifloral honey sample from Skoura (0.92 ± 0.02%). The ash content reflects the mineral composition in honey. It is influenced by soil and botanical origins [12]. On the other hand, Pfund is a parameter that reflects honey color, the analysis of Pfund showed values ranging between 19.67 ± 0.64 mm in Citrus limon honey and 128.74 ± 1.20 mm in Ceratonia siliqua honey. According to the Pfund scale, the color of our honey samples ranged between white color and dark amber [21].Table 2
Physicochemical parameters and minerals content in honey samples.
SamplepHFree acidity (mEq/kg)Lactone acidity (mEq/kg)Total acidity (mEq/kg)Moisture (%)Electrical conductivity (mS/cm)Ash%Pfund (mm)Na (mm/kg)K (mm/kg)Ca (mm/kg)Mg (mm/kg)S14.05 ± 0.03b11.45 ± 0.12hi9.51 ± 0.50cd20.96 ± 0.15j17.20 ± 0.22bc443.00 ± 20.00h0.28 ± 0.01g83.06 ± 0.01j41.64 ± 0.54o270.57 ± 0.85s28.49 ± 0.21o19.76 ± 1.25nS23.71 ± 0.02b21.00 ± 1.00d11.10 ± 0.10cd32.10 ± 0.35e18.60 ± 0.41b472.33 ± 3.22g0.31 ± 0.01ef64.21 ± 0.10m128.16 ± 0.4f286.78 ± 0.42r25.98 ± 1.20o17.44 ± 0.56oS33.59 ± 0.10b20.51 ± 0.50d8.10 ± 0.36f28.61 ± 0.52g17.50 ± 0.50bc444.00 ± 20.00h0.33 ± 0.02ef75.01 ± 0.12l117.89 ± 1.21h667.46 ± 0.78h77.27 ± 1.35i42.36 ± 1.25eS44.30 ± 0.20a14.41 ± 0.90g7.51 ± 0.12f21.92 ± 0.45j17.70 ± 0.23bc620.83 ± 1.53c0.42 ± 0.05d63.43 ± 0.11m154.35 ± 0.54d1014.28 ± 2.36c67.57 ± 0.78k13.55 ± 0.56pS54.07 ± 0.04b12.00 ± 0.15hi6.30 ± 0.22f18.30 ± 0.30k19.80 ± 0.12a242.00 ± 30.00m0.17 ± 0.01gh25.58 ± 0.23q114.74 ± 0.41i424.24 ± 1.25o79.35 ± 1.20h25.76 ± 0.87jS64.20 ± 0.10a13.51 ± 0.60h6.00 ± 0.12f19.51 ± 0.40k19.50 ± 0.52a223.00 ± 40.00n0.20 ± 0.01g33.13 ± 0.13p71.48 ± 0.52l622.95 ± 2.32j86.36 ± 0.36g28.15 ± 0.47hS73.94 ± 0.20a14.45 ± 0.12g7.41 ± 0.52f21.86 ± 0.40j17.10 ± 0.36d577.33 ± 4.22d0.35 ± 0.01e103.77 ± 0.62e114.63 ± 0.21i947.36 ± 3.25e93.21 ± 0.98f35.13 ± 0.36fS83.94 ± 0.40a16.05 ± 0.80g10.11 ± 0.63cd26.16 ± 0.30h18.80 ± 0.20b510.67 ± 1.23e0.55 ± 0.02c78.74 ± 0.14k78.16 ± 1.21k847.21 ± 1.25f167.57 ± 1.25d218.72 ± 1.54aS93.85 ± 0.00b21.51 ± 0.16d8.20 ± 0.10f29.71 ± 0.62f18.70 ± 0.12b762.33 ± 2.33b0.92 ± 0.02a76.19 ± 0.20l153.04 ± 1.54d987.89 ± 1.25d22.12 ± 0.21p34.35 ± 0.54gS103.75 ± 0.10b11.45 ± 0.90hi9.81 ± 0.41cd21.26 ± 0.16j17.10 ± 0.32bc484.33 ± 4.53f0.31 ± 0.01ef102.92 ± 0.35e131.15 ± 0.65e1016.13 ± 4.23c26.00 ± 1.23o28.07 ± 0.88hS113.79 ± 0.08b26.10 ± 0.60bc12.15 ± 0.10c38.25 ± 0.55bc18.80 ± 0.22b948.33 ± 1.63a0.74 ± 0.01b108.91 ± 0.35c174.24 ± 2.87c806.56 ± 1.23g252.43 ± 0.508b60.97 ± 0.58dS124.00 ± 0.01a28.21 ± 0.23b14.00 ± 0.36b42.28 ± 0.30a18.50 ± 0.52b700.00 ± 40.00b0.52 ± 0.04c128.74 ± 1.20a233.92 ± 1.32b1017.98 ± 3.24c80.3 ± 0.69i62.52 ± 0.41cS134.18 ± 0.12a18.15 ± 0.90f8.03 ± 0.09f26.18 ± 0.22hi17.30 ± 0.23bc345.44 ± 3.12k0.31 ± 0.05e40.83 ± 0.52o104.77 ± 0.65j588.92 ± 1.25l59.44 ± 1.25l40.17 ± 0.78eS143.80 ± 0.20b20.56 ± 0.12d6.70 ± 0.23f27.26 ± 0.30h18.74 ± 0.62b231.83 ± 2.21n0.13 ± 0.07gh19.67 ± 0.64r41.66 ± 0.36o331.89 ± 2.32p37.50 ± 1.80n19.44 ± 0.44lS154.08 ± 0.02a21.47 ± 0.40d5.09 ± 0.96g26.56 ± 0.52h19.28 ± 0.20a521.60 ± 2.36e0.42 ± 0.09e56.20 ± 0.12n121.30 ± 0.21g616.38 ± 1.23k77.90 ± 0.36i24.02 ± 0.21jS163.82 ± 0.01b18.05 ± 0.22f7.15 ± 0.52f25.20 ± 0.41hi19.15 ± 0.15a329.70 ± 5.21l0.20 ± 0.01g87.34 ± 0.69i43.30 ± 0.22o494.05 ± 0.23n42.05 ± 0.25m20.24 ± 0.14lS173.99 ± 0.01a28.43 ± 0.12b6.85 ± 0.23f35.28 ± 0.30d18.15 ± 0.52b489.99 ± 5.21f0.21 ± 0.06g97.24 ± 0.54f112.75 ± 0.25i303.67 ± 1.24q43.97 ± 0.87m17.56 ± 0.25mS183.94 ± 0.02a20.04 ± 0.25de11.37 ± 0.63c31.41 ± 0.40e17.75 ± 0.63bc428.89 ± 6.12i0.29 ± 0.04g119.65 ± 0.41b78.42 ± 1.22k585.99 ± 1.23l172.42 ± 1.25c36.80 ± 0.87fS194.14 ± 0.05a27.40 ± 0.63b7.70 ± 0.52f35.10 ± 0.23d20.90 ± 0.15a913.51 ± 6.21b0.49 ± 0.04d103.75 ± 0.25e264.56 ± 0.45a1299.56 ± 3.25a283.13 ± 2.36a77.23 ± 0.59bS203.76 ± 0.06b22.05 ± 0.41d17.50 ± 0.22a39.55 ± 0.60b19.67 ± 0.34a466.67 ± 3.21g0.20 ± 0.05g89.13 ± 0.74h64.56 ± 0.98m521.26 ± 0.98m121.47 ± 1.54e26.42 ± 0.54iS213.82 ± 0.02b30.81 ± 1.20a8.45 ± 0.40e39.26 ± 0.52b17.32 ± 0.42bc372.04 ± 2.35j0.48 ± 0.05d94.18 ± 0.98g103.29 ± 1.02j1105.58 ± 2.36b76.26 ± 0.65i30.10 ± 0.86hS223.94 ± 0.04a18.70 ± 0.17f7.25 ± 0.54f25.95 ± 0.45hi19.00 ± 0.56a582.35 ± 4.21d0.32 ± 0.01e106.63 ± 1.21d51.60 ± 0.45n632.20 ± 3.20i71.85 ± 0.89j21.45 ± 0.32lMean3.9419.838.9228.7618.33505.010.3779.92113.62699.5090.5736.37SD0.185.762.907.070.92195.950.1930.2858.05296.2470.2824.51All values are expressed as means of triplicate determinations ± SD. Values in the same column followed by the same letter are not significantly different by Tukey’s multiple range test (p<0.05).


### 3.2. Minerals Content

The mineral contents in honey samples were summarized in Table2; the predominant metal in the majority of honey samples was potassium followed by sodium, calcium, and then magnesium. The same results were obtained by Bouhlali et al. [22] for eleven Moroccan honey from various floral origins.TheArbutus unedo honey sample presented the highest potassium content of 1299.56 ± 3.25 mg/kg, and Salvia rosmarinus honey presented the lowest content (270.57 ± 0.85 mg/kg). The contents in sodium ranged between 41.64 ± 0.54 mg/kg in Salvia rosmarinus honey and 264.56 ± 0.45 mg/kg in Arbutus unedo honey, and the contents in calcium ranged between 22.12 ± 0.21 in multifloral honey sample from Skoura and 283.13 ± 2.36 mg/kg in Arbutus unedo honey, and the contents in magnesium ranged between 13.55 ± 0.56 mg/kg in Ziziphus lotus honey from Marmoucha and 218.72 ± 1.54 mg/kg in multifloral from Elmers. The mineral content in honey is influenced by soil composition, botanical origin, climatic conditions, and seasonal variations [23].Honey is a good source of trace elements that are essential for the proper functioning of the body [24]. Many studies have shown the pharmacological effect of dietary minerals; for instance, it was proven that potassium plays a crucial role in endothelial and cardiovascular function and can reduce blood pressure [25]. Similarly, calcium dietary intake is very important for the good health of the skeleton, the function of skeletal muscles, and nerve conduction [26, 27]. In addition, it was proven in a study conducted by Kh et al. [28] that magnesium supplementation can prevent blood pressure elevation and significantly reduce platelet aggregation.


### 3.3. Polyphenols, Flavones, Flavonols, and Antioxidant Activities

Polyphenols are secondary metabolites widely present in the plant kingdom and known for their pharmacological properties, such as antioxidant, anti-inflammatory, immunomodulatory, and antidiabetic effects [29, 30]. The analysis of polyphenols in honey samples revealed a range between 13.70 ± 0.30 mg GAE/100 g found in the Citrus sinensis honey sample (S5) and 246.20 ± 10.40 mg GAE/100 g in Ceratonia siliqua honey (S12). Flavones and flavonols ranged between 0.70 ± 0.10 mg QE/100 g in Citrus sinensis honey sample (S5) and 31.70 ± 1.70 mg QE/100 g Ceratonia siliqua honey (S12) (Table 3). Our results are higher than those obtained by Petretto et al. [31] for seven commercial Moroccan honey samples from different floral origins and lower than those found in monofloral honey (Bupleurum spinosum) collected from Moroccan Middle Atlas [12]. Moreover, the antioxidant compounds’ rate in honey is affected by several factors such as its floral source, geographical origin, honey maturation processing, handling, and storage [32]. In addition, we observed that honey samples from the same botanical origin (S2 and S17: Ammi visnaga), (S4 and S15: Ziziphus lotus), (S10, S11, S12: Ceratonia siliqua) have different content in antioxidant compounds; the same results were found by Laaroussi et al. [12]. This is explained by the presence of secondary pollen and nectar from other floral sources [33].Table 3
Bioactive compounds and antioxidant activities (DPPH, ABTS, and RP) of honey samples.
SamplePhenols (mg GAE/100 g)Flavones and flavonols (mg QE/100 g)DPPH (IC50 = mg/ml)ABTS (IC50 = mg/ml)RP (EC50 = mg/ml)S189.00 ± 1.40de10.10 ± 0.20e31.20 ± 1.60e11.17 ± 0.30g3.00 ± 0.02hijS263.20 ± 9.60g7.80 ± 1.30ef47.10 ± 2.20d8.76 ± 0.40i3.20 ± 0.01hijS389.10 ± 2.40de11.50 ± 0.30e36.60 ± 1.70e12.49 ± 0.10f5.40 ± 0.40fS4171.70 ± 0.20c23.30 ± 0.00c9.30 ± 1.00h2.30 ± 0.10l1.60 ± 0.01kS513.70 ± 0.30o0.70 ± 0.10h93.40 ± 6.10a45.16 ± 0.70b20.10 ± 0.40aS615.60 ± 1.20n1.00 ± 0.10h62.00 ± 1.50c16.08 ± 0.30c14.20 ± 0.50dS7168.40 ± 1.40c22.90 ± 0.20c36.70 ± 0.40e10.13 ± 0.20g2.10 ± 0.60kS839.40 ± 2.40l4.50 ± 1.10g56.10 ± 8.80c13.43 ± 0.30e12.80 ± 0.02eS9104.10 ± 3.50d13.80 ± 0.60d36.60 ± 3.10e6.88 ± 0.20j3.80 ± 0.30hS1096.30 ± 1.00d9.10 ± 1.00ef15.10 ± 0.70f2.66 ± 0.10l4.00 ± 0.01hS11199.20 ± 8.70b25.90 ± 1.10b14.20 ± 0.20f5.06 ± 0.20k2.00 ± 0.01kS12246.20 ± 10.40a31.70 ± 1.70a12.50 ± 0.50f2.50 ± 0.20l1.90 ± 0.01kS1377.40 ± 0.16def8.80 ± 0.05ef16.29 ± 0.29f14.90 ± 0.13d7.95 ± 0.04gS1417.28 ± 0.71m1.57 ± 0.04h77.65 ± 1.52b58.96 ± 0.92a17.58 ± 0.34bS1573.72 ± 0.84def8.40 ± 0.08ef16.07 ± 0.16f9.39 ± 0.10gh5.54 ± 0.08fS1657.43 ± 0.12k4.23 ± 0.04g30.92 ± 1.34e10.12 ± 0.21g14.62 ± 0.09cS1775.31 ± 0.70def9.76 ± 0.12ef15.30 ± 0.20f3.27 ± 0.25l3.78 ± 0.01hS18102.11 ± 1.03d12.70 ± 0.10d9.52 ± 0.13fg2.93 ± 0.10l3.46 ± 0.01hiS1983.72 ± 1.43def10.11 ± 0.09e10.39 ± 0.11fg5.27 ± 0.32k4.25 ± 0.02hS2098.01 ± 0.28d10.34 ± 0.23e13.85 ± 0.4f2.94 ± 0.163.78 ± 0.01hS2197.56 ± 1.26d11.70 ± 0.22e15.27 ± 0.34f3.67 ± 0.20l4.31 ± 0.01hS2278.20 ± 0.60g9.76 ± 0.08e18.59 ± 0.56f4.70 ± 0.61k7.34 ± 0.03gAll values are expressed as means of triplicate determinations ± SD. Values in the same column followed by the same letter are not significantly different by Tukey’s multiple range test (p<0.05).Concerning the antioxidant activity, the DPPH method is largely used to test the free radical scavenging ability of various samples [20, 34]. The DPPH IC50 values (the concentration with scavenging activity of 50%) of Moroccan honey samples showed significant differences among analyzed samples and ranged between 9.3 ± 1.0 mg/ml in Ziziphus lotus honey sample and 93.40 ± 6.10 mg/ml in Citrus sinensis honey sample.The second method used to evaluate the antioxidant activity of honey samples was the ABTS cation radical assay (ABTS•+). It is based on the interaction between the nitrogen atom of ABTS and hydrogen donating antioxidant; this reaction produces the decolorization of the solution [34]. The IC50 of the ABTS test ranged between 58.96 ± 0.92 mg/ml in the Citrus limon honey sample and 2.30 ± 0.10 mg/ml in the Ziziphus lotus honey sample (Table 3). Furthermore, the reducing power on the ferric ion of honey samples was analyzed; this method is based on the reduction of ferric ion Fe (III) to ferrous ion Fe (II) by the antioxidant compound. The reaction was visualized by the formation of Perl’s Prussian blue complex with maximum absorption at 700 nm [35]. The results of reducing power showed a range of EC50 between 1.60 ± 0.01 mg/ml obtained by Ziziphus lotus honey from Marmoucha and 17.58 ± 0.34 mg/ml obtained by Citrus limon honey from Khnichat.


### 3.4. Proteins and Apalbumin 1 Content

Proteins are one of the minor compounds found in honey. Their percentages vary according to the honeybee origin. The honey produced byApis mellifera contains a range between 0.6% and 1.6% of proteins while the honey produced by Apis cerana contains a range between 0.1% and 3.3% of proteins [36]. It was reported that the amount of protein from bee secretions in honey is higher than that of protein from plants [37]. Among these proteins, there are enzymes responsible for the transformation of nectar components into honey such as glucose-6-oxidase, invertase, and diastase [9]. Furthermore, the major royal jelly protein “apalbumin 1” is one of the most abundant proteins in honey originating from bee secretions, it is an important criterion for honey quality examination, and its detection in honey samples is an indicator of the authenticity and no adulteration [9]. In addition, this protein of 55 kDa had a wide range of biological properties and health-promoting functions such as the immunostimulation effect, the increase of TNF-α release by macrophages, and the antihypertensive activity [38–40].The results of total protein and apalbumin-1 content in honey are summarized in Table4 and Figure 2. According to the results obtained, the range of protein content in honey samples was between 0.021% found in Salvia rosmarinus honey and 0.412% found in Bupleurum spinosum honey. These results are in line with those obtained in other studies [17, 36].Table 4
Total protein and apalbumin1 contents of the honey samples.
The botanical origin of honeyCodeTotal proteinApalbumin1μg/g of honey% of honeyμg/g of honey%of total protein% of honeySalvia rosmarinusS1212.00o0.021o27.40t12.92m0.013mAmmi visnagaS22215.40d0.222d242.60n10.95n0.011mnEuphorbia resiniferaS31245.00i0.125i185.70o14.92k0.015lZiziphus lotusS41245.00i0.125i338.90i27.22fg0.027gCitrus sinensisS5460.80j0.046n137.70r29.88f0.030fCitrus sinensisS6623.00l0.062m131.00s21.03h0.021hMultifloralS72041.60e0.204e277.40l13.58kl0.014mMultifloralS8897.80k0.090k175.90p19.59i0.020hMultifloralS91518.80g0.152g426.70e28.09f0.028gCeratonia siliquaS101659.90f0.166f302.20j18.21j0.018iCeratonia siliquaS112422.40c0.242c295.20k12.19m0.012mCeratonia siliquaS122052.40e0.205e258.10h12.57m0.013iLavandula angustifoliaS13249.60n0.025o157.44q63.08a0.016jCitrus limonS14446.00m0.045n260.05m58.31b0.026hkZiziphus lotusS15840.80j0.084l404.06f48.06d0.040dSilybum marianumS161630.50f0.163f357.89h21.95h0.036eAmmi visnagaS172582.00c0.258c368.91g14.29k0.037eCapparis spinosaS182010.90e0.201e688.34b34.23e0.069bArbutus unedoS191032.40j0.103j545.40c52.83c0.055cBupleurum spinosumS204121.20a0.412a787.73a19.11i0.079aThymus vulgarisS213326.20b0.333b422.87d12.71m0.042dOriganum vulgareS221357.90h0.136h790.82a58.24b0.079aValues in the same column followed by the same letter are not significantly different by Tukey’s multiple range test (p<0.05).Figure 2
Electrophoretic profile of the proteins of honey samples, 12% SDS-PAGE, Coomassie Brilliant Blue staining. Lines RJ: royal jelly (1 mg/ml); Lines 1–22: honey samples (50%); Line Apa: apalbumin 1.For apalbumin 1 analysis, the highest content was shown inOriganum vulgare honey (790.82 μg/g), it represents 58.24% of total protein presented in this honey, followed by Bupleurum spinosum honey (787.73 μg/g), and it represents 19.11% of total protein. Whereas the lowest content of apalbumin 1 was found in Salvia rosmarinus honey (212 μg/g), it represents 12.92% of total protein. As shown in Table 4, the percentage of apalbumin 1 in honey samples ranged between 0.011% and 0.079%; these results are higher than those obtained for honey samples from Italia and Slovakia ranging between 0.08% and 0.03% [17].Overall, the obtained results showed that all analyzed samples contain an important amount of apalbumin 1, the detection of this protein in honey is a criterion of honey’s good quality, and it is an indicator of authenticity and no adulteration [9].


### 3.5. Correlation and Multivariate Analysis

The correlation test between the studied parameters is presented in Table5. The results revealed significant Pearson correlations (p≤0.05) among the different parameters. The antioxidant compounds (polyphenols, flavones, and flavonols) correlate positively with total acidity, ash, Pfund, and the content of sodium and potassium. Polyphenols content correlates negatively with ABTS, DPPH, and RP tests (r = −0.420, r = −0.595, and r = −0.746, respectively). On the other hand, flavones and flavonols correlate negatively with DPPH (r = −0.578) and RP (r = −0.737). RP test correlates positively with DPPH and ABTS (r = 0.553 and r = 0.799, respectively). Total acidity correlates positively with free acidity (r = 0.917) and lactonic acidity (r = 0.614). Electrical conductivity correlates positively with total acidity (r = 0.459), ash (r = 0.765), Pfund (r = 0.611), sodium (r = 0.738), potassium (r = 0.738), calcium (r = 0.562), polyphenols (r = 0.678), and flavones and flavonols (r = 0.685). The content in apalbumin 1 correlates positively with total acidity (r = 0.431), Pfund (r = 0.471), and the content in protein (r = 0.554) while the protein content correlates negatively with pH (r = −0.480) and positively with free acidity (r = 0.554), lactone acidity (r = 0.629), total acidity (r = 0.710), and Pfund (r = 0.570). Furthermore, a principal component analysis (PCA) was applied to the obtained results (Figure 3). PCA is one of the techniques most used for performing multivariate analysis; it is characterized by low difficulty and rapid analysis [41]. Figure 3(a) represents the PCA for physicochemical analysis of the studied honey samples; the PC1 and PC2 explained a variance of 56%. The first component (PC1) explained 39.274% and represents in its positive part all parameters studied except pH that exists in the negative part. The second component (PC2) explained 16.726% and represents in its positive part pH, moisture, Mg, Ca, K, Na, ash, and electrical conductivity, while in the negative part we found Pfund, free acidity, total acidity, and lactonic acidity. The honey samples S8, S9, S7, S4, and S15 shared characteristics regarding pH, moisture, Mg, Ca, K, Na, ash, and electrical conductivity, whereas the honey samples S21 and S18 shared the characteristics for Pfund, free acidity, total acidity, and lactonic acidity.Table 5
Pearson correlation coefficients between the analyzed parameters.
ParameterspHFALATAMoistureECAshPfundNaKCaMgPolyphenolsFlavones and flavonolsABTSDPPHRPApalbumin1ProteinpH10.3130.0830.1240.4670.013−0.055−0.2450.1830.1390.1200.0100.0170.0370.030−0.0890.084−0.145−0.480∗FA−0.22510.2490.917∗∗∗0.2000.4220.3460.4060.452∗0.2430.3310.0390.3330.333−0.536∗−0.107−0.3570.3780.554∗∗LA0.0830.24910.614∗∗−0.0120.2810.1530.493∗0.1360.0560.2580.1900.457∗0.415−0.221−0.365−0.428∗0.3020.629∗∗TA0.1240.917∗∗∗0.614∗∗10.1580.459∗0.3450.533∗0.425∗0.2210.3750.1100.460∗0.442∗−0.527∗−0.238−0.467∗0.431∗0.710∗∗Moisture0.4670.200−0.0120.15810.1780.013−0.1770.218−0.0010.446∗0.163−0.313−0.3160.1710.2630.4010.297−0.051EC0.0130.4220.2810.459∗0.17810.765∗∗0.611∗∗0.738∗∗0.598∗∗0.562∗∗0.2490.678∗∗0.685∗∗−0.338−0.575∗∗−0.648∗∗0.2830.236Ash−0.0550.3460.1530.3450.0130.765∗∗10.3470.561∗∗0.629∗∗0.3180.383∗∗0.485∗0.527∗−0.076−0.501∗∗−0.442∗0.0450.122Pfund−0.2450.4060.493∗0.533∗−0.1770.611∗∗0.34710.3420.432∗0.3570.1560.668∗∗0.626∗∗−0.509∗−0.624∗−0.683∗∗0.471∗0.570∗∗Na0.1830.452∗0.1360.425∗0.2180.738∗∗0.561∗∗0.34210.670∗∗0.445∗0.1560.572∗∗0.581∗∗−0.135−0.455∗−0.453∗−0.0020.098K0.1390.2430.0560.221−0.0010.598∗∗0.629∗∗0.432∗0.670∗∗10.4220.3290.472∗0.501∗−0.222−0.408−0.3770.1800.159Ca0.1200.3310.2580.3750.446∗0.562∗∗0.3180.3570.445∗0.42210.503∗0.2120.230−0.142−0.213−0.1390.2890.137Mg0.0100.0390.1900.1100.1630.2490.383∗∗0.1560.1560.3290.503∗1−0.025−0.0040.221−0.1540.129−0.149−0.117Polyphenols0.0170.3330.457∗0.460∗−0.3130.678∗∗0.485∗0.668∗∗0.572∗∗0.472∗0.212−0.02510.990∗∗−0.420∗−0.595∗−0.746∗∗0.1110.406Flavones and flavonols0.0370.3330.4150.442∗−0.3160.685∗∗0.527∗0.626∗∗0.581∗∗0.501∗0.230−0.0040.990∗∗1−0.373−0.578∗∗−0.737∗∗0.0970.373ABTS0.030−0.536∗−0.221−0.527∗0.171−0.338−0.076−0.509−0.135−0.222−0.1420.221−0.420∗−0.37310.2180.553∗∗−0.581∗∗−0.398DPPH−0.089−0.107−0.365−0.2380.263−0.575∗∗−0.501∗∗−0.624∗−0.455∗−0.408−0.213−0.154−0.595∗−0.578∗∗0.21810.799∗∗∗−0.127−0.341RP0.084−0.357−0.428∗−0.467∗0.401−0.648∗∗−0.442∗−0.683∗∗−0.453∗−0.377−0.1390.129−0.746∗∗−0.737∗∗0.553∗∗0.799∗∗∗1−0.288−0.497∗Apalbumin1−0.1450.3780.3020.431∗0.2970.2830.0450.471∗−0.0020.1800.289−0.1490.1110.097−0.581∗∗−0.127−0.28810.554∗∗Protein−0.480∗0.554∗∗0.629∗∗0.710∗∗−0.0510.2360.1220.570∗∗0.0980.1590.137−0.1170.4060.373−0.398−0.341−0.497∗0.554∗∗1Figure 3
Principal component analysis (PCA). (a) Biplots for physicochemical analysis of honey samples; (b) biplots for antioxidant activities (ABTS, RP, and DPPH), polyphenols, flavones/flavonols, protein, and apalbumin 1 of honey samples. LA: lactonic acidity; TA; total acidity; FA: free acidity; EC: electrical conductivity.
(a)(b)Figure3(b) represents the PCA for polyphenols, flavones/flavonols, the antioxidant activities (DPPH, ABTS, and RP) and the content in total protein and apalbumin 1. The two principal components (PC1 and PC2) explained a variance of 76.973%. The first component explained 55.936% and represents in its positive part apalbumin 1, proteins, polyphenols, and flavones/flavonols, and in its negative part, we found the antioxidant test (DPPH, ABTS, and RP), whereas the second component explained 21.037% and represents in its positive part proteins, apalbumin 1, DPPH, and ABTS and in its negative part polyphenols, flavones/flavonols, and ABTS. S4, S11, and S7 honey samples shared the characteristics for polyphenols and flavones/flavonols content while S15, S17, S18, S19, S20, S21, and S22 honey samples shared the characteristics for apalbumin 1 and proteins. It was suggested that the difference in the physicochemical parameters, minerals content, and bioactive molecules of honey is related to the soil composition, botanical origin, and climatic conditions [32]. In sum, this study supplied new information about the antioxidant activity, protein, and apalbumin contents in Moroccan honey obtained from different botanical origins.


## 3.1. Physicochemical Parameters

Twenty-two honey samples of different botanical and geographical origins of Morocco were analyzed for physicochemical parameters. The results presented in Table2 showed that the pH values ranged between 3.59 ± 0.10 in Euphorbia resinifera honey from Tiznit and 4.30 ± 0.20 in Ziziphus lotus honey from Marmoucha; these values are within the range recommended by Council Directive 2001/110/EC [18]. Free acidity ranged between 11.45 ± 0.12 mEq/kg in Salvia rosmarinus honey and 30.81 ± 1.20 mEq/kg in Thymus vulgaris from Timhdit, lactone acidity ranged between 5.09 ± 0.96 mEq/kg in Ziziphus lotus honey and 17.50 ± 0.22 mEq/kg in Bupleurum spinosum honey, and total acidity ranged between 18.30 ± 0.30 mEq/kg in Citrus sinensis honey and 42.28 ± 0.30 mEq/kg in Ceratonia siliqua honey. These values are in line with those recommended by Codex Alimentarius Commission [19]. For moisture values, all honey samples are below the maximum limit (20%) except the Arbutus unedo honey sample from Tetouan which was slightly higher (20.90 ± 0.15%). The recommended value of moisture in honey indicates its maturity and reflects the desirable density [20]. For electrical conductivity, all honey samples are bellowing the maximum limit allowed (800 μS/cm) except Ceratonia siliqua honey from Taounate (948.33 ± 1.63 μS/cm) and Arbutus unedo honey from Tetouan (913.51 ± 6.21 μS/cm). Similarly, the ash contents in all honey samples are below the limits (0.6%) except for the multifloral honey sample from Skoura (0.92 ± 0.02%). The ash content reflects the mineral composition in honey. It is influenced by soil and botanical origins [12]. On the other hand, Pfund is a parameter that reflects honey color, the analysis of Pfund showed values ranging between 19.67 ± 0.64 mm in Citrus limon honey and 128.74 ± 1.20 mm in Ceratonia siliqua honey. According to the Pfund scale, the color of our honey samples ranged between white color and dark amber [21].Table 2
Physicochemical parameters and minerals content in honey samples.
SamplepHFree acidity (mEq/kg)Lactone acidity (mEq/kg)Total acidity (mEq/kg)Moisture (%)Electrical conductivity (mS/cm)Ash%Pfund (mm)Na (mm/kg)K (mm/kg)Ca (mm/kg)Mg (mm/kg)S14.05 ± 0.03b11.45 ± 0.12hi9.51 ± 0.50cd20.96 ± 0.15j17.20 ± 0.22bc443.00 ± 20.00h0.28 ± 0.01g83.06 ± 0.01j41.64 ± 0.54o270.57 ± 0.85s28.49 ± 0.21o19.76 ± 1.25nS23.71 ± 0.02b21.00 ± 1.00d11.10 ± 0.10cd32.10 ± 0.35e18.60 ± 0.41b472.33 ± 3.22g0.31 ± 0.01ef64.21 ± 0.10m128.16 ± 0.4f286.78 ± 0.42r25.98 ± 1.20o17.44 ± 0.56oS33.59 ± 0.10b20.51 ± 0.50d8.10 ± 0.36f28.61 ± 0.52g17.50 ± 0.50bc444.00 ± 20.00h0.33 ± 0.02ef75.01 ± 0.12l117.89 ± 1.21h667.46 ± 0.78h77.27 ± 1.35i42.36 ± 1.25eS44.30 ± 0.20a14.41 ± 0.90g7.51 ± 0.12f21.92 ± 0.45j17.70 ± 0.23bc620.83 ± 1.53c0.42 ± 0.05d63.43 ± 0.11m154.35 ± 0.54d1014.28 ± 2.36c67.57 ± 0.78k13.55 ± 0.56pS54.07 ± 0.04b12.00 ± 0.15hi6.30 ± 0.22f18.30 ± 0.30k19.80 ± 0.12a242.00 ± 30.00m0.17 ± 0.01gh25.58 ± 0.23q114.74 ± 0.41i424.24 ± 1.25o79.35 ± 1.20h25.76 ± 0.87jS64.20 ± 0.10a13.51 ± 0.60h6.00 ± 0.12f19.51 ± 0.40k19.50 ± 0.52a223.00 ± 40.00n0.20 ± 0.01g33.13 ± 0.13p71.48 ± 0.52l622.95 ± 2.32j86.36 ± 0.36g28.15 ± 0.47hS73.94 ± 0.20a14.45 ± 0.12g7.41 ± 0.52f21.86 ± 0.40j17.10 ± 0.36d577.33 ± 4.22d0.35 ± 0.01e103.77 ± 0.62e114.63 ± 0.21i947.36 ± 3.25e93.21 ± 0.98f35.13 ± 0.36fS83.94 ± 0.40a16.05 ± 0.80g10.11 ± 0.63cd26.16 ± 0.30h18.80 ± 0.20b510.67 ± 1.23e0.55 ± 0.02c78.74 ± 0.14k78.16 ± 1.21k847.21 ± 1.25f167.57 ± 1.25d218.72 ± 1.54aS93.85 ± 0.00b21.51 ± 0.16d8.20 ± 0.10f29.71 ± 0.62f18.70 ± 0.12b762.33 ± 2.33b0.92 ± 0.02a76.19 ± 0.20l153.04 ± 1.54d987.89 ± 1.25d22.12 ± 0.21p34.35 ± 0.54gS103.75 ± 0.10b11.45 ± 0.90hi9.81 ± 0.41cd21.26 ± 0.16j17.10 ± 0.32bc484.33 ± 4.53f0.31 ± 0.01ef102.92 ± 0.35e131.15 ± 0.65e1016.13 ± 4.23c26.00 ± 1.23o28.07 ± 0.88hS113.79 ± 0.08b26.10 ± 0.60bc12.15 ± 0.10c38.25 ± 0.55bc18.80 ± 0.22b948.33 ± 1.63a0.74 ± 0.01b108.91 ± 0.35c174.24 ± 2.87c806.56 ± 1.23g252.43 ± 0.508b60.97 ± 0.58dS124.00 ± 0.01a28.21 ± 0.23b14.00 ± 0.36b42.28 ± 0.30a18.50 ± 0.52b700.00 ± 40.00b0.52 ± 0.04c128.74 ± 1.20a233.92 ± 1.32b1017.98 ± 3.24c80.3 ± 0.69i62.52 ± 0.41cS134.18 ± 0.12a18.15 ± 0.90f8.03 ± 0.09f26.18 ± 0.22hi17.30 ± 0.23bc345.44 ± 3.12k0.31 ± 0.05e40.83 ± 0.52o104.77 ± 0.65j588.92 ± 1.25l59.44 ± 1.25l40.17 ± 0.78eS143.80 ± 0.20b20.56 ± 0.12d6.70 ± 0.23f27.26 ± 0.30h18.74 ± 0.62b231.83 ± 2.21n0.13 ± 0.07gh19.67 ± 0.64r41.66 ± 0.36o331.89 ± 2.32p37.50 ± 1.80n19.44 ± 0.44lS154.08 ± 0.02a21.47 ± 0.40d5.09 ± 0.96g26.56 ± 0.52h19.28 ± 0.20a521.60 ± 2.36e0.42 ± 0.09e56.20 ± 0.12n121.30 ± 0.21g616.38 ± 1.23k77.90 ± 0.36i24.02 ± 0.21jS163.82 ± 0.01b18.05 ± 0.22f7.15 ± 0.52f25.20 ± 0.41hi19.15 ± 0.15a329.70 ± 5.21l0.20 ± 0.01g87.34 ± 0.69i43.30 ± 0.22o494.05 ± 0.23n42.05 ± 0.25m20.24 ± 0.14lS173.99 ± 0.01a28.43 ± 0.12b6.85 ± 0.23f35.28 ± 0.30d18.15 ± 0.52b489.99 ± 5.21f0.21 ± 0.06g97.24 ± 0.54f112.75 ± 0.25i303.67 ± 1.24q43.97 ± 0.87m17.56 ± 0.25mS183.94 ± 0.02a20.04 ± 0.25de11.37 ± 0.63c31.41 ± 0.40e17.75 ± 0.63bc428.89 ± 6.12i0.29 ± 0.04g119.65 ± 0.41b78.42 ± 1.22k585.99 ± 1.23l172.42 ± 1.25c36.80 ± 0.87fS194.14 ± 0.05a27.40 ± 0.63b7.70 ± 0.52f35.10 ± 0.23d20.90 ± 0.15a913.51 ± 6.21b0.49 ± 0.04d103.75 ± 0.25e264.56 ± 0.45a1299.56 ± 3.25a283.13 ± 2.36a77.23 ± 0.59bS203.76 ± 0.06b22.05 ± 0.41d17.50 ± 0.22a39.55 ± 0.60b19.67 ± 0.34a466.67 ± 3.21g0.20 ± 0.05g89.13 ± 0.74h64.56 ± 0.98m521.26 ± 0.98m121.47 ± 1.54e26.42 ± 0.54iS213.82 ± 0.02b30.81 ± 1.20a8.45 ± 0.40e39.26 ± 0.52b17.32 ± 0.42bc372.04 ± 2.35j0.48 ± 0.05d94.18 ± 0.98g103.29 ± 1.02j1105.58 ± 2.36b76.26 ± 0.65i30.10 ± 0.86hS223.94 ± 0.04a18.70 ± 0.17f7.25 ± 0.54f25.95 ± 0.45hi19.00 ± 0.56a582.35 ± 4.21d0.32 ± 0.01e106.63 ± 1.21d51.60 ± 0.45n632.20 ± 3.20i71.85 ± 0.89j21.45 ± 0.32lMean3.9419.838.9228.7618.33505.010.3779.92113.62699.5090.5736.37SD0.185.762.907.070.92195.950.1930.2858.05296.2470.2824.51All values are expressed as means of triplicate determinations ± SD. Values in the same column followed by the same letter are not significantly different by Tukey’s multiple range test (p<0.05).


## 3.2. Minerals Content

The mineral contents in honey samples were summarized in Table2; the predominant metal in the majority of honey samples was potassium followed by sodium, calcium, and then magnesium. The same results were obtained by Bouhlali et al. [22] for eleven Moroccan honey from various floral origins.TheArbutus unedo honey sample presented the highest potassium content of 1299.56 ± 3.25 mg/kg, and Salvia rosmarinus honey presented the lowest content (270.57 ± 0.85 mg/kg). The contents in sodium ranged between 41.64 ± 0.54 mg/kg in Salvia rosmarinus honey and 264.56 ± 0.45 mg/kg in Arbutus unedo honey, and the contents in calcium ranged between 22.12 ± 0.21 in multifloral honey sample from Skoura and 283.13 ± 2.36 mg/kg in Arbutus unedo honey, and the contents in magnesium ranged between 13.55 ± 0.56 mg/kg in Ziziphus lotus honey from Marmoucha and 218.72 ± 1.54 mg/kg in multifloral from Elmers. The mineral content in honey is influenced by soil composition, botanical origin, climatic conditions, and seasonal variations [23].Honey is a good source of trace elements that are essential for the proper functioning of the body [24]. Many studies have shown the pharmacological effect of dietary minerals; for instance, it was proven that potassium plays a crucial role in endothelial and cardiovascular function and can reduce blood pressure [25]. Similarly, calcium dietary intake is very important for the good health of the skeleton, the function of skeletal muscles, and nerve conduction [26, 27]. In addition, it was proven in a study conducted by Kh et al. [28] that magnesium supplementation can prevent blood pressure elevation and significantly reduce platelet aggregation.


## 3.3. Polyphenols, Flavones, Flavonols, and Antioxidant Activities

Polyphenols are secondary metabolites widely present in the plant kingdom and known for their pharmacological properties, such as antioxidant, anti-inflammatory, immunomodulatory, and antidiabetic effects [29, 30]. The analysis of polyphenols in honey samples revealed a range between 13.70 ± 0.30 mg GAE/100 g found in the Citrus sinensis honey sample (S5) and 246.20 ± 10.40 mg GAE/100 g in Ceratonia siliqua honey (S12). Flavones and flavonols ranged between 0.70 ± 0.10 mg QE/100 g in Citrus sinensis honey sample (S5) and 31.70 ± 1.70 mg QE/100 g Ceratonia siliqua honey (S12) (Table 3). Our results are higher than those obtained by Petretto et al. [31] for seven commercial Moroccan honey samples from different floral origins and lower than those found in monofloral honey (Bupleurum spinosum) collected from Moroccan Middle Atlas [12]. Moreover, the antioxidant compounds’ rate in honey is affected by several factors such as its floral source, geographical origin, honey maturation processing, handling, and storage [32]. In addition, we observed that honey samples from the same botanical origin (S2 and S17: Ammi visnaga), (S4 and S15: Ziziphus lotus), (S10, S11, S12: Ceratonia siliqua) have different content in antioxidant compounds; the same results were found by Laaroussi et al. [12]. This is explained by the presence of secondary pollen and nectar from other floral sources [33].Table 3
Bioactive compounds and antioxidant activities (DPPH, ABTS, and RP) of honey samples.
SamplePhenols (mg GAE/100 g)Flavones and flavonols (mg QE/100 g)DPPH (IC50 = mg/ml)ABTS (IC50 = mg/ml)RP (EC50 = mg/ml)S189.00 ± 1.40de10.10 ± 0.20e31.20 ± 1.60e11.17 ± 0.30g3.00 ± 0.02hijS263.20 ± 9.60g7.80 ± 1.30ef47.10 ± 2.20d8.76 ± 0.40i3.20 ± 0.01hijS389.10 ± 2.40de11.50 ± 0.30e36.60 ± 1.70e12.49 ± 0.10f5.40 ± 0.40fS4171.70 ± 0.20c23.30 ± 0.00c9.30 ± 1.00h2.30 ± 0.10l1.60 ± 0.01kS513.70 ± 0.30o0.70 ± 0.10h93.40 ± 6.10a45.16 ± 0.70b20.10 ± 0.40aS615.60 ± 1.20n1.00 ± 0.10h62.00 ± 1.50c16.08 ± 0.30c14.20 ± 0.50dS7168.40 ± 1.40c22.90 ± 0.20c36.70 ± 0.40e10.13 ± 0.20g2.10 ± 0.60kS839.40 ± 2.40l4.50 ± 1.10g56.10 ± 8.80c13.43 ± 0.30e12.80 ± 0.02eS9104.10 ± 3.50d13.80 ± 0.60d36.60 ± 3.10e6.88 ± 0.20j3.80 ± 0.30hS1096.30 ± 1.00d9.10 ± 1.00ef15.10 ± 0.70f2.66 ± 0.10l4.00 ± 0.01hS11199.20 ± 8.70b25.90 ± 1.10b14.20 ± 0.20f5.06 ± 0.20k2.00 ± 0.01kS12246.20 ± 10.40a31.70 ± 1.70a12.50 ± 0.50f2.50 ± 0.20l1.90 ± 0.01kS1377.40 ± 0.16def8.80 ± 0.05ef16.29 ± 0.29f14.90 ± 0.13d7.95 ± 0.04gS1417.28 ± 0.71m1.57 ± 0.04h77.65 ± 1.52b58.96 ± 0.92a17.58 ± 0.34bS1573.72 ± 0.84def8.40 ± 0.08ef16.07 ± 0.16f9.39 ± 0.10gh5.54 ± 0.08fS1657.43 ± 0.12k4.23 ± 0.04g30.92 ± 1.34e10.12 ± 0.21g14.62 ± 0.09cS1775.31 ± 0.70def9.76 ± 0.12ef15.30 ± 0.20f3.27 ± 0.25l3.78 ± 0.01hS18102.11 ± 1.03d12.70 ± 0.10d9.52 ± 0.13fg2.93 ± 0.10l3.46 ± 0.01hiS1983.72 ± 1.43def10.11 ± 0.09e10.39 ± 0.11fg5.27 ± 0.32k4.25 ± 0.02hS2098.01 ± 0.28d10.34 ± 0.23e13.85 ± 0.4f2.94 ± 0.163.78 ± 0.01hS2197.56 ± 1.26d11.70 ± 0.22e15.27 ± 0.34f3.67 ± 0.20l4.31 ± 0.01hS2278.20 ± 0.60g9.76 ± 0.08e18.59 ± 0.56f4.70 ± 0.61k7.34 ± 0.03gAll values are expressed as means of triplicate determinations ± SD. Values in the same column followed by the same letter are not significantly different by Tukey’s multiple range test (p<0.05).Concerning the antioxidant activity, the DPPH method is largely used to test the free radical scavenging ability of various samples [20, 34]. The DPPH IC50 values (the concentration with scavenging activity of 50%) of Moroccan honey samples showed significant differences among analyzed samples and ranged between 9.3 ± 1.0 mg/ml in Ziziphus lotus honey sample and 93.40 ± 6.10 mg/ml in Citrus sinensis honey sample.The second method used to evaluate the antioxidant activity of honey samples was the ABTS cation radical assay (ABTS•+). It is based on the interaction between the nitrogen atom of ABTS and hydrogen donating antioxidant; this reaction produces the decolorization of the solution [34]. The IC50 of the ABTS test ranged between 58.96 ± 0.92 mg/ml in the Citrus limon honey sample and 2.30 ± 0.10 mg/ml in the Ziziphus lotus honey sample (Table 3). Furthermore, the reducing power on the ferric ion of honey samples was analyzed; this method is based on the reduction of ferric ion Fe (III) to ferrous ion Fe (II) by the antioxidant compound. The reaction was visualized by the formation of Perl’s Prussian blue complex with maximum absorption at 700 nm [35]. The results of reducing power showed a range of EC50 between 1.60 ± 0.01 mg/ml obtained by Ziziphus lotus honey from Marmoucha and 17.58 ± 0.34 mg/ml obtained by Citrus limon honey from Khnichat.


## 3.4. Proteins and Apalbumin 1 Content

Proteins are one of the minor compounds found in honey. Their percentages vary according to the honeybee origin. The honey produced byApis mellifera contains a range between 0.6% and 1.6% of proteins while the honey produced by Apis cerana contains a range between 0.1% and 3.3% of proteins [36]. It was reported that the amount of protein from bee secretions in honey is higher than that of protein from plants [37]. Among these proteins, there are enzymes responsible for the transformation of nectar components into honey such as glucose-6-oxidase, invertase, and diastase [9]. Furthermore, the major royal jelly protein “apalbumin 1” is one of the most abundant proteins in honey originating from bee secretions, it is an important criterion for honey quality examination, and its detection in honey samples is an indicator of the authenticity and no adulteration [9]. In addition, this protein of 55 kDa had a wide range of biological properties and health-promoting functions such as the immunostimulation effect, the increase of TNF-α release by macrophages, and the antihypertensive activity [38–40].The results of total protein and apalbumin-1 content in honey are summarized in Table4 and Figure 2. According to the results obtained, the range of protein content in honey samples was between 0.021% found in Salvia rosmarinus honey and 0.412% found in Bupleurum spinosum honey. These results are in line with those obtained in other studies [17, 36].Table 4
Total protein and apalbumin1 contents of the honey samples.
The botanical origin of honeyCodeTotal proteinApalbumin1μg/g of honey% of honeyμg/g of honey%of total protein% of honeySalvia rosmarinusS1212.00o0.021o27.40t12.92m0.013mAmmi visnagaS22215.40d0.222d242.60n10.95n0.011mnEuphorbia resiniferaS31245.00i0.125i185.70o14.92k0.015lZiziphus lotusS41245.00i0.125i338.90i27.22fg0.027gCitrus sinensisS5460.80j0.046n137.70r29.88f0.030fCitrus sinensisS6623.00l0.062m131.00s21.03h0.021hMultifloralS72041.60e0.204e277.40l13.58kl0.014mMultifloralS8897.80k0.090k175.90p19.59i0.020hMultifloralS91518.80g0.152g426.70e28.09f0.028gCeratonia siliquaS101659.90f0.166f302.20j18.21j0.018iCeratonia siliquaS112422.40c0.242c295.20k12.19m0.012mCeratonia siliquaS122052.40e0.205e258.10h12.57m0.013iLavandula angustifoliaS13249.60n0.025o157.44q63.08a0.016jCitrus limonS14446.00m0.045n260.05m58.31b0.026hkZiziphus lotusS15840.80j0.084l404.06f48.06d0.040dSilybum marianumS161630.50f0.163f357.89h21.95h0.036eAmmi visnagaS172582.00c0.258c368.91g14.29k0.037eCapparis spinosaS182010.90e0.201e688.34b34.23e0.069bArbutus unedoS191032.40j0.103j545.40c52.83c0.055cBupleurum spinosumS204121.20a0.412a787.73a19.11i0.079aThymus vulgarisS213326.20b0.333b422.87d12.71m0.042dOriganum vulgareS221357.90h0.136h790.82a58.24b0.079aValues in the same column followed by the same letter are not significantly different by Tukey’s multiple range test (p<0.05).Figure 2
Electrophoretic profile of the proteins of honey samples, 12% SDS-PAGE, Coomassie Brilliant Blue staining. Lines RJ: royal jelly (1 mg/ml); Lines 1–22: honey samples (50%); Line Apa: apalbumin 1.For apalbumin 1 analysis, the highest content was shown inOriganum vulgare honey (790.82 μg/g), it represents 58.24% of total protein presented in this honey, followed by Bupleurum spinosum honey (787.73 μg/g), and it represents 19.11% of total protein. Whereas the lowest content of apalbumin 1 was found in Salvia rosmarinus honey (212 μg/g), it represents 12.92% of total protein. As shown in Table 4, the percentage of apalbumin 1 in honey samples ranged between 0.011% and 0.079%; these results are higher than those obtained for honey samples from Italia and Slovakia ranging between 0.08% and 0.03% [17].Overall, the obtained results showed that all analyzed samples contain an important amount of apalbumin 1, the detection of this protein in honey is a criterion of honey’s good quality, and it is an indicator of authenticity and no adulteration [9].


## 3.5. Correlation and Multivariate Analysis

The correlation test between the studied parameters is presented in Table5. The results revealed significant Pearson correlations (p≤0.05) among the different parameters. The antioxidant compounds (polyphenols, flavones, and flavonols) correlate positively with total acidity, ash, Pfund, and the content of sodium and potassium. Polyphenols content correlates negatively with ABTS, DPPH, and RP tests (r = −0.420, r = −0.595, and r = −0.746, respectively). On the other hand, flavones and flavonols correlate negatively with DPPH (r = −0.578) and RP (r = −0.737). RP test correlates positively with DPPH and ABTS (r = 0.553 and r = 0.799, respectively). Total acidity correlates positively with free acidity (r = 0.917) and lactonic acidity (r = 0.614). Electrical conductivity correlates positively with total acidity (r = 0.459), ash (r = 0.765), Pfund (r = 0.611), sodium (r = 0.738), potassium (r = 0.738), calcium (r = 0.562), polyphenols (r = 0.678), and flavones and flavonols (r = 0.685). The content in apalbumin 1 correlates positively with total acidity (r = 0.431), Pfund (r = 0.471), and the content in protein (r = 0.554) while the protein content correlates negatively with pH (r = −0.480) and positively with free acidity (r = 0.554), lactone acidity (r = 0.629), total acidity (r = 0.710), and Pfund (r = 0.570). Furthermore, a principal component analysis (PCA) was applied to the obtained results (Figure 3). PCA is one of the techniques most used for performing multivariate analysis; it is characterized by low difficulty and rapid analysis [41]. Figure 3(a) represents the PCA for physicochemical analysis of the studied honey samples; the PC1 and PC2 explained a variance of 56%. The first component (PC1) explained 39.274% and represents in its positive part all parameters studied except pH that exists in the negative part. The second component (PC2) explained 16.726% and represents in its positive part pH, moisture, Mg, Ca, K, Na, ash, and electrical conductivity, while in the negative part we found Pfund, free acidity, total acidity, and lactonic acidity. The honey samples S8, S9, S7, S4, and S15 shared characteristics regarding pH, moisture, Mg, Ca, K, Na, ash, and electrical conductivity, whereas the honey samples S21 and S18 shared the characteristics for Pfund, free acidity, total acidity, and lactonic acidity.Table 5
Pearson correlation coefficients between the analyzed parameters.
ParameterspHFALATAMoistureECAshPfundNaKCaMgPolyphenolsFlavones and flavonolsABTSDPPHRPApalbumin1ProteinpH10.3130.0830.1240.4670.013−0.055−0.2450.1830.1390.1200.0100.0170.0370.030−0.0890.084−0.145−0.480∗FA−0.22510.2490.917∗∗∗0.2000.4220.3460.4060.452∗0.2430.3310.0390.3330.333−0.536∗−0.107−0.3570.3780.554∗∗LA0.0830.24910.614∗∗−0.0120.2810.1530.493∗0.1360.0560.2580.1900.457∗0.415−0.221−0.365−0.428∗0.3020.629∗∗TA0.1240.917∗∗∗0.614∗∗10.1580.459∗0.3450.533∗0.425∗0.2210.3750.1100.460∗0.442∗−0.527∗−0.238−0.467∗0.431∗0.710∗∗Moisture0.4670.200−0.0120.15810.1780.013−0.1770.218−0.0010.446∗0.163−0.313−0.3160.1710.2630.4010.297−0.051EC0.0130.4220.2810.459∗0.17810.765∗∗0.611∗∗0.738∗∗0.598∗∗0.562∗∗0.2490.678∗∗0.685∗∗−0.338−0.575∗∗−0.648∗∗0.2830.236Ash−0.0550.3460.1530.3450.0130.765∗∗10.3470.561∗∗0.629∗∗0.3180.383∗∗0.485∗0.527∗−0.076−0.501∗∗−0.442∗0.0450.122Pfund−0.2450.4060.493∗0.533∗−0.1770.611∗∗0.34710.3420.432∗0.3570.1560.668∗∗0.626∗∗−0.509∗−0.624∗−0.683∗∗0.471∗0.570∗∗Na0.1830.452∗0.1360.425∗0.2180.738∗∗0.561∗∗0.34210.670∗∗0.445∗0.1560.572∗∗0.581∗∗−0.135−0.455∗−0.453∗−0.0020.098K0.1390.2430.0560.221−0.0010.598∗∗0.629∗∗0.432∗0.670∗∗10.4220.3290.472∗0.501∗−0.222−0.408−0.3770.1800.159Ca0.1200.3310.2580.3750.446∗0.562∗∗0.3180.3570.445∗0.42210.503∗0.2120.230−0.142−0.213−0.1390.2890.137Mg0.0100.0390.1900.1100.1630.2490.383∗∗0.1560.1560.3290.503∗1−0.025−0.0040.221−0.1540.129−0.149−0.117Polyphenols0.0170.3330.457∗0.460∗−0.3130.678∗∗0.485∗0.668∗∗0.572∗∗0.472∗0.212−0.02510.990∗∗−0.420∗−0.595∗−0.746∗∗0.1110.406Flavones and flavonols0.0370.3330.4150.442∗−0.3160.685∗∗0.527∗0.626∗∗0.581∗∗0.501∗0.230−0.0040.990∗∗1−0.373−0.578∗∗−0.737∗∗0.0970.373ABTS0.030−0.536∗−0.221−0.527∗0.171−0.338−0.076−0.509−0.135−0.222−0.1420.221−0.420∗−0.37310.2180.553∗∗−0.581∗∗−0.398DPPH−0.089−0.107−0.365−0.2380.263−0.575∗∗−0.501∗∗−0.624∗−0.455∗−0.408−0.213−0.154−0.595∗−0.578∗∗0.21810.799∗∗∗−0.127−0.341RP0.084−0.357−0.428∗−0.467∗0.401−0.648∗∗−0.442∗−0.683∗∗−0.453∗−0.377−0.1390.129−0.746∗∗−0.737∗∗0.553∗∗0.799∗∗∗1−0.288−0.497∗Apalbumin1−0.1450.3780.3020.431∗0.2970.2830.0450.471∗−0.0020.1800.289−0.1490.1110.097−0.581∗∗−0.127−0.28810.554∗∗Protein−0.480∗0.554∗∗0.629∗∗0.710∗∗−0.0510.2360.1220.570∗∗0.0980.1590.137−0.1170.4060.373−0.398−0.341−0.497∗0.554∗∗1Figure 3
Principal component analysis (PCA). (a) Biplots for physicochemical analysis of honey samples; (b) biplots for antioxidant activities (ABTS, RP, and DPPH), polyphenols, flavones/flavonols, protein, and apalbumin 1 of honey samples. LA: lactonic acidity; TA; total acidity; FA: free acidity; EC: electrical conductivity.
(a)(b)Figure3(b) represents the PCA for polyphenols, flavones/flavonols, the antioxidant activities (DPPH, ABTS, and RP) and the content in total protein and apalbumin 1. The two principal components (PC1 and PC2) explained a variance of 76.973%. The first component explained 55.936% and represents in its positive part apalbumin 1, proteins, polyphenols, and flavones/flavonols, and in its negative part, we found the antioxidant test (DPPH, ABTS, and RP), whereas the second component explained 21.037% and represents in its positive part proteins, apalbumin 1, DPPH, and ABTS and in its negative part polyphenols, flavones/flavonols, and ABTS. S4, S11, and S7 honey samples shared the characteristics for polyphenols and flavones/flavonols content while S15, S17, S18, S19, S20, S21, and S22 honey samples shared the characteristics for apalbumin 1 and proteins. It was suggested that the difference in the physicochemical parameters, minerals content, and bioactive molecules of honey is related to the soil composition, botanical origin, and climatic conditions [32]. In sum, this study supplied new information about the antioxidant activity, protein, and apalbumin contents in Moroccan honey obtained from different botanical origins.


## 4. Conclusions

Our results assert that Moroccan honey samples are rich in bioactive compounds such as polyphenols, flavones/flavonols, and proteins and are endowed with great antioxidant activities. Nevertheless, the detection of the major royal jelly protein apalbumin 1 in all analyzed samples confirms the quality and no adulteration of Moroccan honey.



---

*Source: 6045792-2022-09-20.xml*