|Year : 2020 | Volume
| Issue : 2 | Page : 136-141
Molecular identification and phylogenetic relationship of Demodex mites based on mitochondrial 16S rDNA
Pawan Prasher1, Dolly Baghra2, Drishtant Singh3, Sharad Thakur3, Navpreet Kaur Gill2, Anup Kumar Kesavan3
1 Sri Guru Ram Das Institute of Medical Science and Research, Amritsar, Punjab, India
2 Department of Zoology and Environmental Sciences, Punjabi University, Patiala, Punjab, India
3 Molecular Microbiology Lab, Department of Molecular Biology and Biochemistry, Guru Nanak Dev University, Amritsar, Punjab, India
|Date of Submission||11-Nov-2019|
|Date of Acceptance||06-Jan-2020|
|Date of Web Publication||23-Jan-2021|
Anup Kumar Kesavan
Department of Molecular Biology and Biochemistry, Guru Nanak Dev University, Amritsar, Punjab
| Abstract|| |
Background and Objective: Demodex mites are tiny parasites that live around hair follicles of mammals. The two main species of Demodex i.e. Demodex folliculorum and Demodex brevis present in humans are found near the hair follicles of eyes. The present study was to understand the presence of Demodex mites in people suffering from blepharitis in Amritsar, Punjab.
Material and Methods: Demodex mites samples present in blepharitis patients were isolated from the eyelashes. DNA was isolated from three mites and used for PCR amplification of mitochondrial (mt) 16S rDNA. The amplified PCR product were purified and used for molecular identification.
Results: The amplified mt16s rDNA product was sequenced and subjected to BLAST search in the NCBI database for molecular identification. The identified mite belongs to Demodex folliculorum species. The phylogenetic tree constructed by using mt16s rDNA sequence suggests that D. folliculorum is closer to D. canis than to D. brevis.
Conclusion: All the three isolates belong to D. folliculorum and the mitochondrial DNA 16S rDNA partial sequence is applicable for phylogenetic relationship analysis.
Keywords: Demodex blepharitis, Demodex mites, mitochondrial 16S rDNA, phylogenetic tree
|How to cite this article:|
Prasher P, Baghra D, Singh D, Thakur S, Gill NK, Kesavan AK. Molecular identification and phylogenetic relationship of Demodex mites based on mitochondrial 16S rDNA. Trop Parasitol 2020;10:136-41
|How to cite this URL:|
Prasher P, Baghra D, Singh D, Thakur S, Gill NK, Kesavan AK. Molecular identification and phylogenetic relationship of Demodex mites based on mitochondrial 16S rDNA. Trop Parasitol [serial online] 2020 [cited 2021 Mar 9];10:136-41. Available from: https://www.tropicalparasitology.org/text.asp?2020/10/2/136/307799
| Introduction|| |
Demodex mites belong to the class: Arachnida; order: Prostigmata; family: Demodicidae; and genus: Demodex. Simon first described the genus Demodex in 1842. Since then a total of 140 Demodex species or subspecies have been identified. Demodex mites generally infest hair follicles, sebaceous glands, meibomian glands, ceruminous glands, and internal organs of 11 orders of mammals, including dog, sheep, cat, pig, and mouse. It is generally considered that Demodex is a host-specific obligate parasite. Two or more Demodex species might simultaneously parasitize the same host. Cross infection between humans and animals has also been reported.,
The traditional classification of the Demodex species mainly depends on the morphological characteristics or the phenotype, but it has certain limitations that induce certain difficulties and indeterminacy. Therefore, it is necessary to study the Demodex mites phylogenetic relationship at the molecular level. A research on the Demodex mites has been limited by the three factors. First, less attention has been paid to Demodex because its pathogenicity is not certain. Second, DNA extraction of Demodex is troublesome as mites are tiny with thick chitin which is difficult to rupture. Third, it is difficult to obtain a large number of standard samples of mites because these mites are species specific to their host, and it cannot be maintained or cultured in vitro.,,
The increased number of studies have shown a strong association of Demodex infestation with various dermatitis and ocular problems such as rosacea, blepharitis, and acne.,,,,, Therefore, interest in Demodex has been raised, and studies at molecular level are increasing.,
Various molecular techniques have been applied for the classification and phylogenetic studies of the Demodex mites. Zhou and Cheng, analyzed the polymorphism of Demodex folliculorum and Demodex brevis by the random amplified polymorphic method (RAPD). Toops et al. isolated genomic DNA from Demodex canis and suggested that the Demodex sequences obtained through oligonucleotide primer amplification can be useful in analyzing the phylogenetic relationship and disease progression. Ravera et al. detected D. canis DNA 166 bp chitin synthase gene using the real-time polymerase chain reaction (PCR). Similarly, Zhao and Wu analyzed D. folliculorum, D. brevis, and D. canis by RAPD sequence characterized amplified region marker.
Mitochondrial rDNA, being maternally inherited, with fast evolution and relatively simple and nontissue-specific gene structure, is useful for studying the phylogeny of organisms. In general, cytochrome c oxidase subunit genes and protein-coding genes (CO1, COII, NDI, and NDS) are useful markers for the phylogenetic analysis of closely related species, sub species, and population from different geographical regions, whereas 12S rDNA and 16S rDNA genes are suitable for the phylogenetic analysis at the genera and species level.
In the present study, 16S rDNA partial sequences of three isolates of D. folliculorum have been amplified, sequenced, and then analyzed with 17 D. folliculorum isolates (2 isolates from Spain, 5 isolates from Iran, and 10 isolates from China), 4 isolates of D. brevis from China, and 6 isolates of D. canis (1 isolate from USA, 2 isolates from China, and 3 isolates from Spain) to discuss the phylogenetic relationship in mites. Genetic distance was computed, and phylogenetic tree was reconstructed using MEGA 6.0 (www.megasoftware.net).
| Materials and Methods|| |
Mite collection method
Samples of the D. folliculorum were obtained from the eyelashes of three patients diagnosed with Demodex blepharitis. The diagnosis of blepharitis was established on the slit-lamp examination with patients who presented with the symptoms of ocular discomfort and exhibited characteristic cylindrical dandruff like deposits on the base of eyelashes [Figure 1]. Demodex mites were identified using the method given by Coston. Four eyelashes were epilated from each patient and placed on a glass slide and were observed under the light microscope to detect the presence of mites [Figure 2]. The mites were preserved in the Tris–EDTA solution and subjected for molecular analysis.
|Figure 1: Representative picture of a patient with blepharitis showing characteristic cylindrical dandruff like deposits on the bases of eyelashes|
Click here to view
|Figure 2: Representative picture of Demodex mites seen in relation to the base of hair follicles as seen on light microscopy|
Click here to view
Genomic DNA isolation of mites
Mites were suspended in 1 ml of Tris-EDTA. Samples were centrifuged for 5 min at 3000 rpm, and the DNA was extracted using the method given by Toops et al. The primers used for 16S rDNA amplification were forward 5'CTGTGCTAAGGYAGCGAAGTC3' and reverse 5'TCAAAWGCCAACAKCGAGGTAA3' (GenBank no. JF783994).
Polymerase chain reaction
The following conditions were used for PCR amplification of 16S rDNA sequences. Taq buffer 2 μl, dNTPs 1.8 μl, primer forward 0.4 μl, primer reverse 0.4 μl, MgCl2 0.8 μl, DNA polymerase (Taq) 0.2 μl, MQ water 10.4 μl, and DNA 4 μl to make a final volume of 20 μl. The amplification condition for 16S rDNA was 94°C for 3 min, denaturation at 94°C for 30 s, annealing at 52°C for 30 s, elongation at 72°C for 30 s, and final extension at 72°C for 10 min and followed by 35 cycles. PCR fragments were detected by the 1.5% agarose gel electrophoresis, 1 μl of EtBr is added in the solution, along with the molecular marker (100 bp) ladder at 100 V for 1 h, and then visualized under ultraviolet light in a Gel documentation system (Protein Simple-Alpha Imager Mini system) [Figure 3]. The amplified products were stored at 4°C. PCR products were purified using the Gel/PCR DNA fragments extraction kit (IBI Scientific) according to the manufacturer's instructions. DNA fragments were sequenced in both directions using the same primers used in amplification steps.
|Figure 3: Polymerase chain reaction amplification of 16S rDNA from Demodex species. Lane 1: DNA Ladder 100 bp. Lane 2, 3, 4: polymerase chain reaction products of three Demodex species|
Click here to view
Sequence alignments for 16S rDNA sequences of the Demodex isolates obtained in the present study were aligned with DNA sequences of Demodex species obtained from the GenBank by using MUSCLE. Maximum composite likelihood estimate of the pattern of nucleotide substitution, nucleotide composition, neighborhood joining tree, maximum likelihood tree, variable sites, and parsim-info sites were analyzed using MEGA.
| Results|| |
The mitochondrial 16S rDNA PCR products of the three D. folliculorum species from Amritsar, India, were successfully amplified [Figure 3] and sequenced (Accession number KX151160.1, KX151161.1, KX151162.1). These three DNA sequences were aligned with 17 sequences of D. folliculorum which was retrieved from the GenBank [Table 1] to reconstruct the phylogenetic tree by maximum likely hood method, neighborhood, and nucleotide composition in MEGA 6.0 to confirm the reliability of phylogenetic relationship.
|Table 1: Information of 16S rDNA sequences of 27 species of Demodex species retrieved from GenBank|
Click here to view
The alignment of the three sequences of 16S rDNA partial sequences obtained in the present study and the 17 sequences of the D. folliculorum retrieved from the GenBank using MEGA resulted in a total of 300 base pairs including gaps. The MEGA (version 7) analysis presented a (96%) 288/299 sequence identity with the 15 D. folliculorum isolates. It presented a (98%) 243/249 sequence identity with one sequence of D. folliculorum. The overall number of the conserved site in the three isolates was 210/212, variable sites 2/212, and no parsim-info. The average content of T, C, A, and G in D. follculorum was 43.5%, 7.3%, 27.7%, and 21.6% respectively, and the nucleotide frequencies were biased toward A+T, averaging 71.2% [Table 2].
|Table 2: Nucleotide composition of 16S rDNA gene of three isolates of Demodex folliculorum|
Click here to view
Nucleotide pattern and protein translations
In the present study, many silent mutations were observed, and most of them were at the third codon position. The nucleotide sequence and the proportion (percentage) of a particular nucleotide (A, C, T, and G) in the DNA strand are important in the phylogenetic analysis and describe a specific evolutionary trend. In the present study, as shown in [Table 3], the average nucleotide frequencies of A, C, G, and T in D. folliculorum are 30.70% (A), 41.42% (T/U), 7.26% (C), and 20.61% (G). The transition/transversion rate ratios are k1 = 1.938 (purines) and k2 = 25.043 (pyrimidines). The overall transition/transversion bias is R = 3.507. The analysis involved 30 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 190 positions in the final dataset. Evolutionary analyses were conducted in MEGA 6.0.
|Table 3: Maximum composite likelihood estimate of the pattern of nucleotide substitution|
Click here to view
The two phylogenetic trees (NJ and ML) showed similar topological structures. The three isolates of the three phenotypes with finger-like terminus clustered as one branch with high node support, and they form the sister clade indicating that three phenotypes were of the Demodex species. D. canis was closer to D. folliculorum than to D. brevis. D. folliculorum species isolated from China were gathered with the three D. folliculorum isolates from Amritsar, India, and then clustered with D. canis and with the D. brevis [Figure 4] and [Figure 5].
| Discussion|| |
The current study showed that the three isolates have a conserved mitochondrial 16S rDNA sequence with some variations among the same species. These belong to the D. folliculorum species based on the morphological characteristics and sequencing. In addition, phylogenetic studies revealed that all three isolates clustered with other sequences of D. folliculorum, retrieved from GenBank database, with minor single-nucleotide polymorphisms variations. The three isolates were more in A+T nucleotide frequencies than G+C nucleotide frequencies and the estimated transition/transversions bias (R) is 0.93. These findings are in agreement with other studies that showed molecular analysis using several mitochondrial gene fragments.
De Rojas et al. analyzed the Demodex mitochondrial 16S rDNA and showed that the A+T nucleotide frequencies were higher than the G+C, which is in consistent with Rhinonyssid mites. Herbert et al., Tsao and Yeh, and Liu et al. stated that the average genetic distance among the Demodex species is larger than the intraspecific distance. Kumar et al. and Zhou et al. observed the transition/transversions among the three Demodex species are smaller, which indicate that the base substitutions have reached the interspecies level. Zhao and Wu found that the differentiation between D. folliculorum and D. canis is smaller than that between D. folliculorum and D. brevis. Hence, the mitochondrial 16S rDNA partial sequence can be used to identify the difference among the three Demodex species. In the present study, it has been observed that D. folliculorum is more closely related to D. canis than to D. brevis.
De Rojas et al. found the synthetically analytic results of the divergence, genetic distance, and transition/transversions of the two geographic isolates Spain and China, it can be found that the mitochondrial 16S rDNA partial sequence cannot discriminate the difference between the two geographic D. folliculorum isolates. Moreover, 16S rDNA fragments of three D. folliculorum (face) from Xi'an, China and six D. folliculorum (eyelids) from Spain are completely identical and belong to one haplotype, which suggest that neither mitochondrial 16S rDNA partial sequence cannot discriminate the Demodex mites from Chinese faces and Spanish eyelids nor it can tell the difference between the D. folliculorum from the distinct body localization from Spain.
From the above studies, two explanations may arise: one reason is that the mitochondrial 16S rDNA partial sequences were highly conserved and cannot be used for the interspecies identification and second, D. folliculorum has only one species, and subspecies differentiation happened due to the changes of parasitic position and geographic condition. During present studies, it has been observed that three isolates of Demodex show 99% query cover with the D. folliculorum isolates from China.
Zhao and Wu studied the phylogenetic trees inside the genus Demodex, high bootstrap value of 100% supports the presence of three sister branches, where two species, D. folliculorum and D. canis, clustered together (bootstrap 96%–100%), but the third one, D. brevis, seem to be more distant (bootstrap 98%–100%). The phylogenetic tree proved that at the molecular level, the D. folliculorum, D. canis, and D. brevis are three distinct species, and D. folliculorum is closer to D. canis than to D. brevis, which is in consistence with their morphological classification. This may due to the reason: both D. folliculorum and D. brevis are a human parasite, but at a different location, former in the hair follicles and later in the sebaceous glands, where the environment is different. D. canis mites are present in the hair follicles of the dog where the environment is similar to that of the human hair follicle. It might be a similar environment that results into the similar phenotype and genotype of D. folliculorum and D. canis. In conclusion, all three isolates belong to D. folliculorum, and the mitochondrial DNA 16S rDNA partial sequence is applicable for phylogenetic relationship analysis at the lower level, i.e., genus and species but not for interspecies identification.
| Conclusion|| |
Demodex mites present in the patients of blepharitis from Amritsar belongs to D. folliculorum species. The mitochondrial DNA 16S rDNA parital sequence was applicable for phylogenetic relationship analysis at lower level but not for interspecies identification.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Zhao YE, Wu LP. Phylogenetic relationships in Demodex
mites (Acari: Demodicidae
) based on mitochondrial 16S rDNA partial sequences. Parasitol Res 2012;111:1113-21.
Zhao YE, Ma JX, Hu L, Wu LP, De Rojas M. Discrimination between Demodex folliculorum
) isolates from China and Spain based on mitochondrial cox1 sequences. J Zhejiang Univ Sci B 2013;14:829-36.
Morsy TA, el Okbi MM, el-Said AM, Arafa MA, Sabry AH. Demodex
(follicular mite) infesting a boy and his pet dog. J Egypt Soc Parasitol 1995;25:509-12.
Wang YP, Li P, Bing G. A case report of human dermatitis caused by Canine Demodex
. J N Bethune Univ Med Sci 1998;24:265.
Zhao YE, Guo N, Wu LP. Influence of temperature and medium on viability of Demodex folliculorum
and Demodex brevis
). Exp Appl Acarol 2011;54:421-5.
Zhao YE, Guo N, Xun M, Xu JR, Wang M, Wang DL. Sociodemographic characteristics and risk factor analysis of Demodex
infestation (Acari: Demodicidae
). J Zhejiang Univ Sci B 2011;12:998-1007.
Zhao YE, Wu LP, Hu L, Xu JR. Association of blepharitis with Demodex
: A meta-analysis. Ophthalmic Epidemiol 2012;19:95-102.
Moravvej H, Dehghan-Mangabadi M, Abbasian MR, Meshkat-Razavi G. Association of rosacea with demodicosis. Arch Iran Med 2007;10:199-203.
Bonamigo RR, Bakos L, Cartell A, Edelweiss MI. Factors associated with rosacea in population samples of Southern Brazil: Analysis of case-control studies. Bras Dermatol 2008;83:419-24.
Hsu CK, Hsu MM, Lee JY. Demodicosis: A clinicopathological study. J Am Acad Dermatol 2009;60:453-62.
Zhao YE, Wu LP, Peng Y, Cheng H. Retrospective analysis of the association between Demodex
infestation and rosacea. Arch Dermatol 2010;146:896-902.
Zhao YE, Hu L, Wu LP, Ma JX. A meta-analysis of association between acne vulgaris and Demodex
infestation. J Zhejiang Univ Sci B 2012;13:192-202.
Zhao YE, Wu LP, Hu L, Xu Y, Wang ZH, Liu WY. Sequencing for complete rDNA sequences (18S, ITS1, 5.8S, ITS2, and 28S rDNA) of Demodex
and phylogenetic analysis of Acari based on 18S and 28S rDNA. Parasitol Res 2012;111:2109-14.
Zhao YE, Wu LP. RAPD-SCAR marker and genetic relationship analysis of three Demodex
species (Acari: Demodicidae
). Parasitol Res 2012;110:2395-402.
Zhao YE, Guo N, Wu LP. The effect of temperature on the viability of Demodex folliculorum
and Demodex brevis
. Parasitol Res 2009;105:1623-8.
Toops E, Blagburn B, Lenaghan S, Kennis R, MacDonald J, Dykstra C. Extraction and characterization of DNA from Demodex canis
. J Appl Res Vet Med 2010;8:31-43.
Ravera I, Altet L, Francino O, Bardagí M, Sánchez A, Ferrer L. Development of a real-time PCR to detect Demodex canis
DNA in different tissue samples. Parasitol Res 2011;108:305-8.
Xu QG, Hua BZ. Application of mtDNA in phylogenetic analysis of insects. J Northwest Sci Tech Univ Agric For 2001;29:S79-83.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol 2013;30:2725-9.
Coston TO. Demodex folliculorum
blepharitis. Trans Am Ophthalmol Soc 1967;65:361-92.
Ferreira D, Sastre N, Ravera I, Altet L, Francino O, Bardagí M, et al
. Identification of a third feline Demodex
species through partial sequencing of the 16S rDNA and frequency of Demodex
species in 74 cats using a PCR assay. Vet Dermatol 2015;26:239-e53.
De Rojas M, Mora MD, Ubeda JM, Cutillas C, Navajas M, Guevara DC. Phylogenetic relationships in rhinonyssid mites (Acari: Rhinonyssidae) based on mitochondrial 16S rDNA sequences. Exp Appl Acarol 2001;25:957-67.
Hebert PD, Stoeckle MY, Zemlak TS, Francis CM. Identification of Birds through DNA Barcodes. PLoS Biol 2004;2:e312.
Tsao WC, Yeh WB. DNA-based discrimination of subspecies of swallowtail butterflies (Lepidoptera: Papilioninae) from Taiwan. Zool Stud 2008;47:633-43.
Liu SF, Chen LL, Dai FQ, Zhuang ZM. Application of DNA barcoding gene COI for classificating family sciagenidae. Oceanol Limnol Sin 2010;41:223-31.
Kumar S, Nei M, Dudley J, Tamura K. MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform 2008;9:299-306.
Zhou JY, Zhang Q, Tang YL, Yu FY, Zhao S. On phylogenetic relationships of Teraponidae in coastal waters of China. Mar Fish 2010;32:351-5.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5]
[Table 1], [Table 2], [Table 3]