New COI-COII mtDNA Region Haplotypes in the Endemic Honey Bees Apis mellifera intermissa and Apis mellifera sahariensis (Hymenoptera: Apidae) in Algeria.

The practice of beekeeping in Algeria is of great cultural, social, and economic importance. However, the importation of non-local subspecies reported by beekeepers has disrupted the natural geographical distribution area and the genetic diversity of the native honey bees. To assess the genetic diversity of A. m. intermissa and A. m. sahariensis, and their relationships with African and European subspecies, the COI-COII intergenic region was analyzed in 335 individuals, 68 sampled in Algeria, 71 in Europe, Madagascar, and the South West Indian Ocean archipelagos, and 196 sequences recovered from GenBank. The results show the presence of the A lineage exclusively in Algerian samples with the identification of 24 haplotypes of which 16 are described for the first time. These haplotypes were found to be shared by both subspecies, with A74 being the most common haplotype in the population studied. The sequence comparison indicates the existence of three polymorphisms of the COI-COII marker: P0Q, P0QQ, and P0QQQ. One new haplotype was identified in the M lineage in samples from France. No evidence of genetic introgression within the Algerian honey bee population was detected. These data enhance our knowledge of the genetic diversity and emphasize the importance of protecting these local subspecies.

Natural clines and human management impact the genetic structure of Algerian honey bee populations.

BACKGROUND: The Algerian honey bee population is composed of two described subspecies A. m. intermissa and A. m. sahariensis, of which little is known regarding population genomics, both in terms of genetic differentiation and of possible contamination by exogenous stock. Moreover, the phenotypic differences between the two subspecies are expected to translate into genetic differences and possible adaptation to heat and drought in A. m. sahariensis. To shed light on the structure of this population and to integrate these two subspecies in the growing dataset of available haploid drone sequences, we performed whole-genome sequencing of 151 haploid drones. RESULTS: Integrated analysis of our drone sequences with a similar dataset of European reference populations did not detect any significant admixture in the Algerian honey bees. Interestingly, most of the genetic variation was not found between the A. m. intermissa and A. m. sahariensis subspecies; instead, two main genetic clusters were found along an East-West axis. We found that the correlation between genetic and geographic distances was higher in the Western cluster and that close-family relationships were mostly detected in the Eastern cluster, sometimes at long distances. In addition, we selected a panel of 96 ancestry-informative markers to decide whether a sampled bee is Algerian or not, and tested this panel in simulated cases of admixture. CONCLUSIONS: The differences between the two main genetic clusters suggest differential breeding management between eastern and western Algeria, with greater exchange of genetic material over long distances in the east. The lack of detected admixture events suggests that, unlike what is seen in many places worldwide, imports of queens from foreign countries do not seem to have occurred on a large scale in Algeria, a finding that is relevant for conservation purposes. In addition, the proposed panel of 96 markers was found effective to distinguish Algerian from European honey bees. Therefore, we conclude that applying this approach to other taxa is promising, in particular when genetic differentiation is difficult to capture.