A simple score for estimating the long-term risk of fracture in patients using oral glucocorticoids
From the 1Utrecht Institute for Pharmaceutical Sciences, Utrecht University, Utrecht, the Netherlands, 2MRC Epidemiology Resource Centre, University of Southampton, Southampton, UK, 3Procter & Gamble Pharmaceuticals, Egham, UK, 4Department of Rheumatology, University Hospital, Maastricht, Netherlands, 5Limburg University Center, Diepenbeek, Belgium, and 6Erasmus University Medical Center, Rotterdam, Netherlands
Received 29 September 2004 and in revised form 12 January 2005
| Summary |
|---|
|
|
|---|
Background: Previous analyses of risk factors for glucocorticoid (GC)-induced osteoporosis have focused on the estimation of relative rather than absolute fracture probability.
Aim: To estimate risk scores for the individual probability of fracture in GC users.
Design: Retrospective data analysis.
Methods: We evaluated all patients aged 40 years or older with a prescription for oral GCs in the General Practice Research Database (GPRD), which comprises the computerized medical records of around 7 million UK subjects. Individual risk factors for osteoporotic fractures were identified, and combined in a predictive model for 10-year absolute fracture risk.
Results: Of 191 752 oral GC users aged
40 years, 7412 experienced an osteoporotic fracture. Several characteristics independently contributed to the fracture risk score (GC therapy, age, gender, fall history, fracture history, body mass index, smoking, previous diagnoses, use of medication, recent hospitalization and indication for GC treatment). Scores of 30, 40 and 50 corresponded to absolute 5-year fracture risks of 6.2%, 15.3% and 35.2%, respectively. A woman aged 65 years with RA, low BMI, and a previous history of fracture and falls, who used 15 mg GC daily (total risk score 54) would have a 5-year fracture risk of 47% (a man with similar history, 30.1%). Short-term use of high-dose GC therapy (
30 mg) was associated with only a small increased risk of osteoporotic fracture (RR 1.21, 95%CI 1.041.42) in patients with a history of GC use.
Discussion: This risk score helps to predict an individual's risk of fracture during GC use. Decisions about bone protection treatment could be based on long-term risks of fracture.
| Introduction |
|---|
|
|
|---|
Oral glucocorticoids (GCs), also known as oral corticosteroids, are potent anti-inflammatory agents, and are widely used for the treatment of a variety of inflammatory and allergic disorders, such as rheumatoid arthritis or asthma. GC-induced osteoporosis was identified more than 60 years ago, when Cushing first described the tendency of patients with excess endogenous GC to develop fractures.1 Therapy with oral GC can lead to rapid loss of bone mineral density (BMD) and to an increased risk of fracture. A large epidemiological study reported dose-related increased rates of fracture in GC users,2,3 with similar findings in several smaller studies.48 Guidelines for its prevention and therapy of GC-induced osteoporosis have now been developed.9,10
The epidemiological information on the risks of fracture in oral GC users now concern relative risks (RR) of fracture, comparing the rate of fracture in GC users to that of control patients. There is limited information available to calculate the absolute fracture risk of an individual, specific for their age, sex and GC dose. Further, the relative contributions of GC use, bone-related risk factors such as fracture history, and extraskeletal risk factors such as the risk for falls, remain uncertain. In cardiovascular medicine, risk scores have been successfully developed to predict the long-term outcome of coronary heart disease in individuals. Thus, information on age, sex, serum cholesterol, blood pressure and other clinical characteristics can be used to estimate the absolute 5-year risk of death from cardiovascular disease.11 Active intervention is then only considered for patients whose long-term risk is above a certain threshold. It has been proposed that treatment decisions should be based primarily on absolute long-term risks of fracture, rather than on relative risks.12,13 The primary objective of this study was to estimate the long-term risks of fracture in patients using oral GCs by developing a predictive model.
| Methods |
|---|
|
|
|---|
Information for the study was obtained from the General Practice Research Database (GPRD), which comprises the computerized medical records of general practitioners. General practitioners (GP) play a key role in the UK health care system, as they are responsible for primary health care and specialist referrals. Patients are semi-permanently affiliated to a practice, which centralizes the medical information from the GPs, specialist referrals and hospitalizations. The data recorded in the GPRD include demographic information, prescription details, clinical events, preventive care provided, specialist referrals, hospital admissions and their major outcomes.14
Study population
Details of the methods used in the investigation and of the overall fracture rates (relative to a control group of non-users) are available elsewhere.2,3 The study population consisted of all patients aged 40 years or older with a prescription for oral GCs during the period of GPRD data collection. They were followed from the first GC prescription at age 40 years or older to the end of GPRD data collection (GPRD data collection started in 1987 and ended, for this study, in December 1997). The total follow-up period for each study patient was classified into periods of current and past GC exposure. As a patient's use of oral GCs could change over time, patients could move between these exposure periods over time. Each current exposure period started at the date of the prescribing of an oral GC; the duration of this current exposure period was based on the expected duration of GC use plus 3 months. In case of a repeat GC prescription within a current exposure period, this period was extended with the duration of use of the repeat prescription plus 3 months. In case of overlap between two prescriptions (i.e., a repeat prescription given within the duration of use for a previous prescription), the overlap days were added to that of the repeat prescription. The maximum period of current exposure following a GC prescription was set at 6 months. Daily and prior cumulative GC dose were assessed at the start of each current exposure period. The daily dose of each GC prescription was based on the prescribed daily dose by the GP, as obtained from the written dosage instructions in the medical records. Using this daily dose and the information on the total amount of GC prescribed to the patient, the expected duration of GC use was estimated. In case of missing data, the median expected duration of use (based on data from patients of similar age and sex) was used. The period of past exposure was the time period from GC discontinuation (end of current exposure period) up to a new oral GC prescription or up to the end of data collection (whichever end-point came first). The past-exposure period was divided into periods of 6 months.
Cases were patients who had a clinical osteoporotic fracture during follow-up. These were defined as fractures that occurred at one of the following sites: radius/ulna, humerus, rib, femur/hip, pelvis or vertebrae. An earlier validation study found that the vertebral fractures as recorded in GPRD mostly concerned vertebral fractures that were clinically symptomatic and were confirmed radiographically.15 Systematic morphometry of vertebral fractures was not routinely done by the GP. The history of any type of fracture prior to the first GC prescription was determined. Also, the occurrence of osteoporotic fracture at any other site during follow-up was noted. In order to exclude fractures that occurred at the same time but were recorded at a different time in the medical record, any fractures that occurred in the prior 3 months period were not included.
The other risk factors considered in the study included a recorded history of falls in the previous 6 months. The body mass index (BMI) and smoking history recorded around the time of the first GC prescription were also assessed, where available. The analysis also evaluated the number of risk factors, based on the presence of diseases and drugs that have been associated with an increased risk of fracture in a previous GPRD study.16 These included a history of anaemia, dementia, and cerebrovascular disease, or prescriptions in the previous 6 months for anticonvulsants, antiarrythmics, hypnotics/anxiolytics, antidepressants, or anti-Parkinsonian drugs. Information on risk factors was collected at the date of each GC prescription or at the date of starting each 6-monthly period of past exposure.
The indication for GC treatment was based on the illnesses recorded in the medical records in the period of time from 6 months before the first prescription up to the last GC prescription. The following diseases were noted: respiratory disease (ICD9 490496), rheumatoid arthritis (RA) (714), polymyalgia rheumatica (725), other connective tissue disorders (710713, 715719), non-infectious enteritis and colitis (555558), dermatitis (692), other inflammatory skin disorders (690, 691, 693698), urticaria (708), facial nerve disorders (351), other peripheral nervous system disorders (350, 352359) and polyarteritis (446). Hospitalization for any of the diseases indicated for GC treatment in the 12 months before was also measured. Risk factors were included as categorical variables (present vs. absent). With respect to the fracture risks in different GC indications, patients were classified according to GC indication, and each indication was then compared to patients with respiratory disease.
Statistical methods
Cox proportional hazards models were used to estimate the long-term risks of fracture. For each set of patient characteristics, the Cox model allows calculation of an individual's probability of fracture (i.e. survivor function). We first fitted regression models with daily and cumulative GC exposure (square-root-transformed) and an interaction term between daily and cumulative GC exposure. Past GC exposure was considered the reference category. We additionally investigated models with quadratic and cubic GC exposure variables, to address any non-linear relationship between GC exposure and fracture risk. Then regression models were fitted with the exposure variables, risk factors and the various GC indications. Backward regression used a significance level of 0.05. For age, sex, the risk factors and GC indications, we also investigated possible statistical interactions with GC exposure. Only the interaction with age was strong enough to merit inclusion into the final Cox model. The beta-coefficients of this Cox model (the exponentials of which constitute the RRs) were converted into integer risk scores. The value of each integer is the rounded sum of the predictors of the Cox model multiplied by 10. Because of the time-dependent exposure variables, the risk score of a patient was averaged over the total follow-up period. The 5-year risk of fracture was then estimated using these scores. This score represents the probability of fracturing, conditional on patient survival. The 10-year risk of fracture was based on extrapolating the 5-year risk; the average hazard rate was calculated (log of risk) and then applied over 10 years. Various methods were used to test the fitting of the Cox models.17 This included the testing of the proportional hazards assumption. We also compared the observed 5-year probability of fracture (based on the Kaplan-Meier estimate) to the probability predicted by the Cox model. This was done by dividing the study population into 10 groups, based on the predicted probability of fracture. The observed and predicted probabilities for fracture were then compared.
An analysis was conducted estimating the long-term risks of fracture for populations with different underlying fracture rates. As vertebral fractures are under-diagnosed in GPRD,15 we used the morphometric vertebral fracture rates from the European Prospective Osteoporosis Study (EPOS)18 and compared them to the GPRD rates. The log of the ratio of these rates was added to the risk score. In effect, this analysis replaced the underlying baseline fracture rate of the GPRD with that from another population. The underlying assumption was that the effects of GC exposure and of risk factors are similar between the different populations (i.e. RRs of GC exposure are identical).
Validation of predictive model
We also validated the risk score in another dataset, The Health Improvement Network (THIN) Research Database, which contains computerized medical records of patients at general practices in the UK, similar to the GPRD. The study population consisted of all women aged 50 years or older who were registered at one of the THIN practices, and who were prescribed an oral GC after 1997. It included 33 330 women, with 1489 clinical osteoporotic fractures.
| Results |
|---|
|
|
|---|
The study population consisted of 191 752 patients who were aged
40 years and prescribed an oral GC. The patients received a mean of 7.7 GC prescriptions (median 2). The total follow-up period (mean 2.5 years per person) was divided into periods of current and past GC exposure. About 59.5% of the total follow-up was classified as past exposure. There were 7412 patients with a clinical osteoporotic fracture (2144 hip and 1269 clinical vertebral fractures).
Table 1 shows the RRs of fracture for age, sex, risk factors and GC indication. Strong risk factors for fracture included age, sex, low BMI, fall and fracture history. Patients with RA had an increased risk of fractures, compared to those with respiratory disease. GC users with a low BMI and fall history had a five-fold increased risk of hip fracture (RR 5.37, 95%CI 2.3812.09). For GC users with both fracture and fall history, the hip fracture risk was similarly increased (RR 5.40, 95%CI 3.618.08). The RRs for fracture stratified by daily and cumulative GC dose are shown in Table 2. Fracture risk was increased for all doses, with the exception of short-term use of GC at a daily dose <2.5 mg/day. In patients with higher daily dose (
30 mg/day) and low cumulative dose (
1 g), GCs were often used for only short periods of time (average GC duration 10 days). This intermittent high-dose GC therapy was associated with only a small increased risk of clinical osteoporotic fracture. The age- and sex-adjusted RR for fracture was 1.20 (95%CI 0.981.46) in first-time GC users and 1.21 (95%CI 1.041.42) in patients with a high dose of GC at least 3 months after the end of prior GC use. These figures were 2.38 (95%CI 1.523.73) and 1.50 (95%CI 0.972.31) for clinical vertebral fractures and 0.86 (95%CI 0.571.30) and 0.78 (95%CI 0.541.14) for hip fractures, respectively.
|
|
Table 3 shows the risk score for various patient characteristics, with an example calculation shown in Table 4. The 5-year risks for a clinical osteoporotic fracture for patients with total scores of 30, 40, and 50 were 6.2% (95%CI 6.06.4), 15.3% (95%CI 14.915.7), and 35.2% (95%CI 33.736.6%), respectively (Figure 1). The median risk score at start of GC therapy (excluding GC use) was 19 for women aged 4049 years (75% percentile: 21) and 13 (15) for men of this age. These figures were 28 (31) for women and 22 (24) for men aged 6069 years, and 37 (41) for women and 31 (32) for men aged 8089 years.
|
|
|
Receiver operator characteristic (ROC) curves, which plot sensitivity (true positive rate) vs. (1 specificity) (false positive rate) and the areas under the ROC curve were estimated. For clinical osteoporotic fractures, the area under the ROC curve was 69.7%; for hip fracture, 78.1%; and for clinical vertebral fracture, 74.8%.
When applying the risk score as developed in GPRD to another dataset (THIN), the risk score equally differentiated between high- and low-risk patients (in the original GPRD population, the RR for an increase of 10 in the risk score was 2.63 among elderly women; in THIN, it was 2.74). These figures were 2.60 and 2.48 for clinical vertebral fractures, and 2.64 and 2.58 for hip fractures, respectively.
| Discussion |
|---|
|
|
|---|
We have developed a risk score which provides an easily applicable method of estimating a patient's individual risk of fracture, based on routinely available clinical information. Our data suggest that the long-term risks of fracture can be substantial in GC users; for a woman aged 65 years using 7.5 mg prednisolone daily, there is a one in four likelihood of suffering a clinical osteoporotic fracture over 10 years.
The risk score as developed in this study allows calculation of the long-term risk of fracture in patients using oral GCs. As in cardiovascular medicine, such information on individual risks can be useful in identifying patients who require active intervention and monitoring. Current approaches of intervention thresholds based on (e.g.) age, may result in the over-treatment of elderly patients with below-average risks and under-treatment of younger patients with above-average risks. Further study is required to determine the threshold of individual risk at which intervention will become cost-effective in patients using oral GCs.
In postmenopausal osteoporosis, intervention thresholds are currently determined largely by assessment of BMD or a prior history of fracture. However, fracture probability varies substantially with age at any given level of BMD. Older people have much higher fracture risks than younger people, even when their BMDs are similar.12 For GC users, there are even stronger grounds for using BMD as one of several risk indicators, rather than as the single determinant for intervention. There is growing evidence that the increased risk of fracture in GC users is not only related to a reduction in amount of bone (BMD), but also to the quality of bone. Preliminary results of an animal study indicated that bone quality (osteocyte viability) was an important determinant of bone strength, independent of BMD.19 A meta-analysis showed that fractures occur at much higher rates in GC users than expected on the basis of BMD differences,10 and in a comparison of GC users and non-users, GC users had considerably higher fracture risks than non-users at similar levels of BMD.20 Several treatment guidelines now recommend that BMD is only measured in patients with an intermediate, and not high, risk of fracture. Further studies are warranted to examine the value of measuring BMD in this context.
In addition to the effects on bone quality, oral GC therapy may affect the risk of fracture through extraskeletal mechanisms. We found that patients on GC had an increased risk of falling,4 and fall history contributed to our risk score. Although previous studies examining the relationship between GC use and fracture risk have not directly examined the frequency of falling, a British case-control study did suggest that GC users had increased levels of frailty, physical inactivity, and immobility.21 And oral GC therapy may lead to muscle weakness. Thus, a history of falls (a well-documented risk factor for subsequent falls22), is also an important risk factor for fractures in oral GC users.
The daily dose in this study was based on the dosage prescribed by the general practitioner, and could therefore vary over time. In a previous study that was used to calculate RRs of fracture in the study population, the daily dose was derived from a single average (calculated from the total number of tablets prescribed and length of follow-up).2,3 The use of a time-varying measure of exposure may be preferable. Recent research found that fracture risk increases rapidly in GC users and that daily, but not cumulative, GC dose is a strong predictor for incident fracture in GC users.20 Nevertheless, two other studies have reported that cumulative GC dose was a stronger predictor of fracture than daily dose.6,7 The main limitation of these two studies was their reliance on prevalent fracture: patients with high cumulative exposure will have more prevalent fractures even if the incidence of new cases is not increasing with larger cumulative dose, as their time on treatment and disease duration will be longer. Prevalence is a function of underlying incidence rate of fracture as well as of duration of therapy or disease. We believe that it is inappropriate to make inferences on the relative contributions of daily and cumulative dose on the basis of prevalent cases only.
The risk estimates in this study are based on historical data. It has been argued that such data can only be an estimate for prospective prediction of fracture risk, because populations and circumstances are continuously changing.23 At the time of this study, the awareness of fracture as a side-effect of GC therapy was probably limited, as reflected by the low prevalence of use of bone-active medication.24 This may have improved over time, with better information now available on skeletally protective diet or exercise. Also, the risk estimates in this study may not be generalized to other populations, as GC effects may vary. It would be more appropriate to view the risk estimates in this study as a tool to improve the prediction of GC-related fractures, rather than as definitive estimates applicable to every patient.
Our findings are based on a complex mathematical model. We evaluated the key underlying assumptions used, and its overall predictive capacity performed well. However, we did not evaluate all possible interactions between the risk factors and for certain risk factor combinations; the model may therefore have over- or under-estimated fracture risks. Another limitation was that we did not have information on all risk factors for fracture (such as BMD, exercise or diet), which would improve the accuracy of prediction for an individual patient. We also did not include in the model use of bone-protective medications, such as bisphosphonates and hormone replacement. Use of such medication prior to a fracture was infrequent in the study period,24 and the estimates for the risk score did not vary substantially after excluding patients who used these drugs.
In conclusion, patients using GC may have a substantial absolute long-term risk of fracture. A simple risk score based on GC dose and indication, patient BMI, fracture and fall history, and history of other disease and concomitant drug use, can help to quantify this absolute individual risk. This score can also be used to target preventative or investigative action to patients with higher long-term risks.
| Acknowledgments |
|---|
The study was supported by a grant from Procter & Gamble Pharmaceuticals Limited.
| Footnotes |
|---|
Address correspondence to Professor C. Cooper, MRC Epidemiology Resource Centre, Southampton General Hospital, Southampton SO16 6YD. e-mail: cc{at}mrc.soton.ac.uk
| References |
|---|
|
|
|---|
1. Cushing H. The basophilic adenomas of the pituitary body and their clinical manifestations (pituitary basophilism). Bull Johns Hopkins Hosp 1932; 1:13792.
2. van Staa TP, Leufkens HGM, Abenhaim L, Zhang B, Cooper C. Use of oral corticosteroids and risk of fractures. J Bone Miner Res 2000; 15:9931000.[CrossRef][Web of Science][Medline]
3. van Staa TP, Leufkens HGM, Abenhaim L, Zhang B, Cooper C. Fractures and oral corticosteroids: relationship to daily and cumulative dose. Rheumatol 2000; 39:13839.
4. Naganathan V, Jones G, Nash P, Nicholson G, Eisman J, Sambrook PN. Vertebral fracture risk with long-term corticosteroid therapy. Arch Intern Med 2000; 160:291722.
5. Cooper C, Coupland C, Mitchell M. Rheumatoid arthritis, corticosteroid therapy and hip fracture. Ann Rheum Dis 1995; 54:4952.
6. Walsh LJ, Wong CA, Oborne J, Cooper S, Lewis SA, Pringle M, Hubbard R, Tattersfield AE. Adverse effects of oral corticosteroids in relation to dose in patients with lung disease. Thorax 2001; 56:27984.
7. Dykman TR, Gluck OS, Murphy WA, Hahn TJ, Hahn BH. Evaluation of factors associated with glucocorticoid-induced osteopenia in patients with rheumatic diseases. Arthritis Rheum 1985; 28:3618.[Web of Science][Medline]
8. van Staa TP, Leufkens HGM, Cooper C. The epidemiology of corticosteroid-induced osteoporosis: a meta-analysis. Osteoporos Int 2002; 13:77787.[CrossRef][Web of Science][Medline]
9. American College of Rheumatology Ad Hoc Committee on Glucocorticoid-Induced Osteoporosis. Recommendations for the prevention and treatment of glucocorticoid-induced osteoporosis. Arthritis Rheum 2001; 44:1496503.[CrossRef][Web of Science][Medline]
10. Glucocorticoid-induced osteoporosis: guidelines for prevention and treatment. London, Royal College of Physicians, 2002.
11. Pocock SJ, McCormack V, Gueyffier F, Boutitie F, Fagard RH, Boissel J-P. A score for predicting risk of death from cardiovascular disease in adults with raised blood pressure, based on individual data from randomized controlled trials. Br Med J 2001; 323:7581.
12. Kanis JA, Johnell O, Oden A, Dawson A, de Laet C, Jonsson B. Ten year probabilities of osteoporotic fractures according to BMD and diagnostic thresholds. Osteoporos Int 2001; 12:98995.[CrossRef][Web of Science][Medline]
13. Kanis JA, Johnell O, Oden A, de Laet C, Oglesby A, Jonsson B. Intervention thresholds for osteoporosis. Bone 2002; 31:2631.[Medline]
14. Walley T, Mantgani A. The UK General Practice Research Database. Lancet 1997; 350:10979.[CrossRef][Web of Science][Medline]
15. van Staa TP, Abenhaim L, Cooper C, Zhang B, Leufkens HGM. The use of a large pharmacoepidemiological database to study exposure to oral corticosteroids and risk of fractures: validation of study population and results. Pharmacoepi Drug Saf 2000; 9:35966.[CrossRef]
16. van Staa TP, Leufkens HGM, Cooper C. Utility of medical and drug history in fracture risk prediction among men and women. Bone 2002; 31:50814.[Medline]
17. Harrell HE, Lee KL, Mark DB. Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med 1996; 15:36187.[CrossRef][Web of Science][Medline]
18. Felsenberg D, Silman AJ, Lunt M, et al. for the European Prospective Osteoporosis Study (EPOS) group. Incidence of vertebral fracture in Europe: results from the European prospective Osteoporosis Study (EPOS). J Bone Miner Res 2002; 17:71624.[CrossRef][Web of Science][Medline]
19. Weinstein RS, Powers CC, Parfitt AM, Manolagas SC. Preservation of osteocyte viability by bisphosphonates contributes to bone strength in glucocorticoid-treated mice independently of BMD: an unappreciated determinant of bone strength. J Bone Miner Res 2002; 17(Suppl. 1):S156.
20. van Staa TP, Laan RF, Barton I, Cohen S, Reid DM, Cooper C. Predictors and bone density threshold for vertebral fracture in patients using oral glucocorticoids. Arthritis Rheum 2003; 48:32249.[CrossRef][Web of Science][Medline]
21. Cooper C, Coupland C, Mitchell M. Rheumatoid arthritis, corticosteroid therapy and hip fracture. Ann Rheum Dis 1995; 54:4952.
22. Schwartz AV, Capezuti E, Grisso JA. Falls as risk factors for fractures. In: Marcus R, Feldman D, Kelsey J, eds. Osteoporosis, 2nd edn. San Diego, Academic Press, 2001:58598.
23. Tunstall-Pedoe H. Absolute is inappropriate for quantitative risk estimation. Br Med J 2000; 320:723.
24. van Staa TP, Cooper C, Abenhaim L, Begaud B, Zhang B, Leufkens HGM. Utilisation of oral corticosteroids in the United Kingdom. Q J Med 2000; 93:10511.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
S. Parikh, H. Mogun, J. Avorn, and D. H. Solomon Osteoporosis Medication Use in Nursing Home Patients With Fractures in 1 US State Arch Intern Med, May 26, 2008; 168(10): 1111 - 1115. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. P. van Staa, P. Geusens, B. Zhang, H. G. M. Leufkens, A. Boonen, and C. Cooper Individual fracture risk and the cost-effectiveness of bisphosphonates in patients using oral glucocorticoids Rheumatology, March 1, 2007; 46(3): 460 - 466. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. G. Raisz Screening for Osteoporosis N. Engl. J. Med., July 14, 2005; 353(2): 164 - 171. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



