Occurrence of multiple honeybee viruses in the ectoparasitic mites Varroa spp. in Apis cerana colonies

Israeli Acute Paralysis Virus Is an Emerging Pathogen Contributing to Brood Disease of Apis cerana

Larval mortality is the primary symptom of diseased Apis cerana colonies, often attributed to sacbrood virus (SBV) and Melissococcus plutonius. However, the impact of other common honeybee viruses is frequently overlooked, and their pathogenicity to A. cerana remains poorly understood. To investigate the causes of the increasing disease incidence in A. cerana brood, we conducted an epidemiological survey, collecting 70 samples from 19 sites across nine provinces in China. Furthermore, we examined the pathogenicity of Israeli acute paralysis virus (IAPV) in A. cerana brood through artificial inoculation experiments. Our results demonstrate that, besides SBV and M. plutonius, the infection rate and viral load of IAPV in diseased brood are significantly high. Brood artificially inoculated with high concentrations of IAPV exhibited a significant increase in mortality and displayed clinical symptoms similar to those observed in naturally infected colonies. Moreover, a limited resistance to IAPV was observed in A. cerana brood, with some individuals able to restrict viral proliferation. Our study highlights the previously unrecognized pathogenicity of IAPV to A. cerana brood, demonstrating that IAPV poses a significant threat similar to SBV and M. plutonius. We emphasize that IAPV should be recognized as an emerging pathogen causing brood disease in A. cerana and managed accordingly in beekeeping practices.

First detection and prevalence of Apis mellifera filamentous virus in Apis mellifera and Varroa destructor in the Republic of Korea

Apis mellifera filamentous virus (AmFV) is a double-stranded DNA virus that infects Apis mellifera bees. To our knowledge, this is the first comprehensive study aiming to detect and analyse the genetic diversity and prevalence of AmFV in Korean honeybee colonies. Phylogenetic analysis based on baculovirus repeat open reading frame-N gene (Bro-N) sequences revealed that AmFV isolates from the Republic of Korea (ROK) fell into two distinct lineages, with genetic origins in Switzerland and China, with nucleotide similarities of 98.3% and 98.2%, respectively. Our prevalence analysis demonstrated a noteworthy infection rate of AmFV in 545 honeybee colonies, reaching 33.09% in 2022 and increasing to 44.90% by 2023. Intriguingly, we also detected AmFV in Varroa destructor mites, highlighting their potential role as vectors and carriers of AmFV. The presence of AmFV was correlated with an increased infection rate of sacbrood virus, deformed wing virus, Lake Sinai virus 2, black queen cell virus, and Nosema ceranae in honeybee colonies. These findings provide valuable insight into the prevalence and potential transmission mechanisms of AmFV in honeybee colonies in the ROK. The results of this study may be instrumental in the effective management of viral infections in honeybee apiaries.

Biotic and abiotic stresses on honeybee health

Honeybees are the most critical pollinators providing key ecosystem services that underpin crop production and sustainable agriculture. Amidst a backdrop of rapid global change, this eusocial insect encounters a succession of stressors during nesting, foraging, and pollination. Ectoparasitic mites, together with vectored viruses, have been recognized as central biotic threats to honeybee health, while the spread of invasive giant hornets and small hive beetles also increasingly threatens colonies worldwide. Cocktails of agrochemicals, including acaricides used for mite treatment, and other pollutants of the environment have been widely documented to affect bee health in various ways. Additionally, expanding urbanization, climate change, and agricultural intensification often result in the destruction or fragmentation of flower-rich bee habitats. The anthropogenic pressures exerted by beekeeping management practices affect the natural selection and evolution of honeybees, and colony translocations facilitate alien species invasion and disease transmission. In this review, the multiple biotic and abiotic threats and their interactions that potentially undermine bee colony health are discussed, while taking into consideration the sensitivity, large foraging area, dense network among related nestmates, and social behaviors of honeybees.

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.

Molecular Detection and Phylogenetic Analysis of Deformed Wing Virus and Sacbrood Virus Isolated from Pollen

Among many pathogens and pests, honey bee viruses are known as one of the most common cause of diseases in honey bee colonies. In this study, we demonstrate that pollen grains and bee bread are potential sources of viral DNA. We extracted DNA from 3 types of pollen samples: directly provided by beekeepers (n = 12), purchased from trade markets (n = 5), and obtained from honeycombs (bee bread, n = 10). The extracted DNA was used for molecular detection (RT-PCR analysis) of six of the most widely distributed honey bee viruses: deformed wing virus, sacbrood virus, acute bee paralysis virus, black queen cell virus, Kashmir bee virus, Israeli acute paralysis virus, and chronic bee paralysis virus. We successfully managed to establish only the deformed wing virus (DWV) and the sacbrood virus (SBV), with different distribution frequencies depending on the territory of the country. The phylogenetic analyses of Bulgarian isolates were performed with the most similar sequences available in molecular databases from other countries. Phylogenies of Bulgarian viral strains demonstrated genetically heterogeneous populations of DWV and relatively homogenous populations of SBV. In conclusion, the results obtained from the current study have shown that pollen is a valuable source for molecular detection of honey bee pathogens. This allows epidemiological monitoring of honey bee diseases at a regional and a national level.

Genetic Relationships and Signatures of Adaptation to the Climatic Conditions in Populations of Apis cerana Based on the Polymorphism of the Gene Vitellogenin

Apis cerana and Apis mellifera are important honey bee species in Asia. A. cerana populations are distributed from a cold, sharply continental climate in the north to a hot, subtropical climate in the south. Due to the Sacbrood virus, almost all A. cerana populations in Asia have declined significantly in recent decades and have recovered over the past five years. This could lead to a shift in the gene pool of local A. cerana populations that could affect their sustainability and adaptation. It was assumed that adaptation of honey bees could be observed by comparative analysis of the sequences of genes involved in development, labor division, and caste differentiation, such as the gene Vitellogenin VG. The VG gene nucleotide sequences were used to assess the genetic structure and signatures of adaptation of local populations of A. cerana from Korea, Russia, Japan, Nepal, and China. A. mellifera samples from India and Poland were used as the outgroup. The signatures of adaptive selection were found in the local population of A. cerana using VG gene sequence analysis based on Jukes−Cantor genetic distances, cluster analysis, dN/dS ratio evaluation, and Tajima’s D neutrality test. Based on analysis of the VG gene sequences, Apis cerana koreana subspecies in the Korean Peninsula were subdivided into three groups in accordance with their geographic localization from north to south. The VG gene sequences are acceptable tools to study the sustainability and adaptation of A. cerana populations.

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.

Recovery and genetic characterization of black queen cell virus

Black queen cell virus (BQCV) is a severe threat to the honeybee (Apis mellifera) worldwide. Although several BQCV strains have been reported in China, the molecular basis for BQCV pathogenicity has not been well understood. Thus, a reverse genetic system of BQCV is required for studying viral replication and its pathogenic mechanism. Here, the complete genome sequence of BQCV was obtained from honeybees using reverse transcription PCR (RT-PCR), namely a BQCV China-GS1 strain (KY741959). Then, a phylogenetic tree was built to analyse the genetic relationships among BQCV strains from different regions. Our results showed that the BQCV China-GS1 contained two ORFs, consistent with the known reference strains, except for the BQCV China-JL1 strain (KP119603). Furthermore, the infectious clone of BQCV was constructed based on BQCV China-GS1 using a low copy vector pACYC177 and gene recombination. Due to the lack of culture cells for bee viruses, we infected the healthy bees with infectious clone of BQCV, and the rescued BQCV resulted in the recovery of recombinant virus, which induced higher mortality than those of the control group. Immune response after inoculated with BQCV further confirmed that the infectious clone of BQCV caused the cellular and humoral immune response of honeybee (A. mellifera). In conclusion, the full nucleotide sequence of BQCV China-GS1 strain was determined, and the infectious clone of BQCV was constructed in this study. These data will improve the understanding of pathogenesis and the host immune responses to viral infection.

A derived honey bee stock confers resistance to Varroa destructor and associated viral transmission

The ectoparasite Varroa destructor is the greatest threat to managed honey bee (Apis mellifera) colonies globally. Despite significant efforts, novel treatments to control the mite and its vectored pathogens have shown limited efficacy, as the host remains naïve. A prospective solution lies in the development of Varroa-resistant honey bee stocks, but a paucity of rigorous selection data restricts widespread adoption. Here, we characterise the parasite and viral dynamics of a Varroa-resistant honey bee stock, designated ‘Pol-line’, using a large-scale longitudinal study. Results demonstrate markedly reduced Varroa levels in this stock, diminished titres of three major viruses (DWV-A, DWV-B, and CBPV), and a two-fold increase in survival. Levels of a fourth virus that is not associated with Varroa—BQCV—do not differ between stocks, supporting a disruption of the transmission pathway. Further, we show that when decoupled from the influence of Varroa levels, viral titres do not constitute strong independent predictors of colony mortality risk. These findings highlight the need for a reassessment of Varroa etiology, and suggest that derived stocks represent a tractable solution to the Varroa pandemic.

Prevalence of honeybee viruses in Apis mellifera in Gifu prefecture of Japan

Viral infection damages honeybee colony health. Viruses can be carried by queen bees and apicultural production materials when imported from foreign countries. We investigated seven honeybee viruses in worker bees (Apis mellifera) from 26 healthy apiaries in Gifu, Japan between 2018 and 2019. Black queen cell virus (BQCV) was detected in 23 (88.5%) of the apiaries, followed by Israeli acute paralysis virus (42.3%), deformed wing virus (DWV) (38.5%), and sacbrood virus (3.8%). In phylogenetic analysis, BQCV and DWV in Gifu were related to those in China and South Korea. Additionally, a high prevalence of BQCV was observed among worker bees in BQCV-positive colonies. Therefore, BQCV horizontal transmission among worker bees may contribute to the high prevalence of BQCV in Gifu.

Tolerance of Honey Bees to Varroa Mite in the Absence of Deformed Wing Virus

The global spread of the parasitic mite Varroa destructor has emphasized the significance of viruses as pathogens of honey bee (Apis mellifera) populations. In particular, the association of deformed wing virus (DWV) with V. destructor and its devastating effect on honey bee colonies has led to that virus now becoming one of the most well-studied insect viruses. However, there has been no opportunity to examine the effects of Varroa mites without the influence of DWV. In Papua New Guinea (PNG), the sister species, V. jacobsoni, has emerged through a host-shift to reproduce on the local A. mellifera population. After initial colony losses, beekeepers have maintained colonies without chemicals for more than a decade, suggesting that this bee population has an unknown mite tolerance mechanism. Using high throughput sequencing (HTS) and target PCR detection, we investigated whether the viral landscape of the PNG honey bee population is the underlying factor responsible for mite tolerance. We found A. mellifera and A. cerana from PNG and nearby Solomon Islands were predominantly infected by sacbrood virus (SBV), black queen cell virus (BQCV) and Lake Sinai viruses (LSV), with no evidence for any DWV strains. V. jacobsoni was infected by several viral homologs to recently discovered V. destructor viruses, but Varroa jacobsoni rhabdovirus-1 (ARV-1 homolog) was the only virus detected in both mites and honey bees. We conclude from these findings that A. mellifera in PNG may tolerate V. jacobsoni because the damage from parasitism is significantly reduced without DWV. This study also provides further evidence that DWV does not exist as a covert infection in all honey bee populations, and remaining free of this serious viral pathogen can have important implications for bee health outcomes in the face of Varroa.