Acessibilidade / Reportar erro

Visceral leishmaniasis lethality in Brazil: an exploratory analysis of associated demographic and socioeconomic factors

Abstract

INTRODUCTION:

It is believed that delays in diagnosis and treatment of Visceral Leishmaniasis (VL) contribute significantly to the burden of VL lethality in Brazil.

METHODS:

This study included several parts: a descriptive cross-sectional study of the individual characteristics of deaths from disease; a descriptive ecological study of the spatial distribution of deaths from disease; and an ecological analytical study to evaluate the association between disease lethality rates and the demographic, socioeconomic, and health indicators. The study population comprised all cases diagnosed throughout the country per the National Disease Notification System (SINAN) and the total number of disease deaths recorded in the Mortality Information System (SIM) from 2007 to 2012.

RESULTS:

Of the 223 deaths from disease captured by pairing the databases, 59.1% were reported as "death from other causes". There were significant associations between VL lethality rate and municipalities with the highest proportion of vulnerable individuals (rate ratio (RR)=1.18, 95% confidence interval (CI): 1.01-1.27), with VL lower incidence rate (RR=0.62, 95% CI: 0.58-0.67) and a higher incidence rate of Acquired Immune Deficiency Syndrome (AIDS) (RR=1.20, 95% CI: 1.17-1.51).

CONCLUSIONS:

Linking the SINAN and SIM databases allowed the inclusion of 14% of otherwise underreported deaths from VL for the study period, showing that this method is useful for the surveillance of VL-related deaths. The size of the municipal population, proportion of the vulnerable population, incidence of disease, and the incidence of AIDS were associated with municipal lethality rates related to VL in Brazil.

Keywords:
Visceral leishmaniasis; Acquired Immune Deficiency Syndrome; Mortality

INTRODUCTION

Visceral leishmaniasis (VL) is an infectious disease caused by protozoans of the genus Leishmania. In Brazil, VL is associated with infection by Leishmania infantum (syn. L. chagasi) infections, and presents through a zoonotic cycle in which the domestic dog constitutes the main reservoir of the parasite11. Belo VS, Werneck GL, Barbosa DS, Simoes TC, Nascimento BW,da 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):1-12.. The vectors responsible for transmission of the disease in the country are Lutzomyia longipalpis and Lutzomyia cruzi22. Missawa NA, Veloso MAE, Maciel GBML, Michalsky EM, Dias ES. Evidência de transmissão de leishmaniose visceral porLutzomyiacruzino município de Jaciara, Estado de Mato Grosso. Rev Soc Bras Med Trop. 2012;44(1):76-8.. About 90% of new cases are estimated to be concentrated in seven countries, and in 2017 a total of 20,792 cases were reported in these countries33. World Health Organization. Leishmaniasis - Epidemiological Situation, 2019. (Update 2019 Jan 15; cited 2019 Set 10). Alvaiable from: Alvaiable from: https://www.who.int/leishmaniasis/burden/en/
https://www.who.int/leishmaniasis/burden...
. In the Americas, most cases of VL occur in 12 Latin American countries, and 90% of these occur in Brazil44. Barbosa IR, Costa ÍCC (2013). Aspectos clínicos e epidemiológicos da leishmaniose visceral em menores de 15 anos no estado do Rio Grande do Norte, Brasil. Sci Med (Porto Alegre). 2013. 23(1): 5-11..

The Brazilian Ministry of Health has identified reduction of VL lethality as the goal of the Visceral Leishmaniasis Surveillance and Control Program (Programa de Vigilância e Controle da Leishmaniose Visceral (PVC-VL))55. Ministério da Saúde, Manual de Leishmaniose Visceral: recomendações clínicas para redução da letalidade. 1st ed. Brasil: Ministério da Saúde; 2011.78p.. The guidelines described to achieve this goal include early diagnosis and treatment, serological surveys in dogs, vector control, and health education66. Coura-Vital, WA, VEM, Reis, IA, Amancio, FF, Reis, AB, Carneiro M. Prognostic Factors and scoring system for death from visceral leishmaniasis: An historical cohort study in Brazil. PLoS Negl Trop Dis . 2014;8(12):1-12.. In fact, it is believed that late diagnosis and treatment of the disease contribute significantly to the burden of VL lethality in Brazil77. Romero GAS, Boelaert M. Control of visceral leishmaniasis in Latin America a systematic review. PLoS Negl Trop Dis . 2010;4(1):1-17..

In recent years, there has been some progress in the search for solutions to combat VL, especially with regard to diagnosis and treatment. However, few studies have addressed issues related to prognosis, with disease lethality as the outcome88. Belo VS, Struchiner CJ, Barbosa DS, Nascimento BW, Horta MA, da Silva ES, Werneck GL. Risk factors for adverse prognosis and death in American visceral leishmaniasis: a meta-analysis. PLoS Negl Trop Dis . 2014;8(7):1-9..

The fatality rate for VL in Brazil has increased gradually in recent years, from 3.6% of diagnosed cases in 1994 to 8.5% in 2003. This rate decreased between 2003 and 2007, but starting from 2008, the lethality rate again increased to reach 7.1% in 201299. Alvarenga DG, Escalda PMF, Costa ASVM, Duenhas MTF. Leishmaniose visceral: estudo retrospectivo de fatores associados à letalidade. Rev Soc Bras Med Trop . 2010;43(2):194-7.. In Campo Grande (the capital of the state of Mato Grosso do Sul), where diagnostic tracking of human cases of VL began in 2001, the mortality rate reached 7% with a specific mortality rate of 9.1 per 100,000 inhabitants. These statistics demonstrate the severity of the problem1010. Oliveira JM, Fernandes AC, Dorval MEC, Alves TP, Fernandes TD, Oshiro ET. Mortalidade por leishmaniose visceral: aspectos clínicos e laboratoriais. Rev Soc Bras Med Trop . 2010;43(2):188-93.. A more critical situation was observed in the capital city of Belo Horizonte, which experienced a lethality rate of 23.6% in 2009, surpassing the national average of the same year1111. 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):1-9..

In addition, VL lethality observed in people living with Human Immunodeficiency Virus (HIV)/Acquired Immune Deficiency Syndrome (AIDS) (PLHA) has been receiving increasing attention1212. Lindoso JA, Cota GF, da Cruz AM, Goto H, Maia-Elkhoury ANS, Romero GAS, et al. Visceral leishmaniasis and HIV coinfection in Latin America. PLoS Negl Trop Dis . 2014;8(9):e3136.. The estimated lethality rate of VL/AIDS co-infection in Brazil between 2001 and 2010 was approximately 25%, which corresponds to three times the lethality rate observed in non-coinfected individuals1313. Leite de Sousa-Gomes M, Romero GAS, Werneck GL. Visceral leishmaniasis and HIV/AIDS in Brazil: Are we aware enough? PLoS Negl Trop Dis . 2017;11(9):1-13..

Because it is a neglected disease which occurs mainly in tropical countries, the burden of VL falls disproportionately on the poorest segments of the global population. Lack of access to healthcare is associated with late diagnosis and treatment, and with increased VL morbidity and mortality1414. Alvar J, Yactayo S, Bern C. Leishmaniasis and poverty. Trends Parasitol. 2006;22(12):552-7.. Public investment in care activities, combined with the most effective surveillance tools, could contribute to reducing the burden and magnitude of VL44. Barbosa IR, Costa ÍCC (2013). Aspectos clínicos e epidemiológicos da leishmaniose visceral em menores de 15 anos no estado do Rio Grande do Norte, Brasil. Sci Med (Porto Alegre). 2013. 23(1): 5-11..

In view of the above, it is believed that increasing knowledge of the factors associated with VL lethality can help improve the currently recommended public health policies, making it possible to increase the number of epidemiological surveillance programs in municipalities and states, thereby avoiding VL-related deaths. Therefore, the objective of this study was to identify socioeconomic, demographic, and health factors of the municipalities associated with lethality due to VL in Brazil from 2007 to 2012.

METHODS

Study design

The study was divided into three stages: a descriptive cross-sectional study of the individual characteristics of deaths from VL; a descriptive ecological study of the spatial distribution of deaths from VL; and an ecological analytical study to evaluate the association between VL lethality rates and the demographic, socioeconomic, and health indicators of Brazilian municipalities.

Study scope and population

The study population comprised all VL cases diagnosed throughout the country per the National Disease Notification System (SINAN) and the total number of VL deaths recorded in the Mortality Information System (SIM) from 2007 to 2012. For the ecological approach, all 5,569 Brazilian municipalities were included.

Definition of incident case and death from VL

A case of VL was considered to be any case recorded in SINAN for which the variable "final classification" was classified as "confirmed."

Deaths were identified in two ways:

1) cases of VL recorded in SINAN for which the variable “progression” was classified as "death from VL"; and

2) cases of VL recorded in SINAN (without information on progression to death) paired with notification in the SIM during the same period, with the underlying cause of death in the latter system being VL, bleeding, sepsis, or HIV. For the HIV variable we chose the ecological approach considering that the status of this variable was not available for a relevant number of cases in the SINAN database.

To link the SINAN and SIM databases for the identification of VL deaths not recorded in SINAN, probabilistic pairing was applied using Python Brasil with a bloom filter1515. Rothenberg C, Macapuna C, Verdi F, Magalhães M, Zahemszky A. Data center networking with in-packet bloom filters. Computer Networks. 2010;55(6):1364-78.. This program automatically identifies records while considering parameters which are probabilistically estimated. The following variables were used in the comparison fields: patient name, mother's name, and date of birth.

After identifying the records, the notifications which could not be submitted to the pairing process between the banks because they exhibited some inconsistencies among the comparison variables were evaluated by three observers. For these cases, the following pre-established pairing criteria were used: date of being reported into SINAN; date of birth of the individual and the mother; date of death in the SINAN and SIM databases; and municipality of residence. After this evaluation, cases with consistency in the variable data defined above were considered to be deaths.

Spatial distribution of deaths from VL

For the spatial distribution of deaths from VL, the municipality of residence identified in the notification records was considered as an aggregated level of the spatial units. Thematic maps were created using the TerraView 4.2.2 software, and the municipal layers were obtained from the censuses of the Brazilian Institute of Geography and Statistics (IBGE).

Socioeconomic, demographic, and health variables of municipalities

The variables described in Box 1 were explored. Socioeconomic variables were obtained from the Atlas of Human Development in Brazil, published in 2013 by the United Nations Development Program1616. Programa das Nações Unidas. Atlas do Desenvolvimento Humano dos Municípios. Brasil; 2012 (update2014 April 22; cited 2014 Jul 4). Avaiable from Avaiable from http://www.pnud.org.br/IDH/Default.aspx?indiceAccordion=1&li=li_AtlasMunicipios
http://www.pnud.org.br/IDH/Default.aspx?...
. This atlas sources data and indicators from the last three demographic censuses conducted by the Brazilian Institute of Geography and Statistics1717. Brazilian Institute of Geography and Statistics. Estimativas populacionais para os municípios brasileiros. 2013. (Update 2014 May15;cited 2014 Jul 4). Available from: Available from: https://ww2.ibge.gov.br/home/estatistica/populacao/estimativa2013/default.shtm
https://ww2.ibge.gov.br/home/estatistica...
. Population data were obtained from the Brazilian Institute of Geography and Statistics1717. Brazilian Institute of Geography and Statistics. Estimativas populacionais para os municípios brasileiros. 2013. (Update 2014 May15;cited 2014 Jul 4). Available from: Available from: https://ww2.ibge.gov.br/home/estatistica/populacao/estimativa2013/default.shtm
https://ww2.ibge.gov.br/home/estatistica...
.

Healthcare-related variables were obtained from the National Registry of Health Establishments (CNES) and the Department of Primary Health Care (DAB) of the Secretariat of Health Care (SAS), Ministry of Health. Disease-related variables were obtained from SINAN and the Information Technology Department of the Unified Health System1818. Datasus Ministério da Saúde. Informações de saúde. Brasil; 2014 [updated 2014 May 15; cited 2014 Jul 8]. Alvaiable from: Alvaiable from: http://www2.datasus.gov.br/DATASUS/index.php?area=01
http://www2.datasus.gov.br/DATASUS/index...
.

Statistical analysis

In the descriptive step of data analysis, the absolute and relative frequencies of the variables characterizing the cases of and deaths from VL in the evaluated period were calculated. Maps were created for the analysis of the spatial distribution of cases and deaths from VL, and the municipality of residence of the notification records was used for spatial unit aggregation.

In the analytical stage, the ecological associations between the characteristics of the Brazilian municipalities and the VL lethality rate were evaluated. For this purpose, the individual effects of categorical and continuous independent variables associated with the lethality rate were evaluated by a zero-inflated Poisson model with a link function log (λ)1919. Cameron AC, PK Trivedi. Regression Analysis of Count Data. 2nd ed. Cambridge University; 2013. 566 p.,2020. Khan A, Ullah S, Nitz J. Statistical modelling of falls count data with excess zeros. Inj Prev. 2011;17(4):266-70., using the statistical software Stata® 11.2. The estimates adjusted using the zero-inflated Poisson model were obtained as follows. A Pearson correlation matrix was first developed to identify collinearity between the independent variables of interest. Variables with correlation coefficients r≥0.8 were considered collinear, and only one of the variables was selected for the model. This choice was based on parameters such as the ease of measurement of the variable and its potential for being understood by decision-makers. Next, some continuous variables were categorized based on the behavior observed in the preliminary bivariate analysis.

Finally, multivariate models were constructed. Variables exhibiting a statistically significant association with the lethality rate, defined by a p-value < 0.20, were included. Multivariate analysis was conducted using a hierarchical approach with the model incorporating two levels: the distal, including demographic variables; and the proximal, including health indicators. The backward method was used on each of the hierarchical levels to select the independent variables.

Distal variables were added first to the multivariable model, and retained as long as they remained significantly (P<0.05).

Ethical aspects

The project safeguarded the privacy of the included participants by maintaining confidentiality of patient names in VL cases reported in the SIM and SINAN databases. All ethical precepts of human research were implemented. The project was approved by the Research Ethics Committee of the School of Health Sciences of the University of Brasília (CEP/FS-UnB) on March 13, 2014 (CAAE n. 22229613.5.0000.0030). The research was conducted in accordance with the principles of Declaration of Helsinki.

RESULTS

In the period from 2007 to 2012, 1,368 deaths from VL in Brazil were notified in SINAN. After pairing the SINAN and SIM databases, 223 additional deaths were identified which were not indicated as “death from VL” cases in the “progression” field of the SINAN notification form. Thus, the total number of deaths identified for the period studied was 1,591. In 2010, the highest number of deaths from VL (n=48, 21.5%) was recorded, captured by pairing the SINAN and SIM databases (Table 1). Of the 223 deaths from VL captured by pairing the databases, 19.7% (44/223) had "cure" recorded in the “progression” field of the SINAN form; 59.1% (132/223) had "death from other causes" recorded; 14.7% (33/223) reported "not [being] informed"; and 6.2% reported "transferred" (14/223).

TABLE 1:
Distribution of visceral leishmaniasis deaths identified in the National Disease Notification System (SINAN) and through the linkage of this database with the Mortality Information System (SIM) - Brazil, 2007 to 2012.

Figure 1 shows the distribution of total deaths from VL (n=1,591) identified in Brazil between 2007 and 2012. It is possible to observe a concentration of deaths from VL in the Northeast region, though there is a distribution of VL-related deaths throughout the entire country (Figure 1).

FIGURE 1:
Distribution of deaths from Visceral Leishmaniasis in 2007 to 2012.

Deaths from VL in Brazil predominantly affected men (66.8%) during the period from 2007 to 2012. Analysis by age group in the same period showed a predominance of adults aged 35 to 49 years, representing 21% of deaths (Table 2). The accumulated deaths in this period were 1,591, with an annual average of 260 deaths; the year 2007 had the lowest recorded number of deaths (n=234), whereas 2011 showed the highest number of recorded deaths (n=284).

The variables selected to evaluate association with the VL lethality rate in Brazilian municipalities during the 2007-2012 period are described in Table 3.

TABLE 2:
Proportional distribution of deaths from visceral leishmaniasis according to age group. Brazil, 2007-2012.
TABLE 3:
Description of the variables selected to evaluate association with the lethality rate from visceral leishmaniasis in Brazilian municipalities during the 2007-2012 perQiod.

In the crude (bivariate) analysis, the VL lethality rate was significantly associated with small municipalities (municipal population >20,000 and <100,000 inhabitants (p<0.001)) including a higher proportion of vulnerable population (p<0.001), with a lower incidence rate of VL per 100,000 inhabitants (p<0.001), and with a higher incidence of AIDS per 100,000 inhabitants (p<0.001). In the adjusted analysis, small municipalities (population 20,000 to <50,000) were associated with higher rates of VL lethality (rate ratio (RR)=1.90; 95%CI: 1.68-2.15 for population size 20,000 to <50,000 and RR=1.43, 95%CI: 1.29-1.60 for population size 50,000 to <100,000) compared to municipalities with a population of <20,000 inhabitants. There were also significant associations between the VL lethality rate and municipalities with the highest proportion of vulnerable individuals (RR=1.18, 95%CI: 1.01-1.27), with a lower incidence rate of VL (RR=0.62, 95%CI: 0.58-0.67), and a higher incidence rate of AIDS per 100,000 inhabitants (RR=1.20, 95% CI: 1.17-1.51) (Table 4).

TABLE 4:
Crude and adjusted analysis of factors associated with the rate of lethality due to visceral leishmaniasis in Brazil, 2007 to 2012.

DISCUSSION

Reducing the VL lethality rate is currently the main goal of the control program designed by the Brazilian Ministry of Health2121. Ministério da Saúde. Manual de Vigilância e Controle da Leishmaniose Visceral. 1st ed. Brasil; 2003.120 p.. The lethality of VL can be explained by a range of factors, such as late diagnosis and treatment, toxicity of drugs available for treatment, comorbidities, poor healthcare quality, and socioeconomic factors related to the individual and the environment2222. Madalosso G, Fortaleza CM, Ribeiro AF, Cruz LL, Nogueira PA, Lindoso JAL. American Visceral Leishmaniasis: Factors Associated with Lethality in the State of São Paulo, Brazil. J Trop Med J. 2012;2012:281572.. Timely and accurate diagnosis remains a challenge in Brazil, as in other affected countries, where the disease is still treated only on the basis of clinical suspicion. The use of tests with good accuracy which are easy to perform, such as rapid tests employing recombinant proteins, has opened a space for improving the decision-making process for health professionals, regarding the provision of timely treatment2323. Machado de Assis TS, Rabello A, Werneck GL. Predictive models for the diagnostic of human visceral leishmaniasis in Brazil. PLoS Negl Trop Dis . 2012;6(2):1-7.,11. Belo VS, Werneck GL, Barbosa DS, Simoes TC, Nascimento BW,da 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):1-12.,2424. Ritmeijer K, Melaku Y, Mueller M, Kipngetich S, O’keeffe C, Davidson RN. Evaluation of a new recombinant K39 rapid diagnostic test for Sudanese visceral leishmaniasis. Am J Trop Med Hyg. 2006;74(1):76-80..

The results of the present study point to a higher lethality rate in small municipalities (20,000 to fewer than 100,000 inhabitants). This result could be explained in part by the instability of lethality rates in municipalities of this size, where few cases are diagnosed. In the literature, it is possible to observe an association between lethality and larger municipalities which have a more complex healthcare structure and attract more severe cases, thus selecting for individuals with higher risk of death2525. Queiroz MJA, Alves JGB, Correia JB. Leishmaniose visceral: características clínico-epidemiológicas em crianças de área endêmica. J Pediatr (Rio J). 2004;80(2): 141-6.. In these cases, the place of death and not the place of residence of the patient is considered in the analysis/notification of the event. This phenomenon could explain the lower lethality rate in municipalities with fewer than 20,000 inhabitants relative to municipalities with 20,000 to 100,000 inhabitants, which would most likely select for patients with less severe VL. However, the trend in lethality rates in municipalities with more than 100,000 inhabitants still remains to be explained. Because the disease generally has a subacute or chronic course, even considering that the quality of and access to healthcare is precarious in municipalities with fewer than 20,000 inhabitants, there would be an opportunity to seek care in larger and denser neighboring municipalities. This may improve the lethality rate of the smaller municipalities because the deaths would then be recorded in the municipality where they actually occurred and not in the municipality where the patient resided. Other factors may include migration of rural populations to the periphery of large urban centers with more than 100,000 inhabitants, coupled with poor housing, uncontrolled deforestation, and a high number of infected dogs, which would contribute to the increase in VL cases and consequently to death from the disease2626. Martins-Melo FR, Lima MS, Ramos AN, Alencar CH, Heukelbach J. Mortality and case fatality due to visceral leishmaniasis in Brazil: A nationwide analysis of epidemiology, trends and spatial patterns. PLoS One. 2014;9(4):1-14.. In contrast, the spread of the zoonosis to areas with a large population density could negatively affect lethality, as a result of the lack of knowledge about the disease, both among the most vulnerable population and among the healthcare teams responsible for providing services in these places. In this scenario, the development of a healthcare network with qualified referral centers which can quickly accumulate the experience necessary for proper management of the most serious cases and associated complications, becomes essential.

The observed historical trend suggests the existence of a population vulnerable to the disease in the age groups corresponding to the extremes of life: individuals younger than 1 year old and those older than 50 years. In the present study, it was observed that on average, one-fifth of the population studied in this period met this profile, and an association of lethality with the proportion of vulnerable population was found. Factors potentially related to this greater vulnerability of individuals in these age groups may be the time between the onset of first symptoms and physician consultation, malnutrition, immune problems, and the presence of other associated pathologies which may play an important role in the genesis of an unfavorable outcome99. Alvarenga DG, Escalda PMF, Costa ASVM, Duenhas MTF. Leishmaniose visceral: estudo retrospectivo de fatores associados à letalidade. Rev Soc Bras Med Trop . 2010;43(2):194-7.,2727. 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.,2828. Costa DL, Rocha RL, Chaves EBF, Batista VGV, Costa HL, Costa CHN. Predicting death from kala-azar: construction, development, and validation of a score set and accompanying software. Rev Soc Bras Med Trop . 2016;49(6):728-74.,2929. Fu Q, Li SZ, Wu WP, Hou YY, Zhang S, Feng Y, et al. Endemic characteristics of infantile visceral leishmaniasis in the People’s Republic of China. Parasit Vectors. 2013;6(143):2-8..

The ratio of doctors to number of inhabitants and coverage by the family healthcare plan was not significantly associated with VL lethality in this study. In a study which evaluated the determinants of the causes of mortality, a negative association between infant mortality and the number of doctors was identified3030. Boing AF, Boing AC. Mortalidade infantil por causas evitáveis no Brasil: um estudo ecológico no período 2000-2002. Cad Saude Publica. 2008;24(2):447-55.. It is important to note that in the present analysis, all medical specialties were included, and a large proportion of these professionals do not participate directly in the diagnosis and treatment of VL, which may explain the observed result.

Extreme poverty was not significantly associated with VL in bivariate analysis; however, some studies have found that income may be associated with VL. We found that in 47% of the interviewed cases, the income was below the minimum wage, set at BRL120.00 in 1997. In another study, it was identified that increased income may be associated with a decrease in the occurrence of VL11. Belo VS, Werneck GL, Barbosa DS, Simoes TC, Nascimento BW,da 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):1-12.,3131. Cesse EÂP, Carvalho EF, Andrade PP, Ramalho WM, Luna L. Organização do espaço urbano e expansão do calazar. Rev Bras Saude Mater Infant. 2001;1(2):167-76..

None of the indicators directly related to poverty was associated with VL lethality in the period studied. Given that VL has been classically associated with poverty, the observed result may be explained by the municipalities which present with significant extreme poverty and a low human development index (HDI), combined with a long history of the disease in the community. In this scenario, despite the difficulties faced by the local health services, there would be accumulation of knowledge and skills for the management of patients with VL, including the timely identification of the need for referral of the most severe patients to referral centers with greater technological density, which could reduce disease lethality. As a counterpoint, there would be municipalities with a lower proportion of population in extreme poverty and higher HDI, where the disease was introduced recently; such conditions may worsen the quality of case management due to the inexperience of health teams, resulting in a negative impact on lethality. We identified areas at high risk of VL based on socioeconomic indicators and remote sensing data, and concluded that even with the complexity involved in this approach, it was possible to observe a relationship between disease spread and socioeconomic and environmental factors3232. Almeida AS, Werneck GL. Prediction of high-risk areas for visceral leishmaniasis using socioeconomic indicators and remote sensing data. Int J Health Geogr. 2014. 13(13):3-7..

The incidence of AIDS was associated with a lethal outcome. Carvalho and Cols. (2013) found an association between VL/HIV coinfection and death3333. Araújo VEM, Morais MHF, Reis IA, Rabello A, Carneiro M. Early Clinical Manifestations Associated with Death from Visceral Leishmaniasis. PLoS Negl Trop Dis . 2012;6(2):1-9.. The epidemiological profiles of AIDS and VL have changed in Brazil; data suggest that the interiorization of HIV infection simultaneously with the urbanization of VL may be responsible for the greater exposure of the population to both infections3434. Carvalho FL, Aires DLS, Segunda ZF, Azevedo MPSC, Silva de RC, Corrêa RGF, et al. Perfil epidemiológico dos indivíduos HIV positivo e coinfecção HIV-LV em um serviço de referência em São Luís, MA, Brasil. Cien Saude Colet. 2013;18(5):1305-12.. It was also reported that the municipalities with the highest number of notifications of VL/HIV coinfection cases were among the 20 municipalities with the highest incidence of AIDS3535. Sousa-Gomes ML, Maia-Elkhoury AN, Pelissari DM, Junior FEFL, Sena JM, Cechinel MP. Coinfecção Leishmania-HIV no Brasil: aspectos epidemiológicos, clínicos e laboratoriais. Epidemiol Serv Saude. 2011;20(4):519-26..

In addition to the pathophysiological factors of the disease which may contribute to an unfavorable outcome in co-infected patients, the use of the recommended medications for the treatment of VL may also contribute to death. The deleterious effect of the use of pentavalent antimonials on PLHA is already known. A cohort study reported that HIV infection is an independent predictor of poor outcome at six months after VL diagnosis3636. Cota GF, de Sousa MR, de Mendonça AL, Patrocinio A, Assunção LS, de Faria SR, et al. Leishmania-HIV co-infection: clinical presentation and outcomes in an urban area in Brazil. PLoS Negl Trop Dis . 2014;8(4):1-7..

The methodology used in the present study has the limitations intrinsic to ecological studies whose conclusions are limited to the aggregated data and do not allow the establishment of causal relationships at the individual level. Our choice of using the HIV incidence rate instead of the individual HIV status data because of the unavailability of these data in the SINAN records is an example of that limitation. Another relevant limitation refers to the use of secondary data, which may not accurately represent real life conditions. In the present study, the linkage between the databases was not directed toward detecting the underreporting of VL cases but, rather, to correct classification errors of VL cases recorded as survivors when they had actually died of VL. In this sense, the denominator of the number of VL cases was limited to SINAN records, which may have been affected by underreporting. The socioeconomic and demographic variables used were not specifically obtained to try to understand the process of illness and death by VL; therefore, they could have suffered from measurement bias because they were primarily collected for other purposes.

In conclusion, in this study, the linkage of the SINAN and SIM databases allowed for the inclusion of 14% of otherwise underreported deaths from VL for the study period, showing that it is a useful tool for the surveillance of this event. It was found that the size of the municipal population, the proportion of the vulnerable population, the incidence of VL, and the incidence of AIDS were associated with municipal lethality rates from VL in Brazil during the study period. These variables may help healthcare managers intervene to reduce the vulnerability of municipalities with these characteristics, to the fatal outcomes of VL.

REFERENCES

  • 1
    Belo VS, Werneck GL, Barbosa DS, Simoes TC, Nascimento BW,da 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):1-12.
  • 2
    Missawa NA, Veloso MAE, Maciel GBML, Michalsky EM, Dias ES. Evidência de transmissão de leishmaniose visceral porLutzomyiacruzino município de Jaciara, Estado de Mato Grosso. Rev Soc Bras Med Trop. 2012;44(1):76-8.
  • 3
    World Health Organization. Leishmaniasis - Epidemiological Situation, 2019. (Update 2019 Jan 15; cited 2019 Set 10). Alvaiable from: Alvaiable from: https://www.who.int/leishmaniasis/burden/en/
    » https://www.who.int/leishmaniasis/burden/en/
  • 4
    Barbosa IR, Costa ÍCC (2013). Aspectos clínicos e epidemiológicos da leishmaniose visceral em menores de 15 anos no estado do Rio Grande do Norte, Brasil. Sci Med (Porto Alegre). 2013. 23(1): 5-11.
  • 5
    Ministério da Saúde, Manual de Leishmaniose Visceral: recomendações clínicas para redução da letalidade. 1st ed. Brasil: Ministério da Saúde; 2011.78p.
  • 6
    Coura-Vital, WA, VEM, Reis, IA, Amancio, FF, Reis, AB, Carneiro M. Prognostic Factors and scoring system for death from visceral leishmaniasis: An historical cohort study in Brazil. PLoS Negl Trop Dis . 2014;8(12):1-12.
  • 7
    Romero GAS, Boelaert M. Control of visceral leishmaniasis in Latin America a systematic review. PLoS Negl Trop Dis . 2010;4(1):1-17.
  • 8
    Belo VS, Struchiner CJ, Barbosa DS, Nascimento BW, Horta MA, da Silva ES, Werneck GL. Risk factors for adverse prognosis and death in American visceral leishmaniasis: a meta-analysis. PLoS Negl Trop Dis . 2014;8(7):1-9.
  • 9
    Alvarenga DG, Escalda PMF, Costa ASVM, Duenhas MTF. Leishmaniose visceral: estudo retrospectivo de fatores associados à letalidade. Rev Soc Bras Med Trop . 2010;43(2):194-7.
  • 10
    Oliveira JM, Fernandes AC, Dorval MEC, Alves TP, Fernandes TD, Oshiro ET. Mortalidade por leishmaniose visceral: aspectos clínicos e laboratoriais. Rev Soc Bras Med Trop . 2010;43(2):188-93.
  • 11
    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):1-9.
  • 12
    Lindoso JA, Cota GF, da Cruz AM, Goto H, Maia-Elkhoury ANS, Romero GAS, et al. Visceral leishmaniasis and HIV coinfection in Latin America. PLoS Negl Trop Dis . 2014;8(9):e3136.
  • 13
    Leite de Sousa-Gomes M, Romero GAS, Werneck GL. Visceral leishmaniasis and HIV/AIDS in Brazil: Are we aware enough? PLoS Negl Trop Dis . 2017;11(9):1-13.
  • 14
    Alvar J, Yactayo S, Bern C. Leishmaniasis and poverty. Trends Parasitol. 2006;22(12):552-7.
  • 15
    Rothenberg C, Macapuna C, Verdi F, Magalhães M, Zahemszky A. Data center networking with in-packet bloom filters. Computer Networks. 2010;55(6):1364-78.
  • 16
    Programa das Nações Unidas. Atlas do Desenvolvimento Humano dos Municípios. Brasil; 2012 (update2014 April 22; cited 2014 Jul 4). Avaiable from Avaiable from http://www.pnud.org.br/IDH/Default.aspx?indiceAccordion=1&li=li_AtlasMunicipios
    » http://www.pnud.org.br/IDH/Default.aspx?indiceAccordion=1&li=li_AtlasMunicipios
  • 17
    Brazilian Institute of Geography and Statistics. Estimativas populacionais para os municípios brasileiros. 2013. (Update 2014 May15;cited 2014 Jul 4). Available from: Available from: https://ww2.ibge.gov.br/home/estatistica/populacao/estimativa2013/default.shtm
    » https://ww2.ibge.gov.br/home/estatistica/populacao/estimativa2013/default.shtm
  • 18
    Datasus Ministério da Saúde. Informações de saúde. Brasil; 2014 [updated 2014 May 15; cited 2014 Jul 8]. Alvaiable from: Alvaiable from: http://www2.datasus.gov.br/DATASUS/index.php?area=01
    » http://www2.datasus.gov.br/DATASUS/index.php?area=01
  • 19
    Cameron AC, PK Trivedi. Regression Analysis of Count Data. 2nd ed. Cambridge University; 2013. 566 p.
  • 20
    Khan A, Ullah S, Nitz J. Statistical modelling of falls count data with excess zeros. Inj Prev. 2011;17(4):266-70.
  • 21
    Ministério da Saúde. Manual de Vigilância e Controle da Leishmaniose Visceral. 1st ed. Brasil; 2003.120 p.
  • 22
    Madalosso G, Fortaleza CM, Ribeiro AF, Cruz LL, Nogueira PA, Lindoso JAL. American Visceral Leishmaniasis: Factors Associated with Lethality in the State of São Paulo, Brazil. J Trop Med J. 2012;2012:281572.
  • 23
    Machado de Assis TS, Rabello A, Werneck GL. Predictive models for the diagnostic of human visceral leishmaniasis in Brazil. PLoS Negl Trop Dis . 2012;6(2):1-7.
  • 24
    Ritmeijer K, Melaku Y, Mueller M, Kipngetich S, O’keeffe C, Davidson RN. Evaluation of a new recombinant K39 rapid diagnostic test for Sudanese visceral leishmaniasis. Am J Trop Med Hyg. 2006;74(1):76-80.
  • 25
    Queiroz MJA, Alves JGB, Correia JB. Leishmaniose visceral: características clínico-epidemiológicas em crianças de área endêmica. J Pediatr (Rio J). 2004;80(2): 141-6.
  • 26
    Martins-Melo FR, Lima MS, Ramos AN, Alencar CH, Heukelbach J. Mortality and case fatality due to visceral leishmaniasis in Brazil: A nationwide analysis of epidemiology, trends and spatial patterns. PLoS One. 2014;9(4):1-14.
  • 27
    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.
  • 28
    Costa DL, Rocha RL, Chaves EBF, Batista VGV, Costa HL, Costa CHN. Predicting death from kala-azar: construction, development, and validation of a score set and accompanying software. Rev Soc Bras Med Trop . 2016;49(6):728-74.
  • 29
    Fu Q, Li SZ, Wu WP, Hou YY, Zhang S, Feng Y, et al. Endemic characteristics of infantile visceral leishmaniasis in the People’s Republic of China. Parasit Vectors. 2013;6(143):2-8.
  • 30
    Boing AF, Boing AC. Mortalidade infantil por causas evitáveis no Brasil: um estudo ecológico no período 2000-2002. Cad Saude Publica. 2008;24(2):447-55.
  • 31
    Cesse EÂP, Carvalho EF, Andrade PP, Ramalho WM, Luna L. Organização do espaço urbano e expansão do calazar. Rev Bras Saude Mater Infant. 2001;1(2):167-76.
  • 32
    Almeida AS, Werneck GL. Prediction of high-risk areas for visceral leishmaniasis using socioeconomic indicators and remote sensing data. Int J Health Geogr. 2014. 13(13):3-7.
  • 33
    Araújo VEM, Morais MHF, Reis IA, Rabello A, Carneiro M. Early Clinical Manifestations Associated with Death from Visceral Leishmaniasis. PLoS Negl Trop Dis . 2012;6(2):1-9.
  • 34
    Carvalho FL, Aires DLS, Segunda ZF, Azevedo MPSC, Silva de RC, Corrêa RGF, et al. Perfil epidemiológico dos indivíduos HIV positivo e coinfecção HIV-LV em um serviço de referência em São Luís, MA, Brasil. Cien Saude Colet. 2013;18(5):1305-12.
  • 35
    Sousa-Gomes ML, Maia-Elkhoury AN, Pelissari DM, Junior FEFL, Sena JM, Cechinel MP. Coinfecção Leishmania-HIV no Brasil: aspectos epidemiológicos, clínicos e laboratoriais. Epidemiol Serv Saude. 2011;20(4):519-26.
  • 36
    Cota GF, de Sousa MR, de Mendonça AL, Patrocinio A, Assunção LS, de Faria SR, et al. Leishmania-HIV co-infection: clinical presentation and outcomes in an urban area in Brazil. PLoS Negl Trop Dis . 2014;8(4):1-7.

Publication Dates

  • Publication in this collection
    11 Sept 2020
  • Date of issue
    2020

History

  • Received
    07 May 2019
  • Accepted
    17 June 2020
Sociedade Brasileira de Medicina Tropical - SBMT Caixa Postal 118, 38001-970 Uberaba MG Brazil, Tel.: +55 34 3318-5255 / +55 34 3318-5636/ +55 34 3318-5287, http://rsbmt.org.br/ - Uberaba - MG - Brazil
E-mail: rsbmt@uftm.edu.br