Research design and methods
Cost-effectiveness analysis
Cost-effectiveness analysis consists of comparing the incremental costs to deliver an intervention with the long-term benefits derived from it. In the case of the diabetes screening programme, the costs to deliver the intervention were compared with the incremental difference in healthcare expenditures and quality-adjusted life years (QALYs) between those receiving the intervention, and an otherwise identical cohort. The intervention includes screening and diagnosis, followed by standard care consisting of pharmaceutical and lifestyle change recommendations. The control group remains unaware of their diabetes status until they develop complications, at which point clinical diagnosis occurs and treatment begins (figure 1). Thus, we examined the incremental health-related benefits of the programme relative to the costs to deliver it.
Our approach is summarised by the following equation:
The first term of the numerator in the equation depends on the number of additional contacts screened and the cost of treating them if they screen positive for diabetes, and to a much smaller degree, the cost of negative screening results. The second term of the numerator in1 depends on the costs of medical care averted or delayed in the Markov chain described below.
Screening parameters
We obtained data on the prevalence of diabetes among TB contacts from the Texas–Mexico border as part of an ongoing TB research programme.3 10 12 Because we do not have sufficient data on this population, we used estimated sensitivity (68%) and specificity (89%) of screening for diabetes from a universal screening programme in Brazil, in which 47% of patients were in a fasting state for a blood glucose screening, and 53% were not, with cut-offs for a positive result of 100 mg/dL and 140 mg/dL, respectively.13 A patient with a positive result would undergo a diagnostic oral glucose tolerance test (OGTT), which we assume has 100% accuracy and will detect a false positive from the screening stage. These data points were used to calculate the true-positive, true-negative, false-positive and false-negative rates for diabetes screening, which are all inputs to the model. These rates are relevant to the cost-effectiveness of the screening programme, because only true positive results will yield any benefit in terms of QALYs and averted costs. Screening those without diabetes (true negatives and false positives) as well as failing to detect those with diabetes (false negatives) would not prevent or delay complications but would incur the same screening costs as a true positive. Figure 1 represents the potential paths through a screening programme for those with and without diabetes.
Utilities for QALY calculations
Utilities for each health state are from Coffey17 and are summarised in online supplemental table S1.1. In cases where a patient has progressed through multiple complications on the same pathway (clinical nephropathy and end stage renal disease), only the most severe QALY decrement from that pathway is used. In cases where a patient has developed complications along multiple pathways, the QALY decrements are additive.
Transition probabilities of disease states
Baseline transition probabilities of disease states were based on the literature, primarily on data from the United Kingdom Prospective Diabetes Study (UKPDS) and the Framingham Heart Study. The probabilities are based on data from individuals with glycated hemoglobin (HbA1c) values of 9.0% on average at the time of diabetes diagnosis, where probabilities will be generally lower on reductions in HbA1c and vice versa.18 Baseline transition probabilities and hazard ratios for changes in HbA1c are presented in online supplemental table S1.2.
Epidemiologic data
Demographic data, prevalence of cardiovascular risk factors (BMI, smoking, self-reported hypertension) and of new and known diabetes were based on data from close contacts of new patients with TB identified in TB clinics in South Texas and adjacent Mexican border communities as part of research studies described previously (table 1).10 12 19
Economic data
Three types of costs were considered: diabetes screening and diagnosis, diabetes treatment and diabetes complications. Costs associated with blood glucose screening and OGTT are based on the 2020 Q4 Medicare Clinical Laboratory Fee Schedule.20 Costs of diabetes management are based on recent national literature and include two office visits, medications (non-insulin, insulin or both) and other diabetes equipment and supplies.21 22 We assumed that 10% of the cohort took insulin alone, 67% non-insulin treatments only (oral anti-diabetic drugs or diet) and 23% both. These assumptions were based on unpublished data from a Texas-based cohort of diabetes patients. The cost of diabetes-related complications included the cost of nephropathy, neuropathy, retinopathy, coronary heart disease (CHD) and stroke. Event-year (year 1) and annual costs (years 2 and beyond, until death) of complications considered in the model are based on Ward23 and updated to 2021 dollars using the Medical Care Index component of the Consumer Price Index (online supplemental table S1.3). All future costs were discounted to the current year at a rate of 3% per year.
Projecting healthcare
We used the Diabetes Cost-Effectiveness Model (CDC-RTI Model) that was jointly developed by the Centers for Disease Control and Prevention and Research Triangle Institute to build a Microsoft Excel-based simplified version and updated key parameters (all costs inputs are derived from more recent sources, and characteristics of TB contact population are used to adjust all transition probabilities) to reflect the population served and the impact of this specific intervention.14 15 24 The CDC-RTI Model is a deterministic Markov simulation model of disease progression for diabetes that projects the cost-effectiveness of interventions designed to impact diabetes-related outcomes. It stimulates the development of diabetes-related microvascular (ie, nephropathy, neuropathy and retinopathy) and macrovascular (ie, CHD and stroke) complications along five separate pathways. Model outcomes include numbers and rates of complications, deaths, cost of care and QALYs. Markov models are state-transition models whereby progression between disease states is governed by transition probabilities that depend on risk factors. In our model, each model cycle represents 1 year. At the beginning of each cycle, an individual can experience a transition along each complication pathway or remain in the same state, and they cannot experience another transition until the following cycle. Relevant risk factors in this case include glycaemic level (HbA1c), blood pressure, cholesterol, smoking status and disease duration. Outcomes are projected from the year of diagnosis to either death or age 95, effectively a lifetime horizon. The primary mechanism through which the intervention impacts health outcomes in our model is through its impact on HbA1c, using data from the literature on rates of complications at different HbA1c levels.18 Our previous research on diabetes screening among TB contacts from the Texas–Mexico border showed that the mean HbA1c at diagnosis of diabetes was 7.6% (table 1). We assume that clinical diagnosis of diabetes occurs 9–12 years after onset of initial metabolic abnormalities and decline of pancreatic beta cell function (mean: 10.5 years).25 We further assumed that opportunistic screening would cut that in half to 5.25 years, under the assumption that screening is equally likely to detect metabolic abnormalities (ie, pre-diabetes) at any preclinical stage.26 Prior to clinical diagnosis, HbA1c is assumed to increase by 0.2 percentage points per year to 8.7% at clinical diagnosis, comparable to the value found in the UKPDS (9.0%)27 and that used by the CDC/RTI group (8.9%).15 We conservatively estimated that HbA1c would decrease by 1.0 percentage points in the first 12 months, then increasing by 0.156 percentage points per year following diagnosis (for comparison, the UKPDS diet-only group saw an average HbA1c reduction of 2.1 percentage points in the first 12 months). The cost of the programme is calculated from a health system perspective, meaning that it does not differentiate between costs paid by different parties (government, insurance company, patient, etc). A 3% discount rate was applied to both QALYs and healthcare costs.28 The full technical report on the model is available in online supplemental file 2.
Sensitivity analysis
Because our analysis projects outcomes into the future using estimates for several key data points, sensitivity analyses were carried out to test the effects of varying key parameters on the base model findings. Univariate and probabilistic sensitivity analyses were performed. For the univariate analysis, key parameters were adjusted by 20% and the low and high estimates were plugged into the model individually while holding all other parameters constant at base model levels to determine the extent to which they impacted the results. For the probabilistic analysis, random values were drawn from selected probability distributions (online supplemental table S2.16), as a Monte Carlo simulation with 1000 iterations of 1000 patients. The distributions were chosen either based on the literature, or when sources could not be found, so that the bounds of the 95% CI were 20% away from the mean.
Data and resource availability
The dataset analysed in the current study are available from the corresponding author on reasonable request.
Patient and public involvement
Patients or the public were not involved in the design, or conduct, or reporting or dissemination plans of our research.