Introduction
In 2019, tuberculosis (TB) was the leading cause of death from a single infectious disease in the world, causing more deaths than HIV/AIDS and malaria combined [1]. Slow progress has been made in reducing TB mortality, and the World Health Organization's (WHO) End TB strategy target of achieving a 90% reduction in TB deaths between 2015 and 2030 [1] is unlikely to be met.
In the WHO Africa (AFRO) region, TB treatment success rates (cured or completed treatment) in 2017 were 86% for HIV-negative people and 78% for HIV-positive people [1]. Despite these improvements from 2007 (when overall TB treatment success was estimated to be 79% in Africa), a considerable fraction of people with TB symptoms are delayed in accessing TB diagnosis and care at health facilities [1, 2]. The reasons for delayed access to treatment are complex and multi-layered, but include health care-seeking behaviours, clinical and geographical factors, and suboptimal care quality in health facilities [Reference Heunis3].
Distance to health facilities is a well-recognised access barrier to prompt diagnosis [Reference Mayer4–Reference O'Donnell7]. Individuals who live far from health facilities are at an increased risk of unfavourable health outcomes [Reference Mayer4–Reference Momenyan6] and at an increased risk of death on TB treatment [Reference Tripathy8–Reference Cai10]. However, there is limited evidence about associations between distance and death on TB treatment in urban African settings.
Using prospectively collected data from people initiating TB treatment at health facilities in urban Blantyre, Malawi in the era of high antiretroviral therapy (ART) coverage for HIV, we hypothesised that people initiating TB care at facilities at greater distances from their homes might be at greater risk of death on TB treatment compared to people who lived nearer to clinics. Additionally, as accurate measurement of clinic distance may be challenging under routine programmatic conditions, we compared distance measurement using two approaches that could be used by health planners and epidemiologists: Cartesian distance and shortest road network distance.
Methods
Study setting and design
Blantyre District is located in the southern region of Malawi. Blantyre city, in the centre of the district, has large areas of densely-populated informal settlements. Blantyre's 2018 census population was 1 264 304 [11] and HIV prevalence was 18% [12].
Blantyre enhanced TB monitoring and evaluation
In Malawi, patients diagnosed with TB register to receive treatment at primary health care centres and hospitals. In Blantyre, TB registration clinics include one referral hospital (Queen Elizabeth Central Hospital (QECH)), one large private church-supported clinic (Mlambe), three private clinics and seven government public primary health care clinics [Reference Makwero13]. QECH is a tertiary referral hospital for the whole southern region of Malawi and offers inpatient care including TB diagnosis and treatment [Reference Makwero13].
In a joint project between the Malawi-Liverpool-Wellcome Trust Clinical Research Programme, Blantyre District Health Office and the Malawi National TB Control Programme (NTP), TB Officers (a cadre of health worker employed by the Ministry of Health of Malawi) received training to strengthen the TB surveillance system in Blantyre. Since 2015, TB Officers stationed at TB treatment clinics in Blantyre provided TB treatment, HIV testing and linkage to treatment in accordance with Malawi guidelines. They further recorded TB registration data into NTP TB registers and additionally recorded individual-level clinical, socio-demographic and household data using an electronic data collection application (ePaL), which we have previously developed and validated (18). Using ePaL, TB officers obtained global positioning satellite (GPS) coordinates for the location of TB patients’ households (18). A spot sputum sample was also collected from all patients starting TB treatment (if patients were able to produce sputum) that underwent smear microscopy and mycobacteria growth indicator tube culture. On a quarterly basis, study and TB programme registers were reconciled.
TB outcomes were defined according to mutually-exclusive WHO TB treatment outcome definition guidelines [14], with patients classified as either cured, completed treatment, treatment failed, confirmed died, lost to follow-up or not evaluated (typically because still on treatment). The outcome of interest for this analysis was confirmed death on TB treatment; we did not follow-up patients to confirm death. The analysis was limited to study participants that were registered for treatment from 1 January 2015 to 30 December 2018.
Statistical analysis
Baseline characteristics
We compared the characteristics of participants, recorded at TB treatment registration, between those who died during TB treatment and those who were alive at the end of treatment. For categorical variables, we calculated percentages and used the χ 2 test for comparisons; for continuous variables, we calculated medians (and interquartile ranges (IQR)), mean (and standard deviation) and compared between groups using the Kruskal–Wallis test.
Distance estimation
We estimated the distance from study participants’ households to their TB treatment initiation clinic using two approaches (Fig. 1). The first approach was to estimate the distance based on a ‘straight line’ distance (Cartesian distance). In the second approach, we used Blantyre urban road network downloaded from OpenStreetMap (OpenStreetMap Foundation) to calculate the shortest road network distance using the stplanr R package [Reference Lovelace and Ellison15].
Statistical models
To estimate the causal relationship between clinic distance and TB case fatality, we constructed a directed acyclic graph (Fig. 2). The minimum adjustment set of confounders identified from the DAG were: sex, age, HIV, hospital admission (QECH) and household wealth score. We included a term for the interaction between household to clinic distance and the clinic at which TB registration occurred (QECH, the city's referral hospital, vs. other primary health care centres). This was because QECH is a large referral hospital known to have a patient population with more advanced TB disease, more complicated TB disease and patients that are likely to have travelled a longer distance compared to patients that attend other clinics.
Using asset ownership data, we created a household wealth score variable, calculated using a proxy means test developed for urban populations from the Malawi Integrated Household Survey [Reference Payongayong16]. Individuals for whom components needed to calculate their wealth score were missing; their wealth score was imputed using multiple interval imputation, with the observed wealth score providing the lower bound.
We constructed Bayesian multi-level logistic regression models (Equation 1) to investigate the association between participants’ household to clinic distance and their risk of death, adjusting for confounders and an interaction between household to clinic distance and clinic of registration. Weakly regularising priors, refined by inspecting prior predictive plots, were assigned to model intercepts and slopes. Model convergence was assessed by visual inspection of trace plots, effective sample numbers and Gelman–Rubin statistics.
where Dij indicates death on TB treatment for the i th patient in clinic j, sex indicates sex (male or female), age indicates age in years, hivstatus indicates HIV status (positive or negative), wealthscore indicates household wealth score, distance indicates distance to clinic in kilometres, and QECH indicates whether participant registered for treatment at QECH or the other primary health centres, distance×QECH is the interaction effect between distance and QECH and clinic of registration indicates registering for treatment at the j th TB treatment registration health centre.
Four thousand samples with a warm up of 600 iterations were drawn from each model posterior distribution using Markov Chain Monte Carlo methods implemented within Stan via the brms R package [Reference Bürkner17]. Posterior means and 95% credible intervals (CI) were calculated on the log-odds scale and were exponentiated to give odds ratios. In sensitivity analysis, we refitted models recoding participants who were lost to follow-up or transferred out during TB treatment as having died since it was possible that some of these patients had actually died [Reference Squire18]. Additionally, we restricted the analysis to those with microbiologically-confirmed TB. To investigate the predictive accuracy of the two distance estimates, we used models that just had distance and the interaction effect of distance and registering for TB treatment at QECH. We plotted the difference in the predicted probability of death for the 100 patients with the greatest difference in distance from the two measurement methods. Analysis was conducted using R version 3.5.2 (R Foundation for Statistical Computing, Vienna).
Ethical considerations
Ethical approval was granted by the London School of Hygiene and Tropical Medicine (16228) and the College of Medicine, University of Malawi Research Ethics Committee (P.12/18/2556). Participants provided oral consent to participate in TB surveillance with a waiver for written consent granted by both research ethics committees.
Results
Baseline characteristics
A total of 4461/5199 (85.8%) patients that initiated TB treatment at the 12 study clinics were included in this analysis. In total, 64/4461 (1.4%) had a missing TB outcome and were excluded from the primary analysis (Appendix Table A1). Among participants that had treatment outcomes, 479/4397 (10.9%) died while taking TB treatment, and 258/4397 (5.9%) were reported as lost to follow-up or transferred out (Table 1). The percentage of participants who died did not substantially differ across the study years (P = 0.313).
HIV, human immunodeficiency virus; km, kilometres; (Q1, Q3) interquartile range; SD, standard deviation; TB, tuberculosis.
a Wealth score, household wealth developed using household asset ownership.
More deaths were reported among participants who initiated TB treatment at QECH (315/2170, 14.5%) than in the other 11 TB clinics (164/2227, 7.4%, P < 0.001; Table 1). Similar percentages of men (303/2789, 10.9%) and women (176/1608, 10.9%, P = 0.934) died during TB treatment. The median age of participants who died was older (median: 37 years, IQR: 30.0–45.0) compared to those who were alive at the end of TB treatment (35 years, IQR: 28.0–41.0, P < 0.001). A higher proportion of deaths occurred in patients who were HIV-positive (400/2981, 13.4%) compared to HIV-negative (79/1416, 5.6%, P < 0.001). Participants who did not have a microbiologically-confirmed TB (318/2382, 13.4%) were considerably more likely to die compared to participants with either sputum smear, Xpert or culture-positive disease (161/2015, 8.0%, P < 0.001).
The median road network household to clinic distance (median: 5.2 km, IQR: 3.2–7.7) was consistently higher than the Cartesian household to clinic distance estimates (median: 3.8 km, IQR: 2.1–5.9). Both household to clinic distance estimates showed that patients who died on TB treatment lived further from clinics where they initiated TB treatment than those that were alive at the end of TB treatment. For the road network distance measure, those who died on TB treatment had greater clinic distances (median: 5.8 km, IQR: 4.1–9.0) compared to participants who were still alive at TB treatment completion (median: 5.1 km, IQR: 3.1–7.5, P < 0.001). Likewise, those who died on TB treatment had a longer Cartesian distance (median: 4.4 km, IQR: 2.9–6.8) vs. participants who were still alive at the end of TB treatment (median: 3.8 km, IQR: 2.0–5.7, P < 0.001). Patients who did not live within geographically-mapped areas of the city were excluded 738/5199 (14.2%) because we did not have data to ascertain their household location, and so could not calculate their household to clinic distance (Appendix Table A1).
Unadjusted analysis of the association of household to clinic distance and risk of death
TB patients who lived further from the clinic at which they registered for treatment had higher odds of death compared to those who lived nearer, whether measured by the road network distance method (odds ratio (OR) 1.05 per km increase, 95% CI 1.03–1.08) or the Cartesian distance method (OR 1.06 per km increase, 95% CI 1.03–1.09; Table 2).
km, Kilometres; OR, odds ratio.
The analysis was done for all the data at once (N = 4397) using Bayesian multi-level logistic regression models.
a,b Adjusted for sex (male vs. female), age in years, HIV status (HIV-positive vs. HIV-negative), wealth score the household wealth score, TB treatment registration at Queen Elizabeth Central Hospital (QECH) vs. the other health care facilities, distance household to clinic distance in kilometres, distance (km)×QECH interaction effect of QECH and household to clinic distance and a term for the random intercept of clinic of registration.
c Age in years, household to clinic distance and household wealth score were centred by subtracting the mean.
d Wealth score of household wealth developed using household asset ownership.
e,f Adjusted for TB treatment registration at QECH vs. other health care facilities, household to clinic distance, and the interaction effect of QECH and household to clinic distance.
Unadjusted analysis for confounders
In unadjusted analysis, each one year increase in age was associated with a 2% increase (OR 1.02, 95% CI 1.01–1.02) in the odds of death. Being HIV-positive was associated with a 2.5-times (OR vs. HIV-negative status 2.53, 95% CI 2.0–3.25) increase in the odds of death on TB treatment. TB patients who registered at QECH had two times higher (OR 2.10, 95% CI 1.74–2.56) odds of death compared to people registering at the other primary clinics. Each one-unit increase in wealth score was associated with a 40% (OR 1.40, 95% CI 0.98–2.04) increase in the odds of death.
Adjusted analysis for household to clinic distance models
Greater household to clinic distance (road network) had a higher odds of death in TB patients registered at QECH (OR 1.07 per km increase, 95% CI 0.99–1.16), but was not associated with odds of death in TB patients registered at other clinics (OR 0.98 per km increase, 95% CI 0.92–1.03). The average risk of death varied between clinics by 6% (95% CI 4–11%) holding categorical variables at their base level and continuous variables at their average value, with a standard deviation of 0.67 (0.32–1.18; Table 2).
For the model containing Cartesian distance to clinic, model coefficients and uncertainty bounds were very similar to those from the road distance model, with distance from household to clinic of treatment registration associated with higher odds of death in TB patients registered at QECH (OR 1.09 per km increase, 95% CI 1.0–1.21), but was not associated with odds of death on TB treatment in patients registered at other clinics (OR 0.97 per km increase, 95% CI 0.89–1.04).
Comparison of model predicted risk of death among 100 participants with biggest difference between their road network distance and Cartesian distance estimates
The majority of the 100 TB patients with the greatest difference in distance between the two household to clinic distance measurements registered for treatment at QECH or Mlambe private clinic (Fig. 3). For these participants who had the greatest measurement difference, there was, as expected, a trend of increased difference in probability of death with a greater difference between the two distance measures. The greatest differences were seen at Mlambe private clinic, reflecting greater distances travelled by road. However, the maximum difference in the predicted probability of death over these 100 patients was only 0.79%.
Sensitivity analysis
When we reclassified TB patients who were lost to follow-up or had transferred out during TB treatment as having died, model estimates were very similar to our primary analysis (Table 3). Analysis restricted to participants with microbiologically-confirmed disease did not alter our model conclusion (Table 4).
km, Kilometres; OR, odds ratio.
The analysis was done for all the data at once (N = 4397) using Bayesian multi-level logistic regression models with participants that defaulted and those that transferred out (N = 258) recoded as died.
a,b Adjusted for sex (male vs. female), age in years, HIV status (HIV-positive vs. HIV-negative), wealth score the household wealth score, TB treatment registration at Queen Elizabeth Central Hospital (QECH) vs. other health care facilities, distance household to clinic distance in kilometres, distance (km)×QECH interaction effect of QECH and household to clinic distance and a term for the random intercept of clinic of registration.
c Age in years, household to clinic distance and household wealth score were centred by subtracting the mean.
d Wealth score of household wealth developed using household asset ownership.
e,f Adjusted for TB treatment registration at QECH vs. other health care facilities, household to clinic distance, and the interaction effect of QECH and household to clinic distance.
km, Kilometres; OR, odds ratio.
The analysis was done using N = 2015 participants that had a microbiologically confirmed TB diagnosis using Bayesian multi-level logistic regression models.
a,b Adjusted for sex (male vs. female), age in years, HIV status (HIV-positive vs. HIV-negative), wealth score the household wealth score, TB treatment registration at Queen Elizabeth Central Hospital (QECH) vs. the other health care facilities, household to clinic distance in kilometres, distance (km)×QECH interaction effect of QECH and household to clinic distance and a term for the random intercept of clinic of registration.
c Age in years, household to clinic distance and household wealth score were centred by subtracting the mean.
d Wealth score of household wealth developed using household asset ownership.
e,f Adjusted for TB treatment registration at QECH vs. other health care facilities, household to clinic distance, and the interaction effect of QECH and household to clinic distance
Discussion
Our primary hypothesis was that greater household to clinic distance would be associated with an increased risk of death for people starting treatment for TB in urban Blantyre, Malawi. The results of the pre-specified unadjusted analysis showed clinic distance was a risk factor for death after starting TB treatment. However, in adjusted analysis based on our causal graph, greater clinic distance was a risk factor for death only for participants initiating care at the tertiary referral hospital (QECH). Treatment initiation at the city's tertiary hospital, older age and HIV-positive status were important predictors of death whilst taking TB treatment. We additionally found moderate variation in the risk of death between treatment clinics that was not explained by TB patient clinical and socio-demographic characteristics, indicating that clinic quality of care might be an important residual determinant of treatment outcomes.
We used two different measures (Cartesian distance, and shortest road network distance) of household to clinic distance that could be applied by public health planners and epidemiologists, recognising that road network distance may be considerably more challenging to estimate under routine programmatic conditions. Our analysis found that models using either distance measurement reached very similar conclusions, and therefore we recommend that in low-resource settings, the simpler measurement (Cartesian distance) should likely be more appropriate for routine programme use.
For TB patients that were treated at the city's tertiary hospital (QECH), greater household to clinic distance might be a proxy for lack of access to quality health care service [Reference Tripathy8, Reference Laux and Patil19]. Individuals who live in the peri-urban areas far from the centre of Blantyre have limited resources for transport, and may struggle to access timely tertiary-level health care compared to those living nearer to the centre of the city [Reference MacPherson20]. Delayed presentation for care will then result in admission to hospital with more advanced disease and at a higher risk of death [5, Reference Virenfeldt21].
To reduce the high risk of death on TB treatment, prompt diagnosis and treatment for people presenting to health facilities is required [Reference Virenfeldt21, 22]. At health facilities in high HIV-prevalence settings, TB symptoms are common (up to 60% reporting at least one of the WHO four cardinal TB symptoms) [Reference Chihota23]. Therefore, novel screening approaches – including triage testing using high sensitivity initial tests such as chest X-ray or C-reactive protein, followed by highly specific tests such as Xpert for those who triage test positive, require evaluation [5, Reference van't Hoog24].
QECH is a tertiary referral hospital, and most patients treated for TB at QECH will have previously sought diagnosis – likely on multiple occasions before referral – at a primary health centre [Reference Chikovore25]. Consequently, earlier diagnosis in primary health care could reduce hospital admission and potentially reduce mortality [1]. We need to improve the TB diagnostic capacity of primary health facilities in the outskirts of the city and ensure that primary health care facilities can make prompt referrals for hospital care where required [1, 5, Reference Tripathy8].
In addition, there is a need to evaluate community interventions of TB screening in areas further away from the centre of the city (QECH) since prevalent cases identified in community interventions are usually diagnosed at an earlier stage of disease [1, 5]. Our previous work has shown that TB case notifications decline with greater distance from the centre of Blantyre city; this might be because of lack of access to quality health care, including TB diagnosis, at areas in the periphery of the city [Reference MacPherson20]. Health promotion activities to promote early treatment seeking, and quality improvement activities with health workers to support early diagnosis and reduce variation in practice between clinics could contribute to reduced TB case fatality, whilst strengthening universal healthcare provision within health facilities [1, 5].
Malawi has made tremendous progress towards achieving the WHO 90-90-90 HIV targets [26]: 90% of individuals with HIV are aware of their HIV status, 79% of individuals with a diagnosis are on ART and 72% of individuals on ART treatment have viral suppression [27]. Nevertheless, in this study, HIV-positive status remained strongly associated with an increased risk of death during TB treatment. Nearly all TB patients in this study were aware of their HIV-positive status, and 89% were taking ART at TB treatment initiation. Persistently high case fatality for HIV-positive people despite high ART coverage suggests that HIV-positive people have severe acute illness [Reference Kwan and Ernst28], have profound delays in TB diagnosis and treatment initiation [Reference Nogueira29, Reference Gesesew30], or have virological failure to HIV treatment, which itself confers a high risk of death. In a study among HIV-positive people taking ART who were admitted to hospital in Zomba District, Malawi (32%) had virological treatment failure, and resistance to first-line ART drugs was near-universal [Reference Gupta-Wright31]. Improvements in viral load monitoring and rapid treatment changes linked to supportive adherence interventions could reduce the number of patients that present with severe disease at the start of TB treatment and save more lives [Reference Gupta-Wright31, Reference Craik32]. Additionally, implementation of point-of-care HIV viral load monitoring for people living with HIV who are admitted to hospital could identify treatment failure earlier. Following WHO recommendations, the Malawi National HIV Programme has recommended that all people living with HIV take a Dolutegravir containing regimen, including by switching treatment regimens for those already taking [Reference Gupta-Wright31, 33]. It will be important to investigate the effect of this change on TB case fatality among patients treated at primary clinics and in hospital (QECH) [33].
A key strength of this study was that we used prospectively collected data at TB registration through our citywide enhanced TB surveillance system. Our routine monitoring and evaluation show consistently very high agreement with national TB treatment registers. We additionally captured TB patients’ household GPS coordinates at the point of registration for TB treatment; most previous studies have attempted to retrospectively geolocate patients’ households using physical addresses, a method that is prone to error [Reference Cegielski34, Reference Robsky35]. Our analysis of household to clinic distance was analysed at a continuous scale rather than on a transformed categorical scale with arbitrary value cut-off points [Reference Royston, Altman and Sauerbrei36, Reference Altman and Royston37]. To investigate the causal relationship between household to clinic distance and risk of death, we selected confounding variables for adjustment by constructing a directed acyclic graph [Reference Textor38]. Sensitivity analysis restricted to TB patients with microbiologically-confirmed TB and where patients who were lost to follow-up or transferred out to another clinic did not alter our findings. However, timing of death was not collected meaning survival analysis could not be done. A small number of TB patients may not have been able to accurately geolocate their household, or may have deliberately misidentified their household; monthly quality assurance checking of a 5% random sample of households attempted to mitigate this. While Appendix Table A1 indicates that TB patients registering for treatment who were residents in the areas of the city not mapped by our GPS had some different characteristics compared to those living in the mapped areas. Nevertheless, the main objective of this analysis was to predict TB case fatality using data obtained at treatment registration and there is no reason to believe that the underlying relationships should differ between those included and excluded.
In conclusion, in prespecified multilevel modelling using citywide enhanced surveillance data linked to our satellite GPS location system, we found that household to clinic distance was associated with increased risk of death for a patient starting TB treatment at the city's central hospital. This could be explained by barriers in accessing prompt diagnosis and treatment, and by variable quality of care at primary health facilities. HIV-positive status and older age remain important risk factors for death. Therefore, interventions that improve access to TB diagnosis through community-based active case finding, and improved quality of health facility TB screening, and prioritise HIV-positive people with potentially high levels of viral failure for TB screening and ART optimisation are required to reduce the unacceptably high case fatality.
Financial support
This work was supported by two grants from Wellcome Trust (E.L.C. grant number WT200901/Z/16/Z) and (P.M. grant number 206575). J.R.C is funded by MRC programme grant MC_UU_00004/07.
Conflict of interest
None.
Data availability statement
Data will be made available upon request.