Phylogenetic structure of Leishmania tropica in the new endemic focus Birjand in East Iran in comparison to other Iranian endemic regions

Iran has been identified being among the countries with the highest number of cutaneous leishmaniasis (CL) cases. South Khorasan province in East Iran is an emerging focus of CL. Species identification of sixty clinical samples by ITS1 PCR-RFLP presented evidence for the dominance of Leishmania tropica (90%) in this region. Analysis of the ITS1 sequence of 19 L. tropica isolates revealed seven closely related sequence types. In addition, ITS1 sequences available in GenBank from other Iranian regions were compiled for comparison with the studied isolates. Iranian L. tropica was distributed in two main clusters. All East Iranian sequence types were grouped with strains from foci from Southeast and Central regions in cluster A, showing highly similar sequences. The highest similarity was observed between most L. tropica from East and all isolates from Southeast regions and from Savojbolagh county in Central Iran. Southwest L. tropica was shown to be paraphyletic as the isolates were distributed in both clusters A and B. All Northeastern L. tropica were part of cluster B, however they showed significant heterogeneity and were distributed in different subclusters. Distribution of L. tropica populations was to some extent congruent with genetic lineages of Phlebotomus sergenti in Iran and may be an evidence for parasite-vector co-evolution. Southeast-East L. tropica was also similar to strains from Herat province in Afghanistan at the East border of Iran. This is the first comprehensive study on population structure of L. tropica in Iran that provides a guideline for appropriate sampling for further molecular based epidemiological studies.