Deformed wing virus: replication and viral load in mites (Varroa destructor)

Representative honey bee viruses do not replicate in the small hive beetle, Aethina tumida Murray

The small hive beetle (SHB), Aethina tumida Murray, is an invasive pest of the honey bee and causes significant damage through the consumption of colony resources and brood. Two assumptions related to honey bee virus transmission have been made about SHB: first, that SHB vectors honey bee viruses and second, that these viruses replicate in SHB based on the detection of both positive and negative strand viral genomic RNA within SHB. To clarify the role of SHB in virus transmission, we sought to address whether selected honey bee viruses replicate in SHB. Sequences derived from five honey bee viruses were identified in the transcriptomes of field-caught SHB from the U.S., but not in those of lab-reared SHB, suggesting that these viruses do not replicate in SHB. To elucidate whether the representative viruses, Israeli acute paralysis virus (IAPV; Dicistroviridae) and Deformed wing virus (DWV; Iflaviridae) replicate in SHB, we tested for replication in vitro in an SHB-derived cell line (BCIRL-AtumEN-1129-D6). Following treatment of the cell line with viral particles or viral RNA, the number of virus genomes was monitored by reverse transcription quantitative PCR (RT-qPCR). In contrast to the positive control, IAPV and DWV RNA levels steadily decreased over a period of 8 days. Collectively, these results from bioinformatic observations and in vitro experiments indicate that IAPV and DWV do not replicate in SHB. These results are consistent with the host specificity of most insect viruses within a single insect order and indicate that while SHB may serve as a mechanical vector of honey bee viruses within and between hives, this insect does not serve as a biological vector for these honey bee viruses.

Sacbrood viruses and select Lake Sinai virus variants dominated Apis mellifera colonies symptomatic for European foulbrood

European foulbrood (EFB) is a prevalent disease in the European honey bee (Apis mellifera) in the United States, which can lead to colony decline and collapse. The bacterial components of EFB are well-studied, but the diversity of viral infections within infected colonies has not been explored. In this study, we use meta-transcriptomics sequencing of 12 honey bee hives, symptomatic (+, n = 6) and asymptomatic (-, n = 6) for EFB, to investigate viral infection associated with the disease. We assembled 41 viral genomes, belonging to three families (Iflaviridae, Dicistroviridae, and Sinhaliviridae), all previously reported in honey bees, including Lake Sinai virus, deformed wing virus, sacbrood virus, Black queen cell virus, and Israeli acute paralysis virus. In colonies with severe EFB, we observed a higher occurrence of viral genomes (34 genomes) in contrast to fewer recovered from healthy colonies (seven genomes) and a complete absence of Dicistroviridae genomes.We observed specific Lake Sinai virus clades associated exclusively with EFB + or EFB - colonies, in addition to EFB-afflicted colonies that exhibited an increase in relative abundance of sacbrood viruses. Multivariate analyses highlighted that a combination of site and EFB disease status influenced RNA virome composition, while EFB status alone did not significantly impact it, presenting a challenge for comparisons between colonies kept in different yards. These findings contribute to the understanding of viral dynamics in honey bee colonies compromised by EFB and underscore the need for future investigations to consider viral composition when investigating EFB.IMPORTANCEThis study on the viromes of honey bee colonies affected by European foulbrood (EFB) sheds light on the dynamics of viral populations in bee colonies in the context of a prevalent bacterial brood disease. The identification of distinct Lake Sinai virus and sacbrood virus clades associated with colonies affected by severe EFB suggests a potential connection between viral composition and disease status, emphasizing the need for further investigation into the role of viruses during EFB infection. The observed increase in sacbrood viruses during EFB infection suggests a potential viral dysbiosis, with potential implications for honey bee brood health. These findings contribute valuable insights related to beekeeping practices, offering a foundation for future research aimed at understanding and mitigating the impact of bacterial and viral infection in commercial honey bee operations and the management of EFB.

Development of a multiplex real-time quantitative reverse-transcription polymerase chain reaction for the detection of four bee viruses

Viruses in the families Dicistroviridae and Iflaviridae are among the main threats to western honey bees (Apis mellifera) and native bee species. Polymerase chain reaction (PCR) is the gold standard for pathogen detection in bees. However, high throughput screening for bee virus infections in singleplex PCR reactions is cumbersome and limited by the high quantities of sample RNA required. Thus, the development of a sensitive and specific multiplex PCR detection method for screening for multiple viruses simultaneously is necessary. Here, we report the development of a one-step multiplex reverse-transcription quantitative polymerase chain reaction (RT-qPCR) assay to detect four viruses commonly encountered in pollinator species. The optimized multiplex RT-qPCR protocol described in this study allows simultaneous detection of two dicistroviruses (Israeli acute paralysis virus and Black queen cell virus) and two iflaviruses (Sacbrood virus and Deformed wing virus) with high efficiency and specificity comparable to singleplex detection assays. This assay provides a broad range of detection and quantification, and the results of virus quantification in this study are similar to those performed in other studies using singleplex detection assays. This method will be particularly useful for data generation from small-bodied insect species that yield low amounts of RNA.

Oxalic Acid Treatment: Short-Term Effects on Enzyme Activities, Vitellogenin Content, and Residual Oxalic Acid Content in House Bees, Apis mellifera L

Honeybees (Apis mellifera L.) have to face many challenges, including Varroa destructor infestation, associated with viral transmission. Oxalic acid is one of the most common treatments against Varroa. Little is known about the physiological effects of oxalic acid, especially those on honeybees' immune systems. In this study, the short-term effects (0-96 h) of oxalic acid treatment on the immune system components (i.e., glucose oxidase, phenoloxidase, glutathione S-transferase, catalase activities, and vitellogenin contents) of house bees were preliminarily investigated. Oxalic acid contents of bee bodies and haemolymphs were also measured. The results confirm that oxalic acid is constitutively present in bee haemolymphs and its concentration is not affected by treatment. At 6 h after the treatment, a maximum peak of oxalic acid content was detected on bees' bodies, which gradually decreased after that until physiological levels were reached at 48 h. In the immune system, the oxalic acid treatment determined a peak in glucose oxidase activity at 48 h, indicating a potential defence response and an increase in vitellogenin content at 24 h. No significant changes were recorded in phenoloxidase, glutathione S-transferase, and catalase activities. These results suggest a time-dependent response to oxalic acid, with potential immune system activation in treated bees.

Ecological and social factors influence interspecific pathogens occurrence among bees

The interspecific transmission of pathogens can occur frequently in the environment. Among wild bees, the main spillover cases are caused by pathogens associated with Apis mellifera, whose colonies can act as reservoirs. Due to the limited availability of data in Italy, it is challenging to accurately assess the impact and implications of this phenomenon on the wild bee populations. In this study, a total of 3372 bees were sampled from 11 Italian regions within the BeeNet project, evaluating the prevalence and the abundance of the major honey bee pathogens (DWV, BQCV, ABPV, CBPV, KBV, Nosema ceranae, Ascosphaera apis, Crithidia mellificae, Lotmaria passim, Crithidia bombi). The 68.4% of samples were positive for at least one pathogen. DWV, BQCV, N. ceranae and CBPV showed the highest prevalence and abundance values, confirming them as the most prevalent pathogens spread in the environment. For these pathogens, Andrena, Bombus, Eucera and Seladonia showed the highest mean prevalence and abundance values. Generally, time trends showed a prevalence and abundance decrease from April to July. In order to predict the risk of infection among wild bees, statistical models were developed. A low influence of apiary density on pathogen occurrence was observed, while meteorological conditions and agricultural management showed a greater impact on pathogen persistence in the environment. Social and biological traits of wild bees also contributed to defining a higher risk of infection for bivoltine, communal, mining and oligolectic bees. Out of all the samples tested, 40.5% were co-infected with two or more pathogens. In some cases, individuals were simultaneously infected with up to five different pathogens. It is essential to increase knowledge about the transmission of pathogens among wild bees to understand dynamics, impact and effects on pollinator populations. Implementing concrete plans for the conservation of wild bee species is important to ensure the health of wild and human-managed bees within a One-Health perspective.

Distribution of infectious and parasitic agents among three sentinel bee species across European agricultural landscapes

Infectious and parasitic agents (IPAs) and their associated diseases are major environmental stressors that jeopardize bee health, both alone and in interaction with other stressors. Their impact on pollinator communities can be assessed by studying multiple sentinel bee species. Here, we analysed the field exposure of three sentinel managed bee species (Apis mellifera, Bombus terrestris and Osmia bicornis) to 11 IPAs (six RNA viruses, two bacteria, three microsporidia). The sentinel bees were deployed at 128 sites in eight European countries adjacent to either oilseed rape fields or apple orchards during crop bloom. Adult bees of each species were sampled before their placement and after crop bloom. The IPAs were detected and quantified using a harmonised, high-throughput and semi-automatized qPCR workflow. We describe differences among bee species in IPA profiles (richness, diversity, detection frequencies, loads and their change upon field exposure, and exposure risk), with no clear patterns related to the country or focal crop. Our results suggest that the most frequent IPAs in adult bees are more appropriate for assessing the bees' IPA exposure risk. We also report positive correlations of IPA loads supporting the potential IPA transmission among sentinels, suggesting careful consideration should be taken when introducing managed pollinators in ecologically sensitive environments.

Integrated Pest Management Strategies to Control Varroa Mites and Their Effect on Viral Loads in Honey Bee Colonies

Honey bee viruses in combination with varroa mite are very damaging for honey bee colonies worldwide. There are no effective methods to control the viral load in honey bee colonies except regular and effective control of mites. Integrated Pest Management strategies are required to effectively control mites with veterinary medicines based on organic compounds. We evaluated the effect of two brood interruption techniques, queen caging (QC) and trapping comb (TC), followed by an oxalic acid treatment, on the mite fall, colony strength, and viral load of Deformed Wing Virus (DWV) and Acute Bee Paralysis Virus (ABPV). In this paper, we report the data obtained in two experimental sites, in Slovenia and Italy, in terms of the varroacide efficacy, colony strength, and viral load. The number of adult bees after the adoption of the two techniques showed similar decreasing trends in both locations. The viral load of Acute Bee Paralysis Virus did not show any significant reduction after 25 days, reported as the number of Real-Time PCR cycles needed to detect the virus. The viral load of DWV also did not show a significant reduction after 25 days. The acaricidal efficacy of the applied protocols was high in both experimental groups and in both apiaries. Both the queen caging and trapping comb techniques, followed by an oxalic acid treatment, can be considered effective varroa treatment strategies, but further studies should be carried out to evaluate the long-term effects on viral loads to plan the Integrated Pest Management strategy with the right timing before wintering.

Susceptible and infectious states for both vector and host in a dynamic pathogen-vector-host system

Deformed wing virus (DWV) is a resurgent insect pathogen of honeybees that is efficiently transmitted by vectors and through host social contact. Continual transmission of DWV between hosts and vectors is required to maintain the pathogen within the population, and this vector-host-pathogen system offers unique disease transmission dynamics for pathogen maintenance between vectors and a social host. In a series of experiments, we measured vector-vector, host-host and host-vector transmission routes and show how these maintain DWV in honeybee populations. We found co-infestations on shared hosts allowed for movement of DWV from mite to mite. Additionally, two social behaviours of the honeybee, trophallaxis and cannibalization of pupae, provide routes for horizontal transmission from bee to bee. Circulation of the virus solely among hosts through communicable modes provides a reservoir of DWV for naïve Varroa to acquire and subsequently vector the pathogen. Our findings illustrate the importance of community transmission between hosts and vector transmission. We use these results to highlight the key avenues used by DWV during maintenance and infection and point to similarities with a handful of other infectious diseases of zoonotic and medical importance.

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.

Effects of Deformed Wing Virus-Targeting dsRNA on Viral Loads in Bees Parasitised and Non-Parasitised by Varroa destructor

The Varroa destructor mite is a devastating parasite of honey bees; however the negative effects of varroa parasitism are exacerbated by its role as an efficient vector of the honey bee pathogen, Deformed wing virus (DWV). While no direct treatment for DWV infection is available for beekeepers to use on their hives, RNA interference (RNAi) has been widely explored as a possible biopesticide approach for a range of pests and pathogens. This study tested the effectiveness of three DWV-specific dsRNA sequences to lower DWV loads and symptoms in honey bees reared from larvae in laboratory mini-hives containing bees and varroa. The effects of DWV-dsRNA treatment on bees parasitised and non-parasitised by varroa mites during development were investigated. Additionally, the impact of DWV-dsRNA on viral loads and gene expression in brood-parasitising mites was assessed using RNA-sequencing. Bees parasitised during development had significantly higher DWV levels compared to non-parasitised bees. However, DWV-dsRNA did not significantly reduce DWV loads or symptoms in mini-hive reared bees, possibly due to sequence divergence between the DWV variants present in bees and varroa and the specific DWV-dsRNA sequences used. Varroa mites from DWV-dsRNA treated mini-hives did not show evidence of an elevated RNAi response or significant difference in DWV levels. Overall, our findings show that RNAi is not always successful, and multiple factors including pathogen diversity and transmission route may impact its efficiency.

Age-related response to mite parasitization and viral infection in the honey bee suggests a trade-off between growth and immunity

Host age at parasites' exposure is often neglected in studies on host-parasite interactions despite the important implications for epidemiology. Here we compared the impact of the parasitic mite Varroa destructor, and the associated pathogenic virus DWV on different life stages of their host, the western honey bee Apis mellifera. The pre-imaginal stages of the honey bee proved to be more susceptible to mite parasitization and viral infection than adults. The higher viral load in mite-infested bees and DWV genotype do not appear to be the drivers of the observed difference which, instead, seems to be related to the immune-competence of the host. These results support the existence of a trade-off between immunity and growth, making the pupa, which is involved in the highly energy-demanding process of metamorphosis, more susceptible to parasites and pathogens. This may have important implications for the evolution of the parasite's virulence and in turn for honey bee health. Our results highlight the important role of host's age and life stage at exposure in epidemiological modelling. Furthermore, our study could unravel new aspects of the complex honey bee-Varroa relationship to be addressed for a sustainable management of this parasite.

Prevalence of Adult Honey Bee (Apis mellifera L.) Pests and Pathogens in the Five Beekeeping Regions of Mexico

Mexico is a major honey producer, but not much information exists about the health status of honey bees (Apis mellifera L.) in the country. This study was conducted to determine the sanitary status of adult honey bees in Mexico's five beekeeping regions. Samples from 369 apiaries were diagnosed to identify pathogens such as Varroa destructor, which was quantified, Acarapis woodi, Nosema spp., and five viruses. Colonies were also inspected for the presence of the small hive beetle (SHB), Aethina tumida. Varroa destructor was found in 83.5% of the apiaries, with the Pacific Coast region having the highest prevalence (>95%) and rates (4.5% ± 0.6). Acarapis woodi was detected in only one apiary from the Pacific Coast, whereas Nosema spp. were prevalent in 48.5% of the apiaries, with the highest and lowest frequencies in the Yucatan Peninsula and North regions (64.6% and 10.2%, respectively). For viruses, deformed wing virus (DWV) was detected in 26.1% of the apiaries, with the highest frequency in the Pacific Coast region (44.7%). Israeli acute paralysis virus (IAPV) was diagnosed in 3.2% of the samples and sacbrood bee virus (SBV) in 23.3% of them, with the highest frequency in the High Plateau region (36.4%). Chronic bee paralysis and Kashmir bee viruses were not detected. SHB prevalence was 25.2% nationwide, with the highest frequency in the Yucatan Peninsula (39.2%). This study shows that the most common parasites of adult honey bees in Mexico are V. destructor and Nosema spp., and that the most prevalent virus is DWV, whereas SHB is highly prevalent in the Yucatan Peninsula. This information could be useful to design disease control strategies for honey bee colonies in different regions of Mexico.

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.

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.

A Novel Variant of Deformed Wing Virus (DWV) from the Invasive Honeybee Apis florea (Apidae, Hymenoptera) and Its Ectoparasite Euvarroa sinhai (Acarina, Mesostigmata) in Taiwan

The invasion of Apis florea in Taiwan was first recorded in 2017. The deformed wing virus (DWV) has been identified as a common bee virus in apiculture around the world. Ectoparasitic mites are the main DWV vector for horizontal transmission. However, there are few studies about the ectoparasitic mite of Euvarroa sinhai, which has been found in A. florea. In this study, the prevalence of DWV among four hosts, including A. florea, Apis mellifera, E. sinhai, and Varroa destructor, was determined. The results showed that a high DWV-A prevalence rate in A. florea, ranging from 69.2% to 94.4%, was detected. Additionally, the genome of DWV isolates was sequenced and subjected to phylogenetic analysis based on the complete polyprotein sequence. Furthermore, isolates from A. florea and E. sinhai both formed a monophyletic group for the DWV-A lineage, and the sequence identity was 88% between the isolates and DWV-A reference strains. As noted above, two isolates could be the novel DWV strain. It cannot be excluded that novel DWV strains could pose an indirect threat to sympatric species, such as A. mellifera and Apis cerana.

Molecular and kinetic properties of three acetylcholinesterases in the Varroa mite, Varroa destructor

The Varroa mite, Varroa destructor, poses one of the most serious threats to honey bees worldwide. Although coumaphos, an anticholinesterase pesticide, is widely used for varroa mite control, little information is available on the properties of Varroa mite acetylcholinesterases (VdAChEs). In this study, three putative VdAChEs were annotated and named VdAChE1, VdAChE2, and VdAChE3. All VdAChEs possessed most of the functionally important signature domains, suggesting that they are catalytically active. Phylogenetic analysis revealed that VdAChE1 was clustered into a clade containing most arthropod AChE1s, whereas VdAChE2 and VdAChE3 formed a unique clade with other arachnid AChEs. VdAChE1 was determined to be membrane-anchored, but both VdAChE2 and VdAChE3 are soluble, as judged by electrophoresis in conjunction with western blotting. Tissue-specific transcription profiling revealed that VdAChE1 was most predominantly expressed in the synganglion. In contrast, VdAChE2 was most predominantly expressed in the legs and cuticle. VdAChE3 showed negligible expression levels in all the tissues examined. In a kinetic analysis using recombinant VdAChEs, VdAChE1 exhibited the highest catalytic efficiency, followed by VdAChE2 and VdAChE3. Inhibition experiments revealed that VdAChE1 was most sensitive to all tested inhibitors. Taken together, VdAChE1 appears to be the major synaptic enzyme with a more toxicological relevance, whereas VdAChE2 is involved in other noncatalytic functions, including chemical defense against xenobiotics. Current findings contribute to a more detailed understanding of the evolutionary and functional traits of VdAChEs and to the design of novel anticholinesterase varroacides.

Viral species differentially influence macronutrient preferences based on honey bee genotype

Food quantity and macronutrients contribute to honey bee health and colony survival by mediating immune responses. We determined if this held true for bees injected with chronic bee paralysis virus (CBPV) and deformed wing virus (DWV), two common honey bee ssRNA viruses. Pollen-substitute diet and syrup consumption rates and macronutrient preferences of two Varroa-resistant stocks (Pol-Line and Russian bees) were compared to Varroa-susceptible Italian bees. Bee stocks varied in consumption, where Italian bees consumed more than Pol-Line and Russian bees. However, the protein: lipid (P:L) ratios of diet consumed by the Italian and Russian bees was greater than that of the Pol-Line bees. Treatment had different effects on consumption based on the virus injected. CBPV was positively correlated with syrup consumption, while DWV was not correlated with consumption. P:L ratios of consumed diet were significantly impacted by the interaction of bee stock and treatment, with the trends differing between CBPV and DWV. Variation in macronutrient preferences based on viral species may indicate differences in energetic costs associated with immune responses to infections impacting different systems. Further, virus species interacted with bee genotype, indicating different mechanisms of viral resistance or tolerance among honey bee genotypes.

The Epidemiological Situation of the Managed Honey Bee (Apis mellifera) Colonies in the Italian Region Emilia-Romagna

The recent decades witnessed the collapse of honey bee colonies at a global level. The major drivers of this collapse include both individual and synergic pathogen actions, threatening the colonies’ survival. The need to define the epidemiological pattern of the pathogens that are involved has led to the establishment of monitoring programs in many countries, Italy included. In this framework, the health status of managed honey bees in the Emilia–Romagna region (northern Italy) was assessed, throughout the year 2021, on workers from 31 apiaries to investigate the presence of major known and emerging honey bee pathogens. The prevalence and abundance of DWV, KBV, ABPV, CBPV, Nosema ceranae, and trypanosomatids (Lotmaria passim, Crithidia mellificae, Crithidia bombi) were assessed by molecular methods. The most prevalent pathogen was DWV, followed by CBPV and N. ceranae. Trypanosomatids were not found in any of the samples. Pathogens had different peaks in abundance over the months, showing seasonal trends that were related to the dynamics of both bee colonies and Varroa destructor infestation. For some of the pathogens, a weak but significant correlation was observed between abundance and geographical longitude. The information obtained in this study increases our understanding of the epidemiological situation of bee colonies in Emilia–Romagna and helps us to implement better disease prevention and improved territorial management of honey bee health.

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.

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.

Genotype, but Not Climate, Affects the Resistance of Honey Bees (Apis mellifera) to Viral Infections and to the Mite Varroa destructor

This study was conducted to analyze the effect of genotype and climate on the resistance of honey bee (Apis mellifera) colonies to parasitic and viral diseases. The prevalence and intensity of parasitism by Varroa destructor, or infection by Nosema spp., and four honey bee viruses were determined in 365 colonies of predominantly European or African ancestry (descendants of A. m. scutellata) in subtropical and temperate regions of Mexico. Varroa destructor was the most prevalent parasite (95%), whilst N. ceranae was the least prevalent parasite (15%). Deformed wing virus (DWV) and black queen cell virus (BQCV) were the only viruses detected, at frequencies of 38% and 66%, respectively. Varroa destructor was significantly more prevalent in colonies of European ancestry (p < 0.05), and the intensity of parasitism by V. destructor or infection by DWV and BQCV was also significantly higher in colonies of European descent than in African descent colonies (p < 0.01), although no genotype−parasite associations were found for N. ceranae. Additionally, significant and positive correlations were found between V. destructor and DWV levels, and the abundance of these pathogens was negatively correlated with the African ancestry of colonies (p < 0.01). However, there were no significant effects of environment on parasitism or infection intensity for the colonies of both genotypes. Therefore, it is concluded that the genotype of honey bee colonies, but not climate, influences their resistance to DWV, BQCV, and V. destructor.

Honey Bee Genetic Stock Determines Deformed Wing Virus Symptom Severity but not Viral Load or Dissemination Following Pupal Exposure

Honey bees exposed to Varroa mites incur substantial physical damage in addition to potential exposure to vectored viruses such as Deformed wing virus (DWV) that exists as three master variants (DWV-A, DWV-B, and DWV-C) and recombinants. Although mite-resistant bees have been primarily bred to mitigate the impacts of Varroa mites, mite resistance may be associated with increased tolerance or resistance to the vectored viruses. The goal of our study is to determine if five honey bee stocks (Carniolan, Italian, Pol-Line, Russian, and Saskatraz) differ in their resistance or tolerance to DWV based on prior breeding for mite resistance. We injected white-eyed pupae with a sublethal dose (105) of DWV or exposed them to mites and then evaluated DWV levels and dissemination and morphological symptoms upon adult emergence. While we found no evidence of DWV resistance across stocks (i.e., similar rates of viral replication and dissemination), we observed that some stocks exhibited reduced symptom severity suggestive of differential tolerance. However, DWV tolerance was not consistent across mite-resistant stocks as Russian bees were most tolerant, while Pol-Line exhibited the most severe symptoms. DWV variants A and B exhibited differential dissemination patterns that interacted significantly with the treatment group but not bee stock. Furthermore, elevated DWV-B levels reduced adult emergence time, while both DWV variants were associated with symptom likelihood and severity. These data indicate that the genetic differences underlying bee resistance to Varroa mites are not necessarily correlated with DWV tolerance and may interact differentially with DWV variants, highlighting the need for further work on mechanisms of tolerance and bee stock-specific physiological interactions with pathogen variants.

Viral communities in the parasite Varroa destructor and in colonies of their honey bee host (Apis mellifera) in New Zealand

The parasitic mite Varroa destructor is a leading cause of mortality for Western honey bee (Apis mellifera) colonies around the globe. We sought to confirm the presence and likely introduction of only one V. destructor haplotype in New Zealand, and describe the viral community within both V. destructor mites and the bees that they parasitise. A 1232 bp fragment from mitochondrial gene regions suggests the likely introduction of only one V. destructor haplotype to New Zealand. Seventeen viruses were found in bees. The most prevalent and abundant was the Deformed wing virus A (DWV-A) strain, which explained 95.0% of the variation in the viral community of bees. Black queen cell virus, Sacbrood virus, and Varroa destructor virus 2 (VDV-2) played secondary roles. DWV-B and the Israeli acute paralysis virus appeared absent from New Zealand. Ten viruses were observed in V. destructor, with > 99.9% of viral reads from DWV-A and VDV-2. Substantially more variation in viral loads was observed in bees compared to mites. Where high levels of VDV-2 occurred in mites, reduced DWV-A occurred in both the mites and the bees co-occurring within the same hive. Where there were high loads of DWV-A in mites, there were typically high viral loads in bees.

Colony-Level Viral Load Influences Collective Foraging in Honey Bees

Nutrition is an important component of social insect colony health especially in the face of stressors such as parasitism and viral infections. Honey bees are known to preferentially select nectar and pollen based on macronutrient and phytochemical contents and in response to pathogen loads. However, given that honey bees live in colonies, collective foraging decisions may be impacted directly by forager infection status but also by colony health. This field experiment was conducted to determine if honey bee viral infections are correlated with pollen and nectar foraging and if these associations are impacted more by colony or forager infection. By comparing regressions with and without forager and colony variables and through structural equation models, we were able to determine the relative contributions of colony and forager virus loads on forager decisions. We found that foragers had higher numbers and levels of BQCV and CBPV but lower levels of DWV viruses than their respective colonies. Overall, individuals appeared to forage based a combination of their own and colony health but with greater weight given to colony metrics. Colony parasitism by Varroa mites, positively correlated with both forager and colony DWV-B levels, was negatively associated with nectar weight. Further, colony DWV-B levels were negatively associated with individually foraged pollen protein: lipid ratios but positively correlated with nectar weight and sugar content. This study shows that both colony and forager health can simultaneously mediate individual foraging decisions and that the importance of viral infections and parasite levels varies with foraging metrics. Overall, this work highlights the continued need to explore the interactions of disease, nutrition, and genetics in social interactions and structures.

An Insight Into the microRNA Profile of the Ectoparasitic Mite Varroa destructor (Acari: Varroidae), the Primary Vector of Honey Bee Deformed Wing Virus

The remarkably adaptive mite Varroa destructor is the most important honey bee ectoparasite. Varroa mites are competent vectors of deformed wing virus (DWV), and the Varroa-virus complex is a major determinant of annual honey bee colony mortality and collapse. MicroRNAs (miRNAs) are 22-24 nucleotide non-coding RNAs produced by all plants and animals and some viruses that influence biological processes through post-transcriptional regulation of gene expression. Knowledge of miRNAs and their function in mite biology remains limited. Here we constructed small RNA libraries from male and female V. destructor using Illumina's small RNA-Seq platform. A total of 101,913,208 and 91,904,732 small RNA reads (>18 nucleotides) from male and female mites were analyzed using the miRDeep2 algorithm. A conservative approach predicted 306 miRNAs, 18 of which were upregulated and 13 downregulated in female V. destructor compared with males. Quantitative real-time PCR validated the expression of selected differentially-expressed female Varroa miRNAs. This dataset provides a list of potential miRNA targets involved in regulating vital Varroa biological processes and paves the way for developing strategies to target Varroa and their viruses.

Context-Dependent Viral Transgenerational Immune Priming in Honey Bees (Hymenoptera: Apidae)

Transgenerational immune priming is the process of increased resistance to infection in offspring due to parental pathogen exposure. Honey bees (Apis mellifera L. (Hymenoptera: Apidae)) are hosts to multiple pathogens, and this complex immune function could help protect against overwhelming infection. Honey bees have demonstrated transgenerational immune priming for the bacterial pathogen Paenibacillus larvae; however, evidence for viral transgenerational immune priming is lacking across insects in general. Here we test for the presence of transgenerational immune priming in honey bees with Deformed wing virus (DWV) by injecting pupae from DWV-exposed queens and measuring virus titer and immune gene expression. Our data suggest that there is evidence for viral transgenerational immune priming in honey bees, but it is highly context-dependent based on route of maternal exposure and potentially host genetics or epigenetic factors.

Standard Methods for Dissection of Varroa destructor Females

Varroa destructor (Anderson and Trueman) is known as a major pest of Apis mellifera L, especially in the Northern Hemisphere where its effects can be deleterious. As an obligate parasite, this mite relies entirely on its host to reproduce and complete its cycle. Studies focusing on isolated organs are needed to better comprehend this organism. To conduct such targeted molecular or physiological studies, the dissection of V. destructor mites is crucial as it allows the extraction of specific organs. Here, we propose a technical article showing detailed steps of females V. destructor dissection, illustrated with pictures and videos. These illustrated guidelines will represent a helpful tool to go further in V. destructor research.

Host Genotype and Tissue Type Determine DWV Infection Intensity

Varroa mite-vectored viruses such as Deformed wing virus (DWV) are of great concern for honey bee health as they can cause disease in individuals and increase colony mortality. Two genotypes of DWV (A and B) are prevalent in the United States and may have differential virulence and pathogenicity. Honey bee genetic stocks bred to resist Varroa mites also exhibit differential infection responses to the Varroa mite-vectored viruses. The goal of this project was to determine if interactions between host genotype could influence the overall infection levels and dissemination of DWV within honey bees. To do this, we injected DWV isolated from symptomatic adult bees into mite-free, newly emerged adult bees from five genetic stocks with varying levels of resistance to Varroa mites. We measured DWV-A and DWV-B dissemination among tissues chosen based on relevance to general health outcomes for 10 days. Injury from sham injections did not increase DWV-A levels but did increase DWV-B infections. DWV injection increased both DWV-A and DWV-B levels over time with significant host stock interactions. While we did not observe any differences in viral dissemination among host stocks, we found differences in virus genotype dissemination to different body parts. DWV-A exhibited the highest initial levels in heads and legs while the highest initial levels of DWV-B were found in heads and abdomens. These interactions underscore the need to evaluate viral genotype and tissue specificity in conjunction with host genotype, particularly when the host has been selected for traits relative to virus-vector and virus resistance.

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.

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.

Study of Morphological Features in Pre-Imaginal Honey Bee Impaired by Varroa destructor by Means of Computer Tomography

Simple Summary

The Varroa mite (Varroa destructor) is the most important natural pest of the honey bee, Apis mellifera, worldwide. The extent to which impairments in honey bees occur concomitantly upon infestation by this parasite greatly varies. Inter alia, the Varroa mite causes developmental disorders mediated by deformed wing virus in this host. Although there is a plethora of information regarding the consequences of this parasitism in the fully developed stage, data concerning the pre-imaginal honey bee stage inside the comb are rather scarce. In this study, morphological differences in the main body parts of the honey bee during the development stages of both intact and parasitized larvae were measured inside the comb by means of computed tomography. The images obtained reveal a visualization of the harmful effects of the Varroa mite on the pre-imaginal host. Our results demonstrate that the deformation of certain body parts was due to the presence of the parasite. Deformity, as the most conspicuous sign of infestation, is coupled with a decrease in the total-body size and abdomen size together with a disproportionate ratio of different body parts. In summary, information on the impairment of honey bee development triggered by the Varroa mite gives the opportunity to assess the damage caused by this serious pest, which occurs latently in honey bees.

Abstract

The honey bee (Apis mellifera L. 1778) is an essential element in maintaining the diversity of the biosphere and food production. One of its most important parasites is Varroa destructor, Anderson and Trueman, 2000, which plays a role in the vectoring of deformed wing virus (DWV) in honey bee colonies. Our aim was to measure the potential morphometric changes in the pre-imaginal stage of A. mellifera caused by varroosis by means of computed tomography, hence supplying evidence for the presumable role that V. destructor plays as a virus vector. Based on our results, the developmental disorders in honey bees that ensued during the pre-imaginal stages were evident. The total-body length and abdomen length of parasitized specimens were shorter than those of their intact companions. In addition, the calculated quotients of the total-body/abdomen, head/thorax, and head/abdomen in parasitized samples were significantly altered upon infestation. In our view, these phenotypical disorders can also be traced to viral infection mediated by parasitism, which was confirmed by reverse transcriptase polymerase chain reaction (RT-PCR) analysis. Capitalizing on a non-destructive method, our study reveals the deformation of the honey bee due to mite parasitism and the intermediary role this pest plays in viral infection, inside the brood cell.

Transcriptome-level assessment of the impact of deformed wing virus on honey bee larvae

Deformed wing virus (DWV) prevalence is high in honey bee (Apis mellifera) populations. The virus infects honey bees through vertical and horizontal transmission, leading to behavioural changes, wing deformity, and early mortality. To better understand the impacts of viral infection in the larval stage of honey bees, artificially reared honey bee larvae were infected with DWV (1.55 × 10¹⁰ copies/per larva). No significant mortality occurred in infected honey bee larvae, while the survival rates decreased significantly at the pupal stage. Examination of DWV replication revealed that viral replication began at 2 days post inoculation (d.p.i.), increased dramatically to 4 d.p.i., and then continuously increased in the pupal stage. To better understand the impact of DWV on the larval stage, DWV-infected and control groups were subjected to transcriptomic analysis at 4 d.p.i. Two hundred fifty-five differentially expressed genes (DEGs) (fold change ≥ 2 or ≤ -2) were identified. Of these DEGs, 168 genes were downregulated, and 87 genes were upregulated. Gene Ontology (GO) analysis showed that 141 DEGs (55.3%) were categorized into molecular functions, cellular components and biological processes. One hundred eleven genes (38 upregulated and 73 downregulated) were annotated by KO (KEGG Orthology) pathway mapping and involved metabolic pathways, biosynthesis of secondary metabolites and glycine, serine and threonine metabolism pathways. Validation of DEGs was performed, and the related gene expression levels showed a similar tendency to the DEG predictions at 4 d.p.i.; cell wall integrity and stress response component 1 (wsc1), cuticular protein and myo-inositol 2-dehydrogenase (iolG) were significantly upregulated, and small conductance calcium-activated potassium channel protein (SK) was significantly downregulated at 4 d.p.i. Related gene expression levels at different d.p.i. revealed that these DEGs were significantly regulated from the larval stage to the pupal stage, indicating the potential impacts of gene expression levels from the larval to the pupal stages. Taken together, DWV infection in the honey bee larval stage potentially influences the gene expression levels from larvae to pupae and reduces the survival rate of the pupal stage. This information emphasizes the consequences of DWV prevalence in honey bee larvae for apiculture.

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.

Replicative Deformed Wing Virus Found in the Head of Adults from Symptomatic Commercial Bumblebee (Bombus terrestris) Colonies

The deformed wing virus (DWV) is one of the most common honey bee pathogens. The virus may also be detected in other insect species, including Bombus terrestris adults from wild and managed colonies. In this study, individuals of all stages, castes, and sexes were sampled from three commercial colonies exhibiting the presence of deformed workers and analysed for the presence of DWV. Adults (deformed individuals, gynes, workers, males) had their head exscinded from the rest of the body and the two parts were analysed separately by RT-PCR. Juvenile stages (pupae, larvae, and eggs) were analysed undissected. All individuals tested positive for replicative DWV, but deformed adults showed a higher number of copies compared to asymptomatic individuals. Moreover, they showed viral infection in their heads. Sequence analysis indicated that the obtained DWV amplicons belonged to a strain isolated in the United Kingdom. Further studies are needed to characterize the specific DWV target organs in the bumblebees. The result of this study indicates the evidence of DWV infection in B. terrestris specimens that could cause wing deformities, suggesting a relationship between the deformities and the virus localization in the head. Further studies are needed to define if a specific organ could be a target in symptomatic bumblebees.

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.

The Effect of Covert and Overt Infections on Disease Dynamics in Honey-Bee Colonies

Viral diseases of honey bees are important economically and ecologically and have been widely modelled. The models reflect the fact that, in contrast to the typical case for vertebrates, invertebrates cannot acquire immunity to a viral disease, so they are of SIS or (more often) SI type. Very often, these diseases may be transmitted vertically as well as horizontally, by vectors as well as directly, and through the environment, although models do not generally reflect all these transmission mechanisms. Here, we shall consider an important additional complication the consequences of which have yet to be fully explored in a model, namely that both infected honey bees and their vectors may best be described using more than one infection class. For honey bees, we consider three infection classes. Covert infections occur when bees have the virus under control, such that they do not display symptoms of the disease, and are minimally or not at all affected by it. Acutely overtly infected bees often exhibit severe symptoms and have a greatly curtailed lifespan. Chronically overtly infected bees typically have milder symptoms and a moderately shortened lifespan. For the vector, we consider just two infection classes which are covert infected and overt infected as has been observed in deformed-wing virus (DWV) vectored by varroa mites. Using this structure, we explore the impact of spontaneous transition of both mites and bees from a covertly to an overtly infected state, which is also a novel element in modelling viral diseases of honey bees made possible by including the different infected classes. The dynamics of these diseases are unsurprisingly rather different from the dynamics of a standard SI or SIS disease. In this paper, we highlight how our compartmental structure for infection in honey bees and their vectors impact the disease dynamics observed, concentrating in particular on DWV vectored by varroa mites. If there is no spontaneous transition, then a basic reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0$$\end{document}R0 exists. We derive a condition for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0>1$$\end{document}R0>1 that reflects the complexities of the system, with components for vertical and for direct and vector-mediated horizontal transmission, using the directed graph of the next-generation matrix of the system. Such a condition has never previously been derived for a honey-bee–mite–virus system. When spontaneous transitions do occur, then \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0$$\end{document}R0 no longer exists, but we introduce a modification of the analysis that allows us to determine whether (i) the disease remains largely covert or (ii) a substantial outbreak of overt disease occurs.

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.

Detection of Lotmaria passim, Crithidia mellificae and Replicative Forms of Deformed Wing Virus and Kashmir Bee Virus in the Small Hive Beetle (Aethina tumida)

Knowledge regarding the honey bee pathogens borne by invasive bee pests remains scarce. This investigation aimed to assess the presence in Aethina tumida (small hive beetle, SHB) adults of honey bee pathogens belonging to the following groups: (i) bacteria (Paenibacillus larvae and Melissococcus plutonius), (ii) trypanosomatids (Lotmaria passim and Crithidia mellificae), and (iii) viruses (black queen cell virus, Kashmir bee virus, deformed wing virus, slow paralysis virus, sacbrood virus, Israeli acute paralysis virus, acute bee paralysis virus, chronic bee paralysis virus). Specimens were collected from free-flying colonies in Gainesville (Florida, USA) in summer 2017. The results of the molecular analysis show the presence of L. passim, C. mellificae, and replicative forms of deformed wing virus (DWV) and Kashmir bee virus (KBV). Replicative forms of KBV have not previously been reported. These results support the hypothesis of pathogen spillover between managed honey bees and the SHB, and these dynamics require further investigation.

Genetic variability and phylogenetic analysis among strains of deformed wing virus infesting honey bees and other organisms

Various viruses can infect honey bees, but deformed wing virus (DWV) is considered the most dangerous virus to them and has role in the sudden decline of bee colonies. This virus has different strains; however, there are no available studies to compare the characteristics of these strains utilizing bioinformatics. In this study, 27 strains of deformed wing virus were analyzed based on their sequences and their genetic relationships. Also, some primers were designed and tested to identify their ability to separate DWV strains. The percentages range from 28.99% to 29.63%, 22.28% to 22.78%, 15.73% to 16.28%, and 31.71% to 32.86% for nucleotides A, G, C, and T, respectively in all strains. The numbers of polymorphic sites as well as nucleotide diversity were highly similar in all strains. Statistical analyses generally showed the absence of high variations between sequences. Also, the phylogenetic tree classified strains into three groups. The network between strains of each group was established and discussed based on their geographical locations. Two groups contained strains from USA and Europe while one group contained strains from Asia. Rapid variations and mutations in the sequences of DWV were suggested. Notably, genetic studies on DWV are lacking in some geographical regions. The variations between strains detected in honey bees and other organisms were discussed. Four primers were designed and tested beside two reference primers. One of the designed primers showed the best results in binding with all DWV strains except one.

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.

Transcriptomic Responses of the Honey Bee Brain to Infection with Deformed Wing Virus

Managed colonies of European honey bees (Apis mellifera) are under threat from Varroa destructor mite infestation and infection with viruses vectored by mites. In particular, deformed wing virus (DWV) is a common viral pathogen infecting honey bees worldwide that has been shown to induce behavioral changes including precocious foraging and reduced associative learning. We investigated how DWV infection of bees affects the transcriptomic response of the brain. The transcriptomes of individual brains were analyzed using RNA-Seq after experimental infection of newly emerged adult bees with DWV. Two analytical methods were used to identify differentially expressed genes from the ~15,000 genes in the Apis mellifera genome. The 269 genes that had increased expression in DWV infected brains included genes involved in innate immunity such as antimicrobial peptides (AMPs), Ago2, and Dicer. Single bee brain NMR metabolomics methodology was developed for this work and indicates that proline is strongly elevated in DWV infected brains, consistent with the increased presence of the AMPs abaecin and apidaecin. The 1361 genes with reduced expression levels includes genes involved in cellular communication including G-protein coupled, tyrosine kinase, and ion-channel regulated signaling pathways. The number and function of the downregulated genes suggest that DWV has a major impact on neuron signaling that could explain DWV related behavioral changes.

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.

The Role of Pathogen Dynamics and Immune Gene Expression in the Survival of Feral Honey Bees

Studies of the ecoimmunology of feral organisms can provide valuable insight into how host–pathogen dynamics change as organisms transition from human-managed conditions back into the wild. Honey bees (Apis mellifera Linnaeus) offer an ideal system to investigate these questions as colonies of these social insects often escape management and establish in the wild. While managed honey bee colonies have low probability of survival in the absence of disease treatments, feral colonies commonly survive in the wild, where pathogen pressures are expected to be higher due to the absence of disease treatments. Here, we investigate the role of pathogen infections [Deformed wing virus (DWV), Black queen cell virus (BQCV), and Nosema ceranae] and immune gene expression (defensin-1, hymenoptaecin, pgrp-lc, pgrp-s2, argonaute-2, vago) in the survival of feral and managed honey bee colonies. We surveyed a total of 25 pairs of feral and managed colonies over a 2-year period (2017–2018), recorded overwintering survival, and measured pathogen levels and immune gene expression using quantitative polymerase chain reaction (qPCR). Our results showed that feral colonies had higher levels of DWV but it was variable over time compared to managed colonies. Higher pathogen levels were associated with increased immune gene expression, with feral colonies showing higher expression in five out of the six examined immune genes for at least one sampling period. Further analysis revealed that differential expression of the genes hymenoptaecin and vago increased the odds of overwintering survival in managed and feral colonies. Our results revealed that feral colonies express immune genes at higher levels in response to high pathogen burdens, providing evidence for the role of feralization in altering pathogen landscapes and host immune responses.

The novel insecticides flupyradifurone and sulfoxaflor do not act synergistically with viral pathogens in reducing honey bee (Apis mellifera) survival but sulfoxaflor modulates host immunocompetence

The decline of insect pollinators threatens global food security. A major potential cause of decline is considered to be the interaction between environmental stressors, particularly between exposure to pesticides and pathogens. To explore pesticide-pathogen interactions in an important pollinator insect, the honey bee, we used two new nicotinic acetylcholine receptor agonist insecticides (nACHRs), flupyradifurone (FPF) and sulfoxaflor (SULF), at sublethal and field-realistic doses in a fully crossed experimental design with three common viral honey bee pathogens, Black queen cell virus (BQCV) and Deformed wing virus (DWV) genotypes A and B. Through laboratory experiments in which treatments were administered singly or in combination to individual insects, we recorded harmful effects of FPF and pathogens on honey bee survival and immune gene expression. Though we found no evidence of synergistic interactions among stressors on either honey bee survival or viral load, the combined treatment SULF and DWV-B led to a synergistic upregulation of dicer-like gene expression. We conclude that common viral pathogens pose a major threat to honey bees, while co-exposure to these novel nACHR insecticides does not significantly exacerbate viral impacts on host survival in the laboratory.

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.

Selective Breeding for Low and High Varroa destructor Growth in Honey Bee (Apis mellifera) Colonies: Initial Results of Two Generations

Simple Summary

The mite Varroa destructor is considered the most damaging parasite of honey bees worldwide. Beekeepers use synthetic chemical products to control mite infestations in colonies, but the parasites soon develop resistance to them, which compromises their control. One alternative control strategy is the development of Varroa-resistant honey bees. Therefore, a breeding program was initiated to select for lower and higher rates of Varroa-population growth (LVG and HVG, respectively) and deformed wing virus (DWV) levels, which is transmitted by the mites. After two years of bidirectional selection, LVG colonies had Varroa population increases over the summer of 1.7 fold compared to 9.6 fold for HVG colonies. Additionally, HVG colonies had higher mite infestation rates in adult bees compared to LVG colonies. DWV presence and levels were higher in HVG colonies than in LVG colonies and winter mortality rates were 26% and 14% for the HVG and LVG bee types, respectively. The results of this study thus far indicate that selection for LVG may result in colonies with lower Varroa infestation rates, lower prevalence, and levels of DWV and higher colony winter survivorship. Future work will focus on determining mechanisms responsible for genetic differences and in identifying genes associated with Varroa-resistance in honey bees.

Abstract

After two years of bidirectional selection for low and high rates of Varroa destructor population growth (LVG and HVG, respectively) in honey bee (Apis mellifera) colonies in Ontario, Canada, significant differences between the two genotypes were observed. LVG colonies had V. destructor population increases over the summer of 1.7 fold compared to 9.6 fold for HVG colonies by Generation 2. Additionally, HVG colonies had significantly higher mite infestation rates in adult bees compared to LVG colonies for both selected generations. DWV prevalence and levels were significantly higher in HVG colonies than in LVG colonies in Generation 1 but not in Generation 2. Winter mortality rates of Generation 1 colonies were significantly different at 26% and 14% for the HVG and LVG genotypes, respectively. The results of this study thus far indicate that selection for LVG may result in colonies with lower V. destructor infestation rates, lower prevalence, and levels of DWV and higher colony winter survivorship. Future work will focus on determining what mechanisms are responsible for the genotypic differences, estimating genetic parameters, and molecular analyses of the genotypes to identify candidate genes associated with resistance to V. destructor and DWV that could potentially be used for marker-assisted selection.

Evaluation and comparison of the effects of three insect growth regulators on honey bee queen oviposition and egg eclosion

Honey bees (Apis mellifera) are highly valued pollinators that help to ensure national food security in the United States, but reports of heavy annual losses to managed colonies have caused concerns and prompted investigations into the causes of colony losses. One factor that can negatively affect honey bee health and survival is agrochemical exposure. Investigations into the sublethal effects of agrochemicals on important metrics of colony health such as reproduction and queen fecundity has been limited by the availability of targeted methods to study honey bee queens. This work investigates the effects of three insect growth regulators (IGR), a class of agrochemicals known to target pathways involved in insect reproduction, on honey bee queen oviposition, egg hatching, and worker hypopharyngeal development in order to quantify their effects on the fecundity of mated queens. The reported results demonstrate that none of the IGRs affected oviposition, but all three affected egg eclosion. Worker bees consuming methoxyfenozide had significantly larger hypopharyngeal glands at two weeks of age than bees not fed this compound. The results suggest that although IGRs may not exhibit direct toxic effects on adult honey bees, they can affect larval eclosion from eggs and the physiology of workers, which may contribute to colony population declines over time.

Metagenomic Approach with the NetoVIR Enrichment Protocol Reveals Virus Diversity within Ethiopian Honey Bees (Apis mellifera simensis)

Metagenomics studies have accelerated the discovery of novel or divergent viruses of the honey bee. However, most of these studies predominantly focused on RNA viruses, and many suffer from the relatively low abundance of viral nucleic acids in the samples (i.e., compared to that of the host). Here, we explored the virome of the Ethiopian honey bee, Apis mellifera simensis, using an unbiased metagenomic approach in which the next-generation sequencing step was preceded by an enrichment protocol for viral particles. Our study revealed five well-known bee viruses and 25 atypical virus species, most of which have never been found in A. mellifera before. The viruses belong to Iflaviridae, Dicistroviridae, Secoviridae, Partitiviridae, Parvoviridae, Potyviridae, and taxonomically unclassified families. Fifteen of these atypical viruses were most likely plant-specific, and the remaining ten were presumed to be insect-specific. Apis mellifera filamentous virus (AmFV) was found in one sampling site out of 10. Two samples contained high read counts of a virus similar to Diatraea saccharales densovirus (DsDNV), which is a virus that causes high mortality in the sugarcane borer. AmFV and the DsDNV-like virus were the only DNA viruses found. Three viruses that primarily infect Drosophila spp. were also discovered: La Jolla virus (LJV), Kilifi virus (KiV), and Thika virus. Our study suggests that phoretic varroa mites are involved in the transmission of LJV and KiV and that both viruses replicate in mites and adult bees. We also found an overwhelming dominance of the deformed wing virus type B variant, which fits well with the apparently harmless infestation by Varroa destructor. It was suggested that Ethiopian bees have developed tolerance against virus infections as the result of natural selection.