Genetic variability and pyrethroid susceptibility of the parasitic honey bee mite Varroa destructor (Acari: Varroidae) in Iran

The status of pyrethroid resistance mutation frequencies in Varroa destructor populations in the most important beekeeping areas of Türkiye

The Varroa destructor (hereafter referred to as Varroa) is a major pest of honeybees that is generally controlled using pyrethroid-based acaricides. However, resistance to these insecticides has become a growing problem, driven by the acquisition of knockdown resistance (kdr) mutations in the mite's voltage-gated sodium channel (vgsc) gene. Resistance mutations in the vgsc gene, such as the L925V mutation, can confer resistance to pyrethroids like flumethrin and tau-fluvalinate. Monitoring genotypic resistance through molecular mutation screening is crucial for tracking and mitigating resistance spread. In this study, the frequency of resistance mutations in the vgsc was examined using a Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) approach in Varroa populations sampled throughout the Mediterranean, Aegean, and Black Sea regions of Türkiye. Considering all the samples analyzed, the results demonstrated a mean resistance allele frequency of 83.29%, indicating a relatively high frequency of resistant alleles. We observed 94.58%, 85.71%, and 69.58% resistant allele frequencies in populations sampled from the Mediterranean, Aegean, and Black Sea regions, respectively, in our study. The results of our investigation demonstrated substantial regional variations in the frequencies of resistant alleles among Varroa populations throughout Türkiye, with notably elevated resistance levels observed in the Mediterranean and Aegean regions. Due to the significant resistance mutation frequency differences between both provinces and regions, long-term monitoring of resistance alleles and the planning of regional control strategies are required for effective control of this pest.

Detection of mutations in the voltage-gated sodium channel gene of Varroa destructor (Acari: Varroidae) populations

The honeybee (Apis mellifera) ectoparasite, Varroa destructor, is one of the most important honeybee pests worldwide. Acaricides, including the pyrethroids (tau-fluvalinate, flumethrin), and organophosphate (coumaphos) have been applied to control this mite within apiaries, still the long-term, constant, and excessive use of these products has led to the development of resistance in many populations. Three different mutations (L925V, L925I, L925M) at position 925 and one mutation (M918L) at 918 position of the V. destructor voltage-gated sodium channel (VGSC) have been associated with the resistance to these compounds. In the present study, we examined the presence of resistance mutations in the VGSC gene, encoding the target of pyrethroids, in the V. destructor population collected from the Kayseri and Sivas provinces of Central Anatolia Region of Türkiye. A total of 200 V. destructor samples were collected from 20 apiaries in two provinces throughout 2023. To investigate the mutations in Varroa samples, the domain II region of the VGSC gene was amplified using PCR and sequenced. The nucleotide sequencing of the IIS4-IIS5 linker region of the VGSC gene revealed one amino acid change at position 925: a leucine to isoleucine substitution (L925I). No mutations at other positions were identified. Homozygous resistant alleles were detected in 20 (40 %) of the sequenced 50 samples in the study areas. However, we detected the homozygous sensitive allele (wild-type allele, L925) in the remaining samples (50/30, 60 %). The result shows that this status may indicate a problem for Varroa control in the future. Thus, alternative acaricide with a mode of action of different pyrethroids should be considered in the control of V. destructor populations in these provinces.

Geographical distribution of pyrethroid resistance mutations in Varroa destructor across Türkiye and a European overview

Varroa destructor Anderson & Trueman (Acari: Varroidae) is of paramount significance in modern beekeeping, with infestations presenting a primary challenge that directly influences colony health, productivity, and overall apicultural sustainability. In order to control this mite, many beekeepers rely on a limited number of approved synthetic acaricides, including the pyrethroids tau-fluvalinate, flumethrin and organophosphate coumaphos. However, the excessive use of these substances has led to the widespread development of resistance in various beekeeping areas globally. In the present study, the occurrence of resistance mutations in the voltage-gated sodium channel (VGSC) and acetylcholinesterase (AChE), the target-site of pyrethroids and coumaphos, respectively, was examined in Varroa populations collected throughout the southeastern and eastern Anatolia regions of Türkiye. All Varroa samples belonged to the Korean haplotype, and a very low genetic distance was observed based on cytochrome c oxidase subunit I (COI) gene sequences. No amino acid substitutions were determined at the key residues of AChE. On the other hand, three amino acid substitutions, (L925V/I/M), previously associated with pyrethroid resistance, were identified in nearly 80% of the Turkish populations. Importantly, L925M, the dominant mutation in the USA, was detected in Turkish Varroa populations for the first time. To gain a more comprehensive perspective, we conducted a systematic analysis of the distribution of pyrethroid resistance mutations across Europe, based on the previously reported data. Varroa populations from Mediterranean countries such as Türkiye, Spain, and Greece exhibited the highest frequency of resistance mutation. Revealing the occurrence and geographical distribution of pyrethroid resistance mutations in V. destructor populations across the country will enhance the development of more efficient strategies for mite management.

Mutations associated with pyrethroid resistance in Varroa mite, a parasite of honey bees, are widespread across the United States

BACKGROUND

Managed honey bees are key pollinators of many crops and play an essential role in the United States food production. For more than ten years, beekeepers in the United States have been reporting high rates of colony losses. One of the drivers of these losses is the parasitic mite Varroa destructor. Maintaining healthy honey bee colonies in the United States is dependent on a successful control of this mite. The pyrethroid tau-fluvalinate (Apistan®) was among the first synthetic varroacides registered in the United States. With over 20 years of use, mites resistant to Apistan® have emerged, and so it is unsurprising that treatment failures have been reported. Resistance to tau-fluvalinate in US mite populations is associated with point mutations at position 925 of the voltage-gated sodium channel.

RESULTS

Here, we have generated a distribution map of pyrethroid resistance alleles in Varroa samples collected from US apiaries in 2016 and 2017, using a high throughput allelic discrimination assay based on TaqMan®. Our results evidence that knockdown resistance (kdr)-type mutations are widely distributed in Varroa populations across the country showing high variability among apiaries. We used these data to predict the phenotype of the mites in the case of treatments with pyrethroids.

CONCLUSION

We highlight the relevance of monitoring the resistance in mite populations to achieve an efficient control of this pest. We also put forward the benefits of implementing this methodology to provide data for designing pest management programs aiming to control Varroa. © 2021 Society of Chemical Industry.

Reduced proinsecticide activation by cytochrome P450 confers coumaphos resistance in the major bee parasite Varroa destructor

Varroa destructor is one of the main problems in modern beekeeping. Highly selective acaricides with low toxicity to bees are used internationally to control this mite. One of the key acaricides is the organophosphorus (OP) proinsecticide coumaphos, that becomes toxic after enzymatic activation inside Varroa We show here that mites from the island Andros (AN-CR) exhibit high levels of coumaphos resistance. Resistance is not mediated by decreased coumaphos uptake, target-site resistance, or increased detoxification. Reduced proinsecticide activation by a cytochrome P450 enzyme was the main resistance mechanism, a powerful and rarely encountered evolutionary solution to insecticide selection pressure. After treatment with sublethal doses of [¹⁴C] coumaphos, susceptible mite extracts had substantial amounts of coroxon, the activated metabolite of coumaphos, while resistant mites had only trace amounts. This indicates a suppression of the P450 (CYP)-mediated activation step in the AN-CR mites. Bioassays with coroxon to bypass the activation step showed that resistance was dramatically reduced. There are 26 CYPs present in the V. destructor genome. Transcriptome analysis revealed overexpression in resistant mites of CYP4DP24 and underexpression of CYP3012A6 and CYP4EP4 RNA interference of CYP4EP4 in the susceptible population, to mimic underexpression seen in the resistant mites, prevented coumaphos activation and decreased coumaphos toxicity.

Genetic analysis and screening of pyrethroid resistance mutations in Varroa destructor populations from Turkey

Varroa destructor is the most common ectoparasite of the Western honey bee (Apis mellifera L.) worldwide and poses a serious threat to bee health. Synthetic acaricides, particularly pyrethroids, are frequently used to control Varroa mites. However, long-term and repeated use of synthetic pyrethroids has led to the development of resistance. In this study, we report on the presence of resistance mutations in the voltage-gated sodium channel in V. destructor populations from Turkish beekeeping areas. Two resistance mutations, L925V and L925I, that were previously associated with pyrethroid resistance, were found in more than 75% of the populations. A general correlation between the presence of mutations and the history of acaricide usage was observed for the sampled hives. In addition, we show there is only a low genetic distance among the sampled V. destructor populations, based on the analysis of three mitochondrial genes: cytochrome b (cytb), ATP synthase subunit 6 (atp6), and cytochrome c oxidase subunit III (cox3). Revealing the presence and geographical distribution of pyrethroid resistance mutations in V. destructor populations from Turkish apiaries will contribute to create more effective mite management programmes.

Genetic variation and heteroplasmy of Varroa destructor inferred from ND4 mtDNA sequences

Varroa destructor, a parasitic mite of the western honey bee, Apis mellifera L., is a serious threat to colonies and beekeeping worldwide. Population genetics studies of the mite have provided information on two mitochondrial haplotypes infecting honey bee colonies, named K and J (after Korea and Japan, respectively, where they were originally identified). On the American continent, the K haplotype is much more prevalent, with the J haplotype only detected in some areas of Brazil. The aims of the present study were to assess the genetic diversity of V. destructor populations in the major beekeeping region of Argentina and to evaluate the presence of heteroplasmy at the nucleotide level. Phoretic mites were collected from managed A. mellifera colonies in ten localities, and four mitochondrial DNA (mtDNA) regions (COXI, ND4, ND4L, and ND5) were analyzed. Based on cytochrome oxidase subunit I (COXI) sequencing, exclusively the K haplotype of V. destructor was detected. Furthermore, two sub-haplotypes (KArg-N1 and KArg-N2) were identified from a variation in ND4 sequences and the frequency of these sub-haplotypes was found to significantly correlate with geographical latitude. The occurrence of site heteroplasmy was also evident for this gene. Therefore, ND4 appears to be a sensitive marker for detecting genetic variability in mite populations. Site heteroplasmy emerges as a phenomenon that could be relatively frequent in V. destructor.

Coexistence of genetically different Varroa destructor in Apis mellifera colonies

The aim of this study was to investigate the genetic diversity of Varroa destructor parasitizing Apis mellifera colonies and to test for possible host-parasite association at the mitochondrial DNA (mtDNA) level. Six A. mellifera haplotypes (including a novel C2aa) and five haplotypes of V. destructor were detected in 29 analyzed colonies from eight sampling sites in Serbia. We revealed the presence of the K and S1 haplotypes as well as KS1 and KP1 heteroplasmic mite individuals in all localities, while the P1 haplotype was only found in four sampling sites. Significant differences in V. destructor genetic diversity were found at both apiary and colony levels, with mite haplotypes coexisting in almost all tested colonies. In addition, a significant correlation between the number of analyzed mites per colony and the number of identified V. destructor haplotypes was observed. However, no significant host-parasite relationship was found, suggesting that mites bearing different haplotypes as well as those heteroplasmic individuals are well adapted to the host, A. mellifera, independently of the identified haplotype present in each colony. Our results will contribute to future population and biogeographic studies concerning V. destructor infesting A. mellifera, as well as to better understanding their host-parasite relationship.