Abstract

It is an ecological study that analyzed the time trend of visceral leishmaniasis incidence rates in Brazil using segmented time regression by joinpoints. There was a decreasing incidence rate of this disease in the country with an average annual percent change (AAPC) of -5 (CI95%: -9.1; -0.6) and a reduction of 1.69 cases/100 thousand inhabitants in 2007, and 0.91/100 thousand inhabitants in 2020. The Central-West region showed the highest reduction percent (AAPC: -9.1; CI95%: -13.8; -4.3), followed by the Southeast region (AAPC: -8.7; -14.6; -2.5). The North and South regions showed the largest number of joinpoints in the time series. The highest incidences were recorded in the male population, however, stable (AAPC: 2.14; CI95%: -8.3; 0). In the age group analysis, the trend was decreasing for the groups from 0 to 4 years old (AAPC: -7.7; CI95%: -12.6; -2.4), 5 to 9 years old (AAPC: -7.3; CI95%: -13.6; -0,4) and 10 to 14 years old (AAPC: -5.5; CI95%: -10.3; -0.3). It was found that although Visceral Leishmaniasis is an endemic disease in Brazil, there was a decrease in its incidence rate from 2007 to 2020.

Key words:
Visceral leishmaniasis; Epidemiological monitoring; Temporal distribution; Time series study; Ecological studies

Resumo

Trata-se de um estudo ecológico que analisou a tendência temporal das taxas de incidência de leishmaniose visceral no Brasil mediante regressão temporal segmentada por pontos de inflexão. Observou-se tendência de decréscimo na taxa de incidência dessa patologia no país, com variação variação percentual média anual (average annual percent change - AAPC) de -5 (IC95%: -9,1; -0,6) e redução de 1,69 casos/100 mil habitantes em 2007, para 0,91/100 mil habitantes em 2020. A região Centro-Oeste apresentou a maior redução do AAPC (AAPC: -9,1; IC95%: -13,8; -4,3), seguida da região Sudeste (AAPC: -8,7; -14,6; -2,5). As regiões Norte e Sul apresentaram o maior número de pontos de inflexão (joinpoint) na série temporal. As maiores incidências foram registradas na população masculina, porém com tendência estacionária (AAPC: 2,14; IC95%: -8,3; 0). Na análise por faixa etária, a tendência foi decrescente nos grupos de 0 a 4 anos (AAPC: -7,7; IC95%: -12,6; -2,4), 5 a 9 anos (AAPC: -7,3; IC95%: -13,6; -0,4) e de 10 a 14 anos (AAPC: -5,5; IC95%: -10,3; -0,3). Verificou-se que, apesar de a leishmaniose visceral se tratar de uma doença endêmica no Brasil, houve declínio na sua taxa de incidência no período de 2007 a 2020.

Palavras-chave:
Leishmaniose visceral; Monitoramento epidemiológico; Distribuição temporal; Estudos de séries temporais; Estudos ecológicos

Introduction

Visceral leishmaniasis (VL) or kala-azar is a zoonosis of worldwide distribution caused by protozoa of the genus Leishmania that are transmitted by sandflies¹1 Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, Janin J, Boer M. Leishmaniasis worldwide and global estimates of its incidence. PLoS One 2012; 7(5):e35671.^,22 Jones CM, Welburn SC. Leishmaniasis beyond East Africa. Front Vet Sci 2021; 8:618766.. The main urban reservoirs of VL are dogs³3 Ministério da Saúde (MS). Guia de vigilância epidemiológica Brasília: MS; 2019..

This disease is considered a problem of Public Health of international importance, especially affecting tropical countries, whose socio-environmental conditions favor the development and proliferation of Lutzomyia longipalpis, the main vector ofLeishmania chagasi, the causative agent of VL in Latin America²2 Jones CM, Welburn SC. Leishmaniasis beyond East Africa. Front Vet Sci 2021; 8:618766.

3 Ministério da Saúde (MS). Guia de vigilância epidemiológica Brasília: MS; 2019.

4 World Health Organization (WHO). Status of endemicity of visceral leishmaniasis worldwide [Internet]. 2018. [cited 2022 maio 7]. Available from: https://www.who.int/leishmaniasis/burden/GHO_VL_2018.pdf
https://www.who.int/leishmaniasis/burden... ^-55 Cruz CSS, Barbosa DS, Oliveira VC, Cardoso DT, Guimarães NS, Carneiro M. Factors associated with human visceral leishmaniasis cases during urban epidemics in Brazil: a systematic review. Parasitology 2021; 148(6):639-647..

It is a disease of wide distribution in Asia, Europe, the Middle East, Africa, and the Americas, with the highest incidence in rural and suburban areas⁶6 Safavi M, Eshaghi H, Hajihassani Z. Visceral Leishmaniasis: Kala-azar. Diagn Cytopathol 2021; 49(3):446-448.^,77 Bi K, Chen Y, Zhao S, Kuang Y, John Wu CH. Current visceral leishmaniasis research: a research review to inspire future study. Biomed Res Int 2018; 2018:9872095 ..

VL is the most severe form of leishmaniasis, being potentially fatal due to its multisystem involvement and the tropism of its agentS for hematopoietic organs, such as the spleen, liver, and bone marrow³3 Ministério da Saúde (MS). Guia de vigilância epidemiológica Brasília: MS; 2019.^,88 Graepp-Fontoura I, Barbosa DS, Nascimento LFC, Fontoura VM, Ferreira AGN, Santos FAAS, Sousa BS, Santos FS, Santos-Neto M, Santos LH, Silva ALA. Epidemiological aspects and spatial patterns of human visceral leishmaniasis in Brazil. Parasitology 2020; 147(14):1665-1677., and affects more frequently socially vulnerable populations and individuals with immunosuppression conditions, such as malnourished children, patients with HIV/Aids infection, elderly people with comorbidities, among other conditions that impair the host’s immune response³3 Ministério da Saúde (MS). Guia de vigilância epidemiológica Brasília: MS; 2019.^,44 World Health Organization (WHO). Status of endemicity of visceral leishmaniasis worldwide [Internet]. 2018. [cited 2022 maio 7]. Available from: https://www.who.int/leishmaniasis/burden/GHO_VL_2018.pdf
https://www.who.int/leishmaniasis/burden... ^,99 Araújo VEM, Pinheiro LC, Almeida MCM, Menezes FC, Morais MHF, Reis IA, Assunção RM, Carneiro M. Relative risk of visceral leishmaniasis in Brazil: a spatial analysis in urban area. PLoS Negl Trop Dis 2013; 7(11):e2540..

It is estimated that 400 thousand of new cases of this disease occur annually across the globe¹1 Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, Janin J, Boer M. Leishmaniasis worldwide and global estimates of its incidence. PLoS One 2012; 7(5):e35671.^,1010 Bezerra JMT, Araújo VEM, Barbosa DS, Martins-Melo FR, Werneck GL, Carneiro M. Burden of leishmaniasis in Brazil and federated units, 1990-2016: findings from Global Burden of Disease Study 2016. PLoS Negl Trop Dis 2018; 12(9):e0006697.. In the Americas, 13 countries are endemic, and from 2001 to 2020, 67,922 new cases were recorded, with an average of 3,400 occurrences per year. In 2020, of the total number of new cases in the American continent, 97% were reported by Brazil and the others by Argentina, Bolivia, Colombia, Paraguay, Venezuela, and Uruguay⁴4 World Health Organization (WHO). Status of endemicity of visceral leishmaniasis worldwide [Internet]. 2018. [cited 2022 maio 7]. Available from: https://www.who.int/leishmaniasis/burden/GHO_VL_2018.pdf
https://www.who.int/leishmaniasis/burden... ^,1111 Organização Pan-Americana da Saúde (OPAS). Leishmanioses: Informe Epidemiológico das Américas [internet]. 2021. [cited 2022 maio 7]. Available from: https://iris.paho.org/handle/10665.2/55386
https://iris.paho.org/handle/10665.2/553... .

In the last three decades, however, many aspects of the VL ecoepidemiology were reviewed and enabled the development of surveillance measures by the control bodies, such as the incorporation of actions and activities in the Visceral Leishmaniasis Surveillance and Control Program (VLSCP). Since 2011, the Regional Leishmaniasis Program (RLP) by the Pan American Health Organization (PAHO)¹¹11 Organização Pan-Americana da Saúde (OPAS). Leishmanioses: Informe Epidemiológico das Américas [internet]. 2021. [cited 2022 maio 7]. Available from: https://iris.paho.org/handle/10665.2/55386
https://iris.paho.org/handle/10665.2/553... has been monitoring the behavior of this disease in the Americas region with planning and control measures.

The identification of territorial units of greater epidemiological importance is a factor that needs to be considered in the planning of VL promotion, prevention, surveillance, and control measures by health agencies, especially in endemic regions, as is the case of the national territory. Despite this reality, no recent study was identified in Brazil on the behavior of the disease over the years in its different regions, Federation units, and according to sociodemographic characteristics. Therefore, this study sought to analyze the trends of VL incidence in the country from 2007 to 2020 by major national regions and Federation units, and according to the gender and age group of those affected.

Method

This is an ecological time series study to analyze the trends of VL incidence rates in the population of Brazil, from 2007 to 2020.

The study scenario is Brazil, the largest nation in South America in terms of land area and population. Information on new VL cases (CID 10 - B55.0) was collected from the Ministry of Health’s Information System on Notifiable Diseases (SINAN Net) database, available on the website of the Department of Informatics of the Unified Health System (DATASUS; https://datasus.saude.gov.br/informacoes-de-saude-tabnet/).

The collection of data took place on March 20, 2022. SINAN Net is a system that aims to collect, transmit, and disseminate data routinely generated by the Epidemiological Surveillance System of the three levels of government in Brazil, through a computerized network, to support the investigation process and provide subsidies for the analysis of epidemiological surveillance information on compulsory notifiable diseases, as it is the case with VL¹²12 Ministério da Saúde (MS). Informações de saúde [Internet]. 2022. [acessado 2022 mar 20]. Disponível em: https://datasus.saude.gov.br/informacoes-de-saude-tabnet/
https://datasus.saude.gov.br/informacoes... .

Data on the population residing in the country by gender and group age were obtained through demographic information from the Brazilian Institute of Geography and Statistics (IBGE), which is information based on censuses and the inter-census population estimates (www.ibge.gov.br/estatisticas/sociais/populacao).

VL cases were tabulated using DATASUS Tabnet and Microsoft Excel Office 2013 was used to organize the database, build indicators, and create tables.

VL incidence rates were calculated as the number of new VL new cases reported to SINAN divided by the population residing in the location and year, multiplied by 100 thousand inhabitants. The rates were calculated for the country and according to regions North, Northeast, Central-West, Southeast, and South; federation units (FU - North Region: Rondônia, Acre, Amazonas, Roraima, Pará, Amapá, Tocantins; Northeast Region: Maranhão, Piauí, Ceará, Rio Grande do Norte, Paraíba, Pernambuco, Alagoas, Sergipe, Bahia; Central-West Region: Mato Grosso do Sul, Mato Grosso, Goiás, Distrito Federal, Southeast Region: Minas Gerais, Espírito Santo, Rio de Janeiro, São Paulo; South Region: Paraná, Santa Catarina, Rio Grande do Sul), gender, and age groups (in years, younger than 4; 5 to 9; 10 to 14; 15 to 19; 20 to 39; 40 to 59; 60 to 79; and 80 or older).

Then, a joinpoint regression analysis model was used for the time analysis. This model verifies if one line with multiple segments is statistically better to describe the time evolution of a data set than a straight line, or with fewer segments. In this way, it allows identify the trend of the indicator (increasing, stable, or decreasing) and the points where there is a change in this trend, allowing us to know the annual percent change (APC) and that of the entire period, called average annual percent change (AAPC)¹³13 Kim HJ, Fay MP, Feuer EJ, Midthune DN. Permutation tests for joinpoint regression with applications to cancer rates. Stat Med 2000; 19(3):335-351.^,1414 National Cancer Institute. Joinpoint Help Manual, 4.8.0.1 [Internet]. 2020. [accessed 2022 maio 7]. Available from: https://surveillance.cancer.gov/help/joinpoint/
https://surveillance.cancer.gov/help/joi... .

The incidence rates were considered a dependent variable (Y) and calendar years as an independent variable (X). For each trend detected, 95% confidence intervals (CI95%) and p values (5% significance level) were considered. Non-significant p values resulted in a stability trend (accepting the null hypothesis that the incidence rates have not changed over the years). On the other hand, significant p values resulted in a classification of growth trend (positive APC or AAPC) and decreasing trend (negative APC or AAPC)¹³13 Kim HJ, Fay MP, Feuer EJ, Midthune DN. Permutation tests for joinpoint regression with applications to cancer rates. Stat Med 2000; 19(3):335-351.^,1414 National Cancer Institute. Joinpoint Help Manual, 4.8.0.1 [Internet]. 2020. [accessed 2022 maio 7]. Available from: https://surveillance.cancer.gov/help/joinpoint/
https://surveillance.cancer.gov/help/joi... .

All tests were performed on Joinpoint Regression Program, version 4.9.0.1, and to minimize the effect of eventual autocorrelations, the option “fit an autocorrelated errors model based on the data” available on the software was used.

Results

Between 2007 and 2020, 45,270 new cases of VL were recorded, with a predominance in the Northeast region (52.78%). The male gender was responsible for 64.13% of the occurrences. As for age, the most affected group was from 1 to 4 years old (34.33%). The year with the highest incidence of VL was 2017, with 1.98 cases/per 100 thousand inhabitants.

The joinpoint regression analysis model pointed out, for the VL incidence rate in Brazil, a time behavior with two APCs: the first segment, stationary between 2007 and 2018 (APC:-0.5; CI95%: -2.6; 1.6; p > 0.05); the second also between 2018 and 2020 (APC: -26; CI95%: -45.9; 1.4; p > 0.05) (Figure 1). When analyzing the total period, however, there was a significant decreasing trend in the incidence of the disease in the national territory, with an average annual percent change (AAPC) of -5 (CI95%: -9.1; -0.6; p < 0.05), and a reduction from 1.69 cases/100 thousand inhabitants in 2007 to 0.91/100 thousand inhabitants in 2020 (Table 1).

In the time analysis, two major administrative regions presented decreasing trend, highlighting the Central-West region with the largest average annual percent change (AAPC: -9.1; CI95%: -13.8; -4.3; p < 0.05) and Southeast (AAPC: -8.7; -14.6; -2.5; p < 0.05) with the lowest. There was an increasing trend in the South region (AAPC: 21.7; CI95%: 5.7; 40.2; p < 0.05) and stability in the North (AAPC: -8.6; CI95%: -20.4; 5.0; p > 0.05) and Northeast regions (AAPC: -4.2; CI95%: -9.4; 1.4; p > 0.05) (Table 1). The regression model also showed that the North and Southeast regions were the ones that presented the highest number of segments in the time series (joinpoints), representing the greatest oscillation in rated over the years (Figure 1).

Concerning the analysis by federation unit (FU), the time modeling showed that five regions and the Federal District presented decreasing trend, with one from the North region, one from the Northeast, one from the Southeast, and tho from the Central-West (two states and the Federal District). On the other hand, 16 presented stationary patterns, with four in the North, seven in the Northeast, two in the Southeast, two in the South, and 1 in the Central-West; and four showed an increasing trend, one in each region (except for the Central-West region) (Table 1).

From 2007 to 2020, the autochthony of the disease was confirmed in 26 federation units, distributed in five Brazilian regions, with the Northeast region being responsible for the largest number of cases in the country (n = 23,893; 44.6%). Only the state of Acre, in the North region, did not present VL cases during the studied period. The state of Tocantins, in turn, concentrated the highest incidence of VL throughout the data series.

Higher VL incidence rates were observed in the male population. In this group, the incidence trend was stationary (AAPC: 2.14; CI95%: -8.3; 0); in females, in turn, there was a decrease (AAPC: -6.3; CI95%: -20.8; -1.5) (Table 2). When stratifying the indicator according to age group, a decreasing trend was observed in the groups from 0 to 4 years old (AAPC: -7.7; CI95%: -12.6; -2.4; p < 0.05), 5 to 9 years old (AAPC: -7.3; CI95%: -13.6; -0.4; p < 0.05), and from 10 to 14 years old (AAPC: -5.5; CI95%: -10.3; -0.3; p < 0.05). In the other strata, the trends present a stationary behavior (Table 2).

Discussion

In Latin America, more than 95% of VL cases occur in Brazil, most of them in the Northeast region¹1 Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, Janin J, Boer M. Leishmaniasis worldwide and global estimates of its incidence. PLoS One 2012; 7(5):e35671.^,1111 Organização Pan-Americana da Saúde (OPAS). Leishmanioses: Informe Epidemiológico das Américas [internet]. 2021. [cited 2022 maio 7]. Available from: https://iris.paho.org/handle/10665.2/55386
https://iris.paho.org/handle/10665.2/553... . This scenario reflects that VL is a relevant public health problem in the national territory and makes the country a priority area for measures to monitor and control the disease, either through epidemiological or environmental surveillance actions.

In 2020, the lowest number of VL occurrences of the time series was recorded in the national territory, a phenomenon that was significantly reflected in the incidence of this disease (n = 1,933; 0.91 cases/100 thousand inhabitants). This significant decrease in the incidence of VL cases may represent an indication of the reduction in the force of morbidity and in the magnitude of the endemicity of this zoonosis nationwide, formerly treated as a disease of increasing incidence¹⁵15 Gontijo CMF, Melo MN. Leishmaniose visceral no Brasil: quadro atual, desafios e perspectivas. Rev Bras Epidemiol 2004; 7(3):338-349.. And it may also be one of the effects generated by the COVID-19 pandemic¹⁶16 Miguel DC, Brioschi MBC, Rosa LB, Minori K, Grazzia N. The impact of COVID-19 on neglected parasitic diseases: what to expect? Trends Parasitol 2021; 37(8):694-697.^,1717 Bertollo DMB, Soares MMCN. Impacto da pandemia de COVID-19 nas ações de vigilância e controle da infecção por leishmaniose visceral. Brazilian J Infect Dis 2022; 26:102288..

In analyzing the data, it is observed that there was a formation of an inflection point (joinpoint) in the time series of the incidence rate for the national territory in 2018. This behavior was influenced by the reduction in the number of cases at the level of states and regions, especially in the Northeast, which concentrates the greatest burden of the disease in the country and which, according to the data series, presented stability in seven of its Federation Units, decrease in Ceará, and increase only in the state of Alagoas.

Two studies carried out in Ceará indicated heterogeneity of VL in that state. In the first one, VL presented a trend of increase and expansion for the Central-West and East regions between 2003-2017, and a decrease from 2016 on¹⁸18 Cavalcante FRA, Cavalcante KKS, Florencio CMGD, Moreno JO, Correia FGS, Alencar CH. Human visceral leishmaniasis: epidemiological, temporal and spatial aspects in Northeast Brazil, 2003-2017. Rev Inst Med Trop Sao Paulo 2020; 62:e12.. In the second one, which also used segmented time regression, the global incidence rate of VL showed an increasing trend in the capital between 2007-2017 and a decrease for the interval 2015-2017¹⁹19 Almeida CP, Cavalcante FRA, Moreno JO, Florêncio CMGD, Cavalcante KKS, Alencar CH. Visceral leishmaniasis: temporal and spatial distribution in Fortaleza, Ceará State, Brazil, 2007-2017. Epidemiol Serv Saude 2020; 29(5):e2019422.. These diverse findings possibly result from time frames and trend analysis adopted in the studies.

Another study, which had the state of Alagoas as its analysis space, confirms a progressive evolution of VL incidence rates in the Federation Unit, showing geographic expansion, a trend of zoonosis persistence in areas that already had reported cases, and overall strengthening of the endemic status of the disease in that territory between 2007 and 2018²⁰20 Braz BMDA, Silva RBS, Lins SC, Silva DRX, Ramalho WM, Melo MA. Demographic and spatial study of visceral leishmaniasis in the state of Alagoas, Brazil, during 2007-2018. Rev Soc Bras Med Trop 2021; 54:e06102020..

A historical series study of VL in the Northeast region developed by Batista et al.²¹21 Batista FMA, Sousa RA, Aguiar BGA, Ibiapina AB, Araújo Albuquerque LP, Mendonça VJ, Costa CHN. Perfil epidemiológico e tendência temporal da leishmaniose visceral: Piauí, Brasil, 2008 a 2018. Cad Saude Publica 2021; 37(11):e00340320. highlighted the cyclical nature of the disease. The authors identified epidemic peaks in 1983-1984, 1993-1994, 2003-2004, and 2013-2014 in Piauí and a general trend of stability for VL incidence, mortality, and lethality rates between 1971 to 2018 in the state. Another study, in agreement with our findings, also indicated a decreasing trend of VL incidence for Tocantins²²22 Rodrigues MGM, Viana JA, Pereira EG, Doutor B. Epidemiological analysis of cases of visceral and human tegumentary leishmaniosis in the state of Tocantins in the years of 2009 to 2019. Brazilian J Dev 2021; 7(9):87507-87528..

Some national and Latin American studies also point to this cyclical behavior and the heterogeneous nature of VL transmission, which has varied between high and low incidences¹⁹19 Almeida CP, Cavalcante FRA, Moreno JO, Florêncio CMGD, Cavalcante KKS, Alencar CH. Visceral leishmaniasis: temporal and spatial distribution in Fortaleza, Ceará State, Brazil, 2007-2017. Epidemiol Serv Saude 2020; 29(5):e2019422.^,2323 Góes MAO, Melo CM, Jeraldo VLS. Série temporal da leishmaniose visceral em aracaju, estado de sergipe, Brasil (1999 a 2008): aspectos humanos e caninos. Rev Bras Epidemiol 2012; 15(2):298-307.

24 Belo VS, Werneck GL, Barbosa DS, Simões TC, Nascimento BWL, Silva ES, Struchiner CJ. Factors associated with visceral leishmaniasis in the americas: a systematic review and meta-analysis. PLoS Negl Trop Dis 2013; 7(4):e2182.^-2525 Bermudi PMM, Guirado MM, Rodas LAC, Dibo MR, Chiaravalloti-Neto F. Spatio-temporal analysis of the occurrence of human visceral leishmaniasis in Araçatuba, State of São Paulo, Brazil. Rev Soc Bras Med Trop 2018; 51(4):452-460.. In a study developed by Servadio et al.²⁶26 Servadio JL, Machado G, Alvarez J, Lima FEF, Alves RV, Convertino M. Information differences across spatial resolutions and scales for disease surveillance and analysis: the case of visceral leishmaniasis in Brazil. PLoS One 2020; 15(7):e0235920., the asymmetry of VL distribution in the country is also revealed, with emphasis on the few states that followed the same trend of national incidence.

In this study, the trend for the Southeast regions and the state of São Paulo was decreasing. In agreement, a survey that evaluated the space-time distribution of VL in the western region of that state, from 2004 to 2018, also indicated a decreasing trend in the number of cases of the disease, especially after 2010²⁷27 Rancan EA, Chagas EFB, Sperança MA, Carvalho VCL, Martins LPA, Suzuki RB. Spatio-temporal distribution of human American visceral leishmaniasis in the Western region of Sao Paulo State, from 2004 to 2018. Rev Inst Med Trop Sao Paulo 2020; 62:e80.. The authors of the study suggest that this phenomenon may be related to VLSCP settlement in that region²⁷27 Rancan EA, Chagas EFB, Sperança MA, Carvalho VCL, Martins LPA, Suzuki RB. Spatio-temporal distribution of human American visceral leishmaniasis in the Western region of Sao Paulo State, from 2004 to 2018. Rev Inst Med Trop Sao Paulo 2020; 62:e80..

During the time series, the state of Tocantins concentrated the highest burden of the disease in the country. Reis et al.²⁸28 Reis LL, Silva Balieiro AA, Fonseca FR, Gonçalves MJF. Leishmaniose visceral e sua relação com fatores climáticos e ambientais no Estado do Tocantins, Brasil, 2007 a 2014. Cad Saude Publica 2019; 35(1):e00047018. point out that such behavior takes place mostly in the north and south regions of that state, being influenced by high values of annual precipitation and humidity, in addition to high levels of vegetation and night temperature - environmental conditions favorable vector reproduction and to wild reservoirs of VL and that end up favoring the increase of the disease in this territory.

The differences between the behavior of national and federation units’ incidence rates suggest that different factors can influence the risk of VL among states and regions of the country. Rocha et al.²⁹29 Rocha MAN, Matos-Rocha TJ, Ribeiro CMB, Abreu SRO. Epidemiological aspects of human and canine visceral leishmaniasis in State of Alagoas, Northeast, Brazil. Braz J Biol 2018; 78(4):609-614., Braz et al. ^20, and Ribeiro et al.³⁰30 Ribeiro CJN, Santos AD, Lima SVMA, Silva ER, Ribeiro BVS, Duque AM, Peixoto MVS, Santos PL, Oliveira IM, Lipscomb MW, Araújo KCGM, Moura TR. Space-time risk cluster of visceral leishmaniasis in Brazilian endemic region with high social vulnerability: an ecological time series study. PLoS Negl Trop Dis 2021; 15(1):e0009006. agree that the behavior of this zoonosis in the regions of the country is directly associated with the environmental and geographic conditions of the territory, where situations of high temperature and humidity predominate, which enable a good adaption and reproduction of the vector; and precarious socio-environmental issues such as lack of basic sanitation and water supply, poor housing conditions, accumulation of waste from the peridomicile both in an urban and rural area, the intense population of vectors and canine cases of the disease, presence of livestock close to residences and proximity of households to green areas, factors that directly interfere in the adaptation and reproduction of L. longipalpis.

It is worth noting in this area that the trends of incidence rates to remain stable, with no significant variations, in most federation units, may indicate the fragility of prevention and control measures, either due to lack of effectiveness or implementation difficulties in the VLSCP actions, or possible failures in the surveillance actions of this disease. This reality added to the high occurrence and burden of this disease still creates major challenges for health services, especially for the unified system, responsible for VL surveillance and control policies nationwide. On this topic, Werneck³¹31 Werneck GL. Controle da leishmaniose visceral no Brasil: o fim de um ciclo? Cad Saude Publica 2016; 32(6):S0102-311X2016000600201. highlights that despite the efforts and resources committed to the full functioning of VLSCP, the perception is consolidated that the actions aimed at reducing the transmission of VL have not had the desired effect of reducing the incidence of this disease.

Throughout the series analyzed, there was a higher incidence of cases among males and children under 5 years of age, with a record of stability for the occurrence in males and decreasing trend of cases in females and also in the group under 5 years of age. The specialized literature mentions that the greater the number of cases in a region, the greater the vulnerability of infected children, and when affecting the group of male adults or elderly⁶6 Safavi M, Eshaghi H, Hajihassani Z. Visceral Leishmaniasis: Kala-azar. Diagn Cytopathol 2021; 49(3):446-448.^,1010 Bezerra JMT, Araújo VEM, Barbosa DS, Martins-Melo FR, Werneck GL, Carneiro M. Burden of leishmaniasis in Brazil and federated units, 1990-2016: findings from Global Burden of Disease Study 2016. PLoS Negl Trop Dis 2018; 12(9):e0006697.^,3030 Ribeiro CJN, Santos AD, Lima SVMA, Silva ER, Ribeiro BVS, Duque AM, Peixoto MVS, Santos PL, Oliveira IM, Lipscomb MW, Araújo KCGM, Moura TR. Space-time risk cluster of visceral leishmaniasis in Brazilian endemic region with high social vulnerability: an ecological time series study. PLoS Negl Trop Dis 2021; 15(1):e0009006..

In time, it is pertinent to consider that this work has all the limitations inherent to the fact that it is an ecological study, namely: the possibility of ecological fallacy, chance correlations, and difficulties in controlling confounding factors. It also has limitations arising from the use of secondary data, which, as is known, are subject to underreporting, eventual errors due to incompleteness, precision or inconsistency, or even underdiagnosis, a common phenomenon across the country¹1 Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, Janin J, Boer M. Leishmaniasis worldwide and global estimates of its incidence. PLoS One 2012; 7(5):e35671..

Finally, it is important to have in mind that the data presented here, especially those concerning the last two years of the time series, may represent one of the indirect consequences generated by the COVID-19 pandemic¹⁶16 Miguel DC, Brioschi MBC, Rosa LB, Minori K, Grazzia N. The impact of COVID-19 on neglected parasitic diseases: what to expect? Trends Parasitol 2021; 37(8):694-697.^,1717 Bertollo DMB, Soares MMCN. Impacto da pandemia de COVID-19 nas ações de vigilância e controle da infecção por leishmaniose visceral. Brazilian J Infect Dis 2022; 26:102288.. Therefore, the number of cases and consequently, the VL incidence rate in the series analyzed may be higher than what was observed, wither because of the recognized impact of the SARS-CoV-2 public health emergency in the fight against neglected tropical diseases, including the VL, or throughout the series, due to deficiencies in health services, especially in locations where assistance and surveillance actions have operational problems linked to the chronic lack of inputs and financial, material, or personnel resources. It should be noted, however, that the results of this study were based on official sources of information, the same ones adopted for decision-making regarding public policies to combat and monitor VL nationwide.

Conclusion

Despite VL being an endemic disease in the national territory, there was a statistically significant decrease in the general incidence rate of this anthropozoonosis in Brazil for the period studied. It is important to consider, however, that this disease remains an important public health problem in the country and its regions, because, despite the verified trends of reduction in the general rate for the State, it is concomitantly observed inequalities in the epidemiological behavior of this indicator at the level of the political-administrative units, with an important portion of Federation units with incidence rates that remain on an increasing trend or in a stationary behavior. In this way, this study reaffirms the need to enhance effective measures to promote VL control in the country, with emphasis on the territories with the highest burden of this disease.

References

Publication Dates

History