INTRODUCTION
Visceral leishmaniasis (VL) is classified as a Neglected Tropical Climate Disease (NTD)1, characterized by acute, subacute or chronic evolution, with more than 90% of untreated cases leading to death2. Among NTDs, VL placed third in terms of mortalities in Brazil from 2000 to 2011, causing more deaths than dengue or malaria3. VL has the highest years of lost life (YLL) due to premature death, as well as a high mortality rate among NTD, according to studies of global burden of diseases4,5.
Visceral leishmaniasis is present on five continents and more than 90% of cases occur in only 6 countries: India, Bangladesh, Sudan, South Sudan, Brazil, and Ethiopia6. In Brazil, VL is expanding territorially, with autochthonous cases in 25% of Brazilian municipalities in 21 of 26 states7. VL is characterized by epidemiological changes, especially the urbanization of the disease, with human and canine cases reported in medium and large-sized cities and also expansion to other Brazilian macroregions8.
Prior studies have presented data regarding the epidemiology of VL, however the analyses were performed based on municipality9,10 or state11,12 and based on hospital data13,14. One study examined the mortality of VL in Brazil3 with another focused on the prediction of VL in Brazil15. This study analyzed the incidence rate of VL throughout Brazil, using the notifications sent to the Information System for Notifiable Diseases [Sistema de Informação de Agravos de Notificação (SINAN)] from 2001 to 2014. Considering that VL is an obligatory notifiable disease16, these data represent all of the cases diagnosed in the country. Therefore, this large, comprehensive data set enabled a description of changes in the epidemiology of VL in Brazil, between the periods from 2001 to 2006 and from 2007 to 2014.
METHODS
Study population
A descriptive study was undertaken on all new human cases of VL reported to the SINAN of the Brazilian Ministry of Health from 2001 to 2014, analyzing the municipality of patient residence. The database was subgrouped into two periods, 2001-2006 and 2007-2014, due to the changes in SINAN in the year 2007. No patient identification information of the subjects was provided during this study, so Human ethics approval was not required. Notifications in relation to relapses and transferred or unconfirmed VL cases were excluded. Estimates regarding population data were obtained from the Federal Audit Court (TCU), which were used to determine quotas of the Municipal Participation Fund (MPF); data was available from the Department of Information of the Unified Health System [Departamento de Informática do Sistema Único de Saúde (DATASUS)], according to the Brazilian Institute of Geography and Statistics [Instituto Brasileiro de Geografia e Estatística (IBGE)]17.
Analysis of regional incidence rates and spatial distribution
The variables examined in this analysis were derived from the notification records18,19. The analysis was stratified by sex, based on age group, years of education, type of living area, macroregion, evolution of cases and clinical manifestations.
The incidence rates of VL were calculated for each state and region of Brazil by using the number of cases as the numerator and the number of people in the total population as the denominator and multiplying the outcome by 100,000 inhabitants. For each period studied (2001-2006 and 2007-2014), the population data was aggregated in the denominator. In this way, the numerator contains the sum of cases, and the denominator contains the sum of the population.
The crude incidence rates were calculated for the periods 2001-2006 and 2007-2014. Sociodemographic data, clinical data, and evolution of the cases were analyzed across the two time periods and presented by number and percentage.
The variation of the crude rate of VL per 100,000 inhabitants between the two time periods was calculated as follows: ((x1-x2)/x2 * 100), where x1 = crude rate in VL in the period from 2007 to 2014 and x2 = crude rate in VL in the period from 2001 to 2006. This rate variation was compared to the average VL disease incidence rate variation in Brazil. The value was multiplied by 100 and therefore represented the percentage variation of the rate. The decision to group and compare 6 and 8 years, 2001 to 2006 and 2007 to 2014, respectively, does not affect the result interpretation, as a previous test found interpolation in the confidence interval of the mean incidence rates between the groups.
The spatial distribution of VL was evaluated by mapping incidence rates per municipality and comparing distributions between the periods 2001-2006 and 2007-2014. For the mapping, QGIS® software version 2.16 was used. The categorization of rates by mapping, as shown in Table 1, followed the cutoff suggested by the Pan American Health Organization20.
Years | ||||
---|---|---|---|---|
States and macroregions | 2001-2006 | 2007-2014 | ||
crude rate | cases | crude rate | Cases | |
Rondônia* | 0.02 | 2 | 0.1 | 12 |
Amazonas* | 0.03 | 5 | 0.04 | 11 |
Roraima | 2.5 | 56 | 2.4 | 85 |
Pará | 4.7 | 1,888 | 4.1 | 2,495 |
Amapá* | 0.1 | 3 | 0.1 | 3 |
Tocantins | 15.2 | 1,143 | 26.2 | 2,885 |
North | 3.7 | 3,097 | 4.3 | 5,491 |
Maranhão | 9.7 | 3,454 | 7.3 | 3,844 |
Piauí | 8.6 | 1,516 | 6.8 | 1,705 |
Ceará | 4.4 | 2,064 | 6.2 | 4,257 |
Rio Grande do Norte | 2.7 | 480 | 2.8 | 713 |
Paraíba | 1.2 | 249 | 0.9 | 284 |
Pernambuco | 1.5 | 736 | 1 | 681 |
Alagoas | 3.2 | 569 | 1 | 262 |
Sergipe | 2 | 226 | 2.7 | 454 |
Bahia | 2.8 | 2,269 | 2.3 | 2,683 |
Northeast | 3.9 | 11,563 | 3.5 | 14,883 |
Minas Gerais | 2.1 | 2,360 | 2.3 | 3,612 |
Espírito Santo | 0.11 | 21 | 0.1 | 28 |
Rio de Janeiro | 0.02 | 21 | 0.03 | 33 |
São Paulo | 0.4 | 985 | 0.5 | 1,699 |
Southeast | 0.7 | 3,387 | 0.8 | 5,372 |
Paraná* | 0.024 | 15 | 0.019 | 17 |
Santa Catarina* | 0.006 | 2 | 0.014 | 7 |
Rio Grande do Sul | 0.003 | 2 | 0.02 | 19 |
South | 0.01 | 19 | 0.02 | 43 |
Mato Grosso do Sul | 8.4 | 1,108 | 9.2 | 1,796 |
Mato Grosso | 0.6 | 102 | 1.5 | 365 |
Goiás | 0.5 | 157 | 0.6 | 297 |
Distrito Federal | 0.2 | 33 | 0.5 | 99 |
Central-West | 1.9 | 1,400 | 2.2 | 2,557 |
Brazil | 1.808 | 19,466 | 1.833 | 28,346 |
Note: The crude rates, which are depicted in bold and italic, indicate high rates (23.05- 46.44), while the lines in bold are considered averages (9.56-23.04). *No reported autochthonous human cases in the study period. The difference between reports and the total number of cases (47 cases) can be explained by missing data of the municipality.
Data analysis was performed using Stata 9.0 software (StataCorp). Data preparation was conducted in a spreadsheet using Microsoft Excel 2010 (Microsoft Corp, Redmond, WA, USA).
RESULTS
From 2001 to 2014, 47,859 new cases of VL in Brazil were reported to SINAN, with an annual average of 3418.6 cases [standard deviation (SD) = 397.0] and an average incidence rate of 1.8 per 100,000 inhabitants. The lowest rate was recorded in 2002 (1.5/100,000 inhabitants) and the highest in 2011 (2.0/100,000 inhabitants). There was an increase in the number of municipalities reporting VL cases (Figure 1), especially in the interior area of Brazil, which was previously an unregistered area. The percentage of municipalities reporting VL cases ranged from 11.7% in 2002 to 16.8% in 2014 (data not shown). Some states experienced a change in VL epidemiology and started to register vector or canine and human cases (Figure 1).
The number of VL cases in the Southern macroregion of Brazil increased; the region changed from a disease-free area in the first period (2001-2006) to an area with autochthonous canine and human VL cases in the second period (2007-2014). The Rondônia and Amapá States recorded autochthonous canine and vector cases, respectively. Only the Amazonas and Acre states did not have records of vector, canine, and human autochthonous cases. Cases reported in Amazonas were considered not to be of autochthonous origin. These data were added in the manuscript to reinforce the need for VL surveillance, but were not the object of primary data collection. However, a literature review was used to gather this information and to make a relationship with the findings in this study.
Aside from the Northeast region, all regions in Brazil showed an increase in the crude rate of VL in the second study period (2007-2014). The Northern macroregion had the highest rate of VL, due to cases in the state of Tocantins, which had the highest number of new cases in the country (Table 1). The highest rates in each macroregion were in the Mato Grosso do Sul, Minas Gerais, Tocantins, and Maranhão States, which are located in the Central-West, Southeast, North and Northeast macroregions, respectively.
The mean percentage rate variation of VL in Brazil between the first and second period had increased by 1.38%. Compared to the average percentage rate variation in Brazil, there were a number of states that experienced an increase of more than 100%, and a few states that experienced a decrease (Figure 2).
The profile of VL cases is outlined in Table 2. This data showed a predominance in males, the under 9 year age group, less educated members of the population or people who were not yet in school, and those living in urban areas. While the rate of patients treated for VL decreased in the second period, there was an increase in the number of deaths, especially in males. A predominance of treatment abandonment was observed in males.
Characteristics | 2001-2006 (n = 19,496) | 2007-2014 (n = 28,363) | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
male | female | total | male | female | Total | |||||||
n | % | n | % | n | % | n | % | n | % | n | % | |
Years group 0 to 4 | 3,517 | 50.3 | 3,477 | 49.7 | 6,994 | 35.9 | 4,697 | 50.9 | 4,522 | 49.1 | 9,219 | 32.5 |
5 to 9 | 1,998 | 54.7 | 1,654 | 45.3 | 3,652 | 18.7 | 1,437 | 54.3 | 1,209 | 45.7 | 2,646 | 9.3 |
10 to 19 | 1,679 | 59.6 | 1,138 | 40.4 | 2,817 | 14.4 | 1,596 | 63 | 938 | 37 | 2,534 | 8.9 |
20 to 39 | 2,741 | 75.3 | 901 | 24.7 | 3,642 | 18.7 | 3,959 | 76.4 | 1,221 | 23.6 | 5,180 | 18.3 |
≥ 40 | 1,766 | 73.9 | 625 | 26.1 | 2,391 | 12.3 | 6,120 | 69.7 | 2,664 | 30.3 | 8,784 | 31.0 |
Years of education | ||||||||||||
none | 665 | 71 | 272 | 29 | 937 | 4.8 | 504 | 73 | 186 | 27 | 690 | 2.4 |
1 to 8 | 3,024 | 67.7 | 1,441 | 32.3 | 4,465 | 22.9 | 5,554 | 71.8 | 2,185 | 28.2 | 7,739 | 27.3 |
≥ 9 | 757 | 68.8 | 344 | 31.2 | 1,101 | 5.7 | 1,401 | 69.7 | 610 | 30.3 | 2,011 | 7.1 |
not apply or ignored | 7,255 | 55.8 | 5,738 | 44.2 | 12,993 | 66.6 | 10,350 | 57.7 | 7,573 | 42.3 | 17,923 | 63.2 |
Type living area | ||||||||||||
urban | 7,745 | 60.8 | 4,988 | 39.2 | 12,733 | 65.3 | 12,851 | 62.9 | 7,586 | 37.1 | 20,437 | 72.1 |
non-urban | 3,956 | 58.5 | 2,807 | 41.5 | 6,763 | 34.7 | 4,958 | 62.6 | 2,968 | 37.4 | 7,926 | 27.9 |
Macroregion | ||||||||||||
Northeast | 6,845 | 59.8 | 4,605 | 40.2 | 11,450 | 58.7 | 9,395 | 63.6 | 5,388 | 36.4 | 14,783 | 52.1 |
North | 1,827 | 58.5 | 1,297 | 41.5 | 3,124 | 16 | 3,246 | 59.3 | 2,232 | 40.7 | 5,478 | 19.3 |
Southeast | 2,036 | 61.4 | 1,280 | 38.6 | 3,316 | 17 | 3,353 | 62.9 | 1,979 | 37.1 | 5,332 | 18.8 |
Central-West | 983 | 61.8 | 608 | 38.2 | 1,591 | 8.2 | 1,779 | 65.4 | 943 | 34.6 | 2,722 | 9.6 |
South | 10 | 66.7 | 5 | 33.3 | 15 | 0.1 | 36 | 75 | 12 | 25 | 48 | 0.2 |
Evolution of cases | ||||||||||||
cure | 9,237 | 59.8 | 6,202 | 40.2 | 15,439 | 79.2 | 12,564 | 62.3 | 7,602 | 37.7 | 20,166 | 71.1 |
death | 870 | 61.8 | 539 | 38.3 | 1,409 | 7.2 | 1,541 | 66.9 | 763 | 33.1 | 2,304 | 8.1 |
ignored | 586 | 58.1 | 423 | 41.9 | 1,009 | 5.2 | – | - | - | - | - | - |
abandonment | – | - | - | - | - | - | 124 | 70.5 | 52 | 29.5 | 176 | 0.6 |
death by VL | – | – | – | – | – | – | 1,168 | 65.7 | 611 | 34.3 | 1,779 | 77.2 |
death by other causes | – | – | – | – | – | – | 373 | 71 | 152 | 29 | 525 | 22.8 |
Clinical manifestations | ||||||||||||
fever | 10,800 | 60.1 | 7,159 | 39.9 | 17,959 | 92.1 | 16,310 | 62.7 | 9,695 | 37.3 | 26,005 | 91.7 |
weakness | 8,242 | 60.7 | 5,333 | 39.3 | 13,575 | 69.6 | 14,023 | 63.5 | 8,053 | 36.5 | 22,076 | 77.8 |
weight loss | 7,964 | 61.6 | 4,957 | 38.4 | 12,921 | 66.3 | 12,831 | 65.2 | 6,860 | 34.8 | 19,691 | 69.4 |
cough or diarrhea | 4,956 | 60.2 | 3,271 | 39.8 | 8,227 | 42.2 | 7,803 | 63.2 | 4,541 | 36.8 | 12,344 | 43.5 |
splenomegaly | 9,751 | 59.9 | 6,531 | 40.1 | 16,282 | 83.5 | 13,538 | 62.5 | 8,117 | 37.5 | 21,655 | 76.3 |
hepatomegaly | 8,563 | 60 | 5,713 | 40 | 14,277 | 73.2 | 11,721 | 62.4 | 7,070 | 37.6 | 18,791 | 66.3 |
HIV coinfection | 261 | 71.7 | 103 | 28.3 | 364 | 1.9 | 1,178 | 76.7 | 357 | 23.3 | 1,535 | 5.4 |
VL: visceral leishmaniasis; HIV: human immunodeficiency virus; -: absence of the variable on the report form. Note: The percentages of the variables by sex were calculated per row and the total percentages were calculated per column. The percentages of the variable death by VL and death by other causes were calculated from the total number of deaths, which included death by all causes, because the databases were not the same in the two periods. The number of evolution cases did not reach the total number of cases in both periods due to missing data. From 2007 to 2014, the percentage of evolution cases was not 100%, because there were 7% of cases with the outcome transfer.
The incidence of visceral leishmaniasis coinfection with human immunodeficiency virus (VL-HIV) was highest in the second period (2007-2014), although there was missing data for this variable, as well as information on other clinical manifestations (Table 2).
DISCUSSION
This study identified geographical changes in the incidence of VL cases, which was characterized by an expansion to previously disease-free areas; mainly to the interior of Brazil (Figure 1), but there is also a strong urbanization component (Table 2). Despite the stability of the incidence rate in Brazil in general (Figure 2), the increase in the number of municipalities that commenced reporting cases of this important neglected disease deserves attention, as it confirms the territorial expansion.
The expansion over time (Figure 1), as shown by VL-free areas in the first period (2001-2006) that reported vector, canine, or human cases in the second period (2007-2014), can be explained by simple adaptation of sandflies to varying temperatures and to the peridomiciles21, migratory movement of people with VL-contaminated dogs22, and to locals at the borders who have reported the disease23. Some states that were previously considered disease-free started to report the first autochthonous canine cases; for example Rondônia24,25 in 2010 and Paraná26 and Santa Catarina27,28 in 2012. These areas require active epidemiologic surveillance, because canine cases precede the human cases21.
The State of Amapá has had no reported human autochthonous cases to date, although it is located on the border with Pará, which has records of human VL cases dating back to 193429 and is part of the Guianan Ecoregion Complex (GEC), with autochthonous human cases in Venezuela and Northern Brazil30. However, this state reported the presence of the sandfly Lutzomyia longipalpis in 2013 for the first time31, which increases the potential risk of the disease in humans at this site, as also pointed out in State of Rondonia32. In the State of Rio Grande do Sul, vector and canine cases were identified in 2008, and the first autochthonous human case was identified in 2009 in the municipality of São Borja33, which probably originated in the provinces of Argentina that border this municipality23,34.
The changes in VL epidemiology are also reflected in the extensive urbanization of cities, as the disease is predominantly associated with urban areas (Table 2). Previously, VL was considered a rural disease, however in the early 1980s VL epidemics were recorded in an urban environment in Brazil8. The switch to urban areas is corroborated by the current study, in which approximately 70% of the cases were residents of urban areas. Although it is not possible to state that urban transmission is different from that in rural areas, some factors that might be involved in the process of urbanization of VL in Brazil are the environmental modifications caused by anthropic action, caused by migratory movements and nonplanned urban occupation together with poor sanitation35. In addition, the main vector of VL, the sandfly Lutzomyia longipalpis, has adapted to the peridomicile, especially in the presence of domestic animals such as dogs18.
The stability of the crude incidence rate in Brazil between the two time periods indicates that, even with an increase in the number of reported cases, there is no increase in the incidence rate, when population growth is considered. Therefore, the incidence of VL had increased, especially when analyzed per macroregion and state separately. In the 1990s, approximately 90% of the reported cases of VL were located in the Northeastern macroregion. With the spread of the disease to other regions, this situation has changed and a decrease in the number of VL cases has been observed in the Northeast21, whereas the North reported an increase in cases (Table 1). This expansion of the incidence rate is unrelated to the decision to subgroup the database into two time periods with an unequal number of years, considering that if an equal number of years were compared in the groups, a similar mean and interpolation in the confidence interval was found (data not shown). With respect to the spatial distribution, it is evident that municipalities that were previously free from the disease had changed their status.
Although the State of Maranhão reported the highest rate in the Northeast region (Table 1) and has a history of high incidence, overall there was a reduction in incidence reported. This may be due to the surveillance efforts in the Northeast, or to the emigration of the population36. For example, the State of Piauí previously had the highest incidence in the Northeast37 and now has a reported reduction, even though it is still the state with the second highest incidence in the Northeast. This expansion of the incidence in VL cases to other regions has occurred without a sufficient amount of time to organize health services for diagnoses, perform appropriate follow-up of the cases, and to train health professionals, who without the correct knowledge might recommend an inappropriate treatment regimen for patients with VL2,38.
The North of Tocantins State (Figure 1) had a remarkably high VL incidence rate when the geographical distribution was analyzed. This is mainly due to the extensive migratory flow of people from the Northeast, who live in poor housing conditions that lack urban infrastructure, sanitation and essential public services, such as garbage collection, health care and education39. These demographic and social problems are associated with the environmental impact caused by the deforestation in this state40. This situation contributed to the epidemiologic situation identified in Tocantins, which deserves greater efforts to be made by the health services, both in the diagnosis and treatment of cases, as well as towards adequate surveillance services for this particular population.
Of all of the characteristics of VL, the predominance of the disease in males requires attention. The frequency of this disease in men increased with age (Table 2). Physiological factors are the most likely cause for the increased risk in males, indicating that from a certain age, sex hormones and the immune system in men result in a higher susceptibility to infection and disease41.
An increase of the disease in adults older than 40 years is noted in the second period of this study (Table 2), which can also be attributed to HIV coinfection42–44. From 2008, the number of adult patients exceeded the number of children with VL in Ceará, however, VL-HIV coinfection was predominant in the 20-39 year age group12. This phenomenon has also occurred in Southern Europe, North Africa, and Western and Central Asia. Since the beginning of HIV infections and increased use of immunosuppressants for transplantation and chemotherapy, approximately half of the VL cases in Europe are adults45. Therefore, HIV infection should also be considered in our study. Although a poorly performed routine HIV test was done, the data show an increase in the incidence of VL-HIV coinfection (Table 2).
The percentage of patients that were cured (Table 2) is lower than the number suggested by the Pan American Health Organization, which advocates that at least 95% of patients treated for VL are cured46. A study that was performed in the City of Bauru showed that 90.3% of the treated VL patients were cured9. The low number of patients that were cured may reflect the performance of the VL control program and the records in the information system. It is not known if the patients were cured or if there was a problem with the information system, since the evolution of the cases showed an increase of missing data; from 7% in the period 2001-2006 to 13.2% in the period 2007-2014 (data not shown). This is a very serious problem as the information on the percentage of patients that were cured is critically important and related to the capacity of the health services to perform early diagnosis and to have disposal resources such as materials, laboratory, medicines and trained professionals to give the correct treatment to cases. This low percentage of patients that were cured indicates that it is necessary to evaluate the data record for completeness and quality47, as well as to have correct follow-up of the cases, in order to avoid abandoning of the treatment.
Failure to complete the evolution of the case, especially the cure information, is a concerning issue, because the absence of a cure contributes to unfavorable outcomes, such as abandonment and death. In addition to the increase in deaths in the second period, this may be an underestimate due to missing data. Although some patients survive even when they are not cured, they may have subclinical disease and with a return to illness in case of decreased immunity48, and may worsen if there is coinfection, especially with HIV49,50 or malaria, which has already been described in Brazil51 and Africa52.
VL should be suspected when a patient presents with fever and splenomegaly that might be associated with hepatomegaly18. These symptoms, which characterize the initial phase of the disease, were the most frequent in this study, together with weakness and weight loss. However, the latter are also observed in other infectious diseases18, which may cause confusion and delay the diagnosis of VL, thereby compromising the condition of the patient with malnutrition, bleeding and other bacterial infections that can lead to death. Therefore, there is a need for trained health professionals38 and sensitizing health teams to recognize this important neglected disease in addition to the installation of adequate infrastructure for prompt laboratory diagnosis. In addition, health service structuring is needed in order to optimize epidemiologic surveillance, as well as vector control measures and inclusion of new methods of disease control in dogs.