Allele Frequencies of Genetic Variants Associated with Varroa Drone Brood Resistance (DBR) in Apis mellifera Subspecies across the European Continent.

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.

Population genetics of ectoparasitic mites suggest arms race with honeybee hosts.

The ectoparasitic mite, Varroa destructor, is the most severe biotic threat to honeybees (Apis mellifera) globally, usually causing colony death within a few years without treatments. While it is known that a few A. mellifera populations survive mite infestations by means of natural selection, the possible role of mite adaptations remains unclear. To investigate potential changes in mite populations in response to host adaptations, the genetic structure of V. destructor in the mite-resistant A. mellifera population on Gotland, Sweden, was studied. Spatio-temporal genetic changes were assessed by comparing mites collected in these colonies, as well as from neighboring mite-susceptible colonies, in historic (2009) and current (2017/2018) samples. The results show significant changes in the genetic structure of the mite populations during the time frame of this study. These changes were more pronounced in the V. destructor population infesting the mite-resistant honeybee colonies than in the mite-susceptible colonies. These results suggest that V. destructor populations are reciprocating, in a coevolutionary arms race, to the selection pressure induced by their honeybee host. Our data reveal exciting new insights into host-parasite interactions between A. mellifera and its major parasite.

Association of Varroa destructor females in multiply infested cells of the honeybee Apis mellifera.

The genetic diversity of Varroa destructor (Anderson & Trueman) is limited outside its natural range due to population bottlenecks and its propensity to inbreed. In light of the arms race between V. destructor and its honeybee (Apis mellifera L.) host, any mechanism enhancing population admixture of the mite may be favored. One way that admixture can occur is when two genetically dissimilar mites coinvade a brood cell, with the progeny of the foundresses admixing. We determined the relatedness of 393 pairs of V. destructor foundresses, each pair collected from a single bee brood cell (n = five colonies). We used six microsatellites to identify the genotypes of mites coinvading a cell and calculated the frequency of pairs with different or the same genotypes. We found no deviation from random coinvasion, but the frequency of cells infested by mites with different genotypes was high. This rate of recombination, coupled with a high transmission rate of mites, homogenized the allelic pool of mites within the apiary.

Seasonal cycle of inbreeding and recombination of the parasitic mite Varroa destructor in honeybee colonies and its implications for the selection of acaricide resistance.

Varroa destructor is the most devastating parasite of the Western honeybee, Apis mellifera. In the light of the arm race opposing the host and its parasite, the population dynamics and genetic diversity of these organisms are key parameters. However, the life cycle of V. destructor is characterized by extreme inbreeding due to full sibling mating in the host brood cells. We here present an equation reflecting the evolution of inbreeding in such a clonal system, and compare our predictions with empirical data based on the analysis of seven microsatellite markers. This comparison revealed that the mites perform essentially incestuous mating in the beginning of the brood season. However, this pattern changes with the development of mite infestation. Despite the fact that the overall level of genetic diversity of the mites remained low through the season, multiple inbred lineages were identified in the mites we sampled in June. As a response to the decrease of brood availability and the increase of the parasite population in parallel in the colonies, these lineages recombined towards the end of the season as mites co-infest brood cells. Our results suggest that the ratio of the number of mite per brood cell in the colony determines the genetic structure of the populations of V. destructor. This intracolonial population dynamics has great relevance for the selection of acaricide resistance in V. destructor. If chemical treatments occur before the recombination phase, inbreeding will greatly enhance the fixation of resistance alleles at the colony level.

Host Specificity in the Honeybee Parasitic Mite, Varroa spp. in Apis mellifera and Apis cerana.

The ectoparasitic mite Varroa destructor is a major global threat to the Western honeybee Apis mellifera. This mite was originally a parasite of A. cerana in Asia but managed to spill over into colonies of A. mellifera which had been introduced to this continent for honey production. To date, only two almost clonal types of V. destructor from Korea and Japan have been detected in A. mellifera colonies. However, since both A. mellifera and A. cerana colonies are kept in close proximity throughout Asia, not only new spill overs but also spill backs of highly virulent types may be possible, with unpredictable consequences for both honeybee species. We studied the dispersal and hybridisation potential of Varroa from sympatric colonies of the two hosts in Northern Vietnam and the Philippines using mitochondrial and microsatellite DNA markers. We found a very distinct mtDNA haplotype equally invading both A. mellifera and A. cerana in the Philippines. In contrast, we observed a complete reproductive isolation of various Vietnamese Varroa populations in A. mellifera and A. cerana colonies even if kept in the same apiaries. In light of this variance in host specificity, the adaptation of the mite to its hosts seems to have generated much more genetic diversity than previously recognised and the Varroa species complex may include substantial cryptic speciation.

Extensive population admixture on drone congregation areas of the giant honeybee, Apis dorsata (Fabricius, 1793).

The giant honeybee Apis dorsata often forms dense colony aggregations which can include up to 200 often closely related nests in the same location, setting the stage for inbred matings. Yet, like in all other Apis species, A. dorsata queens mate in mid-air on lek like drone congregation areas (DCAs) where large numbers of males gather in flight. We here report how the drone composition of A. dorsata DCAs facilitates outbreeding, taking into the account both spatial (three DCAs) and temporal (subsequent sampling days) dynamics. We compared the drones' genotypes at ten microsatellite DNA markers with those of the queen genotypes of six drone-producing colonies located close to the DCAs (Tenom, Sabah, Malaysia). None of 430 sampled drones originated from any of these nearby colonies. Moreover, we estimated that 141 unidentified colonies were contributing to the three DCAs. Most of these colonies were participating multiple times in the different locations and/or during the consecutive days of sampling. The drones sampled in the DCAs could be attributed to six subpopulations. These were all admixed in all DCA samples, increasing the effective population size an order of magnitude and preventing matings between potentially related queens and drones.