Eumenis autonoe belonging to a lepidopteran family Nymphalidae (superfamily Papilionoidea) is an endangered species in Korea. Historically, the species was distributed in Europe and Asian region including a wide region in Korean peninsula. However, in Korean peninsula, the species is found only in two isolate dregions: South in a remote island Jeju, where altitude is higher than1, 400 meter on Halla Mt. and North in far northern Korean peninsula around Mt. Bekdu. In this study, we report the complete mitochondrial genome of the endangered E. autonoe collected from Mt. Halla. The 15,489-bp long E. autonoe genome has a typical gene content found in animal mitochondrial genomes and contains the gene arrangement identical to all other sequenced lepidopteran insects, which differs from the most common type found in insects, as the result of the movement of tRNAMet to a position 5’-upstream of tRNAIle. As seen in many other lepidopteran insects, no typical ATN codon for COI gene is available. Thus, we tentatively designated the CGA (arginine) found at the beginning of the COI gene, as has been suggested for lepidopteran COI starter. The intergenic spacer sequence located between tRNASer (UCN) and ND1 of E. autonoe mitogenome also contains the ATACTAA motif which is conserved in all sequenced lepidopteran species. The 678-bp long A+T-rich region, which is longest in sequenced lepidopteran insects contains ten identical tandem repeats composed of 27 bp plus one 13-bp long identical incomplete final repeat. Such repeat sequence is rare in the lepidopteran mitogenomes known so far. The E. autonoe A+T-rich region also contains a poly-T stretch located at the end of the region as 19 bp and also contains the downstream conserved motif ATAGA that were previously suggested to serve as a structural signal for minor-strand mtDNA replication. Phylogenetic analysis using the concatenated 13 amino acid sequences of PCGs among available six lepidopteran superfamilies (Tortricoidea, Pyraloidea, Papilionoidea, Bombycoidea, Geometroidea, and Noctuoidea) rooted with three dipteran species with BI and ML analyses supported the following topology: ((((Bombycoidea + Geometroidea +Noctuoidea) + Papilionoidea) + Pyraloidea) + Tortricoidea). Within Papilionoidea, a closer relationship between Lycaenidae and Pieridae, excluding Nymphalidae was observed. Further fruitful information will be available after more analysis is done.

The complete nucleotide sequences of the mitochondrial genome (mitogenome) from the white-spotted flower chafer, Protaetia brevitarsis (Coleoptera: Cetoniidae) was determined. The 20,319-bp long circular genome is the longest among the completely sequenced arthropods. This extraordinary length of the genome stemmed from 5,654-bp long A+T-rich region composed of twenty-eight 117-bp tandem repeats, seven 82-bp tandem repeats, and each two 19-bp and 38-bp tandem repeats. The P. brevitarsis contains a typical gene complement, order, and arrangement identical to most common type found in insects. The P. brevitarsis COI gene does not have typical ATN codon. Thus, we also designated it as AAC (asparagine), which is found in the start context of all sequenced Polyphaga within Coleoptera. All tRNAs showed stable canonical clover-leaf structure of other mt tRNAs, except for tRNASer (AGN), DHU arm of which could not form stable stem-loop structure. The 5bp-long motif sequence (TAGTA) that has been suggested to be the possible binding site for the transcription termination peptide for the major-strand also was found betweent RNASer (UCN) and ND1, as have been detected in all sequenced coleopteran insects.

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 15,389-bp long complete mitogenome of the endangered red-spotted apollo butterfly, Parnassius bremeri (Lepidoptera: Papilionidae) was determined. This genome has a gene arrangement identical to those of all other sequenced lepidopteran insects, which have the gene order of tRNAMet, tRNAIle, and tRNAGln at the beginning. Due to the uncertainty the start codon for COI gene in insect has been discussed extensively. We propose the CGA sequence as the start codon for COI gene in lepidopteran insects, based on complete mitogenome sequences of lepidopteran insects including our P. bremerii and additional sequences of the COI start region from a diverse taxonomic range of lepidopteran species (a total of 51 species belonging to 15 families). As has been suggested in other sequenced lepidopteran insects the 18 bp-long poly-T stretch and the downstream conserved motif ATAGA that were previously suggested to serve as a structural signal for minor-strand mtDNA replication also was found at the 3’-end region of the P. bremerii A+T-rich region. In an extensive search to find out tRNA-like structure in the A+T-rich region, each one tRNATrp-like sequence and tRNALeu (UUR)-like sequence were found in the P. bremeri A+T-rich region, and most of other sequenced lepidopteran insects were shown to have tRNA-like structure within the A+T-rich region, thereby indicating that such feature is frequent in the lepidopteran A+T-rich region. Phylogenetic analysis using the concatenated 13 amino acid sequences and nucleotide sequences of PCGs of the four macrolepidopteran suferfamilies together with Tortricoidea and Pyraloidea well recovered a monophyly of Papilionoidea and a monophyly of Bombycoidea. However, Geometroidea and Noctuoidea were unexpectedly clustered as one group and placed this group to the sister group to Bombycoidea, instead of Papilionoidea in most analyses.

Invertebrate mitochondrial genome contains 13 protein-coding genes and major start codons for them are ATA (Met) and ATG (Met). However, alternative start codons such as ATT (Ile), ATC (Ile), TTG (Leu), and GTG (Val) also have been suggested from a diverse organism. Approximately 120 complete mitochondrial genome reported showed that the start codon for COI gene evidences an array of diverse designation of COI start codon such as typical ATN, tetranucleotide TTAG and ATAA, newly proposed AAT and AAC and so on. In the case of Lepidoptera, many completely sequenced species showed no typical start codon at the start context of COI and even within the neighboring tRNATyr. In order to clarify, we newly sequenced the beginning context of COI gene, encompassing the neighboring tRNATyr and start region of COI gene from 39 species belonging to eight lepidopteran families. We found the newly sequenced 39 species and 14 available complete lepidopteran mitochondrial genomes all possessed CGA (arginine), which is the first non-overlapping in-frame codon in COI gene. Furthermore, this CGA is highly well aligned in terms of both nucleotide and amin o acid sequences with neighboring region. Thus, the CGA (arginine) may be synapomorphic character for Lepidoptera, functionally constrained. We, therefore, propose the CGA sequence as the start codon for COI gene in lepidopteran insects.

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.

To understand geographic genetic variation of the species and relationships among populations of the bumble bee, Bombus ardens, that is utilized as green house pollinator we expanded our investigation by sequencing somewhat longer mitochondrial DNA (mtDNA) fragment, covering some uninvestigated regions within the species distribution, and analyzing the sequence data in terms of population genetic structure. For the purpose of study, a portion of mitochondrial COI gene, corresponding to "DNA Barcode" region (658 bp) was sequenced from 160 individuals of B. ardens collected over 15 localities in Korea. The sequence data revealed overall relatively low genetic diversity within species, with a maximum sequence divergence of 0.3%. Geographically, one haplotype (BARBA01) was found in all localities surveyed, with the frequency of 91% (145 among 160 individuals), whereas other haplotypes were found in a locality mostly as a single individual, suggesting that haplotype distribution can be summarized as coexistence of widespread, one dominant haplotype and regionally restricted, other haplotypes. Overall, very high rate of per generation female migration (Nm = 4.6 ~ infinite) and very low level of geographic substitution (FST = 0 ~ 0.099) among localities were characteristic. Although some populations were genetically subdivided from the remaining localities in the hierarchical analysis, there was regional polarity on this subdivision. Taken together with gene flow estimates, the nature of genetic divergence of the bumble bee populations is characterized as one that possessing low genetic diversity, high gene flow, and wide spread of one dominant haplotype, consistent with the previous finding. To have further detailed information of this valuable genetic resource, further longer and variable molecular portion is under investigating.

We determined the complete mitochondrial genome (mitogenome) of the jewel beetle, Chrysochroa fulgidissima (Coleoptera: Buprestidae) from two overlapping fragments and subsequent sub fragments. The 15,592-bp long C. fulgidissima mitogenome contains gene arrangement and content identical to the most common arrangement found in insects. Most individual C. fulgidissima mitochondrial (mt) genes were well within the range found in the respective genes of other insects. The 875-bp A+T-rich region is shortest among the coleopteran mitogenomes sequenced in their entirety. The region is interesting in that it contains several stem-and-loop structures and tRNA-like structure found in the A+T-rich regions of other insect mitogenomes. As seen in other insect motogenomes the start codon of C. fulgidissima COI gene also is unusual because no typical start codon is available. Three of the 13 protein-coding genes have incomplete termination codon T or TA. All tRNA formed stable stem-and-loop structure, except for tRNASer(AGN), the DHU arm of which formed a simple loop as seen in many other metazoan mt tRNASer(AGN).

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.

The turnip aphid is a worldwide pest, damaging mainly to crucifers. In order to understand the life parameters of Lipaphis erysimi for the eventual goal of control, the developmental periods, survival rates, lifespan, and fecundity of the species were investigated under five temperature regimes (15℃ - 35℃). Furthermore, the efficacy of several environment-friendly agricultural materials (EFAMs) that are on the market was subjected to test in order to obtain further accurate information. The developmental period of the turnip aphid nymph was longest at 15ºC as 16.9 days, shortened as temperature goes up to 25ºC (5.4 days), and then somewhat increased at 30ºC (5.9 days), suggesting that the most efficient temperature for nymphal development could be around 25ºC. Mortality of the nymphal turnip aphid was obvious at 35ºC, whereas it was minimal at other temperature schemes. The longevity of adults shortened as temperature goes up to 30ºC. In particular, the maximum lifespan for adults continued for 55 days at 15ºC, but shortened to 21 days at 30ºC. The total fecundity per day was 35.7 at 15ºC, 81 at 20ºC, 64.2 at 25ºC, and 6.6 individuals at 30ºC, showing the highest fecundity at 20ºC. After the turnip aphids were successfully stabilized in indoor environment the insecticidal activity was tested and mortality was determined 12, 24, 36, and 48 hrs after EFAMs are treated. Several on-the-market EFAMs showed more than 90% of insecticidal activity within 24 hrs or 48 hrs, but a few showed less than 90% activity, signifying importance of selection of proper EFAMs.

In the process of bone remodeling, mineral phase of bone is dissolved by osteoclasts, resulting in elevation of calcium concentration in micro-environment. This study was performed to explore the effect of high extracellular calcium (Cα²+e) on mineralized nodule formation and on the expression of progressive ankylosis (Ank), plasma cell membrane glycoprotein-1 (PC-1) and osteopontin by primary cultured mouse calvarial cells. Osteoblastic differentiation and mineralized nodule formation was induced by culture of mouse calvarial cells in osteoblast differentiation medium containing ascorbic acid and β-glycerophosphate. Although Ank, PC-1 and osteopontin are well known inhibitors of mineralization, expression of these genes were induced at the later stage of osteoblast differentiation during when expression of osteocalcin, a late marker gene of osteoblast differentiation, was induced and mineralization was actively progressing. High Cα²+e(10 mM) treatment highly enhanced mRNA expression of Ank, PC-1 and osteopontin in the late stage of osteoblast differentiation but not in the early stage. Inhibition of p44/42 MAPK activation but not that of protein kinase C suppressed high Cα²+e- induced expression of Ank, PC-1 and osteopontin. When high Cα²+e (5 mM or 10 mM) was present in culture medium during when mineral deposition was actively progressing, matrix calcifiation was significantly increased by high Cα²+e. This stimulatory effect was abolished by pyrophosphate (5 mM) or levamisole (0.1-0.5 mM), an alkaline phosphatase inhibitor. In addition, probenecid (2mM), an inhibitor of Ank, suppressed matrix calcification in both control and high Cα²+e- treated group, suggesting the possible role of Ank in matrix calcification by osteoblasts. Taken together, these results showed that high Cα²+e stimulates expression of Ank, PC-1 and osteopontin as well as matrix calcification in late differentiation stage of osteoblasts and that p44/42 MAPK activation is involved in high Cα²+e- induced expression of Ank, PC-1 and osteopontin.

Dlx3 is a homeodomain protein and is known to play a role in development and differentiation of many tissues. Deletion of four base pairs in DLX3 (NT3198) is causally related to tricho-dento-osseous (TDO) syndrome (OMIM #190320), a genetic disorder manifested by taurodontism, hair abnormalities, and increased bone density in the cranium. The molecular mechanisms that explain the phenotypic characteristics of TDO syndrome have not been clearly determined. In this study, we examined phenotypic characteristics of wild type DLX3(wtDlx3) and 4-BP DEL DLX3 (TDO mtDlx3) in C2C12 cells. To investigate how wtDlx3 and TDO mtDlx3 differentially regulate osteoblastic differentiation, reporter assays were performed by using luciferase reporters containing the promoters of alkaline phosphatase, bone sialoprotein or osteocalcin. Both wtDlx3 and TDO mtDlx3 enhanced significantly all the reporter activities but the effect of mtDlx3 was much weaker than that of wtDlx3. In spite of these differences in reporter activity, electrophoretic mobility shift assay showed that both wtDlx3 and TDO mtDlx3 formed similar amounts of DNA binding complexes with Dlx3 binding consensus sequence or with ALP promoter oligonucleotide bearing the Dlx3 binding core sequence. TDO mtDlx3 exhibits a longer half-life than wtDlx3 and it corresponds to PESTfind analysis result showing that potential PEST sequence was missed in carboxy terminal of TDO mtDlx3. In addition, co-immunoprecipitation demonstrated that TDO mtDlx3 binds to Msx2 more strongly than wtDlx3. Taken together, though TDO mtDlx3 acted as a weaker transcriptional activator than wtDlx3 in osteoblastic cells, there is possibility that during in vivo osteoblast differentiation TDO mtDlx3 may antagonize transcriptional repressor activity of Msx2 more effectively and for longer period than wtDlx3, resulting in enhancement of osteoblast differentiation.

Dlx3 is a homeodomain protein and is known to playa role in development and differentiation of many tissues. Deletion of four base pairs in DLX3 (NT3198) is causally related to tricho-dento-osseous (TDO) syndrome (OMIM # 190320), a genetic disorder manifested by taurodontism, hair abnormalities, and increased bone density in the cranium. Although the observed defects of TDO syndrome involves bone, little is known about the role of Dlx3 in bone remodeling process. In this study, we examined the effect of wild type DLX3 (wtDlx3) expression on osteoclast differentiation and compared it with that of 4-BP DEL DLX3 (TDO mtDlx3). To examine whether Dlx3 is expressed during RANKL-induced osteoclast differentiation, RAW264.7 cells were cultured in the presence of receptor activator of nuclear factor-B ligand (RANKL). Dlx3 protein level increased slightly after RANKL treatment for 1 day and peaked when the fusion of prefusion osteoclasts actively progressed. When wtDlx3 and TDO mtDlx3 were overexpressed in RAW264.7 cells, they enhanced RANKL-induced osteoclastogenesis and the expression of osteoclast differentiation marker genes such as calcitonin receptor, vitronectin receptor and cathepsin K. Since osteoclast differentiation is critically regulated by the balance between RANKL and osteoprotegerin (OPG), we examined the effect of Dlx3 overexpression on expression of RANKL and OPG in C2C12 cells in the presence of bone morphogenetic protein 2. Overexpression of wtDlx3 enhanced RANKL mRNA expression while slightly suppressed OPG expression. However, TDO mtDlx3 did not exert significant effects. This result suggests that inability of TDO mtDlx3 to regulate expression of RANKL and OPG may contribute to increased bone density in TDO syndrome patients. Taken together, it is suggested that Dlx3 playa role as a positive regulator of osteoclast differentiation via up-regulation of osteoclast differentiation-associated genes in osteoclasts, as well as via increasing the ratio of RANKL to OPG in osteoblastic cells.

Tooth loss in elderly is mainly caused by alveolar bone loss via severe periodontitis. Although the severity of periodontitis is known to be affected by age, the aging process or the genetic changes during the aging of periodontal tissue cells are not well characterized. In this study, we investigated the effect of in vitro aging on the change of gene expression pattern in periodontal fibroblasts. Gingival fibroblasts (GF) and periodontal ligament fibroblasts (PDL) were obtained from two young patients and replicative senescence was induced by sequential subcultivation. When more than 90% cells were positively stained with senescence-associated β-galactosidase, those cells were regarded as aged cells. In aged GF and PDL, the level of phosphorylated retinoblastoma (RB) and p16INK4A protein was significantly decreased and increased, respectively. However, the protein level of p53 and p21, well known senescence-inducing genes, did not increase in aged GF and PDL. Although P27Kip1 and p15INK4B, another cyclin-dependent kinase inhibitors, were reported to be involved in replicative senescence of human cells, they were decreased in aged GF and PDL. Because senescent cells showed flattened and enlarged cell shape and are known to have increased focal adhesion, we examined the protein level of several integrins. Aged GF and PDL showed increased protein level of integrin α2, αu, and β1. When the gene expression profiles of actively proliferating young cells and aged cells were compared by cDNA microarray of 3,063 genes and were confirmed by reverse transcription-polymerase chain reaction, 7 genes and 15 genes were significantly and commonly increased and decreased, respectively, in aged GF and PDL. Among them, included are the genes that were known to be involved in the regulation of cell cycle, gene transcription, or integrin signaling. The change of gene expression pattern in GF and PDL was minimally similar to that of oral keratinocyte. These results suggest that p16INK4A/RB might be involved in replicative senescence of periodontal fibroblasts and the change of gene expression profile during aging process is cell type specific.