Arthropod phylogenetics in light of three novel millipede (myriapoda: diplopoda) mitochondrial genomes with comments on the appropriateness of mitochondrial genome sequence data for inferring deep level relationships

The First Complete Mitochondrial Genome of the Genus Litostrophus: Insights into the Rearrangement and Evolution of Mitochondrial Genomes in Diplopoda

This study presents the complete mitochondrial genome (mitogenome) of Litostrophus scaber, which is the first mitogenome of the genus Litostrophus. The mitogenome is a circular molecule with a length of 15,081 bp. The proportion of adenine and thymine (A + T) was 69.25%. The gene ND4L used TGA as the initiation codon, while the other PCGs utilized ATN (A, T, G, C) as the initiation codons. More than half of the PCGs used T as an incomplete termination codon. The transcription direction of the L. scaber mitogenome matched Spirobolus bungii, in contrast to most millipedes. Novel rearrangements were found in the L. scaber mitogenome: trnQ -trnC and trnL1- trnP underwent short-distance translocations and the gene block rrnS-rrnL-ND1 moved to a position between ND4 and ND5, resulting in the formation of a novel gene order. The phylogenetic analysis showed that L. scaber is most closely related to S. bungii, followed by Narceus magnum. These findings enhance our understanding of the rearrangement and evolution of Diplopoda mitogenomes.

A Rearrangement of the Mitochondrial Genes of Centipedes (Arthropoda, Myriapoda) with a Phylogenetic Analysis

Due to the limitations of taxon sampling and differences in results from the available data, the phylogenetic relationships of the Myriapoda remain contentious. Therefore, we try to reconstruct and analyze the phylogenetic relationships within the Myriapoda by examining mitochondrial genomes (the mitogenome). In this study, typical circular mitogenomes of Mecistocephalus marmoratus and Scolopendra subspinipes were sequenced by Sanger sequencing; they were 15,279 bp and 14,637 bp in length, respectively, and a control region and 37 typical mitochondrial genes were annotated in the sequences. The results showed that all 13 PCGs started with ATN codons and ended with TAR codons or a single T; what is interesting is that the gene orders of M. marmoratus have been extensively rearranged compared with most Myriapoda. Thus, we propose a simple duplication/loss model to explain the extensively rearranged genes of M. marmoratus, hoping to provide insights into mitogenome rearrangement events in Myriapoda. In addition, our mitogenomic phylogenetic analyses showed that the main myriapod groups are monophyletic and supported the combination of the Pauropoda and Diplopoda to form the Dignatha. Within the Chilopoda, we suggest that Scutigeromorpha is a sister group to the Lithobiomorpha, Geophilomorpha, and Scolopendromorpha. We also identified a close relationship between the Lithobiomorpha and Geophilomorpha. The results also indicate that the mitogenome can be used as an effective mechanism to understand the phylogenetic relationships within Myriapoda.

Complete mitochondrial genomes from museum specimens clarify millipede evolution in the Eastern Arc Mountains

The Eastern Arc Mountains in Tanzania represent a hotspot for biological diversity of global importance. The level of endemism is high, and Eastern Arc biodiversity has been studied extensively in vertebrates and invertebrates, including millipedes. However, millipede evolution is vastly understudied at the molecular level. Therefore, we used next-generation ‘shotgun’ sequencing to obtain mitochondrial genome sequences of 26 museum specimens, representing six genera and 12 millipede species found across the Eastern Arc Mountains. Bayesian and maximum likelihood methods yielded consistent topologies with high node support, confirming a high level of congruence between molecular and morphological analyses. The only exception was a Tropostreptus sigmatospinus individual from Zanzibar, which was placed outside an otherwise monophyletic group consisting of mainland individuals of the same assumed species. For two species with a distribution across several mountain blocks (Tropostreptus sigmatospinus and Tropostreptus hamatus), each mountain population represents a distinct monophyletic lineage. In contrast, we also observe that distinct species exist sympatrically in the same montane forests, indicative of older speciation events that are not defined by current forest distribution. Our results are important for understanding speciation mechanisms in montane rain forests and highlight that ethanol-preserved invertebrates exhibit a tremendous potential for genomic analyses.

The Complete Mitochondrial Genome of Spirobolus bungii (Diplopoda, Spirobolidae): The First Sequence for the Genus Spirobolus

Millipedes (Diplopoda) comprise one of the most important groups of large soil arthropods in terrestrial ecosystems; however, their phylogenetic relationships are poorly understood. Herein, the mitochondrial genome (mitogenome) of Spirobolus bungii was sequenced and annotated, which was 14,879 bp in size and included 37 typical mitochondrial genes (13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs)). Most of the 13 PCGs had ATN (AT/A/T/G) as the start codon except for COX1, which used CGA, and most PCGs ended with the T end codon. By comparing the gene arrangements of the mitogenomes among Diplopoda species, rearrangement occurred between and within orders. In contrast to Narceus annularus, the mitogenome genes of S. bungii had consistent orders but were transcribed in completely opposite directions, which was a novel finding in Spirobolidae. Moreover, the phylogenetic relationships within Diplopoda, which were based on the sequences of 13 PCGs, showed that S. bungii was clustered with N. annularus, followed by Abacion magmun. This indicated that there might be a close relationship between Callipodida and Spirobolida. These results could contribute to further studies on the genetics and evolutionary processes of S. bungii and other Diplopoda species.

Novel mitochondrial gene rearrangements pattern in the millipede Polydesmus sp. GZCS-2019 and phylogenetic analysis of the Myriapoda

The subphylum Myriapoda included four extant classes (Chilopoda, Symphyla, Diplopoda, and Pauropoda). Due to the limitation of taxon sampling, the phylogenetic relationships within Myriapoda remained contentious, especially for Diplopoda. Herein, we determined the complete mitochondrial genome of Polydesmus sp. GZCS-2019 (Myriapoda: Polydesmida) and the mitochondrial genomes are circular molecules of 15,036 bp, with all genes encoded on + strand. The A+T content is 66.1%, making the chain asymmetric, and exhibits negative AT-skew (-0.236). Several genes rearrangements were detected and we propose a new rearrangement model: "TD (N\R) L + C" based on the genome-scale duplication + (non-random/random) loss + recombination. Phylogenetic analyses demonstrated that Chilopoda and Symphyla both were monophyletic group, whereas Pauropoda was embedded in Diplopoda to form the Dignatha. Divergence time showed the first split of Myriapoda occurred between the Chilopoda and other classes (Wenlock period of Silurian). We combine phylogenetic analysis, divergence time, and gene arrangement to yield valuable insights into the evolutionary history and classification relationship of Myriapoda and these results support a monophyletic Progoneata and the relationship (Chilopoda + (Symphyla + (Diplopoda + Pauropoda))) within myriapod. Our results help to better explain the gene rearrangement events of the invertebrate mitogenome and lay the foundation for further phylogenetic study of Myriapoda.

The complete mitochondrial genome of the millipede Epanerchodus koreanus Verhoeff, 1937 collected in limestone cave of Korea (Polydesmidae: Polydesmida)

We have determined the second mitochondrial genome of Epanerchodus koreanus Verhoeff, 1937 collected in limestone cave of Korea. The circular mitochondrial genome of E. koreanus is 15,581 bp long. It includes 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes. Its gene order was different from the rest three Polydesmida mitochondrial genomes, resulted from relocation of tRNAs, rRNAs, and ND1. The base composition was AT-biased (75.1%). Phylogenetic trees displayed phylogenetic relationship, which is congruent to previous study, except Sphaerotheriidae sp. clustering with Helminthomorpha.

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.

Complete mitochondrial genome of Sellanucheza jaegeri Golovatch, 2013 by next generation sequencing (Polydesmida: Paradoxosomatidae) and phylogenetic analysis in Diplopoda

The mitogenome of Sellanucheza jaegeri was 15,623 bp long, revealed the same gene order to that of typical Polydesmida. Both the tRNA^Ser(AGN) and tRNA^Ser(UCN) lacked the DHU arms. No tandem repeat was found in two control regions. Phylogenetic analysis indicated that Sphaerotheriida was so antiquity that divided out earlier than others. We supported that Polydesmida had a relatively systematic affinity between Julida and Playtdesmida, and suggested that the interordinal phylogenetic relationships within Diplopoda should be further investigated.

Complete mitochondrial genomes of two flat-backed millipedes by next-generation sequencing (Diplopoda, Polydesmida)

A lack of mitochondrial genome data from myriapods is hampering progress across genetic, systematic, phylogenetic and evolutionary studies. Here, the complete mitochondrial genomes of two millipedes, Asiomorphacoarctata Saussure, 1860 (Diplopoda: Polydesmida: Paradoxosomatidae) and Xystodesmus sp. (Diplopoda: Polydesmida: Xystodesmidae) were assembled with high coverage using Illumina sequencing data. The mitochondrial genomes of the two newly sequenced species are circular molecules of 15,644 bp and 15,791 bp, within which the typical mitochondrial genome complement of 13 protein-coding genes, 22 tRNAs and two ribosomal RNA genes could be identified. The mitochondrial genome of Asiomorphacoarctata is the first complete sequence in the family Paradoxosomatidae (Diplopoda: Polydesmida) and the gene order of the two flat-backed millipedes is novel among known myriapod mitochondrial genomes. Unique translocations have occurred, including inversion of one half of the two genomes with respect to other millipede genomes. Inversion of the entire side of a genome (trnF-nad5-trnH-nad4-nad4L, trnP, nad1-trnL2-trnL1-rrnL-trnV-rrnS, trnQ, trnC and trnY) could constitute a common event in the order Polydesmida. Last, our phylogenetic analyses recovered the monophyletic Progoneata, subphylum Myriapoda and four internal classes.

Phylogeny of Syndermata (syn. Rotifera): Mitochondrial gene order verifies epizoic Seisonidea as sister to endoparasitic Acanthocephala within monophyletic Hemirotifera

A monophyletic origin of endoparasitic thorny-headed worms (Acanthocephala) and wheel-animals (Rotifera) is widely accepted. However, the phylogeny inside the clade, be it called Syndermata or Rotifera, has lacked validation by mitochondrial (mt) data. Herein, we present the first mt genome of the key taxon Seison and report conflicting results of phylogenetic analyses: while mt sequence-based topologies showed monophyletic Lemniscea (Bdelloidea+Acanthocephala), gene order analyses supported monophyly of Pararotatoria (Seisonidea+Acanthocephala) and Hemirotifera (Bdelloidea+Pararotatoria). Sequence-based analyses obviously suffered from substitution saturation, compositional bias, and branch length heterogeneity; however, we observed no compromising effects in gene order analyses. Moreover, gene order-based topologies were robust to changes in coding (genes vs. gene pairs, two-state vs. multistate, aligned vs. non-aligned), tree reconstruction methods, and the treatment of the two monogonont mt genomes. Thus, mt gene order verifies seisonids as sister to acanthocephalans within monophyletic Hemirotifera, while deviating results of sequence-based analyses reflect artificial signal. This conclusion implies that the complex life cycle of extant acanthocephalans evolved from a free-living state, as retained by most monogononts and bdelloids, via an epizoic state with a simple life cycle, as shown by seisonids. Hence, Acanthocephala represent a rare example where ancestral transitional stages have counterparts amongst the closest relatives.

A New Dimension in Documenting New Species: High-Detail Imaging for Myriapod Taxonomy and First 3D Cybertype of a New Millipede Species (Diplopoda, Julida, Julidae)

We review the state-of-the-art approaches currently applied in myriapod taxonomy, and we describe, for the first time, a new species of millipede (Ommatoiulus avatar n. sp., family Julidae) using high-resolution X-ray microtomography (microCT) as a substantive adjunct to traditional morphological examination. We present 3D models of the holotype and paratype specimens and discuss the potential of this non-destructive technique in documenting new species of millipedes and other organisms. The microCT data have been uploaded to an open repository (Dryad) to serve as the first actual millipede cybertypes to be published.

Genome of the Rusty Millipede, Trigoniulus corallinus, Illuminates Diplopod, Myriapod, and Arthropod Evolution

The increasing availability of genomic information from the Arthropoda continues to revolutionize our understanding of the biology of this most diverse animal phylum. However, our sampling of arthropod diversity remains uneven, and key clade such as the Myriapoda are severely underrepresented. Here we present the genome of the cosmopolitanly distributed Rusty Millipede Trigoniulus corallinus, which represents the first diplopod genome to be published, and the second example from the Myriapoda as a whole. This genomic resource contains the majority of core eukaryotic genes (94.3%), and key transcription factor classes that were thought to be lost in the Ecdysozoa. Mitochondrial genome and gene family (transcription factor, Dscam, circadian clock-driving protein, odorant receptor cassette, bioactive compound, and cuticular protein) analyses were also carried out to shed light on their states in the Diplopoda and Myriapoda. The ready availability of T. corallinus recommends it as a new model for evolutionary developmental biology, and the data set described here will be of widespread utility in investigating myriapod and arthropod genomics and evolution.

Evolutionary implications of dipluran hexamerins

Hexamerin, as a member of the highly conserved arthropod hemocyanin superfamily, has been shown to be a good marker for the phylogenetic study of insects. However, few studies have been conducted on hexamerins in basal hexapods. The first Diplura hexamerin CspHex1 was reported only recently (Pick and Burmester, 2009). Remarkably, CspHex1 was suggested to have evolved from hexapod hemocyanin subunit type 2, which is very different from all insect hexamerins originated from hexapod hemocyanin subunit type 1. Does this finding suggest double or even multiple origins of hexamerins in Hexapoda? To find more evidence on the evolution of dipluran hexamerins, eight putative hexamerin gene sequences were obtained from three dipluran species, as were three hemocyanin genes from two collembolan species. Unexpectedly, after adding the new sequences into the phylogenetic analyses, all dipluran hexamerins including CspHex1 grouped together and as sister to the insect hexamerins, with high likelihood and Bayesian support. Our analysis supports a single origin of the hexamerins in Hexapoda, and suggests the close relationship between Diplura and Insecta. In addition, our study indicates that a relatively comprehensive taxa sampling is essential to solve some problems in phylogenetic reconstruction.

Revival of forgotten characters and modern imaging techniques help to produce a robust phylogeny of the Diplopoda (Arthropoda, Myriapoda)

The external and internal anatomy of millipedes (Diplopoda) is poorly known compared to some of the other myriapod and arthropod groups. Due to both language barriers, which hindered the assessment of the character-rich older literature, and non-phylogenetic thinking, our knowledge of morphological characters useful for phylogenetic work diminished over the last decades. Here, a new character matrix with 64 characters, mainly derived from old literature data, is used to reconstruct a phylogeny of Diplopoda. As a tool to further our knowledge about the morphology of the different millipede orders, we show how micro-computer tomography (μCT) can be used to assess and illustrate specific parts of the Platydesmida. With the advent of μCT it is now possible to analyse many taxa and characters in a comparatively short time. A focus is put on potential phylogenetically useful characters. Our results support a Verhoeffian classification of the Diplopoda: Polyxenida + Chilognatha. Pentazonia are the sistergroup to the Helminthomorpha. Colobognatha form the sistergroup to Eugnatha, the latter split into monophyletic Juliformia and Polydesmida + Nematophora.

Ordinal-level phylogenomics of the arthropod class Diplopoda (millipedes) based on an analysis of 221 nuclear protein-coding loci generated using next-generation sequence analyses

BACKGROUND

The ancient and diverse, yet understudied arthropod class Diplopoda, the millipedes, has a muddled taxonomic history. Despite having a cosmopolitan distribution and a number of unique and interesting characteristics, the group has received relatively little attention; interest in millipede systematics is low compared to taxa of comparable diversity. The existing classification of the group comprises 16 orders. Past attempts to reconstruct millipede phylogenies have suffered from a paucity of characters and included too few taxa to confidently resolve relationships and make formal nomenclatural changes. Herein, we reconstruct an ordinal-level phylogeny for the class Diplopoda using the largest character set ever assembled for the group.

METHODS

Transcriptomic sequences were obtained from exemplar taxa representing much of the diversity of millipede orders using second-generation (i.e., next-generation or high-throughput) sequencing. These data were subject to rigorous orthology selection and phylogenetic dataset optimization and then used to reconstruct phylogenies employing Bayesian inference and maximum likelihood optimality criteria. Ancestral reconstructions of sperm transfer appendage development (gonopods), presence of lateral defense secretion pores (ozopores), and presence of spinnerets were considered. The timings of major millipede lineage divergence points were estimated.

RESULTS

The resulting phylogeny differed from the existing classifications in a number of fundamental ways. Our phylogeny includes a grouping that has never been described (Juliformia+Merocheta+Stemmiulida), and the ancestral reconstructions suggest caution with respect to using spinnerets as a unifying characteristic for the Nematophora. Our results are shown to have significantly stronger support than previous hypotheses given our data. Our efforts represent the first step toward obtaining a well-supported and robust phylogeny of the Diplopoda that can be used to answer many questions concerning the evolution of this ancient and diverse animal group.