Honey Bee Viruses

Nutrition, pesticide exposure, and virus infection interact to produce context-dependent effects in honey bees (Apis mellifera)

Declines in pollinator health are frequently hypothesized to be the combined result of multiple interacting biotic and abiotic stressors; namely, nutritional limitations, pesticide exposure, and infection with pathogens and parasites. Despite this hypothesis, most studies examining stressor interactions have been constrained to two concurrent factors, limiting our understanding of multi-stressor dynamics. Using honey bees as a model, we addressed this gap by studying how variable diet, field-realistic levels of multiple pesticides, and virus infection interact to affect survival, infection intensity, and immune and detoxification gene expression. Although we found evidence that agrochemical exposure (a field-derived mixture of chlorpyrifos and two fungicides) can exacerbate infection and increase virus-induced mortality, this result was nutritionally-dependent, only occurring when bees were provided artificial pollen. Provisioning with naturally-collected polyfloral pollen inverted the effect, reducing virus-induced mortality and suggesting a hormetic response. To test if the response was pesticide specific, we repeated our experiment with a pyrethroid (lambda-cyhalothrin) and a neonicotinoid (thiamethoxam), finding variable results. Finally, to understand the underpinnings of these effects, we measured viral load and expression of important immune and detoxification genes. Together, our results show that multi-stressor interactions are complex and highly context-dependent, but have great potential to affect bee health and physiology.

Molecular Detection and Phylogenetic Relationships of Honey Bee-Associated Viruses in Bee Products

In the last few years, the isolation and amplification of DNA or RNA from the environment (eDNA/eRNA) has proven to be an alternative and non-invasive approach for molecular identification of pathogens and pests in beekeeping. We have recently demonstrated that bee pollen and bee bread represent suitable biological material for the molecular identification of viral RNA. In the present study, we extracted total RNA from different bee products (pollen, n = 25; bee bread, n = 17; and royal jelly, n = 15). All the samples were tested for the presence of six of the most common honey bee-associated viruses-Deformed wing virus (DWV), Acute bee paralysis virus (ABPV), Chronic bee paralysis virus (CBPV), Sacbrood virus (SBV), Kashmir bee virus (KBV), and Black queen cell virus (BQCV)-using a reverse transcription polymerase chain reaction (RT-PCR). We successfully detected six records of DWV (10.5%, 6/57), four of ABPV (7.0%, 4/57), three of Israeli acute paralysis virus (IAPV) (5.3%, 3/57), and two of BQCV (3.5%, 2/57). Using ABPV primers, we also successfully detected the presence of IAPV. The obtained viral sequences were analyzed for phylogenetic relationships with the highly similar sequences (megablast) available in the GenBank database. The Bulgarian DWV isolates revealed a high homology level with strains from Syria and Turkey. Moreover, we successfully detected a DWV strain B for the first time in Bulgaria. In contrast to DWV, the ABPV isolates formed a separate clade in the phylogenetic tree. BQCV was closely grouped with Russian isolates, while Bulgarian IAPV formed its own clade and included a strain from China. In conclusion, the present study demonstrated that eRNA can be successfully used for molecular detection of honey bee-associated viruses in bee products. The method can assist the monitoring of the health status of honey bee colonies at the local, regional, and even national levels.

Do Varroa destructor (Acari: Varroidae) mite flows between Apis mellifera (Hymenoptera: Apidae) colonies bias colony infestation evaluation for resistance selection?

Since the global invasion of the ectoparasitic mite Varroa destructor (Anderson and Trueman), selection of mite-resistant honey bee (Apis mellifera L.) colonies appears challenging and has to date not broadly reduced colony mortality. The low published estimated heritability values for mite infestation levels could explain the limited genetic progresses obtained so far. We hypothesize that intercolonial horizontal mite transmission could differentially affect the single colonies located in a given apiary and therefore invisibly bias colony infestation phenotypes. This bias may be lower in regions with lower colony density, providing suitable conditions to set up evaluation apiaries. To verify these hypotheses, we monitored mite infestation and reinvasion in experimental colonies, as well as infestation in neighboring colonies belonging to beekeepers in three areas with variable colony densities in the canton of Bern, Switzerland during three consecutive beekeeping seasons. Mite immigration fluctuated between apiaries and years and significantly contributed to colony infestation level. Depending on apiary and year, 17-48% of the mites present in the experimental colonies at the time of the summer oxalic acid final treatment potentially derived from mite immigration that had occurred since mid-spring. Mite immigration was not linked to local colony density or the infestation levels of beekeepers' colonies located within 2 km. Our results do not prove that apiaries for colony evaluation should necessarily be established in areas with low colony density. However, they highlight the high impact of beekeeping management practices on mite colony infestation levels.

Oxalic Acid Treatment: Short-Term Effects on Enzyme Activities, Vitellogenin Content, and Residual Oxalic Acid Content in House Bees, Apis mellifera L

Honeybees (Apis mellifera L.) have to face many challenges, including Varroa destructor infestation, associated with viral transmission. Oxalic acid is one of the most common treatments against Varroa. Little is known about the physiological effects of oxalic acid, especially those on honeybees' immune systems. In this study, the short-term effects (0-96 h) of oxalic acid treatment on the immune system components (i.e., glucose oxidase, phenoloxidase, glutathione S-transferase, catalase activities, and vitellogenin contents) of house bees were preliminarily investigated. Oxalic acid contents of bee bodies and haemolymphs were also measured. The results confirm that oxalic acid is constitutively present in bee haemolymphs and its concentration is not affected by treatment. At 6 h after the treatment, a maximum peak of oxalic acid content was detected on bees' bodies, which gradually decreased after that until physiological levels were reached at 48 h. In the immune system, the oxalic acid treatment determined a peak in glucose oxidase activity at 48 h, indicating a potential defence response and an increase in vitellogenin content at 24 h. No significant changes were recorded in phenoloxidase, glutathione S-transferase, and catalase activities. These results suggest a time-dependent response to oxalic acid, with potential immune system activation in treated bees.

Eggs sampling as an effective tool for identifying the incidence of viruses in honey bees involved in artificial queen rearing

The Carniolan honey bee (Apis mellifera carnica) plays an essential role in crop pollination, environment diversity, and the production of honey bee products. However, the health of individual honey bees and their colonies is under pressure due to multiple stressors, including viruses as a significant threat to bees. Monitoring various virus infections could be a crucial selection tool during queen rearing. In the present study, samples from all developmental stages (eggs, larvae, pupae, and queens) were screened for the incidence of seven viruses during queen rearing in Slovenia. The screening of a total of 108 samples from five queen breeders was performed by the RT-qPCR assays. The results showed that the highest incidence was observed for black queen cell virus (BQCV), Lake Sinai virus 3 (LSV3), deformed wing virus B (DWV-B), and sacbrood virus (SBV). The highest viral load was detected in queens (6.07 log10 copies/queen) and larvae (5.50 log10 copies/larva) for BQCV, followed by SBV in larvae (5.47 log10 copies/larva). When comparing all the honey bee developmental stages, the eggs exhibited general screening for virus incidence and load in queen mother colonies. The results suggest that analyzing eggs is a good indicator of resilience to virus infection during queen development.

Apis mellifera filamentous virus from a honey bee gut microbiome survey in Hungary

In Hungary, as part of a nationwide, climatically balanced survey for a next-generation sequencing-based study of the honey bee (Apis mellifera) gut microbiome, repeated sampling was carried out during the honey production season (March and May 2019). Among other findings, the presence of Apis mellifera filamentous virus (AmFV) was detected in all samples, some at very high levels. AmFV-derived reads were more abundant in the March samples than in the May samples. In March, a higher abundance of AmFV-originated reads was identified in samples collected from warmer areas compared to those collected from cooler areas. A lower proportion of AmFV-derived reads were identified in samples collected in March from the wetter areas than those collected from the drier areas. AmFV-read abundance in samples collected in May showed no significant differences between groups based on either environmental temperature or precipitation. The AmFV abundance correlated negatively with Bartonella apihabitans, Bartonella choladocola, and positively with Frischella perrara, Gilliamella apicola, Gilliamella sp. ESL0443, Lactobacillus apis, Lactobacillus kullabergensis, Lactobacillus sp. IBH004. De novo metagenome assembly of four samples resulted in almost the complete AmFV genome. According to phylogenetic analysis based on DNA polymerase, the Hungarian strains are closest to the strain CH-05 isolated in Switzerland.

Seasonal trends of the ABPV, KBV, and IAPV complex in Italian managed honey bee (Apis mellifera L.) colonies

Acute bee paralysis virus (ABPV), Kashmir bee virus (KBV), and Israeli acute paralysis virus (IAPV) usually persist as covert infections in honey bee colonies. They can cause rapid bee mortality in cases of severe infection, often associated with high Varroa destructor infestation, by which they are transmitted. In various countries, these viruses have been associated with colony collapse. Despite their potential danger, these viruses are often disregarded, and little information is available on their occurrence in many countries, including Italy. In 2021, 370 apiaries representing all of the Italian regions were investigated in four different months (June, September, November, and March) for the presence of ABPV, KBV, and IAPV. IAPV was not found in any of the apiaries investigated, whereas 16.45% and 0.67% of the samples tested positive for ABPV and KBV, respectively. Most ABPV cases occurred in late summer-autumn in both northern and southern regions. We observed a scattered pattern of KBV-positive colonies that did not allow any seasonal or regional trends to be discerned. Differences observed among regions and months were potentially related to the dynamics of varroa infestation, viral genetic variations, and different climatic conditions resulting in variations in bee behaviour. This study improves our understanding of the circulation of bee viruses and will contribute to better disease prevention and preservation of bee health.

Plants affect the horizontal transmission of a new densovirus infecting the green peach aphid Myzus persicae by modulating honeydew production

In a tritrophic context of plant-insect-entomopathogen, plants play important roles in modulating the interaction of insects and their pathogenic viruses. Currently, the influence of plants on the transmission of insect viruses has been mainly studied on baculoviruses and some RNA viruses, whereas the impact of plants on other insect viruses is largely unknown. Here, we identified a new densovirus infecting the green peach aphid Myzus persicae and tested whether and how host plants influence the transmission of the aphid densovirus. The complete single-stranded DNA genome of the virus, M. persicae densovirus 2, is 5 727 nt and contains inverted terminal repeats. Transcription and phylogenetic analysis indicated that the virus was distinct from other a few identified aphid densoviruses. The virus abundance was detected highly in the intestinal tract of aphids, compared with the lower level of it in other tissues including head, embryo, and epidermis. Cabbage and pepper plants had no obvious effect on the vertical transmission and saliva-mediated horizontal transmission of the virus. However, the honeydew-mediated horizontal transmission among aphids highly depended on host plants (65% on cabbages versus 17% on peppers). Although the virus concentration in the honeydew produced by aphids between 2 plants was similar, the honeydew production of the infected aphids reared on peppers was dramatically reduced. Taken together, our results provide evidence that plants influence the horizontal transmission of a new densovirus in an aphid population by modulating honeydew secretion of aphids, suggesting plants may manipulate the spread of an aphid-pathogenic densovirus in nature.

Evaluation of the Potential Effect of Postbiotics Obtained from Honey Bees against Varroa destructor and Their Combination with Other Organic Products

The Varroa destructor mite infests Apis mellifera colonies and causes significant harm. Traditional treatments have become less effective because of mite resistance development and can also generate residues inside beehives. This study aimed to gauge the efficacy of a beehive-derived postbiotic in reducing V. destructor viability and to explore its synergies with organic compounds. Four lactic acid bacteria (LAB) species, Leuconostoc mesenteroides, Lactobacillus helsingborgensis, Bacillus velezensis, and Apilactobacillus kunkeei, were isolated and tested in a postbiotic form (preparations of inanimate microorganisms and/or their components) via bioassays. L. mesenteroides, L. helsingborgensis, and B. velezensis notably reduced the mite viability compared to the control, and they were further tested together as a single postbiotic product (POS). Further bioassays were performed to assess the impact of the POS and its combinations with oxalic acid and oregano essential oil. The simple products and combinations (POS/Oregano, POS/Oxalic, Oregano/Oxalic, and POS/Oregano/Oxalic) decreased the mite viability. The most effective were the oxalic acid combinations (POS/Oregano/Oxalic, Oxalic/Oregano, POS/Oxalic), showing significant improvements compared to the individual products. These findings highlight the potential of combining organic products as a vital strategy for controlling V. destructor infection. This study suggests that these combinations could serve as essential tools for combating the impact of mites on bee colonies.

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.

Establishment and application of bioassay- and molecular marker-based methods for monitoring fluvalinate resistance of Varroa mites

The Varroa mite, Varroa destructor, is an ectoparasite that infests honey bees. The extensive use of acaricides, including fluvalinate, has led to the emergence of resistance in Varroa mite populations worldwide. This study's objective is to monitor fluvalinate resistance in field populations of Varroa mites in Korea through both bioassay-based and molecular marker-based methods. To achieve this, a residual contact vial (RCV) bioassay was established for on-site resistance monitoring. A diagnostic dose of 200 ppm was determined based on the bioassay using a putative susceptible population. In the RCV bioassay, early mortality evaluation was effective for accurately discriminating mites with the knockdown resistance (kdr) genotype, while late evaluation was useful for distinguishing mites with additional resistance factors. The RCV bioassay of 14 field mite populations collected in 2021 indicated potential resistance development in four populations. As an alternative approach, quantitative sequencing was employed to assess the frequency of the L925I/M mutation in the voltage-gated sodium channel (VGSC), associated with fluvalinate kdr trait. While the mutation was absent in 2020 Varroa mite populations, it emerged in 2021, increased in frequency in 2022, and became nearly widespread across the country by 2023. This recent emergence and rapid spread of fluvalinate resistance within a span of three years demonstrate the Varroa mite's significant potential for developing resistance. This situation further underscores the urgent need to replace fluvalinate with alternative acaricides. A few novel VGSC mutations potentially involved in resistance were identified. Potential factors driving the rapid expansion of resistance were further discussed.

The threat of honey bee RNA viruses to yellow-legged hornets: Insights from cross-species transmission events

Viral diseases are a significant challenge in beekeeping, and recent studies have unveiled a potential link between these diseases and the yellow-legged hornets (Vespa velutina), notorious predators of honey bees. However, it remains unclear whether virus diseases are commonly shared between honey bees and hornets or are merely sporadic cross-species transmission events. To address this knowledge gap, we conducted a study utilizing hornet-keeping practices in Yunnan, Southwest China. Our findings demonstrate that deformed wing virus (DWV-A) and Israeli acute paralysis virus (IAPV) can be transmitted from honey bees to yellow-legged hornets. We detected virus replication in various hornet stages, including pupae with IAPV infections, indicating the similarities between infected hornet and honey bee stages. Furthermore, we observed signs and infection intensities of DWV-A and IAPV comparable to those in honey bees. While different polymorphisms were found in the virus isolates from yellow-legged hornets, the sequences remain similar to honey bee counterparts. While our findings suggest that DWV-A and IAPV behave like common diseases, we observed a natural elimination of the viruses in hornet colonies, with minimal alterations in viral sequences. Consequently, these events appear to be cross-species transmission from honey bees, with yellow-legged hornets acting as potential incidental hosts. Further investigations of virus monitoring in hornets promise valuable insights into the disease ecology of bee-infecting viruses.

Mediating a host cell signaling pathway linked to overwinter mortality offers a promising therapeutic approach for improving bee health

INTRODUCTION

Honey bees provides valuable pollination services for world food crops and wild flowering plants which are habitats of many animal species and remove carbon dioxide from the atmosphere, a powerful tool in the fight against climate change. Nevertheless, the honey bee population has been declining and the majority of colony losses occur during the winter.

OBJECTIVES

The goal of this study was to understand the mechanisms underlying overwinter colony losses and develop novel therapeutic strategies for improving bee health.

METHODS

First, pathogen prevalence in overwintering bees were screened between 2015 and 2018. Second, RNA sequencing (RNA-Seq) for transcriptional profiling of overwintering honey bees was conducted and qRT-PCR was performed to confirm the results of the differential expression of selected genes. Lastly, laboratory bioassays were conducted to measure the effects of cold challenges on bee survivorship and stress responses and to assess the effect of a novel medication for alleviating cold stress in honey bees.

RESULTS

We identified that sirtuin signaling pathway is the most significantly enriched pathway among the down-regulated differentially expressed genes (DEGs) in overwintering diseased bees. Moreover, we showed that the expression of SIRT1 gene, a major sirtuin that regulates energy and immune metabolism, was significantly downregulated in bees merely exposed to cold challenges, linking cold stress with altered gene expression of SIRT1. Furthermore, we demonstrated that activation of SIRT1 gene expression by SRT1720, an activator of SIRT1 expression, could improve the physiology and extend the lifespan of cold-stressed bees.

CONCLUSION

Our study suggests that increased energy consumption of overwintering bees for maintaining hive temperature reduces the allocation of energy toward immune functions, thus making the overwintering bees more susceptible to disease infections and leading to high winter colony losses. The novel information gained from this study provides a promising avenue for the development of therapeutic strategies for mitigating colony losses, both overwinter and annually.

Signatures of adaptive decreased virulence of deformed wing virus in an isolated population of wild honeybees (Apis mellifera)

Understanding the ecological and evolutionary processes that drive host-pathogen interactions is critical for combating epidemics and conserving species. The Varroa destructor mite and deformed wing virus (DWV) are two synergistic threats to Western honeybee (Apis mellifera) populations across the globe. Distinct honeybee populations have been found to self-sustain despite Varroa infestations, including colonies within the Arnot Forest outside Ithaca, NY, USA. We hypothesized that in these bee populations, DWV has been selected to produce an avirulent infection phenotype, allowing for the persistence of both host and disease-causing agents. To investigate this, we assessed the titre of viruses in bees from the Arnot Forest and managed apiaries, and assessed genomic variation and virulence differences between DWV isolates. Across groups, we found viral abundance was similar, but DWV genotypes were distinct. We also found that infections with isolates from the Arnot Forest resulted in higher survival and lower rates of symptomatic deformed wings, compared to analogous isolates from managed colonies, providing preliminary evidence to support the hypothesis of adaptive decreased viral virulence. Overall, this multi-level investigation of virus genotype and phenotype indicates that host ecological context can be a significant driver of viral evolution and host-pathogen interactions in honeybees.

Can the examination of different types of hive samples be a non-invasive method for detection and quantification of viruses in honey bee (Apis mellifera L.) colonies?

Introduction

Honey bee viruses have been shown to negatively affect the vigour and longevity of European honey bees (Apis mellifera L). In the present work, beehive materials were tested for their potential to serve as non-invasive samples for honey bee virus detection.

Material and Methods

Honey, pollen, hive debris, hive grid smears and forager honey bees were collected from 24 hives at four locations in the Czech Republic. Deformed wing virus (DWV), acute bee paralysis virus (ABPV), sacbrood virus (SBV) and black queen cell virus (BQCV) were detected using a reverse transcription PCR (RT-PCR) and real-time quantitative RT-PCR and the results for bees and alternative materials compared.

Results

All forager bee samples contained DWV, BQCV and SBV and 54.2% had ABPV. When comparing beehive materials to bees, the most promising results were obtained from honey and pollen samples, with BQCV and SBV detected in all honey samples and ABPV in 12.5%. Detection of SBV was achieved in 91.6% of pollen samples, detection of BQCV in 87.5% and detection of DWW in 75%. The results for debris and smears were less consistent with the viral profile of the forager samples.

Conclusion

The best candidate materials for honey bee virus detection in a non-invasive technique are honey and pollen.

Impacts of seasonality and parasitism on honey bee population dynamics

The honeybee plays an extremely important role in ecosystem stability and diversity and in the production of bee pollinated crops. Honey bees and other pollinators are under threat from the combined effects of nutritional stress, parasitism, pesticides, and climate change that impact the timing, duration, and variability of seasonal events. To understand how parasitism and seasonality influence honey bee colonies separately and interactively, we developed a non-autonomous nonlinear honeybee-parasite interaction differential equation model that incorporates seasonality into the egg-laying rate of the queen. Our theoretical results show that parasitism negatively impacts the honey bee population either by decreasing colony size or destabilizing population dynamics through supercritical or subcritical Hopf-bifurcations depending on conditions. Our bifurcation analysis and simulations suggest that seasonality alone may have positive or negative impacts on the survival of honey bee colonies. More specifically, our study indicates that (1) the timing of the maximum egg-laying rate seems to determine when seasonality has positive or negative impacts; and (2) when the period of seasonality is large it can lead to the colony collapsing. Our study further suggests that the synergistic influences of parasitism and seasonality can lead to complicated dynamics that may positively and negatively impact the honey bee colony's survival. Our work partially uncovers the intrinsic effects of climate change and parasites, which potentially provide essential insights into how best to maintain or improve a honey bee colony's health.

Endoparasitoid lifestyle promotes endogenization and domestication of dsDNA viruses

The accidental endogenization of viral elements within eukaryotic genomes can occasionally provide significant evolutionary benefits, giving rise to their long-term retention, that is, to viral domestication. For instance, in some endoparasitoid wasps (whose immature stages develop inside their hosts), the membrane-fusion property of double-stranded DNA viruses have been repeatedly domesticated following ancestral endogenizations. The endogenized genes provide female wasps with a delivery tool to inject virulence factors that are essential to the developmental success of their offspring. Because all known cases of viral domestication involve endoparasitic wasps, we hypothesized that this lifestyle, relying on a close interaction between individuals, may have promoted the endogenization and domestication of viruses. By analyzing the composition of 124 Hymenoptera genomes, spread over the diversity of this clade and including free-living, ecto, and endoparasitoid species, we tested this hypothesis. Our analysis first revealed that double-stranded DNA viruses, in comparison with other viral genomic structures (ssDNA, dsRNA, ssRNA), are more often endogenized and domesticated (that is, retained by selection) than expected from their estimated abundance in insect viral communities. Second, our analysis indicates that the rate at which dsDNA viruses are endogenized is higher in endoparasitoids than in ectoparasitoids or free-living hymenopterans, which also translates into more frequent events of domestication. Hence, these results are consistent with the hypothesis that the endoparasitoid lifestyle has facilitated the endogenization of dsDNA viruses, in turn, increasing the opportunities of domestications that now play a central role in the biology of many endoparasitoid lineages.

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.

First molecular detection of the presence of honey bee viruses in insects, Varroa destructor mites, and pollinated plants in an isolated region of Armenia

Background and Aim

Recently, viral diseases of honey bees (Apis mellifera) have presented an increasing threat to beekeeping. This study aimed to examine the presence of honey bee viruses in Apis and non-Apis bee species, the mite Varroa destructor, and pollinated plants in Armenia.

Materials and Methods

Sampling was performed in Tavush Province, in the northeast of the Republic of Armenia, from August to November 2019. Overall, 200 A. mellifera bees, 50 V. destructor mites, and 20 wasps were collected (corresponding to three bees, five mites, and 2-11 wasps in each investigated sample) and homogenized for RNA isolation and detection of viruses. Ten pollinated plants were taken from each plant, and 2 g of each sample was used for homogenization. In each investigated case Apis mellifera, Varroa destructor, Vespula germanica and plants received percentages of the virus presence.

Results

Six important honey bee viruses (acute bee paralysis virus [ABPV], deformed wing virus [DWV], A. mellifera norovirus [ANV], Lake Sinai virus-2 [LSV-2], Big Sioux River virus [BSRV], and A. mellifera filamentous virus [AmFV]) were detected in samples by polymerase chain reaction. Our results showed that DWV, ANV, and ABPV were the most common viruses in honey bees. All viruses were detected in wasps, but LSV-2 and ANV were present in almost all samples.

Conclusion

Our results showed that almost all viruses were present in V. destructor. Although ANV is very common in honey bees, it did not appear in any mite samples. Our study indicates that viruses typically associated with honey bees were also actively infecting wasps. Our data suggest that the survival of viruses in plants can be an important source of seasonal transmission of viruses to bees. In addition, pollinated plants can potentially serve as reservoirs for honey bee viruses.

Nationwide genomic surveillance reveals the prevalence and evolution of honeybee viruses in China

BACKGROUND

The economic and environmental value of honeybees has been severely challenged in recent years by the collapse of their colonies worldwide, often caused by outbreaks of infectious diseases. However, our understanding of the diversity, prevalence, and transmission of honeybee viruses is largely obscure due to a lack of large-scale and longitudinal genomic surveillance on a global scale.

RESULTS

We report the meta-transcriptomic sequencing of nearly 2000 samples of the two most important economic and widely maintained honeybee species, as well as an associated ectoparasite mite, collected across China during 2016-2019. We document the natural diversity and evolution of honeybee viruses in China, providing evidence that multiple viruses commonly co-circulate within individual bee colonies. We also expanded the genomic data for 12 important honeybee viruses and revealed novel genetic variants and lineages associated with China. We identified more than 23 novel viruses from the honeybee and mite viromes, with some exhibiting ongoing replication in their respective hosts. Together, these data provide additional support to the idea that mites are an important reservoir and spill-over host for honeybee viruses.

CONCLUSIONS

Our data show that honeybee viruses are more widespread, prevalent, and genetically diverse than previously realized. The information provided is important in mitigating viral infectious diseases in honeybees, in turn helping to maintain sustainable productive agriculture on a global scale. Video Abstract.

Metagenomic analysis of viromes in honey bee colonies (Apis mellifera; Hymenoptera: Apidae) after mass disappearance in Korea

After the nationwide, massive winter losses of honey bees in Korea during the winter of 2021, samplings were conducted from live honey bees in colonies and dead honey bees nearby colonies in the same bee-farms in six regions in Korea. Each sample was subjected to virome analysis using high-throughput sequencing technology. The number of viral reads was the lowest in the live honey bee group sample with 370,503 reads and the highest in the dead honey bee group sample with 42,659,622 reads. Viral contigs were matched with the viral genomes of the black queen cell virus, deformed wing virus, Israeli acute paralysis virus, and sacbrood virus, all of which have been previously reported in Korea. However, Apis rhabdovirus 5, bee macula-like virus, Varroa orthomyxovirus-1, Hubei partiti-like virus 34, Lake Sinai virus 2, 3, and 4, and the Ditton virus, were also discovered in this study, which are the first records in Korea. Plant viral sequences resembling those of Arabidopsis latent virus 1, and a novel viral sequence was also discovered. In the present study 55 complete viral genome sequences were identified. This study is the first virome analysis of domestic honey bees and provides the latest information on the diversity of honey bee viruses in Korea.

First Detection of Honeybee Pathogenic Viruses in Butterflies

Several pathogens are important causes of the observed pollinator decline, some of which could be transmitted between different pollinator species. To determine whether honeybee viruses can be transmitted to butterflies, a total of 120 butterflies were sampled at four locations in Slovenia. At each location, butterflies from three families (Pieridae, Nymphalidae, Hesperiidae/Lycenidae) and Carniolan honeybees (Apis mellifera carnica) were collected. The RNA of six honeybee viruses, i.e., acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), deformed wing virus A (DWV-A), Sacbrood bee virus (SBV), and Lake Sinai virus 3 (LSV3), was detected by a specific quantitative method (RT-PCR). The presence of ABPV, BQCV, LSV3, and SBV was detected in both butterflies and honeybees. All butterfly and bee samples were negative for CBPV, while DWV-A was detected only in honeybees. The viral load in the positive butterfly samples was much lower than in the positive bee samples, which could indicate that butterflies are passive carriers of bee viruses. The percentage of positive butterfly samples was higher when the butterflies were collected at sampling sites with a higher density of apiaries. Therefore, we believe that infected bees are a necessary condition for the presence of viruses in cohabiting butterflies. This is the first study on the presence of pathogenic bee viruses in butterflies.

Transmission of deformed wing virus between Varroa destructor foundresses, mite offspring and infested honey bees

Background

Varroa destructor is the major ectoparasite of the western honey bee (Apis mellifera). Through both its parasitic life-cycle and its role as a vector of viral pathogens, it can cause major damage to honey bee colonies. The deformed wing virus (DWV) is the most common virus transmitted by this ectoparasite, and the mite is correlated to increased viral prevalence and viral loads in infested colonies. DWV variants A and B (DWV-A and DWV-B, respectively) are the two major DWV variants, and they differ both in their virulence and transmission dynamics.

Methods

We studied the transmission of DWV between bees, parasitic mites and their offspring by quantifying DWV loads in bees and mites collected in in vitro and in situ environments. In vitro, we artificially transmitted DWV-A to mites and quantified both DWV-A and DWV-B in mites and bees. In situ, we measured the natural presence of DWV-B in bees, mites and mites’ offspring.

Results

Bee and mite viral loads were correlated, and mites carrying both variants were associated with higher mortality of the infected host. Mite infestation increased the DWV-B loads and decreased the DWV-A loads in our laboratory conditions. In situ, viral quantification in the mite offspring showed that, after an initially non-infected egg stage, the DWV-B loads were more closely correlated with the foundress (mother) mites than with the bee hosts.

Conclusions

The association between mites and DWV-B was highlighted in this study. The parasitic history of a mite directly impacts its DWV infection potential during the rest of its life-cycle (in terms of variant and viral loads). Regarding the mite’s progeny, we hypothesize that the route of contamination is likely through the feeding site rather than by vertical transmission, although further studies are needed to confirm this hypothesis.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05463-9.

Sacbrood Virus: A Growing Threat to Honeybees and Wild Pollinators

Sacbrood virus (SBV) is one of the many viruses that infect both the Western honeybee (Apis mellifera) and the Eastern honeybee (Apis cerana). Recently, the interspecies transmission of SBV has been discovered, especially among wild pollinators. This newly discovered evolutionary occurrence regarding SBV indicates a much wider host range than previously believed, causing further concern about the future sustainability of agriculture and the resilience of ecosystems. Over the past few decades, vast numbers of studies have been undertaken concerning SBV infection in honeybees, and remarkable progress has been made in our understanding of the epidemiology, pathogenesis, transmission, and manifestations of SBV infection in honeybees and other pollinators. Meanwhile, some methods, including Chinese medicine, have been established to control and prevent sacbrood disease in A. cerana in Asian countries. In this review, we summarize the existing knowledge of SBV and address the gaps in the knowledge within the existing literature in the hope of providing future directions for the research and development of management strategies for controlling the spread of this deadly disease.

The Epidemiological Situation of the Managed Honey Bee (Apis mellifera) Colonies in the Italian Region Emilia-Romagna

The recent decades witnessed the collapse of honey bee colonies at a global level. The major drivers of this collapse include both individual and synergic pathogen actions, threatening the colonies’ survival. The need to define the epidemiological pattern of the pathogens that are involved has led to the establishment of monitoring programs in many countries, Italy included. In this framework, the health status of managed honey bees in the Emilia–Romagna region (northern Italy) was assessed, throughout the year 2021, on workers from 31 apiaries to investigate the presence of major known and emerging honey bee pathogens. The prevalence and abundance of DWV, KBV, ABPV, CBPV, Nosema ceranae, and trypanosomatids (Lotmaria passim, Crithidia mellificae, Crithidia bombi) were assessed by molecular methods. The most prevalent pathogen was DWV, followed by CBPV and N. ceranae. Trypanosomatids were not found in any of the samples. Pathogens had different peaks in abundance over the months, showing seasonal trends that were related to the dynamics of both bee colonies and Varroa destructor infestation. For some of the pathogens, a weak but significant correlation was observed between abundance and geographical longitude. The information obtained in this study increases our understanding of the epidemiological situation of bee colonies in Emilia–Romagna and helps us to implement better disease prevention and improved territorial management of honey bee health.

Colony-Level Viral Load Influences Collective Foraging in Honey Bees

Nutrition is an important component of social insect colony health especially in the face of stressors such as parasitism and viral infections. Honey bees are known to preferentially select nectar and pollen based on macronutrient and phytochemical contents and in response to pathogen loads. However, given that honey bees live in colonies, collective foraging decisions may be impacted directly by forager infection status but also by colony health. This field experiment was conducted to determine if honey bee viral infections are correlated with pollen and nectar foraging and if these associations are impacted more by colony or forager infection. By comparing regressions with and without forager and colony variables and through structural equation models, we were able to determine the relative contributions of colony and forager virus loads on forager decisions. We found that foragers had higher numbers and levels of BQCV and CBPV but lower levels of DWV viruses than their respective colonies. Overall, individuals appeared to forage based a combination of their own and colony health but with greater weight given to colony metrics. Colony parasitism by Varroa mites, positively correlated with both forager and colony DWV-B levels, was negatively associated with nectar weight. Further, colony DWV-B levels were negatively associated with individually foraged pollen protein: lipid ratios but positively correlated with nectar weight and sugar content. This study shows that both colony and forager health can simultaneously mediate individual foraging decisions and that the importance of viral infections and parasite levels varies with foraging metrics. Overall, this work highlights the continued need to explore the interactions of disease, nutrition, and genetics in social interactions and structures.

Molecular Detection and Differentiation of Arthropod, Fungal, Protozoan, Bacterial and Viral Pathogens of Honeybees

The honeybee Apis mellifera is highly appreciated worldwide because of its products, but also as it is a pollinator of crops and wild plants. The beehive is vulnerable to infections due to arthropods, fungi, protozoa, bacteria and/or viruses that manage to by-pass the individual and social immune mechanisms of bees. Due to the close proximity of bees in the beehive and their foraging habits, infections easily spread within and between beehives. Moreover, international trade of bees has caused the global spread of infections, several of which result in significant losses for apiculture. Only in a few cases can infections be diagnosed with the naked eye, by direct observation of the pathogen in the case of some arthropods, or by pathogen-associated distinctive traits. Development of molecular methods based on the amplification and analysis of one or more genes or genomic segments has brought significant progress to the study of bee pathogens, allowing for: (i) the precise and sensitive identification of the infectious agent; (ii) the analysis of co-infections; (iii) the description of novel species; (iv) associations between geno- and pheno-types and (v) population structure studies. Sequencing of bee pathogen genomes has allowed for the identification of new molecular targets and the development of specific genotypification strategies.

Presence of Known and Emerging Honey Bee Pathogens in Apiaries of Veneto Region (Northeast of Italy) during Spring 2020 and 2021

A progressive honey bee population decline has been reported worldwide during the last decades, and it could be attributed to several causes, in particular to the presence of pathogens and parasites that can act individually or in synergy. The health status of nine apiaries located in different areas of the Veneto region (northeast of Italy) was assessed for two consecutive years (2020 and 2021) in spring, during the resumption of honey bee activity, for determining the presence of known (Nosema spp., Varroa mite and viruses) and less known or emerging pathogens (Lotmaria passim and Crithidia mellificae) in honey bees. After honey bees sampling from each of the nine apiaries, Nosema apis, Nosema ceranae, L. passim, C. mellificae, ABPV, CBPV, IAPV, KBV, BQCV, SBV, DWV-A, DWV-B and V. destructor were investigated either by microscopic observation or PCR protocols. The viruses BQCV, SBV, CBPV followed by N. ceranae and L. passim were the most prevalent pathogens, and many of the investigated hives, despite asymptomatic, had different degrees of co-infection. This study aimed to highlight, during the resumption of honey bee activity in spring, the prevalence and spreading in the regional territory of different honey bee pathogens, which could alone or synergistically alter the homeostasis of bees colonies. The information gathered would increase our knowledge about the presence of these microorganisms and parasites in the territory and could contribute to improve beekeepers practice.

Virome Analysis Reveals Diverse and Divergent RNA Viruses in Wild Insect Pollinators in Beijing, China

Insect pollinators provide major pollination services for wild plants and crops. Honeybee viruses can cause serious damage to honeybee colonies. However, viruses of other wild pollinating insects have yet to be fully explored. In the present study, we used RNA sequencing to investigate the viral diversity of 50 species of wild pollinating insects. A total of 3 pathogenic honeybee viruses, 8 previously reported viruses, and 26 novel viruses were identified in sequenced samples. Among these, 7 novel viruses were shown to be closely related to honeybee pathogenic viruses, and 4 were determined to have potential pathogenicity for their hosts. The viruses detected in wild insect pollinators were mainly from the order Picornavirales and the families Orthomyxoviridae, Sinhaliviridae, Rhabdoviridae, and Flaviviridae. Our study expanded the species range of known insect pollinator viruses, contributing to future efforts to protect economic honeybees and wild pollinating insects.

Research progress on VARROA DESTRUCTOR

Varroa destructor is the most serious threat for western honeybees worldwide. It parasitizes mostly on honeybee's pupae and larvae, resulting in impaired development, decreased immunity, residual wings, loss of flight ability, and even direct death. At the same time, it can cause a series of diseases, greatly infesting bee colonies. This paper introduces the classification, distribution and biological characteristics of Varroa destructor, and discusses its transmission mode, harm form and common control measures.

Warmer winters are associated with lower levels of the cryoprotectant glycerol, a slower decrease in vitellogenin expression and reduced virus infections in winter honeybees

Within the context of climate change, winter temperatures at high latitudes are predicted to rise faster than summer temperatures. This phenomenon is expected to negatively affect the diapause performance and survival of insects, since they largely rely on low temperatures to lower their metabolism and preserve energy. However, some insects like honeybees, remain relatively active during the winter and elevate their metabolic rate to produce endothermic heat when temperatures drop. Warming winters are thus expected to improve overwintering performance of honeybees. In order to verify this hypothesis, for two consecutive years, we exposed honeybee colonies to either a mild or cold winter. We then monitored the influence of wintering conditions on several parameters of honeybee overwintering physiology, such as levels of the cryoprotectant glycerol, expression levels of immune and antioxidant genes, and genes encoding multifunctional proteins, including vitellogenin, which promotes bee longevity. Winter conditions had no effect on the expression of antioxidant genes, and genes related to immunity were not consistently affected. However, mild winters were consistently associated with a lower investment in glycerol synthesis and a higher expression of fat body genes, especially apidaecin and vitellogenin. Finally, while we found that viral loads generally decreased through the winter, this trend was more pronounced under mild winter conditions. In conclusion, and without considering how warming temperatures might affect other aspects of honeybee biology before overwintering, our data suggest that warming temperatures will likely benefit honeybee vitality by notably reducing their viral loads over the winter.

Factors affecting virus prevalence in honey bees in the Pacific-Northwest, USA

Global efforts to assess honey bee health show viruses are major stressors that undermine colony performance. Identifying factors that affect virus incidence, such as management practices and landscape context, could aid in slowing virus transmission. Here we surveyed viruses in honey bees from 86 sites in the Pacific Northwest, USA, and tested effects of regional bee density, movement associated with commercial pollination, julian date, and hive management on virus prevalence. We also explored patterns of virus co-occurrence and spatial autocorrelation to identify whether local transmission was a primary driver of pathogen distribution. Our surveys found widespread prevalence of Deformed wing virus (DWV), Sacbrood virus (SBV), and Black queen cell virus (BQCV). BQCV and SBV were most prolific in commercial apiaries, while Chronic bee paralysis virus (CPBV) was more common in hobbyist apiaries than commercial apiaries. DWV was most common in urban landscapes and was best predicted by mite prevalence and julian date, while the incidence of both SBV and BQCV were best predicted by regional apiary density. We did not find evidence of additional spatial autocorrelation for any viruses, although high co-occurrence suggests parallel transmission patterns. Our results support the importance of mite management in slowing virus spread and suggest that greater bee density increases transmission. Our study provides support that viruses are widespread in honey bees and connects known mechanisms of virus transmission to the distribution of pathogens observed across the Pacific Northwest.

Genetic Variation in Antimicrobial Activity of Honey Bee (Apis mellifera) Seminal Fluid

Honey bees can host a remarkably large number of different parasites and pathogens, and some are known drivers of recent declines in wild and managed bee populations. Here, we studied the interactions between the fungal pathogen Nosema apis and seminal fluid of the Western honey bee (Apis mellifera). Honey bee seminal fluid contains multiple antimicrobial molecules that kill N. apis spores and we therefore hypothesized that antimicrobial activities of seminal fluid are genetically driven by interactions between honey bee genotype and different N. apis strains/ecotypes, with the virulence of a strain depending on the genotype of their honey bee hosts. Among the antimicrobials, chitinases have been found in honey bee seminal fluid and have the predicted N. apis killing capabilities. We measured chitinase activity in the seminal fluid of eight different colonies. Our results indicate that multiple chitinases are present in seminal fluid, with activity significantly differing between genotypes. We therefore pooled equal numbers of N. apis spores from eight different colonies and exposed subsamples to seminal fluid samples from each of the colonies. We infected males from each colony with seminal fluid exposed spore samples and quantified N. apis infections after 6 days. We found that host colony had a stronger effect compared to seminal fluid treatment, and significantly affected host mortality, infection intensity and parasite prevalence. We also found a significant effect of treatment, as well as a treatment × colony interaction when our data were analyzed ignoring cage as a blocking factor. Our findings provide evidence that N. apis-honey bee interactions are driven by genotypic effects, which could be used in the future for breeding purposes of disease resistant or tolerant honey bee stock.

Simulated vector transmission differentially influences dynamics of two viral variants of deformed wing virus in honey bees (Apis mellifera)

Understanding how vectors alter the interactions between viruses and their hosts is a fundamental question in virology and disease ecology. In honey bees, transmission of deformed wing virus (DWV) by parasitic Varroa mites has been associated with elevated disease and host mortality, and Varroa transmission has been hypothesized to lead to increased viral titres or select for more virulent variants. Here, we mimicked Varroa transmission by serially passaging a mixed population of two DWV variants, A and B, by injection through in vitro reared honey bee pupae and tracking these viral populations through five passages. The DWV-A and DWV-B variant proportions shifted dynamically through passaging, with DWV-B outcompeting DWV-A after one passage, but levels of both variants becoming equivalent by Passage 5. Sequencing analysis revealed a dominant, recombinant DWV-B strain (DWV-A derived 5′ IRES region with the rest of the genome DWV-B), with low nucleotide diversity that decreased through passaging. DWV-A populations had higher nucleotide diversity compared to DWV-B, but this also decreased through passaging. Selection signatures were found across functional regions of the DWV-A and DWV-B genomes, including amino acid mutations in the putative capsid protein region. Simulated vector transmission differentially impacted two closely related viral variants which could influence viral interactions with the host, demonstrating surprising plasticity in vector-host-viral dynamics.

Review on mathematical modeling of honeybee population dynamics

Honeybees have an irreplaceable position in agricultural production and the stabilization of natural ecosystems. Unfortunately, honeybee populations have been declining globally. Parasites, diseases, poor nutrition, pesticides, and climate changes contribute greatly to the global crisis of honeybee colony losses. Mathematical models have been used to provide useful insights on potential factors and important processes for improving the survival rate of colonies. In this review, we present various mathematical tractable models from different aspects: 1) simple bee-only models with features such as age segmentation, food collection, and nutrient absorption; 2) models of bees with other species such as parasites and/or pathogens; and 3) models of bees affected by pesticide exposure. We aim to review those mathematical models to emphasize the power of mathematical modeling in helping us understand honeybee population dynamics and its related ecological communities. We also provide a review of computational models such as VARROAPOP and BEEHAVE that describe the bee population dynamics in environments that include factors such as temperature, rainfall, light, distance and quality of food, and their effects on colony growth and survival. In addition, we propose a future outlook on important directions regarding mathematical modeling of honeybees. We particularly encourage collaborations between mathematicians and biologists so that mathematical models could be more useful through validation with experimental data.

Honey bees increase social distancing when facing the ectoparasite Varroa destructor

Social distancing in response to infectious diseases is a strategy exhibited by human and nonhuman animals to counteract the spread of pathogens and/or parasites. Honey bee (Apis mellifera) colonies are ideal models to study this behavior because of the compartmentalized structure of these societies, evolved under exposure to parasite pressure and the need to ensure efficient functioning. Here, by using a combination of spatial and behavioral approaches, we investigated whether the presence of the ectoparasite mite Varroa destructor induces changes in the social organization of A. mellifera colonies that could reduce the spread of the parasite. Our results demonstrated that honey bees react to the intrusion of V. destructor by modifying space use and social interactions to increase the social distancing between young (nurses) and old (foragers) cohorts of bees. These findings strongly suggest a behavioral strategy not previously reported in honey bees to limit the intracolony parasite transmission.

A simple method for ex vivo honey bee cell culture capable of in vitro gene expression analysis

Cultured cells are a very powerful tool for investigating biological events in vitro; therefore, cell lines have been established not only in model insect species, but also in non-model species. However, there are few reports on the establishment of stable cell lines and development of systems to introduce genes into the cultured cells of the honey bee (Apis mellifera). We describe a simple ex vivo cell culture system for the honey bee. Hemocyte cells obtained from third and fourth instar larvae were cultured in commercial Grace’s insect medium or MGM-450 insect medium for more than two weeks maintaining a normal morphology without deterioration. After an expression plasmid vector bearing the enhanced green fluorescent protein (egfp) gene driven by the immediate early 2 (IE2) viral promoter was transfected into cells, EGFP fluorescence was detected in cells for more than one week from one day after transfection. Furthermore, double-stranded RNA corresponding to a part of the egfp gene was successfully introduced into cells and interfered with egfp gene expression. A convenient and reproducible method for an ex vivo cell culture that is fully practicable for gene expression assays was established for the honey bee.

Bee Products: A Representation of Biodiversity, Sustainability, and Health

Biodiversity strengthens the productivity of any ecosystem (agricultural land, forest, lake, etc.). The loss of biodiversity contributes to food and energy insecurity; increases vulnerability to natural disasters, such as floods or tropical storms; and decreases the quality of both life and health. Wild and managed bees play a key role in maintaining the biodiversity and in the recovery and restoration of degraded habitats. The novelty character of this perspective is to give an updated representation of bee products’ biodiversity, sustainability, and health relationship. The role of bees as bioindicators, their importance in the conservation of biodiversity, their ecosystem services, and the variety of the bee products are described herein. An overview of the main components of bee products, their biological potentials, and health is highlighted and detailed as follows: (i) nutritional value of bee products, (ii) bioactive profile of bee products and the related beneficial properties; (iii) focus on honey and health through a literature quantitative analysis, and (iv) bee products explored through databases. Moreover, as an example of the interconnection between health, biodiversity, and sustainability, a case study, namely the “Cellulose Park”, realized in Rome (Italy), is presented here. This case study highlights how bee activities can be used to assess and track changes in the quality of agricultural ecosystems—hive products could be valid indicators of the quality and health of the surrounding environment, as well as the changes induced by the biotic and abiotic factors that impact the sustainability of agricultural production and biodiversity conservation in peri-urban areas.

Pathogens Spillover from Honey Bees to Other Arthropods

Honey bees, and pollinators in general, play a major role in the health of ecosystems. There is a consensus about the steady decrease in pollinator populations, which raises global ecological concern. Several drivers are implicated in this threat. Among them, honey bee pathogens are transmitted to other arthropods populations, including wild and managed pollinators. The western honey bee, Apis mellifera, is quasi-globally spread. This successful species acted as and, in some cases, became a maintenance host for pathogens. This systematic review collects and summarizes spillover cases having in common Apis mellifera as the mainteinance host and some of its pathogens. The reports are grouped by final host species and condition, year, and geographic area of detection and the co-occurrence in the same host. A total of eighty-one articles in the time frame 1960-2021 were included. The reported spillover cases cover a wide range of hymenopteran host species, generally living in close contact with or sharing the same environmental resources as the honey bees. They also involve non-hymenopteran arthropods, like spiders and roaches, which are either likely or unlikely to live in close proximity to honey bees. Specific studies should consider host-dependent pathogen modifications and effects on involved host species. Both the plasticity of bee pathogens and the ecological consequences of spillover suggest a holistic approach to bee health and the implementation of a One Health approach.

A Case Report of Chronic Stress in Honey Bee Colonies Induced by Pathogens and Acaricide Residues

In this case report, we analyze the possible causes of the poor health status of a professional Apis mellifera iberiensis apiary located in Gajanejos (Guadalajara, Spain). Several factors that potentially favor colony collapse were identified, including Nosema ceranae infection, alone or in combination with other factors (e.g., BQCV and DWV infection), and the accumulation of acaricides commonly used to control Varroa destructor in the beebread (coumaphos and tau-fluvalinate). Based on the levels of residues, the average toxic unit estimated for the apiary suggests a possible increase in vulnerability to infection by N. ceranae due to the presence of high levels of acaricides and the unusual climatic conditions of the year of the collapse event. These data highlight the importance of evaluating these factors in future monitoring programs, as well as the need to adopt adequate preventive measures as part of national and international welfare programs aimed at guaranteeing the health and fitness of bees.

Investigating Virus-Host Interactions in Cultured Primary Honey Bee Cells

Honey bee (Apis mellifera) health is impacted by viral infections at the colony, individual bee, and cellular levels. To investigate honey bee antiviral defense mechanisms at the cellular level we further developed the use of cultured primary cells, derived from either larvae or pupae, and demonstrated that these cells could be infected with a panel of viruses, including common honey bee infecting viruses (i.e., sacbrood virus (SBV) and deformed wing virus (DWV)) and an insect model virus, Flock House virus (FHV). Virus abundances were quantified over the course of infection. The production of infectious virions in cultured honey bee pupal cells was demonstrated by determining that naïve cells became infected after the transfer of deformed wing virus or Flock House virus from infected cell cultures. Initial characterization of the honey bee antiviral immune responses at the cellular level indicated that there were virus-specific responses, which included increased expression of bee antiviral protein-1 (GenBank: MF116383) in SBV-infected pupal cells and increased expression of argonaute-2 and dicer-like in FHV-infected hemocytes and pupal cells. Additional studies are required to further elucidate virus-specific honey bee antiviral defense mechanisms. The continued use of cultured primary honey bee cells for studies that involve multiple viruses will address this knowledge gap.

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.

Reduced proinsecticide activation by cytochrome P450 confers coumaphos resistance in the major bee parasite Varroa destructor

Varroa destructor is one of the main problems in modern beekeeping. Highly selective acaricides with low toxicity to bees are used internationally to control this mite. One of the key acaricides is the organophosphorus (OP) proinsecticide coumaphos, that becomes toxic after enzymatic activation inside Varroa We show here that mites from the island Andros (AN-CR) exhibit high levels of coumaphos resistance. Resistance is not mediated by decreased coumaphos uptake, target-site resistance, or increased detoxification. Reduced proinsecticide activation by a cytochrome P450 enzyme was the main resistance mechanism, a powerful and rarely encountered evolutionary solution to insecticide selection pressure. After treatment with sublethal doses of [¹⁴C] coumaphos, susceptible mite extracts had substantial amounts of coroxon, the activated metabolite of coumaphos, while resistant mites had only trace amounts. This indicates a suppression of the P450 (CYP)-mediated activation step in the AN-CR mites. Bioassays with coroxon to bypass the activation step showed that resistance was dramatically reduced. There are 26 CYPs present in the V. destructor genome. Transcriptome analysis revealed overexpression in resistant mites of CYP4DP24 and underexpression of CYP3012A6 and CYP4EP4 RNA interference of CYP4EP4 in the susceptible population, to mimic underexpression seen in the resistant mites, prevented coumaphos activation and decreased coumaphos toxicity.

Breeding for Virus Resistance and Its Effects on Deformed Wing Virus Infection Patterns in Honey Bee Queens

Viruses, and in particular the deformed wing virus (DWV), are considered as one of the main antagonists of honey bee health. The 'suppressed in ovo virus infection' trait (SOV) described for the first time that control of a virus infection can be achieved from genetically inherited traits and that the virus state of the eggs is indicative for this. This research aims to explore the effect of the SOV trait on DWV infections in queens descending from both SOV-positive (QDS+) and SOV-negative (QDS-) queens. Twenty QDS+ and QDS- were reared from each time four queens in the same starter-finisher colony. From each queen the head, thorax, ovaries, spermatheca, guts and eviscerated abdomen were dissected and screened for the presence of the DWV-A and DWV-B genotype using qRT-PCR. Queens descending from SOV-positive queens showed significant lower infection loads for DWV-A and DWV-B as well as a lower number of infected tissues for DWV-A. Surprisingly, differences were less expressed in the reproductive tissues, the ovaries and spermatheca. These results confirm that selection on the SOV trait is associated with increased virus resistance across viral genotypes and that this selection drives DWV towards an increased tissue specificity for the reproductive tissues. Further research is needed to explore the mechanisms underlying the interaction between the antiviral response and DWV.

Transcriptome Profiling Reveals a Novel Mechanism of Antiviral Immunity Upon Sacbrood Virus Infection in Honey Bee Larvae (Apis cerana)

The honey bee is one of the most important pollinators in the agricultural system and is responsible for pollinating a third of all food we eat. Sacbrood virus (SBV) is a member of the virus family Iflaviridae and affects honey bee larvae and causes particularly devastating disease in the Asian honey bees, Apis cerana. Chinese Sacbrood virus (CSBV) is a geographic strain of SBV identified in China and has resulted in mass death of honey bees in China in recent years. However, the molecular mechanism underlying SBV infection in the Asian honey bee has remained unelucidated. In this present study, we employed high throughput next-generation sequencing technology to study the host transcriptional responses to CSBV infection in A. cerana larvae, and were able to identify genome-wide differentially expressed genes associated with the viral infection. Our study identified 2,534 differentially expressed genes (DEGs) involved in host innate immunity including Toll and immune deficiency (IMD) pathways, RNA interference (RNAi) pathway, endocytosis, etc. Notably, the expression of genes encoding antimicrobial peptides (abaecin, apidaecin, hymenoptaecin, and defensin) and core components of RNAi such as Dicer-like and Ago2 were found to be significantly upregulated in CSBV infected larvae. Most importantly, the expression of Sirtuin target genes, a family of signaling proteins involved in metabolic regulation, apoptosis, and intracellular signaling was found to be changed, providing the first evidence of the involvement of Sirtuin signaling pathway in insects’ immune response to a virus infection. The results obtained from this study provide novel insights into the molecular mechanism and immune responses involved in CSBV infection, which in turn will contribute to the development of diagnostics and treatment for the diseases in honey bees.

Africanized honey bees in Colombia exhibit high prevalence but low level of infestation of Varroa mites and low prevalence of pathogenic viruses

The global spread of the ectoparasitic mite Varroa destructor has promoted the spread and virulence of highly infectious honey bee viruses. This phenomenon is considered the leading cause for the increased number of colony losses experienced by the mite-susceptible European honey bee populations in the Northern hemisphere. Most of the honey bee populations in Central and South America are Africanized honey bees (AHBs), which are considered more resistant to Varroa compared to European honey bees. However, the relationship between Varroa levels and the spread of honey bee viruses in AHBs remains unknown. In this study, we determined Varroa prevalence and infestation levels as well as the prevalence of seven major honey bee viruses in AHBs from three regions of Colombia. We found that although Varroa exhibited high prevalence (92%), its infestation levels were low (4.5%) considering that these populations never received acaricide treatments. We also detected four viruses in the three regions analyzed, but all colonies were asymptomatic, and virus prevalence was considerably lower than those found in other countries with higher rates of mite-associated colony loss (DWV 19.88%, BQCV 17.39%, SBV 23.4%, ABPV 10.56%). Our findings indicate that AHBs possess a natural resistance to Varroa that does not prevent the spread of this parasite among their population, but restrains mite population growth and suppresses the prevalence and pathogenicity of mite-associated viruses.

Directed evolution of Metarhizium fungus improves its biocontrol efficacy against Varroa mites in honey bee colonies

Entomopathogenic fungi show great promise as pesticides in terms of their relatively high target specificity, low non-target toxicity, and low residual effects in agricultural fields and the environment. However, they also frequently have characteristics that limit their use, especially concerning tolerances to temperature, ultraviolet radiation, or other abiotic factors. The devastating ectoparasite of honey bees, Varroa destructor, is susceptible to entomopathogenic fungi, but the relatively warm temperatures inside honey bee hives have prevented these fungi from becoming effective control measures. Using a combination of traditional selection and directed evolution techniques developed for this system, new strains of Metarhizium brunneum were created that survived, germinated, and grew better at bee hive temperatures (35 °C). Field tests with full-sized honey bee colonies confirmed that the new strain JH1078 is more virulent against Varroa mites and controls the pest comparable to current treatments. These results indicate that entomopathogenic fungi are evolutionarily labile and capable of playing a larger role in modern pest management practices.

Virion structure and in vitro genome release mechanism of dicistrovirus Kashmir bee virus

Infections of Kashmir bee virus (KBV) are lethal for honeybees and have been associated with colony collapse disorder. KBV and closely related viruses contribute to the ongoing decline in the number of honeybee colonies in North America, Europe, Australia, and other parts of the world. Despite the economic and ecological impact of KBV, its structure and infection process remain unknown. Here we present the structure of the virion of KBV determined to a resolution of 2.8 Å. We show that the exposure of KBV to acidic pH induces a reduction in inter-pentamer contacts within capsids and the reorganization of its RNA genome from a uniform distribution to regions of high and low density. Capsids of KBV crack into pieces at acidic pH, resulting in the formation of open particles lacking pentamers of capsid proteins. The large openings of capsids enable the rapid release of genomes and thus limit the probability of their degradation by RNases. The opening of capsids may be a shared mechanism for the genome release of viruses from the family Dicistroviridae ImportanceThe western honeybee (Apis mellifera) is indispensable for maintaining agricultural productivity as well as the abundance and diversity of wild flowering plants. However, bees suffer from environmental pollution, parasites, and pathogens, including viruses. Outbreaks of virus infections cause the deaths of individual honeybees as well as collapses of whole colonies. Kashmir bee virus has been associated with colony collapse disorder in the US, and no cure of the disease is currently available. Here we report the structure of an infectious particle of Kashmir bee virus and show how its protein capsid opens to release the genome. Our structural characterization of the infection process determined that therapeutic compounds stabilizing contacts between pentamers of capsid proteins could prevent the genome release of the virus.