- Research article
- Open Access
Early trajectories of skin thickening are associated with severity and mortality in systemic sclerosis
Arthritis Research & Therapy volume 22, Article number: 30 (2020)
Systemic sclerosis (SSc) is a severe and highly heterogeneous disease. The modified Rodnan skin score (mRSS) is a widely used tool for the assessment of the extent and degree of skin thickness. This study aimed to identify the classes of patients with early similar skin thickening trajectories without any a priori assumptions and study their associations with organ involvement and survival.
From the French SSc national cohort, patients with a disease duration of less than 2 years at inclusion and with at least 2 mRSS available within the first 4 years of follow-up were enrolled. Classes of patients with similar mRSS trajectories were identified based on a latent class mixed model. The clinical characteristics and survival rate were compared between the obtained classes.
A total of 198 patients fulfilled the inclusion criteria, with a total of 641 mRSS available. The median disease duration and follow-up were 0.8 (interquartile range 0.4; 1.2) and 6.3 (3.8; 8.9) years, respectively. Individual trajectories of mRSS were highly heterogeneous between patients. Models with 1–6 latent classes of trajectories were sequentially assessed, and the 5-class model represented the best fit to data. Each class was characterized by a unique global trajectory of mRSS. The median disease duration did not differ significantly between classes. Baseline organ involvement was more frequent in classes with significant change over time (classes 2–5) than in class 1 (low baseline mRSS without significant change over time). Using Cox regression, we observed a progressively increasing risk of death from classes 1 to 5.
Early identification of clinical phenotype based on skin thickening trajectories could predict morbi-mortality in SSc. This study suggested that mRSS trajectories characterization might be pivotal for clinical practice and future trial designs.
Systemic sclerosis (SSc) is a chronic connective tissue disease characterized by widespread fibrosis of the skin and/or internal organs . Among the hallmarks of SSc, skin thickening is one of the pivotal symptoms and is used in routine practice to classify patients within the subsets—limited cutaneous SSc (lcSSc) and diffuse cutaneous SSc (dcSSc) . Modified Rodnan skin score (mRSS) is a semiquantitative score, ranging from 0 (normal) to 3 (severe), used to evaluate the skin thickness in 17 different cutaneous sites (for a total score from 0 to 51), and is correlated with histological skin thickness . mRSS is a validated clinical instrument  often used as primary or secondary outcomes in clinical trials .
Few studies have focused on the evolution of mRSS over time. Using latent trajectory modeling, Shand et al.  divided the early dcSSc patients (less than 2 years after disease onset) into three subgroups: “low baseline/improvers,” “high baseline/improvers,” and “high baseline/non-improvers.” Survival was associated with the subsequent mRSS trajectories, albeit only 68% of patients included could be included in one of these 3 groups. Perera et al.  described five skin thickening profiles in anti-topoisomerase I antibodies positive SSc patients early in the course of the disease. The 3 subgroups of dcSSc patients based on their skin thickening progression rate (STPR) were as follows: rapid (≥ 40 units per year), intermediate (15–40 units per year), and slow (≤ 15 units per year). The two subgroups of lcSSc patients were as follows: one where the skin clinical phenotype subsequently became diffuse, and the remaining one was limited throughout the follow-up. These studies underlined the well-known heterogeneity of the skin thickening evolution in SSc, and the complexity of the relations between mRSS at baseline and the skin thickening course (improvers/non-improvers and STPR). Some predictive factors of mRSS progression (defined as > 5 units and ≥ 25% increment in mRSS at 1 year follow-up) have been identified such as tendon friction rubs, joint synovitis, mRSS at baseline ≤ 22/51, disease duration < 15 months, and antibody status [8,9,10]. Nevertheless, it remains difficult to accurately predict the trajectory of mRSS in a given patient, which might limit the homogeneity and thus comparability of patients included in clinical trials . Deciphering the skin thickening heterogeneity is therefore of utmost importance considering the wide use of mRSS as the primary outcome and the three recent negative clinical trials to prove a benefice on skin thickening in SSc [11,12,13]. We herein aimed to identify the early mRSS longitudinal trajectories in SSc patients from the prospective French SSc national database without any a priori assumptions and to examine their associations with organ involvement and survival.
Patients and methods
Study population and inclusion criteria
The French SSc national database is a multicenter observational study conducted in 42 French hospital centers. Adult SSc patients under standard care were enrolled consecutively since 2010. Data were retrospectively collected before 2010 and then prospectively collected using a standardized form recorded on an online database (Clean Web®) with plausibility checks. In accordance with the French legislation, the database has received ethical approval from CCTIRS (approval no. 13.145; Advisory Committee on Information Processing in Material Research in the Field of Health). Data protection complied with the requirements of the National Information science and Liberties Commission and recorded under no. 914607. All included patients provided informed consent. Patients could also be included in other cohorts (e.g., European Scleroderma Trials and Research group) in a non-competitive way. Data were extracted in August 2015, and patients were eligible for the present analysis if (i) the ACR 1980 preliminary classification criteria  and/or 2013-ACR/EULAR SSc classification criteria  were fulfilled, (ii) inclusion visit occurred less than 2 years after the first onset of non-Raynaud phenomenon (RP) symptom, and (iii) the baseline mRSS and at least 1 mRSS during follow-up were available.
Data collection and variables
Baseline was defined as the date of inclusion in the database and follow-up as the time between the inclusion and the last available visit at the time of the extraction. Data on demographics, dates of first RP and first non-RP symptom, cutaneous subset, telangiectasia, calcinosis, and autoantibody status (anti-centromere [ACA], anti-topoisomerase I [ATA], anti-RNA polymerase III [anti-RNAP3], anti-U1 RNP, anti-PM/Scl, and other autoantibodies) were collected. Disease duration was defined as the time between the first non-RP symptom by patient report and inclusion visit. Organ involvements were defined by the occurrence of clinical events at baseline—skin involvement: mRSS and STPR defined as the mRSS at baseline visit divided by the disease duration (in years) [7, 16]; joint involvement: arthritis, arthralgia, friction rubs, or synovitis; muscle involvement: myalgia, myositis, or rhabdomyolysis; lung involvement: interstitial lung disease (ILD) diagnosed on high-resolution computerized tomography or chest X-ray, forced vital capacity (FVC % predicted value), and diffusing capacity of the lung for carbon monoxide (DLCO % predicted value); heart involvement: arrhythmia or conduction block or systolic dysfunction (left ventricular fraction ejection ≤ 45% of predicted value) or pericardial effusion; pulmonary hypertension (PH): mean pulmonary arterial pressure measured by right heart catheterization > 25 mmHg at rest; gastrointestinal tract (GIT) involvement: esophageal reflux, dysmotility, constipation, diarrhea, signs of bacterial overgrowth and/or malabsorption, abnormal esophageal manometry, and/or endoscopy test; digital ulcer (DU): history or active DU, digital tip, pitting scar, or digital ischemia; and scleroderma renal crisis (SRC): defined by new onset hypertension (≥ 150/85 mmHg) associated with a decrease in renal function defined by a decrement of at least 10% in the estimated glomerular filtration rate. C-reactive protein elevation was defined as a C-reactive protein level of > 6 mg/L. All immunosuppressive drugs were recorded during the follow-up. Death was also recorded.
The primary objective of this study was to delineate groups of patients according to their skin thickening trajectories (classes) as measured by mRSS over time using latent class mixed models (LCMM) [17, 18]. LCMM assumes that the population is divided into a finite number of groups called latent classes. Each latent class is characterized by a specific mean trajectory, which is described by a class-specific linear mixed model. In a given latent class, the individual trajectories of patients are close to each other, while the individual trajectories of patients of different classes tend to be dissimilar. In LCMM, latent classes correspond to an unknown categorical variable which is identified from data using a multinomial logistic model, and trajectories of mRSS are analyzed using the mixed model with random coefficients to take into account individual trajectories. The random effects (linear or quadratic) are determined from the analysis of residuals according to Verbeke and Molenberghs . To identify the number of classes, several LCMMs are performed. Each model predicts the shape of the trajectory of each class, estimates the probability for each individual of class membership and assigns each of them to the class for which the likelihood is the highest. Time 0 was defined by the date of baseline mRSS recorded. Trajectories were censored after 4 years of follow-up because of substantial missing records after this duration. To determine the best number of latent classes that represented the heterogeneity of developmental trajectories, we considered both formal statistical criteria (such as Bayes information criterion (BIC)) and model adequacy. A low value of BIC and average posterior probabilities of class membership greater than 0.7–0.8 correspond to the better model [19, 20]. A sensitivity analysis was performed with disease duration as an adjustment factor in the different LCMM.
Continuous variables were expressed as mean ± standard deviation (SD) or median and interquartile range (IQR) to describe classes, and as mean with 95% confidence interval (CI) to describe trajectory shapes. The normality of distribution was checked graphically and using the Shapiro-Wilk test. Categorical variables were expressed as frequencies and percentages. Comparisons of classes were performed using the analysis of variance or the Kruskal-Wallis test for quantitative variables and the Fisher’s exact test or the chi-square test for categorical variables. In the case of significant results, pairwise comparisons were performed, and a Bonferroni correction was applied. The survival rate in every class was estimated using the Kaplan-Meier method and compared using the Cox regression adjusted for age and sex. All statistical tests were performed at 2-tailed α level of 0.05. All data analyses were performed using the SAS software version 9.4 (SAS Institute Inc., Cary, NC, USA) and the R software (Package LCMM).
Of the 2063 patients in the database, 611 had a disease duration of less than 2 years at the time of inclusion (median disease duration (IQR): 0.7 (0.3; 1.2) years). Among them, 198 patients with an available baseline mRSS and at least 1 mRSS obtained within the first 4 years of follow-up were included (Fig. 1). The mean age of included patients was 51.1 ± 14.3 years, and majority of them were White patients; the male to female ratio was 3:1. The median disease duration at baseline and follow-up were 0.8 (IQR 0.4; 1.2) years and 6.3 (3.8; 8.9) years, respectively. The proportion of dcSSc patients was 49.7%. Nearly 95.0% of patients were positive for antinuclear antibodies (ACA 28.3%, ATA 55.9%, and anti-RNAP3 5.3%). The median baseline mRSS was 8 (2; 18) (Table 1). The proportion of male patients with dcSSc, ATA, and anti-RNAP3 was higher in the group included than in those with ≤ 1 mRSS (Additional file 1). No significant difference was observed in the mortality rate between the 2 groups (log-rank test, p = 0.40). A total of 641 mRSS values were available, and 55.0% of the patients had at least 3 mRSS records (Additional file 2).
Model fit evaluation
Individual trajectories of 198 patients included are presented in Fig. 2 and showed a notable heterogeneity between patients. Models with 1 to 6 latent classes were sequentially performed (Additional file 3). The 5-class model had the lowest value of BIC index, which suggests that it represented the best fit to data (Fig. 2, Additional file 4). The averages of posterior probabilities of belonging to a class (indicating that the modeled trajectories gathered individuals with similar patterns of skin change and distinguished the aforementioned individuals from those with dissimilar patterns of skin change) were 0.96, 0.88, 0.92, 0.95, and 0.93, respectively, for classes 1 to 5 (Additional file 5, Additional file 6). The median disease duration did not differ significantly between classes (p = 0.21). A sensitivity analysis with disease duration as an adjustment factor yielded similar results and confirmed that the 5-class model best fitted the data (Additional files 7, 8, and 9).
Demographics and clinical characteristics of the 5 mRSS trajectory classes (Table 2, Fig. 3, Additional file 10)
Class 1 was characterized by a low baseline mRSS (mean mRSS 4.1 [95% CI 3.2; 5.0]) with no significant change over time (mean mRSS at 1 year, 2 years, 3 years, and 4 years: 4.6 [3.3; 5.9], 5.1 [3.7; 6.4], 5.6 [4.3; 6.9], and 6.2 [3.8; 8.6], respectively). This class included 117 patients with lcSSc (82.6%), which primarily affects White women with ACA (42.3%) or ATA (42.3%). Almost all ACA-positive patients (95%) were assigned to this class. At baseline, one-third of the patients had ILD, joint, GIT, and DU involvements. The median STPR was 3.9 (IQR 1.2; 9.3) units/year.
Class 2 slightly improved from a mean baseline mRSS of 20.8 [95% CI 19.0; 22.5] to a mean mRSS at 1 year, 2 years, 3 years, and 4 years of 18.4 [16.3; 20.5], 16.8 [14.5; 19.0], 15.7 [13.6; 17.9], and 15.4 [11.0; 19.8], respectively. This class comprised 43 patients composed of White (87.5%) women (65.1%) with dcSSc (97.6%) associated with ATA (75.0%). Joint, DU, GIT, and ILD involvements were common at baseline. SRC was found in 7 patients (22.6%). The median STPR was 21.8 (IQR 16.0; 35.9) units/year.
Class 3 was characterized by a 2-step trajectory with a low baseline mRSS (mean mRSS: 8.7 [95% CI 6.0; 11.5]) rapidly increasing to a mean estimated peak mRSS of 23.2 [18.8; 27.6] at 2.3 years of follow-up, then followed by an improvement (mean mRSS at 4 years: 14.5 [8.4; 20.7]). Three of them were Black patients, 11 had dcSSc, and 2 had lcSSc. ATA was common. More than two-thirds of them had joint, DU, GIT, and ILD involvements. The median STPR was 7.8 (IQR 5.7; 13.4) units/year.
Class 4 was characterized by a 2-step trajectory with a mean baseline mRSS of 25.1 [95% CI 22.6; 27.6], which is rapidly increasing to a mean estimated peak mRSS of 41.6 [37.2; 46.0] at 2.2 years of follow-up, then followed by an improvement (mean mRSS at 4 years: 29.5 [22.7; 36.2]). This class was composed of 13 patients including 6 men and 3 Black patients. ATA, joint, GIT, and ILD involvements were frequent at baseline. The median STPR was 38.5 (IQR: 29.0; 131.3) units/year.
Class 5 was characterized by a mean baseline mRSS of 35.1 [95% CI 32.2; 37.9] subsequently improving (mean mRSS at 1 year, 2 years, 3 years, and 4 years: 27.6 [23.3; 31.8], 21.5 [17.2; 25.7], 16.8 [10.1; 23.5], and 13.5 [0; 29.5], respectively). All 12 patients had dcSSc mainly associated with ATA. At baseline, most of them had joint and DU involvements. A third was affected by ILD for whom the median FVC and DLCO were 54.0% (IQR: 53.0; 92.0) and 44.0% (41.0; 53.0), respectively (Additional file 10). The median STPR was 34.4 (IQR 27.6; 74.9) units/year.
Steroids and immunosuppressive treatments were more frequently used in classes 2 to 5 than in class 1 (from 82.5% to 100% versus 53.3% in class 1; p < .001) during the follow-up. No significant difference was noted between classes 2 to 5 in terms of steroids, methotrexate, azathioprine, cyclophosphamide, and rituximab. Mycophenolate mofetil (MMF) was used more often in class 5 (90.0%) than in classes 2 (44.1%; p = 0.013) and 3 (41.7%; p = 0.031). No significant difference was found between classes 4 (72.7%) and 5 (90.0%).
Kaplan-Meier curves are shown in Fig. 4. Survival was different according to the trajectory classes (p = 0.025). Using Cox regression analysis (Additional file 11), we observed a progressive increase in the risk of death from classes 2 to 5 compared with class 1 (reference): hazard ratio (HR) (class 2) = 1.35 [95% CI 0.33; 5.46], HR (class 3) = 2.99 [0.74; 12.07], HR (class 4) = 4.05 [1.09; 15.13], and HR (class 5) = 5.85 [1.63; 21.03]. The results were similar after adjusting for age and sex (Fig. 4, Additional file 12).
The main results were as follows: (i) LCMM identified without any a priori assumptions five distinct trajectories of mRSS during the follow-up in early SSc patients under standard care (inclusion within 2 years of the first non-RP symptom), (ii) the mRSS trajectory classes were associated with different organ involvement and survival, and (iii) patients with high baseline mRSS or those who peaked after baseline had the highest severity.
The natural evolution of skin thickening is very heterogeneous, yet it is generally accepted that it tends to worsen at the beginning of dcSSc to a maximum, which usually occurs over the first 2–3 years after disease onset, then followed by improvement at the advanced stage [10, 21,22,23]. However, the data related to this are limited and complex. Modeling the evolution of mRSS has only been performed in a few studies. Shand et al.  classified 131/192 (68%) early dcSSc patients into 3 subgroups using latent trajectory modeling over the first 2 years of follow-up. Those three subgroups, “low baseline (mean mRSS: 20 ± 6)/improvers,” “high baseline (mean mRSS: 42 ± 8)/non-improvers,” and “high baseline (mean mRSS 35 ± 7)/improvers,” had similar trajectories than classes 2, 4, and 5. The survival rate in the “high baseline/non-improvers” (similar to the class 4 in our study) subgroup was significantly worse than that in the 2 other subgroups. In our study, we observed that classes 4 and 5 had the worst survival. Moreover, we found 2 additional trajectories: one mainly composed of lcSSc (class 1), and another one characterized by a 2-step trajectory with a low baseline mRSS rapidly increasing to a mean estimated peak mRSS of 23.2 [95% CI 18.8; 27.6] before improving (class 3). We may have captured this last trajectory due to the following reasons: (i) we included patients with lcSSc, and (ii) the very short disease duration and the longer follow-up (4 years) in our study allowed us to discriminate these individual trajectories with different patterns of early skin change. Using a prespecified definition of progressive skin disease (increase in mRSS of > 5 points and ≥ 25% from baseline), Maurer et al.  identified several independent factors associated with skin thickening progression: baseline mRSS of ≤ 22/51, low baseline STPR, and disease duration of ≤ 15 months. The best prediction model of worsening performed correctly in only 44.4% of the cases, suggesting that it is still not easy to accurately predict the skin thickening progression. To our knowledge, no study has attempted to model the evolution of mRSS without any a priori assumptions or prespecified definition.
Our approach identified five original distinct trajectories over time meeting the fittest formal statistical criteria, model adequacy, and clinical relevance of discriminated trajectories. These five classes of trajectories were distinguished by intrinsic characteristics such as baseline mRSS, trajectory slopes of worsening or improvement, and mRSS peaks. The 5-class model remained the best model to decipher the global heterogeneity with disease duration as an adjustment factor. Class 1 had low mRSS and STPR values at baseline, less of organ involvements at baseline, and better survival, as reported in lcSSc [24,25,26]. Most patients with ACA (95%) were assigned to class 1. However, 42% of patients in class 1 had ATA, which is a higher proportion than that usually observed in lcSSc [24, 27, 28]. In addition to lcSSc, there were 20 patients (17%) with dcSSc associated with ATA and whose median baseline mRSS was 9 (IQR 6.5; 9.5). Low mRSS values have already been reported in dcSSc [29,30,31] and some limitations on the current classification have been identified [28, 32], especially when the forearms are involved. Their assignation to class 1 was probably related to the modeled trajectory shape of class 1 that was fitter to their individual trajectories than the other modeled trajectories classes.
Organ involvement was more frequent in classes 2 to 5 than in class 1. As expected, the use of immunosuppressive drugs was also more frequent in classes 2 to 5 than in class 1. With regard to survival, we found an increasing risk of death from classes 2 to 5 compared with class 1, especially in trajectory classes with a high baseline mRSS (class 5) or an intermediate baseline mRSS that peaked after baseline (class 4). These 2 classes also shared a rapid STPR at baseline (median > 30 units/years) compared with other classes. In our study, patients with ATA were present in all five trajectory classes while we noted that lcSSc patients were predominantly in class 1 and those with dcSSc were mainly in classes 2 to 5. This heterogeneity is well known and has been confirmed by several recent works highlighting the importance of autoantibodies, cutaneous subset, and disease duration to predict organ complications and survival [7, 23, 28]. In particular, it appears relevant to distinguish especially ATA-associated lcSSc from ATA-associated dcSSc [33, 34]. Recently, Wu et al. reported that in dcSSc, skin progression that occurred within 1 year was independently associated with forced vital capacity decline and all-cause death , and Zheng et al. reported worst disease outcomes in early dcSSc patients with worsening skin score (≤ 3 years), whose data were extracted from the Canadian Scleroderma Research Group database . Taken together, these results highlight the need to consider the baseline mRSS and the early changes in skin thickening (worsening/improving) in patient risk stratification.
Our study has the following strengths: we used prospective data from a French multicenter cohort on SSc including patients with early disease duration from the first non-RP symptom (≤ 2 years). We were able to identify five trajectory classes using an approach without any a priori assumptions. Moreover, the mRSS changes that occurred in the trajectory classes 2 to 5 seemed clinically relevant. Indeed, they were higher than the minimal clinically important difference in mRSS (≥ 3–4 units) estimated in the Scleroderma Lung Studies . Finally, we showed under standard care that patients with high baseline mRSS or those whose mRSS peaked after baseline had the highest severity.
The potential limitations include firstly potential biases in relation to inclusion criteria. Skin involvement, which is an important outcome in dcSSc, could have been more frequently assessed in dcSSc patients explaining a higher proportion of men with dcSSc, ATA, and anti-RNAP3 among included patients than among excluded patients. However, mortality did not seem different between the two groups (Additional file 1). We cannot exclude having missed out the most severe patients who would have died before the second record. Second, we chose to include patients with a disease duration of ≤ 2 years. A longer disease duration might have enabled a larger study size but would have increased the proportion of patients in which the mRSS has already been reached (usually during the first 2–3 years of the disease) . A shorter disease duration could also have been more relevant but would have reduced the sample size preventing a robust modeling. Third, the disease duration appeared longer in class 5 compared to class 4, but modeling with disease duration as adjustment factor confirmed that the 5-class model outperformed the other models, suggesting that the 5-class model remained as the optimal model. The small size of some classes (reflecting that patients assigned to them were rarely encountered in clinical practice) and the unknown causes of death may limit the interpretation of data. Fourth, mRSS was not recorded by the same practitioner for each patient leading to potential inter-observer variability. The long duration of the study and the turnover of the investigators could have led to a high variability in mRSS. However, all participating centers are SSc referral centers in which physicians attend regular mRSS assessment training, which might have reduced this limitation [38, 39]. Fifth, the start and end dates for each immunosuppressive drug in parallel to mRSS assessments were lacking in the database. Therefore, the influence of immunosuppression could not be considered for the trajectory modeling. Finally, our findings should be confirmed on an external validation cohort.
This study identified five distinct mRSS trajectories in early SSc patients. The mRSS trajectory classes were associated with organ involvement and survival. Early identification of clinical phenotype based on skin thickening trajectories could predict morbi-mortality in SSc and influence clinical management.
Availability of data and materials
The data that support the findings of this study are available on separate scientific request (contact Prof. Eric Hachulla, Department of Internal Medicine, University of Lille, Lille, France; firstname.lastname@example.org), but restrictions apply to the availability of these data, which were used under license for the current study and, hence, are not publicly available.
Anti-RNA polymerase III antibodies
Anti-topoisomerase I antibodies
Bayes information criterion
Diffuse systemic sclerosis
Diffusing capacity of the lung for carbon monoxide (% predicted value)
Forced vital capacity (% predicted value)
Interstitial lung disease
Latent class mixed model
Limited cutaneous systemic sclerosis
Modified Rodnan skin score
Scleroderma renal crisis
Skin thickening progression rate
Hachulla E, Launay D. Diagnosis and classification of systemic sclerosis. Clin Rev Allergy Immunol. 2011;40(2):78–83. https://doi.org/10.1007/s12016-010-8198-y.
LeRoy EC, Medsger TA. Criteria for the classification of early systemic sclerosis. J Rheumatol. 2001;28(7):1573–6.
Furst DE, Clements PJ, Steen VD, Medsger TA, Masi AT, D'Angelo WA, et al. The modified Rodnan skin score is an accurate reflection of skin biopsy thickness in systemic sclerosis. J Rheumatol. 1998;25(1):84–8.
Kumánovics G, Péntek M, Bae S, Opris D, Khanna D, Furst DE, et al. Assessment of skin involvement in systemic sclerosis. Rheumatology (Oxford). 2017;56(suppl_5):v53–66. https://doi.org/10.1093/rheumatology/kex202.
Khanna D, Furst DE, Clements PJ, Allanore Y, Baron M, Czirjak L, et al. Standardization of the modified Rodnan skin score for use in clinical trials of systemic sclerosis. J Scleroderma Relat Disord. 2017;2(1):11–8. https://doi.org/10.5301/jsrd.5000231.
Shand L, Lunt M, Nihtyanova S, Hoseini M, Silman A, Black CM, et al. Relationship between change in skin score and disease outcome in diffuse cutaneous systemic sclerosis: application of a latent linear trajectory model. Arthritis Rheum. 2007;56(7):2422–31.
Perera A, Fertig N, Lucas M, Rodriguez-Reyna TS, Hu P, Steen VD, et al. Clinical subsets, skin thickness progression rate, and serum antibody levels in systemic sclerosis patients with anti–topoisomerase I antibody. Arthritis Rheum. 2007;56(8):2740–6.
Avouac J, Walker UA, Hachulla E, Riemekasten G, Cuomo G, Carreira PE, et al. Joint and tendon involvement predict disease progression in systemic sclerosis: a EUSTAR prospective study. Ann Rheum Dis. 2016;75:103–9. https://doi.org/10.1136/annrheumdis-2014-205295.
Maurer B, Graf N, Michel BA, Müller-Ladner U, Czirják L, Denton CP, et al. Prediction of worsening of skin fibrosis in patients with diffuse cutaneous systemic sclerosis using the EUSTAR database. Ann Rheum Dis. 2015;74(6):1124–31. https://doi.org/10.1136/annrheumdis-2014-205226.
Herrick AL, Peytrignet S, Lunt M, Pan X, Hesselstrand R, Mouthon L, et al. Patterns and predictors of skin score change in early diffuse systemic sclerosis from the European Scleroderma Observational Study. Ann Rheum Dis. 2018;77(4):563–70. https://doi.org/10.1136/annrheumdis-2017-211912.
Distler O, Allanore Y, Denton C, Kuwana M, Matucci-Cerinic M, Pope J, et al. Riociguat in Patients with Early Diffuse Cutaneous Systemic Sclerosis: a Randomized, Double-Blind, Placebo-Controlled Phase IIb Study (RISE-SSc). Arthritis Rheumatol. 2018;70(suppl 10) https://acrabstracts.org/abstract/riociguat-in-patients-with-early-diffuse-cutaneous-systemic-sclerosis-a-randomized-double-blind-placebo-controlled-phase-iib-study-rise-ssc/. Accessed 9 Nov 2018.
Khanna D, Lin C, Kuwana M, Allanore Y, Batalov A, Butrimiene I, et al. Efficacy and safety of tocilizumab for the treatment of systemic sclerosis: results from a phase 3 randomized controlled trial [abstract]. Arthritis Rheumatol. 2018;70(suppl 10) https://acrabstracts.org/abstract/efficacy-and-safety-of-tocilizumab-for-the-treatment-of-systemic-sclerosis-results-from-a-phase-3-randomized-controlled-trial/. Accessed 9 Nov 2018.
Khanna D, Spino C, Johnson S, Chung L, Whitfield ML, et al. Abatacept in early diffuse cutaneous systemic sclerosis: results of a phase II investigator-initiated, multicenter, double-blind, randomized, placebo-controlled trial. Arthritis Rheumatol. 2020;72(1):125–36. https://doi.org/10.1002/art.41055 Epub 2019 Dec 10.
Masi AT. Subcommittee For Scleroderma Criteria of the American Rheumatism Association Diagnostic and Therapeutic Criteria Committee. Preliminary criteria for the classification of systemic sclerosis (scleroderma). Arthritis Rheum. 1980;23:581–90. https://doi.org/10.1002/art.1780230510.
Hoogen F, Khanna D, Fransen J, Johnson SR, Baron M, Tyndall A, et al. 2013 classification criteria for systemic sclerosis: an American College of Rheumatology/European League Against Rheumatism Collaborative Initiative. Arthritis Rheum. 2013;65(11):2737–47. https://doi.org/10.1002/art.38098.
Domsic RT, Rodriguez-Reyna T, Lucas M, Fertig N, Medsger TA. Skin thickness progression rate: a predictor of mortality and early internal organ involvement in diffuse scleroderma. Ann Rheum Dis. 2010;70:104–9. https://doi.org/10.1136/ard.2009.127621.
McCulloch CE, Lin H, Slate EH, Turnbull BW. Discovering subpopulation structure with latent class mixed models. Stat Med. 2002;21(3):417–29.
Verbeke G, Molenberghs G. Linear mixed models for longitudinal data. New-York: Springer-Verlag; 2000. https://doi.org/10.1007/978-1-4419-0300-6.
Schreiber JB. Latent class analysis: an example for reporting results. Res Social Adm Pharm. 2017;13(6):1196–201. https://doi.org/10.1016/j.sapharm.2016.11.011.
Van de Schoot R, Sijbrandij M, Winter SD, Depaoli S, Vermunt JK, van de Schoot R, et al. The GRoLTS-Checklist: guidelines for reporting on latent trajectory studies. Struct Equ Model Multidiscip J. 2017;24:451–67. https://doi.org/10.1080/10705511.2016.1247646.
Clements P, Lachenbruch P, Furst D, Paulus H. The course of skin involvement in systemic sclerosis over three years in a trial of chlorambucil versus placebo. Arthritis Rheum. 1993;36(11):1575–9.
Medsger TA Jr. Natural history of systemic sclerosis and the assessment of disease activity, severity, functional status, and psychologic well-being. Rheum Dis Clin North Am. 2003;29(2):255–73. Elsevier. https://doi.org/10.1016/S0889-857X(03)00023-1.
Nihtyanova SI, Sari A, Harvey JC, Leslie A, Derrett-Smith EC, et al. Using autoantibodies and cutaneous subset to develop outcome-based disease classification in systemic sclerosis. Arthritis Rheumatol. 2019. https://doi.org/10.1002/art.41153 [Epub ahead of print].
Walker UA, Tyndall A, Czirják L, Denton C, Farge-Bancel D, et al. Clinical risk assessment of organ manifestations in systemic sclerosis: a report from the EULAR scleroderma trials and research group database. Ann Rheum Dis. 2007;66(6):754–63. https://doi.org/10.1136/ard.2006.062901.
Tyndall AJ, Bannert B, Vonk M, Airo P, Cozzi F, et al. Causes and risk factors for death in systemic sclerosis: a study from the EULAR scleroderma trials and research (EUSTAR) database. Ann Rheum Dis. 2010;69(10):1809–15. https://doi.org/10.1136/ard.2009.114264.
Frantz C, Huscher D, Avouac J, Hachulla E, Balbir-Gurman A, et al. Outcomes of limited cutaneous systemic sclerosis patients: results on more than 12,000 patients from the EUSTAR database. Autoimmun Rev. 2020;19(2):102452. https://doi.org/10.1016/j.autrev.2019.102452.
Meier FMP, Frommer KW, Dinser R, Walker UA, Czirjak L, et al. Update on the profile of the EUSTAR cohort: an analysis of the EULAR scleroderma trials and research group database. Ann Rheum Dis. 2012;71(8):1355–60. https://doi.org/10.1136/annrheumdis-2011-200742.
Sobanski V, Giovannelli J, Allanore Y, Riemekasten G, Cozzi F, Distler O, et al. Phenotypes determined by cluster analysis and their survival in the prospective EUSTAR cohort of patients with systemic sclerosis. Arthritis Rheumatol. 2019;71(9):1553–70. https://doi.org/10.1002/art.40906.
Herrick AL, Pan X, Peytrignet S, Lunt M, Hesselstrand R, Mouthon L, et al. Treatment outcome in early diffuse cutaneous systemic sclerosis: the European Scleroderma Observational Study (ESOS). Ann Rheum Dis. 2017;76(7):1207–18. https://doi.org/10.1136/annrheumdis-2016-210503.
Merkel PA, Silliman NP, Clements PJ, Denton CP, Furst DE, Mayes MD, et al. Patterns and predictors of change in outcome measures in clinical trials in scleroderma: an individual patient meta-analysis of 629 subjects with diffuse cutaneous systemic sclerosis. Arthritis Rheum. 2012;64(10):3420–9. https://doi.org/10.1002/art.34427.
Amjadi S, Maranian P, Furst DE, Clements PJ, Wong WK, Postlethwaite AE, et al. Course of the modified Rodnan skin thickness score in systemic sclerosis clinical trials: analysis of three large multicenter, double-blind, randomized controlled trials. Arthritis Rheum. 2009;60(8):2490–8. https://doi.org/10.1002/art.24681.
Johnson SR, Soowamber ML, Fransen J, Khanna D, Van Den Hoogen F, Baron M, et al. There is a need for new systemic sclerosis subset criteria. A content analytic approach. Scand J Rheumatol. 2018;47:62–70. https://doi.org/10.1080/03009742.2017.1299793.
Srivastava N, Hudson M, Tatibouet S, Wang M, Baron M, et al. Thinking outside the box--the associations with cutaneous involvement and autoantibody status in systemic sclerosis are not always what we expect. Semin Arthritis Rheum. 2015;45(2):184–9. https://doi.org/10.1016/j.semarthrit.2015.04.009.
Kranenburg P, van den Hombergh WMT, Knaapen-Hans HKA, van den Hoogen FHJ, Fransen J, et al. Survival and organ involvement in patients with limited cutaneous systemic sclerosis and anti-topoisomerase-I antibodies: determined by skin subtype or auto-antibody subtype? A long-term follow-up study. Rheumatology (Oxford). 2016;55(11):2001–8. https://doi.org/10.1093/rheumatology/kew298.
Wu W, Jordan S, Graf N, de Oliveira PJ, Curram J, Allanore Y, et al. Progressive skin fibrosis is associated with a decline in lung function and worse survival in patients with diffuse cutaneous systemic sclerosis in the European Scleroderma Trials and Research (EUSTAR) cohort. Ann Rheum Dis. 2019;78(5):648–56. https://doi.org/10.1136/annrheumdis-2018-213455.
Zheng B, Nevskaya T, Baxter CA, Ramey DR, Pope JE, et al. Changes in skin score in early diffuse cutaneous systemic sclerosis are associated with changes in global disease severity. Rheumatology (Oxford). 2019. https://doi.org/10.1093/rheumatology/kez299 [Epub ahead of print].
Khanna D, Clements PJ, Volkmann ER, Wilhalme H, Tseng C-H, Furst DE, et al. Minimal clinically important differences for the modified Rodnan skin score: results from the scleroderma lung studies (SLS-I and SLS-II). Arthritis Res Ther. 2019;21(1):23. https://doi.org/10.1186/s13075-019-1809-y.
Czirják L, Nagy Z, Aringer M, Riemekasten G, Matucci-Cerinic M, Furst DE, et al. The EUSTAR model for teaching and implementing the modified Rodnan skin score in systemic sclerosis. Ann Rheum Dis. 2007;66(7):966–9. https://doi.org/10.1136/ard.2006.066530.
Park JW, Ahn GY, Kim J-W, Park ES, Kang J-H, Chang SH, et al. Impact of EUSTAR standardized training on accuracy of modified Rodnan skin score in patients with systemic sclerosis. Int J Rheum Dis. 2019;22(1):96–102. https://doi.org/10.1111/1756-185X.13433.
The authors thank the French Autoimmune and Autoinflammatory Rare Diseases Network (FAI2R).
French National Scleroderma Cohort coauthors: Zahir Amoura, Paris; Olivier Aumaitre, Clermont-Ferrand; Eric Auxenfants, Roubaix; Marie-Hélène Balquet, Lens; Cristina Belizna, Angers; Boris Bienvenu, Marseille; Emmanuel Chatelus, Strasbourg; Robin Dhote, Avicenne; Yves Frances, Marseille; Jean-Baptiste Gaultier, St-Etienne; Bernard Imbert, Grenoble; Jean-Emmanuel Kahn, Paris; Gilles Kaplanski, Marseille; Pierre Kieffer, Mulhouse; Noémie Le Gouellec, Valenciennes; Philippe Guilpain, Montpellier; Olivier Lidove, Paris; Nadine Magy-Bertrand, Besançon; François Maurier, Metz; Thomas Papo, Paris; Jean-Loup Pennaforte, Reims; Jacques Pouchot, Paris; Vivianne Queyrel, Nice; Denis Wahl, Nancy.
Ethics approval and consent to participate
In accordance with the French legislation, the database has received ethical approval from CCTIRS (approval no. 13.145; Advisory Committee on Information Processing in Material Research in the Field of Health). Data protection complied with the requirements of the National Information science and Liberties Commission and was recorded under no. 914607. Patients provided informed consent before entering the database.
Consent for publication
David Launay reports grants and personal fees from Actelion, Shire, CSL Behring, GSK, and Roche; personal fees from Boehringer Ingelheim; and grants from Octapharma outside the submitted work. Christian Agard reports grants and personal fees from Roche-Chugai and Boehringer Ingelheim outside the submitted work. Yannick Allanore reports grants from Inventiva and Sanofi; personal fees from Bayer, Boehringer Ingelheim, Chemomab, Inventiva, Roche, and Sanofi. Nadine Magy-Bertrand reports grants from GSK outside the submitted work. Olivier Lidove reports grants and personal fees from Amicus, Genzyme/Sanofi, and Shire HGT outside the submitted work. Francois Maurier reports personal fees from Actelion, and Chugaï; non-financial support from Amgen, Janssen, GSK, and Novartis outside the submitted work. Eric Hachulla reports grants from CSL Behring, Grifols and Octapharma; personal fees from Actelion, Bayer, Boehringer Ingelheim, GSK, Roche-Chugai, and Sanofi-Genzyme; grants from Actelion, GSK, and Roche-Chugai outside of the submitted work. Vincent Sobanski reports personal fees from Grifols; grants from Actelion, CSL Behring, Grifols, GSK, Octapharma, Pfizer, and Shire outside the submitted work. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Demographics and disease characteristics of patients with less than 2 years of the first non-RP symptom (n = 611)
Number of mRSS available in patients included in LCMM.
The 6 different LCMM tested
Model fit evaluation information for each LCMM tested
Averages of posterior probabilities of belonging to a class in each LCMM tested
Averages of posterior probabilities of belonging to a class in the 5-class LCMM
Sensitivity analysis: the 6 different LCMM with disease duration as adjustment factor
Sensitivity analysis: model fit evaluation information for each LCMM with disease duration as adjustment factor
Sensitivity analysis: averages of posterior probabilities of belonging to a class in each LCMM with disease duration as adjustment factor
Pulmonary function test values in patients with interstitial lung disease in the 5-class model (n = 72)
Survival analyses using Cox regression analysis without adjustment for age and sex in the 5-class LCMM
Survival analyses using Cox regression analysis adjusted for age and sex in the 5-class LCMM
About this article
Cite this article
Ledoult, E., Launay, D., Béhal, H. et al. Early trajectories of skin thickening are associated with severity and mortality in systemic sclerosis. Arthritis Res Ther 22, 30 (2020). https://doi.org/10.1186/s13075-020-2113-6
- Systemic sclerosis
- Modified Rodnan skin score
- Skin thickening trajectories
- Clinical heterogeneity