Characterization of the complete mitochondrial genome of Bombyx mori strain H9 (Lepidoptera: Bombycidae)

The complete mitochondrial genome (mitogenome) of Bombyx mori strain H9 (Lepidoptera: Bombycidae) is 15,670 base pairs (bp) in length, encoding 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The nucleotide composition of the genome is highly A + T biased, accounting for 81.31%, with a slightly positive AT skewness (0.059). The arrangement of 13 PCGs is similar to that of other sequenced lepidopterans. All the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is proposed by the TTAG sequence as observed in other lepidopterans. Unlike the other PCGs, the cox1 and cytochrome c oxidase subunit 2 (cox2) genes have incomplete stop codons consisting of just a T. All tRNAs have typical structures of insect mitochondrial tRNAs, which is different from other sequenced lepidopterans. The structure of A + T-rich region is similar to that of other sequenced lepidopterans, including non-repetitive sequences, the ATAGA binding domain, a 18 bp poly-T stretch and a poly-A element upstream of transfer RNA M (trnM) gene. Phylogenetic analysis shows that the domesticated silkmoth B. mori originated from the Chinese Bombyx mandarina.     

The complete mitogenome of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) and comparison with other lepidopteran insects

The complete mitochondrial genome (mitogenome) of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) was determined to be 15,653 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a A + T-rich region. It has the typical gene organization and order of mitogenomes from lepidopteran insects. The AT skew of this mitogenome was slightly positive and the nucleotide composition was also biased toward A + T nucleotides (81.31%). All PCGs were initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which was initiated by CGA. The cox1 and cox2 genes had incomplete stop codons consisting of just a T. All the tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome was 495 bp in length and consisted of several features common to the lepidopteras. Phylogenetic analysis showed that the B. mori Dazao was close to Bombycidae.     

The complete mitochondrial genome of Gomphocerus tibetanus Uvarov, 1935 (Orthoptera: Acrididae: Gomphocerinae)

The complete nucleotide sequence of the mitochondrial genome (mitogenome) of Gomphocerus tibetanus Uvarov, 1935 (Orthoptera: Acrididae: Gomphocerinae) was determined. It is 15,571 bp in length and contains 74.8% A + T. All Gomphocerus tibetanus protein-coding sequences start with a typical ATN codon. The usual termination codons (TAA and TAG) were found from 13 PCGs except COI and COII which took incomplete codon T as termination codons. All tRNA genes could be folded into the typical cloverleaf secondary structure, except tRNA^Ser(AGN) lacking of dihydrouridine (D) arm. The sizes of the large and small ribosomal RNA genes are 1313 and 822 bp, respectively. The A + T content of the A + T-rich region is 82.3%. A preliminary analysis on characteristics of Gomphocerinae mitogenome was made by comparision among three Gomphocerinae mitogenomes and Locusta migratoria.     

The complete mitochondrial genome of fall armyworm <Emphasis Type="Italic">Spodoptera frugiperda</Emphasis> (Lepidoptera:Noctuidae)

The mitochondrial genome (mitogenome) can provide important information for understanding molecular evolution and phylogenetic analyses. The complete mitogenome of Spodoptera frugiperda (Lepidoptera:Noctuidae) was determined to be 15,365 bp in length and has the typical gene order found in Noctuidae mitogenomes, it includes 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a A+T-rich region. The nucleotide composition was biased toward A+T nucleotides (81.09 %) and the AT skew of this mitogenome was slightly positive (0.004). All PCGs were initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which was initiated by CGA. Eight of the 13 PCGs have the incomplete termination codon, T or TA. All the tRNA genes displayed the typical clover-leaf structure of mitochondrial tRNAs, with the exception of trnS1 (AGN). The A+T-rich region was 328 bp in length and consisted of several features common to the Noctuidae insects. Phylogenetic analysis showed that the S. frugiperda was within the Noctuidae.     

The complete mitochondrial genomes of three grasshoppers, Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis (Orthoptera: Pamphagidae)

The complete mitogenomes of Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis are 15,660 bp, 15,657 bp and 15,661 bp in size, respectively. All three mitogenomes contain a standard set of 13 protein - coding genes, 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A + T-rich region in the same order as those of the other analysed caeliferan species, including the rearrangement of trnAsp and trnLys. The putative initiation codon for the cox1 gene in the three species is CCG. The long polythymine stretch (T-stretch) in the A + T-rich region of the three species is not adjacent to the trnIle but inside the stem–loop sequence in the majority strand. The mitogenomes of F. helanshanensis and P. rubimarginis have higher overall similarities. The characterization of the three mitogenomes will enrich our knowledge on the Pamphagidae mitogenome. The phylogenetic analyses indicated that within the Caelifera, Pyrgomorphoidea is a sister group to Acridoidea. The species from the Pamphagidae form a monophyletic group, as is the case for Acrididae. Furthermore, the two families cluster as sister groups, supporting the monophyly of Acridoidea. The relationships among eight acridid subfamilies were (Cyrtacanthacridinae + (Calliptaminae + (Catantopinae + (Oxyinae + (Melanopline + (Acridinae + (Oedipodinae + Gomphocerinae))))))).     

Characterization of the complete mitochondrial genomes of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae)

The complete mitochondrial genomes (mitogenomes) of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae) were determined and analyzed. The circular genomes were 15,388 bp long for C. medinalis and 15,395 bp long for C. suppressalis. Both mitogenomes contained 37 genes, with gene order similar to that of other lepidopterans. Notably, 12 protein-coding genes (PCGs) utilized the standard ATN, but the cox1 gene used CGA as the initiation codon; the cox1, cox2, and nad4 genes in the two mitogenomes had the truncated termination codons T, T, and TA, respectively, but the nad5 gene was found to use T as the termination codon only in the C. medinalis mitogenome. Additionally, the codon distribution and Relative Synonymous Codon Usage of the 13 PCGs in the C. medinalis mitogenome were very different from those in other pyralid moth mitogenomes. Most of the tRNA genes had typical cloverleaf secondary structures. However, the dihydrouridine (DHU) arm of the trnS1(AGN) gene did not form a stable stem-loop structure. Forty-nine helices in six domains, and 33 helices in three domains were present in the secondary structures of the rrnL and rrnS genes of the two mitogenomes, respectively. There were four major intergenic spacers, except for the A+T-rich region, spanning at least 12 bp in the two mitogenomes. The A+T-rich region contained an 'ATAGT(A)'-like motif followed by a poly-T stretch in the two mitogenomes. In addition, there were a potential stem-loop structure, a duplicated 25-bp repeat element, and a microsatellite '(TA)(13)' observed in the A+T-rich region of the C. medinalis mitogenome. A poly-T motif, a duplicated 31-bp repeat element, and a 19-bp triplication were found in the C. suppressalis mitogenome. However, there are many differences in the A+T-rich regions between the C. suppressalis mitogenome sequence in the present study and previous reports. Finally, the phylogenetic relationships of these insects were reconstructed based on amino acid sequences of mitochondrial 13 PCGs using Bayesian inference and maximum likelihood methods. These molecular-based phylogenies support the traditional morphologically based view of relationships within the Pyralidae.     

The complete mitogenome of Apocheima cinerarius (Lepidoptera: Geometridae: Ennominae) and comparison with that of other lepidopteran insects

The complete mitochondrial genome (mitogenome) of a female flightless geometrid moth Apocheima cinerarius was found to be 15,722 bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a control region. The A + T content of the complete mitogenome is 80.83%. The AT skew value ([A − T] / [A + T]) is 0.027. The 13 PCGs of the mitogenome start with typical ATN codons, except for cox1 with the start codon CGA. All the tRNA genes have typical cloverleaf secondary structures, except for trnSer(AGN). The secondary structures of rrnL and rrnS were predicted. Six structural domains including conserved regions (IV, V) and variable regions (I, II, III, VI) were identified in the secondary structure of rrnL. The secondary structure of rrnS consists of 3 structural domains. The control region of A. cinerarius begins with conserved motifs of “ATAGA” + 19-bp poly T. It also contains a microsatellite-like (TA)₂₆, a stem-and-loop structure, and a poly-A stretch. Phylogenetic analysis showed that Geometroidea is more closely related to Bombycoidea than to Noctuoidea. A. cinerarius is more closely related to Biston panterinaria than to Phthonandria atrilineata, which is in accordance with the conventional morphology-based classification.     

The complete mitochondrial genome of Sasakia funebris (Leech) (Lepidoptera: Nymphalidae) and comparison with other Apaturinae insects

Sasakia funebris, a member of the lepidopteran family, Nymphalidae (superfamily Papilionoidea) is a rare species and is found only in some areas of South China. In this study, the 15,233 bp long complete mitochondrial genome of S. funebris was determined, and harbors the gene arrangement identical to all other sequenced lepidopteran insects. The nucleotide composition of the genome is highly A + T biased, accounting for 81.2%. All protein-coding genes (PCGs) start with typical ATN codons, except for COI which begins with the CGA codon. All tRNAs have a typical clover-leaf secondary structure, except for tRNA^Ser(AGN), the dihydrouridine (DHU) arm of which forms a simple loop. The S. funebris A + T-rich region of 370 bp contains several features common to the Lepidoptera insects, including the motif ATAGA followed by a 19 bp poly-T stretch, and two tandem repeats consisting of 18 bp repeat units and 14 bp repeat units. The phylogenetic analyses of Apaturinae based on mitogenome sequences showed: (S. funebris + Sasakia charonda) + (Apatura metis + Apatura ilia). This result is consistent with the morphological classification.     

Complete mitochondrial genome of the hemp borer,Grapholita delineana (Lepidoptera: Tortricidae): Gene variability and phylogeny among Grapholita

The mitochondrial genome (mitogenome) has been extensively used in phylogenetics and species-level evolutionary investigations. The lepidopteran family Tortricidae (leaf-roller moths), including the genus Grapholita, contains numerous species of economic importance, but for the majority of Grapholita species, their mitogenomes remain poorly studied. Here, we sequence and annotate the full mitogenome of Grapholita delineana, an important pest of hemp worldwide and compare it with the mitogenomes of two congeneric species available from GenBank. The G. delineana mitogenome is 15,599 bp long, including 37 typical mitochondrial genes and an A + T-rich region. Gene content, order and orientation are identical to other reported tortricid mitogenomes. Analyses of nucleotide diversity, Ka/Ks, genetic distance and number of variable sites together suggest that nad6 is the fastest-evolving gene among the mitochondrial PCGs of Grapholita. Our analyses indicate that Grapholita, as presently defined, is not monophyletic, confirming previous morphological and multiple-gene studies, using mitogenomic evidence. Our study provides information on comparative mitogenomics of Tortricidae especially Grapholita.     

The complete nucleotide sequence of the mitochondrial genome of <Emphasis Type="Italic">Phthonandria atrilineata</Emphasis> (Lepidoptera: Geometridae)

Using long-polymerase chain reaction (Long-PCR) method, we determined the complete nucleotide sequence of the mitochondrial genome (mitogenome) of Phthonandria atrilineata. The complete mtDNA from P. atrilineata was 15,499 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The P. atrilineata genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. The nucleotide composition of P. atrilineata mitogenome was biased toward A + T nucleotides (81.02%), and the 13 PCGs show different A + T contents that range from 73.25% (cox1) to 92.12% (atp8). Phthonandria had the canonical set of 22 tRNA genes, that fold in the typical cloverleaf structure described for metazoan mt tRNAs, with the unique exception of trnS(AGN). The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, which confirmed that P. atrilineata is most closely related to the superfamily Bombycoidea.     

Complete mitochondrial genome of Spilosoma lubricipedum (Noctuoidea: Erebidae) and implications for phylogeny of noctuid insects

Mitochondrial genomes (mitogenomes) have been widely used for studies on phylogenetic relationships and molecular evolutionary biology. Here, the complete mitogenome sequence of Spilosoma lubricipedum (Noctuoidea: Erebidae: Arctiinae) was determined (total length 15,375 bp) and phylogenetic analyses S. lubricipedum were inferred from available noctuid sequence data. The mitogenome of S. lubricipedum was found to be highly A + T-biased (81.39%) and exhibited negative AT- and GC-skews. All 13 protein-coding genes (PCGs) were initiated by ATN codons, except for cox1 with CGA. All tRNAs exhibited typical clover-leaf secondary structures, except for trnS1. The gene order of the S. lubricipedum mitogenome was trnM-trnI-trnQ-nad2. The A + T-rich region of S. lubricipedum contained several conservative features common to noctuid insects. Phylogenetic analysis within Noctuoidea was carried out based on mitochondrial data. Results showed that S. lubricipedum belonged to Erebidae and the Noctuoidea insects could be divided into five well-supported families (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))).     

The complete mitochondrial genome of the wild silkworm moth, Actias selene

The complete mitochondrial genome (mitogenome) of Actias selene (Lepidoptera: Saturniidae) was determined to be 15,236 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The arrangement of 13 PCGs was similar to that of other sequenced lepidopterans. The AT skew of the mitogenome of A. selene was slightly negative, indicating a higher number of T compared to A nucleotides. The nucleotide composition of the mitogenome of A. selene was also biased toward A+T nucleotides (78.91%). All PCGs were initiated by ATN codons, except for the gene encoding cytochrome c oxidase subunit 1 (cox1), which may be initiated by the TTAG, as observed in other lepidopterans. Three genes, including cox1, cox2, and nad5, had incomplete stop codons consisting of just a T. With an exception for trnS1(AGN), all the other tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A+T-rich region of the mitogenome of A. selene was 339 bp in length, and contains several features common to the Lepidopteras, including non-repetitive sequences, a conserved structure combining the motif ATAGA and an 18-bp poly-T stretch and a poly-A element upstream of trnM gene. Phylogenetic analysis showed that A. selene was close to Saturniidae.     

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 Eriocheir japonica sinensis (Decapoda: Varunidae) and its phylogenetic analysis

The complete mitochondrial genome (mitogenome) can provide novel insights into understanding the mechanisms underlying mitogenome evolution. In this study, the complete mitogenome of Eriocheir japonica sinensis (Decapoda: Varunidae) was determined to be 16,378 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a D-loop region. The AT skew of the E. j. sinensis mitogenome was slightly negative (−0.016), indicating a higher number of T compared with A nucleotides. The nucleotide composition of the mitogenome was also biased toward A + T nucleotides (71.6%). All PCGs were initiated by ATN codons. Eight of the 13 PCGs harbored the incomplete termination codon by T, or TA. All other tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The D-loop region of the E. j. sinensis mitogenome was 918 bp in length. Based on 13 PCGs, phylogenetic analysis confirmed the placement of E. j. sinensis within the Varunidae.     

Characterization of the complete mitochondrial genome of Tryporyza incertulas, in comparison with seven other Pyraloidea moths

The complete 15,223-bp mitochondrial genome (mitogenome) of Tryporyza incertulas (Walker) (Lepidoptera: Pyraloidea: Crambidae) was determined, characterized and compared with seven other species of superfamily Pyraloidea. The order of 37 genes was typical of insect mitochondrial DNA sequences described to date. Compared with other moths of Pyraloidea, the A + T biased (77.0%) of T. incertulas was the lowest. Eleven protein-coding genes (PCGs) utilized the standard ATN, but cox1 used CGA and nad4 used AAT as the initiation codons. Ten protein-coding genes had the common stop codon TAA, except nad3 having TAG as the stop codon, and cox2, nad4 using T, TA as the incomplete stop codons, respectively. All of the tRNA genes had typical cloverleaf secondary structures except trnS1(AGN), in which the dihydrouridine (DHU) arm did not form a stable stem-loop structure. There was a spacer between trnQ and nad2, which was common in Lepidoptera moths. A 6-bp motif ‘ATACTA’ between trnS2(UCN) and nad1, a 7-bp motif “AGC(T)CTTA” between trnW and trnC and a 6-bp motif “ATGATA” of overlapping region between atp8 and atp6 were found in Pyraloidea moths. The A + T-rich region contained an ‘ATAGT(A)’-like motif followed by a poly-T stretch. In addition, two potential stem-loop structures, a duplicated 19-bp repeat element, and two microsatellites ‘(TA)₁₂’ and ‘(TA)₉’ were observed in the A + T-rich region of T. incertulas mitogenome. Finally, the phylogenetic relationships of Pyraloidea species were constructed based on amino acid sequences of 13 PCGs of mitogenomes using Bayesian inference (BI) and maximum likelihood (ML) methods. These molecular-based phylogenies supported the morphological classification on relationships within Pyraloidea species.     

Fourteen complete mitochondrial genomes of butterflies from the genus Lethe (Lepidoptera,Nymphalidae, Satyrinae) with mitogenome-based phylogenetic analysis

The mitochondrial genome (mitogenome) can help us understand the phylogenetic relationships within the genus Lethe and the subfamily Satyrinae. In this study, we sequenced the complete mitogenomes of 14 Lethe species, which range in size from 15,225 to 15,271 bp, with both 37 genes (13 PCGs, 22 tRNAs, 2 rRNAs) and a noncoding A + T-rich region. The gene arrangement and orientation is similar to typical mitogenomes of Lepidoptera. The Ka/Ks ratio shows that cox1 has the slowest evolutionary rate. The secondary structure of trnN lacks the Pseudouracil loop (TψC loop) in most Lethe species. The inferred phylogenetic analyses show that Lethe is a well-supported monophyletic group, and reveal 2 major clades within the genus Lethe, which is consistent with previous morphological classifications.     

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.     

Complete mitochondrial genome of the stonefly Cryptoperla stilifera Sivec (Plecoptera: Peltoperlidae) and the phylogeny of Polyneopteran insects

We present the complete mitogenome of a stonefly, Cryptoperla stilifera Sivec (Plecoptera; Peltoperlidae). The mitogenome was a circular molecule consisting of 15,633 nucleotides, 37 genes and a A + T-rich region. C. stilifera mitogenome was similar to Pteronarcys princeps mitogenome (Plecoptera; Pteronarcyidae). All transfer RNA genes (tRNAs) had typical cloverleaf secondary structures except for trnSer (AGN), where the stem-loop structure of the dihydrouridine (DHU) arm was missing. The A + T-rich region of C. stilifera had two stem-loops and each had two interlink. Three conserved sequence blocks (CSBs) were present in the A + T-rich regions of C. stilifera, Peltoperla tarteri and Peltoperla arcuata. Moreover, many polynucleotide stretches (Poly N, N = A, T and C) in the A + T-rich region of C. stilifera Phylogenetic relationships of Polyneopteran species were constructed based on the nucleotide sequences of 13 protein coding genes (PCGs). Both maximum likelihood (ML) and Bayesian inference (BI) analyses supported Grylloblattodea as the sister group to Plecoptera + Dermaptera and Embiidina and Phasmatodea as sister groups.     

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.