In Korea, the Asian honey bee (Apis cerana) and the European honey bee (Apis mellifera) (Hymenoptera: Apidae) are the two most common honey bee species. These two closely related species are known to have different sensitivity levels to various insecticides due to millennia of exposure to different pests and pesticides. It is reported that A. cerana is known to be more sensitive to several insecticides, such as amitraz, fenitrothion, and fipronil, than A. melllifera. Multiple studies investigated toxicological responses and related CYPome in A. mellifera, but little is known in A. cerana. The goal of this study is to elucidate the underlying mechanism of different toxicological responses between two bee species, with an emphasis on cytochrome P450 (P450), a significant enzyme involved in metabolic activities. The differences in basal P450 expression patterns were investigated by comparing the relative expression levels of P450 orthologs in several dissected organisms of each species. To compare the sensitivity against major insecticides, lethal doses of major insecticides relevant to both honey bee species were assessed by topical and oral ingestion bioassays. The determined sublethal doses of insecticides were applied to honey bees, and the inducibility of P450s was investigated by comparing the expression patterns of multiple P450s. From these results, this study eventually attempts to compare the toxicological differences between two Apis species with differences in induced cytochrome P450 expression levels.

새롭게 육성된 낭충봉아부패병 저항성 신품종 토종벌(Apis cerana koreana) 과 기존 농가에서 관행적으로 사육되는 토종벌 사이의 형태학적 차 이를 육안으로 확연하게 구분하는 것은 어렵지만 본 연구에서는 신품종 토종벌(A. c. koreana) 을 기존 토종벌(A. c. koreana) 품종 및 계통 간 형태학적 비교를 통해 신품종 만의 특성을 결정할 수 있는 표현형 정보를 제공하였다. 신품종 토종벌(A. c. koreana)의 외부형질을 이용한 품종 특성은 22가지의 형태학적 특성을 기하학적, 형태학적 분석 방법을 적용하고 토종벌(A. c. koreana)의 로얄젤리 생산량, 일벌, 여왕벌, 수벌의 특성을 비교 분석하였다. 본 연구 결과, 신품종 토종벌(A. c. koreana)은 기존 토종벌과 앞날개의 길이에 차이를 보였으며, 중국의 동양종꿀벌(A. cerana)과 비교한 결과, 일벌은 몸무게, 혀의 길이, 앞날개의 길이 등의 값이 높았다. 또한, 신품종 토종벌(A. c. koreana)은 A. cerana indica 보다 두 가지 부위에서 형태학적인 차이를 보였다. 그리고 신품종 토종벌(A. c. koreana)은 로얄젤리를 다른 품종과 비교하여 많이 분비하여 봉군의 발육에 긍정적인 영향을 끼쳤다. 따라서 본 연구결과는 신규 육성 토종벌(A. cerana)에 대한 형태학적 분석 방법을 이용하여 품종을 분류하는데 도움이 될 것으로 기대한다.

To control an external parasitic mite, a honey bee line possessing high hygienic behavior (HHB) against an external parasitic mite, Varroa destructor, has been bred in South Korea and an assessment method has been necessitated to diagnose HHB line from the low hygienic behavior (LHB) line. Thus, in this study, we developed single nucleotide polymorphism (SNP) markers from whole genome sequencing of each 20 worker bees from HHB and LHB lines of A. mellifera ligustica (Hymenoptera: Apidae). An average of 319,445,977 sequence reads was mapped to the known A. mellifera reference genome (an average of 87.46%). In 2,316,128 and 3,266,756 SNPs from each HHB and LHB line, an average of 93.6% and was located in the intergenic spacers and introns, whereas, the remaining 6.4% was located in the genic region, respectively. Among them 20 SNPs that were fixed at each line possessing within-individual homozygosity were selected and each four SNPs were used to diagnose the two honey bee lines either by typical PCR-restriction fragment length polymorphism method or allele-specific PCR. The remaining six SNPs had the size difference, enabling relatively easy differentiation between the two honey bee lines on typical agarose gel and another remaining six SNPs only has sequence difference including SNP sites. Thus, these SNP markers can be used to diagnose the honey bee line with HHB from LHB line against V. destructor.

Bumblebee, Bombus ardens ardens (Apidae: Hymenopera), is an important resource for pollination that is most widely distributed in Korea. This study utilized microsatellite markers for investigation of genetic diversity and geographic relevance of the B. a. ardens populations in Korea. Through Next Generation Sequencing analysis, we identified 10 microsatellite markers and genotyped for 107 individuals of B. ardens collected from 10 populations. At each locus the number of alleles ranged from 10 to 23; the observed and expected heterozygosities ranged from 0.8909 to 1.0000 and 0.6641 to 0.8422, respectively; and inbreeding coefficient(FIS) ranged from –0.5053 to –0.0891. Significant deviation from the Hardy–Weinberg equilibrium was not observed at any locus. Population structure analysis indicated that there are three genetic groups in Korea with each Jeongseon and Ulleung-do composed of different gene pool from the remaining other populations. Similarly, Principal coordinates analysis also showed the same pattern. FST and RST analyses showed that each Jeongseon and Ulleung-do population had a significant genetic distance from other populations. Considering these results, genetic isolation of Ulleung-do may be explained by “Oceanic island” status and Jeongseon, which showed the positive FIS (0.069) and genetic isolation may be caused by its location on the east side of Baekdudaegan and by on-going inbreeding with a small population size.

We developed single nucleotide polymorphism (SNP) markers and are establishing diagnostic systems to distinguish disease resistance- and susceptible-strains of honey bees using the SNPs. For development of SNP markers, whole genome was sequenced each from 20 individuals of “disease resistance-strain” and “susceptible-strain” of Apis mellifera ligustica using the Illumina HiSeq 2000 sequencer. Approximately, 344 and 294 million sequence reads were mapped to the honeybee reference assembly (Amel_4.5) for each strain, respectively. Among the total 2,246,428 SNPs yielded, 33 were found to be fixed between the two strains with all homozygosity. Sixteen of them were casually amplified and sequenced from randomly selected each 10 individual of honey bees from each strain and presented strain specific SNPs. These ten SNPs were used to diagnose the two strains either by original size difference, caused by indel-accompanying SNP, typical PCR-RFLP, or AS PCR.

We sequenced complete mitochondrial genomes (mt genomes) of four Apis species, such as A. cerana, A. dorsata, A. laboriosa, and A. mellifera ligustica. The A. m. ligustica used in this study was originated from United States and selected for a high hygienic behavior during several years of breeding program. Considering current phylogenetic hypothesis of Apis the positive relationships between phylogeny and gene rearrangement was not supported. Phylogenetic reconstruction using available Apoidea species has shown the relationships ((Apidae + Colletidae) + Crabronidae). Within Apis, a strong support for three groups was obtained: A. florea and A. andreniformis as a basal group to the other A. cerana and A. mellifera group and A. dorsata and A. laboriosa group. The test to find out suitable single mt gene for the phylogeny at the genus level and below-genus level of Apis provided an identical topology from ND4L and ND6 in the ML analysis and ND3 and NF4L in the BI analysis to that from whole mt genome, with relatively strong support, indicating that these individual genes can potentially be utilized for within-species level phylogeny of Apis.

We sequenced 17,329 bp of mitochondrial genome (mitogenome) of the black dwarf honey bee, Apis andreniformis (Hymenoptera: Apidae), that lacks ~200 bp of the A+T-rich region for the completion of the genomic sequence. The gene arrangement of A. andreniformis mitogenome is identical to that of A. cerana. However, the genome contains 5 additional tRNALeu(CUN) located 4 copies between tRNAMet and tRNAGln, and 1 copy between tRNAGln and tRNAAla, along with the typical sets of genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs) including regular tRNALeu(CUN) and the A+T-rich region (at least 923 bp). Only 1 copy of tRNALeu(CUN) differed by 1 bp from other 4 copies of tRNALeu(CUN). Each additional tRNALeu(CUN) is followed by nearly identical 68-bp long repeat sequence (95.6% identity). All 13 protein coding genes have typical start codons found in insect mitochondrial PCGs (2 ATA, 9 ATT, and 2 ATG).

The bumblebee, Bombus ignitus (Hymenoptera: Apidae), is a valuable natural resource that is widely utilized for greenhouse pollination in South Korea. Understanding the magnitude of genetic diversity and geographic relationships is of fundamental importance for long term preservation and utilization. As a first step, we sequenced a partial COI gene of mitochondrial DNA (mtDNA) corresponding to the “DNA barcode” region and the complete internal transcribed spacer 2 (ITS2) of nuclear ribosomal DNA from 88 individuals collected in nine South Korean localities. The complete ITS2 sequences were longest among known insects, ranging in size from 2,034 bp ~ 2,052 bp, harboring two duplicated 112-bp long repeats. The 658-bp long mtDNA sequences provided only six haplotypes with a maximum sequence divergence of 0.61% (4 bp), whereas the ITS sequences provided 84 sequence types with a maximum sequence divergence of 1.02% (21 sites). The combination of the current COI data with those of published data suggest that the B. ignitus in South Korea and China are genetically a large group, but those in Japan can be roughly separated into another group. Overall, a very high per generation migration ratio, a very low level of genetic fixation, and no discernable hierarchical population were found to exist among the South Korean populations of B. ignitus, which suggests panmixia. This finding is consistent with our understanding of the dispersal capability of the species.