Honey Bee Infecting Lake Sinai Viruses

Honey bee viruses in the yellow-legged hornet Vespa velutina (Lepelieter 1836): Prevalence, loads, and detection of replicative DWV and LSV forms

Apiaries in Galicia, northwestern Spain, are currently facing the invasive alien species Vespa velutina, which is well established in the region. The pressure on honey bee colonies is high, resulting in both economic and ecological losses. Honey bee colonies also face the challenge of viruses, which are becoming increasingly diverse. In recent years, honey bee viruses have been spreading across taxonomic groups beyond Apoidea, infecting the Vespoidea superfamily. This cross-species spillover has raised concerns in the scientific community due to the potential risk of viruses spreading in ecosystems. Currently, there is a lack of knowledge on this topic, and further research is needed to address this issue. This study employed qPCR and sequencing to investigate the prevalence, loads, and presence of replicative forms of important honey bee viruses in V. velutina individuals collected from 11 apiaries in Galicia. All V. velutina individuals tested positive for DWV, BQCV, AKI complex (ABPV, KBV, and IAPV), or LSV but not for CBPV. DWV showed the highest prevalence (97.0 %) and loads, with both DWV-A (67.4 %) and DWV-B (32.6 %) being detected. The AKI complex (46.3 %) and LSV (43.3 %) were also common, whereas BQCV (11.9 %) was rarer. LSV is detected for the first time in V. velutina. LSV-2 was the dominant strain (82.1 %), and two less frequent (17.9 %) unknown strains were also detected. All 44 screened V. velutina samples carried the replicative form of DWV, and six of these also carried the replicative form of LSV, raising for the first time the possibility of co-infection in the hornet. The detection of honey bee viruses in V. velutina, and the ability of these viruses to spread to other species, may indicate a potential risk of spillover in the apiaries.

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.

Identification and pathogen detection of a Neocypholaelaps species (Acari: Mesostigmata: Ameroseiidae) from beehives in the Republic of Korea

In this study, we identified a new strain of the genus Neocypholaelaps from the beehives of Apis mellifera colonies in the Republic of Korea (ROK). The Neocypholaelap sp. KOR23 mites were collected from the hives of honeybee apiaries in Wonju, Gangwon-do, in May 2023. Morphological and molecular analyses based on 18S and 28S rRNA gene regions conclusively identified that these mites belong to the genus Neocypholaelaps, closely resembling Neocypholaelaps sp. APGD-2010 that was first isolated from the United States. The presence of 9 of 25 honeybee pathogens in these mite samples suggests that Neocypholaelaps sp. KOR23 mite may act as an intermediate vector and carrier of honeybee diseases. The identification of various honeybee pathogens within this mite highlights their significance in disease transmission among honeybee colonies. This comprehensive study provides valuable insights into the taxonomy and implications of these mites for bee health management and pathogen dissemination.

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.

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.

WIV, a protein domain found in a wide number of arthropod viruses, which probably facilitates infection

The most powerful approach to detect distant homologues of a protein is based on structure prediction and comparison. Yet this approach is still inapplicable to many viral proteins. Therefore, we applied a powerful sequence-based procedure to identify distant homologues of viral proteins. It relies on three principles: (1) traces of sequence similarity can persist beyond the significance cutoff of homology detection programmes; (2) candidate homologues can be identified among proteins with weak sequence similarity to the query by using 'contextual' information, e.g. taxonomy or type of host infected; (3) these candidate homologues can be validated using highly sensitive profile-profile comparison. As a test case, this approach was applied to a protein without known homologues, encoded by ORF4 of Lake Sinai viruses (which infect bees). We discovered that the ORF4 protein contains a domain that has homologues in proteins from >20 taxa of viruses infecting arthropods. We called this domain 'widespread, intriguing, versatile' (WIV), because it is found in proteins with a wide variety of functions and within varied domain contexts. For example, WIV is found in the NSs protein of tospoviruses, a global threat to food security, which infect plants as well as their arthropod vectors; in the RNA2 ORF1-encoded protein of chronic bee paralysis virus, a widespread virus of bees; and in various proteins of cypoviruses, which infect the silkworm Bombyx mori. Structural modelling with AlphaFold indicated that the WIV domain has a previously unknown fold, and bibliographical evidence suggests that it facilitates infection of arthropods.

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.

Quantitative trait loci mapping for survival of virus infection and virus levels in honey bees

Israeli acute paralysis virus (IAPV) is a highly virulent, Varroa-vectored virus that is of global concern for honey bee health. Little is known about the genetic basis of honey bees to withstand infection with IAPV or other viruses. We set up and analyzed a backcross between preselected honey bee colonies of low and high IAPV susceptibility to identify quantitative trait loci (QTL) associated with IAPV susceptibility. Experimentally inoculated adult worker bees were surveyed for survival and selectively sampled for QTL analysis based on SNPs identified by whole-genome resequencing and composite interval mapping. Additionally, natural titers of other viruses were quantified in the abdomen of these workers via qPCR and also used for QTL mapping. In addition to the full dataset, we analyzed distinct subpopulations of susceptible and non-susceptible workers separately. These subpopulations are distinguished by a single, suggestive QTL on chromosome 6, but we identified numerous other QTL for different abdominal virus titers, particularly in the subpopulation that was not susceptible to IAPV. The pronounced QTL differences between the susceptible and non-susceptible subpopulations indicate either an interaction between IAPV infection and the bees' interaction with other viruses or heterogeneity among workers of a single cohort that manifests itself as IAPV susceptibility and results in distinct subgroups that differ in their interaction with other viruses. Furthermore, our results indicate that low susceptibility of honey bees to viruses can be caused by both, virus tolerance and virus resistance. QTL were partially overlapping among different viruses, indicating a mixture of shared and specific processes that control viruses. Some functional candidate genes are located in the QTL intervals, but their genomic co-localization with numerous genes of unknown function delegates any definite characterization of the underlying molecular mechanisms to future studies.

Exploring viral infections in honey bee colonies: insights from a metagenomic study in southern Brazil

The decline in honey bee colonies in different parts of the world in recent years is due to different reasons, such as agricultural practices, climate changes, the use of chemical insecticides, and pests and diseases. Viral infections are one of the main causes leading to honey bee population declines, which have a major economic impact due to honey production and pollination. To investigate the presence of viruses in bees in southern Brazil, we used a metagenomic approach to sequence adults' samples of concentrated extracts from Apis mellifera collected in fifteen apiaries of six municipalities in the Rio Grande do Sul state, Brazil, between 2016 and 2017. High-throughput sequencing (HTS) of these samples resulted in the identification of eight previously known viruses (Apis rhabdovirus 1 (ARV-1), Acute bee paralysis virus (ABPV), Aphid lethal paralysis virus (ALPV), Black queen cell virus (BQCV), Bee Macula-like virus (BeeMLV), Deformed wing virus (DWV), Lake Sinai Virus NE (LSV), and Varroa destructor virus 3 (VDV-3)) and a thogotovirus isolate. This thogotovirus shares high amino acid identities in five of the six segments with Varroa orthomyxovirus 1, VOV-1 (98.36 to 99.34% identity). In contrast, segment 4, which codes for the main glycoprotein (GP), has no identity with VOV-1, as observed for the other segments, but shares an amino acid identity of 34-38% with other glycoproteins of viruses from the Orthomyxoviridae family. In addition, the putative thogotovirus GP also shows amino acid identities ranging from 33 to 41% with the major glycoprotein (GP64) of insect viruses of the Baculoviridae family. To our knowledge, this is the second report of a thogotovirus found in bees and given this information, this thogotovirus isolate was tentatively named Apis thogotovirus 1 (ATHOV-1). The detection of multiple viruses in bees is important to better understand the complex interactions between viruses and their hosts. By understanding these interactions, better strategies for managing viral infections in bees and protecting their populations can be developed.

Virome Analysis Provides an Insight into the Viral Community of Chinese Mitten Crab Eriocheir sinensis

Recent advances in viromics have led to the discovery of a great diversity of RNA viruses and the identification of a large number of viral pathogens. A systematic exploration of viruses in Chinese mitten crab (Eriocheir sinensis), one of the most important aquatic commercial species, is still lacking. Here, we characterized the RNA viromes of asymptomatic, milky disease (MD)-affected, and hepatopancreatic necrosis syndrome (HPNS)-affected Chinese mitten crabs collected from 3 regions in China. In total, we identified 31 RNA viruses belonging to 11 orders, 22 of which were first reported here. By comparing viral composition between samples, we observed high variation in viral communities across regions, with most of the viral species being region-specific. We proposed to establish several novel viral families or genera based on the phylogenetic relationships and genome structures of viruses discovered in this study, expanding our knowledge of viral diversity in brachyuran crustaceans. IMPORTANCE High-throughput sequencing and meta-transcriptomic analysis provide us with an efficient tool to discover unknown viruses and explore the composition of viral communities in specific species. In this study, we investigated viromes in asymptomatic and diseased Chinese mitten crabs collected from three distant locations. We observed high regional variation in the composition of viral species, highlighting the importance of multi-location sampling. In addition, we classified several novel and ICTV-unclassified viruses based on their genome structures and phylogenetic relationships, providing a new perspective on current viral taxa.

Area Wide Monitoring of Plant and Honey Bee (Apis mellifera) Viruses in Blueberry (Vaccinium corymbosum) Agroecosystems Facilitated by Honey Bee Pollination

Healthy agroecosystems are dependent on a complex web of factors and inter-species interactions. Flowers are hubs for pathogen transmission, including the horizontal or vertical transmission of plant-viruses and the horizontal transmission of bee-viruses. Pollination by the European honey bee (Apis mellifera) is critical for industrial fruit production, but bees can also vector viruses and other pathogens between individuals. Here, we utilized commercial honey bee pollination services in blueberry (Vaccinium corymbosum) farms for a metagenomics-based bee and plant virus monitoring system. Following RNA sequencing, viruses were identified by mapping reads to a reference sequence database through the bioinformatics portal Virtool. In total, 29 unique plant viral species were found at two blueberry farms in British Columbia (BC). Nine viruses were identified at one site in Ontario (ON), five of which were not identified in BC. Ilarviruses blueberry shock virus (BlShV) and prune dwarf virus (PDV) were the most frequently detected viruses in BC but absent in ON, while nepoviruses tomato ringspot virus and tobacco ringspot virus were common in ON but absent in BC. BlShV coat protein (CP) nucleotide sequences were nearly identical in all samples, while PDV CP sequences were more diverse, suggesting multiple strains of PDV circulating at this site. Ten bee-infecting viruses were identified, with black queen cell virus frequently detected in ON and BC. Area-wide bee-mediated pathogen monitoring can provide new insights into the diversity of viruses present in, and the health of, bee-pollination ecosystems. This approach can be limited by a short sampling season, biased towards pollen-transmitted viruses, and the plant material collected by bees can be very diverse. This can obscure the origin of some viruses, but bee-mediated virus monitoring can be an effective preliminary monitoring approach.

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.

Omics-based analysis of honey bee (Apis mellifera) response to Varroa sp. parasitisation and associated factors reveals changes impairing winter bee generation

The extensive annual loss of honey bees (Apis mellifera L.) represents a global problem affecting agriculture and biodiversity. The parasitic mite Varroa destructor, associated with viral co-infections, plays a key role in this loss. Despite years of intensive research, the complex mechanisms of Varroa - honey bee interaction are still not fully defined. Therefore, this study employed a unique combination of transcriptomic, proteomic, metabolomic, and functional analyses to reveal new details about the effect of Varroa mites and naturally associated factors, including viruses, on honey bees. We focused on the differences between Varroa parasitised and unparasitised ten-day-old worker bees collected before overwintering from the same set of colonies reared without anti-mite treatment. Supplementary comparison to honey bees collected from colonies with standard anti-Varroa treatment can provide further insights into the effect of a pyrethroid flumethrin. Analysis of the honey bees exposed to mite parasitisation revealed alterations in the transcriptome and proteome related to immunity, oxidative stress, olfactory recognition, metabolism of sphingolipids, and RNA regulatory mechanisms. The immune response and sphingolipid metabolism were strongly activated, whereas olfactory recognition and oxidative stress pathways were inhibited in Varroa parasitised honey bees compared to unparasitised ones. Moreover, metabolomic analysis confirmed the depletion of nutrients and energy stores, resulting in a generally disrupted metabolism in the parasitised workers. The combined omics-based analysis conducted on strictly parasitised bees revealed the key molecular components and mechanisms underlying the detrimental effects of Varroa sp. and its associated pathogens. This study provides the theoretical basis and interlinked datasets for further research on honey bee response to biological threats and the development of efficient control strategies against Varroa mites.

Vector-virus interaction affects viral loads and co-occurrence

BACKGROUND

Vector-borne viral diseases threaten human and wildlife worldwide. Vectors are often viewed as a passive syringe injecting the virus. However, to survive, replicate and spread, viruses must manipulate vector biology. While most vector-borne viral research focuses on vectors transmitting a single virus, in reality, vectors often carry diverse viruses. Yet how viruses affect the vectors remains poorly understood. Here, we focused on the varroa mite (Varroa destructor), an emergent parasite that can carry over 20 honey bee viruses, and has been responsible for colony collapses worldwide, as well as changes in global viral populations. Co-evolution of the varroa and the viral community makes it possible to investigate whether viruses affect vector gene expression and whether these interactions affect viral epidemiology.

RESULTS

Using a large set of available varroa transcriptomes, we identified how abundances of individual viruses affect the vector's transcriptional network. We found no evidence of competition between viruses, but rather that some virus abundances are positively correlated. Furthermore, viruses that are found together interact with the vector's gene co-expression modules in similar ways, suggesting that interactions with the vector affect viral epidemiology. We experimentally validated this observation by silencing candidate genes using RNAi and found that the reduction in varroa gene expression was accompanied by a change in viral load.

CONCLUSIONS

Combined, the meta-transcriptomic analysis and experimental results shed light on the mechanism by which viruses interact with each other and with their vector to shape the disease course.

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.

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.

The Virome of Healthy Honey Bee Colonies: Ubiquitous Occurrence of Known and New Viruses in Bee Populations

The Western honey bee, Apis mellifera , is a vital part of our ecosystem as well as cultural heritage. Annual colony losses endanger beekeeping.

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.

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.

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.

Global similarity, and some key differences, in the metagenomes of Swedish varroa-surviving and varroa-susceptible honeybees

There is increasing evidence that honeybees (Apis mellifera L.) can adapt naturally to survive Varroa destructor, the primary cause of colony mortality world-wide. Most of the adaptive traits of naturally varroa-surviving honeybees concern varroa reproduction. Here we investigate whether factors in the honeybee metagenome also contribute to this survival. The quantitative and qualitative composition of the bacterial and viral metagenome fluctuated greatly during the active season, but with little overall difference between varroa-surviving and varroa-susceptible colonies. The main exceptions were Bartonella apis and sacbrood virus, particularly during early spring and autumn. Bombella apis was also strongly associated with early and late season, though equally for all colonies. All three affect colony protein management and metabolism. Lake Sinai virus was more abundant in varroa-surviving colonies during the summer. Lake Sinai virus and deformed wing virus also showed a tendency towards seasonal genetic change, but without any distinction between varroa-surviving and varroa-susceptible colonies. Whether the changes in these taxa contribute to survival or reflect demographic differences between the colonies (or both) remains unclear.

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.

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.

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.

Honey Bee Health in Maine Wild Blueberry Production

Simple Summary

Wild blueberry is an important native North American crop that requires insect pollination. Migratory western honey bee colonies constitute the majority of commercial bees brought into Maine for pollination of wild blueberry. Currently, many stressors impact the western honey bee in the US. We designed a two-year monitoring study (2014 and 2015) to assess the potential health of honey bee colonies hired for pollination services in wild blueberry fields. We monitored the colony health of nine hive locations (three hives/location) in 2014 and nine locations (five hives/location) in 2015 during bloom (May–June). Queen health status, colony strength, rate of population increase, and pesticide residues on pollen, wax, and honey bee workers were measured. In addition, each hive was sampled to assess levels of mite parasites, viruses, and Microsporidian and Trypanosome pathogens. Different patterns in colony health were observed over the two years. Factors predicting colony growth rate over both years were Varroa mite infestation and risk due to pollen pesticide residues during bloom. In addition, recently discovered parasites and pathogens were already observed in most of the colonies suggesting that parasites and diseases spread rapidly and become established quickly in commercial honey bee colonies.

Abstract

A two-year study was conducted in Maine wild blueberry fields (Vaccinium angustifolium Aiton) on the health of migratory honey bee colonies in 2014 and 2015. In each year, three or five colonies were monitored at each of nine wild blueberry field locations during bloom (mid-May until mid-June). Colony health was measured by assessing colony strength during wild blueberry bloom. Potential factors that might affect colony health were queen failure or supersedure; pesticide residues on trapped pollen, wax comb, and bee bread; and parasites and pathogens. We found that Varroa mite and pesticide residues on trapped pollen were significant predictors of colony health measured as the rate of change in the amount of sealed brood during bloom. These two factors explained 71% of the variance in colony health over the two years. Pesticide exposure was different in each year as were pathogen prevalence and incidence. We detected high prevalence and abundance of two recently discovered pathogens and one recently discovered parasite, the trypanosome Lotmaria passim Schwartz, the Sinai virus, and the phorid fly, Apocephalus borealis Brues.

Differential Viral Distribution Patterns in Reproductive Tissues of Apis mellifera and Apis cerana Drones

Honeybee drones are male bees that mate with virgin queens during the mating flight, consequently transferring their genes to offspring. Therefore, the health of drones affects the overall fitness of the offspring and ultimately the survivability of the colony. Honeybee viruses are considered to be a major threat to the health of honeybees. In the present study, we demonstrated the pattern of common honeybee viruses in various tissues of drones in the western honeybee, Apis mellifera, and the eastern honeybee, Apis cerana. Drones were collected during the mating flight and analyzed using quantitative real-time (qRT-PCR) to detect the presence of seven honeybee viruses. The qRT-PCR result revealed that three honeybee viruses, namely Black Queen Cell Virus (BQCV), Deformed Wing Virus (DWV), and Chinese Sacbrood Virus (CSBV), were detected in the reproductive tissues of A. mellifera and A. cerana drones. The results from qRT-PCR showed that the Israeli Acute Paralysis Virus (IAPV) was only detected in A. mellifera drone body tissues. Moreover, the prevalence of DWV and BQCV in the drones collected from A. mellifera colonies was significantly higher than that of A. cerana. In addition, virus multiple infections were higher in A. mellifera drones compared to those in A. cerana. CSBV was found predominantly in the reproductive tissues of A. cerana drones. This study is the first report describing the presence of the CSBV in reproductive tissues of A. mellifera drones. Our results may reflect the preference of honeybee viruses in honeybee species and may provide a piece of interesting evidence for understanding the virus transmission in A. cerana.

Viruses that affect Argentinian honey bees (Apis mellifera)

Beekeeping is a widespread activity in Argentina, mainly producing honey that has gained both national and international recognition. There are more than 3,000,000 hives in the country, mainly concentrated in Buenos Aires Province (approximately 1,000,000 hives). In recent decades, worrying rates of hive loss have been observed in many countries around the world. In Latin America, the estimated loss of hives is between 13% (Peru and Ecuador) and 53% (Chile). Argentina had annual losses of 34% for the period of October 1, 2016 to October 1, 2017. The causes of these losses are not clear but probably involve multiple stressors that can act simultaneously. One of the main causes of loss of bee colonies worldwide is infestation by the ectoparasitic mite Varroa destructor in combination with viral infections. To date, 10 viruses have been detected that affect honey bees (Apis mellifera) in Argentina. Of these, deformed wing virus, sacbrood virus, acute bee paralysis virus, chronic bee paralysis virus, and Israeli acute bee paralysis can be transmitted by mites. Deformed wing virus and the AIK complex are the viruses most often associated with loss of hives worldwide. Considering that bee viruses have been detected in Argentina in several hymenopteran and non-hymenopteran insects, these hosts could act as important natural reservoirs for viruses and play an important role in their dispersal in the environment. Further studies to investigate the different mechanisms by which viruses spread in the environment will enable us to develop various strategies for the control of infected colonies and the spread of viruses in the habitat where they are found.

Metatranscriptome Analysis of Sympatric Bee Species Identifies Bee Virus Variants and a New Virus, Andrena-Associated Bee Virus-1

Bees are important plant pollinators in agricultural and natural ecosystems. High average annual losses of honey bee (Apis mellifera) colonies in some parts of the world, and regional population declines of some mining bee species (Andrena spp.), are attributed to multiple factors including habitat loss, lack of quality forage, insecticide exposure, and pathogens, including viruses. While research has primarily focused on viruses in honey bees, many of these viruses have a broad host range. It is therefore important to apply a community level approach in studying the epidemiology of bee viruses. We utilized high-throughput sequencing to evaluate viral diversity and viral sharing in sympatric, co-foraging bees in the context of habitat type. Variants of four common viruses (i.e., black queen cell virus, deformed wing virus, Lake Sinai virus 2, and Lake Sinai virus NE) were identified in honey bee and mining bee samples, and the high degree of nucleotide identity in the virus consensus sequences obtained from both taxa indicates virus sharing. We discovered a unique bipartite + ssRNA Tombo-like virus, Andrena-associated bee virus-1 (AnBV-1). AnBV-1 infects mining bees, honey bees, and primary honey bee pupal cells maintained in culture. AnBV-1 prevalence and abundance was greater in mining bees than in honey bees. Statistical modeling that examined the roles of ecological factors, including floral diversity and abundance, indicated that AnBV-1 infection prevalence in honey bees was greater in habitats with low floral diversity and abundance, and that interspecific virus transmission is strongly modulated by the floral community in the habitat. These results suggest that land management strategies that aim to enhance floral diversity and abundance may reduce AnBV-1 spread between co-foraging bees.

Persistent effects of management history on honeybee colony virus abundances

Infectious diseases are a major threat to both managed and wild pollinators. One key question is how the movement or transplantation of honeybee colonies under different management regimes affects honeybee disease epidemiology. We opportunistically examined any persistent effect of colony management history following relocation by characterising the virus abundances of honeybee colonies from three management histories, representing different management histories: feral, low-intensity management, and high-intensity "industrial" management. The colonies had been maintained for one year under the same approximate 'common garden' condition. Colonies in this observational study differed in their virus abundances according to management history, with the feral population history showing qualitatively different viral abundance patterns compared to colonies from the two managed population management histories; for example, higher abundance of sacbrood virus but lower abundances of various paralysis viruses. Colonies from the high-intensity management history exhibited higher viral abundances for all viruses than colonies from the low-intensity management history. Our results provide evidence that management history has persistent impacts on honeybee disease epidemiology, suggesting that apicultural intensification could be majorly impacting on pollinator health, justifying much more substantial investigation.

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.

Colony-Level Effects of Amygdalin on Honeybees and Their Microbes

Amygdalin, a cyanogenic glycoside, is found in the nectar and pollen of almond trees, as well as in a variety of other crops, such as cherries, nectarines, apples and others. It is inevitable that western honeybees (Apis mellifera) consistently consume amygdalin during almond pollination season because almond crops are almost exclusively pollinated by honeybees. This study tests the effects of a field-relevant concentration of amygdalin on honeybee microbes and the activities of key honeybee genes. We executed a two-month field trial providing sucrose solutions with or without amygdalin ad libitum to free-flying honeybee colonies. We collected adult worker bees at four time points and used RNA sequencing technology and our HoloBee database to assess global changes in microbes and honeybee transcripts. Our hypothesis was that amygdalin will negatively affect bee microbes and possibly immune gene regulation. Using a log2 fold-change cutoff at two and intraday comparisons, we show no large change of bacterial counts, fungal counts or key bee immune gene transcripts, due to amygdalin treatment in relation to the control. However, relatively large titer decreases in the amygdalin treatment relative to the control were found for several viruses. Chronic bee paralysis virus levels had a sharp decrease (-14.4) with titers then remaining less than the control, Black queen cell virus titers were lower at three time points (<-2) and Deformed wing virus titers were lower at two time points (<-6) in amygdalin-fed compared to sucrose-fed colonies. Titers of Lotmaria passim were lower in the treatment group at three of the four dates (<-4). In contrast, Sacbrood virus had two dates with relative increases in its titers (>2). Overall, viral titers appeared to fluctuate more so than bacteria, as observed by highly inconstant patterns between treatment and control and throughout the season. Our results suggest that amygdalin consumption may reduce several honeybee viruses without affecting other microbes or colony-level expression of immune genes.

Longitudinal monitoring of honey bee colonies reveals dynamic nature of virus abundance and indicates a negative impact of Lake Sinai virus 2 on colony health

Honey bees (Apis mellifera) are important pollinators of plants, including those that produce nut, fruit, and vegetable crops. Therefore, high annual losses of managed honey bee colonies in the United States and many other countries threaten global agriculture. Honey bee colony deaths have been associated with multiple abiotic and biotic factors, including pathogens, but the impact of virus infections on honey bee colony population size and survival are not well understood. To further investigate seasonal patterns of pathogen presence and abundance and the impact of viruses on honey bee colony health, commercially managed colonies involved in the 2016 California almond pollination event were monitored for one year. At each sample date, colony health and pathogen burden were assessed. Data from this 50-colony cohort study illustrate the dynamic nature of honey bee colony health and the temporal patterns of virus infection. Black queen cell virus, deformed wing virus, sacbrood virus, and the Lake Sinai viruses were the most readily detected viruses in honey bee samples obtained throughout the year. Analyses of virus prevalence and abundance revealed pathogen-specific trends including the overall increase in deformed wing virus abundance from summer to fall, while the levels of Lake Sinai virus 2 (LSV2) decreased over the same time period. Though virus prevalence and abundance varied in individual colonies, analyses of the overall trends reveal correlation with sample date. Total virus abundance increased from November 2015 (post-honey harvest) to the end of the almond pollination event in March 2016, which coincides with spring increase in colony population size. Peak total virus abundance occurred in late fall (August and October 2016), which correlated with the time period when the majority of colonies died. Honey bee colonies with larger populations harbored less LSV2 than weaker colonies with smaller populations, suggesting an inverse relationship between colony health and LSV2 abundance. Together, data from this and other longitudinal studies at the colony level are forming a better understanding of the impact of viruses on honey bee colony losses.

An updated genetic marker for detection of Lake Sinai Virus and metagenetic applications

Background

Lake Sinai Viruses (LSV) are common RNA viruses of honey bees (Apis mellifera) that frequently reach high abundance but are not linked to overt disease. LSVs are genetically heterogeneous and collectively widespread, but despite frequent detection in surveys, the ecological and geographic factors structuring their distribution in A. mellifera are not understood. Even less is known about their distribution in other species. Better understanding of LSV prevalence and ecology have been hampered by high sequence diversity within the LSV clade.

Methods

Here we report a new polymerase chain reaction (PCR) assay that is compatible with currently known lineages with minimal primer degeneracy, producing an expected 365 bp amplicon suitable for end-point PCR and metagenetic sequencing. Using the Illumina MiSeq platform, we performed pilot metagenetic assessments of three sample sets, each representing a distinct variable that might structure LSV diversity (geography, tissue, and species).

Results

The first sample set in our pilot assessment compared cDNA pools from managed A. mellifera hives in California (n = 8) and Maryland (n = 6) that had previously been evaluated for LSV2, confirming that the primers co-amplify divergent lineages in real-world samples. The second sample set included cDNA pools derived from different tissues (thorax vs. abdomen, n = 24 paired samples), collected from managed A. mellifera hives in North Dakota. End-point detection of LSV frequently differed between the two tissue types; LSV metagenetic composition was similar in one pair of sequenced samples but divergent in a second pair. Overall, LSV1 and intermediate lineages were common in these samples whereas variants clustering with LSV2 were rare. The third sample set included cDNA from individual pollinator specimens collected from diverse landscapes in the vicinity of Lincoln, Nebraska. We detected LSV in the bee Halictus ligatus (four of 63 specimens tested, 6.3%) at a similar rate as A. mellifera (nine of 115 specimens, 7.8%), but only one H. ligatus sequencing library yielded sufficient data for compositional analysis. Sequenced samples often contained multiple divergent LSV lineages, including individual specimens. While these studies were exploratory rather than statistically powerful tests of hypotheses, they illustrate the utility of high-throughput sequencing for understanding LSV transmission within and among species.

The Pathogen Profile of a Honey Bee Queen Does Not Reflect That of Her Workers

Throughout a honey bee queen's lifetime, she is tended to by her worker daughters, who feed and groom her. Such interactions provide possible horizontal transmission routes for pathogens from the workers to the queen, and as such a queen's pathogen profile may be representative of the workers within a colony. To explore this further, we investigated known honey bee pathogen co-occurrence, as well as pathogen transmission from workers to queens. Queens from 42 colonies were removed from their source hives and exchanged into a second, unrelated foster colony. Worker samples were taken from the source colony on the day of queen exchange and the queens were collected 24 days after introduction. All samples were screened for Nosema spp., Trypanosome spp., acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), Israeli acute paralysis virus (IAPV), Lake Sinai virus (LSV), and deformed wing virus master variants (DWV-A, B, and C) using RT-qPCR. The data show that LSV, Nosema, and DWV-B were the most abundant pathogens in colonies. All workers (n = 42) were LSV-positive, 88% were Nosema-positive, whilst pathogen loads were low (<1 × 106 genome equivalents per pooled worker sample). All queens (n = 39) were negative for both LSV and Nosema. We found no evidence of DWV transmission occurring from worker to queen when comparing queens to foster colonies, despite DWV being present in both queens and workers. Honey bee pathogen presence and diversity in queens cannot be revealed from screening workers, nor were pathogens successfully transmitted to the queen.

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

The complete genome of Lake Sinai virus 3 (LSV3) was sequenced by the Ion Torrent next-generation sequencing (NGS) technology from an archive sample of honey bees collected in 2010. This strain M92/2010 is the first complete genome sequence of LSV lineage 3. From October 2016 to December 2017, 56 honey bee samples from 32 different locations and 41 bumble bee samples from five different locations were collected. These samples were tested using a specific reverse transcriptase-polymerase chain reaction (RT-PCR) method; 75.92% of honey bee samples and 17.07% of bumble bee samples were LSV-positive with the RT-PCR method. Phylogenetic comparison of 557-base pair-long RNA-dependent RNA polymerase (RdRp) genome region of selected 23 positive samples of honey bees and three positive bumble bee samples identified three different LSV lineages: LSV1, LSV2, and LSV3. The LSV3 lineage was confirmed for the first time in Slovenia in 2010, and the same strain was later detected in several locations within the country. The LSV strains detected in bumble bees are from 98.6 to 99.4% identical to LSV strains detected among honey bees in the same territory.

Bee Viruses: Routes of Infection in Hymenoptera

Numerous studies have recently reported on the discovery of bee viruses in different arthropod species and their possible transmission routes, vastly increasing our understanding of these viruses and their distribution. Here, we review the current literature on the recent advances in understanding the transmission of viruses, both on the presence of bee viruses in Apis and non-Apis bee species and on the discovery of previously unknown bee viruses. The natural transmission of bee viruses will be discussed among different bee species and other insects. Finally, the research potential of in vivo (host organisms) and in vitro (cell lines) serial passages of bee viruses is discussed, from the perspective of the host-virus landscape changes and potential transmission routes for emerging bee virus infections.

The Heat Shock Response in the Western Honey Bee (Apis mellifera) is Antiviral

Honey bees (Apismellifera) are an agriculturally important pollinator species that live in easily managed social groups (i.e., colonies). Unfortunately, annual losses of honey bee colonies in many parts of the world have reached unsustainable levels. Multiple abiotic and biotic stressors, including viruses, are associated with individual honey bee and colony mortality. Honey bees have evolved several antiviral defense mechanisms including conserved immune pathways (e.g., Toll, Imd, JAK/STAT) and dsRNA-triggered responses including RNA interference and a non-sequence specific dsRNA-mediated response. In addition, transcriptome analyses of virus-infected honey bees implicate an antiviral role of stress response pathways, including the heat shock response. Herein, we demonstrate that the heat shock response is antiviral in honey bees. Specifically, heat-shocked honey bees (i.e., 42 °C for 4 h) had reduced levels of the model virus, Sindbis-GFP, compared with bees maintained at a constant temperature. Virus-infection and/or heat shock resulted in differential expression of six heat shock protein encoding genes and three immune genes, many of which are positively correlated. The heat shock protein encoding and immune gene transcriptional responses observed in virus-infected bees were not completely recapitulated by administration of double stranded RNA (dsRNA), a virus-associated molecular pattern, indicating that additional virus-host interactions are involved in triggering antiviral stress response pathways.

Floral Traits Predict Frequency of Defecation on Flowers by Foraging Bumble Bees

Flowers may become inoculated with pathogens that can infect bees and other critical pollinators, but the mechanisms of inoculation remain unclear. During foraging, bees may regurgitate or defecate directly onto flower parts, which could inoculate flowers with pollinator pathogens and lead to subsequent disease transmission to floral visitors. We tested if captive eastern bumble bees (Bombus impatiens) (Cresson) (Hymenoptera: Apidae) defecate on floral surfaces during foraging and if flower shape played a role in the probability of defecation and the quantity of feces deposited on floral surfaces. Captive Bombus impatiens were fed a solution of fluorescent dye and sucrose, then allowed to forage freely on flowers of a variety of shapes in a flight cage. Flowers were then examined under ultraviolet light for fluorescing fecal matter. We found that bumble bees did defecate on floral surfaces during foraging and that composite flowers with a large area of disk flowers were the most likely to have feces on them. Our results point to defecation by bumble bees during foraging as a potential mechanism for inoculation of flowers with pollinator pathogens and suggest that flower shape could play a significant role in inoculation.

The virome of an endangered stingless bee suffering from annual mortality in southern Brazil

Meliponiculture - the management of stingless bee colonies - is an expanding activity in Brazil with economic, social and environmental potential. However, unlike in apiculture, the pathogens that impact on meliponiculture remain largely unknown. In southern Brazil, every year at the end of the summer, managed colonies of the stingless bee Melipona quadrifasciata manifest a syndrome that eventually leads to collapse. Here we characterize the M. quadrifasciata virome using high-throughput sequencing, with the aim of identifying potentially pathogenic viruses, and test whether they are related to the syndrome outbreaks. Two paired viromes are explored, one from healthy bees and another from unhealthy ones. Each virome is built from metagenomes assembled from sequencing reads derived either from RNA or DNA. A total of 40 621 reads map to viral contigs of the unhealthy bees' metagenomes, whereas only 11 reads map to contigs identified as viruses of healthy bees. The viruses showing the largest copy numbers in the virome of unhealthy bees belong to the family Dicistroviridae - common pathogenic honeybee viruses - as well as Parvoviridae and Circoviridae, which have never been reported as being pathogenic in insects. Our analyses indicate that they represent seven novel viruses associated with stingless bees. PCR-based detection of these viruses in individual bees (healthy or unhealthy) from three different localities revealed a statistically significant association between viral infection and symptom manifestation in one meliponary. We conclude that although viral infections may contribute to colony collapses in the annual syndrome in some meliponaries, viruses spread opportunistically during the outbreak, perhaps due to colony weakness.

Multiple Virus Infections in Western Honeybee (Apis mellifera L.) Ejaculate Used for Instrumental Insemination

Instrumental insemination of Apis mellifera L. queens is a widely employed technique used in honeybee breeding that enables the effective control of mating. However, drone semen represents a potential source of honeybee viruses. In this study, 43 semen doses collected from apparently healthy drones, and consequently used in instrumental insemination, were analysed using PCR or RT-PCR to detect the presence of viral genome of 11 honeybee viruses. In 91% of samples, viral infection was detected. The survey revealed genomes of five viruses, namely Deformed wing virus (DWV), Acute bee paralysis virus (ABPV), Black queen cell virus (BQCV), Sacbrood virus (SBV), and A. mellifera filamentous virus (AmFV) in 84%, 19%, 14%, 2%, and 67% of samples, respectively. Single infection (30% of samples) as well as multiple infection (61% of samples) of two, three or four pathogens were also evaluated. To the best of our knowledge, this is the first study describing the presence of the BQCV and SBV genome sequence in drone ejaculate. Phylogenetic analysis of BQCV partial helicase gene sequence revealed the high similarity of nucleotide sequence of described Czech strains, which varied from 91.4% to 99.6%. The findings of our study indicate the possibility of venereal transmission of BQCV and SBV.

Relative abundance and molecular evolution of Lake Sinai Virus (Sinaivirus) clades

Lake Sinai Viruses (Sinaivirus) are commonly detected in honey bees (Apis mellifera) but no disease phenotypes or fitness consequences have yet been demonstrated. This viral group is genetically diverse, lacks obvious geographic structure, and multiple lineages can co-infect individual bees. While phylogenetic analyses have been performed, the molecular evolution of LSV has not been studied extensively. Here, I use LSV isolates from GenBank as well as contigs assembled from honey bee Sequence Read Archive (SRA) accessions to better understand the evolutionary history of these viruses. For each ORF, substitution rate variation, codon usage, and tests of positive selection were evaluated. Outlier regions of high or low diversity were sought with sliding window analysis and the role of recombination in creating LSV diversity was explored. Phylogenetic analysis consistently identified two large clusters of sequences that correspond to the current LSV1 and LSV2 nomenclature, however lineages sister to LSV1 were the most frequently detected in honey bee SRA accessions. Different expression levels among ORFs suggested the occurrence of subgenomic transcripts. ORF1 and RNA-dependent RNA polymerase had higher evolutionary rates than the capsid and ORF4. A hypervariable region of the ORF1 protein-coding sequence was identified that had reduced selective constraint, but a site-based model of positive selection was not significantly more likely than a neutral model for any ORF. The only significant recombination signals detected between LSV1 and LSV2 initiated within this hypervariable region, but assumptions of the test (single-frame coding and independence of substitution rate by site) were violated. LSV codon usage differed strikingly from that of honey bees and other common honey-bee viruses, suggesting LSV is not strongly co-evolved with that host. LSV codon usage was significantly correlated with that of Varroa destructor, however, despite the relatively weak codon bias exhibited by the latter. While codon usage between the LSV1 and LSV2 clusters was similar for three ORFs, ORF4 codon usage was uncorrelated between these clades, implying rapid divergence of codon use for this ORF only. Phylogenetic placement and relative abundance of LSV isolates reconstructed from SRA accessions suggest that detection biases may be over-representing LSV1 and LSV2 in public databases relative to their sister lineages.

Bee Viruses: Ecology, Pathogenicity, and Impacts

Bees-including solitary, social, wild, and managed species-are key pollinators of flowering plant species, including nearly three-quarters of global food crops. Their ecological importance, coupled with increased annual losses of managed honey bees and declines in populations of key wild species, has focused attention on the factors that adversely affect bee health, including viral pathogens. Genomic approaches have dramatically expanded understanding of the diversity of viruses that infect bees, the complexity of their transmission routes-including intergenus transmission-and the diversity of strategies bees have evolved to combat virus infections, with RNA-mediated responses playing a prominent role. Moreover, the impacts of viruses on their hosts are exacerbated by the other major stressors bee populations face, including parasites, poor nutrition, and exposure to chemicals. Unraveling the complex relationships between viruses and their bee hosts will lead to improved understanding of viral ecology and management strategies that support better bee health.

Temporal changes in the viromes of Swedish Varroa-resistant and Varroa-susceptible honeybee populations

The parasitic mite, Varroa destructor, in combination with the viruses it vectors, is the main cause for global colony losses of the European honeybee, Apis mellifera. However, an isolated honeybee population established in 1999 on the Island of Gotland, Sweden has naturally acquired resistance to the mite, and has survived without mite control treatment for more than 18 years. A recent study has shown that this mite resistant (MR) population also appears to be resistant to Black queen cell virus (BQCV) and Sacbrood virus (SBV) and tolerant to Deformed wing virus (DWV), relative to nearby mite susceptible (MS) honeybee populations. In this study, RNA sequencing was employed to corroborate these previous findings and identify other viral factors that may play a role in the enhanced survival of this mite resistant honeybee population. Two additional honeybee-infecting viruses, Apis rhabdovirus-1 (ARV-1) and Lake Sinai virus (LSV), were identified and near-complete genomes of these two viruses were obtained. Phylogenetic analyses of the assembled virus sequences revealed consistent separation between the MR and MS honeybee populations, although it is unclear whether this is due to pre-existing differences between the viruses in the two populations when they were established, and isolated, or due to virus genetic adaptation towards reduced virulence in the MR population, to promote colony survival. Reverse transcription quantitative polymerase chain reaction(RT-qPCR) analyses show higher ARV and LSV titres in MS colonies compared to MR colonies, gradually increasing from summer to autumn 2009, and reaching maximum titres in the following spring 2010. While the DWV and BQCV titres in MR colonies increased between autumn 2009 and spring 2010, the SBV practically disappeared entirely by spring 2010. Possible explanations for the apparent virus tolerance or resistance in the Gotland mite-resistant honeybee population are discussed.

Extracts of Polypore Mushroom Mycelia Reduce Viruses in Honey Bees

Waves of highly infectious viruses sweeping through global honey bee populations have contributed to recent declines in honey bee health. Bees have been observed foraging on mushroom mycelium, suggesting that they may be deriving medicinal or nutritional value from fungi. Fungi are known to produce a wide array of chemicals with antimicrobial activity, including compounds active against bacteria, other fungi, or viruses. We tested extracts from the mycelium of multiple polypore fungal species known to have antiviral properties. Extracts from amadou (Fomes) and reishi (Ganoderma) fungi reduced the levels of honey bee deformed wing virus (DWV) and Lake Sinai virus (LSV) in a dose-dependent manner. In field trials, colonies fed Ganoderma resinaceum extract exhibited a 79-fold reduction in DWV and a 45,000-fold reduction in LSV compared to control colonies. These findings indicate honey bees may gain health benefits from fungi and their antimicrobial compounds.

Honey Bee and Bumble Bee Antiviral Defense

Bees are important plant pollinators in both natural and agricultural ecosystems. Managed and wild bees have experienced high average annual colony losses, population declines, and local extinctions in many geographic regions. Multiple factors, including virus infections, impact bee health and longevity. The majority of bee-infecting viruses are positive-sense single-stranded RNA viruses. Bee-infecting viruses often cause asymptomatic infections but may also cause paralysis, deformity or death. The severity of infection is governed by bee host immune responses and influenced by additional biotic and abiotic factors. Herein, we highlight studies that have contributed to the current understanding of antiviral defense in bees, including the Western honey bee (Apis mellifera), the Eastern honey bee (Apis cerana) and bumble bee species (Bombus spp.). Bee antiviral defense mechanisms include RNA interference (RNAi), endocytosis, melanization, encapsulation, autophagy and conserved immune pathways including Jak/STAT (Janus kinase/signal transducer and activator of transcription), JNK (c-Jun N-terminal kinase), MAPK (mitogen-activated protein kinases) and the NF-κB mediated Toll and Imd (immune deficiency) pathways. Studies in Dipteran insects, including the model organism Drosophila melanogaster and pathogen-transmitting mosquitos, provide the framework for understanding bee antiviral defense. However, there are notable differences such as the more prominent role of a non-sequence specific, dsRNA-triggered, virus limiting response in honey bees and bumble bees. This virus-limiting response in bees is akin to pathways in a range of organisms including other invertebrates (i.e., oysters, shrimp and sand flies), as well as the mammalian interferon response. Current and future research aimed at elucidating bee antiviral defense mechanisms may lead to development of strategies that mitigate bee losses, while expanding our understanding of insect antiviral defense and the potential evolutionary relationship between sociality and immune function.

Metagenomic analysis of Varroa-free Australian honey bees (Apis mellifera) shows a diverse Picornavirales virome

The viral landscape of the honey bee (Apismellifera) has changed as a consequence of the global spread of the parasitic mite Varroa destructor and accompanying virulent strains of the iflavirus deformed wing virus (DWV), which the mite vectors. The presence of DWV in honey bee populations is known to influence the occurrence of other viruses, suggesting that the current known virome of A. mellifera may be undercharacterized. Here we tested this hypothesis by examining the honey bee virome in Australia, which is uniquely free of parasitic mites or DWV. Using a high-throughput sequencing (HTS) approach, we examined the RNA virome from nine pools of A. mellifera across Australia. In addition to previously reported honey bee viruses, several other insect viruses were detected, including strains related to aphid lethal paralysis virus (ALPV) and Rhopalosiphum padi virus (RhPV), which have recently been identified as infecting honey bees in the USA, as well as several other viruses recently found in Drosophila spp. A further 42 putative novel insect virus genomes spanning the order Picornavirales were assembled, which significantly increases the known viral diversity in A. mellifera. Among these novel genomes, we identified several that were similar (but different) to key A. mellifera viruses, such as DWV, that warrant further investigation. We propose that A. mellifera may be preferentially infected with viruses of the order Picornavirales and that a diverse population of these viruses may be representative of a Varroa-free landscape.