British Columbia beekeeping revenues and costs: survey data and profit modeling.

British Columbia beekeepers, like many beekeepers around the world, are currently facing declines in honey bee health and high overwinter colony losses. To better understand the economics and the cycle of yearly colony loss and replacement of this critical agricultural industry, we collected and analyzed survey data on beekeeping costs and returns. Forty British Columbia beekeepers provided details about revenue sources, variable costs, capital costs, and investments. Ten surveyed beekeepers managed between 1 and 9 colonies, 10 managed between 10 and 39 colonies, 9 managed between 40 and 100 colonies, 5 managed between 101 and 299 colonies, 3 managed between 300 and 699 colonies, and 3 managed 700 colonies or more. The data was used to calculate beekeeping profit and to parameterize a model that explores the economic impact of colony loss rates and replacement strategies. Survey results show that when the data is aggregated, revenues exceed costs for beekeeping operations in British Columbia with a per colony profit of $56.92 or $0.87 per pound of honey produced. Surveyed operations with fewer than 100 colonies have negative profits, while operations with 100-299 colonies have positive profits. Surveyed operations in the Cariboo, North Coast, and Okanagan regions have the highest profits while surveyed operations in the Peace region have the lowest profits. Profit modeling shows that replacing losses with packages generates lower profit than replacing losses with split colonies. Our modeling shows that operations that diversify their revenue to include bee sales and commercial pollination accrue higher profits and can withstand higher winter loss rates.

Honey bee protein atlas at organ-level resolution.

Genome sequencing has provided us with gene lists but cannot tell us where and how their encoded products work together to support life. Complex organisms rely on differential expression of subsets of genes/proteins in organs and tissues, and, in concert, evolved to their present state as they function together to improve an organism's overall reproductive fitness. Proteomics studies of individual organs help us understand their basic functions, but this reductionist approach misses the larger context of the whole organism. This problem could be circumvented if all the organs in an organism were comprehensively studied by the same methodology and analyzed together. Using honey bees (Apis mellifera L.) as a model system, we report here an initial whole proteome of a complex organism, measuring 29 different organ/tissue types among the three honey bee castes: queen, drone, and worker. The data reveal that, e.g., workers have a heightened capacity to deal with environmental toxins and queens have a far more robust pheromone detection system than their nestmates. The data also suggest that workers altruistically sacrifice not only their own reproductive capacity but also their immune potential in favor of their queen. Finally, organ-level resolution of protein expression offers a systematic insight into how organs may have developed.

Higher prevalence of sacbrood virus in Apis mellifera (Hymenoptera: Apidae) colonies after pollinating highbush blueberries.

Highbush blueberry pollination depends on managed honey bees (Apis mellifera) L. for adequate fruit sets; however, beekeepers have raised concerns about the poor health of colonies after pollinating this crop. Postulated causes include agrochemical exposure, nutritional deficits, and interactions with parasites and pathogens, particularly Melisococcus plutonius [(ex. White) Bailey and Collins, Lactobacillales: Enterococcaceae], the causal agent of European foulbrood disease, but other pathogens could be involved. To broadly investigate common honey bee pathogens in relation to blueberry pollination, we sampled adult honey bees from colonies at time points corresponding to before (t1), during (t2), at the end (t3), and after (t4) highbush blueberry pollination in British Columbia, Canada, across 2 years (2020 and 2021). Nine viruses, as well as M. plutonius, Vairimorpha ceranae, and V. apis [Tokarev et al., Microsporidia: Nosematidae; formerly Nosema ceranae (Fries et al.) and N. apis (Zander)], were detected by PCR and compared among colonies located near and far from blueberry fields. We found a significant interactive effect of time and blueberry proximity on the multivariate pathogen community, mainly due to differences at t4 (corresponding to ~6 wk after the beginning of the pollination period). Post hoc comparisons of pathogens in near and far groups at t4 showed that detections of sacbrood virus (SBV), which was significantly higher in the near group, not M. plutonius, was the primary driver. Further research is needed to determine if the association of SBV with highbush blueberry pollination is contributing to the health decline that beekeepers observe after pollinating this crop.

Honey bee stressor networks are complex and dependent on crop and region.

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.

Effects of dialkoxybenzenes against Varroa destructor and identification of 1-allyloxy-4-propoxybenzene as a promising acaricide candidate.

The honey bee is responsible for pollination of a large proportion of crop plants, but the health of honey bee populations has been challenged by the parasitic mite Varroa destructor. Mite infestation is the main cause of colony losses during the winter months, which causes significant economic challenges in apiculture. Treatments have been developed to control the spread of varroa. However, many of these treatments are no longer effective due to acaricide resistance. In a search of varroa-active compounds, we tested the effect of dialkoxybenzenes on the mite. A structure-activity relationship revealed that 1-allyloxy-4-propoxybenzene is most active of a series of dialkoxybenzenes tested. We found that three compounds (1-allyloxy-4-propoxybenzene, 1,4-diallyloxybenzene and 1,4-dipropoxybenzene) cause paralysis and death of adult varroa mites, whereas the previously discovered compound, 1,3-diethoxybenzene, which alters host choice of adult mites in certain conditions, did not cause paralysis. Since paralysis can be caused by inhibition of acetylcholinesterase (AChE), a ubiquitous enzyme in the nervous system of animals, we tested dialkoxybenzenes on human, honey bee and varroa AChE. These tests revealed that 1-allyloxy-4-propoxybenzene had no effects on AChE, which leads us to conclude that 1-allyloxy-4-propoxybenzene does not exert its paralytic effect on mites through AChE. In addition to paralysis, the most active compounds affected the ability of the mites to find and remain at the abdomen of host bees provided during assays. A test of 1-allyloxy-4-propoxybenzene in the field, during the autumn of 2019 in two locations, showed that this compound has promise in the treatment of varroa infestations.

Impacts of COVID-19 on Canadian Beekeeping: Survey Results and a Profitability Analysis.

To gauge the impact of COVID-19 on the Canadian beekeeping sector, we conducted a survey of over 200 beekeepers in the fall of 2020. Our survey results show Canadian beekeepers faced two major challenges: 1) disrupted importation of honey bees (Hymenoptera: Apidae) (queen and bulk bees) that maintain populations; and 2) disrupted arrival of temporary foreign workers (TFWs). Disruptions in the arrival of bees and labor resulted in fewer colonies and less colony management, culminating in higher costs and lower productivity. Using the survey data, we develop a profitability analysis to estimate the impact of these disruptions on colony profit. Our results suggest that a disruption in either foreign worker or bee arrival allows beekeepers to compensate and while colony profits are lower, they remain positive. When both honey bee and foreign workers arrivals are disrupted for a beekeeper, even when the beekeeper experiences less significant colony health and cost impacts, a colony with a single pollination contract is no longer profitable, and a colony with two pollination contracts has significantly reduced profitability. As COVID-19 disruptions from 2020 and into 2021 become more significant to long-term colony health and more costly to a beekeeping operation, economic losses could threaten the industry's viability as well as the sustainability of pollination-dependent crop sectors across the country. The economic and agricultural impacts from the COVID-19 pandemic have exposed a vulnerability within Canada's beekeeping industry stemming from its dependency on imported labor and bees. Travel disruptions and border closures pose an ongoing threat to Canadian agriculture and apiculture in 2021 and highlight the need for Canada's beekeeping industry to strengthen domestic supply chains to minimize future risks.

Integrative Genomics Reveals the Genetics and Evolution of the Honey Bee's Social Immune System.

Social organisms combat pathogens through individual innate immune responses or through social immunity-behaviors among individuals that limit pathogen transmission within groups. Although we have a relatively detailed understanding of the genetics and evolution of the innate immune system of animals, we know little about social immunity. Addressing this knowledge gap is crucial for understanding how life-history traits influence immunity, and identifying if trade-offs exist between innate and social immunity. Hygienic behavior in the Western honey bee, Apis mellifera, provides an excellent model for investigating the genetics and evolution of social immunity in animals. This heritable, colony-level behavior is performed by nurse bees when they detect and remove infected or dead brood from the colony. We sequenced 125 haploid genomes from two artificially selected highly hygienic populations and a baseline unselected population. Genomic contrasts allowed us to identify a minimum of 73 genes tentatively associated with hygienic behavior. Many genes were within previously discovered QTLs associated with hygienic behavior and were predictive of hygienic behavior within the unselected population. These genes were often involved in neuronal development and sensory perception in solitary insects. We found that genes associated with hygienic behavior have evidence of positive selection within honey bees (Apis), supporting the hypothesis that social immunity contributes to fitness. Our results indicate that genes influencing developmental neurobiology and behavior in solitary insects may have been co-opted to give rise to a novel and adaptive social immune phenotype in honey bees.

A death pheromone, oleic acid, triggers hygienic behavior in honey bees (Apis mellifera L.).

Eusocial insects live in teeming societies with thousands of their kin. In this crowded environment, workers combat disease by removing or burying their dead or diseased nestmates. For honey bees, we found that hygienic brood-removal behavior is triggered by two odorants - ß-ocimene and oleic acid - which are released from brood upon freeze-killing. ß-ocimene is a co-opted pheromone that normally signals larval food-begging, whereas oleic acid is a conserved necromone across arthropod taxa. Interestingly, the odorant blend can induce hygienic behavior more consistently than either odorant alone. We suggest that the volatile ß-ocimene flags hygienic workers' attention, while oleic acid is the death cue, triggering removal. Bees with high hygienicity detect and remove brood with these odorants faster than bees with low hygienicity, and both molecules are strong ligands for hygienic behavior-associated odorant binding proteins (OBP16 and OBP18). Odorants that induce low levels of hygienic behavior, however, are weak ligands for these OBPs. We are therefore beginning to paint a picture of the molecular mechanism behind this complex behavior, using odorants associated with freeze-killed brood as a model.

Peptide biomarkers used for the selective breeding of a complex polygenic trait in honey bees.

We present a novel way to select for highly polygenic traits. For millennia, humans have used observable phenotypes to selectively breed stronger or more productive livestock and crops. Selection on genotype, using single-nucleotide polymorphisms (SNPs) and genome profiling, is also now applied broadly in livestock breeding programs; however, selection on protein/peptide or mRNA expression markers has not yet been proven useful. Here we demonstrate the utility of protein markers to select for disease-resistant hygienic behavior in the European honey bee (Apis mellifera L.). Robust, mechanistically-linked protein expression markers, by integrating cis- and trans- effects from many genomic loci, may overcome limitations of genomic markers to allow for selection. After three generations of selection, the resulting marker-selected stock outperformed an unselected benchmark stock in terms of hygienic behavior, and had improved survival when challenged with a bacterial disease or a parasitic mite, similar to bees selected using a phenotype-based assessment for this trait. This is the first demonstration of the efficacy of protein markers for industrial selective breeding in any agricultural species, plant or animal.

Honey bee (Hymenoptera: Apidae) distribution and potential for supplementary pollination in commercial tomato greenhouses during winter.

This study examined the use of honey bees, Apis mellifera L., to supplement bumble bee, Bombus spp., pollination in commercial tomato, Lycopersicon esculentum Miller, greenhouses in Western Canada. Honey bee colonies were brought into greenhouses already containing bumble bees and left for 1 wk to acclimatize. The following week, counts of honey and bumble bees foraging and flying throughout the greenhouse were conducted three times per day, and tomato flowers open during honey bee pollination were marked for later fruit harvest. The same counts and flower-marking also were done before and after the presence of honey bees to determine the background level of bumble bee pollination. Overall, tomato size was not affected by the addition of honey bees, but in one greenhouse significantly larger tomatoes were produced with honey bees present compared with bumble bees alone. In that greenhouse, honey bee foraging was greater than in the other greenhouses. Honey bees generally foraged within 100 m of their colony in all greenhouses. Our study invites further research to examine the use of honey bees with reduced levels of bumble bees, or as sole pollinators of greenhouse tomatoes. We also make specific recommendations for how honey bees can best be managed in greenhouses.

Worker honey bee ovary development: seasonal variation and the influence of larval and adult nutrition.

We examined the effect of larval and adult nutrition on worker honey bee (Apis mellifera L.) ovary development. Workers were fed high or low-pollen diets as larvae, and high or low-protein diets as adults. Workers fed low-protein diets at both life stages had the lowest levels of ovary development, followed by those fed high-protein diets as larvae and low- quality diets as adults, and then those fed diets poor in protein as larvae but high as adults. Workers fed high-protein diets at both life stages had the highest levels of ovary development. The increases in ovary development due to improved dietary protein in the larval and adult life stages were additive. Adult diet also had an effect on body mass. The results demonstrate that both carry-over of larval reserves and nutrients acquired in the adult life stage are important to ovary development in worker honey bees. Carry-over from larval development, however, appears to be less important to adult fecundity than is adult nutrition. Seasonal trends in worker ovary development and mass were examined throughout the brood rearing season. Worker ovary development was lowest in spring, highest in mid-summer, and intermediate in fall.

New components of the honey bee (Apis mellifera L.) queen retinue pheromone.

The honey bee queen produces pheromones that function in both releaser and primer roles such as attracting a retinue of workers around her, attracting drones on mating flights, preventing workers from reproducing at the individual (worker egg-laying) and colony (swarming) level, and regulating several other aspects of colony functioning. The queen mandibular pheromone (QMP), consisting of five synergistic components, is the only pheromone chemically identified in the honey bee (Apis mellifera L.) queen, but this pheromone does not fully duplicate the pheromonal activity of a full queen extract. To identify the remaining unknown compounds for retinue attraction, honey bee colonies were selectively bred to have low response to synthetic QMP and high response to a queen extract in a laboratory retinue bioassay. Workers from these colonies were then used in the bioassay to guide the isolation and identification of the remaining active components. Four new compounds were identified from several glandular sources that account for the majority of the difference in retinue attraction between synthetic QMP and queen extract: methyl (Z)-octadec-9-enoate (methyl oleate), (E)-3-(4-hydroxy-3-methoxyphenyl)-prop-2-en-1-ol (coniferyl alcohol), hexadecan-1-ol, and (Z9,Z12,Z15)-octadeca-9,12,15-trienoic acid (linolenic acid). These compounds were inactive alone or in combination, and they only elicited attraction in the presence of QMP. There was still unidentified activity remaining in the queen extract. The queen therefore produces a synergistic, multiglandular pheromone blend of at least nine compounds for retinue attraction, the most complex pheromone blend known for inducing a single behavior in any organism.