The effect of climate variables on the incidence of cutaneous leishmaniasis in Isfahan, Central Iran

In recent years, there have been considerable changes in the distribution of diseases that are potentially tied to ongoing climate variability. The aim of this study was to investigate the association between the incidence of cutaneous leishmaniasis (CL) and climatic factors in an Iranian city (Isfahan), which had the highest incidence of CL in the country. CL incidence and meteorological data were acquired from April 2010 to March 2017 (108 months) for Isfahan City. Univariate and multivariate seasonal autoregressive integrated moving average (SARIMA), generalized additive models (GAM), and generalized additive mixed models (GAMM) were used to identify the association between CL cases and meteorological variables, and forecast CL incidence. AIC, BIC, and residual tests were used to test the goodness of fit of SARIMA models; and R-2 was used for GAM/GAMM. 6798 CL cases were recorded during this time. The incidence had a seasonal pattern and the highest number of cases was recorded from August to October. In univariate SARIMA, (1,0,1) (0,1,1)(12) was the best fit for predicting CL incidence (AIC=8.09, BIC=8.32). Time series regression (1,0,1) (0,1,1)(12) showed that monthly mean humidity after 4-month lag was inversely related to CL incidence (AIC=8.53, BIC=8.66). GAMM results showed that average temperature with 2-month lag, average relative humidity with 3-month lag, monthly cumulative rainfall with 1-month lag, and monthly sunshine hours with 1-month lag were related to CL incidence (R-2=0.94). The impact of meteorological variables on the incidence of CL is not linear and GAM models that include non-linear structures are a better fit for prediction. In Isfahan, Iran, meteorological variables can greatly predict the incidence of CL, and these variables can be used for predicting outbreaks.