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

Unique brood ester profile in a Varroa destructor resistant population of European honey bee (Apis mellifera)

Varroa destructor is one of the greatest threats to Apis mellifera worldwide and if left untreated will kill a colony in less than three years. A Varroa-resistant population from Gotland, Sweden, has managed to survive for 25 years with little to no Varroa treatment by reducing the mite's reproductive success. The underlying mechanisms of this trait is currently not known, though previous research indicates that it is the honey bee brood, and not adult bee influence, that contributes to this phenotype. As the mite's own reproduction is synchronized with the brood's development though the interception of brood pheromones, it is possible that a change in pheromone profile would disrupt the mite's reproductive timing. To investigate this, we characterized the brood ester pheromone (BEP) profile of our resistant Gotland population compared to a non-resistant control. This was done by extracting and analyzing key cuticular compounds of the BEP using gas chromatography. A significant difference was found immediately after brood capping, indicating a divergence in their pheromonal production at this time point. This is an important step to understanding the mechanisms of the Gotland population's Varroa-resistance and contributes to our global understanding of Varroa destructor infestation and survival.

Ecological and social factors influence interspecific pathogens occurrence among bees

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

Life-history stage determines the diet of ectoparasitic mites on their honey bee hosts

Ectoparasitic mites of the genera Varroa and Tropilaelaps have evolved to exclusively exploit honey bees as food sources during alternating dispersal and reproductive life history stages. Here we show that the primary food source utilized by Varroa destructor depends on the host life history stage. While feeding on adult bees, dispersing V. destructor feed on the abdominal membranes to access to the fat body as reported previously. However, when V. destructor feed on honey bee pupae during their reproductive stage, they primarily consume hemolymph, indicated by wound analysis, preferential transfer of biostains, and a proteomic comparison between parasite and host tissues. Biostaining and proteomic results were paralleled by corresponding findings in Tropilaelaps mercedesae, a mite that only feeds on brood and has a strongly reduced dispersal stage. Metabolomic profiling of V. destructor corroborates differences between the diet of the dispersing adults and reproductive foundresses. The proteome and metabolome differences between reproductive and dispersing V. destructor suggest that the hemolymph diet coincides with amino acid metabolism and protein synthesis in the foundresses while the metabolism of non-reproductive adults is tuned to lipid metabolism. Thus, we demonstrate within-host dietary specialization of ectoparasitic mites that coincides with life history of hosts and parasites.

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

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

Reduced male fertility of an Antarctic mite following extreme heat stress could prompt localized population declines

Climate change is leading to substantial global thermal changes, which are particularly pronounced in polar regions. Therefore, it is important to examine the impact of heat stress on the reproduction of polar terrestrial arthropods, specifically, how brief extreme events may alter survival. We observed that sublethal heat stress reduces male fecundity in an Antarctic mite, yielding females that produced fewer viable eggs. Females and males collected from microhabitats with high temperatures showed a similar reduction in fertility. This impact is temporary, as indicated by recovery of male fecundity following return to cooler, stable conditions. The diminished fecundity is likely due to a drastic reduction in the expression of male-associated factors that occur in tandem with a substantial increase in the expression of heat shock proteins. Cross-mating between mites from different sites confirmed that heat-exposed populations have impaired male fertility. However, the negative impacts are transient as the effect on fertility declines with recovery time under less stressful conditions. Modeling indicated that heat stress is likely to reduce population growth and that short bouts of non-lethal heat stress could have substantial reproductive effects on local populations of Antarctic arthropods.

Host brood traits, independent of adult behaviours, reduce Varroa destructor mite reproduction in resistant honeybee populations

The ectoparasitic mite Varroa destructor is an invasive species of Western honey bees (Apis mellifera) and the largest pathogenic threat to their health world-wide. Its successful invasion and expansion is related to its ability to exploit the worker brood for reproduction, which results in an exponential population growth rate in the new host. With invasion of the mite, wild honeybee populations have been nearly eradicated from Europe and North America, and the survival of managed honeybee populations relies on mite population control treatments. However, there are a few documented honeybee populations surviving extended periods without control treatments due to adapted host traits that directly impact Varroa mite fitness. The aim of this study was to investigate if Varroa mite reproductive success was affected by traits of adult bee behaviours or by traits of the worker brood, in three mite-resistant honey bee populations from Sweden, France and Norway. The mite's reproductive success was measured and compared in broods that were either exposed to, or excluded from, adult bee access. Mite-resistant bee populations were also compared with a local mite-susceptible population, as a control group. Our results show that mite reproductive success rates and mite fecundity in the three mite-resistant populations were significantly different from the control population, with the French and Swedish populations having significantly lower reproductive rates than the Norwegian population. When comparing mite reproduction in exposed or excluded brood treatments, no differences were observed, regardless of population. This result clearly demonstrates that Varroa mite reproductive success can be suppressed by traits of the brood, independent of adult worker bees.

The adverse impact on lifespan, immunity, and forage behavior of worker bees (Apis mellifera Linnaeus 1758) after exposure to flumethrin

Several pyrethroids (such as flumethrin and fluvalinate) with low toxicity to honey bees and comparable high toxicity to mites are used worldwide as acaricides. However, flumethrin has been used for a long time in colonies to control Varroa destructor and the honey bees might be exposed to flumethrin cumulatively, which could affect the health of honey bee colonies. This study evaluated the potential adverse effects of direct flumethrin exposure on worker bees under laboratory and colony conditions. Under laboratory conditions, downregulation of genes related to immune was observed when worker bees were exposed to flumethrin above 1/16 LD₅₀; at levels above 1/8 LD₅₀, olfactory learning was impaired, and genes related to learning memory were downregulated; and at >1/4 LD₅₀, their lifespan was shortened. Monitoring with radio frequency identification (RFID) revealed that worker bees in a colony exposed to flumethrin above 1/8 LD₅₀ had a shortened lifespan and reduced foraging ability. When worker bees are exposed to >1/4 LD₅₀ of flumethrin, it can lead to excessive rest day behavior. These results indicate that applying flumethrin in colonies may pose a severe health risk to honey bees and reveal the urgent need to develop non-toxic and highly effective acaricides.

Developing Strategies to Help Bee Colony Resilience in Changing Environments

Climate change, loss of plant biodiversity, burdens caused by new pathogens, predators, and toxins due to human disturbance and activity are significant causes of the loss of bee colonies and wild bees. The aim of this review is to highlight some possible strategies that could help develop bee resilience in facing their changing environments. Scientists underline the importance of the links between nutrition, microbiota, and immune and neuroendocrine stress resistance of bees. Nutrition with special care for plant-derived molecules may play a major role in bee colony health. Studies have highlighted the importance of pollen, essential oils, plant resins, and leaves or fungi as sources of fundamental nutrients for the development and longevity of a honeybee colony. The microbiota is also considered as a key factor in bee physiology and a cornerstone between nutrition, metabolism, growth, health, and pathogen resistance. Another stressor is the varroa mite parasite. This parasite is a major concern for beekeepers and needs specific strategies to reduce its severe impact on honeybees. Here we discuss how helping bees to thrive, especially through changing environments, is of great concern for beekeepers and scientists.

Ornamental tobacco floral nectar is a rich source of antimicrobial peptides

Although floral nectar is a rich source of nutrients, it is rarely infected by microorganisms. Defense molecules such as proteins have been identified in this fluid, but defense peptides have been largely overlooked. Thus, the aim of this study was to perform an extensive peptidomic analysis of the ornamental tobacco floral nectar to seek peptides involved in nectar defense. Using LC-MS/MS, 793 peptides were sequenced and characterized. After extensive bioinformatics analysis, six peptides were selected for further characterization, synthesis, and evaluation of their antimicrobial properties against phytopathogenic fungi and bacteria. All six peptides had antimicrobial activity to some extent. However, the activity varied by peptide concentration and microorganism tested. An analysis of the action mechanism revealed damage in the cell membrane induced by peptides. The results show that floral nectar is rich in peptides and that, together with proteins and hydrogen peroxide, they contribute to plant defense against microorganisms during pollination.

Effects of Niemann-Pick type C2 (NPC2) gene transcripts silencing on behavior of Varroa destructor and molecular changes in the putative olfactory gene networks

To understand the role of two Niemann-Pick type C2 (NPC2) transcripts, Vd40090 (NP1) and Vd74517 (NP5), in the chemosensing pathway of Varroa destructor, we evaluated the impact of NP5 silencing on mites behavior and compared the effect of silencing of either transcripts on the interaction between chemosensory transcripts. In contrast to silencing NP1, which reduced feeding and reproduction by the mite (Nganso et al., 2021), silencing of NP5 reduced significantly the host reaching ability, but it did not affect the feeding on nurse bee. However, silencing of either transcript changed dramatically the co-expression patterns among the putative chemosensory genes, binding proteins and receptors. The results suggest the role of gustatory receptors in the detection of long-range chemical cues in the chemosensory cascade of the Varroa mite.

Integrative Analysis of lncRNA-mRNA Co-expression Provides Novel Insights Into the Regulation of Developmental Transitions in Female Varroa destructor

Varroa destructor is a major pathogenic driver of the Western honeybee colony losses globally. Understanding the developmental regulation of V. destructor is critical to develop effective control measures. Development is a complex biological process regulated by numerous genes and long non-coding RNAs (lncRNAs); however, the underlying regulation of lncRNAs in the development of V. destructor remains unknown. In this study, we analyzed the RNA sequencing (RNA-Seq) data derived from the four stages of female V. destructor in the reproductive phase (i.e., egg, protonymph, deutonymph, and adult). The identified differentially expressed mRNAs and lncRNAs exhibited a stage-specific pattern during developmental transitions. Further functional enrichment established that fat digestion and absorption, ATP-binding cassette (ABC) transporters, mitogen-activated protein kinase (MAPK) signaling pathway, and ubiquitin-proteasome pathway play key roles in the maturation of female V. destructor. Moreover, the lncRNAs and mRNAs of some pivotal genes were significantly upregulated at the deutonymph stage, such as cuticle protein 65/6.4/63/38 and mucin 5AC, suggesting that deutonymph is the key stage of metamorphosis development and pathogen resistance acquisition for female V. destructor. Our study provides novel insights into a foundational understanding of V. destructor biology.

Towards disrupting Varroa -honey bee chemosensing: A focus on a Niemann-Pick type C2 transcript

We focused our study on the 12 recently identified putative odorant carrier proteins in the ectoparasitic mite, Varroa destructor. Here we show, via an exclusion of the chemosensory appendages (forelegs and gnathosoma) that transcripts of five of the 12 genes were significantly lower, suggesting that they are likely involved in carrying host volatiles. Specifically, three transcripts were found to be foreleg-specific while the other two transcripts were expressed in both the forelegs and gnathosoma. We focused on one of the highly expressed and foreleg-specific transcript Vd40090, which encodes a Niemann-Pick disease protein type C2 (NPC2) protein. Effects of dsRNA-mediated silencing of Vd40090 were first measured by quantifying the transcript levels of genes that encode other putative odorant carrier proteins as well as reproduction related proteins. In addition, the impact of silencing on mites behaviour and survival was tested. Silencing of Vd40090 effectively disrupted Varroa host selection, acceptance and feeding and significantly impaired the expression of genes that regulate its reproduction in brood cells, resulting in reduced reproduction and survival.

A toolkit for studying Varroa genomics and transcriptomics: preservation, extraction, and sequencing library preparation

BACKGROUND

The honey bee parasite, Varroa destructor, is a leading cause of honey bee population declines. In addition to being an obligate ectoparasitic mite, Varroa carries several viruses that infect honey bees and act as the proximal cause of colony collapses. Nevertheless, until recently, studies of Varroa have been limited by the paucity of genomic tools. Lab- and field-based methods exploiting such methods are still nascent. This study developed a set of methods for preserving Varroa DNA and RNA from the field to the lab and processing them into sequencing libraries. We performed preservation experiments in which Varroa mites were immersed in TRIzol, RNAlater, and absolute ethanol for preservation periods up to 21 days post-treatment to assess DNA and RNA integrity.

RESULTS

For both DNA and RNA, mites preserved in TRIzol and RNAlater at room temperature degraded within 10 days post-treatment. Mites preserved in ethanol at room temperature and 4 °C remained intact through 21 days. Varroa mite DNA and RNA libraries were created and sequenced for ethanol preserved samples, 15 and 21 days post-treatment. All DNA sequences mapped to the V. destructor genome at above 95% on average, while RNA sequences mapped to V. destructor, but also sometimes to high levels of the deformed-wing virus and to various organisms.

CONCLUSIONS

Ethanolic preservation of field-collected mites is inexpensive and simple, and allows them to be shipped and processed successfully in the lab for a wide variety of sequencing applications. It appears to preserve RNA from both Varroa and at least some of the viruses it vectors.

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

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

Genome-Wide Identification of Long Non-coding RNAs in the Gravid Ectoparasite Varroa destructor

Long non-coding RNAs (lncRNAs) emerge as critical regulators across a wide variety of biological functions in living organisms. However, to date, no systematic characterization of lncRNAs has been investigated in the ectoparasitic mite Varroa destructor, the most severe biotic threat to honey bees worldwide. Here, we performed an initial genome-wide identification of lncRNAs in V. destructor via high-throughput sequencing technology and reported, for the first time, the transcriptomic landscape of lncRNAs in the devastating parasite. By means of a lncRNA identification pipeline, 6,645 novel lncRNA transcripts, encoded by 3,897 gene loci, were identified, including 2,066 sense lncRNAs, 2,772 lincRNAs, and 1,807 lncNATs. Compared with protein-coding mRNAs, V. destructor lncRNAs are shorter in terms of full length, as well as of the ORF length, contain less exons, and express at lower level. GO term and KEGG pathway enrichment analyses of the lncRNA target genes demonstrated that these predicted lncRNAs may be potentially responsible for the regulatory functions of cellular and biological progresses in the reproductive phase of V. destructor. To our knowledge, this is the first catalog of lncRNA profile in the parasitiformes species, providing a valuable resource for genetic and genomic studies. Understanding the characteristics and features of lncRNAs in V. destructor would promote sustainable parasite control.

Significant transcriptional changes in mature daughter Varroa destructor mites during infestation of different developmental stages of honeybees

BACKGROUND

Varroa destructor is considered a major cause of honeybee (Apis mellifera) colony losses worldwide. Although V. destructor mites exhibit preference behavior for certain honeybee lifecycle stages, the mechanism underlying host finding and preference remains largely unknown.

RESULTS

By using a de novo transcriptome assembly strategy, we sequenced the mature daughter V. destructor mite transcriptome during infestation of different stages of honeybees (brood cells, newly emerged bees and adult bees). A total of 132 779 unigenes were obtained with an average length of 2745 bp and N50 of 5706 bp. About 63.1% of the transcriptome could be annotated based on sequence homology to the predatory mite Metaseiulus occidentalis proteins. Expression analysis revealed that mature daughter mites had distinct transcriptome profiles after infestation of different honeybee stages, and that the majority of the differentially expressed genes (DEGs) of mite infesting adult honeybees were down-regulated compared to that infesting the sealed brood cells. Gene ontology and KEGG pathway enrichment analyses showed that a large number of DEGs were involved in cellular process and metabolic process, suggesting that Varroa mites undergo metabolic adjustment to accommodate the cellular, molecular and/or immune response of the honeybees. Interestingly, in adult honeybees, some mite DEGs involved in neurotransmitter biosynthesis and transport were identified and their levels of expression were validated by quantitative polymerase chain reaction (qPCR).

CONCLUSION

These results provide evidence for transcriptional reprogramming in mature daughter Varroa mites during infestation of honeybees, which may be relevant to understanding the mechanism underpinning adaptation and preference behavior of these mites for honeybees. © 2020 Society of Chemical Industry.

A high quality method for hemolymph collection from honeybee larvae

The drastic decline of bees is associated with several factors, including the immune system suppression due to the increased exposure to pesticides. A widely used method to evaluate these effects on these insects' immune systems is the counting of circulating hemocytes in the hemolymph. However, the extraction of hemolymph from larvae is quite difficult, and the collected material is frequently contaminated with other tissues and gastrointestinal fluids, which complicates counting. Therefore, the present work established a high quality and easily reproducible method of extracting hemolymph from honeybee larvae (Apis mellifera), the extraction with ophthalmic scissors. Extraction methods with the following tools also were tested: 30G needle, fine-tipped forceps, hypodermic syringe, and capillaries tubes. The hemolymph was obtained via an incision on the larvae’s right side for all methods, except for the extraction with ophthalmic scissors, in which the hemolymph was extracted from the head region. To assess the purity of the collected material, turbidity analyses of the samples using a turbidimeter were proposed, tested, and evaluated. The results showed that the use of ophthalmic scissors provided the clearest samples and was free from contamination. A reference range between 22,432.35 and 24,504.87 NTU (nephelometric turbidity units) was established, in which the collected samples may be considered of high quality and free from contamination.

Varroa destructor: A Complex Parasite, Crippling Honey Bees Worldwide

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

Divergent evolutionary trajectories following speciation in two ectoparasitic honey bee mites

Multispecies host-parasite evolution is common, but how parasites evolve after speciating remains poorly understood. Shared evolutionary history and physiology may propel species along similar evolutionary trajectories whereas pursuing different strategies can reduce competition. We test these scenarios in the economically important association between honey bees and ectoparasitic mites by sequencing the genomes of the sister mite species Varroa destructor and Varroa jacobsoni. These genomes were closely related, with 99.7% sequence identity. Among the 9,628 orthologous genes, 4.8% showed signs of positive selection in at least one species. Divergent selective trajectories were discovered in conserved chemosensory gene families (IGR, SNMP), and Halloween genes (CYP) involved in moulting and reproduction. However, there was little overlap in these gene sets and associated GO terms, indicating different selective regimes operating on each of the parasites. Based on our findings, we suggest that species-specific strategies may be needed to combat evolving parasite communities.

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.

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

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

A gene for resistance to the Varroa mite (Acari) in honey bee (Apis mellifera) pupae

Social insect colonies possess a range of defences which protect them against highly virulent parasites and colony collapse. The host-parasite interaction between honey bees (Apis mellifera) and the mite Varroa destructor is unusual, as honey bee colonies are relatively poorly defended against this parasite. The interaction has existed since the mid-20th Century, when Varroa switched host to parasitize A. mellifera. The combination of a virulent parasite and relatively naïve host means that, without acaricides, honey bee colonies typically die within 3 years of Varroa infestation. A consequence of acaricide use has been a reduced selective pressure for the evolution of Varroa resistance in honey bee colonies. However, in the past 20 years, several natural-selection-based breeding programmes have resulted in the evolution of Varroa-resistant populations. In these populations, the inhibition of Varroa's reproduction is a common trait. Using a high-density genome-wide association analysis in a Varroa-resistant honey bee population, we identify an ecdysone-induced gene significantly linked to resistance. Ecdysone both initiates metamorphosis in insects and reproduction in Varroa. Previously, using a less dense genetic map and a quantitative trait loci analysis, we have identified Ecdysone-related genes at resistance loci in an independently evolved resistant population. Varroa cannot biosynthesize ecdysone but can acquire it from its diet. Using qPCR, we are able to link the expression of ecdysone-linked resistance genes to Varroa's meals and reproduction. If Varroa co-opts pupal compounds to initiate and time its own reproduction, mutations in the host's ecdysone pathway may represent a key selection tool for honey bee resistance and breeding.