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 Teinopalpus aureus guangxiensis (Lepidoptera: Papilionidae) and related phylogenetic analyses

In this study, we sequenced the complete mitochondrial genome of Teinopalpus aureus guangxiensis (Lepidoptera: Papilionidae), which is considered as an endemic species in China. It is listed as a vulnerable species by International Union for Conservation of Nature and Natural Resources Red List and also a first class endangered species in China. The complete mtDNA from T. aureus guangxiensis was 15,235 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The T. aureus guangxiensis genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. All PCGs of T. aureus guangxiensis mitogenome start with a typical ATN codon and terminate in the common stop codon TAA, except that ND1 gene uses TTA, ND3 gene uses ATT, and ND4 and ND4L gene use TAA. The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, and phylogenetic results confirmed that Nymphalidae, Lycaenidae, Papilionidae, Pieridae are monophyletic clades.     

The complete mitochondrial genome of Calappa bilineata: The first representative from the family Calappidae and its phylogenetic position within Brachyura

In this study, we determined the complete mitogenome sequence of Calappa bilineata, which is the first mitogenome of Calappidae up to now. The total length is 15,606 bp and includes 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs and one control region. The genome composition is highly A + T biased (68.7%), and exhibits a negative AT-skew (−0.010) and GC-skew (−0.267). As with other invertebrate mitogenomes, the PCGs start with the standard ATN and stop with the standard TAN codons or incomplete T. Phylogenetic analysis showed that C. bilineata was most closely related to Matuta planipes (Matutidae), and these two species formed a sister clade, constituting a Calappoidea group and forming a sister clade with part of Eriphioidea. The existence of the polyphyletic families raised doubts over the traditional classification system. These results will help to better understand the features of the C. bilineata mitogenome and lay foundation for further evolutionary relationships within Brachyura.     

The complete mitochondrial genome of <Emphasis Type="Italic">Damora sagana</Emphasis> and phylogenetic analyses of the family Nymphalidae

The monotypic genus Damora (Nymphalidae, Heliconiinae) contains a single species, Damora sagana, which is widely distributed across southern China. Herein, its complete mitogenome was sequenced to further understand lepidopteran mitogenome characteristics, reconstruct the nymphalid family phylogeny, and infer the subdivision of Heliconiinae species. The circular mitogenome was 15,151 bp long, abundant in A and T, and comprised of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and one control region with a gene arrangement typical of lepidopteran mitogenomes. ATN codons initiated all PCGs, except cytochrome c oxidase subunit 1 (COX1), which was initiated by a CGA sequence as has been observed in other lepidopterans. Three PCGs (COX1, COX2 and ND4) employed a single T termination signal, whereas others had the typical complete termination codon (TAA). All tRNA genes were folded into the typical cloverleaf structure except for tRNA-Ser (AGN). The A+T-rich region included the conserved motif ‘ATAGA’ followed by a 17 bp poly-T stretch, which was also observed in tribe Argynnini mitogenomes. A phylogenetic tree was constructed via multiple methods using the 13 PCGs data of D. sagana and other available mitogenomes of nymphalid species. All three phylogenetic trees yielded the same topology. These results were consistent with those from previous studies of most major nymphalid groups, except those regarding tribe subdivision in certain subfamilies such as Argynnini + (Acraeini + Heliconiini) for Heliconiine. Furthermore, our analyses identified that the genus Cethosia was grouped with the genus Acraea composing the tribe Acraeini with strong support.     

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.     

The Complete Mitochondrial Genome of the Beet Webworm,Spoladea recurvalis (Lepidoptera: Crambidae) and Its Phylogenetic Implications

The complete mitochondrial genome (mitogenome) of the beet webworm, Spoladea recurvalis has been sequenced. The circular genome is 15,273 bp in size, encoding 13 protein-coding genes (PCGs), two rRNA genes, and 22 tRNA genes and containing a control region with gene order and orientation identical to that of other ditrysian lepidopteran mitogenomes. The nucleotide composition of the mitogenome shows a high A+T content of 80.9%, and the AT skewness is slightly negative (-0.023). All PCGs start with the typical ATN codons, except for COX1, which may start with the CGA codon. Nine of 13 PCGs have the common stop codon TAA; however, COX1, COX2 and ND5 utilize the T nucleotide and ND4 utilizes TA nucleotides as incomplete termination codons. All tRNAs genes are folded into the typical cloverleaf structure of mitochondrial tRNAs, except for the tRNA^Ser(AGY) gene, in which the DHU arm fails to form a stable stem-loop structure. A total of 157 bp intergenic spacers are scattered in 17 regions. The overlapping sequences are 42 bp in total and found in eight different locations. The 329 bp AT-rich region is comprised of non-repetitive sequences, including the motif ATAG, which is followed by a 14 bp poly-T stretch, a (AT₁₁ microsatellite-like repeat, which is adjacent to the motif ATTTA, and a 9 bp poly-A, which is immediately upstream from the tRNA^Met gene. Phylogenetic analysis, based on 13 PCGs and 13 PCGs+2 rRNAs using Bayesian inference and Maximum likelihood methods, show that the classification position of Pyraloidea is inconsistent with the traditional classification. Hesperioidea is placed within the Papilionoidea rather than as a sister group to it. The Pyraloidea is placed within the Macrolepidoptera with other superfamilies instead of the Papilionoidea.     

Characterization of the complete mitogenome of Trachylophus sinensis (Coleoptera: Cerambycidae: Cerambycinae), the type species of Trachylophus and its phylogenetic implications

Complete mitochondrial genomes (mitogenomes) have long been proved as reliable markers for phylogenetic reconstruction among diverse animal groups, especially benefited from recent rapid development of sequencing techniques. However, the mitogenomes of many important clades remain poorly represented, which restricted the understanding of macroscale evolutionary history of these groups. Here, we sequenced and characterized the complete mitogenome of Trachylophus sinensis, a type species of the Trachylophus genus, which also represents the first sequenced mitogenome in this genus. The complete circular mitogenome was 15,746 bp in length, containing 37 typical genes and one noncoding AT-rich control region. The nucleotide composition of the mitogenome was highly A + T biased, accounting for 70.07 % of the whole mitogenome with a slightly positive AT skewness (0.106). The 13 Protein coding genes (PCGs) used ATN as their start codons, except nad1 which used TTG. All tRNA genes were predicted with a characteristic cloverleaf secondary structure except trnS1^(AGN), whose dihydrouridine (DHU) arm was replaced by a simple loop. Phylogenetic analyses recovered Cerambycinae as a monophyletic group with high node supports and the sister relationship between T. sinensis and Nadezhdiella cantori. However, we found that deeper nodes showed not strong support, which may be caused by limited taxa sampling in our study. More mitogenomes should be sequenced representing various taxonomic levels, especially closely related species, which will enhance our understanding of phylogenetic relationships among Cerambycinae.     

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.     

Complete mitochondrial genome of the black‐tailed hornet,Vespa ducalis (Hymenoptera: Vespidae): Genomic comparisons in Vespoidea

We sequenced the complete mitochondrial genome (mitogenome) of the black‐tailed hornet, Vespa ducalis (Hymenoptera: Vespidae). The genome was 15,779‐bp long and contained typical sets of genes [13 protein‐coding genes (PCGs), 22 tRNAs, and 2 rRNAs]. The V. ducalis A + T‐rich region was 166‐bp long and was the shortest of all sequenced Vespoidea genomes, including Vespa. The genome was highly biased toward A/T nucleotides—80.1 % in the whole genome, 77.8 % in PCGs, 83.4–85.6 % in RNAs, and 92.8 % in the A + T‐rich region. These values are well within the typical range for genes and regions of Vespoidea mitogenomes. Start and stop codons in several Vespa species—including V. ducalis—were diversified, despite these species belonging to the same genus. In comparison with the ancestral mitogenomes, Vespa mitogenomes—including that of V. ducalis—showed substantial gene rearrangement; however, we detected no gene rearrangement among Vespa species. We conducted phylogenetic reconstruction based on concatenated sequences of 13 PCGs and two rRNAs (12,755 bp ) in available species of Vespoidea—21 species in six subfamilies in two families (Vespidae and Formicidae). The Bayesian inference and maximum likelihood (ML) methods revealed that each family formed strong monophyletic groups [Bayesian posterior probability (BPP) = 1; ML, 100 %]. Moreover, V. ducalis and V. mandarinia formed a strong sister group (BPP = 1; ML, 94 %).     

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.     

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.     

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.     

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.     

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

Mitochondrial genome of <Emphasis Type="Italic">Pteronotus personatus</Emphasis> (Chiroptera: Mormoopidae): comparison with selected bats and phylogenetic considerations

We described the complete mitochondrial genome (mitogenome) of the Wagner’s mustached bat, Pteronotus personatus, a species belonging to the family Mormoopidae, and compared it with other published mitogenomes of bats (Chiroptera). The mitogenome of P. personatus was 16,570 bp long and contained a typically conserved structure including 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one control region (D-loop). Most of the genes were encoded on the H-strand, except for eight tRNA and the ND6 genes. The order of protein-coding and rRNA genes was highly conserved in all mitogenomes. All protein-coding genes started with an ATG codon, except for ND2, ND3, and ND5, which initiated with ATA, and terminated with the typical stop codon TAA/TAG or the codon AGA. Phylogenetic trees constructed using Maximum Parsimony, Maximum Likelihood, and Bayesian inference methods showed an identical topology and indicated the monophyly of different families of bats (Mormoopidae, Phyllostomidae, Vespertilionidae, Rhinolophidae, and Pteropopidae) and the existence of two major clades corresponding to the suborders Yangochiroptera and Yinpterochiroptera. The mitogenome sequence provided here will be useful for further phylogenetic analyses and population genetic studies in mormoopid bats.     

Mitogenome of Knodus borki (Cypriniformes: Characidae): genomic characterization and phylogenetic analysis

Background

The taxonomic status of Knodu in the family Characidae is not yet clear. This study aimed to address this by sequencing and annotating Knodu borki Zarske, 2008.

Materials and results

K. borki Zarske, 2008 was sequenced using a Hiseq platform and the complete mitogenome was assembled in SPAdes v3.15.2 and SOAPdenovo2 v.2.01. The mitogenome of K. borki from Guangzhou, the first sequenced species of the genus Knodu, is 16,837 bp in length and contains 13 protein-coding genes (PCGs), two ribosomal (r) RNAs, 22 transfer (t) RNAs, and one D-loop. Among these 37 genes, 28 are encoded by the heavy strand, while nine are encoded by the light strand. Twenty-one of the tRNAs can form typical cloverleaf secondary structures, except tRNA-Ser1, which lacks dihydrouridine arms. All PCGs have the same start codon (ATG), with the exception of COI (GTG). Four PCGs (ND1, ATP8, ND4L, and ND5) have TAA as the stop codon, ND6 has TAG as the stop codon, COI has AGG as the stop codon, and the remaining seven genes have incomplete stop codons of TA-/T-(ND2, COII, COIII, ND3, ND4, and Cyt b as T-, ATP6 as TA-). Phylogenetic analysis showed that K. borki belongs to the family Characidae.

Conclusions

Our findings demonstrate that K. borki belongs to the family Characidae, due to consistency with the morphological identification. This study provides molecular information for further research on the phylogeny of the genus Knodus and for analyses of the taxonomic status of Characidae.

     

Complete mitochondrial genome of <Emphasis Type="Italic">Ampittia dioscorides</Emphasis> (Lepidoptera: Hesperiidae) and its phylogenetic analysis

The complete mitochondrial genome of Ampittia dioscorides (Lepidoptera: Hesperiidae) was determined. The sequenced genome is a circular molecule of 15313 bp, containing 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an A + T-rich region. The gene arrangements and transcribing directions are identical to those in most of the reported lepidopteran mitogenomes. The base composition of the whole genome and genes or regions are also similar to those in other lepidopteran species. All the PCGs are initiated by typical ATN codons; the exception being COI, which begins with a CGA codon. Eight genes (ND2, ATPase8, ATPase6, COIII, ND5, ND4L, ND6, and Cytb) end with a TAA stop codon, and two genes (ND1 and ND3) end with TAG. The remaining three genes (COI and COII, which end with TA-, and ND4, which ends with T-) have incomplete stop codons. All tRNAs have the typical clover-leaf structure of mitochondrial tRNAs, with the exception of tRNA^Ser(AGY). On the basis of the concatenated nucleotide and amino acid sequences of the 13 PCGs and wingless gene of 22 butterfly species, maximum parsimony (MP) and Bayesian inference (BI) trees were constructed, respectively. Both MP and BI trees had the same topological structure: ((((Nymphalidae + Danaidae) + Lycaenidae) + Pieridae) + Papilionidae) + Hesperiidae). The results provide support for Hesperiidae as a superfamily-level taxon.     

Characterization of the complete mitochondrial genome of Drawida gisti (Metagynophora,Moniligastridae) and comparison with other Metagynophora species

Here, the complete mitochondrial genome (mitogenome) of Drawida gisti was sequenced and compared with the mitogenomes of other Metagynophora species. The circular mitogenome was 14,648 bp in length and contained two ribosomal RNA genes (rRNAs), 13 protein-coding genes (PCGs), and 22 transfer RNA genes (tRNAs). The types of constitutive genes and the direction of the coding strand that appeared in Drawida mitogenome were identical to those observed in other Metagynophora species, except for a missing lengthy non-coding region. The conservative relationships between Drawida species were supported by the overall analyses of 13 PCGs, two rRNAs, and 22 tRNAs. A comparison of the Metagynophora mitogenomes revealed that the ATP8 gene possessed the highest polymorphism among the 13 PCGs and two rRNAs. Phylogenetic analysis suggested that the Moniligastridae contained Drawida, which is a primitive Metagynophora group. Our study provides a step forward toward elucidating the evolutionary linkages within Drawida and even Metagynophora.