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.     

Characterization of the complete mitochondrial genome of Diaphania pyloalis (Lepidoptera: Pyralididae)

The complete mitochondrial genome (mitogenome) of Diaphania pyloalis (Lepidoptera: Pyralididae) was determined to be 15,298 bp and has the typical gene organization of mitogenomes from lepidopteran insects. It consists of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The A + T content of this mitogenome is 80.83% and the AT skew is slightly positive. All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which is initiated by CGA. Only the cox2 gene has an incomplete stop codon consisting of just a T. All the tRNA genes display a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome is 332 bp in length, including several common features found in lepidopteran mitogenomes. Phylogenetic analysis showed that the D. pyloalis is close to Pyralididae.     

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.     

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

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.     

The first mitochondrial genome for Brahmin moths (Brahmaeidae) and implications for the higher phylogeny of Bombycoidea

Bombycoidea comprises 10 families and 4723 species, and the phylogenetic relationships among families are still in debate. In this study, we have determined the complete mitochondrial genome (mitogenome) of Brahmaea porphyria. The 15,429-bp mitogenome contains a common set of 37 mitochondrial genes including 13 protein-coding genes, 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and an inferred control region, and shares the conserved gene rearrangement (trnM-trnI-trnQ) in most ditrysian mitogenomes. Moreover, we analysed the secondary structure for all the tRNA genes of B. porphyria and the preference of codon usage in the PCGs of B. porphyria. The putative 373-bp control region (CR) possesses three types of conserved elements, including ATAGA, Ploy-T stretch, and microsatellite-like elements. A phylogenetic analysis among available Bombycoidea mitogenomes using the concatenated 37 mitochondrial genes appears to support the hypothesis of (Sphingidae+Bombycidae)+Saturniidae and the relatively basal phylogenetic position of Brahmaeidae within Bombycoidea.     

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

Complete mitochondrial genome sequence of bighead croaker Collichthys niveatus (Perciformes, Sciaenidae): a mitogenomic perspective on the phylogenetic relationships of Pseudosciaeniae

The monophyly and phylogenetic relationships of Pseudosciaeniae have long been controversial. Here we describe the mitochondrial genome (mitogenome) sequence of Collichthys niveatus. It is a circular double-stranded DNA molecule of 16,450 base pairs (bp) in length with a standard set of 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), 13 protein-coding genes as well as a non-coding control region. The mitogenome of C. niveatus shared common features with those of other bony fishes in terms of gene arrangement, base composition, and tRNA structures. The C. niveatus mitogenome exhibited pronounced strand-specific asymmetry in nucleotide composition, which was also reflected in the codon usage of genes oriented in opposite directions. Contrary to the typical structure of the control region, the central conserved blocks (CSB-D, -E, and -F) could not be detected in C. niveatus mitogenome. Phylogenetic analysis based on whole mitogenome sequences provided strong support for the monophyly of Pseudosciaeniae, and sister-group relationships of C. niveatus + Collichthys lucidus and Larimichthys crocea + Larimichthys polyactis, which was consistent with the traditional taxonomy. Unexpected divergence was found in two C. niveatus mitogenomes and several hypotheses were proposed to explain this observation including misidentification and introgressive hybridization between C. niveatus and L. polyactis, and polyphyletic origin of C. niveatus. We considered species misidentification to be the main hypothesis. However, additional data is essential to test these proposed hypotheses.     

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 Biston panterinaria (Lepidoptera: Geometridae), with phylogenetic utility of mitochondrial genome in the Lepidoptera

The complete mitochondrial genome (mitogenome) of the Chinese pistacia looper Biston panterinaria was sequenced and annotated (15,517 bp). It contains the typical 37 genes of animal mitogenomes and a high A + T content (79.5%). All protein coding genes (PCGs) use standard ATN initiation codons except for cytochrome c oxidase 1 (COX1) with CGA. Eleven PCGs use a common stop codon of TAA or TAG, whereas COX2 and NADH dehydrogenase 4 (ND4) use a single T. All transfer RNA (tRNA) genes have the typical clover-leaf structure with the exception of tRNA^Ser(AGN). We reconstructed a preliminary mitochondrial phylogeny of six ditrysian superfamilies and performed comparative analyses of inference methods (Bayesian Inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP)), dataset compositions (including and excluding 3rd codon positions), and alignment methods (Muscle, Clustal W, and MAFFT). Our analyses indicated that inference methods and dataset compositions more significantly affected the phylogenetic results than alignment methods. BI analysis consistently revealed uncontroversial relationships with all dataset compositions. By contrast, ML analysis failed to reconstruct stable phylogeny at two nodes, whereas MP analysis had more difficulties in the tree resolution and nodal support. Distinct from most previous studies, our analyses revealed that Geometroidea had a closer lineage relationship with Bombycoidea than Noctuoidea. Similar to previous molecular studies, our analyses revealed that Hesperiidae were nested in the Papilionoidea clade, providing further evidence to the previous concept that Papilionoidea was paraphyletic, and none of the butterflies were associated with the Macroheterocera.     

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 mitogenome sequence of the Japanese oak silkmoth, <Emphasis Type="Italic">Antheraea yamamai</Emphasis> (Lepidoptera: Saturniidae)

The 15,338-bp long complete mitochondrial genome (mitogenome) of the Japanese oak silkmoth, Antheraea yamamai (Lepidoptera: Saturniidae) was determined. This genome has a gene arrangement identical to those of all other sequenced lepidopteran insects, but differs from the most common type, as the result of the movement of tRNA^Met to a position 5′-upstream of tRNA^Ile. No typical start codon of the A. yamamai COI gene is available. Instead, a tetranucleotide, TTAG, which is found at the beginning context of all sequenced lepidopteran insects was tentatively designated as the start codon for A. yamamai COI gene. Three of the 13 protein-coding genes (PCGs) harbor the incomplete termination codon, T or TA. All tRNAs formed stable stem-and-loop structures, with the exception of tRNA^Ser(AGN), the DHU arm of which formed a simple loop as has been observed in many other metazoan mt tRNA^Ser(AGN). The 334-bp long A + T-rich region is noteworthy in that it harbors tRNA-like structures, as has also been seen in the A + T-rich regions of other insect mitogenomes. Phylogenetic analyses of the available species of Bombycoidea, Pyraloidea, and Tortricidea bolstered the current morphology-based hypothesis that Bombycoidea and Pyraloidea are monophyletic (Obtectomera). As has been previously suggested, Bombycidae (Bombyx mori and B. mandarina) and Saturniidae (A. yamamai and Caligula boisduvalii) formed a reciprocal monophyletic group.     

The complete mitochondrial genome of the sycamore lace bug Corythucha ciliata (Hemiptera: Tingidae)

The complete mitochondrial genome of the sycamore lace bug, Corythucha ciliata, was sequenced in this study. It represents the first sequenced mitogenome of family Tingidae in Heteroptera. The mitogenome of C. ciliata is 15,257 bp and contains 37 genes including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and a large non-coding region. Gene arrangement, nucleotide content, codon usage, and amino acid composition and asymmetry indicate a high degree of conservation with six other species of Cimicomorpha. The 13 PCGs initiated with ATN as the start codon and terminated with TAA, TA or T as stop codon. The evolutionary rate of each PCG was different, among which ATP8 showed the highest rate while ATP6 indicated the lowest rate. The 22 tRNAs genes apparently fold into a typical cloverleaf structure; however, the anticodon (TTC) of trnSer (AGN) differs from other Heteropteran insects. Secondary structure modeling of rRNA genes revealed similarity to other insects, except for two incomplete helices (H1648 and H2735) in lrRNA. The predicted secondary structure of lrRNA indicates 45 helices in six domains, whereas srRNA has 27 helices in three domains. Three potential stem–loops and two tandem repeats (–TCTAAT–) were identified in the A+T-rich region. Phylogenetic analysis indicated that C. ciliata is a sister group to other Heteroptera species based on analysis of the 13 PCGs.     

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.     

Characterization of the complete mitochondrial genome of Cryptotermes domesticus (Blattodea: Kalotermitidae): Genome description and phylogenetic implications

The complete mitochondrial genome of Cryptotermes domesticus (Haviland) was sequenced and annotated to study its characteristics and the phylogenetic relationship of C. domesticus to other termite species. The mitogenome of C. domesticus is a circular, close, and double-stranded molecule with a length of 15,655 bp. The sequenced mitogenome contains 37 typical genes, which are highly conserved in gene size, organization, and codon usage. Transfer RNA genes (tRNAs) also have typical secondary structures. All of the 13 protein-coding genes (PCGs) start with an ATN codon, except for nad4, which starts with GTG and terminates with the terminal codon TAA and TAG or the incomplete form T-- (cox2 and nad5). Most tRNAs have a typical cloverleaf structure, except for trnS1, in which this form is replaced by a simple loop and lacks the dihydrouridine (DHU) arm. The nucleotide diversity (Pi) and nonsynonymous (Ka)/synonymous (Ks) mutation rate ratios indicate that nad1, cox1, and cox3 are the most conserved genes, and that cox1 has the lowest rate of evolution. In addition, an 89 bp repeated sequence was found in the A + T-rich region. Phylogenetic analysis was performed using Bayesian inference (BI) and maximum likelihood (ML) methods based on 13 PCGs, and the monophyly of Kalotermitidae was supported.