Honey bee viruses are serious pathogens that can cause poor colony health and productivity. We analyzed a multi-year longitudinal dataset of abundances of nine honey bee viruses (deformed wing virus A, deformed wing virus B, black queen cell virus, sacbrood virus, Lake Sinai virus, Kashmir bee virus, acute bee paralysis virus, chronic bee paralysis virus, and Israeli acute paralysis virus) in colonies located across Canada to describe broad trends in virus intensity and occurrence among regions and years. We also tested climatic variables (temperature, wind speed, and precipitation) as predictors in an effort to understand possible drivers underlying seasonal patterns in viral prevalence. Temperature was a significant positive predictor of the total number of viruses per sample, which was highest in British Columbia (mean = 5.0). Lake Sinai virus (LSV) was the most prevalent overall (at 89%) and had the highest infection intensity, at an average of 3.9 x 108 copies per bee. Acute bee paralysis virus was the least prevalent virus (at 4.7%) and had the lowest infection intensity (1.9 x 105 copies per bee). Surprisingly, including Varroa abundance as a covariate did not significantly improve model fit for any virus. All viruses, except Kashmir bee virus, varied by region, and one or more climatic variables were significant predictors for six of the nine viruses. Although climatic effects were often inconsistent among individual viruses, we show that climatic variables can be better predictors of virus intensity and occurrence than Varroa mite abundance, at least when infestation rates are low.

Honey bees play a major role in crop pollination but have experienced declining health throughout most of the globe. Despite decades of research on key honey bee stressors (e.g., parasitic Varroa destructor mites and viruses), researchers cannot fully explain or predict colony mortality, potentially because it is caused by exposure to multiple interacting stressors in the field. Understanding which honey bee stressors co-occur and have the potential to interact is therefore of profound importance. Here, we used the emerging field of systems theory to characterize the stressor networks found in honey bee colonies after they were placed in fields containing economically valuable crops across Canada. Honey bee stressor networks were often highly complex, with hundreds of potential interactions between stressors. Their placement in crops for the pollination season generally exposed colonies to more complex stressor networks, with an average of 23 stressors and 307 interactions. We discovered that the most influential stressors in a network—those that substantively impacted network architecture—are not currently addressed by beekeepers. Finally, the stressor networks showed substantial divergence among crop systems from different regions, which is consistent with the knowledge that some crops (e.g., highbush blueberry) are traditionally riskier to honey bees than others. Our approach sheds light on the stressor networks that honey bees encounter in the field and underscores the importance of considering interactions among stressors. Clearly, addressing and managing these issues will require solutions that are tailored to specific crops and regions and their associated stressor networks.

Recent declines in the health of honey bee colonies used for crop pollination pose a considerable threat to global food security. Foraging by honey bee workers represents the primary route of exposure to a plethora of toxins and pathogens known to affect bee health, but it remains unclear how foraging preferences impact colony-level patterns of stressor exposure. Resolving this knowledge gap is crucial for enhancing the health of honey bees and the agricultural systems that rely on them for pollination. To address this, we carried out a national-scale experiment encompassing 456 Canadian honey bee colonies to first characterize pollen foraging preferences in relation to major crops and then explore how foraging behavior influences patterns of stressor exposure. We used a metagenetic approach to quantify honey bee dietary breadth and found that bees display distinct foraging preferences that vary substantially relative to crop type and proximity, and the breadth of foraging interactions can be used to predict the abundance and diversity of stressors a colony is exposed to. Foraging on diverse plant communities was associated with increased exposure to pathogens, while the opposite was associated with increased exposure to xenobiotics. Our work provides the first large-scale empirical evidence that pollen foraging behavior plays an influential role in determining exposure to dichotomous stressor syndromes in honey bees.

Improving our understanding of how climate influences honey bee parasites and pathogens is critical as weather patterns continue to shift under climate change. While the prevalence of diseases vary according to regional and seasonal patterns, the influence of specific climatic predictors has rarely been formally assessed. To address this gap, we analyzed how occurrence and intensity of three prominent honey bee disease agents (Varroa destructor ― hereon Varroa ― Melissococcus plutonius, and Vairimorpha spp.) varied according to regional, temporal, and climatic factors in honey bee colonies across five Canadian provinces that were sampled at three time points. We found strong regional effects for all disease agents, with consistently high Varroa intensity and infestation probabilities and high M. plutonius infection probabilities in British Columbia, and year-dependent regional patterns of Vairimorpha spp. spore counts. Increasing wind speed and precipitation were linked to lower Varroa infestation probabilities, whereas warmer temperatures were linked to higher infestation probabilities. Analysis of an independent dataset shows that these trends for Varroa are consistent within a similar date range, but temperature is the strongest climatic predictor of season-long patterns. Vairimorpha spp. intensity decreased over the course of the summer, with the lowest spore counts found at later dates when temperatures were warm. Vairimorpha spp. intensity increased with wind speed and precipitation, consistent with inclement weather limiting defecation flights. Probability of M. plutonius infection generally increased across the spring and summer, and was also positively associated with inclement weather. These data contribute to building a larger dataset of honey bee disease agent occurrence that is needed in order to predict how epidemiology may change in our future climate.

Abstract Honey bees are efficient pollinators of flowering plants, aiding in the plant reproductive cycle and acting as vehicles for evolutionary processes. Their role as agents of selection and drivers of gene flow is instrumental to the structure of plant populations, but historically, our understanding of their influence has been limited to predominantly insect‐dispersed flowering species. Recent metagenetic work has provided evidence that honey bees also forage on pollen from anemophilous species, suggesting that their role as vectors for transmission of plant genetic material is not confined to groups designated as entomophilous, and leading us to ask: could honey bees act as dispersal agents for non‐flowering plant taxa? Using an extensive pollen metabarcoding dataset from Canada, we discovered that honey bees may serve as dispersal agents for an array of sporophytes (Anchistea, Claytosmunda, Dryopteris, Osmunda, Osmundastrum, Equisetum) and bryophytes (Funaria, Orthotrichum, Sphagnum, Ulota). Our findings also suggest that honey bees may occasionally act as vectors for the dispersal of aquatic phototrophs, specifically Coccomyxa and Protosiphon, species of green algae. Our work has shed light on the broad resource‐access patterns that guide plant‐pollinator interactions and suggests that bees could act as vectors of gene flow, and potentially even agents of selection, across Plantae.

Alternatives to the antibiotic fumagillin for the control of Nosema ceranae, a gut parasite of the honey bee, are needed. The prebiotics eugenol, chitosan, and naringenin and the probiotic Protexin® (Enterococcus faecium) provided in sugar syrup or protein patty either in spring or fall were evaluated for their effects on N. ceranae infection, colony population, honey yield and winter survivorship using field colonies. In the first year, spring treatments with eugenol, naringenin, and Protexin® significantly reduced N. ceranae infection and increased honey production, while Protexin® also increased adult bee populations and chitosan was ineffective. Fall treatments increased survivorship and decreased N. ceranae infection the following spring. In the second year, selected compounds were further tested with a larger number of colonies per treatment and only protein patty used in the spring and sugar syrup in the fall. Protexin® and naringenin significantly decreased N. ceranae infections and increased the population of adult bees after spring treatment, but did not affect honey yields. There were no differences between treatments for colony winter mortality, but surviving colonies that had been treated with Protexin® and naringenin were significantly more populated and had lower N. ceranae spore counts than control, non-treated colonies. Protexin® and naringenin were the most promising candidates for controlling N. ceranae and promoting honey bee populations, warranting further investigation. Future research should investigate the optimal colony dose and treatment frequency to maximize colony health.

Nosema ceranae is a microsporidian fungus that parasitizes the midgut epithelial cells of honey bees, Apis mellifera. Due to the role that midgut microorganisms play in bee health and immunity, food supplementation with prebiotics and probiotics may assist in the control of N. ceranae. The dietary fiber prebiotics acacia gum, inulin, and fructooligosaccharides, as well as the commercial probiotics Vetafarm Probotic, Protexin Concentrate single-strain (Enterococcus faecium), and Protexin Concentrate multi-strain (Lactobacillus acidophilus, L. plantarum, L. rhamnosus, L. delbrueckii, Bifidobacterium bifidum, Streptococcus salivarius, and E. faecium) were tested for their effect on N. ceranae spore loads and honey bee survivorship. Bees kept in cages were inoculated with N. ceranae spores and single-dose treatments were administered in sugar syrup. Acacia gum caused the greatest reduction in N. ceranae spore numbers (67%) but also significantly increased bee mortality (62.2%). However, Protexin Concentrate single-strain gave similarly reduced spore numbers (59%) without affecting the mortality. In a second experiment, multiple doses of the probiotics revealed significantly reduced spore numbers with 2.50 mg/mL Vetafarm Probotic, and 0.25, 1.25, and 2.50 mg/mL Protexin Concentrate single-strain. Mortality was also significantly reduced with 1.25 mg/mL Protexin Concentrate single-strain. N. ceranae-inoculated bees fed 3.75 mg/mL Vetafarm Probotic had higher survival than N. ceranae-inoculated bees, which was similar to that of non-inoculated bees, while N. ceranae-inoculated bees fed 2.50 mg/mL Protexin Concentrate single-strain, had significantly higher survival than both N. ceranae-inoculated and non-inoculated bees. Protexin Concentrate single-strain is promising as it can reduce N. ceranae proliferation and increase bee survivorship of infected bees, even compared to healthy, non-infected bees.

Nosema ceranae is a microsporidian parasite that causes nosemosis in the honey bee (Apis mellifera). As alternatives to the antibiotic fumagillin, ten nutraceuticals (oregano oil, thymol, carvacrol, trans-cinnmaldehyde, tetrahydrocurcumin, sulforaphane, naringenin, embelin, allyl sulfide, hydroxytyrosol) and two immuno-stimulatory compounds (chitosan, poly I:C) were examined for controlling N. ceranae infections. Caged bees were inoculated with N. ceranae spores, and treatments were administered in sugar syrup. Only two compounds did not significantly reduce N. ceranae spore counts compared to the infected positive control, but the most effective were sulforaphane from cruciferous vegetables, carvacrol from oregano oil, and naringenin from citrus fruit. When tested at several concentrations, the highest sulforaphane concentration reduced spore counts by 100%, but also caused 100% bee mortality. For carvacrol, the maximum reduction in spore counts was 57% with an intermediate concentration and the maximum bee mortality was 23% with the highest concentration. For naringenin, the maximum reduction in spore counts was 64% with the highest concentration, and the maximum bee mortality was only 15% with an intermediate concentration. In the longevity experiment, naringenin-fed bees lived as long as Nosema-free control bees, both of which lived significantly longer than infected positive control bees. While its antimicrobial properties may be promising, reducing sulforaphane toxicity to bees is necessary before it can be considered as a candidate for controlling N. ceranae. Although further work on formulation is needed with naringenin, its effect on extending longevity in infected bees may give it an additional value as a potential additive for bee feed in honey bee colonies.

After two years of bidirectional selection for low and high rates of Varroa destructor population growth (LVG and HVG, respectively) in honey bee (Apis mellifera) colonies in Ontario, Canada, significant differences between the two genotypes were observed. LVG colonies had V. destructor population increases over the summer of 1.7 fold compared to 9.6 fold for HVG colonies by Generation 2. Additionally, HVG colonies had significantly higher mite infestation rates in adult bees compared to LVG colonies for both selected generations. DWV prevalence and levels were significantly higher in HVG colonies than in LVG colonies in Generation 1 but not in Generation 2. Winter mortality rates of Generation 1 colonies were significantly different at 26% and 14% for the HVG and LVG genotypes, respectively. The results of this study thus far indicate that selection for LVG may result in colonies with lower V. destructor infestation rates, lower prevalence, and levels of DWV and higher colony winter survivorship. Future work will focus on determining what mechanisms are responsible for the genotypic differences, estimating genetic parameters, and molecular analyses of the genotypes to identify candidate genes associated with resistance to V. destructor and DWV that could potentially be used for marker-assisted selection.