The apparent scarcity of available sequence data has greatly impeded evolutionary studies in Acari (mites and ticks). This subclass encompasses over 48,000 species and forms the largest group within the Arachnida. Although mitochondrial genomes are widely utilised for phylogenetic and population genetic studies, only 20 mitochondrial genomes of Acari have been determined, of which only one belongs to the diverse order of the Sarcoptiformes. In this study, we describe the mitochondrial genome of the European house dust mite Dermatophagoides pteronyssinus, the most important member of this largely neglected group.
The mitochondrial genome of D. pteronyssinus is a circular DNA molecule of 14,203 bp. It contains the complete set of 37 genes (13 protein coding genes, 2 rRNA genes and 22 tRNA genes), usually present in metazoan mitochondrial genomes. The mitochondrial gene order differs considerably from that of other Acari mitochondrial genomes. Compared to the mitochondrial genome of Limulus polyphemus, considered as the ancestral arthropod pattern, only 11 of the 38 gene boundaries are conserved. The majority strand has a 72.6% AT-content but a GC-skew of 0.194. This skew is the reverse of that normally observed for typical animal mitochondrial genomes. A microsatellite was detected in a large non-coding region (286 bp), which probably functions as the control region. Almost all tRNA genes lack a T-arm, provoking the formation of canonical cloverleaf tRNA-structures, and both rRNA genes are considerably reduced in size. Finally, the genomic sequence was used to perform a phylogenetic study. Both maximum likelihood and Bayesian inference analysis clustered D. pteronyssinus with Steganacarus magnus, forming a sistergroup of the Trombidiformes.
Although the mitochondrial genome of D. pteronyssinus shares different features with previously characterised Acari mitochondrial genomes, it is unique in many ways. Gene order is extremely rearranged and represents a new pattern within the Acari. Both tRNAs and rRNAs are truncated, corroborating the theory of the functional co-evolution of these molecules. Furthermore, the strong and reversed GC- and AT-skews suggest the inversion of the control region as an evolutionary event. Finally, phylogenetic analysis using concatenated mt gene sequences succeeded in recovering Acari relationships concordant with traditional views of phylogeny of Acari.
Approximately 48,000 mite and tick species (Arthropoda: Chelicerata: Arachnida: Acari) [1, 2] have been described to this day. Since the number of undescribed species is thought to be twenty-fold higher, this subclass is by far the most species-rich group among the Arachnida . Acari diversified 400 million years ago and currently three major lineages are recognised: Opilioacariformes, Acariformes and Parasitiformes [2, 3]. The Acariformes comprise two major groups, the Trombidiformes and Sarcoptiformes [2, 4]. Two of the most prominent members of the Sarcoptiformes are the European house dust mite Dermatophagoides pteronyssinus (Trouessart, 1897) and the American house dust mite Dermatophagoides farinae (Hughes, 1961), both belonging to the family of the Pyroglyphidae (cohort Astigmata). Pyroglyphid mites are typical inhabitants of animal nests. In the human environment, they are mainly found in upholstery, textile floor covers and beddings, where they primarily feed on the skin scale fraction in house dust . About 40 years ago, house dust mites were first recognized as one of the major sources of allergens in house dust . The allergenic proteins are found in high concentrations in mite faeces, which, after drying and pulverizing, become airborne and can be inhaled. The presence of these allergens in sensitive persons is able to cause diseases like asthma, dermatitis and rhinitis [7, 8]. In countries with a temperate climate, 6 to 35 per cent of the population is sensitive to house dust mite-derived allergens .
Complete mitochondrial (mt) genome sequences are becoming increasingly important for effective evolutionary and population studies. Mt genome sequences are not only more informative than shorter sequences of individual genes, they also provide sets of genome-level characters, such as the relative position of different genes, RNA secondary structures and modes of control of replication and transcription [10–13]. However, the applicability of mt genomes as a marker of highly divergent lineages is still controversial [14, 15] and remains to be elucidated . In addition, unravelling mt genomes can be of economic importance as well, since several chemical classes of pesticides target mt proteins. Well-known acaricides like acequinocyl and fluacrypyrim affect mt electron transport through the inhibition of the mt encoded cytochrome b in complex III . Also, the economically important class of METI (Mitochondrial Electron Transfer Inhibitors)-acaricides target the mt complex I, although their exact molecular target has not yet been elucidated. Recently, resistance to the acaricide bifenazate was shown to be caused by mutations in the mt encoded cytochrome b and to have evolved rapidly through a short stage of mt heteroplasmy .
At present, the mt genomes of 20 species belonging to the Acari are available at NCBI (, status January 10, 2009). Most of the submitted sequences have the typical features of metazoan mt genomes. They are circular, between 13 and 20 kb in length, contain a coding region with 37 genes (22 tRNAs, 2 rRNAs and 13 protein coding genes) and a relatively small non-coding region. The latter is mostly AT-rich and fulfils a role in the initiation of replication and transcription [20, 21]. Compared to this typical configuration, the mt genomes of Steganacarus magnus, Metaseiulus occidentalis and Leptotrombidium pallidum show some abnormal features. S. magnus lacks 16 of the 22 tRNAs normally present in mt genomes . M. occidentalis has a unusually large mt genome (24.9 kb) resulting from a duplication event of a large fragment of the codon region. Despite its large size, genes coding for nad6 and nad3 were not found during the initial annotation process . L. pallidum on the other hand has 38 mt genes due to a duplication of the 16S-rRNA .
In this study, we analyse the complete mt genome of a member of the Sarcoptiformes, the European house dust mite D. pteronyssinus, after obtaining the complete sequence using a long PCR approach.
Results and discussion
The mt genome of D. pteronyssinus was amplified, using long PCR, in three overlapping fragments. The final assembled sequence was 14,203 bp [GenBank: EU884425; Fig. 1], making it the fifth smallest sequenced genome within the Acari. Only the mt genomes of Tetranychus urticae (13,103 bp), Leptotrombidium akamushi (13,698 bp), Leptotrombidium deliense (13,731 bp) and S. magnus (13,818 bp) are smaller (Table 1 [25–27]). As non-specific amplification artefacts and incomplete coverage of genes are well-known drawbacks of a PCR approach , we checked the genome size by restriction digest on rolling circle amplified mtDNA (Fig. 2). This approach confirmed the sequence size, considering that the relative mobility of mtDNA restriction fragments can show slight (5–12%) deviations compared to their sequence length . The mt genome of D. pteronyssinus is the first mt sequence of a mite belonging to the Astigmata and is together with the mt genome of S. magnus the only representative from the order of the Sarcoptiformes. Adding this genome to the database resulted in 21 publicly available Acari mtDNA sequences. Twelve belong to species in the superorder of the Parasitiformes whereas nine – among which D. pteronyssinus – belong to species in the superorder of the Acariformes.
Nucleotides composition of completely sequenced mt genomes of Acari and Limulus polyphemus*.
* values were obtained from the corresponding GenBank flat-file in the NCBI database (status January 10, 2009)
a D = Dermatophagoides, Am = Amblyomma, As = Ascoschoengastia, C = Carios, H = Haemaphysalis, I = Ixodes, Le = Leptotrombidium, M = Metaseiulus, O = Ornithodoros, R = Rhipicephalus, S = Steganacarus, T = Tetranychus, U = Unionicola, V = Varroa, W = Walchia, Li = Limulus
b S = Acari superorder (A = Acariformes, P = Parasitiformes)
c PCG = Protein coding genes
d duplications of the control region were also considered
e Ref = References; unpub = unpublished
f GC- and AT-skew for the strand coding for cox1, calculated following 
g for these species the largest non-coding region(s) was/were assumed to be the control region(s)
hL. pallidum has a duplication of 16S-rRNA; in this table the largest 16S-rRNA gene is considered
i only single copy genes were considered for protein coding gene length calculation of M. occidentalis
jL. polyphemus belongs to the order of the Xiphosura within the class of the Merostomata
All 37 genes present in a standard metazoan mt genome could be identified (Fig. 1). Gene overlap exists between trnD/atp8 (1 bp), trnR/nad3 (17 bp), trnM/trnS2 (12 bp), trnP/trnV (11 bp), trnV/trnK (7 bp), trnW/trnY (1 bp), trnY/nad1 (4 bp), trnI/trnQ (1 bp) and trnL1/trnC (5 bp). No overlap was found between protein coding genes. Small non-coding regions (> 20 bp) are present between trnS2/trnA (25 bp), trnA/trnP (20 bp) and nad1/nad6 (28 bp). A large non-coding region is positioned between trnF and trnS1 (286 bp). Twenty-five genes of the mt genome of D. pteronyssinus are transcribed on the majority strand (J-strand), whereas the others are oriented on the minority strand (N-strand).
The mt genome of the horseshoe crab Limulus polyphemus is considered to represent the ground pattern for arthropod mt genomes [30, 31]. Comparing the D. pteronyssinus genome to this sequence revealed that only 11 of the 38 gene boundaries in L. polyphemus are conserved in D. pteronyssinus (Fig. 3 ). Moreover, by making use of the pattern search function in the Mitome-database (, status January 10, 2009), the mt gene order of D. pteronyssinus appeared to be unique among arthropods. Remarkably, the relative position of trnL2 (between nad1 and 16S-rRNA), which differentiates the Chelicerata, Myriapoda and Onychophora from the Insecta and Crustacea according to Boore [12, 34], is not conserved. However, Boore's hypothesis was based on mt genome data from only 2 Chelicerata that were available in 1998. At present, 41 complete chelicerate mt genomes are available in the NCBI-database (, status January 10, 2009). Out of these, only 29 depict the specific arrangement of trnL2 between nad1 and 16S-rRNA (see additional file 1 for an overview of gene arrangements of chelicerate mt genomes). This illustrates that care should be taken when general rules are deduced from limited datasets.
Mt gene arrangements have already provided strong support toward the resolution of several long-standing controversial phylogenetic relationships . Surprisingly, the mt gene order of D. pteronyssinus differs considerably from that of other mites (see additional file 1). Comparing the D. pteronyssinus mt genome to the mt sequence of the oribatid S. magnus , the closest relative of D. pteronyssinus, revealed that only 6 of the 22 gene boundaries in S. magnus are conserved in D. pteronyssinus (Fig. 3).
Extending this analysis to the other Acari mt genomes showed that in several cases the set of neighboured genes that were not separated during the evolution  was greater between members of different superorders (e.g. D. pteronyssinus (Acariformes) and Rhipicephalus sanguineus (Parasitiformes)) than between members of the same superorder (e.g. D. pteronyssinus (Acariformes) and T. urticae (Acariformes)) (Table2). Exclusion of tRNAs in our analysis showed a similar trend, suggesting that protein coding genes were also involved in mt gene rearrangements. These results indicate that mt gene orders seem less useful for deduction of phylogenetic relationships between superorders within the Acari. However, comparing gene order might be more powerful to establish phylogenetic relations within families, as was previously proposed [14, 36]. In the case of the Ixodidae family, it was shown that the division of Prostriata (Ixodes sp.) and Metastriata (R. sanguineus, Amblyomma triguttatum, Haemaphysalis flava) could be linked to mt gene arrangements [37, 38].
Pairwise common interval distance matrix of mt gene orders of Acari*.
D. pteronyssinus (Dp)
Ixodes sp. (I) a
R. sanguineus (Rs) b
V. destructor (Vd)
L. akamushi (La) c
T. urticae (Tu)
W. hayashii (Wh)
U. foilii (Uf)
S. magnus (Sm) d
*bold numbers represent pairwise common interval distances between mt gene orders of Acari (37 genes in total), while italic numbers represent pairwise common interval distances between mt gene orders of Acari without tRNAs (15 genes in total) (L. pallidum and M. occidentalis were excluded from the dataset as the genomes contain duplicated genes whereas other genes are absent). Accession numbers of Acari mt genomes are listed in Table 1.
a similar gene order as L. polyphemus, Ornithodoros sp. and C. capensis
b similar gene order as H. flava and A. triguttatum
c similar gene order as L. deliense
d due to lack of tRNAs only protein coding gene and rRNA order comparison was possible
Base composition and codon usage
The overall AT-content of the mt genome of D. pteronyssinus is 72.6% (Table 1). This is within the range of the average AT-content of Acari mt genomes (74.6 +/- 4.0%). The high AT-content is reflected in the codon usage (Table 3 ) with nucleotides 'A' and 'T' preferred over 'C' and 'G' on the wobble position and the predominant use of codons deficient in 'C' or 'G'. For example, the most frequently used codons are TTT (F) (105 codons per 1000 codons) and TTA (L2) (78 codons per 1000 codons).
Relative synonymous codon usage (RSCU) and number of codons per 1000 codons (NC1000) in the protein coding genes of the mitochondrial genome of D. pteronyssinus.
* RSCU is the number of times a particular codon is observed relative to the number of times a codon would be observed in the absence of any codon usage bias .
Metazoan mt genomes usually present a clear strand bias in nucleotide composition . This is probably due to asymmetric patterns of mutations during transcription and replication when one strand remains transiently in a single-stranded state, making it more vulnerable to DNA damage . However, in the case of mtDNA-replication, this hypothesis is not without controversy [42–45]. The strand bias in nucleotide composition can be measured as GC- and AT-skews ((G%-C%)/(G%+C%) and (A%-T%)/(A%+T%), respectively) . The overall GC- and AT-skews of the J-strand of the D. pteronyssinus mt genome are 0.194 and -0.199, respectively. These are the most extreme values encountered within mite mt genomes up till now (Table 1) and they are reversed compared to the usual strand biases of metazoan mtDNA (negative GC-skew and positive AT-skew for the J-strand). Moreover, a positive GC-skew for mite mt genomes seems to be rare since at present, it was only encountered in Varroa destructor. Although hypothetical, it could be the result of a strand swap of the control region . This region contains all initiation sites for transcription  and an inversion of the control region is expected to produce a global reversal of asymmetric mutational constraints in the mtDNA, resulting with time in a complete reversal of strand compositional bias . The asymmetrical directional mutation pressure is also reflected in the codon usage of genes oriented in opposite directions . Whereas NNG and NNU codons are preferred over NNA and NNC codons on the J-strand, genes on the N-strand show the exact opposite trend (see additional file 2 for an across-strand (N and J) comparison of frequencies of codons ending with the same nucleotide).
Protein coding genes
Nine proteins are encoded by genes on the J-strand (cox1, cox2, cox3, atp6, atp8, nad3, nad4, nad4L, nad5), while four are encoded by genes on the N-strand (nad1, nad6, nad2, cytB). The total length (10,826 bp) and AT-content (71.61%) of the protein-coding genes are within the range of values typical for Acari (10,639.0 +/- 272.0 bp; 74.0 +/- 4.0%, respectively) (Table 1). Compared to other mite mt proteins, cox1, cox2 and cytB are best conserved. On the other hand, atp8, nad6 and nad4L showed lowest similarity values (see additional file 3 for the average identity and similarity % of mt proteins of D. pteronyssinus).
Start and stop codons were determined based on alignments with the corresponding genes and proteins of other mite species. In the case of stop codons, we could also benefit from available expressed sequence tags (ESTs) of D. pteronyssinus (n = 1797) and D. farinae (n = 1735) (Fig. 1) . As for other metazoan mt proteins, unorthodox initiation codons are used  (see additional file 4 for start and stop codons of protein coding genes of Acari mt genomes). Eight genes (cox2, atp6, cox3, nad3, nad6, nad4L, nad4, cytB) use the standard ATG start codon, 3 genes (cox1, nad1, nad2) start with ATA and nad5 initiates with ATT. atp8 most likely starts with codon TTG.
Eleven genes employ a complete translation termination codon, either TAG (cox1, cox3) or TAA (cox2, atp8, atp6, nad1, nad3, nad6, nad4L, nad5, cytB). With the exception of nad3, atp8 and nad4L, D. pteronyssinus ESTs for all these genes confirmed the position of the stop codon (Fig. 1, see additional file 5 for sequences of Dermatophagoides ESTs covering the mt genome of D. pteronyssinus). Berthier et al.  showed that the adjacent genes, nad4L/nad4 and atp8/atp6, were transcribed and translated as a bicistronic mRNA in the model organism Drosophila melanogaster. However, as no ESTs were found that aligned with the nad4L/nad4 and atp8/atp6 gene boundaries, it could not be confirmed whether this was also the case for D. pteronyssinus. Despite its efficiency, the use of sequence alignments to determine the position of stop codons resulted in several cases in overlapping genes. For example, based on a highly conserved tryptophan at the C-terminal end of Acari nad3 proteins, a stop codon was positioned despite the resulting 17 bp overlap with trnaR. The two remaining genes (nad2 and nad4) are likely equipped with a truncated stop codon (T). Polyadenylation of the mRNA is needed in these cases to form a fully functional TAA stop codon . Although speculative, ESTs of D. farinae confirm the truncated stop of nad4 (Fig. 1, see additional file 5).
Fourteen tRNAs are encoded on the J-strand and 8 on the N-strand (Fig. 1). Secondary structures were predicted for all tRNAs (Fig. 4). With the exception of trnS1 (UCU instead of GCU in L. pallidum) and trnP (UGG instead of AGG in S. magnus), all anticodon sequences were identical to those of L. pallidum and S. magnus, the only acariform mites for which tRNA secondary structures have been reported [22, 24]. Usually, T is in the first anticodon position for tRNAs that recognise either four-fold degenerate codon families or NNR codons. G is usually in this position only to specifically recognize NNY-codons . Except for trnM, all of the D. pteronyssinus mt tRNAs follow this pattern. trnM has the anticodon CAT (to recognise both ATG and ATA), which is the case for almost all animal mt systems  (Fig. 4).
Only one tRNA lacks the D-arm: trnS1, as is common for most metazoans. With the exception of trnC, trnV and trnS1, all tRNAs have T-arm variable loops (TV replacement loops) instead of the T-arm. Similar structures were found for tRNAs of L. pallidum  and S. magnus . The absence of the T-arm is a typical feature for tRNAs of Chelicerata belonging to the orders of the Araneae, Scorpiones and Thelyphonida. However, other taxa within the Chelicerata (Amblypygi, Opiliones, Ricinulei, Solifugae and ticks) possess typical metazoan cloverleaf tRNAs . Masta and Boore  suggested a multi-step evolutionary process in an attempt to understand how so many tRNAs in these chelicerate groups could lose their T-arm. According to this speculative theory, changes in mt ribosomes, resulting in the fact that the loss of arms from tRNAs was tolerated [53, 54], and/or changes in specific elongation factors [55–57] are considered as a first step in this process.
Only 7 of the 22 tRNAs have a completely matched 7 bp acceptor stem (trnG, trnW, trnH, trnF, trnE, trnL1 and trnL2). A maximum of 3 mismatches in this stem is found in trnR. In contrast, almost all tRNAs (18) possess a completely matched 5 bp anticodon stem. trnC, trnS1 and trnN have a single mismatch whereas trnY has two mismatches in this stem. All tRNAs, except trnL2, have a symmetric anticodon loop consisting of 2 bp up- and 2 bp downstream of the 3 bp anticodon. The anticodon loop of trnL2 consists of 2 nucleotides preceding the anticodon and 3 nucleotides immediately following it. This kind of aberrant anticodon loops have also been reported for the two-humped camel Camelus bactrianus ferus (trnS2)  and the scorpion Mesobuthus gibbosus (trnH and trnN) . As mentioned before, sequences of some tRNAs overlap with neighbouring genes. The extreme examples are trnR, trnS2 and trnV. trnR overlaps with the adjacent gene nad3 on the same strand for 17 bp at its 3'-end whereas trnS2 overlaps with the adjacent gene trnM on the same strand for 12 bp at its 3'-end. trnV overlaps with the adjacent gene trnP on the opposite strand for 11 bp at its 3'-end and with trnK on the opposite strand for 7 bp at its 5'-start. Despite these overlaps, we consider these genes not likely to be pseudogenes. First of all, their sequence is relatively well conserved when compared to corresponding genes of other Acari. Secondly, besides sequence conservation they depict a conserved secondary structure. Thirdly, an EST [GenBank: CB284825] of the related species D. farinae was found corresponding to the region covering trnR, trnM and trnS2 of D. pteronyssinus indicating that the genes are expressed (see additional file 6 for an alignment of trnR, trnM and trnS2 of D. pteronyssinus with an EST of D. farinae). Finally, and most importantly, stem mismatches and sequence overlap are not uncommon for mt tRNAs of arachnids [13, 60], and are probably repaired by a post-transcriptional editing process [54, 61].
The largest non-coding region (286 bp) is flanked by trnF and trnS1. It is highly enriched in AT (91.61%) and can form stable stem-loop secondary structures. Based on these features, it possibly functions as a control region [20, 62]. With the exception of T. urticae (95.45%), it has the highest AT-content of all Acari mt control regions (Table 1). The position of the non-coding region differs from most insect and arachnid mt genomes, where the region is mostly located in close proximity to 12S-rRNA (, see additional file 1).
Based on the sequence pattern, the control region can be subdivided in a repeat region and a stem-loop region. The first region (11,491–11,528 bp) contains several AT-repeats. In order to verify the exact number of repeats we resequenced this region. For this purpose, two flanking primers, Dp-Ms-F and Dp-Ms-R, were synthesised spanning approximately 700 bp. The PCR product was cloned and ten independent clones were sequenced. This revealed that the number of AT-repeats varied between 7 to 28, suggesting that this domain can be considered as a microsatellite . This is remarkable as a mt microsatellite was never reported before for species belonging to the Chelicerata. Also in metazoan mtDNA such microsatellites are rare and have, to our knowledge, only been reported for butterflies , a dragonfish, Scleropages formosus , a bat, Myotis bechsteinii , a turtle, Pelomedusa subrufa  and several seal species [68–70].
The second region (11,529–11,768 bp) holds two short palindromic sequences, TACAT and ATGTA, which are conserved in mt genomes of mammals  and fishes [65, 72]. They can form a stable stem-loop structure (Fig. 5-A2), which might be involved as a recognition site for the arrest of J-strand synthesis . Near this region other stem-loop structures could be folded (Fig. 5-A) but none of them had flanking sequences similar to those that are conserved in the control region of the mt genome of insects  and metastriate ticks .
As described before, four other stretches of non-coding nucleotides were found outside the control region. These short sequences can fold into stable stem-loop structures (Fig. 5-B, C, D) which may function as splicing recognition sites during processing of the transcripts .
12S-rRNA and 16S-rRNA are located on the J-strand. This does not coincide with their position in most Chelicerata where they are located on the N-strand (see additional file 1). The AT-contents of both genes are comparable (72.9% and 76.1% for the 12S- and 16S-RNA, respectively) and are within the range of rRNAs of other Acari (76.5 +/- 4.2%; 78.0 +/- 4.1%, respectively). The sizes of the rRNAs (665 bp and 1078 bp) are slightly larger than those of other acariform mite rRNAs (626.0 +/- 29.9 bp and 1018.5 +/- 21.3 bp) but are shorter than those found in the Parasitiformes (706.0 +/- 17.5 bp and 1207.3 +/- 31.4 bp) (Table 1).
The 12S-rRNA and 16S-rRNA genes of Leptotrombidium species (Acariformes: Trombidiformes: Trombiculidae) are 23.4% and 23.5% shorter than their counterparts in Drosophila yakuba. This substantial reduction is mainly caused by the loss of stem-loop structures at the 5'-end of the rRNA genes . To identify whether similar domains are absent in the rRNAs of D. pteronyssinus, we constructed their secondary structures (Fig. 6 [75, 76]). This revealed that the D. pteronyssinus 12S-rRNA indeed lacks similar stem-loops as L. pallidum, compared to D. yakuba. The structure also revealed 1 additional stem-loop (stem-loop 1) not present in 12S-rRNA of L. pallidum. Like in L. pallidum, one stem-loop replaces three stem-loops (24, 25 and 26) whereas another replaces a region of four stem-loops (39, 40, 41 and 42) of the D. yakuba 12S-rRNA . Based on the modelled structure in combination with an alignment of other acariform 12S-rRNAs, the greatest sequence conservation was found in the loop region of stem-loops 21 and 27 and the region between stem-loops 48 and 50.
In analogy to the 16S-rRNA gene of L. pallidum, the main deletions of the D. pteronyssinus 16S-rRNA are located at the 5'-end. With the exception of D19, all stem-loops of L. pallidum are present in D. pteronyssinus. We also discovered three additional stem-loops (C1, E2 and E19) which are absent in the 16S-rRNA of L. pallidum. The 3'-end of the 16S-rRNA structure is best conserved compared to other acariform 16S-rRNAs. This is in agreement with the idea that this region is the main component of the peptidyl-transferase centre, and as such most vulnerable to mutations . Recently, the 12S-rRNA and 16S-rRNA secondary structures of S. magnus have been published . The 12S-rRNA structure of S. magnus has 5 extra stem-loops (2, 4, 5, 40 and 42) compared to the one of D. pteronyssinus whereas the 16S-rRNA lacks 6 stem-loops (D4, D16, E1, E2, E19 and G9) and has 5 stem-loops (cd1, D1, D17, D19, G13) not present in the 16S-rRNA of D. pteronyssinus.
It is still an open question how relatively well-conserved structures such as rRNAs can dramatically decrease in size while remaining functional. Wolstenholme et al.  and Masta  suggested a correlation between the occurrence of truncated rRNAs (compared to Drosophila) and the loss of the T-arm in tRNAs. The coincidence of short rRNAs and missing T-arms in tRNAs was also observed in S. magnus, L. pallidum and D. pteronyssinus. Other acariform mites like T. urticae, Ascoschoengastia sp. and Walchia hayashii also exhibit short rRNAs (Table 1) and the prediction of their tRNA secondary structures could further support this hypothesis. However, examples contradicting this hypothesis also occur e.g. pulmonate gastropods with tRNAs lacking T-arms have no truncated rRNAs. Therefore, it remains possible that truncation of both tRNAs and rRNA genes only reflects an independent trend towards minimisation of the mt genome as suggested by Yamazaki et al. .
A phylogenetic tree was constructed based on nucleotide and amino acid sequences from all mt protein coding genes of Acari. The ILD-test  indicated a significant incongruence (P = 0.01) among data set partitions for nucleotide alignments and low congruence (P = 0.07) among data set partitions for amino acid alignments. A considerable debate exists on the utility of this test [79–84]. However, the principle of Kluge  implies that all data should always be included in a combined analysis for any phylogenetic problem and therefore we combined data partitions for both amino acid and nucleotide alignments for phylogenetic analysis. A maximum parsimony (MP) analysis based on nucleotide alignments (data not shown) grouped V. destructor (Parasitiformes) within the Acariformes, close to D. pteronyssinus. This is in contrast with the generally accepted view on the phylogeny of the Acariformes and Parasitiformes [2, 3]. As mentioned before, V. destructor and D. pteronyssinus both have a reversal of asymmetrical mutation pattern. When such reversals occurred independently, D. pteronyssinus and V. destructor could have acquired a similar base composition and as a consequence group together due to the long-branch attraction (LBA) phenomenon [40, 86]. Model-based methods such as maximum likelihood (ML) and Bayesian inference (BI) are less sensitive to LBA [40, 87] and were for this reason considered for phylogenetic analysis.
ML and BI analysis performed on the amino acid data set resolved trees with an identical topology (Fig. 7-A) in which D. pteronyssinus clusters with S. magnus, forming a sistergroup of the Trombidiformes. This is in agreement with the most recent views on the classification of the Acariformes [2–4]. The nucleotide data set resulted in similar trees, confirming the evolutionary position of D. pteronyssinus (Fig. 7-B). The only major inconsistency over the trees was the position of T. urticae. Although this species is generally considered as a member of the Trombidiformes [2–4], it was clustered with the sarcoptiform mites D. pteronyssinus and S. magnus in the trees based on the nucleotide dataset. (Fig. 7-B). However, the position in the different trees is questionable as it is supported by low bootstrap values/Bayesian posterior probabilities (Fig. 7-A/B). Adding additional mt genome data from closely related taxa of T. urticae and from taxa located between T. urticae and Trombiculidae would probably position T. urticae with higher support values within the Trombidiformes.
In the trees based on the nucleotide dataset, H. flava is, compared to A. triguttatum, evolutionary closer related to R. sanguineus while in the trees based on the amino acid dataset this is the opposite. However, as the clustering of H. flava and R. sanguineus is in agreement with the most recent views on the classification of the Ixodida [88, 89], we consider the nucleotide topology as the most correct one. Murrell et al.  considers the Parasitiformes to be paraphyletic with respect to the Opilioacariformes, but as there are no complete mt genomes of Opilioacariformes available, we were not able to verify this hypothesis.
This is the first description of a complete mt genome of a species belonging to the Astigmata, a cohort within the Sarcoptiformes. Although the length, gene and AT-content are similar to other Acari mtDNA, the mt genome of D. pteronyssinus exhibits some interesting features. The gene order of D. pteronyssinus is completely different from that of other Acari mt genomes. Gene order comparison indicated that mt gene orders seem less useful for deduction of phylogenetic relationships between superorders within the Acari. GC- and AT-skews of the J-strand were very large and reversed as compared to those found in most metazoan mtDNA.
Compared to parasitiform mites, both D. pteronyssinus rRNAs were considerably shorter and almost all transfer RNAs lacked the T-arm. It would be interesting to investigate whether the occurrence of truncated rRNAs and the loss of the T-arm in tRNAs are correlated or just a trend toward minimisation of the mt genome. Finally, phylogenetic analysis using concatenated mt gene sequences succeeded in recovering Acari relationships concordant with traditional views of phylogeny of Acari.
Upon arrival in the laboratory, mites were identified as D. pteronyssinus by J. Witters (ILVO, Belgium) and F. Th. M. Spieksma (Laboratory of Aerobiology, LUMC, The Netherlands) using morphological characteristics. To back up this identification, molecular techniques were applied. For this purpose DNA was extracted and used as a template for PCR. Primers 12SID-F and 12SID-R (see additional file 7 for primer sequences) successfully amplified a 316 bp fragment. BLASTn searches against non-redundant nucleotide sequences using the amplified fragment as query resulted in a perfect match with a mt 12S-rRNA sequence of D. pteronyssinus [GenBank: AF529911].
Mite strain, mass rearing and isolation
The initial D. pteronyssinus culture was provided by D. Bylemans (Janssen Pharmaceutica, Belgium). Mites were cultured on a 1:1 mixture of Premium Gold (Vitacraft, Germany) and beard shavings at 75% R.H., 25°C and permanent dark conditions [91, 92]. Mites were isolated from the colony using a modified heat-escape technique [93, 94]. Briefly, mite cultures were transferred to small plastic petri dishes (75 mm in diameter, 28 mm high) with a lid on top. These dishes were placed in the dark on a hot plate set at 45°C (Bekso, Belgium). After 15–20 minutes the mites moved away from the heat source, formed groups on the lid of the petri dish and could be collected using a fine hair brush.
Approximately 1000 D. pteronyssinus mites were collected in an Eppendorf tube and were ground in 800 μl SDS-lysis buffer (400 mM NaCl, 200 mM TRIS, 10 mM EDTA, 2% SDS) using a small sterile plastic pestle (Eppendorf, Germany). After incubation for 30 min at 60°C under continuous rotation, a standard phenol-chloroform extraction was performed . Total genomic DNA was precipitated with 0.7 volumes of isopropanol at 4°C for 1 hour, centrifuged for 45 minutes at 21,000 × g and washed with 70% ethanol. Precipitated DNA was resolved in 50 μl 0.1 M Tris pH 8.2.
Standard PCR (amplicon < 500 bp) was performed in 50 μl volumes (38.5 μl double-distilled water; 5 μl buffer; 2 mM MgCl2; 0.2 mM dNTP-mix; 0.2 μM of each primer; 1 μl template DNA and 0.5 μl Taq polymerase (Invitrogen, Belgium). PCR conditions were as follows: 2' 94°C, 35 × (20" 92°C, 30" 53°C, 1' 72°C) and 2' 72°C. The annealing time was extended to 1 minute and the primer concentration was increased to 2 μM when degenerate primers were used. Long PCR (amplicon > 500 bp) was performed with the Expand Long Range Kit (Roche, Switzerland) in 50 μl volumes (28.5 μl double-distilled water; 10 μl buffer; 0.5 mM dNTP-mix; 0.3 μM of each primer; 4 μl 100% DMSO; 1 μl template DNA and 1 μl enzyme-mix). PCR conditions were: 2' 94°C, 10 × (10" 92°C, 20" at a temperature that varies depending on the primers, 1'/kb 58°C), 25 × (10" 92°C, 20" at a temperature that varies depending on the primers, 1'/kb 58°C with 20" added for every consecutive cycle) and 7' 58°C. All PCR products were separated by electrophoresis on a 1% agarose gel and visualised by EtBr staining. Fragments (amplicon < 1000 bp) of interest were excised from gel, purified with the QIAquick PCR Purification Kit (Qiagen, Belgium) and cloned into the pGEM-T vector (Promega, Belgium). After heat-shock transformation of E. coli (DH5α) cells, plasmid DNA was obtained by miniprep and inserted fragments were sequenced with SP6 and T7-primers. Long PCR products were sequenced by primer-walking. All sequencing reactions were performed by AGOWA sequencing service.
Amplification of the mt genome
Primers COXI-F and 12S-R, based on partial D. pteronyssinus cox1 and 12S-rRNA sequences [GenBank: AY525570 and AF529911, respectively] (see additional file 7 for primer sequences), successfully amplified a 4.6 kb sequence of the mt genome of D. pteronyssinus. Degenerate primers CYTB-F-Deg and CYTB-R-Deg (see additional file 7 for primer sequences), designed on conserved regions of Acari cytB, amplified a partial cytB sequence from D. pteronyssinus. A specific primer COXI-R, designed from the 3' end of the 4.5 kb sequence in combination with the primer CYTB-F, designed from the partial cytB sequence, successfully amplified a 2.2 kb sequence. Another primer CYTB-R, designed from the 5'-end of this 2.2 kb sequence, in combination with the primer 12S-F successfully amplified a 8.6 kb sequence, making the mt genome sequence complete.
Annotation and bioinformatics analysis
The complete genomic sequence was assembled and annotated using VectorNTI (Invitrogen, Belgium) according to Masta and Boore . Open reading frames (ORFs) were identified with the program Getorf from the EMBOSS-package . The obtained ORFs were used as query in BLASTp  searches against the non-redundant protein database at NCBI. Two large non-protein-coding regions were candidates for the rRNAs (16S and 12S respectively). The boundaries were identified based on alignments and secondary structures of rRNA genes of other mite species. Sixteen of the 22 tRNAs were identified by tRNA-scan SE  with a cove cutoff score of 0.1 and the tRNA-model set to "nematode mito". The remaining tRNAs (trnM, trnV, trnY, trnS1, trnI, trnC) were determined in the unannotated regions by sequence similarity to tRNAs of other mite species. In order to obtain additional information on mt gene boundaries, BLASTn  searches of D. pteronyssinus tRNA, rRNA and protein encoding nucleotide sequences were carried out against ESTs  restricted to Dermatophagoides sequences (n = 3532). ESTs with statistically significant matches (E-value cutoff: 0.1) were collected, checked for vector contamination and aligned by Clustal W  as implemented in BioEdit 7.0.1  against the appropriate nucleotide sequence of D. pteronyssinus. MatGAT 2.02 was used to calculate similarity and identity values  of mt proteins. The identification of gene subsets that appear consecutively in different genomes was performed by common interval distance analysis using CREx  (see additional file 8 for input data of the CREx program).
Construction of secondary structures of RNAs and non-coding regions
Secondary structures of tRNAs were determined following the method of Masta and Boore . Secondary structures of tRNAs were drawn with CorelDraw 12.0 (Corel Corporation, Canada). The rRNA genes of D. pteronyssinus were aligned with those of other Acariformes and conserved areas were identified. These regions were mapped on the published structures of L. pallidum rRNA . Regions lacking significant homology were folded using Mfold . Secondary structures of rRNAs were drawn using the RnaViz2 program  and afterwards modified with CorelDraw 12.0 (Corel Corporation, Canada). Secondary structures of non-coding regions were folded using Mfold . When multiple secondary structures were possible, the most stable (lowest free energy (-ΔG)) one was preferred. Drawing and editing of these structures was done in a similar way as for rRNA secondary structures.
Rolling circle amplification and restriction enzyme digestion
Extraction and rolling circle amplification of the mtDNA of D. pteronyssinus was done according to Van Leeuwen et al. . Rolling circle amplified mtDNA was digested with two enzymes (XmnI and EcoRI; New England Biolabs) following the manufacturer's instructions. Restriction digests were fractionated by agarose gel electrophoresis as described before.
Sequence data were obtained from 21 Acari species (for GenBank accession numbers see Table 1) and two outgroup taxa (Limulus polyphemus [GenBank: NC_003057] and Locusta migratoria [Genbank: NC_001712]). Only mite species with a completely sequenced mt genome were selected. Alignments from all mt protein-coding genes were used in phylogenetic analysis. Amino acid sequences and nucleotide sequences were aligned by Clustal W  as implemented in BioEdit 7.0.1 . The nucleotide alignment was generated based on the protein alignment using codon alignment. Ambiguously aligned parts were omitted from the analysis by making use of Gblocks 0.91b , with default block parameters except for changing "allowed gap positions" to "with half". Abascal et al.  recently presented evidence that some insects and ticks use a modified mitochondrial code, with AGG coding for lysine rather than serine as in the standard invertebrate mitochondrial code. As 10 out of 20 Acari species in our dataset are ticks all positions aligning to AGG codons in the final amino acid alignment were removed.
For the nucleotide alignments the "codons" option was used in Gblocks 0.91b . Due to the results of a saturation analysis  on single codon positions, implemented in DAMBE 4.2.13 , third codon positions were eliminated from the nucleotide alignment. An incongruence length difference test (ILD-test)  as implemented in PAUP* (version 4.0b10; ) was used to assess congruence among gene partitions.
Model selection was done with ProtTest 1.4  for amino acid sequences and with Modeltest 3.7  for nucleotide sequences. According to the Akaike information criterion, the mtART+G+I+F model was optimum for phylogenetic analysis with amino acid alignments and the GTR+I+G model was optimal for analysis with nucleotide alignments.
Two different analyses were performed. (1) Maximum likelihood (ML) analysis was performed using Treefinder , bootstrapping with 1000 pseudoreplicates (2) Bayesian inference (BI) was done with MrBayes 3.1.2 . As the mtART model is not implemented in the current version of MrBayes, the mtREV+G+I model was used for phylogenetic analysis with the amino acid alignment. Four chains ran for 1,000,000 generations, while tree sampling was done every 100 generations. Burnin was calculated when the average standard deviation of split frequencies had declined to < 0.01. The remaining trees were used to calculate Bayesian posterior probabilities (BPP).
small (12S) rRNA subunit (gene)
large (16S) rRNA subunit (gene)
atp6 and 8:
ATPase subunit 6 and 8
cytochrome oxidase subunits I-III
dihydrouridine-arm of a tRNA secondary structure
expressed sequence tag
large non-coding region/control region
nad1-6 and nad4L:
NADH dehydrogenase subunits 1–6 and 4L
open reading frame
polymerase chain reaction
TΨC-arm of a tRNA secondary structure
(where X is replaced by a one letter amino acid code for the corresponding amino acid, with L1 = CUN
L2 = UUR S1 = AGN S2 = UCN),:
TVL is a Postdoctoral Fellow of the Research Foundation Flanders (Belgium). BV had a postdoctoral grant from Ghent University (BOF-01P04105).
Department of Crop Protection, Faculty of Bioscience Engineering, Ghent University
Department of Molecular Biotechnology, Faculty of Bioscience Engineering, Ghent University
VIB Department of Plant Systems Biology, Ghent University
Harvey MS: The neglected cousins: What do we know about the smaller Arachnid orders?J Arachnol 2002, 30:357–372.
Dunlop JA, Alberti G: The affinities of mites and ticks: a review.J Zool Syst Evol Res 2008, 46:1–18.
O'Connor BM: Phylogenetic relationships among higher taxa in the Acariformes, with particular reference to the Astigmata.Acarology VI(Edited by: Griffiths DA, Bowman CE). Chichester: Ellis–Horwood Ltd 1984, 1:19–27.
Spieksma FTM: Domestic mites from an acarologic perspectives.Allergy 1997, 52:360–368.PubMed
Voorhorst R, Spieksma FT, Varekamp H, Leupen MJ, Lyklema AW: House–dust mite (Dermatophagoides pteronyssinus) and allergens it produces. Identity with house–dust allergen.J Allergy 1967, 39:325–339.
van Bronswijk JEMA: House dust biology. The Netherlands: NIB Publishers 1981.
Arlian LG, Platts–Mills TAE: The biology of dust mites and the remediation of mite allergens in allergic disease.J Allergy Clin Immunol 2001, 107:S406–S413.PubMed
Janson C, Anto J, Burney P, Chinn S, de Marco R, Heinrich J, Jarvis D, Kuenzli N, Leynaert B, Luczynska C, Neukirch F, Svanes C, Sunyer J, Wjst M: The European Community Respiratory Health Survey: what are the main results so far?Eur Respir J 2001, 18:598–611.PubMed
Dowton M, Castro LR, Austin AD: Mitochondrial gene rearrangements as phylogenetic characters in the invertebrates: the examination of genome 'morphology'.Invertebr Syst 2002, 16:345–356.
Boore JL, Macey JR, Medina M: Sequencing and comparing whole mitochondrial genomes of animals.Molecular Evolution: Producing the Biochemical Data, Part B San Diego: Elsevier Academic Press Inc 2005, 395:311–348.
Boore JL: The use of genome–level characters for phylogenetic reconstruction.Trends Ecol Evol 2006, 21:439–446.PubMed
Masta SE, Boore JL: Parallel evolution of truncated transfer RNA genes in arachnid mitochondrial genomes.Mol Biol Evol 2008, 25:949–959.PubMed
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:534–557.
Castro LR, Dowton M: Mitochondrial genomes in the Hymenoptera and their utility as phylogenetic markers.Syst Entomol 2007, 32:60–69.
Dekeyser MA: Acaricide mode of action.Pest Manag Sci 2005, 61:103–110.PubMed
Van Leeuwen T, Vanholme B, Van Pottelberge S, Van Nieuwenhuyse P, Nauen R, Tirry L, Denholm I: Mitochondrial heteroplasmy and the evolution of insecticide resistance: Non–Mendelian inheritance in action.Proc Natl Acad Sci USA 2008, 105:5980–5985.PubMedPubMed Central
Domes K, Maraun M, Scheu S, Cameron SL: The complete mitochondrial genome of the sexual oribatid mite Steganacarus magnus: genome rearrangements and loss of tRNAs.BMC Genomics 2008, 9:532.PubMedPubMed Central
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.PubMed
Shao RF, 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.PubMed
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.PubMed
Shao RF, Barker SC, Mitani H, Aoki Y, Fukunaga M: Evolution of duplicate control regions in the mitochondrial genomes of metazoa: A case study with Australasian Ixodes ticks.Mol Biol Evol 2005, 22:620–629.PubMed
Mitani H, Talbert A, Fukunaga M: New World relapsing fever Borrelia found in Ornithodoros porcinus ticks in central Tanzania.Microbiol Immunol 2004, 48:501–505.PubMed
Gorrochotegui–Escalante N, Black WC: Amplifying whole insect genomes with multiple displacement amplification.Insect Mol Biol 2003, 12:195–200.PubMed
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.PubMed
Lavrov DV, Boore JL, Brown WM: The complete mitochondrial DNA sequence of the horseshoe crab Limulus polyphemus.Mol Biol Evol 2000, 17:813–824.PubMed
Fahrein K, Talarico G, Braband A, Podsiadlowski L: The complete mitochondrial genome of Pseudocellus pearsei (Chelicerata: Ricinulei) and a comparison of mitochondrial gene rearrangements in Arachnida.BMC Genomics 2007, 8:386.PubMedPubMed Central
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.PubMedPubMed Central
Boore JL, Lavrov DV, Brown WM: Gene translocation links insects and crustaceans.Nature 1998, 392:667–668.PubMed
Berard S, Bergeron A, Chauve C, Paul C: Perfect sorting by reversals is not always difficult.IEEE/ACM Trans Comput Biol Bioinform 2007, 4:4–16.PubMed
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:1293–1304.PubMed
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.PubMed
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.PubMed
Sharp PM, Cowe E, Higgins DG, Shields DC, Wolfe KH, Wright F: Codon usage patterns in Escherichia coli, Bacillus subtilis, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Drosophila melanogaster and Homo sapiens – A review on the considerable within species diversity.Nucleic Acids Res 1988, 16:8207–8211.PubMedPubMed Central
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.PubMed
Francino MP, Ochman H: Strand asymmetries in DNA evolution.Trends Genet 1997, 13:240–245.PubMed
Clayton DA: Replication and transcription of vertebrate mitochondrial DNA.Annu Rev Cell Biol 1991, 7:453–478.PubMed
Yang MY, Bowmaker M, Reyes A, Vergani L, Angeli P, Gringeri E, Jacobs HT, Holt IJ: Biased incorporation of ribonucleotides on the mitochondrial L–strand accounts for apparent strand–asymmetric DNA replication.Cell 2002, 111:495–505.PubMed
Yasukawa T, Yang MY, Jacobs HT, Holt IJ: A bidirectional origin of replication maps to the major noncoding region of human mitochondrial DNA.Mol Cell 2005, 18:651–662.PubMed
Brown TA, Cecconi C, Tkachuk AN, Bustamante C, Clayton DA: Replication of mitochondrial DNA occurs by strand displacement with alternative light–strand origins, not via a strand–coupled mechanism.Genes Dev 2005, 19:2466–2476.PubMedPubMed Central
Perna NT, Kocher TD: Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes.J Mol Evol 1995, 41:353–358.PubMed
Taanman JW: The mitochondrial genome: structure, transcription, translation and replication.Biochim Biophys Acta, Bioenerg 1999, 1410:103–123.
Carapelli A, Comandi S, Convey P, Nardi F, Frati F: The complete mitochondrial genome of the antartic springtail Cryptopygus antarcticus (Hexapoda: Collembola).BMC Genomics 2008, 9:315.PubMedPubMed Central
Gissi C, Pesole G: Transcript mapping and genome annotation of ascidian mtDNA using EST data.Genome Res 2003, 13:2203–2212.PubMedPubMed Central
Berthier F, Renaud M, Alziari S, Durand R: RNA mapping on Drosophila mitochondrial DNA – precursors and template strands.Nucleic Acids Res 1986, 14:4519–4533.PubMedPubMed Central
Ojala D, Montoya J, Attardi G: Transfer–RNA punctuation model of RNA processing in human mitochondria.Nature 1981, 290:470–474.PubMed
Boore JL: The complete sequence of the mitochondrial genome of Nautilus macromphalus (Mollusca: Cephalopoda).BMC Genomics 2006, 7:182.PubMedPubMed Central
Wolstenholme DR, Macfarlane JL, Okimoto R, Clary DO, Wahleithner JA: Bizarre transfer–RNAs inferred from DNA–sequences of mitochondrial genomes of nematode worms.Proc Natl Acad Sci USA 1987, 84:1324–1328.PubMedPubMed Central
Masta SE: Mitochondrial sequence evolution in spiders: Intraspecific variation in tRNAs lacking the T psi C arm.Mol Biol Evol 2000, 17:1091–1100.PubMed
Ohtsuki T, Watanabe Y, Takemoto C, Kawai G, Ueda T, Kita K, Kojima S, Kaziro Y, Nyborg J, Watanabe K: An "elongated" translation elongation factor Tu for truncated tRNAs in nematode mitochondria.J Biol Chem 2001, 276:21571–21577.PubMed
Arita M, Suematsu T, Osanai A, Inaba T, Kamiya H, Kita K, Sisido M, Watanabe Y, Ohtsuki T: An evolutionary 'intermediate state' of mitochondrial translation systems found in Trichinella species of parasitic nematodes: co–evolution of tRNA and EF–Tu.Nucleic Acids Res 2006, 34:5291–5299.PubMedPubMed Central
Ohtsuki T, Watanabe Y: T–armless tRNAs and elongated elongation factor Tu.IUBMB Life 2007, 59:68–75.PubMed
Cui P, Ji R, Ding F, Qi D, Gao HW, Meng H, Yu J, Hu SN, Zhang HP: A complete mitochondrial genome sequence of the wild two–humped camel (Camelus bactrianus ferus): an evolutionary history of camelidae.BMC Genomics 2007, 8:241.PubMedPubMed Central
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.PubMed
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.PubMed
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:13738–13742.PubMedPubMed Central
Zhang DX, Hewitt GM: Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies.Biochem Syst Ecol 1997, 25:99–120.
Goldstein DB, Schlotterer C: Microsatellites: Evolution and Applications. USA: Oxford University Press 1999.
Cameron SL, Whiting MF: The complete mitochondrial genome of the tobacco hornworm, Manduca sexta, (Insecta: Lepidoptera: Sphingidae), and an examination of mitochondrial gene variability within butterflies and moths.Gene 2008, 408:112–123.PubMed
Yue GH, Liew WC, Orban L: The complete mitochondrial genome of a basal teleost, the Asian arowana (Scleropages formosus, Osteoglossidae).BMC Genomics 2006, 7:242.PubMedPubMed Central
Mayer F, Kerth G: Microsatellite evolution in the mitochondrial genome of Bechstein's bat (Myotis bechsteinii).J Mol Evol 2005, 61:408–416.PubMed
Zardoya R, Meyer A: Cloning and characterization of a microsatellite in the mitochondrial control region of the African side–necked turtle, Pelomedusa subrufa.Gene 1998, 216:149–153.PubMed
Arnason U, Johnsson E: The complete mitochondrial DNA sequence of the harbor seal, Phoca vitulina.J Mol Evol 1992, 34:493–505.PubMed
Arnason U, Gullberg A, Johnsson E, Ledje C: The nucleotide sequence of the mitochondrial DNA molecule of the gray seal Halichoerus grypus, and a comparison with mitochondrial sequences of other true seals.J Mol Evol 1993, 37:323–330.PubMed
Hoelzel AR, Hancock JM, Dover GA: Generations of VNTRS and heteroplasmy by sequence turnover in the mitochondrial control region of 2 elephant seal species.J Mol Evol 1993, 37:190–197.PubMed
Saccone C, Pesole G, Sbisa E: The main regulatory region of mammalian mitochondrial DNA – Structure function model and evolutionary pattern.J Mol Evol 1991, 33:83–91.PubMed
Zardoya R, Meyer A: The complete nucleotide sequence of the mitochondrial genome of the lungfish (Protopterus dolloi) supports its phylogenetic position as a close relative of land vertebrates.Genetics 1996, 142:1249–1263.PubMedPubMed Central
He Y, Jones J, Armstrong M, Lamberti F, Moens M: The mitochondrial genome of Xiphinema americanum sensu stricto (Nematoda: Enoplea): Considerable economization in the length and structural features of encoded genes.J Mol Evol 2005, 61:819–833.PubMed
Shao RF, 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.PubMed
De Rijk P, Robbrecht E, De Hoog S, Caers A, Peer Y, De Wachter R: Database on the structure of large subunit RNA.Nucleic Acids Res 1999, 27:174–178.PubMedPubMed Central
Peer Y, Caers A, De Rijk P, De Wachter R: Database on the structure of small ribosomal subunit RNA.Nucleic Acids Res 1998, 26:179–182.PubMedPubMed Central
Yamazaki N, Ueshima R, Terrett JA, Yokobori S, Kaifu M, Segawa R, Kobayashi T, Numachi K, Ueda T, Nishikawa K, Watanabe K, Thomas RH: Evolution of pulmonate gastropod mitochondrial genomes: Comparisons of gene organizations of Euhadra, Cepaea and Albinaria and implications of unusual tRNA secondary structures.Genetics 1997, 145:749–758.PubMedPubMed Central
Farris JS, Kallersjo M, Kluge AG, Bult C: Constructing a significance test for incongruence.Syst Biol 1995, 44:570–572.
Cunningham CW: Can three incongruence tests predict when data should be combined?Mol Biol Evol 1997, 14:733–740.PubMed
Dolphin K, Belshaw R, Orme CDL, Quicke DLJ: Noise and incongruence: Interpreting results of the incongruence length difference test.Mol Phylogenet Evol 2000, 17:401–406.PubMed
Yoder AD, Irwin JA, Payseur BA: Failure of the ILD to determine data combinability for slow loris phylogeny.Syst Biol 2001, 50:408–424.PubMed
Barker FK, Lutzoni FM: The utility of the incongruence length difference test.Syst Biol 2002, 51:625–637.PubMed
Dowton M, Austin AD: Increased congruence does not necessarily indicate increased phylogenetic accuracy – The behavior of the incongruence length difference test in mixed–model analyses.Syst Biol 2002, 51:19–31.PubMed
Hipp AL, Hall JC, Sytsma KJ: Congruence versus phylogenetic accuracy: Revisiting the incongruence length difference test.Syst Biol 2004, 53:81–89.PubMed
Kluge AG: A concern for evidence and a phylogenetic hypothesis of relationships among Epicrates (Boidae, Serpentes).Syst Zool 1989, 38:7–25.
Felsenstein J: Cases in which parsimony or compatibility methods will be positvely misleading.Syst Zool 1978, 27:401–410.
Cunningham CW, Zhu H, Hillis DM: Best–fit maximum–likelihood models for phylogenetic inference: Empirical tests with known phylogenies.Evolution 1998, 52:978–987.
Barker SC, Murrell A: Phylogeny, evolution and historical zoogeography of ticks: a review of recent progress.Exp Appl Acarol 2002, 28:55–68.PubMed
Klompen H, Lekvelshvili M, Black WC: Phylogeny of parasitiform mites (Acari) based on rRNA.Mol Phylogenet Evol 2007, 43:936–951.PubMed
Murrell A, Dobson SJ, Walter DE, Campbell NJH, Shao RF, Barker SC: Relationships among the three major lineages of the Acari (Arthropoda: Arachnida) inferred from small subunit rRNA: paraphyly of the parasitiformes with respect to the opilioacariformes and relative rates of nucleotide substitution.Invertebr Syst 2005, 19:383–389.
Miyamoto J, Ishii A, Sasa M: Successful method for mass–culture of house dust mite, Dermatophagoides pteronyssinus (Trouessart, 1897).Jpn J Exp Med 1975, 45:133–138.PubMed
Kalpaklioglu AF, Ferizli AG, Misirligil Z, Demirel YS, Gurbuz L: The effectiveness of benzyl benzoate and different chemicals as acaricides.Allergy 1996, 51:164–170.PubMed
Bischoff E: Méthodes actuelles de quantification des acariens dans l'habitat.Rev Fr Allergol 1988, 28:115–122.
Hart B: Ecology and biology of allergenic mites.Mites and Allergic disease(Edited by: Guerin). Belgium: Groeninge Press 1990, 135–152.
Sambrook J, Russel D: Molecular Cloning: a Laboratory Manual2 Edition United states: Cold Spring Harbor Laboratory Press 1987.
Rice P, Longden I, Bleasby A: EMBOSS: The European Molecular Biology Open Software Suite.Trends Genet 2000, 16:276–277.PubMed
Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI–BLAST: a new generation of protein database search programs.Nucleic Acids Res 1997, 25:3389–3402.PubMedPubMed Central
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.PubMedPubMed Central
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 Res 1994, 22:4673–4680.PubMedPubMed Central
Hall T: BioEdit: a user–friendly biological sequence alignment–editor and analysis program for Windows 95/98/NT.Nucleic Acids Symp Ser 1999, 41:95–98.
Campanella JJ, Bitincka L, Smalley J: MatGAT: An application that generates similarity/identity matrices using protein or DNA sequences.BMC Bioinformatics 2003, 4:29.PubMedPubMed Central
Bernt M, Merkle D, Ramsch K, Fritzsch G, Perseke M, Bernhard D, Schlegel M, Stadler PF, Middendorf M: CREx: inferring genomic rearrangements based on common intervals.Bioinformatics 2007, 23:2957–2958.PubMed
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction.Nucleic Acids Res 2003, 31:3406–3415.PubMedPubMed Central
De Rijk P, Wuyts J, De Wachter R: RnaViz 2: an improved representation of RNA secondary structure.Bioinformatics 2003, 19:299–300.PubMed
Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis.Mol Biol Evol 2000, 17:540–552.PubMed
Abascal F, Posada D, Knight RD, Zardoya R: Parallel evolution of the genetic code in arthropod mitochondrial genomes.PLoS Biol 2006, 4:711–718.
Xia XH, Xie Z, Salemi M, Chen L, Wang Y: An index of substitution saturation and its application.Mol Phylogenet Evol 2003, 26:1–7.PubMed
Xia X, Xie Z: DAMBE: Software package for data analysis in molecular biology and evolution.J Hered 2001, 92:371–373.PubMed
Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods), version 4.0b10. Sunderland, Massachusetts: Sinauer Asociates 2003.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best–fit models of protein evolution.Bioinformatics 2005, 21:2104–2105.PubMed
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution.Bioinformatics 1998, 14:817–818.PubMed
Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics.BMC Evol Biol 2004, 4:18.PubMedPubMed Central
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.