Known and novel viruses in Belgian honey bees: yearly differences, spatial clustering, and associations with overwintering loss.

In recent years, managed honey bee colonies have been suffering from an increasing number of biotic and abiotic stressors, resulting in numerous losses of colonies worldwide. A pan-European study, EPILOBEE, estimated the colony loss in Belgium to be 32.4% in 2012 and 14.8% in 2013. In the current study, absolute viral loads of four known honey bee viruses (DWV-A, DWV-B, AmFV, and BMLV) and three novel putative honey bee viruses (Apis orthomyxovirus 1, apthili virus, and apparli virus) were determined in 300 Flemish honey bee samples, and associations with winter survival were determined. This revealed that, in addition to the known influence of DWV-A and DWV-B on colony health, one of the newly described viruses (apthili virus) shows a strong yearly difference and is also associated with winter survival. Furthermore, all scrutinized viruses revealed significant spatial clustering patterns, implying that despite the limited surface area of Flanders, local virus transmission is paramount. The vast majority of samples were positive for at least one of the seven investigated viruses, and up to 20% of samples were positive for at least one of the three novel viruses. One of those three, Apis orthomyxovirus 1, was shown to be a genuine honey bee-infecting virus, able to infect all developmental stages of the honey bee, as well as the Varroa destructor mite. These results shed light on the most prevalent viruses in Belgium and their roles in the winter survival of honey bee colonies. IMPORTANCE: The western honey bee (Apis mellifera) is a highly effective pollinator of flowering plants, including many crops, which gives honey bees an outstanding importance both ecologically and economically. Alarmingly high annual loss rates of managed honey bee colonies are a growing concern for beekeepers and scientists and have prompted a significant research effort toward bee health. Several detrimental factors have been identified, such as varroa mite infestation and disease from various bacterial and viral agents, but annual differences are often not elucidated. In this study, we utilize the viral metagenomic survey of the EPILOBEE project, a European research program for bee health, to elaborate on the most abundant bee viruses of Flanders. We complement the existing metagenomic data with absolute viral loads and their spatial and temporal distributions. Furthermore, we identify Apis orthomyxovirus 1 as a potentially emerging pathogen, as we find evidence for its active replication honey bees.

Comprehensive Approach to Phenotype Varroa destructor Reproduction in Honey Bee Drone Brood and Its Correlation with Decreased Mite Reproduction (DMR).

The mechanisms of action behind decreased mite reproduction (DMR) are still unknown, but current hypotheses state that DMR is the result of brood-intrinsic and/or external disturbances in the V. destructor-honey bee pupa signal interactions. For accurate and precise DMR phenotyping, sufficient single infested honey bee brood cells are required (e.g., 35), which requires extensive labor and time and may exclude many samples not reaching the threshold. We defined a new comprehensive trait called the 'mean V. destructor reproduction rate' (mVR), which describes the mean number of offspring mites per infested cell in the sample while compensating for the reduced number of offspring with increasing multiple infested cells. We found a significant correlation between mVR and DMR, allowing for an estimation of DMR based on the mVR only. When the mVR was calculated with 10 infested cells, we found an average variation in mVR of 16.8%. For the same variation in DMR determination, 40 single infested cells are required. This broader look at V. destructor resistance phenotyping can improve the applicability and effectiveness of traits related to V. destructor reproduction in honey bee breeding programs.

Allele Frequencies of Genetic Variants Associated with Varroa Drone Brood Resistance (DBR) in Apis mellifera Subspecies across the European Continent.

Implementation of marker-assisted selection (MAS) in modern beekeeping would improve sustainability, especially in breeding programs aiming for resilience against the parasitic mite Varroa destructor. Selecting honey bee colonies for natural resistance traits, such as brood-intrinsic suppression of varroa mite reproduction, reduces the use of chemical acaricides while respecting local adaptation. In 2019, eight genomic variants associated with varroa non-reproduction in drone brood were discovered in a single colony from the Amsterdam Water Dune population in the Netherlands. Recently, a new study tested the applicability of these eight genetic variants for the same phenotype on a population-wide scale in Flanders, Belgium. As the properties of some variants varied between the two studies, one hypothesized that the difference in genetic ancestry of the sampled colonies may underly these contribution shifts. In order to frame this, we determined the allele frequencies of the eight genetic variants in more than 360 Apis mellifera colonies across the European continent and found that variant type allele frequencies of these variants are primarily related to the A. mellifera subspecies or phylogenetic honey bee lineage. Our results confirm that population-specific genetic markers should always be evaluated in a new population prior to using them in MAS programs.

A field realistic model to assess the effects of pesticides residues and adulterants on honey bee gene expression.

While studies on the sublethal effects of chemical residues in beeswax on adult honey bees are increasing, the study protocols assessing the impacts on honey bee brood in realistic conditions still need to be investigated. Moreover, little is known about the residue's effect on gene expression in honey bee brood. This study reports the effects of chlorpyriphos-ethyl, acrinathrin and stearin worker pupae exposure through contaminated or adulterated beeswax on the gene expression of some key health indicators, using a novel in vivo realistic model. Larvae were reared in acrinathrin (12.5, 25, 10 and 100 ppb) and chlorpyriphos-ethyl (5, 10, 500 and 5000 ppb) contaminated or stearin adulterated beeswax (3, 4, 5, 6 and 9%) in newly formed colonies to reduce the influence of external factors. On day 11, mortality rates were assessed. Honey bee pupae were extracted from the comb after 19 days of rearing and were analysed for the gene expression profile of four genes involved in the immune response to pathogens and environmental stress factors (Imd, dorsal, domeless and defensin), and two genes involved in detoxifications mechanisms (CYP6AS14 and CYP9Q3). We found no linear relation between the increase in the pesticide concentrations and the brood mortality rates, unlike stearin where an increase in stearin percentage led to an exponential increase in brood mortality. The immune system of pupae raised in acrinathrin contaminated wax was triggered and the expression of CYP6AS14 was significantly upregulated (exposure to 12.5 and 25 ppb). Almost all expression levels of the tested immune and detoxification genes were down-regulated when pupae were exposed to chlorpyrifos-contaminated wax. The exposure to stearin triggered the immune system and detoxification system of the pupae. The identification of substance-specific response factors might ultimately serve to identify molecules that are safer for bees and the ecosystem's health.

Eggs sampling as an effective tool for identifying the incidence of viruses in honey bees involved in artificial queen rearing.

The Carniolan honey bee (Apis mellifera carnica) plays an essential role in crop pollination, environment diversity, and the production of honey bee products. However, the health of individual honey bees and their colonies is under pressure due to multiple stressors, including viruses as a significant threat to bees. Monitoring various virus infections could be a crucial selection tool during queen rearing. In the present study, samples from all developmental stages (eggs, larvae, pupae, and queens) were screened for the incidence of seven viruses during queen rearing in Slovenia. The screening of a total of 108 samples from five queen breeders was performed by the RT-qPCR assays. The results showed that the highest incidence was observed for black queen cell virus (BQCV), Lake Sinai virus 3 (LSV3), deformed wing virus B (DWV-B), and sacbrood virus (SBV). The highest viral load was detected in queens (6.07 log10 copies/queen) and larvae (5.50 log10 copies/larva) for BQCV, followed by SBV in larvae (5.47 log10 copies/larva). When comparing all the honey bee developmental stages, the eggs exhibited general screening for virus incidence and load in queen mother colonies. The results suggest that analyzing eggs is a good indicator of resilience to virus infection during queen development.

Population-wide modelling reveals prospects of marker-assisted selection for parasitic mite resistance in honey bees.

In 2019, a joint eight-variant model was published in which eight single nucleotide polymorphisms (SNPs) in seven Apis mellifera genes were associated with Varroa destructor drone brood resistance (DBR, i.e. mite non-reproduction in drone brood). As this model was derived from only one Darwinian Black Bee Box colony, it could not directly be applied on a population-overarching scale in the northern part of Belgium (Flanders), where beekeepers prefer the carnica subspecies. To determine whether these eight SNPs remained associated with the DBR trait on a Flemish colony-broad scope, we performed population-wide modelling through sampling of various A. mellifera carnica colonies, DBR scoring of Varroa-infested drone brood and variant genotyping. Novel eight-variant modelling was performed and the classification performance of the eight SNPs was evaluated. Besides, we built a reduced three-variant model retaining only three genetic variants and found that this model classified 76% of the phenotyped drones correctly. To examine the spread of beneficial alleles and predict the DBR probability distribution in Flanders, we determined the allelic frequencies of the three variants in 292 A. mellifera carnica queens. As such, this research reveals prospects of marker-assisted selection for Varroa drone brood resistance in honeybees.

Bridging the Gap between Field Experiments and Machine Learning: The EC H2020 B-GOOD Project as a Case Study towards Automated Predictive Health Monitoring of Honey Bee Colonies.

Honey bee colonies have great societal and economic importance. The main challenge that beekeepers face is keeping bee colonies healthy under ever-changing environmental conditions. In the past two decades, beekeepers that manage colonies of Western honey bees (Apis mellifera) have become increasingly concerned by the presence of parasites and pathogens affecting the bees, the reduction in pollen and nectar availability, and the colonies' exposure to pesticides, among others. Hence, beekeepers need to know the health condition of their colonies and how to keep them alive and thriving, which creates a need for a new holistic data collection method to harmonize the flow of information from various sources that can be linked at the colony level for different health determinants, such as bee colony, environmental, socioeconomic, and genetic statuses. For this purpose, we have developed and implemented the B-GOOD (Giving Beekeeping Guidance by computational-assisted Decision Making) project as a case study to categorize the colony's health condition and find a Health Status Index (HSI). Using a 3-tier setup guided by work plans and standardized protocols, we have collected data from inside the colonies (amount of brood, disease load, honey harvest, etc.) and from their environment (floral resource availability). Most of the project's data was automatically collected by the BEEP Base Sensor System. This continuous stream of data served as the basis to determine and validate an algorithm to calculate the HSI using machine learning. In this article, we share our insights on this holistic methodology and also highlight the importance of using a standardized data language to increase the compatibility between different current and future studies. We argue that the combined management of big data will be an essential building block in the development of targeted guidance for beekeepers and for the future of sustainable beekeeping.

Shift in virus composition in honeybees (Apis mellifera) following worldwide invasion by the parasitic mite and virus vector Varroa destructor.

Invasive vectors can induce dramatic changes in disease epidemiology. While viral emergence following geographical range expansion of a vector is well known, the influence a vector can have at the level of the host's pathobiome is less well understood. Taking advantage of the formerly heterogeneous spatial distribution of the ectoparasitic mite Varroa destructor that acts as potent virus vector among honeybees Apis mellifera, we investigated the impact of its recent global spread on the viral community of honeybees in a retrospective study of historical samples. We hypothesized that the vector has had an effect on the epidemiology of several bee viruses, potentially altering their transmissibility and/or virulence, and consequently their prevalence, abundance, or both. To test this, we quantified the prevalence and loads of 14 viruses from honeybee samples collected in mite-free and mite-infested populations in four independent geographical regions. The presence of the mite dramatically increased the prevalence and load of deformed wing virus, a cause of unsustainably high colony losses. In addition, several other viruses became more prevalent or were found at higher load in mite-infested areas, including viruses not known to be actively varroa-transmitted, but which may increase opportunistically in varroa-parasitized bees.

Expression of Molecular Markers of Resilience against Varroa destructor and Bee Viruses in Ethiopian Honey Bees (Apis mellifera simensis) Focussing on Olfactory Sensing and the RNA Interference Machinery.

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.

Virus Prevalence in Egg Samples Collected from Naturally Selected and Traditionally Managed Honey Bee Colonies across Europe.

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.

Breeding for Virus Resistance and Its Effects on Deformed Wing Virus Infection Patterns in Honey Bee Queens.

Viruses, and in particular the deformed wing virus (DWV), are considered as one of the main antagonists of honey bee health. The 'suppressed in ovo virus infection' trait (SOV) described for the first time that control of a virus infection can be achieved from genetically inherited traits and that the virus state of the eggs is indicative for this. This research aims to explore the effect of the SOV trait on DWV infections in queens descending from both SOV-positive (QDS+) and SOV-negative (QDS-) queens. Twenty QDS+ and QDS- were reared from each time four queens in the same starter-finisher colony. From each queen the head, thorax, ovaries, spermatheca, guts and eviscerated abdomen were dissected and screened for the presence of the DWV-A and DWV-B genotype using qRT-PCR. Queens descending from SOV-positive queens showed significant lower infection loads for DWV-A and DWV-B as well as a lower number of infected tissues for DWV-A. Surprisingly, differences were less expressed in the reproductive tissues, the ovaries and spermatheca. These results confirm that selection on the SOV trait is associated with increased virus resistance across viral genotypes and that this selection drives DWV towards an increased tissue specificity for the reproductive tissues. Further research is needed to explore the mechanisms underlying the interaction between the antiviral response and DWV.

qPCR assays with dual-labeled probes for genotyping honey bee variants associated with varroa resistance.

BACKGROUND: The varroa mite is one of the main causes of honey bee mortality. An important mechanism by which honey bees increase their resistance against this mite is the expression of suppressed mite reproduction. This trait describes the physiological inability of mites to produce viable offspring and was found associated with eight genomic variants in previous research. RESULTS: This paper presents the development and validation of high-throughput qPCR assays with dual-labeled probes for discriminating these eight single-nucleotide variants. Amplicon sequences used for assay validation revealed additional variants in the primer/probe binding sites in four out of the eight assays. As for two of these the additional variants interfered with the genotyping outcome supplementary primers and/or probes were developed. Inclusion of these primers and probes in the assay mixes allowed for the correct genotyping of all eight variants of interest within our bee population. CONCLUSION: These outcomes underline the importance of checking for interfering variants in designing qPCR assays. Ultimately, the availability of this assay allows genotyping for the suppressed mite reproduction trait and paves the way for marker assisted selection in breeding programs.

Metagenomic Approach with the NetoVIR Enrichment Protocol Reveals Virus Diversity within Ethiopian Honey Bees (Apis mellifera simensis).

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.

Heritability estimates of the novel trait 'suppressed in ovo virus infection' in honey bees (Apis mellifera).

Honey bees are under pressure due to abnormal high colony death rates, especially during the winter. The infestation by the Varroa destructor mite and the viruses that this ectoparasite transmits are generally considered as the bees' most important biological threats. Almost all efforts to remedy this dual infection have so far focused on the control of the Varroa mite alone and not on the viruses it transmits. In the present study, the sanitary control of breeding queens was conducted on eggs taken from drone brood for 4 consecutive years (2015-2018). The screening was performed on the sideline of an ongoing breeding program, which allowed us to estimate the heritabilities of the virus status of the eggs. We used the term 'suppressed in ovo virus infection' (SOV) for this novel trait and found moderate heritabilities for the presence of several viruses simultaneously and for the presence of single viral species. Colonies that expressed the SOV trait seemed to be more resilient to virus infections as a whole with fewer and less severe Deformed wing virus infections in most developmental stages, especially in the male caste. The implementation of this novel trait into breeding programs is recommended.

Honey bee predisposition of resistance to ubiquitous mite infestations.

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.

Characterization of the honeybee venom proteins C1q-like protein and PVF1 and their allergenic potential.

Honeybee (Apis mellifera) venom (HBV) represents an ideal model to study the role of particular venom components in allergic reactions in sensitized individuals as well as in the eusociality of Hymenoptera species. The aim of this study was to further characterize the HBV components C1q-like protein (C1q) and PDGF/VEGF-like factor 1 (PVF1). C1q and PVF1 were produced as recombinant proteins in insect cells. Their allergenic properties were examined by determining the level of specific IgE antibodies in the sera of HBV-allergic patients (n = 26) as well as by their capacity to activate patients' basophils (n = 11). Moreover, the transcript heterogeneity of PVF1 was analyzed. It could be demonstrated that at least three PVF1 variants are present in the venom gland, which all result from alternative splicing of one transcript. Additionally, recombinant C1q and PVF1 from Spodoptera frugiperda insect cells exhibited specific IgE reactivity with approximately 38.5% of sera of HBV-allergic patients. Interestingly, both proteins were unable to activate basophils of the patients, questioning their role in the context of clinically relevant sensitization. Recombinant C1q and PVF1 can build the basis for a deeper understanding of the molecular mechanisms of Hymenoptera venoms. Moreover, the conflicting results between IgE sensitization and lack of basophil activation, might in the future contribute to the identification of factors that determine the allergenic potential of proteins.

Nosema neumanni n. sp. (Microsporidia, Nosematidae), a new microsporidian parasite of honeybees, Apis mellifera in Uganda.

The microsporidium Nosema neumanni n. sp., a new parasite of the honeybee Apis mellifera is described based on its ultra-structural and molecular characteristics. Structures resembling microsporidian spores were found by microscopic examination of honeybees from Uganda. Molecular confirmation failed when PCR primers specific for Nosema apis and Nosema ceranae were used, but was successful with primers covering the whole family of Nosematidae. We performed transmission electron microscopy and found typical microsporidian spores which were smaller (length: 2.36±0.14µm and width: 1.78±0.06µm; n=6) and had fewer polar filament coils (10-12) when compared to those of known species infecting honeybees. The entire 16S SSU rRNA region was amplified, cloned and sequenced and was found to be unique with the highest resemblance (97% identity) to N. apis. The incidence of N. neumanni n. sp. in Ugandan honeybees was found to be much higher than of the two other Nosema species.

Covert deformed wing virus infections have long-term deleterious effects on honeybee foraging and survival.

Several studies have suggested that covert stressors can contribute to bee colony declines. Here we provide a novel case study and show using radiofrequency identification tracking technology that covert deformed wing virus (DWV) infections in adult honeybee workers seriously impact long-term foraging and survival under natural foraging conditions. In particular, our experiments show that adult workers injected with low doses of DWV experienced increased mortality rates, that DWV caused workers to start foraging at a premature age, and that the virus reduced the workers' total activity span as foragers. Altogether, these results demonstrate that covert DWV infections have strongly deleterious effects on honeybee foraging and survival. These results are consistent with previous studies that suggested DWV to be an important contributor to the ongoing bee declines in Europe and the USA. Overall, our study underlines the strong impact that covert pathogen infections can have on individual and group-level performance in bees.

Stress indicator gene expression profiles, colony dynamics and tissue development of honey bees exposed to sub-lethal doses of imidacloprid in laboratory and field experiments.

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.

Unbiased RNA Shotgun Metagenomics in Social and Solitary Wild Bees Detects Associations with Eukaryote Parasites and New Viruses.

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.