Skip to main content

Finding gaps in TB notifications: spatial analysis of geographical patterns of TB notifications, associations with TB program efforts and social determinants of TB risk in Bangladesh, Nepal and Pakistan



In order to effectively combat Tuberculosis, resources to diagnose and treat TB should be allocated effectively to the areas and population that need them. Although a wealth of subnational data on TB is routinely collected to support local planning, it is often underutilized. Therefore, this study uses spatial analytical techniques and profiling to understand and identify factors underlying spatial variation in TB case notification rates (CNR) in Bangladesh, Nepal and Pakistan for better TB program planning.


Spatial analytical techniques and profiling was used to identify subnational patterns of TB CNRs at the district level in Bangladesh (N = 64, 2015), Nepal (N = 75, 2014) and Pakistan (N = 142, 2015). A multivariable linear regression analysis was performed to assess the association between subnational CNR and demographic and health indicators associated with TB burden and indicators of TB programme efforts. To correct for spatial dependencies of the observations, the residuals of the multivariable models were tested for unexplained spatial autocorrelation. Spatial autocorrelation among the residuals was adjusted for by fitting a simultaneous autoregressive model (SAR).


Spatial clustering of TB CNRs was observed in all three countries. In Bangladesh, TB CNR were found significantly associated with testing rate (0.06%, p < 0.001), test positivity rate (14.44%, p < 0.001), proportion of bacteriologically confirmed cases (− 1.33%, p < 0.001) and population density (4.5*10–3%, p < 0.01). In Nepal, TB CNR were associated with population sex ratio (1.54%, p < 0.01), facility density (− 0.19%, p < 0.05) and treatment success rate (− 3.68%, p < 0.001). Finally, TB CNR in Pakistan were found significantly associated with testing rate (0.08%, p < 0.001), positivity rate (4.29, p < 0.001), proportion of bacteriologically confirmed cases (− 1.45, p < 0.001), vaccination coverage (1.17%, p < 0.001) and facility density (20.41%, p < 0.001).


Subnational TB CNRs are more likely reflective of TB programme efforts and access to healthcare than TB burden. TB CNRs are better used for monitoring and evaluation of TB control efforts than the TB epidemic. Using spatial analytical techniques and profiling can help identify areas where TB is underreported. Applying these techniques routinely in the surveillance facilitates the use of TB CNRs in program planning.

Peer Review reports


Tuberculosis (TB) is an infectious respiratory disease which affects millions of people all over the world, but overwhelmingly affects the most vulnerable, hard to reach and socioeconomically disadvantaged people [1] . Despite efforts to drastically reduce the burden of TB, it remains within the top ten causes of deaths worldwide and outranks HIV/AIDS as one of the leading causes of death from an infectious disease [2]. In 2017, an estimated 10.0 million people fell ill from TB. However, only 6.4 million were notified to the national authorities. This means that an estimated 3.6 million people with TB were either not detected by the health system or not notified to the local authorities, and therefore missed by the formal health systems [2]. Many of these missing people with TB do not receive the care they need, leaving them vulnerable to develop severe and potentially fatal infections as well as being a potential source of transmission to those around them [1,2,3].

There are several barriers that TB patients may encounter which causes them to not be diagnosed or notified by the health systems. These barriers are well described by Uplekar et al. and others [4,5,6,7,8,9,10,11]. First of all, the patient needs to recognize the symptoms, but due to misperceptions, lack of knowledge or even social or internalized stigma, the patient might fail to do so [4,5,6]. Secondly, the patient needs to seek healthcare which can be compromised by distance to health facilities or transportation costs, loss of wages, costs for diagnosis and treatment and the perception of poor quality services [4, 7, 8]. Thirdly, health workers might fail to recognize the symptoms due to lack of training or lack of human resources [4, 9]. Fourth, the TB patient must be diagnosed as such, which can be complicated by insensitive screening and diagnostic tests, delay between testing and diagnosis that leads to loss to follow-up, or the inability for the individual with suspected TB to produce sputum [4, 10]. Fifth, the patient with TB might not initiate treatment due to direct and indirect treatment costs, geographical access to TB services, inadequate knowledge of the importance of timely treatment and stigma [11]{Citation}. Finally, due to poor knowledge of reporting procedures of the health care provider or poor engagement with the private health care sector, TB patients started on treatment might not be notified to the authorities [11].

TB case notification rates (CNR) are geographically heterogeneous [12,13,14,15]. Geographic variations in the presence of risk populations, TB transmission and burden are potential drivers of geographic variations in TB CNR. In Portugal for example, clusters of TB were associated with higher HIV/AIDS incidence, household crowding and incarceration [16]. Whereas in Brazil, higher rates of TB were associated with poor economic conditions, non-white population, urbanization and household crowding [17].

However, TB CNR are also a function of TB control efforts (i.e. the extent to which the health system effectively reaches out, diagnoses and notifies people with TB) and are therefore not necessarily an appropriate indicator of TB incidence. TB CNR can be only used as an indicator of TB incidence in places with a strong control program. A strong control program should take into account the subnational variation in TB burden and tailor its control efforts to the local epidemic [18].

It is our hypothesis that TB CNRs are only reflective of TB incidence in places where the TB control efforts are tailored to the local epidemic. Nonetheless, TB CNRs are often used as a proxy for TB incidence in the absence of a more reliable indicator of TB burden. If TB CNRs are used for the allocation of resources, we need a better understanding of what is driving the case notification rates. Therefore, the objective of this study is to gain better understanding of the drivers of subnational variations in TB notification rates in three South Asian countries with intermediate to high TB burden; Bangladesh, Nepal and Pakistan. The aim is to use spatial analytical techniques and profiling to gain better understanding of sociodemographic, health and programmatic indicators that underlie spatial variation of TB case notification rates.



This study uses aggregated subnational data from Bangladesh, Nepal and Pakistan. These countries were purposely selected based on the availability of recent subnational data on administrative areas (i.e. spatial data), TB notification, TB program efforts and indicators of demography and health; as well as their geographical proximity and thus similarities in geographic context. The study focusses on TB notification data from 2014 and 2015.


Data on TB case notification and programmatic indicators routinely collected by the National Tuberculosis Programs (NTPs) of Bangladesh (2015), Nepal (2014) and Pakistan (2015) were obtained for the most recent and complete year. Additional data on demography and health were derived from publicly available reports, such as the Demographic and Health Survey (DHS), Multiple Indicator and Cluster Survey (MICS), population census reports and statistical yearbooks [19,20,21,22,23,24,25,26,27,28,29,30,31] . These reports were collected online from governmental websites such as the USAID and the Bureau of Statistics of each country [32,33,34,35] . Only reports with subnational data were included and reports prior to 2010 were excluded from this study. One dataset was made per country, which included all available data on district level or higher, leading to three datasets: one dataset covering 64 administrative units of Bangladesh, one dataset covering 75 administrative units of Nepal and one dataset covering 142 administrative units of Pakistan.

To be able to map and visualise TB indicators, spatial data files of the respective administrative boundaries for each of the three countries were obtained from online spatial data repositories [36,37,38]. As subnational boundaries and administrative units change throughout the years, adjustments to the spatial data files were made (i.e. merging of 2 administrative units into 1) to ensure that the spatial data is the same as the TB reporting units.


Subnational, annual TB case notification rates were used as the primary outcome variable for all analyses. It is defined as the number of reported cases per 100,000 population. It includes all reported cases, independent of the type of TB, the diagnostic method that had been used, or the patient’s TB history [39].

The covariates that were used for this study were selected foremost because of their known association to TB burden or TB programme performance, but also on the availability of these indicators in public reports. The covariates that were used for this study can be grouped into four broader themes. An overview of all covariates and their definitions can be found in Table 1.

  1. (1)

    Access to healthcare; the TB case notification rate tends to be higher in areas with better access to healthcare [40,41,42].

  2. (2)

    Socioeconomic status (SES); TB is associated with the socially and economically disadvantaged [17, 43].

  3. (3)

    Demography and key populations; Certain demographic and socio-economic key populations are known to have an increased risk of developing clinical TB, these include: miners, migrants, males, the malnourished and the elderly [44,45,46,47,48].

  4. (4)

    Quality of TB treatment and diagnostic services; The ability of the healthcare system to detect and treat TB patients is associated to various performance indicators [48, 49].

Table 1 Defintions of covariates

In order to link the covariate data to the units in the spatial data files, every geographical unit was assigned an unique identifier in the form of a numerical code. The same code was given to the spatial information belonging to that geographical unit. This way, the compiled subnational data were merged with the spatial information and subsequently visualized in GIS software.

Statistical analysis

Spatial autocorrelation

The global univariable Moran’s I was used to test for the presence of spatial clustering (spatial autocorrelation) in the TB case notification rate at subnational level. In addition, Local Indicators of Spatial Autocorrelation (LISA) were calculated in order to identify and locate clusters of districts with a relatively high or low TB case notification rate. A first order queen contiguity connectivity matrix was used to define neighbouring districts. These statistics give valuable information about the independency (i.e. absence of spatial dependency) of observations, an important assumption in regression analyses [50].

Univariable and multivariable analyses

Generalized linear models (GLM) with log transformed outcome variable (TB CNR) were fitted to the data for the univariable and multivariable analyses - separately for each country - based on the fit and nature of the data, as disease rates often follow a log-normal distribution [51]. A poisson and negative binomial model were also considered as these are frequently used to model disease rates, but the lognormal model provided a better fit for all three countries. First, univariable analyses were conducted for each covariate to assess the strength and direction of the association. Next, all variables were included in a multivariable GLM Finally, simultaneous autoregressive (SAR) models were fitted to the data to account for unexplained spatial autocorrelation in the residuals of the multivariable GLM. This model uses a spatially correlated error structure based on a contiguity weights matrix taking into account only directly adjacent spatial neighbours (i.e. shared borders). The reported coefficients are exponentiated and reflect percentage change in TB CNR. All models are reported with corresponding 95% confidence interval (CI) and p- values, as well as the global Moran’s I statistic with corresponding p-value for the residuals of the multivariable models.

Pearson’s correlation coefficients were calculated for all covariates before multivariable modelling. If two variables were correlated (i.e. correlation coefficient exceeding 0.7), one of the variables was excluded from further analysis. The variable to be excluded was determined by the number of correlations, where the variable with highest number of correlations was excluded. In case both variables had equal numbers of correlations than the variable with a non-significant result in the univariable analysis was excluded.

Data analyses were performed using GeoDa version for the assessment of spatial autocorrelation, RStudio version 1.0.143 for all other statistical analyses and QGIS version 2.18.4 for the geographical visualization of the data [52,53,54].


Descriptive analysis

A total of 205,98; 37,025 and 326,152 people were diagnosed with TB and notified in Bangladesh (2015), Nepal (2014) and Pakistan (2015) respectively. Subnational TB case notifications per hundred thousand population ranged from 50 to 187 in Bangladesh, from 31 to 227 in Nepal and from 9 to 689 in Pakistan (Table 2, 3 and 4).

Table 2 Descriptive statistics for Bangladeshi districts
Table 3 Descriptive statistics for Nepali districts
Table 4 Descriptive statistics for Pakistani districts

Spatial clustering of TB case notification rates was observed in Nepal (Moran’s I: 0.52, p-value < 0.001, Fig. 1). In the northern mountainous region of Nepal clusters of low notification rates (cold spots) were observed whereas in the lower “Terai” region in the south of Nepal several clusters of high notification rates (hot spots) were identified (Fig. 2) Moderate spatial clustering of TB case notification rates were observed in Bangladesh (Moran’s I: 0.23, p-value < 0.01, Fig. 3), with one larger hot spot in the eastern division Sylhet and two small cold spots in Rajshahi and Dhaka divisions (Fig. 4).

Fig. 1
figure 1

Moran’s I plot showing the strength of spatial autocorrelation of TB CNR (2014), Nepal

Fig. 2
figure 2

Map of Nepal showing the geographical variations and spatial clustering of the TB CNR 2014

Fig. 3
figure 3

Moran’s I plot showing the strength of spatial autocorrelation of TB CNR (2015), Bangladesh

Fig. 4
figure 4

Map of Bangladesh showing the geographical variations and spatial clustering of the TB CNR 2015

Moderate spatial clustering was also observed for Pakistan (Moran’s I: 0.23, p-value < 0.001, Fig. 5), with a large cold spot covering almost the entire province of Balochistan and one major hot spot in Punjab province (Fig. 6).

Fig. 5
figure 5

Moran’s I plot showing the strength of spatial autocorrelation of TB CNR (2015), Pakistan

Fig. 6
figure 6

Map of Pakistan showing the geographical variations and spatial clustering of the TB CNR 2015

Statistical results

Tables 5, 6 and 7 show the results of the univariable, non-spatial multivariable and spatial multivariable models of the natural logarithm of the TB case notification rate for Bangladesh, Nepal and Pakistan.

Table 5 Results of the univariable and multivariable analyses for Bangladesh and global Moran’s I of the residuals of the multivariable model
Table 6 Results of the univariable and multivariable analyses for Nepal and global Moran’s I of the residuals of the multivariable models
Table 7 Results of the univariable and multivariable analyses for Pakistan and global Moran’s I of the residuals of the multivariable models

Based on the correlation matrices (Additional file 1) and the univariable analyses, the following covariates were excluded prior to multivariable modelling: facility density and migrant population (Bangladesh), miners and under 5 mortality rate (Nepal), and literacy rate (Pakistan).


The spatial GLM model for Bangladesh showed that TB CNR is positively associated to population density (β = 4.5e-3; 95%CI: 1.7e-3 - 7.3e-3), testing rate (β = 0.06, CI: 0.06–0.07) and positivity rate (β = 14.44, CI: 12.44–16.48). Furthermore, the model suggests that an increase in the proportion of bacteriologically confirmed cases is associated with a decrease in TB case notification rate of 1.33% (CI: − 1.65 - -1.01). The model fully accounted for spatial autocorrelation in TB CNR (Moran’s I: − 0.07, p-value = 0.28), therefore no SAR model was fitted.


The residuals of the non-spatial multivariable model for Nepal were found to spatially auto correlate (Moran’s I: 0.17, p-value < 0.01) and a spatial model correcting for the special dependencies was computed. The spatial multivariable model shows that facility density is inversely associated to the TB CNR, where an increase of one unit in facility density is associated with a decrease of the TB CNR of 0.19% (CI: − 0.37 - -0.01). Furthermore, an increase of one unit in sex ratio is associated with an increase in TB CNR of 1.54% (CI: 0.40–2.69). Treatment success rate (β:-3.68, CI:-5.27 - -2.08) and stunting (β:-1.88, CI: − 3.15 - -0.60), are both inversely associated to TB CNR in the spatial model. The spatial model accounted for all spatial autocorrelation (Moran’s I: − 0.02, p-value = 0.44).


According to the spatial multivariable model TB CNR (which was necessary with a Moran’s I: 0.11 after the GLM) was found positively associated with testing rate (β = 0.08, CI: 0.05–0.10), test positivity rate (β = 4.29, CI: 2.73–5.87), proportion bacteriologically diagnosed (β = − 1.45, CI: − 2.02 – − 0.88), vaccination coverage (β = 1.17, CI: 0.61–1.72) and facility density (β = 20.41, CI: 8.42–33.72). The SAR model fully accounted for the spatial autocorrelation of the TB CNR (Moran’s I: − 0.04, p-value = 0.14).


From TB prevalence studies we know that there is a gap between TB notification and TB burden. Onozaki et al. assessed national TB prevalence surveys in Asia from 1990 to 2012 and found the TB prevalence to be twice as high as the number of notified cases [55]. Yet TB notification rates are often used as a measure of TB incidence, making it imperative that we have a better understanding of the drivers of TB notification rates.

The data and results presented in this paper show that TB case notifications across the three countries analysed are spatially heterogeneous and spatially clustered. The positive association between population density in Bangladesh and sex ratio in Nepal suggest that part of the variation in TB CNR can be explained by proxies for TB risk. In Bangladesh an increase in population density is positively associated with TB CNR (B:4.5e-3, 95%CI: 1.7e-3 – 7.3e-3), a higher population density is indicative of more crowding which increases transmission of TB [56]. Crowding has been associated with increased risk of TB in Bangladesh in both adults and children [57, 58]. Furthermore, a positive association between sex ratio and TB CNR in Nepal (B:1.54, 95%CI:0.40–2.69) shows that more TB patients are diagnosed in districts with more men. This is line with findings of the TB prevalence survey, which found a much higher prevalence in men [59].

However, the models also suggest that part of the variation in TB CNR can be explained by programmatic factors. TB CNR is inversely associated with the proportion of TB patients with a bacteriologically confirmation in both Bangladesh (B:-1.33, 95%CI:-.65 - -1.01) and Pakistan (B:-1.45, 95%CI: − 2.02 - -0.88). One explanation is that more sensitive diagnostic methods can result in a reduction of clinically diagnosed pulmonary TB, which may cause a reduction in the overall TB CNR [60].

The positive association between testing rate and TB CNR in Bangladesh (B:0.06, 95%CI:0.06–0.07) and Pakistan (B:0.08, 95%CI:0.05–0.10) suggests that more testing yields higher notification. In a well-adjusted system, the level of testing is direct response to the local TB burden. But the positive association between test positivity rates in both countries (Bangladesh: B:14.44, 95%CI:12.44–16.48, Pakistan B:4.29, 95%CI:2.73–5.87) suggests otherwise. When the testing rate increases, one would expect the positivity rate to decrease or remain stable. The increasing test positivity rate suggests that the current level of testing does not meet the local need for testing or that testing efforts may be targeted towards populations more likely to suffer from TB.

The positive association between TB CNR with facility density (B:20.41, 95%CI:8.42–33.72) and vaccination coverage (B:1.17, 95%CI:0.61–1.72) in Pakistan substantiate the results above. Facility density and vaccination coverage are both indicative of better access to health care and improved uptake of health services [61].

In contrast with the findings above, facility density is inversely associated with TB CNR in Nepal (B:-0.19, 95%CI:-0.37 - -0.01). One explanation is that in areas with very low population numbers (such as the hilly and mountainous regions in Nepal), the facility density is likely to be very high. Likewise, in highly urbanized areas (e.g. Kathmandu) the facility to population ratio can be very low due to the high population denominator. Although the placement of health facilities takes into account population density and the unmet need of the population, the facility to population ratio calculated on district level does not reflect this and might therefore not be the right metric to assess access to healthcare on a district level.

The inverse association between treatment success rate and TB CNR in Nepal (B:-3.68, 95%CI:-5.27 - -2.08) suggests that higher notification rates negatively affect case management, possibly due a higher burden on the health system to follow-up on TB patients or to provide the required medication.

The TB CNR decreases with 1.88% (95%CI: − 3.15 - -0.60) with increasing prevalence of stunting (Nepal). Stunting is strongly associated with lower socioeconomic status and poor health, which increases the risk of TB, but in this case may reflect lower access to healthcare. This is in line with findings from the TB prevalence survey, where the proportion of TB patients not seeking care was higher among the poor [59].

Although the models largely agree with one another, we see differences in the associations between sociodemographic and access indicators with the CNR. This finding underscore that the results are to some extent influenced by country context, specifically health care infrastructure and TB epidemiology.

Strengths and limitations

We were able to combine data from multiple sources, hence capturing different dimensions of the same phenomenon as well as minimizing inadequacies that occur in data from a single source. Furthermore, the models are congruent in what they suggest and the strength and direction of the associations are also consistent between the models. These similarities between independent models decrease the likelihood of the observed associations being the result of chance. In addition, most of the spatial clustering is adjusted for in the final models.

The data that were used in this study were derived from various sources and were therefore available on different administrative levels and for different years. However, we expect demographic indicators such as poverty, migration and population size to remain stable over the course of three to 5 years or to grow proportionally over time. Although natural disasters and conflict are known to influence demographic indicators, these changes are often not reflected in district level population statistics as presented by national statistical offices which also cover a broader time frame. However, the possibility remains that data from different years does not accurately reflect the situation in the year from which TB data were available. Furthermore, some data were only available on a higher administrative level (e.g. province or regional). These data could not be used to reflect the district-level variations that this study is trying to address. In addition, these data points are not independent from one another, which increases the risk of a type II error. Finally, subnational data on HIV or TB prevalence were not available for this study and could therefore not be included.


The results give clear indications of spatial clustering of the tuberculosis case notification rates in Bangladesh, Nepal and Pakistan. Where this is a result that is often found by other research concerning spatial epidemiology of TB, most of these studies attribute this to social and demographic indicators and neglect the influence that TB program efforts might have.

The result of this study show that notification of TB is mainly associated with access to healthcare and TB program efforts. This is not necessarily a problem in case the local efforts are a direct response to the TB burden, in which TB notification rates can be used as a proxy for TB incidence. However, the associations that were found do not suggest that the TB program efforts are a response to TB burden. In fact, they suggest that TB notification is dependent on programmatic response such as the ability to test, diagnose and treat individuals, but also the ability of patients to access health care. Hence, TB notifications should not be used as a proxy for TB incidence.

However, TB notifications are a great source of information if they are interpreted in the context of the local health system. As such, assessing changes over time in the geographical distribution of TB notification rates can be useful to monitor changes in policy, interventions or programmatic efforts. Spatial analytical techniques and profiling allows for the identification of spatial outliers and local inconsistencies which can be indicative of TB under notification. This valuable information can be used to prioritize areas which require further supervision and tailor interventions according to their local needs in an effort to diagnose and successfully treat the missing people with TB.

Availability of data and materials

Part of the data used for this study comes from publicly available documents such as the Demographic and Health Survey or reports from National Statistical Offices. All publicly available documents used for this study can be found in the references of this article.

All TB data that support the finding of this study are available from the National TB programs from which they were retrieved but restrictions apply to the availability of these data, which were 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 National TB programs.



Confidence interval


Case notification rate


Demographic and Health Survey


Global Moran’s I


Geographical Information System


Generalized Linear Model


Human Immunodeficiency Virus


Local Indicators of Spatial Autocorrelation


Multiple Indicator Cluster Survey


National TB Program


Simultaneous Autoregressive model


Socio economic status




  1. Stop TB Partnership. The paradigm shift 2016–2020: Global Plan to End TB [Internet]. Geneva: Stop TB Partnership; 2015 [cited 2018 Jan 10]. Available from:

  2. World Health Organization. Global tuberculosis report 2018 2018.

    Google Scholar 

  3. Knechel NA. Tuberculosis: pathophysiology, clinical features, and diagnosis. Crit Care Nurse. 2009 Apr 1;29(2):34–43.

    Article  Google Scholar 

  4. Uplekar M, Creswell J, Ottmani S-E, Weil D, Sahu S, Lönnroth K. Programmatic approaches to screening for active tuberculosis. Int J Tuberc Lung Dis Off J Int Union Tuberc Lung Dis. 2013 Oct;17(10):1248–56.

    Article  CAS  Google Scholar 

  5. Hadley M, Maher D. Community involvement in tuberculosis control: lessons from other health care programmes. Int J Tuberc Lung Dis. 2000;4(5):401–8.

    CAS  PubMed  Google Scholar 

  6. Somma D, Thomas BE, Karim F, Kemp J, Aria N, Auer C, et al. Gender and socio-cultural determinants of TB-related stigma in Bangladesh, India, Malawi and Colombia [special section on gender and TB]. Int J Tuberc Lung Dis. 2008;12(7):856–66.

    CAS  PubMed  Google Scholar 

  7. Abebe G, Deribew A, Apers L, Woldemichael K, Shiffa J, Tesfaye M, et al. Knowledge, health seeking behavior and perceived stigma towards tuberculosis among tuberculosis suspects in a rural community in Southwest Ethiopia. PLoS One. 2010;5(10):e13339.

    Article  Google Scholar 

  8. Needham DM, Bowman D, Foster SD, Godfrey-Faussett P. Patient care seeking barriers and tuberculosis programme reform: a qualitative study. Health Policy. 2004;67(1):93–106.

    Article  Google Scholar 

  9. Floyd K, Blanc L, Raviglione M, Lee J-W. Resources required for global tuberculosis control. Science. 2002;295(5562):2040–1.

    Article  CAS  Google Scholar 

  10. Joseph HA, Shrestha-Kuwahara R, Lowry D, Lambert LA, Panlilio AL, Raucher BG, et al. Factors influencing health care workers’ adherence to work site tuberculosis screening and treatment policies. Am J Infect Control. 2004;32(8):456–61.

    Article  Google Scholar 

  11. de Vries SG, Cremers AL, Heuvelings CC, Greve PF, Visser BJ, Bélard S, et al. Barriers and facilitators to the uptake of tuberculosis diagnostic and treatment services by hard-to-reach populations in countries of low and medium tuberculosis incidence: a systematic review of qualitative literature. Lancet Infect Dis [Internet]. [cited 2017 Apr 4]; Available from:

  12. Nanoo A, Izu A, Ismail NA, Ihekweazu C, Abubakar I, Mametja D, et al. Nationwide and regional incidence of microbiologically confirmed pulmonary tuberculosis in South Africa, 2004–12: a time series analysis. Lancet Infect Dis. 2015;15(9):1066–76.

    Article  Google Scholar 

  13. Tiwari N, Adhikari CMS, Tewari A, Kandpal V. Investigation of geo-spatial hotspots for the occurrence of tuberculosis in Almora district, India, using GIS and spatial scan statistic. Int J Health Geogr. 2006;5(1):33.

    Article  Google Scholar 

  14. Kakchapati S, Yotthanoo S, Choonpradup C. Modeling tuberculosis incidence in Nepal. Asian Biomed. 2010;4:2.

    Article  Google Scholar 

  15. Wang T, Xue F, Chen Y, Ma Y, Liu Y. The spatial epidemiology of tuberculosis in Linyi City, China, 2005–2010. BMC Public Health. 2012;12(1):885.

    Article  Google Scholar 

  16. Couceiro L, Santana P, Nunes C. Pulmonary tuberculosis and risk factors in Portugal: a spatial analysis. Int J Tuberc Lung Dis. 2011;15(11):1445–55.

    Article  CAS  Google Scholar 

  17. Maciel ELN, Pan W, Dietze R, Peres RL, Vinhas SA, Ribeiro FK, et al. Spatial patterns of pulmonary tuberculosis incidence and their relationship to socio-economic status in Vitoria, Brazil. Int J Tuberc Lung Dis. 2010;14(11):1395–402.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Theron G, Jenkins HE, Cobelens F, Abubakar I, Khan AJ, Cohen T, et al. Data for action: collection and use of local data to end tuberculosis. Lancet. 2015;386(10010):2324–33.

    Article  Google Scholar 

  19. NIPORT/Bangladesh NI of PR and T-, Associates M and, International ICF. Bangladesh Demographic and Health Survey 2014. 2016 [cited 2018 Jan 11]; Available from:

  20. Bangladesh EPI Coverage Evaluation Survey 2015. Mohakhali, Bangladesh: Center for Social and Market Research; 2016.

  21. Department of Health Services. Annual Report 2013/2014. Kathmandu: Ministry of Health and Population; 2015.

    Google Scholar 

  22. Elderly Population in bangladesh: Current Features and Guture Perspectives. Dhaka, Bangladesh: Bangladesh Bureau of Statistics; 2015. (Population Monograph). Report No.: 4 2015.

    Google Scholar 

  23. FATA Development Indicators household survey 2013–14 [internet]. Peshawar, Pakistan: Bureau of Statistics; 2015 [cited 2017 Jul 20]. Available from:

  24. Naveed A, Wood G, Ghaus MU. Geography of Poverty in Pakistan– 2008-09 to 2012–13: Distribution, Trends and Explanations. Think Asia [Internet]. 2016 Jun 30 [cited 2018 Jan 11]; Available from:

  25. Literacy Trends in Pakistan. Islamabad, Pakistan: UNESCO Office;.

  26. Nepal Demographic and Health Survey 2011 [Internet]. Kathmandu, Nepal: Ministry of Health and Population, New ERA, and ICF International, Calverton, Maryland: Ministry of Health and Population (MOHP) [Nepal], New ERA, and ICF International Inc; 2012 [cited 2017 Jul 20]. Available from:

  27. Sharma P, Guha-Khasnobis B, Khanal DR. Nepal Human Development Report 2014. Kathmandu: Government of Nepal National Planning Commission; 2014.

    Google Scholar 

  28. Nepal Living Standards Survey 2010/2011 [Internet]. Kathmandu, Nepal: Central bureau of Statistics; 2011 [cited 2017 Jul 20]. Available from:

  29. Rahman MN, Kamal GM, Rahman MM, Sultana S, Abdullah HP, Parvin A, Mondal SR, Alam MD. Population and Housing Census 2011 - Socio-Economic and Demographic Report. Dhaka: Bangladesh Bureau of Statistics; 2012. (National Series). Report No.: 4.

    Google Scholar 

  30. Progotir Pathey Multiple Indicator Cluster Survey 2012-2013, Key Findings [Internet]. Dhaka, Bangladesh: Bangladesh Bureau of Statistics and UNICEF Bangladesh; 2014 [cited 2017 Jul 20]. Available from:

  31. Social Demography. Katmandu, Nepal: Central Bureau of Statistics; 2014. (Population Monograph of Nepal). Report No.: 2.

  32. Bangladesh Bureau of Statistics-Government of the People’s Republic of Bangladesh | বাংলাদেশ পরিসংখ্যান ব্যুরো-গণপ্রজাতন্ত্রী বাংলাদেশ সরকার [Internet]. [cited 2018 Jan 10]. Available from:

  33. Central Bureau of Statistics [Internet]. [cited 2018 Jan 10]. Available from:

  34. Pakistan Bureau of Statistics [Internet]. [cited 2018 Jan 10]. Available from:

  35. The DHS Program - Quality information to plan, monitor and improve population, health, and nutrition programs [Internet]. [cited 2018 Jan 10]. Available from:

  36. Global Administrative Areas | Boundaries without limits [Internet]. [cited 2018 Jan 10]. Available from:

  37. DIVA-GIS | free, simple & effective [Internet]. [cited 2018 Jan 10]. Available from:

  38. Team G. GeoNetwork opensource portal to spatial data and information [Internet]. 2007 [cited 2018 Jan 11]. Available from:

  39. Organization WH. Definitions and reporting framework for tuberculosis–2013 revision. 2013;.

    Google Scholar 

  40. Levesque J-F, Harris MF, Russell G. Patient-centred access to health care: conceptualising access at the interface of health systems and populations. Int J Equity Health. 2013 Mar 11;12:18.

    Article  Google Scholar 

  41. Dangisso MH, Datiko DG, Lindtjørn B. Accessibility to tuberculosis control services and tuberculosis programme performance in southern Ethiopia. Glob Health Action. 2015;8 Available from:

  42. Datar A, Mukherji A, Sood N. Health infrastructure & immunization coverage in rural India. Indian J Med Res. 2007;125(1):31–42.

    PubMed  Google Scholar 

  43. Lönnroth K, Jaramillo E, Williams BG, Dye C, Raviglione M. Drivers of tuberculosis epidemics: the role of risk factors and social determinants. Soc Sci Med. 2009 Jun 1;68(12):2240–6.

    Article  Google Scholar 

  44. Horton KC, MacPherson P, Houben RMGJ, White RG, Corbett EL. Sex differences in tuberculosis burden and notifications in low- and middle-income countries: a systematic review and meta-analysis. PLoS Med. 2016 Sep 6;13(9):e1002119.

    Article  Google Scholar 

  45. WHO | Moderate malnutrition [Internet]. WHO. [cited 2017 Jun 16]. Available from:

  46. Dhavan P, Dias HM, Creswell J, Weil D. An overview of tuberculosis and migration. Int J Tuberc Lung Dis Off J Int Union Tuberc Lung Dis. 2017;21(6):610–23.

    Article  CAS  Google Scholar 

  47. Ross MH. Occupational respiratory disease in mining. Occup Med. 2004 Aug 1;54(5):304–10.

    Article  CAS  Google Scholar 

  48. World Development Indicators 2007. World Bank Publications; 2007. 430 p.

  49. McCarthy KD, Metchock B, Kanphukiew A, Monkongdee P, Sinthuwattanawibool C, Tasaneeyapan T, et al. Monitoring the performance of mycobacteriology laboratories: a proposal for standardized indicators. Int J Tuberc Lung Dis. 2008 Sep 1;12(9):1015–20.

    CAS  PubMed  Google Scholar 

  50. Twisk JWR. Inleiding in de toegepaste biostatistiek; 2014.

    Google Scholar 

  51. Limpert E, Stahel WA, Abbt M. Log-normal distributions across the sciences: keys and clues: on the charms of statistics, and how mechanical models resembling gambling machines offer a link to a handy way to characterize log-normal distributions, which can provide deeper insight into variability and probability—normal or log-normal: that is the question. AIBS Bull. 2001;51(5):341–52.

    Google Scholar 

  52. QGIS D. QGIS geographic information system. Open source geospatial Foundation project. 2015.

  53. Team Rs. RStudio: integrated development for R. Boston MA URL Httpwww Rstudio Com: RStudio Inc; 2015.

    Google Scholar 

  54. Anselin L, Syabri I, Kho Y. GeoDa: an introduction to spatial data analysis. Geogr Anal. 2006;38(1):5–22.

    Article  Google Scholar 

  55. Onozaki I, Law I, Sismanidis C, Zignol M, Glaziou P, Floyd K. National tuberculosis prevalence surveys in Asia, 1990-2012: an overview of results and lessons learned. Tropical Med Int Health. 2015 Sep;20(9):1128–45.

    Article  Google Scholar 

  56. Myers WP, Westenhouse JL, Flood J, Riley LW. An ecological study of tuberculosis transmission in California. Am J Public Health. 2006 Apr;96(4):685–90.

    Article  Google Scholar 

  57. Irfan SD, Faruque MO, Islam MU, Sanjoy SS, Afrin D, Hossain A. Socio-demographic determinants of adult tuberculosis: a matched case-control study in Bangladesh. Am J Infect Dis. 2017 Mar 1;13(3):32–7.

    Article  Google Scholar 

  58. Karim M, Rahman M, Mamun S, Alam M, Akhter S. Risk factors of childhood tuberculosis: a case control study from rural Bangladesh. WHO South-East Asia J Public Health. 2012;1(1):76.

    Article  Google Scholar 

  59. NATIONAL TB PREVALENCE SURVEY (2018–19) - FACTSHEET. National Tuberculosis Control Centre;.

  60. Creswell J, Rai B, Wali R, Sudrungrot S, Adhikari LM, Pant R, et al. Introducing new tuberculosis diagnostics: the impact of Xpert® MTB/RIF testing on case notifications in Nepal [Internet]. 2015 [cited 2020 Apr 3]. Available from:

  61. Lakew Y, Bekele A, Biadgilign S. Factors influencing full immunization coverage among 12–23 months of age children in Ethiopia: evidence from the national demographic and health survey in 2011. BMC Public Health. 2015 Jul 30;15(1):728.

    Article  Google Scholar 

Download references


Not applicable.


No funding was obtained for this study

Author information

Authors and Affiliations



MG analysed the data. RF, PJ, SCV and AHK participated in interpretation and contextualization of the results. ER, MIB, LB and CM conceived of the study and participated in its design and coordination. All authors have read an approved the final manuscript.

Corresponding author

Correspondence to Margo van Gurp.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

No competing interest to declare.

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.

Correlation Matrices of Bangladesh, Nepal and Pakistan. Three individual correlation matrices for Bangladesh, Nepal and Pakistan.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

van Gurp, M., Rood, E., Fatima, R. et al. Finding gaps in TB notifications: spatial analysis of geographical patterns of TB notifications, associations with TB program efforts and social determinants of TB risk in Bangladesh, Nepal and Pakistan. BMC Infect Dis 20, 490 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: