Environmental risk modelling and potential sand fly vectors of cutaneous leishmaniasis in Chitral district: a leishmanial focal point of mount Tirich Mir, Pakistan

Objective To provide baseline information about suspected vectors and the incidence, distribution and an active zone of transmission for cutaneous leishmaniasis (CL) in Chitral, Pakistan, using GIS tools; and to investigate the role of environmental factors in the disease dynamics. Method Two surveys in 2014 and 2016 as a basis for choropleth and environmental risk mapping. Results A total of 769 captured specimens yielded 14 Phlebotomus and six Sergentomyia species including two potential vectors of CL, i.e. Phlebotomus papatasi and Phlebotomus sergenti. P. papatasi (71%) was dominant, followed by P. sergenti (18%). A choropleth map generated in Arcmap 10.1 based on 1560 CL case reports displayed maximum prevalence (0.92–2.5%) in Ayun, Broz, Charun, Chitral 1 and 2 and Darosh 1 and 2 union councils. An environmental risk map constructed by MaxEnt 3.3.3 defined an active zone of transmission based on leishmaniasis occurrence records (n = 315). The analysis of variable contribution in MaxEnt indicates significance of elevation (54.4%), population density (23.3%) and land use/land cover (6.6%) in CL disease dynamics. Conclusion The probability of CL increases (0.6–1 on logistic scale) in severely deforested areas, in lowland valleys and in regions with high‐population density. Objectif Fournir des informations de base sur les vecteurs soupçonnés et sur l'incidence, la distribution et une zone active de transmission de la Leishmaniose cutanée (LC) à Chitral, au Pakistan, en utilisant des outils du SIG; et étudier le rôle des facteurs environnementaux dans la dynamique de la maladie. Méthode Deux enquêtes ont été menées en 2014 et 2016 comme base pour la carte chloroplète et les risques environnementaux. Résultats 769 spécimens capturés ont révélé 18 espèces Phlebotomus et 4 Sergentomyia incluant deux vecteurs potentiels de la LC, c'est‐à‐dire Phlebotomus papatasi et Phlebotomus sergenti. P. papatasi (71%) était dominant, suivi de P. sergenti (18%). Une carte chloroplète générée dans Arcmap 10.1 sur base de 1.560 rapports de cas de LC a montré une prévalence maximale (0,92 à 2,5%) dans les conseils syndicaux d'Ayun, Broz, Charun, Chital 1 et 2 et Darosh 1 et 2. Une carte de risque environnementale générée sur MaxEnt 3.3.3 a défini une zone de transmission active basée sur des cas enregistrés d'apparition de la leishmaniose (n = 315). L'analyse de la contribution des variables dans MaxEnt indique une importance de l’élévation (54,4%), de la densité de populatio