J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 41 http://jad.tums.ac.ir Published Online: March 18, 2018 Original Article Modeling of Environmental Factors Affecting the Prevalence of Zoonotic and Anthroponotic Cutaneous, and Zoonotic Visceral Leishmaniasis in Foci of Iran: A Remote Sensing and GIS Based Study Abdol Ali Golpayegani 1, 2, 3, Ali Reza Moslem 4, Amir Ahmad Akhavan 1, 5, Azam Zeydabadi 1, 2, *Amir Hossein Mahvi 1, 3, 6, *Ahmad Allah-Abadi7 1School of Public Health, Tehran University of Medical Sciences, Tehran, Iran 2Bam University of Medical Sciences, Bam, Iran 3National Institute of Health Research, Ministry of Health, Tehran, Iran 4Sabzevar University of Medical Sciences, Sabzevar, Iran 5Institute for Environmental Research, Tehran University of Medical Sciences, Tehran, Iran 6Center for Solid Waste Research, Institute for Environmental Research, Tehran University of Medical Sciences, Tehran, Iran 7Leishmaniasis Research Center, Sabzevar University of Medical Sciences, Sabzevar, Iran (Received 31 Jan 2016; accepted 10 Feb 2018) Abstract Background: Leishmaniasis is a re-emerging serious international public health problem, and both visceral and cu- taneous types of leishmaniasis became important endemic diseases in Iran. In this study, the relationships between environmental factors (vegetation and elevation) and the prevalence of diseases have been investigated. Methods: All international and national online databases were searched by terms such as leishmaniasis, incidence, prevalence and other related words attributed to Iran and published until first quarter of 2015. The developed database in Excel, later imported to the ArcMap for spatial analyst and mapping. Afterwards, the software was used for modeling the relationship between the prevalence/incidence and environmental variables (vegetation and elevation) by both linear and nonlinear regression. Results: After mapping the prevalence data from 144 studies, considering non-parametric ANOVA, the tendency of zoonotic visceral leishmaniasis to presence in high elevation and high vegetation was more than Anthroponotic and zoonotic cutaneous leishmaniasis. While linear regression showed weaker results for modeling, however, additive nonparametric regression analysis suggested that 10km buffers for elevation, and 10 as well as 50km buffers for veg- etation could contribute in better fitness in modeling of these variables. Conclusion: The detailed maps for distribution of disease concluded. The nonlinear regression is a reliable predictor of the relationship between environmental factors and disease incidence, although more and wide researchers are needed to confirm it. Keywords: Leishmaniasis, GIS, Remote sensing, Environmental variables, Nonlinear regression Introduction Leishmaniasis –a re-emerging serious in- ternational public health problem- is considered as one of the most neglected tropical diseases. The disease is distributed in new foci due to in- fluence of many risk factors, including envi- ronmental factors (1-4). Almost 200000 to 400000 new cases of visceral leishmaniasis (VL), 0.7 million to 1.3 million new cases of cu- taneous leishmaniasis (CL) and 20000 to 30000 deaths annually" are the estimations of WHO last report about importance of the disease (2). Two dominant types of leishmaniasis there in Iran are: zoonotic visceral or kala-azar (ZVL), the most serious and fatal type of the disease, and both types of cutaneous leishmaniasis namely urban form or anthroponotic cutaneous *Corresponding authors: Dr Amir Hossein Mahvi, E- mail: ahmahvi@yahoo.com, Dr Ahmad Allah-Abadi, E- mail: ahmad_health@yahoo.com http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 42 http://jad.tums.ac.ir Published Online: March 18, 2018 leishmaniasis (ACL) and rural form or zoon- otic cutaneous leishmaniasis (ZCL). Leishmania donovani complex are the causes of ZVL pro- tozoan disease which can lead to patient mor- tality (5). Evidence indicates an increase in the incidence of disease in the Old and New Worlds at the early years of the present cen- tury. In recent decades, ZVL became an im- portant endemic disease in the Northwest (6), South and Southwest (7), and Northeast (8) of Iran, especially in nomads. The majority of ZVL cases (92.8%) in country were found among children with age under 12yr old (9). The most recently, direct agglutination test (DAT) as a simple, affordable and Practical test, has found its way as a special quick and accurate test. Thus it has been recommended and applied widely for clinical diagnosis and seroprevalence studies of ZVL in zoonosis and anthroponotic reservoirs in endemic and nonendemic areas of Iran (10). Due to ecologically dependence of the disease vectors and reservoirs to envi- ronment, outbreaks mostly are common in low- lands rural or suburban areas (less than 600M elevation), with a heavy annual rainfall, a mod- erate temperature (15–38 °C), and almost dense vegetation (2). Countries of Mediterranean coasts (Algeria and Syrian Arab Republic), South Americas (Brazil and Colombia), and Middle East (Is- lamic Republic of Iran and Afghanistan) and Central Asia countries are the most reported foci for CL and the major contribution of CL occur in this area (2). The CL cases in Iran dur- ing the period of 2001 to 2008 dramatically has been increased more than two-fold and also over half of 31 provinces of country have CL endemic foci (11). As reported by Center for Disease Control, the annual incidence of leish- maniasis in the country due to referring to health centers has reached about 20000 cases. While the actual cases of disease may be far more than these reports that may even reach up to five times (12). As for the first time carried by John Snow in 1854 due to diarrheal disease detection in London, spatial modeling as a useful tool in epidemiological studies could be used. To- day, Geographic Information System (GIS) in relation to the processed images from remote sensing (RS) has been used to map the geo- graphical distributions of disease prevalence, factors affecting the transmission and control of disease, and the spatial modeling of envi- ronmental factors that have meaningful im- pacts on disease occurrences (9, 13). Vector control programs are associated with logistic problems, mainly due to poor understanding of vector ecology (14). Expanding inexpensive and effective models of integrated management for better control of CL is a critical target which need a thorough will, alongside epidemiology studies and understanding of the ecology of the disease (11) and the role of social, ecological and specially environmental variables for risk factors such as the vegetation and elevation (two most widespread and accessible environ- mental variables). The relationship between ep- idemiological components of leishmaniasis (in- cidence, prevalence, and seroprevalence) and environmental variables (vegetation, elevation, and/or other factors) can be used in GIS mod- eling and may be useful for secondary uses in- cluding prediction. Therefore, contributing to better planning for control activities and the establishment of early warning systems, incor- porating remotely sensed information in epide- miological studies can help public health pol- icy makers in having a better notion and com- prehensive understanding to find the relation- ships between disease and its environmental risk factors (15). The remotely sensed Normalized Differ- ence Vegetation Index (NDVI) –the predom- inantly index which used for vegetation cov- erage– has found its widespread applications due to fluctuations it has along with meteor- ological and environmental components (14). This index together with the other remotely sensed data (11, 16) like the Digital Elevation Model (DEM), have been used extensively in monitoring some vector-borne diseases in- http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 43 http://jad.tums.ac.ir Published Online: March 18, 2018 cluding leishmaniasis throughout the world. While many studies have surveyed the asso- ciations between the mean of NDVI and DEM with incidence or prevalence rate of a type of leishmaniasis locally, however, they do not study all variables in a vast area. In this study, we systematically reviewed all studies regarding the incidence and/or prev- alence rate of all categories of leishmaniasis investigated in Iran from 1976 to 2015. After mapping the results, using Spatial Analysis and Regression techniques, statistical association between variables have been modeled in detail. Materials and Methods To identify the incidence and prevalence of leishmaniasis in endemic and nonendemic area of ZVL, ZCL, and ACL of Iran, we searched and reviewed the literature based on international database including MEDLINE, ISI, SCOPUS, and Google Scholar, and nation- al databases such as SID, IRANDOC, MagIran, and ISC. Search protocols have been based on using the terms related to all types of leishman- iasis, fauna, vectors, reservoirs, incidence, prev- alence, seroprevalence, Direct Agglutination Test (DAT), and their Persian equivalents with the geographic term "Iran". Overall, 503 articles published from 1976 through Jun 2015 in English, French, and Per- sian were founded in first search from data- bases. Besides, by referring to CV of famous professors, abstract seminars, thesis, non-dig- ital papers, and the reference/citations of arti- cles, 121 additional sources were found. After the elimination of duplicate and inappropriate articles by two separate teams, we consequent- ly included 144 studies that clearly addressed confirmed clinical cases for human seroprev- alence (DAT ≥ 1: 3200)/ incidence of ZVL and active lesion prevalence/total prevalence (active lesion + scar)/incidence of ACL and ZCL in Iran. Disease data were extracted by researchers directly from body text of arti- cles to database tables (Dbase format using the MS Excel software) containing all cities (1246 point) and villages with a population more than 500 people (12034 points), disease information, and descriptive data. Depending on the raw data presented in each article, one or more than one category of epidemiologic data was included in the main table. These numeric values (all in percent) for ACL and ZCL includ- ing active lesion prevalence, total prevalence (active lesion + scar), and incidence and for VCL including seroprevalence and incidence were directly or indirectly extracted using study population. The developed database later im- ported to the ArcMap® software version 9.3 (ESRI, Redlands, USA) software platform for Spatial Analyst and mapping all types of leish- maniasis. R software version 3.2.2 (17) was used for analyzing the relationship between the prevalence/incidence and continues spatial var- iables. Shapefile maps for political subdivisions of provinces, counties, rural districts, and points of cities, towns, and villages with their attribute tables were obtained from National Cartograph- ic Center, Iran. The Moderate Resolution Im- aging Spectroradiometer (MODIS) onboard the Terra satellite, is one of most applied NDVI images source in literature. NDVI calculated from surface reflectance in the near-infrared ρNIR and red ρR portions of the electromag- netic spectrum, using the equation 1 (18): 1) NDVI= (ρNIR–ρR)/ (ρNIR+ρR) The 16d, 250m NDVI raster maps (MOD 13Q1 series, version 5) projected in Geograph- ic Coordinate System WGS_1984, derived from MODIS Reprojection Tool (MRT) (19) for the years 2000 to 2014. Sixteen days maps unified to annual average NDVI in Spatial Analysis, then reclassified to 1–20. Altitude was based on DEM, obtained from the 100m ASTER Global Digital Elevation Model (ASTER GDEM) da- ta center Version 2 (20). These maps then re- projected to Lambert Conformal Conic in ArcMap, and additional Raster Analysis such as reclassification, annual and 15yr average http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 44 http://jad.tums.ac.ir Published Online: March 18, 2018 and zonal statistics were performed on them. Overall, 5km, 10km, 20km, and 50km buffer zones centered on each city/village point re- ported values for leishmaniasis, and their zonal Min, Max, Minority, Majority, Median and Mean statistics were calculated and exported in Dbase tables. Finally, massive joined ta- bles were called in R, and statistical calcula- tions were completed for each. Data analysis CSV tables for any purposes derived from the main disease data table and called in R console. We chose to test the normality of data with One-sample Kolmogorov-Smirnov and Shapiro-Wilk normality test. Next, we mod- eled the relationships between prevalence/ incidence values and the continuous independ- ent variables. Two regression scenarios were considered. First, using simple linear regres- sion models, a stepwise backward elimination procedure was used with the variables, final- ly resulting in the model contains merely the variables which associated with the outcome. Secondly, using generalized additive models -which use spline or local regression- a step- wise backward elimination procedure, was used again to reach the final model. The additive nonparametric regression model is presented in equation 2 (21): 2) y= β0+ m1(x1)+ m2(x2)+ . . .+ mk(xk)+ ɛ Where the partial-regression functions mj are fit using a simple-regression smoother, such as local polynomial regression or smoothing splines. In nonparametric regression, in con- trast, the object is to estimate the regression function directly without specifying its form explicitly. In nonparametric comparing with simple regression, the ultimate goal is to ob- tain the regression function directly from in- ternal computation without any open aspect (22). The level of 5% significance for all tests was considered. Results We found 144 studies that directly in- formed the active lesion/scar prevalence and/or incidence of leishmaniasis in Iran from 1976– 2015 (Table 1). These data considered as a base for construction of 89 ACL points including 15, 19, and 71 urban/rural points for active le- sion, total prevalence (active lesion + scar), and incidence respectively. The same way, data for 420 ZCL points included 67, 68, and 362 points. Cases for 355 ZVL points for the seropreva- lence and incidence were 283 and 82 respec- tively. Figures 1, 2, and 3 show the maps of three types of leishmaniasis in Iran separate- ly. The DEM map and 15yr average of NDVI map also attached on VL and ACL maps re- spectively. Rating in the categories are High, Middle and Low for prevalence data and ob- tained based on expert consensus. Scatter diagrams and normality plot of each independent variable and the disease incidence/ prevalence were created and visually inspected before running regression. Some of these plots traced in Fig. 4. Normal probability of data in our variables has been questioned, then, the normality of data statistically tested. Therefore, using both One-sample Kolmogorov-Smirnov and Shapiro-Wilk normality test, the data are without normality. However, both linear and nonparametric regression fitted on data. To compare of the means of main inde- pendent variables (DEM and NDVI) versus types of leishmaniasis, homogeneity of vari- ances rejected by Fligner-Killeen test, but due to nonparametric Kruskal-Wallis test, we result in significant difference between means of var- iables. Regard to boxplots for comparison of means in Fig. 5, and nonparametric analysis of variances in Table 2, tendency of ZVL to presence in high elevation and high vegetation is more than ACL and ZCL. Similarly, we com- pare the means of types of leishmaniasis oc- currences and epidemiological variables of disease in Table 3. http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 45 http://jad.tums.ac.ir Published Online: March 18, 2018 The regression models were evaluated us- ing R software packages. Tables 4 and 5 gen- erally descript data structure and assessment of the effect of, or relationship between, ex- planatory variables on the response by using multiple linear regression models. To determine the best final fitted models for each group of data, based on P-value, ad- justed R-squared, and AIC (Akaike's 'An In- formation Criterion'), using stepwise back- ward simple linear regression, following model produced: (although second order or in- teraction models had a little more R-squared, due to the complexity were excluded). 3) Total prevalence= 12.14– 0.014 (MEAN$ DEM20)+ 0.005 (MEAN$DEM50)– 1.64 (MIN$NDVI20)+ 1.8 (MEAN$NDVI20) Multiple R-squared: 0.2111, Adjusted R-squared: 0.198 F-statistic: 16.06 on 4 and 240 DF, P: 1.151e–11 Important plots of main model presented in Fig. 6. Regarding weak results of linear re- gression (R-squared located near the zero), eventually, Additive Nonparametric Regres- sion analysis, using GAM (1Generalized Addi- tive Model) command in library “mgcv”, ap plied to find best parameters and buffers for modeling variables of this study. Here two models were fitted: first with the highest per- centage of predictability by equation 4, sec- ond with the convenience practical application (same parameter mean of variables in unique buffer 10km) which came in equation 5. 4) Total prevalence= 7.69+ s (MEAN$DEM10)+ s (MAX$NDVI50) R-sq. (adj)= 0.524 Deviance explained= 54.5% 5) Total prevalence= 7.42+ s (MEAN$DEM10)+ s (MEAN$NDVI10) R-sq. (adj)= 0.371 Deviance explained= 39.6% The “s” function, used in specifying the model formula, indicates that each term is to be fitted with a smoothing spline. The sur- face that fitted on the model and smoothing splines have been shown in Fig. 7 and 8. In equation 5, the 5.43 parameters are used for the MEAN$NDVI10 (mean of 10km buffer for NDVI) term, and 7.98 for the MEAN$DEM10 term, the degrees of freedom for the model is the sum of these plus 1 for the regression constant. Fig. 1. Distribution of anthroponotic cutaneous leishmaniasis city points in Iranian studies from 1976–2013 on 15yr annual average Normalized Difference Vegetation Index map. High: active lesion ≥ 2% or active lesion+scar ≥ 20% or incidence ≥ 2%; Middle: 2% ˃ active lesion ≥ 0.5% or 20% ˃ active lesion+scar ≥ 5% or 2 ˃ incidence ≥ 0.5%, Low: 0.5% ˃ active lesion or 5% ˃ active lesion+scar or 0.5˃ incidence http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 46 http://jad.tums.ac.ir Published Online: March 18, 2018 Table 1. Literatures published for prevalence of leishmaniasis, divided by leishmaniasis types and provinces of Iran Province ACL references ZCL references ZVL references Alborz Fagih-Nayini et al. 2002(23) Ardabil Mazloumi Gavgani et al. 2002(24), Soleimanzadeh et al. 1993(25), Mohebali et al. 2010(26), Edrissian 1996(27), Mohebali et al. 2011(28), Mahami et al. 2006(29), Mohammadi-Kheyrabadi et al. 2004(30), Nadim et al. 1998(31), Arshi et al. 2002(32), Mazloumi Gavgani et al. 2011(33) Azerbaijan, East Mazloumi Gavgani et al. 2002(24), Mohebali et al. 2011(28), Mirsamadi et al. 2003(34), Mazloumi Gavgani et al. 2011(33) Azerbaijan, West Bushehr Yaghoobi-Ershadi et al. 2013(35) Yaghoobi-Ershadi et al. 2013(35), Hamzavi et al. 2000(36) Mohebali et al. 2001(37) Chahar Mahaal and Bakhtiari Fars Ghatee et al. 2013(38) Razmjou et al. 2009(39), Ghatee et al. 2013(38), Rassi et al. 2004(40), Fakoorziba et al. 2011(41), Davami et al. 2010(42), Sharafi et al. 2013(43), Noorpisheh et al. 2013(44), Entezari and Eskandari 2014(45) Fakhar et al. 2008(46), Edrissian et al. 1993(47), Sahabi et al. 1992(48), Fakhar et al. 2006(49) Gilan Golestan Mollalo et al. 2015(50), Rahbarian et al. 2009(51), Sofizadeh et al. 2012(52), Baghaei et al. 2012(53), Mesgarian et al. 2010(54) (Fakhar et al. 2014) Hamadan Hanafi Bojd et al. 2006(55), Nazari 2012(56) Hormozgān Azizi et al. 2012(57) Ilam Asgari Nezhad et al. 2012(58), Shazad et al. 2005(59), Bahrami et al. 2011(60), Yazdanpanah and Rostamianpur 2013(61), Roghani et al. 2012(62) Isfahan Shiee et al. 2012(63), Zahraei-Ramazani et al. 2007(64), Doroodgar et al. 2012(65) Talari et al. 2006(66), Yaghoobi- Ershadi et al. 2001(67), Mohammadian et al. 1999(68), Yaghoobi-Ershadi and Javadian 1995(69), Yaghoobi-Ershadi et al. 2006(70), Nilforoushzadeh et al. http://jad.tums.ac.ir/ https://en.wikipedia.org/wiki/Alborz_Province https://en.wikipedia.org/wiki/Ardabil_Province https://en.wikipedia.org/wiki/East_Azerbaijan_Province https://en.wikipedia.org/wiki/West_Azerbaijan_Province https://en.wikipedia.org/wiki/West_Azerbaijan_Province https://en.wikipedia.org/wiki/Bushehr_Province https://en.wikipedia.org/wiki/Chaharmahal_and_Bakhtiari_Province https://en.wikipedia.org/wiki/Chaharmahal_and_Bakhtiari_Province https://en.wikipedia.org/wiki/Fars_Province https://en.wikipedia.org/wiki/Gilan_Province https://en.wikipedia.org/wiki/Golestan_Province https://en.wikipedia.org/wiki/Hamadan_Province https://en.wikipedia.org/wiki/Hormozgan_Province https://en.wikipedia.org/wiki/Ilam_Province https://en.wikipedia.org/wiki/Isfahan_Province J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 47 http://jad.tums.ac.ir Published Online: March 18, 2018 Table 1. Continued … 2002(71), Doroodgar et al. 2009(72), Ahmadi et al. 2013(73), Ebadi and Hejazi 2003(74), Yaghoobi-Ershadi et al. 1999(75), Doroudgar et al. 1996(76) Kerman Nadim and Aflatoonian 1995(77), Sharifi et al. 2011(78), Sharifi et al. 2015(3), Yaghoobi-Ershadi 1977(79), Mirzaei et al. 2012(80), Aflatoonian et al. 2014(81), Sharifi et al. 1998(82), Sharifi et al. 2011(83), Ghatee et al. 2013(38), Aflatoonian et al. 2013(84), Nadim et al. 1995(85), Sharifi et al. 2015(86), Aflatoonian and Sharifi 2011(87), Aflatoonian and Sharifi 2010(88), Mirzadeh et al. 2008(89), Aflatoonian et al. 2014(90), Aflatoonian and Sharifi 2007(91) Sharifi et al. 2015(3), Akhavan et al. 2007(92), Aghaei-Afshar et al. 2013(93), Khosravi et al. 2013(94), Sharifi et al. 2008(95) Mahmoudvand et al. 2011(96), Barati et al. 2008(97) Kermanshah Hamzavi and Khademi 2015(98), Nazari et al. 2012(99), Hamzavi and Khademi 2013(100) Hamzavi et al. 2012(101) Khorasan, North Alavinia et al. 2009(102) Alavinia et al. 2009(102) Torabi et al. 2007(103) Khorasan, Razavi Sattar Pagheh et al. 2013(104), Moosa-Kazemi et al. 2007(105), Saadabadi et al. 2013(106), Hassanpour et al. 2014(107), Hoseini Farash et al. 2011(108), Khajedaluee et al. 2014(109), Mohajery et al. 2008(110) Yaghoobi-Ershadi et al. 2003(111), Hassanpour et al. 2014(107), Khajedaluee et al. 2014(109), Karimi Zarchi et al. 2004(112), Sahabi et al. 1999(113), Akbari et al. 2014(114) Khorasan, South Karamian et al. 2013(115) Khuzestan Kassiri et al. 2012(116), Kassiri et al. 2012(117), Kassiri et al. 2014(118), Kassiri et al. 2014(119), Vazirianzadeh et al. 2014(120), Kassiri et al. 2013(121), Spotin et al. 2014(122), Ghasemian et al. 2011(123), Kassiri et al. 2013(124), Nejati et al. 2014(125), Behbahani et al. 2012(126), Kassiri et al. 2011(127) Kohgiluyeh and Boyer-Ahmad Sarkari et al. 2010(128), Sarkari et al. 2007(129) Kurdistan http://jad.tums.ac.ir/ https://en.wikipedia.org/wiki/Kerman_Province https://en.wikipedia.org/wiki/Kermanshah_Province https://en.wikipedia.org/wiki/North_Khorasan_Province https://en.wikipedia.org/wiki/North_Khorasan_Province https://en.wikipedia.org/wiki/Razavi_Khorasan_Province https://en.wikipedia.org/wiki/Razavi_Khorasan_Province https://en.wikipedia.org/wiki/South_Khorasan_Province https://en.wikipedia.org/wiki/Khuzestan_Province https://en.wikipedia.org/wiki/Kohgiluyeh_and_Boyer-Ahmad_Province https://en.wikipedia.org/wiki/Kohgiluyeh_and_Boyer-Ahmad_Province https://en.wikipedia.org/wiki/Kurdistan_Province J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 48 http://jad.tums.ac.ir Published Online: March 18, 2018 Lorestan Kheirandish et al. 2013(130) Amraee et al. 2013(131), Ahmadi et al. 2013(132), Kheirandish et al. 2013(130), Chegeni-Sharafi et al. 2011(133) Chegeni-Sharafi et al. 2005(134) Markazi Mazandaran Qazvin Qom Rostami et al. 2013(135), Rassi et al. 2013(12), Akhavan et al. 2003(136), Saghafipour et al. 2013(137), Saghafipour et al. 2012(138) (Rakhshanpour et al. 2014), (Fakhar et al. 2004) Semnan Mohammadi Azni et al. 2010(139), Mohammadi Azni et al. 2011(140), Mohammadi Azni et al. 2010(141), Rafati et al. 2007(142) Sistan and Ba- luchestan Fazaeli et al. 2008(143), Fazaeli et al. 2009(144), Fouladi et al. 2007(145) Tehran Salehzadeh and Seyedi Rashti 1996(146), Hossein Abadi 2004(147) Yazd Yaghoobi-Ershadi et al. 2002(148) Yaghoobi-Ershadi et al. 2004(149), Jaafary et al. 2007(150), Mozafari and Bakhshi Zadeh 2011(151), Hanafi Bojd et al. 2003(152), Yaghoobi-Ershadi et al. 2008(153) Zanjan Multy-provinces Mohebali et al. 2006(6) Table 2. Nonparametric analysis of variances between types of zoonotic, anthroponotic cutaneous, and zoonotic visceral leishmaniasis and environmental variables in Iran (P˂ 0.05), buffer = 5km DEM_MEAN DEM_VARIANCE NDVI_MEAN* NDVI_VARIANCE ACL 1216 213046 8.57 3.93 ZCL 666 385752 10.12 12.16 ZVL 1296 332153 11.93 6.28 *(The range is 1–20 for NDVI) Table 3. Nonparametric Analysis of variances between types of zoonotic, anthroponotic cutaneous, and zoonotic visceral leishmaniasis and epidemiological variables of disease in Iran (P˂ 0.05) MEAN VARIANCE ACL (%) ZCL (%) ZVL (%) ACL ZCL ZVL Active lesion prevalence 1.48 4.42 Na* 3.47 37.97 Na Total prevalence 11.88 25.02 2.69 125.38 622.39 3.73 Incidence 0.76 1.04 1.01 2.62 8.47 1.78 *Na: Not available Table 1. Continued … http://jad.tums.ac.ir/ https://en.wikipedia.org/wiki/Lorestan_Province https://en.wikipedia.org/wiki/Markazi_Province https://en.wikipedia.org/wiki/Mazandaran_Province https://en.wikipedia.org/wiki/Qazvin_Province https://en.wikipedia.org/wiki/Qom_Province https://en.wikipedia.org/wiki/Semnan_Province https://en.wikipedia.org/wiki/Sistan_and_Baluchestan_Province https://en.wikipedia.org/wiki/Sistan_and_Baluchestan_Province https://en.wikipedia.org/wiki/Tehran_Province https://en.wikipedia.org/wiki/Yazd_Province https://en.wikipedia.org/wiki/Zanjan_Province J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 49 http://jad.tums.ac.ir Published Online: March 18, 2018 Table 4. Multiple linear regression models for types of zoonotic, anthroponotic cutaneous, and zoonotic visceral leishmaniasis and Digital Elevation Model subcategory (buffer distance around the cities or villages) DEM 5km_total DEM 5km_ACL DEM 5km_ZVL DEM 5km_ZCL Buffer types Buffer types_ACL Buffer types_ZVL Buffer types_ZCL Active lesion prevalence ~MEAN+MAX AR*: 0.1004 ~MEAN+MIN AR: 0.284 Na*** ~MEAN AR: 0.077 ~MEAN20km AR: 0.035 ns ns ~MEAN5km AR: 0.077 Total prevalence ~MAX AR: 0.1028 ~MEAN AR:0.1649 ~MAX AR:0.0484 ~MEAN+MAX AR: 0.234 ~MEAN20km AR: 0.1178 ns ~MEAN50km AR: 0.029 ~MEAN20km AR: 0.25 Incidence ~MEAN+MAX AR: 0.029 Ns** ~MIN AR:0.2954 ~MAX+MAJORITY AR: 0.077 ~MEAN5km AR: 0.011 ns ~MEAN20km AR: 0.21 ~MEAN5 AR: 0.031 *AR: Adjusted R-squared **Ns: Not significant ***Na: Not available Fig. 2. Distribution of zoonotic cutaneous leishmaniasis rural region points in Iranian studies from 1991–2013. High: active lesion ≥2% or active lesion+scar ≥20% or incidence ≥ 2%, Middle: 2% ˃ active lesion ≥ 0.5% or 20% ˃ active lesion+scar ≥ 5% or 2˃ incidence ≥ 0.5%, Low: 0.5% ˃ active lesion or 5% ˃ active lesion+scar or 0.5 ˃ inci- dence Fig. 3. Distribution of zoonotic visceral leishmaniasis county region points in Iranian studies from 1989–2012 on Digital Elevation Model map. High: DAT ≥ 2% or incidence ≥ 1%, Middle: 2% ˃ DAT ≥ 0.5% or 1 ˃ incidence ≥ 0.1%, Low: 0.5% ˃ DAT or 0.1 ˃ incidence http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 50 http://jad.tums.ac.ir Published Online: March 18, 2018 Table 4. Multiple linear regression models for types of zoonotic, anthroponotic cutaneous, and zoonotic visceral leishmaniasis and Digital Elevation Model subcategory (buffer distance around the cities or villages) DEM 5km_total DEM 5km_ACL DEM 5km_ZVL DEM 5km_ZCL Buffer types Buffer types_ACL Buffer types_ZVL Buffer types_ZCL Active lesion prevalence ~ MEAN+MAX AR*: 0.1004 ~ MEAN+MIN AR: 0.284 Na*** ~ MEAN AR: 0.077 ~ MEAN20km AR: 0.035 ns ns ~ MEAN5km AR: 0.077 Total prevalence ~ MAX AR: 0.1028 ~ MEAN AR:0.1649 ~ MAX AR:0.0484 ~ MEAN+MAX AR: 0.234 ~MEAN20km AR: 0.1178 ns ~MEAN50km AR: 0.029 ~MEAN20km AR: 0.25 Incidence ~MEAN+MAX AR: 0.029 Ns** ~ MIN AR:0.2954 ~ MAX+MAJORITY AR: 0.077 ~ MEAN5km AR: 0.011 ns ~ MEAN20km AR: 0.21 ~ MEAN5 AR: 0.031 *AR: Adjusted R-squared **Ns: Not significant ***Na: Not available 0 500 1000 1500 2000 2500 0 5 1 0 1 5 2 0 2 5 3 0 3 5 plot of Active Ulcer Prevalence and 5kmDEM MEAN of DEM U lc e r P re v a le n c e o f li e s h m a n ia s is -3 -2 -1 0 1 2 3 0 5 0 0 1 0 0 0 2 0 0 0 normal probability plot of Lieshmaniasis Incidence and 5kmDEM mean Theoretical Quantiles S a m p le Q u a n ti le s -3 -2 -1 0 1 2 3 0 5 0 0 1 0 0 0 1 5 0 0 2 0 0 0 normal probability plot of ZCL prevalence and 5kmDEM mean Theoretical Quantiles S a m p le Q u a n ti le s -3 -2 -1 0 1 2 3 0 2 0 4 0 6 0 8 0 1 0 0 normal probability plot of total prevalence and 5kmNDVI mean Theoretical Quantiles S a m p le Q u a n ti le s Fig. 4. Typical plots and normality plots of independent and dependent variables ACL VL ZCL 0 5 0 0 1 0 0 0 2 0 0 0 Boxplot of leishmaniasis types Vs MEAN of DEM leishmaniasis type M e a n o f 5 k m B u ff e r D E M ACL VL ZCL 5 1 0 1 5 2 0 Boxplot of leishmaniasis types Vs MEAN of NDVI leishmaniasis type M e a n o f 5 k m B u ff e r 1 5 y e a r N D V I Fig. 5. Boxplots for compares of means of 5km-buffer Digital Elevation Model and Normalized Difference Vegeta- tion Index ((in the range of 1–20)) for anthroponotic cutaneous leishmaniasis, zoonotic cutaneous leishmaniasis, and zoonotic visceral leishmaniasis in Iran for study period (1976–2013) http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 51 http://jad.tums.ac.ir Published Online: March 18, 2018 -10 0 10 20 30 -2 0 2 0 6 0 Fitted values R e s id u a ls Residuals vs Fitted 900 905 151 -3 -2 -1 0 1 2 3 0 2 4 6 Theoretical Quantiles S ta n d a rd iz e d r e s id u a ls Normal Q-Q 900 905 151 -10 0 10 20 30 0 .0 1 .0 2 .0 Fitted values S ta n d a rd iz e d r e s id u a ls Scale-Location 900 905 151 0.00 0.02 0.04 0.06 0.08 0.10 -2 0 2 4 6 Leverage S ta n d a rd iz e d r e s id u a ls Cook's distance 0.5 Residuals vs Leverage 900 870 880 Fig. 6. Plots of final linear regression model for leishmaniasis and environmental variables (Digital Elevation Model and Normalized Difference Vegetation Index) m ean of 10km D E M 0 1000 2000 3000 m ea n of 1 0k m N D V I 5 10 15 20 fitte d P re v a le n c e -20 0 20 40 gam(tot_Prev~DEM$MEAN10+NDVI$MEAN10) Fig. 7. Fitted surface for the additive nonparametric regression of total prevalence of leishmaniasis on mean of 10km buffer for Digital Elevation Model and Normalized Difference Vegetation Index http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 52 http://jad.tums.ac.ir Published Online: March 18, 2018 Table 5. Multiple linear regression models for disease mode and Normalized Difference Vegetation Index subcate- gory (buffer distance around the cities or villages) NDVI 5km_total NDVI 5km_ACL NDVI 5km_ZVL NDVI 5km_ZCL Buffer types Active lesion prevalence ~MEAN+MAJORITY AR*: 0.0338 **Ns ***Na Ns ~MAX50km AR: 0.093 Total prevalence ~MIN AR: 0.012 Ns ~MEAN AR:0.39 ~ MIN + MI- NORITY AR: 0.3479 ~ MIN5km AR: 0.012 Incidence ~MEAN AR: 0.0085 ~ MINORITY AR: 0.042 Ns ~MEDIAN AR: 0.0259 ~MEAN20km AR: 0.007 *AR: Adjusted R-squared **Ns: Not significant ***Na: Not available 0 500 1000 1500 2000 2500 3000 3500 -3 0 -2 0 -1 0 0 1 0 2 0 3 0 4 0 prevalence of leishmaniasis on mean of 10km buffer Digital Elevation Model MEAN of 10km DEM s (M E A N $ D E M 1 0 ,7 .9 8 ) 5 10 15 -3 0 -2 0 -1 0 0 1 0 2 0 3 0 4 0 prevalence of leishmaniasis on mean of 10km buffer Normalized Difference Vegetation Index Model MEAN of 10km NDVI s (M E A N $ N D V I1 0 , 5 .4 3 ) Fig. 8. Partial-regression functions for the additive regression of total prevalence of leishmaniasis on mean of 10km buffer Digital Elevation Model and Normalized Difference Vegetation Index. The broken lines give point-wise 95- percent confidence envelopes around the fit Discussion Maps based on the actual data for distribu- tion of all leishmaniasis types in Iran, with iden- tifying the detailed variables of environmental factors affecting the disease distribution, are the most important findings of this study. Due to the limited access to official leishmaniasis data of Center for Disease Control, integration all data from both CL types and rural and urban foci in the heart of the county official reports (lack of separation species), tens of kilometers distance between some foci and center of coun- ties often identified as disease points (lack of spatial resolution), unknown types of CL in most of official data, and passive screening in Public Health Centers, we decided to use the real data available in the literatures with a de- tailed focus on disease points. Many of scien- tific epidemiological studies predominantly fo- cus on foci of diseases and have very valuable and verified information about all aspects of http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 53 http://jad.tums.ac.ir Published Online: March 18, 2018 the disease, such as prevalence, incidence, vec- tors, reservoirs, and control of disease data. Maps that contain all of the scientific infor- mation investigate by researchers of a country over several decades could well represent the disease situation in that country. Although the nature of the leishmaniasis itself was not the aim of our study, however, distribution of dis- ease to the breakdown mapping can present more accurate view of disease for researchers. Though vegetation statistically has a sig- nificant relationship with the prevalence of vector-borne diseases, this study (especially for the simple linear regression models), with a weak confirmation, implicitly point to a com- plex and indirect relationship. Therefore, more detailed studies and comprehensive investiga- tions would be needed to prove the theorem. A study on effect of minimum, maximum and mean NDVI on the distribution of the VL vec- tor in district of Bihar, India showed that max and mean variables with an R2 about 0.6 strong- ly associated with the disease (14). Werneck et al. with a mention to the NDVI, emphasized on the impact of multilevel variable modeling on the development of leishmaniasis in Bra- zil (154). Study of environmental factors in- fluencing the distribution of vectors involved in transmission of disease in north-eastern It- aly revealed that areas with high winter NDVI may be related to the survival of the larvae in moist soil (155). According to the results of relationship between CL and NDVI in north- east province of Iran, Golestan, arid and semi- arid regions with low vegetation are the most involved area for CL occurrence (11). In our study, as expected due to their rural and no- madic natures of ZCL and ZVL in Iran, the values for NDVI index increase for ACL, ZCL, and ZVL respectively (Table 2). There is not much variation among them, and all three types are almost inclined to intermediate vegetation (neither too sparse nor too dense). The distinction in the result of numerous studies would be aroused from differences between climate factors, socio-economic variables, in- strumental confounders and differences in the interpretation of the results. Incidence of various types of leishmaniasis in different heights may result from the diver- sity in the reservoirs and vectors of the dis- ease. Whereas ACL and ZVL had an interest for presence in high altitude areas, ZCL have been distributed in lower altitude (Table 2). This result for ZVL has obtained also in Ilam Province spatial modeling, and elevation had a negative impact on CL prevalence (156). Multivariate analysis in a VL study in eastern Sudan found the elevation as important vari- ables, while the primary analysis did not show any correlation with disease incidence (157). However, while DEM can use as a good var- iable for topographic purposes, it may do not has most appropriate function as elevation in all regions (155). The role of environmental factors in the development of disease vectors is obvious. While the CL is the main vector-borne disease in Iran, about 80% of cases are zoonotic and spread in 17 out of 31 provinces (158). Some studies have been modeled the environmen- tal factors for vectors and reservoirs of leish- maniasis. Altitude and land cover was found to be from most important factors for suitabil- ity of niches for sand flies. Areas with 990m were found and 1235m average altitude were related to probability of presence of ZCL and ACL vectors respectively (159). Niche mod- eling of main reservoir hosts in Iran showed that among topography variables, slope has the most contribution in forecasting risk of disease (160). A predictive degree-day model was used for development time, population dynamics and activities of VL vector sandflies in field in northwest Iran (161). Among climate factors, minimum of temperature, mean of humidity, and rainfall had the most impact on ZVL dis- tribution in Golestan, north-east of Iran (162). None of these studies include NDVI as a var- iable. The vegetation coverage has been dis- cussed as a risk factor for VL and CL types http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 54 http://jad.tums.ac.ir Published Online: March 18, 2018 in previous studies (163, 164). The NDVI has not considered for intrinsic specifications it- self but also it covers some other important environmental and ecological variables such as soil types and moisture, humidity, slope degree and even elevation of its ambient land (16). The majority of studies, evaluated the presence/absence of vegetation or have used a cut-off value point for NDVI, but less at- tention has been paid to the sheer numerical modeling. Linear regression results in this study show that none of subcategories of in- dependent variables have a favorable associ- ated with DEM and NDVI and do not fit on the linear model very well (low R-squared). Only some of models when focused on leish- maniasis types, the R-squared reach to 0.2– 0.4. This fact could be predicted from plots associated with variables (Fig. 4). Despite the lack of goodness of fitness of simple linear models, Additive Nonparametric Regression analysis presents a better prediction model for data of this study. These models do not out- put real intercepts for variables, but using their internal formulas, they are able to provide pre- dictive models for input data. As a nonpara- metric regression attributes, no obvious esti- mates were seen for parameters (21), and then to figure out regression results, regards to ne- cessitate of using fitted regression graphics, we used this method as shown in Fig. 7. Da- ta frame used to find predicted values on the regression surface, fitted to data of model to drawn this figure. For more clarification, functions which ex- plained previously that present in inherent con- tent of the additive regression, have been shown in partial regression of our model in Fig. 8. Conclusion Although the results of modeling and pre- dictive power of the models in this study was not great but was somewhat satisfied. Mod- eling environmental factors which affecting ecological disease, need attention to the role of all these factors together. Moreover, inte- gration the scenarios such as finding hotspots in GIS, using statistical logistic regression, more specific factors of diseases such as vector and reservoir species, can be used in modeling the relationships, more meaningful and more clear. However, a good model needs to consider all the factors involved the prevalence and inci- dence of a disease. Acknowledgements The authors would like to thank financial grant supporting this study by National Insti- tute of Health Research, Tehran University of Medical Sciences, Contract No: 241/M/9644 and also by Sabzevar University of Medical Sciences, as well. The NDVI raster maps (MOD13Q1 series) Version 5 and 100m AS- TER Global DEM Version 2 were retrieved from the online MODIS Reprojection Tool and online ASTER GDEM Data product of METI and NASA respectively, courtesy of the NASA EOSDIS Land Processes Distributed Active Archive Center (LP DAAC), USGS/ Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota. The authors declare that there is no conflict of in- terests. References 1. Cerbino Neto J, Werneck GL, Costa CHN (2009) Factors associated with the in- cidence of urban visceral leishmania- sis: an ecological study in Teresina, Piauí State, Brazil. Cadernos de Saúde Pública. 25(7): 1543–1551. 2. WHO (2016) Leishmaniasis fact sheet no. 375. World Health Organization, Geneva, Switzerland. Available at: http://www.who.int/mediacentre/fact sheets/fs375/en/. 3. Sharifi I, Aflatoonian MR, Fekri AR, Hakimi Parizi M, Aghaei Afshar A, http://jad.tums.ac.ir/ http://www.who.int/mediacentre/factsheets/fs375/en/ http://www.who.int/mediacentre/factsheets/fs375/en/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 55 http://jad.tums.ac.ir Published Online: March 18, 2018 Khosravi A, Sharifi F, Aflatoonian B, Khamesipour A, Dowlati Y, Modabber F, Nadim A (2015) A comprehensive review of cutaneous leishmaniasis in Kerman Province, southeastern iran- narrative review article. Iran J Public Health. 44(3): 299–307. 4. Kassiri H, Sharifinia N, Jalilian M, Shemshad K (2012) Epidemiological aspects of cutaneous leishmaniasis in Ilam Province, west of Iran (2000– 2007). Asian Pac J Trop Med. 2: 382– 386. 5. Mohebali M, Edrissian GH, Shirzadi MR, Akhoundi B, Hajjaran H, Zarei Z, Molaei S, Sharifi I, Mamishi S, Mahmoudvand H, Torabi V, Moshfe A, Malmasi A, Motazedian MH, Fakhar M (2011) An observational study on the current distribution of visceral leishmaniasis in different ge- ographical zones of Iran and impli- cation to health policy. Travel Med Infect Dis. 9(2): 67–74. 6. Mohebali M, Edrissian GH, Nadim A, Hajjaran H, Akhoundi B, Hooshmand B, Zarei Z, Arshi S, Mirsamadi N, Naeini Manouchehri K (2006) Ap- plication of direct agglutination test (DAT) for the diagnosis and seroep- ide-miological studies of visceral leish- maniasis in Iran. Iran J Parasitol. 1 (1): 15–25. 7. Ghatee MA, Sharifi I, Haghdoost AA, Kanannejad Z, Taabody Z, Hatam G, Abdollahipanah A (2013) Spatial cor- relations of population and ecological factors with distribution of visceral leishmaniasis cases in southwestern Iran. J Vector Borne Dis. 50(3): 179– 187. 8. Mohebali M (2012) Epidemiological status of visceral leishmaniasis in Iran: ex- periences and review of literature. Int J Clin Exp Pathol. S3: 003. 9. Masimalai P (2014) Remote sensing and Geographic Information Systems (GIS) as the applied public health and en- vironmental epidemiology. Int J Med Sci Public Health 3(12): 1430–1438. 10. Ready PD (2010) Leishmaniasis emer- gence in Europe. Euro Surveill. 15 (10): 19505–19516. 11. Mollalo A, Alimohammadi A, Shah- risvand M, Shirzadi MR, Malek MR (2014) Spatial and statistical analyses of the relations between vegetation cover and incidence of cutaneous leishmaniasis in an endemic prov- ince, northeast of Iran. Asian Pac J Trop Dis. 4(3): 176–180. 12. Rassi Y, Saghafipour A, Abai MR, Oshaghi MA, Mohebali M, Mostafavi R (2013) Determination of Leishmania Parasite Species of Cutaneous Leish- maniasis Using PCR Method in Cen- tral County, Qom Province. Zahedan J Res Med Sci. 15(12): 13–16. 13. Palaniyandi M (2014) GIS for disease surveillance and health information management in India. Geospatial Today. 13(5): 44–46. 14. Bhunia GS, Kesari S, Chatterjee N, Mandal R, Kumar V, Das P (2012) Seasonal relationship between normal- ized difference vegetation index and abundance of the Phlebotomus kala- azar vector in an endemic focus in Bihar, India. Geospatial Health. 7(1): 51–62. 15. Tsegaw T, Gadisa E, Seid A, Abera A, Teshome A, Mulugeta A, Herrero M, Argaw D, Jorge A, Aseffa A (2013) Identification of environmental param- eters and risk mapping of visceral leish- maniasis in Ethiopia by using geo- graphical information systems and a statistical approach. Geospatial Health. 7(2): 299–308. 16. Palaniyandi M, Anand PH, Maniyosai R (2014) Climate, landscape and the en- vironments of visceral leishmaniasis http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 56 http://jad.tums.ac.ir Published Online: March 18, 2018 transmission in India, using remote sensing and GIS. J Geophys Remote Sens. 3(3): 1–6. 17. R Core Team (2015) R: A Language and Environment for Statistical Compu- ting. R Foundation for Statistical Com- puting. Vienna, Austria. 18. Hereher ME (2015) Environmental mon- itoring and change assessment of Tosh- ka lakes in southern Egypt using re- mote sensing. Environ Earth Sci. 73 (7): 3623–3632. 19. Land Processes DAAC NASA (2015) USGS Earth Resources Observation and Science (EROS) Center, MODIS Reprojection Tool (MRT). Sioux Falls, South Dakota. Available at: https://mrtweb.cr.usgs.gov/. 20. NASA and Japan ASTER Program (2015) ASTER Global Digital Elevation Model (ASTER GDEM) data center. Version 2. Available at: http://gdem.ersdac.jspacesystems.or.j p/search.jsp. 21. Fox J, Weisberg S (2010) Nonlinear regression and nonlinear least squares in R: An appendix to an R companion to applied regression. Retrieved from McMaster University: http://socserv.mcmaster.ca/jfox/Book s/Companion/appendix/Appendix- Nonlinear-Regression.pdf. 22. Fox J, Weisberg S (2011) An {R} Com- panion to Applied Regression. Sec- ond ed. Sage, Thousand Oaks, CA. 23. Fagih-Nayini F, Mohebali M, Javadian E (2002) Epidemiology of Visceral leishmaniasis in Kordan-savojbalagh. Pajoohandeh J. 7(2): 9–15. 24. Mazloumi Gavgani AS, Hodjati MH, Mohite H, Davies CR (2002) Effect of insecticide-impregnated dog col- lars on incidence of zoonotic visceral leishmaniasis in Iranian children: a matchedcluster randomised trial. The Lancet. 360(9330): 374–379. 25. Soleimanzadeh G, Edrissian GH, Movah- hed-Danesh AM, Nadim A (1993) Epidemiological aspects of kala-azar in Meshkin-Shahr, Iran: human in- fection. Bull World Health Organ. 71(6): 759–762. 26. Mohebali M, Edrissian GH, Shirzadi MR, Hosseingholizadeh G, Pashaei MH, Ganji A, Zarei Z, Kousha A, Ak- houndi B, Hajjaran H, Malekafzali H (2010) Integrated visceral leishman- iasis surveillance system in primary care for children in Meshkin-Shahr district, north-western Islamic Repub- lic of Iran. Eastern Mediterranean health journal = La revue de sante de la Mediterranee orientale = al-Majal- lah al-sihhiyah li-sharq al-mutawas- sit. 16(10): 1050–1054. 27. Edrissian GH (1996) Visceral leishman- iasis in Iran and the role of serologi- cal tests in diagnosis and epidemi- ological studies. Parasitology for 21st Century ICOPA VIII, Izmir, Turkey: CAB International, pp. 63–78. 28. Mohebali M, Edrissian GH, Shirzadi MR, Hosseingholizadeh Y, Pashaei MH, Ganji A, Zarei Z, Kousha A, Ak- houndi B, Hajjaran H (2011) Estab- lishment and integration of a visceral leishmaniasis surveillance system in the primary health care system and its evaluation, Meshkin-Shahr District, Iran. J Sch Public Health Inst Public Health Res. 9(2): 1–8. 29. Mahami M, Moheb Ali M, Keshavarz H, Hajaran H, Akhoondi B, Zarei Z, Charedar S (2006) A seroepidemio- logical survey of visceral leishman- iasis (KALA-AZAR) in Germi Dis- trict, Ardabil Province. J Sch Public Health Inst Public Health Res. 4(1): 45–55. 30. Mohammadi-Kheyrabadi K, Mohebali M, Mamishi S, Arshi S (2004) Epidemi- ological characteristics of Kala-Azar http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 57 http://jad.tums.ac.ir Published Online: March 18, 2018 in hospitalized patients in Ardebil Prov- ince. J Sch Public Health Inst Public Health Res. 2(2): 11–24. 31. Nadim A, Tahvildar Birooni GH, Mohe- bali M (1998) A report of cutaneous and post kala-azar dermal leishmaniasis in Moghan, Ardabil Province Hakim Res J. 1(1): 57–60. 32. Arshi S, Mohebali M, Akhoundi B, Sadeghi-Bazargani H, Sepehram V, Zarei Z, Hajikhani S, Sezavar SH (2002) Identification of a new en- demic focus of kala-azar and seroepi- demiological study of visceral Leish- mania infection in Ardabil Province. J Sch Public Health Inst Public Health Res. 1: 9–18. 33. Mazloumi Gavgani AS, Malaki Rav- asan N, Mazloumi Gavgani F (2011) Comparison of Nomadic and non- Nomadic lifestyles in transmission of Visceral Leishaniasis. J Gorgan Univ Med Sci. 13(1): 94–100. 34. Mirsamadi N, Mohebali M, Atari MR, Edrisian GH (2003) Serological survey of Visceral leishmaniasis (Kala-Azar) in Azarshahr, Azarbaijan Province, northwest of Iran. Hakim Res J. 6(1): 17–22. 35. Yaghoobi-Ershadi MR, Shahbazi F, Dar- vishi M, Akhavan AA, Jafari R, Khajeian M, Rassi Y, Soleimani H, Shirzadi MR, Hanafi-Bojd AA, Da- rabi H, Arandian MH, Sanei-Dehkor- di A, Heidari M (2013) Molecular ep- idemiological study of cutaneous leish- maniasis in the focus of bushehr City, southwestern Iran. J Arthropod Borne Dis. 7(2): 113–121. 36. Hamzavi V, Mohebali M, Edrissian G, Fourozani A (2000) Epidemiological study of cutaneous leishmaniasis in Dashti and Dashtestani district of Bu- shehr Province. Iran J Public Health. 29: 177–191. 37. Mohebali M, Hamzavi Y, Edrissian GH, Forouzani A (2001) Seroepidemiologi- cal study of visceral leishmaniasis among humans and animal reservoirs in Bushehr Province, Islamic Repub- lic of Iran. Eastern Mediterranean health journal = La revue de sante de la Mediterranee orientale = al-Majal- lah al-sihhiyah li-sharq al-mutawassit. 7(6): 912–917. 38. Ghatee M, Sharifi I, Mirhendi H, Kanan- nejad Z, Hatam G (2013) Investigation of Double-Band Electrophoretic Pattern of ITS-rDNA Region in Iranian Iso- lates of Leishmania tropica. Iran J Par- asitol. 8(2): 264–272. 39. Razmjou S, Hejazy H, Motazedian MH, Baghaei M, Emamy M, Kalantary M (2009) A new focus of zoonotic cuta- neous leishmaniasis in Shiraz, Iran. Trans R Soc Trop Med Hyg. 103(7): 727–730. 40. Rassi Y, Javadian E, Jalali M, Motazedian MH, Vatndoost H (2004) Investigation on Zoonotic Cutaneous Leishmaniasis, Southern Iran. Iran J Public Health. 33(1): 31–35. 41. Fakoorziba MR, Baseri A, Eghbal F, Rezaee S, Azizi K, Moemenbellah- Fard MD (2011) Post-earthquake out- break of cutaneous leishmaniasis in a rural region of southern Iran. Ann Trop Med Parasitol. 105(3): 217–224. 42. Davami MH, Motazedian MH, Sarkari B (2010) The changing profile of cu- taneous leishmaniasis in a focus of the disease in Jahrom District, southern Iran. Ann Trop Med Parasitol. 104 (5): 377–382. 43. Sharafi M, Pezeshki B, Reisi A, Kalantari M, Naghizadeh MM, Dast Manesh S (2013) Detection of Cutaneous Leish- maniasis by PCR in Fasa District in 2012. J Fasa Univ Med Sci. 3(3): 266–270. 44. Noorpisheh S, Naghizadeh MM, Nikrouz L (2013) A study on the life quality http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 58 http://jad.tums.ac.ir Published Online: March 18, 2018 of patients suffering from leishmaniasis. J Fasa Univ Med Sci. 3(2): 155–162. 45. Entezari M, Eskandari F (2014) Rela- tionship between Climatic Factors and the Prevalence of Cutaneous Leishman- iasis in Larestan City J Mil Med. 16 (2): 99–104. 46. Fakhar M, Motazedian MH, Hatam GR, Asgari Q, Kalantari M, Mohebali M (2008) Asymptomatic human carriers of Leishmania infantum: possible res- ervoirs for Mediterranean visceral leish- maniasis in southern Iran. Ann Trop Med Parasitol. 102(7): 577–583. 47. Edrissian GH, Ahanchin AR, Gharachahi AM, Ghorbani M, Nadim A, Ardeha- li S, Hafizi A (1993) Seroepidemio- logical studies of visceral leishman- iasis and search for animal reservoirs in Fars province, southern Iran. Iran J Med Sci. 18(3–4): 99–105. 48. Sahabi Z, Seyedi Rashti MA, Nadim A, Javadian E, Kazemeini M, Abai MR (1992) A Preliminary Report on the Natural Leptomonad Infection of Phlebotomus Majorin an Endemic Focus of Visceral Leishmaniasis (Vl) in Fars Province, South of Iran. Iran J Public Health. 21(1–4): 87–93. 49. Fakhar M, Motazedian M, Asgari G, Mohebali M, Mehrbani D (2006) A new endemic focus of Visceral Leish- maniosis in southern Iran. Armaghane- Danesh. 11(2): 103–113. 50. Mollalo A, Alimohammadi A, Shirzadi MR, Malek MR (2015) Geographic information system-based analysis of the spatial and spatio-temporal distri- bution of zoonotic cutaneous leish- maniasis in Golestan Province, north- east of Iran. Zoonoses Public Health. 62(1): 18–28. 51. Rahbarian N, Mesgarian A, Mahmoudi Rad M, Hajaran H, Shahbazi F, Mesgarian Z, Taghipour N (2009) Identification of Leishmania species isolated from human cutaneous leish- maniasis using PCR method. J Res Health Sci. 9(2): 48–51. 52. Sofizadeh A, Cherabin M, Mehravaran A (2012) Cutaneous leishmaniasis in Gonbad Kavoos, North of Iran (2009– 11): an epidemiological study. J Gor- gan Univ Med Sci. 14(4): 100–106. 53. Baghaei A, Parvizi P, Amirkhani A, Ho- narvar MR, Badiei F (2012) Identifi- cation of Leishmania using micro- scopic and molecular methods in sus- pected patients of Cutaneous Leish- maniasis by targeting ITS-rDNA gene, Golestan Province, Iran (2009–10). J Gorgan Univ Med Sci. 14(3): 72–81. 54. Mesgarian F, Rahbarian N, Hajaran H, Shahbaz F, Mesgarian Z, Taghipour N (2010) Identification of Leishma- nia species isolated from human cu- taneous Leishmaniasis in Gonbad-e- Qabus city using a PCR method dur- ing 2006–2007. J Tehran Univ Med Sci. 68(4): 250–256. 55. Hanafi Bojd AA, Yaghoobi-Ershadi MR, Zamani GH, Barzekar A, Jafari R, Pour Abazari G (2006) Epidemiologic as- pect of Cutaneaus Leishmaniasis in Hajiabad, Hormozgan, Iran (2003). Medical J Hormozgan Univ Med sci. 10(1): 63–70. 56. Nazari M (2012) Cutaneous leishmaniasis in Hamadan, Iran (2004–2010). Zahedan J Res Med Sci. 13(9): 39–42. 57. Azizi K, Soltani A, Alipour H (2012) Molecular detection of Leishmania isolated from cutaneous leishmaniasis patients in Jask County, Hormozgan Province, Southern Iran, 2008. Asian Pac J Trop Med. 5(7): 514–517. 58. Asgari Nezhad H, Borhani M, Norouzi M, Merzaie M (2012) Cutaneous Leish- maniasis in school children in a bor- der area at south-west of Iran. Sci Par- asitol. 13(4): 153–158. http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 59 http://jad.tums.ac.ir Published Online: March 18, 2018 59. Shazad B, Abbaszadeh B, Khamesipour A (2005) Comparison of topical paro- momycin sulfate (twice/day) with in- tralesional meglumine antimoniate for the treatment of cutaneous leishman- iasis caused by L. major. Eur J Der- matol. 15(2): 85–87. 60. Bahrami AM, Delpisheh A, Zamandousti A (2011) Epidemiologic aspects of cutaneous Leishmaniasis and trapeu- tic effect of traditional ointment con- taining onion juice compared with Am- photricin B. J Kermanshah Univ Med Sci. 15(4): 251–259. 61. Yazdanpanah HA, Rostamianpur M (2013) Analysis of Spatial Distribu- tion of Leishmaniasis and its Relation- ship with Climatic Parameters (Case Study: Ilam Province). Bull Env Phar- macol Life Sci. 2(12): 80–86. 62. Roghani A, Yasemi M, Jalilian M, Abdi J, Rezayee K (2012) Epidemiology of cutaneous leishmaniasis in Ilam. Pa- joohesh dar Pezeshki. 36(5): 50–53. 63. Shiee MR, Mohebali M, Doroodgar A, Teimouri A, Afzali H, Shirzadi MR (2012) A molecular and parasitologi- cal survey on cutaneous leishmaniasis patients from historical city of Kashan in Isfahan Province, center of Iran. Asian Pac J Trop Dis. 2(6): 421–425. 64. Zahraei-Ramazani AR, Yaghoobi-Ershadi MR, Mokhtari AR, Akhavan AA, Abdoli H, Arandian MH (2007) An- throponotic Cutaneous Leishmaniasis in Nonendemic Quarters of a Centeral City in Iran. Iran J Public Health. 36 (2): 7–11. 65. Doroodgar A, Sayyah M, Doroodgar M, Mahbobi S, Nemetian M, Rafizadeh S, Rassi Y (2012) Progressive increas- ing of cutaneous leishmaniasis in Kashan district, central of Iran. Asian Pac J Trop Med. 2(4): 260–263. 66. Talari S, Rezvan Talaei R, Shajari G (2006) A clinical and epidemiology survey of cutaneous leishmaniasis in childhood. J Med Sci. 6(4): 693–697. 67. Yaghoobi-Ershadi M, Hanafi-Bojd A, Akhavan A, Zahrai-Ramazani A, Mohebali M (2001) Epidemiological study in a new focus of cutaneous leishmaniosis due to Leishmania ma- jor in Ardestan town, central Iran. Acta Trop. 79(2): 115–121. 68. Mohammadian A, Almasi M, Kheybari H, Karimi F, Tehrani S (1999) Iden- tification of species and characteris- tics of an outbreak of cutaneous leish- maniasis in a new focus of Iran. Bi- omedical J Iran. 3(1): 31–39. 69. Yaghoobi-Ershadi M, Javadian E (1995) leishmaniose cutanée zoonotique au nod d'Ispahan: le point sur l'infection humaine en 1991. Bull Soc Pathol Exot. 88(1): 42–45. 70. Yaghoobi-Ershadi M, Moosa-Kazemi S, Zahraei-Ramazani A, Jalai-Zand A, Akhavan A, Arandian M, Abdoli H, Houshmand B, Nadim A, Hosseini M (2006) Evaluation of deltamethrin-im- pregnated bed nets and curtains for con- trol of zoonotic cutaneous leishmaniasis in. Bull Soc Pathol Exot. 99(1): 43–48. 71. Nilforoushzadeh MA, Tahmores J, Ataei B, Jaafary F (2002) Evaluation of in- tegrated pest management method on controlling zoonotic Cutaneous Leish- maniasis in Emamzadeh Agha Ali Ab- bas, Natanz County, Isfahan Province, 1996–1998. Iran J Dermatol. 6(2): 78– 83. 72. Doroodgar A, Mahbobi S, Nemetian M, Sayyah M, Doroodgar M (2009) An epidemiological study of cutaneous leishmaniasis in Kashan (2007–2008). Koomesh. 10(3): 177–184. 73. Ahmadi NA, Ghafarzadeh M, Jalali Ga- loosang A, Gholami Parizad E (2013) An epidemiological study of Cutane- ous Leishmaniasis with emphasis on incidence rate in Kashan, Isfahan Prov- http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 60 http://jad.tums.ac.ir Published Online: March 18, 2018 ince. Scientific J Ilam Univ Med Sci. 21(2): 1–9 (In Persian). 74. Ebadi M, Hejazi SH (2003) The epide- miological study of cutaneous Leish- maniasis situation in the students of primary school in Isfahan Borkhar re- gion. J Kerman Univ Med Sci. 10(2): 92–98. 75. Yaghoobi-Ershadi M, Hanafi A, Akha- van A, Zahraei Ramazani A, Mohe- bali M (1999) Cutaneous leishmania- sis in Ardestan Town. Hakim Res J. 1(3): 206–215. 76. Doroudgar A, Dehghani R, Hooshyar H, Sayyah M (1996) Study of the prev- alence of the cutaneous Leishmania- sis in south-east part of Kashan. J Kerman Univ Med Sci. 3(2): 80–86. 77. Nadim A, Aflatoonian MR (1995) An- throponotic cutaneous leishmaniasis in the city of Bam, southeast Iran. Iran J Public Health. 24(1–2): 15–24. 78. Sharifi I, Nakhaei N, Aflatoonian M, Parizi MH, Fekri A, Safizadeh H, Shirzadi M, Gooya M, Khamesipour A, Nadim A (2011) Cutaneous leish- maniasis in bam: a comparative evalu- ation of pre- and post-earthquake years (1999–2008). Iran J Public Health. 40(2): 49–56. 79. Yaghoobi-Ershadi MR (1977) Cutene- ous leishmaniasis in Kerman. [MSPH theses]. School of Public Health, Tehran University, Iran. 80. Mirzaei M, Sharifi I, Poursmaelian S (2012) A new focus of anthroponotic cutaneous leishmaniasis and identi- fication of parasite species by nested PCR in Jiroft, Iran. Comp Clin Path. 21(5): 1071–1075. 81. Aflatoonian MR, Sharifi I, Hakimi Parizi M, Fekri AR, Aflatoonian B, Sharifi M, Khosravi A, Khamesipour A, Shar- ifi H (2014) A prospective cohort study of cutaneous leishmaniasis risk and opium addiction in south eastern Iran. PLoS One. 9(2): e89043. 82. Sharifi I, Fekri AR, Aflatonian MR, Nadim A, Nikian Y, Kamesipour A (1998) Cutaneous leishmaniasis in primary school children in the south-eastern Ira- nian city of Bam, 1994–95. Bull World Health Organ. 76(3): 289–293. 83. Sharifi I, Poursmaelian S, Aflatoonian MR, Ardakani RF, Mirzaei M, Fekri AR, Khamesipour A, Parizi MH, Harandi MF (2011) Emergence of a new focus of anthroponotic cutane- ous leishmaniasis due to Leishmania tropica in rural communities of Bam District after the earthquake, Iran. Trop Med Int Health. 16(4): 510–513. 84. Aflatoonian MR, Sharifi I, Poursmaelian S, Hakimi-Parizi M, Ziaali N (2013) The Emergence of Anthroponotic Cu- taneous Leishmaniasis Following the Earthquake in Southern Villages of Bam District, southeastern Iran, 2010. J Arthropod Borne Dis. 7(1): 8–14. 85. Nadim A, Motabar M, Houshmand B, Keyghobadi K, Aflatonian MR (1995) Evaluation of pyrethroid impregnated bednets for control of anthroponotic cutaneous leishmaniasis in Bam (Islam- ic Republic of Iran), WHO/LEISH/95. 37: Geneva. p. 21. 86. Sharifi I, Aflatoonian MR, Aflatoonian B, Kermanizadeh A (2015) The severity of cutaneous leishmaniasis before and after the earthquake in Bam, southeast- ern Iran. J Parasit Dis. 39(4): 741–744. 87. Aflatoonian MR, Sharifi I (2011) The Epidemiology of cutaneous leishman- iasis in the city and suburb of Bam in 2010: Active case-finding, treatment and health education of the school children. Iran J Epidemiol. 7(3): 52–57. 88. Aflatoonian MR, Sharifi I (2010) Preva- lence rate of cutaneous leishmaniasis in Bam District during 20 Years (1988–2007). J Kerman Univ Med Sci. 17(4): 297–306. http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 61 http://jad.tums.ac.ir Published Online: March 18, 2018 89. Mirzadeh A, Hejazi Zadeh B, Mesgar pour B, Golezar A, Halakooiiy Nayini K (2008) The process of development of high-risk centers of cutaneous leish- maniasis in the city of Kerman in 2002–2005 years andevaluation of efec- tive environmental factors by using GIS. Iran J Epidemiol. 4(3 and 4): 5–17. 90. Aflatoonian MR, Sharifi I, Aflatoonian B (2014) The Effects of Prevention Knowledge Compared to the Improve- ment of Housing Conditions in Con- trol of Anthroponotic Cutaneous Leish- maniasis: a Natural Experience of Earth- quake in Bam/Iran. J Kerman Univ Med Sci. 21(3): 247–258. 91. Aflatoonian MR, Sharifi I (2007) Prev- alence of Cutaneous Leishmaniasis in School Children in Bam and Barawat. Iran in 2006. J Kerman Univ Med Sci. 14(2): 82–89. 92. Akhavan AA, Yaghoobi-Ershadi MR, Ha- sibi F, Jafari R, Abdoli H, Arandian MH, Soleimani H, Zahraei-Ramazani AR, Mohebali M, Hajjaran H (2007) Emergence of cutaneous leishmania- sis due to Leishmania major in a new focus of southern Iran. J Arthropod Borne Dis. 1(1): 1–8. 93. Aghaei-Afshar A, Vatandoost H, Sharifi I, Rassi Y, Abai MR, Oshaghi MA, Yaghoobi-Ershadi MR, Rafizadeh S (2013) First determination of impact and outcome indicators following in- door residual spraying (IRS) with del- tamethrin in a new focus of anthropon- otic cutaneous leishmaniasis (ACL) in Iran. Asian Pac J Trop Med. 3(1): 5–9. 94. Khosravi A, Sharifi I, Dortaj E, Aghaei Afshar A, Mostafavi M (2013) The present status of cutaneous leishmania- sis in a recently emerged focus in South-west of Kerman Province, Iran. Iran J Public Health. 42(2): 182–187. 95. Sharifi I, Zemani F, Aflatoonian MR, Fekri AR (2008) A report of cutaneous leish- maniasis epidemic and its probable causative factors in Baft District, Ker- man Province. Iran J Epidemiol. 4(1): 53–58. 96. Mahmoudvand H, Mohebali M, Sharifi I, Keshavarz H, Hajjaran H, Ak- houndi B, Jahanbakhsh S, Zarean M, Javadi A (2011) Epidemiological as- pects of visceral leishmaniasis in Baft District, Kerman Province, Southeast of Iran. Iran J Parasitol. 6(1): 1–11. 97. Barati M, Parizi MHD, Sharifi I (2008) Epidemiological and Clinical aspects of kala-azar in hospitalized children of Kerman Province, during 1991–2006. J Kerman Univ Med Sci. 15(2): 148– 155. 98. Hamzavi Y, Khademi N (2015) Trend of cutaneous leishmaniasis in Kermanshah Province, west of iran from 1990 to 2012. Iran J Parasitol. 10(1): 78–86. 99. Nazari N, Faraji R, Vejdani M, Mekaeili A, Hamzavi Y (2012) Prevalence of Cutaneous leishmaniases in Pations referred to Kermanshah hygienic cen- ters (2006–2008) Zahedan. J Res Med Sci. 14(8): 77–79. 100. Hamzavi Y, Khademi N (2013) The Analytical study of Cutaneous Leish- maniasis in Kermanshah (2011–2012). J Kermanshah Univ Med Sci. 17(9): 582–589. 101. Hamzavi Y, Hamzeh B, Mohebali M, Akhoundi B, Ajhang K, Khademi N, Ghadiri K, Bashiri H, Pajhouhan M (2012) Human visceral leishmaniasis in Kermanshah Province, Western Iran, during 2011–2012. Iran J Parasitol. 7 (4): 49–56. 102. Alavinia S, Arzamani K, Reihani M, Jafari J (2009) Some epidemiological aspects of cutaneous leishmaniasis in Northern Khorasan Province, Iran. J Arthropod Borne Dis. 3(2): 50–54. 103. Torabi V, Mohebali M, Edrissian GH, Keshavarz H, Mohajery M, Hajaran http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 62 http://jad.tums.ac.ir Published Online: March 18, 2018 H, Akhoondi B, Sanaati A, Zarei Z, Delshad A (2007) Seroepidemiology of visceral leishmaniasis by DAT in Bojnord-North Khorasan. Iran J Ep- idemiol. 4(3 and 4): 43–50. 104. Sattar Pagheh A, Fakhar M, Sharif M, Danesh V, Ahmadi Z (2013) Epide- miological Survey of Cutaneous Leish- maniasis due to Leishmania tropica in a New Focus in Khorasan Razavi Prov- ince. J Mazandaran Univ Med Sci. 23(103): 47–53. 105. Moosa-Kazemi SH, Yaghoobi-Ershadir MR, Akhavan AA, Abdoli H, Zahraei- Ramazani AR, Jafari R, Houshmand B, Nadim A, Hosseini M (2007) Deltame- thrin-impregnated bed nets and curtains in an anthroponotic cutaneous leish- maniasis control program in northeast- ern Iran. Ann Saudi Med. 27(1): 6–12. 106. Saadabadi F, Mohajery M, Poostchi E, Shamsian SAA (2013) Identification of Leishmania species causing cuta- neous leishmaniasis using Random Am- plified Polymorphic DNA (RAPD- PCR) in Kharve, Iran. Reports of Bio- chemistry and Molecular Biology. 1 (2): 69–73. 107. Hassanpour K, Aghamollaei H, Golpich M, Amani J, Taheri A, Farnoosh G (2014) Molecular epidemiological study of cutaneous leishmaniasis in the east north of Iran. Asian Pac J Trop Dis. 4: 540–544. 108. Hoseini Farash BR, Mohajery M, Shamsian AA, Rezaie A (2011) Iden- tification of Leishmania species caus- ing cutaneous Leishmaniasis by PCR method in Torghabeh-Shandiz during 2009–2010. Forth Congress of Labor- atory and Clinic, 2011 December 21– 23, Tehran, Iran. 109. Khajedaluee M, Yazdanpanah M, Seyed Nozadi SM, Fata A, Juya MR, Ma- soudi MH, Najaf Najafi M (2014) Epidemiology of cutaneous leishman- iasis in Razavi Khorasan in 2011. Med- ical J Mashhad Univ Med sci. 57(4): 647–654. 110. Mohajery M, Hajaran H, Shamsian A, Tavakol Afshari J, Saadabadi F (2008) Identification of cutaneous leishman- iasis species by RAPD-PCR in Ney- shabur. Medical J Mashhad Univ Med sci. 51(100): 79–86. 111. Yaghoobi-Ershadi MR, Akhavan AA, Zahraei-Ramazani AV, Abai MR, Ebrahimi B, Vafaei-Nezhad R, Hanafi- Bojd AA, Jafari R (2003) Epidemio- logical study in a new focus of cuta- neous leishmaniasis in the Islamic Re- public of Iran. Eastern Mediterranean health journal = La revue de sante de la Mediterranee orientale = al-Majallah al-sihhiyah li-sharq al-mutawassit. 9 (4): 816–826. 112. Karimi Zarchi AA, Mahmodzadeh A, Vatani H, Shirbazo S (2004) Epidemi- ological study of Cutaneous Leishman- iasis in Sarakhs District. J Shahid Se- doghi Univ Med Sci. 12(1): 30–35. 113. Sahabi Z, Mohajery M, Seyedi Rashti M (1999) Fauna of arthropods vec- tors and their medical importance in the free trade zone Sarakhs ferns, sand flies and leishmaniasis. Hakim Res J. 2(3): 177–182. 114. Akbari E, Mayvaneh E, Entezari A, Nazari M (2014) Survey of the Role of Bioclimatic Factors in the Out- break of Cutaneous Leishmaniasis. Iran J Epidemiol. 10(3): 65–74. 115. Karamian M, Bojd F, Sedigh M, Hem- mati M, Saadatjoo A, Barati DA (2013) Molecular identification of cu- taneous leishmaniasis agents in Bir- jand, Iran. J Birjand Univ Med Sci. 20(2): 183–190. 116. Kassiri H, Shemshad K, Kassiri A, Sho- jaee S, Sharifinia N, Shemshad M (2012) Clinical laboratory and epide- miological research on cutaneous leish- http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 63 http://jad.tums.ac.ir Published Online: March 18, 2018 maniasis in the south west of Iran. Ar- chives of Clinical Infectious Diseas- es. 7(4): 128–131. 117. Kassiri H, Shojaei S, Kolahkaji K (2012) Epidemiological status of cutaneous leishmaniasis in Behbahan County, Khuzestan Province, Southwestern Iran. Int J Infect Dis. 16(Supplement 1): e348. 118. Kassiri H, Lotfi M, Farajifard P, Kassiri E (2014) Laboratory diagnosis, clinical manifestations, epidemiological situa- tion and public health importance of cutaneous leishmaniasis in Shushtar County, Southwestern Iran. J Acute Disease. 3(2): 93–98. 119. Kassiri H, Feizhaddad MH, Abdehpanah M (2014) Morbidity, surveillance and epidemiology of scorpion sting, cuta- neous leishmaniasis and pediculosis capitis in Bandar-Mahshahr County, Southwestern Iran. J Acute Disease. 3(3): 194–200. 120. Vazirianzadeh B, Hoseini S, Pour Rezaee S, Gardani H, Amraee K (2014) Prevalence of Cutaneous Leishman- iasis in Ramshir, Iran, An Epidemi- ological Study. International Archives of Health Sciences. 1(1): 41–37. 121. Kassiri H, Shemshad K, Lotfi M, Shemshad M (2013) Relationship Trend Analysis of of Cutaneous Leish- maniasis Prevalence and Climatolog- ical Variables in Shush County, South- West of Iran (2003–2007). Acad J En- tomol. 6(2): 79–84. 122. Spotin A, Rouhani S, Parvizi P (2014) The associations of Leishmania major and Leishmania tropica aspects by focusing their morphological and mo- lecular features on clinical appearances in Khuzestan Province, Iran. Bio Med Res Int. 2014: 913510. 123. Ghasemian M, Maraghi S, Samarbafza- deh AR, Jelowdar A, Kalantari M (2011) The PCR-based detection and identification of the parasites causing human cutaneous leishmaniasis in the Iranian city of Ahvaz. Ann Trop Med Parasitol. 105(3): 209–215. 124. Kassiri H, Shemshad K, Shojaee S (2013) The incidence rate of cutaneous leish- maniasis in Behbahan County of Khuzestan Province, Southwest of Iran. Jundishapur J Microbiol. 6(5): e7045. 125. Nejati J, Mojadam M, Bojd H, Ali A, Keyhani A, Habibi Nodeh F (2014) An epidemiological study of Cuta- neous Leishmaniasis in Andimeshk (2005–2010). Scientific J Ilam Univ Med Sci. 21(7): 94–101. 126. Behbahani A, Ahmadi S, Latifi SM, Sadeghi M (2012) Study of the fre- quency of Cutaneous Leishmaniasis in Omidieh District, Khouzestan Prov- ince, Southwest of Iran (2008–2010). Jundishapur J Health Sci. 4(4): 37–46. 127. Kassiri H, Mortazavi SH, Kazemi S (2011) The epidemiological study of cutaneous leishmaniasis in Khorram- shahr City, Khuzestan Province, South- west of Iran. Jundishapur J Health Sci. 3(2): 11–20. 128. Sarkari B, Pedram N, Mohebali M, Moshfe AA, Zargar MA, Akhoundi B, Shirzadi MR (2010) Seroepidemi- ological study of visceral leishman- iasis in Booyerahmad district, south- west Islamic Republic of Iran. Eastern Mediterranean health journal = La re- vue de sante de la Mediterranee ori- entale = al-Majallah al-sihhiyah li-sharq al-mutawassit. 16(11): 1133–1136. 129. Sarkari B, Moshfe AA, Pedram N, Zargar MA, Yazdanpanah B, Akhoundi B, Hasani Z, Mohebali M (2007) Serologi- cal Study of Visceral Leishmaniasis in Boyer Ahmad Township in 2005. Armaghane-Danesh. 12(2): 69–77. 130. Kheirandish F, Chegeni Sharafi A, Kazemi B, Mohebali M, Sarlak A, Tarahi MJ, Holakouee K, Hajaran H (2013) Iden- http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 64 http://jad.tums.ac.ir Published Online: March 18, 2018 tification of Leishmania species using PCR assay on giemsa-stained slides prepared from cutaneous leishmaniasis patients. Iran J Parasitol. 8(3): 382–388. 131. Amraee K, Rastegar HA, Beiranvand E (2013) An epidemiological study of cutaneous leishmaniasis in Poledokhtar District, Lorestan Province, Southwest- ern of Iran, 2001–2011. Jundishapur J Health Sci. 5(1): 55–62. 132. Ahmadi NA, Modiri M, Mamdohi S (2013) First survey of cutaneous leish- maniasis in Borujerd county, western Islamic Republic of Iran. Eastern Med- iterranean Health J = La revue de sante de la Mediterranee orientale = al-Majal- lah al-sihhiyah li-sharq al-mutawas- sit. 19(10): 847–853. 133. Chegeni-Sharafi A, Amani H, Kayedi MH, Yarahahmadi A, Saki M, Mehrdad M, Nasiri E (2011) Epidemiological survey of cutaneous leishmaniasis in Lorestan Province (Iran) and introduction of dis- ease transmission in new local areas. J Ilam Univ Med Sci. 19(1):54–60. [In Persian]. 134. Chegeni-Sharafi A, Ourmazdi H, Mohe- bali M, Akhlaghi L, Sharafi MM, Akhondi B (2005) Seroepidemiolog- ical study of Visceral Leishmaniasis (human infection) in east Miankooh area in Lorestan Province by Direct Agglutination Test (DAT). Yafte J Med Sci. 7(3 and 4): 31–35. 135. Rostami MN, Saghafipour A, Vesali E (2013) A newly emerged cutaneous leishmaniasis focus in central Iran. Int J Infect Dis. 17(12): e1198–e1206. 136. Akhavan AA, Yaghoobi-Ershadi MR, Mehdipour D, Abdoli H, Farzinnia B, Mohebali M, Hajjaran H (2003) Ep- idemic outbreak of cutaneous leish- maniasis due to leishmania major in Ghanavat Rural District, Qom Prov- ince, Central Iran. Iran J Public Health. 32(4): 35–41. 137. Saghafipour A, Rassi Y, Abai MR, Oshaghi MA, Farzinnia B, Mostafavi R, Karimian F (2013) Outbreak of Zoonotic Cutaneous Leishmaniasis: A Report. Arch Hyg Sci. 2(2): 45–54. 138. Saghafipour A, Amir A, Rassi Y, Mo- stafavi R (2012) Epidemiology of cu- taneous leishmaniasis in Qom. J Qom Univ Med Sci. 6(1): 83–88. 139. Mohammadi Azni S, Nokandeh Z, Khorsandi AA, Sanei Dehkordi AR (2010) Epidemiology of cutaneous leishmaniasis in Damghan. Iran J Mil Med. 12(3): 131–135. 140. Mohammadi Azni s, Rassi Y, Oshaghi MA, Yaghoobi-Ershdi M, Mohebali M, Abai MR, Mohtarami F, Nokandeh Z, Rafizadeh S, Khojami GM (2011) Determination of parasite species of cutaneous leishmaniasis using Nested PCR in Damghan-Iran, during 2008. J Gorgan Univ Med Sci. 13(1): 59–65. 141. Mohammadi Azni S, Nokandeh Z, Sanei Dehkordi AR (2010) Control of rural cutaneous leishmaniasis in Damghan in 2005–2006. Iran J Infect Dis. 15 (48): 29–32. 142. Rafati N, Shapouri Moghadam M, Ghor- bani R (2007) Epidemiological study of cutaneous Leishmaniasis in Damghan (2000–2006). Scientific J Semnan Univ Med Sci. 8(4): 247–253. 143. Fazaeli A, Fouladi B, Hashemi-Shahri SM, Sharifi I (2008) Clinical features of cu- taneous leishmaniasis and direct PCR- based identification of parasite spe- cies in a new focus in southeast of Iran. Iran J Public Health. 37(3): 44–51. 144. Fazaeli A, Fouladi B, Sharifi I (2009) Emergence of cutaneous leishmaniasis in a border area at south-east of Iran: an epidemiological survey. J Vector Borne Dis. 46(1): 36–42. 145. Fouladi B, Sharifi I, Hashemi SS, Morad GH, Sarabandi NA, Ebrahimzadeh A, Fazaeli A (2007) Evaluation of a http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 65 http://jad.tums.ac.ir Published Online: March 18, 2018 direct PCR in comparison with routine microscopy and in vitro culture for diagnosis of cutaneous leishmaniasis. Zahedan J Res Med Sci. 9(3): 9–15. 146. Salehzadeh A, Seyedi Rashti M (1996) Cutaneous leishmaniasis in North of Tehran. Scientific J Hamadan Univ Med Sci. 4(1): 34–39. 147. Abadi M (2004) Survey the situation of anthroponotic Cutaneous leishmaniasis in North-West of Tehran. [PhD dis- sertation]. Tehran Islamic Azad uni- versity, Iran. 148. Yaghoobi-Ershadi MR, Hanafi-Bojd AA, Javadian E, Jafari R, Zahraei-Ramazani AR, Mohebali M (2002) A new focus of cutaneous leishmaniasis caused by Leishmania tropica. Saudi Med J. 23 (3): 291–294. 149. Yaghoobi-Ershadi MR, Jafari R, Hanafi- Bojd AA (2004) A new epidemic focus of zoonotic cutaneous leishmaniasis in central Iran. Ann Saudi Med. 24(2): 98–101. 150. Jaafary R, Mohebali M, Dehghan Dehnavi AR, Soleimani H, Akhavan AA, Haja- ran H, Dehghan shadkam A, Fatahi J (2007) Epidemiological status of cu- taneous leishmaniasis in Bafgh City, Yazd Province 2005. J Shahid Sa- doghi Univ Med Sci. 15(2): 76–83. 151. Mozafari GH, Bakhshi Zadeh F (2011) Analysis the Roll of bioclimatic Fac- tors in Prevalence of Cutaneous Leish- maniasis in Yazd-Ardakan plain Geog- raphy and Development J. 23: 185–202. 152. Hanafi Bojd AA, Dehghani Tafti AA, Jafari R, Ehrampoosh MH (2003) Ef- ficacy of control program of zoonotic cutaneous leishmaniasis in a new focus of the disease. Medical J Mashhad Univ Med sci. 46(80): 19–25. 153. Yaghoobi-Ershadi MR, Marvi Moghadam N, Jaafary R, Akhavan A, Soleimani H, Zahraei Ramazani AR, Arandian MH, Dehghan Dehnavi AR (2008) Study of certain epidemiological as- pects of cutaneous Leishmaniasis in Khatam County, Yazd Province, Iran. J Shahid Sedoghi Univ Med Sci. 15 (4): 47–52. 154. Werneck G, Costa C, Walker A, David J, Wand M, Maguire J (2007) Multi- level modelling of the incidence of visceral leishmaniasis in Teresina, Bra- zil. Epidemiol Infect. 135(2): 195–201. 155. Signorini M, Cassini R, Drigo M, di Re- galbono AF, Pietrobelli M, Montarsi F, Stensgaard AS (2014) Ecological niche model of Phlebotomus pernicio- sus, the main vector of canine leish- maniasis in north-eastern Italy. Geo- spat Health. 9(1): 193–201. 156. Mokhtari M, Miri M, Nikoonahad A, Jalilian A, Naserifar R, Ghaffari HR, Kazembeigi F (2016) Cutaneous leish- maniasis prevalence and morbidity based on environmental factors in Ilam, Iran: Spatial analysis and land use regression models. Acta Trop. 163: 90–97. 157. Elnaiem DEA, Schorscher J, Bendall A, Obsomer V, Osman ME, Mekkawi AM, Connor SJ, Ashford RW, Thom- son MC (2003) Risk mapping of vis- ceral leishmaniasis: the role of local variation in rainfall and altitude on the presence and incidence of Kala-Azar in eastern Sudan. Am J Trop Med Hyg. 68(1): 10–17. 158. Karimi A, Hanafi-Bojd AA, Yaghoobi- Ershadi MR, Akhavan AA, Ghezel- bash Z (2014) Spatial and temporal distributions of phlebotomine sand flies (Diptera: Psychodidae), vectors of leishmaniasis, in Iran. Acta Trop. 132: 131–139. 159. Hanafi-Bojd AA, Yaghoobi-Ershadi MR, Haghdoost AA, Akhavan AA, Rassi Y, Karimi A, Charrahy Z (2015) Modeling the distribution of cutaneous leishman- iasis vectors (Psychodidae: Phlebotom- http://jad.tums.ac.ir/ J Arthropod-Borne Dis, March 2018, 12(1): 41–66 A Golpayegani et al.: Modeling of … 66 http://jad.tums.ac.ir Published Online: March 18, 2018 inae) in Iran: a potential transmission in disease prone areas. J Med Ento- mol. 52(4): 557–565. 160. Gholamrezaei M, Mohebali M, Hanafi- Bojd AA, Sedaghat MM, Shirzadi MR (2016) Ecological Niche Modeling of main reservoir hosts of zoonotic cu- taneous leishmaniasis in Iran. Acta Trop. 160: 44–52. 161. Oshaghi M, Ravasan NM, Javadian E, Rassi Y, Sadraei J, Enayati A, Vatan- doost H, Zare Z, Emami S (2009) Application of predictive degree day model for field development of sand- fly vectors of visceral leishmaniasis in northwest of Iran. J Vector Borne Dis. 46(4): 247–255. 162. Shirzadi MR, Mollalo A, Yaghoobi- Ershadi MR (2015) Dynamic relations between incidence of zoonotic cutane- ous leishmaniasis and climatic factors in Golestan Province, Iran. J Arthro- pod Borne Dis. 9(2): 148–160. 163. Sudhakar S, Srinivas T, Palit A, Kar SK, Battacharya SK (2006) Mapping of risk prone areas of Kala-Azar (Visceral leishmaniasis) in parts of Bihar State, India: an RS and GIS approach. J Vec- tor Borne Dis. 43(3): 115–122. 164. Menezes JA, Ferreira EdC, Andrade-Filho JD, Sousa AMd, Morais MHG, Rocha AMS, Machado-Coelho GLL, Lima FP, Madureira AP, Garcia TC (2015) An In- tegrated Approach Using Spatial Anal- ysis to Study the Risk Factors for Leish- maniasis in Area of Recent Transmis- sion. Bio Med Res Int. 2015: 621854. http://jad.tums.ac.ir/