The mitochondrial genome of a social wasp, Vespula rufa (Linnaeus, 1758) (Hymenoptera: Vespidae)

For an ongoing systematic study of Korean Vespidae, we analyzed the complete mitochondrial genome of a social wasp, Vespula rufa (Linnaeus 1758), from the South Korea. The mitogenome is 17,521 bp in length, comprising 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and two ribosomal RNA (rRNA) genes. The nucleotide composition is 40.5% adenines, 43.0% thymines, 6.0% guanines, and 10.5% cytosines. The GC content is 16.5%. A maximum-likelihood analysis was conducted to understand phylogenetic relationships, based on 13 complete mitogenome sequences of Vespinae species. We recognized that V. rufa is being placed basal within the genus Vespula. The complete mitochondrial genome of V. rufa provides useful genetic information for further studies.

De novo assembly of the complete mitochondrial genome of Mycetophylax simplex Emery, 1888 through organelle targeting revels no substantial expansion of gene spacers, but rather some slightly shorter genes

Mitochondria play a key role in cell biology and have their own genome, residing in a highly oxidative environment that induces faster changes than the nuclear genome. Because of this, mitochondrial markers have been exploited to reconstruct phylogenetic and phylogeographic relationships in studies of adaptation and molecular evolution. In this study, we determined the complete mitogenome of the fungus-farming ant Mycetophylax simplex (Hymenoptera, Formicidae) and conducted a comparative analysis among 29 myrmicine ant mitogenomes. Mycetophylax simplex is an endemic ant that inhabits sand dunes along the southern Atlantic coast. Specifically, the species occur in the ecosystem known as "restinga", within the Atlantic Forest biome. Due to habitat degradation, land use and decline of restinga habitats, the species is considered locally extinct in extremely urban beaches and is listed as vulnerable on the Brazilian Red List (ICMBio). We employed a mitochondrion-targeting approach to obtain the complete mitogenome through high-throughput DNA sequencing technology. This method allowed us to determine the mitogenome with high performance, coverage and low cost. The circular mitogenome has a length of 16,367 base pairs enclosing 37 genes (13 protein-coding genes, 22 tRNAs and 2 rRNAs) along with one control region (CR). All the protein-coding genes begin with a typical ATN codon and end with the canonical stop codons. All tRNAs formed the fully paired acceptor stems and fold into the typical cloverleaf-shaped secondary structures. The gene order is consistent with the shared Myrmicinae structure, and the A + T content of the majority strand is 81.51%. Long intergenic spacers were not found but some gene are slightly shorter. The phylogenetic relationships based on concatenated nucleotide and amino acid sequences of the 13 protein-coding genes, using Maximum Likelihood and Bayesian Inference methods, indicated that mitogenome sequences were useful in resolving higher-level relationship within Formicidae.

Characterization of Complete Mitochondrial Genome and Phylogenetic Analysis of a Nocturnal Wasps-Provespa barthelemyi (Hymenoptera: Vespidae)

Genus Provespa contains nocturnal wasps mainly found in the southeastern region of Asia. There are no complete genome resources available of this genus, which hinders the study of its phylogenetic evolution and the origin of nocturnal behavior in the Vespidae family. Through high-throughput sequencing, we obtained the mitochondrial genome of Provespa barthelemyi (Hymenoptera: Vespidae), which is 17,721 base pairs in length and contains 13 protein-coding genes (PCGs), 22 tRNAs, and two rRNAs. We identified four gene rearrangement events of P. barthelemyi that frequently occur in the Vespidae family. We used Maximum Likelihood (ML) methodologies to construct a phylogenetic tree based on the sequenced mitochondrial genome and the available data of reported species belonging to Vespinae. Our findings confirmed the monophyly of Vespinae. Our study reports the first complete mitochondrial genome of Provespa and compares its characteristics with other mitochondrial genomes in the family Vespidae. This research should shed light on the phylogenetic relationships and ecological characteristics of the Vespidae family.

Novel Gene Rearrangement Pattern in Pachycrepoideus vindemmiae Mitochondrial Genome: New Gene Order in Pteromalidae (Hymenoptera: Chalcidoidea)

The mitochondrial genomes of Muscidifurax similadanacus, M. sinesensilla, Nasonia vitripennis, and Pachycrepoideus vindemmiae were sequenced to better understand the structural evolution of Pteromalidae mitogenomes. These newly sequenced mitogenomes all contained 37 genes. Nucleotide composition was AT-biased and the majority of the protein-coding genes exhibited a negative AT skew. All 13 protein-coding genes (PCGs) initiated with the standard start codon of ATN, excepted for nad1 of N. vitripennis, which started with TTG, and terminated with a typical stop codon TAA/TAG or an incomplete stop codon T. All transfer RNA (tRNA) genes were predicted to fold into the typical clover-leaf secondary structures, except for trnS1, which lacks the DHU arm in all species. In P. vindemmiae, trnR and trnQ lack the DHU arm and TΨC arm, respectively. Although most genes evolved under a strong purifying selection, the Ka/Ks value of the atp8 gene of P. vindemmiae was greater than 1, indicating putative positive selection. A novel transposition of trnR in P. vindemmiae was revealed, which was the first of this kind to be reported in Pteromalidae. Two kinds of datasets (PCG12 and AA) and two inference methods (maximum likelihood and Bayesian inference) were used to reconstruct a phylogenetic hypothesis for the newly sequenced mitogenomes of Pteromalidae and those deposited in GenBank. The topologies obtained recovered the monophyly of the three subfamilies included. Pachyneurinae and Pteromalinae were recovered as sister families, and both appeared sister to Sycophaginae. The pairwise breakpoint distances of mitogenome rearrangements were estimated to infer phylogeny among pteromalid species. The topology obtained was not totally congruent with those reconstructed using the ML and BI methods.

New Dielis species and structural dichotomy of the mitochondrial cox2 gene in Scoliidae wasps

Some mitochondrial protein-coding genes of protists and land plants have split over the course of evolution into complementary genes whose products can form heteromeric complexes that likely substitute for the undivided proteins. One of these genes, cox2, has also been found to have split in animals, specifically in Scoliidae wasps (Hymenoptera: Apocrita) of the genus Dielis (Campsomerini), while maintaining the conventional structure in related Scolia (Scoliini). Here, a hitherto unrecognized Nearctic species of Dielis, D. tejensis, is described based on its phenotype and mtDNA. The mitogenome of D. tejensis sp. nov. differs from that of the sympatric sibling species Dielis plumipes fossulana by the reduced size of the cox2-dividing insert, which, however, still constitutes the fifth part of the mtDNA; an enlarged nad2-trnW intergenic region; the presence of two trnK^ttt paralogues; and other features. Both species of Dielis have a unique insertion of a threonine in COXIIA, predicted to be involved in COXIIA-COXIIB docking, and substitutions of two hydrophobic residues with redox-active cysteines around the Cu_A centre in COXIIB. Importantly, the analysis of mtDNA from another Campsomerini genus, Megacampsomeris, shows that its cox2 gene is also split. The presented data highlight evolutionary processes taking place in hymenopteran mitogenomes that do not fall within the mainstream of animal mitochondrion evolution.

The Mitochondrial Genomes of Two Parasitoid Wasps Protapanteles immunis and Parapanteles hyposidrae (Hymenoptera: Braconidae) with Phylogenetic Implications and Novel Gene Rearrangements

Parapanteles hypsidrae (Wilkinson, 1928) and Protapanteles immunis (Haliday, 1834) are the most important parasitic wasps of Ectropis grisescens Warren and Ectropis obliqua (Prout). We sequenced and annotated the mitochondrial genomes of Pa. hyposidrae and Pr. immunis, which are 17,063 bp and 16,397 bp in length, respectively, and possess 37 mitochondrial genes. We discovered two novel types of gene rearrangement, the local inversion of nad4L in Pa. hyposidrae and the remote inversion of the block cox3-nad3-nad5-nad4 in Pr. immunis, within the mitogenomes of Braconidae. The phylogenetic analysis supported the subfamily Microgastrinae is a monophyletic group, but the tribes Apantelini and Cotesiini within this subfamily are paraphyletic groups.

Mitochondrial genome rearrangements and phylogenomics of the Hymenoptera (Insecta) using an expanded taxon sample

The order Hymenoptera is one of the most species-rich insect orders, with more than 150,000 described extant species. Many hymenopteran insects have very different mitochondrial genome (mitogenome) organizations compared to the putative ancestral organization of insects. In this study, we sequenced 18 mitogenomes of representatives in the order Hymenoptera to increase taxonomic sampling. A total of 475 species were used in phylogenetic analyses, including 18 new mitogenomes and 457 existing mitogenomes. Using a site-heterogeneous model, Bayesian's inference from amino acid data yielded more resolved relationships among Hymenoptera than maximum-likelihood analysis and coalescent-based species analyses. The monophyly of Symphyta was not supported. The Xyeloidea was the earliest branching clade in the Hymenoptera. The Orussoidea was closely related to Apocrita. Within Apocrita, the Parasitoida was non-monophyletic. The monophyly of most Parasitoida superfamilies received strong support. The Proctotrupomorpha clade was supported in Bayesian's analysis. The Apoidea was monophyletic when excluding Ampulex compressa from consideration. The superfamilies Vespoidea and Chrysidoidea were found to be non-monophyletic. Comparisons of mitochondrial gene order revealed a higher frequency of gene rearrangement among lineages with a parasitoid lifestyle, particularly prominent in Chalcidoidea. The degree of gene rearrangement ranked second in specific taxa of Cynipoidea and Ichneumonoidea.

The mitochondrial genome of Chelonus formosanus (Hymenoptera: Braconidae) with novel gene orders and phylogenetic implications

Chelonus formosanus Sonan is an important egg-larval parasitoid of noctuid moths and a potential candidate for understanding interactions between host and parasitoid mediated by polydnavirues (PDVs). We sequenced and annotated the mitochondrial genome of C. formosanus, which is 15,466 bp in length and possesses 38 mitochondrial genes. However, unlike most animal mitochondrial genomes, it contains one extra trnF gene. There are five transfer RNA (tRNA) rearrangement events compared with the ancestral gene order, which is a novel rearrangement type in Hymenoptera for all published mitogenomes so far. Phylogenetic trees supported C. formosanus from the subfamily Cheloninae was closely related to the subfamily Cardiochilinae and Microgastrinae.

First report of complete mitochondrial genome in the subfamily Alleculinae and mitochondrial genome-based phylogenetics in Tenebrionidae (Coleoptera: Tenebrionoidea)

Despite worldwide distribution and rich diversity, the knowledge of mitochondrial genome (mtgenome) characteristics within the family Tenebrionidae (Coleoptera) is still very limited, and phylogenetics remain unresolved for the family. In this study, the mtgenomes for 12 species are newly sequenced and annotated. Three of the species represent the first report of complete mtgenomes from the subfamily Alleculinae. Comparative analysis of 36 tenebrionid mtgenomes finds that gene composition and order are the same as a theoretical ancestral insect mtgenome, and AT bias, length variation, and codon usage are consistent with other reported beetle mtgenomes. Some intergenic overlap and gap sequences may contain phylogenetically informative information, whereas neither the conserved nor repeat sequences in the control region do. The subfamilies Lagriinae, Pimeliinae, Stenochiinae, and Alleculinae are found to be monophyletic, and the subfamilies Diaperinae and Tenebrioninae polyphyletic in our study. Furthermore, Lagriinae is sister to the rest of the subfamilies. At the tribal level, the tribes Lagriini in the subfamily Lagrrinae; Strongyliini in Stenochiinae; Cteniopodini in Alleculinae; and Triboliini, Opatrini, and Amarygmini in Tenebrioninae are monophyletic; while Diaperini in Diaperinae, and Tenebrionini in Tenebrioninae are polyphyletic.

First Comprehensive Analysis of Both Mitochondrial Characteristics and Mitogenome-Based Phylogenetics in the Subfamily Eumeninae (Hymenoptera: Vespidae)

The subfamily Eumeninae plays a significant role in the biological control of agricultural pests. However, the characteristics of eumenine mitogenomes that are important molecular markers for phylogenetics are not clearly revealed. Here, 52 eumenine mitogenomes are newly sequenced and annotated, and the phylogenetic relationships of the subfamily are comprehensively analyzed based on 87 vespid mitogenomes. Through the comparative analysis of the 54 eumenine mitogenomes, the gene compositions of about one half of the 54 species match with ancestral insect mitogenome, and remaining others contain two trnM which are highly similar, with 51.86% (Eumenes tripunctatus) to 90.65% (Pseumenes nigripectus) sequence identities, which is unique among the reported mitogenomes of the family Vespidae. Moreover, the translocation trnL1 upstream of nad1 is a common rearrangement event in all eumenine mitogenomes. The results of phylogenetic analyses support the paraphyly of the subfamily Eumeninae and the tribe Odynerini, respectively, and the monophyly of the tribe Eumenini, and verify that the tribe Zethini is a valid subfamily Zethinae. In this study, the relationships between some genera such as Allorhynchium and Pararrhynchium or the taxonomic status of the subgenera such as Eremodynerus and Dirhynchium are found to be confusing and there should be further inquiry with more samples.

Mitochondrial composition of and diffusion limiting factors of three social wasp genera Polistes, Ropalidia, and parapolybia (Hymenoptera: Vespidae)

Background

Social wasps Polistes, Ropalidia, and Parapolybia, belonging to the subfamily Polistinae, have obviously different distribution patterns, yet the factors leading to this difference remain unknown.

Results

The 17 newly sequenced mitogenomes of Polistes, Ropalidia, and Parapolybia contain 37 genes, and there are obvious differences among the compositions of the three genera. The monophyly of the genus Polistes and a monophyletic Ropalidiini: (Ropalidia + Parapolybia) are concordant with previous morphological analysis of the subfamily Polistinae. Our inferred divergence time demonstrates Polistes (at around 69 Ma) was diverged earlier than Ropalidia and Parapolybia (at around 61 Ma). The rearrangement of both trnY and trnL1 are shared by all the Polistinae. In addition, the unique rearrangement of TDRL derived at 69 Ma is detected in Polistes, and Ropalidia contains a Reversal which may derive at 61 Ma. Hereafter, the possibility is elaborated that Polistes originated in Aisa and then dispersed from Africa to South America, and Polistes and Ropalidia spread from Southeast Asia to Australia. At last, continental drift and Quaternary Ice Ages are inferred to be two main limiting factors in the current distributions of the three genera.

Conclusions

Obvious differences occur in the mitochondrial composition of Polistes, Ropalidia, and Parapolybia. According to the reconstructed time-calibrated framework, it is inquired that the continental drifts and the climate are mainly diffusion limiting factors of the three genera.

Supplementary information

The online version contains supplementary material available at 10.1186/s12862-022-02017-6.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea)

Simple Summary

Papilionoidea is a superfamily of Lepidoptera encompassing about 19,000 species. In the present work, we study the evolution of the structure of the mitogenome of these lepidopterans. The mechanisms generating the eight arrangements known for Papilionoidea were investigated analysing the movements of different mitochondrial genes. Five newly sequenced/assembled mitogenomes were included in our analysis involving more than 600 genomes. We provide new findings that help to understand the evolution of the gene orders MIQGO, IMQGO, 2S1GO, ES1GO and S1NGO in different butterflies. We demonstrate that the evolution of the 2S1GO in Lycaenidae followed a complicated pathway with multiple events of duplication and loss of trnS1 and changes in anticodon. We describe two new gene orders 2FFGO and 4QGO for Ampittia subvittatus (Hesperiidae) and Bhutanitis thaidina (Papilionidae).

Abstract

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

Novel gene rearrangement in the mitochondrial genome of Anastatus fulloi (Hymenoptera Chalcidoidea) and phylogenetic implications for Chalcidoidea

The genus Anastatus comprises a large group of parasitoids, including several biological control agents in agricultural and forest systems. The taxonomy and phylogeny of these species remain controversial. In this study, the mitogenome of A. fulloi Sheng and Wang was sequenced and characterized. The nearly full-length mitogenome of A. fulloi was 15,692 bp, compromising 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and a control region (CR). The total A + T contents were 83.83%, 82.18%, 87.58%, 87.27%, and 82.13% in the whole mitogenome, 13 PCGs, 22 tRNA genes, 2 rRNA genes, and CR, respectively. The mitogenome presented negative AT skews and positive GC skews, except for the CR. Most PCGs were encoded on the heavy strand, started with ATN codons, and ended with TAA codons. Among the 3736 amino acid-encoding codons, TTA (Leu1), CGA (Arg), TCA (Ser2), and TCT (Ser2) were predominant. Most tRNAs had cloverleaf secondary structures, except trnS1, with the absence of a dihydrouridine (DHU) arm. Compared with mitogenomes of the ancestral insect and another parasitoid within Eupelmidae, large-scale rearrangements were found in the mitogenome of A. fulloi, especially inversions and inverse transpositions of tRNA genes. The gene arrangements of parasitoid mitogenomes within Chalcidoidea were variable. A novel gene arrangement was presented in the mitogenome of A. fulloi. Phylogenetic analyses based on the 13 protein-coding genes of 20 parasitoids indicated that the phylogenetic relationship of 6 superfamilies could be presented as Mymaridae + (Eupelmidae + (Encyrtidae + (Trichogrammatidae + (Pteromalidae + Eulophidae)))). This study presents the first mitogenome of the Anastatus genus and offers insights into the identification, taxonomy, and phylogeny of these parasitoids.

Insight into the Phylogenetic Relationships among Three Subfamilies within Heptageniidae (Insecta: Ephemeroptera) along with Low-Temperature Selection Pressure Analyses Using Mitogenomes

We determined 15 complete and two nearly complete mitogenomes of Heptageniidae belonging to three subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae) and six genera (Afronurus, Epeorus, Leucrocuta, Maccaffertium, Stenacron, and Stenonema). Species of Rhithrogeninae and Ecdyonurinae had the same gene rearrangement of CR-I-M-Q-M-ND2, whereas a novel gene rearrangement of CR-I-M-Q-NCR-ND2 was found in Heptageniinae. Non-coding regions (NCRs) of 25-47 bp located between trnA and trnR were observed in all mayflies of Heptageniidae, which may be a synapomorphy for Heptageniidae. Both the BI and ML phylogenetic analyses supported the monophyly of Heptageniidae and its subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae). The phylogenetic results combined with gene rearrangements and NCR locations confirmed the relationship of the subfamilies as (Heptageniinae + (Rhithrogeninae + Ecdyonurinae)). To assess the effects of low-temperature stress on Heptageniidae species from Ottawa, Canada, we found 27 positive selection sites in eight protein-coding genes (PCGs) using the branch-site model. The selection pressure analyses suggested that mitochondrial PCGs underwent positive selection to meet the energy requirements under low-temperature stress.

Novel tRNA gene rearrangements in the mitochondrial genomes of praying mantises (Mantodea: Mantidae): Translocation, duplication and pseudogenization

Gene rearrangements have been found in several mitochondrial genomes of Mantodea, located in the gene blocks CR-I-Q-M-ND2, COX1-K-D-ATP8 and ND3-A-R-N-S-E-F-ND5. We have sequenced one mitogenome of Amelidae (Yersinia mexicana) and six mitogenomes of Mantidae to discuss the mitochondrial gene rearrangement and the phylogenetic relationship within Mantidae. These mitogenomes showed rearrangements of tRNA genes except for Asiadodis yunnanensis and Hierodula zhangi. These novel gene rearrangements of Mantidae were primarily concentrated in the region of CR-I-Q-M-ND2, including gene translocation, duplication and pseudogenization. For the occurrences of these rearrangements, the tandem duplication-random loss (TDRL) model and slipped-strand mispairing model were suitable to explain. Large non-coding regions (LNCRs) located in the region of CR-I-Q-M-ND2 were detected in most Mantidae species, whereas some LNCRs had high similarity to the control region (CR). Both BI and ML phylogenetic analyses supported the monophyly of Mantidae and the paraphyly of Mantinae. The phylogenetic results with the gene order and the location of NCRs acted as forceful evidence that specific gene rearrangements and special LNCRs may be synapomorphies for several groups of mantises.

New Mitochondrial Gene Rearrangement in Psyttalia concolor, P. humilis and P. lounsburyi (Hymenoptera: Braconidae), Three Parasitoid Species of Economic Interest

The family Braconidae consists mostly of specialized parasitoids, some of which hold potential in biocontrol of agricultural pests. Psyttalia concolor, Psyttalia humilis and Psyttalia lounsburyi are parasitoids associated with Bactrocera oleae, a major pest of cultivated olives. The native range of Psyttalia concolor is the Mediterranean, and P. humilis and P. lounsburyi are native to sub-Saharan Africa. This study reports the mitochondrial genomes of the three species, thus laying the foundation for mitogenomic analyses in the genus Psyttalia. Comparative mitogenomics within Braconidae showed a novel gene arrangement in Psyttalia in involving translocation and inversion of transfer RNA genes. The placement of Psyttalia in the subfamily Opiinae was well-supported, and the divergence between Psyttalia and its closest relative (Diachasmimorpha longicaudata) was at ~55 MYA [95% highest posterior density (HPD): 34-83 MYA]. Psyttalia lounsburyi occupied the most basal position among the three Psyttalia, having diverged from the other two species ~11 MYA (95% HPD: 6-17 MYA). Psyttalia concolor and P. humilis were recovered as sister species diverged at ~2 MYA (95% HPD: 1.1-3.6 MYA). This phylogeny combining new sequences and a set of 31 other cyclostomes and non-cyclostomes highlights the importance of a comprehensive taxonomic coverage of Braconidae mitogenomes to overcome the lack of robustness in the placement of several subfamilies.

Six complete mitochondrial genomes of mayflies from three genera of Ephemerellidae (Insecta: Ephemeroptera) with inversion and translocation of trnI rearrangement and their phylogenetic relationships

As a small order of Pterygota (Insecta), Ephemeroptera has almost 3,500 species around the world. Ephemerellidae is a widely distributed common group of Ephemeroptera. However, the relationship among Ephemerellidae, Vietnamellidae and Teloganellidae is still in dispute. In this study, we sequenced six complete mitogenomes of three genera from Ephemerellidae (Insecta: Ephemeroptera): Ephemerella sp. Yunnan-2018, Serratella zapekinae, Serratella sp. Yunnan-2018, Serratella sp. Liaoning-2019, Torleya grandipennis and T. tumiforceps. These mitogenomes were employed to reveal controversial phylogenetic relationships among the Ephemeroptera, with emphasis on the phylogenetic relationships among Ephemerellidae. The lengths of the six mayfly mitogenomes ranged from 15,134 bp to 15,703 bp. Four mitogenomes of Ephemerella sp. Yunnan-2018, Serratella zapekinae, Serratella sp. Yunnan-2018 and Serratella sp. Liaoning-2019 had 22 tRNAs including an inversion and translocation of trnI. By contrast, the mitogenomes of T. tumiforceps and T. grandipennis had 24 tRNAs due to an extra two copies of inversion and translocation of trnI. Within the family Ephemerellidae, disparate gene rearrangement occurred in the mitogenomes of different genera: one copy of inversion and translocation trnI in the genera Ephemerella and Serratella, and three repeat copies of inversion and translocation of trnI in the genus Torleya. A large non-coding region (≥200 bp) between trnS1 (AGN) and trnE was detected in T. grandipennis and T. tumiforceps. Among the phylogenetic relationship of the Ephemeroptera, the monophyly of almost all families except Siphlonuridae was supported by BI and ML analyses. The phylogenetic results indicated that Ephemerellidae was the sister clade to Vietnamellidae whereas Teloganellidae was not a sister clade of Ephemerellidae and Vietnamellidae.

Rearrangement and evolution of mitochondrial genomes in Thysanoptera (Insecta)

Prior to this study, complete mitochondrial genomes from Order Thysanoptera were restricted to a single family, the Thripidae, resulting in a biased view of their evolution. Here we present the sequences for the mitochondrial genomes of four additional thrips species, adding three extra families and an additional subfamily, thus greatly improving taxonomic coverage. Thrips mitochondrial genomes are marked by high rates of gene rearrangement, duplications of the control region and tRNA mutations. Derived features of mitochondrial tRNAs in thrips include gene duplications, anticodon mutations, loss of secondary structures and high gene translocation rates. Duplicated control regions are found in the Aeolothripidae and the 'core' Thripinae clade but do not appear to promote gene rearrangement as previously proposed. Phylogenetic analysis of thrips mitochondrial sequence data supports the monophyly of two suborders, a sister-group relationship between Stenurothripidae and Thripidae, and suggests a novel set of relationships between thripid genera. Ancestral state reconstructions indicate that genome rearrangements are common, with just eight gene blocks conserved between any thrips species and the ancestral insect mitochondrial genome. Conversely, 71 derived rearrangements are shared between at least two species, and 24 of these are unambiguous synapomorphies for clades identified by phylogenetic analysis. While the reconstructed sequence of genome rearrangements among the protein-coding and ribosomal RNA genes could be inferred across the phylogeny, direct inference of phylogeny from rearrangement data in MLGO resulted in a highly discordant set of relationships inconsistent with both sequence-based phylogenies and previous morphological analysis. Given the demonstrated rates of genomic evolution within thrips, extensive sampling is needed to fully understand these phenomena across the order.

The first divergence time estimation of the subfamily Stenogastrinae (Hymenoptera: Vespidae) based on mitochondrial phylogenomics

In this study, the mitochondrial genomes of three Stenogastrinae species, Eustenogaster scitula, Liostenogaster nitidipennis and Parishnogaster mellyi were sequenced and annotated, and a total of 16 vespid mtgenomes are comparatively analyzed. Our results indicate that codon usage bias is mainly dominated by mutational pressure, and affected only slightly by natural selection. Selective pressure analysis of protein-coding genes (PCGs) shows that the highest evolutionary rate is present in NADH complex I, and the lowest in cox1. Compared with the reported mtgenomes of other Vespidae, in Stenogastrinae, trnH is shifted to a new position. Phylogenetic analyses are performed using Bayesian method and Maximum Parsimony. Phylogenetic analysis further confirms that the Stenogastrinae is the sister group of all remaining Vespidae. Divergence time of Stenogastrinae from other Vespidae is estimated at ~ 166 Mya. Our results also support that eusociality evolved twice in the family Vespidae.

The complete mitochondrial genome of Orancistrocerusaterrimusaterrimus and comparative analysis in the family Vespidae (Hymenoptera, Vespidae, Eumeninae)

To date, only one mitochondrial genome (mitogenome) in the Eumeninae has been reported in the world and this is the first report in China. The mitogenome of O.a.aterrimus is 17 972 bp long, and contains 38 genes, including 13 protein coding genes (PCGs), 23 tRNA genes, two rRNA genes, a long non-coding region (NCR), and a control region (CR). The mitogenome has 79.43% A + T content, its 13 PCGs use ATN as the initiation codon except for cox1 using TTG, and nine genes used complete translation termination TAA and four genes have incomplete stop codon T (cox2, cox3, nad4, and cytb). Twenty-two of 23 tRNAs can form the typical cloverleaf secondary structure except for trnS1. The CR is 1 078 bp long with 84.69% A+T content, comprising 28 bp tandem repeat sequences and 13 bp T-strech. There are two gene rearrangements which are an extra trnM2 located between trnQ and nad2 and the trnL2 in the upstream of nad1. Within all rearrangements of these mitogenomes reported in the family Vespidae, the translocation between trnS1 and trnE genes only appears in Vespinae, and the translocation of trnY in Polistinae and Vespinae. The absent codons of 13 PCGs in Polistinae are more than those both in Vespinae and Eumeninae in the family Vespidae. The study reports the complete mitogenome of O.a.aterrimus, compares the characteristics and construct phylogenetic relationships of the mitogenomes in the family Vespidae.

Extensive gene rearrangements in the mitochondrial genomes of two egg parasitoids, Trichogramma japonicum and Trichogramma ostriniae (Hymenoptera: Chalcidoidea: Trichogrammatidae)

Animal mitochondrial genomes usually exhibit conserved gene arrangement across major lineages, while those in the Hymenoptera are known to possess frequent rearrangements, as are those of several other orders of insects. Here, we sequenced two complete mitochondrial genomes of Trichogramma japonicum and Trichogramma ostriniae (Hymenoptera: Chalcidoidea: Trichogrammatidae). In total, 37 mitochondrial genes were identified in both species. The same gene arrangement pattern was found in the two species, with extensive gene rearrangement compared with the ancestral insect mitochondrial genome. Most tRNA genes and all protein-coding genes were encoded on the minority strand. In total, 15 tRNA genes and seven protein-coding genes were rearranged. The rearrangements of cox1 and nad2 as well as most tRNA genes were novel. Phylogenetic analysis based on nucleotide sequences of protein-coding genes and on gene arrangement patterns produced identical topologies that support the relationship of (Agaonidae + Pteromalidae) + Trichogrammatidae in Chalcidoidea. CREx analysis revealed eight rearrangement operations occurred from presumed ancestral gene order of Chalcidoidea to form the derived gene order of Trichogramma. Our study shows that gene rearrangement information in Chalcidoidea can potentially contribute to the phylogeny of Chalcidoidea when more mitochondrial genome sequences are available.

Gene characteristics of the complete mitochondrial genomes of Paratoxodera polyacantha and Toxodera hauseri (Mantodea: Toxoderidae)

The family Toxoderidae (Mantodea) contains an ecologically diverse group of praying mantis species that have in common greatly elongated bodies. In this study, we sequenced and compared the complete mitochondrial genomes of two Toxoderidae species, Paratoxodera polyacantha and Toxodera hauseri, and compared their mitochondrial genome characteristics with another member of the Toxoderidae, Stenotoxodera porioni (KY689118). The lengths of the mitogenomes of T. hauseri and P. polyacantha were 15,616 bp and 15,999 bp, respectively, which is similar to that of S. porioni (15,846 bp). The size of each gene as well as the A+T-rich region and the A+T content of the whole genome were also very similar among the three species as were the protein-coding genes, the A+T content and the codon usages. The mitogenome of T. hauseri had the typical 22 tRNAs, whereas that of P. polyacantha had 26 tRNAs including an extra two copies of trnA-trnR. Intergenic regions of 67 bp and 76 bp were found in T. hauseri and P. polyacantha, respectively, between COX2 and trnK; these can be explained as residues of a tandem duplication/random loss of trnK and trnD. This non-coding region may be synapomorphic for Toxoderidae. In BI and ML analyses, the monophyly of Toxoderidae was supported and P. polyacantha was the sister clade to T. hauseri and S. porioni.

Complete mitochondrial DNA sequence of the alien hornet Vespa velutina (Insecta: Hymenoptera) invading Kyushu Island, Japan

We analyzed the complete mitochondrial genome of the invasive Asian hornet Vespa velutina from Kyushu Island, Japan. The mitochondrial genome of V. velutina was identified as a circular molecule of 16,388 bp. We predicted that the genome contains 13 protein-coding genes (PCGs), 22 tRNA genes, and 2 rRNA genes, along with one A + T-rich control region. The average AT content is 81.68%. Molecular phylogenetic analysis using the 13 mitochondrial PCGs from 11 closely related taxa of Vespidae indicated that the V. velutina invading the Japanese Islands of Kyushu and Tsushima have a common origin.

Complete mitochondrial DNA sequence of the tropical hornet Vespa affinis (Insecta, Hymenoptera)

We analyzed the complete mitochondrial genome of the Asian tropical hornet Vespa affinis from Ishigaki Island, Japan. It consisted of a circular molecule with 19,109 bp, which is larger to other hornet species e.g. V. velutina. We predicted the genome contained 13 protein-coding, 22 tRNA, and two rRNA genes, along with one A + T-rich control region. The repetitive sequences were confirmed at multiple positions in the non-coding genes. The initiation codons ATA was found in one, ATG in seven, and ATT in five genes, while the termination codons TAA and TAG were observed 11 and two genes, respectively. The average AT content of the genome was 78.4%.

The Evolution of Dark Matter in the Mitogenome of Seed Beetles

Animal mitogenomes are generally thought of as being economic and optimized for rapid replication and transcription. We use long-read sequencing technology to assemble the remarkable mitogenomes of four species of seed beetles. These are the largest circular mitogenomes ever assembled in insects, ranging from 24,496 to 26,613 bp in total length, and are exceptional in that some 40% consists of non-coding DNA. The size expansion is due to two very long intergenic spacers (LIGSs), rich in tandem repeats. The two LIGSs are present in all species but vary greatly in length (114-10,408 bp), show very low sequence similarity, divergent tandem repeat motifs, a very high AT content and concerted length evolution. The LIGSs have been retained for at least some 45 my but must have undergone repeated reductions and expansions, despite strong purifying selection on protein coding mtDNA genes. The LIGSs are located in two intergenic sites where a few recent studies of insects have also reported shorter LIGSs (>200 bp). These sites may represent spaces that tolerate neutral repeat array expansions or, alternatively, the LIGSs may function to allow a more economic translational machinery. Mitochondrial respiration in adult seed beetles is based almost exclusively on fatty acids, which reduces the need for building complex I of the oxidative phosphorylation pathway (NADH dehydrogenase). One possibility is thus that the LIGSs may allow depressed transcription of NAD genes. RNA sequencing showed that LIGSs are partly transcribed and transcriptional profiling suggested that all seven mtDNA NAD genes indeed show low levels of transcription and co-regulation of transcription across sexes and tissues.

Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes

To better understand the phylogeny and evolution of mosquitoes, the complete mitochondrial genome (mitogenome) of Anopheles stephensi and An. dirus were sequenced and annotated, and a total of 50 mosquito mitogenomes were comparatively analyzed. The complete mitogenome of An. stephensi and An. dirus is 1,5371 bp and 1,5406 bp long, respectively. The main features of the 50 mosquito mitogenomes are conservative: 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, positive AT-skew and negative GC-skew. The gene order trnA-trnR in ancestral insects is rearranged. All tRNA genes have the typical clover leaf secondary structure but tRNA Ser . The control regions are highly variable in size. PCGs show signals of purifying selection, but evidence for positive selection in ND2, ND4 and ND6 is found. Bayesian and Maximum Likelihood phylogenetic analyses based on all PCG nucleotides produce an identical tree topology and strongly support the monophyly of subgenera Cellia, Anopheles, Keterszia and Nyssorhynchus, the sister relationship of the subgenera Nyssorhynchus and Keterszia, and Cellia and Anopheles. The most recent ancestor of the genus Anopheles and Culicini + Aedini exited ~145 Mya ago. This is the first comprehensive study of mosquito mitogenomes, which are effective for mosquito phylogeny at various taxonomic levels.

Novel insights into mitochondrial gene rearrangement in thrips (Insecta: Thysanoptera) from the grass thrips, Anaphothrips obscurus

We sequenced the mitochondrial (mt) genome of the grass thrips, Anaphothrips obscurus, which is highly rearranged and differs from the four thrips species reported previously in the arrangement of both tRNA genes and a protein-coding gene, nad3, and in the copy number of the control region (CR). We reconstructed the phylogeny of the thrips with mt genome sequences, and used it as a framework to gain insights into mt genome evolution in thrips. It is evident that A. obscurus is less rearranged in mt genome organization than the other four known thrips. nad3 is in its ancestral location in A. obscurus but was translocated in other four thrips. Also, A. obscurus has one CR, which is ancestral to hexapods whereas other thrips have two or three CRs. All of the five thrips whose mt genomes have been sequenced to date are from the subfamily Thripinae, which represents about a quarter of the species richness in the order Thysanoptera. The high variation in mt genome organization observed in a subfamily challenges our knowledge about animal mt genomes. It remains to be investigated why mt genomes evolved so fast in the subfamily Thripinae and how mt genomes evolved in other lineages of thrips.

Nearly complete mitogenome of hairy sawfly, Corynis lateralis (Brullé, 1832) (Hymenoptera: Cimbicidae): rearrangements in the IQM and ARNS1EF gene clusters

The Cimbicidae is a small family of the primitive and relatively less diverse suborder Symphyta (Hymenoptera). Here, nearly complete mitochondrial genome (mitogenome) of hairy sawfly, Corynis lateralis (Hymenoptera: Cimbicidae) was sequenced using next generation sequencing and comparatively analysed with the mitogenome of Trichiosoma anthracinum. The sequenced length of C. lateralis mitogenome was 14,899 bp with an A+T content of 80.60%. All protein coding genes (PCGs) are initiated by ATN codons and all are terminated with TAR or T- stop codon. All tRNA genes preferred usual anticodons. Compared with the inferred insect ancestral mitogenome, two tRNA rearrangements were observed in the IQM and ARNS1EF gene clusters, representing a new event not previously reported in Symphyta. An illicit priming of replication and/or intra/inter-mitochondrial recombination and TDRL seem to be responsible mechanisms for the rearrangement events in these gene clusters. Phylogenetic analyses confirmed the position of Corynis within Cimbicidae and recovered a relationship of Tenthredinoidea + (Cephoidea + Orussoidea) in Symphyta.

Complete mitochondrial DNA sequence of the invasive hornet Vespa velutina (Insecta, Hymenoptera) found in Japan

In this study, we analyzed the complete mitochondrial genome of the invasive Asian hornet Vespa velutina from Japan. The mitochondrial genome of V. velutina was identified as a circular molecule of 16,765 bp, similar to that in other hornet species. It was predicted to contain 13 protein-coding, 20 tRNA, and two rRNA genes, along with one A + T-rich control region. The initiation codons ATC was found in one, ATG in four, ATT in five, and ATA in three genes, while TAA was the termination codon in all these genes. The average AT content of 13 protein-coding was 82%.

Next-Generation Sequencing of Two Mitochondrial Genomes from Family Pompilidae (Hymenoptera: Vespoidea) Reveal Novel Patterns of Gene Arrangement

Animal mitochondrial genomes have provided large and diverse datasets for evolutionary studies. Here, the first two representative mitochondrial genomes from the family Pompilidae (Hymenoptera: Vespoidea) were determined using next-generation sequencing. The sequenced region of these two mitochondrial genomes from the species Auplopus sp. and Agenioideus sp. was 16,746 bp long with an A + T content of 83.12% and 16,596 bp long with an A + T content of 78.64%, respectively. In both species, all of the 37 typical mitochondrial genes were determined. The secondary structure of tRNA genes and rRNA genes were predicted and compared with those of other insects. Atypical trnS1 using abnormal anticodons TCT and lacking D-stem pairings was identified. There were 49 helices belonging to six domains in rrnL and 30 helices belonging to three domains in rrns present. Compared with the ancestral organization, four and two tRNA genes were rearranged in mitochondrial genomes of Auplopus and Agenioideus, respectively. In both species, trnM was shuffled upstream of the trnI-trnQ-trnM cluster, and trnA was translocated from the cluster trnA-trnR-trnN-trnS1-trnE-trnF to the region between nad1 and trnL1, which is novel to the Vespoidea. In Auplopus, the tRNA cluster trnW-trnC-trnY was shuffled to trnW-trnY-trnC. Phylogenetic analysis within Vespoidea revealed that Pompilidae and Mutillidae formed a sister lineage, and then sistered Formicidae. The genomes presented in this study have enriched the knowledge base of molecular markers, which is valuable in respect to studies about the gene rearrangement mechanism, genomic evolutionary processes and phylogeny of Hymenoptera.

Multiple Lines of Evidence from Mitochondrial Genomes Resolve Phylogenetic Relationships of Parasitic Wasps in Braconidae

The rapid increase in the number of mitochondrial genomes in public databases provides opportunities for insect phylogenetic studies; but it also provides challenges because of gene rearrangements and variable substitution rates among both lineages and sites. Typically, phylogenetic studies use mitochondrial sequence data but exclude other features of the mitochondrial genome from analyses. Here, we undertook large-scale sequencing of mitochondrial genomes from a worldwide collection of specimens belonging to Braconidae, one of the largest families of Metazoa. The strand-asymmetry of base composition in the mitochondrial genomes of braconids is reversed, providing evidence for monophyly of the Braconidae. We have reconstructed a backbone phylogeny of the major lineages of Braconidae from gene order of the mitochondrial genomes. Standard phylogenetic analyses of DNA sequences provided strong support for both Cyclostomes and Noncyclostomes. Four subfamily complexes, that is, helconoid, euphoroid, sigalphoid, and microgastroid, within the Noncyclostomes were reconstructed robustly, the first three of which formed a monophyletic group sister to the last one. Aphidiinae was recovered as a lineage sister to other groups of Cyclostomes, while the Ichneutinae was recovered as paraphyletic. Separate analyses of the subdivided groups showed congruent relationships, employing different matrices and methods, for the internal nodes of the Cyclostomes and the microgastroid complex of subfamilies. This research, using multiple lines of evidence from mitochondrial genomes, illustrates multiple uses of mitochondrial genomes for phylogenetic inference in Braconidae.

Back to Water: Signature of Adaptive Evolution in Cetacean Mitochondrial tRNAs

The mitochondrion is the power plant of the eukaryotic cell, and tRNAs are the fundamental components of its translational machinery. In the present paper, the evolution of mitochondrial tRNAs was investigated in the Cetacea, a clade of Cetartiodactyla that retuned to water and thus had to adapt its metabolism to a different medium than that of its mainland ancestors. Our analysis focussed on identifying the factors that influenced the evolution of Cetacea tRNA double-helix elements, which play a pivotal role in the formation of the secondary and tertiary structures of each tRNA and consequently manipulate the whole translation machinery of the mitochondrion. Our analyses showed that the substitution pathways in the stems of different tRNAs were influenced by various factors, determining a molecular evolution that was unique to each of the 22 tRNAs. Our data suggested that the composition, AT-skew, and GC-skew of the tRNA stems were the main factors influencing the substitution process. In particular, the range of variation and the fluctuation of these parameters affected the fate of single tRNAs. Strong heterogeneity was observed among the different species of Cetacea. Finally, it appears that the evolution of mitochondrial tRNAs was also shaped by the environments in which the Cetacean taxa differentiated. This latter effect was particularly evident in toothed whales that either live in freshwater or are deep divers.

Complete mitochondrial genome of Camponotus atrox (Hymenoptera: Formicidae): a new tRNA arrangement in Hymenoptera

We sequenced the complete mitochondrial (mt) genome of Camponotus atrox (Hymenoptera: Formicidae), which is only distributed in Korea. The genome was 16 540 bp in size and contained typical sets of genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs). The C. atrox A+T-rich region, at 1402 bp, was the longest of all sequenced ant genomes and was composed of an identical tandem repeat consisting of six 100-bp copies and one 96-bp copy. A total of 315 bp of intergenic spacer sequence was spread over 23 regions. An alignment of the spacer sequences in ants was largely feasible among congeneric species, and there was substantial sequence divergence, indicating their potential use as molecular markers for congeneric species. The A/T contents at the first and second codon positions of protein-coding genes (PCGs) were similar for ant species, including C. atrox (73.9% vs. 72.3%, on average). With increased taxon sampling among hymenopteran superfamilies, differences in the divergence rates (i.e., the non-synonymous substitution rates) between the suborders Symphyta and Apocrita were detected, consistent with previous results. The C. atrox mt genome had a unique gene arrangement, trnI-trnM-trnQ, at the A+T-rich region and ND2 junction (underline indicates inverted gene). This may have originated from a tandem duplication of trnM-trnI, resulting in trnM-trnI-trnM-trnI-trnQ, and the subsequent loss of the first trnM and second trnI, resulting in trnI-trnM-trnQ.

The mitochondrial genome of the German wasp Vespula germanica (Fabricius, 1793) (Hymenoptera: Vespoidea: Vespidae)

The mitochondrial genome of the German wasp Vespula germanica (Fabricius, 1793) (Hymenoptera: Vespidae) (GenBank accession no. KR703583) was sequenced in the study. It represents the first mitochondrial genome from the genus Vespula. There are totally 163 42 bp in the currently sequenced portion of the genome, containing 13 protein-coding, two rRNA, and 18 tRNA genes and a partial A + T-rich region. Four tRNA genes of trnI, trnQ, trnM and trnY located at the downstream of the A + T-rich region were failed to sequence. At least two rearrangement events occurred in the sequenced region compared with the pupative ancestral arrangement of insects, corresponding to the translocation or remote inversion of tnnY from trnW-trnC-trnY cluster to the region of trnI-trnQ-trnM cluster and translocation of trnL1 from the downstream to the upstream of nad1 gene. All protein-coding genes start with ATN codons. Twelve and one protein-coding genes stop with termination codon TAA and T, respectively. Phylogenetic analysis using the Bayesian method based on all codon positions of the 13 protein-coding genes supports the monophyly of Vespidae and Formicidae. Within the Formicidae, the Myrmicinae and Formicinae form a sister group and then sister to the Dolichoderinae, while within the Vespidae, the Eumeninae sister to the lineage of Vespinae + Polistinae.

Phylogenetic analysis of the complete mitogenome sequence of the raspberry weevil, Aegorhinus superciliosus (Coleoptera: Curculionidae), supports monophyly of the tribe Aterpini

The superfamily Curculionoidea is one of the most diverse groups of insects in the world, including many species which are crop pests. Within this group, the native raspberry weevil, Aegorhinus superciliosus (Guérin, 1830), is an important pest in blueberry and raspberry fields in southern South America. Using a 454 sequencing approach, we sequenced and annotated the mitogenome of A. superciliosus, it, providing the first such information for any species in the tribe Aterpini, subfamily Cyclominae. The assembled mitogenome is a circular DNA molecule 15,121bp in length containing all 37 genes normally found in metazoans. Mitogenome organization and transcriptional orientation in A. superciliosus showed the same pattern that characterizes the suborder Polyphaga. Bayesian and Maximum Likelihood phylogenetic analyses supported the monophyly of the tribe Aterpini and the subfamily Cyclominae, recovering this clade in a sister group relationship with Entiminae and Hyperinae. The monophyly of these three subfamilies defines a critical transition to an ectophagous lifestyle in the larvae, from an ancestrally endophagous larval lifestyle in all other lineages. The sequenced mitogenome of A. superciliosus can provide basic data for future studies investigating population history, molecular systematics, stress ecophysiology and phylogeography.

The mitochondrial genome of Polistes jokahamae and a phylogenetic analysis of the Vespoidea (Insecta: Hymenoptera)

The mitochondrial genome sequence of Polistes jokahamae (Radoszkowski, 1887) (Hymenoptera: Vespidae) (GenBank accession no. KR052468) was sequenced. The current length with partial A + T-rich region of this mitochondrial genome is 16,616 bp. All the typical mitochondrial genes were sequenced except for three tRNAs (trnI, trnQ, and trnY) located between the A + T-rich region and nad2. At least three rearrangement events occurred in the sequenced region compared with the pupative ancestral arrangement of insects, corresponding to the shuffling of trnK and trnD, translocation or remote inversion of tnnY and translocation of trnL1. All protein-coding genes start with ATN codons. Eleven, one, and another one protein-coding genes stop with termination codon TAA, TA, and T, respectively. Phylogenetic analysis using the Bayesian method based on all codon positions of the 13 protein-coding genes supports the monophyly of Vespidae and Formicidae. Within the Formicidae, the Myrmicinae and Formicinae form a sister lineage and then sister to the Dolichoderinae, while within the Vespidae, the Eumeninae is sister to the lineage of Vespinae + Polistinae.

Phylogenetic utility, and variability in structure and content, of complete mitochondrial genomes among genetic lineages of the Hawaiian anchialine shrimp Halocaridina rubra Holthuis 1963 (Atyidae:Decapoda)

The Atyidae are caridean shrimp possessing hair-like setae on their claws and are important contributors to ecological services in tropical and temperate fresh and brackish water ecosystems. Complete mitochondrial genomes have only been reported from five of the 449 species in the family, thus limiting understanding of mitochondrial genome evolution and the phylogenetic utility of complete mitochondrial sequences in the Atyidae. Here, comparative analyses of complete mitochondrial genomes from eight genetic lineages of Halocaridina rubra, an atyid endemic to the anchialine ecosystem of the Hawaiian Archipelago, are presented. Although gene number, order, and orientation were syntenic among genomes, three regions were identified and further quantified where conservation was substantially lower: (1) high length and sequence variability in the tRNA-Lys and tRNA-Asp intergenic region; (2) a 317-bp insertion between the NAD6 and CytB genes confined to a single lineage and representing a partial duplication of CytB; and (3) the putative control region. Phylogenetic analyses utilizing complete mitochondrial sequences provided new insights into relationships among the H. rubra genetic lineages, with the topology of one clade correlating to the geologic sequence of the islands. However, deeper nodes in the phylogeny lacked bootstrap support. Overall, our results from H. rubra suggest intra-specific mitochondrial genomic diversity could be underestimated across the Metazoa since the vast majority of complete genomes are from just a single individual of a species.

Two nearly complete mitogenomes of wheat stem borers, Cephus pygmeus (L.) and Cephus sareptanus Dovnar-Zapolskij (Hymenoptera: Cephidae): an unusual elongation of rrnS gene

Two nearly complete mitochondrial genomes (mitogenomes) of wheat stem borers, Cephus pygmeus and Cephus sareptanus (Hymenoptera: Cephidae), were sequenced, characterised and compared with the previously known mitogenome of Cephus cinctus. The gene orders are mostly conserved, except for translocation of trnM and swapped position of trnI and trnQ. An A+T bias was found, but a deviation from strand asymmetry was also detected on the J strand. All protein coding genes (PCGs) are initiated by ATN codons, except for nad1, nad2 and atp8, and all are terminated with TAA, TA- or T- as a stop codon. The predicted secondary structures of rrnS and rrnL genes are mostly consistent with reported hymenopteran species. However, an unusual elongation in rrnS, not know elsewhere in the order, was discovered in Cephus species. Three autonomous sequences detected in domains I and II are mainly responsible for the length expansions.

Higher-level phylogeny of the Hymenoptera inferred from mitochondrial genomes

Higher-level hymenopteran relationships remain unresolved in both morphological and molecular analyses. In this study, we present the most comprehensive analyses of hymenopteran relationships based on 48 mitochondrial (mt) genomes. One complete and two nearly complete mt genomes representing three hymenopteran superfamilies were newly sequenced. We assessed the influence of inclusion/exclusion of 3rd codon positions, alignment approaches, partition schemes and phylogenetic approaches on topology and nodal support within the Hymenoptera. The results showed that the topologies were sensitive to the variation of dataset and analytical approach. However, some robust and highly supported relationships were recovered: the Ichneumonomorpha was monophyletic; the Trigonalyoidea+Megalyroidea and the Diaprioidea+Chalcidoidea were consistently recovered; the Cynipoidea was generally recovered as the sister group to the Diaprioidea+Chalcidoidea. In addition, the monophyletic Aculeata and Proctotrupomorpha were recovered in some analyses. Several gene rearrangements were detected in each of the three newly sequenced mt genomes. Specifically, the Ibalia leucospoides mt genome harbors a large inversion of a gene block from trnE to trnS2. Inverted, duplicated A+T rich regions were detected in the Ibalia leucospoides mt genome, which probably played an important role during the formation of the large gene block inversion via recombination.

Complete mitochondrial genomes of Ceratobaeus sp. and Idris sp. (Hymenoptera: Scelionidae): shared gene rearrangements as potential phylogenetic markers at the tribal level

We sequenced the complete mitochondrial genomes of two sceliond taxa (Ceratobaeus sp. and Idris sp.). An atypical tRNA-Arg which lacks a D-stem was identified in both taxa, and represents a potentially derived character of sceliond wasps. A number of tRNA genes have rearranged in the two mitochondrial genomes compared with the ancestral organization. Some of these derived genome organizations are shared, and thus have much potential as phylogenetic markers at the tribal level in the subfamily Scelioninae. We test the influence of third codon inclusion/exclusion, alignment methods and partition schemes on the reconstruction of phylogenetic relationships. The results show that inclusion of third codon positions does not appear to be problematic when investigating the phylogeny of closely related taxa. Muscle and PartitionFinder schemes significantly improve the likelihood scores.

The mitochondrial genome of the Vespa bicolor Fabricius (Hymenoptera: Vespidae: Vespinae)

Here, we report the first representative mitochondrial genome of the subfamily Vespinae (Hymenopterea: Vespidae) from the Vespa bicolor Fabricius (GenBank accession No. KJ735511). Nearly complete mitochondrial genome was sequenced with a length of 16,937 bp, including 13 protein-coding genes, 2 rRNA and 20 tRNA genes, as well as a portion of A+T-rich region. Two tRNA genes, i.e. trnI and trnY, were failed to sequence, which were presumed to be located within a region between A+T-rich region and trnM-trnQ-nad2. In the V. bicolor mitochondrial genome, at least three tRNA genes were rearranged. trnY was rearranged to the unsequenced region between A+T-rich region and trnM-trnQ-nad2. trnL1 was rearranged from a location between nad1 and rrnL to the upstream of nad1 gene. trnS1 and trnE were shuffled in the tRNA cluster of trnA-trnR-trnN-trnS1-trnE-trnF. Our study showed that the mitochondrial genomes between Vespinae and Polistinae shared more arrangement pattern than that between Vespinae and Eumeninae.

The mitochondrial genome of the leaf-cutter ant Atta laevigata: a mitogenome with a large number of intergenic spacers

In this paper we describe the nearly complete mitochondrial genome of the leaf-cutter ant Atta laevigata, assembled using transcriptomic libraries from Sanger and Illumina next generation sequencing (NGS), and PCR products. This mitogenome was found to be very large (18,729 bp), given the presence of 30 non-coding intergenic spacers (IGS) spanning 3,808 bp. A portion of the putative control region remained unsequenced. The gene content and organization correspond to that inferred for the ancestral pancrustacea, except for two tRNA gene rearrangements that have been described previously in other ants. The IGS were highly variable in length and dispersed through the mitogenome. This pattern was also found for the other hymenopterans in particular for the monophyletic Apocrita. These spacers with unknown function may be valuable for characterizing genome evolution and distinguishing closely related species and individuals. NGS provided better coverage than Sanger sequencing, especially for tRNA and ribosomal subunit genes, thus facilitating efforts to fill in sequence gaps. The results obtained showed that data from transcriptomic libraries contain valuable information for assembling mitogenomes. The present data also provide a source of molecular markers that will be very important for improving our understanding of genomic evolutionary processes and phylogenetic relationships among hymenopterans.