Demographic, laboratory, and disease-related data
A total of 188 participants, 73 patients with SSc (50 patients -68% - with limited form and 23 patients -32%- with diffuse form) and 115 age- and sex-matched controls, were included in this study. Demographic characteristics of patients and controls, as well as disease-related characteristics of the participants are shown in Table 1. There were no differences between patients and controls regarding body mass index and the presence of hypertension, current smoking, diabetes or obesity. Only waist circumference and current use of statins (34% vs. 11%, p=0.000) were found to be higher in patients with SSc. Many differences were found in the lipid profiles between patients and controls. In this sense, cholesterol, HDL and LDL-cholesterol, and apolipoprotein A1 were found to be lower in SSc patients. In contrast, triglycerides, lipoprotein A, the apoB:A1 ratio, and the atherogenic index revealed higher serum levels in SSc patients. The mean CEC of HDL was significantly lower in SSc patients compared to controls (17±11 vs. 8±3%, p=0.000) when the univariable analysis was performed.
Table 1. Demographics of SSc patients and controls vis-à-vis comorbidities, analytical and disease-related data.
|
|
|
Controls (n=115)
|
SSc patients (n=73)
|
p
|
Demographics
|
|
|
|
|
Female, n (%)
|
107 (93)
|
68 (93)
|
0.98
|
|
Age, years
|
60 ± 6
|
59 ± 1
|
0.40
|
|
BMI, mg/cm2
|
28 ± 5
|
29 ± 6
|
0.38
|
|
Waist circumference, cm
|
91 ± 14
|
96 ± 14
|
0.042
|
|
Systolic pressure, mmHg
|
131 ± 15
|
133 ± 19
|
0.54
|
|
Diastolic pressure, mmHg
|
84 ± 53
|
80 ± 13
|
0.53
|
Comorbidities, n(%)
|
|
|
|
|
Hypertension
|
40 (35)
|
27 (37)
|
0.76
|
|
Current smoking
|
19 (17)
|
14 (19)
|
0.61
|
|
Diabetes
|
7 (6)
|
8 (11)
|
0.22
|
|
BMI > 30
|
30 (26)
|
22 (30)
|
0.55
|
|
Statins
|
13 (11)
|
25 (34)
|
0.000
|
Analytical data
|
|
|
|
|
CRP, mg/dl
|
1.00 (1.00-3.00)
|
2.35 (1.15-4.29)
|
0.96
|
|
Cholesterol, mg/dl
|
223 ± 38
|
198 ± 39
|
0.000
|
|
Triglycerides, mg/dl
|
104 ± 52
|
153 ± 83
|
0.000
|
|
HDL-cholesterol, mg/dl
|
66 ± 18
|
49 ± 13
|
0.000
|
|
LDL-cholesterol, mg/dl
|
136 ± 35
|
111 ± 45
|
0.000
|
|
LDL:HDL cholesterol ratio
|
2.23 ± 0.84
|
2.50 ± 1.62
|
0.23
|
|
Non-HDL cholesterol, mg/dl
|
157 ± 37
|
149 ± 41
|
0.24
|
|
Lipoprotein A, mg/dl
|
16 (9.32)
|
23 (12.65)
|
0.020
|
|
Apolipoprotein A1, mg/dl
|
197 ± 34
|
168 ± 32
|
0.000
|
|
Apolipoprotein B, mg/dl
|
102 ± 22
|
100 ± 26
|
0.67
|
|
Apo B: Apo A ratio
|
0.54 ± 0.16
|
0.63 ± 0.30
|
0.024
|
|
Atherogenic index
|
3.58 ± 0.98
|
4.37 ± 1.74
|
0.003
|
|
Cholesterol efflux capacity, %
|
17 ± 11
|
8 ± 3
|
0.000
|
Systemic sclerosis-related data, n (%)
|
|
|
|
|
Disease duration, years
|
|
9 (4-15)
|
|
|
Modified Rodnan Skin Score, units
|
|
3 (1-6)
|
|
|
Raynaud phenomenon
|
|
62 (85)
|
|
|
Digital ulcers
|
|
8 (11)
|
|
|
Calcinosis
|
|
7 (10)
|
|
|
Arthritis
|
|
15 (21)
|
|
|
Gastric reflux
|
|
29 (40)
|
|
|
Pathological esophageal manometry
|
|
40 (55)
|
|
|
|
Interstitial lung disease
|
|
16 (22)
|
|
|
Pulmonary hypertension
|
|
16 (22)
|
|
|
Anti-centromere antibody positivity
|
|
47 (64)
|
|
|
Anti-Scl70 antibody
|
|
10 (14)
|
|
Treatments
|
|
Current prednisone use,
|
|
15 (21)
|
|
|
Prednisone, mg/day
|
|
|
5 (5-10)
|
|
|
DMARDs, n (%)
|
|
10 (14)
|
|
|
Azathioprine
|
|
4 (5)
|
|
|
Methotrexate,
|
|
2 (3)
|
|
|
Mycophenolate
|
|
2 (3)
|
|
|
Hydroxychloroquine,
|
|
2 (3)
|
|
Data represent patient numbers (%), means±SD or median (IQR) when data were not normally distributed.
BMI: body mass index; CRP: C-reactive protein; DMARD: disease-modifying antirheumatic drug; HDL: high-density lipoprotein; LDL: low-density lipoprotein.
Disease duration was 9 (IQR 4-15) years in SSc patients. Disease durations for the limited and diffuse types were 12 (IQR 6-17) and 8 (IQR 4-15) (p=0.13) years, respectively (see Supplementary Table 1). The mRSS score was 3 (IQR 1-6; minimum-maximum 0-27) for the entire study population (Table 1) and no significant differences were found (p=0.12) when this score was analyzed independently for the limited (3, IQR 0-6) and diffuse (12, IQR 6-17) forms (see Supplementary Table 1). The presence of digital ulcers and calcinosis was reported, respectively, in 11% and 10% of the patients. At the time the study was conducted, about one-fifth of patients (21%) were taking prednisone with a median dose of 5 [IQR 5-10] mg/day. Additionally, at the time of the study, 47 (64%) patients were found to be positive for anti-centromere, and 10 (14%) were positive for anti-Scl70. Disease-modifying anti-rheumatic drugs (DMARDs) use was reported in 14% of the patients, including 5% azathioprine, 3% hydroxychloroquine and 3% methotrexate. Other features related to the disease are shown in Table 1.
Multivariable analysis of the differences in lipid profiles between SSc patients and controls
Table 2 and Figure 1 shows the differences in lipid-related molecules between patients and controls. After adjusting for abdominal circumference and statins (Model 1 in Table 2), most of these molecules showed differences between the two populations. In this sense, total cholesterol, HDL and LDL-cholesterol, apolipoprotein A1, and CEC were downregulated in SSc patients to a statistically significant degree. In contrast, triglycerides, lipoprotein A, the apolipoprotein B:A1 ratio, and the atherogenic index showed higher levels in SSc patients.
Table 2. Multivariable analysis of the differences in lipid profiles between SSc patients and controls.
|
|
|
|
Univariate
|
|
Model
|
|
Model #2
|
|
|
Controls (n=115)
|
SSc patients (n=73)
|
|
model
|
|
#1
|
|
beta coef. (95% CI), p
|
|
Lipid profile
|
|
|
|
|
|
|
|
|
|
Cholesterol, mg/dl
|
223 ± 38
|
198 ± 39
|
|
0.000
|
|
0.001
|
|
-22 (-37- -7)
|
0.004
|
|
Triglycerides, mg/dl
|
104 ± 52
|
153 ± 83
|
|
0.000
|
|
0.001
|
|
24 (2-47)
|
0.033
|
|
HDL cholesterol, mg/dl
|
66 ± 18
|
49 ± 13
|
|
0.000
|
|
0.000
|
|
-1 (-5-3)
|
0.75
|
|
LDL cholesterol, mg/dl
|
136 ± 35
|
111 ± 45
|
|
0.000
|
|
0.001
|
|
|
|
|
LDL:HDL cholesterol ratio
|
2.23 ± 0.84
|
2.50 ± 1.62
|
|
0.23
|
|
0.098
|
|
|
|
|
Non-HDL cholesterol, mg/dl
|
157 ± 37
|
149 ± 41
|
|
0.24
|
|
0.24
|
|
|
|
|
Lipoprotein A, mg/dl
|
16 (9-32)
|
23 (12-65)
|
|
0.020
|
|
0.055
|
|
22 (2-43)
|
0.033
|
|
Apolipoprotein A1, mg/dl
|
197 ± 34
|
168 ± 32
|
|
0.000
|
|
0.000
|
|
-9 (-19-0)
|
0.062
|
|
Apolipoprotein B, mg/dl
|
102 ± 22
|
100 ± 26
|
|
0.67
|
|
0.52
|
|
|
|
|
Apo B: Apo A ratio
|
0.54 ± 0.16
|
0.63 ± 0.30
|
|
0.024
|
|
0.010
|
|
|
|
|
Atherogenic index
|
3.58 ± 0.98
|
4.37 ± 1.74
|
|
0.003
|
|
0.001
|
|
|
|
|
Cholesterol efflux capacity, %
|
17 ± 11
|
8 ± 3
|
|
0.000
|
|
0.000
|
|
-6 (-10- -2)
|
0.002
|
|
Data represent means ± standard deviation; Model #1: Adjusted for waist circumference and statins; Model #2: Adjusted for model #1 + rest of lipid molecules (with a p value < 0.20 in the univariate analysis) other than the one that is compared; Because collinearity LDL cholesterol, LDL:HDL ratio, non-HDL cholesterol, apoB:apoA, and atherogenic index were excluded from the multivariable analyses in Model 2, HDL: high-density lipoprotein; LDL: low-density lipoprotein.
Because lipid-related molecules are interrelated (they share metabolic pathways and distinguishing the effect of one from another is difficult), we performed a multivariable analysis adjusting for demographics and cardiovascular risk factors plus all the lipid-related molecules that were found to be different between patients and controls in Model 1 (Model 2 in Table 2). Because of collinearity, lipid molecules derived from a formula were excluded from the regression model (LDL-cholesterol, LDL:HDL ratio, non-HDL cholesterol, apoB:apoA1, and atherogenic index). Total cholesterol (beta coef. -22 [95%CI -37- -7], p=0.004), triglycerides (beta coef. 24 [95%CI 2-47], p=0.033), and lipoprotein A (beta coef. 22 [95%CI 2-43], p=0.033) maintained their differences between patients and controls. Interestingly, CEC conserved its downregulation in SSc patients after adjusting for other lipid profile-related molecules (beta coefficient -6 [95%CI -10- -2] %, p=0.002). Remarkably, HDL-cholesterol was not different between patients and controls after the multivariable analysis.
Association of demographics, lipid profile, and disease-related data with efflux capacity in SSc patients and controls
Demographics and lipid profiles were, in general, not associated with CEC in either patients or controls. In this sense, only abdominal circumference showed a positive association with CEC in patients (beta coefficient 0.09 [95%CI 0.03-0.14], p=0.002), though not in controls. Similarly, no traditional cardiovascular risk factors were linked to CEC in either population. Regarding lipid profiles, no correlations were identified between the standard lipid profile molecules and CEC. Remarkably, the use of statins was not related to CEC in patients or controls. Lastly, concerning SSc-related data, a negative association between mRSS and CEC was identified (beta coef. -0.19 [95%CI -0.34- -0.04]%, p=0.014) (Table 3). Disease duration was not related to CEC in either the diffuse or limited SSc types. Moreover, the use of calcium channel blockers, PDE5 inhibitors or endothelin receptor antagonists was not related to CEC (data not shown).
Table 3. Relation of demographics, lipid profile, and disease-related data with cholesterol efflux capability in SSc patients and controls.
|
|
|
Cholesterol efflux capacity beta coef. (CI95%), p
|
|
|
|
|
Controls (n=115)
|
Patients (n=73)
|
|
|
Age, years
|
-0.03 (-0.35-0.30), 0.88
|
-0.04 (-0.11-0.04), 0.32
|
|
|
Male
|
-3.34 (-11.34-4.66), 0.41
|
-1.67 (-4.81-1.47), 0.29
|
|
|
Body mass index, kg/m2
|
-0.18 (-0.59-0.24), 0.41
|
0.14 (0.00-0.28), 0.050
|
|
|
Abdominal circumference, cm
|
-0.13 (-0.27-0.02), 0.097
|
0.09 (0.03-0.14), 0.002
|
|
|
Systolic blood pressure, mmHg
|
-0.14 (-0.27-0.00), 0.057
|
0.01 (-0.03-0.06), 0.57
|
|
|
Diastolic blood pressure, mmHg
|
-0.02 (-0.06-0.02), 0.39
|
0.01 (-0.05-0.08), 0.67
|
|
Cardiovascular co-morbidities
|
|
|
|
Smoking
|
4.05 (-1.39-9.50), 0.14
|
1.50 (-0.47-3.47), 0.13
|
|
|
Diabetes
|
0.81 (-7.72-9.35), 0.85
|
1.63 (-0.87-4.12), 0.20
|
|
|
Hypertension
|
-2.75 (-7.00-1.51), 0.20
|
0.41 (-1.25-2.07), 0.62
|
|
|
Statins
|
-2.73 (-9.16-3.69), 0.40
|
0.35 (-1.32-2.02), 0.68
|
|
Analytical and lipid profiles
|
|
|
|
CRP, mg/dl
|
|
|
|
|
Cholesterol, mg/dl
|
0.04 (-0.02-0.09), 0.17
|
0.01 (-0.01-0.03), 0.35
|
|
|
Triglycerides, mg/dl
|
-0.02 (-0.06-0.03), 0.43
|
0.01 (-0.01-0.02), 0.36
|
|
|
HDL cholesterol, mg/dl
|
0.08 (-0.03-0.20), 0.16
|
-0.00 (-0.07-0.07), 0.97
|
|
|
LDL cholesterol, mg/dl
|
0.02 (-0.04-0.08), 0.43
|
0.01 (-0.01-0.03), 0.19
|
|
|
LDL: HDL cholesterol ratio
|
-0.31 (-2.78-2.17), 0.81
|
0.11 (-0.42-0.64), 0.69
|
|
|
Non-HDL cholesterol, mg/dl
|
0.01 (-0.04-0.07), 0.66
|
0.01 (-0.01-0.03), 0.37
|
|
|
Lipoprotein A, mg/dl
|
0.04 (-0.02-0.11), 0.16
|
0.01 (-0.01-0.02), 0.41
|
|
|
Apolipoprotein A1, mg/dl
|
0.03 (-0.03-0.09), 0.33
|
0.01 (-0.02-0.03), 0.62
|
|
|
Apolipoprotein B, mg/dl
|
-0.01 (-0.10-0.09), 0.91
|
0.01 (-0.02-0.04), 0.58
|
|
|
Apo B: Apo A ratio
|
-5.19 (-18.11-7.73), 0.43
|
-0.39 (-3.34-2.56), 0.79
|
|
|
Atherogenic index
|
-0.52 (-2.63-1.60), 0.63
|
0.01 (-0.50-0.52), 0.97
|
|
SS- and treatment-related data
|
|
|
|
log Disease duration, years
|
|
-0.51 (-1.47-0.45), 0.29
|
|
|
Modified Rodnan Skin Score, units
|
|
-0.19 (-0.34- -0.04), 0.014
|
|
|
Raynaud phenomenon
|
|
-0.54 (-2.89-1.80), 0.64
|
|
|
Digital ulcers
|
|
-2.50 (-4.99-0.01), 0.049
|
|
|
Calcinosis
|
|
-0.53 (-3.25-2.18), 0.70
|
|
|
Arthritis
|
|
0.28 (-2.26-1.71), 0.78
|
|
|
Gastric reflux
|
|
0.87 (-0.78-2.53), 0.30
|
|
|
Pathological esophageal manometry
|
|
0.15 (-1.78-2.08), 0.88
|
|
|
Interstitial lung disease
|
|
-1.34 (-3.24-0.55), 0.16
|
|
|
Pulmonary hypertension
|
|
-1.14 (-3.15-0.86), 0.26
|
|
|
Anti-centromere antibody positivity
|
|
-0.34 (-2.08-1.40), 0.70
|
|
|
Anti-Scl70 antibody positivity
|
|
1.80 (-0.49-4.09), 0.12
|
|
|
Current prednisone
|
|
-0.69 (-2.67-1.28), 0.49
|
|
|
log Prednisone
|
|
-2.22 (-6.83-2.40), 0.32
|
|
|
DMARDs
|
|
-0-00 (-2.33-2.33), 0.99
|
|
|
Azathioprine
|
|
-0.33 (-3.95-3.07), 0.80
|
|
|
Methotrexate, n (%)
|
|
0.96 (-3.93-5.86), 0.70
|
|
|
Hydroxychloroquine, n (%)
|
|
-0.85 (-5.74-4.05), 0.73
|
|
ANA: antinuclear antibodies; BMI: body mass index; CRP: C-reactive protein; DMARD: disease-modifying antirheumatic drug; HDL: high-density lipoprotein; LDL: low-density lipoprotein.
Modified Rodnan Skin Score's association with lipid profile molecules and CEC
mRSS was positively related to presence of hypertension, but negatively associated with CEC when the univariable correlations were assessed (Table 4). When the relation of mRSS to these lipid-related molecules was adjusted for traditional CV risk factors, similar results were found. Moreover, when the relation between mRSS and CEC was additionally adjusted for other lipid-related molecules, its significance was conserved (beta coef. -0.21 [95%CI -0.37- -0.05]%, p=0.011) (Figure 2). When CEC was analyzed by SSc forms, a significant difference was found between the limited and diffuse forms (7.98±3.33 vs 9.79±3.52, p=0.035). However, when this difference was analyzed by adjusting for demographic characteristics and cardiovascular comorbidities, no significant difference was found in CEC between the two forms of SSc (p=0.18) (see Supplementary Table 1).
Table 4. Modified Rodnan Skin Score association with lipid profile molecules and CV risk factors in SSc patients.
|
|
|
Modified Rodnan Skin Score
|
|
|
|
|
Pearson's r
|
p
|
|
Model #1
Beta coef. (95% CI), p
|
Model #2
Beta coef. (95% CI), p
|
|
Age, years
|
0.202
|
0.13
|
|
|
|
|
|
Male
|
-0.142
|
0.29
|
|
|
|
|
|
Body mass index, kg/m2
|
-0.028
|
0.84
|
|
|
|
|
|
Abdominal circumference, cm
|
-0.018
|
0.90
|
|
|
|
|
|
Systolic blood pressure, mmHg
|
0.016
|
0.91
|
|
|
|
|
|
Diastolic blood pressure, mmHg
|
-0.174
|
0.20
|
|
|
|
|
Cardiovascular co-morbidities
|
|
|
|
|
|
|
Smoking
|
0.009
|
0.94
|
|
|
|
|
|
Diabetes
|
0.182
|
0.17
|
|
|
|
|
|
Hypertension
|
0.294
|
0.025
|
|
|
|
|
|
Statins
|
0.109
|
0.41
|
|
|
|
|
Analytical and lipid profiles
|
|
|
|
|
|
|
Cholesterol, mg/dl
|
-0.176
|
0.25
|
|
-1 (-3-1), 0.25
|
-0.29 (-1.36-0.78), 0.59
|
|
|
Triglycerides, mg/dl
|
0.150
|
0.33
|
|
|
|
|
|
HDL cholesterol, mg/dl
|
0.031
|
0.84
|
|
|
|
|
|
LDL cholesterol, mg/dl
|
-0.237
|
0.10
|
|
-2 (-4-0), 0.10
|
|
|
|
LDL:HDL cholesterol ratio
|
-0.191
|
0.19
|
|
-0.03 (-0.08-0.02), 0.19
|
|
|
|
Non-HDL cholesterol, mg/dl
|
-0.182
|
0.23
|
|
-1 (-3-1), 0.23
|
|
|
|
Lipoprotein A, mg/dl
|
0.072
|
0.64
|
|
|
|
|
|
Apolipoprotein A1, mg/dl
|
0.168
|
0.27
|
|
0.61 (-0.80-2.03), 0.39
|
|
|
|
Apolipoprotein B, mg/dl
|
-0.220
|
0.15
|
|
-0.91 (-2.14-0.33), 0.15
|
-0.04 (-0.85-0.78), 0.93
|
|
|
Apo B:Apo A ratio
|
-0.280
|
0.062
|
|
-0.00 (-0.02-0.00), 0.062
|
|
|
|
Atherogenic index
|
-0.100
|
0.51
|
|
|
|
|
|
Cholesterol efflux capacity, %
|
-0.323
|
0.014
|
|
-0.19 (-0.34- -0.04), 0.014
|
-0.21 (-0.37- -0.05), 0.011
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
mRSS is considered the independent variable in the multivariable analyses. Model #1: Adjusted for gender, age and hypertension; Model #2: Adjusted for model #1 + rest of lipid molecules (with a p value < 0.20 in the univariate analysis) other than the one that is compared. Because collinearity parameters derived from a formula (LDL cholesterol, LDL:HDL ratio, non-HDL cholesterol, apoB:apoA, and atherogenic index) were excluded from the multivariable analyses in Model 2.
HDL: high-density lipoprotein; LDL: low-density lipoprotein.