Homology differences between complete Sacbrood virus genomes from infected Apis mellifera and Apis cerana honeybees in Korea

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.

Sacbrood viruses cross-infection between Apis cerana and Apis mellifera: Rapid detection, viral dynamics, evolution and spillover risk assessment

Recent outbreaks of sacbrood virus (SBV) have caused serious epizootic disease in Apis cerana populations across Asia including Taiwan. Earlier phylogenetic analyses showed that cross-infection of AcSBV and AmSBV in both A. cerana and A. mellifera seems common, raising a concern of cross-infection intensifying the risk of disease resurgence in A. cerana. In this study, we analyzed the dynamics of cross-infection in three different types of apiaries (A. mellifera-only, A. cerana-only and two species co-cultured apiaries) over one year in Taiwan. Using novel, genotype-specific primer sets, we showed that SBV infection status varies across apiaries: AmSBV-AM and AcSBV-AC were the major genotype in the A. mellifera-only and the A. cerana-only apiaries, respectively, while AmSBV-AC and AcSBV-AC were the dominant genotypes in the co-cultured apiaries. Interestingly, co-cultured apiaries were among the only apiary type that harbored all variants and dual infections (i.e., AC and AM genotype co-infection in a single sample), indicating the interactions between hosts may form a conduit for cross-infection. The cross-infection between the two honey bee species appears to occur in a regular cycle with temporal fluctuation of AmSBV-AC and AcSBV-AC prevalence synchronized to each other in the co-cultured apiaries. Artificial infection of AcSBV in A. mellifera workers showed the suppression of viral replication, suggesting the potential of A. mellifera serving as a AcSBV reservoir that may contribute to virus spillover. Furthermore, the survival rate of A. cerana larvae was significantly reduced after artificial infections of both SBVs, indicating fitness costs of cross-infection on A. cerana and thus a high risk of disease resurgence in co-cultured apiaries. Our field and laboratory data provide baseline information that facilitates understanding of the risk of SBV cross-infection, and highlights the urgent need of SBV monitoring in co-cultured apiaries.

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.

Genomic Sequencing and Comparison of Sacbrood Viruses from Apis cerana and Apis mellifera in Taiwan

Sacbrood virus (SBV) was the first identified bee virus and shown to cause serious epizootic infections in the population of Apis cerana in Taiwan in 2015. Herein, the whole genome sequences of SBVs in A. cerana and A. mellifera were decoded and designated AcSBV-TW and AmSBV-TW, respectively. The whole genomes of AcSBV-TW and AmSBV-TW were 8776 and 8885 bp, respectively, and shared 90% identity. Each viral genome encoded a polyprotein, which consisted of 2841 aa in AcSBV-TW and 2859 aa in AmSBV-TW, and these sequences shared 95% identity. Compared to 54 other SBVs, the structural protein and protease regions showed high variation, while the helicase was the most highly conserved region among SBVs. Moreover, a 17-amino-acid deletion was found in viral protein 1 (VP1) region of AcSBV-TW compared to AmSBV-TW. The phylogenetic analysis based on the polyprotein sequences and partial VP1 region indicated that AcSBV-TW was grouped into the SBV clade with the AC-genotype (17-aa deletion) and was closely related to AmSBV-SDLY and CSBV-FZ, while AmSBV-TW was grouped into the AM-genotype clade but branched independently from other AmSBVs, indicating that the divergent genomic characteristics of AmSBV-TW might be a consequence of geographic distance driving evolution, and AcSBV-TW was closely related to CSBV-FZ, which originated from China. This 17-amino-acid deletion could be found in either AcSBV or AmSBV in Taiwan, indicating cross-infection between the two viruses. Our data revealed geographic and host specificities between SBVs. The amino acid difference in the VP1 region might serve as a molecular marker for describing SBV cross-infection.

Phylogenetic analysis of sacbrood virus structural polyprotein and non-structural RNA dependent RNA polymerase gene: Differences in Turkish strains

Sacbrood virus (SBV) is one of the most damaging viruses in honey bee colonies. Genetic differences among sacbrood viruses detected in honey bees in different locales have been reported in previous studies. The aim of this study was to construct phylogenetic trees based on the structural polyprotein and non-structural RNA dependent RNA polymerase gene regions and to make a molecular characterization of the Tur/Bur/Sac01 and Tur/Bur/Sac02 strains identified in Apis mellifera in Turkey. As a result of the study, the tree based on the structural polyprotein region separated into four lineages: Tur/Bur/Sac01 and Tur/Bur/Sac02 were in the same branch as the Turkish sacbrood virus strains identified in previous studies and formed the Turkish clade. Strains isolated from adjacent geographical areas were in the same clade in this tree. The phylogenetic tree based on the non-structural RNA dependent RNA polymerase gene region divides into two main branches, reflecting host affiliation: Apis cerana and A. mellifera. Strains formed clusters based on their geographic distribution and host affiliation. The Tur/Bur/Sac01 and Tur/Bur/Sac02 strains formed a separate cluster among the European strains. Sacbrood viruses from Turkey were genetically different from SBV strains detected in other countries and in A. cerana.

Impact of Sacbrood Virus on Larval Microbiome of Apis mellifera and Apis cerana

In this study, we examined the impact of Sacbrood virus (SBV), the cause of larval honeybee (Apis mellifera) death, producing a liquefied a larva sac, on the gut bacterial communities on two larval honeybee species, Apis mellifera and Apis cerana. SBV was added into a worker jelly food mixture and bee larvae were grafted into each of the treatment groups for 24 h before DNA/RNA extraction. Confirmation of SBV infection was achieved using quantitative reverse transcription polymerase chain reaction (RT-qPCR) and visual symptomology. The 16S rDNA was sequenced by Illumina sequencing. The results showed the larvae were infected with SBV. The gut communities of infected A. cerana larvae exhibited a dramatic change compared with A. mellifera. In A. mellifera larvae, the Illumina sequencing revealed the proportion of Gilliamella, Snodgrassella and Fructobacillus was not significantly different, whereas in A. cerana, Gilliamella was significantly decreased (from 35.54% to 2.96%), however, with significant increase in Snodgrassella and Fructobacillus. The possibility of cross-infection should be further investigated.

Genetic and phylogenetic analysis of Chinese sacbrood virus isolates from Apis mellifera

Background

Sacbrood virus (SBV) is one of the most pathogenic honeybee viruses that exhibits host specificity and regional variations. The SBV strains that infect the Chinese honeybee Apis cerana are called Chinese SBVs (CSBVs).

Methods

In this study, a CSBV strain named AmCSBV-SDLY-2016 (GenBank accession No. MG733283) infecting A. mellifera was identified by electron microscopy, its protein composition was analyzed by sodium dodecyl sulfate polyacrylamide gel electrophoresis and agar gel immunodiffusion assay, and its nucleotide sequence was identified using a series of reverse-transcription polymerase chain reaction fragments of AmCSBV-SDLY-2016 generated using SBV/CSBV-specific primers. To investigate phylogenetic relationships of the CSBV isolates, a phylogenetic tree of the complete open reading frames (ORF) of the CSBV sequences was constructed using MEGA 6.0; then, the similarity and recombination events among the isolated CSBV strains were analyzed using SimPlot and RDP4 software, respectively.

Results

Sequencing results revealed the complete 8,794-nucleotide long complete genomic RNA of the strain, with a single large ORF (189–8,717) encoding 2,843 amino acids. Comparison of the deduced amino acid sequence with the SBV/CSBV reference sequences deposited in the GenBank database identified helicase, protease, and RNA-dependent RNA polymerase domains; the structural genes were located at the 5′ end, whereas the non-structural genes were found at the 3′ end. Multiple sequence alignment showed that AmCSBV-SDLY-2016 had a 17-amino acid (aa) and a single aa deletion at positions 711–729 and 2,128, respectively, as compared with CSBV-GD-2002, and a 16-aa deletion (positions 711–713 and 715–728) as compared with AmSBV-UK-2000. However, AmCSBV-SDLY-2016 was similar to the CSBV-JLCBS-2014 strain, which infects A. cerana. AmCSBV-SDLY-2016 ORF shared 92.4–97.1% identity with the genomes of other CSBV strains (94.5–97.7% identity for deduced amino acids). AmCSBV-SDLY-2016 was least similar (89.5–90.4% identity) to other SBVs but showed maximum similarity with the previously reported CSBV-FZ-2014 strain. The phylogenetic tree constructed from AmCSBV-SDLY-2016 and 43 previously reported SBV/CSBV sequences indicated that SBV/CSBV strains clustered according to the host species and country of origin; AmCSBV-SDLY-2016 clustered with other previously reported Chinese and Asian strains (AC genotype SBV, as these strains originated from A. cerana) but was separate from the SBV genomes originating from Europe (AM genotype SBV, originating from A. mellifera). A SimPlot graph of SBV genomes confirmed the high variability, especially between the AC genotype SBV and AM genotype SBV. This genomic diversity may reflect the adaptation of SBV to specific hosts, ability of CSBV to cross the species barrier, and the spatial distances that separate CSBVs from other SBVs.

Complete genome sequence of sacbrood virus isolated from Asiatic honey bee Apis cerana indica in India

We determined the complete genome sequence of a sacbrood virus (SBV) infecting Indian honey bee (Apis cerana indica) from Tamil Nadu, India named as AcSBV-IndTN1. The genome of AcSBV-IndTN1 comprised of 8740 nucleotides, encoding a single large ORF containing 2849 amino acids flanked by 5' and 3' untranslated regions. Results of phylogenetic tree analysis based on complete genomes of SBV isolates indicated that the virus isolates from India isolated from the Asiatic honey bee A. cerana (AcSBVs) formed a separate group along with six Vietnam isolates and three Chinese isolates. The AcSBV-IndTN1 isolate showed closer genetic relationship with other isolates from India. The second major group had both AcSBVs and AmSBVs (virus isolated from European honey bee, Apis mellifera SBV) of Korea, China and Vietnam. The third and a distantly related group had AmSBVs of Australia, UK, USA and Korea. The results obtained from phylogenetic analysis were further supported with evolutionary distance analysis. AcSBV-IndTN1 isolate open reading frame had 95-99% amino acid sequence similarity with other Indian isolates and 92-96% with AcSBVs and AmSBVs of other geographical locations. In addition, sequence difference count matrix ranged from 154 to 907 nt among all the SBV isolates. This suggests that the virus isolates have evolved significantly in different geographical locations but isolates on different hosts in a given location/country are closely related. The high similarity in the genome among the AcSBV and AmSBV isolates indicate possible cross-infections and recombination of SBV isolates in Asian continent where both the honey bee species are reared in close proximity. Gene flow between SBV population indicating that an infrequent gene flow occur between them. The pattern of molecular diversity in SBV population revealed that the occurrence of recent population expansion of SBV. To the best of our knowledge this is the first report of the complete nucleotide sequence of AcSBV from Tamil Nadu, India. This study provided an opportunity to establish the molecular evolution of SBV isolates and shall be useful in the development of diagnostics and effective disease control strategies.

Molecular detection and phylogenetic assessment of six honeybee viruses in Apis mellifera L. colonies in Bulgaria

Honey bee colonies suffer from various pathogens, including honey bee viruses. About 24 viruses have been reported so far. However, six of them are considered to cause severe infection which inflicts heavy losses on beekeeping. The aim of this study was to investigate incidence of six honey bee 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) by a reverse transcription polymerase chain reaction (RT-PCR). A total of 250 adult honey bee samples were obtained from 50 colonies from eight apiaries situated in three different parts of the country (South, North and West Bulgaria). The results showed the highest prevalence of DWV followed by SBV and ABPV, and one case of BQCV. A comparison with homology sequences available in GenBank was performed by phylogenetic analysis, and phylogenetic relationships were discussed in the context of newly described genotypes in the uninvestigated South Eastern region of Europe. In conclusion, the present study has been the first to provide sequencing data and phylogenetics analyses of some honey bee viruses in Bulgaria.

Comparative Genomic Analysis for Genetic Variation in Sacbrood Virus of Apis cerana and Apis mellifera Honeybees From Different Regions of Vietnam

Sacbrood virus (SBV) is one of the most common viral infections of honeybees. The entire genome sequence for nine SBV infecting honeybees, Apis cerana and Apis mellifera, in Vietnam, namely AcSBV-Viet1, AcSBV-Viet2, AcSBV-Viet3, AmSBV-Viet4, AcSBV-Viet5, AmSBV-Viet6, AcSBV-Viet7, AcSBV-Viet8, and AcSBV-Viet9, was determined. These sequences were aligned with seven previously reported complete genome sequences of SBV from other countries, and various genomic regions were compared. The Vietnamese SBVs (VN-SBVs) shared 91-99% identity with each other, and shared 89-94% identity with strains from other countries. The open reading frames (ORFs) of the VN-SBV genomes differed greatly from those of SBVs from other countries, especially in their VP1 sequences. The AmSBV-Viet6 and AcSBV-Viet9 genome encodes 17 more amino acids within this region than the other VN-SBVs. In a phylogenetic analysis, the strains AmSBV-Viet4, AcSBV-Viet2, and AcSBV-Viet3 were clustered in group with AmSBV-UK, AmSBV-Kor21, and AmSBV-Kor19 strains. Whereas, the strains AmSBV-Viet6 and AcSBV-Viet7 clustered separately with the AcSBV strains from Korea and AcSBV-VietSBM2. And the strains AcSBV-Viet8, AcSBV-Viet1, AcSBV-Viet5, and AcSBV-Viet9 clustered with the AcSBV-India, AcSBV-Kor and AcSBV-VietSBM2. In a Simplot graph, the VN-SBVs diverged stronger in their ORF regions than in their 5' or 3' untranslated regions. The VN-SBVs possess genetic characteristics which are more similar to the Asian AcSBV strains than to AmSBV-UK strain. Taken together, our data indicate that host specificity, geographic distance, and viral cross-infections between different bee species may explain the genetic diversity among the VN-SBVs in A. cerana and A. mellifera and other SBV strains.

A comparison of biological characteristics of three strains of Chinese sacbrood virus in Apis cerana

We selected and sequenced the entire genomes of three strains of Chinese sacbrood virus (CSBV): LNQY-2008 (isolated in Qingyuan, Liaoning Province), SXYL-2015 (isolated in Yulin, Shanxi Province), and JLCBS-2014 (isolated in Changbaishan, Jilin Province), by VP1 amino acid (aa) analysis. These strains are endemic in China and infect Apis cerana. Nucleotide sequences, deduced amino acid sequences, genetic backgrounds, and other molecular biological characteristics were analysed. We also examined sensitivity of these virus strains to temperature, pH, and organic solvents, as well as to other physicochemical properties. On the basis of these observations, we compared pathogenicity and tested cross-immunogenicity and protective immunity, using antisera raised against each of the three strains. Our results showed that compared with SXYL-2015, LNQY-2008 has a 10-aa deletion and 3-aa deletion (positions 282–291 and 299–301, respectively), whereas JLCBS-2014 has a 17-aa deletion (positions 284–300). However, the three strains showed no obvious differences in physicochemical properties or pathogenicity. Moreover, there was immune cross-reactivity among the antisera raised against the different strains, implying good protective effects of such antisera. The present study should significantly advance the understanding of the pathogenesis of Chinese sacbrood disease, and offers insights into comprehensive prevention and treatment of, as well as possible protection from, the disease by means of an antiserum.