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

Wang, Ah Rha;Jeong, Su Yeon;Jeong, Jun Seong;Kim, Seong Ryul;Choi, Yong Soo;Kim, Iksoo

  • Received : 2015.10.07
  • Accepted : 2015.11.03
  • Published : 2015.12.31


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 $tRNA^{Leu}$ 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.


Asian cavity-nesting honey;Apis cerana;introns;white gene;MRJP9;sequenced diversity


  1. Sihanuntavong D, Sittipraneed S, Klinbunga S (1999) Mitochondrial DNA diversity and population structure of the honey bee, Apis cerana, in Thailand. J Apic Res 38, 211-219.
  2. Smith DR, Hagen RH. (1996) The biogeography of Apis cerana as revealed by mitochondrial DNA sequence data. J Kans Entomol Soc 69, 249-310.
  3. Smith DR, Villafuerte L, Otis G, Palmer MR (2000) Biogeography of Apis cerana F. and A. nigrocincta Smith: Insights from mtDNA studies. Apidologie 31, 265-279.
  4. Smith DR, Warrit N, Hepburn HR (2004) Apis cerana from Myanmar (Burma): unusual distribution of mitochondrial lineages, Apidologie 35, 637-644.
  5. Swofford DL (2002) PAUP*. Phylogenetic Analysis Using Parsimony (*and other method) software ver. 4.10, Sinauer Associates. Sunderland, Massachusetts.
  6. Takahashi J, Yoshida T, Takagi T, Akimoto S, Wood KS, Deowanish S, Hepburn R, Nakamura S, Matsuka M (2007) Geographic variation in the Japanese islands of Apis cerana japonica and in A. cerana populations bordering its geographic range Apidologie 38, 335-340.
  7. Tan K, Warrit N, Smith DR (2007) Mitochondrial DNA diversity of Chinese Apis cerana. Apidologie 38, 238-246.
  8. Warrit N, Smith DR, Lekprayoon C (2006) Genetic subpopulations of Varroa mites and their Apis cerana hosts in Thailand. Apidologie 37, 19-30.
  9. Yu L, Luan PT, Jin W, Ryder OA, Chemnick LG, Davis HA, Zhang YP (2011) Phylogenetic utility of nuclear introns in interfamilial relationships of Caniformia (order Carnivora). Syst Biol 60, 175-187
  10. Zink RM, Barrowclough GF (2008) Mitochondrial DNA under siege in avian phylogeography. Mol Ecol 17, 2107-2121
  11. De la Rúa P, Simon UE, Tilde AC, Moritz RFA, Fuchs S (2000) MtDNA variation in Apis cerana populations from the Philippines. Heredity 84, 124-130.
  12. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10, 564-567.
  13. Katoh K, Misawa K, Kuma K, Miyata T (2002) MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res 30, 3059-3066.
  14. Kawakita A, Ascher JS, Sota T, Kato M, Roubik DW (2008) Phylogenetic analysis of the corbiculate bee tribes based on 12 nuclear protein-coding genes (Hymenoptera: Apoidea: Apidae). Apidologie 39, 163-175.
  15. Kim HK, Lee MY, Lee ML, Choi YS, Hong IP, Kim MS, Lee KG, Jin BR (2010) Genetic diversity of the Apis cerana in Korea by mitochondrial DNA and ITS ribosomal DNA sequence. J Apic 25, 9-16.
  16. Lee JY, Wang AR, Choi YS, Thapa R, Kwon HW and Kim I (2015). Mitochondrial DNA variations in Korean Apis cerana (Hymenoptera: Apidae) and development of another potential marker. Apidologie In Press.
  17. Lo N, Gloag RS, Anderson DL, Oldroyd BP (2010) A molecular phylogeny of the genus Apis suggests that the Giant Honey Bee of the Philippines, A. breviligula Maa, and the Plains Honey Bee of southern India, A. indica Fabricius, are valid species. Syst Entomol 35, 226-233.
  18. Peiren N, de Graaf DC, Vanrobaeys F, Danneels EL, Devreese B, Van Beeumen J, Jacobs FJ (2008) Proteomic analysis of the honey bee worker venom gland focusing on the mechanisms of protection against tissue damage. Toxicon 52, 72-83.
  19. Sang T (2002) Utility of low-copy nuclear gene sequences in plant phylogenetics. Crit Rev Biochem Mol Biol 37, 121-147.
  20. Cornuet JM, Garnery L, Solignac M (1991) Putative origin and function of the intergenic region between COI and COII of Apis mellifera L. mitochondrial DNA. Genetics 128, 393-403.
  21. Crozier RH, Crozier YC (1993) The mitochondrial genome of the honeybee Apis mellifera: complete sequence and genome organization. Genetics 133, 97-117.
  22. Deowanish S, Nakamura J, Matsuka M, Kimura K (1996) mtDNA variation among subspecies of Apis cerana using restriction fragment length polymorphism. Apidologie 27, 407-413.


Grant : 글로벌 식물 헬스케어 전문인력 양성 사업팀