Abstract
Introduction Thise study aims to assess the cost effectiveness of romosozumab versus teriparatide, both sequenced to denosumab, for the treatment of severe postmenopausal osteoporosis at very high risk of fractures in Mexican women.
Methods A Markov model was used to assess the relative cost effectiveness of 1 year of romosozumab versus 2 years of teriparatide, both sequenced to denosumab for a total treatment duration of 5 years. Outcomes for a cohort of women with a mean age of 74 years, a T-score ≤ -2.5 and a previous fragility fracture were simulated over a lifetime horizon. The analysis was conducted from the perspective of the Mexican healthcare system and used a discount rate of 5% per annum. To inform relative fracture incidence, the bone mineral density (BMD) advantage of romosozumab over teriparatide was translated into relative risks of fracture, using relationships provided by a meta-regression of osteoporosis therapy trials. Outcomes were assessed in terms of lifetime costs (2023 Mexican pesos), quality-adjusted life years (QALYs) and life-years gained (LYs).
Results Base case results showed that, compared with teriparatide/ denosumab, romosozumab/ denosumab reduced costs by $51,363 MXN per patient and yielded 0.03 additional QALYs and 0.01 LYs. Scenario analyses and probabilistic sensitivity analyses confirmed that results are robust to uncertainty in model assumptions and inputs.
Conclusions Results show that romosozumab/ denosumab produces greater health benefits at a lower total cost than teriparatide/ denosumab.
Introduction
Osteoporosis is a chronic condition characterized by reduced bone strength, leading to an increased susceptibility to fractures. The risk of fractures due to osteoporosis becomes more prominent with advancing age, particularly among postmenopausal women (1). These fractures have a significant economic impact globally, resulting in substantial direct and indirect costs (2). In Mexico, the direct cost of more than 75,000 fragility fractures in 2010 amounted to $256.2 million US dollars, and it was projected to grow by 41.7% by 2020 (3). The likelihood of experiencing a hip fracture over a lifetime was estimated to be 8.5%. Fragility fractures, which occur due to weakened bone strength from minimal trauma, can result in loss of independence for patients and increased burden for both patients and their caregivers (4). They are also associated with higher risks of disability, hospitalization, and mortality, which can persist for several years, especially in the case of hip fractures (5,6).
Approximately one-third of postmenopausal women will experience a fragility fracture associated with osteoporosis during their lifetime (7). Once a postmenopausal woman experiences her first fracture due to osteoporosis, her risk for subsequent fractures increases significantly, with a more than fivefold higher likelihood of suffering another fracture within one year (8–10). A recent study showed that nearly 18% of individuals experienced a second fragility fracture, and the median time to the occurrence of the second fracture was less than two years (11).
Severe or established osteoporosis is diagnosed when a previous fracture is present and bone mineral density (BMD), measured by dual-energy X-ray absorptiometry (DXA), is -2.5 standard deviations or more below the mean value of a health young adult (12). According to a recent consensus statement in Mexico, teriparatide and denosumab were recommended as treatment options for severe osteoporosis (13). Furthermore, until 2022 teriparatide was the only therapy approved for reimbursement in patients with severe postmenopausal osteoporosis at very high risk of fracture (14). Teriparatide, an anabolic hormone, is recommended for individuals with osteoporosis who have a history of previous hip or spinal fractures, are at a high risk of refracture, and have not responded to bisphosphonate treatments. However, it should be noted that teriparatide has not been proven to prevent hip fractures. Furthermore, in accordance with the Clinical Practice Guidelines in Mexico for the management of osteoporosis in postmenopausal women established by CENETEC (National Center for Excellence in Health Technology), hip fractures remain prevalent and do not exhibit a decline. These fractures are associated with significant morbidity, mortality, and financial burdens (15,16).
Romosozumab is a humanized monoclonal antibody which works by inhibiting the protein sclerostin. It has a dual effect on the bones: it increases bone formation and decreases bone resorption, unlike other bone-forming agents such as teriparatide, where both bone formation and resorption are increased. Compared to bone-forming agents or aniresoptives, romosozumab treatment duration is only one year. Romosozumab has undergone comprehensive investigation to evaluate its potential in reducing the risk of fractures, with a focus on two Phase III studies that incorporated primary and secondary fracture endpoints: the ARCH (17) and FRAME (18) studies. Additionally, a separate study assessed the percentage change in a real bone mineral density (BMD) at the total hip as an outcome measure (19) (STRUCTURE trial). Romosozumab has been shown to increase trabecular and cortical bone mass, improve bone structure and strength, and reduce the incidence of new vertebral fractures by 73% compared to placebo in just one year of treatment. In comparison to teriparatide, romosozumab has demonstrated to double the BMD at key sites: +4.4% in the spine, +2.4% in the hip, and +3% in the femur.
Economic modeling is frequently employed to assess the relative cost-effectiveness of pharmacological treatments for osteoporosis (20–22). Romosozumab has been subject to evaluation in various contexts within these studies (23–25). The findings of these assessments indicate that romosozumab offers a cost-effective alternative compared to teriparatide and alendronate as a first-line treatment option for postmenopausal women. However, the specific evidence regarding the cost-effectiveness of romosozumab in the Mexican healthcare system is currently lacking.
As a result, the objective of this study was to analyze the cost-effectiveness of romosozumab, in comparison with teriparatide, both administered sequentially with denosumab (romosozumab/denosumab versus teriparatide/denosumab), for women with severe postmenopausal osteoporosis who face a particularly high risk of fractures within the framework of the Mexican healthcare system. In Mexico, teriparatide is recommended following the lack of success with bisphosphonates. As a result, the sequential therapy employed did not involve an antiresorptive bisphosphonate; instead, it encompassed a selective modulator of oestrogen, receptors, denosumab, available in the Mexican national basic formulary (26).
Methods
Patient population
The patient population modelled in this analysis consisted of a cohort of postmenopausal osteoporotic women who were at very high risk for future fracture. In line with the ARCH trial, the model included patients with a mean age of 74 years with a femoral neck BMD T-score ≤ ― 2.5 and a history of fragility fracture. Based on previous clinical evidence, it was assumed that 50% of patients had a single previous fracture, and 50% had multiple prior fractures (27).
Time horizon, perspective, and discount rate
Given that osteoporosis is a chronic disease, we used a lifetime time horizon to capture relevant benefits and costs associated with treatment. An annual discount rate of 5% was applied to both costs and future outcomes, in accordance with the Mexican HTA guidelines (14). Alternative discount rates were considered in scenario analyses. The analysis was conducted from the perspective of the public healthcare payer in Mexico.
Comparators
Our model compared two arms:
In the first arm, patients received romosozumab 210 mg monthly for 12 months, and were then sequenced to denosumab, 60 mg every six months (romosozumab/denosumab).
In the second arm, patients received teriparatide 20 μg daily for 24 months, sequenced to denosumab, 60 mg every six months (teriparatide/denosumab)
The choice of these comparators relies on current guidelines for osteoporosis management in Mexico (13). Romosozumab was administered for one year (consistent with the ARCH and FRAME trials and the product label). Teriparatide was administered for 2 years consistent with the clinical guidelines. Patients in both arms were assumed to be treated for a total of 5 years; a commonly recommended duration for pharmacological osteoporosis therapy (28).
Model structure
A Markov state transition model with a 6-month cycle length was used to assess costs, life years (LYs), and quality-adjusted life-years (QALYs) associated with the intervention arms. The structure of the model was based on the model created by the International Osteoporosis foundation. Such a model has been used widely as a basis to provide economic analyses of osteoporosis (29–33). Markov models are especially appropriate to use when: 1) a problem involves a continuous risk over time 2) the timing of events is important 3) events happen more than once, which is the case for our decision-problem (34).
Seven Markov health states were considered (Fig1): baseline (at risk of fracture), clinical vertebral fracture, post clinical vertebral fracture, hip fracture, post hip fracture, “other” fragility fracture, and death. All patients started at the baseline state, and throughout each cycle, they were at risk of sustaining a fracture (hip, vertebral, and other fracture) or death. A year after the event, patients with a hip or vertebral fracture transitioned into the post-hip fracture or post-vertebral fracture state. On the other hand, patients with another type of fracture returned to the “baseline” state.
The structure of the model is hierarchical, and the transition of the patient between the fracture states is based on the severity of the fracture type. Regarding the highest cost and lower utility, the most severe fracture is hip fracture, followed by vertebral fracture. Therefore, once the patient has sustained a vertebral fracture, the patient cannot go back to the baseline, or to the “other” fracture state(s). Similarly, patients that have sustained a hip fracture, cannot go back to the baseline, just move to the post hip state, or die. The reason for having a hierarchical structure is due to sustaining a clinical vertebral or hip fracture affects the cost and utility for the rest of the patient’s life. To rephrase, long-term effects are incurred in hip and clinical vertebral health states (21,24). However, the hierarchical structure that the Markov model holds results in underestimating the fracture incidence. To correct this, lower hierarchy fractures were separately estimated by multiplying the number of subjects in each higher hierarchical state, and with the incidence rate of the lower hierarchy fracture type in the model population.
Model inputs
Fracture risk
The risk of sustaining a fracture in the model depends on three elements:
The risk for an individual in the general population of incurring a fracture,
The increased fracture risk associated with osteoporosis (the relative risk RR), and
A risk reduction, if any, attributed to an intervention.
The general population risk depends on age and sex. The risk of fracture relative to the general population depends on age, bone mineral density and previous fracture. Age-specific general population rates of hip and vertebral fracture were derived from Ettinger et al (35). The incidences of other osteoporotic fractures were derived from Melton et al (36). Other osteoporotic fractures include ribs, pelvis, shaft/distal humerus, proximal humerus, clavicle, scapula, sternum, distal forearm, tibia,1 and fibula fractures. For age-specific incidence, these values were linearly interpolated or extrapolated as required to produce fracture rates for each year of age. General population fracture rates used in the model are shown in Supplementary Table 1.
To estimate fracture risks in untreated patients with severe postmenopausal osteoporosis, general population fracture rates were adjusted for the lower BMD T-scores and higher prevalence of previous fracture in the modelled population. Fracture risks for the patient population were adjusted with RRs of subsequent fracture in patients with a prior vertebral fracture and RRs of subsequent fracture per standard deviation decline in BMD, as described in previous economic evaluations (24,32).
Treatment efficacy
For individuals undergoing treatment, the model incorporated data on the effectiveness of the treatment, specifically in terms of the RRs of experiencing fractures. These RRs were applied over time to the fracture rates observed in untreated very high-risk postmenopausal osteoporotic patients. The efficacy data utilized in the model provided separate RRs for hip fractures, new vertebral fractures, and nonvertebral fractures. The RRs associated with new vertebral fractures were utilized to determine the treatment-specific efficacy in preventing vertebral fractures. Similarly, the RRs of nonvertebral fractures informed the efficacy of the treatment in preventing other types of fractures, including those affecting the wrist and all other non-hip, non-vertebral fractures. The efficacy of preventing hip fractures was directly informed by the treatment-specific RRs of hip fracture.
For patients receiving romosozumab/denosumab, relative fracture incidence compared with no treatment were derived by applying the RRs for romosozumab/denosumab versus denosumab taken from a propensity-score matching (PSM) post-hoc analysis of the FRAME trial (37) to the RRs for denosumab versus placebo from a network meta-analysis of osteoporosis therapies (38) (Table 1). To allow for changing relative efficacy over time, parametric survival functions were fitted to time-to-event data for hip and nonvertebral fractures (FRAME trial). Therefore, RRs were calculated for each 6-month model cycle over the course of the 5-year treatment period. For each arm of the FRAME trial (romosozumab/denosumab and denosumab), model fitting was applied for each fracture type (hip and nonvertebral fractures). The parametric functional form with the lowest Akaike information criterion (AIC) was selected, from exponential, Weibull, log-logistic, log-normal, Gompert, generalized gamma and gamma distributions. Selected regression parameters are shown in Supplementary Table 2. This was not possible for vertebral fractures, given that accurate time-to-event data were not available. Therefore, the RR of vertebral fractures for romosozumab/denosumab versus denosumab were calculated from new vertebral fracture incidence data from the FRAME trial, shown in Table 1. Efficacy estimates for romosozumab/denosumab versus placebo from this exercise (i.e. applying RRs romosozumab/denosumab versus denosumab to RRs for denosumab versus placebo) are shown in Table 1.
For patients in the teriparatide/ denosumab arm, the efficacy in reducing fractures for this arm compared to no treatment was determined by applying the RRs of fracture for teriparatide versus romosozumab to the RRs for romosozumab/ denosumab versus placebo (described in the previous paragraph).
To derive the RRs of fracture for teriparatide versus romosozumab, BMD efficacy results from the STRUCTURE trial (19) were converted to fracture outcomes using the relationship between percentage total hip BMD change from baseline and RRs of hip, vertebral, and nonvertebral fracture on the log scale, provided by the meta-regression conducted by the FNIH Bone Quality Study (39). While the regression equations were not directly reported in the publication, fracture risk reductions associated with a 2%, 4%, and 6% improvement in BMD were available, which were used to reproduce the parameters. The slopes of these equations were used to translate the difference of 3.4% total hip BMD change from baseline for romosozumab versus teriparatide at 12 months in the STRUCTURE trial into RRs of fracture. BMD at the total hip was used in these calculations due to its high predictive value for fractures (39). Resulting RRs of fracture for romosozumab versus teriparatide are shown in Table 1. Full details of their derivation are shown in Supplmentary Table 3. Given that the STRUCTURE trial exclusively offers relative BMD efficacy data for romosozumab versus teriparatide at the 12-months, the model’s base case involves assuming that the fracture reduction advantage of romosozumab over teriparatide endures until patients in both arms transition to denosumab, with the assumed duration being 2 years. Following this period, the model posits equal efficacy between romosozumab/ denosumab and teriparatide/ denosumab. This assumption was explored in scenario analysis.
Finally, the STRUCTURE trial was conducted in subjects previously treated with bisphosphonates. Given that newer osteoporosis therapies are often used as secondary prevention treatments (41), this is likely to be consistent with the characteristics of patients receiving bone-forming agents in clinical practice. However, to assess the cost effectiveness of romosozumab/ denosumab versus teriparatide/ denosumab in a treatment-naïve population, a scenario analysis was conducted where BMD efficacy was taken from a romosozumab phase II trial (2.8% difference in total hip change from baseline at 12 months), conducted in patients who were not bisphosphonate pre-treated (40).
Treatment persistence
In practice, adherence to pharmacological osteoporosis treatments is often less than perfect (42). Nevertheless, by recognizing that the controlled conditions of a randomized trial may not accurately reflect treatment discontinuation in practical settings, the model’s base case assumes full persistence with all treatments. The model then explores the impact of imperfect persistence through scenario analysis.
Treatment offset time
Although anti-fracture efficacy is likely to persist for a period of time (offset time) after a treatment is stopped, there have been very few studies that report offset time and there is a lack of consensus on the duration of offset time (29,43). As an example, some evidence has shown that after 5 years’ treatment with alendronate followed by 5 years’ treatment with placebo, patients’ BMD remained at or above pre-treatment levels (43), indicating that the benefit of pharmacological treatment persists for some time. Therefore, our model assumes that after the 5-year treatment course, fracture reduction benefit declines linearly to 0 over the course of a further 5 years for both arms.
Mortality
General population all-cause mortality was informed by life tables for females in Mexico from the National Institute of Statistics and Geography (INEGI) (44). The model also accounted for the increased risk of mortality following a fracture. Two key assumptions were made regarding mortality following osteoporosis-related fractures: 1) 30% of the excess mortality following a fracture was attributable to the fracture itself, in line with previous analyses (21,32) and 2) the increased risk of mortality after hip and vertebral fractures was assumed to last for 8 years as per previous analyses. This duration of excess mortality only applied to hip and vertebral fractures as other fractures were assumed to only have effects in the first year of fracture. Due to lack of Mexico-specific data related to mortality due to hip, clinical vertebral and other fractures, the default values were derived using data from Sweden (20).
Health-related quality of life
The estimates of utility in the general population were obtained from the US (45). Utility multipliers specific to the type of fracture were applied to this background age- and sex-specific utility value. To account for the health-related quality of life (HRQoL) loss due to fracture, in the first year after hip, vertebral, and other fracture, and for the second and subsequent years after hip and vertebral fracture, utility multipliers were applied to utilities of the general population. Data specific to subsequent “other” fractures were not available. Therefore, the utility multiplier in the first year after “other” fracture was assumed to correspond to that of a distal forearm fracture. These values were taken from Svedbom et al. (46). HRQoL inputs used in the model are shown in Table 2.
Costs
The model included drug acquisition costs, treatment monitoring/administration costs, and direct medical costs due to fracture (Table 2). All costs were expressed in 2023 Mexican pesos (MXN), inflated to 2023 values where required using the Consumer Price Index (CPI) for Mexico (51). Drug acquisition costs were obtained from the manufacturer for romosozumab and from the list prices from the Governmental Electronic System for Purchases (Compranet) (47).
For treatment monitoring/administration costs, the model assumed that patients receiving treatment incurred the cost of a physician visit once a year, and the cost of a BMD measurement every 2 years, as per on the IMSS inpatient and outpatient unit costs (48). The model assumed that all patients treated with romosozumab required a monthly outpatient visit for subcutaneous injection administration (outpatient consultation visit). This cost was obtained by the “Hospital General Gea González” inpatient and outpatient unit costs (49). Fracture morbidity costs included the direct medical cost of fracture (50), applied for the year in which fracture occur.
Analysis
The model estimated total discounted lifetime costs, LYs (total and without fractures) and QALYs for each intervention, with cost-effectiveness assessed in terms of incremental cost-effectiveness ratios (ICERs).
Our base case analysis considered LYs without fractures gained as the effectiveness of the decision-making process. In Mexico, the willingness-to-pay threshold values are typically set at 1 to 3 times the Gross Domestic Product (GDP) per capita, which translates to approximately $200,000 to $600,000 Mexican Pesos per LY gained. Our results were assessed through a probabilistic model with 1,000 simulations, where parameters were varied simultaneously according to distributions representing their uncertainty (52). Cost-effectiveness acceptability curves (CEACs) were derived to summarize the proportion of probabilistic iterations in which each comparator was cost-effective across a range of willingness to pay per QALY-gained thresholds. In addition, sensitivity analyses using the deterministic model were performed to assess the sensitivity of results to changes in individual parameters. Parameters were varied using published confidence intervals or standard errors, where available, and by 20% above and below point estimates where measures of uncertainty were unavailable. Cost-effectiveness was assessed with the incremental net monetary benefit (INMB) for each deterministic sensitivity analysis.
Scenario analyses were conducted to assess the impact of using alternative model assumptions. Scenario analyses around modelling assumptions were considered to account for annual discount rates of 0 or 3%, treatment efficacy rates estimated from parametric models, modified treatment offset time (1 year), increased excess mortality duration, using QALYs as instead of LYs, using a different time horizon (5 years), and reduced duration of fracture reduction benefits associated with romosozumab/ denosumab. We also considered relaxing the assumption of treatment persistence, i.e., treatment discontinuation. Persistence data for denosumab were taken from an analysis of osteoporotic patients in US Medicare Fee-for-Service data reported by Singer et al. (53). Persistence with romosozumab and teriparatide were reported by Chien et al. (54).Cost-effectiveness was estimated across each intervention for each scenario.
Results
Base case results
The total and disaggregated base case model results are presented in Table 3. Romosozumab/ denosumab yielded 7.85 LYs without fractures, and a total cost of $213,925 MXN. Compared to teriparatide/ denosumab, it was associated with a lifetime gain of .02 LYs, and a cost reduction of $46,61 MXN. Except for treatment monitory costs, all other costs (drug and morbidity costs) were higher for those patients treated with teriparatide/ denosumab. Therefore, the improvements in QALYs and cost reductions for romosozumab/ denosumab were driven by a reduction of the expected number of fractures (2542.3 per 1000 patients versus 2572.2 for teriparatide/ denosumab). When considering QALYs, romosozumab/ denosumab mean estimate was 6.09 and teriparatide/ denosumab 6.07. Consequently, teriparatide/ denosumab was dominated by romosozumab/ denosumab independently of the measure of effectiveness.
In our probabilistic sensitivity analyses, 97% our simulations yield lower costs and more LYs without fractures for romosozumab/ denosumab when compared to teriparatide/ denosumab (Figure 2). The probability of romosozumab/ denosumab being cost-effective was 100%, independently of the threshold chosen (Figure 3).
Fx, fracture; Vert, vertebral. Death can occur from any other state. For simplicity, specific transitions are not shown.
LY, life year
LY, life year
Deterministic sensitivity analyses and scenario analyses
Deterministic sensitivity analysis results demonstrated that the model outcomes were the most sensitive to the RR for vertebral fractures between denosumab vs placebo. None of the cost variations proposed in this study had an impact on the model outcomes. Moreover, romosozumab/ denosumab was consistently dominant, as it represented the highest INMB relative to teriparatide/ denosumab for each run of the one-way sensitivity analysis (see Figure 4).
FRAME, FRActure study in postmenopausal woMen with ostEoporosis; BMD, bone mineral density; Fx, fracture; QALY, quality-adjusted life year; RR, relative risk; SD, standard deviation; vert, vertebral
The incremental costs, LYs, and ICER of romosozumab/ denosumab relative to teriparatide/ denosumab for all scenarios are presented in Table 4. All scenarios yielded similar results to the base case, where romosozumab/ denosumab yielded additional LYs without fractures and fewer costs relative to teriparatide/ denosumab (i.e., remained the dominant intervention).
Discussion
This economic evaluation assessed the lifetime cost-effectiveness of 1 year of romosozumab sequenced to 4 years of denosumab versus teriparatide (2 years) sequenced with denosumab (3 years), for the treatment of severe osteoporosis in postmenopausal women in Mexico with a history of osteoporotic fracture and who are at very high risk for future fracture. Our base case and in all scenarios results indicate that romosozumab dominates teriparatide; meaning that it produces a gain in LYs (and QALYs) at a lower cost. This is due to the lower drug acquisition cost of romosozumab and the superior BMD efficacy of romosozumab which, translated to fracture reduction efficacy, produces both a QALY gain and a further cost reduction from avoided fractures. Deterministic and probabilistic sensitivity analyses show that model results are robust to uncertainty in parameters. When comparing romosozumab with teriparatide, romosozumab remains cost effective at a threshold of $600,000 per LY without fractures across all deterministic sensitivity analyses, with probabilistic sensitivity analysis showing a 100% probability that romosozumab is cost effective across all assessed thresholds.
The inclusion of romosozumab within the Mexican public health institutions could contribute to reduce the clinical and economic burden associated with osteoporotic fractures. Moreover, by presenting a simple administration scheme with two injections once a month, in comparison with daily injections of teriparatide, romosozumab facilitates adherence to treatment for osteoporosis, and confers more comfort to the patient.
Our cost-effectiveness analysis model is equivalent of the one that was submitted to the Canadian Agency for Drugs and Technologies in Health (CADTH) (55), and the General Health Council (Consejo de Salubridad General [CSG]) for reimbursement purposes. Although the cost effectiveness of romosozumab has already been estimated in different countries, our findings contribute to Mexico’s pharmacoeconomic literature of osteoporosis. To our knowledge, there has been only one study found that compared romosozumab with teriparatide (24). However, instead of being sequenced with denosumab, Hagino et al study used alendronate agent as the subsequent treatment (24). This specific study was conducted in Japan, and included women with severe postmenopausal osteoporosis that was treated with bisphosphonates as candidates. Hagino et al study found that romosozumab/ alendronate produced greater health benefits at a lower economic cost than teriparatide/ alendronate, which is consistent with our findings.
Limitations and future work
As with all analyses based on economic models, this evaluation has several limitations. The most important one is that this model uses a surrogate outcome – BMD – to estimate relative fracture incidence for romosozumab versus teriparatide, rather than directly observed fracture outcomes. However, as has been demonstrated by the FNIH meta-regression, treatment-induced changes in BMD are a strong predictor of fracture outcomes. Moreover, this method is preferable to estimating relative fracture outcomes via indirect methods which, as previously discussed, would necessitate synthesis of heterogeneous trials, resulting in a high degree of uncertainty.
A further limitation is that direct BMD efficacy data for romosozumab versus teriparatide are only available for up to 12 months in the STRUCTURE and romosozumab phase II trials. This is pertinent because romosozumab and teriparatide are provided for different durations (maximum of 1 year and 2 years, respectively) prior to sequencing to denosumab, meaning there is uncertainty in relative efficacy over time. Trial-based evidence also shows that teriparatide continues to produce an improvement in total hip BMD after 12 months of treatment (57). However, even under the conservative assumption that both regimens are equally efficacious after 12 months, romosozumab/ denosumab produces equivalent health outcomes at a lower cost than teriparatide/ denosumab.
The hierarchical nature of the model structure is also a limitation. Because patients cannot experience less severe fractures once they have had a more severe fracture, the model most likely underestimates the number of fractures in the population over their lifetime. Additionally, there is uncertainty in the duration of treatment “offset time” (the duration of fracture reduction benefit after treatment discontinuation). While the duration of the offset time was based on clinical evidence (43), it is not possible to precisely quantify the duration of the offset time. This assumption was tested through scenario analyses.
Another important limitation is that international data was used where appropriate local data were not available. These included HRQoL, the baseline age of the population, and RRs used to adjust general population fracture rates for fracture history and BMD. Mexican HTA guideline specifies to perform cost-effectiveness analysis as opposed to cost-utility analysis given the lack of HRQoL studies done in the Mexican population. Nevertheless, when using LYs instead of QALYs, romosozumab was still a dominant option.
Despite the limitations, this analysis provides clear evidence of the value-for-money of romosozumab sequenced to denosumab versus teriparatide/ denosumab, with results that are robust to alternative deterministic and scenario analyses. These results would not have been possible without the approach of converting BMD efficacy to RRs of fracture, considering the inherent uncertainty associated with indirect comparisons of fracture outcomes. Results of this evaluation are specific to the Mexican healthcare system perspective. However, it is reasonable to expect romosozumab to be cost effective in any setting where the cost per course of romosozumab/ denosumab is lower than that of teriparatide/ denosumab, since the BMD advantage of romosozumab over teriparatide guarantees a QALY gain and lower total cost. More research is needed to facilitate future cost-effectiveness analyses of bone-forming agents sequenced antiresorptives. Lastly, further research on the effectiveness of treatment in other populations (i.e., male osteoporosis, patients with secondary forms of osteoporosis) could inform the potential to extrapolate the results of this analysis.
Conclusion
This is the first economic model that evaluates the cost effectiveness of romosozumab/ denosumab as treatment for patients with postmenopausal osteoporosis from a Mexican perspective. The results from this cost-effectiveness analysis indicates that the strategy of using romosozumab for 1 year sequenced to 4 years of denosumab produces a greater number of QALYs (and LYs), reduce the number of fractures and lower the total cost when compared to 2 years of teriparatide, sequenced to 3 years of denosumab. According to probabilistic, deterministic, and scenario analysis, there is enough evidence to claim that romosozumab/ alendronate is likely to be the most cost-effective option at any decision-maker threshold. Due to this result, romosozumab should be considered as a reimbursable medication from Mexico’s public drug plan. This will aid in treating women with post-menopausal osteoporosis who has a history and is also at high risk for osteoporotic fractures.
Transparency
Author Contributions
TAM, EY and JPDM contributed to the design and implementation of the research, to the analysis of the results and to the writing of the manuscript; ES and LMC to the writing of the manuscript.
Declaration of financial/other relationships
TAM, LG, EY and LMC are employed by Amgen. JPDM has served as a consultant for Amgen and has received research grants from them.
Funding
This work was financed by Amgen, Mexico City, Mexico.
Availability of data and materials
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Data Availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request