Genetic Properties and Evolution of Asian Honey Bee Apis cerana ussuriensis from Primorsky Krai, Russia

Abstract Apis cerana ussuriensis Ilyasov et al., 2019 is the northernmost subspecies of the Asian honey bee A. cerana Fabricius, 1793, common in the forests of Primorsky krai and Khabarovsk krai as far as 47°54′ N. Genetic studies of this subspecies are of great interest for science and apiculture, since all its adaptive traits were formed under the influence of the natural environment without human interference. We sequenced and annotated the complete mitochondrial DNA (mtDNA) sequences of bees of subspecies Apis cerana ussuriensis Ilyasov et al., 2019 (GenBank accession number AP018450) from Primorsky krai and Apis cerana koreana Ilyasov et al., 2019 (AP018431) from South Korea, as well as six exons of the nuclear DNA (nDNA) vitellogenin VG E2–E7 gene of bee subspecies A. c. ussuriensis, A. c. koreana, A. c. japonica Radoszkowski, 1887, A. c. cerana, and A. c. indica Fabricius, 1798. Cluster analysis of the mtDNA and the nDNA VG gene sequences showed the division of bees into two groups, the southern subspecies A. c. indica and the northern subspecies A. c. ussuriensis, A. c. koreana, A. c. japonica, and A. c. cerana. On the basis of the genetic divergence, we showed that subspecies A. c. ussuriensis was genetically closer to subspecies A. c. japonica, A. c. koreana, and A. c. cerana than to subspecies A. c. indica. Values of genetic divergence (0.80–8.00%) and Jukes–Cantor genetic distance (0.005–0.100) for mtDNA and nDNA VG gene between subspecies A. c. ussuriensis, A. c. koreana, A. c. japonica, A. c. cerana, and A. c. indica are within the range of intraspecific differences between insect subspecies. The estimated time of the emergence of the A. cerana subspecies is from two to one million years ago.