Complete mitochondrial genomes of three skippers in the tribe Aeromachini (Lepidoptera: Hesperiidae: Hesperiinae) and their phylogenetic implications

The mitochondrial genome is now widely used in the study of phylogenetics and molecular evolution due to its maternal inheritance, fast evolutionary rate, and highly conserved gene content. To explore the phylogenetic relationships of the tribe Aeromachini within the subfamily Hesperiinae at the mitochondrial genomic level, we sequenced and annotated the complete mitogenomes of 3 skippers: Ampittia virgata, Halpe nephele, and Onryza maga (new mitogenomes for 2 genera) with a total length of 15,333 bp, 15,291 bp, and 15,381 bp, respectively. The mitogenomes all contain 13 protein‐coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a noncoding A + T‐rich region and are consistent with other lepidopterans in gene order and type. In addition, we reconstructed the phylogenetic trees of Hesperiinae using maximum likelihood (ML) and Bayesian inference (BI) methods based on mitogenomic data. Results show that the tribe Aeromachini in this study robustly constitute a monophyletic group in the subfamily Hesperiinae, with the relationships Coeliadinae + (Euschemoninae + (Pyrginae + ((Eudaminae + Tagiadinae) + (Heteropterinae + ((Trapezitinae + Barcinae) + Hesperiinae))))). Moreover, our study supports the view that Apostictopterus fuliginosus and Barca bicolor should be placed out of the subfamily Hesperiinae.     

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.     

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.     

Mitochondrial genomes of the green macroalga Ulva pertusa (Ulvophyceae,Chlorophyta): novel insights into the evolution of mitogenomes in the Ulvophyceae

To further understand the trends in the evolution of mitochondrial genomes (mitogenomes or mtDNAs) in the Ulvophyceae, the mitogenomes of two separate thalli of Ulva pertusa were sequenced. Two U. pertusa mitogenomes (Up1 and Up2) were 69,333 bp and 64,602 bp in length. These mitogenomes shared two ribosomal RNAs (rRNAs), 28 transfer RNAs (tRNAs), 29 protein‐coding genes, and 12 open reading frames. The 4.7 kb difference in size was attributed to variation in intron content and tandem repeat regions. A total of six introns were present in the smaller U. pertusa mtDNA (Up2), while the larger mtDNA (Up1) had eight. The larger mtDNA had two additional group II introns in two genes (cox1 and cox2) and tandem duplication mutations in noncoding regions. Our results showed the first case of intraspecific variation in chlorophytan mitogenomes and provided further genomic data for the undersampled Ulvophyceae.     

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

Evolutionary relatedness of mackerels of the genus Scomber based on complete mitochondrial genomes: Strong support to the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as distinct species

Mackerels of the genus Scomber are commercially important species, but their taxonomic status is still controversial. Although previous phylogenetic data support the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as separate species, it is only based on the analysis of partial mitochondrial and nuclear DNA sequences. In an attempt to shed light on this relevant issue, we have determined the complete mitochondrial DNA sequence of S. colias, S. japonicus, and Scomber australasicus. The total length of the mitogenomes was 16,568 bp for S. colias and 16,570 bp for both S. japonicus and S. australasicus. All mitogenomes had a gene content (13 protein-coding, 2 rRNAs, and 22 tRNAs) and organization similar to that observed in Scomber scombrus and most other vertebrates. The major noncoding region (control region) ranged between 865 and 866 bp in length and showed the typical conserved blocks. Phylogenetic analyses revealed a monophyletic origin of Scomber species with regard to other scombrid fish. The major finding of this study is that S. colias and S. japonicus were significantly grouped in distinct lineages within Scomber cluster, which phylogenetically constitutes evidence that they may be considered as separate species. Additionally, molecular data here presented provide a useful tool for evolutionary as well as population genetic studies.     

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 and phylogenetic analysis of the complete mitogenome of Allactaga sibirica (Rodentia: Dipodidae)

Allactaga sibirica (Dipodidae) is widely distributed in the northwestern arid regions of China. The complete mitochondrial genome (mitogenome) of A. sibirica was 16,685 bp in length; included 13 protein-coding genes, 2 ribosome genes, 22 tRNA genes, and one control region; and had a structure that was typical of vertebrates. The base composition and codon usage of the mitogenome are described, and the structure of the non-coding sequence in the A. sibirica is reported for the first time. The putative origin of replication for the light strand of A. sibirica was approximately 45 bp long, and was highly conserved in the stem-loop and adjacent sequences. Phylogenetic analyses showed high resolution in each of the main divergent clades within Dipodoidea using mitogenomes data. The results indicated that the Zapodidae group was a representative of very basal taxon in Dipodoidea, and shared a common ancestor with Dipodidae species. Within Dipodidae clade, Allactaginae species was at basal position, and this result was in line with previous molecular systematic and morphological studies. Furthermore, Euchoreutes naso and A. sibirica had a close relationship could implicate a sister-group relationship between Euchoreutinae and Allactaginae. Meanwhile, this work also provided a set of useful data on phylogeny and molecular evolution in Dipodidae species.     

Two novel cricetine mitogenomes: Insight into the mitogenomic characteristics and phylogeny in Cricetinae (Rodentia: Cricetidae)

Both Cricetus cricetus and Phodopus sungorus mitochondrial genomes (mitogenomes) were sequenced and elaborated for the first time in the present study. Their mitogenomes contained 37 genes and showed typical characteristics of the vertebrate mitogenome. Comparative analysis of 10 cricetine mitogenomes indicated that they shared similar characteristics with those of other cricetines in terms of genes arrangement, nucleotide composition, codon usage, tRNA structure, nucleotide skew and the origin of replication of light strand. Phylogenetic relationship of the subfamily Cricetinae was reconstructed using mitogenomes data with the methods of Bayesian Inference and Maximum Likelihood. Phylogenetic analysis indicated that Cricetulus kamensis was at basal position and phylogenetically distant from all other Cricetulus species but had a close relationship with the group of Phodopus, and supported that the genus Urocricetus deserved as a separate genus rank. The phylogenetic status of Tscherskia triton represented a separate clade corresponding to a diversified cricetine lineage (Cricetulus, Allocricetulus, and Cricetus).     

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

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 sequences of <Emphasis Type="Italic">Chelodina rugosa</Emphasis> and <Emphasis Type="Italic">Chelus fimbriata</Emphasis> (Pleurodira: Chelidae): implications of a common absence of initiation sites (O<Subscript>L</Subscript>) in pleurodiran turtles

Within the order Testudines, while phylogenetic analyses have been performed on the suborder Cryptodira with complete mitochondrial genomes (mitogenomes), mitogenomic information from another important suborder Pleurodira has been inadequate. In the present study, complete mitochondrial DNA (mtDNA) sequences of two chelid turtles Chelodina rugosa and Chelus fimbriata were firstly determined, the lengths of which were 16,582 and 16,661 bp respectively. As the typical vertebrate mitogenome, both mtDNAs consist of 13 protein coding genes, 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and a long noncoding region (control region, CR). However, the initiation sites for light-strand replication (O_L), which has been identified in all reported Cryptodire mitogenomes, were not found in the putative position of the two chelid turtles and African helmeted turtle Pelomedusa subrufa. The results suggested that the absence of mitogenomic initiation sites (O_L) could be a characteristic of Pleurodira. Phylogenetic relationships of chelid turtles and other turtles were reconstructed using the reported mitogenomes. Both maximum likelihood (ML) and Bayesian inference (BI) analyses suggested the monophyly of Pleurodira and Cryptodira as well as a sister group relationship between the two chelid turtles with strong statistical support. This phylogenetic framework was also utilized to estimate divergence dates among lineages using relaxed-clock methods combined with fossil evidence. Divergence estimates revealed that genus Chelodina diverged from genus Chelus in Late Cretaceous (~83 million years ago (mya)), and the time is consistent with the vicariance of the fragments which was caused by Gondwana split.     

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.     

Characterization of the complete mitochondrial genome of Japanagallia spinosa and Durgades nigropicta (Hemiptera: Cicadellidae: Megophthalminae)

In this study, we determined and analyzed the complete mitochondrial genomes (mitogenomes) of Japanagallia spinosa and Durgades nigropicta (Hemiptera: Megophthalminae). The circular genome were 15,655 bp long in J. spinosa (GenBank: KY123686) and 15,974 bp long in D. nigropicta (GenBank: KY123687). The J. spinosa and D. nigropicta mitogenomes both contained 37 genes and the gene order was similar to that in other leafhoppers. All of the protein-coding genes started with ATN. In the J. spinosa mitogenome, the nad3, nad4L, and cytb genes used TAG as a stop codon, the atp8 and nad1 genes used TGA, and the cox2 gene used a single T. However, in the D. nigropicta mitogenome, three genes used a single T as the stop codon, whereas the nad3 gene used TAG. We predicted the secondary structures of the rRNAs in J. spinosa and D. nigropicta. The secondary structure of rrnL comprised six domains (domain III is absent in arthropods) with 42 helices and that of rrnS comprised three structural domains with 26 helices. Comparisons of J. spinosa and D. nigropicta detected some differences in H577 and H673. We determined the structural organization of the control regions in the mitogenomes of leafhoppers, where three types of repeat regions were found in most. The phylogenetic relationships between J. spinosa and D. nigropicta with related lineages were reconstructed using Bayesian inference and maximum likelihood analyses. The monophyly of each superfamily considered in this study was confirmed by the clades in the phylogenetic tree. And in this study, Cicadellidae was resolved as monophyletic by the phylogenetic analysis. This mitogenome information for J. spinosa and D. nigropicta could facilitate future studies of mitogenomic diversity and the evolution of related insect lineages.     

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.     

LINE-1 distribution in six rodent genomes follow a species-specific pattern

L1 distribution in mammal’s genomes is yet a huge riddle. However, these repetitive sequences were already found in all chromosomic regions, and in general, they seem to be nonrandomly distributed in the genome. It also seems that after insertion and when they are not deleterious, they are always involved in dynamic processes occurring on that particular chromosomic region. Furthermore, it seems that large-scale genome rearrangements and L1 activity and accumulation are somehow interconnected. In the present study, we analysed L1 genomic distribution in Tateragambiana (Muridae, Gerbillinae), Acomys sp. (Muridae, Deomyinae), Cricetomys sp. (Nesomyidae, Cricetomyinae), Microtusarvalis (Cricetidae, Arvicolinae), Phodopus roborovskii and P. sungorus (Cricetidae, Cricetinae). All the species studied here seems to exhibit a species-specific pattern. Possible mechanisms, and processes involved in L1 distribution and preferential accumulation in certain regions are discussed.     

The mitogenomes of three beetles (Coleoptera: Polyphaga: Cucujiformia): New gene rearrangement and phylogeny

The mitochondrial genome (mitogenome) has been extensively used for studying phylogenetic relationships at different taxonomic levels. Several molecular analyses have been performed, but the phylogenetic relationships among infraorders in Polyphaga have not been well resolved. In this work, three nearly complete mitogenomes of Coleoptera, Sitophilus oryzae, Oryzaephilus surinamensis and Callosobruchus chinensis, were determined. The O. surinamensis and S. oryzae mitogenomes harbor gene content typical of other Polyphaga mitogenomes, while a gene rearrangement (trnQ) was found in the C. chinensis mitogenome. The mitogenomes of these three Coleoptera species each consist of approximately 13 protein-coding genes, 22 tRNA genes, two rRNA genes and one A + T-rich region. Phylogenetic analysis within Polyphaga was carried out based on mitochondrial data. The phylogenetic results within Polyphaga support the basal position of Cyphon sp., which belonged to Scirtoidea, Elateriformia. Within Cucujiformia, monophyletic Curculionoidea, Chrysomeloidea and Tenebrionoidea were confirmed.     

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.     

Two novel mitogenomes of Dipodidae species and phylogeny of Rodentia inferred from the complete mitogenomes

Two novel mitogenomes of Eozapus setchuanus (KJ648495) and Sicista concolor (KJ648496) were reported and their total lengths were 16,630 bp and 16,493 bp, respectively. Both mitogenomes which were analogous to other rodent mitogenomes, contained 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a control region. Specifically, the ND2 gene of S. concolor has three amino acids lesser than that of two other Dipodidae species (E. setchuanus and Jaculus jaculus) due to a premature termination codon in the 3′ end. We detected a tandem repeat cluster of 221 bp and 274 bp in the control region of S. concolor and E. setchuanus, respectively. Along with phylogenetic relationship analysis, we speculated that the tandem repeats in control regions might be common in Dipodinae species. Our phylogenetic analysis using concatenated mitochondrial gene datasets suggested five suborder and 16 family monophyletic groups in 54 rodent taxa sampled and strongly supported a basal position of the squirrel-related clade (PP = 1; BP = 100). Dipodidae had a sister-group relationship with Muroidea, and Sicistinae was in the base of Dipodidae clade. The complete mitochondrial genomes showed high resolution in deep-level phylogenetic relationship reconstructions of Rodentia.     

The complete mitochondrial genome of Tuta absoluta (Lepidoptera: Gelechiidae) and genetic variation in two newly invaded populations in China

Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae) is a devastating invasive pest worldwide, causing severe damage to tomatoes. Recently, it has been recorded in the northwestern and southwestern parts of China. Here, the mitogenomes and genetic variation of two newly invaded T. absoluta populations in Xinjiang and Yunnan, were determined. The results showed that the complete mitogenome size of T. absoluta is 15298 bp for the individual from Xinjiang and 15296 bp for the individual from Yunnan, which were both longer than the reported mitogenome from Spain (15290 bp). The mitogenome sequences of individuals collected from three locations showed high levels of sequence similarity, except for 8 polymorphic sites, which were in genes cox2 (1 site), cox3 (2 sites), cob (1 site), atp6 (1 site), nad1 (2 sites) and nad5 (1 site). Tuta absoluta mitogenomes share many features with other 6 Gelechiidae mitogenomes, except for several differences in the start and stop codons of protein-coding genes and the length of intergenic spacers. Seven partial mitochondrial genes (cox1, cox2, cox3, atp6, cob, nad1, and nad5) were used for genetic variation analysis, and significant population differentiation was found between the two populations based on cox2, atp6, nad1, and nad5. The complete mitogenomes and sensitive mitochondrial gene markers reported here provide useful data for further population genetics study of this pest.