qPCR assays with dual-labeled probes for genotyping honey bee variants associated with varroa resistance
BMC Veterinary Research volume 17, Article number: 179 (2021)
The varroa mite is one of the main causes of honey bee mortality. An important mechanism by which honey bees increase their resistance against this mite is the expression of suppressed mite reproduction. This trait describes the physiological inability of mites to produce viable offspring and was found associated with eight genomic variants in previous research.
This paper presents the development and validation of high-throughput qPCR assays with dual-labeled probes for discriminating these eight single-nucleotide variants. Amplicon sequences used for assay validation revealed additional variants in the primer/probe binding sites in four out of the eight assays. As for two of these the additional variants interfered with the genotyping outcome supplementary primers and/or probes were developed. Inclusion of these primers and probes in the assay mixes allowed for the correct genotyping of all eight variants of interest within our bee population.
These outcomes underline the importance of checking for interfering variants in designing qPCR assays. Ultimately, the availability of this assay allows genotyping for the suppressed mite reproduction trait and paves the way for marker assisted selection in breeding programs.
Since the first occurrence of the Varroa destructor mite in the Western honey bee Apis mellifera, honey bee health has become tightly interwoven with the presence and abundance of this ectoparasitic mite . Originally the varroa mite occurred in South-East Asia where it has a balanced host-parasite relationship with the Asian honey bee Apis cerana . Arriving in Europe around the 80’s  the varroa mite encountered a large pool of susceptible hosts lacking natural resistance [3, 4]. In addition to feeding on the fat body of bees , the mite provides a new transmission pathway for viruses causing rising problems with virus infections and ultimately colony mortality . The initial response of the beekeeping community to control the varroa mite was heavily reliant on chemicals . Although being effective, these chemicals also harm honey bees  and include disadvantages such as the deposition of residues in hive products  and the prevention of co-evolutionary processes to create a stable host-parasite relationship . A long-term solution overcoming these disadvantages is to find and select on honey bee traits linked with varroa resistance or tolerance . Over the last decades several research and breeding programs discovered multiple of these traits [2, 11, 12].
One of the key traits linked with varroa resistance is suppressed mite reproduction (SMR) which describes the non-reproduction of varroa mites in honey bee drone pupae cells [7, 13]. The precise mechanisms behind SMR are still not fully understood. Possible hypotheses are a suppression of the varroa reproduction cycle by lower levels of juvenile hormone , alterations in a gene from the ecdysone pathway  or diminished production of the brood pheromone . Varroa reproduction may also be influenced by variations in the genotype of the mite or in the physiological status of the brood cell [16, 17]. In order to initiate breeding programs on SMR it is important to start by screening honey bee populations for the presence of the trait . As an alternative for performing elaborative phenotypic assays, genotypic information can greatly increase the scale at which local populations can be screened. In addition, genotypic information can provide crucial insight in the mechanisms underlying varroa resistance .
Since the publication of the honey bee genome in 2006  many studies identified quantitative trait loci or single nucleotide variants (SNV) associated with different varroa resistance traits . For SMR, eight single-nucleotide variants were discovered by Broeckx and colleagues  using a novel whole exome sequencing design. Of the variants discovered six were risk associated variants and two were protective variants. The present research describes the design and validation of eight dual-labeled probe based qPCR assays for the high-throughput genotyping of the SMR trait.
Table 1 provides an overview of the primer and probe sequences of the qPCR assays used to genotype the eight SNVs associated with SMR, along with their amplicon lengths and optimal Ta. The resulting amplification plots are shown in Fig. 1. For the assays genotyping SNV 1, SNV 2, SNV 6 and SNV 8 no additional SNVs were present in the primer/probe binding sites in the 92 sequenced worker bees distributed throughout Flanders. For the assay genotyping SNV 3, a rare G > A variant was detected 3 bp downstream SNV 3 in the probe binding site of the Wt-allele in one bee, but it did not influence correct genotyping as the Wt-probe could still specifically bind to the Wt-allele and produce a strong signal, despite its single mismatch with the target. Similarly, for the assay genotyping SNV 7, a rare C > T variant was detected 3 bp upstream SNV 7 in the probe binding site of the Wt-allele in two bees, but did not influence correct genotyping either.
However, for the assay genotyping SNV 5, an abundant C > G variant was detected 3 bp downstream SNV 5 in the probe binding site of the Vt-allele in 57 bees, that prevented binding of the original Vt1-probe to the mismatch-containing Vt-allele and thus the detection of the Vt-allele. Adding the Vt2-probe to the assay, containing that additional variant, solved this problem. Most problematic was the assay genotyping SNV 4. A rare C > T variant (found in one bee) at the fourth last position of the F1-primer prevented amplification of the Vt-allele. In addition, an abundant G > A variant 5 bp downstream SNV 4 in the probe binding site of the Vt-allele (found in 35 bees) prevented the detection of the Vt-allele using the original Vt1-probe. Including the F2-primer and the Vt2-probe to the assay mix (both containing the corresponding additional variant), resulted in correct genotyping.
These results show that additional variants, which are common in honey bees due to their extremely high recombination rate , can lead to false results depending on their nature and position, as they can in all PCR-based techniques. Also, this finding underscores the importance of taking into account known variants from available sequence data or variant databases when designing assays and of checking results for anomalous amplification curves. Haploid samples (males) should always have a signal with only one probe. Diploid samples (females) should always have a typical homozygous signal (steep slope and high RFU-value) with one probe and no signal with the other probe, or a typical heterozygous signal (intermediate slope and intermediate RFU-value) with both probes . Sequencing the complete qPCR amplicon region of samples with anomalous amplification plots with external primers, as described in Broeckx et al. , can identify unknown influencing variants.
Although sequencing is more informative, qPCR assays with dual-labeled probes allow for faster and cheaper high-throughput screening. For screening new populations, it is recommended to first sequence the qPCR amplicon regions in a small representative subset of that population and adjust the qPCR assays based on newly identified variants before using them for high-throughput screening. The assays can be performed on other qPCR platforms and reagents from other suppliers, but we advise to check the specific annealing temperature for every assay experimentally in the lab-specific setup with all positive controls and a no template control (NTC). In case of a lack of control samples, artificial oligos can be ordered containing the correctly orientated primer and probe sequences.
There are multiple options to tune up the assays. Assays can be run in the presence of fluorescent nucleic acid binding dyes such as SYBR green, to include melt curve analysis, useful to detect (non) specific amplification, oligo dimers, null-alleles or primer/probe-related problems during optimization, as described by Van Poucke et al. . Assays performed at the same annealing temperature might be combined per two, if the probes of the second assay are differently labeled, e.g. with 5′-HEX-BHQ1–3′ and 5′-Cy5-BHQ2–3′. Although pooling techniques to detect allele ratios exist , it is difficult to precisely determine Wt/Vt allele ratios when working with more than two haploid drones or multiple diploid worker bees.
The use of genomic markers, such as the eight variants found by Broeckx et al. , in marker-assisted selection (MAS) is a promising method to accelerate the breeding progress on varroa resistance traits . Thus far, only protein markers have been successfully applied in breeding programs in honey bees . Compared to protein markers, genomic markers have the advantage of being independent of expression levels and are considered more stable . The prerequisite however is that the high recombination rate in honey bees  does not cause the breakdown of inter-allele linkages through repeated rounds of meiosis . Further research is ongoing to validate the effect of the eight SNVs on the SMR trait and thus of the applicability of genomic markers in MAS.
Supporting ongoing selective breeding programs with honey bees by applying different ‘omics tools opens new possibilities for better understanding underlying mechanisms and unrolling marker-assisted selection programs . The qPCR assays described in this paper neatly fits in with these future perspectives as it provides a novel laboratory based detection method to genotype honey bee colonies for the presence of the SMR trait.
Ninety-two DNA samples, used to analyze the allelic frequency of the eight SNVs in the Belgian honey bee population via Sanger sequencing in Broeckx et al. , were reused to optimize and validate the qPCR assays. They were isolated from two individual worker bees of the subspecies A. m. carnica from 46 different colonies located throughout Flanders, the northern part of Belgium.
A qPCR genotyping assay with dual-labeled probes was designed for each of the eight SNVs following the strategy described by Van Poucke et al. . Depending on the amplicon sequence, probes were designed on the forward or the reverse strand. Wild type (Wt) probes were labeled with 5′-FAM and 3’-BHQ1, variant type (Vt) probes with 5′-TexasRed and 3′-BHQ2.
Assay optimization and validation
The assays were performed in 10 μl containing 10x KEY buffer, 500 nM of each primer/probe, 800 μM dNTPs, 0.5 U TEMPase Hot Start DNA Polymerase (VWR) and 2 ng DNA on the CFX96 Touch Real-Time PCR Detection System (Bio-Rad). Thermal cycling conditions comprised 1 cycle of 14′40″ at 95 °C (activation Hot Start Polymerase and denaturation DNA) followed by 40 cycles of 20″ at 95 °C (denaturation DNA) and 40″ at the assay specific combined annealing/elongation/signal detection temperature (Table 1). Specific amplicon generation was checked by evaluating the PCR products using agarose gel electrophoresis. Optimal annealing temperature (Ta) was determined by performing gradient PCR and assessing probe specific signals on Wt/Wt (wild type homozygote), Wt/Vt (heterozygote) and Vt/Vt (variant type homozygote) samples. No template controls (NTC) were included to account for possible contaminations. The Sanger sequenced samples were used for validation and also checked for additional SNVs in the primer and probe binding sites. If present, they were evaluated for their influence on the result. For the assays genotyping SNV 4 and SNV 5 additional primers/probes had to be included for correct genotyping.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Suppressed mite reproduction
Single nucleotide variant
Rosenkranz P, Aumeier P, Ziegelmann B. Biology and control of Varroa destructor. J Invertebr Pathol. 2010;103(SUPPL. 1):S96–119. https://doi.org/10.1016/j.jip.2009.07.016.
Büchler R, Berg S, Le Conte Y. Breeding for resistance to Varroa destructor in Europe. Apidologie. 2010;41(3):393–408. https://doi.org/10.1051/apido/2010011.
Fries I, Hansen H, Imdorf A, Rosenkranz P. Swarming in honey bees (Apis mellifera) and Varroa destructor population development in Sweden. Apidologie. 2003;34:389–98.
Kruitwagen A, van Langevelde F, van Dooremalen C, Blacquière T. Naturally selected honey bee (Apis mellifera) colonies resistant to Varroa destructor do not groom more intensively. J Apic Res. 2017;56(4):354–65. https://doi.org/10.1080/00218839.2017.1329797.
Ramsey SD, Ochoa R, Bauchan G, Gulbronson C, Mowery JD, Cohen A, et al. Varroa destructor feeds primarily on honey bee fat body tissue and not hemolymph. Proc Natl Acad Sci U S A. 2019;116(5):1792–801. https://doi.org/10.1073/pnas.1818371116.
Gisder S, Möckel N, Eisenhardt D, Genersch E. In vivo evolution of viral virulence: switching of deformed wing virus between hosts results in virulence changes and sequence shifts. Environ Microbiol. 2018;20(12):4612–28. https://doi.org/10.1111/1462-2920.14481.
Locke B, Le Conte Y, Crauser D, Fries I. Host adaptations reduce the reproductive success of Varroa destructor in two distinct European honey bee populations. Ecol Evol. 2012;2(6):1144–50. https://doi.org/10.1002/ece3.248.
Johnson RM, Dahlgren L, Siegfried BD, Ellis MD. Acaricide, fungicide and drug interactions in honey bees (Apis mellifera). PLoS One. 2013;8:e54092.
Beaurepaire AL, Krieger KJ, Moritz RFA. Seasonal cycle of inbreeding and recombination of the parasitic mite Varroa destructor in honeybee colonies and its implications for the selection of acaricide resistance. Infect Genet Evol. 2017;50:49–54. https://doi.org/10.1016/j.meegid.2017.02.011.
Neumann P, Blacquière T. The Darwin cure for apiculture? Natural selection and managed honeybee health. Evol Appl. 2017;10(3):226–30. https://doi.org/10.1111/eva.12448.
Harbo JR, Hoopingarner RA. Honey bees (Hymenoptera: Apidae) in the United States that express resistance to Varroa jacobsoni (Mesostigmata: Varroidae). J Econ Entomol. 1997;90(4):893–8. https://doi.org/10.1093/jee/90.4.893.
Bienefeld K. Breeding success or genetic diversity in honey bees? Bee World. 2016;93(2):40–4. https://doi.org/10.1080/0005772X.2016.1227547.
Broeckx BJG, De Smet L, Blacquière T, Maebe K, Khalenkow M, Van Poucke M, et al. Honey bee predisposition of resistance to ubiquitous mite infestations. Sci Rep. 2019;9:1–12.
Camazine S. Differential reproduction of the mite, Varroa jacobsoni (Mesostigmata: Varroidae), on Africanized and European honey bees (Hymenoptera: Apidae). Ann Entomol Soc Am. 1986;79(5):801–3. https://doi.org/10.1093/aesa/79.5.801.
Conlon BH, Aurori A, Giurgiu AI, Kefuss J, Dezmirean DS, Moritz RFA, et al. A gene for resistance to the Varroa mite (Acari) in honey bee (Apis mellifera) pupae. Mol Ecol. 2019;28(12):2958–66. https://doi.org/10.1111/mec.15080.
Techer MA, Rane RV, Grau ML, Roberts JMK, Sullivan ST, Liachko I, et al. Divergent evolutionary trajectories following speciation in two ectoparasitic honey bee mites. Commun Biol. 2019;2(1):357. https://doi.org/10.1038/s42003-019-0606-0.
Beaurepaire A, Sann C, Arredondo D, Mondet F, Le Conte Y. Behavioral genetics of the interactions between Apis mellifera and Varroa destructor. Insects. 2019;10(9):299. https://doi.org/10.3390/insects10090299.
Mondet F, Parejo M, Meixner MD, Costa C, Kryger P, Bigio G, et al. Evaluation of suppressed mite reproduction ( SMR ) reveals potential for Varroa resistance in European. Insects. 2020;11:1–17.
Mondet F, Beaurepaire A, McAfee A, Locke B, Alaux C, Blanchard S, et al. Honey bee survival mechanisms against the parasite Varroa destructor: a systematic review of phenotypic and genomic research efforts. Int J Parasitol. 2020;50(6-7):433–47. https://doi.org/10.1016/j.ijpara.2020.03.005.
Weinstock G, Robinson G, Gibbs R. Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006;443:931–49.
Wallberg A, Glémin S, Webster MT. Extreme recombination frequencies shape genome variation and evolution in the honeybee, Apis mellifera. PLoS Genet. 2015;11:1–27.
Van Poucke M, Vandesompele J, Mattheeuws M, Van Zeveren A, Peelman LJ. A dual fluorescent multiprobe assay for prion protein genotyping in sheep. BMC Infect Dis. 2005;5(1):13. https://doi.org/10.1186/1471-2334-5-13.
Van Poucke M, Van Zeveren A, Peelman LJ. Combined FAM-labeled TaqMan probe detection and SYBR green I melting curve analysis in multiprobe qPCR genotyping assays. Biotechniques. 2012;52:81–5.
Chen Y, Bogema DR, Barchia IM, Herron GA. Quantification of the pirimicarb resistance allele frequency in pooled cotton aphid (Aphis gossypii glover) samples by TaqMan SNP genotyping assay. PLoS One. 2014;9(3):e91104. https://doi.org/10.1371/journal.pone.0091104.
Rinderer TE, Harris JW, Hunt GJ, De Guzman LI. Breeding for resistance to Varroa destructor in North America. Apidologie. 2010;41(3):409–24. https://doi.org/10.1051/apido/2010015.
Guarna MM, Hoover SE, Huxter E, Higo H, Moon KM, Domanski D, et al. Peptide biomarkers used for the selective breeding of a complex polygenic trait in honey bees. Sci Rep. 2017;7:1–10.
Liu H. Causes and consequences of crossing-over evidenced via a high-resolution recombinational landscape of the honey bee. Genome Biol. 2015;16:1–16.
We wish to thank local Flemish (Belgium) beekeepers for supplying bee samples, and Linda Impe, Dominique Vander Donckt and Ruben Van Gansbeke for excellent technical assistance.
This work was supported by FOD Volksgezondheid, Veiligheid van de Voedselketen en Leefmilieu (Cel Contractueel Onderzoek) RT19/6336 and by the European Union’s Horizon 2020 research and innovation programme under grant agreement No 817622 (B-GOOD).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Claeys Boúúaert, D., Van Poucke, M., De Smet, L. et al. qPCR assays with dual-labeled probes for genotyping honey bee variants associated with varroa resistance. BMC Vet Res 17, 179 (2021). https://doi.org/10.1186/s12917-021-02886-x
- Honey bee
- Varroa destructor
- Varroa resistance
- Suppressed mite reproduction
- High-throughput DNA test