# Risk score for predicting mortality including urine lipoarabinomannan detection in hospital inpatients with HIV-associated tuberculosis in sub-Saharan Africa: Derivation and external validation cohort study

Ankur Gupta-Wright, Elizabeth L Corbett, Douglas Wilson, Joep van Oosterhout, Keertan Dheda, Helena Huerga, Jonny Peter, Maryline Bonnet, Melanie Alufandika-Moyo, Daniel Grint, et al.

► To cite this version:

Ankur Gupta-Wright, Elizabeth L Corbett, Douglas Wilson, Joep van Oosterhout, Keertan Dheda, et al.. Risk score for predicting mortality including urine lipoarabinomannan detection in hospital inpatients with HIV-associated tuberculosis in sub-Saharan Africa: Derivation and external validation cohort study. PLoS Medicine, 2019, 16 (4), pp.e1002776. 10.1371/journal.pmed.1002776. hal-02474761

HAL Id: hal-02474761

https://hal.science/hal-02474761

Submitted on 11 Feb 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Risk score for predicting mortality including urine lipoarabinomannan detection in hospital inpatients with HIV-associated tuberculosis in sub-Saharan Africa: Derivation and external validation cohort study

Ankur Gupta-Wright<sup>1,2,3*</sup>, Elizabeth L. Corbett<sup>1,2,3</sup>, Douglas Wilson<sup>4</sup>, Joep J. van Oosterhout<sup>5,6</sup>, Keertan Dheda<sup>7,8</sup>, Helena Huerga<sup>9</sup>, Jonny Peter<sup>10</sup>, Maryline Bonnet<sup>9,11</sup>, Melanie Alufandika-Moyo<sup>3,5</sup>, Daniel Grint<sup>1,12</sup>, Stephen D. Lawn<sup>1,2†</sup>, Katherine Fielding<sup>1,12,13</sup>

1 TB Centre, London School of Hygiene & Tropical Medicine, London, United Kingdom, 2 Clinical Research Department, London School of Hygiene & Tropical Medicine, London, United Kingdom, 3 Malawi–Liverpool–Wellcome Trust Clinical Research Programme, College of Medicine, University of Malawi, Blantyre, Malawi, 4 Department of Medicine, Edendale Hospital, University of KwaZulu-Natal, Pietermaritzburg, South Africa, 5 Dignitas International, Zomba, Malawi, 6 Department of Medicine, College of Medicine, University of Malawi, Blantyre, Malawi, 7 Centre for Lung Infection and Immunity, University of Cape Town, Cape Town, South Africa, 8 Institute of Infectious Diseases and Molecular Medicine, University of Cape Town, Cape Town, South Africa, 9 Epicentre, Paris, France, 10 Division of Allergology and Clinical Immunology, Department of Medicine, University of Cape Town, Cape Town, South Africa, 11 Institute of Research for Development (IRD), UMI 233 TransVIHMI–UM–INSERM U1175, Montpellier, France, 12 Department of Infectious Disease Epidemiology, London School of Hygiene & Tropical Medicine, London, United Kingdom, 13 University of the Witwatersrand, Johannesburg, South Africa

† Deceased.

* ankur.gupta-wright@lshtm.ac.uk

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Gupta-Wright A, Corbett EL, Wilson D, van Oosterhout JJ, Dheda K, Huerga H, et al. (2019) Risk score for predicting mortality including urine lipoarabinomannan detection in hospital inpatients with HIV-associated tuberculosis in sub-Saharan Africa: Derivation and external validation cohort study. PLoS Med 16(4): e1002776. https://doi.org/10.1371/journal.pmed.1002776

**Academic Editor:** Mark Hatherill, University of Cape Town, SOUTH AFRICA

**Received:** October 28, 2018

**Accepted:** March 6, 2019

**Published:** April 5, 2019

**Copyright:** © 2019 Gupta-Wright et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Data are publicly available through London School of Hygiene & Tropical Medicine Data Compass repository. https://doi.org/10.17037/DATA.00001026

**Funding:** The funding for the STAMP trial was provided by the Joint Global Health Trials Scheme of the UK Department of Health and Social Care, the Department for International Development, the Global Challenges Research Fund, the Medical

## Abstract

### Background

The prevalence of and mortality from HIV-associated tuberculosis (HIV/TB) in hospital inpatients in Africa remains unacceptably high. Currently, there is a lack of tools to identify those at high risk of early mortality who may benefit from adjunctive interventions. We therefore aimed to develop and validate a simple clinical risk score to predict mortality in high-burden, low-resource settings.

### Methods and findings

A cohort of HIV-positive adults with laboratory-confirmed TB from the STAMP TB screening trial (Malawi and South Africa) was used to derive a clinical risk score using multivariable predictive modelling, considering factors at hospital admission (including urine lipoarabinomannan [LAM] detection) thought to be associated with 2-month mortality. Performance was evaluated internally and then externally validated using independent cohorts from 2 other studies (LAM-RCT and a Médecins Sans Frontières [MSF] cohort) from South Africa, Zambia, Zimbabwe, Tanzania, and Kenya. The derivation cohort included 315 patients
Research Council and Wellcome Trust (MR/
M007375/1). AG-W has received a Royal College of
Physicians London JMGP Fellowship, ELC has
received a Wellcome Trust Fellowship (WT200901/
Z/16/Z). Funding for LAM-RCT was provided by the
European Developing Clinical Trials Partnership the
South African Medical Research Council, and the
South African National Research Foundation for the
TB Neat project. The MSF cohort study at Homa
Bay County Hospital, Kenya was funded by
Medecins Sans Frontières. The funders had no role
in study design, data collection and analysis,
decision to publish, or preparation of the
manuscript.

**Competing interests:** I have read the journal's policy and the authors of this manuscript have the following competing interests: KD has obtained speaker fees at industry-sponsored symposia and non-financial support from Alere in the form of kits and test strips, outside the submitted work. No other authors declare competing interests.

**Abbreviations:** ART, antiretroviral therapy; C-index, concordance index; HIV/TB, HIV-associated tuberculosis; LAM, lipoarabinomannan; MSF, Médecins Sans Frontières; OR, odds ratio; TB, tuberculosis; TB-LAM, Determine TB LAM Ag; Xpert, Xpert MTB/RIF.

enrolled from October 2015 and September 2017. Their median age was 36 years (IQR 30–43), 45.4% were female, median CD4 cell count at admission was 76 cells/µl (IQR 23–206), and 80.2% (210/262) of those who knew they were HIV-positive at hospital admission were taking antiretroviral therapy (ART). Two-month mortality was 30% (94/315), and mortality was associated with the following factors included in the score: age 55 years or older, male sex, being ART experienced, having severe anaemia (haemoglobin < 80 g/l), being unable to walk unaided, and having a positive urinary Determine TB LAM Ag test (Alere). The score identified patients with a 46.4% (95% CI 37.8%–55.2%) mortality risk in the high-risk group compared to 12.5% (95% CI 5.7%–25.4%) in the low-risk group (*p* < 0.001). The odds ratio (OR) for mortality was 6.1 (95% CI 2.4–15.2) in high-risk patients compared to low-risk patients (*p* < 0.001). Discrimination (c-statistic 0.70, 95% CI 0.63–0.76) and calibration (Hosmer-Lemeshow statistic, *p* = 0.78) were good in the derivation cohort, and similar in the external validation cohort (complete cases *n* = 372, c-statistic 0.68 [95% CI 0.61–0.74]). The validation cohort included 644 patients between January 2013 and August 2015. Median age was 36 years, 48.9% were female, and median CD4 count at admission was 61 (IQR 21–145). OR for mortality was 5.3 (95% CI 2.2–9.5) for high compared to low-risk patients (complete cases *n* = 372, *p* < 0.001). The score also predicted patients at higher risk of death both pre- and post-discharge. A simplified score (any 3 or more of the predictors) performed equally well. The main limitations of the scores were their imperfect accuracy, the need for access to urine LAM testing, modest study size, and not measuring all potential predictors of mortality (e.g., tuberculosis drug resistance).

Conclusions

This risk score is capable of identifying patients who could benefit from enhanced clinical
care, follow-up, and/or adjunctive interventions, although further prospective validation stud-
ies are necessary. Given the scale of HIV/TB morbidity and mortality in African hospitals,
better prognostic tools along with interventions could contribute towards global targets to
reduce tuberculosis mortality.

Author summary

Why was this study done?

*   HIV-associated tuberculosis (TB) is very common in hospitals in sub-Saharan Africa, and is a major cause of morbidity and mortality.

*   There is a lack of tools to identify which patients are more likely to die early; therefore, these patients cannot be targeted for more intensive clinical care or other treatments in addition to TB antibiotics and antiretroviral drugs.

*   A new urine TB diagnostic test (detecting a substance called lipoarabinomannan [LAM]) can identify patients at higher risk of dying and, along with other simply measured clinical signs or symptoms, may be useful for predicting which patients are most likely to do poorly.
## What did the researchers do and find?

* We used data from HIV-positive adults admitted to hospital in Malawi and South Africa and diagnosed with TB to develop a simple clinical risk score to identify patients with a 50% chance of dying within 2 months.

* The score included 6 factors that can be measured at admission to hospital, including the results of the urine TB test, presence of anaemia, and some demographic factors.

* We then tested the performance of the risk score using data from studies of different patients from sub-Saharan Africa, and it performed similarly.

* The score was able to identify patients at higher risk of dying during admission to hospital, and after discharge from hospital.

## What do these findings mean?

* This score could be used to identify patients admitted to hospital and diagnosed with HIV-associated TB who may benefit most from more intensive clinical care, additional treatments, and/or closer follow-up after discharge.

* It could also be used as a research tool to study new drugs or strategies aimed at reducing mortality from with HIV-associated TB.

* This is the first such tool to our knowledge in this patient population, and future studies could optimise such predictive tools, in particular if appropriate new interventions and/or diagnostics become available.

## Introduction

Tuberculosis (TB) is the leading infectious disease killer globally, causing an estimated 1.7 million deaths globally in 2017 [1]. This burden lies disproportionately in people living with HIV, who account for approximately 1 in 4 TB deaths. The case fatality rate of HIV-associated TB (HIV/TB) is particularly high in hospitals, estimated at 29% in a recent meta-analysis [2]. This may be an underestimate, given that post-mortem studies from sub-Saharan Africa have demonstrated that a high proportion of HIV-positive deaths in facilities have evidence of undiagnosed TB [3].

Interventional studies aiming to reduce mortality in this patient population have demonstrated mortality reductions with improved TB diagnostics [4,5] and appropriately timed initiation of antiretroviral therapy (ART) [6,7]. However, mortality remains substantial despite these interventions, and adjunctive interventions are likely to be needed to further impact mortality. Currently, predictors for mortality are poorly defined. Being able to identify patients at the highest risk of mortality could inform the development and assessment of new interventions, and also identify which patients would benefit most from interventions beyond TB therapy and appropriately timed ART [8].

Clinical decision tools and risk scores are used widely in clinical practice to simplify the identification of patients at highest risk for poor health outcomes. Predictor scores for
mortality have been developed for HIV-associated cryptococcal meningitis and pneumonia, and are used to guide management in *Pneumocystis jiroveci* pneumonia [9–11]. Although scores have been developed to predict risk of TB disease in various populations, including TB bacteraemia in hospitalised patients [12], to our knowledge no externally validated scores exist to predict outcomes of TB disease among hospitalised patients with HIV [13,14]. Scores developed to predict TB mortality in settings with low HIV prevalence are also of limited use in people living with HIV due to differences in clinical presentation, pathogenesis, and outcomes [15–18]. A recent study from the US developed and internally validated a score to predict mortality in HIV/TB in low-prevalence settings (US), but this would not be applicable to hospitalised patients in Africa given that many of the variables are not routinely available [19].

We have previously shown that detection of lipoarabinomannan (LAM) in the urine of HIV/TB patients using a cheap (approximately US$3) and quick (testing takes 25 minutes) lateral flow assay is independently associated with a 2- to 3-fold increased risk of mortality [20]. We therefore aimed to investigate if urinary LAM detection, along with other clinical variables readily available in high-burden settings, could be used to predict which HIV-positive patients admitted to hospital and diagnosed with TB were at high risk of early mortality, and to externally validate the predictive tool.

## Methods

### Study design and participants for prediction tool development

We used data from the STAMP ('rapid urine-based screening for tuberculosis in HIV-positive patients admitted to hospital in Africa') trial for the clinical risk score derivation [5,21]. The STAMP trial recruited HIV-positive adults (aged 18 years or more), irrespective of symptoms or clinical presentation, who were admitted to medical wards of 2 hospitals in Malawi and South Africa between 26 October 2015 and 19 September 2017. On admission, patients were screened for TB using Xpert MTB/RIF (Xpert; Cepheid) on sputum in both study arms, and Xpert and Determine TB LAM Ag (TB-LAM; Alere) assays on urine in the intervention arm.

Exclusion criteria in the trial were already taking TB treatment and inability to give consent. The clinical teams managing the patients were masked to which TB tests were positive; therefore, management of TB patients should not have differed between arms. The management of HIV/TB in the study hospitals was representative of their local settings and followed local and national guidelines, with no input from the study team (beyond TB diagnostic tests).

Patients diagnosed with TB in the standard-of-care arm had stored urine tested with Xpert and TB-LAM retrospectively. Data were collected at baseline (at or close to admission) on demographics and clinical characteristics, and subsequently on TB investigations and treatment, and clinical events, including death or discharge from hospital. Patients discharged alive were followed up at 2 months by outpatient attendance, home visit, or telephone for vital status. The derivation cohort included all patients (from both trial arms) with laboratory-confirmed TB. The outcome was mortality risk at 2 months after admission. Patients lost to follow-up were assumed alive at 56 days.

### Definitions

Laboratory-confirmed TB was defined as any 1 of a positive smear microscopy, mycobacterial culture, Xpert from any site, or urinary TB-LAM. TB-LAM assay was positive if recorded as 'grade 1' or higher on the manufacturer's (post-2014) reference card. Ability to walk unaided was assessed by healthcare workers (not self-reported by patients), and was equivalent to a Karnofsky functional score below 40 points [22]. WHO danger signs were heart rate > 120 beats per minute, respiratory rate > 30 per minute, temperature > 39°C, and being unable to
walk unaided. 'ART experienced' was defined as receiving ART at the time of enrolment to the study.

## Score derivation

Candidate predictor variables were identified for inclusion in the predictive model based on a priori clinical knowledge, previous literature, and the need for variables to be objective, reproducible, and available in resource-constrained settings [23]. We considered variables known to be associated with mortality in HIV/TB, including age, sex, ART experience, physiological measurements at admission, weight and/or body mass index, CD4 cell count, functional status (being unable to walk unaided), and haemoglobin [24–28]. Time on ART was not considered as not all patients take ART, and because of challenges in accurately ascertaining duration. Where 2 or more predictors were highly correlated, only 1 was selected, to simplify the prognostic model, as inclusion of all would contribute little additional predictive information [23]. Analyses were planned prospectively (see S1 Appendix) except where indicated as post hoc.

Continuous variables were assessed for non-linearity using fractional polynomials, and categorised based on previously established cutoffs (e.g., CD4 cell count and haemoglobin) or associations with mortality (e.g., age and weight, using the *fp plot* command in Stata). Complete case analysis was chosen for the derivation score as few data (<5%) were missing. We first performed univariable analyses assessing the association of each variable with mortality risk using logistic regression. We then used a backward elimination, stepwise approach to create a multivariable predictive model, starting with all candidate variables, and excluding variables sequentially if *p* > 0.1 using likelihood ratio tests and the Akaike information criterion. Given that there were 94 deaths, we did not want to estimate more than 9 candidate predictors (various studies have shown each candidate predictor studied requires a minimum of 10 events) [29]. Interactions were also assessed using likelihood ratio testing. All analyses were done using Stata version 14, and all *p*-values were 2-sided.

Regression coefficients from the final multivariable model were multiplied by the smallest possible constant and then rounded to the nearest integer, and then assigned as 'points' to each variable. The clinical risk score was derived by combining the points based on each patient's characteristics. High-, medium-, and low-risk groups for mortality were then arbitrarily defined after plotting risk score against observed mortality such that the high-risk group accounted for most (>50%) deaths and the low-risk group accounted for as few deaths as possible.

## Risk score evaluation and internal performance

Mortality risk at 2 months and 95% confidence intervals (CIs) were calculated for each risk group, as were odds ratios (ORs) and 95% CI for mortality. In exploratory analyses, inpatient and outpatient (post-discharge) deaths were also compared between risk groups by restricting analyses to deaths occurring during hospital admission or to deaths occurring after discharge in the subset of patients who were discharged alive from hospital. CD4 cell count and TB-LAM grade were also compared between risk groups. Mortality risk was compared between groups using chi-squared tests.

We assessed the model discrimination (ability to differentiate patients who would die within 2 months and those who would survive) by calculating the concordance index (c-statistic) (also known as the area under the receiver operator curve), assuming a c-statistic < 0.6 showed poor discrimination [30]. Model calibration was assessed by plotting the probability of mortality predicted by the model against observed mortality in the derivation dataset using a calibration plot and the Hosmer-Lemeshow test, assuming a *p* < 0.05 indicated poor
calibration. In post hoc analysis, in response to reviewer request and to better understand the
utility of the score, the sensitivity and positive predictive value of the score were calculated.

## External validation

To externally validate the clinical risk score, we used data collected independently from 2 stud-
ies: (1) a multicentre diagnostic clinical trial of adjunctive urine TB-LAM testing in HIV-posi-
tive patients with TB symptoms who were admitted to hospitals in 4 sub-Saharan African
countries (South Africa, Zambia, Zimbabwe, and Tanzania) (LAM-RCT) [31] and (2) a pro-
spective cohort study assessing the diagnostic yield of TB-LAM in HIV-positive patients with
TB symptoms in Kenya (Médecins Sans Frontières [MSF] cohort study) [32]. Patients were
included in the validation cohort if they were adults and had laboratory-confirmed TB (as pre-
viously defined). Patients from the LAM-RCT in the ‘no TB-LAM’ arm were excluded, as were
outpatients (i.e., patients not admitted to hospital) from the MSF cohort study.

The validation cohort sites were all in settings in sub-Saharan Africa with high HIV prevalence and TB incidence, but differed from the derivation cohort in that all patients had at least 1 TB symptom (cough, fever, weight loss, or night sweats). The LAM-RCT recruitment occurred between 1 January 2013 and 2 October 2014, and the MSF cohort recruitment between 22 October 2013 and 20 August 2015. Mortality outcomes were assessed at 2 months in both studies.

The clinical risk score for mortality was calculated by assigning the same ‘points’ to vari-
ables as for the derivation cohort, and the same cutoffs were used to define high-, medium-,
and low-risk groups for mortality. Patients with missing observations were excluded, for a
complete case analysis. However, sensitivity analyses were done for score performance using
multivariate multiple imputation with chained equations for missing data as 42% of patients
had missing data in the validation cohort. Data were assumed to be missing at random, and
were imputed for missing candidate predictor variables using mortality risk, other candidate
predictor variables, and other baseline demographic variables, with 100 imputations.

Evaluation of the score in the validation dataset was done using the same statistical methods
as the internal evaluation, with calculation of mortality risk at 2 months, ORs for mortality,
and survival curves. Discrimination was assessed using the c-statistic, and calibration with a
calibration plot and the Hosmer-Lemeshow test.

The study is reported in concordance with TRIPOD guidance for multivariable prediction
models (see S2 Appendix) [33]. Ethical approval for each of the source studies was obtained
from the relevant ethics committees in the country of data collection and from the trial sponsors
(see S3 Appendix for list of ethics committees). All patients provided informed written consent.

## Results

### Baseline characteristics

Of 506 HIV-positive patients diagnosed with TB in the STAMP trial derivation cohort, 322 had laboratory-confirmed TB. Seven patients were excluded from the complete case analysis for missing data (Fig 1). The median age of TB patients included in the derivation cohort was 36 years (interquartile range [IQR] 30–43), 172 (55%) were men, 53 (17%) were newly diagnosed with HIV, and the median CD4 cell count was 76 cells/µl (IQR 23–206; Table 1). In all, 209 (65%) patients were positive on urine TB-LAM testing, indicating probable disseminated TB disease. Anaemia was common and median haemoglobin was 86 g/l (IQR 67–108). Patients presented with advanced disease: 133 (42%) had 1 or more WHO danger signs, and 71 (23%) were severely disabled or unable to walk unaided.

In the derivation cohort, 94 (30%) patients died within 2 months, with 66 (70%) dying dur-
ing their hospital admission; 29 (31% of deaths) patients died by 1 week, and 52 (55%) by 2
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Study profile. LAM, lipoarabinomannan; MSF, Médecins Sans Frontières; TB, tuberculosis.</figcaption></figure>

https://doi.org/10.1371/journal.pmed.1002776.g001

weeks, after admission. In unadjusted analyses, mortality risk was higher in patients aged 55 years or older, men, ART-experienced patients, those unable to walk, patients with severe anaemia (haemoglobin < 80 g/l), patients with CD4 cell count < 100 cells/µl, and those with positive urine TB-LAM tests (Table 2). Six out of 322 (2%) patients were lost to follow-up after hospital discharge.

## Multivariable model and clinical risk score

The final multivariable logistic regression model for mortality at 2 months included age, sex, ART experience, haemoglobin, functional status (being unable to walk unaided), and urine TB-LAM result (Table 2). For associations of linear continuous variables with mortality see S1 Table. CD4 count and weight were dropped from the final predictor score model as their relationship with mortality was mediated by functional status and urine TB-LAM result. We found no significant interactions between variables in the final model. The c-statistic for the predictive model in the derivation dataset was 0.70 (95% CI 0.63–0.76), showing moderate discrimination. Calibration of the predictive model was good, as shown by the calibration plot (see S1 Fig) and a Hosmer-Lemeshow statistic *p* = 0.78.

The clinical risk score for mortality, based on the regression coefficients, is outlined in Fig 2. Observed and predicted mortality risks for the risk score are reported in S2 Fig. Mortality risk groups were defined as low risk (10 points or fewer), medium risk (11 to 20 points), or high risk (more than 20 points) (Fig 3). Therefore, in the derivation cohort, 48 (15%) patients were deemed low risk, 142 (45%) were deemed medium risk, and 125 (40%) were deemed high risk. Median risk score was 19 (IQR 13–22, range 0–42). Observed mortality risk by 2 months was 12.5% (95% CI 5.7%–25.4%), 21.1% (95% CI 15.1%–28.7%), and 46.4% (95% CI 37.8%–55.2%) in the low-, medium-, and high-risk groups, respectively (*p* < 0.001). ORs for mortality were 6.1 (95% CI 2.4–15.2) in the high-risk group and 1.9 (95% CI 0.7–4.8) in the medium-risk group compared to low-risk patients (*p* < 0.001).

## Simplified clinical risk score for mortality

As the regression coefficients and points in the clinical risk score were similar for all 6 variables, we created a simplified version of the score by assigning each variable within the score 1
<table><caption>Table 1. Baseline characteristics.</caption><thead><tr><th rowspan="2">Characteristic</th><th rowspan="2">Category</th><th colspan="2">Median (IQR) or N (%)</th></tr><tr><th>Derivation dataset (n = 315)</th><th>Validation dataset (n = 644)</th></tr></thead><tbody><tr><td colspan="4"><strong>Demographics</strong></td></tr><tr><td>Age (years)</td><td></td><td>36 (30-43)</td><td>35 (30-42)</td></tr><tr><td>Sex</td><td>Female</td><td>143 (45.4)</td><td>315 (48.9)</td></tr><tr><td rowspan="6">Country of enrolment</td><td>South Africa</td><td>162 (51.4)</td><td>102 (15.8)</td></tr><tr><td>Tanzania</td><td>—</td><td>70 (10.9)</td></tr><tr><td>Zambia</td><td>—</td><td>151 (23.5)</td></tr><tr><td>Zimbabwe</td><td>—</td><td>139 (21.6)</td></tr><tr><td>Kenya</td><td>—</td><td>182 (28.3)</td></tr><tr><td>Malawi</td><td>153 (48.6)</td><td>—</td></tr><tr><td colspan="4"><strong>HIV history</strong></td></tr><tr><td>New HIV diagnosis</td><td>Yes</td><td>53 (16.8)</td><td>—</td></tr><tr><td>Currently taking ART</td><td>Yes</td><td>210 (80.2)</td><td>290 (45.0)</td></tr><tr><td>Time on ART (years)<sup>a</sup></td><td></td><td>1.0 (0.2-4.4)</td><td>0.7 (0.1-3.2)</td></tr><tr><td>CD4 cell count (cells/µl)<sup>b</sup></td><td></td><td>76 (23-206)</td><td>61 (21-145)</td></tr><tr><td colspan="4"><strong>TB history</strong></td></tr><tr><td>Cough</td><td>Yes</td><td>228 (72.4)</td><td>601 (93.5)</td></tr><tr><td>Fever</td><td>Yes</td><td>223 (70.8)</td><td>562 (87.3)</td></tr><tr><td>Weight loss</td><td>Yes</td><td>286 (90.8)</td><td>595 (95.8)</td></tr><tr><td>Night sweats</td><td>Yes</td><td>165 (52.4)</td><td>531 (82.6)</td></tr><tr><td>WHO TB symptom screen</td><td>Yes</td><td>310 (98.4)</td><td>644 (100)</td></tr><tr><td>Previous history of TB</td><td>Yes</td><td>72 (22.9)</td><td>123 (19.1)</td></tr><tr><td colspan="4"><strong>Clinical presentation</strong></td></tr><tr><td>Weight (kg)<sup>c</sup></td><td></td><td>50 (42-57)</td><td>49 (43-55)</td></tr><tr><td>BMI<sup>d</sup></td><td></td><td>19.1 (16.2-21.0)</td><td>17.6 (15.9-20.3)</td></tr><tr><td>Heart rate (bpm)</td><td></td><td>104 (90-118)</td><td>102 (90-119)</td></tr><tr><td>Respiratory rate (per minute)</td><td></td><td>22 (20-26)</td><td>24 (22-28)</td></tr><tr><td>Systolic blood pressure (mm Hg)</td><td></td><td>102 (92-116)</td><td>104 (95-116)</td></tr><tr><td>Temperature (°C)</td><td></td><td>36.5 (36.1-37.2)</td><td>37.0 (36.6-38.0)</td></tr><tr><td>Haemoglobin (g/l)<sup>e</sup></td><td></td><td>86 (67-108)</td><td>85 (68-100)</td></tr><tr><td>WHO danger sign<sup>f</sup></td><td>Yes</td><td>133 (42.2)</td><td>399 (62.0)</td></tr><tr><td>Unable to walk unaided<sup>g</sup></td><td>Yes</td><td>71 (22.5)</td><td>262 (40.7)</td></tr><tr><td colspan="4"><strong>TB diagnosis</strong></td></tr><tr><td>Sputum Xpert positive</td><td>Yes</td><td>168 (52.2)</td><td>217 (33.7)</td></tr><tr><td>Sputum smear microscopy positive</td><td>Yes</td><td>—</td><td>211 (32.8)</td></tr><tr><td>TB culture positive (any site)</td><td>Yes</td><td>—</td><td>388 (60.3)</td></tr><tr><td>Urine LAM positive</td><td>Yes</td><td>209 (64.9)</td><td>424 (65.8)</td></tr><tr><td>Chest X-ray suggestive of TB</td><td>Yes</td><td>107 (33.2)</td><td>336 (52.2)</td></tr><tr><td colspan="4"><strong>Outcome</strong></td></tr><tr><td>Died by 2 months</td><td>Yes</td><td>94 (29.8)</td><td>147 (22.8)</td></tr></tbody></table>

Sputum smear and TB culture were not routinely performed in the STAMP trial (derivation dataset). Missing data are for the validation dataset only.

ART, antiretroviral therapy; BMI, body mass index; bpm, beats per minute; IQR, interquartile range; LAM, lipoarabinomannan; TB, tuberculosis; WHO, World Health Organization; Xpert, Xpert MTB/RIF.

<footnote>aTime on ART missing for 19 (3%) patients.</footnote>

<footnote>bCD4 cell count missing for 27 (4%) patients.</footnote>

<footnote>cWeight missing for 75 (12%) patients.</footnote>

<footnote>dBMI missing for 90 (14%) patients.</footnote>

<footnote>eHaemoglobin missing for 272 (42%) patients.</footnote>

<footnote>fOne of heart rate > 120 bpm, respiratory rate > 30 per minute, temperature > 39°C, or unable to walk unaided.</footnote>

<footnote>gAbility to walk unaided was assessed by healthcare worker and not self-reported.</footnote>
<table><caption>Table 2. Univariable and multivariable logistic regression analysis of factors associated with mortality in the derivation cohort (n = 315).</caption><thead><tr><th rowspan="2">Characteristic</th><th rowspan="2">Category</th><th rowspan="2">Died<sup>a</sup> (n = 94)</th><th colspan="2">Univariable</th><th colspan="2">Multivariable</th><th rowspan="2">Regression (β) coefficient</th></tr><tr><th>OR (95% CI)</th><th>p-Value</th><th>OR (95% CI)</th><th>p-Value</th></tr></thead><tbody><tr><td colspan="8"><strong>Demographics</strong></td></tr><tr><td rowspan="2">Age</td><td>&lt;55 years</td><td>82 (28.4)</td><td>1 (ref)</td><td rowspan="2">0.067</td><td>1 (ref)</td><td rowspan="2">0.10</td><td rowspan="2">0.710</td></tr><tr><td>≥55 years</td><td>12 (46.2)</td><td>2.2 (1.0-4.9)</td><td>2.0 (0.9-4.9)</td></tr><tr><td rowspan="2">Sex</td><td>Female</td><td>32 (22.4)</td><td>1 (ref)</td><td rowspan="2">0.012</td><td>1 (ref)</td><td rowspan="2">0.001</td><td rowspan="2">0.923</td></tr><tr><td>Male</td><td>62 (36.0)</td><td>2.0 (1.2-3.2)</td><td>2.5 (1.5-4.3)</td></tr><tr><td colspan="8"><strong>HIV infection</strong></td></tr><tr><td rowspan="2">ART experienced</td><td>No</td><td>18 (20.9)</td><td>1 (ref)</td><td rowspan="2">0.024</td><td>1 (ref)</td><td rowspan="2">0.048</td><td rowspan="2">0.621</td></tr><tr><td>Yes</td><td>76 (33.2)</td><td>1.9 (1.1-3.4)</td><td>1.9 (1.0-3.5)</td></tr><tr><td rowspan="2">CD4 cell count<sup>#</sup></td><td>≥100 cells/μl</td><td>31 (23.9)</td><td>1 (ref)</td><td rowspan="2">0.040</td><td>—</td><td rowspan="2">—</td><td rowspan="2">—</td></tr><tr><td>&lt;100 cells/μl</td><td>62 (33.9)</td><td>1.7 (1.0-2.8)</td><td>—</td></tr><tr><td colspan="8"><strong>Clinical presentation</strong></td></tr><tr><td rowspan="2">WHO danger sign</td><td>No</td><td>50 (27.5)</td><td>1 (ref)</td><td rowspan="2">0.185</td><td>—</td><td rowspan="2">—</td><td rowspan="2">—</td></tr><tr><td>Yes</td><td>48 (34.8)</td><td>1.4 (0.9-2.2)</td><td>—</td></tr><tr><td rowspan="4">Weight</td><td>&lt;35 kg</td><td>10 (43.5)</td><td>3.3 (1.2-8.9)</td><td rowspan="4">0.054</td><td>—</td><td rowspan="4">—</td><td rowspan="4">—</td></tr><tr><td>35-60 kg</td><td>73 (31.5)</td><td>1.8 (0.9-3.5)</td><td>—</td></tr><tr><td>&gt;60 kg</td><td>11 (18.3)</td><td>1 (ref)</td><td>—</td></tr><tr><td></td><td></td><td></td><td>—</td><td>—</td></tr><tr><td rowspan="2">Haemoglobin (g/l)</td><td>≥80 g/l</td><td>44 (23.7)</td><td>1 (ref)</td><td rowspan="2">0.003</td><td>1 (ref)</td><td rowspan="2">0.008</td><td rowspan="2">0.703</td></tr><tr><td>&lt;80 g/l</td><td>50 (38.8)</td><td>2.0 (1.3-3.3)</td><td>2.0 (1.2-3.4)</td></tr><tr><td rowspan="2">Unable to walk unaided</td><td>No</td><td>64 (26.2)</td><td>1 (ref)</td><td rowspan="2">0.004</td><td>1 (ref)</td><td rowspan="2">0.022</td><td rowspan="2">0.689</td></tr><tr><td>Yes</td><td>30 (42.3)</td><td>2.2 (1.3-3.8)</td><td>2.0 (1.1-3.6)</td></tr><tr><td colspan="8"><strong>TB diagnosis</strong></td></tr><tr><td rowspan="2">Urine LAM positive</td><td>No</td><td>24 (22.6)</td><td>1 (ref)</td><td rowspan="2">0.044</td><td>1 (ref)</td><td rowspan="2">0.040</td><td rowspan="2">0.603</td></tr><tr><td>Yes</td><td>70 (33.5)</td><td>1.7 (1.0-2.9)</td><td>1.8 (1.0-3.2)</td></tr></tbody></table>

The constant (intercept) was -2.8. *p*-Values were calculated by likelihood ratio tests. There was no evidence of interaction between urine LAM positivity, being unable to walk, and haemoglobin < 80 g/l in the multivariable model (likelihood ratio test *p*-values all >0.1). Weight and being unable to walk were strongly associated.

point if present (age 55 years or over, male sex, ART experienced, severe anaemia, being unable to walk unaided, or urine TB-LAM positive; see S3 Fig). A high mortality risk was defined as 3 or more points, medium risk as 2 points, and low risk as 0 or 1 point.

In the derivation cohort, patients with 3 or more points (high risk) had a mortality of 40.0% (70/175, 95% CI 33.0%–47.5%), compared to 19.6% (18/92, 95% CI 12.6%–29.0%) mortality in those with 2 points and 12.5% (6/48, 95% CI 5.7%–25.4%) in those with 0 or 1 point (*p* < 0.001) (see Fig 3 and S2 Table). The sensitivity of the risk score for mortality in was 0.75 (the score correctly identified 70/94 deaths), and the positive predictive value was 0.4.

The clinical risk score was useful in predicting deaths that occurred during inpatient admission (50 [28.6%, 95% CI 22.3%–35.8%] in the high-risk group compared to 6 [10.9%, 95% CI 5.9%–19.1%] in the low-risk group, *p* = 0.001) as well as deaths occurring after discharge (20 [16.0%, 95% CI 10.5%–23.6%] in the high-risk group compared to 0 [0%] in the low-risk group, *p* = 0.015; Fig 4). More patients in the high-risk group were TB-LAM positive and had higher grades of positive result, but CD4 cell count did not differ by risk group (see S4 Fig). Survival curves by risk group are presented in Fig 5.

<footnote><sup>a</sup>Data are number of patients in category who died (%).</footnote>

<footnote>ART, antiretroviral therapy; LAM, lipoarabinomannan; OR, odds ratio; TB, tuberculosis.</footnote>
**Demographic factors:**

1. Is the patient male? Yes: add 9 points

2. Is the patient aged 55 years or older? Yes: add 7 points

**HIV factors:**

3. Is the patient currently taking antiretroviral therapy? Yes: add 6 points

**Clinical presentation and TB diagnosis:**

4. Is the patient unable to walk unaided? Yes: add 7 points

5. Does the patient have severe anaemia (haemoglobin <80g/l)? Yes: add 7 points

6. Is the patient positive on urine TB-LAM testing? Yes: add 6 points

**Total points (min 0, max 42):**

**Fig 2. Risk score calculation to predict mortality.** TB, tuberculosis; TB-LAM, Determine TB LAM Ag.

https://doi.org/10.1371/journal.pmed.1002776.g002

**External validation**

The external validation cohort included 644 HIV-positive patients with laboratory-confirmed TB, of whom 372 (58%) patients had no missing data for the risk score and were therefore included in the complete case analysis (Fig 1). Baseline characteristics were similar between cohorts, although fewer patients reported taking ART and more patients presented with severe functional impairment and 1 or more WHO danger signs in the validation cohort (Table 1). A similar proportion of patients were positive on urine TB-LAM testing (65% in the derivation cohort compared to 66% in the validation cohort). Mortality at 2 months was lower in the validation cohort (22.8%) compared to the derivation cohort (29.8%). Loss to follow-up was 4% in the validation cohort (15/372).

In complete case analysis (n = 372), the observed mortality risks in the validation cohort were 8.7% (95% CI 4.6%–16.0%) in the low-risk group, 19.1% (95% CI 13.2%–26.8%) in the medium-risk group, and 35.5% (95% CI 27.9%–43.9%) in the high-risk group (see Fig 4). Median risk score was 16 (IQR 10–22, range 0–42). The ORs for mortality by risk group were similar to those in the derivation cohort (5.8 [95% CI 2.7–12.3] for the high-risk group and 2.5 [95% CI 1.1–5.5] for the medium-risk group compared to the low-risk group). The risk score was also useful in predicting both inpatient and post-discharge deaths (the high-risk group had a 20% risk of post-discharge death compared to 5% in the low-risk group). The simplified risk score performed similarly to the full score in the validation cohort.

The predictive model had similar calibration and discrimination in the validation cohort as in the derivation cohort: the c-statistic was 0.68 (95% CI 0.61–0.74; see S3 Table), and the Hosmer-Lemeshow statistic had *p* = 0.13 (see S5 Fig for the calibration plot). In a sensitivity analysis using multiple imputation for missing data in the validation dataset (n = 644), the c-statistic for the predictive model was 0.64 (95% CI 0.60–0.69), and the Hosmer-Lemeshow statistic had *p* = 0.67. ORs for mortality were 5.3 (95% CI 2.2–9.5) for the high-risk group and 2.1 (95% CI 1.0–4.6) for the medium-risk group compared to low-risk patients (*p* < 0.001).

**Discussion**

In this study, we developed and externally validated a pragmatic clinical risk score to predict early mortality in HIV-positive patients admitted to hospital and diagnosed with laboratory-confirmed TB. Our score used 6 clinical and laboratory factors that could be readily collected at admission to hospital in settings with high HIV and TB burden. The score was able to
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig 3.** Distribution of risk scores and mortality in the derivation dataset. Distribution of risk scores for mortality stratified by outcome at 2 months (stacked bar chart) and mortality risk (percent, shown by blue line) for (A) the full risk score (based on the regression coefficients) and (B) the simplified risk score. Mortality risks and absolute numbers in each category are presented in S2 Table.</figcaption></figure>

https://doi.org/10.1371/journal.pmed.1002776.g003

categorise patients into 3 risk groups. One-third of the high-risk group died during hospital admission, and almost 50% had died by 2 months. A simplified ‘3 of 6 predictors’ version of the score performed similarly. This is the first study to our knowledge to derive and externally validate a risk score to predict mortality in this patient population.

We found older age, being male, being ART experienced, having severe anaemia, being severely disabled or unable to walk unaided, and being urine TB-LAM positive were all risk
A)

<figure><img src="image_7.png" /></figure>

B)

<figure><img src="image_8.png" /></figure>

C)

<figure><img src="image_9.png" /></figure>
**Fig 4.** Observed mortality risk by risk score category in the derivation and validation cohorts. Observed mortality risk (A) at 56 days, (B) during inpatient stay, and (C) post-discharge in the derivation and validation cohorts, stratified by risk score category (derivation and validation cohorts) and simplified risk score category (derivation cohort). Numbers on bars represent absolute mortality risk; error bars represent 95% confidence intervals. For the full risk score, low risk was defined as 10 points or fewer, medium risk as 11 to 20 points, and high risk as more than 20 points. For the simplified risk score, the low-risk group had a predictor score of 0 or 1 point, the medium-risk group had a predictor score of 2 points, and the high-risk group had predictor score of ≥3 points. *p*-Values based on chi-squared tests between groups for derivation and validation cohorts, respectively, are (A) *p* < 0.001 and *p* < 0.001, (B) *p* = 0.001 and *p* = 0.001, and (C) *p* = 0.015 and *p* = 0.003.

https://doi.org/10.1371/journal.pmed.1002776.g004

factors for mortality. These factors have been established as being associated with outcome in HIV/TB in previous studies [24–28], and most likely reflect more advanced HIV-related immunosuppression, late presentation to healthcare services, and/or poorer underlying physiological reserve. Positive urine diagnostic tests (including LAM detection and *Mycobacterium tuberculosis* nucleic acid detection) in the context of HIV infection are thought to represent haematogenously disseminated renal TB with high mycobacterial burden, which may explain why it is associated with a worse prognosis [34]. Interestingly, clinical signs and symptoms (such as WHO danger signs) were not predictive of mortality.

In contrast to previously published data, ART-experienced patients had a higher mortality risk in our study [28,35]. This likely reflects a high burden of unrecognised ART failure, due to either poor adherence or drug resistance among patients admitted to hospital. Another potential cause is immune reconstitution inflammatory syndrome (IRIS) in patients who have recently started ART. The relationship between ART and mortality is likely to be more complex, representing different groups of patients with different mortality risks, but for this pragmatic tool we have not been able to explore this further. CD4 cell count, which has been previously shown to be associated with mortality in HIV-positive patients, dropped out of our final multivariable predictive model due to mediation by other variables. Furthermore, in the era of test and treat for HIV and use of quantitative HIV viral load for monitoring, CD4 testing services are being scaled back, and are often not available in resource-limited settings.

Mechanisms and causes of mortality in advanced HIV/TB are still not well understood. Co-pathologies, including other opportunistic infections and bacterial pneumonia or sepsis, are commonly detected post-mortem [36,37]. High-risk patients could be prioritised for screening for co-infections, for example using cryptococcal antigen point-of-care tests, or empirical prophylactic treatment with antibacterial agents, an approach that has been shown to reduce mortality in advanced HIV infection [38].

Whilst this clinical risk score can identify patients with the highest risk of mortality, there remains an absence of proven interventions (beyond TB therapy and appropriately timed ART) to reduce mortality in this population. Therefore, we propose this score could be used as a clinical tool to alert clinicians to patients at high risk of mortality who should be reviewed before discharge and/or flagged for early clinical follow-up in settings where urine TB-LAM scale-up is occurring. The score could also be used as a research tool to aid evaluation of intensified or optimised TB treatment regimens or adjunctive interventions aimed at reducing high mortality in this population.

Possible interventions include rapid viral load testing with ART adherence support and early switching for those with virological failure. Host-directed therapies, which target host immune responses, are in clinical trials for TB, including some specifically for HIV/TB [39,40]. Patients identified as being at highest risk for mortality could also be offered more intensive monitoring or supportive care, for example better management of severe anaemia [41], although optimal strategies of supportive care are not clear [42]. Enhanced treatment and prophylaxis for co-infections have been shown to reduce early mortality in patients with
<figure><img src="image_10.png" /><figcaption>(A)</figcaption></figure>

<table><caption>(B)</caption><thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr><tr><th>Number at risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th></tr></thead><tbody><tr><td>49</td><td>46</td><td>45</td><td>42</td><td>41</td><td>41</td><td>41</td><td>41</td><td>41</td><td>41</td><td>41</td><td>41</td><td>41</td><td>41</td><td>41</td></tr><tr><td>92</td><td>87</td><td>83</td><td>80</td><td>78</td><td>78</td><td>77</td><td>76</td><td>76</td><td>73</td><td>73</td><td>76</td><td>76</td><td>73</td><td>73</td></tr><tr><td>181</td><td>154</td><td>133</td><td>125</td><td>119</td><td>113</td><td>109</td><td>105</td><td>105</td><td>103</td><td>103</td><td>105</td><td>105</td><td>103</td><td>103</td></tr></tbody></table>

<figure><img src="image_11.png" /></figure>

<table><thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr><tr><th>Number at risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th><th>Low risk</th><th>Medium risk</th><th>High risk</th></tr></thead><tbody><tr><td>103</td><td>102</td><td>97</td><td>96</td><td>93</td><td>93</td><td>93</td><td>93</td><td>93</td><td>91</td><td>91</td><td>93</td><td>91</td><td>91</td><td>91</td></tr><tr><td>131</td><td>122</td><td>115</td><td>110</td><td>107</td><td>106</td><td>103</td><td>101</td><td>101</td><td>100</td><td>100</td><td>101</td><td>101</td><td>100</td><td>100</td></tr><tr><td>138</td><td>117</td><td>104</td><td>100</td><td>96</td><td>91</td><td>86</td><td>84</td><td>84</td><td>81</td><td>81</td><td>84</td><td>84</td><td>81</td><td>81</td></tr></tbody></table>
**Fig 5.** Survival curves stratified by clinical risk score category in the derivation dataset. Survival curves and risk tables with number at risk for the (A) derivation cohort and (B) validation cohort stratified by risk group using the simplified clinical risk score. The low-risk group (blue line) had a risk score of 0 or 1 point, the medium-risk group (red line) had a risk score of 2 points, and the high-risk group (green line) had a risk score of ≥3 points. Log-rank test *p* < 0.001.

https://doi.org/10.1371/journal.pmed.1002776.g005

advanced HIV initiating ART [38], and may also benefit those with HIV/TB disease. Interventions will likely need to be instituted rapidly after TB diagnosis to alter outcomes.

The risk score was able to highlight patients at highest risk of death post-discharge, in addition to those at high risk of death during hospitalisation, and could be used to prevent too early discharges. Enhanced community support, including home visits, has been shown to reduce mortality after starting ART in advanced HIV [43], and could have a similar impact for HIV/TB patients. Current services in high-burden settings take a public health approach to service delivery, whereas prognostic risk scores can identify patients suitable for differentiated care [8].

The main aim of this risk score was to detect patients at high risk of early mortality who may benefit from interventions in addition to TB treatment. Although the discrimination of the model was not perfect, the sensitivity of the simplified score was 75%; the score did not identify 25% of patients who died within 2 months, and such patients would still receive standard-of-care management of HIV/TB. Proposed interventions to reduce mortality would have limited adverse events, so those deemed as 'high risk' by the score but surviving to 2 months are unlikely to come to significant harm from such interventions. However, if adjunctive interventions are found to reduce early mortality, better predictive biomarkers or more accurate predictive tools would allow more efficient use of resources through targeting of patients.

Limitations of our study include the potential for selection bias. In the STAMP trial standard-of-care arm, only patients started on TB treatment for clinical/radiological criteria or following a positive sputum Xpert result had stored urine retrieved for TB testing. Patients with otherwise undiagnosed TB who would have been urine test positive if they had been tested were not included in this study. Patients unable to provide consent, mostly due to being severely unwell and having altered consciousness, were also excluded. Although our risk score did not have optimal discrimination and calibration, performance was adequate and similar to that of other prognostic scores widely used in clinical practice (e.g., the Framingham cardiovascular risk score) [15,44]. Performance may have been reduced by categorising continuous variables for simplicity. TB drug resistance was not a predictor of mortality in this cohort; however, prevalence of rifampicin resistance was low in these settings. Not all established risk factors for mortality were characterised, leaving potential to improve on performance. Future studies could assess more detailed markers of physiology, as well as social and more distal risk factors.

Whilst the score is pragmatic and its constituent factors are widely available in hospitals in African regions with high HIV and TB burdens, it does rely on access to the TB-LAM lateral flow assay. There is now good evidence to support mortality reductions with the use of TB-LAM in HIV-positive patients admitted to hospital [4,5], and its use as a screening test has been incorporated into the latest guidelines in Malawi and South Africa. The assay has also been scaled up nationally in eSwatini, Kenya, and Uganda [45]. Missing data were common in the validation cohort. However, sensitivity analyses using multiple imputation gave similar results as the complete case analysis. We assumed patients lost to follow-up were alive at 2 months, although only 2% in the derivation cohort and 4% in the validation cohort were not followed up after hospital discharge. Our cohort did not include patients treated for TB without a positive diagnostic test, which remains common in HIV-positive patients admitted to hospital, and this patient group may be an important group for whom to apply risk
stratification and predictive scores. The biomarkers studied are imperfect predictors of mortality, and further research is needed to focus on better biomarkers to predict outcome.

Strengths of this study include that the derivation cohort and the LAM-RCT external validation cohort were nested within randomised controlled trials. Our predictive model had similar discrimination and calibration in the validation cohort, and was able to identify groups of patients with similarly increased odds of mortality. This was despite the validation cohort being from geographically distinct locations, collected at different times by different investigators, and with a lower overall mortality risk at 2 months. The factors required for the score can be obtained rapidly after admission.

In conclusion, we have developed and externally validated a clinical risk score capable of identifying, among patients admitted to hospital in settings with high HIV/TB burden, those with the highest risk of early mortality. This score could be a useful clinical and research tool, and could prove beneficial in identifying patients who would gain most from adjunctive interventions to reduce mortality. Further work to assess the impact of such risk scores, and to identify which interventions could potentially reduce mortality, is urgently needed if ambitious global targets to reduce TB mortality are to be met by 2025.

**Supporting information**

**S1 Appendix. Statistical analysis plan.** Prospective statistical analysis plan (version 1.2, 4 June 2017).
(PDF)

**S2 Appendix. TRIPOD checklist.**
(PDF)

**S3 Appendix. List of ethics committees that provided approval.**
(PDF)

**S1 Fig. Performance of clinical risk score in derivation cohort.** (A) Receiver operator curve of the predictive model: area under the curve = 0.70 (95% CI 0.63–0.76). (B) Calibration plot of observed probability of mortality plotted against predicted probability of mortality by the risk score multivariable regression model, with variables grouped into deciles based on predicted probability, and 95% CIs. Black dashed line shows perfect prediction. Hosmer-Lemeshow statistic *p* = 0.78.
(TIF)

**S2 Fig. Observed and predicted mortality for risk score values in the derivation cohort.**
The size of the blue circles representing observed mortality risk is proportional to the number of patients with that score. Predicted mortality risk is represented by the green line/triangles.
(TIF)

**S3 Fig. Risk score calculation for the simplified risk score to predict mortality.**
(TIF)

**S4 Fig. Urine TB-LAM grade and CD4 cell strata stratified by clinical risk score for the derivation cohort (n = 315).**
(TIF)

**S5 Fig. Calibration plot for the predictive model in the external validation dataset (n = 372).** Plot shows the observed compared to expected probability of risk for the external validation cohort as deciles based on risk score, with 95% CIs. Hosmer-Lemeshow statistic *p* = 0.13. *c*-Statistic (or area under the receiver operator curve) was 0.68 (95% CI 0.61–0.74). Dotted line
represents perfect prediction.

(TIF)

**S1 Table. Univariable analysis of continuous variables and associations with mortality.**

(PDF)

**S2 Table. Risk score and mortality data.** Number of patients surviving and patients dying, and observed and predicted mortality risk for (A) the full clinical risk score (based on the regression coefficients) and (B) the simplified risk score.

(PDF)

**S3 Table. C-statistic, 95% confidence intervals, and Hosmer-Lemeshow test for final model and risk scores in derivation and validation cohorts.**

(PDF)

**Acknowledgments**

The authors would like to acknowledge study staff and participants from the STAMP trial, the LAM-RCT, and the MSF study at Homa Bay County Hospital, Kenya.

**Author Contributions**

**Conceptualization:** Ankur Gupta-Wright, Elizabeth L. Corbett, Stephen D. Lawn, Katherine Fielding.

**Data curation:** Ankur Gupta-Wright, Keertan Dheda, Helena Huerga, Jonny Peter, Daniel Grint, Katherine Fielding.

**Formal analysis:** Ankur Gupta-Wright, Katherine Fielding.

**Funding acquisition:** Elizabeth L. Corbett, Stephen D. Lawn, Katherine Fielding.

**Investigation:** Ankur Gupta-Wright, Douglas Wilson, Joep J. van Oosterhout, Keertan Dheda, Helena Huerga, Jonny Peter, Maryline Bonnet, Melanie Alufandika-Moyo, Daniel Grint, Stephen D. Lawn, Katherine Fielding.

**Methodology:** Ankur Gupta-Wright, Elizabeth L. Corbett, Katherine Fielding.

**Project administration:** Ankur Gupta-Wright, Douglas Wilson, Joep J. van Oosterhout, Helena Huerga, Jonny Peter, Melanie Alufandika-Moyo.

**Supervision:** Ankur Gupta-Wright, Elizabeth L. Corbett, Douglas Wilson, Joep J. van Oosterhout, Keertan Dheda, Helena Huerga, Jonny Peter, Maryline Bonnet, Katherine Fielding.

**Writing – original draft:** Ankur Gupta-Wright, Elizabeth L. Corbett, Katherine Fielding.

**Writing – review & editing:** Ankur Gupta-Wright, Elizabeth L. Corbett, Douglas Wilson, Joep J. van Oosterhout, Keertan Dheda, Helena Huerga, Jonny Peter, Maryline Bonnet, Melanie Alufandika-Moyo, Daniel Grint, Katherine Fielding.

**References**

1. World Health Organization. Global tuberculosis report 2017. Geneva: World Health Organization; 2017.

2. Ford N, Matteelli A, Shubber Z, Hermans S, Meintjes G, Grinsztejn B, et al. TB as a cause of hospitalization and in-hospital mortality among people living with HIV worldwide: a systematic review and meta-analysis. J Int AIDS Soc. 2016; 19:20714. https://doi.org/10.7448/IAS.19.1.20714 PMID: 26765347
3. Gupta RK, Lucas SB, Fielding KL, Lawn SD. Prevalence of tuberculosis in post-mortem studies of HIV-infected adults and children in resource-limited settings: a systematic review and meta-analysis. AIDS. 2015; 29:1987-2002. https://doi.org/10.1097/QAD.0000000000000802 PMID: 26266773

4. Godfrey-Faussett P, Ayles H. Why are people living with HIV still dying of tuberculosis? Lancet. 2016; 387:1141-3. https://doi.org/10.1016/S0140-6736(16)00699-1 PMID: 27025319

5. Gupta-Wright A, Corbett EL, van Oosterhout JJ, Wilson D, Grint D, Alufandika-Moyo M, et al. Rapid urine-based screening for tuberculosis in HIV-positive patients admitted to hospital in Africa (STAMP): a pragmatic, multicentre, parallel-group, double-blind, randomised controlled trial. Lancet. 2018; 392:292-301. https://doi.org/10.1016/S0140-6736(18)31267-4 PMID: 30032978

6. Blanc F-X, Sok T, Laureillard D, Borand L, Rekacewicz C, Nerrienet E, et al. Earlier versus later start of antiretroviral therapy in HIV-infected adults with tuberculosis. N Engl J Med. 2011; 365:1471-81. https://doi.org/10.1056/NEJMoa1013911 PMID: 22010913

7. Uthman OA, Okwundu C, Gbenga K, Volmink J, Dowdy D, Zumla A, et al. Optimal timing of antiretroviral therapy initiation for HIV-infected adults with newly diagnosed pulmonary tuberculosis: a systematic review and meta-analysis. Ann Intern Med. 2015; 163:32-9. https://doi.org/10.7326/M14-2979 PMID: 26148280

8. McNairy ML, Abrams EJ, Rabkin M, El-Sadr WM. Clinical decision tools are needed to identify HIV-positive patients at high risk for poor outcomes after initiation of antiretroviral therapy. PLoS Med. 2017; 14(4):e1002278. https://doi.org/10.1371/journal.pmed.1002278 PMID: 28419097

9. Jarvis JN, Bicanic T, Loyse A, Namarika D, Jackson A, Nussbaum JC, et al. Determinants of mortality in a combined cohort of 501 patients with HIV-associated cryptococcal meningitis: Implications for improving outcomes. Clin Infect Dis. 2014; 58:736-45. https://doi.org/10.1093/cid/cit794 PMID: 24319084

10. Koss CA, Jarlsberg LG, Den Boon S, Cattamanchi A, Davis JL, Worodria W, et al. A clinical predictor score for 30-day mortality among HIV-infected adults hospitalized with pneumonia in Uganda. PLoS ONE. 2015; 10(5):e0126591. https://doi.org/10.1371/journal.pone.0126591 PMID: 25962069

11. Armstrong-James D, Copas AJ, Walzer PD, Edwards SG, Miller RF. A prognostic scoring tool for identification of patients at high and low risk of death from HIV-associated Pneumocystis jirovecii pneumonia. Int J STD AIDS. 2011; 22:628-34. https://doi.org/10.1258/ijsa.2011.011040 PMID: 22096046

12. Jacob ST, Pavlinac PB, Nakiyingi L, Banura P, Baeten JM, Morgan K, et al. Mycobacterium tuberculosis bacteremia in a cohort of HIV-infected patients hospitalized with severe sepsis in Uganda—high frequency, low clinical suspicion [corrected] and derivation of a clinical prediction score. PLoS ONE. 2013; 8(8):e70305. https://doi.org/10.1371/journal.pone.0070305 PMID: 23940557

13. Saunders MJ, Wingfield T, Tovar MA, Baldwin MR, Datta S, Zevallos K, et al. A score to predict and stratify risk of tuberculosis in adult contacts of tuberculosis index cases: a prospective derivation and external validation cohort study. Lancet Infect Dis. 2017; 17:1190-9. https://doi.org/10.1016/S1473-3099(17)30447-4 PMID: 28827142

14. Van Wyk SS, Lin HH, Claassens MM. A systematic review of prediction models for prevalent pulmonary tuberculosis in adults. Int J Tuberc Lung Dis. 2017; 21:405-11. https://doi.org/10.5588/ijtld.16.0059 PMID: 28284255

15. Liao Y, McGee DL, Cooper RS, Sutkowski MB. How generalizable are coronary risk prediction models? Comparison of Framingham and two national cohorts. Am Heart J. 1999; 137:837-45. PMID: 10220632

16. Bastos HN, Osório NS, Castro AG, Ramos A, Carvalho T, Meira L, et al. A prediction rule to stratify mortality risk of patients with pulmonary tuberculosis. PLoS ONE. 2016; 11(9):e0162797. https://doi.org/10.1371/journal.pone.0162797 PMID: 27636095

17. Nguyen DT, Graviss EA. Development and validation of a prognostic score to predict tuberculosis mortality. J Infect. 2018; 77:283-90. https://doi.org/10.1016/j.jinf.2018.02.009 PMID: 29649520

18. Pefura-Yone EW, Balkissou AD, Poka-Mayap V, Fatime-Abaicho HK, Enono-Edende PT, Kengne AP. Development and validation of a prognostic score during tuberculosis treatment. BMC Infect Dis. 2017; 17:251. https://doi.org/10.1186/s12879-017-2309-9 PMID: 28388895

19. Nguyen DT, Jenkins HE, Graviss EA. Prognostic score to predict mortality during TB treatment in TB/HIV co-infected patients. PLoS ONE. 2018; 13(4):e0196022. https://doi.org/10.1371/journal.pone.0196022 PMID: 29659636

20. Gupta-Wright A, Peters JA, Flach C, Lawn SD. Detection of lipoarabinomannan (LAM) in urine is an independent predictor of mortality risk in patients receiving treatment for HIV-associated tuberculosis in sub-Saharan Africa: a systematic review and meta-analysis. BMC Med. 2016; 14:53. https://doi.org/10.1186/s12916-016-0603-9 PMID: 27007773

21. Gupta-Wright A, Fielding KL, van Oosterhout JJ, Wilson DK, Corbett EL, Flach C, et al. Rapid urine-based screening for tuberculosis to reduce AIDS-related mortality in hospitalized patients in Africa (the
STAMP trial): study protocol for a randomised controlled trial. BMC Infect Dis. 2016; 16:501. https://doi.org/10.1186/s12879-016-1837-z PMID: 27659507

22. Karnofsky D. Nitrogen mustards in the treatment of neoplastic disease. Adv Intern Med. 1950; 4:1-75. PMID: 14770961

23. Royston P, Moons KGM, Altman DG, Vergouwe Y. Prognosis and prognostic research: developing a prognostic model. BMJ. 2009; 338:1373-7.

24. Kyeyune R, den Boon S, Cattamanchi A, Davis JL, Worodria W, Yoo SD, et al. Causes of early mortality in HIV-infected TB suspects in an East African referral hospital. J Acquir Immune Defic Syndr. 2010; 55:446-50. PMID: 21105258

25. Crump JA, Ramadhani HO, Morrissey AB, Saganda W, Mwako MS, Yang L-Y, et al. Bacteremic dis-seminated tuberculosis in sub-saharan Africa: a prospective cohort study. Clin Infect Dis. 2012; 55:242-50. https://doi.org/10.1093/cid/cis409 PMID: 22511551

26. Marcy O, Laureillard D, Madec Y, Chan S, Mayaud C, Borand L, et al. Causes and determinants of mor-tality in HIV-infected adults with tuberculosis: an analysis from the CAMELIA ANRS 1295-CIPRA KH001 randomized trial. Clin Infect Dis. 2014; 59:435-45. https://doi.org/10.1093/cid/ciu283 PMID: 24759827

27. Subbarao S, Wilkinson KA, Van Halsema CL, Rao SS, Boyles T, Utay NS, et al. Raised venous lactate and markers of intestinal translocation are associated with mortality amongst in-patients with HIV-asso- ciated TB in rural South Africa. J Acquir Immune Defic Syndr. 2015; 70:406-13. https://doi.org/10.1097/ QAI.0000000000000763 PMID: 26186506

28. Talbot E, Munseri P, Teixeira P, Matee M, Bakari M, Lahey T, et al. Test characteristics of urinary lipoar- abinomannan and predictors of mortality among hospitalized HIV-infected tuberculosis suspects in Tan-zania. PLoS ONE. 2012; 7(3):e32876. https://doi.org/10.1371/journal.pone.0032876 PMID: 22412939

29. Moons KGM, Royston P, Vergouwe Y, Grobbee DE, Altman DG. Prognosis and prognostic research: what, why, and how? BMJ. 2009; 338:1317-20.

30. Altman DG, Vergouwe Y, Royston P, Moons KGM. Prognosis and prognostic research: validating a prognostic model. BMJ. 2009; 338:1432-5.

31. Peter JG, Zijenah LS, Chanda D, Clowes P, Lesosky M, Gina P, et al. Effect on mortality of point-of-care, urine-based lipoarabinomannan testing to guide tuberculosis treatment initiation in HIV-positive hospital inpatients: a pragmatic, parallel-group, multicountry, open-label, randomised controlled trial. Lancet. 2016; 387:1187-97. https://doi.org/10.1016/S0140-6736(15)01092-2 PMID: 26970721

32. Huerga H, Ferlazzo G, Bevilacqua P, Kirubi B, Ardizzoni E, Wanjala S, et al. Incremental yield of includ-ing determine-TB LAM assay in diagnostic algorithms for hospitalized and ambulatory HIV-positive patients in Kenya. PLoS ONE. 2017; 12(1): e0170976. https://doi.org/10.1371/journal.pone.0170976

33. van Schaik S, O'Brien BC, Almeida SA, Adler SR. Interprofessional training ward pilot phase evaluation project report. Med Educ. 2014; 48:583-92. https://doi.org/10.1111/medu.12424 PMID: 24807434

34. Lawn SD, Gupta-Wright A. Detection of lipoarabinomannan (LAM) in urine is indicative of disseminated tuberculosis with renal involvement in patients living with HIV and advanced immunodeficiency: evi-dence and implications. Trans R Soc Trop Med Hyg. 2016; 110:180-5. https://doi.org/10.1093/trstmh/ trw008 PMID: 26884498

35. Kirenga BJ, Levin J, Ayakaka I, Worodria W, Reilly N, Mumbowa F. et al. Treatment outcomes of new tuberculosis patients hospitalized in Kampala, Uganda: a prospective cohort study. PLoS ONE. 2014; 9 (3):e90614. https://doi.org/10.1371/journal.pone.0090614 PMID: 24608875

36. Bates M, Mudenda V, Shibemba A, Kaluwaji J, Tembo J, Kabwe M, et al. Burden of tuberculosis at post mortem in inpatients at a tertiary referral centre in sub-Saharan Africa: a prospective descriptive autopsy study. Lancet Infect Dis. 2015; 15:544-51. https://doi.org/10.1016/S1473-3099(15)70058-7 PMID: 25765217

37. Karat AS, Omar T, Von Gottberg A, Tlali M, Chihota VN, Churchyard GJ, et al. Autopsy prevalence of tuberculosis and other potentially treatable infections among adults with advanced HIV enrolled in out-patient care in South Africa. PLoS ONE. 2016; 11(11): e0166158. https://doi.org/10.1371/journal.pone.0166158 PMID: 27829072

38. Hakim J, Musiime V, Szubert AJ, Mallewa J, Siika A, Agutu C, et al. Enhanced prophylaxis plus antire-troviral therapy for advanced HIV infection in Africa. N Engl J Med. 2017; 377:233-45. https://doi.org/ 10.1056/NEJMoa1615822 PMID: 28723333

39. Wallis RS, Hafner R. Advancing host-directed therapy for tuberculosis. Nat Rev Immunol. 2015; 15:255-63. https://doi.org/10.1038/nri3813 PMID: 25765201

40. Kolloli A, Subbian S. Host-directed therapeutic strategies for tuberculosis. Front Med. 2017; 4:171.

41. Kerkhoff AD, Lawn SD, Schutz C, Burton R, Boulle A, Cobelens FJ, et al. Anemia, blood transfusion requirements and mortality risk in human immunodeficiency virus-infected adults requiring acute
medical admission to hospital in South Africa. Open Forum Infect Dis. 2015; 2:ofv173. https://doi.org/10.1093/ofid/ofv173 PMID: 26730391

42. Andrews B, Semler MW, Muchemwa L, Kelly P, Lakhi S, Heimburger DC, et al. Effect of an early resuscitation protocol on in-hospital mortality among adults with sepsis and hypotension: a randomized clinical trial. JAMA. 2017; 318:1233–40. https://doi.org/10.1001/jama.2017.10913 PMID: 28973227

43. Mfinanga S, Chanda D, Kivuyo SL, Guinness L, Bottomley C, Simms V, et al. Cryptococcal meningitis screening and community-based early adherence support in people with advanced HIV infection starting antiretroviral therapy in Tanzania and Zambia: an open-label, randomised controlled trial. Lancet. 2015; 385:2173–82. https://doi.org/10.1016/S0140-6736(15)60164-7 PMID: 25765698

44. Moons KGM, Altman DG, Vergouwe Y, Royston P. Prognosis and prognostic research: application and impact of prognostic models in clinical practice. BMJ. 2009; 338:1487–90.

45. Treatment Action Group. Activists call on countries and donors to immediately scale up use of life-saving TB LAM test. New York: Treatment Action Group; 2018[cited 2019 Jan 5]. Available from: http://www.treatmentactiongroup.org/content/activists-call-countries-and-donors-immediately-scale-use-life-saving-tb-lam-test.