1. Introduction
Despite its size, Uruguay has variations in climate1, geology2, soil3 and vegetation4. These variations would allow the production of different honeys. In addition, the country has a national system of protected areas (SNAP) where native vegetation occurs without competition from crops and their associated weeds, and where agrotoxic contamination is not expected. However, honey produced in Uruguay is exported without differentiating as a commodity. Nevertheless, honey is a food product that many countries want to regulate and also differentiate in its botanical or geographical origin. Honey is valorized by its phenolic compounds5, mineral content6, its content of particular sugars7 and also by the positive effect of its intestinal absorption8, and to prevent gastric affections9. Concerning the origin, honey is characterized by the botanical origin, searching to differentiate between floral or extra floral and among the floral, mono or multiflorals10)(11. Some honey such as Manuka honey has an added-value because it came from a native New Zealand myrtaceae (Leptospermum scoparium) and it is associated with medicinal properties12. Geographical origin13 enables to avoid honey coming from contaminated regions as much as to ensure consumers that honey comes from natural or protected areas5. To determine the honey origin many analytical procedures were developed concerning flavonoids or sugars profile, or trace element and mineral content14)(15)(16. Geographical origin was not specified in previous works about honey quality in Uruguay17)(18)(19. The importance of honey as an exportable product and the opportunity to enhance protected areas for this production make it necessary to investigate honey from these areas. In these areas the risk of appearance of contaminants such as glyphosate is minimized. In the last years glyphosate was detected and stopped honey exportation. In this sense, the objective of this work was to characterize Uruguayan honey from three protected areas and one not protected area throughout the year by physicochemical parameters, nutritional compounds and palynology, looking for possible markers of geographical or botanical origin. Also, glyphosate presence was evaluated.
2. Materials and methods
2.1 Location of apiaries
The locations of the apiaries chosen allowed to show both a latitude effect (higher temperature and rainfall in the north) and a longitude effect (similarity to Argentinean flora in the west or similarity to Brazilian flora in the east). Three of the selected sites (north, west and east) are Protected Areas belonging to the Uruguayan National System of Protected Areas (SNAP in Spanish), while one of them is not (south). The georeferentiation, soil types and location names are presented in Table 1. Local bees characterized by the predominance of hybrids were used20. The rainfall and temperature data of each production season are presented in Table 2. Oxalic acid was used against Varroa destructor control, except in the East.
(*)SNAP: Sistema Nacional de Áreas protegidas. (**) Soil names according United States Department of Agriculture (USDA). 1 Hives lost by Varroa destructor 2) Only a harvest was obtained due to Climatic adverse conditions.
2.2 Extraction of samples of honey
In each region 10 hives (Langstroth) were used, and unstretched commercial wax was placed in the center (frames 4 and 5) of the first box on the brood chamber. Two samples were taken from each hive. One hundred percent capped pieces of honeycomb were cut with disposable plastic knives. These combs were placed in labeled sterile jars and stored at -20°C. Honey was extracted from these combs into the laboratory. For this reason, hydroxymethyl furfural (HMF) was not measured21. Combs were extracted in two harvests, at spring-summer (December-February; harvest 1) and at summer-autumn (March-May; harvest 2). 66 samples were obtained (Table 1).
2.3 Palynological analysis
Twenty g of honey were dissolved with distilled water, according to Louveaux22. After centrifugation residual extract was observed with a microscope at 400x magnification. Six hundred pollen grains were counted to determine the relative abundance of each taxa23. Moreover, the presence of honeydew elements (HDE) was determined (e.g. algae, spores or hyphae fungi and microcrystals), as present or not.
2.4 Electrical conductivity (EC), pH and moisture
The EC and pH of 20% (w/v) honey solution was measured at 20ºC using a conductivity meter CON 10 (Oakton, USA) and a pHmeter JENWAY 3305 (U.K.). The honey solution was prepared in milli-Q® water (conductivity < 10 μS/cm). The results of conductivity were expressed as milli Siemens per centimeter (mS/cm). Moisture content was measured by refractometry according to the AOAC method24 with a manual refractometer ATAGO-MASTER 3 M (58° to 90° Brix), using a drop of honey of each hive. Data are expressed as moisture percentage (%).
2.5 Ca, Mg, Na and K content
The total Ca, Mg, Na and K were determined according to Paul25, using a flame atomic absorption spectrophotometer (PinAAcleTM 500, Perkin Elmer, USA). Briefly, 2 grams of honey were dissolved on a hot plate with a solution of HNO3 1M (Merck 65% p.p.a., distilled by sub boiling) and HCl 6 M (Merck, p.p.a.) in a Erlenmeyer under trap of steam during an hour at sub boiling (< 75°C), and were taken to volume of 25 ml with Milli-Q® water (18 MΩ.cm). Determinations were carried up with air-acetylene flame, (10-2.5, l/min). A multielement lamp, Lumina (Perkin Elmer, USA) was used for Ca (422.7 wavelength) and Mg (285.2), whereas Na and K were measured with emission (589 and 404.4 wavelength, respectively). For each analyte a standard curve was prepared from stock solutions of Ca, Mg, Na y K of 1000 mg/l (CertipurR, Merck, Germany). A blank (only the acids) was run with the samples. LQ was 1.386, 0.176, 0.110 and 7.826 mg/l for Ca, Mg, Na and k, respectively. Data were expressed in milligrams of each mineral by kilogram of honey (mg/kg).
2.6 Sugars content
Sugars were determined by high performance liquid chromatography (HPLC). One gram of honey of each sample was dissolved with 20 ml of deionized H2O and vortexed to homogenize. Solution was filtered previously to inject in the HPLC, with PVDF filters of 13 mm of diameter and 0.45µm of porous (Merck Millipore). A HPLC Prominence, LC-20A series, Shimadzu Corporation (Japan) was used, equipped with IR differential detector (RID-20A), a low-pressure gradient valve (LC-20AT), enabling 4-solvent gradient elution, auto sampler (SIL-20AC HT), oven (CTO-10AS VP) and degassing unit (DGU-20A 5R). Data was processed using LabSolutions software. A column Luna Omega SUGAR 100 Å, 3 µm, 250 x 4.6 mm (Phenomenex, USA) thermostatized at 40ºC was used. Mobile phase was Acetonitrile: H2O (80:20) with a flux of a 1ml.min-1. RID detector was thermostatized at 40ºC. Total time of each run was 30 min. To quantify sugars the following external standards from Sigma-Aldrich, USA, were used: D(-)-Fructose, D(+/-)- Sucrose (99.5%), D(+/-)-Galactose (≥99%), D(+/-)-Mannose (≥99%). For Glucose and isomaltulose, D(+/-) -Glucose anhydride AR® ACS (97-102%, Macron Fine Chemicals, USA) and Isomaltulose (94.5%, USP Reference Standard, USA) were utilized. Calibration curves were made for each sugar. Data were expressed as g of sugar by 100 g of honey (%). Limit of Quantification (LOQ) ≥ 0.01%.
2.7 Glyphosate content
Glyphosate determination was carried up with HPLC coupled with a mass spectrometer (LC-MS/MS), in the QSI Laboratory (Bremen Deutschland/Germany). For each region and harvest date samples were grouped to obtain 50 grams of honey. All regions and harvest dates had two groups, except for East 2 due to having only 7 hives to sample. Data were expressed as mg/kg. Limit of Detection was 0.005 mg/kg.
2.8 Statistical analysis
Data were analyzed by analysis of variance (ANOVA), comparing seven groups of honey (East 1 and 2, North 1 and 2, South 1 and 2, and West), and post hoc Tukey-Kramer test. The normality of residuals was adjusted by graphic methods (QQplot) and the Levene's test was used to test the homogeneity of the variance. Principal Component Analysis (PCA) was used to search the associations between physicochemical variables of the different honey and regions. The data was standardized. Electric conductivity, pH, sugar and minerals content were analyzed to group the different types of honey and the region considering two harvests. The relation between pollen content, Galactose and Mannose was determined by linear regression. Software R was used26.
3. Results
3.1 Palynology
It was possible to determine the family, genus or species of 37 types of pollen. Results are shown by region and harvest date in Tables 3, 4, 5, 6, 7, 8 and 9. The 66 honeys underwent a principal component analysis that is represented in Figure 1. West honeys show a great variation and are distributed in 3 quadrants. On the other hand, the other honeys are associated with one of the quadrants, varying with the harvest date. Eastern honeys are associated with a single quadrant regardless of the harvest date. Some regions present similar honeys despite the distance (South 1 and 2 and North 2). North has different honeys in different seasons.
3.2 Mineral content, conductivity, pH and moisture
In Table 10 mean values of calcium, potassium, magnesium, sodium and the sum of them are presented. All samples analyzed had a higher mineral content than previously reported by others27)(28)(29. In Table 11 mean values of conductivity, moisture and pH are presented. Despite floral origin, some honey samples showed high conductivity. The moisture values are within what is accepted30. The pH values are similar to other values found for Uruguayan honey17.
3.3 Sugar content
In Table 12 mean values of Glucose (G), Fructose (F), Sucrose (S), Isomaltulose (I), Mannose (M) and Galactose are presented. A quadratic relation (Figure 2) was obtained between the content of Mannose and the percentage of pollen from Echium plantagineum in the first harvest of the south (R2= 0.64, y = 0.029 + 0.155x+ 0.0068x²). A linear relation (Figure 3) was detected between Galactose content and percentage of pollen from Eucalyptus sp. in the west region (R2=0.81 y = 0.351 + 0.077 x).
Values are means ± SD; values within a column not sharing a common letter are significantly different (Tukey-Kramer < 0.05).
Values are means ± SD; values within a column not sharing a common letter are significantly different (Tukey-Kramer < 0.05.
3.4 Glyphosate content
Glyphosate was detected in south and west regions (Table 13). In the South, above the maximum limit of residual admissible by the European Union (0.05 mg/kg).
Values are means ± SD; values within a column not sharing a common letter are significantly different (Tukey-Kramer < 0.05)
4. Discussion
4.1 Palynology and PCA
The great dispersion in different quadrants of western honeys can be botanically explained. In the west we found 4 monofloral samples of Myrcianthes, 2 of Parkinsonia aculeata, 1 of Trifolium pratense. The remaining 3 are multifloral with different percentages of different pollen. The differences in vegetation and the soils where these plant species are found originate both groups of north honeys that could explain the differences between North 1 and 2; Baccharis is found in shallow soils (North 1) while Lithraea, Myrcianthes and Gleditsia in deep soils (North 2).
The two southern honeys are close in the PCA and present minor differences. The vegetation that originates is in the same soils but one is multifloral honey (South 1), while the other is monofloral Lotus honey (South 2). Some of the analyzed variables have a great impact on the results (e.g., Ca and Mg in East, Na in South 1). Species not previously cited from Uruguay honeys were found, like Parkinsonia aculeata, Tripodanthus acutifolius and Lithraea brasiliensis18)(31. This would be explained by the fact that samples were taken from different regions of the country in different seasons. Unlike other works17)(18)(19)(31, this one collected honey from different parts of the country for a year and directly from the hives. In addition, direct extraction from the honeycomb ensures that different or lesser pollen appears and is not diluted. These species generate monofloral honeys, which marks the preference of bees for them.
4.2 Moisture, pH and minerals
Values within world ranges were found for humidity, pH and sugars. High conductivity values (> 0.8 milliS / cm) were found for some honeys. These values could indicate an extrafloral origin. However, although some samples have HDE, the HDE/pollen ratio is less than 323. In addition, the sum of fructose plus glucose is greater than 60%, complying with the European standard for floral honeys. In the region (Argentina) other works report honeys with high conductivity, although they relate it to honeydews32. The reported values correspond to other studies from Uruguay17)(19. Also, high mineral content was found in the sum of the four minerals analyzed. Potassium is the dominant mineral as in other honeys from Apis6)(27)(28) or from Meliponas33.This mineral exceeds the average of 742 mg / kg cited28. The soils of Uruguay present variable potassium contents associated with the content of clays and minerals that originate the soil34. Differences in Ca values were found in the east more than in other regions, higher than the average calcium found in Brazil6, and in other countries of the world29. In Europe, only 4 values higher than the average for the eastern region are reported28. Also in the east are high values of Mg, although Mg is reported as a mineral with greater variation in its values29. The soils and waters of this region have high Ca and Mg contents explained by geology35. In the south, a higher sodium content was observed on both dates, without statistical differences. In contrast to the Ca, Na would not be a mineral that allows discrimination origin16, but in Uruguay, proximity to the sea and the absence of barriers to the sea winds could explain the values found in the south apiary. Sodium could be considered a “pollutant”, which comes through the air as heavy metals or pesticide molecules36. This could be valid for the entire coastal region of the Río de la Plata and the Atlantic Ocean.
4.3 Sugars
Sugar values are among the range reported in many studies around the world37)(38)(39. In all honeys the sum of F+G is more than 60%. This supports the floral origin of the honeys despite the presence of HDE. Eastern region shows higher fructose values than the other regions, but within the world ranges40)(41. The relationship found between Mannose and Equium R2 = 0.64 in the south was not found in the other group, where there is Equium pollen (North 2). This could be explained by the different behavior of the species on different soils10. In the case of the relationship between Galactose and Eucalyptus R2 = 0.81, the Eucalyptus species level should be reached to propose its use as a marker of botanical origin. Both of them are introduced species, Echium associated with crops and Eucalyptus in commercial forestry. Indeed, both species could not be associated with a specific region.
4.4 Glyphosate
The presence of glyphosate, even in protected areas, is explained by human activities. Agriculture in the south and west appears to be responsible for the pollution levels. Honey samples with Glycine max, Lotus, Trifolium pratense and Trifolium repens pollen are evidence of this intensive production of forage and crops in the south. Also, the presence of pollen from weeds associated with crops such as Senecio, Cirsium, Echium and Ammi, that are controlled with glyphosate. The levels found in South 2 would be explained by the accumulation of applications during the summer in soybeans, and by the killing of old pastures to implant winter crops. Low levels below the maximum admissible were also detected in the west protected area. Honey samples with Lotus and Trifolium pratense pollen are evidence of this intensive production of forage and crops in the west. There were soybean crops in the western region, the absence of pollen in the samples is due to the harvest date prior to flowering. A later harvest of summer honey could have a higher accumulation of glyphosate. The protected areas, with a limited surface, do not prevent bees from going out of bounds and work in agricultural areas. To avoid this problem, the location of the apiaries inside the protected area is crucial. A right management is necessary in this sense. It is interesting to note that in east and north protected areas, located near fields of pastoral animal production (cattle and sheep) and forest activities, this contaminant was not detected.
5. Conclusions
Different honeys were found in different regions and at different seasons. Values variables are between the ranges reported in the world. However, the high value of conductivity (> 0.8 mS/cm) is not explained by the presence of honey dew. The high mineral content of honeys must continue to be investigated since it can be an element of the soils of Uruguay. Different monofloral honeys of unreported species were found. Calcium could be used as a marker for the east region associated with soil type. In the same way sodium could be a marker for the region of the maritime coast but more honeys should be studied. Mannose was associated with Echium plantagineum in southern samples and could be suggested as a botanic marker in this region. More studies should be carried out to determine the species of eucalyptus that correlates with galactose. Further investigations are necessary in order to relate other sugars in honey with native species. Honey production in protected areas is not enough to ensure the absence of glyphosate.