The mitochondrial genome of Sinentomon erythranum(Arthropoda: Hexapoda: Protura): an example of highly divergent evolution
- Wan-Jun Chen†1,
- Yun Bu†1,
- Antonio Carapelli2,
- Romano Dallai2,
- Sheng Li1,
- Wen-Ying Yin1 and
- Yun-Xia Luan1Email author
© Chen et al; licensee BioMed Central Ltd. 2011
Received: 28 April 2011
Accepted: 27 August 2011
Published: 27 August 2011
The phylogenetic position of the Protura, traditionally considered the most basal hexapod group, is disputed because it has many unique morphological characters compared with other hexapods. Although mitochondrial genome information has been used extensively in phylogenetic studies, such information is not available for the Protura. This has impeded phylogenetic studies on this taxon, as well as the evolution of the arthropod mitochondrial genome.
In this study, the mitochondrial genome of Sinentomon erythranum was sequenced, as the first proturan species to be reported. The genome contains a number of special features that differ from those of other hexapods and arthropods. As a very small arthropod mitochondrial genome, its 14,491 nucleotides encode 37 typical mitochondrial genes. Compared with other metazoan mtDNA, it has the most biased nucleotide composition with T = 52.4%, an extreme and reversed AT-skew of -0.351 and a GC-skew of 0.350. Two tandemly repeated regions occur in the A+T-rich region, and both could form stable stem-loop structures. Eighteen of the 22 tRNAs are greatly reduced in size with truncated secondary structures. The gene order is novel among available arthropod mitochondrial genomes. Rearrangements have involved in not only small tRNA genes, but also PCGs (protein-coding genes) and ribosome RNA genes. A large block of genes has experienced inversion and another nearby block has been reshuffled, which can be explained by the tandem duplication and random loss model. The most remarkable finding is that trnL2(UUR) is not located between cox1 and cox2 as observed in most hexapod and crustacean groups, but is between rrnL and nad1 as in the ancestral arthropod ground pattern. The "cox1-cox2" pattern was further confirmed in three more representative proturan species. The phylogenetic analyses based on the amino acid sequences of 13 mitochondrial PCGs suggest S. erythranum failed to group with other hexapod groups.
The mitochondrial genome of S. erythranum shows many different features from other hexapod and arthropod mitochondrial genomes. It underwent highly divergent evolution. The "cox1-cox2" pattern probably represents the ancestral state for all proturan mitogenomes, and suggests a long evolutionary history for the Protura.
The Protura is a group of mysterious soil-dwelling micro-arthropods (usually 0.5-2.0 mm in length), first described by Silvestri in 1907 . Traditionally, it was regarded as a basal hexapod group, but it owns many unique and primitive morphological characteristics compared with other hexapods. For example, they lack antennae and wings, the foretarsus are enlarged with many sensilla serving the role of antennae, eyes and tentorium are absent, they have anamorphic post-embryonic development, and they have 12 abdominal segments (instead of 11) . The proturan spermatozoan has a variable number of doublet microtubules (9-16), with no accessory or central microtubules. It is different from those of other hexapods, but similar to the sperm of sea spider (Arthropoda: Pycnogonida). This probably reflects a high diversification rate, or a lengthy evolution [3–5]. Historically, there were many controversies about the relationship of proturans to other hexapods, and their evolutionary position in the Arthropoda [2, 3, 6–9]. This is because proturans are understudied, being so small and rare, making them difficult to collect, identify, culture and experiment on [2, 10, 11].
The higher-level phylogeny of the major arthropod groups (Chelicerata, Myriapoda, Crustacea and Hexapoda) continues to be a matter of debate despite extensive research based on phylogenetic analysis and genetic data [12–14]. Almost all molecular analyses strongly support the Pancrustacea hypothesis: crustaceans, instead of myriapods, are the closest relatives of the hexapods [15–18]. The Hexapoda (Insecta s. lat.), which includes four groups, Protura, Collembola, Diplura and Insecta (Insecta s. str.), was traditionally considered a monophyletic lineage based on the synapomorphies of body segments, six legs on the thorax, and adaptation to the terrestrial environment. The monophyly of the Insecta has been well established by morphological and molecular studies [8, 10, 17, 18], but the monophyly of the Hexapoda is less certain [17, 19]. Three basal hexapod groups (Protura, Collembola and Diplura) show many different features from insects according to morphology [10, 20] and ultrastructure of spermatozoa . The mitogenomic data of basal hexapod collembolans and diplurans reject the monophyly of Hexapoda, and suggest that some crustaceans are more closely related to the Insecta than Collembola and Diplura [17, 19, 21]. However, recent studies based on EST data and nuclear genes (18S and 28S ribosomal RNA genes, nuclear protein-coding sequences) support the monophyly of the Hexapoda [12, 13, 18].
The arthropod mitochondrial genome is a single circular DNA molecule encoding 13 proteins, 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one A+T-rich region for the control of replication and transcription of the mtDNA. It is used extensively for studying phylogenetic relationships at various taxonomic levels. Unlike nuclear molecular markers, mtDNA is of maternal inheritance, and does not experience intermolecular genetic recombination. In addition, the mitochondrial gene order can provide additional phylogenetic information, since rearrangements appear to be generally rare events, and most mitochondrial gene arrangements often remain unchanged over a long evolutionary period . Mitogenomic data also strongly support the Pancrustacea hypothesis [14, 17, 23], especially with the evidence of the gene order [16, 24]. The gene trnL2 (UUR) is located between rrnL and nad1 in the ancestral arthropod ground pattern, but is translocated to the position between cox1 and cox2 in Pancrustacea . It has been considered a distinctive synapomorphic character for crustaceans and hexapods. The mitochondrial genomes of basal hexapod Collembola  and Diplura  also agree with the "cox1-trnL2-cox2" pattern. So far, no mitochondrial genome information is available for the Protura. This has impeded comprehensive discussions on the evolution of the arthropod mitochondrial genome, and the validity of using mtDNA to study the phylogeny of the Hexapoda [27–29].
In this study, we sequenced the complete mitochondrial genome of Sinentomon erythranum (Protura: Sinentomata: Sinentomidae), to describe the molecular features of the proturan mitochondrial genome, to judge how these evolved, and to see if it has any phylogenetic information, which may help resolve the discrepancy on the monophyly of the Hexapoda between mitochondrial and nuclear DNA markers.
Results and Discussion
General description of the mitochondrial genome of S. erythranum
Annotation table for the mitochondrial genome of S. erythranum
Figure 2B shows the nucleotide composition, AT-skew and GC-skew for each of the 13 PCGs and two rRNA genes of the mitochondrion of S. erythranum. Cox1 has the lowest AT content (70.2%) and atp8 has the highest AT content (85.6%). The AT content of these 15 genes does not fluctuate far from the overall average AT content (77.6%). Nad3 has the most negative AT-skew (-0.685), and nad4 and rrnS share the least extreme AT-skew (-0.204). The AT-skew values of the adjacent genes nad5, nad4 and nad4L are less extreme than in other adjacent genes, and all three of these genes are encoded by the minority strand, so it seems that some constraints shaped the genome that evolved under a strong directional mutation pressure (Figure 2B) .
Gene rearrangements and possible evolutionary mechanisms
Rearrangements 1 and 2: the translocation of trnF may be an independent event, and this kind of minor rearrangement is very common in mtDNA [42, 43]. The trnP changed its coding strand from N to J during its "long range" translocation, and this situation is rarely reported.
Rearrangements 3 and 4: The tandem duplication and random loss (TDRL) model is a popular hypothesis for explaining many mtDNA gene rearrangements [44–46]. Here, it can readily explain the reshuffling of tRNAs in the region from trnI to trnC (rearrangement 4 in Figure 5), although it does not explain the gene inversion (rearrangement 3 in Figure 5). For that inversion, the implication is strong that the gene block "rrnS-V-rrnL-trnL2-trnL1-nad1" was locally reversed as a whole. Gene inversions are probably the result of intra-molecular recombination, which can not only rearrange parts of the genome but also invert them at the same time. In the mitogenomic sequence of S. erythranum, both gene relocation and inversion must have occurred, although it is uncertain which of these two processes dominated. Here, we have some new thoughts. For the TDRL model, gene duplication is necessary, which can be achieved by replication slippage in single stranded templates. At the same time, a loop must be produced by slippage, so it is possible for the loop to perform intra-molecular recombination simultaneously . Namely, the reshuffling of tRNAs and local inversion of a gene block may happen together in a stepwise rearrangement process. We further checked available mitochondrial genomes, and found that recombination involving PCGs has rarely occurred in hexapods, except in some lice whose mitochondrial genomes were extensively shuffled .
Rearrangement 5: it is not easy to explain the translocation of the A+T-rich region. There is a hint of an orientation change of replication due to the nucleotide-bias change from the majority type (AT-skew and GC-skew) (Figure 2A), but it is hard to explain it as a consequence of the inversion of gene block "rrnS-V-rrnL-trnL2-trnL1-nad1".
Position of trnL2(UUR) and its phylogenetic implications
The Protura has three groups: Acerentomata, Sinentomata and Eosentomata. Besides S. erythranum, a member of the Sinentomata, we also sequenced the cox1/cox2 region (about 1.4 kb) from Baculentulus tianmushanensis of Acerentomata (GenBank accession HQ416715), Eosentomon nivocolum of Eosentomata (GenBank accession HQ416716), and Zhongguohentomon piligeroum of Eosentomata (GenBank accession HQ416714). They all agree with the cox1-cox2 pattern and have no intervening trnL2. In addition, cox1 is the exact neighbor to cox2 with no nucleotide between them in S. erythranum, B. tianmushanensis and E. nivocolum, and only four intergenic nucleotides in Z. piligeroum. Therefore, based on the available data, we believe it is more reasonable to conclude that the ancestral state is the cox1-cox2 pattern for all proturan mtDNAs.
The "cox1-trnL2-cox2" pattern occurs in almost all hexapods. We compared all published data of arthropod mitogenomes (available until January 16, 2011), and found only eight of 226 mtDNAs of Insecta are not consistent with this pattern (Figure 6 and Additional File 2), but they are clearly secondary mtDNA rearrangements or with multiple trnL2 copies. Five of them are from the Hemiptera, three parasitic lice from the Phthiraptera (Bothriometopus macrocnemis, C. bidentatus compar and Heterodoxus macropus) [52, 53], one bark louse from the Psocoptera (Lepidopsocid sp. RS-2001) and one species from the Thysanoptera (Thrips imaginis) . Their mitochondrial gene arrangements are reshuffled rigorously. The other three exceptions are from the Hymenoptera (Vanhornia eucnemidarum, Abispa ephippium and Diadegma semiclausum) . It was noticed that in Hymenoptera, tRNA rearrangements (termed minor rearrangements) are very common, especially in the hot-spot areas . In Abispa ephippium, trnL2 has four copies, but is still located between cox1 and cox2 . However, most hemipteran and hymenopteran mtDNAs are still consistent with the cox1-trnL2-cox2 pattern. In Crustacea, only nine of 60 mitochondrial genomes are not consistent with the cox1-trnL2-cox2 pattern (Additional File 2). In addition, only seven of 53 mitochondrial genomes from the Chelicerata are not consistent with the cox1-cox2 pattern (Additional File 2), and all eight reported mitochondrial genomes from the Myriapoda are consistent with the cox1-cox2 pattern (Figure 6).
These statistics reflect the fact that translocation of trnL2 out of the cox1/cox2 junction has rarely happened within Pancrustacea lineage, and no case of the cox1-trnL2-cox2 pattern was detected within Myriapoda and Chelicerata lineages, whose trnL2 tends to stay between rrnL and nad1. This information leads to a single plausible scenario of the ancestral state being cox1-trnL2-cox2 in the Hexapoda, but the proturan mitochondrial genomes likely retain the ancestral state of the Arthropoda, the cox1-cox2 pattern. This seems to cast new doubt on the monophyly of Hexapoda. The Protura probably has a very ancient origin and a long evolutionary history, with distant affinity to other hexapods, evolving even earlier than other pancrustaceans. However, we cannot exclude the possibility of the secondary reversion to the primitive arthropod condition in the proturan ancestor since our gene sequence is so highly divergent. In this case, the mtDNA of S. erythranum provides a remarkable example of secondary reversion.
Phylogenetic position of Protura
Since the position of trnL2 cast doubt on the relationship between the Protura and other hexapods, it is important to verify it with a phylogenetic tree. As revealed in Figure 2A, the base composition of S. erythranum is so different from that of most arthropod mitochondrial genomes, long-branch attraction (LBA) can be expected. Translating the PCGs into amino acid sequences is an effective method of dealing with the problem caused by base compositional heterogeneity in tree reconstruction [14, 17, 56], so we performed all phylogenetic analyses on conceptually translated amino acid data of 13 mitochondrial PCGs using maximum likelihood and Bayesian inference methods.
In our trees (Figure 7), the clade of Diplura and Collembola is sister to Insecta, although the bootstrap value is relatively low. It is different from previous studies based on mitochondrial gene sequences of diplurans and collembolans, which suggested that some crustaceans are more closely related to Insecta than Collembola and Diplura . More arthropod taxa are needed to further discuss this problem.
The unusual long-branch length indicates that the S. erythranum mitochondrial genomes are evolving rapidly. The population of soil-dwelling proturans is usually very small. Mutations may accumulate faster in such organisms due to the slow rate of gene flow. This also seems true for nematodes, parasitic lice and mites, in which high levels of genome diversity are commonly detected. The study on the mitochondrial genome of two diplurans also reveals that high genetic divergence existed in the morphologically uniform taxa .
Whether the Protura is a real hexapod group or not has been debated for a long time . The Protura have many unique morphological characters compared with other hexapods: 1) they have no eyes and no antennae; 2) they have abdominal legs on abdominal segments 1-3; 3) they have no caudal cerci but have a telson tail, which is common in crustaceans but absent in other hexapods [1–3]; 4) the axoneme of flagellated spermatozoa lacks central microtubules, which is similar to the condition in pycnogonid spermatozoa ; 5) the serosa (embryonic membrane) of proturans retains the ability to differentiate into a tergum or definitive dorsal closure during embryonic development, which is similar to crustaceans and myriapods, but different from other hexapods. Based on information from embryonic development, Machida (2006) proposed that the Protura may have a much longer evolutionary history than previously thought . However, a few recent studies based on EST data and rRNA genes have presented relatively robust evidence supporting the monophyly of Hexapoda and Pancrustacea (although only one proturan species was included in these studies) [12, 18].
Although the mitochondrial genome sequence of S. erythranum is unique, with little phylogenetic affinity to the insects, we cannot equate this to the evolutionary history of the Protura. Mitochondrial genome data alone are not enough to unambiguously resolve the relationships of Protura, Diplura, Collembola and Insecta. It is necessary to understand the limits and applicability of these data . Our sequence data showed many unique molecular features, which can provide valuable information for studying problems of mitochondrial genome evolution, for example, the mechanisms of mitochondrial gene rearrangements, truncation of tRNA secondary structures, and nucleotide frequency bias. Understanding these fundamental biology problems should be helpful in phylogenetic analyses when using mitochondrial genomic data.
This is the first report of a complete mitochondrial genome from the Protura. With highly divergent evolution, their mtDNA has many different features to that of other hexapods, including nucleotide-frequency bias, gene order, and tRNA secondary structure. Therefore, it is a valuable example to study the mechanism of mitochondrial gene evolution and rearrangement in the Arthropoda.
Our study suggests that proturan mtDNAs do not agree with the "cox1-trnL2-cox2" pattern, which was thought to be an important character shared by hexapod and crustacean groups. It may be a result of secondary reversion due to extensive rapid and divergent evolution, but also may suggests that the Protura have a long evolutionary history, and do not have a close affinity to hexapods and crustaceans. S. erythranum did not group with other hexapods in our phylogenetic trees, and its extreme long-branch implies that its mtDNA underwent highly divergent evolution. More evidence is needed to verify this hypothesis and to solve the conflict between the studies on mitochondrial and nuclear gene markers.
mtDNA sequencing of S. erythranum
Specimens of S. erythranum were collected from Tianping Mountain (Jiangsu Province, China). The total DNA of one individual was extracted with the commercial kit Wizard SV Genomic Purification System (Promega), and then was used as the template for PCR amplifications. Initially, two small fragments of cox1 and cob were amplified using two universal primer pairs of LCO1490/HCO2198  and CobF424/CobR876 , respectively, and the PCR products were sequenced directly by the amplification primers. Four primers were designed according to these obtained sequences for two long PCR amplifications encompassing the cox1/cob (~9 kb) and cob/cox1 (~6 kb) fragments, respectively. These primers were SI-C1-J320 (CTGGTTGAACTGTTTATCCTCCTC)/SI-Cb-N239 (ATAAGGATGAAAACTAACCCTATCA), and SI-Cb-J181 (GTTCTTCTAATCCTTTAGGAGTTGG)/SI-C1-N343 (GAGGAGGATAAACAGTTCAACCAG). Long PCRs were generated with LA Taq (Takara, Dalian, China) under the following two-step conditions: 35 cycles of 96°C for 2 min and 68°C for 10 min, followed by incubation at 68°C for 10 min. The 9 kb and 6 kb products were mixed together after gel-purification, and then sequenced with the shotgun sequencing approach as described by Masta and Boore (2004) . The sequencing service was from Shanghai Majorbio Biotech Co., Ltd. Two contigs were assembled by Phred/Phrap [59, 60] from the shotgun sequencing readings, guaranteed to have 10 times coverage for both contigs. More specific primers were designed for PCR amplifications to bridge two remaining gaps (primers available on request). All PCR products were then cloned and then sequenced by an ABI 3730 automated DNA sequencer. A consensus sequence was assembled from all the contigs using Seqman in the DNAStar software package (DNASTAR Inc., Madison, WI) .
Gene annotation and secondary structure prediction
The sequence was submitted in Fasta format to the web-based software DOGMA (Dual Organellar Genome Annotator)  for primary annotation. BLAST searches were done on NCBI Blast Entrez databases to ensure the identity of PCGs and rRNA genes. To identify the tRNA genes in the genome, we used the annotation obtained by DOGMA (with the COVE threshold for tRNAs set to 7(low)), and further used tRNAscan-SE via the web interface and the "Nematode Mito" settings for the COVE program . The ARWEN (version 1.2) program was also used by the web interface with the "mtmam" option switched off . Finally, the tRNAs were determined by comparing the secondary structures suggested by these different programs. Tandemly repetitive sequences in the A+T-rich region were determined both manually and by using the Tandem Repeats Finder . The putative minimum-free-energy structures of TRRs were given by RNAfold WebServer in the Vienna RNA Websuite .
Sequence determination of cox1/cox2junction region
In order to find if trnL2 lay outside of cox1 and cox2, not only in the Sinentomata but also in the other proturan groups, we amplified and sequenced the cox1/cox2 junction (about 1.4 kb) of B. tianmushanensis (Acerentomata: Berberentomidae), E. nivocolum (Eosentomata: Eosentomidae) and Z. piligeroum (Eosentomata: Eosentomidae) using the universal primer pair C1-HCO-J and C2-B-3665 . We followed the above-mentioned methods to annotate these genes.
Statistical comparison of strand asymmetry and of trnL2positions of arthropod mtDNAs
We retrieved the nucleotide sequences and DNA compositions for all 359 published arthropod mtDNAs (before January 16, 2011) from the Mitome database  or NCBI Organelle Genome Resources. Strand asymmetry represents strand compositional bias, usually reflected by the AT skew = (A-T)/(A+T) and GC-skew = (G-C)/(G+C) [32, 68].
We further checked the position of trnL2 in all 359 available arthropod mtDNAs. For the pancrustacean groups, we checked whether each mtDNA agreed with the typical patterns of cox1-trnL2-cox2 and rrnL-trnL1-nad1; then, we did the same for the other arthropods, the myriapods and chelicerates, which typically have the different pattern of cox1-cox2 and rrnL-trnL1-trnL2-nad1 .
First, we choose 24 Panarthropoda representatives (Additional File 3) for phylogenetic tree construction based on previous studies [14, 17], including three groups with the similar base composition to S. erythranum (negative AT-skew and positive GC-skew, Additional File 4), in order to see if S. erythranum will group with them because of LBA. Then, we reconstructed the phylogenetic trees after removing these three taxa, focusing on the relationship of S. erythranum and other hexapods. The onychophoran Opisthopatus cinctipes was defined as the outgroup in our analyses.
The nucleotide sequences of each PCG were retro-aligned using DAMBE, version 5.1.1 . The 13 amino acid data were concatenated as an alignment of 3819 positions after individually aligned, and then, 2520 aligned characters for 24 taxa and 2616 aligned characters for 21 taxa were retained respectively after Gblocks screening with default settings . The best model "mtREV24+G+I+F" was selected using MEGA 5.0 . We carried out ML searches with RAxML through the web portal http://phylobench.vital-it.ch/raxml-bb/index.php. Bayesian analysis was performed using MrBayes (version 3.1.2), with mtRev+I+G model . Four Markov chains were run for 1,000,000 generations, and sampled every 100 generations to yield a posterior probability distribution of 10,000 trees. The first 2,000 trees were discarded as burn-in. The standard deviation of split frequencies was lower than 0.01 in 21 taxa dataset analysis, but we failed to obtain a meaningful convergence for the 24 taxa dataset.
We are grateful to Jon Mallatt from Washington State University and Yonggang Yao from the Kunming Institute of Zoology, CAS who kindly offered critical comments and valuable suggestions on the manuscript. We also thank two anonymous referees and the Associate Editor for helpful comments on the manuscript. This work was funded by grants from National Natural Science Foundation of China (No. 30870282, 31071911), and Innovative Program of The Chinese Academy of Sciences (No. KSCX2-YW-Z-0930).
- Silvestri F: Descrizione di un novo genere d'insetti apterigoti rappresentante di un novo ordine. Boll Lab Zool Portici. 1907, 1: 296-311.Google Scholar
- Yin WY: Fauna Sinica Arthropoda: Protura. 1999, Beijing: Science Press, in Chinese with English summaryGoogle Scholar
- Yin WY: A new idea on phylogeny of Protura with approach to its origin and systematic position. Sci Sin Ser B. 1984, 27: 149-160.Google Scholar
- Yin WY, Xue LZ: Comparative Spermatology of Protura and Its Significance on Proturan Systematics. Sci China Ser B. 1993, 36 (5): 575-586.Google Scholar
- Dallai R, Mercati D, Bu Y, Yin YW, Callaini G, Riparbelli MG: The spermatogenesis and sperm structure of Acerentomon microrhinus (Protura, Hexapoda) with considerations on the phylogenetic position of the taxon. Zoomorphology. 2010, 129 (1): 61-80. 10.1007/s00435-009-0100-1.View ArticleGoogle Scholar
- Hennig W: Insect phylogeny. 1981, New York: John Wiley & SonsGoogle Scholar
- Dallai R: Are Protura really insects?. The Early Evolution of Metazoa and the Significance of Problematic Taxa. Edited by: Simonetta AM, Morris SC. 1991, Cambridge: The Cambridge University Press, 263-269.Google Scholar
- Luan YX, Mallatt JM, Xie RD, Yang YM, Yin WY: The phylogenetic positions of three basal-hexapod groups (Protura, Diplura, and Collembola) based on ribosomal RNA gene sequences. Mol Biol Evol. 2005, 22 (7): 1579-1592. 10.1093/molbev/msi148.View ArticlePubMedGoogle Scholar
- Machida R: Evidence from embryology for reconstructing the relationships of hexapod basal clades. Arthropod Systematics & Phylogeny. 2006, 64 (1): 95-104.Google Scholar
- Bitsch C, Bitsch J: Phylogenetic relationships of basal hexapods among the mandibulate arthropods: a cladistic analysis based on comparative morphological characters. Zoolog Sci. 2004, 33 (6): 511-550. 10.1111/j.0300-3256.2004.00162.x.View ArticleGoogle Scholar
- Szeptycki A: Catalogue of the world Protura. Acta Zoologica Cracoviensia. 2007, 50B (1):
- Meusemann K, von Reumont BM, Simon S, Roeding F, Strauss S, Kuck P, Ebersberger I, Walzl M, Pass G, Breuers S, et al: A phylogenomic approach to resolve the arthropod tree of life. Mol Biol Evol. 2010, 27 (11): 2451-2464. 10.1093/molbev/msq130.View ArticlePubMedGoogle Scholar
- Regier JC, Shultz JW, Zwick A, Hussey A, Ball B, Wetzer R, Martin JW, Cunningham CW: Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences. Nature. 2010, 463 (7284): 1079-1083. 10.1038/nature08742.View ArticlePubMedGoogle Scholar
- Rota-Stabelli O, Kayal E, Gleeson D, Daub J, Boore JL, Telford MJ, Pisani D, Blaxter M, Lavrov DV: Ecdysozoan mitogenomics: evidence for a common origin of the legged invertebrates, the Panarthropoda. Genome Biol Evol. 2010, 2: 425-440. 10.1093/gbe/evq030.View ArticlePubMedPubMed CentralGoogle Scholar
- Friedrich M, Tautz D: Ribosomal DNA phylogeny of the major extant arthropod classes and the evolution of myriapods. Nature. 1995, 376 (6536): 165-167. 10.1038/376165a0.View ArticlePubMedGoogle Scholar
- Boore JL, Lavrov DV, Brown WM: Gene translocation links insects and crustaceans. Nature. 1998, 392 (6677): 667-668. 10.1038/33577.View ArticlePubMedGoogle Scholar
- Carapelli A, Lio P, Nardi F, van der Wath E, Frati F: Phylogenetic analysis of mitochondrial protein coding genes confirms the reciprocal paraphyly of Hexapoda and Crustacea. BMC Evol Biol. 2007, 7 (Suppl 2): S8-10.1186/1471-2148-7-S2-S8.View ArticlePubMedPubMed CentralGoogle Scholar
- Mallatt J, Craig CW, Yoder MJ: Nearly complete rRNA genes assembled from across the metazoan animals: Effects of more taxa, a structure-based alignment, and paired-sites evolutionary models on phylogeny reconstruction. Mol Phylogenet Evol. 2010, 55 (1): 1-17. 10.1016/j.ympev.2009.09.028.View ArticlePubMedGoogle Scholar
- Nardi F, Spinsanti G, Boore JL, Carapelli A, Dallai R, Frati F: Hexapod origins: monophyletic or paraphyletic?. Science. 2003, 299 (5614): 1887-1889. 10.1126/science.1078607.View ArticlePubMedGoogle Scholar
- Manton SM, Harding MJP: The evolution of arthropodan locomotory mechanisms. Part 10. Locomotory habits, morphology and evolution of the hexapod classes. Zool J Linn Soc-Lond. 1972, 51: 203-400. 10.1111/j.1096-3642.1972.tb02550.x.View ArticleGoogle Scholar
- Cook CE, Yue QY, Akam M: Mitochondrial genomes suggest that hexapods and crustaceans are mutually paraphyletic. P Roy Soc B-Biol Sci. 2005, 272 (1569): 1295-1304. 10.1098/rspb.2004.3042.View ArticleGoogle Scholar
- Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27 (8): 1767-1780. 10.1093/nar/27.8.1767.View ArticlePubMedPubMed CentralGoogle Scholar
- Negrisolo E, Minelli A, Valle G: The mitochondrial genome of the house centipede scutigera and the monophyly versus paraphyly of myriapods. Mol Biol Evol. 2004, 21 (4): 770-780. 10.1093/molbev/msh078.View ArticlePubMedGoogle Scholar
- Boore JL, Collins TM, Stanton D, Daehler LL, Brown WM: Deducing the pattern of arthropod phylogeny from mitochondrial DNA rearrangements. Nature. 1995, 376 (6536): 163-165. 10.1038/376163a0.View ArticlePubMedGoogle Scholar
- Nardi F, Carapelli A, Fanciulli PP, Dallai R, Frati F: The complete mitochondrial DNA sequence of the basal hexapod Tetrodontophora bielanensis: evidence for heteroplasmy and tRNA translocations. Mol Biol Evol. 2001, 18 (7): 1293-1304.View ArticlePubMedGoogle Scholar
- Podsiadlowski L, Carapelli A, Nardi F, Dallai R, Koch M, Boore JL, Frati F: The mitochondrial genomes of Campodea fragilis and Campodea lubbocki (Hexapoda: Diplura): High genetic divergence in a morphologically uniform taxon. Gene. 2006, 381: 49-61.View ArticlePubMedGoogle Scholar
- Cameron SL, Miller KB, D'Haese CA, Whiting MF, Barker SC: Mitochondrial genome data alone are not enough to unambiguously resolve the relationships of Entognatha, Insecta and Crustacea sensu lato (Arthropoda). Cladistics. 2004, 20 (6): 534-557. 10.1111/j.1096-0031.2004.00040.x.View ArticleGoogle Scholar
- Glenner H, Thomsen PF, Hebsgaard MB, Sorensen MV, Willerslev E: Evolution. The origin of insects. Science. 2006, 314 (5807): 1883-1884. 10.1126/science.1129844.View ArticlePubMedGoogle Scholar
- Budd GE, Telford MJ: The origin and evolution of arthropods. Nature. 2009, 457 (7231): 812-817. 10.1038/nature07890.View ArticlePubMedGoogle Scholar
- Simon C, Buckley TR, Frati F, Stewart JB, Beckenbach AT: Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA. Annu Rev Ecol Evol S. 2006, 37: 545-579. 10.1146/annurev.ecolsys.37.091305.110018.View ArticleGoogle Scholar
- Ojala D, Montoya J, Attardi G: tRNA punctuation model of RNA processing in human mitochondria. Nature. 1981, 290 (5806): 470-474. 10.1038/290470a0.View ArticlePubMedGoogle Scholar
- Wei SJ, Shi M, Chen XX, Sharkey MJ, van Achterberg C, Ye GY, He JH: New Views on Strand Asymmetry in Insect Mitochondrial Genomes. Plos One. 2010, 5 (9):
- Hassanin A, Leger N, Deutsch J: Evidence for multiple reversals of asymmetric mutational constraints during the evolution of the mitochondrial genome of metazoa, and consequences for phylogenetic inferences. Syst Biol. 2005, 54 (2): 277-298. 10.1080/10635150590947843.View ArticlePubMedGoogle Scholar
- Zhang DX, Hewitt GM: Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol. 1997, 25 (2): 99-120. 10.1016/S0305-1978(96)00042-7.View ArticleGoogle Scholar
- Wolstenholme DR, Macfarlane JL, Okimoto R, Clary DO, Wahleithner JA: Bizarre tRNAs inferred from DNA sequences of mitochondrial genomes of nematode worms. Proc Natl Acad Sci USA. 1987, 84 (5): 1324-1328. 10.1073/pnas.84.5.1324.View ArticlePubMedPubMed CentralGoogle Scholar
- Garey JR, Wolstenholme DR: Platyhelminth mitochondrial DNA: evidence for early evolutionary origin of a tRNA(serAGN) that contains a dihydrouridine arm replacement loop, and of serine-specifying AGA and AGG codons. J Mol Evol. 1989, 28 (5): 374-387. 10.1007/BF02603072.View ArticlePubMedGoogle Scholar
- Lavrov DV, Brown WM, Boore JL: A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus. Proc Natl Acad Sci USA. 2000, 97 (25): 13738-13742. 10.1073/pnas.250402997.View ArticlePubMedPubMed CentralGoogle Scholar
- Segovia R, Pett W, Trewick S, Lavrov DV: Extensive and evolutionarily persistent mitochondrial tRNA editing in velvet worms (phylum Onychophora). Mol Biol Evol. 2011Google Scholar
- Masta SE, Boore JL: The complete mitochondrial genome sequence of the spider Habronattus oregonensis reveals rearranged and extremely truncated tRNAs. Mol Biol Evol. 2004, 21 (5): 893-902. 10.1093/molbev/msh096.View ArticlePubMedGoogle Scholar
- Klimov PB, OConnor BM: Improved tRNA prediction in the American house dust mite reveals widespread occurrence of extremely short minimal tRNAs in acariform mites. Bmc Genomics. 2009, 10: 598-10.1186/1471-2164-10-598.View ArticlePubMedPubMed CentralGoogle Scholar
- Beckenbach AT, Joy JB: Evolution of the mitochondrial genomes of gall midges (Diptera: Cecidomyiidae): rearrangement and severe truncation of tRNA genes. Genome Biol Evol. 2009, 1: 278-287.View ArticlePubMedPubMed CentralGoogle Scholar
- Moritz C, Dowling TE, Brown WM: Evolution of Animal Mitochondrial-DNA - Relevance for Population Biology and Systematics. Annu Rev Ecol Syst. 1987, 18: 269-292. 10.1146/annurev.es.18.110187.001413.View ArticleGoogle Scholar
- Negrisolo E, Minelli A, Valle G: Extensive gene order rearrangement in the mitochondrial genome of the centipede Scutigera coleoptrata. J Mol Evol. 2004, 58 (4): 413-423. 10.1007/s00239-003-2563-x.View ArticlePubMedGoogle Scholar
- Macey JR, Larson A, Ananjeva NB, Fang ZL, Papenfuss TJ: Two novel gene orders and the role of light-strand replication in rearrangement of the vertebrate mitochondrial genome. Mol Biol Evol. 1997, 14 (1): 91-104.View ArticlePubMedGoogle Scholar
- Lavrov DV, Boore JL, Brown WM: Complete mtDNA sequences of two millipedes suggest a new model for mitochondrial gene rearrangements: Duplication and nonrandom loss. Mol Biol Evol. 2002, 19 (2): 163-169.View ArticlePubMedGoogle Scholar
- San Mauro D, Gower DJ, Zardoya R, Wilkinson M: A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome. Mol Biol Evol. 2006, 23 (1): 227-234.View ArticlePubMedGoogle Scholar
- Lunt DH, Hyman BC: Animal mitochondrial DNA recombination. Nature. 1997, 387 (6630): 247-247. 10.1038/387247a0.View ArticlePubMedGoogle Scholar
- Cameron SL, Dowton M, Castro LR, Ruberu K, Whiting MF, Austin AD, Diement K, Stevens J: Mitochondrial genome organization and phylogeny of two vespid wasps. Genome. 2008, 51 (10): 800-808. 10.1139/G08-066.View ArticlePubMedGoogle Scholar
- Rawlings TA, Collins TM, Bieler R: Changing identities: tRNA duplication and remolding within animal mitochondrial genomes. Proc Natl Acad Sci USA. 2003, 100 (26): 15700-15705. 10.1073/pnas.2535036100.View ArticlePubMedPubMed CentralGoogle Scholar
- Covacin C, Shao R, Cameron S, Barker SC: Extraordinary number of gene rearrangements in the mitochondrial genomes of lice (Phthiraptera: Insecta). Insect Mol Biol. 2006, 15 (1): 63-68. 10.1111/j.1365-2583.2005.00608.x.View ArticlePubMedGoogle Scholar
- Kilpert F, Podsiadlowski L: The Australian fresh water isopod (Phreatoicidea: Isopoda) allows insights into the early mitogenomic evolution of isopods. Comp Biochem Phys D. 2010, 5 (1): 36-44.Google Scholar
- Shao RF, Barker SC: The highly rearranged mitochondrial genome of the plague thrips, Thrips imaginis (Insecta: thysanoptera): Convergence of two novel gene boundaries and an extraordinary arrangement of rRNA genes. Mol Biol Evol. 2003, 20 (3): 362-370. 10.1093/molbev/msg045.View ArticlePubMedGoogle Scholar
- Shao R, Dowton M, Murrell A, Barker SC: Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects. Mol Biol Evol. 2003, 20 (10): 1612-1619. 10.1093/molbev/msg176.View ArticlePubMedGoogle Scholar
- Shao R, Campbell NJ, Schmidt ER, Barker SC: Increased rate of gene rearrangement in the mitochondrial genomes of three orders of hemipteroid insects. Mol Biol Evol. 2001, 18 (9): 1828-1832.View ArticlePubMedGoogle Scholar
- Dowton M, Castro LR, Campbell SL, Bargon SD, Austin AD: Frequent mitochondrial gene rearrangements at the hymenopteran nad3-nad5 junction. J Mol Evol. 2003, 56 (5): 517-526. 10.1007/s00239-002-2420-3.View ArticlePubMedGoogle Scholar
- Sheffield NC, Song H, Cameron SL, Whiting MF: Nonstationary evolution and compositional heterogeneity in beetle mitochondrial phylogenomics. Syst Biol. 2009, 58 (4): 381-394. 10.1093/sysbio/syp037.View ArticlePubMedGoogle Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994, 3 (5): 294-299.PubMedGoogle Scholar
- Boore JL, Brown WM: Mitochondrial genomes of Galathealinum, Helobdella, and Platynereis: sequence and gene arrangement comparisons indicate that Pogonophora is not a phylum and Annelida and Arthropoda are not sister taxa. Mol Biol Evol. 2000, 17 (1): 87-106.View ArticlePubMedGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8 (3): 175-185.View ArticlePubMedGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.View ArticlePubMedGoogle Scholar
- Burland TG: DNASTAR's Lasergene sequence analysis software. Methods Mol Biol. 2000, 132: 71-91.PubMedGoogle Scholar
- Wyman SK, Jansen RK, Boore JL: Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004, 20 (17): 3252-3255. 10.1093/bioinformatics/bth352.View ArticlePubMedGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25 (5): 955-964. 10.1093/nar/25.5.955.View ArticlePubMedPubMed CentralGoogle Scholar
- Laslett D, Canback B: ARWEN: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences. Bioinformatics. 2008, 24 (2): 172-175. 10.1093/bioinformatics/btm573.View ArticlePubMedGoogle Scholar
- Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27 (2): 573-580. 10.1093/nar/27.2.573.View ArticlePubMedPubMed CentralGoogle Scholar
- Gruber AR, Lorenz R, Bernhart SH, Neubock R, Hofacker IL: The Vienna RNA websuite. Nucleic Acids Res. 2008, 36 (Web Server): W70-74. 10.1093/nar/gkn188.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee YS, Oh J, Kim YU, Kim N, Yang S, Hwang UW: Mitome: dynamic and interactive database for comparative mitochondrial genomics in metazoan animals. Nucleic Acids Res. 2008, 36: D938-D942.View ArticlePubMedPubMed CentralGoogle Scholar
- Perna NT, Kocher TD: Patterns of Nucleotide Composition at Fourfold Degenerate Sites of Animal Mitochondrial Genomes. J Mol Evol. 1995, 41 (3): 353-358. 10.1007/BF01215182.View ArticlePubMedGoogle Scholar
- Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92 (4): 371-373. 10.1093/jhered/92.4.371.View ArticlePubMedGoogle Scholar
- Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17 (4): 540-552.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol.
- Stamatakis A, Hoover P, Rougemont J: A Rapid Bootstrap Algorithm for the RAxML Web Servers. Syst Biol. 2008, 57 (5): 758-771. 10.1080/10635150802429642.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17 (8): 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.