The complete mitochondrial genome of the rice moth, Corcyra cephalonica

Comparative mitochondrial genomic analysis provides new insights into the evolution of the subfamily Lamiinae (Coleoptera: Cerambycidae)

The genus Monochamus within the subfamily Lamiinae is the main vector of Bursaphelenchus xylophilus, which causes pine wilt disease and induces substantial economic and ecological losses. Only three complete mitochondrial genomes of the genus Monochamus have been sequenced to date, and no comparative mitochondrial genomic studies of Lamiinae have been conducted. Here, the mitochondrial genomes of two Monochamus species, M. saltuarius and M. urussovi, were newly sequenced and annotated. The composition and order of genes in the mitochondrial genomes of Monochamus species are conserved. All transfer RNAs exhibit the typical clover-leaf secondary structure, with the exception of trnS1. Similar to other longhorn beetles, Lamiinae mitochondrial genomes have an A + T bias. All 13 protein-coding genes have experienced purifying selection, and tandem repeat sequences are abundant in the A + T-rich region. Phylogenetic analyses revealed congruent topologies among trees inferred from the five datasets, with the monophyly of Acanthocinini, Agapanthiini, Batocerini, Dorcaschematini, Pteropliini, and Saperdini receiving high support. The findings of this study enhance our understanding of mitochondrial genome evolution and will provide a basis for future studies of population genetics and phylogenetic investigations in this group.

Characterization of four mitochondrial genomes of Crambidae (Lepidoptera, Pyraloidea) and phylogenetic implications

Loxostege turbidalis, Loxostege aeruginalis, Pyrausta despicata, and Crambus perlellus belong to Crambidae, Pyraloidea. Their mitochondrial genomes (mitogenomes) were successfully sequenced. The mitogenomes of L. turbidalis, L. aeruginalis, P. despicata, and C. perlellus are 15 240 bp, 15 339 bp, 15 389 bp, and 15 440 bp. The four mitogenomes all have a typical insect mitochondrial gene order, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and one A + T rich region (control region). The PCGs are initiated by the typical ATN codons, except CGA for the cox1 gene. Most PCGs terminate with common codon TAA or TAG, the incomplete codon T is found as the stop codon for cox2, nad4, and nad5. Most tRNA genes exhibit typical cloverleaf structure, except trnS1 (AGN) lacking the dihydrouridine (DHU) arm. The secondary structure of rRNA of four mitogenomes were predicted. Poly-T structure and micro-satellite regions are conserved in control regions. The phylogenetic analyses based on 13 PCGs showed the relationships of subfamilies in Pyraloidea. Pyralidae, and Crambidae are monophyletic, respectively. Pyralidae comprises four subfamilies, which form the following topology with high support values: (Galleriinae + ((Pyralinae + Epipaschiinae)+ Phycitinae)). Crambidae includes seven subfamilies and is divided into two lineages. Pyraustinae and Spilomelinae are sister groups of each other, and form the "PS clade." Other five subfamilies (Crambinae, Acentropinae, Scopariinae, Schoenobiinae, and Glaphyriinae) form the "non-PS clade" in the Bayesian inference tree. However, Schoenobiinae is not grouped with the other four subfamilies and located at the base of Crambidae in two maximum likelihood trees.

Comparative mitogenomic and evolutionary analysis of Lycaenidae (Insecta: Lepidoptera): Potential association with high-altitude adaptation

Harsh environments (e.g., hypoxia and cold temperatures) of the Qinghai-Tibetan Plateau have a substantial influence on adaptive evolution in various species. Some species in Lycaenidae, a large and widely distributed family of butterflies, are adapted to the Qinghai-Tibetan Plateau. Here, we sequenced four mitogenomes of two lycaenid species in the Qinghai-Tibetan Plateau and performed a detailed comparative mitogenomic analysis including nine other lycaenid mitogenomes (nine species) to explore the molecular basis of high-altitude adaptation. Based on mitogenomic data, Bayesian inference, and maximum likelihood methods, we recovered a lycaenid phylogeny of [Curetinae + (Aphnaeinae + (Lycaeninae + (Theclinae + Polyommatinae)))]. The gene content, gene arrangement, base composition, codon usage, and transfer RNA genes (sequence and structure) were highly conserved within Lycaenidae. TrnS1 not only lacked the dihydrouridine arm but also showed anticodon and copy number diversity. The ratios of non-synonymous substitutions to synonymous substitutions of 13 protein-coding genes (PCGs) were less than 1.0, indicating that all PCGs evolved under purifying selection. However, signals of positive selection were detected in cox1 in the two Qinghai-Tibetan Plateau lycaenid species, indicating that this gene may be associated with high-altitude adaptation. Three large non-coding regions, i.e., rrnS-trnM (control region), trnQ-nad2, and trnS2-nad1, were found in the mitogenomes of all lycaenid species. Conserved motifs in three non-coding regions (trnE-trnF, trnS1-trnE, and trnP-nad6) and long sequences in two non-coding regions (nad6-cob and cob-trnS2) were detected in the Qinghai-Tibetan Plateau lycaenid species, suggesting that these non-coding regions were involved in high-altitude adaptation. In addition to the characterization of Lycaenidae mitogenomes, this study highlights the importance of both PCGs and non-coding regions in high-altitude adaptation.

Nine Mitochondrial Genomes of the Pyraloidea and Their Phylogenetic Implications (Lepidoptera)

Simple Summary

The Pyraloidea is a large superfamily of Lepidoptera in species composition. To date, the higher-level phylogenetic relationships in this group remain unresolved, and many taxa, with taxonomic positions historically established by morphological characters, need to be confirmed through sequencing of DNA, including mitochondrial genome sequences (mitogenomes). Here, we newly generated nine complete mitogenomes for Pyraloidea that shared identical gene content, and arrangements that are typical of Lepidoptera. The current phylogenetic results confirmed previous multilocus studies, indicating the effectiveness of mitogenomes for inference of Pyraloidea higher-level relationships. Unexpectedly, Orybina Snellen was robustly placed as basal to the remaining Pyralidae taxa, rather than nested in the Pyralinae of Pyralidae as morphologically defined and placed. Our results bring a greater understanding to Pyraloidea phylogeny, and highlight the necessity of sequencing more pyraloid taxa to reevaluate their phylogenetic positions.

Abstract

The Pyraloidea is one of the species-rich superfamilies of Lepidoptera and contains numerous economically important pest species that cause great loss in crop production. Here, we sequenced and annotated nine complete mitogenomes for Pyraloidea, and further performed various phylogenetic analyses, to improve our understanding of mitogenomic evolution and phylogeny of this superfamily. The nine mitogenomes were circular, double-stranded molecules, with the lengths ranging from 15,214 bp to 15,422 bp, which are comparable to other reported pyraloid mitogenomes in size. Gene content and arrangement were highly conserved and are typical of Lepidoptera. Based on the hitherto most extensive mitogenomic sampling, our various resulting trees showed generally congruent topologies among pyraloid subfamilies, which are almost in accordance with previous multilocus studies, indicating the suitability of mitogenomes in inferring high-level relationships of Pyraloidea. However, nodes linking subfamilies in the “non-PS clade” were not completely resolved in terms of unstable topologies or low supports, and future investigations are needed with increased taxon sampling and molecular data. Unexpectedly, Orybina Snellen, represented in a molecular phylogenetic investigation for the first time, was robustly placed as basal to the remaining Pyralidae taxa across our analyses, rather than nested in Pyralinae of Pyralidae as morphologically defined. This novel finding highlights the need to reevaluate Orybina monophyly and its phylogenetic position by incorporating additional molecular and morphological evidence.

The First Mitogenomes of the Subfamily Odontiinae (Lepidoptera, Crambidae) and Phylogenetic Analysis of Pyraloidea

Simple Summary

The Odontiinae is a small group in the Pyraloidea comprised of 388 species in 88 genera, but externally, these moths are diverse, including heterogeneous maculation and a size range from 9 to 50 mm in total wingspan. The monophyly of Pyraloidea and the two families (Pyralidae and Crambidae) is well supported by phylogenetic analyses based on morphology and molecular data of multiple nuclear genes. However, only a few mito-phylogenetic analyses have been conducted and no mitogenome of Odontiinae species has been reported. Three complete mitogenomes of odontiine species were sequenced and analyzed for the first time herein. The results showed that Odontiinae mitogenomes shared similar genomic characters with other Pyraloidea. The phylogenetic analyses based on 13 PCGs of mitogenomes confirmed the monophyly of Odontiinae and its position within Crambidae.

Abstract

The complete mitochondrial genomes of three species of Odontiinae were newly sequenced: Dausara latiterminalis Yoshiyasu, Heortia vitessoides (Moore), and Pseudonoorda nigropunctalis (Hampson). These circular and double-stranded mitogenomes vary from 15,084 bp to 15,237 bp in size, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs) and an A + T-rich region. The nucleotide composition indicated a strong A/T bias. Most PCGs are initiated with an ATN codon and terminated by a codon of TAR. All tRNAs could be folded into the clover-leaf structure with the exception of trnS1 (AGN), in which the dihydrouridine (DHU) arm formed a simple loop, and the motif ‘ATAG’ and ‘ATTTA’ in the A + T-rich region was also founded. The phylogenomic analyses covering Odontiinae + 11 subfamilies of Pyraloidea were conducted. Similar topologies were generated from both Bayesian inference (BI) and maximum likelihood (ML) analyses based on the nucleotide and amino acid sequence data. There was some discrepancy in the sister-group relationship of Odontiinae and Glaphyriinae, and the relationships among the subfamilies in the ‘CAMMSS clade’ of the Crambidae. The results of this study suggest that mitogenomic data are useful for resolving the deep-level relationships of Pyraloidea and the topologies generated from amino acid data might be more realistic and reliable. Moreover, more mitogenomic taxon sampling and larger scale analyses with more genes or a combination of mitogenomic and nuclear genes are needed to reconstruct a comprehensive framework of the pyraloid phylogeny.

Characterization of the complete mitochondrial genome of Orthaga olivacea Warre (Lepidoptera Pyralidae) and comparison with other Lepidopteran insects

Orthaga olivacea Warre (Lepidoptera: Pyralidae) is an important agricultural pest of camphor trees (Cinnamomum camphora). To further supplement the known genome-level features of related species, the complete mitochondrial genome of Orthaga olivacea is amplified, sequenced, annotated, analyzed, and compared with 58 other species of Lepidopteran. The complete sequence is 15,174 bp, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a putative control region. Base composition is biased toward adenine and thymine (79.02% A+T) and A+T skew are slightly negative. Twelve of the 13 PCGs use typical ATN start codons. The exception is cytochrome oxidase 1 (cox1) that utilizes a CGA initiation codon. Nine PCGs have standard termination codon (TAA); others have incomplete stop codons, a single T or TA nucleotide. All the tRNA genes have the typical clover-leaf secondary structure, except for trnS(AGN), in which dihydrouridine (DHU) arm fails to form a stable stem-loop structure. The A+T-rich region (293 bp) contains a typical Lepidopter motifs 'ATAGA' followed by a 17 bp poly-T stretch, and a microsatellite-like (AT)13 repeat. Codon usage analysis revealed that Asn, Ile, Leu2, Lys, Tyr and Phe were the most frequently used amino acids, while Cys was the least utilized. Phylogenetic analysis suggested that among sequenced lepidopteran mitochondrial genomes, Orthaga olivacea Warre was most closely related to Hypsopygia regina, and confirmed that Orthaga olivacea Warre belongs to the Pyralidae family.

The first mitochondrial genome for Phaudidae (Lepidoptera) with phylogenetic analyses of Zygaenoidea

Phauda flammans Walker belongs to Phaudidae (Lepidoptera), which is a holometabolous and leaf-eating pest that harms trees. So far, there is no mitochondrial (mt) genome reported of Phaudidae. Herein, we sequenced and annotated the complete mt genome of P. flammans representing the first mt genome of Phaudidae and predicted the secondary structures of its RNAs in this study. This mt genome is 15470 bp long consisting of 13 protein coding genes (PCGs), 22 tRNAs, 2 rRNAs and the control region, which are usually conserved in insects. Most PCGs used the standard ATN start codons and complete TAA/TAG termination codons. Almost all of tRNA genes exhibited cloverleaf secondary structures except that the dihydorouridine (DHU) arm of tRNA^Ser(AGN) was absent. The phylogenetic analyses using both Bayesian inference (BI) and maximum likelihood (ML) methods all supported that Phaudidae was a single family being the sister group to Zygaenidae. More mt genomes are needed to better understand the phylogenetic relationships within Zygaenoidea in the future.

Gene characteristics of the complete mitochondrial genomes of Paratoxodera polyacantha and Toxodera hauseri (Mantodea: Toxoderidae)

The family Toxoderidae (Mantodea) contains an ecologically diverse group of praying mantis species that have in common greatly elongated bodies. In this study, we sequenced and compared the complete mitochondrial genomes of two Toxoderidae species, Paratoxodera polyacantha and Toxodera hauseri, and compared their mitochondrial genome characteristics with another member of the Toxoderidae, Stenotoxodera porioni (KY689118). The lengths of the mitogenomes of T. hauseri and P. polyacantha were 15,616 bp and 15,999 bp, respectively, which is similar to that of S. porioni (15,846 bp). The size of each gene as well as the A+T-rich region and the A+T content of the whole genome were also very similar among the three species as were the protein-coding genes, the A+T content and the codon usages. The mitogenome of T. hauseri had the typical 22 tRNAs, whereas that of P. polyacantha had 26 tRNAs including an extra two copies of trnA-trnR. Intergenic regions of 67 bp and 76 bp were found in T. hauseri and P. polyacantha, respectively, between COX2 and trnK; these can be explained as residues of a tandem duplication/random loss of trnK and trnD. This non-coding region may be synapomorphic for Toxoderidae. In BI and ML analyses, the monophyly of Toxoderidae was supported and P. polyacantha was the sister clade to T. hauseri and S. porioni.

Complete mitochondrial genome of the soybean leaffolder, Omiodes indicata (Lepidoptera: Pyraloidea: Crambidae), and phylogenetic analysis for Pyraloidea

Complete mitochondrial genome (mitogenome) of the Omiodes indicata was sequenced and characterized. The circular mitogenome is 15,367bp long, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A+T-rich region. Nucleotide composition is highly biased toward A+T nucleotides (81.6%). All 13 PCGs initiate with canonical start codon (ATN), except for cox1 that initiates with CGA. All tRNAs have a typical clover-leaf structure, except for trnS1 (AGN) in which the base pairs of the dihydrouridine (DHU) arm are reduced. In O. indicata, the motifs "ATGATAA" and "ATACTAA" between atp8 and atp6, trnS2 and nad1, respectively, and the motifs "ATAG" and "ATTTA" in the A+T-rich region can be identified. Comparative phylogenetic analyses based on four datasets show that the dataset including all coding positions of 13 PCGs exhibit the highest informativeness in resolving higher phylogeny of Pyraloidea. Bayesian inference (BI) and maximum likelihood (ML) analyses yield generally well-supported phylogenetic relationships among the eleven pyraloid subfamilies involved. However, the relationships among the five subfamilies (Acentropinae, Crambinae, Glaphyriinae, Schoenobiinae and Scopariinae) in ML analysis are ambiguous, which might be resolved by ample sampling in future mitogenome-based phylogenetic studies.

Complete sequence and characterization of the Ectropis oblique mitochondrial genome and its phylogenetic implications

In the present study, we sequenced the entire mitochondrial genome (mitogenome) of Ectropis oblique using PCR amplification and sequencing methods The entire mitogenome is 15,356bp long, including 13 protein coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and A+T rich element. The base composition and gene arrangement are identical to those of other lepidopterans. Of these 37 genes, twenty three are resided on heavy strand, while fourteen located on light strand. The newly sequenced mitogenome displayed a biased Adenine/thymine 84.94% usage versus guanine/cytosine. The AT skewness is 0.030 and the GC skewness is -0.180. Twelve out of 13 PCGs initiated with canonical start codon (ATN), while cox1 started with CGA. The A+T rich element of E. oblique is 363bp long and contains of many features common to lepidopteran insects, including the 'ATAGA' motif, a 22bp poly (T) stretch and a microsatellite-like (AT)₈ element upstream of trnM. Further, this mitogenome comprises 15 intergenic spacer and overlapping regions. Phylogenetic analyses showed that E. oblique belongs to Geometridae, and that the monophyly of Lepidoptera superfamilies is well supported.

Comparative mitochondrial genome analysis of Daphnis nerii and other lepidopteran insects reveals conserved mitochondrial genome organization and phylogenetic relationships

In the present study, the complete sequence of the mitochondrial genome (mitogenome) of Daphnis nerii (Lepidoptera: Sphingidae) is described. The mitogenome (15,247 bp) of D.nerii encodes13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and an adenine (A) + thymine (T)-rich region. Its gene complement and order is similar to that of other sequenced lepidopterans. The 12 PCGs initiated by ATN codons except for cytochrome c oxidase subunit 1 (cox1) gene that is seemingly initiated by the CGA codon as documented in other insect mitogenomes. Four of the 13 PCGs have the incomplete termination codon T, while the remainder terminated with the canonical stop codon. This mitogenome has six major intergenic spacers, with the exception of A+T-rich region, spanning at least 10 bp. The A+T-rich region is 351 bp long, and contains some conserved regions, including 'ATAGA' motif followed by a 17 bp poly-T stretch, a microsatellite-like element (AT)9 and also a poly-A element. Phylogenetic analyses based on 13 PCGs using maximum likelihood (ML) and Bayesian inference (BI) revealed that D. nerii resides in the Sphingidae family.

The mitochondrial genome of the lepidopteran host cadaver (Thitarodes sp.) of Ophiocordyceps sinensis and related phylogenetic analysis

To understand the phylogeny of the host insect (Thitarodes sp.) of the fungus Ophiocordyceps sinensis, we sequenced, annotated and characterized the complete mitochondrial (mt) genome of the host cadaver of a natural O. sinensis. Further, we compared the Thitarodes sp. mt genome with those of the other 7 sequenced Hepialidae and examined the phylogenetic relationships using a constructed Maximum Likelihood (ML) phylogenetic tree and mt genomic features (genetic distances and intergenic spacers). The mt genome is a circular molecule of 16,280bp in length with a high A+T content (81.20%) and contains 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and an AT-rich region. The gene arrangement is identical to the ancestral arrangement but differs from those of other lepidopteran mt genomes because of the arrangement of tRNA genes. The tRNA region, which is located between the AT-rich region and nad2, is trnI/trnQ/trnM (IQM) in Thitarodes sp., rather than the trnM/trnI/trnQ (MIQ) of the Lepidoptera-specific rearrangement. All PCGs begin with the canonical start codons ATN or NTG, except for cox1, which starts with CGA. Most PCGs terminate with the typical stop codon TAA, although some have an incomplete stop codon (T). The 1473bp AT-rich region is located between the rrnS (12S rRNA) and trnI, which is the longest sequenced in a Thitarodes mt genome to date, containing nine 112bp copies and one partial copy of a 55bp sequence. The results derived from the phylogenetic tree, the genetic distances and the intergenic spacers of the mt genome show that the host insect of O. sinensis belongs to the Thitarodes, while Endoclita signifer and Napialus hunanensis form a relatively distinct lineage from Thitarodes. The sequence and full annotation of this moth mt genome will provide more molecular information about the Exoporia within the Lepidoptera, and the clarification of its phylogeny will improve the management of this insect resource and the conservation and sustainable use of this endangered medicinal species in China.

Mitochondrial genome characterization of Tecia solanivora (Lepidoptera: Gelechiidae) and its phylogenetic relationship with other lepidopteran insects

The complete mitogenome of the potato tuber moth Tecia solanivora (Lepidoptera: Gelechiidae) was sequenced, annotated, characterized and compared with 140 species of the order Lepidoptera. The circular genome is 15,251 bp, containing 37 genes (13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A+T-rich region). The gene arrangement was identical to other lepidopteran mitogenomes but different from the ancestral arrangement found in most insects for the tRNA-Met gene (A+T-region, tRNA-I, tRNA-Q, tRNA-M). The mitogenome of T. solanivora is highly A+T-biased (78.2%) and exhibits negative AT- and GC-skews. All PCGs are initiated by canonical ATN start codons, except for Cytochrome Oxidase subunit 1 (COI), which is initiated by CGA. Most PCGs have a complete typical stop codon (TAA). Only NAD1 has a TAG stop codon and the COII and NAD5 genes have an incomplete stop codon consisting of just a T. The A+T-rich region is 332 bp long and contains common features found in lepidopteran mitogenomes, including the 'ATAGA' motif, a 17 bp poly (T) stretch and a (AT)8 element preceded by the 'ATTTA' motif. Other tandem repeats like (TAA)4 and (TAT)7 were found, as well as (T)6 and (A)10 mononucleotide repeat elements. Finally, this mitogenome has 20 intergenic spacer regions. The phylogenetic relationship of T. solanivora with 28 other lepidopteran families (12 superfamilies) showed that taxonomic classification by morphological features coincides with the inferred phylogeny. Thus, the Gelechiidae family represents a monophyletic group, suggesting that T. solanivora and Pectinophora gossypiella have a recent common ancestor.

The complete mitochondrial genome of Plodia interpunctella (Lepidoptera: Pyralidae) and comparison with other Pyraloidea insects

The mitochondrial (mt) genome can provide important information for the understanding of phylogenetic relationships. The complete mt genome of Plodia interpunctella (Lepidoptera: Pyralidae) has been sequenced. The circular genome is 15 287 bp in size, encoding 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The AT skew of this mt genome is slightly negative, and the nucleotide composition is biased toward A+T nucleotides (80.15%). All PCGs start with the typical ATN (ATA, ATC, ATG, and ATT) codons, except for the cox1 gene which may start with the CGA codon. Four of the 13 PCGs harbor the incomplete termination codon T or TA. All the tRNA genes are folded into the typical clover-leaf structure of mitochondrial tRNA, except for trnS1 (AGN) in which the DHU arm fails to form a stable stem-loop structure. The overlapping sequences are 35 bp in total and are found in seven different locations. A total of 240 bp of intergenic spacers are scattered in 16 regions. The control region of the mt genome is 327 bp in length and consisted of several features common to the sequenced lepidopteran insects. Phylogenetic analysis based on 13 PCGs using the Maximum Likelihood method shows that the placement of P. interpunctella was within the Pyralidae.

The complete mitochondrial genomes of three bristletails (Insecta: Archaeognatha): the paraphyly of Machilidae and insights into archaeognathan phylogeny

The order Archaeognatha was an ancient group of Hexapoda and was considered as the most primitive of living insects. Two extant families (Meinertellidae and Machilidae) consisted of approximately 500 species. This study determined 3 complete mitochondrial genomes and 2 nearly complete mitochondrial genome sequences of the bristletail. The size of the 5 mitochondrial genome sequences of bristletail were relatively modest, containing 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. The gene orders were identical to that of Drosophila yakuba and most bristletail species suggesting a conserved genome evolution within the Archaeognatha. In order to estimate archaeognathan evolutionary relationships, phylogenetic analyses were conducted using concatenated nucleotide sequences of 13 protein-coding genes, with four different computational algorithms (NJ, MP, ML and BI). Based on the results, the monophyly of the family Machilidae was challenged by both datasets (W12 and G12 datasets). The relationships among archaeognathan subfamilies seemed to be tangled and the subfamily Machilinae was also believed to be a paraphyletic group in our study.

The complete mitochondrial genome of Acleris fimbriana (Lepidoptera: Tortricidae)

The yellow tortrix, Acleris fimbriana belongs to Tortricidae in Lepidoptera. We described the complete mitogenome of A. fimbriana, which is typical circular duplex molecules and 15,933 bp in length containing the standard metazoan set of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and an A + T-rich region with macro-repeat sequences. All the inferred tRNA secondary structures show the common cloverleaf pattern, with the exception of trnS1(AGN) which lacks the DHU arm. The A. fimbriana mitochondrial genome has the same gene order with other lepidopterans.