- Open Access
The complete mitochondrial genome of the scab mite Psoroptes cuniculi (Arthropoda: Arachnida) provides insights into Acari phylogeny
Parasites & Vectorsvolume 7, Article number: 340 (2014)
Limited available sequence information has greatly impeded population genetics, phylogenetics and systematics studies in the subclass Acari (mites and ticks). Mitochondrial (mt) DNA is well known to provide genetic markers for investigations in these areas, but complete mt genomic data have been lacking for many Acari species. Herein, we present the complete mt genome of the scab mite Psoroptes cuniculi.
P. cuniculi was collected from a naturally infected New Zealand white rabbit from China and identified by morphological criteria. The complete mt genome of P. cuniculi was amplified by PCR and then sequenced. The relationships of this scab mite with selected members of the Acari were assessed by phylogenetic analysis of concatenated amino acid sequence datasets by Bayesian inference (BI), maximum likelihood (ML) and maximum parsimony (MP).
This mt genome (14,247 bp) is circular and consists of 37 genes, including 13 genes for proteins, 22 genes for tRNA, 2 genes for rRNA. The gene arrangement in mt genome of P. cuniculi is the same as those of Dermatophagoides farinae (Pyroglyphidae) and Aleuroglyphus ovatus (Acaridae), but distinct from those of Steganacarus magnus (Steganacaridae) and Panonychus citri (Tetranychidae). Phylogenetic analyses using concatenated amino acid sequences of 12 protein-coding genes, with three different computational algorithms (BI, ML and MP), showed the division of subclass Acari into two superorders, supported the monophylies of the both superorders Parasitiformes and Acariformes; and the three orders Ixodida and Mesostigmata and Astigmata, but rejected the monophyly of the order Prostigmata.
The mt genome of P. cuniculi represents the first mt genome of any member of the family Psoroptidae. Analysis of mt genome sequences in the present study has provided new insights into the phylogenetic relationships among several major lineages of Acari species.
The metazoan mitochondrial (mt) genome possess a circular double strand DNA that varies in size from 15 to 20 kb, generally encoding for 36–37 genes: 12–13 protein-coding genes (PCGs), two ribosomal RNAs (rRNA) genes and 22 transfer RNAs (tRNA) [1, 2]. Mt genomes have been extensively used as genetic markers in molecular phylogenetic studies due to it having several useful properties (i.e., haploidy, compactness, maternal inheritance, relatively high mutation rates, and the lack of recombination) . To date, there are over 4000 complete mt genome sequences of metazoans available in GenBank, including some from parasites [4–10]. In spite of the availability of advanced DNA technologies, however, there is still major gap in our knowledge of mt genomes of the subclass Acari.
Currently, the Acari encompasses over 55,000 species (ticks and mites) and forms the largest group among the class Arachnida . Many species of ticks and mites threaten the health of plants, animals and humans on a global scale, causing major economic losses and significant public health problems [12, 13]. Despite their abundance, medical, veterinary and economic significance, the genetics, epidemiology, biology and phylogenetic relationships in this group remain poorly understood. Acari were traditionally classified into three monophyletic lineages: superorders Parasitiformes (=Anactinotrichida), Opiloacariformes and Acariformes (=Actinotrichida) . However, phylogenetic relationships among Acari have been revised in the past using nuclear rDNA gene sequences, indicating Opiloacariformes + Parasitiformes as the second monophyletic lineage of Acari [15, 16]. Recently, mt genome sequences have been used to understand the phylogenetic relationships among Acari. Dermauw et al. inferred the Acari phylogeny with the mt genome sequences of 21 species. Liu et al. and Burger et al. inferred the phylogeny of ticks with the mt genome sequences of 16 and 21 species, respectively. A difference among mt genome phylogenies is on the relationship within the order Ixodida. Dermauw et al.  and Burger et al.  showed strong support that Metastriate + Prostriate are sister to Ornithodorinae. However, Liu et al. showed strong support that Ornithodorinae is more closely related to Prostriate, and Metastriate is sister to Ornithodorinae + Prostriate. More recently, Burger et al. showed that Metastriate + Prostriate are sister to Ornithodorinae based on nuclear rDNA and mt gene sequences. The systematics and classification of the Acari are controversial in these analyses due to limited or underrepresented taxon sampling. So, the phylogenetic relationships inferred from nuclear and mt sequences need to be retested with more taxa from a wide range of lineages.
In the present study, we sequenced the complete mt genome of the scab mite P. cuniculi in the family Psoroptidae of the order Astigmata. P. cuniculi is a common worldwide parasite of rabbits, causing considerable economic loss due to weight loss, decreased food consumption, vestibular dysfunction, meningitis and death . Also, we inferred the phylogenetic relationships with the concatenated mt amino acid sequences of P. cuniculi and 50 other Acari species that have been sequenced to date.
This study was approved by the Animal Ethics Committee of The Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (Permit code. LVRIAEC2013-006). The New Zealand white rabbit, from which P. cuniculi was collected, was handled in accordance with good animal practices required by the Animal Ethics Procedures and Guidelines of the People’s Republic of China.
Parasites and DNA extraction
P. cuniculi samples were collected from a naturally infected New Zealand white rabbit in Sichuan Province, China. The adult mites were isolated from ear cerumens of the rabbit, according to the previous descriptions . The presence of mites was detected by light microscopic examination, and identification of P. cuniculi was conducted by morphological criteria and site of predilection , then fixed in 70% alcohol and stored at −20°C until use. Total genomic DNA was isolated from these specimens using sodium dodecyl sulphate/proteinase K treatment, followed by spin-column purification (WizardSV Genomic DNA Purification System, Promega). In brief, P. cuniculi collected from the same host were ground in a mortar and subsequently suspended in 275 μL lysis buffer (200 μL Nuclei Lysis Solution, 50 μL EDTA, 20 μL proteinase K and 5 μL RNase A Solution). After vortexing, the sample was incubated at 55°C for 16–18 h, and the DNA was extracted by shaking with 250 μL Wizard SV Lysis Buffer. After centrifugation, the DNA was washed, eluted and stored.
PCR amplification and sequencing
Nine pairs of PCR primers were used to amplify overlapping segments of the complete mt genome of P. cuniculi as shown in Table 1. Briefly, the segment of cox 1 was amplified using the primers COIF/COIR . The other segments were amplified using designed primers (Table 1) based on sequences well conserved in the Dermatophagoides farinae and D. pteronyssinus[17, 23]. All PCR reactions (25 μL) were performed in 2.5 mM dNTP mixture, 2.5 μL 10 × Ex Taq buffer (Mg2+ Plus), 10 pmol of each primer, 1 U Ex Taq polymerase (Takara), 1 μL of DNA sample in a thermocycler (Eppendorf, Germany) under the following conditions: 94°C for 5 min (initial denaturation); then followed by 35 cycles of 94°C for 30 s (denaturation), 36 ~ 60°C for 30 s (annealing), and 68°C for 1 ~ 4.5 min (extension) according to the product length; with a final extension step at 68°C for 10 min. Each PCR reaction yielded a single band as detected in a 1% (W/V) agarose gel upon ethidium-bromide staining (not shown). PCR products were subsequently sent to Invitrogen Biotechnology Company (Shanghai, China) for sequencing using a primer-walking strategy.
Sequences were assembled manually and aligned against the complete mt genome sequences of D. farinae available using the computer program MAFFT 7.122  to identify gene boundaries. Each gene was translated into amino acid sequence using the invertebrate mitochondrial genetic code in MEGA 5 , and aligned based on its amino acid sequence using default settings. The translation initiation and termination codons were identified to avoid gene overlap and to optimize the similarity between the gene lengths of closely-related mite mt genomes. For analyzing tRNA genes, the program tRNAscan-SE  and the program ARWEN  were used to detect tRNA and infer their secondary structure. For tRNAscan-SE, the following parameters were changed: Source: = Nematode Mito, Search Mode = tRNAscan > Cove, Genetic Code = Invertebrate Mito, Cove score cutoff = 0.1. EufindtRNA search parameters = relaxed. Putative secondary structures of 17 tRNA genes were identified, the remaining five tRNA genes (tRNA-SerUCN, tRNA-Ile, tRNA-Pro, tRNA-Val and tRNA-Ala) were identified by recognizing potential secondary structures and anticodon sequences by eye, and two rRNA genes were predicted by comparison with those of other mites previously reported [17, 23].
For comparative purposes, amino acid sequences predicted from published mt genomes of 50 Acari species were also included in the present analysis, using the scorpion Centruroides limpidus (GenBank accession number NC_006896) as an outgroup . The 12 amino acid sequences (not ATP8) were single aligned using MAFFT 7.122 and then concatenated, and ambiguously aligned regions were excluded using Gblocks 0.91b (doc) with the default parameters  using the options for a less stringent selection. Phylogenetic analyses were conducted using three methods: Bayesian inference (BI), Maximum likelihood (ML) and Maximum parsimony (MP). The MtArt + I + G + F model of amino acid evolution was selected as the most suitable model of evolution by ProtTest 2.4  based on the Akaike information criterion (AIC). As MtArt model is not implemented in the current version of MrBayes, an alternative model, MtREV, was used in BI and four chains (three heated and one cold) were run simultaneously for the Monte Carlo Markov Chain. Two independent runs for 1,000,000 metropolis-coupled MCMC generations, sampling a tree every 100 generation in MrBayes 3.1.1 ; the first 2,500 trees represented burn-in and the remaining trees were used to calculate Bayesian posterior probabilities (Bpp). The analysis was performed until the potential scale reduction factor approached 1 and the average standard deviation of split frequencies was less than 0.01. ML analysis was performed with PhyML 3.0  using the subtree pruning and regrafting (SPR) method with a BioNJ starting tree, and the MtArt model of amino acid substitution with proportion of invariant sites (I) and gamma distribution (G) parameters estimated from the data with four discretized substitution rate classes, the middle of which was estimated using the median. Bootstrap frequency (Bf) was calculated using 100 bootstrap replicates. MP analysis was conducted using PAUP 4.0 Beta 10 program , with indels treated as missing character states; 1,000 random additional searches were performed using TBR. Bf was calculated using 1,000 bootstrap replicates, and 100 random taxon additions in PAUP. Phylograms were drawn using the program FigTree v.1.4.
General features of the mt genome of P. cuniculi
The complete mt genome of P. cuniculi was 14,247 bp in length (Figure 1). This size is within the range of other mite mt genomes (Table 2). The sequences have been deposited in GenBank under the accession number KJ957822. This circular mt genome contains 13 PCGs (cox 1-3, nad 1-6, nad 4L, cyt b, atp 6 and atp 8), 22 tRNA genes, two rRNA genes and one D-loop region (Table 3). The genes are transcribed in both different directions. Except for four PCGs (nad 1, nad 6, nad 2 and cyt b) and seven tRNA genes (tRNA-Cys, tRNA-Thr, tRNA-SerAGN, tRNA-Gln, tRNA- Ile, tRNA-Glu and tRNA-LeuCUN) encoded on the minority strand (N-strand), all other genes were encoded on the majority strand (J-strand) (Figure 1). The gene arrangement is the same as those of D. farinae (Pyroglyphidae) and Aleuroglyphus ovatus (Acaridae), but distinct from those of Steganacarus magnus (Steganacaridae) and Panonychus citri (Tetranychidae). The nucleotide composition of the complete mt genome of P. cuniculi is A = 29.4%, T = 41.7%, G = 16.8% and C = 12.1%, with a typically high A + T content of 71.1% within the range of values found in D. farinae (71.4%) and D. pteronyssinus (72.6%) (Table 2).
The length of PCGs of P. cuniculi was in the following order: nad 5 > cox 1 > nad 4 > cyt b > nad 1 > nad 2 > cox 3 > cox 2 > atp 6 > nad 6 > nad 3 > nad 4L > atp 8 (Table 3). A total of 3588 amino acids are encoded in the mt genome of P. cuniculi. In this mt genome, 4 genes (cox 1, nad 3, atp 8 and nad 6) use ATA, 6 genes (cox 2, cox 3, atp 6, nad 4, nad 4L and cyt b) use ATG, 2 genes (nad 1, nad 5) use ATT and 1 gene (nad 2) uses TTG as start codon, respectively (Table 3). All genes have complete termination codon except for cyt b, cox 2, nad 3, nad 5, nad 2 and nad 4 genes which use abbreviated stop codon T or TA, 5 genes (cox 1, cox 3, nad 1, nad 4L and atp 6) use TAA and 2 genes (atp 8 and nad 6) use TAG as termination codon, respectively (Table 3). The rrn S of P. cuniculi is located between tRNA-Asn and tRNA-Val, and rrn L is located between tRNA-Val and tRNA-Trp. The length of the rrn S gene is 658 bp and rrn L gene is 1024 bp (Table 3). The A + T contents of the rrn S and rrn L are 72.4% and 72.2%, respectively. As the mt genomes of Acariform mites can contain short minimal tRNA genes , the identification of tRNA genes can sometimes be challenging. A total of 22 tRNA sequences were identified in the P. cuniculi mt genome, ranging from 50–62 bp (Table 3). Almost all tRNA genes are atypical lacking either the T- or D-loop. Their predicted secondary structures (Figure 2) are similar to that of D. farinae and D. pteronyssinus. As in typical arthropod mtDNA, there is only a small D-loop region between tRNA genes. The D-loop of the P. cuniculi mt genome is located between the tRNA-Phe and tRNA-SerAGN (Table 3). The size of D-loop is 382 bp, and the A + T content is 89.8%. This region contains several dinucleotide [AT]5–40 repeats, which are conserved sequences in mt genomes of D. farinae and D. pteronyssinus[17, 23]. In addition, this region holds a short palindromic sequence, TACAT, which is usually a conserved motif that is often found in mt genomes of insects, including D. farinae and D. pteronyssinus[17, 23].
Of the 51 Acari species included in the phylogenetic analyses in this study, 35 species belonged to the superorder Parasitiformes while 16 belonged to the superorder Acariformes. Both superorders Parasitiformes and Acariformes were monophyletic in all of the trees inferred by the BI, ML and MP methods. The monophyly of the superorder Parasitiformes was strongly supported with a Bpp of 1 in BI analyses (Figure 3), a Bf of 81% in ML analyses (Figure 4), and a Bf of 82% in MP analyses (Figure 5). The monophyly of the superorder Acariformes was strongly supported in BI and ML analyses (Bpp = 1, Figure 3; Bf = 100%, Figure 4), and was moderately supported in MP analysis (Bf = 97%, Figure 5). The 35 species of ticks in the superorder Parasitiformes included in this study were from two orders: Ixodida (32 species) and Mesostigmata (3 species). The monophyly of the order Ixodida was strongly supported in BI and ML analyses (Bpp = 1, Figure 3; Bf = 95%, Figure 4), but was weakly supported in MP analysis (Bf = 46%, Figure 5). The Mesostigmata were monophyletic with strong support in BI analysis (Bpp = 1, Figure 3), and was moderately supported in ML analysis (Bf = 81%, Figure 4), and weakly supported in MP analysis (Bf = 56%, Figure 5). Of the 16 species of mites in the superorder Acariformes, 4 species were from the order Astigmata, one species from the order Oribatida and 11 species from the order Prostigmata. The Astigmata was monophyletic with strong support in all of the three phylogenetic analyses (Bpp = 1, Figure 3; Bf = 100%, Figure 4; Bf = 100%, Figure 5). However, the Prostigmata was not monophyletic in all of the three phylogenetic analyses in the present study (Figures 3, 4 and 5). The Prostigmata was paraphyletic with respect to the family Tetranychidae. Four species from the Tetranychidae of the Prostigmata were more closely related to families Steganacaridae (order Oribatida), Acaridae, Psoroptidae and Pyroglyphidae (order Astigmata) than they were to the other two families (Trombiculidae and Unionicolidae) of the order Prostigmata. The close relationship between these species of the Tetranychidae and Steganacaridae was strongly supported in BI analysis (Bpp = 1, Figure 3), and was moderately supported in ML and MP analyses (Bf = 60%, Figure 4; Bf = 74%, Figure 5). In addition to the Prostigmata, four families of the orders Astigmata and Oribatida were also represented in our analyses: Pyroglyphidae (2 species), Psoroptidae (1 species), Acaridae (1 species) and Steganacaridae (1 species). The Pyroglyphidae was monophyletic with strong support in all of the three phylogenetic analyses in this study (Bpp = 1, Figure 3; Bf = 94%, Figure 4; Bf = 91%, Figure 5). Of the 32 species of ticks in the Ixodida, 21 species were from the family Ixodidae (hard ticks), 10 species from the family Argasidae (soft ticks) and one species from the family Nuttalliellidae. The two families Ixodidae and Argasidae were monophyletic with strong support in BI analysis (Bpp = 1, Figure 3), ML analysis (Bf = 100%, Figure 4) and MP analysis (Bf > 99%, Figure 5). The four species of the genus Amblyomma within Ixodidae were included in this study. These results indicated that Amblyomma is paraphyletic with strong support in BI analysis (Figure 3), and is moderately supported in ML and MP analyses (Figures 3 and 4).
For two decades, there has been considerable debate surrounding the systematics of Acari. The subclass Acari were initially recognized as three lineages (Acariformes, Parasitiformes and Opiloacariformes), however, some authors later demonstrated that there are two distinct groups: Acariformes and Parasitiformes [34–37], and Opiloacariformes as an independent lineage were rejected as they represent basal Parasitiformes . Unfortunately, as there are no complete mt genomes of Opilioacariformes available, we were not able to verify this hypothesis. Dabert et al.  confirmed the monophyly of Acariformes based on the nuclear small subunit rRNA gene (18S rDNA) and the mt cytochrome c oxidase subunit I (COX1). Dermauw et al. tested the hypothesis of the Acari phylogeny with the mt genome sequences of 21 species, and showed strong support of the monophylies of the superorders Parasitiformes and Acariformes, represented by 12 species from the Parasitiformes and 9 species from the Acariformes. This study, however, could not reject the monophyly of the order Prostigmata. In the present study, our analyses support the division of the subclass Acari into two superorders Parasitiformes and Acariformes , however, the order Prostigmata was paraphyletic in all of the three phylogenetic analyses in this study. Although morphological characteristics remain an important source of phylogenetic signal in organismal evolution, morphological characters and molecular characters vary among orders and major groups of Acari. Prostigmata phylogeny based on morphological characteristics shows that the order Prostigmata has an apparent monophyly and is sister to orders Endeostigmata + Oribatids . However, molecular data provide some support that the order Prostigmata was paraphyletic in the current study. However, we should be cautious that the phylogenetic signal may be influenced by long-branch attraction artifacts, this phenomenon may challenge the analyses of phylogeny of Acariform mites .
The superorder Parasitiformes is composed of four orders: Opilioacarida, Holothyrida, Ixodida, and Mesostigmata, and the superorder Acariformes is divided into two suborders: Sarcoptiformes (consisting of Endeostigmata, Oribatida and Astigmata) and Trombidiformes (Prostigmata or Actinedida) . Domes et al. (2007) showed that the order Astigmata is located on the outside of the order Oribatida, usually as sister-group of the order Endeostigmata . However, Dabert et al. speculated that order Oribatida is paraphyletic with respect to order Astigmata. In the present study, although our results have shown that the orders Oribatida and Astigmata are sister groups, only one species within the order Oribatida was included in the present study.
In the present study, our results provide strong support that families Ixodidae and Argasidae are sister taxa, consistent with those of previous studies [4, 42]. Here, we will not discuss the phylogenetic relationships of other families within the superorder Parasitiformes (eg., Nuttalliellidae, Phytoseiidae, Varroidae and Ologamasidae) due to limited data. Phylogenetic analyses based on mt genome sequences in the present study strongly support a close relationship between families Tetranychidae and Steganacaridae + Acaridae + Psoroptidae + Pyroglyphidae to the exclusion of families Trombiculidae and Unionicolidae. Apparently, additional markers for phylogenetic inference are required to resolve the controversy on the order Prostigmata between phylogeny based on mt genome sequences and nuclear SSU rRNA sequences . The monophyly of the order Oribatida will also need further testing with expanded taxon sampling as only one species, Steganacarus magnus, from the Steganacaridae was included in the present study. No species from the superorder Opilioacariformes and order Holothyrida was included in our analyses. Therefore, expanding taxon sampling from these lineages of Acari is clearly the next step for phylogenetic studies of Acari using mtDNA. In addition, phylogenetic analysis in the present study was based only on mtDNA sequences, so we believe it is still necessary to employ nuclear genomic sequences to provide additional evidence for phylogenetic analyses and genome evolution of the Acari.
The present study determined the mt genome of the mite species P. cuniculi of animal health significance, and represents the first mt genome study of any member of the family Psoroptidae. Phylogenetic analyses of the mt genome sequences of the P. cuniculi species, in combination with 50 other Acari species, support the monophylies of the superorders Parasitiformes and Acariformes; the orders Ixodida, Mesostigmata and Astigmata, however, reject the monophyly of the order Prostigmata. Our results indicate the need to test hypotheses of Acari phylogeny with more taxa and different molecular markers.
Wolstenholme DR: Animal mitochondrial DNA, structure and evolution. Int Rev Cytol. 1992, 141: 173-216.
Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27: 1767-1780. 10.1093/nar/27.8.1767.
Tao M, You CP, Zhao RR, Liu SJ, Zhang ZH, Zhang C, Liu Y: Animal mitochondria: evolution, function, and disease. Curr Mol Med. 2014, 14: 115-124. 10.2174/15665240113136660081.
Burger TD, Shao R, Beati L, Miller H, Barker SC: Phylogenetic analysis of ticks (Acari: Ixodida) using mitochondrial genomes and nuclear rRNA genes indicates that the genus Amblyomma is polyphyletic. Mol Phylogenet Evol. 2012, 64: 45-55. 10.1016/j.ympev.2012.03.004.
Jabbar A, Mohandas N, Jex AR, Gasser RB: The mitochondrial genome of Protostrongylus rufescens-implications for population and systematic studies. Parasit Vectors. 2013, 6: 263-271. 10.1186/1756-3305-6-263.
Cameron SL: Insect mitochondrial genomics: implications for evolution for evolution and phylogeny. Annu Rev Entomol. 2014, 59: 95-117. 10.1146/annurev-ento-011613-162007.
Liu GH, Gasser RB, Su A, Nejsum P, Peng L, Lin RQ, Li MW, Xu MJ, Zhu XQ: Clear genetic distinctiveness between human- and pig-derived Trichuris based on analyses of mitochondrial datasets. PLoS Negl Trop Dis. 2012, 6: e1539-10.1371/journal.pntd.0001539.
Liu GH, Chen F, Chen YZ, Song HQ, Lin RQ, Zhou DH, Zhu XQ: Complete mitochondrial genome sequence data provides genetic evidence that the brown dog tick Rhipicephalus sanguineus (Acari: Ixodidae) represents a species complex. Int J Biol Sci. 2013, 9: 361-369. 10.7150/ijbs.6081.
Liu GH, Shao R, Li JY, Zhou DH, Li H, Zhu XQ: The complete mitochondrial genomes of three parasitic nematodes of birds: a unique gene order and insights into nematode phylogeny. BMC Genomics. 2013, 14: 414-426. 10.1186/1471-2164-14-414.
Liu GH, Gasser RB, Young ND, Song HQ, Ai L, Zhu XQ: Complete mitochondrial genomes of the ‘intermediate form’ of Fasciola and Fasciola gigantica, and their comparison with F. hepatica. Parasit Vectors. 2014, 7: 150-159. 10.1186/1756-3305-7-150.
Walter DE, Proctor HC: Systematic and Morphological Survey. Mites: Ecology, Evolution & Behaviour. Edited by: Walter DE, Proctor HC. 2013, New York: Springer, 39-68.
McCoy KD, Léger E, Dietrich M: Host specialization in ticks and transmission of tick-borne diseases: a review. Front Cell Infect Microbiol. 2013, 3: 57-68.
Javed S, Khan F, Ramirez-Fort M, Tyring SK: Bites and mites: prevention and protection of vector-borne disease. Curr Opin Pediatr. 2013, 25: 488-491. 10.1097/MOP.0b013e328362c4ab.
Dunlop JA, Alberti G: The affinities of mites and ticks: a review. J Zool Syst Evol Res. 2008, 46: 1-18.
Murrell A, Dobson SJ, Walter DE, Campbell NJH, Shao R, 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. Invert Syst. 2005, 19: 383-389. 10.1071/IS05027.
Klompen H, Lekveishvili M, Black WC: Phylogeny of parasitiform mites (Acari) based on rRNA. Mol Phylogenet Evol. 2007, 43: 936-951. 10.1016/j.ympev.2006.10.024.
Dermauw W, Van Leeuwen T, Vanholme B, Tirry L: The complete mitochondrial genome of the house dust mite Dermatophagoides pteronyssinus (Trouessart): a novel gene arrangement among arthropods. BMC Genomics. 2009, 10: 107-127. 10.1186/1471-2164-10-107.
Burger TD, Shao R, Labruna MB, Barker SC: Molecular phylogeny of soft ticks (Ixodida: Argasidae) inferred from mitochondrial genome and nuclear rRNA sequences. Ticks Tick Borne Dis. 2014, 5: 195-207. 10.1016/j.ttbdis.2013.10.009.
Cutler SL: Ectopic Psoroptes cuniculi infestation in a pet rabbit. J Small Anim Pract. 1998, 39: 86-87. 10.1111/j.1748-5827.1998.tb03599.x.
Zheng WP, Zhang RH, Wu XH, Ren Y, Nong X, Gu XB, Wang SH, Peng XR, Yang GY: Evaluating troponin C from Psoroptes cuniculi as a diagnostic antigen for a dot-ELISA assay to diagnose mite infestations in rabbits. Parasite Immunol. 2014, 36: 53-59. 10.1111/pim.12076.
Sweatman GK: On the life history and validity of the species in Psoroptes, a genus of mange mites. Can J Zool. 1958, 36: 905-929. 10.1139/z58-078.
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994, 3: 294-299.
Klimov PB, OConnor BM: Improved tRNA prediction in the American house dust mite reveals widespread occurrence of extremely short minimal tRNAs in acariform mites. BMC Genomics. 2009, 10: 598-616. 10.1186/1471-2164-10-598.
Katoh K, Standley DM: MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013, 30: 772-780. 10.1093/molbev/mst010.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA 5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
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.0955.
Laslett D, Canback B: Arwen: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences. Bioinformatics. 2008, 24: 172-175. 10.1093/bioinformatics/btm573.
Dávila S, Piñero D, Bustos P, Cevallos MA, Dávila 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.
Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007, 56: 564-577. 10.1080/10635150701472164.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony, Version 4.0b10. 2002, Sunderland: Sinauer Associates
Alberti G: Zur Feinstruktur der Spermien und Spermiocytogenese der Milben (Acari) I. Anactinotrichida. Zool Jb Anat. 1980, 104: 77-138.
Alberti G: Zur Feinstruktur der Spermien und Spermiocytogenese der Milben (Acari) II. Actinotrichida. Zool Jb Anat. 1980, 104: 144-203.
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. Acarology. Volume 1. Edited by: Griffiths DA, Bowman CE. 1984, Chichester: Ellis Horwood Ltd, 28-62.
Bernini F: Current ideas on the phylogeny and the adaptive radiations of Acarida. Boll Zool. 1986, 53: 279-313. 10.1080/11250008609355512.
Dabert M, Witalinskib W, Kazmierskic A, Olszanowskid Z, Dabertc J: Molecular phylogeny of acariform mites (Acari, Arachnida): Strong conflict between phylogenetic signal and long-branch attraction artifacts. Mol Phylogenet Evol. 2010, 56: 222-241. 10.1016/j.ympev.2009.12.020.
Jeffrey WS: A phylogenetic analysis of the arachnid orders based on morphological characters. Zool J Linn Soc. 2007, 150: 221-265. 10.1111/j.1096-3642.2007.00284.x.
Lindquist EE, Krantz GW, Walter DE: Classification. A Manual of Acarology. Edited by: Krantz GE, Walter DE. 2009, Lubbock: Texas Tech University Press, 97-103. 3
Domes K, Althammer M, Norton RA, Scheu S, Maraun M: The phylogenetic relationship between Astigmata and Oribatida (Acari) as indicated by molecular markers. Exp Appl Acarol. 2007, 42: 159-171. 10.1007/s10493-007-9088-8.
Mans BJ, Klerk D, Pienaar R, de Castro MH, Latif AA: The Mitochondrial genomes of Nuttalliella namaqua (Ixodoidea: Nuttalliellidae) and Argas africolumbae (Ixodoidae: Argasidae): estimation of divergence dates for the major tick lineages and reconstruction of ancestral blood-feeding characters. PLoS One. 2012, 7: e49461-10.1371/journal.pone.0049461.
This work was supported by the International Science & Technology Cooperation Program of China (Grant No. 2013DFA31840), the Science Fund for Creative Research Groups of Gansu Province (Grant No. 1210RJIA006) and Sichuan Provincial Department of Education (Grant No. 14ZB0003).
The authors declare that they have no competing interests.
XBG, GHL, HQS, TYL and GYY performed the experiments, analyzed the data and drafted parts of the manuscript. XQZ conceived and funded the study, revised and edited the manuscript. All authors read and approved the final manuscript.