Next Article in Journal
Effects of Four Strains of Actinomycetes on the Content of Terpenoids in Baijiu
Next Article in Special Issue
Antioxidants Discovery for Differentiation of Monofloral Stingless Bee Honeys Using Ambient Mass Spectrometry and Metabolomics Approaches
Previous Article in Journal
College Students’ Preferences for Milk Tea: Results from a Choice Experiment
Previous Article in Special Issue
Influence of the Physical State of Two Monofloral Honeys on Sensory Properties and Consumer Satisfaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multivariate Statistical Approach for the Discrimination of Honey Samples from Galicia (NW Spain) Using Physicochemical and Pollen Parameters

by
Olga Escuredo
1,*,
María Shantal Rodríguez-Flores
1,
Montserrat Míguez
2 and
María Carmen Seijo
1
1
Department of Vegetal Biology and Soil Sciences, Faculty of Sciences, University of Vigo, As Lagoas, 32004 Ourense, Spain
2
Department Analytical and Food Chemistry, Faculty of Sciences, University of Vigo, As Lagoas, 32004 Ourense, Spain
*
Author to whom correspondence should be addressed.
Submission received: 28 February 2023 / Revised: 29 March 2023 / Accepted: 30 March 2023 / Published: 1 April 2023

Abstract

:
Raw honey is a food with a close relation to the territory in which it is produced because of factors such as soil conditions, weather patterns, and plant communities living in the area together. Furthermore, beekeeping management affects the properties of honey. Protected Geographical Indication Miel de Galicia protects the honey produced in Galicia (Northwest Spain). Various types of honeys (362 samples) from this geographical area were analyzed using chemometric techniques. Principal component analysis was favorable to analyzing the physicochemical and pollen variables with the greatest weight in the differentiation of honey. The linear discriminant analysis correctly classified 89.8% of the samples according to the botanical origin using main pollen spectra and physicochemical attributes (moisture, pH, electrical conductivity, diastase content, phenols, flavonoids, and color). Regarding unifloral honey, blackberry, eucalyptus, and heather honeys were correctly grouped, while five chestnut honeys and fourteen samples of honeydew honeys were misclassified. The chestnut and honeydew honeys have similar physicochemical properties and frequently similar pollen spectra profiles complicating the differentiation. Experimental evidence suggests the potential of multivariate statistics in the characterization of honey of the same geographical origin. Therefore, the classification results were good, with electrical conductivity, total phenol content, total flavonoid content and dominant pollens Eucalyptus, Erica, Rubus and Castanea sativa as the variables of higher importance in the differentiation of botanical origin of honeys.

1. Introduction

Honey is considered a natural and animal product because it is produced by honey bees, but the main sources are plant secretions such as nectar (known as a blossom or floral honey) or honeydew harvested on the living parts of plants (honeydew honey). In the hive, the bees transform this mixture of different substances composed mainly of sugars, water, proteins, and other compounds such as enzymes, organic acids, vitamins, and phenols in honey [1,2,3]. Glucose and fructose are the carbohydrates present in the greatest amount and contribute mainly to the energy value and physical characteristics of honey. Other minor constituents of honey are proteins, enzymes (such as diastase, invertase, and glucose-oxidase), amino acids and organic acids (proline and gluconic acid, as the most abundant, respectively), vitamins (mainly vitamin C and group B), and minerals (K, Ca, Na and P, among others) [1,4]. Honey contains diverse phytochemical substances biosynthesized by plants that have important antioxidant activity [5]. Among them, phenol acids and flavonoids play a significant role in their bio-functional properties, as well as in their antioxidant and anti-inflammatory activities [6,7]. Some findings also confirmed the contribution of some aliphatic acids from royal jelly in the antibacterial action of honey [8]. Hence, honey is recognized as one of the most popular natural functional foods, which is attributed to many healthy properties. Therefore, this bee product is used in a wide variety of nutritional and medicinal products and cosmetic applications [2,3,9]. Moreover, it has been used as a food ingredient in interesting culinary arts for its attractive physical and organoleptic qualities [4].
The quantitative variability of each chemical compound in the composition of bee products defines the versatility of its physical, sensory, and functional characteristics [6]. The chemical constituents of honey, as well as the viscosity of its texture, are variable and primarily depend on the floral source, but some external factors, such as seasonal and environmental factors and the processing method, are also influential [1,7,10,11]. One of the most affected properties is color because the nectar origin and its composition contribute to the wide diversity perceptions, from almost water white to black, that can be seen according to the honey type. Thus, some compounds are related to a darker hue of honey, and this is the case of minerals, phenol acids and flavonoids [12,13,14]. The importance of the color of honey lies in consumer preferences and the choice for one type of honey or another, and consequently, it is closely linked to the commercial value of the product. Generally, lighter honeys fetch higher prices on the market, but these preferences vary by country and region. In recent years, dark honey types such as honeydew honey have been gaining the market in countries such as Germany, Switzerland, Greece, and Turkey [4], as well as in Spain. Hence, the importance of their study relating this honey type with the area of production and the source of honeydew.
The relationship of the honey composition with its origin is a relevant issue since it is the tool used to assess the authenticity of honey [7,9]. Several researchers showed differences in the characteristics of honey according to their floral and geographical origins [1,11,14,15,16,17,18]. In Spain, the production of a wide range of honeys with different physicochemical and antioxidant characteristics has been reported depending on the production area [1,12,14,19,20,21,22], with a distinguishable unifloral gap between the north and the south of the country. The most representative unifloral honeys in south-central Spain are orange blossom honey (Citrus aurantium), holm oak or oak honey (Quercus ilex), lavender honey (Lavandula angustifolia), sunflower honey (Helianthus annuus), rosemary honey (Rosmarinus officinalis), and thyme honey (Thymus). In the northern area of the Iberian Peninsula, chestnut honey (Castanea sativa), blackberry honey (Rubus), eucalyptus honey (Eucalyptus), and honeydew honey (Quercus pyrenaica) are produced. Heather honey obtained from a different Erica species or from Calluna vulgaris can also be found throughout the Spanish territory.
The first honey quality criteria observed by consumers are the color and the density or degree of crystallization. Therefore, they are the determining parameters in the selection of the purchase. In this sense, beekeepers, retailers, consumers, and food safety regulatory offices are interested in knowing the quality of honey based on botanical origin and its particular physicochemical properties [1]. Unfortunately, honey products are often adulterated and mislabeled. Honey can be adulterated with cheaper and poorer quality honey, sugar solutions or syrups with high fructose content. The composition of honey can also be manipulated by feeding the bees with sugars and syrups, producing artificial honey [23]. Sometimes, the labelling is changed to a specific botanical or geographical origin to increase the commercial value in the market. For this reason, the adulteration of honey has become a global problem, which in addition to affecting the consumer, undermines the efforts of many beekeepers in the struggle to produce excellent products with qualities worthy of recognition. Honey authentication, regarding its genuine botanical and geographical origins, as well as the detection of any adulteration, is essential in order to protect consumer health and avoid unfair competition that could create a destabilized market [9].
The Protected Geographical Indication (PGI) Miel de Galicia is a designation of origin registered in the European Union since 2007. It covered a territory of about 29,500 km2, with traditional apiculture mainly based on non-migratory apiaries, producing more than 500 t of honey under the label [24]. This quality designation is crucial for the development of beekeeping in the area and the increase of the economic value of honey. This is also in addition to the close contribution to the maintenance of biodiversity and rural environments. Therefore, this type of research can collaborate with the strategy of linking honey and territory throughout the botanical origin and key physicochemical parameters. The present study is the first work that includes such a large number of samples from the same geographical territory for this purpose. There are several investigations that typify the physicochemical qualities of European honeys depending on the botanical origin [1,4,15,16,18,19,20,22,25,26,27,28,29]. However, the reliability of the results depends on the large extent of the number of samples treated with similar characteristics, which allows the successful valorization of unifloral honey.
Currently, tools that handle a large set of data and resulting variables are being sought to identify common physicochemical patterns in the honey. Multivariate analysis has the advantage of extracting information from complex data sets using mathematical and statistical techniques. Among them, the principal component analysis (PCA), the cluster analysis (CA) or the linear discriminant analysis (LDA) offers the possibility of analyzing a complex food matrix such as honey and making a possible classification by botanical origin [4,18,20,21,22,27,28,29].
The aim of this study was to characterize and classify honeys collected in a specific Atlantic region according to their botanical origin. Firstly, a description of the palynological and physicochemical characteristics of the set of honey samples from the geographic area is detailed. Secondly, it is intended to classify honey regarding the botanical origin using a multivariate statistical treatment applied to pollen variables and physicochemical quality parameters such as color, total phenol, and flavonoid content. The similarities and the weight of the variables analyzed according to different statistical treatments are evaluated for the subsequent interpretation of the classification rate of the samples according to their botanical origin.

2. Materials and Methods

2.1. Geographical Origin of Honey Samples

The study was carried out with 362 samples of honey collected in different localities of Galicia (Northwest Spain). Sampling was carried out in collaboration with beekeepers and technicians of beekeeping associations as well as PGI Miel de Galicia during eight harvest seasons. The samples were deposited in glass containers and transported to the laboratory in the Faculty of Sciences (Ourense, Spain). The quality parameters were carried out upon arrival at the laboratory, and the remaining samples were stored frozen until further analytical determinations. The analytical determinations were performed in duplicate.
The honeys were collected without classification of their botanical origin. The palynological analysis and the physicochemical parameters (detailed below) were the basis for classifying the honey samples according to their botanical origin and subsequent statistical treatment.

2.2. Melissopalynological Analysis of Honey

The melissopalynological analysis of honey was performed by extracting the sediment from the samples. A total of ten grams of sample were dissolved in 50 mL of water until it had been completely dissolved. After this, the solution was centrifuged at 4500 rpm (3373 g) twice, for 10 min and 5 min, respectively [5]. Then, the supernatant was discarded. For the quantitative analysis, the volume of the sediment was completed until 5 mL, after stirring 10 μL of the sediment, was deposited in a slide over a heat plate. When the aliquot dried up, a drop of glycerol–gelatine stained with fuchsine was added to the sample and then covered with the coverslip for counting. In the case of qualitative analysis, new centrifugation at the same conditions was performed, and the supernatant was discarded. An aliquot of 100 μL of the sediment was deposited on the slide and covered as commented before. The procedure was carried out by duplicate. The identification and counting of the pollen grains were performed in both subsamples under optical microscopy (Olympus BX50 microscope, Olympus Corp., Tokyo, Japan). The minimum number of pollen grains counted for qualitative analysis was 700 pollen grains, and the results of each pollen were expressed as percentage over the total number of pollen grains counted. For quantitative analysis, all the pollen grains in both subsamples were counted, and results were expressed as number of pollen grains by gram of honey.

2.3. Determination of Quality Parameters: Moisture, pH and Electrical Conductivity

The moisture content of the honey was determined with an ABBE URA-2WAJ-325 digital refractometer (Auxilab S.L., Navarra, Spain) using the refractive index values from the Chataway table at 20 °C. The pH was measured directly on honey sample solutions in bi-distilled water (0.2 g/mL) using a pH meter (Crison micro pH 2001; Crison Instruments S.A., Barcelona, Spain). The electrical conductivity (EC) was measured on the same honey solution with a portable conductivity meter (Knick Portamess 913 Conductivity, Beuckestr, Berlin), expressing the results as mS/cm. Schade method was used to determine diastase activity of honeys [30]. It was calculated based on the hydrolysis rate of the starch solution by α-amylase present in a honey buffer solution at 40 °C. The amount of converted starch in the honey solution was analyzed using a UV-VIs spectrophotometer (Jenway 6305 UV-Visible spectrophotometer, Staffordshire, UK) at an absorbance of 660 nm. Measurements were taken at various time intervals until an absorbance of less than 0.235 was reached. Finally, the diastase activity was calculated as the diastase number or grams of hydrolyzed starch per hour per 100 g of honey.

2.4. Determination of Color

The color of the honey was determined with a HANNA Honey Color C221 colorimeter (HANNA C221 Honey Color Analyzer, Rhode Island, RI, USA). The fluidity of the honeys was previously treated to the measurement of the color for the correct reading. The treatment consisted of heating the slightly fluid or crystallized honeys no more than 45 °C in a thermostatic bath [13]. After a short break (for the total elimination of possible bubbles), approximately 4 mL of honey sample was introduced into a smooth plastic vial. Glycerin was used to calibrate the instrument. Finally, the color was expressed in mm, according to Pfund scale.

2.5. Determination of Total Phenol and Flavonoid Concentration

The total phenol content (TPC) and the total flavonoid content (TFC) were determined by spectrophotometric techniques according to reference methods proposed by Singleton et al. [31] and Arvouet-Grand et al. [32], respectively. For the determination of TPC, solutions of honey samples (0.1 g/mL) were prepared. These solutions were mixed with Folin–Ciocalteu reagent and calcium carbonate solutions, and the absorbance at 765 nm using a UV-Vis spectrophotometer (Jenway 6305, Staffordshire, UK) was measured. A calibration curve was obtained using gallic acid solutions (0.01–0.50 mg/mL) as a reference standard to quantify TPC. The methodology for the determination of TFC starts with preparing the honey sample solutions (0.33 g/mL). Then, a volume of aluminum chloride solution was added, and the absorbance was measured against a blank at 425 nm with a UV-Vis spectrophotometer (Jenway 6305, Staffordshire, UK). For the quantification of flavonoids, a curve with quercetin (0.002–0.01 mg/mL) as a reference standard was used. Finally, TPC and TFC were expressed in mg/100 g honey.

2.6. Statistical Analysis

The statistical treatments were carried out with SPSS Statistic 23.0 (IBM SPSS Statistics, Armonk, NY, USA) and Statgraphics Centurion 17.0 for Windows (Statgraphics Technologies, Inc., The Plains, VA, USA). Multivariate techniques were applied as association tools, searching for common patterns and relationships in masses of data. Principal component analysis (PCA) and cluster analysis (CA) were performed to reduce the amplitude of the data matrix and to establish significant relationships between palynological and physicochemical variables of honey. The statistical results were represented graphically with precise representations that integrate the interrelation of the significant elements between the main pollen types and the physicochemical parameters. Based on the pollen profile and physicochemical parameters, the honeys into four unifloral honey groups were classified: chestnut, blackberry, heather, and eucalyptus (with 52, 56, 36 and 33 samples, respectively), a honeydew honey group (with 53 samples) and the multifloral group (with 132 samples).
The classification of honey samples was checked through linear discriminant analysis (LDA) based on certain similarities in their physicochemical and botanical characteristics. Therefore, LDA was tested to quantify the probability of belonging to one type of honey or another. Finally, with the objective of comparing the groups of samples classified according to the results of the multivariate treatment, an analysis of variance (ANOVA) was performed using the Bonferroni test (p < 0.05).

3. Results

3.1. Representation of Botanical Diversity in Galician Honeys

The variability of families and pollen types identified in all the samples are classified as 52 and 111, respectively. The pollen types with the highest representation in the honey samples were Rubus, Castanea sativa, Cytisus type, Erica, Eucalyptus, Trifolium type, Quercus, and Echium (in more than 60% of the honeys) (Table 1).
The pollen grains with a higher mean percentage were Castanea sativa (42.9%), Rubus (23.2%), Eucalyptus (13.4%), Erica (7.9%), and Cytisus type (5.3%). Eucalyptus, Castanea sativa, and Rubus had a maximum value above 90%. Finally, the diversity in the pollen profile of the honey produced in Galicia was reflected in the quantitative analysis performed in the sediment of the samples. The counted number of pollen grains had an average value of 20,879 grains/g, with a wide standard deviation of 21,398 grains/g.

3.2. Physicochemical Characteristics of Honeys

Descriptive analyses for the results of physicochemical parameters and color are shown in Table 2. The moisture content of the studied honeys showed a mean value of 17.74%, with a range between 14.4% and 21.2%. The minimum values of pH and EC were 3.29 and 0.22 ms/cm, and the maximum values were 5.14 and 1.65 ms/cm, respectively. In terms of mean values, 4.24 for pH and 0.76 mS/cm for EC were obtained. Diastase content had a low mean value (21.31), with a range between 6.14 and 44.04. The mean values for TPC and TFC were 116.43 mg/100 g and 6.72 mg/100 g, respectively. However, these compounds presented a great variability in the set of samples, with a range between 33.91 and 254.5 mg/100 g for TPC and 1.28 and 16.7 mg/100 g for TFC. Finally, the honeys ranged from an amber color to dark color, from 150 mm to 36 mm on the Pfund scale.

3.3. Distribution of Honeys According to Botanical Origin, Physicochemical Parameters and Multivariate Classification Techniques

PCA was the multivariate technique used to simplify the large database matrix, and at the same time, it allowed us to show the relationships between the physicochemical and pollen variables. A total of four components were extracted that explain 76.14% of the variability of the original data (Table 3). The first two components accounted for more than 50% of the data variability. The variables with the greatest weight in the first component of PCA were pollen variables, such as Eucalyptus and Castanea sativa, and the physicochemical variables EC, TPC, and TFC. In the second component, Erica, Rubus, moisture content, color, and pH had the highest weight in the analysis.
The graphical representation of these two components shows the distribution of the variables included in the analysis and the distribution of the honey samples according to them (Figure 1). The honeys were categorized by a number from 1 to 6 due to their botanical origin (1: chestnut, 2: blackberry, 3: eucalyptus, 4: heather, 5: honeydew, 6: multifloral) for better visualization of the distribution. Thus, the honeys with a high percentage of Castanea sativa, high EC, pH, higher enzyme content and higher TFC were located on the left side of the quadrant. TPC variable was located between the honeys with a high content of Castanea sativa and Erica. At the same time, the samples with higher Erica pollen were those with higher moisture content. The Rubus pollen variable is placed on the opposite side of Erica. Finally, a group of samples was located together with the Eucalyptus variable on the opposite side of Castanea sativa pollen, EC, pH, color, and diastase content. Therefore, the honeys located on the right side of the quadrant are characterized by being the clearest, with lower EC and pH, lower diastase content, and lower TPC and TFC, mainly honeys which predominate Eucalyptus pollen. The darker honeys had the highest electrical conductivity, flavonoid content, and frequently large percentage of Castanea pollen.
CA was applied to analyze the groups of homogeneous honeys based on the physicochemical variables and main pollens (Figure 2). The results showed a good grouping of honeys of heather and eucalyptus (clusters E and D, respectively). Cluster C grouped blackberry samples with some honeydew and multifloral honeys. Cluster B included chestnut honeys and some honeydew honeys. Finally, cluster A (the largest group of samples, with 102 honeys) included samples of chestnut, honeydew, blackberry and multifloral.

3.4. Classification Rate of Honeys According to Botanical Origin

Based on the results of the pollen profile and the physicochemical parameters, 177 honeys were classified as unifloral honeys (chestnut, blackberry, eucalyptus, and heather), 53 as honeydew honeys and the remaining samples were grouped as multifloral honey. This previous classification was considered in the LDA for the discrimination of honeys based on botanical origin. LDA results were satisfactory, with five discriminant functions (Table 4). The first three functions optimally separated the samples with a percentage of the relative variance of the data higher than 90% and a canonical correlation greater than 0.80 with eigenvalue values above 2.2. Results of LDA showed that two statistically significant discriminant functions are formed (Wilks Lambda = 0.01, Chi-Square = 1777.24, degrees of freedom = 55, p < 0.05 for the first function, and Wilks Lambda = 0.03, Chi-Square = 1211.01, degrees of freedom = 40, p < 0.05 for the second, respectively). Low values of Wilks Lambda (close to 0) indicated high discriminant power because the mean of the explanatory variables included in the analysis is different between the groups (honey type), mainly in the first two discriminant functions.
The classification of honeys by LDA resulted in satisfactory, correctly classifying 89.8% of all the samples. Blackberry, eucalyptus, and heather honeys were properly classified (100%). However, five samples of chestnut honeys (9.6%), 14 samples of honeydew honeys (22.4%) and 18 samples of multifloral honeys (13.6%) were misclassified (Table 5). Figure 3 shows the distribution of the samples by honey type, with a clear differentiation of heather and eucalyptus honeys with respect to the others, which is due to their different physicochemical qualities. The first two functions extracted by LDA are dominated by the variables: Eucalyptus (with a standardized discriminant function coefficient of 1.07), Erica (−0.88), Rubus (0.37), color (0.32), and TPC (0.31).

3.5. Physicochemical and Melissopalynological Characterization of Unifloral and Honeydew Honeys

The characterization of the unifloral, honeydew and multifloral honeys according to the main palynological and physicochemical characteristics through the multivariate treatment are detailed in Table 6. An analysis of variance (ANOVA) according to the Bonferroni test was applied to evaluate the significant differences and similarities by honey type. The most outstanding characteristics of honey type are detailed below.
Based on these results, chestnut honeys are characterized by a mean percentage of Castanea sativa of 76.8% (significantly different from the other characterized honey types, p < 0.05) and with a confidence limit of 95% above 70%. The counted pollen grains were significantly higher in average terms than eucalyptus, heather and multifloral honeys. EC and pH of chestnut honeys had significantly higher mean values (1.02 mS/cm and 4.5, respectively) with respect to unifloral honeys (p < 0.05) but similar to honeydew honey (1.14 mS/cm and 4.5, respectively). The diastase content was intermediate, with a mean value significantly lower (23.5) than honeydew honey (29.4) and significantly higher than eucalyptus honey (14.6) (p < 0.05). Regarding the quantification of TPC and TFC, the chestnut honeys were characterized by a significantly higher content with respect to the other unifloral honeys (p < 0.05), except for heather honey (with similar mean values) and honeydew honey that have higher mean values (p < 0.05). Finally, the color of this group of honeys is dark, with a mean value of 128 mm, similar to heather and honeydew honeys (117 and 142 mm, respectively).
The blackberry honeys were characterized by a mean percentage of Rubus of 56.7% and a lower limit (95% confidence level) above 53%. Some physicochemical characteristics that differentiate them from other honeys of the area were pH, which is significantly higher (4.3, p < 0.05) and TFC, with a significantly lower mean value (6.1 mg/100 g, p < 0.05). In addition, TPC was significantly lower than chestnut, heather and honeydew honeys (p < 0.05). On the other hand, the EC of blackberry honey (0.69 mS/cm) was significantly lower with respect to chestnut and honeydew honeys (mean values above 1.0 mS/cm) and significantly higher with eucalyptus honey (0.51 mS/cm).
The eucalyptus honeys were characterized by a mean value of 73.8% for Eucalyptus pollen, with a minimum value at a 95% confidence level of 70.1%. These are honeys with the lightest color (77 mm), lower pH, EC, and TFC, and are significantly different to the other groups of honeys studied (p < 0.05). pH was similar to heather honeys (4.1) but significantly lower than the other unifloral honeys (p < 0.05).
The mean percentage of Erica in the heather honeys was 36.7%, with confidence limits between 33.1% and 40.3%. The counted pollen grains were significantly lower than chestnut and blackberry honeys (p < 0.05). Despite having a lower percentage of dominant pollen representation than other unifloral honeys, it has well-marked physicochemical properties. It highlights the color, with a significantly higher average value (117 mm) with respect to eucalyptus and blackberry and significantly lower with honeydew honeys (p < 0.05). TPC and TFC presented significantly higher mean values (143.1 mg/100 g and 7.3 mg/100 g, respectively), similar to chestnut honeys (p < 0.05).
The honeydew honey had more pollen diversity, with mean values significantly different in Rubus and Castanea sativa compared to the other types of unifloral honey (p < 0.05), but with lower values. This type of honey presented physicochemical parameters statistically different to eucalyptus, heather, and blackberry honey for pH, EC, and color (with significantly higher values, p < 0.05). However, the values were similar to those of chestnut honey, hence the possibility of finding honeys based on nectar secretion obtained from chestnut and contributions of honeydew that give particular properties to these samples.
Finally, multifloral honeys were composed of samples with a heterogeneous pollen profile. In most of these honeys, there was not a predominant pollen type in the pollen spectra, or if there were one, it was generally Castanea sativa. Rubus, Eucalyptus or Cytisus type usually appeared with values over 10%. Regarding the physicochemical characteristics, they differed in terms of botanical contributions. Although generally, they had lower pH, EC and moisture content and an intermediate enzymatic and polyphenolic content (TPC and TFC). Regarding the color scale, it ranged from dark amber to light amber, with a mean value slightly greater than eucalyptus honey.

4. Discussion

Honey, as a valuable natural product of the bee, offers substantial nutritional, therapeutic and medicinal benefits attributed to its botanical origin resulting from its complex chemical composition [2,29]. Some analytical methods are standardized for honey, but an exhaustive analysis of a set of analytes and physical properties is required for its correct characterization [9,23]. Experts from the scientific community and the beekeeping sector emphasized that the labelling with respect to certain botanical or geographical origins cannot be conducted based on a single group of chemical markers but rather on a combination of several [9]. Hence, the characterization of honeys based on their physicochemical characteristics supported by a palynological analysis, sensory analysis, and biologically active compounds contributes to supporting the demand of consumers, regulatory councils, and the beekeeping sector.
Considering the honey types studied from the Northwest of Spain, chestnut, heather, and honeydew honeys had the highest TPC and TFC in comparison to eucalyptus and blackberry honeys. With respect to the quality parameters, chestnut and honeydew honeys had the highest values of EC, pH, diastase, and color by the Pfund scale, unlike heather and eucalyptus honeys which had the lowest EC and color. In previous research on the composition of honeys collected in the Atlantic area of the Iberian Peninsula, the contribution of Erica pollen to the content of polyphenolic compounds was reported [1,26], and that derives in a high content of these compounds in unifloral heather honeys, as found in the present study. The physicochemical pattern of Galician chestnut honeys was common to chestnut honeys produced in other countries [25,33], although higher EC, TPC and lower color by CIELab were found for Croatian honeys [27], and lower TPC for Italian chestnut honeys [18]. Chestnut honeys from Tenerife Island are characterized by higher EC and pH than Galician honeys as a consequence of the specific edaphoclimatic characteristics of this territory and the abundance of endemic plants in the Islands [20]. Galician heather honeys were analogous to Portuguese heather honeys [26,33]. In fact, most are obtained from the same species standing out as Erica umbellata and E. arborea in soil conditions comparable to those of the Galician community. However, Algerian heather honeys showed slightly higher EC, TFC and color measured by the Pfund scale [11,34]. In the case of eucalyptus honey, the predominant specie for Galician honey production is Eucalyptus globulus, which flowers in winter and the early spring, in contrast to other European areas where summer production obtained from species such as E. camaldulensis is relevant. Concretely, an important production of unifloral eucalyptus honey was documented in Italy, Portugal, and other areas of Spain [25]. The physicochemical properties were similar to these honeys produced in Portugal because of similarities in the ecosystems formed by eucalyptus trees [33]. However, Italian eucalyptus honeys had lower TPC [18] than Galician honeys. In the case of Uruguayan eucalyptus honeys, higher EC, pH, and color were determined [35], while for Algerian eucalyptus honey, similar physicochemical characteristics were reported [34].
Honeydew honey usually presents higher content of bioactive compounds such as phenolic acids, flavonoids, proteins, and amino acids compared to blossom honeys, as documented in honeydew honeys from other geographical origins [2,12,18,27,29,34]. The difficulties in the discrimination of honeydew honey have been referenced because there are various sources of honeydew depending on the plant and the secretion itself by sucking insects or by the plant [17,21,25]. As indicated by Vasić et al. [17], the description of honeydew honeys without specific botanical attribution could explain the variability of physicochemical results in this group of honeys. Currently, in Spain, the types of honeydew honey identified are from holm oak (Quercus ilex), green oak (Quercus pyrenaica), and other oak (Quercus sp.) [21,22,36]. Therefore, in the specific case of this type of honey, a combination of chemical, physical, organoleptic, and statistical data could contribute to its discrimination.
The influence of physical and chemical indicators and the botanical origin in the color of honey is known, contributing to the diversity of its commercial assortment. According to Szabó et al. [37], the color of the honeys is attributed to the predominant plant species in them. However, it is a parameter strongly affected by biogeography because it determines the different plant communities. This justifies that unifloral honeys from the same plant origin present differences in physicochemical properties, as already mentioned for unifloral chestnut honeys produced in Canary Island [20], Croatia [27], or Italia [18]. Each biogeographical area has particular soil conditions, weather patterns, and plant communities that contribute to the differentiation of small nuances in honey. Hence the importance of geographical indications corroborated with the pollen profile in samples. Some specific chemical compounds may be involved in color attributes. This is the case with some phytochemicals, such as flavonoids, which are considered substances with a major effect on chromatic parameters [12,14]. These compounds are present in nectar and honeydew, therefore, are transferred to the hive and become part of the final product as honey [5,10]. In this context, the polyphenol content of dark-colored honeys (such as heather, chestnut or honeydew) [1,4,13,35,38] is usually greater than light samples (such as citrus, eucalyptus or acacia, for example) [1,4,18,23]. The role of polyphenols and flavonoids in plants is not only related to sensory properties such as color or odor or, for example, bitterness; their importance in chemical defense and oxidative processes is a crucial key [39]. Consequently, dark honeys show higher antioxidant activity [1,4,5,15,18,23,38].
At the same time, handling such a large matrix of data requires complex data management. The application of advanced multivariate chemometric techniques contributes to analyzing and extracting information from the dataset [23,40]. Therefore, chemometrics helps to reduce the complexity of large chemical data sets, offering better understandings, simplifications, explanations, and accuracy in results.
In recent years, statistical techniques combined with traditional analytical techniques, as well as melissopalynology, have been proven successful in making decisions about differentiation criteria in the complex matrix of honey [1,5,15,21,22,40]. Attempts at assessing botanical or geographic origins are made based on the physicochemical and antioxidant properties of honeys or their chemical composition with the use of multivariate techniques. Especially, PCA, CA and LDA are the best-known and most used techniques in the classification of honeys with satisfactory results [5,18,20,21,27,28,29,38]. Some Spanish honeys from Canary Island (fennel, chestnut, retama, and tajinaste honeys) were correctly classified (95.1%) by PCA-LDA, choosing the physicochemical variables (EC, acidity, pH, color, proline, diastase, invertase, fructose, glucose, trehalose, and melezitose) [20]. The quality parameters such as moisture, EC and pH were the variables that better discriminated Uruguayan honeys, with more than 80% of the samples correctly classified according to their floral origin (pasture, Citrus and Baccharis) by PCA-LDA [35]. The research carried out by Tarapatskyy et al. [29] showed that the specific content of phenolic acids, minerals, proline, and sugar, in combination with chemometric analysis (PCA, CA, and LDA), can successfully differentiate Polish honey samples according to their botanical origin (lime, buckwheat, and pine honeydew), as a preliminary verification of samples before performing pollen analysis. The application of an LDA model succeeded in classifying the Italian unifloral honeys (acacia, orange, honeydew, chestnut, strawberry tree, sulla, eucalyptus, dandelion, and linden), as they greatly differed in the polyphenol content and color attributes due to their botanical origin [18]. Eight physicochemical parameters (L*, a*, total dissolved solids, salinity, moisture, free acidity, total acidity, and dissolved solids/total acidity ratio) were enough for classifying Egyptian honeys such as clover or citrus honey. LDA function is classified correctly at a rate higher than 90% [41]. Akbari et al. [23] provided a 97% classification rate with PCA-DA for Iranian honeys from thyme, jujube, coriander, barberry, acacia, and alfalfa. PCA and LDA identified as significant variables main pollen data and physicochemical variables (free acidity, reducing sugars, and moisture) to build a discriminatory model with a cumulative variance of 90%, and correctly classifying three different groups of Argentine honeys (Eucalyptus, Salix humboldtiana, and Baccharis) [40]. Other mathematical treatments covering PCA and LDA models included a sugars profile for the discrimination of honeys from Southern Italy (Calabria region) by the botanical origin [28]. Some antioxidant properties, minerals, and color had the highest discriminating power for cherry, apple, saffron, and wild bush honeys collected in India, with a successful classification of 100% by LDA-CA [38]. In the case of the samples from Northwest Spain, the variables that better contributed to the discrimination between chestnut and honeydew honeys with PCA-LDA (97.6% of samples) were moisture, diastase, CIELab coordinates for color, flavonoids, radical scavenging activity, Mg, Na, fructose, turanose, maltose, trehalose, and main pollen variables [21]. The variables with the greatest discriminatory power using LDA for unifloral honeys of Citrus and Eucalyptus from southern Spain were water activity and EC [19]. Chemometrics provided satisfactory results for the classification of honey samples covered by PGI Miel de Galicia of the present study. All the blackberry, eucalyptus and heather honeys were properly grouped. In the case of chestnut and honeydew honeys, some samples were misclassified (5 and 14, respectively). It should be highlighted that chestnut honey and honeydew honey presented similar fingerprints regarding the routine physicochemical parameters, only a slightly high EC and color in honeydew honeys are marked, but the natural variations for each type of honeys complicate a full differentiation [21]. Other analyses, such as sensorial analysis, the polyphenol or volatile profile, could contribute to the discrimination. Choosing the most suitable chemical compounds for the discrimination of honey samples based on botanical origin is a complicated task and requires long professional experience in analytical techniques and knowledge of the composition and origin of this complex matrix. Some studies demonstrated the strong relationship between the floral origin of honey with physicochemical parameters, the profile and quantity of bioactive compounds (polyphenols, minerals) and color attributes. Therefore, these variables can be used as a simple approach for the discrimination of floral origin and, at the same time, as a preliminary evaluation of the antioxidant properties of honey.

5. Conclusions

Chemometrics contributed to the discrimination of the botanical origin of honeys produced in Galicia (NW Spain). The statistical classification methods on a large number of honey samples from different floral origins were tested, including unifloral honeys collected in PGI Miel de Galicia. The variables with the greatest weight in the differentiation of honeys based on PCA were EC, TPC, TFC, and dominant pollens (Eucalyptus, Erica, Rubus, and Castanea sativa). LDA classified the honeys in the function of the botanical origin, simplifying the interpretation of data from samples of the same geographical area. LDA showed better results graphically in the differentiation of the honey groups than CA. However, the differentiation of honeydew honey with chestnut honey from Galicia can sometimes be difficult due to having similar qualities in pH, EC, enzyme content, TPC, and TFC. TPC and TFC were higher for dark-colored honey (chestnut, heather, and honeydew) compared to light-colored honey (blackberry and eucalyptus). Therefore, the application of multivariate techniques helps to characterize honeys according to their botanical origin linked to a differentiating quality and to the geographical territory.

Author Contributions

O.E., M.S.R.-F., M.M. and M.C.S. conceived and designed the methodology. O.E., M.S.R.-F. and M.C.S. analyzed data. O.E., M.S.R.-F., M.M. and M.C.S. writing original draft preparation. O.E. and M.C.S. writing review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Xunta de Galicia (2014/2020-FEADER 2018/054B).

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors want to thank the collaboration of the beekeepers and the technician of the Regulation Council PGI Miel de Galicia for their support in the collection of samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Escuredo, O.; Míguez, M.; Fernández-González, M.; Seijo, M.C. Nutritional value and antioxidant activity of honeys produced in a European Atlantic area. Food Chem. 2013, 138, 851–856. [Google Scholar] [CrossRef]
  2. Seraglio, S.K.T.; Silva, B.; Bergamo, G.; Brugnerotto, P.; Gonzaga, L.V.; Fett, R.; Costa, A.C.O. An overview of physicochemical characteristics and health-promoting properties of honeydew honey. Food Res. Int. 2019, 119, 44–66. [Google Scholar] [CrossRef] [PubMed]
  3. Nikhat, S.; Fazil, M. History, phytochemistry, experimental pharmacology and clinical uses of honey: A comprehensive review with special reference to Unani medicine. J. Ethnopharmacol. 2022, 282, 114614. [Google Scholar] [CrossRef] [PubMed]
  4. Machado, A.A.; Almeida-Muradian, L.B.D.; Sancho, M.T.; Pascual-Maté, A. Composition and properties of Apis mellifera honey: A review. J. Apic. Res. 2018, 57, 5–37. [Google Scholar] [CrossRef]
  5. Rodríguez-Flores, M.S.; Escuredo, O.; Seijo, M.C. Characterization and antioxidant capacity of sweet chestnut honey produced in North-West Spain. J. Apic. Sci. 2016, 60, 19. [Google Scholar] [CrossRef] [Green Version]
  6. Becerril-Sánchez, A.L.; Quintero-Salazar, B.; Dublán-García, O.; Escalona-Buendía, H.B. Phenolic compounds in honey and their relationship with antioxidant activity, botanical origin, and color. Antioxidants 2021, 10, 1700. [Google Scholar] [CrossRef]
  7. Valverde, S.; Ares, A.M.; Elmore, J.S.; Bernal, J. Recent trends in the analysis of honey constituents. Food Chem. 2022, 387, 132920. [Google Scholar] [CrossRef] [PubMed]
  8. Isidorow, W.; Witkowski, S.; Iwaniuk, P.; Zambrzycka, M.; Swiecicka, I. Royal jelly aliphatic acids contribute to antimicrobial activity of honey. J. Apic. Sci. 2018, 62, 111–123. [Google Scholar] [CrossRef] [Green Version]
  9. Trifković, J.; Andrić, F.; Ristovojević, P.; Guzelmeric, E.; Yesilada, E. Analytical methods in tracing honey authenticity. J. AOAC Int. 2017, 100, 827–839. [Google Scholar] [CrossRef] [PubMed]
  10. Da Silva, P.M.; Gauche, C.; Gonzaga, L.V.; Costa, A.C.O. Honey: Chemical composition, stability and authenticity. Food Chem. 2016, 196, 309–323. [Google Scholar] [CrossRef]
  11. Ghorab, A.; Rodríguez-Flores, M.S.; Nakib, R.; Escuredo, O.; Haderbache, L.; Bekdouche, F.; Seijo, M.C. Sensorial, Melissopalynological and Physico-Chemical Characteristics of Honey from Babors Kabylia’s Region (Algeria). Foods 2021, 10, 225. [Google Scholar] [CrossRef] [PubMed]
  12. Vela, L.; De Lorenzo, C.; Pérez, R.A. Antioxidant capacity of Spanish honeys and its correlation with polyphenol content and other physicochemical properties. J. Sci. Food Agric. 2007, 87, 1069–1075. [Google Scholar] [CrossRef]
  13. Escuredo, O.; Rodríguez-Flores, M.S.; Rojo-Martínez, S.; Seijo, M.C. Contribution to the chromatic characterization of unifloral honeys from Galicia (NW Spain). Foods 2019, 8, 233. [Google Scholar] [CrossRef] [Green Version]
  14. Combarros-Fuertes, P.; Valencia-Barrera, R.M.; Estevinho, L.M.; Dias, L.G.; Castro, J.M.; Tornadijo, M.E.; Fresno, J.M. Spanish honeys with quality brand: A multivariate approach to physicochemical parameters, microbiological quality, and floral origin. J. Apic. Res. 2019, 58, 92–103. [Google Scholar] [CrossRef] [Green Version]
  15. Dobre, I.; Escuredo, O.; Rodriguez-Flores, S.; Seijo, M.C. Evaluation of several Romanian honeys based on their palynological and biochemical profiles. Int. J. Food Prop. 2014, 17, 1850–1860. [Google Scholar] [CrossRef]
  16. Karabagias, I.K.; Badeka, A.V.; Kontakos, S.; Karabournioti, S.; Kontominas, M.G. Botanical discrimination of Greek unifloral honeys with physico-chemical and chemometric analyses. Food Chem. 2014, 165, 181–190. [Google Scholar] [CrossRef]
  17. Vasić, V.; Gašić, U.; Stanković, D.; Lušić, D.; Vukić-Lušić, D.; Milojković-Opsenica, D.; Tešić, Ž.; Trifković, J. Towards better quality criteria of European honeydew honey: Phenolic profile and antioxidant capacity. Food Chem. 2019, 274, 629–641. [Google Scholar] [CrossRef] [Green Version]
  18. Preti, R.; Tarola, A.M. Chemometric evaluation of the antioxidant properties and phenolic compounds in Italian honeys as markers of floral origin. Eur. Food Res. Technol. 2022, 248, 991–1002. [Google Scholar] [CrossRef]
  19. Serrano, S.; Villarejo, M.; Espejo, R.; Jodral, M. Chemical and physical parameters of Andalusian honey: Classification of Citrus and Eucalyptus honeys by discriminant analysis. Food Chem. 2004, 87, 619–625. [Google Scholar] [CrossRef]
  20. Bentabol-Manzanares, A.; Hernández-García, Z.; Rodríguez-Galdón, B.; Rodríguez-Rodríguez, E.M.; Díaz-Romero, C. Physicochemical characteristics and pollen spectrum of monofloral honeys from Tenerife, Spain. Food Chem. 2017, 228, 441–446. [Google Scholar] [CrossRef]
  21. Rodríguez-Flores, M.S.; Escuredo, O.; Míguez, M.; Seijo, M.C. Differentiation of oak honeydew and chestnut honeys from the same geographical origin using chemometric methods. Food Chem. 2019, 297, 124979. [Google Scholar] [CrossRef] [PubMed]
  22. Seijo, M.C.; Escuredo, O.; Rodríguez-Flores, M.S. Physicochemical properties and pollen profile of oak honeydew and evergreen oak honeydew honeys from Spain: A comparative study. Foods 2019, 8, 126. [Google Scholar] [CrossRef] [Green Version]
  23. Akbari, E.; Baigbabaei, A.; Shahidi, M. Determination of the floral origin of honey based on its phenolic profile and physicochemical properties coupled with chemometrics. Int. J. Food Prop. 2020, 23, 506–519. [Google Scholar] [CrossRef] [Green Version]
  24. MAGRAMA. El Sector Apícola Español en 2020: Principales Magnitudes e Indicadores Económicos. Informes del Sector Apícola Español. Available online: https://www.mapa.gob.es/es/ganaderia/temas/produccion-y-mercados-ganaderos/indicadoreseconomicossectorapicola2020.pdf (accessed on 15 December 2022).
  25. Persano-Oddo, L.; Piro, R. Main European unifloral honeys: Descriptive sheets. Apidologie 2004, 35, S38–S81. [Google Scholar] [CrossRef]
  26. Estevinho, L.M.; Feás, X.; Seijas, J.A.; Vázquez-Tato, M.P. Organic honey from Trás-Os-Montes region (Portugal): Chemical, palynological, microbiological and bioactive compounds characterization. Food Chem. Toxicol. 2012, 50, 258–264. [Google Scholar] [CrossRef] [PubMed]
  27. Flanjak, I.; Kenjerić, D.; Bubalo, D.; Primorac, L. Characterisation of selected Croatian honey types based on the combination of antioxidant capacity, quality parameters, and chemometrics. Eur. Food Res. Technol. 2016, 242, 467–475. [Google Scholar] [CrossRef]
  28. Carabetta, S.; Di Sanzo, R.; Campone, L.; Fuda, S.; Rastrelli, L.; Russo, M. High-performance anion exchange chromatography with Pulsed Amperometric Detection (HPAEC–PAD) and chemometrics for geographical and floral authentication of honeys from Southern Italy (Calabria region). Foods 2020, 9, 1625. [Google Scholar] [CrossRef]
  29. Tarapatskyy, M.; Sowa, P.; Zaguła, G.; Dżugan, M.; Puchalski, C. Assessment of the botanical origin of polish honeys based on physicochemical properties and bioactive components with chemometric analysis. Molecules 2021, 26, 4801. [Google Scholar] [CrossRef]
  30. Bogdanov, S.; Martin, P.; Lullmann, C. Harmonized methods of the international honey commission. Apidologie 1997, 1–59. Available online: https://www.ihc-platform.net/ihcmethods2009.pdf (accessed on 20 July 2019).
  31. Singleton, V.L.; Orthofer, R.; Lamuela-Raventos, R.M. Analysis of total phenols and other oxidation substrates and antioxidants by means of Folin-Ciocalteu reagent. Meth. Enzymol. 1999, 299, 152–178. [Google Scholar]
  32. Arvouet-Grand, A.; Vennat, B.; Pourrat, A.; Legret, P. Standardization of propolis extract and identification of principal constituents. J. Pharm. Belg. 1994, 49, 462–468. [Google Scholar] [PubMed]
  33. Karabagias, I.K.; Maia, M.; Karabagias, V.K.; Gatzias, I.; Badeka, A.V. Characterization of eucalyptus, chestnut and heather honeys from Portugal using multi-parameter analysis and chemo-calculus. Foods 2018, 7, 194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Homrani, M.; Escuredo, O.; Rodríguez-Flores, M.S.; Fatiha, D.; Mohammed, B.; Homrani, A.; Seijo, M.C. Botanical origin, pollen profile, and physicochemical properties of Algerian honey from different bioclimatic areas. Foods 2020, 9, 938. [Google Scholar] [CrossRef] [PubMed]
  35. Corbella, E.; Cozzolino, D. Classification of the Floral Origin of Uruguayan Honeys by Chemical and Physical Characteristics Combined with Chemometrics. LWT-Food Sci. Technol. 2006, 39, 534–539. [Google Scholar] [CrossRef]
  36. Pita-Calvo, C.; Vázquez, M. Honeydew honeys: A Review on the characterization and authentication of botanical and geographical origins. J. Agric. Food Chem. 2018, 66, 2523–2537. [Google Scholar] [CrossRef]
  37. Szabó, R.T.; Mézes, M.; Szalai, T.; Zajácz, E.; Kovács-Weber, M. Colour identification of honey and methodical development of its instrumental measuring. Columella J. Agric. Environ. Sci. 2016, 3, 29–36. [Google Scholar] [CrossRef]
  38. Nayik, N.; Nanda, V. A chemometric approach to evaluate the phenolic compounds, antioxidant activity and mineral content of different unifloral honey types from Kashmir, India. LWT-Food Sci. Technol. 2016, 74, 504–513. [Google Scholar] [CrossRef]
  39. Singh, S.; Kaur, I.; Kariyat, R. The multifunctional roles of polyphenols in plant-herbivore interactions. Int. J. Mol. Sci. 2021, 22, 1442. [Google Scholar] [CrossRef]
  40. Rios, F.; Sanchez, A.C.; Lobo, M.; Lupo, L.; Coelho, I.; Castanheira, I.; Samman, N. A Chemometric Approach: Characterization of Quality and Authenticity of Artisanal Honeys from Argentina. J. Chemom. 2014, 28, 834–843. [Google Scholar] [CrossRef]
  41. Karabagias, I.K.; Karabournioti, S. Discrimination of clover and citrus honeys from Egypt according to floral type using easily assessable physicochemical parameters and discriminant analysis: An external validation of the chemometric approach. Foods 2018, 7, 70. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Biplot of the first two components and distribution of honey samples extracted by PCA. EC: electrical conductivity; TPC: total phenol content; TFC: total flavonoid content. Numerical nomenclature of honeys according to botanical origin: 1, chestnut; 2, blackberry; 3, eucalyptus; 4, heather; 5, honeydew; 6, multifloral.
Figure 1. Biplot of the first two components and distribution of honey samples extracted by PCA. EC: electrical conductivity; TPC: total phenol content; TFC: total flavonoid content. Numerical nomenclature of honeys according to botanical origin: 1, chestnut; 2, blackberry; 3, eucalyptus; 4, heather; 5, honeydew; 6, multifloral.
Foods 12 01493 g001
Figure 2. Honey samples distributed in five groups (cluster A, B, C, D and E) according to CA. Numerical nomenclature of honeys according to botanical origin: 1, chestnut; 2, blackberry; 3, eucalyptus; 4, heather; 5, honeydew; 6, multifloral.
Figure 2. Honey samples distributed in five groups (cluster A, B, C, D and E) according to CA. Numerical nomenclature of honeys according to botanical origin: 1, chestnut; 2, blackberry; 3, eucalyptus; 4, heather; 5, honeydew; 6, multifloral.
Foods 12 01493 g002
Figure 3. Plot of the representation of the first two discriminant functions and distribution of samples by honey type.
Figure 3. Plot of the representation of the first two discriminant functions and distribution of samples by honey type.
Foods 12 01493 g003
Table 1. Descriptive analysis of the families and pollen types found in more than 30% of honey samples (Rep.). Distribution of the samples according to the frequency classes (D: dominant pollen; A: accompanying pollen; I: important pollen; R: frequent pollen; P: minor pollen). Max.: maximum.
Table 1. Descriptive analysis of the families and pollen types found in more than 30% of honey samples (Rep.). Distribution of the samples according to the frequency classes (D: dominant pollen; A: accompanying pollen; I: important pollen; R: frequent pollen; P: minor pollen). Max.: maximum.
FamilyPollen TypeRep.(%)MeanSDMax.D > 45A (15–45)I (3–15)R (1–3)P (0–1)
RosaceaeRubus99.223.219.891.318.037.630.47.75.5
FagaceaeCastanea sativa98.642.924.892.545.638.19.73.61.7
FabaceaeCytisus type97.85.36.648.60.36.645.929.016.0
EricaceaeErica96.77.911.668.61.914.935.622.421.8
MyrtaceaeEucalyptus81.513.423.894.811.910.517.116.925.1
FabaceaeTrifolium type73.80.81.515.4-0.36.416.950.3
FagaceaeQuercus71.80.71.927.8-0.34.714.951.9
BoraginaceaeEchium67.10.91.813.1--9.413.843.9
SalicaceaeSalix59.70.72.327.4-0.84.110.544.2
PlantaginaceaePlantago53.60.20.67.4--0.64.448.6
PoaceaePoaceae53.30.20.53.3--0.85.846.7
RosaceaeCrataegus monogyna type44.50.31.117.5-0.31.96.435.9
RosaceaePrunus type42.00.20.43.9--0.82.538.7
ApiaceaeConium maculatum type41.40.30.77.0--1.45.234.8
BrassicaceaeBrassica type40.10.20.55.5--0.63.336.2
RhamnaceaeFrangula alnus35.90.41.821.8-0.61.46.427.6
CampanulaceaeCampanula type35.10.20.43.9--0.63.331.2
ResedaceaeSesamoides34.30.20.916.0-0.3-3.330.7
ScrophulariaceaeScrophularia type32.30.10.22.9---1.131.2
Table 2. Descriptive analysis of the physicochemical parameters determined in the set of honey samples. SD: standard deviation, EC: electrical conductivity; TPC: total phenol content; TFC: total flavonoid content.
Table 2. Descriptive analysis of the physicochemical parameters determined in the set of honey samples. SD: standard deviation, EC: electrical conductivity; TPC: total phenol content; TFC: total flavonoid content.
MeanSDMinimumMaximum
Moisture (%)17.741.0214.4021.20
pH4.240.353.295.15
EC (mS/cm)0.760.280.221.65
Diastase content21.318.796.1444.04
TPC (mg/100 g)116.4344.3933.91254.50
TFC (mg/100 g)6.722.721.2816.70
Color (mm)1053036150
Table 3. Results of components extracted and weights of pollen and physicochemical variables by PCA. EC: electrical conductivity; TPC: total phenol content; TFC: total flavonoid content.
Table 3. Results of components extracted and weights of pollen and physicochemical variables by PCA. EC: electrical conductivity; TPC: total phenol content; TFC: total flavonoid content.
ComponentsC1C2C3C4
Eigenvalue4.031.881.291.18
Variance (%)36.6317.0811.7310.69
Cumulative variance (%)36.6353.7165.4576.14
Component weights of pollen variables
Rubus−0.02−0.380.720.07
Castanea sativa−0.30−0.08−0.18−0.54
Erica−0.010.600.230.11
Eucalyptus0.310.06−0.540.34
Component weights of physicochemical variables
Moisture−0.030.410.06−0.55
pH−0.29−0.34−0.150.21
EC−0.44−0.10−0.19−0.03
Diastase content−0.31−0.15−0.19−0.20
TPC−0.430.090.020.28
TFC−0.420.170.080.22
Color−0.290.370.030.27
Table 4. Results of the statistics according to LDA. DF: degrees of freedom.
Table 4. Results of the statistics according to LDA. DF: degrees of freedom.
Discriminant FunctionEigenValueRelative PercentageCanonical CorrelationWilks LambdaChi-SquareDFp
13.9837.270.890.011777.2455<0.001
23.4832.540.880.031211.0140<0.001
32.2921.420.830.14682.527<0.001
40.76.520.640.47262.6516<0.001
50.242.260.440.8176.317<0.001
Table 5. Number of samples and correct percentage of classification (in brackets) by honey type according to the LDA.
Table 5. Number of samples and correct percentage of classification (in brackets) by honey type according to the LDA.
Predicted Honey Type (%)
Honey TypenChestnutBlackberryEucalyptusHeatherHoneydewMultifloral
Chestnut5247 (90.4)0--2 (3.8)3 (5.8)
Blackberry56056 (100)--00
Eucalyptus330033 (100)000
Heather3600036 (100)00
Honeydew537 (13.2)5 (9.4)0039 (73.6)2 (3.8)
Multifloral1323 (2.3)4 (3.0)6 (4.5)3 (2.3)2 (1.5)114 (86.4)
Table 6. Descriptive analysis of the botanical and physicochemical characteristics by honey type based on LDA. * The numbers show significant differences in the means for each honey type (1: chestnut, 2: blackberry, 3: eucalyptus, 4: heather, 5: honeydew and 6: multifloral) according to the Bonferroni test (p < 0.05). SD: standard deviation, PG/g: number of pollen grains per gram of honey; EC: electrical conductivity; TPC: total phenol content; TFC: total flavonoid content.
Table 6. Descriptive analysis of the botanical and physicochemical characteristics by honey type based on LDA. * The numbers show significant differences in the means for each honey type (1: chestnut, 2: blackberry, 3: eucalyptus, 4: heather, 5: honeydew and 6: multifloral) according to the Bonferroni test (p < 0.05). SD: standard deviation, PG/g: number of pollen grains per gram of honey; EC: electrical conductivity; TPC: total phenol content; TFC: total flavonoid content.
MeanSDLower Limit than 95%Upper Limit than 95%ANOVA *
Chestnut honey (n = 57)
Rubus (%)12.66.710.814.42, 3, 5, 6
Castanea sativa (%)76.88.074.678.92, 3, 4, 5, 6
Erica (%)4.03.73.05.04
Eucalyptus (%)1.32.50.61.93, 6
PG/g30,23923,50523,88536,5933, 4, 6
Moisture (%)18.21.017.918.52, 3, 5, 6
pH4.50.44.44.62, 3, 4, 6
EC (mS/cm)1.020.210.971.082, 3, 4, 6
Diastase content23.57.221.625.43, 5
TPC (mg/100 g)122.829.4115.0130.62, 3, 5
TFC (mg/100 g)8.22.27.68.82, 3, 5, 6
Color (mm Pfund)128241221352, 3, 6
Blackberry honey (n = 65)
Rubus (%)56.711.453.959.51, 3, 4, 5, 6
Castanea sativa (%)26.613.923.130.01, 3, 5, 6
Erica (%)3.34.52.14.44, 6
Eucalyptus (%)1.62.50.92.23, 6
PG/g26,57129,84718,72334,4184
Moisture (%)17.41.017.117.61, 4
pH4.30.34.24.41, 3, 4, 5, 6
EC (mS/cm)0.690.250.630.751, 3, 5
Diastase content19.77.018.021.53, 5
TPC (mg/100 g)94.731.786.8102.51, 4, 5
TFC (mg/100 g)6.12.15.66.71, 3, 4, 5, 6
Color (mm Pfund)9630891041, 3, 4, 5
Eucalyptus honey (n = 39)
Rubus (%)2.22.61.33.01, 2, 4, 5, 6
Castanea sativa (%)8.37.95.710.91, 2, 4, 5, 6
Erica (%)3.03.31.94.04, 6
Eucalyptus (%)73.811.470.177.51, 2, 4, 5, 6
PG/g16,37110,94612,12620,6151
Moisture (%)17.40.917.117.71, 4
pH4.10.34.04.21, 2, 5
EC (mS/cm)0.510.100.480.541, 2, 4, 5, 6
Diastase content14.67.012.416.91, 2, 4, 5, 6
TPC (mg/100 g)83.738.271.396.11, 4, 5, 6
TFC (mg/100 g)4.61.24.25.01, 2, 4, 5
Color (mm Pfund)771971831, 2, 4, 5, 6
Heather honey (n = 39)
Rubus (%)10.98.78.113.72, 3, 5, 6
Castanea sativa (%)28.414.723.633.11, 3, 5, 6
Erica (%)36.711.033.140.31, 2, 3, 5, 6
Eucalyptus (%)5.69.52.68.73, 6
PG/g11,04511,028742014,6701, 2
Moisture (%)18.51.218.118.92, 3, 5, 6
pH4.00.23.94.11, 2, 5
EC (mS/cm)0.680.170.620.731, 3, 5
Diastase content20.18.217.422.83, 5
TPC (mg/100 g)143.149.6127.0159.22, 3, 6
TFC (mg/100 g)7.32.06.78.02, 3, 5, 6
Color (mm Pfund)117211101242, 3, 5, 6
Honeydew honey (n = 43)
Rubus (%)26.912.723.030.81, 2, 3, 4, 6
Castanea sativa (%)53.619.147.759.41, 2, 3, 4
Erica (%)3.53.62.44.64
Eucalyptus (%)1.02.30.31.73, 6
PG/g18,53114,36513,35223,710
Moisture (%)17.40.917.117.61, 4
pH4.50.24.44.52, 3, 4, 6
EC (mS/cm)1.140.201.071.202, 3, 4, 6
Diastase content29.47.927.031.81, 2, 3, 4, 6
TPC (mg/100 g)166.239.2154.1178.31, 2, 3, 6
TFC (mg/100 g)11.22.210.511.81, 2, 3, 4, 6
Color (mm Pfund)142151371462, 3, 4, 6
Multifloral honey (n = 119)
Rubus (%)19.511.717.421.61, 2, 3, 4, 5
Castanea sativa (%)47.817.044.750.91, 2, 3, 4
Erica (%)6.26.15.17.32, 3, 4
Eucalyptus (%)12.814.810.115.51, 2, 3, 4, 5
PG/g18,15318,86814,39021,9161
Moisture (%)17.70.817.617.91, 4
pH4.10.34.14.21, 2, 5
EC (mS/cm)0.660.190.620.691, 3, 5
Diastase content20.89.119.122.43, 5
TPC (mg/100 g)109.337.9102.4116.23, 4, 5
TFC (mg/100 g)5.21.45.05.51, 2, 4, 5
Color (mm Pfund)912187941, 3, 4, 5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Escuredo, O.; Rodríguez-Flores, M.S.; Míguez, M.; Seijo, M.C. Multivariate Statistical Approach for the Discrimination of Honey Samples from Galicia (NW Spain) Using Physicochemical and Pollen Parameters. Foods 2023, 12, 1493. https://0-doi-org.brum.beds.ac.uk/10.3390/foods12071493

AMA Style

Escuredo O, Rodríguez-Flores MS, Míguez M, Seijo MC. Multivariate Statistical Approach for the Discrimination of Honey Samples from Galicia (NW Spain) Using Physicochemical and Pollen Parameters. Foods. 2023; 12(7):1493. https://0-doi-org.brum.beds.ac.uk/10.3390/foods12071493

Chicago/Turabian Style

Escuredo, Olga, María Shantal Rodríguez-Flores, Montserrat Míguez, and María Carmen Seijo. 2023. "Multivariate Statistical Approach for the Discrimination of Honey Samples from Galicia (NW Spain) Using Physicochemical and Pollen Parameters" Foods 12, no. 7: 1493. https://0-doi-org.brum.beds.ac.uk/10.3390/foods12071493

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop