Next Article in Journal
Liver Injury in Favipiravir-Treated COVID-19 Patients: Retrospective Single-Center Cohort Study
Next Article in Special Issue
Ecological Niche Modelling Approaches: Challenges and Applications in Vector-Borne Diseases
Previous Article in Journal
Strongyloides and COVID-19: Challenges and Opportunities for Future Research
Previous Article in Special Issue
Epidemiological Survey on Tick-Borne Pathogens with Zoonotic Potential in Dog Populations of Southern Ethiopia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distribution and Risk of Cutaneous Leishmaniasis in Khyber Pakhtunkhwa, Pakistan

1
Department of Zoology, Abdul Wali Khan University, Mardan 23300, Khyber Pakhtunkhwa, Pakistan
2
Institute of Environmental and Occupational Health Sciences, College of Public Health, National Taiwan University, Taipei 100025, Taiwan
3
Facultad de Medicina Veterinaria y Zootecnia, Campus de Ciencias Biologicas y Agropecuarias, Universidad Autonoma de Yucatán, Km 15.5 Carretera Mérida–Xmatkuil, Merida 97100, Yucatan, Mexico
4
Department of Botany/Zoology, Bacha Khan University, Charsadda 24420, Khyber Pakhtunkhwa, Pakistan
5
Department of Public Health, College of Public Health, National Taiwan University, Taipei 100025, Taiwan
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Trop. Med. Infect. Dis. 2023, 8(2), 128; https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed8020128
Submission received: 18 January 2023 / Revised: 16 February 2023 / Accepted: 17 February 2023 / Published: 20 February 2023

Abstract

:
Cutaneous leishmaniasis (CL) is a zoonotic infection caused by obligate intracellular protozoa of the genus Leishmania. This study aimed to investigate CL in Khyber Pakhtunkhwa, Pakistan and to estimate the risk of epidemics. Clinico-epidemiological data of 3188 CL patients were collected from health facilities in 2021. Risk factors were analyzed using the chi-square test. ArcGIS V.10.7.1 was applied for spatial analysis. The association between CL occurrence and climatic variables was examined by Bayesian geostatistical analysis. The clinical data revealed males or individuals younger than 20 years old were more affected. Most patients presented with a single lesion, and the face was the most attacked body part. CL was prevalent in the southern region in winter. A proportional symbol map, a choropleth map, and a digital elevation model map were built to show the distribution of CL. Focal transmission was predicted by inverse distance weighting interpolation. Cluster and outlier analysis identified clusters in Bannu, Dir Lower, and Mardan, and hotspot analysis suggested Bannu as a high-risk foci. Bayesian geostatistical analysis indicated that increasing precipitation and temperature as well as low altitudes were associated with CL infection. The study has provided important information for public health sectors to develop intervention strategies for future CL epidemics.

1. Introduction

Leishmaniasis is a zoonotic infection caused by obligate intracellular protozoa of the genus Leishmania. Over 20 species of Leishmania parasites have been linked to diseases in humans. Natural transmission of leishmaniasis is carried out by sandflies of the genus Phlebotomus in the Old World and Lutzomyia in the New World [1,2]. Hyraxes, canids, and rodents are common reservoirs, but some species of Leishmania are found to be anthroponotic [3]. After being injected into the host’s skin, Leishmania escapes host immune responses by infecting and proliferating in phagocytes, especially macrophages, and results in varied clinical manifestations with diverse prognoses [4,5]. Four clinical forms have been categorized, namely visceral leishmaniasis (VL), cutaneous leishmaniasis (CL), mucocutaneous, and post kala-azar dermal leishmaniasis (PKDL). Leishmaniasis is endemic in many tropical and subtropical regions, affecting more than 12 million people in at least 88 countries and leading to a morbidity and mortality loss of 2.4 million disability-adjusted life-years (DALYs) and nearly 70,000 deaths [6,7,8]. According to the World Health Organization (WHO), the estimated annual incidence of CL and VL have ranged from 600,000 to 1 million and 50,000 to 90,000, respectively. About 90% of both cutaneous and visceral leishmaniases have been reported from developing countries such as Afghanistan, Bangladesh, Brazil, Saudi Arabia, Syria, and Peru [9].
CL is the most common form of leishmaniasis and is characterized by one or several ulcers or nodules developing at the sites of the infectious sandfly bites. CL has been known in human history for centuries. It was mentioned in the Old World in the first century and in Ecuador and Peru in the New World from 400–900 AD [10,11,12]. Illness resembling CL was called ‘Aleppo evil’ in the Middle East and ‘Dehli boil’ in the Indian subcontinent in the 16th and 17th centuries, when protozoa were identified in skin lesions [13,14]. Different names, including Rose of Jericho, forest yaws, Baghdad sore, Uta, and Chiclero’s ulcer have been used to describe CL lesions [14]. CL continued to spread via the aid of travel and habitat expansion. Risk factors such as urbanization, migration, population surge and displacement, drug resistance, human modification of the environment, and new agricultural practices have been associated with the emergence of new foci. Infection caused by new species has also been documented [15,16,17,18].
In Pakistan, approximately 21,700–35,700 cases of CL are reported annually, and epidemics have occurred in Punjab (Multan), Baluchistan, and Khyber Pakhtunkhwa [19]. CL has been widespread in communities with low-income and in neglected areas due to their limited access to health care [20]. Moreover, the situation of leishmaniasis in Pakistan has been exacerbated. Nowadays, both VL and CL have been reported in numerous regions of the country, including Khyber Pakhtunkhwa [21]. Located in the north-west, Khyber Pakhtunkhwa has been one of the provinces most affected by CL [3]. The inflow of Afghan refugees resulted in the spread of CL to formerly non-endemic areas [22]. CL has also become an issue faced by the peacekeeping forces and soldiers posted in Afghanistan and the Federally Administered Tribal Areas (FATA) of Pakistan [23,24].
Geographic information systems (GIS) is a computer-based tool which has been widely applied to health science research [25]. As global positioning system (GPS) mobile devices have become popular, GIS and remote sensing have been increasingly adopted to determine the spatial epidemiology of vector-borne diseases, including malaria, dengue fever, leishmaniasis, etc. [26,27,28,29]. The distribution and spread of vector-borne diseases largely depend on the environmental requirement of vectors. Temperature, rainfall, urbanization, vegetation, and human society are known factors influencing vector behaviors [30]. Commercial GIS tools such ArcGIS incorporating statistical models in the software can not only be used for disease mapping for risk assessment but also can be used to estimate vector abundance and incidence by integrating data of climatic, demographic, and environmental factors [31,32]. Hence, control programs can be designed and implemented targeting a specific region based on local vector dynamics and their interactions with the hosts [33,34]. However, sandfly activities have not been constantly monitored in Pakistan despite the fact that a few studies have reported the environmental risk factors affecting sandfly populations [35,36,37]. Bayesian geostatistical analysis would provide a flexible approach to the availability of data, demonstrated by the identification of critical variables for disease transmission of malaria, schistosomiasis, and leishmaniasis [38,39,40,41]. The probabilistic likelihood-based framework of the Bayesian approach allows assimilation of available information, and the use of Markov chain Monte Carlo methods for data augmentation offers the chance to treat missing data [42].
The epidemiology of CL is changing in Pakistan, especially in the border lands in Khyber Pakhtunkhwa, and CL has emerged as a public health concern in past decades. However, lack of disease surveillance and vector data have limited the effectiveness of control efforts. This study aimed to illustrate a risk map of CL in Khyber Pakhtunkhwa. The epidemiology of CL was investigated based on clinical data. The spatial distribution of cases was analyzed by ArcGIS software, and the transmission efficacy, active zones, and coming epidemics throughout the study area were explored. Bayesian geostatistical analysis was applied to determine the environmental risk factors. The results provide an overview of CL in Khyber Pakhtunkhwa and would be helpful in developing control strategies for the disease in the region.

2. Materials and Methods

2.1. Study Area

Khyber Pakhtunkhwa province, formerly known as the North-West Frontier Province (NWFP), has an area of 101,741 km2 and a population of about 35.5 million, comprising 50.7% males and 49.7% of females. Sharing a border with Afghanistan and tribal areas, it consists of 35 administrative districts (Figure 1). Earthquakes attack Khyber Pakhtunkhwa province frequently because of its location in the weak tectonic zone. In addition, floods from the Indus River happen during the monsoon season every year. This province has experienced many shattering floods in the last twenty years. Of 22 serious floods recorded during 1950 to 2014, the flood in 2010 was the most devastating, affecting millions of people and households in the province [8]. In other parts of the province, particularly the mountainous areas in the north, floods occur due to torrents, landslides, glacial lake outbursts, and rapid glacial run. Forestry, mining, and agriculture are the province’s primary economic activities. Approximately 78% of marble production in Pakistan is from the province. The gross domestic product (GDP) in 2021–2022 was estimated at 1071 US dollar per capita.

2.2. Ethical Statement

The study has been approved by the Ethical Committee of Chemical and Life Section, Department of Zoology, Abdul Wali Khan University Mardan, Pakistan under the approval no. AWKUM-20053310.

2.3. Data Collection and Processing

The logic and procedure of the study is shown in Supplementary Figure S1. The clinical-epidemiological data of microscopically confirmed CL patients were obtained from health facilities of health departments across the province throughout 2021. CL patients were diagnosed at local facilities and treated and followed-up systematically. The treatment courses of CL were excluded from this study due to incompleteness of information. Other data were entered into the spreadsheets of Microsoft Excel (Windows version, 2016) for descriptive analysis. Coordinates of the districts were obtained from Google Earth Pro (version 7.3), and the average incidence of CL in each district was aligned respectively. ArcGIS (version 10.7.1) was used for spatial risk and statistical analysis. A proportional symbol map was illustrated according to the cases of risk groups in districts. A digital elevation model (DEM) was generated using the geostatistics from the link http://www.diva-gis.org/datadown (accessed on 27 January 2023), which were light detection and ranging (LiDAR) and interferometric synthetic aperture radar (IfSAR, only for Alaska) data with a grid of 10 m2 resolution. The exact elevation of each district was obtained by computing DEM raster output data.

2.4. Choropleth Map, Inverse Distance Weighting (IDW) Interpolation

A choropleth map was built by representing the number of CL cases with continuous colors in lack of fidelity among districts. Model accuracy was validated by field statistics [43]. The IDW approach was used for spatial interpolation to evaluate CL transmission. IDW is an estimation method in which a linear combination of known values of sampling points is used to predict unknown values of non-sampling points with corresponding weighted values of inverse distance [44,45]. The calculation was carried out with the settings of neighboring type, 0.5 smoothing factor, and zero angles of spatial pattern.

2.5. Cluster and Outlier Analysis and Hotspot Analysis

The spatial statistics tools suite in ArcGIS was used. Spatial clusters and outliers were detected by calculating a local Moran’s I value, a z-score, a pseudo p-value, and a code representing the cluster type for each significant feature (Anselin Local Moran’s I). Statistically significant spatial clusters of high values (hot spots) and low values (cold spots) were identified via hot spot analysis (Getis-Ord Gi*). Fixed distance band was used to ensure each feature had at least one neighbor. The Euclidean distance technique was applied to calculate the distance of the neighboring zone around the indices. Spatial weights were not row standardized [43].

2.6. Bayesian Geostatistical Analysis

Climate data comprising 19 bioclimatic variables which represented varied temperature and precipitation conditions with 1 (30 s) to 340 (10 min) km2 resolution for the years 1970–2000 were obtained from WorldClim Global Climate Data (https://www.worldclim.org/data/worldclim21.html, accessed on 23 January 2023). The data were organized in Microsoft excel sheets and R-codes were created against each variable using R-software (windows version 4.2.2) for Bayesian geostatistical analysis [41,46].

2.7. Statistical Analysis

The data of CL patients were inputted to spreadsheets of Microsoft Excel (Windows version, 2016) and analyzed using SPSS (v. 26, IBM Corp. in Armonk, NY, USA). Differences between groups of each risk factor were compared using chi-square tests, a nonparametric method, after testing for a normal distribution. The observed difference was considered to be statistically significant if the p-value was less than 0.05. Associations between variables were examined using the chi-square test of independence with a significance level of 0.05. All variables were treated as categorical.

3. Results

A total of 3188 microscopically confirmed CL cases were reported in Khyber Pakhtunkhwa province in Pakistan in 2021. Sixty three percent of them were males (2003/3188, 62.8%). Male patients were more abundant in all age groups and in all studied districts (p < 0.001) (Figure 1). The age of patients ranged from 1 to 82 years, but the majority were younger than 20 years old (1800/3188, 56.5%) (p < 0.001). In fact, the number of patients was observed to be inversely correlated to age groups (Table 1). Regarding clinical presentation of CL, one to five lesions were observed at different body parts of the patients. Most patients had a single lesion (1540/3188, 48.3%) (p < 0.001), but there were 552 patients who had more than three lesions (552/3188, 17.3%) (Table 1). The face was most frequently attacked (1151/3188, 36.1%) (p < 0.001), followed by the hand (1068/3188, 33.5%), foot (633/3188, 19.9%), and multiple sites (336/3188, 10.5%). The dry form of lesion (2746/3188, 86.1%) was more common than the wet form (442/3188, 13.9%) (p < 0.001). The type of lesion was associated with season (p = 0.004) and district (p < 0.001).
Temporally, the incidence gradually increased from October (n = 212) and reached a peak in February (n = 643), which was followed by a sharp decrease in March (n = 313). Overall, most of infections occurred in winter (1443/3188, 45.26%), while fewer cases were found in autumn (446/3188, 13.98%) (p < 0.001) (Table 1). The seasonal trend was particularly apparent in the southern region (Table 2). However, in northern districts no patient was identified in Chitral in summer, although five male patients were reported in autumn. In Dir Lower, the fewest CL cases were recorded in spring (2/421, 0.5%), in contrast to Malakand, which had the highest number of infections in spring (107/250, 42.8%). The seasonal difference of cases in Swat was not significant. In the central region, there were almost as many patients in winter as in summer in Mardan. CL was least seen in Noshehra but most prevalent in Peshawar in spring. Charsadda did not report any patients in summer and autumn. In addition, the gender of patients was associated with season (p < 0.001). The sex ratios were 1.4, 1.8, 1.9, and 2.4 males/female for winter, spring, summer, and autumn, respectively. Spatially, CL was more prevalent in southern region (1501/3188, 47.08%) compared with northern (865/3188, 27.13%) and central regions (822/3188, 25.78%) (p < 0.001) (Table 1). The highest incidence of CL was observed in the city of Kaka Khal (n = 370) in the Bannu district, followed by the city of Domain (n = 263) and the city of Khal (n = 232) in the Dir Lower district. In total, 633 (633/3188, 19.9%) patients were from the Bannu district, 421 (421/3188, 13.2%) patients were from the Dir Lower district, and 411 (411/3188, 12.9%) patients were from the Mardan district.
CL was observed in 24 cities in 12 districts of Khyber Pakhtunkhwa, showing transmission of CL across the study area. The incidence of CL in each city was tagged to provide a visual concept of prevalence across the province (Figure 2A). The choropleth map showed the darkest color in Bannu in the southern region, highlighting the highest average case number in the district. The next dark color filled in Dir Lower in the northern region, Mardan in the central region, and Lakki Marwat in the southern region. The lightest color was used for Swat in the northern region, representing the fewest average cases (Figure 2B). The DEM map revealed that CL was found at altitudes of 125 to 7632 m, but incidence decreased with high altitudes. Cases were mostly distributed at low altitudes in the south (Figure 2C). An elevation range of 946 to 1977 m was identified to be associated with the highest occurrence of Leishmania infection.
IDW interpolation was applied to analyze the transmission patterns and therefore predict the future epidemic risk of CL. Threat focus areas were presumed to be located in the vicinity of the current CL-prevalent districts, including Bannu (IDW Spatial interpolation: 330.44–369.74), Dir Lower (IDW Spatial interpolation: 291.13–330.44), and Mardan (IDW Spatial interpolation: 251.83–291.13) (Figure 3A). The results of IDW showed high and low endemic areas across the province, with the highest value of focal statistics of 367.972 and the lowest value of 16.7361 (Figure 3B). Cluster and outlier analysis indicated Bannu (99% confidence interval, p < 0.05), Dir lower (95% confidence interval, p < 0.05), and Mardan (95% confidence interval, p < 0.05) as high cluster districts, as the other districts were not significant (Figure 3B). Hotspot analysis identified Bannu and Tank as significantly high clusters (hotspots) and Shangla as a low cluster district (coldspot) (Figure 3C).
Bayesian geostatistical analysis showed high temperature, low altitude, and annual precipitation were associated with the predicted highest number of cases in Kaka Khel (Bannu, p = 0.001), Khal (Dir Lower, p = 0.001), and Toro (Mardan, p = 0.001). Humidity and varied seasonal precipitation were associated with CL less affected areas such as Swat and Chitral.

4. Discussion

CL has been prevalent and is still expanding its frontier in Pakistan, including Khyber Pakhtunkhwa, a northern province bordering Afghanistan and Baluchistan [35,47]. One of the major outbreaks in that area occurred in the Timer Camp for Afghan refugees in 1997, indicating the spread of the disease due to the migration of the population [22]. However, transmission and geological expansion of CL were rather difficult to identify because of the diversity of disease incubation periods. This cross-sectional study was conducted to explore the spatial and temporal distribution of CL in Khyber Pakhtunkhwa. Clinico-epidemiological data of 3188 patients were collected from local health facilities throughout 2021. The seasonal occurrence of CL was analyzed statistically, and the geographical distribution, focal transmission, and hot spots were investigated using GIS tools.
Of 3188 patients reported in 2021, nearly 20% (633/3188, 19.9%) came from Bannu. This may be due to the high temperature and population density in Bannu compared with other districts in the province. Higher incidence of CL has been observed in areas with high population density and unhygienic environments [15]. Unsatisfactory housing conditions, such as cracked mud walls, dampness, and darkness as an indicator of poverty also increased the risk of leishmaniasis [48,49]. Dir Lower was found to be the second most affected district in the study. One of the reasons could be attributed to the inflow of Afghan refugees since 1970, but CL became even widespread among local inhabitants in combination with the abundance of sandflies in the district [50,51,52,53,54]. Moreover, the northern region was famous for various types of tourism and business activities. Migration of people between non-endemic and endemic areas would play a key role in the spread of CL [55]. Similar differences in prevalence between districts have been demonstrated in previous reports [5,17,32].
Sandfly species of Phlebotomus and Sergentomyia genera are abundant in Khyber Pakhtunkhwa [56,57,58,59,60]. Although there have been no conclusive reports to clarify the pathogen-vector association between Leishmania and sandflies in Pakistan, Phlebotomus papatasi and Phlebotomus sergenti, the two most common species in several districts, are considered to be vectors of CL, while Phlebotomus major and Phlebotomus hindustanicus are suspected vectors of VL [56]. Studies in Iran have identified Ph. papatasi as a vector of Leishmania major [61,62], and Ph. sergenti has been shown to transmit Leishmania tropica in Morocco [63]. Phlebotomus caucasicus and Phlebotomus salehi were also indicated as potential vectors of CL in Middle East countries [64,65,66]. Other sandflies found in Khyber Pakhtunkhwa included Phlebotomus salengensis, Phlebotomus andrejevi, Phlebotomus kazeroni, Phlebotomus bergeroti, Phlebotomus ansari, Phlebotomus alexandri, Sergentomyia babu, Sergentomyia dentata, Sergentomyia baghdadis, Sergentomyia bailyi, Sergentomyia hospittii, Sergentomyia montana, Sergentomyia grekovi, Sergentomyia hodgsoni, Sergentomyia turkistnica, Sergentomyia tiberidis, Sergentomyia theodori, Sergentomyia sumbarica, Sergentomyia dreyfussi turkestanica, and Sergentomyia fallax afghanica [57,58,59]. On the other hand, records about reservoir hosts of CL are scarce in Pakistan. Great gerbils (Rhombomys opimus) were the primary natural hosts of Leishmania spp. in many countries in Central Asia [66]. The role great gerbils and other small mammals, e.g., Meriones erythrourus, played in the transmission cycle of Leishmania in Khyber Pakhtunkhwa remains to be investigated.
In our study, males appeared to be more affected than females. This may be due to the fact that females cover themselves and remain indoors most of the time in Pakistan. Males were more exposed to the bites of sandflies as they went for labor with their faces, shoulders, and hands uncovered. The association between season and gender further supported the hypothesis. The sex ratios of CL patients were 1.4, 1.8, 1.9, and 2.4 males/females for winter, spring, summer, and autumn, respectively (853/590, 426/240, 411/222, and 313/133). The low temperature in winter made both females and males stay indoors at night, when sandflies were seeking bloodmeals actively, while the warmer weather in summer would allow male workers to go outdoors after sunset. CL infection was observed across all ages, but children and young people under the age of 20 years were most vulnerable. Playing outside in the shade of trees or near humid ground surfaces where sandflies bred placed children at high risk of vector bites. Their immature immune systems also led to increasing susceptibility to infection [67,68,69].
Our data indicated that CL was more prevalent in winter, especially in the southern region, although other studies identified a high incidence in spring [32,70,71]. The occurrence of CL was in fact associated with the fluctuation of sandflies in the rainy season [72]. In the northern region, frontal cloud bands result in precipitation in winter. Heavy thunderstorms sometimes happen in Chitral in spring. Dir Lower is one of the wettest districts in Pakistan with an annual rainfall of 1473.2 mm, of which 400 mm are brought by the monsoon in summer. Almost twice that amount takes place during December to April, showing a bimodal rainfall regime. CL infection was mostly reported in winter (n = 221) and summer (n = 131), echoing the precipitation pattern.
Clinically, most of the patients in the study had a single facial lesion, consistent with previous findings [73,74]. The face was the body part most exposed to sandfly bites. The situation was worsened if people did not carefully take preventive measures, such as wearing repellents, using bed nets, and avoiding outdoor activities during the time sandflies were active [48]. Lesions over five centimeters were noticed on seven patients. The size of lesion was associated with limited access to medical care [48,75,76]. As a result, patients with low socio-economic status suffering from CL might have to live with facial disfigurements in developing countries. On the other hand, we also found the type of lesion was significantly associated with temporal and spatial factors (season and districts, p = 0.004 and p < 0.001, respectively). The wet type of lesion was chiefly caused by L. major and had a shorter incubation period (two weeks). This type of CL, also referred as zoonotic CL, was identified to be prevalent at low elevations in arid and semi-arid areas in Pakistan [3]. The dry form, caused by L. tropica and also called anthroponotic CL, was transmitted at high elevations (500–2837 m) [3]. The ratio of patients with dry to wet form CL was lowest in autumn (360/86, 4.2), while it was higher in winter and spring (1260/183 and 582/84, 6.9 and 6.9, respectively). One limitation of the study was the lack of species identification for the etiological agents causing the observed lesions due to limited resources in the local health facilities and the study design. However, previous reports have shown that L. tropica and L. major were responsible for most CL cases in Khyber Pakhtunkhwa [3,23,48]. Other pathogens, such as Leishmania infantum, have been detected in patients, but the ratios were quite low [47]. The Leishmania species causing CL in Khyber Pakhtunkhwa and their dynamics will be further investigated in the future.
The spatial pattern of CL patients was presented with a choropleth map. The continuous color represented the abundance of patients varied in districts. There were more patients reported in Bannu, in contrast to fewer patients reported in Swat and Charsadda. Visualization of case distribution would allow public health sectors to easily assess regional severity of the infectious disease [77]. Meanwhile, the vertical occurrence of CL was illustrated by a DEM map. CL patients were distributed between the elevations of 125 to 7632 m, although most of the patients lived at 946 to 1977 m. The vertical dissemination of CL would be restricted by the activities of sandflies, which prefer warm weather at low altitudes.
The transmission of CL was predicted using the IDW method based on the hypothesis that the chance of being infected decreased with growing distance to an existing case, as sandflies are weak fliers. The CL occurrence of neighborhood locations close to the cases was estimated by spatial interpolation, as both locations shared similar characteristics, including land use, lifestyle of the inhabitants, and breeding sources of sandflies. The interpolation could provide public health sectors with a reference for disease prevention [77]. Risk of CL was further assessed using cluster and outlier analysis as well as hotspot analysis. Apparent clusters were identified in areas around the current CL-prevalent districts including Bannu and Dir Lower. Hotspot analysis also identified Bannu as an area of high clustering (hotspots). Previous studies have described household clusters and climate as the main factors affecting transmission of CL [78,79]. Bayesian geostatistical analysis was then applied to evaluate the effects of environmental and climatic variables on the spatial distribution of the disease. The results indicated a high incidence in the southern region was associated with low altitudes, high temperature, and increasing annual precipitation. Warm weather and precipitation have been shown to facilitate the transmission of CL, as the development of both sandflies and Leishmania accelerate at high temperatures, and precipitation results in more suitable breeding sites for sandflies [41].
Although the WHO has taken action in an attempt to control CL in endemic areas, CL remains a major public health problem in Pakistan [80,81]. GIS has been extensively applied to describe spatial patterns and shows great promise in research and surveillance of diseases. This study utilized GIS tools to illustrate the horizontal and vertical distribution of CL in Khyber Pakhtunkhwa, Pakistan. Transmission, risk, and hotspots of the disease were evaluated and presented visually. Clinical presentations as well as risk factors including gender, age, season, and region were analyzed statistically. Additionally, Bayesian inference was used in geostatistics to explore the relationship between CL occurrence and climatic variables. The results could help public health sectors to establish prevention and control strategies for upcoming CL epidemics. For instance, health education especially targeting risk groups could be provided. Combined with population data, intervention could be prioritized in areas with the highest disease risk, e.g., the Bannu district, and resources could be distributed to those in need. Efforts for disease prevention and vector control could be made before the epidemic season of each district specifically according to the results of temporal analysis. Public health workers should help to raise the awareness of clinicians and residents in hotspots, thus allowing infected patients to receive treatment earlier. When climate factors favor the growth of vectors, special attention has to be paid to surveillance, and sufficient medical supplies must be prepared for potential outbreaks. Moreover, from a worldwide perspective, this study could be beneficial for global CL control.

5. Conclusions

The study presented an overview of CL in Khyber Pakhtunkhwa in Pakistan. The majority of patients were males or individuals under the age of 20 years. Infection mostly occurred in the southern region in winter. GIS tools were used to visualize the spatial patterns of the disease, and high clusters were found in Bannu, Dir Lower, and Mardan. Identification of risk areas would be helpful for public health sectors to predict future epidemics and implement corresponding interventions.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/tropicalmed8020128/s1, Figure S1: Flowchart of the study.

Author Contributions

Conceptualization: W.U., A.K., K.-H.T., and S.N.; methodology: W.U. and N.N.; software: W.U. and A.K.; validation: K.-H.T., T.-Y.Y., Y.-F.T., and A.K.; formal analysis: T.-Y.Y. and R.I.R.-V.; investigation: W.U.; resources: K.-H.T. and T.-Y.Y.; data curation: T.-Y.Y., S.N., Y.-F.T., and N.N.; writing—original draft preparation: W.U. and T.-Y.Y.; writing—review and editing: K.-H.T. and T.-Y.Y.; visualization: A.K. and R.I.R.-V.; supervision: S.N.; project administration: S.N and N.N.; funding acquisition: K.-H.T. and T.-Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partially supported by the Vector-borne Infectious Diseases Laboratory, College of Public Health, National Taiwan University, Taipei, Taiwan (NTUCPH-2022-SEA-01).

Institutional Review Board Statement

The study was approved under the approval no. AWKUM-20053310 by the Ethical Committee of Chemical and Life Section, Department of Zoology, Abdul Wali Khan University Mardan, Pakistan.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the local district headquarter hospitals for their administrative and technical support in sample and data collection.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Mandell, G.L.; Bennett, J.E.; Dolin, R. Principles and Practice of Infectious Diseases, 6th ed.; Elsevier Churchill Livingstone: Philadelphia, PA, USA, 2005; pp. 1864–1890. [Google Scholar]
  2. Malani, P.N. Mandell, Douglas, and Bennett’s principles and practice of infectious diseases. JAMA 2010, 304, 2067–2071. [Google Scholar] [CrossRef]
  3. Khan, K.; Khan, N.H.; Wahid, S. Systematic review of leishmaniasis in Pakistan: Evaluating spatial distribution and risk factors. J. Parasitol. 2021, 107, 630–638. [Google Scholar] [CrossRef]
  4. Arevalo, J.; Ramirez, L.; Adaui, V.; Zimic, M.; Tulliano, G.; Miranda-Verástegui, C.; Lazo, M.; Loayza-Muro, R.; De Doncker, S.; Maurer, A.; et al. Influence of Leishmania (Viannia) species on the response to antimonial treatment in patients with American tegumentary leishmaniasis. J. Infect. Dis. 2007, 195, 1846–1851. [Google Scholar] [CrossRef] [Green Version]
  5. Demirel, R.; Erdogan, S. Determination of high risk regions of cutaneous leishmaniasis in Turkey using spatial analysis. Turkiye Parazitol. Derg. 2009, 33, 8–14. [Google Scholar]
  6. Pourmohammadi, B.; Motazedian, M.; Hatam, G.; Kalantari, M.; Habibi, P.; Sarkari, B. Comparison of three methods for diagnosis of cutaneous leishmaniasis. Iran. J. Parasitol. 2010, 5, 1–8. [Google Scholar]
  7. Davies, C.R.; Kaye, P.; Croft, S.L.; Sundar, S. Leishmaniasis: New approaches to disease control. BMJ 2003, 326, 377–382. [Google Scholar] [CrossRef] [Green Version]
  8. Reithinger, R.; Dujardin, J.C.; Louzir, H.; Pirmez, C.; Alexander, B.; Brooker, S. Cutaneous leishmaniasis. Lancet Infect. Dis. 2007, 7, 581–596. [Google Scholar] [CrossRef] [Green Version]
  9. Kumar, R.; Bumb, R.A.; Ansari, N.A.; Mehta, R.D.; Salotra, P. Cutaneous leishmaniasis caused by Leishmania tropica in Bikaner, India: Parasite identification and characterization using molecular and immunologic tools. Am. J. Trop. Med. Hyg. 2007, 76, 896–901. [Google Scholar] [CrossRef]
  10. Lainson, R.; Ryan, L.; Shaw, J.J. Infective stages of Leishmania in the sandfly vector and some observations on the mechanism of transmission. Mem. Inst. Oswaldo Cruz 1987, 82, 421–424. [Google Scholar] [CrossRef]
  11. Shah, A.A.; Funk, W.C.; Ghalambor, C.K. Thermal acclimation ability varies in temperate and tropical aquatic insects from different elevations. Integr. Comp. Biol. 2017, 57, 977–987. [Google Scholar] [CrossRef] [Green Version]
  12. Peter, H.; Hand, I. Patterns of patient-spouse interaction in agoraphobics: Assessment by Camberwell Family Interview (CFI) and impact on outcome of self-exposure treatment. In Panic and Phobias, 2nd ed.; Springer: Cham, Switzerland, 1988; pp. 240–251. [Google Scholar]
  13. ul Bari, A. Chronology of cutaneous leishmaniasis: An overview of the history of the disease. J. Pak. Assoc. Dermatol. 2006, 16, 24–27. [Google Scholar]
  14. Grevelink, S.A.; Lerner, E. Leishmaniasis. J. Am. Acad. Dermatol. 1996, 34, 80121–80126. [Google Scholar] [CrossRef] [PubMed]
  15. Ashford, R.W. The leishmaniases as emerging and reemerging zoonoses. Int. J. Parasitol. 2000, 30, 1269–1281. [Google Scholar] [CrossRef] [PubMed]
  16. Daszak, P.; Cunningham, A.A.; Hyatt, A.D. Anthropogenic environmental change and the emergence of infectious diseases in wildlife. Acta Trop. 2001, 78, 103–116. [Google Scholar] [CrossRef] [PubMed]
  17. Patz, J.A.; Graczyk, T.K.; Geller, N.; Vittor, A.Y. Effects of environmental change on emerging parasitic diseases. Int. J. Parasitol. 2000, 30, 1395–1405. [Google Scholar] [CrossRef] [Green Version]
  18. Jeddi-Tehrani, M.; Torabi, R.; Zarnani, A.H.; Mohammadzadeh, A.; Arefi, S.; Zeraati, H.; Akhondi, M.M.; Chamani-Tabriz, L.; Idali, F.; Emami, S.; et al. Analysis of plasminogen activator inhibitor-1, integrin beta3, beta fibrinogen, and methylenetetrahydrofolate reductase polymorphisms in Iranian women with recurrent pregnancy loss. Am. J. Reprod. Immunol. 2011, 66, 149–156. [Google Scholar] [CrossRef]
  19. Alvar, J.; Vélez, I.D.; Bern, C.; Herrero, M.; Desjeux, P.; Cano, J.; Jannin, J.; den Boer, M.; WHO Leishmaniasis Control Team. Leishmaniasis worldwide and global estimates of its incidence. PLoS ONE 2012, 7, e35671. [Google Scholar] [CrossRef]
  20. Kassi, M.; Kassi, M.; Afghan, A.K.; Rehman, R.; Kasi, P.M. Marring leishmaniasis: The stigmatization and the impact of cutaneous leishmaniasis in Pakistan and Afghanistan. PLoS Negl. Trop. Dis. 2008, 2, e259. [Google Scholar] [CrossRef] [Green Version]
  21. Saifullah Khan, S.; Shahab, M.H.; Naz, T. Visceral leishmaniasis in a resident of Swat, Khyber Pakhtunkhwa, Pakistan, presenting to Civil Hospital Karachi: A case report. Cureus 2019, 11, e6059. [Google Scholar] [CrossRef]
  22. Rowland, M.; Munir, A.; Durrani, N.; Noyes, H.; Reyburn, H. An outbreak of cutaneous leishmaniasis in an Afghan refugee settlement in north-west Pakistan. Trans. R. Soc. Trop. Med. Hyg. 1999, 93, 133–136. [Google Scholar] [CrossRef]
  23. Hussain, M.; Munir, S.; Khan, T.A.; Khan, A.; Ayaz, S.; Jamal, M.A.; Ahmed, I.; Aziz, S.; Watany, N.; Kasbari, M. Epidemiology of cutaneous leishmaniasis outbreak, Waziristan, Pakistan. Emerg. Infect. Dis. 2018, 24, 159–161. [Google Scholar] [CrossRef] [Green Version]
  24. Hussain, M.; Munir, S.; Jamal, M.A.; Ayaz, S.; Akhoundi, M.; Mohamed, K. Epidemic outbreak of anthroponotic cutaneous leishmaniasis in Kohat District, Khyber Pakhtunkhwa, Pakistan. Acta Trop. 2017, 172, 147–155. [Google Scholar] [CrossRef]
  25. Musa, G.J.; Chiang, P.H.; Sylk, T.; Bavley, R.; Keating, W.; Lakew, B.; Tsou, H.C.; Hoven, C.W. Use of GIS mapping as a public health tool-from cholera to cancer. Health Serv. Insights 2013, 6, 111–116. [Google Scholar] [CrossRef]
  26. Gwitira, I.; Mukonoweshuro, M.; Mapako, G.; Shekede, M.D.; Chirenda, J.; Mberikunashe, J. Spatial and spatio-temporal analysis of malaria cases in Zimbabwe. Infect. Dis. Poverty. 2020, 9, 146. [Google Scholar] [CrossRef]
  27. Duncombe, J.; Clements, A.; Hu, W.; Weinstein, P.; Ritchie, S.; Espino, F.E. Geographical information systems for dengue surveillance. Am. J. Trop. Med. Hyg. 2012, 86, 753–755. [Google Scholar] [CrossRef] [Green Version]
  28. Bavia, M.E.; Carneiro, D.D.; Gurgel Hda, C.; Madureira Filho, C.; Barbosa, M.G. Remote sensing and geographic information systems and risk of American visceral leishmaniasis in Bahia, Brazil. Parassitologia 2005, 47, 165–169. [Google Scholar]
  29. Malone, J.B.; Bergquist, R.; Martins, M.; Luvall, J.C. Use of geospatial surveillance and response systems for vector-borne diseases in the elimination phase. Trop. Med. Infect. Dis. 2019, 4, 15. [Google Scholar] [CrossRef] [Green Version]
  30. Kuleshov, Y.; Wei, Y.; Inape, K.; Liu, G.J. Spatio-temporal distribution of vector borne diseases in Australia and Papua New Guinea vis-à-vis climatic factors. J. Vector Borne Dis. 2022, 59, 115–126. [Google Scholar] [CrossRef]
  31. Karunaweera, N.D.; Ginige, S.; Senanayake, S.; Silva, H.; Manamperi, N.; Samaranayake, N.; Siriwardana, Y.; Gamage, D.; Senerath, U.; Zhou, G. Spatial epidemiologic trends and hotspots of leishmaniasis, Sri Lanka, 2001–2018. Emerg. Infect. Dis. 2020, 26, 1–10. [Google Scholar] [CrossRef] [Green Version]
  32. Abedi-Astaneh, F.; Hajjaran, H.; Yaghoobi-Ershadi, M.R.; Hanafi-Bojd, A.A.; Mohebali, M.; Shirzadi, M.R.; Rassi, Y.; Akhavan, A.A.; Mahmoudi, B. Risk mapping and situational analysis of cutaneous leishmaniasis in an endemic area of Central Iran: A GIS-based survey. PLoS ONE 2016, 11, e0161317. [Google Scholar] [CrossRef]
  33. Eisen, L.; Eisen, R.J. Using geographic information systems and decision support systems for the prediction, prevention, and control of vector-borne diseases. Annu. Rev. Entomol. 2011, 56, 41–61. [Google Scholar] [CrossRef] [PubMed]
  34. Ostad, M.; Shirian, S.; Pishro, F.; Abbasi, T.; Ai, A.; Azimi, F. Control of cutaneous leishmaniasis using geographic information systems from 2010 to 2014 in Khuzestan Province, Iran. PLoS ONE 2016, 11, e0159546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zaidi, F.; Fatima, S.H.; Jan, T.; Fatima, M.; Ali, A.; Khisroon, M.; Adnan, M.; Rasheed, S.B. Environmental risk modelling and potential sand fly vectors of cutaneous leishmaniasis in Chitral district: A leishmanial focal point of mount Tirich Mir, Pakistan. Trop. Med. Int. Health 2017, 22, 1130–1140. [Google Scholar] [CrossRef]
  36. Khan, K.; Khan, N.H.; Anwar, F.; Ullah, I.; Badshah, N.; Irfan, I.; Iqbal, K.; Shah, I.; Aziz, S.T.; Khan, M.S.N.; et al. Characterization of sand fly breeding sites in district Malakand, Khyber Pakhtunkhwa, Pakistan, and evaluation of risk factors for cutaneous leishmaniasis in the region. Zoonoses Public Health 2022, 69, 33–45. [Google Scholar] [CrossRef]
  37. Wahid, S.; Khan, K.; Khan, N.H. Sand fly (Diptera: Psychodidae) species diversity, habitat preferences, and ecological aspects of distribution in Bajaur District, Khyber Pakhtunkhwa, Pakistan. J. Med. Entomol. 2020, 57, 1432–1439. [Google Scholar] [CrossRef] [PubMed]
  38. Johnson, L.R.; Ben-Horin, T.; Lafferty, K.D.; McNally, A.; Mordecai, E.; Paaijmans, K.P.; Pawar, S.; Ryan, S.J. Understanding uncertainty in temperature effects on vector-borne disease: A Bayesian approach. Ecology 2015, 96, 203–213. [Google Scholar] [CrossRef] [Green Version]
  39. Adigun, A.B.; Gajere, E.N.; Oresanya, O.; Vounatsou, P. Malaria risk in Nigeria: Bayesian geostatistical modelling of 2010 malaria indicator survey data. Malar. J. 2015, 14, 156. [Google Scholar] [CrossRef] [Green Version]
  40. Schur, N.; Hürlimann, E.; Garba, A.; Traoré, M.S.; Ndir, O.; Ratard, R.C.; Tchuem Tchuenté, L.A.; Kristensen, T.K.; Utzinger, J.; Vounatsou, P. Geostatistical model-based estimates of Schistosomiasis prevalence among individuals aged ≤ 20 years in West Africa. PLoS Negl. Trop. Dis. 2011, 5, e1194. [Google Scholar] [CrossRef] [Green Version]
  41. Karagiannis-Voules, D.A.; Scholte, R.G.; Guimarães, L.H.; Utzinger, J.; Vounatsou, P. Bayesian geostatistical modeling of leishmaniasis incidence in Brazil. PLoS Negl. Trop. Dis. 2013, 7, e2213. [Google Scholar] [CrossRef] [Green Version]
  42. Jewell, C.P.; Brown, R.G. Bayesian data assimilation provides rapid decision support for vector-borne diseases. J. R. Soc. Interface 2015, 12, 20150367. [Google Scholar] [CrossRef]
  43. Pakzad, R.; Dabbagh-Moghaddam, A.; Mohebali, M.; Safiri, S.; Barati, M. Spatio-temporal analysis of cutaneous leishmaniasis using geographic information system among Iranian Army Units and its comparison with the general population of Iran during 2005–2014. J. Parasit. Dis. 2017, 41, 1114–1122. [Google Scholar] [CrossRef] [PubMed]
  44. Messina, J.P.; Taylor, S.M.; Meshnick, S.R.; Linke, A.M.; Tshefu, A.K.; Atua, B.; Mwandagalirwa, K.; Emch, M. Population, behavioural and environmental drivers of malaria prevalence in the Democratic Republic of Congo. Malar. J. 2011, 10, 161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Moeletsi, M.E.; Shabalala, Z.P.; de Nysschen, G.; Walker, S. Evaluation of an inverse distance weighting method for patching daily and dekadal rainfall over the free state province, South Africa. Water SA 2016, 42, 466–474. [Google Scholar] [CrossRef] [Green Version]
  46. Becker, R.; Wilks, A.R. Maps: Draw Geographical Maps. Version 3.3.0 2017. Available online: https://cran.r-project.org/web/packages/maps/index.html (accessed on 21 November 2022).
  47. Khan, N.H.; Bari, A.U.; Hashim, R.; Khan, I.; Muneer, A.; Shah, A.; Wahid, S.; Yardley, V.; O’Neil, B.; Sutherland, C.J. Cutaneous leishmaniasis in Khyber Pakhtunkhwa province of Pakistan: Clinical diversity and species-level diagnosis. Am. J. Trop. Med. Hyg. 2016, 95, 1106–1114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Zeb, I.; Ali, A.; Nawab, J.; Khan, M.Q.; Kamil, A.; Tsai, K.H. Cutaneous leishmaniasis in male schoolchildren in the upper and lower Dir districts of Khyber Pakhtunkhwa, and a review of previous record in Pakistan. Acta Trop. 2020, 209, 105578. [Google Scholar] [CrossRef]
  49. Calderon-Anyosa, R.; Galvez-Petzoldt, C.; Garcia, P.J.; Carcamo, C.P. Housing characteristics and leishmaniasis: A systematic review. Am. J. Trop. Med. Hyg. 2018, 99, 1547–1554. [Google Scholar] [CrossRef] [Green Version]
  50. Kolaczinski, J.; Brooker, S.; Reyburn, H.; Rowland, M. Epidemiology of anthroponotic cutaneous leishmaniasis in Afghan refugee camps in northwest Pakistan. Trans. R. Soc. Trop. Med. Hyg. 2004, 98, 373–378. [Google Scholar] [CrossRef]
  51. Ali, A.; Rehman, T.U.; Qureshi, N.A.; Rahman, H.U. New endemic focus of cutaneous leishmaniasis in Pakistan and future epidemics threats. Asian Pac. J. Trop. Dis. 2016, 6, 155–159. [Google Scholar] [CrossRef]
  52. Khan, K.; Wahid, S.; Shah, A.; Ali, N. First report of sand fly (Diptera: Psycodidae: Phlebotomine) immature males and various reproductive stages of females in Upper and Lower Dir Districts, Khyber Pakhtunkhwa Pakistan. Pak. J. Zool. 2016, 48, 1581–1584. [Google Scholar]
  53. Khan, K.; Wahid, S.; Khan, N.H.; Ali, N. Potential resting and breeding sites of sand flies (Diptera: Psychodidae) and their habitat characteristics in leishmaniasis foci of Dir districts, Khyber Pakhtunkhwa, Pakistan. J. Med. Entomol. 2017, 54, 1390–1396. [Google Scholar] [CrossRef]
  54. Zeb, I.; Qureshi, N.A.; Shaheen, N.; Zafar, M.I.; Ali, A.; Hamid, A.; Shah, S.A.A.; Ashraf, A. Spatiotemporal patterns of cutaneous leishmaniasis in the district Upper and Lower Dir, Khyber Pakhtunkhwa, Pakistan: A GIS-based spatial approaches. Acta Trop. 2021, 217, 105861. [Google Scholar] [CrossRef] [PubMed]
  55. Torres-Guerrero, E.; Quintanilla-Cedillo, M.R.; Ruiz-Esmenjaud, J.; Arenas, R. Leishmaniasis: A review. F1000Research 2017, 6, 750. [Google Scholar] [CrossRef] [PubMed]
  56. Khan, K.; Wahid, S.; Khan, N.H. Habitat characterization of sand fly vectors of leishmaniasis in Khyber Pakhtunkhwa, Pakistan. Acta Trop. 2019, 199, 105147. [Google Scholar] [CrossRef]
  57. Ali, N.; Khan, K.; Wahid, S.; Khan, N.H.; Shah, S.U. Species composition and activity patterns of sand flies (Psycodidae: Phlebotomine) in four tehsils of Dir Districts, Khyber Pakhtunkhwa, Pakistan. Acta Trop. 2016, 156, 144–151. [Google Scholar] [CrossRef] [PubMed]
  58. Ali, N.; Ullah, A.; Wahid, S.; Khisroon, M.; Rasheed, S.B. Distribution, species composition and relative abundances of sandflies in North Waziristan Agency, Pakistan. Med. Vet. Entomol. 2016, 30, 89–94. [Google Scholar] [CrossRef]
  59. Khan, K.; Sajjad, M.; Wahid, S.; Gul, M.; Khan, L.; Ullah, H.; Rahman, Y.; Khan, D.; Khan, K.; Khan, M.Y.; et al. Bionomics of the unexplored sand flies fauna of District Mohmand, Khyber Pakhtunkhwa, Pakistan: Assessing risk factors associated with cutaneous leishmaniasis. Trans. R. Soc. Trop. Med. Hyg. 2022, 116, 832–844. [Google Scholar] [CrossRef]
  60. Durrani, A.Z.; Durrani, H.Z.; Kamal, N. Prevalence of Leishmania in sand fly in Pakistan. Pakistan J. Zool. 2012, 44, 61–65. [Google Scholar]
  61. Parvizi, P.; Mauricio, I.; Aransay, A.M.; Miles, M.A.; Ready, P.D. First detection of Leishmania major in peridomestic Phlebotomus papatasi from Isfahan province, Iran: Comparison of nested PCR of nuclear ITS ribosomal DNA and semi-nested PCR of minicircle kinetoplast DNA. Acta Trop. 2005, 93, 75–83. [Google Scholar] [CrossRef]
  62. Yaghoobi-Ershadi, M.R.; Javadian, E.; Tahvildare-Bidruni, G.H. Leishmania major MON-26 isolated from naturally infected Phlebotomus papatasi (Diptera: Psychodidae) in Isfahan Province, Iran. Acta Trop. 1995, 59, 279–282. [Google Scholar] [CrossRef]
  63. Es-Sette, N.; Ajaoud, M.; Bichaud, L.; Hamdi, S.; Mellouki, F.; Charrel, R.N.; Lemrani, M. Phlebotomus sergenti a common vector of Leishmania tropica and Toscana virus in Morocco. J. Vector. Borne Dis. 2014, 51, 86–90. [Google Scholar]
  64. Azizi, K.; Fakoorziba, M.R.; Jalali, M.; Moemenbellah-Fard, M.D. First molecular detection of Leishmania major within naturally infected Phlebotomus salehi from a zoonotic cutaneous leishmaniasis focus in southern Iran. Trop. Biomed. 2012, 29, 1–8. [Google Scholar] [PubMed]
  65. Zahraei-Ramazani, A.; Saghafipour, A.; Mehdi Sedaghat, M.; Absavaran, A.; Azarm, A. Molecular identification of Phlebotomus caucasicus and Phlebotomus mongolensis (Diptera: Psychodidae) in a hyperendemic area of zoonotic cutaneous leishmaniasis in Iran. J. Med. Entomol. 2017, 54, 1525–1530. [Google Scholar] [CrossRef] [PubMed]
  66. Yurchenko, V.; Chistyakov, D.S.; Akhmadishina, L.V.; Lukashev, A.N.; Sádlová, J.; Strelkova, M.V. Revisiting epidemiology of leishmaniasis in central Asia: Lessons learnt. Parasitology 2023, 150, 129–136. [Google Scholar] [CrossRef] [PubMed]
  67. Yemisen, M.; Ulas, Y.; Celik, H.; Aksoy, N. Epidemiological and clinical characteristics of 7172 patients with cutaneous leishmaniasis in Sanliurfa, between 2001 and 2008. Int. J. Dermatol. 2012, 51, 300–304. [Google Scholar] [CrossRef] [PubMed]
  68. Bennis, I.; de Brouwere, V.; Ameur, B.; El Idrissi Laamrani, A.; Chichaoui, S.; Hamid, S.; Boelaert, M. Control of cutaneous leishmaniasis caused by Leishmania major in south-eastern Morocco. Trop. Med. Int. Health 2015, 20, 1297–1305. [Google Scholar] [CrossRef] [PubMed]
  69. Haouas, N.; Amer, O.; Ishankyty, A.; Alazmi, A.; Ishankyty, I. Profile and geographical distribution of reported cutaneous leishmaniasis cases in Northwestern Saudi Arabia, from 2010 to 2013. Asian Pac. J. Trop. Med. 2015, 8, 287–291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Hayat, U.; Ayaz, S.; Akram, M.; Khattak, A.M. Cutaneous leishmaniasis: Its prevalence and role of PCR in its detection. J. Islamabad Med. Dent. Coll. 2015, 4, 15–18. [Google Scholar]
  71. Khosravani, M.; Moemenbellah-Fard, M.D.; Sharafi, M.; Rafat-Panah, A. Epidemiologic profile of oriental sore caused by Leishmania parasites in a new endemic focus of cutaneous leishmaniasis, southern Iran. J. Parasit. Dis. 2016, 40, 1077–1081. [Google Scholar] [CrossRef] [Green Version]
  72. Mann, R.S.; Kaufman, P.E. The seasonal abundance of phlebotomine sand flies, Lutzomyia species in Florida. J. Am. Mosq. Control Assoc. 2010, 26, 10–17. [Google Scholar] [CrossRef]
  73. Khattak, A.S.; Khan, S.; Shams, S.; Saqalain, M.; Ahmad, J.; Bibi, A.; Ahmad, M.; Noreen, S.; Hussain, M. Cutaneous leishmaniasis in Karak, Pakistan: Report of an outbreak and comparison of diagnostic techniques. Afr. J. Biotechnol. 2011, 10, 9908–9910. [Google Scholar] [CrossRef]
  74. Fakhar, M.; Karamian, M.; Ghatee, M.A.; Taylor, W.R.; Pazoki Ghohe, H.; Rasooli, S.A. Distribution pattern of anthroponotic cutaneous leishmaniasis caused by Leishmania tropica in Western Afghanistan during 2013–2014. Acta Trop. 2017, 176, 22–28. [Google Scholar] [CrossRef]
  75. Machado, P.; Araújo, C.; Da Silva, A.T.; Almeida, R.P.; D’Oliveira, A., Jr.; Bittencourt, A.; Carvalho, E.M. Failure of early treatment of cutaneous leishmaniasis in preventing the development of an ulcer. Clin. Infect. Dis. 2002, 34, e69–e73. [Google Scholar] [CrossRef] [PubMed]
  76. Sandanayaka, R.; Kahawita, I.; Gamage, A.; Siribaddana, S.; Agampodi, S. Emergence of cutaneous leishmaniasis in Polonnaruwa, Sri Lanka 2008–2011. Trop. Med. Int. Health. 2014, 19, 140–145. [Google Scholar] [CrossRef] [Green Version]
  77. Osei, F.B.; Stein, A. Spatial variation and hot-spots of district level diarrhea incidences in Ghana: 2010–2014. BMC Public Health 2017, 17, 617. [Google Scholar] [CrossRef] [PubMed]
  78. Bhunia, G.S.; Kesari, S.; Chatterjee, N.; Kumar, V.; Das, P. Spatial and temporal variation and hotspot detection of kala-azar disease in Vaishali district (Bihar), India. BMC Infect. Dis. 2013, 13, 64. [Google Scholar] [CrossRef] [Green Version]
  79. Chaves, L.F.; Pascual, M. Climate cycles and forecasts of cutaneous leishmaniasis, a nonstationary vector-borne disease. PLoS Med. 2006, 3, e295. [Google Scholar] [CrossRef]
  80. Ministry of Health, Pakistan; World Health Organization; Health Net International. Guidelines for the Treatment and Prevention of Cutaneous Leishmaniasis in Pakistan; World Health Organization: Geneva, Switzerland, 2002. [Google Scholar]
  81. Akram, A.; Khan, H.A.; Qadir, A.; Sabir, A.M. A cross-sectional survey of knowledge, attitude and practices related to cutaneous leishmaniasis and sand flies in Punjab, Pakistan. PLoS ONE 2015, 10, e0130929. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A proportional symbol map of the study area. Location of Khyber Pakhtunkhwa province is shown in the upper right corner. Gender-wise distribution of total cases reported during 2021 is scaled in colored squares.
Figure 1. A proportional symbol map of the study area. Location of Khyber Pakhtunkhwa province is shown in the upper right corner. Gender-wise distribution of total cases reported during 2021 is scaled in colored squares.
Tropicalmed 08 00128 g001
Figure 2. (A) Cutaneous leishmaniasis (CL) cases in each city, showing the highest and lowest incidences across the province; (B) a choropleth map showing district-based CL cases; (C) a digital elevation model (DEM) map illustrating the occurrence of CL at different altitudes.
Figure 2. (A) Cutaneous leishmaniasis (CL) cases in each city, showing the highest and lowest incidences across the province; (B) a choropleth map showing district-based CL cases; (C) a digital elevation model (DEM) map illustrating the occurrence of CL at different altitudes.
Tropicalmed 08 00128 g002
Figure 3. (A) Cutaneous leishmaniasis (CL) transmission analyzed by inverse distance weighting (IDW) interpolation; (B) focal statistics with cluster and outlier analysis; (C) hotspot analysis.
Figure 3. (A) Cutaneous leishmaniasis (CL) transmission analyzed by inverse distance weighting (IDW) interpolation; (B) focal statistics with cluster and outlier analysis; (C) hotspot analysis.
Tropicalmed 08 00128 g003
Table 1. Characteristics of cutaneous leishmaniasis patients in Khyber Pakhtunkhwa, Pakistan, 2021.
Table 1. Characteristics of cutaneous leishmaniasis patients in Khyber Pakhtunkhwa, Pakistan, 2021.
VariableDescriptionCase no.Percentage (%)p Value
GenderMale200362.8<0.001
Female118537.2
Age group1–20180056.5<0.001
21–4087227.4
41–6040912.8
>601073.4
Season *Winter144345.3<0.001
Spring66620.9
Summer63319.9
Autumn44614.0
Lesion siteFace115136.1<0.001
Hand106833.5
Foot63319.9
Multiple33610.5
Lesion typeDry274686.1<0.001
Wet44213.9
Lesion numberOne154048.3<0.001
Two62619.6
Three47014.7
>Three55217.3
Northern regionChitral1504.7<0.001
Dir Lower42113.2
Malakand2507.8
Swat441.4
Central regionMardan41112.9<0.001
Charsadda642.0
Noshehra2287.2
Peshawar1193.7
Southern regionBannu63319.9<0.001
DI Khan2156.7
Lakki Marwat35411.1
Tank2999.4
* Winter: December ~ February; spring: March ~ May; summer: June ~ September; autumn: October ~ November.
Table 2. Seasonal and gender-wise prevalence of CL in Khyber Pakhtunkhwa, 2021.
Table 2. Seasonal and gender-wise prevalence of CL in Khyber Pakhtunkhwa, 2021.
DistrictSeasonn%p ValueFemaleMale
n%p Valuen%p Value
BannuWinter21734.3<0.0018938.9< 0.00112831.7<0.001
Spring17227.2 6729.3 10526.0
Summer17026.9 5122.3 11929.5
Autumn7411.7 229.6 5212.9
CharsaddaWinter3757.80.211741.20.4673063.80.058
Spring2742.2 1058.8 1736.2
Summer00.0 00.0 00.0
Autumn00.0 00.0 00.0
ChitralWinter10268.0<0.0011477.80.0188866.7<0.001
Spring4328.7 422.2 3929.5
Summer00.0 00.0 00.0
Autumn53.3 00.0 53.8
DI KhanWinter9644.7< 0.0014147.1<0.0015543.0<0.001
Spring4721.9 1618.4 3124.2
Summer2813.0 1213.8 1612.5
Autumn4420.5 1820.7 2620.3
Dir LowerWinter22152.5< 0.0019555.2<0.00112650.6<0.001
Spring20.5 21.2 00.0
Summer13131.1 5733.1 7429.7
Autumn6715.9 1810.5 4919.7
Lakki MarwatWinter13939.3<0.0014939.20.0029039.3<0.001
Spring7822.0 3024.0 4821.0
Summer6418.1 2217.6 4218.3
Autumn7320.6 2419.2 4921.4
MalakandWinter6927.6<0.0013633.0< 0.0013323.4<0.001
Spring10742.8 4339.4 6445.4
Summer4317.2 1311.9 3021.3
Autumn3112.4 1715.6 149.9
MardanWinter13432.6<0.0015636.6< 0.0017830.2<0.001
Spring8921.7 4026.1 4919.0
Summer13031.6 4328.1 8733.7
Autumn5814.1 149.2 4417.1
NoshehraWinter8637.7<0.0014043.50.0014633.8<0.001
Spring3414.9 2021.7 1410.3
Summer4218.4 1617.4 2619.1
Autumn6628.9 1617.4 5036.8
PeshawarWinter3025.2<0.001266.70.5642824.1<0.001
Spring5747.9 133.3 5648.3
Summer119.2 00.0 119.5
Autumn2117.6 00.0 2118.1
SwatWinter1329.50.436624.00.706736.80.443
Spring1022.7 728.0 315.8
Summer1431.8 832.0 631.6
Autumn715.9 416.0 315.8
TankWinter299100.0 155100.0 144100.0
Spring00.0 00.0 00.0
Summer00.0 00.0 00.0
Autumn00.0 00 00.0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ullah, W.; Yen, T.-Y.; Niaz, S.; Nasreen, N.; Tsai, Y.-F.; Rodriguez-Vivas, R.I.; Khan, A.; Tsai, K.-H. Distribution and Risk of Cutaneous Leishmaniasis in Khyber Pakhtunkhwa, Pakistan. Trop. Med. Infect. Dis. 2023, 8, 128. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed8020128

AMA Style

Ullah W, Yen T-Y, Niaz S, Nasreen N, Tsai Y-F, Rodriguez-Vivas RI, Khan A, Tsai K-H. Distribution and Risk of Cutaneous Leishmaniasis in Khyber Pakhtunkhwa, Pakistan. Tropical Medicine and Infectious Disease. 2023; 8(2):128. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed8020128

Chicago/Turabian Style

Ullah, Wasia, Tsai-Ying Yen, Sadaf Niaz, Nasreen Nasreen, Yu-Feng Tsai, Roger Ivan Rodriguez-Vivas, Adil Khan, and Kun-Hsien Tsai. 2023. "Distribution and Risk of Cutaneous Leishmaniasis in Khyber Pakhtunkhwa, Pakistan" Tropical Medicine and Infectious Disease 8, no. 2: 128. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed8020128

Article Metrics

Back to TopTop