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.

Multiple Locus Variable number of tandem repeat Analysis: A molecular genotyping tool for Paenibacillus larvae.

American Foulbrood, caused by Paenibacillus larvae, is the most severe bacterial disease of honey bees (Apis mellifera). To perform genotyping of P. larvae in an epidemiological context, there is a need of a fast and cheap method with a high resolution. Here, we propose Multiple Locus Variable number of tandem repeat Analysis (MLVA). MLVA has been used for typing a collection of 209 P. larvae strains from which 23 different MLVA types could be identified. Moreover, the developed methodology not only permits the identification of the four Enterobacterial Repetitive Intergenic Consensus (ERIC) genotypes, but allows also a discriminatory subdivision of the most dominant ERIC type I and ERIC type II genotypes. A biogeographical study has been conducted showing a significant correlation between MLVA genotype and the geographical region where it was isolated.

Factors influencing the prevalence and infestation levels of Varroa destructor in honeybee colonies in two highland agro-ecological zones of Uganda.

Varroa mites are ecto-parasites of honeybees and are a threat to the beekeeping industry. We identified the haplotype of Varroa mites and evaluated potential factors that influence their prevalence and infestation levels in the eastern and western highland agro-ecological zones of Uganda. This was done by collecting samples of adult worker bees between December 2014 and September 2015 in two sampling moments. Samples of bees were screened for Varroa using the ethanol wash method and the mites were identified by molecular techniques. All DNA sequences obtained from sampled mite populations in the two zones were 100 % identical to the Korean Haplotype (AF106899). Mean mite prevalence in the apiaries was 40 and 53 % for the western and eastern zones, respectively, during the first sampling. Over the second sampling, mean mite prevalence increased considerably in the western (59 %) but not in the eastern (51 %) zone. Factors that were associated with Varroa mite infestation levels include altitude, nature of apiary slope and apiary management practices during the first sampling. Our results further showed that Varroa mites were spreading from lower to higher elevations. Feral colonies were also infested with Varroa mites at infestation levels not significantly different from those in managed colonies. Colony productivity and strength were not correlated to mite infestation levels. We recommend a long-term Varroa mite monitoring strategy in areas of varying landscape and land use factors for a clear understanding of possible changes in mite infestation levels among African honeybees for informed decision making.

Differential diagnosis of the honey bee trypanosomatids Crithidia mellificae and Lotmaria passim.

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.

Vertical transmission of honey bee viruses in a Belgian queen breeding program.

BACKGROUND: The Member States of European Union are encouraged to improve the general conditions for the production and marketing of apicultural products. In Belgium, programmes on the restocking of honey bee hives have run for many years. Overall, the success ratio of this queen breeding programme has been only around 50%. To tackle this low efficacy, we organized sanitary controls of the breeding queens in 2012 and 2014. RESULTS: We found a high quantity of viruses, with more than 75% of the egg samples being infected with at least one virus. The most abundant viruses were Deformed Wing Virus and Sacbrood Virus (≥40%), although Lake Sinai Virus and Acute Bee Paralysis Virus were also occasionally detected (between 10-30%). In addition, Aphid Lethal Paralysis Virus strain Brookings, Black Queen Cell Virus, Chronic Bee Paralysis Virus and Varroa destructor Macula-like Virus occurred at very low prevalences (≤5%). Remarkably, we found Apis mellifera carnica bees to be less infected with Deformed Wing Virus than Buckfast bees (p < 0.01), and also found them to have a lower average total number of infecting viruses (p < 0.001). This is a significant finding, given that Deformed Wing Virus has earlier been shown to be a contributory factor to winter mortality and Colony Collapse Disorder. Moreover, negative-strand detection of Sacbrood Virus in eggs was demonstrated for the first time. CONCLUSIONS: High pathogen loads were observed in this sanitary control program. We documented for the first time vertical transmission of some viruses, as well as significant differences between two honey bee races in being affected by Deformed Wing Virus. Nevertheless, we could not demonstrate a correlation between the presence of viruses and queen breeding efficacies.