Gene-knockdown in the honey bee mite Varroa destructor by a non-invasive approach: studies on a glutathione S-transferase
© Campbell et al; licensee BioMed Central Ltd. 2010
Received: 26 July 2010
Accepted: 16 August 2010
Published: 16 August 2010
The parasitic mite Varroa destructor is considered the major pest of the European honey bee (Apis mellifera) and responsible for declines in honey bee populations worldwide. Exploiting the full potential of gene sequences becoming available for V. destructor requires adaptation of modern molecular biology approaches to this non-model organism. Using a mu-class glutathione S-transferase (Vd GST-mu1) as a candidate gene we investigated the feasibility of gene knockdown in V. destructor by double-stranded RNA-interference (dsRNAi).
Intra-haemocoelic injection of dsRNA-Vd GST-mu1 resulted in 97% reduction in Vd GST-mu1 transcript levels 48 h post-injection compared to mites injected with a bolus of irrelevant dsRNA (LacZ). This gene suppression was maintained to, at least, 72 h. Total GST catalytic activity was reduced by 54% in Vd GST-mu1 gene knockdown mites demonstrating the knockdown was effective at the translation step as well as the transcription steps. Although near total gene knockdown was achieved by intra-haemocoelic injection, only half of such treated mites survived this traumatic method of dsRNA administration and less invasive methods were assessed. V. destructor immersed overnight in 0.9% NaCl solution containing dsRNA exhibited excellent reduction in Vd GST-mu1 transcript levels (87% compared to mites immersed in dsRNA-LacZ). Importantly, mites undergoing the immersion approach had greatly improved survival (75-80%) over 72 h, approaching that of mites not undergoing any treatment.
Our findings on V. destructor are the first report of gene knockdown in any mite species and demonstrate that the small size of such organisms is not a major impediment to applying gene knockdown approaches to the study of such parasitic pests. The immersion in dsRNA solution method provides an easy, inexpensive, relatively high throughput method of gene silencing suitable for studies in V. destructor, other small mites and immature stages of ticks.
The European honey bee, Apis mellifera, is vital to the pollination of agricultural and wild plants . There is widespread concern about the worldwide decline in the abundance of A. mellifera. The ectoparasitic Varroa mite (Varroa destructor) is the most important pest of A. mellifera and plays a central role to honey bee losses . V. destructor originally parasitized the Asian bee (A. cerana) where it nearly exclusively parasitized the male bees (drones), thus making little impact on the bee colony. However, in the 1950's V. destructor shifted to the European honey bee (A. mellifera) upon which it parasitizes both the drones and female bees (workers). This shift in parasitized caste is significant because the workers make up the bulk of the adult bee population within a colony . V. destructor entered mainland Europe in the 1970's, the USA in 1987 and the UK in 1992 and subsequently has been associated with the loss of millions of colonies . The mite causes damage by feeding on the haemolymph of both the developing bee within brood cells and the adult bee. Mites also transmit a variety of pathogens, most notably deformed wing virus. In terms of both the number of enterprises affected and the impact of global food production, varroasis is arguably the most serious disease of livestock in any species. Previous control of V. destructor by chemical treatment is increasingly ineffective due to the development of widespread resistance in mites to the limited available acaricides .
There is an urgent need to harness modern molecular techniques for research into the biology and, ultimately, the control of this non-model organism, V. destructor. RNA interference (RNAi) is a gene silencing technique that is becoming an ever more powerful tool in investigating the functional role of specific genes that may be potential targets for chemotherapeutic intervention. The RNAi mechanism involves the in vivo production of small interfering RNA molecules (siRNAs) from larger introduced double-stranded RNA (dsRNA). siRNA molecules target and destroy specific mRNA, silencing the target gene at the post-transcriptional stage. RNAi has been used successfully to investigate genes in a range of invertebrate species to date including gut proteins in tsetse flies ; moulting genes in Tribolium castaneum, aquaporins in aphids , a putative prostaglandin E synthase in sealice  and extensively in ticks (reviewed in ). In most cases, administration of dsRNA to the invertebrate target is achieved by intrahaemocoelic injection, though immersion in solutions containing the dsRNA has been employed for nematodes and sea lice.
In this study, we assessed the feasibility of gene knockdown in V. destructor using dsRNA and determined the speed and persistence of the gene knockdown effect. Additionally, we examined less invasive methods of dsRNA administration to increase ease of experimentation and improve mite survivability. The chosen target gene was a glutathione S-transferase of the mu class. GSTs were selected for study as they have been implicated in acaricide resistance in ticks [11, 12] and scabies mites [13, 14] and were readily amenable to activity assays to assess the effect of gene knockdown at the translational step.
Bioinformatics and phylogeny
Tissue distribution of three V. destructor GSTs
Effectiveness and persistence of Vd GST-mu1 knockdown in V. destructor
Effectiveness and mortality rates of different routes of dsRNA administration
Mortality of V. destructor and the degree of gene knockdown under different RNAi experimental conditions.
Survival at 48 h (%)
Vd GSTMu knockdown (%)a
dsRNA in 0.9% NaCl
79 (n = 14)
0% (n = 3)
dsRNA in 0.2% Triton-X100
0 (n = 6)
0 (n = 15)
0 (n = 15)
80 (n = 15)
dsRNA in 0.9% NaCl
76 (n = 21)
87 ± 0.7% (3)
50 (n = 62)
96 ± 3.7% (4)
Phenotype of Vd GST-mu1 knockdown mites
GSTs are a large family of multifunctional enzymes  that catalyse the conjugation of reduced glutathione to the electrophilic centres of lipophilic compounds rendering them water soluble, less toxic products that can be rapidly excreted from the organism's body. As such, GSTs play an important role in phase II detoxification of both endogenous and xenobiotic compounds such as pesticides. Specifically, GSTs have been shown to be involved in acaricide resistance in ticks [11, 12], scabies mites [13, 14] and plant mites . In the present study we demonstrated the presence and tissue expression of two mu- and one kappa-class GST in V. destructor. The largest sequence fragment termed, Vd GST-mu1 was 100% identical to the sequence present in an ongoing pyrosequencing project (Dr. Jay Evans, personal comm.) performed on V. destructor samples from the USA. Vd GST-mu1 was present in all mite tissues assessed and was, by far, the most abundant of the three GSTs examined. Vd GST-mu1 contained conserved glutathione and substrate binding residues and motifs characteristic of mu-class GSTs. However, both the mu-loop and SNAIL/TRAIL motifs were notably different from those of other acarine muGSTs. Interestingly, Vd GST-mu1 was more similar to tick muGSTs than that of mites.
Vd GST-mu1 was used as a test gene to assess the feasibility of dsRNAi gene knockdown in V. destructor. Essentially total suppression of Vd GST-mu1 expression was achieved within 48 h after an intrahaemocoelic bolus injection of dsRNA and this suppression was maintained until at least 72 h. RNAi by injected dsRNA is an extensively employed approach for research in ticks (reviewed in ) and total gene knockdown can be maintained for up to, at least, 9 days . However, V. destructor are considerably smaller (1 mm length × 1.5 mm width) than adult ticks and injection was only possible by pulled-glass capillary needles connected to a relatively sophisticated micro-injector system, rather than using manually operated micro-syringes as is the case for ticks and most insects. Even though the micro-injector system allowed us to demonstrate effective dsRNAi in V. destructor, survival of mites undergoing this very invasive and traumatic procedure was only 50% over a 48 h period.
Increasing the osmotic concentration of the immersion solution by using 0.9% NaCl allowed the mites to remain under solution overnight at 4°C with 75-80% survival and achieved excellent gene silencing. It is not known what is the site of dsRNA entry into the mite, but this is of considerable interest and importance. We presumed entry to be either by ingestion or though the tracheal system via the spiracles. Adult female V. destructor enter the prepupal cells 24 - 48 h before the cell is capped and they remain hidden from vigilant nurse bees by submerging themselves in the liquid brood food  with the peritremes surrounding the spiracles protruding from the liquid enabling the mites to breathe [19, 20]. It may be possible to improve survivability and, possibly increase gene knockdown efficacy, if the osmolality of the immersion solution were increased to that of liquid brood food that is presumed to be greater than that of physiological saline (0.9% NaCl). Indeed, further optimisation of this approach could be attempted, if thought necessary, by altering immersion temperature and time, but the method presented here is a good starting point.
In conclusion, we have established dsRNA knockdown in the critically important bee parasite, V. destructor, using a mu-class GST as the candidate gene. A simple, inexpensive, high throughput method of dsRNA administration was established with low mortality rates for V. destructor. This immersion approach may be applied to other small acarine pests such as scabies mites, sheep scab mites, poultry mites and plant pest mites. Additionally, the approach can be employed in ticks where the immersion approach is more convenient than injections for adult ticks but can also be applied to the smaller larval and nymphal stages which are presently not included in such gene knockdown studies. The degree of knockdown was essentially complete and persisted for several days and was apparent at both the transcript and protein level. The dsRNAi approach will allow researchers to make full use of the V. destructor gene sequences becoming available to investigate biological questions in this important parasite, notably to develop control strategies in the era when current suitable acaricides are no longer effective due to widespread resistance. Peculiarly to Varroa, because of the relatively close relatedness between the parasite and host (Varroa and bee) it is especially difficult to develop pesticides effective against the Varroa, but innocuous to the bee. In this respect, the specificity of dsRNAi may allow targeting of V. destructor genes as a control strategy, as suggested for other insect pests .
Mite collection and husbandry
Varroa destructor mites were collected from capped brood cells of Apis mellifera from hives in York, England that had purposefully been untreated for Varroa control. Prior to harvesting mites, the bee frames were maintained in Aberdeen at 27ºC and 40% relative humidity with a 15.5 h : 8.5 h, light:dark regime. Mites were either removed from newly emerged adult bees or from pupae following uncapping of the brood cell and extraction of the pupae. Mites were collected, bisected and placed in RNA ZR extraction buffer (Zymo Research, Orange, California, USA) and stored at -80 ºC prior to initial total RNA isolation.
RT-PCR of GSTs from V. destructor
Six adult (pharate and mobile) and 2 deutonymph mites were homogenized in 600 μl extraction buffer and total RNA prepared using a mini RNA Isolation II Kit (Zymo Research), as per manufacturer's instructions. Eluted RNA was co-precipitated with glycogen in 95% ethanol and resuspended in 10 μl of DEPC-treated water. After isolation, 2 μg total RNA was DNase treated with 2 μl (2U) RQ1-DNase (Promega, Southampton, UK) and 2 μl RQ1 buffer and incubated at 37°C for 30 min. DNase-treated total RNA (2 μg) was incubated at 70°C with 0.5 μg of oligo d(T)15 (Promega) in a total volume of 10 μl for 5 min. Material was snap-chilled on ice for 5 min prior to the addition of 5 μl 5×RT buffer, 1 μl dNTPs (25 mM each), 0.5 μl Bioscript-reverse transcriptase and DEPC water to 25 μl. The reaction was incubated at 42°C for 60 min prior to arrest by heating to 70°C for 5 min.
A V. destructor database from an ongoing pyrosequencing project (Dr. Jay Evans, pers. comm.) was mined for putative V. destructor GST sequences by BLASTx analysis using GST sequences from the deer tick, Ixodes scapularis. Of 14 putative V. destructor GST partial sequences identified, we chose three to evaluate further: two mu-class GSTs (Vd GST-mu1 and Vd GST-mu2) and one kappa class (Vd GST-kap). These Varroa transcripts were isolated by RT-PCR using primers based on these sequences. PCR reactions consisted of 1 μl cDNA template, 5 μl 10 × reaction buffer, 2 μl 50 mM MgCl2, 1 μl dNTPs (25 mM each), 1 μl each primer (10 mM each), 0.5 μl (1.25 U) Taq (Bioline, London,UK) and DEPC-treated water to give a 50 μl total volume. Primer pairs used were Vd GST-mu-F1 (CAAGTTAAATCTTGCGTTTCC) and Vd GST-mu1-R1 (CTGTCCGTCTTCATGTATGC), Vd GST-mu2-F1 (ACGATTTATCCGTTTTGACG) and Vd GST-mu2-R2 (CCAACTGATGTGTCGTTCC), Vd GST-Kap-F1 (TACTGGTGGTCGTTTCAGG) and Vd GST-Kap-R1 (CTTCATAGGCCAAGAGATGC), generating products of 494, 396 and 343 bp, respectively. PCR cycling conditions were as follows: 1 cycle of 5 min at 94°C, followed by 35 cycles of 1 min at 94°C, 1 min at 53°C and 45 s at 72°C followed by a final extension time of 15 min at 72°C. Products were visualised on an agarose gel and specific bands excised and cloned into pCR4-TOPO TA vector (Invitrogen). Purified plasmids were submitted for sequencing to Eurofins MWG (Ebersberg, Germany) from flanking T7 and T3 promoter regions to confirm the original sequence data.
The GST class of the different Varroa sequences was determined by database searches and phylogenetic analysis. For Vd GST-mu1, protein sequences from various acari were extracted from both non-redundant and EST databases at GenBank by tBLASTx. Additionally, the most similar GST sequences in four model insect species for which complete genomes are available were also obtained. Protein sequences were aligned in CLUSTALW2 . These sequences were used to estimate phylogeny with neighbour-joining, minimum evolution and maximum parsimony methods using MEGA version 4.1 . Phylogenetic trees were constructed with 10,000 bootstrap replicates.
Expression of GST transcripts in tissues of V. destructor
V. destructor mites, were removed from brood cells, as described above, and dissected under ice-cold dissection buffer (20 mM TRIS, 5 mM EDTA, 0.9% NaCl, pH 7.4). Malpighian tubules, gut and synganglia and were dissected from 20, 25 and 98 mites respectively, washed in fresh ice-cold dissection buffer and stored in 50 μl RNAlater (Sigma, Poole, UK) at -80 °C. Prior to RNA extraction, an additional 450 μl dissection buffer was added to sample tubes and centrifuged at 10,000 g for 15 min at 4°C. Supernatant was removed and tissue washed with fresh dissection buffer before a final spin at 10,000 g for 15 min at 4°C. Supernatant was discarded and 600 μl RNA ZR extraction buffer added to each tissue sample. Total RNA extraction, DNase treatment and reverse transcription were as described above. Vd GST-mu1, Vd GST-mu2 and Vd GST-kap transcripts in tissue-specific cDNA were assayed by PCR. PCR protocol and cycling conditions were carried out using the Vd GST-mu1, Vd GST-mu2 and Vd GST-kap primers, as described above over 30 cycles. Gel loading was normalized to V. destructor actin (Vd ActinF, CATCACCATTGGTAACGAG and Vd ActinR, CGATCCAGACGGAATACTT) generating a fragment of approximately 195 bp.
Preparation of dsRNA
Vd GST-mu1-dsRNA was prepared using a BLOCK-iT RNAi TOPO transcription kit (Invitrogen), according to the manufacturer's instructions. LacZ-dsRNA was prepared and used as a negative control. Briefly, PCR was carried out as described above using adult female V. destructor cDNA in conjunction with Vd GST-mu1 specific primers (Vd GST-mu1 F1/R1), or with control LacZ-plasmid and LacZ specific primers (LacZ-F2, ACCAGAAGCGGTGCCGGAAA and LacZ-R2, CCACAGCGGTGGTTCGGAT). Products were resolved on an agarose gel, excised and purified using a Qiagen gel extraction kit (Qiagen, Crawley, UK). TOPO-T7 linker was ligated to Vd GST-mu1 and LacZ reactions before a secondary PCR was carried out to gain sense and antisense templates. T7-RNA polymerase was used in transcription reactions with target templates to generate sense and antisense RNA. Finally, RNA strands were annealed and the resultant dsRNA purified and quantified in a micro-spectrophotometer (Nanodrop Technology Ltd). dsRNA was ethanol precipitated and resuspended in DEPC-treated water to a working concentration of 2.5 μg/μl and stored at -80°C.
Protocol of dsRNA injection and soaking
Adult female V. destructor were removed from capped brood cells along with associated bee larvae. Microinjections were carried out using pulled glass capillary needles in conjunction with a Harvard micro-injector system. Mites were placed on double-sided tape ventral side up, and injected with 20 nl (2.5 μg/μl) of either Vd GST-mu1-dsRNA or LacZ-dsRNA in either the soft tissue proximal to the anal region and postcoxal plate, or in the coxa IV region, as indicated in Figure 7. Needles were left in each mite for 1 - 2 min to reduce the expulsion of fluid from the wound and withdrawn slowly. Mites were left for 1 - 2 min to allow the injection site to "seal" then returned to Petri dishes containing 1 bee larvae per 4 mites. Dead or unhealthy looking mites were removed after 1 hour and mortality was monitored over 72 h in LacZ-dsRNA, Vd GSTmu1-dsRNA and non-injected mites.
To assess non-invasive techniques for dsRNA delivery, mites were either completely immersed in dsRNA or were exposed to a droplet of dsRNA on their ventral carapace. For soaking experiments, adult mites were removed from capped brood cells and placed in 500 μl microfuge tubes containing 20 μl Vd GST-mu1-dsRNA or LacZ-dsRNA (2.5 μg/μl) supplemented with either nothing, 0.9% NaCl, 0.2% Triton-X100 or both. Mites were soaked at 4ºC overnight before being removed, dried and placed in Petri dishes at 27ºC, 95% relative humidity with bee larvae. Alternatively, a sample of mites was exposed to dsRNA by attaching them to double-sided tape and placing a 1 μl drop of Vd GST-mu1-dsRNA or LacZ-dsRNA (2.5 μl/μg) supplemented with either nothing, 0.9% NaCl, 0.2% Triton-X100 or both on the ventral carapace. Mortality was monitored for 48 h prior to collection and validation of knockdown.
Validation of RNAi
To validate RNAi in injected, soaked and droplet-exposed adults, the total RNA was extracted from individual mites 48 h post-treatment. In addition, persistence of the RNAi effect was measured in injected mites by harvesting at 18, 24, 48 and 72 h. Total RNA was extracted from mites using mini-RNA isolation kit II, prior to DNase-treatment and reverse transcription, as described previously. PCR was carried out using either Vd GST-mu1-dsRNA or LacZ-dsRNA treated sample cDNA in conjunction with primers specific for actin or Vd GST-mu1 using primers and cycling conditions, as described above. Products were visualized on an agarose gel normalized to actin loading. To assess persistence and approximate % knockdown, ImageJ software was used to carry out semi-quantitative densitometric analysis on gel images.
Glutathione S-transferase activity in RNAi silenced mites
An enzymatic assay was used to assess the effect of silencing Vd GST-mu1 transcript on GST enzyme activity in dsRNA-treated mites. GST activity was measured by using monochlorobimane (MCB) as the substrate, as previously described for scabies mites  and plant mites , but here we adapted it for cuvettes rather than microplates. MCB forms a stable fluorescent conjugate with reduced glutathione when catalyzed with GST. Mites injected with LacZ-dsRNA or Vd GST-mu1-dsRNA were harvested 48 h post-treatment and frozen at -80ºC. Mites were grouped into pairs and homogenized by milling 2 × 30 s, with a single 3 mm tungsten bead in 175 μl of ice cold 0.05 M Tris buffer pH7.5 in a mechanical tissue disruptor (TissueLyser, Qiagen). After milling, samples were disrupted further in a sonicating water bath for 3 × 20 s with chilling in between bursts then centrifuged at 10,000 g for 5 min at 4ºC. The supernatant was retained, subjected to three freeze thaw-cycles and sonicated for a further 2 × 20 s. Samples were then centrifuged at 10,000 g for 5 min at 4ºC and 150 μl of supernatant from each sample used in the GST assay. Reduced glutathione was dissolved in 0.05 M Tris pH 7.5 and MCB dissolved in methanol and added to each reaction to a final concentration of 3 mM and 0.3 mM, respectively, in a total reaction volume of 200 μl. Mite homogenate and substrate were incubated at 37ºC for 10 min with fluorescence of GSH-bimane adduct measured every 2 s using a F4500 fluorimeter (Hitachi) with emission wavelength at 465 nm and excitation wavelength at 390 nm. All results presented were corrected for non-enzymatic fluorescence obtained from samples substituting 150 μl 0.05 M Tris for the mite material.
We are indebted to Dr. Jay Evans and Dr. Scott Cornman, USDA-ARS Bee Research Laboratories, Beltsville, MD, USA for access to GST sequences from their Varroa sequencing project prior to its collation, analysis and annotation and publication . This project was sponsored by the Scottish Beekeepers Association. We gratefully acknowledge funding from a Genesis Faraday SPARK Award, funded as part of a Scottish Government SEEKIT project; the Department for Environment, Food and Rural Affairs and the Welsh Assembly Government.
- Klein A-, Vaissière BE, Cane JH, Steffan-Dewenter I, Cunningham SA, Kremen C, Tscharntke T: Importance of pollinators in changing landscapes for world crops. Proceedings of the Royal Society B: Biological Sciences. 2007, 274: 303-313. 10.1098/rspb.2006.3721.PubMed CentralView ArticlePubMedGoogle Scholar
- van Engelsdorp D, Meixner MD: A historical review of managed honey bee populations in Europe and the United States and the factors that may affect them. J Invertebr Pathol. 2010, 103 (SUPPL. 1): S80-S95. 10.1016/j.jip.2009.06.011.View ArticleGoogle Scholar
- Martin SJ: The role of varroa and viral pathogens in the collapse of honey bee colonies: A modelling approach. J Appl Ecol. 2001, 38: 1082-1093. 10.1046/j.1365-2664.2001.00662.x.View ArticleGoogle Scholar
- Winston ML: The honey bee colony: Life history. The hive and the honey bee. Edited by: Graham JM. 1992, Hamilton Illinois: Dadant & Sons, 10Google Scholar
- Martin SJ: Acaricide (pyrethroid) resistance in Varroa destructor. Bee World. 2004, 85: 67-69.View ArticleGoogle Scholar
- Walshe DP, Lehane SM, Lehane MJ, Haines LR: Prolonged gene knockdown in the tsetse fly Glossina by feeding double stranded RNA. Insect Mol Biol. 2009, 18: 11-19. 10.1111/j.1365-2583.2008.00839.x.View ArticlePubMedGoogle Scholar
- Minakuchi C, Namiki T, Shinoda T: Krüppel homolog 1, an early juvenile hormone-response gene downstream of Methoprene-tolerant, mediates its anti-metamorphic action in the red flour beetle Tribolium castaneum. Dev Biol. 2009, 325: 341-350. 10.1016/j.ydbio.2008.10.016.View ArticlePubMedGoogle Scholar
- Shakesby AJ, Wallace IS, Isaacs HV, Pritchard J, Roberts DM, Douglas AE: A water-specific aquaporin involved in aphid osmoregulation. Insect Biochem Mol Biol. 2009, 39: 1-10. 10.1016/j.ibmb.2008.08.008.View ArticlePubMedGoogle Scholar
- Campbell EM, Pert CC, Bowman AS: RNA-interference methods for gene-knockdown in the sea louse, Lepeophtheirus salmonis: Studies on a putative prostaglandin E synthase. Parasitology. 2009, 136: 867-874. 10.1017/S0031182009990357.View ArticlePubMedGoogle Scholar
- de la Fuente J, Kocan KM, Almazán C, Blouin EF: RNA interference for the study and genetic manipulation of ticks. Trends Parasitol. 2007, 23: 427-433. 10.1016/j.pt.2007.07.002.View ArticlePubMedGoogle Scholar
- Enayati AA, Asgarian F, Amouei A, Sharif M, Mortazavi H, Boujhmehrani H, Hemingway J: Pyrethroid insecticide resistance in Rhipicephalus bursa (Acari, Ixodidae). Pestic Biochem Physiol. 2010, 97: 243-248. 10.1016/j.pestbp.2010.03.003.View ArticleGoogle Scholar
- Enayati AA, Asgarian F, Sharif M, Boujhmehrani H, Amouei A, Vahedi N, Boudaghi B, Piazak N, Hemingway J: Propetamphos resistance in Rhipicephalus bursa (Acari, Ixodidae). Vet Parasitol. 2009, 162: 135-141. 10.1016/j.vetpar.2009.02.005.View ArticlePubMedGoogle Scholar
- Molin EU, Mattsson JG: Effect of acaricides on the activity of glutathione transferases from the parasitic mite Sarcoptes scabiei. Parasitology. 2008, 135: 115-123. 10.1017/S0031182007003472.View ArticlePubMedGoogle Scholar
- Mounsey KE, Pasay CJ, Arlian LG, Morgan MS, Holt DC, Currie BJ, Walton SF, McCarthy JS: Increased transcription of glutathione S-transferases in acaricide exposed scabies mites. Parasites and Vectors. 2010, 3: 43-10.1186/1756-3305-3-43.PubMed CentralView ArticlePubMedGoogle Scholar
- Sheehan D, Meade G, Foley VM, Dowd CA: Structure, function and evolution of glutathione transferases: Implications for classification of non-mammalian members of an ancient enzyme superfamily. Biochem J. 2001, 360: 1-16. 10.1042/0264-6021:3600001.PubMed CentralView ArticlePubMedGoogle Scholar
- Stumpf N, Nauen R: Biochemical markers linked to abamectin resistance in Tetranychus urticae (Acari: Tetranychidae). Pestic Biochem Physiol. 2002, 72: 111-121. 10.1006/pest.2001.2583.View ArticleGoogle Scholar
- Campbell EM, Burdin M, Hoppler S, Bowman AS: Role of an aquaporin in the sheep tick Ixodes ricinus: Assessment as a potential control target. Int J Parasitol. 2010, 40: 15-23. 10.1016/j.ijpara.2009.06.010.View ArticlePubMedGoogle Scholar
- Sammataro D, Gerson U, Needham G: Parasitic mites of honey bees: Life history, implications, and impact. Annu Rev Entomol. 2000, 45: 519-548.Google Scholar
- Bruce WA, Delfinado-Baker M, Vincenta DL: Comparative morphology of the peritremes of Varroa and Euvarroa (Varroidae), parasites of honey bees (Apidae). Internat J Acarol. 1997, 23: 13-20. 10.1080/01647959708684114.View ArticleGoogle Scholar
- Donzé G, Guerin PM: Time-activity budgets and space structuring by the different life stages of Varroa jacobsoni in capped brood of the honey bee, Apis mellifera. J Insect Behav. 1997, 10: 371-393. 10.1007/BF02765605.View ArticleGoogle Scholar
- Kulkarni MM, Booker M, Silver SJ, Friedman A, Hong P, Perrimon N, Mathey-Prevot B: Evidence of off-target effects associated with long dsRNAs in Drosophila melanogaster cell-based assays. Nature Methods. 2006, 3: 833-838.PubMedGoogle Scholar
- Whyard S, Singh AD, Wong S: Ingested double-stranded RNAs can act as species-specific insecticides. Insect Biochem Mol Biol. 2009, 39: 824-832. 10.1016/j.ibmb.2009.09.007.View ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, Mcgettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Pasay C, Arlian L, Morgan M, Gunning R, Rossiter L, Holt D, Walton S, Beckham S, McCarthy J: The effect of insecticide synergists on the response of scabies mites to pyrethroid acaricides. PLoS Neglected Tropical Diseases. 2009, 3 (1): e354-10.1371/journal.pntd.0000354. Epub 2009 Jan 6PubMed CentralView ArticlePubMedGoogle Scholar
- Nauen R, Stumpf N: Fluorometric microplate assay to measure glutathione S-transferase activity in insects and mites using monochlorobimane. Anal Biochem. 2002, 303: 194-198. 10.1006/abio.2002.5578.View ArticlePubMedGoogle Scholar
- Cornman RC, Schatz MC, Johnston S, Chen YP, Pettis J, Hunt G, Bourgeois L, Elsik C, Grozinger CM, Evans JD: Genomic survey of the ectoparasitic mite Varroa destructor, a major pest of the honey bee Apis mellifera. BMC Genomics.
- Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ: Jalview Version 2-A multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25: 1189-1191. 10.1093/bioinformatics/btp033.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.