- Open Access
Exploiting genetic polymorphisms in metabolic enzymes for rapid screening of Leishmania infantum genotypes
Parasites & Vectorsvolume 11, Article number: 572 (2018)
Leishmania infantum is the aetiological agent of visceral leishmaniasis (VL) and cutaneous leishmaniasis (CL). Numerous strains and/or zymodemes have been identified and characterized by multilocus enzyme electrophoresis (MLEE). MLEE is considered the reference method for L. infantum parasite typing and it is based upon enzyme electrophoretic mobility analysis from promastigote cultures. However, the MLEE technique is cumbersome, time-consuming and does not detect silent genetic mutations or nucleotide changes that give rise to amino acid changes that do not alter electrophoretic mobility. As a result of these difficulties, many DNA-based typing methods have been developed over the past few years. However, relative to the enzymes utilized in MLEE analysis, we observed a shortage of DNA sequences available in the GenBank database or an absolute lack of sequences belonging to specific zymodemes. The aims of the present study were to (i) implement the number of sequences coding for metabolic enzymes used in MLEE; (ii) identify polymorphisms that characterize L. infantum zymodemes most prevalent in the Mediterranean basin; and (iii) exploit these polymorphisms to develop a rapid screening test that would give results comparable with existing MLEE typing.
Partial sequences of seven metabolic enzyme genes (malic enzyme, 6-phosphogluconate dehydrogenase, mitochondrial isocitrate dehydrogenase, glucose-6-phosphate isomerase, glucose-6-phosphate dehydrogenase, phosphoglucomutase and mannose phosphate isomerase) were obtained from 11 L. infantum strains. The comparison of these sequences with those obtained from GenBank allowed for the identification of a few polymorphisms that could distinguish several zymodemes. In particular, the polymorphism 390T>G in the malic enzyme gene has been exploited to develop a high-resolution melt (HRM)-based assay to rapidly differentiate the genotype 390T, associated with zymodemes MON-1, MON-72 and MON-201, evidencing a partial agreement between genotyping results and MLEE. The assay has been successfully applied to L. infantum clinical isolates and clinical samples.
A HRM-based assay for rapid identification of genotypes associated with the most common L. infantum zymodemes in the Mediterranean basin has been developed and its potential application in epidemiological research for L. infantum population screening, without parasite isolation and culturing, has been demonstrated.
The taxonomy of the genus Leishmania is complex and has been revised several times in light of parasite biology and biochemistry [1, 2]. Currently, the genus Leishmania includes four subgenera: Leishmania, Viannia, Sauroleishmania and L. enriettii complex. Each subgenus includes several species . Leishmania infantum, the aetiological agent of visceral leishmaniasis (VL) and cutaneous leishmaniasis (CL), belongs to subgenus Leishmania .
Multilocus enzyme electrophoresis (MLEE) is the isoenzymatic analysis presently considered the reference method for parasite typing by the World Health Organization . This technique, developed at the Centre for Leishmaniasis of Montpellier (France), is based on the electrophoretic mobility of several enzymes obtained from promastigote culture, as described by Rioux et al. . Up to 15 enzymes can be considered in this analysis: malate dehydrogenase, malic enzyme, isocitrate dehydrogenase (NADP), phosphogluconate dehydrogenase, glucose-6-phosphate dehydrogenase, glutamate dehydrogenase, cytochrome-b5 reductase, purine-nucleoside phosphorylase 1 and 2, glutamic-oxaloacetic transaminase 1 and 2, phosphoglucomutase, fumarase, mannose-6-phosphate isomerase, and glucose-6-phosphate isomerase . The comparison of isoenzyme mobility with a reference strain identified over 300 zymodemes, also termed MON.
Concerning L. infantum, 45 zymodemes have been found, 36 recognized as real zymodemes and 9 as variants . In dogs in the Mediterranean area, 12 zymodemes have been found: MON-1, MON-24, MON-34, MON-72, MON-77, MON-80, MON-98, MON-105, MON-108, MON-199, MON-199 variant NP1130 and MON-281. With the exception of MON-105, these zymodemes have also all been found in humans . Other zymodemes, such as MON-11, MON-27, MON-28, MON-29 MON-33 and MON-189, have been isolated only from humans . Leishmania infantum MON-1 is the most frequent zymodeme in both humans and dogs in the Mediterranean basin  and in the New World . It is present in 30 countries and represents about 70% of all identified strains. In particular, MON-1 represents 88% of L. infantum zymodemes in the south of France , 96.7% in Portugal  and 45–58% in Spain [12, 13]. Moreover, in HIV-infected subjects, MON-1 was found in about 73% of co-infections, while in immunocompetent patients it was found in about 90% of VL cases and 20% of CL cases . Concerning canine infections in Italy, Gramiccia et al.  showed a high prevalence of MON-1 zymodeme (91%) in 497 canine clinical isolates, with the remaining percentage composed almost exclusively of MON-72 zymodeme. The homogeneity of the isoenzymes identified in the dog population does not fully reflect the isoenzyme diversity in humans; therefore, the role of canine population as the only reservoir for all L. infantum zymodemes remains unclear. In light of this epidemiological evidence, it is advisable that a more extensive typing of L. infantum be carried out. However, MLEE presents several disadvantages: it is time-consuming, technically demanding, and requires bulk cultures of parasites. Moreover, due to its inability to detect nucleotide substitutions that do not change the amino acid composition, its discriminatory power is poor. Therefore, molecular approaches such as multilocus microsatellite typing (MLMT) [16, 17] and multilocus sequence typing (MLST) [18,19,20] have been introduced for strain discrimination. MLST results were generally in agreement with MLEE, though increased resolution was obtained and some key discrepancies were found [18,19,20,21]. For example, L. infantum zymodeme MON-1 has been shown to be genetically heterogeneous. On the other hand, it was also shown that strains classified as different zymodemes can present the same genotype [18, 20]. However, the low number of typed strains and the lack of consensus on the marker genes, as well as the lack of MLST database for Leishmania, have hindered the implementation of MLST as a reference typing method . Until now, most of the epidemiological data available in the literature is based on MLEE classification, and many DNA-based typing methods correlate data with existing MLEE classification [23, 24].
In this study, we explored the sequences of seven L. infantum genes encoding enzymes considered in MLEE analysis with the aim of (i) enriching the actual databases; (ii) identifying polymorphisms that characterize L. infantum zymodemes most prevalent in the Mediterranean basin; and (iii) developing a rapid screening test that would give results that could be associated with existing MLEE data, overcoming the drawbacks of this laborious technique.
Leishmania infantum DNA
Leishmania infantum strains or isolates used in this study are listed in Table 1. The Chelex-purified DNA from promastigotes of L. infantum strains or isolates was obtained from the National Reference Center for Leishmaniasis (C.Re.Na.L.), Istituto Zooprofilattico Sperimentale della Sicilia (Palermo, Italy). The DNA was quantified using a Qubit fluorometer (Life Technologies, Carlsbad, USA).
Identification of molecular targets and PCR amplification
A preliminary in silico analysis was performed in the GenBank sequence database to search for genes encoding enzymes used in MLEE for L. infantum. These genes were selected not only to gain information on genetic diversity, but also to connect the genotyping results to the existing MLEE data. Six genes were initially selected, based on previous studies performed on L. infantum that allowed the retrieval of available polymorphic sequences from GenBank, i.e. malic enzyme (me), 6-phosphogluconate dehydrogenase (pgd), mitochondrial isocitrate dehydrogenase (icd), glucose-6-phosphate isomerase (gpi), glucose-6-phosphate dehydrogenase (g6pdh) and mannose phosphate isomerase (mpi). Moreover, the phosphoglucomutase (pgm) gene was also taken into consideration in an attempt to enrich the existing database. Primers were designed with Primer-BLAST  using L. infantum MHOM/FR/78/LEM75 or MCAN/AR/10/MDP1 sequences as reference. The primers are listed in Table 2 and their position on the seven target genes are shown in Additional file 1: Figure S1.
Conventional PCR was performed on a total volume of 50 μl with 1–2 μl of template, 200 μM dNTP, 2.5 mM MgCl2, 200 nM of each primer and 1 unit of Hot-Rescue DNA polymerase (Diatheva s.r.l., Fano, Italy). The amplification was performed in a PCR GeneAmp 2700 thermocycler (Applied Biosystems, Foster City, USA). The thermal cycling profile was as follows: 94 °C for 7 min; followed by 35 cycles at 94 °C for 30 s, 60 °C for 20 s or 60 s (depending on amplicon length) and 72 °C for 20 s; with a final extension at 72 °C for 5 min. Each sample was amplified in duplicate. Amplified fragments were analyzed by agarose gel electrophoresis and visualized with GelRed DNA stain (Biotium, Fremont, USA). The gels were visualized under UV light using a gel doc apparatus (Bio-Rad, Hercules, USA). A 100-bp double-stranded DNA ladder or ΦX174 DNA/BSU/ HaeIII/marker 9 (MBI Fermentas, Waltham, USA) was included on the gels as a size standard.
PCR product sequencing and phylogenetic analysis
The amplification products of me, pgd, icd, gpi, g6pdh, mpi and pgm genes for all the strains/isolates are indicated in Table 1 (with the exception of isolate 1 and isolate 2) were purified using a Minelute PCR purification kit (Qiagen, Hilden, Germany) and directly sequenced. DNA sequencing was performed using the BigDye Terminator v. 1.1 Cycle Sequencing Kit on an ABI PRISM 310 Genetic Analyzer (Applied Biosystems). Bases with Phred values < 20 were checked by visual analysis of electropherograms and aligned with BioEdit Sequence Alignment Editor  using default options. Heterozygosity was considered to be present when direct sequencing of PCR product yielded a similar peak at the same site. The sequences were deposited in GenBank with the following accession numbers: MF375413-MF375423 (me), MF280205-MF280215 (pgd), MF347625-MF347635 (icd), MF288905-MF288915 (gpi), MF479731-MF479741 (g6pdh), MF462101-MF462111 (mpi) and MF347614-MF347624 (pgm).
For each housekeeping gene, allele numbers were assigned to unique sequences and a genotype was determined as the combination of the six genes selected for analysis (Additional file 2: Table S1). In total, 54 isolates/strains belonging to 22 different MON, were analyzed; 23 of them had sequences available for me, pgd, icd, gpi, g6pdh and mpi genes. For each of these isolates/strains, we concatenated those sequences to obtain a 3578 bp-long sequence. A maximum likelihood tree was constructed from the concatenated nucleotide dataset using PhyML 3.0 . The best-fit substitution model was determined by the Akaike information criterion (AIC) using Smart Model Selection . The optimal model of evolution was GTR+I with proportion of invariable sites of 0.989 and one category for substitution rate. Bootstrap values were calculated from 100 replications. The tree was visualized using iTOL .
Real-time PCR and high resolution melting (HRM) analysis
A 111 bp internal region of the me gene encompassing the single nucleotide polymorphism (SNP) 390T>G was amplified with primers listed in Table 2. Real-time PCR (named qPCR-MEint) was carried out in a 25 μl volume with 1 μl template DNA and SYBR green PCR master mix (Diatheva s.r.l.) containing 1U Taq Polymerase and 200 nM of each primer. The PCR reactions were performed in a Rotor-Gene 6000 instrument (Corbett Life Science, Mortlake, Australia). The amplification profile was: 94 °C for 10 min; followed by 33 cycles at 94 °C for 30 s, 60 °C for 20 s and 72 °C for 20 s. The reactions were performed in duplicate. After amplification, the high-resolution melting (HRM) analysis was performed over the range 85–95 °C, rising by 0.1 °C/s and waiting for 2 s at each temperature. Raw HRM curves were normalized by the Rotorgene 6000 v.1.7 software. Difference graphs of the normalized curves were obtained using MHOM/FR/78/LEM75 and MHOM/DZ/82/LIPA59 strains as reference curves for genotypes 390T and 390G, respectively.
The qPCR-MEint and HRM analysis was also tested using DNA extracted from 4 clinical samples (whole blood, buffy coat, conjunctival swabs, bone marrow aspirate) as described previously [30, 31]. To ensure adequate sensitivity with these samples, a pre-amplification step (10 cycles) was performed and 2 μl of pre-amplified mixture were used as template in the qPCR described above.
The L. infantum isolates/strains gave single products of expected size with the primer pairs listed in Table 2. The PCR products were sequenced and the sequences were aligned with those obtained from GenBank. Partial sequences of me, pgd, icd, gpi, g6pdh, mpi and pgm genes were obtained from 11 strains/isolates of L. infantum. Seven represented MON-1 zymodeme and the remaining four represented MON-24, MON-29, MON-72 and MON-201 (Table 1). Notably, the following L. infantum sequences from zymodemes not previously reported in the databases were obtained: me sequences from MON-24, MON-72 and MON-201; a pgd sequence from MON-201; icd sequences from MON-24, MON-72 and MON-201; gpi sequences from MON-24, MON-72 and MON-201; g6pdh sequences from MON-24, MON-72 and MON-201; mpi sequences from MON-24, MON-72 and MON-201; and pgm sequences from MON-1, MON-24, MON-29, MON-72 and MON-201.
Concerning the me gene, 13 SNPs were evidenced (Table 3). Nine were silent and four changed the amino acid sequence (V59M; I133S; V330I; E563E/D). These data confirmed previously reported findings obtained with L. donovani complex . All sequences from MON-1 strains were identical, except for the MHOM/TN/80/IPT1 strain. Notably, the genotype heterogeneity of this strain was previously confirmed on the gpi sequence . It is also noteworthy that the two MON-29 strains (MHOM/ES/82/BCN1 and MHOM/FR/1996/LEM3249) showed differences at nucleotide positions 327, 329 and 507, confirming the higher discriminatory power of a sequencing approach compared to MLEE. The pgd gene sequence showed only the silent SNP 678A>G in 11 strains (Table 4). Concerning the icd gene, it is of relevance that two silent SNPs (204T>C, 369T>C) were evidenced in strain MHOM/DZ/82/LIPA59 (MON-24) and three silent SNPs (150G>A, 369T>C, 1038C>A) in strain MHOM/SD/1997/LEM3472 (MON-267) (Table 5). The polymorphisms 204T>C and 1038C>A appeared unique for MON-24 and MON-267, respectively. The gpi gene sequence showed 1503G>T/K and 1831G>A/R SNPs in 4 and 5 strains, respectively (Table 6). The g6pdh gene showed 8 SNPs (Table 7). The mpi gene showed 5 SNPs (Table 8), with SNP 486G>A unique for MON-136. Finally, concerning the pgm gene, no differences were found among the zymodemes examined, including the two sequences available in the database.
The concatenated partial sequences for me, pgd, icd, gpi, g6pdh, mpi from 23 strains were used for phylogenetic tree construction (Fig. 1, Additional file 2: Table S1). The tree showed a 100% bootstrap support for the L. infantum strains and a 75% support for the clade including all strains of MON-1, 72, 201 zymodemes. Within this clade, the MHOM/TN/80/IPT1 strain differed from the bulk of closely related strains (bootstrap support value of 66%) (Fig. 1). This made it possible to identify the most common genotypes in Europe using a simple screening method (see below).
Exploiting sequence polymorphisms for typing
As shown in Fig. 1, the clade including zymodemes MON-1, 72, 201 could be discriminated from all other zymodemes with the exception of the MHOM/TN/80/IPT1 strain. The sequence comparison analysis showed that this group could be discriminated by exploiting the polymorphism 390T>G in the me gene sequence (Table 3). In fact, MON-1, 72, 201 present a T in position 390 (genotype 390T), while the other zymodemes present a G (genotype 390G). Due to this observation, we developed a HRM-based assay to monitor the SNP at position 390. To this end, two internal primers were designed upstream and downstream nucleotide 390, and used for a new qPCR assay (qPCR-MEint) (Additional file 1: Figure S1a). The qPCR-MEint showed a linear correlation between the log of DNA concentration (from 2 to 2 × 10-5 ng/reaction) and Ct values, with a reaction efficiency of 90% (Fig. 2).
HRM analysis allowed us to distinguish amplicons with genotype 390G (i.e. MHOM/TN/80/IPT1 and MHOM/DZ/82/LIPA59 strains) from amplicons with genotype 390T (i.e. MHOM/FR/78/LEM75 and MHOM/IT/86/ISS218 strains) (Fig. 3a). Furthermore, using the PCR product from strains MHOM/DZ/82/LIPA59 and MHOM/FR/78/LEM75 as reference, two L. infantum clinical isolates (isolate 1 and 2)  were typed with qPCR-MEint followed by HRM analysis. The HRM analysis tool of Rotor-Gene 6000 software allowed for the assignment of amplicons of clinical isolate 1 and 2 to genotype 390G and genotype 390T, respectively, with a confidence > 95% in at least one replicate (Fig. 3b). These results have been confirmed by PCR product sequencing (Additional file 3: Figure S2).
Application of HRM analysis in clinical samples
To test the feasibility of the HRM-based typing approach without parasite cultivation, the qPCR-MEint followed by HRM analysis was performed using DNA extracted from canine clinical samples [31, 33] and a human clinical sample (blood), using PCR product from strains MHOM/DZ/82/LIPA59 and MHOM/FR/78/LEM75 as reference. The human peripheral venous blood sample was obtained during routine diagnosis of a patient with VL by venipuncture of the upper limb. The HRM analysis allowed for the assignment of all samples to group 390T (Table 9). The results were confirmed by PCR product sequencing (Additional file 4: Figure S3).
Epidemiological studies examining the dynamic balance between host, vector, and pathogen populations is an important aspect of Leishmania infections, and the importance of reservoirs alternative to the canine population (e.g. lagomorphs, wolves) is currently being investigated [34, 35]. For example, in Italy, L. infantum zymodemes retrieved in infected dogs are primarily represented by MON-1 and MON-72, while infection in humans is caused by a more heterogeneous zymodeme population, suggesting that the canine population is not the only reservoir for all L. infantum zymodemes .
In this study, we investigated the genetic diversity within seven L. infantum genes, aiming to provide information that could aid in the development of tools for fast genetic characterization and epidemiological screening of L. infantum in mammalian host and/or vectors. In particular, we investigated whether it was possible to distinguish the genotype associated with the most common zymodemes in the Mediterranean area (i.e. MON-1 and MON-72), from the genotype associated with all other zymodemes, therefore limiting the number of samples needing costly and time-consuming characterizations. To this end, we collected and aligned all L. infantum gene sequences of seven enzymes (me, pgd, icd, gpi, g6pdh, mpi, pgm) available in GenBank to identify the regions that were more polymorphic. These regions were amplified in 11 other L. infantum strains available in our laboratory to increase the number of sequences to compare. The phylogenetic analysis of the concatenated sequences showed that the group including zymodemes MON-1, 72, 201 clustered separately. Zymodeme MON-201 has been found only in some HIV-positive individuals [12, 36] and in one case of feline leishmaniasis ; therefore, its retrieval may be considered rare. Notably, the genetic distance of MON-81 and MON-267 from other L. infantum zymodemes was confirmed (Fig. 1, Additional file 2: Table S1). In fact, MON-81 and MON-267, initially considered L. infantum, have been successively assigned to L. donovani .
Remarkably, the silent polymorphism at position 390 (T/G) of the me gene distinguishes strains of MON-1, 72, 201 (genotype 390T) from strains of other zymodemes (genotype 390G) (Table 3). This could be important for epidemiological applications, as MON-1 is the most common zymodeme of L. infantum in Europe. The only exception was the reference strain MHOM/TN/80/IPT1, which presents a guanine at nucleotide 390. This may be due to the geographical origin of this strain (Tunisia) compared to the European strains/isolates. It is also noteworthy that, among the MON-1 strains, only strain MHOM/TN/80/IPT1 cannot be distinguished from the MON-98 strains (considering that only me sequences were available for this MON) (Table 3). This could be in agreement with results reported by Haralambous et al. , which showed that some MON-1 and MON-98 zymodemes were not distinguished with a PCR-based typing method using the K26 antigen-coding gene as target sequence. Notably, the use of a real-time PCR-based assay with internal primers followed by HRM analysis allowed for the quick differentiation of sequences with the polymorphysm 390T>G not only in L. infantum isolates/strains, but also in clinical samples, enabling the potential application of this approach in population screening without parasite isolation and culturing.
Many marker sequences can be taken into consideration for genotyping purposes. The most suitable marker sequence depends on the level of resolution needed: highly variable markers can distinguish closely related genotypes but may not faithfully represent more distant relationships, while moderately variable markers will not be able to evidence the differences between closely related samples. The assay that we propose would not be able to distinguish closely related genotypes; therefore, it could be used to differentiate between relapse and re-infection only in cases where different 390T/G genotypes were found. However, it could be useful for epidemiological screening in regions where MLEE data are available, to rapidly identify zymodemes different from MON-1, 72 or 201 (e.g. for rapid screening of L. infantum genotypes in mammalian hosts). Once attributed a genotype (390T or 390G), other tests (e.g. MLEE or MLST) can be performed, depending on the level of resolution needed and on the classification system with which the results are to be compared.
Other zymodemes (i.e. MON-24, MON-78 and MON-199) could be identified by exploiting the unique SNP found across the seven gene sequences. In particular MON-24, which is another common zymodeme [12, 38], could be differentiated by exploiting the polymorphisms 204T>C in the icd gene (Table 5) by using a HRM-based assay similar to the assay described here.
A possible limitation of this study is the unavailability of some zymodeme sequences in the database. Moreover, in many cases, only one strain with sequence information per zymodeme was available. This is important to keep in mind, since strains of the same zymodeme could present different genotypes (e.g. MON-29 in me, icd and g6pdh genes) (Additional file 2: Table S1).
A total of 77 new sequences of seven genes encoding for metabolic enzymes in L. infantum isolates/strains have been produced and deposited in the GenBank database. The analysis of these sequences, together with sequences available in the database, allowed for the identification of genetic polymorphisms exploitable to differentiate the most common L. infantum zymodemes in the Mediterranean basin. In particular, a HRM-based assay aimed to differentiate the genotype 390T and 390G in the me gene was developed. The genotype 390T correlated with zymodemes MON-1, 72, 201, allowing for the rapid identification of the majority of L. infantum genotypes. Once a parasite is attributed to genotype 390T or 390G, other tests (e.g. MLEE, MLST or MLMT) can be performed if more detailed information is needed. This assay has been successfully applied to clinical samples, demonstrating its potential applicability in investigating the role of other mammalian hosts in epidemiological screening.
Quantitative real-time PCR
Multilocus enzyme electrophoresis
Multilocus sequence typing
Multilocus microsatellite typing
Akhoundi M, Kuhls K, Cannet A, Votýpka J, Marty P, Delaunay P, et al. A historical overview of the classification, evolution, and dispersion of Leishmania parasites and sandflies. PLoS Negl Trop Dis. 2016;10:e0004349.
Kostygov AY, Yurchenko V. Revised classification of the subfamily Leishmaniinae (Trypanosomatidae). Folia Parasitol (Praha). 2017;64:020.
Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, et al. Leishmaniasis worldwide and global estimates of its incidence. PLoS One. 2012;7:e35671.
World Health Organization. Control of the leishmaniases. World Health Organ Tech Rep Ser. 2010;949:xii–xiii 1-186.
Rioux JA, Lanotte G, Serres E, Pratlong F, Bastien P, Perieres J. Taxonomy of Leishmania. Use of isoenzymes. Suggestions for a new classification. Ann Parasitol Hum Comp. 1990;65:111–25.
Aït-Oudhia K, Harrat Z, Benikhlef R, Dedet JP, Pratlong F. Canine Leishmania infantum enzymatic polymorphism: a review including 1023 strains of the Mediterranean area, with special reference to Algeria. Acta Trop. 2011;118:80–6.
Millán J, Ferroglio E, Solano-Gallego L. Role of wildlife in the epidemiology of Leishmania infantum infection in Europe. Parasitol Res. 2014;113:2005–14.
Pratlong F, Lami P, Ravel C, Balard Y, Dereure J, Serres G, et al. Geographical distribution and epidemiological features of Old World Leishmania infantum and Leishmania donovani foci, based on the isoenzyme analysis of 2277 strains. Parasitology. 2013;140:423–34.
Kuhls K, Alam MZ, Cupolillo E, Ferreira GEM, Mauricio IL, Oddone R, et al. Comparative microsatellite typing of New World Leishmania infantum reveals low heterogeneity among populations and its recent Old World origin. PLoS Negl Trop Dis. 2011;5:e1155.
Pratlong F, Rioux J-A, Marty P, Faraut-Gambarelli F, Dereure J, Lanotte G, et al. Isoenzymatic analysis of 712 strains of Leishmania infantum in the south of France and relationship of enzymatic polymorphism to clinical and epidemiological features. J Clin Microbiol. 2004;42:4077–82.
Campino L, Pratlong F, Abranches P, Rioux J-A, Santos-Gomes G, Alves-Pires C, et al. Leishmaniasis in Portugal: enzyme polymorphism of Leishmania infantum based on the identification of 213 strains. Trop Med Int Heal. 2006;11:1708–14.
Chicharro C, Jiménez MI, Alvar J. Iso-enzymatic variability of Leishmania infantum in Spain. Ann Trop Med Parasitol. 2003;97(Suppl. 1):57–64.
Martin-Sanchez J, Gramiccia M, Di Muccio T, Ludovisi A, Morillas-Márquez F. Isoenzymatic polymorphism of Leishmania infantum in southern Spain. Trans R Soc Trop Med Hyg. 2004;98:228–32.
Alvar J, Cañavate C, Gutiérrez-Solar B, Jiménez M, Laguna F, López-Vélez R, et al. Leishmania and human immunodeficiency virus coinfection: the first 10 years. Clin Microbiol Rev. 1997;10:298–319.
Gramiccia M, Scalone A, Di Muccio T, Orsini S, Fiorentino E, Gradoni L. The burden of visceral leishmaniasis in Italy from 1982 to 2012: a retrospective analysis of the multi-annual epidemic that occurred from 1989 to 2009. Euro Surveill. 2013;18:20535.
Botilde Y, Laurent T, Quispe Tintaya W, Chicharro C, Cañavate C, Cruz I, et al. Comparison of molecular markers for strain typing of Leishmania infantum. Infect Genet Evol. 2006;6:440–6.
Cortes S, Maurício IL, Kuhls K, Nunes M, Lopes C, Marcos M, et al. Genetic diversity evaluation on Portuguese Leishmania infantum strains by multilocus microsatellite typing. Infect Genet Evol. 2014;26:20–31.
Mauricio IL, Yeo M, Baghaei M, Doto D, Pratlong F, Zemanova E, et al. Towards multilocus sequence typing of the Leishmania donovani complex: resolving genotypes and haplotypes for five polymorphic metabolic enzymes (ASAT, GPI, NH1, NH2, PGD). Int J Parasitol. 2006;36:757–69.
Zemanová E, Jirků M, Mauricio IL, Horák A, Miles MA, Lukeš J. The Leishmania donovani complex: genotypes of five metabolic enzymes (ICD, ME, MPI, G6PDH, and FH), new targets for multilocus sequence typing. Int J Parasitol. 2007;37:149–60.
El Baidouri F, Diancourt L, Berry V, Chevenet F, Pratlong F, Marty P, et al. Genetic structure and evolution of the Leishmania genus in Africa and Eurasia: what does MLSA tell us. PLoS Negl Trop Dis. 2013;7:e2255.
Schonian G, Kuhls K, Mauricio IL. Molecular approaches for a better understanding of the epidemiology and population genetics of Leishmania. Parasitology. 2011;138:405–25.
Herrera G, Hernández C, Ayala MS, Flórez C, Teherán AA, Ramírez JD. Evaluation of a Multilocus Sequence Typing (MLST) scheme for Leishmania (Viannia) braziliensis and Leishmania (Viannia) panamensis in Colombia. Parasit Vectors. 2017;10:236.
Chicharro C, Llanes-Acevedo P, García E, Nieto J, Moreno J, Cruz I. Molecular typing of Leishmania infantum isolates from a leishmaniasis outbreak in Madrid, Spain, 2009 to 2012. Euro Surveill. 2013;18:20545.
Haralambous C, Antoniou M, Pratlong F, Dedet J-PP, Soteriadou K. Development of a molecular assay specific for the Leishmania donovani complex that discriminates L. donovani/Leishmania infantum zymodemes: a useful tool for typing MON-1. Diagn Microbiol Infect Dis. 2008;60:33–42.
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13:134.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.
Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.
Lefort V, Longueville J-E, Gascuel O. SMS: Smart Model Selection in PhyML. Mol Biol Evol. 2017;34:2422–4.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44:W242–5.
Ceccarelli M, Galluzzi L, Migliazzo A, Magnani M. Detection and characterization of Leishmania (Leishmania) and Leishmania (Viannia) by SYBR green-based real-time PCR and high resolution melt analysis targeting kinetoplast minicircle DNA. PLoS One. 2014;9:e88845.
Ceccarelli M, Galluzzi L, Diotallevi A, Gasparini E, Migliazzo A, Magnani M. The relevance of molecular diagnosis in a dog vaccinated against leishmaniasis. Vet Med Anim Sci. 2016;4:4.
Galluzzi L, Diotallevi A, De Santi M, Ceccarelli M, Vitale F, Brandi G, et al. Leishmania infantum induces mild unfolded protein response in infected macrophages. PLoS One. 2016;11:e0168339.
Ceccarelli M, Galluzzi L, Sisti D, Bianchi B, Magnani M. Application of qPCR in conjunctival swab samples for the evaluation of canine leishmaniasis in borderline cases or disease relapse and correlation with clinical parameters. Parasit Vectors. 2014;7:460.
Miró G, Müller A, Montoya A, Checa R, Marino V, Marino E, et al. Epidemiological role of dogs since the human leishmaniosis outbreak in Madrid. Parasit Vectors. 2017;10:209.
Oleaga A, Zanet S, Espí A, Pegoraro de Macedo MR, Gortázar C, Ferroglio E. Leishmania in wolves in northern Spain: a spreading zoonosis evidenced by wildlife sanitary surveillance. Vet Parasitol. 2018;255:26–31.
Gramiccia M. The identification and variability of the parasites causing leishmaniasis in HIV-positive patients in Italy. Ann Trop Med Parasitol. 2003;97(Suppl. 1):65–73.
Pennisi M-G, Cardoso L, Baneth G, Bourdeau P, Koutinas A, Miró G, et al. LeishVet update and recommendations on feline leishmaniosis. Parasit Vectors. 2015;8:302.
Haralambous C, Dakkak A, Pratlong F, Dedet J-P, Soteriadou K. First detection and genetic typing of Leishmania infantum MON-24 in a dog from the Moroccan Mediterranean coast: genetic diversity of MON-24. Acta Trop. 2007;103:69–79.
We thank Dr Federica Bruno and Dr Germano Castelli from IZS Palermo for providing Leishmania spp. DNA. We thank Dr Emanule Gasparini and Dr Barbara Bianchi from the Santa Teresa Veterinary Clinic, Fano (PU) Italy, for providing canine clinical samples. We also thank Dr Benedetta Canovari and Professor Francesco Barchiesi from Ospedali Riuniti Marche Nord for providing a blood sample from a human VL patient, and Dr Elizabeth Tremmel for English revision.
This work was partially supported by the Italian Ministry of Health (grant number Ricerca corrente 2013 IZS SI 3/14 RC). The funder had no role in study design, data collection and analysis, interpretation of the data, or writing the manuscript.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its Additional files.
Approval of the study was obtained on 31st July 2012 from the Ethical Committee for Animal Experiments of the University of Urbino (CESA) [study title “Diagnosi biomolecolare della leishmaniosi attraverso l’uso di campioni clinici non invasivi e loro utilizzo per il monitoraggio terapeutico” (Prot. CESA 2/2012)].
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Reference sequences and position of primers used in this study. a L. infantum strain MHOM/FR/1978/LEM75 cytosolic NADP-malic enzyme (me) gene, partial cds. (GenBank: DQ449701.1). Primers amplifying 5' region and 3' region are boxed and underlined, respectively. Internal primers are in bold. b L. infantum gene for phosphogluconate dehydrogenase (decarboxylating), strain MHOM/FR/1978/LEM75 (GenBank: AM157139.1). c L. infantum strain MHOM/FR/1978/LEM75 mitochondrial isocitrate dehydrogenase (icd) gene, complete cds (GenBank: DQ449672.1). d L. infantum gpi gene for glucose-6-phosphate isomerase, strain MHOM/FR/1978/LEM75 (GenBank: AJ620617.1). e L. infantum strain MHOM/FR/1978/LEM75 glucose-6-phosphate dehydrogenase (g6pdh) gene, complete cds (GenBank:DQ449770.1). f L. infantum strain MHOM/FR/1978/LEM75 mannose phosphate isomerase (mpi) gene, complete cds (GenBank: DQ449737.1). g L. infantum isolate MCAN/AR/10/MDP1 phosphoglucomutase gene, partial cds (GenBank: KJ643214.1). (DOCX 20 kb)
Table S1. Summary of allelic profile and genotype of L. infantum zymodemes considered in this study. (DOCX 18 kb)
Figure S2. Partial sequences of qPCR-MEint amplification products from reference strains and clinical isolates. Electropherograms encompassing polymorphic nucleotide at position 390 (arrows) are represented. (PDF 72 kb)
Figure S3. Partial sequences of qPCR-MEint amplification products from clinical samples. Electropherograms encompassing polymorphic nucleotide at position 390 (arrows) are represented. (PDF 570 kb)