Acaricide Treatment Affects Viral Dynamics in Varroa destructor-Infested Honey Bee Colonies via both Host Physiology and Mite Control

ABSTRACT Honey bee (Apis mellifera) colonies are declining, and a number of stressors have been identified that affect, alone or in combination, the health of honey bees. The ectoparasitic mite Varroa destructor, honey bee viruses that are often closely associated with the mite, and pesticides used to control the mite population form a complex system of stressors that may affect honey bee health in different ways. During an acaricide treatment using Apistan (plastic strips coated with tau-fluvalinate), we analyzed the infection dynamics of deformed wing virus (DWV), sacbrood virus (SBV), and black queen cell virus (BQCV) in adult bees, mite-infested pupae, their associated Varroa mites, and uninfested pupae, comparing these to similar samples from untreated control colonies. Titers of DWV increased initially with the onset of the acaricide application and then slightly decreased progressively coinciding with the removal of the Varroa mite infestation. This initial increase in DWV titers suggests a physiological effect of tau-fluvalinate on the host's susceptibility to viral infection. DWV titers in adult bees and uninfested pupae remained higher in treated colonies than in untreated colonies. The titers of SBV and BQCV did not show any direct relationship with mite infestation and showed a variety of possible effects of the acaricide treatment. The results indicate that other factors besides Varroa mite infestation may be important to the development and maintenance of damaging DWV titers in colonies. Possible biochemical explanations for the observed synergistic effects between tau-fluvalinate and virus infections are discussed.

R ecent declines in honey bee abundance in the United States and Europe have led to an increased gap between supply and demand in pollination services, causing pressure both economically in agricultural crop production (18) and ecologically on natural plant biodiversity (4). Today, no clear explanation for these declines exists. However, a variety of stressors affecting honey bee health, such as viruses (13), microsporidia (23), and pesticides (51), have been identified as potential causes of such colony losses.
The ectoparasitic mite Varroa destructor is a well-known pest and has become the most significant economic threat to apiculture on a global scale (30). By feeding on hemolymph, the mite is a significant stressor on honey bee health, causing a variety of physical and physiological effects for both individual bees and the colony (30,41). In addition, the indirect effects of Varroa, primarily caused by the viral infections vectored by the mite, are the most devastating and can lead to severe disease and mortality at the individual or colony level (2,5,10,15,33,40,44). These viruses usually persist as covert infections in the colony without the vectoring capacity of Varroa (37). Besides Varroa-mediated transmission, viruses can also be transmitted horizontally within the colony through trophallaxis, contact, feces, and salivary gland secretions and vertically from infected queens to their progeny (10,12,14,54). There are over 18 identified and characterized honey bee viruses, many of which are suggested to be associated with mite infestation to various degrees (2,9,12,33,34,40,44,49). Most notably, deformed wing virus (DWV) has become the most prevalent virus associated with Varroa mite infestation (15,20), causing wing deformities in developing pupae that result in flightless adults that die shortly after emerging (2,35).
Exponential mite population growth can lead to increased virus transmission, causing overt viral infections ultimately resulting in colony mortality within 2 to 3 years (5,33). In order to prevent colony mortality, management strategies are imple-mented by beekeepers to reduce the mite population within the colony and thereby limit the transmission opportunities for potentially lethal virus infections. A common Varroa control treatment is Apistan (plastic strips coated with tau-fluvalinate), a waxsoluble, nonvolatile, potent pyrethroid acaricide with 98% to 100% efficacy against Varroa mite infestation (41,52) except where fluvalinate resistance is present (25,52). Although taufluvalinate had previously been considered to have low toxicity for bees (38), recent studies have shown it to be a stressor on colony health, particularly affecting the development and performance of queens (22). In addition, tau-fluvalinate has been found to act synergistically with sublethal doses of other Varroa-targeted pesticides administered to honey bee colonies (28). Recently, it was demonstrated that the neonicotinoid pesticide imidacloprid interacted synergistically with Nosema sp., a honey bee pathogen that is otherwise relatively benign in infectivity (1). If pesticides are an additional stress to honey bee colonies and can influence host susceptibility to the infectivity of certain pathogens, what effect might be seen on the susceptibility of honey bees to virus infections when using a pesticide against their virus vector, the Varroa mite?
The relationship between Varroa and honey bee viruses has been well documented in laboratory studies and field surveys (3,6,9,34,35,37,49,50), and the outcomes of virus infections and Varroa mite population growth rates and transmission dynamics have been modeled (33,47). It has also been demonstrated that remedial Varroa control treatment can reverse virus titer accumulation in colonies (3,34). However, these studies looked at longterm virus titers after the completion of the acaricide treatment. The present study investigates the dynamics of virus infections during the acaricide treatment, aiming to identify any effects of tau-fluvalinate (Apistan) treatment on virus titers in adult bees, Varroa-infested and uninfested pupae, and Varroa mites.

MATERIALS AND METHODS
Experimental design. In July 2008, nine single-brood-chamber honey bee colonies were established from splits containing equal proportions of adult bees, brood, honey, pollen, and mite infestations in Uppsala, Sweden. Colonies with high mite infestations were chosen under the assumption that higher mite levels would cause higher viral titers and, therefore, that a decrease in titers due to changing mite infestation would be more easily recognized. On 8 August 2008, six of the colonies were treated with Apistan (plastic strips coated with tau-fluvalinate), while the remaining three colonies served as an untreated control group.
Data and sample collection. An initial set of measurements and samples were taken from the colonies prior to the acaricide treatment, while subsequent sampling continued approximately once a week for the duration of the treatment, giving a total of six data sets per colony. Population estimations of the adult bees and brood were made using the Liebefeld estimation method (27). The phoretic Varroa mite infestation rates were determined by washing samples of approximately 200 bees with soapy water to dislodge the mites (16). The brood Varroa mite infestation rates were determined by opening 100 randomly selected pupal cells (17). Samples of adult bees, mite-infested pupae, the mites infesting these pupae, and uninfested pupae were collected and stored at Ϫ20°C until subsequent virus analysis (11).
RNA extraction. Bulk samples of 30 adult bees were placed with 5 ml nuclease-free water in a plastic mesh bag, flash frozen with liquid nitrogen, and homogenized with a pestle. Bulk samples of 5 infested pupae were separated from their mites and homogenized in 2.5 ml RLT buffer (Qiagen) with a micropestle. The corresponding 5 mites were pooled and homogenized mechanically in 0.5 ml RLT buffer with glass beads using the FastPrep system (Thermo Electron Corporation Sweden AB). Bulk samples of 5 uninfested pupae were homogenized similarly to the infested pupae. Total RNA was extracted by a Qiacube automated extraction robot (Qiagen) from 100 l of each crude extract using the RNEasy protocol for plant tissues (Qiagen) and stored as two 50-l aliquots at Ϫ80°C.
RT-qPCR assays. To determine which viruses were present in the colonies at the start of the experiment, adult bees, infested pupae, and uninfested pupae collected from the first sampling date were screened by reverse transcription-quantitative PCR (RT-qPCR) for 10 honey bee viruses: acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), DWV, Israeli acute paralysis virus (IAPV), Kashmir bee virus (KBV), sacbrood virus (SBV), slow bee paralysis virus (SBPV), Varroa destructor virus 1 (VDV-1), and Varroa destructor macula-like virus (VdMLV) (see Table S1 in the supplemental material). Only those viruses detected in the initial screen, DWV, BQCV, and SBV, were analyzed in the remaining samples. DWV, as mentioned above, is widespread and has a well-characterized association with Varroa mites, while BQCV and SBV are both virulent viruses at the individual level affecting honey bee brood, but less virulent at the colony level and with unclear associations with Varroa mites. Each sample was also assayed for Apis mellifera or Varroa ␤-actin mRNA, using intron-spanning primers (see Table S1 in the supplemental material). ␤-Actin is a common internal reference RNA (32), used here to normalize the RT-qPCR data for differences between samples in the quantity and quality of the RNA. All assays had excellent linearity (r 2 ) (see Table S1 in the supplemental material) between log 10 [template] and cycle number, covering the full range of experimental titers, as well as highly acceptable PCR efficiencies (E) (see Table S1 in the supplemental material).
The amounts of DWV, BQCV, SBV, and ␤-actin in each sample were determined using the Bio-Rad iScript One-Step RT-qPCR Kit with SYBR Green as the detection chemistry, 96-well optical qPCR plates, and the Bio-Rad Chromo4 thermocycler. The assays were performed in 20-l volumes containing 3.0 l RNA, 0.2 l of each primer, and 0.4 l iScript. All primers were designed to have the same optimal annealing temperature so that all four assays for each sample can be run concurrently in the same plate, together with five positive controls and one negative H 2 O control for each assay. The positive controls were obtained from 10-fold serial dilutions of a purified PCR product of known concentration and covered 6 orders of magnitude. The following reaction conditions were used: 10 min at 50°C for cDNA synthesis plus 5 min at 95°C for reverse transcription inactivation and Taq polymerase activation, followed by 35 cycles of 10 s at 95°C for denaturation and 30 s at 58°C for primer annealing, extension, and data collection. The amplification reaction was followed by a melting curve analysis to determine the specificity of the amplification products by incubating them for 60 s at 95°C and 60 s at 55°C and then reading the fluorescence at 0.5°C increments from 55°C to 95°C. Nonspecific amplification products melt and lose fluorescence at a distinctly lower melting temperature (T m ) than the specific "true" reaction products (8).
RT-qPCR data conversion. The RT-qPCR data were first screened for the presence of specific target PCR product, as determined by the melting curve analyses. The quantification cycle (Cq) values (8) of all confirmed target amplifications were calculated with Bio-Rad Opticon Monitor software version 3.1 following a reaction baseline subtraction using the Global Minimum Trend option and with the fluorescence threshold set uniformly at 0.05 for all plates. All Cq values coincided with the logarithmic phase of the amplifications. For each target RNA, the Cq values of the external dilution standards of all RT-qPCR runs were pooled into a single linear regression analysis of Cq value onto log 10 [template]. These regression equations were then used to estimate the absolute amounts of virus and ␤-actin RNA in each reaction. These were subsequently converted to estimated amounts per bee through the different reaction and extraction dilution factors. The regression slopes were also used to calculate the amplification efficiencies (E) of the different assays as follows: E assay ϭ 10 Ϫ1/slope (39).
Data transformations. The virus titers were exponentially distributed, spanning between 4 and 7 orders of magnitude, depending on the virus. This broad titer distribution is a logical consequence of the largely exponential nature of virus replication and proliferation (7). A log transformation was therefore used to render the data suitable for parametric analyses (45). Since it is not possible to log transform zero values (i.e., instances of no virus amplification), a constant value is usually added to the data to allow the inclusion of zero values in the analyses (45). For these experiments, this constant value was determined by assuming a hypothetical Cq value of 41 for the virus assays and converting these to theoretical titers as indicated above.
Data normalization. The data set was partitioned into 12 subsets (3 viruses ϫ 4 sample types) for data normalization and statistical analysis. This was a necessary consequence of the separate ␤-actin titer adjustments for the different sample types and the distinct data and error distributions of the three viruses ( Fig. 1). For each subset, the virus titer and mite infestation rate data for each colony were normalized to the average pretreatment value for all colonies to avoid statistically significant effects purely due to natural, preexisting differences between the colonies.
Statistical analysis. The effect of acaricide treatment on the normalized virus titers in each of the four sample types (adults, infested pupae, uninfested pupae, and Varroa mites) was analyzed by mixed-effects models as implemented by SAS 9.1 for Windows and included the normalized adult and pupal mite infestation rates as independent explanatory variables. The models included a linear repeated-measures factor to account for the bias inherent in sampling the same colonies repeatedly over time, with the covariance structure selected based on the Aikaike's information criteria (AIC) value (31). Assumptions of normality were verified by analysis of residuals and equality of variance (31).
Correlation analysis was used to assess the relationships across the 12 data subsets, i.e., between the three viruses in each sample type, between the different sample types for each virus, or between the virus titers and the mite infestation rates in pupae or adults, using the original, prenormalized data. The correlations involving the mite infestation rates included only data from the pretreated and untreated colonies to avoid biasing the natural global correlation between virus titers and the mite infestation rates by the rapid artificial removal of the mite population from the treated colonies.

Mite and bee population development.
Although the colonies were randomly allocated to the treated and untreated groups, the initial mite infestation rates were higher in the treated group than in the untreated group for both the brood and adult infestation rates (Fig. 2). The acaricide treatment significantly reduced the mite infestation rates (df in parentheses) in the brood [F (10,35) ϭ 15.04; P Ͻ 0.0001] and on adult bees [F (10,35) ϭ 8.27; P Ͻ 0.0001] but did not eliminate mite infestation completely (Fig. 2). A very low adult mite infestation rate persisted in all the treated colonies after the second week of treatment, even though mites were no longer detectable in the brood in all but one of the acaricide-treated colonies after the third week of treatment (Fig. 2). The numbers of adult bees and brood in the untreated colonies were significantly reduced during the study to nearly half their initial numbers compared to the treated colonies, whose populations remained relatively constant [adult bees, F (10,33) ϭ 7.09, P Ͻ 0.0001; brood, F (10,33) ϭ 2.22, P ϭ 0.0413]. Clinical symptoms of DWV were observed with increasing frequency in the untreated colonies, while in the treated colonies, clinical DWV symptoms decreased. The presence of clinically deformed bees, which have very short life expectancy, was likely the cause of the reduced adult bee population in the untreated group, also considering that the mite infestation rates of adult bees significantly affected the numbers of adult bees [F (1, 33) ϭ 5.38; P ϭ 0.0267]. All untreated colonies died before or during the following winter, while all treated colonies survived.
Assays. There were clear differences between the four sample types (infested pupae, uninfested pupae, adults, and mites) in the total estimated levels of ␤-actin mRNA per bee or mite, covering a 6-fold range, although there were no differences between treated and untreated colonies in this regard ( Fig. 1) (t tests, P Ͼ 0.05). The amount of virus in each sample was therefore adjusted to the average ␤-actin levels for that particular sample type, effectively using ␤-actin as a molecular marker for the quantity and quality of the RNA sample. These values were then multiplied by the different experimental dilution factors to arrive at the estimated number of virus genome equivalents per sample. A negative result worth noting is the apparent absence of all viruses other than DWV, BQCV, and SBV, including the absence of VDV-1, a close relative of DWV that is specifically associated with V. destructor (15).
Virus titer analyses. The statistical effects of the acaricide treatment on the normalized virus titers over the 6-week sampling period, and how much of this is due to the variable effect of treatment with time or can be explained by the mite infestation rates, are summarized in Table 1. The most significant effects of acaricide treatment are seen in titers in adult bees and mites for DWV, in titers in infested pupae for BQCV, and in titers in both infested and uninfested pupae for SBV. The mite infestation rates (both adult and pupal) affected only the DWV titers of adult bees (Table 1).
Relative to pretreatment values, an initial increase in the average DWV titers in the treated colonies was observed in miteinfested and uninfested pupae and adult bees, which was not observed in the untreated colonies. The untreated colonies experienced a decrease in DWV titers relative to the initial titers (Fig. 3, DWV) in mite-infested pupae and adult bees that did not rise above the titers in the treated colonies until the 4th week of treatment. The DWV titers in uninfested pupae in the untreated colonies were never higher than in the treated colonies. The most pronounced treatment effect is the reduction of DWV titers in mite samples from the treated group by the 4th and 5th weeks of treatment (Fig. 3). Although this is based on a single colony, since mites were undetectable in all other treated colonies by week 4, this reduction with treatment is a true observation for that colony and not an artifact of using group means. The significance in the mixed-effects model analysis for this drastic reduction in the mite DWV titer is reflected in both the treatment and treatment-date effects.
BQCV and SBV titers were consistently higher in adult bees of the acaricide-treated colonies than in those of the untreated control colonies (Fig. 3). However, relative to the pretreatment titers, the fluctuations in both BQCV and SBV titers in the different samples with increasing treatment duration appear random rather than progressive (Fig. 3), although they were large enough to render most of the corresponding treatment-date interactions significant (Table 1).
In all samples, the DWV titers were several orders of magnitude higher than the BQCV and SBV titers. While DWV titers in infested pupae were 2 to 3 orders of magnitude higher than in uninfested pupae, there was no difference between infested and uninfested pupae for the BQCV or SBV titers (Fig. 3). DWV titers for adult bees were slightly lower than for the mite-infested pupae but considerably higher than for the uninfested pupae (Fig. 3, DWV). The BQCV and SBV titers in adult bees were both slightly higher than for either infested or uninfested pupae (Fig. 3, BQCV and SBV).
Correlation analyses. The relationships between the 12 data subsets were assessed through correlation analyses of the original prenormalized data. The titers of the three viruses within each sample type were generally not correlated, i.e., the viruses appear to act independently of each other ( Table 2). The only exception is in the infested pupae, where SBV varies positively with both DWV and BQCV. The DWV titers in infested pupae and their mites were highly correlated, but not the BQCV or the SBV titers, suggesting that during pupal infestation, Varroa has a more active relationship with DWV than with BQCV or SBV (Table 2). This observation is indirectly supported by the significant correlation between infested and uninfested pupae for the BQCV and SBV titers, which appear to be unaffected by mite infestation, but not for the DWV titers, where mite infestation has a large impact ( Fig. 2 and Table 2).
Of special interest are the relationships between the virus titers in adult bees and those in the pupae and mites, since the pupae will become part of the adult population, while the adults can influence the pupal titers indirectly through oral transmission to the larvae. For DWV, all these correlations were positive and highly significant, but this is true for only one of the SBV correlations (between adults and infested pupae) and not at all for any of the BQCV correlations (Table 2). There was excellent correspondence between the two mite infestation rate estimates. Of the viruses, only the DWV titers are correlated with the mite infestation rates in infested pupae, their mites, and adult bees. These correlations for the pupal infestation rates are superior to those for the adult infestation rates (Table 2), which is largely due to the more even distribution of the pupal infestation rate data across the correlation. There was no relationship between the mite infestation rates and the BQCV or SBV titers in any of the samples ( Table 2).

Viral dynamics during acaricide treatment.
Given the strong influence of Varroa infestation and Varroa-vectored transmission on DWV titers in honeybee colonies (6,33,42,48,53), the expectation before this study was that removal of the mites using an acaricide treatment should result in a progressive reduction of DWV titers throughout the colony with a trend similar to the mite infestation rate, as was previously shown (34). On the contrary, our results showed a significant increase in DWV titers coinciding with the onset of the tau-fluvalinate treatment application in adult honey bees, mite-infested pupae, and uninfested pupae. A slight progressive decrease in DWV titers in adult bees and pupae was seen in treated colonies after this initial increase. However, the DWV titers in adult bees and uninfested pupae never fell below those in the untreated colonies during the treatment period. The DWV titers in mite-infested pupae eventually dropped below the titers in the untreated colonies, but only in the final week of the treatment. These results illustrate the short-term colony level changes in virus titers during the actual acaricide treatment rather than the long-term dynamics of virus titers in either treated or untreated colonies.
The initial increase of DWV titers in adult bees and pupae of the acaricide-treated colonies coincides with the most potent chemical effect of the treatment and could therefore potentially be a consequence of debilitating direct effects of tau-fluvalinate on honey bee physiology and/or immune system responses, causing increased host susceptibility to DWV infection. The subsequent decrease in DWV titers after the initial increase is at least partly due to the removal of Varroa-mediated transmission, as implied by the strong correlation between DWV titers and Varroa infestation rates in both adults and pupae. However, it may also be an effect of the adaptation of the host immune system to the presence of the tau-fluvalinate.
Although the primary effect of the acaricide treatment is the removal of mites from the colony, tau-fluvalinate is a potent chemical pesticide that can have distinct physiological effects on exposed bees (22). Synergistic interactions between chemicals and infectious pathogens causing reduced fitness in honey bee colonies have been demonstrated before (1). At the cellular and biochemical levels, fluvalinate interferes with the voltage-gated sodium ion transport channels (36,41). These channels regulate the osmotic potential of the cell and are therefore also frequently a target for virus infections that use osmotic pressure to burst the cell to release newly formed virus particles (26,29,46). This could be a possible explanation for the initial increase in DWV titers seen in the bees. The sharp reduction in DWV titers in those mites surviving the fluvalinate treatment could have a similar origin. A variable proportion of Varroa mites are themselves also hosts for DWV (53), and such mites have highly elevated DWV titers compared to those mites that do not replicate DWV (21). DWVreplicating mites may therefore suffer enhanced mortality due to synergistic interactions between fluvalinate and DWV at the sodium transport channels so that only mites that do not replicate DWV, with drastically reduced DWV titers, have a chance of surviving the fluvalinate treatment. These are speculative arguments, but certainly plausible. Mechanisms behind the possible direct effects of the tau-fluvalinate treatment on host immunity and pathogen virulence in both bees and mites, as well as those of other chemical treatments administered to honey bee colonies, need further investigation.
With respect to BQCV and SBV infection, the acaricide treatment seemed to have a random influence rather than a progressive one. Although for these viruses there are also slight increase in titers following the initial treatment application in most sample types, it is much less pronounced than for DWV. In adult bees, BQCV and SBV titers remained relatively constant in colonies with the tau-fluvalinate acaricide treatment, while those in the untreated colonies decreased, something also seen for the SBV titers in infested pupae. The fluctuations in BQCV and SBV titers over time and between the treated and untreated colonies were large enough to produce statistically significant treatment-date effects for infested and uninfested pupae in the mixed-effects model analyses, but it is unclear how these can be interpreted, since they did not show a clear progressive increase or decrease throughout the experiment.
Varroa-virus relationships. The strong positive correlation between DWV titers in infested pupae, adult bees, and Varroa mites and the colony pupal mite infestation; the statistically significant explanatory effect of colony mite infestation rates on DWV titers in adult bees; and the elevated DWV titers in infested pupae compared to uninfested pupae are all observations that complement numerous prior reports demonstrating the active re-lationship between DWV and Varroa (2,6,15,21,53). However, the unexpected observation in these experiments was that DWV titers remained high in our treated colonies even after the removal of the mites. Notably, DWV titers also remained high in the uninfested pupae throughout the experiment, even though these pupae did not have any contact with Varroa mites. These observa- tions suggest that high-titer DWV infections may be able to persist for a considerable time after the removal of Varroa through alternative transmission routes, such as oral or vertical transmission. This supports observations by Highfield et al. (24) that winter colony mortality due to DWV infections is partly independent of Varroa mite infestation levels. The second unexpected observation was a significant initial increase in the DWV titer immediately following acaricide application in both adults and pupae, which subsequently decreased progressively. Highfield et al. (24) also implied the existence of an additional, unknown factor interacting with DWV that contributed to colony mortality. Perhaps the results of the present study offer a hint as to the nature of this unknown factor, considering that the colonies of the Highfield et al. study (24) were treated with acaricides before the winter and before their experiment, began. The incidental observation that clinical symptoms, i.e., deformed bees, disappeared in the treated colonies while those in the untreated colonies increased, even though there was very little net change in DWV titer between treated and untreated colonies over the study period, reflects the very particular circumstances that lead to wing deformities. Although deformed bees usually have higher DWV titers than nondeformed bees, it is not the titers alone that result in deformities but rather how and when the virus is transmitted. While DWV is directly responsible for these wing deformities (35), they only appear in natural circumstances through Varroa-mediated DWV transmission during the pupal stage (6,21,53). Varroa-mediated DWV transmission among adult bees also results in high DWV titers but does not result in deformities (10,19,35,48). In addition, deformed bees are ex-pelled more rapidly from stronger colonies and may therefore also simply be more easily detected in the untreated colonies than in the treated ones. During an acaricide treatment, the adult bee population is increasingly derived from uninfested pupae (with lower DWV titers) as the pupal mite infestation rate drops, resulting in a rapid decrease in deformed emerging bees, even with a similar overall colony level DWV titer. The gradual decrease in DWV titers following mite treatment is due to a progressive turnover of the bee population with increasingly healthier bees rather than an instant cleansing with the removal of the mite vector alone. Just as it takes more than one season to build up to lethal DWV infection levels with Varroa mite infestation, it takes more than 5 weeks to bring them back down again through the reverse process of eliminating the mite. A decrease in DWV titers may therefore be more pronounced over a longer period than was covered by this study.
BQCV fails to replicate after injection into adult bees, but it is able to multiply after injection into pupae (40), and SBV has been detected in mites and their saliva (43), so a potential for Varroamediated transmission of BQCV and SBV exists. In this study, however, mite infestation had no effect or explanatory influence on BQCV titers in pupae or adult bees, and none of the correlations involving BQCV titers in infested pupae, their mites, or the mite infestation rates were significant. The mere detection of BQCV in the Varroa samples is inconclusive with regard to the ability of Varroa to transmit BQCV, since RT-qPCR detection does not distinguish between infectious particles and fragmented noninfectious BQCV RNA. The same conclusion, with the same arguments, can be reached for SBV, which, moreover, could not even be detected in any of the mites examined in this study. However, the probability of successful transmission depends on the amount of virus transmitted, and the BQCV and SBV titers were very low in all sample types. Although it is safe to conclude that in these experiments Varroa-mediated BQCV or SBV transmission did not occur, we cannot rule out the possibility that such transmission could occur with higher titers of these viruses. Our results are consistent with other studies that have been unable to demonstrate an epidemiological relationship between the Varroa mites and BQCV or SBV (19,40,43,49,50). DWV, BQCV, and SBV each react differently to Varroa mite infestation, despite similarities in particle and genome structures (12). Therefore, knowledge of virus epidemiology and, more specifically, virus-vector interactions is important in order to implement effective techniques to manage different virus infections. This complex system of stressors, including Varroa mites, the viruses associated with Varroa, and the variety of harsh treatments used to control Varroa, needs further investigation, especially for implementing effective honey bee health control strategies. One obvious line of research resulting from these experiments includes measuring the effects of tau-fluvalinate on virus titers in bees directly, using Varroa-free colonies, to better differentiate between the effects of Varroa transmission and acaricide treatment on colony DWV titers. Such experiments are currently being planned.