Phylogeny of Arthropoda inferred from mitochondrial sequences: strategies for limiting the misleading effects of multiple changes in pattern and rates of substitution

Drivers of interlineage variability in mitogenomic evolutionary rates in Platyhelminthes

Studies of forces driving interlineage variability in the evolutionary rates (both sequence and architecture) of mitochondrial genomes often produce contradictory results. Flatworms (Platyhelminthes) exhibit the fastest-evolving mitogenomic sequences among all bilaterian phyla. To test the effects of multiple factors previously associated with different aspects of mitogenomic evolution, we used mitogenomes of 223 flatworm species, phylogenetic multilevel regression models, and causal inference. Thermic host environment (endothermic vs. ectothermic) had nonsignificant impacts on both sequence evolution and mitogenomic size. Mitogenomic gene order rearrangements (GORR) were mostly positively correlated with mitogenomic size (R2 ≈ 20-30%). Longevity was not (negatively) correlated with sequence evolution in flatworms. The predominantly free-living "turbellaria" exhibited much shorter branches and faster-evolving mitogenomic architecture than parasitic Neodermata. As a result, "parasitism" had a strong explanatory power on the branch length variability (>90%), and there was a negative correlation between GORR and branch length. However, the stem branch of Neodermata comprised 63.6% of the total average branch length. This evolutionary period was also marked by a high rate of gene order rearrangements in the ancestral Neodermata. We discuss how this period of rapid evolution deep in the evolutionary history may have decoupled sequence evolution rates from longevity and GORR, and overestimated the explanatory power of "parasitism". This study shows that impacts of variables often vary across lineages, and stresses the importance accounting for the episodic nature of evolutionary patterns in studies of mitogenomic evolution.

Arthropod Phylotranscriptomics With a Special Focus on the Basal Phylogeny of the Myriapoda

Arthropoda represents the most diverse animal phylum, but clarifying the phylogenetic relationships among arthropod taxa remains challenging given the numerous arthropod lineages that diverged over a short period of time. In order to resolve the most controversial aspects of deep arthropod phylogeny, focusing on the Myriapoda, we conducted phylogenetic analyses based on ten super-matrices comprised of 751 to 1,233 orthologous genes across 64 representative arthropod species, including 28 transcriptomes that were newly generated in this study. Our findings provide unambiguous support for the monophyly of the higher arthropod taxa, Chelicerata, Mandibulata, Myriapoda, Pancrustacea, and Hexapoda, while the Crustacea are paraphyletic, with the class Remipedia supported as the lineage most closely related to hexapods. Within the Hexapoda, our results largely affirm previously proposed phylogenetic relationships among deep hexapod lineages, except that the Paraneoptera (Hemiptera, Thysanoptera, and Psocodea) was recovered as a monophyletic lineage in some analyses. The results corroborated the recently proposed phylogenetic framework of the four myriapod classes, wherein Symphyla and Pauropoda, as well as Chilopoda and Diplopoda, are each proposed to be sister taxa. The findings provide important insights into understanding the phylogeny and evolution of arthropods.

The mitochondrial genome of Bottapotamon fukienense (Brachiura: Potamidae) is fragmented into two chromosomes

BACKGROUND

China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.

RESULTS

Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.

CONCLUSIONS

This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

Dinoflagellate chloroplasts as a model for extreme genome reduction and fragmentation in organelles - The COCOA principle for gene retention

The genomes of peridinin-containing dinoflagellate chloroplasts have a very unusual organisation. These genomes are highly fragmented and greatly reduced, with most of the usual complement of chloroplast genes relocated to the nucleus. Dinoflagellate chloroplasts highlight evolutionary changes that are found to varying extents in a number of other organelle genomes. These include the chloroplast genome of the green alga Boodlea and other Cladophorales, and the mitochondrial genomes of blood-sucking and chewing lice, the parasitic plant Rhopalocnemis phalloides, the red alga Rhodosorus marinus and other members of the Stylonematophyceae, diplonemid flagellates, and some Cnidaria. Consideration of the coding content of the remnant chloroplast genomes indicates that organelles may preferentially retain genes for proteins important in initiating assembly of complexes, and the same is largely true for mitochondria. We propose a new principle, of CO-location for COntrol of Assembly (COCOA), indicating the importance of retaining these genes in the organelle. This adds to, but does not invalidate, the existing hypotheses of the multisubunit completion principle, CO-location for Redox Regulation (CORR) and Control by Epistasy of Synthesis (CES).

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Mitochondrial genome comparison and phylogenetic position of Fannia pusio among the Calyptratae flies

Fannia pusio, the chicken dung fly species, remains unexplored despite its forensic, sanitary, and veterinary importance in the Nearctic and Neotropical regions. In this study, we obtained the complete mitochondrial genome of Fannia pusio for the first time using next-generation sequencing. We compared it with previously published mitogenomes of the genus from the Palearctic region, and its phylogenetic position was studied based on the concatenated protein-coding genes (PCGs) dataset of Calyptratae flies. The circular mitochondrial genome of F. pusio is 16,176 bp in length, with a high A + T content (78.3%), whose gene synteny, codon usage analysis, and amino acid frequency are similar to previously reported Fannia mitogenomes. All PCGs underwent purifying selection except the nad2 gene. Interspecific K2P distances of PCGs of Fannia yielded an average of 12.4% (8.1%-21.1%). The Fannia genus is monophyletic and closely related to Muscidae based on molecular data. Further taxonomic sampling is required to deep into the phylogenetic relationships of the originally proposed species-groups and subgroups within the genus. These results provide a valuable dataset for studying the mitochondrial genome evolution and a resource for the taxonomy and systematics of Fannia.

Climate-driven mitochondrial selection in lacertid lizards

The mitochondrion, which is an intracellular organelle responsible for most of the energy-producing pathways, can have its genome targeted for climate-driven selection. However, climate-driven mitochondrial selection remains a sparsely studied area in reptiles. Here, we reported the complete mitochondrial genome sequence of a lacertid lizard (Takydromus intermedius) and used mitogenomes from 54 species of lacertid lizards to study their phylogenetic relationships and to identify the mitochondrial genes under positive selection by climate. The length of the complete mitochondrial genome sequence of T. intermedius was 17,713 bp, which was within the range of lengths (17,224-18,943) ever reported for Takydromus species. The arrangement of mitochondrial genes in T. intermedius was the same as in other congeneric species. The 54 lacertid species could be divided into three geographically and climatically different clades. We identified three mitochondrial genes (ATP6, ATP8, and ND3) under positive selection by climate, and found that isothermality, temperature seasonality, precipitation of wettest month, and precipitation seasonality were the most important climatic variables contributing to the gene selection.

The extraordinary rearrangement of mitochondrial genome of the wheat pest, Aptinothrips stylifer and the mitochondrial phylogeny of Thripidae (Thysanoptera)

The mitochondrial gene order in Thysanoptera is notably distinct and highly rearranged, with each species exhibiting its own unique arrangement. To elucidate the relationship between gene rearrangements and phylogeny, the complete mitochondrial genome (mitogenome) of the wheat pest, Aptinothrips stylifer, was sequenced and assembled, spanning a total length of 16,033 bp. Compared with the ancestral arthropod mitogenome, significant rearrangement differences were evident in A. stylifer, whereas the gene order between A. stylifer and Anaphothrips obscurus was similar. Phylogenetic trees were reconstructed based on all 13 protein-coding gene sequences using Bayesian inference and maximum-likelihood methods, both yielding similar topological structures. Notably, A. stylifer was robustly clustered with A. obscurus, affirming its classification within Anaphothrips genus group. This exemplifies the potential correlation between gene rearrangements and phylogeny in the Thripidae family. Additionally, the mitogenome of A. stylifer exhibited several atypical features, including: (1) Three putative control regions (CRs) in close proximity, with CR2 and CR3 displaying partial similarity, and CR1 differing in base composition; (2) Two transfer RNAs (tRNAs), trnS1 and trnV, lacking the DHU arm; (3) Two ribosomal RNA (rRNA) genes inverted and positioned distant from each other; (4) Negative AT and GC skew (AT skew = -0.001, GC skew = -0.077); (5) One transposition (nad6), one inverse transposition (trnQ), four inversions (trnF, trnH, trnC, and gene block nad1-trnL1-rrnL-trnV-rrnS), and four tandem duplication random loss events; and (6) Two protein-coding genes, nad2 and atp8, terminated with an incomplete stop codon "T".

Mitogenomic evolutionary rates in bilateria are influenced by parasitic lifestyle and locomotory capacity

The evidence that parasitic animals exhibit elevated mitogenomic evolutionary rates is inconsistent and limited to Arthropoda. Similarly, the evidence that mitogenomic evolution is faster in species with low locomotory capacity is limited to a handful of animal lineages. We hypothesised that these two variables are associated and that locomotory capacity is a major underlying factor driving the elevated rates in parasites. Here, we study the evolutionary rates of mitogenomes of 10,906 bilaterian species classified according to their locomotory capacity and parasitic/free-living life history. In Bilateria, evolutionary rates were by far the highest in endoparasites, much lower in ectoparasites with reduced locomotory capacity and free-living lineages with low locomotory capacity, followed by parasitoids, ectoparasites with high locomotory capacity, and finally micropredatory and free-living lineages. The life history categorisation (parasitism) explained ≈45%, locomotory capacity categorisation explained ≈39%, and together they explained ≈56% of the total variability in evolutionary rates of mitochondrial protein-coding genes in Bilateria. Our findings suggest that these two variables play major roles in calibrating the mitogenomic molecular clock in bilaterian animals.

The First Complete Mitochondrial Genome of the Genus Pachycondyla (Formicidae, Ponerinae) and Insights into the Phylogeny of Ants

Ants are the standout group among eusocial insects in terms of their exceptional species richness and ecological dominance. The phylogenetic relationships among the group remain elusive. Mitochondrial genome sequences, as a kind of molecular marker, have been widely utilized in the phylogenetic analysis of insects. However, the number of ant mitogenomes published is still very limited. In this study, we utilized next-generation sequencing to determine the complete mitogenome of Pachycondyla annamita (Formicidae, Ponerinae). This is the first mitogenome from the genus Pachycondyla. Two gene rearrangements were identified in the mitogenome, the transposition of trnQ and trnM and the transposition of trnV and rrnS. The secondary structures of tRNAs were predicted. The tRNA genes trnR and trnS1 lacked the dihydrouridine (DHU) arm, and the trnE lacked the TΨC (T) arm. Phylogenetic analyses of the mitochondrial protein-coding genes under maximum likelihood (ML) and Bayesian inference (BI) criteria resulted in conflicting hypotheses. BI analysis using amino acid data with the site-heterogeneous mixture model produced a tree topology congruent with previous studies. The Formicidae was subdivided into two main clades, namely the "poneroid" clade and the "formicoid" clade. A sister group relationship between Myrmicinae and Formicinae was recovered within the "formicoid" clade.

Phylogeny of sea spiders (Arthropoda: Pycnogonida) inferred from mitochondrial genome and 18S ribosomal RNA gene sequences

The phylogeny of sea spiders has been debated for more than a century. Despite several molecular studies in the last twenty years, interfamilial relationships remain uncertain. In the present study, relationships within Pycnogonida are examined in the light of a new dataset composed of 160 mitochondrial genomes (including 152 new sequences) and 130 18S rRNA gene sequences (including 120 new sequences), from 141 sea spider morphospecies representing 26 genera and 9 families. Node congruence between mitochondrial and nuclear markers was analysed to identify the most reliable relationships. We also reanalysed a multilocus dataset previously published and showed that the high percentages of missing data make phylogenetic conclusions difficult and uncertain. Our results support the monophyly of most families currently accepted, except Callipallenidae and Nymphonidae, the monophyly of the superfamilies Ammotheoidea (Ammotheidae + Pallenopsidae), Nymphonoidea (Nymphonidae + Callipallenidae), Phoxichilidioidea (Phoxichilidiidae + Endeidae) and Colossendeoidea (Colossendeidae + Pycnogonidae + Rhynchothoracidae), and the sister-group relationship between Ammotheoidea and Phoxichilidioidea. We discuss the morphological evolution of sea spiders, identifying homoplastic characters and possible synapomorphies. We also discuss the palaeontological and phylogenetic arguments supporting either a radiation of sea spiders prior to Jurassic or a progressive diversification from Ordovician or Cambrian.

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

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

The mitochondrial genome of Heterosentis pseudobagri (Wang & Zhang, 1987) Pichelin & Cribb, 1999 reveals novel aspects of tRNA genes evolution in Acanthocephala

BACKGROUND

Acanthocephala is a clade of obligate endoparasites whose mitochondrial genomes (mitogenomes) and evolution remain relatively poorly understood. Previous studies reported that atp8 is lacking from acanthocephalan mitogenomes, and that tRNA genes often have nonstandard structures. Heterosentis pseudobagri (Arhythmacanthidae) is an acanthocephalan fish endoparasite for which no molecular data are currently available, and biological information is unavailable in the English language. Furthermore, there are currently no mitogenomes available for Arhythmacanthidae.

METHODS

We sequenced its mitogenome and transcriptome, and conducted comparative mitogenomic analyses with almost all available acanthocephalan mitogenomes.

RESULTS

The mitogenome had all genes encoded on the same strand and unique gene order in the dataset. Among the 12 protein-coding genes, several genes were highly divergent and annotated with difficulty. Moreover, several tRNA genes could not be identified automatically, so we had to identify them manually via a detailed comparison with orthologues. As common in acanthocephalans, some tRNAs lacked either the TWC arm or the DHU arm, but in several cases, we annotated tRNA genes only on the basis of the conserved narrow central segment comprising the anticodon, while the flanking 5' and 3' ends did not exhibit any resemblance to orthologues and they could not be folded into a tRNA secondary structure. We corroborated that these are not sequencing artefacts by assembling the mitogenome from transcriptomic data. Although this phenomenon was not observed in previous studies, our comparative analyses revealed the existence of highly divergent tRNAs in multiple acanthocephalan lineages.

CONCLUSIONS

These findings indicate either that multiple tRNA genes are non-functional or that (some) tRNA genes in (some) acanthocephalans might undergo extensive posttranscriptional tRNA processing which restores them to more conventional structures. It is necessary to sequence mitogenomes from yet unrepresented lineages and further explore the unusual patterns of tRNA evolution in Acanthocephala.

Integrative taxonomy of the subfamily Orbioninae (Crustacea: Isopoda) based on mitochondrial and nuclear data with evidence that supports Epicaridea as a suborder

Epicaridea is a group of isopods with high morphological diversity, reduction and loss of characters, and strong sexual dimorphism due to their parasitic lifestyles but their systematics is not well understood. Despite the use of nuclear and mitochondrial genes to test the phylogeny of many invertebrate groups, few molecular data from epicarideans are known, especially from the subfamily Orbioninae. Species in this group are obligate penaeoid shrimp parasites and the lack molecular data has hampered studies on the phylogeny of Orbioninae. To rectify this, mitochondrial and nuclear genes of 9 orbionine species are sequenced here. Compared to the isopod ground pattern, the sequences of orbionines seem to be more plastic near the control region and major translocations are located between rrns and cob. A phylogenetic analysis based on three data sets showed strong support for a monophyletic Orbioninae and that Epicaridea should be accepted at the rank of a suborder within Isopoda. The monophyly of Parapenaeon and Orbione is in doubt based on morphological and molecular data. The genus Parapenaeon is revised and a new genus Aparapenaeon is erected for Parapenaeon japonica and three closely related species.

Complete mitochondrial genome of Idea leuconoe (Lepidoptera: Danaidae) and related phylogenetic analyses

In the present study, we first sequenced and determined the complete mitochondrial genome (mitogenome) of the giant Danaidae butterfly, Idea leuconoe (Lepidoptera: Danaidae). The mitogenome was a typical closed, circular, double-stranded DNA molecule of 15,278 bp length (GenBank accession number: KR815449), similar to the metazoan mitogenomes containing 37 genes and one A + T-rich region. All the protein-coding genes (PCGs) were initiated with a typical ATN codon. Seven genes (COII, ATP6, COIII, nad4, nad4L, cytb, and nad1) adopted the standard ATG start codon, but the remaining six genes were initiated with ATA. All the 13 PCGs harbored complete termination codons (TAA). The overlap nucleotides ATGATAA were conserved for the ATP8/ATP6 gene. The largest intergenic spacer was located between trnGln and nad2, a common finding in Lepidoptera butterflies. All the transfer RNA genes in the I. leuconoe mitogenome could be folded into typical clover-leaf secondary structures, except for trnSer (AGN) that lacked a dihydrouridine arm. The control region with 94.8% A + T content was 444 bp in length and located between rrnS and trnMet. Finally, the phylogenetic relationships obtained using the maximum likelihood and Bayesian methods confirmed two well-supported phylogenetic trees of Danaidae, Papilionidae, and Nymphalidae from the order Lepidoptera, which were consistent with the traditional morphological classification. Results provided additional information for butterfly phylogenetic analysis and insights into the evolution of genomes.

Comparative Mitogenomics of Flesh Flies: Implications for Phylogeny

Flesh flies (Diptera: Sarcophagidae) represent a rapid radiation belonging to the Calyptratae. With more than 3000 known species, they are extraordinarily diverse in terms of their breeding habits and are therefore of particular importance in human and veterinary medicine, forensics, and ecology. To better comprehend the phylogenetic relationships and evolutionary characteristics of the Sarcophagidae, we sequenced the complete mitochondrial genomes of five species of flesh flies and performed mitogenomic comparisons amongst the three subfamilies. The mitochondrial genomes match the hypothetical condition of the insect ancestor in terms of gene content and gene arrangement. The evolutionary rates of the subfamilies of Sarcophagidae differ significantly, with Miltogramminae exhibiting a higher rate than the other two subfamilies. The monophyly of the Sarcophagidae and each subfamily is strongly supported by phylogenetic analysis, with the subfamily-level relationship inferred as (Sarcophaginae, (Miltogramminae, Paramacronychiinae)). This study suggests that phylogenetic analysis based on mitochondrial genomes may not be appropriate for rapidly evolving groups such as Miltogramminae and that the third-codon positions could play a considerable role in reconstructing the phylogeny of Sarcophagidae. The protein-coding genes ND2 and ND6 have the potential to be employed as DNA markers for species identification and delimitation in flesh flies.

The complete mitochondrial genome of Platygaster robiniae (Hymenoptera: Platygastridae): A novel tRNA secondary structure, gene rearrangements and phylogenetic implications

Platygaster robiniae is economically important as a highly specific parasitoid of the invasive pest Obolodiplosis robiniae which was introduced into the Euro-Asia region in the last decade. Despite being a critical and specific parasitoid of the invasive pest O. robiniae and its use as an effective biocontrol agent, the absence of sequence information from P. robiniae have limited its genetic applications for pest management in forests. Mitochondrial (mt) genomes generally contain abundant nucleotide information and thus are helpful for understanding species history. Here, we sequenced the complete mt genome of P. robiniae using next generation sequencing, and annotated 13 protein-coding, 22 tRNA, and 2 rRNA genes and a 702 bp noncoding region. Comparative analysis indicated that this mt genome has a normal A + T content and codons use, however possessed both the expected and unique rearrangements. Ten tRNAs at four gene blocks COII-ATP8, COIII-ND3, ND3-ND5 and the A + T-rich region-ND2 were rearranged, including gene shuffles, transpositions and inversions. Notably, two genes tRNA^Ser(UCN) and tRNA^Leu(CUN) had undergone long-range inversions, which is the first record of this rearrangement type in the superfamily Platygastroidea. The D-loops of both tRNA^Ile and tRNA^Leu(CUN) were absent from the tRNA secondary structure, which has not been reported from hymenopteran previously. Phylogenetic analysis based with the maximum likelihood and Bayesian methods showed that P. robiniae grouped with other species of Platygastridae, and that the superfamily Platygastridea is sister to the other Proctotrupomorpha superfamilies. Our tree strongly supports the monophyly of the five superfamilies of Proctotrupomorpha. This study discovered some unique characters of P. robiniae, and contributes to our understanding of genome rearrangements in the order Hymenoptera.

•

Assemble the complete mitochondria genome of the obligate parasitoid (Platygaster robiniae) of Obolodiplosis robiniae.

•

Discover some unique tRNA secondary structures.

•

Discover some unique rearrangements of Platygaster robiniae and Platygastroidea.

•

The abundance rearrangements in the order Hymenoptera.

Complete Mitochondrial Genome of Malenka flexura (Plecoptera: Nemouridae) and Phylogenetic Analysis

The genus-level relationships within the subfamily Amphinemurinae have been controversial, although attempts have been made based on morphology and limited molecular data. With the establishment of two new genera, the phylogenetic relationships within Amphinemurinae should be re-examined. In this study, the complete mitochondrial genome (mitogenome) of Malenka flexura of the genus Malenka was firstly sequenced and analyzed. The phylogeny of Amphinemurinae was also reconstructed using 13 proteincoding genes (PCGs) from previously published stoneflies. This mitogenome was 15,744 bp long and encoded the typical 37 genes, as well as a putative control region. The gene arrangement of M. flexura mitogenome is identical with the putative ancestral mitogenome in Drosophila yakuba. Most PCGs used standard ATN as start codons and TAA/TAG as termination codons. All tRNA genes exhibited the typical cloverleaf secondary structure, except for tRNA^Ser(AGN), whose dihydrouridine (DHU) arm was lacking. Some structural elements in the control region were founded, such as tandem repeat regions, stemloop structures, polyN stretch and microsatellite structure, etc. Phylogenetic analyses of sequenced Amphinemurinae mitogenomes unsupported the sister relationship of Amphinemura and Malenka. Finally, the phylogenetic analyses inferred a relationship within Amphinemurinae: Amphinemura + (Malenka + (Protonemura + (Indonemoura + (Sphaeronemoura + Mesonemoura)))).

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda)

Background

Within the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes.

Results

P. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture.

Conclusions

Longidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08607-4.

Describing novel mitochondrial genomes of Antarctic amphipods

To date, only one mitogenome from an Antarctic amphipod has been published. Here, novel complete mitochondrial genomes (mitogenomes) of two morphospecies are assembled, namely, Charcotia amundseni and Eusirus giganteus. For the latter species, we have assembled two mitogenomes from different genetic clades of this species. The lengths of Eusirus and Charcotia mitogenomes range from 15,534 to 15,619 base pairs and their mitogenomes are composed of 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and 1 putative control region CR. Some tRNAs display aberrant structures suggesting that minimalization is also ongoing in amphipod mitogenomes. The novel mitogenomes of the two Antarctic species have features distinguishing them from other amphipod mitogenomes such as a lower AT-richness in the whole mitogenomes and a negative GC- skew in both strands of protein coding genes. The genetically most variable mitochondrial regions of amphipods are nad6 and atp8, while cox1 shows low nucleotide diversity among closely and more distantly related species. In comparison to the pancrustacean mitochondrial ground pattern, E. giganteus shows a translocation of the nad1 gene, while cytb and nad6 genes are translocated in C. amundseni. Phylogenetic analysis based on mitogenomes illustrates that Eusirus and Charcotia cluster together with other species belonging to the same amphipod superfamilies. In the absence of reference nuclear genomes, mitogenomes can be useful to develop markers for studying population genetics or evolutionary relationships at higher taxonomic levels.

A Mitochondrial Genome Phylogeny of Cleridae (Coleoptera, Cleroidea)

The predaceous beetle family Cleridae includes a large and widely distributed rapid radiation, which is vital for the ecosystem. Despite its important role, a number of problems remain to be solved regarding the phylogenetic inter-relationships, the timing of divergence, and the mitochondrial biology. Mitochondrial genomes have been widely used to reconstruct phylogenies of various insect groups, but never introduced to Cleridae until now. Here, we generated 18 mitochondrial genomes to address these issues, which are all novel to the family. In addition to phylogenomic analysis, we have leveraged our new sources to study the mitochondrial biology in terms of nucleotide composition, codon usage and substitutional rate, to understand how these vital cellular components may have contributed to the divergence of the Cleridae. Our results recovered Korynetinae sister to the remaining clerids, and the calde of Clerinae+Hydnocerinae is indicated more related to Tillinae. A time-calibrated phylogeny estimated the earliest divergence time of Cleridae was soon after the origin of the family, not later than 160.18 Mya (95% HPD: 158.18–162.07 Mya) during the mid-Jurassic. This is the first mitochondrial genome-based phylogenetic study of the Cleridae that covers nearly all subfamily members, which provides an alternative evidence for reconstructing the phylogenetic relationships.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Saturnia japonica

The complete mitochondrial genome (mitogenome) of Saturnia japonica (Lepidoptera: Saturniidae) was sequenced and annotated. It is a circular molecule of 15, 376 bp, composed of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNA), and an adenine (A) + thymine (T)-rich region. All protein-coding genes (PCGs) are initiated by the ATN codon except for cytochrome c oxidase subunit 1 (cox1) gene that is seemingly initiated by the CGA codon. Except for cox2 and nad4, which were terminated by incomplete stop codon T or TA, the rest were terminated by canonical stop codon TAA. The A + T-rich region is high conservative, including 'ATAGA' motif followed by a 19 bp poly-T stretch, a microsatellite-like element (AT)₉ and also a poly-A element, with a total length of 332 bp. The Asn codon was the most frequently used codon, followed by Ile, Leu2, Lys, Met, Phe, and Tyr, while Cys was the least frequently used codon. Phylogenetic relationships analysis based on the 13 PCGs by using maximum likelihood (ML) and neighbor Joining (NJ) revealed that S. japonica belongs to the Saturniidae family. In this study, the annotation and characteristics of the mitogenome of S. japonica were resolved for the first time, which laid a foundation for species classification and the molecular evolution of Lepidoptera: Saturniidae.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

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

Simple Summary

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

Abstract

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

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.

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea)

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson–Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes

The order Carnivora, which currently includes 296 species classified into 16 families, is distributed across all continents. The phylogeny and the timing of diversification of members of the order are still a matter of debate. Here, complete mitochondrial genomes were analysed to reconstruct the phylogenetic relationships and to estimate divergence times among species of Carnivora. We assembled 51 new mitogenomes from 13 families, and aligned them with available mitogenomes by selecting only those showing more than 1% of nucleotide divergence and excluding those suspected to be of low-quality or from misidentified taxa. Our final alignment included 220 taxa representing 2,442 mitogenomes. Our analyses led to a robust resolution of suprafamilial and intrafamilial relationships. We identified 21 fossil calibration points to estimate a molecular timescale for carnivorans. According to our divergence time estimates, crown carnivorans appeared during or just after the Early Eocene Climatic Optimum; all major groups of Caniformia (Cynoidea/Arctoidea; Ursidae; Musteloidea/Pinnipedia) diverged from each other during the Eocene, while all major groups of Feliformia (Nandiniidae; Feloidea; Viverroidea) diversified more recently during the Oligocene, with a basal divergence of Nandinia at the Eocene/Oligocene transition; intrafamilial divergences occurred during the Miocene, except for the Procyonidae, as Potos separated from other genera during the Oligocene.

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.

Disrupted architecture and fast evolution of the mitochondrial genome of Argeia pugettensis (Isopoda): implications for speciation and fitness

BACKGROUND

Argeia pugettensis is an isopod species that parasitizes other crustaceans. Its huge native geographic range spans the Pacific from China to California, but molecular data are available only for a handful of specimens from North-American populations. We sequenced and characterised the complete mitogenome of a specimen collected in the Yellow Sea.

RESULTS

It exhibited a barcode (cox1) similarity level of only 87-89% with North-American populations, which is unusually low for conspecifics. Its mitogenome is among the largest in isopods (≈16.5 Kbp), mostly due to a large duplicated palindromic genomic segment (2 Kbp) comprising three genes. However, it lost a segment comprising three genes, nad4L-trnP-nad6, and many genes exhibited highly divergent sequences in comparison to isopod orthologues, including numerous mutations, deletions and insertions. Phylogenetic and selection analyses corroborated that this is one of the handful of most rapidly evolving available isopod mitogenomes, and that it evolves under highly relaxed selection constraints (as opposed to positive selection). However, its nuclear 18S gene is highly conserved, which suggests that rapid evolution is limited to its mitochondrial genome. The cox1 sequence analysis indicates that elevated mitogenomic evolutionary rates are not shared by North-American conspecifics, which suggests a breakdown of cox1 barcoding in this species.

CONCLUSIONS

A highly architecturally disrupted mitogenome and decoupling of mitochondrial and nuclear rates would normally be expected to have strong negative impacts on the fitness of the organism, so the existence of this lineage is a puzzling evolutionary question. Additional studies are needed to assess the phylogenetic breadth of this disrupted mitochondrial architecture and its impact on fitness.

Evolution of tRNA gene rearrangement in the mitochondrial genome of ichneumonoid wasps (Hymenoptera: Ichneumonoidea)

Gene rearrangements in the mitochondrial genome (mt genome) are common in certain insect groups and can be an informative character for phylogenetic reconstruction. However, knowledge of the mechanism and biases of gene rearrangement in insect mt genomes is still limited. With an accelerated rate of gene rearrangements, Hymenoptera is an important group for mt genome rearrangements diversity and for understanding the gene rearrangement evolution in mt genomes. Here, we sequenced the complete mt genome of Aphidius gifuensis and analyzed the evolution of tRNA gene rearrangements in the mt genomes of ichneumonoid wasps. Two control regions were detected in A. gifuensis and most of the tRNA rearrangement events occurred around these control regions. tRNA gene rearrangements occurred in almost all of the sequenced mt genomes of Ichneumonoidea and the gene block CR-trnI-trnQ-trnM-ND2-trnW-trnC-trnY was the main hot spot of gene rearrangement. Mapped over the backbone phylogeny of Ichneumonoidea, we found that the inversion and translocation of both trnI and trnM is likely a synapomorphic rearrangement in Braconidae. Our study also demonstrated that the gene block CR-trnI-trnQ-trnM-ND2-trnW-trnC-trnY was important for inferring the gene rearrangement dynamics in Ichneumonoidea.

The conserved mitochondrial genome of the jewel beetle (Coleoptera: Buprestidae) and its phylogenetic implications for the suborder Polyphaga

In this study, we sequenced the mitochondrial (mt) genome of Agrilus mali (Coleoptera: Buprestidae) using next-generation sequencing, and accordingly annotated 13 protein-coding, 22 tRNA, and 2 rRNA genes and a 1458-bp non-coding region. Comparative analysis indicated that the mt genome of A. mali is relatively conserved, with a typical gene content and order identical to those of other coleopterans. However, the newly sequenced mt genome is characterized by a relatively higher A + T content compared with that of other species within the family Buprestidae. Phylogenetic analysis based on Bayesian inference revealed that the evolutionary relationship among the six infraorders of the suborder Polyphaga is (Scirtiformia + (Elateriformia + ((Scarabaeiformia + Staphyliniformia) + (Bostrichiformia + (Cucujiformia))))). However, the topology indicated that the family Buprestidae is a sister group to other Polyphaga infraorders, excluding Scirtiformia as a monophyly, and thus the monophyly of Elateriformia was not supported. This study not only presents the mt genome of a species in the family Buprestidae and a comparative analysis of jewel beetles but also examines the contribution of mt genomes in elucidating phylogenetic relationships within the suborder Polyphaga of Coleoptera.

Major Revisions in Arthropod Phylogeny Through Improved Supermatrix, With Support for Two Possible Waves of Land Invasion by Chelicerates

Deep phylogeny involving arthropod lineages is difficult to recover because the erosion of phylogenetic signals over time leads to unreliable multiple sequence alignment (MSA) and subsequent phylogenetic reconstruction. One way to alleviate the problem is to assemble a large number of gene sequences to compensate for the weakness in each individual gene. Such an approach has led to many robustly supported but contradictory phylogenies. A close examination shows that the supermatrix approach often suffers from two shortcomings. The first is that MSA is rarely checked for reliability and, as will be illustrated, can be poor. The second is that, to alleviate the problem of homoplasy at the third codon position of protein-coding genes due to convergent evolution of nucleotide frequencies, phylogeneticists may remove or degenerate the third codon position but may do it improperly and introduce new biases. We performed extensive reanalysis of one of such "big data" sets to highlight these two problems, and demonstrated the power and benefits of correcting or alleviating these problems. Our results support a new group with Xiphosura and Arachnopulmonata (Tetrapulmonata + Scorpiones) as sister taxa. This favors a new hypothesis in which the ancestor of Xiphosura and the extinct Eurypterida (sea scorpions, of which many later forms lived in brackish or freshwater) returned to the sea after the initial chelicerate invasion of land. Our phylogeny is supported even with the original data but processed with a new "principled" codon degeneration. We also show that removing the 1673 codon sites with both AGN and UCN codons (encoding serine) in our alignment can partially reconcile discrepancies between nucleotide-based and AA-based tree, partly because two sequences, one with AGN and the other with UCN, would be identical at the amino acid level but quite different at the nucleotide level.

Architectural instability, inverted skews and mitochondrial phylogenomics of Isopoda: outgroup choice affects the long-branch attraction artefacts

The majority strand of mitochondrial genomes of crustaceans usually exhibits negative GC skews. Most isopods exhibit an inversed strand asymmetry, believed to be a consequence of an inversion of the replication origin (ROI). Recently, we proposed that an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families resulted in a double-inverted skew (negative GC), and that taxa with homoplastic skews cluster together in phylogenetic analyses (long-branch attraction, LBA). Herein, we further explore these hypotheses, for which we sequenced the mitogenome of Asotana magnifica (Cymothoidae), and tested whether our conclusions were biased by poor taxon sampling and inclusion of outgroups. (1) The new mitogenome also exhibits a double-inverted skew, which supports the hypothesis of an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families. (2) It exhibits a unique gene order, which corroborates that isopods possess exceptionally destabilized mitogenomic architecture. (3) Improved taxonomic sampling failed to resolve skew-driven phylogenetic artefacts. (4) The use of a single outgroup exacerbated the LBA, whereas both the use of a large number of outgroups and complete exclusion of outgroups ameliorated it.

The first mitochondrial genome for Phaudidae (Lepidoptera) with phylogenetic analyses of Zygaenoidea

Phauda flammans Walker belongs to Phaudidae (Lepidoptera), which is a holometabolous and leaf-eating pest that harms trees. So far, there is no mitochondrial (mt) genome reported of Phaudidae. Herein, we sequenced and annotated the complete mt genome of P. flammans representing the first mt genome of Phaudidae and predicted the secondary structures of its RNAs in this study. This mt genome is 15470 bp long consisting of 13 protein coding genes (PCGs), 22 tRNAs, 2 rRNAs and the control region, which are usually conserved in insects. Most PCGs used the standard ATN start codons and complete TAA/TAG termination codons. Almost all of tRNA genes exhibited cloverleaf secondary structures except that the dihydorouridine (DHU) arm of tRNA^Ser(AGN) was absent. The phylogenetic analyses using both Bayesian inference (BI) and maximum likelihood (ML) methods all supported that Phaudidae was a single family being the sister group to Zygaenidae. More mt genomes are needed to better understand the phylogenetic relationships within Zygaenoidea in the future.

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 gene arrangement and phylogeny using mitochondrial genomes in spiders (Arachnida: Araneae)

The complete mitochondrial genome (mitogenome) of Cheiracanthium triviale was sequenced for the first time. The 14,595 bp C. triviale mitogenome contained 37 genes (13 protein coding genes, 2 ribosomal RNAs, 22 transfer RNAs) and one control region. The mitogenome of Dysdera silvatica which was available at NCBI GenBank was annotated. The mitogenome of C. triviale was compared with 43 previously sequenced spider species to observe the gene arrangements, control region and phylogeny. TreeREx analysis identified 19 mitochondrial gene rearrangements (11 transposition, 6 inversion, 2 inverse transposition) in spiders as compared with the putative ancestral gene order and lead to form new gene boundaries: trnQ-trnA, trnA-trnM for Loxosceles similis; nad3-trnS1, trnE-trnL2, trnL2-trnA, trnN-trnF for Agelena silvatica; trnN-trnE, trnE-trnA, trnR-trnF, nad4L-trnW, trnW-trnP for Carrhotus xanthogramma; trnQ-trnW, trnW-trnG, trnG-trnM for Tetragnatha nitens. Our study revealed that the gene rearrangement in spiders with putative ancestor is accelerated in Araneomorphae as compared to Mygalomorphae. Phylogenetic analysis of spiders using mitochondrial sequence data supports the monophyly of two infraorders, and sister relationship of Cheiracanthiidae with Selenopidae and Salticidae. The systematic position of the Cheiracanthium species always a controversial issue as this taxa was placed in different families by different authors.

Disentangling the interplay of positive and negative selection forces that shaped mitochondrial genomes of Gammarus pisinnus and Gammarus lacustris

We hypothesized that the mitogenome of Gammarus lacustris (GL), native to the Qinghai-Tibet Plateau, might exhibit genetic adaptations to the extreme environmental conditions associated with high altitudes (greater than 3000 m). To test this, we also sequenced the mitogenome of Gammarus pisinnus (GP), whose native range is close to the Tibetan plateau, but at a much lower altitude (200-1500 m). The two mitogenomes exhibited conserved mitochondrial architecture, but low identity between genes (55% atp8 to 76.1% cox1). Standard (homogeneous) phylogenetic models resolved Gammaridae as paraphyletic, but 'heterogeneous' CAT-GTR model as monophyletic. In indirect support of our working hypothesis, GL, GP and Gammarus fossarum exhibit evidence of episodic diversifying selection within the studied Gammaroidea dataset. The mitogenome of GL generally evolves under a strong purifying selection, whereas GP evolves under directional (especially pronounced in atp8) and/or relaxed selection. This is surprising, as GP does not inhabit a unique ecological niche compared to other gammarids. We propose that this rapid evolution of the GP mitogenome may be a reflection of its relatively recent speciation and heightened non-adaptive (putatively metabolic rate-driven) mutational pressures. To test these hypotheses, we urge sequencing mitogenomes of remaining Gammarus species populating the same geographical range as GP.

Mitochondrial Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda

The phylogeny of Isopoda, a speciose order of crustaceans, remains unresolved, with different data sets (morphological, nuclear, mitochondrial) often producing starkly incongruent phylogenetic hypotheses. We hypothesized that extreme diversity in their life histories might be causing compositional heterogeneity/heterotachy in their mitochondrial genomes, and compromising the phylogenetic reconstruction. We tested the effects of different data sets (mitochondrial, nuclear, nucleotides, amino acids, concatenated genes, individual genes, gene orders), phylogenetic algorithms (assuming data homogeneity, heterogeneity, and heterotachy), and partitioning; and found that almost all of them produced unique topologies. As we also found that mitogenomes of Asellota and two Cymothoida families (Cymothoidae and Corallanidae) possess inversed base (GC) skew patterns in comparison to other isopods, we concluded that inverted skews cause long-branch attraction phylogenetic artifacts between these taxa. These asymmetrical skews are most likely driven by multiple independent inversions of origin of replication (i.e., nonadaptive mutational pressures). Although the PhyloBayes CAT-GTR algorithm managed to attenuate some of these artifacts (and outperform partitioning), mitochondrial data have limited applicability for reconstructing the phylogeny of Isopoda. Regardless of this, our analyses allowed us to propose solutions to some unresolved phylogenetic debates, and support Asellota are the most likely candidate for the basal isopod branch. As our findings show that architectural rearrangements might produce major compositional biases even on relatively short evolutionary timescales, the implications are that proving the suitability of data via composition skew analyses should be a prerequisite for every study that aims to use mitochondrial data for phylogenetic reconstruction, even among closely related taxa.

Characterization of the mitochondrial genome of an ancient amphipod Halice sp. MT-2017 (Pardaliscidae) from 10,908 m in the Mariana Trench

Small amphipods (Halice sp. MT-2017) with body length <1 cm were collected from the Challenger Deep (~10,920 m below sea level). The divergence time of their lineage was approximately 109 Mya, making this group ancient compared to others under study. The mitochondrial genome of Halice sp. shared the usual gene components of metazoans, comprising 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs). The arrangement of these genes, however, differed greatly from that of other amphipods. Of the 15 genes that were rearranged with respect to the pancrustacean gene pattern, 12 genes (2 PCGs, 2 rRNAs, and 8 tRNAs) were both translocated and strand-reversed. In contrast, the mitochondrial genomes in other amphipods never show so many reordered genes, and in most instances, only tRNAs were involved in strand-reversion-coupled translocation. Other characteristics, including reversed strand nucleotide composition bias, relatively higher composition of non-polar amino acids, and lower evolutionary rate, were also identified. Interestingly, the latter two features were shared with another hadal amphipod, Hirondellea gigas, suggesting their possible associations with the adaptation to deep-sea extreme habitats. Overall, our data provided a useful resource for future studies on the evolutionary and adaptive mechanisms of hadal faunas.

The first complete mitochondrial genome of marigold pest thrips, Neohydatothrips samayunkur (Sericothripinae) and comparative analysis

Complete mitogenomes from the order Thysanoptera are limited to representatives of the subfamily Thripinae. Therefore, in the present study, we sequenced the mitochondrial genome of Neohydatothrips samayunkur (15,295 bp), a member of subfamily Sericothripinae. The genome possesses the canonical 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs) as well as two putative control regions (CRs). The majority strand was 77.42% A + T content, and 22.58% G + C with weakly positive AT skew (0.04) and negative GC skew (-0.03). The majority of PCGs start with ATN codons as observed in other insect mitochondrial genomes. The GCG codon (Alanine) was not used in N. samayunkur. Most tRNAs have the typical cloverleaf secondary structure, however the DHU stem and loop were absent in trnV and trnS1, while the TΨC loop was absent in trnR and trnT. The two putative control regions (CR1 and CR2) show 99% sequence similarity indicated a possible duplication, and shared 57 bp repeats were identified. N. samayunkur showed extensive gene rearrangements, with 11 PCGs, 22 tRNAs, and two rRNAs translocated when compared to the ancestral insect. The gene trnL2 was separated from the 'trnL2-cox2' gene block, which is a conserved, ancestral gene order found in all previously sequenced thrips mitogenomes. Both maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees resulted in similar topologies. The phylogenetic position of N. samayunkur indicates that subfamily Sericothripinae is sister to subfamily Thripinae. More molecular data from different taxonomic groups is needed to understand thrips phylogeny and evolution.

The complete mitochondrial genome of Melon thrips, Thrips palmi (Thripinae): Comparative analysis

The melon thrips, Thrips palmi is a serious pest and vector for plant viruses on a wide range of economically important crops. DNA barcoding evidenced the presence of cryptic diversity in T. palmi and that warrants exhaustive molecular studies. Our present study is on decoding the first complete mitochondrial genome of T. palmi (15,333 bp) through next-generation sequencing (NGS). The T. palmi mt genome contains 37 genes, including 13 Protein coding genes (PCGs), two ribosomal RNA (rRNAs), 22 transfer RNA (tRNAs), and two control regions (CRs). The majority strand of T. palmi revealed 78.29% A+T content, and 21.72% G+C content with positive AT skew (0.09) and negative GC skew (-0.06). The ATN initiation codons were observed in 12 PCGs except for cox1 which have unique start codon (TTG). The relative synonymous codon usage (RSCU) analysis revealed Phe, Leu, Ile, Tyr, Asn, Lys and Met were the most frequently used amino acids in all PCGs. The codon (CGG) which is assigned to Arginine in most insects but absent in T. palmi. The Ka/Ks ratio ranges from 0.078 in cox1 to 0.913 in atp8. We observed the typical cloverleaf secondary structure in most of the tRNA genes with a few exceptions; absence of DHU stem and loop in trnV and trnS, absence of DHU loop in trnE, lack of T-arm and loop in trnN. The T. palmi gene order (GO) was compared with ancestral GO and observed an extensive gene arrangement in PCGs, tRNAs and rRNAs. The cox2 gene was separated from the gene block 'cox2-trnL2' in T. palmi as compared with the other thrips mt genomes, including ancestor GO. Further, the nad1, trnQ, trnC, trnL1, trnV, trnF, rrnS, and rrnL were inversely transpositioned in T. palmi GO. The gene blocks 'trnQ-trnS2-trnD' and 'trnN-trnE-trnS1-trnL1' seems to be genus specific. The T. palmi mt genome contained 24 intergenic spacer regions and 12 overlapping regions. The 62 bp of CR2 shows the similarity with CR1 indicating a possible duplication. The occurrence of multiple CRs in thrips mt genomes seems to be a derived trait which needs further investigation. Although, the study depicted extensive gene rearrangements in T. palmi mt genome, but the negative GC skew reflects only strand asymmetry. Both the ML and BI phylogenetic trees revealed the close relationships of Thrips with Scirtothrips as compared to Frankliniella. Thus, more mt genomes of the diverse thrips species are required to understand the in-depth phylogenetic and evolutionary relationships.

Complete mitochodrial genome of the crab spider Ebrechtella tricuspidata (Araneae: Thomisidae): A novel tRNA rearrangement and phylogenetic implications for Araneae

Mitochondrial genomes are widely used for phylogenetic and phylogeographic analyses among arthropods, but there is a lack of sufficient mitochondrial genome sequence data for spiders. Herein, we sequenced and characterized the complete mitochondrial genome of a crab spider Ebrechtella tricuspidata (Araneae: Thomisidae). The circular mitochondrial genome is 14,352 bp long, including a standard set of 37 genes and an A + T-rich region. Nucleotide composition is highly biased toward A + T nucleotides (77.3%). A novel gene order rearrangement is detected by a tRNA (trnL1) translocation. Tandem repeats are not identified in the A + T-rich region. Most of the tRNAs are greatly reduced in size and cannot be folded into typical cloverleaf-shaped secondary structures. The phylogenetic analysis confirms that the mitochondrial genome sequences are useful in resolving higher-level relationship of Araneae. Overall, our data present in this study will elevate our knowledge on the architecture and evolution of spider mitochondrial genome.

Species delimitation and mitogenome phylogenetics in the subterranean genus Pseudoniphargus (Crustacea: Amphipoda)

The amphi-Atlantic distributions exhibited by many thalassoid stygobiont (obligate subterranean) crustaceans have been explained by fragmentation by plate tectonics of ancestral shallow water marine populations. The amphipod stygobiont genus Pseudoniphargus is distributed across the Mediterranean region but also in the North Atlantic archipelagos of Bermuda, Azores, Madeira and the Canaries. We used species delimitation methods and mitogenome phylogenetic analyses to clarify the species diversity and evolutionary relationships within the genus and timing their diversification. Analyses included samples from the Iberian Peninsula, northern Morocco, the Balearic, Canarian, Azores and Madeira archipelagoes plus Bermuda. In most instances, morphological and molecular-based species delimitation analyses yielded consistent results. Notwithstanding, in a few cases either incipient speciation with no involvement of detectable morphological divergence or species crypticism were the most plausible explanations for the disagreement found between morphological and molecular species delimitations. Phylogenetic analyses based on a robust calibrated mitochondrial tree suggested that Pseudoniphargus lineages have a younger age than for other thalassoid amphipods displaying a disjunct distribution embracing both sides of the Atlantic Ocean. A major split within the family was estimated to occur at the Paleocene, when a lineage from Northern Iberian Peninsula diverged from the rest of pseudoniphargids. Species diversification in the peri-Mediterranean area was deduced to occur in early Miocene to Tortonian times, while in the Atlantic islands it started in the Pliocene. Our results show that the current distribution pattern of Pseudoniphargus resulted from a complex admix of relatively ancient vicariance events and several episodes of long- distance dispersal.

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.

Sequencing of complete mitochondrial genomes confirms synonymization of Hyalomma asiaticum asiaticum and kozlovi, and advances phylogenetic hypotheses for the Ixodidae

Phylogeny of hard ticks (Ixodidae) remains unresolved. Mitochondrial genomes (mitogenomes) are increasingly used to resolve phylogenetic controversies, but remain unavailable for the entire large Hyalomma genus. Hyalomma asiaticum is a parasitic tick distributed throughout the Asia. As a result of great morphological variability, two subspecies have been recognised historically; until a morphological data-based synonymization was proposed. However, this hypothesis was never tested using molecular data. Therefore, objectives of this study were to: 1. sequence the first Hyalomma mitogenome; 2. scrutinise the proposed synonymization using molecular data, i.e. complete mitogenomes of both subspecies: H. a. asiaticum and kozlovi; 3. conduct phylogenomic and comparative analyses of all available Ixodidae mitogenomes. Results corroborate the proposed synonymization: the two mitogenomes are almost identical (99.6%). Genomic features of both mitogenomes are standard for Metastriata; which includes the presence of two control regions and all three "Tick-Box" motifs. Gene order and strand distribution are perfectly conserved for the entire Metastriata group. Suspecting compositional biases, we conducted phylogenetic analyses (29 almost complete mitogenomes) using homogeneous and heterogeneous (CAT) models of substitution. The results were congruent, apart from the deep-level topology of prostriate ticks (Ixodes): the homogeneous model produced a monophyletic Ixodes, but the CAT model produced a paraphyletic Ixodes (and thereby Prostriata), divided into Australasian and non-Australasian clades. This topology implies that all metastriate ticks have evolved from the ancestor of the non-Australian branch of prostriate ticks. Metastriata was divided into three clades: 1. Amblyomminae and Rhipicephalinae (Rhipicephalus, Hyalomma, Dermacentor); 2. Haemaphysalinae and Bothriocrotoninae, plus Amblyomma sphenodonti; 3. Amblyomma elaphense, basal to all Metastriata. We conclude that mitogenomes have the potential to resolve the long-standing debate about the evolutionary history of ticks, but heterogeneous evolutionary models should be used to alleviate the effects of compositional heterogeneity on deep-level relationships.

The mitochondrial genome of the oribatid mite Paraleius leontonychus: new insights into tRNA evolution and phylogenetic relationships in acariform mites

Bilaterian mitochondrial (mt) genomes are circular molecules that typically contain 37 genes. To date, only a single complete mitogenome sequence is available for the species-rich sarcoptiform mite order Oribatida. We sequenced the mitogenome of Paraleius leontonychus, another species of this suborder. It is 14,186 bp long and contains 35 genes, including only 20 tRNAs, lacking tRNA^Gly and tRNA^Tyr. Re-annotation of the mitogenome of Steganacarus magnus increased the number of mt tRNAs for this species to 12. As typical for acariform mites, many tRNAs are highly truncated in both oribatid species. The total number of tRNAs and the number of tRNAs with a complete cloverleaf-like structure in P. leontonychus, however, clearly exceeds the numbers previously reported for Sarcoptiformes. This indicates, contrary to what has been previously assumed, that reduction of tRNAs is not a general characteristic for sarcoptiform mites. Compared to other Sarcoptiformes, the two oribatid species have the least rearranged mt genome with respect to the pattern observed in Limulus polyphemus, a basal arachnid species. Phylogenetic analysis of the newly sequenced mt genome and previously published data on other acariform mites confirms paraphyly of the Oribatida and an origin of the Astigmata within the Oribatida.