Skip to main content

Risk-prone territories for spreading tuberculosis, temporal trends and their determinants in a high burden city from São Paulo State, Brazil



To identify risk-prone areas for the spread of tuberculosis, analyze spatial variation and temporal trends of the disease in these areas and identify their determinants in a high burden city.


An ecological study was carried out in Ribeirão Preto, São Paulo, Brazil. The population was composed of pulmonary tuberculosis cases reported in the Tuberculosis Patient Control System between 2006 and 2017. Seasonal Trend Decomposition using the Loess decomposition method was used. Spatial and spatiotemporal scanning statistics were applied to identify risk areas. Spatial Variation in Temporal Trends (SVTT) was used to detect risk-prone territories with changes in the temporal trend. Finally, Pearson's Chi-square test was performed to identify factors associated with the epidemiological situation in the municipality.


Between 2006 and 2017, 1760 cases of pulmonary tuberculosis were reported in the municipality. With spatial scanning, four groups of clusters were identified with relative risks (RR) from 0.19 to 0.52, 1.73, 2.07, and 2.68 to 2.72. With the space–time scan, four clusters were also identified with RR of 0.13 (2008–2013), 1.94 (2010–2015), 2.34 (2006 to 2011), and 2.84 (2014–2017). With the SVTT, a cluster was identified with RR 0.11, an internal time trend of growth (+ 0.09%/year), and an external time trend of decrease (− 0.06%/year). Finally, three risk factors and three protective factors that are associated with the epidemiological situation in the municipality were identified, being: race/brown color (OR: 1.26), without education (OR: 1.71), retired (OR: 1.35), 15 years or more of study (OR: 0.73), not having HIV (OR: 0.55) and not having diabetes (OR: 0.35).


The importance of using spatial analysis tools in identifying areas that should be prioritized for TB control is highlighted, and greater attention is necessary for individuals who fit the profile indicated as “at risk” for the disease.

Peer Review reports


Tuberculosis is one of the oldest diseases in the world and is still the cause of illness for millions of people each year. It is estimated that in 2020 there were approximately 10 million cases of tuberculosis in the world, with 56% of these cases affecting males, 33% females, and 11% children under 15 years of age [1]. In addition to these alarming statistics, the WHO still estimates 1.3 million deaths from tuberculosis and more than 214,000 deaths caused by co-infection with tuberculosis and human immunodeficiency virus (TB-HIV) [1].

Brazil remains among the 30 countries with a high burden for TB and TB-HIV co-infection, being considered a priority for disease control globally by the WHO [2]. In 2020, Brazil recorded about 67,000 new TB cases, presenting an incidence coefficient of 31.6 cases per 100 thousand inhabitants, also, in 2019, about 4500 deaths from the disease, was presented with a coefficient of mortality of 2.2 deaths per 100 thousand inhabitants [2].

It is known that the incidence of any disease changes over time and that the time trend also varies according to the geographic region. Thus, the importance of monitoring the emerging spatial patterns and temporal risk trends for tuberculosis is emphasized to provide additional information which helps prevent the disease, implement control measures and address new health risks within a geographical region [3].

In scientific literature, many studies have analyzed the behavior of the disease in a territory under study and sought to identify areas of risk. However, few studies have addressed the temporal behavior of tuberculosis and concomitantly placed variations in the temporal trends in each territory at risk. Therefore, this study aimed to identify risk-prone areas for the spread of tuberculosis, analyze spatial variation and temporal trends in these areas, and their determinants in a municipality with high endemicity due to the disease.


Research design and scenario

An ecological study [4] was carried out in Ribeirão Preto, a city in the interior of the state of São Paulo (SP)—Brazil. Located 314 km from the capital, Ribeirão Preto has an area of approximately 650 km2 and a high population density of 995.3 inhabitants/km2. It also had an estimated population of 711,825 inhabitants in 2020, of which 99.7% live in urban areas [5].

The unit of ecological analysis used in the study is the census sector. Ribeirão Preto is composed of 972 census sectors, of which 956 were considered in the present study because they are urban areas of the municipality.


The study population was composed of cases of pulmonary tuberculosis reported through the Tuberculosis Control System (TBWeb) [6] from 2006 to 2017. Data were obtained from the Municipal Tuberculosis Control Program of the Ribeirão Preto Municipal Secretariat. Only pulmonary tuberculosis cases were selected since it is the unique contagious form of the disease and with chances of spreading in the territory.

It was adopted as a selection criterion that the notification was carried out between 2006 and 2017, with only one registration per person, the most current registration being selected if there was more than one entry in the system and only residents in urban area of the city of Ribeirão Preto were included. It is noteworthy that only pulmonary tuberculosis records were considered so that extrapulmonary or concomitant forms (pulmonary and extrapulmonary forms together) were excluded as part of the exclusion criteria.

Analysis plan

Time series analysis

Initially, monthly time series of tuberculosis cases were constructed, spanning between January 2006 to December 2017. To verify the behavior of the time series over the study period and also its trend, the decomposition method called Seasonal Trend Decomposition using Loess (STL) was used, based on a locally weighted regression [7]. This analysis was performed using RStudio software through the forecast package [8].

Identification of clusters

The georeferencing of pulmonary tuberculosis cases were performed using the Google Earth Pro® software in order to obtain the geographical coordinates (latitude and longitude) of the residential addresses of the notified cases.

In order to identify areas at higher risk for pulmonary tuberculosis, the spatial analysis technique called scanning statistics, developed by Kulldorff and Nagarwalla [9] was used.

It is considered as a null hypothesis if there is no high or low-risk cluster, which means, the entire population has the same probability of contracting pulmonary tuberculosis, regardless of its location; while the alternative hypothesis predicts the existence of clusters that are areas in which the population would be more or less likely to contract the disease [10].

Unlike the purely spatial scan that is based on circles, in the space–time scan, cylinders are created around each centroid, incorporating time as a variable of interest, it is possible to verify the existence of clusters in a given area, and also prove that in a specific period of time, there was a greater or lesser proportion of cases when compared to the other areas that were analyzed [11].

Still referring to the analysis of cluster detection, the SVTT technique was also performed, which differs from the other analyses presented by calculating the temporal trend of the clusters [3]. The time trend is calculated inside and outside the scan circle, we call the internal temporal trend (ITT) the change in the time trend of the event within a cluster and the external temporal trend (ETT) is called the trend of all other areas that do not belong to the cluster in question. Therefore, what is statistically significant in this analysis are the temporal trends and not the cluster formation as in spatial and space–time scanning [3, 11].

SVTT is considered a null hypothesis when there is no difference in the temporal trends in the analyzed areas while we have as an alternative hypothesis the occurrence that the temporal trends are statistically different.

In addition, the relative risk (RR) and 95% confidence interval (95%CI) of each cluster was calculated, allowing the comparison of information in different areas, with the exception of SVTT, because, as explained above, it is emphasized that what is significant in this analysis are ITT and ETT, so that the RR of the identified cluster may not be within the CI. Clusters with p < 0.05 were considered statistically significant. The analyses were performed using SaTScan software 9.3, and thematic maps were created using ArcGis software 10.5.

Descriptive and association analysis

In order to identify factors associated with the epidemiological situation in the municipality, an exploratory analysis was carried out (absolute and relative frequencies), and then the association of these variables with the fact of living in a risk area was identified by means of Pearson's chi-square analysis (χ2) and for the variables that were statistically significant (p < 0.05), the Odds Ratio (OR) and 95%CI were calculated using the IBM SPSS version 25 software.

We emphasize that the χ2 test will be used to identify whether there is an association between the variables and the OR to quantify this association and classify whether it characterizes risk or protection.

We also emphasize that the dichotomized dependent variable was having tuberculosis and residing in an identified cluster (0 and 1) and this variable was crossed with all the independent variables that were also dichotomized (0 and 1) and that category analyzed was considered as a reference. Only for variables with two categories (HIV, diabetes, alcoholism, mental disease, drug addiction and smoking) we adopted as a reference category not having the comorbidity or not using substances.

By observing the association between two independent dichotomous categorical variables, the Chi-square test can be used to generate a traditional p-value that verifies the existence or not of association. However, within the applied statistics, the chi-square p-value is of little value due to the loss of precision, accuracy and variance that comes with categorical variables and, in this way, the use of the Odds Ratio calculation is indicated with a confidence interval of 95% for a more accurate measurement.

Additionally, a binary logistic regression was also conducted because it is a more robust analysis. Initially, from the variables present in the TBWeb notification form, the variables that could explain the variable of interest were chosen. All selected variables were dichotomized (0 and 1) and the dependent variable (having TB and residing in the risk cluster) was also dichotomized.

Then, the selected variables were tested for multicollinearity based on the variance inflation factor (VIF), and those with an index greater than ten were excluded [12] (Additional file 1). After the selection process of the independent variables was completed, logistic regression was conducted using the RStudio 4.0.4 program.

The backward selection was used to insert the variables in the model. This method incorporates all the variables in the model and one by one it is removed and the individual contributions of the variables to this model are investigated, and the worst performing variable is eliminated, comparing the complete model with the reduced model, by removing such variable. After exhausting all the possibilities of the analysis, the best model was chosen based on the lowest values of the Akaike Information Criterion (AIC) [13]. It is also noteworthy that for the final model with the best comparison parameter, the Odds Ratio (OR) and their respective 95%CI were calculated.


Between 2006 and 2017, 2,259 tuberculosis cases were reported in Ribeirão Preto, of which 1760 (77.9%) were pulmonary tuberculosis. The minimum age of the cases was two months, and the maximum was 102 years old, with an average of 42 years old and a standard deviation of 17.3.

Figure 1 shows the behavior of the pulmonary tuberculosis time series and its time trend over the study period. It is possible to observe the presence of peaks and falls in specific periods of the years.

Fig. 1
figure 1

Series and time trend of pulmonary tuberculosis cases, Ribeirão Preto-SP, Brazil (2006–2017)

In the spatial analysis stage, 112 cases were excluded because they lived in a rural area, the address field on the notification form was not valid (wrong, blank, incomplete address or address of government agencies—hospitals, health units, prisons), making georeferencing impossible. Therefore, 1648 cases (93.6%) were georeferenced and integrated into the following analyses.

With the application of the purely spatial scanning technique, it was possible to identify four groups of statistically significant clusters (p < 0.01) (Fig. 2). Spatial cluster 1 (SC1), considered as protection for the event, presented RR: 0.19–0.52 (95% CI: 0.09–0.69), composed of 238 census sectors in the eastern, western, and central districts of the municipality. Spatial cluster 2 (SC2) with RR: 1.73 (95%CI: 1.49–1.99) composed 60 census sectors in the west district, two census sectors in the central district, and two census sectors in the north district.

Fig. 2
figure 2

Areas of spatial risk for the occurrence of pulmonary tuberculosis, Ribeirão Preto, SP, Brazil (2006–2017)

Spatial cluster 3 (SC3), presented RR: 2.07 (95%CI: 1.70–2.32), was composed of 52 census sectors in the north district and two census sectors in the central district, and finally, spatial cluster 4 (SC4) with RR: 2.68 to 2.72 (95%CI: 2.00–3.77), was composed of 24 census sectors in the southern district.

With the space–time scan, it was possible to identify four statistically significant clusters (p < 0.01), one being protective and three at risk for pulmonary tuberculosis (Fig. 3).

Fig. 3
figure 3

Space–time risk areas for the occurrence of pulmonary tuberculosis, Ribeirão Preto, SP, Brazil (2006–2017)

The space–time cluster 1 (STC1), considered protection for the event, presented RR: 0.13 (95%CI: 0.03–0.26) between 2008 and 2013, was composed of 87 census sectors in the eastern district. The space–time cluster 2 (STC2) presented RR: 1.94 (95%CI: 1.77–2.11) between 2010 and 2015; it was composed of 58 census sectors in the west, north, and central district.

The space–time cluster 3 (STC3), with RR: 2.34 (95%CI: 1.67–2.94) and the period from 2006 to 2011, was composed of 44 census sectors in the northern district and the space–time cluster 4 (STC4) (p < 0.01) presented RR: 2.84 (95%CI: 2.50–3.90) between the years 2014 and 2017, it was composed of 20 census sectors in the southern district.

With the application of SVTT, it was found that the cases of pulmonary tuberculosis in Ribeirão Preto showed an average decrease of 0.12% per year, making the findings presented in Fig. 1 in which periods of peaks and falls were more evident. Still, with this analysis, it was possible to identify a cluster with a statistically significant variation in the temporal trend (p < 0.01) (Fig. 4).

Fig. 4
figure 4

Areas with spatial variation in temporal trends for the occurrence of pulmonary tuberculosis, Ribeirão Preto, SP, Brazil (2006–2017)

The cluster with spatial variation in temporal trends 1 (SVTT1) has RR: 0.11 and was composed of 13 census sectors in the eastern district of the municipality. The cluster, classified as protection for the event, presented an ITT of growth (+ 0.09%/year), while the ETT indicated a decrease (–0.06%/year).

Table 1 shows the clinical-epidemiological profile of all pulmonary tuberculosis cases reported in Ribeirão Preto from 2006 to 2017 and the clinical-epidemiological profile of patients of pulmonary tuberculosis resident in the clusters identified with the techniques of the scan applied.

Table 1 Clinical and epidemiological profile of pulmonary tuberculosis cases, Ribeirão Preto, SP, Brazil (2006–2017)

From the analysis of Pearson's chi-square analysis (χ2) and OR, it was possible to identify two risk factors and three protection factors for the analyzed event. In the identified areas, mixed-race people (brown) are 1.77 times more likely and people without any studies are 1.16 times more likely to get sick from pulmonary tuberculosis than people living in other areas of the municipality.

Higher education (15 years or more of study) was identified as a protective factor for tuberculosis in the identified areas (OR: 0.73). Also classified as protection against the event are HIV-negative people who do not have diabetes (OR: 0.55 and 0.35, respectively).

From the logistic regression analysis, the variables that showed multicollinearity (VIF value > 10) were excluded (Supplementary material), as follows: age (15–59 years and 60 years or older) and sex, the other variables entered the logistic regression stage and the final model was composed of nine variables. Through OR, it was possible to identify three risk factors and two protective factors associated with the epidemiological situation in the municipality (Table 2). In the identified areas, people of race/color brown are 1.26 times more likely to fall ill with pulmonary tuberculosis and also people without any education (OR: 1.71).

Table 2 Factors associated with the epidemiological situation in Ribeirão Preto, SP, Brazil (2006–2017)

Higher education (8–11 years of study and 15 years or more of study) was identified as a protective factor against tuberculosis in the identified areas (OR: 0.72 and 0.26, respectively).


The study has evidenced the risk-prone territories for TB spread and the changes in the clusters over the years. When analyzing the time series of pulmonary tuberculosis cases between 2006 and 2017, it was possible to observe a series of undulations in its temporal trend, which was later classified through the SVTT, which indicated an average decrease of 0.12% per year of pulmonary tuberculosis in Ribeirão Preto, SP, Brazil.

With the analysis of purely spatial scanning, it was possible to verify the formation of clusters in areas that can be considered at risk for the occurrence and transmissibility of pulmonary tuberculosis.

SC4, classified with the highest RR (2.68–2.72) (95%CI: 2.00–3.77), was located in the southern district where the largest subnormal agglomerate (slums) of residents is found in the municipality, and it is noteworthy that Family Health Units have not yet been implemented in this area [14, 15]. SC3 (RR: 2.07) (95%CI: 1.70–2.32), located mostly in the northern district, is the region with the highest population density in the municipality and concentrates the largest number of residents per residence. It is worthy of note that there is a settlement of rural workers in this district, and this region also concentrates the most significant number of subnormal agglomerates (slums) in the city [15].

The SC2 (RR: 1.73) (95%CI: 1.49–1.99) identified in the western district of the municipality has one of the lowest municipal human development rates in the city, with the population receiving mostly two minimum wages. It is the region with the largest number of health units and has the highest percentage of exclusive users of the Unified Health System (SUS) in the city [15].

Finally, SC1 was classified as protection for the event (RR: 0.19–0.52) (95% CI: 0.09–0.69) and located in the eastern, western, and central districts. It is noteworthy that the center of Ribeirão Preto has the lowest density of inhabitants per residence and presents mainly commercial properties. In the eastern community, the population has the highest municipal income and the highest education level. It is noteworthy that this region has areas of expansion of high-standard condominiums. However, there are also population clusters here with lower socioeconomic resources that deserve attention because although this population is located in a more economically favored region in the municipality, they do not have the same resources and there is risk for other people who live in this area and thus can cause outbreak of tuberculosis or other infectious diseases [14, 15].

With the use of the space–time scanning statistics, four statistically significant clusters were identified; one for protection against illness from pulmonary tuberculosis (STC1; RR: 0.13) (95%CI: 0.03–0.26) between the years 2008 and 2013 in the eastern district, which may indicate that during this period there were improvements in surveillance services in health, such as awareness campaigns and active search for respiratory symptoms or probably there were fewer diagnoses of the disease in this period or no notification of new cases.

In contrast, STC2 (RR: 1.94) (95%CI: 1.77–2.11) in the period from 2010 to 2015, mainly in the western region, and STC3 (RR: 2.34) (95%CI: 1.67–2.94) in the period from 2006 to 2011 in the north district were classified as risk, which may indicate that there was intense activity of active search for respiratory symptoms, and therefore, there were more notifications than other regions of the municipality.

STC4 (RR: 2.84) (95%CI: 2.50–3.90) was identified in the southern district between the years 2014 to 2017. In 2014, the diagnosis was implemented in the municipality of Ribeirão Preto using the rapid molecular test for tuberculosis (RMT-TB) performed through GeneXpert®MTB/RIF. The system is commonly used for the detection of Mycobacterium tuberculosis and resistance to rifampicin (the primary drug used in the treatment) automatically [16, 17].

The municipality has two rapid diagnostic devices, one in the west and the other in the south, and studies [16,17,18,19] carried out to compare the RMT-TB with the tests commonly performed showed that the sensitivity of the RMT-TB for a sputum sample was 72.5% (for three samples it reached 90%) and the specificity reached 99%. Therefore, the hypothesis is raised that the high number of cases diagnosed in this region may be related to the implementation of RMT-TB in the municipality.

Through the use of the SVTT, a cluster with RR: 0.11 (event protection) was identified in the eastern region of the municipality, which, in the analysis of purely spatial scanning statistics, was identified as protection for the event (RR: 0.19–0.52). The SVTT1 cluster indicated a specific area in the east with growth ITT (0.09%/year), while outside that area, the identified time trend is decreasing (0.06%/year).

It is important that municipal managers focus on strategies to combat tuberculosis in this region with the identification of respiratory symptoms and their communicants. In addition, aiming at improving the epidemiological indicators of tuberculosis, a very necessary and effective tool is the periodic training of health professionals not only in the biological sense but also to understand the epidemiological and social reality of the territory in which they are inserted, so that they can understand the relationship between tuberculosis and the social determinants of health in order to fully look at the individual [20].

Geography and its sub-areas of knowledge, including geoprocessing, have great power to support decision-making based on the analysis of cases of a given disease, uniting their socio-spatial characteristics so that it is possible to highlight the disparities between population groups and also the delimitation of risk and/or spatial protection areas for a given event of interest.

Over the years with the evolution of technology, several studies have sought to understand the behavior of tuberculosis in space in Brazil [21,22,23,24], but few studies used temporal [25, 26] or Spatio-temporal [27, 28] approaches, and the numbers of studies that combined both techniques are even smaller [29, 30]. Therefore, we encourage new studies to be carried out using spatial and temporal processes.

In the scientific literature, there are a large number of studies that have analyzed the behavior of the disease in the territory, and that sought to identify areas of risk. However, few studies have addressed the temporal behavior of tuberculosis or sought to identify variations in the temporal trends in the risk territories.

With the association analysis, it was identified that, in the risk areas classified in the municipality, people of brown race/color were 1.26 (95%CI: 1.10–1.60) times more likely to get sick with pulmonary tuberculosis, corroborating the findings of the χ2 analysis, which indicated OR = 1.77 (95%CI: 1.01–1.86) for this population. Other studies have shown that people of mixed race/color are less likely to be cured [31] and that the mortality rate from tuberculosis is growing in this population [32].

No plausible biological justification or relationship was found in the literature to justify this difference between brown race/color and other classifications, but the construction of Brazilian Society [33] and the social determinants of health [34] must be taken into account.

It was also identified that people without any education had 1.71 (95%CI: 1.53–1.95) times more chance, and people with high schooling (8–11 years of study and 15 years or more of study) have less chance (OR: 0.72 (95%CI: 0.53–0.96) and 0.26 (95%CI: 0.07–0.91), respectively) of falling ill with pulmonary tuberculosis in high-risk areas in the municipality, corroborating the results of the χ2 analysis, which indicated OR = 1.16 (95%CI: 1.03–1.85) for people without schooling and OR = 0.73 (95%CI: 0.55–0.98) for people with 15 years of schooling or more. In a literature review[35] conducted with the aim of relating the level of education to infection with Mycobacterium tuberculosis, it was identified that the educational level of people with tuberculosis is directly related to income, highlighting the relationship between tuberculosis and social conditions of life.

Additionally, the χ2 analysis also identified not having HIV or diabetes as protection for the event (OR: 0.55 (95%CI: 0.08–0.63) and 0.35 (95%CI: 0.09–0.73), respectively). It is known that HIV is an aggressive virus that attacks the immune system, and, in this way, the chances of contracting opportunistic infections such as tuberculosis increase significantly. The impact of tuberculosis co-infection with HIV in the body is bidirectional; as the HIV virus grows in the body, the risk and progression of other infections, such as tuberculosis, also increases. In this way, tuberculosis delays the recovery of CD4 T-cells that are also destroyed by HIV, which increases the progression of the disease to AIDS and, consequently, the deaths from tuberculosis in PLHIV [36,37,38].

The same is true for people with TB and diabetes, so TB makes glycemic control difficult, and in turn, high blood glucose makes TB worse [39, 40]. It was identified in studies [39, 40] that people with diabetes had a 2.44–8.33 times greater chance of developing tuberculosis than those without the disease. This is because people with diabetes have decreased cellular and humoral immunity, in addition to hyperglycemia and cellular insulinopenia, which have effects on the function of macrophages and lymphocytes, thus increasing the chance of infections [41,42,43].

Concerning the limitations of this study, it was an ecological study; the so-called ecological fallacy stands out in that because variables are used at the aggregate level, the results may not represent associations at the individual level. It is also worth mentioning the use of secondary data sources, which may include incomplete data or typos.

We would also like to emphasize the importance and also encourage further studies using different methodologies to further explore tuberculosis and also its latent and extrapulmonary forms.


In view of the above, using statistical techniques and spatial analysis to understand the behavior of the disease over time and the location of areas in the municipality in which tuberculosis is a problem and using the characteristics of its surroundings to estimate the population at risk can assist managers in making assertive decisions, so that it becomes easier to understand the process about the chain of transmission of the disease and the entire context in which that population is inserted.

Availability of data and materials

The data that support the findings of this study are available from the Municipal Health Secretariat of Ribeirão Preto, but restrictions were applied to the availability of these data used under license for the current study, and so are not publicly available. Data are, however, available from the authors upon reasonable request and with permission of the Municipal Health Secretariat of Ribeirão Preto.



95% Confidence interval


External temporal trend


Human immunodeficiency virus


Institute of Geography and Statistics


Internal temporal trend



Km2 :

Square kilometers


Odds ratio


People living with the human immunodeficiency virus


Rapid molecular test for tuberculosis


Relative risk


Spatial cluster 1


Spatial cluster 2


Spatial cluster 3


Spatial cluster 4


Information System for Notifiable Diseases


São Paulo


Space–time cluster 1


Space–time cluster 2


Space–time cluster 3


Space–time cluster 4


Seasonal Trend Decomposition using Loess


Unified Health System


Spatial variation in temporal trends


Cluster with spatial variation in temporal trends 1


Co-infection with tuberculosis and human immunodeficiency virus


Tuberculosis Control system


World Health Organization


  1. World Health Organization (WHO). Global Tuberculosis Report 2021. Geneva: World Health Organization; 2019. Accessed 12 Mar 2022.

  2. Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Boletim Epidemiológico - Tuberculose 2021, 44p. Accessed 12 Mar 2022.

  3. Moraga P, Kulldorff M. Detection of spatial variations in temporal trends with a quadratic function. Stat Methods Med Res. 2016;25(4):1422–37.

    Article  Google Scholar 

  4. Kleinbaum DG, Kupper LL, Morgenstern H. Epidemiologic research: principles and quantitative Methods. 1st ed. USA: Wiley; 1982.

    Google Scholar 

  5. Instituto Brasileiro de Geografia e Estatística. Mapas: bases e referenciais. Rio de Janeiro: IBGE; 2010. Accessed 12 Mar 2022.

  6. Galesi VMN. Data on tuberculosis in the state of São Paulo, Brazil. Rev Saúde Pública. 2007;41(1):121.

    Article  Google Scholar 

  7. Cleveland R, Cleveland W, McRee JE. Seasonal-trend decomposition procedure based on LOESS. J Offic Stat. 1990;6:3–73.

    Google Scholar 

  8. Hyndman R, Athanasopoulos G, Bergmeir C, Caceres G, Chhay L, O'Hara-Wild M et al. Package ‘forecast’: forecasting functions for time series and linear models. Accessed 27 Oct 2020.

  9. Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Stat Med. 1995;14(8):799–810.

    CAS  Article  Google Scholar 

  10. Kulldorff M. SaTScan User Guide, 2018. Accessed 27 Oct 2020.

  11. Jaisankar R, Kesavan J. A study on spatial variations in temporal trends of dengue incidences In Tamil Nadu, India. Int J Sci Technol Res. 2019;8(09).

  12. Zuur AF, Ieno EM, Elphick CS. A protocol for data exploration to avoid common statistical problems. Methods Ecol Evol. 2010;1(1):3–14.

    Article  Google Scholar 

  13. Wagh YS, Kamalja KK. Zero-inflated models and estimation in zero-inflated Poisson distribution. Commun Stat Simul Comput. 2018;47(8):2248–65.

    Article  Google Scholar 

  14. Berra TZ, Bruce ATI, Alves YM, Campoy LT, Arroyo LH, Crispim JA, et al. Related factors, time trend and spatial association of abandonment of treatment for tuberculosis in Ribeirão Preto-SP. Rev Eletr Enferm. 2020;22:1–10.

    Article  Google Scholar 

  15. Ribeirão Preto. Secretaria Municipal de Saúde, Departamento de Vigilância em Saúde e Planejamento. Fatores de risco relacionados à saúde da população residente na zona urbana de Ribeirão Preto (SP) 2008–2011. Ribeirão Preto (SP). Accessed 27 Oct 2020.

  16. BRASIL. Ministério da Saúde. 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: Ministério da Saúde, 2015. Accessed 27 Oct 2020.

  17. Lima TM, Belotti NCU, Nardi SMT, Pedro HSP. GeneXpert MTB/RIF assay for diagnosis of tuberculosis. Rev Pan-Amaz Saude. 2017;8(2):67–78.

    Google Scholar 

  18. Shah W. To determine diagnostic accuracy of gene xpert and sputum Ziehl-Neelsen staining taking sputum culture as gold standard. Eur Respir J. 2016;48:PA2779.

    Article  Google Scholar 

  19. 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.

    Article  Google Scholar 

  20. Macedo SM, Andrade RPS, Souza CRBA, Andrade ASS, Villa TCS, Pinto ESG. Estratégias para capacitação ao cuidado em tuberculose. Cogitare Enferm. 2016;21(3):01–8.

    Article  Google Scholar 

  21. Prado-Junior J, Medronho RA. Spatial analysis of tuberculosis cure in primary care in Rio de Janeiro, Brazil. BMC Public Health. 2021;21(1):1841.

    Article  Google Scholar 

  22. Santos BA, Cruz RPS, Lima SVMA, Santos AD, Duque AM, Araújo KCGM, et al. Tuberculosis among children and adolescents: an epidemiological and spatial analysis in the state of Sergipe, Brazil, 2001–2017. Cien Saude Colet. 2020;25(8):2939–48.

    Article  Google Scholar 

  23. Mendes MS, Oliveira ALS, Pimentel LMLM, Fiqueiredo TMRM, Schindler HC. Spatial analysis of tuberculosis in children under 15 years of age and socioeconomic risk: an ecological study in Paraíba, Brazil, 2007–2016. Epidemiol Serv Saude. 2021;30(3): e20201038.

    Article  Google Scholar 

  24. Andrade HLP, Ramos ACV, Crispim JA, Santos Neto M, Arroyo LH, Arcêncio RA. Spatial analysis of risk areas for the development of tuberculosis and treatment outcomes. Rev Bras Enferm. 2021;74(2): e20200564.

    Article  Google Scholar 

  25. Souza CDF, Paiva JPS, Silva LF, Leal TC, Magalhães MAFM. Trends in tuberculosis mortality in Brazil (1990–2015): joinpoint analysis. J Bras Pneumol. 2019;45(2): e20180393.

    Article  Google Scholar 

  26. Berra TZ, Gomes D, Ramos ACV, Alves YM, Bruce ATI, Arroyo LH, et al. Effectiveness and trend forecasting of tuberculosis diagnosis after the introduction of GeneXpert in a city in south-eastern Brazil. PLoS ONE. 2021;16(5): e0252375.

    CAS  Article  Google Scholar 

  27. Roza DL, Caccia-Bava MCGG, Martinez EZ. Spatio-temporal patterns of tuberculosis incidence in Ribeirão Preto, State of São Paulo, southeast Brazil, and their relationship with social vulnerability: a Bayesian analysis. Rev Soc Bras Med Trop. 2012;45(5):607–15.

    Article  Google Scholar 

  28. Santos JA, Santos DT, Arcêncio RA, Nunes C. Space–time clustering and temporal trends of hospitalizations due to pulmonary tuberculosis: potential strategy for assessing health care policies. Eur J Public Health. 2021;31(1):57–62.

    Article  Google Scholar 

  29. Giacomet CL, Santos MS, Berra TZ, Alves YM, Alves LS, Costa FBP, et al. Temporal trend of tuberculosis incidence and its spatial distribution in Macapá-Amapá. Rev Saude Publica. 2021;55:96.

    Article  Google Scholar 

  30. Arcêncio RA, Berra TZ, Terena NFM, Rocha MP, Alecrim TFA, Kihara FMS, et al. Spatial clustering and temporal trend analysis of international migrants diagnosed with tuberculosis in Brazil. PLoS ONE. 2021;16(6): e0252712.

    Article  Google Scholar 

  31. Santos JN, Sales CMM, Maciel EL. Factors associated with cure when treating tuberculosis in the state of Rio de Janeiro, Brazil, 2011–2014. Epidemiol Serv Saude. 2018;27(3): e2017464.

    PubMed  Google Scholar 

  32. Arcoverde MAM, Berra TZ, Alves LS, Santos DT, Belchior AS, Ramos ACV, et al. How do social-economic differences in urban areas affect tuberculosis mortality in a city in the tri-border region of Brazil, Paraguay and Argentina. BMC Public Health. 2018;18:795.

    Article  Google Scholar 

  33. Maio MC. Raça, doença e saúde pública no Brasil: um debate sobre o pensamento higienista do século XIX. In: Maio MC, Santos RV, editors. Raça como questão: história, ciência e identidades no Brasil. Rio de Janeiro: Fiocruz; 2010. p. 51–82.

    Chapter  Google Scholar 

  34. World Health Organization (WHO). A conceptual framework for action on the social determinants of health. Geneva: World Health Organization; 2010. Accessed 12 Mar 2022.

  35. Rodrigues MW, Mello AGNC. Tuberculose e escolaridade: Uma revisão da literatura. Rev Int Apoyo Inclusión, Logopedia, Sociedad y Multiculturalidad. 2018;4(2).

  36. Bell LCK, Noursadeghi M. Pathogenesis of HIV-1 and Mycobacterium tuberculosis co-infection. Nat Rev Microbiol. 2018;16(2):80–90.

    CAS  Article  Google Scholar 

  37. Meintjes G, Brust JCM, Nuttall J, Maartens G. Management of active tuberculosis in adults with HIV. Lancet HIV. 2019;6(7):e463–74.

    Article  Google Scholar 

  38. Tornheim JA, Dooley KE. Tuberculosis associated with HIV infection. Microbiol Spectr. 2017;5(1)

  39. Silva DR, Muñoz-Torrico M, Duarte R, Galvão T, Bonini EH, Arbex FF, et al. Fatores de risco para tuberculose: diabetes, tabagismo, álcool e uso de outras drogas. J Bras Pneumol. 2018;44(2):145–52.

    Article  Google Scholar 

  40. Pereira SM, Araújo GS, Santos CAST, Oliveira MG, Barreto ML. Associação entre diabetes e tuberculose: estudo caso controle. Rev Saúde Pública. 2016;50:82.

    PubMed  PubMed Central  Google Scholar 

  41. Marais BJ, Lönnroth K, Lawn SD, Migliori GB, Mwaba P, Glaziou P, et al. Tuberculosis comorbidity with communicable and non-communicable diseases: integrating health services and control efforts. Lancet Infect Dis. 2013;13(5):436–48.

    Article  Google Scholar 

  42. Ayelign B, Negash M, Genetu M, Wondmagegn T, Shibabaw T. Immunological impacts of diabetes on the susceptibility of Mycobacterium tuberculosis. J Immunol Res. 2019;2019:6196532.

    Article  Google Scholar 

  43. Restrepo BI. Diabetes and tuberculosis. Microbiol Spectr. 2016;4(6):10.

    Article  Google Scholar 

Download references


Municipal Health Secretariat of Ribeirão Preto. Municipal Tuberculosis Control Program of Ribeirão Preto


This work was supported by the São Paulo State Research Support Foundation (FAPESP) (FAPESP Process no. 2018/03700-7).

Author information

Authors and Affiliations



TZB, ACVR, and RAA participated in the conception, planning, analysis, interpretation, and writing of the work; LHA, FMD, JAC, YMA, FBPC, MSS, and LSA participated in the report of the work; Author RCF, AAM, and DG participated on interpretation and writing of the work. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Thaís Zamboni Berra.

Ethics declarations

Ethics approval and consent to participate

In compliance with Resolution 499/2012 of the National Health Council, the study was approved by the Research Ethics Committee of Nursing College of Ribeirão Preto, University of São Paulo, under the Certificate of Ethical Appraisal number 87696318.3.0000.5393 issued on April 30, 2019. Consent to participate was not applicable because we worked with secondary data of cases diagnosed with tuberculosis and reported on TBWeb (tuberculosis case notification system used in the state of São Paulo). We declare that all applied methods were carried out following all guidelines and precepts for their validity.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Test for multicollinearity based on the variance inflation factor (VIF).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Berra, T.Z., Ramos, A.C.V., Arroyo, L.H. et al. Risk-prone territories for spreading tuberculosis, temporal trends and their determinants in a high burden city from São Paulo State, Brazil. BMC Infect Dis 22, 515 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Tuberculosis
  • Spatial analysis
  • Temporal trend