The spotted-wing drosophila (SWD), Drosophila suzukii (Diptera: Drosophilidae), is an economically damaging pest that feeds on most thin-skinned fruits. In this study, we sequenced portions of the mitochondrial (mt) COI and ND4 genes from a total of 195 individuals collected mainly from Korea. A total of 139 haplotypes were obtained from the concatenated COI and ND4 sequences. A dataset combining GenBank sequences with our own data identified a total of 94 worldwide COI haplotypes with a maximum sequence divergence of 5.433% (32 bp). A rough estimate of genetic diversity in each country showed higher diversity in ancestral distributional ranges, but the invasion over Asian countries seems to have been substantial because haplotype diversity was only 2.35-3.97-fold lower in the USA, Canada, and Italy than that in the populations ancestral ranges.

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 lycanid butterfly, Shijimiaeoides divina (Lepidoptera: Lycaenidae) is listed as the second-degree endangered wild animal in Korea from 2012. The 15,259-bp long complete mitochondrial genome of the species consisted of a typical set of genes (13 protein-coding genes, two rRNA genes and 22 tRNA genes) and one major non-coding A+T-rich region, with the typical arrangement found in majority of Lepidoptera. The 15,259-bp long S. divina mitogenome is well within the range found in Lycaenidae and has typical sets of 37 genes and a major non-coding A+T-rich region as 379 bp. As other lycanid butterflies S. divina COI also started with CGA. The gene arrangement of S. divina is identical to that of the Ditrysia in Lepidoptera that has the order trnM-trnI-trnQ (underline for inverted gene) between the A+T-rich region and ND2. Comparison of the skewness between the PCGs encoded in major and minor strand indicates a substantial difference between them in GC skewness (0.261 ~ 0.340 in minor strand vs. -0.081 ~ -0.115 in major strand). The 151-bp intergenic spacer sequence of the S. divina mitogenome is spread over 16 regions ranging in size from 1-53 bp. The longest one (53 bp) located between trnQ and ND2 shows substantially high sequence homology to neighboring ND2 may indicating the origination of the region by a partial duplication of the ND2 gene. One of the unusual features of the S. divina mitogenome is the presence of a trnK-like sequence that is encoded at the major strand of the genome in the A+T-rich region.

We have determined the complete mitochondrial genome of the yellow-spotted long horned beetle, Psacothea hilaris (Coleoptera: Cerambycidae), an endangered insect species in Korea. The 15,856-bp long P. hilaris mitogenome harbors gene content typical of the animal mitogenome and a gene arrangement identical to the most common type found in insect mitogenomes. As with all other sequenced coleopteran species, the 5-bp long TAGTA motif was also detected in the intergenic space sequence located between tRNASer (UCN) and ND1 of P. hilaris. The 1,190-bp long non-coding A+T-rich region harbors an unusual series of seven identical repeat sequences of 57-bp in length and several stretches of sequences with the potential to form stem-and-loop structures. Furthermore, it contains one tRNAArg-like sequence and one tRNALys-likes equence. Phylogenetic analysis among available coleopteran mitogenomes using the concatenated amino acid sequences of PCGs appear to support the sister group relationship of the suborder Polyphaga to all remaining suborders, including Adephaga, Myxophaga, and Archostemata. Among the two available infraorders in Polyphaga, a monophyletic Cucujiformia was confirmed, with the placement of Cleroidea as the basal lineage for Cucujiformia. On the other hand, the infraorder Elateriformia was not identified as monophyletic, thereby indicating that Scirtoidea and Buprestoidea are the basal lineages for Cucujiformia and the remaining Elateriformia.

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 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.