Abstract
Implementation of marker-assisted selection (MAS) in modern beekeeping would improve sustainability, especially in breeding programs aiming for resilience against the parasitic mite Varroa destructor. Selecting honey bee colonies for natural resistance traits, such as brood-intrinsic suppression of varroa mite reproduction, reduces the use of chemical acaricides while respecting local adaptation. In 2019, eight genomic variants associated with varroa non-reproduction in drone brood were discovered in a single colony from the Amsterdam Water Dune population in the Netherlands. Recently, a new study tested the applicability of these eight genetic variants for the same phenotype on a population-wide scale in Flanders, Belgium. As the properties of some variants varied between the two studies, one hypothesized that the difference in genetic ancestry of the sampled colonies may underly these contribution shifts. In order to frame this, we determined the allele frequencies of the eight genetic variants in more than 360 Apis mellifera colonies across the European continent and found that variant type allele frequencies of these variants are primarily related to the A. mellifera subspecies or phylogenetic honey bee lineage. Our results confirm that population-specific genetic markers should always be evaluated in a new population prior to using them in MAS programs.