Next Article in Journal
The Biological Effects of Forsythia Leaves Containing the Cyclic AMP Phosphodiesterase 4 Inhibitor Phillyrin
Next Article in Special Issue
Bee Bread Can Alleviate Lipid Abnormalities and Impaired Bone Morphology in Obese Zucker Diabetic Rats
Previous Article in Journal
Role of pH on Nanostructured SERS Active Substrates for Detection of Organic Dyes
Previous Article in Special Issue
Propolis, Bee Honey, and Their Components Protect against Coronavirus Disease 2019 (COVID-19): A Review of In Silico, In Vitro, and Clinical Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cluster Analysis Classification of Honey from Two Different Climatic Zones Based on Selected Physicochemical and of Microbiological Parameters

1
Department of Food Hygiene and Quality Management, Institute of Human Nutrition, Warsaw University of Life Sciences in Warsaw—SGGW, Nowoursynowska str. 166, 02-787 Warsaw, Poland
2
Apiculture Division, Institute of Animal Sciences, Warsaw University of Life Sciences in Warsaw—SGGW, 02-787 Warsaw, Poland
3
Apiculture Division in Puławy, The Research Institute of Horticulture, 24-100 Puławy, Poland
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(8), 2361; https://doi.org/10.3390/molecules26082361
Submission received: 23 March 2021 / Revised: 14 April 2021 / Accepted: 16 April 2021 / Published: 19 April 2021
(This article belongs to the Special Issue Bee Products: Recent Progress in Health Benefits Studies)

Abstract

:
The geographical origin of honey affects its composition, which is of key importance for the health-promoting properties and safety of the product. European regulations clearly define the physicochemical requirements for honey that determine the microbiological quality. On the other hand, legislation abolishes microbiological criteria. In the study 40 honey samples originating from two different climatic zones were analyzed. The water content, pH, water activity analysis and the microbiological quality of honey samples have been tested using the reference plate method (total viable count, yeast and molds, lactic acid bacteria, Bacillus spp.). The cluster classification showed that total viable count of bacteria could be used as a measure alternative to the count of Bacillus spp. and 70% of honeys from the tropical climate zone had different microbiological quality than honeys from the temperate climate zone but still under the level 3.0 log cfu/g. The study has revealed that geographical origin of honey may significantly affect the quality and safety of honey. It was considered that water content can be the most informative and handy marker of the microbiological quality of honeys. Analysis of lactic acid bacteria showed temperate climate zone honeys as a source of beneficial bacteria in the diet.

1. Introduction

The properties of honey result from their antimicrobial, antioxidant, enzymatic dietary and sensory properties, as well as prebiotic effects and the presence of probiotic bacteria [1,2]. Some strains of lactic acid bacteria have a beneficial effect on the human body by normalizing the microbiota of the gastrointestinal tract and the body’s antiallergic response [3,4]. The analysis of the presence of the 16S rRNA gene of Lactobacillus bacteria confirmed the presence of these bacteria in a large percentage of bee products, respectively: in honey—90.9%; in pollen—70.6%; in propolis—83.9%. Lactobacillus kunkeei was the dominant species (98%) of Lactobacillus spp. in bee products [5]. It was also proved that the Lactobacillus kunkeei YB38 strain promoted the production of IgA antibodies in humans. In vitro studies showed that strains YB83 and YB38 present in bee pollen increased the production of IgA antibodies in Peyer cells in mice and showed mitogenic activity. Depending on the species and strain of Lactobacillus, the immunomodulatory effects may vary, and these strains can safely improve the immune responses of the human body [5]. It has been shown that lactic acid bacteria have a beneficial effect on health with a concentration of at least 109 cfu/g. However, confirmation of the probiotic effect requires detailed research [6].
Among the major components of honey, the next can be reported: sugars, water, ni-trogenous substances, proteins, organic acids and polyphenols. Honey sugars are a complex mixture of carbohydrates dominated by fructose and glucose [7,8,9,10,11,12].
The composition of honey and the content of compounds influencing the properties of honey depends on many factors: the botanical origin of honey, environmental and climatic factors, as well as the honey extraction process [11,13,14].
The microbiota of honey also results on its chemical composition depending on the botanical, environmental and climatic factors. Beekeepers from the tropics, where the vegetation period is long, can harvest large quantities of honey almost all year round. However, high ambient temperature and humidity levels impede the aging, while the environment promotes the growth of bacteria and yeasts. Typically, excess water is removed from honey through heating, which leads to thermal stress, consequently, deterioration of health-promoting properties and succulence. On the other hand, in the moderate climate, the harvesting is conducted a few times, at sufficient intervals, over an approximately half-year period which favors the natural ripening process [9]. An element of the health quality of honey is microbiological safety. Ripe honeys are dominated by aerobic bacteria, while yeast and molds are in the minority. Identification of microbiome of unripe honey showed the presence of Gram-positive bacteria mainly Lactobacillus (1.0 × 102–1.2 × 103/g) and genera: Bacillus, Staphylococcus and Enterococcus. Gram-negative are represented by Gluconobacter bacteria (6.0 × 102–7.0 × 103/g), and intrinsic biota of bees which in 75% consists of Enterobacteriaceae (E. coli, Salmonella, Shigella, Klebsiella, Proteus, Serratia, Citrobacter and Edwardsiella, Erwinia) [15,16]. Secondary sources of microbes in honey, are likely to be the same as for other foods (humans, animals, water, soil, air and/or processing facilities or equipment, such as honey harvesting or storage containers) and, consequently, appropriate standards of hygiene must be applied in all operations involving honey handling [17]. Clostridium spp., Corynebacterium spp., Bacillus spp. and Pseudomonas spp. are bacteria commonly found in soil. Brochotrix spp., Citrobacter spp., Enterobacter spp., Lactobacillus spp., Lactococcus spp., Pediococcus spp., Listeria spp. and Flavobacterium are found in plants and plant products. Air and dust are important source of Bacillus spp., Clostridium spp. and Micrococcus spp. species. Saccharomycses and Torula yeasts can be found in high-moisture sugars [15,16]. The ripe honeys contain mostly molds and different species of sporulating bacteria Clostridium and Bacillus (because of their possibility of development in conditions of limited water availability; aw 0.57–0.62) and osmophilic yeast (Schizosaccharomyces, Hansenula, Torula, Pichia, Nemaiospora, Schwanniomyces and Rhodotorula). The nectar honey types are dominated by the Saccharomyces yeasts, represented by a dozen or so strains and aerobic and anaerobic bacterial spores [15]. On the other hand, another safety issue identified in honey are chemical pollution, the source of which are improper practices in agribusiness. Many authors indicated that the presence of pesticides in honeys can be serious problem [8].
In European countries there is a lack of the national microbiological guidelines for assessing the microbiological quality of honey and bee products, abolished by European Union regulations. The microbiological quality of honey is assessed according to the criteria adopted for food and feed in accordance with Regulation (EC) No 178/2002 of the European Parliament and of the Council [18] applicable to all stages of production, processing and distribution of food and feed as well as in Food Code (Codex Alimentarius) where the microbiological quality of honey relates to the criteria established in accordance with the principles and guidelines for the determination and application of microbiological and hygiene criteria related to food [19]. The quality of the honey is defined by physicochemical criteria well compiled by European Directive 2014/63/EU, while the microbiological aspects are ignored [20]. Under this Directive also there is no requirement for disclosing the country of origin on a honey package [21]. Where the product is a blend originating from more than one EU and/or non-EU country, it may be designated as an “blend of EU honey”, “blend of non-EU honey” or “blend of EU and non-EU honey” [22]. Because of the heavy reliance of the product quality on its origin, the consumer should not be “geographically misguided”. The modern consumer has the possibility to choose products from a very wide range and markets in many retail countries, and the role of information provided with the product is growing. When making a choice, the consumer can take into account the origin of the product, its safety and health quality [23].
The aim of the research was to evaluate and classify honey samples from different climatic zones on the basis physicochemical properties and the evaluation of their microbiological quality and safety. Most of the research publications concerns physicochemical properties of honey in the light of the EU Directive [21]. Honey microbiota research is mainly focused on the analysis of the occurrence of toxynogenic Clostridium botulinum bacteria [12,24,25,26]. Only the small number of current scientific publications on the microbiological and hygienic quality of honeys are available. To our best knowledge these are the first honey study that present a microbiological assessment in comparison to the physicochemical requirements.

2. Results

2.1. Pollen Analysis

The analyzed honey samples had been declared by the beekeepers as varietals honeys, based on organoleptic features as well as bee forage (available to bees flowering surrounding plants) (Table 1). However, the pollen samples analyze showed that (in all samples predominant pollen > 45% was not identified), so all the samples were considered as multifloral honeys [27,28,29].

2.2. Physicochemical Analysis

The physicochemical results of honey samples analysis were presented in detail in Table 2. The pH value of Polish and Thai honey samples ranged from 3.38 to 4.43 and 3.44 to 4.90, respectively. The Polish honey samples were characterized by very balanced water activity. The water activity of Polish honey samples ranges from 0.501 to 0.578. The median value of aw (0.535) differed from the average value (0.534) by 0.001 which is confirmed by the even distribution of this parameter in the tested samples. Thai honeys samples aw values were found significantly higher (p < 0,05) than Polish one; the results ranged from 0.553 to 0.673, the average and median value of aw were 0.605 and 0.609, respectively. Moreover, most of the Polish honey samples contained up to 20% of water. This limit was exceeded in samples P1, P8 and P16 (slightly, by 0.5–0.9%) and in sample P3 (by 2%). For the Thai honeys, thirteen of the twenty samples were above the limit (by more than 6% for samples T12 and T20).

2.3. Microbiological Analysis

The results of microbiological analysis were shown at Figure 1. The bacteria of Salmonella spp. was not detected in any of the tested honey samples. The mean values obtained for the total number of mesophilic aerobic bacteria (TVC) in the examined honey samples were 0.98 log cfu/g and 1.13 log cfu/g for Polish and Thai honeys, respectively. The yeasts and molds were found in 11 of Thai and 9 of Polish honeys samples, the mean count of them was respectively 1.24 and 0.11 log cfu/g. Thai honey samples (30%) contained a small amount of lactic acid bacteria (LAB), not exceed 0.97 log cfu/g. The Polish honey samples (50%) contains the population of LAB ranged from 0.30 log cfu/g up to 1.75 log cfu/g. In the case of mesophilic spores of the Bacillus spp. bacterium, the mean value obtained in the Polish honeys samples was 0.82 log cfu/g, and for Thai honey samples the mean value was 0.98 log cfu/g.
Figure 2 shows the relationship between water content and the number of yeasts and molds only in samples where yeasts and molds have been detected. In the case of three samples of Polish honeys, the excessive water content did not affect the number of yeasts and molds, which not exceeding 0.5 log cfu/g, while in six samples of Thai honeys, the increased water content caused high number of yeast and molds 1.8–3.01 log cfu/g. The correlation between the number of yeast and molds and water content was 0.575 and 0.428 for Polish and Thai honey samples respectively, it has also clearly seen that in some cases water content can affect yeast and molds dynamic growth.
Figure 3 and Figure 4 shows the results of cluster analysis classification (variable and cases, respectively) of the Polish and Thai honeys based on results of four microbiological analysis (TVC, count of yeast and molds, LAB, Bacillus spp.).
Based on an analysis of the classification of variables for the Polish and Thai honeys, the counts of mesophilic aerobic microorganisms and the Bacillus spp. form mutually correlated pairs and the results can be used interchangeably for microbiological quality evaluation of the honeys irrespective of the cluster merging method. The variable—LAB represented a distinct cluster, slightly correlated to the variables: TVC and the Bacillus spp., particularly for the Thai honeys. The variable—yeasts and molds represented a different cluster (Figure 3).
Figure 4 shows the classification of cases of the Polish and Thai honeys. The analysis identified two main groups (Cluster I and Cluster II) of honeys with a distance of 20 between cluster centers, which attests to a significant variation of the features of the honeys. Ten Thai honey samples, varying significantly in their microbiological contamination levels, were classified as belonging to the first major cluster (CI). Similar to each other in this group were the honeys: T13 and T14 as well as T1 and T19. The second major cluster (CII) was formed from two smaller clusters marked as C3 and C4, with a binding length of approximately nine. The C3 groups of honey with similar microbial contamination as in group (CI) included 11 honey samples: 4 Thai and 7 Polish. The most similar in terms of microbiological quality in this group were multiflorous Polish honeys P14 and P15, Thai honeys constituted 36.4% in this group (T4, T10, T15, T16). The last, third, cluster (C3) consisting 19 samples featured the smallest variation in their microbiological quality including contamination. Thai honeys (T2, T3, T6, T12, T17, T20) represented 31.5% of the samples from this cluster. The classification of honey originating from the two climatic zones demonstrated significant differences −70% of the tropical honeys (C1 and 4 samples of C3) had different microbiological quality than honey from temperate climate zone.

3. Discussion

Pollen analysis of honey is widely used to verify the claimed geographic and floral origin of honey samples [30]. Beekeepers do not always perform pollen analysis due to cost, lack of the necessary equipment and knowledge. Based on the literature data, it is known that beekeepers declare the type of honey on the basis of the organoleptic characteristics and properties of the bee forage [31,32]. On the other hand, honey with the presence of one type of pollen (predominant pollen) at the level above 45% [27,28,29] can be considered as varietal honey. In case of the present study, future, more detailed analyzed should be performed to determinate the predominant pollen and type of tested honey samples.
The data obtained from physicochemical study indicated that honey samples originating from the two different climatic zones (transitional between the continental moderate and oceanic moderate one vs. monsoon tropical) differ in case of pH, water activity and water content. In the present study, the pH value of the honeys of Polish and Thai samples (Table 2) were 3.3–4.9, which is consistent with other authors findings [15,16,33]. None of the normative acts imposes a minimum or maximum pH or water activity values. The pH value depends on the acids present in the honey (e.g., acetic, butyric, citric, formic, gluconic, lactic and malic). It was found, that the gluconic acid present in the honey in the range of 0.23–0.98% plays the greatest role in antimicrobial activity of honey [7].
The water content decreases as the water evaporates during the honey maturation period. In the tropical climate zone, air humidity is high, which interferes with the natural ripening process. Due to the high sugar content (especially fructose and glucose) of honey, the osmotic pressure of honey is usually high leading to low water activity (aw) in reported range 0.56–0.62. The limiting water activity for growth of osmotolerant yeast is about 0.61–0.62 and much other microorganisms. Knowledge of water activity of honey is also needed to predict moisture exchange with the environment, since water activity is the driving force behind water transfer from/to honey [16,34]. Water activity (aw) and water content (%) of Polish honey samples was lower than Thai honey samples. The average and median value of Polish samples were identical, 0.535 and 0.534, respectively; in Thai samples increased value of aw (average 0.605 and median 0.609 could lead to an increase in yeasts and molds count.
Moreover, most of the Polish honey samples contained up to 20% of water, which is in line with EU Directive [21]. This limit was exceeded in four Polish samples (P1, P3, P8 and P16) and in the thirteen samples of Thai honeys, which indicates the unripe of Thai honeys [17,21,27].
If the water content is high, the presence of a single mold cell may be sufficient to initiate a fermentation, causing adverse sensory changes. When yeasts and molds grow in honey, they break down monosaccharides which indicates a fermentation. As a result of this process, ethyl alcohol and carbon dioxide are formed. This is followed by a honey acidification due to the growth of bacteria, which cause the formation of acetic acid and non-volatile acids that have a characteristic aroma [35,36,37]. In the research of Madras–Majewska et al., 2016 one of three Thailand multifloral honey samples was negatively evaluated with respect to sensory attributes because of the tart, burning and sharp taste. The changes were connected with above 20% of water content (21.7%) and high count of yeasts and molds (more than 3.5 log cfu/g) [9].
In the present study the slight correlation between the number of yeast and molds and water content also was found. According to Snowdon and Cliver, 1996, yeast and molds activity is limited by the content of free water, causing that honey from humid regions is more prone to be contaminated with osmophilic yeasts [15]. Piana et al., 1991 reported yeast growth only in honey samples with a water activity < 0.65 [38]. The presence of yeasts and molds is also often caused by cross-contamination from product handling equipment [17].
Currently, the situation with the honey legislation complicates even more if we consider that some countries issues are national provisions, decisions and guidelines filling the gap in European and International legislation, despite the EU recommendations. Most of them set limits to define the physicochemical, organoleptic, microbiological and microscopic characteristics of monofloral honey, some have provisions regarding the country where the honey has been harvested and several others are differentiated from the set criteria. The differences among the national provisions enhance the difficulties of the applicability of honey regulations and make the necessity of uniformity of honey legislation [22]. In the repealed by EU low Polish legislation established in relation to bee products which are not honeys (pollen, propolis, bee bread), the microbiological criteria for the total number of mesophilic aerobic microorganisms were set at maximum level 5 × 104 (4.7 log cfu/g); not more than 5 × 102 (2.7 log cfu/g) of yeast and molds; not more than 100 (2 log cfu/g) of Staphylococcus aureus and aerobe Bacillus cereus [39]. Mexico, the major world honey exporter countries use no mandatory standard (Mexico NMX-036-Normex-2006) establishing the accepted presence of no more than 1000 cfu/g (3.0 log cfu/g) of non-pathogenic bacteria and up to 100 cfu/g (2.0 log cfu/g) of yeast and molds. Similarly, in Argentina, commercial honeys are ruled by Código Alimentario Argentino (CAA) and Mercado Común del Sur rules. The maximum level of microorganisms allowed by this legislation for molds and yeast with trading purposes is 100 cfu/g of honey. Likewise, the legislation does not allow the presence of Salmonella and Shigella bacteria or total coliforms in honey [16,17].
In the present study yeasts and molds were found in 11 of Thai honeys samples (55%) and the mean count of them was 1.24 log cfu/g. Eight from eleven positive honey samples (T5, T8, T9, T11, T13, T14, T18 and T19) were above the national microbiological limit (2.0 log cfu/g) [20]. The growth of the microorganisms was probably caused by elevated values of water activity of in range 0.569–0.647 and water contents of five of these eight samples (T5, T9, T11, T13 and T14) contained an above-the-limit twenty percent quantity of water (21.65–25.07%). Gomes et al., 2010 reported yeast and mold count on a similar level with Thai honey with average count of 1.53 log cfu/g [10]. Iurlina et al., 2005 showed count yeast and molds equal to or lower than 2.67 log cfu/g (57% of honey samples) [33]. Similar level of mold 2.0–2.2 log cfu/g (17% of samples) and yeast 2.07–2.99 log cfu/g examined separately showed in Mexico honey samples [16]. In the present study yeasts and mold were found in 45% and the mean count of them was 0.26 log cfu/g in the Polish honey samples. Only three samples contained some small number of yeasts and molds (0.1–0.56 log cfu/g) and above-the-limit of water content 20.5–22.0%.
The mean values obtained for the total number of mesophilic aerobic bacteria in the tested honey samples were 0.98 and 1.13 log cfu/g for Polish and Thai honeys, respectively. Very similar results found Gomes et al., 2010, who have found that all tested honey samples were contaminated <1.0 log cfu/g of TVC except one sample which was contaminated on the level 1.30 log cfu/g [10]. Iurlina et al., 2005 reported higher level of contamination (average 2.38 log cfu/g), but still in range of no mandatory Mexico NMX-036-Normex-2006 standard [16,33]. According to Mexico standard (no more than 1000 cfu/g of nonpathogenic bacteria and up to 100 cfu/g of yeast and molds) 40.5% of honey samples from Mexico did not meet the specification in the case of aerobic mesophilic bacteria [16].
We have found that the number of LAB were higher in the Polish samples in comparison to the Thai honey. Lactic acid bacteria are recognized as safe and play important role in the preservation of the product since many of them have the ability to produce antimicrobial agents such as organic acids and bacteriocins that can inhibit or destroy pathogenic bacteria although this characteristic may depending on the type of LAB. Very few (2%) of Mexico honey samples contain LAB and only 15.79% of the samples contained more than 2 log cfu/g [16]. Earlier reports also detected the presence of some vegetative, non-spore forming lactic acid bacteria in raw honey [40,41,42,43]. Lactobacillus, Lactococcus and Leuconostoc genera were isolated from flowers, plant surfaces and plant associated products suggesting that LAB which present in honey may come from plant sources and the bees [44].
However, in the case of mesophilic spores of the Bacillus spp. bacterium, in the present study the mean value obtained in the case of 18 (90%) Polish honeys samples was 0.82 log cfu/g, for 16 (80%) of the 20 Thai honeys samples the mean value was 0.98 log cfu/g (Figure 2).
Moreover, other authors indicated that some opportunistic spore-forming bacteria, molds and yeasts are typically found in honey, often at low numbers, while spores can persist indefinitely [43]. Spores from the Bacillus genus are regularly found in honey. Iurlina et al., 2005 and Iurlina, et al., 2006 showed 23–27% Bacillus spp. positive honey samples which were identified as B. cerus, B. pumilus and B. laterosporus [33,45]. Some authors reported that the potential toxigenic effects of Bacillus were achieved with 104 spores per g of honey [7]. Consumption of foods that contain more than 104 spores or vegetative cells B. cereus per gram may results in food poisoning. B cereus can produce several toxins, most important being the emetic toxin and the HBL enterotoxin complex. Investigation of honey samples from Argentina, conducted various methods, found B. cereus in 27% samples and other species of Bacillus 14% of honey samples [46]. In another investigation the most frequently isolated from honey samples was Bacillus amylioliquefacien [20].
The microbiota of ripening honey can also pose a threat to human health, due to the viable from of pathogens as: Klebsiella spp., Salmonella spp. and Shigella spp. [47]. In the present study, the Salmonella spp. was not detected in any of the tested honey samples, which is in accordance with Regulation (EC) No 2073/2005 [48]. The same results were obtained by other authors when examining honey samples [10,16,33] and the samples coming from honey drums [17]. However, up to date no vegetative forms of disease-causing bacterial species have been found in honey. Bacteria do not replicate in honey and as such high numbers of vegetative bacteria could indicate recent contamination from a secondary source [15].

4. Materials and Methods

4.1. Research Material

The research material consisted of ripe honey samples purchased directly from beekeepers in apiaries in order to avoidance of manipulative contamination research material. The time from acquisition from the apiary to analysis did not exceed 3 months. Each sample was stored at room temperature in darkness for future study. Twenty Polish samples were representing a climate transitional between the continental moderate and oceanic moderate one (P1–P20). Honeys purchased in Thailand represented the monsoon tropical climate zone (T1–T20). In Table 1 characteristic of honey samples was shown, according to producer declaration.

4.2. Methods of Analysis

4.2.1. Pollen Analysis

Conformity of the honey samples was determined by palynologic microscopic observations consistent with the method recommended by the International Commission for Bee Botany and by the International Honey Commission Variety and conformity with the Polish standard [28,29].

4.2.2. Physicochemical Analysis

The pH value of honey was determined using Lab 860 pH meter (SI Analytics GmbH, Schott Instruments, Germany). Samples of honey were prepared by dissolving 10 g of honey in 75 mL of distilled water carbon dioxide free [28].
The water activity was measured at 25 ± 0.2 °C using Aqualab TE series 4 analyzer (Decagon Devices, Pullman, Washington, DC, USA), in a temperature stable sampling environment, calibrated with saturated salt solutions in the aw range of 0.40–0.70. This device operates based on electronic dew-point measurements. AquaLab analyzer continues the analyses of water activity until the difference of three consecutive measurements is less than 0.0005 aw [34].
Water content in undiluted honeys was established using PAL-22S refractometer (ConbestCo, Kraków, Poland)—three/four measurements have been performed.

4.2.3. Microbial Counts Analysis

For each sample a mass of 10 g of honey was weighed aseptically and homogenized for 1 min with 90 mL of Buffered Peptone Water (BPW) in Stomacher 400 instrument (IUL Instruments, Königswinter, Germany). To perform 10-fold dilutions BPW has been used. The microbiological quality of honey samples was tested using the reference plate method, performing analysis on two parallel plates. The number of microorganisms was expressed as colony forming units per gram of honey (cfu/g).
Determination of TVC was performed using the ISO standard [49] with Plate Count Agar (PCA Biokar Diagnostics, Austria). The plates were incubated at 30 °C for 48–72 h days. The counts of yeast and mold were determined on Rose Bengal Chloramphenicol Agar (RBC Agar, Becton Dickinson and Co., Franklin Lakes, NJ, USA), after incubation at 25 °C for 5 days [50]. The number of mesophilic cells of LAB was determined with a pour plate method using MRS agar (de Man, Rogosa and Sharpe, Biokar Diagnostic, Wien, Austria). Plates were incubated anaerobically at the temperature 30 °C for up to 72 h [51]. The presence of mesophilic Bacillus spp. bacteria, determined using Mannitol Egg Yolk Agar supplemented with Polymyxin—Mossel (MYP) Agar (Bio-Rad, Watford, UK), after incubation at 37 °C for 24–48 h [52].

4.2.4. Bacterial Detection Analysis

Detection of Salmonella species in 25 g was performed instrumentally using BacTrac 4300 screening impedance method (SyLabGeräte GmbH, Purkersdorf, Austria). Analysis has been performed using BiMedia 201C Salmonella culture medium (modif. Rappaport-Vassiliadis media), measurement cells were incubated at 40 °C for 24 h, threshold 10% for the E-value (Electrode-value, impedance around electrode).

4.2.5. Data Analysis

The cluster analysis method and the Statistica 13.3 PL software were used to classify the results of microbiological analyses of honey [53]. The distance between clusters was measured by Euclidean distance function, while the Ward method was used to bind the clusters. The Ward method uses the assumptions of variance analysis and aims to minimize the sum of deviations within clusters. As a result of joining cluster pairs, the pair that gives the cluster with the minimum differentiation is chosen. Error sum of squares (ESS) is a measure of the difference to the mean value. The test t-student was used to assess differences between the samples (α = 0.05) and the linear coefficient of Pearson to assess the correlation between the variables.

5. Conclusions

The study has revealed that geographical origin of honey may significantly affect the quality and safety of honey especially from different climatic zones. Our study indicated that the Polish and the Thai honey samples, although considered as multiflower, were differ in case of count yeast and mold, which was affected by water content. It was considered that water content can be the most informative and handy marker of the microbiological quality of honeys. The classification of honeys samples found that the total number of viable microorganisms could be used as a measure alternative to the count of Bacilus spp., which is important hint for the legislator. Moreover, 70% of the tropical honey samples had different microbiological quality than honey from temperate climate zone. Analysis of lactic acid bacteria showed temperate climate zone honeys as a source of beneficial bacteria in the diet.
Due to the lack of EU and World regulations on the presence of microorganisms in honey, the attention should be given to the origin of the honey. The microbiological criteria used in the legislations of different countries are the result of good production and hygiene practices applied by honey producers and constitute an indispensable reference point in assessing the quality of honey and show that microbiological control of honey is needed in order to ensure the safety of consumers.

Author Contributions

Conceptualization, E.R. and D.Z.; methodology, E.R.; formal analysis, D.T., E.R. and A.Ł.; resources, B.M.-M.; writing—original draft preparation, E.R.; writing—review and editing, D.Z.; visualization, E.R.; funding acquisition, B.M.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Polish Ministry of Science and Higher Education within funds of the Faculty of Human Nutrition and Consumer Sciences, Warsaw University of Life Sciences (WULS), for scientific research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data and figures generated or used during the study appear in the submitted article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Nguyen, H.T.L.; Panyoyai, N.; Kasapis, S.; Pang, E.; Mantri, N. Honey and its role in relieving multiple facets of atherosclerosis. Nutrients 2019, 11, 167. [Google Scholar] [CrossRef] [Green Version]
  2. Lashani, E.; Davoodabadi, A.; Dallal, M.M.S. Some probiotic properties of Lactobacillus species isolated from honey and their antimicrobial activity against foodborne pathogens. Vet. Res. Forum 2020, 11, 121–126. [Google Scholar]
  3. Esawy, M.A.; Awad, G.E.; Ahmed, E.F.; Danial, E.N.; Mansour, N.M. Evaluation of honey as new reservoir for probiotic bacteria. Adv. Food Sci. 2012, 34, 72–81. [Google Scholar]
  4. Begum, S.B.; Roobia, R.R.; Karthikeyan, M.; Murugappan, R.M. Validation of nutraceutical properties of honey and probiotic potential of its innate microflora. LWT Food Sci. Technol. 2015, 60, 743–750. [Google Scholar] [CrossRef]
  5. Asama, T.; Arima, T.-H.; Gomi, T.; Keishi, T.; Tani, H.; Kimura, Y.; Tatefuji, T.; Hashimoto, K. Lactobacillus kunkeei YB38 from honeybee products enhances IgA production in healthy adults. J. Appl. Microbiol. 2015, 119, 818–826. [Google Scholar] [CrossRef]
  6. Hill, C.; Guarner, F.; Reid, G.; Gibson, G.R.; Merenstein, D.J.; Pot, B.; Morelli, L.; Canani, R.B.; Flint, H.J.; Salminen, S.; et al. Expert consensus document: The International Scientific Association for Probiotics and Prebiotics consensus statement on the scope and appropriate use of the term probiotic. Nat. Rev. Gastroenterol. Hepatol. 2014, 11, 506–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Saranraj, P.; Sivasakthi, S. Comprehensive review on honey: Biochemical and medicinal properties. J. Acad. Res. 2018, 6, 165–181. [Google Scholar]
  8. Lignor, M.; Bukowska, M.; Ratiu, I.A.; Gadzała-Kopciuch, R.; Buszewski, B. Determination of neonicotinoids in honey samples originated from Poland and other world countries. Molecules 2020, 25, 5817. [Google Scholar] [CrossRef] [PubMed]
  9. Madras-Majewska, B.; Rosiak, E.; Jaworska, D.; Kulesza, K.; Wasiak-Zys, G.; Teper, D. Comparison of selected quality characteristics of domestic and Thailand multifloral honeys. Vet. Med. 2016, 72, 593–656. [Google Scholar] [CrossRef] [Green Version]
  10. Gomes, S.; Dias, L.G.; Moreira, L.L.; Rodrigues, P.; Estevinho, L. Physicochemical, microbiological and antimicrobial properties of commercial honeys from Portugal. Food Chem. Toxicol. 2010, 48, 544–548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Ratiu, I.A.; Al-Suod, H.; Bukowska, M.; Ligor, M.; Buszewski, B. Correlation sudy of honey regarding their physicochemical properties and sugars and cyclitols content. Molecules 2020, 25, 34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Finola, M.S.; Lasagno, M.C.; Marioli, J.M. Microbiological and chemical characterization of honeys from central Argentina. Food Chem. 2007, 100, 1649–1653. [Google Scholar] [CrossRef]
  13. Gheldof, N.; Wang, X.H.; Engeseth, N.J. Identification and quantification of antioxidant components of honeys from various floral sources. J. Agric. Food Chem. 2002, 50, 5870–5877. [Google Scholar] [CrossRef] [PubMed]
  14. Vásquez, A.; Forsgren, E.; Fries, I.; Paxton, R.J.; Flaberg, E.; Szekely, L.; Olofsson, T.C. Symbionts as major modulators of insect health: Lactic acid bacteria and honeybees. PLoS ONE 2012, 7, e33188. [Google Scholar] [CrossRef]
  15. Snowdon, J.A.; Cliver, D.O. Microorganisms in honey. Int. J. Food Microbiol. 1996, 31, 1–26. [Google Scholar] [CrossRef]
  16. Vázquez-Quiñones, C.R.; Moreno-Terrazas, R.; Natividad-Bonifacio, I.; Quiñones-Ramírez, E.I.; Vázquez-Salinas, C. Microbiological assessment of honey in México. Rev. Argent. Microbiol. 2018, 50, 75–80. [Google Scholar] [CrossRef]
  17. Fernandez, L.A.; Ghilardi, C.; Hoffmann, B.; Busso, C.; Gallez, L.M. Microbiological quality of honey from the Pampas Region (Argentina) throughout the extraction process. Rev. Argent. Microbiol. 2017, 49, 55–61. [Google Scholar]
  18. European Comission. Regulation (EC) No 178/2002 of the European Parliament and of the Council of 28 January 2002 Laying Down the General Principles and Requirements of Food Law, Establishing the European Food Safety Authority and Laying Down Procedures in Matters of Food Safety; European Comission: Brussels, Belgium, 2002; pp. 1–24. [Google Scholar]
  19. WHO; FAO. Codex 2019. Codex Alimentarius Standard for Honey; CXS 12-1981; FAO: Italy, Rome, 1981. [Google Scholar]
  20. Sinacori, M.; Francesca, N.; Alfonso, A.; Cruciata, M.; Sannino, C.; Settanni, L.; Moschetti, G. Cultivable microorganisms associated with honeys of different geographical and botanical origin. Food Microbiol. 2014, 38, 284–294. [Google Scholar] [CrossRef] [Green Version]
  21. European Comission. Directive 2014/63/EU of the European Parliament and of the Council of 15 May 2014 Amending Council Directive 2001/110/EC Relating to Honey; European Comission: Brussels, Belgium, 2014; pp. 1–5. [Google Scholar]
  22. Thrasyvoulou, A.; Tananaki, C.; Goras, G.; Karazafiris, E.; Liolios, V.; Kanalis, D.; Gounari, S. Legislation of honey criteria and standards. J. Apic. Res. 2018, 15, 88–96. [Google Scholar] [CrossRef]
  23. Grunert, K.G.; Aachmann, K. Consumer reactions to the use of EU quality labels on food products: A review of the literature. Food Control 2016, 59, 178–187. [Google Scholar] [CrossRef]
  24. Nevas, M.; Hielm, S.; Lindstrom, M.; Horn, H.; Koivulehto, K.; Korkeala, H. High prevalence of Clostridium botulinum types A and B in honey samples detected by polymerase chain reaction. Int. J. Food Microbiol. 2002, 72, 45–52. [Google Scholar] [CrossRef]
  25. Nevas, M.; Lindstrom, M.; Hautamaki, K.; Puoskari, S.; Korkeala, H. Prevalence and diversity of Clostridium botulinum types A, B, E and F in honey produced in the Nordic countries. Int. J. Food Microbiol. 2005, 105, 145–151. [Google Scholar] [CrossRef] [PubMed]
  26. Küplülü, Ö.; Göncüoğlu, M.; Özdemir, H.; Koluman, A. Incidence of Clostridium botulinum spores in honey in Turkey. Food Control 2006, 17, 222–224. [Google Scholar] [CrossRef]
  27. Ministry of Agriculture and Rural Development of Poland. Regulation of the Minister of Agriculture and Rural Development of 29 May 2015 [Rozporządzenie Ministra Rolnictwa i Rozwoju Wsi z Dnia 29 Maja 2015 Zmieniające Rozporządzenie w Sprawie Szczegółowych Wymagań w Zakresie Jakości Handlowej Miodu]; Ministry of Agriculture and Rural Development of Poland: Warsaw, Poland, 2015; p. 850. (In Polish)
  28. Louveaux, J.; Maurizio, A.; Vorwohl, G. Methods of Melissopalynology. Bee World 1978, 59, 139–157. [Google Scholar] [CrossRef]
  29. Ministry of Agriculture and Rural Development of Poland. Regulation of the Minister of Agriculture and Rural Development of January 14 2009 [Rozporządzenie Ministra Rolnictwa i Rozwoju Wsi z Dnia 14 Stycznia 2009 r. w Sprawie Metod Analiz Związanych z Dokonywaniem Oceny Miodu]; Ministry of Agriculture and Rural Development of Poland: Warsaw, Poland, 2009; p. 94. (In Polish)
  30. Sniderman, J.K.; Matley, K.A.; Haberle, S.G.; Cantrill, D.J. Pollen analysis of Australian honey. PLoS ONE 2018, 13, e0197545. [Google Scholar] [CrossRef] [Green Version]
  31. Piekut, J.; Witkowska, A.; Borawska, M.; Hejft, R. Próba zastosowania analizy spektrofotometrycznej w bliskiej podczerwieni do rozróżniania miodów. Bromatol. Chem. Toksykol. 2000, 33, 73–78. [Google Scholar]
  32. Puścion-Jakubik, A.; Borawska, M. Varietal bee honey—Principal and accompanying pollen as a basis for their classification. Probl. Hig. Epidemiol. 2016, 97, 275–278. [Google Scholar]
  33. Iurlina, M.O.; Fritz, R. Characterization of microorganisms in Argentinian honeys from different sources. Int. J. Food Microbiol. 2005, 105, 297–304. [Google Scholar] [CrossRef] [PubMed]
  34. Zamora, M.C.; Chirife, J.; Roldan, D. On the nature of the relationship between water activity and % moisture in honey. Food Control 2006, 17, 642–647. [Google Scholar] [CrossRef] [Green Version]
  35. Guau, A.; Akalakou, C.; Degrande, A.; Biloso, A. Determinants of consumer preferences for honey in the Democratic Republic of Congo. J. Food Prod. Market. 2014, 20, 476–490. [Google Scholar] [CrossRef]
  36. Ismaiel, S.; Al-Kahtani, S.; Adgaba, N.; Al-Ghamdi, A.A.; Zulail, A. Factors that affect consumption patterns and market demands for honey in the Kingdom of Saudi Arabia. Food Nutr. Sci. 2014, 5, 1725–1737. [Google Scholar] [CrossRef] [Green Version]
  37. Pocol, C.B.; Bolboac, S.D. Perceptions and trends related to the consumption of honey. A case study of North-West Romania. Int. J. Consum. Stud. 2013, 37, 642–649. [Google Scholar] [CrossRef]
  38. Piana, M.L.; Poda, G.; Cesaroni, D.; Cuetti, L.; Bucci, M.A.; Gotti, P. Research on microbial characteristics of honey samples of Udine province. Riv. Stor. Ital. 1991, 20, 293–301. [Google Scholar]
  39. Ministry of Health of Poland. Regulation of the Minister of Health of January 13 2003 [Rozporządzenie Ministra Zdrowia z Dnia 13 Stycznia 2003 r. w Sprawie Maksymalnych Poziomów Zanieczyszczeń Chemicznych i Biologicznych, Które Mogą Znajdować Się w Żywności, Składnikach Żywności, Dozwolonych Substancjach Dodatkowych, Substancjach Pomagających w Przetwarzaniu Albo na Powierzchni Żywności]; Ministry of Health of Poland: Warsaw, Poland, 2003; p. 326. (In Polish)
  40. Forsgren, E.; Olofsson, T.C.; Vasquez, A.; Fries, I. Novel lactic acid bacteria inhibiting Paenibacillus larvae in honey bee larvae. Apidologie 2009, 41, 99–108. [Google Scholar] [CrossRef] [Green Version]
  41. Hosny, I.M.; El-Ghani, S.; Nadir, A.S. Nutrient composition and microbiological quality of three unifloral honeys with emphasis on processing of honey probiotic yoghurt. Glob. Vet. 2009, 3, 107–112. [Google Scholar]
  42. Bahiru, B.T.; Mehari, T.M.; Ashenafi, M. Yeast and lactic acid flora of an indigenous Ethiopian honey wine. Variations within and between production units. Food Microbiol. 2006, 23, 277–282. [Google Scholar] [CrossRef]
  43. Ruiz-Argueso, T.; Rodriguez-Navarro, A. Microbiology of ripening honey. Appl. Environ. Microbiol. 1975, 30, 893–896. [Google Scholar] [CrossRef]
  44. Tavaria, F.K.; Dahl, S.; Carballo, F.J.; Malcata, F.X. Amino acid catabolism and generation of volatiles by lactic acid bacteria. J. Dairy Sci. 2002, 85, 2462–2470. [Google Scholar] [CrossRef]
  45. Iurlina, M.O.; Saiz, A.I.; Fuselli, S.R.; Fritz, R. Prevalence of Bacillus spp. in different food products collected in Argentina. LWT Food Sci. Technol. 2006, 39, 105–110. [Google Scholar] [CrossRef]
  46. Lopez, A.C.; Alippi, A.M. Phenothypic and genotyphic diversity of Bacillus cereus isolates recovered from honey. Int. J. Food Microbiol. 2007, 117, 175–184. [Google Scholar] [CrossRef] [PubMed]
  47. Grabowski, N.T.; Klein, G. Microbiology and foodborne pathogens in honey. Crit. Rev. Food Sci. Nutr. 2017, 57, 1852–1862. [Google Scholar] [PubMed]
  48. European Comission. Commission Regulation (EC) No 2073/2005 of 15 November 2005 on microbiological criteria for foodstuffs. Off. J. Eur. Union. 2005, L338, 1–26. [Google Scholar]
  49. ISO. ISO 4833-1-2:2013. Microbiology of Food and Animal Feeding Stuffs—Horizontal Method for the Enumeration of Microorganisms—Part 1: Colony Count at 30 Degrees C by the Pour Plate Technique; ISO: Geneva, Switzerland, 2020. [Google Scholar]
  50. ISO. ISO 21527-2:2008. Microbiology of Food and Animal Feeding Stuffs—Horizontal Method for the Enumeration of Yeasts and Molds. Part 2: Colony Count Technique in Products with Water Activity less than or Equal to 0,95; ISO: Geneva, Switzerland, 2018. [Google Scholar]
  51. ISO. ISO 15214:1998. Microbiology of Food and Animal Feeding Stuffs—Horizontal Method for the Enumeration of Mesophilic Lactic Acid Bacteri–Colony-Count Technique at 30 Degrees C; ISO: Geneva, Switzerland, 2021. [Google Scholar]
  52. ISO. ISO 7932:2004. Microbiology of Food and Animal Feeding Stuffs—Horizontal Method for the Enumeration of Presumptive Bacillus cereus—Colony-Count Technique at 30 Degrees C; ISO: Geneva, Switzerland, 2004. [Google Scholar]
  53. StatSoft Inc. STATISTICA (Data Analysis Software System), version 13.3; StatSoft: Tulsa, OK, USA, 2011; Available online: www.statsoft.com (accessed on 31 March 2021).
Figure 1. Results of microbiological analyses of Polish and Thai honeys.
Figure 1. Results of microbiological analyses of Polish and Thai honeys.
Molecules 26 02361 g001
Figure 2. Water content [%] and the number of yeasts and molds [log cfu/g] in Polish (P) and Thai (T) honeys, in which yeast and molds were detected.
Figure 2. Water content [%] and the number of yeasts and molds [log cfu/g] in Polish (P) and Thai (T) honeys, in which yeast and molds were detected.
Molecules 26 02361 g002
Figure 3. Dendrogram of hierarchical cluster analysis of microbiological variables determined by the Ward method in Thai honeys (a) and Polish honeys (b).
Figure 3. Dendrogram of hierarchical cluster analysis of microbiological variables determined by the Ward method in Thai honeys (a) and Polish honeys (b).
Molecules 26 02361 g003
Figure 4. Dendrogram of hierarchical cluster analysis of Polish and Thai honeys cases determined by the Ward method.
Figure 4. Dendrogram of hierarchical cluster analysis of Polish and Thai honeys cases determined by the Ward method.
Molecules 26 02361 g004
Table 1. Analyzed honey samples according to beekeepers’ declarations.
Table 1. Analyzed honey samples according to beekeepers’ declarations.
Sample No.Producer DeclarationSample No.Producer Declaration
P1goldenrod nectar honeyT1coffee nectar honey
P2raspberry nectar honeyT2coffee nectar honey
P3rapeseed nectar honeyT3coffee nectar honey
P4rapeseed nectar honeyT4coffee nectar honey
P5linden nectar honeyT5coffee nectar honey
P6linden nectar honeyT6wild forest nectar honey
P7multiflorous nectar honeyT7wild forest nectar honey
P8multiflorous nectar honeyT8wild forest nectar honey
P9multiflorous nectar honeyT9wild forest nectar honey
P10buckwheat nectar honeyT10longan nectar honey
P11multiflorous nectar honeyT11longan nectar honey
P12dandelion nectar honeyT12longan nectar honey
P13forest nectar honeyT13longan nectar honey
P14multiflorous nectar honeyT14longan nectar honey
P15raspberry nectar honeyT15longan nectar honey
P16buckwheat nectar honeyT16longan nectar honey
P17multiflorous nectar honeyT17longan nectar honey
P18acacia nectar honeyT18lychee nectar honey
P19acacia nectar honeyT19lychee nectar honey
P20forest nectar honeyT20lychee nectar honey
Table 2. Physicochemical properties of Polish and Thai honey samples.
Table 2. Physicochemical properties of Polish and Thai honey samples.
Polish and Thai Honey Sample
No.
of Sample
pH ± SDWater Content
[%] ± SD
Water
Activity ± SD
P14.07 ± 0.0220.90 ± 0.040.566 ± 0.00
P23.47 ± 0.0217.20 ± 0.050.544 ± 0.00
P33.38 ± 0.0422.00 ± 0.010.549 ± 0.00
P43.71 ± 0.0516.70 ± 0.050.503 ± 0.00
P54.12 ± 0.0319.70 ± 0.000.521 ± 0.00
P64.15 ± 0.1116.70 ± 0.020.533 ± 0.01
P74.12 ± 0.0118.80 ± 0.020.526 ± 0.00
P84.43 ± 0.0320.50 ± 0.000.543 ± 0.00
P93.92 ± 0.0918.30 ± 0.100.547 ± 0.02
P103.86 ± 0.0519.10 ± 0.060.559 ± 0.01
P113.43 ± 0.0619.20 ± 0.080.538 ± 0.01
P124.25 ± 0.1016.00 ± 0.100.501 ± 0.00
P134.22 ± 0.1216.90 ± 0.090.525 ± 0.00
P143.53 ± 0.0915.10 ± 0.100.516 ± 0.00
P153.78 ± 0.0519.60 ± 0.000.532 ± 0.01
P163.54 ± 0.0720.70 ± 0.090.578 ± 0.00
P173.62 ± 0.1119.00 ± 0.050.537 ± 0.00
P183.66 ± 0.1017.60 ± 0.050.513 ± 0.00
P193.72 ± 0.0117.40 ± 0.800.558 ± 0.00
P203.52 ± 0.0816.60 ± 0.200.505 ± 0.00
Average value3.75 ± 0.0518.55 ± 0.050.535 ± 0.00
Median value3.82 ± 0.0618.40 ± 0.090.534 ± 0.00
T14.67 ± 0.0121.80 ± 0.070.594 ± 0.00
T23.44 ± 0.0321.70 ± 0.190.643 ± 0.00
T34.18 ± 0.0317.31 ± 0.110.572 ± 0.01
T44.53 ± 0.0021.33 ± 0.090.621 ± 0.00
T54.90 ± 0.0921.65 ± 0.030.629 ± 0.02
T64.16 ± 0.0117.86 ± 0.160.575 ± 0.00
T74.53 ± 0.0020.87 ± 0.120.633 ± 0.00
T83.50 ± 0.0417.14 ± 0.100.567 ± 0.01
T93.67 ± 0.0325.07 ± 0.160.656 ± 0.01
T103.73 ± 0.1124.21 ± 0.180.673 ± 0.00
T113.87 ± 0.0621.70 ± 0.220.608 ± 0.00
T123.82 ± 0.0826.20 ± 0.310.657 ± 0.01
T133.72 ± 0.1023.31 ± 0.240.647 ± 0.02
T143.83 ± 0.1422.20 ± 0.150.596 ± 0.00
T153.66 ± 0.2020.50 ± 0.110.592 ± 0.00
T163.69 ± 0.0817.40 ± 0.090.582 ± 0.00
T174.62 ± 0.0316.00 ± 0.090.577 ± 0.02
T184.17 ± 0.1017.60 ± 0.100.553 ± 0.00
T193.90 ± 0.1219.50 ± 0.040.602 ± 0.00
T203.64 ± 0.0126.41 ± 0.200.612 ± 0.01
Average value3.85 ± 0.0521.49 ± 0.110.605 ± 0.00
Median value4.01 ± 0.0620.98 ± 0.130.609 ± 0.00
P—Polish honey samples; T—Thai honey samples. SD-standard deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rosiak, E.; Madras-Majewska, B.; Teper, D.; Łepecka, A.; Zielińska, D. Cluster Analysis Classification of Honey from Two Different Climatic Zones Based on Selected Physicochemical and of Microbiological Parameters. Molecules 2021, 26, 2361. https://doi.org/10.3390/molecules26082361

AMA Style

Rosiak E, Madras-Majewska B, Teper D, Łepecka A, Zielińska D. Cluster Analysis Classification of Honey from Two Different Climatic Zones Based on Selected Physicochemical and of Microbiological Parameters. Molecules. 2021; 26(8):2361. https://doi.org/10.3390/molecules26082361

Chicago/Turabian Style

Rosiak, Elżbieta, Beata Madras-Majewska, Dariusz Teper, Anna Łepecka, and Dorota Zielińska. 2021. "Cluster Analysis Classification of Honey from Two Different Climatic Zones Based on Selected Physicochemical and of Microbiological Parameters" Molecules 26, no. 8: 2361. https://doi.org/10.3390/molecules26082361

APA Style

Rosiak, E., Madras-Majewska, B., Teper, D., Łepecka, A., & Zielińska, D. (2021). Cluster Analysis Classification of Honey from Two Different Climatic Zones Based on Selected Physicochemical and of Microbiological Parameters. Molecules, 26(8), 2361. https://doi.org/10.3390/molecules26082361

Article Metrics

Back to TopTop