BackgroundPhlebotomus papatasi is considered the primary vector of Leishmania major parasites that cause zoonotic cutaneous leishmaniasis (ZCL) in the Middle East and North Africa. Phlebotomus papatasi populations have been studied extensively, revealing the existence of different genetic populations and subpopulations over its large distribution range. Genetic diversity and population structure analysis using transcriptome microsatellite markers is important to uncover the vector distribution dynamics, essential for controlling ZCL in endemic areas.MethodsIn this study, we investigated the level of genetic variation using expressed sequence tag-derived simple sequence repeats (EST-SSRs) among field and colony P. papatasi samples collected from 25 different locations in 11 countries. A total of 302 P. papatasi sand fly individuals were analyzed, including at least 10 flies from each region.ResultsThe analysis revealed a high-level population structure expressed by five distinct populations A through E, with moderate genetic differentiation among all populations. These genetic differences in expressed genes may enable P. papatasi to adapt to different environmental conditions along its distribution range and likely affect dispersal.ConclusionsElucidating the population structuring of P. papatasi is essential to L. major containment efforts in endemic countries. Moreover, the level of genetic variation among these populations may improve our understanding of Leishmania–sand fly interactions and contribute to the efforts of vaccine development based on P. papatasi salivary proteins.Graphical