In Tunisia, cases of zoonotic cutaneous leishmaniasis caused by Leishmania major are increasing and spreading from the south-west to new areas in the center. To improve the current knowledge on L. major evolution and population dynamics, we performed multi-locus microsatellite typing of human isolates from Tunisian governorates where the disease is endemic (Gafsa, Kairouan and Sidi Bouzid governorates) and collected during two periods: 1991–1992 and 2008–2012. Analysis (F-statistics and Bayesian model-based approach) of the genotyping results of isolates collected in Sidi Bouzid in 1991–1992 and 2008–2012 shows that, over two decades, in the same area, Leishmania parasites evolved by generating genetically differentiated populations. The genetic patterns of 2008–2012 isolates from the three governorates indicate that L. major populations did not spread gradually from the south to the center of Tunisia, according to a geographical gradient, suggesting that human activities might be the source of the disease expansion. The genotype analysis also suggests previous (Bayesian model-based approach) and current (F-statistics) flows of genotypes between governorates and districts. Human activities as well as reservoir dynamics and the effects of environmental changes could explain how the disease progresses. This study provides new insights into the evolution and spread of L. major in Tunisia that might improve our understanding of the parasite flow between geographically and temporally distinct populations.
Author Summary In Tunisia, zoonotic cutaneous leishmaniasis (ZCL) constitutes a significant public health problem. Since 1884, the Gafsa, Kairouan and Sidi Bouzid governorates are the most endemic areas of ZCL. This study used a multi-locus microsatellite typing approach to study the evolution and the population dynamics of Leishmania major in Tunisia. Within the same area, in twenty years, parasite populations evolved by producing a genetically differentiated population, probably better adapted to the ecosystem. In agreement with the reported human cases of ZCL, the genetic data on samples from the three governorates shows that the disease did not spread according to a geographical gradient. Furthermore, L. major flows seem to still occur between governorates and neighboring districts. This study suggests that environmental changes, human activities and reservoir systems have influenced the spread and evolution of L. major populations. Our findings provide important knowledge on the epidemiology of L. major in Tunisia and might help understanding why the disease is still spreading from the south to the center, despite the control measures that have been put into place.