Abstract

Objectives

To characterize patients with positive anti-topoisomerase I (ATA) in lcSSc.

Methods

SSc patients enrolled in the EUSTAR cohort with a disease duration of ≤3 years at database entry were considered. We assessed the risk of major organ involvement in the following groups: ATA-lcSSc vs ACA-lcSSc and vs ANA without specificity (ANA)-lcSSc, and ATA-lcSSc vs ATA-dcSSc. Cox regression models with time-dependent covariates were performed with the following outcomes: new-onset interstitial lung disease (ILD), ILD progression [forced vital capacity (FVC) decline ≥10% and ≥5% vs values at ILD diagnosis), primary myocardial involvement (PMI), pulmonary hypertension (PH), any organ involvement and all-cause mortality.

Results

We included 1252 patients [194 ATA-lcSSc (15.5%)], with 7.7 years (s.d. 3.5) of follow-up. ILD risk was higher in ATA-lcSSc vs ACA- and ANA-lcSSc and similar to ATA-dcSSc, although with less frequent restrictive lung disease. The risk of FVC decline ≥10% (35% of ATA-lcSSc) was lower in ATA-lcSSc than in ATA-dcSSc, whereas FVC decline ≥5% occurs similarly between ATA-lcSSc (58% of patients) and other SSc subsets, including ATA-dcSSc. The risk of PMI was similar in ATA-lcSSc and ANA-lcSSc but lower than in ACA-lcSSc; no difference in PH and mortality risk was observed among lcSSc subsets. The risk of any organ involvement, PMI and PH was lower and the mortality tended to be lower in ATA-lcSSc vs ATA-dcSSc.

Conclusion

ATA-lcSSc patients have a high risk of ILD, albeit with a lower risk of progression compared with ATA-dcSSc, supporting careful screening for ILD in this subgroup.

Rheumatology key messages
  • ATA-lcSSc patients have a high risk of ILD, both at baseline and during follow-up.

  • ATA-lcSSc patients constitute a frequent SSc subset in which careful screening for ILD is recommended.

  • ILD progression is less frequent in ATA-lcSSc than ATA-dcSSc with an FVC decline ≥10%.

Introduction

SSc is a connective tissue disease characterized by an extremely wide heterogeneity, ranging from mild and stable to life-threatening forms [1, 2]. Besides the European Scleroderma Trial and Research (EUSTAR) activity score [3], reliable markers of disease activity and robust predictive factors are lacking in SSc. Therefore, improving risk stratification is of the utmost importance [4]. Patients affected with SSc are classified into limited (lcSSc) and diffuse (dcSSc) cutaneous forms [5]. However, recent evidence suggests that this classification may fail to fully capture the heterogeneity of the scleroderma spectrum [6]. Notably, there is renewed interest in the potential role of SSc-specific autoantibodies in predicting organ involvement [7] thanks to their high reproducibility and worldwide availability. Typically lcSSc patients have positive anticentromere (ACA), whereas those with dcSSc have positive anti-topoisomerase I (ATA). However, by exploring the ‘many faces of scleroderma’ back in 1988, Steen etal. [8] highlighted a consistent population of SSc patients with positive ATA but stable lcSSc, hypothesizing that it may be an ‘aborted form of dcSSc’ [9]. Although some cohort studies have more recently reported that this population may represent 11–23% of the lcSSc forms [10–12], very few have investigated the clinical spectrum of ATA-lcSSc patients, thus their disease presentation and course remain largely unknown [13, 14].

Interstitial lung disease (ILD) is the leading cause of death in SSc to date [15, 16]. The disease course in SSc-ILD is extremely variable, ranging from mild asymptomatic to rapidly progressive [17]. Interestingly, there have recently been some major achievements in recent randomized controlled trials (RCTs) leading to the labelling of two drugs for SSc-ILD: nintedanib and tocilizumab [18–21]. It must be pointed out that in the SENSCIS trial (NCT02597933), not only patients with dcSSc were enrolled, but also patients with lcSSc, with similar efficacy of nintedanib irrespective of the cutaneous subset, highlighting that some lcSSc may have significant ILD [22]. Risk stratification for SSc-ILD development and progression is a milestone towards early diagnosis and tailored treatments to improve outcome [23–25]. Our aim was to investigate the clinical features and risks of major organ involvement, including ILD, in patients with lcSSc and positive ATA in a large international multicentre cohort.

Methods

The study was conducted on prospectively collected data in the large multicentre EUSTAR database through July 2019. The structure of the online database, the collected data set and the definitions of clinical variables have been previously reported [26, 27]. Inclusion criteria were patients fulfilling the ACR 1980 and/or ACR/EULAR 2013 classification criteria for SSc [28], disease duration (from the first non-Raynaud’s sign/symptom) ≤3 years at database entry and known and stable cutaneous form during the first 3 years (i.e. lcSSc → lcSSc). In order to avoid the potential misclassification of some dcSSc patients as lcSSc in a very early stage of SSc, only patients with an available follow-up ≥3 years after database entry were included. Exclusion criteria were patients presenting multiple autoantibody reactivities and patients with unknown or positive anti-RNA polymerase III and anti-ENA antibodies other than ACA and ATA (i.e. anti-U1RNP and anti-PM/Scl). This was to avoid the inclusion of patients with multiple autoantibody positivity and those not tested for less-frequent but relevant SSc-associated antibodies (e.g. anti-RNA polymerase III).

Demographic and clinical data were recorded at baseline (i.e. database entry). Multivariable Cox regression models with time-dependent covariates were performed to assess the risk of developing the following outcomes during follow-up:

  • ILD: evidence of lung fibrosis on standard X-ray and/or high-resolution CT (HRCT)

  • Scleroderma renal crisis (SRC): according to the EUSTAR definition [29]

  • Primary myocardial involvement (PMI): ventricular arrhythmia or conduction defects by electrocardiogram or systolic dysfunction (i.e. ejection fraction <52% for males, <54% for females) [30] or diastolic dysfunction or significant pericardial effusion (>10 mm and diffuse) by transthoracic echocardiography (TTE)

  • Pulmonary hypertension (PH): systolic pulmonary arterial pressure (PAP) >45 mmHg by TTE and/or a mean PAP ≥25 mmHg on right heart catheterization (RHC)

  • Any major organ involvement: ILD and/or PMI and/or PH and/or SRC

  • Mortality: death from all-cause mortality.

Outside of Cox regression models, gastrointestinal (GI) involvement was defined as gastroesophageal reflux disease or stomach symptoms or malabsorption or intestinal symptoms or paralytic ileus.

As secondary endpoints, two adjunctive models were performed to assess the risk of ILD progression, defined by two different cut-offs of Forced vital capacity (FVC) decline:

  • ILD progression (FVC ≥10%): FVC decline ≥10% from ILD onset (defined as the first imaging evidence of ILD) to the last follow-up.

  • ILD progression (FVC ≥5 %): FVC decline ≥5% from ILD onset (defined as the first imaging evidence of ILD) to the last follow-up.

Restrictive lung disease was defined as FVC <80% of predicted value [31].

For each model, if FVC at ILD onset was not available, the closer FVC measure to ILD onset (within 1 year) was considered as the baseline FVC. The period of observation after ILD onset (i.e. ILD follow-up) was calculated from the baseline FVC to the first occurrence of progression or to the last available follow-up in patients without progression. Each Cox regression model was adjusted for age, sex and confounders.

Our SSc population was subdivided into four subsets: ATA-lcSSc, ACA-lcSSc, lcSSc with ANA positive but without anti-ENA specificity (ANA-lcSSc) and ATA-dcSSc. For baseline comparisons and Cox regression models, the analysis was performed in two groups of patients:

  • Group 1: lcSSc, by comparing patients with ATA-lcSSc vs ACA-lcSSc and vs ANA-lcSSc

  • Group 2: ATA-positive SSc, by comparing ATA-lcSSc vs ATA-dcSSc patients.

Ethics

Approval of the local ethics committee was obtained by each participating EUSTAR centre and informed consent was given by each registered patient.

Statistical analysis

For continuous variables, mean (s.d.) or median and interquartile range (IQR) are shown; counts and percentages are used for categorical variables. Group comparisons of baseline variables were performed using the Mann–Whitney test for continuous variables or the χ2 test for categorical variables. For groupwise comparisons of ACA-lcSSc and ANA-lcSSc vs ATA-lcSSc, P-values were adjusted for two parallel tests. Crude group comparisons of time until respective events of interest were conducted using Kaplan–Meier analysis. For multivariable time-to-event analysis, missing values were imputed 10 times with predictive mean matching and fixed random seed (packages mice and mitools); only cases with any available follow-up information were analysed. Multivariable Cox regression models allowing for time-dependent variables were used to investigate the influence of sets of variables, preselected according to medical considerations [32], on the events of interest. The starting point was the baseline visit for models assessing the risk of organ involvement and death and the first evidence of ILD for models evaluating ILD progression. The risk of event was expressed as the hazard ratio (HR) and its 95% confidence interval (CI). SPSS Statistics version 24 (IBM, Armonk, NY, USA) and R version 4.0.3 (R Foundation for Statistical Computing, Vienna, Austria) were used for analysis.

Results

Study population

Among the 16 831 total adult SSc patients, 4817 had a disease duration at database entry of ≤3 years; among them, 1949 were followed up for at least 3 years. In the latter group of patients, 133 (6.8%) transitioned from lcSSc to dcSSc (according to LeRoy et al. [5]) and were therefore excluded. Notably, the proportion of lcSSc patients treated with immunosuppressants was higher in those who transitioned to dcSSc (49.6%) compared with lcSSc who did not (32.3%) (P = 0.0003). After applying the other exclusion criteria (Supplementary Fig. S1, available at Rheumatology online), 1252 patients [1042 (83.2% female)] were enrolled; 95.5% of patients were Caucasian, 827 (66.1%) had lcSSc and 425 (33.9%) had dcSSc, the median age was 52 years (IQR 41.4–61.8) and disease duration at database entry was 1.3 years (IQR 0.7–2.2). Among patients with lcSSc, 194 (23.5%) had positive ATA, 537 (64.9%) had positive ACA and 96 (11.6%) had positive ANA without anti-ENA specificity. Among patients with positive ATA (n = 619), 194 (31.3%) had lcSSc and 425 (68.7%) had dcSSc (Tables 1 and 2).

Table 1

Demographic and clinical characteristics of patients with the limited cutaneous form at baseline

CharacteristicsPatients with data available, nTotal lcSSc (N = 827)ACA-lcSSc (n = 537)P-valueATA-lcSSc (n = 194)P-valueANA-lcSSc (n = 96)
Female, n (%)827740 (89.5)498 (92.7)<0.0001161 (83)0.76681 (84.4)
Age at RP onset, years82045.6 (15.2)46.9 (15.3)0.28444.2 (15.2)0.42346.3 (14.5)
Time from RP to SSc, years, median (IQR)8201.9 (0.1–7.9)2.5 (0.3–9.9)<0.00010.7 (0–3)0.9780.8 (0–4)
Age at SSc onset, years82751.4 (13.7)53 (13.3)<0.000146.9 (14.2)0.04351.1 (13.3)
Age, years82752.8 (13.7)54.4 (13.3)<0.000148.4 (14.3)0.0552.4 (13.3)
Disease duration, years8271.4 (0.9)1.4 (0.9)1.0001.4 (0.9)1.0001.4 (0.8)
BMI45525.4 (4.9)25.6 (5.1)1.00025.5 (4.7)1.00025.1 (4.4)
ESR >20 mm/ha, n (%)396113 (28.5)65 (26)0.0836 (35.3)0.34412 (27.3)
CRP elevation, n (%)38158 (13.2)25 (8.8)<0.000125 (23.6)0.6658 (16.7)
Hba, g/dl6413.2 (1.3)13.2 (1.3)1.00013.3 (1.3)0.92813.7 (1.0)
mRSS7524.8 (4.2)4.3 (4)0.0015.8 (4.8)1.0005.1 (3.9)
Puffy fingers (ever), n (%)451207 (45.9)123 (42.4)0.02562 (56.4)0.23622 (43.1)
DUs ever, n (%)820219 (26.7)141 (26.5)0.89359 (30.7)1.00019 (20.0)
DUs current, n (%)820200 (24.4)130 (24.4)0.75953 (27.6)0.31517 (17.9)
Telangiectasiaa, n (%)13867 (48.6)49 (55.1)0.03010 (30.3)0.3608 (50)
Joint synovitis, n (%)82197 (11.8)46 (8.6)0.00233 (17.2)1.00018 (18.8)
TFRs, n (%)81925 (3.1)9 (1.7)0.00811 (5.7)1.0006 (6.3)
Muscle weakness, n (%)821110 (13.4)52 (9.8)<0.000141 (21.1)0.98517 (17.7)
CK elevation, n (%)78355 (7.0)18 (3.6)<0.000125 (13.3)1.00012 (13.3)
GI involvement, n (%)826492 (59.6)338 (63.1)0.015101 (52.1)1.00053 (55.2)
GERD, n (%)824463 (56.2)323 (60.5)0.00593 (47.9)1.00047 (49.0)
Stomach symptoms, n (%)822123 (15)74 (13.9)0.92231 (16.1)1.00018 (18.8)
Intestinal symptoms, n (%)824145 (17.6)94 (17.6)1.00033 (17)1.00018 (18.8)
SRC, n (%)8237 (0.9)2 (0.4)0.5812 (1)0.3953 (3.1)
Serum creatininea, mg/dl, median (IQR)650.7 (0.7–0.9)0.8 (0.7–0.9)1.0000.7 (0.7–0.8)0.0540.65 (0.55–0.7)
Proteinuria, n (%)78122 (2.8)8 (1.6)0.01110 (5.4)1.0004 (4.4)
ILD, n (%)784163 (20.8)52 (10.3)<0.000187 (46.3)0.00424 (27.0)
Restrictive LD, n (%)78397 (12.4)43 (8.4)<0.000139 (21.4)0.62815 (16.3)
FVC, % predicted, median (IQR)417103 (91–115)106 (95–118)<0.000195 (84.5–106)0.357101 (85–114)
TLCa, % predicted, median (IQR)28799 (88–110)100 (90–111)0.01293 (81–107)0.79997 (82–110)
DLCO, % predicted, median (IQR)65677 (64–89)78 (66–90)0.00672 (59–86)1.00071 (58–90)
FEV1a, % predicted, median (IQR)28999 (86–109)100 (90–112)0.00694 (83–103)0.56598.5 (87–105)
PHa, n (%)32114 (4.4)11 (5.1)0.3501 (1.4)0.4392 (5.6)
sPAPa, mmHg31128.4 (10.7)28.5 (11.2)1.00027.7 (6.7)1.00029.3 (13.2)
PMI, n (%)827154 (18.1)100 (18.6)1.00036 (18.6)1.00018 (18.8)
PMI (<45%), n (%)812147 (17.9)95(18.1)1.00034(17.8)1.00018 (18.8)
Ejection fractiona, %39162.7 (5.8)62.7 (5.7)1.00062.5 (6.7)1.00063.5 (4.5)
Cardiac blocks, n (%)76952 (6.8)26 (5.2)0.10417 (9.3)1.0009 (10.3)
Diastolic dysfunction, n (%)773109 (14.1)77 (15.4)0.63222 (12.3)1.00010 (10.8)
Pericardial effusiona, n (%)39815 (3.8)9 (3.5)1.0004 (4.1)0.4834 (8.7)
Any severe, n (%)827293 (35.4)149 (27.7)<0.0001104 (53.6)0.11140 (41.7)
EScSG activity index 2001, median (IQR)8270.5 (0–1)0.5 (0–0.5)<0.00010.5 (0–1.5)0.1260.5 (0–1)
EScSG activity index 2016a, median (IQR)642.7 (1.3–3.7)2.4 (0.7–3.4)0.2633.5 (1.5–5)0.2602.4 (1.8–2.6)
EScSG 2001 active, n (%)79631 (3.7)13 (2.4)0.00216 (8.2)0.2142 (2.1)
EScSG 2016 activea, n (%)6435 (54.7)18 (48.6)0.18014 (73.7)0.1893 (37.5)
Corticosteroids, n (%)827174 (21)73 (13.6)<0.000163 (32.5)0.46438 (39.6)
Prednisone equivalent, mg/day, median (IQR)8275 (3.5–5)5 (2.5–5)0.1335 (4–6)0.5675 (3.4–5)
Immunosuppressants, n (%)827267 (32.3)124 (23.1)<0.000190 (46.4)0.31553 (55.2)
CharacteristicsPatients with data available, nTotal lcSSc (N = 827)ACA-lcSSc (n = 537)P-valueATA-lcSSc (n = 194)P-valueANA-lcSSc (n = 96)
Female, n (%)827740 (89.5)498 (92.7)<0.0001161 (83)0.76681 (84.4)
Age at RP onset, years82045.6 (15.2)46.9 (15.3)0.28444.2 (15.2)0.42346.3 (14.5)
Time from RP to SSc, years, median (IQR)8201.9 (0.1–7.9)2.5 (0.3–9.9)<0.00010.7 (0–3)0.9780.8 (0–4)
Age at SSc onset, years82751.4 (13.7)53 (13.3)<0.000146.9 (14.2)0.04351.1 (13.3)
Age, years82752.8 (13.7)54.4 (13.3)<0.000148.4 (14.3)0.0552.4 (13.3)
Disease duration, years8271.4 (0.9)1.4 (0.9)1.0001.4 (0.9)1.0001.4 (0.8)
BMI45525.4 (4.9)25.6 (5.1)1.00025.5 (4.7)1.00025.1 (4.4)
ESR >20 mm/ha, n (%)396113 (28.5)65 (26)0.0836 (35.3)0.34412 (27.3)
CRP elevation, n (%)38158 (13.2)25 (8.8)<0.000125 (23.6)0.6658 (16.7)
Hba, g/dl6413.2 (1.3)13.2 (1.3)1.00013.3 (1.3)0.92813.7 (1.0)
mRSS7524.8 (4.2)4.3 (4)0.0015.8 (4.8)1.0005.1 (3.9)
Puffy fingers (ever), n (%)451207 (45.9)123 (42.4)0.02562 (56.4)0.23622 (43.1)
DUs ever, n (%)820219 (26.7)141 (26.5)0.89359 (30.7)1.00019 (20.0)
DUs current, n (%)820200 (24.4)130 (24.4)0.75953 (27.6)0.31517 (17.9)
Telangiectasiaa, n (%)13867 (48.6)49 (55.1)0.03010 (30.3)0.3608 (50)
Joint synovitis, n (%)82197 (11.8)46 (8.6)0.00233 (17.2)1.00018 (18.8)
TFRs, n (%)81925 (3.1)9 (1.7)0.00811 (5.7)1.0006 (6.3)
Muscle weakness, n (%)821110 (13.4)52 (9.8)<0.000141 (21.1)0.98517 (17.7)
CK elevation, n (%)78355 (7.0)18 (3.6)<0.000125 (13.3)1.00012 (13.3)
GI involvement, n (%)826492 (59.6)338 (63.1)0.015101 (52.1)1.00053 (55.2)
GERD, n (%)824463 (56.2)323 (60.5)0.00593 (47.9)1.00047 (49.0)
Stomach symptoms, n (%)822123 (15)74 (13.9)0.92231 (16.1)1.00018 (18.8)
Intestinal symptoms, n (%)824145 (17.6)94 (17.6)1.00033 (17)1.00018 (18.8)
SRC, n (%)8237 (0.9)2 (0.4)0.5812 (1)0.3953 (3.1)
Serum creatininea, mg/dl, median (IQR)650.7 (0.7–0.9)0.8 (0.7–0.9)1.0000.7 (0.7–0.8)0.0540.65 (0.55–0.7)
Proteinuria, n (%)78122 (2.8)8 (1.6)0.01110 (5.4)1.0004 (4.4)
ILD, n (%)784163 (20.8)52 (10.3)<0.000187 (46.3)0.00424 (27.0)
Restrictive LD, n (%)78397 (12.4)43 (8.4)<0.000139 (21.4)0.62815 (16.3)
FVC, % predicted, median (IQR)417103 (91–115)106 (95–118)<0.000195 (84.5–106)0.357101 (85–114)
TLCa, % predicted, median (IQR)28799 (88–110)100 (90–111)0.01293 (81–107)0.79997 (82–110)
DLCO, % predicted, median (IQR)65677 (64–89)78 (66–90)0.00672 (59–86)1.00071 (58–90)
FEV1a, % predicted, median (IQR)28999 (86–109)100 (90–112)0.00694 (83–103)0.56598.5 (87–105)
PHa, n (%)32114 (4.4)11 (5.1)0.3501 (1.4)0.4392 (5.6)
sPAPa, mmHg31128.4 (10.7)28.5 (11.2)1.00027.7 (6.7)1.00029.3 (13.2)
PMI, n (%)827154 (18.1)100 (18.6)1.00036 (18.6)1.00018 (18.8)
PMI (<45%), n (%)812147 (17.9)95(18.1)1.00034(17.8)1.00018 (18.8)
Ejection fractiona, %39162.7 (5.8)62.7 (5.7)1.00062.5 (6.7)1.00063.5 (4.5)
Cardiac blocks, n (%)76952 (6.8)26 (5.2)0.10417 (9.3)1.0009 (10.3)
Diastolic dysfunction, n (%)773109 (14.1)77 (15.4)0.63222 (12.3)1.00010 (10.8)
Pericardial effusiona, n (%)39815 (3.8)9 (3.5)1.0004 (4.1)0.4834 (8.7)
Any severe, n (%)827293 (35.4)149 (27.7)<0.0001104 (53.6)0.11140 (41.7)
EScSG activity index 2001, median (IQR)8270.5 (0–1)0.5 (0–0.5)<0.00010.5 (0–1.5)0.1260.5 (0–1)
EScSG activity index 2016a, median (IQR)642.7 (1.3–3.7)2.4 (0.7–3.4)0.2633.5 (1.5–5)0.2602.4 (1.8–2.6)
EScSG 2001 active, n (%)79631 (3.7)13 (2.4)0.00216 (8.2)0.2142 (2.1)
EScSG 2016 activea, n (%)6435 (54.7)18 (48.6)0.18014 (73.7)0.1893 (37.5)
Corticosteroids, n (%)827174 (21)73 (13.6)<0.000163 (32.5)0.46438 (39.6)
Prednisone equivalent, mg/day, median (IQR)8275 (3.5–5)5 (2.5–5)0.1335 (4–6)0.5675 (3.4–5)
Immunosuppressants, n (%)827267 (32.3)124 (23.1)<0.000190 (46.4)0.31553 (55.2)

Values are presented as mean (s.d.) unless stated otherwise.

a

Data available in less than half of patients.

Significant P-values (<0.05) are in bold.

CK: creatine kinase; DLCO: diffusing capacity for carbon monoxide; DUs: digital ulcers; EScSG: European Scleroderma Study Group; FEV1: forced expiratory volume in the first second; FVC: forced vital capacity; GI: gastrointestinal; GERD: gastroesophageal reflux disease; Hb: haemoglobin; ILD: interstitial lung disease; LD: lung disease; sPAP: systolic pulmonary artery pressure; immunosuppressants: CYC, SSZ, HCQ, MTX, LEF, AZA, MMF, ciclosporin A, d-penicillamine, rituximab, imatinib, anti-TNF inhibitors, tocilizumab or abatacept; mRSS: modified Rodnan Skin Score; PH: pulmonary hypertension; PMI: primary myocardial involvement; SRC: scleroderma renal crisis; TFRs: tendon friction rubs; TLC: total lung capacity.

Table 1

Demographic and clinical characteristics of patients with the limited cutaneous form at baseline

CharacteristicsPatients with data available, nTotal lcSSc (N = 827)ACA-lcSSc (n = 537)P-valueATA-lcSSc (n = 194)P-valueANA-lcSSc (n = 96)
Female, n (%)827740 (89.5)498 (92.7)<0.0001161 (83)0.76681 (84.4)
Age at RP onset, years82045.6 (15.2)46.9 (15.3)0.28444.2 (15.2)0.42346.3 (14.5)
Time from RP to SSc, years, median (IQR)8201.9 (0.1–7.9)2.5 (0.3–9.9)<0.00010.7 (0–3)0.9780.8 (0–4)
Age at SSc onset, years82751.4 (13.7)53 (13.3)<0.000146.9 (14.2)0.04351.1 (13.3)
Age, years82752.8 (13.7)54.4 (13.3)<0.000148.4 (14.3)0.0552.4 (13.3)
Disease duration, years8271.4 (0.9)1.4 (0.9)1.0001.4 (0.9)1.0001.4 (0.8)
BMI45525.4 (4.9)25.6 (5.1)1.00025.5 (4.7)1.00025.1 (4.4)
ESR >20 mm/ha, n (%)396113 (28.5)65 (26)0.0836 (35.3)0.34412 (27.3)
CRP elevation, n (%)38158 (13.2)25 (8.8)<0.000125 (23.6)0.6658 (16.7)
Hba, g/dl6413.2 (1.3)13.2 (1.3)1.00013.3 (1.3)0.92813.7 (1.0)
mRSS7524.8 (4.2)4.3 (4)0.0015.8 (4.8)1.0005.1 (3.9)
Puffy fingers (ever), n (%)451207 (45.9)123 (42.4)0.02562 (56.4)0.23622 (43.1)
DUs ever, n (%)820219 (26.7)141 (26.5)0.89359 (30.7)1.00019 (20.0)
DUs current, n (%)820200 (24.4)130 (24.4)0.75953 (27.6)0.31517 (17.9)
Telangiectasiaa, n (%)13867 (48.6)49 (55.1)0.03010 (30.3)0.3608 (50)
Joint synovitis, n (%)82197 (11.8)46 (8.6)0.00233 (17.2)1.00018 (18.8)
TFRs, n (%)81925 (3.1)9 (1.7)0.00811 (5.7)1.0006 (6.3)
Muscle weakness, n (%)821110 (13.4)52 (9.8)<0.000141 (21.1)0.98517 (17.7)
CK elevation, n (%)78355 (7.0)18 (3.6)<0.000125 (13.3)1.00012 (13.3)
GI involvement, n (%)826492 (59.6)338 (63.1)0.015101 (52.1)1.00053 (55.2)
GERD, n (%)824463 (56.2)323 (60.5)0.00593 (47.9)1.00047 (49.0)
Stomach symptoms, n (%)822123 (15)74 (13.9)0.92231 (16.1)1.00018 (18.8)
Intestinal symptoms, n (%)824145 (17.6)94 (17.6)1.00033 (17)1.00018 (18.8)
SRC, n (%)8237 (0.9)2 (0.4)0.5812 (1)0.3953 (3.1)
Serum creatininea, mg/dl, median (IQR)650.7 (0.7–0.9)0.8 (0.7–0.9)1.0000.7 (0.7–0.8)0.0540.65 (0.55–0.7)
Proteinuria, n (%)78122 (2.8)8 (1.6)0.01110 (5.4)1.0004 (4.4)
ILD, n (%)784163 (20.8)52 (10.3)<0.000187 (46.3)0.00424 (27.0)
Restrictive LD, n (%)78397 (12.4)43 (8.4)<0.000139 (21.4)0.62815 (16.3)
FVC, % predicted, median (IQR)417103 (91–115)106 (95–118)<0.000195 (84.5–106)0.357101 (85–114)
TLCa, % predicted, median (IQR)28799 (88–110)100 (90–111)0.01293 (81–107)0.79997 (82–110)
DLCO, % predicted, median (IQR)65677 (64–89)78 (66–90)0.00672 (59–86)1.00071 (58–90)
FEV1a, % predicted, median (IQR)28999 (86–109)100 (90–112)0.00694 (83–103)0.56598.5 (87–105)
PHa, n (%)32114 (4.4)11 (5.1)0.3501 (1.4)0.4392 (5.6)
sPAPa, mmHg31128.4 (10.7)28.5 (11.2)1.00027.7 (6.7)1.00029.3 (13.2)
PMI, n (%)827154 (18.1)100 (18.6)1.00036 (18.6)1.00018 (18.8)
PMI (<45%), n (%)812147 (17.9)95(18.1)1.00034(17.8)1.00018 (18.8)
Ejection fractiona, %39162.7 (5.8)62.7 (5.7)1.00062.5 (6.7)1.00063.5 (4.5)
Cardiac blocks, n (%)76952 (6.8)26 (5.2)0.10417 (9.3)1.0009 (10.3)
Diastolic dysfunction, n (%)773109 (14.1)77 (15.4)0.63222 (12.3)1.00010 (10.8)
Pericardial effusiona, n (%)39815 (3.8)9 (3.5)1.0004 (4.1)0.4834 (8.7)
Any severe, n (%)827293 (35.4)149 (27.7)<0.0001104 (53.6)0.11140 (41.7)
EScSG activity index 2001, median (IQR)8270.5 (0–1)0.5 (0–0.5)<0.00010.5 (0–1.5)0.1260.5 (0–1)
EScSG activity index 2016a, median (IQR)642.7 (1.3–3.7)2.4 (0.7–3.4)0.2633.5 (1.5–5)0.2602.4 (1.8–2.6)
EScSG 2001 active, n (%)79631 (3.7)13 (2.4)0.00216 (8.2)0.2142 (2.1)
EScSG 2016 activea, n (%)6435 (54.7)18 (48.6)0.18014 (73.7)0.1893 (37.5)
Corticosteroids, n (%)827174 (21)73 (13.6)<0.000163 (32.5)0.46438 (39.6)
Prednisone equivalent, mg/day, median (IQR)8275 (3.5–5)5 (2.5–5)0.1335 (4–6)0.5675 (3.4–5)
Immunosuppressants, n (%)827267 (32.3)124 (23.1)<0.000190 (46.4)0.31553 (55.2)
CharacteristicsPatients with data available, nTotal lcSSc (N = 827)ACA-lcSSc (n = 537)P-valueATA-lcSSc (n = 194)P-valueANA-lcSSc (n = 96)
Female, n (%)827740 (89.5)498 (92.7)<0.0001161 (83)0.76681 (84.4)
Age at RP onset, years82045.6 (15.2)46.9 (15.3)0.28444.2 (15.2)0.42346.3 (14.5)
Time from RP to SSc, years, median (IQR)8201.9 (0.1–7.9)2.5 (0.3–9.9)<0.00010.7 (0–3)0.9780.8 (0–4)
Age at SSc onset, years82751.4 (13.7)53 (13.3)<0.000146.9 (14.2)0.04351.1 (13.3)
Age, years82752.8 (13.7)54.4 (13.3)<0.000148.4 (14.3)0.0552.4 (13.3)
Disease duration, years8271.4 (0.9)1.4 (0.9)1.0001.4 (0.9)1.0001.4 (0.8)
BMI45525.4 (4.9)25.6 (5.1)1.00025.5 (4.7)1.00025.1 (4.4)
ESR >20 mm/ha, n (%)396113 (28.5)65 (26)0.0836 (35.3)0.34412 (27.3)
CRP elevation, n (%)38158 (13.2)25 (8.8)<0.000125 (23.6)0.6658 (16.7)
Hba, g/dl6413.2 (1.3)13.2 (1.3)1.00013.3 (1.3)0.92813.7 (1.0)
mRSS7524.8 (4.2)4.3 (4)0.0015.8 (4.8)1.0005.1 (3.9)
Puffy fingers (ever), n (%)451207 (45.9)123 (42.4)0.02562 (56.4)0.23622 (43.1)
DUs ever, n (%)820219 (26.7)141 (26.5)0.89359 (30.7)1.00019 (20.0)
DUs current, n (%)820200 (24.4)130 (24.4)0.75953 (27.6)0.31517 (17.9)
Telangiectasiaa, n (%)13867 (48.6)49 (55.1)0.03010 (30.3)0.3608 (50)
Joint synovitis, n (%)82197 (11.8)46 (8.6)0.00233 (17.2)1.00018 (18.8)
TFRs, n (%)81925 (3.1)9 (1.7)0.00811 (5.7)1.0006 (6.3)
Muscle weakness, n (%)821110 (13.4)52 (9.8)<0.000141 (21.1)0.98517 (17.7)
CK elevation, n (%)78355 (7.0)18 (3.6)<0.000125 (13.3)1.00012 (13.3)
GI involvement, n (%)826492 (59.6)338 (63.1)0.015101 (52.1)1.00053 (55.2)
GERD, n (%)824463 (56.2)323 (60.5)0.00593 (47.9)1.00047 (49.0)
Stomach symptoms, n (%)822123 (15)74 (13.9)0.92231 (16.1)1.00018 (18.8)
Intestinal symptoms, n (%)824145 (17.6)94 (17.6)1.00033 (17)1.00018 (18.8)
SRC, n (%)8237 (0.9)2 (0.4)0.5812 (1)0.3953 (3.1)
Serum creatininea, mg/dl, median (IQR)650.7 (0.7–0.9)0.8 (0.7–0.9)1.0000.7 (0.7–0.8)0.0540.65 (0.55–0.7)
Proteinuria, n (%)78122 (2.8)8 (1.6)0.01110 (5.4)1.0004 (4.4)
ILD, n (%)784163 (20.8)52 (10.3)<0.000187 (46.3)0.00424 (27.0)
Restrictive LD, n (%)78397 (12.4)43 (8.4)<0.000139 (21.4)0.62815 (16.3)
FVC, % predicted, median (IQR)417103 (91–115)106 (95–118)<0.000195 (84.5–106)0.357101 (85–114)
TLCa, % predicted, median (IQR)28799 (88–110)100 (90–111)0.01293 (81–107)0.79997 (82–110)
DLCO, % predicted, median (IQR)65677 (64–89)78 (66–90)0.00672 (59–86)1.00071 (58–90)
FEV1a, % predicted, median (IQR)28999 (86–109)100 (90–112)0.00694 (83–103)0.56598.5 (87–105)
PHa, n (%)32114 (4.4)11 (5.1)0.3501 (1.4)0.4392 (5.6)
sPAPa, mmHg31128.4 (10.7)28.5 (11.2)1.00027.7 (6.7)1.00029.3 (13.2)
PMI, n (%)827154 (18.1)100 (18.6)1.00036 (18.6)1.00018 (18.8)
PMI (<45%), n (%)812147 (17.9)95(18.1)1.00034(17.8)1.00018 (18.8)
Ejection fractiona, %39162.7 (5.8)62.7 (5.7)1.00062.5 (6.7)1.00063.5 (4.5)
Cardiac blocks, n (%)76952 (6.8)26 (5.2)0.10417 (9.3)1.0009 (10.3)
Diastolic dysfunction, n (%)773109 (14.1)77 (15.4)0.63222 (12.3)1.00010 (10.8)
Pericardial effusiona, n (%)39815 (3.8)9 (3.5)1.0004 (4.1)0.4834 (8.7)
Any severe, n (%)827293 (35.4)149 (27.7)<0.0001104 (53.6)0.11140 (41.7)
EScSG activity index 2001, median (IQR)8270.5 (0–1)0.5 (0–0.5)<0.00010.5 (0–1.5)0.1260.5 (0–1)
EScSG activity index 2016a, median (IQR)642.7 (1.3–3.7)2.4 (0.7–3.4)0.2633.5 (1.5–5)0.2602.4 (1.8–2.6)
EScSG 2001 active, n (%)79631 (3.7)13 (2.4)0.00216 (8.2)0.2142 (2.1)
EScSG 2016 activea, n (%)6435 (54.7)18 (48.6)0.18014 (73.7)0.1893 (37.5)
Corticosteroids, n (%)827174 (21)73 (13.6)<0.000163 (32.5)0.46438 (39.6)
Prednisone equivalent, mg/day, median (IQR)8275 (3.5–5)5 (2.5–5)0.1335 (4–6)0.5675 (3.4–5)
Immunosuppressants, n (%)827267 (32.3)124 (23.1)<0.000190 (46.4)0.31553 (55.2)

Values are presented as mean (s.d.) unless stated otherwise.

a

Data available in less than half of patients.

Significant P-values (<0.05) are in bold.

CK: creatine kinase; DLCO: diffusing capacity for carbon monoxide; DUs: digital ulcers; EScSG: European Scleroderma Study Group; FEV1: forced expiratory volume in the first second; FVC: forced vital capacity; GI: gastrointestinal; GERD: gastroesophageal reflux disease; Hb: haemoglobin; ILD: interstitial lung disease; LD: lung disease; sPAP: systolic pulmonary artery pressure; immunosuppressants: CYC, SSZ, HCQ, MTX, LEF, AZA, MMF, ciclosporin A, d-penicillamine, rituximab, imatinib, anti-TNF inhibitors, tocilizumab or abatacept; mRSS: modified Rodnan Skin Score; PH: pulmonary hypertension; PMI: primary myocardial involvement; SRC: scleroderma renal crisis; TFRs: tendon friction rubs; TLC: total lung capacity.

Table 2

Demographic and clinical characteristics of ATA-positive patients at baseline

CharacteristicsPatients with data available, nTotal ATA (N = 619)ATA-dcSSc (n = 425)P-valueATA-lcSSc (n = 194)
Female, n (%)619463 (74.8)302 (71.1)0.003161 (83)
Age at RP onset, years61545.3 (34.5–55.6)45.1 (35.4–55.3)0.83245.4 (31.4–56)
Time from RP to SSc, years6150.3 (0–2)0.1 (0–1.2)<0.00010.7 (0–3)
Age at SSc onset, years61947.2 (36.9–57)46.8 (36.9–56.3)1.00047.8 (37.3–57.8)
Age, years61948.9 (38.6–58.6)48.2 (38.6–58.2)1.00049.7 (38.4–59.3)
Disease duration, years, mean (s.d.)6191.4 (0.8)1.5 (0.8)1.0001.4 (0.9)
BMI, mean (s.d.)35724.7 (4.4)24.3 (4.2)0.03025.5 (4.7)
ESR >20 mm/h, n (%)330147 (44.5)111 (48.7)0.04736 (35.3)
CRP elevation, n (%)337109 (32.3)84 (36.4)0.04025 (23.6)
Hba, g/dl, mean (s.d.)5713.4 (0.2)13.4 (0.2)1.00013.3 (0.3)
mRSS58612.5 (5–20)16 (11–23)<0.00015 (2–8)
Puffy fingers (ever), n (%)336115 (34.2)79 (35)1.00036 (32.7)
DUs ever, n (%)616212 (34.4)153 (36)0.42659 (30.7)
DUs current, n (%)616190 (30.8)137 (32.3)0.48353 (27.6)
Teleangiectasiaa, n (%)10946 (42.2)36 (47.4)0.19510 (30.3)
Joint synovitis, n (%)615145 (23.6)112 (26.5)0.02433 (17.2)
TFRs, n (%)610112 (18.4)102 (24.4)<0.000111(5.7)
Muscle weakness, n (%)615145 (23.6)104 (24.7)0.66541 (21.1)
CK elevation, n (%)58083 (14.3)58 (14.8)1.00025 (13.3)
GI involvement, n (%)619365 (59)264 (62.1)0.037101 (52.1)
GERD, n (%)618344 (55.7)251 (59.2)0.01893 (47.9)
Stomach symptoms, n (%)617119 (19.3)88 (20.8)0.34231 (16.1)
Intestinal symptoms, n (%)619101 (16.3)68 (16.0)1.00033 (17.0)
SRC, n (%)61810 (1.6)8 (1.9)0.8682 (1.0)
Serum creatininea, mg/dl, mean (s.d.)580.83 (0.73)0.72 (0.16)0.2091.02 (1.19)
Proteinuria, n (%)58948 (8.1)38 (9.4)0.18910 (5.4)
ILD, n (%)598290 (48.5)203 (49.5)0.92587 (46.3)
Restrictive LD, n (%)579202 (34.9)163 (41.1)<0.000139 (21.4)
FVC, % predicted31790 (77–101)86 (75–99)<0.000195 (84.5–106)
TLCa, % predicted20190 (76–100)86 (73–99)0.00893 (81–107)
DLCO, % predicted48868 (54–82)64 (52–80)0.00172 (59–86)
FEV1a, % predicted23289 (78–100)87 (76–99)0.01594 (83–103)
PHa, n (%)2439 (3.7)8 (4.7)0.4471 (1.4)
sPAPa, mmHg, mean (s.d.)23328.7 (8.4)29.1 (8.9)0.38927.7 (6.7)
PMI, n (%)619134 (21.6)98 (23.1)0.41436 (18.6)
PMI (<45%), n (%)607127 (20.9)93 (22.4)0.40034 (17.8)
Ejection fractiona, %, mean (s.d.)29162.7 (6.2)62.8 (6)1.00062.5 (6.7)
Cardiac blocks, n (%)58157 (9.8)40 (10.1)1.00017 (9.3)
Diastolic dysfunction, n (%)57187 (15.2)65 (16.6)0.37122 (12.3)
Pericardial effusiona, n (%)30115 (5)11 (5.4)1.0004 (4.1)
Any severe, n (%)619351 (56.7)247 (58.1)0.587104 (53.6)
EScSG activity index 20016191 (0.5–2.5)1.5 (0.5–3)<0.00010.5 (0–1.5)
EScSG activity index 2016a524 (3.1–5.4)4.4 (3.8–5.6)0.0393.5 (1.5–5)
EScSG 2001 active, n (%)619103 (16.6)87 (20.5)<0.000116 (8.2)
EScSG 2016 activea, n (%)5244 (84.6)30 (90.9)0.09714 (73.7)
Corticosteroids, n (%)619228 (36.8)165 (38.8)0.25763 (32.5)
Prednisone equivalent, mg/day6195 (5–7.5)5 (5–7.5)0.0595 (4–6)
Immunosuppressants, n (%)619332 (53.6)242 (56.9)0.02990 (46.4)
CharacteristicsPatients with data available, nTotal ATA (N = 619)ATA-dcSSc (n = 425)P-valueATA-lcSSc (n = 194)
Female, n (%)619463 (74.8)302 (71.1)0.003161 (83)
Age at RP onset, years61545.3 (34.5–55.6)45.1 (35.4–55.3)0.83245.4 (31.4–56)
Time from RP to SSc, years6150.3 (0–2)0.1 (0–1.2)<0.00010.7 (0–3)
Age at SSc onset, years61947.2 (36.9–57)46.8 (36.9–56.3)1.00047.8 (37.3–57.8)
Age, years61948.9 (38.6–58.6)48.2 (38.6–58.2)1.00049.7 (38.4–59.3)
Disease duration, years, mean (s.d.)6191.4 (0.8)1.5 (0.8)1.0001.4 (0.9)
BMI, mean (s.d.)35724.7 (4.4)24.3 (4.2)0.03025.5 (4.7)
ESR >20 mm/h, n (%)330147 (44.5)111 (48.7)0.04736 (35.3)
CRP elevation, n (%)337109 (32.3)84 (36.4)0.04025 (23.6)
Hba, g/dl, mean (s.d.)5713.4 (0.2)13.4 (0.2)1.00013.3 (0.3)
mRSS58612.5 (5–20)16 (11–23)<0.00015 (2–8)
Puffy fingers (ever), n (%)336115 (34.2)79 (35)1.00036 (32.7)
DUs ever, n (%)616212 (34.4)153 (36)0.42659 (30.7)
DUs current, n (%)616190 (30.8)137 (32.3)0.48353 (27.6)
Teleangiectasiaa, n (%)10946 (42.2)36 (47.4)0.19510 (30.3)
Joint synovitis, n (%)615145 (23.6)112 (26.5)0.02433 (17.2)
TFRs, n (%)610112 (18.4)102 (24.4)<0.000111(5.7)
Muscle weakness, n (%)615145 (23.6)104 (24.7)0.66541 (21.1)
CK elevation, n (%)58083 (14.3)58 (14.8)1.00025 (13.3)
GI involvement, n (%)619365 (59)264 (62.1)0.037101 (52.1)
GERD, n (%)618344 (55.7)251 (59.2)0.01893 (47.9)
Stomach symptoms, n (%)617119 (19.3)88 (20.8)0.34231 (16.1)
Intestinal symptoms, n (%)619101 (16.3)68 (16.0)1.00033 (17.0)
SRC, n (%)61810 (1.6)8 (1.9)0.8682 (1.0)
Serum creatininea, mg/dl, mean (s.d.)580.83 (0.73)0.72 (0.16)0.2091.02 (1.19)
Proteinuria, n (%)58948 (8.1)38 (9.4)0.18910 (5.4)
ILD, n (%)598290 (48.5)203 (49.5)0.92587 (46.3)
Restrictive LD, n (%)579202 (34.9)163 (41.1)<0.000139 (21.4)
FVC, % predicted31790 (77–101)86 (75–99)<0.000195 (84.5–106)
TLCa, % predicted20190 (76–100)86 (73–99)0.00893 (81–107)
DLCO, % predicted48868 (54–82)64 (52–80)0.00172 (59–86)
FEV1a, % predicted23289 (78–100)87 (76–99)0.01594 (83–103)
PHa, n (%)2439 (3.7)8 (4.7)0.4471 (1.4)
sPAPa, mmHg, mean (s.d.)23328.7 (8.4)29.1 (8.9)0.38927.7 (6.7)
PMI, n (%)619134 (21.6)98 (23.1)0.41436 (18.6)
PMI (<45%), n (%)607127 (20.9)93 (22.4)0.40034 (17.8)
Ejection fractiona, %, mean (s.d.)29162.7 (6.2)62.8 (6)1.00062.5 (6.7)
Cardiac blocks, n (%)58157 (9.8)40 (10.1)1.00017 (9.3)
Diastolic dysfunction, n (%)57187 (15.2)65 (16.6)0.37122 (12.3)
Pericardial effusiona, n (%)30115 (5)11 (5.4)1.0004 (4.1)
Any severe, n (%)619351 (56.7)247 (58.1)0.587104 (53.6)
EScSG activity index 20016191 (0.5–2.5)1.5 (0.5–3)<0.00010.5 (0–1.5)
EScSG activity index 2016a524 (3.1–5.4)4.4 (3.8–5.6)0.0393.5 (1.5–5)
EScSG 2001 active, n (%)619103 (16.6)87 (20.5)<0.000116 (8.2)
EScSG 2016 activea, n (%)5244 (84.6)30 (90.9)0.09714 (73.7)
Corticosteroids, n (%)619228 (36.8)165 (38.8)0.25763 (32.5)
Prednisone equivalent, mg/day6195 (5–7.5)5 (5–7.5)0.0595 (4–6)
Immunosuppressants, n (%)619332 (53.6)242 (56.9)0.02990 (46.4)

Values are presented as median (IQR) unless stated otherwise.

a

Data available in less than half of patients.

Significant P-values (<0.05) are in bold.

CK: creatine kinase; DUs: digital ulcers; EScSG: European Scleroderma Study Group; FEV1: forced expiratory volume in the first second; GERD: gastroesophageal reflux disease; Hb: haemoglobin; LD: lung disease; sPAP: systolic pulmonary artery pressure; immunosuppressants: CYC, SSZ, HCQ, MTX, LEF, AZA, MMF, ciclosporin A, d-penicillamine, rituximab, imatinib, anti-TNF inhibitors, tocilizumab or abatacept; DLCO: diffusing capacity for carbon monoxide; FVC: forced vital capacity; GI: gastrointestinal; ILD: interstitial lung disease; mRSS: modified Rodnan Skin Score; PH: pulmonary hypertension; PMI: primary myocardial involvement; SRC: scleroderma renal crisis; TFRs: tendon friction rubs; TLC: total lung capacity. 

Table 2

Demographic and clinical characteristics of ATA-positive patients at baseline

CharacteristicsPatients with data available, nTotal ATA (N = 619)ATA-dcSSc (n = 425)P-valueATA-lcSSc (n = 194)
Female, n (%)619463 (74.8)302 (71.1)0.003161 (83)
Age at RP onset, years61545.3 (34.5–55.6)45.1 (35.4–55.3)0.83245.4 (31.4–56)
Time from RP to SSc, years6150.3 (0–2)0.1 (0–1.2)<0.00010.7 (0–3)
Age at SSc onset, years61947.2 (36.9–57)46.8 (36.9–56.3)1.00047.8 (37.3–57.8)
Age, years61948.9 (38.6–58.6)48.2 (38.6–58.2)1.00049.7 (38.4–59.3)
Disease duration, years, mean (s.d.)6191.4 (0.8)1.5 (0.8)1.0001.4 (0.9)
BMI, mean (s.d.)35724.7 (4.4)24.3 (4.2)0.03025.5 (4.7)
ESR >20 mm/h, n (%)330147 (44.5)111 (48.7)0.04736 (35.3)
CRP elevation, n (%)337109 (32.3)84 (36.4)0.04025 (23.6)
Hba, g/dl, mean (s.d.)5713.4 (0.2)13.4 (0.2)1.00013.3 (0.3)
mRSS58612.5 (5–20)16 (11–23)<0.00015 (2–8)
Puffy fingers (ever), n (%)336115 (34.2)79 (35)1.00036 (32.7)
DUs ever, n (%)616212 (34.4)153 (36)0.42659 (30.7)
DUs current, n (%)616190 (30.8)137 (32.3)0.48353 (27.6)
Teleangiectasiaa, n (%)10946 (42.2)36 (47.4)0.19510 (30.3)
Joint synovitis, n (%)615145 (23.6)112 (26.5)0.02433 (17.2)
TFRs, n (%)610112 (18.4)102 (24.4)<0.000111(5.7)
Muscle weakness, n (%)615145 (23.6)104 (24.7)0.66541 (21.1)
CK elevation, n (%)58083 (14.3)58 (14.8)1.00025 (13.3)
GI involvement, n (%)619365 (59)264 (62.1)0.037101 (52.1)
GERD, n (%)618344 (55.7)251 (59.2)0.01893 (47.9)
Stomach symptoms, n (%)617119 (19.3)88 (20.8)0.34231 (16.1)
Intestinal symptoms, n (%)619101 (16.3)68 (16.0)1.00033 (17.0)
SRC, n (%)61810 (1.6)8 (1.9)0.8682 (1.0)
Serum creatininea, mg/dl, mean (s.d.)580.83 (0.73)0.72 (0.16)0.2091.02 (1.19)
Proteinuria, n (%)58948 (8.1)38 (9.4)0.18910 (5.4)
ILD, n (%)598290 (48.5)203 (49.5)0.92587 (46.3)
Restrictive LD, n (%)579202 (34.9)163 (41.1)<0.000139 (21.4)
FVC, % predicted31790 (77–101)86 (75–99)<0.000195 (84.5–106)
TLCa, % predicted20190 (76–100)86 (73–99)0.00893 (81–107)
DLCO, % predicted48868 (54–82)64 (52–80)0.00172 (59–86)
FEV1a, % predicted23289 (78–100)87 (76–99)0.01594 (83–103)
PHa, n (%)2439 (3.7)8 (4.7)0.4471 (1.4)
sPAPa, mmHg, mean (s.d.)23328.7 (8.4)29.1 (8.9)0.38927.7 (6.7)
PMI, n (%)619134 (21.6)98 (23.1)0.41436 (18.6)
PMI (<45%), n (%)607127 (20.9)93 (22.4)0.40034 (17.8)
Ejection fractiona, %, mean (s.d.)29162.7 (6.2)62.8 (6)1.00062.5 (6.7)
Cardiac blocks, n (%)58157 (9.8)40 (10.1)1.00017 (9.3)
Diastolic dysfunction, n (%)57187 (15.2)65 (16.6)0.37122 (12.3)
Pericardial effusiona, n (%)30115 (5)11 (5.4)1.0004 (4.1)
Any severe, n (%)619351 (56.7)247 (58.1)0.587104 (53.6)
EScSG activity index 20016191 (0.5–2.5)1.5 (0.5–3)<0.00010.5 (0–1.5)
EScSG activity index 2016a524 (3.1–5.4)4.4 (3.8–5.6)0.0393.5 (1.5–5)
EScSG 2001 active, n (%)619103 (16.6)87 (20.5)<0.000116 (8.2)
EScSG 2016 activea, n (%)5244 (84.6)30 (90.9)0.09714 (73.7)
Corticosteroids, n (%)619228 (36.8)165 (38.8)0.25763 (32.5)
Prednisone equivalent, mg/day6195 (5–7.5)5 (5–7.5)0.0595 (4–6)
Immunosuppressants, n (%)619332 (53.6)242 (56.9)0.02990 (46.4)
CharacteristicsPatients with data available, nTotal ATA (N = 619)ATA-dcSSc (n = 425)P-valueATA-lcSSc (n = 194)
Female, n (%)619463 (74.8)302 (71.1)0.003161 (83)
Age at RP onset, years61545.3 (34.5–55.6)45.1 (35.4–55.3)0.83245.4 (31.4–56)
Time from RP to SSc, years6150.3 (0–2)0.1 (0–1.2)<0.00010.7 (0–3)
Age at SSc onset, years61947.2 (36.9–57)46.8 (36.9–56.3)1.00047.8 (37.3–57.8)
Age, years61948.9 (38.6–58.6)48.2 (38.6–58.2)1.00049.7 (38.4–59.3)
Disease duration, years, mean (s.d.)6191.4 (0.8)1.5 (0.8)1.0001.4 (0.9)
BMI, mean (s.d.)35724.7 (4.4)24.3 (4.2)0.03025.5 (4.7)
ESR >20 mm/h, n (%)330147 (44.5)111 (48.7)0.04736 (35.3)
CRP elevation, n (%)337109 (32.3)84 (36.4)0.04025 (23.6)
Hba, g/dl, mean (s.d.)5713.4 (0.2)13.4 (0.2)1.00013.3 (0.3)
mRSS58612.5 (5–20)16 (11–23)<0.00015 (2–8)
Puffy fingers (ever), n (%)336115 (34.2)79 (35)1.00036 (32.7)
DUs ever, n (%)616212 (34.4)153 (36)0.42659 (30.7)
DUs current, n (%)616190 (30.8)137 (32.3)0.48353 (27.6)
Teleangiectasiaa, n (%)10946 (42.2)36 (47.4)0.19510 (30.3)
Joint synovitis, n (%)615145 (23.6)112 (26.5)0.02433 (17.2)
TFRs, n (%)610112 (18.4)102 (24.4)<0.000111(5.7)
Muscle weakness, n (%)615145 (23.6)104 (24.7)0.66541 (21.1)
CK elevation, n (%)58083 (14.3)58 (14.8)1.00025 (13.3)
GI involvement, n (%)619365 (59)264 (62.1)0.037101 (52.1)
GERD, n (%)618344 (55.7)251 (59.2)0.01893 (47.9)
Stomach symptoms, n (%)617119 (19.3)88 (20.8)0.34231 (16.1)
Intestinal symptoms, n (%)619101 (16.3)68 (16.0)1.00033 (17.0)
SRC, n (%)61810 (1.6)8 (1.9)0.8682 (1.0)
Serum creatininea, mg/dl, mean (s.d.)580.83 (0.73)0.72 (0.16)0.2091.02 (1.19)
Proteinuria, n (%)58948 (8.1)38 (9.4)0.18910 (5.4)
ILD, n (%)598290 (48.5)203 (49.5)0.92587 (46.3)
Restrictive LD, n (%)579202 (34.9)163 (41.1)<0.000139 (21.4)
FVC, % predicted31790 (77–101)86 (75–99)<0.000195 (84.5–106)
TLCa, % predicted20190 (76–100)86 (73–99)0.00893 (81–107)
DLCO, % predicted48868 (54–82)64 (52–80)0.00172 (59–86)
FEV1a, % predicted23289 (78–100)87 (76–99)0.01594 (83–103)
PHa, n (%)2439 (3.7)8 (4.7)0.4471 (1.4)
sPAPa, mmHg, mean (s.d.)23328.7 (8.4)29.1 (8.9)0.38927.7 (6.7)
PMI, n (%)619134 (21.6)98 (23.1)0.41436 (18.6)
PMI (<45%), n (%)607127 (20.9)93 (22.4)0.40034 (17.8)
Ejection fractiona, %, mean (s.d.)29162.7 (6.2)62.8 (6)1.00062.5 (6.7)
Cardiac blocks, n (%)58157 (9.8)40 (10.1)1.00017 (9.3)
Diastolic dysfunction, n (%)57187 (15.2)65 (16.6)0.37122 (12.3)
Pericardial effusiona, n (%)30115 (5)11 (5.4)1.0004 (4.1)
Any severe, n (%)619351 (56.7)247 (58.1)0.587104 (53.6)
EScSG activity index 20016191 (0.5–2.5)1.5 (0.5–3)<0.00010.5 (0–1.5)
EScSG activity index 2016a524 (3.1–5.4)4.4 (3.8–5.6)0.0393.5 (1.5–5)
EScSG 2001 active, n (%)619103 (16.6)87 (20.5)<0.000116 (8.2)
EScSG 2016 activea, n (%)5244 (84.6)30 (90.9)0.09714 (73.7)
Corticosteroids, n (%)619228 (36.8)165 (38.8)0.25763 (32.5)
Prednisone equivalent, mg/day6195 (5–7.5)5 (5–7.5)0.0595 (4–6)
Immunosuppressants, n (%)619332 (53.6)242 (56.9)0.02990 (46.4)

Values are presented as median (IQR) unless stated otherwise.

a

Data available in less than half of patients.

Significant P-values (<0.05) are in bold.

CK: creatine kinase; DUs: digital ulcers; EScSG: European Scleroderma Study Group; FEV1: forced expiratory volume in the first second; GERD: gastroesophageal reflux disease; Hb: haemoglobin; LD: lung disease; sPAP: systolic pulmonary artery pressure; immunosuppressants: CYC, SSZ, HCQ, MTX, LEF, AZA, MMF, ciclosporin A, d-penicillamine, rituximab, imatinib, anti-TNF inhibitors, tocilizumab or abatacept; DLCO: diffusing capacity for carbon monoxide; FVC: forced vital capacity; GI: gastrointestinal; ILD: interstitial lung disease; mRSS: modified Rodnan Skin Score; PH: pulmonary hypertension; PMI: primary myocardial involvement; SRC: scleroderma renal crisis; TFRs: tendon friction rubs; TLC: total lung capacity. 

Demographic and clinical characteristics at baseline (Tables 1 and 2)

Group 1 (patients with lcSSc). ATA-lcSSc patients were more frequently male with a shorter time from RP to SSc onset (P < 0.0001) vs ACA-lcSSc, but similar to ANA-lcSSc patients (Table 1). In ATA-lcSSc patients, the median value of the modified Rodnan skin score (mRSS) was consistent with lcSSc at baseline [ 5 (IQR 2–8)] and at the end of follow-up [4 (IQR 0–6)]. At baseline, our group of interest had a higher frequency of increased CRP, articular tendon involvement [i.e. joints synovitis and tendon friction rubs (TFRs)] and muscular impairment (i.e. muscle weakness and elevated creatine kinase levels) than ACA-lcSSc patients but similar to ANA-lcSSc patients. The frequency of SRC, PH and PMI was similar among groups. In contrast, ILD occurred more frequently in ATA-lcSSc (46.3%) than in ANA-lcSSc (P = 0.004) and ACA-lcSSc (P < 0.0001). Restrictive lung disease was more frequent in ATA-lcSSc than in ACA-lcSSc patients (P < 0.0001) but was similar to ANA-lcSSc patients. The frequency of any major organ involvement was higher in ATA-lcSSc compared with ACA-lcSSc (P < 0.0001). The use of both glucocorticoids and immunosuppressants was higher in ATA-lcSSc than in ACA-lcSSc patients (P < 0.0001, for both), but was similar in ATA-lcSSc and ANA-lcSSc patients. Among ATA-lcSSc patients, methotrexate was the most frequently prescribed drug (28% of patients taking immunosuppressants), followed by azathioprine (19%) and mycophenolate mofetil (13%).

Group 2 (ATA-positive SSc patients). The percentage of females was higher and the time from RP to SSc onset was longer in ATA-lcSSc than in ATA-dcSSc patients (Table 2). An increase in serum inflammatory biomarkers, joint synovitis and TFRs was less frequently observed in ATA-lcSSc than in ATA-dcSSc patients, whereas the frequency of muscular impairment was similar in the two groups. At baseline, the frequency of SRC, PMI, PH and ILD (46.3% in lcSSc vs 49.5% in dcSSc) were similar between groups. Nevertheless, restrictive lung disease was less frequent in patients with ATA-lcSSc vs ATA-dcSSc (P < 0.0001), with higher values of FVC, total lung capacity (TLC) and diffusion capacity for carbon monoxide (DLCO) in the former. The occurrence of any major organ involvement was similar in both groups even though immunosuppressants (P = 0.029) were prescribed less frequently in lcSSc vs dcSSc patients.

Multivariable Cox regression models

The mean follow-up in the whole study population was 7.7 years (s.d. 3.5) and similar in all four groups, ranging from 7.2 years (s.d. 3.1) in ATA-lcSSc to 8.1 years (s.d. 3.9) in ACA-lcSSc, with a disease duration at last available follow-up ranging from 8.7 years (s.d. 3.2) to 9.5 (s.d. 4). A total of 136 patients (10.9%) died during follow-up: 16 (8.2%) ATA-lcSSc, 53 (9.9%) ACA-lcSSc, 6 (6.3%) ANA-lcSSc and 61 (14.4%) ATA-dcSSc. According to the clinical practice of expert centres, at least one HRCT was performed in 86% of patients with ATA-lcSSc, 86% with ATA-dcSSc and 74% with ACA-lcSSc and ANA-lcSSc.

Survival curves are reported in Fig. 1 and multivariable Cox regression models for the risk of developing major organ involvement and death are reported in Tables 3 and 4.

Kaplan–Meier survival curves for risk of major organ involvement and death in the four subsets
Fig. 1

Kaplan–Meier survival curves for risk of major organ involvement and death in the four subsets

Survival curves indicating the risk of (A) ILD), (B) PMI, (C) PH, (D) any major organ involvement and (E) death in the four subsets: ATA-lcSSc, ACA-lcSSc, ANA-lcSSc and ATA-dcSSc.

Table 3

Multivariable Cox regression models for the risk of major organ involvement in the limited forms

VariablesHRCIP-value
ILD (n = 854, 178 with an event)
vs ACA-lcSSc4.553.16, 6.29<0.0001
 ANA-lcSSc2.561.51, 4.27<0.0001
 Male1.100.56, 1.710.735
 Age1.021.01, 1.04<0.0001
 Smoke ever0.930.59, 1.450.612
PMI (n = 860, 336 with an event)
vs ACA-lcSSc1.371.05, 1.790.020
 ANA-lcSSc1.140.78, 1.690.479
 Male1.741.23, 2.460.001
 Age1.021.01, 1.03<0.0001
 Arterial hypertension1.250.96, 1.640.098
 BMI1.010.98, 1.050.387
PH (n = 856, 120 with an event)
vs ACA-lcSSc1.060.65, 1.750.803
 ANA-lcSSc0.940.49, 1.830.866
 Male0.820.41, 1.650.058
 Age1.051.03, 1.07<0.0001
 Arterial hypertension1.771.19, 2.640.005
Any major organ involvement (n = 860, 356 with an event)
vs ACA-lcSSc1.881.41, 2.490.0006
 ANA-lcSSc1.761.16, 2.690.008
 Male1.440.99, 2.080.072
 Age1.021.01, 1.03<0.0001
 Arterial hypertension1.190.90, 1.580.220
Death (n = 860, 75 with an event)
vs ACA-lcSSc1.350.75, 2.380.311
 ANA-lcSSc1.630.62, 4.160.314
 Age1.081.06, 1.11<0.0001
 Male1.610.75, 3.450.208
VariablesHRCIP-value
ILD (n = 854, 178 with an event)
vs ACA-lcSSc4.553.16, 6.29<0.0001
 ANA-lcSSc2.561.51, 4.27<0.0001
 Male1.100.56, 1.710.735
 Age1.021.01, 1.04<0.0001
 Smoke ever0.930.59, 1.450.612
PMI (n = 860, 336 with an event)
vs ACA-lcSSc1.371.05, 1.790.020
 ANA-lcSSc1.140.78, 1.690.479
 Male1.741.23, 2.460.001
 Age1.021.01, 1.03<0.0001
 Arterial hypertension1.250.96, 1.640.098
 BMI1.010.98, 1.050.387
PH (n = 856, 120 with an event)
vs ACA-lcSSc1.060.65, 1.750.803
 ANA-lcSSc0.940.49, 1.830.866
 Male0.820.41, 1.650.058
 Age1.051.03, 1.07<0.0001
 Arterial hypertension1.771.19, 2.640.005
Any major organ involvement (n = 860, 356 with an event)
vs ACA-lcSSc1.881.41, 2.490.0006
 ANA-lcSSc1.761.16, 2.690.008
 Male1.440.99, 2.080.072
 Age1.021.01, 1.03<0.0001
 Arterial hypertension1.190.90, 1.580.220
Death (n = 860, 75 with an event)
vs ACA-lcSSc1.350.75, 2.380.311
 ANA-lcSSc1.630.62, 4.160.314
 Age1.081.06, 1.11<0.0001
 Male1.610.75, 3.450.208

Patients with an event already at baseline were excluded.

Significant P values (<0.05) are in bold. CI: confidence interval; HR: hazard ratio; ILD: interstitial lung disease; PMI: primary myocardial involvement; PH: pulmonary hypertension.

Table 3

Multivariable Cox regression models for the risk of major organ involvement in the limited forms

VariablesHRCIP-value
ILD (n = 854, 178 with an event)
vs ACA-lcSSc4.553.16, 6.29<0.0001
 ANA-lcSSc2.561.51, 4.27<0.0001
 Male1.100.56, 1.710.735
 Age1.021.01, 1.04<0.0001
 Smoke ever0.930.59, 1.450.612
PMI (n = 860, 336 with an event)
vs ACA-lcSSc1.371.05, 1.790.020
 ANA-lcSSc1.140.78, 1.690.479
 Male1.741.23, 2.460.001
 Age1.021.01, 1.03<0.0001
 Arterial hypertension1.250.96, 1.640.098
 BMI1.010.98, 1.050.387
PH (n = 856, 120 with an event)
vs ACA-lcSSc1.060.65, 1.750.803
 ANA-lcSSc0.940.49, 1.830.866
 Male0.820.41, 1.650.058
 Age1.051.03, 1.07<0.0001
 Arterial hypertension1.771.19, 2.640.005
Any major organ involvement (n = 860, 356 with an event)
vs ACA-lcSSc1.881.41, 2.490.0006
 ANA-lcSSc1.761.16, 2.690.008
 Male1.440.99, 2.080.072
 Age1.021.01, 1.03<0.0001
 Arterial hypertension1.190.90, 1.580.220
Death (n = 860, 75 with an event)
vs ACA-lcSSc1.350.75, 2.380.311
 ANA-lcSSc1.630.62, 4.160.314
 Age1.081.06, 1.11<0.0001
 Male1.610.75, 3.450.208
VariablesHRCIP-value
ILD (n = 854, 178 with an event)
vs ACA-lcSSc4.553.16, 6.29<0.0001
 ANA-lcSSc2.561.51, 4.27<0.0001
 Male1.100.56, 1.710.735
 Age1.021.01, 1.04<0.0001
 Smoke ever0.930.59, 1.450.612
PMI (n = 860, 336 with an event)
vs ACA-lcSSc1.371.05, 1.790.020
 ANA-lcSSc1.140.78, 1.690.479
 Male1.741.23, 2.460.001
 Age1.021.01, 1.03<0.0001
 Arterial hypertension1.250.96, 1.640.098
 BMI1.010.98, 1.050.387
PH (n = 856, 120 with an event)
vs ACA-lcSSc1.060.65, 1.750.803
 ANA-lcSSc0.940.49, 1.830.866
 Male0.820.41, 1.650.058
 Age1.051.03, 1.07<0.0001
 Arterial hypertension1.771.19, 2.640.005
Any major organ involvement (n = 860, 356 with an event)
vs ACA-lcSSc1.881.41, 2.490.0006
 ANA-lcSSc1.761.16, 2.690.008
 Male1.440.99, 2.080.072
 Age1.021.01, 1.03<0.0001
 Arterial hypertension1.190.90, 1.580.220
Death (n = 860, 75 with an event)
vs ACA-lcSSc1.350.75, 2.380.311
 ANA-lcSSc1.630.62, 4.160.314
 Age1.081.06, 1.11<0.0001
 Male1.610.75, 3.450.208

Patients with an event already at baseline were excluded.

Significant P values (<0.05) are in bold. CI: confidence interval; HR: hazard ratio; ILD: interstitial lung disease; PMI: primary myocardial involvement; PH: pulmonary hypertension.

Table 4

Multivariable Cox regression models for the risk of major organ involvement in ATA positive patients

VariablesHRCIP-value
ILD (n = 617, 192 with an event)
vs ATA-dcSSc0.750.54, 1.050.092
 Male0.950.64, 1.390.781
 Age1.011.0, 1.020.05
 Smoke ever0.980.57, 1.690.954
PMI (n = 619, 272 with an event)
vs ATA-dcSSc0.700.54, 0.920.011
 Male1.130.85, 1.490.409
 Age0.980.97, 0.99<0.0001
 Arterial hypertension1.050.76, 1.480.739
 BMI1.030.99, 1.060.075
PH (n = 615, 116 with an event)
vs ATA-dcSSc0.430.27, 0.69<0.0001
 Male1.320.86, 2.010.196
 Age1.041.03, 1.06<0.0001
 Arterial hypertension1.180.75, 1.850.748
Any major organ involvement (n = 619, 214 with an event)
vs ATA-dcSSc0.660.49, 0.880.006
 Male1.270.90, 1.780.171
 Age1.011.003, 1.020.008
 Arterial hypertension0.820.51, 1.290.391
Death (n = 619, 77 with an event)
vs ATA-dcSSc0.570.33, 1.010.053
 Age1.051.03, 1.07<0.0001
 Male1.81.11, 2.910.01
VariablesHRCIP-value
ILD (n = 617, 192 with an event)
vs ATA-dcSSc0.750.54, 1.050.092
 Male0.950.64, 1.390.781
 Age1.011.0, 1.020.05
 Smoke ever0.980.57, 1.690.954
PMI (n = 619, 272 with an event)
vs ATA-dcSSc0.700.54, 0.920.011
 Male1.130.85, 1.490.409
 Age0.980.97, 0.99<0.0001
 Arterial hypertension1.050.76, 1.480.739
 BMI1.030.99, 1.060.075
PH (n = 615, 116 with an event)
vs ATA-dcSSc0.430.27, 0.69<0.0001
 Male1.320.86, 2.010.196
 Age1.041.03, 1.06<0.0001
 Arterial hypertension1.180.75, 1.850.748
Any major organ involvement (n = 619, 214 with an event)
vs ATA-dcSSc0.660.49, 0.880.006
 Male1.270.90, 1.780.171
 Age1.011.003, 1.020.008
 Arterial hypertension0.820.51, 1.290.391
Death (n = 619, 77 with an event)
vs ATA-dcSSc0.570.33, 1.010.053
 Age1.051.03, 1.07<0.0001
 Male1.81.11, 2.910.01

Patients with an event already at baseline were excluded.

Significant P-values (<0.05) are in bold. CI: confidence interval; HR: hazard ratio; ILD: interstitial lung disease; PMI: primary myocardial involvement; PH: pulmonary hypertension.

Table 4

Multivariable Cox regression models for the risk of major organ involvement in ATA positive patients

VariablesHRCIP-value
ILD (n = 617, 192 with an event)
vs ATA-dcSSc0.750.54, 1.050.092
 Male0.950.64, 1.390.781
 Age1.011.0, 1.020.05
 Smoke ever0.980.57, 1.690.954
PMI (n = 619, 272 with an event)
vs ATA-dcSSc0.700.54, 0.920.011
 Male1.130.85, 1.490.409
 Age0.980.97, 0.99<0.0001
 Arterial hypertension1.050.76, 1.480.739
 BMI1.030.99, 1.060.075
PH (n = 615, 116 with an event)
vs ATA-dcSSc0.430.27, 0.69<0.0001
 Male1.320.86, 2.010.196
 Age1.041.03, 1.06<0.0001
 Arterial hypertension1.180.75, 1.850.748
Any major organ involvement (n = 619, 214 with an event)
vs ATA-dcSSc0.660.49, 0.880.006
 Male1.270.90, 1.780.171
 Age1.011.003, 1.020.008
 Arterial hypertension0.820.51, 1.290.391
Death (n = 619, 77 with an event)
vs ATA-dcSSc0.570.33, 1.010.053
 Age1.051.03, 1.07<0.0001
 Male1.81.11, 2.910.01
VariablesHRCIP-value
ILD (n = 617, 192 with an event)
vs ATA-dcSSc0.750.54, 1.050.092
 Male0.950.64, 1.390.781
 Age1.011.0, 1.020.05
 Smoke ever0.980.57, 1.690.954
PMI (n = 619, 272 with an event)
vs ATA-dcSSc0.700.54, 0.920.011
 Male1.130.85, 1.490.409
 Age0.980.97, 0.99<0.0001
 Arterial hypertension1.050.76, 1.480.739
 BMI1.030.99, 1.060.075
PH (n = 615, 116 with an event)
vs ATA-dcSSc0.430.27, 0.69<0.0001
 Male1.320.86, 2.010.196
 Age1.041.03, 1.06<0.0001
 Arterial hypertension1.180.75, 1.850.748
Any major organ involvement (n = 619, 214 with an event)
vs ATA-dcSSc0.660.49, 0.880.006
 Male1.270.90, 1.780.171
 Age1.011.003, 1.020.008
 Arterial hypertension0.820.51, 1.290.391
Death (n = 619, 77 with an event)
vs ATA-dcSSc0.570.33, 1.010.053
 Age1.051.03, 1.07<0.0001
 Male1.81.11, 2.910.01

Patients with an event already at baseline were excluded.

Significant P-values (<0.05) are in bold. CI: confidence interval; HR: hazard ratio; ILD: interstitial lung disease; PMI: primary myocardial involvement; PH: pulmonary hypertension.

ILD development. A total of 304 of 950 patients without ILD at baseline developed ILD during follow-up. ILD risk was nearly 5-fold higher in ATA-lcSSc patients than in ACA-lcSSc and 2.5-fold higher than in ANA-lcSSc (Table 3). In contrast, it was similar between ATA-lcSSc and ATA-dcSSc (Table 4).

ILD progression (Table 5; Supplementary Table S1 and S2, available at Rheumatology online). Among the 194 ATA-lcSSc patients, 143 (74%) had ILD ‘ever’ (i.e. at baseline or during follow-up). After the exclusion of patients with no available FVC at the ILD diagnosis or with only one FVC measurement after the ILD diagnosis, data on progression were available in 83 cases (58%). Among them, 29/83 (35%) had an FVC decline ≥10% and 48/83 (58%) had an FVC decline ≥5% (Supplementary Table S1, available at Rheumatology online).

Table 5

Multivariable Cox regression models for the risk of ILD progression

VariablesHRCIP-value
FVC decline ≥10%
Group 1: n = 196, 81 with an event
Group 2: n = 289, 138 with an event
vs ACA-lcSSc0.630.37, 1.080.093
 ANA-lcSSc0.860.42, 1.790.678
 Male0.930.43, 1.990.856
 Age1.010.98, 1.010.537
 Smoke ever1.170.60, 2.260.635
 FVC_01.021.01, 1.030.004
 IS_therapy_lung2.121.23, 3.650.007
vs ATA-dcSSc0.610.40, 0.950.028
 Age1.000.99, 1.020.296
 Male1.090.71, 1.660.696
 Smoke ever1.320.77, 2.260.301
 FVC_00.990.98, 1.000.394
 IS_therapy_lung1.160.82, 1.650.402
FVC decline ≥5%
Group 1: n = 196, 112 with an event
Group 2: n = 289, 189 with an event
vs ACA-lcSSc0.860.55, 1.350.500
 ANA-lcSSc1.250.68, 2.340.469
 Male0.950.52, 1.760.954
 Age0.990.98, 1.010.894
 Smoke ever1.030.54, 1.990.915
 FVC_01.010.99, 1.020.054
 IS_therapy_lung1.290.80, 2.070.286
vs ATA-dcSSc0.730.51, 1.010.060
 Male1.120.77, 1.600.564
 Age1.000.99, 1.010.760
 Smoke ever0.980.56, 1.750.963
 FVC_00.990.99, 1.010.756
 IS_therapy_lung0.890.65, 1.220.474
VariablesHRCIP-value
FVC decline ≥10%
Group 1: n = 196, 81 with an event
Group 2: n = 289, 138 with an event
vs ACA-lcSSc0.630.37, 1.080.093
 ANA-lcSSc0.860.42, 1.790.678
 Male0.930.43, 1.990.856
 Age1.010.98, 1.010.537
 Smoke ever1.170.60, 2.260.635
 FVC_01.021.01, 1.030.004
 IS_therapy_lung2.121.23, 3.650.007
vs ATA-dcSSc0.610.40, 0.950.028
 Age1.000.99, 1.020.296
 Male1.090.71, 1.660.696
 Smoke ever1.320.77, 2.260.301
 FVC_00.990.98, 1.000.394
 IS_therapy_lung1.160.82, 1.650.402
FVC decline ≥5%
Group 1: n = 196, 112 with an event
Group 2: n = 289, 189 with an event
vs ACA-lcSSc0.860.55, 1.350.500
 ANA-lcSSc1.250.68, 2.340.469
 Male0.950.52, 1.760.954
 Age0.990.98, 1.010.894
 Smoke ever1.030.54, 1.990.915
 FVC_01.010.99, 1.020.054
 IS_therapy_lung1.290.80, 2.070.286
vs ATA-dcSSc0.730.51, 1.010.060
 Male1.120.77, 1.600.564
 Age1.000.99, 1.010.760
 Smoke ever0.980.56, 1.750.963
 FVC_00.990.99, 1.010.756
 IS_therapy_lung0.890.65, 1.220.474

FVC_0: FVC at ILD onset; IS_therapy_lung: immunosuppressants with efficacy on lung involvement (i.e. MMF, AZA, tocilizumab, rituximab, CYC).

Significant P-values (<0.05) are in bold. CI: confidence interval; FVC: forced vital capacity; HR: hazard ratio.

Table 5

Multivariable Cox regression models for the risk of ILD progression

VariablesHRCIP-value
FVC decline ≥10%
Group 1: n = 196, 81 with an event
Group 2: n = 289, 138 with an event
vs ACA-lcSSc0.630.37, 1.080.093
 ANA-lcSSc0.860.42, 1.790.678
 Male0.930.43, 1.990.856
 Age1.010.98, 1.010.537
 Smoke ever1.170.60, 2.260.635
 FVC_01.021.01, 1.030.004
 IS_therapy_lung2.121.23, 3.650.007
vs ATA-dcSSc0.610.40, 0.950.028
 Age1.000.99, 1.020.296
 Male1.090.71, 1.660.696
 Smoke ever1.320.77, 2.260.301
 FVC_00.990.98, 1.000.394
 IS_therapy_lung1.160.82, 1.650.402
FVC decline ≥5%
Group 1: n = 196, 112 with an event
Group 2: n = 289, 189 with an event
vs ACA-lcSSc0.860.55, 1.350.500
 ANA-lcSSc1.250.68, 2.340.469
 Male0.950.52, 1.760.954
 Age0.990.98, 1.010.894
 Smoke ever1.030.54, 1.990.915
 FVC_01.010.99, 1.020.054
 IS_therapy_lung1.290.80, 2.070.286
vs ATA-dcSSc0.730.51, 1.010.060
 Male1.120.77, 1.600.564
 Age1.000.99, 1.010.760
 Smoke ever0.980.56, 1.750.963
 FVC_00.990.99, 1.010.756
 IS_therapy_lung0.890.65, 1.220.474
VariablesHRCIP-value
FVC decline ≥10%
Group 1: n = 196, 81 with an event
Group 2: n = 289, 138 with an event
vs ACA-lcSSc0.630.37, 1.080.093
 ANA-lcSSc0.860.42, 1.790.678
 Male0.930.43, 1.990.856
 Age1.010.98, 1.010.537
 Smoke ever1.170.60, 2.260.635
 FVC_01.021.01, 1.030.004
 IS_therapy_lung2.121.23, 3.650.007
vs ATA-dcSSc0.610.40, 0.950.028
 Age1.000.99, 1.020.296
 Male1.090.71, 1.660.696
 Smoke ever1.320.77, 2.260.301
 FVC_00.990.98, 1.000.394
 IS_therapy_lung1.160.82, 1.650.402
FVC decline ≥5%
Group 1: n = 196, 112 with an event
Group 2: n = 289, 189 with an event
vs ACA-lcSSc0.860.55, 1.350.500
 ANA-lcSSc1.250.68, 2.340.469
 Male0.950.52, 1.760.954
 Age0.990.98, 1.010.894
 Smoke ever1.030.54, 1.990.915
 FVC_01.010.99, 1.020.054
 IS_therapy_lung1.290.80, 2.070.286
vs ATA-dcSSc0.730.51, 1.010.060
 Male1.120.77, 1.600.564
 Age1.000.99, 1.010.760
 Smoke ever0.980.56, 1.750.963
 FVC_00.990.99, 1.010.756
 IS_therapy_lung0.890.65, 1.220.474

FVC_0: FVC at ILD onset; IS_therapy_lung: immunosuppressants with efficacy on lung involvement (i.e. MMF, AZA, tocilizumab, rituximab, CYC).

Significant P-values (<0.05) are in bold. CI: confidence interval; FVC: forced vital capacity; HR: hazard ratio.

The risk of FVC decline ≥10% was similar in ATA-lcSSc and other limited forms but lower in ATA-lcSSc than in ATA-dcSSc patients. The risk of FVC decline ≥5% was similar between ATA-lcSSc and all the other subsets, including ATA-dcSSc patients (Table 5, Supplementary Table S2, available at Rheumatology online).

In all patients, the period of observation after ILD onset was 3.6 years (IQR 2–5.5) for an FVC decline ≥10% and 2.7 years (IQR 1.3–4.6) for an FVC decline ≥5%.

PMI. A total of 515 of 1000 patients without PMI at baseline developed PMI during follow-up. The risk of PMI in ATA-lcSSc was higher than in ACA-lcSSc but similar to that in ANA-lcSSc (Table 3). Patients with ATA-lcSSc had a lower risk of PMI vs ATA-dcSSc (Table 4).

PH. A total of 213 of 1118 patients without PH at baseline developed PH during follow-up. ATA-lcSSc patients had a risk of PH similar to both ACA-lcSSc and ANA-lcSSc patients but lower than ATA-dcSSc patients (Tables 3 and 4).

Any major organ involvement. A total of 486 of 712 patients without any major organ involvement at baseline developed it during follow-up. The risk of any major organ involvement in ATA-lcSSc was higher than in ACA-lcSSc and ANA-lcSSc (Table 3) and lower than in ATA-dcSSc (Table 4).

All-cause mortality (136 events). The risk of death from all causes in ATA-lcSSc patients was similar to that observed in ACA-lcSSc and ANA-lcSSc patients and tended to be lower than in ATA-dcSSc patients [HR 0.57 (95% CI 0.33–1.01), P = 0.053].

Discussion

Given the wide heterogenicity of SSc, several efforts have been made in the last few years to improve patient stratification, with a renewed interest in the potential role of SSc-specific autoantibodies in predicting outcome [1, 6, 33].

ATA has been traditionally associated with diffuse skin involvement, aggressive disease and poor prognosis. However, recent reports have highlighted a wide spectrum of disease severity among ATA-positive patients [34]. This might be due to the higher sensitivity of ACR/EULAR 2013 criteria in identifying less severe SSc forms and the different titre and Ig classes of ATA. Indeed, it has been observed that high titres of IgG and the presence of IgM ATA are associated with a more progressive disease phenotype [35, 36].

Although ATA-positive lcSSc patients may be considered as a non-classic SSc subset, they are reported in ∼10–11% of recent SSc cohorts [33]. We investigated the largest number of patients with ATA-lcSSc and early disease at baseline. This allowed us a more precise characterization and risk stratification of the ATA-lcSSc clinical phenotype over time compared with previous studies [13, 37]. We compared ATA-lcSSc patients with those affected with other limited forms (ACA-lcSSc and ANA-lcSSc) and with the ‘classical’ ATA-dcSSc subset. Our subclassification differs from previous studies [13, 14] where ATA- and ACA-positive or negative lcSSc and dcSSc patients were considered irrespective of other antibody specificities (e.g. anti-RNA polymerase III, anti-PM/Scl). Given the emerging prognostic role of some autoantibodies in SSc [38, 39] and the potential coexistence of ACA and ATA in the same patient [40], the inclusion or exclusion of a specific autoantibody reactivity might have a great influence on the patient’s characterization.

In contrast to Kranenburg et al. [14], we found a higher frequency of ILD at baseline in ATA-lcSSc vs other lcSSc groups, highlighting a very close association between ATA and ILD since the early stages of SSc. Even during follow-up, the risk of ILD development was similar in ATA-lcSSc patients and ATA-dcSSc, suggesting that all ATA-positive patients should be closely monitored for ILD over time, regardless of the skin subset, as has also been recommended by the recent European consensus on SSc-ILD [21].

In our study population, ATA-lcSSc appears to have a more inflammatory phenotype (i.e. elevated CRP, musculoskeletal involvement) than ACA-lcSSc, but similar to ANA-lcSSc-patients, where few data are available in the literature. These findings might help to identify lcSSc patients with ILD who may benefit from biologics. In fact, tocilizumab is more effective in ATA-positive vs ATA-negative dcSSc patients, raising the idea of investigating this drug in ATA-lcSSc and ANA-lcSSc patients [41].

The risk of PMI in ATA-lcSSc was between that of ACA-lcSSc and ATA-dcSSc, which may be expected since both ATA and dcSSc have been recognized as risk factors for PMI-SSc [42]. Probably due to the different definitions adopted, the PMI risk was higher in our ATA-lcSSc population than in other cohorts [33]. Similar to Kranenburg et al. [14] and Nythyanova et al. [33], we found that the risk of PH was similar among the limited forms irrespective of antibody specificity. On the other hand, the higher risk of PH in ATA-dcSSc vs ATA-lcSSc may be due to the more frequent PMI in the former, given the high prevalence of post-capillary PH vs PAH in SSc [43].

The overall survival rate in our ATA-lcSSc patients was similar to that of the other limited forms, confirming that ATA-positive patients may not invariably have a bad prognosis, as recently highlighted [33, 34]. Compared with dcSSc patients, we found a 15 year survival of 81% in ATA-lcSSc vs 64% in ATA-dcSSc. Interestingly, in the study by Perera et al. [35], published in 2007, 10 year mortality was good in ATA-lcSSc (86%) but remarkably poor in ATA-dcSSc (40%). Thus our results are probably due to the therapeutic advances in dcSSc management over the past 2 decades.

One of the most intriguing observations is the good survival rate of ATA-lcSSc patients despite their high frequency of ILD [14]. Due to the great heterogenicity of SSc-ILD severity and progression, and the evidence that nintedanib reduces the rate of FVC decline in patients with SSc-ILD and chronic progressive fibrosing ILD, the identification of an ILD progressive phenotype has become a priority in the research agenda [23]. To the best of our knowledge, this is the first study assessing the risk of progressive ILD in patients with ATA-lcSSc. At baseline, restrictive lung disease was less common in ATA-lcSSc vs ATA-dcSSc and FVC was higher in the former at the time of ILD diagnosis, suggesting a milder lung involvement in ATA-lcSSc at ILD onset. ATA-lcSSc patients had a risk of ILD progression based on an FVC decline ≥10% lower than that of ATA-dcSSc and similar to the other limited forms. Conflicting data have been reported on the role of ATA in predicting lung disease progression, with most studies evaluating FVC decline at 1 year [44]. A recent EUSTAR study [23] showed that a high skin score, but not positive ATA, was an independent predictor of FVC decline over 5 years. It has been argued that the cut-off ≥10% of FVC decline was derived from studies on idiopathic pulmonary fibrosis and may not be appropriate for SSc-ILD [45, 46]. By adopting the 5% cut-off, we identified twice as many cases of progressive ILD in ATA-lcSSc (59%), a proportion similar to that observed in ATA-dcSSc.

Besides ATA-lcSSc, our study provided new insights into the characterization and risk stratification of ANA-lcSSc patients: >40% had ILD over time and the severity of ILD as well as the risk of PH and PMI were similar to those of ATA-lcSSc. The EUSTAR group has recently developed recommendations for the annual monitoring of SSc, providing a common minimal core set of investigations for all SSc patients [47]. Due to the great heterogeneity of the disease, our findings could help to define a more tailored follow-up in different SSc subgroups.

Among the several strengths of our study, we used a large multicentre SSc cohort and included only patients with early disease at database entry, well-characterized skin involvement and autoantibody profile. The follow-up of at least 3 years allows us to minimize any misclassification of the cutaneous form, since the transition from lcSSc to dcSSc has been reported in the first 3 years of SSc in the great majority of cases [14, 48]; moreover, the percentage of transitions was low in our cohort, in line with other recent reports [46]. Unlike previous studies, we included several important confounders in the models for the risk of organ involvement, such as systemic arterial hypertension in the PMI model, given the importance of considering common cardiovascular comorbidities in the definition of SSc-PMI [49]. As a limitation of the study, the mortality risk may have been slightly underestimated, in particular in the ATA-dcSSc cohort; however, the early mortality in dcSSc patients is limited nowadays: 6% at 3 years in a recent SSc cohort [48]. Being a multicentre study, the autoantibody assays were not all done in the same laboratory, but all EUSTAR centres use standardized and internationally validated assays. Given the paucity of RHC data in the database, the presence of PH was also defined by TTE; however, the cut-off we adopted for a PH diagnosis was shown to be highly correlated with RHC [50]. Finally, since we considered FVC within 1 year from the ILD diagnosis as baseline, we cannot exclude that a few patients might have already progressed at that point. However, this reflects the fact that it is often not possible to perform pulmonary function tests soon after the radiological diagnosis in real-life settings. We did not consider severe GI involvement in the model for any major organ involvement, due to the lack of data to define it.

In conclusion, we provide a better understanding of the clinical phenotype and risk of major organ involvement in ATA-lcSSc patients. Due to the observed high risk of ILD, albeit with a lower risk of progression compared with ATA-dcSSc, our study supports careful screening for ILD in ATA-lcSSc both at baseline and during follow-up. A significant risk of cardiopulmonary involvement was also highlighted in ANA-lcSSc patients that deserves more investigations and may suggest close monitoring in these patients as well. Our study shed more light on the many faces of SSc and provides insights towards a more precise risk stratification of SSc patients.

Acknowledgements

EUSTAR collaborators: Giovanna Cuomo (Naples, Italy), Gianluca Moroncini (Ancona, Italy), Jiri Stork (Prague, Czech Republic), Fiorenzo Iannone (Bari, Italy), Ulrich Walker (Basel, Switzerland), Eugenia Bertoldo (Verona, Italy), Dorota Krasowska (Lublin, Poland), Maria João Salvador (Coimbra, Portugal), Mohammed Tikly (Johannesburg, South Africa), Eric Hachulla (Lille, France), Valeria Riccieri (Rome, Italy), Ami Sha (Baltimore, MD, USA), Ana Maria Gheorghiu (Bucharest, Romania), Cord Sunderkötter (Münster, Germany), Francesca Ingegnoli (Milan, Italy), Luc Mouthon (Paris, France), Vanessa Smith (Gent, Belgium), Francesco Paolo Cantatore (Foggia, Italy), Kilian Eyerich (Munich, Germany), Piotr Wiland (Wroclaw, Poland), Marie Vanthuyne (Brussels, Belgium), Branimir Anic (Zagreb, Croatia), Maria Üprus (Tallin, Estonia), Brigitte Granel (Marseille, France), Alessandra Vacca (Monserrato, Italy), Cristina-Mihaela Tanaseanu (Bucharest, Romania), Paloma García de la Peña Lefebvre (Madrid, Spain), Jean Sibilia (Strasbourg, France), Ira Litinsky (Tel Aviv, Israel), Lesley Ann Saketkoo (New Orleans, LA, USA), Eduardo Kerzberg (Buenos Aires, Argentina), Massimiliano Limonta (Bergamo, Italy), Doron Rimar (Haifa, Israel), Petros Sfikakis (Athens, Greece), Maurizio Cutolo (Genova, Italy), Patricia E. Carreira (Madrid, Spain), Rosario Foti (Catania, Italy), Srdan Novak (Rijeka, Croatia), Michele Iudici (Geneva, Switzerland), Mislav Radic (Split, Croatia), Raffaele Pellerito (Torino, Italy), Carlo Francesco Selmi Rozzano (Milan, Italy), Lidia P. Ananieva (Moscow, Russia), Gabriela Szücs (Debrecen, Hungary), Carlos de la Puente (Madrid, Spain), Ruxandra Maria Ionescu (Bucharest, Romania), Jörg Distler (Erlangen, Germany), Maria Rosa Pozzi (Monza, Italy), Juan Jose Alegre-Sancho (Valencia, Spain), Kristine Herrmann (Dresden, Germany), Ellen De Langhe (Leuven, Belgium), Sule Yavuz Altunizade (Istanbul, Turkey), Carolina de Souza Müller (Curitiba, Brazil), Svetlana Agachi (Chisinau, Republic of Moldova), Douglas Veale (Dublin, Ireland), Esthela Loyo (Santiago, Dominican Republic), Mengtao Li (Beijing, China), Edoardo Rosato (Rome, Italy), Britta Maurer (Bern, Switzerland), Ivan Castellví (Barcelona, Spain), François Spertini (Lausanne, Switzerland), Kamal Solanki (Hamilton, New Zealand), Nicoletta Del Papa (Milan, Italy), Gerard Espinosa (Barcelona, Spain), László Czirják (Pecs, Hungary), Bernard Coleiro (Balzan, Malta), Dominique Farge Bancel (Paris, France), Raffaele Pellerito (Torino, Italy), Christopher Denton (London, UK), Nemanja Damjanov (Belgrade, Serbia and Montenegro), Jörg Henes (Tübingen, Germany), Vera Ortiz Santamaria Granollers (Barcelona, Spain), Michaela Kohm (Frankfurt am Main, Germany), Bojana Stamenkovic (Niska Banja, Serbia and Montenegro).

Funding: No specific funding was received from any bodies in the public, commercial or not-for-profit sectors to carry out the work described in this article.

Disclosure statement: E.Z. has received speaker honoraria or consultancy fees from Actelion Pharmaceuticals, GlaxoSmithKline and Boehringer Ingelheim. J.A. had consultancy relationships or has received research funding from Actelion, Boehringer Ingelheim, Roche, Pfizer, Sanofi and Bristol-Myers Squibb. P.A. has received personal fees from Bristol-Myers Squibb and Boehringer Ingelheim and non-financial support from CSL Behring, SOBI, Janssen, Roche and Sanofi. A.M.H.V. had grant/research support from Boehringer Ingelheim and received speaker and personal fees from Boehringer Ingelheim, Roche, Actelion, MSD, Medscape and Bayer. O.D. had consultancy relationships and/or has received research funding from AbbVie, Acceleron Pharma, Amgen, AnaMar, Arxx Therapeutics, Baecon Discovery, Blade Therapeutics, Bayer, Boehringer Ingelheim, ChemomAb, Corbus Pharmaceuticals, CSL Behring, Galapagos NV, Glenmark Pharmaceuticals, GlaxoSmithKline, Horizon (Curzion) Pharmaceuticals, Inventiva, iQvia, Italfarmaco, iQone, Kymera Therapeutics, Eli Lilly, Medac, Medscape, Mitsubishi Tanabe Pharma, MSD, Novartis, Pfizer, Roche, Sanofi, Serodapharm, Topadur, Target Bioscience and UCB and a patent issued ‘mir-29 for the treatment of systemic sclerosis’ (US8247389, EP2331143). Y.A. had consultancy relationships and/or has received research funding from Actelion, Bayer, Boehringer Ingelheim, Genentech/Roche, Inventiva, Medsenic and Sanofi. The remaining authors have no competing interests to declare.

Data availability statement

The data that support the findings of this study are available upon reasonable request.

Supplementary data

Supplementary data are available at Rheumatology online.

References

1

Allanore
Y
,
Simms
R
,
Distler
O
et al.
Systemic sclerosis
.
Nat Rev Dis Prim
2015
;
1
:
15002
.

2

Zanatta
E
,
Codullo
V
,
Avouac
J
,
Allanore
Y.
Systemic sclerosis: recent insight in clinical management
.
Joint Bone Spine
2020
;
87
:
293
9
.

3

Fasano
S
,
Riccardi
A
,
Messiniti
V
et al.
Revised European Scleroderma Trials and Research Group Activity Index is the best predictor of short-term severity accrual
.
Ann Rheum Dis
2019
;
78
:
1681
5
.

4

van den Hombergh
WM
,
Knaapen-Hans
HK
,
van den Hoogen
FH
et al.
Prediction of organ involvement and survival in systemic sclerosis patients in the first 5 years from diagnosis
.
J Sclerod Related Disord
2020
;
5
:
57
65
.

5

LeRoy
EC
,
Black
C
,
Fleischmajer
R
et al.
Scleroderma (systemic sclerosis): classification, subsets and pathogenesis
.
J Rheumatol
1988
;
15
:
202
5
.

6

Sobanski
V
,
Giovannelli
J
,
Allanore
Y
et al.
Phenotypes determined by cluster analysis and their survival in the prospective European Scleroderma Trials and Research cohort of patients with systemic sclerosis
.
Arthritis Rheumatol
2019
;
71
:
1553
70
.

7

Mahler
M
,
Hudson
M
,
Bentow
C
et al.
Autoantibodies to stratify systemic sclerosis patients into clinically actionable subsets
.
Autoimmun Rev
2020
;
19
:
102583
.

8

Steen
VD
,
Powell
DL
,
Medsger
TA
Jr.
Clinical correlations and prognosis based on serum autoantibodies in patients with systemic sclerosis
.
Arthritis Rheum
1988
;
31
:
196
203
.

9

Steen
VD.
The many faces of scleroderma
.
Rheum Dis Clin North Am
2008
;
34
:
1
15
.

10

Diab
S
,
Dostrovsky
N
,
Hudson
M
et al.
Systemic sclerosis sine scleroderma: a multicener study of 1417 subjects
.
J Rheumatol
2014
;
41
:
2179
85
.

11

Igusa
T
,
Hummers
LK
,
Visvanathan
K
et al.
Autoantibodies and scleroderma phenotype define subgroups at high-risk and low-risk for cancer
.
Ann Rheum Dis
2018
;
77
:
1180
7
.

12

Walker
UA
,
Tyndall
A
,
Czirják
L
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
:
754
63
.

13

Srivastava
N
,
Hudson
M
,
Tatibouet
S
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
:
184
9
.

14

Kranenburg
P
,
van den Hombergh
WM
,
Knaapen-Hans
HK
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
:
2001
8
.

15

Elhai
M
,
Meune
C
,
Boubaya
M
et al.
Mapping and predicting mortality from systemic sclerosis
.
Ann Rheum Dis
2017
;
76
:
1897
905
.

16

Hoffmann-Vold
AM
,
Fretheim
H
,
Halse
AK
et al.
Tracking impact of interstitial lung disease in systemic sclerosis in a complete nationwide cohort
.
Am J Respir Crit Care Med
2019
;
200
:
1258
66
.

17

Vonk
MC
,
Walker
UA
,
Volkmann
ER
et al.
Natural variability in the disease course of SSc-ILD: implications for treatment
.
Eur Respir Rev
2021
;
30
:
200340
.

18

Khanna
D
,
Tashkin
DP
,
Denton
CP
et al.
Ongoing clinical trials and treatment options for patients with systemic sclerosis-associated interstitial lung disease
.
Rheumatology (Oxford)
2019
;
58
:
567
79
.

19

Distler
O
,
Gahlemann
M
,
Maher
TM.
Nintedanib for systemic sclerosis–associated interstitial lung disease
.
N Engl J Med
2019
;
381
:
1596
7
.

20

Khanna
D
,
Lin
CJF
,
Furst
DE
et al.
Tocilizumab in systemic sclerosis: a randomised, double-blind, placebo-controlled, phase 3 trial
.
Lancet Respir Med
2020
;
8
:
963
74
.

21

Hoffmann-Vold
AM
,
Maher
TM
,
Philpot
EE
et al.
The identification and management of interstitial lung disease in systemic sclerosis: evidence-based European consensus statements
.
Lancet Rheumatol
2020
;
2
:
e71
83
.

22

Kuwana
M
,
Allanore
Y
,
Denton
CP
et al.
Nintedanib in patients with systemic sclerosis-associated interstitial lung disease: subgroup analyses by autoantibody status and skin score
.
Arthritis Rheumatol
2022
;
74
:
518
26
.

23

Hoffmann-Vold
AM
,
Allanore
Y
,
Alves
M
et al.
Progressive interstitial lung disease in patients with systemic sclerosis-associated interstitial lung disease in the EUSTAR database
.
Ann Rheum Dis
2021
;
80
:
219
27
.

24

Distler
O
,
Assassi
S
,
Cottin
V
et al.
Predictors of progression in systemic sclerosis patients with interstitial lung disease
.
Eur Respir J
2020
;
55
:
1902026
.

25

Roofeh
D
,
Distler
O
,
Allanore
Y
et al.
Treatment of systemic sclerosis–associated interstitial lung disease: lessons from clinical trials
.
J Sclerod Related Disord
2020
;
5(Suppl 2
):
61
71
.

26

Meier
FM
,
Frommer
KW
,
Dinser
R
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
:
1355
60
.

27

Benvenuti
F
,
Zanatta
E
,
Avouac
J
et al.
The power of the EUSTAR cohort: key findings to date and implications for management of systemic sclerosis patients
.
Expert Rev Clin Immunol
2020
;
16
:
1065
74
.

28

van den Hoogen
F
,
Khanna
D
,
Fransen
J
et al.
2013 classification criteria for systemic sclerosis: an American College of Rheumatology/European League Against Rheumatism collaborative initiative
.
Ann Rheum Dis
2013
;
72
:
1747
55
.

29

Bütikofer
L
,
Varisco
PA
,
Distler
O
et al.
ACE inhibitors in SSc patients display a risk factor for scleroderma renal crisis-a EUSTAR analysis
.
Arthritis Res Ther
2020
;
22
:
59
.

30

Lang
RM
,
Badano
LP
,
Mor-Avi
V
et al.
Recommendations for cardiac chamber quantification by echocardiography in adults: an update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging
.
J Am Soc Echocardiogr
2015
;
28
:
1
39.e14
.

31

Pellegrino
R
,
Viegi
G
,
Brusasco
V
et al.
Interpretative strategies for lung function tests
.
Eur Respir J
2005
;
26
:
948
68
.

32

Lydersen
S.
Statistical review: frequently given comments
.
Ann Rheum Dis
2015
;
74
:
323
5
.

33

Nihtyanova
SI
,
Sari
A
,
Harvey
JC
et al.
Using autoantibodies and cutaneous subset to develop outcome-based disease classification in systemic sclerosis
.
Arthritis Rheumatol
2020
;
72
:
465
76
.

34

Boonstra
M
,
Ninaber
MK
,
Ajmone Marsan
N
et al.
Prognostic properties of anti-topoisomerase antibodies in patients identified by the ACR/EULAR 2013 systemic sclerosis criteria
.
Rheumatology (Oxford)
2019
:
58
:
730
2
.

35

Perera
A
,
Fertig
N
,
Lucas
M
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
:
2740
6
.

36

Boonstra
M
,
Bakker
JA
,
Grummels
A
et al.
Association of anti-topoisomerase I antibodies of the IgM isotype with disease progression in anti-topoisomerase I-positive systemic sclerosis
.
Arthritis Rheumatol
2020
;
72
:
1897
904
.

37

Iniesta Arandia
N
,
Espinosa
G
,
Tolosa Vilella
C
et al.
Serodiscordant patients with systemic sclerosis: when antibody does not correspond to skin involvement
.
Clin Exp Rheumatol
2020
;
38(Suppl 125
):
106
14
.

38

Lazzaroni
M-G
,
Cavazzana
I
,
Colombo
E
et al.
Malignancies in patients with anti-RNA polymerase III antibodies and systemic sclerosis: analysis of the EULAR scleroderma trials and research cohort and possible recommendations for screening
.
J Rheumatol
2017
;
44
:
639
47
.

39

Lazzaroni
MG
,
Marasco
E
,
Campochiaro
C
et al.
The clinical phenotype of Systemic Sclerosis patients with anti-PM/Scl antibodies: results from the EUSTAR cohort
.
Rheumatology (Oxford)
2021
;
60
:
5028
41
.

40

Heijnen
IA
,
Foocharoen
C
,
Bannert
B
et al.
Clinical significance of coexisting antitopoisomerase I and anticentromere antibodies in patients with systemic sclerosis: a EUSTAR group-based study
.
Clin Exp Rheumatol
2013
;
31
:
96
102
.

41

Suleman
Y
,
Clark
KEN
,
AR
C
et al.
Real-world experience of tocilizumab in systemic sclerosis: potential benefit on lung function for anti-topoisomerase-positive patients
.
Rheumatology (Oxford)
2021
;
60
:
3945
6
.

42

Bissell
L-A
,
Md Yusof
MY
,
Buch
MH.
Primary myocardial disease in scleroderma—a comprehensive review of the literature to inform the UK Systemic Sclerosis Study Group cardiac working group
.
Rheumatology
2016
;
56
:
882
95
.

43

Fox
BD
,
Shimony
A
,
Langleben
D
et al.
High prevalence of occult left heart disease in scleroderma-pulmonary hypertension
.
Eur Respir J
2013
;
42
:
1083
91
.

44

Volkmann
ER
,
Fischer
A.
Update on morbidity and mortality in systemic sclerosis–related interstitial lung disease
.
J Sclerod Related Disord
2021
;
6
:
11
20
.

45

Kafaja
S
,
Clements
PJ
,
Wilhalme
H
et al.
Reliability and minimal clinically important differences of forced vital capacity: results from the Scleroderma Lung Studies (SLS-I and SLS-II)
.
Am J Respir Crit Care Med
2018
;
197
:
644
52
.

46

Frantz
C
,
Huscher
D
,
Avouac
J
et al.
Outcomes of limited cutaneous systemic sclerosis patients: results on more than 12,000 patients from the EUSTAR database
.
Autoimmun Rev
2020
;
19
:
102452
.

47

Hoffmann-Vold
AM
,
Distler
O
,
Murray
B
et al.
Setting the international standard for longitudinal follow-up of patients with systemic sclerosis: a Delphi-based expert consensus on core clinical features
.
RMD Open
2019
;
5
:
e000826
.

48

Jaafar
S
,
Lescoat
A
,
Huang
S
et al.
Clinical characteristics, visceral involvement, and mortality in at-risk or early diffuse systemic sclerosis: a longitudinal analysis of an observational prospective multicenter US cohort
. Arthritis Res Ther
2021
;
23
:
170
.

49

Bruni
C
,
Buch
MH
,
De Luca
G
et al.
Primary systemic sclerosis heart involvement: a systematic literature review and preliminary data-driven, consensus-based WSF/HFA definition
.
J Sclerod Related Disord
2022
;
7
:
24
32
.

50

Hsu
VM
,
Moreyra
AE
,
Wilson
AC
et al.
Assessment of pulmonary arterial hypertension in patients with systemic sclerosis: comparison of noninvasive tests with results of right-heart catheterization
.
J Rheumatol
2008
;
35
:
458
65
.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://dbpia.nl.go.kr/pages/standard-publication-reuse-rights)

Supplementary data

Comments

0 Comments
Submit a comment
You have entered an invalid code
Thank you for submitting a comment on this article. Your comment will be reviewed and published at the journal's discretion. Please check for further notifications by email.