ABSTRACT:

Objective

To analyze the spatiotemporal pattern of tuberculosis incidence and its associated factors.

Methods:

Ecological study, which used tuberculosis notifications during 2001-2017 and as units of analysis the municipalities of the state of Ceará of the Notifiable Diseases Information System. Time pattern analysis techniques and geographically weighted regression were used.

Results:

The gross incidence rate in the state reached a peak of 226.1/100 thousand inhabitants with a significant decrease of 1.9% per year (95%CI −3.0–−0.7). The main clusters were identified in Fortaleza and its metropolitan region, as well as in the Sobral region. The socioeconomic indicators that were associated with the incidence were: occupation in the agricultural sector and services, population in households with bathrooms and running water and vulnerable to poverty.

Conclusion:

There was a significant decrease in the incidence of the disease in the state of Ceará. The spatial clusters were mostly located in areas with high population and the indicators most related to the formation of spatial clusters were related to employment/income, housing and vulnerability.

Keywords:
Tuberculosis; Incidence; Spatial analysis; Epidemiology

INTRODUCTION

Tuberculosis (TB) is an infectious disease of importance to global public health. A report by the World Health Organization (WHO) points out that about 10 million new cases of TB have been reported on a global scale. Of these, more than one million people died, making the disease one of the ten leading causes of death worldwide¹1 World Health Organization. Global Tuberculosis Report 2019. Geneva: World Health Organization; 2020. [cited on Aug 1, 2020]. Available at: https://apps.who.int/iris/bitstream/handle/10665/329368/9789241565714-eng.pdf
https://apps.who.int/iris/bitstream/hand... . In addition, Brazil was the only country on the American continent listed among the 30 with a high burden of TB and occupied the first position on the continent, with 32% of cases²2 Aridja UM, Gallo LG, Oliveira AFM, Silva AWM, Duarte EC. Casos de tuberculose com notificação após o óbito no Brasil, 2014: um estudo descritivo com base nos dados de vigilância. Epidemiol Serv Saude 2020; 29(5): e2020060. https://doi.org/10.1590/S1679-49742020000500014
https://doi.org/10.1590/S1679-4974202000... ,³3 Brasil. Ministério da Saúde. Secretária de Vigilância em Saúde. Departamento de vigilância das Doenças Transmissíveis. Manual de recomendações para o controle da tuberculose no Brasil. Brasília: Ministério da Saúde; 2019. [cited on Aug 10, 2020]. Available at: https://sbpt.org.br/portal/wp-content/uploads/2019/06/manual_recomendacoes_tb_2ed_atualizada_8maio19.pdf
https://sbpt.org.br/portal/wp-content/up... .

Although Brazil is among the countries with the most cases in the world, from 2011 to 2016 there was a decrease in this number, with a lower coefficient presented in 2016, with 34.4 cases per 100,000 inhabitants. On the other hand, in the period 2017-2019, there is an increase in incidence in the country⁴4 Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Tuberculose 2021. Boletim Epidemiológico; 2021 (número especial). [cited on Feb 09, 2022]. Available at: https://www.gov.br/saude/pt-br/media/pdf/2021/marco/24/boletim-tuberculose-2021_24.03
https://www.gov.br/saude/pt-br/media/pdf... . In Ceará, the decrease occurred in the period from 2003 to 2017, with an average reduction of 1.9% per year. As a result, the state had the lowest incidence ever recorded: 34.1 cases per 100,000 inhabitants in 2017⁵5 Souza J, Oliveira KS, Ávila TT, Quadros SR, Zilly A, Sobrinho RAS. Incidência da tuberculose e a correlação entre a realização do tratamento e a cura. Rev Enferm UFSM 2018; 8(4): 637-48. https://doi.org/10.5902/2179769230534
https://doi.org/10.5902/2179769230534... .

TB is associated with some population groups considered more susceptible to illnesses⁶6 Silva PHS, Cirilo SSV, Sousa Junior SC, Cruz VT, Correira RS, Santos AF, et al. Aspectos sociodemográficos e clínicoepidemiológicos da tuberculose em um município do nordeste brasileiro. Revista Eletrônica Acervo Saúde 2020; (12):5: e1916. https://doi.org/10.5123/s1679-49742017000200014
https://doi.org/10.5123/s1679-4974201700... . Furthermore, studies indicate that TB is also related to socioeconomic conditions observed in so-called developing countries, such as situations of poverty, malnutrition, poor sanitary conditions, and high population density⁷7 Pereira AGL, Escosteguy CC, Valencia LIO, Magalhães MAFM, Medronho RA. Análise espacial de casos de tuberculose e associação com fatores socioeconômicos: uma experiência no município do Rio de Janeiro. Cad Saúde Colet 2018; (26)2: 203-10. https://doi.org/10.1590/1414-462x201800020013
https://doi.org/10.1590/1414-462x2018000... .

In view of the above, the importance of using epidemiological methods of spatial analysis in health⁸8 Martins-Melo FR, Lima MS, Ramos Jr 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): e93770. https://doi.org/10.1371/journal.pone.0093770
https://doi.org/10.1371/journal.pone.009... for accurate knowledge of places that present clusters of high rates of TB is highlighted. Thus, this study aimed to analyze the spatiotemporal pattern of TB incidence and its associated factors.

METHODS

This is an ecological study, which used the municipalities of the state of Ceará as units of analysis. It has 184 municipalities and the eighth largest population in the country, estimated at 9.13 million inhabitants. Of the 27 states of the federation, Ceará has the 11^th highest population density, with 56.76 inhabitants per km²2 Aridja UM, Gallo LG, Oliveira AFM, Silva AWM, Duarte EC. Casos de tuberculose com notificação após o óbito no Brasil, 2014: um estudo descritivo com base nos dados de vigilância. Epidemiol Serv Saude 2020; 29(5): e2020060. https://doi.org/10.1590/S1679-49742020000500014
https://doi.org/10.1590/S1679-4974202000... and a nominal monthly per capita income of R$ 850⁹9 Instituto Brasileiro de Geografia e Estatística. Pesquisa de orçamentos familiares 2017-2018: primeiros resultados 2019. [cited on Oct 02, 2020]. Available at: https://www.ibge.gov.br/estatisticas/sociais/populacao/24786-pesquisa-de-orcamentos-familiares-2.html?edicao=25578&t=publicacoes
https://www.ibge.gov.br/estatisticas/soc... .

The data in this study are of the secondary type and were obtained from the Secretary of Health of the state of Ceará (Secretaria da Saúde do Estado do Ceará – SESA), after prior authorization to access the information. The records contain information on TB from the Notifiable Diseases Information System (Sistema de Informação de Agravos de Notificação – SINAN).

This investigation comprises all cases reported as “new TB cases” during the period from January 1^st, 2001 to December 31^st, 2017 in the state of Ceará. Thus, it had a population of 60,614 new TB cases during the period. The mean incidence of TB in the period studied was calculated for each municipality using the average number of cases in the numerator and the population of the central year (2009) as the denominator, multiplied by 100 thousand inhabitants. Population data from the municipalities of Ceará were obtained from the website of the Brazilian Institute of Geography and Statistics (Instituto Brasileiro de Geografia e Estatística – IBGE).

Initially, the characteristics of the cases were described, including the number of missing data in each variable. After that, the analysis of the temporal trend of the incidence of TB was conducted through regression by inflection points. This analysis indicates whether one or more line segments should be added in a linear regression to indicate any change in time trend, which refutes the null hypothesis that no points should be added. Its results are possible to estimate the annual percentage change (APC) of the studied trend, as well as its 95% confidence interval (95%CI) and its statistical significance¹⁰10 Sousa GJB, Garces TS, Pereira MLD, Moreira TMM, Silveira GMD. Temporal pattern of tuberculosis cure, mortality, and treatment abandonment in Brazilian capitals. Rev Latino-Am Enfermagem 2019; 27: e3218. https://doi.org/10.1590/1518-8345.3019.3218
https://doi.org/10.1590/1518-8345.3019.3... .

In their results, significantly positive APC indicate a growth trend of the phenomenon studied, and negative APC indicate a decrease in the trend. Non-significant APC values indicate a stationary trend. Thus, the model was adjusted assuming that the number of inflection points could vary from zero to two, which means a straight line or a line with two or three segments. The 5% significance level was established in all trend analyses of all times¹⁰10 Sousa GJB, Garces TS, Pereira MLD, Moreira TMM, Silveira GMD. Temporal pattern of tuberculosis cure, mortality, and treatment abandonment in Brazilian capitals. Rev Latino-Am Enfermagem 2019; 27: e3218. https://doi.org/10.1590/1518-8345.3019.3218
https://doi.org/10.1590/1518-8345.3019.3... .

For the spatial analysis, the thematic map of the average gross incidence of TB in the municipalities of Ceará was initially created, and then these rates were smoothed using the local empirical Bayesian method to reduce the instabilities caused by the gross rates. The application of this method is necessary since it generates rates closer to reality, as it considers not only the rate value of a given municipality, but weights it in relation to those that border it through a spatial proximity matrix. For the construction of the matrix, the criterion of contiguity was taken into account, assigning the value 1 to municipalities that have common borders and 0 to municipalities that do not share borders.

For the identification of spatial clusters of TB, two different methods were used, though with the same purpose. The first concerns the spatial autocorrelation function, through the global and local Moran index. The global index was used to test the spatial dependence hypothesis and provide an overall measure of association for the entire study area. Once the presence of global spatial autocorrelation was verified, the local Moran index (Local Index Spatial Analysis — LISA) was applied to verify the presence of spatial aggregates and quantify the degree of spatial association in each municipality of the sample set, considering it if p<0.05.

The results of the local Moran index are presented using the Moran and the Lisa maps. The first one allows you to graphically visualize the degree of similarity between neighbors and is represented by four quadrants: in red are the municipalities with high rates and that are close to municipalities with equally high rates (high/high spatial pattern); in green are the municipalities that have low rates and that are surrounded by municipalities that also have low rates (low/low spatial pattern); the municipalities in yellow (high/low spatial pattern) and blue (low/high spatial pattern) represent areas of epidemiological transition and that exhibit high and low rates, but are very close to municipalities that have low and high rates, respectively⁸8 Martins-Melo FR, Lima MS, Ramos Jr 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): e93770. https://doi.org/10.1371/journal.pone.0093770
https://doi.org/10.1371/journal.pone.009... .

Aiming at detecting spatial clusters, as well as areas that present greater risk for TB, a purely spatial scan analysis, called Scan statistic, was carried out. For the identification of clusters, the discrete Poisson model was adopted, which considers that the number of events in each area is proportional to the size of the population⁸8 Martins-Melo FR, Lima MS, Ramos Jr 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): e93770. https://doi.org/10.1371/journal.pone.0093770
https://doi.org/10.1371/journal.pone.009... , as well as the following requirements: no geographical overlap of clusters, maximum cluster size equal to 50.0% of the exposed population, clusters in a circular shape and 999 replications. The Scan statistic also made it possible to calculate the relative risk (RR) in each municipality in Ceará; those with values >1 have a relative risk for TB that is higher than the risk for Ceará as a whole.

To identify the influence of socioeconomic indicators on TB incidence, the ordinary least squares (OLS) regression method was used. The indicators used were: Theil L index, Gini index, unemployment rate, employed with complete high school education, employed in the agricultural sector, employed in the mining sector, employed in the manufacturing industry, employed in the construction sector, employed in the commerce sector, employed in the service sector, population in households with toilet and running water, population in households with density >2, municipal human development index (MHDI), illiteracy rate, and vulnerable to poverty. All these indicators were taken from the Atlas of Human Development in Brazil¹¹11 Atlas do Desenvolvimento Humano no Brasil. Acervo. Entenda o que é o projeto Atlas, o conceito de desenvolvimento humano, o IDHM e sua metodologia [cited on Feb 10, 2022]. Available at: http://www.atlasbrasil.org.br/acervo/biblioteca
http://www.atlasbrasil.org.br/acervo/bib... .

Initially, these variables were entered into a multiple linear regression (OLS) model. The variables were selected using backward selection, while the presence of multicollinearity was assessed using the variance inflation factor (VIF). Those that remained in the final OLS regression model (p<0.05) were also inserted into a geographically weighted regression (GWR) model, since it is able to analyze phenomena that vary according to the area in which they are studied. The GWR generated a regression coefficient for each of the 184 municipalities in the state, as a social indicator can act as a risk factor in a certain territory and a protective factor in another¹²12 Charlton M, Fotheringham S, Brunsdon C. Geographically weighted regression. United Kingdon: ESRC National Centre for Research Methods; 2005.. The result of the GWR regression was presented in the format of two thematic maps for each socioeconomic indicator: one for the value of the regression coefficient and another that represents the statistical significance, considering p<0.05.

The calculation of the local empirical Bayesian rate and the spatial autocorrelation test was performed using the TerraView 4.2.2^® software, while the purely spatial scan analysis was performed using the SaTScan 9.6^® software. The non-spatial regression OLS was performed in the software Stata 12^® and the spatial regression GWR, in the software GWR4.0.9^®. All maps were produced in QGIS 2.4.17^® software.

The project of which this study is a part was sent for ethical consideration by the State University of Ceará, being approved under protocol No. 2.687.046.

RESULTS

From 2001 to 2017, 60,614 new TB cases were reported in Ceará. Table 1 shows that the cases had a median age of 38 years (interquartile range — IQR 25–52) and were mainly male (62.1%; n=37,649). Most self-declared brown (74.2%; n=38,859), had less than high school education (78.4%; n=34,101), and lived in urban areas (84.1%; n=49,134).

The mean incidence of the disease in the period was 28.6/100 thousand inhabitants, with the lowest incidence occurring in 2017 (34.1 cases per 100 thousand inhabitants) and the highest in 2003 (61.2 cases per 100 thousand inhabitants). The analysis of the temporal pattern of TB incidence in the studied period showed a significant decrease of 1.9% per year (95%CI −3.0–−0.7; p=0.003). Furthermore, the regression did not insert any inflection points, which means that there was no change in trend over time, which is explained only by a line segment (data not shown in tables or graphics).

In Figure 1A, it is possible to observe that all municipalities in Ceará registered at least one case of TB during the period 2001–2017, since there are no municipalities with a rate equal to zero. The gross incidence rate in the state reached a maximum value of 226.1/100 thousand inhabitants. With the smoothing of the gross rates by the local empirical Bayesian method (Figure 1B), it is possible to perceive a more apparent spatial pattern, with the aggregation of municipalities with higher rates close to Fortaleza and Sobral. It is noteworthy that the municipality of Itaitinga, in the metropolitan region of Fortaleza, had the highest gross (226.1/100 thousand inhab.) and Bayesian (168.1/100 thousand inhab.) rates of TB.

In the spatial autocorrelation analysis, the global Moran index was equal to 0.32 (p=0.01), indicating positive spatial autocorrelation. Through this index, spatial clusters can be identified (Figure 1C), in which the high-high pattern of distribution is found in Fortaleza and its metropolitan region, as well as in the region of Sobral, an important municipality in the countryside of the state. The main low-low clusters are concentrated in the Sertão Cearense and Sul Cearense mesoregions. The Lisa map (Figure 1D) shows the intensity of statistical significance of the spatial clusters pointed out by the Moran map.

Using the Scan technique, it was possible to identify that the population of Itaitinga has 4.8 times more risk of contracting TB when compared to the population of other municipalities in Ceará (Figure 1E). The technique also made it possible to identify four statistically significant (p<0.05) clusters of TB in the state. The most likely cluster, that is, the one that has the lowest probability of having occurred by chance (in red), was composed of Fortaleza and five municipalities in its metropolitan region (Caucaia, Maranguape, Maracanaú, Eusébio, and Itaitinga), whose RR was of 2.75 (Figure 1F). The characterization of the other clusters can be found in Table 2.

Table 3 shows the significant socioeconomic indicators in the OLS and GWR regression models. In the first model, the value of the Moran index (I) was 0.077 (p=0.04) and it could explain 27.0% of the outcome variation. The significant value indicates that the regression residuals have a spatial pattern and that the incidence can be explained by the neighbors of the neighbors. When applying the GWR regression, 32.0% of the model was explained and I=0.04.

Table 3 presents the final OLS non-spatial regression and GWR spatial regression models for the factors associated with TB incidence in Ceará. The socioeconomic indicators that were associated with the incidence of TB were: employed in the agricultural sector, employed in the service sector, population in households with bathroom and piped water, and population vulnerable to poverty. These did not show multicollinearity, as the global VIF was equal to 2.2, which indicates a low correlation between predictors.

Figure 2 shows the thematic maps derived from the GWR regression results. It was found that, throughout the state, there was a significant negative relationship (p<0.01) between the percentage of people employed in the agricultural sector and the incidence of TB (Figures A and B). Likewise, in the north, northwest and Sertão Cearense mesoregions, a significant negative correlation (p<0.05) was also found between the percentage of people employed in the service sector and the incidence of TB (Figures C and D). In addition, an increase in the incidence of TB was observed as the percentage of people living in houses with bathroom and piped water increased (p<0.05), mainly in northern Ceará (Figures E and F). Finally, a positive correlation (p<0.05) was observed between the percentage of those vulnerable to poverty and the incidence of TB, especially in municipalities far from the capital and metropolitan region (Figures G and H).

DISCUSSION

The present study identified the predominance of male individuals, of mixed race, with low education, and residing in the urban area. In addition, a decrease in the incidence of TB was demonstrated in the period 2001–2017. On the other hand, the spatial analysis explained the formation of clusters in populous municipalities in the state of Ceará. Finally, the indicators that showed a relationship with the incidence of TB were: employed in the agricultural sector, employed in the service sector, population in households with bathroom and piped water, and population vulnerable to poverty.

The characterization of the population of this study shows that the disease remains in a specific segment. It is identified that, regardless of the region studied, illness from TB has social aspects in common throughout Brazil¹³13 Mendes AM, Bastos JL, Bresan D, Leite MS. Situação epidemiológica da tuberculose no Rio Grande do Sul: uma análise com base nos dados do Sinan entre 2003 e 2012 com foco nos povos indígenas. Rev Bras Epidemiol 2016; (19)3: 658-69. https://doi.org/10.1590/1980-5497201600030015
https://doi.org/10.1590/1980-54972016000... –¹⁷17 Santos BO, Brito TVR, Mesquita CR, Guimarães RJPS, Leão LA, Rocha MP. Análise espaço-temporal da incidência de tuberculose na atenção primária. Pará Res Med J 2017; 1(2): e21. https://doi.org/10.4322/prmj.2017.021
https://doi.org/10.4322/prmj.2017.021... . These findings help to perpetuate an image of TB linked to neglected and now stigmatized segments, which lack effective actions to eliminate the disease.

As for the temporal pattern, it was identified that there was a significant decrease in the incidence of TB in the period 2001-2017. This drop may be associated with the implementation of the Millennium Goals, which enabled the development of national programs to reduce the incidence of TB. However, this result is currently far from the recommended by the goals of the Sustainable Development Goals (SDG) and the End of TB strategy (End TB). These seek a 95% decline in the incidence of the disease¹⁸18 Uplekar M, Weil D, Lonnroth K, Jaramillo E, Lienhardt C, Dias HM, et al. WHO's new end TB strategy. Lancet 2015; 385(9979): 1799-801. https://doi.org/10.1016/S0140-6736(15)60570-0
https://doi.org/10.1016/S0140-6736(15)60... ,¹⁹19 Trajman A, Saraceni V, Durovni B. Sustainable development goals and tuberculosis in Brazil: challenges and potentialities. Cad Saúde Pública 2018; 34(6): e00030318. https://doi.org/10.1590/0102-311X00030318
https://doi.org/10.1590/0102-311X0003031... . Thus, it is pointed out that actions such as health education, the active search for contacts, the screening of people living with the human immunodeficiency virus (HIV/Aids) and the treatment of cases of latent TB infection are fundamental to achieve the goals set²⁰20 Santos AP, Silva DR, Mello FCQ. In the time of strategies to end tuberculosis, prevention is better than treatment. J Bras Pneumol 2020; 46(2): e20200017. https://doi.org/10.36416/1806-3756/e20200017
https://doi.org/10.36416/1806-3756/e2020... .

Regarding the analysis of the spatial pattern of TB incidence, spatial clusters located in Fortaleza and in municipalities in its metropolitan region were revealed, with the highest risk observed in Itaitinga, a municipality that houses a prison in its territory. It is known that this factor may have influenced the formation of clusters, since populations deprived of liberty are at 31 times greater risk of becoming ill with TB than the general population¹⁹19 Trajman A, Saraceni V, Durovni B. Sustainable development goals and tuberculosis in Brazil: challenges and potentialities. Cad Saúde Pública 2018; 34(6): e00030318. https://doi.org/10.1590/0102-311X00030318
https://doi.org/10.1590/0102-311X0003031... . In addition, the most likely cluster is located in the state capital and in municipalities in its metropolitan region, which shows the characteristic of urban agglomerations, constant transit of people between municipalities and areas with greater poverty.

As for the indicators associated with the incidence of the disease, this study found that it decreases when the proportion of employees in the agricultural and service sectors increases. This is due to the fact that the state is one of the most economically diversified in the Brazilian Northeast²¹21 Ceará. Instituto de Pesquisa e Estratégia Econômica do Ceará. Indicadores econômicos do Ceará 2020. Fortaleza: IPECE; 2020. [cited on Oct 25, 2020]. Available at: https://www.ipece.ce.gov.br/wp-content/uploads/sites/45/2020/12/Indicadores_Economicos2020.pdf
https://www.ipece.ce.gov.br/wp-content/u... . Thus, it is believed that these activities showed a significant relationship as they are two of the most important in the state, employing a large portion of the population of Ceará²¹21 Ceará. Instituto de Pesquisa e Estratégia Econômica do Ceará. Indicadores econômicos do Ceará 2020. Fortaleza: IPECE; 2020. [cited on Oct 25, 2020]. Available at: https://www.ipece.ce.gov.br/wp-content/uploads/sites/45/2020/12/Indicadores_Economicos2020.pdf
https://www.ipece.ce.gov.br/wp-content/u... . Thus, the inverse relationship of being employed in this sector helps to highlight the importance of work/income in the incidence of TB.

In addition, this study observed that the incidence of TB was positively associated with the percentage indicator of the population living in households with a bathroom and running water. This relationship proved to be significant in regions where there are large urban agglomerations, including the capital Fortaleza and Sobral, which is an important municipality in the interior of the state. Through this finding, it can be hypothesized that, since the disease is located in large centers, it is also possible to find households with greater access to toilets and running water. Thus, these two variables may present a correlation; however, it is important to emphasize that the greater social exclusion, the result of the strong stratification between the classes of society, is one of the main points to be verified regarding the social determinants of health, facilitating the transmission of TB in the community²²22 Munayco CV, Mújica OJ, León FX, del Granado M, Espinal MA. Social determinants and inequalities in tuberculosis incidence in Latin America and the Caribbean. Rev Panam Salud Publica 2015; 38(3): 177-85. PMID: 26757995,²³23 Pelissari DM, Diaz-Quijano FA. Household crowding as a potential mediator of socioeconomic determinants of tuberculosis incidence in Brazil. PLoS One 2017; 12(4): e0176116. https://doi.org/10.1371/journal.pone.0176116
https://doi.org/10.1371/journal.pone.017... .

Finally, the incidence of TB also showed a relationship with the proportion of people vulnerable to poverty. Thus, this population is more prone to exposure to the bacillus, especially in the periphery or in places of agglomeration and home density in municipalities¹⁴14 Moraes MFV, Corrêa RGCF, Coutinho NPS, Caldas AJM, Silva TC, Santos KCB, et al. Perfil epidemiológico de casos de tuberculose em um município prioritário no estado do Maranhão. Rev Pesq Saúde 2017; (18)3: 147-50.. In this sense, it is understood that this variable can mediate the relationship between other variables with the incidence of TB, such as the lack of basic health services, malnutrition, inadequate nutrition, drug abuse, lack of access to information, and low education²⁴24 Moreira ASR, Kritski AL, Carvalho ACC. Determinantes sociais da saúde e custos catastróficos associados ao diagnóstico e tratamento da tuberculose. J Bras Pneumol 2020; 46(5): e20200015. https://doi.org/10.36416/1806-3756/e20200015
https://doi.org/10.36416/1806-3756/e2020... . However, it is important to emphasize that these are just hypotheses to be tested in future studies, especially those that consider the individual level in order to avoid the ecological fallacy.

In this context, given that TB is a disease related to poverty, income transfer programs are effective in reducing social inequalities and the consequent incidence of diseases. In the case of Brazil, the Bolsa Família Program can be mentioned as a successful strategy in this regard.

Although the mechanism by which the cash transfer contributes to the reduction of the incidence of TB has not been investigated, the impact on the social determinants that influence the number of new cases is undeniable²⁵25 Nery JS, Rodrigues LC, Rasella D, Aquino D, Barreira D, Torrens AW, et al. Effect of Brazil's conditional cash transfer programme on tuberculosis incidence. Int J Tuberc Lung Dis 2017; 21(7): 790-6. https://doi.org/10.5588/ijtld.16.0599
https://doi.org/10.5588/ijtld.16.0599... . Therefore, even knowing that the program alone is not enough to contribute to the reduction of the incidence of TB, an approach that promotes the reduction of poverty and the expansion of social protection would accelerate progress toward the goals recommended by the End TB strategy²⁵25 Nery JS, Rodrigues LC, Rasella D, Aquino D, Barreira D, Torrens AW, et al. Effect of Brazil's conditional cash transfer programme on tuberculosis incidence. Int J Tuberc Lung Dis 2017; 21(7): 790-6. https://doi.org/10.5588/ijtld.16.0599
https://doi.org/10.5588/ijtld.16.0599... ,²⁶26 Carter DJ, Glaziou P, Lönnroth K, Siroka A, Floyd K, Weil D, et al. The impact of social protection and poverty elimination on global tuberculosis incidence: a statistical modelling analysis of sustainable development goal 1. Lancet Glob Health 2018; 6(5): e514-22. https://doi.org/10.1016/S2214-109X(18)30195-5
https://doi.org/10.1016/S2214-109X(18)30... .

This study has some limitations. First, the ecological design makes it impossible to establish relationships at an individual level. In this case, results must be considered only at the municipal level to avoid the ecological fallacy (inference about individual phenomena based on aggregated observations). In addition, the use of secondary databases presents incompleteness and variable information quality, which can make it difficult to interpret the results. It is known that the systems used in data collection do not have a primary character for research, but were adapted to this reality. In this sense, variables that could be considered important in this study were not included because this type of research can only use those that are maintained by the system. Even so, the results of this research are important for the epidemiological characterization of the incidence of TB and associated factors. Finally, it is suggested to carry out studies like the present one, with other designs or methods of analysis, to estimate the same effect at an individual level and strengthen the causal hypothesis.

The results of this investigation show the spatiotemporal dynamics of TB incidence in recent years. The variables related to those employed in the agricultural sector, employed in the service sector, population in households with bathroom and piped water, and population vulnerable to poverty were associated with the incidence of the disease. These findings denote the economic and social character related to the disease. Thus, focusing on these areas can help break the chain of TB transmission.

REFERENCES

Publication Dates

History