Skip to main content

Beehives possess their own distinct microbiomes



Honeybees use plant material to manufacture their own food. These insect pollinators visit flowers repeatedly to collect nectar and pollen, which are shared with other hive bees to produce honey and beebread. While producing these products, beehives accumulate a considerable number of microbes, including bacteria that derive from plants and different parts of the honeybees’ body. Whether bacteria form similar communities amongst beehives, even if located in close proximity, is an ecologically important question that has been addressed in this study. Specific ecological factors such as the surrounding environment and the beekeeping methods used can shape the microbiome of the beehive as a whole, and eventually influence the health of the honeybees and their ecosystem.


We conducted 16S rRNA meta-taxonomic analysis on honey and beebread samples that were collected from 15 apiaries in the southeast of England to quantify the bacteria associated with different beehives. We observed that honeybee products carry a significant variety of bacterial groups that comprise bee commensals, environmental bacteria and symbionts and pathogens of plants and animals. Remarkably, this bacterial diversity differs not only amongst apiaries, but also between the beehives of the same apiary. In particular, the levels of the bee commensals varied significantly, and their fluctuations correlated with the presence of different environmental bacteria and various apiculture practices.


Our results show that every hive possesses their own distinct microbiome and that this very defined fingerprint is affected by multiple factors such as the nectar and pollen gathered from local plants, the management of the apiaries and the bacterial communities living around the beehives. Based on our findings, we suggest that the microbiome of beehives could be used as a valuable biosensor informing of the health of the honeybees and their surrounding environment.


Honeybees (Apis mellifera) use plant material to produce honey and beebread [1, 2]. Honey is made in the stomach of the adult workers, where the nectar collected from flowers is digested before regurgitation. Beebread is the collected pollen mixed with the young workers’ saliva. Both products are then further processed by microbes, including fermentative bacteria and yeasts that are thought to be involved in the crucial step of preservation [2,3,4,5]. Whether these microbes derive from the bees, or the environment is a very intriguing question that the scientific community has only addressed very recently [6, 7].

Recent studies have reported that the composition of the microbial community found in honey is dependent on the variety of floral nectars used by the bees [2, 8]. The nectar seems to contribute more significantly to species richness and microbial abundance than the honeybee gut [9, 10]. This microbial divergence is even more obvious in the beebread, where most of the microbes present in the pollen originates from the soil and phyllosphere [3, 11]. Furthermore, the microbes of nectar and pollen can be transferred by bees from plant to plant, from (or to) other insects and pollinators, and also shared with house bees within the same beehive, including beneficial bacteria [5, 12, 13] and pathogens [14, 15]. To date, we know that bees may accumulate a significant variety of different bacteria in their beehives; however, most of the studies have employed isolated beehives, primarily focusing on the microbiome of the pollen (or nectar) or the bee [3, 6, 8, 10]. No previous investigations have delved into the bacterial communities present in different beehives and how their fluctuations mirror the environment where honeybees forage and live.

In this study we have used samples of honey and beebread from beehives representing same or different apiaries to determine their bacterial profile. We postulate that the microbiome of beehives reflects the specific ecosystem where they are located, including the beekeeping methods used, and that honeybee products can be used as pooled samples to elucidate the bacterial species present in that ecosystem.


Experimental design, sampling, and DNA extraction

To test our hypothesis, we used 16S rDNA metataxonomy to characterise and compare the bacterial diversity present in samples of honey and beebread (referred as to pollen henceforward) that were collected from 15 apiaries in southeast England. Sample collection took place between mid-June and mid-August and targeted several habitats and soils in 4 different counties, as well as different beehives, some of which were located within the same apiary (same postcode) (see sampling in the Additional file 1: Excel File). Samples were collected directly from honeycomb frames using sterile swab tubes for pollen and containers for honey from which ten grams of each were put into sterile tubes that were immersed in liquid nitrogen and stored at − 80 °C. To ensure full representation of the whole beehive, samples were collected from different parts of the honeycomb. One hundred milligrams of the frozen samples were then ground manually using a mortar and pestle under sterile conditions and DNA was extracted using the BIO101 FastDNA® SPIN Kit for Soil in conjunction with a FastPrep® cell disrupter instrument (Bio101, ThermoSavant, Qbiogene) as we have previously reported [16]. We know that this kit results in enhanced extraction of DNA from both Gram-positive and Gram-negative bacteria, and therefore a realistic representation of complex microbial environments.

16S rRNA gene sequencing and taxonomy

DNA was quantified and quality-assured with a Thermo Scientific™ Nanodrop, and sequenced using the Ion Torrent PGM sequencer as previously described [17]. The sequencing process targeted the V1–V2 variable region of the bacterial 16S rRNA gene as the length of this region match with the coverage capacity of the sequencer. The V1–V2 amplicons were generated in triplicates, pooled to minimize the effect of PCR bias, and subjected to sequencing using the Ion PGM Template OT2 400 and Ion PGM Hi-Q Sequencing kits (Life Technologies Ltd, Paisley, UK). In total, we sequenced 39 samples from the 15 apiaries, including 24 honey and 15 pollen samples representing at least 1 beehive from each apiary and, in some cases, between 1 and 4 beehives from the same apiary (see sampling in the Additional file 1: Excel File). Therefore, sample IDs were designated with a number (1–15) to identify the apiary, followed by H or P and then A, B, C or D to indicate the product -honey or pollen- and the different hives of the apiary, respectively. Samples where A, B, C or D is not indicated correspond to apiaries with only 1 beehive. Following sequencing we used The CD-HIT-OTU pipeline to remove low quality sequences, pyrosequencing errors and chimeras [18], with the resulting sequences clustered into Operational Taxonomic Units (OTUs) at 97% identity. OTUs were then taxonomically classified down to genus rank against the Greengenes 16S rRNA gene database (13.5) using MOTHUR [19].

Statistical analysis

We performed calculations of alpha and beta diversity at genus level using the phyloseq Bioconductor package in R [20]. Alpha diversity was calculated to assess and compare the richness and evenness between apiaries and samples of honey and pollen based on the observed OTUs and the Inverse Simpson (InvSimpon) index. Different OTUS were converted to percentage of total reads and subjected to ANOVA with Tukey–Kramer post-hoc analysis for multiple comparisons with a confidence level of at least 95% (p < 0.05) using GenStat [21]. Results were illustrated using box and whisker plots. For the beta diversity calculations, we employed distance matrices of Jaccard and Bray–Curtis to identify compositional dissimilarities amongst apiaries based on the presence/absence and abundance of microbial communities, respectively. These dissimilarities were visualized on Principal Coordinates Analysis (PCoA) plots, and then analysed using permutation tests for Adonis (Permanova) and homogeneity of multivariate dispersions (Betadispersion) followed by Tukey–Kramer multiple comparisons with a confidence level of at least 95% (p < 0.05). Permanova and Betadispersion let us confirm significant differences in the composition and spread of the bacterial communities not only between apiaries but also between hives within the same apiary.

Taxonomical analysis

To complement the alpha and beta diversity analysis described above, we used ClustVis [22] to cluster the OTUs from all the beehives based on the type of sample (honey and pollen) and the soil habitat. The soil types were identified using the postcode on the Cranfield University Soil and Agrifood Institute Soilscapes tool ( Despite being in different locations (postcodes) some apiaries share the same habitat and potentially a similar ecosystem around them, introducing a supplementary variable that could help group samples. The resulting clusters were visualized at phylum, class and order level using a heatmap. The identification of the OTUs at species level was finally performed by carrying out a BLASTN search of the representative sequences of each OTU against the NR database. A match was considered significant if it had greater or equal to 98% sequence identity and 100% coverage of the query sequence. Species that passed these filters were then classified as either bee symbionts, invertebrate symbionts, vertebrate symbionts, environmental bacteria, or pathogens by reference to the NCBI scientific literature.

Results and discussion

Sequencing data and analysis

The sequencing of the V1-V2 amplicons from the 39 samples generated 2.2 million reads, with an average of 57,500 reads/sample and length of 300 bp, and a total of 90 potential OTUs (see OTU’s raw data in the Additional file 1: Excel File). The read coverage was nearly 100% for all samples and, on average, 36% (min < 1%; max 99%) of the reads derived from each sample were non-bacterial, generally representing matches to chloroplasts or mitochondria of plants. The details of the target and off target OTUs are also indicated in the Additional file 1: Excel file. Although rarefaction curves showed that the sequencing depths of the 39 samples were adequate to maintain all of them in our analyses (Additional file 2: Fig. S1), we decided to exclude the OTUs that were only classified as plant as well as those with a bacterial identity % less than 97 and a very low confidence at phylum level (highlighted in red within the OUT’s raw data). Three pollen samples (6P; 10PA and 15P) with fewer than 2000 reads from OTUs classified as bacteria but containing very high levels of off-target matches were also ruled out. This resulted in 74 OTUs and 36 samples with an average of 29,868 reads per sample (see OUT’s normalized data in the Additional file 1: Excel File), from which all OTU counts were scaled to the minimum sample size for normalization and subsequent comparative analysis. The sequences of the initial 39 samples are available in the NCBI database under bioproject number PRJEB45401.

Bacterial diversity in beehive samples

The variety of bacterial communities present in all samples was first estimated using the two following alpha diversity measures: observed OTUs and InvSimpson index (Fig. 1). When samples were grouped by the type of honeybee product, no significant differences were observed between honey and pollen samples, either in terms of richness or evenness (Fig. 1A). Similar results were obtained using the apiaries as a grouping factor (Fig. 1B), specially with regards to the abundance of OTUs. The number of OTUs identified in all samples was not significantly different amongst the apiaries. However, estimations with the InvSimpon index revealed that apiary 13 display the highest diversity. With the exception of apiaries 5, 11, 14 and 15, the relative abundance of OTUs found in apiary 13 is significantly higher than that of the remainder apiaries. This first evidence of difference between apiaries led us to investigate the dissimilarities between the type of bacterial communities present in the samples.

Fig. 1
figure 1

The observed-OTUs richness and inverse Simpson index of bacterial communities found in beehives and grouped by the type of sample (A) and different apiaries (B). Tukey–Kramer post-hoc analysis was carried out for multiple comparisons (*, p < 0.05)

To visualize the distance in the bacterial OTUs found in honey and pollen samples isolated from different apiaries, we generated Principal Coordinates Analysis (PCoA) plots based on Jaccard and Bray–Curtis dissimilarity metrics (Fig. 2). The plots show that samples did not group by product type (honey vs pollen) or apiary. In fact, the beta diversity Adonis test confirmed this different compositional dissimilarity between the apiaries with a confidence level of higher than 99% (p = 0.0056). In terms of distribution, the differences were even significantly higher (p < 0.001), especially when the OTUs of apiaries 3, 4 and 14 were pairwise compared with those of 5, 6, 9, 10 and 11 (p < 0.05). In this respect, the different location of the apiaries could explain these variabilities; however, it is worth highlighting that apiaries 3, 4 and 11 are in close proximity, with approximately 4 miles of distance between the 3 of them. On average, honeybees forage distances no longer than 5 miles [23]. Furthermore, when those apiaries represented by more than 2 beehives were compared (9, 10 and 11), we also observed significant differences in their bacterial communities (p < 0.05). Taken together, our data suggests that there might be no consistent bacterial fingerprint for beehives, even when samples of honey and pollen are taken from the same apiary. These divergences were confirmed following a comprehensive taxonomical analysis described below.

Fig. 2
figure 2

PCoA plots of the Jaccard (A) and Bray–Curtis (B) dissimilarities for bacterial communities (OTUs) found in samples of honey and pollen that were collected from the 15 different apiaries selected in this study. Samples with the same colour derive from the same apiary

Taxonomical analysis

The classification of the OTUs into different taxonomic ranks resulted in the identification of 5 different phyla across samples, with Firmicutes and Proteobacteria the most abundant, for which the number of different classes and orders detected were 7 and 19, respectively (Additional file 1: Excel file and Additional file 2: Fig. S2). Sequences from Gammaproteobacteria and Bacilli were very frequent, and within these classes, the orders Pseudomonadales, Enterobacterales and Lactobacillales were the most predominant. Of particular note is the fact that seven unclassified bacterial communities, potentially derived from animal faeces and tissues, were found in nearly all samples, with a very high relative proportion in sample 13H. Moreover, chloroplast sequences were also identified within honey and pollen samples, likely as a result of the bacterial origin of this organelle (Additional file 1: Excel file). Recent studies have reported that plant chloroplast sequences are very prevalent on honeybee products and can help to determine the foraging patterns of the bees [24]. In this study we observed matches to chloroplasts belonged to different plant species, including Adenophora stricta, Citrullus lanatus, Fagus sylvatica, Quercus fenchengensis, Raphanus sativus and Salix paraflabellaris (“Off target” OTUs in the Additional file 1: Excel file). However, it should be noted that 16S rDNA techniques do not give the best resolution for distinguishing different plants, with the RuBisCo large subunit (rbcL) and maturase K (matK) genes being better biomarkers [25].

The clustering of OTUs at phylum, class and order level across the different samples and the various soil type habitats of the apiaries gave us, once again, a very undefined bacterial profile amongst the beehives (Additional file 2: Fig. S2). Samples of honey and pollen did not group at all, and although more similarity was observed using the soil variable, especially when samples derive from the same apiary, instances of variability were confirmed amongst the beehive samples, with no clear clusters grouped by soil-type habitat. To delve into the reasons of this irregular microbiome structure, we carried out a further OTU analysis at lower taxonomical levels. This analysis revealed the presence of 40 genera (Additional file 1: Excel file), from which we identified 44 different species representing different bacterial communities as referred to their ecological features (Fig. 3). Although most of the samples were dominated by bee symbionts and environmental bacteria (Fig. 3), we also detected symbionts of invertebrates and vertebrates, as well as potential pathogens of plants and humans such as Enterococcus faecalis, Lonsdalea britannica, Pseudomonas syringae, Staphylococcus aureus, Xanthomonas campestris and Yersinia mollaretii [26,27,28]. Similar to the lack of microbiome consistency discussed above, the abundance and distribution of OTUs within the different bacterial communities varied amongst the different beehive samples, with no clear core microbiome defining honey and pollen and the different apiaries that the samples represent (Fig. 3 and Additional file 1: Excel file). The only few species that were detected in all samples of honey and pollen included the plant endophyte Cutibacterium acnes [29] and the two bee symbionts Lactobacillus kunkeei and Parasaccharibacter apium [30]. Furthermore, we observed some contradictions to the general agreement that pollen carries more environmental bacteria than bee commensals, when compared to honey, and vice versa [9, 11]. For instance, the bee symbionts Arsenophonus nasoniae, Gilliamella apicola, Lactobacillus apis, and Snodgrassella alvi [31] were more frequently found in our pollen samples, while the plant and water associated bacteria Lactococcus lactis [32], Pelomonas puraquae [33] and Pseudomonas graminis [34] showed a higher prevalence in honey. As expected, the number of honey samples populated by bacteria previously found in nectar, such as Acinetobacter boissieri [35] and Fructobacillus fructosus [36], was higher than that of pollen.

Fig. 3
figure 3

Bacterial communities of bee symbionts (orange), invertebrate symbionts (pink), vertebrate symbionts (yellow), environmental bacteria (green), and pathogens (red) found in honey (A) and pollen (B) samples. Other bacteria from which only the genus was identified are indicated in blue, while those unclassified are represented in grey. The species-level analysis of the identified OTUs revealed: (i) Bee symbionts (orange) isolated from honey, pollen and honeybees, including Arsenophonus nasoniae, Bartonella apis, Bombilactobacillus mellis, Frischella perrara, Gilliamella apicola, Lactobacillus kunkeei, L. helsingborgensis, L. apis, Parasaccharibacter apium, Snodgrassella alvi, and Spiroplasma melliferum; (ii) Invertebrate symbionts (pink) found in other insects and nematods, including Commensalibacter intestine, Moraxella osloensis, Photorhabdus kayaii and Serratia symbiotica; (iii) Vertebrate symbionts (yellow) found in the skin and gut of birds, mammals and humans, including Acinetobacter pullicarnis, Haemophilus parainfluenzae, Lactobacillus salivarius, and Microbacterium hominis; (iv) Environmental bacteria (green) found in water, soil, plants, seeds, fruits, food and animal faeces, some of which may cause infections in plants and animals, such as Acinetobacter boissieri, A. chinensis, A. junii, Bacillus thuringiensis, Brevundimonas diminuta, B. mediterranea, Burkholderia cepacia, Cutibacterium acnes, Fructobacillus fructosus, F. tropaeoli, Lactococcus lactis, Leuconostoc mesenteroides, Methyloversatilis discipulorum, Neokomagataea tanensis, Pantoea vagans, P. agglomerans, Pelomonas puraquae, Pseudomonas fluorescens, P. graminis, and Zymobacter palmae; (v) Pathogens (red) that cause diseases in plants, animals and humans, including Enterococcus faecalis, Lonsdalea britannica, Pseudomonas syringae, Staphylococcus aureus, Xanthomonas campestris, and Yersinia mollaretii; and (vi) other bacteria (blue) representing vertebrate symbionts and environmental bacteria, including Acinetobacter, Erwinia, Fibrobacter, Mycoplasma, Pantoea, Prevotella, Ralstonia, and Undibacterium

Bacterial community of bee symbionts

Finally, we investigated why the community of bee symbionts fluctuate amongst the beehives (Fig. 3) and whether those fluctuations may be dependent on the presence and abundance of certain environmental bacteria and the apiculture methods used by the beekeepers. To answer these questions, we took advantage of the sequencing data from the 24 honey samples to obtain the total percentage of bee symbionts and environmental bacteria present in each of the samples. We also calculated the relative proportion of the different species and families found within both groups of bacteria. As indicated in Table 1, there was a clear correlation between the abundance of symbiotic and environmental bacteria and the use of antimicrobials in the beehives. On the one hand, the microbiome of apiaries that did not receive any antimicrobial treatment (e.g. amitraz, fumidil B, oxalic/formic acid and/or thymol), has an overrepresentation of the very well know honeybee gut commensal L. kunkeei (samples from apiaries 1, 2, 8 and 12) and a low percentage of bacteria commonly found in plants, including species of the Erwiniacea (Pantoea agglomerans and P. vagans) and Pseudomonadaceae (Pseudomonas fluorescens and P. graminis). In contrast, beehives exposed to antimicrobial treatments due to previous infestations and/or infections (Chalkbrood, Nosema apis, sacbrood, Varroa, and/or wax moth), showed a much lower percentage of symbiotic bacteria, but a more diverse community comprising not only A. kunkeei but also Parasaccharibacter apium, Arsenophonus nasoniae, Serratia symbiotica and Gilliamella apicola (samples from apiaries 7, 9, 11, 14, 15). This reduced level of bee commensals shifted significantly in favour of some environmental bacteria belonging to the families of Propionibacteriaceae (Cutibacterium acnes), Lactobacillaceae (Fructobacillus fructosus, F. tropaeola, Lactobacillus salivarius, Leuconostoc mesenteroides) and Streptococcaceae (Lactococcus lactis). In beehives where antimicrobial treatments were used as a prophylactic measure, the abundance of bee symbionts ranged from high to moderate or even low (samples from apiaries 5, 10, 13), but the structure of the bacterial communities was similar to that of the beehives suffering from infestations and/or infections and subsequent antimicrobial therapies.

Table 1 The percentage of bee symbionts and environmental bacteria and their corresponding species and families that were detected in honey samples representing 24 beehives where different antimicrobials were used as a prophylactic measure or therapeutic treatment due to previous infestations or infections

Here, we have shown that changes in the microbiome of the beehive, in particular those affecting the bee symbionts, associate with specific apiculture methods. This finding contradicts a recent study reporting that different hive practices do not influence the bee gut microbiome [37]. However, we must consider that the authors of this study used preparations of whole bees and dissected bee samples instead of samples originated from beehives. On the other hand, our data indicate that the use of antimicrobials could be an underlying factor leading to more susceptibility to infestations or infections in the beehive. Previous studies have shown that some species of symbionts are more abundant in bees from pathogen-infected colonies [38, 39]. For instance, L. kunkeei antagonizes bee pathogens [40], and this is the main symbiotic species that we have found in antibiotic-free beehives. However, the eventual presence of other species such Bartonella apis, Frischella perrara, G. apicola, and Snodgrasella alvi, which have also been detected in our samples, may be crucial to maintaining the health of the colony [41]. Additionally, we have observed that the presence of certain environmental bacteria correlates with the levels of bee symbionts, suggesting the influence of microbes that, in principle, derive from plants, water or soil, in the conformation of the beehive microbiome. On this point, samples derived from apiaries 3, 4 and 6 are a clear example of how spore-forming bacteria, in particular Bacillus thuringiensis may outcompete symbiotic bacteria and overpopulate the beehive, as recently reported [42].


Our exploratory study shows that honeybee products carry a significant diversity of bacterial species, particularly from the bees, plants and the environment; and also that there is an inconsistent microbial pattern, not only between honey and pollen, but also among samples collected from the same apiary and/or habitat. In agreement with very recent studies, we have confirmed that the beehive microbiome is defined by multiple environmental and ecological factors, such as the soil, habitat, local plants and bee forage [43, 44]; and most importantly, our results suggest, for the first time, that every beehive possesses their own distinct bacteria. Furthermore, we have also described that the presence of certain environmental bacteria seems to influence the levels at which the most prevalent bee symbionts are detected. Whether apiaries have received any antimicrobial treatment may also affect the type of bacteria present in the beehives. We thus postulate that the use of antibiotics and the incidental existence of some environmental bacteria in the beehive defines the prevalence and abundance of honeybee gut symbionts. Depending on how beneficial these symbionts are, they might then promote or jeopardize the health of the colony.

In conclusion, we propose that the DNA present in honey and pollen could inform us of microbial changes indicative of the health of the beehive ecosystem, including not only social bees and plants but also solitary bees and other animals living within that ecosystem. Further ecological studies that include a comprehensive monitoring on the beehive microbiome throughout different seasons and years as well as additional sampling of potential sources of bacteria from the environment would be important to verify the origins and dynamics of the microbial communities present in the beehives.

Availability of data and materials

The sequencing data is available at under bio-project number PRJEB45401.


  1. Ball DW. The chemical composition of honey. J Chem Educ. 2007;84(10):1643–6.

    Article  CAS  Google Scholar 

  2. Anderson KE, Sheehan TH, Mott BM, Maes P, Snyder L, Schwan MR, et al. Microbial ecology of the hive and pollination landscape: bacterial associates from floral nectar, the alimentary tract and stored food of honey bees (Apis mellifera). PLoS ONE. 2013;8(12):e83125.

    Article  CAS  Google Scholar 

  3. Anderson KE, Carroll MJ, Sheehan T, Lanan MC, Mott BM, Maes P, et al. Hive-stored pollen of honey bees: many lines of evidence are consistent with pollen preservation, not nutrient conversion. Mol Ecol. 2014;23(23):5904–17.

    Article  CAS  Google Scholar 

  4. Olofsson TC, Vasquez A. Detection and identification of a novel lactic acid bacterial flora within the honey stomach of the honeybee Apis mellifera. Curr Microbiol. 2008;57(4):356–63.

    Article  CAS  Google Scholar 

  5. Vasquez A, Forsgren E, Fries I, Paxton RJ, Flaberg E, Szekely L, et al. Symbionts as major modulators of insect health: Lactic Acid Bacteria and honeybees. PLoS ONE. 2012;7(3):e33188.

    Article  CAS  Google Scholar 

  6. Martinson VG, Danforth BN, Minckley RL, Rueppell O, Tingek S, Moran NA. A simple and distinctive microbiota associated with honey bees and bumble bees. Mol Ecol. 2011;20(3):619–28.

    Article  Google Scholar 

  7. Saraiva MA, Zemolin AP, Franco JL, Boldo JT, Stefenon VM, Triplett EW, et al. Relationship between honeybee nutrition and their microbial communities. Antonie Van Leeuwenhoek. 2015;107(4):921–33.

    Article  CAS  Google Scholar 

  8. Sinacori M, Francesca N, Alfonzo A, Cruciata M, Sannino C, Settanni L, et al. Cultivable microorganisms associated with honeys of different geographical and botanical origin. Food Microbiol. 2014;38:284–94.

    Article  CAS  Google Scholar 

  9. Moran NA, Hansen AK, Powell JE, Sabree ZL. Distinctive gut microbiota of honey bees assessed using deep sampling from individual worker bees. PLoS ONE. 2012.

    Article  Google Scholar 

  10. Kwong WK, Moran NA. Gut microbial communities of social bees. Nat Rev Microbiol. 2016;14(6):374–84.

    Article  CAS  Google Scholar 

  11. Corby-Harris V, Maes P, Anderson KE. The bacterial communities associated with honey bee (Apis mellifera) foragers. PLoS ONE. 2014;9(4):e95056.

    Article  CAS  Google Scholar 

  12. Santos-Beneit F, Ceniceros A, Nikolaou A, Salas JA, Gutierrez-Merino J. Identification of antimicrobial compounds in two Streptomyces sp. strains isolated from beehives. Front Microbiol. 2022;13:742168.

    Article  Google Scholar 

  13. Grubbs KJ, May DS, Sardina JA, Dermenjian RK, Wyche TP, Pinto-Tomas AA, et al. Pollen Streptomyces produce antibiotic that inhibits the honey bee pathogen Paenibacillus larvae. Front Microbio. 2021.

    Article  Google Scholar 

  14. Card SD, Pearson MN, Clover GRG. Plant pathogens transmitted by pollen. Aust Plant Path. 2007;36(5):455–61.

    Article  Google Scholar 

  15. Gruter C, Moore H, Firmin N, Helantera H, Ratnieks FL. Flower constancy in honey bee workers (Apis mellifera) depends on ecologically realistic rewards. J Exp Biol. 2011;214(Pt 8):1397–402.

    Article  Google Scholar 

  16. Huws SA, Kim EJ, Cameron SJ, Girdwood SE, Davies L, Tweed J, et al. Characterization of the rumen lipidome and microbiome of steers fed a diet supplemented with flax and echium oil. Microb Biotechnol. 2015;8(2):331–41.

    Article  CAS  Google Scholar 

  17. Wilkinson TJ, Cowan AA, Vallin HE, Onime LA, Oyama LB, Cameron SJ, et al. Characterization of the microbiome along the gastrointestinal tract of growing turkeys. Front Microbiol. 2017;8:1089.

    Article  Google Scholar 

  18. Li W, Fu L, Niu B, Wu S, Wooley J. Ultrafast clustering algorithms for metagenomic sequence analysis. Brief Bioinform. 2012;13(6):656–68.

    Article  Google Scholar 

  19. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.

    Article  CAS  Google Scholar 

  20. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8(4):e61217.

    Article  CAS  Google Scholar 

  21. Hilbe JM. GenStat 9: a review. Am Stat. 2007;61(3):269–73.

    Article  Google Scholar 

  22. Metsalu T, Vilo J. ClustVis: a web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 2015;43(W1):W566–70.

    Article  CAS  Google Scholar 

  23. Hagler JR, Mueller S, Teuber LR, Machtley SA, Van Deynze A. Foraging range of honey bees, Apis mellifera, in alfalfa seed production fields. J Insect Sci. 2011;11:144.

    Article  Google Scholar 

  24. Jones L, Brennan GL, Lowe A, Creer S, Ford CR, de Vere N. Shifts in honeybee foraging reveal historical changes in floral resources. Commun Biol. 2021;4(1):37.

    Article  Google Scholar 

  25. Group CPW. A DNA barcode for land plants. Proc Natl Acad Sci U S A. 2009;106(31):12794–7.

    Article  Google Scholar 

  26. Mansfield J, Genin S, Magori S, Citovsky V, Sriariyanum M, Ronald P, et al. Top 10 plant pathogenic bacteria in molecular plant pathology. Mol Plant Pathol. 2012;13(6):614–29.

    Article  Google Scholar 

  27. Li Y, Xue H, Guo LM, Koltay A, Palacio-Bielsa A, Chang JP, et al. Elevation of three subspecies of Lonsdalea quercina to species level: Lonsdalea britannica sp nov., Lonsdalea iberica sp nov and Lonsdalea populi sp nov. Int J Syst Evol Micr. 2017;67(11):4680–4.

    Article  CAS  Google Scholar 

  28. Cantas L, Suer K. Review: the important bacterial zoonoses in “one health” concept. Front Public Health. 2014;2:144.

    Article  Google Scholar 

  29. Campisano A, Ometto L, Compant S, Pancher M, Antonielli L, Yousaf S, et al. Interkingdom transfer of the acne-causing agent, Propionibacterium acnes, from human to grapevine. Mol Biol Evol. 2014;31(5):1059–65.

    Article  CAS  Google Scholar 

  30. Drew GC, Budge GE, Frost CL, Neumann P, Siozios S, Yanez O, et al. Transitions in symbiosis: evidence for environmental acquisition and social transmission within a clade of heritable symbionts. ISME J. 2021;15(10):2956–68.

    Article  Google Scholar 

  31. Moran NA. Genomics of the honey bee microbiome. Curr Opin Insect Sci. 2015;10:22–8.

    Article  Google Scholar 

  32. Song AA, In LLA, Lim SHE, Rahim RA. A review on Lactococcus lactis: from food to factory. Microb Cell Fact. 2017;16(1):55.

    Article  CAS  Google Scholar 

  33. Gomila M, Bowien B, Falsen E, Moore ERB, Lalucat J. Description of Pelomonas aquatica sp. nov. and Pelomonas puraquae sp. nov., isolated from industrial and haemodialysis water. Int J Syst Evol Microbiol. 2007;57(Pt 11):2629–35.

    Article  CAS  Google Scholar 

  34. Behrendt U, Ulrich A, Schumann P, Erler W, Burghardt J, Seyfarth W. A taxonomic study of bacteria isolated from grasses: a proposed new species Pseudomonas graminis sp. nov. Int J Syst Bacteriol. 1999;49 Pt 1:297–308.

    Article  CAS  Google Scholar 

  35. Alvarez-Perez S, Lievens B, Jacquemyn H, Herrera CM. Acinetobacter nectaris sp. nov. and Acinetobacter boissieri sp. nov., isolated from floral nectar of wild Mediterranean insect-pollinated plants. Int J Syst Evol Microbiol. 2013;63(Pt 4):1532–9.

    Article  Google Scholar 

  36. Endo A, Maeno S, Tanizawa Y, Kneifel W, Arita M, Dicks L, et al. Fructophilic Lactic Acid Bacteria, a unique group of fructose-fermenting microbes. Appl Environ Microbiol. 2018.

    Article  Google Scholar 

  37. Subotic S, Boddicker AM, Nguyen VM, Rivers J, Briles CE, Mosier AC. Honey bee microbiome associated with different hive and sample types over a honey production season. PLoS ONE. 2019;14(11):e0223834.

    Article  CAS  Google Scholar 

  38. Rubanov A, Russell KA, Rothman JA, Nieh JC, McFrederick QS. Intensity of Nosema ceranae infection is associated with specific honey bee gut bacteria and weakly associated with gut microbiome structure. Sci Rep. 2019;9(1):3820.

    Article  CAS  Google Scholar 

  39. Floyd AS, Mott BM, Maes P, Copeland DC, McFrederick QS, Anderson KE. Microbial ecology of European Foul Brood disease in the honey bee (Apis mellifera): Towards a microbiome understanding of disease susceptibility. Insects. 2020.

    Article  Google Scholar 

  40. Arredondo D, Castelli L, Porrini MP, Garrido PM, Eguaras MJ, Zunino P, et al. Lactobacillus kunkeei strains decreased the infection by honey bee pathogens Paenibacillus larvae and Nosema ceranae. Benef Microbes. 2018;9(2):279–90.

    Article  CAS  Google Scholar 

  41. Raymann K, Moran NA. The role of the gut microbiome in health and disease of adult honey bee workers. Curr Opin Insect Sci. 2018;26:97–104.

    Article  Google Scholar 

  42. Kim DR, Cho G, Jeon CW, Weller DM, Thomashow LS, Paulitz TC, et al. A mutualistic interaction between Streptomyces bacteria, strawberry plants and pollinating bees. Nat Commun. 2019;10(1):4802.

    Article  CAS  Google Scholar 

  43. Donkersley P, Rhodes G, Pickup RW, Jones KC, Wilson K. Bacterial communities associated with honeybee food stores are correlated with land use. Ecol Evol. 2018;8(10):4743–56.

    Article  Google Scholar 

  44. Donkersley P, Rhodes G, Pickup RW, Jones KC, Power EF, Wright GA, et al. Nutritional composition of honey bee food stores vary with floral composition. Oecologia. 2017;185(4):749–61.

    Article  Google Scholar 

Download references


We would like to express our gratitude to the Surrey Beekeepers Association and beekeepers in Buckinghamshire, Essex and Surrey for helping us with sampling. We also thank Jaime Gonzalez Gutierrez de la Concha for his scientific advice and helpful discussions.


This study has been funded by University of Surrey start-up funds.

Author information

Authors and Affiliations



The first and corresponding authors LS and JGM planned and performed experiments, carried out data analysis and prepared and edited the manuscript. TW conducted technical experiments and contributed to data analysis and interpretation. RA and YR aided in experimentation and data analysis. CJC and SH designed experiments, helped with data interpretation and aided in preparing and editing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jorge Gutierrez-Merino.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors have read and approved the final manuscript.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Sampling data, OTU's raw data, OTU's normalized data, offtarget OTU's raw counts and target and offtarget OTU's data.

Additional file 2: Fig S1.

Rarefaction curves illustrating the sequencing depth of bacterial communities (identified OTUs) in the 39 samples of honey and pollen collected from the 15 apiaries of this study. Curves were generated in R using the package ggrare and plotted with ggplot2. Fig S2. Clustering of OTUs from our beehive samples provides a proportional representation of microbiomes at 3 taxonomical levels (phylum-class-order) across different samples of honey and pollen and soil type habitats. The IDs for the samples are indicated at the bottom, where the numbers designate the apiary (postcode location) and letters H or P denote the product type -honey or pollen, respectively. If H or P are followed by another letter (A-B-C-D) indicates different beehives within the same apiary. Different soils include: 5, Herb-rich chalk and limestone pastures, lime-rich deciduous woodlands; 6, Neutral and acid pastures and deciduous woodlands, acid communities such as bracken and gorse in the uplands; 7, Base-rich pastures and deciduous woodlands; 8, Wide range of pasture and woodland types; 14, Mostly lowland dry heath communities; 15, Mixed dry and wet lowland heath communities; 18, Grassland and arable some woodland; and 22: Arable grassland and woodland. The heatmap was generated using ClustVis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Santorelli, L.A., Wilkinson, T., Abdulmalik, R. et al. Beehives possess their own distinct microbiomes. Environmental Microbiome 18, 1 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: