Solitary bees and honey bees from a neighbouring apiary were screened for a broad set of putative pathogens including protists, fungi, spiroplasmas and viruses. Most sampled bees appeared to be infected with multiple parasites. Interestingly, viruses exclusively known from honey bees such as Apis mellifera Filamentous Virus and Varroa destructor Macula-like Virus were also discovered in solitary bees. A microsporidium found in Andrena vaga showed most resemblance to Nosema thomsoni. Our results suggest that bee hives represent a putative source of pathogens for other pollinators. Similarly, solitary bees may act as a reservoir of honey bee pathogens.

The health of honey bees is threatened by multiple factors, including viruses and parasites. We screened 557 honey bee (Apis mellifera) colonies from 155 beekeepers distributed all over Belgium to determine the prevalence of seven widespread viruses and two parasites (Varroa sp. and Nosema sp.). Deformed wing virus B (DWV-B), black queen cell virus (BQCV), and sacbrood virus (SBV) were highly prevalent and detected by real-time RT-PCR in more than 95% of the colonies. Acute bee paralysis virus (ABPV), chronic bee paralysis virus (CBPV) and deformed wing virus A (DWV-A) were prevalent to a lower extent (between 18 and 29%). Most viruses were only present at low or moderate viral loads. Nevertheless, about 50% of the colonies harbored at least one virus at high viral load (>107 genome copies/bee). Varroa mites and Nosema sp. were found in 81.5% and 59.7% of the honey bee colonies, respectively, and all Nosema were identified as Nosema ceranae by real time PCR. Interestingly, we found a significant correlation between the number of Varroa mites and DWV-B viral load. To determine the combined effect of these and other factors on honey bee health in Belgium, a follow up of colonies over multiple years is necessary.

Venoms of invertebrates contain an enormous diversity of proteins, peptides, and other classes of substances. Insect venoms are characterized by a large interspecific variation resulting in extended lists of venom compounds. The venom composition of several hymenopterans also shows different intraspecific variation. For instance, venom from different honeybee castes, more specifically queens and workers, shows quantitative and qualitative variation, while the environment, like seasonal changes, also proves to be an important factor. The present study aimed at an in-depth analysis of the intraspecific variation in the honeybee venom proteome. In summer workers, the recent list of venom proteins resulted from merging combinatorial peptide ligand library sample pretreatment and targeted tandem mass spectrometry realized with a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS/MS). Now, the same technique was used to determine the venom proteome of queens and winter bees, enabling us to compare it with that of summer bees. In total, 34 putative venom toxins were found, of which two were never described in honeybee venoms before. Venom from winter workers did not contain toxins that were not present in queens or summer workers, while winter worker venom lacked the allergen Api m 12, also known as vitellogenin. Venom from queen bees, on the other hand, was lacking six of the 34 venom toxins compared to worker bees, while it contained two new venom toxins, in particularly serine proteinase stubble and antithrombin-III. Although people are hardly stung by honeybees during winter or by queen bees, these newly identified toxins should be taken into account in the characterization of a putative allergic response against Apis mellifera stings.

The eusocial societies of honeybees, where the queen is the only fertile female among tens of thousands sterile worker bees, have intrigued scientists for centuries. The proximate factors, which cause the inhibition of worker bee ovaries, remain largely unknown; as are the factors which cause the activation of worker ovaries upon the loss of queen and brood in the colony. In an attempt to reveal key players in the regulatory network, we made a proteomic comparison of hemolymph profiles of workers with completely activated ovaries vs. rudimentary ovaries. An unexpected finding of this study is the correlation between age matched worker sterility and the enrichment of Picorna-like virus proteins. Fertile workers, on the other hand, show the upregulation of potential components of the immune system. It remains to be investigated whether viral infections contribute to worker sterility directly or are the result of a weaker immune system of sterile workers.

Metagenomics studies have accelerated the discovery of novel or divergent viruses of the honey bee. However, most of these studies predominantly focused on RNA viruses, and many suffer from the relatively low abundance of viral nucleic acids in the samples (i.e., compared to that of the host). Here, we explored the virome of the Ethiopian honey bee, Apis mellifera simensis, using an unbiased metagenomic approach in which the next-generation sequencing step was preceded by an enrichment protocol for viral particles. Our study revealed five well-known bee viruses and 25 atypical virus species, most of which have never been found in A. mellifera before. The viruses belong to Iflaviridae, Dicistroviridae, Secoviridae, Partitiviridae, Parvoviridae, Potyviridae, and taxonomically unclassified families. Fifteen of these atypical viruses were most likely plant-specific, and the remaining ten were presumed to be insect-specific. Apis mellifera filamentous virus (AmFV) was found in one sampling site out of 10. Two samples contained high read counts of a virus similar to Diatraea saccharales densovirus (DsDNV), which is a virus that causes high mortality in the sugarcane borer. AmFV and the DsDNV-like virus were the only DNA viruses found. Three viruses that primarily infect Drosophila spp. were also discovered: La Jolla virus (LJV), Kilifi virus (KiV), and Thika virus. Our study suggests that phoretic varroa mites are involved in the transmission of LJV and KiV and that both viruses replicate in mites and adult bees. We also found an overwhelming dominance of the deformed wing virus type B variant, which fits well with the apparently harmless infestation by Varroa destructor. It was suggested that Ethiopian bees have developed tolerance against virus infections as the result of natural selection.

As scientists continue to pursue various 'omics-based research, there is a need for high quality data for the most fundamental 'omics of all: genomics. The bacterium Paenibacillus larvae is the causative agent of the honey bee disease American foulbrood. If untreated, it can lead to the demise of an entire hive; the highly social nature of bees also leads to easy disease spread, between both individuals and colonies. Biologists have studied this organism since the early 1900s, and a century later, the molecular mechanism of infection remains elusive. Transcriptomics and proteomics, because of their ability to analyze multiple genes and proteins in a high-throughput manner, may be very helpful to its study. However, the power of these methodologies is severely limited without a complete genome; we undertake to address that deficiency here.

The diversity of eukaryote organisms and viruses associated with wild bees remains poorly characterized in contrast to the well-documented pathosphere of the western honey bee, Apis mellifera. Using a deliberate RNA shotgun metagenomic sequencing strategy in combination with a dedicated bioinformatics workflow, we identified the (micro-)organisms and viruses associated with two bumble bee hosts, Bombus terrestris and Bombus pascuorum, and two solitary bee hosts, Osmia cornuta and Andrena vaga. Ion Torrent semiconductor sequencing generated approximately 3.8 million high quality reads. The most significant eukaryote associations were two protozoan, Apicystis bombi and Crithidia bombi, and one nematode parasite Sphaerularia bombi in bumble bees. The trypanosome protozoan C. bombi was also found in the solitary bee O. cornuta. Next to the identification of three honey bee viruses Black queen cell virus, Sacbrood virus and Varroa destructor virus-1 and four plant viruses, we describe two novel RNA viruses Scaldis River bee virus (SRBV) and Ganda bee virus (GABV) based on their partial genomic sequences. The novel viruses belong to the class of negative-sense RNA viruses, SRBV is related to the order Mononegavirales whereas GABV is related to the family Bunyaviridae. The potential biological role of both viruses in bees is discussed in the context of recent advances in the field of arthropod viruses. Further, fragmentary sequence evidence for other undescribed viruses is presented, among which a nudivirus in O. cornuta and an unclassified virus related to Chronic bee paralysis virus in B. terrestris. Our findings extend the current knowledge of wild bee parasites in general and addsto the growing evidence of unexplored arthropod viruses in valuable insects.

Honey bees (Apis mellifera) are susceptible to a wide range of pathogens, including a broad set of viruses. Recently, next-generation sequencing has expanded the list of viruses with, for instance, two strains of Lake Sinai Virus. Soon after its discovery in the USA, LSV was also discovered in other countries and in other hosts. In the present study, we assemble four almost complete LSV genomes, and show that there is remarkable sequence heterogeneity based on the Orf1, RNA-dependent RNA polymerase and capsid protein sequences in comparison to the previously identified LSV 1 and 2 strains. Phylogenetic analyses of LSV sequences obtained from single honey bee specimens further revealed that up to three distinctive clades could be present in a single bee. Such superinfections have not previously been identified for other honey bee viruses. In a search for the putative routes of LSV transmission, we were able to demonstrate the presence of LSV in pollen pellets and in Varroa destructor mites. However, negative-strand analyses demonstrated that the virus only actively replicates in honey bees and mason bees (Osmia cornuta) and not in Varroa mites.

In this study, different context-dependent effects of imidacloprid exposure on the honey bee response were studied. Honey bees were exposed to different concentrations of imidacloprid during a time period of 40 days. Next to these variables, a laboratory-field comparison was conducted. The influence of the chronic exposure on gene expression levels was determined using an in-house developed microarray targeting different immunity-related and detoxification genes to determine stress-related gene expression changes. Increased levels of the detoxification genes encoding, CYP9Q3 and CYT P450, were detected in imidacloprid-exposed honey bees. The different context-dependent effects of imidacloprid exposure on honey bees were confirmed physiologically by decreased hypopharyngeal gland sizes. Honey bees exposed to imidacloprid in laboratory cages showed a general immunosuppression and no detoxification mechanisms were triggered significantly, while honey bees in-field showed a resilient response with an immune stimulation at later time points. However, the treated colonies had a brood and population decline tendency after the first brood cycle in the field. In conclusion, this study highlighted the different context-dependent effects of imidacloprid exposure on the honey bee response. These findings warn for possible pitfalls concerning the generalization of results based on specific experiments with short exposure times. The increased levels of CYT P450 and CYP9Q3 combined with an immune response reaction can be used as markers for bees which are exposed to pesticides in the field.

Varroa destructor mites and the viruses it vectors are two major factors leading to high losses of honey bees (Apis mellifera) colonies worldwide. However, honey bees in some African countries show resilience to varroa infestation and/or virus infections, although little is known about the mechanisms underlying this resilience. In this study, we investigated the expression profiles of some key molecular markers involved in olfactory sensing and RNA interference, as these processes may contribute to the bees' resilience to varroa infestation and virus infection, respectively. We found significantly higher gene expression of the odorant binding protein, OBP14, in the antennae of Ethiopian bees compared to Belgian bees. This result suggests the potential of OBP14 as a molecular marker of resilience to mite infestation. Scanning electron microscopy showed no significant differences in the antennal sensilla occurrence and distribution, suggesting that resilience arises from molecular processes rather than morphological adaptations. In addition, seven RNAi genes were upregulated in the Ethiopian honey bees and three of them-Dicer-Drosha, Argonaute 2, and TRBP2-were positively correlated with the viral load. We can conclude that the antiviral immune response was triggered when bees were experiencing severe viral infection and that this might contribute to the bees' resilience to viruses.

Bees are associated with a remarkable diversity of microorganisms, including unicellular parasites, bacteria, fungi, and viruses. The application of next-generation sequencing approaches enables the identification of this rich species composition as well as the discovery of previously unknown associations. Using high-throughput polyadenylated ribonucleic acid (RNA) sequencing, we investigated the metatranscriptome of eight wild bee species (Andrena cineraria, Andrena fulva, Andrena haemorrhoa, Bombus terrestris, Bombus cryptarum, Bombus pascuorum, Osmia bicornis, and Osmia cornuta) sampled from four different localities in Belgium. Across the RNA sequencing libraries, 88-99% of the taxonomically informative reads were of the host transcriptome. Four viruses with homology to insect pathogens were found including two RNA viruses (belonging to the families Iflaviridae and Tymoviridae that harbor already viruses of honey bees), a double stranded DNA virus (family Nudiviridae) and a single stranded DNA virus (family Parvoviridae). In addition, we found genomic sequences of 11 unclassified arthropod viruses (related to negeviruses, sobemoviruses, totiviruses, rhabdoviruses, and mononegaviruses), seven plant pathogenic viruses, and one fungal virus. Interestingly, nege-like viruses appear to be widespread, host-specific, and capable of attaining high copy numbers inside bees. Next to viruses, three novel parasite associations were discovered in wild bees, including Crithidia pragensis and a tubulinosematid and a neogregarine parasite. Yeasts of the genus Metschnikowia were identified in solitary bees. This study gives a glimpse of the microorganisms and viruses associated with social and solitary wild bees and demonstrates that their diversity exceeds by far the subset of species first discovered in honey bees.

Monitoring virus infections can be an important selection tool in honey bee breeding. A recent study pointed towards an association between the virus-free status of eggs and an increased virus resistance to deformed wing virus (DWV) at the colony level. In this study, eggs from both naturally surviving and traditionally managed colonies from across Europe were screened for the prevalence of different viruses. Screenings were performed using the phenotyping protocol of the 'suppressed in ovo virus infection' trait but with qPCR instead of end-point PCR and a primer set that covers all DWV genotypes. Of the 213 screened samples, 109 were infected with DWV, 54 were infected with black queen cell virus (BQCV), 3 were infected with the sacbrood virus, and 2 were infected with the acute bee paralyses virus. It was demonstrated that incidences of the vertical transmission of DWV were more frequent in naturally surviving than in traditionally managed colonies, although the virus loads in the eggs remained the same. When comparing virus infections with queen age, older queens showed significantly lower infection loads of DWV in both traditionally managed and naturally surviving colonies, as well as reduced DWV infection frequencies in traditionally managed colonies. We determined that the detection frequencies of DWV and BQCV in honey bee eggs were lower in samples obtained in the spring than in those collected in the summer, indicating that vertical transmission may be lower in spring. Together, these patterns in vertical transmission show that honey bee queens have the potential to reduce the degree of vertical transmission over time.

Trypanosomatids infecting honey bees have been poorly studied with molecular methods until recently. After the description of Crithidia mellificae (Langridge and McGhee, 1967) it took about forty years until molecular data for honey bee trypanosomatids became available and were used to identify and describe a new trypanosomatid species from honey bees, Lotmaria passim (Evans and Schwarz, 2014). However, an easy method to distinguish them without sequencing is not yet available. Research on the related bumble bee parasites Crithidia bombi and Crithidia expoeki revealed a fragment length polymorphism in the internal transcribed spacer 1 (ITS1), which enabled species discrimination. In search of fragment length polymorphisms for differential diagnostics in honey bee trypanosomatids, we studied honey bee trypanosomatid cell cultures of C. mellificae and L. passim. This research resulted in the identification of fragment length polymorphisms in ITS1 and ITS1-2 markers, which enabled us to develop a diagnostic method to differentiate both honey bee trypanosomatid species without the need for sequencing. However, the amplification success of the ITS1 marker depends probably on the trypanosomatid infection level. Further investigation confirmed that L. passim is the dominant species in Belgium, Japan and Switzerland. We found C. mellificae only rarely in Belgian honey bee samples, but not in honey bee samples from other countries. C. mellificae was also detected in mason bees (Osmia bicornis and Osmia cornuta) besides in honey bees. Further, the characterization and comparison of additional markers from L. passim strain SF (published as C. mellificae strain SF) and a Belgian honey bee sample revealed very low divergence in the 18S rRNA, ITS1-2, 28S rRNA and cytochrome b sequences. Nevertheless, a variable stretch was observed in the gp63 virulence factor.

Honey bees (Apis mellifera) produce an enormous economic value through their pollination activities and play a central role in the biodiversity of entire ecosystems. Recent efforts have revealed the substantial influence that the gut microbiota exert on bee development, food digestion, and homeostasis in general. In this study, deep sequencing was used to characterize prokaryotic viral communities associated with honey bees, which was a blind spot in research up until now. The vast majority of the prokaryotic viral populations are novel at the genus level, and most of the encoded proteins comprise unknown functions. Nevertheless, genomes of bacteriophages were predicted to infect nearly every major bee-gut bacterium, and functional annotation and auxiliary metabolic gene discovery imply the potential to influence microbial metabolism. Furthermore, undiscovered genes involved in the synthesis of secondary metabolic biosynthetic gene clusters reflect a wealth of previously untapped enzymatic resources hidden in the bee bacteriophage community.

Host-parasite co-evolution history is lacking when parasites switch to novel hosts. This was the case for Western honey bees (Apis mellifera) when the ectoparasitic mite, Varroa destructor, switched hosts from Eastern honey bees (Apis cerana). This mite has since become the most severe biological threat to A. mellifera worldwide. However, some A. mellifera populations are known to survive infestations, largely by suppressing mite population growth. One known mechanism is suppressed mite reproduction (SMR), but the underlying genetics are poorly understood. Here, we take advantage of haploid drones, originating from one queen from the Netherlands that developed Varroa-resistance, whole exome sequencing and elastic-net regression to identify genetic variants associated with SMR in resistant honeybees. An eight variants model predicted 88% of the phenotypes correctly and identified six risk and two protective variants. Reproducing and non-reproducing mites could not be distinguished using DNA microsatellites, which is in agreement with the hypothesis that it is not the parasite but the host that adapted itself. Our results suggest that the brood pheromone-dependent mite oogenesis is disrupted in resistant hosts. The identified genetic markers have a considerable potential to contribute to a sustainable global apiculture.

Bombus atratus bumblebees from Colombia that were caught in the wild and from breeding programs were screened for a broad set of bee pathogens. We discovered for the first time Lake Sinai Virus and confirmed the infection by other common viruses. The prevalence of Apicystis bombi, Crithidia bombi and Nosema ceranae was remarkably high. According to other studies the former two could have been co-introduced in South America with exotic bumble bees as Bombus terrestris or Bombus ruderatus. Given the fact that none of these species occur in Colombia, our data puts a new light on the spread of these pathogens over the South American continent.

The worldwide spread of diseases is considered a major threat to biodiversity and a possible driver of the decline of pollinator populations, particularly when novel species or strains of parasites emerge. Previous studies have suggested that populations of introduced European honeybee (Apis mellifera) and bumblebee species (Bombus terrestris and Bombus ruderatus) in Argentina share the neogregarine parasite Apicystis bombi with the native bumblebee (Bombus dahlbomii). In this study we investigated whether A. bombi is acting as an emergent parasite in the non-native populations. Specifically, we asked whether A. bombi, recently identified in Argentina, was introduced by European, non-native bees. Using ITS1 and ITS2 to assess the parasite's intraspecific genetic variation in bees from Argentina and Europe, we found a largely unstructured parasite population, with only 15% of the genetic variation being explained by geographic location. The most abundant haplotype in Argentina (found in all 9 specimens of non-native species) was identical to the most abundant haplotype in Europe (found in 6 out of 8 specimens). Similarly, there was no evidence of structuring by host species, with this factor explaining only 17% of the genetic variation. Interestingly, parasites in native Bombus ephippiatus from Mexico were genetically distant from the Argentine and European samples, suggesting that sufficient variability does exist in the ITS region to identify continent-level genetic structure in the parasite. Thus, the data suggest that A. bombi from Argentina and Europe share a common, relatively recent origin. Although our data did not provide information on the direction of transfer, the absence of genetic structure across space and host species suggests that A. bombi may be acting as an emergent infectious disease across bee taxa and continents.

A diverse supply of pollen is an important factor for honey bee health, but information about the pollen diversity available to colonies at the landscape scale is largely missing. In this COLOSS study, beekeeper citizen scientists sampled and analyzed the diversity of pollen collected by honey bee colonies. As a simple measure of diversity, beekeepers determined the number of colors found in pollen samples that were collected in a coordinated and standardized way. Altogether, 750 beekeepers from 28 different regions from 24 countries participated in the two-year study and collected and analyzed almost 18,000 pollen samples. Pollen samples contained approximately six different colors in total throughout the sampling period, of which four colors were abundant. We ran generalized linear mixed models to test for possible effects of diverse factors such as collection, i.e., whether a minimum amount of pollen was collected or not, and habitat type on the number of colors found in pollen samples. To identify habitat effects on pollen diversity, beekeepers' descriptions of the surrounding landscape and CORINE land cover classes were investigated in two different models, which both showed that both the total number and the rare number of colors in pollen samples were positively affected by 'urban' habitats or 'artificial surfaces', respectively. This citizen science study underlines the importance of the habitat for pollen diversity for bees and suggests higher diversity in urban areas.

Since the last decade, unusually high honey bee colony losses have been reported mainly in North-America and Europe. Here, we report on a comprehensive bee pathogen screening in Belgium covering 363 bee colonies that were screened for 18 known disease-causing pathogens and correlate their incidence in summer with subsequent winter mortality. Our analyses demonstrate that, in addition to Varroa destructor, the presence of the trypanosomatid parasite Crithidia mellificae and the microsporidian parasite Nosema ceranae in summer are also predictive markers of winter mortality, with a negative synergy being observed between the two in terms of their effects on colony mortality. Furthermore, we document the first occurrence of a parasitizing phorid fly in Europe, identify a new fourth strain of Lake Sinai Virus (LSV), and confirm the presence of other little reported pathogens such as Apicystis bombi, Aphid Lethal Paralysis Virus (ALPV), Spiroplasma apis, Spiroplasma melliferum and Varroa destructor Macula-like Virus (VdMLV). Finally, we provide evidence that ALPV and VdMLV replicate in honey bees and show that viruses of the LSV complex and Black Queen Cell Virus tend to non-randomly co-occur together. We also noticed a significant correlation between the number of pathogen species and colony losses. Overall, our results contribute significantly to our understanding of honey bee diseases and the likely causes of their current decline in Europe.

The beekeeping sector is suffering from the detrimental effects of climate change, both directly and indirectly. Despite numerous studies conducted on this subject, large-scale research incorporating stakeholders' and beekeepers' perspectives has remained elusive. This study aims to bridge this gap by assessing the extent to which stakeholders involved in the European beekeeping sector and European beekeepers perceive and experience the impacts of climate change on their operations, and whether they had to adapt their practices accordingly. To this end, a mixed-methods study including in-depth stakeholder interviews (n = 41) and a pan-European beekeeper survey (n = 844) was completed within the frame of the EU-funded H2020-project B-GOOD. The development of the beekeeper survey was informed by insights from literature and the stakeholder interviews. The results highlighted significant regional disparities in the perceived impacts of climate change, with beekeepers in Southern European regions expressing more negative outlooks, while Northern European beekeepers reported more favourable experiences. Furthermore, survey analysis revealed beekeepers who were classified as 'heavily impacted' by climate change. These beekeepers reported lower average honey yields, higher colony winter loss rates and a stronger perceived contribution of honey bees to pollination and biodiversity, underscoring climate change's detrimental impacts on the beekeeping sector. Multinomial logistic regression revealed determinants of the likelihood of beekeepers being classified as 'heavily impacted' by climate change. This analysis indicates that Southern European beekeepers experienced a 10-fold likelihood of being classified as heavily impacted by climate change compared to Northern European beekeepers. Other significant factors distinguishing 'winners' and 'losers' were self-reported level of professionalism as a beekeeper (ranging from pure hobbyist to fully professional, Odds Ratio (OR) = 1.31), number of years active in beekeeping (OR = 1.02), availability of floral resources throughout the bee season (OR = 0.78), beehives located in a forested environment (OR = 1.34), and the presence of local policy measures addressing climate change-related challenges (OR = 0.76).