- Open Access
High resolution melting analysis and detection of Leishmania resistance: the role of multi drug resistance 1 gene
Genes and Environment volume 43, Article number: 36 (2021)
Pentavalent antimonial compounds are currently used to treat leishmaniasis and resistance to these drugs is a serious problem. Multidrug resistance protein is an efflux pump of the cell membrane that expels foreign compounds. This study designed to evaluate the mutations in the multi-drug resistance 1 (MDR1) gene, in biopsy specimens of Leishmania tropica, with high resolution melting (HRM) method. In this experimental study, genomic DNA was extracted from 130 patients with skin leishmaniasis. Then, nucleotide changes were investigated throughout the gene using HRM and sequencing methods. The samples categorized in 5 groups by differences in the melting temperature (Tm).
The nucleotide changes analysis showed that 61% of the samples of different groups that were unresponsive to drug had mutations in the MDR1 gene, which were also confirmed by the sequencing method. These mutations can be one of the factors responsible for non-responsiveness to the treatment.
According to the findings, it seems that mutation in MDR1 gene could be responsible for drug resistance to pentavalent antimonial compounds. Furthermore, HRM method can be used to diagnose drug resistance in leishmaniasis. It is also recommended that further studies be done regarding the importance of drug resistance in the leishmania affected patients.
Cutaneous leishmaniasis is one of the major diseases in tropical and subtropical areas  which is caused by parasites of Leishmania tropica and Leishmania major [2,3,4,5]. It is a health problem in 88 countries and its controlling and preventing is one of the top priorities of the international community, including the World Health Organization [6,7,8]. Currently, pentavalent antimony compounds, including glucantime and sodium stibogluconate (SSG), are used to treat leishmaniasis and through various direct and indirect mechanisms may eliminate this parasite. Besides inhibition of the activity of enzymes involving in the fatty acid oxidation pathway [9, 10] and stimulation of macrophages to produce lethal molecules such as nitric oxide (NO) and reactive oxygen species (ROS) , these compounds via creating the flow of intracellular thiol (decrease of parasitic thiol), the disturbance in the oxidation potential and by inhibiting trypanothione reductase lead to oxidative stress-dependent parasite death .
Studies on the failure of a pentavalent antimony treatment have recently been reported in several countries . Since uncovering the molecular mechanisms involved in drug resistance will help to treat patients, efforts to identify these mechanisms are ongoing.
Recent our experiment evidenced that expression of multidrug resistance protein 1 (MDR1) or CD243 is different in various forms of cutaneous leishmaniasis . This ATP-dependent drug efflux protein pumps several substances out of cells and can reduce drug accumulation in resistant cells . About 65% of the Mediterranean leishmania strains have amplified MDR1 genes . In addition, it has been made clear that overexpression of this gene has been associated with an increased resistance of leishmania to miltefosine . Overall, the mutation of this gene appears to play a key role in the phenomenon of drug resistance [16,17,18]. Therefore, screening of patients with this mutation has an important role in the treatment of patients with leishmaniasis. Due to the limitations of gold standard method (sequencing), trying to identify simple, feasible, and high-performance methods is inevitable. One of these methods is high resolution melting analysis (HRM). This new, accurate, and affordable method is a post-PCR analysis of genetic variation that performed in a sealed tube [19, 20].
Given the importance of MDR1 gene mutations in the drug resistance of leishmania and the lack of similar study, in the present experiment, we evaluated the capacity of HRM assay for detection of MDR1 gene mutation in genomic DNA of tissue specimens of patients with skin leishmaniasis.
Materials and methods
Tissue specimens from 130 patients (chemoresistant samples: 61% (80 samples) and chemosensitive samples: 39% (50 samples)) with cutaneous leishmaniasis along with one standard sample (sensitive strain of leishmania tropica; MHOM/IR/10/175) were obtained from Dermatology Department of Afzalipour Hospital in Kerman, Iran.
Genomic DNA was isolated by using QIAamp DNA Mini Kit (QIAGEN, Germany). We used the manufacturer’s protocol for extraction. DNA concentration was determined by Nano Drop® at a wavelength of 260 nm. Also, the absorption ratio of a pure DNA at 260 and 280 nm was calculated (A260/280).
Specific primers were designed to amplify the MDR1 encoded region of the MDR gene. For this purpose, we aligned the sequence of HM854717.1 through HM854725.1 (partial sequence) with the whole sequence of the MDR1 (U63320.1). The forward primer sequence was selected from the nucleotides before the conserved points of the 10 mentioned sequences (HM854717.1-HM854725.1 and whole sequence of the MDR1). The reverse primer is selected among sequences that located after the conserved points (160 nucleotides after the forward primer). The designed primers had the following sequence and were synthesized by TibMolbiol (Berlin, Germany):
Real-time PCR was performed using 10 μL of Type-it Master Mix (QIAGEN, Germany), 2 μL of DNA, 0.7 μL of each of forward and reverse primers (10 pmol) and 6.6 μL of nuclease-free water in a total volume of 20 μL. Thermal cycling conditions included an initial activation step at 95 °C for 5 min followed by 40 cycles including a denaturation step at 94 °C for 20 s and a combined annealing/elongation step at 60 °C for 30 s. The reaction took place in rotor gene 6000 (QIAGEN, Germany). For HRM analysis, the PCR products were melted by warming up the temperature from 40 °C to 95 °C (0.007 °C.s− 1) and 20 fluorescence acquisitions were recorded in each temperature. The normalized graph and the normalized temperature-shifted difference graph (difference graph) from the gene scanning analysis were used to analyze the data. HRM data were analyzed using rotor gene 6000 (QIAGEN, Germany) software. Based on the melting temperature (Tm), studied samples were categorized in five groups. Namely, group 1 was the standard and treatment-sensitive samples and each of the other groups also illustrated samples with similar melting temperature.
Agarose gel electrophoresis
After amplification, the amplicon was electrophoresis in agarose gel. The gel (1.5%) was prepared in Tris/Borate/EDTA (TBE) buffer and intercalating dye (SYBR Safe) added to gel solution to make the DNA visible under UV. After casting the gel, creating wells, and sinking the gel in TBE buffer, PCR products were loaded into wells. The DNA ladder was used to determine the sizes of sample bands.
To confirm nucleotide changes and compare the nucleotide differences among samples of different groups, some of the samples (10% of each group) were amplified by MDR1-R and MDR1-F primers and bilateral sequencing was accomplished by Bioneer Company in South Korea. The results of sequencing were studied by the choromas pro.v.2.1 software and sequence alignment was carried out. Alignment diagram was drawn by the use of Vector NTI advance 11 software.
Amplification and electrophoresis
The amplification was confirmed by conventional PCR and gel electrophoresis. As illustrated in Fig. 1, all samples had 158 bp amplification products.
Real-time PCR and HRM results
Corresponding primers amplified a 158 bp product. In HRM analysis, all samples had melting peaks with certain Tm. Different Tm can be attributed to genomic differences in the studied area. To further differentiate, a high-resolution difference curve was plotted by subtracting the melting curve of each sample from the baseline (reference) curve. Accordingly, samples with similar curves were placed in separate groups (Fig. 2).
To compare the amplified products, some of the samples were sent for sequencing. Next, the nucleotide sequence of different groups was subjected to alignment with reference sequence. In this method, each nucleotide of a sample is compared with all nucleotide of the reference sample and any similarity is shown by a dark dot. Since the highest difference was obtained between groups 1 and 3, so, the Fig. 3 is representing their corresponding alignment chart. As mentioned before, different melting temperature can be attributed to genomic differences in the studied area. As presented in Fig. 4, some insertion and mutation were detected in the sequence of MDR1 gene of most of the patients that we studied.
Currently, chemotherapy using the pentavalent antimonials (such as glucantime) is applied more than other methods in the treatment of leishmaniasis. However, the use of glucantime in patients with cutaneous leishmaniasis has some side effects on vital organs such as kidney, liver, heart, and blood [18, 21,22,23,24]. Studies conducted in different parts of the world showed drug resistance to glucantime in some people with leishmaniasis [24,25,26] and had caused a lot of concern after about 65 years of usage. Although the resistance mechanisms are poorly understood, different factors such as genetic factors, proteinaceous factors, enzymatic factors, and alteration of drug uptake and accumulation appear to be effective in causing this condition [17, 18, 25,26,27,28,29].
ATP-binding cassette (ABC) transporters such as P-glycoprotein (P-gp) are transmembrane proteins, which is encoded by the MDR1 gene. P-gp is an ATP-dependent drug efflux pump which can reduce the drug’s accumulation in resistant cells and often increases the drug resistance [30,31,32]. The MDR1 gene of leishmania is 37% similar to human MDR gene and located on chromosome 34. Heretofore, the association of increased resistance to several drugs with overexpression of this gene has been proven [14, 33,34,35]. In line with these studies, our experiment in 2019 showed that the expression of this gene is different in various forms of the disease (acute, chronic lupoid, and chronic non-lupoid) and can be considered as an important factor in the drug resistance of leishmania parasite .
Although sequencing is considered as a gold standard approach for detection of mutations, the major obstacles are overpriced instrumentation, inaccessibility, and need to striking technical expertise. Therefore, and according to the importance of drug resistance, modern alternative techniques can provide an opportunity for leishmania control and treatment. Flow cytometry, RFLP-PCR, and HRM are among these techniques. For instance, Alizadeh and colleagues by using RFLP-PCR method illustrated that about 11% of the collected samples of cutaneous leishmaniasis from some endemic areas of Iran (Yazd, Mashhad, and Kashan provinces) had mutations in MDR1 gene . Tsirigotakis et al. also evaluated drug resistance in leishmania isolates using flow cytometry . Since Rhodamine-123 is an established substrate for P-gp 170 and its efflux rate is largely dependent on the number of efflux pumps, they introduced the efflux rate of this dye as an indicator of drug resistance.
Previously, HRM assay with direct sequencing confirmation possibility was reported as a rapid and inexpensive method for detection of mutations in pathological . Furthermore, this test can be used as a pre-screening method to decide whether sequencing should be done on a sample. To discriminate between wild-type and mutant DNA, the genomic DNA is amplified with specific primers and in the presence of a fluorescent intercalating dye. Then, the fluorescent melting curve of amplicons is analyzed . Regarding validation, Heideman and co-workers challenged the performance of HRM assay for detection of mutations in genomic DNAs of 68 non-small cell lung cancer specimens . They evidenced that this assay had a sensitivity of ≤5% of mutated DNA in a background of wild-type DNA and could be applied for screening of mutations.
In the present experiment, we evaluated the capacity of HRM assay for detection of MDR1 gene mutation in genomic DNA of tissue specimens of 130 patients with skin leishmaniasis. Subsequent sequencing of the HRM products was applied to confirm the presence of the nucleotide alteration.
Figure 1 illustrated that MDR1 gene was amplified in our samples. Based on the Tm, studied samples were categorized in five groups. Group 1 actually refers to standard sample and samples of treatment-sensitive patients. Each of the other groups also illustrates samples with similar melting temperature (Fig. 2). Our results showed that 61% of the specimens had mutations in the MDR1 gene region. To confirm findings and compare the amplified products, sequencing and then alignment analysis was applied (Figs. 3 and 4). The analysis showed that the type of mutations in each group is similar but different from the mutations in the other groups. Interestingly, no discrepancy was found between HRM and sequencing findings as well as triplicate analysis confirmed the reproducibility of HRM assay.
Totally, the present experiment established that HRM assay is sensitive and suitable enough to screen mutations of MDR1 gene. Furthermore, rapidity, simplicity, feasibility, and cost-effectiveness are other advantages that support HRM use to detect mutations. In addition, unlike methods such as RFLP-PCR, HRM assay detects mutations anywhere in the amplicon and does not require knowing the exact location of the mutation. Nevertheless, this assay could not identify the specific nucleotide change. It should be noted that this inability is predictable for a screening test. As a suggestion, additional experiments for more validation of HRM assay, as a screening test for detection of resistant patients, as well as determination of the prevalence of genetic factors affecting the resistance of antimony compounds in endemic areas of cutaneous leishmaniasis are necessary.
The present study has potential limitations. The lack of information about the duration, amount, and type of drugs which were used to treat the studied patients can be considered as the most important ones. It should be noted that due to the overlap of melting curves, identifying single-point mutations whose melting temperatures do not differ much can be challenging.
Availability of data and materials
Please contact corresponding author (S.D) for data requests.
Multi-drug resistance 1
High resolution melting
Reactive oxygen species
Polymerase chain reaction
Kubba R, Al-Gindan Y, El-Hassan A, Omer A. Clinical diagnosis of cutaneous leishmaniasis (oriental sore). J Am Acad Dermatol. 1987;16(6):1183–9. https://doi.org/10.1016/S0190-9622(87)70155-8.
Mirzaei A, Rouhani S, Taherkhani H, Farahmand M, Kazemi B, Hedayati M, et al. Isolation and detection of Leishmania species among naturally infected Rhombomis opimus, a reservoir host of zoonotic cutaneous leishmaniasis in Turkemen Sahara, North East of Iran. Exp Parasitol. 2011;129(4):375–80.
Alborzi A, Pouladfar GR, Fakhar M, Motazedian MH, Hatam GR, Kadivar MR. Isolation of Leishmania tropica from a patient with visceral leishmaniasis and disseminated cutaneous leishmaniasis, southern Iran. Am J Trop Med Hyg. 2008;79(3):435–7. https://doi.org/10.4269/ajtmh.2008.79.435.
Karamian M, Motazedian MH, Mehrabani D, Gholami K. Leishmania major infection in a patient with visceral leishmaniasis: treatment with amphotericin B. Parasitol Res. 2007;101(5):1431–4. https://doi.org/10.1007/s00436-007-0649-x.
Fekrisoofiabadi M, Fekri M, Moradabadi A, Vahidi R, Khaleghi M, Ram M, et al. Evaluation of MDR1 and MRPA genes expression in different types of dry cutaneous leishmaniasis. BMC Res Notes. 2019;12(1):1–4.
Solbach W, Laskay T. The host response to Leishmania infection. Adv Immunol. 1999;74:275–317.
Handman E. Leishmaniasis: current status of vaccine development. Clin Microbiol Rev. 2001;14(2):229–43. https://doi.org/10.1128/CMR.14.2.229-243.2001.
Moradabadi A, Farsinejad A, Fekri SM. Fast method for diagnosis of leishmania by PCR and FLASH PCR. J Arak Univ Med Sci. 2017;19(11):79–86.
Sundar S, More DK, Singh MK, Singh VP, Sharma S, Makharia A, et al. Failure of pentavalent antimony in visceral leishmaniasis in India: report from the center of the Indian epidemic. Clin Infect Dis. 2000;31(4):1104–7. https://doi.org/10.1086/318121.
Mousavi Z, Yazdani Z, Moradabadi A, Hoseinpourkasgari F, Hassanshahi G. Role of some members of chemokine/cytokine network in the pathogenesis of thalassemia and sickle cell hemoglobinopathies: a mini review. Exp Hematol Oncol. 2019;8(1):1–6.
Frézard F, Demicheli C, Ribeiro RR. Pentavalent antimonials: new perspectives for old drugs. Molecules. 2009;14(7):2317–36. https://doi.org/10.3390/molecules14072317.
Vanaerschot M, Decuypere S, Berg M, Roy S, Dujardin J-C. Drug-resistant microorganisms with a higher fitness–can medicines boost pathogens? Crit Rev Microbiol. 2013;39(4):384–94. https://doi.org/10.3109/1040841X.2012.716818.
Neves JI. Avaliação da expressão do gene MDR1 (Glicoproteína-P) e atividade de efluxo em células do sangue periférico de pacientes sob tratamento da tuberculose multirresistente; 2013.
Castanys-Munoz E, Pérez-Victoria JM, Gamarro F, Castanys S. Characterization of an ABCG-like transporter from the protozoan parasite Leishmania with a role in drug resistance and transbilayer lipid movement. Antimicrob Agents Chemother. 2008;52(10):3573–9. https://doi.org/10.1128/AAC.00587-08.
Tiwari N, Gedda MR, Tiwari VK, Singh SP, Singh RK. Limitations of current therapeutic options, possible drug targets and scope of natural products in control of leishmaniasis. Mini-Rev Med Chem. 2018;18(1):26–41. https://doi.org/10.2174/1389557517666170425105129.
Rai K, Cuypers B, Bhattarai NR, Uranw S, Berg M, Ostyn B, et al. Relapse after treatment with miltefosine for visceral leishmaniasis is associated with increased infectivity of the infecting Leishmania donovani strain. MBio. 2013;4(5).
Messaritakis I, Christodoulou V, Mazeris A, Koutala E, Vlahou A, Papadogiorgaki S, et al. Drug resistance in natural isolates of Leishmania donovani sl promastigotes is dependent of Pgp170 expression. PLoS One. 2013;8(6):e65467. https://doi.org/10.1371/journal.pone.0065467.
Fekri-SoofiAbadi M, Fekri M, Vahidi R, Shamsi-Meymandi S, Dabiri D, Dabiri S. Ability of real-time PCR for differential diagnosis of various forms of cutaneous leishmaniasis: a comparative study with histopathology. BMC Res Notes. 2019;12(1):1–5.
Moradabadi A, Farsinejad A, Khansarinejad B, Fatemi A. Development of a high resolution melting analysis assay for rapid identification of JAK2 V617F missense mutation and its validation. Exp Hematol. 2019;8(1):1–7.
Moradabadi A, Fatemi A, Noroozi-Aghideh A. Analysis of the reannealing-instead of melting-curve in the detection of JAK2 V617F mutation by HRM method. J Blood Med. 2019;10:235–41. https://doi.org/10.2147/JBM.S204222.
Croft SL, Coombs GH. Leishmaniasis–current chemotherapy and recent advances in the search for novel drugs. Trends Parasitol. 2003;19(11):502–8. https://doi.org/10.1016/j.pt.2003.09.008.
Demicheli C, Ochoa R, da Silva JB, Falcão CA, Rossi-Bergmann B, de Melo AL, et al. Oral delivery of meglumine antimoniate-β-cyclodextrin complex for treatment of leishmaniasis. Antimicrob Agents Chemother. 2004;48(1):100–3. https://doi.org/10.1128/AAC.48.1.100-103.2004.
Sadr SF, Talari SA. Effect of Glucantime on electrocardiogram of patients suffering from cutaneous Leishmaniasis. KAUMS J (FEYZ). 1998;2(1):13–20.
Talari SA, Vakili Z, Emami AH. Effect of Glucantime on blood parameters in patients with cutaneous Leishmaniasis. KAUMS J (FEYZ). 1999;3(2):17–22.
Croft SL, Sundar S, Fairlamb AH. Drug resistance in leishmaniasis. Clin Microbiol Rev. 2006;19(1):111–26. https://doi.org/10.1128/CMR.19.1.111-126.2006.
Hadighi R, Mohebali M, Boucher P, Hajjaran H, Khamesipour A, Ouellette M. Unresponsiveness to Glucantime treatment in Iranian cutaneous leishmaniasis due to drug-resistant Leishmania tropica parasites. PLoS Med. 2006;3(5):e162. https://doi.org/10.1371/journal.pmed.0030162.
Kristensen LS, Wojdacz TK, Thestrup BB, Wiuf C, Hager H, Hansen LL. Quality assessment of DNA derived from up to 30 years old formalin fixed paraffin embedded (FFPE) tissue for PCR-based methylation analysis using SMART-MSP and MS-HRM. BMC Cancer. 2009;9(1):453. https://doi.org/10.1186/1471-2407-9-453.
Liew M, Pryor R, Palais R, Meadows C, Erali M, Lyon E, et al. Genotyping of single-nucleotide polymorphisms by high-resolution melting of small amplicons. Clin Chem. 2004;50(7):1156–64. https://doi.org/10.1373/clinchem.2004.032136.
Croft SL. Monitoring drug resistance in leishmaniasis. Tropical Med Int Health. 2001;6(11):899–905. https://doi.org/10.1046/j.1365-3156.2001.00754.x.
Coelho AC, Tosi LR, Cotrim PC. Mapping of a Leishmania major gene/locus that confers pentamidine resistance by deletion and insertion of transposable element. Rev Inst Med Trop Sao Paulo. 2004;46(2):109–12. https://doi.org/10.1590/S0036-46652004000200010.
Coelho AC, Beverley SM, Cotrim PC. Functional genetic identification of PRP1, an ABC transporter superfamily member conferring pentamidine resistance in Leishmania major. Mol Biochem Parasitol. 2003;130(2):83–90. https://doi.org/10.1016/S0166-6851(03)00162-2.
Légaré D, Richard D, Mukhopadhyay R, Stierhof Y-D, Rosen BP, Haimeur A, et al. The Leishmania ATP-binding cassette protein PGPA is an intracellular metal-thiol transporter ATPase. J Biol Chem. 2001;276(28):26301–7. https://doi.org/10.1074/jbc.M102351200.
Ponte-Sucre A, Gamarro F, Dujardin J-C, Barrett MP, López-Vélez R, García-Hernández R, et al. Drug resistance and treatment failure in leishmaniasis: a 21st century challenge. PLoS Negl Trop Dis. 2017;11(12):e0006052. https://doi.org/10.1371/journal.pntd.0006052.
Henderson DM, Sifri CD, Rodgers M, Wirth DF, Hendrickson N, Ullman B. Multidrug resistance in Leishmania donovani is conferred by amplification of a gene homologous to the mammalian mdr1 gene. Mol Cell Biol. 1992;12(6):2855–65. https://doi.org/10.1128/mcb.12.6.2855-2865.1992.
Pérez-Victoria JM, Bavchvarov BI, Torrecillas IR, Martínez-García M, López-Martín C, Campillo M, et al. Sitamaquine overcomes ABC-mediated resistance to miltefosine and antimony in Leishmania. Antimicrob Agents Chemother. 2011;55(8):3838–44. https://doi.org/10.1128/AAC.00065-11.
Alizadeh R, Hooshyar H, Bandehpor M, Arbabi M, Kazemi F, Talari A, et al. Detection of drug resistance gene in cutaneous leishmaniasis by PCR in some endemic areas of Iran. Iran Red Crescent Med J. 2011;13(12):863–7.
Tsirigotakis N, Christodoulou V, Ntais P, Mazeris A, Koutala E, Messaritakis I, et al. Geographical distribution of MDR1 expression in Leishmania isolates, from Greece and Cyprus, measured by the rhodamine-123 efflux potential of the isolates, using flow cytometry. Am J Trop Med Hyg. 2016;94(5):987–92. https://doi.org/10.4269/ajtmh.15-0658.
Heideman D, Thunnissen F, Doeleman M, Kramer D, Verheul H, Smit E, et al. A panel of high resolution melting (HRM) technology-based assays with direct sequencing possibility for effective mutation screening of EGFR and K-ras genes. Anal Cell Pathol. 2009;31(5):329–33. https://doi.org/10.1155/2009/465850.
We would like to thank the staff of the Pathology Department of Kerman University of Medical Sciences and Stem Cell Research Center who assisted in performing this study.
No funding sources used in this study.
Ethics approval and consent to participate
The study approved in Ethical Committee of Kerman University of Medical Sciences. The Ethical Approval Code is IR.KMU.REC.1397.260.
Consent for publication
All authors consent for publication in BMC Genes and Environment journal.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Abadi, M.F.S., Moradabadi, A., Vahidi, R. et al. High resolution melting analysis and detection of Leishmania resistance: the role of multi drug resistance 1 gene. Genes and Environ 43, 36 (2021). https://doi.org/10.1186/s41021-021-00210-5
- Unresponsive leishmaniasis
- HRM method
- Pentavalent antimonial compounds
- Multi-drug resistance 1 gene