Intron sequence diversity of the asian cavity-nesting honey bee, Apis cerana (Hymenoptera: Apidae)

The Asian cavity-nesting honeybee, Apis cerana (Hymenoptera: Apidae), has been extensively studied for its biogeography and genetic diversity, but the molecules utilized in past studies were mainly ~90 bp long mitochondrial non-coding sequences, located between tRNALeu and COII. Thus, additional molecular markers may enrich our understanding of the biogeography and genetic diversity of this valuable bee species. In this study, we reviewed the public genome database to find introns of cDNA sequences, with the assumption that these introns may have less evolutionary constraints. The six introns selected were subjected to preliminary tests. Thereafter, two introns, titled White gene and MRJP9 gene, were selected. Sequencing of 552 clones from 184 individual bees showed a total of 222 and 141 sequence types in the White gene and MRJP9 gene introns, respectively. The sequence divergence ranged from 0.6% to 7.9% and from 0.26% to 17.6% in the White gene and the MRJP9 introns, respectively, indicating higher sequence divergence in both introns. Analysis of population genetic diversity for 16 populations originating from Korea, China, Vietnam, and Thailand shows that nucleotide diversity (π) ranges from 0.003117 to 0.025837 and from 0.016541 to 0.052468 in the White gene and MRJP9 introns, respectively. The highest π was found in a Vietnamese population for both intron sequences, whereas the nine Korean populations showed moderate to low sequence divergence. Considering the variability and diversity, these intron sequences can be useful as non-mitochondrial DNA-based molecular markers for future studies of population genetics.