Study of Genetic Variation of Leishmania major Based on Internal Transcribed Spacer 1 (ITS1) in Chabahar, Iran

authors:

avatar Mansour Dabirzadeh ORCID 1 , * , avatar Mohammad Hashemi 1 , avatar Yahya Maroufi 1

Department of Parasitology and Mycology, Zabol University of Medical Sciences, Zabol, IR Iran

how to cite: Dabirzadeh M, Hashemi M, Maroufi Y. Study of Genetic Variation of Leishmania major Based on Internal Transcribed Spacer 1 (ITS1) in Chabahar, Iran. Jundishapur J Microbiol. 2016;9(6):e60076. https://doi.org/10.5812/jjm.33498.

Abstract

Background:

Zoonotic cutaneous leishmaniasis (ZCL) is polymorphic disease that may show various clinical manifestations.

Objectives:

This study investigates the determination of genetic variation within the species of Leishmania major isolates from new cases in Chabahar, a port city in Southeast Iran (situated at the Iran-Pakistan border). Migration in this region indicates that leishmaniasis is spreading gradually, and a new micro-habitat focus appears each year.

Materials and Methods:

A variety of nucleic acid detection methods that target both DNA and RNA have been developed. The restriction fragment length polymorphism analysis of amplified internal transcribed spacer 1 with polymerase chain reaction (ITS1-RFLP PCR) assay is a multipurpose tool for the diagnosis of Leishmania from clinical samples and for enabling the determination of the infecting Leishmania species. The goal of this study was the identification of species based on ITS1-RFLP in the ribosomal operon of L. major from clinically different forms of ZCL amplified by PCR, followed by the digestion of the PCR product with restriction enzymes. The profiles were observed and visualized in agarose gel under UV light. We used direct smears to identify the parasites. While taking the smear, samples were collected for culture or direct PCR. We used the PCR-RFLP assay of the ITS1 genes for direct identification of Leishmania species in 24 out of 33 suspected patients. PCR-ITS1 amplification was done on the 24 samples confirmed by culture via growth and parasitological methods.

Results:

Of the 24 isolates, 21 had 350 bp bands (87.5%) and three had 450 bp bands (12.5%). After using the restriction enzyme, banding patterns including fragments of 210 and 140 bp for L. major were detected in 19 cases.

Conclusions:

The L. major species causing ZCL in Chabahar have limited genetic variation. There seems to be little manifestation of diversity between these lesions as a new focus of disease, and new micro-habitats for the disease are appearing in parts of this region.

1. Background

Zoonotic cutaneous leishmaniasis (ZCL) is polymorphic disease with various clinical manifestations (1, 2). Leishmaniasis, a group of parasitic infections, has a worldwide distribution (3); the world health organization estimates a prevalence of approximately 12 million cases, with an annual mortality rate of 60,000 and an at-risk population of approximately 350 million (4). Approximately 90% of CL cases occur in just seven countries, including Iran (5). Several molecular targets for diagnostic PCR testing have been evaluated in Leishmania, including minicircle kinetoplast DNA (kDNA) (5, 6), the mini-exon (spliced leader RNA) gene (7), gp63 polymerase chain reaction analysis of restriction fragment length polymorphism (PCR-RFLP) (8), and the internal transcribed spacer (ITS) (1, 3, 9, 10). PCR is a new alternative to existing diagnostic procedures, such as direct smear of parasites from clinical specimens or by cultivation, with microscopic examination.

2. Objectives

In the present study, as described previously (1, 11, 12), we used the RFLP analysis of amplified ITS1 in ribosomal operons to investigate the parasite that causes ZCL, and the genetic variations among L. major isolates from Chabahar, a port city in Southeast Iran (situated at the Iran-Pakistan border). In this region, the manifestation of ZCL as a new focus of disease was in doubt, but our findings showed that the species is L. major with variable genes. The disease was correlated with the clinical manifestations of ZCL in Chabahar, so this is a new major public health problem.

3. Materials and Methods

3.1. Study Area and Population

A variety of nucleic acid detection methods that target both DNA and RNA have been developed (1). The ITS1-RFLP PCR assay is a multipurpose tool for the diagnosis of Leishmania from clinical samples, and it enables determination of the infecting species (3). The procedure was performed as described by Schonian et al. The study was approved by the institutional ethics committee, and all patients signed informed consent forms.

We evaluated the PCR-RFLP assay of the ITS1 genes for direct identification of Leishmania species in 24 out of 33 suspected patients. These patients were referred for diagnosis to the central laboratory of Chabahar (25° 17’ N.70° 38’E; 13,162 km2) for evaluation of skin diseases. The average annual temperature and humidity in Chabahar are 36.4°C and 75.9%, respectively. The population migration into and emigration out of this port have important medical implications.

The clinical samples were taken from patients suspected to have CL in different endemic areas of Chabahar, including Negor and other districts and villages. Each diagnosis was confirmed by the demonstration of amastigotes on the slit smear and/or the presence of flagellated promastigotes in Novy-MacNeal-Nicolle (NNN) medium from skin samples of 24 patients with typical and atypical lesions.

For the parasitological investigation, three skin samples were taken from the edges of lesions using sterile surgical blades. Two were used for direct smears and were stained with giemsa, and the third was inoculated into sterile screw-tap tubes that contained blood agar slanted in NNN medium, then incubated at 25 ± 1°C. The isolated promastigotes were subcultured in RPMI-1640 medium supplemented with 15% fetal bovine serum (FBS), 100 μg/mL streptomycin, and 100 U/mL penicillin. Finally, 24 positive-culture isolates were selected for molecular examination.

3.2. Parasite Samples

The diagnosis of leishmaniasis was confirmed by the presence of amastigotes in slit smears and/or the presence of flagellated promastigotes in NNN medium from typical and atypical lesions. For parasitological investigations, three skin samples were collected from the edges of the lesions using sterile surgical blades. Two of the samples were used for direct smears and were stained with giemsa, and the third was inoculated into sterile screw-top tubes containing blood agar slanted in NNN medium, then incubated at 25 ± 1°C. The isolated promastigotes were subcultured in RPMI-1640 medium supplemented with 15% FBS. They were washed twice by resuspension in buffered saline, then centrifuged. The pellets were stored at -70°C until further processing.

Parasites from a 15 mL mid-logarithmic phase of the bulk culture were harvested by centrifugation (700 × g for 20 min at 4°C) and washed three times in ice-cold sterile PBS (pH 7.2).

3.3. DNA Preparation

For species identification, the DNA isolated from amastigotes in ulcers and promastigotes in cultures was extracted using a commercial extraction kit (High Pure Template DNA Preparation Kit, Roche, Germany), according to the manufacturer’s instructions.

Leishmania reference strains were obtained from Isfahan and Tehran universities. The target ITS1-DNA was amplified in a reaction mixture that consisted of 1 μL of LITSR: 5’-CTTG GATCATTTTCCGATG-3’ and 1 μL of L5.8S 5’-TGA TAC CAC TTA TCG CAT T-3’ (12). The reaction was carried out with the PCR-Ready Supreme mix in 47 μL of total reaction, 15 μL of deionized water, 5 μL of target DNA, and 25 μL of master mix. Amplification products were separated in 1% agar gel and visualized under ultraviolet light after staining with ethidium bromide.

Leishmania. major (MHOM/IR/15/ER) and L. tropica (MHM/SU/79/K27) (13) were used as positive controls. Distilled water instead of the DNA template was used as the negative control. The designed program for Leishmania DNA amplification was an initial denaturation cycle at 95°C for 5 minutes; 30 cycles with cycle 1 at 95°C for 20 seconds (denaturation), cycle 2 at 53°C for 30 seconds (annealing), cycle 3 at 72°C for 1 minute (extension), and a final extension at 72°C for 7 minutes.

3.4. Restriction Fragment Length Polymorphism

The PCR products from ITS1-DNA PCR were digested using the HaeIII restriction endonuclease enzyme and the Taq1 restriction enzyme, as recommended by the manufacture (New England BioLabs, Inc.). The ITS1 sequence was amplified from all isolates, and reference stains were digested separately with restriction enzymes, according to the manufacturer’s instructions.

4. Results

We used the PCR-RFLP assay of ITS1 genes for the direct identification of Leishmania species in 24 out of 33 patients suspected to have CL. These 24 patients, who were positive on smears and cultures for L. donovani bodies demonstrated microscopically, were selected for the study. The patient characteristics are shown in Table 1, including gender, age, number of lesions, ulcer duration, and distribution of ulcers on the body.

Table 1.

Main Characteristics of 24 CL-Positive Patients

Data SamplesAge, yGenderNumber of UlcersUlcer duration, mUlcer LocationLeishmania spp.Genotype Pattern
19M22HandL. majorLmA
212M11LegL. majorLmA
38M33Hand and legL. majorLmA
45F32Hand and legL. majorLmA
510M43Hand, leg, and trunkL. majorLmA
615F32LegL. majorLmA
718M22LegL. majorLmA
88F22NeckFlagellates-
925M33Hand and faceL. majorLmA
106M11HandFlagellates-
1123F22HandL. majorLmA
1217F13HandL. majorLmA
139M22Hand and legL. majorLmA
1417M22Hand and faceL. majorLmA
1529F24HandL. majorLmA
167M32Hand, leg, and trunkL. majorLmA
1721M24HandL. majorLmA
1812F21HandL. majorLmB
199M42Hand and legL. majorLmA
2015M23LegCrithidia-
216F11HandL. majorLmB
2228M32HandL. majorLmA
2317M22Hand and legL. majorLmA
247M32FaceL. majorLmA

With the use of LITSR/L5.8S primer for PCR amplification, single products of the expected size of 350 bp of L. major were detected in 21 patients (Figure 1). In three isolates, a nonspecific 450 bp band was seen. RFLP was performed on the ITS1 region in the ribosomal operons of 24 isolates of L. major. After using the restriction enzyme HaeIII, banding patterns, including fragments of 140 bp and 210 bp bands, were observed in 19 cases (LmA). In three cases with this enzyme, 300 bp and 150 bp bands (Figure 2) were observed, probably flagellates. No digestion was observed in two cases.

Electrophoresis Results of ITS1-PCR from Smears
M, molecular size marker (50 bp); lane 1, standard L. major (MRHO/IR/75/ER); lane 2, standard L. tropica (MHOM/IR/99/YAZ1); lane 3, PCR products (350 bp) as a standard; lane 4, PCR products (450 bp; different from standard); lanes 5 - 10, PCR products; lane 11, negative control.
Restriction Enzyme Digestion Profile of Amplified ITS1 Region with the HaeIII Restriction Enzyme
M, molecular marker (50 bp); lane 1, standard L. major; lane 2, standard L. tropica; lanes 3 - 9, PCR products; lane 10, the only sample different from the normal samples and standard patterns, probably flagellates; lane 11, negative control.

Digestion of ITS1 from ribosomal DNA with Taq1 enzymes enabled us to identify two different patterns. Using the Taq1 restriction enzyme, banding patterns including 150 bp and 200 bp were observed in 19 cases, same as the standard of L. major (genotype pattern LmA). In three isolates, digestion with this enzyme showed 300 bp and 150 bp banding patterns, dissimilar to the L. major reference strain. In two cases, no digestion was observed (LmB) (Figures 3 and 4).

Restriction Enzyme Digestion Profile of Amplified ITS1 Region from L. major with the TaqI Enzyme
M, molecular marker (50 bp); lane 1, standard L. major; lane 2, L. tropica; lanes 3 - 9, PCR products; lane 10, the only sample different from the normal samples and standard patterns, probably flagellates; lane 11, negative control.
Restriction Enzyme Digestion Profile of Amplified ITS1 Region from L. major with the TaqI Enzyme
M, molecular marker (50 bp); lane 1, standard L. major; lane 2, L. tropica; lane 3, PCR products; lane 4, PCR products; lane 5, PCR products; lane 6, the only sample different from the normal samples and standard patterns, probably flagellates; lane 7, PCR products; lanes 8 - 10, PCR products; lane 11, negative control.

5. Discussion

Leishmaniasis is a major societal health problem that mostly affects populations where essential health services are not easily accessible and the detection of CL is based on clinical characteristics (14). The present study shows that CL is caused by L. major in Chabahar, Iran. Parasitological methods are the gold standard for the diagnosis of CL, with a high dependence on the number of parasites in the samples and a requirement of technical skill for sampling. The sensitivity of these methods varies from 27% to 85% for the diagnosis of leishmaniasis, which is their main disadvantage (15). In addition, contamination of the cultured medium is an occasional problem; such material must be discarded.

Molecular techniques have proved to be sensitive and powerful tools for detecting Leishmania directly in clinical samples, as well as for parasite characterization using PCR (1). When PCR is appropriately applied, genetic information can be obtained (16). Some studies on Viannia isolates from different hosts and geographic areas have found high levels of intra- and interspecific variation. Schonian et al. established a diagnostic ITS1 PCR-RFLP method using the restriction enzyme HaeIII for leishmaniasis, combining high sensitivity for detecting Leishmania directly in clinical materials and the ability to identify all medically relevant species groups (12). In our study, using the HaeIII restriction enzyme, the restriction analysis of the amplified ITS1 ribosomal DNA revealed two different schizodeme patterns (LmA and LmB) in L. major.

On the other hand, some researchers have found that HaeIII was not sufficient to distinguish all species of Viannia (16, 17). We applied more restriction enzymes in the diagnosis of Leishmania species using ITS1-RFLP PCR. We found that HaeIII, Rsa, and Hinf1 can produce two profile bandings; Alu1 produces three; Dde1 and Scrf1 produce four; and Taq1 produces five (1). In this study, similar to the studies of de Almeida et al. (18), Doudi et al. (19), and Mirzaie et al. (20), two profile bandings (220 bp and 140 bp) in HaeIII digestion were detected, similar to the standard. The most prevalent genotypes related to isolates in Chabahar were LmA. Taq1 digestion showed two patterns with 150 bp and 200 bp bands. In contrary, in a study in Yazd (20) on molecular identification using Taq1 restriction enzymes, four patterns were produced. This situation must be related to special genotypes in some regions. It should be emphasized that leishmaniasis in Chabahar is new, and genetic variation is limited.

The results of the present study showed that CL infections were more common in males (48.4%) than in females (24.2%) (Table 2), probably because males spend more time outdoors than females do in this region. Our results showed relatively good concordance between parasitological results, cultures, and PCR. Some studies have shown that PCR methods achieved positive results in more than 90% of CL cases (13). The Crithidias are able to endosymbiotically live within some parasites related to the Trypanosomatidae family and sub-families of Leishmania, and move to different hosts along with them. Some researchers believe that the only way to isolate Crithidia from Leishmania is through kDNA and/or ribosomal RNA genes. It has been suggested that by extracting kDNA and/or rDNA, then amplifying them with PCR using the PCR-RFLP technique with the HaeIII enzyme and isolation of bands at 350 bp for Leishmania and at 450 bp for Crithidia, these two parasites can be isolated from each other (21). However, ITS1-RFLP PCR showed good sensitivity, and two ITS1 analyses showed the same results, as in a recent study by Doudi et al. (21).

Table 2.

PCR Diagnostic Ratio of CL in Chabahar According to Sex

SamplesFemale, No. (%)Male, No. (%)
Positive8 (24.2)14 (48.4)
Negative2 (6.06)7 (21.2)
Total10 (30.26)23 (69.6)

Thus, the results of the present study confirm that in the border area of Sistan Va Baluchestan, on the peninsula of Chabahar, L. major is the causative agent of CL, with genotypic variation in this area.

References

  • 1.

    Dabirzadeh M, Mirmohammad Sadeghi H, Baghaie M, Hejazi H. Genetic polymorphism of Leishmania major in two hyper endemic regions of Iran revealed by PPIP-PCR and ITS- RFLP. Arch Iran Med. 2012;15(3):151-6. [PubMed ID: 22369303].

  • 2.

    Baghaei M. Intraspecific variation in Leishmania major isolated from different forms of zoonotic cutaneous Leishmaniasis. Iran J Med Sci. 2015;30(2).

  • 3.

    Monroy-Ostria A, Nasereddin A, Monteon VM, Guzman-Bracho C, Jaffe CL. ITS1 PCR-RFLP Diagnosis and Characterization of Leishmania in Clinical Samples and Strains from Cases of Human Cutaneous Leishmaniasis in States of the Mexican Southeast. Interdiscip Perspect Infect Dis. 2014;2014:6072-87. [PubMed ID: 25104958]. https://doi.org/10.1155/2014/607287.

  • 4.

    Monteiro Tinel JM, Benevides MF, Frutuoso MS, Rocha CF, Arruda FV, Vasconcelos MA, et al. A lectin from Dioclea violacea Interacts with midgut surface of Lutzomyia migonei, unlike its homologues, Cratylia floribunda lectin and Canavalia gladiata lectin. Sci World J. 2014;2014:1-7. [PubMed ID: 25431778]. https://doi.org/10.1155/2014/239208.

  • 5.

    Gramiccia M, Gradoni L. The current status of zoonotic leishmaniases and approaches to disease control. Int J Parasitol. 2005;35(11-2):1169-80. [PubMed ID: 16162348].

  • 6.

    Gibson LJ, McCann SH, Brewster S, Barker DC. A kinetoplast DNA minicircle probe (B18) specific for the Leishmania subgenus Viannia. Trans R Soc Trop Med Hyg. 2002;96 Suppl 1:S71-3. [PubMed ID: 12055854].

  • 7.

    Souto RP, Fernandes O, Macedo AM, Campbell DA, Zingales B. DNA markers define two major phylogenetic lineages of Trypanosoma cruzi. Mol Biochem Parasitol. 1996;83(2):141-52. [PubMed ID: 9027747].

  • 8.

    Victoir K, Banuls AL, Arevalo J, Llanos-Cuentas A, Hamers R, Noel S, et al. The gp63 gene locus, a target for genetic characterization of Leishmania belonging to subgenus Viannia. Parasitology. 1998;117 ( Pt 1):1-13. [PubMed ID: 9695095].

  • 9.

    Bensoussan E, Nasereddin A, Jonas F, Schnur LF, Jaffe CL. Comparison of PCR assays for diagnosis of cutaneous leishmaniasis. J Clin Microbiol. 2006;44(4):1435-9. [PubMed ID: 16597873]. https://doi.org/10.1128/JCM.44.4.1435-1439.2006.

  • 10.

    Reithinger R, Dujardin JC. Molecular diagnosis of leishmaniasis: current status and future applications. J Clin Microbiol. 2007;45(1):21-5. [PubMed ID: 17093038]. https://doi.org/10.1128/JCM.02029-06.

  • 11.

    Cupolillo E, Grimaldi Junior G, Momen H, Beverley SM. Intergenic region typing (IRT): a rapid molecular approach to the characterization and evolution of Leishmania. Mol Biochem Parasitol. 1995;73(1-2):145-55. [PubMed ID: 8577322].

  • 12.

    Schonian G, Nasereddin A, Dinse N, Schweynoch C, Schallig HD, Presber W, et al. PCR diagnosis and characterization of Leishmania in local and imported clinical samples. Diagn Microbiol Infect Dis. 2003;47(1):349-58. [PubMed ID: 12967749].

  • 13.

    Tashakori M, Kuhls K, Al-Jawabreh A, Mauricio IL, Schonian G, Farajnia S, et al. Leishmania major: genetic heterogeneity of Iranian isolates by single-strand conformation polymorphism and sequence analysis of ribosomal DNA internal transcribed spacer. Acta Trop. 2006;98(1):52-8. [PubMed ID: 16513079]. https://doi.org/10.1016/j.actatropica.2006.01.010.

  • 14.

    Khan SN, Ayaz S, Khan S, Jan AH, Attaullah S. Molecular epidemiology of cutaneous leishmaniasis in Southern belt of Khyber Pakhtunkhwa, Pakistan. IJOART. 2013;2(4):472-81.

  • 15.

    Hajjaran H, Vasigheh F, Mohebali M, Rezaei S, Mamishi S, Charedar S. Direct diagnosis of Leishmania species on serosity materials punctured from cutaneous leishmaniasis patients using PCR-RFLP. J Clin Lab Anal. 2011;25(1):20-4. [PubMed ID: 21254238]. https://doi.org/10.1002/jcla.20377.

  • 16.

    Belli A, Rodriguez B, Aviles H, Harris E. Simplified polymerase chain reaction detection of new world Leishmania in clinical specimens of cutaneous leishmaniasis. Am J Trop Med Hyg. 1998;58(1):102-9. [PubMed ID: 9452300].

  • 17.

    Spanakos G, Piperaki ET, Menounos PG, Tegos N, Flemetakis A, Vakalis NC. Detection and species identification of Old World Leishmania in clinical samples using a PCR-based method. Trans R Soc Trop Med Hyg. 2008;102(1):46-53. [PubMed ID: 17669452]. https://doi.org/10.1016/j.trstmh.2007.05.019.

  • 18.

    de Almeida ME, Steurer FJ, Koru O, Herwaldt BL, Pieniazek NJ, da Silva AJ. Identification of Leishmania spp. by molecular amplification and DNA sequencing analysis of a fragment of rRNA internal transcribed spacer 2. J Clin Microbiol. 2011;49(9):3143-9. [PubMed ID: 21752983]. https://doi.org/10.1128/JCM.01177-11.

  • 19.

    Doudi M, Hejazi SH, Razavi MR, Narimani M, Khanjani S, Eslami G. Comparative molecular epidemiology of Leishmania major and Leishmania tropica by PCR-RFLP technique in hyper endemic cities of Isfahan and Bam, Iran. Med Sci Monit. 2010;16(11):CR530-5. [PubMed ID: 20980956].

  • 20.

    Mirzaie F, Eslami G, Yosefi MH, Pestehchian N. Molecular identification of Leishmania isolates obtained from patients suspected as having cutaneous leishmaniasis referred to reference laboratories from Yazd province in central Iran. Adv Biomed Res. 2013;2:92. [PubMed ID: 24524036]. https://doi.org/10.4103/2277-9175.122525.

  • 21.

    Doudi M, Setorki M, Eslami G, Hejazi SH. The Molecular Epidemiology of Crithidia Strains in Isfahan City and Surrounding Areas, Iran [in Persian]. J Isfahan Med Sch. 2014;32(275).