Population structure analysis of Phlebotomus papatasi populations using transcriptome microsatellites: possible implications for leishmaniasis control and vaccine development

Abstract

Abstract Background Phlebotomus 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. Methods In 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. Results The 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. Conclusions Elucidating 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 Abstract

Keywords

• Phlebotomus papatasi
• Microsatellites
• Population structure analysis
• Expressed sequence tags