Deformed Wing Virus Loads Research Articles

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

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

Learning performance and GABAergic pathway link to deformed wing virus in the mushroom bodies of naturally infected honey bees.

Viral infections can be detrimental to the foraging ability of the western honey bee, Apis mellifera. The deformed wing virus (DWV) is the most common honey bee virus and has been proposed as a possible cause of learning and memory impairment. However, evidence for this phenomenon so far has come from artificially infected bees, while less is known about the implications of natural infections with the virus. Using the proboscis extension reflex (PER), we uncovered no significant association between a simple associative learning task and natural DWV load. However, when assessed through a reversal associative learning assay, bees with higher DWV load performed better in the reversal learning phase. DWV is able to replicate in the honey bee mushroom bodies, where the GABAergic signalling pathway has an antagonistic effect on associative learning but is crucial for reversal learning. Hence, we assessed the pattern of expression of several GABA-related genes in bees with different learning responses. Intriguingly, mushroom body expression of selected genes was positively correlated with DWV load, but only for bees with good reversal learning performance. We hypothesise that DWV might improve olfactory learning performance by enhancing the GABAergic inhibition of responses to unrewarded stimuli, which is consistent with the behavioural patterns that we observed. However, at higher disease burdens, which might be induced by an artificial infection or by a severe, natural Varroa infestation, this DWV-associated increase in GABA signalling could impair associative learning as previously reported by other studies.

Parasite and virus dynamics in the honeybee Apis mellifera unicolor on a tropical island recently invaded by Varroa destructor

In La Réunion, the established honeybee subspecies Apis mellifera unicolor, an endemic subspecies of African lineage, is facing considerable challenges. Since the introduction of the Varroa destructor mite in 2017 high colony losses have been recorded. We investigated the dynamics of V. destructor and two viruses, the Deformed Wing Virus (DWV), known to be transmitted by the mite, and the Chronic Bee Paralysis Virus (CBPV), in A. m. unicolor. Colonies from two apiaries located at 300 and 900 m a.s.l were monitored twice for one year without any acaricide treatment. The brood area, V. destructor infestation rates, DWV and CBPV prevalence and load were recorded monthly. A. m. unicolor maintained brood rearing throughout the year. Varroa destructor infestation resulted in high colony mortality (up to 85 %) and high phoretic mite rates (up to 52 mites per hundred bees). The establishment of DWV in colonies occurred after that of V. destructor and the mite infestation rate had a significant effect on the virus prevalence and load. CBPV appeared only transiently throughout the surveys. The data showed that, in tropical colonies with permanent brood rearing, V. destructor and DWV can reach high levels, but are still subject to seasonal variations that appear to be influenced by environmental conditions. This suggests that beekeeping practices could be adapted by favouring sites and periods for transhumance or acaricide treatment.

Sampling methodology and storage management effect on deformed wing virus (DWV) quantification using real-time PCR

Deformed wing virus (DWV) is one of the most important viruses affecting Apis mellifera. A valid method for estimating the viral load at the honey bee colony level is of utmost importance. In this study, different samplings of bees were performed, in order to determine the appropriate sample size and collection point in the hive. A similar trend was observed at both sampling points, and subsequently, the DWV load was quantified in samples collected from the entire hive area. Whether 10 or 50 bees were sampled from each frame, the viral titre showed a similar dispersion and the lowest coefficient of variance was recorded in 10 bees/frame samples. According to our results, each one of the tested sampling methods can be used to quantify the viral load at the colony level, eliminating the need to sample 50 bees per frame due to the large number of bees. The effect of post-collection sample management, specifically the combination of storage temperature and duration, on RNA extraction yield, purity and virus quantification results was also investigated. The DWV titre determination was successful for all groups. It seems that improper storage affects the results of the viral quantification, however, significant differences were not observed among the treatments. Undoubtedly, immediate sample preservation in extremely low temperatures is the most appropriate method for virus analyses.

Transmission of deformed wing virus between Varroa destructor foundresses, mite offspring and infested honey bees

BackgroundVarroa destructor is the major ectoparasite of the western honey bee (Apis mellifera). Through both its parasitic life-cycle and its role as a vector of viral pathogens, it can cause major damage to honey bee colonies. The deformed wing virus (DWV) is the most common virus transmitted by this ectoparasite, and the mite is correlated to increased viral prevalence and viral loads in infested colonies. DWV variants A and B (DWV-A and DWV-B, respectively) are the two major DWV variants, and they differ both in their virulence and transmission dynamics.MethodsWe studied the transmission of DWV between bees, parasitic mites and their offspring by quantifying DWV loads in bees and mites collected in in vitro and in situ environments. In vitro, we artificially transmitted DWV-A to mites and quantified both DWV-A and DWV-B in mites and bees. In situ, we measured the natural presence of DWV-B in bees, mites and mites’ offspring.ResultsBee and mite viral loads were correlated, and mites carrying both variants were associated with higher mortality of the infected host. Mite infestation increased the DWV-B loads and decreased the DWV-A loads in our laboratory conditions. In situ, viral quantification in the mite offspring showed that, after an initially non-infected egg stage, the DWV-B loads were more closely correlated with the foundress (mother) mites than with the bee hosts.ConclusionsThe association between mites and DWV-B was highlighted in this study. The parasitic history of a mite directly impacts its DWV infection potential during the rest of its life-cycle (in terms of variant and viral loads). Regarding the mite’s progeny, we hypothesize that the route of contamination is likely through the feeding site rather than by vertical transmission, although further studies are needed to confirm this hypothesis.Graphical

The Beneficial Effect of Pollen on Varroa Infested Bees Depends on Its Influence on Behavioral Maturation Genes.

Honey bees collect nectar and pollen to fulfill their nutritional demands. In particular, pollen can influence longevity, the development of hypopharyngeal glands, and immune-competence of bees. Pollen can also mitigate the deleterious effects caused by the parasitic mite Varroa destructor and related deformed wing virus (DWV) infections. It has been shown that V. destructor accelerates the physiological and behavioral maturation of honey bees by influencing the interaction between two core physiological factors, Vitellogenin and juvenile hormone. In this study, we test the hypothesis that the beneficial effects of pollen on Varroa-infested bees are related to the hormonal control underpinning behavioral maturation. By analyzing the expression of genes associated to behavioral maturation in pollen-fed mite-infested bees, we show that treatment with pollen increases the lifespan of mite-infested bees by reversing the faster maturation induced by the parasite at the gene expression level. As expected, from the different immune-competence of nurse and forager bees, the lifespan extension triggered by pollen is also correlated with a positive influence of antimicrobial peptide gene expression and DWV load, further reinforcing the beneficial effect of pollen. This study lay the groundwork for future analyses of the underlying evolutionary processes and applications to improve bee health.

Pupal cannibalism by worker honey bees contributes to the spread of deformed wing virus

Transmission routes impact pathogen virulence and genetics, therefore comprehensive knowledge of these routes and their contribution to pathogen circulation is essential for understanding host–pathogen interactions and designing control strategies. Deformed wing virus (DWV), a principal viral pathogen of honey bees associated with increased honey bee mortality and colony losses, became highly virulent with the spread of its vector, the ectoparasitic mite Varroa destructor. Reproduction of Varroa mites occurs in capped brood cells and mite-infested pupae from these cells usually have high levels of DWV. The removal of mite-infested pupae by worker bees, Varroa Sensitive Hygiene (VSH), leads to cannibalization of pupae with high DWV loads, thereby offering an alternative route for virus transmission. We used genetically tagged DWV to investigate virus transmission to and between worker bees following pupal cannibalisation under experimental conditions. We demonstrated that cannibalization of DWV-infected pupae resulted in high levels of this virus in worker bees and that the acquired virus was then transmitted between bees via trophallaxis, allowing circulation of Varroa-vectored DWV variants without the mites. Despite the known benefits of hygienic behaviour, it is possible that higher levels of VSH activity may result in increased transmission of DWV via cannibalism and trophallaxis.

Effects of Deformed Wing Virus Infection on Expressions of Immune- and Apoptosis-Related Genes in Western Honeybees (Apis mellifera).

Simple SummaryDeformed wing virus (DWV) is a widespread viral pathogen in western honeybees. This work investigated the effects of DWV on honeybee and host defense by analyzing the immune and apoptosis gene expressions. The results revealed that the mortality rates and deformed wing symptoms in newly emerged adult bees increased at high initial concentration of DWV injection. Moreover, abnormality of wings in DWV injected bees did not correlate with the virus levels nor the immune and apoptosis responses. After DWV injection, the immune response was activated in pupae and newly emerged adult bees. However, apoptosis genes were stimulated only in newly emerged adult bees but suppressed in honeybee pupae. The apoptosis gene suppression during the initial phase of infection may promote the DWV survival and parasitism in the honeybee host.Honeybees are globally threatened by several pathogens, especially deformed wing virus (DWV), as the presence of DWV in western honeybees is indicative of colony loss. The high mortality rate is further exacerbated by the lack of effective treatment, and therefore understanding the immune and apoptosis responses could pave an avenue for the treatment method. In this study, DWV was directly injected into the white-eyed pupae stage of western honeybees (Apis mellifera). The DWV loads and selected gene responses were monitored using the real-time PCR technique. The results showed that honeybee pupae that were injected with the highest concentration of viral loads showed a significantly higher mortality rate than the control groups. Deformed wings could be observed in newly emerged adult bees when the infected bees harbored high levels of viral loads. However, the numbers of viral loads in both normal and crippled wing groups were not significantly different. DWV-injected honeybee pupae with 104 and 107 copy numbers per bee groups showed similar viral loads after 48 h until newly emerged adult bees. Levels of gene expression including immune genes (defensin, abaecin, and hymenoptaecin) and apoptosis genes (buffy, p53, Apaf1, caspase3-like, caspase8-like, and caspase9-like) were analyzed after DWV infection. The expressions of immune and apoptosis genes were significantly different in infected bees compared to those of the control groups. In the pupae stage, the immune genes were activated by injecting DWV (defensin and hymenoptaecin) or Escherichia coli (defensin, abaecin, and hymenoptaecin), a positive control. On the contrary, the expression of apoptosis-related genes (buffy, caspase3-like, caspase8-like, and caspase9-like genes) was suppressed at 96 h post-infection. In DWV-infected newly emerged adult bees, abaecin, hymenoptaecin, Apaf1, and caspase8-like genes were upregulated. However, these genes were not significantly different between the normal and crippled wing bees. Our results suggested that DWV could activate the humoral immunity in honeybees and that honeybee hosts may be able to protect themselves from the virus infection through immune responses. Apoptosis gene expressions were upregulated in newly emerged adult bees by the virus, however, they were downregulated during the initial phase of viral infection.

Success of Selection in Terms of Hygienic Behavior in Struggle Against DWV and Varroa

Abstract Many studies have reported sudden and large-scale losses in honey bee colonies. There is no precise explanation for this happening, but it may be caused by the interaction of pathogenic viruses and ectoparasite Varroa destructor. Deformed-wing virus (DWV) is common in Apis mellifera L. and has been implicated in worldwide Varroa-associated colony losses. Hygienic behavior may be one of the ways that limit V. destructor infestation level and DWV load. The purpose of this research was to determine whether the selection of higher hygienic behavior is effective against V. destructor and DWV. Workers and pupae of A. m. anatoliaca were collected from colonies showing four generations of hygienic behaviodr, and non-hygienics were used as a control. Loads of DWV were investigated by one-step RT-qPCR, and the bottom-board method was used to allow mites count. The prevalence of DWV in pupae were higher in hygienic colonies (69.0%, average viral load 234.38) than control colonies (60.8% average viral load 937.70). However, the prevalence of worker bees infected with DWV was lower in hygienic colonies (41.3%, average viral load 181627.69) than in control's (66.0%, average viral load 241982.35). Although the averages of the Varroa counts were significantly different for both hygienic (28.92) and control colonies (108.90), we could not find any correlation between DWV loads of both workers and pupae in this study. When all these results are evaluated, the selection performed in terms of hygienic behavior can said to have been said successful against DWV.

Interactions Between Thiamethoxam and Deformed Wing Virus Can Drastically Impair Flight Behavior of Honey Bees.

Exposure to multiple stress factors is believed to contribute to honey bee colony decline. However, little is known about how co-exposure to stress factors can alter the survival and behavior of free-living honey bees in colony conditions. We therefore studied the potential interaction between a neonicotinoid pesticide, thiamethoxam, and a highly prevalent honey bee pathogen, Deformed wing virus (DWV). For this purpose, tagged bees were exposed to DWV by feeding or injection, and/or to field-relevant doses of thiamethoxam, then left in colonies equipped with optical bee counters to monitor flight activity. DWV loads and the expression of immune genes were quantified. A reduction in vitellogenin expression level was observed in DWV-injected bees and was associated with precocious onset of foraging. Combined exposure to DWV and thiamethoxam did not result in higher DWV loads compared to bees only exposed to DWV, but induced precocious foraging, increased the risk of not returning to the hive after the first flight, and decreased survival when compared to single stress exposures. We therefore provided the first evidence for deleterious interactions between DWV and thiamethoxam in natural conditions.

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.

Validation of quantitative real-time RT-PCR assays for the detection of six honeybee viruses

Acute bee paralysis virus (ABPV), Black queen cell virus (BQCV), Chronic bee paralysis virus (CBPV), Deformed wing virus (DWV), Sacbrood virus (SBV) and Varroa destructor virus 1 (VDV1) are the six main honeybee viruses reported in Europe. We assessed the accuracy (trueness and precision) of reverse transcriptase quantitative TaqMan® PCR methods (RT-qPCR) for quantifying ABPV, BQCV, DWV, VDV1 and SBV loads. Once the systematic bias in quantitative results had been corrected (overestimation in ABPV and BQCV quantification and underestimation in that of SBV and VDV1), measurements were taken to determine the viral load ranges for which quantification uncertainty was below ± 1 log10 equivalent of genome copies per bee (hereafter reported as genome copies/bee). The accuracy range of RT-qPCR was found to be between 6.4 and 10.4 log10 genome copies/bee for ABPV, between 3.0 and 10.0 log10 genome copies/bee for BQCV, between 2.4 and 10.4 log10 genome copies/bee for DWV and between 3.4 and 10.4 log10 genome copies/bee for SBV. Outside these ranges, the results’ uncertainty is higher. VDV1 RT-qPCR accuracy was outside accuracy limits for all viral loads. Using these RT-qPCR methods, we quantified viral loads in naturally-infected honeybees. The viral load distribution and clinical signs reported with the honeybee samples allowed us to define a threshold that could be used to differentiate between covert and overt infections. These methods will be useful in diagnosing the main viral infections impairing honeybee health.

Deformed Wing Virus Infection in Honey Bees (Apis mellifera L.).

Deformed wing virus (DWV) is a single-stranded RNA virus of honey bees (Apis mellifera L.) transmitted by the parasitic mite Varroa destructor. Although DWV represents a major threat to honey bee health worldwide, the pathological basis of DWV infection is not well documented. The objective of this study was to investigate clinicopathological and histological aspects of natural DWV infection in honey bee workers. Emergence of worker honey bees was observed in 5 colonies that were clinically affected with DWV and the newly emerged bees were collected for histopathology. DWV-affected bees were 2 times slower to emerge and had 30% higher mortality compared to clinically normal bees. Hypopharyngeal glands in bees with DWV were hypoplastic, with fewer intracytoplasmic secretory vesicles; cells affected by apoptosis were observed more frequently. Mandibular glands were hypoplastic and were lined by cuboidal epithelium in severely affected bees compared to tall columnar epithelium in nonaffected bees. The DWV load was on average 1.7 × 106 times higher (P < .001) in the severely affected workers compared to aged-matched sister honey bee workers that were not affected by deformed wing disease based on gross examination. Thus, DWV infection is associated with prolonged emergence, increased mortality during emergence, and hypoplasia of hypopharyngeal and mandibular glands in newly emerged worker honey bees in addition to previously reported deformed wing abnormalities.

Occurrence of deformed wing virus variants in the stingless bee Melipona subnitida and honey bee Apis mellifera populations in Brazil.

The global spread of the parasitic Varroa mite has introduced a new bee to the bee horizontal transmission route for several RNA viruses that bypasses existing barriers in honey bees. From among these viruses, deformed wing virus (DWV) is now among the most widespread insect pathogens in the world. Brazilian stingless bees are a diverse group often managed in close proximity to honey bees. Therefore, we investigated the prevalence and load of DWV in 21 stingless bee (Melipona subnitida) and 26 honey bee (Apis mellifera) colonies from Brazil. DWV was detected in all colonies with DWV-A and DWV-C dominating in M. subnitida, while DWV-A dominated in A. mellifera. Average total viral loads per bee were 8.8E+07 and 6.3E+07 in M. subnitida and A. mellifera, respectively, which are much lower than DWV levels (>1E+10) found in honey bees in the northern hemisphere. In colonies introduced 30 years ago to the remote island of Fernando de Noronha, the DWV load was low (<1E+03) in honey bees but we detected higher loads (1.6E+08) in all M. subnitida colonies on the island. This may suggest that minimal, if any, viral transmission of DWV from stingless bees to honey bees has occurred on this island. Furthermore, the ubiquitous presence of the DWV-C variant in M. subnitida colonies, and its rarity in A. mellifera, may again suggest that limited viral exchange between these two species is occurring.

A real-time PCR method for quantification of the total and major variant strains of the deformed wing virus.

European honey bees (Apis mellifera) are critically important to global food production by virtue of their pollination services but are severely threatened by deformed wing virus (DWV) especially in the presence of the external parasite Varroa destructor. DWV exists as many viral strains with the two major variants (DWV-A and DWV-B) varying in virulence. A single plasmid standard was constructed containing three sections for the specific determination of DWV-A (VP2 capsid region), DWV-B (IRES) and a conserved region suitable for total DWV (helicase region). The assays were confirmed as specific and discriminatory with limits of detections of 25, 25 and 50 genome equivalents for DWV-A, DWV-B and total-DWV, respectively. The methods were successfully tested on Apis mellifera and V. destructor samples with varying DWV profiles. The new method determined a more accurate total DWV titre in samples with substantial DWV-B than the method currently described in the COLOSS Beebook. The proposed assays could be utilized for the screening of large quantities of bee material for both a total DWV load overview along with more detailed investigations into DWV-A and DWV-B profiles.

Colony adaptive response to simulated heat waves and consequences at the individual level in honeybees (Apis mellifera)

Since climate change is expected to bring more severe and frequent extreme weather events such as heat waves, assessing the physiological and behavioural sensitivity of organisms to temperature becomes a priority. We therefore investigated the responses of honeybees, an important insect pollinator, to simulated heat waves (SHW). Honeybees are known to maintain strict brood thermoregulation, but the consequences at the colony and individual levels remain poorly understood. For the first time, we quantified and modelled colony real-time activity and found a 70% increase in foraging activity with SHW, which was likely due to the recruitment of previously inactive bees. Pollen and nectar foraging was not impacted, but an increase in water foragers was observed at the expense of empty bees. Contrary to individual energetic resources, vitellogenin levels increased with SHW, probably to protect bees against oxidative stress. Finally, though immune functions were not altered, we observed a significant decrease in deformed wing virus loads with SHW. In conclusion, we demonstrated that honeybees could remarkably adapt to heat waves without a cost at the individual level and on resource flow. However, the recruitment of backup foraging forces might be costly by lowering the colony buffering capacity against additional environmental pressures.

Oldest Varroa tolerant honey bee population provides insight into the origins of the global decline of honey bees

The ecto-parasitic mite Varroa destructor has transformed the previously inconsequential Deformed Wing Virus (DWV) into the most important honey bee viral pathogen responsible for the death of millions of colonies worldwide. Naturally, DWV persists as a low level covert infection transmitted between nest-mates. It has long been speculated that Varroa via immunosuppression of the bees, activate a covert infection into an overt one. Here we show that despite Varroa feeding on a population of 20–40 colonies for over 30 years on the remote island of Fernando de Noronha, Brazil no such activation has occurred and DWV loads have remained at borderline levels of detection. This supports the alternative theory that for a new vector borne viral transmission cycle to start, an outbreak of an overt infection must first occur within the host. Therefore, we predict that this honey bee population is a ticking time-bomb, protected by its isolated position and small population size. This unique association between mite and bee persists due to the evolution of low Varroa reproduction rates. So the population is not adapted to tolerate Varroa and DWV, rather the viral quasispecies has simply not yet evolved the necessary mutations to produce a virulent variant.

A Comparison of Deformed Wing Virus in Deformed and Asymptomatic Honey Bees.

Deformed wing virus (DWV) in association with Varroa destructor is currently attributed to being responsible for colony collapse in the western honey bee (Apis mellifera). The appearance of deformed individuals within an infested colony has long been associated with colony losses. However, it is unknown why only a fraction of DWV positive bees develop deformed wings. This study concerns two small studies comparing deformed and non-deformed bees. In Brazil, asymptomatic bees (no wing deformity) that had been parasitised by Varroa as pupae had higher DWV loads than non-parasitised bees. However, we found no greater bilateral asymmetry in wing morphology due to DWV titres or parasitisation. As expected, using RT-qPCR, deformed bees were found to contain the highest viral loads. In a separate study, next generation sequencing (NGS) was applied to compare the entire DWV genomes from paired symptomatic and asymptomatic bees from three colonies on two different Hawaiian islands. This revealed no consistent differences between DWV genomes from deformed or asymptomatic bees, with the greatest variation seen between locations, not phenotypes. All samples, except one, were dominated by DWV type A. This small-scale study suggests that there is no unique genetic variant associated with wing deformity; but that many DWV variants have the potential to cause deformity.

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.

Varroa destructor (Mesostigmata: Varroidae) Parasitism and Climate Differentially Influence the Prevalence, Levels, and Overt Infections of Deformed Wing Virus in Honey Bees (Hymenoptera: Apidae).

The prevalence and loads of deformed wing virus (DWV) between honey bee (Apis mellifera L.) colonies from a tropical and a temperate environment were compared. The interaction between these environments and the mite Varroa destructor in relation to DWV prevalence, levels, and overt infections, was also analyzed. V. destructor rates were determined, and samples of mites, adult bees, brood parasitized with varroa mites and brood not infested by mites were analyzed. DWV was detected in 100% of the mites and its prevalence and loads in honey bees were significantly higher in colonies from the temperate climate than in colonies from the tropical climate. Significant interactions were found between climate and type of sample, with the highest levels of DWV found in varroa-parasitized brood from temperate climate colonies. Additionally, overt infections were observed only in the temperate climate. Varroa parasitism and DWV loads in bees from colonies with overt infections were significantly higher than in bees from colonies with covert infections. These results suggest that interactions between climate, V. destructor, and possibly other factors, may play a significant role in the prevalence and levels of DWV in honey bee colonies, as well as in the development of overt infections. Several hypotheses are discussed to explain these results.