Does treatment with bisphosphonates protect against fractures in real life? The HUNT study, Norway.

Bisphosphonates reduce fractures in randomized controlled trials (RCT); however, there is less information from real life. In our population including 14,990 women and 13,239 men, use of bisphosphonates reduced risk of fractures in hip and forearm in women. The magnitude of the effect was comparable to results from RCT.
INTRODUCTION: The objective was to examine if treatment with bisphosphonates (BPs) was associated with reduced risk of fractures in the hip and forearm in women and men in the general population.
METHODS: In a cohort study based on data from the third wave of the population-based HUNT Study (HUNT3), the fracture registry in Nord-Trøndelag, and the Norwegian Prescription Database, 14,990 women and 13,239 men 50-85 years were followed from the date of participating in HUNT3 (2006-2008) until the date of first fracture in the hip or forearm, death, or end of study (31 December 2012). Hazard ratios with 95% confidence intervals for hip and forearm fracture according to use of BPs were estimated using Cox proportional hazards models with time-dependent exposure. Adjustment for individual FRAX® fracture risk assessment scores was included.
RESULTS: BPs, predominantly alendronate, were used by 9.4% of the women and 1.5% of the men. During a median of 5.2 years of follow-up, 265 women and 133 men had a hip fracture, and 662 women and 127 men had a forearm fracture. Compared with non-users of BPs, the hazard ratios with 95% confidence interval for a fracture among users of BPs adjusted for age and FRAX® were 0.67 (0.52-0.86) for women and 1.13 (0.50-2.57) for men. Among users of glucocorticoids, the corresponding figures were 0.35 (0.19-0.66) and 1.16 (0.33-4.09), respectively.
CONCLUSIONS: Use of BPs was associated with reduced risk of fractures in hip and forearm in women, and the magnitude of effect is comparable to results from RCTs.

Introduction

Osteoporosis is a skeletal disease characterized by low-bone mineral density (BMD) and microarchitectural deterioration of bone tissue, leading to bone fragility and increased fracture risk [1]. Over the last decades, several anti-osteoporotic drugs (AODs) have emerged, including inhibitors of bone resorption such as bisphosphonates (BPs), raloxifene, and denosumab, as well as teriparatide, romosozumab, and abaloparatide which stimulate bone formation [2, 3].

Treatment with AODs has in randomized clinical trials (RCTs) been found to reduce the relative fracture risk for vertebral and non-vertebral fractures by 40–80% and 20–60%, respectively [4-6]. Oral BPs, preferentially alendronate, are the most widely used worldwide because they are effective and inexpensive [3]. Most RCTs have included women with postmenopausal osteoporosis, excluding cases with secondary osteoporosis and those with several diagnoses. In contrast, patients treated in clinical practice are more heterogenous and include frail individuals with comorbidities [7, 8].

Reyes et al. observed that half of the patients who received alendronate treatment had secondary or severe osteoporosis [9]. A few RCTs have addressed secondary osteoporosis, preferentially glucocorticoid (GC)-induced osteoporosis (GIOP). A meta-analysis including 10 studies and 1002 patients with GIOP showed that alendronate increased BMD of the lumbar spine and femoral neck without a corresponding effect on fracture risk [10]. In contrast, a Cochrane review concluded that BPs reduced vertebral fractures and prevented bone loss both in the lumbar spine and femoral neck. An effect on non-vertebral fractures was not found; however, according to the authors the effect of non-vertebral fracture was of low-certainty evidence [11].

In real life, adherence to oral BPs is generally low compared to that in RCTs [7, 8], and premature cessation of treatment is particularly common in men [12]. Discontinuation of BPs may be due to polypharmacy or side effects [3]. The adherence is also dependent on the route of administration of the drug. These factors may contribute to a reduced external validity of results from RCT.

There is little information regarding the effect of BPs in the general population. Danish and Swedish registry studies have reported 30% risk reduction in hip fractures in women aged 50–94 years [7, 13, 14]. Similarly, a large population-based study from the USA showed a 21% reduced risk of any clinical fracture in users of BPs, when comparing 1-year fracture incidence before and after initiating treatment [8].

Although the incidence of osteoporotic fractures in Norway is among the highest in the world, the prescription of BPs has been relatively low compared to some other European countries; a little less than 5% among all women > 40 years [15-18]. In the HUNT study, prescription of AODs (whereof BPs contributed to 98%) was low even among individuals with high fracture risk assessed by the Fracture Risk Assessment Tool, FRAX® (www.FRAX.com). FRAX® predicts the 10-year absolute risk of hip fracture as well as a major osteoporotic fracture. Only 25% of the women and 17% of the men with > 20% risk of a major osteoporotic fracture were treated with AODs [16]. Moreover, merely 11% of the women and 3% of the men used AODs the first year after a forearm fracture [15].

The objective of this study was to examine if use of BPs was associated with reduced risk of fractures in the hip and forearm in women and men aged 50–85 years in a Norwegian population.

Methods

Data sources

Participants from the third survey of the HUNT Study (HUNT3) comprised the study population. Their data were linked to the fracture registry of Nord-Trøndelag and the Norwegian Prescription Database (NorPD) via the national personal identification number.

The Trøndelag Health Study (HUNT study)

HUNT3 was performed from 2006 to 2008 in the county of Nord-Trøndelag, located in mid-Norway. The geographic, demographic, and occupational structure is considered fairly representative of the country as a whole [19]. All individuals aged 20 years and older residing in the county were invited to participate.

In the current study, we included the age group 50–85 years. Of the 43,760 invited, 28,692 (65.6%) responded, completed a comprehensive questionnaire, and underwent a short clinical examination at the screening station. Of these, 231 were excluded due to lack of data on height (N = 213) or weight (N = 220) and 26 used other AODs than BPs, leaving 14,990 women and 13,239 men (Fig. 1).

Fig. 1

Flowchart of the included subjects

The fracture registry of Nord-Trøndelag

The registry provides validated information on forearm and hip fracture in subjects older than 16 years, treated, or followed up from 15 August 1995 to 31 December 2012 at the two hospitals located in Nord-Trøndelag County. Data were retrieved through the electronic discharge registers, namely The Patient Administrative System (PAS) based on diagnoses coded according to the International Classification of Diseases (ICD), as well as surgical procedures coded according to the NOMESCO Classification of Surgical Procedures (NCSP). A fracture was defined when (1) the ICD code was accompanied by a medical record confirmation of hip fracture or (2) a fracture was diagnosed by X-ray.

Fractures due to metastatic disease were excluded. Details about the classification and validity of this fracture information have been published previously [15, 20, 21].

Norwegian Prescription Database (NorPD)

Data on prescriptions of BPs and oral GCs were collected from the NorPD established 01 January 2004. NorPD contains information on all prescribed drugs that are dispensed at all pharmacies in Norway to individual patients in ambulatory care. Drugs prescribed to patients in hospital or other institutions, are not registered in NorPD [22]. Month of death was retrieved form the NorPD. Information regarding other types of AODs, including denosumab (M05BX04), raloxifene (G03XC01), and teriparatide (H05AA02), were also collected from the NorPD as use of these was an exclusion criterion from the study.

Exposure: bisphosphonates

For each individual, exposure to BPs, based on filled prescriptions with Anatomical Therapeutic Chemical (ATC) code M05BA with all subgroups, was registered from 1 year before participation in HUNT3 and until either fracture, death, or end of study (31 December 2012).

Time-dependent exposure of BPs was applied to estimate the association between such use and risk of fractures. Individuals who filled at least two prescriptions for bisphosphonates within a 6-month period at any time during follow-up were defined as exposed from the date of filling the first prescription until 6 months after estimated end of drug supply. An exception from this was use of zoledronic acid where the inclusion criterion was one infusion per year.

Outcome: fracture

Data on first forearm or hip fracture after the HUNT3 examination until study end were obtained from the fracture registry of Nord-Trøndelag. The ICD 10 codes included to identify hip fractures were S72.0–2 and 9, and for forearm fractures S52.0–S52.9.

Covariates

Fracture risk was calculated by use of the FRAX® (www.shef.ac.uk/FRAX) (for Norway) [23]. Individual FRAX scores in terms of predicted 10-year risk of hip fracture (FRAXHIP) and major osteoporotic fracture (FRAXMOF) without BMD were calculated based on information from the HUNT study and the NorPD; details are previously described [24].

The following variables were included in the FRAX calculation: gender, age, body mass index (BMI), self-reported previous fracture, parental hip fracture, current smoking, use of alcohol (units per week), use of oral GCs, rheumatoid arthritis (RA), and other causes of secondary osteoporosis, including diabetes type 1, premature menopause, and hyperthyroidism. We did not have any information regarding malabsorption, chronic liver disease, or osteogenesis imperfecta. As the fracture register was established in 1995 and only contains fractures in the hip and wrist, self-reported fractures at the hip, spine, or wrist after the age of 40 years in HUNT3 were included.

There were no self-reported data on use of oral GCs in HUNT3, so this information was retrieved from the NorPD (ATC code H02A and H02B). In line with the guidelines of FRAX, use of at least 5 mg GCs for more than 3 months prior to the inclusion in HUNT3 or current use were included in the FRAX calculation [24] (www.FRAX.com).

All other information regarding FRAX was collected from HUNT3 [24].

Statistics

Descriptive data are given as means with standard deviations (SD) for continuous data and numbers and percentages for categorical data, and 95% confidence intervals (CIs) for counts were calculated by the continuity-corrected score interval method [25].

Cox proportional hazards model with time-dependent exposure to BPs was used to assess the association between BPs and fracture rates. Here, for each time increment (months) an individual would be classified as unexposed or exposed according to the description in the “Exposure: bisphosphonates” section. Thus, an individual may contribute with unexposed time during the first part of follow-up, switch to exposed when fulfilling the exposure criteria, and possibly become unexposed at a later point in time. The study population was followed from the date of participating in HUNT3 (the index date) until date of first fracture, death, or end of study (31 December 2012), whichever came first.

We present three models: the first adjusted for age, the second adjusted for FRAXMOF without BMD, and the third adjusted for both age and FRAXMOF without BMD for analyses of the composite endpoint of hip- and forearm fracture. FRAXHIP without BMD was used in the models with hip fracture as outcome.

The proportional hazards assumptions were checked by visual inspection of log minus log plots.

For descriptive purposes medication possession ratio (MPR) was calculated as a measure of compliance. MPR is defined as the proportion of a time period where medication supply is available and was calculated as the total amount defined daily dose (DDD) prescribed/(last fill date − first fill date + DDD of the last prescription). Adherence was defined as MPR ≥ 80% [26, 27]. The estimate for MPR is only descriptive and not included in our statistical models.

In our main analyses, we assessed the composite endpoint of first hip or forearm fracture in women and men separately. In addition, we also assessed hip and forearm fractures separately.

Sensitivity analyses

The impact of BPs was assessed in the following subgroups among women: secondary osteoporosis, those using GCs, and those with high risk for fractures classified as having a FRAXMOF score ≥ 20%.

As the effect of BPs is dependent on the duration of use [3], we also applied models with different assumptions on duration of exposure to BPs after the last filled prescription. In addition to our primary analysis, assuming that the effect on bone would last for 6 months after termination of drug supply, we also performed separate analyses replacing this assumption with 3 and 9 months, as well as censoring 6 months after stopping BP treatment.

In additional sensitivity analyses, we estimated propensity scores for BP exposure by logistic regression including age, BMI, previous fractures, osteoporosis among parents, smoking, daily milk intake, alcohol consumption, exposure to glucocorticoids, rheumatoid arthritis, premature menopause, hyperthyroidism, and diabetes type 1 as covariates. We compared the distribution of propensity scores between users and never users of PPIs and performed a stratified analysis using 6 strata as well as separate analyses in each stratum (women) to assess whether the association seemed to differ between strata.

All analyses were done using STATA 14.1 (StataCorp. 2015. Stata Statistical Software: Release 14. College Station, TX: StataCorp LP).

Ethics

Participants in HUNT3 gave written informed consent for use of their data in research including linkage to named registries, such as NorPD. The study was approved by the Regional Committee for Medical and Health Research Ethics in Central Norway (2012/1906/REK). Linkage of databases was approved by the Norwegian Data Protection Authority.

Results

Use of bisphosphonates

The study includes 14,990 women and 13,239 men (Fig. 1). Alendronate was by far the most prescribed BPs, used by 98% of the women and 99% of the men. Only 26 individuals were excluded in the analyses due to use of other AODs (Fig. 1 and Supplementary Table 1).

A total of 1402 women (9.4%) and 201 men (1.5%) met the criteria of filling at least two prescriptions of peroral BPs during a 6-month period or infusion of zoledronic acid once a year and were included in the analysis. Four women and one man had been treated with zoledronic acid as the only BPs, while 57 women and 7 men used zoledronic acid after side effects from peroral BPs (Supplementary Table 1). We classified 251 women and 72 men who filled only one prescription for peroral BPs as non-users.

Those treated with BPs were substantially older than those who were not, with a mean age at baseline of 70.2 and 68.6 years for women and men, respectively. Among those treated with BPs, women reported more previous fractures (44.5% vs. 26.9%) and osteoporosis among parents than men (29.8% vs. 16.9%). A larger proportion of men than women using BPs had secondary osteoporosis (63.7% vs. 50.4%), and of these, 120 men and 441 women were treated with GCs (59.7% vs. 31.5%) (Table 1).

Table 1

Baseline values stratified for use of bisphosphonates1 in the follow-up period

	Women, N = 14,990		Men, N = 13,239
	Exposed to bisphosphonatesN = 1402	Not exposed to bisphosphonatesN = 13,588	Exposed to bisphosphonatesN = 201	Not exposed to bisphosphonatesN = 13,038
Age, years mean (SD)	70.2 (8.7)	63.2 (8.9)	68.6 (9.0)	63.6 (898)
BMI, mg/kg2 mean (SD)	26.3 (4.5)	27.7 (4.7)	26.7 (4.0)	27.7 (3.7)
FRAX® calculated 10-years risk of major osteoporotic fracture, mean (SD)	26.3 (15.6)	12.3 (9.5)	12.2 (8.9)	6.3 (3.7)
Previous fractures, N (%)b,c	624 (44.5)	1980 (14.6)	54 (26.9)	1276 (9.8)
Osteoporosis among parents, N (%)b	418 (29.8)	1513 (11.1)	34 (16.9)	795 (6.1)
Current smoker, N (%)b	307(21.9)	3204 (23.6)	40 (19.9)	2785 (21.4)
More than 2 glass milk/day, N (%)b	896 (63.9)	8190 (60.3)	84 (41.8)	4929 (37.8)
Alcohol consumptionb
≥ 2 units/day, N (%)	30 (2.2)	372 (2.8)	7 (3.5)	1121 (8.7)
≥ 3 units/day, N (%)a	0	4 (0)	0	28 (0.2)
Secondary osteoporosis, N (%)d	706 (50.4)	3656 (26.9)	128 (63.7)	1951 (15.0)
Exposed to glucocorticoids, N (%)e	441 (31.5)	1587(11.7)	120 (59.7)	1417(10.9)
Rheumatoid arthritis, N (%)b	188 (13.4)	672 (4.9)	33 (16.4)	375 (2.9)
Premature menopause, N (%)b	259 (8.5)	1591 (11.7)	–	–
Hyperthyroidism, N (%)b	57 (4.1)	400 (2.9)	1 (0.5)	127 (1.0)
Diabetes type 1, N (%)b	15 (1.1)	62 (0.5)	4 (2.0)	86 (0.7)

aUse of bisphosphonates were defined as at least two prescriptions of bisphosphonates were filled during a 6 months period or infusion of zoledronic acid once a year

bSelf-reported

cPrevious fractures in hip, wrist, or spine after 40 years

dSecondary osteoporosis defined as at least one of the following: exposed to glucocorticoids (min 90 DDD or 3 months treatment), rheumatoid arthritis, menopause or surgical removal of ovaries before 45 years, hyperthyroidism, or diabetes mellitus type 1. The number of patients with of secondary osteoporosis is less than the sum of the different conditions, as 602 women and 145 men had more than one criterion for secondary osteoporosis

eUsed the drug during the follow-up period, data were retrieved from the Norwegian prescription register

The mean duration of use of BPs was 52.2 months (SD 27.8) among women and 42.9 months (SD 26.9) for men. Regarding adherence among the included participants with at least two filled prescriptions, 79.8% of the women and 77.1% of the men had an MPR ≥ 80 during follow-up.

Fracture risk according to use of bisphosphonates

During a median of 5.2 years follow-up, 265 women and 133 men had a first hip fracture, and 662 women and 127 men had a first forearm fracture. Among women with fractures, 94 were treated with BPs and 828 were not. The respective numbers for men were 9 and 250 (Table 2). Adjusted only for age, the HR for hip/forearm fracture among women using BPs was 0.90 (95% CI 0.71–1.14) (Table 2). However, after adjustment for individual FRAX scores, women had a HR for hip/ forearm fractures of 0.67 (95% CI 0.52–0.86) and 0.68 (0.46–1.02) for hip fractures solely. Only 9 men using BPs experienced fractures, and no difference in the fracture risk was shown between users and non-users (Table 2).

Table 2

Hazard ratios for the association between use of bisphosphonates and fractures in the 3 different models

	N	Fractures	Fractures/1000 individuals	Model 1, adjusted for age	Model 2, adjusted for FRAX	Model 3, adjusted for age and FRAX
Hip and forearm fractures
Women
Bisphosphonates no	13,588	828	60.9 (57.0–65.1)	Reference	Reference	Reference
Bisphosphonates yes	1402	94	67.0 (54.8–81.7)	0.90 (0.71–1.14)	0.64 (0.50–0.82)	0.67 (0.52–0.86)
Men
Bisphosphonates no	13,038	251	19.2 (16.9–21.7)	Reference	Reference	Reference
Bisphosphonates yes	201	9	44.7 (22.0–86.0)	2.00 (0.94–4.25)	0.80 (0.34–1.88)	1.13 (0.50–2.58)
Hip fractures
Women
Bisphosphonates no	13,588	230	16.9 (14.9–19.3)	Reference	Reference	Reference
Bisphosphonates yes	1402	35	25.0 (17.7–34.9)	0.87 (0.59–1.27)	0.64 (0.42–0.98)	0.68 (0.46–1.02)
Men
Bisphosphonates no	13,038	128	9.8 (8.2–11.7)	Reference	Reference	Reference
Bisphosphonates yes	201	5	24.9 (9.2–60.3)	1.33 (0.42–4.20)	0.31 (0.09–1.17)	0.78 (0.22–2.68)
Forearm fractures
Women
Bisphosphonates no	13,588	601	44.2 (40.9–47.8)	Reference	Reference	Reference
Bisphosphonates yes	1402	61	43.5 (33.7–55.9)	0.96 (0.72–1.29)	0.87 (0.64–1.18)	0.87 (0.64–1.18)
Men
Bisphosphonates no	13,038	123	9.4 (7.8–11.3)	Reference	Reference	Reference
Bisphosphonates yes	201	4	19.9 (6.4–53.5)	3.14 (1.15–8.55)	2.33 (0.80–6.81)	2.21 (0.73–6.63)

In women with secondary osteoporosis, the HR for fractures of the hip/forearm was 0.54 (95% CI 0.37–0.82), hip fractures 0.54 (95% CI 0.30–0.99), and forearm fractures 0.62 (95% CI 0.37–1.02) (Table 3a). In women using GCs (Table 3b), the HR was 0.35 (95% CI 0.19–0.66) for fractures in hip/forearm, 0.33 (95% CI 0.13–0.90) for hip fractures, and 0.45 (95% CI 0.21–0.98) for forearm fractures. When examining women with high risk for fracture defined as FRAX ≥ 20 (N = 3001), the HR for fracture in hip/forearm was 0.59 (0.44–0.80) (data not shown).

Table 3

Sensitivity analyses including (a) women with all types of secondary osteoporosis and (b) women and men using glucocorticoids

	N	Fractures	Fractures/1000 individuals	Model 1, adjusted for age	Model 2, adjusted for age and FRAX
a. Women with all types of secondary osteoporosis
Hip and forearm fractures
Bisphosphonates no	3656	211	57.7 (50.5–65.9)	Reference	Reference
Bisphosphonates yes	706	40	56.7 (41.3–77.0)	0.75 (0.51–1.10)	0.54 (0.37–0.82)
Hip fractures
Bisphosphonates no	3656	73	20.0 (15.7–25.2)	Reference	Reference
Bisphosphonates yes	706	17	24.1 (14.5–39.1)	0.75 (0.42–1.34)	0.54 (0.30–0.99)
Forearm fractures
Bisphosphonates no	3656	139	38.0 (32.1–44.9)	Reference	Reference
Bisphosphonates yes	706	25	35.4 (23.5–52.6)	0.81 (0.50–1.32)	0.62 (0.37–1.02)
b. Women exposed to glucocorticoids
Hip and forearm fractures
Bisphosphonates no	1587	81	51.0 (41.0–63.3)	Reference	Reference
Bisphosphonates yes	441	17	38.5 (23.3–62.2)	0.48 (0.26–0.88)	0.35 (0.19–0.66)
Hip fractures,
Bisphosphonates no	1587	33	20.8 (14.6--29.4)	Reference	Reference
Bisphosphonates yes	441	6	13.6 (5.5–30.9)	0.64 (0.42–0.98)	0.33 (0.13–0.90)
Forearm fractures, women
Bisphosphonates no	1587	48	30.2 (22.6–40.2)	Reference	Reference
Bisphosphonates yes	441	12	27.2 (14.8–48.4)	0.57 (0.27–1.21)	0.45 (0.21–0.98)

There were only minor differences when accounting for potential variation in assumed duration of the antiresorptive effect of BPs after end of supply. When duration of exposure was assumed to last 3 months after end of drug supply, the HR among women (adjusted for age and FRAXMOF without BMD) was 0.63 (0.49–0.81) When the exposure was assumed to last 9 months, the respective values were 0.68 (0.53–0.87). In the model censoring BPs users 6 months after stopping BPs treatment, the HR were 0.73 (0.58–0.93).

In propensity score stratified models, estimated associations between use of bisphosphonates and fractures were very similar to those estimated from the model adjusted for age and FRAX. For women, the HR among users of BPs was 0.75 (0.58–0.95) and men 1.38 (0.63–4.05).

Discussion

In this population-based study conducted between 2006 and 2012, the risk of fracture in the forearm or hip was 33% lower in women using BPs when taking into account individual fracture risk assessment (FRAX®) scores. This is in agreement with results from RCTs, and the fact that different sensitivity analyses broadly gave the same results seems reassuring. The association was most pronounced for hip fractures. A risk reduction could not be shown for men, probably due to lack of power as only 201 men were BP users with only nine fractures occurring among these during follow-up. Notably, 98% and 99% of women and men, respectively, were treated with oral BPs, preferentially alendronate. About 50% of the women using BPs had secondary osteoporosis and a third of these had GIOP. Interestingly, women with secondary osteoporosis seemed to benefit the most, exhibiting a relative risk reduction of 46%, and those with GIOP a reduction by 65%.

It is important to complement RCTs with observational studies addressing the effect of BPs in daily clinical practice. The clinical picture differs between those participating in RCTs and patients in real life, the latter exhibiting comorbidity and more severe osteoporosis [9]. The drug adherence also seems to be poorer in clinical practice. Accordingly, the anti-fracture efficacy seen in RCTs does not necessarily translate to real life.

A crude comparison of the fracture rates demonstrates that BP users had a higher fracture rate than non-users. This mirrors the appropriate targeting of BPs to individuals at substantially elevated fracture risk, with BP users being older, three times as likely to have sustained prior fractures and more likely to suffer from secondary osteoporosis. Hence, in the present study, women prescribed BPs had a risk for a major osteoporotic fracture of 26% according to FRAX (without measured BMD), a base risk more than twice that of untreated women in the study. When effects are adjusted for this difference, it becomes evident that the observed fracture rates represent a substantial reduction of risk in those who received BP treatment.

The estimated reduction of fracture risk was somewhat lower than in the initial RCTs of oral BPs. The results are not necessarily comparable, due to strict inclusion and exclusion criteria in RCTs and potential residual confounding in our observational study. In the pivotal fracture intervention trial (FIT) of alendronate including women with existing vertebral fractures, fractures at the hip and wrist were reduced by 51 and 48%, respectively [28]. In the subgroups of women with a very high risk for fracture defined as FRAX ≥ 20, the risk reduction was estimated to 41%. This emphasizes the importance of treating women at high risk. Women with high FRAX often have comorbidities that exclude them from participation in RCTs [29]. Our results are comparable to those found in studies from Denmark and Sweden [14, 30].

In total, among those treated with BPs, 50% of women and 63% of men had secondary osteoporosis, where GIOP accounted for 31% and 59%, respectively. Given potential differences in the effect of BPs on fracture reduction in primary and secondary osteoporosis, a subanalysis of these groups was performed. Among women with secondary osteoporosis, the fracture reduction in our study was 46%, and among users of GCs as high as 65%. There is limited evidence for an anti-fracture effect of BPs for most types of secondary osteoporosis, except for GIOP and osteoporosis induced by hypogonadism in men. Notably, most studies were not powered to assess fracture risk reduction. In an RCT addressing GIOP, a trend towards risk reduction of vertebral fractures was observed after 1 year of alendronate treatment. In the open-label extension study at 2 years, the risk of vertebral fracture was reduced by 90%. The study did not have sufficient power to demonstrate an effect on non-vertebral fractures [31, 32].

The reason for the more pronounced fracture reduction among GCs users may be attributed to the severity of their osteoporosis. GCs exert a multitude of adverse effects on bone and muscle contributing to a high risk for fractures. Due to impairment of bone quality, these patients experience fractures at a higher BMD than non-users. GCs compromise inhibition of bone formation with a concomitant stimulation of bone resorption including increased apoptosis of mature osteoblasts and osteocytes; impaired differentiation of osteoblasts; and an increased lifespan of osteoclasts [33]. GCs further inhibit the intestinal calcium absorption. In addition, the diseases treated by GCs may also increase fracture risk themselves through inflammatory mediators, inappropriate diet, and inactivity [33-35].

Oral BPs, and preferentially alendronate, are the recommended initial treatment for osteoporosis worldwide due to their efficacy, beneficial safety profile, and low cost. In line with this, alendronate was used by the majority of the study subjects, whereas only 69 individuals (4 % of the BP users) were treated with iv zoledronic acid. Zoledronic acid is mainly given in a hospital setting in Norway—and thus, not registered in NorPD. Hence, the frequency of use of this drug is probably underestimated. It should be recalled that the number of patients treated with zoledronic acid was much lower at the time the study was conducted than today. Moreover, zoledronic acid is also used in the treatment of hypercalcemia and as adjuvant treatment to patients with cancer. Other treatment options than BPs were used by very few.

Regarding adherence among the included participants with at least two filled prescriptions, 80% of the women and 77% of the men had an MPR ≥ 80 during follow-up. However, 15% of the women and 26% of the men filled only one prescription or used oral BP for less than 6 months (counted as unexposed in this study), illustrating that there is still a problem regarding adherence. Nevertheless, the adherence observed in our study is satisfactory compared to most studies on oral BPs [26].

We may assume that the compliance and adherence is better in RCTs with close follow-up than in the general patient population. Moreover, the route of administration will impact the adherence. Most probably the adherence is better for zoledronic acid that is given intravenously once yearly, and for denosumab administered subcutaneously every 6 months than for alendronate orally once a week.

The strength of our study is the population-based design and the large registers. The HUNT study also includes substantial information regarding risk for fractures, which reduces residual confounding. An additional strength is that we had access to calculated FRAX scores based on data collected in the HUNT study and the NorPD in all participants, although without BMD. Further, we have used time-dependent drug exposure to evaluate the effect of BPs during the follow-up period, thus avoiding immortal time bias, which refers to a period of follow-up during which, by design, death, or the study outcome cannot occur [36].

An important limitation of any observational study is that there may be unmeasured confounders which may lead to biased estimates [7]. In our study, there are obvious differences between individuals exposed and not exposed to BPs, and such potential confounders are not necessarily captured by age and FRAX. However, the fact that estimates from a sensitivity analysis using a stratified propensity score model based on a large number of potential confounders led to broadly similar results is reassuring. Furthermore, our results are comparable to a review including RCTs addressing the effect of BPs on non-vertebral fractures. In a total of five RCTs including 30.000 postmenopausal women, the HR for any non-vertebral fracture was 0.76 (0.70–0.81) among users of BPs [37].

Not all types of fractures could be included in the study as the fracture register only includes fractures of the hip and forearm. Further, calcium supplements with or without vitamin D were not included in the analyses as they are available without prescription in Norway. Among men, there were both few users of BPs as well as few fractures, with a resulting low statistical power. As noted above, the use of zoledronic acid may be underestimated since drugs given in hospital are not included in the NorPD. However, any uncaptured zoledronic acid use in the nominally untreated population would only conservatively bias the treatment effect we report here. As in all pharmacoepidemiological studies based on filled prescriptions, there is the possibility that individuals may accumulate old drugs. In general, in such studies there is no way to ensure that individuals actually use medication as prescribed. At last we only have information regarding prescription of BPs from 2004 (when NorPD was established) and cannot distinguish between treatment naïve and historical users.

The conclusion of this population-based study is that women treated with BPs displayed a reduced risk of fractures at the hip and forearm comparable to results from RCTs. The relative risk reduction was higher among women with secondary osteoporosis compared to primary osteoporosis. The most pronounced relative risk reduction was seen in women with glucocorticoid-induced osteoporosis.

(DOCX 15 kb)