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.

Post-invasion selection acts on standing genetic variation despite a severe founding bottleneck.

Invasive populations often have lower genetic diversity relative to the native-range populations from which they derive.1,2 Despite this, many biological invaders succeed in their new environments, in part due to rapid adaptation.3,4,5,6 Therefore, the role of genetic bottlenecks in constraining the adaptation of invaders is debated.7,8,9,10 Here, we use whole-genome resequencing of samples from a 10-year time-series dataset, representing the natural invasion of the Asian honey bee (Apis cerana) in Australia, to investigate natural selection occurring in the aftermath of a founding event. We find that Australia's A. cerana population was founded by as few as one colony, whose arrival was followed by a period of rapid population expansion associated with an increase of rare variants.11 The bottleneck resulted in a steep loss of overall genetic diversity, yet we nevertheless detected loci with signatures of positive selection during the first years post-invasion. When we investigated the origin of alleles under selection, we found that selection acted primarily on the variation introduced by founders and not on the variants that arose post-invasion by mutation. In all, our data highlight that selection on standing genetic variation can occur in the early years post-invasion, even where founding bottlenecks are severe.

Validating a multi-locus metabarcoding approach for characterizing mixed-pollen samples.

BACKGROUND: The mutualistic interaction between entomophilous plants and pollinators is fundamental to the structure of most terrestrial ecosystems. The sensitive nature of this relationship has been disrupted by anthropogenic modifications to natural landscapes, warranting development of new methods for exploring this trophic interaction. Characterizing the composition of pollen collected by pollinators, e.g. Apis mellifera, is a common means of exploring this relationship, but traditional methods of microscopic pollen assessment are laborious and limited in their scope. The development of pollen metabarcoding as a method of rapidly characterizing the abundance and diversity of pollen within mixed samples presents a new frontier for this type of work, but metabarcoding may have limitations, and validation is warranted before any suite of primers can be confidently used in a research program. We set out to evaluate the utility of an integrative approach, using a set of established primers (ITS2 and rbcL) versus melissopalynological analysis for characterizing 27 mixed-pollen samples from agricultural sites across Canada. RESULTS: Both individual markers performed well relative to melissopalynology at the family level with decreases in the strength of correlation and linear model fits at the genus level. Integrating data from both markers together via a multi-locus approach provided the best rank-based correlation between metagenetic and melissopalynological data at both the genus (ρ = 0.659; p < 0.001) and family level (ρ = 0.830; p < 0.001). Species accumulation curves indicated that, after controlling for sampling effort, melissopalynological characterization provides similar or higher species richness estimates than either marker. The higher number of plant species discovered via the metabarcoding approach simply reflects the vastly greater sampling effort in comparison to melissopalynology. CONCLUSIONS: Pollen metabarcoding performed well at characterizing the composition of mixed pollen samples relative to a traditional melissopalynological approach. Limitations to the quantitative application of this method can be addressed by adopting a multi-locus approach that integrates information from multiple markers.

Environmental metagenetics unveil novel plant-pollinator interactions.

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.

Genetics of tolerance in honeybees to the neonicotinoid clothianidin.

The effects of neonicotinoid insecticides (NNIs) on honeybee health are intensely debated, with numerous studies showing negative effects of exposure, while others report no such effects. We carried out experiments to study the genetic and molecular basis of NNI tolerance in honeybees, which may underlie the discrepancies observed in the literature. We discovered that worker survival post-exposure to an acute oral dose of clothianidin is heritable (H 2 = 37.8%). Tolerance to clothianidin was not associated with differences in the expression of detoxification enzymes in our experiments. Instead, mutations in the primary neonicotinoid detoxification genes CYP9Q1 and CYP9Q3 were strongly associated with worker survival post-clothianidin exposure. In some instances, the strong association between CYP9Q haplotypes and worker survival was associated with the protein's predicted binding affinity for clothianidin. Our findings have implications regarding future toxicological studies utilizing honeybees as a model pollinator.

Bees in the six: Determinants of bumblebee habitat quality in urban landscapes.

With growing urbanization, it is becoming increasingly important to design cities in a manner that sustains and enhances biodiversity and ecosystem services. Native bees are critical pollinators that have experienced substantive declines over the past several decades. These declines have captured the attention of the public, particularly urbanites, prompting a large interest in protecting pollinators and their habitats in cities across North America and Europe. Unfortunately, we currently lack research about specific features of urban environments that can enhance the fitness of pollinators. We carried out an intensive study of Bombus impatiens, the Common Eastern Bumblebee, in the city of Toronto (Canada's largest city), to better understand landscape parameters that provide high-quality habitat for this species and likely other generalist bees. We divided the city into 270 grid cells and sampled a large number of worker bees, which were then genotyped at twelve hypervariable microsatellite loci. The genetic data allowed us to quantify the effective number of colonies and foraging distance for bumblebees in our study area. We then asked how the city's landscape and human population demography and income are associated with the availability of high-quality habitat for B. impatiens. Several aspects of Toronto's landscape influenced colony density and foraging range. Urbanization had a clear effect on both colony density and foraging distance of workers. On the other hand, functional (i.e., not cosmetic) green space was often associated with higher quality habitats for bumblebees. Our study suggests several planning strategies to enhance habitat quality for bumblebees and other pollinators in cities.

Phenomic analysis of the honey bee pathogen-web and its dynamics on colony productivity, health and social immunity behaviors.

Many pathogens and parasites have evolved to overwhelm and suppress their host's immune system. Nevertheless, the interactive effects of these agents on colony productivity and wintering success have been relatively unexplored, particularly in large-scale phenomic studies. As a defense mechanism, honey bees have evolved remarkable social behaviors to defend against pathogen and parasite challenges, which reduce the impact of disease and improve colony health. To investigate the complex role of pathogens, parasites and social immunity behaviors in relation to colony productivity and outcomes, we extensively studied colonies at several locations across Canada for two years. In 2016 and 2017, colonies founded with 1-year-old queens of diverse genetic origin were evaluated, which represented a generalized subset of the Canadian bee population. During each experimental year (May through April), we collected phenotypic data and sampled colonies for pathogen analysis in a standardized manner. Measures included: colony size and productivity (colony weight, cluster size, honey production, and sealed brood population), social immunity traits (hygienic behavior, instantaneous mite population growth rate, and grooming behavior), as well as quantification of gut parasites (Nosema spp., and Lotmaria passim), viruses (DWV-A, DWV-B, BQCV and SBV) and external parasites (Varroa destructor). Our goal was to examine: 1) correlations between pathogens and colony phenotypes; 2) the dynamics of pathogens and parasites on colony phenotypes and productivity traits; and 3) the effects of social immunity behaviors on colony pathogen load. Our results show that colonies expressing high levels of some social immunity behaviors were associated with low levels of pathogens/parasites, including viruses, Nosema spp., and V. destructor. In addition, we determined that elevated viral and Nosema spp. levels were associated with low levels of colony productivity, and that five out of six pathogenic factors measured were negatively associated with colony size and weight in both fall and spring periods. Finally, this study also provides information about the incidence and abundance of pathogens, colony phenotypes, and further disentangles their inter-correlation, so as to better understand drivers of honey bee colony health and productivity.

Thrice out of Asia and the adaptive radiation of the western honey bee.

The origin of the western honey bee Apis mellifera has been intensely debated. Addressing this knowledge gap is essential for understanding the evolution and genetics of one of the world's most important pollinators. By analyzing 251 genomes from 18 native subspecies, we found support for an Asian origin of honey bees with at least three expansions leading to African and European lineages. The adaptive radiation of honey bees involved selection on a few genomic "hotspots." We found 145 genes with independent signatures of selection across all bee lineages, and these genes were highly associated with worker traits. Our results indicate that a core set of genes associated with worker and colony traits facilitated the adaptive radiation of honey bees across their vast distribution.

Practical Applications of Genomics in Managing Honey bee Health.

The honey bee Apis mellifera is a model organism for sociogenomics and one of the most important managed pollinators. High mortalities experienced by honey bee colonies over the past several decades are expected to have a substantive effect on crop pollination and global food security. These threats and the availability of a growing number of genomic resources for the honey bee have motivated research on how genetics and genomics can be practically applied to manage bee health. The authors review 3 such applications: (1) Certification of bee lineages using single-polymorphism markers; (2) breeding bees using marker-assisted selection; (3) diagnosing honey bee stressors using biomarkers.

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.

Conservation genomics reveals pesticide and pathogen exposure in the declining bumble bee Bombus terricola.

In recent years, many pollinators have experienced large population declines, which threaten food security and the stability of natural ecosystems. Bumble bees are particularly important because their ability to "buzz" pollinate and tolerate cooler conditions make them critical pollinators for certain plants and regions. Here, we apply a conservation genomics approach to study the vulnerable Bombus terricola. We sequenced RNA from 30 worker abdomens, 18 of which were collected from agricultural sites and 12 of which were collected from nonagricultural sites. We found transcriptional signatures associated with exposure to insecticides, with gene expression patterns suggesting that bumble bees were exposed to neonicotinoids and/or fipronil-two compounds known to negatively impact bees. We also found transcriptional signatures associated with pathogen infections. In addition to the transcriptomic analysis, we carried out a metatranscriptomic analysis and detected five pathogens in the abdomens of workers, three of which are common in managed honey bee and bumble bee colonies. Our conservation genomics study provides functional support for the role of pesticides and pathogen spillover in the decline of B. terricola. We demonstrate that conservation genomics is an invaluable tool which allows researchers to quantify the effects of multiple stressors that impact pollinator populations in the wild.

Adaptive, caste-specific changes to recombination rates in a thelytokous honeybee population.

The ability to clone oneself has clear benefits-no need for mate hunting or dilution of one's genome in offspring. It is therefore unsurprising that some populations of haplo-diploid social insects have evolved thelytokous parthenogenesis-the virgin birth of a female. But thelytokous parthenogenesis has a downside: the loss of heterozygosity (LoH) as a consequence of genetic recombination. LoH in haplo-diploid insects can be highly deleterious because female sex determination often relies on heterozygosity at sex-determining loci. The two female castes of the Cape honeybee, Apis mellifera capensis, differ in their mode of reproduction. While workers always reproduce thelytokously, queens always mate and reproduce sexually. For workers, it is important to reduce the frequency of recombination so as to not produce offspring that are homozygous. Here, we ask whether recombination rates differ between Cape workers and Cape queens that we experimentally manipulated to reproduce thelytokously. We tested our hypothesis that Cape workers have evolved mechanisms that restrain genetic recombination, whereas queens have no need for such mechanisms because they reproduce sexually. Using a combination of microsatellite genotyping and whole-genome sequencing we find that a reduction in recombination is confined to workers only.

Searching beyond the streetlight: Neonicotinoid exposure alters the neurogenomic state of worker honey bees.

Neonicotinoid insecticides have been implicated in honey bee declines, with many studies showing that sublethal exposure impacts bee behaviors such as foraging, learning, and memory. Despite the large number of ecotoxicological studies carried out to date, most focus on a handful of worker phenotypes leading to a "streetlight effect" where the a priori choice of phenotypes to measure may influence the results and conclusions arising from the studies. This bias can be overcome with the use of toxicological transcriptomics, where changes in gene expression can provide a more objective view of how pesticides alter animal traits. Here, we used RNA sequencing to examine the changes in neurogenomic states of nurse and forager honey bees that were naturally exposed to neonicotinoids in the field and artificially exposed to neonicotinoids in a controlled experiment. We found that neonicotinoid exposure influenced the neurogenomic state of foragers and nurses in different ways; foragers experienced shifts in expression of genes involved in cognition and development, while nurses experienced shifts in expression of genes involved in metabolism. Our study suggests that neonicotinoids influence nurse and forager bees in a different manner. We also found no to minimal overlap in the differentially expressed genes in our study and in previously published studies, which might help reconcile the seemingly contradictory results often reported in the neonicotinoid literature.

Eusociality influences the strength of negative selection on insect genomes.

While much of the focus of sociobiology concerns identifying genomic changes that influence social behaviour, we know little about the consequences of social behaviour on genome evolution. It has been hypothesized that social evolution can influence the strength of negative selection via two mechanisms. First, division of labour can influence the efficiency of negative selection in a caste-specific manner; indirect negative selection on worker traits is theoretically expected to be weaker than direct selection on queen traits. Second, increasing social complexity is expected to lead to relaxed negative selection because of its influence on effective population size. We tested these two hypotheses by estimating the strength of negative selection in honeybees, bumblebees, paper wasps, fire ants and six other insects that span the range of social complexity. We found no consistent evidence that negative selection was significantly stronger on queen-biased genes relative to worker-biased genes. However, we found strong evidence that increased social complexity reduced the efficiency of negative selection. Our study clearly illustrates how changes in behaviour can influence patterns of genome evolution by modulating the strength of natural selection.

Defense Response in Brazilian Honey Bees (Apis mellifera scutellata × spp.) Is Underpinned by Complex Patterns of Admixture.

In 1957, an invasive and highly defensive honey bee began to spread across Brazil. In the previous year, Brazilian researchers hoped to produce a subtropical-adapted honey bee by crossing local commercial honey bees (of European origin) with a South African honey bee subspecies (Apis mellifera scutellata; an A-lineage honey bee subspecies). The resulting cross-African hybrid honey bees (AHBs)-escaped from their enclosure and spread through the Americas. Today, AHB is the most common honey bee from Northern Argentina to the Southern United States. AHBs are much more likely to sting nest intruders than managed European-derived honey bee colonies. Previous studies have explored how genetic variation contributes to differences in defense response between European-derived honey bee and AHB. Although this work demonstrated very strong genetic effects on defense response, they have yet to pinpoint which genes influence variation in defense response within AHBs, specifically. We quantified defense response for 116 colonies in Brazil and performed pooled sequencing on the most phenotypically divergent samples. We identified 65 loci containing 322 genes that were significantly associated with defense response. Loci were strongly associated with metabolic function, consistent with previous functional genomic analyses of this phenotype. Additionally, defense-associated loci had nonrandom and unexpected patterns of admixture. Defense response was not simply the product of more A-lineage honey bee ancestry as previously assumed, but rather an interaction between A-lineage and European alleles. Our results suggest that a combination of A-lineage and European alleles play roles in defensive behavior in AHBs.

Honey Bee Queen Production: Canadian Costing Case Study and Profitability Analysis.

The decline in managed honey bee (Hymenoptera: Apidae) colony health worldwide has had a significant impact on the beekeeping industry. To mitigate colony losses, beekeepers in Canada and around the world introduce queens into replacement colonies; however, Canada's short queen rearing season has historically limited the production of early season queens. As a result, Canadian beekeepers rely on the importation of foreign bees, particularly queens from warmer climates. Importing a large proportion of (often mal-adapted) queens each year creates a dependency on foreign bee sources, putting beekeeping, and pollination sectors at risk in the event of border closures, transportation issues, and other restrictions as is currently happening due to the 2020 Covid-19 pandemic. Although traditional Canadian queen production is unable to fully meet early season demand, increasing domestic queen production to meet mid- and later season demand would reduce Canada's dependency. As well, on-going studies exploring the potential for overwintering queens in Canada may offer a strategy to have early season domestic queens available. Increasing the local supply of queens could provide Canadian beekeepers, farmers, and consumers with a greater level of agricultural stability and food security. Our study is the first rigorous analysis of the economic feasibility of queen production. We present the costs of queen production for three Canadian operations over two years. Our results show that it can be profitable for a beekeeping operation in Canada to produce queen cells and mated queens and could be one viable strategy to increase the sustainability of the beekeeping industry.

Developmental plasticity shapes social traits and selection in a facultatively eusocial bee.

Developmental plasticity generates phenotypic variation, but how it contributes to evolutionary change is unclear. Phenotypes of individuals in caste-based (eusocial) societies are particularly sensitive to developmental processes, and the evolutionary origins of eusociality may be rooted in developmental plasticity of ancestral forms. We used an integrative genomics approach to evaluate the relationships among developmental plasticity, molecular evolution, and social behavior in a bee species (Megalopta genalis) that expresses flexible sociality, and thus provides a window into the factors that may have been important at the evolutionary origins of eusociality. We find that differences in social behavior are derived from genes that also regulate sex differentiation and metamorphosis. Positive selection on social traits is influenced by the function of these genes in development. We further identify evidence that social polyphenisms may become encoded in the genome via genetic changes in regulatory regions, specifically in transcription factor binding sites. Taken together, our results provide evidence that developmental plasticity provides the substrate for evolutionary novelty and shapes the selective landscape for molecular evolution in a major evolutionary innovation: Eusociality.

A Single Gene Causes Thelytokous Parthenogenesis, the Defining Feature of the Cape Honeybee Apis mellifera capensis.

In honeybees, the ability of workers to produce daughters asexually, i.e., thelytokous parthenogenesis, is restricted to a single subspecies inhabiting the Cape region of South Africa, Apis mellifera capensis. Thelytoky has unleashed new selective pressures and the evolution of traits such as social parasitism, invasiveness, and social cancer. Thelytoky arises from an abnormal meiosis that results in the fusion of two maternal pronuclei, restoring diploidy in newly laid eggs. The genetic basis underlying thelytoky is disputed. To resolve this controversy, we generated a backcross between thelytokous A. m. capensis and non-thelytokous A. m. scutellata from the neighboring population and looked for evidence of genetic markers that co-segregated with thelytokous reproduction in 49 backcross females. We found that markers associated with the gene GB45239 on chromosome 11, including non-synonymous variants, showed consistent co-segregation with thelytoky, whereas no other region did so. Alleles associated with thelytoky were present in all A. m. capensis genomes examined but were absent from all other honeybees worldwide including A. m. scutellata. GB45239 is derived in A. m. capensis and has a putative role in chromosome segregation. It is expressed in ovaries and is downregulated in thelytokous bees, likely because of polymorphisms in the promoter region. Our study reveals how mutations affecting the sequence and/or expression of a single gene can change the reproductive mode of a population.

Paternally-biased gene expression follows kin-selected predictions in female honey bee embryos.

The Kinship Theory of Genomic Imprinting (KTGI) posits that, in species where females mate with multiple males, there is selection for a male to enhance the reproductive success of his offspring at the expense of other males and his mating partner. Reciprocal crosses between honey bee subspecies show parent-of-origin effects for reproductive traits, suggesting that males modify the expression of genes related to female function in their female offspring. This effect is likely to be greater in the Cape honey bee (Apis mellifera capensis), because a male's daughters have the unique ability to produce female offspring that can develop into reproductive workers or the next queen without mating. We generated reciprocal crosses between Capensis and another subspecies and used RNA-seq to identify transcripts that are over- or underexpressed in the embryos, depending on the parental origin of the gene. As predicted, 21 genes showed expression bias towards the Capensis father's allele in colonies with a Capensis father, with no such bias in the reciprocal cross. A further six genes showed a consistent bias towards expression of the father's allele across all eight colonies examined, regardless of the direction of the cross. Consistent with predictions of the KTGI, six of the 21 genes are associated with female reproduction. No gene consistently showed overexpression of the maternal allele.