1Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 ORIGINAL RESEARCH Sex-­specific­difference­in­cardiac­ function­in­patients­with­systemic­ sclerosis:­association­with­ cardiovascular­outcomes Tea Gegenava ,1,2 Federico Fortuni,1,3 Nina Marijn van Leeuwen ,4 Anders H Tennoe,5 Anna- Maria Hoffmann- Vold ,5 Ruxandra Jurcut,6 Adrian Giuca,6 Laura Groseanu,7,8 Felix Tanner,9 Oliver Distler ,10 Jeroen J Bax,1,11 Jeska De Vries- Bouwstra,4 Nina Ajmone Marsan1 To cite: Gegenava T, Fortuni F, van Leeuwen NM, et al. Sex- specific difference in cardiac function in patients with systemic sclerosis: association with cardiovascular outcomes. RMD Open 2023;9:e003380. doi:10.1136/ rmdopen-2023-003380 ► Additional supplemental material is published online only. To view, please visit the journal online (http:// dx. doi. org/ 10. 1136/ rmdopen- 2023- 003380). Received 8 June 2023 Accepted 4 September 2023 For numbered affiliations see end of article. Correspondence to Dr Nina Ajmone Marsan; n. ajmone@ lumc. nl Systemic sclerosis © Author(s) (or their employer(s)) 2023. Re- use permitted under CC BY- NC. No commercial re- use. See rights and permissions. Published by BMJ. ABSTRACT Background Cardiovascular involvement is one of the leading causes of mortality in systemic sclerosis (SSc) and is reported to be higher in men as compared with women. However, the cause of this difference is largely unknown. The objective of this study was to assess sex differences in echocardiographic characteristics, including left ventricular global longitudinal strain (LV GLS), as a potential explanation of sex differences in outcomes. Methods A total of 746 patients with SSc from four centres, including 628 (84%, 54±13 years) women and 118 (16%, 55±15 years) men, were evaluated with standard and advanced echocardiographic examinations. The independent association of the echocardiographic parameters with the combined endpoint of cardiovascular events- hospitalisation/death was evaluated. Results Men and women with SSc showed significant differences in disease characteristics and cardiac function. After adjusting for the most important clinical characteristics, while LV ejection fraction and diastolic function were not significantly different anymore, men still presented with more impaired LV GLS as compared with women (−19% (IQR −20% to −17%) vs −21% (IQR: −22% to −19%), p<0.001). After a median follow- up of 48 months (IQR: 26–80), the combined endpoint occurred in 182 patients. Men with SSc experienced higher cumulative rates of cardiovascular events- hospitalisation/mortality (χ2=8.648; Log- rank=0.003), and sex differences were maintained after adjusting for clinical confounders, but neutralised when matching the groups for LV GLS. Conclusion In patients with SSc, male sex is associated with worse cardiovascular outcomes even after adjusting for important clinical characteristics. LV GLS was more impaired in men as compared with women and potentially explains the sex difference in cardiovascular outcomes. INTRODUCTION Systemic sclerosis (SSc) is a chronic auto- immune disease characterised by progres- sive fibrosis, microvascular changes and multiorgan involvement. According to previ- ously published data, 7%–40% of patients with SSc have cardiovascular involvement,1–4 which also showed an association with poor WHAT IS ALREADY KNOWN ON THIS TOPIC ⇒ Cardiovascular involvement is highly prevalent in patients with systemic sclerosis (SSc) and is asso- ciated with poor prognosis. Sex has been also asso- ciated with the prognosis in patients with SSc, with men being at higher risk of death and cardiovascular complications. WHAT THIS STUDY ADDS ⇒ Men and women with SSc showed differences in terms of disease characteristics and cardiac func- tion. Left ventricular global longitudinal strain (LV GLS) was more impaired in men as compared with women. After adjusting for clinical characteristics, while LV ejection fraction and diastolic function were not significantly different, men still presented with more impaired LV GLS as compared with women. Men with SSc experienced higher cumulative rates of cardiovascular outcomes. Sex difference in out- comes was maintained after adjusting for clinical characteristics, but neutralised when matching the groups for LV GLS. HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY ⇒ In patients with SSc, early assessment of subtle myocardial changes and cardiac involvement by LV GLS may improve risk stratification and define the need for further diagnostic assessment, closer follow- up and eventually specific treatment in both men and women; in fact, the lack of interaction be- tween sex and LV GLS suggested an equally high prognostic role of LV GLS for both sexes among pa- tients with SSc. However, taking into account sex- specific clinical and echocardiographic parameters may help to better define factors associated with outcomes. Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 2 Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 RMD Open prognosis.5 Sex has also been associated with adverse events in patients with SSc, men being at higher risk of death and cardiovascular complications.1 6 Echocardiography is the first- line imaging technique to identify left ventricular (LV) systolic or diastolic dysfunc- tion in these patients; however, diagnosis of cardiac involvement remains challenging, especially at an early stage, as conventional echocardiographic measures lack sensitivity to detect subtle myocardial dysfunction. Previous studies have demonstrated that 2D speckle tracking strain echocardiography (STE) is a reliable and sensitive tool for detecting subtle LV dysfunction in various cardiovascular diseases,7 8 and LV global longitu- dinal strain (LV GLS) has been introduced in the most recent recommendations for chamber quantification.9 In patients with SSc, LV GLS has shown an association with lower functional capacity, occurrence of arrhythmias7 and all- cause mortality and cardiovascular events.10 However, sex differences in standard and novel echocardiographic measures have never been explored in these patients as well as the impact of these differences on prognosis and cardiovascular risk profile. The aim of this study was therefore to identify poten- tial differences in clinical and echocardiographic param- eters, including LV GLS, between men and women in a large, well- characterised multicentre cohort of patients with SSc and to test their potential association with all- cause mortality and cardiovascular events. METHODS Patient population Consecutive patients with SSc referred for a specifi- cally designed multidisciplinary healthcare programme at the Department of Rheumatology of the Leiden University Medical Centre (The Netherlands), Univer- sity Hospital of Zurich (Switzerland), Oslo University Hospital (Norway) and the University of Medicine and Pharmacy ‘Carol Davila’, Bucharest (Romania) were included.11 12 Patients were diagnosed according to the classification criteria American College of Rheumatology (ACR) 2013.13 A baseline echocardiogram, defined as the first one available, was selected and corresponding visit at the abovementioned departments was defined as a baseline visit. Only patients with a minimum of 1 year follow- up were included in this analysis. Furthermore, to avoid the potential confounding effect on the assess- ment of cardiac function, patients with previous myocar- dial infarction or severe valvular heart disease were excluded. Written informed consent was obtained at the time of inclusion in the multidisciplinary healthcare programmes. Additional informed consent was acquired when patients were contacted by telephone to evaluate the occurrence of the pre- specified endpoints. Baseline clinical variables Disease- related characteristics at baseline included SSc subtype according to LeRoy et al,14 modified Rodnan Skin Score (mRSS), presence of digital ulcers, arthritis, proximal muscle weakness, myositis, calcinosis and pitting scars.15 General medical history, medications and cardiovascular risk factors were also noted. Labora- tory testing was performed systematically and included creatine phosphokinase (CK), erythrocyte sedimenta- tion rate, renal function (estimated glomerular filtra- tion rate (eGFR)) and N- terminal pro- brain natriuretic peptide (NT- proBNP). Spirometry was performed in all patients according to the American Thoracic Society and the European Respiratory Society guidelines,16 17 and included the percentage of predicted values for forced vital capacity (FVC), forced expiratory volume in 1 s (FEV1) and single breath diffusing capacity of carbon monoxide (DLCO- SB). The presence of lung involve- ment was assessed using thoracic high- resolution CT (HRCT) and defined by abnormal findings including interstitial fibrosis in patients with available HRCT images.18 19 To assess the percentage of predicted peak oxygen consumption (VO2max), cardiopulmonary exer- cise testing was performed using an electrically ramped cycle ergometer according to the guidelines.20 Finally, (24- hour) Holter ECG was performed. The presence of a bundle branch block, ventricular arrhythmias including premature ventricular contraction >100/day and sustained ventricular tachycardias, and supraventricular arrhythmias including atrial fibrillation, atrial flutter and premature atrial contraction when >7% were recorded and considered abnormal.21 Conventional transthoracic echocardiography Transthoracic echocardiography was performed in the left lateral decubitus position, using a commercially avail- able system (Vivid 7 and E9; General Electric- Vingmed, Horten, Norway) and a 3.5- MHz or 5MS transducer. Standard M- mode and 2D, colour, pulsed wave and continuous wave Doppler images were acquired and digitally stored in cine- loop format. Off- line analysis was performed using EchoPAC (version 112.0.1; GE Medical Systems, Horten, Norway). LV end- diastolic volume (LVEDV), LV end- systolic volume (LVESV), LV ejection fraction (LVEF) and LV dimensions, including end- diastolic diameter, interventricular septum thickness and posterior wall thickness were measured according to current guidelines.9 LV diastolic function was assessed according to current recommendations based on a multiparametric approach and included: peak early (E) and late (A) diastolic velocities measured on pulsed wave Doppler recordings of the trans- mitral flow,22 Left Atrial Volume Index measured before mitral valve opening on the apical four- chamber and two- chamber views according to the biplane Simpson’s method and indexed to the patient’s body surface area. Additionally, e′ with tissue Doppler imaging at the lateral and septal sides of the mitral annulus in the four- chamber view was meas- ured and averaged; E- wave/e′ ratio was then calculated.22 Systolic pulmonary arterial pressure was estimated by the right ventricular systolic pressure which was calculated by Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 3Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 Systemic sclerosis adding the peak tricuspid regurgitation gradient to the right atrial pressure which was estimated by the inferior vena cava diameter and degree of respiratory collapse.23 Right ventricular function was evaluated by measuring the tricuspid annular plane systolic excursion.24 The presence of pericardial effusion was also recorded. Two-dimensional STE LV GLS was measured as previously described by using the four- chamber, two- chamber and three- chamber views to obtain measurements of all LV segments.10 25 26 LV GLS was automatically calculated as the average peak systolic strain of 17 LV segments (figure 1). Negative values are used to express LV GLS since this represents the short- ening of the myocardial wall, a more negative value represents better myocardial deformation. Follow-up and endpoint The endpoint of this study was defined as all- cause mortality or hospitalisations for cardiovascular reasons which included pacemaker/cardioverter defibrillator implantations, heart failure, myocardial infarction, arrhythmias, angina pectoris, percutaneous coronary interventions, stroke and peripheral ischaemic disease. Endpoint data were obtained by reviewing the electronic information system and retrieval of survival status through the municipal civil registries (according to availability). If the last available follow- up data were dated more than 1 year ago or incomplete, patients and their general prac- titioners were contacted by telephone for more recent information (in selected centres). Statistical analysis Data analysis was performed using SPSS software V.23.0 (IBM, Armonk, New York, USA). Continuous variables are presented as mean±SD or as median with IQRs according to presence of normal distribution. Categor- ical data are presented as frequencies and percentages. Continuous variables were compared using one- way analysis of variance (ANOVA), applying the Bonfer- roni’s post- hoc analysis, or the Kruskal- Wallis one- way ANOVA. Categorical variables were compared with the χ2 test. Propensity matching score was used to match men and women groups based on age, disease duration (since Raynaud), subtype of SSc (diffuse cutaneous SSc (dcSSc)), presence of interstitial lung fibrosis, DLCO- SB and NT- proBNP, and in a second step also adding LV GLS (these variables were chosen as signifi- cantly different among men and women, and also based on clinical relevance). The matching tolerance was set at 0.05. Cumulative event rates for all- cause mortality or cardiovascular hospitalisation were estimated with the Kaplan- Meier curves and compared with log- rank tests. P values <0.05 were considered significant. The associ- ation between LV GLS and cardiovascular outcomes, as well as interaction between LV GLS and sex was evalu- ated with Cox regression analyses. To assess the incre- mental value of LV GLS, we compared χ2 of different multivariate cox- regression models which included the variables selected for the propensity score matching and other variables which were significantly different between men and women. RESULTS Baseline clinical characteristics A total of 746 patients with SSc from four different centres were included, of which 628 (84%, 55±14 years) were women and 118 (16%, 54±13 years) were men. Base- line clinical characteristics are shown in table 1. Men and women patients with SSc showed several differences in terms of disease characteristics: greater mRSS, higher prevalence of dcSSc than limited cutaneous SSc (lcSSc), myositis and interstitial lung disease (ILD), higher CK and higher eGFR were observed in men as compared with women (table 1). In addition, disease duration was significantly longer in women as compared with men and Figure 1 Example of assessment of LV GLS by speckle tracking strain echocardiography in a woman with SSc (panel A) versus a man with SSc (panel B). Color- coded Bull’s eye plots and the curves for longitudinal strain (panel A, LV GLS=−17%, panel B, LV GLS=−12%) are displayed; the strain curves are showed per segment and averaged among the segments (dotted line for the three- chamber, four- chamber and two- chamber apical views). LV GLS, left ventricular global longitudinal strain; SSc, systemic sclerosis. Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 4 Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 RMD Open Table 1 Baseline clinical characteristics in the overall population and divided by men and women Baseline clinical characteristics Total (n=746) Men (n=118) Women (n=628) P value Age (years), mean±SD 55±14 54±13 55±14 0.351 Body surface area (m2), mean±SD 1.77±0.19 1.95±0.19 1.73±0.17 <0.001 Diffuse cutaneous SSc, n (%) 192 (30) 52 (48) 140 (26) <0.001 Time since Raynaud (years), median (IQR)- till first echo 10 (4–19) 7 (3–14) 10 (4–19) 0.006 Time since non- Raynaud (years), median (IQR)- till first echo 5 (2–10) 4 (1–9) 5 (2–11) 0.059 Disease characteristics, n (%) mRSS >15 94 (16) 26 (26) 68 (13) 0.002 Digital ulcers 226 (32) 39 (35) 187 (32) 0.267 Pitting scars 227 (38) 45 (45) 182 (36) 0.060 Proximal muscle weakness 39 (9) 8 (12) 31 (9) 0.240 Myositis 12 (3) 7 (13) 5 (2) <0.001 Calcinosis 166 (28) 18 (19) 148 (29) 0.097 Lung fibrosis 332 (49) 67 (63) 265 (46) 0.002 Medical history and cardiovascular risk factors, n (%) Hypertension 151 (21) 24 (21) 127 (21) 0.537 Diabetes 16 (4) 4 (6) 12 (3) 0.234 (History of) smoking 346 (50) 76 (65) 270 (47) <0.001 Coronary artery disease 67 (9) 14 (12) 53 (9) 0.172 Prior/current pericarditis 55 (7) 9 (8) 46 (7) 0.513 Holter ECG abnormalities (prior or current arrhythmias) 125 (31) 25 (40) 100 (29) 0.052 History of renal crisis 15 (4) 4 (6) 11 (3) 0.188 Laboratory tests Creatine phosphokinase (U/L), median (IQR) 84 (61–122) 115 (79–200) 79 (58–111) <0.001 eGFR (mL/min/1.73 m2), mean±SD 89±27 103±31 87±26 <0.001 NT- proBNP (ng/L), median (IQR) 60 (19–143) 37 (11–98) 62 (22–150) 0.005 ANA/ANF (n=397) 380 (96) 55 (89) 325 (97) 0.009 Pulmonary function tests, mean±SD FVC % of predicted 96±22 86±20 98±22 <0.001 FEV1 % of predicted 90±20 82±20 91±20 <0.001 DLCO- SB % of predicted 66±21 60±22 67±20 0.001 VO2max % of predicted, n=284 87±24 72±23 90±23 <0.001 Mean heart rate, beats per minute 79±10 80±11 79±10 0.242 Cardiovascular medication, n (%) ACEi/ARB 268 (37) 56 (48) 212 (35) 0.004 Beta- blocker 126 (18) 21 (18) 105 (17) 0.480 Ca2+ channel blocker 431 (60) 77 (66) 354 (59) 0.077 Diuretics 184 (26) 30 (26) 154 (26) 0.520 Immunosuppressive medication, n (%) Corticosteroids 63 (15) 16 (25) 47 (14) 0.020 Cyclophosphamide 10 (2) 1 (1) 9 (3) 0.511 Methotrexate, n=414 54 (13) 7 (11) 47 (14) 0.627 Azathioprine, n=414 21 (5) 6 (9) 15 (4) 0.097 p value <0.05 is considered significant. ACEi, ACE inhibitor; ANA, antinuclear antibodies; ANF, antinuclear factor; ARB, angiotensin II receptor antagonist; DLCO- SB, diffusing capacity for carbon monoxide single breath; eGFR, estimated glomerular filtration rate; FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity; mRSS, modified Rodnan skin score; NT- proBNP, N- terminal pro- brain natriuretic peptide; SSc, systemic sclerosis; VO2max, maximal oxygen uptake. Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 5Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 Systemic sclerosis positivity for antinuclear antibodies/antinuclear factor was higher in women (table 1). Women and men were comparable in terms of cardio- vascular disease risk factors, although smoking was more prominent in men. Pulmonary functional tests (FVC, FEV, DLCO- SB, VO2max) were more impaired in men and NT- proBNP was higher in women. ACE/angiotensin II receptor antagonists and corticosteroids were more frequently prescribed in men (table 1). Baseline echocardiographic characteristics The baseline echocardiographic characteristics of the overall population are presented in table 2. When comparing standard echocardiographic measures between men and women patients with SSc (particu- larly LV systolic function), we observed that men were characterised by larger LV indexed volumes (LVEDV: 54 (44–62 mL/m2) vs 47 (40–55 mL/m2), p<0.001; LVESV: 21 (17–25 mL/m2) vs 18 (14–22 mL/m2), p<0.001) and lower LVEF (men SSc 59±8% vs women SSc 60±7%, p=0.035). By the STE analysis, LV GLS was more preserved in women (−21% (IQR: −22% to −19%) as compared with men (−19% (IQR −20% to −17%), p<0.001); in turn, there were no significant differences between men and women in terms of LV diastolic func- tion parameters (table 2). Considering the significant difference between men and women for important clinical characteristics, a propensity matching score was applied to explore whether the sex differences in the echocardiographic parameters were maintained also after adjusting for these confounders. The matching was performed according to age, disease duration since Raynaud, presence of dcSSc, interstitial lung fibrosis, DLCO- SB and NT- proBNP; after matching (n=140 patients), LV GLS still showed a significant difference between men and women, whereas LV volumes and LVEF did not (table 3). Sex differences in outcomes in patients with SSc After a median follow- up of 48 months (IQR: 26–80), the combined endpoint occurred in 182 (24.4%) patients. Kaplan- Meier survival curves showed that men experi- enced higher cumulative rates of cardiovascular events- hospitalisation/death as compared with women (χ2 8.648; Log rank 0.003, figure 2). When using propensity score matching to match men and women according to the abovementioned clinical characteristics (age, disease dura- tion since Raynaud, type of SSc, lung fibrosis, DLCO- SB and NT- proBNP, n=140 patients), men still experienced higher cumulative rates of cardiovascular events/death as compared with women (χ2 7.211; Log rank 0.007, figure 3A). Sex difference in outcome was neutralised when matching the groups according to the LV GLS on top of the abovementioned clinical characteristics (n=112 patients, χ2 0.474; Log rank 0.491, figure 3B). Table 2 Baseline echocardiographic characteristics Total (n=746) Men (n=118) Women (n=628) P value LVEDV (mL), median (IQR) 85 (69–101) 104 (89–118) 82 (57–96) <0.001 LVEDV/BSA (mL/m2) 48 (40–56) 54 (44–62) 47 (40–55) <0.001 LVESV (m), median (IQR) 32 (25–41) 41 (34–50) 31 (25–39) <0.001 LVESV/BSA (mL/m2) 19 (15–23) 21(17–25) 18 (14–22) <0.001 LVEF %, mean±SD 60±7 59±8 60±7 0.035 IVST (mm), mean±SD 9±2 10±2 9±2 <0.001 LVEDd (mm), mean±SD 48±5 51±5 47±5 <0.001 PWT (mm), mean±SD 8±2 10±2 9±2 <0.001 LAVI (mL), median (IQR) 31 (25–38) 32 (25–38) 30 (25–37) 0.502 E/A ratio, median (IQR) 1 (0.9–1.3) 1 (0.9–1.00) 1 (0.9–1.00) 0.852 e′ (cm/s), mean±SD 9±3 9±3 9±3 0.912 E/e′ ratio, median (IQR) 9 (7–11) 8 (6–11) 9 (7–11) 0.865 sPAP (mm Hg), mean±SD 32±17 32±14 32±17 0.988 sPAP >35 mm Hg, n (%) 97 (23) 16(29) 81(22) 0.219 TAPSE (mm), mean±SD 23±4 23±4 22±4 0.163 Pericardial effusion, n (%) 55 (7) 11 (9) 44 (7) 0.379 LV GLS %, median (IQR)* −20 (−22 to −18) −19 (−20 to −17) −21 (−22 to −19) 0.001 p value <0.05 is considered significant. BSA, body surface area; IVST, interventricular septum thickness; LAVI, Left Atrial Volume Index; LVEDd, left ventricular end- diastolic dimension; LVEDV, left ventricular end- diastolic volume; LVEF, left ventricular ejection fraction; LVESV, left ventricular end- systolic volume; LV GLS, left ventricular global longitudinal strain; PWT, posterior wall thickness; sPAP, systolic pulmonary artery pressure; TAPSE, tricuspid annular plane systolic excursion. Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 6 Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 RMD Open Table 3 Echocardiographic findings in men and women patients with SSc after propensity score matching (based on age, disease duration since onset of Raynaud, SSc type, interstitial lung fibrosis, DLCO- SB, NT- proBNP) Baseline echocardiographic characteristics Propensity matched patients Total (n=140) Men (n=70) Women (n=70) P value LVEDV/BSA (mL/m2) 49 (41–57) 53 (42–60) 46 (40–55) 0.059 LVESV/BSA (mL/m2) 18 (15–23) 19 (17–24) 18 (14–23) 0.059 LVEF %, mean±SD 60±7 60±8 61±6 0.226 LAVI (mL), median (IQR) 27 (22–34) 26 (20–33) 27 (23–35) 0.583 E/A ratio, median (IQR) 1 (0.9–1.3) 1 (0.9–1.3) 1 (0.9–1.2) 0.812 E/e′ ratio, median (IQR) 8 (6–10) 7 (5–8) 9 (6–10) 0.132 sPAP (mm Hg), mean±SD 34±17 33±15 35±20 0.656 TAPSE (mm), mean±SD 23±4 24±4 22±3 0.059 LV GLS %, median (IQR) −20 (−21 to −18) −19 (−20 to −18) −20 (−22 to −18) 0.003 p value <0.05 is considered significant. BSA, body surface area; DLCO- SB, diffusing capacity for carbon monoxide single breath; LAVI, Left Atrial Volume Index; LVEDV, left ventricular end- diastolic volume; LVEF, left ventricular ejection fraction; LVESV, left ventricular end- systolic volume; LV GLS, left ventricular global longitudinal strain; NT- proBNP, N- terminal pro- brain natriuretic peptide; sPAP, systolic pulmonary artery pressure; SSc, systemic sclerosis; TAPSE, tricuspid annular plane systolic excursion. Figure 2 Survival function (for cardiovascular hospitalisation/death) in women and men patients with systemic sclerosis (SSc). Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 7Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 Systemic sclerosis In addition, LV GLS showed an association with outcomes in the overall group (HR: 1.163; 95% CI: 1.103 to 1.226, p<0.001) and no interaction between sex and LV GLS was detected (HR: 0.936; 95% CI: 0.811 to 1.082, p=0.373), indicating an equally high prognostic role of LV GLS for both sex in the SSc population. We also performed an additional Cox regression analysis in the overall populations including the variables of the propensity score matching and sex. When adding LV GLS to a multivariable model including all these variables (basal model), it significantly increased the predictivity of the model; in turn, adding any of the other variables which were different between men and women (but still not included in the propensity matching analysis) did not increase the predictivity of the model significantly, confirming the importance of GLS to predict prognosis and possibly its role in the higher cardiovascular risk of men (online supplemental figures 1a–e). DISCUSSION The main findings of the present study, which included a large cohort of patients with SSc from multiple centres, can be summarised as follows: (1) after matching for the various clinical and SSc- related characteristics, men with SSc presented with more impaired LV GLS as compared with women, but no differences in LVEF or diastolic function parameters were observed and (2) LV GLS was significantly associated with cardiovascular outcome both in men and women but the sex difference in outcome was no longer observed when matching the groups according to LV GLS on top of the clinical characteristics. Differences in LV function between men and women with SSc In patients with SSc, cardiovascular system involvement is observed both in men and women, but previous studies have shown that LV dysfunction (as defined by an impaired LVEF) is more prevalent in men. Similarly in our study, men showed lower LVEF and, as shown for the first time, also more impaired LV GLS as compared with women. However, it has not been demonstrated so far whether this sex difference in cardiac involvement remains significant after adjusting for important char- acteristics which have been shown to be associated with LV dysfunction, such as age, disease duration and type of disease, digital ulcers, muscle and lung involvement.5 Like in other cohorts, also in our study men patients with SSc had more severe disease with more frequent dcSSc, SSc- associated ILD and more impaired DLCO. In general, higher burden of organ fibrosis formation may explain the higher prevalence of LV dysfunction in men.27–29 Also the dcSSc subtype has been previously associated with cardiovascular involvement.30 31 In the study by Ferri et al,2 patients with dcSSc presented more frequently with cardiac involvement at the time of diag- nosis as compared with patients with lcSSc. An increased risk of cardiac disease in dcSSc was also reported in the German Network for Systemic Scleroderma cohort.30 Sex differences in LV dysfunction should be therefore corrected for these important factors, and in our study a more impaired LV GLS (but not LVEF) was confirmed in men as compared with women even after adjusting for the presence of dcSSc and ILD, both reflecting more severe fibrotic disease. Figure 3 Survival function in women and men patients with SSc in panel A after adjusting for age, disease duration (since Raynaud), type of SSc, lung fibrosis, DLCO- SB and NT- proBNP, in panel B after adjusting for age, disease duration (since Raynaud), type of SSc, lung fibrosis, DLCO- SB, NT- proBNP and LV GLS. DLCO- SB, diffusing capacity for carbon monoxide single breath; LV GLS, left ventricular global longitudinal strain; NT- proBNP, N- terminal pro- brain natriuretic peptide; SSc, systemic sclerosis. Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 8 Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 RMD Open In our cohort, no significant sex differences were observed in digital ulcers or proximal muscle weak- ness; furthermore, NT- proBNP levels were higher in women as compared with men. In the study by Lau et al, which extensively explored sex differences in circulating biomarkers among different cardiovascular diseases, NT- proBNP showed typically higher values in women.32 This might explain the observation in our patient cohort (higher NT- proBNP level in women) as well, although LVEF values were still relatively preserved. However, also after correction for NT- proBNP values, the difference in LV GLS between men and women remained significant. Clinical and echocardiographic associates of cardiovascular outcomes in men and women with SSc Cardiovascular events are one of the major causes of mortality in patients with SSc.1 33 Although several studies have shown that male sex is associated with worse cardi- ovascular outcome/mortality in general in SSc,5 6 it has never been shown whether this difference was related with the prevalence of LV dysfunction, particularly when assessed with advanced echocardiographic meas- ures such as LV GLS. In SSc, LV GLS has emerged as an important marker of subtle myocardial involvement (due to myocardial fibrosis related to microangiopathy and microvascular injury, myocardial stunning due to repeated focal ischemia caused by abnormal vasoreac- tivity or peri- myocarditis) and dysfunction, even when LVEF is preserved.5 34 In the study by van Wijngaarden et al, LV GLS showed to be independently associated with cardiovascular outcomes in a large cohort of patient with SSc together with NT- proBNP and DLCO.10 In the current study, including a large cohort of patients from multiple centres, we confirmed that men were characterised by a higher incidence of cardiovas- cular events/mortality and LV GLS was associated with the outcome both in men and women. Interestingly, sex differences in outcomes were maintained also after correcting for important clinical characteristics but were neutralised after matching for LV GLS, suggesting there- fore an important role of LV GLS in the higher incidence of cardiovascular events/mortality in men. GLS might be the underlying problem by which the survival difference between men and women is caused. Additional studies aiming at demonstrating the underlying pathophysi- ology of LV GLS impairment in SSc are therefore highly warranted as well as characterisation of sex differences and potentially changes over time. CLINICAL IMPLICATIONS In patients with SSc, the presence of impaired LVEF and cardiovascular symptoms is of important prognostic value but identifies an advanced stage of the cardiac involve- ment, whereas risk stratification should help to identify patients at higher risk prior to developing advanced disease for which earlier treatment can still potentially change the natural course and prognosis of the disease. Early assessment of subtle myocardial changes and cardiac involvement by LV GLS may improve risk strat- ification and define the need for further diagnostic assessment (eg, cardiac MRI, right heart catheterisation, Holter ECG monitoring), closer follow- up and eventu- ally specific treatment in both men and women; in fact, lack of interaction between sex and LV GLS suggested an equally high prognostic role of LV GLS for both sexes among patients with SSc. However, considering sex- specific clinical and echocardiographic parameters may help to better define factors associated with outcomes of the disease and therefore risk- stratifying and managing these patients. In particular, normal values of LV GLS in men patients with SSc may reassure the treating physi- cian of the absence of cardiac involvement in this group considered at high risk for cardiovascular complications. LIMITATIONS The limitations of this study include: (1) specific anal- yses for myocardial ischemia were not systematically performed; however, most of the patients underwent cardiopulmonary exercise testing without showing signs of myocardial ischemia making therefore its presence unlikely; (2) treatment strategies were not taken into account during follow- up and potential association with cardiac function could not be evaluated; however, consid- ering the absence of homogeneous recommendations, therapeutic changes were also not included in previously published prognostic models; (3) the relatively small number of men patients with SSc still gives some statis- tical power limitation but propensity score matching was used for group comparison; (4) echocardiograms were not read in a core lab but were directly entered in the database in each centre that, nevertheless, has great experience in analysing standard and advanced echocar- diography and (5) significant differences between men and women for important clinical characteristics were observed, and not all of those variables were included in the propensity score matching; however, we have also built different models to further test the incremental prognostic value of LV GLS on top of those variables observed as different between men and women. CONCLUSIONS In SSc, male sex is associated with worse cardiovascular outcomes even after adjusting for important clinical char- acteristics. LV GLS was more impaired in men patients with SSc, as compared with women and this may poten- tially explain the sex differences in cardiovascular events. Using LV GLS for early assessment of myocardial involve- ment for both sexes may improve risk stratification and surveillance in patients with SSc. Further research that aims at elucidating the cause of LV GLS impairment in SSc is highly recommended, as this might lead to improved and targeted treatment strategies. Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 9Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 Systemic sclerosis Author affiliations 1Department of Cardiology, Leiden University Medical Center, Leiden, The Netherlands 2Internal Medicine, Tbilisi State Medical University, Tbilisi, Georgia 3Department of Cardiology, Ospedale Nuovo San Giovanni Battista, Foligno, Umbria, Italy 4Department of Rheumatology, Leiden University Medical Center, Leiden, The Netherlands 5Department of Rheumatology, Oslo University Hospital, Oslo, Norway 6Department of Cardiology, Emergency Institute for Cardiovascular Diseases Prof C C Iliescu, Bucuresti, Romania 7Department of Internal Medicine Rheumatology, University of Medicine and Pharmacy "Carol Davila", Bucharest, Romania 8Internal Medicine and Rheumatology, Santa Maria Clinical Hospital, Bucharest, Romania 9Department of Cardiology, University Hospital Zurich, Zurich, Switzerland 10Department of Rheumatology, University Hospital Zurich, Zurich, Switzerland 11Department of Cardiology, Turku University Finland, Turku, Finland Twitter Tea Gegenava @TeaGegenava Acknowledgements Abstract entitled “Sex- specific difference in cardiac function in patients with systemic sclerosis: association with cardiovascular outcomes” has been presented at European Society of Cardiology Congress (ESC) 2021. Contributors TG and NAM contributed to conception and design of the study; analysis and interpretation of data; drafting of the manuscript and final approval of the manuscript submitted. TG and NAM accept full responsibility for the work and/ or the conduct of the study, had access to the data, and controlled the decision to publish. NAM is a guarantor. FF and NMvL contributed to analysis and interpretation of data; drafting of the manuscript and final approval of the manuscript submitted. AHT, RJ and AG contributed to interpretation of data; drafting of the manuscript and final approval of the manuscript submitted. A- MH- V, FT, OD, JJB, LG and JDV- B contributed to critical revision of the manuscript for important intellectual content and final approval of the manuscript submitted. Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not- for- profit sectors. Competing interests 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. JJB received speaker fees from Abbot Vascular and Edwards Lifesciences. NAM received speaker fees from GE Healthcare, Philips Ultrasound and Abbott Vascular; also research grant from Alnylam and Pfizer. JDV- B received consulting fees from Abbvie, Janssen and Boehringer Ingelheim, and received research grants from Roche, Galapagos and Janssen; all payments were made to her institution. AHT has received consulting fees or other remuneration from GlaxoSmithKline (GSK) and Actelion. A- MH- V has received unrestricted research grants from Boehringer Ingelheim, research funding from Janssen, consultancy and speaker fees from Arxx, Boehringer Ingelheim, Janssen, Lilly, Medscape, Merck Sharp & Dohme and Roche. OD has/had consultancy relationship with and/or has received research funding from and/or has served as a speaker for the following companies in the area of potential treatments for systemic sclerosis and its complications in the last three calendar years: 4P- Pharma, Abbvie, Acceleron, Alcimed, Altavant, Amgen, AnaMar, Arxx, AstraZeneca, Blade, Bayer, Boehringer Ingelheim, Corbus, CSL Behring, Galderma, Galapagos, Glenmark, Gossamer, Horizon, Janssen, Kymera, Lupin, Medscape, Merck, Miltenyi Biotec, Mitsubishi Tanabe, Novartis, Prometheus, Redxpharma, Roivant and Topadur. Patent issued “mir- 29 for the treatment of systemic sclerosis” (US8247389, EP2331143). Patient consent for publication Consent obtained directly from patient(s). Ethics approval This study involves human participants and was approved by the METC institutional and other local ethical committees (numbers of the protocols approved by METC P09.003/SH/sh; REU 036/SH/sh; REU 043/SH/shB; 19.008/KB/ kb). Participants gave informed consent to participate in the study before taking part. Provenance and peer review Not commissioned; externally peer reviewed. Data availability statement Data are available upon reasonable request. Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer- reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise. Open access This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY- NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non- commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non- commercial. See: http://creativecommons.org/licenses/by-nc/4.0/. ORCID iDs Tea Gegenava http://orcid.org/0000-0003-2560-169X Nina Marijn van Leeuwen http://orcid.org/0000-0003-4228-012X Anna- Maria Hoffmann- Vold http://orcid.org/0000-0001-6467-7422 Oliver Distler http://orcid.org/0000-0002-0546-8310 REFERENCES 1 Tyndall AJ, Bannert B, Vonk M, et al. Causes and risk factors for death in systemic sclerosis: a study from the EULAR scleroderma trials and research (EUSTAR) database. Ann Rheum Dis 2010;69:1809–15. 2 Ferri C, Valentini G, Cozzi F, et al. Systemic sclerosis: demographic, clinical, and serologic features and survival in 1,012 Italian patients. Medicine (Baltimore) 2002;81:139–53. 3 Ngian G- S, Sahhar J, Proudman SM, et al. Prevalence of coronary heart disease and cardiovascular risk factors in a national cross- sectional cohort study of systemic sclerosis. Ann Rheum Dis 2012;71:1980–3. 4 Tennøe AH, Murbræch K, Andreassen JC, et al. Left ventricular diastolic dysfunction predicts mortality in patients with systemic sclerosis. J Am Coll Cardiol 2018;72:1804–13. 5 Allanore Y, Meune C, Vonk MC, et al. Prevalence and factors associated with left ventricular dysfunction in the EULAR scleroderma trial and research group (EUSTAR) database of patients with systemic sclerosis. Ann Rheum Dis 2010;69:218–21. 6 Liem SIE, Boonstra M, le Cessie S, et al. Sex- specific risk of anti- topoisomerase antibodies on mortality and disease severity in systemic sclerosis: 10- year analysis of the Leiden CCISS and EUSTAR cohorts. The Lancet Rheumatology 2022;4:e699–709. 7 Yiu KH, Schouffoer AA, Marsan NA, et al. Left ventricular dysfunction assessed by speckle- tracking strain analysis in patients with systemic sclerosis: relationship to functional capacity and ventricular arrhythmias. Arthritis Rheum 2011;63:3969–78. 8 Gegenava T, Gegenava M, Steup- Beekman GM, et al. Left ventricular systolic function in patients with systemic lupus erythematosus and its association with cardiovascular events. J Am Soc Echocardiogr 2020;33:1116–22. 9 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. 10 van Wijngaarden SE, Boonstra M, Bloem B, et al. Clinical and echocardiographic associates of all- cause mortality and cardiovascular outcomes in patients with systemic sclerosis. JACC Cardiovasc Imaging 2019;12:2273–6. 11 Schouffoer AA, Ninaber MK, Beaart- van de Voorde LJJ, et al. Randomized comparison of a multidisciplinary team care program with usual care in patients with systemic sclerosis. Arthritis Care Res (Hoboken) 2011;63:909–17. 12 Meijs J, Schouffoer AA, Ajmone Marsan N, et al. Therapeutic and diagnostic outcomes of a standardised, comprehensive care pathway for patients with systemic sclerosis. RMD Open 2016;2:e000159. 13 van den Hoogen F, Khanna D, Fransen J, et al. Classification criteria for systemic sclerosis: an American r/European League against rheumatism collaborative initiative. Ann Rheum Dis 2013;72:1747–55. 14 LeRoy EC, Medsger TA. Criteria for the classification of early systemic sclerosis. J Rheumatol 2001;28:1573–6. 15 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. 16 Pellegrino R, Viegi G, Brusasco V, et al. Interpretative strategies for lung function tests. Eur Respir J 2005;26:948–68. 17 Graham BL, Brusasco V, Burgos F, et al. Executive summary: 2017 ERS/ATS standards for single- breath carbon monoxide uptake in the lung. Eur Respir J 2017;49:16E0016. Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from 10 Gegenava T, et al. RMD Open 2023;9:e003380. doi:10.1136/rmdopen-2023-003380 RMD Open 18 Hansell DM, Bankier AA, MacMahon H, et al. Fleischner society: glossary of terms for thoracic imaging. Radiology 2008;246:697–722. 19 Yiu KH, Ninaber MK, Kroft LJ, et al. Impact of pulmonary fibrosis and elevated pulmonary pressures on right ventricular function in patients with systemic sclerosis. Rheumatology (Oxford) 2016;55:504–12. 20 Ross RM. ATS/ACCP statement on cardiopulmonary exercise testing. Am J Respir Crit Care Med 2003;167:1451; author reply 1451. 21 Vacca A, Meune C, Gordon J, et al. Cardiac arrhythmias and conduction defects in systemic sclerosis. Rheumatology (Oxford) 2014;53:1172–7. 22 Nagueh SF, Smiseth OA, Appleton CP, et al. Recommendations for the evaluation of left ventricular diastolic function by echocardiography: an update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging. Eur Heart J Cardiovasc Imaging 2016;17:1321–60. 23 Kircher BJ, Himelman RB, Schiller NB. Noninvasive estimation of right atrial pressure from the Inspiratory collapse of the inferior vena cava. Am J Cardiol 1990;66:493–6. 24 Rudski LG, Lai WW, Afilalo J, et al. Guidelines for the echocardiographic assessment of the right heart in adults: a report from the American Society of Echocardiography endorsed by the European Association of Echocardiography, a registered branch of the European Society of Cardiology, and the Canadian Society of Echocardiography. J Am Soc Echocardiogr 2010;23:685–713; 25 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. Eur Heart J Cardiovasc Imaging 2015;16:233–70. 26 Delgado V, Ypenburg C, van Bommel RJ, et al. Assessment of left ventricular dyssynchrony by speckle tracking strain imaging comparison between longitudinal, circumferential, and radial strain in cardiac resynchronization therapy. J Am Coll Cardiol 2008;51:1944–52. 27 Perelas A, Silver RM, Arrossi AV, et al. Systemic sclerosis- associated interstitial lung disease. Lancet Respir Med 2020;8:304–20. 28 Guillén- Del- Castillo A, Simeón- Aznar CP, Callejas- Moraga EL, et al. Quantitative videocapillaroscopy correlates with functional respiratory parameters: a clue for vasculopathy as a pathogenic mechanism for lung injury in systemic sclerosis. Arthritis Res Ther 2018;20:281. 29 Mayes MD, Lacey JV Jr, Beebe- Dimmer J, et al. Prevalence, incidence, survival, and disease characteristics of systemic sclerosis in a large US population. Arthritis Rheum 2003;48:2246–55. 30 Hunzelmann N, Genth E, Krieg T, et al. The Registry of the German network for systemic scleroderma: frequency of disease subsets and patterns of organ involvement. Rheumatology (Oxford) 2008;47:1185–92. 31 Ostojić P, Damjanov N. Different clinical features in patients with limited and diffuse cutaneous systemic sclerosis. Clin Rheumatol 2006;25:453–7. 32 Lau ES, Paniagua SM, Guseh JS, et al. Sex differences in circulating biomarkers of cardiovascular disease. J Am Coll Cardiol 2019;74:1543–53. 33 Tennøe AH, Murbræch K, Didriksen H, et al. Serum markers of cardiac complications in a systemic sclerosis cohort. Sci Rep 2022;12:4661. 34 Tzelepis GE, Kelekis NL, Plastiras SC, et al. Pattern and distribution of myocardial fibrosis in systemic sclerosis: a delayed enhanced magnetic resonance imaging study. Arthritis Rheum 2007;56:3827–36. Tiedekuntakirjasto. Protected by copyright. o n D ecem ber 11, 2023 at Turun Yliop Laaketieteellinen http://rm dopen.bmj.com/ R M D O pen: first published as 10.1136/rm dopen-2023-003380 on 10 Novem ber 2023. Downloaded from