- Open Access
Evaluation of a Multilocus Sequence Typing (MLST) scheme for Leishmania (Viannia) braziliensis and Leishmania (Viannia) panamensis in Colombia
Parasites & Vectorsvolume 10, Article number: 236 (2017)
Leishmaniases are parasitic vector-borne diseases affecting more than 12 million people in 98 countries. In Colombia, leishmaniasis is widespread and the most common clinical manifestation is cutaneous, mainly caused by L. panamensis and L. braziliensis. Currently, the genetic diversity of these species in Colombia is unknown. To address this, we applied molecular techniques for their characterization, using multilocus sequence typing (MLST) to explore the genetic variability and phylodynamics of the disease.
Seven previously described genetic markers were selected highlighting the implementation of a mitochondrial marker. Markers were applied to 163 samples from isolates obtained between 1980 and 2001.
The identification of the samples showed an excellent correlation with typing tests previously applied (MLEE, monoclonal antibodies). Isolates of L. braziliensis showed greater genetic diversity than L. panamensis, and a greater number of diploid sequence types (DSTs). In addition, the geographical distribution of DSTs for each species were obtained through georeferencing maps.
To our knowldge, this study represents the first description of the genetic variability of L. panamensis in Colombia and South America, and is the first to propose a scheme of MLST for epidemiological surveillance of leishmaniasis in the country.
Leishmaniases are a group of diseases caused by Leishmania parasites and transmitted to humans and other mammals by the bite of psychodid insects . Currently, leishmaniases are ranked among the most neglected diseases constituting a serious public health problem worldwide because of high morbidity and the fact that the areas where most cases occur are usually rural or suburban, hindering the diagnosis and treatment of infected patients. The clinical presentations of the disease are cutaneous leishmaniasis (CL), mucocutaneous leishmaniasis (MCL) and visceral leishmaniasis (VL) [1, 2]. These pathologies affect more than 12 million people in 98 countries and is estimated that 350 million individuals are at risk of acquiring the infection, with an incidence of 1.3 million of new cases per year [3, 4]. In the Americas, CL and MCL are endemic in 18 countries with a distribution of cases from Mexico to Argentina, where Brazil, Colombia and Peru report the highest percentage of cases . In Colombia, it is estimated that the population at risk of infection is 12,277,606 people for CL and 6,795,047 people for MCL respectively .
The genus Leishmania includes more than 20 species grouped in complexes constantly under review and strong debate . In Colombia, L. panamensis, L. braziliensis, L. lainsoni, L. guyanensis, L. infantum chagasi, L. equatoriensis, L. mexicana, L. amazonensis and L. colombiensis have been reported, with L. panamensis and L. braziliensis as the most prevalent and mainly associated to CL [7, 8]. However, many studies in the country do not report the species involved in the transmission, which has hindered the association with the clinical presentation, profile of susceptibility to treatment and epidemiological management of outbreaks [9,10,11]. In this sense, molecular biology tools have been shown to be very useful, not only for diagnosis, but also for the identification of species . In addition to multilocus enzyme electrophoresis (MLEE), considered the gold standard test , other techniques such as multilocus sequence typing (MLST) have been used as a valuable tool for the characterization of the L. donovani complex [14, 15]. However, the lack of databases for Leishmania, added to the low number of typed strains and the lack of consensus on the markers has hindered the implementation of this technique as a reference method. Parallel to MLST, multilocus sequence analysis (MLSA) has been employed in the study of Leishmania braziliensis for analysis of genetic diversity and outbreaks surveillance [16, 17]. This methodology, widely used for several pathogens [18,19,20], has become an important tool for the study of Leishmania [16, 17]. Nevertheless, few studies have been reported in South America (Brazil and Argentina). In Colombia, no MLST or MLSA have been applied resulting in the absence of information concerning the transmission dynamics of leishmaniasis in the country.
Therefore, the aims of this study were to standardize and evaluate an MLST scheme based on seven markers in L. (Viannia) panamensis and L. (Viannia) braziliensis strains from different geographical regions of Colombia. To determine the circulating diploid sequence types (DSTs) and strengthen the epidemiological surveillance in the country; and finally, to provide further data for the understanding of the genetic diversity of Leishmania parasites in South America.
In this study, we included a total of 163 isolates of L. panamensis (105 isolates) and L. braziliensis (58 isolates) from human cutaneous leishmaniasis cases obtained by the National Institute of Health of Colombia (NIHC) as part of epidemiological surveillance between 1980 and 2001. Isolates were retrieved from patients from 20 departments (Antioquia, Bolivar, Boyaca, Caldas, Caquetá, Casanare, Chocó, Córdoba, Cundinamarca, Guainía, Guaviare, Magdalena, Meta, Norte de Santander, Putumayo, Risaralda, Santander, Sucre, Tolima and Vichada).
Parasite isolation, DNA extraction and typing
Punch biopsies were triturated in sterile Ten Broeck homogenizers containing phosphate buffered saline (PBS), gentamicin (40 μg/ml), and 5-fluorocytosine (500 μg/ml). The resultant tissue suspension was inoculated directly into two tubes of NNN medium. Leishmania panamensis and L. braziliensis species identification was performed for each isolate by direct Sanger sequencing of cytochrome b, heat shock protein 70 (HSP70) and MLEE as described by Ramirez et al. . DNA extraction was performed using the Quick-DNA™ Universal Kit (Zymo Research, Orange, CA, USA) following the manufacturer's protocol.
Marker selection, amplification by PCR and sequencing
Amplification was performed based on specific primers for gene fragments coding for metabolic enzymes such as glucose 6-phosphate dehydrogenase (g6pdh, EC 18.104.22.168), phosphoglucomutase (pgm, EC 22.214.171.124), mannose phosphate isomerase (mpi, EC 126.96.36.199), alanino-transferase (alat, EC 188.8.131.52), aspartate aminotransferase (asat, EC 184.108.40.206), isocitrate dehydrogenase (icd, EC 220.127.116.11) and cytochrome b (cytb) molecule as previously reported [15, 16, 21] (Table 1). The PCR reactions were performed in a final volume of 25 μl containing 12.5 μl GoTaq® Green Master Mix (Promega, Madison, WI, USA), 0.5 μM of each primer and 1.25 μl of DNA (<250 ng). Amplification was conducted in a Thermal Cycler S1000™ (Bio-Rad Laboratories, Inc., Hercules, CA, USA) with an initial denaturation at 95 °C for 5 min, followed by 40 cycles of denaturation at 95 °C for 1 min, annealing at 60 °C for 1 min and extension at 72 °C for 1 min. Finally, an extension at 72 °C was performed for 10 min. Amplification and amplicon size was verified by electrophoresis on agarose gel stained with SYBR Safe DNA Gel Stain (Life Technologies, Carlsband, CA, USA) and a molecular weight marker (Promega) adding 2 μl of each PCR product. The remaining 23 μl were purified using the PCR kit ExoSAP-IT® Product Cleanup (Affymetrix, Santa Clara, CA, USA). Sequencing of both strands was carried out in the Unit of Sequencing and Genomic Analysis of the National Institute of Health, Colombia (NIHC) using the Big Dye Terminator v3.1 kit and ABI 3130 Genetic Analyzer (Applied Biosystems, California, USA).
Sequence assembly and genetic diversity indices
The assembly of the consensus sequences using the forward and reverse reads, alignment and sequence editing was performed using MEGA 5.0 . Additionally, electropherograms were revised in order to detect ambiguous sites (heterozygosity) evidenced by two overlapped peaks in one position. Reconstruction of haplotypes was performed using the algorithm PHASE, through DnaSP v.5.0 . The indices of genetic diversity of the sequences (total number of mutations Eta, haplotype diversity Hd, genetic diversity ɵ and π and the rate of synonymous and non-synonymous substitutions dN/dS) were calculated for each marker, each species and by geographical region (Andean, Orinoquia, Caribbean, Pacific and Amazon), using DnaSP v.5.0 . Nucleotide diversity (π) is defined as the average number of nucleotide differences per site between two sequences. The θ index is defined as an indicator of mutation rate per nucleotide site per generation. Eta (h) is defined as the total number of mutations and S the number of segregating (polymorphic) sites. Comparison charts of haplotype diversity for each gene were constructed for each species using Microsoft Excel 2013. Additionally, using the same tool, sample distributions were constructed by geographical region.
Allelic profile coding, geographical and temporal distribution of DSTs
MLSTest (https://mlstest.codeplex.com/) bioinformatic tool was used for DST assignment, based on the allelic profile of each sample using the seven gene fragments of the MLST scheme herein proposed. The DST is defined as the combination of the alleles from the different gene fragments used for an MLST scheme. Similarly, the number of alleles and polymorphic nucleotide numbers were calculated . In addition, a map of DSTs distribution of alleles for each species, based on the georeferenced coordinates of each of the samples was constructed for L. panamensis and L. braziliensis using the ArcGIS10.3 program (http://www.esri.com/ArcGIS10.3). Finally, we assessed the temporal variation of the different DSTs for each species across the 21 years of sampling.
Phylogenetic analysis and clonal complexes determination
The final set of sequences from each individual gene fragment and the concatenated alignment were evaluated in ModelTest 3.7 where the most appropriate evolutionary model was selected based on the AIC (Akaike information criterion). A maximum composite likelihood (MCL) analysis using a Tamura-3 parameter model and the Neighbour-Joining algorithm was run in RAxML 7.2.5 on the CIPRES project (Cyberinfrastructure for Phylogenetic Research) portal 2.0 servers. To evaluate the robustness of the nodes in the resulting phylogenetic trees, 1,000 bootstrap replicates were performed. The final trees were rooted with sequences for Leishmania donovani, Leishmania guyanensis, Leishmania tropica, Leishmania major, Leishmania amazonensis and Leishmania mexicana, retrieved from TriTrypDB. DSTs obtained based on the sequences of the seven genes for each species were subjected to eBURST program (http://eburst.mlst.net) to presumably determine the evolutionary relationships among the samples analyzed. eBURST predicts descent from a genotype "founder" for all other genotypes, establishing the groups of related strains which share certain number of identical alleles with other members (Clonal Complex).
Assesment of clonality via allelic plot
The population genetics of Leishmania is complex. Therefore, we decided to calculate the rate of topology incongruence as a measure of clonality vs sexuality across the seven markers employed. We constructed MCL phylogenetic trees from each gene fragment and the allelic diversity was measured in terms of ‘allelic types’ supported as nodes with bootstrap values over 80% and colored appropriately depending on the allelic type within L. panamensis (red) and L. braziliensis (blue); this allelic diversity patterns was used to construct an allelic plot to determine the extent of allelic repertoire exchange at inter and intra species level. Finally, based on the combination of different shades of blue and/or red (different number of alleles supported by bootstrap values equal to or over 80%), the number of multilocus genotypes (MLG) were calculated manually. Here MLGs are defined as the combination of statistically supported alleles based on phylogenetic reconstruction using maximum likelihood estimations.
Genetic diversity of the samples analyzed
The seven markers evaluated in the present study were successfully amplified in all the samples tested. Genetic diversity calculations showed a greater degree of intraspecific variation in L. braziliensis than in L. panamensis. However, out of the seven markers evaluated only cytochrome b showed greater diversity in L. panamensis than in L. braziliensis (Table 2). Nucleotide diversity indices π and θ in L. panamensis were on average 0.016499 and 0.00815, respectively, while in L. braziliensis were 0.02452 and 0.03337, respectively. Fragments of the tested genes showed between 6 and 80 polymorphic sites in L. panamensis (80 only in cytb) and between 5 and 21 polymorphic sites in L. braziliensis. On average, the ratio of the rate of nonsynonymous and synonymous substitutions was higher in L. panamensis (dN/dS = 0.1705) than in L. braziliensis (dN/dS = 0.1261). The pgm and alat genes in L. panamensis showed the highest values of haplotype diversity, followed by asat and g6pdh genes, whereas genes in L. braziliensis showed more haplotype diversity across mpi, asat and cytb (Additional file 1: Figure S1).
Genetic diversity by geographical regions
When comparing the genetic diversity by geographical region (Fig. 1), it was found that L. panamensis had a high degree of haplotype diversity among samples from the regions of the Amazon and Orinoco (Hd = 1 in both regions, N = 9 and N = 12, respectively) while samples from the Andean region showed less diversity (Hd = 0.118, N = 36). When calculating nucleotide diversity indices π and θ the trend remains only in the Andean region, while in the Orinoco region is reversed. Meanwhile, samples of L. braziliensis showed greater haplotype diversity in the Andean region (Hd = 0.564, N = 36) and Amazon region (Hd = 0.534, N = 9), while those from the Caribbean region showed a lack of haplotype diversity (H = 0, N = 1). As for nucleotide diversity indices π and θ, the samples from the Andean region showed greater diversity (0.019 and 0.033, respectively), followed by the samples from the Amazon region (0.00451 and 0.00342, respectively).
Geographical and temporal distribution of DSTs
The results from MLSTest showed a total of seven DSTs for L. braziliensis and five DSTs for L. panamensis (due to the lack of a DST database, the DSTs detected were included in Additional file 2: Tables S1, S2). Based on the DSTs obtained and georeferencing tool data, distribution maps were constructed, showing the distribution of six of the seven DSTs of L. braziliensis in the Andean region, with major density in three departments (Antioquia, Cundinamarca and Norte de Santander). DST 2 is uniformly distributed in four of the five regions of Colombia, while the DST 3 is restricted to the Orinoco region; and DSTs 6 and 7 are found only in the department of Norte de Santander (Fig. 2a). In contrast, the distribution of DSTs of L. panamensis showed that four of the five DSTs of this species are circulating in the Andean region, with a trend to cluster in the central and western departments of the country. DST 1 was found circulating in the Andean, Caribbean, Orinoco and Pacific regions, while the DST 5 is restricted to the Pacific region (Fig. 2b). We determined the temporal distribution of the DSTs detected (Fig. 3). No association between the DSTs and the year of isolation was observed. Two patterns were depicted: (i) in L. panamensis, DST1 was constant across the years but DST 5 only appears in 1996 (Fig. 3a); (ii) in L. braziliensis DST 6 and DST 7 only appear in 1985 (Fig. 3b).
Phylogenetic analyses and clonal complexes
Phylogenetic trees were constructed for each gene fragment including one final concatenated tree (Fig. 4a). The analysis of L. braziliensis clonal complex evidenced the presence of two main complexes followed by another five smaller groups. In contrast, L. panamensis had only one major clonal complex composed of a predominant genotype with few (4) relatives (Fig. 4b).
Assesment of clonality via allelic plot
Phylogenetic trees were constructed for each gene fragment in order to depict an allelic plot that showed the incongruences in the topologies as a measure of clonality vs sexuality (Fig. 5). Overall, the allelic plot for L. panamensis showed conservation in the topologies, with few swapping events and exchange of alleles at intra- and interspecies level (Fig. 5). In contrast, L. braziliensis showed gross exchange of alleles, mainly presented in the genes encoding g6pdh and asat. Also, the number of alleles was higher in asat than in the other genes tested. In this regard, the genes coding for isocitrate dehydrogenase and phosphoglucomutase in L. braziliensis showed conservation in the topologies among the tested samples, while all other genes are events of alleles swapping at intra- and interspecific level with presence of L. panamensis alleles. Finally, when calculating the MLGs, a total of 16 MLGs were detected for L. braziliensis and 8 MLGs for L. panamensis.
Despite the many existing tools for classifying Leishmania species, its taxonomy remains contentious. Its complex biology and differing use and interpretation of the biological, biochemical and genetic markers used by the scientific community have not allowed a consensus to be reached, leading some to suggest a need for complete genome sequencing [25,26,27]. Precisely, this lack of consensus has prevented modification of the current gold standard for typing (MLEE) despite difficulties with its implementation and complex analysis. However, molecular techniques such as MLST are considered useful tools for studies of outbreaks and improving the understanding of the taxonomy of Leishmania [25, 26, 28]. Despite the benefits of MLST, few studies have been published, most of them focused on the study of Old World Leishmania species and only a couple in South America (Brazil and Argentina) [16, 17, 21]. To our knowledge, this is the first time MLST has been applied to Leishmania spp. in Colombia. Herein, we propose a scheme with seven previously reported different markers for typing the two most prevalent species (L. panamensis and L. braziliensis) (Table 1) [15,16,17, 21, 29]. Among the markers employed, the use of cytb gene is highlighted due to its specificity and sensitivity including its character of uniparental inheritance [30,31,32,33].
Leishmania panamensis showed lower diversity than L. braziliensis, this was clearly depicted by comparing rates of haplotype diversity of different genes (Table 2). The greatest value of diversity in L. panamensis was observed in pgm and alat markers (Hd: 0.234), which was matched only by the diversity index of alat in L. braziliensis (Hd: 0.223) (Table 2; Additional file 1: Figure S1). Nevertheless, when comparing the genetic variability of L. braziliensis with previously reported data from Argentina and Brazil, the latter exhibit a much greater diversity than in the present study. As shown by Marco et al. , haplotype diversity indices for alat and pgm markers were 0.89 and 0.921, respectively, compared with 0.223 and 0.324 of this study. In this sense, the haplotype diversity values for g6pdh, mpi and icd were closer to those reported by Boite et al.  (Table 2; Additional file 1: Figure S1). The differences between the values reported in each study are attributed not only to the fact that genes are different, but also in the study of Marco et al. estimates of genetic diversity were performed based on different species of the subgenus Viannia, resulting in an increase in these values of diversity . The other fact is the vast diversity of L. braziliensis across South America. Van der Auwera et al.  suggested the presence of three different groups within the same species based on the heat-shock protein 70 (HSP70), demonstrating the essential role that geography and ecology can play as a major driver of genetic diversity. However, the genetic variability within L. braziliensis will only be resolved when a consensus MLST scheme is applied to isolates across South America. It is noteworthy that g6pdh marker showed more diversity in L. braziliensis than in L. panamensis. This is possibly due to rearrangements present in chromosomes 20 and 34 in L. braziliensis where this gene is located [34, 35].
Regarding the geographical distribution of DSTs, it is noted that L. braziliensis DST 1 is concentrated in the center of the country (Fig. 2), wheres DTS 2 is widely dispersed. The heterogeneous distribution covers the Caribbean, Andean, Orinoco and Amazon regions. DSTs 4 and 5 were found only in the department of Antioquia. This geographical restriction to a single department was also evident for the DSTs 6 and 7. In contrast, the L. panamensis DST 1 is distributed in the Pacific, Andean, Caribbean and Orinoco, while the DST 3 was found in central Colombia (Cundinamarca department). The DSTs 2 and 4 were distributed in departments of the Andean region, while the DST 5 was restricted to the department of Choco (Fig. 2a, b). When these dispersions are compared with the geographical distribution of Leishmania vectors in Colombia reported by Ferro et al. , we observed that L. braziliensis DSTs 6 and 7 are restricted to the department of Norte de Santander overlapping with the distribution of Pintomyia spinicrassa in that department; and no other DSTs were found in the other departments (Cundinamarca and Boyacá) . Likewise, a slight overlap of L. panamensis DST 3 was observed with the vector P. ovallesi in the department of Cundinamarca , indicating that this sandfly species might be a potential vector of L. panamensis . This suggests a possible selection by the vector of some genetic variants. However, MLST studies in vector isolates are needed to confirm this hypothesis. Regarding the temporal variation, no strict associations were observed based on the DSTs; this can be explained by the strict selection of isolates in this study, a cross-sectional study is required to determine the true temporal distribution of DSTs in Colombia (Fig. 3).
The phylogenetic reconstruction of the concatenated alignment (Fig. 4a) demonstrated the genetic relatedness of the DSTs depicted. The MLST scheme had the power and resolution to readily discriminate L. braziliensis and L. panamensis, including L. guyanensis. Different authors have proposed that L. guyanensis and L. panamensis might be the same species, but the MLST scheme herein proposed was able to discriminate these two species. However, more L. guyanensis isolates are required to fully corroborate this statement. The calculation of the number of clonal complexes confirmed the high diversity of L. braziliensis compared to L. panamensis (Fig. 3b). Our results reinforce the need to establish a consensus MLST scheme for a proper epidemiological surveillance of cutaneous leishmaniasis, not only in Colombia but at a continental level.
The study of the genetic diversity and population structure of Leishmania is complex due to the absence of accurate population genetics tools that correlate with the biology of the protozoan parasites. The type of reproduction of Leishmania is controversial, some authors have demonstrated genetic exchange in Leishmania  but others suggest preponderant clonal evolution . Recently, whole genome sequencing demonstrated the capacity of hybridization in Leishmania favouring a sexual propagation theory in protozoan parasites . Our results are consistent with previous results that highlight the ability of Leishmania for genetic exchange, where different alleles of L. panamensis were observed in the genes of L. braziliensis and vice versa (Fig. 5). Clear events of intra- an inter-species allelic exchange were observed. This phenomenon of alleles swapping/exchange was more frequently observed in L. braziliensis. This has been also demonstrated in Brazil where L. braziliensis showed a high level of genetic exchange . Herein, we showed the heightened potential of Leishmania braziliensis for genetic exchange and the reduced capacity for genetic exchange observed in L. panamensis. Such instances of genetic exchange can impact the distribution of DSTs, since shifting of alleles may generate novel DSTs across time and space. However, that is exactly the benefit of MLST as typing tool in molecular studies since it could be pivotal to understanding emerging strains and outbreaks arising from genetic exchange across endemic regions.
To our knowledge, this is the first study of MLST in Leishmania species from Colombia, which enables a robust and reliable strategy to identify the two most prevalent species in the country. The phylogenetic analyses suggest its applicability to other New World species of Leishmania, but a set of well-characterized isolates of these species are required to fulfill this statement. It is also the first study to determine intraspecific genetic diversity of L. panamensis, and the high genetic diversity of L. braziliensis in Colombia. The usefulness of MLST is demonstrated and we encourage the Colombian government to establish a database that can be used for epidemiological surveillance. This will afford establishment of measures for prevention and control; towards early detection of drug resistant and virulent strains that pose substantial threats to the susceptible population that suffer from cutaneous leishmaniasis.
- cytb :
Diploid sequence type
Glucose 6-phosphate dehydrogenase
Instituto Nacional de Salud
Multilocus enzyme electrophoresis
Multilocus sequence analysis
Multilocus sequence typing
Mannose phosphate isomerase
WHO, World Health Organization. Control of the leishmaniases. Report of a meeting of the WHO Expert Committee on the Control of Leishmaniases. 2010.
Assembly WH. The World Health Assembly Resolution (WHA60.13) on the “Control of leishmaniasis”. Geneva, Switzerland. 2007:5. http://www.who.int/neglected_diseases/mediacentre/WHA_60.13_Eng.pdf.
Desjeux P. Leishmaniasis: current situation and new perspectives. Comp Immunol Microbiol Infect Dis. 2004;275:305–18.
Alvar J, Velez ID, Bern C, Herrero M, Desjeux P, Cano J, et al. Leishmaniasis worldwide and global estimates of its incidence. PLoS One. 2012;75, e35671.
Pan American Health Organization. Leishmaniases: Epidemiological report of the Americas. Report N° 3 - July. 2015.
Schonian G, Kuhls K, Mauricio IL. Molecular approaches for a better understanding of the epidemiology and population genetics of Leishmania. Parasitology. 2011;1384:405–25.
Ramirez JD, Hernandez C, Leon CM, Ayala MS, Florez C, Gonzalez C. Taxonomy, diversity, temporal and geographical distribution of cutaneous leishmaniasis in Colombia: A retrospective study. Sci Rep. 2016;6:28266.
Ovalle CE, Porras L, Rey M, Rios M, Camargo YC. Geographic distribution of Leishmania species isolated from patients at the National Institute of Dermatology Federico Lleras Acosta E.S.E., 1995–2005. Biomedica. 2006;26 Suppl 1:145–51.
Blum J, Desjeux P, Schwartz E, Beck B, Hatz C. Treatment of cutaneous leishmaniasis among travellers. J Antimicrob Chemother. 2004;532:158–66.
Kato H, Caceres AG, Gomez EA, Mimori T, Uezato H, Marco JD, et al. Molecular mass screening to incriminate sand fly vectors of Andean-type cutaneous leishmaniasis in Ecuador and Peru. Am J Trop Med Hyg. 2008;795:719–21.
Sobrino R, Ferroglio E, Oleaga A, Romano A, Millan J, Revilla M, et al. Characterization of widespread canine leishmaniasis among wild carnivores from Spain. Vet Parasitol. 2008;1553–4:198–203.
Reithinger R, Dujardin JC. Molecular diagnosis of leishmaniasis: current status and future applications. J Clin Microbiol. 2007;451:21–5.
Bañuls A-L, Hide M, Prugnolle F. Leishmania and the leishmaniases: a parasite genetic update and advances in taxonomy, epidemiology and pathogenicity in humans. Adv Parasitol. 2007;64:1–458.
Zemanova E, Jirku M, Mauricio IL, Horak A, Miles MA, Lukes 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;372:149–60.
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;367:757–69.
Boite MC, Mauricio IL, Miles MA, Cupolillo E. New insights on taxonomy, phylogeny and population genetics of Leishmania (Viannia) parasites based on multilocus sequence analysis. PLoS Negl Trop Dis. 2012;611, e1888.
Marlow MA, Boite MC, Ferreira GE, Steindel M, Cupolillo E. Multilocus sequence analysis for Leishmania braziliensis outbreak investigation. PLoS Negl Trop Dis. 2014;82, e2695.
Odds FC. Molecular phylogenetics and epidemiology of Candida albicans. Future Microbiol. 2010;51:67–79.
Naser SM, Thompson FL, Hoste B, Gevers D, Dawyndt P, Vancanneyt M, et al. Application of multilocus sequence analysis (MLSA) for rapid identification of Enterococcus species based on rpoA and pheS genes. Microbiology. 2005;151Pt 7:2141–50.
Martens M, Dawyndt P, Coopman R, Gillis M, De Vos P, Willems A. Advantages of multilocus sequence analysis for taxonomic studies: a case study using 10 housekeeping genes in the genus Ensifer (including former Sinorhizobium). Int J Syst Evol Microbiol. 2008;58Pt 1:200–14.
Marco JD, Barroso PA, Locatelli FM, Cajal SP, Hoyos CL, Nevot MC, et al. Multilocus sequence typing approach for a broader range of species of Leishmania genus: describing parasite diversity in Argentina. Inf Genet Evol. 2015;30:308–17.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;2810:2731–9.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;2511:1451–2.
Tomasini N, Lauthier JJ, Llewellyn MS, Diosque P. MLSTest: novel software for multi-locus sequence data analysis in eukaryotic organisms. Infect Genet Evol. 2013;20:188–96.
Schonian G, Mauricio I, Cupolillo E. Is it time to revise the nomenclature of Leishmania? Trends Parasitol. 2010;2610:466–9.
Van der Auwera G, Fraga J, Montalvo AM, Dujardin JC. Leishmania taxonomy up for promotion? Trends Parasitol. 2011;272:49–50.
Downing T, Imamura H, Decuypere S, Clark TG, Coombs GH, Cotton JA, et al. Whole genome sequencing of multiple Leishmania donovani clinical isolates provides insights into population structure and mechanisms of drug resistance. Genome Res. 2011;2112:2143–56.
Van der Auwera G, Maes I, De Doncker S, Ravel C, Cnops L, Van Esbroeck M, et al. Heat-shock protein 70 gene sequencing for Leishmania species typing in European tropical infectious disease clinics. Euro Surveill. 2013;1830:20543.
Lukes J, Mauricio IL, Schonian G, Dujardin JC, Soteriadou K, Dedet JP, et al. Evolutionary and geographical history of the Leishmania donovani complex with a revision of current taxonomy. Proc Natl Acad Sci USA. 2007;10422:9375–80.
Foulet F, Botterel F, Buffet P, Morizot G, Rivollet D, Deniau M, et al. Detection and identification of Leishmania species from clinical specimens by using a real-time PCR assay and sequencing of the cytochrome B gene. J Clin Microbiol. 2007;457:2110–5.
Asato Y, Oshiro M, Myint CK, Yamamoto Y, Kato H, Marco JD, et al. Phylogenic analysis of the genus Leishmania by cytochrome b gene sequencing. Exp Parasitol. 2009;1214:352–61.
Marco JD, Uezato H, Mimori T, Barroso PA, Korenaga M, Nonaka S, et al. Are cytochrome b gene sequencing and polymorphism-specific polymerase chain reaction as reliable as multilocus enzyme electrophoresis for identifying Leishmania spp. from Argentina? Am J Trop Med Hyg. 2006;752:256–60.
Kato H, Gomez EA, Martini-Robles L, Muzzio J, Velez L, Calvopina M, et al. Geographic distribution of Leishmania species in Ecuador based on the cytochrome b gene sequence analysis. PLoS Negl Trop Dis. 2016;107, e0004844.
Britto C, Ravel C, Bastien P, Blaineau C, Pages M, Dedet JP, et al. Conserved linkage groups associated with large-scale chromosomal rearrangements between Old World and New World Leishmania genomes. Gene. 1998;2221:107–17.
Peacock CS, Seeger K, Harris D, Murphy L, Ruiz JC, Quail MA, et al. Comparative genomic analysis of three Leishmania species that cause diverse human disease. Nat Genet. 2007;397:839–47.
Ferro C, Lopez M, Fuya P, Lugo L, Cordovez JM, Gonzalez C. Spatial distribution of sand fly vectors and eco-epidemiology of cutaneous leishmaniasis transmission in Colombia. PLoS One. 2015;1010, e0139391.
Akhoundi M, Kuhls K, Cannet A, Votypka 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;103, e0004349.
Akopyants NS, Kimblin N, Secundino N, Patrick R, Peters N, Lawyer P, et al. Demonstration of genetic exchange during cyclical development of Leishmania in the sand fly vector. Science. 2009;324:265–8.
Tibayrenc M, Ayala FJ. How clonal are Trypanosoma and Leishmania? Trends Parasitol. 2013;29:264–9.
Rogers MB, Downing T, Smith BA, Imamura H, Sanders M, Svobodova M, et al. Genomic confirmation of hybridisation and recent inbreeding in a vector-isolated Leishmania population. PLoS Genet. 2014;10, e1004092.
Kuhls K, Cupolillo E, Silva SO, Schweynoch C, Boite MC, Mello MN, et al. Population structure and evidence for both clonality and recombination among Brazilian strains of the subgenus Leishmania (Viannia). PLoS Negl Trop Dis. 2013;7, e2490.
We thank Camila González from Universidad de Los Andes for the construction of the maps and Jorge Cardenas for his help in the experiments.
We thank the Departamento Administrativo de Ciencia, Tecnología e Innovación ‘COLCIENCIAS’ for funding the Project “Fortalecimiento de la capacidad diagnóstica, de investigación y de vigilancia de enfermedades transmisibles emergentes y reemergentes en Colombia” grant number 757–13 and funding GH as “Joven Investigador” convocatoria 706 del 2015.
Availability of data and materials
All the data are included in the article and its additional files.
JDR, GH, CH and AT analyzed the data and wrote the manuscript. MA and CF analyzed the data and collected the samples. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Leishmania parasites were isolated from human patients performing methods that were carried out in accordance with the approved guidelines by the ethics and technical scientific committee of National Institute of Health in Bogotá, Colombia (All experimental protocols were approved by this institute's committee). Written informed consent was obtained from all subjects. All subjects that provided written informed consent were adults.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Comparison of haplotype diversity (Hd) by species. a L. panamensis. b L. braziliensis. (PPTX 126 kb)
List of allelic profile by each gene fragment for the assignment of diploid sequence types in L. braziliensis. Table S2. List of allelic profile by each gene fragment for the assignment of diploid sequence types in L. panamensis. (XLSX 14 kb)