Transcriptome analysis of the synganglion from the honey bee mite, Varroa destructor and RNAi knockdown of neural peptide targets

Vector-virus interaction affects viral loads and co-occurrence

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

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.

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.

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.

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: 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.

Identification of Immune Response to Sacbrood Virus Infection in Apis cerana Under Natural Condition

Chinese sacbrood virus (CSBV) is a serious threat to eastern honeybees (Apis cerana), especially larvae. However, the pathological mechanism of this deadly disease remains unclear. Here, we employed mRNA and small RNA (sRNA) transcriptome approach to investigate the microRNAs (miRNAs) and small interfering RNAs (siRNAs) expression changes of A. cerana larvae infected with CSBV under natural condition. We found that serine proteases involved in immune response were down-regulated, while the expression of siRNAs targeted to serine proteases were up-regulated. In addition, CSBV infection also affected the expression of larvae cuticle proteins such as larval cuticle proteins A1A and A3A, resulting in increased susceptibility to CSBV infection. Together, our results provide insights into sRNAs that they are likely to be involved in regulating honeybee immune response.

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.

Varroa destructor: how does it harm Apis mellifera honey bees and what can be done about it?

Since its migration from the Asian honey bee (Apis cerana) to the European honey bee (Apis mellifera), the ectoparasitic mite Varroa destructor has emerged as a major issue for beekeeping worldwide. Due to a short history of coevolution, the host–parasite relationship between A. mellifera and V. destructor is unbalanced, with honey bees suffering infestation effects at the individual, colony and population levels. Several control solutions have been developed to tackle the colony and production losses due to Varroa, but the burden caused by the mite in combination with other biotic and abiotic factors continues to increase, weakening the beekeeping industry. In this synthetic review, we highlight the main advances made between 2015 and 2020 on V. destructor biology and its impact on the health of the honey bee, A. mellifera. We also describe the main control solutions that are currently available to fight the mite and place a special focus on new methodological developments, which point to integrated pest management strategies for the control of Varroa in honey bee colonies.

Veterinary Diagnostic Approach of Common Virus Diseases in Adult Honeybees

Veterinarians are educated in prevention, diagnosis and treatment of diseases in various vertebrate species. As they are familiar with multifactorial health problems in single animals as well as in herd health management, their knowledge and skills can be beneficial for the beekeepers and honeybee health. However, in education and in practice, honeybees are not a common species for most veterinarians and the typical veterinary diagnostic methods such as blood sampling or auscultation are not applicable to the superorganism honeybee. Honeybee colonies may be affected by various biotic and abiotic factors. Among the infectious agents, RNA-viruses build the largest group, causing covert and overt infections in honeybee colonies which may lead to colony losses. Veterinarians could and should play a more substantial role in managing honeybee health—not limited to cases of notifiable diseases and official hygiene controls. This review discusses the veterinary diagnostic approach to adult bee examination with a special focus on diagnosis and differential diagnosis of the common virus diseases Acute Bee Paralysis Virus (ABPV)-Kashmir Bee Virus (KBV)-Israeli Acute Paralysis Virus (IAPV)-Complex, Chronic Bee Paralysis Virus (CBPV) and Deformed Wing Virus (DWV), as well as coinfections like Varroa spp. and Nosema spp.

Bee Viruses: Routes of Infection in Hymenoptera

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

Varroa destructor: A Complex Parasite, Crippling Honey Bees Worldwide

The parasitic mite, Varroa destructor, has shaken the beekeeping and pollination industries since its spread from its native host, the Asian honey bee (Apis cerana), to the naïve European honey bee (Apis mellifera) used commercially for pollination and honey production around the globe. Varroa is the greatest threat to honey bee health. Worrying observations include increasing acaricide resistance in the varroa population and sinking economic treatment thresholds, suggesting that the mites or their vectored viruses are becoming more virulent. Highly infested weak colonies facilitate mite dispersal and disease transmission to stronger and healthier colonies. Here, we review recent developments in the biology, pathology, and management of varroa, and integrate older knowledge that is less well known.

Green Bees: Reverse Genetic Analysis of Deformed Wing Virus Transmission, Replication, and Tropism

Environmental and agricultural pollination services by honey bees, Apis mellifera, and honey production are compromised by high levels of annual colony losses globally. The majority are associated with disease caused by deformed wing virus (DWV), a positive-strand RNA virus, exacerbated by the ectoparasitic mite Varroa destructor. To improve honey bee health, a better understanding of virus transmission and pathogenesis is needed which requires the development of tools to study virus replication, transmission, and localisation. We report the use of reverse genetic (RG) systems for the predominant genetically distinct variants of DWV to address these questions. All RG-recovered viruses replicate within 24 h post-inoculation of pupae and could recapitulate the characteristic symptoms of DWV disease upon eclosion. Larvae were significantly less susceptible but could be infected orally and subsequently developed disease. Using genetically tagged RG DWV and an in vitro Varroa feeding system, we demonstrate virus replication in the mite by accumulation of tagged negative-strand viral replication intermediates. We additionally apply a modified DWV genome expressing a fluorescent reporter protein for direct in vivo observation of virus distribution in injected pupae or fed larvae. Using this, we demonstrate extensive sites of virus replication in a range of pupal tissues and organs and in the nascent wing buds in larvae fed high levels of virus, indicative of a direct association between virus replication and pathogenesis. These studies provide insights into virus replication kinetics, tropism, transmission, and pathogenesis, and produce new tools to help develop the understanding needed to control DWV-mediated colony losses.

Diversity and Global Distribution of Viruses of the Western Honey Bee, Apis mellifera

In the past centuries, viruses have benefited from globalization to spread across the globe, infecting new host species and populations. A growing number of viruses have been documented in the western honey bee, Apis mellifera. Several of these contribute significantly to honey bee colony losses. This review synthetizes the knowledge of the diversity and distribution of honey-bee-infecting viruses, including recent data from high-throughput sequencing (HTS). After presenting the diversity of viruses and their corresponding symptoms, we surveyed the scientific literature for the prevalence of these pathogens across the globe. The geographical distribution shows that the most prevalent viruses (deformed wing virus, sacbrood virus, black queen cell virus and acute paralysis complex) are also the most widely distributed. We discuss the ecological drivers that influence the distribution of these pathogens in worldwide honey bee populations. Besides the natural transmission routes and the resulting temporal dynamics, global trade contributes to their dissemination. As recent evidence shows that these viruses are often multihost pathogens, their spread is a risk for both the beekeeping industry and the pollination services provided by managed and wild pollinators.

The honeybee (Apis mellifera) developmental state shapes the genetic composition of the deformed wing virus-A quasispecies during serial transmission

The main biological threat to the western honeybee (Apis mellifera) is the parasitic mite Varroa destructor, largely because it vectors lethal epidemics of honeybee viruses that, in the absence of this mite, are relatively innocuous. The severe pathology is a direct consequence of excessive virus titres caused by this novel transmission route. However, little is known about how the virus adapts genetically during transmission and whether this influences the pathology. Here, we show that upon injection into honeybee pupae, the deformed wing virus type-A (DWV-A) quasispecies undergoes a rapid, extensive expansion of its sequence space, followed by strong negative selection towards a uniform, common shape by the time the pupae have completed their development, with no difference between symptomatic and asymptomatic adults in either DWV titre or genetic composition. This suggests that the physiological and molecular environment during pupal development has a strong, conservative influence on shaping the DWV-A quasispecies in emerging adults. There was furthermore no evidence of any progressive adaptation of the DWV-A quasispecies to serial intra-abdominal injection, simulating mite transmission, despite the generation of ample variation immediately following each transmission, suggesting that the virus either had already adapted to transmission by injection, or was unaffected by it.

Dynamic evolution in the key honey bee pathogen deformed wing virus: Novel insights into virulence and competition using reverse genetics

The impacts of invertebrate RNA virus population dynamics on virulence and infection outcomes are poorly understood. Deformed wing virus (DWV), the main viral pathogen of honey bees, negatively impacts bee health, which can lead to colony death. Despite previous reports on the reduction of DWV diversity following the arrival of the parasitic mite Varroa destructor, the key DWV vector, we found high genetic diversity of DWV in infested United States honey bee colonies. Phylogenetic analysis showed that divergent US DWV genotypes are of monophyletic origin and were likely generated as a result of diversification after a genetic bottleneck. To investigate the population dynamics of this divergent DWV, we designed a series of novel infectious cDNA clones corresponding to coexisting DWV genotypes, thereby devising a reverse-genetics system for an invertebrate RNA virus quasispecies. Equal replication rates were observed for all clone-derived DWV variants in single infections. Surprisingly, individual clones replicated to the same high levels as their mixtures and even the parental highly diverse natural DWV population, suggesting that complementation between genotypes was not required to replicate to high levels. Mixed clone–derived infections showed a lack of strong competitive exclusion, suggesting that the DWV genotypes were adapted to coexist. Mutational and recombination events were observed across clone progeny, providing new insights into the forces that drive and constrain virus diversification. Accordingly, our results suggest that Varroa influences DWV dynamics by causing an initial selective sweep, which is followed by virus diversification fueled by negative frequency-dependent selection for new genotypes. We suggest that this selection might reflect the ability of rare lineages to evade host defenses, specifically antiviral RNA interference (RNAi). In support of this hypothesis, we show that RNAi induced against one DWV strain is less effective against an alternate strain from the same population.

Deformed wing virus, a key pathogen of honey bees, shows rapid diversification after genetic bottlenecks; a novel reverse-genetic system provides insights into the forces that shape virus diversity, suggesting that virus quasi-species diversification may be driven by selection of genotypes capable of evading host RNAi defences.

Deformed wing virus type A, a major honey bee pathogen, is vectored by the mite Varroa destructor in a non-propagative manner

Honey bees, the primary managed insect pollinator, suffer considerable losses due to Deformed wing virus (DWV), an RNA virus vectored by the mite Varroa destructor. Mite vectoring has resulted in the emergence of virulent DWV variants. The basis for such changes in DWV is poorly understood. Most importantly, it remains unclear whether replication of DWV occurs in the mite. In this study, we exposed Varroa mites to DWV type A via feeding on artificially infected honey bees. A significant, 357-fold increase in DWV load was observed in these mites after 2 days. However, after 8 additional days of passage on honey bee pupae with low viral loads, the DWV load dropped by 29-fold. This decrease significantly reduced the mites’ ability to transmit DWV to honey bees. Notably, negative-strand DWV RNA, which could indicate viral replication, was detected only in mites collected from pupae with high DWV levels but not in the passaged mites. We also found that Varroa mites contain honey bee mRNAs, consistent with the acquisition of honey bee cells which would additionally contain DWV replication complexes with negative-strand DWV RNA. We propose that transmission of DWV type A by Varroa mites occurs in a non-propagative manner.

Deformed Wing Virus in Honeybees and Other Insects

Deformed wing virus (DWV) has become the most well-known, widespread, and intensively studied insect pathogen in the world. Although DWV was previously present in honeybee populations, the arrival and global spread of a new vector, the ectoparasitic mite Varroa destructor, has dramatically altered DWV epidemiology. DWV is now the most prevalent virus in honeybees, with a minimum average of 55% of colonies/apiaries infected across 32 countries. Additionally, DWV has been detected in 65 arthropod species spanning eight insect orders and three orders of Arachnida. Here, we describe the significant progress that has been made in elucidating the capsid structure of the virus, understanding its ever-expanding host range, and tracking the constantly evolving DWV genome and formation of recombinants. The construction of molecular clones, working with DWV in cell lines, and the development of immunohistochemistry methods will all help the community to move forward. Identifying the tissues in which DWV variants are replicating and understanding the impact of DWV in non-honeybee hosts are major new goals.

Varroa chemosensory proteins: some are conserved across Arthropoda but others are arachnid specific

The tight synchronization between the life cycle of the obligatory parasitic mite Varroa destructor (Varroa) and its host, the honeybee, is mediated by honeybee chemical stimuli. These stimuli are mainly perceived by a pit organ located on the distal part of the mite's foreleg. In the present study, we searched for Varroa chemosensory molecular components by comparing transcriptomic and proteomic profiles between forelegs from different physiological stages, and rear legs. In general, a comparative transcriptomic analysis showed a clear separation of the expression profiles between the rear legs and the three groups of forelegs (phoretic, reproductive and tray-collected mites). Most of the differentially expressed transcripts and proteins in the mite's foreleg were previously uncharacterized. Using a conserved domain approach, we identified 45 transcripts with known chemosensory domains belonging to seven chemosensory protein families, of which 14 were significantly upregulated in the mite's forelegs when compared to rear legs. These are soluble and membrane bound proteins, including the somewhat ignored receptors of degenerin/epithelial Na+ channels and transient receptor potentials. Phylogenetic clustering and expression profiles of the putative chemosensory proteins suggest their role in chemosensation and shed light on the evolution of these proteins in Chelicerata.

Genetic Diversity of Deformed Wing Virus From Apis mellifera carnica (Hymenoptera: Apidae) and Varroa Mite (Mesostigmata: Varroidae)

Deformed wing virus (DWV) is one of the most widespread viruses that infect honey bee colonies. The route of infection is directly through contaminated food, feces, and air, or indirectly through the varroa mite, which acts as a vector. Positive DWV samples were obtained from Carniolan honey bee (Apis mellifera carnica) colonies and of varroa mites from the whole territory of Slovenia during a survey between 2007 and 2014. Nucleotide sequences of 471 nucleotides for the L protein gene and 573 nucleotides for the helicase gene were compared. High genetic diversity was observed among these Slovenian Carniolan honey bee DWV field samples, as well as with almost all the strains previously found in other countries in Europe. Phylogenetic analyses in two regions of the viral genome show that several of the DWV strains obtained from honey bees and varroa are genetically very closely related, confirming the important role of varroa in the transmission of DWV. Identification of closely related sequences also confirmed that the same strains of DWV have been successfully transmitted between various honey bee colonies and apiaries. It has also been established that simultaneous infection, in one apiary, of honey bees with two or more different strains of DWV is quite frequent. This is phylogenetic study that compares honey bee and varroa DWV strains from Carniolan honeybees.

Genes important for survival or reproduction in Varroa destructor identified by RNAi

The Varroa mite, (Varroa destructor), is the worst threat to honey bee health worldwide. To explore the possibility of using RNA interference to control this pest, we determined the effects of knocking down various genes on Varroa mite survival and reproduction. Double-stranded RNA (dsRNA) of six candidate genes (Da, Pros26S, RpL8, RpL11, RpP0 and RpS13) were synthesized and each injected into Varroa mites, then mite survival and reproduction were assessed. Injection of dsRNA for Da (Daughterless) and Pros26S (Gene for proteasome 26S subunit adenosine triphosphatase) caused a significant reduction in mite survival, with 3.57% ± 1.94% and 30.03% ± 11.43% mites surviving at 72 h post-injection (hpi), respectively. Control mites injected with green fluorescent protein (GFP)-dsRNA showed survival rates of 81.95% ± 5.03% and 82.36 ± 2.81%, respectively. Injections of dsRNA for four other genes (RpL8, RpL11, RpP0 and RpS13) did not affect survival significantly, enabling us to assess their effect on Varroa mite reproduction. The number of female offspring per mite was significantly reduced for mites injected with dsRNA of each of these four genes compared to their GFP-dsRNA controls. Knockdown of the target genes was verified by real-time polymerase chain reaction for two genes important for reproduction (RpL8, RpL11) and one gene important for survival (Pros26S). In conclusion, through RNA interference, we have discovered two genes important for mite survival and four genes important for mite reproduction. These genes could be explored as possible targets for the control of Varroa destructor in the future.

In vivo evolution of viral virulence: switching of deformed wing virus between hosts results in virulence changes and sequence shifts

The health of the Western honey bee is threatened by a global epidemic of deformed wing virus (DWV) infections driven by the ectoparasitic mite Varroa destructor acting as mechanical and biological virus vector. Three different variants of DWV, DWV-A, -B and -C exist. Virulence differences between these variants and their relation to V. destructor are still controversially discussed. We performed laboratory experiments to analyze the virulence of DWV directly isolated from crippled bees (DWV_P0 ) or after one additional passage in bee pupae (DWV_P1 ). We demonstrated that DWV_P0 was more virulent than DWV_P1 for pupae, when pupal mortality was taken as virulence marker, and for adult bees, when neurotropism and cognitive impairment were taken as virulence markers. Phylogenetic analysis supported that DWV exists as quasispecies and showed that DWV_P0 clustered with DWV-B and DWV_P1 with DWV-A when the phylogeny was based on the master sequences of the RNA-dependent RNA polymerase but not so when it was based on the VP3 region master sequences. We propose that switching of DWV between the bee and the mite host is accompanied by changes in viral sequence, tissue tropism and virulence and that the RNA-dependent RNA polymerase is involved in determining host range and virulence.

Beyond insects: current status and achievements of RNA interference in mite pests and future perspectives

Mites comprise a group of key agricultural pests on a wide range of crops. They cause harm through feeding on the plant and transferring dangerous pathogens, and the rapid evolution of pesticide resistance in mites highlights the need for novel control methods. Currently, RNA interference (RNAi) shows great potential for insect pest control. Here, we review the literature regarding RNAi in mite pests. We discuss different target genes and RNAi efficiency in various mite species, a promising Varroa control program using RNAi, the synergy of RNAi with plant defense mechanisms and microorganisms, and current understanding of systemic movement of double-stranded RNA (dsRNA). On the basis of this evidence, we can conclude that there is clear potential for application of RNAi-based mite control, but further research on several aspects of RNAi in mites is needed, including: (i) the factors influencing RNAi efficiency, (ii) the mechanism of environmental RNAi and cross-kingdom dsRNA trafficking, (iii) the mechanism of possible systemic and parental RNAi, and (iv) non-target effects, specifically in predatory mites, which should be considered during RNAi target selection. © 2018 Society of Chemical Industry.

RNA interference in the cat flea, Ctenocephalides felis: Approaches for sustained gene knockdown and evidence of involvement of Dicer-2 and Argonaute2

Effective RNA interference (RNAi) methods have been developed in many pest species, enabling exploration of gene function. Until now RNAi had not been attempted in the cat flea, Ctenocephalides felis, although the development of RNAi approaches would open up potential avenues for control of this important pest. This study aimed to establish if an RNAi response occurs in adult C. felis upon exposure to double-stranded RNA (dsRNA), which administration methods for dsRNA delivery could bring about effective gene knockdown and to investigate dynamics of any RNAi response. Knockdown of 80% of GSTσ was achieved by intrahaemoceolic microinjection of dsGSTσ but this invasive technique was associated with relatively high mortality rates. Immersing C. felis in dsGSTσ or dsDicer-2 overnight resulted in 65% knockdown of GSTσ or 60% of Dicer-2, respectively, and the degree of knockdown was not improved by increasing the dsRNA concentration in the bathing solution. Unexpectedly, the greatest degree of knockdown was achieved with the continuous administration of dsRNA in whole blood via a membrane feeding system, resulting in 96% knockdown of GSTσ within 2 days and sustained up to, at least, 7 days. Thus, unlike in many other species, the gut nucleases do not impair the RNAi response to ingested dsRNA in C. felis. A modest, but significant, upregulation of Dicer-2 and Argonaute2 was detectable 3 h after exposure to exogenous dsRNA, implicating the short-interfering RNA pathway. To our knowledge this study represents the first demonstration of experimentally induced RNAi in the cat flea as well as giving insight into how the gene knockdown response progresses.

Diverse Factors Affecting Efficiency of RNAi in Honey Bee Viruses

Infection and transmission of honey bee viruses pose a serious threat to the pollination services of crops and wild plants, which plays a vital role in agricultural economy and ecology. RNA interference (RNAi) is an effective defense mechanism against commonly occurring viral infections of animals and plants. However, recent studies indicate that the effects of RNAi on the honey bee can induce additional impacts and might not always be effective in suppressing the virus. Moreover, the RNAi responses differed in relation to the developmental stage of the insect and the target tissue used, even though the same method of delivery was used. These results indicate that further analysis and field experiments should be performed to characterize the varying effectiveness of RNAi-based methods for treating honey bee viral infections. In this review, we provide an overview of the current knowledge and the recent progress in RNAi-based anti-viral treatments for honey bees, focusing in particular highlight the role of the dsRNA-delivery method used and its effect on RNAi efficiency and demonstrate the potential practical value of this tool for controlling the virus. We conclude studying the gene function and disease control of honey bee by RNAi technology requires a complex consideration from physiology, genetics to environment.

Transcriptome profiling of the honeybee parasite Varroa destructor provides new biological insights into the mite adult life cycle

Background

The parasite Varroa destructor represents a significant threat to honeybee colonies. Indeed, development of Varroa infestation within colonies, if left untreated, often leads to the death of the colony. Although its impact on bees has been extensively studied, less is known about its biology and the functional processes governing its adult life cycle and adaptation to its host. We therefore developed a full life cycle transcriptomic catalogue in adult Varroa females and included pairwise comparisons with males, artificially-reared and non-reproducing females (10 life cycle stages and conditions in total).

Results

Extensive remodeling of the Varroa transcriptome was observed, with an upregulation of energetic and chitin metabolic processes during the initial and final phases of the life cycle (e.g. phoretic and post-oviposition stages), whereas during reproductive stages in brood cells genes showing functions related to transcriptional regulation were overexpressed. Several neurotransmitter and neuropeptide receptors involved in behavioural regulation, as well as active compounds of salivary glands, were also expressed at a higher level outside the reproductive stages. No difference was detected between artificially-reared phoretic females and their counterparts in colonies, or between females who failed to reproduce and females who successfully reproduced, indicating that phoretic individuals can be reared outside host colonies without impacting their physiology and that mechanisms underlying reproductive failure occur before oogenesis.

Conclusions

We discuss how these new findings reveal the remarkable adaptation of Varroa to its host biology and notably to the switch from living on adults to reproducing in sealed brood cells. By spanning the entire adult life cycle, our work captures the dynamic changes in the parasite gene expression and serves as a unique resource for deciphering Varroa biology and identifying new targets for mite control.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4668-z) contains supplementary material, which is available to authorized users.

Chemosensing of honeybee parasite, Varroa destructor: Transcriptomic analysis

Chemosensing is a primary sense in nature, however little is known about its mechanism in Chelicerata. As a model organism we used the mite Varroa destructor, a key parasite of honeybees. Here we describe a transcriptomic analysis of two physiological stages for the Varroa foreleg, the site of primary olfactory organ. The transcriptomic analysis revealed transcripts of chemosensory related genes belonging to several groups. These include Niemann-Pick disease protein, type C2 (NPC2), gustatory receptors (GRs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and odorant binding proteins (OBP). However, no insect odorant receptors (ORs) and odorant co-receptors (ORcos) were found. In addition, we identified a homolog of the most ancient IR co-receptor, IR25a, in Varroa as well as in other members of Acari. High expression of this transcript in the mite's forelegs, while not detectable in the other pairs of legs, suggests a function for this IR25a-like in Varroa chemosensing.

A Varroa destructor protein atlas reveals molecular underpinnings of developmental transitions and sexual differentiation

Varroa destructor is the most economically damaging honey bee pest, weakening colonies by simultaneously parasitizing bees and transmitting harmful viruses. Despite these impacts on honey bee health, surprisingly little is known about its fundamental molecular biology. Here, we present a Varroa protein atlas crossing all major developmental stages (egg, protonymph, deutonymph, and adult) for both male and female mites as a web-based interactive tool (http://foster.nce.ubc.ca/varroa/index.html). We used intensity-based label-free quantitation to find 1,433 differentially expressed proteins across developmental stages. Enzymes for processing carbohydrates and amino acids were among many of these differences as well as proteins involved in cuticle formation. Lipid transport involving vitellogenin was the most significantly enriched biological process in the foundress (reproductive female) and young mites. In addition, we found that 101 proteins were sexually regulated and functional enrichment analysis suggests that chromatin remodeling may be a key feature of sex determination. In a proteogenomic effort, we identified 519 protein-coding regions, 301 of which were supported by two or more peptides and 169 of which were differentially expressed. Overall, this work provides a first-of-its-kind interrogation of the patterns of protein expression that govern the Varroa life cycle and the tools we have developed will support further research on this threatening honey bee pest.

A feeding protocol for delivery of agents to assess development in Varroa mites

A novel feeding protocol for delivery of bio-active agents to Varroa mites was developed by providing mites with honey bee larva hemolymph supplemented with cultured insect cells and selected materials delivered on a fibrous cotton substrate. Mites were starved, fed on treated hemolymph to deliver selected agents and then returned to bee larvae. Transcript levels of two reference genes, actin and glyceraldehyde 3-phosphate dehydrogenase (GAPDH), as well as for nine selected genes involved in reproductive processes showed that the starvation and feeding protocol periods did not pose a high level of stress to the mites as transcript levels remained comparable between phoretic mites and those completing the protocol. The feeding protocol was used to deliver molecules such as hormone analogs or plasmids. Mites fed with Tebufenozide, an ecdysone analog, had higher transcript levels of shade than untreated or solvent treated mites. In order to extend this feeding protocol, cultured insect cells were incorporated to a final ratio of 1 part cells and 2 parts hemolymph. Although supplementation with Bombyx mori Bm5 cells increased the amount of hemolymph consumed per mite, there was a significant decrease in the percentage of mites that fed and survived. On the other hand, Drosophila melanogaster S2 cells reduced significantly the percentage of mites that fed and survived as well as the amount of hemolymph consumed. The feeding protocol provides a dynamic platform with which to challenge the Varroa mite to establish efficacy of control agents for this devastating honey bee pest.

A Toolbox for Quantitative Gene Expression in Varroa destructor: RNA Degradation in Field Samples and Systematic Analysis of Reference Gene Stability

Varroa destructor is the major pest of Apis mellifera and contributes to the global honey bee health crisis threatening food security. Developing new control strategies to combat Varroa will require the application of molecular biology, including gene expression studies by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Both high quality RNA samples and suitable stable internal reference genes are required for accurate gene expression studies. In this study, ten candidate genes (succinate dehydrogenase (SDHA), NADH dehydrogenase (NADH), large ribsosmal subunit, TATA-binding protein, glyceraldehyde-3-phosphate dehydrogenase, 18S rRNA (18S), heat-shock protein 90 (HSP90), cyclophilin, α-tubulin, actin), were evaluated for their suitability as normalization genes using the geNorm, Normfinder, BestKeeper, and comparative ΔCq algorithims. Our study proposes the use of no more than two of the four most stable reference genes (NADH, 18S, SDHA and HSP90) in Varroa gene expression studies. These four genes remain stable in phoretic and reproductive stage Varroa and are unaffected by Deformed wing virus load. When used for determining changes in vitellogenin gene expression, the signal-to-noise ratio (SNR) for the relatively unstable genes actin and α-tubulin was much lower than for the stable gene combinations (NADH + HSP90 +18S; NADH + HSP90; or NADH). Using both electropherograms and RT-qPCR for short and long amplicons as quality controls, we demonstrate that high quality RNA can be recovered from Varroa up to 10 days later stored at ambient temperature if collected into RNAlater and provided the body is pierced. This protocol allows the exchange of Varroa samples between international collaborators and field sample collectors without requiring frozen collection or shipping. Our results make important contributions to gene expression studies in Varroa by proposing a validated sampling protocol to obtain high quality Varroa RNA and the validation of suitable reference genes for expression studies in this globally important pest.