Skip to main content
  • Research article
  • Open access
  • Published:

Disparities in access to diagnosis and care in Blantyre, Malawi, identified through enhanced tuberculosis surveillance and spatial analysis

Abstract

Background

A sizeable fraction of tuberculosis (TB) cases go undiagnosed. By analysing data from enhanced demographic, microbiological and geospatial surveillance of TB registrations, we aimed to identify modifiable predictors of inequitable access to diagnosis and care.

Methods

Governmental community health workers (CHW) enumerated all households in 315 catchment areas during October–December 2015. From January 2015, government TB Officers routinely implemented enhanced TB surveillance at all public and private TB treatment registration centres within Blantyre (18 clinics in total). This included collection from registering TB patients of demographic and clinical characteristics, a single sputum sample for TB microscopy and culture, and geolocation of place of residence using an electronic satellite map application. We estimated catchment area annual TB case notification rates (CNRs), stratified by microbiological status. To identify population and area-level factors predictive of CHW catchment area TB case notification rates, we constructed Bayesian spatially autocorrelated regression models with Poisson response distributions. Worldpop data were used to estimate poverty.

Results

In total, the 315 CHW catchment areas comprised 753,489 people (range 162 to 13,066 people/catchment area). Between 2015 and 2017, 6077 TB cases (61% male; 99% HIV tested; 67% HIV positive; 55% culture confirmed) were geolocated, with 3723 (61%) resident within a CHW catchment area. In adjusted models, greater distance to the nearest TB registration clinic was negatively correlated with TB CNRs, which halved for every 3.2-fold (95% CI 2.24–5.21) increase in distance. Poverty, which increased with distance from clinics, was negatively correlated with TB CNRs; a 23% increase (95% CI 17–34%) in the mean percentage of the population living on less than US$2 per day corresponded to a halving of the TB case notification rates.

Conclusions

Using enhanced surveillance of TB cases in Blantyre, we show an ecological relationship consistent with an ‘inverse care law’ whereby poorer neighbourhoods and those furthest from TB clinics have lower relative CNRs. If confirmed as low case detection, then pro-poor strategies to facilitate equitable access to TB diagnosis and treatment are required.

Peer Review reports

Background

A strong socio-economic gradient between poverty and incidence of tuberculosis (TB) has been demonstrated at both population and individual levels [1, 2], with the poorest in society having the greatest risk of developing active TB disease [3]. Poor housing quality, malnutrition, HIV and unhealthy behaviours such as smoking appear to be the major intermediary determinants [4]. People living in poverty also face the greatest difficulties in accessing TB diagnosis and treatment [5] and frequently experience catastrophic care-seeking costs [6]. The World Health Organization (WHO) End-TB strategy has called for ambitious policies to ensure all TB cases are detected and treated and to ensure that no person affected by TB faces socio-economic ruin [7]. However, in sub-Saharan Africa, a substantial ‘detection gap’ exists, with an estimated 51% of prevalent TB cases remaining undiagnosed in the community in 2016 [1].

There is likely substantial heterogeneity in the determinants of incident active TB and in access to TB diagnosis care within cities [8], meaning that preventive interventions and improved access to high-quality TB services targeted towards local neighbourhoods with the greatest need could plausibly result in improvements in TB control [9]. However, identifying these ‘hot-spots’ of inequitable access to TB care depends on integrating high-quality epidemiological and microbiological surveillance with geospatial technologies. Lower case notification rates from certain areas may reflect either worse access to care or lower TB burden. Where lower notification rates coincide with higher levels of known risk factors for TB (e.g. poverty), they may instead reflect barriers to access rather than a relative absence of TB.

We report on the implementation and evaluation of an enhanced TB surveillance programme in Blantyre, Malawi, where nearly one-in-five adults are HIV positive [10]. We aimed to use high-quality spatially resolved surveillance data to support public health practitioners and policymakers to identify the most important and modifiable barriers to TB diagnosis and treatment. We hypothesised that neighbourhood characteristics such as poverty rates, age and sex structure, and distance from TB registration clinics might be correlated with TB case notification rates, and so guide the implementation of interventions to improve access to TB diagnosis and prevention.

Methods

Study site and population

Blantyre is a major urban centre located in the south of Malawi. Blantyre is divided into non-overlapping geographical areas, each served by locally resident community health workers (CHW) known as Health Surveillance Assistants. CWHs are employed by the Ministry of Health of Malawi and serve as a link between health facilities and communities.

Enumeration of Health Surveillance Assistant areas

To estimate population denominators for the study, we defined boundaries around 315 CHW catchment areas, excluding business and industrial areas and the most affluent areas of the city. Through community engagement activities, and with the support of the Blantyre District Health Officer, we provided training in census enumeration to CHWs in Blantyre. Between 10 October 2015 and 30 December 2015, CHWs, accompanied by study Research Assistants, undertook a circumferential walk around their catchment area to record sets of boundary coordinates using global positioning satellite devices. Where there was overlap between two contiguous CHW catchment areas, the circumferential walks were repeated and boundaries revised to ensure that no areas overlapped.

CHWs and Research Assistants subsequently undertook a door-to-door census from all dwellings within each CHW catchment area. We defined a dwelling to be a physical structure within which members of one or more households had slept in the previous week. We excluded dwellings that were derelict or abandoned. We defined a household to be a person living alone or a group of people living together who shared meals together and who may have been related or unrelated. The head of each household (or, in their absence, the next most senior household member) was interviewed using a structured census questionnaire that recorded the number of households and people residing within the dwelling. Five per cent of enumerated households were randomly selected, revisited and re-enumerated by a separate Research Assistant team for quality control purposes.

Enhanced surveillance of tuberculosis case notifications

Within Blantyre’s administrative border, there are a total of 18 health facilities where TB Officers employed by the Ministry of Health register TB cases and initiate treatment for TB. As TB diagnostic and treatment centres are not evenly spread through Blantyre, geographical proximity to diagnostic services may be an important predictor of accessing care.

At all health facilities in the city, patients are investigated for TB by health workers, usually with an initial symptom screen and subsequent examination of a diagnostic sample, most commonly sputum smear microscopy. Once a decision to initiate TB treatment has been made by a health worker—either on the basis of a positive microbiological sample or on clinical decision—the patient is referred to their nearest TB registration clinic, which will usually be located within the same health facility. Here, Ministry of Health TB Officers initiate TB treatment and register the patient using the National TB Register, recording demographic and clinical characteristics. In this study, we implemented enhanced TB surveillance at the point of TB registration.

We identified all TB Officers who registered TB cases at the 18 health facilities in Blantyre (one public and two private hospitals, 13 public and one private primary health care centres, and one prison clinic), and provided study-specific training in enhanced TB surveillance activities. In addition to registering TB cases using the National TB Programme register, TB Officers were asked to complete an enhanced surveillance form that recorded registration centre, age, sex, any sputum smear results taken as part of routine care (positive, negative or not done), TB treatment category (using standard Ministry of Health/WHO categories), TB classification (pulmonary or extrapulmonary TB) and HIV status (positive, negative or unknown). Monthly, study coordinators held joint meetings with all TB Officers and research team members where clinic TB Registers were reconciled with enhanced surveillance data.

Geolocation of TB cases

Following registration of a TB case, TB Officers used a satellite map application (electronic patient locator [ePAL]) that we had previously developed and evaluated to geolocate the physical location of each case’s dwelling. ePAL and its precursor—the Mapbook—have been described before [11, 12]. In brief, we purchased high-resolution satellite maps from the European Space Agency, and CHWs geotagged important points of reference within their catchment areas, including schools, clinics, bars, roads, churches and other physical landmarks. Throughout the study period, points of reference were periodically added-to and updated. Within the ePAL application map screen, TB patients searched for points of interest known to be near their dwelling using a text matching algorithm, then scrolled and zoomed to identify their dwelling; a subsequent long-press on the map screen recorded a set of coordinates from the selected dwelling location.

Microbiological surveillance

TB Officers asked all registering adult TB cases (regardless of TB category or classification) to submit a single sputum sample for smear and culture. Sputum samples were collected daily from each registration centre and transported in cooler boxes to the TB Research Laboratory at the College of Medicine, University of Malawi, for fluorescence microscopy and MGIT culture. Patients with positive results requiring additional clinical input were traced and linked to care.

Statistical analysis

We took the total population denominators of each CHW catchment area from the study census enumeration. CHW catchment area population densities were calculated and plotted in people per square kilometre. Using the ‘st_distance’ function in the R ‘sf’ package [13], we estimated the shortest direct Cartesian distance between the centroid of each CHW catchment area and the nearest TB registration clinic, without accounting for interceding physical structures. We overlaid 1-km gridded raster images containing data on the percentage of the population living on less than US $2 per day obtained from the Worldpop Project [14] onto CHW catchment area boundaries; where CHW catchment areas were covered by more than one 1-km grid square, we took the mean of the percentages.

We used the ‘st_within’ function in the R ‘sf’ package [13] to classify the CHW catchment area of residence of each case. Where notified TB cases were resident within a study-mapped CWH catchment area, we classified them as CHW catchment area residents; where TB cases were resident in an area of the city not mapped by study activities, or in another part of the country, we classified them as non-CHW catchment area resident. We compared the characteristics of TB cases that resided within and outside of study CHW catchment areas using proportions and means. Using census denominators, we estimated CHW catchment area TB case notification rates (CNRs) per 100,000 per year and disaggregated by TB microbiological status (smear positive or culture positive for Mycobacterium tuberculosis on study enhanced surveillance sample vs. all TB cases).

Using the CHW catchment area as the unit of analysis, we undertook two Bayesian conditional autocorrelated regression modelling analyses with Poisson response distributions [15]: the first to estimate all TB CNRs and the second to estimate TB CNRs that were microbiologically confirmed by the study laboratory as part of enhanced surveillance activities. We included within both models linear predictors for the log10(Poisson rate), including log10 population density, mean number of people per household, log10 Cartesian distance in metres from the nearest TB registration centre, percentage of population living on less than US $2 per day, adult male-to-female ratio, the percentage of the population that was aged 15 years or older and—as a marker of late presentation for diagnosis—the ratio of sputum smear-positive to smear-negative notified cases (on smears done by the routine health system; where it was not possible to estimate a ratio, either because no cases were notified from that CHW catchment area or there were no smear-positive cases, we inputted a value of 1). We additionally included an offset term for log10 population size of each CHW catchment area so as to model per capita rates. We included a diffuse Gaussian prior with mean 0 and standard deviation 10 on each covariate and the population-averaged intercept.

The spatial conditional autoregressive structure [16] used to account for spatial correlations in TB CNRs was a Besag-York-Mollié prior using a nearest-neighbours matrix method where the two catchment areas with the shortest Cartesian distance between CHW centroids were classified as neighbours. To determine the optimal nearest-neighbours correlation matrix, we plotted network graphs of CHW catchment areas, varying from k = 0 (no spatial autocorrelation prior) to k = 6 (each CHW catchment area was correlated with five other CHW catchment areas). To select the final models, we constructed for each analysis seven Bayesian regression models including priors k = 0 to k = 6, compared their widely applicable information criteria (WAIC) statistics [17] and selected the model with the lowest WAIC value.

Samples from posterior distributions were drawn using Markov chain Monte Carlo (MCMC) methods, specifically the non-U-turn sampler (NUTS) compiled in Stan using the ‘brms’ version 2.4.0 package [18] running in R version 3.5.1 (R Foundation for Statistical Computing, Vienna). For both analyses (all TB CNRs and microbiologically confirmed TB CNRs), we ran three chains of 4000 iterations with 1000 warm-up samples. We checked model convergence by estimating Gelman-Rubin statistics, inspecting trace plots, plotting residuals and plotting the geographical distribution of model-fitted to observed case notifications. We estimated posterior means and 95% credible intervals and plotted sample posterior distributions to examine the strength of evidence for covariate correlation with TB CNRs per CHW catchment area.

Results

Population characteristics

A total of 753,489 people resident within 315 CHW catchment areas were enumerated (Fig. 1). CHW populations ranged from 162 people to 13,066 people, with a mean of 2392 (standard deviation 1504) and median of 2043 (25th percentile 1338; 75th percentile 3035). The mean catchment area male-to-female adult ratio was 1.04 and ranged from 0.83 to 1.35 (Additional file 1: Figure S1). Overall, 61.5% of catchment area residents were aged 15 years or older, with cluster level percentages of adults ranging from 44.4 to 78.0%. The mean percentage of the population living in poverty (< $US2 per day) was 32% (range 13 to 79%).

Fig. 1
figure 1

Health Surveillance Assistant boundaries and populations 2015, Blantyre, Malawi. Population counts from study census conducted from October to December 2015. Boundaries are community health worker catchment areas, recorded by circumferential walks recording global positioning satellite coordinates by community health worker study team during October to December 2015. Black triangles are TB registration clinics. White areas in the centre of the map are mountainous areas or business districts with few residents that were not enumerated in the census or mapped by CHWs. The black triangle in the far north of map is a health centre with a TB registration clinic located outside of Blantyre District that may be used by Blantyre residents; to increase accuracy of case notification rates, we captured TB registrations at this clinic

TB case notifications

Between January 2015 and December 2017, 7799 TB cases were registered for TB treatment at the 18 TB registration clinics. Overall, 6077 (78%) registering TB cases were captured by the ePAL geolocation system. Of these, 3723 cases were resident within a CHW catchment area and 2354 were not: 1722 residing in a District outside Blantyre city, and 632 resided inside Blantyre city but outside of a study CHW catchment area or in an area of the city that had not been mapped by study activities and where enhanced surveillance GPS coordinates could not be obtained.

Characteristics of CHW-resident and non-CHW-resident TB cases were broadly similar, but with some important differences (Table 1). Overall, 61% (3727/6077) of TB cases were male, 60% (3629/6076) had pulmonary TB and 60% (1536/4089) were recorded as being sputum smear positive in the clinic TB register. Ascertainment of HIV status by TB Officers was very high, with HIV status recorded in > 99% of cases. Overall, 67% (3964/5915) of TB cases were HIV positive. However, compared to non-CHW area-resident cases, CWH area-resident cases were more likely to be male (64% vs. 58%), have pulmonary TB (62% vs. 55%), be sputum smear positive on routine clinic sample (62% vs. 56%) and on study laboratory sample (45% vs. 39%) and be culture positive on study laboratory sample (59% vs. 49%).

Table 1 Baseline characteristics of TB cases recorded in enhanced surveillance system in Blantyre

Sputum samples taken for TB culture as part of enhanced TB surveillance activities were successfully collected and reported for 71% (4306/6077) of cases overall and 73% (2170/3723) of CHW catchment area-resident cases. Overall, 59% of TB cultures done among CHW catchment area-resident cases were positive, with all speciated as Mycobacterium tuberculosis. A further 103 cases were sputum smear positive, but culture negative, meaning that for this analysis, 54% (1700/3723) of CHW catchment area-resident cases were classified as microbiologically confirmed.

TB case notification rates

Annual CHW catchment area TB CNRs ranged from 0 to 986 per 100,000 population per year. Across all 3 years, CNRs were highest in the densely populated areas of central urban Blantyre (Fig. 2). Catchment areas in the poorer and less densely populated outer suburbs recorded very few TB case notifications.

Fig. 2
figure 2

Annual Health Surveillance Assistant TB case notification rates, Blantyre, Malawi: 2015–2017. a All TB cases. b Microbiologically confirmed TB cases. Case notification rates per 100,000 population per year. Boundaries are community health worker catchment areas. Black triangles are TB registration clinics. White areas in the centre of the map are mountainous areas or business districts with few residents that were not enumerated in the census. The black triangle in the far north of the map is a health centre with a TB registration clinic located outside of Blantyre District that may be used by Blantyre residents; to increase accuracy of case notification rates, we captured TB registrations at this clinic

Factors associated with TB case notifications

We ran seven Bayesian regression analyses to determine the optimal spatial autocorrelation prior for inclusion in each final regression analysis (all TB CNRs and microbiologically confirmed TB CNRs)—Additional file 2: Figure S2 and Additional file 3: Figure S3. For analysis 1 (all TB cases) and analysis 2 (microbiologically confirmed TB cases), a k = 6 and k = 4 nearest-neighbour correlation prior respectively provided the lowest WAIC statistic, and there was evidence for better fit when compared to models without spatial autocorrelation priors (Additional file 4: Figure S4).

Both final regression analyses (analysis 1: all TB case notification rates; analysis 2: microbiologically confirmed TB case notification rates) showed evidence of convergence with stable mixing across chains. Gelman-Rubin statistics were 1.00 for all covariates, and residual plots showed evidence of good model fit (Additional file 4: Figure S4). Posterior probability distributions for fully adjusted covariate effects are shown in Table 2, and marginal effect plots are shown in Fig. 3.

Table 2 Posterior distributions of population-level covariates in a Bayesian spatial model for predicting tuberculosis case notification rates in Blantyre, Malawi: 2015–2017
Fig. 3
figure 3

Marginal effects of covariates on TB case notification rates, Blantyre, 2015–2017. Red: analysis 1 (all TB); blue: analysis 2 (microbiologically confirmed TB). Estimated by fitting a Bayesian spatial regression model with Poisson response, a k nearest-neighbours conditional spatial autocorrelation prior (with k = 6 for analysis 1 and k = 4 for analysis 2), with linear terms fitted for community health worker catchment area log10 population density, adult M:F ratio, mean number of people per household, log10 Cartesian distance from geographical centroid to the nearest TB clinic, percentage of population aged 15 years or older, mean percentage living on less than US $2 per day, offset term for log10 population size, and with weakly informative prior on the population-level effects intercept (Gaussian: mean = 0, sd = 10), and predictor intercept (Gaussian, mean = 0, sd = 10). Lines are medians, estimated from 100 draws from posterior samples. Shaded areas are 95% credible intervals. Each marginal distribution estimated by holding all other variables constant at their median

Adjusting for the effects of other covariates, we found strong evidence for a negative correlation between CHW catchment area TB CNRs and greater Cartesian distance between the catchment area centroid and the nearest TB clinic (analysis 1: relative rate [RR] per log10 increase in clinic distance 0.60, 95% credible interval [CI] 0.42–0.86); analysis 2: RR 0.55, 95% CI 0.36–0.84). For analyses 1 and 2, the TB case notification rate halved for every 3.2-fold (95% CI 2.24–5.21) and 3.5-fold (95% CI 2.28–6.86) increase in distance respectively.

We additionally found strong evidence for a negative correlation between TB CNRs and the percentage of CHW catchment area population that were living in poverty (analysis 1: RR 0.97, 95% CI 0.96–0.98; analysis 2: RR 0.97, 95% CI 0.95–0.98). In adjusted analysis, for analysis 1, a 23% (95% CI 17–34%) increase in the mean percentage of the population living on less than US$2 per day corresponded to a halving of the TB CNRs. The effect was similar for analysis 2 (23%, 95% CI 14–34%).

For both analyses, there was evidence in support of a correlation between higher TB CNRs and lower CHW catchment area male-to-female adult ratios (i.e. where a greater percentage of the adult population was female). For analysis 1, a 1.8% (95% CI 1.03–7.07%) decrease in the male-to-female adult ratio was associated with a 10% reduction in the CNR for all TB cases. For analysis 2 (microbiologically confirmed cases), each 1.6% (95% CI 0.91–7.72%) decrease in the male-to-female adult ratio was associated with a 10% reduction in the CNR.

The sputum smear positive to negative ratio (from samples collected for TB diagnosis by the routine health system) showed mixed results. For analysis 1 (all TB cases), the TB CNR halved for every 2.9% (95% CI 1.61–22.76%) decrease in the smear positive to negative ratio. In contrast, for analysis 2 (microbiologically confirmed TB), there was little evidence of an association between TB CNRs and the smear positive to negative ratio (RR 1.19, 95% CI 0.94–1.47).

There was little evidence to support correlations between TB CNRs and population density, mean number of people per household or the proportion of the population aged 15 years or older.

Discussion

Through enhanced TB surveillance in Blantyre, Malawi, we found evidence supporting the ‘inverse case law’, with adjusted TB case notification rate ratios for the poorest neighbourhoods, and those furthest from TB registration clinics, being substantially lower than for less poor neighbourhoods and for those in closer proximity to TB diagnosis and care centres. Evidence was consistently strong in models investigating microbiologically confirmed TB disease and all TB cases, supporting our hypothesis that poor access to TB diagnosis and treatment is predominately determined by poverty, and disproportionately affects the most disadvantaged in society. If a high burden of undiagnosed TB is confirmed through future prevalence surveys, then policymakers should strongly consider prioritising the implementation of pro-poor interventions to improve access to TB services in these underserved neighbourhoods. This could include expanding or relocating primary health care centres with comprehensive TB and HIV services [19], periodic mobile TB screening camps [20] or targeted community-wide active case-finding interventions [21, 22]. There may also be a potential impact from wider development interventions, not specific to TB.

Cities—and in particular informal settlements—provide conditions highly conducive to the spread of infectious disease [23], although it is clear that there is considerable social and geographical heterogeneity in the prevalence and incidence of infections, including TB [2, 9]. The prevalence of undiagnosed active TB is known to be high in many urban Africa cities [21, 22]. However, understanding of how and where to best implement TB control strategies in resource-limited settings has been hampered by a lack of census data and TB case notification data disaggregated by age, sex and HIV status and at lower-level geospatial boundaries.

In Blantyre, we successfully implemented enhanced TB surveillance in an extremely densely populated setting, including routine sputum sample collection for mycobacterial culture and geolocation using an innovative electronic application. We believe this approach could be widely applied in other similar settings and would be beneficial to public health authorities, providing neighbourhood-level data on underserved areas to inform TB care and prevention strategies. Critical to achieving the high levels of case-capture were the partnerships we established with key stakeholders, including the Malawi National TB Programme and Blantyre District Health Office, and with the Blantyre TB Officers and community health workers who undertook most of the enhanced surveillance activities.

We additionally found strong evidence that neighbourhoods with lower male-to-female adult ratios (i.e. a higher proportion of women) had higher TB CNRs even after adjusting for other important determinants of access to care. Although we cannot speculate about the causal nature of the relationship between neighbourhood sex ratios and TB case notification rates, if confirmed by prevalence survey data, this suggests that lower case notification rates may, in part, be driven by male barrier to TB diagnosis. Community TB prevalence and case notifications are known to be substantially higher for men than women [24], perhaps driven by gendered patterns of social mixing and barriers in care-seeking [25]. Areas in Blantyre with the highest male-to-female ratios are centred around extremely densely populated trading areas, where young men migrate for work. For such individuals, intense economic pressure may mean that seeking care for symptoms of TB is not a priority, a hypothesis that is supported by our previous qualitative work with symptomatic men in the city [25]. Conversely, residents of neighbourhoods where sex ratios were equal, or where women outnumber men, may be relatively more stable, with social stability a prerequisite for prioritising individual or household health.

As a further potential marker of poor access and late presentation for care that could be investigated using routine programmatic surveillance data, we investigated associations between sputum smear positive to negative notification ratios (on samples collected by the routine TB programme) and neighbourhood TB case notification rates. Smear positivity is a marker of infectivity that increases with diagnostic delay but, once present, is considerably easier for routine programmes to diagnose than smear-negative TB; thus, greater diagnostic access barriers are likely to be reflected by a higher proportion of smear-positive patients. Here, we show that TB case notification rates halved with every 2.4-fold increase in the sputum smear positive to negative ratio. This supports our hypothesis that late and/or under-diagnosis of TB from the poorest, more distant neighbourhoods was explained by poor access to diagnosis rather than by less prevalent disease. However, confirming this would need neighbourhood-level data from TB prevalence surveys (and in particular prevalence-to-notification ratios), ideally combined with local data on other intermediary determinants.

There were a number of strengths and limitations to this study. Overall, we successfully recorded clinical and demographic data from 78% (6077/7799) of all registering TB cases at all 18 health facilities in urban Blantyre. We additionally completed microbiological surveillance by sputum culture from 71% of captured cases, and geolocated 70% of all TB cases registered at city health centres. We concentrated mapping activities on the largest and most deprived neighbourhoods of the city and excluded business districts and the most affluent neighbourhoods. Thus, we will have systematic under-sampled TB patients from middle-class homes. Reflecting these different underlying populations, non-CHW catchment area-resident cases did have different characteristics compared to CHW-resident cases. Additionally, because we could only collect geospatial data from TB cases who were resident within neighbourhoods mapped by our study activities, we are unable to include TB cases resident outwith mapped areas in our regression analysis. There therefore may be important predictors of access to TB diagnosis and care among populations excluded from the present analysis that we were unable to identify. Surveillance data were collected by routine Ministry of Health employees (TB Officers) using the electronic GPS application and not by research staff. Being programmatic data, we anticipate some imprecision in geolocating TB patients. The household enumeration data was limited to basic data on age and sex, and so we are unable to estimate correlations with important intermediary determinants—notably HIV prevalence; the relationship with sputum smear positivity rates at the population level remains uncertain without detailed HIV prevalence data. We relied on data from the Worldpop Project from for estimating neighbourhood-level poverty. TB cultures were not done on 27% of CHW-resident TB patients who did not submit sputum, including 17% of pulmonary TB cases: this is likely to include a bias against less gravely ill patients too sick to produce sputum, or started on TB treatment out-of-hours while inpatients. Overall, 58% of notified TB cases initiated on TB treatment by the National TB Programme in study health facilities were microbiologically confirmed by the study laboratory, reflecting the programmatic and diagnostic challenges in confirming TB disease in low-resource settings. Finally, without supporting epidemiological data, our hypothesis relating to under-diagnosis in the more remote neighbourhoods of Blantyre remains speculative.

Future priorities include investigation of the prevalence of undiagnosed TB in neighbourhoods with low CNRs but high poverty and poor access to care indicators, where TB may be underdiagnosed: if so, then the prevalence of undiagnosed infectious TB may in fact be very high in these areas. Community-level TB interventions including active case finding have tended to focus on high-density informal residential areas, guided by high CNRs and prevalence survey data that show urban settings to be high risk for TB [26]. Even within densely populated urban settings such as Blantyre, there may be substantial heterogeneity in TB prevalence and risk factors for transmission, meaning that the potential effects of community-wide interventions may become diluted. We note that across 3 years studied, there were consistent patterns of high TB case notifications arising from the same set of CHW catchment areas. It is possible that case-finding interventions directed towards areas with the greatest need (both in terms of risk factors for TB and in likely under-diagnosis) could have greater and quicker impact on controlling TB than broader population-based strategies. However, further epidemiological and modelling studies, as well as trials, are required.

Conclusions

In summary, we found evidence from statistical analysis of enhanced city-wide surveillance data that the poorest neighbourhoods of Blantyre, and those farthest from TB diagnosis and treatment clinics, had substantially lower relative TB CNRs. This remained true after adjusting for important epidemiological risk factors for TB, implying inequity in access to diagnosis and care. The innovative georeferenced surveillance system described here could be implemented widely in low-resource settings and potentially provide high-quality, real-time epidemiological and microbiological data to direct TB control efforts. In particular, the impact on TB burden of facilitated access to diagnosis and spatially directed active TB case-finding interventions should be examined further.

Abbreviations

CHW:

Community health worker

CNR:

Case notification rate

ePAL:

Electronic patient locator

HIV:

Human immunodeficiency virus

MCMC:

Markov chain Monte Carlo

NUTS:

Non-U-turn sampler

TB:

Tuberculosis

WAIC:

Widely applicable information criteria

WHO:

World Health Organization

References

  1. World Health Organization. Global tuberculosis report 2017. World Health Organization Geneva: 2017 9789241565394 924156539X.

  2. Yates TA, Ayles H, Leacy FP, Schaap AB, Boccia D, Beyers N, et al. Socio-economic gradients in prevalent tuberculosis in Zambia and the Western Cape of South Africa. Tropical Med Int Health. 2018;23(4):375–90. https://doi.org/10.1111/tmi.13038 Epub 2018/02/13. PubMed PMID: 29432669.

    Article  Google Scholar 

  3. Lonnroth 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;68(12):2240–6. https://doi.org/10.1016/j.socscimed.2009.03.041 Epub 2009/04/28. PubMed PMID: 19394122.

    Article  PubMed  Google Scholar 

  4. Pedrazzoli D, Boccia D, Dodd PJ, Lonnroth K, Dowdy DW, Siroka A, et al. Modelling the social and structural determinants of tuberculosis: opportunities and challenges. Int J Tuberc Lung Dis. 2017;21(9):957–64. https://doi.org/10.5588/ijtld.16.0906 Epub 2017/08/23. PubMed PMID: 28826444; PubMed Central PMCID: PMCPMC5566999.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Storla DG, Yimer S, Bjune GA. A systematic review of delay in the diagnosis and treatment of tuberculosis. BMC Public Health. 2008;8:15. https://doi.org/10.1186/1471-2458-8-15 Epub 2008/01/16. PubMed PMID: 18194573; PubMed Central PMCID: PMCPMC2265684.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Tanimura T, Jaramillo E, Weil D, Raviglione M, Lonnroth K. Financial burden for tuberculosis patients in low- and middle-income countries: a systematic review. Eur Respir J. 2014;43(6):1763–75. https://doi.org/10.1183/09031936.00193413 Epub 2014/02/15. PubMed PMID: 24525439; PubMed Central PMCID: PMCPMC4040181.

    Article  PubMed  PubMed Central  Google Scholar 

  7. World Health Organization. END TB Strategy. Geneva: World Health Organization; 2016.

    Google Scholar 

  8. Middelkoop K, Bekker LG, Mathema B, Shashkina E, Kurepina N, Whitelaw A, et al. Molecular epidemiology of Mycobacterium tuberculosis in a South African community with high HIV prevalence. J Infect Dis. 2009;200(8):1207–11. https://doi.org/10.1086/605930 Epub 2009/09/22. PubMed PMID: 19764885; PubMed Central PMCID: PMCPMC2932637.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Dowdy DW, Golub JE, Chaisson RE, Saraceni V. Heterogeneity in tuberculosis transmission and the role of geographic hotspots in propagating epidemics. Proc Natl Acad Sci USA. 2012;109(24):9557–62. https://doi.org/10.1073/pnas.1203517109 Epub 2012/05/31. PubMed PMID: 22645356; PubMed Central PMCID: PMCPMC3386125.

    Article  PubMed  Google Scholar 

  10. ICAP, Columbia University, USA and Ministry of Health of Malawi. 2016. https://phia.icap.columbia.edu/wp-content/uploads/2016/09/MALAWI-Factsheet.FIN_.pdf.

  11. MacPherson P, Choko AT, Webb EL, Thindwa D, Squire SB, Sambakunsi R, et al. Development and validation of a global positioning system-based “map book” system for categorizing cluster residency status of community members living in high-density urban slums in Blantyre, Malawi. Am J Epidemiol. 2013;177(10):1143–7. https://doi.org/10.1093/aje/kws376 Epub 2013/04/17. PubMed PMID: 23589586; PubMed Central PMCID: PMCPMC3649638.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Harris RC. Informing development strategies for new tuber- culosis vaccines: mathematical modelling and novel epidemiological tools. PhD thesis. UK: London School of Hygiene and Tropical Medicine; 2018.

    Google Scholar 

  13. Pebesma E. sf: simple features for R. 2018. https://github.com/r-spatial/sf.

  14. Tatem A, Gething P, Pezzulo C, Weiss D, Bhatt S. Development of high-resolution gridded poverty surfaces, report for the Bill and Melinda Gates Foundation. 2014. http://www.worldpop.org.uk/resources/docs/Poverty-mapping-report.pdf.

  15. Gschlößl S, Czado C. Modelling count data with overdispersion and spatial effects. Stat Pap. 2008;49(3):531.

    Article  Google Scholar 

  16. Best N, Richardson S, Thomson A. A comparison of Bayesian spatial models for disease mapping. Stat Methods Med Res. 2005;14(1):35–59. https://doi.org/10.1191/0962280205sm388oa Epub 2005/02/05. PubMed PMID: 15690999.

    Article  PubMed  Google Scholar 

  17. Gelman A, Hwang J, Vehtari A. Understanding predictive information criteria for Bayesian models. Stat Comput. 2014;24(6):997–1016.

    Article  Google Scholar 

  18. Bürkner P-C. brms: an R package for Bayesian multilevel models using Stan. J Stat Softw. 2018;80(1):1–28.

    Google Scholar 

  19. Uplekar M, Weil D, Lonnroth K, Jaramillo E, Lienhardt C, Dias HM, et al. WHO’s new end TB strategy. Lancet. 2015;385(9979):1799–801. https://doi.org/10.1016/s0140-6736(15)60570-0 Epub 2015/03/31. PubMed PMID: 25814376.

    Article  PubMed  Google Scholar 

  20. Dodd PJ, White RG, Corbett EL. Periodic active case finding for TB: when to look? PLoS One. 2011;6(12):e29130. https://doi.org/10.1371/journal.pone.0029130 Epub 2012/01/05. PubMed PMID: 22216182; PubMed Central PMCID: PMCPMC3245256.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Corbett EL, Bandason T, Duong T, Dauya E, Makamure B, Churchyard GJ, et al. Comparison of two active case-finding strategies for community-based diagnosis of symptomatic smear-positive tuberculosis and control of infectious tuberculosis in Harare, Zimbabwe (DETECTB): a cluster-randomised trial. Lancet. 2010;376(9748):1244–53. https://doi.org/10.1016/s0140-6736(10)61425-0 Epub 2010/10/07. PubMed PMID: 20923715; PubMed Central PMCID: PMCPMC2956882.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Ayles H, Muyoyeta M, Du Toit E, Schaap A, Floyd S, Simwinga M, et al. Effect of household and community interventions on the burden of tuberculosis in southern Africa: the ZAMSTAR community-randomised trial. Lancet. 2013;382(9899):1183–94. https://doi.org/10.1016/s0140-6736(13)61131-9 Epub 2013/08/07. PubMed PMID: 23915882.

    Article  PubMed  Google Scholar 

  23. Ezeh A, Oyebode O, Satterthwaite D, Chen YF, Ndugwa R, Sartori J, et al. The history, geography, and sociology of slums and the health problems of people who live in slums. Lancet. 2017;389(10068):547–58. https://doi.org/10.1016/s0140-6736(16)31650-6 Epub 2016/10/21. PubMed PMID: 27760703.

    Article  PubMed  Google Scholar 

  24. Horton KC, MacPherson P, Houben RM, 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;13(9):e1002119. https://doi.org/10.1371/journal.pmed.1002119 Epub 2016/09/07. PubMed PMID: 27598345; PubMed Central PMCID: PMCPMC5012571.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Chikovore J, Hart G, Kumwenda M, Chipungu G, Desmond N, Corbett EL. TB and HIV stigma compounded by threatened masculinity: implications for TB health-care seeking in Malawi. Int J Tuberc Lung Dis. 2017;21(11):26–33. https://doi.org/10.5588/ijtld.16.0925 Epub 2017/10/14. PubMed PMID: 29025482.

    Article  CAS  PubMed  Google Scholar 

  26. Corbett EL, MacPherson P. Tuberculosis screening in high human immunodeficiency virus prevalence settings: turning promise into reality. Int J Tuberc Lung Dis. 2013;17(9):1125–38. https://doi.org/10.5588/ijtld.13.0117 Epub 2013/08/10. PubMed PMID: 23928165; PubMed Central PMCID: PMCPMC3966511.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the Wellcome Trust [206575 to PM, 200901 to ELC]. It was additional supported by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement [MR/K012126/1] to ELW. ELW is also part of the EDCTP2 programme supported by the European Union. PJD was supported by the UK Medical Research Council (MR/P022081/1). TC was supported by US National Institutes of Health (U54 GM088558). RH was funded by the UK Medical Research Council (MRC) under the London School of Hygiene and Tropical Medicine UK Medical Research Council vaccines scholarship program. The funders had no role in study design; in the collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the paper for publication.

Availability of data and materials

The datasets and code supporting the conclusions of this article are available in the GitHub repository (https://github.com/petermacp/BlantyreTBEpi), or installed as an R package by running ‘devtools::install_github(“petermacp/BlantyreTBEpi”).

Author information

Authors and Affiliations

Authors

Contributions

PM, ELC, RH and ELW designed the study. MK, VKP and RH contributed to the data collection and management. PM, MK and PJD contributed to the data analysis. PM, PJD, ELW, TC, ELC and MN contributed to the data interpretation. PM wrote the first draft. PM, MK, MN, ATC, VKP, ELW, PJD, TC, RH and ELC reviewed, contributed to, and approved the final draft. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Peter MacPherson.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was granted for the study by the research ethics committees of the College of Medicine, University of Malawi, and the London School of Hygiene and Tropical Medicine. All TB cases provided verbal informed consent for collection of clinical, demographic, and geospatial data. The Research Ethics Committees granted waivers of written informed consent because surveillance was administered through the routine public health service.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Figure S1. Characteristics of Health Surveillance Assistant catchment areas, Blantyre, Malawi. A) Total population; B) population density per km2; C) Mean number of people per household; D) Mean proportion of population living in poverty; E) Cartesian distance from catchment area centroid to nearest TB registration centre; F) Ratio of male-to-female adults G) Proportion of population aged 15 years or older H) Sputum smear positive to negative ratio, on sample from routine clinic. Poverty estimates from Worldpop (www.worldpop.org.uk). Boundaries are Health Surveillance Assistant catchment areas. Black triangles are TB registration clinic. White areas in centre of maps are mountainous areas or business districts with few residents that were not enumerated in the census. The black triangle in the far north of maps is a health centre with a TB registration clinic located outside of Blantyre District that may be used by Blantyre residents; to increase accuracy of case notification rates, we captured TB registrations at this clinic. (PDF 933 kb)

Additional file 2:

Figure S2. Nearest neighbour structure of Blantyre Health Surveillance Assistant catchment areas. Blue lines indicate the nearest k neighbours to each 315 Health Surveillance Assistant catchment area, calculated as the Cartesian distance between pairs of Health Surveillance Assistant catchment area centroids. (PDF 37 kb)

Additional file 3:

Figure S3. Sensitivity analysis. X-axes are relative rates per unit increase in variable. Model k0 is model without spatial autocorrelation prior. Models k1 to k6 are models with k = 1 to 6 nearest neighbour matrix spatial autocorrelation priors. k nearest neighbours calculated from Cartesian distances between pairs of Health Surveillance Assistant catchment area centroids (Additional file 2: Figure S2). WAIC: Widely applicable information criteria based on the posterior distribution of each model. (PDF 14 kb)

Additional file 4:

Figure S4. Predicted and observed TB cases in 315 Blantyre Health Surveillance Assistant catchment areas. A: All TB cases. B: Microbiologically confirmed TB cases. X-axes are Health Surveillance catchment areas, ordered in decreasing frequency by number of TB cases registered through enhanced TB surveillance. Red points are observed numbers of cases, blue points are predicted values, estimated by the Bayesian spatial regression models. Transparency of blue fitted points are inversely proportional to absolute distance from observed value. Predicted numbers of TB cases estimated by fitting a Bayesian spatial regression model with Poisson response, a k = 6 nearest-neighbours conditional spatial autocorrelation structure, with linear terms fitted for health surveillance assistant catchment area log10 total population, log10 population density, adult M:F ratio, mean number of people per household, log10 Cartesian distance from geographical centroid to the nearest TB clinic, proportion of population aged 15 years or older, mean proportion living on less than US $2 per day, sputum smear positive to negative ratio, offset term for log10 HSA population size, and with weakly informative prior on the population-level effects intercept (Gaussian: mean = 0, sd = 10), and predictor intercept (Gaussian, mean = 0, sd = 10). (PDF 109 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

MacPherson, P., Khundi, M., Nliwasa, M. et al. Disparities in access to diagnosis and care in Blantyre, Malawi, identified through enhanced tuberculosis surveillance and spatial analysis. BMC Med 17, 21 (2019). https://doi.org/10.1186/s12916-019-1260-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-019-1260-6

Keywords