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.

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

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human-mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)-based genome scan in honey bees collected across an environmental gradient in Iberia and used four FST -based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross-validated by four FST -based methods, and 8.9% (34 of 383) were cross-validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human-mediated environmental changes.