Temporal trend of tuberculosis incidence and its spatial distribution in Macapá – Amapá

Clóvis Luciano Giacomet Marcio Souza Santos Thaís Zamboni Berra Yan Mathias Alves Luana Seles Alves Fernanda Bruzadelli Paulino da Costa Antonio Carlos Vieira Ramos Juliane de Almeida Crispim Aline Aparecida Monroe Ione Carvalho Pinto Regina Célia Fiorati Marcos Augusto Moraes Arcoverde Dulce Gomes Giselle Lima de Freitas Mellina Yamamura Ricardo Alexandre Arcêncio About the authors

ABSTRACT

OBJECTIVE:

To evaluate the temporal trend of tuberculosis incidence after the implementation of the rapid molecular test (RMT-TB), to identify whether tuberculosis presents seasonal variation and to classify the territory according to case density and risk areas in Macapá, Amapá.

METHODS:

Ecological study of tuberculosis cases registered in the Sistema de Informação de Agravos de Notificação (SINAN – Information System for Notifiable Diseases) between 2001 and 2017. We used the Prais-Winsten test to classify the temporal trend of incidence and the interrupted time series to identify changes in the temporal trend before and after the implementation of the rapid molecular test, and to verify seasonality in the municipality. The Kernel estimator was used to classify case density and scan statistics to identify areas of tuberculosis risk.

RESULTS:

A total of 1,730 cases were identified, with a decreasing temporal trend of tuberculosis incidence (−0.27% per month, 95%CI −0.13 to −0.41). The time series showed no change in level after the implementation of the GeneXpert®MTB/RIF molecular test; however, the incidence increased in the post-test period (+2.09% per month, 95%CI 0.92 to 3.27). Regarding the seasonal variation, it showed growth (+13.7%/month, 95%CI 4.71 to 23.87) from December to June, the rainy season – called amazon winter season –, and decrease (−9.21% per month, CI95% −1.37 to −16.63) in the other periods. We classified areas with high density of cases in the Central and Northern districts using Kernel and identified three protection clusters, SC1 (RR = 0.07), SC2 (RR = 0.23) and SC3 (RR = 0.36), and a high-risk cluster, SC4 (RR = 1.47), with the scan statistics.

CONCLUSION:

The temporal trend of tuberculosis incidence was decreasing in the time series; however, detection increased after the introduction of RMT-TB, and tuberculosis showed seasonal behavior. The case distribution was heterogeneous, with a tendency to concentrate in vulnerable and risk territories, evidencing a pattern of disease inequality in the territory.

DESCRIPTORS:
Tuberculosis, epidemiology; Space-Time Clustering; Spatial Analysis; Ecological Studies

INTRODUCTION

According to the World Health Organization (WHO), tuberculosis is one of the top 10 causes of death in the world and the leading cause of death by a single infectious agent in people living with the human immunodeficiency virus (HIV), thus, remaining a worldwide public health problem. In 2019, about 10 million people became ill with tuberculosis worldwide and about 1.4 million people died as a result of the disease11 World Health Organization. Global Tuberculosis Report 2020. Geneva (CH): WHO; 2020 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/336069/9789240013131-eng.pdf
https://apps.who.int/iris/bitstream/hand...
.

In Brazil, in 2019, 73,864 new cases of tuberculosis were diagnosed (incidence of 35.0 cases/100,000 inhabitants) and note that the Northern region of the country is the most affected by the disease22 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde. Bol Epidemiol Tuberculose 2020. 2020 [cited 2021 Mar 15].];N° Espec:1-40. Available from: http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
http://www.aids.gov.br/pt-br/pub/2020/bo...
. In the state of Amapá, in 2019, 296 new cases of the disease were diagnosed in 14 of its 16 municipalities, Macapá being the city with the highest number of cases, with a recorded increase of 23.9% compared to 201833 Amapá [Estado], Superintendência de Vigilância em Saúde, Unidade de Doenças Transmissíveis. Bol Epidemiol Tuberculose. 2020 [cited 2021 Mar 15]. Available from: https://editor.amapa.gov.br/arquivos_portais/publicacoes/SVS_c20f41e26fed90da418341d2d2135a3a.pdf
https://editor.amapa.gov.br/arquivos_por...
.

For its eradication, wider-ranging actions are necessary, such as the active search for respiratory symptoms aiming at the early diagnosis of cases and the rapid start of treatment, to prevent the tuberculosis-infected person from remaining in the bacilliferous phase and transmitting the disease to other people and communities44 Ministério da Saúde (BR), Secretaria 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. 2. ed. atual. Brasília, DF: 2019 [cited 2021 Mar 15]. Available from: https://bvsms.saude.gov.br/bvs/publicacoes/manual_recomendacoes_controle_tuberculose_brasil_2_ed.pdf
https://bvsms.saude.gov.br/bvs/publicaco...
.

Regarding diagnosis, the rapid molecular test for tuberculosis (RMT-TB), performed by using the GeneXpert® MTB/IfR system, consists of an automated examination, which reduces the time of diagnosis result to a maximum of 2 hours and saves on costs, including those from hospitalizations resulting from the disease evolution and the late diagnosis. Hospitalization costs due to tuberculosis increase by up to one hundred times in relation to those of outpatient treatment55 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde, Departamento de Vigilância das Doenças Transmissíveis. Rede de Teste Rápido para Tuberculose no Brasil: primeiro ano da implantação. Brasília, DF; 2015 [cited 2021 Mar 15]. Available from: http://portalarquivos.saude.gov.br/images/pdf/2016/janeiro/19/rtr-tb-15jan16-isbn-web.pdf
http://portalarquivos.saude.gov.br/image...
.

With faster and highly sensitive and specific – respectively 98% and 88%55 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde, Departamento de Vigilância das Doenças Transmissíveis. Rede de Teste Rápido para Tuberculose no Brasil: primeiro ano da implantação. Brasília, DF; 2015 [cited 2021 Mar 15]. Available from: http://portalarquivos.saude.gov.br/images/pdf/2016/janeiro/19/rtr-tb-15jan16-isbn-web.pdf
http://portalarquivos.saude.gov.br/image...
– results, treatment can be started early, without the need to wait for confirmation by culture (which can take up to 60 days).

The WHO validated the RMT-TB results in 2010 after which the diagnostic technology was recommended and adopted in several health systems worldwide, including in Brazil, where the National Commission for the Incorporation of Technologies of the Unified Health System (CONITEC-SUS) and the National Health Surveillance Agency (ANVISA) approved the use of the test in 2013. The Ministry of Health acquired 160 Cepheid laboratory equipment and strategically distributed them throughout the country55 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde, Departamento de Vigilância das Doenças Transmissíveis. Rede de Teste Rápido para Tuberculose no Brasil: primeiro ano da implantação. Brasília, DF; 2015 [cited 2021 Mar 15]. Available from: http://portalarquivos.saude.gov.br/images/pdf/2016/janeiro/19/rtr-tb-15jan16-isbn-web.pdf
http://portalarquivos.saude.gov.br/image...
.

We hypothesized that the RMT-TB, due to its high sensitivity, may have modified the number of diagnosed cases and contributed to the elucidation of cases, specifically paucibacillary cases, previously undetected by bacilloscopy. In other words, we investigated whether RMT-TB may have changed the epidemiological panorama of tuberculosis.

In addition to the increase in cases, another hypothesis of this study is the seasonal variation of the disease, a phenomenon that has been observed in studies conducted in several countries66 Fares A. Seasonality of tuberculosis. J Glob Infect Dis. 2011;3(1):46-55. https://doi.org/10.4103/0974-777X.77296
https://doi.org/10.4103/0974-777X.77296...
, but little explored in Brazil. Another study77 Fundação de Vigilância em Saúde do Amazonas (BR). Situação epidemiológica da Síndrome Respiratória Aguda Grave no Estado do Amazonas. Bol Epidemiol. 2021 [cited 2021 Mar 15];8(11). Available from: https://www.fvs.am.gov.br/media/publicacao/Boletim_Epidemiol%C3%B3gico_N%C2%BA_11.pdf
https://www.fvs.am.gov.br/media/publicac...
showed that respiratory diseases increased in the Amazon in the winter periods, when rainfall in the region increases, therefore, the hypothesis arises that tuberculosis rates may also increase.

Considering the above, our study sought to evaluate the temporal trend of tuberculosis incidence after the implementation of RMT-TB, identify whether it presents seasonal variation and classify the territory according to case density and risk areas for the disease in a municipality of eastern Amazonia.

METHODS

Study Design

Ecological study88 Rothman KJ, Greenland S, Lash TL. Modern epidemiology. 3. ed. Philadelphia, PA: Lippincott Williams & Wilkins; 2008. conducted in Macapá, capital of the Brazillian state of Amapá, which is part of the Eastern Amazon and borders French Guiana.

Place of Study

The city of Macapá has an area of approximately 6,563.849 km², a population density of 62.14hab/km² and an estimated population of 512,902 people in 2020. It is administratively divided into four districts: Northern, Western, Central and Southern99 Instituto Brasileiro de Geografia e Estatística. Censo 2010: cidades: Macapá. Rio de Janeiro: IBGE; 2016 [cited 2021 Mar 15]. Available from: https://cidades.ibge.gov.br/brasil/ap/macapa/panorama
https://cidades.ibge.gov.br/brasil/ap/ma...
.

The diagnosis of tuberculosis cases in Macapá follows the protocol of the Ministry of Health55 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde, Departamento de Vigilância das Doenças Transmissíveis. Rede de Teste Rápido para Tuberculose no Brasil: primeiro ano da implantação. Brasília, DF; 2015 [cited 2021 Mar 15]. Available from: http://portalarquivos.saude.gov.br/images/pdf/2016/janeiro/19/rtr-tb-15jan16-isbn-web.pdf
http://portalarquivos.saude.gov.br/image...
. Thus, all basic health units (UBS) receive suspected cases that are later referred to the state reference unit, called the Reference Center for Tropical Diseases (CRDT). It is a state-wide centralized service, which monitors cases and dispenses medicines.

The state also has the Central Laboratory (LACEN), which performs laboratory tests, including those for tuberculosis diagnosis. The only GeneXpert® MTB/RIF equipment, a reference for the diagnosis of tuberculosis, in Amapá is in the LACEN. We emphasize that the RMT-TB was implemented in September 2014, which was considered the cut-off year in the time series55 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde, Departamento de Vigilância das Doenças Transmissíveis. Rede de Teste Rápido para Tuberculose no Brasil: primeiro ano da implantação. Brasília, DF; 2015 [cited 2021 Mar 15]. Available from: http://portalarquivos.saude.gov.br/images/pdf/2016/janeiro/19/rtr-tb-15jan16-isbn-web.pdf
http://portalarquivos.saude.gov.br/image...
.

Population

The study's population consisted of tuberculosis cases reported in the Sistema de Informação de Agravos de Notificação (SINAN – Information System for Notifiable Diseases) between 2001 and 2017.

Analysis Plan

Initially, to characterize the profile of tuberculosis cases reported in the municipality under study, the absolute and relative frequencies of the variables presented in the notification form were calculated using the IBM SPSS Statistics software version 25.

To classify the temporal trend of tuberculosis incidence, the rate was calculated monthly, considering the absolute number of cases in the numerator and the population in the denominator, with multiplication factor per 1000 inhabitants. The incidence rates were converted into logarithms (log10) to stabilize variance over time1010 Antunes JLF, Cardoso MRA. Uso da análise de séries temporais em estudos epidemiológicos. Epidemiol Serv Saude. 2015;24(3):565-76. https://doi.org/10.5123/S1679-49742015000300024
https://doi.org/10.5123/S1679-4974201500...
. The graphs of case distribution and of estimated rates throughout the time series were made using the RStudio software.

The Prais-Winsten self-regression method was used in the STATA software to classify the tuberculosis incidence temporal trend into increasing, decreasing or stationary and, if increasing or decreasing, the monthly percent change (MPC) and its respective 95% confidence interval (95%CI) was calculated1010 Antunes JLF, Cardoso MRA. Uso da análise de séries temporais em estudos epidemiológicos. Epidemiol Serv Saude. 2015;24(3):565-76. https://doi.org/10.5123/S1679-49742015000300024
https://doi.org/10.5123/S1679-4974201500...
.

The interrupted time series (ITS) was used to verify whether the temporal trend of tuberculosis incidence changed after the implementation of RMT-TB and whether the disease presents seasonal variation. It is defined as the most effective resource to evaluate the impact of an intervention, allowing to verify whether there is immediate influence (change of level) and/or progressive influence (trend change) on the values of the series1111 Wagner AK, Soumerai SB, Zhang F, Ross-Degnan D. Segmented regression analysis of interrupted time series studies in medication use research. J Clin Pharm Ther. 2002;27(4):299-309. https://doi.org/10.1046/j.1365-2710.2002.00430.x
https://doi.org/10.1046/j.1365-2710.2002...
.

For the application of the ITS, the previously calculated logarithmic monthly incidence rates (log10) were considered. The change in level was called “intervention” and the progressive change, “post-intervention”1111 Wagner AK, Soumerai SB, Zhang F, Ross-Degnan D. Segmented regression analysis of interrupted time series studies in medication use research. J Clin Pharm Ther. 2002;27(4):299-309. https://doi.org/10.1046/j.1365-2710.2002.00430.x
https://doi.org/10.1046/j.1365-2710.2002...
. Trigonometric functions sine and cosine were also incorporated into the analyses to verify whether tuberculosis presents seasonal variation1010 Antunes JLF, Cardoso MRA. Uso da análise de séries temporais em estudos epidemiológicos. Epidemiol Serv Saude. 2015;24(3):565-76. https://doi.org/10.5123/S1679-49742015000300024
https://doi.org/10.5123/S1679-4974201500...
. STATA software version 14 was used for ITS.

In the first stage, the inclusion criteria was all reported cases of tuberculosis from residents in the urban area of Macapá. For cases notified more than once in the system, the most up to date notification was considered. Note that, in the analyzed period, all cases listed in the study were from new notifications in SINAN, therefore, the incidence was estimated.

The Kernel density estimator was used to classify the density of cases. Initially, performing the geocoding of the cases by the geographic coordinates (latitude and longitude) of the residential addresses of the reported cases – obtained using the Google Earth Pro tool – was necessary.

The georeferencing of the cases included in the study was performed by the ArcGis 10.5 software, considering the urban census tracts of the municipality of Macapá as an analysis unit and the residing population by census tract, as a configuration for analysis standardization.

Note that, in the spatial analysis phase, cases of homeless people, individuals deprived of liberty, indigenous people or people whose notification address was from municipal agencies (hospitals or health units) were excluded due to the impossibility of obtaining the geographic coordinates of residential addresses.

The Kernel intensity estimator was used from the ArcGis 10.5 software. It is an exploratory interpolation method that defines circular areas of influence around occurrence points of a phenomenon, generating a surface density to identify vulnerable areas1212 Environmental Systems Research Institute. ArcGIS for Desktop. How Kernel Density works. West Redlands, CA: Esri; c2016 [cited 2021 Mar 15]. Available from: https://desktop.arcgis.com/en/arcmap/10.3/tools/spatial-analyst-toolbox/how-kernel-density-works.htm
https://desktop.arcgis.com/en/arcmap/10....
. Thus, considering a radius of 1,000 meters1313 Oliveira U, Brescovit AD, Santos AJ. Delimiting areas of endemism through Kernel Interpolation. PLoS One. 2015;10(1):e0116673. https://doi.org/10.1371/journal.pone.0116673
https://doi.org/10.1371/journal.pone.011...
,142 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde. Bol Epidemiol Tuberculose 2020. 2020 [cited 2021 Mar 15].];N° Espec:1-40. Available from: http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
http://www.aids.gov.br/pt-br/pub/2020/bo...
, thematic maps of the tuberculosis case density distribution were also generated in the ArcGIS 10.5 software.

Finally, to identify the risk areas for tuberculosis in Macapá, the scan statistic1515 Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Stat Med. 1995;14(8):799-810. https://doi.org/10.1002/sim.4780140809
https://doi.org/10.1002/sim.4780140809...
was used, in which a circle of variable radius around the centroide of each unit of analysis, which were the census tracts in this study, is used to search for spatial clusters. In this technique, the number of observed and expected cases within each circle is calculated beforehand, until all centroides are tested. When the observed value in the area inside the circle is greater or smaller than expected, the circle is called a cluster1616 Kulldorff M. SaTScanTM manual do usuário para versão 9.4: versão do original traduzido para o Português. Pelini, ACG tradutor. São Paulo; 2016 [cited 2021 Mar 15]. Available from: https://www.satscan.org/SaTScan_TM_Manual_do_Usu%C3%A1rio_v9.4_Portugues.pdf
https://www.satscan.org/SaTScan_TM_Manua...
.

In these analyses, the following characteristics were adopted: Poisson discrete model; circular-shaped agglomerates; no geographical overlap of the agglomerates; 999 replications in Monte Carlo simulation; and exposed population size stipulated by the Gini coefficient, in which the number of cases is compared to the data of the base population and the expected number of cases in each census tract is proportional to the size of the population at risk1616 Kulldorff M. SaTScanTM manual do usuário para versão 9.4: versão do original traduzido para o Português. Pelini, ACG tradutor. São Paulo; 2016 [cited 2021 Mar 15]. Available from: https://www.satscan.org/SaTScan_TM_Manual_do_Usu%C3%A1rio_v9.4_Portugues.pdf
https://www.satscan.org/SaTScan_TM_Manua...
,172 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde. Bol Epidemiol Tuberculose 2020. 2020 [cited 2021 Mar 15].];N° Espec:1-40. Available from: http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
http://www.aids.gov.br/pt-br/pub/2020/bo...
. The relative risk (RR) and confidence interval (95%CI) of each cluster were also calculated. Note that clusters with p < 0.05 were considered statistically significant.

The study was approved by the Research Ethics Committee of the Escola de Enfermagem de Ribeirão Preto (CAAE No. 23043019.2.0000.5393).

RESULTS

Between 2001 and 2017, 1,730 cases of tuberculosis were reported in Macapá, with a minimum age of one year and a maximum of 89 years.

Most cases were female (59.3%), aged between 31 and 59 years (47.6%), mixed race (67.9%) and incomplete elementary school education (46.6%). The pulmonary form is predominant (85.1%).

Table 1 shows the high frequency of unanswered and/or ignored information regarding schooling (16.4%), TB-HIV co-infection (50%), TB-Diabetes (41.8%), alcoholism (39.7%) and mental illness (41.7%), despite the relevance of this information for the clinical management of cases and conducts.

Table 1
Sociodemographic and epidemiological clinical profile of tuberculosis (TB) cases reported in Macapá, state of Amapá, Brazil (2001 to 2017).

Figure 1 shows the monthly time series of the number of cases and incidence of tuberculosis from 2001 to 2017, in which positive and negative peaks occur throughout the study period. We concluded that the series begins with higher numbers, suffers a decrease and ends in an increase.

Figure 1
Time series of the number of cases and tuberculosis incidence in Macapá, Amapá, Brazil (2001 to 2017).

Table 2 shows the results of the Prais-Winsten and ITS techniques. The temporal trend of tuberculosis incidence in Macapá was classified as decreasing at a −0.27%/month rate (95%CI −0.41 to −0.13).

Table 2
Temporal trend, impact of diagnosis by rapid molecular testing and seasonal variation in the incidence of tuberculosis cases, Macapá, Amapá, Brazil (2001 to 2017).

We saw no change of level (intervention) in the time series after the implementation of RMT-TB. However, we classified the period after the implantation of RMT-TB as increasing, with a 2.09%/month rate (95%CI 3.27 to 0.92), indicating that the incidence of tuberculosis increased progressively after the implementation of the test.

Regarding seasonality, the disease in Macapá showed a growth of 13.7%/month (95%CI 23.87 to 4.71) in the months between December and June (sine), a period that coincides with the rains of the so-called Amazonian winter. We saw a decrease of 9.21%/month (95%CI −1.37 to −16.63) in the remaining months (cosine), which are warmer and dryer seasons.

Of the 1,730 tuberculosis cases identified, 59 (3.41%) were excluded for lack of address. Of the remaining 1,671 cases, 1,475 (88.2%) had their geographic coordinates identified, thus georeferenced.

With the Kernel estimator, we observed areas classified as very high density of cases in the Central (neighborhoods of Laguinho and Perpétuo Socorro) and Northern (neighborhoods of Pacoval and Cidade Nova) districts, ranging from 70.67 to 88.33 cases/km². Figure 2 shows that areas classified as high density of cases are in the Southern (neighborhood of Buritizal), Central (neighborhoods of Laguinho and Perpétuo Socorro) and Northern (neighborhoods of Pacoval and Cidade Nova) districts, ranging from 53.01 to 70.66 cases/km2.

Figure 2
Geographic location of tuberculosis cases reported in Macapá-AP, Brazil (2001 to 2017).

Using the scan statistics, we identified four spatial clusters (SC) at risk for tuberculosis in Macapá (Figure 3), considering the parameter of 40% of the exposed population, according to the Gini index.

Figure 3
Areas of spatial risk for tuberculosis in Macapá, Amapá, Brazil (2001 to 2017).

The SC1, considered a protection cluster (RR = 0.07; 95%CI 0.01 to 0.48), consisted of six census tracts (p = 0.03) of the Southern district, with a population of 3,799 inhabitants, one observed case and 14 expected cases.

The SC2, another protection cluster (RR = 0.23; 95%CI 0.10 to 0.52), consisted of 10 census tracts (p = 0.03) of the Southern district, with a population of 6,578 inhabitants, six observed cases and 25 expected cases.

The SC3, also considered protective (RR = 0.36; 95%CI 0.20 to 0.59), consisted of 11 census tracts of the Western district (p = 0.04), with 10,172 inhabitants, 14 observed cases and 38 expected cases.

Finally, SC4 (p < 0.01), considered a high-risk cluster for tuberculosis (RR = 1.47; 95%CI 1.39 to 1.72), consisted of 146 census tracts of the Northern, Southern and Central districts, with a population of 122,356 inhabitants, 620 observed cases and 487 expected cases.

DISCUSSION

We managed to trace the sociodemographic and clinical-epidemiological profile of tuberculosis cases in Macapá, Amapá. By the calculation of the monthly incidence in a 16-year time series, we could classify the temporal trend and identify changes in the series after RMT-TB implementation, besides showing that tuberculosis has seasonal variation. Finally, we could classify areas according to the case density and identify risk and protection areas in Macapá, Amapá.

Most people affected by tuberculosis in Macapá were female, which differs from the findings of the literature that shows cases mostly in males – explained by the fact that men tend to minimize their self-care and be more exposed to risk factors for the disease1818 Hoshino H, Uchimura K, Yamauchi Y. [Comparison of TB incidence of young and middle age groups between urban/suburban prefectures and other prefectures]. Kekkaku. 2009;84(1):1-8. Japanese. https://doi.org/10.11400/kekkaku.84.1
https://doi.org/10.11400/kekkaku.84.1...
.

The predominance of women in this case may be related to the fact that women are heads of families, undervalued in the labor market and often work triple hours, which are stress factors that can be considered a risk for illness1818 Hoshino H, Uchimura K, Yamauchi Y. [Comparison of TB incidence of young and middle age groups between urban/suburban prefectures and other prefectures]. Kekkaku. 2009;84(1):1-8. Japanese. https://doi.org/10.11400/kekkaku.84.1
https://doi.org/10.11400/kekkaku.84.1...
. The gender inequality, inherited from local cultural issues of the riverside peoples in which women are seen as adjunct and the man as a provider, persists to this day and is another stress factor.

Most cases were detected in adults in an economically active age group, a phase of life in which individuals tend to face greater agglomerations, due to work and daily activities, increasing the chances of illness1919 Fontes GJF, Silva TG, Sousa JCM, Feitosa ANA, Silva ML, Bezerra ALD, et al. Perfil epidemiológico da tuberculose no Brasil no período de 2012 a 2016. Rev Bras Educ Saude. 2019;9(1):19-26. https://doi.org/10.18378/rebes.v9i1.6376
https://doi.org/10.18378/rebes.v9i1.6376...
.

Mixed race people were more prone to the disease in the study, which can be justified by being the predominant race/color in the municipality. Almost 52% of the Macapá population declares itself as mixed race and, in this study, the self-declared as mixed race are 67.9%, which would justify the high rates of tuberculosis in this population.

Most reported cases have incomplete elementary school, which is an indirect indicator of social conditions caused by the population's social context. This data represents an obstacle to TB control, since these individuals may have difficulties in understanding the disease and the guidelines for its treatment and prevention2020 Santos JN, Sales CMM, Prado TN, Maciel EL. Fatores associados à cura no tratamento da tuberculose no estado do Rio de Janeiro, 2011-2014. Epidemiol Serv Saude. 2018;27(3):e2017464. https://doi.org/10.5123/s1679-49742018000300015
https://doi.org/10.5123/s1679-4974201800...
.

Considering this panorama, the health professional must be able to clarify these situations in complementary ways, seeking ways to increase the effectiveness of TB control and promote a holistic view of health care2121 Coriolano-Marinus MWL, Queiroga BAM, Ruiz-Moreno L, Lima LS. Comunicação nas práticas em saúde: revisão integrativa da literatura. Saude Soc. 2014;23(4):1356-69. https://doi.org/10.1590/S0104-12902014000400019
https://doi.org/10.1590/S0104-1290201400...
.

We observed the predominance of the pulmonary clinical form of the disease, which was expected, since it is its most common form11 World Health Organization. Global Tuberculosis Report 2020. Geneva (CH): WHO; 2020 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/336069/9789240013131-eng.pdf
https://apps.who.int/iris/bitstream/hand...
,22 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde. Bol Epidemiol Tuberculose 2020. 2020 [cited 2021 Mar 15].];N° Espec:1-40. Available from: http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
http://www.aids.gov.br/pt-br/pub/2020/bo...
. Most cases reached a cure (75.3%), however, this rate was far below the WHO recommendation of 85%11 World Health Organization. Global Tuberculosis Report 2020. Geneva (CH): WHO; 2020 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/336069/9789240013131-eng.pdf
https://apps.who.int/iris/bitstream/hand...
. Our findings point to the high number of dropouts (12.5%), above that established by the WHO (5%), which is a predictor for the development of resistant tuberculosis, whose treatment is longer and more costly for the health system2020 Santos JN, Sales CMM, Prado TN, Maciel EL. Fatores associados à cura no tratamento da tuberculose no estado do Rio de Janeiro, 2011-2014. Epidemiol Serv Saude. 2018;27(3):e2017464. https://doi.org/10.5123/s1679-49742018000300015
https://doi.org/10.5123/s1679-4974201800...
.

Similarly, we observed a high percentage of uninformed variables, such as the comorbidities HIV, diabetes, alcoholism and mental illnesses; which is information of the utmost relevance for the clinical management of cases. Some assumptions for non-records are the non-prioritization of this information in the clinic, the incipient qualification of the teams on the importance of health surveillance based on the quality of the data, and the lack of evaluation and feedback to the teams about the data produced in the units in which the diagnoses were made.

In the studied period, the incidence of tuberculosis was decreasing, a pattern that follows the worldwide trend11 World Health Organization. Global Tuberculosis Report 2020. Geneva (CH): WHO; 2020 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/336069/9789240013131-eng.pdf
https://apps.who.int/iris/bitstream/hand...
,222 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde. Bol Epidemiol Tuberculose 2020. 2020 [cited 2021 Mar 15].];N° Espec:1-40. Available from: http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
http://www.aids.gov.br/pt-br/pub/2020/bo...
. However, this decline occurs at a speed below the WHO expectations11 World Health Organization. Global Tuberculosis Report 2020. Geneva (CH): WHO; 2020 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/336069/9789240013131-eng.pdf
https://apps.who.int/iris/bitstream/hand...
, thus it is unlikely that the disease will be fully erradicated by 2050.

With ITS, we identified no change of level after the implementation of RMT-TB, however, the period of “post-intervention” was classified as increasing, confirming the hypothesis that the test, thanks to its high sensitivity2323 Boehme CC, Nabeta P, Hillemann D, Nicol MP, Shenai S, Krapp F, et al. Rapid molecular detection of tuberculosis and rifampin resistance. N Engl J Med. 2010;363(11):1005-15. https://doi.org/10.1056/NEJMoa0907847
https://doi.org/10.1056/NEJMoa0907847...
252 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde. Bol Epidemiol Tuberculose 2020. 2020 [cited 2021 Mar 15].];N° Espec:1-40. Available from: http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
http://www.aids.gov.br/pt-br/pub/2020/bo...
, may have influenced the increase of cases in the municipality.

Bacilloscopy remains in the Unified Health System (SUS) as a method for monthly follow-up of treatment and RMT-TB is performed only for initial diagnosis. The culture, considered the gold standard for the diagnosis, takes around 45 to 60 days, due to slow replication of the bacillus, which impedes its use as the only strategy for early diagnosis and timely therapy44 Ministério da Saúde (BR), Secretaria 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. 2. ed. atual. Brasília, DF: 2019 [cited 2021 Mar 15]. Available from: https://bvsms.saude.gov.br/bvs/publicacoes/manual_recomendacoes_controle_tuberculose_brasil_2_ed.pdf
https://bvsms.saude.gov.br/bvs/publicaco...
,52 Ministério da Saúde (BR), Secretaria de Vigilância em Saúde. Bol Epidemiol Tuberculose 2020. 2020 [cited 2021 Mar 15].];N° Espec:1-40. Available from: http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
http://www.aids.gov.br/pt-br/pub/2020/bo...
.

However, the use of RMT-TB does not dispense with the culture, which should always be considered, so that its data are compared with those of RMT-TB (subject to false positives regarding the sensitivity or rifampicin resistance) and drug-resistant strains are identified2424 Lima TM, Belotti NCU, Nardi SMT, Pedro HSP. Teste rápido molecular GeneXpert MTB/RIF para diagnóstico da tuberculose. Rev Pan Amaz Saude 2017;8(2):67-78..

The main strategy for tuberculosis control is the early detection of cases – especially bacilliferous pulmonary cases, which have greater epidemiological importance – aiming at maintaining the transmission chain of the disease, followed by adequate treatment and complete cure2626 Spagnolo LML, Tomberg JO, Vieira DA, Gonzales RIC. Detecção da tuberculose: fluxo dos sintomáticos respiratórios e resultados alcançados. Rev Bras Enferm. 2018;71(5):2543-51. https://doi.org/10.1590/0034-7167-2017-0457
https://doi.org/10.1590/0034-7167-2017-0...
. According to the WHO2727 World Health Organization. Global Tuberculosis Report 2016. Geneva (CH): WHO; 2019 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/250441/9789241565394-eng.pdf?sequence=1&isAllowed=y
https://apps.who.int/iris/bitstream/hand...
, about 43 million deaths were avoided between the years 2000 and 2014 by the use of early diagnosis and appropriate treatment. Thus, strategies and public policies with this purpose can, in the medium or long term, reduce disease rates. This means reliable means of diagnosis are essential.

Our study also evidenced the seasonality of TB, with significant oscillations of cases and incidence throughout the year, thus confirming one of our initial hypotheses. The findings showed that the incidence increased between December and June, coinciding with the rainy season. Respiratory diseases in the Eastern Amazon77 Fundação de Vigilância em Saúde do Amazonas (BR). Situação epidemiológica da Síndrome Respiratória Aguda Grave no Estado do Amazonas. Bol Epidemiol. 2021 [cited 2021 Mar 15];8(11). Available from: https://www.fvs.am.gov.br/media/publicacao/Boletim_Epidemiol%C3%B3gico_N%C2%BA_11.pdf
https://www.fvs.am.gov.br/media/publicac...
showed this behavior, but it had not yet been evidenced for TB.

This can be attributed to the need for people to stay on overcrowded places with little ventilation in these periods. In warmer times, on the other hand, the temporal trend decreased, since people seek more airy places, which may decrease the contagion with the bacillus causing TB.

The use of Kernel and scan statistics made possible to identify that the highest densities of cases and risk areas (SC4) are in the Central, Northern and Southern districts and in the regions that concentrate the lowest Human Development Indexes of the municipality, which, added to the existence of bridge areas (stilt houses), to the absence of sanitation, to the clusters of population, and to the large number of informal workers corroborate the spread of the disease.

Near those locations, other areas considered as protection (SC1, SC2 and SC3) were identified, which leads to the assumption that these regions have undiagnosed cases. Therefore, underreporting is a possibility, since social vulnerability is also notable in these regions. This result should be cautiously analyzed and serves as a warning to health surveillance teams.

The use of strategies such as income transfer and/or compensatory policies, to alleviate social inequality and its deleterious effects on the vulnerable population is critical2828 Moreira ASR, Kinski AL, Carvalho ACC. Social determinants of health and catastrophic costs associated with the diagnosis and treatment of tuberculosis. J Bras Pneumol. 2020;46(5):e20200015. https://doi.org/10.36416/1806-3756/e20200015
https://doi.org/10.36416/1806-3756/e2020...
.

Among the limitations of this study, the ecological fallacy stands out2929 Veras RP, Almeida Filho N, Barreto ML, Veras RP, Barata RB. Teoria epidemiológica hoje: fundamentos, interfaces, tendências. Rio de Janeiro: Editora FIOCRUZ; 1998 [cited 2021 Mar 15]. (Séries Epidemiológicas, N° 2). Available from: http://books.scielo.org/id/5btwk
http://books.scielo.org/id/5btwk...
, meaning the data are analyzed at the aggregate level and the results cannot be interpreted at the individual level. Also note the use of secondary data, which may contain typing errors and missing information, which may interfere with the analyses3030 Melo GBT, Valongueiro S. Incompleteness of Mortality Information System records on deaths from external causes in Pernambuco, Brazil, 2000-2002 and 2008-2010. Epidemiol Serv Saude. 2015;24(4):651-60. https://doi.org/10.5123/S1679-4974201500040000
https://doi.org/10.5123/S1679-4974201500...
. Another limitation was the lack of the socioeconomic variable income, which could be used as a covariate in the adjustment of the temporal model.

The use of two distinct approaches – time series and spatial analyses – to understand the behavior of tuberculosis in the eastern Amazon makes this study relevant. The evidence of the influence of RMT-TB on TB incidence rates and of its seasonal behavior in the Eastern Amazon adds to its originality, making this a possible reference for future studies.

This will also help orient public policies, given that differences in of TB load between the territories are evidenced, with social vulnerability as one of its probable explanations.

  • Funding:Fundação de Amparo à Pesquisa do Estado de São Paulo (Fapesp - Doctorate grant Process No. 2018/03700-7). Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq Research Productivity Grant – Process 304483/2018-4 – PQ Modality – Level 1C). Coordenação de Aperfeicoamento de Pessoal de Nível Superior (Capes – Financing Code 001 – PROEX Program – Aid No. 0524/2018, Process No.: 23038.003643/2018-11).

REFERENCES

  • 1
    World Health Organization. Global Tuberculosis Report 2020. Geneva (CH): WHO; 2020 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/336069/9789240013131-eng.pdf
    » https://apps.who.int/iris/bitstream/handle/10665/336069/9789240013131-eng.pdf
  • 2
    Ministério da Saúde (BR), Secretaria de Vigilância em Saúde. Bol Epidemiol Tuberculose 2020. 2020 [cited 2021 Mar 15].];N° Espec:1-40. Available from: http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
    » http://www.aids.gov.br/pt-br/pub/2020/boletim-epidemiologico-de-turbeculose-2020
  • 3
    Amapá [Estado], Superintendência de Vigilância em Saúde, Unidade de Doenças Transmissíveis. Bol Epidemiol Tuberculose. 2020 [cited 2021 Mar 15]. Available from: https://editor.amapa.gov.br/arquivos_portais/publicacoes/SVS_c20f41e26fed90da418341d2d2135a3a.pdf
    » https://editor.amapa.gov.br/arquivos_portais/publicacoes/SVS_c20f41e26fed90da418341d2d2135a3a.pdf
  • 4
    Ministério da Saúde (BR), Secretaria 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. 2. ed. atual. Brasília, DF: 2019 [cited 2021 Mar 15]. Available from: https://bvsms.saude.gov.br/bvs/publicacoes/manual_recomendacoes_controle_tuberculose_brasil_2_ed.pdf
    » https://bvsms.saude.gov.br/bvs/publicacoes/manual_recomendacoes_controle_tuberculose_brasil_2_ed.pdf
  • 5
    Ministério da Saúde (BR), Secretaria de Vigilância em Saúde, Departamento de Vigilância das Doenças Transmissíveis. Rede de Teste Rápido para Tuberculose no Brasil: primeiro ano da implantação. Brasília, DF; 2015 [cited 2021 Mar 15]. Available from: http://portalarquivos.saude.gov.br/images/pdf/2016/janeiro/19/rtr-tb-15jan16-isbn-web.pdf
    » http://portalarquivos.saude.gov.br/images/pdf/2016/janeiro/19/rtr-tb-15jan16-isbn-web.pdf
  • 6
    Fares A. Seasonality of tuberculosis. J Glob Infect Dis. 2011;3(1):46-55. https://doi.org/10.4103/0974-777X.77296
    » https://doi.org/10.4103/0974-777X.77296
  • 7
    Fundação de Vigilância em Saúde do Amazonas (BR). Situação epidemiológica da Síndrome Respiratória Aguda Grave no Estado do Amazonas. Bol Epidemiol. 2021 [cited 2021 Mar 15];8(11). Available from: https://www.fvs.am.gov.br/media/publicacao/Boletim_Epidemiol%C3%B3gico_N%C2%BA_11.pdf
    » https://www.fvs.am.gov.br/media/publicacao/Boletim_Epidemiol%C3%B3gico_N%C2%BA_11.pdf
  • 8
    Rothman KJ, Greenland S, Lash TL. Modern epidemiology. 3. ed. Philadelphia, PA: Lippincott Williams & Wilkins; 2008.
  • 9
    Instituto Brasileiro de Geografia e Estatística. Censo 2010: cidades: Macapá. Rio de Janeiro: IBGE; 2016 [cited 2021 Mar 15]. Available from: https://cidades.ibge.gov.br/brasil/ap/macapa/panorama
    » https://cidades.ibge.gov.br/brasil/ap/macapa/panorama
  • 10
    Antunes JLF, Cardoso MRA. Uso da análise de séries temporais em estudos epidemiológicos. Epidemiol Serv Saude. 2015;24(3):565-76. https://doi.org/10.5123/S1679-49742015000300024
    » https://doi.org/10.5123/S1679-49742015000300024
  • 11
    Wagner AK, Soumerai SB, Zhang F, Ross-Degnan D. Segmented regression analysis of interrupted time series studies in medication use research. J Clin Pharm Ther. 2002;27(4):299-309. https://doi.org/10.1046/j.1365-2710.2002.00430.x
    » https://doi.org/10.1046/j.1365-2710.2002.00430.x
  • 12
    Environmental Systems Research Institute. ArcGIS for Desktop. How Kernel Density works. West Redlands, CA: Esri; c2016 [cited 2021 Mar 15]. Available from: https://desktop.arcgis.com/en/arcmap/10.3/tools/spatial-analyst-toolbox/how-kernel-density-works.htm
    » https://desktop.arcgis.com/en/arcmap/10.3/tools/spatial-analyst-toolbox/how-kernel-density-works.htm
  • 13
    Oliveira U, Brescovit AD, Santos AJ. Delimiting areas of endemism through Kernel Interpolation. PLoS One. 2015;10(1):e0116673. https://doi.org/10.1371/journal.pone.0116673
    » https://doi.org/10.1371/journal.pone.0116673
  • 14
    Yamamura M, Freitas IM, Santos-Neto M, Chiaravalloti-Neto F, Popolin MAP, Arroyo LH, et al. Análise espacial das internações evitáveis por tuberculose em Ribeirão Preto, SP (2006-2012). Rev Saude Publica. 2016;50:20. https://doi.org/10.1590/S1518-8787.2016050006049
    » https://doi.org/10.1590/S1518-8787.2016050006049
  • 15
    Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Stat Med. 1995;14(8):799-810. https://doi.org/10.1002/sim.4780140809
    » https://doi.org/10.1002/sim.4780140809
  • 16
    Kulldorff M. SaTScanTM manual do usuário para versão 9.4: versão do original traduzido para o Português. Pelini, ACG tradutor. São Paulo; 2016 [cited 2021 Mar 15]. Available from: https://www.satscan.org/SaTScan_TM_Manual_do_Usu%C3%A1rio_v9.4_Portugues.pdf
    » https://www.satscan.org/SaTScan_TM_Manual_do_Usu%C3%A1rio_v9.4_Portugues.pdf
  • 17
    Han J, Zhu L, Kulldorff M, Hostovich S, Stinchcomb DG, Tatalovich Z, et al. Using Gini coefficient to determining optimal cluster reporting sizes for spatial scan statistics. Int J Health Geogr. 2016;15(1):27. https://doi.org/10.1186/s12942-016-0056-6
    » https://doi.org/10.1186/s12942-016-0056-6
  • 18
    Hoshino H, Uchimura K, Yamauchi Y. [Comparison of TB incidence of young and middle age groups between urban/suburban prefectures and other prefectures]. Kekkaku. 2009;84(1):1-8. Japanese. https://doi.org/10.11400/kekkaku.84.1
    » https://doi.org/10.11400/kekkaku.84.1
  • 19
    Fontes GJF, Silva TG, Sousa JCM, Feitosa ANA, Silva ML, Bezerra ALD, et al. Perfil epidemiológico da tuberculose no Brasil no período de 2012 a 2016. Rev Bras Educ Saude. 2019;9(1):19-26. https://doi.org/10.18378/rebes.v9i1.6376
    » https://doi.org/10.18378/rebes.v9i1.6376
  • 20
    Santos JN, Sales CMM, Prado TN, Maciel EL. Fatores associados à cura no tratamento da tuberculose no estado do Rio de Janeiro, 2011-2014. Epidemiol Serv Saude. 2018;27(3):e2017464. https://doi.org/10.5123/s1679-49742018000300015
    » https://doi.org/10.5123/s1679-49742018000300015
  • 21
    Coriolano-Marinus MWL, Queiroga BAM, Ruiz-Moreno L, Lima LS. Comunicação nas práticas em saúde: revisão integrativa da literatura. Saude Soc. 2014;23(4):1356-69. https://doi.org/10.1590/S0104-12902014000400019
    » https://doi.org/10.1590/S0104-12902014000400019
  • 22
    World Health Organization. Global Tuberculosis Report 2019. Geneva (CH): WHO; 2019 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/329368/9789241565714-eng.pdf
    » https://apps.who.int/iris/bitstream/handle/10665/329368/9789241565714-eng.pdf
  • 23
    Boehme CC, Nabeta P, Hillemann D, Nicol MP, Shenai S, Krapp F, et al. Rapid molecular detection of tuberculosis and rifampin resistance. N Engl J Med. 2010;363(11):1005-15. https://doi.org/10.1056/NEJMoa0907847
    » https://doi.org/10.1056/NEJMoa0907847
  • 24
    Lima TM, Belotti NCU, Nardi SMT, Pedro HSP. Teste rápido molecular GeneXpert MTB/RIF para diagnóstico da tuberculose. Rev Pan Amaz Saude 2017;8(2):67-78.
  • 25
    Pandey P, Pant ND, Rijal KR, Shrestha B, Kattel S, Banjara MR, et al. Diagnostic accuracy of GeneXpert MTB/RIF assay in comparison to conventional drug susceptibility testing method for the diagnosis of multidrug-resistant tuberculosis. PLoS One. 2017;12(1):e0169798. https://doi.org/10.1371/journal.pone.0169798
    » https://doi.org/10.1371/journal.pone.0169798
  • 26
    Spagnolo LML, Tomberg JO, Vieira DA, Gonzales RIC. Detecção da tuberculose: fluxo dos sintomáticos respiratórios e resultados alcançados. Rev Bras Enferm. 2018;71(5):2543-51. https://doi.org/10.1590/0034-7167-2017-0457
    » https://doi.org/10.1590/0034-7167-2017-0457
  • 27
    World Health Organization. Global Tuberculosis Report 2016. Geneva (CH): WHO; 2019 [cited 2021 Mar 15]. Available from: https://apps.who.int/iris/bitstream/handle/10665/250441/9789241565394-eng.pdf?sequence=1&isAllowed=y
    » https://apps.who.int/iris/bitstream/handle/10665/250441/9789241565394-eng.pdf?sequence=1&isAllowed=y
  • 28
    Moreira ASR, Kinski AL, Carvalho ACC. Social determinants of health and catastrophic costs associated with the diagnosis and treatment of tuberculosis. J Bras Pneumol. 2020;46(5):e20200015. https://doi.org/10.36416/1806-3756/e20200015
    » https://doi.org/10.36416/1806-3756/e20200015
  • 29
    Veras RP, Almeida Filho N, Barreto ML, Veras RP, Barata RB. Teoria epidemiológica hoje: fundamentos, interfaces, tendências. Rio de Janeiro: Editora FIOCRUZ; 1998 [cited 2021 Mar 15]. (Séries Epidemiológicas, N° 2). Available from: http://books.scielo.org/id/5btwk
    » http://books.scielo.org/id/5btwk
  • 30
    Melo GBT, Valongueiro S. Incompleteness of Mortality Information System records on deaths from external causes in Pernambuco, Brazil, 2000-2002 and 2008-2010. Epidemiol Serv Saude. 2015;24(4):651-60. https://doi.org/10.5123/S1679-4974201500040000
    » https://doi.org/10.5123/S1679-4974201500040000

Publication Dates

  • Publication in this collection
    01 Dec 2021
  • Date of issue
    2021

History

  • Received
    21 Dec 2020
  • Accepted
    19 Mar 2021
Faculdade de Saúde Pública da Universidade de São Paulo São Paulo - SP - Brazil
E-mail: revsp@org.usp.br