Hostname: page-component-cd9895bd7-7cvxr Total loading time: 0 Render date: 2024-12-27T08:38:53.942Z Has data issue: false hasContentIssue false

Clinical, health systems and neighbourhood determinants of tuberculosis case fatality in urban Blantyre, Malawi: a multilevel epidemiological analysis of enhanced surveillance data

Published online by Cambridge University Press:  02 September 2021

McEwen Khundi*
Affiliation:
Malawi-Liverpool-Wellcome Trust Clinical Research Programme, Blantyre, Malawi London School of Hygiene and Tropical Medicine, London, UK
Peter MacPherson
Affiliation:
Malawi-Liverpool-Wellcome Trust Clinical Research Programme, Blantyre, Malawi London School of Hygiene and Tropical Medicine, London, UK Department of Clinical Sciences, Liverpool School of Tropical Medicine, Liverpool, UK
Helena R. A. Feasey
Affiliation:
Malawi-Liverpool-Wellcome Trust Clinical Research Programme, Blantyre, Malawi London School of Hygiene and Tropical Medicine, London, UK
Rebeca Nzawa Soko
Affiliation:
Malawi-Liverpool-Wellcome Trust Clinical Research Programme, Blantyre, Malawi London School of Hygiene and Tropical Medicine, London, UK
Marriott Nliwasa
Affiliation:
Malawi-Liverpool-Wellcome Trust Clinical Research Programme, Blantyre, Malawi Helse Nord TB Initiative, College of Medicine, University of Malawi, Blantyre, Malawi
Elizabeth L. Corbett
Affiliation:
Malawi-Liverpool-Wellcome Trust Clinical Research Programme, Blantyre, Malawi London School of Hygiene and Tropical Medicine, London, UK
James R. Carpenter
Affiliation:
London School of Hygiene and Tropical Medicine, London, UK
*
Author for correspondence: McEwen Khundi, E-mail: mcewenkhundi@gmail.com
Rights & Permissions [Opens in a new window]

Abstract

We investigated whether household to clinic distance was a risk factor for death on tuberculosis (TB) treatment in Malawi. Using enhanced TB surveillance data, we recorded all TB treatment initiations and outcomes between 2015 and 2018. Household locations were geolocated, and distances were measured by a straight line or shortest road network. We constructed Bayesian multi-level logistic regression models to investigate associations between distance and case fatality. A total of 479/4397 (10.9%) TB patients died. Greater distance was associated with higher (odds ratio (OR) 1.07 per kilometre (km) increase, 95% credible interval (CI) 0.99–1.16) odds of death in TB patients registered at the referral hospital, but not among TB patients registered at primary clinics (OR 0.98 per km increase, 95% CI 0.92–1.03). Age (OR 1.02 per year increase, 95% CI 1.01–1.02) and HIV-positive status (OR 2.21, 95% CI 1.73–2.85) were also associated with higher odds of death. Model estimates were similar for both distance measures. Distance was a risk factor for death among patients at the main referral hospital, likely due to delayed diagnosis and suboptimal healthcare access. To reduce mortality, targeted community TB screening interventions for TB disease and HIV, and expansion of novel sensitive diagnostic tests are required.

Type
Original Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2021. Published by Cambridge University Press

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 Mayer4Reference O'Donnell7]. Individuals who live far from health facilities are at an increased risk of unfavourable health outcomes [Reference Mayer4Reference Momenyan6] and at an increased risk of death on TB treatment [Reference Tripathy8Reference 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].

Fig. 1. Illustration of the two methods for the measurement of household to clinic distance: Cartesian distance (2.7 km), shortest road network distance (4.5 km). Note: The patient place of residence is a randomly generated point for illustration and does not correspond to any patient in the dataset.

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.

Fig. 2. Directed acyclic graph (causal diagram). Illustrating the relationship between distance to TB clinic and risk of death on TB treatment and other covariates. The variables sex, age, HIV, Queens Elizabeth hospital registration (QECH) vs. registration at other clinics and poverty were selected as the minimum adjustment set.

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.

(1)$$\matrix{D_{ij}\sim Binomial( {1, \;p_{ij}} ) \cr logit( {\,p_{ij}} ) = \alpha _{clinic\;of\;registration\_j} + \beta _{sex\_ij} + + \beta _{age\_ij} + \cr \beta _{wealthscore\_ij\;} + \beta _{distance\_ij\;} + \beta _{QECH\_ij} + \beta _{distance_{ij}\ast QECH\_ij} \cr {\rm PRIORS} \cr \beta _{sex}\sim Normal( {0, \;\;0.5} ) \cr \beta _{age}\sim Normal( {0, \;\;0.5} ) \cr \beta _{hivstatus}\sim Normal( {0, \;\;1} ) \cr \beta _{wealthscore}\sim Normal( {0, \;\;0.5} ) \cr \beta _{distance}\sim Normal( {0, \;\;0.5} ) \cr \beta _{{\rm QECH}}\sim Normal( {0, \;\;0.5} ) \cr \beta _{distance\ast {\rm \;QECH}\;}\sim Normal( {0, \;\;0.5} ) \cr \beta _{clinic\;of\;registration\_j\;}\sim \alpha + u_j \cr \alpha \;\sim \;Normal( {0, \;\;0.5} ) \cr u_j\sim Normal( {0, \;\sigma_\alpha } ) \cr \sigma _\alpha \sim HalfNormal( {0, \;\;0.5} ) }$$

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

Table 1. Characteristics for notified tuberculosis cases in urban Blantyre Malawi, 2015–2018

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

Table 2. Statistical model results for the main analysis, notified TB cases in urban Blantyre, Malawi from 2015–2018

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

Fig. 3. Plot of difference between shortest road network and Cartesian distance, vs. fitted probability of death, for the 100 notified TB cases with the largest distance differences. q, Queen Elizabeth Central Hospital; m, Mlambe Private Hospital; l, Limbe Health Centre; g, Bangwe Health Centre; z, Zingwangwa Health Centre; a, Blantyre Adventist Hospital; c, Chilomoni Health Centre.

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

Table 3. Statistical model results for the sensitivity analysis: patients with loss to follow-up or transfer out treatment status recoded as having died, in urban Blantyre, Malawi, 2015–2016

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.

Table 4. Statistical model results for the sensitivity analysis: restricted to patients with microbiologically-confirmed tuberculosis disease only in urban Blantyre Malawi, 2015–2018

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.

Appendix

Table A1. Characteristics for notified Tuberculosis cases in Urban Blantyre Malawi, 2015–2018. Comparing participants who were included in the study versus those who were not part of the study but were from Blantyre.

Footnotes

Abbreviations: HIV, Human immunodeficiency virus; (Q1, Q3) interquartile range; SD standard deviation; TB tuberculosis

References

World Health Organization (2019) Global Tuberculosis Report 2018. WHO. Geneva.Google Scholar
World Health Organization (2010) Global Tuberculosis Control Report 2009. Geneva.Google Scholar
Heunis, JC, et al. (2017) Risk factors for mortality in TB patients: a 10-year electronic record review in a South African province. BMC Public Health 17, 38.Google Scholar
Mayer, CM, et al. (2019) Distance to clinic is a barrier to PrEP uptake and visit attendance in a community in rural Uganda. Journal of the International AIDS Society 22, e25276.CrossRefGoogle Scholar
Wold Health Organisation (2013) Systematic Screening for Active Tuberculosis. Geneva.Google Scholar
Momenyan, S, et al. (2018) Spatial inequalities and predictors of HIV/AIDS mortality risk in Hamadan, Iran: a retrospective cohort study. Epidemiology and Health 40, e2018038.CrossRefGoogle ScholarPubMed
O'Donnell, O (2007) Access to health care in developing countries: breaking down demand side barriers. Cadernos de Saude Publica 23, 28202834.CrossRefGoogle ScholarPubMed
Tripathy, JP, et al. (2013) Is physical access an impediment to tuberculosis diagnosis and treatment? A study from a rural district in North India. Public Health Action 3, 235239.CrossRefGoogle ScholarPubMed
Fluegge, K, et al. (2018) Impact of geographic distance on appraisal delay for active TB treatment seeking in Uganda: a network analysis of the Kawempe Community Health Cohort Study. BMC Public Health 18, 798.Google ScholarPubMed
Cai, J, et al. (2015) Factors associated with patient and provider delays for tuberculosis diagnosis and treatment in Asia: a systematic review and meta-analysis. PLoS ONE 10, e0120088. Published online: 25 March 2015. doi: 10.1371/journal.pone.0120088.Google ScholarPubMed
National Statistical Office of Malawi (NSO) (2018) Malawi Population Census Preliminary Report. Zomba.Google Scholar
Malawi Ministry of Health. Malawi Population-Based HIV Impact Assessment (MPHIA) 2015–2016: Final Report. Lilongwe, Ministry of Health. October 2018.Google Scholar
Makwero, MT (2018) Delivery of primary health care in Malawi. African Journal of Primary Health Care and Family Medicine 10, 1799.Google ScholarPubMed
World Health Organization (WHO) (2014) TB case and treatment outcome definitions. World Health Organization; Published online: 2014.Google Scholar
Lovelace, R and Ellison, R (2018) stplanr: A Package for Transport Planning. The R Journal 10, 723.CrossRefGoogle Scholar
Payongayong, E, et al. (2006) Simple household poverty assessment models for Malawi: proxy means test from the 1997–98 Malawi Integrated Household Survey. Lilongwe, Malawi: National Statistical Office, Government of Malawi.Google Scholar
Bürkner, P-C (2017) brms: an R package for Bayesian multilevel models using stan. Journal of Statistical Software 80, 128.CrossRefGoogle Scholar
Squire, SB, et al. (2005) ‘Lost’ smear-positive pulmonary tuberculosis cases: where are they and why did we lose them? International Journal of Tuberculosis and Lung Disease 9, 2531.Google Scholar
Secretary of Jan Swasthya Sahyog, Laux, TS and Patil, S (2018) Predictors of tuberculosis treatment outcomes among a retrospective cohort in rural, Central India. Journal of Clinical Tuberculosis and Other Mycobacterial Diseases 12, 4147.Google Scholar
MacPherson, P, et al. (2019) Disparities in access to diagnosis and care in Blantyre, Malawi, identified through enhanced tuberculosis surveillance and spatial analysis. BMC Medicine 17, 21.CrossRefGoogle ScholarPubMed
Virenfeldt, J, et al. (2014) Treatment delay affects clinical severity of tuberculosis: a longitudinal cohort study. BMJ Open 4, 4818.CrossRefGoogle ScholarPubMed
World Health Organization (2013) Systematic screening for active tuberculosis. Principles and recommendations. Geneva.Google Scholar
Chihota, VN, et al. (2015) Missed opportunities for TB investigation in primary care clinics in South Africa: experience from the XTEND Trial. PLoS ONE 10, e0138149.CrossRefGoogle ScholarPubMed
van't Hoog, AH, et al. (2012) Screening strategies for tuberculosis prevalence surveys: the value of chest radiography and symptoms. PLoS ONE 7, e38691.Google ScholarPubMed
Chikovore, J, et al. (2015) ‘For a mere cough, men must just chew Conjex, gain strength, and continue working’: the provider construction and tuberculosis care-seeking implications in Blantyre, Malawi. Global Health Action 8, e26292. Published online: 2015. doi: 10.3402/gha.v8.26292.CrossRefGoogle ScholarPubMed
United Nations Joint Programme on HIV/AIDS (UNAIDS) (2020) Joint United Nations Programme on HIV/AIDS. UNAIDS data 2020. Geneva, Switzerland; 436.Google Scholar
United Nations Joint Programme on HIV/AIDS (UNAIDS). Malawi country HIV/AIDS factsheet 2019. Available at https://www.unaids.org/en/regionscountries/countries/malawi (Accessed 29 October 2020).Google Scholar
Kwan, C and Ernst, JD (2011) HIV and tuberculosis: a deadly human syndemic. Clinical Microbiology Reviews 24, 351376.CrossRefGoogle ScholarPubMed
Nogueira, BMF, et al. (2018) Factors associated with tuberculosis treatment delay in patients co-infected with HIV in a high prevalence area in Brazil. PLoS ONE 13, e0195409. Published online: 1 April 2018. doi: 10.1371/journal.pone.0195409.Google Scholar
Gesesew, H, et al. (2016) The prevalence and associated factors for delayed presentation for HIV care among tuberculosis/HIV co-infected patients in Southwest Ethiopia: a retrospective observational cohort. Infectious Diseases of Poverty 5, 96.Google ScholarPubMed
Gupta-Wright, A, et al. (2020) Virological failure, HIV-1 drug resistance, and early mortality in adults admitted to hospital in Malawi: an observational cohort study. The Lancet HIV8 7, e620e628.CrossRefGoogle Scholar
Craik, A, et al. (2016) Challenges with targeted viral load testing for medical inpatients at Queen Elizabeth Central Hospital in Blantyre, Malawi. Malawi Medical Journal 28, 179181.CrossRefGoogle Scholar
World Health Organisation (2019) Update of recommendations on first-and second-line antiretroviral regimens. Geneva.Google Scholar
Cegielski, JP, et al. (2013) Eliminating tuberculosis one neighborhood at a time. American Journal of Public Health 103, 12921300.CrossRefGoogle ScholarPubMed
Robsky, KO, et al. (2020) Is distance associated with tuberculosis treatment outcomes? A retrospective cohort study in Kampala, Uganda. BMC Infectious Diseases 20, 406.CrossRefGoogle ScholarPubMed
Royston, P, Altman, DG and Sauerbrei, W (2006) Dichotomizing continuous predictors in multiple regression: a bad idea. Statistics in Medicine 25, 127141.CrossRefGoogle ScholarPubMed
Altman, DG and Royston, P (2006) The cost of dichotomising continuous variables. British Medical Journal 332, 1080.CrossRefGoogle ScholarPubMed
Textor, J, et al. (2016) Robust causal inference using directed acyclic graphs: the R package ‘dagitty’. International Journal of Epidemiology 45, 18871894.Google Scholar
Figure 0

Fig. 1. Illustration of the two methods for the measurement of household to clinic distance: Cartesian distance (2.7 km), shortest road network distance (4.5 km). Note: The patient place of residence is a randomly generated point for illustration and does not correspond to any patient in the dataset.

Figure 1

Fig. 2. Directed acyclic graph (causal diagram). Illustrating the relationship between distance to TB clinic and risk of death on TB treatment and other covariates. The variables sex, age, HIV, Queens Elizabeth hospital registration (QECH) vs. registration at other clinics and poverty were selected as the minimum adjustment set.

Figure 2

Table 1. Characteristics for notified tuberculosis cases in urban Blantyre Malawi, 2015–2018

Figure 3

Table 2. Statistical model results for the main analysis, notified TB cases in urban Blantyre, Malawi from 2015–2018

Figure 4

Fig. 3. Plot of difference between shortest road network and Cartesian distance, vs. fitted probability of death, for the 100 notified TB cases with the largest distance differences. q, Queen Elizabeth Central Hospital; m, Mlambe Private Hospital; l, Limbe Health Centre; g, Bangwe Health Centre; z, Zingwangwa Health Centre; a, Blantyre Adventist Hospital; c, Chilomoni Health Centre.

Figure 5

Table 3. Statistical model results for the sensitivity analysis: patients with loss to follow-up or transfer out treatment status recoded as having died, in urban Blantyre, Malawi, 2015–2016

Figure 6

Table 4. Statistical model results for the sensitivity analysis: restricted to patients with microbiologically-confirmed tuberculosis disease only in urban Blantyre Malawi, 2015–2018

Figure 7

Table A1. Characteristics for notified Tuberculosis cases in Urban Blantyre Malawi, 2015–2018. Comparing participants who were included in the study versus those who were not part of the study but were from Blantyre.