Insect mitochondrial genomes (mitogenomes) contain a conserved set of 37 genes for an extensive diversity of lineages. Previously reported dictyopteran mitogenomes share this conserved mitochondrial gene arrangement, although surprisingly little is known about the mitogenome of Mantodea. We sequenced eight mantodean mitogenomes including the first representatives of two families: Hymenopodidae and Liturgusidae. Only two of these genomes retain the typical insect gene arrangement. In three Liturgusidae species, the trnM genes have translocated. Four species of mantis (Creobroter gemmata, Mantis religiosa, Statilia sp., and Theopompa sp.-HN) have multiple identical tandem duplication of trnR, and Statilia sp. additionally includes five extra duplicate trnW. These extra trnR and trnW in Statilia sp. are erratically arranged and form another novel gene order. Interestingly, the extra trnW is converted from trnR by the process of point mutation at anticodon, which is the first case of tRNA reassignment for an insect. Furthermore, no significant differences were observed amongst mantodean mitogenomes with variable copies of tRNA according to comparative analysis of codon usage. Combined with phylogenetic analysis, the characteristics of tRNA only possess limited phylogenetic information in this research. Nevertheless, these features of gene rearrangement, duplication, and reassignment provide valuable information toward understanding mitogenome evolution in insects.

Leishmaniasis is a worldwide but neglected disease of humans and animal transmitted by sand flies, vectors that also transmit other important diseases. Mitochondrial genomes contain abundant information for population genetic and phylogenetic studies, important in disease management. However, the available mitochondrial sequences of these crucial vectors are limited, emphasizing the need for developing more mitochondrial genetic markers. The complete mitochondrial genome of Phlebotomus chinensis was amplified in eight fragments and sequenced using primer walking. The mitochondrial genome of Phlebotomus papatasi was reconstructed from whole-genome sequencing data available on Genbank. The phylogenetic relationship of 24 selected representatives of Diptera was deduced from codon positions 1 and 2 for 13 protein coding genes, using Bayesian inference (BI) and maximum likelihood (ML) methods. We provide the first Phlebotomus (P. chinensis and P. papatasi) mitochondrial genomes. Both genomes contain 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and an A + T-rich region. The gene order of Phlebotomus mitochondrial genomes is identical with the ancestral gene order of insect. Phylogenetic analyses demonstrated that Psychodidae and Tanyderidae are sister taxa. Potential markers for population genetic study of Phlebotomus species were also revealed. The generated mitochondrial genomes of P. chinensis and P. papatasi represent a useful resource for comparative genomic studies and provide valuable future markers for the population genetic study of these important Leishmania vectors. Our results also preliminary demonstrate the phylogenetic placement of Psychodidae based on their mitochondrial genomes.

Paragyrodactylus Gvosdev and Martechov, 1953, a viviparous genus of ectoparasite within the Gyrodactylidae, contains three nominal species all of which infect Asian river loaches. The group is suspected to be a basal lineage within Gyrodactylus Nordmann, 1832 sensu lato although this remains unclear. Further molecular study, beyond characterization of the standard Internal Transcribed Spacer region, is needed to clarify the evolutionary relationships within the family and the placement of this genus. The mitochondrial genome of Paragyrodactylus variegatus You, King, Ye and Cone, 2014 was amplified in six parts from a single worm, sequenced using primer walking, annotated and analyzed using bioinformatic tools. The mitochondrial genome of P. variegatus is 14,517 bp, containing 12 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a major non-coding region (NCR). The overall A + T content of the mitochondrial genome is 76.3%, which is higher than all reported mitochondrial genomes of monogeneans. All of the 22 tRNAs have the typical cloverleaf secondary structure, except tRNACys, tRNASer1 and tRNASer2 that lack the dihydrouridine (DHU) arm. There are six domains (domain III is absent) and three domains in the inferred secondary structures of the large ribosomal subunit (rrnL) and small ribosomal subunit (rrnS), respectively. The NCR includes six 40 bp tandem repeat units and has the double identical poly-T stretches, stem-loop structure and some surrounding structure elements. The gene order (tRNAGln, tRNAMet and NCR) differs in arrangement compared to the mitochondrial genomes reported from Gyrodactylus spp. The Duplication and Random Loss Model and Recombination Model together are the most plausible explanations for the variation in gene order. Both morphological characters and characteristics of the mitochondrial genome support Paragyrodactylus as a distinct genus from Gyrodactylus. Considering their specific distribution and known hosts, we believe that Paragyrodactylus is a relict freshwater lineage of viviparous monogenean isolated in the high plateaus of central Asia on closely related river loaches.

The complete mitochondrial genome sequence is determined for Paracymoriza prodigalis (Leech, 1889). The 15,326 bp circular molecule possesses a gene organization and order identical to other sequenced Pyraloidea mitochondrial genomes. All tRNAs have the typical clover-leaf structure except for tRNA^Ser(AGN), which lacks the dihydrouridine (DHU) arm. The A+T-rich region of 343 bp includes the features common to the Lepidoptera, including the ‘ATAGA’ followed by an 19-bp poly-T stretch, but the tandem repeat sequences often appearing in available insects are not found. Phylogenetic relationships of eight subfamilies of 14 Pyraloidea species were constructed based on 13 PCGs of mitochondrial genomes using Bayesian inference (BI) and maximum likelihood (ML) methods. These phylogenies of the subfamilies within Pyraloidea accord well with morphological phylogenetic analysis except for the position of Schoenobiinae.

In the present study, we report five complete and one nearly complete mitochondrial genomes of the Pyraloidea including the first representatives from the Pyralinae (Pyralidae) and Glaphyriinae (Crambidae). We also conduct a comparative analysis of mitogenomic features of this group. Our results show that Pyraloidea mitogenomes evolved under a common trend found in lepidopteran mitogenomes and share several typical genomic characters. The extra conserved blocks are identified in the Pyraloidea control region, and diverse missing codons formed another unique trait within Pyraloidea mitogenome. Furthermore, we reconstruct the mitogenomic phylogeny of Pyraloidea and confirm the phylogenetic position of Pyralinae and Glaphyriinae within the Pyraloidea using BI and ML method based on multiple mitochondrial datasets.