Complete nucleotide sequence and organization of the mitogenome of the red-spotted apollo butterfly, Parnassius bremeri (Lepidoptera: Papilionidae) and comparison with other lepidopteran insects

The 15,389-bp long complete mitogenome of the endangered red-spotted apollo butterfly, Parnassius bremeri (Lepidoptera: Papilionidae) was determined in this study. The start codon for the COI gene in insects has been extensively discussed, and has long remained a matter of some controversy. Herein, we propose that the CGA (arginine) sequence functions as the start codon for the COI gene in lepidopteran insects, on the basis of complete mitogenome sequences of lepidopteran insects, including P. bremeri, as well as additional sequences of the COI start region from a diverse taxonomic range of lepidopteran species (a total of 53 species from 15 families). In our extensive search for a tRNA-like structure in the A+T-rich region, one tRNA^Trp-like sequence and one tRNA^Leu (UUR)-like sequence were detected in the P. bremeri A+T-rich region, and one or more tRNA-like structures were detected in the A+T-rich region of the majority of other sequenced lepidopteran insects, thereby indicating that such features occur frequently in the lepidopteran mitogenomes. Phylogenetic analysis using the concatenated 13 amino acid sequences and nucleotide sequences of PCGs of the four macrolepidopteran superfamilies together with the Tortricoidea and Pyraloidea resulted in the successful recovery of a monophyly of Papilionoidea and a monophyly of Bombycoidea. However, the Geometroidea were unexpectedly identified as a sister group of the Bombycoidea, rather than the Papilionoidea.     

Complete nucleotide sequence and organization of the mitogenome of the silk moth Caligula boisduvalii (Lepidoptera: Saturniidae) and comparison with other lepidopteran insects

The 15,360-bp long complete mitogenome of Caligula boisduvalii possesses a gene arrangement and content identical to other completely sequenced lepidopteran mitogenomes, but different from the common arrangement found in most insect order, as the result of the movement of tRNA(Met) to a position 5'-upstream of tRNA Ile. The 330-bp A+T-rich region is apparently capable of forming a stem-and-loop structure, which harbors the conserved flanking sequences at both ends. Dissimilar to what has been seen in other sequenced lepidopteran insects, the initiation codon for C. boisduvalii COI appears to be TTG, which is a rare, but apparently possible initiation codon. The ATP8, ATP6, ND4L, and ND6 genes, which neighbor another PCG at their 3' end, all harbored potential sequences for the formation of a hairpin structure. This is suggestive of the importance of such structures for the precise cleavage of the mRNA of mature PCGs. Phylogenetic analyses of available sequenced species of Bombycoidea, Pyraloidea, and Tortricidea supported the morphology-based current hypothesis that Bombycoidea and Pyraloidea are monophyletic (Obtectomera). As previously suggested, Bombycidae (Bombyx mori and B. mandarina) and Saturniidae (Antheraea pernyi and C. boisduvalii) formed a reciprocal monophyletic group.     

The complete mitochondrial genome of the butterfly <Emphasis Type="Italic">Apatura metis</Emphasis> (Lepidoptera: Nymphalidae)

As an important pest in the Slender Leaved Willow (Salix alba), Apatura metis is called Freyer’s purple emperor, and its mitochondrial genome is 15,236 bp long. The encoded genes for 22 tRNA genes, two ribosomal RNA (rrnL and rrnS) genes, and 13 protein-coding genes (PCGs), and a control region in the A. metis mitochondria are highly homologous to other lepidopteran species. The mitochondrial genome of A. metis is biased toward a high A + T content (A + T = 80.5%). All protein-coding genes, except for COI begins with the CGA codon as observed in other lepidopterans, start with a typical ATN initiation codon. All tRNAs show the classic clover-leaf structure, except that the dihydrouridine (DHU) arm of tRNA ^Ser(AGN) forms a simple loop. The A. metis A + T-rich region contains some conserved structures including a structure combining the motif ‘ATAGA’ and 19 bp poly (T) stretch, which is similar to those found in other lepidopteran mitogenomes. The phylogenetic analyses of lepidopterans based on mitogenomes sequences demonstrate that each of the six superfamilies is monophyletic, and the relationship among them is (((Noctuoidea + (Geometroidea + Bombycoidea)) + Pyraloidea) + Papilionoidea) + Tortricoidea. In Papilionoidea group, our conclusion argues that ((Lycaenidae + Pieridae) + Nymphalidae) + Papilionidae.     

Complete mitochondrial genome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae) and the phylogenetic relationship of Saturniidae species

Mitochondrial genome (mitogenome) can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. In this study, we present the complete mitogenome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae), a well-known silk-producing and ornamental insect with the largest wing surface area of all moths. The mitogenome of A. atlas is a circular molecule of 15,282 bp long, and its nucleotide composition shows heavily biased towards As and Ts, accounting for 79.30%. This genome comprises 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A + T-rich region. It is of note that this genome exhibits a slightly positive AT skew, which is different from the other known Saturniidae species. All PCGs are initiated by ATN codons, except for COI with CGA instead. Only six PCGs use a common stop codon of TAA or TAG, whereas the remaining seven use an incomplete termination codon T or TA. All tRNAs have the typical clover-leaf structure, with an exception for tRNA^Ser(AGN). The A. atlas A + T-rich region contains non-repetitive sequences, but harbors several features common to the Bombycoidea insects. The phylogenetic relationships based on Maximum Likelihood method provide a well-supported outline of Saturniidae, which is in accordance with the traditional morphological classification and recent molecular works.     

The complete mitochondrial genome of <Emphasis Type="Italic">Spilonota lechriaspis</Emphasis> Meyrick (Lepidoptera: Tortricidae)

We determined the nucleotide sequence of the mitochondrial genome (mtgenome) of Spilonota lechriaspis Meyrick (Lepidoptera: Tortricidae). The entire closed circular molecule is 15,368 bp and contains 37 genes with the typical gene complement and order for lepidopteran mtgenomes. All tRNAs except tRNA^Ser(AGN) can be folded into the typical cloverleaf secondary structures. The protein-coding genes (PCGs) have typical mitochondrial start codons, with the exception of COI, which uses the unusual CGA one as is found in all other Lepidoptera sequenced to date. In addition, six of 13 PCGs harbor the incomplete termination codons, a single T. The A + T-rich region contains some conserved structures that are similar to those found in other lepidopteran mtgenomes, including a structure combining the motif ‘ATAGA’, a 19-bp poly(T) stretch and three microsatellite (AT)_n elements which are part of larger 122+ bp macrorepeats. This is the first report of macrorepeats in a lepidopteran mtgenome.     

The mitochondrial genome of the butterfly Papilio xuthus (Lepidoptera: Papilionidae) and related phylogenetic analyses

The nearly complete mitochondrial genome of the butterfly Papilio xuthus (Lepidoptera: Papilionidae) was sequenced for its nucleotide sequence of 13,964 bp. The genome has a typical gene order identical to other lepidopteran species. All tRNAs showed same stable canonical clover-leaf structure as those of other insects, except for tRNA^Ser (AGN), in which the dihydrouracil arm (DHU arm) could not form stable stem–loop structure. Anomalous initiation codons have been observed for the cox1 gene, where the ATTACG hexa-nucleotide was believed to be involved in the initiation signaling. Twelve mitochondrial protein-coding gene sequence data were used to infer the phylogenetic relationships among the insect orders. Even though the number of insect orders represented by complete mitochondrial genomes is still limited, several well-established relationships are evident in the phylogenetic analysis of the complete sequences. Monophyly of the Homometabola was not supported in this paper. Phylogenetic analyses of the available species of Bombycoidea, Pyraloidea, Papilionoidea and Tortricidea bolstered the current morphology-based hypothesis that Bombycoidea, Pyraloidea and Papilionoidea are monophyletic (Obtectomera). Bombycoidea (Bombyx mandarina and Antheraea pernyi) and Papilionoidea (P. xuthus and Coreana raphaelis) formed a sister group.     

Complete mitochondrial genomes of five skippers (Lepidoptera: Hesperiidae) and phylogenetic reconstruction of Lepidoptera

We sequenced mitogenomes of five skippers (family Hesperiidae, Lepidoptera) to obtain further insight into the characteristics of butterfly mitogenomes and performed phylogenetic reconstruction using all available gene sequences (PCGs, rRNAs, and tRNAs) from 85 species (20 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 both strands, consistently length intergenic spacer sequence between tRNA^Gln and ND2 (64–87 bp), conserved ATACTAA motif between tRNA^{Ser (UCN)} and ND1, and characteristic features of the A + T-rich region (the ATAGA motif, 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. All species had the ancestral arrangement tRNA^Asn/tRNA^{Ser (AGN)}, instead of the rearrangement tRNA^{Ser (AGN)}/tRNA^Asn, found in another skipper species (Erynnis). Phylogenetic analyses using all available genes (PCGs, rRNAS, and tRNAs) yielded the consensus superfamilial relationships ((((((Bombycoidea + Noctuoidea + Geometroidea) + Pyraloidea) + Papilionoidea) + Tortricoidea) + Yponomeutoidea) + Hepialoidea), confirming the validity of Macroheterocera (Bombycoidea, Noctuoidea, and Geometroidea in this study) and its sister relationship to Pyraloidea. Within Rhopalocera (butterflies and skippers) the familial relationships (Papilionidae + (Hesperiidae + (Pieridae + ((Lycaenidae + Riodinidae) + Nymphalidae)))) were strongly supported in all analyses (0.98–1 by BI and 96–100 by ML methods), rendering invalid the superfamily status for Hesperioidea. On the other hand, current mitogenome-based phylogeny did not find consistent superfamilial relationships among Noctuoidea, Geometroidea, and Bombycoidea and the familial relationships within Bombycoidea between analyses, requiring further taxon sampling in future studies.     

The complete mitochondrial genome of Paracymoriza prodigalis (Leech, 1889) (Lepidoptera), with a preliminary phylogenetic analysis of Pyraloidea

The complete mitochondrial genome sequence is determined for Paracymoriza prodigalis (Leech, 1889). The 15,326 bp circular molecule possesses a gene organization and order identical to other sequenced Pyraloidea mitochondrial genomes. All tRNAs have the typical clover-leaf structure except for tRNA^Ser(AGN), which lacks the dihydrouridine (DHU) arm. The A+T-rich region of 343 bp includes the features common to the Lepidoptera, including the ‘ATAGA’ followed by an 19-bp poly-T stretch, but the tandem repeat sequences often appearing in available insects are not found. Phylogenetic relationships of eight subfamilies of 14 Pyraloidea species were constructed based on 13 PCGs of mitochondrial genomes using Bayesian inference (BI) and maximum likelihood (ML) methods. These phylogenies of the subfamilies within Pyraloidea accord well with morphological phylogenetic analysis except for the position of Schoenobiinae.     

武铠蛱蝶线粒体基因组全序列测定和分析

【目的】了解闪蛱蝶亚科属间及种间的分子系统进化关系。【方法】采用PCR步移法对武铠蛱蝶 Chitoria ulupi 线粒体基因组全序列进行了测定和分析。基于线粒体基因组13个蛋白质编码基因的核苷酸序列构建了38种鳞翅目昆虫的系统发育树。【结果】分析结果表明,武铠蛱蝶线粒体基因组全长15 279 bp,包括13个蛋白质编码基因、22个tRNA基因、2个rRNA基因和一段长度为391 bp的A+T富含区,基因排列顺序与其他已知近缘种昆虫相同。武铠蛱蝶线粒体基因组中存在很高的A+T含量（79.9％）。13个蛋白质编码基因中,COII以TTG作为起始密码子,COI以CGA作为起始密码子外,其余均为昆虫典型的起始密码子ATN。COII和ND4基因使用了不完全终止密码子T,其余基因均以典型的TAA为终止密码子。在所测得的22个tRNA基因中,除tRNA^Ser(AGN)缺少DHU臂外,其余tRNA均能形成典型的三叶草结构。与其他多数鳞翅目昆虫一样,武铠蛱蝶的A+T富含区中有一段由ATAGAA引导的保守的多聚T结构,长度为21 bp,并散布着一些长短不一的串联重复单元。系统发育树结果显示,总科级别的系统发育关系为：卷蛾总科＋（凤蝶总科＋（螟蛾总科＋(夜蛾总科＋蚕蛾总科＋尺蛾总科)））;在蛱蝶科物种中,武铠蛱蝶与猫蛱蝶Timelaea maculate 亲缘关系最近。【结论】基于分子标记构建的鳞翅目昆虫系统发育关系与传统的形态学分类结果基本一致。     

The complete mitochondrial genome of <Emphasis Type="Italic">Deracantha onos</Emphasis> (Orthoptera: Bradyporidae)

The complete mitochondrial genome 15,650 bp in size of the Deracantha onos has been determined. The gene content, base composition and codon usage of D. onos are coincident to typical hexapods mitochondrial genomes. Genes arrangement of D. onos is identical to Gryllotalpa orientalis, Ruspolia dubia and Anabrus simplex, in that the relative locations of tRNA^Lys and tRNA^Asp was different to that of Locusta migratoria. All tRNAs could be folded into the typical cloverleaf secondary structure, excluding tRNA^Ser(AGN) which forms another structure according to the Steinberg–Cedergren tertiary structure. Sequence analysis of the A + T-rich region with Dot-plot did not find any conspicuous repeat clusters. Two poly-thymine (poly-T) nucleotide stretches of 20 bp and 11 bp in size, which may involved in the recognition of replication origin, were found on the H-strand and L-strand in the A + T-rich region of the D. onos mitogenome, respectively. One open reading frame (ORF) 87 amino acids in size was found on the H-strand, but Protein Blast searches analysis indicated that it was a nonfunctional ORF.     

The complete mitochondrial genome of <Emphasis Type="Italic">Thrinchus schrenkii</Emphasis> (Orthoptera: Caelifera,Acridoidea, Pamphagidae)

Complete nucleotide sequences of mitochondrial genome (mitogenome) of Thrinchus schrenkii (Orthoptera: Acridoidea: Pamphagidae) were determined. It is 15672 bp in length and contains 71.2% A + T. All T. schrenkii protein-coding sequences except for the cytochrome oxidase subunit I (COI) start with a typical ATN codon. Instead, CCG, which is a rare but possible initiation codon, is located at the initiation context of COI. The usual termination codons (TAA and TAG) were found from 12 PCGs. However, the ND5 had incomplete termination codon (T). All tRNA genes could be folded into the typical cloverleaf secondary structure, excluding tRNA ^Ser(AGN) which forms another structure according to the Steinberg–Cedergren tertiary structure. The sizes of the large and small ribosomal RNA genes are 1319 and 848 bp, respectively. The A + T content of the A + T-rich region is 78.7%, which is the lowest among the known mitogenome of Acridoidea.     

Complete nucleotide sequence and organization of the mitochondrial genome of eri-silkworm,Samia cynthia ricini (Lepidoptera: Saturniidae)

Samia cynthia ricini is a commercial silk-producing insect that is now reared year-round in Korea, with the expectation of being utilized for diverse purposes. In this report, we present the complete mitochondrial genome (mitogenome) of S. c. ricini. The 15,384-bp long S. cynthia ricini mitogenome was amplified into 26 short fragments using three long overlapping fragments using primers designed from reported lepidopteran mitogenome sequences. The genome comprises 37 genes (13 protein-coding genes, two rRNA genes, and 22 tRNA genes), and one large non-coding region termed the A + T-rich region. The A/T content of the third codon position was 91.7%, which was 18.8% and 21.6% higher than those of first and second codon positions, respectively. The high A/T content in the genome is reflected in codon usage, accounting for 39.5% of A/T-composed codons (TTA, ATT, TTT, and ATA). Unlike a previous report on the start codon for the COI gene, the S. c. ricini COI gene commences with a typical ATT codon. A total of 221 bp of non-coding sequences are dispersed in 17 regions, ranging in size from 1 to 54 bp, which comprise 1.4% of the total genome. One of the non-coding sequence located between tRNA^Gln and ND2 (54 bp) has 77% sequence homology with the 5′-sequence of the neighboring ND2 gene, suggesting partial duplication of the sequence during evolution. The 361-bp long A + T-rich region contains an 18 bp-long poly-T stretch, ATAGA motif, ATTTA element, microsatellite-like A/T sequence, poly-A stretch and one tRNA-like sequence, as typically found in Lepidoptera including Bombycoidea.     

The complete mitochondrial genome of the Ailanthus silkmoth, Samia cynthia cynthia (Lepidoptera: Saturniidae)

The complete mitochondrial genome (mitogenome) of the Ailanthus silkmoth, Samia cynthia cynthia (Lepidoptera: Saturniidae) was determined. The circular genome is 15,345 bp long, and presents a typical gene organization and order for sequenced mitogenomes of Bombycidea species. The nucleotide composition of the genome is highly A+T biased, accounting for 79.86%. The AT skew of the genome is slightly negative, indicating the occurrence of more Ts than As, as found in other Saturniidae species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI and COII, which are tentatively designated by CGA and GTG, respectively, as observed in other insects. Four of 13 PCGs, including COI, COII, ATP6, and ND3, harbor the incomplete termination codons, T or TA. With an exception for tRNA^Ser(AGN), all other tRNAs can form a typical clover-leaf structure of mitochondrial tRNA. The 359 bp A+T-rich region of S. c. cynthia contains non-repetitive sequences, but harbors several features common to the Bombycidea insects, including the motif ATAGA followed by a poly-T stretch of 19 bp, a microsatellite-like (AT)₇ element preceded by the ATTTA motif, and a poly-A element upstream tRNA^Met. The phylogenetic analyses support the morphology-based current hypothesis that Bombycidae and Saturniidae are monophyletic. Our result confirms that Saturniini and Attacini form a reciprocal monophyletic group within Saturniidae.     

Complete mitochondrial genome sequence of the yellow-spotted long-horned beetle <Emphasis Type="Italic">Psacothea hilaris</Emphasis> (Coleoptera: Cerambycidae) and phylogenetic analysis among coleopteran insects

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 tRNA^Ser(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 tRNA^Arg-like sequence and one tRNA^Lys-like sequence. 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 mitochondrial genome of the leafminer <Emphasis Type="Italic">Liriomyza trifolii</Emphasis> (Diptera: Agromyzidae)

Liriomyza trifolii (Diptera: Agromyzidae) is one of the most economically significant pests in the world. In this paper we present sequence data for the complete mitochondrial genome of L. trifolii. The circular genome is 16,141 bp long and contains one encoding region including 37 genes and one non-coding A+T-rich region. Gene numbers and organization is similar to that of the typical insect mitochondrial genomes except that two additional tRNA genes are found in the A+T-rich region (tRNA^Thr and tRNA^Leu(UUR)). All of the protein initiation codons are ATN, except ND1 which begins with GTG and COI which is initiated by the quadruplet ATCA. The 22 tRNA anticodons of L. trifolii match those observed in Drosophila yakuba, and all of tRNAs form the typical cloverleaf structure except for tRNA^Ser(AGN), which has lost the DHU-arm. The A+T-rich region of L. trifolii also contains two previously noted Diperan features—a highly conserved polyT stretch and a (TA)_n stretch.     

Characterization of the complete mitochondrial genome of the giant silkworm moth, Eriogyna pyretorum (Lepidoptera: Saturniidae)

The complete mitochondrial genome (mitogenome) of Eriogyna pyretorum (Lepidoptera: Saturniidae) was determined as being composed of 15,327 base pairs (bp), including 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The arrangement of the PCGs is the same as that found in the other sequenced lepidopteran. The AT skewness for the E. pyretorum mitogenome is slightly negative (-0.031), indicating the occurrence of more Ts than As. The nucleotide composition of the E. pyretorum mitogenome is also biased toward A + T nucleotides (80.82%). All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 and 2 (cox1 and cox2). Two of the 13 PCGs harbor the incomplete termination codon by T. All tRNA genes have a typical clover-leaf structure of mitochondrial tRNA, with the exception of trnS1(AGN) and trnS2(UCN). Phylogenetic analysis among the available lepidopteran species supports the current morphology-based hypothesis that Bombycoidea, Geometroidea, Notodontidea, Papilionoidea and Pyraloidea are monophyletic. As has been previously suggested, Bombycidae (Bombyx mori and Bombyx mandarina), Sphingoidae (Manduca sexta) and Saturniidae (Antheraea pernyi, Antheraea yamamai, E. pyretorum and Caligula boisduvalii) formed a group.     

The complete mitochondrial genome of the fall webworm,Hyphantria cunea (Lepidoptera: Arctiidae)

The complete mitochondrial genome (mitogenome) of the fall webworm, Hyphantria cunea (Lepidoptera: Arctiidae) was determined. The genome is a circular molecule 15 481 bp long. It presents a typical gene organization and order for completely sequenced lepidopteran mitogenomes, but differs from the insect ancestral type for the placement of tRNA^Met. The nucleotide composition of the genome is also highly A + T biased, accounting for 80.38%, with a slightly positive AT skewness (0.010), indicating the occurrence of more As than Ts, as found in the Noctuoidea species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI, which is tentatively designated by the CGA codon as observed in other lepidopterans. Four of 13 PCGs harbor the incomplete termination codon, T or TA. All tRNAs have a typical clover-leaf structure of mitochondrial tRNAs, except for tRNA^Ser(AGN), the DHU arm of which could not form a stable stem-loop structure. The intergenic spacer sequence between tRNA^Ser(AGN) and ND1 also contains the ATACTAA motif, which is conserved across the Lepidoptera order. The H. cunea A+T-rich region of 357 bp is comprised of non-repetitive sequences, but harbors several features common to the Lepidoptera insects, including the motif ATAGA followed by an 18 bp poly-T stretch, a microsatellite-like (AT)₈ element preceded by the ATTTA motif, an 11 bp poly-A present immediately upstream tRNA^Met. The phylogenetic analyses support the view that the H. cunea is closerly related to the Lymantria dispar than Ochrogaster lunifer, and support the hypothesis that Noctuoidea (H. cunea, L. dispar, and O. lunifer) and Geometroidea (Phthonandria atrilineata) are monophyletic. However, in the phylogenetic trees based on mitogenome sequences among the lepidopteran superfamilies, Papillonoidea (Artogeia melete, Acraea issoria, and Coreana raphaelis) joined basally within the monophyly of Lepidoptera, which is different to the traditional classification.     

The complete mitochondrial genome of the yellow coaster, <Emphasis Type="Italic">Acraea</Emphasis> <Emphasis Type="Italic">issoria</Emphasis> (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini): sequence,gene organization and a unique tRNA translocation event

In this paper, the complete mitochondrial genome of Acraea issoria (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini) is reported; a circular molecule of 15,245 bp in size. For A. issoria, genes are arranged in the same order and orientation as the complete sequenced mitochondrial genomes of the other lepidopteran species, except for the presence of an extra copy of tRNA^Ile(AUR)b in the control region. All protein-coding genes of A. issoria mitogenome start with a typical ATN codon and terminate in the common stop codon TAA, except that COI gene uses TTG as its initial codon and terminates in a single T residue. All tRNA genes possess the typical clover leaf secondary structure except for tRNA^Ser(AGN), which has a simple loop with the absence of the DHU stem. The sequence, organization and other features including nucleotide composition and codon usage of this mitochondrial genome were also reported and compared with those of other sequenced lepidopterans mitochondrial genomes. There are some short microsatellite-like repeat regions (e.g., (TA)₉, polyA and polyT) scattered in the control region, however, the conspicuous macro-repeats units commonly found in other insect species are absent.     

The complete mitochondrial genome of the leafminer Liriomyza sativae (Diptera: Agromyzidae): great difference in the A+T-rich region compared to Liriomyza trifolii

The complete mitochondrial genome sequence of Liriomyza sativae Blanchard (15,551 bp) was determined and analyzed in this study. The circular genome contained 37 genes including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and an A + T-rich region. The initiation codons of COI and ND1 were ‘ATCA’ and ‘GTG’, respectively. ND2 gene used the truncated termination codon ‘T’. All the tRNA genes had the typical cloverleaf secondary structures except for tRNA^Ser(AGN) gene, which was found with the absence of a DHU arm. In addition, a tRNA-like secondary structure (tRNA^Met) was found in the A + T-rich region. The great difference was that the length of L. sativae A + T-rich region was 597 bp shorter than that of Liriomyza trifolii (Burgess). Meanwhile, some minor differences such as ‘TATA’ block were also observed in L. sativae in contrast to ‘TACA’ block in L. trifolii. There were also some essential structure elements such as ‘TATA’ block, ‘G(A)_nT’ block, poly-T stretch and stem-and-loop structure in the A + T-rich region of L. sativae mitochondrial genome.     

The complete mitochondrial genome of the endangered Apollo butterfly,Parnassius apollo (Lepidoptera: Papilionidae) and its comparison to other Papilionidae species

The Apollo butterfly, Parnassius apollo is a representative species of the butterfly subfamily Parnassiinae. This charming species is one of the most endangered butterfly species in the world. In this study, we sequenced its complete mitochondrial genome (mitogenome), with the aim of accumulating genetic information for further studies of population genetics and mitogenome evolution in the Papilionidae. The 15,404-bp long mitogenome harbors a typical set of 37 genes and is the largest butterfly mitogenome determined, except for Papilio maraho (16,094 bp). Like many other sequenced lepidopteran species, one tRNA^Trp-like and one tRNA^Leu(UUR)-like sequences were detected in the AT-rich region. A total of 164 bp of non-coding sequences are dispersed in 14 regions throughout the genome. The longest intergenic spacer (68 bp) is located between tRNA^Ser(AGN) and tRNA^Glu, and is the largest spacer at this location among Papilionidae species. This spacer may have resulted from an 8-fold repetition of a TTTCTTCT motif or a 4-fold repetition of a CTTTATTT motif.