The availability of powerful high-throughput genomic tools, combined with genome scans, has helped identifying genes and genetic changes responsible for environmental adaptation in many organisms, including the honeybee. Here, we resequenced 87 whole genomes of the honeybee native to Iberia and used conceptually different selection methods (Samβada, LFMM, PCAdapt, iHs) together with in sillico protein modelling to search for selection footprints along environmental gradients. We found 670 outlier SNPs, most of which associated with precipitation, longitude and latitude. Over 88.7% SNPs laid outside exons and there was a significant enrichment in regions adjacent to exons and UTRs. Enrichment was also detected in exonic regions. Furthermore, in silico protein modelling suggests that several non-synonymous SNPs are likely direct targets of selection, as they lead to amino acid replacement...

The tRNAleu-COX2 intergenic region of the mitochondrial DNA has been used for assessing diversity in honey bee (Apis mellifera L.) populations worldwide. However, differential mutation rates in different partitions of the mitogenome may produce incongruent results. In this study, we sequenced 123 mitogenomes of 7 subspecies from lineages A, M, and C. This allowed generating a comprehensive dataset to investigate the phylogenetic and phylogeographic congruence among the mitogenome, individual genes, and the tRNAleu-COX2 region. We showed that the diversity patterns inferred from the tRNAleu-COX2 marker are not fully paralleled by those obtained with the mitogenome and the individual genes; while the three lineages are supported by these, the African sub-lineages and the haplotypes are not. Thus, conclusions drawn from the tRNAleu-COX2 region need to be taken with caution and this marker m...

A large-scale survey of the Iberian honey bee (Apis mellifera iberiensis) diversity patterns, using sequence data of the tRNAleu-cox2 mitochondrial DNA (mtDNA) region, demonstrates that earlier studies based on the DraI 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 AI, 13 to AII, and 31 to AIII. 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 a...