Species traits, landscape quality and floral resource overlap with honeybees determine virus transmission in plant-pollinator networks.

Emerging infectious diseases pose a threat to pollinators. Virus transmission among pollinators via flowers may be reinforced by anthropogenic land-use change and concomitant alteration of plant-pollinator interactions. Here, we examine how species' traits and roles in flower-visitation networks and landscape-scale factors drive key honeybee viruses-black queen cell virus (BQCV) and deformed wing virus-in 19 wild bee and hoverfly species, across 12 landscapes varying in pollinator-friendly (flower-rich) habitat. Viral loads were on average more than ten times higher in managed honeybees than in wild pollinators. Viral loads in wild pollinators were higher when floral resource use overlapped with honeybees, suggesting these as reservoir hosts, and increased with pollinator abundance and viral loads in honeybees. Viral prevalence decreased with the amount of pollinator-friendly habitat in a landscape, which was partly driven by reduced floral resource overlap with honeybees. Black queen cell virus loads decreased with a wild pollinator's centrality in the network and the proportion of visited dish-shaped flowers. Our findings highlight the complex interplay of resource overlap with honeybees, species traits and roles in flower-visitation networks and flower-rich pollinator habitat shaping virus transmission.

Inhibition of Renin Expression Is Regulated by an Epigenetic Switch From an Active to a Poised State.

BACKGROUND: Renin-expressing cells are myoendocrine cells crucial for the maintenance of homeostasis. Renin is regulated by cAMP, p300 (histone acetyltransferase p300)/CBP (CREB-binding protein), and Brd4 (bromodomain-containing protein 4) proteins and associated pathways. However, the specific regulatory changes that occur following inhibition of these pathways are not clear. METHODS: We treated As4.1 cells (tumoral cells derived from mouse juxtaglomerular cells that constitutively express renin) with 3 inhibitors that target different factors required for renin transcription: H-89-dihydrochloride, PKA (protein kinase A) inhibitor; JQ1, Brd4 bromodomain inhibitor; and A-485, p300/CBP inhibitor. We performed assay for transposase-accessible chromatin with sequencing (ATAC-seq), single-cell RNA sequencing, cleavage under targets and tagmentation (CUT&Tag), and chromatin immunoprecipitation sequencing for H3K27ac (acetylation of lysine 27 of the histone H3 protein) and p300 binding on biological replicates of treated and control As4.1 cells. RESULTS: In response to each inhibitor, Ren1 expression was significantly reduced and reversible upon washout. Chromatin accessibility at the Ren1 locus did not markedly change but was globally reduced at distal elements. Inhibition of PKA led to significant reductions in H3K27ac and p300 binding specifically within the Ren1 super-enhancer region. Further, we identified enriched TF (transcription factor) motifs shared across each inhibitory treatment. Finally, we identified a set of 9 genes with putative roles across each of the 3 renin regulatory pathways and observed that each displayed differentially accessible chromatin, gene expression, H3K27ac, and p300 binding at their respective loci. CONCLUSIONS: Inhibition of renin expression in cells that constitutively synthesize and release renin is regulated by an epigenetic switch from an active to poised state associated with decreased cell-cell communication and an epithelial-mesenchymal transition. This work highlights and helps define the factors necessary for renin cells to alternate between myoendocrine and contractile phenotypes.

Eco-evolutionary processes shaping floral nectar sugar composition.

Floral nectar sugar composition is assumed to reflect the nutritional demands and foraging behaviour of pollinators, but the relative contributions of evolutionary and abiotic factors to nectar sugar composition remain largely unknown across the angiosperms. We compiled a comprehensive dataset on nectar sugar composition for 414 insect-pollinated plant species across central Europe, along with phylogeny, paleoclimate, flower morphology, and pollinator dietary demands, to disentangle their relative effects. We found that phylogeny was strongly related with nectar sucrose content, which increased with the phylogenetic age of plant families, but even more strongly with historic global surface temperature. Nectar sugar composition was also defined by floral morphology, though it was not related to our functional measure of pollinator dietary demands. However, specialist pollinators of current plant-pollinator networks predominantly visited plant species with sucrose-rich nectar. Our results suggest that both physiological mechanisms related to plant water balance and evolutionary effects related to paleoclimatic changes have shaped floral nectar sugar composition during the radiation and specialisation of plants and pollinators. As a consequence, the high velocity of current climate change may affect plant-pollinator interaction networks due to a conflicting combination of immediate physiological responses and phylogenetic conservatism.

Allele Frequencies of Genetic Variants Associated with Varroa Drone Brood Resistance (DBR) in Apis mellifera Subspecies across the European Continent.

Implementation of marker-assisted selection (MAS) in modern beekeeping would improve sustainability, especially in breeding programs aiming for resilience against the parasitic mite Varroa destructor. Selecting honey bee colonies for natural resistance traits, such as brood-intrinsic suppression of varroa mite reproduction, reduces the use of chemical acaricides while respecting local adaptation. In 2019, eight genomic variants associated with varroa non-reproduction in drone brood were discovered in a single colony from the Amsterdam Water Dune population in the Netherlands. Recently, a new study tested the applicability of these eight genetic variants for the same phenotype on a population-wide scale in Flanders, Belgium. As the properties of some variants varied between the two studies, one hypothesized that the difference in genetic ancestry of the sampled colonies may underly these contribution shifts. In order to frame this, we determined the allele frequencies of the eight genetic variants in more than 360 Apis mellifera colonies across the European continent and found that variant type allele frequencies of these variants are primarily related to the A. mellifera subspecies or phylogenetic honey bee lineage. Our results confirm that population-specific genetic markers should always be evaluated in a new population prior to using them in MAS programs.

Host Barriers Limit Viral Spread in a Spillover Host: A Study of Deformed Wing Virus in the Bumblebee Bombus terrestris.

The transmission of pathogens from reservoir to recipient host species, termed pathogen spillover, can profoundly impact plant, animal, and public health. However, why some pathogens lead to disease emergence in a novel species while others fail to establish or do not elicit disease is often poorly understood. There is strong evidence that deformed wing virus (DWV), an (+)ssRNA virus, spills over from its reservoir host, the honeybee Apis mellifera, into the bumblebee Bombus terrestris. However, the low impact of DWV on B. terrestris in laboratory experiments suggests host barriers to virus spread in this recipient host. To investigate potential host barriers, we followed the spread of DWV genotype B (DWV-B) through a host's body using RT-PCR after experimental transmission to bumblebees in comparison to honeybees. Inoculation was per os, mimicking food-borne transmission, or by injection into the bee's haemocoel, mimicking vector-based transmission. In honeybees, DWV-B was present in both honeybee faeces and haemolymph within 3 days of inoculation per os or by injection. In contrast, DWV-B was not detected in B. terrestris haemolymph after inoculation per os, suggesting a gut barrier that hinders DWV-B's spread through the body of a B. terrestris. DWV-B was, however, detected in B. terrestris faeces after injection and feeding, albeit at a lower abundance than that observed for A. mellifera, suggesting that B. terrestris sheds less DWV-B than A. mellifera in faeces when infected. Barriers to viral spread in B. terrestris following oral infection may limit DWV's impact on this spillover host and reduce its contribution to the community epidemiology of DWV.

A new exposure protocol adapted for wild bees reveals species-specific impacts of the sulfoximine insecticide sulfoxaflor.

Wild bees are crucial pollinators of flowering plants and concerns are rising about their decline associated with pesticide use. Interspecific variation in wild bee response to pesticide exposure is expected to be related to variation in their morphology, physiology, and ecology, though there are still important knowledge gaps in its understanding. Pesticide risk assessments have largely focussed on the Western honey bee sensitivity considering it protective enough for wild bees. Recently, guidelines for Bombus terrestris and Osmia bicornis testing have been developed but are not yet implemented at a global scale in pesticide risk assessments. Here, we developed and tested a new simplified method of pesticide exposure on wild bee species collected from the field in Belgium. Enough specimens of nine species survived in a laboratory setting and were exposed to oral and topical acute doses of a sulfoximine insecticide. Our results confirm significant variability among wild bee species. We show that Osmia cornuta is more sensitive to sulfoxaflor than B. terrestris, whereas Bombus hypnorum is less sensitive. We propose hypotheses on the mechanisms explaining interspecific variations in sensitivity to pesticides. Future pesticide risk assessments of wild bees will require further refinement of protocols for their controlled housing and exposure.

Bridging the Gap between Field Experiments and Machine Learning: The EC H2020 B-GOOD Project as a Case Study towards Automated Predictive Health Monitoring of Honey Bee Colonies.

Honey bee colonies have great societal and economic importance. The main challenge that beekeepers face is keeping bee colonies healthy under ever-changing environmental conditions. In the past two decades, beekeepers that manage colonies of Western honey bees (Apis mellifera) have become increasingly concerned by the presence of parasites and pathogens affecting the bees, the reduction in pollen and nectar availability, and the colonies' exposure to pesticides, among others. Hence, beekeepers need to know the health condition of their colonies and how to keep them alive and thriving, which creates a need for a new holistic data collection method to harmonize the flow of information from various sources that can be linked at the colony level for different health determinants, such as bee colony, environmental, socioeconomic, and genetic statuses. For this purpose, we have developed and implemented the B-GOOD (Giving Beekeeping Guidance by computational-assisted Decision Making) project as a case study to categorize the colony's health condition and find a Health Status Index (HSI). Using a 3-tier setup guided by work plans and standardized protocols, we have collected data from inside the colonies (amount of brood, disease load, honey harvest, etc.) and from their environment (floral resource availability). Most of the project's data was automatically collected by the BEEP Base Sensor System. This continuous stream of data served as the basis to determine and validate an algorithm to calculate the HSI using machine learning. In this article, we share our insights on this holistic methodology and also highlight the importance of using a standardized data language to increase the compatibility between different current and future studies. We argue that the combined management of big data will be an essential building block in the development of targeted guidance for beekeepers and for the future of sustainable beekeeping.

Shift in virus composition in honeybees (Apis mellifera) following worldwide invasion by the parasitic mite and virus vector Varroa destructor.

Invasive vectors can induce dramatic changes in disease epidemiology. While viral emergence following geographical range expansion of a vector is well known, the influence a vector can have at the level of the host's pathobiome is less well understood. Taking advantage of the formerly heterogeneous spatial distribution of the ectoparasitic mite Varroa destructor that acts as potent virus vector among honeybees Apis mellifera, we investigated the impact of its recent global spread on the viral community of honeybees in a retrospective study of historical samples. We hypothesized that the vector has had an effect on the epidemiology of several bee viruses, potentially altering their transmissibility and/or virulence, and consequently their prevalence, abundance, or both. To test this, we quantified the prevalence and loads of 14 viruses from honeybee samples collected in mite-free and mite-infested populations in four independent geographical regions. The presence of the mite dramatically increased the prevalence and load of deformed wing virus, a cause of unsustainably high colony losses. In addition, several other viruses became more prevalent or were found at higher load in mite-infested areas, including viruses not known to be actively varroa-transmitted, but which may increase opportunistically in varroa-parasitized bees.

Mason bees and honey bees synergistically enhance fruit set in sweet cherry orchards.

Mason bees (Osmia spp.) are efficient fruit tree pollinators that can be encouraged to occupy and breed in artificial nesting material. In sweet cherry orchards, they are occasionally used as an alternative managed pollinator as a replacement for or in addition to honey bees (Apis mellifera). Yet, the lack of practical guidelines on management practices, for example optimal stocking rates, for both mason bee nesting material and honey bees might compromise pollination service provision. In this study, we assessed the relationship between stocking rates (honey bee hives and mason bee nesting material) and the abundance of honey bees and mason bees in 17 sweet cherry (Prunus avium) orchards in Central Germany. We furthermore performed a pollination experiment to explore the interactive effect of mason bees and honey bees on sweet cherry fruit set. In the orchards, both honey bee and mason bee abundance increased with increasing stocking rates of hives or nesting material, respectively. Honey bee abundance increased linearly with stocking rates. In contrast, mason bee abundance asymptoted at 2-3 nesting boxes per ha, beyond which more boxes resulted in little increase in visitation rate. Our pollination experiment demonstrated that orchards were pollen limited, with only 28% of insect-pollinated flowers setting fruit versus 39% of optimally hand-pollinated flowers. Honey bees and mason bees enhanced sweet cherry fruit set, but only when both were present and not when either was present alone in an orchard. Our findings demonstrate that offering nesting material for mason bees and employing honey bee hives can enhance bee abundance in sweet cherry orchards. By increasing honey bee abundance in combination with enhanced mason bee abundance, farmers can substantially boost fruit set and potentially sweet cherry yield. To enhance pollination services, farmers should consider the benefits of increasing pollinator biodiversity as an immediate benefit to improve crop yields.

Transcriptomic Responses Underlying the High Virulence of Black Queen Cell Virus and Sacbrood Virus following a Change in Their Mode of Transmission in Honey Bees (Apis mellifera).

BACKGROUND: Over the last two decades, honey bees (Apis mellifera) have suffered high rates of colony losses that have been attributed to a variety of factors, chief among which are viral pathogens, such as deformed wing virus (DWV), whose virulence has increased because of vector-based transmission by the invasive, ectoparasitic varroa mite (Varroa destructor). A shift in the experimental mode of transmission of the black queen cell virus (BQCV) and sacbrood virus (SBV) from fecal/food-oral (direct horizontal) to vector-mediated (indirect horizontal) transmission also results in high virulence and viral titers in pupal and adult honey bees. Agricultural pesticides represent another factor that acts independently or in interaction with pathogens, and they are also thought to cause colony loss. Understanding the molecular mechanisms underlying the higher virulence following a vector-based mode of transmission provides deeper insight into honey bee colony losses, as does determining whether or not host-pathogen interactions are modulated by exposure to pesticides. METHODS: Through an experimental design with controlled laboratory, we investigated the effects of the modes of transmission of BQCV and SBV (feeding vs. vector-mediated via injection) alone or in combination with chronic exposure to sublethal and field-realistic concentrations of flupyradifurone (FPF), a novel agricultural insecticide, on honey bee survival and transcription responses by using high-throughput RNA sequencing (RNA-seq) analysis. RESULTS: Co-exposure to viruses via feeding (VF) or injection (VI) and FPF insecticide had no statistically significant interactive effect on their survival compared to, respectively, VF or VI treatments alone. Transcriptomic analysis revealed a distinct difference in the gene expression profiles of bees inoculated with viruses via injection (VI) and exposed to FPF insecticide (VI+FPF). The number of differentially expressed genes (DEGs) at log2 (fold-change) > 2.0 in VI bees (136 genes) or/and VI+FPF insecticide (282 genes) was very high compared to that of VF bees (8 genes) or the VF+FPF insecticide treatment (15 genes). Of these DEGs, the expression in VI and VI+FPF bees of some immune-related genes, such as those for antimicrobial peptides, Ago2, and Dicer, was induced. In short, several genes encoding odorant binding proteins, chemosensory proteins, odor receptors, honey bee venom peptides, and vitellogenin were downregulated in VI and VI+FPF bees. CONCLUSIONS: Given the importance of these suppressed genes in honey bees' innate immunity, eicosanoid biosynthesis, and olfactory associative function, their inhibition because of the change in the mode of infection with BQCV and SBV to vector-mediated transmission (injection into haemocoel) could explain the high virulence observed in these viruses when they were experimentally injected into hosts. These changes may help explain why other viruses, such as DWV, represent such a threat to colony survival when transmitted by varroa mites.

The presence of identical deformed wing virus sequence variants in co-occurring Apis species in Northern Thailand may represent a potential epidemiological threat to native honey bees of Southeast Asia.

Widespread native honey bee species in South and East Asia (Apis cerana, Apis dorsata and Apis florea) and the imported western honey bee (Apis mellifera) share habitats and potentially also share pathogens. Chief among the threats facing A. mellifera in Europe and North America is deformed wing virus (DWV), including its two principal genotypes: A and B (DWV-A and DWV-B respectively). Though DWV-A has been recorded in Asia's native Apis species, it is not known if DWV-B, or both DWV-A and DWV-B, are currently widespread in Asia and, if so, whether viral transmission is primarily intraspecific or interspecific. This study aims to fill these knowledge gaps by (i) determining the DWV genotype in four co-occurring Apis host species using qPCR and (ii) inferring viral transmission between them using nucleotide sequences of DWV from Apis host species collected at three independent localities in Northern Thailand. We found DWV-A and -B in all four Apis species, the exotic A. mellifera and the native A. cerana, A. dorsata and A. florea. That DWV-A sequences were identical across Apis species at the same locality, with a similar pattern for DWV-B sequences, suggests that DWV's epidemiology is largely driven by ongoing interspecific transmission (spillover) of DWV across co-occurring native and exotic Apis species. Both genotypes of DWV represent a serious threat to Asia's exotic and native honey bee species.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees.

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Nationwide Screening for Bee Viruses in Apis mellifera Colonies in Egypt.

Honey bees are essential for crop and wild plant pollination. However, many countries have reported high annual colony losses caused by multiple possible stressors. Diseases, particularly those caused by viruses, are a major cause of colony losses. However, little is known about the prevalence of honey bee pathogens, particularly virus prevalence, in Egyptian honey bees. To address this shortfall, we determined the prevalence of widespread bee viruses in honey bee colonies in Egypt-whether it is affected by geography, the season, or infestation with Varroa destructor (varroa) mites. Honey bee worker samples were collected from 18 geographical regions across Egypt during two seasons: winter and summer of 2021. Three apiaries were chosen in each region, and a pooled sample of 150 worker bees was collected from five colonies in each apiary then screened by qPCR for 10 viral targets: acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), deformed wing virus (DWV) genotypes A (DWV-A), B (DWV-B) and D (Egyptian bee virus), Israeli acute paralysis virus (IAPV), Kashmir bee virus (KBV), sacbrood virus (SBV), and slow bee paralysis virus (SBPV). Our results revealed that DWV-A was the most prevalent virus, followed by BQCV and ABPV; the DWV genotype now spreading across the world, DWV-B, was not detected. There was no difference in varroa infestation rates as well as virus prevalence between winter and summer. However, colonies infected with BQCV had a significantly higher varroa count (adjusted p < 0.05) in the winter season, indicating that there is a seasonal association between the intensity of infestation by varroa and the presence of this virus. We provide data on the current virus prevalence in Egypt, which could assist in the protection of Egypt's beekeeping industry. Moreover, our study aids in the systematic assessment of the global honey bee virome by filling a knowledge gap about the prevalence of honey bee viruses in Egypt.

Characterization Factors to Assess Land Use Impacts on Pollinator Abundance in Life Cycle Assessment.

While wild pollinators play a key role in global food production, their assessment is currently missing from the most commonly used environmental impact assessment method, Life Cycle Assessment (LCA). This is mainly due to constraints in data availability and compatibility with LCA inventories. To target this gap, relative pollinator abundance estimates were obtained with the use of a Delphi assessment, during which 25 experts, covering 16 nationalities and 45 countries of expertise, provided scores for low, typical, and high expected abundance associated with 24 land use categories. Based on these estimates, this study presents a set of globally generic characterization factors (CFs) that allows translating land use into relative impacts to wild pollinator abundance. The associated uncertainty of the CFs is presented along with an illustrative case to demonstrate the applicability in LCA studies. The CFs based on estimates that reached consensus during the Delphi assessment are recommended as readily applicable and allow key differences among land use types to be distinguished. The resulting CFs are proposed as the first step for incorporating pollinator impacts in LCA studies, exemplifying the use of expert elicitation methods as a useful tool to fill data gaps that constrain the characterization of key environmental impacts.

Chronic Exposure to Polystyrene Microplastic Fragments Has No Effect on Honey Bee Survival, but Reduces Feeding Rate and Body Weight.

Microplastics (MPs), in the form of fragments and fibers, were recently found in honey samples collected in Ecuador as well as in honey bees collected from Denmark and China. However, little is known about how MPs impact bee health. To fill this knowledge gap, we investigated the potential toxicity of irregularly shaped polystyrene (PS)-MP fragments on honey bee health. In the first experiment of its kind with honey bees, we chronically exposed bees with a well-established gut microbiome to small (27 ± 17 µm) or large (93 ± 25 µm) PS-MP fragments at varying concentrations (1, 10, 100 µg mL-1) for 14 days. Bee mortality, food consumption, and body weight were all studied. We found that chronic exposure to PS-MP fragments has no effect on honey bee survival, but reduced the feeding rate and body weight, particularly at 10 µg PS-MP fragments per mL, which may have long-term consequences for honey bee health. The findings of this study could assist in the risk assessment of MPs on pollinator health.

The Insecticide Imidacloprid Decreases Nannotrigona Stingless Bee Survival and Food Consumption and Modulates the Expression of Detoxification and Immune-Related Genes.

Stingless bees are ecologically and economically important species in the tropics and subtropics, but there has been little research on the characterization of detoxification systems and immune responses within them. This is critical for understanding their responses to, and defenses against, a variety of environmental stresses, including agrochemicals. Therefore, we studied the detoxification and immune responses of a stingless bee, Nanotrigona perilampoides, which is an important stingless bee that is widely distributed throughout Mexico, including urban areas, and has the potential to be used in commercial pollination. We first determined the LC50 of the neonicotinoid insecticide imidacloprid for foragers of N. perilampoides, then chronically exposed bees for 10 days to imidacloprid at two field-realistic concentrations, LC10 (0.45 ng/µL) or LC20 (0.74 ng/µL), which are respectively 2.7 and 1.3-fold lower than the residues of imidacloprid that have been found in honey (6 ng/g) in central Mexico. We found that exposing N. perilampoides stingless bees to imidacloprid at these concentrations markedly reduced bee survival and food consumption, revealing the great sensitivity of this stingless bee to the insecticide in comparison to honey bees. The expression of detoxification (GSTD1) and immune-related genes (abaecin, defensin1, and hymenopteacin) in N. perilampoides also changed over time in response to imidacloprid. Gene expression was always lower in bees after 8 days of exposure to imidacloprid (LC10 or LC20) than it was after 4 days. Our results demonstrate that N. perilampoides stingless bees are extremely sensitive to imidacloprid, even at low concentrations, and provide greater insight into how stingless bees respond to pesticide toxicity. This is the first study of its kind to look at detoxification systems and immune responses in Mexican stingless bees, an ecologically and economically important taxon.

Intra-specific variation in sensitivity of Bombus terrestris and Osmia bicornis to three pesticides.

There is growing evidence that pesticides may be among the causes of worldwide bee declines, which has resulted in repeated calls for their increased scrutiny in regulatory assessments. One recurring concern is that the current frameworks may be biased towards assessing risks to the honey bee. This paradigm requires extrapolating toxicity information across bee species. Most research effort has therefore focused on quantifying differences in sensitivity across species. However, our understanding of how responses to pesticides may vary within a species is still very poor. Here we take the first steps towards filling this knowledge gap by comparing acute, lethal hazards in sexes and castes of the eusocial bee Bombus terrestris and in sexes of the solitary bee Osmia bicornis after oral and contact exposure to the pesticides sulfoxaflor, Amistar (azoxystrobin) and glyphosate. We show that sensitivity towards pesticides varies significantly both within and across species. Bee weight was a meaningful predictor of pesticide susceptibility. However, weight could not fully explain the observed differences, which suggests the existence of unexplored mechanisms regulating pesticide sensitivity across bee sexes and castes. Our data show that intra-specific responses are an overlooked yet important aspect of the risk assessment of pesticides in bees.

Bees under interactive stressors: the novel insecticides flupyradifurone and sulfoxaflor along with the fungicide azoxystrobin disrupt the gut microbiota of honey bees and increase opportunistic bacterial pathogens.

The gut microbiome plays an important role in bee health and disease. But it can be disrupted by pesticides and in-hive chemicals, putting honey bee health in danger. We used a controlled and fully crossed laboratory experimental design to test the effects of a 10-day period of chronic exposure to field-realistic sublethal concentrations of two nicotinic acetylcholine receptor agonist insecticides (nACHRs), namely flupyradifurone (FPF) and sulfoxaflor (Sulf), and a fungicide, azoxystrobin (Azoxy), individually and in combination, on the survival of individual honey bee workers and the composition of their gut microbiota (fungal and bacterial diversity). Metabarcoding was used to examine the gut microbiota on days 0, 5, and 10 of pesticide exposure to determine how the microbial response varies over time; to do so, the fungal ITS2 fragment and the V4 region of the bacterial 16S rRNA were targeted. We found that FPF has a negative impact on honey bee survival, but interactive (additive or synergistic) effects between either insecticide and the fungicide on honey bee survival were not statistically significant. Pesticide treatments significantly impacted the microbial community composition. The fungicide Azoxy substantially reduced the Shannon diversity of fungi after chronic exposure for 10 days. The relative abundance of the top 10 genera of the bee gut microbiota was also differentially affected by the fungicide, insecticides, and fungicide-insecticide combinations. Gut microbiota dysbiosis was associated with an increase in the relative abundance of opportunistic pathogens such as Serratia spp. (e.g. S. marcescens), which can have devastating consequences for host health such as increased susceptibility to infection and reduced lifespan. Our findings raise concerns about the long-term impact of novel nACHR insecticides, particularly FPF, on pollinator health and recommend a novel methodology for a refined risk assessment that includes the potential effects of agrochemicals on the gut microbiome of bees.

Bumble bee colony health and performance vary widely across the urban ecosystem.

Urbanization is a global phenomenon that can affect fitness and could challenge the persistence of most species, including wild bee pollinators. Yet, how and which environmental features affect bee health and fitness within the urban ecosystem remain unclear. Here, we placed experimental Bombus terrestris colonies in sites spanning from the edge into a city's core to investigate bumble bee parasitism, foraging behaviour, energetic stress, colony growth and reproductive output. In each site, ambient temperature was recorded, the availability of floral resources was evaluated and landscape heterogeneity was characterized using land-cover maps. We found that Bombus terrestris parasitism levels increased across the season in line with colony growth but were negatively related to the proportion of impervious surfaces surrounding a site. Bombus terrestris foraging trip duration decreased with increasing ecotones (edge density) but, conversely, increased in sites with honey bee hives present. Energetic stress was evaluated as lowered trehalose titre in the haemolymph of returning foragers; stress increased with the proportion of impervious surfaces. Furthermore, our analyses identified ambient temperature to be a strong predictor of Bombus terrestris colony performance in that high ambient temperature reduced colony growth and indirectly the production of sexual offspring (gynes). Our results highlight the importance of ecotones as well as minimizing the intensity of urbanization and urban honey bee beekeeping for bumble bee colony health and foraging behaviour. They also point to the importance of microclimate (i.e. temperature) for bumble bee colony performance and suggest that increasing temperatures could have a negative impact in slowing colony weight gain, and indirectly in reducing colony reproduction.