The spotted-wing drosophila (SWD), Drosophila suzukii (Diptera: Drosophilidae), was originally observed in a few Asian countries, but is now found even in North America and Europe. Genetic information on geographic variation and relationship may broaden our understanding of origin and migration. As a first step, in this study, a portion of mitochondrial COI gene was sequenced to understand genetic relationship and diversity in Korea. Sequencing of 104 individuals provided 57 haplotypes, with the maximum sequence divergence of 1.5%, suggesting high haplotype diversity and moderate sequence divergence. Comparison to GenBankregistered D. suzukii haplotypes (possibly from Spain, Portugal, USA) has shown 100% sequence identity to most of the haplotypes found in this study, but two USA sequences were found to be independent haplotypes, with the sequence divergence ranging from 0.5% ~ 1.4% from our samples in the 553-bp comparison. Phylogenetically, no separable group was found, but, population genetically, the only Chinese population, Sandong, was significantly differentiated (p < 0.05) from all Korean populations, without sharing any haplotype. Among 28 pairwise comparisons of Korean populations only two comparisons showed a significant genetic differentiation, indicating that no population in Korea is completely isolated. Geographically, one haplotype (SWDBA06) was relatively widespread (five among nine localities) and a few haplotypes were found in more than one locality, but most haplotypes were restricted in a locality as a single individual. Overall, high rate of per generation female migration (Nm = 0.75 ~ infinite) and low level of geographic separation (FST=0~0.40) among localities were characteristic. Current data is limited mainly to Korean localities, thus, an expanded study may provide further scrutinized analysis for the fly.

Worldwide studies on Apis cerana variation for biogeography and genetic diversity depended largely on a 86~93 bp-long mitochondrial non-coding region (internal spacer region) located between tRNALeu and COII (named as NC2), possibly due to higher variability among available markers. In order to incorporate the A. cerana occurring in South Korea into world extensive data, we also sequenced the NC2 from 118 A. cerana samples collected over nine Korean localities and 66 A. cerana samples over seven Asian localities, such as China, Vietnam, and Thailand. These data were combined with preexisting world data to scrutinize genetic relationships of A. cerana in South Korea to outside distributional range. Sequencing of 184 samples provided a total of ten haplotypes: five from Korea, six from China, one from Vietnam, and two from Thailand. Among them eight were new, whereas two were previously reported ones. Phylogenetic analysis of A. cerana NC2 haplotypes so far found including ours has confirmed the presence of four major groups of A. cerana (Asian mainland group, Sundaland group, Palawan group, and Luzon-Mindahnao group) and all haplotypes found in this study also were included in the Asian mainland group. In order to find further variable regions that can be used as sequence-based marker several mitochondrial non-coding regions and nuclear intron regions are in the middle of testing.

Honeybees (Apis mellifera) adapted themselves to different geographical and climatical conditions since they had been introduced in Korea. Beekeepers have tried to breed valuable lineages with artificial insemination or conventional mating techniques. However, evaluation of breeding resultants still relies on timeconsuming observation data. Genetical characterization of breeds has proven its usefulness to preserve genetic resources of livestock. In recent years, microsatellites are most commonly used to evaluate population structures and diversities of living organisms in that the characteristics of locus specificity, rich polymorphism, abundant and random distribution over the genome, and their co-dominant inheritance. Determining classic genetic distances using neutral, highly polymorphic microsatellite markers is a reliable method to investigate genetic relationships and breed differentiation. This methodology can be used to establish preservation priorities for livestock breeds. The final aim of this study was to develop potent markers for assessing genetic structure of lineages after artificial insemination. In this study, the genetic structure of ten microsatellite markers were sequenced or analysed with polymerase chain reaction for eleven European honeybee populations. The results may help to develop reliable microsatellite markers for more efficient preservation strategies of valuable honeybee breeds.

The rice leaf roller, Cnaphalocroci smedinalis Guenée (Lepidoptera: Pyralidae) is a leaf-feeding pest of rice world-widely distributed. For better understanding of the pest insect, geographic sequence variation of the species were performed using the mitochondrial A+T-rich region, with the samples collected from seven Korean and six Chinese localities. A total of 94 haplotypes obtained from 187 individuals showed the length variation, ranging from 339 bp to 348 bp. The maximum divergence of 4.57% appears to evidence a substantial sequence variation, indicating the applicability of this molecular marker to the study of geographic variation. Overall, a high per generation migration ratio (Nm = 3.67742 ~ infinite), a low level of genetic fixation (FST = 0 ~ 0.11969), and no discernable isolated population were noted in the most C. medinalis populations. AMOVA analysis to find out allocation of genetic variability of C. medinalis populations has shown allocation of majority of variation to the within-population, rather than among-populations and between-region, suggesting that the C. medinalis populations in both China and Korea are largely well connected. This result is consistent with current knowledge of the dispersal ability of the species. The structure analysis of the A+T-rich region has shown that the typical structural elements found in other lepidopteran insects also is well preserved in the C. medinalis A+T-rich region (i.e., a poly-T stretch and a microsatellite-like A/T repeat).

The bumblebee, Bombus ignitus (Hymenoptera: Apidae), is a valuable natural resource that is one of the most notably utilized for greenhouse pollination in Korea. In order to understand the nature of genetic relationships, gene flow, and population structure of the species we sequenced a partial COI gene of mitochondrial DNA (mtDNA) corresponding to “animal barcode” region and the complete internal transcribed spacer 2 (ITS2) of the nuclear ribosomal DNA (nrDNA) collected from Korean localities. Although the 658-bp long mtDNA sequence provided only six haplotypes with the maximum sequence divergence of 0.61% (4 bp), the ITS sequences provided 84 sequence types with the maximum sequence divergence of 1.02% (21 sites), confirming better applicability of the ITS sequences to the study of intraspecific variation. The complete ITS2 sequences of B. ignitus were shown to be longest among known insects, ranging in size from 2,034 bp ~ 2,052 bp, harboring two duplicated repeats. Overall, a very high per generation migration ratio, a very low level of genetic fixation, and no discernable hierarchical population/ population group were noted to exist among populations of B. ignitus on the basis of both molecules, thus suggesting that the B. ignitus populations on the Korean peninsula are panmictic, which is consistent with our understanding of the dispersal capability of the species

In an effort to gain a better understanding of the nature of the population genetic structure of the pest insect, Spodoptera litura (Lepidoptera: Noctuidae), tobacco cutworms were collected from six Korean and five Chinese localities and their mitochondrial A+T-rich region and nuclear internal transcribed spacer 2 (ITS2) regions were cloned and sequenced. A total of 106 A+T-rich region haplotypes and 92 ITS2 sequence types were obtained from 158 individuals. Compared to the preliminary mitochondrial COI gene sequence data that provided all identical sequence types in all samples, the A+T-rich region and nuclear ITS2 are highly variable. Overall, a low level of genetic fixation in the A+T-rich region (FST = 0 ~ 0.02965) and nuclear ITS2 (FST = 0 ~ 0.34491), and no discernable isolated population was noted to exist among most S. litura populations. The presence of a majority of within-population variations, rather than variations among populations or between China and Korea, indicates that the S. litura populations are composed of heterogeneous individuals. The fixation index at hierarchical level of among regions shows no significance, suggesting that the S. litura populations in both Korea and China are profoundly interrelated with one another. This finding is in accord with the current knowledge that S. litura has sufficient flight capacity for dispersal.

In order to understand the nature of population genetic structure of the notorious pest insect, Spodoptera litura (Lepidoptera: Noctuidae), a total of 159 individuals of the tobacco cutworm collected from six Korean and five Chinese localities were cloned and sequenced their mitochondrial A+T-rich region. A total of 90 haplotypes ranged 324 bp ~ 372 bp in length and 0.30% ~ 4.85% in sequence divergence. Compared to the preliminary mitochondrial COI gene sequence data that provided nearly identical sequences in all samples, the A+T-rich region provided substantially high sequence variation, validating the applicability of this molecule to the study of intraspecific variation. Overall, a very high per generation migration ratio (Nm = 5.91 ~ infinite), a very low level of genetic fixation (FST = 0 ~ 0.077796), and no discernable isolated population were noted to exist among populations of S. litura, although some significant genetic differentiations were found between several pairs of populations. These results collectively suggest that the S. litura populations in both Korea and China are very well interrelated to each other. The structural analysis revealed that the S. litura A+T-rich region harbored two stretches of the [TA(A)]n sequence, the poly-A stretch, and a stem-and-loop structure that are well found in other lepidopteran A+T-rich region. Furthermore, the region harbors the 13 ~ 19-bp long T stretch and 5-bp long motif that have been suggested to function as a possible recognition site for the initiation of replication of the minor strand of mtDNA also were found.

The bumblebee, Bombus ignitus (Hymenoptera: Apidae), is a valuable natural resource that is one of the most notably utilized for greenhouse pollination in Korea. In order to understand the nature of genetic relationships, gene flow, and population structure of the species we sequenced a partial COI gene of mitochondrial DNA (mtDNA) corresponding to “animal barcode” region and the complete internal transcribed spacer 2 (ITS2) of the nuclear ribosomal DNA (nrDNA) with the individuals collected from Korean localities. Although the 658-bp long mtDNA sequence provided only six haplotypes with the maximum sequence divergence of 0.456% (3 bp, from 91 individuals), but the ITS sequences provided 33 sequence types with the maximum sequence divergence of 0.78% (16 bp, from 35 individuals), confirming better applicability of the ITS sequences to the study of intraspecific variation. The ITS2 sequences of B. ignitus were shown to be the longest among known insects, ranging in size from 2,034 bp ~ 2,045 bp and harbored two duplicated repeats, indicating unusual structure of B. ignitus ITS2 sequences compared with other insect ITS sequences. Overall, a very high per generation migration ratio, a very low level of genetic fixation, and no discernable hierarchical population/population group were noted to exist among populations of B. ignitus on the basis of both molecules, thus suggesting that the B. ignitus populations on the Korean peninsula are panmictic, which is consistent with our understanding of the dispersal capability.

The bumblebees, Bombus species are valuable natural resources being utilized for greenhouse pollination. Low level of genetic variation of Bombus species has been reported previously. In this study, we sequenced complete internal transcribed spacer 2 (ITS2) of the nuclear rDNA from 100 individuals of B. ardens collected from seven localities in Korean peninsula. The ITS2 sequence of B. ardens is longest, ranging in size from 1,940 bp-1,954 bp among known in insects, which ranges approximately from 241 bp-1,728 bp. The ITS2 sequences have ~51% of C+C content and contain each two 27 bp repeats, 20 bp repeats, 33 bp repeats, and 34 bp repeats at the beginning. Such repeats were not found in other insects. Uncorrected pairwise distance among 96 sequences that were obtained from 100 individuals revealed a maximum sequence divergence of 1.03%. Genetic diversity (π) of each population ranged from 0.007801 to 0.009627, and the lowest diversity was obtained from islet population of Ulleungdo, indicating possibly small, isolation of the population. Significant level of genetic distance was only found when Ulleungdo population was compared to two other mainland populations. Except for this, overall, a very high rate of per generation migration ratio (Nm=7.1-infinite) and a very low level of genetic fixation (FST=0-0.06546) were detected between pairs of localities. Analysis of hierarchical relationships among localities consistently revealed no statistically significant structure among populations. Taken these together, the B. ardens populations on the Korean peninsula are panmictic this is consistent with our understanding of the dispersal capability.

The bumblebee, Bombus ardens, is a valuable natural resource, and is most notably utilized for greenhouse pollination. In order to gain a greater understanding of the population genetic structure and the genetic diversity of this species, we sequenced a portion of the mitochondrial COI gene corresponding to the “DNA Barcode” region (658 bp) from 160 individuals collected over 15 Korean localities. Uncorrected pairwise distances among the eight haplotypes suggested low intraspecific genetic diversity, with a maximum sequence divergence of 0.3%. Such a low level of intraspecific genetic diversity was further reflected in local populations, particularly to islet populations, such as those of Youngheungdo, Jakyakdo, and Ulleugdo, which had zero genetic diversity. Geographically, one haplotype (BARBA01) was widespread and dominant, with a frequency of 90.6% (145 among 160 individuals). Other haplotypes were restricted to one to three localities and had low frequency. Overall, a very high rate of per generation female migration ratio (Nm= 4.6 to infinite) and a very low level of genetic fixation (FST= 0 to 0.099) were detected between pairs of localities, suggesting that the B. ardens population sonthe Korean peninsula are panmictic, which is consistent with our understanding of their dispersal capability.

We analyzed a portion of mitochondrial COI gene sequences (658 bp) to investigate the genetic diversity and geographic variation of the swallowtail butterfly, Papilioxuthus L., and the cabbage butterfly, Pieris rapae (Lepidoptera: Papilionidae). P. xuthus showed a moderate level of sequence divergence (0.91% at maximum) in 15 haplotypes, whereas P. rapae showed a moderate to high level of sequence divergence (1.67% at maximum) in 30 haplotypes, compared with other relevant studies. Analyses of population genetic structure showed that most populations are not genetically differentiated in both species. The distribution pattern of both species appears to be consistent with category IV of the phylogeographic pattern sensu Avise (Avise et al. 1987): a phylogenetic continuity, an absence of regional isolation of mtDNA clones, and extensive distribution of close clones. The observed pattern of genetic diversity and geographic variation of the two butterfly species seems to reflect the abundant habitats, abundant host plants, and flying abilities in connection with the lack of historical biogeographic barriers.

The mason bee (Osmia cornifrons Radoszkowsky) is an excellent pollinator of apple. To understand geographic genetic variation of the species and relationships among populations sequenced a portion of mitochondrial COI gene, which corresponds to “DNA Barcode” region (658 bp) from 81 individuals collected over eight localities in Korea. The sequence data were used to investigate genetic diversity within populations and species, geographic variation within species, phylogeographic relationship among populations, and phylogenetic relationship among haplotypes. Summarized, overall moderate to low genetic diversity within populations and species was characteristic, concordant with the high potential to disperse of O. cornifrons in Korea. Although two populations were genetically subdivided from the remaining localities, no clear regional bias was observed. Overall, high rate of gene flow among localities and low FST was characteristic considering other relevant studies that investigated population genetic structure of other insects occurring in Korean peninsula.

To understand geographic genetic variation of the species and relationships among populations of the bumble bee, Bombus ardens, we sequenced a portion of mitochondrial COI gene, which corresponds to "DNA Barcode" region (658 bp) from 101 B. ardens individuals collected over 11 localities in Korea. The sequence data were used to investigate genetic diversity within populations and species, geographic variation within species, phylogeographic relationship among populations, and phylogenetic relationship among haplotypes. A total of nine haplortpes were found, but they very close to each other (a maximum sequence divergence of 0. 304%). Summarized, overall moderate to low genetic diversity within populations and species was characteristic, concordant with the high potential to disperse of B. ardens in Korea. There was no clear regional subdivision was observed and relatively high rate of gene flow among localities and low FST was characteristic.

본 연구는 젖소(Holstein종)의 mtDNA D-loop 영역 염기변이 다형성과 경제형질간의 관련성을 분석하기 위하여 수행하였다. 젖소의 mtDNA D-loop 영역에서 단일염기의 치환에 의해 총 35개의 polymorphic site가 확인되었다. 그중 주요 Polymorphic site의 염기변이 빈도는 106, 169, 16057, 16231 및 16255 지역에서 높은 빈도의 염기치환이 검출되었다. 169 지역은 A가 G로 염기치환이 일어났고 그 빈도는 0.555로 높게 나타났으며 염기치환에 의한 산유량 효과는 1259.6 kg(P<0.05)로 나타났다. 유지방과의 관계를 분석한 결과 16118 지역은 A가 G로 치환되는 빈도가 0.02로 검출되었으며, 16135 지역은 T가 C로 치환되는 빈도가 0.02로 검출되었고, 16302 지역은 G가 A로 치환되는 빈도가 0.04로 검출되었다. 이 지역들이 유지방에 미치는 변이 효과는 - 156, - 118.15 및 - 67.77 kg(P<0.1)으로 나타났다. 본 연구에서 검출한 젖소 mtDNA내 D-loop 영역의 염기서열 변이 빈도와 유량, 유지방과의 연관성 분석 결과 등은 젖소집단의 유전적 변이성 추정과 좀 더 다양한 경제형질과의 관련성 분석으로 다양한 유전적 지표인자 발굴에 도움이 될 것이며 이를 통해 분자유전학적 기법을 이용한 젖소의 육종전략을 확립하는데 기초 자료가 되는 것은 물론 분자육종학적인 연구에 기초 자료로서 유용하게 활용할 수 있을 것으로 기대된다.

한국에 분포하고 있는 개구리 속에 대한 염기서열을 결정하고 상호 비교하여 종간 유전적 변이 정도를 밝히고자 한국산 개구리 속 6종과 일본산 개구리 속 1종에 대한 미토콘드리아 165 rDNA를 분석하였으며, Gene-bank에 수록된 일본산 산개구리류 3종도 함께 비교 분석하여 총 437 bp의 염기서열을 결정하였다. 산개구리류 7종에 대한 similarity는 91.3∼97.3％이며, 참개구리류는 96.1∼97.3％로 나타났다. 또한, 참개구리류와 옴

국내 늦반딧불이(Pyrocoelia rufa)이 미토콘드리아 DNA중 COI 유전자 일부(403 bp)의 염기서열을 결정, 집단내 유전적 다양도, 지역적 변이, 계통유전적 관련에 대한 분석을 실시하였다. 남해, 부산, 무주, 용인 등 우리나라 4개 지역으로부터 채집된 총 26개체로부터 7개의 mtDNA haplotype을 얻었으며 이들의 변이는 0.2~1.2%이었다. 도심인 부산에서 채집한 늦반딧불이는 다른 산림 및 농업의 채집지역과 달리 하나의 haplotype으로 고정되어 있어 도시화에 따른 집단의 병목현상과 서식처 파편화가 심각했음을 보여주었다. 그러나 우리나라 최대의 반딧불이 서식처이자 보호지역으로 지정된 무주로부터 4개의 haplotype을 얻었으며 이들의 최대염기 치환율은 1.0%로 가장 높은 집단내 유전적 다양도를 나타내었다. 근해의 섬인 남해의 늦반딧불이는 상대적으로 낮은 haplotype 다양도(H=0.25)와 계통유전적으로 이질적인 haplotype(PR7)의 존재로 요약되었는데 이는 비교적 가까운 과거의 한반도 생물지리 역사 및 유전자 이동에 의해 나타난 현상으로 설명하였다. 계층적 유전분석 결과 무주-용인과 부산-남해 그룹의 형성은 역사적으로 두 그룹사이에 유전자 이동에 반한 장벽의 존재 가능성을 제시한다고 설명하였다.

It is important for radish to have late flowering characteristics especially in the case of spring and winter cultivars. To understand late flowering characteristics of radish at the molecular level the flowering time genes of two radish lines (NH-JS1 and NH-JS2) with different flowering time were compared by re-sequencing their genomes. There were a total of 872,587 SNPs and 194,637 INDELs between the two lines. The SNP density of each chromosome was relatively uniform throughout, but the region with low SNP density was found at the end of R3 and the middle of R9. To compare the flowering time genes of the two lines, we first looked for the flowering time genes in radish using Arabidopsis thaliana flowering time genes. As a result, homologs of radish were found for most flowering time genes, but FRIGIDA was not found. Among 224 radish flowering time gene-homologs found, 97 genes showed more than one sequence difference (SNP or INDEL) between the two lines, and 127 genes had no difference. In particular, no sequence differences were found in FT, CO, and FLC, core flowering time control genes. Rs350520 (FVE), Rs193800 (CURLY LEAF) and Rs255320 (ATX1) with more than 100 sequence variations were expected to have a significant effect on flowering time difference between the two lines. These results will be of great help in understanding the flowering timing difference between the two lines at the molecular level.