Baseline clinical characteristics
One hundred and fourteen participants from the registry were included in the study. Two participants who underwent HSCT and four participants who did not undergo HSCT were excluded from the final analysis due to poor image quality for STE analysis. A majority (90%) of participants had diffuse SSc and nearly all (99%) were on immunosuppressant therapy prior to HSCT (Table 1). The median (IQR) time between diagnosis of SSc and HSCT was 2.7 (1.5-6.4) years. Compared to HSCT recipients, participants in the non-HSCT cohort were older (59±6 versus 51±11 years, p<0.01), had higher systolic blood pressure, and lower frequency of seropositivity to antinuclear antibody and anti-topoisomerase. Patients in the non-HSCT group were more likely to have a history of coronary artery disease, diabetes mellitus, hyperlipidemia, and were taking concurrent cardiovascular medications, including ACEi/ARB, diuretics, beta-blockers and phosphodiesterase inhibitors.
Among patients who received HSCT, the mean pulmonary artery pressure (mPAP) measured prior to HSCT was 18.4±4.4 mmHg, and 27 (31%) participants had a mPAP greater than 20 mmHg (Supplemental Table 2). Four patients had a mPAP greater than 25 mmHg but were deemed to have post-capillary pulmonary hypertension with elevated left-sided filling pressures and normal pulmonary vascular resistance, and thus were considered appropriate for HSCT.
Table 1. Baseline clinical characteristics of 108 study participants by HSCT status.
Clinical characteristics, No. (%)
|
HSCT
(N=88)
|
Non-HSCT
(N=20)
|
P-value1
|
Age, years, mean±SD
|
51±11
|
59±6
|
<0.01
|
Female gender
|
66 (75)
|
19 (95)
|
0.05
|
Race
|
|
|
0.89
|
Caucasian
|
67 (76)
|
14 (70)
|
|
African American
|
10 (11)
|
3 (15)
|
|
Other
|
11 (13)
|
3 (15)
|
|
Body mass index, kg/m2, mean±SD
|
24±5
|
25±5
|
0.41
|
Heart rate, beats/min, mean±SD
|
84±14
|
86±15
|
0.56
|
Systolic blood pressure, mmHg, mean±SD
|
110±15
|
118±14
|
0.03
|
Diastolic blood pressure, mmHg, mean±SD
|
69±10
|
70±11
|
0.64
|
Hemoglobin, g/dL, mean±SD
|
12.0±1.7
|
12.5±1.5
|
0.20
|
BNP, pg/mL, median (IQR)
|
35 (20-70)
|
42 (22-97)
|
0.27
|
mRSS, median (IQR)
|
21 (13-34)
|
14 (10-25)
|
0.05
|
SSc characteristics and Seropositivity
|
|
|
|
Diffuse cutaneous SSc
|
79 (90)
|
17 (85)
|
0.54
|
Antinuclear antibody
|
86 (98)
|
16 (80)
|
<0.01
|
Anti-RNA polymerase III
|
30 (34)
|
9 (45)
|
0.36
|
Anti-topoisomerase I
|
32 (36)
|
2 (10)
|
0.02
|
Anti-centromere
|
5 (6)
|
0 (0)
|
0.27
|
Comorbidities/Cardiovascular risk factors
|
|
|
|
Coronary artery disease
|
2 (2)
|
3 (15)
|
0.02
|
Systemic arterial hypertension
|
9 (10)
|
5 (25)
|
0.08
|
Diabetes mellitus
|
1 (1)
|
2 (10)
|
0.03
|
Hyperlipidemia
|
7 (8)
|
8 (40)
|
<0.01
|
Smoking history
|
18 (20)
|
6 (30)
|
0.35
|
Medications
|
|
|
|
Calcium channel blockers
|
26 (30)
|
10 (50)
|
0.08
|
ACEi/ARB
|
18 (20)
|
19 (95)
|
<0.01
|
Diuretics
|
11 (13)
|
8 (40)
|
<0.01
|
Beta-blockers
|
3 (3)
|
5 (25)
|
<0.01
|
Phosphodiesterase inhibitors
|
13 (15)
|
8 (40)
|
0.01
|
Prostacyclin receptor agonists
|
1 (1)
|
0 (0)
|
0.63
|
Endothelin receptor antagonists
|
2 (2)
|
2 (10)
|
0.10
|
Mycophenolate mofetil
|
70 (80)
|
17 (85)
|
0.58
|
Prednisone
|
59 (67)
|
9 (45)
|
0.07
|
Methotrexate
|
47 (53)
|
3 (15)
|
<0.01
|
Cyclophosphamide
|
25 (28)
|
2 (10)
|
0.09
|
Hydroxychloroquine
|
16 (18)
|
3 (15)
|
0.74
|
Biologics or other
|
16 (18)
|
1 (5)
|
0.14
|
Intravenous Immunoglobulin
|
11 (13)
|
2 (10)
|
0.76
|
Azathioprine
|
11 (13)
|
2 (10)
|
0.76
|
Rituximab
|
14 (16)
|
1 (5)
|
0.20
|
D-penicillamine
|
1 (1)
|
0 (0)
|
0.63
|
Abbreviations: ACEi: angiotensin-converting enzyme inhibitors; ARB: aldosterone receptor antagonists; BNP: brain natriuretic peptide; HSCT: hematopoietic stem cell transplant; mRSS: modified Rodnan skin score; RNA: ribonucleic acid; SSc: systemic sclerosis.
1Two-sample t-test, Wilcoxon sum rank test, or chi-squared test as appropriate.
Baseline echocardiographic characteristics
The pre-HSCT echocardiogram was performed a median (IQR) of 2 (1-4) months before HSCT. The median (IQR) time between diagnosis of SSc and baseline echocardiogram was 2.4 (1.3-6) years for the HSCT group, and 0.9 (0.2-1.6) years for the non-HSCT group. At baseline, all mean values of conventional echocardiography measures were within the normal range for both the HSCT and non-HSCT groups, and there were no significant differences in these measures between the groups (Table 2). In the HSCT group, 56% of these patients had abnormal LV GLS, 22% had abnormal LV GCS, 48% had abnormal RV FWS, and 70% had abnormal LA reservoir strain. In the non-HSCT group, 55% had abnormal LV GLS, 30% had abnormal LV GCS, 30% had abnormal RV FWS, and 35% had abnormal LA reservoir strain.
Follow-up echocardiographic findings
The follow-up echocardiographic evaluation was performed a median of 12 (6-14) months after HSCT. Among those patients who did not receive HSCT, the follow-up echocardiographic evaluation was performed a median of 12 (10-13) months after baseline. At follow-up, significant differences were noted in several conventional echocardiographic measures in both groups (Table 2). However, the absolute difference was small, and all follow-up measures remained in the normal range.
Among strain parameters, both RVGLS and RVFWS significantly improved post-HSCT. Regionally, the improvement was observed in the mid (20.4±9.5% vs 23.7±8.0%, p=0.04) and apical segments (15.3±8.6% vs 20.9±9.0%, p<0.01) of the RV free wall, but not in the basal segment. LA reservoir strain also improved (35.9±8.7% vs 47.8±11.4%, p<0.01) and LA stiffness index decreased from 0.24 to 0.18 (p<0.01). No significant changes were observed in LVGLS or LVGCS after HSCT. In the patients who did not receive HSCT, no significant changes in strain parameters were observed.
Table 2. Echocardiographic and clinical measures of 108 study participants by HSCT status.
Echocardiography variables, mean±SD
|
HSCT (N=88)
|
Non-HSCT (N=20)
|
|
Baseline
|
Follow-Up
|
P-value1
|
Baseline
|
Follow-Up
|
P-value2
|
P-value3
|
Left Ventricle
|
|
|
|
|
|
|
|
End-diastolic diameter index, cm/m2
|
2.51±0.30
|
2.47±0.32
|
0.14
|
2.46±0.25
|
2.54±0.32
|
0.28
|
0.52
|
IVS thickness, cm
|
0.89±0.18
|
0.90±0.18
|
0.68
|
0.96±0.25
|
0.93±0.19
|
0.62
|
0.18
|
Posterior wall thickness, cm
|
0.92±0.16
|
0.93±0.17
|
0.51
|
0.92±0.22
|
0.93±0.18
|
0.87
|
0.97
|
Mass index, g/m2
|
74.7±22.0
|
74.1±21.7
|
0.82
|
75.5±25.3
|
77.3±24.2
|
0.60
|
0.88
|
Ejection Fraction, %
|
61.7±5.3
|
61.0±6.9
|
0.37
|
62.6±5.2
|
59.9±5.5
|
0.02
|
0.49
|
GLS, %
|
18.7±4.4
|
19.0±3.4
|
0.61
|
19.8±3.8
|
19.2±3.8
|
0.34
|
0.31
|
GCS, %
|
26.4±7.1
|
25.7±6.8
|
0.49
|
24.5±7.7
|
25.2±6.1
|
0.54
|
0.29
|
Left Atrium
|
|
|
|
|
|
|
|
Volume index, ml/m2
|
24.7±8.3
|
24.2±7.2
|
0.64
|
24.8±6.6
|
25.2±10.9
|
0.88
|
0.95
|
Septal mitral annular e’ velocity, cm/sec
|
10.1±3.0
|
9.3±2.8
|
0.01
|
9.2±2.0
|
8.4±2.2
|
0.25
|
0.22
|
Lateral mitral annular e’ velocity, cm/sec
|
12.2±3.4
|
11.5±3.4
|
0.03
|
11.8±2.5
|
10.4±2.3
|
0.09
|
0.64
|
Transmitral Doppler E/A ratio
|
1.4±0.5
|
1.2±0.4
|
0.01
|
1.3±0.3
|
1.4±0.4
|
0.66
|
0.57
|
E/e’
|
8.0±2.4
|
8.3±2.6
|
0.30
|
8.2±1.6
|
10.3±3.1
|
0.04
|
0.73
|
Stiffness index
|
0.24±0.12
|
0.18±0.08
|
<0.01
|
0.21±0.07
|
0.29±0.19
|
0.06
|
0.21
|
Reservoir strain, %
|
35.9±8.7
|
47.8±11.4
|
<0.01
|
43.2±12.5
|
40.9±11.3
|
0.31
|
<0.01
|
Right Atrium
|
|
|
|
|
|
|
|
Right atrial area, cm2
|
13.3±3.8
|
13.4±4.5
|
0.78
|
12.6±4.1
|
13.6±3.8
|
0.20
|
0.45
|
Right Ventricle
|
|
|
|
|
|
|
|
Basal diameter, cm
|
3.6±1.1
|
3.4±0.6
|
0.13
|
3.5±0.7
|
3.7±0.7
|
0.36
|
0.78
|
End-diastolic area, cm2
|
16.7±5.5
|
17.5±4.9
|
0.01
|
18.1±5.6
|
19.2±4.4
|
0.31
|
0.34
|
End-systolic area, cm2
|
9.6±3.9
|
9.7±3.7
|
0.21
|
10.9±3.4
|
11.7±3.6
|
0.26
|
0.16
|
Fractional area change, %
|
43.0±11.8
|
44.8±11.0
|
0.25
|
39.2±8.4
|
39.9±9.7
|
0.79
|
0.17
|
TAPSE, cm
|
2.2±0.4
|
2.1±0.4
|
<0.01
|
2.4±0.6
|
2.3±0.6
|
0.26
|
0.13
|
Tricuspid annular S’ velocity, cm/sec
|
12.8±1.9
|
12.0±2.2
|
<0.01
|
12.2±2.4
|
13.0±2.9
|
0.91
|
0.27
|
TR velocity, cm/sec
|
2.4±0.4
|
2.3± 0.6
|
0.31
|
2.6±0.5
|
2.9±0.8
|
0.22
|
0.05
|
Pericardial effusion, No. (%)
|
12 (14)
|
10 (11)
|
0.62
|
2 (10)
|
4 (20)
|
0.63
|
0.66
|
GLS, %
|
18.1±3.9
|
20.0±4.5
|
<0.01
|
19.9±4.5
|
19.5±4.9
|
0.73
|
0.07
|
Free Wall Strain, %
|
20.7±5.3
|
23.2±5.6
|
<0.01
|
23.2±5.4
|
23.1±5.3
|
0.97
|
0.07
|
Free Wall Strain Segments
|
|
|
|
|
|
|
|
Basal, %
|
27.9±9.4
|
28.0 ± 7.9
|
0.94
|
29.0±7.7
|
28.5±9.7
|
0.76
|
0.63
|
Mid, %
|
20.4±9.5
|
23.7 ± 8.0
|
0.04
|
23.0±7.0
|
22.2±5.3
|
0.62
|
0.27
|
Apical, %
|
15.3±8.6
|
20.9 ± 9.0
|
<0.01
|
15.7±8.7
|
18.0±10.0
|
0.42
|
0.89
|
Clinical outcome measures
|
|
|
|
|
|
|
|
mRSS, median (IQR)
|
21 (13-34)
|
9 (4-20)
|
<0.01
|
14 (10-25)
|
10 (8-23)
|
0.01
|
0.01
|
Abbreviations: A: late diastolic mitral inflow velocity; E: early diastolic mitral inflow velocity; e’: early diastolic mitral annulus velocity; FWS: free wall strain; GCS: global circumferential strain; GLS: global longitudinal strain; IVS: interventricular septum; mRSS: modified Rodnan skin score; TAPSE: tricuspid annular plane systolic excursion; TR: tricuspid regurgitation.
1Paired t-test, McNemar’s test or Wilcoxon sum rank test as appropriate comparing baseline and follow-up measures of the HSCT group.
2Paired t-test, McNemar’s test or Wilcoxon sum rank test as appropriate comparing baseline and follow-up measures of the non-HSCT group.
3Two-sample t-test or chi-squared test as appropriate comparing baseline measures of the HSCT and non-HSCT group.
Comparison of echocardiographic changes between HSCT and non-HSCT patients
Among participants who received HSCT, there was a 9.70 (95% CI [3.80, 15.61], p<0.01) unit increase in LA reservoir strain at follow-up compared to those who did not receive HSCT after controlling for baseline measures. After additional adjustment for age, sex, race, and comorbidities, HSCT therapy remained significantly associated with an improvement in LA reservoir strain at follow-up. There were no significant between-group differences observed for changes in RV and LV strain (Table 3 and Figure 2).
Strain parameters and clinical outcomes
Among the HSCT recipients, mRSS decreased, reflecting an improvement in symptom burden (21 [13-34] at baseline to 9 [4-20] follow up, p<0.01). Among participants who did not undergo HSCT and were only treated with conventional immunosuppressant therapies, mRSS was lower at baseline compared to the HSCT recipients, but still improved at follow up (14 [10-25] at baseline to 10 [8-23] follow up, p=0.01, Table 2). However, the between-group analysis demonstrated that there was an improvement of a 6 (95% CI [2.36, 9.61], p<0.01) unit reduction in mRSS at follow up in participants who received HSCT compared to those who did not after adjustment for baseline mRSS (Table 3). There was no significant correlation between change in strain parameters and change in mRSS (Supplemental Table 3).
Table 3. Association between treatment status and change in mRSS or strain parameters.
|
Baseline Strain
|
HSCT Treatment
|
β coefficient (95% CI)
|
P-value
|
β coefficient (95% CI)
|
P-value
|
ANCOVA1
|
|
|
|
|
mRSS
|
-0.34 (-0.46, -0.22)
|
<0.01
|
-5.98 (-9.61, -2.36)
|
<0.01
|
LVGLS
|
-0.61 (-0.76, -0.47)
|
<0.01
|
0.25 (-1.31, 1.80)
|
0.75
|
LVGCS
|
-0.64 (-0.80, -0.47)
|
<0.01
|
-0.15 (-3.23, 2.93)
|
0.92
|
RVGLS
|
-0.68 (-0.90, -0.47)
|
<0.01
|
1.17 (-1.05, 3.39)
|
0.30
|
RVFWS
|
-0.72 (-0.92, -0.52)
|
<0.01
|
0.81 (-1.91, 3.52)
|
0.56
|
LA reservoir strain
|
-0.67 (-0.91, -0.44)
|
<0.01
|
9.70 (3.80, 15.61)
|
<0.01
|
Demographic Adjusted2
|
|
|
|
|
mRSS
|
-0.33 (-0.46, -0.20)
|
<0.01
|
-6.38 (-11.16, -1.60)
|
0.01
|
LVGLS
|
-0.63 (-0.78, -0.47)
|
<0.01
|
0.52 (-1.33, 2.37)
|
0.58
|
LVGCS
|
-0.65 (-0.83, -0.48)
|
<0.01
|
1.44 (-2.15, 5.02)
|
0.43
|
RVGLS
|
-0.67 (-0.91, -0.44)
|
<0.01
|
1.17 (-1.43, 3.77)
|
0.37
|
RVFWS
|
-0.70 (-0.91, -0.49)
|
<0.01
|
1.22 (-1.97, 4.40)
|
0.45
|
LA reservoir strain
|
-0.66 (-0.90, -0.42)
|
<0.01
|
8.03 (1.15, 14.90)
|
0.02
|
Abbreviations: FWS: free wall strain; GCS: global circumferential strain; GLS: global longitudinal strain; LA: left atrium; LV: left ventricle; mRSS: modified Rodnan skin score; RV: right ventricle.
1ANCOVA models evaluate the effect of HSCT on the change in strain adjusted for baseline measures.
2Demographic adjusted models were additionally adjusted for age, gender, race and comorbidities.