Evolution and structural conservation of the control region of insect mitochondrial DNA

A Comparative Analysis of Mitogenomes in Species of the Tapinoma nigerrimum Complex and Other Species of the Genus Tapinoma (Formicidae, Dolichoderinae)

Using next-generation sequencing data, the complete mitogenomes of six species from the genus Tapinoma were assembled. This study explores the mitochondrial genomes of Tapinoma species, among them the five species from the Tapinoma nigerrimum complex, comparing them with each other and with other species from Dolichoderinae subfamily to understand their evolutionary relationships and evolution. Tapinoma mitochondrial genomes contain the typical set of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and the A + T-rich control region. A phylogenetic analysis using the protein-coding gene sequences from available Dolichoderinae mitogenomes supports the monophyletic nature of the genus Tapinoma, with the T. nigerrimum complex forming a well-supported clade. Key findings include genetic traits unique to the T. nigerrimum complex, such as a start codon in the atp8 gene and a complete stop codon in cox1, distinguishing them from other Tapinoma species. Additionally, a gene rearrangement involving tRNA-Trp, tRNA-Cys, and tRNA-Tyr was found exclusively in the Tapinoma species, suggesting a potential phylogenetic marker for the genus.

Characterization, Codon Usage Pattern and Phylogenetic Implications of the Waterlily Aphid Rhopalosiphum nymphaeae (Hemiptera: Aphididae) Mitochondrial Genome

The water lily aphid, Rhopalosiphum nymphaeae, is the only known aphid that can live in both terrestrial and aquatic conditions. In this study, the complete mitochondrial genome of R. nymphaeae was generated using Illumina sequencing technology. The typical circular DNA mitochondrial genome of R. nymphaeae is 15,772 bp in length, with a high A+T content (84.34%). It contains 37 coding genes (13 protein-coding genes, 22 transport RNAs, and two ribosomal RNAs) and two non-coding regions (one control region and one repeat region). Enc-plot, PR2-bias, and neutrality plot analysis indicated that the codon usage of the protein-coding genes is mainly affected by natural selection. The evolution rate analysis (the ratio of nonsynonymous to synonymous, Ka/Ks) indicated that all the PCGs in R. nymphaeae are under a strong purifying selection. The control region has conserved structure elements, and two types of tandem repeat units exist. The length and sequence of the aphid-unique repeat region has high similarity with closely related species. Phylogenetic analyses determined by both maximum likelihood and Bayesian inference support the monophyly of Aphidinae, Aphidini, Aphidina, and Rhopalosiphina. However, the monophyly of the genera in Rhopalosiphina, such as Rhopalosiphum, is still not resolved. This study may help us to understand the phylogenetic relationship of aphids, and much more aphid data are needed in future studies.

The Complete Mitochondrial Genomes of Penthe kochi (Coleoptera: Tetratomidae) with Its Phylogenetic Implications

To explore the mitogenome characteristics of Tetratomidae and the phylogenetic position of this family in Tenebrionoidea, the mitogenome of Penthe kochi Mařan, 1940 was sequenced, annotated, and analyzed. The P. kochi mitogenome is consistent with Tenebrionoidea species in gene length, genomic organization, codon usage, and secondary structures of transfer genes (tRNAs). Most protein-coding genes (PCGs) originate with a typical ATN start codon, except nad1 and nad3, which start with TTG. In total, 10 PCGs are terminated with complete stop codon TAA and TAG, while cox1, cox2, and nad 4 contain an incomplete stop codon T-. Among the 13 PCGs, nad2 (Pi = 0.282) has the most diverse nucleotide composition, and cox2 is the most conserved gene with the lowest value (Pi = 0.154). The Ka/Ks ratio of cox1 (0.076) and cox2 (0.124) has a lower value. All the tRNAs can be folded in a typical clover-leaf secondary structure, except trnS1, which lacked a dihydrouridine arm. And phylogenetic analyses were performed based on 13 PCGs using the Bayesian inference (BI) method. The results showed that the clade of Tenebrionoidea was well separated from the outgroups, and Tetratomidae and Mycetophagidae were not well resolved. Phylogenetic analyses with more mitogenome samplings are needed to resolve the phylogeny of Tenebrionoidea.

The first report of complete mitogenomes of two endangered species of genus Propomacrus (Coleoptera: Scarabaeidae: Euchirinae) and phylogenetic implications

To understand the mitochondrial genome structure of two endangered and long-armed scarab beetles, Propomacrus davidi and Propomacrus bimucronatus, their complete mitogenomes were sequenced for the first time in this study. The complete mitogenomes of P. davidi and P. bimucronatus were 18, 042 bp and 18, 104 bp in length, respectively. The gene orders of their mitogenomes were highly consistent with other Coleopteran species, and the typical ATN was used as the start codon in most protein coding genes. The incomplete stop codon T was used in cox1, cox2, and nad5, and TAN was used as a complete stop codon in most protein coding genes. All predicted tRNAs could form a typical cloverleaf secondary structure, except that trnS1 lacked the dihydrouridine arm. Based on the maximum likelihood and the Bayesian inference methods, phylogenetic trees of 50 species were reconstructed. The results showed that P. davidi, P. bimucronatus, Cheirotonus jansoni and Cheirotonus gestroi clustered in the same branch, and were the most closely related. The results supported that subfamily Euchirinae is a monophyletic group of Scarabaeidae, which was consistent with the morphological classification. These molecular data enriched the complete mitogenome database of Euchirinae, and improved our understanding of the phylogenetic relationship and evolutionary characteristics of these two endangered species.

The Fast Evolution of the Stenobothrini Grasshoppers (Orthoptera, Acrididae, and Gomphocerinae) Revealed by an Analysis of the Control Region of mtDNA, with an Emphasis on the Stenobothrus eurasius Group

The two cryptic grasshopper species of the genus Stenobothrus, S. eurasius and S. hyalosuperficies, demonstrate different acoustic behavior despite a strong similarity in morphology. A hybridization between these species is possible in the contact zone; however, there are little molecular data about the relationships of these species. The analysis of the mtDNA control region (CR) reveals that haplotypes of S. hyalosuperficies have more in common with the more distant Stenobothrus species than with the closely related S. eurasius. In the contact zone, S. eurasius has mt-haplotypes shared with S. hyalosuperficies, which might indicate an introgression of mtDNA from S. hyalosuperficies to the S. eurasius gene pool. We also analyze the structure and evolutionary rate of the mtDNA CR for the Stenobothrus genus and estimate the time of divergence of the species within the genus. The phylogenetic tree of the tribe Stenobothrini reconstructed with either the CR or COI gave the same four groups. The phylogenetic tree of the Stenobothrus genus has a star-like topology with each mtDNA haplotype found in any analyzed species, except S. eurasius, which forms a separate branch. The maximum degree of incomplete lineage sorting can demonstrate either ancestral polymorphism or introgression.

Characterization and Molecular Phylogenetic Analysis of Subfamily Erebinae (Lepidoptera: Noctuoidea: Erebidae) Using Five Complete Mitochondrial Genomes

In this study, the complete mitogenomes of Sympis rufibasis, Lacera noctilio, Oxyodes scrobiculata, Mocis undata, and Artena dotata were newly sequenced to bring up-to-date the database using the next-generation sequencing methods. The gene order of all sequenced mitogenomes was identical consisting of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a non-coding A+T-rich region, which were common to other Lepidopteran insects. All protein-coding genes (PCGs) initiated with a canonical ATN codon and ended with TAN or an incomplete stop codon, single T. The A+T-rich region of S. rufibasis, L. noctilio, O. scrobiculata, M. undata, and A. dotata are 406 bp, 462 bp, 372 bp, 410 bp, and 406 bp long, respectively, containing number of characteristics that are distinctive to Noctuoidea moths. We analyzed concatenated amino acid sequences of protein-coding genes not including rRNAs, using Maximum Likelihood and Bayesian Inference methods. The phylogenetic analyses indicated that the tribe relationships within Erebinae were reconstructed as (Sypnini+((Erebini 1+Poaphilini 1)+((Euclidiini+Catocalini+(Hypopyrini+Erebini 2))+((Hulodini+(Poaphilini 2+Ophiusini))))). Phylogenetic analyses supported and confirmed the monophyly of the subfamilies' relationships as follows: (Hypeninae+Lymantriinae)+((Scoliopterginae+((Calpinae+Erebinae)+((Herminiinae+Aganainae)+Arctiinae)))) within Erebidae.

Novel insight into lepidopteran phylogenetics from the mitochondrial genome of the apple fruit moth of the family Argyresthiidae

BACKGROUND

The order Lepidoptera has an abundance of species, including both agriculturally beneficial and detrimental insects. Molecular data has been used to investigate the phylogenetic relationships of major subdivisions in Lepidoptera, which has enhanced our understanding of the evolutionary relationships at the family and superfamily levels. However, the phylogenetic placement of many superfamilies and/or families in this order is still unknown. In this study, we determine the systematic status of the family Argyresthiidae within Lepidoptera and explore its phylogenetic affinities and implications for the evolution of the order. We describe the first mitochondrial (mt) genome from a member of Argyresthiidae, the apple fruit moth Argyresthia conjugella. The insect is an important pest on apples in Fennoscandia, as it switches hosts when the main host fails to produce crops.

RESULTS

The mt genome of A. conjugella contains 16,044 bp and encodes all 37 genes commonly found in insect mt genomes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a large control region (1101 bp). The nucleotide composition was extremely AT-rich (82%). All detected PCGs (13) began with an ATN codon and terminated with a TAA stop codon, except the start codon in cox1 is ATT. All 22 tRNAs had cloverleaf secondary structures, except trnS1, where one of the dihydrouridine (DHU) arms is missing, reflecting potential differences in gene expression. When compared to the mt genomes of 507 other Lepidoptera representing 18 superfamilies and 42 families, phylogenomic analyses found that A. conjugella had the closest relationship with the Plutellidae family (Yponomeutoidea-super family). We also detected a sister relationship between Yponomeutoidea and the superfamily Tineidae.

CONCLUSIONS

Our results underline the potential importance of mt genomes in comparative genomic analyses of Lepidoptera species and provide valuable evolutionary insight across the tree of Lepidoptera species.

The complete mitochondrial genome of Echinolaelaps fukienensis provide insights into phylogeny and rearrangement in the superfamily Dermanyssoidea

BACKGROUND

Echinolaelaps fukienensis is the dominant mite species parasitic on the body surface of the genus Niviventer. The mitochondrial genome (mitogenome) has its own independent genetic material and genetic system, and is now widely used in population genetics, genealogical biogeography, phylogeny and molecular evolution studies. Species diversity of the superfamily Dermanyssoidea is very rich, but its mitogenomes AT content is high, and it is difficult to amplify the complete mitogenome by routine PCR. To date, we have only obtained the mitogenomes of 6 species, scarcity on sequence data has greatly impeded the studies in the superfamily Dermanyssoidea.

METHODS

Echinolaelaps fukienensis were collected in 2019 from the body surface of Niviventer confucianus (Rodentia, Muridae) in Yunnan Province. The E. fukienensis mitogenome was determined and analyzed for the first time using the Illumina Novoseq 6000 platform. Phylogenetic analyses of the superfamily Dermanyssoidea were conducted based on the entire mitogenome sequences.

RESULTS

The E. fukienensis mitogenome was 14,402 bp, which is known the smallest genome of the superfamily Dermanyssoidea, encoding a total of 37 genes, including 13 PCGs, 22 tRNAs, 2 rRNAs and 1 control region. Most protein-coding genes use ATN as the start codon and TAN as the stop codon. AT and GC skew of atp8 genes in E. fukienensis were both 0. The average length of 22 tRNA genes of E. fukienensis was 64 bp, and secondary structures of tRNAs showed base mismatches and missing D-arms in many places. Compared with gene arrangement pattern of the hypothetical ancestor of arthropods, the E. fukienensis mitogenome shows a novel arrangement pattern. Phylogenetic tree supported the monophyly of the superfamily Dermanyssoidea. Echinolaelaps fukienensis being the least genetic distant (0.2762) and most closely related to Varroa destructor.

CONCLUSIONS

This study analyzed comprehensive the structure and evolution of the E. fukienensis mitogenome for the first time, enriches molecular data of the genus Echinolaelaps, which will contribute to further understand phylogeny and rearrangement patterns of the superfamily Dermanyssoidea.

The analyses of the complete mitochondrial genomes of three crabs revealed novel gene rearrangements and phylogenetic relationships of Brachyura

BACKGROUND

Brachyura crab is the largest branch of Decapoda crustacean. Phylogenetic relationships within Brachyura remain controversial to be investigated. The mitochondrial genome (mitogenome) is an important molecular marker for studying the phylogenetic relationships of Brachyura.

METHODS AND RESULTS

To understand the phylogeny of Brachyura, the three complete mitogenomes from Charybdis annulata, Leptodius exaratus, and Spider crab were sequenced and annotated. Their full length was 15,747, 15,716, and 16,608 bp long, respectively. The first two crabs both contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. However, Spider crab contained 13 PCGs, two rRNA genes, 25 tRNA genes and a control region. The mitogenomes of each of the three crabs exhibited high AT content (67.8%, 69.1%, and 70.8%), with negative AT skews (-0.014, - 0.028, and - 0.017) and GC skews (-0.269, - 0.286, and - 0.341). The gene order of C. annulata was identical to the ancestor of Brachyura. Compared with the ancestor of Brachyura, L. exaratus exhibited the gene rearrangements of Val (V)-rrnS-control region, and Spider crab had the four copies of Lys (K). Phylogenetic analyses indicated that C. annulata belonged to Portunidae family, Portunoidea superfamilies, L. exaratus belonged to Xanthidae family, Xanthoidea superfamilies, and Spider crab belonged to Mithracidae family, Majoidea superfamilies. Phylogenetic analyses showed that the two species (Somanniathelphusa boyangensis and Huananpotamon lichuanense) belonging to the Potamoidea were sister groups to the Thoracotremata, thus supporting the conclusion that Heterotremata is polyphyletic.

CONCLUSION

The results of this study enriched the crab mitogenome database and enabled us to better understand the phylogenetic relationships of Brachyura.

New Mitogenomes of the Green Lacewing Tribe Ankylopterygini (Neuroptera: Chrysopidae: Chrysopinae) and Phylogenetic Implications of Chrysopidae

Chrysopidae (green lacewings) are a cosmopolitan and species-rich family of Neuroptera, with remarkable significance of biological control against various agricultural and forestry pests. However, the phylogenetic position of Chrysopidae in Neuroptera and the internal relationships within the family remain equivocal among previous studies based on different types of data and sampling. Here we sequenced the mitochondrial genomes (mitogenomes) of two species of the genus Ankylopteryx in the chrysopine tribe Ankylopterygini for the first time. The characteristics of these mitogenomes were analyzed in comparison with other green lacewing mitogenomes published to date. In the phylogeny herein reconstructed based on mitogenomes, Chrysopinae were recovered as the sister group to Apochrysinae + Nothochrysinae. Within the subfamily of Chrysopinae, Nothancylini were recovered as the sister group to (Leucochrysini + Belonopterygini) + (Ankylopterygini + Chrysopini). The divergence time estimation suggested an Early Cretaceous initial divergence within the extant Chrysopidae. Within Chrysopinae, the four tribes except Nothancylini diverged around mid-Cretaceous.

Novel tRNA Gene Rearrangements in the Mitochondrial Genomes of Poneroid Ants and Phylogenetic Implication of Paraponerinae (Hymenoptera: Formicidae)

Ants (Formicidae) are the most diverse eusocial insects in Hymenoptera, distributed across 17 extant subfamilies grouped into 3 major clades, the Formicoid, Leptanilloid, and Poneroid. While the mitogenomes of Formicoid ants have been well studied, there is a lack of published data on the mitogenomes of Poneroid ants, which requires further characterization. In this study, we first present three complete mitogenomes of Poneroid ants: Paraponera clavata, the only extant species from the subfamily Paraponerinae, and two species (Harpegnathos venator and Buniapone amblyops) from the Ponerinae subfamily. Notable novel gene rearrangements were observed in the new mitogenomes, located in the gene blocks CR-trnM-trnI-trnQ-ND2, COX1-trnK-trnD-ATP8, and ND3-trnA-trnR-trnN-trnS1-trnE-trnF-ND5. We reported the duplication of tRNA genes for the first time in Formicidae. An extra trnQ gene was identified in H. venator. These gene rearrangements could be explained by the tandem duplication/random loss (TDRL) model and the slipped-strand mispairing model. Additionally, one large duplicated region containing tandem repeats was identified in the control region of P. clavata. Phylogenetic analyses based on protein-coding genes and rRNA genes via maximum likelihood and Bayes methods supported the monophyly of the Poneroid clade and the sister group relationship between the subfamilies Paraponerinae and Amblyoponinae. However, caution is advised in interpreting the positions of Paraponerinae due to the potential artifact of long-branch attraction.

The complete mitochondrial genome of Trigonisca nataliae (Hymenoptera, Apidae) assemblage reveals heteroplasmy in the control region

The evolution of mitochondrial genomes in the stingless bees is surprisingly dynamic, making them a model system to understand mitogenome structure, function, and evolution. Out of the seven mitogenomes available in this group, five exhibit atypical characteristics, including extreme rearrangements, rapid evolution and complete mitogenome duplication. To further explore the mitogenome diversity in these bees, we utilized isolated mtDNA and Illumina sequencing to assemble the complete mitogenome of Trigonisca nataliae, a species found in Northern Brazil. The mitogenome of T. nataliae was highly conserved in gene content and structure when compared to Melipona species but diverged in the control region (CR). Using PCR amplification, cloning and Sanger sequencing, six different CR haplotypes, varying in size and content, were recovery. These findings indicate that heteroplasmy, where different mitochondrial haplotypes coexist within individuals, occurs in T. nataliae. Consequently, we argue that heteroplasmy might indeed be a common phenomenon in bees that could be associated with variations in mitogenome size and challenges encountered during the assembly process.

Characterization of the complete mitochondrial genome of the longhorn beetle, Batocerahorsfieldi (Coleoptera, Cerambycidae) and its phylogenetic analysis with suitable longhorn beetles

Mitochondrial genome analysis is an important tool for studying insect phylogenetics. The longhorn beetle, Batocerahorsfieldi, is a significant pest in timber, economic and protection forests. This study determined the mitochondrial genome of B.horsfieldi and compared it with the mitochondrial genomes of other Cerambycidae with the aim of exploring the phylogenetic status of the pest and the evolutionary relationships among some Cerambycidae subgroups. The complete mitochondrial genome of B.horsfieldi was sequenced by the Illumina HiSeq platform. The mitochondrial genome was aligned and compared with the existing mitochondrial genomes of Batoceralineolata and B.rubus in GenBank (MF521888, MW629558, OM161963, respectively). The secondary structure of transfer RNA (tRNA) was predicted using tRNAScan-SE server v.1.21 and MITOS WebSever. Thirteen protein-coding genes (PCGs) and two ribosomal RNA gene sequences of 21 longhorn beetles, including B.horsfieldi, plus two outgroups, Dryopsernesti (Dryopidae) and Heterocerusparallelus (Heteroceridae), were analyzed. The phylogenetic tree was constructed using maximum likelihood and Bayesian inference methods. In this study, we successfully obtained the complete mitochondrial genome of B.horsfieldi for the first time, which is 15 425 bp in length. It contains 37 genes and an A + T-rich region, arranged in the same order as the recognized ancestor of longhorn beetles. The genome of B.horsfieldi is composed of 33.12% A bases, 41.64% T bases, 12.08% C bases, and 13.16% G bases. The structure, nucleotide composition, and codon usage of the new mitochondrial genome are not significantly different from other longhorn mitochondrial genomes. Phylogenetic analyses revealed that Cerambycidae formed a highly supported single clade, and Vesperidae was either clustered with Cerambycidae or formed a separate clade. Interestingly, B.horsfieldi, B.rubus and B.lineolata were clustered with Monochamus and Anoplophora species in both analyses, with high node support. Additionally, the VesperidaeSpiniphilusspinicornis and Vesperussanzi and the 19 Cerambycidae species formed a sister clade in the Bayesian analysis. Our results have produced new complete mitogenomic data, which will provide information for future phylogenetic and taxonomic research, and provide a foundation for future relevant research.

The First Complete Mitochondrial Genome of Genus Isocapnia (Plecoptera: Capniidae) and Phylogenetic Assignment of Superfamily Nemouroidea

Capniidae are a family of stoneflies, also known as snow flies, who emerge in winter. The phylogeny of Capniidae is widely accepted to be based on morphological analysis. Until now, only five Capniidae mitochondrial genomes have been sequenced so far. In addition, sampling is required to determine an accurate phylogenetic association because the generic classification of this family is still controversial and needs to be investigated further. In this study, the first mitogenome of genus Isocapnia was sequenced with a length of 16,200 bp and contained 37 genes, including a control region, two rRNAs, 22 tRNAs, and 13 PCGs, respectively. Twelve PCGs originated with the common start codon ATN (ATG, ATA, or ATT), while nad5 used GTG. Eleven PCGs had TAN (TAA or TAG) as their last codon; however, cox1 and nad5 had T as their final codon due to a shortened termination codon. All tRNA genes demonstrated the cloverleaf structure, which is distinctive for metazoans excluding the tRNASer1 (AGN) that missed the dihydrouridine arm. A Phylogenetic analysis of the superfamily Nemouroidea was constructed using thirteen PCGs from 32 formerly sequenced Plecoptera species. The Bayesian inference and maximum likelihood phylogeny tree structures derived similar results across the thirteen PCGs. Our findings strongly supported Leuctridae + ((Capniidae + Taeniopterygidae) + (Nemouridae + Notonemouridae)). Ultimately, the best well-supported generic phylogenetic relationship within Capniidae is as follows; (Isocapnia + (Capnia + Zwicknia) + (Apteroperla + Mesocapnia)). These findings will enable us to better understand the evolutionary relationships within the superfamily Nemouroidea and the generic classification and mitogenome structure of the family Capniidae.

Phylogenomics including the newly sequenced mitogenomes of two moths (Noctuoidea, Erebidae) reveals Ischyja manlia (incertae sedis) as a member of subfamily Erebinae

We sequenced the mitogenomes of two Erebid species, namely Ischyja manlia (Cramer, 1776) and Rusicada privata (Walker, 1865) to analyse the phylogenetic relationship and to establish the taxonomic position of incertae sedis members of the family Erebidae. The two circular genomes of I. manlia and R. privata were 15,879 bp and 15,563 bp long, respectively. The gene order was identical, containing 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an A + T-rich region. The nucleotide compositions of the A + T-rich region of both mitogenomes were similar: 80.65% for R. privata, and 81.09% for I. manlia. The AT skew and GC skew were slightly positive in I. manlia and negative in R. privata. In I. manlia and R. privata, except for cox1 which started with CGA and TTG codons, all the other 12 PCGs started with ATN codon. The A + T-rich regions of I. manlia and R. privata were 433 and 476 bp long, respectively, and contained common characteristics of Noctuoidea moths. At present, Ischyja is treated as Erebinae incertae sedis. However, phylogenetic analysis conducted in the present study reveals that the genus Ischyja is most likely to be a member of the subfamily Erebinae.

Comparative genome and phylogenetic analysis revealed the complex mitochondrial genome and phylogenetic position of Conopomorpha sinensis Bradley

Conopomorpha sinensis Bradley is a destructive pest that causes severe economic damage to litchi and longan. Previous C. sinensis research has focused on population life tables, oviposition selectivity, pest population prediction, and control technology. However, there are few studies on its mitogenome and phylogenetic evolution. In this study, we sequenced the whole mitogenome of C. sinensis by the third-generation sequencing, and analyzed the characteristics of its mitogenome by comparative genome. The complete mitogenome of C. sinensis is a typical circular and double-stranded structure. The ENC-plot analyses revealed that natural selection could affect the information of codon bias of the protein-coding genes in the mitogenome of C. sinensis in the evolutionary process. Compared with 12 other Tineoidea species, the trnA-trnF gene cluster of tRNA in the C. sinensis mitogenome appears to have a new arrangement pattern. This new arrangement has not been found in other Tineoidea or other Lepidoptera, which needs further exploration. Meanwhile, a long AT repeated sequence was inserted between trnR and trnA, trnE and trnF, ND1 and trnS in the mitogenome of C. sinensis, and the reason for this sequence remains to be further studied. Furthermore, the results of phylogenetic analysis showed that the litchi fruit borer belonged to Gracillariidae, and Gracillariidae was monophyletic. The results will contribute to an improved understanding of the complex mitogenome and phylogeny of C. sinensis. It also will provide a molecular basis for further research on the genetic diversity and population differentiation of C. sinensis.

Features and evolution of control regions in leafroller moths (Lepidoptera: Tortricidae) inferred from mitochondrial genomes and phylogeny

The control region (CR) of the mitochondrial genome (mitogenome) represents a major noncoding fragment with several special structural features that are thought to be responsible for the initiation of mitogenome transcription and replication. However, few studies have revealed the evolutionary patterns of CR in the phylogenetic context. Here, we explain the characteristics and evolution of CR in Tortricidae, inferred from a mitogenome-based phylogeny. The first complete mitogenomes of the genera Meiligma and Matsumuraeses were sequenced. Both mitogenomes are double-stranded circular DNA molecules with lengths of 15,675 bp and 15,330 bp, respectively. Phylogenetic analyses derived from 13 protein-coding genes and two ribosomal RNAs showed that most tribes, including subfamilies Olethreutinae and Tortricinae, were recovered as monophyletic clades, similar to previous studies based on morphological or nuclear data. Moreover, comprehensive comparative analyses of the structural organization and role of tandem replications on the length variation and high AT content of CR sequences were conducted. The results reveal a significant positive correlation between the total length and AT content of tandem repeats and whole CR sequences in Tortricidae. The structural organization in CR sequences is diverse, even between closely related tribes, which demonstrates the plasticity of the mitochondrial DNA molecule in Tortricidae.

Analysis of the complete mitogenome of Daphnia galeata from the Han River, South Korea: structure comparison and control region evolution

Daphnia galeata is an important plankton in aquatic ecosystems. As a widely distributed species, D. galeata has been found throughout the Holarctic region. Understanding the genetic diversity and evolution of D. galeata requires the accumulation of genetic information from different locations. Even though the mitochondrial genome (mitogenome) sequence of D. galeata has already been reported, little is known about the evolution of its mitochondrial control region. In this study, D. galeata samples were collected from the Han River on the Korean Peninsula and its partial nd2 gene was sequenced for haplotype network analysis. This analysis showed that four clades of D. galeata were present in the Holarctic region. Moreover, the D. galeata examined in this study belonged to clade D and was specific to South Korea. The mitogenome of D. galeata from the Han River showed similar gene content and structure compared to sequences reported from Japan. Furthermore, the structure of control region of the Han River was similar to those of Japanese clones and differed substantially from European clone. Finally, a phylogenetic analysis based on the amino acid sequences of 13 protein-coding genes (PCGs) indicated that D. galeata from the Han River formed a cluster with clones collected from Lakes Kasumigaura, Shirakaba, and Kizaki in Japan. The differences in control region structure and stem and loop structure reflect the different evolutionary directions of the mitogenomes from Asian and European clones. These findings improve our understanding of the mitogenome structure and genetic diversity of D. galeata.

First mitochondrial genome of subfamily Julodinae (Coleoptera, Buprestidae) with its phylogenetic implications

Complete mitochondrial genomes of three species of the family Buprestidae were sequenced, annotated, and analyzed in this study. To explore the mitogenome features of the subfamily Julodinae and verify its phylogenetic position, the complete mitogenome of Julodisvariolaris was sequenced and annotated. The complete mitogenomes of Ptosimachinensis and Chalcophorajaponica were also provided for the phylogenetic analyses within Buprestidae. Compared to the known mitogenomes of Buprestidae species varied from 15,499 bp to 16,771 bp in length, three newly sequenced mitogenomes were medium length (15,759-16,227 bp). These mitogenomes were encoded 37 typical mitochondrial genes. Among the three studied mitogenomes, Leu2 (L2), Ser2 (S2), and Pro (P) were the three most frequently encoded amino acids. Within the Buprestidae, the heterogeneity in sequence divergences of Agrilinae was highest, whereas the sequence homogeneity of Chrysochroinae was highest. Moreover, phylogenetic analyses were performed based on nucleotide matrix (13 PCGs + 2 rRNAs) among the available sequenced species of Buprestidae using Bayesian Inference and Maximum Likelihood methods. The results showed that the Julodinae was closely related to the subfamily Polycestinae. Meanwhile, the genera Melanophila, Dicerca, and Coomaniella were included in Buprestinae, which was inconsistent with the current classification system of Buprestidae. These results could contribute to further studies on genetic diversity and phylogeny of Buprestidae.

Complete Mitochondrial Genome of Scolytoplatypodini Species (Coleoptera: Curculionidae: Scolytinae) and Phylogenetic Implications

The complete mitochondrial genomes (mitogenomes) of beetles in the tribe Scolytoplatypodini (genus Scolytoplatypus) were sequenced and annotated. These included Scolytoplatypus raja (15,324 bp), Scolytoplatypus sinensis (15,394 bp), Scolytoplatypus skyliuae (15,167 bp), and Scolytoplatypus wugongshanensis (15,267 bp). The four mitogenomes contained 37 typical genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The gene orientation and arrangement of the four mitogenomes were similar to other Coleoptera mitogenomes. PCGs mostly started with ATN and terminated with TAA. The Ka/Ks ratio of 13 PCGs in the four species revealed that cox1 had the slowest evolutionary rate and atp8 and nad6 had a higher evolutionary rate. All tRNAs had typical cloverleaf secondary structures, but trnS1 lacked dihydrouridine arm. Partial tRNAs lost the discriminator nucleotide. The trnY did not possess the discriminator nucleotide and also lost three bases, showing a special amino-acyl arm. Bayesian inference (BI) and maximum likelihood (ML) methods were conducted for phylogenetic analyses using 13 PCGs. Scolytoplatypodini was clustered with Hylurgini and Hylastini, and the monophyly of Scolytoplatypodini was supported. The four newly sequenced mitogenomes increase understanding of the evolutionary relationships of Scolytoplatypodini and other Scolytinae species.

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.

Characterization of complete mitochondrial genomes of five noctuoid moths (Lepidoptera: Noctuoidea) and their phylogenetic implications

In this study, the complete mitochondrial genomes of five noctuoid moth species, namely, Paectes cristatrix, Hemichloridia euprepia, Trigonodes hyppasia, Mythimna unipuncta, and Risoba obstructa, were newly sequenced and characterized. Based on our investigations, the five mitogenomes showed sizes of 15409, 15619, 15654, 15403, and 15233 bp, respectively, and encoded 37 genes (13 protein-coding, 2 ribosomal RNA, and 22 transfer RNA genes) and a control region. All protein-coding genes (PCGs) in these mitogenomes had typical ATN start codons except for cox1 and nad1. The non-canonical codon CGA was the starter in H. euprepia and M. unipuncta and TTG was the starter in T. hyppasia for cox1. The unusual codon GTG was the starter in R. obstructa for nad1. All tRNA genes showed distinctive putative cloverleaf structures except trnS1 (AGN), which lacked the dihydrouridine (DHU) stem. Bayesian inference (BI) and maximum likelihood (ML) analyses were employed based on 13 PCGs using 159 taxa among the five families of the superfamily Noctuoidea, and the results showed the monophyly of this superfamily with well-supported values. Phylogenetic analyses showed that the newly sequenced endemic species Hemichloridia euprepia clustered separately along with the subfamilies Cuculliinae and Acronictinae with fairly supported values.

The first complete mitochondrial genome of Matsucoccidae (Hemiptera, Coccoidea) and implications for its phylogenetic position

The mitochondrial genome (mitogenome) has been extensively used to better understand the phylogenetic relationships within the hemipteran suborder Sternorrhyncha, but sequenced mitogenomes remain unavailable for the entire family Matsucoccidae to date. To address this, here we sequenced the complete mitogenome of Matsucoccusmatsumurae; the first for this family. The mitogenome is 15,360 bp in size and comprises the typical set of 37 mitochondrial genes and a large non-coding region (AT-rich region). Gene order, nucleotide composition and codon usage of protein-coding genes (PCGs) of M.matsumurae differ considerably from those of the other two sequenced Coccidae species. All PCGs were initiated by the ATN start codons and ended with the TAA/G or single T-- stop codons. Nine transfer RNA genes could be folded into typical clover-leaf secondary structures. The length and AT content of the ribosomal RNA genes are highly conserved in the Coccoidea mitogenomes. In contrast, the AT-rich control region is highly variable in size and in the number of tandem repeats. The sliding window analysis showed that the cox1 gene is the most conserved amongst the 13 PCGs, while the ratios of non-synonymous to synonymous substitution rates indicated that the evolution of this mitogenome has been dominated by positive selection. Phylogenetic analyses, based on nucleotide sequence data of 37 mitochondrial genes and amino acid sequence data of 13 PCGs using Bayesian Inference and Maximum Likelihood methods, showed that Matsucoccidae diverged before the Coccidae.

Mitochondrial genomes of Sternochetus species (Coleoptera: Curculionidae) and the phylogenetic implications

The three weevil species, Sternochetus gravis, S. mangiferae, and S. olivieri, have all been reported to be serious pests of mango fruits. Morphology, biology, and various management approaches of these economically important weevils have been well studied. However, no mitochondrial genomes have been reported from the genus Sternochetus. Herein, we assembled mitogenomes of all the three Sternochetus species to reveal their mitogenomic characteristics. A DNA library of 350 bp insert size was constructed and sequenced in Illumina's HiSeq 6000 platform with a pair-end 150 bp sequencing strategy by Novogene. The sequence reads were assembled using GetOrganelle v1.7.1 and the genes were annotated by Geneious Prime 2021.0.3 and MITOS Web Server. Coupled with 61 published mitogenomes from 13 subfamilies of Curculionidae, we reconstructed phylogenetic trees to resolve evolutionary relationships of these closely related species and also examined subfamily-level classification among Curculionidae. All three mitogenomes are double-stranded circular molecules with 22 transfer RNA genes, 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 1 noncoding control region as in other insects. Higher interspecific nucleotide divergence (about 10%) of 13 PCGs indicated these three Sternochetus species diverged a long time ago. Phylogenetic analyses using both maximum likelihood and Bayesian inference methods showed that Sternochetus falls into the basal clade of Cryptorhynchini, a tribe in the subfamily Molytinae. The relationship of S. olivieri as a sister species to S. gravis + S. mangiferae was strongly supported. The monophyly of Cryptorhynchini was also well supported whereas Molytinae was suggested to be a polyphyletic group.

Comparison of mitogenomes of three Petalocephala species (Hemiptera: Cicadellidae: Ledrinae) and their phylogenetic analysis

Ledrinae is a unique group of leafhoppers with a distinct appearance. Petalocephala is the largest Ledrinae genus that is difficult to identify except by dissecting the male genitals. To date, research on Ledrinae is relatively less compared with other leafhoppers. Therefore, to better understand this group, we sequenced and analyzed three complete Petalocephala mitochondrial genomes. We comparatively analyzed these general Petalocephala genomic features (including size, AT content, AT/GC skew, 13 protein-coding gene nucleotide compositions, etc.), and predicted 22 transfer RNA secondary structures. We obtained highly consistent phylogenetic results within Cicadellidae based on mitogenomic data using the maximum likelihood and Bayesian methods. Our results showed that all subfamilies were monophyletic and had a high node support rate, and there was a sister group relationship between Ledrinae and all other leafhopper groups. Furthermore, treehoppers were found to originate from leafhoppers and showed sister group relationships with Megophthalminae. Within Ledrinae, all phylogenetic trees supporting phylogenetic relationships were as follows: ([P. dicondylica + P. gongshanensis] + [Tituria pyramidata + [Ledra auditura + P. gongshanensis]]) Based on the complete mitogenome phylogenetic analysis and the comparison of morphological characteristics, we propose that Petalocephala is not monophyletic.

Genome Organization and Comparative Evolutionary Mitochondriomics of Brown Planthopper, Nilaparvata lugens Biotype 4 Using Next Generation Sequencing (NGS)

Nilaparvata lugens is the main rice pest in India. Until now, the Indian N. lugens mitochondrial genome has not been sequenced, which is a very important basis for population genetics and phylogenetic evolution studies. An attempt was made to sequence two examples of the whole mitochondrial genome of N. lugens biotype 4 from the Indian population for the first time. The mitogenomes of N. lugens are 16,072 and 16,081 bp long with 77.50% and 77.45% A + T contents, respectively, for both of the samples. The mitochondrial genome of N. lugens contains 37 genes, including 13 protein-coding genes (PCGs) (cox1-3, atp6, atp8, nad1-6, nad4l, and cob), 22 transfer RNA genes, and two ribosomal RNA (rrnS and rrnL) subunits genes, which are typical of metazoan mitogenomes. However, both samples of N. lugens mitogenome in the present study retained one extra copy of the trnC gene. Additionally, we also found 93 bp lengths for the atp8 gene in both of the samples, which were 60–70 bp less than that of the other sequenced mitogenomes of hemipteran insects. The phylogenetic analysis of the 19 delphacids mitogenome dataset yielded two identical topologies when rooted with Ugyops sp. in one clade, and the remaining species formed another clade with P. maidis and M. muiri being sisters to the remaining species. Further, the genus Nilaparvata formed a separate subclade with the other genera (Sogatella, Laodelphax, Changeondelphax, and Unkanodes) of Delphacidae. Additionally, the relationship among the biotypes of N. lugens was recovered as the present study samples (biotype-4) were separated from the three biotypes reported earlier. The present study provides the reference mitogenome for N. lugens biotype 4 that may be utilized for biotype differentiation and molecular-aspect-based future studies of N. lugens.

Comparative analysis of mitochondrial genomes among the family Peltoperlidae (Plecoptera: Systellognatha) and phylogenetic implications

Nowadays, the position of Peltoperlidae in Systellognatha has been resolved based on morphological analyses. However, there are different opinions based on molecular data. To date, only three peltoperlid mitogenomes are available, and more sampling is needed to obtain precise phylogenetic relationships. In this study, we obtained the complete mitogenomes of Cryptoperla kawasawai (15,832 bp) and Peltoperlopsis sagittata (15,756 bp). Our results show that gene content, gene order, DmTTF binding site, nucleotide composition, codon usage, ribonucleic acid (RNA) structure, and structural elements in the control region are highly conserved in peltoperlids. Heatmap analysis of codon usage shows that the AT-rich codons UUA, AUU, UUU, and AUA were commonly used codons in the Peltoperlidae. Evolutionary rate analyses of protein-coding genes reveal that different genes have been subject to different rates of molecular evolution correlated with the GC content. All tRNA genes in peltoperlid mitogenomes have a canonical cloverleaf secondary structure except for trnS1, whose dihydrouridine arm simply forms a loop. The control region of the family has several distinct structural characteristics and has the potential to serve as effective phylogenetic markers. Phylogenetic analyses support the monophyly of Perloidea, but the monophyly of Pteronarcyoidea is still not supported. The Peltoperlidae is placed as the earliest branch within the Systellognatha, and the estimated phylogenetic relationship is: Peltoperlidae + {(Styloperlidae + Pteronarcyidae) + [Perlidae + (Chloroperlidae + Perlodidae)]}. Our results provide new insight into the phylogeny of this group.

Complete mitochondrial genomes of four deep-sea echinoids: conserved mitogenome organization and new insights into the phylogeny and evolution of Echinoidea

Echinoids are an important component in benthic marine environments, which occur at all depths from the shallow-water hard substrates to abyssal depths. To date, the phylogeny of the sea urchins and the macro-evolutionary processes of deep-sea and shallow water groups have not yet been fully resolved. In the present study, we sequenced the complete mitochondrial genomes (mitogenomes) of four deep-sea sea urchins (Echinoidea), which were the first representatives of the orders Aspidodiadematoida, Pedinoida and Echinothurioida, respectively. The gene content and arrangement were highly conserved in echinoid mitogenomes. The tRNA-Ser ^AGY with DHU arm was detected in the newly sequenced echinoid mitogenomes, representing an ancestral structure of tRNA-Ser ^AGY. No difference was found between deep-sea and shallow water groups in terms of base composition and codon usage. The phylogenetic analysis showed that all the orders except Spatangoida were monophyletic. The basal position of Cidaroida was supported. The closest relationship of Scutelloida and Echinolampadoida was confirmed. Our phylogenetic analysis shed new light on the position of Arbacioida, which supported that Arbacioida was most related with the irregular sea urchins instead of Stomopneustoida. The position Aspidodiadematoida (((Aspidodiadematoida + Pedinoida) + Echinothurioida) + Diadematoida) revealed by mitogenomic data discredited the hypothesis based on morphological evidences. The macro-evolutionary pattern revealed no simple onshore-offshore or an opposite hypothesis. But the basal position of the deep-sea lineages indicated the important role of deep sea in generating the current diversity of the class Echinoidea.

Complete Nucleotide Sequence of the Mitogenome of Tapinoma ibericum (Hymenoptera: Formicidae: Dolichoderinae), Gene Organization and Phylogenetics Implications for the Dolichoderinae Subfamily

The ant Tapinoma ibericum Santschi, 1925 is native to the Iberian Peninsula. This species, as well as other species from the Tapinoma nigerrimum complex, could form supercolonies that make these species potentially invasive and could give rise to pests. Recently a mature colony from this species has been found in the Isle of Wight (United Kingdom). Mitogenomes have been used to study the taxonomy, biogeography and genetics of species, improving the development of strategies against pest invasion. However, the number of available mitogenomes from the subfamily Dolichoderinae is still scarce and only two of these mitogenomes belong to Tapinoma species. Herein, the complete mitogenome of T. ibericum is presented in order to increase the molecular information of the genus. The T. ibericum mitogenome, retrieved by Next-Generation Sequencing data, is 15,715 bp in length. It contains the typical set of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNAs and the A + T-rich control region. Comparisons of the T. ibericum mitogenome with other dolichoderine mitogenomes revealed the existence of four gene rearrangements in relation with the ancestral insect mitogenome. One of these rearrangements, involving the tRNA-Ile, tRNA-Gln and tRNA-Met genes, was found in most of the analyzed ant mitogenomes. Probably this rearrangement was an ancestral or plesiomorphic character in Formicidae. Interestingly, another rearrangement that affects to tRNA-Trp, tRNA-Cys and tRNA-Tyr genes was found only in Tapinoma species. This change could be a synapomorphic character for the genus Tapinoma, and could be used as a phylogenetic marker. Additionally, a phylogenetic analysis was performed using the protein-coding gene sequences from available Dolichoderinae mitogenomes, as well as mitogenomes from representative species from other Formicidae subfamilies. Results support the monophyletic nature of the genus Tapinoma placing it within the same clade as the rest of Dolichoderinae species.

Mitochondrial Genome of Strophopteryx fasciata (Plecoptera: Taeniopterygidae), with a Phylogenetic Analysis of Nemouroidea

Taeniopterygidae is a medium-sized family of stoneflies. The phylogeny of Taeniopterygidae was widely accepted based on the morphological analyses. However, there are different opinions based on molecular data. To date, only two taeniopterygid mitochondrial genomes (mitogenomes) were available, and more sampling is needed to obtain precise phylogenetic relationships. In this research, the Strophopteryx fasciata mitogenome was sequenced and analyzed. The complete mitogenome of S. fasciata was 15,527 bp in length and contained 37 genes and a non-coding control region. Among taeniopterygid mitogenomes, the length variation was minimal in protein-coding genes (PCGs), transfer RNA genes (tRNAs) and ribosomal RNA genes (rRNAs), but very different in the control region. Similar to mitogenomes of other taeniopterygid species, the S. fasciata mitogenome was consistently AT biased and displayed positive AT- and negative GC-skews of the whole mitogenome. Most PCGs used ATN as the start codon and TAA/TAG as the stop codon. The stop codons were far less variable than the start codons in taeniopterygid mitogenomes. All Ka/Ks ratios were less than 1, indicating the presence of purifying selection in these genes. The secondary structures of transfer and ribosomal RNA genes of S. fasciata mitogenome are highly conserved with other taeniopterygid species. In the control region of the S. fasciata mitogenome, some essential elements (tandem repeats, stem–loop structures, and poly−N stretch, etc.) were observed. Two phylogenetic trees were inferred from Bayesian inference (BI) and Maximum Likelihood (ML) methods generated the identical topology across the PCGR dataset. The relationships of five families in Nemouroidea were recovered as Leuctridae + ((Capniidae + Taeniopterygidae) + (Nemouridae + Notonemouridae)). These results will help us understand the mitogenome structure of taeniopterygid species and the evolutionary relationship within Plecoptera.

The Complete Mitochondrial Genome of the Chicken Body Louse, Menacanthus cornutus, and Evolutionary Patterns of Extensive Gene Rearrangements in the Mitochondrial Genomes of Amblycera (Psocodea: Phthiraptera)

Animal mitochondrial (mt) genomes are typically double-strand circular DNA molecules, but diverse structural variations have been widely found in multiple groups. In parasitic lice (Phthiraptera), the structure of mt genomes varies remarkably across all five suborders. In this study, we reported the complete mt genome of a chicken body louse, Menacanthus cornutus, which has a typical single circular mt chromosome and drastic mt gene rearrangements. This mt genome is 15,693 bp in length, consisting of 13 protein-coding genes, 23 tRNA genes, 2 rRNA genes, and a control region. A comparison with a typical insect mt genome suggested that two highly similar trnM are present in the mt genome of M. cornutus. Moreover, almost every single gene was rearranged, and over half of mt genes were inverted. Phylogenetic analyses inferred from the mt genome sequences supported the monophyly and position of Amblycera. Mapped over the phylogenetic relationships of Amblycera, we identified two inversion events for the conserved gene blocks in Boopidae and Menoponidae. The inverted ND4L-ND4 was likely a synapomorphic rearrangement in Menoponidae. Our study demonstrated the importance of sequencing mt genomes for additional taxa to uncover the mechanism underlying the structural evolution of the mt genome in parasitic lice.

Comparative analysis of four complete mitogenomes from hoverfly genus Eristalinus with phylogenetic implications

The genus Eristalinus is widely distributed globally. Four complete mitochondrial genomes (i.e., mitogenomes) of Eristalinus were sequenced and analyzed in this study: Eristalinus viridis (Coquillett, 1898), E. quinquestriatus (Fabricius, 1781), E. tarsalis (Macquart, 1855), and E. sp. Within these four sequenced mitogenomes, most protein-coding genes (ND2, CO1, COX2, COX3, ND3, ND5, ND4, ND4L, ND6, and Cytb) began with a typical ATN (T/C/G/A) start codon and ended with a stop codon TAA or incomplete T, whereas ND1 began with the start codon TTG. ND3 ended with TAG. The secondary tRNA structure was that of a typical cloverleaf, and only the tRNA-Ser1 lacked a DHU arm. Three and five domains appeared in the 12S and 16S rRNA secondary structures, respectively. The phylogenetic relationships among the four Eristalinus species combined with the published mitogenomes of Syrphidae were reconstructed using the maximum likelihood and Bayesian inference methods, which support the monophyly of the subfamily Syrphinae but do not support that of the subfamily Eristalinae. Of note, Eristalini and Syrphini are monophyletic groups. The mitogenomes of E. viridis, E. quinquestriatus, E. sp., and E. tarsalis are useful for determining the phylogenetic relationships and evolution of Syrphidae.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

Mitogenomes of Nine Asian Skipper Genera and Their Phylogenetic Position (Lepidoptera: Hesperiidae: Pyrginae)

In this study, complete mitochondrial genomes of nine species representing three tribes in the subfamily Pyrginae sensu lato were newly sequenced. The mitogenomes are closed double-stranded circular molecules, with the length ranging from 15,232 bp to 15,559 bp, which all encode 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region. The orientation and gene order of these nine mitogenomes are identical to the inferred ancestral arrangement of insects. All PCGs exhibit the typical start codon ATN except for cox1 (using CGA) and cox2 (using TTG) in Mooreana trichoneura. Most of the PCGs terminate with a TAA stop codon, while cox1, cox2, nad4, and nad5 end with the incomplete codon single T. For the different datasets, we found that the one comprising all 37 genes of the mitogenome produced the highest nodal support, indicating that the inclusion of RNAs improves the phylogenetic signal. This study re-confirmed the status of Capila, Pseudocoladenia, and Sarangesa; namely, Capila belongs to the tribe Tagiadini, and Pseudocoladenia and Sarangesa to the tribe Celaenorrhini. Diagnostic characters distinguishing the two tribes, the length of the forewing cell and labial palpi, are no longer significant. Two populations of Pseudocoladenia dan fabia from China and Myanmar and P. dan dhyana from Thailand are confirmed as conspecific.

The First Eight Mitogenomes of Leaf-Mining Dactylispa Beetles (Coleoptera: Chrysomelidae: Cassidinae) Shed New Light on Subgenus Relationships

The taxonomic classification of Dactylispa, a large genus of leaf-mining beetles, is problematic because it is currently based on morphology alone. Here, the first eight mitochondrial genomes of Dactylispa species, which were used to construct the first molecular phylogenies of this genus, are reported. The lengths of the eight mitogenomes range from 17,189 bp to 20,363 bp. All of the mitochondrial genomes include 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and 1 A + T-rich region. According to the nonsynonymous/synonymous mutation ratio (Ka/Ks) of all PCGs, the highest and the lowest evolutionary rates were found for atp8 and cox1, respectively, which is a common phenomenon among animals. According to relative synonymous codon usage, UUA(L) has the highest frequency. With two Gonophorini species as the outgroup, mitogenome-based phylogenetic trees of the eight Dactylispa species were constructed using maximum likelihood (ML) and Bayesian inference (BI) methods based on the PCGs, tRNAs, and rRNAs. Two DNA-based phylogenomic inferences and one protein-based phylogenomic inference support the delimitation of the subgenera Dactylispa s. str. and Platypriella as proposed in the system of Chen et al. (1986). However, the subgenus Triplispa is not recovered as monophyletic. The placement of Triplispa species requires further verification and testing with more species. We also found that both adult body shape and host plant relationship might explain the subgeneric relationships among Dactylispa beetles to a certain degree.

Three Complete Mitochondrial Genomes of Erotylidae (Coleoptera: Cucujoidea) with Higher Phylogenetic Analysis

Simple Summary

Erotylid beetles are phytophagous and mycophagous. Phylogenetic studies on this family were mainly based on morphological characters or several gene fragments. Research on the mitochondrial genome of Erotylidae is rare. Therefore, we sequenced and analyzed three complete mt genomes of Erotylinae with a comparative mt genomic analysis of Erotylinae and Languriinae for the first time to reveal mitochondrial genome characterizations and reconstruct phylogenetic relationships of this group. The comparative analyses showed the mt genome characterizations of Erotylinae are similar to Languriinae. These results provided a comprehensive framework and worthy information for the future research of this family.

Abstract

The family Erotylidae belongs to the superfamily Cucujoidea, which are phytophagous and mycophagous. So far, only two representative complete mitochondrial (mt) genomes of Erotylidae have been sequenced. Mitochondrial genomes of Tritoma metasobrina, Neotriplax arisana, and Episcapha opaca, which all belong to the subfamily Erotylinae, were sequenced using next-generation sequencing technology to better understand the diversity of mt genomes of Erotylidae. A comparative mt genomic analysis was conducted on the three sequenced representatives of Erotylinae and Languriinae sp. (Languriinae). The size of the complete mt genome of the 4 species ranged from 15,581 bp to 16,502 bp in length, including 37 genes (13 protein-coding genes, 22 transfer RNAs, and 2 ribosomal RNAs) and the control region. The arrangements of their mt genomes are highly consistent with other Coleoptera species. The start codons of two PCGs (ND1 and ND5) and the stop codons of one PCG (ATP8) were illustrated differences between Languriinae sp. and the other three species. All tRNAs of these 4 species exhibited cloverleaf secondary structures except that the dihydorouridine (DHU) arm of tRNA^Ser(AGN) was absent. The phylogenetic analyses using both Bayesian inference (BI) and maximum likelihood (ML) methods all supported that Erotylidae as monophyletic. Erotylinae was monophyletic being the sister group to Xenocelinae. Languriinae was closely related to ‘Erotylinae-Xenocelinae’. Our results recovered Languriinae nested within Erotylidae.

The First Mitogenomes of the Subfamily Odontiinae (Lepidoptera, Crambidae) and Phylogenetic Analysis of Pyraloidea

Simple Summary

The Odontiinae is a small group in the Pyraloidea comprised of 388 species in 88 genera, but externally, these moths are diverse, including heterogeneous maculation and a size range from 9 to 50 mm in total wingspan. The monophyly of Pyraloidea and the two families (Pyralidae and Crambidae) is well supported by phylogenetic analyses based on morphology and molecular data of multiple nuclear genes. However, only a few mito-phylogenetic analyses have been conducted and no mitogenome of Odontiinae species has been reported. Three complete mitogenomes of odontiine species were sequenced and analyzed for the first time herein. The results showed that Odontiinae mitogenomes shared similar genomic characters with other Pyraloidea. The phylogenetic analyses based on 13 PCGs of mitogenomes confirmed the monophyly of Odontiinae and its position within Crambidae.

Abstract

The complete mitochondrial genomes of three species of Odontiinae were newly sequenced: Dausara latiterminalis Yoshiyasu, Heortia vitessoides (Moore), and Pseudonoorda nigropunctalis (Hampson). These circular and double-stranded mitogenomes vary from 15,084 bp to 15,237 bp in size, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs) and an A + T-rich region. The nucleotide composition indicated a strong A/T bias. Most PCGs are initiated with an ATN codon and terminated by a codon of TAR. All tRNAs could be folded into the clover-leaf structure with the exception of trnS1 (AGN), in which the dihydrouridine (DHU) arm formed a simple loop, and the motif ‘ATAG’ and ‘ATTTA’ in the A + T-rich region was also founded. The phylogenomic analyses covering Odontiinae + 11 subfamilies of Pyraloidea were conducted. Similar topologies were generated from both Bayesian inference (BI) and maximum likelihood (ML) analyses based on the nucleotide and amino acid sequence data. There was some discrepancy in the sister-group relationship of Odontiinae and Glaphyriinae, and the relationships among the subfamilies in the ‘CAMMSS clade’ of the Crambidae. The results of this study suggest that mitogenomic data are useful for resolving the deep-level relationships of Pyraloidea and the topologies generated from amino acid data might be more realistic and reliable. Moreover, more mitogenomic taxon sampling and larger scale analyses with more genes or a combination of mitogenomic and nuclear genes are needed to reconstruct a comprehensive framework of the pyraloid phylogeny.

Mitogenome and phylogenetic analysis of typhlocybine leafhoppers (Hemiptera: Cicadellidae)

Mitogenomes have been widely used to estimate phylogenetic relationships among insects and provide data useful for augmenting traditional morphological characters in delimiting species. Here, complete mitogenome sequences of two closely related typhlocybine leafhoppers, Cassianeuracassiae (Ahmed, 1970) and C. bimaculata Dworakowska, 1984, were obtained and found to be 15,423 bp and 14,597 bp in length, respectively. The gene order was found to be similar to other published leafhopper mitogenomes, but the control region of C.bimaculata is the shortest among known leafhoppers and lacks tandem repeats. Phylogenetic analysis of 13 protein-coding genes (PCGs), the first and second codons of 13 PCGs, 13 PCGs and two rRNAs formed three well-supported tree topologies. The topologies of phylogenetic trees inferred from three datasets were almost identical, which was consistent with previous molecular phylogenies of this group. Comparative morphological study of the ovipositors revealed several characters potentially useful for diagnosing genera and resolving their phylogenetic relationships. Phylogenetic analysis of these and other morphological characters yielded a tree that is mostly consistent with the tree obtained from analysis of mitogenome sequences. In both molecular and morphological phylogenenies, Typhlocybini and Zyginellini clustered into one clade, but neither was recovered as monophyletic.

First mitogenome of subfamily Langiinae (Lepidoptera: Sphingidae) with its phylogenetic implications

To date, a relatively complete classification of Sphingidae (Lepidoptera) has been generated, but the phylogeny of the family remains need to be fully resolved. Some phylogenetic relationships within Sphingidae still remains uncertain, especially the taxonomic status of the subfamily Langiinae and its sole included genus and species, Langia zenzeroides. To begin to address this problem, we generated nine new complete mitochondrial genomes, including that of Langia, and together with that of Theretra oldenlandiae from our previous study and 25 other Sphingidae mitogenomes downloaded from GenBank, analyzed the phylogenetic relationships of Sphingidae and investigated the mitogenomic differences among members of the Langiinae, Sphinginae, Smerinthinae and Macroglossinae. The mitogenomes of Sphingidae varied from 14995 bp to 15669 bp in length. The gene order of all newly sequenced mitogenomes was identical, containing 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and the A + T-rich region. Nucleotide composition was A + T biased, and all the protein-coding genes exhibited a positive AT-skew, which was reflected in the nucleotide composition, codon, and amino acid usage. The A + T-rich region was comprised of nonrepetitive sequences, which contained regulatory elements related to the control of replication and transcription. We analyzed concatenated gene sequences, with third codon positions of protein coding genes and rRNAs excluded, using Maximum Likelihood and Bayesian Inference techniques. All four currently recognized subfamilies were recovered as monophyletic but in contrast to the most recent studies, our preferred tree placed Langiinae as the first subfamily to diverge within Sphingidae rather as sister to Smerinthinae + Sphinginae. Our results also support the removal of the genus Barbourion from the smerinthine tribe Ambulycini to an unresolved position in "Smerinthinae incertae sedis".

Unique Duplication of trnN in Odontoptilum angulatum (Lepidoptera: Pyrginae) and Phylogeny within Hesperiidae

To explore the variation and relationship between gene rearrangement and phylogenetic effectiveness of mitogenomes among lineages of the diversification of the tribe Tagiadini in the subfamily Pyrginae, we sequenced the complete mitogenome of Odontoptilum angulatum. The genome is 15,361 bp with the typical 37 genes, a large AT-rich region and an additional trnN (trnN2), which is completely identical to trnN (sequence similarity: 100%). The gene order differs from the typical Lepidoptera-specific arrangement and is unique to Hesperiidae. The presence of a "pseudo-trnS1" in the non-coding region between trnN1 and trnN2 supports the hypothesis that the presence of an extra trnN can be explained by the tandem duplication-random loss (TDRL) model. Regarding the phylogenetic analyses, we found that the dataset comprising all 37 genes produced the highest node support, as well as a monophyly of Pyrginae, indicating that the inclusion of RNAs improves the phylogenetic signal. Relationships among the subfamilies in Hesperiidae were also in general agreement with the results of previous studies. The monophyly of Tagiadini is strongly supported. Our study provides a new orientation for application of compositional and mutational biases of mitogenomes in phylogenetic analysis of Tagiadini and even all Hesperiidae based on larger taxon sampling in the future.

First Complete Mitochondrial Genome of Melyridae (Coleoptera, Cleroidea): Genome Description and Phylogenetic Implications

To explore the characteristics of the mitogenome of Melyridae and reveal phylogenetic relationships, the mitogenome of Cordylepherus sp. was sequenced and annotated. This is the first time a complete mitochondrial genome has beengenerated in this family. Consistent with previous observations of Cleroidea species, the mitogenome of Cordylepherus sp. is highly conserved in gene size, organization and codon usage, and secondary structures oftRNAs. All protein-coding genes(PCGs) initiate withthe standard start codon ATN,except ND1, which starts with TTG, and terminate with the complete stop codons of TAA and TAG, or incomplete forms, TA- and T-. MosttRNAs have the typical clover-leaf structure, except trnS1 (Ser, AGN), whose dihydrouridine (DHU) armis reduced. In the A+T-rich region, three types of tandem repeat sequence units are found, including a 115 bp sequence tandemly repeated twice, a 16 bp sequence tandemly repeated three timeswith a partial third repeat and a 10 bp sequence tandemly repeated seven times. Phylogenetic analyses based on 13 protein-coding genes by both Bayesian inference (BI) and maximum likelihood (ML) methods suggest that Melyridae sensu lato is polyphyletic, and Dasytinae and Malchiinae are supported as independent families.

Characterization of the Complete Mitochondrial Genomes from Two Nitidulid Pests with Phylogenetic Implications

The complete mitochondrial genomes of Xenostrongylusvariegatus and Epuraea sp. were sequenced and analyzed. The total genome lengths are 17,657 and 16,641 bp, with an A+T content of 77.2% and 76.4%, respectively. Each mitochondrial genome consists of 37 coding genes and a non-coding (AT-rich) region. All protein-coding genes (PCGs) start with the standard start codon, ATN, and end with complete stop codons, TAA and TAG, or an incomplete stop codon, T. All tRNAs can be folded into the typical clover-leaf secondary structure, with the exception of trnS1 in both species with a reduced dihydrouridine (DHU) arm. The AT-rich region has tandem repeats differing in both number and length. Genetic distance and Ka/Ks analyses show that nad6 has a higher variability and more rapid evolutionary rate than other PCGs. Both maximum likelihood and Bayesian inference phylogenetic analyses based on 13 PCGs and 2 ribosome DNAs (rDNAs) agree with the previous phylogenies in supporting the Nitidulidae monophyly and the sister-group relationship of Kateretidae + (Monotomidae + Nitidulidae).

Novel tRNA gene rearrangements in the mitochondrial genome of Camarochiloides weiweii (Hemiptera: Pachynomidae)

Pachynomidae is a poorly known heteropteran group currently considered as a distinct family closely related to Reduviidae. Taxonomic position of Pachynomidae has been debated for a long time. Knowledge about mitochondrial genome (mitogenome) of this family also remain limited. In the present study, we describe the complete mitogenome of Camarochiloides weiweii, the first mitogenome sequenced for this family. This mitogenome consists of 13 protein-coding genes, two rRNA genes and 21 tRNA genes. Two tRNA gene rearrangements in this genome, including the translocation of trnA and trnR as well as the transposition of trnH, were observed and could be explained by the tandem duplication-random loss (TDRL) model. Bayesian and maximum likelihood analyses based on sequences of three different mitogenomic datasets recovered the monophyly of Reduvioidea, and the sister relationship between Reduvioidea and the remaining Cimicomorpha lineages. This study enlightens the future application of mitogenomic data in the phylogenetic study of Cimicomorpha based on the large-scale sampling.

Characterization of the Complete Mitochondrial Genomes of Two Species with Preliminary Investigation on Phylogenetic Status of Zyginellini (Hemiptera: Cicadellidae: Typhlocybinae)

To explore the characteristics of mitogenomes and reveal phylogenetic relationships of the tribes of Zyginellini and Typhlocybini in Typhlocybinae, mitogenomes of two species of the Zyginellini, Parazyginella tiani and Limassolla sp., were sequenced. Mitogenomes of both species contain 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a large non-coding region (A + T-rich region). These characteristics are similar to other Membracoidea mitogenomes. All PCGs initiate with the standard start codon of ATN and terminate with the complete stop codon of TAA/G or with an incomplete T codon. All tRNAs have the typical clover-leaf structure, except trnS1 which has a reduced DHU arm and the acceptor stem of trnR is 5 or 6 bp in some species, an unusual feature here reported for the first time in Typhlocybinae. The A + T-rich region is highly variable in length and in numbers of tandem repeats present. Our analyses indicate that nad6 and atp6 exhibit higher evolutionary rates compared to other PCGs. Phylogenetic analyses by both maximum likelihood and Bayesian methods based on 13 protein-coding genes of 12 species of Typhlocybinae suggest that Zyginellini are paraphyletic with respect to Typhlocybini.

Characterization of the Complete Mitochondrial Genome of Epicauta impressicornis (Coleoptera: Meloidae) and Its Phylogenetic Implications for the Infraorder Cucujiformia

The complete mitochondrial genome (mitogenome) of Epicauta impressicornis Pic (Coleoptera: Meloidae) was determined. The circular genome is 15,713-bp long, and encodes 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region (CR). The 13 PCGs start with the typical ATN codon and terminate with the typical stop codon TAA (ND2, ND4L, ND6, ATP6, ATP8, and CYTB), TAG (ND1 and ND3), and T- (COX1, COX2, COX3, ND4, and ND5). The two rRNA genes (rrn12S and rrn16S) are encoded on the minority strand. All tRNAs genes except trnS1 (AGN) are predicted to fold into the typical cloverleaf structure. The longest overlap (10 bp) is observed between ATP8 and ATP6. CR mainly harbors a conserved poly-T stretch (15 bp), a short repeat unit (17 bp), some universal microsatellite-like repeats, and a canonical poly-A tail. Phylogenetic analysis using Bayesian inferences and maximum likelihood based on nucleotide and corresponding amino acid sequences of the 13 PCGs showed that E. impressicornis is closely related to E. chinensis, this relationship is and supported within Cucujiformia belonging to Meloidae (Tenebrionoidea). Our results further confirmed the monophyly of Tenebrionoidea, Lymexyloidea, Curculionoidea, Chrysomeloidea, Cucujoidea, Coccinelloidea, and Cleroidea within Cucujiformia, and revealed the sister relationships of (Cleroidea + Coccinelloidea), (Lymexyloidea + Tenebrionoidea), and ((Chrysomeloidea + Cucujoidea) + Curculionoidea). We believe that the complete mitogenome of E. impressicornis will contribute to further studies on molecular bases for the classification and phylogeny of Meloidae or even Cucujiformia.

Characteristic and variability of five complete aphid mitochondrial genomes: Aphis fabae mordvilkoi, Aphis craccivora, Myzus persicae, Therioaphis tenera and Appendiseta robiniae (Hemiptera; Sternorrhyncha; Aphididae)

The complete mitochondrial genomes of aphids Aphis fabae mordvilkoi, A. craccivora, Myzus persicae from Aphidinae as well as Therioaphis tenera and Appendiseta robiniae from Calaphidinae were sequenced and compared with the genomes of other aphid species. A. fabae mordvilkoi, Th. tenera and A. robiniae mitogenomes were sequenced and analyzed for the first time. The annotation of A. craccivora and M. persicae were corrected compared to what was previously published. According to our data there is no translocation of tRNA-Tyr gene in A. craccivora mitogenome and this aphid species has an ancestral type of mitochondrial gene order. A + T content in all 5 mitogenomes was higher than 80%. A + T content in the Th. tenera CR was 59.5% which is untypically low. CRs of all 5 studied mitogenomes had 2 conserved motifs at their ends and extended G + C rich region. A. craccivora, M. persicae and Th. tenera had large tandem repeats inside the CRs. Detailed molecular analysis of all 5 aphid mitochondrial genomes showed the importance of a deep understanding of the molecular organization of all the functional regions of the mitochondrial DNA, which helps to avoid mistakes during genome annotation.

The complete mitochondrial genome of stag beetle Lucanus cervus (Coleoptera: Lucanidae) and phylogenetic analysis

Background

The stag beetle Lucanus cervus (Coleoptera: Lucanidae) is widely distributed in Europe. Habitat loss and fragmentation has led to significant reductions in numbers of this species. In this study, we sequenced the complete mitochondrial genome of L. cervus and reconstructed phylogenetic relationships among Lucanidae using complete mitochondrial genome sequences.

Methods

Raw data sequences were generated by the next generation sequencing using Illumina platform from genomic DNA of L. cervus. The mitochondrial genome was assembled by IDBA and annotated by MITOS. The aligned sequences of mitochondrial genes were partitioned using PartitionFinder 2. Phylogenetic relationships among 19 stag beetle species were constructed using Maximum Likelihood (ML) method implemented in IQ-TREE web server and Bayesian method implemented in PhyloBayes MPI 1.5a. Three scarab beetles were used as outgroups.

Results

The complete mitochondrial genome of L. cervus is 20,109 bp in length, comprising 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNAs and a control region. The A + T content is 69.93% for the majority strand. All protein-coding genes start with the typical ATN initiation codons except for cox1, which uses AAT. Phylogenetic analyses based on ML and Bayesian methods shown consistent topologies among Lucanidae.

Characterization of the Complete Mitochondrial Genome of Harpalus sinicus and Its Implications for Phylogenetic Analyses

In this study, we report the complete mitochondrial genome of Harpalus sinicus (occasionally named as the Chinese ground beetle) which is the first mitochondrial genome for Harpalus. The mitogenome is 16,521 bp in length, comprising 37 genes, and a control region. The A + T content of the mitogenome is as high as 80.6%. A mitochondrial origins of light-strand replication (OL)-like region is found firstly in the insect mitogenome, which can form a stem-loop hairpin structure. Thirteen protein-coding genes (PCGs) share high homology, and all of them are under purifying selection. All tRNA genes (tRNAs) can be folded into the classic cloverleaf secondary structures except tRNA-Ser (GCU), which lacks a dihydrouridine (DHU) stem. The secondary structure of two ribosomal RNA genes (rRNAs) is predicted based on previous insect models. Twelve types of tandem repeats and two stem-loop structures are detected in the control region, and two stem-loop structures may be involved in the initiation of replication and transcription. Additionally, phylogenetic analyses based on mitogenomes suggest that Harpalus is an independent lineage in Carabidae, and is closely related to four genera (Abax, Amara, Stomis, and Pterostichus). In general, this study provides meaningful genetic information for Harpalus sinicus and new insights into the phylogenetic relationships within the Carabidae.