Is it an ant or a butterfly? Convergent evolution in the mitochondrial gene order of Hymenoptera and Lepidoptera

A Mysterious Asian Firefly Genus, Oculogryphus Jeng, Engel & Yang (Coleoptera, Lampyridae): The First Complete Mitochondrial Genome and Its Phylogenetic Implications

The firefly genus Oculogryphus Jeng, Engel & Yang, 2007 is a rare-species group endemic to Asia. Since its establishment, its position has been controversial but never rigorously tested. To address this perplexing issue, we are the first to present the complete mitochondrial sequence of Oculogryphus, using the material of O. chenghoiyanae Yiu & Jeng, 2018 determined through a comprehensive morphological identification. Our analyses demonstrate that its mitogenome exhibits similar characteristics to that of Stenocladius, including a rearranged gene order between trnC and trnW, and a long intergenic spacer (702 bp) between the two rearranged genes, within which six remnants (29 bp) of trnW were identified. Further, we incorporated this sequence into phylogenetic analyses of Lampyridae based on different molecular markers and datasets using ML and BI analyses. The results consistently place Oculogryphus within the same clade as Stenocladius in all topologies, and the gene rearrangement is a synapomorphy for this clade. It suggests that Oculogryphus should be classified together with Stenocladius in the subfamily Ototretinae at the moment. This study provides molecular evidence confirming the close relationship between Oculogryphus and Stenocladius and discovers a new phylogenetic marker helpful in clarifying the monophyly of Ototretinae, which also sheds a new light on firefly evolution.

Comparative mitogenomic analysis provides evolutionary insights into Formica (Hymenoptera: Formicidae)

Formica is a large genus in the family Formicidae with high diversity in its distribution, morphology, and physiology. To better understand evolutionary characteristics of Formica, the complete mitochondrial genomes (mitogenomes) of two Formica species were determined and a comparative mitogenomic analysis for this genus was performed. The two newly sequenced Formica mitogenomes each included 37 typical mitochondrial genes and a large non-coding region (putative control region), as observed in other Formica mitogenomes. Base composition, gene order, codon usage, and tRNA secondary structure were well conserved among Formica species, whereas diversity in sequence size and structural characteristics was observed in control regions. We also observed several conserved motifs in the intergenic spacer regions. These conserved genomic features may be related to mitochondrial function and their highly conserved physiological constraints, while the diversity of the control regions may be associated with adaptive evolution among heterogenous habitats. A negative AT-skew value on the majority chain was presented in each of Formica mitogenomes, indicating a reversal of strand asymmetry in base composition. Strong codon usage bias was observed in Formica mitogenomes, which was predominantly determined by nucleotide composition. All 13 mitochondrial protein-coding genes of Formica species exhibited molecular signatures of purifying selection, as indicated by the ratio of non-synonymous substitutions to synonymous substitutions being less than 1 for each protein-coding gene. Phylogenetic analyses based on mitogenomic data obtained fairly consistent phylogenetic relationships, except for two Formica species that had unstable phylogenetic positions, indicating mitogenomic data are useful for constructing phylogenies of ants. Beyond characterizing two additional Formica mitogenomes, this study also provided some key evolutionary insights into Formica.

Phylogenomics resolves the higher-level phylogeny of herbivorous eriophyoid mites (Acariformes: Eriophyoidea)

BACKGROUND

Eriophyoid mites (Eriophyoidea) are among the largest groups in the Acariformes; they are strictly phytophagous. The higher-level phylogeny of eriophyoid mites, however, remains unresolved due to the limited number of available morphological characters-some of them are homoplastic. Nevertheless, the eriophyoid mites sequenced to date showed highly variable mitochondrial (mt) gene orders, which could potentially be useful for resolving the higher-level phylogenetic relationships.

RESULTS

Here, we sequenced and compared the complete mt genomes of 153 eriophyoid mite species, which showed 54 patterns of rearranged mt gene orders relative to that of the hypothetical ancestor of arthropods. The shared derived mt gene clusters support the monophyly of eriophyoid mites (Eriophyoidea) as a whole and the monophylies of six clades within Eriophyoidea. These monophyletic groups and their relationships were largely supported in the phylogenetic trees inferred from mt genome sequences as well. Our molecular dating results showed that Eriophyoidea originated in the Triassic and diversified in the Cretaceous, coinciding with the diversification of angiosperms.

CONCLUSIONS

This study reveals multiple molecular synapomorphies (i.e. shared derived mt gene clusters) at different levels (i.e. family, subfamily or tribe level) from the complete mt genomes of 153 eriophyoid mite species. We demonstrated the use of derived mt gene clusters in unveiling the higher-level phylogeny of eriophyoid mites, and underlines the origin of these mites and their co-diversification with angiosperms.

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.

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.

Phylomitogenomics elucidates the evolution of symbiosis in Thoracotremata (Decapoda: Cryptochiridae, Pinnotheridae, Varunidae)

Background

Thoracotremata belong to the large group of "true" crabs (infraorder Brachyura), and they exhibit a wide range of physiological and morphological adaptations to living in terrestrial, freshwater and marine habitats. Moreover, the clade comprises various symbiotic taxa (Aphanodactylidae, Cryptochiridae, Pinnotheridae, some Varunidae) that are specialised in living with invertebrate hosts, but the evolutionary history of these symbiotic crabs is still partially unresolved.

Methods

Here we assembled and characterised the complete mitochondrial genomes (hereafter mitogenomes) of three gall crab species (Cryptochiridae): Kroppcarcinus siderastreicola, Opecarcinus hypostegus and Troglocarcinus corallicola. A phylogenetic tree of the Thoracotremata was reconstructed using 13 protein-coding genes and two ribosomal RNA genes retrieved from three new gall crab mitogenomes and a further 72 available thoracotreme mitogenomes. Furthermore, we applied a comparative analysis to characterise mitochondrial gene order arrangement, and performed a selection analysis to test for selective pressure of the protein-coding genes in symbiotic Cryptochiridae, Pinnotheridae, and Varunidae (Asthenognathus inaequipes and Tritodynamia horvathi).

Results

The results of the phylogenetic reconstruction confirm the monophyly of Cryptochiridae, which clustered separately from the Pinnotheridae. The latter clustered at the base of the tree with robust branch values. The symbiotic varunids A. inaequipes and T. horvathi clustered together in a clade with free-living Varunidae species, highlighting that symbiosis in the Thoracotremata evolved independently on multiple occasions. Different gene orders were detected in symbionts and free-living species when compared with the ancestral brachyuran gene order. Lastly, the selective pressure analysis detected two positively selected sites in the nad6 gene of Cryptochiridae, but the evidence for positive selection in Pinnotheridae and A. inaequipes and T. horvathi was weak. Adaptive evolution of mitochondrial protein-coding genes is perhaps related to the presumably higher energetic demands of a symbiotic lifestyle.

The complete mitochondrial genome of Eulaelaps huzhuensis (Mesostigmata: Haemogamasidae)

Some mites of the family Haemogamasidae can transmit a variety of zoonotic diseases and have important public health and safety implications. Currently, however, little attention has been paid to molecular data of Haemogamasidae species, limiting our understanding of their evolutionary and phylogenetic relationships. In this study, the complete mitochondrial genome of Eulaelaps huzhuensis was determined for the first time, and its genomic information was analyzed in detail. The mitochondrial genome of E. huzhuensis is 14,872 bp in length with 37 genes and two control regions. The base composition showed a distinct AT preference. Twelve protein-coding genes have a typical ATN as the start codon, and three protein-coding genes have incomplete stop codons. During the folding of tRNA genes, a total of 30 mismatches occurred, and three tRNA genes had an atypical cloverleaf secondary structure. The order of the E. huzhuensis mitochondrial genome arrangement is a new type of rearrangement in Mesostigmata. The phylogenetic analysis confirmed that the family Haemogamasidae is a monophyletic branch and does not belong to a subfamily of the Laelapidae. Our results lay the foundation for subsequent studies on the phylogeny and evolutionary history of the family Haemogamasidae.

Evolution of mitogenomic gene order in Orthoptera

Mitochondrial gene order has contributed to the elucidation of evolutionary relationships in several animal groups. It generally has found its application as a phylogenetic marker for deep nodes. Yet, in Orthoptera limited research has been performed on the gene order, although the group represents one of the oldest insect orders. We performed a comprehensive study on mitochondrial genome rearrangements (MTRs) within Orthoptera in the context of mitogenomic sequence-based phylogeny. We used 280 published mitogenome sequences from 256 species, including three outgroup species, to reconstruct a molecular phylogeny. Using a heuristic approach, we assigned MTR scenarios to the edges of the phylogenetic tree and reconstructed ancestral gene orders to identify possible synapomorphies in Orthoptera. We found all types of MTRs in our dataset: inversions, transpositions, inverse transpositions, and tandem-duplication/random loss events (TDRL). Most of the suggested MTRs were in single and unrelated species. Out of five MTRs which were unique in subgroups of Orthoptera, we suggest four of them to be synapomorphies; those were in the infraorder Acrididea, in the tribe Holochlorini, in the subfamily Pseudophyllinae, and in the two families Phalangopsidae and Gryllidae or their common ancestor (leading to the relationship ((Phalangopsidae + Gryllidae) + Trigonidiidae)). However, similar MTRs have been found in distant insect lineages. Our findings suggest convergent evolution of specific mitochondrial gene orders in several species, deviant from the evolution of the mitogenome DNA sequence. As most MTRs were detected at terminal nodes, a phylogenetic inference of deeper nodes based on MTRs is not supported. Hence, the marker does not seem to aid resolving the phylogeny of Orthoptera, but adds further evidence for the complex evolution of the whole group, especially at the genetic and genomic levels. The results indicate a high demand for more research on patterns and underlying mechanisms of MTR events in Orthoptera.

Comparative analysis of the mitochondrial genome of Dermacentor steini from different regions in China

Ticks are a group of blood-sucking ectoparasites that play an important role in human health and livestock production development as vectors of zoonotic diseases. The phylogenetic tree of single genes cannot accurately reflect the true kinship between species. Based on the complete mitochondrial genome analysis one can help to elucidate the phylogenetic relationships among species. In this study, the complete mitochondrial genome of Dermacentor steini (isolate Longyan) was sequenced and compared with the mitochondrial genes of 3 other Chinese isolates (Nanchang, Jinhua and Yingtan). In Dermacentor steini 4 isolates had identical or similar mitochondrial genome lengths and an overall variation of 0.76% between sequences. All nucleotide compositions showed a distinct AT preference. The most common initiation and stop codons were ATG and TAA, respectively. Fewer base mismatches were found in the tRNA gene of D. steini (isolate Longyan), and the vicinity of the control region and tRNA gene was a hot rearrangement region of the genus Dermacentor. Maximum likelihood trees and Bayesian trees indicate that D. steini is most closely related to Dermacentor auratus. The results enrich the mitochondrial genomic data of species in the genus Dermacentor and provide novel insights for further studies on the phylogeographic classification and molecular evolution of ticks.

A novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

Characterization and Phylogenetic Implication of Complete Mitochondrial Genome of the Medicinal Ant Formica sinae (Hymenoptera: Formicidae): Genomic Comparisons in Formicidae

Ants (Hymenoptera: Formicidae) are among the largest and most widespread families of terrestrial insects and are valuable to medical and ecological investigations. The mitochondrial genome has been widely used as a reliable genetic marker for species identification and phylogenetic analyses. To further understand the mitogenome-level characteristics of the congeneric Formicidae species, the complete mitogenome of Formica sinae (Hymenoptera: Formicidae) was sequenced, annotated, and compared with other 48 Formicidae species. The results showed that gene composition, content, and codon usage were conserved. The complete mitochondrial genome of F. sinae was 17,432 bp, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and one control region located between rrnS and trnM, which was 1,256 bp long, the longest of all sequenced species. Gene rearrangement was not detected in Formica species (Hymenoptera: Formicidae). All PCGs of F. sinae were initiated with ATN codons and terminated with the TAA codon. The overall nucleotide composition of F. sinae was AT-biased (83.51%), being 80.58% in PCGs, 86.68% in tRNAs, 87.10% in rRNAs, and 88.70% in the control region. Phylogenetic analyses indicated that each subfamily formed a strongly monophyletic group. Furthermore, F. sinae clustered with Formica fusca (Hymenoptera: Formicidae) and Formica selysi (Hymenoptera: Formicidae). This work enhances the genetic data of Formicidae and contributes to our understanding of their phylogenic relationship, evolution, and utilization.

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.

The complete mitochondrial genome of Aphaenogaster japonica (Forel, 1911) (Hymenoptera: Formicidae)

Aphaenogaster japonica (Forel, 1911) is an omnivorous ant that is widely distributed in eastern Asia. The mitochondrial genome of A. japonica reported here was 18,607 bp in length, consisting of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a control region. The base composition was AT biased (the GC ratio is 18.9%). With A. japonica added, we obtain weak evidence that the sister group of the Stenammini group, including Aphaenogaster, is the Myrmicini group. Therefore, the Stenammini and Myrmicini groups may be not a robust monophyletic group, unlike the previous results based on the complete mitochondrial genome.

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.

Phylogenetic relationship and characterization of the complete mitochondrial genome of Camponotus japonicus (Hymenoptera: Formicoidea: Formicidae)

Camponotus japonicus Mayr, 1866 is a widespread and open-field formicine ant species in east Asia. In this study, we sequenced and analyzed the complete mitochondrial genome (mitogenome) of C. japonicus. This mitogenome was 16,422 bp long and encoded 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs) and two ribosomal RNA unit genes (rRNAs). Gene order was conserved and identical to most other previously sequenced Formicidae. All PCGs of C. japonicus have the conventional start codon for invertebrate mitochondrial PCGs (five ATT, four ATA and four ATG). All PCGs terminated with the stop codon TAA. The whole mitogenome exhibited heavy AT nucleotide bias (77.8%). Phylogenetic analysis positioned C. japonicus in a well-supported clade with C. atrox and C. concavus, and C. japonicus was more closely related to C. atrox than to C. concavus. Within Formicinae, the topology (Colobopsis + (Polyrhachis + Camponotus) + (Formica + (Nylanderia + Lasius) + (Anoplolepis + Acropyga))) was recovered.

Highly diversified mitochondrial genomes provide new evidence for interordinal relationships in the Arachnida

Arachnida is an exceptionally diverse class in the Arthropoda, consisting of 20 orders and playing crucial roles in the terrestrial ecosystems. However, their interordinal relationships have been debated for over a century. Rearranged or highly rearranged mitochondrial genomes (mitogenomes) were consistently found in this class, but their various extent in different lineages and efficiency for resolving arachnid phylogenies are unclear. Here, we reconstructed phylogenetic trees using mitogenome sequences of 290 arachnid species to decipher interordinal relationships as well as diversification through time. Our results recovered monophyly of ten orders (i.e. Amblypygi, Araneae, Ixodida, Mesostigmata, Opiliones, Pseudoscorpiones, Ricinulei, Sarcoptiformes, Scorpiones and Solifugae), while rejecting monophyly of the Trombidiformes due to the unstable position of the Eriophyoidea. The monophyly of Acari (subclass) was rejected, possibly due to the long-branch attraction of the Pseudoscorpiones. The monophyly of Arachnida was further rejected because the Xiphosura nested within arachnid orders with unstable positions. Mitogenomes that are highly rearranged in mites but less rearranged or conserved in the remaining lineages point to their exceptional diversification in mite orders; however, shared derived mitochondrial (mt) gene clusters were found within superfamilies rather than interorders, confusing phylogenetic signals in arachnid interordinal relationships. Molecular dating results show that arachnid orders have ancient origins, ranging from the Ordovician to the Carboniferous, yet have significantly diversified since the Cretaceous in orders Araneae, Mesostigmata, Sarcoptiformes, and Trombidiformes. By summarizing previously resolved key positions of some orders, we propose a plausible arachnid tree of life. Our results underline a more precise framework for interordinal phylogeny in the Arachnida and provide new insights into their ancient evolution.

The complete mitochondrial genome of Pheidole nodus (Smith, 1874) (Hymenoptera: Formicidae)

Pheidole nodus (Smith, 1874) belongs to a famously hyperdiverse and ecologically dominant ant genus. The mitochondrial genome of P. nodus is 15,579 bp in length, and the overall base composition is 78.6% AT. It includes 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNAs, and a control region. Phylogenetic trees show that P. nodus is more closely related to Wasmannia than to Atta. These sequence data will play an important role in the investigation of the phylogenetic relationships and taxonomy of the group Attini.

The complete mitochondrial genome of Tetramorium tsushimae (Emery, 1925) (Hymenoptera: Formicidae)

Tetramorium tsushimae (Emery, 1925) is an omnivorous ant species native to East Asia and has been introduced to North America. The mitochondrial genome of T. tsushimae is 19,207 bp in length with an A + T content of 81.3% and includes 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNAs, and a control region. This sequence data would play an important role in the investigation of mitochondrial evolution of the subfamily Myrmicinae.

The complete mitochondrial genome of a parthenogenetic ant Monomorium triviale (Hymenoptera: Formicidae)

Monomorium is one of the most species-rich yet taxonomically problematic ant genus. An East Asian species, M. triviale Wheeler, W.M., 1906, is reproduced by obligate thelytokous parthenogenesis and performs strict reproductive division of labor. We sequenced the M. triviale mitogenome using next-generation sequencing methods. The circular mitogenome of M. triviale was 16,290 bp in length, consisting of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and a single non-coding region of 568 bp. The base composition was AT-biased (82%). Gene order rearrangements were detected and likely to be unique to the genus Monomorium. We announce the M. triviale mitogenome as additional genomic resources for elucidating phylogenetic and taxonomic problems of Monomorium and comparative genomics of parthenogenetic ant species.

How are the mitochondrial genomes reorganized in Hexapoda? Differential evolution and the first report of convergences within Hexapoda

The evolution of the Hexapoda mitochondrial genome has been the focus of several genetic and evolutionary studies over the last decades. However, they have concentrated on certain taxonomic orders of economic or health importance. The recent increase of mitochondrial genomes sequencing of diverse taxonomic orders generates an important opportunity to clarify the evolution of this group of organisms. However, there is no comparative study that investigates the evolution of the Hexapoda mitochondrial genome. In order to verify the level of rearrangement and the mitochondrial genome evolution, we performed a comparative genomic analysis of the Hexapoda mitochondrial genome available in the NCBI database. Using a combination of bioinformatics methods to carefully examine the mitochondrial gene rearrangements in 1198 Hexapoda species belonging to 32 taxonomic orders, we determined that there is a great variation in the rate of rearrangement by gene and by taxonomic order. A higher rate of genetic reassortment is observed in Phthiraptera, Thysanoptera, Protura, and Hymenoptera; compared to other taxonomic orders. Twenty-four events of convergence in the genetic order between different taxonomic orders were determined, most of them not previously reported; which proves the great evolutionary dynamics within Hexapoda.

Complete Genome Sequence of a Blochmannia Endosymbiont of Colobopsis nipponica

Blochmannia endosymbionts (Gammaproteobacteria) live in bacteriocytes, which are specialized cells found in the genus Camponotus and its neighbor genera. In this announcement, we describe the complete genome sequence of the Blochmannia endosymbiont of Colobopsis nipponica, which originated from a colony collected in the Republic of Korea.

Eusociality Shapes Convergent Patterns of Molecular Evolution across Mitochondrial Genomes of Snapping Shrimps

Eusociality is a highly conspicuous and ecologically impactful behavioral syndrome that has evolved independently across multiple animal lineages. So far, comparative genomic analyses of advanced sociality have been mostly limited to insects. Here, we study the only clade of animals known to exhibit eusociality in the marine realm-lineages of socially diverse snapping shrimps in the genus Synalpheus. To investigate the molecular impact of sociality, we assembled the mitochondrial genomes of eight Synalpheus species that represent three independent origins of eusociality and analyzed patterns of molecular evolution in protein-coding genes. Synonymous substitution rates are lower and potential signals of relaxed purifying selection are higher in eusocial relative to noneusocial taxa. Our results suggest that mitochondrial genome evolution was shaped by eusociality-linked traits-extended generation times and reduced effective population sizes that are hallmarks of advanced animal societies. This is the first direct evidence of eusociality impacting genome evolution in marine taxa. Our results also strongly support the idea that eusociality can shape genome evolution through profound changes in life history and demography.

Comparative mitogenomics and phylogenetics of the stinging wasps (Hymenoptera: Aculeata)

The stinging wasps (Hymenoptera: Aculeata) include diverse groups such as vespid wasps, ants and bees. Phylogenetic relationships among major lineages of stinging wasps have been inferred from molecular and morphological data. However, the genomic features of the mitochondrial genomes and their phylogenetic utility remain to be explored. In this study, we determined 23 mitochondrial genomes from the Aculeata. Four Mutillidae species showed relatively low A + T content compared to other species of the Aculeata (69.7%-77.4%). Eleven out of 44 species, mainly from the Chrysididae and the Pompilidae, showed reversals of GC skews. Gene rearrangements occurred across the species. Patterns of tRNA rearrangement were conserved in some groups, including the Chrysididae, Bethylidae, Pompilidae, Scolioidea and Vespoidea. Rearrangement of protein-coding genes were found in 12 out of 44 species of the Aculeata, including all four species from the Chrysididae, both species from the Bethylidae, one species from the Dryinidae, all three Scolioidea species and two Apoidea species. Phylogenetic inference showed a long branch in species with unusual genomic features, such as in the Mutillidae and Bethylidae. By excluding these species, we found paraphyly of the Chrysidoidea and a sister group relationship between the Formicoidea and Vespoidea. These results improve our understanding of the evolution of mitochondrial genomes in the Aculeata and, in general, the evolution across this subclade.

Complete mitochondrial genome of the acrobat ant Crematogaster teranishii Santschi, 1930 (Formicidae; Hymenoptera)

The genus Crematogaster is a diverse group of ants found around the world. We have completed the mitochondrial genome of Crematogaster teranishii, which is the first mitochondrial genome of the genus. The mitochondrial genome is 17,442 bp long and 20.3% in GC ratio, which is similar to those of other ants. It contains 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and a control region with same gene order to other myrmicine species. The intergenic region between nad3 and trnA was unusually long compared to other ant species. Phylogenetic analysis showed that C. teranishii was closely related to other members of tribe Crematogastrini.

Complete mitochondrial genome of the jet ant Lasius spathepus Wheeler, W.M., 1910 (Formicidae; Hymenoptera)

The genus Lasius is a conspicuous and popular genus of ants found in the Holarctic regions. We have completed the mitochondrial genome of Lasius spathepus as the first mitochondrial genome of Lasius. The mitochondrial genome is 18,951 bp long, which is the sixth longest ant mitochondrial genome known to science. It contains 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a control region in a gene order shared with other species of subfamily Formicinae. The control region is 2,147 bp long, longest of all ants. Phylogenetic analysis shows L. spathepus groups with Nylanderia flavipes of the same tribe Lasiini.

Not Frozen in the Ice: Large and Dynamic Rearrangements in the Mitochondrial Genomes of the Antarctic Fish

The vertebrate mitochondrial genomes generally present a typical gene order. Exceptions are uncommon and important to study the genetic mechanisms of gene order rearrangements and their consequences on phylogenetic output and mitochondrial function. Antarctic notothenioid fish carry some peculiar rearrangements of the mitochondrial gene order. In this first systematic study of 28 species, we analyzed known and undescribed mitochondrial genome rearrangements for a total of eight different gene orders within the notothenioid fish. Our reconstructions suggest that transpositions, duplications, and inversion of multiple genes are the most likely mechanisms of rearrangement in notothenioid mitochondrial genomes. In Trematominae, we documented an extremely rare inversion of a large genomic segment of 5,300 bp that partially affected the gene compositional bias but not the phylogenetic output. The genomic region delimited by nad5 and trnF, close to the area of the Control Region, was identified as the hot spot of variation in Antarctic fish mitochondrial genomes. Analyzing the sequence of several intergenic spacers and mapping the arrangements on a newly generated phylogeny showed that the entire history of the Antarctic notothenioids is characterized by multiple, relatively rapid, events of disruption of the gene order. We hypothesized that a pre-existing genomic flexibility of the ancestor of the Antarctic notothenioids may have generated a precondition for gene order rearrangement, and the pressure of purifying selection could have worked for a rapid restoration of the mitochondrial functionality and compactness after each event of rearrangement.

Quantification and evolution of mitochondrial genome rearrangement in Amphibians

BACKGROUND

Rearrangement is an important topic in the research of amphibian mitochondrial genomes ("mitogenomes" hereafter), whose causes and mechanisms remain enigmatic. Globally examining mitogenome rearrangements and uncovering their characteristics can contribute to a better understanding of mitogenome evolution.

RESULTS

Here we systematically investigated mitogenome arrangements of 232 amphibians including four newly sequenced Dicroglossidae mitogenomes. The results showed that our new sequenced mitogenomes all possessed a trnM tandem duplication, which was not exclusive to Dicroglossidae. By merging the same arrangements, the mitogenomes of ~ 80% species belonged to the four major patterns, the major two of which were typical vertebrate arrangement and typical neobatrachian arrangement. Using qMGR for calculating rearrangement frequency (RF) (%), we found that the control region (CR) (RF = 45.04) and trnL2 (RF = 38.79) were the two most frequently rearranged components. Forty-seven point eight percentage of amphibians possessed rearranged mitogenomes including all neobatrachians and their distribution was significantly clustered in the phylogenetic trees (p < 0.001). In addition, we argued that the typical neobatrachian arrangement may have appeared in the Late Jurassic according to possible occurrence time estimation.

CONCLUSION

It was the first global census of amphibian mitogenome arrangements from the perspective of quantity statistics, which helped us to systematically understand the type, distribution, frequency and phylogenetic characteristics of these rearrangements.

Comparative Mitochondrial Analysis of Cnaphalocrocis exigua (Lepidoptera: Crambidae) and Its Close Relative C. medinalis

Rice leaffolders are important pests on rice in Asia, Oceania, and Africa, causing serious loss to rice production. There are two main rice leaffolders in China, namely Cnaphalocrocis medinalis (Guenée) and C. exigua (Butler) with the former having the ability of long-distance migration. To reveal the differences in the mitochondrial genomes (mitogenome) between them, we compared the completed mitogenome of C. exigua with three C. medinalis individuals. Although phylogenetic analysis based on the mitogenomic data strongly supported the close relationship between these two species, many differences were still being revealed. The results showed that the mitogenome of C. exigua was shorter in length (15,262 bp) and slight lower in AT content than that of C. medinalis. Except for the different start codons of nad3 and nad6 gene, we also found the cox1 gene had a typical start codon 'ATG' which suggested that the starting position of this gene must be reconsidered in the entire superfamily Pyraloidea. All tRNAs have a typical clover-leaf structure, except for the dihydrouridine (DHU) stem losing of trnS1, which has the atypical anticondon 'TCT' instead of 'GCT' in C. medinalis and most Pyraloidea species. Two intergenic regions (between trnY and cox1, nad3 and trnA) featured by AT repeats were only found in C. medinalis and even rarely appeared in reported Pyraloidea species. Furthermore, regardless of interspecific comparison or intraspecific comparison of these two species, protein coding genes, especially the atp8 genes, had quite different evolutionary rates.

The complete mitochondrial genome of the endangered Assam Roofed Turtle, Pangshura sylhetensis (Testudines: Geoemydidae): Genomic features and phylogeny

The Assam Roofed Turtle, Pangshura sylhetensis is an endangered and least studied species endemic to India and Bangladesh. The present study decodes the first complete mitochondrial genome of P. sylhetensis (16,568 bp) by using next-generation sequencing. The assembly encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). Most of the genes were encoded on the majority strand, except NADH dehydrogenase subunit 6 (nad6) and eight tRNAs. All PCGs start with an ATG initiation codon, except for Cytochrome oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 5 (nad5), which both start with GTG codon. The study also found the typical cloverleaf secondary structures in most of the predicted tRNA structures, except for serine (trnS1) which lacks of conventional DHU arm and loop. Both Bayesian and maximum-likelihood phylogenetic inference using 13 concatenated PCGs demonstrated strong support for the monophyly of all 52 Testudines species within their respective families and revealed Batagur trivittata as the nearest neighbor of P. sylhetensis. The mitogenomic phylogeny with other amniotes is congruent with previous research, supporting the sister relationship of Testudines and Archosaurians (birds and crocodilians). Additionally, the mitochondrial Gene Order (GO) analysis indicated plesiomorphy with the typical vertebrate GO in most of the Testudines species.

Mitochondrial and Innate Immunity Transcriptomes from Spodoptera frugiperda Larvae Infected with the Spodoptera frugiperda Ascovirus

Ascoviruses are large, enveloped DNA viruses that induce remarkable changes in cellular architecture during which the cell is partitioned into numerous vesicles for viral replication. Previous studies have shown that these vesicles arise from a process resembling apoptosis yet which differs after nuclear lysis in that mitochondria are not degraded but are modified by the virus, changing in size, shape, and motility. Moreover, infection does not provoke an obvious innate immune response. Thus, we used in vivo RNA sequencing to determine whether infection by the Spodoptera frugiperda ascovirus 1a (SfAV-1a) modified expression of host mitochondrial, cytoskeletal, and innate immunity genes. We show that transcripts from many mitochondrial genes were similar to those from uninfected controls, whereas others increased slightly during vesicle formation, including those for ATP6, ATP8 synthase, and NADH dehydrogenase subunits, supporting electron microscopy (EM) data that these organelles were conserved for virus replication. Transcripts from 58 of 106 cytoskeletal genes studied increased or decreased more than 2-fold postinfection. More than half coded for mitochondrial motor proteins. Similar increases occurred for innate immunity transcripts and their negative regulators, including those for Toll, melanization, and phagocytosis pathways. However, those for many antimicrobial peptides, such as moricin, increased more than 20-fold. In addition, transcripts for gloverin-3, spod_x_tox, Hdd23, and lebocin, also antimicrobial, increased more than 20-fold. Interestingly, a phenoloxidase inhibitor transcript increased 12-fold, apparently to interfere with melanization. SfAV-1a destroys most fat body cells by 7 days postinfection, so innate immunity gene transcripts apparently occur in remaining cells in this tissue and possibly other major tissues, namely, epidermis and tracheal matrix.IMPORTANCE Ascoviruses are large DNA viruses that infect insects, inducing a cellular pathology that resembles apoptosis but which differs by causing enormous cellular hypertrophy followed by cleavage of the cell into numerous viral vesicles for replication. Previous EM studies suggest that mitochondria are important for vesicle formation. Transcriptome analyses of Spodoptera frugiperda larvae infected with SfAV-1a showed that mitochondrial transcripts were similar to those from uninfected controls or increased slightly during vesicle formation, especially for ATP6, ATP8 synthase, and NADH dehydrogenase subunits. This pattern resembles that for chronic disease-inducing viruses, which conserve mitochondria, differing markedly from viruses causing short-term viral diseases, which degrade mitochondrial DNA. Though mitochondrial transcript increases were low, our results demonstrate that SfAV-1a alters host mitochondrial expression more than any other virus. Regarding innate immunity, although SfAV-1a destroys most fat body cells, certain immunity genes were highly upregulated (greater than 20-fold), suggesting that these transcripts may originate from other tissues.

Mitochondrial genomes of genus Atta (Formicidae: Myrmicinae) reveal high gene organization and giant intergenic spacers

The ants of the genus Atta are considered important pests to agriculture in the Americas, although Atta species are also important contributors to ecosystem functions in the various habitats in which they occur. The aim of this study was to assemble four complete mitochondrial genomes of the genus Atta, construct the phylogenomic tree, and analyze the gene content, order, and organization. The mitogenomes of A. colombica, A. opaciceps, A. texana, and A. sexdens rubropilosa comprise 18,392, 19,257, 19,709, and 19,748 bp, respectively. The four Atta mitogenomes showed the charactistics typical of those of insects, with 13 protein-coding genes, 22 tRNAs, and 2 rRNAs, with genes displayed in the conventional order. Analysis for intergenic spacer regions showed that Atta intergenic spacers are larger than those of the outgroups. Phylogenomic analyses using partial cytochrome oxidase I gene sequences showed similar topologies to previous phylogenetic analyses, with high clade support values. We conclude that Atta mitogenomes are characterized by high conservation in gene order and have giant intergenic spacers in the genus Atta.

The complete mitochondrial genome of Aphaenogaster famelica (Smith, 1874) (Hymenoptera: Formicidae)

Aphaenogaster famelica (Smith, 1874) is a slender ant found in Korea. We have completed the mitochondrial genome of A. famelica of which length is 19,464 bp, the third longest ant mitogenome so far. It was AT-biased (GC ratio is 18.0%) and includes 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a control region. The extra bases of A. famelica are placed in intergenic regions each between trnV-trnM and trnI-trnQ. Gene order of A. famelica is identical to other Myrmicinae ants. Phylogenetic trees show that tribe Stenammini and tribe Myrmicini are closely related.

The complete mitochondrial genome of Nylanderia flavipes (Smith, 1874) (Hymenoptera: Formicidae)

Nylanderia flavipes (Smith, 1874) is a Formicine ant found in East Asia. We have completed mitochondrial genome of N. flavipes of which length is 16,687 bp including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNAs, and a control region. The base composition was AT-biased (GC ratio is 18.9%). Gene order of N. flavipes presents a unique inversion of trnP in comparison to Camponotus species. Our phylogenetic trees together with phylogenomic analysis suggest that the gene order rearrangements occurred independently in Camponotus and N. flavipes.

MitoZ: a toolkit for animal mitochondrial genome assembly, annotation and visualization

Mitochondrial genome (mitogenome) plays important roles in evolutionary and ecological studies. It becomes routine to utilize multiple genes on mitogenome or the entire mitogenomes to investigate phylogeny and biodiversity of focal groups with the onset of High Throughput Sequencing (HTS) technologies. We developed a mitogenome toolkit MitoZ, consisting of independent modules of de novo assembly, findMitoScaf (find Mitochondrial Scaffolds), annotation and visualization, that can generate mitogenome assembly together with annotation and visualization results from HTS raw reads. We evaluated its performance using a total of 50 samples of which mitogenomes are publicly available. The results showed that MitoZ can recover more full-length mitogenomes with higher accuracy compared to the other available mitogenome assemblers. Overall, MitoZ provides a one-click solution to construct the annotated mitogenome from HTS raw data and will facilitate large scale ecological and evolutionary studies. MitoZ is free open source software distributed under GPLv3 license and available at https://github.com/linzhi2013/MitoZ.

Homoplasy or plesiomorphy? Reconstruction of the evolutionary history of mitochondrial gene order rearrangements in the subphylum Neodermata

Recent mitogenomic studies have exposed a gene order (GO) shared by two classes, four orders and 31 species ('common GO') within the flatworm subphylum Neodermata. There are two possible hypotheses for this phenomenon: convergent evolution (homoplasy) or shared ancestry (plesiomorphy). To test those, we conducted a meta-analysis on all available mitogenomes to infer the evolutionary history of GO in Neodermata. To improve the resolution, we added a newly sequenced mitogenome that exhibited the common GO, Euryhaliotrema johni (Ancyrocephalinae), to the dataset. Phylogenetic analyses conducted on two datasets (nucleotides of all 36 genes and amino acid sequences of 12 protein coding genes) and four algorithms (MrBayes, RAxML, IQ-TREE and PhyloBayes) produced topology instability towards the tips, so ancestral GO reconstructions were conducted using TreeREx and MLGO programs using all eight obtained topologies, plus three unique topologies from previous studies. The results consistently supported the second hypothesis, resolving the common GO as a plesiomorphic ancestral GO for Neodermata, Cestoda, Monopisthocotylea, Cestoda + Trematoda and Cestoda + Trematoda + Monopisthocotylea. This allowed us to trace the evolutionary GO scenarios from each common ancestor to its descendants amongst the Monogenea and Cestoda classes, and propose that the common GO was most likely retained throughout all of the common ancestors, leading to the extant species possessing the common GO. Neodermatan phylogeny inferred from GOs was largely incongruent with all 11 topologies described above, but it did support the mitogenomic dataset in resolving Polyopisthocotylea as the earliest neodermatan branch. Although highly derived GOs might be of some use in resolving isolated taxonomic and phylogenetic uncertainties, we conclude that, due to the discontinuous nature of their evolution, they tend to produce artefactual phylogenetic relationships, which makes them unsuitable for phylogenetic reconstruction in Neodermata. Wider and denser sampling of neodermatan mitogenomic sequences will be needed to infer the evolutionary pathways leading to the observed diversity of GOs with confidence.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Mitochondrial genome reorganization provides insights into the relationship between oribatid mites and astigmatid mites (Acari: Sarcoptiformes: Oribatida)

Oribatida s.l. represents one of the most species-rich mite lineages, including two recognized groups: oribatid mites (Oribatida s.s., non-astigmatan oribatids) and astigmatid mites (Astigmata). However, the relationship between these two groups has been debated. Here, we sequenced the complete mitochondrial (mt) genome of one oribatid mite and one astigmatid mite, retrieved complete mt genomes of three oribatid mites, and compared them with two other oribatid mites and 12 astigmatid mites sequenced previously. We find that gene orders in the mt genomes of both oribatid mites and astigmatid mites are rearranged relative to the hypothetical ancestral arrangement of the arthropods. Based on the shared derived gene clusters in each mt genome group, rearranged mt genomes are roughly divided into two groups corresponding to each mite group (oribatid mites or astigmatid mites). Phylogenetic results show that Astigmata nested in Oribatida. The monophyly of Astigmata is recovered, while paraphyly of Oribatida s.s. is observed. Our results show that rearranged gene orders in the mt genomes characterize various lineages of oribatid mites and astigmatid mites, and have potential phylogenetic information for resolving the high-level (cohort or supercohort) phylogeny of Oribatida.

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.

Mitochondrial Genome Fragmentation Unites the Parasitic Lice of Eutherian Mammals

Organelle genome fragmentation has been found in a wide range of eukaryotic lineages; however, its use in phylogenetic reconstruction has not been demonstrated. We explored the use of mitochondrial (mt) genome fragmentation in resolving the controversial suborder-level phylogeny of parasitic lice (order Phthiraptera). There are approximately 5000 species of parasitic lice in four suborders (Amblycera, Ischnocera, Rhynchophthirina, and Anoplura), which infest mammals and birds. The phylogenetic relationships among these suborders are unresolved despite decades of studies. We sequenced the mt genomes of eight species of parasitic lice and compared them with 17 other species of parasitic lice sequenced previously. We found that the typical single-chromosome mt genome is retained in the lice of birds but fragmented into many minichromosomes in the lice of eutherian mammals. The shared derived feature of mt genome fragmentation unites the eutherian mammal lice of Ischnocera (family Trichodectidae) with Anoplura and Rhynchophthirina to the exclusion of the bird lice of Ischnocera (family Philopteridae). The novel clade, namely Mitodivisia, is also supported by phylogenetic analysis of mt genome and cox1 gene sequences. Our results demonstrate, for the first time, that organelle genome fragmentation is informative for resolving controversial high-level phylogenies.

Accessible molecular phylogenomics at no cost: obtaining 14 new mitogenomes for the ant subfamily Pseudomyrmecinae from public data

The advent of Next Generation Sequencing has reduced sequencing costs and increased genomic projects from a huge amount of organismal taxa, generating an unprecedented amount of genomic datasets publicly available. Often, only a tiny fraction of outstanding relevance of the genomic data produced by researchers is used in their works. This fact allows the data generated to be recycled in further projects worldwide. The assembly of complete mitogenomes is frequently overlooked though it is useful to understand evolutionary relationships among taxa, especially those presenting poor mtDNA sampling at the level of genera and families. This is exactly the case for ants (Hymenoptera:Formicidae) and more specifically for the subfamily Pseudomyrmecinae, a group of arboreal ants with several cases of convergent coevolution without any complete mitochondrial sequence available. In this work, we assembled, annotated and performed comparative genomics analyses of 14 new complete mitochondria from Pseudomyrmecinae species relying solely on public datasets available from the Sequence Read Archive (SRA). We used all complete mitogenomes available for ants to study the gene order conservation and also to generate two phylogenetic trees using both (i) concatenated set of 13 mitochondrial genes and (ii) the whole mitochondrial sequences. Even though the tree topologies diverged subtly from each other (and from previous studies), our results confirm several known relationships and generate new evidences for sister clade classification inside Pseudomyrmecinae clade. We also performed a synteny analysis for Formicidae and identified possible sites in which nucleotidic insertions happened in mitogenomes of pseudomyrmecine ants. Using a data mining/bioinformatics approach, the current work increased the number of complete mitochondrial genomes available for ants from 15 to 29, demonstrating the unique potential of public databases for mitogenomics studies. The wide applications of mitogenomes in research and presence of mitochondrial data in different public dataset types makes the "no budget mitogenomics" approach ideal for comprehensive molecular studies, especially for subsampled taxa.

Information from the mitochondrial genomes of two egg parasitoids, Gonatocerus sp. and Telenomus sp., reveals a controversial phylogenetic relationship between Mymaridae and Scelionidae

The taxonomic status and phylogenetic affinities of Mymaridae and Scelionidae are controversial, based on similarities between these families in the characteristics of adults, larvae, and eggs. In this study, we sequenced the mitochondrial (mt) genomes of representatives from these two families and found that the derived secondary structure of tRNA-Arg was the same in each family due to the absence of the D-stem. The segment of "cox1 trnL₂cox2 trnK trnD atp8 atp6 cox3" in Gonatocerus sp. (Mymaridae) is conserved and distinct from those of four other species of Chalcidoidea but similar to that in Proctotrupoidea and Platygastroidea. However, phylogenetic analysis indicated that Gonatocerus sp. was sister group to other species of Chalcidoidea. Comparisons based on complete gene orders may be more useful in a phylogenetic and systematic context, as different branches may exhibit partially homoplastic gene orders.

The first two mitochondrial genomes of the family Aphelinidae with novel gene orders and phylogenetic implications

Chalcidoidea is one of the most diverse group in Hymenoptera by possessing striking mitochondrial gene arrangement. By using next generation sequencing method, the first two nearly complete mitochondrial genomes in the family Aphelinidae (Insecta, Hymenopetra, Chalcidoidea) were obtained in this study. Almost all previously sequenced mitochondrial genome of Chalcidoidea species have a large inversion including six genes (atp6-atp8-trnD-trnK-cox2-trnL2-cox1) as compared with ancestral mitochondrial genome, but these two Encarsia mitochondrial genomes had a large inversion including nine genes (nad3-trnG-atp6-atp8-trnD-trnK-cox2-trnL2-cox1), which was only congruent with the species in the genus Nasonia. Moreover, we found that one shuffling changes (trnD and trnK) happened in the species E. obtusiclava but not in another species E. formosa within the same genus, of which such shuffling within the same genus at this region was only detected in Polisters within Insecta. Phylogenetic analysis displayed that different data matrix (13PCG+ 2 rRNA or 13 PCG) and inference methods (BI or ML) indicate the identical topology with high nodal supports that Aphelinidae formed a sister group with (Trichogrammatidae + Aganoidae) and the monophyly of Pteramalidae. Our results also indicated the validity of assembling and feasibility of next-generation technology to obtain the mitochondrial genomes of parasitic Hymenoptera.

The complete mitochondrial genome of yellow crazy ant, Anoplolepis gracilipes (Hymenoptera: Formicidae)

The yellow crazy ant Anoplolepis gracilipes is an invasive species that threatens biodiversity in introduced ecosystems. We sequenced the A. gracilipes mitogenome using next-generation sequencing methods. The circular mitogenome of A. gracilipes was 16,943 bp included 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and a single large non-coding region of 893 bp. The base composition was AT-biased (72%). Three genomic rearrangements compared to ancestral insects were found. Phylogenetic analysis based on the concatenated nucleotide sequences of the 13 protein-coding genes supports A. gracilipes belonging to the Formicinae subfamily. We announce the A. gracilipes mitogenome as a DNA reference for further population genetic, phylogenetic, and evolutionary analyses.

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.