Predictive modeling of sand fly distribution incriminated in the transmission of Leishmania (Viannia) braziliensis and the incidence of Cutaneous Leishmaniasis in the state of Paraná, Brazil.

Southern Brazil concentrates a considerable number of cases of cutaneous leishmaniasis reported since 1980, and Paraná is the state that most records CL cases in the region. The main sand fly species incriminated as vectors of Leishmania (Viannia) braziliensis (Vianna,1911) are Migonemyia (Migonemyia) migonei (França, 1920), Nyssomyia (Nyssomyia) neivai (Pinto, 1926) and Nyssomyia (Nyssomyia) whitmani (Antunes & Coutinho, 1936). In this study, we evaluated areas with climatic suitability for the distribution of these vectors and correlated these data with CL incidence in the state. The occurrence points of Mg. migonei, Ny. neivai, and Ny. whitmani were extracted from a literature review and field data. For CL analysis in the state of Paraná, data were obtained from the Informatics Department of the Unified Health System of Brazil (DATASUS), covering the period from 2001 to 2019. The layers of bioclimatic variables from the WorldClim database were used in the study. Species distribution modeling was developed using the MaxEnt Software version 3.4.4. ArcGIS software version 10.5 was used to develop suitability maps and the graphical representation of disease incidence. The AUC values were acceptable for all models (> 0,8). Bioclimatic variables BIO13 and BIO14 were the most influential in the distribution of Mg. migonei, while BIO19 and BIO6 were the variables that most influenced the distribution of Ny. neivai, and Ny. whitmani was most influenced by variables BIO5 and BIO9. During 19 years, 4992 cases of CL were reported in the state by 286 municipalities (71,6%). Northern Paraná showed the highest number of areas with very high and high climatic suitability for the occurrence of these species, coinciding with the highest number of CL cases. The modeling tools allowed analyzing the association between climatic variables and the geographical distribution of CL in the state. Moreover, they provided a better understanding of the climatic conditions related to the distribution of different species, favoring the monitoring of risk areas, the implementation of preventive measures, risk awareness, early and accurate diagnosis, and consequent timely treatment.