SummaryThe recognition that the Dark European honey bee, Apis mellifera mellifera, is increasingly threatened in its native range has led to the establishment of conservation programmes and protected areas throughout western Europe. Previous molecular surveys showed that, despite management strategies to preserve the genetic integrity of A. m. mellifera, protected populations had a measurable component of their gene pool derived from commercial C-lineage honey bees. Here we used both sequence data from the tRNA leu -cox2 intergenic mtDNA region and a genome-wide scan, with over 1183 single nucleotide polymorphisms (SNPs), to assess genetic diversity and introgression levels in several protected populations of A. m. mellifera, which were then compared with samples collected from unprotected populations. MtDNA analysis of the protected populations revealed a single colony bearing a foreign haplotype, whereas SNPs showed varying levels of introgression ranging from virtually zero in Norway to about 14% in Denmark. Introgression overall was higher in unprotected (30%) than in protected populations (8%), and is reflected in larger SNP diversity levels of the former, although opposite diversity levels were observed for mtDNA. These results suggest that, despite controlled breeding, some protected populations still require adjustments to the management strategies to further purge foreign alleles, which can be identified by SNPs. 270Pinto et al.

The natural distribution of the honeybee (Apis mellifera L.) has been changed by humans in recent decades to such an extent that the formerly widest-spread European subspecies, Apis mellifera mellifera, is threatened by extinction through introgression from highly divergent commercial strains in large tracts of its range. Conservation efforts for A. m. mellifera are underway in multiple European countries requiring reliable and cost-efficient molecular tools to identify purebred colonies. Here, we developed four ancestry-informative SNP assays for high sample throughput genotyping using the iPLEX Mass Array system. Our customized assays were tested on DNA from individual and pooled, haploid and diploid honeybee samples extracted from different tissues using a diverse range of protocols. The assays had a high genotyping success rate and yielded accurate genotypes. Performance assessed against whole-genome data showed that individual assays behaved well, although the most accurate introgression estimates were obtained for the four assays combined (117 SNPs). The best compromise between accuracy and genotyping costs was achieved when combining two assays (62 SNPs). We provide a ready-to-use cost-effective tool for accurate molecular identification and estimation of introgression levels to more effectively monitor and manage A. m. mellifera conservatories.

Beekeeping activities, especially queen trading, have shaped the distribution of honey bee (Apis mellifera) subspecies in Europe, and have resulted in extensive introductions of two eastern European C-lineage subspecies (A. m. ligustica and A. m. carnica) into the native range of the M-lineage A. m. mellifera subspecies in Western Europe. As a consequence, replacement and gene flow between native and commercial populations have occurred at varying levels across western European populations. Genetic identification and introgression analysis using molecular markers is an important tool for management and conservation of honey bee subspecies. Previous studies have monitored introgression by using microsatellite, PCR-RFLP markers and most recently, high density assays using single nucleotide polymorphism (SNP) markers. While the latter are almost prohibitively expensive, the information gained to date can be exploited to create a reduced panel containing the most ancestry-informative markers (AIMs) for those purposes with very little loss of information. The objective of this study was to design reduced panels of AIMs to verify the origin of A. m. mellifera individuals and to provide accurate estimates of the level of C-lineage introgression into their genome. The discriminant power of the SNPs using a variety of metrics and approaches including the Weir & Cockerham’s FST, an FST-based outlier test, Delta, informativeness (In), and PCA was evaluated. This study shows that reduced AIMs panels assign individuals to the correct origin and calculates the admixture level with a high degree of accuracy. These panels provide an essential tool in Europe for genetic stock identification and estimation of admixture levels which can assist management strategies and monitor honey bee conservation programs.

A large-scale survey of the Iberian honey bee (Apis mellifera iberiensis) diversity patterns, using sequence data of the tRNA leu-cox2 mitochondrial DNA (mtDNA) region, demonstrates that earlier studies based on the Dra I test missed significant components of genetic variation. Based on results from this survey, existing haplotype names were revised and updated following a nomenclature system established earlier and extended herein for the intergenic region. A more complete picture of the complex diversity patterns of IHBs is revealed that includes 164 novel haplotypes, 113 belonging to lineage A and 51 to lineage M and within lineage A and 69 novel haplotypes that belong to sub-lineage A I , 13 to A II , and 31 to A III. Within lineage M, two novel haplotypes show a striking architecture with features of lineages A and M, which based on sequence comparisons and relationships among haplotypes are seemingly ancestral. These data expand our knowledge of the complex architecture of the tRNA leu-cox2 intergenic region in Apis mellifera and re-emphasizes the importance of Iberia as a source of honey bee mtDNA diversity. Iberian honey bee / Apis mellifera intermissa / ancestral haplotype M / Dra I test