Abstract
Objective. To evaluate bone quality in patients with systemic lupus erythematosus (SLE) who were undergoing longterm glucocorticoid (GC) therapy, and to focus on the correlation between bone quality and organ damage.
Methods. Seventy-eight female patients with SLE and organ damage taking longterm GC, and 72 age-matched SLE patients without damage taking longterm GC were recruited for study. Clinical variables of interest included disease activity, cumulative organ damage (by Systemic Lupus International Collaborating Clinics/American College of Rheumatology Damage Index; SDI), major organ involvement (musculoskeletal damage and neuropsychiatric damage, etc.), and use of medication. Areal bone mineral density (aBMD) was measured by dual-energy X-ray absorptiometry. Bone geometry, volumetric BMD (vBMD), microarchitecture, and biomechanical properties were measured by high-resolution peripheral quantitative computed tomography (HR-pQCT).
Results. Patients were mean age of 45 years (SD 10) and 54% were postmenopausal. The median SDI score of the cohort was 1 (interquartile range 1−2, range 1−5). Compared with patients without damage, the prevalence of osteopenia at either total hip or lumbar spine was significantly higher, and there were trends of deterioration of bone geometry, vBMD, microarchitecture, and biomechanical properties in patients with organ damage. Potential risk factors for bone quality in patients with damage were screened by univariate analysis. During multiple regression analysis, SDI was the only clinical variable consistently associated with deterioration of vBMD and microarchitecture.
Conclusion. Cumulative organ damage consistently correlated with deterioration of vBMD and bone microarchitecture in SLE patients with damage on longterm GC therapy. HR-pQCT provides an insight into the underlying mechanism of bone loss in SLE.
- BONE MICROARCHITECTURE
- SYSTEMIC LUPUS ERYTHEMATOSUS
- ORGAN DAMAGE
- HIGH-RESOLUTION PERIPHERAL QUANTITATIVE COMPUTED TOMOGRAPHY
Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with a broad range of clinical and laboratory presentations involving almost all organs, including the neuropsychiatric, cardiovascular, and musculoskeletal systems. Women with SLE have higher risk of osteoporosis and fracture than the general population as a result of exposure to direct or indirect risk factors, such as treatment with glucocorticoid (GC), sun avoidance, or premature menopause caused by the disease itself or its treatment with cytotoxic drugs1. Prevalence of osteoporosis in Chinese patients with SLE ranged from 4% to 6% in premenopausal patients2 and 3% to 48% in postmenopausal patients3, depending on different measurement sites.
Many studies using dual-energy X-ray absorptiometry (DEXA) have shown that cumulative organ damage of SLE is a risk factor for decreased areal bone mineral density (aBMD)4,5. As a 2-dimensional projectional imaging technique, DEXA measures integral aBMD of cortical and trabecular bone and is confounded by bone geometry, and thus cannot truly reflect bone quality including microarchitecture and bone strength6. With advances in high-resolution imaging techniques, interest in investigating bone microarchitecture as an important factor in development of osteoporosis has increased. High-resolution peripheral quantitative computed tomography (HR-pQCT) capable of achieving an isotropic voxel size of 80 μm is now available for assessment of trabecular and cortical microarchitecture at the distal radius and tibia7. Compared with DEXA, HR-pQCT has the ability to distinguish cortical and trabecular bone with 3-dimensional images, higher resolution, and lower radiation exposure, as well as shorter scanning time. HR-pQCT has been validated against microcomputed tomography, which is usually considered the “gold standard” to assess bone microarchitecture8. Further, it is feasible to apply finite-element analysis based on the HR-pQCT images to predict biomechanical properties (bone strength)9. In our previous study, we found that the discriminatory power for vertebral fracture was greater for radial volumetric BMD (vBMD) and microarchitecture by HR-pQCT than that for central aBMD in patients with SLE10.
In this cross-sectional study, we aimed to (1) evaluate various bone quality measurements, including aBMD by DEXA and bone geometric, densitometric, microarchitectural, and biomechanical measures, by HR-pQCT in female Chinese patients with SLE and organ damage who were undergoing longterm GC therapy, compared with age-matched SLE patients with no organ damage; and (2) investigate the risk factors for deterioration of vBMD and microarchitectural and biomechanical variables in SLE patients with organ damage.
MATERIALS AND METHODS
Patients
Consecutive female patients with SLE who regularly attended the outpatient rheumatology clinic at Prince of Wales Hospital or Alice Ho Miu Ling Nethersole Hospital in Hong Kong were screened during September 2010 and July 2011. Patients were eligible if (1) they fulfilled the 1997 American College of Rheumatology (ACR) revised criteria for the classification of SLE11; (2) they were receiving prednisone at least 5 mg per day continuously for at least 1 year prior to study entry; and (3) they had at least 1 manifestation of organ damage defined as Systemic Lupus International Collaborating Clinics/ACR Damage Index (SDI) ≥ 112. Patients were ineligible if (1) they had a known metabolic disorder that could affect bone metabolism; (2) they were receiving treatment that could affect bone metabolism, including antiresorptive drugs, thyroid or parathyroid hormone, hormonal replacement therapy; or (3) they were pregnant or breastfeeding. A total of 180 patients were screened, and 78 fulfilled our criteria. As a comparative group, 72 age-matched (± 5 yrs) female patients with SLE fulfilling the recruiting criteria but with no organ damage (SDI = 0) were chosen. The study was approved by the Joint Chinese University of Hong Kong-New Territories East Cluster Clinical Research Ethics Committee. Written informed consent was obtained from all participants.
Data collection and clinical measurement
Demographic data collected from the study interview included age, body weight, body height, menstrual status (age and duration of menopause if applicable), smoking habit, aerobic exercise habits for at least 40 min per week, history of fracture after age 25 years, and family history of osteoporosis and fracture. Body mass index (BMI) was calculated for each participant. Disease duration, comorbidities, use of immunosuppressants (mycophenolate mofetil, azathioprine, cyclosporine A, and cyclophosphamide) and supplements (multivitamins and active form of vitamin D) in the preceding year, and use of GC (duration of use, current dose, cumulative dose, and average dose) during the whole disease course up to entry to the study were collected from medical records. Disease activity was assessed by a validated disease activity measure, the SLE Disease Activity Index (SLEDAI)13, and organ damage by the SDI12. It contains 24 descriptors in 9 organ systems, including clinical and laboratory variables, and measures disease activity of the preceding 10 days. The total SLEDAI score ranges from 0 (no activity) to 105 (maximum activity). The SDI is designed to assess organ damage due to the disease itself or the treatment in 12 organ systems occurring since the onset of SLE and presenting for at least 6 months. The total SDI score ranges from 0 (no damage) to 49 (maximum damage).
25-hydroxy vitamin D [25(OH)D] measurement
Serum levels of 25(OH)D were measured with a commercial ELISA kit (Immunodiagnostic Systems Inc.). All samples were collected at the time of the interview. Vitamin D deficiency was defined as serum level < 25 nmol/l (10 ng/ml), and insufficiency as < 75 nmol/l (30 ng/ml)14.
aBMD measured by DEXA
aBMD measures of left hip, lumbar spine (L1-L4, anteroposterior view), and nondominant ultradistal radius were performed using the same DEXA equipment in all participants (model 4500A; Hologic). Results were expressed in g/cm2 and T scores calculated with reference to general population norms of Hong Kong15.
Bone geometry, vBMD, and microarchitecture measurements
Bone geometry, vBMD, and microarchitecture were measured using a 3D HR-pQCT device (XtremeCT) at the nondominant distal radius (dominant distal radius was scanned in 1 patient because of a history of fracture at distal radius). This system uses a 2D detector array in combination with a 0.08-mm point-focus X-ray tube, enabling simultaneous acquisition of a stack of parallel CT slices with a nominal resolution (voxel size) of 82 μm. Details of image acquisition and analysis have been described7. The arm of the patient was fixed during the examination in an anatomically formed carbon fiber shell. An anteroposterior scout view was used to define the measurement region. A reference line was manually placed at the endplate of the radius. At the distal radius, 110 CT slices were obtained, thus delivering a 3D representation of roughly 9 mm in the axial direction. The first CT slice was 9.5 mm proximal to the reference line for the distal radius. The entire volume of interest was automatically separated into a cortical and trabecular region using a threshold-based algorithm. The threshold used to discriminate cortical from trabecular bone was set to one-third of the apparent cortical bone density value. Mean cortical thickness was defined as the mean cortical volume divided by the outer bone surface. The outcome variables used in our analyses included area for total bone, cortical, trabecular and estimated cross-sectional area; volumetric bone density (mg hydroxyapatite/cm3) for total average (D100), trabecular (Dtrab), and cortical bone (Dcomp); cortical thickness (Ct.Th, μm); and trabecular bone volume fraction (BV/TV, percentage), trabecular thickness (Tb.Th, mm), trabecular number (Tb.N, mm−1), and trabecular separation (Tb.Sp, mm). Derivations and definitions of trabecular variables were described in a previous study7. The precision of HR-pQCT at the distal radius demonstrates a coefficient of variation of 0.7% to 1.5% for density and 0.9% to 4.4% for trabecular architecture measurements7.
Finite-element analysis
The HR-pQCT images were then subjected to 3D finite-element analysis using Image Processing Language software (Scanco Medical) to derive further measures of bone biomechanical properties (bone strength). All images were segmented and only the voxels that represented bone tissue were converted to equal-size brick elements, resulting in models with 1 to 8 million elements. All elements were given a Young modulus of 10 GPa and a Poisson ratio of 0.316. A uniaxial compression in the longitudinal direction, which approximates to the common loading condition for Colle’s fracture of radius, was applied to the distal surface of the radius, to simulate a fall from standing height on the outstretched hand for determining bone biomechanical properties, i.e., axial stiffness and strength of radius. The variables used in our study included bone stiffness (kN/mm), apparent modulus (N/mm2), and estimated failure load (N), computed based on a criterion developed by Pistoia, et al17. With this criterion, it is assumed that failure will occur if the bone tissue equivalent strain exceeds 0.7% for > 2% of the bone tissue. Based on the results of the simulated compression tests at 1% strain, a scaling factor is determined such that the calculated bone tissue equivalent strain meets this criterion. The estimated failure load is then calculated by multiplying the total reaction force calculated from the compression test simulation by this scaling factor.
Statistical analysis
Analyses were performed using SPSS for Windows, version 13.0 (SPSS Inc., Chicago, IL, USA). Results were expressed as number (percentage) for categorical data, mean ± SD for normally distributed data, or median and interquartile range (IQR) for non-normally distributed data. Comparisons between continuous variables were performed by independent-samples t test and between categorical variables by chi-square test. Potential variables associated with bone quality and biomechanical measures were examined first by univariate tests and subsequently by multiple linear regression analysis. Tested variables included age, body height, body weight, BMI, menstrual status, disease duration, SLEDAI, SDI, duration of GC use, average GC dose, cumulative GC dose, use of immunosuppressants, and level of 25(OH)D. Univariate correlations for continuous variables were examined by Pearson’s or Spearman’s correlation for normally or non-normally distributed data, respectively. For categorical variables, Student’s t test or Mann-Whitney U test was used, depending on the data distribution. Univariate analysis was used as a screening tool to select the most potential variables in the multiple regression analysis. Variables with p < 0.05 in the univariate analysis were then entered into the multiple linear regression analysis to determine the independent explanatory factors for deterioration of bone quality.
RESULTS
Demographic and clinical characteristics and use of medications
SLE patients with organ damage had a mean age of 45.3 years, and over half were postmenopausal, with a mean duration of menopause of 4 years. Patients with damage had a median disease duration of 13.6 years and mild disease activity at the time of study entry. The median SDI score was 1, ranging from 1 to 5. Neuropsychiatric damage (25.6%), ocular damage (24.4%), pulmonary damage (14.1%), and musculoskeletal damage (14.1%) were the most common deficits comprising the SDI. Arthritis and nephritis were the most common features involving major organs in our cohort. Over half the patients with damage had arthritis (51.3%) and nephritis (66.7%). Serum level of 25(OH)D was low, with 98.7% of patients with damage diagnosed as having vitamin D insufficiency, but none as vitamin D deficiency. Table 2 shows the medication use of the study cohort. Patients with damage had a median duration of GC use of 10.4 years, with current dose 5 mg and cumulative dose 22.3 g. Over half the patients with damage had received multivitamins and a certain proportion were treated with immunosuppressant in the preceding year.
Demographic and clinical characteristics of the study cohort. Results are expressed as mean ± SD or median (interquartile range) unless otherwise indicated.
Use of medications by study cohort in the preceding year. Results are expressed as median (interquartile range) unless otherwise indicated.
Compared with patients with damage, SLE patients with no organ damage had significantly shorter disease duration (p = 0.01). Other variables related to demographic and clinical characteristics and use of medications were not significantly different between these 2 groups.
Bone density, microarchitecture, and biomechanical properties by DEXA and HR-pQCT
Table 3 gives data for aBMD measured by DEXA, bone geometry, vBMD, microarchitecture, and biomechanical properties assessed by HRCT and finite-element analysis in SLE patients with organ damage compared to patients without damage. Among patients with organ damage, 38 (48.7%) were osteopenic and 9 (11.5%) osteoporotic at either total hip or lumbar spine. The corresponding prevalence in patients without damage was 23 (31.9%) and 4 (5.6%), respectively. The prevalence of osteopenia at either total hip or lumbar spine was significantly higher in patients with organ damage (p = 0.037). Patients with damage had trends of deteriorated bone geometry, BMD, microarchitecture, and biomechanical properties compared to patients without damage, and the decrease in aBMD at total hip and trabecular thickness at distal radius achieved a significant difference.
Bone geometric, densitometric, microarchitectural, and biomechanical variables by the HR-pQCT. Results are expressed as mean ± SD.
Risk factors for deterioration of bone geometry, vBMD, microarchitecture, and biomechanical properties
Univariate analysis (Tables 4 and 5) was used to screen for selection of the most potential variables that might relate to deterioration of bone quality for the subsequent multiple regression analysis. There was no difference in bone quality between patients with and those without other damage such as neuropsychiatric, ocular, pulmonary, or musculoskeletal damage. As well, no difference of bone quality was found between patients with and those without treatment with immunosuppressant in the preceding year (data not shown).
Univariate analysis: correlation between HR-pQCT measures and demographic and clinical variables (continuous variables). Results are expressed as correlation coefficient (r).
Relationship between HR-pQCT measures and demographic and clinical variables (categorical variables). Results are expressed as mean ± SD.
Variables with p < 0.05, such as age, body height, body weight, SLEDAI, SDI, postmenopausal status, nephritis, and arthritis, were then entered into the multiple linear regression analysis to determine the independent explanatory factors for deterioration of bone quality (Tables 6A, 6B). During multiple linear regression, body height was found to be independently associated with bone geometry (total area, trabecular area, and cross-sectional area), while body weight was associated with cortical area, bone stiffness, and estimated failure load. The influences of age and postmenopausal status were limited to Tb.Sp and apparent modulus, respectively. Cumulative organ damage (SDI) was the only clinical variable that correlated consistently to deterioration of cortical area, vBMD (D100, Dcomp, Dtrab, Dmeta, and Dinn), and bone microarchitecture (Ct.Th, BV/TV, Tb.N, and Tb.Th). However, nephritis seemed to correlate positively with total area and trabecular area. SLEDAI or arthritis did not influence any bone quality during the multiple regression analysis (data not shown).
Multiple analysis for independent variables associated with deterioration of bone quality (HR-pQCT measures). Only significant variables are shown. Results are expressed as coefficients (ß) and 95% CI and p value.
Multiple analysis for independent variables associated with deterioration of bone quality (HR-pQCT measures). Only significant variables are shown. Results are expressed as coefficients (ß) and 95% CI and p value.
DISCUSSION
The use of high-resolution imaging techniques such as HR-pQCT allows us to investigate the microarchitectural and biomechanical properties of peripheral bone, which are emerging determinants of bone quality and fracture risks7. We previously showed that densitometric and microarchitectural measures assessed by HR-pQCT had greater discriminatory power for fracture in patients with SLE18. In this study, we evaluated bone geometry, vBMD, microarchitecture, and biomechanical properties by HR-pQCT in SLE patients with organ damage, compared with a group with no damage, and investigated factors correlated with bone quality assessed by HR-pQCT in SLE patients with organ damage. We found that increase in organ damage was independently correlated with densitometric and microarchitectural measures in patients with organ damage.
Bone loss and fracture in SLE are multifactorial. The traditional risk factors associated with reduced aBMD in SLE include female sex, older age, white or Asian race, lower BMI, early menopause, and a personal or maternal history of fracture and others. SLE-specific risk factors include active disease status, cumulative organ damage, sun avoidance inducing vitamin D deficiency, therapy-related premature menopause, and chronic GC use19,20. It has been documented that organ damage in SLE was significantly associated with reduced aBMD21,22. In our study, we selected SLE patients with at least 1 expression of organ damage, who were at a higher risk of bone loss. This is the first study to investigate the correlation between organ damage in SLE and vBMD, bone microarchitecture, and biomechanical properties. We found that increasing organ damage was the only clinical variable consistently correlated with deterioration of vBMD and bone microarchitecture after adjustment for other risk factors. The correlations between organ damage and vBMD and bone microarchitecture assessed by HR-pQCT demonstrate the potential use of this modality to further our knowledge of the underlying mechanism of bone deterioration in SLE.
Organ damage is a summary of longterm irreversible changes related to the disease itself and to related therapy12, both of which could be responsible for impairment of bone quality in SLE. As a systemic autoimmune inflammatory disorder, the chronic elevated levels of inflammatory cytokines in active SLE, such as tumor necrosis factor-α (TNF-α), interleukin 1 (IL-1), and IL-6, could increase the rate of bone resorption through activation of osteoclastogenesis23, leading to loss of bone density and deterioration of bone microarchitecture. The effect of SLE disease itself on bone quality is supported by the finding that premenopausal patients with early disease and no prior use of GC still had lower aBMD at the spine24. On the other hand, studies have shown that organ damage was one of the risk factors for lower aBMD among SLE patients who were on GC therapy5,21,22,25,26. GC is the mainstay of therapy for SLE, with the dose that is prescribed primarily depending on disease severity and distribution of organ damage. However, numerous adverse side effects may result from GC use. To evaluate the association between GC use and organ damage in patients with SLE, Petri, et al investigated patients with SLE from the Hopkins Lupus Cohort and found that cumulative GC dose was significantly associated with many manifestations of organ damage such as osteoporotic fracture, coronary artery disease, cataracts, and avascular necrosis; among these, osteoporotic fracture was the most strongly influenced27,28. The mechanism of GC-induced bone loss is mainly due to its transient increase of osteoclast lifespan and its continual induction of apoptosis of osteoblasts and osteocytes and suppression of osteoblastogenesis, leading to uncoupled resorption rate and formation rate29. Longterm GC therapy significantly and promptly increases bone loss and fracture risk30. Therefore, the predominant role of organ damage (by SDI) in the deterioration of bone quality might at least in part reflect the negative effect of GC on bone.
Few studies have found a relationship between disease activity or damage and aBMD. Dhillon and colleagues found no correlation between disease activity, measured by British Isles Lupus Assessment Group scores, and lumbar spine aBMD31. A similar finding was reported by Hansen, et al in a prospective study examining consecutive measurements of aBMD of 37 SLE patients, in which disease activity was measured by SLEDAI, and no correlation was found at baseline or at first and second year followup32. In our study, disease activity (by SLEDAI) had no significant association with bone quality or strength. This finding is not unexpected, since the majority of disease activity measures are limited to assessment of disease status during a short period of time, restricting its utility in predicting longterm outcomes such as bone quality, particularly in a condition such as SLE, which usually has a fluctuating disease course.
Photosensitivity is a common feature of SLE that is usually associated with avoidance of ultraviolet light, which decreases synthesis of vitamin D. Low serum 25(OH)D level predisposes the patient to secondary hyperparathyroidism, an additional risk for osteoporosis in SLE33. There is a high prevalence of hypovitaminosis D in South Asia34. The high prevalence of insufficiency of serum vitamin D in our cohort was comparable with that in women of child-bearing age in Hong Kong35, although our patients were older (mean age 45 vs 28 years). We did not find significant correlations between serum 25(OH)D and bone quality. This could be partly explained by the fact that 70% of our patients received vitamin D supplements (multivitamin or active vitamin D) in the preceding year. Disease activity, found to be inversely associated with level of vitamin D3, was relatively low in our cohort. Since patients were recruited over a period of almost a year, the influence of seasonal change of vitamin D level should be considered. Finally, the effects of GC therapy on vitamin D metabolism remain inconclusive30.
Our study has several limitations. First, it was a cross-sectional study. A longitudinal study is needed to determine the cause-effect relationship between organ damage and bone quality. Second, the sample size was relatively small and the number of patients with symptomatic fracture was small, and we did not have data related to asymptomatic vertebral fracture. Therefore, the correlation between organ damage and fracture could not be examined. Since there are no data related to precision of HR-pQCT in patients with SLE, the precision we noted was referred from the study of 256 healthy women by Boutroy, et al7. However, hand positioning and imaging analysis were acceptable in the patients with SLE. Serum levels of 25(OH)D were measured in different seasons because of the long period of patient recruitment, which limited the prediction ability and interpretation of the influence of 25(OH)D on bone. It would be more reliable to measure the mean 25(OH)D level during a certain time period to accurately determine its effect on a longterm cumulative variable such as bone quality. The lack of information for serum creatinine and glomerular filtration rate in the cohort at time of study entry prevented investigation of correlations between bone quality and renal function. Similarly, the lack of data for erythrocyte sedimentation rate prohibited study of the relationship between bone quality and inflammation. Finally, we could not demonstrate the respective contributions of disease-induced damage and GC therapy-induced damage on bone deterioration, because all patients were receiving GC; GC remains the mainstay of treatment for SLE in Hong Kong, thus a GC-naive group is likely to be a very small sample.
Cumulative organ damage was negatively correlated with vBMD and bone microarchitecture in female patients with SLE with organ damage receiving longterm GC. HR-pQCT provides insight into the underlying mechanism of bone loss in SLE that may be due to disease damage and partly to the chronic use of GC. Therapies without deleterious effect on bone that can effectively control disease activity and at the same time prevent organ damage are needed. Future longitudinal studies, and studies in particular including GC-naive SLE populations in multiple regions, are needed to monitor bone quality changes during the course of treatment in SLE and to distinguish the respective contribution of disease itself and GC therapy on bone deterioration.
Footnotes
-
Supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (Project no. CUHK471010).
- Accepted for publication July 3, 2012.