Cost-effectiveness analysis of five drugs for treating postmenopausal women in the United States with osteoporosis and a very high fracture risk.

PURPOSE: Five strategies were recommended by the American Association of Clinical Endocrinologists/American College of Endocrinology (AACE/ACE) guidelines for the treatment of postmenopausal osteoporosis (PMO) patients with a very high fracture risk. We aimed to assess their cost-effectiveness in the United States (US).
METHODS: A microsimulation Markov model was created to compare the cost-effectiveness of five treatment strategies, including zoledronate, denosumab, abaloparatide, teriparatide, and romosozumab in PMO patients with a recent fracture from the healthcare perspective of the US. The data used in the model were obtained from published studies or online resources. Base-case analysis, one-way deterministic sensitivity analysis (DSA) and probability sensitivity analysis (PSA) were conducted for 65-, 70-, 75-, and 80-year-old patients.
RESULTS: In base case, at 65 years, zoledronate was the cheapest strategy. The incremental cost-effectiveness ratios (ICER, which represent incremental costs per QALY gained) of denosumab, teriparatide, abaloparatide, and romosozumab against zoledronate were $13,020/QALY (quality-adjusted years), $477,331 /QALY, $176,287/QALY, and $98,953/QALY, respectively. Under a willing-to-pay (WTP, which means the highest price a consumer will pay for one unit of a good of service) threshold of $150,000/QALY, denosumab and romosozumab were cost-effective against zoledronate. The PSA results showed that denosumab was the most cost-effective option with WTP thresholds of $50,000/QALY, $100,000/QALY and $150,000/QALY. The results were similar in other age groups. The DSA results indicated that the most common parameters that have important influence on the outcome were drug persistence, incidence of adverse events, the efficacy of drugs on hip fractures and the cost of the drug. CONCLUSION AND RELEVANCE: Among PMO patients with a very high fracture risk in the US, zoledronate is the cheapest strategy and denosumab is the most cost-effective choice among these five strategies.

Introduction

Osteoporosis is a disease characterized by low bone mass and microarchitectural deterioration of bone [1], which has become a significant public health issue. There are about 6–11% and 32–46% adults, age 50 years and older, in the US who have, respectively, osteoporosis or low bone mass [2], and it is estimated that more than 3 million osteoporotic fractures will occur in 2025, causing more than $25 billion in Medicare costs [3].

In 2020, the AACE/ACE recommended that the treatment of osteoporosis should be stratified by the patient’s fracture risk. For patients with a very high fracture risk (e.g., patients with a recent fracture, fractures while on approved osteoporosis therapy or multiple fractures), zoledronate, denosumab, abaloparatide, teriparatide, or romosozumab are recommended for initial treatment [4], as these drugs are more effective in increasing bone mineral density (BMD) and reducing the incidence of fractures [5, 6]. As the guidelines did not establish priorities for drug selection, and it is unrealistic to conduct a head-to-head trial that includes all five drugs to compare differences in their efficacy to reduce fracture risk, clinicians may be uncertain when deciding on a patient’s specific treatment regimen.

The objective of this study was to assess the cost-effectiveness of zoledronate, denosumab, teriparatide, abaloparatide, and romosozumab for treating postmenopausal women in the US with osteoporosis and a very high fracture risk. It aimed to provide a reference for selecting drug treatments for this population, and help promote the rational allocation of medical resources to reduce the national financial burden.

We chose patients with recent fractures as a representative group of those with a very high fracture risk, as it is widely recognized that fracture risk decreases gradually over time since the last fracture occurred [7, 8]. Moreover, a number of studies have confirmed the importance of a recent fracture for imminent fracture risk [9, 10]. It is of great individual value to treat patients with osteoporosis who have had a recent fracture and this will help in reducing osteoporotic fracture incidence.

Methods

As osteoporotic fractures are ongoing over time, a microsimulation Markov model was constructed to model the progression of it [11, 12]. This report followed the Consolidated Health Economic Evaluation Reporting Standards (CHEERS) statement [13] and the guidelines for the conduct of economic evaluations of osteoporosis [14]. We adopted a 3% discount rate for both costs and effectiveness, in line with US guidelines for economic evaluations [15].

Target population

The population was postmenopausal women with osteoporosis of various ages (65, 70, 75, 80 years old) who had a recent fracture (including hip, vertebral, wrist or other osteoporotic fracture).

Model structure

We constructed a microsimulation Markov model (Fig. 1) consisting of six health states: ‘At risk of fracture’, ‘Hip fracture’, ‘Clinical vertebral fracture’, ‘Wrist fracture’, ‘Other osteoporotic fracture’, and ‘Death’ from the US health care perspective. The horizon of the model was lifetime to capture the long-term effect of treatment. Patients were followed from the age they entered the model until 100 years old or death. The cycle length was 6 months [16, 17]. All the patients were supposed to start the cycle with a state based on their fracture site. Every 6 months, the patient would come to a state transition, based on the transition probabilities between states.

Fig. 1

Model structure. Patients in any states could transit to ‘Death’

Fracture probabilities and mortality

Patients’ fracture risk was modeled based on three factors: the fracture risk of the general population, the relative risk of patients compared with the general population, and the risk reduced by treatment.

The annual incidence rates of hip, vertebral and wrist fractures in the general population were based on data used to adjust the incidence of fractures of the US version of FRAX [18]; other osteoporotic fracture rates were calculated based on other published research [3]. Then, we calculated the incidence of fractures among postmenopausal women with osteoporosis without prior fractures, using the relative risk (RR) described by Mori et al. [19]. The fracture incidence of patients with a prior fracture was calculated according to a Swedish study, which considered the impact of the number of prior fractures and the time since the last fracture on the risk of subsequent fractures at the same time [20]. The RR for any fracture in that study was used for wrist and other osteoporotic fractures here.

Background mortality was obtained from the 2018 US life table [21]. The mortality rate of patients without hip fracture was equal to the background rate of mortality. When a hip fracture occurred, there would be excess mortality during the year of the fracture and subsequent years [22]. As comorbidities play a significant role in excess mortality after a hip fracture [23, 24], we assumed 25% of excess mortality was directly caused by fractures. Increased mortality following clinical vertebral fractures has been found to be very similar than those of a hip fracture [25, 26]. A scenario analysis that also considered excess mortality caused by clinical vertebral fracture was conducted.

Treatment and efficacy

Five strategies, which were recommended by the AACE/ACE guidelines [4], were compared in the model: (a) a 5 mg intravenous injection of zoledronate for 6 years (ZOLE); (b) a 60 mg subcutaneous injection of denosumab every 6 months for 10 years, followed by alendronate, po, weekly for 10 years (DENO/ALN); (c) an 80 μg subcutaneous injection of abaloparatide daily for 2 years, followed by 70 mg of alendronate, po, weekly, for 10 years (ABLA/ALN); (d) a subcutaneous injection of 20 μg of teriparatide daily for 2 years, followed by 70 mg of alendronate, po, weekly for 10 years (TPTD/ALN); and (e) a 210 mg subcutaneous injection of romosozumab monthly for 1 year, followed by 70 mg of alendronate, po, weekly for 10 years (ROMO/ALN).

Treatment efficacy was obtained from a meta-analysis and data used in prior models [26-28]. The efficacy of non-vertebral fractures was used for ‘Other osteoporotic fracture’ for all the drugs [14], as well as wrist fractures for zoledronate and romosozumab because other data were not available [16].

After treatment was completed, efficacy decreased linearly to zero over a period of time equal to the length of treatment [19, 29]. As BMD returned to baseline level within 1 years after discontinuation of denosumab [30], the offset time of denosumab was set to 1 years. For DENO/ALN, ABAL/ALN, TPTD/ALN, and ROMO/ALN, we assumed that efficacy of denosumab, teriparatide, abaloparatide, and romosozumab would be maintained during sequential therapy [26, 31, 32], and would start to decline after treatment was completed or the patient dropped out of treatment.

Persistence

Persistence refers to the duration of time from the initiation to the discontinuation of therapy [33]. Persistence data were obtained from published research [34, 35]. Because of a lack of long-term persistence data for zoledronate, denosumab, and alendronate, we assumed the persistence of patients who completed 2 or 5years of treatment would continue until the end of the prescribed treatment.

Adverse events

The adverse events that were considered in the model included gastrointestinal events for alendronate, acute-phase reaction for zoledronate, cellulitis for denosumab and hypercalcemia for both teriparatide and abaloparatide. Other adverse events, such as subtrochanteric femoral fracture, were excluded because of their very low incidence [36], and they were supposed to have little influence on the results. The incidence of adverse events was obtained from published studies [37-39].

Costs

Only direct costs were included in the model. The costs in the model were the costs of drugs, fractures and administration and management of adverse events. We obtained drug prices from average sale prices or red book [40, 41].The fracture costs were obtained from published Medicare claims data [42].The costs of managing adverse events were referred to previously published economic analysis models [32, 37, 43]. It was assumed that 12% of patients with a hip fracture required long-term care and the cost directly attributable to fractures was 25% [29]. Patients receiving treatment were assumed to make a physician visit and have a DXA test every year, as guidelines recommend [4]. For the costs of physician visits (Current Procedural Terminology [CPT] codes 99,213) and DXA scans (CPT codes 77,080), we used the allowable charges based on national payment amounts from the 2021 Medicare Physician Fee Schedule [44].A renal function test was required for patients treated with zoledronate ever year [43, 44]. For patients taking alendronate, 0·041 and 0·021 extra physician visits were required during the first 6 months and each subsequent cycle, and a proton pump inhibitor was prescribed at each visit [37].

All costs were reported in US dollars and converted to 2020 US dollars based on the Consumer Price Index (CPI).

Utility

Background utility was equal to that of the health of US women based on the EQ-5D [45]. There was a disutility in the first and subsequent years after hip and vertebral fractures [46], and when a gastrointestinal event [37], acute-phase reaction [47] or cellulitis [37] occurred. The disutility of hip and vertebra fracture were obtained from a systematic review [46]. The disutility of gastrointestinal event [37] and acute-phase reaction [47] were obtained from published economic analysis. Utility of cellulitis were measured using time trade-off methods [37, 48].

For wrist and other osteoporotic fractures, we only assumed a disutility during the first year of fracture [46]. We did not consider a disutility associated with hypercalcemia, as per a previous study [32].

Analysis

Consistent with the recommendations for conducting economic evaluations [14], we conducted base case, one-way deterministic sensitivity analysis and probabilistic sensitivity analysis for each age group, and quality-adjusted life years (QALY) were used as the outcome. Costs and QALY were calculated for each strategy. DSA was conducted on drug, fracture, cost of the administration and management of adverse events, background utility and disutility caused by fracture and adverse events, excess mortality after fractures, discount rate, treatment effects, and persistence. In the DSA, all the variables varied over a plausible range, which were obtained from 95% confidence intervals or a variance of 20% from base-case values. PSA was conducted using distributions to assess the influence of the joint uncertainty surrounding the model variables. Most of the parameters were contained in the PSA. The range and distribution of the variables are presented in Table 1.

Table 1

Key model parameters

	Value	Range	Distribution	References
Annual fracture incidence of general population (per1000 person)
Hip fracture, age 65–69	2.03	N/A	N/A	[18]
Hip fracture, age 70–74	3.94	N/A	N/A	[18]
Hip fracture, age 75–79	7.93	N/A	N/A	[18]
Hip fracture, age 80–84	14.47	N/A	N/A	[18]
Hip fracture, age 85 +	26.05	N/A	N/A	[18]
Clinical vertebral fracture, age 65–69	2.33	N/A	N/A	[18]
Clinical vertebral fracture, age 70–74	4.73	N/A	N/A	[18]
Clinical vertebral fracture, age 75–79	5.23	N/A	N/A	[18]
Clinical vertebral fracture, age 80–84	6.22	N/A	N/A	[18]
Clinical vertebral fracture, age 85 +	10.95	N/A	N/A	[18]
Wrist fracture, age 65–69	8.22	N/A	N/A	[18]
Wrist fracture, age 70–74	8.24	N/A	N/A	[18]
Wrist fracture, age 75–79	8.35	N/A	N/A	[18]
Wrist fracture, age 80–84	8.70	N/A	N/A	[18]
Wrist fracture, age 85 +	8.49	N/A	N/A	[18]
Other osteoporotic fracture, age 65–69	9.00	N/A	N/A	[3]
Other osteoporotic fracture, age 70–74	12.82	N/A	N/A	[3]
Other osteoporotic fracture, age 75–79	18.87	N/A	N/A	[3]
Other osteoporotic fracture, age 80–84	23.86	N/A	N/A	[3]
Other osteoporotic fracture, age 85 +	30.74	N/A	N/A	[3]
Mortality of general population
Age 65–100	a	N/A	N/A	[21]
Excess mortality after a fracture, RR
Hip fracture, first year	2.87	2.52–3.27	Gamma	[22]
Hip fracture, subsequent years	1.73	1.56–1.90	Gamma	[22]
Proportion of excess mortality caused by fracture after fracture	0.25	± 20%b	Gamma	[23, 24]
RR of fractures for osteoporotic individuals without fracture compared with general population
Hip fracture, age 65–69	3.91	3.28–4.56	Gamma	[19]
Hip fracture, age 70–74	3.13	2.80–3.47	Gamma	[19]
Hip fracture, age 75–79	2.60	2.39–2.82	Gamma	[19]
Hip fracture, age 80–84	2.04	1.91–2.17	Gamma	[19]
Hip fracture, age 85 +	1.92	1.78–2.05	Gamma	[19]
Clinical vertebral fracture, age 65–69	2.59	1.19–4.27	Gamma	[19]
Clinical vertebral fracture, age 70–79	2.15	1.15–3.15	Gamma	[19]
Clinical vertebral fracture, age 80 +	1.82	1.12–2.41	Gamma	[19]
Wrist fracture, age 65–69	1.78	1.78–2.19	Gamma	[19]
Wrist fracture, age 70–79	1.60	1.60–1.88	Gamma	[19]
Wrist fracture, age 80 +	1.45	1.45–1.64	Gamma	[19]
Other osteoporotic fracture, age 65–69	2.19	1.78–2.59	Gamma	[19]
Other osteoporotic fracture, age 70–79	1.88	1.60–2.15	Gamma	[19]
Other osteoporotic fracture, age80 +	1.64	1.45–1.82	Gamma	[19]
Risk of refracture after prior fracture	c	N/A	N/A	[20]
Utility
General population utility
Age 65–69	0.811	0.800–0.822	Beta	[45]
Age 70–79	0.771	0.758–0.784	Beta	[45]
Age 80 +	0.724	0.701–0.747	Beta	[45]
Fracture-related utility (multiplier)
Hip fracture, first year	0.797	0.77–0.825	Beta	[46]
Hip fracture, subsequent years	0.899	0.885–0.91	Beta	[46]
Clinical vertebral fracture, first year	0.720	0.66–0.775	Beta	[46]
Clinical vertebral fracture, subsequent years	0.931	0.916–0.946	Beta	[46]
Wrist fracture, first year	0.940	0.91–0.96	Beta	[46]
Other osteoporotic fracture, first year	0.910	0.88–0.94	Beta	[46]
Adverse event-related disutility (multiplier)
GI	0.910	0.89–0.96	Beta	[37]
APR	0.35	± 20%b	Beta	[47]
Cellulitis	0.820	0.79–0.85	Beta	[37, 48]
Treatment efficacy, RR
Alendronate
Hip fracture	0.61	0.42–0.90	Beta	[27]
Clinical vertebral fracture	0.57	0.45–0.71	Beta	[27]
Wrist fracture	0.82	0.25–2.53	Beta	[28]
Other osteoporotic fracture	0.84	0.74–0.94	Beta	[27]
Zoledronate
Hip fracture	0.60	0.45–0.81	Beta	[27]
Clinical vertebral fracture	0.38	0.25–0.58	Beta	[27]
Wrist fracture	0.79	0.67–0.94	Beta	[27]
Other osteoporotic fracture	0.79	0.67–0.94	Beta	[27]
Denosumab
Hip fracture	0.56	0.35–0.90	Beta	[27]
Clinical vertebral fracture	0.32	0.22–0.45	Beta	[27]
Wrist fracture	0.840	0.17–4.00	Beta	[28]
Other osteoporotic fracture	0.80	0.67–0.96	Beta	[27]
Teriparatide
Hip fracture	0.64	0.25–1.68	Beta	[27]
Clinical vertebral fracture	0.27	0.19–0.38	Beta	[27]
Wrist fracture	0.24	0.02–2.01	Beta	[28]
Other osteoporotic fracture	0.62	0.47–0.80	Beta	[27]
Abaloparatide
Hip fracture	0.24	0.01–4.84	Beta	[27]
Clinical vertebral fracture	0.14	0.05–0.42	Beta	[27]
Wrist fracture	0.63	0.41–0.98	Beta	[26]
Other osteoporotic fracture	0.51	0.29–0.87	Beta	[27]
Romosozumab
Hip fracture	0.44	0.24–0.79	Beta	[27]
Clinical vertebral fracture	0.33	0.22–0.49	Beta	[27]
Wrist fracture	0.67	0.53–0.86	Beta	[27]
Other osteoporotic fracture	0.67	0.53–0.86	Beta	[27]
Persistence (%)
12 months
Zoledronate	100	80–100b	Beta	[34]
Denosumab	70.5	± 20%b	Beta	[34]
Teriparatide	63.4	± 20%b	Beta	[34]
Abaloparatide	63.4	± 20%b	Beta	[34]
Romosozumab	80	± 20%b	Beta	[35]
24 months
Zoledronate	33.9	± 20%b	Beta	[34]
Denosumab	41.2	± 20%b	Beta	[34]
Teriparatide	40.8	± 20%b	Beta	[34]
Abaloparatide	39.8	± 20%b	Beta	[34]
Alendronate
12 months	54	± 20%b	Beta	[35]
24 months	39	± 20%b	Beta	[35]
36 months	30	± 20%b	Beta	[35]
48 months	22	± 20%b	Beta	[35]
60 months	15	± 20%b	Beta	[35]
Cost
Treatment cost
Alendronate (annual)	120	± 20%b	Gamma	[41]
Zoledronate (annual)	600d	± 20%b	Gamma	[59]
Romosozumab (annual)	23,547	± 20%b	Gamma	[40]
Denosumab (annual)	2542	± 20%b	Gamma	[40]
Teriparatide (annual)	47,256	± 20%b	Gamma	[41]
Abaloparatide (annual)	22,224	± 20%b	Gamma	[41]
Fracture cost
Hip fracture	29,986	25,677–42,913	Gamma	[42]
Clinical vertebral fracture	8325	5775–15,975	Gamma	[42]
Wrist fracture	4577	2543–10,674	Gamma	[42]
Other osteoporotic fracture	14,144	10,086–26,314	Gamma	[42]
Annual long-term care after hip fracture	2577	0–5154	Gamma	[60]
Proportion of the long-term cost attributable to fracture	25%			[29]
Adverse effects cost
PPI	2·41	1.81–3.01	Gamma	[37]
APR	106	± 20%b	Gamma	[43]
Cellulitis (daily)	1458	1449–1466	Gamma	[37]
Hypercalcemia	208	166–249	Gamma	[32]
Management cost
Renal monitoring fee	26	± 20%b	Gamma	[43, 44]
DXA (CPT code 77,080)	38·73	± 20%b	Gamma	[44]
Physician visit (CPT code 99,213)	92·47	± 20%b	Gamma	[44]
Annual incidence of adverse events
APR (first year)	0·42	± 20%b	Gamma	[38]
APR (second year)	0·12	± 20%b	Gamma	[38]
APR (third year)	0·06	± 20%b	Gamma	[38]
Cellulitis	0·0013	± 20%b	Gamma	[37]
Hypercalcemia (teriparatide)	0·04	± 20%b	Gamma	[39]
Hypercalcemia (abaloparatide)	0·034	± 20%b	Gamma	[39]
Discount rate	0·03	0–0.05	Gamma	[15]

RR relative risk, GI gastrointestinal event, APR acute-phase reaction, PPI proton-pump inhibitor, DXA dual-emission X-ray absorptiometry, CPT current procedural terminology, N/A not applicable

aRefer to reference #21

bSensitivity values 20% lower and 20% higher compared with base case value, which was based on assumption

cRefer to reference #20 and Online Resource 1(Table 3–5)

dCost of zoledronate consisted drug cost and intravenous injection fee

The incremental cost-effectiveness ratio (ICER) was calculated from the difference in average costs divided by the difference in the average effectiveness of two strategies. And we used willing-to-pay thresholds of $50,000, $100,000 and $150,000 in PSA.

Analyses were conducted using TreeAge Pro 2021 software (TreeAge, Williamstown, Massachusetts) and R software (version 4·1·0, http://www.r-project.org).

Results

Model validation

The model predicted that 99·5%, 99·5%, 99·4% and 99·3% of people died before age 100 in a starting age of 65, 70, 75 and 80 years old, separately. It was close to the US life table [21]. Consistent with the incidence of new osteoporotic fractures among Medicare beneficiaries in the US [3, 4], for patients with a recent fracture, 17% were hip fractures, 23% were vertebral fractures, 13% were wrist fractures, and 47% were other osteoporotic fractures. According to previous researches, 2.4% and 10.2% of women ≥ 65 years had a subsequent hip fracture within 1 and 5 years, respectively, following their prior clinical fracture regardless the site of the prior fracture. [49] And the 25-year incidence rates of subsequent hip were 23.8% [50]. Our model predicted that without an intervention, the probabilities for 65, 70, 75 and 80 year-old patients to have a hip fracture within 1 year following prior fracture were 1.9%, 2.5%, 4.0% and 4.6%, respectively. And the probability to have a hip fracture within 5 years were 6.3%, 8.0%, 11.4% and 13.3%, respectively. Lifetime probability to have a hip fracture following prior fracture in 65, 70, 75 and 80 year-old patients with recent fracture was 17.80%, 19.80%, 22.40% and 23.50%. Overall, refracture rates predicted by our model were similar to epidemiological data.

Base case analysis

The results of the base-case analysis of 65-, 70-, 75- and 80-year-old patients are shown in Table 2. For the 65-year-old patients, zoledronate had the lowest cost and utility. Compared with zoledronate, the ICER of denosumab/alendronate, romosozumab/alendronate, abaloparatide/alendronate and teriparatide/alendronate were $13,020, $98,953, $176,287 and $477,331 per QALY gains, respectively. Teriparatide/alendronate was absolutely dominated by abaloparatide/alendronate, which means it had lower effectiveness, but higher costs. Extended domination means that the ICER for a given alternative strategy is higher than that of the next more effective alternative. Romosozumab/alendronate exhibited extended domination because, when compared with denosumab/alendronate, it had a higher ICER than abaloparatide/alendronate did. The results were similar in other age groups.

Table 2

The results of the base case analysis

	Cost ($)	QALYs	ICER($/QALY)
Aged 65 y
ZOLE	31,084	10.31	Comparator	Undominated
DENO/ALN	32,367	10.31	13,020	Undominated
ROMO/ALN	40,827	10.38	98,953	Abs. dominated
ABAL/ALN	60,767	10.36	176,287	Undominated
TPTD/ALN	99,179	10.31	477,331	Abs. dominated
Aged 70 y
ZOLE	28,236	8.43	Comparator	Undominated
DENO/ALN	29,520	8.53	12,772	Undominated
ROMO/ALN	37,975	8.54	88,438	Ext. dominated
ABAL/ALN	56,171	8.60	163,444	Undominated
TPTD/ALN	93,030	8.57	449,882	Abs. dominated
Aged 75 y
ZOLE	27,209	6.72	Comparator	Undominated
DENO/ALN	28,919	6.81	19,424	Undominated
ROMO/ALN	37,172	6.82	104,959	Ext. dominated
ABAL/ALN	52,986	6.86	184,384	Undominated
TPTD/ALN	86,795	6.84	515,168	Abs. dominated
Aged 80 y
ZOLE	21,924	5.06	Comparator	Undominated
DENO/ALN	23,398	5.13	20,793	Undominated
ROMO/ALN	31,771	5.14	124,770	Abs. dominated
ABAL/ALN	42,959	5.18	182,721	Undominated
TPTD/ALN	71,606	5.16	513,241	Abs. dominated

ZOLE zoledronate, DENO denosumab, ALN alendronate, ABAL abaloparatide, TPTD teriparatide, ROMO romosozumab, Abs. absolute, Ext. extended

Sensitivity analysis

DSA was performed separately to compare undominated strategies against zoledronate. In 65-year-old patients, for the comparison of denosumab/alendronate against zoledronate (Fig. 2), the parameters that had a significant effect on the outcome were the price of zoledronate, the persistence of denosumab and zoledronate. For the comparison of abaloparatide/alendronate against zoledronate, the parameters that had a significant effect on the outcome were the price of abaloparatide, the efficacy of abaloparatide on hip fractures and persistence of zoledronate. The results for patients in the 70-, 75- and 80-year-old age groups are shown in Online Resource 1 (Figs. 4–6).

Fig. 2

Results of deterministic sensitivity analysis, age 65 years. Tornado plots shows the lower and upper values for the ICER of denosumab/alendronate strategy (a, dotted line = $13,020/QALY) and abaloparatide/alendronate strategy (b, dotted line = $176,287/QALY) to zoledronate strategy

The scenario analysis that considered excess mortality caused by hip as well as clinical vertebral fracture showed very similar results to those of base case analysis. For 65-, 70-, 75- and 80-year-old age group, zoledronate had the lowest cost and utility. Teriparatide/alendronate was absolutely dominated by abaloparatide/alendronate. And romosozumab/alendronate exhibited extended domination. (Online Resource 1 Table 6). When dominated strategies were eliminated from the PSA, the results showed, with a 100% probability, that denosumab/alendronate was the most cost-effective strategy at three WTP thresholds ($50,000/QALY, $100,000/QALY and $150,000/QALY) (Fig. 3).

Fig. 3

Result of probability sensitivity analysis comparing zoledronate, denosumab/alendronate and abaloparatide/alendronate in 65 (a), 70 (b), 75 (c) and 80 (d) years old. Circle (blue) represents for zoledronate, square(red) represents for denosumab/alendronate and cross (yellow) represents for abaloparatide/alendronate

Discussion

In this study, we evaluated the cost-effectiveness of zoledronate, denosumab, abaloparatide, teriparatide, and romosozumab in the treatment of postmenopausal women with osteoporosis and a recent fracture. The base-case analysis results indicate that zoledronate has the lowest cost and utility in 65-, 70-, 75-, and 80-year-old patients. The cost-effectiveness of strategies when compared with zoledronate differed under different WTP thresholds. Under a WTP threshold of $150,000/QALY, denosumab and romosozumab were cost-effective compared with zoledronate. But when WTP thresholds was lowered to $100,000/QALY or $50,000/QALY, only denosumab was cost-effective compared with zoledronate. The DSA results indicate that drug persistence, incidence of acute-phase reactions, the efficacy of drugs on hip fractures and the cost of the drugs are the most common parameters that have important influences on outcomes. The PSA results indicate that the probability of denosumab being the most cost-effective strategy is 100% in all age groups under different WTP thresholds.

Some previous US studies that analyzed the cost-effectiveness of osteoporosis treatment included the same drugs we investigated in our study. Stuart et al. [51] evaluated the cost-effectiveness of denosumab versus zoledronate, teriparatide and other osteoporotic treatments in 78 year-old men who had a BMD T-score of − 2·12 and a vertebral prevalence of 23%. The results showed that denosumab was dominant over zoledronate and teriparatide. Parthan et al. [37], which evaluated the cost-effectiveness of denosumab versus oral bisphosphonates for postmenopausal osteoporosis, suggested that denosumab was cost-effective or dominant compared with other strategies. A study by Mori and colleagues compared two strategies, including alendronate alone and teriparatide followed by alendronate [29]. The results indicated that for women with osteoporosis who had a previous vertebral fracture, the ICERs of sequential teriparatide/alendronate over alendronate varied from $288,200/QALY to $441,700/QALY from ages 65 to 80, which all exceeded the WTP threshold of $150,000/QALY. Another study that compared the cost-effectiveness of abaloparatide and teriparatide in postmenopausal women with osteoporosis, which used discrete-event simulation [32], showed that abaloparatide was a dominant treatment strategy over teriparatide and yielded an ICER of $18,891/QALY against alendronate.

It should be noted that the results of the cost-effectiveness analysis are affected by the choice of target population and the strategies being compared, as well as the model used in the analyses. Hence, our findings cannot be directly compared with the findings of previous studies. However, to some extent, they are consistent with previous findings. Our results also agree that the cost-effectiveness of denosumab may be better compared with other strategies, and that teriparatide seems to be a strategy that is not cost-effective because of its high cost and being dominated by abaloparatide.

This is the first study to investigate the cost-effectiveness of romosozumab in the US, and it shows the US costs and QALYs of 1-year treatment of romosozumab followed by 10-year treatment of alendronate in 70 year-old women with osteoporosis with a recent fracture was $37,975 and 8·54 QALY. Research in Japan showed that in 78-year-old women with a T score ≤  − 2·5 and a previous fragility fracture, the costs and QALY for a 1 year of romosozumab followed by 4 years of alendronate were, respectively, $35,178 and 7·290 QALY [52]. In Sweden, the costs and QALY were $60,396 and 8·547 in 74-year-old women with a recent major osteoporotic fracture [35].

In line with the summary of a systematic review of cost-effectiveness analyses of osteoporosis [53], the DSA results of our study confirmed that drug effect, drug cost, and medication persistence were key drivers of cost effectiveness. Our results also indicate that the incidence of adverse events during treatment and its influence on utility play an important role in outcomes. In previous studies, one of the most frequently unreported recommendations from the guidelines for economic evaluations of osteoporosis is adverse events [53]. Our study further illustrated the importance of including adverse events in cost-effectiveness analyses.

There are a few limitations associated with our study. First, because epidemiological data and WTP thresholds vary across countries, our findings should be generalized to other countries, conservatively. Second, we used efficacy of non-vertebral fracture for other osteoporosis fracture and vertebral fracture for clinical vertebral fracture. But we included the efficacy in DSA and PSA, and the results showed that denosumab always was the most cost-effective strategy when parameters vary. Third, although most of the data in our model were extracted from US data sets, some parameters were obtained from other countries, such as the 5-year persistence of alendronate and the RR of fractures, which varies with the number of prior fractures and the stage after the last fracture. Besides, the RR of subsequent fracture after prior fracture did not adjust for BMD which might cause fracture risk being over-evaluated. Fourth, the long-term persistence of zoledronate, denosumab and alendronate was assumed in this study. This assumption was based on research indicating that dropout rates are highest shortly after the initiation of treatment and remain stable for more than 5 years after that [54, 55]. The results should be updated when better data are available. At last, we used alendronate as a sequential treatment because it is the most widely used drug for the treatment of osteoporosis and it is also the drug that has been most frequently investigated in studies of the cost-effectiveness of sequential therapy [26, 35, 52]. However, in clinical practice, patients might choose other drugs. Thus, cost-effectiveness analyses using drugs other than alendronate as sequential therapy should be conducted in the future. Finally, we chose patients with a recent fracture to represent patients with a very high fracture risk, so the results may differ for patients meeting other standards.

Despite of these limitations, our study has notable strengths. This is the first study to investigate the cost-effectiveness of treatment strategies for US patients with osteoporosis who have a very high fracture risk. In addition, the strategies we compared were developed according to the recommendation of the AACE/ACE guidelines that zoledronate, denosumab, abaloparatide, teriparatide, and romosozumab can be considered as initial options for patients in a very high fracture risk, and that these patients should be treated with zoledronate for 6 years or receive a sequential treatment that consists of 10 years of bisphosphonate, p.o., following 10 years of denosumab, 2 years of abaloparatide or teriparatide or 1 year of romosozumab [4]. Also, we considered sequential treatment for denosumab for the first time, as recommended by the AACE/ACE guidelines. In addition, to the best of our knowledge, this is first study to assess the cost-effectiveness of romosozumab in the US. Finally, we used and referenced the guidelines for economic evaluations in osteoporosis, which helps to improve the transparency and quality of the study [14].

It should be noted that the pandemic of Coronavirus disease 19 (COVID-19) has brought some impact on the management of osteoporosis [56-58]. For example, the prevalence of non-adherence of denosumab was higher during the COVID-19 lockdown than before [56]. Since medicine persistence is one of the key drivers of cost effectiveness. The cost-effective of osteoporosis drug treatment might be overestimated during COVID-19 period. Therefore, results should be interpreted cautiously under this special situation.

In conclusion, this study provides economic results comparing zoledronate, denosumab, teriparatide, abaloparatide, and romosozumab in the US, which indicate zoledronate is the cheapest strategy. When compared with zoledronate, the cost-effectiveness of other strategies was dependent on the WTP threshold. Denosumab is the most cost-effective option when all the strategies are compared together under WTP thresholds of $50,000/QALY, $100,000/QALY and $150,000/QALY.

Below is the link to the electronic supplementary material.

Supplementary file1 (DOCX 1418 KB)