Abstract

The phenomenon of antibiotic resistance stands as a paramount health challenge in the contemporary era. Within a One Health approach, it becomes crucial to effectively track the dissemination of antibiotic resistance, not only within humans and animals but also within the environment. To investigate the environment, the honey bee (Apis mellifera) has emerged as a prominent environmental bioindicator due to its social, behavioral, and morphological features. The objective of this study was to describe the antimicrobial resistance (AMR) patterns of bacterial isolates from the body surface and the gut of honey bees sampled from 33 colonies throughout the Emilia-Romagna region (Italy). A total of 608 strains were examined for 19 distinct antimicrobial compounds from various classes, and the results showed that more than 50% of the isolates for eight out of nine provinces showed characteristics of nonsusceptibility toward amoxicillin and penicillin, and, generally, 98.19% of isolated strains were considered AMR and 74.67% exhibited multidrug resistance (MDR) characteristics, more frequent in Gram-negative strains (87.74%) than in Gram-positive ones (60.34%). Additionally, a significant correlation with a lower prevalence of MDR bacteria was demonstrated for one province (Ferrara, odds ratio (OR) = 3.33, (1.67; 6.64), ). In conclusion, this study provides evidence for the utility of A. mellifera colonies as bioindicators for MDR bacteria, enabling their characterization and distribution at a geographical level. Additional investigations are required to further explore the potential role of honey bees as bioindicators for antimicrobial-resistant bacteria, particularly in terms of their association with environmental characteristics.

1. Introduction

Antimicrobial resistance (AMR) is a phenomenon born with antimicrobial drug discovery and use, whose exponential growth is racing with improper antimicrobial use both in human and veterinary medicine [1, 2]. The main institutions have established surveillance systems for monitoring the AMR spread in a distinct pathogen bacteria group; in Europe, it is represented by the European Antimicrobial Resistance Surveillance Network (EARS-Net) (https://www.ecdc.europa.eu/en/about-us/networks/disease-networks-and-laboratory-networks/ears-net-data, accessed June 2023). However, all these monitoring systems are limited in the sanitary contest without considering external factors, and within a One Health approach, it can be necessary to evaluate the possible spread of AMR and multidrug resistance (MDR) bacteria in the environment. To date, some studies have found AMR and antimicrobial resistance genes (ARGs) in different environments, such as rivers, embankments, and ponds, probably due to fecal contamination [3, 4]. Farming and husbandry areas can also be a potential reservoir of AMR and ARGs because soil bacteria can be influenced by animal manure derived from animals treated with antimicrobics [57]. Urban contexts may also be a possible source of transmission of AMR or ARGs when considering water treatment plants [8] or sanitary structures [9]. For this reason, the AMR may be regarded as an environmental factor in biomonitoring. Honey bee (Apis mellifera L.) colonies are commonly used as a bioindicator, thanks to the individual morphology and behavioral colony characteristics. Individually, they present a body surface covered by bristles and hairs, whereby they catch the pollen and intercept other types of particles during the flight. Concerning the behavior, they present a huge member population and an elevated number of flights of foragers per day, with a flight radius of around 1.5 km [1012]. Considering each bee as a microsampler, the total amount of foragers bees are able to perform a representative sample of the explored environment [11]. Besides, the honey bee colonies can be managed in any type of environment, including marginal, rural, agricultural, urban, and industrialized areas [10]. All these characteristics have made honey bee colonies useful for several environmental biomonitoring plans through the analysis of bee products and/or their bodies [1316] to detect environmental contaminants, such as heavy metals, pesticides, and microplastics [15, 1725], environmental pollutant, such as atmospheric particulate matters [2628], and, lately, plants, animals, human pathogens [15, 2933], AMR and ARGs [3439]. Honey bees can acquire AMR bacteria and ARGs during the foraging flight, especially during the ingestion of nectar and water. Besides, bacteria on flowers can be picked up by the fimbriae on their bodies during flight activity [10, 11, 27].

This study aims to describe the AMR patterns of bacterial strains isolated from the body surface and the gastrointestinal tract of honey bees, sampled from 33 colonies throughout the Emilia-Romagna region (Italy).

2. Materials and Methods

2.1. Sampling

The study was conducted as a part of the BeeNet project [40], including a network of 33 apiaries located in all provinces of the Emilia-Romagna (Supplementary 1). The sampling was performed in four different periods of the year, namely November 2021, March 2022, June 2022, and September 2022. All apiaries were sampled during each monthly sampling period. For each apiary, the province of origin was considered. Three hives were chosen for each apiary, and 10 returning foragers have been sampled for each hive [41, 42]. The 10 foragers have been placed in sterile 25 ml tubes type falcon and stored at 4 ± 3°C until analysis.

2.2. Bacterial Isolation and Identification

The bacterial isolation was made from both the body surface and the gastrointestinal tract. For bacterial isolation from the body surface, the 10 foragers belonging to the same colony were put in a 25 ml falcon tube with 6 ml of sterile physiologic solution and vortexed for 10 s [43, 44]. The consequent suspension was streaked by a 10 μl sterile disposable loop onto UTIC (Condalab, Madrid, Spain) agar medium and incubated aerobically at 37 ± 1°C for 24 hr; for the gastrointestinal tract’s bacterial isolation, all the 10 foragers used were dissected and for each forager, the ventriculus (small intestine and midgut) was extracted and put in 2 ml sterile microtubes with 1.5 ml of a sterile physiologic solution to create a pool with all the 10 guts. Subsequently, the pool was ground with a microbiological micro pestle and then streaked onto the UTIC medium as described above and incubated equally. After the incubation, the bacterial growth was evaluated, and the different colonies with different colors were isolated in purity onto the UTIC agar medium [45]. For the identification of isolated bacteria, each bacterial colony was subcultured in Tryptone Soy Agar (TSA) (Oxoid, Basington, UK) medium. The fresh colonies were identified using Matrix-Assisted Laser Desorption/Ionization Time-Of-Flight Mass Spectrometry (MALDI-TOF MS) (Biotyper, Bruker Inc., USA). All the colonies were identified at level species with a minimal score of 1.80 using Bruker Biotyper version 3.0 software.

2.3. Antimicrobial Susceptibility Testing (AST)

All the identified strains were tested for antimicrobial susceptibility to 18 antimicrobials according to the agar diffusion method described by the Clinical and Laboratory Standard Institute [46]. The 19 antimicrobials tested (see Table 1) are amoxicillin (25 μg), amoxicillin/clavulanic acid (20/10 μg), ampicillin (10 μg), aztreonam (30 μg), cefotaxime (30 μg), cefoxitin (30 μg), cephalothin (30 μg), chloramphenicol (30 μg), doxycycline (30 μg), enrofloxacin (5 μg), erythromycin (15 μg), gentamicin (10 μg), imipenem (10 μg), nalidixic acid (30 μg), penicillin (10 U), streptomycin (10 μg), tetracycline (30 μg), trimethoprim/sulfamethoxazole (1.25/23.75 μg), and vancomycin (30 μg). The antimicrobial susceptibility was evaluated according to clinical breakpoints provided by the European Committee on Antimicrobial Susceptibility Testing (EUCAST) or CSLI [46, 47]; the clinical breakpoints used were reported in Table 1. According to the National Reference Laboratory for Antimicrobial Resistance (Reg. 882/2004/CE) “Istituto Zooprofilattico Sperimentale del Lazio e della Toscana (IZS),” the intrinsic resistance of each bacterial species was evaluated and excluded from results (https://www.izslt.it/crab/wp-content/uploads/sites/8/2018/08/Tabelle-Resistenze-intrinseche-in-batteri-di-interesse-veterinario.pdf, accessed June 2023). For AST interpretation, the strains were divided into “susceptible” and “non-susceptible,” where the “non-susceptible” category included resistant and intermediate isolates. Following the definition given by Magiorakos et al. [48], isolates that were nonsusceptible to at least one antimicrobial drug were considered AMR, and isolates that were not susceptible to at least one antimicrobial drug from at least three different antimicrobial classes were considered MDR.

2.4. Statistical Analysis

For each isolate, data about the province of sampling (BO, RA, FC, FE, RN, MO, PC, PR, and RE) were extrapolated. Nonsusceptibility percentages for each tested antimicrobial were calculated by dividing the number of nonsusceptible strains by the number of tested strains. AMR and MDR percentages were evaluated by dividing the number of AMR/MDR strains by the number of total strains. The association between the nonsusceptibility percentages and province was tested with the Chi-squared test. The alpha risk was set to 0.05. A multivariate logistic regression was performed to assess the relation between AMR/MDR and province. Data were checked for multicollinearity with the Belsley–Kuh–Welsch technique. Heteroskedasticity and normality of residuals were assessed respectively by the White test and the Shapiro–Wilk test. A p-value < 0.05 was considered statistically significant. For the evaluation of the correlation between antimicrobials found to be statistically significant distributed by province and the livestock density, the Shapiro–Wilk tests were used. Pearson’s coefficient was used to assess the correlation, judged very strong from 1 to 0.9, strong from 0.9 to 0.7, moderate from 0.7 to 0.5, low from 0.5 to 0.3, and poor from 0.3 to 0. The alpha risk was set to 0.05. Statistical analysis was performed with EasyMedStat (version 3.24; https://www.easymedstat.com). Bovine and swine populations were the only animal species considered, since they represent the most frequent type of livestock in the considered area (National Database of Italian Zootechnical Registry (BDN) https://www.vetinfo.it/j6_statistiche/index.html#/, accessed June 2023).

3. Results

3.1. Bacteria Isolation and Identification

A total of 608 strains, belonging to 29 genera and 79 species, were isolated. In detail, 84 strains (13.82%) were isolated from the body surfaces of foragers, while 524 (86.18%) were from the gut. Three-hundred-eight strains (52.3%) were classified as Gram-negative bacteria, while 290 (47.7%) as Gram-positive. All genera and species identified are reported in Supplementary 2. Considering the sampling period, a total of 97, 138, 207, and 166 bacteria were isolated, respectively, in November 2021, March 2022, June 2022, and September 2022, as shown in Figure 1.

The higher percentage of the strain isolated belonged to the genus Bacillus (27%), isolated in all provinces; the other two bacteria genera isolated in all provinces were Klebsiella spp., Pantoea spp., and Enterobacter spp., which represented 8.39%, 8.22%, and 15.95% of the total bacteria isolated, respectively (Figure 2).

3.2. Antimicrobial Susceptibility Testing and Multivariate Analysis

AST results and percentages for each tested drug of the total number of isolated strains and other categories (Gram-positive, Gram-negative, isolation from the gut or external surface) are shown in Table 2. Five hundred ninety-seven isolates (98.19%) were considered AMR, while 454 (74.67%) were considered MDR (Figure 3).

The highest percentages of nonsusceptibility strains resulted in amoxicillin and penicillin with percentages of 62.99% and 62.34%, respectively, followed by erythromycin (59.11%), ampicillin (54.95%), and aztreonam (54.82%); these results were in line with the categories considered (Gram-positive, Gram-negative, isolation from the gut or external surface) except for Gram-positives, where the only nonsusceptibility percentages above 50% were for cefotaxime and aztreonam (Table 2).

The nonsusceptibility percentages distribution of isolates by sampling province are shown in Table 3. Among all the molecules tested, amoxicillin and penicillin resulted in higher percentages of nonsusceptibility (> 50% of strains isolated) in all the provinces but Ferrara, where the corresponding percentages were 48.75% for amoxicillin and 46.15% for penicillin. Among all provinces, the highest percentage for amoxicillin was in Modena (71.79%); instead, the highest percentage for penicillin was in Reggio Emilia (69.39%). Regarding ampicillin and aztreonam, the nonsusceptibility percentages of isolates were above 50% in six out of nine provinces (except Ferrara, Modena, and Reggio Emilia for aztreonam and Ravenna, Forlì-Cesena, and Ferrara for ampicillin). Considering erythromycin, the strains presented nonsusceptibility percentages above 50% in seven provinces except for Ferrara and Rimini, and the highest value was recorded in the Modena province, with 84.21% of isolates found to be nonsusceptible (Table 3). In the statistical analysis between the province and the single-drug nonsusceptibility, the p-value was considered statistically significant for FOX (), AX (), CTX (), DO (), TE (), and KF (). In multivariate analysis, FE province (odds ratio (OR) = 3.33, (1.67; 6.64), ) was associated with lower percentages of MDR.

A strong positive correlation was found between TE nonsusceptibility percentages and bovine density (ρ = 0.74; r2 = 0.542; ) and between DO percentages and swine density (ρ = 0.72; r2 = 0.517; ). A strong positive correlation was found between the percentage of nonsusceptibility to at least one of the two tested tetracyclines and swine (ρ = 0.77; r2 = 0.591; ) and bovine density (ρ = 0.8; r2 = 0.634; ). These results are shown in Figure 4.

4. Discussion

Honey bees are widely considered excellent bioindicators due to their distinctive characteristics [13, 30, 49]. Notably, for AMR, bees have demonstrated their potential as specific indicators [34, 50], providing insights into the epidemiological distribution of ARGs.

In this study, 608 bacterial strains were isolated from bees collected from apiaries distributed across all nine Emilia-Romagna region provinces. These strains were subsequently identified; some identified species, such as Escherichia coli, Klebsiella pneumoniae, Enterococcus faecium, and Staphylococcus aureus, are widely recognized as commensal organisms but also have the potential to be pathogenic and cause infectious diseases in both animals and humans [5155]. These species are included in the list of “priority pathogens” redacted by the World Health Organization [56], namely bacteria for which new antimicrobials are urgently needed. The possibility of comparing their susceptibility to antimicrobials between the clinical and environmental settings, using bee colonies as bioindicators, may provide insights into the indirect effects of antimicrobial use in healthcare on the environment. Furthermore, it could potentially serve as an indication for improving the management of antimicrobial treatments for specific bacterial species of interest. The results of this investigation indicated also that isolated strains genera and species change during the sampling period. According to previous studies, the bee microbiota and their colonizing bacteria strains change seasonally, probably due to different ontogenic and climate conditions, including diets and environments, which influence the growth percentages of bacteria [31, 33, 43, 44, 57, 58].

All isolated strains were tested for their susceptibility to different classes of antimicrobials, and the results revealed alarmingly high percentages of nonsusceptibility, which can be considered indicative of the regional AMR situation.

The results indicate that, in nearly all of the different antimicrobials considered, amoxicillin and penicillin exhibited higher percentages of nonsusceptibility among the isolated strains. This finding can be attributed to the fact that the penicillin class is extensively used in Italy, both in human medicine and, based on preliminary data from the Italian Drug Association, in animal husbandry as well (https://www.aifa.gov.it/documents/20142/1853258, accessed June 2023). Due to the incompleteness of these data, making an accurate comparison is challenging. Nevertheless, it still provides some insight into the use of honey bee colonies as indicators for the spread of antimicrobial-resistant bacteria to specific antimicrobials.

The resistance percentages observed against antimicrobials classified as the highest priority for human medicine by the World Health Organization are of particular concern (https://www.who.int/publications/i/item/9789241515528, accessed June 2023), such as quinolones, carbapenems, penicillins with beta-lactamase inhibitors, macrolides, and others [59]. Given that the use of antibiotics is prohibited for bees in Europe (https://eur-lex.europa.eu/legal-content/IT/ALL/?uri=celex:32004L0028, accessed June 2023), it can be excluded the possibility that the isolated bacteria have experienced selective pressure from direct drug treatments, enabling them to acquire nonsusceptible traits. Instead, these traits likely stem from surrounding environmental contamination, particularly through resources such as water and pollen. Notably, no significant differences were observed in terms of genus/species variability and AMR patterns between samples collected from the external surface and the gut, suggesting that both sites serve as valid indicators. It is worth noting that the gut provides a more favorable physiological environment for bacterial growth and colonization compared to the external surface, which was reflected in a higher number of isolates obtained from the gastrointestinal tract [49].

Generally, Gram-negative bacteria exhibited higher resistance percentages, likely attributed to their propensity for acquiring mobile resistance elements, such as plasmids, as compared to Gram-positive bacteria [60]. However, this evaluation is based on phenotypic qualitative assessment, and further genotypic analysis is necessary to confirm this hypothesis, and it should be acknowledged as a limitation of the study.

Another important aspect to consider is the variation in resistance percentages across different provinces, which underscores the potential of honey bees as bioindicators for detecting geographic differences. Although the limited number of apiaries sampled does not allow to conclude that the differences between provinces are truly representative of the real AMR situation, it can be speculated that the results from each apiary somehow reflect the area in which they are located and their sum can be considered a partial representation of the AMR situation in the province. Interestingly, one particular province (FE) exhibited statistically significant lower percentages of MDR bacteria. It is currently impossible to determine the exact origin of this specific situation, and we are unable to formulate a specific hypothesis regarding its cause. Analyzing the land use and weather conditions during the sampling period could potentially provide valuable insights and help identify explanatory variables. Indeed, previous studies have demonstrated that the urban environment, characterized by human activities, serves as a reservoir for AMR due to factors like water disposal systems and pollution [9, 61, 62]. Additionally, specific molecules, including cephalosporins (CTX, FOX, KF), tetracyclines (DO, TE), and penicillins (AX), were found to be statistically associated with provincial disparities. Once again, it is not possible to determine the exact reason for this disparity. A possibility could be to conduct a cross-analysis that might reveal a correlation between the total consumption of antimicrobials (including human and animal use) in the different provinces and the results obtained in this study. However, such data are not publicly available, and the only aspect we can emphasize is the fact that the utilization of bee colonies highlights a geographical disparity for certain antimicrobials and that these findings are important in evaluating their potential use as bioindicators. The honey bees have long been recognized as reservoirs for various tetracycline resistance genes [63, 64]. Furthermore, as shown in Figure 4 and the statistical analysis, the comparison of geographical distributions of resistance to TE and DO reveals some correlation with the density of swine and bovine farms in the respective areas (measured in animals per km2). This correlation has to be considered in relation to the use of this antimicrobial class in livestock. In Italy, tetracyclines are frequently used in livestock animals (they represent 23.2% of total antimicrobials sales in 2021) (https://www.salute.gov.it/imgs/C_17_pubblicazioni_3281_allegato.pdf, accessed June 2023) in both swine and bovine farms (19% and 28%, respectively). Although data about tetracyclines consumption in livestock for each province were not available, it can be speculated that the antimicrobial use is proportional to the animal density for each province and that the nonsusceptibility percentages toward tetracyclines found in bees could reflect their use in farms. Considering cephalosporins, the distribution among provinces was found to have statistical significance, but no strong correlation with bovine or swine density was identified. This lack of correlation could be attributed to the limited usage of cephalosporins in livestock, accounting for only 0.2% of total sales in 2021. As a result, it is not possible to establish a direct connection with animal husbandry. However, the influence of urban and industrial settings with high human activities can be hypothesized as a contributing factor.

In general, it can be argued that honey bees act as location-dependent bioindicators, reflecting the selective pressure exerted by antibiotics in localized areas. This phenomenon is attributed to the foraging and drinking activities of honey bees, as well as the horizontal transfer of resistance genes within their microbiota [36, 65].

5. Conclusions

This study indicates another use of honey bee colonies as bioindicators for human, animal, and environmental health. The observed high percentages of nonsusceptibility highlight the importance of considering this insect species within a One Health approach. Honey bee colonies can be used as a readily accessible and sensitive indicator for evaluating the geographic distribution of certain antimicrobic resistances. Specifically, the statistical association between MDR percentages and province suggests a potential use in terms of AMR location-dependent bioindicators, especially for some antimicrobials such as tetracyclines and cephalosporins, for which honey bees seemed to show a specific sensitivity. The potential link between AMR percentages in honey bees and the antimicrobial use in both humans and animals and its use for surveillance purposes still needs to be better investigated, especially in the broader context of a One Health perspective.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Conceptualization was done by G.C.; sampling was made by G.C., I.R., and S.A.; methodology was done by G.C. and S.P.; investigation was done by G.C., I.R., R.S., and S.P.; data analysis was done by G.C., I.R., R.S., L.Z., S.A., L.B., A.N., and S.P.; writing—original draft preparation was done by G.C., I.R., R.S., and S.P.; writing—review and editing was done by G.C., L.Z., S.A., L.B., A.N., and S.P.; funding acquisition was made by L.B. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This study was supported by the project BeeNet (Italian National Fund under FEASR 2014–2020) from the Italian Ministry of Agriculture, Food Sovereignty and Forestry (MASAF). The authors are grateful to Amanda Dettori, Irene Guerra, Valeria Caringi and Vittorio Capano (CREA Research Centre for Agriculture and Environment) for their effort in the sampling activity.

Supplementary Materials

Supplementary 1. Geographical characteristics of investigated apiaries.

Supplementary 2. List of bacteria species isolated from honey bees’ surface and ventriculus samples.