Martin Kamler

Variation in honey bee gut microbial diversity affected by ontogenetic stage, age and geographic location

Social honey bees, Apis mellifera, host a set of distinct microbiota, which is similar across the continents and various honey bee species. Some of these bacteria, such as lactobacilli, have been linked to immunity and defence against pathogens. Pathogen defence is crucial, particularly in larval stages, as many pathogens affect the brood. However, information on larval microbiota is conflicting. Seven developmental stages and drones were sampled from 3 colonies at each of the 4 geographic locations of A. mellifera carnica, and the samples were maintained separately for analysis. We analysed the variation and abundance of important bacterial groups and taxa in the collected bees. Major bacterial groups were evaluated over the entire life of honey bee individuals, where digestive tracts of same aged bees were sampled in the course of time. The results showed that the microbial tract of 6-day-old 5th instar larvae were nearly equally rich in total microbial counts per total digestive tract weight as foraging bees, showing a high percentage of various lactobacilli (Firmicutes) and Gilliamella apicola (Gammaproteobacteria 1). However, during pupation, microbial counts were significantly reduced but recovered quickly by 6 days post-emergence. Between emergence and day 6, imago reached the highest counts of Firmicutes and Gammaproteobacteria, which then gradually declined with bee age. Redundancy analysis conducted using denaturing gradient gel electrophoresis identified bacterial species that were characteristic of each developmental stage. The results suggest that 3-day 4th instar larvae contain low microbial counts that increase 2-fold by day 6 and then decrease during pupation. Microbial succession of the imago begins soon after emergence. We found that bacterial counts do not show only yearly cycles within a colony, but vary on the individual level. Sampling and pooling adult bees or 6th day larvae may lead to high errors and variability, as both of these stages may be undergoing dynamic succession.

In-depth proteomic analysis of Varroa destructor: Detection of DWV-complex, ABPV, VdMLV and honeybee proteins in the mite.

We investigated pathogens in the parasitic honeybee mite Varroa destructor using nanoLC-MS/MS (TripleTOF) and 2D-E-MS/MS proteomics approaches supplemented with affinity-chromatography to concentrate trace target proteins. Peptides were detected from the currently uncharacterized Varroa destructor Macula-like virus (VdMLV), the deformed wing virus (DWV)-complex and the acute bee paralysis virus (ABPV). Peptide alignments revealed detection of complete structural DWV-complex block VP2-VP1-VP3, VDV-1 helicase and single-amino-acid substitution A/K/Q in VP1, the ABPV structural block VP1-VP4-VP2-VP3 including uncleaved VP4/VP2, and VdMLV coat protein. Isoforms of viral structural proteins of highest abundance were localized via 2D-E. The presence of all types of capsid/coat proteins of a particular virus suggested the presence of virions in Varroa. Also, matches between the MWs of viral structural proteins on 2D-E and their theoretical MWs indicated that viruses were not digested. The absence/scarce detection of non-structural proteins compared with high-abundance structural proteins suggest that the viruses did not replicate in the mite; hence, virions accumulate in the Varroa gut via hemolymph feeding. Hemolymph feeding also resulted in the detection of a variety of honeybee proteins. The advantages of MS-based proteomics for pathogen detection, false-positive pathogen detection, virus replication, posttranslational modifications, and the presence of honeybee proteins in Varroa are discussed.

Honeybee ( Apis mellifera )-associated bacterial community affected by American foulbrood: detection of Paenibacillus larvae via microbiome analysis

Honeybee (Apis mellifera L.) workers act as passive vectors of Paenibacillus larvae spores, which cause the quarantine disease American foulbrood (AFB). We assessed the relative proportions of P. larvae within the honeybee microbiome using metabarcoding analysis of the 16 S rRNA gene. The microbiome was analyzed in workers outside of the AFB zone (control - AFB0), in workers from asymptomatic colonies in an AFB apiary (AFB1), and in workers from colonies exhibiting clinical AFB symptoms (AFB2). The microbiome was processed for the entire community and for a cut-off microbiome comprising pathogenic/environmental bacteria following the removal of core bacterial sequences; varroosis levels were considered in the statistical analysis. No correlation was observed between AFB status and varroosis level, but AFB influenced the worker bee bacterial community, primarily the pathogenic/environmental bacteria. There was no significant difference in the relative abundance of P. larvae between the AFB1 and AFB0 colonies, but we did observe a 9-fold increase in P. larvae abundance in AFB2 relative to the abundance in AFB1. The relative sequence numbers of Citrobacter freundii and Hafnia alvei were higher in AFB2 and AFB1 than in AFB0, whereas Enterococcus faecalis, Klebsiella oxytoca, Spiroplasma melliferum and Morganella morganii were more abundant in AFB0 and AFB1 than in AFB2.

Comparison of Varroa destructor and Worker Honeybee Microbiota Within Hives Indicates Shared Bacteria

The ectoparasitic mite Varroa destructor is a major pest of the honeybee Apis mellifera. In a previous study, bacteria were found in the guts of mites collected from winter beehive debris and were identified using Sanger sequencing of their 16S rRNA genes. In this study, community comparison and diversity analyses were performed to examine the microbiota of honeybees and mites at the population level. The microbiota of the mites and honeybees in 26 colonies in seven apiaries in Czechia was studied. Between 10 and 50 Varroa females were collected from the bottom board, and 10 worker bees were removed from the peripheral comb of the same beehive. Both bees and mites were surface sterilized. Analysis of the 16S rRNA gene libraries revealed significant differences in the Varroa and honeybee microbiota. The Varroa microbiota was less diverse than was the honeybee microbiota, and the relative abundances of bacterial taxa in the mite and bee microbiota differed. The Varroa mites, but not the honeybees, were found to be inhabited by Diplorickettsia. The relative abundance of Arsenophonus, Morganella, Spiroplasma, Enterococcus, and Pseudomonas was higher in Varroa than in honeybees, and the Diplorickettsia symbiont detected in this study is specific to Varroa mites. The results demonstrated that there are shared bacteria between Varroa and honeybee populations but that these bacteria occur in different relative proportions in the honeybee and mite bacteriomes. These results support the suggestion of bacterial transfer via mites, although only some of the transferred bacteria may be harmful.

Bacterial community associated with worker honeybees (Apis mellifera) affected by European foulbrood

Background Melissococcus plutonius is an entomopathogenic bacterium that causes European foulbrood (EFB), a honeybee (Apis mellifera L.) disease that necessitates quarantine in some countries. In Czechia, positive evidence of EFB was absent for almost 40 years, until an outbreak in the Krkonose Mountains National Park in 2015. This occurrence of EFB gave us the opportunity to study the epizootiology of EFB by focusing on the microbiome of honeybee workers, which act as vectors of honeybee diseases within and between colonies. Methods The study included worker bees collected from brood combs of colonies (i) with no signs of EFB (EFB0), (ii) without clinical symptoms but located at an apiary showing clinical signs of EFB (EFB1), and (iii) with clinical symptoms of EFB (EFB2). In total, 49 samples from 27 honeybee colonies were included in the dataset evaluated in this study. Each biological sample consisted of 10 surface-sterilized worker bees processed for DNA extraction. All subjects were analyzed using conventional PCR and by metabarcoding analysis based on the 16S rRNA gene V1–V3 region, as performed through Illumina MiSeq amplicon sequencing. Results The bees from EFB2 colonies with clinical symptoms exhibited a 75-fold-higher incidence of M. plutonius than those from EFB1 asymptomatic colonies. Melissococcus plutonius was identified in all EFB1 colonies as well as in some of the control colonies. The proportions of Fructobacillus fructosus, Lactobacillus kunkeei, Gilliamella apicola, Frischella perrara, and Bifidobacterium coryneforme were higher in EFB2 than in EFB1, whereas Lactobacillus mellis was significantly higher in EFB2 than in EFB0. Snodgrassella alvi and L. melliventris, L. helsingborgensis and, L. kullabergensis exhibited higher proportion in EFB1 than in EFB2 and EFB0. The occurrence of Bartonella apis and Commensalibacter intestini were higher in EFB0 than in EFB2 and EFB1. Enterococcus faecalis incidence was highest in EFB2. Conclusions High-throughput Illumina sequencing permitted a semi-quantitative analysis of the presence of M. plutonius within the honeybee worker microbiome. The results of this study indicate that worker bees from EFB-diseased colonies are capable of transmitting M. plutonius due to the greatly increased incidence of the pathogen. The presence of M. plutonius sequences in control colonies supports the hypothesis that this pathogen exists in an enzootic state. The bacterial groups synergic to both the colonies with clinical signs of EFB and the EFB-asymptomatic colonies could be candidates for probiotics. This study confirms that E. faecalis is a secondary invader to M. plutonius; however, other putative secondary invaders were not identified in this study.

Spatio-temporal dynamics of Varroa destructor resistance to tau-fluvalinate in Czechia, associated with L925V sodium channel point mutation: Spread of tau-fluvalinate resistance in Varroa

BACKGROUND Extensive application of pyrethroids to control Varroa destructor, an invasive mite devastating bee colonies, has resulted in a global spread of resistant mite populations. In this study, we analyzed the spatio-temporal dynamics of resistant V. destructor populations in Czechia, stemming from the L925V mutation. Mites were collected during 2011-2018 directly or from winter beeswax debris, and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and densitometry was used to detect the L925V mutation. RESULTS Pooled samples of 10 mites were classified, based on their PCR-RFLP patterns, as tau-fluvalinate-sensitive (56%), resistant (9%), or mixed (35%), with the latter including sensitive and resistant homo- and heterozygotes. We identified two zones with higher frequencies of resistance, one in southern Moravia and the other in Bohemia. The mutant populations were evenly distributed throughout the monitored districts, with a few temporal and spatial local fluctuations. The greatest increase in resistance was observed in 2016, following massive losses of bee colonies in the winter of 2015. This event appeared to be closely associated with fluctuations in resistant mite populations and their dispersion. CONCLUSION Two outbreaks of resistance were detected in Czechia; however, the amount of applied tau-fluvalinate was not correlated with the frequency of resistance in mites. There was no remarkable increase in mite resistance in 2011-2018, although the use of tau-fluvalinate increased 40-fold between 2011 and 2015. PCR-RFLP analysis, performed on mites present in beeswax debris, is a suitable method for monitoring the L925V mutation in V. destructor.