- Research article
- Open Access
RAPD-PCR reveals genetic polymorphism among Leishmania major strains from Tunisian patients
BMC Infectious Diseases volume 15, Article number: 269 (2015)
Zoonotic cutaneous leishmaniasis caused by Leishmania (L.) major is endemoepidemic in the Center and South of Tunisia. The clinical course of the disease varies widely among different patients and geographic regions. Although genetic diversity in L. major parasites has been suggested as a potential factor influencing their pathogenic variability, little information on genetic polymorphism among L. major strains is available in the literature. This work aimed to estimate the genetic variability within different isolates of L. major.
Our sample comprised 39 isolates (confirmed as L. major by restriction fragment length polymorphism typing) from patients experiencing the same clinical manifestations but living in different regions of Tunisia where L. major is endemic. Random amplified polymorphic DNA (RAPD) PCR marker polymorphism was estimated by calculating Nei and Li’s genetic distances and by an analysis of molecular variance (AMOVA).
Analysis of the genetic diversity among the isolates revealed a high level of polymorphism (43 %) among them. AMOVA indicated that the highest variability (99 %) existed within the study regions.
Our results revealed a heterogeneous genetic profile for L. major with similar clinical manifestations occurring within the different geographical regions. Additional L. major isolates from patients, insect vectors, and reservoir hosts from different endemic foci should be collected for further analysis.
Zoonotic cutaneous leishmaniasis (ZCL) caused by L. major is a vector borne disease. In Tunisia, the disease is distributed mainly in the central and southern parts of the country where it is considered a major public health problem.
Various clinical presentations have been reported in patients, varying from asymptomatic infections to ulcerative lesions of the skin, which last for many months and may cause disfigurement and stigma. The clinical pleomorphism of ZCL is associated with multiple factors such as variability in the species or subspecies of Leishmania , host genetic factors or immune status , host-parasite interactions , factors related to environmental upheaval depending on the distribution and density of the respective vectors and reservoirs hosts, transmission pressure, parasite density  and the presence of other pathogens . However, because all Leishmania parasites (including L. major) are vector-borne obligate endoparasites, adaptation to the internal microenvironment of their long-term residence in the macrophages of their mammalian reservoir hosts, is considered the most significant factor influencing the evolution of this genus . Several studies have been conducted to identify and type naturally occurring Leishmania parasites to attempt to correlate their genetic diversity with their clinical features such as virulence, pathogenicity and drug resistance. This information is crucial for showing the biology of this parasite [6–8].
The present study aimed to investigate genetic variability among L. major isolates from patients who displayed similar clinical features. To achieve this we used the genetic markers generated by the random amplified polymorphic DNA (RAPD) technique.
Parasite strains and DNA extraction
In total, 39 isolates were used in this study. Isolates were collected between 2007 and 2009 from areas in central and southern Tunisia where ZCL is endemic. Each isolate was obtained from a different human subject who had a typical nodulo-ulcerative skin lesion that had healed within few months of appearance after using only local care. Only individuals aged between 5 and 65 years who gave their written informed consent (or their parents or legal guardians consent in case of minors) were enrolled.
The Leishmania isolates were divided into five geographical groups comprising those from the Sidi Bouzid Governorate, the Kairouan Governorate and the Gafsa Governorate, and represented by the three delegations, M’dhila, Metlaoui and Gafsa Center (Table 1). Species identification of all the isolates as L. major was achieved by restriction fragment length polymorphism (RFLP) analysis. Isoenzyme typing of one isolate (MHOM/TN/2009/S566) was conducted at the Leishmania Reference Center of Montpellier, France. This isolate exhibited the MON-25 zymodeme, the most common zymodeme of L. major in Tunisia  and was considered a L. major reference strain.
For genomic DNA extraction, promastigotes corresponding to each isolate were harvested on the sixth day of culture in RPMI 1640 medium supplemented with 20 % heat inactivated fetal bovine serum (Invitrogen, Carlsbad, CA, USA), 2 mM glutamine, 100 U/ml penicillin/50 and U/ml streptomycin at 26 °C. A pellet of 2 × 108 parasites of each isolate was stored at −20 °C. DNA was extracted from each pellet using the QIAamp® DNA Mini Kit (QIAGEN, Germany) according to the manufacturer’s instructions. DNA concentrations were estimated by spectrophotometry by reading the absorbance at 260 nm. Scanning spectrophotometry and agarose gel electrophoresis were used to judge the purity and integrity of the DNA.
The 16 decameric primers (Invitrogen) described previously [10, 11] were used for this study. All primers used were resuspended in TE buffer, stored at −20 °C, and 10 mM (10pmol/ml) working solutions were prepared. The same person did the RAPD screening and amplification using a single protocol.
Four primers yielded patterns that were nondistinct for all of the isolates and were, therefore, excluded from the study. The excluded primers were A1 (CAGGCCCTTC), AB1-12 (CCTTGACGCA), AB1-14 (TTCCCCCGCT) and AB1-18 (CCACAGCAGT). Hence, the results presented herein are from only 12 of the primers (Table 2).
The RAPD mixtures were processed in 50 μl reactions containing 0.2 mM of each dNTP, 0.3 μM primer, 1 U of Taq DNA polymerase (Invitrogen), 10 x buffer (supplied with the enzyme; 10 mM Tris-HCl, 1.5 mM MgCl2, 50 mM KCl, pH 8.3), and 20 ng of genomic DNA. Amplification was performed in a thermocycler (PXE 0.5 Thermal Cycler, Thermo Electron Corporation, Waltham, USA) using 1 cycle at 94 °C for 2 min, followed by 35 cycles at 94 °C for 1 min, 36 °C for 1 min, 72 °C for 2 min, and final extension at 72 °C for 10 min. A positive control (containing ADN extracted from Leishmania culture reference strains MHOM/TN/2009/S566) and a negative control to detect contamination (containing water and no DNA) were included in each PCR run. The amplified products were analyzed on 2 % agarose gels run in 1 × TBE buffer stained with SYBR Safe (Invitrogen), visualized under UV light, and photographed (GenoView, VWR, France). A 1 kb DNA ladder (Invitrogen) was used as a molecular marker.
Genetic similarity between the isolates was evaluated through simple association, and Nei and Li (1979) genetic distances were calculated using RAPD Distance 1.04 software . The binary matrix was built pair-wise, and the presence or absence of a RAPD band was scored 1 or 0, respectively [13, 14]. The genetic distance dendrogram was constructed using MVSP software package (Multi Variate Statistical Package 3.1) . Analysis of molecular variance (AMOVA) using WINAMOVA 1.55 , with the significance levels of the variance components based on 1000 permutations, was applied to quantify the degree of variation between isolates from the same or different collection sites.
Ethical clearance for this research study was obtained from the ethical committee of the Pasteur Institute of Tunis, Tunisia.
The genetic variability of 39 L. major isolates from five endemic foci for L. major (Sidi Bouzid, Kairouan and Gafsa represented by the three delegations, Mdhila, Metlaoui and Gafsa Center) was investigated by RAPD. The consistency of the amplified DNA bands at the same or a different electrophoretic position for L. major from the different endemic areas as well as from the same endemic area were analyzed. A total of 53 RAPD fragments were generated whose size ranged from 0.1 to 2.3 kb. The reproducibility of the results was confirmed by repeating the RAPD-PCR analysis with 39 randomly selected isolates. The profile of the isolates in the reproducibility study was similar to their profiles in the original study. The number of reproducible bands produced by each primer ranged from 1 to 10 with an average of 4.5 polymorphic bands per primer.
The RAPD profiles are illustrated in Figs. 1a and b. The greatest number of bands (10) was produced by primer 329, while the A5 primer produced only one band. The resulting band pattern shows that the S124 isolate from Metlaoui had a completely distinct profile from the other isolates, yielding two major products of 0.1–2.3 kb, thereby differentiating it from L. major DNA.
The polymorphic band percentages ranged from 16 to a maximum of 68 with an average polymorphism level of 43 %. The polymorphism levels varied according to the primers used and the geographic origin of the isolates (Table 3). The most polymorphic RAPD markers were identified with AB01-01 and 329 primers with a percentage of 68 %, while the AB1-09 primer produced the lowest polymorphism (16 %). On average, for all primers, isolates from Metlaoui produced the highest percentage (91 %) of polymorphic bands, while the lowest percentage was seen among the Kairouan isolates (23 %).
The genetic identity between the pairs of individuals ranged from 0.59 to 1.00, as based on the criteria of Nei and Li (1979). The dendrogram constructed from Nei’s genetic distances showed two obvious clusters (Fig. 2). The S124 strain belonging to the Metlaoui area appeared on its own, while all the other strains formed a cluster. The major cluster was divided further into two branches with one consisting of the S004 strain from M’dhila while the second divided the remaining strains into several sub-branches. We did not observe any specific clusters that corresponded with the geographical origin of isolates.
The hierarchical analysis of the molecular diversity in the isolates using AMOVA was performed with genetic distances between pairs to analyze the partitioning of RAPD variation in L. major isolates among and within populations belonging to the different areas (Table 4).
AMOVA analysis within and among the study areas revealed that 98.9 % of the total variation resided between isolates within foci, and only 1.1 % within isolates between foci (Table 4).
Despite the fairly low structuring between foci, the phi-statistics (ΦST) between the population pairs were variable. The largest ΦST value (0.189) was observed between the Gafsa Center and Kairouan regions, while the lowest value (0.006) was between M’dhila and Kairouan.
ZCL caused by L. major in Tunisia has various clinical manifestations. Herein, we used RAPD-PCR to study genetic diversity among 39 L. major isolates from the central and southern parts of Tunisia where the parasite is endemic. The twelve polymorphic primers we used revealed a high level of polymorphism (43 %) among the L. major isolates collected from patients with similar clinical manifestations of the disease.
Despite many studies suggesting that the genetic variability observed among L. major isolates is correlated with different clinical manifestations [6, 7, 17] or responses to antileishmanial treatments , our findings suggest that genetic diversity can also be found among L. major isolates that cause the same clinical picture. Obviously, the severity of the disease is influenced by the pathogenesis of the parasite species or strain and the host immune responses against it . In fact, a series of complex and not fully understood interactions that occur between species-specific virulence factors and genetically determined cell-mediated immunity contributes to disease susceptibility, as shown in murine models [2, 19–21], but the contribution of each component within a particular region is unknown in humans. Therefore, this raises the question of whether variation in the pathogenicity of the parasite and the clinical pattern of ZCL may be related to genetic heterogeneity in the parasite. However, to address this question it is first necessary to characterize genetic polymorphism in the parasite.
Furthermore, considerable heterogeneity was detected among isolates belonging to the same endemic area, with percentages ranging from 23 to 91 %. In agreement with this result, a previous study in Iran showed substantial variability between some isolates from the same regions using the same RAPD technique . This result might be explained by considerable gene flow among isolates belonging to the same area. Genetic exchange may result from the relatively old evolutionary history of L. major in the study sites. Additionally, the zoonotic transmission cycles and the number of reservoir hosts and vectors are known to play an important role in sustaining genetic diversity. Transmission of L. major in Tunisia is associated with one major vector (Phlebotomus papatasi), but the existence of different species of wild rodent (Psammomys obesus, Meriones shawi, Meriones libycus) serving as animal reservoir hosts, and sometimes their coexistence within a small geographic area seems to support this hypothesis . L. major strains isolated from both sandflies and rodents should be analyzed for a reliable characterization of L. major Tunisian strains.
Although our study isolates from within the same area exhibited genetic polymorphism, lower genetic diversity was found between populations of L. major based on their collection sites. This was first noted in the dendrogram built from the RAPD-PCR data based on Nei’s genetic distances. The dendrogram showed no significant clustering of the isolates according to their geographical origin, a finding that was confirmed later by the AMOVA analysis, which indicated that only 1.1 % of the observed variation occurred among the groups, while 98.9 % occurred among isolates within the groups. This finding most likely reflects the unbalanced and low number of isolates per area sampled rather than a lack of correlation between genetic polymorphism and the geographical origin of the parasite isolates, as has been noted before . Investigating this type of correlation was beyond the scope of the present work, and should be tested using an appropriate sampling strategy with a larger sample size per group to reveal potential genetic diversity. The underestimation of genetic variability between regions caused by sampling differences was also supported by the fact that the highest polymorphism percentages were obtained within the Metlaoui and M’dhila regions, which included the highest numbers of isolates among the five study regions.
Some authors have failed to establish an association between genetic variability in their isolates and their geographic distribution when too few strains and/or poorly discriminating markers were used, while others explained the lack of focus-specific genotypes by the small size of the country and the frequent migration between foci [25, 26].
However, other studies have identified geographic structuring of L. major populations using the RAPD technique , or other PCR-based methods targeting different nuclear multi-copy sequences or antigen-coding genes, followed by RFLP, permissively primed intergenic polymorphic-polymerase chain reaction, single-stranded conformation polymorphism analysis, or multilocus microsatellite typing [25, 27, 28]. Geographical differences have been attributed by some authors to the existence of different populations of vectors and/or to different species of rodent hosts [25, 27–29]. Thus, Leishmania strains isolated from different stages of the parasite life cycle belonging to the same or different geographic origins need to be investigated to explain the mechanisms involved in the evolution of L. major.
Both AB01-01 and 329 primers produced considerable genetic variability between L. major isolates from different regions of Tunisia as well as within the same regions. These primers were chosen because of their reproducibility and ability to reveal intraspecies variation.
The RAPD technique was used for identification and differentiation of old world species of Leishmania [10, 25, 27–29]. For L. major, these authors reported that the RAPD technique amplified a 0.45 kb product with the A7 primer and generated a 2.7–1.9 kb product with the A10 primer. Nevertheless, our results suggest the validity of the RAPD technique when used for the identification and characterization of Leishmania species, especially L. major. The RAPD technique amplified a 0.2–2.3 kb product with the A7 primer and generated a 0.2–1 kb product with the A10 primer. The resulting band pattern shows that the S124 isolate from Metlaoui produced different amplicons from the primers.
The first case of cutaneous leishmaniasis in Tunisia was described more than 120 years ago  in the region of Gafsa. Since this date, this province has been known as an active focus of ZCL. However, in 2012, twenty-four human cases of cutaneous leishmaniasis caused by L. killicki were reported for the first time in Metlaoui district . This important set of data about the new L. killicki focus in Metlaoui, and the mixed transmission cycle in this province, allows speculation for S124 about the existence of hybrids and gene flow between genetically different populations.
This is the first time that such a large number of this species isolated from humans have been included from different endemic foci in Tunisia. Here, we have used the RAPD technique as a tool to identify genetic variation within Leishmania species. Even though our results have shown that DNA fingerprinting and RAPD analysis were able to discriminate between different L major isolates of closely related strains, highly polymorphic markers using fast evolving microsatellites need to be developed to achieve the highest discrimination level among strains. According to these markers, microsatellites are hypervariable, genetically neutral, co-dominant and are, therefore, ideal to reliably address the role that genetic polymorphism plays in the epidemiology of Leishmania disease at a fine scale .
Our results suggest that isolates of Tunisian L. major are highly polymorphic, despite inducing similar clinical manifestations and belonging to the same geographical origin. This finding requires confirmation using a method with a higher discriminatory capacity and with larger sample size to decipher more fully the pathogenesis of the parasite and show the mechanisms by which it circulates in a restricted geographical area. In future research it would also be worth collecting samples from patients, vectors, and reservoir hosts in the main areas where the parasite is endemic to help improve our understanding of the molecular epidemiology of ZCL in Tunisia.
Zoonotic Cutaneous Leishmaniasis
Random amplified polymorphic DNA
Analysis of molecular variance
Masmoudi A, Ayadi N, Boudaya S, Meziou TJ, Mseddi M, Marrekchi S, et al. [Clinical polymorphism of cutaneous leishmaniasis in centre and south of Tunisia]. Bull Soc Pathol Exot. 2007;100:36–40.
Solbach W, Laskay T. The host response to Leishmania infection. AdvImmunol. 2000;74:275–317.
Ashford RW. The leishmaniases as emerging and reemerging zoonoses. Int J Parasitol. 2000;30:1269–81.
Eswarappa SM, Estrela S, Brown SP. Within-host dynamics of multi-species infections: facilitation, competition and virulence. PLoS One. 2012;7:e38730.
Waki K, Dutta S, Ray D, Kolli BK, Akman L, Kawazu S, et al. Transmembrane molecules for phylogenetic analyses of pathogenic protists: Leishmania-specific informative sites in hydrophilic loops of trans- endoplasmic reticulum N-acetylglucosamine-1-phosphate transferase. Eukaryot Cell. 2007;6:198–210.
Baghaei M. Intraspecific variation in Leishmania major isolated from different forms of zoonotic Cutaneous Leishmaniasis. Iran J Med Sci. 2005;30:51–4.
Oryan A, Shirian S, Tabandeh MR, Hatam GR, Randau G, Daneshbod Y. Genetic diversity of Leishmania major strains isolated from different clinical forms of cutaneous leishmaniasis in southern Iran based on minicircle kDNA. Infect Genet Evol. 2013;19:226–31.
Cupolillo E, Momen H, Grimaldi Jr G. Genetic diversity in natural populations of New World Leishmania. Mem Inst Oswaldo Cruz. 1998;93:663–8.
Haouas N, Gorcii M, Chargui N, Aoun K, Bouratbine A, Messaadi Akrout F, et al. Leishmaniasis in central and southern Tunisia: current geographical distribution of zymodemes. Parasite. 2007;14:239–46.
Ikram G, Dellagi K, Ismail RB. Random amplified polymorphic DNA technique for identification and differentiation of old world Leishmania species. Am J Trop Med Hyg. 2002;66:152–6.
Hajjaran H, Mohebali M, Razavi MR, Rezaei S, Kazemi B, Edrissian GH, et al. Identification of Leishmania species isolated from human cutaneous leishmaniasis, using random amplified polymorphic DNA, (RAPD-PCR). Iranian Journal of Public Health. 2004;33:8–15.
Armstrong JS, Gibbs AJ, Peakall R, Weiler G. The RAPDistance package, ver. 1.03. Research School of Biological Sciences, Australian National University: Canberra, Australia; 1994.
Mattnera J, Zawkoa G, Rossettoc M, Kraussb SL, Dixonb KW, Sivasithamparama K. Conservation genetics and implications for restoration of Hemigenia exilis (Lamiaceae), a serpentine endemic from Western Australia. Biol Conserv. 2002;107:37–45.
Mekuria G, Collins G, Sedgley M. Genetic diversity within an isolated olive (Olea europaea L.) population in relation to feral spread. Sci Hort. 2002;94:91–105.
Kovach W. MVSP-Multivariate Statistical Package (software) version 3.1. Pentraeth, UK: Kovach Computing Services; 1999.
Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131:479–91.
van Thiel PP, van Gool T, Faber WR, Leenstra T, Kager PA, Bart A. Variation in clinical presentation and genotype of causative Leishmania major strain in cutaneous leishmaniasis in north and south Afghanistan. Am J Trop Med Hyg. 2011;85:60–3.
Garcia AL, Kindt A, Quispe-Tintaya KW, Bermudez H, Llanos A, Arevalo J, et al. American tegumentary leishmaniasis: antigen-gene polymorphism, taxonomy and clinical pleomorphism. Infect Genet Evol. 2005;5:109–16.
Blackwell JM. Genetic susceptibility to leishmanial infections: studies in mice and man. Parasitology. 1996;112(Suppl):S67–74.
Reiner SL, Locksley RM. The regulation of immunity to Leishmania major. Annu Rev Immunol. 1995;13:151–77.
Russo DM, Barral-Netto M, Barral A, Reed SG. Human T-cell responses in Leishmania infections. Prog Clin Parasitol. 1993;3:119–44.
Mahmoudzadeh-Niknam H, Ajdary S, Riazi-Rad F, Mirzadegan E, Rezaeian A, Khaze V, et al. Molecular epidemiology of cutaneous leishmaniasis and heterogeneity of Leishmania major strains in Iran. Trop Med Int Health. 2012;17:1335–44.
Ghawar W, Toumi A, Snoussi MA, Chlif S, Zaatour A, Boukthir A, et al. Leishmania major infection among Psammomys obesus and Meriones shawi: reservoirs of zoonotic cutaneous leishmaniasis in Sidi Bouzid(central Tunisia). Vector Borne Zoonotic Dis. 2011;11:1561–8.
Fitzpatrick BM. Power and sample size for nested analysis of molecular variance. Mol Ecol. 2009;18:3961–6.
Al-Jawabreh A, Diezmann S, Muller M, Wirth T, Schnur LF, Strelkova MV, et al. Identification of geographically distributed sub-populations of Leishmania (Leishmania) major by microsatellite analysis. BMC Evol Biol. 2008;8:183.
Seridi N, Amro A, Kuhls K, Belkaid M, Zidane C, Al-Jawabreh A, et al. Genetic polymorphism of Algerian Leishmania infantum strains revealed by multilocus microsatellite analysis. Microbes Infect. 2008;10:1309–15.
Mahnaz T, Al-Jawabreh A, Kuhls K, Schonian G. Multilocus microsatellite typing shows three different genetic clusters of Leishmania major in Iran. Microbes Infect. 2011;13:937–42.
Hamarsheh O, Presber W, Yaghoobi-Ershadi MR, Amro A, Al-Jawabreh A, Sawalha S, et al. Population structure and geographical subdivision of the Leishmania major vector Phlebotomus papatasi as revealed by microsatellite variation. Med Vet Entomol. 2009;23:69–77.
Elfari M, Schnur LF, Strelkova MV, Eisenberger CL, Jacobson RL, Greenblatt CL, et al. Genetic and biological diversity among populations of Leishmania major from Central Asia, the Middle East and Africa. Microbes Infect. 2005;7:93–103.
Deperet C, Boinet E. Du bouton de Gafsa au camp de Sathonay. Archives Médicales et Pharmaceutiques Militaires. 1884;3:296–302.
Jaouadi K, Depaquit J, Haouas N, Chaara D, Gorcii M, Chargui N, et al. Twenty-four new human cases of cutaneous leishmaniasis due to Leishmania killicki in Metlaoui, southwestern Tunisia. Probable role of Phlebotomus sergenti in the transmission. Acta Trop. 2012;122:276–83.
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.
We thank Dr Francine Pratlong from the National Reference Center for Leishmaniasis of Montpellier for her technical contribution to typing the Leishmania parasites. We also thank Dr Chokri Massoud from INSAT for help in statistical analysis. The work was supported by the Laboratory of Transmission, Control and Immunobiology of Infections (LR11IPT02), Pasteur Institute of Tunis.
The authors declare that they have no competing interests.
RY carried out the molecular genetic studies, JB performed the statistical analysis and drafted the manuscript, WG participated in discussions about the analysis and helped to draft the manuscript, KJ and SC participated in the design and molecular analysis of the study, AZ and ABS conceived and coordinated the study. All authors have read and approved the final manuscript.
Rihab Yazidi and Jihene Bettaieb contributed equally to this work.
About this article
Cite this article
Yazidi, R., Bettaieb, J., Ghawar, W. et al. RAPD-PCR reveals genetic polymorphism among Leishmania major strains from Tunisian patients. BMC Infect Dis 15, 269 (2015). https://doi.org/10.1186/s12879-015-1010-0
- Leishmania major
- Genetic Polymorphism