Article Text

Extended report
Prediction of improvement in skin fibrosis in diffuse cutaneous systemic sclerosis: a EUSTAR analysis
  1. Rucsandra Dobrota1,2,
  2. Britta Maurer1,
  3. Nicole Graf3,
  4. Suzana Jordan1,
  5. Carina Mihai2,
  6. Otylia Kowal-Bielecka4,
  7. Yannick Allanore5,
  8. Oliver Distler1
  9. on behalf of EUSTAR coauthors
    1. 1Division of Rheumatology, University Hospital Zurich, Zurich, Switzerland
    2. 2Department of Internal Medicine and Rheumatology, Cantacuzino Hospital, Carol Davila University of Medicine and Pharmacy, Bucharest, Romania
    3. 3Graf biostatistics, Winterthur, Switzerland
    4. 4Department of Rheumatology and Internal Medicine, Medical University of Bialystok, Bialystok, Poland
    5. 5Department of Rheumatology, University Paris Descartes and Cochin Hospital, Paris, France
    1. Correspondence to Dr Oliver Distler, Division of Rheumatology, University Hospital Zurich, Gloriastrasse 25, Zürich 8091, Switzerland; Oliver.Distler{at}

    Statistics from


    Systemic sclerosis (SSc) is a highly heterogeneous disease, making clinical management of SSc and patient selection for clinical trials challenging. For skin fibrosis, the modified Rodnan skin score (mRSS) is the most widely used measure in clinical practice and it is also the most frequent primary end point in clinical trials.1 ,2 Identifying predictors of change in mRSS over time is therefore of much interest for risk-tailored clinical management, as well as for clinical trial design to enrich for patients with worsening skin fibrosis.3

    In a recent, large-scale, observational study on the European Scleroderma Trials And Research (EUSTAR) database, we identified short disease duration (<15 months), joint synovitis and low baseline mRSS (>7/51 and ≤22/51) as independent predictors of skin progression in patients with diffuse cutaneous SSc (dcSSc).4 This provided an evidence-based tool for the improved identification of patients at risk for progressive skin involvement and is also valuable for cohort enrichment in clinical trials on skin fibrosis.

    While the identification of factors predicting progression has been improved, little is known about predictors of regression of skin fibrosis in patients with SSc. Regression of mRSS has long been identified as a characteristic feature of the natural history of skin fibrosis in patients with dcSSc. Regression is thought to occur after the early active phase of the disease has stabilised. However, the time to peak skin score varies widely in patients, leading to a highly heterogeneous pattern of regressing, progressing and stable patients even in early disease stages.5 ,6

    From a therapeutic perspective, there is general agreement that prevention of progression in patients with active skin fibrosis is one of the major treatment goals. However, it is much less established whether a therapeutic benefit can be achieved for patients who are already likely to show improvement of skin fibrosis under standard of care. In such a population the benefit/risk ratio of any treatment would have to be accurately assessed. Therefore, in order to identify those patients who would benefit most from therapeutic interventions, it is important to identify patients with progressive skin fibrosis, and to be aware of patients with skin regression. For trial design, it is important to enrich for progressive patients, and to exclude regressing patients under standard of care to increase the likelihood of identifying treatment effects.

    So far, previous attempts to identify predictors of change in mRSS have been largely inconclusive,4 and patients recruited for clinical studies targeting skin fibrosis often show spontaneous regression of mRSS.4 ,5 Thus, the objective of our study was to provide evidence-based predictors of skin improvement in patients with dcSSc using the EUSTAR database.


    Patients and study design

    The longitudinally followed EUSTAR cohort was analysed for this observational study. The whole EUSTAR data set, consisting of 12 274 patients at the time of the first data export (20 February 2015), was considered.

    The following inclusion criteria were used for cohort selection: fulfilment of American College of Rheumatology 1980 classification criteria, dcSSc, mRSS ≥7 at the first visit (baseline) and available data for mRSS at 12±2 months follow-up.

    Patients with dcSSc were identified according to LeRoy et al7 or, in case of missing values for the LeRoy criteria, by diffuse skin involvement at any visit. The minimum mRSS ≥7 was chosen because it reflects the lowest value classifiable as dcSSc, thus allowing the inclusion of patients with dcSSc with less severe to extensive skin fibrosis. The 1 year follow-up has been shown adequate for capturing significant changes in mRSS and is often used in clinical trials in skin fibrosis in SSc.2

    The clinical data in EUSTAR are prospectively collected in a multicentre approach following a standardised protocol.8–10 Regular training courses in skin scoring are organised by EUSTAR and all centres are advised to have the same examiner assessing the skin score in individual patients at follow-up visits.4 Quality indicators for data from the registry include regular external monitoring of large centres and plausibility checks on key items with written requests to centres for clarification. Ethics approval has been obtained from the respective local ethics committees by all participating EUSTAR centres.

    Statistical analysis

    The primary end point, improvement of skin fibrosis, was defined as a decrease in mRSS of >5 points and ≥25% within 1 year. These thresholds were chosen according to the minimal clinically important difference.11 Similarly, progression of skin fibrosis was defined as an increase in mRSS of >5 points and ≥25% within 1 year as used previously.4 For definitions of the clinical variables see the online supplement.

    A subanalysis using receiver operating characteristic analysis was performed with skin improvement, and, respectively, skin progression as the state variable, in order to explore the relationship between different baseline mRSS cut-off points and the proportion of regressors and progressors included in the cohort.

    Candidate predictors for skin improvement were selected based on nominal group technique by SSc experts (OD, YA, OK-B, CM, RD), who were asked to suggest clinically meaningful variables with face validity (see the online supplement). All parameters suggested by the experts (see online supplementary table S1) with <50% missing data were considered for the analysis. As a first step, a multivariable logistic regression model including all selected 11 parameters was run after single conditional mean imputation of the data. Baseline mRSS was centred at 7 points as all patients had baseline mRSS ≥7. Because baseline mRSS did not behave linearly, a quadratic term for baseline mRSS was included in the model. The Wald statistics (see online supplementary table S3) showed that the effects were very far from being significant (p value >0.7) for joint contractures and diffusing capacity of the lung for carbon monoxide ≥70%. Thus, these parameters were excluded from further models. The interaction between disease duration and baseline mRSS was also tested, but proved to be insignificant, meaning that the effect of baseline mRSS on regression of SSc did not depend on disease duration. The model after single imputation is shown in the online supplementary table S4. Single imputation was done for validation as it is not possible to validate models with multiply imputed data. Thus, bootstrap with 100 repetitions was used to validate the model (see online supplementary table S5). However, as multiple imputation provides more trustworthy estimates and ORs than single imputation, the final logistic regression model from the multiply imputed data set is presented.

    The statistical analysis was performed by the biostatistician (NG) using R V.3.1.0 (see the online supplement).


    Study population

    A total of 919 patients with dcSSc who met the inclusion criteria were analysed. Of these, 218/919 (24%) patients showed skin improvement over 1-year follow-up. The patients' demographic and clinical characteristics are shown in table 1.

    Table 1

    Description of the study cohort at baseline (total N=919)

    Multivariable analysis

    The candidate predictors selected after nominal group technique and exclusion of parameters with higher missing values and finally included in the multivariable analysis are shown in box 1.

    Box 1

    Candidate predictors of skin improvement selected for the analysis


    • Baseline mRSS

    • Disease duration

    • ANA positive

    • Anti-Scl70 positive

    • Tendon friction rubs

    • Proteinuria

    • Conduction blocks

    • Abnormal diastolic function

    • Fibrosis on chest X-ray

    • DLCO ≥70%

    • Joint contractures

    ANA, antinuclear antibodies; Anti-Scl70 antibodies, antitopoisomerase I antibodies; DLCO, diffusing capacity of the lung for carbon monoxide; mRSS, modified Rodnan skin score.

    The prediction model for skin improvement after single imputation is shown in online supplementary table S4. The most significant predictor was baseline mRSS (p<0.001). Other significant predictors were absence of tendon friction rubs and negative Scl-70 antibodies. The performance parameters of the model before and after validation are shown in online supplementary table S5.

    Multiple imputation was used to fit the model and to obtain SEs. The final model with multiply imputed data is shown in table 2.

    Table 2

    Final logistic regression model for prediction of skin improvement at 1 year

    High baseline mRSS remained the strongest predictor of skin improvement (p<0.001). The model for example indicates that the risk for regression within 12 months is more than doubled (OR=2.316) for a patient with a baseline mRSS of 22 points, in comparison to a patient with a baseline mRSS of 14 (all other parameters being equal). Furthermore, absence of tendon friction rubs significantly predicted skin improvement. Absence of anti-Scl-70 antibodies, which was also a significant predictor in the model after single imputation (see online supplementary table S3), only retained a trend in the model from the multiply imputed data set.

    Baseline mRSS as predictor of the pattern of skin change in dcSSc over 1 year

    The observation that high baseline mRSS was the strongest predictor of skin improvement complements our previous findings indicating low baseline mRSS as an important predictor of skin worsening.4 This was also confirmed in the current cohort: the 95/919 (10%) patients with dcSSc who showed skin progression within 1 year had lower baseline mRSS (p<0.001). Baseline mRSS is thus a predictor of change in skin score after 1 year, patients with lower skin scores being prone to progress and those with higher skin scores prone to improve within the next 12 months (see online supplementary figure S1).

    Having in mind the optimisation of cohort enrichment with maximal number of progressive patients and minimal numbers of regressive patients, we checked for the optimal mRSS cut-off (figure 1). In this cohort, an upper baseline mRSS cut-off value of 18 points performed best, including the highest proportion of progressors (78.9%) and the lowest proportion of regressors (35.3%, figure 1).

    Figure 1

    (A) Percentage of progressors and regressors per baseline modified Rodnan skin score (mRSS) range. Patients with lower skin score are more likely to progress, whereas those with higher skin scores are much more likely to regress. (B) Sensitivity for progression and regression depending on different cut-off values for baseline mRSS.

    If translated into clinical study design, this suggests a baseline mRSS between 7 and 18 as an inclusion criterion, raising questions for feasibility of recruiting patients. Thus, we next analysed the impact of higher cut-offs for upper baseline mRSS on the proportion of progressive and regressive patients. This analysis showed that a baseline mRSS between 18 and 25 would still allow identifying a reasonably high rate of progressors over regressors, whereas for skin scores higher than 25 a considerable drop of the included progressors and a dramatic increase in the percentage of regressors was observed (table 3).

    Table 3

    Proportion of progressors and regressors depending on different cut-offs for baseline mRSS


    Patients with improvement of skin fibrosis under standard of care are less likely to benefit from therapeutic interventions than patients with progressive skin fibrosis. In this large EUSTAR analysis of 919 patients with dcSSc with clinically derived real life data, we have identified parameters which can predict improvement of skin fibrosis over a 12 month observation period, the strongest being high baseline mRSS.

    To our knowledge, this is the first report of an evidence-based model for the prediction of skin improvement in a non-selected cohort of patients with dcSSc. In a previous study, Steen et al specifically focused on improvement of skin thickening in a cohort of 278 patients with early dcSSc. In this study, independent predictors for skin improvement could not be identified (except for D-Penicillamine use, which is however contradictory to the negative results of the dedicated randomised controlled trial).13–15 Potential explanations for this lack of predictors include the lower sample size, and the higher baseline mRSS in this study. The higher baseline mRSS might also explain the higher rate of improvers in this study (63% vs 22% in our EUSTAR analysis).

    A key message resulting from our study is the important role of baseline mRSS to predict either progression or regression of skin fibrosis. This also supports our previous EUSTAR analysis on worsening of skin fibrosis.4 In a pooled analysis of patients with dcSSc from seven multicentre clinical trials, mRSS at baseline had a weak negative correlation with any change in mRSS.6 Furthermore, in a recent analysis from the Canadian cohort, baseline mRSS was the only baseline parameter significantly associated with a significant change in mRSS at follow-up (defined as difference in 8 points).16 These data underline the value of mRSS for cohort enrichment in clinical trials. Our study also provides evidence-derived data on different thresholds of baseline mRSS and their performance to enrich for progressors over regressors in clinical trials (table 3). Thus, the optimal baseline mRSS cut-off for a specific study can be chosen from these data taking into account feasibility versus optimised cohort enrichment.

    Another important aspect to consider is the natural regression to the mean phenomenon: the more extreme the skin score values in the study population at baseline, the more likely they are to decrease towards the mean at follow-up, thus not necessarily reflecting treatment response. The regression to the mean is most likely the statistical effect that explains the selection of baseline mRSS as a strong predictor of worsening and regression, respectively.

    While our study addresses a very large SSc cohort with multiple data quality controls and external data monitoring, it also has several limitations. It has the natural drawbacks of registry data, such as missing data. However, we have addressed this issue by applying acknowledged imputation methods to compensate for missing data. Nonetheless, some of the candidate predictors had too much missing data which did not allow a trustable imputation, hence we could not include them in the analysis. The inclusion threshold of mRSS>7 aimed at identifying patients with true dcSSc with minimum involvement of distal upper extremities, but is, nonetheless, somewhat arbitrary. Our data also have to be confirmed in other cohorts with different baseline characteristics, for example, with higher prevalence of anti-RNA polymerase III antibodies. Additional information from in-between visits (eg, at 3 months, 6 months) as well as health assessment questionnaire data17 could bring additional valuable information on the course of skin fibrosis. Moreover, it has to be mentioned that the final prediction model only explains about 16% of the variation in skin fibrosis regression (see online supplementary table S4), indicating that other yet unknown factors such as, for example, biomarkers have an important role in determining the improvement of skin fibrosis in dcSSc. Further, modifications on how to measure the mRSS (eg, averaging vs maximising skin thickness over a certain skin area) will have great impact on the specific baseline mRSS values and have to be considered when these data are applied to clinical trials.

    In conclusion, our study provides novel evidence-based data for cohort enrichment in clinical trials on skin fibrosis in patients with dcSSc. These data further support a lower baseline mRSS as inclusion criteria to optimise the ratio of progressors over regressors for recruitment into clinical trials. Other significant predictors of improvement with potential application to clinical trials resulting from these data are absence of tendon friction rubs, and, potentially, negativity for anti-Scl70 antibodies.


    The authors thank Nicole Schneider for excellent administration and data entry to the EUSTAR cohort and Terkia Medkour for data management of the EUSTAR database.


    View Abstract

    Supplementary materials

    • Lay summary

      This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Linked Articles