- Research article
- Open Access
The complete mitochondrial genome of Pseudocellus pearsei (Chelicerata: Ricinulei) and a comparison of mitochondrial gene rearrangements in Arachnida
BMC Genomics volume 8, Article number: 386 (2007)
Mitochondrial genomes are widely utilized for phylogenetic and population genetic analyses among animals. In addition to sequence data the mitochondrial gene order and RNA secondary structure data are used in phylogenetic analyses. Arachnid phylogeny is still highly debated and there is a lack of sufficient sequence data for many taxa. Ricinulei (hooded tickspiders) are a morphologically distinct clade of arachnids with uncertain phylogenetic affinities.
The first complete mitochondrial DNA genome of a member of the Ricinulei, Pseudocellus pearsei (Arachnida: Ricinulei) was sequenced using a PCR-based approach. The mitochondrial genome is a typical circular duplex DNA molecule with a size of 15,099 bp, showing the complete set of genes usually present in bilaterian mitochondrial genomes. Five tRNA genes (trnW, trnY, trnN, trnL(CUN), trnV) show different relative positions compared to other Chelicerata (e.g. Limulus polyphemus, Ixodes spp.). We propose that two events led to this derived gene order: (1) a tandem duplication followed by random deletion and (2) an independent translocation of trnN. Most of the inferred tRNA secondary structures show the common cloverleaf pattern except tRNA-Glu where the TψC-arm is missing. In phylogenetic analyses (maximum likelihood, maximum parsimony, Bayesian inference) using concatenated amino acid and nucleotide sequences of protein-coding genes the basal relationships of arachnid orders remain unresolved.
Phylogenetic analyses (ML, MP, BI) of arachnid mitochondrial genomes fail to resolve interordinal relationships of Arachnida and remain in a preliminary stage because there is still a lack of mitogenomic data from important taxa such as Opiliones and Pseudoscorpiones. Gene order varies considerably within Arachnida – only eight out of 23 species have retained the putative arthropod ground pattern. Some gene order changes are valuable characters in phylogenetic analysis of intraordinal relationships, e.g. in Acari.
Due to their bacterial origin [1, 2] mitochondria have retained a circular DNA double-helix, which in animals is sized between 12–30 kb. This is only a small part of the original bacterial chromosome, the majority was eliminated or transferred to the nucleus . The mitochondrial DNA of Bilateria typically contains 37 genes and one AT-rich non-coding part, which putatively bears regulatory elements for transcription and translation and is therefore referred to as the mitochondrial control region . In general the genes encode 13 protein subunits necessary for oxidative phosphorylation (atp6+8, cob, cox1–3, nad1–6 and nad4L), 22 transfer RNAs and two rRNAs (rrnS and rrnL) . Except for the control region, mtDNA possesses only few non-coding sections between genes, even gene overlaps are common. E.g., in many species the last seven nucleotides of atp8 are also the first seven nucleotides of atp6. A similar overlap is often seen on the boundary between nad4 and nad4L. As a consequence, rearrangements in mitochondrial genomes most often disrupt genes and thus are deleterious – a possible reason for the stability of mitochondrial gene order .
Mitochondrial genomes have been proven useful for phylogenetic analyses . Nucleotide or amino acid sequences as well as rearrangements in mitochondrial gene order are used as phylogenetic markers [6, 7]. Gene rearrangements are considered to be valuable characters, because it is very unlikely that closely related taxa exhibit homoplastic translocations . In addition the secondary structure of encoded tRNAs  and changes in the mitochondrial genetic code  have also been used as characters in phylogenetic analysis.
Ricinulei (hooded tickspiders) are a small order of arachnids, comprising 3 genera with 55 described species [11–14]. They are predatory animals that live in humid caves or leaf litter of tropical regions [15–17]. Species of Ricinoides occur in West Africa whereas species of Cryptocellus and Pseudocellus live in Central and South America . Ricinuleids have body lengths of 3 to 10 mm  and their cuticle is strongly sclerotized and extraordinarily thick . Several peculiarities characterize ricinuleids – a moveable hood (cucullus) in front of the prosoma covering the mouthparts, two jointed chelicerae, chelate pedipalps, elongated second legs, a tarsal copulatory organ on the third pair of legs of adult males, a locking mechanism between pro- and opisthosoma, which can be unlocked during mating and egg-laying, a 6-legged larvae, the lack of distinct eyes, and tracheal respiration [20–25]. According to morphological studies and some combined morphological and molecular analyses, ricinuleids are often considered to be the sister group of Acari [26–29], together forming the Acaromorpha. This clade is characterized by a unique post-embryonic development: a hexapodal larva followed by three octapod nymphal instars . A gnathosoma with medially fused palpal coxae is another unique character of Acaromporpha. However the definition of a "gnathosoma" varies and its presence in Ricinulei is questioned by some authors [e.g. ]. The analysis of Shultz  obtained seven additional "homoplasious synapomorphies" supporting Acaromorpha (= these characters are not exclusively found in Acari and Ricinulei). Van der Hammen [31, 32] however has questioned the monophyly of Acari and placed Ricinulei within the Acari as the sister group to Anactinotrichida (Parasitiformes s. str. + Opilioacariformes) . Dunlop  suggested Ricinulei as sister group to the extinct Trigonotarbida which are collectively the sister of the extant Tetrapulmonata (Araneae, Uropygi, Amblypygi).
Because available sequence data from nuclear or mitochondrial genes is very limited for Ricinulei, recent molecular studies of arthropod phylogeny seldom included ricinuleids [e.g. ]. With a combined data set (93 morphological characters and 18S + 28S rRNA sequences) Wheeler & Hayashi  placed Ricinulei as sister group to Acari. In another combined analysis by Giribet et al. (253 morphological characters, 18S and 28S rRNA, ) Ricinulei appear as the sister group of Tetrapulmonata (Araneae + Amblypygi + Uropygi). When fossil taxa were included, a close relationship between the fossil Trigonotarbida and Ricinulei was recovered and both together were the sister group to Tetrapulmonata. This phylogenetic hypothesis corresponds with that of Dunlop , but Giribet et al. mentioned the instability of this relationship: In their analysis of only the molecular data set Ricinulei is early branching within the arachnid tree as the sister group to all remaining arachnids .
In this study, we present the first complete mitochondrial genome of a member of Ricinulei, the hooded tickspider Pseudocellus pearsei (Chamberlin & Ivie, 1938). The sequence data is used to unveil phylogenetic relationships between Arachnida. Furthermore, the gene order of mitochondrial genomes from all available arachnid species is compared in order to reconstruct the events leading to derived genome arrangements and to evaluate the phylogenetic significance of gene translocations within the Arachnida.
Results and Discussion
Genome organization, gene order and non-coding parts
The generation of overlapping PCR fragments and subsequent sequencing demonstrated that the mitochondrial genome of P. pearsei is a typical circular DNA molecule with a length of 15099 bp [GenBank:EU024483]. All 37 genes usually present in bilaterian mitochondrial genomes have been identified (Fig. 1, Table 1). With the exception of five translocated tRNAs, gene order is similar to that of the horseshoe crab, Limulus polyphemus, which is considered to represent the putative ground pattern of the Arthropoda [38, 39]. The tRNAs (trnW, trnY, trnN, trnL(CUN), trnV) changed their position to a new location between trnM and nad2. In this part of the genome there are five non-coding regions, which are ranging in size from 87 to 250 bp. The three largest non-coding regions are located between rrnS and trnI (250 bp), trnY and trnN (183 bp) and trnV and nad2 (169 bp). The shorter ones are situated between trnI and trnQ (100 bp) and trnM and trnW (87 bp). The longest non-coding region (250 bp) is flanked by rrnSI and trn, and based on its similarity to other arthropods [5, 40, 41] can be identified as the putative mitochondrial control region. Part of this region is capable of folding into a hairpin-like formation (Fig. 2) with a loop consisting of 10 nucleotides and a stem composed of 21 paired nucleotides (five mismatches). Furthermore conserved motifs occur in the flanking sequences around the stem-loop structure: a TATA motif appears in the 5'-flanking sequence whereas the motif GA(A)T is found in the 3'-flanking sequence (Fig. 2). Both motifs are also present in flanking sequences of other arthropods, e.g. in metastriate ticks , crustaceans , and insects and are presumed to have functional significance in transcription and/or replication . The other larger non-coding regions do not bear similar hairpin-like structures.
Besides the above mentioned regions only smaller non-coding regions of 2–23 bp are found in the mitochondrial genome. Gene overlaps occur between cox1 and cox2 (1 bp), nad4 and nad4L (4 bp), atp6 and atp8 (7 bp), and in four cases tRNA-genes are overlapping each other (Fig. 1, Table 1).
Protein-coding genes and nucleotide composition
All of the 13 identified protein-coding genes begin with one of the common start codons for mtDNA ATG, ATA or ATT (Table 1). Out of these 13 protein-coding genes, six show incomplete stop codons (cox2, cob and nad1, 3, 4, 5). In cox2 the stop codon is truncated and the gene terminates with TA, in the other five protein-coding genes only a single thymine serves as partial stop codon. Similar structural events have also been described for the mtDNA of other species where post-transcriptional polyadenylation completes a T or TA partial stop codon into a functional one .
The CG-skew (= (%C-%G)/(%C+%G)) of mitochondrial genes is a good indicator of the strand specific nucleotide frequency bias [46, 47]. In P. pearsei the CG-skew is positive in all (+)strand encoded genes and negative in all (-)strand encoded genes (Table 1). We defined the (+)strand as the strand bearing the majority of coding sequence. The strand specific bias found in P. pearsei occurs in most other arthropods, while a reversal of that bias has been reported for only a few species [43, 46, 48–50]. Table 2 shows the CG-skews of third codon positions of cox1 for all chelicerates and outgroup taxa included in our phylogenetic analyses. We chose cox1 because this gene is found on (+)strand in all species examined. A reversal of CG-skew is seen in the two scorpions, in three Araneae (but not in Heptathela) and in the parasitiform mite Varroa.
Secondary structure of transfer RNAs
The mitochondrial genome of P. pearsei bears all of the 22 tRNAs commonly found in metazoan mtDNA (Fig. 1, Table 1). Except for tRNA-Glu, all tRNAs possess the typical cloverleaf secondary structure, though the TψC stem is shortened in several tRNAs (Fig. 3). The TψC-arm of tRNA-Glu is entirely absent. It is shortened to a single pair of nucleotides in tRNA-Met and tRNA-Phe and it is composed of just two paired bases in tRNA-Gly, tRNA-His, tRNA-Ile, tRNA-Leu (CUN), tRNA-Ser (AGN) and tRNA-Thr.
We performed phylogenetic analyses with two different data sets: concatenated amino acid and concatenated nucleotide sequences (without third codon positions) from all protein-coding genes. Topologies with best likelihood scores from maximum likelihood (ML) analysis are shown (Fig. 4). Topologies of the 50% majority rule consensus trees from ML bootstrapping and Bayesian Inference (BI) differ only slightly from the best topologies of ML analysis: BI with the amino acid dataset resulted in a basal split between Oltacola and the remainder of chelicerates including Limulus, but that node is not supported by bootstrapping or BI; in ML bootstrapping with the nucleotide dataset there was no resolution between Ricinulei, Scorpiones, Araneae and Acari (ML bootstrap <50%). In all performed analyses, good support was found for monophyly of Scorpiones, Opisthothelae (= all Araneae except Mesothelae, here represented by Heptathela), Ixodes, Leptotrombidium, Ornithodorinae (Carios + Ornithodoros), Ixodidae (Ixodes + Amblyomma + Haemaphysalis + Rhipicephalus), Metastriata (Amblyomma + Haemaphysalis + Rhipicephalus), and Dermanyssina (Metaseiulus + Varroa). Monophyly of Acari-Parasitiformes was also well supported by most analyses, except for the maximum parsimony (MP) bootstrapping (56%). Monophyly of Acari was recovered in the best topology of the nucleotide ML analysis and moderately supported by ML bootstrapping (80%), but not by BI or MP bootstrapping. Acari and Ricinulei form sister groups in the best ML topology in the analysis of the nucleotide alignment, but this clade found support only by BI (0.98), not by ML or MP bootstrapping. In contrast analyses of the amino acid dataset resulted in paraphyletic Acari, as Acariformes (Leptotrombidium) form a clade with Ricinulei, and both together form the sister group to Araneae. Both clades are well supported only by BI (Acariformes + Ricinulei: 1.0; Acariformes + Ricinulei + Araneae: 0.99).
Arachnid monophyly is only supported by ML analysis with the amino acid dataset (bootstrap: 74%). In all other analyses the three basal chelicerate branches – Limulus (Xiphosura), Oltacola (Solifugae), and the remainder of arachnids – are not well resolved. Monophyly of Arachnida excl. Solifugae (in this case Scorpiones + Araneae + Acari + Ricinulei) is well supported in ML and BI analyses of the amino acid alignment and weakly supported by ML analysis of the nucleotide alignment.
The classical, morphology based view of the phylogenetic position of Ricinulei is as the sister group to Acari [26–29]. Because our analysis does not resolve the basal relationships of arachnid orders well, it is not in conflict with actual morphology based analyses of arachnid interrelationships [27, 36]. A sister group relation between Acari and Ricinulei is neither well supported, nor refuted by our results.
The phylogenetic analysis of mitochondrial genome data in chelicerate phylogenetics suffers from several problems. First, the currently incomplete taxon sampling (complete lack of data for Opiliones, Pseudoscorpiones, Palpigradi, Uropygi, and Amblypygi) highlights the preliminary nature of this analysis. While the problem of taxon sampling will be overcome in the near future, other problems lie in the nature of the sequence data itself. Several taxa show a reversed strand bias of nucleotide composition (Table 2), which complicates phylogenetic analyses [46, 49]. This may strongly affect the phylogenetic position of Scorpiones in analyses with mitochondrial datsets . Another problem is the heterogeneous substitution rate among arachnids. Our analysis demonstrates great variability in branch lengths, with very long branches for some Acari species (Leptotrombidium, Varroa, Metaseiulus), and very short branches, e.g. from Limulus and Oltacola.
Mitochondrial gene order variation in Arachnida
Apart from five translocated tRNA-genes, the mitochondrial gene order of P. pearsei does not differ from the putative arthropod ground pattern (Figs. 5, 6). A minimum of two events is required to lead to the derived gene order of P. pearsei: a tandem duplication/random deletion event and a single tRNA gene transposition. Tandem duplication and random loss of genes is widely accepted as a mode of genome shuffling in mitochondria . In the present case a segment ranging from trnL to trnY (ground pattern) can be inferred to be involved in such a duplication event. Subsequently the first copy may have lost trnL, trnV, nad1, and trnC, the second copy may have lost rrnL, rrnS, trnI, trnQ, trnM, trnW, and trnY, as well as the control region. Besides the control region (between rrnS and trnI) there are four larger non-coding regions (87–183 bp) present between trnI and nad2 (Fig. 1, Table 1). This is a further hint towards a tandem duplication random deletion event, although there are no apparent sequence homologies to any of the lost genes. In addition, trnN is located in a novel relative position which is best explained by the transposition of this single gene, the second inferred event.
The comparison of mitochondrial gene order from Arachnida reveals a great variation (Fig. 5). Only eight out of 23 species have retained the arthropod ground pattern as represented by Limulus . We present a parsimonious scenario of gene order changes and mapped these events on a phylogenetic tree (Fig. 6). Three different modes leading to a change in gene order were assumed: (1) tandem duplication random deletion events , (2) inversions, and (3) transposition of single tRNAs. Although the latter mode is currently not well explained by a molecular model, we have assumed that single tRNA transpositions have occurred in all cases in which tandem duplication and random deletion seems to be an implausible explanation, e.g. when only one tRNA gene found a new position distant from its original position, or on the other strand. The tandem duplication random deletion model was proposed when more than one gene was involved and/or non-coding sequence was found between the genes involved. The tree we have used for mapping the events is the best scoring topology from the ML analysis of the amino acid alignment (Fig. 4, lower tree), because this analysis provided bootstrap support for most of the branches. However, regarding gene order changes there is also no character conflict (= homoplastic characters/events) with the other topology obtained from the nucleotide alignment.
None of the gene order characters is capable to resolve interordinal relationships of the taxa studied (Fig. 6). The arthropod ground pattern  is retained in some parasitiform mites (e.g. Ixodes , Carios and Ornithodoros ) (Figs. 5, 6), and in the spider Heptathela (Araneae: Mesothelae ). The mitochondrial genomes of Araneae-Opisthothelae (Ornithoctonus , Nephila, [GenBank:NC_008063], Habronattus ) share seven translocated tRNA genes. A subsequent translocation of trnI is found in Nephila and Habronattus . The camel spider Oltacola gomezi (Solifugae) has a mitochondrial genome arrangement almost unaltered from the arthropod ground pattern [GenBank:EU024482], only trnP was translocated to a new relative position. The scorpions Centruroides  and Mesobuthus  share a translocated and inverted trnQ and both lack trnD.
Several independently derived mitochondrial gene orders are found in Acari. Metastriata (Parasitiformes: Ixodidae) show a derived gene order probably caused by a tandem duplication random deletion event and an additional tRNA translocation (Rhipicephalus , Haemaphysalis , Amblyomma, [GenBank:NC_005963]). The honeybee mite Varroa shows translocation of six tRNA genes . In addition the reversed arrangement of rrnS and the control region in this species is hypothesized to be the result of a tandem duplication and random deletion event.
Metaseiulus (Parasitiformes: Mesostimata) possesses the most derived mitochondrial gene order among arachnids, probably due to multiple rearrangements . This genome is the largest (25 kb) within the Chelicerata due to the presence of duplicate and triplicate regions. The duplicate region contains 18 genes plus a copy of the control region and partial trnL(UUR) sequence while the triplicate region comprises only the control region sequence and partial trnL(UUR) sequence (Fig. 5). The mitochondrial genome of Metaseiulus is also remarkable for the absence of nad3 and nad6, because no other chelicerate is known to have lost any of the protein coding genes. Due to the magnitude of these changes, it is difficult to reliably reconstruct the events leading to the gene order found in Metaseiulus which is unique amongst arachnids.
The same is true for the genus Leptotrombidium (Acariformes) [58, 59]. The three species share a common, derived arrangement with secondary rearrangements in L. pallidum. The secondarily derived gene order of L. pallidum (translocation of trnQ, duplication of rrnL and the presence of four large noncoding regions) is considered to have evolved by a combination of tandem duplication with random deletion plus intergenomic recombination between several genes with subsequent gene conversion .
It is apparent that the two taxa showing the most complicated genome rearrangement history (Leptotrombidium and Metaseiulus) exhibit also the longest branches (= highest substitution rates) in phylogenetic analysis (Fig. 4, 6).
Reversals of strand bias in nucleotide frequency
A reversal of nucleotide strand bias was detected in six species (Table 2), most probably due to three independent events: (1) in the Varroa mite , (2) in Scorpiones and (3) in Araneae-Opisthothelae (Ornithoctonus , Nephila, [GenBank:NC_008063], Habronattus ). These changes are also mapped on the tree in Fig. 6 (characters 2, 4, and 20). The comparatively low substitution rate in the mitochondrial genome of Heptathela (Araneae: Mesothelae; Fig. 6) argues against a re-reversal of nucleotide bias within Araneae and suggests that the reversal is in fact a synapomorphy of the Opisthothelae (Fig. 6).
The first complete mitochondrial sequence of a hooded tickspider (P. pearsei, Ricinulei, Arachnida) reveals a typical circular duplex DNA molecule with a compact gene organisation as found in other bilaterians. In comparison to the putative arthropod ground pattern we observed a derived gene order with five tRNA genes found in different relative positions compared to the gene order of Limulus. Probably two events led to the derived gene order of Pseudocellus: (1) a duplication and random deletion event may be responsible for the translocation of four of these tRNA genes and (2) an additional translocation of the trnN. The putative mitochondrial control region is situated in the same position as in other arthropod mitochondrial genomes. Part of the putative control region can be folded into a characteristic stem-loop structure with conserved flanking sequences as found in other arthropods. All tRNAs, except tRNA-Glu, can be folded into a typical cloverleaf secondary structure.
Aligments of nucleotide and amino acid sequences from mitochondrial protein-coding genes were used in a phylogenetic analysis of arachnid relationships. In the best scoring topologies from ML analyses, Ricinulei appear as sister taxon to Acari (nucleotide alignment) or Acariformes (amino acid alignment), but in both cases without sufficient support from bootstrapping or Bayesian inference. Within the Acari, monophyly of Parasitiformes is well supported, while Acariformes are represented by only a single genus (Leptotrombidium). Monophyly of Acari was only recovered in the analysis of the nucleotide alignment and found strong support by BI and only weak support by ML/MP bootstrapping.
Because of highly divergent substitution rates amongst arachnid species, phylogenetic analyses may be generally biased due to long-branch attraction. Another complicating factor is the reversal of nucleotide bias, which has occurred independently in the arachnid clades Scorpiones, Acari (Varroa) and Araneae (Opisthothelae). Due to the lack of support for interordinal relations within Arachnida, the phylogenetic analysis of mitochondrial genome data reveals no strong conflict with recent morphological analyses.
Investigation of mitochondrial gene rearrangements across the range of taxa for which complete mitochondrial genomes are available on GenBank, reveals high variability in gene order within arachnids. Among arachnids, 15 of the 23 species investigated have derived features in gene order. When mapped onto the tree, the hypotheses of events leading to gene order changes are not in conflict with our phylogenetic analyses of sequence data. None of the hypothesized events is useful as a phylogenetic character to resolve interordinal relationships within arachnids, but some of these characters are promising in analyses of intraordinal relationships. Especially in Acari the comparison of mitochondrial genomes from a larger taxon sampling would be promising for the detection of gene order changes, which will be valuable in a phylogenetic analysis. Altogether, further mitogenomic data of a broader taxon sampling are necessary, especially from Opiliones, Pseudoscorpiones, Palpigradi, Uropygi and Amblypygi.
Samples and DNA extraction
P. pearsei specimens were collected by Gerd Alberti in the Gruta Sabac-Ha, near Merida, province of Yucatán, México (20°10'18"N, 89°16'03"W). DNA extraction was done with two legs obtained from a freshly killed animal. We used "DNeasy blood and tissue kit" (Qiagen, Hilden, Germany) according to the manufacturers protocol, except for reduction of the elution volume to 50 μl.
PCR and sequencing
Initial PCRs were done with a published primer set designed for crustacean mitochondrial genomes . All PCRs were performed on an Eppendorf Mastercycler and Mastercycler gradient. The Eppendorf 5-prime-Taq kit (Eppendorf, Germany) was used in 50 μl volumes (5 μl buffer; 1 μl dNTP mix, 10 μM; 0.25 μl Taq polymerase; 1 μl template DNA, 40,75 μl water, 1 μl primer mix, 10 μM each). PCR conditions were: initial denaturation (94°C, 1 min), 40 cycles of denaturation (94°C, 30 sec), annealing (50°C, 30 sec), and elongation (68°C, 1 min), followed by a final elongation step (68°C, 1 min). Successful PCR amplification was obtained with six primer pairs: S8(cox1, about 650 bp), S10(cox1, about 1050 bp), S15(cox2-atp6, about 830 bp), S20(cox3-trnF, about 950 bp), S37(nd4l-cob, about 1600 bp), and S46(rrnL-rrnS, about 1000 bp), see  for primer sequences. PCR products were visualized on 1% agarose gels. Before sequencing PCR products were purified using the Bluematrix DNA purification kit (EURx, Gdansk, Poland). Except for S8 and S10 there was no overlap between these initial sequences. Species specific primer pairs were designed to bridge the five gaps between these initial sequences (see Table 3 for primer information). Long PCRs were performed with Takara LA Taq kit (Takara) in 50 μl volumes (5 μl Buffer, Mg+; 8 μl dNTP mix, 2.5 μM; 1 μl DNA; 1 μl primer mix, 10 μM each; 0.5 μl Takara LA Taq; 32 μl water). PCR conditions were: initial denaturation (94°C, 1 min), 40 cycles of denaturation (94°C, 30 sec), annealing (primer specific temperature see Table 3, 1 min), and elongation (68°C, 3 min), followed by a final elongation step (68°C, 2 min). Long PCR products were sequenced by primer walking. All sequencing was performed on a CEQ 8000 capillary sequencer (Beckmann-Coulter) using CEQ DCTS kits (Beckmann-Coulter) with 10 μl reaction volumes (4 μl DCTS master mix, 1 μl primer, 10 μM, 1–5 μl DNA, 0–4 μl water). Sequencing reaction was done with 30 cycles of denaturation (94°C, 30 sec), annealing (primer specific temperature, 30 sec) and elongation (60°C, 2 min).
Sequence assemblage and annotation
Primary sequence analysis was performed with the CEQ software (quality check). Sequence assemblage was done with Bioedit 7.0.1 . Protein-coding and ribosomal RNA genes and gene boundaries were identified by BLAST search and in comparison with alignments from other chelicerate species. Genomic position and secondary structure from 20 out of 22 transfer RNAs were identified by tRNA-scan SE , the remaining two by eye inspection of the regions under suspect. Sequence data was deposited at NCBI database [GenBank:EU024483].
Alignments from all protein-coding genes were used in phylogenetic analysis. Alignments of amino acid sequences and nucleotide sequences were performed with ClustalW , as implemented in Bioedit 7.0.1 . Nucleotide sequences were translated to amino acid sequence prior to the multiple alignment and translated back afterwards (making use of the "toggle" function of Bioedit). Sequence data was obtained from 24 chelicerate species and three outgroup taxa (for a list of taxa with GenBank accession numbers see Table 2). Ambiguously aligned regions were omitted using Gblocks ver. 0.91b , using default settings, except for changing "allowed gap positions" to "with half". The final amino acid alignment consisted of 2786 amino acids. For the nucleotide dataset the "codons" option was used, so that only complete codons (not single nucleotides) were omitted from the alignment by Gblocks. A saturation analysis  was performed with DAMBE 4.2.13  using subsets of the nucleotide dataset representing first, second and third codon positions. Due to significant saturation, third codon positions were omitted from the final nucleotide alignment (6490 nucleotides).
Maximum parsimony (MP) analysis was performed with PAUP* ver. 4.0b10 , 1,000 bootstrap replicates were performed, each with 10 replicates with random addition of taxa. Maximum likelihood (ML) analysis was done with Treefinder . Model selection was done with ProtTest ver. 1.3 . According to the Akaike information criterion, the mtArt+G+I model was optimal for ML analysis. The tree with the best likelihood value was chosen out of 10 independent analyses with random taxon order. Bootstrapping was performed with 100 pseudoreplicates. Bayesian inference (BI) with MrBayes ver. 3.1.2  was done with the model mtRev+G+I, as the mtArt model is not implemented in the current version of MrBayes. Eight chains ran for 1,000,000 generations, while tree sampling was done every 1,000 generations; we examined likelihood values of the sampled trees and decided to omit the first 100 trees as burn-in, while the remaining 900 were used to calculate Bayesian posterior probabilities (BPP).
- atp6 and 8:
genes encoding ATPase subunit 6 and 8
Bayesian posterior probability
genes encoding cytochrome oxidase subunits I-III
- cob :
gene encoding cytochrome b
- nad1-6 and nad4L:
genes encoding NADH dehydrogenase subunits 1–6 and 4L
polymerase chain reaction
- rrnL :
large (16S) rRNA subunit (gene)
- rrnS :
small (12S) rRNA subunit (gene)
tandem duplication and random deletion
- transfer RNA:
tRNA-Xxx (where Xxx is replaced by three letter amino acid code of the corresponding amino acid)
- tRNA gene:
trnX (where X is replaced by one letter amino acid code of the corresponding amino acid)
T-loop and T-stem of a tRNA secondary structure.
Altmann R: Die Elementarorganismen und ihre Beziehungen zu den Zellen. 1890, Leipzig, Veit & comp
Margulis L: Origin of eukaryotic cellsevidence and research implications for a theory of the origin and evolution of microbial, plant, and animal cells on the Precambrian earth. 1970, New Haven, Yale University Press,
Adams KL, Palmer JD: Evolution of mitochondrial gene content: gene loss and transfer to the nucleus. Mol Phylogenet Evol. 2003, 29: 380-395. 10.1016/S1055-7903(03)00194-5.
Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27: 1767-1780. 10.1093/nar/27.8.1767.
Wolstenholme DR: Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992, 141: 173-216.
Boore JL, Lavrov DV, Brown WM: Gene translocation links insects and crustaceans. Nature. 1998, 392: 667-668. 10.1038/33577.
Boore JL, Collins TM, Stanton D, Daehler LL, Brown WM: Deducing the pattern of arthropod phylogeny from mitochondrial DNA rearrangements. Nature London. 1995, 376: 163-165. 10.1038/376163a0.
Dowton M, Castro LR, Austin AD: Mitochondrial gene rearrangements as phylogenetic characters in the invertebrates: The examination of genome 'morphology'. Invertebrate Systematics. 2002, 16: 345-356. 10.1071/IS02003.
Macey JR, Schulte JA, Larson A: Evolution and phylogenetic information content of mitochondrial genomic structural features illustrated with acrodont lizards. Syst Biol. 2000, 49: 257-277. 10.1080/10635159950173843.
Castresana J, Feldmaier-Fuchs G, Paabo S: Codon reassignment and amino acid composition in hemichordate mitochondria. Proc Natl Acad Sci U S A. 1998, 95: 3703-3707. 10.1073/pnas.95.7.3703.
Bonaldo AB, Pinto-da-Rocha R: On a new species of Cryptocellus from the Brasilian Amazonas (Arachnida, Ricinulei). Revista Ibérica de Aracnología. 2003, 7: 103-108.
Cokendolpher J, Enríquez T: A new species and records of Pseudocellus (Arachnida: Ricinulei: Ricinoididae) from caves in Yucatán, Mexico and Belize. Texas Memorial Museum, Speleological Monographs. 2004, 6: 95-99.
Harvey MS: Catalogue of the Smaller Arachnid Orders of the World: Amblypygi, Uropygi, Schizomida, Palpigradi, Ricinulei and Solifugae. 2003, CSIRO Publishing, 1-400.
Pinto-da-Rocha R, Bonaldo AB: A new species of Cryptocellus (Arachnida, Ricinulei) from Oriental Amazonia. Zootaxa. 2007, 1386: 47-51.
Adis JU, Platnick NI, de Morais JW, Rodrigues JMG: On the abundance and ecology of Ricinulei (Arachnida) from Central Amazonia, Brazil. Journal of the New York Entomological Society. 1989, 97: 133-140.
Cooke JAL: Observations on the biology of Ricinulei (Arachnida) with description of two new species of Cryptocellus. Journal of Zoology, London. 1967, 151: 31-42.
Mitchell RW: Population size and dispersion and species associations of a Mexican cavernicole ricinuleid (Arachnida). Ciencia, Mexico. 1970, 27: 63-74.
Moritz M: Überklasse Chelicerata, Fühlerlose. Lehrbuch der speziellen Zoologie/ Band I: Wirbellose Tiere/ 4 Teil: Arthropoda (ohne Insecta). Edited by: Gruner HE. 1993, Stuttgart, Jena, New York, Gustav Fischer Verlag, 73-442.
Kennaugh JH: An examination of the cuticle of three species of Ricinulei (Arachnida). Journal of Zoology, London. 1968, 156: 393-404.
Ax P: Multicellular animals. A new approach to the phylogenetic order in nature. 1996, Berlin, Heidelberg, New York, Springer, 1-225.
Hansen HJ, Sørensen W: The order Ricinulei. On two Orders of Arachnida. Edited by: Hansen HJ and Sorensen W. 1904, Cambridge, Cambridge University Press, 114-157.
Legg G: The external morphology of a new species of ricinuleid (Arachnida) from Sierra Leone. Zoological Journal of the Linnean Society. 1976, 59: 1-58.
Legg G: Sperm transfer and mating in Ricinoides hanseni (Ricinulei: Arachnida). Journal of Zoology, London. 1977, 182: 51-61.
Selden PA, Dunlop JA: Fossil taxa and relationships of chelicerates. Arthropod Fossils and Phylogeny. Edited by: Edgecombe G. 1998, Columbia University Press, 303-331.
Pittard K, Mitchell RW: Comparative morphology of the life stages of Cryptocellus pelaezi (Arachnida, Ricinulei). Grad Stud Texas Tech Univ. 1972, 1: 3-77.
Weygoldt P, Paulus HF: Studies on the Morphology, Taxonomy and Phylogeny of the Chelicerata. 2. Cladograms and Development of Chelicerata. Zeitschrift fur Zoologische Systematik und Evolutionsforschung. 1979, 17: 177-200.
Shultz JW: A phylogenetic analysis of the arachnid orders based on morphological characters. Zoological Journal of the Linnean Society. 2007, 150: 221-265. 10.1111/j.1096-3642.2007.00284.x.
Lindquist EE: Current theories on the evolution of major groups of Acari and on their relationships with other groups of Arachnida, with consequent implications for their classification. Edited by: D.A.Griffiths CEB. 1984, New York, John Wiley, 28-62. 1
Shultz JW: Evolutionary morphology and phylogeny of Arachnida. Cladistics. 1990, 6: 1-38.
Van der Hammen L: An introduction to comparative arachnology. 1989, The Hague, SPB Academic Publishing
Van der Hammen L: Comparative studies in Chelicerata II. Epimerata (Palpigradi and Actinotrichida). Zool Verh Leiden. 1982, 196: 1-70.
Van der Hammen L: A revised classification of the mites (Arachnidea, Acarida) with diagnoses, a key, and notes on phylogeny. Zool Meded Leiden. 1972, 47: 273-292.
Van der Hammen L: Comparative studies in Chelicerata I. The Cryptognomae (Ricinulei, Architarbi and Anactinotrichida). Zool Verh Leiden. 1979, E. J. Brill, 174: 1-62.
Dunlop JA: Evidence for a sister group relationship between Ricinulei and Trigonotarbida. Bulletin of the British Arachnological Society. 1996, 10: 193-204.
Mallatt J, Giribet G: Further use of nearly complete 28S and 18S rRNA genes to classify Ecdysozoa: 37 more arthropods and a kinorhynch. Mol Phylogenet Evol. 2006, 40: 772-794. 10.1016/j.ympev.2006.04.021.
Wheeler WC, Hayashi CY: The phylogeny of the extant chelicerate orders. Cladistics. 1998, 14: 173-192. 10.1111/j.1096-0031.1998.tb00331.x.
Giribet G, Edgecombe GD, Wheeler WC, Babbitt C: Phylogeny and systematic position of Opiliones: a combined analysis of chelicerate relationships using morphological and molecular data. Cladistics. 2002, 18: 5-70.
Staton JL, Daehler LL, Brown WM: Mitochondrial gene arrangement of the horseshoe crab Limulus polyphemus L.: conservation of major features among arthropod classes. Mol Biol Evol. 1997, 14: 867-874.
Lavrov DV, Boore JL, Brown WM: The complete mitochondrial DNA sequence of the horseshoe crab Limulus polyphemus. Mol Biol Evol. 2000, 17: 813-824.
Zhang DX, Hewitt GM: Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochemical Systematics and Ecology. 1997, 25: 99-120. 10.1016/S0305-1978(96)00042-7.
Qiu Y, Song D, Zhou K, Sun H: The mitochondrial sequences of Heptathela hangzhouensis and Ornithoctonus huwena reveal unique gene arrangements and atypical tRNAs. J Mol Evol. 2005, 60: 57-71. 10.1007/s00239-004-0010-2.
Black WC, Roehrdanz RL: Mitochondrial gene order is not conserved in arthropods: prostriate and metastriate tick mitochondrial genomes. Mol Biol Evol. 1998, 15: 1772-1785.
Kilpert F, Podsiadlowski L: The complete mitochondrial genome of the common sea slater, Ligia oceanica (Crustacea, Isopoda) bears a novel gene order and unusual control region features. BMC Genomics. 2006, 7: 241-10.1186/1471-2164-7-241.
Zhang DX, Szymura JM, Hewitt GM: Evolution and structural conservation of the control region of insect mitochondrial DNA. J Mol Evol. 1995, 40: 382-391. 10.1007/BF00164024.
Ojala D, Montoya J, Attardi G: tRNA punctuation model of RNA processing in human mitochondria. Nature. 1981, 290: 470-474. 10.1038/290470a0.
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: 277-298. 10.1080/10635150590947843.
Perna NT, Kocher TD: Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995, 41: 353-358. 10.1007/BF01215182.
Podsiadlowski L, Bartolomaeus T: Major rearrangements characterize the mitochondrial genome of the isopod Idotea baltica (Crustacea: Peracarida). Mol Phylogenet Evol. 2006, 40: 893-899. 10.1016/j.ympev.2006.04.008.
Hassanin A: Phylogeny of Arthropoda inferred from mitochondrial sequences: strategies for limiting the misleading effects of multiple changes in pattern and rates of substitution. Mol Phylogenet Evol. 2006, 38: 100-116. 10.1016/j.ympev.2005.09.012.
Navajas M, Le Conte Y, Solignac M, Cros-Arteil S, Cornuet JM: The complete sequence of the mitochondrial genome of the honeybee ectoparasite mite Varroa destructor (Acari: Mesostigmata). Mol Biol Evol. 2002, 19: 2313-2317.
Jones M, Gantenbein B, Fet V, Blaxter M: The effect of model choice on phylogenetic inference using mitochondrial sequence data: Lessons from the scorpions. Molecular Phylogenetics and Evolution. 2007, 43: 583-595. 10.1016/j.ympev.2006.11.017.
Moritz C, Dowling TE, Brown WM: Evolution of animal mitochondrial-DNA - relevance for population biology and systematics. Annual Review of Ecology and Systematics. 1987, 18: 269-292. 10.1146/annurev.es.18.110187.001413.
Shao R, Aoki Y, Mitani H, Tabuchi N, Barker SC, Fukunaga M: The mitochondrial genomes of soft ticks have an arrangement of genes that has remained unchanged for over 400 million years. Insect Mol Biol. 2004, 13: 219-224. 10.1111/j.0962-1075.2004.00447.x.
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: 893-902. 10.1093/molbev/msh096.
Davila S, Pinero D, Bustos P, Cevallos MA, Davila G: The mitochondrial genome sequence of the scorpion Centruroides limpidus (Karsch 1879) (Chelicerata; Arachnida). Gene. 2005, 360: 92-102. 10.1016/j.gene.2005.06.008.
Choi EH, Park SJ, Jang KH, Hwang W: Complete mitochondrial genome of a Chinese scorpion Mesobuthus martensii (Chelicerata, Scorpiones, Buthidae). DNA Seq. 2007, 18: 459-471.
Jeyaprakash A, Hoy MA: The mitochondrial genome of the predatory mite Metaseiulus occidentalis (Arthropoda : Chelicerata : Acari : Phytoseiidae) is unexpectedly large and contains several novel features. Gene. 2007, 391: 264-274. 10.1016/j.gene.2007.01.012.
Shao R, Mitani H, Barker SC, Takahashi M, Fukunaga M: Novel mitochondrial gene content and gene arrangement indicate illegitimate inter-mtDNA recombination in the chigger mite, Leptotrombidium pallidum. J Mol Evol. 2005, 60: 764-773. 10.1007/s00239-004-0226-1.
Shao R, Barker SC, Mitani H, Takahashi M, Fukunaga M: Molecular mechanisms for the variation of mitochondrial gene content and gene arrangement among chigger mites of the genus Leptotrombidium (Acari: Acariformes). J Mol Evol. 2006, 63: 251-261. 10.1007/s00239-005-0196-y.
Yamauchi MM, Miya MU, Machida RJ, Nishida M: PCR-based approach for sequencing mitochondrial genomes of decapod crustaceans, with a practical example from kuruma prawn (Marsupenaeus japonicus). Mar Biotechnol (NY). 2004, 6: 419-429. 10.1007/s10126-003-0036-2.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964. 10.1093/nar/25.5.955.
Thompson JD, Higgins DG, Gibson TJ: Clustal-W - Improving the Sensitivity of Progressive Multiple Sequence Alignment Through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice. Nucleic Acids Research. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17: 540-552.
Xia X, Xie Z, Salemi M, Chen L, Wang Y: An index of substitution saturation and its application. Molecular Phylogenetics and Evolution. 2003, 26: 1-7. 10.1016/S1055-7903(02)00326-3.
Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92: 371-373. 10.1093/jhered/92.4.371.
Swofford DL: Paup - A Computer-Program for Phylogenetic Inference Using Maximum Parsimony. Journal of General Physiology. 1993, 102: A9-A9.
Jobb G: Treefinder. Version of Feb. 2007. Munich, Germany. 2007, [http://www.treefinder.de]
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.
The authors thank Thomas Bartolomaeus (Freie Universität Berlin, Germany) for his kind support through all stages of this work. We are also grateful to Gerd Alberti (Universität Greifswald), Gerhard Scholtz (Humboldt-Universität Berlin), and J.G. Palacios-Vargas (México) for further support. GT wishes to express his thanks to Tanja Müther and Thomas Schierl for their kind hospitality. We also thank three anonymous reviewers for improvement of the manuscript. This study was supported by German Research Foundation (DFG) grants Ba 1520/10-1 (to LP) and DFG Scho 442/8-1 (to AB), both from Priority Program 1174 (Deep Metazoan Phylogeny), and Al 138/13-1 (to GT).
The author(s) declares that there are no competing interests.
GT sampled and determined animals and isolated DNA. KF, GT and LP performed PCRs, KF and LP did all sequencing, sequence analysis and annotation. KF, AB and LP performed phylogenetic analyses. KF and LP wrote the manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Fahrein, K., Talarico, G., Braband, A. et al. The complete mitochondrial genome of Pseudocellus pearsei (Chelicerata: Ricinulei) and a comparison of mitochondrial gene rearrangements in Arachnida. BMC Genomics 8, 386 (2007). https://doi.org/10.1186/1471-2164-8-386
- Mitochondrial Genome
- Bayesian Inference
- Gene Order
- Sister Group
- tRNA Gene