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.

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.

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.

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.

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.

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.

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.

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.

A library preparation optimized for metagenomics of RNA viruses

Our understanding of the viral communities associated to animals has not yet reached the level attained on the bacteriome. This situation is due to, among others, technical challenges in adapting metagenomics using high-throughput sequencing to the study of RNA viromes in animals. Although important developments have been achieved in most steps of viral metagenomics, there is yet a key step that has received little attention: the library preparation. This situation differs from bacteriome studies in which developments in library preparation have largely contributed to the democratisation of metagenomics. Here, we present a library preparation optimized for metagenomics of RNA viruses from insect vectors of viral diseases. The library design allows a simple PCR-based preparation, such as those routinely used in bacterial metabarcoding, that is adapted to shotgun sequencing as required in viral metagenomics. We first optimized our library preparation using mock viral communities and then validated a full metagenomic approach incorporating our preparation in two pilot studies with field-caught insect vectors; one including a comparison with a published metagenomic protocol. Our approach provided a fold increase in virus-like sequences compared to other studies, and nearly-full genomes from new virus species. Moreover, our results suggested conserved trends in virome composition within a population of a mosquito species. Finally, the sensitivity of our approach was compared to a commercial diagnostic PCR for the detection of an arbovirus in field-caught insect vectors. Our approach could facilitate studies on viral communities from animals and the democratization of metagenomics in community ecology of viruses.

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.

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.

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.

Virome Analysis of Aphid Populations That Infest the Barley Field: The Discovery of Two Novel Groups of Nege/Kita-Like Viruses and Other Novel RNA Viruses

Aphids (order Hemiptera) are important insect pests of crops and are also vectors of many plant viruses. However, little is known about aphid-infecting viruses, particularly their diversity and relationship to plant viruses. To investigate the aphid viromes, we performed deep sequencing analyses of the aphid transcriptomes from infested barley plants in a field in Japan. We discovered virus-like sequences related to nege/kita-, flavi-, tombus-, phenui-, mononega-, narna-, chryso-, partiti-, and luteoviruses. Using RT-PCR and sequence analyses, we determined almost complete sequences of seven nege/kitavirus-like virus genomes; one of which was a variant of the Wuhan house centipede virus (WHCV-1). The other six seem to belong to four novel viruses distantly related to Wuhan insect virus 9 (WhIV-9) or Hubei nege-like virus 4 (HVLV-4). We designated the four viruses as barley aphid RNA virus 1 to 4 (BARV-1 to -4). Moreover, some nege/kitavirus-like sequences were found by searches on the transcriptome shotgun assembly (TSA) libraries of arthropods and plants. Phylogenetic analyses showed that BARV-1 forms a clade with WHCV-1 and HVLV-4, whereas BARV-2 to -4 clustered with WhIV-9 and an aphid virus, Aphis glycines virus 3. Both virus groups (tentatively designated as Centivirus and Aphiglyvirus, respectively), together with arthropod virus-like TSAs, fill the phylogenetic gaps between the negeviruses and kitaviruses lineages. We also characterized the flavi/jingmen-like and tombus-like virus sequences as well as other RNA viruses, including six putative novel viruses, designated as barley aphid RNA viruses 5 to 10. Interestingly, we also discovered that some aphid-associated viruses, including nege/kita-like viruses, were present in different aphid species, raising a speculation that these viruses might be distributed across different aphid species with plants being the reservoirs. This study provides novel information on the diversity and spread of nege/kitavirus-related viruses and other RNA viruses that are associated with aphids.

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.

The detection of honey bee (Apis mellifera)-associated viruses in ants

Interspecies virus transmission involving economically important pollinators, including honey bees (Apis mellifera), has recently sparked research interests regarding pollinator health. Given that ants are common pests within apiaries in the southern U.S., the goals of this study were to (1) survey ants found within or near managed honey bee colonies, (2) document what interactions are occurring between ant pests and managed honey bees, and 3) determine if any of six commonly occurring honey bee-associated viruses were present in ants collected from within or far from apiaries. Ants belonging to 14 genera were observed interacting with managed honey bee colonies in multiple ways, most commonly by robbing sugar resources from within hives. We detected at least one virus in 89% of the ant samples collected from apiary sites (n = 57) and in 15% of ant samples collected at non-apiary sites (n = 20). We found that none of these ant samples tested positive for the replication of Deformed wing virus, Black queen cell virus, or Israeli acute paralysis virus, however. Future studies looking at possible virus transmission between ants and bees could determine whether ants can be considered mechanical vectors of honey bee-associated viruses, making them a potential threat to pollinator health.

Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids

Since the discovery that honey bee viruses play a role in colony decline, researchers have made major breakthroughs in understanding viral pathology and infection processes in honey bees. Work on virus transmission patterns and virus vectors, such as the mite Varroa destructor, has prompted intense efforts to manage honey bee health. However, little is known about the occurrence of honey bee viruses in bee predators, such as vespids. In this study, we characterized the occurrence of 11 honey bee viruses in five vespid species and one wasp from four provinces in China and two vespid species from four locations in France. The results showed that all the species from China carried certain honey bee viruses, notably Apis mellifera filamentous virus (AmFV), Deformed wing virus (DWV), and Israeli acute paralysis virus (IAPV); furthermore, in some vespid colonies, more than three different viruses were identified. In France, DWV was the most common virus; Sacbrood virus (SBV) and Black queen cell virus (BQCV) were observed in one and two samples, respectively. Phylogenetic analyses of IAPV and BQCV sequences indicated that most of the IAPV sequences belonged to a single group, while the BQCV sequences belonged to several groups. Additionally, our study is the first to detect Lake Sinai virus (LSV) in a hornet from China. Our findings can guide further research into the origin and transmission of honey bee viruses in Vespidae, a taxon of ecological, and potentially epidemiological, relevance.

Viruses in the Invasive Hornet Vespa velutina

The Asian yellow-legged hornet Vespa velutina nigrithorax, a major predator of honeybees, is spreading in Europe in part due to a lack of efficient control methods. In this study, as a first step to identify biological control agents, we characterized viral RNA sequences present in asymptomatic or symptomatic hornets. Among 19 detected viruses, the honey bee virus Deformed wing virus-B was predominant in all the samples, particularly in muscles from the symptomatic hornet, suggesting a putative cause of the deformed wing symptom. Interestingly, two new viruses closely related to Acyrthosiphon pisumvirus and Himetobi Pvirus and viruses typically associated with honey bees, Acute bee paralysis virus and Black queen cell virus, were detected in the brain and muscles, and may correspond to the circulation and possible replication forms of these viruses in the hornet. Aphid lethal paralysis virus, Bee Macula-like virus, and Moku virus, which are known to infect honey bees, were also identified in the gut virus metagenome of hornets. Therefore, our study underlined the urgent need to study the host range of these newly discovered viruses in hornets to determine whether they represent a new threat for honey bees or a hope for the biocontrol of V. velutina.

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.

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.

Emerging Viruses in Bees: From Molecules to Ecology

Emerging infectious diseases arise as a result of novel interactions between populations of hosts and pathogens, and can threaten the health and wellbeing of the entire spectrum of biodiversity. Bees and their viruses are a case in point. However, detailed knowledge of the ecological factors and evolutionary forces that drive disease emergence in bees and other host-pathogen communities is surprisingly lacking. In this review, we build on the fundamental insight that viruses evolve and adapt over timescales that overlap with host ecology. At the same time, we integrate the role of host community ecology, including community structure and composition, biodiversity loss, and human-driven disturbance, all of which represent significant factors in bee virus ecology. Both of these evolutionary and ecological perspectives represent major advances but, in most cases, it remains unclear how evolutionary forces actually operate across different biological scales (e.g., from cell to ecosystem). We present a molecule-to-ecology framework to help address these issues, emphasizing the role of molecular mechanisms as key bottom-up drivers of change at higher ecological scales. We consider the bee-virus system to be an ideal one in which to apply this framework. Unlike many other animal models, bees constitute a well characterized and accessible multispecies assemblage, whose populations and interspecific interactions can be experimentally manipulated and monitored in high resolution across space and time to provide robust tests of prevailing theory.

Recently identified bee viruses and their impact on bee pollinators

Bees are agriculturally and ecologically important plant pollinators. Recent high annual losses of honey bee colonies, and reduced populations of native and wild bees in some geographic locations, may impact the availability of affordable food crops and the diversity and abundance of native and wild plant species. Multiple factors including viral infections affect pollinator health. The majority of well-characterized bee viruses are picorna-like RNA viruses, which may be maintained as covert infections or cause symptomatic infections or death. Next generation sequencing technologies have been utilized to identify additional bee-infecting viruses including the Lake Sinai viruses and Rhabdoviruses. In addition, sequence data is instrumental for defining specific viral strains and characterizing associated pathogenicity, such as the recent characterization of Deformed wing virus master variants (DWV-A, DWV-B, and DWV-C) and their impact on bee health.