Complete sequence of a picorna-like virus of the genus Iflavirus replicating in the mite Varroa destructor

Entomophthovirus: an insect-derived iflavirus that infects a behavior-manipulating fungal pathogen of dipterans

We report a virus infecting Entomophthora muscae, a behavior-manipulating fungal pathogen of dipterans. The virus, which we name Berkeley Entomophthovirus, is a positive-strand RNA virus in the iflaviridae family of capsid-forming viruses, which are mostly known to infect insects. The viral RNA is expressed at high levels in fungal cells in vitro and during in vivo infections of Drosophila melanogaster, and virus particles can be seen intracellularly in E. muscae. This virus, of which we find two closely related variants in our culture of E. muscae, is also closely related to three different viruses reported from metagenomic surveys, two of which were isolated from wild dipterans, and a third isolated from wild ticks. By analyzing sequencing data from these earlier reports, we find abundant reads aligning to E. muscae specifically in the samples from which viral reads were sequenced. These data establish a wide and perhaps obligate association with E. muscae in the wild, consistent with our laboratory data that E. muscae is the host for these closely related viruses. Because of this, we propose the name Entomophthovirus (EV) for this group of highly related virus variants. As other members of the iflaviridae have been reported to cause behavioral changes in insects, we speculate on the possibility that EV plays a role in the behavioral manipulation of flies infected with E. muscae.

Molecular Characterization and Phylogenetic Analysis of Honeybee (Apis mellifera) Mite-Borne Pathogen DWV-A and DWV-B Isolated from Lithuania

Deformed wing virus (DWV) is known as one of the main viruses that affect honeybees' health all around the world. The virus has two widespread genotypes, DWV-A and DWV-B (VDV-1), transmitted mainly by V. destructor mites. In this study, we collected honeycombs with covered broods from 73 apiaries in eight Lithuanian regions and initially investigated the prevalence of V. destructor mites. Mites were collected from May to the end of July in 2021 from 124 hives. The prevalence of V. destructor infestations in beehives reached 30% and 63% in investigated apiaries. The presence of DWV-A and DWV-B pathogens in mites and broods was examined by RT-qPCR targeting the CRPV-capsid region. The molecular characterization of the virus in mite samples was based on sequence analysis of the RNA-dependent RNA polymerase (RdRp) region. In addition, leader polypeptide (LP), structural protein (Vp3), Helicase, and RdRp genes were used for phylogenetic characterization of dual infection. The prevalences of DWV-B in mites and broods were 56.5% and 31.5%, respectively, while DWV-A was detected in 12.9% of mite samples and 24.7% of brood samples. Some of the examined mite samples harboured dual virus infections. Our findings showed that bee colonies from the same apiary were not always infected by the same viruses. Some bee colonies were virus-free, while others were highly infected. Phylogenetic analysis of 21 sequences demonstrated the presence of highly variable DWV-B and DWV-A genotypes in Lithuania and possible recombinant variants of the virus. This study represents the first molecular characterization of mite-borne pathogens hosted by honeybees (Apis mellifera) in Lithuania.

Three picorna-like viruses found associated with the spider mite, Tetranychus truncatus (Acari: Tetranychidae)

Herbivorous arthropods, such as mites and insects, host a variety of microorganisms that significantly influence their ecology and evolution. While insect viruses have been extensively studied, our understanding of the diversity and composition of mite viromes and the interactions with mite hosts remains limited. The Asian spider mite, Tetranychus truncatus Ehara (Acari: Tetranychidae), a major agricultural pest, has not yet been reported to harbor any viruses. Here, using publicly available RNA-Seq data, we identified and characterized three picorna-like viruses associated with T. truncatus: Tetranychus truncatus-associated iflavirus 1 (TtAIV-1), Tetranychus truncatus-associated picorna-like virus 1 (TtAV-1), and Tetranychus truncatus-associated picorna-like virus 2 (TtAV-2). TtAIV-1 has a typical Iflaviridae genome structure with a single ORF, representing the first iflavirus associated with the Tetranychus genus. TtAV-1 and TtAV-2 exhibit bicistronic arrangements similar to dicistroviruses and other picorna-like viruses, with complex secondary structures in their non-coding regions. Phylogenetic analysis places TtAIV-1 within Iflaviridae, possibly as a new species, while TtAV-1 and TtAV-2 form distinct clades within unclassified picorna-like viruses, suggesting new families within Picornavirales. We analyzed in silico the presence and abundance of these viruses in T. truncatus across four bioproject SRAs, mostly finding them co-associated, with viral reads reaching up to 30% of total reads. Their presence and abundance varied by mite treatment and origin, with no significant impact from Wolbachia infection or abamectin exposure, although TtAV-2 was absent in abamectin-treated mites. Temperature influenced virus abundance, and variations were observed among Chinese mite populations based on geography and host plant association. Our findings offer insights into picorna-like virus diversity and dynamics in T. truncatus, revealing potential roles in mite biology and suggesting applications for mite population control, thereby enhancing agricultural productivity and food security.

Screening and Characterization of a New Iflavirus Virus in the Fruit Tree Pest Pyrops candelaria

Pyrops candelaria is one of the common pests of fruit trees, but the research on the pathogenic microorganisms it may carry is very limited. Therefore, it is essential to reveal the pathogenic microbes it carries and their potential hazards. This study found a new virus from the transcriptome of P. candelaria, which was first reported in P. candelaria and named PyCaV (Pyrops candelaria associated virus). RACE and bioinformatics assay revealed that the full length of PyCaV is 10,855 bp with the polyA tail, containing a single open-reading frame (ORF) encoding a polyprotein consisting of 3171 amino acid (aa). The virus has a typical iflavirus structure, including two rhv domains, an RNA helicase domain (HEL), a 3C cysteine protease domain (Pro), and an RNA-dependent RNA polymerase domain (RdRp). Further phylogenetic analysis revealed that this virus belongs to family Iflaviridae and sequence alignments analysis suggested PyCaV is a new member in an unassigned genus of family Iflaviridae. Further in-depth analysis of the virus infection showed that PyCaV is distributed throughout the whole P. candelaria, including its head, chest, and abdomen, but more PyCaV was identified in the chest. The distribution of PyCaV in different parts of P. candelaria was further explored, which showed that more PyCaV was detected in its piercing-sucking mouthparts and chest viscera. Statistical analysis showed that the PyCaV infection was affected by time and location.

Molecular Detection and Phylogenetic Relationships of Honey Bee-Associated Viruses in Bee Products

In the last few years, the isolation and amplification of DNA or RNA from the environment (eDNA/eRNA) has proven to be an alternative and non-invasive approach for molecular identification of pathogens and pests in beekeeping. We have recently demonstrated that bee pollen and bee bread represent suitable biological material for the molecular identification of viral RNA. In the present study, we extracted total RNA from different bee products (pollen, n = 25; bee bread, n = 17; and royal jelly, n = 15). All the samples were tested for the presence of six of the most common honey bee-associated viruses-Deformed wing virus (DWV), Acute bee paralysis virus (ABPV), Chronic bee paralysis virus (CBPV), Sacbrood virus (SBV), Kashmir bee virus (KBV), and Black queen cell virus (BQCV)-using a reverse transcription polymerase chain reaction (RT-PCR). We successfully detected six records of DWV (10.5%, 6/57), four of ABPV (7.0%, 4/57), three of Israeli acute paralysis virus (IAPV) (5.3%, 3/57), and two of BQCV (3.5%, 2/57). Using ABPV primers, we also successfully detected the presence of IAPV. The obtained viral sequences were analyzed for phylogenetic relationships with the highly similar sequences (megablast) available in the GenBank database. The Bulgarian DWV isolates revealed a high homology level with strains from Syria and Turkey. Moreover, we successfully detected a DWV strain B for the first time in Bulgaria. In contrast to DWV, the ABPV isolates formed a separate clade in the phylogenetic tree. BQCV was closely grouped with Russian isolates, while Bulgarian IAPV formed its own clade and included a strain from China. In conclusion, the present study demonstrated that eRNA can be successfully used for molecular detection of honey bee-associated viruses in bee products. The method can assist the monitoring of the health status of honey bee colonies at the local, regional, and even national levels.

Deformed wing virus genotypes A and B do not elicit immunologically different responses in naïve honey bee hosts

Iflavirus aladeformis (Picornavirales: Iflaviridae), commonly known as deformed wing virus(DWV), in association with Varroa destructor Anderson and Trueman (Mesostigmata: Varroidae), is a leading factor associated with honey bee (Apis mellifera L. [Hymenoptera: Apidae]) deaths. The virus and mite have a near global distribution, making it difficult to separate the effect of one from the other. The prevalence of two main DWV genotypes (DWV-A and DWV-B) has changed over time, leading to the possibility that the two strains elicit a different immune response by the host. Here, we use a honey bee population naïve to both the mite and the virus to investigate if honey bees show a different immunological response to DWV genotypes. We examined the expression of 19 immune genes by reverse transcription quantitative PCR (RT-qPCR) and analysed small RNA after experimental injection with DWV-A and DWV-B. We found no evidence that DWV-A and DWV-B elicit different immune responses in honey bees. RNA interference genes were up-regulated during DWV infection, and small interfering RNA (siRNA) responses were proportional to viral loads yet did not inhibit DWV accumulation. The siRNA response towards DWV was weaker than the response to another honey bee pathogen, Triatovirus nigereginacellulae (Picornavirales: Dicistroviridae; black queen cell virus), suggesting that DWV is comparatively better at evading host antiviral defences. There was no evidence for the production of virus-derived Piwi-interacting RNAs (piRNAs) in response to DWV. In contrast to previous studies, and in the absence of V. destructor, we found no evidence that DWV has an immunosuppressive effect. Overall, our results advance our understanding of the immunological effect that DWV in isolation elicits in honey bees.

Origins, diversity, and adaptive evolution of DWV in the honey bees of the Azores: the impact of the invasive mite Varroa destructor

Deformed wing virus (DWV) is a honey bee virus, whose emergence from relative obscurity is driven by the recent host-switch, adaptation, and global dispersal of the ectoparasitic mite Varroa destructor (a highly efficient vector of DWV) to reproduction on honey bees (Apis mellifera). Our study examines how varroa affects the continuing evolution of DWV, using the Azores archipelago, where varroa is present on only three out of the eight Islands, as a natural experimental system for comparing different evolutionary conditions and trajectories. We combined qPCR of 494 honey bee colonies sampled across the archipelago with amplicon deep sequencing to reveal how the DWV genetic landscape is altered by varroa. Two of the varroa-free Islands were also free of DWV, while a further two Islands were intriguingly dominated by the rare DWV-C major variant. The other four Islands, including the three varroa-infested Islands, were dominated by the common DWV-A and DWV-B variants. The varroa-infested Islands had, as expected, an elevated DWV prevalence relative to the uninfested Islands, but not elevated DWV loads, due the relatively high prevalence and loads of DWV-C on the varroa-free Islands. This establishes the Azores as a stable refuge for DWV-C and provides the most convincing evidence to date that at least some major strains of DWV may be capable of not just surviving, but actually thriving in honey bees in the absence of varroa-mediated transmission. We did not detect any change in DWV genetic diversity associated with island varroa status but did find a positive association of DWV diversity with virus load, irrespective of island varroa status.

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.

Characterization of a Molecular Clone of Deformed Wing Virus B

Honey bees (Apis mellifera) play a crucial role in agriculture through their pollination activities. However, they have faced significant health challenges over the past decades that can limit colony performance and even lead to collapse. A primary culprit is the parasitic mite Varroa destructor, known for transmitting harmful bee viruses. Among these viruses is deformed wing virus (DWV), which impacts bee pupae during their development, resulting in either pupal demise or in the emergence of crippled adult bees. In this study, we focused on DWV master variant B. DWV-B prevalence has risen sharply in recent decades and appears to be outcompeting variant A of DWV. We generated a molecular clone of a typical DWV-B strain to compare it with our established DWV-A clone, examining RNA replication, protein expression, and virulence. Initially, we analyzed the genome using RACE-PCR and RT-PCR techniques. Subsequently, we conducted full-genome RT-PCR and inserted the complete viral cDNA into a bacterial plasmid backbone. Phylogenetic comparisons with available full-length sequences were performed, followed by functional analyses using a live bee pupae model. Upon the transfection of in vitro-transcribed RNA, bee pupae exhibited symptoms of DWV infection, with detectable viral protein expression and stable RNA replication observed in subsequent virus passages. The DWV-B clone displayed a lower virulence compared to the DWV-A clone after the transfection of synthetic RNA, as evidenced by a reduced pupal mortality rate of only 20% compared to 80% in the case of DWV-A and a lack of malformations in 50% of the emerging bees. Comparable results were observed in experiments with low infection doses of the passaged virus clones. In these tests, 90% of bees infected with DWV-B showed no clinical symptoms, while 100% of pupae infected with DWV-A died. However, at high infection doses, both DWV-A and DWV-B caused mortality rates exceeding 90%. Taken together, we have generated an authentic virus clone of DWV-B and characterized it in animal experiments.

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.

Virus replication in the honey bee parasite, Varroa destructor

IMPORTANCE

The parasitic mite Varroa destructor is a significant driver of worldwide colony losses of our most important commercial pollinator, the Western honey bee Apis mellifera. Declines in honey bee health are frequently attributed to the viruses that mites vector to honey bees, yet whether mites passively transmit viruses as a mechanical vector or actively participate in viral amplification and facilitate replication of honey bee viruses is debated. Our work investigating the antiviral RNA interference response in V. destructor demonstrates that key viruses associated with honey bee declines actively replicate in mites, indicating that they are biological vectors, and the host range of bee-associated viruses extends to their parasites, which could impact virus evolution, pathogenicity, and spread.

Processing of the 3C/D Region of the Deformed Wing Virus (DWV)

The deformed wing virus (DWV) belongs to the genus Iflavirus and the family Iflaviridae within the order Picornavirales. It is an important pathogen of the Western honey bee, Apis mellifera, causing major losses among honey bee colonies in association with the ectoparasitic mite Varroa destructor. Although DWV is one of the best-studied insect viruses, the mechanisms of viral replication and polyprotein processing have been poorly studied in the past. We investigated the processing of the protease-polymerase region at the C-terminus of the polyprotein in more detail using recombinant expression, novel serological reagents, and virus clone mutagenesis. Edman degradation of purified maturated polypeptides uncovered the C- and N-termini of the mature 3C-like (3CL) protease and RNA-dependent RNA polymerase (3DL, RdRp), respectively. Autocatalytic processing of the recombinant DWV 3CL protease occurred at P1 Q2118 and P1' G2119 (KPQ/GST) as well as P1 Q2393 and P1' S2394 (HAQ/SPS) cleavage sites. New monoclonal antibodies (Mab) detected the mature 3CL protease with an apparent molecular mass of 32 kDa, mature 3DL with an apparent molecular mass of 55 kDa as well as a dominant 3CDL precursor of 90 kDa in DWV infected honey bee pupae. The observed pattern corresponds well to data obtained via recombinant expression and N-terminal sequencing. Finally, we were able to show that 3CL protease activity and availability of the specific protease cleavage sites are essential for viral replication, protein synthesis, and establishment of infection using our molecular clone of DWV-A.

Establishment and Application of CRISPR-Cas12a-Based Recombinase Polymerase Amplification and a Lateral Flow Dipstick and Fluorescence for the Detection and Distinction of Deformed Wing Virus Types A and B

Deformed wing virus (DWV) is one of the important pathogens of the honey bee (Apis mellifera), which consists of three master variants: types A, B, and C. Among them, DWV types A (DWV-A) and B (DWV-B) are the most prevalent variants in honey bee colonies and have been linked to colony decline. DWV-A and DWV-B have different virulence, but it is difficult to distinguish them via traditional methods. In this study, we established a visual detection assay for DWV-A and DWV-B using recombinase polymerase amplification (RPA) and a lateral flow dipstick (LFD) coupled with the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein (Cas) 12a fluorescence system (RPA-CRISPR-Cas12a-LFD). The limit of detection of this system was ~6.5 × 100 and 6.2 × 101 copies/μL for DWV-A and DWV-B, respectively. The assays were specific and non-cross-reactive against other bee viruses, and the results could be visualized within 1 h. The assays were validated by extracting cDNA from 36 clinical samples of bees that were suspected to be infected with DWV. The findings were consistent with those of traditional reverse transcription-quantitative polymerase chain reaction, and the RPA-CRISPR-Cas12a assay showed the specific, sensitive, simple, and appropriate detection of DWV-A and DWV-B. This method can facilitate the visual and qualitative detection of DWV-A and DWV-B as well as the monitoring of different subtypes, thereby providing potentially better control and preventing current and future DWV outbreaks.

The presence of identical deformed wing virus sequence variants in co-occurring Apis species in Northern Thailand may represent a potential epidemiological threat to native honey bees of Southeast Asia

Widespread native honey bee species in South and East Asia (Apis cerana, Apis dorsata and Apis florea) and the imported western honey bee (Apis mellifera) share habitats and potentially also share pathogens. Chief among the threats facing A. mellifera in Europe and North America is deformed wing virus (DWV), including its two principal genotypes: A and B (DWV-A and DWV-B respectively). Though DWV-A has been recorded in Asia's native Apis species, it is not known if DWV-B, or both DWV-A and DWV-B, are currently widespread in Asia and, if so, whether viral transmission is primarily intraspecific or interspecific. This study aims to fill these knowledge gaps by (i) determining the DWV genotype in four co-occurring Apis host species using qPCR and (ii) inferring viral transmission between them using nucleotide sequences of DWV from Apis host species collected at three independent localities in Northern Thailand. We found DWV-A and -B in all four Apis species, the exotic A. mellifera and the native A. cerana, A. dorsata and A. florea. That DWV-A sequences were identical across Apis species at the same locality, with a similar pattern for DWV-B sequences, suggests that DWV's epidemiology is largely driven by ongoing interspecific transmission (spillover) of DWV across co-occurring native and exotic Apis species. Both genotypes of DWV represent a serious threat to Asia's exotic and native honey bee species.

Viral Co-Infections and Antiviral Immunity in Honey Bees

Over the past few decades, honey bees have been facing an increasing number of stressors. Beyond individual stress factors, the synergies between them have been identified as a key factor in the observed increase in colony mortality. However, these interactions are numerous and complex and call for further research. Here, in line with our need for a systemic understanding of the threats that they pose to bee health, we review the interactions between honey bee viruses. As viruses are obligate parasites, the interactions between them not only depend on the viruses themselves but also on the immune responses of honey bees. Thus, we first summarise our current knowledge of the antiviral immunity of honey bees. We then review the interactions between specific pathogenic viruses and their interactions with their host. Finally, we draw hypotheses from the current literature and suggest directions for future research.

Significant, but not biologically relevant: Nosema ceranae infections and winter losses of honey bee colonies

The Western honey bee Apis mellifera, which provides about 90% of commercial pollination, is under threat from diverse abiotic and biotic factors. The ectoparasitic mite Varroa destructor vectoring deformed wing virus (DWV) has been identified as the main biotic contributor to honey bee colony losses worldwide, while the role of the microsporidium Nosema ceranae is still controversially discussed. In an attempt to solve this controversy, we statistically analyzed a unique data set on honey bee colony health collected from a cohort of honey bee colonies over 15 years and comprising more than 3000 data sets on mite infestation levels, Nosema spp. infections, and winter losses. Multivariate statistical analysis confirms that V. destructor is the major cause of colony winter losses. Although N. ceranae infections are also statistically significantly correlated with colony losses, determination of the effect size reveals that N. ceranae infections are of no or low biological relevance.

Characterization of Iflavirus in the Red Flour Beetle, Tribolium castaneum (Coleoptera; Tenebrionidae)

Iflavirus is a group of viruses distributed mainly in arthropod species. We surveyed Tribolium castaneum iflavirus (TcIV) in different laboratory strains and in Sequence Read Archives (SRA) in GenBank. TcIV is highly specific to only T. castaneum and is not found in seven other Tenebrionid species, including the closely related species T. freemani. The same strains from different laboratories and different strains displayed largely different degrees of infections in the examination of 50 different lines by using Taqman-based quantitative PCR. We found that ~63% (27 out of 43 strains) of T. castaneum strains in different laboratories are positive for TcIV PCR with large degrees of variation, in the range of seven orders of magnitude, indicating that the TcIV is highly fluctuating depending on the rearing conditions. The TcIV was prevalent in the nervous system with low levels found in the gonad and gut. The transovarial transmission was supported in the experiment with surface-sterilized eggs. Interestingly, TcIV infection did not show observable pathogenicity. TcIV offers an opportunity to study the interaction between the virus and the immune system of this model beetle species.

Nationwide Screening for Bee Viruses in Apis mellifera Colonies in Egypt

Honey bees are essential for crop and wild plant pollination. However, many countries have reported high annual colony losses caused by multiple possible stressors. Diseases, particularly those caused by viruses, are a major cause of colony losses. However, little is known about the prevalence of honey bee pathogens, particularly virus prevalence, in Egyptian honey bees. To address this shortfall, we determined the prevalence of widespread bee viruses in honey bee colonies in Egypt-whether it is affected by geography, the season, or infestation with Varroa destructor (varroa) mites. Honey bee worker samples were collected from 18 geographical regions across Egypt during two seasons: winter and summer of 2021. Three apiaries were chosen in each region, and a pooled sample of 150 worker bees was collected from five colonies in each apiary then screened by qPCR for 10 viral targets: acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), deformed wing virus (DWV) genotypes A (DWV-A), B (DWV-B) and D (Egyptian bee virus), Israeli acute paralysis virus (IAPV), Kashmir bee virus (KBV), sacbrood virus (SBV), and slow bee paralysis virus (SBPV). Our results revealed that DWV-A was the most prevalent virus, followed by BQCV and ABPV; the DWV genotype now spreading across the world, DWV-B, was not detected. There was no difference in varroa infestation rates as well as virus prevalence between winter and summer. However, colonies infected with BQCV had a significantly higher varroa count (adjusted p < 0.05) in the winter season, indicating that there is a seasonal association between the intensity of infestation by varroa and the presence of this virus. We provide data on the current virus prevalence in Egypt, which could assist in the protection of Egypt's beekeeping industry. Moreover, our study aids in the systematic assessment of the global honey bee virome by filling a knowledge gap about the prevalence of honey bee viruses in Egypt.

The vectoring competence of the mite Varroa destructor for deformed wing virus of honey bees is dynamic and affects survival of the mite

The ectoparasitic mite, Varroa destructor and the viruses it vectors, including types A and B of Deformed wing virus (DWV), pose a major threat to honey bees, Apis mellifera. Analysis of 256 mites collected from the same set of field colonies on five occasions from May to October 2021 showed that less than a half of them, 39.8% (95% confidence interval (CI): 34.0 - 46.0%), were able to induce a high (overt) level DWV infection with more than 109 viral genomes per bee in the pupa after 6 days of feeding, with both DWV-A and DWV-B being vectored at similar rates. To investigate the effect of the phoretic (or dispersal) stage on adult bees on the mites' ability to vector DWV, the mites from two collection events were divided into two groups, one of which was tested immediately for their infectiveness, and the other was kept with adult worker bees in cages for 12 days prior to testing their infectiveness. We found that while 39.2% (95% CI: 30.0 - 49.1%) of the immediately tested mites induced overt-level infections, 12-day passage on adult bees significantly increased the infectiousness to 89.8% (95% CI: 79.2 - 95.6%). It is likely that Varroa mites that survive brood interruptions in field colonies are increasingly infectious. The mite lifespan was affected by the DWV type it transmitted to pupae. The mites, which induced high DWV-B but not DWV-A infection had an average lifespan of 15.5 days (95% CI: 11.8 - 19.2 days), which was significantly shorter than those of the mites which induced high DWV-A but not DWV-B infection, with an average lifespan of 24.3 days (95% CI: 20.2 - 28.5), or the mites which did not induce high levels of DWV-A or DWV-B, with an average survival of 21.2 days (95% CI: 19.0 - 23.5 days). The mites which transmitted high levels of both DWV-A and DWV-B had an intermediate average survival of 20.5 days (95% CI: 15.1 - 25.9 days). The negative impact of DWV-B on mite survival could be a consequence of the ability of DWV-B, but not DWV-A to replicate in Varroa.

Epidemiology of a major honey bee pathogen, deformed wing virus: potential worldwide replacement of genotype A by genotype B

The western honey bee (Apis mellifera) is of major economic and ecological importance, with elevated rates of colony losses in temperate regions over the last two decades thought to be largely caused by the exotic ectoparasitic mite Varroa destructor and deformed wing virus (DWV), which the mite transmits. DWV currently exists as two main genotypes: the formerly widespread DWV-A and the more recently described and rapidly expanding DWV-B. It is an excellent system to understand viral evolution and the replacement of one viral variant by another. Here we synthesise published results on the distribution and prevalence of DWV-A and -B over the period 2008-2021 and present novel data for Germany, Italy and the UK to suggest that (i) DWV-B has rapidly expanded worldwide since its first description in 2004 and (ii) that it is potentially replacing DWV-A. Both genotypes are also found in wild bee species. Based on a simple mathematical model, we suggest that interference between viral genotypes when co-infecting the same host is key to understanding their epidemiology. We finally discuss the consequences of genotype replacement for beekeeping and for wild pollinator species.

Quantitative and Qualitative Changes in the Deformed Wing Virus Population in Honey Bees Associated with the Introduction or Removal of Varroa destructor

Varroa destructor is an ectoparasitic mite associated with significant losses of honeybee colonies globally. The mite vectors a range of pathogenic viruses, the most important of which is the Deformed wing virus (DWV). In the absence of Varroa, DWV exists as a low-level, highly diverse virus population. However, when transmitted by Varroa, certain variants become highly elevated, and may become near-clonal and cause symptomatic infections. Mite transmission between colonies can occur when parasitised workers drift from or rob adjacent hives. These activities can result in elevated mite levels, but the resulting change in the DWV population, the primary determinant of winter colony losses, has not been determined. In reciprocal studies, we investigated the influence of the removal of mites, or their acquisition, on the DWV population. When mites were removed from heavily infested colonies, there was a striking and rapid reduction in virus load. Conversely, siting Varroa-naïve colonies in a mite-infested apiary resulted in the acquisition of mites and concomitant changes in the virus population. We observed both near-clonal and highly divergent virus populations regardless of titre, suggesting changes were stochastic and colony-specific. Our findings have implications for the outcome of strategies in areas with total or patchy implementation of Varroa control plans.

Interaction between Thiamethoxam and Deformed Wing Virus Type A on Wing Characteristics and Expression of Immune and Apoptosis Genes in Apis mellifera

Simple Summary

Honey bees are key pollinators in agricultural crops. Today, honey bee colonies in decline are a global concern as a result of various stressors, including pesticides, pathogens, honey bee health, and parasites. A healthy honey bee colony refers to colonies that are not exposed to biotic and abiotic stressors. In this study, we examine how thiamethoxam (pesticide) and deformed wing virus type A (DWV-A) interact in effects on honey bee health. The results revealed that the honey bees were infected with DWV-A and were additionally exposed to thiamethoxam, showing effects that increased the mortality rate, and crippled wings in newly emerged adult honey bees. Moreover, the exposure to thiamethoxam and DWV-A injection resulted in induced expression of immune genes (hymenoptaecin gene) while downregulation of two apoptosis genes (caspase8-like, caspase9-like genes). The impact interaction of pesticide and DWV-A have on the expression of apoptosis genes can directly affect viral susceptibility in the honey bee host.

Abstract

Honey bees are economically important insects for crop pollination. They play a significant role as pollinators of wild plants and agricultural crops and produce economical products, such as honey, royal jelly, wax, pollen, propolis, and venom. Despite their ecological and economical importance, the global honey bee population is in decline due to factors including pathogens, parasites, intensive agriculture, and pesticides. Moreover, these factors may be interlinked and exacerbate the loss of honey bees. This study aimed to investigate the interaction between a pesticide, thiamethoxam, and deformed wing virus type A (DWV-A) to honey bees and the effects on survival rate, wing characteristics, and expression of immune and apoptosis genes in Apis mellifera. We described the potential interaction between thiamethoxam and DWV-A on honey bee wing characteristics, DWV-A loads, and the expressions of immune (defensin, abaecin, and hymenoptaecin) and apoptosis genes (buffy, apaf1, caspase3-like, caspase8-like, and caspase9-like). Honey bee larvae were fed with three different thiamethoxam doses (0.001, 1.4, and 14.3 ng/µL of the diet). Then, thiamethoxam-treated white-eyed pupae were injected with 10⁷ copy numbers/honey bee of the DWV-A genome. The interaction between thiamethoxam and DWV-A caused a high mortality rate, crippled wings in newly emerged adult honey bees (100%), and resulted in induced expression of hymenoptaecin gene compared to the control group, while downregulation of caspase8-like, caspase9-like genes compared to the DWV injection group. Therefore, the potential interaction between thiamethoxam and DWV-A might have a deleterious effect on honey bee lifespan. The results from this study could be used as a tool to combat DWV-A infection and mitigate pesticide usage to alleviate the decrease in the honey bee population.

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

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

Cold case: The disappearance of Egypt bee virus, a fourth distinct master strain of deformed wing virus linked to honeybee mortality in 1970's Egypt

In 1977, a sample of diseased adult honeybees (Apis mellifera) from Egypt was found to contain large amounts of a previously unknown virus, Egypt bee virus, which was subsequently shown to be serologically related to deformed wing virus (DWV). By sequencing the original isolate, we demonstrate that Egypt bee virus is in fact a fourth unique, major variant of DWV (DWV-D): more closely related to DWV-C than to either DWV-A or DWV-B. DWV-A and DWV-B are the most common DWV variants worldwide due to their close relationship and transmission by Varroa destructor. However, we could not find any trace of DWV-D in several hundred RNA sequencing libraries from a worldwide selection of honeybee, varroa and bumblebee samples. This means that DWV-D has either become extinct, been replaced by other DWV variants better adapted to varroa-mediated transmission, or persists only in a narrow geographic or host range, isolated from common bee and beekeeping trade routes.

Haplotype Analysis of Varroa destructor and Deformed Wing Virus Using Long Reads

As a phoretic parasite and virus vector, the mite Varroa destructor and the associated Deformed wing virus (DWV) form a lethal combination to the honey bee, Apis mellifera. Routine acaricide treatment has been reported to reduce the diversity of mites and select for tolerance against these treatments. Further, different DWV strains face selective pressures when transmitted via mites. In this study, the haplotypes of Varroa mites and associated DWV variants were quantified using long reads. A single haplotype dominated the mite mitochondrial gene cytochrome oxidase subunit I, reflecting an ancient bottleneck. However, highly polymorphic genes were present across the mite genome, suggesting the diversity of mites could be actively maintained at a regional level. DWV detected in both mites and honey bees show a dominant variant with only a few low-frequency alternate haplotypes. The relative abundances of DWV haplotypes isolated from honey bees and mites were highly consistent, suggesting that some variants are favored by ongoing selection.

First come, first served: superinfection exclusion in Deformed wing virus is dependent upon sequence identity and not the order of virus acquisition

Deformed wing virus (DWV) is the most important globally distributed pathogen of honey bees and, when vectored by the ectoparasite Varroa destructor, is associated with high levels of colony losses. Divergent DWV types may differ in their pathogenicity and are reported to exhibit superinfection exclusion upon sequential infections, an inevitability in a Varroa-infested colony. We used a reverse genetic approach to investigate competition and interactions between genetically distinct or related virus strains, analysing viral load over time, tissue distribution with reporter gene-expressing viruses and recombination between virus variants. Transient competition occurred irrespective of the order of virus acquisition, indicating no directionality or dominance. Over longer periods, the ability to compete with a pre-existing infection correlated with the genetic divergence of the inoculae. Genetic recombination was observed throughout the DWV genome with recombinants accounting for ~2% of the population as determined by deep sequencing. We propose that superinfection exclusion, if it occurs at all, is a consequence of a cross-reactive RNAi response to the viruses involved, explaining the lack of dominance of one virus type over another. A better understanding of the consequences of dual- and superinfection will inform development of cross-protective honey bee vaccines and landscape-scale DWV transmission and evolution.

Screening of Honey Bee Pathogens in the Czech Republic and Their Prevalence in Various Habitats

Western honey bee (Apis mellifera) is one of the most important pollinators in the world. Thus, a recent honey bee health decline and frequent honey bee mass losses have drawn attention and concern. Honey bee fitness is primarily reduced by pathogens, parasites, and viral load, exposure to pesticides and their residues, and inadequate nutrition from both the quality and amount of food resources. This study evaluated the prevalence of the most common honey bee pathogens and viruses in different habitats across the Czech Republic. The agroecosystems, urban ecosystems, and national park were chosen for sampling from 250 colonies in 50 apiaries. Surprisingly, the most prevalent honey bee pathogens belong to the family Trypanosomatidae including Lotmaria passim and Crithidia mellificae. As expected, the most prevalent viruses were DWV, followed by ABPV. Additionally, the occurrence of DWV-B and DWV-C were correlated with honey bee colony mortality. From the habitat point of view, most pathogens occurred in the town habitat, less in the agroecosystem and least in the national park. The opposite trend was observed in the occurrence of viruses. However, the prevalence of viruses was not affected by habitat.

Occurrence and Phylogenetic Analysis of DWV in Stingless Bee (Apidae sp.) in China: A Case Report

Honey bees play a vital role in providing pollination services for agricultural crops and wild flowering plants. However, the spillover risk of their pathogens to other pollinators or wild insects is becoming a cause for concern. There is some evidence that stingless bees can carry honey bee viruses, but little is known about the presence of honey bee viruses in stingless bees in China. Here, we investigate the occurrence of major honey bee pathogens including bacteria, fungi, and viruses in stingless bees (Apidae: sp.). Our results show that the stingless bees (Apidae: sp.) were mainly infected with DWV-A, but no DWV-B and DWV-C. Phylogenetic analysis on fragments of lp, RdRp, and VP3 of DWV-A indicated that genetic variation in VP3 might an important indicator for host-specific viruses, but it requires further study. Our results indicated that DWV-A is not only the major strain of virus currently circulating in managed bee colonies in China and globally, but in stingless bee species as a whole.

Simulated vector transmission differentially influences dynamics of two viral variants of deformed wing virus in honey bees (Apis mellifera)

Understanding how vectors alter the interactions between viruses and their hosts is a fundamental question in virology and disease ecology. In honey bees, transmission of deformed wing virus (DWV) by parasitic Varroa mites has been associated with elevated disease and host mortality, and Varroa transmission has been hypothesized to lead to increased viral titres or select for more virulent variants. Here, we mimicked Varroa transmission by serially passaging a mixed population of two DWV variants, A and B, by injection through in vitro reared honey bee pupae and tracking these viral populations through five passages. The DWV-A and DWV-B variant proportions shifted dynamically through passaging, with DWV-B outcompeting DWV-A after one passage, but levels of both variants becoming equivalent by Passage 5. Sequencing analysis revealed a dominant, recombinant DWV-B strain (DWV-A derived 5' IRES region with the rest of the genome DWV-B), with low nucleotide diversity that decreased through passaging. DWV-A populations had higher nucleotide diversity compared to DWV-B, but this also decreased through passaging. Selection signatures were found across functional regions of the DWV-A and DWV-B genomes, including amino acid mutations in the putative capsid protein region. Simulated vector transmission differentially impacted two closely related viral variants which could influence viral interactions with the host, demonstrating surprising plasticity in vector-host-viral dynamics.

Real-time monitoring of deformed wing virus-infected bee foraging behavior following histone deacetylase inhibitor treatment

Impairment in the learning/memory behavior of bees is responsible for the massive disappearance of bee populations and its consequent agricultural economic losses. Such impairment might be because of o both pesticide exposure and pathogen infection, with a key contributor deformed wing virus (DWV). The present study found that sodium butyrate (NaB) significantly increased survival and reversed the learning/memory impairment of DWV-infected bees. A next-generation sequencing analysis showed that NaB affected the expression of genes involved in glycolytic processes and memory formation, which were suppressed by DWV infection. In addition, we performed a large-scale movement tracking experiment by using a wireless sensor network-based automatic real-time monitoring system and confirmed that NaB could improve the homing ability of DWV-infected bees. In short, we demonstrated the mechanism of how epigenetic regulation can resume the memory function of honeybees and suggest strategies for applying NaB to reduce the incidence of colony losses.

Variant A of the Deformed Wings Virus Alters the Olfactory Sensitivity and the Expression of Odorant Binding Proteins on Antennas of Apis mellifera

Simple Summary

Honey bees, Apis melllifera, are the most commonly managed bee in the world for pollination services. However, worldwide continuous colony losses have been reported for almost a decade. One factor of these losses is associated to pathogens being the virus one of the most important problems in honey bee health. One of the known viruses that affect the honey bee population is deformed wing virus (DWV). DWV causes physical malformation and behavioral disturbances, but also, this virus can be found in the antenna affecting the anatomical integrity of infected areas, which could compromise normal antennal functioning associated to aroma perception. Thus, we evaluate olfactory sensitivity and the expression of antenna-specific odorant-binding proteins (OBP) genes in honey bees inoculated with variant A of the DWV. We performed olfactory sensitivity analysis using the essential oils Eucalyptus globulus and Mentha piperita, but also, and molecular analysis of gene expression of nine OBPs. We found that the high level of replication of DWV-A in the antennae decreased the olfactory sensitivity and led to a down-regulation of some OBPs in middle- and forager-age worker bees. Thus, DWV-A infection in adults of honey bees could compromise volatile compound recognition inside the hive and outside the hive.

Abstract

Insects have a highly sensitive sense of smell, allowing them to perform complex behaviors, such as foraging and peer recognition. Their sense of smell is based on the recognition of ligands and is mainly coordinated by odorant-binding proteins (OBPs). In Apis mellifera, behavior can be affected by different pathogens, including deformed wing virus (DWV) and its variants. In particular, it has been shown that variant A of DWV (DWV-A) is capable of altering the ultra-cellular structure associated with olfactory activity. In this study was evaluated olfactory sensitivity and the expression of OBP genes in honey bees inoculated with DWV-A. Electroantennographic analyses (EAG) were carried out to determine the olfactory sensitivity to the essential oils Eucalyptus globulus and Mentha piperita. The expression of nine antenna-specific OBP genes and DWV-A load in inoculated bees was also quantified by qPCR. We observed an inverse relationship between viral load and olfactory sensitivity and the expression of some OBP proteins. Thus, high viral loads reduced olfactory sensitivity to essential oils and the gene expression of the OBP2, OBP5, OBP11, and OBP12 proteins on the antennas of middle- and forager-age bees. These results suggest that DWV-A could have negative effects on the processes of aroma perception by worker bees, affecting their performance in tasks carried out in and outside the colony.

Varroa destructor from the Laboratory to the Field: Control, Biocontrol and IPM Perspectives-A Review

Varroa destructor is a real challenger for beekeepers and scientists: fragile out of the hive, tenacious inside a bee colony. From all the research done on the topic, we have learned that a better understanding of this organism in its relationship with the bee but also for itself is necessary. Its biology relies mostly on semiochemicals for reproduction, nutrition, or orientation. Many treatments have been developed over the years based on hard or soft acaricides or even on biocontrol techniques. To date, no real sustainable solution exists to reduce the pressure of the mite without creating resistances or harming honeybees. Consequently, the development of alternative disruptive tools against the parasitic life cycle remains open. It requires the combination of both laboratory and field results through a holistic approach based on health biomarkers. Here, we advocate for a more integrative vision of V. destructor research, where in vitro and field studies are more systematically compared and compiled. Therefore, after a brief state-of-the-art about the mite's life cycle, we discuss what has been done and what can be done from the laboratory to the field against V. destructor through an integrative approach.

Lipidomic Profiling Reveals Distinct Differences in Sphingolipids Metabolic Pathway between Healthy Apis cerana cerana larvae and Chinese Sacbrood Disease

Chinese sacbrood disease (CSD), which is caused by Chinese sacbrood virus (CSBV), is a major viral disease in Apis cerana cerana larvae. Analysis of lipid composition is critical to the study of CSBV replication. The host lipidome profiling during CSBV infection has not been conducted. This paper identified the lipidome of the CSBV-larvae interaction through high-resolution mass spectrometry. A total of 2164 lipids were detected and divided into 20 categories. Comparison of lipidome between healthy and CSBV infected-larvae showed that 266 lipid species were altered by CSBV infection. Furthermore, qRT-PCR showed that various sphingolipid enzymes and the contents of sphingolipids in the larvae were increased, indicating that sphingolipids may be important for CSBV infection. Importantly, Cer (d14:1 + hO/21:0 + O), DG (41:0e), PE (18:0e/18:3), SM (d20:0/19:1), SM (d37:1), TG (16:0/18:1/18:3), TG (18:1/20:4/21:0) and TG (43:7) were significantly altered in both CSBV_24 h vs. CK_24 h and CSBV_48 h vs. CK_48 h. Moreover, TG (39:6), which was increased by more than 10-fold, could be used as a biomarker for the early detection of CSD. This study provides evidence that global lipidome homeostasis in A. c. cerana larvae is remodeled after CSBV infection. Detailed studies in the future may improve the understanding of the relationship between the sphingolipid pathway and CSBV replication.

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.

Deformed wing virus variant shift from 2010 to 2016 in managed and feral UK honey bee colonies

Deformed wing virus (DWV) has been linked to the global decline of honey bees. DWV exists as three master variants (DWV-A, DWV-B, and DWV-C), each with differing outcomes for the honey bee host. Research in the USA showed a shift from DWV-A to DWV-B between 2010 to 2016 in honey bee colonies. Likewise, in the UK, a small study in 2007 found only DWV-A, whereas in 2016, DWV-B was the most prevalent variant. This suggests a shift from DWV-A to DWV-B might have occurred in the UK between 2007 and 2016. To investigate this further, data from samples collected in 2009/10 (n = 46) were compared to existing data from 2016 (n = 42). These samples also allowed a comparison of DWV variants between Varroa-untreated (feral) and Varroa-treated (managed) colonies. The results revealed that, in the UK, DWV-A was far more prevalent in 2009/10 (87%) than in 2016 (43%). In contrast, DWV-B was less prevalent in 2009/10 (76%) than in 2016 (93%). Regardless if colonies had been treated for Varroa (managed) or not (feral), the same trend from DWV-A to DWV-B occurred. Overall, the results reveal a decrease in DWV-A and an increase in DWV-B in UK colonies.

Adapted tolerance to virus infections in four geographically distinct Varroa destructor-resistant honeybee populations

The ectoparasitic mite, Varroa destructor, is unarguably the leading cause of honeybee (Apis mellifera) mortality worldwide through its role as a vector for lethal viruses, in particular, strains of the Deformed wing virus (DWV) and Acute bee paralysis virus (ABPV) complexes. Several honeybee populations across Europe have well-documented adaptations of mite-resistant traits but little is known about host adaptations towards the virus infections vectored by the mite. The aim of this study was to assess and compare the possible contribution of adapted virus tolerance and/or resistance to the enhanced survival of four well-documented mite-resistant honeybee populations from Norway, Sweden, The Netherlands and France, in relation to unselected mite-susceptible honeybees. Caged adult bees and laboratory reared larvae, from colonies of these four populations, were inoculated with DWV and ABPV in a series of feeding infection experiments, while control groups received virus-free food. Virus infections were monitored using RT-qPCR assays in individuals sampled over a time course. In both adults and larvae the DWV and ABPV infection dynamics were nearly identical in all groups, but all mite-resistant honeybee populations had significantly higher survival rates compared to the mite-susceptible honeybees. These results suggest that adapted virus tolerance is an important component of survival mechanisms.

Breeding for Virus Resistance and Its Effects on Deformed Wing Virus Infection Patterns in Honey Bee Queens

Viruses, and in particular the deformed wing virus (DWV), are considered as one of the main antagonists of honey bee health. The 'suppressed in ovo virus infection' trait (SOV) described for the first time that control of a virus infection can be achieved from genetically inherited traits and that the virus state of the eggs is indicative for this. This research aims to explore the effect of the SOV trait on DWV infections in queens descending from both SOV-positive (QDS+) and SOV-negative (QDS-) queens. Twenty QDS+ and QDS- were reared from each time four queens in the same starter-finisher colony. From each queen the head, thorax, ovaries, spermatheca, guts and eviscerated abdomen were dissected and screened for the presence of the DWV-A and DWV-B genotype using qRT-PCR. Queens descending from SOV-positive queens showed significant lower infection loads for DWV-A and DWV-B as well as a lower number of infected tissues for DWV-A. Surprisingly, differences were less expressed in the reproductive tissues, the ovaries and spermatheca. These results confirm that selection on the SOV trait is associated with increased virus resistance across viral genotypes and that this selection drives DWV towards an increased tissue specificity for the reproductive tissues. Further research is needed to explore the mechanisms underlying the interaction between the antiviral response and DWV.

A New Strain of Virus Discovered in China Specific to the Parasitic Mite Varroa destructor Poses a Potential Threat to Honey Bees

The ectoparasitic mite, Varroa destructor, feeds directly on honey bees and serves as a vector for transmitting viruses among them. The Varroa mite causes relatively little damage to its natural host, the Eastern honey bee (Apis cerana) but it is the most devastating pest for the Western honey bee (Apis mellifera). Using Illumina HiSeq sequencing technology, we conducted a metatranscriptome analysis of the microbial community associated with Varroa mites. This study led to the identification of a new Chinese strain of Varroa destructor virus-2 (VDV-2), which is a member of the Iflaviridae family and was previously reported to be specific to Varroa mites. A subsequent epidemiological investigation of Chinese strain of VDV-2 (VDV-2-China) showed that the virus was highly prevalent among Varroa populations and was not identified in any of the adult workers from both A. mellifera and A.cerana colonies distributed in six provinces in China, clearly indicating that VDV-2-China is predominantly a Varroa-adapted virus. While A. mellifera worker pupae exposed to less than two Varroa mites tested negative for VDV-2-China, VDV-2-China was detected in 12.5% of the A. mellifera worker pupae that were parasitized by more than 10 Varroa mites, bringing into play the possibility of a new scenario where VDV-2 could be transmitted to the honey bees during heavy Varroa infestations. Bioassay for the VDV-2-China infectivity showed that A. cerana was not a permissive host for VDV-2-China, yet A. mellifera could be a biological host that supports VDV-2-China's replication. The different replication dynamics of the virus between the two host species reflect their variation in terms of susceptibility to the virus infection, posing a potential threat to the health of the Western honey bee. The information gained from this study contributes to the knowledge concerning genetic variabilities and evolutionary dynamics of Varroa-borne viruses, thereby enhancing our understanding of underlying molecular mechanisms governing honey bee Varroosis.

Adaptation to vector-based transmission in a honeybee virus

Global pollinator declines as a result of emerging infectious diseases are of major concern. Managed honeybees Apis mellifera are susceptible to numerous parasites and pathogens, many of which appear to be transmissible to sympatric non-Apis taxa. The ectoparasitic mite Varroa destructor is considered to be the most significant threat to honeybees due to its role in vectoring RNA viruses, particularly Deformed wing virus (DWV). Vector transmission of DWV has resulted in the accumulation of high viral loads in honeybees and is often associated with colony death. DWV has two main genotypes, A and B. DWV-A was more prevalent during the initial phase of V. destructor establishment. In recent years, the global prevalence of DWV-B has increased, suggesting that DWV-B is better adapted to vector transmission than DWV-A. We aimed to determine the role vector transmission plays in DWV genotype prevalence at a colony level. We experimentally increased or decreased the number of V. destructor mites in honeybee colonies, and tracked DWV-A and DWV-B loads over a period of 10 months. Our results show that the two DWV genotypes differ in their response to mite numbers. DWV-A accumulation in honeybees was positively correlated with mite numbers yet DWV-A was largely undetected in the absence of the mite. In contrast, colonies had high loads of DWV-B even when mite numbers were low. DWV-B loads persisted in miticide-treated colonies, indicating that this genotype has a competitive advantage over DWV-A irrespective of mite numbers. Our findings suggest that the global increase in DWV-B prevalence is not driven by selective pressure by the vector. Rather, DWV-B is able to persist in colonies at higher viral loads relative to DWV-A in the presence and absence of V. destructor. The interplay between V. destructor and DWV genotypes within honeybee colonies may have broad consequences upon viral diversity in sympatric taxa as a result of spillover.

Pollinator community species richness dilutes prevalence of multiple viruses within multiple host species

Most pathogens are embedded in complex communities composed of multiple interacting hosts, but we are still learning how community-level factors, such as host diversity, abundance, and composition, contribute to pathogen spread for many host-pathogen systems. Evaluating relationships among multiple pathogens and hosts may clarify whether particular host or pathogen traits consistently drive links between community factors and pathogen prevalence. Pollinators are a good system to test how community composition influences pathogen spread because pollinator communities are extremely variable and contain several multi-host pathogens transmitted on shared floral resources. We conducted a field survey of four pollinator species to test the prevalence of three RNA viruses (deformed wing virus, black queen cell virus, and sacbrood virus) among pollinator communities with variable species richness, abundance, and composition. All three viruses showed a similar pattern of prevalence among hosts. Apis mellifera and Bombus impatiens had significantly higher viral prevalence than Lasioglossum spp. and Eucera pruinosa. In each species, lower virus prevalence was most strongly linked with greater pollinator community species richness. In contrast, pollinator abundance, species-specific pollinator abundance, and community composition were not associated with virus prevalence. Our results support a consistent dilution effect for multiple viruses and host species. Pollinators in species-rich communities had lower viral prevalence than pollinators from species-poor communities, when accounting for differences in pollinator abundance. Species-rich communities likely had lower viral prevalence because species-rich communities contained more native bee species likely to be poor viral hosts than species-poor communities, and all communities contained the highly competent hosts A. mellifera and B. impatiens. Interestingly, the strength of the dilution effect was not consistent among hosts. Instead, host species with low viral prevalence exhibited weaker dilution effects compared to hosts with high viral prevalence. Therefore, host species susceptibility and competence for each virus may contribute to variation in the strength of dilution effects. This study expands biodiversity-disease studies to the pollinator-virus system, finding consistent evidence of the dilution effect among multiple similar pathogens that infect "replicate" host communities.

Virus Prospecting in Crickets-Discovery and Strain Divergence of a Novel Iflavirus in Wild and Cultivated Acheta domesticus

Orthopteran insects have high reproductive rates leading to boom-bust population dynamics with high local densities that are ideal for short, episodic disease epidemics. Viruses are particularly well suited for such host population dynamics, due to their supreme ability to adapt to changing transmission criteria. However, very little is known about the viruses of Orthopteran insects. Since Orthopterans are increasingly reared commercially, for animal feed and human consumption, there is a risk that viruses naturally associated with these insects can adapt to commercial rearing conditions, and cause disease. We therefore explored the virome of the house cricket Acheta domesticus, which is both part of the natural Swedish landscape and reared commercially for the pet feed market. Only 1% of the faecal RNA and DNA from wild-caught A. domesticus consisted of viruses. These included both known and novel viruses associated with crickets/insects, their bacterial-fungal microbiome, or their plant food. Relatively abundant among these viral Operational Taxonomic Units (OTUs) was a novel Iflavirus, tentatively named Acheta domesticus Iflavirus (AdIV). Quantitative analyses showed that AdIV was also abundant in frass and insect samples from commercially reared crickets. Interestingly, the wild and commercial AdIV strains had short, extremely divergent variation hotspots throughout the genome, which may indicate specific adaptation to their hosts’ distinct rearing environments.

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.

Direct Evidence for Infection of Varroa destructor Mites with the Bee-Pathogenic Deformed Wing Virus Variant B - but Not Variant A - via Fluorescence-in situ-Hybridization Analysis

Deformed wing virus (DWV) is a bee pathogenic, single- and positive-stranded RNA virus that has been involved in severe honey bee colony losses worldwide. DWV, when transmitted horizontally or vertically from bee to bee, causes mainly covert infections not associated with any visible symptoms or damage. Overt infections occur after vectorial transmission of DWV to the developing bee pupae through the ectoparasitic mite Varroa destructor Symptoms of overt infections are pupal death, bees emerging with deformed wings and shortened abdomens, or cognitive impairment due to brain infection. So far, three variants of DWV, DWV-A, DWV-B, and DWV-C, have been described. While it is widely accepted that V. destructor acts as vector of DWV, the question of whether the mite only functions as a mechanical vector or whether DWV can infect the mite thus using it as a biological vector is hotly debated, because in the literature data can be found that support both hypotheses. In order to settle this scientific dispute, we analyzed putatively DWV-infected mites with a newly established protocol for fluorescence-in situ-hybridization of mites and demonstrated DWV-specific signals inside mite cells. We provide compelling and direct evidence that DWV-B infects the intestinal epithelium and the salivary glands of V. destructor In contrast, no evidence for DWV-A infecting mite cells was found. Our data are key to understanding the pathobiology of DWV, the mite's role as a biological DWV vector and the quasispecies dynamics of this RNA virus when switching between insect and arachnid host species.IMPORTANCE Deformed wing virus (DWV) is a bee pathogenic, originally rather benign, single- and positive-stranded RNA virus. Only the vectorial transmission of this virus to honey bees by the ectoparasitic mite Varroa destructor leads to fatal or symptomatic infections of individuals, usually followed by collapse of the entire colony. Studies on whether the mite only acts as a mechanical virus vector or whether DWV can infect the mite and thus use it as a biological vector have led to disparate results. In our study using fluorescence-in situ-hybridization we provide compelling and direct evidence that at least the DWV-B variant infects the gut epithelium and the salivary glands of V. destructor Hence, the host range of DWV includes both, bees (Insecta) and mites (Arachnida). Our data contribute to a better understanding of the triangular relationship between honey bees, V. destructor and DWV and the evolution of virulence in this viral bee pathogen.

Deformed wing virus type a and b in managed honeybee colonies of Argentina

Apis mellifera is infected by more than 24 virus species worldwide, mainly positive-sense, single-stranded RNA viruses of the Dicistroviridae and Iflaviridae families. Among the viruses that infect honey bees, Deformed wing virus is the most prevalent and is present as three master variants DWV-A, B, and C. Given that the ectoparasitic mite Varroa destructor vectors these virus variants, recombination events between them are expected, and variants and their recombinants can co-exist in mites and honeybees at the same time. In this study, we detect, through RT-qPCR, the presence of DWV-A and B in the same samples of adult bees from colonies of Argentina. Total RNA was extracted from pools of ten adult bees from 45 apiaries distributed across the main beekeeping Provinces of Argentina (Buenos Aires, Santa Fe, Córdoba, Santiago del Estero, Río Negro, and Mendoza); then RT-qPCR reactions were performed to detect DWV-A and B, with specific primer pairs. After the amplifications, PCR products (204 and 660 bp amplicons for DWV-B, and ~250 bp for DWV-A) were purified and sequenced to verify that they corresponded to reported sequences, analyzing them using the Blast software. Of the 45 samples analyzed by RT-qPCR, over 90% were infected with DWV-A and 47% were also positive for DWV-B, where it was found in high prevalence specifically in colonies of A. mellifera of the Buenos Aires Province. Future studies will determine the impact of this type of the virus and its ability to recombine with the other DWV types in the apiaries of our country.

A Novel Iflavirus Was Discovered in Green Rice Leafhopper Nephotettix cincticeps and Its Proliferation Was Inhibited by Infection of Rice Dwarf Virus

The green rice leafhopper, Nephotettix cincticeps (Hemiptera: Cicadellidae), is a key insect vector transmitting rice dwarf virus (RDV) that causes rice dwarf disease. We discovered a novel iflavirus from the transcriptomes of N. cincticeps and named it as Nephotettix cincticeps positive-stranded RNA virus-1 (NcPSRV-1). The viral genome consists of 10,524 nucleotides excluding the poly(A) tail and contains one predicted open reading frame encoding a polyprotein of 3,192 amino acids, flanked by 5' and 3' untranslated regions. NcPSRV-1 has a typical iflavirus genome arrangement and is clustered with the members of the family Iflaviridae in the phylogenetic analysis. NcPSRV-1 was detected in all tested tissues and life stages of N. cincticeps and could be transmitted horizontally and vertically. Moreover, NcPSRV-1 had high prevalence in the laboratory populations and was widely spread in field populations of N. cincticeps. NcPSRV-1 could also infect the two-striped leafhopper, Nephotettix apicalis, at a 3.33% infection rate, but was absent in the zigzag leafhopper, Recilia dorsalis, and rice Oryza sativa variety TN1. The infection of RDV altered the viral load and infection rate of NcPSRV-1 in N. cincticeps, for which it seems that RDV has an antagonistic effect on NcPSRV-1 infection in the host.

Deformed wing virus: using reverse genetics to tackle unanswered questions about the most important viral pathogen of honey bees

Deformed wing virus (DWV) is the most important viral pathogen of honey bees. It usually causes asymptomatic infections but, when vectored by the ectoparasitic mite Varroa destructor, it is responsible for the majority of overwintering colony losses globally. Although DWV was discovered four decades ago, research has been hampered by the absence of an in vitro cell culture system or the ability to culture pure stocks of the virus. The recent developments of reverse genetic systems for DWV go some way to addressing these limitations. They will allow the investigation of specific questions about strain variation, host tropism and pathogenesis to be answered, and are already being exploited to study tissue tropism and replication in Varroa and non-Apis pollinators. Three areas neatly illustrate the advances possible with reverse genetic approaches: (i) strain variation and recombination, in which reverse genetics has highlighted similarities rather than differences between virus strains; (ii) analysis of replication kinetics in both honey bees and Varroa, in studies that likely explain the near clonality of virus populations often reported; and (iii) pathogen spillover to non-Apis pollinators, using genetically tagged viruses to accurately monitor replication and infection.

OneHealth implications of infectious diseases of wild and managed bees

The OneHealth approach aims to further our understanding of the drivers of human, animal and environmental health, and, ultimately, to improve them by combining approaches and knowledge from medicine, biology and fields beyond. Wild and managed bees are essential pollinators of crops and wild flowers. Their health thus directly impacts on human and environmental health. At the same time, these bee species represent highly amenable and relevant model organisms for a OneHealth approach that aims to study fundamental epidemiological questions. In this review, we focus on how infectious diseases of wild and managed bees can be used as a OneHealth model system, informing fundamental questions on ecological immunology and disease transmission, while addressing how this knowledge can be used to tackle the issues facing pollinator health.

RNAseq of Deformed Wing Virus and Other Honey Bee-Associated Viruses in Eight Insect Taxa with or without Varroa Infestation

The global spread of a parasitic mite (Varroa destructor) has resulted in Deformed wing virus (DWV), a previously rare pathogen, now dominating the viromes in honey bees and contributing to large-scale honey bee colony losses. DWV can be found in diverse insect taxa and has been implicated in spilling over from honey bees into associated ("apiary") and other ("non-apiary") insects. Here we generated next generation sequence data from 127 insect samples belonging to diverse taxa collected from Hawaiian islands with and without Varroa to identify whether the mite has indirectly affected the viral landscapes of key insect taxa across bees, wasps, flies and ants. Our data showed that, while Varroa was associated with a dramatic increase in abundance of (predominantly recombinant) DWV in honey bees (and no other honey bee-associated RNA virus), this change was not seen in any other taxa sampled. Honey bees share their environment with other insect populations and exist as a homogenous group, frequently sharing common viruses, albeit at low levels. Our data suggest that the threat of Varroa to increase viral load in an apiary does not automatically translate to an increase in virus load in other insects living in the wider community.