Abstract

Aims

Epicardial adipose tissue (EAT) has emerged as a mediator between systemic inflammatory disorders and cardiovascular disease, and may therefore play a role in the pathophysiology of cardiac involvement in systemic sclerosis (SSc). The aim of this study was to assess the correlation between EAT and left ventricular (LV) function, and to determine the prognostic value of EAT in patients with SSc.

Methods and results

Consecutive patients with SSc who underwent non-contrast thorax computed tomography and echocardiography were included. EAT mass was quantified using dedicated software. The study endpoint was all-cause mortality. A total of 230 SSc patients (age 53 ± 15 years, 14% male) were included. The median value of EAT mass was 67 g (interquartile range: 45–101 g). Patients with increased EAT mass (≥67 g) showed more impaired LV diastolic function as compared with patients with less EAT mass (<67 g), and even after adjusting for age and comorbidities, EAT mass was independently associated with LV diastolic function parameters. During a median follow-up of 8 years, 42 deaths occurred. Kaplan–Meier analysis showed that patients with increased EAT mass had higher all-cause mortality rate as compared with patients with less EAT mass (29% vs. 7%; P < 0.001). In the multivariable analysis, EAT was independently associated with all-cause mortality after adjusting for important covariates (HR: 1.006; 95% CI: 1.001–1.010).

Conclusion

In patients with SSc, EAT is independently associated with LV diastolic dysfunction and higher mortality rate.

Introduction

Systemic sclerosis (SSc) is an autoimmune connective tissue disease characterized by microvascular injury and progressive fibrosis of skin and visceral organs.1 The heart is frequently involved in SSc, with cardiovascular abnormalities being one of the major causes of mortality in these patients.2 Cardiac manifestations typically include myocarditis, heart failure with systolic and/or diastolic dysfunction, pulmonary arterial hypertension, pericardial involvement, and arrhythmias.3 Particularly, previous studies indicated that left ventricular (LV) diastolic dysfunction is relatively common in patients with SSc, ranging from 18% to 62%,4,5 and is an important predictor of mortality.6 Potential substrates leading to LV diastolic dysfunction are the presence of diffuse and focal myocardial fibrosis, possibly secondary to inflammation or to ischaemia from microcirculation impairment.7–9 Also, accumulation of epicardial adipose tissue (EAT) has been associated with local inflammation that may involve the adjacent cardiac tissues, causing microvascular dysfunction and fibrosis.10–12 When the process adjoins the LV, it impairs myocardial distensibility, increasing LV stiffness and filling pressures and leading to LV diastolic dysfunction.8,13 There is already growing evidence that EAT is related to cardiovascular risk factors, coronary artery disease, atrial fibrillation, and heart failure with preserved ejection fraction (HFpEF).14–16 We therefore hypothesized that EAT extent, as assessed by computed tomography (CT), may be related to cardiac involvement in patients with SSc and particularly to LV diastolic dysfunction. Aims of the present study were two-fold: (i) to explore the association between EAT and cardiac dysfunction in patients with SSc; and (ii) to determine the prognostic value of EAT in these patients.

Methods

Study population

Consecutive SSc patients referred to Leiden University Medical Center (Leiden, The Netherlands) and enrolled in the Leiden Comprehensive Care in Systemic Sclerosis (CCISS) Cohort (REU 043/SH/sh) were included in this study from April 2009 to May 2015.17,18 Patients were diagnosed according to the American College of Rheumatology/the European League Against Rheumatism classification criteria for SSc.19 Multidisciplinary team care and organ screening including thorax CT and echocardiography were performed to assess interstitial lung disease (ILD) and cardiac involvement for all patients, who gave written informed consent for this assessment as according to the approval of the Leiden University Medical Center Medic Ethical Committee. However, only patients who underwent echocardiography within one month from the initial thorax CT were enrolled. In addition, patients with missing baseline clinical data and uninterpretable thorax CT images were excluded. Also, patients with significant valvular heart disease and pericardial involvement were excluded to avoid interference with the EAT measures.

Disease-related characteristics included among others the modified Rodnan skin score (mRSS), SSc subtype (limited or diffuse cutaneous), the time since diagnosis, and the time since Raynaud. Complete laboratory testing was performed, including renal function, N-terminal pro-brain natriuretic peptide (NT-proBNP), C-reactive protein (CRP), and creatine phosphokinase (CPK). Pulmonary function was assessed by performing spirometry tests according to the American Thoracic Society/European Respiratory Society recommendations,20,21 and included maximum O2 uptake during exercise (VO2max) and diffusion capacity of the lung for carbon monoxide (DLCO). Pulmonary hypertension (PAH) was diagnosed according to 2015 ESC/ERS guidelines for the diagnosis and treatment of pulmonary hypertension.22

Echocardiography

Transthoracic echocardiography was performed using a commercially available ultrasound system equipped with 3.5 MHz transducers (Vivid 7 and E9, GE-Vingmed, Horten, Norway). Standard M-mode and two-dimensional, colour, pulsed, and continuous wave Doppler images were acquired. The images were digitally stored and retrospectively measured (EchoPAC Version 203.0.1, GE Medical Systems, Horten, Norway). LV volumes were measured on the apical four- and two-chamber views and LV ejection fraction (LVEF) was calculated using the Simpson biplane method.23 The left atrial volume was measured on the apical four- and two-chamber views according to the method of discs and indexed to body surface area [left atrial volume index (LAVI)]. LV internal diameters, interventricular septal thickness (IVST), and posterior wall thickness (PWT) were measured at end diastole (and end systole for LV diameter). LV mass was calculated based on the recommended formula and indexed to body surface area [left ventricular mass index (LVMI)].23 The peak velocity of the tricuspid regurgitation (TR) jet was derived from continuous wave Doppler. Inspiratory collapsibility and the inferior vena cava diameter were measured to determine right atrial pressure and calculate pulmonary pressure. Pulsed wave Doppler of mitral inflow was used to measure the peak early diastolic (E) wave and late diastolic (A) wave. The ratio of peak early diastolic wave-to-late diastolic wave was calculated. The septal and lateral peak early diastolic mitral annular velocities (Eʹ) were measured using tissue Doppler imaging in the apical four-chamber view. Subsequently, the ratio between E and Eʹ was calculated to assess LV filling pressures.

CT acquisition

Non-contrast thorax CT scans were acquired using a 64-detector row helical scanner (Aquilion 64, Toshiba Medical Systems, Toshiba Medical Systems, Otawara, Japan) or a 320-detector row volumetric scanner (Aquilion ONE, Toshiba Medical Systems, Otawara, Japan). The CT datasets were reconstructed at a slice thickness of either 2.5 or 3.0 mm in hybrid iterative reconstruction technique (iDose5, Philips Healthcare, Best, The Netherlands).

Quantification of EAT mass

EAT was measured in gram (mass) using the MASS software (Leiden University Medical Center, Leiden, The Netherlands) as described before.24 Briefly, a cardiac cross-sectional view of thorax non-contrast CT was acquired from the reconstructed four-chamber views with 2 mm slice thickness. The limits of the heart were defined as the pulmonary artery bifurcation (superior limit) to the posterior descending coronary artery (inferior limit). The pericardium was manually traced from the pulmonary artery bifurcation until the posterior descending coronary artery. The EAT was automatically identified by the software as tissue with HU between −195 and −45 within the region of pericardium24 (Figure 1A). To define the intra- and inter-observer variability, the EAT mass measurements were repeated for 10 randomly selected patients > 4 weeks apart by the same observer (X.W. blinded to the echocardiographic results) on the same thorax CT images and by a second independent observer (R.M.).

EAT mass assessment. (A) Example of quantification of EAT mass. Manual tracing was performed from the pulmonary artery bifurcation (superior limit) to the PDA (inferior limit), and the adipose tissue within the pericardium was automatically recognized by the software as tissue with Hounsfield Units (HU) between −195 and −45. (B) Histograms of frequencies of EAT mass. Median value of EAT mass in study cohort was 67 g. EAT, epicardial adipose tissue; PDA, posterior descending coronary artery.
Figure 1

EAT mass assessment. (A) Example of quantification of EAT mass. Manual tracing was performed from the pulmonary artery bifurcation (superior limit) to the PDA (inferior limit), and the adipose tissue within the pericardium was automatically recognized by the software as tissue with Hounsfield Units (HU) between −195 and −45. (B) Histograms of frequencies of EAT mass. Median value of EAT mass in study cohort was 67 g. EAT, epicardial adipose tissue; PDA, posterior descending coronary artery.

Follow-up and study endpoint

The study endpoint was all-cause mortality. Medical records review and survival data were obtained from the hospital information systems (EPD-Vision and EZIS; Leiden University Medical Center, Leiden, The Netherlands) or the Social Security Death Index. Follow-up began from the date that CT was performed.

Statistical analysis

Normally distributed continuous variables were presented as mean ± standard deviation; and non-normally distributed continuous variables were presented as median values and interquartile range (IQR). The Student’s t-test was used for comparison of normally distributed continuous variables, while the Mann–Whitney U test was used for non-normally distributed continuous variables. Categorical data were presented as frequencies with percentages and were compared using the χ2 test. For the reproducibility of EAT measurements, the intra-class correlation coefficient (ICC) was calculated to assess the inter- and intra-observer variability. Correlation between EAT and echocardiographic cardiac functional parameters was evaluated in multivariable linear regression analysis. The cumulative survival rate was calculated using Kaplan–Meier analysis, and comparison between groups was performed using the log-rank test. Cox regression analyses were used to investigate the association between clinical, echocardiographic, and CT variables and all-cause mortality. Variables with a P-value < 0.05 at the univariable analysis were included in the Cox proportional hazards backward stepwise model. The hazard ratio (HR) and 95% confidence interval were reported. The incremental value of EAT over baseline clinical variables for the risk of all-cause mortality was assessed using the likelihood ratio test. All statistical tests were two-sided, and a P-value < 0.05 was considered statistically significant. Statistical analysis was performed using SPSS version 25.0 (IBM Corporation, Armonk, NY, USA) and R version 4.1.1 (R Foundation for Statistical Computing, Vienna, Austria).

Results

Patient characteristics

A total of 230 patients [mean age 53 ± 15 years and 33 (14%) males] with SSc were enrolled. In the overall population, the median value of EAT mass was 67 g (IQR: 45 to 101 g) (Figure 1B), by which the subjects were therefore dichotomized. Other baseline characteristics of the study cohort are summarized in Table 1. Patients with increased EAT (≥67 g) were older as compared with patients with less EAT (<67 g), had more frequently PAH, and had also significantly higher prevalence of comorbidities including diabetes and hypertension. Moreover, patients with increased EAT had significantly higher values of CRP and NT-proBNP, lower mean estimated glomerular filtration rate (eGFR), and more impaired DLCO as compared with patients with less EAT. Finally, cardiovascular medications were more frequently used in patients with increased EAT, including angiotensin-converting enzyme inhibitors or angiotensin II receptor blockers, calcium channel blockers, beta-blockers, and diuretics. These results suggest that patients with higher EAT had also more progressed disease stage (including cardiac, renal, and pulmonary involvement) and more comorbidities, with a worse cardiovascular burden (diabetes and hypertension).

Table 1

Baseline clinical characteristics of the total SSc population and divided based on the median value of EAT mass

Total populationEAT ≥ 67 gEAT < 67 gP-value
(n = 230)(n = 115)(n = 115)
Clinical variables
 Age, y53 ± 1560 ± 1145 ± 14<0.001
 Male, n (%)33(14.3)24(20.9)9(7.8)0.005
 BMI, kg/m224.6 ± 4.326.0 ± 4.423.1 ± 3.6<0.001
 Diffuse SSc, n (%)51(23.1)24(21.6)27(24.5)0.458
 HTN, n (%)58(25.2)39(33.9)19(16.5)0.002
 DM, n (%)11(4.8)9(7.8)2(1.7)0.031
 CAD, n (%)8(3.5)6(5.2)1(0.8)0.150
 AF, n (%)7(3.2)6(5.5)1(0.9)0.055
 History of smoking, n (%)110(47.8)61(53.0)49(42.6)0.113
 mRSS5 ± 86 ± 94 ± 50.124
 ILD, n (%)91(41.2)52(46.8)39(35.5)0.085
 PAH, n (%)8(3.6)7(6.3)1(0.9)0.032
 Time since Raynaud, y13 ± 1415 ± 1511 ± 120.013
 Time since diagnosis, y5 ± 76 ± 84 ± 60.014
Laboratory findings
 eGFR, mL/min/1.73 m288(72–103)82(69–96)95(79–106)<0.001
 NT-proBNP, ng/L95.8(53.6–170.5)107.6(59.2–252.4)79.0(49.1–137.6)0.005
 CRP, mg/L0(0–5)0(0–6)0(0–0)0.002
 CPK, U/L83(61–118)93(60–134)79(61–109)0.239
 ANA-positive, n (%)211(95.5)105(94.6)106(96.4)0.423
 ACA-positive, n (%)87(39.4)44(39.6)43(39.1)0.933
 ATA-positive, n (%)54(24.4)23(20.7)31(28.3)0.197
 DLCO%64.7 ± 19.159.5 ± 18.570.0 ± 18.2<0.001
 VO2max%87.0 ± 23.785.1 ± 22.488.9 ± 24.90.262
Immunosuppressive medication, n (%)
 Corticosteroids31(13.5)18(15.7)13(11.3)0.334
 Cyclophosphamide3(1.3)2(1.7)1(0.9)0.561
 Methotrexate29(12.6)18(15.7)11(9.6)0.164
 Azathioprine11(4.8)9(7.8)2(1.7)0.031
Cardiovascular medications, n (%)
 ACEI/ARB85(37.0)57(49.6)28(24.3)<0.001
 CCB108(47.0)64(55.7)44(38.3)0.008
 Beta-blockers21(9.1)17(14.8)4(3.5)0.003
 Diuretics30(13.0)22(19.1)8(7.0)0.006
Total populationEAT ≥ 67 gEAT < 67 gP-value
(n = 230)(n = 115)(n = 115)
Clinical variables
 Age, y53 ± 1560 ± 1145 ± 14<0.001
 Male, n (%)33(14.3)24(20.9)9(7.8)0.005
 BMI, kg/m224.6 ± 4.326.0 ± 4.423.1 ± 3.6<0.001
 Diffuse SSc, n (%)51(23.1)24(21.6)27(24.5)0.458
 HTN, n (%)58(25.2)39(33.9)19(16.5)0.002
 DM, n (%)11(4.8)9(7.8)2(1.7)0.031
 CAD, n (%)8(3.5)6(5.2)1(0.8)0.150
 AF, n (%)7(3.2)6(5.5)1(0.9)0.055
 History of smoking, n (%)110(47.8)61(53.0)49(42.6)0.113
 mRSS5 ± 86 ± 94 ± 50.124
 ILD, n (%)91(41.2)52(46.8)39(35.5)0.085
 PAH, n (%)8(3.6)7(6.3)1(0.9)0.032
 Time since Raynaud, y13 ± 1415 ± 1511 ± 120.013
 Time since diagnosis, y5 ± 76 ± 84 ± 60.014
Laboratory findings
 eGFR, mL/min/1.73 m288(72–103)82(69–96)95(79–106)<0.001
 NT-proBNP, ng/L95.8(53.6–170.5)107.6(59.2–252.4)79.0(49.1–137.6)0.005
 CRP, mg/L0(0–5)0(0–6)0(0–0)0.002
 CPK, U/L83(61–118)93(60–134)79(61–109)0.239
 ANA-positive, n (%)211(95.5)105(94.6)106(96.4)0.423
 ACA-positive, n (%)87(39.4)44(39.6)43(39.1)0.933
 ATA-positive, n (%)54(24.4)23(20.7)31(28.3)0.197
 DLCO%64.7 ± 19.159.5 ± 18.570.0 ± 18.2<0.001
 VO2max%87.0 ± 23.785.1 ± 22.488.9 ± 24.90.262
Immunosuppressive medication, n (%)
 Corticosteroids31(13.5)18(15.7)13(11.3)0.334
 Cyclophosphamide3(1.3)2(1.7)1(0.9)0.561
 Methotrexate29(12.6)18(15.7)11(9.6)0.164
 Azathioprine11(4.8)9(7.8)2(1.7)0.031
Cardiovascular medications, n (%)
 ACEI/ARB85(37.0)57(49.6)28(24.3)<0.001
 CCB108(47.0)64(55.7)44(38.3)0.008
 Beta-blockers21(9.1)17(14.8)4(3.5)0.003
 Diuretics30(13.0)22(19.1)8(7.0)0.006

Values are mean ± standard deviation if normally distributed and median (interquartile range) if not normally distributed.

ACA, anti-centromere antibodies; ACEI, angiotensin-converting enzyme inhibitor; AF, atrial fibrillation; ANA, anti-nuclear antibodies; ARB, angiotensin II receptor blocker; ATA, anti-topoisomerase antibodies; BMI, body mass index; CAD, coronary artery disease; CCB, calcium channel blocker; CPK, creatine kinase; CRP, C-reactive protein; DLCO, diffusion capacity of the lung for carbon monoxide; DM, diabetes mellitus; eGFR, estimated glomerular filtration rate; HTN, hypertension; ILD, interstitial lung disease; MRSS, modified Rodnan skin score; NT-proBNP, N-terminal pro-brain natriuretic peptide; PAH, pulmonary artery hypertension; VO2max, maximum O2 uptake during exercise.

Table 1

Baseline clinical characteristics of the total SSc population and divided based on the median value of EAT mass

Total populationEAT ≥ 67 gEAT < 67 gP-value
(n = 230)(n = 115)(n = 115)
Clinical variables
 Age, y53 ± 1560 ± 1145 ± 14<0.001
 Male, n (%)33(14.3)24(20.9)9(7.8)0.005
 BMI, kg/m224.6 ± 4.326.0 ± 4.423.1 ± 3.6<0.001
 Diffuse SSc, n (%)51(23.1)24(21.6)27(24.5)0.458
 HTN, n (%)58(25.2)39(33.9)19(16.5)0.002
 DM, n (%)11(4.8)9(7.8)2(1.7)0.031
 CAD, n (%)8(3.5)6(5.2)1(0.8)0.150
 AF, n (%)7(3.2)6(5.5)1(0.9)0.055
 History of smoking, n (%)110(47.8)61(53.0)49(42.6)0.113
 mRSS5 ± 86 ± 94 ± 50.124
 ILD, n (%)91(41.2)52(46.8)39(35.5)0.085
 PAH, n (%)8(3.6)7(6.3)1(0.9)0.032
 Time since Raynaud, y13 ± 1415 ± 1511 ± 120.013
 Time since diagnosis, y5 ± 76 ± 84 ± 60.014
Laboratory findings
 eGFR, mL/min/1.73 m288(72–103)82(69–96)95(79–106)<0.001
 NT-proBNP, ng/L95.8(53.6–170.5)107.6(59.2–252.4)79.0(49.1–137.6)0.005
 CRP, mg/L0(0–5)0(0–6)0(0–0)0.002
 CPK, U/L83(61–118)93(60–134)79(61–109)0.239
 ANA-positive, n (%)211(95.5)105(94.6)106(96.4)0.423
 ACA-positive, n (%)87(39.4)44(39.6)43(39.1)0.933
 ATA-positive, n (%)54(24.4)23(20.7)31(28.3)0.197
 DLCO%64.7 ± 19.159.5 ± 18.570.0 ± 18.2<0.001
 VO2max%87.0 ± 23.785.1 ± 22.488.9 ± 24.90.262
Immunosuppressive medication, n (%)
 Corticosteroids31(13.5)18(15.7)13(11.3)0.334
 Cyclophosphamide3(1.3)2(1.7)1(0.9)0.561
 Methotrexate29(12.6)18(15.7)11(9.6)0.164
 Azathioprine11(4.8)9(7.8)2(1.7)0.031
Cardiovascular medications, n (%)
 ACEI/ARB85(37.0)57(49.6)28(24.3)<0.001
 CCB108(47.0)64(55.7)44(38.3)0.008
 Beta-blockers21(9.1)17(14.8)4(3.5)0.003
 Diuretics30(13.0)22(19.1)8(7.0)0.006
Total populationEAT ≥ 67 gEAT < 67 gP-value
(n = 230)(n = 115)(n = 115)
Clinical variables
 Age, y53 ± 1560 ± 1145 ± 14<0.001
 Male, n (%)33(14.3)24(20.9)9(7.8)0.005
 BMI, kg/m224.6 ± 4.326.0 ± 4.423.1 ± 3.6<0.001
 Diffuse SSc, n (%)51(23.1)24(21.6)27(24.5)0.458
 HTN, n (%)58(25.2)39(33.9)19(16.5)0.002
 DM, n (%)11(4.8)9(7.8)2(1.7)0.031
 CAD, n (%)8(3.5)6(5.2)1(0.8)0.150
 AF, n (%)7(3.2)6(5.5)1(0.9)0.055
 History of smoking, n (%)110(47.8)61(53.0)49(42.6)0.113
 mRSS5 ± 86 ± 94 ± 50.124
 ILD, n (%)91(41.2)52(46.8)39(35.5)0.085
 PAH, n (%)8(3.6)7(6.3)1(0.9)0.032
 Time since Raynaud, y13 ± 1415 ± 1511 ± 120.013
 Time since diagnosis, y5 ± 76 ± 84 ± 60.014
Laboratory findings
 eGFR, mL/min/1.73 m288(72–103)82(69–96)95(79–106)<0.001
 NT-proBNP, ng/L95.8(53.6–170.5)107.6(59.2–252.4)79.0(49.1–137.6)0.005
 CRP, mg/L0(0–5)0(0–6)0(0–0)0.002
 CPK, U/L83(61–118)93(60–134)79(61–109)0.239
 ANA-positive, n (%)211(95.5)105(94.6)106(96.4)0.423
 ACA-positive, n (%)87(39.4)44(39.6)43(39.1)0.933
 ATA-positive, n (%)54(24.4)23(20.7)31(28.3)0.197
 DLCO%64.7 ± 19.159.5 ± 18.570.0 ± 18.2<0.001
 VO2max%87.0 ± 23.785.1 ± 22.488.9 ± 24.90.262
Immunosuppressive medication, n (%)
 Corticosteroids31(13.5)18(15.7)13(11.3)0.334
 Cyclophosphamide3(1.3)2(1.7)1(0.9)0.561
 Methotrexate29(12.6)18(15.7)11(9.6)0.164
 Azathioprine11(4.8)9(7.8)2(1.7)0.031
Cardiovascular medications, n (%)
 ACEI/ARB85(37.0)57(49.6)28(24.3)<0.001
 CCB108(47.0)64(55.7)44(38.3)0.008
 Beta-blockers21(9.1)17(14.8)4(3.5)0.003
 Diuretics30(13.0)22(19.1)8(7.0)0.006

Values are mean ± standard deviation if normally distributed and median (interquartile range) if not normally distributed.

ACA, anti-centromere antibodies; ACEI, angiotensin-converting enzyme inhibitor; AF, atrial fibrillation; ANA, anti-nuclear antibodies; ARB, angiotensin II receptor blocker; ATA, anti-topoisomerase antibodies; BMI, body mass index; CAD, coronary artery disease; CCB, calcium channel blocker; CPK, creatine kinase; CRP, C-reactive protein; DLCO, diffusion capacity of the lung for carbon monoxide; DM, diabetes mellitus; eGFR, estimated glomerular filtration rate; HTN, hypertension; ILD, interstitial lung disease; MRSS, modified Rodnan skin score; NT-proBNP, N-terminal pro-brain natriuretic peptide; PAH, pulmonary artery hypertension; VO2max, maximum O2 uptake during exercise.

Echocardiography

The echocardiographic variables of the study population are displayed in Table 2. Patients with increased EAT had significantly larger IVST and PWT as compared with those with less EAT. Furthermore, patients with increase EAT had larger LAVI, larger LVMI, reduced mean Eʹ, increased E/Eʹ, and higher TR peak velocity compared with those with less EAT. No significant differences in LV end-diastolic volume, LV end-systolic volume, and LVEF were observed between the two groups.

Table 2

Baseline echocardiographic and CT characteristics in the total population and divided based on the median value of EAT

Total populationEAT ≥ 67 gEAT < 67 gP-value
(n = 230)(n = 115)(n = 115)
Echocardiographic parameters
 LVEDV, mL83.5 ± 24.382.3 ± 26.284.7 ± 22.70.562
 LVESV, mL32.6 ± 13.632.2 ± 15.032.9 ± 12.20.694
 LVEF, %61.7 ± 6.861.7 ± 7.061.6 ± 7.00.921
 LVMI, g/m2148.2 ± 42.3156.4 ± 39.6140.0 ± 43.40.003
 IVST, mm8.9 ± 1.89.2 ± 1.78.7 ± 1.90.048
 PWT, mm8.9 ± 1.69.1 ± 1.68.7 ± 1.70.037
 LAVI, mL/m223.9 ± 9.427.0 ± 9.920.8 ± 7.9<0.001
 E/A1.1 ± 0.41.0 ± 0.31.3 ± 0.4<0.001
 Eʹ, cm/s9.0 ± 2.97.5 ± 2.210.5 ± 2.8<0.001
 E/E′9.5 ± 3.910.8 ± 4.38.2 ± 2.8<0.001
 TRV, m/s2.4 ± 0.42.5 ± 0.42.3 ± 0.3<0.001
CT parameter
 EAT, g67 (45–101)100 (82–122)45 (35–56)<0.001
Total populationEAT ≥ 67 gEAT < 67 gP-value
(n = 230)(n = 115)(n = 115)
Echocardiographic parameters
 LVEDV, mL83.5 ± 24.382.3 ± 26.284.7 ± 22.70.562
 LVESV, mL32.6 ± 13.632.2 ± 15.032.9 ± 12.20.694
 LVEF, %61.7 ± 6.861.7 ± 7.061.6 ± 7.00.921
 LVMI, g/m2148.2 ± 42.3156.4 ± 39.6140.0 ± 43.40.003
 IVST, mm8.9 ± 1.89.2 ± 1.78.7 ± 1.90.048
 PWT, mm8.9 ± 1.69.1 ± 1.68.7 ± 1.70.037
 LAVI, mL/m223.9 ± 9.427.0 ± 9.920.8 ± 7.9<0.001
 E/A1.1 ± 0.41.0 ± 0.31.3 ± 0.4<0.001
 Eʹ, cm/s9.0 ± 2.97.5 ± 2.210.5 ± 2.8<0.001
 E/E′9.5 ± 3.910.8 ± 4.38.2 ± 2.8<0.001
 TRV, m/s2.4 ± 0.42.5 ± 0.42.3 ± 0.3<0.001
CT parameter
 EAT, g67 (45–101)100 (82–122)45 (35–56)<0.001

Values are mean ± standard deviation if normally distributed and median (interquartile range) if abnormally distributed.

CT, computed tomography; E, peak early diastolic mitral flow velocity; E′, peak early diastolic mitral annular tissue velocity; EAT, epicardial adipose tissue; IVST, interventricular septum thickness; LAVI, left atrial volume index; LV, left ventricular; LVEDV, left ventricle end-diastolic volume; LVEF, left ventricular ejection fraction; LVESV, left ventricular end-systolic volume; LVMI, left ventricular mass index; PWT, posterior LV wall thickness; TRV, tricuspid regurgitation velocity.

Table 2

Baseline echocardiographic and CT characteristics in the total population and divided based on the median value of EAT

Total populationEAT ≥ 67 gEAT < 67 gP-value
(n = 230)(n = 115)(n = 115)
Echocardiographic parameters
 LVEDV, mL83.5 ± 24.382.3 ± 26.284.7 ± 22.70.562
 LVESV, mL32.6 ± 13.632.2 ± 15.032.9 ± 12.20.694
 LVEF, %61.7 ± 6.861.7 ± 7.061.6 ± 7.00.921
 LVMI, g/m2148.2 ± 42.3156.4 ± 39.6140.0 ± 43.40.003
 IVST, mm8.9 ± 1.89.2 ± 1.78.7 ± 1.90.048
 PWT, mm8.9 ± 1.69.1 ± 1.68.7 ± 1.70.037
 LAVI, mL/m223.9 ± 9.427.0 ± 9.920.8 ± 7.9<0.001
 E/A1.1 ± 0.41.0 ± 0.31.3 ± 0.4<0.001
 Eʹ, cm/s9.0 ± 2.97.5 ± 2.210.5 ± 2.8<0.001
 E/E′9.5 ± 3.910.8 ± 4.38.2 ± 2.8<0.001
 TRV, m/s2.4 ± 0.42.5 ± 0.42.3 ± 0.3<0.001
CT parameter
 EAT, g67 (45–101)100 (82–122)45 (35–56)<0.001
Total populationEAT ≥ 67 gEAT < 67 gP-value
(n = 230)(n = 115)(n = 115)
Echocardiographic parameters
 LVEDV, mL83.5 ± 24.382.3 ± 26.284.7 ± 22.70.562
 LVESV, mL32.6 ± 13.632.2 ± 15.032.9 ± 12.20.694
 LVEF, %61.7 ± 6.861.7 ± 7.061.6 ± 7.00.921
 LVMI, g/m2148.2 ± 42.3156.4 ± 39.6140.0 ± 43.40.003
 IVST, mm8.9 ± 1.89.2 ± 1.78.7 ± 1.90.048
 PWT, mm8.9 ± 1.69.1 ± 1.68.7 ± 1.70.037
 LAVI, mL/m223.9 ± 9.427.0 ± 9.920.8 ± 7.9<0.001
 E/A1.1 ± 0.41.0 ± 0.31.3 ± 0.4<0.001
 Eʹ, cm/s9.0 ± 2.97.5 ± 2.210.5 ± 2.8<0.001
 E/E′9.5 ± 3.910.8 ± 4.38.2 ± 2.8<0.001
 TRV, m/s2.4 ± 0.42.5 ± 0.42.3 ± 0.3<0.001
CT parameter
 EAT, g67 (45–101)100 (82–122)45 (35–56)<0.001

Values are mean ± standard deviation if normally distributed and median (interquartile range) if abnormally distributed.

CT, computed tomography; E, peak early diastolic mitral flow velocity; E′, peak early diastolic mitral annular tissue velocity; EAT, epicardial adipose tissue; IVST, interventricular septum thickness; LAVI, left atrial volume index; LV, left ventricular; LVEDV, left ventricle end-diastolic volume; LVEF, left ventricular ejection fraction; LVESV, left ventricular end-systolic volume; LVMI, left ventricular mass index; PWT, posterior LV wall thickness; TRV, tricuspid regurgitation velocity.

Association between EAT and cardiac function

Univariable linear regression analyses showed that EAT was significantly associated to LV diastolic function parameters including LAVI, LVMI, Eʹ, E/Eʹ, and TR peak velocity, while there was no significant association between EAT and LVEF (see Supplementary data online, Table S1).

At the multivariable linear regression analyses, EAT showed an independent association with LV diastolic function parameters including LAVI (B = 0.031, P = 0.025), LVMI (B = 0.139, P = 0.036), E/E′ (B = 0.025, P < 0.001), and Eʹ (B = −0.012, P < 0.001) after correction for clinical characteristics including age, sex, hypertension, diabetes, and eGFR, indicating that increased EAT is associated with impaired LV diastolic function (Table 3).

Table 3

Multivariable linear regression analysis for EAT and echocardiographic parameters of LV diastolic function (therefore corrected for demographics and comorbidities)

LAVILVMIE/E′TRV
B coeffP-valueB coeffP-valueB coeffP-valueB coeffP-valueB coeffP-value
Age0.261<0.0010.3590.1010.079<0.001−0.111<0.0010.158<0.001
Male0.6910.68128.139<0.001−0.7800.222−0.0280.9440.8510.615
DM0.2710.9239.0680.9190.9950.327−0.7050.2910.0850.974
HTN−2.9080.034−0.3880.9520.5320.303−0.2730.412−2.2620.075
eGFR0.0140.5720.1050.001−0.0010.9210.0070.247−0.0170.461
EAT0.0310.0250.1390.0360.025<0.001−0.012<0.0010.0250.112
LAVILVMIE/E′TRV
B coeffP-valueB coeffP-valueB coeffP-valueB coeffP-valueB coeffP-value
Age0.261<0.0010.3590.1010.079<0.001−0.111<0.0010.158<0.001
Male0.6910.68128.139<0.001−0.7800.222−0.0280.9440.8510.615
DM0.2710.9239.0680.9190.9950.327−0.7050.2910.0850.974
HTN−2.9080.034−0.3880.9520.5320.303−0.2730.412−2.2620.075
eGFR0.0140.5720.1050.001−0.0010.9210.0070.247−0.0170.461
EAT0.0310.0250.1390.0360.025<0.001−0.012<0.0010.0250.112

DM, diabetes mellitus; E, peak early diastolic mitral flow velocity; E′, peak early diastolic mitral annular tissue velocity; EAT, epicardial adipose tissue; eGFR, estimated glomerular filtration rate; HTN, hypertension; LAVI, left atrial volume index; LVMI, left ventricular mass index; TRV, tricuspid regurgitation velocity.

Table 3

Multivariable linear regression analysis for EAT and echocardiographic parameters of LV diastolic function (therefore corrected for demographics and comorbidities)

LAVILVMIE/E′TRV
B coeffP-valueB coeffP-valueB coeffP-valueB coeffP-valueB coeffP-value
Age0.261<0.0010.3590.1010.079<0.001−0.111<0.0010.158<0.001
Male0.6910.68128.139<0.001−0.7800.222−0.0280.9440.8510.615
DM0.2710.9239.0680.9190.9950.327−0.7050.2910.0850.974
HTN−2.9080.034−0.3880.9520.5320.303−0.2730.412−2.2620.075
eGFR0.0140.5720.1050.001−0.0010.9210.0070.247−0.0170.461
EAT0.0310.0250.1390.0360.025<0.001−0.012<0.0010.0250.112
LAVILVMIE/E′TRV
B coeffP-valueB coeffP-valueB coeffP-valueB coeffP-valueB coeffP-value
Age0.261<0.0010.3590.1010.079<0.001−0.111<0.0010.158<0.001
Male0.6910.68128.139<0.001−0.7800.222−0.0280.9440.8510.615
DM0.2710.9239.0680.9190.9950.327−0.7050.2910.0850.974
HTN−2.9080.034−0.3880.9520.5320.303−0.2730.412−2.2620.075
eGFR0.0140.5720.1050.001−0.0010.9210.0070.247−0.0170.461
EAT0.0310.0250.1390.0360.025<0.001−0.012<0.0010.0250.112

DM, diabetes mellitus; E, peak early diastolic mitral flow velocity; E′, peak early diastolic mitral annular tissue velocity; EAT, epicardial adipose tissue; eGFR, estimated glomerular filtration rate; HTN, hypertension; LAVI, left atrial volume index; LVMI, left ventricular mass index; TRV, tricuspid regurgitation velocity.

EAT and outcome

During a median follow-up of 8 years (IQR 7–10 years), 42 patients (18%) died. Survival analysis demonstrated that patients with increased EAT (≥67 g) experienced significantly higher rates of all-cause mortality (29.6% vs. 7.0%; P < 0.001) (Figure 2), compared with patients with less EAT (<67 g).

Kaplan–Meier curves for all-mortality. Time to death according to baseline EAT < 67 g (less) and EAT ≥ 67 g (more). The lighter shaded colour around the curves indicates a 95% confidence interval. EAT, epicardial adipose tissue.
Figure 2

Kaplan–Meier curves for all-mortality. Time to death according to baseline EAT < 67 g (less) and EAT ≥ 67 g (more). The lighter shaded colour around the curves indicates a 95% confidence interval. EAT, epicardial adipose tissue.

Univariable cox regression analyses demonstrated significant associations between all-cause mortality and age, ILD, DLCO, CRP, NT-proBNP, and EAT mass (Table 4). Because of the total number of events, five variables including age, DLCO (preferred over ILD to represent pulmonary involvement), NT-proBNP, CRP, and EAT mass were ultimately included in the multivariable Cox proportional hazards backward stepwise model. In this analysis, age, DLCO, and EAT mass were independently associated with all-cause mortality (Table 4). Of note, when considering a combined endpoint of all-cause mortality, heart failure hospitalization, myocardial infarction, and arrhythmias, at the multivariable analysis EAT remained also independently associated with the endpoint [HR 1.005 (1.001–1.009), P = 0.008] together with NT-proBNP and DLCO (data not shown in the tables).

Table 4

Univariable and multivariable analyses to evaluate the association between EAT and all-cause mortality

VariableUnivariable analysisMultivariable analysis
HR (95% CI)P-valueHR 95% CIP-value
Age (per one year increase)1.072 (1.044–1.101)<0.0011.047 (1.017–1.079)0.002
Male (yes/no)1.877 (0.922–3.821)0.082
Diffuse SSc (yes/no)1.785 (0.913–3.489)0.090
PAH (yes/no)1.436 (0.990–2.084)0.056
ILD (yes/no)3.738 (1.901–7.331)<0.001
Time since Raynaud (per one year increase)1.001 (0.979–1.023)0.918
Autoantibodies
0Reference
10.949 (0.251–3.578)0.938
20.556 (0.153–2.023)0.373
31.379 (0.411–4.629)0.603
CRP (per one unit increase)1.046 (1.024–1.069)<0.001NS
eGFR (per one unit increase)0.990 (0.978–1.001)0.087
DLCO (per one unit increase)0.975 (0.968–0.981)<0.0010.978 (0.971–0.985)<0.001
LVEF (per one unit increase)0.976 (0.937–1.017)0.249
NT-proBNP (per one unit increase)1.121 (1.083–10160)<0.001NS
EAT (per gram increase)1.009 (1.006–1.012)<0.0011.006 (1.001–1.010)0.010
VariableUnivariable analysisMultivariable analysis
HR (95% CI)P-valueHR 95% CIP-value
Age (per one year increase)1.072 (1.044–1.101)<0.0011.047 (1.017–1.079)0.002
Male (yes/no)1.877 (0.922–3.821)0.082
Diffuse SSc (yes/no)1.785 (0.913–3.489)0.090
PAH (yes/no)1.436 (0.990–2.084)0.056
ILD (yes/no)3.738 (1.901–7.331)<0.001
Time since Raynaud (per one year increase)1.001 (0.979–1.023)0.918
Autoantibodies
0Reference
10.949 (0.251–3.578)0.938
20.556 (0.153–2.023)0.373
31.379 (0.411–4.629)0.603
CRP (per one unit increase)1.046 (1.024–1.069)<0.001NS
eGFR (per one unit increase)0.990 (0.978–1.001)0.087
DLCO (per one unit increase)0.975 (0.968–0.981)<0.0010.978 (0.971–0.985)<0.001
LVEF (per one unit increase)0.976 (0.937–1.017)0.249
NT-proBNP (per one unit increase)1.121 (1.083–10160)<0.001NS
EAT (per gram increase)1.009 (1.006–1.012)<0.0011.006 (1.001–1.010)0.010

Autoantibodies is a categorical variable (0 = ANA−, 1 = ATA+, 2 = ACA+, and 3 = other).

ACA, anti-centromere antibodies; ANA, anti-nuclear antibodies; ATA, anti-topoisomerase antibodies; DLCO, diffusion capacity of the lung for carbon monoxide; DM, diabetes mellitus; EAT, epicardial adipose tissue; E′, peak early diastolic mitral annular tissue velocity; eGFR, estimated glomerular filtration rate; HTN, hypertension; ILD, interstitial lung disease; LVEF, left ventricular ejection fraction; NT-proBNP, N-terminal pro-brain natriuretic peptide; PAH, pulmonary artery hypertension; SSc, systemic sclerosis.

Table 4

Univariable and multivariable analyses to evaluate the association between EAT and all-cause mortality

VariableUnivariable analysisMultivariable analysis
HR (95% CI)P-valueHR 95% CIP-value
Age (per one year increase)1.072 (1.044–1.101)<0.0011.047 (1.017–1.079)0.002
Male (yes/no)1.877 (0.922–3.821)0.082
Diffuse SSc (yes/no)1.785 (0.913–3.489)0.090
PAH (yes/no)1.436 (0.990–2.084)0.056
ILD (yes/no)3.738 (1.901–7.331)<0.001
Time since Raynaud (per one year increase)1.001 (0.979–1.023)0.918
Autoantibodies
0Reference
10.949 (0.251–3.578)0.938
20.556 (0.153–2.023)0.373
31.379 (0.411–4.629)0.603
CRP (per one unit increase)1.046 (1.024–1.069)<0.001NS
eGFR (per one unit increase)0.990 (0.978–1.001)0.087
DLCO (per one unit increase)0.975 (0.968–0.981)<0.0010.978 (0.971–0.985)<0.001
LVEF (per one unit increase)0.976 (0.937–1.017)0.249
NT-proBNP (per one unit increase)1.121 (1.083–10160)<0.001NS
EAT (per gram increase)1.009 (1.006–1.012)<0.0011.006 (1.001–1.010)0.010
VariableUnivariable analysisMultivariable analysis
HR (95% CI)P-valueHR 95% CIP-value
Age (per one year increase)1.072 (1.044–1.101)<0.0011.047 (1.017–1.079)0.002
Male (yes/no)1.877 (0.922–3.821)0.082
Diffuse SSc (yes/no)1.785 (0.913–3.489)0.090
PAH (yes/no)1.436 (0.990–2.084)0.056
ILD (yes/no)3.738 (1.901–7.331)<0.001
Time since Raynaud (per one year increase)1.001 (0.979–1.023)0.918
Autoantibodies
0Reference
10.949 (0.251–3.578)0.938
20.556 (0.153–2.023)0.373
31.379 (0.411–4.629)0.603
CRP (per one unit increase)1.046 (1.024–1.069)<0.001NS
eGFR (per one unit increase)0.990 (0.978–1.001)0.087
DLCO (per one unit increase)0.975 (0.968–0.981)<0.0010.978 (0.971–0.985)<0.001
LVEF (per one unit increase)0.976 (0.937–1.017)0.249
NT-proBNP (per one unit increase)1.121 (1.083–10160)<0.001NS
EAT (per gram increase)1.009 (1.006–1.012)<0.0011.006 (1.001–1.010)0.010

Autoantibodies is a categorical variable (0 = ANA−, 1 = ATA+, 2 = ACA+, and 3 = other).

ACA, anti-centromere antibodies; ANA, anti-nuclear antibodies; ATA, anti-topoisomerase antibodies; DLCO, diffusion capacity of the lung for carbon monoxide; DM, diabetes mellitus; EAT, epicardial adipose tissue; E′, peak early diastolic mitral annular tissue velocity; eGFR, estimated glomerular filtration rate; HTN, hypertension; ILD, interstitial lung disease; LVEF, left ventricular ejection fraction; NT-proBNP, N-terminal pro-brain natriuretic peptide; PAH, pulmonary artery hypertension; SSc, systemic sclerosis.

Incremental prognostic value of EAT

The addition of EAT to a Cox model for all-cause mortality that included age and DLCO resulted in a significantly improved model fit as assessed by using the likelihood ratio test (P = 0.019).

Reproducibility

Measurement of EAT mass in 10 randomly selected patients showed excellent intra- and inter-observer reproducibility. The inter-observer ICC was 0.993 (95% CI 0.815–0.999), and the intra-observer ICC was 0.997 (95% CI 0.963–0.999) for EAT measurement.

Discussion

The present study examined a large cohort of SSc patients, exploring for the first time the association between EAT, LV function, and clinical outcomes. The extent of EAT measured by non-contrast CT was significantly and independently associated with LV diastolic dysfunction and with higher mortality rates.

Association between EAT, severity of SSc, and LV diastolic dysfunction

The current study showed that patients with increased EAT had higher levels of CRP, longer time to diagnosis (or to Raynaud), higher levels of NT-proBNP, lower eGFR values, and DLCO, suggesting an association between EAT and severity of SSc and further organ damage.

Previous studies showed that several systemic inflammatory diseases, including rheumatoid arthritis, systemic lupus erythematosus, and psoriasis, are accompanied by an increase in EAT mass, which is proportional to the duration and severity of the underlying disease.25–27 In SSc, increased EAT volume was also suggested to be associated with the presence and severity of the disease.28 In particular, the activation of aldosterone, leptin, and neprilysin observed in systemic inflammatory disease may facilitate the accumulation and dysfunction of EAT.8 In addition to systemic inflammation, numerous metabolic disorders (such as diabetes, insulin resistance, and metabolic syndrome) are accompanied by the expansion and inflammation of visceral adipose tissue29–32; therefore, assessment of the EAT extent should take into account also the presence of these comorbidities.

On the other hand, because of its endocrine and paracrine activity, secreting pro-inflammatory and anti-inflammatory cytokines and chemokines, EAT may play a role in the pathophysiology of cardiac damage. Previous studies have suggested EAT to be related to coronary atherosclerosis, microcirculatory dysfunction, myocardial fibrosis, and LV hypertrophy,33–39 and therefore to be involved in the development of coronary artery disease, atrial fibrillation, and most importantly HFpEF.14–16 This study shows that EAT is associated with echocardiographic parameters of LV diastolic function including LAVI, LVMI, Eʹ, and E/Eʹ, and independently of other risk factors and comorbidities. It therefore supports the hypothesis that EAT plays a role in the cardiovascular involvement in SSc, and mainly in the development of LV diastolic dysfunction (overt LV systolic dysfunction was not observed in this cohort).

EAT and prognosis

Many chronic systemic inflammatory disorders are accompanied by an increased risk of heart failure, particularly HFpEF,24,40–43 and this risk is independent of macrovascular coronary heart disease. Particularly in SSc, presence of LV diastolic dysfunction has been shown to be of independent prognostic value.6 As above-mentioned, EAT, by playing a role in the development of myocardial fibrosis, LV hypertrophy, and impaired distensibility, may be therefore an important marker of cardiac involvement. In addition, EAT may also help general risk stratification of SSc patients by reflecting the severity of the disease and of multiorgan involvement. To the best of our knowledge, this study is the first demonstrating an association between EAT and clinical outcome in a large SSc cohort, importantly after correction for clinical characteristics and LV function. Considering that most patients with SSc undergo a non-contrast thorax CT, EAT measure could be easily implemented in the routine assessment of these patients and have impact in their management. Importantly, EAT could also represent a novel therapeutic target: because of its rapid metabolism and potential for modification, targeting EAT with drugs such as glucagon-like peptide 1 and sodium-glucose cotransporter-2 inhibitor agonists may reshape the landscape of pharmacotherapy, providing cardiovascular protection beyond just SSc management.

Study limitations

This study has all of the inherent limitations associated with a single centre, retrospective study design. However, the systematic and standardized evaluation of all the patients, according to a prospective cohort study design, supports the robustness of the findings; of note, only patients with the CT scan and the echocardiography performed at the same visit were included and patients with inadequate CT quality were excluded, which may have introduced a selection bias. Finally, a direct demonstration of the EAT inflammatory activity could not be provided neither by imaging or laboratory tests. The role that EAT plays in pathogenesis of cardiac involvement in SSc requires therefore future investigation.

Conclusion

In patients with SSc, EAT is associated with LV diastolic dysfunction and higher mortality rate, independently of traditional risk factors. The mechanisms underlying this association warranted to be explored.

Supplementary data

Supplementary data is available at European Heart Journal - Imaging Methods and Practice online.

Acknowledgements

The authors deeply thank all the staff members of the non-invasive cardiac imaging laboratories, the Department of Cardiology and the Department of Rheumatology at the Leiden University Medical Center for their effort in this article.

Consent: All patients gave written informed consent to this registry in accordance with the approval of the Medical Ethics Committee of Leiden University Medical Center.

Funding

S.C.B. received funding from European Society of Cardiology (ESC Research Grant App000080404). X.W.is supported by a research grant from the University of Turku.

Data availability

The data underlying this article will be shared on reasonable request to the corresponding author.

Lead author biography

Nina Ajmone Marsan, MD, PhD, FESC

Assistant Professor, senior cardiologist, expert in valvular heart disease and cardiovascular imaging, Leiden University Medical Center, The Netherlands

graphicDr Nina Ajmone Marsan is a senior cardiologist, specializing in non-invasive cardiovascular imaging and currently working at the Leiden University Medical Center (LUMC, The Netherlands), where she is Assistant Professor and supervises several Dutch and international PhD students. Her main research topic is the clinical applications of imaging in cardiovascular disease with special focus on valvular heart disease (particularly for patients referred for trans-catheter interventions), cardiomyopathies, and systemic disease. She is also the principal investigator of translational research projects in mitral and aortic valve disease and the coordinator of the LUMC Heart Valve Center. She is in the board of EACVI as Vice president elect 2022–2024 and as Chair of the Task Forces and Shared Interest Groups (2022–2024). She is Associate Editor of the European Heart Journal of the European Heart Journal Cardiovascular Imaging of the Journal of the American Society of Echocardiography. She is in the Congress Program Committee of the ESC 2023 and course director of the EACVI TOE teaching course, of the PCR Madrid Imaging Valve and of the ‘3D echocardiography: meet the experts’ course.

References

1

Denton
 
CP
,
Khanna
 
D
.
Systemic sclerosis
.
Lancet
 
2017
;
390
:
1685
99
.

2

Elhai
 
M
,
Meune
 
C
,
Boubaya
 
M
,
Avouac
 
J
,
Hachulla
 
E
,
Balbir-Gurman
 
A
 et al.  
Mapping and predicting mortality from systemic sclerosis
.
Ann Rheum Dis
 
2020
;
79
:
1608
15
.

3

Butt
 
SA
,
Jeppesen
 
JL
,
Torp-Pedersen
 
C
,
Sam
 
F
,
Gislason
 
GH
,
Jacobsen
 
S
 et al.  
Cardiovascular manifestations of systemic sclerosis: a Danish nationwide cohort study
.
J Am Heart Assoc
 
2019
;
8
:
e013405
.

4

Meune
 
C
,
Avouac
 
J
,
Wahbi
 
K
,
Cabanes
 
L
,
Wipff
 
J
,
Mouthon
 
L
 et al.  
Cardiac involvement in systemic sclerosis assessed by tissue-Doppler echocardiography during routine care: a controlled study of 100 consecutive patients
.
Arthritis Rheum
 
2008
;
58
:
1803
9
.

5

de Groote
 
P
,
Gressin
 
V
,
Hachulla
 
E
,
Carpentier
 
P
,
Guillevin
 
L
,
Kahan
 
A
 et al.  
Evaluation of cardiac abnormalities by Doppler echocardiography in a large nationwide multicentric cohort of patients with systemic sclerosis
.
Ann Rheum Dis
 
2008
;
67
:
31
6
.

6

Tennøe
 
AH
,
Murbræch
 
K
,
Andreassen
 
JC
,
Fretheim
 
H
,
Garen
 
T
,
Gude
 
E
 et al.  
Left ventricular diastolic dysfunction predicts mortality in patients with systemic sclerosis
.
J Am Coll Cardiol
 
2018
;
72
:
1804
13
.

7

Rodríguez-Reyna
 
TS
,
Morelos-Guzman
 
M
,
Hernández-Reyes
 
P
,
Montero-Duarte
 
K
,
Martínez-Reyes
 
C
,
Reyes-Utrera
 
C
 et al.  
Assessment of myocardial fibrosis and microvascular damage in systemic sclerosis by magnetic resonance imaging and coronary angiotomography
.
Rheumatology
 
2015
;
54
:
647
54
.

8

Packer
 
M
,
Lam
 
CS
,
Lund
 
LH
,
Maurer
 
MS
,
Borlaug
 
BA
.
Characterization of the inflammatory-metabolic phenotype of heart failure with a preserved ejection fraction: a hypothesis to explain influence of sex on the evolution and potential treatment of the disease
.
Eur J Heart Fail
 
2020
;
22
:
1551
67
.

9

Wong
 
CY
,
O’Moore-Sullivan
 
T
,
Leano
 
R
,
Byrne
 
N
,
Beller
 
E
,
Marwick
 
TH
.
Alterations of left ventricular myocardial characteristics associated with obesity
.
Circulation
 
2004
;
110
:
3081
7
.

10

Packer
 
M
.
Epicardial adipose tissue may mediate deleterious effects of obesity and inflammation on the myocardium
.
J Am Coll Cardiol
 
2018
;
71
:
2360
72
.

11

Mohammed
 
SF
,
Hussain
 
S
,
Mirzoyev
 
SA
,
Edwards
 
WD
,
Maleszewski
 
JJ
,
Redfield
 
MM
.
Coronary microvascular rarefaction and myocardial fibrosis in heart failure with preserved ejection fraction
.
Circulation
 
2015
;
131
:
550
9
.

12

Paulus
 
WJ
,
Tschöpe
 
C
.
A novel paradigm for heart failure with preserved ejection fraction: comorbidities drive myocardial dysfunction and remodeling through coronary microvascular endothelial inflammation
.
J Am Coll Cardiol
 
2013
;
62
:
263
71
.

13

Obokata
 
M
,
Reddy
 
YN
,
Pislaru
 
SV
,
Melenovsky
 
V
,
Borlaug
 
BA
.
Evidence supporting the existence of a distinct obese phenotype of heart failure with preserved ejection fraction
.
Circulation
 
2017
;
136
:
6
19
.

14

Mahabadi
 
AA
,
Massaro
 
JM
,
Rosito
 
GA
,
Levy
 
D
,
Murabito
 
JM
,
Wolf
 
PA
 et al.  
Association of pericardial fat, intrathoracic fat, and visceral abdominal fat with cardiovascular disease burden: the Framingham Heart Study
.
Eur Heart J
 
2009
;
30
:
850
6
.

15

Al Chekakie
 
MO
,
Welles
 
CC
,
Metoyer
 
R
,
Ibrahim
 
A
,
Shapira
 
AR
,
Cytron
 
J
 et al.  
Pericardial fat is independently associated with human atrial fibrillation
.
J Am Coll Cardiol
 
2010
;
56
:
784
8
.

16

van Woerden
 
G
,
Gorter
 
TM
,
Westenbrink
 
BD
,
Willems
 
TP
,
van Veldhuisen
 
DJ
,
Rienstra
 
M
.
Epicardial fat in heart failure patients with mid-range and preserved ejection fraction
.
Eur J Heart Fail
 
2018
;
20
:
1559
66
.

17

Meijs
 
J
,
Schouffoer
 
AA
,
Marsan
 
NA
,
Kroft
 
LJ
,
Stijnen
 
T
,
Ninaber
 
MK
 et al.  
Therapeutic and diagnostic outcomes of a standardised, comprehensive care pathway for patients with systemic sclerosis
.
RMD open
 
2016
;
2
:
e000159
.

18

Schouffoer
 
A
,
Ninaber
 
M
,
Beaart-Van De Voorde
 
L
,
Van Der Giesen
 
F
,
De Jong
 
Z
,
Stolk
 
J
 et al.  
Randomized comparison of a multidisciplinary team care program with usual care in patients with systemic sclerosis
.
Arthritis Care Res
 
2011
;
63
:
909
17
.

19

Van Den Hoogen
 
F
,
Khanna
 
D
,
Fransen
 
J
,
Johnson
 
SR
,
Baron
 
M
,
Tyndall
 
A
 et al.  
2013 classification criteria for systemic sclerosis: an American College of Rheumatology/European League against Rheumatism collaborative initiative
.
Arthritis Rheum
 
2013
;
65
:
2737
47
.

20

Macintyre
 
N
,
Crapo
 
R
,
Viegi
 
G
,
Johnson
 
D
,
Van der Grinten
 
C
,
Brusasco
 
V
 et al.  
Standardisation of the single-breath determination of carbon monoxide uptake in the lung
.
Eur Respir J
 
2005
;
26
:
720
35
.

21

Wanger
 
J
,
Clausen
 
J
,
Coates
 
A
,
Pedersen
 
O
,
Brusasco
 
V
,
Burgos
 
F
 et al.  
Standardisation of the measurement of lung volumes
.
Eur Respir J
 
2005
;
26
:
511
22
.

22

Galiè
 
N
,
Humbert
 
M
,
Vachiery
 
J-L
,
Gibbs
 
S
,
Lang
 
I
,
Torbicki
 
A
 et al.  
2015 ESC/ERS guidelines for the diagnosis and treatment of pulmonary hypertension: the joint task force for the diagnosis and treatment of pulmonary hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS): endorsed by: Association for European Paediatric and Congenital Cardiology (AEPC), International Society for Heart and Lung Transplantation (ISHLT)
.
Eur Heart J
 
2016
;
37
:
67
119
.

23

Lang
 
RM
,
Badano
 
LP
,
Mor-Avi
 
V
,
Afilalo
 
J
,
Armstrong
 
A
,
Ernande
 
L
 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
.
Eur Heart J Cardiovasc Imaging
 
2015
;
16
:
233
71
.

24

van Rosendael
 
AR
,
Dimitriu-Leen
 
AC
,
van Rosendael
 
PJ
,
Leung
 
M
,
Smit
 
JM
,
Saraste
 
A
 et al.  
Association between posterior left atrial adipose tissue mass and atrial fibrillation
.
Circ Arrhythm Electrophysiol
 
2017
;
10
:
e004614
.

25

Uysal
 
F
,
Akbal
 
E
,
Akbal
 
A
,
Cevizci
 
S
,
Arık
 
K
,
Gazi
 
E
.
Epicardial adipose tissue is increased in patients with inflammatory bowel disease
.
J Ultrasound Med
 
2016
;
35
:
1859
64
.

26

Alpaydın
 
S
,
Buyukterzi
 
Z
,
Akkurt
 
HE
,
Yılmaz
 
H
.
Impaired left ventricular diastolic functions and thickened epicardial adipose tissue in rheumatoid arthritis patients is correlated with DAS-28 score
.
Acta Cardiol Sin
 
2017
;
33
:
182
.

27

Lipson
 
A
,
Alexopoulos
 
N
,
Hartlage
 
GR
,
Arepalli
 
C
,
Oeser
 
A
,
Bian
 
A
 et al.  
Epicardial adipose tissue is increased in patients with systemic lupus erythematosus
.
Atherosclerosis
 
2012
;
223
:
389
93
.

28

Long
 
BD
,
Stojanovska
 
J
,
Brown
 
RK
,
Attili
 
AK
,
Jackson
 
EA
,
Ognenovski
 
V
.
Increased epicardial fat volume is independently associated with the presence and severity of systemic sclerosis
.
Acad Radiol
 
2017
;
24
:
1473
81
.

29

Fontes-Carvalho
 
R
,
Ladeiras-Lopes
 
R
,
Bettencourt
 
P
,
Leite-Moreira
 
A
,
Azevedo
 
A
.
Diastolic dysfunction in the diabetic continuum: association with insulin resistance, metabolic syndrome and type 2 diabetes
.
Cardiovasc Diabetol
 
2015
;
14
:
4
9
.

30

Figtree
 
GA
,
Rådholm
 
K
,
Barrett
 
TD
,
Perkovic
 
V
,
Mahaffey
 
KW
,
de Zeeuw
 
D
 et al.  
Effects of canagliflozin on heart failure outcomes associated with preserved and reduced ejection fraction in type 2 diabetes mellitus: results from the CANVAS Program
.
Circulation
 
2019
;
139
:
2591
3
.

31

Kurnicka
 
K
,
Domienik-Karłowicz
 
J
,
Lichodziejewska
 
B
,
Bielecki
 
M
,
Kozłowska
 
M
,
Goliszek
 
S
 et al.  
Improvement of left ventricular diastolic function and left heart morphology in young women with morbid obesity six months after bariatric surgery
.
Cardiol J
 
2018
;
25
:
97
105
.

32

Sorop
 
O
,
Olver
 
TD
,
van de Wouw
 
J
,
Heinonen
 
I
,
van Duin
 
RW
,
Duncker
 
DJ
 et al.  
The microcirculation: a key player in obesity-associated cardiovascular disease
.
Cardiovasc Res
 
2017
;
113
:
1035
45
.

33

Venteclef
 
N
,
Guglielmi
 
V
,
Balse
 
E
,
Gaborit
 
B
,
Cotillard
 
A
,
Atassi
 
F
 et al.  
Human epicardial adipose tissue induces fibrosis of the atrial myocardium through the secretion of adipo-fibrokines
.
Eur Heart J
 
2015
;
36
:
795
805
.

34

Alam
 
MS
,
Green
 
R
,
de Kemp
 
R
,
Beanlands
 
RS
,
Chow
 
BJ
.
Epicardial adipose tissue thickness as a predictor of impaired microvascular function in patients with non-obstructive coronary artery disease
.
J Nucl Cardiol
 
2013
;
20
:
804
12
.

35

Bakkum
 
M
,
Danad
 
I
,
Romijn
 
M
,
Stuijfzand
 
W
,
Leonora
 
R
,
Tulevski
 
I
 et al.  
The impact of obesity on the relationship between epicardial adipose tissue, left ventricular mass and coronary microvascular function
.
Eur J Nucl Med Mol Imaging
 
2015
;
42
:
1562
73
.

36

Wu
 
CK
,
Lee
 
JK
,
Hsu
 
JC
,
Su
 
MYM
,
Wu
 
YF
,
Lin
 
TT
 et al.  
Myocardial adipose deposition and the development of heart failure with preserved ejection fraction
.
Eur J Heart Fail
 
2020
;
22
:
445
54
.

37

Ng
 
AC
,
Strudwick
 
M
,
van der Geest
 
RJ
,
Ng
 
AC
,
Gillinder
 
L
,
Goo
 
SY
 et al.  
Impact of epicardial adipose tissue, left ventricular myocardial fat content, and interstitial fibrosis on myocardial contractile function
.
Circ Cardiovasc Imaging
 
2018
;
11
:
e007372
.

38

Abe
 
I
,
Teshima
 
Y
,
Kondo
 
H
,
Kaku
 
H
,
Kira
 
S
,
Ikebe
 
Y
 et al.  
Association of fibrotic remodeling and cytokines/chemokines content in epicardial adipose tissue with atrial myocardial fibrosis in patients with atrial fibrillation
.
Heart Rhythm
 
2018
;
15
:
1717
27
.

39

Şeker
 
T
,
Türkoğlu
 
C
,
Harbalıoğlu
 
H
,
Gür
 
M
.
Epicardial fat thickness is associated with abnormal left ventricle geometry in newly diagnosed hypertension
.
Acta Cardiol Sin
 
2018
;
34
:
280
.

40

Davis
 
JM
 III,
Roger
 
VL
,
Crowson
 
CS
,
Kremers
 
HM
,
Therneau
 
TM
,
Gabriel
 
SE
.
The presentation and outcome of heart failure in patients with rheumatoid arthritis differs from that in the general population
.
Arthritis Rheum
 
2008
;
58
:
2603
11
.

41

Davis
 
JM
 III,
Lin
 
G
,
Oh
 
JK
,
Crowson
 
CS
,
Achenbach
 
SJ
,
Therneau
 
TM
 et al.  
Five-year changes in cardiac structure and function in patients with rheumatoid arthritis compared with the general population
.
Int J Cardiol
 
2017
;
240
:
379
85
.

42

Kim
 
CH
,
Al-Kindi
 
SG
,
Jandali
 
B
,
Askari
 
AD
,
Zacharias
 
M
,
Oliveira
 
GH
.
Incidence and risk of heart failure in systemic lupus erythematosus
.
Heart
 
2017
;
103
:
227
33
.

43

Elnady
 
BM
,
Abdelghafar
 
ASM
,
El Shazly Abdul Khalik
 
MM
,
Algethami
 
A
,
Al-Otaibi
 
MDA
,
Al-Otaibi
 
ME
.
The implication of tissue Doppler echocardiography and cardiopulmonary exercise in early detection of cardiac dysfunction in systemic lupus erythematosus patients
.
Eur J Rheumatol
 
2016
;
3
:
109
17
.

Author notes

Conflict of interest: The Department of Cardiology of the Leiden University Medical Center received unrestricted research grants from Abbott Vascular, Alnylam, Bayer, Bioventrix, Biotronik, Boston Scientific, Edwards Lifesciences, GE Healthcare, and Medtronic. V.D. received speaker fees from Abbott Vascular, Edwards Lifesciences, GE Healthcare, MSD, Novartis, and Medtronic. J.J.B. received speaker fees from Abbot Vascular and Edwards Lifesciences. N.A.M. received speaker fees from GE Healthcare, Philips Ultrasound, and Abbott Vascular, and research grant from Alnylam. J.K.D.V.-B. received consulting fees from Abbvie, Janssen, and Boehringer Ingelheim, and received research grants from ReumaNederland, NVLE, Roche, Galapagos, and Janssen. All payments were made to her institution. The remaining authors have nothing to disclose.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Supplementary data