Spatio-temporal distribution of human American visceral leishmaniasis in the Western region of Sao Paulo State, from 2004 to 2018

ABSTRACT

Leishmaniasis comprises a group of zoonotic diseases caused by protozoa belonging to the Leishmania genus, noting that the visceral form is the most severe and lethal, if untreated. Nowadays visceral leishmaniasis is widespread in Brazil and the Adamantina microregion, located in the west of Sao Paulo State, has been affected by Human American Visceral Leishmaniasis (HAVL) since 2004. We evaluated the epidemiological profile of HAVL in the Adamantina microregion through a Geographic Information Systems (GIS) and established its incidence rate by location and time. Notified cases were provided by the Sao Paulo State Epidemiological Surveillance Center. Home addresses of patients who tested positive to HAVL were converted into geographic coordinates through the Google Geocoding Application Programming Interface submitted to ArcMap 10.5 System for georeferencing. Kernel spatial analyses were performed to obtain the incidence distribution and the total area involvement rate. From 2004 to 2018, 325 cases of HAVL were diagnosed in 11 of the 12 municipalities belonging to the of Adamantina microregion. The disease has disseminated to the Northwest and East-Southeast directions, taking place along the Comandante Joao Ribeiro de Barros highway, with higher incidences rates in the municipalities where the highway passes. HAVL incidence was higher in children aged between 0 to 9 years and in the elderly; there was no difference in relation to sex and the majority of cases were located in urban areas. The determination of the epidemiological profile and the the spread of disease patterns can indicate possible areas of vulnerability, in order to contribute to the management and prevention of the disease through a strategic resources optimization.

Visceral leishmaniasis; Kernel density; Incidence; Risk factors; Epidemiology; Georeferencing

INTRODUCTION

Leishmaniasis, a group of infectious diseases caused by protozoa of the genus Leishmania, can affect humans, domestic and wild animals around the world ¹ . These protozoa are transmitted by small insects belonging to genus Phlebotomus (in the Old World) and Lutzomyia (in the New World - Americas). Clinical manifestation of human leishmaniasis are divided into visceral and tegumentary forms depending on the infectious Leishmania species, specific geographical distributions, and the host’s immunological status ¹ .

Visceral leishmaniasis is the most severe form of the disease, being lethal in 95% of untreated cases ² . This disease is part of the group of neglected infectious diseases worldwide, as they mainly affect low socioeconomic populations with restricted access to health services, therefore being more vulnerable ³ .

The first reference that visceral leishmaniasis would be endemic in Brazil was made by Penna in 1934, who found amastigotes in biopsies of 41 fatal cases with suspected Yellow Fever, living in eight States of the Northeastern, North and Midwest regions ⁴ . Subsequently, the Northeast region became endemic for human American visceral leishmaniasis (HAVL) in the country, and the first epidemic dates back to 1950 in Ceara State, presenting with persistent high lethality that has been recently verified, and in 1980 the disease gained expressiveness with the dispersion to other influential States of this region ⁴ - ⁷ . In this period, the Northeast region accounted for 93% of the total cases in Brazil, which were concentrated in rural areas in line with the wild behavior of the transmission vectors ⁶ .

The proximity of rural and urban areas associated with the adaptation of phlebotomines species to urban settings resulted in cases of HAVL in the city of Terezina located in Piaui State ⁶ . In Pernambuco State, the essentially rural disease paradigm shift was observed in the 1990s, with consequent expansion of the disease to new areas ⁸ , ⁹ . The dissemination of phlebotomine vectors, especially Lutzomyia longipalps and the migration of the population along with domestic dogs infected by Leishmania infantum resulted in the dissemination of HAVL to others regions of the North, Midwest and Southeast of Brazil ¹⁰ , ¹¹ , decreasing the representativeness of the number of cases in the Northeast that reported 43.1% of them in the country, in 2012 ¹² . Currently, the HAVL occurs in 22 Brazilian States distributed among the five major Brazilian regions, presenting with social, geographical and climate particularities ¹³ .

In the West of Sao Paulo State, the first cases of HAVL were reported in urban areas in 1999, then quickly disseminating to neighboring municipalities following the spread and the transmission potential of the vector Lutzomyia longipalpis sandfly, found in this region of the State ¹⁴ , ¹⁵ .

Besides the vector density and its dissemination capacity, another major factor that may increase the complexity of transmission in this region is the proximity between municipalities ¹⁶ . Furthermore, changes in the economical activity leading to an increasing flow of goods and people through the highways in Western municipalities of the State are also suggested as causes for the vector and the HAVL increments ¹⁵ .

The city of Adamantina is located in a microregion affected by HAVL since 2004, and it was the source of dissemination of the disease ¹⁷ , due to its proximity to a secondary highway. In this way, building an epidemiological profile based on spatial and temporal analyses, with determination of incidence rates, in addition to associated environmental and socioeconomical characteristics, is a crucial mechanism in the development of disease control strategies ¹⁸ . In this sense, this work aimed to conceive a spatial and temporal distribution of HAVL in the municipalities of the Adamantina microregion, since its onset in 2004, until 2018, through georeferencing and then to verify the pattern of disease distribution, as well as the possible areas of vulnerability in order to contributing to the management and prevention of the disease through strategic resources optimization.

Figure 1
Study area showing the Adamantina microregion, its geographic location in Sao Paulo State and in Brazil; highway Comandante Joao Ribeiro de Barros (SP-294); the HDI of each municipality is highlighted in red.

Figure 2
Geographic distribution of Human American Visceral Leishmaniasis cases in the microregion of Adamantina per year, from 2004 to 2018. Confirmed cases are expressed as incidences and represented by the Kernel density.

METHODS

Study area

As part of the municipalities of the Adamantina microregion, which is inserted in the Presidente Prudente mesoregion, Western region of the Sao Paulo State, are members of the Epidemiological Surveillance Group (GVE), region XIX - Marilia: Pacaembu, Florida Paulista, Adamantina, Mariapolis, Lucelia, Pracinha, Inubia Paulista, Sagres, Osvaldo Cruz, Parapua, Salmourao and Rinopolis. This microregion covers an area of 3,445,789 km ² and 149,066 inhabitants, based on the last census of the Brazilian Institute of Geography and Statistics (IBGE) ¹⁹ in 2010 and a range of HDI - Human Development Index - mostly within the higher range (from 0.700 to 0.799). This census was also used to record each municipality population. Annual rainfall data of the municipality of Adamantina was used for the meteorological referencing of the microregion, obtained by the Integrated Center of Agrometeorological Information of the Agronomic Institute of the Secretariat of Agriculture and Supply – CIIAGRO ²⁰ .

Human LV cases

This research is an observational, longitudinal and retrospective quantitative study. Excel files, provided by the Epidemiological Surveillance Center (CVE) “Prof. Alexandre Vranjac” of the Sao Paulo State Health Department, were used as collection instruments, referring to the HAVL, between 2004 to 2018, which were stratified and used according to the variables: date of notification, municipality code of residence, age, sex, residence address, neighborhood and classification of the diagnosis. Thereafter, the column data referring to the addresses of the respective confirmed cases were converted into geographic coordinates using the Google Geocoding API (Application Programming Interface). The ordered pairs of coordinates were added to the notification numbers, which allowed the georeferencing with a high precision, by plotting the cases on a surface previously delimited by the municipal limits using the ArcMap software (version 10.5, ESRI, California, USA). In order to generate a surface that approximates the studied phenomenon in a realistic way, the plotted confirmed cases were treated using the Kernel density phenomenon smoothing methodology for geographic representation purposes. The radius of influence that defines the degree of vicinity of the points was established in order to generate a 1000 square-meters area. All data treated with the Kernel tool had in its ratio the corresponding population of each municipality.

Statistical analysis

Qualitative variables were described by the distribution of absolute (N) and relative (%) frequencies. The proportion of cases in total and by region was analyzed by the incidence rate, which was calculated by dividing the number of confirmed cases by the aggregate of exposed population in the period and results were expressed per 100,000 inhabitants. The 95% confidence interval (95% CI) for the incidence rate was calculated using the equation p ± 1.96 × √p × q / n, where “p” represents the incidence rate, “q” represents the complementary incidence rate (q = 1 - p), and “n” represents the sample size. Significant differences between the incidence rates were considered when there was no overlapping of the 95% CI. The normality distribution was analyzed using the Shapiro-Wilk test. The relationship between quantitative variables was analyzed using the Pearson’s correlation test. The percentage of variation of the dependent variable (y-axis) explained by the variation of the independent variable (x-axis) was analyzed by the linear R2. The level of significance adopted was 5% (p ≤ 0.05). The data were analyzed using the SPSS software (version 22.0 IBM, Armonk, New York, USA).

Ethical approval

This research was approved by the Human Research Ethical Committee of Marilia Medical School, process nº 21430413.6.0000.5413.

RESULTS

During the analyzed period, 325 cases of HAVL were diagnosed in the municipalities of the Adamantina microregion, as shown in Table 1 . Of these cases, 316 (97.23%) had complete addresses information and were georeferenced. Nine cases (2.77%) were from rural areas, and among them, only three were correctly georeferenced. The highest absolute incidence was observed in the municipality of Adamantina with 108 confirmed cases, followed by 50 cases in the municipality of Osvaldo Cruz. During the studied period, no positive cases were registered in the municipality of Pracinha.

Table 1
Incidence of Human American Visceral Leishmaniasis per year in the municipalities of the microregion of Adamantina, Sao Paulo State, from 2004 to 2018.

Year	Case number	Incidence per 100.000 inhabitants (IC95%)
2004	2	1.3	(0.4 / 4.9)
2005	23	15.4	(10.3 / 23.1)
2006	28	18.8	(13.0 / 27.1)
2007	24	16.1	(10.8 / 24.0)
2008	51	34.2	(26.0 / 45.0)
2009	36	24.2	(17.4 / 33.4)
2010	41	27.5	(20.3 / 37.3)
2011	24	16.1	(10.8 / 24.0)
2012	22	14.8	(9.7 / 22.3)
2013	17	11.4	(7.1 / 18.3)
2014	16	10.7	(6.6 / 17.4)
2015	11	7.4	(4.1 / 13.2)
2016	12	8.1	(4.6 / 14.1)
2017	11	7.4	(4.1 / 13.2)
2018	7	4.7	(2.3 / 9.7)
Total	325	14.5	(13.0 / 16.2)

The temporal distribution of the microregion HAVL cases per year and their respective incidence calculations showed two peaks with a statistically significant difference in relation to the previous year, according to the established confidence interval, being one in 2005, with 15.4 cases per 100,000 inhabitants; and the other in 2008, with 34.2 cases per 100,000 inhabitants ( Table 1 ). It is noteworthy that 2008 was the year with the highest recorded incidence. However, in 2005, the municipality of Adamantina deserves special attention due to the incidence attributed to the microregion, as it concentrated around 18 confirmed cases, in absolute numbers, corresponding to an incidence of approximately 50 cases per 100,000 inhabitants. In 2008, the municipality of Pacaembu deserves to be mentioned, due to 14 confirmed cases in absolute numbers and an incidence of 100 cases per 100,000 inhabitants.

However, as of 2010, there was a decrease in the incidence of HAVL autochthonous cases in the microregion, and in 2018, only seven cases were reported in three municipalities, corresponding the lowest incidence in the period, of 4.7 cases per 100,000 inhabitants ( Table 1 ).

Figure 3 shows the expansion of human cases of HAVL in the microregion, which began in 2004 with the finding of amastigote forms in bone marrow aspirates of two patients living in the municipality of Adamantina. The following year, positive cases were detected in the cities of Lucelia and Rinopolis, 6 km and 48 km from Adamantina, respectively. In subsequent years, the disease was established, in the Northwest and East-Southeast directions, with cases in 2006 recorded in the municipality of Florida Paulista, in 2007 in the municipality of Pacaembu and, in 2008, in the municipalities of Inubia Paulista, Osvaldo Cruz and Parapua.

Figure 3
Geographic distribution of Human American Visceral Leishmaniasis cases in the Adamantina microregion from 2004 to 2018, expressed by the Kernel density ratio.

The density raster indicates a predominance of confirmed cases in urban areas in relation to rural areas, as shown in Figure 3 . Municipalities that are crossed by the Comandante Joao Ribeiro de Barros highway (SP – 294), such as Florida Paulista, Adamantina, Pacaembu and Lucelia, had the highest incidence rates, in that order. No statistically significant difference was observed for the incidence of HAVL between the municipalities of Adamantina, Florida Paulista, Lucelia, Pacaembu and Inubia Paulista, which are interconnected by the Comandante Joao Ribeiro de Barros highway. However, the municipalities of Adamantina, Florida Paulista and Pacaembu showed a statistically significant (95% CI) difference regarding the incidence of confirmed cases of HAVL in relation to the municipalities of Rinopolis, Osvaldo Cruz and Parapua. The municipality of Mariapolis also had a low incidence rate, however, it showed a statistically significant difference only in the municipalities of Adamantina and Florida Paulista.

The municipalities of Adamantina and Florida Paulista had a higher incidence and statistical significance in the microregion ( Table 2 ), whereas, Rinopolis and Parapua had lower incidences and statistical significance. Nevertheless, given the high amplitude of confidence intervals of the municipalities with the lowest populations, the municipalities of Osvaldo Cruz, Inubia Paulista, Lucelia and Parapua did not present a statistically significant difference in relation to Mariapolis and Rinopolis. In relation to sex, only in the municipality of Osvaldo Cruz a significant difference was observed. In the microregion, this difference could not be established, emphasizing that the incidence calculations were performed based on the population according to the 2010 IBGE sense.

Table 2
Relationship between the incidence of Human American Visceral Leishmaniasis by municipality, sex and the total of the microregion of Adamantina, Sao Paulo State, from 2004 to 2018.

Municipalities	Female		Male		Total
	Case Nº	Incidence*	Case Nº	Incidence*	Case Nº	Incidence*
Adamantina	45	17.2	63	25.6	108	21.3
		(12.9 / 23.0)		(20.0 / 32.7)		(17.6 / 25.7)
Florida Paulista	21	24.2	24	22.7	45	23.3
		(15.8 / 36.9)		(15.2 / 33.7)		(17.5 / 31.2)
Inubia Paulista	1	3.7	4	14.7	5	9.2
		(0.7 / 20.9)		(5.7 / 37.9)		(3.9 / 21.5)
Lucelia	20	14.3	24	15.2	44	14.8
		(9.3 / 22.1)		(10.2 / 22.5)		(11.0 / 19.8)
Mariapolis	1	3.5	3	9.9	4	6.8
		(0.6 / 19.9)		(3.4 / 29.1)		(2.6 / 17.5)
Osvaldo Cruz	15	6.5	35	15.0	50	10.8
		(3.9 / 10.7)		(10.8 / 20.9)		(8.2 / 14.2)
Pacaembu	13	15.9	27	23.2	40	20.2
		(9.3 / 27.2)		(16.0 / 33.8)		(14.8 / 27.5)
Rinopolis	4	5.5	4	5.3	8	5.4
		(2.1 / 14.0)		(2.1 / 13.6)		(2.7 / 10.6)
Sagres	1	5.5	1	5.6	2	5.6
		(1.0 / 31.4)		(1.0 / 31.7)		(1.5 / 20.3)
Salmourao	5	14.2	4	10.8	9	12.5
		(6.1 / 33.3)		(4.2 / 27.7)		(6.6 / 23.7)
Parapua	6	7.5	4	4.8	10	6.1
		(3.5 / 16.4)		(1.9 / 12.4)		(3.3 / 11.3)
Microregion	132	12.3	193	16.6	325	14.5
		(10.4 / 14.6)		(14.4 / 19.1)		(13.0 / 16.2)

• The Municipality of Pracinha is not described in Table 2 because it did not present confirmed cases during the studied period; *Incidence per 100,000 inhabitants (95% CI)

The Pearson’s correlation coefficient revealed no influence of annual rainfall index on the HVAL incidence, as seen in Figure 4 . Regarding the demographic analysis, it was observed that the age group from 0 to 9 years old had a higher incidence compared to the others, with 47.4 cases per 100,000 inhabitants, followed by the age group above 59 years, with an incidence of 14.7 per 100,000 inhabitants. The first presented a statistical significance difference in relation to all other groups, and the second presented a statistical significant difference in relation to the group from 20 to 59 years old ( Table 2 ).

Figure 4
Relationship between rainfall indexes and Human American Visceral Leishmaniasis incidence in Adamantina, Sao Paulo State, and Adamantina, Sao Paulo State microregion, from 2004 to 2015. The Pearson’s correlation coefficient (r); p value calculated by the Pearson’s correlation test. Linear R2 = percentage of variation of the dependent variable (y-axis) explained by the variation of the independent variable (x-axis).

DISCUSSION

In 2003, dogs infected by visceral leishmaniasis were reported in Adamantina microregion for the first time in the municipality of Lucelia, then, one-year latter, in 2004, the municipality of Adamantina presented the first confirmed human case, corroborating that infection in dogs precede the occurrence of HAVL, in concordance with its previously described classical dispersion pattern, as well as their importance from an epidemiological point of view ¹¹ , ²¹ , ²² .

Since then, HAVL has spread to 11 of the 12 municipalities in this area, reaching the highest number of notifications in 2008, with an incidence of 34.2 cases per 100,000 inhabitants. Although the municipalities of this microregion are relatively similar in terms of HDI, and also regarding the climate and altitude, all aspects that show importance within the scope of risk factors for HALV transmission, were not relevant in this study when compared to geographic factors ¹⁴ .

The georeferencing of HAVL cases in the Adamantina microregion argues in favor of a correlation between municipalities with higher disease incidence and the proximity to the secondary Highway Comandante Joao Ribeiro de Barros (SP-294), suggesting that this highway has an important role in HAVL dispersion.

Although some studies ¹⁴ , ¹⁵ , ²³ have pointed Marechal Rondon highway as the main pathway of spreading the disease due to its intense flow of people and goods, departing from Aracatuba and heading towards Bauru, in the present study it was observed that the secondary roads were equally important in the dispersion of the disease in the Western Sao Paulo region, since the first reports of human cases in the city of Bauru in 2003 ²⁴ and in the municipality of Adamantina in 2004. Hence, our results suggest that the spread of HAVL in this West region of Sao Paulo State followed a centrifugal dissemination pattern, as observed by Toledo et al . ²⁵ in Araguaina, Tocantins State.

The spatial analysis showed that the spread of HAVL along the SP-294 road occurred in the East-Southeast direction, with progressive dispersion of autochthonous cases. Although other municipalities, which are not part of the Adamantina microregion, but belong to GVE-XIX and are interconnected by this highway were not the subject of this study, it was observed that, in 2010, the first autochthonous case of HAVL was notified in the municipality of Bastos; in 2012, in the municipality of Tupa; in 2013, in Quintana and Pompeia; in 2014 in the city of Marilia; and, finally, in 2019, in the municipality of Garca (CVE).

The comparison of GVE-XIX (Marilia), which belongs to the Adamantina microregion, with the other neighboring GVEs, such as GVE-XI (Aracatuba), GVE-XV (Bauru), GVE-XXI (Presidente Prudente) and GVE-XXII (Presidente Venceslau), found a lower incidence in the first municipality in relation to the others. However, considering the number of reported cases, as well as the populations, the analysis of the Adamantina microregion found higher incidence rates, both within its GVE and in relation to the borderline municipalities, demonstrating the importance of this microregion as an incipient endemic region in Sao Paulo State. This fact probably occurred due to the proximity between the municipalities interconnected by SP-294, facilitating the intense flow of residents, as the most affected municipalities that are crossed by it presented statistically significant differences in the incidence of cases when compared with the municipalities of Mariapolis, Rinopolis and Pracinha, and the latter did not present confirmed cases of the disease. In addition, our Kernel density raster follows the SP-294 direction, showing the most affected areas in the accumulated period.

An elevated rainfall index is an important environmental factor of HAVL as it significantly increases the population of sandflies ²⁶ . Despite that, our data revealed no correlation, either positive or negative, between rainfall indexes and HAVL incidences, decreasing the most relevant environmental factor related to the climate involved in the disease, both in propagation and in perpetuation.

Another factor that has probably influenced the high incidence of the disease in this microregion is the sociodemographic characteristics of the municipalities. They are small, with an average of 12,422 inhabitants that preserve rural characteristics and habits, most residences have fruit trees and livestock, such as chickens and ducks that are used as natural predators of the endemic scorpions in the region ²⁷ . Thus, there is an ideal environment for the development and dispersal of the vector insect, similar to the situation reported by Vieira et al . ²⁸ in the municipality of Birigui.

As mentioned by other authors ²⁹ , ³⁰ , it was found that most cases of HAVL were concentrated in the urban perimeter of the municipalities and, despite a higher prevalence of the disease in males in others regions of Sao Paulo State ²³ , ³¹ , no significant difference between genders was observed in the studied microregion, corroborating the results obtained by Vieira et al . ²⁸ in Birigui, Sao Paulo State, and Moura et al . ³² in the island of Sao Luis, Maranhao State. The HAVL equal incidence between genders can be attributed to the predominance of urban cases, an environment in which there is no tendency for men to be more exposed to a harsh environment than women ²⁸ , unlike the situation in rural areas ³³ .

Regarding the age group, a higher incidence was found in children up to 9 years old and in adults over 59 years old, but the highest lethality was observed in adults between 20 and 59 years old. According to Queiroz et al . ³⁴ , the greater the number of cases in a region, the higher the vulnerability of infected children, who have an immature cellular immune response, in the same way of the elderly, who also show an immune response decline with advancing age, beyond comorbidities, such as chronic diseases that worsen the prognosis in these individuals and justify the increased incidence in both age ranges ³⁵ , ³⁶ .

Although there are reports of HAVL seasonality ²³ , this research could not detect this phenomenon, reported by Lopes et al . ³⁷ in Belo Horizonte, Minas Gerais State. Among some possibilities, we highlight the short period of time since the onset of the disease in this microregion.

After the establishment of HAVL prevention and control measures, a decrease in the number of cases has been observed since 2010, what may have occurred in association with the “Visceral Leishmaniasis Control Program”, which classifies the regions according to areas with confirmed HAVL cases and silent ones. They can also be classified as vulnerable, as they are contiguous municipalities with others with confirmed cases of HAVL and / or are interconnected by the same important road axis between endemic municipalities ³⁸ . Given the geographic distribution shown in this research, municipalities in the West of Sao Paulo State contiguous to all those cities, or presenting with confirmed cases, as well as those connected by SP-294 and that have not yet presented cases of HAVL, fall into this category of vulnerability, requiring preventive control strategies ²⁸ , ³⁷ , ³⁹ .

Furthermore, the implementation of prevention and control measures at the municipal level are hampered by the intermittency of the measures, sometimes due to financial factors, sometimes due to dengue epidemics, as these situations depend on the same team and municipal resources used in health problems. In addition, the prevalence of American visceral leishmaniasis-seropositive dogs whose tutors resist to deliver them to euthanasia and various other reasons such as the difficulties in the environmental management, also jeopardize the implementation of preventive and control HAVL measures ¹¹ , ²² .

However, it can still be presumed that, in general, the control measures employed in the municipalities of the Adamantina microregion were competent after the year 2010, given the subsequent decrease in HAVL incidence rates, as observed by Costa et al . ²² in Aracatuba, between 2007 and 2009. At the same time, studies that identify the places with the highest incidences and risk factors, as well as the characterization of the epidemiological regional characteristics can significantly contribute to the diagnosis-treatment process ¹⁸ , ²⁸ .

In view of the above, the georeferencing of cases in this microregion can assist in the identification of potentially high-risk regions for the occurrence of the disease, as a better understanding of the disease and its expansion process, under a space-time characterization, facilitates the adoption of measures aimed at controlling and preventing. Thus, understanding the spatial and temporal dissemination of HAVL and associate it to notorious environmental characteristics including the geographical space, culture, habits and socioeconomic conditions, enables to construct regional epidemiological profiles in order to develop rational control strategies.

REFERENCES

• ¹ Ashford RW. The leishmaniases as emerging and reemerging zoonoses. Int J Parasitol. 2000;30:1269-81.
• ² Herwaldt BL. Leishmaniasis. Lancet (London, England). 1999;354:1191-9.
• ³ Organização Pan-Americana da Saúde. Organização Mundial da Saúde. Leishmanioses: informe epidemiológico das Américas. [cited 2020 Sep 15]. Available from: https://iris.paho.org/bitstream/handle/10665.2/50505/2019-cde-leish-informe-epi-das-americas.pdf?sequence=2&isAllowed=y
  » https://iris.paho.org/bitstream/handle/10665.2/50505/2019-cde-leish-informe-epi-das-americas.pdf?sequence=2&isAllowed=y
• ⁴ Deane LM, Deane MP. Visceral leishmaniasis in Brazil: geographical distribution and transmission. Rev Inst Med Trop Sao Paulo. 1962;4:198-212.
• ⁵ Costa JM, Viana GM, Saldanha AC, Nascimento MD, Alvim AC, Burattini MN, et al. Leishmaniose Visceral no Estado do Maranhão, Brasil: a evolução de uma epidemia. 1995;11:321-4.
• ⁶ Costa CH, Pereira HF, Araújo MV. Epidemia de leishmaniose visceral no Estado do Piauí, Brasil, 1980-1986. Rev Saude Publica. 1990;24:361-72.
• ⁷ Rodrigues AC, Melo AC, Júnior AD, Franco SO, Rondon F, Bevilaqua CM. Epidemiologia da leishmaniose visceral no município de Fortaleza, Ceará. Pesq Vet Bras. 2017;37:1119-24.
• ⁸ Dantas-Torres F. Situação atual da epidemiologia da leishmaniose visceral em Pernambuco. Rev Saude Publica. 2006;40:537-41.
• ⁹ Dantas-Torres F, Brandão-Filho SP. Expansão geográfica da leishmaniose visceral no Estado de Pernambuco. Rev Soc Bras Med Trop. 2006;39:352-6.
• ¹⁰ Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Departamento de Vigilância Epidemiológica. Guia de vigilância epidemiológica. 7ª ed. Brasília: Ministério da Saúde; 2009.
• ¹¹ von Zuben AP, Donalísio MR. Dificuldades na execução das diretrizes do Programa de Vigilância e Controle da Leishmaniose Visceral em grandes municípios brasileiros. Cad Saude Publica. 2016;32:e00087415.
• ¹² Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Coordenação-Geral de Desenvolvimento da Epidemiologia em Serviços. Guia de vigilância em saúde. 2ª ed. Brasília: Ministério da Saúde; 2017.
• ¹³ Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde.Vigilância em saúde no Brasil 2003/2019: da criação da Secretaria de Vigilância em Saúde aos dias atuais. Bol Epidemiol. 2019;N Esp:1-154.
• ¹⁴ Sevá AP, Mao L, Galvis-Ovallos F, Lima JM, Valle D. Risk analysis and prediction of visceral leishmaniasis dispersion in São Paulo State, Brazil. PLoS Negl Trop Dis. 2017;11:e0005353.
• ¹⁵ Casanova C, Colla-Jacques FE, Hamilton JG, Brazil RP, Shaw JJ. Distribution of Lutzomyia longipalpis chemotype populations in São Paulo state, Brazil. PLoS Negl Trop Dis. 2015;9:e0003620.
• ¹⁶ Cardim MF, Rodas LA, Dibo MR, Guirado MM, Oliveira AM, Chiaravalloti Neto, F. Introduction and expansion of human American visceral leishmaniasis in the state of São Paulo, Brazil, 1999-2011. Rev Saude Publica. 2013;47:691-700.
• ¹⁷ Cardim MF, Vieira CP, Chiaravalloti-Neto, F. Spatial and spatiotemporal occurrence of human visceral leishmaniasis in Adamantina, State of São Paulo, Brazil. Rev Soc Bras Med Trop. 2015;48:716-23.
• ¹⁸ Barcellos CC, Sabroza PC, Peiter P, Iñiguez Rojas L. Organização espacial, saúde e qualidade de vida: análise espacial e uso de indicadores na avaliação de situações de saúde. Inf Epidemiol SUS. 2002;11:129-38.
• ¹⁹ Instituto Brasileiro de Geografia e Estatística. Cidades@. [cited 2020 Sep 15]. Available from: https://cidades.ibge.gov.br/
  » https://cidades.ibge.gov.br/
• ²⁰ São Paulo. Governo do Estado. Centro Integrado de Informações Agrometeorológicas. Balanço hídrico semanal. [cited 2020 Sep 15]. Available from: http://www.ciiagro.sp.gov.br/ciiagroonline/Listagens/BH/LBalancoHidricoLocal.asp
  » http://www.ciiagro.sp.gov.br/ciiagroonline/Listagens/BH/LBalancoHidricoLocal.asp
• ²¹ Camargo-Neves V. A leishmaniose visceral americana no estado de São Paulo: situação atual. Bol Epidemiol Paulista. 2007;4:12-4.
• ²² Costa DN, Bermudi PM, Rodas LA, Nunes CM, Hiramoto RM, Tolezano JE, et al. Human visceral leishmaniasis and relationship with vector and canine control measures. Rev Saude Publica. 2018;52:92.
• ²³ Cardim MF, Guirado MM, Dibo MR, Chiaravalloti Neto F. Leishmaniose visceral no estado de São Paulo, Brasil: análise espacial e espaço-temporal. Rev Saude Publica. 2016;50:48.
• ²⁴ Souza VA, Cortez LR, Dias RA, Amaku M, Ferreira Neto JS, Kuroda RB, et al. Space-time cluster analysis of American visceral leishmaniasis in Bauru, São Paulo State, Brazil. Cad. Saude Publica. 2012;28:1949-64.
• ²⁵ Toledo CR, Almeida AS, Chaves SA, Sabroza PC, Toledo LM, Caldas JP. Vulnerability to the transmission of human visceral leishmaniasis in a Brazilian urban area. Rev Saude Publica. 2017;51:49.
• ²⁶ Macedo IT, Bevilaqua CM, Morais N, Sousa L, Linhares F, Amóra SS, et al. Sazonalidade de flebotomíneos em área endêmica de leishmaniose visceral no município de Sobral, Ceará, Brasil. Cien Anim. 2008;18:67-74.
• ²⁷ Martins KP, Garcia DA, Cortezi AM, Gomes DE. Escorpionismo: revisão de literatura. Rev Cien. 2018;1.
• ²⁸ Vieira CP, Oliveira AM, Rodas LA, Dibo MR, Guirado MM, Chiaravalloti Neto F. Temporal, spatial and spatiotemporal analysis of the occurrence of visceral leishmaniasis in humans in the City of Birigui, State of São Paulo, from 1999 to 2012. Rev Soc Bras Med Trop. 2014;47:350-8.
• ²⁹ Alves WA, Fonseca DS. Leishmaniose visceral humana: estudo do perfil clínico-epidemiológico na região leste de Minas Gerais, Brasil. J Health Biol Sci. 2018;6:133-9.
• ³⁰ Ortiz RC, Anversa LJ. Epidemiologia da leishmaniose visceral em Bauru, São Paulo, no período de 2004 a 2012: um estudo descritivo. Epidemiol Serv Saude. 2015;24:97-104.
• ³¹ Madalosso G, Fortaleza CM, Ribeiro AF, Cruz LL, Nogueira PA, Lindoso JA. American visceral leishmaniasis: factors associated with lethality in the state of São Paulo, Brazil. J Trop Med. 2012;2012:281572.
• ³² Moura GS, Santos AM, Aquino DM, Silva AA, Caldas AJ. Factors associated with asymptomatic infection in family members and neighbors of patients with visceral leishmaniasis. Cad Saude Publica. 2012;28:2306-14.
• ³³ Brazuna JC, Silva EA, Brazuna JM, Domingos IH, Chaves N, Honer MR, et al. Profile and geographic distribution of reported cases of visceral leishmaniasis in Campo Grande, State of Mato Grosso do Sul, Brazil, from 2002 to 2009. Rev Soc Bras Med Trop. 2012;45:601-6.
• ³⁴ Queiroz MJ, Alves JG, Correia JB. Visceral leishmaniasis: clinical and epidemiological features of children in an endemic area. J Pediatr. 2004;80:141-6.
• ³⁵ Oliveira JM, Fernandes AC, Dorval ME, Alves TP, Fernandes TD, Oshiro ET, et al. Mortalidade por leishmaniose visceral: aspectos clínicos e laboratoriais. Rev Soc Bras Med Trop. 2010;43:188-93.
• ³⁶ Santos MA, Rodrigues SL, Nascimento AL, Rodrigues JS, Góes MA. Leishmaniose Visceral: características clínico-epidemiológicas de casos e óbitos no estado de Sergipe. Rev Epidemiol Controle Infecç. 2018;8:428-34.
• ³⁷ Lopes EG, Magalhães DF, Silva JA, Haddad JP, Moreira EC. Distribuição temporal e espacial da leishmaniose visceral em humanos e cães em Belo Horizonte-MG, 1993 a 2007. Arq Bras Med Vet Zootec. 2010;62:1062-71.
• ³⁸ Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Departamento de Vigilância Epidemiológica. Manual de vigilância e controle da leishmaniose visceral. Brasília: Ministério da Saúde; 2014.
• ³⁹ Camargo-Neves VL, Katz G, Rodas LA, Poletto DW, Lage LC, Spínola RM, et al. Utilização de ferramentas de análise espacial na vigilância epidemiológica de leishmaniose visceral americana-Araçatuba, São Paulo, Brasil, 1998-1999. Cad Saude Publica. 2001;17:1263-7.

• FUNDING
  This study was supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), grant Nº 2016/14514-4.

Publication Dates

• Publication in this collection
  30 Oct 2020
• Date of issue
  2020

History

• Received
  20 May 2020
• Accepted
  15 Sept 2020