The effect of geo-climatic determinants on the distribution of cutaneous leishmaniasis in a recently emerging focus in eastern Iran

Cutaneous leishmaniasis (CL) has been reported in recent years in South Khorasan Province, a desert region of eastern Iran, where the main species is Leishmania tropica. Little is known of the influence of geography and climate on its distribution, and so this study was conducted to determine geo-climatic factors by using geographic information system. The home addresses of patients with CL patients who were diagnosed and notified from 2009 to 2017 were retrieved from the provincial health center and registered on the village/town/city point layer. The effects of mean annual rainfall (MAR) and mean annual humidity (MAH), mean annual temperature (MAT), maximum annual temperature (MaxMAT), minimum annual temperature (MinMAT), mean annual number of high-velocity wind days (MAWD), mean annual frosty days (MAFD) and snowy days (MASD), elevation, soil type and land cover on CL distribution were examined. The geographical analysis was done using ArcMap software, and univariate and multivariate binary logistic regression were applied to determine the factors associated with CL. A total of 332 CL patients were identified: 197 (59.3%) male and 135 (40.7%) female. Their mean age was 29.3 ± 2.1 years, with age ranging from 10 months to 98 years. CL patients came from a total of 86 villages/towns/cities. By multivariate analysis, the independent factors associated with increased CL were urban setting (OR = 52.102), agricultural land cover (OR = 3.048), and MAWD (OR = 1.004). Elevation was a protective factor only in the univariate analysis (OR = 0.999). Soil type, MAH, MAT, MinMAT, MaxMAT, and MAFD did not influence CL distribution in eastern Iran. The major risk zones for CL in eastern Iran were urban and agricultural areas with a higher number of windy days at lower altitudes. Control strategies to reduce human vector contact should be focused in these settings.