| Literature DB >> 33989317 |
Nicola Foster1,2,3, Lucy Cunnama1, Kerrigan McCarthy4,5, Lebogang Ramma6, Mariana Siapka3, Edina Sinanovic1, Gavin Churchyard5,7, Katherine Fielding3,5, Alison D Grant3,5,8, Susan Cleary1.
Abstract
BACKGROUND: In South Africa, replacing smear microscopy with Xpert-MTB/RIF (Xpert) for tuberculosis diagnosis did not reduce mortality and was cost-neutral. The unchanged mortality has been attributed to suboptimal Xpert implementation. We developed a mathematical model to explore how complementary investments may improve cost-effectiveness of the tuberculosis diagnostic algorithm.Entities:
Mesh:
Year: 2021 PMID: 33989317 PMCID: PMC8121360 DOI: 10.1371/journal.pone.0251547
Source DB: PubMed Journal: PLoS One ISSN: 1932-6203 Impact factor: 3.752
Fig 1Simplified schematic of the model.
The figure is a simplified representation of the model, with the circular boxes representing health states and the square boxes representing intermittent states used to model shorter time step process. The model structure is presented in more detail in S1 Text. “A” refers to the decision-making process from a negative test result to starting treatment without bacteriological confirmation; “B” represents the decision to continue testing for tuberculosis (negative pathway) in those with a negative test result; “C” is the behaviour around starting treatment after a positive test results; “D” refers to the decision (based on an interpretation of the further diagnostic tests) to start treatment; and “E” refers to the decision to start tuberculosis (TB) treatment after being ‘out of care’. * The model structures following the tuberculosis test are replicated for each of the six patient types, those HIV negative (with and without tuberculosis), HIV positive not on antiretroviral therapy (with and without tuberculosis), and those HIV positive on antiretroviral therapy (with and without tuberculosis). ** The treatments states are replicated for drug-sensitive and multi-drug resistant tuberculosis treatment.
Summary of parameters and distributions.
| Definition | Mean and stratification | Distribution | Comments. References are listed as name of first author, year (Reference). |
|---|---|---|---|
| Gender | 59.9% female | Represents trial population. Churchyard. 2015 [ | |
| Age (IQR) | 37 (29–48) years | Fixed | Represents trial population. Churchyard. 2015 [ |
| Initial population disease characteristics | HIVneg 0.314 (0.030); | Dirichlet | From trial population. Churchyard 2015 [ |
| HIVpos 0.531 (0.015); | |||
| ART 0.155 (0.005) | |||
| CD4 count in those with HIV (IQR) | 315 (192–480) cells/μL | Represents the microscopy arm of the trial population. Churchyard. 2015 [ | |
| True TB prevalence (includes bacteriologically confirmed -, clinical—and undiagnosed TB) in the microscopy arm of the study. | 13.0% | Fixed | Estimated from XTEND trial and model calibration. Churchyard. 2015 [ |
| Proportion of patients diagnosed with drug-resistant TB, any diagnosis | 4.0% (8/195) | Represents trial population. Churchyard. 2015 [ | |
| Proportion of patients starting MDR-TB treatment | 2.0% (3/195) | Represents what was observed in the XTEND trial. Churchyard. 2015 [ | |
| Probability of a positive Xpert test result if symptomatic and able to provide a sputum sample, mean (standard deviation) | HIVneg 0.077 (0.03); | Dirichlet | Estimated from XTEND trial. Churchyard. 2015 [ |
| HIVpos 0.132 (0.05); | |||
| ART 0.135 (0.03) | |||
| Probability of TB if patient had a positive test result | HIVneg 0.877; | Fixed | Estimated based on GX sensitivity 0.86 in HIVneg; 0.79 in HIVpos, 0.94 for Rif resistance, and GX specificity of 0.99 in HIVneg, HIVpos, 0.98 for Rif resistance. Steingart 2014 [ |
| HIVpos 0.936; | |||
| ART 0.938 | |||
| Probability of TB if patient had a negative test result | HIVneg 0.012; | Fixed | Unobserved parameter, estimated from model calibration. Based on GX sensitivity 0.86 in HIVneg; 0.79 in HIVpos, 0.94 for Rif resistance, and GX specificity of 0.99 in HIVneg, HIVpos, 0.98 for Rif resistance. Steingart 2014 [ |
| HIVpos 0.113; | |||
| ART 0.114 | |||
| Probability of starting treatment within 30 days of a positive test result, mean (standard deviation) | HIVneg 0.882 (0.325); | Dirichlet | Estimated from XTEND trial. Churchyard et al. 2015 [ |
| HIVpos 0.802 (0.400); | |||
| ART 0.944 (0.236) | |||
| Probability of starting treatment within one month of a negative test result without further diagnostic tests | HIVneg_TB 0.535; HIVneg 0.002; HIVpos_TB 0.072; HIVpos 0.009; ART_TB 0.017; ART 0.003 | Fixed | Probability of starting treatment was estimated from XTEND trial, whether this clinical decision was correct (treatment started in those with TB vs those without) was estimated through model calibration. Churchyard et al. 2015 [ |
| Probability of receiving further investigations after a negative test result | HIVpos_TB 0.041; HIVpos 0.041; ART_TB 0.073; ART 0.073 | Fixed | Estimated from XTEND trial. Churchyard. 2015 [ |
| Probability of starting TB treatment after further diagnostic tests | HIVpos_TB 0.212; HIVpos 0.027; ART_TB 0.217; ART 0.037 | Fixed | Estimated from XTEND trial and the model calibration. Churchyard. 2015 [ |
| Probability of starting TB treatment from ‘out of care’, by month: from all who do not start TB treatment within one month of the diagnostic test | |||
| Month 2 | HIVneg_TB 0.928; HIVneg 0.005; HIVpos_TB 0.164; HIVpos 0.000; ART_TB 0.100; ART 0.000 | Fixed | Curve estimated from XTEND trial. Assume that the behaviour from out of care remains the same. Churchyard. 2015 [ |
| Month 3 | HIVneg_TB 0.756; HIVneg 0.000; HIVpos_TB 0.066; HIVpos 0.000; ART_TB 0.207; ART 0.000 | Fixed | Curve estimated from XTEND. Assume that the behaviour from out of care remains the same. Churchyard. 2015 [ |
| Month 4 | HIVneg_TB 0.000; HIVneg 0.005; HIVpos_TB 0.146; HIVpos 0.000; ART_TB 0.148; ART 0.000 | Fixed | Curve estimated from XTEND. Assume that the behaviour from out of care remains the same. Churchyard. 2015 [ |
| Month 5 | HIVneg_TB 0.000; HIVneg 0.015; HIVpos_TB 0.064; HIVpos 0.000; ART_TB 0.000; ART 0.000 | Fixed | Curve estimated from XTEND trial. Assume that the behaviour from out of care remains the same. Churchyard. 2015 [ |
| Month 6 | HIVneg_TB 0.000; HIVneg 0.010; HIVpos_TB 0.060; HIVpos 0.000; ART_TB 0.000; ART 0.000 | Fixed | Curve estimated from XTEND trial. Assume that the behaviour from out of care remains the same. Churchyard. 2015 [ |
| Probability of starting MDR-TB treatment if diagnosed with MDR-TB | HIVneg 0.025; HIVpos 0.019; ART 0.000 | Dirichlet | Estimated from XTEND trial. Churchyard. 2015 [ |
| Probability of drug sensitive TB regimen started if TB treatment started | HIVneg 0.952; HIVpos 0.969; ART_TB 0.834 | Dirichlet | Estimated from XTEND trial. Churchyard. 2015 [ |
| Probability of MDR-TB regimen started if TB treatment started, mean (standard deviation) | HIVneg 0.039 (0.208); HIVpos 0.023 (0.002); ART 0.000 (0.000); | Dirichlet | Estimated from XTEND trial. Churchyard. 2015 [ |
| Average life expectancy at birth, South Africa | 63 years | Fixed | From the rapid mortality surveillance report 2014. Assumes that HIVpos patients who are on ART when they enter the model would have the same life expectancy as the general population (varied in the sensitivity analysis). HIV specific mortality considered in model through probabilities. Dorrington. 2015 [ |
| All-cause mortality in those without TB, monthly, mean (standard deviation) | HIVneg 0.001 (0.0005); HIVpos 0.002 (0.000); ART 0.001 (0.001) | Dirichlet | From Statistics South Africa report (P0309.3), mortality and causes of death in South Africa: findings from death notification [ |
| Standardised mortality ratio for all-cause mortality in patients post-TB treatment | 3.76 | Fixed | Increased all-cause mortality in those with a previous episode of TB [ |
| Monthly mortality if living with TB, not currently receiving treatment, mean (standard deviation) | HIVneg 0.018 (0.020); HIVpos 0.132 (0.005); ART 0.039 (0.005) | Changes over time | Based on Tiemersma. 2011 [ |
| Monthly mortality on treatment for those with TB, mean (standard deviation) | HIVneg 0.002 (0.001); HIVpos 0.046 (0.002); ART 0.006 (0.003) | Changes over time | Andrews 2012 [ |
| Disability weights, mean (standard deviation) | HIVneg_TB 0.331 (0.057); HIVpos_TB 0.399 (0.070); HIVpos 0.221 (0.041); ART 0.053 (0.011); ART_TB 0.331 (0.057) | Beta | Salomon. 2015 [ |
| The disability weight is a factor reflecting the severity of disease. | |||
| Microscopy, mean (standard deviation) | $6.30 ($1.34) | Gamma | Cunnama 2016 [ |
| Xpert, mean (standard deviation) | $16.90 ($6.10) | Gamma | Cunnama 2016 [ |
| Sputum liquid culture, mean (standard deviation) | $12.90 ($2.26) | Gamma | Cunnama 2016 [ |
| Digital radiograph, mean (standard deviation) | $15.17 ($7.74) | Gamma | Foster. Unpublished. |
| First-line drug sensitivity test, mean (standard deviation) | $20.30 ($7.28) | Gamma | Cunnama 2016 [ |
| Second-line drug sensitivity test, mean (standard deviation) | $25.10 ($20.22) | Gamma | Cunnama 2016 [ |
| Provider cost of clinic visit for initial diagnosis and monitoring | $8.63 | Fixed | Vassall 2017 [ |
| Provider cost of clinic visit for treatment | $3.89 | Fixed | Vassall 2017 [ |
| Patient cost of clinic visit | $2.90 | Fixed | Foster 2015 [ |
| Guardian cost per clinic visit | $10.04 | Fixed | Foster 2015 [ |
| Cost of caregiver per day | $0.69 | Fixed | Foster 2015 [ |
| Resource use along the diagnostic pathway | Detailed input available from | Gamma | Estimated by disease progression. Reported in Vassall 2017 [ |
| Provider cost of drug sensitive TB treatment, episode | $192.99 | Fixed | Estimated based on patient movements through care observed in the trial. Reported in Vassall 2017 [ |
| Provider cost of multi-drug resistant TB treatment, episode | $10 802.66 | Fixed | Estimated based on patient movements through care observed in the trial. Reported in Vassall 2017 [ |
| Patient cost of drug sensitive TB treatment, episode | Cost of accessing care associated $459.16; | Time-dependent functions | Foster 2015 [ |
| Cost of illness $135.94 | |||
| Patient cost of multi-drug resistant TB treatment, episode | Cost of accessing care associated $3 592.27; | Time-dependent functions | Foster 2015 [ |
| Cost of illness $2 442.03 | |||
In the Table, a fixed distribution refers to a distribution one where no uncertainty interval is estimated in keeping with calibration practice in complex models. Furthermore, IQR = interquartile range; TB = tuberculosis; MDR-TB = multi-drug resistant tuberculosis; Xpert = Xpert MTB/RIF; HIVpos = individuals HIV positive not yet started on anti-retroviral therapy; HIVpos_TB = individuals HIV positive with tuberculosis; ART = individuals HIV positive started on anti-retroviral therapy; ART_TB = individuals HIV positive on anti-retroviral therapy with tuberculosis.
Summary of the investment scenarios modelled.
| Investment | Model implementation | Parameter, events or resource changes | Assumptions | |
|---|---|---|---|---|
| Reduction in initial LTFU (in | All patients with positive TB test results start treatment within one month of testing—simulating a point-of-care or a track-and-trace scenario with active follow-up of people with a positive TB test result. Synergies with investment in a community health worker programme. | ptxfpos = 1—pMort_m1 (stratified by HIV and TB status) | Probability of starting treatment after positive (in month 1), from: | Monthly conditional probabilities of starting treatment from ‘out of care’ were estimated from the trial in the base scenario (reported in |
| The probability of starting treatment from a positive test result was the remainder of all patients in that state after those who would die in that month had been subtracted. The mortality rate was stratified by HIV and TB status. | HIVneg: 0.882 to 1; | |||
| HIVneg_TB: 0.882 to 1; | ||||
| HIVpos: 0.802 to 1; | ||||
| HIVpos_TB: 0.802 to 1; | ||||
| ART: 0.944 to 1; | ||||
| ART_TB: 0.944 to 1 | ||||
| Empirical treatment from negative test result (in | The ability of healthcare workers to correctly act based on continued clinical symptoms, on the same day as the results visit (by giving TB treatment to those with test negative TB expressed as the sensitivity and specificity of that decision). This was based on the behaviour estimated from the microscopy arm of the model calibration and was applied to behaviour after a negative Xpert test result. | pnegpathfeg = 0 | Probability of the | Given the differences in health care worker behaviour after a microscopy test compared to a Xpert test result observed in the XTEND trial, we use the transition probabilities estimated from the microscopy arm of the trial [ |
| ptreatfneg = value estimated from reported behaviour in the control arm of the XTEND study [ | HIVpos: 0.027 to 0.000 | |||
| Assumed that all have at least one visit to a public health clinic (and associated costs) after a negative test result for treatment initiation. | HIVpos_TB: 0.212 to 0.000 | |||
| ART: 0.037 to 0.000 | ||||
| ART_TB: 0.217 to 0.000 | ||||
| Probability of | ||||
| HIVneg: 0.002 to 0.040 | ||||
| HIVneg_TB:0.054 to 0.270 | ||||
| HIVpos: 0.009 to 0.180 | ||||
| HIVpos_TB: 0.072 to 0.360 | ||||
| ART: 0.003 to 0.060 | ||||
| ART_TB: 0.017 to 0.090 | ||||
| Improvements in the test-negative pathway (in | HIV-positive people with negative test results get further investigations (radiograph/culture) for TB, and a proportion are started on TB treatment, simulating additional investment in improving access to further diagnostic tests. | ptreatfneg = 0 | Probability of | Similar to the previous scenario, we model a healthcare worker behaviour change scenario based on the difference in observed behaviour between the microscopy and Xpert arms of the study. This scenario simulates a situation where there is an increase in the proportion of patients who receive further investigations after a negative test result. Therefore, we reduced all empirical treatment to 0 and all eligible patients received a radiograph as part of the negative pathway. |
| pnegpathfneg = 1 (stratified by HIV and TB status) | HIVneg: 0.002 to 0.000 | |||
| treatfnegpath = 0.10 (no TB); 0.80 (with TB) | HIVneg_TB:0.054 to 0.000 | |||
| The probability of starting treatment is shifted from following a negative test result to the decision to order further diagnostic tests. The probability of starting treatment after the negative pathway was 10% in those without TB, and 80% in those with TB. | HIVpos: 0.009 to 0.000 | |||
| Assumed that every person will accumulate two visits to the public clinic during the negative pathway, and that each person getting further tests will get at least one radiograph. | HIVpos_TB: 0.072 to 0.000 | |||
| ART: 0.003 to 0.000 | ||||
| ART_TB: 0.017 to 0.000 | ||||
| Probability of the | ||||
| HIVpos: 0.041 to 0.900 | ||||
| HIVpos_TB: 0.041 to 0.900 | ||||
| ART: 0.073 to 0.900 | ||||
| ART_TB: 0.073 to 0.900 | ||||
| Probability of | ||||
| HIVpos: 0.027 to 0.100 | ||||
| HIVpos_TB: 0.212 to 0.800 | ||||
| ART: 0.037 to 0.100 | ||||
| ART_TB: 0.217 to 0.800 |
In the Table, the individual characteristics of the patients are labelled as HIVneg for people who are HIV negative and don’t have tuberculosis; HIVneg_TB for people who are HIV negative and have been diagnosed with tuberculosis; HIVpos for people who are HIV positive and don’t have tuberculosis; HIVpos_TB for people who are HIV positive and have been diagnosed with tuberculosis; ART represents the individuals who are HIV positive, on anti-retroviral therapy and don’t have tuberculosis; and ART_TB represents the individuals who are HIV positive, on anti-retroviral therapy and have been diagnosed with tuberculosis.
Costs (US$), outcomes and ICERs over three years (36 one-month cycles) in a cohort with an estimated TB prevalence of 13%.
| Status quo and five investment scenarios | TB service costs per symptomatic individual (US$) | Outcomes per symptomatic individual | ICERs: compared against the status quo | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| In cohort of 10 000, true TB treated (range) | Provider costs | Societal costs | DALYs and DALYs averted | Deaths and Deaths averted | Provider cost/ DALY averted (95% UI) | Societal cost/ DALY averted (95% UI) | Provider cost/ death averted (95% UI) | Societal cost/ death averted (95% UI) | |||||
| Total (95% UI) | Incr change from base (%) | Total (95% UI) | Incr % change (range) | Total DALYs (95% UI) | Incr DALYs averted % change (range) | Total deaths (95% UI) | Incr deaths averted % change (range) | (95% UI) | (95% UI) | (95% UI) | (95% UI) | ||
| Xpert (status quo) | 940 | 89.66 | --- | 169.94 | --- | 4.72 | --- | 0.133 | --- | --- | --- | --- | --- |
| (920; 960) | (87; 92) | (167; 173) | (4.6; 4.8) | (0.129; 0.136) | |||||||||
| Xpert plus reduction in initial LTFU (iLTFU) | 1010 | 92.42 | 2.76 | 178.19 | 8.25 | 4.56 | 0.16 | 0.128 | 0.005 | 17.42 | 51.86 | 601.40 | 1790.50 |
| (C and E) | (990; 1030) | (90; 95) | 3% | (175; 181) | 5% | (4.4; 4.7) | 3% | (0.125; 0.132) | 4% | (2.2; 117.6) | (18.5; 271.0) | (75.1; 3806.7) | (644; 8774) |
| Xpert plus treatment from negative (TfN) | 1140 | 110.78 | 21.12 | 256.36 | 86.42 | 4.04 | 0.68 | 0.115 | 0.018 | 31.40 | 128.45 | 1180.00 | 4826.60 |
| (A and E) | (1120; 1160) | (109; 113) | 24% | (253; 260) | 51% | (3.9; 4.1) | 14% | (0.112; 0.118) | 14% | (24.6; 40.5) | (107.2; 157.0) | (905.9; 1567.3) | (3939.0; 6079.2) |
| Xpert plus reduction in initial LTFU, and treatment from negative (iLTFU_TfN) | 1210 | 113.55 | 23.89 | 264.60 | 94.66 | 3.88 | 0.84 | 0.110 | 0.023 | 28.73 | 113.82 | 1061.70 | 4205.70 |
| (A, C and E) | (1190; 1230) | (111; 116) | 27% | (261; 268) | 56% | (3.8; 4.0) | 18% | (0.107; 0.113) | 17% | (23.5; 35.3) | (98.3; 133.3) | (853.3; 1333.3) | (3569.1; 5035.1) |
| Xpert plus improvements in the negative pathway (NP) | 1420 | 141.01 | 51.35 | 278.87 | 108.93 | 3.42 | 1.30 | 0.096 | 0.037 | 39.70 | 84.19 | 1387.70 | 2943.10 |
| (B, D and E) | (1390; 1450) | (139; 143) | 57% | (274; 284) | 64% | (3.3; 3.5) | 28% | (0.093; 0.099) | 28% | (35.2; 44.9) | (75.0; 94.8) | (1225.6; 1576.9) | (2608.4; 3334.0) |
| Xpert plus reduction in initial LTFU, and improvements in the negative pathway (iLTFU_NP) | 1480 | 142.99 | 53.33 | 285.97 | 116.03 | 3.28 | 1.44 | 0.092 | 0.041 | 37.02 | 80.55 | 1292.97 | 2813.00 |
| (B, C, D and E) | (1460; 1510) | (141; 145) | 59% | (281; 291) | 68% | (3.2; 3.4) | 31% | (0.089; 0.094) | 31% | (33.3; 41.3) | (72.8; 89.4) | (1155.8; 1449.9) | (2527.7; 3139.4) |
In the Table, Incr is the incremental change in costs or effectiveness from the base case. The base case in this analysis which represents the current status quo, Xpert as observed in the intervention arm of the XTEND study; dominant: less costly and more effective; dominated: more costly and less effective; The 95% uncertainty interval (UI) is shown in parentheses; ICER: Incremental cost-effectiveness ratio; DALYs: Disability Adjusted Life Years. In the scenario column, the capital letters refer to the decision points upon which the investment scenario acts, as shown in Fig 1.
Fig 2Societal service-level costs (US$) per symptomatic person per episode.
In the Figure, the cost of accessing care (Access) includes out of pocket and time costs incurred by patients and caregivers when accessing care; the cost of illness (Illness) includes the cost of caregiver’s time, the cost of patient’s time when unable to work as well as loan interest, assets sold and the cost of nutritional supplements. Xpert referes to the Xpert baseline; iLTFU (Xpert plus iLTFU) = additional investment to reduce pre-treatment loss-to-follow-up; TfN (Xpert plus TfN) = supporting clinical diagnosis of tuberculosis after a negative test results; Np (Xpert plus NP) = improving access to further tuberculosis diagnostic tests following a negative test result.
Fig 3Provider cost-effectiveness acceptability frontiers (CEAF) at various levels of transaction costs.
Where iLTFU refers to Xpert plus a reduction in initial loss to follow up scenario; TfN refers to the scenario modelling Xpert plus treatment from negative; Np refers to Xpert plus improvements in the negative pathway. The cost-effectiveness acceptability frontier (CEAF) expressing the uncertainty around the cost-effectiveness of investments, by showing which strategy is economically preferred at a range of cost-effectiveness thresholds (on the x-axis). The base case scenario for each of these comparisons is Xpert MTB/RIF, as observed in the XTEND trial. The graph is a plot of the proportion of individual runs that would be cost-effective for each intervention (y-axis) while restricting the options to only those that would be the most cost-effective (optimal) investment for at least one individual, against a range of cost-effectiveness thresholds (x-axis). As the threshold increase, the preferred option changes, the switch point being where the incremental cost-effectiveness ratio (ICER) value increases beyond the threshold [62]. The analysis is repeated at a range of transaction costs per transaction, thereby varying the costs needed to be invested to facilitate systems level change in line with the investment strategy.
Fig 4(A-C) Results from the univariate sensitivity analyses, showing the ten parameters with the greatest influence on the (A) provider cost, (B) the societal costs, and the (C) effectiveness (DALYs) of the base case (Xpert). The full results for these analyses are presented S1 Text. In each one-way analysis, one parameter was varied by a factor of 10 from the mean to produce the low and high estimates, with all other parameters kept constant. Where DALYs are disability adjusted life years; Prov refers provider; and Soc is societal. DS treatment is drug-sensitive treatment. MDR refers to multi-drug resistant tuberculosis.