Acessibilidade / Reportar erro

Two Leishmania species separation targeting the ITS-rDNA and Cyt b genes by developing and evaluating HRM- qPCR

ABSTRACT

Background:

Incidence of Cutaneous Leishmaniasis as an infectious and neglected disease is increasing, for the diagnosis of which several traditional methods and conventional PCR techniques have been developed, employing different genes for species identification.

Methods:

Leishmania parasites were sampled, DNA was extracted, and new specific and sensitive primers were designed. Two ITS-rDNA and Cyt b genes were targeted by qPCR using the High- Resolution Melting method to identify Leishmania parasites. The standard curves were drawn, compared, and identified by high-resolution melting curve analysis.

Results:

Melting temperature and Cycle of Threshold of ITS-rDNA was higher than Cyt b but Cyt b was more sensitive than ITS-rDNA when Leishmania major and Leishmania tropica were analyzed and evaluated. By aligning melt curves, normalizing fluorescence curves, and difference plotting melt curves, each Leishmania species was distinguished easily. L. major and L. tropica were separated at 83.6 °C and 84.7 °C, respectively, with less than 0.9 °C of temperature difference. Developing sensitivity and specificity of real-time PCR based on EvaGreen could detect DNA concentration to less than one pmol.

Conclusions:

Precise identification of Leishmania parasites is crucial for strategies of disease control. Real-time PCR using EvaGreen provides rapid, highly sensitive, and specific detection of parasite’s DNA. The modified High-Resolution Melting could determine unique curves and was able to detect single nucleotide polymorphisms according to small differences in the nucleotide content of Leishmania parasites.

Keywords:
Cyt b; ITS-rDNA; PCR; Real-time PCR; HRM

INTRODUCTION

Leishmaniasis, one of the neglected tropical and subtropical diseases, has been considered based on the medical importance; besides, visceral leishmaniasis is potentially lethal11. Fotouhi-Ardakani R, Ghafari SM, Ready PD, Parvizi P. Developing, modifying, and validating a TaqMan Real-Time PCR technique for accurate identification of Leishmania parasites causing most leishmaniasis in Iran. Front Cell Infect Microbiol. 2021;11:e731595. Available from: https://doi.org/10.3389/fcimb.2021.731595
https://doi.org/10.3389/fcimb.2021.73159...

2. Fotouhi-Ardakani R, Dabiri S, Ajdari S, Alimohammadian MH, AlaeeNovin E, Taleshi N, et al. Assessment of nuclear and mitochondrial genes in precise identification and analysis of genetic polymorphisms for the evaluation of Leishmania parasites. Infect Genet Evol. 2016;46:33-41.
-33. Müller KE, Zampieri RA, Aoki J, Muxel SM, Nerland AH, Floeter-Winter LM. Amino acid permease 3 (AAP3) coding sequence as a target for Leishmania identification and diagnosis of leishmaniases using high resolution melting analysis. Parasit Vectors. 2018;11:421.. From all known Leishmania species, some are medically important. In Cutaneous Leishmaniasis (CL) the lesions caused by different Leishmania species, which are usually large but painless unless a secondary bacterial or fungal infection occurs. The skin manifestations of the disease in final stages can usually lead to obstruction or degradation of the pharynx, larynx, and nose44. Zampieri RA, Laranjeira-Silva MF, Muxel SM, Stocco de Lima AC, Jon Shaw J, Floeter-Winter LM. High Resolution Melting analysis targeting hsp70 as a fast and eficient method for the discrimination of Leishmania species. PLoS Negl Trop Dis. 2016;10(2): e0004485..

CL includes ZCL (Zoonotic Cutaneous Leishmaniasis) and ACL (Anthroponotic Cutaneous Leishmaniasis). These important tropical diseases have many impressive effects on the human health. The mammalian host(s) and human were infected by different Leishmania parasites in the Old World.

The vertebrate animals were bitten by female sand flies and thus developed ZCL. Humans were involved in Leishmania parasites and they were involved ACL by humans oppositely55. Zahirnia AH, Bordbar A, Ebrahimi S, Spotin A, Mohammadi S, Ghafari SM, et al. Predominance of Leishmania major and rare occurrence of Leishmania tropica with haplotype variability at the center of Iran. Braz J Infect Dis. 2018;22(4):278-87.,66. Ready PD. Biology of phlebotomine sand flies as vectors of disease agents. Annu Rev Entomol. 2013;58:227-50..

In general, precise identification of infectious agents of leishmaniasis and their clinical manifestations are needed in reservoir hosts, vectors, and humans. Well known and practical methods for the detection and recognition of Leishmania species were employed. The Giemsa staining, clinical indications, vector estimations, in vivo pathogenicity, growth in media, and isoenzyme were used for distinguishing between different Leishmania parasites. Later PCR, RFLP (Restriction Fragment Length Polymorphism), sequencing, and phylogenetic analysis were applied11. Fotouhi-Ardakani R, Ghafari SM, Ready PD, Parvizi P. Developing, modifying, and validating a TaqMan Real-Time PCR technique for accurate identification of Leishmania parasites causing most leishmaniasis in Iran. Front Cell Infect Microbiol. 2021;11:e731595. Available from: https://doi.org/10.3389/fcimb.2021.731595
https://doi.org/10.3389/fcimb.2021.73159...
,77. Bates PA, Depaquit J, Galati EA, Kamhawi SH, Maroli M, Ann McDowellet M, et al. Recent advances in phlebotomine sand fly research related to leishmaniasis control. Parasit Vectors . 2015;8(131):11-21..

High- Resolution Melting (HRM) is an efficient and cost-effective procedure that eliminates the risk of laboratory contaminations. Real-time PCR using HRM can firmly identify most well-known Leishmania species by employing housekeeping nuclear ITS-rDNA and mitochondrial Cyt b genes33. Müller KE, Zampieri RA, Aoki J, Muxel SM, Nerland AH, Floeter-Winter LM. Amino acid permease 3 (AAP3) coding sequence as a target for Leishmania identification and diagnosis of leishmaniases using high resolution melting analysis. Parasit Vectors. 2018;11:421.,88. Odell ID, Cloud JL, Seipp M, Wittwer CT. Rapid species identification within the Mycobacterium chelonae-abscessus group by high-resolution melting analysis of hsp65 PCR products. Am J Clin Pathol. 2005;123(1):96-101.

9. Pangasa A, Jex AR, Campbell BE, Bott NJ, Whipp M, Hogg G, et al. High resolution melting-curve (HRM) analysis for the diagnosis of cryptosporidiosis in humans. Mol Cell Probes. 2009;23(1):10-5.
-1010. Chua KH, Ng CC, Lee PC, Lee PC, Lian Lim YA, Lau TP, et al. Development of high resolution melting analysis for the diagnosis of human malaria. Cell Rep Phys Sci. 2015;5:e15671.. This method is inexpensive, rapid, and has high-performance detection; it requires no electrophoresis analysis, reducing risk of contamination as well as significantly decreasing workload, time, and cost. HRM is a relatively new analysis technique that allows direct identification of PCR amplicons. In addition, advantage of qPCR over conventional PCR is that initial DNA concentration is determined with high precision and sensitivity. The results using HRM-qPCR can be qualitative and/or quantitative1111. Kaltenboeck B, Wang C. Advances in real-time PCR: application to clinical laboratory diagnostics. Adv Clin Chem. 2005;40:219-59.,1212. Wang W, Zijlstra RT, Gänzle MG. Identification and quantification of virulence factors of enterotoxigenic Escherichia coli by high-resolution melting curve quantitative PCR. BMC Microbiol. 2017;17(1):114..

ITS-rDNA is a housekeeping gene and Cyt b is an enzymatic and expression gene.

For precise identification of Leishmania parasites, two housekeeping ITS-rDNA and mitochondrial Cyt b genes were targeted and HRM-qPCR assay was modified and developed. This was an important purpose of this research. Accuracy, speed, and especially affordability as well as low cost, were the secondary purposes of this study. Recently, it has been shown that the colors of EvaGreen and SYTO rather than SYBR Green are more specific and more susceptible for measuring DNA quantification in qPCR; therefore, this color was used in this investigation1313. Ahmad AI, Ghasemi JB. New unsymmetrical cyanine dyes for real-time thermal cycling. Anal Bioanal Chem. 2007;389(3):983-8.,1414. Eischeid AC. SYTO dyes and EvaGreen outperform SYBR Green in real-time PCR. BMC Res Notes. 2011;4:263..

METHODS

Sampling of Leishmania isolates and DNA extraction

Different Leishmania parasites from different locations of Iran were employed including four strains that were represented in more than 160 collected samples. Following extraction and sequencing, four haplotypes were found and three standard strains of Leishmania from other places (WHO1 (MHOM/SU/73/5ASKH), K27.11 (MHOM/SU/74/K27) from Soviet Union and FR05 (MHOM/IL/80/fredlin) from Italy) were sampled with the DNA extracted for precise recognition of the Leishmania parasites using the ITS-rDNA and Cyt b genes (Supplementary Figure 1), (Table 1).

TABLE 1:
Leishmania species identification using PCR applying ITS-rDNA and Cyt b genes and separation by HRM-qPCR.

Clinical specimens were isolated and smeared according to the Declaration of Helsinki, Ethical Principles for Medical Research Involving Human Subjects protocol, and Ethical Committee Guidelines of Pasteur Institute of Iran1515. World Health Organization (WHO). Manual of basic techniques for a health laboratory. 2nd ed. WHO Publications Centre: USA, Geneva; 2003.,1616. WMA, Declaration of Helsinki, ethical principles for medical research involving human subjects. 2013.. The fixation of smears was followed by Giemsa staining, so that amastigotes were seen using light microscope (×1000).

Culturing samples in Novy-MacNeal-Nicolle (NNN) medium were carried out, then incubation was done and samples were checked for the observation of the presence and growth of promastigotes. The isolation of promastigotes was followed by transferring to RPMI-1640 media and the counting of the parasites in a Neubauer chamber (about 10 million promastigotes per/ml).

DNA was extracted using a tissue extraction kit (Takapozist, Iran). The promastigotes were washed using phosphate buffer after one week of culturing. Proteinase K (20 μl) was added to the parasites (200 μl). Samples were subjected to 56 ºC hot bath. Absolute ethanol was added. Columns were washed with releasing buffer to obtain purified extracted DNA.

Designing new primers and DNA amplification of Leishmania parasites

Cyt b is mitochondrial; it is an enzymatic or expression gene. It has a 20-50 copy number. The infection was detected using forward primer (LCBF1F) (5' GGTGTAGGTTTTAGTTTAGG 3') and reverse primer (RLCBR2)(5'CTACAATAAACAAATCATAATATACAATT 3')1717. Luyo-Acero G E, Uezato H, Oshiro M, Takei K, Kariya K, Katakura K, et al. Sequence variation of the cytochrome b gene of various human infecting members of the genus Leishmania and their phylogeny. Parasitology. 2004;128(Pt 5):483-91.. After cloning the fragment, two species were differentiated with new designed primers using the CLC DNA workbench 5.5 software (CLC, Bio-QIAGEN, Aarhus, Denmark). For the newly designed primers, following criteria were considered: 1- the primer length has to be between 11-30 bp; 2- the hybridization temperature of the primers has to be between 58-60 °C; 3- the candidate region has to be specific for Leishmania genus; 4- the prediction of the temperature regarding the mutations has to be intended for detecting region. Then newly designed primers were used in HRM-qPCR: forward primer (5' ACTCACGGCCTCTAGGAATGA3') and reverse primer (5' AGAATTTCACCTCTGACGCCCCAGT3').

ITS-rDNA is a chromosomal, housekeeping gene. It has 40-200 copy numbers. Using forward primer (ITS1F) (5’ GCAGCTGGATCATTTTCC 3’) and reverse primer ITS2R4 (5’ ATATGCAGAAGAGAGAGGAGGC 3’), 480 bp fragment was amplified and the infection was detected1818. Parvizi P. First detection of Leishmania major in peridomestic Phlebotomus papatasi from Isfahan province, Iran: comparison of nested PCR of nuclear ITS ribosomal DNA and semi-nested PCR of minicircle kinetoplast DNA. Acta Trop. 2005;93(1):75-83.. Then another set of primers was used in HRM-qPCR: (ITSmF1) (5’ GTGTGGAAGCCAAGAGGAGG 3’) and reverse primer (ITSmR2) (5’ GCAAGCACCAGAGAGGAGT 3’); the size of the fragment (microsatellites of ITS-rDNA gene) was 160 bp, including primers1919. Parvizi P, Ready PD. Nested PCRs and sequencing of nuclear ITS-rDNA fragments detect three Leishmania species of gerbils in sandflies from Iranian foci of zoonotic cutaneous leishmaniasis. Trop Med Int Health. 2008;13(9):1159-71..

For these newly- designed primers, the similarity, sensitivity, specificity, proximity, and coverage percentages were evaluated with L. major, L. Infantum, L. turanica, L. tropica, L. donovani, and L. gerbili (Figure 1). For each PCR reaction, a total volume of 25 μl was utilized containing 0.4 pmol of each primer, 1X PCR buffer (750 mM of Tris-HCl, pH: 8.8 at 25 °C, 200 mM of (NH4)2SO4 and 0.1% Tween-20), 2.5 mM of MgCl2, 0.5 mM of deoxy ribonucleotide triphosphate (dNTP), 1 U of Taq polymerase (CinnaClon Co., Iran), and 10 ng/reaction of each DNA. For the newly designed primers, PCR program was as follows: initial denaturation at 94 °C for 4 min, followed by 35 cycles at 94 °C for 30 s, annealing for each gene (Cyt b: at 50 °C for 60 s; ITS-rDNA: at 58 °C for 45 s), and extension at 72 °C for 1 min, as well as final extension at 72 °C for 5 min. After purification, PCR products were sent for sequencing by an ABI PRISM™ 310 automated sequencer (Applied Biosystems, USA), using reverse and forward primers. Sequences were edited and aligned via Sequencher 3.1.1 and the CLC Main Workbench 5.5 software (CLC, Bio-QIAGEN, Aarhus, Denmark).

FIGURE 1:
Comparison among cover percentage of studied species in this research and other species from GenBank using the CLC Main Workbench 5.5 software.

ITS-rDNA and Cyt b genes amplification by Real- time PCR

Cyt b was amplified applying newly designed primers using Real-time PCR. ITS-rDNA was amplified according to the protocols of Parvizi and Ready, 2008. The mapping genes of these primers and time-temperature schedule were determined using the CLC DNA workbench 5.5 software (CLC, Bio-QIAGEN, Aarhus, Denmark). The protocols were modified and performed to detect Leishmania species employing the HRM method by Rotor-Gene 6000 Series Software 1.7. The modified protocol of HRM was followed: first, holding 95 °C for 15 min by using 4 µl of 5x HOT FIREPol® EvaGreen® qPCR Mix Plus (Solis Bio Dyne, Estonia). In the first stage, DNA denaturation of template and activation of polymerase enzyme were carried out at 95 °C for 15 s. In the second stage, DNA replication was performed in 45 cycles, at 55 °C for 30 s and at 72 °C for 30 s, with final extension cycle at 72 °C for 300 s. The separation curve or melting curve were generated at 95 °C for 60 s or 60 °C for 60 s. Finally, to obtain HRM from 65 to 95 °C, with every 0.1 °C rise, two seconds were paused and finally all of these steps were optimized. GC content, length, sequence, and heterozygosity of DNA affect the melting characteristics of a PCR product. PCR products were controlled by saturated EvaGreen color and were fluoresced in the presence of two-strand of DNA.

The standard samples for qPCR and re-sequencing of ambiguous sites were prepared according to Ghafari et al, 2020 with some modifications2020. Ghafari SM, Fotouhi-Ardakani R, Parvizi P. Designing and developing a high- resolution melting technique for accurate identification of Leishmania species by targeting amino acid permease 3 and cytochrome oxidase II genes using real-time PCR and in silico genetic evaluation. Acta Trop. 2020;211:105626..

Absolute quantification qPCR and improving HRM

Leishmania parasites were identified using newly designed primers by high-resolution melting-curve analysis. Amplification of Leishmania DNA for real-time PCR was set for a total volume of 20 μl containing 0.7 µl of each designed forward/reveres primers, EvaGreen 5x (Yekta tajhiz azma), 12.6 µl sterile DNase/RNase-free water, and 10 ng/reaction of each DNA sample. Functionality and sensitivity of HRM assay were developed to distinguish and recognize the Leishmania species. The results of HRM dissociation were compared with the standard PCR data. Melting temperature (Tm) values for each sample were calculated by independent standards from GenBank using the CLC Main Workbench.v5.5 software (CLC, Bio-QIAGEN, Aarhus, Denmark).

RESULTS

Leishmania identification by employing ITS-rDNA and Cyt b genes

Two Leishmania species from Iranian sources and some Leishmania reference strains were precisely identified by employing ITS-rDNA and Cyt b genes (Table 1).

The sensitivity and specificity of each gene were identified and evaluated to accurately identify Leishmania species (L. major and L. tropica).

Tm and Cycle of Threshold (CT) of ITS-rDNA and Cyt b genes of L. major and L. tropica from Iran and some regions of the world were analyzed and evaluated. According to the results, Cyt b gene was more sensitive than ITS-rDNA. Tm for Cyt b gene was lower than ITS-rDNA and could identify the species DNA. The length of ITS-rDNA fragment was 160 bp, meanwhile, that of Cyt b fragment was 866 bp.

For both Leishmania species, amplifying, genotyping, and sequencing of amplified product were performed using newly designed primers. HRM assay managed to detect both species of Leishmania and could determine similar Leishmania species in the New and Old World. Each melting profile (keeping the same chemical environment with similar reagents and DNA concentrations) was considered for consistency examination. Six consecutive dilutions of DNA template and sequencing dilution reactions were prepared to optimize real-time PCR assay. Duplicates of consecutive dilutions with a DNA-free reaction template (negative control) were prepared to find the specific temperature for separating each Leishmania species using Cyt b. At the end of the reaction, the standard curve was plotted. Temperature was drawn on a horizontal axis and CT on a vertical axis. After proliferation of DNA dilutions, L. major and L. tropica were traced to the standard curve (Figure 2, Table 2).

FIGURE 2:
Derivative melt curve analysis of ITS2 gene in different species of isolated Leishmania parasites using the HRM assay.

TABLE 2:
Serial dilution of plasmid copy numbers (106-101) of two Leishmania parasites using Cyt b to find specific temperature for separating of each Leishmania species.

The melting characteristics of amplicons were assessed for Leishmania species by plotting two different curves (Figure 3A,B). In this research, aligned melt curve, the normalized fluorescence curve, and difference plot melt curve were produced to easily distinguish each Leishmania species. Tm for L. major and L. tropica was 84.7 °C and 83.6 °C, respectively, which were separated by less than 0.9 °C of temperature difference.

FIGURE 3:
Normalized melting profiles obtained with the HRM assays for Leishmania major and Leishmania tropica. Each species was tested in triplicate. A.1: normalized melting curves for both species; A.2: a difference plot for both species.

The sensitivity of HRM based on serial dilution and drawing standard curve

Different dilutions of DNA template were prepared for validation, accuracy, and sensitivity of HRM assay. High sensitivity of HRM could withdraw minimum amounts of DNA to draw appropriate melting curve. Serial dilutions of Leishmania were experienced on three consecutive days with similar reactions. Identical DNA concentrations were repeated with the Cyt b gene. The melting curves for serial dilutions were compared. Identical DNAs were analyzed and melt patterns were appropriated for both species in different days (Figure 2).

CONCLUSIONS

Truthful, quick, and delicate diagnostic procedures are essential for detecting and characterizing causative CL for accurate treatment, precise prognosis, and understanding epidemiology of various species of Leishmania parasites1515. World Health Organization (WHO). Manual of basic techniques for a health laboratory. 2nd ed. WHO Publications Centre: USA, Geneva; 2003.,2121. Rojas-Jaimes J, Rojas-Palomino N, Pence J, Lescano AG. Leishmania species in biopsies of patients with different clinical manifestations identified by high resolution melting and nested PCR in an endemic district in Peru. Parasite Epidemiol Control. 2019;3:e00095.,2222. Safeukui I, Millet P, Boucher S, Melinard L, Fregeville F, Receveur MC, et al. Evaluation of FRET real-time PCR assay for rapid detection and differentiation of Plasmodium species in returning travelers and migrants. Malar J. 2008;7:70..

Phylogenetic and evolutionary analysis should investigate different clinical lesions for avoiding misdiagnosis of leishmaniasis, and helping investigators with the development of molecular epidemiological investigations that can be appropriate for extirpating Leishmania parasites in diverse regions2323. Asfaram S, Fakhar M, Mirani N, Derakhshani-Niya M, Valadan R, Ziaei Hezarjaribi H, et al. HRM-PCR is an accurate and sensitive technique for the diagnosis of cutaneous leishmaniasis as compared with conventional PCR. Acta Parasitol. 2020;65(2):310-6. Available from: https://doi.org/10.2478/s11686-019-00154-5.
https://doi.org/10.2478/s11686-019-00154...
,2424. Talmi-Frank D, Nasereddin A, Schnur LF, Schönian G, Ozensoy Töz S, Jaffe CL, et al. Detection and identification of Old World Leishmania by high resolution melt analysis. PLoS Negl Trop Dis . 2010;4(1):e581..

Thus, designing of primers for distinguishing causative agents of leishmaniasis is necessary to find an application with high sensitivity, efficiency, and accuracy22. Fotouhi-Ardakani R, Dabiri S, Ajdari S, Alimohammadian MH, AlaeeNovin E, Taleshi N, et al. Assessment of nuclear and mitochondrial genes in precise identification and analysis of genetic polymorphisms for the evaluation of Leishmania parasites. Infect Genet Evol. 2016;46:33-41.,55. Zahirnia AH, Bordbar A, Ebrahimi S, Spotin A, Mohammadi S, Ghafari SM, et al. Predominance of Leishmania major and rare occurrence of Leishmania tropica with haplotype variability at the center of Iran. Braz J Infect Dis. 2018;22(4):278-87.,1919. Parvizi P, Ready PD. Nested PCRs and sequencing of nuclear ITS-rDNA fragments detect three Leishmania species of gerbils in sandflies from Iranian foci of zoonotic cutaneous leishmaniasis. Trop Med Int Health. 2008;13(9):1159-71.. In healthcare and research centers of endemic foci of the disease, which have a high number of patients and patrons, accuracy and speed are of great importance77. Bates PA, Depaquit J, Galati EA, Kamhawi SH, Maroli M, Ann McDowellet M, et al. Recent advances in phlebotomine sand fly research related to leishmaniasis control. Parasit Vectors . 2015;8(131):11-21.,1111. Kaltenboeck B, Wang C. Advances in real-time PCR: application to clinical laboratory diagnostics. Adv Clin Chem. 2005;40:219-59.,2424. Talmi-Frank D, Nasereddin A, Schnur LF, Schönian G, Ozensoy Töz S, Jaffe CL, et al. Detection and identification of Old World Leishmania by high resolution melt analysis. PLoS Negl Trop Dis . 2010;4(1):e581.. Isoenzyme analyses was a gold standard method for typifying Leishmania species, but recently, PCR, RFLP, sequencing, and phylogenetic analysis, and now, HRM method using qPCR are more acceptable. No distinct extensively accepted method is available to identify Leishmania species and strains today. A valuable molecular technique is needed to be developed to separate Leishmania parasites in primary infections11. Fotouhi-Ardakani R, Ghafari SM, Ready PD, Parvizi P. Developing, modifying, and validating a TaqMan Real-Time PCR technique for accurate identification of Leishmania parasites causing most leishmaniasis in Iran. Front Cell Infect Microbiol. 2021;11:e731595. Available from: https://doi.org/10.3389/fcimb.2021.731595
https://doi.org/10.3389/fcimb.2021.73159...
,44. Zampieri RA, Laranjeira-Silva MF, Muxel SM, Stocco de Lima AC, Jon Shaw J, Floeter-Winter LM. High Resolution Melting analysis targeting hsp70 as a fast and eficient method for the discrimination of Leishmania species. PLoS Negl Trop Dis. 2016;10(2): e0004485.,77. Bates PA, Depaquit J, Galati EA, Kamhawi SH, Maroli M, Ann McDowellet M, et al. Recent advances in phlebotomine sand fly research related to leishmaniasis control. Parasit Vectors . 2015;8(131):11-21..

Real-time PCR (qPCR) is convenient, fast, simple, and cost-effective and is widely used in researches for accurate identification and quantification of pathogens. qPCR with high accuracy, sensitivity, specificity, reliability, and availability has been employed for gene scanning, genotyping, and recognition of dissimilar species of Leishmania2525. Taylor SC, Nadeau K, Abbasi M, Lachance C, Nguyen M, Fenrich J. The ultimate qPCR experiment: producing publication quality, reproducible data the first time. Trends Biotechnol Res. 2019;37(7):761-74.,2626. Antonia AL, Wang L, Ko DC. A real-time PCR assay for quantification of parasite burden in murine models of leishmaniasis. Peer Community J. 2018;6:e5905.. HRM technique has been recorded in recent years to recognize and classify parasitic diseases including leishmaniasis2121. Rojas-Jaimes J, Rojas-Palomino N, Pence J, Lescano AG. Leishmania species in biopsies of patients with different clinical manifestations identified by high resolution melting and nested PCR in an endemic district in Peru. Parasite Epidemiol Control. 2019;3:e00095.,2323. Asfaram S, Fakhar M, Mirani N, Derakhshani-Niya M, Valadan R, Ziaei Hezarjaribi H, et al. HRM-PCR is an accurate and sensitive technique for the diagnosis of cutaneous leishmaniasis as compared with conventional PCR. Acta Parasitol. 2020;65(2):310-6. Available from: https://doi.org/10.2478/s11686-019-00154-5.
https://doi.org/10.2478/s11686-019-00154...
,2727. Hernandez C, Alvarez C, Gonzalez C, Stella Ayala M, Maritza León C, David Ramírez J. Identification of six New World Leishmania species through the implementation of a High-Resolution Melting (HRM) genotyping assay. Parasit Vectors . 2014;7:501-8.. HRM can determine unique curves and detect single nucleotide polymorphisms according to a few differences in the nucleotide content of Leishmania parasites2828. Ahuja K, Vats A, Beg MA, Kariyawasam K, Chaudhury A, Chatterjee M, et al. High resolution melting based method for rapid discriminatory diagnosis of co-infecting Leptomonas seymouri in Leishmania donovani induced leishmaniasis. Parasitol Int. 2020;75:102047..

Also, previously New World Leishmania have been identified and analyzed by qPCR and HRM2929. Galluzzi L, Ceccarelli M, Diotallevi A, Menotta M, Magnani M. Real-time PCR applications for diagnosis of leishmaniasis. Parasit Vectors . 2018;11(1):273.. HRM is a new approach to real-time PCR; working with HRM requires experience and skills to interpret the results33. Müller KE, Zampieri RA, Aoki J, Muxel SM, Nerland AH, Floeter-Winter LM. Amino acid permease 3 (AAP3) coding sequence as a target for Leishmania identification and diagnosis of leishmaniases using high resolution melting analysis. Parasit Vectors. 2018;11:421.. The third generation EvaGreen is a double-stranded DNA-attached color, and one of its characteristics is the ability to inhibit PCR in comparison with SYBR green I. EvaGreen is also suitable for HRM3030. Monis PT, Giglio S, Saint CP. Comparison of SYTO9 and SYBR Green I for real301 timepolymerase chain reaction and investigation of the effect of dye concentration on amplification and DNA melting curve analysis. Anal Biochem. 2005;340(1):24-34.,3131. Mao F, Leung WY, Xin X. Characterization of EvaGreen and the implication of its physicochemical properties for qPCR applications. BMC Biotechnol. 2007;7:76..

HRM was employed to identify the L. major, L. tropica, L. donovani, and L. infantum species by targeting kDNA minicircle, ITS-rRNA, SSU rDNA, g6pd, and hsp70 genes; these have been identified in Eurasia and South America33. Müller KE, Zampieri RA, Aoki J, Muxel SM, Nerland AH, Floeter-Winter LM. Amino acid permease 3 (AAP3) coding sequence as a target for Leishmania identification and diagnosis of leishmaniases using high resolution melting analysis. Parasit Vectors. 2018;11:421.,2424. Talmi-Frank D, Nasereddin A, Schnur LF, Schönian G, Ozensoy Töz S, Jaffe CL, et al. Detection and identification of Old World Leishmania by high resolution melt analysis. PLoS Negl Trop Dis . 2010;4(1):e581.. By targeting parasite ITS-rDNA and Cyt b genes, the most famous Leishmania species were confidently identified using modified and developed HRM as an efficient and cost-effective method22. Fotouhi-Ardakani R, Dabiri S, Ajdari S, Alimohammadian MH, AlaeeNovin E, Taleshi N, et al. Assessment of nuclear and mitochondrial genes in precise identification and analysis of genetic polymorphisms for the evaluation of Leishmania parasites. Infect Genet Evol. 2016;46:33-41.,88. Odell ID, Cloud JL, Seipp M, Wittwer CT. Rapid species identification within the Mycobacterium chelonae-abscessus group by high-resolution melting analysis of hsp65 PCR products. Am J Clin Pathol. 2005;123(1):96-101..

In this research, Cyt b and ITS-rDNA genes were targeted, and HRM technique was modified to be more trustworthy, strong, slight, and obvious, to allow the correct genotyping of Leishmania species2222. Safeukui I, Millet P, Boucher S, Melinard L, Fregeville F, Receveur MC, et al. Evaluation of FRET real-time PCR assay for rapid detection and differentiation of Plasmodium species in returning travelers and migrants. Malar J. 2008;7:70.,3232. Chelbi H, Essid R, Jelassi R, Bouzekri N, Zidi I, Ben Salah H, et al. High-resolution melting-curve (HRM) analysis for C. meleagridis identification in stool samples. Microb Pathog. 2018;115:332-7.. In this method, new design of dedicated primers has been one of the most important parts of the work. New designed primers (more economical matter and efficiency with high sensitivity in patients) were considered as advantages of this method. Melting curves provide different temperatures with new primers designed to isolate Leishmania species and strains using the ITS-rDNA gene with short repeat (STR) microsatellite employing the HRM method. The prediction results were corrected using CLC software. The altered regions could establish a desirable temperature difference in species and strains of Leishmania parasites1919. Parvizi P, Ready PD. Nested PCRs and sequencing of nuclear ITS-rDNA fragments detect three Leishmania species of gerbils in sandflies from Iranian foci of zoonotic cutaneous leishmaniasis. Trop Med Int Health. 2008;13(9):1159-71.,2727. Hernandez C, Alvarez C, Gonzalez C, Stella Ayala M, Maritza León C, David Ramírez J. Identification of six New World Leishmania species through the implementation of a High-Resolution Melting (HRM) genotyping assay. Parasit Vectors . 2014;7:501-8..

Both genes were sequenced to recognize L. major and L. tropica and had clear sequences of readings. If a mutation occurs in the housekeeping genes (ITS-rDNA), it leads to create a new species of Leishmania. These genes are in charge of conserving primary cellular function and their expression is crucial for cell survival3333. Peacock CS, Seeger K, Harris D, Murphy L, Ruiz JC, A Quail M, et al. Comparative genomic analysis of three Leishmania species that cause diverse human disease. Nat Genet. 2007;39(7):839-47.. The sequences of our fragment were 160 bp in the ITS-rDNA gene, which were analyzed and evaluated against control samples.

To validate ITS-rDNA gene, the results were compared with sequences of NCBI database showing L. major (accession number: HQ591339) with a similarity of 96% and L. tropica (accession number: KT 972151.1) with a similarity of 86.59%. The sensitivity and specificity of primers were designed by BLAST program whereby ITS-rDNA sequences revealed similarity of about 100% with both L. major and L. tropica. Melting curve analysis was also used to distinguish these parasites2222. Safeukui I, Millet P, Boucher S, Melinard L, Fregeville F, Receveur MC, et al. Evaluation of FRET real-time PCR assay for rapid detection and differentiation of Plasmodium species in returning travelers and migrants. Malar J. 2008;7:70.,3232. Chelbi H, Essid R, Jelassi R, Bouzekri N, Zidi I, Ben Salah H, et al. High-resolution melting-curve (HRM) analysis for C. meleagridis identification in stool samples. Microb Pathog. 2018;115:332-7.. In expression Cyt b, evolution occurred while no evolution was observed in ITS-rDNA.

Accordingly, we can conclude that HRM assay has many advantages over other genotype methods, including low costs as well as homogeniety, accuracy, and fast speed in a closed tube. Due to its ease of use and rapidness, it is suitable for detection as a quick and low-cost method. HRM is a relatively new technique that provides a clear feature of PCR amplicons. In this method, the changes in fluorescence light that are emitted to DNA at conversion time of a two-stranded DNA into a single-stranded DNA, are measured and can be distinguished from single-nucleotide polymorphism (SNP)3131. Mao F, Leung WY, Xin X. Characterization of EvaGreen and the implication of its physicochemical properties for qPCR applications. BMC Biotechnol. 2007;7:76..

Our experimental results revealed that the HRM method using EvaGreen and applying ITS-rDNA and Cyt b genes is not only definite and slight, but also repeatable and accurate, and can be intended as a technique for the early recognition and approval of a Leishmania infection.

ACKNOWLEDGMENTS

The authors thank staff of Molecular Systematics Laboratory, Parasitology Department, Pasteur Institute of Iran for their assistance.

REFERENCES

  • 1
    Fotouhi-Ardakani R, Ghafari SM, Ready PD, Parvizi P. Developing, modifying, and validating a TaqMan Real-Time PCR technique for accurate identification of Leishmania parasites causing most leishmaniasis in Iran. Front Cell Infect Microbiol. 2021;11:e731595. Available from: https://doi.org/10.3389/fcimb.2021.731595
    » https://doi.org/10.3389/fcimb.2021.731595
  • 2
    Fotouhi-Ardakani R, Dabiri S, Ajdari S, Alimohammadian MH, AlaeeNovin E, Taleshi N, et al. Assessment of nuclear and mitochondrial genes in precise identification and analysis of genetic polymorphisms for the evaluation of Leishmania parasites. Infect Genet Evol. 2016;46:33-41.
  • 3
    Müller KE, Zampieri RA, Aoki J, Muxel SM, Nerland AH, Floeter-Winter LM. Amino acid permease 3 (AAP3) coding sequence as a target for Leishmania identification and diagnosis of leishmaniases using high resolution melting analysis. Parasit Vectors. 2018;11:421.
  • 4
    Zampieri RA, Laranjeira-Silva MF, Muxel SM, Stocco de Lima AC, Jon Shaw J, Floeter-Winter LM. High Resolution Melting analysis targeting hsp70 as a fast and eficient method for the discrimination of Leishmania species. PLoS Negl Trop Dis. 2016;10(2): e0004485.
  • 5
    Zahirnia AH, Bordbar A, Ebrahimi S, Spotin A, Mohammadi S, Ghafari SM, et al. Predominance of Leishmania major and rare occurrence of Leishmania tropica with haplotype variability at the center of Iran. Braz J Infect Dis. 2018;22(4):278-87.
  • 6
    Ready PD. Biology of phlebotomine sand flies as vectors of disease agents. Annu Rev Entomol. 2013;58:227-50.
  • 7
    Bates PA, Depaquit J, Galati EA, Kamhawi SH, Maroli M, Ann McDowellet M, et al. Recent advances in phlebotomine sand fly research related to leishmaniasis control. Parasit Vectors . 2015;8(131):11-21.
  • 8
    Odell ID, Cloud JL, Seipp M, Wittwer CT. Rapid species identification within the Mycobacterium chelonae-abscessus group by high-resolution melting analysis of hsp65 PCR products. Am J Clin Pathol. 2005;123(1):96-101.
  • 9
    Pangasa A, Jex AR, Campbell BE, Bott NJ, Whipp M, Hogg G, et al. High resolution melting-curve (HRM) analysis for the diagnosis of cryptosporidiosis in humans. Mol Cell Probes. 2009;23(1):10-5.
  • 10
    Chua KH, Ng CC, Lee PC, Lee PC, Lian Lim YA, Lau TP, et al. Development of high resolution melting analysis for the diagnosis of human malaria. Cell Rep Phys Sci. 2015;5:e15671.
  • 11
    Kaltenboeck B, Wang C. Advances in real-time PCR: application to clinical laboratory diagnostics. Adv Clin Chem. 2005;40:219-59.
  • 12
    Wang W, Zijlstra RT, Gänzle MG. Identification and quantification of virulence factors of enterotoxigenic Escherichia coli by high-resolution melting curve quantitative PCR. BMC Microbiol. 2017;17(1):114.
  • 13
    Ahmad AI, Ghasemi JB. New unsymmetrical cyanine dyes for real-time thermal cycling. Anal Bioanal Chem. 2007;389(3):983-8.
  • 14
    Eischeid AC. SYTO dyes and EvaGreen outperform SYBR Green in real-time PCR. BMC Res Notes. 2011;4:263.
  • 15
    World Health Organization (WHO). Manual of basic techniques for a health laboratory. 2nd ed. WHO Publications Centre: USA, Geneva; 2003.
  • 16
    WMA, Declaration of Helsinki, ethical principles for medical research involving human subjects. 2013.
  • 17
    Luyo-Acero G E, Uezato H, Oshiro M, Takei K, Kariya K, Katakura K, et al. Sequence variation of the cytochrome b gene of various human infecting members of the genus Leishmania and their phylogeny. Parasitology. 2004;128(Pt 5):483-91.
  • 18
    Parvizi P. First detection of Leishmania major in peridomestic Phlebotomus papatasi from Isfahan province, Iran: comparison of nested PCR of nuclear ITS ribosomal DNA and semi-nested PCR of minicircle kinetoplast DNA. Acta Trop. 2005;93(1):75-83.
  • 19
    Parvizi P, Ready PD. Nested PCRs and sequencing of nuclear ITS-rDNA fragments detect three Leishmania species of gerbils in sandflies from Iranian foci of zoonotic cutaneous leishmaniasis. Trop Med Int Health. 2008;13(9):1159-71.
  • 20
    Ghafari SM, Fotouhi-Ardakani R, Parvizi P. Designing and developing a high- resolution melting technique for accurate identification of Leishmania species by targeting amino acid permease 3 and cytochrome oxidase II genes using real-time PCR and in silico genetic evaluation. Acta Trop. 2020;211:105626.
  • 21
    Rojas-Jaimes J, Rojas-Palomino N, Pence J, Lescano AG. Leishmania species in biopsies of patients with different clinical manifestations identified by high resolution melting and nested PCR in an endemic district in Peru. Parasite Epidemiol Control. 2019;3:e00095.
  • 22
    Safeukui I, Millet P, Boucher S, Melinard L, Fregeville F, Receveur MC, et al. Evaluation of FRET real-time PCR assay for rapid detection and differentiation of Plasmodium species in returning travelers and migrants. Malar J. 2008;7:70.
  • 23
    Asfaram S, Fakhar M, Mirani N, Derakhshani-Niya M, Valadan R, Ziaei Hezarjaribi H, et al. HRM-PCR is an accurate and sensitive technique for the diagnosis of cutaneous leishmaniasis as compared with conventional PCR. Acta Parasitol. 2020;65(2):310-6. Available from: https://doi.org/10.2478/s11686-019-00154-5.
    » https://doi.org/10.2478/s11686-019-00154-5
  • 24
    Talmi-Frank D, Nasereddin A, Schnur LF, Schönian G, Ozensoy Töz S, Jaffe CL, et al. Detection and identification of Old World Leishmania by high resolution melt analysis. PLoS Negl Trop Dis . 2010;4(1):e581.
  • 25
    Taylor SC, Nadeau K, Abbasi M, Lachance C, Nguyen M, Fenrich J. The ultimate qPCR experiment: producing publication quality, reproducible data the first time. Trends Biotechnol Res. 2019;37(7):761-74.
  • 26
    Antonia AL, Wang L, Ko DC. A real-time PCR assay for quantification of parasite burden in murine models of leishmaniasis. Peer Community J. 2018;6:e5905.
  • 27
    Hernandez C, Alvarez C, Gonzalez C, Stella Ayala M, Maritza León C, David Ramírez J. Identification of six New World Leishmania species through the implementation of a High-Resolution Melting (HRM) genotyping assay. Parasit Vectors . 2014;7:501-8.
  • 28
    Ahuja K, Vats A, Beg MA, Kariyawasam K, Chaudhury A, Chatterjee M, et al. High resolution melting based method for rapid discriminatory diagnosis of co-infecting Leptomonas seymouri in Leishmania donovani induced leishmaniasis. Parasitol Int. 2020;75:102047.
  • 29
    Galluzzi L, Ceccarelli M, Diotallevi A, Menotta M, Magnani M. Real-time PCR applications for diagnosis of leishmaniasis. Parasit Vectors . 2018;11(1):273.
  • 30
    Monis PT, Giglio S, Saint CP. Comparison of SYTO9 and SYBR Green I for real301 timepolymerase chain reaction and investigation of the effect of dye concentration on amplification and DNA melting curve analysis. Anal Biochem. 2005;340(1):24-34.
  • 31
    Mao F, Leung WY, Xin X. Characterization of EvaGreen and the implication of its physicochemical properties for qPCR applications. BMC Biotechnol. 2007;7:76.
  • 32
    Chelbi H, Essid R, Jelassi R, Bouzekri N, Zidi I, Ben Salah H, et al. High-resolution melting-curve (HRM) analysis for C. meleagridis identification in stool samples. Microb Pathog. 2018;115:332-7.
  • 33
    Peacock CS, Seeger K, Harris D, Murphy L, Ruiz JC, A Quail M, et al. Comparative genomic analysis of three Leishmania species that cause diverse human disease. Nat Genet. 2007;39(7):839-47.
  • Financial Support: This work was supported by Pasteur Institute of Iran (grant numbers 973) and the Iran National Science Foundation (INSF) (grant number 95840608).

SUPPLEMENTARY MATERIAL

SUPPLEMENTARY FIGURE 1:
Graphical abstract highlighting important methodical aspects of this research.

Publication Dates

  • Publication in this collection
    16 Dec 2022
  • Date of issue
    2022

History

  • Received
    12 Apr 2022
  • Accepted
    24 Aug 2022
Sociedade Brasileira de Medicina Tropical - SBMT Caixa Postal 118, 38001-970 Uberaba MG Brazil, Tel.: +55 34 3318-5255 / +55 34 3318-5636/ +55 34 3318-5287, http://rsbmt.org.br/ - Uberaba - MG - Brazil
E-mail: rsbmt@uftm.edu.br