Background and aims We studied damage accrual and factors determining development and progression of damage in an international cohort of systemic lupus erythematosus (SLE) patients.
Methods The Systemic Lupus International Collaborating Clinics (SLICC) Inception Cohort recruited patients within 15 months of developing four or more 1997 American College of Rheumatology (ACR) criteria for SLE; the SLICC/ACR damage index (SDI) was measured annually. We assessed relative rates of transition using maximum likelihood estimation in a multistate model. The Kaplan–Meier method estimated the probabilities for time to first increase in SDI score and Cox regression analysis was used to assess mortality.
Results We recruited 1722 patients; mean (SD) age 35.0 (13.4) years at cohort entry. Patients with damage at enrolment were more likely to have further worsening of SDI (SDI 0 vs ≥1; p<0.001). Age, USA African race/ethnicity, SLEDAI-2K score, steroid use and hypertension were associated with transition from no damage to damage, and increase(s) in pre-existing damage. Male gender (relative transition rates (95% CI) 1.48 (1.06 to 2.08)) and USA Caucasian race/ethnicity (1.63 (1.08 to 2.47)) were associated with SDI 0 to ≥1 transitions; Asian race/ethnicity patients had lower rates of new damage (0.60 (0.39 to 0.93)). Antimalarial use was associated with lower rates of increases in pre-existing damage (0.63 (0.44 to 0.89)). Damage was associated with future mortality (HR (95% CI) 1.46 (1.18 to 1.81) per SDI point).
Conclusions Damage in SLE predicts future damage accrual and mortality. We identified several potentially modifiable risk factors for damage accrual; an integrated strategy to address these may improve long-term outcomes.
- Systemic Lupus Erythematosus
- Outcomes research
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 3.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/3.0/
Statistics from Altmetric.com
Systemic lupus erythematosus (SLE) is a multisystem autoimmune disease in which adverse long-term outcomes remain a major challenge. In assessing patients with SLE, three disease dimensions are considered in formal outcomes studies: inflammatory disease activity, organ damage and health related quality of life (HRQOL).1 Damage is principally assessed using the Systemic Lupus International Collaborating Clinics (SLICC)/American College of Rheumatology (ACR) Damage Index (SDI), which has been extensively validated.2 ,3 SDI items represent irreversible damage that has occurred after the diagnosis of SLE. However, an item does not have to be attributable to lupus.2 As a general rule, items should be present for at least 6 months and once recorded in the SDI they are permanent such that the score cannot decrease.
The mean SDI tends to increase over time,4 and in time the majority of SLE patients will accrue damage. The SDI also predicts future mortality.5–7 It is therefore important to understand factors related to the development of damage. To date, studies have mainly focused on established SLE cohorts from a single centre or region.8–10 A number of factors have been associated with higher SDI scores, including older age at SLE onset,11 ,12 Hispanic and African ancestry race/ethnicity,8 ,13 ,14 chronic disease activity12 ,15 and major flares.16 Steroid exposure also predicts future damage, especially late-onset damage.8 ,16 ,17 There is also accumulating evidence that antimalarials (AMs) exert a protective role even after adjusting for their propensity for use in milder disease.18
We aimed to study damage accrual over the early years of follow-up in patients enrolled into the SLICC Inception Cohort. We focused on the rate of accrual of total damage as well as the impact of demographic, racial/ethnic and geographical variables. We also assessed the contribution of disease-related factors, therapy, co-morbidities and serological biomarkers to damage accrual. Finally, we sought to determine the relationship between damage and HRQOL, as well as future mortality.
SLICC Inception Cohort study
SLICC comprises 31 centres from 11 countries in North America, Latin America, Europe and Asia. An inception cohort was recruited from 2000 to 2011. Data were submitted to the co-ordinating centre at the University of Toronto at enrolment, and patients were reviewed annually. Laboratory tests necessary to evaluate disease activity and damage parameters were performed locally. The study was approved by the Institutional Research Ethics Boards of participating centres in accordance with the Declaration of Helsinki's guidelines for research in humans.
Patients and clinical assessments
Patients were enrolled within 15 months of recognition of four or more 1997 ACR classification criteria for SLE.19 We included patients who either had two study visits or had died after their first study visit, that is, patients who had two data points to model statistically. There were no specific exclusion criteria other than failing to meet four ACR criteria and it being >15 months since diagnosis. We noted demographic features including age, gender, race/ethnicity, geographical region and years in post-secondary education. We also noted the number of ACR criteria fulfilled by the baseline visit. At each visit we also assessed the Systemic Lupus Erythematosus Disease Activity Index 2000 (SLEDAI-2K)20 and the SDI.21 At each visit we recorded whether the patient was taking corticosteroids (yes or no). In addition to steroid use, we also recorded whether patients were taking AMs, immunosuppressives (ISs), or both AMs and ISs. Co-morbidities (recorded at each visit) included in our analysis were diabetes mellitus (physician confirmed diagnosis) and hypertension (systolic blood pressure >140 mm Hg and/or diastolic blood pressure >90 mm Hg and/or taking antihypertensive medications). Baseline serological markers included antibodies to double-stranded DNA and C3 and C4 complement (in local clinical laboratories at each centre). Antibodies to cardiolipin, β-2-glycoprotein I and the lupus anti-coagulant were measured at a central laboratory at the Oklahoma Medical Research Foundation as previously described.22 HRQOL was assessed using the Medical Outcomes Survey Short-Form 36 (SF-36). All patients provided written informed consent.
Simple descriptive statistics were used to summarise enrolment data. The SDI scores are discrete values that are observed over time for each patient. Thus, we used a multistate model for transitions among damage states, defined by the SDI scores. Specifically, at each visit a patient is assigned to one of the damage states according to their current SDI score. Since there are relatively few transitions to states 5–11, we merged these states into one state indicating an SDI score ≥5. We employed a multistate model with seven states, shown pictorially below:
Patients may only show deterioration in the SDI damage index over time and patients may die at any time during observation.
Multistate models allow estimation of the transition rates between these observed states, and these transition rates can be modelled as a function of explanatory variables (both time-independent and time-dependent). If λij(t) denotes the transition intensity from state i to state j at time t, then λij(t) may be modelled as follows: where λ0ij(t) denotes a baseline state i to state j transition intensity at time t, T denotes matrix transpose, and X(t) is a vector of explanatory variables with associated explanatory variable effects on the state i to state j intensity denoted by βij.
In this work, we assume constant baseline transition intensities for the relevant transitions. That is, λ0ij(t)=λ0ij for all t, for the state i to state j transitions. Maximum likelihood estimation is used to estimate the unknown model parameters (λ0ij and βij) for each transition in the model together with their associated SEs.
The correlation among the states of a patient at the different assessment visits is directly modelled through the Markov assumption that the future evolution of a patient's damage process depends only on his/her current state and not on his/her previous history. Initial modelling was based on a proportional hazards assumption and assumed common explanatory variable effects across selected transition rates. Notably we assumed that transition rates and explanatory variable effects between damage states where SDI ≥1 were the same and also that transitions to death from these higher SDI states were equal. Explanatory variable effects are reported as relative rates of transition, together with corresponding 95% CIs, obtained using maximum likelihood estimation. Age at diagnosis was standardised as (age in years−34.5)/13.4. For disease activity we report effects corresponding to 3-point increments in SLEDAI-2K.The Kaplan–Meier method was used to estimate the probabilities relating to time until first worsening of SDI score.
For the modelling of HRQOL outcomes, we fitted linear models using generalised estimating equations (GEEs) to account for the correlation among observations over time within each patient. Standard Cox regression analysis was performed with patient survival as the outcome and functions of damage scores over time as explanatory variables.
There were a total of 1722 patients in the SLICC Inception Cohort up to September 2011 and the mean (SD) number of follow-up visits was 4.25 (2.72). Demographic and disease-related factors at clinic entry are summarised in table 1. At baseline, 671 patients had their SDI scores recorded as they had more than 6 months of disease; of these 671 patients, 130 (19.4%) had an SDI score of 1 or more at the baseline visit. For this group (figure 1A), the overall estimate of the probability of SDI first worsening at a time greater than 6 years since clinic entry is approximately 0.58 (figure 1A). Put another way, of the 348 patients who were observed 6 years after clinic entry, 178 (51.1%) had at least one item of damage by that time point. When we stratified these 671 patients by baseline SDI score, those with initial damage were significantly more likely to have further worsening of the SDI at each follow-up visit (p<0.001) (figure 1B).
There were 1502 (1337 female and 165 male) patients who had at least two clinic visits or who had one clinic visit and subsequently died. In this group the predicted probability of remaining with the same damage score over a 5-year period was conditional on the pre-existing SDI score. The probability of death also increased with higher SDI scores (table 2).
Influence of age, gender, race/ethnicity and geographical region
Increasing age and male gender both had a significant influence on the probability of damage accrual. Higher standardised age increased the risk of future damage, especially in those with no current damage, and the influence of age was non-linear. The significant effect of (standardised age)2 suggests that the effect of increases in age is most marked in older patients. Therefore, assuming all other covariates have the same values for each age at diagnosis compared to a patient aged 35.4 years (the mean age of diagnosis), the relative transition rate was 1.58 for a patient aged 50 years. Compared to a 50-year-old, the relative transition rate was 2.51 for a 60-year-old, and for patients between 60 and 70 years of age the transition rate increased by a factor of 4.52. We also found that the effects of race/ethnicity and location of study sites were not independent (Pearson's χ2=2096.775, 16 df; p<0.001) (data on file) and that both had a significant impact on damage accrual. We therefore combined these into new variables (tables 3 and 4). Compared to Caucasians in Europe or Canada, USA patients of African ancestry had a higher risk of moving from no damage to damage and also of progressing from baseline damage to higher damage (relative transition rates (RTR) (95% CI) 1.99 (1.33 to 2.96) and 2.55 (1.92, 3.40), respectively), while Asians had lower transition rates (0.66 (0.43, 0.99)). Hispanic patients in Mexico also had a higher risk of progressing from baseline damage to higher damage (RTR (95% CI) 1.36 (1.02 to 1.83)) (tables 3 and 4).
Clinical, therapeutic and serological factors associated with development and/or progression of damage
Corticosteroid use, immunosuppressive use, SLEDAI-2K score and hypertension were all significantly associated with both the development of damage in patients free of damage at baseline, as well as progression of damage in patients with baseline damage (tables 3 and 4). The number of ACR criteria at enrolment was also associated with higher transitions from SDI 0 to ≥1 (RTR (95% CI) 1.19 (1.06 to 1.34)) as was IS use (1.69 (1.08, 2.63)). In addition, AM use was associated with a reduced transition rate to higher damage (0.60 (0.42, 0.84)) (tables 3 and 4).
Multivariate, multistate models for both transitions confirmed that age, USA patients of African ancestry, SLEDAI-2K score, steroid use and hypertension were predictive in both models of damage accrual (tables 3 and 4). For transition from SDI 0 to ≥1, male gender (RTR (95% CI) 1.48 (1.06 to 2.08)) and USA Caucasian race/ethnicity (RTR (95% CI) 1.63 (1.08 to 2.47)) were also associated with damage, while patients of Asian ethnicity had lower rates of transition (RTR (95% CI) 0.60 (0.39 to 0.93)). For transitions from SDI ≥1 to higher damage, patients taking AMs also had lower rates of transition (RTR (95% CI) 0.63 (0.44 to 0.89)). We found no evidence that baseline autoantibody status influenced damage accrual. We also noted a significant interaction between SLEDAI-2K and steroid use for transitions from SDI 0 to ≥1 (SLEDAI/3×(Corticosteroids = Yes); RTR (95% CI) 1.33 (1.02 to 1.74)). This suggests that the association between disease activity and transition to damage is greater for those patients taking corticosteroids.
In a secondary analysis we also assessed the influence of having ‘active renal disease’ during follow-up. We found that transition to higher damage states was greater in those with active renal disease (RTR (95% CI) SDI 0 to ≥1: 1.62 (1.10 to 2.38) and SDI ≥1 to higher damage: 1.66 (1.28 to 2.15), respectively).
Influence of SDI on HRQOL
The physical component domains of the SF-36 were more influenced by damage than the mental health components (figure 2). In a regression analysis with GEEs, the Physical Component Summary score (PCS) declined steadily with increased damage (table 5). We also found that PCS values were most influenced by recent changes in damage (coefficient (95% CI) −1.36 (−1.99 to −0.73) per SDI unit) and to a lesser extent by pre-existing damage (coefficient (95% CI) −0.39 (−0.36 to −0.58)).
Influence of SDI on mortality risk
To date there have been 41 deaths in the cohort. Using a Cox proportional hazards model with SDI score classed as a numerical variable rather than a factor, the SDI score was associated with an increased HR of 1.46 (95% CI 1.18 to 1.81) for mortality. A generalised likelihood ratio test of this model against a model with SDI score stratified by factors produced a test statistic of 14.25 (p=0.007 when compared to the quantiles of a χ2 distribution on 4 df). This suggests that the level of damage has a significant effect on mortality but there is not a simple relationship with SDI score (i.e. log-linear).
In a large international SLE inception cohort we have observed a steady accrual of damage over time and as has been reported by others, patients with damage are more likely to develop further damage over time and are also at higher risk of future mortality.5 ,7 ,9 ,14 We also found that damage has a significant effect on physical functioning. This steady accrual of damage has also been reported by other groups.10 ,14 ,23 ,24 We found no evidence of a plateau effect in our early cohort; studies that have suggested a plateau effect have tended to follow a cohort for more than 10 years.12 These observations are of clinical importance as the SDI is relatively easy to administer with some training in routine clinical settings and clearly identifies lupus patients at particularly increased risk of future adverse health outcomes.
Damage, especially a recent increase in the SDI, had a significant influence on physical functioning. Patients with recent damage may experience the maximal physical disability soon after acquiring the damage item. Over time, this may be ameliorated by physical adaptation or by corrective interventions. For example, patients with cataracts may later have lens replacement, patients with osteonecrosis may have joint arthroplasty, and patients with stroke are likely to rehabilitate over time.
Our data also allowed us to estimate the probabilities of developing future damage based on the patient's current SDI score. Such data allow us to consider how the SDI may be used as a clinical trial endpoint. For example, the estimated probability of remaining damage-free at 2 years is 0.844 if a patient has no damage at baseline. If there is one unit of damage, the estimated probability is 0.664 (similar for more than one). If we want to detect a doubling of the odds of remaining damage-free with an intervention, then in the first case the sample size needed (test at 5% and 80% power) would be approximately 670 (335 per group) and 349 (175 per group) in the second scenario.
We found a non-linear effect of age, with the effect of increases in age being most marked in older patients. Certain damage items such as cataracts, stroke and osteoporosis are, of course, more common with increasing age in the general population. Therefore, there may be a greater sensitivity to the additional effects of SLE and drug adverse effects with increasing age due to reduced organ reserve. There were also important differences among subsets of patients according to race/ethnicity and location. USA patients of African ancestry have an increased risk of damage accrual and USA Caucasians were also more likely to develop new damage. Of note, Asians had a lower risk of developing damage. There are a number of explanations for these findings, including differences in the clinical phenotype and/or its severity across different racial/ethnic groups. Response to therapy may also vary in different racial/ethnic groups25 ,26 and socio-economic factors and access to healthcare may also contribute. We used post-secondary education as a surrogate for socio-economic status, which was not significant in our models. Other more direct measures of socio-economic status were not collected but may have helped address this question more fully.
A number of similar factors drove development of new damage and/or progression of existing damage. Levels of disease activity, use of corticosteroids and hypertension all significantly influenced damage accrual. The significant interaction between disease activity and steroid therapy on new damage suggests that both act together to enhance the development of irreversible organ changes. Conversely, AMs were associated with reduced progression of damage, particularly in patients with baseline damage. These are all potentially modifiable risk factors. A multidimensional approach to damage prevention may therefore be needed and components of this would include better suppression of disease activity, minimising/avoiding corticosteroid use, use of AMs from an early stage and close control of hypertension. Also, if a novel therapy for SLE could achieve better disease control and steroid-sparing/avoidance, this ‘double benefit’ may translate to significant effects on damage accrual; indeed the interaction of inflammation and steroid therapy we found suggests there may be major gains in reducing future damage by such an approach.
Our study has a number of strengths. This is a large international inception cohort from diverse racial/ethnic backgrounds and geographical locations which has helped us understand how these factors influence damage development. We could also estimate probabilities of damage progression to help inform the use of the SDI as a clinical trial endpoint.
There are some limitations to our study. Patients were followed annually so it is difficult to fully model disease activity and therapeutic exposures. We also lack data on psychosocial factors which may influence damage progression. Sundaramurthy et al previously demonstrated that locus of control and time orientation were strong predictors of future damage.27 Our cohort, followed at a number of major teaching centres, may represent a lower estimate of damage accrual rates than that seen in general rheumatology practice; conversely, the tertiary referral case mix in many centres may act in the opposite way to influence our estimates. Finally, while our multivariate modelling suggests an independent effect of steroids and AMs on certain outcomes, we cannot exclude the possibility of residual confounding and that unmeasured factors may also influence the use of these agents in SLE patients.
In conclusion, we describe a steady increase in damage over time in SLE patients, with pre-existing damage being an important predictor of future damage accrual. We have also identified a number of modifiable risk factors that can be targeted as an integrated strategy. Overall, the SDI may therefore act in a way analogous to an erosion score in rheumatoid arthritis and could also act as a valid intermediate surrogate outcome for future mortality in SLE clinical trials.
We are grateful for the generous donation of our patients’ time and the dedication of all the fellows, research coordinators and research assistants in the SLICC network to the completion of this work.
Handling editor Tore K Kvien
Correction notice This article has been corrected since it was published Online First. Figure 1 has been corrected.
Funding INB is a National Institute for Health Research (NIHR) Senior Investigator and is supported by Arthritis Research UK, The Manchester Academic Health Science Centre, the NIHR Biomedical Research Unit Funding Scheme, the NIHR Manchester Wellcome Trust Clinical Research Facility and the Manchester Biomedical Research Centre. AGO'K and VF were supported by MRC (UK) funding U105261167. MBU was supported by the Canadian Institutes of Health Research (grant MOP-49529), the Lupus Foundation of Ontario, the Ontario Lupus Association, Lupus UK, the Lupus Foundation of America, the Lupus Alliance of Western New York, the Conn Smythe Foundation, the Lupus Flare Foundation, and the Tolfo family of Toronto, Ontario, Canada. S-CB was supported by the Korea Healthcare Technology R&D Project, Ministry for Health and Welfare, Republic of Korea (A120404). CG was supported by Lupus UK, the Arthritis Research Campaign, and the NIHR/Wellcome Trust Clinical Research Facility, Birmingham, UK. GSA was supported by the University of Alabama at Birmingham (NIH/NIAMS grant P60-AR-48095). AEC was supported by the Singer Family Fund for Lupus Research; she is also a Fonds de la Recherche en Sante du Quebec National Scholar. SB was supported by the Fonds de la Recherche en Sante du Quebec Jeune Chercheure and the McGill University Health Centre Research Institute; she is also a recipient of a Canadian Institutes of Health Research Junior Investigator Award and a Canadian Arthritis Network Scholar Award. DAI and AR were supported by the National Institute for Health Research, University College London Hospitals Biomedical Research Centre. PRF is a Distinguished Senior Investigator of The Arthritis Society with additional support from the Arthritis Centre of Excellence, University of Toronto. JGH was supported by a grant from the Canadian Institutes of Health Research (grant MOP-86526). DDG was supported by the Canadian Institutes of Health Research. MP was supported by the Hopkins Lupus Cohort (grant AR-43727) and the Johns Hopkins University General Clinical Research Center (grant M01-RR-00052). MAD was supported by the University of North Carolina Chapel Hill General Clinical Research Center (RR00046). ON and GKS were supported by grants from the Medical Faculty at Lund University, Swedish Combine Projects, the Crafoord Foundation and Lund University Hospital. RRG was supported by the NIH (grants UL 1RR 025741, K24-AR-02318 and P60-AR-48098). DLK was supported by the NIH (grants UL1 RR029882 and P60 AR062755). We acknowledge funding support to SLICC for this study from GlaxoSmithKline and Human Genome Sciences. The views expressed in this publication are those of the authors and not necessarily those of the NHS, the National Institute for Health Research or the Department of Health.
Competing interests None.
Ethics approval IRBs and Ethics Committees at all participating centres approved this study.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement The study data have been archived and are available on request. No additional unpublished data are available.