Real-world effectiveness of osteoporosis treatment in the oldest old.

We studied effectiveness of osteoporosis treatment in women older than 80 years, who often are not included in clinical trials. Treatments were as effective on bone density and fractures as in younger women.
INTRODUCTION: To study real-world effectiveness of osteoporosis treatment on BMD and fractures in the oldest old women (≥ 80 years) compared with women (60-79 years) in the clinical setting using Swedish health register data.
METHODS: National registers and data from DXA machines were used to study effectiveness of all available osteoporosis treatments in women 60-79 and ≥ 80 years using three approaches: (1) Total Hip BMD change up to 8 years after treatment start; (2) fracture incidence where patients served as their own controls, comparing the first 3 months after treatment start with the subsequent 12 months; and (3) comparison of fracture incidence post-fracture in women ≥ 80 years treated with osteoporosis treatment or calcium/vitamin D.
RESULTS: Analysis 1: Total Hip BMD increased by up to 6.7% and 7.7% in women 60-79 and ≥ 80 years old, respectively. The mean increase in BMD was 1.1%-units per year in both age groups. Analysis 2: Relative to the 3-month baseline, fracture incidence decreased during the subsequent 12 months of treatment. Incidence rate ratios were estimated at 0.65, 0.74, 0.29, and 0.81 for any, hip, vertebral, and non-hip-non-vertebral fracture, respectively. Analysis 3: A 24-month incidence of any fracture in women ≥ 80 years given post-fracture osteoporosis treatment was lower (HR = 0.78) than in women given calcium/vitamin D, but treatment allocation was not random, with lower mortality (HR = 0.51) in patients receiving OP treatment.
CONCLUSIONS: Osteoporosis medication in women > 80 years in clinical practice likely works, and the magnitude of effect is similar to what was estimated in younger women. The choice between osteoporosis treatment and calcium/vitamin D after fracture in women ≥ 80 years is not random but appears associated with the patient's health status and presence of vertebral fractures, rather than the known risk profile of sustaining a fracture at a high age.

Introduction

The fracture incidence rate among women aged 80 and older (80+), sometimes referred to as the “oldest old”, is considerably higher than among younger women [1-3]. Fractures in older individuals also entail more serious consequences for the patient and higher costs for society compared with fractures in younger individuals [4, 5]. With a growing elderly population, the incidence of fractures will increase over the coming decades, and the number of hip fractures, generally associated with high impact on morbidity and mortality, is estimated to double in the first half of this century [6-8]. The efficacy of the available pharmacological agents for the treatment of osteoporosis in increasing bone strength and reducing osteoporotic fracture risk is well-established, although this evidence is generally better for the prevention of vertebral fractures than that for non-vertebral and hip fractures. The bulk of the evidence of efficacy rests on RCTs in postmenopausal women between the ages of 50 and 80 years and the evidence of anti-osteoporotic efficacy in the oldest old have come primarily from subgroup analyses. Nonetheless, evidence of efficacy in the oldest old is supported by randomized studies, reviewed in an ESCEO consensus article by Rizzoli and colleagues [9]. However, patients treated in clinical practice are more heterogeneous, less adherent to treatment, may have comorbidities or other characteristics excluding them from RCTs, and may be treated for both shorter and longer periods than typically recommended. It is therefore of importance to study if the effectiveness reported in RCTs can be generalized to clinical practice. Notable studies of effectiveness in older patients have been published by Abelson et al. [10] and Axelsson et al. [11] who both reported risk reductions with alendronate, but with varying designs and definitions, as is often the case in real-world effectiveness studies. Osteoporosis treatments are frequently found to be under-prescribed in many countries [12], also in women who have sustained a fracture. Eleven percent of treatment-naïve Swedish women 70–79 years old received treatment within 12 months after a fracture, whereas the corresponding figure for women > 80 years is 6% [13]. Moreover, only 4% of women ≥ 80 years had undergone a DXA within 6 months after fracture [data on file] which may indicate that thorough risk assessment is unusual and that treatment decisions partly are based on other factors than those traditionally governing fracture risk. The reasons for the poor treatment provision are not clear but could relate to that osteoporosis sometimes is regarded to be a natural part of ageing, that other concomitant diseases should be prioritized, and that expected remaining life expectancy is too short to reap the benefits of fracture prevention.

The objective was to study the real-world effectiveness of osteoporosis treatment in the oldest old women (≥ 80 years) compared with younger women (60–79 years) in the clinical setting using both BMD and fracture as outcomes.

Methods

Source data, participants, and study design

The study used the Nordic Osteoporosis Research Dataset (NORD), which contains > 1 million patients with fracture, > 400,000 with osteoporosis treatment, and > 55,000 with BMD measurements. All data were linked via a unique patient ID and spans from 2000 to 2016. Data on diagnoses, pharmaceuticals, and mortality were extracted from the National Patient Register, the Swedish Prescribed Drug Register, and the Cause of Death Register, respectively. A DXA population was selected from three databases containing BMD data (Malmö, Uppsala and Linköping) and included women with at least one BMD measurement between the dates 1 July 2007 and 31 December 2015 while being ≥ 50 years old. The research was approved by the Stockholm ethics review board and the data holders.

This was a retrospective, observational cohort study using three distinct study cohorts drawn from NORD to study the effectiveness of osteoporosis therapy in the oldest old (≥ 80) compared with younger populations (60–79 years) and, in one of the analyses, with subjects receiving calcium/vitamin D. The three different approaches described below were chosen to complement each other to account for the use of non-randomized retrospective data, which are limited by bias to some extent. Women with a diagnosis of Paget’s disease at any time point during the study period or with a primary diagnosis of secondary malignant neoplasm of bone and bone marrow at any time point during the individual look-back period or follow-up were excluded. Identification of fractures was based on primary ICD-10 codes in inpatient and outpatient hospital care, including ER visits and outpatient referrals from primary care. Fractures coded for the same body site as a previous fracture were disregarded and counted as a follow-up visit if coded in the outpatient setting. Inpatient fractures at the same body site were counted as a rehospitalization related to the previous fracture, unless ≥ 6 months had passed from the previous fracture. When adjusted, statistical models for fracture outcomes were controlled for age, Charlson comorbidity index, glucocorticoid use, pre-packaged drug dispensing, index fracture type, exposure to drugs increasing the risk of falls, and use of proton pump inhibitors. Models with TH BMD as outcome were adjusted for baseline T-score ≤ − 2.5 SD, glucocorticoid use, and use of proton pump inhibitors.

Analysis 1: T-score and percent change in total hip BMD in treated women 60–79 and 80 years or older

The analysis of T-score and percent change in BMD over time in women who started osteoporosis treatment was performed on total hip BMD because this measurement was found superior to lumbar spine for predicting fracture risk in a study of risk factors performed on the same dataset [14]. Women who had started an osteoporosis treatment1 and undergone at least two total hip BMD measurements were followed for up to 8 years from the dispensing of the first prescription (index). The first BMD measurement should have taken place no longer than 12 months before or after treatment start. To capture the real-world treatment effect on BMD, patients who discontinued treatment were not removed from observation at the time of discontinuation, and switching to other treatments during follow-up was allowed. Changes in BMD over time was also analysed in a subset of patients that were still persistent to treatment at the time of the DXA scan. Patients were defined as non-persistent if they did not fill a new prescription within 60 days from when their previous dispensing should have been consumed. Consumption of medication surplus was allowed. Follow-up was dived into 12-month windows and each patient had a baseline BMD at t0 and between one and six repeat measurements. Total hip BMD was measured on Lunar DPX and Hologic Discovery. Data was standardized and converted using the International DXA Standardization Committee (IDSC) [15].

Percent change in total hip BMD was used as outcome, and each follow-up year containing a BMD measurement was used in the estimation. Patient characteristics were measured during 24 months prior to index treatment. An unadjusted population averaged generalized linear panel data model was used to estimate the change per year while accounting for contribution of data by different patients at different time points and correlation of within-subject measurements.

Analysis 2: effectiveness in women 60–79 vs. 80+ using patients as their own control

Observational studies have examined the effectiveness of osteoporosis treatments for reducing clinical fractures. The designs of these observational studies include comparisons between patients with or without bisphosphonate use [16, 17] and compliant or not compliant with bisphosphonate use [18-20]. A key limitation in interpreting any of these comparisons is uncertainty whether known or unknown differences in baseline fracture risk between patient populations account for some of the estimated effect. Abelson and colleagues [10] proposed a method where the first 3 months after treatment start (index time point) serves as a proxy for the untreated baseline risk, thereby allowing patients to serve as their own control. Using this method, we estimated change in fracture incidence as incidence rate ratios (IRRs) in the oldest old compared with younger women using hip, vertebral, and non-hip-non-vertebral (NHNV) as fracture outcomes. Included patients had started a treatment during July 2007–August 2013 after a treatment free period of at least 24 months and with at least 4 months of follow-up. Patients were required to be persistent to treatment during their entire included follow-up, but switching between osteoporosis treatments during follow-up was allowed. Follow-up was limited to 15 months (3 + 12) to avoid capturing time-dependent changes not caused by the treatment. Patient characteristics were measured during 24 months before index. An unadjusted multiple failure Poisson survival model was chosen for fracture outcomes to avoid that women who fracture during the first 3 months were removed from the later period (4–15 months after treatment).

Analysis 3: effectiveness of osteoporosis treatment vs. calcium/vitamin D therapy after clinical fracture in women 80 years or older

Resting on the hypothesis that in women 80 years or older the group-level treatment allocation after fracture is mainly random, or dependent on other factors not directly associated with fracture risk, a comparison of fracture risk and mortality was made between treatment-naive women who had received at least one prescription of osteoporosis treatment after fracture (intention-to-treat, ITT group) and treatment-naïve women who instead had received therapy with calcium/vitamin D (> 73,200 IU during the first 6 months of treatment). Treatment with osteoporosis medication or calcium/vitamin D had to start between July 2007 and August 2013 and within 6 months after the fracture. Follow-up was limited to 24 months after treatment initiation (the index time point). Switching to other treatments was allowed after 6 months of follow-up in the ITT group and calcium/vitamin D group. Risk of any fracture was compared between treatment cohorts using an adjusted standard Cox model and an adjusted Fine-Gray competing risk model, accounting for death as a competing event. The sample did not allow analysis of separate outcome fracture types.

Results

Patient characteristics

Patient characteristics (Table 1) for the cohorts and subgroups studied in analyses 1 and 2 indicated that there were differences in comorbidities and concomitant medication between the age groups, as expected. Further, prior hip and vertebral fractures were more frequent, and BMD was lower in older women. Patient characteristics for the women in analysis 3 show that prior vertebral fracture was clearly relatively more frequent than other fracture types among women treated with osteoporosis medication than those who received calcium/vitamin D. There was also a pattern that women who had received calcium/vitamin D were slightly older, had more comorbidities, and were more than twice as likely to receive pre-packaged medication, which is indicative of special living conditions (e.g. nursing homes or extra care housing), home help, or other forms for assistance needs (Tables 2 and 3).

Table 1

Patient characteristics (measured during 24 months look-back from index, unless otherwise specified)

	Analysis 1: Total hip BMD over time in patients starting treatment		Analysis 2: Fracture incidence in patients starting treatment used as own control (0–3 months vs. 4–15 months)		Analysis 3: Fracture incidence after starting osteoporosis treatment vs. calcium/vitamin D after fracture in women 80+
	60–79 years	80+	60–79 years	80+	Osteoporosis treatment	Calcium/vitamin D
n	1907	254	81,039	35,542	5272	20,533
Follow-up in years (mean,( min-max))	3.5 (0.5–8.4)	2.9 (0.6–6.9)	0.8 (0.3–1.3)	0.8 (0.3–1.3)	1.6 (0.0–2.0)	1.5 (0.0–2.0)
Age at index (mean)	69.0	83.0	70.3	84.4	84.4	86.7
Prior fracture, any (%)	25.1	29.9	25.0	32.8	100.0	100.0
Hip (%)	3.9	11.0	4.7	11.1	31.3	43.2
Vertebral (%)	4.1	5.5	3.9	8.9	32.4	14.0
Non-hip, non-vertebral (%)	17.9	15.4	17.5	16.1	41.5	45.3
Previous OP treatment (%)	0.0	0.0	0.0	0.0	0.0	0.0
Diagnose related to secondary osteoporosis (%)	13.3	15.4	15.0	15.2	18.0	19.5
Systemic glucocorticoid use (%)	18.4	14.6	33.1	31.5	17.0	15.8
Charlson-Quan comorbidity index = 0 (%)	68.3	53.9	65.2	58.2	50.4	42.4
Charlson-Quan comorbidity index = 1–2 (%)	25.9	33.1	27.3	32.9	39.0	44.3
Charlson-Quan comorbidity index > 2 (%)	5.8	13.0	7.5	8.9	10.6	13.3
≥6 months calcium/vitamin D during last 24 months (%)	33.2	37.8	30.1	36.1	31.1	63.0
Treatment with proton pump inhibitors (%)	30.8	34.3	38.3	39.9	36.7	36.8
Pre-packaged drug dispensing (ApoDos) (%)	1.2	8.3	2.9	14.0	17.3	38.8
Post-index treatment1 (%)
Alendronate	82.9	82.3	88.7	88.0	88.6	3.2
Risedronate	12.5	6.3	7.9	6.1	5.5	0.2
Denosumab	4.7	11.4	2.3	3.8	4.3	0.4
Zoledronate	6.6	12.2	2.5	2.5	2.9	0.3
Raloxifene	0.6	0.0	0.6	0.3	0.4	0.0
Teriparatide	2.8	2.4	0.4	0.2	0.3	0.0
Number of hip DXA scans during follow-up (mean)	2.4	2.2	N/A	N/A	N/A	N/A
Baseline total hip BMD (mean g/cm3 (SD))2	0.56 (0.17)	0.50 (0.17)	N/A	N/A	N/A	N/A
Baseline total hip T-score (SD)2	− 2.06 (0.86)	− 2.45 (0.92)	N/A	N/A	N/A	N/A

1Switching allowed during follow-up

2Within 12 months before or after treatment start

Table 2

Total hip T-score (SD) and percent change from baseline in total hip BMD (g/cm3), all patients with ≥ 2 DXA measurements

	60–79 years				80+ years
Years from1st DXA	Patients with measurement	Mean T-score (SD) [CI95]	Percent change in TH BMD [CI95]	% treatment persistent	Patients with measurement	Mean T-score (SD) [CI95]	Percent change in TH BMD [CI95]	% treatment persistent
Baseline	1,907	−2.06 [−2.10, 2.02]			254	−2.45 [−2.56, −2.34]
1st	53	−2.01 [−2.06, 1.96]	3.2 [1.1, 5.3]	81%	8	−2.48 [−2.76, −2.21]	2.0 [−4.2, 8.1]	75%
2nd	453	−1.90 [−1.92, 1.88]	5.7 [5.0, 6.5]	73%	69	−2.38 [−2.48, −2.28]	3.3 [1.0, 5.6]	59%
3rd	810	−1.87 [−1.88, 1.85]	6.7 [6.1, 7.3]	66%	106	−2.29 [−2.37, −2.21]	5.9 [3.9, 7.8]	78%
4th	426	−1.88 [−1.90, 1.86]	6.0 [5.2, 6.8]	53%	57	−2.28 [−2.40, −2.17]	5.5 [3.0, 8.0]	60%
5th	356	−1.89 [−1.91, 1.87]	5.3 [4.5, 6.2]	42%	25	−2.27 [−2.44, −2.10]	6.5 [3.0, 10]	48%
6th	264	−1.90 [−1.93, 1.88]	5.3 [4.3, 6.3]	35%	19	−2.28 [−2.47, −2.08]	7.7 [3.7, 12]	47%
7th	144	−1.86 [−1.89, 1.83]	6.7 [5.4, 8.0]	37%	10	−2.21 [−2.46, −1.95]	6.3 [0.8, 12]	50%
8th	59	−1.92 [−1.97, 1.88]	5.0 [3.0, 6.9]	12%	0	–	–	–

Table 3

Total hip T-score (SD) and percent change from baseline in total hip BMD (g/cm3), persistent patients with ≥ 2 DXA measurements

Years from 1st DXA	60–79		80+
	Patients with measurement	Percent change in TH BMD [CI95]	Patients with measurement	Percent change in TH BMD [CI95]
Baseline	1152		166
1st	43	3.9 [1.8, 5.9]	6	1.0 [−5.1, 7.0]
2nd	330	7.5 [6.7, 8.3]	41	4.5 [1.9, 7.0]
3rd	531	8.9 [8.2, 9.6]	83	7.1 [5.2, 9.1]
4th	227	9.1 [8.1, 10]	34	8.6 [5.8, 11.3]
5th	148	7.9 [6.8, 9.1]	12	7.2 [2.8, 11.5]
6th	93	9.5 [8.1, 10.9]	9	18.8 [13.8, 23.7]
7th	53	9.7 [7.9, 11.6]	5	19.5 [12.9, 26.2]
8th	7	21.1 [16, 26.1]	–	–

Analysis 1: TH BMD over time in patients starting treatment

Repeated DXA scans of total hip (TH) BMD were available for 2161 women who had started an osteoporosis treatment. The majority (74%) had only two measurements, 21% had three, and 5% underwent ≥ 4 DXA scans, and the most frequent time point for a repeat measurement was during the third year after the baseline DXA. The linear percent change in TH BMD was estimated at + 1.13% units [1.03–1.22] per year, and there was no statistically significant difference between age groups (p = 0.309). The results did not change when running a fully adjusted model (+ 1.05% units [CI95 0.93–1.17] per year). Percent change in TH BMD was also associated (p < 0.05) with a baseline TH T-score ≤ − 2.5 SD (+ 0.6% units per year) and use of glucocorticoids (− 0.4% per year). The use of proton pump inhibitors (PPIs) was not significantly associated with percent change per year in TH BMD.

When only including women persistent to treatment during follow-up, the TH BMD increase was higher, with percent increase after 3 years reaching 8.9% and 7.1% in women 60–79 and 80+, respectively. This is to be compared with 6.7 and 5.9% units in the whole population with repeated DXA scans. The adjusted linear percent change in TH BMD was estimated at + 1.95% units [CI95 1.79–2.11] per year, and there was no statistically significant difference between age groups (p = 0.636). TH BMD increases of ~ 20% units were observed in women persistent to treatment for 6–8 years in both age groups, but these results should be interpreted with caution given the small patient numbers and risk of selection bias.

Analysis 2: fracture incidence in treated patients used as own control (0–3 months vs. 4–15 months)

Analysis of fracture incidence during 0–3 months after treatment start as a proxy for untreated fracture risk (Table 4) indicated that IRRs for change in fracture risk during treatment were similar across age groups but different across fracture types. There were no statistically significant differences between the age groups.

Table 4

Incidence during 0–3 months and 4–15 months after treatment start and corresponding IRRs

Outcome and age group	No. of fractures	Incidence rate (%) 0–3 months [CI95]	Incidence rate (%) 4–15 months [CI95]	IRR [CI95]*
Any fracture
60–79 years	2992	6.08 [5.74, 6.44]	3.83 [3.66, 4.01]	0.63 [0.58, 0.68]
80+ years	2219	10.41 [9.74, 11.13]	7.18 [6.81, 7.57]	0.69 [0.63, 0.75]
All ages	5211	7.39 [7.07, 7.72]	4.78 [4.62, 4.95]	0.65 [0.61, 0.68]
Hip fracture
60–79 years	474	0.90 [0.78, 1.05]	0.62 [0.56, 0.70]	0.69 [0.57, 0.84]
80+ years	635	2.67 [2.34, 3.05]	2.16 [1.96, 2.37]	0.81 [0.68, 0.95]
All ages	1109	1.43 [1.30, 1.58]	1.06 [0.99, 1.14]	0.74 [0.65, 0.84]
Vertebral fracture
60–79 years	583	1.84 [1.66, 2.04]	0.48 [0.42, 0.54]	0.26 [0.22, 0.31]
80+ years	431	2.92 [2.58, 3.32]	0.98 [0.85, 1.14]	0.34 [0.28, 0.41]
All ages	1014	2.17 [2.00, 2.35]	0.62 [0.57, 0.69]	0.29 [0.25, 0.33]
NHNV** fracture
60–79 years	1935	3.34 [3.09, 3.61]	2.70 [2.56, 2.85]	0.81 [0.73, 0.89]
80+ years	1153	4.82 [4.37, 5.32]	3.94 [3.67, 4.24]	0.82 [0.72, 0.93]
All ages	3088	3.79 [3.56, 4.02]	3.06 [2.93, 3.19]	0.81 [0.75, 0.87]

*Incidence rate ratios < 1 means that incidence was higher during the first 3 months after treatment initiation

**Non-hip/non-vertebral

Analysis 3: fracture incidence and mortality associated with osteoporosis treatment vs. calcium/vitamin D after clinical fracture in women 80 years or older

Comparing the two groups (Fig. 1) showed that patients receiving osteoporosis treatment (ITT group) had lower mortality and lower fracture risk than patients who started calcium/vitamin D. The adjusted Cox model with any fracture as the outcome estimated a HR of 0.77 [0.71–0.83], and the competing risk model estimated a HR of 0.78 [0.72–0.85], indicating that the incorporation of mortality as a competing risk should not change the interpretation. The lower crude mortality (HR 0.51 [0.47–0.55]) in patients receiving osteoporosis treatment indicated that treatment allocation after fracture in the age 80+ strata was not random. Excluding patients with hip fracture as their index fracture did not change this interpretation (HR 0.58).

Fig. 1

Incidence of any fracture (a) and mortality (b) in treatment-naïve women 80 years or older receiving osteoporosis treatment or calcium/vitamin D after a clinical fracture

Discussion

Pharmaceutical prevention of secondary fractures in the women older than 80 years is rare when compared with younger women, and there are likely several reasons for this. Uncertainty regarding treatment effectiveness, concerns about polypharmacy, frailty and limited life expectancy, views on natural ageing, and reduced patient self-advocacy in the elderly are factors that partly may explain this discrepancy. This study aimed at investigating the real-world effectiveness of osteoporosis treatment in the oldest old women (80+) compared with younger women (60–79 years).

The use of non-randomized retrospective data to estimate effectiveness can never replace randomized trials, because the analysis will always be limited by some degree of bias. However, real-world effectiveness can provide an important complement in that large samples efficiently can be acquired to allow for the study of comparative effectiveness, subgroups, interactions, or rare safety events. Possibly even more important is that studies can be performed for the actual context and population in which the treatments are used in clinical practice. We choose to analyse treatment of the oldest old vs. younger women from several aspects to allow them to complement and validate each other.

Analysis 1: percent change in BMD in treated women 60–79 and 80 years or older

Percent change in TH BMD was largely the same in both age groups, and the mean percent change in TH BMD over time in treated patients was very similar to what has been reported in an RCT for up to 10 years by Bone and colleagues [21]. Apart from the lack of randomization and a placebo control, the approaches differed in that our analysis was performed on a sample with repeated BMD measurements in patients starting treatment (≥ 2 DXA scans). A natural consequence of this is that follow-up BMD was estimated in different individuals at different time points. There is a risk of bias in that the number and timing of BMD measurements may be associated with the measurement results, but no such clear pattern was observed. It was beyond the scope of this study to analyse BMD development with individual treatments separately, but it may be worthwhile to revisit this in future research.

BMD change in patients fully persistent to treatment during follow-up was generally higher, as expected. Like for the entire BMD population, there was no statistically significant difference of the linear BMD change per year between women 60–79 years and 80+ years. In agreement with previous studies of persistence to osteoporosis treatments and outcomes [22, 23], these results further reinforce the importance of keeping patients on treatment for as long as it is clinically motivated.

Analysis 2: incidence in treated patients used as own control (0–3 months vs. 4–15 months)

Using a method originally proposed by Abelson [10], we could show that fracture incidences were higher during the first 3 months after treatment initiation compared with the following 12 months. This pattern was consistent across both age groups and for all three fracture type outcomes analysed. IRRs were similar between age groups with IRRs (all ages combined) of 0.65, 0.74, 0.34, and 0.81 for any fracture, hip fracture, vertebral fracture, and NHNV, respectively. These estimates are comparable with the results reported by Abelson et al. [10] as well as estimates from the RCT meta-analysis of bisphosphonate conducted by NICE, which reported HRs of 0.85, 0.45, and 0.79 for hip fracture, vertebral fracture, and non-vertebral fracture, respectively [24]. However, such comparisons alone cannot validate the method, which rests on the assumption that the first 3 months of treatment are representative of a patient’s untreated risk [25] and that treatment is fully effective during the following 12 months. The approach is useful since it allows for the comparison of groups with different characteristics by using patients as their own controls. More typical comparisons when groups or products are compared with each other are often subject to indication bias, where the reasons for an individual’s group allocation are associated with the outcome variable. Such bias can be addressed by using propensity scores or multivariate adjustment, but unobserved confounding may make the interpretation of results challenging [26].

Analysis 3: fracture incidence and mortality associated with osteoporosis treatment vs. calcium/vitamin D after clinical fracture in women 80 years or older

This analysis clearly showed that treatment choice after fracture in the oldest old is not random. This despite that fracture risk assessment with a BMD measurement and/or treatment is rare in this group, and most Swedish patients sustaining a fracture after the age of 80 years will have a fracture risk warranting treatment, irrespective of the presence of other risk factors. For example, the 10-year probability of a major osteoporotic fracture in a Swedish 80-year-old woman with a previous fracture, normal BMI, and no other risk factors is estimated at 39% by the FRAX algorithm (https://www.sheffield.ac.uk/FRAX). Moreover, our previous research on the impact of risk factors in the oldest old indicated that factors, other than total hip BMD and previous fracture, added limited value to the fracture risk assessment (NB: data on fall risk, smoking and heredity was not available) [14]. Calcium/vitamin D therapy is not recommended for fracture prevention in Swedish patients < 80 years but has been shown in a meta-analysis to confer moderate risk reductions in patients > 80 years, with pooled RRs ranging from 0.95 to 0.84, depending on fracture type [27]. Besides being the standard comparator in placebo-controlled RCTs in osteoporosis, calcium/vitamin D was chosen as a comparator because these patients also have been assessed to need some type of treatment after their fracture. Nonetheless, some patient characteristics and mortality were clearly different when comparing patients receiving alendronate and calcium/vitamin D after fracture. Patients receiving only calcium/vitamin D had more often had a hip fracture and less often a vertebral fracture as the index fracture, had more comorbidities, and were more than twice as likely to receive pre-packaged drug dispensing, which indicates a higher need for assistance, like home care or different types special living accommodations. Mortality was twice as high in patients receiving calcium/vitamin D, which to some extent may be a result of avoided fractures but is still suggesting that treatment decisions in this population may focus more on morbidity than on fracture risk. Omitting hip fracture patients in a sensitivity analysis did not substantially change this interpretation (Mortality HR 0.58), indicating that the differences in mortality only partly can be explained by the observed differences in the distribution of index fracture types. In the adjusted multivariate fracture analysis, osteoporosis medication was associated with a 23% lower risk of any fracture compared with calcium/vitamin D, but this estimate should be interpreted with caution given the differences in patient profiles. The exact same risk reduction was estimated by Axelsson et al. [11] who studied alendronate vs. no treatment in a similar fashion and population but with untreated matched controls instead of calcium/vitamin D. The more balanced groups achieved by Axelsson and colleagues with matching reduced the mortality difference between the groups (HR 0.88) but instead limited the possibility to study how these populations can be characterized in the clinical setting. Irrespective of whether the approach employed in the present study is useful for measuring effectiveness, it is noteworthy how many patients > 80 years are treated only with calcium/vitamin D after fracture but who possibly would benefit from osteoporosis treatment (~ 4 times the osteoporosis treatment group size).

Limitations

When treatment effectiveness is studied in the real-world setting, it is by design intended that patients should be treated without strict criteria or frequent visits, adherence is not compulsory, and switching between treatments occurs. Apart from this, the analyses are limited by several aspects in the context of drawing inference regarding effectiveness. The absence of randomization will introduce indication bias, where unobserved characteristics that are unevenly distributed between study groups may be associated with the outcome. Furthermore, zoledronate has over time more often been distributed directly via the hospitals instead of pharmacies. Because only drugs dispensed at pharmacies are captured in the Prescribed Drugs Register, it is likely that some patients who appear to have been prescribed only calcium/vitamin D also were prescribed zoledronate during the follow-up period, which should be factored into the interpretation of the association with fracture risk found in analysis 3.

Conclusions

Our results indicate that osteoporosis medication works in women older than 80 years and that the magnitude of effect is similar to what we estimated in younger women and to what has been reported in randomized trials. Using patients as their own control appears to be a useful method for eliminating indication bias and measuring treatment effect in osteoporosis. The choice between osteoporosis treatment and calcium/vitamin D after fracture in Swedish women ≥ 80 years is not random but appears associated with the patient’s health status and presence of vertebral fractures, rather than the known risk profile of sustaining a fracture at a high age.

(DOCX 16 kb)