Leishmaniasis is a public health problem worldwide. We aimed to predict ecological niche models (ENMs) for visceral (VL) and cutaneous (CL) leishmaniasis and the sand flies involved in the transmission of leishmaniasis in São Paulo, Brazil. Phlebotomine sand flies were collected between 1985 and 2015. ENMs were created for each sand fly species using Maximum Entropy Species Distribution Modeling software, and 20 climatic variables were determined. Nyssomyia intermedia (Lutz & Neiva, 1912) and Lutzomyia longipalpis (Lutz & Neiva, 1912), the primary vectors involved in CL and VL, displayed the highest suitability across the various regions, climates, and topographies. L. longipalpis was found in the border of Paraná an area currently free of VL. The variables with the greatest impact were temperature seasonality, precipitation, and altitude. Co-presence of multiple sand fly species was observed in the cuestas and coastal areas along the border of Paraná and in the western basalt areas along the border of Mato Grosso do Sul. Human CL and VL were found in 475 of 546 (86.7%) and 106 of 645 (16.4%) of municipalities, respectively. Niche overlap between N. intermedia and L. longipalpis was found with 9208 human cases of CL and 2952 cases of VL. ENMs demonstrated that each phlebotomine sand fly species has a unique geographic distribution pattern, and the occurrence of the primary vectors of CL and VL overlapped. These data can be used by public authorities to monitor the dispersion and expansion of CL and VL vectors in São Paulo state.
Keywords: Ecological niche model; Lutzomya longipalpis; Nyssomyia intermedia; co-presence.