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.

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.

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.

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.

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.

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.

The speciation continuum: population structure, gene flow, and maternal ancestry in the Simulium arcticum complex (Diptera: Simuliidae).

Comparative analyses of populations at different stages of divergence can yield insights into the process of speciation. We assess population structure, gene flow, and maternal ancestry at five locations containing sympatric members of the Simulium arcticum complex at different stages of chromosome divergence. We analyze both nuclear and mitochondrial DNA markers, including 11 microsatellite loci, as well as COI, COII, cytb, and ND4 gene sequences. Simulium negativum, representing the later stages of divergence, shows both nuclear and mitochondrial differentiation when compared with allopatric and sympatric chromosomal forms, as well as both low contemporary and historical gene flow in sympatry. At intermediate stages of chromosome divergence, populations differ at nuclear, but not mitochondrial, loci in allopatry and sympatry. In one comparison of intermediate stage chromosomal forms (S. arcticum sensu stricto and S. apricarium), populations demonstrate low contemporary, but higher historical, gene flow in sympatry. In a second sympatric comparison (S. arcticum s. s. and S. brevicercum), both contemporary and historical gene flow are high. All analyses of sympatric populations at the earliest stages of chromosome divergence demonstrate panmixia; yet, some nuclear differentiation in allopatry is apparent. These findings suggest that molecular divergence is tracking chromosome divergence along a chromosomally-defined continuum of speciation in black flies.

Good species behaving badly: Non-monophyly of black fly sibling species in the Simulium arcticum complex (Diptera: Simuliidae).

Mitochondrial based phylogenetic reconstructions often show deviations from species-level monophyly. We used the Simulium arcticum species complex (Diptera: Simuliidae) as a model system for interpreting non-monophyly in light of chromosomal data supporting species status of siblings. For cytogenetic identification of morphologically indistinguishable black fly sibling species, larvae must be preserved in Carnoy's solution, a fixative known to degrade DNA. Consequently, we reconstructed phylogenetic relationships based on 12S, COII, cyt b, and ITS-1 gene sequences obtained from larvae sampled from presumed taxon-pure localities. As species composition at 'taxon-pure' sites may have changed at the time of sampling, we performed a second study that aimed to: (1) assess phylogenetic relationships among cytologically verified members of the S. arcticum species complex using COI and COII gene sequences; (2) determine whether useable genetic information could be gleaned from Carnoy's fixed specimens; and (3) determine the extent to which Carnoy's fixative degrades DNA over time. We consistently obtained genetic data from material stored in Carnoy's solution for two to three months. Genetic analysis of samples fixed in Carnoy's solution for up to six years indicates that larvae preserved for a maximum of five years can provide useable information for molecular analysis. Our preliminary and cytologically confirmed phylogenetic analyses demonstrate that mitochondrial DNA fails to resolve species-level monophyly of chromosomally distinct S. arcticum taxa. As results of analyses based on cytologically verified larvae mirror those of our preliminary study, we rule out imperfect taxonomy as the reason for species-level non-monophyly. Although we cannot confidently reject either inadequate phylogenetic information or incomplete lineage sorting as the cause of non-monophyly, the sharing of alleles between sympatric siblings suggests introgressive hybridization between taxa. We conclude that the patterns present in the S. arcticum phylogeny likely represent the initial stages of chromosome based sibling speciation.