Exploration of Multi-Gene DNA Barcode Markers to Reveal the Broad Genetic Diversity of Field Ticks (Acari: Ixodidae) in a Tropical Environment of Hainan Island, China.

Ticks are hematophagous arthropods and obligate ectoparasites of humans and other animals. This study focused on the molecular discrimination of ticks in the tropical environment of Hainan according to multi-gene DNA barcode markers with the expectation of accurately distinguishing species. A total of 420 ticks, including 49 adult ticks, 203 nymphal ticks, and 168 larval ticks, were collected in the field, and the 49 adult ticks were identified as Rhipicephalus turanicus, Dermacentor marginatus, and Haemaphysalis longicornis. The mitochondrial 16S rRNA, ribosomal 28S rRNA D2, and ribosomal internal transcribed spacer 2 (ITS2) regions were used as DNA barcode markers to discriminate species. According to basic local alignment search tool analysis against the GenBank database, 16S rRNA positively identified ticks in the Rhipicephalus, Dermacentor, and Haemaphysalis genera; the 28S rRNA D2 region identified ticks in the Rhipicephalus and Dermacentor genera; and ITS2 identified ticks as D. marginatus. Pairwise sequence comparisons based on these three regions were visualized with a Sequence Demarcation Tool matrix. Substitution saturation tests using data analysis and molecular biology and evolution revealed little substitution saturation (Iss < Iss.c, p < 0.05) in the 16S rRNA region for the Haemaphysalis genus; 28S rRNA D2 region for the Rhipicephalus, Dermacentor, and Haemaphysalis genera; and ITS2 region for the Rhipicephalus and Dermacentor genera. Distinctive sequences for which it is difficult to obtain good matches with the sequences available in GenBank exist in the ticks of Hainan. Future studies should obtain complementary sequences to refine and update the database for the molecular characterization of ticks.

Establishment of purification method for prokaryotic expression of Serpin gene for Dermatophagoides farinae.

This study aimed to develop an effective method for the expression and purification of the Dermatophagoides farinae serpin protein and to establish an experimental foundation for elucidating its role in the temperature stress response. The total RNA of D. farinae was extracted, and specific primers were designed for serpin amplification. Serpin was joined with pET32a vector and transformed into BL21 (DE3) cells. Expression of recombinant proteins was induced. Proteins were extracted by enzymatic lysis or enzymatic lysis combined with ultrasonication. Recombinant proteins were purified by Ni-NTA method. SDS-PAGE was conducted to evaluate protein expression, extraction, and purification efficiency. Agarose gel electrophoresis and sequencing analysis showed that the amplified serpin open reading frame was 1284 bp, encoding a hydrophilic and stable protein with a relative molecular weight of 48.30 kD. SDS-PAGE demonstrated that there was a specific band at 55-70 kD, which was consistent with the predicted size of the recombinant pET32a-Serpin protein. Enzymatic lysis combined with 30% ultrasonic power promoted the release of soluble protein more effectively than enzymatic lysis alone. 16 °C for 4 h was optimal for inducing expression. The optimal imidazole concentrations for washing non-His-tagged protein and eluting His-tagged protein were determined to be 20 mM and 200 mM, respectively. In this study, A prokaryotic expression and purification system for the D. farinae serpin protein was successfully established, providing a technical reference for functional gene research in mites at the protein level.

De novo transcriptome sequencing and functional annotation of Demodex canis.

Demodex canis is an important parasitic mite causing skin lesions in dogs. However, molecular studies on this species are limited because of the lack of transcriptomic data. To obtain functional genes of D. canis, mites were collected and RNA was extracted for transcriptome sequencing and functional annotation. Coding sequences (CDSs) of important functional genes were screened and identified using bioinformatics strategies. Six complete CDSs were amplified by specific primers, cloned and sequenced. Results demonstrated that the quantity of the RNA extracted from 100 mites was 12.8 ng and the RNA integrity number was 5.4, indicating that it could be used to construct a cDNA library. Transcriptome sequencing yielded 263,619 unigenes, of which 151,048 (57.3%) were functionally annotated. Using bioinformatics strategies, 58 total CDSs were identified as homologous to those of closely related mite species, including 16 types of allergen genes, 13 types of protease genes, and nine types of motion-related genes. The verified CDSs were almost identical to the CDSs of unigenes. Phylogenetic analyses of tropomyosin further revealed that the verified CDSs clustered successively with those of Demodex species and Tetranychus species in Raphignathae, which agreed with the morphological classification, demonstrating the reliability of the transcriptomic data. In conclusion, this study is the first to successfully sequence transcriptome of D. canis, perform functional annotation, and verify CDSs, which will provide ample data for further studies on the functional genes and molecular pathogenic mechanism of D. canis.

Molecular identification and DNA barcode screening of acaroid mites in ground flour dust.

Molecular identification of acaroid mites is difficult because of the scarcity of molecular data in GenBank. Here, acaroid mites collected from ground flour dust in Xi'an, China, were preliminarily morphologically classified/grouped. Universal primers were then designed to amplify and screen suitable DNA barcodes for identifying these mites. Sixty mite samples were morphologically classified into six groups. Groups 1-2 were identified to Dermatophagoides farinae and Tyrophagus putrescentiae, while Groups 3-6 were not identified to the species level. ITS2 exhibited higher efficiency in molecular identification in comparison with COI, 12S, and 16S. Groups 1-6 were identified as D. farinae, T. putrescentiae, Suidasia nesbitti, Chortoglyphus arcuatus, Lepidoglyphus destructor, and Gohieria sp., respectively. The phylogenetic results were consistent with the morphological classification. Group 6 was further identified as G. fusca according to the morphology of the reproductive foramen. We conclude that the use of ITS2 and the availability of universal primers provide an ideal DNA barcode for molecular identification of acaroid mites. The use of multiple target genetic markers in conjunction with morphological approaches will improve the accuracy of Acaridida identification.

Temperature stress response: A novel important function of Dermatophagoides farinae allergens.

Dermatophagoides farinae, an important pathogen, has multiple allergens. However, their expression under physiological conditions are not understood. Our previous RNA-seq showed that allergens of D. farinae were up-regulated under temperature stress, implying that they may be involved in stress response. Here, we performed a comprehensive study. qRT-PCR detection indicated that 26 of the 34 allergens showed differential expression. Der f1 had the most abundant basic expression quantity. Der f 28.0201 (HSP70) and Der f3 had the same regulation pattern in 9 highly expressed transcripts, which only up-regulated at 41 °C and 43 °C, but Der f 28.0201 showed stronger regulation than Der f 3 (19.88-fold vs 6.02-fold). Whereas Der f 1, 2, 7, 21, 22, 27, and 30 were up-regulated under both heat and cold stress, and Der f 27 showed the strongest regulation ability among them. Der f 27 showed more significant up-regulation than Der f 28.0201 under heat stress (23.59-fold vs 19.88-fold), and Der f27 had more obvious up-regulation under cold than heat stress (30.70-fold vs 23.59-fold). The expression of Der f 27, 28.0201 and 1, and D. farinae survival rates significantly decreased following RNAi, indicating the upregulation of these allergens under temperature stress conferred thermo-tolerance or cold-tolerance to D. farinae. In this study, we described for the first time that these allergens have temperature-stress response functions. This new scientific discovery has important clinical value for revealing the more frequent and serious allergic diseases caused by D. farinae during the change of seasons.

LSU rDNA D5 region: the DNA barcode for molecular classification and identification of Demodex.

Whether ribosomal genes can be used as DNA barcodes for molecular identification of Demodex (Acariformes: Demodicidae) is unclear. To examine this, Demodex folliculorum, D. brevis, D. canis, and D. caprae were collected for DNA extraction, rDNA fragments amplification, sequencing, and analysis. The V2 and V4 regions of SSU rDNA; D5, D6, and D8 regions of LSU rDNA; and ITS region were obtained from the four morphospecies. BLAST analysis showed that the obtained sequences matched those of Demodex or Aplonobia (Acariformes: Tetranychidae) in Raphignathae. Phylogenetic trees derived from V2, V4, D5, D6, and D8 regions, but not from ITS region, showed that the four species of Demodex clustered independently. Sequence divergence analysis further demonstrated that D5, D6, and D8 regions had obvious barcoding gap between intraspecific and interspecific divergences, with the gap of D5 (16.91%) larger than that of D6 (11.82%) and D8 (4.66%). The V2 and V4 regions did not have a barcoding gap, as the intraspecific and interspecific divergences partially overlapped. For the ITS region, intraspecific and interspecific divergences completely overlapped. These results suggest that the D5, D6, and D8 regions of LSU rDNA, especially D5, are suitable DNA barcodes for Demodex.

Establishing an RNA extraction method from a small number of Demodex mites for transcriptome sequencing.

Demodex is a type of parasitic mite which could cause serious dermatoses in 11 orders of mammals. However, due to the tiny body with thick chitin hard to be ruptured as well as the difficulty in obtaining a large number of mites, the quantity and quality of extracted RNA could hardly satisfied for transcriptome sequencing. This has hampered the research on functional genes and molecular pathogenesis of Demodex for a long time. To solve the problems above, the present study established a new RNA extraction method in combination Azanno method with liquid nitrogen grinding using 16 human and canine Demodex mite samples. The RNA quality detection results of Agilent 2100 Bioanalyzer showed that 8 of 16 RNA samples met the requirements for trace RNA-Seq, with RIN of 5.0-6.5 and RNA quantity of 1.1-16.0 ng. RNA quality was affected by grinding process and parasitic position of Demodex. Enough grinding number (≥2000) in moderate time (≤20 min) was significant for mites' complete rupture and RNA degradation prevention. D. brevis (100%, 3/3) parasitizing in human sebaceous glands had significantly higher RNA qualification rate than D. folliculorum (57.14%, 4/7) parasitizing in human hair follicles. Yet D. canis parasitizing in dog had lower RNA qualification rate (16.67%, 1/6) as mites were embedded in skin tissues and blood clots. It should be pointed out that microplate reader had defects with a lower RNA qualification rate of 6.25% (1/16) unmatched with 2100 Bioanalyzer, reminding that it could be only used as reference in RNA quality evaluation.

Screening of reference genes and quantitative real-time PCR detection and verification in Dermatophagoides farinae under temperature stress.

Dermatophagoides farinae is an important source of indoor allergens that shows strong tolerance to external temperatures. However, the regularity and mechanism of tolerance are still unclear. Based on our previous RNA-seq and annotation of D. farinae under temperature stress, it is planned to identify differentially expressed genes (DEGs) involved in the temperature stress response by quantitative real-time PCR (qRT-PCR). However, the lack of reference genes directly limited the detection and confirmation of DEGs. Accordingly, in this study, we have selected six candidates as reference genes in D. farinae: 60S RP L11, 60S RP L21, α tubulin, GAPDH, Der f Mal f 6, and calreticulin, and evaluated their expression stabilities as affected by heat and cold stresses, using geNorm, NormFinder, BestKeeper, comparative ΔCt and RefFinder methods. Then the expression level of 15 DEGs were detected and verified. geNorm analysis showed that α tubulin and calreticulin were the most stable reference genes under heat stress and cold stress of D. farinae. Similar evaluation results were obtained by NormFinder and BestKeeper, in which 60S RP L21 and α tubulin were the most stable reference genes. By comparative ΔCt method and a comprehensive evaluation of RefFinder, α tubulin was identified as the most ideal reference gene of D. farinae under heat and cold stresses. Furthermore, qRT-PCR detection results of 15 DEGs were almost identical to the RNA-seq results, indicating that α tubulin is stable as a reference gene. This study provided technical support for DEGs expression studies in D. farinae using qRT-PCR.

De novo transcriptome sequencing and differential gene expression analysis of two parasitic human Demodex species.

Demodex are among the tiniest organisms in Acari and are important mammalian parasites. However, differences in pathogenicity between two human parasites, Demodex folliculorum and Demodex brevis, remain unknown. Related genetic studies are limited by RNA extraction difficulties and molecular data deficiencies. In this study, RNA extraction, de novo sequencing, functional annotation, and differential gene expression analyses were performed to compare D. folliculorum and D. brevis. This yielded 67.09 and 65.10 million clean reads, respectively, with similar annotations. Bioinformatics analyses and manual alignments identified 237 coding sequences comprising 48 genes from 29 families, including five important functional classes. Of these, 30 genes from 20 families related to metabolism, motion, detoxification and stress response, and allergic reaction were differentially expressed between the two species. Cathepsin type 1, serine protease inhibitor, arginine kinase, triosephosphate isomerase, muscle-specific protein 20-2, myosin alkaline light chain, troponin C, tropomyosin, and heat shock protein 90 were highly expressed in D. folliculorum, whereas cathepsin type 2, aspartic protease, serine protease, myosin heavy chain type 2, and alpha tubulin type 1C were highly expressed in D. brevis. Verified coding sequences were nearly consistent with unigene clusters. Further, absolute quantification results demonstrated that differentially expressed genes followed the predicted expression trend. Therefore, the first RNA sequencing and functional annotation analysis of two Demodex species was successful. Differential expression of important functional genes is likely implicated in pathogenicity disparities between these two species. Our study provides molecular data and technical support for further studies on human Demodex pathogenicity and functional genes.

De novo RNA-seq and functional annotation of Ornithonyssus bacoti.

Ornithonyssus bacoti (Hirst) (Acari: Macronyssidae) is a vector and reservoir of pathogens causing serious infectious diseases, such as epidemic hemorrhagic fever, endemic typhus, tularemia, and leptospirosis. Its genome and transcriptome data are lacking in public databases. In this study, total RNA was extracted from live O. bacoti to conduct RNA-seq, functional annotation, coding domain sequence (CDS) prediction and simple sequence repeats (SSRs) detection. The results showed that 65.8 million clean reads were generated and assembled into 72,185 unigenes, of which 49.4% were annotated by seven functional databases. 23,121 unigenes were annotated and assigned to 457 species by non-redundant protein sequence database. The BLAST top-two hit species were Metaseiulus occidentalis and Ixodes scapularis. The procedure detected 12,426 SSRs, of which tri- and di-nucleotides were the most abundant types and the representative motifs were AAT/ATT and AC/GT. 26,936 CDS were predicted with a mean length of 711 bp. 87 unigenes of 30 functional genes, which are usually involved in stress responses, drug resistance, movement, metabolism and allergy, were further identified by bioinformatics methods. The unigenes putatively encoding cytochrome P450 proteins were further analyzed phylogenetically. In conclusion, this study completed the RNA-seq and functional annotation of O. bacoti successfully, which provides reliable molecular data for its future studies of gene function and molecular markers.

Association study of Demodex bacteria and facial dermatoses based on DGGE technique.

The role of bacteria is unclear in the facial skin lesions caused by Demodex. To shed some light on this issue, we conducted a case-control study comparing cases with facial dermatoses with controls with healthy skin using denaturing gradient gel electrophoresis (DGGE) technique. The bacterial diversity, composition, and principal component were analyzed for Demodex bacteria and the matched facial skin bacteria. The result of mite examination showed that all 33 cases were infected with Demodex folliculorum (D. f), whereas 16 out of the 30 controls were infected with D. f, and the remaining 14 controls were infected with Demodex brevis (D. b). The diversity analysis showed that only evenness index presented statistical difference between mite bacteria and matched skin bacteria in the cases. The composition analysis showed that the DGGE bands of cases and controls were assigned to 12 taxa of 4 phyla, including Proteobacteria (39.37-52.78%), Firmicutes (2.7-26.77%), Actinobacteria (0-5.71%), and Bacteroidetes (0-2.08%). In cases, the proportion of Staphylococcus in Firmicutes was significantly higher than that in D. f controls and D. b controls, while the proportion of Sphingomonas in Proteobacteria was significantly lower than that in D. f controls. The between-group analysis (BGA) showed that all the banding patterns clustered into three groups, namely, D. f cases, D. f controls, and D. b controls. Our study suggests that the bacteria in Demodex should come from the matched facial skin bacteria. Proteobacteria and Firmicutes are the two main taxa. The increase of Staphylococcus and decrease of Sphingomonas might be associated with the development of facial dermatoses.

DNA barcoding for molecular identification of Demodex based on mitochondrial genes.

There has been no widely accepted DNA barcode for species identification of Demodex. In this study, we attempted to solve this issue. First, mitochondrial cox1-5' and 12S gene fragments of Demodex folloculorum, D. brevis, D. canis, and D. caprae were amplified, cloned, and sequenced for the first time; intra/interspecific divergences were computed and phylogenetic trees were reconstructed. Then, divergence frequency distribution plots of those two gene fragments were drawn together with mtDNA cox1-middle region and 16S obtained in previous studies. Finally, their identification efficiency was evaluated by comparing barcoding gap. Results indicated that 12S had the higher identification efficiency. Specifically, for cox1-5' region of the four Demodex species, intraspecific divergences were less than 2.0%, and interspecific divergences were 21.1-31.0%; for 12S, intraspecific divergences were less than 1.4%, and interspecific divergences were 20.8-26.9%. The phylogenetic trees demonstrated that the four Demodex species clustered separately, and divergence frequency distribution plot showed that the largest intraspecific divergence of 12S (1.4%) was less than cox1-5' region (2.0%), cox1-middle region (3.1%), and 16S (2.8%). The barcoding gap of 12S was 19.4%, larger than cox1-5' region (19.1%), cox1-middle region (11.3%), and 16S (13.0%); the interspecific divergence span of 12S was 6.2%, smaller than cox1-5' region (10.0%), cox1-middle region (14.1%), and 16S (11.4%). Moreover, 12S has a moderate length (517 bp) for sequencing at once. Therefore, we proposed mtDNA 12S was more suitable than cox1 and 16S to be a DNA barcode for classification and identification of Demodex at lower category level.

De novo RNA-Seq and functional annotation of Sarcoptes scabiei canis.

The transcriptomic data of Sarcoptes is still lacking in the public database due to the difficulty in extracting high-quality RNA from tiny mites with thick chitin. In this study, total RNA was extracted from live Sarcoptes mites for quality assessment, RNA-Seq, functional annotation, and coding region (CD) prediction and verification. The results showed that the sample JMQ-lngm was qualified for cDNA library construction. Firstly, Agilent 2100 detection showed that the RNA baseline was smooth and the 18S peak was single. Second, the Illumina platform generated 65.78M clean reads and 20,826 unigenes with 35.43M were assembled, occupying 62.98 % of the 56.26M genome. In total, 15,034 unigenes were annotated in seven functional databases. Finally, 13,122 CDs were detected in the 20,826 unigenes, of which 70 complete CDs were matched with Sarcoptes manually in non-redundant nucleotide (NT). Three CDs with indels ≥10 bp were verified. Those results indicated that peritrophin sequences of JMQ-lngm missed 35 bp during the assembly; the pressure-sensitive sodium channel sequences of all the six Sarcoptes scabiei canis isolates were confirmed to be 90 bp shorter than that of a Sarcoptes scabiei hominis isolate; three introns remained in PH chlorine ion channel gating sequences of JMQ-lngm. Moreover, the allergen gene prediction for JMQ-lngm indicated that 61 unigenes were matched with 19 allergen genes of Dermatophagoides, of which Der 1, Der 3, Der 8, and Der 10 had been confirmed in NT. In conclusion, this study successfully completed the RNA-Seq and functional annotation of S. s. canis for the first time, which provides molecular data for future studies on the identification and pathogenic genes of Sarcoptidae.

Improvement on the extraction method of RNA in mites and its quality test.

To solve the long-existing difficult problems in extracting RNA and constructing a complementary DNA (cDNA) library for trace mites, we conducted a further comparative experiment among three RNA extraction methods (TRIzol method, Omega method, and Azanno method) based on our previous attempts at the construction of cDNA library of mites, with Psoroptes cuniculi still used as the experimental subject. By subsequently decreasing the number of mites, the least number of mites needed for RNA extraction of each method were found by criteria of completeness, concentration, and purity of the extracted RNA. Specific primers were designed according to the allergen Pso c1, Pso c2, and Actin gene sequences of Psoroptes to test the reliability of cDNA library. The results showed that Azanno method needed only 10 mites with sensitivity 204 times higher than previously used TRIzol method and 20 times higher than Omega method; clear RNA band was detected by agarose gel electrophoresis; and ultraviolet spectrophotometer determination showed that RNA concentration, 260/280, and 260/230 were in the range of 102 to 166 ng/µl, 1.83 to 1.99, and 1.49 to 1.72, respectively. Finally, specific primers detection showed that the amplified sequences had 98.33, 98.19, and 99.52% identities with those of P. cuniculi or Psoroptes ovis in GenBank, respectively, indicating that the cDNA library constructed using 10 mites was successful and it could meet the requirements for molecular biology research. Therefore, we concluded that Azanno method was more effective than TRIzol method and Omega method in RNA extraction and cDNA library construction of trace mites.

[Human Demodicosis: Emerging Dermatosis Caused by Demodex].

Demodex mites were discovered 175 years ago. However, the pathogenicity of Demodex remains questioned by clinicians due to the high infection yet low incidence of demodicosis and incomplete correlation between mite loads and clinical symptoms. In the recent decades, the pathogenicity of Demodex has attracted public concern as there are continuous case reports of skin lesions caused by mite infection. The pathogenicity of Demodex remains poorly understood. Demodex infection is suspected to be assocaited with rosacea, but its relation with other types of facial lesions such as papule and pustule is usually misdiagnosed or missed in clinical practice. This review gives an overview on the discovery, typical cases, pathogenesis, and clinical diagnosis and treatment of demodicosis, in the aim of promoting social recognition and attention on demodicosis.

Secondary structure of expansion segment D1 in LSU rDNA from Arachnida and its phylogenetic application in Eriophyoid mites and in Acari.

An increasing number of researchers have applied secondary-structure based multiple alignments of rDNA genes in phylogeny. These studies mostly depended on a few valuable divergent domains in LSU and SSU rDNA. Yet other divergent domains, e.g. D1, were poorly investigated and rarely used. However, these domains might contain additional evolutionary data and play a vital role in DNA-based phylogenetic study. Here, we investigated all available D1 sequences of Arachnida taxa and predicted corresponding secondary structures to help identify homologous positions in the D1 region. Long insertions were found exclusive to Eriophyoidea and folded into three newly proposed helices. Non-Acari taxa were all GC rich. In Acari, most Trombidiformes and all Mesostigmata (Parasitiformes) taxa were AT rich and Ixodida (Parasitiformes) GC rich; however there was no consistent base bias in Sarcoptiformes sequences. For Eriophyoid mites, genera Cecidophyopsis and Aceria were both well supported in MP, NJ, ME and ML tress based on D1 sequences, and clusters of Cecidophyopsis species were identical with former study. This demonstrated that the D1 region could act as a valuable molecular marker in phylogenetic reconstruction of Eriophyoidea. Additionally, D1 has been proven suitable in phylogenetic analysis at the family and genus level in Acari, but not in Opiliones.

Constructing and detecting a cDNA library for mites.

RNA extraction and construction of complementary DNA (cDNA) library for mites have been quite challenging due to difficulties in acquiring tiny living mites and breaking their hard chitin. The present study is to explore a better method to construct cDNA library for mites that will lay the foundation on transcriptome and molecular pathogenesis research. We selected Psoroptes cuniculi as an experimental subject and took the following steps to construct and verify cDNA library. First, we combined liquid nitrogen grinding with TRIzol for total RNA extraction. Then, switching mechanism at 5' end of the RNA transcript (SMART) technique was used to construct full-length cDNA library. To evaluate the quality of cDNA library, the library titer and recombination rate were calculated. The reliability of cDNA library was detected by sequencing and analyzing positive clones and genes amplified by specific primers. The results showed that the RNA concentration was 836 ng/µl and the absorbance ratio at 260/280 nm was 1.82. The library titer was 5.31 × 10(5) plaque-forming unit (PFU)/ml and the recombination rate was 98.21%, indicating that the library was of good quality. In the 33 expressed sequence tags (ESTs) of P. cuniculi, two clones of 1656 and 1658 bp were almost identical with only three variable sites detected, which had an identity of 99.63% with that of Psoroptes ovis, indicating that the cDNA library was reliable. Further detection by specific primers demonstrated that the 553-bp Pso c II gene sequences of P. cuniculi had an identity of 98.56% with those of P. ovis, confirming that the cDNA library was not only reliable but also feasible.

Population identification of Sarcoptes hominis and Sarcoptes canis in China using DNA sequences.

There has been no consistent conclusion on whether Sarcoptes mites parasitizing in humans and animals are the same species. To identify Sarcoptes (S.) hominis and S. canis in China, gDNA was extracted from individual mites (five from patients with scabies and five from dogs with mange) for amplification of rDNA ITS2, mtDNA 16S, and cox1 fragment sequences. Then, the sequences obtained were aligned with those from different hosts and geographical locations retrieved from GenBank and sequence analyses were conducted. Phylogenetic trees based on 317-bp mtDNA cox1 showed five distinctive branches (species) of Sarcoptes mites, four for S. hominis (S. hominis Chinese, S. nr. hominis Chinese, S. hominis Australian, and S. hominis Panamanian) and one for S. animal (S. animal). S. animal included mites from nine animal species, with S. canis China, S. canis Australia, and S. canis USA clustering as a subbranch. Further sequence divergence analysis revealed no overlap between intraspecific (≤ 2.6 %) and interspecific (2.6-10.5 %) divergences in 317-bp mtDNA cox1. However, overlap was detected between intra- and interspecific divergences in 311-bp rDNA ITS2 or 275-bp mtDNA 16S when the divergences exceeded 1.0 %, which resulted in failure in identification of Sarcoptes. The results showed that the 317-bp mtDNA cox1 could be used as a DNA barcode for molecular identification of Sarcoptes mites. In addition, geographical isolation was observed between S. hominis Chinese, S. hominis Australian, and S. hominis Panamanian, but not between all S. canis. S. canis and the other S. animal belonged to the same species.

Molecular identification of four phenotypes of human Demodex in China.

Traditional classification of Demodex mites by hosts and phenotypic characteristics is defective because of environmental influences. In this study, we proposed molecular identification of four phenotypes of two human Demodex species based on mitochondrial cox1 fragments for the first time. Mites collected from sufferers' facial skin were classified into four phenotypes: phenotype A-C with finger-like terminus, and phenotype D with cone-like terminus. The results of molecular data showed that cox1 sequences were all 429 bp. Divergences, genetic distances and transition/transversion ratios among the three phenotypes with finger-like terminus were 0.0-3.0%, 0.000-0.031, and 6/3-5/0, respectively, in line with intraspecific differences. However, those measures between the phenotype with cone-like terminus and phenotypes with finger-like terminus were 19.6-20.5%, 0.256-0.271, and 0.58 (31/53)-0.66 (35/53), respectively, reaching interspecific level. Phylogenetic trees also showed that the three phenotypes with finger-like terminus clustered as one clade, and the phenotype with cone-like terminus formed another one. Therefore, we conclude that mitochondrial cox1 sequence is a good marker for identification of two human Demodex species. Molecular data indicate no subspecies differentiation. Terminus is an effective character for species identification. Mites with finger-like terminus are Demodex folliculorum, and those with cone-like terminus are Demodex brevis.

Molecular identification and phylogenetic study of Demodex caprae.

The DNA barcode has been widely used in species identification and phylogenetic analysis since 2003, but there have been no reports in Demodex. In this study, to obtain an appropriate DNA barcode for Demodex, molecular identification of Demodex caprae based on mitochondrial cox1 was conducted. Firstly, individual adults and eggs of D. caprae were obtained for genomic DNA (gDNA) extraction; Secondly, mitochondrial cox1 fragment was amplified, cloned, and sequenced; Thirdly, cox1 fragments of D. caprae were aligned with those of other Demodex retrieved from GenBank; Finally, the intra- and inter-specific divergences were computed and the phylogenetic trees were reconstructed to analyze phylogenetic relationship in Demodex. Results obtained from seven 429-bp fragments of D. caprae showed that sequence identities were above 99.1% among three adults and four eggs. The intraspecific divergences in D. caprae, Demodex folliculorum, Demodex brevis, and Demodex canis were 0.0-0.9, 0.5-0.9, 0.0-0.2, and 0.0-0.5%, respectively, while the interspecific divergences between D. caprae and D. folliculorum, D. canis, and D. brevis were 20.3-20.9, 21.8-23.0, and 25.0-25.3, respectively. The interspecific divergences were 10 times higher than intraspecific ones, indicating considerable barcoding gap. Furthermore, the phylogenetic trees showed that four Demodex species gathered separately, representing independent species; and Demodex folliculorum gathered with canine Demodex, D. caprae, and D. brevis in sequence. In conclusion, the selected 429-bp mitochondrial cox1 gene is an appropriate DNA barcode for molecular classification, identification, and phylogenetic analysis of Demodex. D. caprae is an independent species and D. folliculorum is closer to D. canis than to D. caprae or D. brevis.