Article Text
Abstract
OBJECTIVES To examine relations between osteoporosis and low bone mass and demographic and clinical variables in patients with rheumatoid arthritis (RA), in an attempt to develop a data driven clinical tool for identification of patients at high risk of osteoporosis.
METHODS All patients were recruited from a county based register and were examined cross sectionally with a variety of clinical and health status measures as well as bone density measures (anteroposterior spine L2-4, total hip, and femoral neck). Associations between osteoporosis (T score ⩽−2.5SD) and low bone mass (T score ⩽−1SD), on the one hand, and demographic and clinical measures, on the other, were examined bivariately and by logistic regression analyses.
RESULTS 394 patients with a mean age of 54.8 years were examined. The percentages having osteoporosis/low bone mass were 16.8/45.8, 14.7/54.5 and 14.7/55.5 in spine L2-4, total hip, and femoral neck, respectively. Osteoporosis and low bone mass were bivariately related to age, body mass index (BMI), disease duration, disease process measures, presence of deformed joints, physical disability, current use of corticosteroids, and history of non-vertebral fracture. In multivariate analyses, age >60 years, low BMI, and current use of corticosteroids were consistently related to osteoporosis and to low bone mass at all sites. The presence of deformed joints was associated with osteoporosis at the total hip, and a history of previous non-vertebral fracture with osteoporosis at the femoral neck. The Modified Health Assessment Questionnaire (MHAQ) ⩾1.5 and non-vertebral fracture were also independently associated with low bone mass at the hip. The logistic regression analyses models could, however, only predict osteoporosis with a sensitivity of about 50–60% and a specificity of 80–90% at the various measurement sites, and low bone mass with a sensitivity and specificity of about 70%.
CONCLUSION Consideration of demographic and disease markers may be of some help in predicting presence of osteoporosis or low bone mass, but a combination of markers cannot be used as a clinical tool with sufficient sensitivity and specificity for the identification of osteoporosis or low bone mass in patients with RA.
- rheumatoid arthritis
- osteoporosis
- risk factors
- epidemiology
Statistics from Altmetric.com
Osteoporosis is recognised as a complication of rheumatoid arthritis (RA) (see review1). Several studies have explored the relation between osteoporosis, on the one hand, and demographic and disease related variables, on the other.2 ,3 Use of corticosteroids has been found to be related to osteoporosis in RA,4 and several studies have also indicated that the development of bone loss (osteoporosis) is related to disease activity, especially early in the disease course.5 ,6
Guidelines have been issued on the management of osteoporosis with corticosteroids.7 ,8 These guidelines have also included directions on which patients to select for diagnostic osteoporosis procedures, including bone mineral densitometry (BMD). Such guidelines have partly been based on published scientific data, partly on expert opinions, and partly on consenting group processes.
Longitudinal data are generally required to make statements about causal relations. However, for the clinician facing individual patients retrospective evaluation of a series of markers of disease activity and severity is usually not available. Thus the clinician must often rely on cross sectional data to identify patients at high risk of osteoporosis being candidates for diagnostic procedures. Recently, the need for a clinical algorithm for osteoporosis in patients with RA has been discussed.9 A preliminary proposal was made to measure bone density in patients with RA fulfilling two out of three of the following criteria: (a) age above 60, (b) immobility, and (c) high disease activity, defined as mean C reactive protein (CRP) above 20 mg/l or persistently raised erythrocyte sedimentation rate (ESR) above 20 mm/1st h.9
Previous studies on the identification of clinical and demographic variables predicting high risk of osteoporosis have been performed in rather limited patient numbers recruited from clinical settings, with a possible overrepresentation of patients with RA using corticosteroids and with active and severe disease. Our present approach was to examine a large number of patients recruited from a county based RA register in an attempt to develop a clinical tool to identify female patients who require diagnostic procedures and treatment.
Materials and methods
SETTING AND INCLUSION CRITERIA
A register of patients with RA has been established in the county of Oslo since 1994.10 Inclusion criteria are a diagnosis of RA11 and a residential address in Oslo. The register is continuously updated with new cases and also withdrawals due to death or new addresses outside Oslo.
The register has been shown to be about 85% complete in the age range 18–79 years, and it has been successful in providing data on incidence,12 prevalence, and disease severity of RA.10
Inclusion criteria for the present study were enrolment in the Oslo RA register, female sex, and age between 18 and 70 years (born 1926 or later). The number of patients with RA fulfilling the inclusion criteria was 721. Of these 721 patients, 394 white subjects were willing to meet for clinical examination and BMD measurement (attendance rate 55%).13 All patients were younger than 70 years when completing the data collection (clinical data and BMD).
COLLECTION OF CLINICAL DATA
The collection of data was performed during a 15 month period between September 1996 and December 1997. Clinical data were collected partly by questionnaire, partly by interview and examination by trained research nurses under the supervision of a rheumatologist.
The following demographic variables were assessed: age, disease duration, body weight (light indoor clothing), and height. From the two last measures the body mass index (BMI) was computed. Comorbidities, exposure to oestradiol (current or not current user), menopause status, and previous oophorectomy were recorded as well as a history of previous non-vertebral fracture. Lifestyle variables included smoking (current, previous, never).
Disease activity was assessed by the following eight variables: 28 tender and 28 swollen joint counts, patients' and investigators' global assessment of disease activity (on a five point categorical scale and 100 mm visual analogue scale (VAS), respectively), pain severity and fatigue on 100 mm VAS, and acute phase reactants by ESR and CRP. The disease activity score (DAS) was computed using 28 joint counts.14
Damage was assessed by the number of deformed joints (score 0–18), and functional health status by the Modified Health Assessment Questionnaire (MHAQ) (range 1–4).15
Use of corticosteroids was assessed as a categorical variable (never user, previous user, current user).
BMD MEASUREMENTS
The BMD measurements of the hip (total hip and femoral neck) and the anteroposterior lumbar spine (L2–4) were performed by three trained technicians using the same dual energyx ray absorptiometry equipment (Lunar Expert, Madison, Wisconsin). The machine was calibrated daily with a spine phantom supplied by Lunar. Measurements of the left hip and spine were done according to standardised procedures.13 In 14 patients no hip measurement was performed as they had had two hip replacements, and in 14 patients the right hip only was measured owing to replacement of severe destruction of the left hip. No spine measurement was performed in one patient.
To ensure uniformity of the results all the scans were reanalysed by one of the experienced technicians after all the data had been collected. The spine phantom precision error calculated as a coefficient of variation (CV%) was 0.9% for the whole measurement period (15 months). The in vivo reproducibility of BMD measurements was assessed from duplicate measurements in 31 healthy female hospital workers (mean age 56.1 years, range 50–66). The CV was 2.2% at the lumbar spine, 1.5% at the total hip, and 1.5% at the femoral neck. The interobserver variation for the three trained technicians performing measurements and analyses was assessed by pairwise comparisons of analyses of 30 randomly selected patients with RA. The precision varied from 0.7 to 1.4% at the spine L2–4, from 0.4 to 0.5% in the total hip, and from 0.5 to 0.8% in the femoral neck.
The BMD measurements in the RA population were compared with a pooled European/US reference population aged 20–69, based on different normal data from Europe and US, containing 9505 women with hip and 10 281 with spine measurements.16-24 The variation of mean BMD values among the geographical sites in this pooled population was about 1.3% (SD/mean).16 The reference database comprised measures of ambulatory subjects generally free from earlier fractures, chronic diseases, and drugs influencing bone metabolism (for example, corticosteroids, anticonvulsant drugs, and thyroxine).16
The BMD data for the pooled European/US reference population were provided by Lunar Corporation, Madison, including the data for T score estimations. For T score estimations the following mean (SD) BMD results for young adults (age 20–39) were used: femoral neck 0.98 (0.12) g/cm2, total hip 1.00 (0.12) g/cm2, and spine L2–4 1.20 (0.12) g/cm2 .
ETHICS AND LEGAL ASPECTS
The local ethic committee approved this study. The Data Inspectorate had approved the register of patients with RA in Oslo.
DATA ANALYSES
Data analysis was performed in the total female patient group (n=394) and, separately, for patients aged over 50 (n=271). All analyses were performed with SPSS, version 8.0.
Osteoporosis was defined as a T score ⩽−2.5SD,25 low bone mass as T score ⩽−1. After computing the number of patients with osteoporosis and low bone mass, bivariate comparisons of demographic, disease activity, and disease severity measures between patients with and without osteoporosis/low bone mass were performed, using two sided t tests (continuous variables) and χ2 tests (counts).
Independent variables to be used in the multivariate analyses were dichotomised as follows: age over and under 60 years, BMI and DAS divided by the two highest and lowest quartiles, presence and absence of deformed joints, MHAQ score ⩾1.5 (indicating physical disability), corticosteroids into current versus not users, and presence versus absence of previous non-vertebral fracture.
The occurrence of such possible predictors of osteoporosis (age ⩾60 years, low BMI, high DAS, deformed joints, MHAQ ⩾1.5, current user of corticosteroids, previous non-vertebral fracture) were compared between patients with and without osteoporosis/low bone mass. These possible predictors were subsequently entered into a logistic regression analysis applying osteoporosis/low bone mass at different sites as dependent variables. The logistic regression analyses were run with the enter procedure.
Results
PATIENTS
Table 1 shows the demographic and clinical characteristics of the 394 female patients with RA finally included in the study, born 1926 and later, and the 271 patients aged between 50 and 70.
Of all the women with RA from the register, born 1926 or later (n=721), there were no statistically significant differences in age, disease duration, and rheumatoid factor positivity between those measured and not measured. A group of 109 patients had clinical examinations, without bone mass measurement. The clinical and disease characteristics, including disease activity and severity measures, in these patients did not differ from those whose BMD was measured. Thus, the final sample of 394 patients with RA seems to be representative of the overall RA population in the county.
BMD MEASURES
The BMD values (mean (SD) in g/cm2) for the whole patient group were 1.10 (0.21) at the spine L2–4, 0.87 (0.16) at the total hip, and 0.85 (0.16) at the femoral neck. Corresponding values in the 50–70 age group were 1.05 (0.20), 0.83 (0.16), and 0.80 (0.15).
The percentages (95% confidence interval (CI)) having osteoporosis (T score ⩽−2.5SD) were in the total patient group 16.8 (13.1 to 20.5) at L2–4, 14.7 (11.1 to 18.3) at the total hip, and 14.7 (11.1 to 18.3) at the femoral neck. Among those aged 50–70 the percentages were as expected higher (23.3 (18.3 to 28.3), 21.1 (16.2 to 26.0), and 20.7 (15.9 to 25.5)).The percentages (95%CI) with low bone mass (T score ⩽−1SD) were 45.8 (40.9 to 50.7), 54.5 (49.5 to 59.5), and 55.5 (50.5 to 60.5) in the spine L2–4, total hip, and femoral neck, respectively.
COMPARISON OF PATIENTS WITH AND WITHOUT OSTEOPOROSIS/LOW BONE MASS
Table 2 gives bivariate comparisons of demographic and disease activity and severity variables between patients with and without osteoporosis. Osteoporosis was statistically associated with increasing age, longer disease duration, lower BMI, increasing disease activity, and use of corticostertoids at all sites of measurement. The computed DAS also discriminated between patients with and without osteoporosis. Osteoporosis was more consistently related to disease activity and severity measures at the hip location than at the spine L2–4 site of measurement (table 2).
Similar patterns of differences were seen for demographic and disease variables between patients with and without low bone mass (data not shown).
RISK FACTORS FOR OSTEOPOROSIS
Tables 3-5 present bivariate comparisons and logistic regression analyses of the associations between osteoporosis and categorised risk factors (age, BMI, joint deformity, disability, DAS, current use of corticosteroids, and previous non-vertebral fracture). Table 6 presents the multivarate relations between low bone mass and the categorised risk factors.
Patients with osteoporosis differed for all tested variables (table 4) from those without osteoporosis of the total hip, whereas group differences were not found for physical disability (MHAQ ⩾1.5) between patients with and without osteoporosis at the spine L2–4 (table 3) and femoral neck (table 5). The same trends were found for patients aged 50–70, but owing to the lower sample size not all differences reached statistical significance (data not shown).
Age over 60, low BMI, and current use of corticosteroids were the most consistent risk factors in multivariate analyses for osteoporosis and low bone mass (tables 3-6). The presence of deformed joints was associated with osteoporosis at the total hip, and a history of previous non-vertebral fracture with osteoporosis at the femoral neck (tables 4 and 5). MHAQ ⩾1.5 and non-vertebral fracture were independently associated with low bone mass at the hip.
The classification tables showed that models generally reached a sensitivity between 50 and 60% and specificity between 80 and 90% for osteoporosis (tables 3-5) and a sensitivity and specificity of about 70% for low bone mass (table 6).
In addition to the model examining the main effects, an expanded model which included all first order interactions was fitted. The overall test for all the interactions turned out to be non-significant. Consequently, the main effects model had an adequate fit.
Discussion
This study confirms that patients with RA with osteoporosis and low bone mass generally have a more active and severe disease than patients without osteoporosis. Thus markers of disease activity and severity as well as treatment with corticosteroids indicate an increased risk of osteoporosis. However, this study also highlights that these associations with clinical markers are weak, as increasing age and low BMI were the strongest and most consistent predictors of osteoporosis at all sites, both in the total patient group and in the analyses of patients aged 50–70. Age and BMI are also consistent predictors in postmenopausal osteoporosis,26 suggesting that these demographic variables may be more important at a group level than the disease factors.
However, clinical markers are also important, even if the logistic regression models used had rather low sensitivity in detecting patients with osteoporosis. The most important clinical factor was the current use of corticosteroids, and then the presence of deformed joints, indicating structural joint damage. In this cross sectional approach the DAS and the MHAQ score did not turn out to be independent predictors of osteoporosis, and the MHAQ score only reached borderline statistical significance at one measurement site. Analyses were also performed with low bone mass as the dependent variable because the threshold of the T score ⩽−1SD may be as important when advocating treatment, especially in patients using corticosteroids.8
The results can be used to compute algorithms of some help in clinical practice. The predicted probabilities for a combination of clinical variables can be calculated by adding the parameter estimates from any combinations. For example, the results in table 4 indicate that the odds of having total hip osteoporosis for women in the oldest age group and lowest BMI group with inactive and not severe disease is:
- logit
- (constant, age ⩾60, BMI two lowest quartiles, not deformed joints, MHAQ ⩽1.5, DAS two lowest quartiles, no use of corticosteroids, and no previous non-vertebral fracture) = logit (−5.63 + 2.58 + 1.82 + 0 + 0 + 0 + 0 + 0),
- while
- the odds for having osteoporosis for women with the same age and BMI characteristics, but with severe and active disease and previous non-vertebral fracture is:
- logit
- (constant, age ⩾60, BMI two lowest quartiles, deformed joints, MHAQ ⩾1.5, DAS two highest quartiles, current use of corticosteroids, and previous non-vertebral fracture) = logit (−5.63 + 2.58 + 1.82 + 0.75 + 0.42 + 0.01 + 1.00 + 0.67).
The log odds ratio (OR) for women over 60 and with low BMI with active and severe disease and with previous non-vertebral fracture compared with those without active disease and with no previous non-vertebral fracture is:
- (logit
- (−5.63 + 2.58 + 1.82 + 0.75 + 0.42 + 0.01 + 1.00 + 0.67)) − (logit (−5.63 + 2.58 + 1.82 + 0 + 0 + 0 + 0 + 0)) = 1.62 −(−1.23) = 2.85 .
Thus the OR for having total hip osteoporosis for those with compared with those without active and severe disease and a previous history of non-vertebral fracture is: e2.85 = 17.29, where e is the base of the natural logarithms, approximately 2.718.
It can be argued that other disease activity measures should have been included in the model. A lot of exploratory analyses were done before we decided to use the present model, which incorporates various dimensions of measures that are recommended for longitudinal observational studies27 and also may be used in the routine clinical care of patients with RA.
For corticosteroids we chose to use current instead of ever use of corticosteroids as the independent variable as corticosteroid induced bone loss may be reversed when corticosteroid treatment is stopped.28 In the exploratory analyses we found that acute phase reactants, number of swollen joints, and rheumatoid factor were not independent predictors. From the group comparisons it was also shown that only marginal differences were seen for these variables between the patients with and without osteoporosis (table 2). This result partially contradicts the results of a previous study showing that the mean ESR measured during the six months before the BMD measurement was negatively associated with hip BMD.29 We had expected a stronger effect of the MHAQ level on osteoporosis, as this previous Dutch study also found an association between physical disability and BMD.29 The level of MHAQ is, however, closely related to age, and age came out as the strongest overall predictor.
Previous studies on osteoporosis in RA have mainly focused on selected groups of patients recruited from referral centres.2 ,3 ,30 ,31 In our study we had a different approach, examining patients from a county based RA register. This register has been examined for its completeness and includes about 85% of the total underlying patient population,10 and prevalence data have shown a twofold increase in osteoporosis in patients with RA.13 Furthermore, the patients examined did not differ significantly from the patients on the register who were eligible for this study. Thus our results are relevant for all patients with RA as a group, but different results might have been found if similar analyses had been performed in selected patients with active and severe disease, or in a selected group of patients who were not receiving hormone replacement therapy. It should be noted that more than one quarter of the patients in this study were current users of oestradiol. However, the results were similar when we repeated the analyses in a subsample excluding users of oestradiol (data not shown).
The effect of inflammatory markers on osteoporosis might also have been different if longitudinal data had been analysed.6 ,29 For radiographic progression it is well known that the level of acute phase reactants correlates with structural damage.32 This correlation has been shown to be even stronger if, for example, a series of CRP measurements are related to radiographic damage.33
The important clinical outcome of low bone mass is symptomatic fracture. The incidence of hip and vertebral fractures is increased in patients with RA, but the relation between bone mass and subsequent fractures remains to be established.34 It is not certain whether data from patients with postmenopausal osteoporosis, showing a double risk of fracture for each SD reduction in BMD,35 ,36 is also applicable to patients with RA. In our study a history of non-vertebral fracture was used as a possible predictor, showing an independent relation with osteoporosis and low bone mass in the hip, but not in the lumbar spine.
Our aim was to develop a clinical tool which would be helpful in the identification of patients with RA who were likely to develop osteoporosis. We found that osteoporosis and low bone mass in RA are more consistenly associated with age, BMI, and current use of corticosteroids than disease variables. Osteoporosis could be predicted with an average sensitivity of only 50–60% and a specificity of 80–90% when the predictors were combined into logistic regression models. This indicates that consideration of risk factors may be helpful in clinical practice, but that many patients with osteoporosis will be missed if only those with known risk factors, including high disease activity and severity, are considered.
Acknowledgments
We thank Petter Mowinchel for statistical advice, our technicians, Ingerid Müller, Sidsel Arnkværn, and Espen Haavardsholm, for expert assistance with the BMD measurements, and research secretary Kirsten Mossin for keeping the Oslo Rheumatoid Arthritis Register updated.
References
Footnotes
This study was funded by grants from The Research Council of Norway, Lions Club International MD 104 Norway, the Norwegian Women Public Health Association, Trygve Gythfeldt and Wife's Legacy, Grethe Harbitz's Legacy, and Marie and Else Mustad's Legacy.