In order to understand evolutionary characteristics of gene rearrangement in Lepidoptera, we collected all available complete mitogenome (mitogenome) sequences registered in GenBank (274 mitogenomes from 44 families in 23 superfamilies as of August 6, 2015). It turned up six rearrangements that differ from the arrangement of ancestral insects, including that of the gelechioid Mesophleps albilinella that we sequenced in this study. The M. albilinella mitogenome has a unique gene arrangement among the Gelechioidea: ARNESF (the underline signifies an inverted gene) at the ND3 and ND5 junction, as opposed to the ARNSEF that is found in ancestral insects. Most of the rearrangements can be explained by the tandem duplication-random loss model, but inversion, which requires recombination, is also found in two cases, including M. albilinella. Excluding the MIQ rearrangement at the A+T-rich region and ND2 junction, which is found in nearly all Ditrysia, most of the remaining rearrangements found in Lepidoptera appear to be independently derived in that they are automorphic at several taxonomic scales. Current mitogenomic data are limited, particularly for congeneric data. Thus, future research focused on congenerics could clarify evolutionary independency at the generic level also.

The Gelechioidea is the second most species-rich group of Lepidoptera, but only limited number of mitochondrial genome (mitogenome) sequences is available. Thus, we sequenced the complete mitogenome of a gelechioid Hieromantis kurokoi (Lepidoptera: Stathmopodidae) to use the data for future study for the higher phylogeny of Ditrysia in Lepidoptera. The arrangement of the genome was identical to typical one found in Ditrysia (trnM-trnI-trnQ) (underline for inverted gene). The COI began with CGA, which has been designated as the start codon for majority of lepidopteran species, whereas other protein-coding genes (PCGs) began with the typical ATN codon. The 360-bp long A+T-rich region harbored the conserved sequence blocks Phylogenetic analysis using the 13 PCGs both by Bayesian inference (BI) and Maximum-likelihood (ML) methods indicated that H. kurokoi belonging to the family Stathmopodidae grouped together with within-familial species Atrijuglans hetaohei with the highest nodal support (BI, 1.0; ML, 100%).

The spotted-wing drosophila (SWD), Drosophila suzukii (Diptera: Drosophilidae), originally distributed across a few Asian countries, including South Korea, has invaded North America and Europe, but is absent from Australia. In order to export the South Korean grape cultivar Campbell Early to Australia, its potential to serve as oviposition and development medium for SWD must first be determined. In this study, we determined the oviposition and development potential of SWD on Campbell Early, after elucidating the SWD life cycle and establishing an artificial diet-based mass-culturing system. An investigation of the life cycle under five temperature regimes (16, 19, 22, 25, and 28°C) showed that the durations of the egg, larval, and adult stages were shorten when temperature was increased from 16, 19, 22, 25, and 28°C, but pupal duration was shortest at 25°C and extended again at 28°C. A test of oviposition and development potential of SWD on Campbell Early grape clusters showed oviposition of 30.8 ± 6.8 eggs per cluster of injured grapes and 157.7 ± 16.2 eggs on a culture dish of artificial diet. However, in a similar experiment using uninjured grape clusters, only a single egg was deposited on the grape skin, which soon dried. In light of these results, newly harvested grapes left at vineyards during daily harvests are unlikely to serve as an oviposition and development medium for SWD, as long as the grapes remain uninjured.

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.

We newly sequenced mitogenomes of five skippers belonging to Lepidoptera to obtain further insight into characteristics of butterfly mitogenomes and performed phylogenetic reconstruction using all available gene sequences (PCGs, rRNAs, and tRNAs) from 85 species in 19 families in eight superfamilies. The general genomic features found in the butterflies also were found in the five skippers: a high A/T composition (79.3% - 80.9%), dominant usage of TAA stop codon, similar skewness pattern in various levels, consistently long intergenic spacer sequence between tRNAGln - ND2 (64-87 bp), the ATACTAA motif betweent RNASer(UCN) and ND1, and characteristic features of the A+T-rich region (the motif ATAGA, varying length of poly-T stretch, and poly-A stretch). The start codon for COI was CGA in four skippers as typical, but Lobocla bifasciatus evidently possessed canonical ATG as start codon. Phylogenetic analyses mainly yielded the consensus superfamilial relationships ((((((Bombycoidea + Noctuoidea + Geometroidea) + Pyraloidea) + Papilionoidea) + Tortricoidea) + Yponomeutoidea) + Hepialoidea) with a high support for most nodes, confirming the validity of Macroheterocera and its sister relationship to Pyraloidea. Within Rhopalocera the familial relationships (Papilionidae + (Hesperiidae + (Pieridae + ((Lycaenidae + Riodinidae) + Nymphalidae))) were strongly supported, confirming invalidity of the superfamily Hesperioidea. On the other hand, superfamilial relationships among Noctuoidea, Geometroidea, and Bombycoidea and the familial relationships among Saturniidae, Sphingidae, and Bombycidae were dubious, requiring further representative taxon sampling.

Glyphodes quadrimaculalis (Lepidoptera: Crambidae) feed on a root tuber of Cynanchum wilfordii (Asclepiadaceae) that is one of the most famous traditional medicines in Korea. The genus Glyphodes includes ~130 species distributed worldwide, so the complete mitochondrial genome (mitogenome) would be helpful for bio-identification, biogeographic studies, and multigene-based phylogeny. The 15,255-bp long G. quadrimaculalis genome comprises 37 typical genes and one large non-coding region, with the typical arrangement found in Lepidoptera. Of the 13 protein coding genes (PCGs), 12 begin with typical start codons found in insect mitochondrial PCGs, but the COI gene starts with atypical CGA. One of the noteworthy features of the genome includes the presence of a 51-bp long non-coding space sequence located between tRNAGln and ND2 that reveals high sequence homology (71.4%) to the neighboring ND2 gene, indicating the origin of the region by partial duplication of the ND2 gene.

Hipparchia autonoe belongs to the family Nymphalidae (Lepidoptera) and is designated as an endangered insect and national monument in Korea. It only inhabits a very restricted area on Mt. Halla but is widely distributed in several Asian countries including Mongolia. A previous study conducted to understand the genetic relationship between Mt. Halla and Mongolian H. autonoe for conservation purposes suffered from a limited number of samples. Therefore, we sequenced the DNA barcode region of an additional 36 H. autonoe individuals, combined them with previous data from 19 individuals, and performed phylogenetic and population genetic analyses. Furthermore, the internal transcribed spacer 2 (ITS2) region was also sequenced from the 36 samples as a nuclear DNA marker. The existence of independent haplotypes, sequence types, and significant FST estimates (p < 0.05) between Mt. Halla and Mongolian populations indicated hampered gene flow between the populations. Nevertheless, an absence of a reciprocal monophyletic group in Mt. Halla and Mongolian populations by cytochrome oxidase subunit I gene- and ITS2-based phylogeny suggests that the genetic isolation of the Mt. Halla population from the Mongolian populations seemed not large enough to consider them independent genetic entities.

Lepidoptera is one of the largest insect orders, but the phylogenetic relationships within this order, have yet to be completely described. One of the unresolved relationships includes the monophyly of Papilionoidea in relationship with the monotypic superfamily Hesperioidea. We newly sequenced five hesperid mitochondrial genomes (mitogenomes), representing four subfamilies: Pyrginae (Daimio tethys and Lobocla bifasciatus), Coeliadinae (Choaspes benjaminii), and Hesperiinae (Potanthus flavus), and Heteropterinae (Carterocephalus silvicola). Along with these newly sequenced hesperid genomes phylogenetic analysis was conducted with all available lepidopteran mitogenomes including three reported species of Hesperiidae that consisted of ~70 species in ten lepidopteran superfamilies. The test for the effect of optimization schemes, such as exclusion and inclusion of third codon position of 13 PCGs, other genes (22 tRNAs and two rRNAs), and with and without partitions also was performed. Majority of datasets consistently placed the monophyletic Hesperiidae the sister to ((Pieridae + Lycaenidae) + Nymphalidae), placing another true butterfly family Papilionidae as the basal lineage of this group, presenting the relationships (Papilionidae + (Hesperiidae + ((Pieridae + Lycaenidae) + Nymphalidae))). Consistent to previous result, Pyraloidea was placed as the sister to ((Bombycoidea + Geometroidea) + Noctuoidea), placing the Macrolepidoptera as non-monophyletic group.

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.

In the present study, the 17,694-bp long complete mitochondrial genome (mitogenome) of the dwarf honey bee, Apis florea (Hymenoptera: Apidae), is described with an emphasis on the noteworthy triplicated tRNAser(AGN) region and an extraordinary long A+T-rich region with repeat regions. The gene arrangement of A. florea mitogenome is identical to that of A. mellifera, but has triplicated tRNASer(AGN), each of which contains the precedent 44 bp-long and following another 64 bp-long repeats plus one complete first repeat abutting to tRNAMet. A total of 1,610-bp long two repeat regions in 1,987 bp-long A+T-rich region is composed of nearly identical 141 ~ 219-bp long five tandem repeats and 50 ~ 52-bp long 12 tandem repeats that are encompassed by three non-repeat sequences. One of the common interpretations for such repeat sequence is slipped-strand mispairing and unequal crossing-over events during DNA replication.

We previously reported Pear Pest Forecasting Management System (PPFMS) for the Improvement of pass ratio of Korean exporting pears. It is consisted of regular field forecasting by pear farmers, meteorological data obtained by automatic weather station (AWS), an internet web page (http://pearpest.jnu.ac.kr/) as information collecting and providing ground, and information providing service. Currently, we are expanding this system to the area, Cheonan and Ansung, where pear orchards are organized into exportation-specific group. Further, the information obtained from field forecasting and AWS were up-loaded to under-constructing upgraded webpage (http://www.kpear.kr), with several pest/disease-related information. We hope this pest forecasting management system increases the pass ratio of Korean exporting pears throughout establishment of farmer-oriented forecasting, inspiring farmers’ effort for the prevention and forecasting of diseases and pests occurring at pear orchards.

The black-veined white, Aporia crataegi (Lepidoptera: Papilionoidea), is nearly extinct in South Korea, although substantial numbers of dried specimens are available. One of the common practices for such species is to launch re-introduction program after proper amount of genetic information are analyzed from donor and donee populations. In this study, we sequenced complete mitochondrial genome (mitogenome) of A. crataegi to design species-specific primers for subsequent population works and to further understand the mitogenome evolution in lepodiopteran Papilionoidea. The 15,140-bp long A. crataegi mitogenome that has typical sets of 37 genes is smallest among true butterfly species with overall slightly smaller size in genes and regions throughout the genome. Arrangement of the genome is identical to those of other lepidopteran mitogenomes, in which tRNA cluster located between the A+T-rich region and ND2 gene is translocated into tRNAMet, tRNAIle, and tRNAGln from ancestral arrangement, tRNAIle, and tRNAGln, tRNAMet. The A/T content of the genome at 81.3% is the highest in Pieridae, but lower than that of lycaenid species (81.7% ~ 82.7%) The high A/T content in the genome is also reflected in codon usage, accounting for 41.69% of A/T-composed codons (TTA, ATT, TTT, and ATA). Unlikely the diversified or modified usage of anticodon for tRNASer(AGN) the species of Pieridae including A. crataegi all unanimously have GCT that has been hypothesized as ancestral for Lepidoptera. A total of 111 bp of non-coding sequences are dispersed in 13 regions, ranging in size from 1–49 bp. Among them relatively longer ones (≥ 16 bp) all have relatively higher sequence identity to other regions of the genome, suggesting partial duplication of the sequences during A. crataegi evolution. As has been reported in some species of Lepidoptera, the A. crataegi A+T-region also has typically found conserved sequences (e.g., poly-T stretch, ATAGA motif, ATTTA element, microsatellite-like A/T sequence, and poly-A stretch) and one tRNA-like sequence, and this feature was commonly found in true butterfly species.

The phylogenetic relationships among the Nymphalidae (Lepidoptera: Papilionoidea) have been controversial in several perspective. The present study sequenced a total of ~ 3,500 bp from cytochrome oxidase subunit I (COI), 16S ribosomal RNA (16S rRNA), and elongation factor-1 alpha (EF-1α) in 80 nymphalid species belonging to seven subfamilies (Linmenitidinae, Heliconiinae, Nymphalinae, Apaturinae, Libytheinae, Satyrinae, and Danainae), along with those of six lycaenid species as outgroups. Phylogenetic analyses via Bayesian Inference (BI) and Maximum Likelihood (ML) algorithms concordantly supported the subfamilial relationships of (((((Linmenitidinae + Heliconiinae) + (Nymphalinae + Apaturinae)) + Libytheinae) + Satyrinae) + Danainae), with high nodal support for monophyletic subfamilies and tribes. This result is largely consistent with a previous study performed with a substantially large sequence information and morphological characters, except for the position of Libytheinae that has previously been placed as the sister to all reminder of Nymphalidae.

The Samia cynthia ricini (Lepidoptera: Saturniidae) is a commercial silk-producing insect belonging to an insect family Saturniidae in Bombycoidea. The species that has presumably been originated in India, is distributed in India, China, and Japan. Unlikely domestic silkworm the prime host plant for the species is a castor-oil plant (Ricinus communis in Euphorbiaceae). Recently, the eri-silkworm also is reared in Korea and is expected to be utilized for a diverse purpose. In this report, we present the complete mitochondrial genome of the species with the emphasis of a few major characteristics. The 15,384-bp long S. cynthia ricini (Lepidoptera: Saturniidae) mitochondrial genome was amplified into three long overlapping fragments (from COI ~ ND4, ND5 ~ lrRNA, and lrRNA ~ COI) and subsequent several short fragments using the long fragments as temperate. The primers for both long and short fragments were designed solely for lepidopteran genomes, without any species-specific primers. As a usual the genome is composed of 37 genes: 13 protein-coding genes (PCGs), two rRNA genes, and 22 tRNA genes, and one large non-coding region termed the A+T-rich region. Arrangement of the genome is identical to those of other lepidopteran mitochondrial genome, but this differs from the common arrangement found in a diverse insect order, by the movement of tRNAMet to a position 5’- up stream of tRNAIle. Unlikely previous report on the start codon for COI gene in Lepidoptera S. cynthia ricini COI gene starts with typical ATT codon located between tRNATyr and the beginning region of COI gene. The 22 tRNAs that are interspersed throughout the mitogenome ranged in length from 62 to 71 bp. All tRNAs but tRNASer(AGN) were shown to be folded into the expected cloverleaf secondary structures. More detailed structural and phylogenetic analyses among Bombycidae and Saturniidae in connection with other families in the Bombycoidea will be performed soon

The leaf beetle, Chrysolina aurichalcea (Coleoptera: Chysomelidae), is a pest damaging plants of Compositae. In order to understand the genetic diversity and geographic variation of the species we sequenced a portion of mitochondrial COI gene (658 bp) and complete nuclear internal transcribed spacer 2 (ITS2) collected from seven Korean localities. A total of 18 haplotypes (BARCA01 ~ BARCA18), with the maximum sequence divergence of 3.04% (20 bp) were obtained from COI gene sequence, whereas 17 sequence types (ITS2CA01 ~ ITS2CA17), with the maximum sequence divergence of 2.013% (9 bp) were obtained from ITS2, indicating substantially larger sequence divergence in mitochondrial gene sequence. Phylogenetically, the mitochondrial DNA has shown several haplotypes formed independent groups with substantially high node support (≥ 90%), whereas no such grouping was evidenced for ITS2, indicating different behaviors of the two molecules. Such difference may reflect a diverse dynamics of the species such as biogeographic history, mating behaviors, and also possibly different mode of inheritance of the two molecules, but requires further scrutinized examination of the dataset. In terms of population genetic perspective, overall no population subdivision was detected from both molecules, except for locality 7 (Eocheong islet) from mitochondrial DNA. As more scrutinized analysis is performed, further fruitful inference on the geographic contour of the species might be available.