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

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.

Keywords: Phlebotomus papatasi; Expressed sequence tags; Microsatellites; Population structure analysis.

Fig. 1

Study locations where P. papatasi sand flies were collected and provided for this project

Fig. 2

Analysis of the population structure of 302 P. papatasi sand fly individuals from 11 countries and genotyped using 14 microsatellites. A The output of PopCluster run from K = 1 to K = 9. B The optimal K which is represented by a peak of DLK2 over K = 5. C The total number of individuals in y axes assigned to P. papatasi individuals from each country. D The total number of individual flies per country assigned to the resolved populations at K = 5