Complex population structure and haplotype patterns in the Western European honey bee from sequencing a large panel of haploid drones.

Honey bee subspecies originate from specific geographical areas in Africa, Europe and the Middle East, and beekeepers interested in specific phenotypes have imported genetic material to regions outside of the bees' original range for use either in pure lines or controlled crosses. Moreover, imported drones are present in the environment and mate naturally with queens from the local subspecies. The resulting admixture complicates population genetics analyses, and population stratification can be a major problem for association studies. To better understand Western European honey bee populations, we produced a whole genome sequence and single nucleotide polymorphism (SNP) genotype data set from 870 haploid drones and demonstrate its utility for the identification of nine genetic backgrounds and various degrees of admixture in a subset of 629 samples. Five backgrounds identified correspond to subspecies, two to isolated populations on islands and two to managed populations. We also highlight several large haplotype blocks, some of which coincide with the position of centromeres. The largest is 3.6 Mb long and represents 21% of chromosome 11, with two major haplotypes corresponding to the two dominant genetic backgrounds identified. This large naturally phased data set is available as a single vcf file that can now serve as a reference for subsequent populations genomics studies in the honey bee, such as (i) selecting individuals of verified homogeneous genetic backgrounds as references, (ii) imputing genotypes from a lower-density data set generated by an SNP-chip or by low-pass sequencing, or (iii) selecting SNPs compatible with the requirements of genotyping chips.

Population Structure and Diversity in European Honey Bees (Apismellifera L.)-An Empirical Comparison of Pool and Individual Whole-Genome Sequencing.

BACKGROUND: Whole-genome sequencing has become routine for population genetic studies. Sequencing of individuals provides maximal data but is rather expensive and fewer samples can be studied. In contrast, sequencing a pool of samples (pool-seq) can provide sufficient data, while presenting less of an economic challenge. Few studies have compared the two approaches to infer population genetic structure and diversity in real datasets. Here, we apply individual sequencing (ind-seq) and pool-seq to the study of Western honey bees (Apis mellifera). METHODS: We collected honey bee workers that belonged to 14 populations, including 13 subspecies, totaling 1347 colonies, who were individually (139 individuals) and pool-sequenced (14 pools). We compared allele frequencies, genetic diversity estimates, and population structure as inferred by the two approaches. RESULTS: Pool-seq and ind-seq revealed near identical population structure and genetic diversities, albeit at different costs. While pool-seq provides genome-wide polymorphism data at considerably lower costs, ind-seq can provide additional information, including the identification of population substructures, hybridization, or individual outliers. CONCLUSIONS: If costs are not the limiting factor, we recommend using ind-seq, as population genetic structure can be inferred similarly well, with the advantage gained from individual genetic information. Not least, it also significantly reduces the effort required for the collection of numerous samples and their further processing in the laboratory.

Resistance and Vulnerability of Honeybee (Apis mellifera) Gut Bacteria to Commonly Used Pesticides.

Agricultural and apicultural practices expose honeybees to a range of pesticides that have the potential to negatively affect their physiology, neurobiology, and behavior. Accumulating evidence suggests that these effects extend to the honeybee gut microbiome, which serves important functions for honeybee health. Here we test the potential effects of the pesticides thiacloprid, acetamiprid, and oxalic acid on the gut microbiota of honeybees, first in direct in vitro inhibition assays and secondly in an in vivo caged bee experiment to test if exposure leads to gut microbiota community changes. We found that thiacloprid did not inhibit the honeybee core gut bacteria in vitro, nor did it affect overall community composition or richness in vivo. Acetamiprid did also not inhibit bacterial growth in vitro, but it did affect community structure within bees. The eight bacterial genera tested showed variable levels of susceptibility to oxalic acid in vitro. In vivo, treatment with this pesticide reduced amplicon sequence variant (ASV) richness and affected gut microbiome composition, with most marked impact on the common crop bacteria Lactobacillus kunkeei and the genus Bombella. We conducted network analyses which captured known associations between bacterial members and illustrated the sensitivity of the microbiome to environmental stressors. Our findings point to risks of honeybee exposure to oxalic acid, which has been deemed safe for use in treatment against Varroa mites in honeybee colonies, and we advocate for more extensive assessment of the long-term effects that it may have on honeybee health.

Pollination and Plant Reproductive Success of Two Ploidy Levels in Red Clover (Trifolium pratense L.).

Plant reproduction in red clover requires cross-fertilization via insect pollination. However, the influences of visitation rate and timing on maximizing ovule utilization are yet to be determined. We aimed to study the influences of visitation rate, flowering stage, and self-incompatibility on reproductive success. We applied hand and honey bee pollination in the study of eight red clover cultivars with two ploidy levels released between 1964 and 2001. In hand pollination, increasing the visitation rates (from 10 to 80 pollinated florets per flower head) increased the seed number per flower head but reduced the seed number per pollinated floret. Different flowering stages (early, middle, and full flowering) did not influence the seed number per pollinated floret significantly. There was a marked difference in reproductive success depending on the ploidy level, with 0.52 seeds per pollinated floret in diploid and 0.16 in tetraploid cultivars. During the cultivar release history, seed number per pollinated floret seemed to decrease in diploid cultivars, whereas it increased in tetraploids. In honey bee pollination, diploid cultivars had more two-seeded florets than tetraploids. Different visitation rates and the stochastic nature of pollen transfer resulted in difficulties when the plant reproductive success between hand and bee pollination was compared. A maximum of 0.27 seeds per pollinated floret were produced in hand pollination compared to the 0.34 in honey bee pollination. In spite of this, hand pollination provided a valuable method for studying the pollination biology and reproduction of red clover. Future studies may employ hand pollination to unravel further aspects of the low reproductive success with the future perspective of improving seed number per pollinated floret in tetraploid red clover.

Evaluating competition for forage plants between honey bees and wild bees in Denmark.

A recurrent concern in nature conservation is the potential competition for forage plants between wild bees and managed honey bees. Specifically, that the highly sophisticated system of recruitment and large perennial colonies of honey bees quickly exhaust forage resources leading to the local extirpation of wild bees. However, different species of bees show different preferences for forage plants. We here summarize known forage plants for honey bees and wild bee species at national scale in Denmark. Our focus is on floral resources shared by honey bees and wild bees, with an emphasis on both threatened wild bee species and foraging specialist species. Across all 292 known bee species from Denmark, a total of 410 plant genera were recorded as forage plants. These included 294 plant genera visited by honey bees and 292 plant genera visited by different species of wild bees. Honey bees and wild bees share 176 plant genera in Denmark. Comparing the pairwise niche overlap for individual bee species, no significant relationship was found between their overlap and forage specialization or conservation status. Network analysis of the bee-plant interactions placed honey bees aside from most other bee species, specifically the module containing the honey bee had fewer links to any other modules, while the remaining modules were more highly inter-connected. Despite the lack of predictive relationship from the pairwise niche overlap, data for individual species could be summarized. Consequently, we have identified a set of operational parameters that, based on a high foraging overlap (>70%) and unfavorable conservation status (Vulnerable+Endangered+Critically Endangered), can guide both conservation actions and land management decisions in proximity to known or suspected populations of these species.

Metabolomics unveils the influence of dietary phytochemicals on residual pesticide concentrations in honey bees.

The losses of honey bee colonies and declines of other insect pollinators have been associated with negative effects of pesticides. Honey bees as well as other pollinators are nectar and pollen foragers and thus are exposed to an extensive range of phytochemicals. Understanding the synergistic, additive, and antagonistic effects of plant secondary metabolites and pesticides in honey bees may help to protect honey bee colonies against agrochemicals. In this study, we used untargeted metabolomics to investigate the impact of dietary phytochemical composition on the residual concentration of three pesticides: imidacloprid, tau-fluvalinate and tebuconazole in exposed honey bees. Honey bees were given different diets based on pollen or nectar from four plants: Reseda odorata, Borago officinalis, Phacelia tanacetifolia, and Trifolium repens for two days. Thereafter, they were orally exposed to 10 ng/bee imidacloprid or contact-exposed to 0.9 µg/bee tau-fluvalinate or 5 µg/bee tebuconazole. After 1 h of oral exposure or 24 h of contact exposure, the honey bees were anaesthetised with CO2, sacrificed by freezing, extracted with a validated QuEChERS method, and residual pesticide concentrations were determined by LC-QTRAP-MS/MS. The phytochemical composition in the given diets were profiled with an UHPLC-Q Exactive-MS/MS. The results revealed that the dietary phytochemical composition has a noteworthy influence on the concentration of residual pesticides in honey bees. The correlation coefficient analysis demonstrated that flavonoids have a reducing effect on the residual concentration of imidacloprid and tau-fluvalinate in honey bees. The results also highlighted that exposure to imidacloprid impaired the metabolism of sugars in honey bees. Exploiting flavonoid-rich plants may protect honey bees against pesticides and hold promise as forage plants in future beekeeping.

Authoritative subspecies diagnosis tool for European honey bees based on ancestry informative SNPs.

BACKGROUND: With numerous endemic subspecies representing four of its five evolutionary lineages, Europe holds a large fraction of Apis mellifera genetic diversity. This diversity and the natural distribution range have been altered by anthropogenic factors. The conservation of this natural heritage relies on the availability of accurate tools for subspecies diagnosis. Based on pool-sequence data from 2145 worker bees representing 22 populations sampled across Europe, we employed two highly discriminative approaches (PCA and FST) to select the most informative SNPs for ancestry inference. RESULTS: Using a supervised machine learning (ML) approach and a set of 3896 genotyped individuals, we could show that the 4094 selected single nucleotide polymorphisms (SNPs) provide an accurate prediction of ancestry inference in European honey bees. The best ML model was Linear Support Vector Classifier (Linear SVC) which correctly assigned most individuals to one of the 14 subspecies or different genetic origins with a mean accuracy of 96.2% ± 0.8 SD. A total of 3.8% of test individuals were misclassified, most probably due to limited differentiation between the subspecies caused by close geographical proximity, or human interference of genetic integrity of reference subspecies, or a combination thereof. CONCLUSIONS: The diagnostic tool presented here will contribute to a sustainable conservation and support breeding activities in order to preserve the genetic heritage of European honey bees.

LC-MS/MS Quantification Reveals Ample Gut Uptake and Metabolization of Dietary Phytochemicals in Honey Bees (Apis mellifera).

The honey bee pollen/nectar diet is rich in bioactive phytochemicals and recent studies have demonstrated the potential of phytochemicals to influence honey bee disease resistance. To unravel the role of dietary phytochemicals in honey bee health it is essential to understand phytochemical uptake, bioavailability, and metabolism but presently limited knowledge exists. With this study we aim to build a knowledge foundation. For 5 days, we continuously fed honey bees on eight individual phytochemicals and measured the concentrations in whole and dissected bees by HPLC-MS/MS. Ample phytochemical metabolization was observed, and only 6-30% of the consumed quantities were recovered. Clear differences in metabolization rates were evident, with atropine, aucubin, and triptolide displaying significantly slower metabolism. Phytochemical gut uptake was also demonstrated, and oral bioavailability was 4-31%, with the highest percentages observed for amygdalin, triptolide, and aucubin. We conclude that differences in the chemical properties and structure impact phytochemical uptake and metabolism.

Dietary quercetin impacts the concentration of pesticides in honey bees.

Honey bees are important pollinators and are subject to numerous stressors, such as changing floral resources, parasites, and agrochemical exposure. Pesticide exposure has been linked to the decline in the global honey bee population. We have limited knowledge of the metabolic pathways and synergistic effects of xenobiotics in bees. Quercetin is one of the most abundant phytochemicals in plants and is therefore abundant in the honey bee diet. Quercetin can upregulate the detoxification system in honey bees; however, it is still unknown to what extent quercetin ingestion can reduce the content of absorbed pesticides. In this study, we investigated the effect of dietary quercetin on the contents of three pesticides in honey bees: imidacloprid (insecticide), tebuconazole (fungicide), and tau-fluvalinate (insecticide and acaricide). Bees were divided into two main groups and fed either quercetin-sucrose paste or only sucrose for 72 h. Thereafter, they were orally exposed to ∼10 ng/bee imidacloprid or contact-exposed to ∼0.9 µg/bee tau-fluvalinate or ∼5.2 µg/bee tebuconazole. After 1 h of oral exposure or 24 h of contact exposure, the bees were anaesthetised with CO2, sacrificed by freezing, and extracted with a validated QuEChERS method. Subsequently, the concentrations of the three pesticides and quercetin in the bees were determined with a triple quadrupole tandem mass spectrometer coupled to an HPLC system. No significant effect on the concentration of tebuconazole or tau-fluvalinate was observed in bees fed quercetin. Intake of quercetin led to a reduction in the concentration of imidacloprid in honey bees. Quercetin-rich plants may be exploited in future beekeeping.

Evaluation of Suppressed Mite Reproduction (SMR) Reveals Potential for Varroa Resistance in European Honey Bees (Apis mellifera L.).

In the fight against the Varroa destructor mite, selective breeding of honey bee (Apis mellifera L.) populations that are resistant to the parasitic mite stands as a sustainable solution. Selection initiatives indicate that using the suppressed mite reproduction (SMR) trait as a selection criterion is a suitable tool to breed such resistant bee populations. We conducted a large European experiment to evaluate the SMR trait in different populations of honey bees spread over 13 different countries, and representing different honey bee genotypes with their local mite parasites. The first goal was to standardize and validate the SMR evaluation method, and then to compare the SMR trait between the different populations. Simulation results indicate that it is necessary to examine at least 35 single-infested cells to reliably estimate the SMR score of any given colony. Several colonies from our dataset display high SMR scores indicating that this trait is present within the European honey bee populations. The trait is highly variable between colonies and some countries, but no major differences could be identified between countries for a given genotype, or between genotypes in different countries. This study shows the potential to increase selective breeding efforts of V. destructor resistant populations.

High sample throughput genotyping for estimating C-lineage introgression in the dark honeybee: an accurate and cost-effective SNP-based tool.

The natural distribution of the honeybee (Apis mellifera L.) has been changed by humans in recent decades to such an extent that the formerly widest-spread European subspecies, Apis mellifera mellifera, is threatened by extinction through introgression from highly divergent commercial strains in large tracts of its range. Conservation efforts for A. m. mellifera are underway in multiple European countries requiring reliable and cost-efficient molecular tools to identify purebred colonies. Here, we developed four ancestry-informative SNP assays for high sample throughput genotyping using the iPLEX Mass Array system. Our customized assays were tested on DNA from individual and pooled, haploid and diploid honeybee samples extracted from different tissues using a diverse range of protocols. The assays had a high genotyping success rate and yielded accurate genotypes. Performance assessed against whole-genome data showed that individual assays behaved well, although the most accurate introgression estimates were obtained for the four assays combined (117 SNPs). The best compromise between accuracy and genotyping costs was achieved when combining two assays (62 SNPs). We provide a ready-to-use cost-effective tool for accurate molecular identification and estimation of introgression levels to more effectively monitor and manage A. m. mellifera conservatories.

Quantitative patterns of vertical transmission of deformed wing virus in honey bees.

Deformed wing virus (DWV) is an important pathogen in a broad range of insects, including honey bees. Concordant with the spread of Varroa, DWV is present in the majority of honey bee colonies and can result in either low-level infections with asymptomatic bees that nonetheless exhibit increased colony loss under stress, or high-level infections with acute effects on bee health and viability. DWV can be transmitted vertically or horizontally and evidence suggests that horizontal transmission via Varroa is associated with acute symptomatic infections. Vertical transmission also occurs and is presumably important for the maintenance of DWV in honey bee populations. To further our understanding the vertical transmission of DWV through queens, we performed three experiments: we studied the quantitative effectiveness of vertical transmission, surveyed the prevalence of successful egg infection under commercial conditions, and distinguished among three possible mechanisms of transmission. We find that queen-infection level predicts the DWV titers in their eggs, although the transmission is not very efficient. Our quantitative assessment of DWV demonstrates that eggs in 1/3 of the colonies are infected with DWV and highly infected eggs are rare in newly-installed spring colonies. Additionally, our results indicate that DWV transmission occurs predominantly by virus adhering to the surface of eggs (transovum) rather than intracellularly. Our combined results suggest that the queens' DWV vectoring capacity in practice is not as high as its theoretical potential. Thus, DWV transmission by honey bee queens is part of the DWV epidemic with relevant practical implications, which should be further studied.

Infection dynamics of Nosema ceranae in honey bee midgut and host cell apoptosis.

Nosema ceranae is an intracellular microsporidian parasite that infects epithelial cells of the honey bee (Apis mellifera) midgut. Previous studies have shown that Nosema may alter cell renewal and apoptosis in honey bees. We found that the amount of apoptotic cells progressively declines from the anterior towards posterior regions of the midgut in Nosema-infected sensitive bees. There was no such pattern in the infected Nosema tolerant honey bees and controls. These data provide additional evidence that N. ceranae appears to alter apoptosis in its host cells for its own advantage.

Trueness and precision of the real-time RT-PCR method for quantifying the chronic bee paralysis virus genome in bee homogenates evaluated by a comparative inter-laboratory study.

The Chronic bee paralysis virus (CBPV) is the aetiological agent of chronic bee paralysis, a contagious disease associated with nervous disorders in adult honeybees leading to massive mortalities in front of the hives. Some of the clinical signs frequently reported, such as trembling, may be confused with intoxication syndromes. Therefore, laboratory diagnosis using real-time PCR to quantify CBPV loads is used to confirm disease. Clinical signs of chronic paralysis are usually associated with viral loads higher than 108 copies of CBPV genome copies per bee (8 log10 CBPV/bee). This threshold is used by the European Union Reference Laboratory for Bee Health to diagnose the disease. In 2015, the accuracy of measurements of three CBPV loads (5, 8 and 9 log10 CBPV/bee) was assessed through an inter-laboratory study. Twenty-one participants, including 16 European National Reference Laboratories, received 13 homogenates of CBPV-infected bees adjusted to the three loads. Participants were requested to use the method usually employed for routine diagnosis. The quantitative results (n=270) were analysed according to international standards NF ISO 13528 (2015) and NF ISO 5725-2 (1994). The standard deviations of measurement reproducibility (SR) were 0.83, 1.06 and 1.16 at viral loads 5, 8 and 9 log10 CBPV/bee, respectively. The inter-laboratory confidence of viral quantification (+/- 1.96SR) at the diagnostic threshold (8 log10 CBPV/bee) was+/- 2.08 log10 CBPV/bee. These results highlight the need to take into account the confidence of measurements in epidemiological studies using results from different laboratories. Considering this confidence, viral loads over 6 log10 CBPV/bee may be considered to indicate probable cases of chronic paralysis.

Differential proteomics reveals novel insights into Nosema-honey bee interactions.

Host manipulation is a common strategy by parasites to reduce host defense responses, enhance development, host exploitation, reproduction and, ultimately, transmission success. As these parasitic modifications can reduce host fitness, increased selection pressure may result in reciprocal adaptations of the host. Whereas the majority of studies on host manipulation have explored resistance against parasites (i.e. ability to prevent or limit an infection), data describing tolerance mechanisms (i.e. ability to limit harm of an infection) are scarce. By comparing differential protein abundance, we provide evidence of host-parasite interactions in the midgut proteomes of N. ceranae-infected and uninfected honey bees from both Nosema-tolerant and Nosema-sensitive lineages. We identified 16 proteins out of 661 protein spots that were differentially abundant between experimental groups. In general, infections of Nosema resulted in an up-regulation of the bee's energy metabolism. Additionally, we identified 8 proteins that were differentially abundant between tolerant and sensitive honey bees regardless of the Nosema infection. Those proteins were linked to metabolism, response to oxidative stress and apoptosis. In addition to bee proteins, we also identified 3 Nosema ceranae proteins. Interestingly, abundance of two of these Nosema proteins were significantly higher in infected Nosema-sensitive honeybees relative to the infected Nosema-tolerant lineage. This may provide a novel candidate for studying the molecular interplay between N. ceranae and its honey bee host in more detail.

Label-Free Differential Proteomics and Quantification of Exoenzymes from Isolates of the Entomopathogenic Fungus Beauveria bassiana.

Beauveria bassiana is an entomopathogenic fungus that grows both in vivo and in vitro. In vivo it can colonize live insect hosts, and tissue digestion occurs by secreted hydrolytic exoenzymes. It can also colonize dead insect tissue provided this is free from competing microorganisms. Depending on whether the host is alive or dead the expression (quality/quantity) of the exoenzymes may vary. We have grown several isolates of B. bassiana in shaking flasks for 120 h at 25 °C in order to evaluate the maximal exoenzyme production using two diet regimes. As sole carbon, nitrogen, and phosphate sources we used 1% shrimp chitin and either 0.5% w/v of dead intact American cockroach (Periplaneta americana) or their isolated cuticles. This is the first report of a differential proteomics of B. bassiana exoenzymes performed by label-free nano-LC MS/MS. Total proteolytic enzyme activity was mainly due to Pr1A or Pr1B depending on the isolate and the diet regime. The most differentially secreted enzymes were: the cuticle-degrading subtilisin Pr1A, GH13 alpha-glycosidase, glucan endo-1,3-beta-glucosidase, subtilisin-like proteinase Spm1, lipase 1, beta-1,3 exoglucanase, and endo-1,3-beta-glucosidase. Among the B. bassiana isolates analyzed, Bb 678 and Bb BG were the most active in Pr1A secretion.

Deformed wing virus can be transmitted during natural mating in honey bees and infect the queens.

Deformed wing virus is an important contributor to honey bee colony losses. Frequently queen failure is reported as a cause for colony loss. Here we examine whether sexual transmission during multiple matings of queens is a possible way of virus infection in queens. In an environment with high prevalence of deformed wing virus, queens (n = 30) were trapped upon their return from natural mating flights. The last drone's endophallus (n = 29), if present, was removed from the mated queens for deformed wing virus quantification, leading to the detection of high-level infection in 3 endophalli. After oviposition, viral quantification revealed that seven of the 30 queens had high-level deformed wing virus infections, in all tissues, including the semen stored in the spermathecae. Two groups of either unmated queens (n = 8) with induced egg laying, or queens (n = 12) mated in isolation with drones showing comparatively low deformed wing virus infections served as control. None of the control queens exhibited high-level viral infections. Our results demonstrate that deformed wing virus infected drones are competitive to mate and able to transmit the virus along with semen, which occasionally leads to queen infections. Virus transmission to queens during mating may be common and can contribute noticeably to queen failure.

Nosema spp. infections cause no energetic stress in tolerant honeybees.

Host-pathogen coevolution leads to reciprocal adaptations, allowing pathogens to increase host exploitation or hosts to minimise costs of infection. As pathogen resistance is often associated with considerable costs, tolerance may be an evolutionary alternative. Here, we examined the effect of two closely related and highly host dependent intracellular gut pathogens, Nosema apis and Nosema ceranae, on the energetic state in Nosema tolerant and sensitive honeybees facing the infection. We quantified the three major haemolymph carbohydrates fructose, glucose, and trehalose using high-performance liquid chromatography (HPLC) as a measure for host energetic state. Trehalose levels in the haemolymph were negatively associated with N. apis infection intensity and with N. ceranae infection regardless of the infection intensity in sensitive honeybees. Nevertheless, there was no such association in Nosema spp. infected tolerant honeybees. These findings suggest that energy availability in tolerant honeybees was not compromised by the infection. This result obtained at the individual level may also have implications at the colony level where workers in spite of a Nosema infection can still perform as well as healthy bees, maintaining colony efficiency and productivity.

Parasite infection accelerates age polyethism in young honey bees.

Honey bees (Apis mellifera) are important pollinators and their health is threatened worldwide by persistent exposure to a wide range of factors including pesticides, poor nutrition, and pathogens. Nosema ceranae is a ubiquitous microsporidian associated with high colony mortality. We used lab micro-colonies of honey bees and video analyses to track the effects of N. ceranae infection and exposure on a range of individual and social behaviours in young adult bees. We provide detailed data showing that N. ceranae infection significantly accelerated the age polyethism of young bees, causing them to exhibit behaviours typical of older bees. Bees with high N. ceranae spore counts had significantly increased walking rates and decreased attraction to queen mandibular pheromone. Infected bees also exhibited higher rates of trophallaxis (food exchange), potentially reflecting parasite manipulation to increase colony infection. However, reduction in queen contacts could help bees limit the spread of infection. Such accelerated age polyethism may provide a form of behavioural immunity, particularly if it is elicited by a wide variety of pathogens.

Four Categories of Viral Infection Describe the Health Status of Honey Bee Colonies.

Honey bee virus prevalence data are an essential prerequisite for managing epidemic events in a population. A survey study was carried out for seven viruses in colonies representing a healthy Danish honey bee population. In addition, colonies from apiaries with high level Varroa infestation or high level of winter mortality were also surveyed. Results from RT-qPCR showed a considerable difference of virus levels between healthy and sick colonies. In the group of healthy colonies, no virus was detected in 36% of cases, while at least one virus was found in each of the sick colonies. Virus titers varied among the samples, and multiple virus infections were common in both groups with a high prevalence of Sacbrood virus (SBV), Black queen cell virus (BQCV) and Deformed wing virus (DWV). Based on the distribution of virus titers, we established four categories of infection: samples free of virus (C = 0), samples with low virus titer (estimated number of virus copies 0 < C < 103), samples with medium virus titer (103 ≤ C < 107) and samples with high virus titer (C ≥ 107). This allowed us to statistically compare virus levels in healthy and sick colonies. Using categories to communicate virus diagnosis results to beekeepers may help them to reach an informed decision on management strategies to prevent further spread of viruses among colonies.