First Complete Genome of Lake Sinai Virus Lineage 3 and Genetic Diversity of Lake Sinai Virus Strains From Honey Bees and Bumble Bees

Varroa destructor shapes the unique viral landscapes of the honey bee populations of the Azores archipelago

The worldwide dispersal of the ectoparasitic mite Varroa destructor from its Asian origins has fundamentally transformed the relationship of the honey bee (Apis mellifera) with several of its viruses, via changes in transmission and/or host immunosuppression. The extent to which honey bee-virus relationships change after Varroa invasion is poorly understood for most viruses, in part because there are few places in the world with several geographically close but completely isolated honey bee populations that either have, or have not, been exposed long-term to Varroa, allowing for separate ecological, epidemiological, and adaptive relationships to develop between honey bees and their viruses, in relation to the mite's presence or absence. The Azores is one such place, as it contains islands with and without the mite. Here, we combined qPCR with meta-amplicon deep sequencing to uncover the relationship between Varroa presence, and the prevalence, load, diversity, and phylogeographic structure of eight honey bee viruses screened across the archipelago. Four viruses were not detected on any island (ABPV-Acute bee paralysis virus, KBV-Kashmir bee virus, IAPV-Israeli acute bee paralysis virus, BeeMLV-Bee macula-like virus); one (SBV-Sacbrood virus) was detected only on mite-infested islands; one (CBPV-Chronic bee paralysis virus) occurred on some islands, and two (BQCV-Black queen cell virus, LSV-Lake Sinai virus,) were present on every single island. This multi-virus screening builds upon a parallel survey of Deformed wing virus (DWV) strains that uncovered a remarkably heterogeneous viral landscape featuring Varroa-infested islands dominated by DWV-A and -B, Varroa-free islands naïve to DWV, and a refuge of the rare DWV-C dominating the easternmost Varroa-free islands. While all four detected viruses investigated here were affected by Varroa for one or two parameters (usually prevalence and/or the Richness component of ASV diversity), the strongest effect was observed for the multi-strain LSV. Varroa unambiguously led to elevated prevalence, load, and diversity (Richness and Shannon Index) of LSV, with these results largely shaped by LSV-2, a major LSV strain. Unprecedented insights into the mite-virus relationship were further gained from implementing a phylogeographic approach. In addition to enabling the identification of a novel LSV strain that dominated the unique viral landscape of the easternmost islands, this approach, in combination with the recovered diversity patterns, strongly suggests that Varroa is driving the evolutionary change of LSV in the Azores. This study greatly advances the current understanding of the effect of Varroa on the epidemiology and adaptive evolution of these less-studied viruses, whose relationship with Varroa has thus far been poorly defined.

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.

Prevalence and genome features of lake sinai virus isolated from Apis mellifera in the Republic of Korea

Lake Sinai Virus (LSV) is an emerging pathogen known to affect the honeybee (Apis mellifera). However, its prevalence and genomic characteristics in the Republic of Korea (ROK) remain unexplored. This study aimed to assess the prevalence of and analyze the LSVs by examining 266 honeybee samples from the ROK. Our findings revealed that LSV exhibited the highest infection rate among the pathogens observed in Korean apiaries, particularly during the reported period of severe winter loss (SWL) in A. mellifera apiaries in 2022. Three LSV genotypes- 2, 3, and 4 -were identified using RNA-dependent RNA polymerase gene analysis. Importantly, the infection rates of LSV2 (65.2%) and LSV3 (73.3%) were significantly higher in colonies experiencing SWL than in those experiencing normal winter loss (NWL) (p < 0.03). Furthermore, this study provides the first near-complete genome sequences of the Korean LSV2, LSV3, and LSV4 strains, comprising 5,759, 6,040, and 5,985 nt, respectively. Phylogenetic analysis based on these near-complete genome sequences demonstrated a close relationship between LSVs in the ROK and China. The high LSV infection rate in colonies experiencing a heightened mortality rate during winter suggests that this pathogen might contribute to SWL in ROK. Moreover, the genomic characteristic information on LSVs in this study holds immense potential for epidemiological information and the selection of specific genes suitable for preventing and treating LSV, including the promising utilization of RNA interference medicine in the future.

Preliminary Survey of Pathogens in the Asian Honey Bee (Apis cerana) in Thailand

Widespread parasites, along with emerging threats, globalization, and climate change, have greatly affected honey bees' health, leading to colony losses worldwide. In this study, we investigated the detection of biotic stressors (i.e., viruses, microsporidian, bacteria, and fungi) in Apis cerana by surveying the colonies across different regions of Thailand (Chiang Mai in the north, Nong Khai and Khon Kaen in the northeast, and Chumphon and Surat Thani in the south, in addition to the Samui and Pha-ngan islands). In this study, we detected ABPV, BQCV, LSV, and Nosema ceranae in A. cerana samples through RT-PCR. ABPV was only detected from the samples of Chiang Mai, whereas we found BQCV only in those from Chumphon. LSV was detected only in the samples from the Samui and Pha-ngan islands, where historically no managed bees are known. Nosema ceranae was found in all of the regions except for Nong Khai and Khon Kaen in northeastern Thailand. Paenibacillus larvae and Ascosphaera apis were not detected in any of the A. cerana samples in this survey. The phylogenetic tree analysis of the pathogens provided insights into the pathogens' movements and their distribution ranges across different landscapes, indicating the flow of pathogens among the honey bees. Here, we describe the presence of emerging pathogens in the Asian honey bee as a valuable step in our understanding of these pathogens in terms of the decline in eastern honey bee populations.

Structures of honeybee-infecting Lake Sinai virus reveal domain functions and capsid assembly with dynamic motions

Understanding the structural diversity of honeybee-infecting viruses is critical to maintain pollinator health and manage the spread of diseases in ecology and agriculture. We determine cryo-EM structures of T = 4 and T = 3 capsids of virus-like particles (VLPs) of Lake Sinai virus (LSV) 2 and delta-N48 LSV1, belonging to tetraviruses, at resolutions of 2.3-2.6 Å in various pH environments. Structural analysis shows that the LSV2 capsid protein (CP) structural features, particularly the protruding domain and C-arm, differ from those of other tetraviruses. The anchor loop on the central β-barrel domain interacts with the neighboring subunit to stabilize homo-trimeric capsomeres during assembly. Delta-N48 LSV1 CP interacts with ssRNA via the rigid helix α1', α1'-α1 loop, β-barrel domain, and C-arm. Cryo-EM reconstructions, combined with X-ray crystallographic and small-angle scattering analyses, indicate that pH affects capsid conformations by regulating reversible dynamic particle motions and sizes of LSV2 VLPs. C-arms exist in all LSV2 and delta-N48 LSV1 VLPs across varied pH conditions, indicating that autoproteolysis cleavage is not required for LSV maturation. The observed linear domino-scaffold structures of various lengths, made up of trapezoid-shape capsomeres, provide a basis for icosahedral T = 4 and T = 3 architecture assemblies. These findings advance understanding of honeybee-infecting viruses that can cause Colony Collapse Disorder.

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.

Two Distinct Genomic Lineages of Sinaivirus Detected in Guyanese Africanized Honey Bees

Over the past decade or so, PCR-based screening programs have reported that Africanized honey bees (AHB) are also hosts to viruses commonly found in European honey bees. Very little is known about the genomic variants found in AHB. Here, we present two distinct lineages of sinaiviruses in AHB.

First detection of Lake Sinai virus in the Czech Republic: a potential member of a new species

Lake Sinai virus (LSV) is one of over 20 honey bee viruses. Variants of LSV have been classified as members of two officially recognised species, Lake Sinai virus 1 and Lake Sinai virus 2. However, there are currently a limited number of whole-genome sequences, and the genetic variability of the virus indicates that additional species may need to be established. Extracted nucleic acid of 209 honey bee samples was screened by PCR for 11 honey bee viruses. LSV was the third most abundant virus (36.9% of positive samples), after Apis mellifera filamentous virus (72.2%) and deformed wing virus (52.5%). LSV-positive samples were analyzed further by PCR with primers targeting the region encoding the viral RNA-dependent RNA polymerase. Subsequently, the PCR products were sequenced, and the resulting sequences were used for a first round of phylogenetic analysis. Based on those results, several isolates were selected for whole-genome sequencing, and the complete genome sequences were used for additional phylogenetic analysis. The results indicated the presence of at least three genetically distinct groups of LSV in the Czech Republic, the most prevalent one being related to LSV 2 but too dissimilar to be considered a member of the same species. Two sequences of a major LSV strain cluster native to the Czech Republic were determined, representing the first Czech LSV strains published to date.

First detection of Lake Sinai virus in honeybees (Apis mellifera) and wild arthropods in Japan

Lake Sinai virus (LSV), an RNA virus, is suspected to be associated with poor health in honeybees (Apis mellifera). We examined LSV in 26 specimens of healthy honeybees and 44 specimens of wild arthropods in the Gifu Prefecture, Japan. LSV was found more frequently in honeybee specimens (11/26, 42.3%) than in wild arthropod specimens (1/44, 2.3%) (P<0.01). Phylogenetic and nucleotide sequence analysis revealed two lineages: LSV3 in honeybees, and LSV4 in both honeybees and wild arthropods. To our knowledge, this is the first report of LSV prevalence in honeybees and wild arthropods in Japan.

The Comparison of Honeybee Viral Loads for Six Honeybee Viruses (ABPV, BQCV, CBPV, DWV, LSV3 and SBV) in Healthy and Clinically Affected Honeybees with TaqMan Quantitative Real-Time RT-PCR Assays

The viral loads of acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), deformed wing virus (DWV), Lake Sinai virus 3 (LSV3), and sacbrood bee virus (SBV) were determined in samples with the use of quantitative TaqMan real-time reverse transcription and polymerase chain reaction (RT-qPCR). A total of 108 samples of healthy adult honeybees from four differently located apiaries and samples of honeybees showing different clinical signs of viral infections from 89 apiaries were collected throughout Slovenia. The aim of this study was to discover correlations between viral loads and clinical signs in adult honeybees and confirm previously set threshold viral load levels between healthy and clinically affected honeybees. Within this study, two new RT-qPCR assays for quantification of LSV3 and SBV were developed. Statistically significant differences in viral loads of positive samples were identified between healthy and clinically affected honeybees for ABPV, CBPV, DWV, and SBV, while for BQCV and LSV3, no statistical differences were observed between both groups. Despite high detected LSV3 prevalence and viral loads around 6.00 log₁₀ viral copies/bee, this lineage probably has a limited impact on the health status of honeybee colonies. The determined viral loads between 3.94 log₁₀ and 13.17 log₁₀ in positive samples for six viruses, collected over 10 consecutive months, including winter, present additional information of high viral load variations in healthy honeybee colonies.

Determination of Genetically Identical Strains of Four Honeybee Viruses in Bumblebee Positive Samples

In recent years, there has been growing evidence that certain types of honeybee viruses could be transmitted between different pollinators. Within a voluntary monitoring programme, 180 honeybee samples (Apis mellifera carnica) were collected from affected apiaries between 2007 and 2018. Also from August 2017 to August 2018, a total 148 samples of healthy bumblebees (Bombus lapidarius, B. pascuorum, B. terrestris, B. lucorum, B. hortorum, B. sylvarum, B. humilis) were collected at four different locations in Slovenia, and all samples were tested by using RT-PCR methods for six honeybee viruses. Direct sequencing of a total 158 positive samples (acute bee paralysis virus (ABPV n = 33), black queen cell virus (BQCV n = 75), sacbrood bee virus (SBV n = 25) and Lake Sinai virus (LSV n = 25)) was performed from obtained RT-PCR products. The genetic comparison of identified positive samples of bumblebees and detected honeybee field strains of ABPV, BQCV, SBV, and LSV demonstrated 98.74% to 100% nucleotide identity between both species. This study not only provides evidence that honeybees and bumblebees are infected with genetically identical or closely related viral strains of four endemically present honeybee viruses but also detected a high diversity of circulating strains in bumblebees, similar as was observed among honeybees. Important new genetic data for endemic strains circulating in honeybees and bumblebees in Slovenia are presented.