Introduction

Systemic sclerosis (SSc) is a rare, autoimmune disease characterized by diffuse fibrosis, inflammation and vasculopathy and remains one of the most complex and lethal diseases among connective tissue disorders1. The severity of organ involvement predicts the outcome in SSc patients2,3.

Pulmonary complications, mainly pulmonary arterial hypertension (PAH) and interstitial lung disease (ILD), remain the leading cause of mortality among SSc patients1. Consequently, screening and early diagnosis of both the disease and organ involvement are crucial for prognosis.

The establishment of risk stratification methods is essential to improve the survival through rapid medical and therapeutic interventions. The majority of available scoring models are yet to be validated1. Previous studies have already reported the prognostic meaning of impaired oxygen saturation (SaO2) during exercise for survival and disease progression among patients with SSc-ILD4,5. SaO2 may be influenced by both cardiac and pulmonary comorbidities. Furthermore, anemia has been shown to be a prognostic predictor of survival among patients with SSc6. A combination of these parameters might therefore be beneficial for risk stratification and prognosis.

Trudzinski et al. showed that apart from white blood cell count, oxygenated hemoglobin (OxyHem), the product of hemoglobin concentration and fractional oxygen saturation (SaO2) in arterial capillary blood gas analysis, is a suitable and easy to estimate novel biomarker for the prognosis of patients with chronic obstructive pulmonary disease (COPD)7. OxyHem measured in patients with stable COPD, but also during acute COPD exacerbations was independently associated with higher mortality among those patients8. OxyHem in arterial blood may also be a novel biomarker for the prognosis in patients with SSc.

Hence, the aim of this study was to analyze the prognostic value of OxyHem in SSc patients referred for screening for pulmonary hypertension (PH). Furthermore, the study investigated the association of OxyHem with the development of pulmonary vascular disease (PVD) among these patients.

Results

Baseline clinical and demographic characteristics

A total of 280 patients with SSc were screened for pulmonary hypertension. From this cohort, 13 patients met exclusion criteria, leading to a final study group of 267 patients (Fig. 1, Table 1). The mean age of the patients was 59.8 ± 13.4 years, 82% were female, 73.8% had limited cutaneous SSc (lcSSc) and 26.2% diffuse cutaneous SSc (dcSSc). The mean duration of SSc at baseline was 9.0 ± 9.3 years. Of these, 70 had significant ILD (26.2%) and 141 (52.8%) showed a mean pulmonary arterial pressure (mPAP) ≥ 21 mmHg (Fig. 1) while 126 had no signs of PH (41.6%). Of the 141 patients with PH 19 presented with PH due to left heart disease (Fig. 1). Patients with PAH received targeted PAH medication according to the contemporary guidelines9. Out of 267 patients, 46.4% were in World Health Organization functional class (WHO-FC) II, 30.3% in WHO-FC III and 2.6% in WHO-FC IV. The mean 6-min walking distance (6MWD) was 427.7 ± 102.5 m. Comorbidities were identified in 166 patients: 101 (37.8%) suffered from arterial hypertension, 58 (21.7%) had coronary heart disease, one (0.4%) had reported gastric antral vascular ectasia, three (1%) had lymphoma in the past history and 41 (15.4%) were obese, having a body mass index (BMI) > 30 kg/m2.

Figure 1
figure 1

Study flow-chart. The flow-chart of the study cohort gives reasons for exclusion from the analysis set, distribution of patients according to OxyHem and to hemodynamic criteria.

Table 1 Clinical characteristics of the study cohort.

Clinical parameters characterizing patients with low OxyHem

The median OxyHem in this cohort was 12.89 g/dl (interquartile range: 151.8). Receiver operating characteristic (ROC) analyses for OxyHem were in line with the cutoff from the literature7 of 12.5 g/dl (ROC area under the curve P = 0.004). Overall, 94 SSc patients had an OxyHem ≤ 12.5 g/dl and 173 SSc patients OxyHem values above 12.5 g/dl (Table 2). Compared to patients with OxyHem > 12.5 g/dl, patients with low OxyHem were mainly female (n = 87, 92.6% vs. n = 132, 76.3%, P = 0.001), had mostly lcSSc (n = 75, 79.8% vs. n = 122, 70.5%, P = 0.011), a more impaired iron metabolism in terms of lower ferritin levels (75.3 ± 83.6 vs. 126.7 ± 150.6 µg/l, P = 0.004) in normal range, mean corpuscular volume of erythrocytes (MCV) (88.1 ± 8.7 vs. 90.3 ± 5.2 fl, P = 0.027) and mean corpuscular hemoglobin in erythrocytes (29.1 ± 3.0 vs. 30.1 ± 2.4 pg, P = 0.006). Moreover, they were clinically more limited in terms of WHO-FC (P = 0.001), shorter 6MWD (408.6 ± 101.3 vs. 437.8 ± 102.0 m, P = 0.036) and greater hemodynamic impairment with a higher pulmonary vascular resistance (PVR) (3.2 ± 2.6 vs. 2.4 ± 1.9 Wood Units (WU), P = 0.018) and mPAP (25.6 ± 11.7 vs. 21.6 ± 9.3 mmHg, P = 0.005) (Table 2). There was no significant difference in N-terminal pro brain natriuretic peptide (NTproBNP) between the two groups. Out of 94 patients with OxyHem ≤ 12.5 g/dl, 29 suffered from ILD (23.7%). In 59 patients with low OxyHem, mPAP was ≥ 21 mmHg (Fig. 1). Concomitant ILD was diagnosed in 70 patients, out of which 50 showed an mPAP ≥ 21 mmHg. Elevation of pulmonary arterial wedge pressure (PAWP) > 15 mmHg was detected in 22 patients, 19 showed an mPAP ≥ 21 mmHg (Fig. 1). Six patients had PH due to both left heart and lung disease. Pulmonary function test (PFT) parameters were significantly more impaired among patients with lower OxyHem including maximal vital capacity (VCmax), forced expiratory volume in one second (FEV1) and total lung capacity (TLC). In the laboratory analysis C-reactive protein (CRP) was higher in the low OxyHem group (P = 0.020) and these patients had significantly lower hemoglobin (Hb) (P < 0.0001). There was no difference in renal function between the two groups (Table 2).

Table 2 Clinical characteristics of the OxyHem subgroups.

Survival in SSc patients

Within the observation period of 3.2 ± 2.6 (median 2.7; interquartile range 3.5) years 47 patients died (17.6%): 22 due to PH (46.8%), 4 due to PH and/or ILD (8.5%), 3 due to ILD (6.4%), 7 due to cancer (14.9%), one due to pulmonary sepsis (2.1%) and 3 due to left heart disease (6.4%). Two patients underwent lung transplantation due to PH and ILD. In seven cases the cause of death remained unknown (14.9%). The mean estimated overall survival was 7.7 ± 0.3 standard error of the mean years (95% confidence interval 7.1 to 8.2) from baseline. Survival was similar in male and female patients (P = 0.624) and between patients with dcSSc and lcSSc (P = 0.694).

In the univariable Kaplan–Meier analysis, OxyHem ≤ 12.5 g/dl (P = 0.046), age at baseline of ≥ 60 years (P = 0.002), the presence of ILD (P < 0.0001), high PVR ≥ 2 WU (P = 0.002) and low DLCO ≤ 65% predicted (P < 0.001) were significantly associated with worse survival (Table 3). White blood cell (WBC) count ≥ 10,000/ml was analyzed together with OxyHem in the initially described COPD cohort7. It was also a significant prognostic parameter for survival in our cohort (P = 0.02).

Table 3 Uni- and multivariable Cox-regression analysis for survival.

Patients with low OxyHem ≤ 12.5 g/dl showed significantly worse survival than patients with OxyHem > 12.5 g/dl with 1-, 3-, 5-year survival of 98.7%, 82.9% and 71.2% versus 99.3%, 95.7% and 80.3%, respectively (Kaplan–Meier P = 0.046, Fig. 2). In the multivariable Cox-regression analysis, low DLCO, the presence of ILD and older age were identified as independent prognostic predictors for survival (combined P < 0.0001, Fig. 3).

Figure 2
figure 2

Kaplan–Meier analysis of OxyHem. Patients with an OxyHem value equal and lower than 12.5 g/dl had a significantly shorter survival in our study cohort.

Figure 3
figure 3

Kaplan–Meier analysis of multivariable risk set. A multivariable risk set including DLCO ≤ 65% predicted, ILD and age ≥ 60 years was used for risk stratification.

OxyHem in pulmonary vascular disease

At baseline, patients with PH (mPAP ≥ 21 mmHg) had the highest rate of low OxyHem (41.8%), compared to patients with normal hemodynamics (27.8% low OxyHem; chi-square P = 0.016). Low OxyHem at follow-up was though not associated with hemodynamic status (P = 0.562) (Table 4). PH at baseline could be predicted with 41.8% sensitivity and 72.7% specificity by low OxyHem. When combining low DLCO and low OxyHem, prediction of PH at baseline had a sensitivity of 76.1% and a specificity of 50.6%. In comparison, step 1 of the DETECT algorithm10 reached 85.3% sensitivity and 37.0% specificity. In step 2 of the DETECT algorithm10, 79.4% sensitivity and 73.6% specificity could be reached for PH at baseline.

Table 4 Frequency of low OxyHem in patients with PH or PAH.

When using the three independent risk factors DLCO ≤ 65% predicted, age ≥ 60 years and presence of ILD, occurrence of PH with mPAP ≥ 21 mmHg during follow-up could be predicted with a sensitivity of 91.1% and a specificity of 34.6% when at least one of the three risk factors was met. The combination of DLCO ≤ 65% predicted and OxyHem ≤ 12.5 g/dl had a sensitivity of 58.7% and a specificity of 72.9% to predict PH during follow-up when at least one of the two risk factors was met (prediction of 44 out of 57; true negative 35 out of 66).

Discussion

This study detected for the first time that SSc patients with an OxyHem level ≤ 12.5 g/dl have a significantly worse survival than patients with an OxyHem level > 12.5 g/dl (P = 0.046). ILD, low DLCO, older age were robust prognostic factors associated with survival in the multivariable analysis. Furthermore, low OxyHem in combination with low DLCO at baseline were predictors of PH at baseline and during follow-up (sensitivity 76.1% and 58.7%, specificity of 50.6% and 72.9%, respectively). To our knowledge, this is the first study indicating the impact of OxyHem as predictor of survival among SSc patients.

OxyHem and SSc prognosis

Lower OxyHem was identified in patients with a shorter survival, highlighting the impact of oxygenation and the significance of oxygen transport among patients suffering from extended fibrosis and vasculopathy, that typically characterize SSc11. Patients with low OxyHem revealed more often clinical characteristics with worse survival in previously described risk stratification models, such as greater WHO-FC12, elevated CRP1, presence of ILD1, lower VCmax and DLCO13 as well as lower Hb and shorter 6MWD12,13. Patients with lcSSc and female patients had significantly more often a low OxyHem level in our cohort. As these patients would not be suspected to be at risk for an impaired outcome, this underlines the importance of evaluating patients also at low risk1,12,13,14,15,16.

SSc and oxygenation

The importance of oxygen saturation after exercise as a predictor of mortality among SSc-ILD patients was previously evaluated5. Exercise induced desaturation was also shown to be predictor of mortality among patients at risk for PH in the PHAROS registry, though SaO2 at rest did not differ between the groups17 and was found to be associated with 3-year-survival in the DIBOSA study18. Wu et al. reported that desaturation < 94% in 6MWD combined with the presence of arthritis was associated with disease progression in SSc-ILD4. In contrast, for the calculation of OxyHem, resting oxygen saturation can be used, which makes the parameter easier to determine and less exhausting for the patient compared to assessment of oxygen saturation during exertion.

PVD

In our cohort, OxyHem was associated with pulmonary vasculopathy only at baseline but not at follow-up. This might be caused by initiation of PAH targeted therapy after diagnosis. Hoeper et al. reported that oxygenation was not a prognostic marker for survival among prevalent patients with idiopathic PAH under therapy, which supports our results19. Recently, Valentin et al. (2021) showed that oxygen desaturation (SaO2) ≥ 3% at the time of the first follow-up after one year in combination with low DLCO was significantly associated with worse prognosis of PAH20, matching our results linking low OxyHem and low DLCO with early PVD. In this cohort, male sex did not correlate with survival consistent with results from the literature6,21,22.

Strengths and limitations of the study

This study provides important insights into prognostic predictors and stratification models in SSc patients screened for PH. The screening for PH was performed in an expert center including a thorough diagnostic workup, leading to a well characterized study cohort. OxyHem offers a parameter, which is pathophysiologically important in different organ involvements in SSc such as ILD, PH and cardiac diseases. While it was a significant predictor in the univariable analysis it was outperformed in the multivariable model. This could be due to the limited sample size or the strong effect of ILD on prognosis. Since we only included extensive ILD in our assessment, this could have also biased the results. One of the main limitations of this study is its retrospective character leading to a smaller data set due to missing parameters. Early prognostic markers need to be investigated in future, prospective studies. We also had no control cohort of healthy individuals or other well characterized patients to compare our results to. However, the prognostic value of OxyHem has already been described in COPD patients7. Furthermore, other organ manifestations of SSc, such as gastrointestinal manifestations and renal crisis, could not be considered due to missing data. Although the cohort size of our study is large for such a rare disease, sample size of subgroups was too small to differentiate the impact of OxyHem on survival separately in ILD and PH. Another limitation is that in our cohort blood gas analysis was performed with arterialized capillary blood from the earlobe, which might not be included in standard screening procedures by rheumatologists. However, this method has already shown good correlation with arterial puncture23. The meaning of OxyHem during the course of the disease including follow-up data should be addressed in future studies. Finally, our cohort included only SSc patients referred for screening for PH, which could have led to biased data.

Conclusions

The results of this study revealed that OxyHem may be a new biomarker which could be used for risk stratification in SSc and identification of patients with early PVD and PH, respectively. OxyHem is easy to calculate in daily routine and values ≤ 12.5 g/dl were significantly associated with a worse survival. Moreover, the combination of DLCO ≤ 65% predicted and low OxyHem could predict the development of early PVD in this cohort. Further prospective studies are needed to confirm the results.

Methods

Study population

SSc patients with dcSSc or lcSSc 24 referred to our center for PH at the Thoraxklinik Heidelberg gGmbH at Heidelberg University Hospital, Germany for screening from their rheumatologists between 2008 and 2020 were retrospectively analyzed in this single center, cross-sectional cohort study. All patients fulfilled the SSc classification criteria of the American College of Rheumatology/European League against Rheumatism25 and were followed yearly according to the current recommendations9. Individuals were excluded, if they were underaged, had rheumatic diseases other than SSc, had missing data for the calculation of OxyHem or were unable to provide informed consent. Data of part of this cohort has already been analyzed and published before26. The ethics committee of the Medical Faculty of Heidelberg University (internal number S-126/2021) had no objection against the conduct of this study. All patients provided written informed consent for the analysis of the data. The study complied with the Declaration of Helsinki in its current version.

Study design

Data from routinely performed examinations were extracted from the patient's files at the first presentation in the PH center. Patients were referred from their rheumatologists for screening for PH due to unexplained breathlessness and/or suspicion for PH. A detailed clinical work-up was performed in all patients. The routinely performed screening included medical history, recording of comorbidities such as arterial hypertension, coronary heart disease or obesity (BMI > 30 kg/m2), physical examination, electrocardiogram, PFT with DLCO, blood gas analysis obtained from arterialized capillary blood from earlobe, echocardiography at rest and during exercise, high resolution computed tomography scan of the lungs (HRCT) to investigate the presence of any pulmonary manifestations or diseases, WHO-FC, 6MWD test under standardized conditions27, laboratory examinations including blood count tests, NTproBNP, CRP, analysis for iron deficiency, renal function including glomerular filtration rate, and right heart catheterization at rest and/or during exercise.

The presence of cardiac disease, such as coronary heart disease, and pulmonary manifestations, such as PH or ILD, was documented and, in case of pathology, furtherly examined. ILD was only considered in case of extensive ILD defined as the presence of significant interstitial fibrosis in the HRCT of the lungs (> 20% of the parenchyma) or restrictive pattern in PFT (VCmax < 70% predicted with a normal ratio of FEV1 to VC) in case of missing quantification of interstitial fibrosis in the HRCT28. Manifest PAH was diagnosed according to valid guidelines at the time of the study (mPAP ≥ 25 mmHg, PAWP ≤ 15 mmHg and PVR > 3 WU) measured by RHC9. Signs of early PVD were defined as mPAP ≥ 21 mmHg. PH due to left heart disease was defined by an mPAP ≥ 21 mmHg with a PAWP > 15 mmHg29. PH due to lung disease was based on lung function parameters and in case of ILD HRCT imaging30.

Statistical methods

Statistical analyses were performed by a medical statistician (NB). Data are presented as mean ± standard deviation (SD) with 95% confidence interval of the mean. Frequency data are presented as n and %. The concentration of OxyHem, (g/dl) was calculated at baseline by the multiplication of Hb concentration in blood with the SaO2 measured during blood gas analysis as described before7.

Survival analysis was performed using the Kaplan-Meier and Cox-regression analysis. The first screening visit date was set as baseline for the survival analysis. Death was defined as death due to any cause. Thresholds for dichotomization of parameters were either obtained by the literature of known prognostic predictors or by ROC in case of unknown threshold or in need of threshold adjustment. Clinical parameters were compared with the two-sided student’s t-test for independent samples for patients above versus below the OxyHem threshold. Qualitative data was analyzed using the chi-square test, respectively. This study is of exploratory nature. P-values < 0.05 were considered as statistically significant and were not adapted for multiple testing.

Variables were evaluated by univariable analysis to identify their impact on survival. Parameters were then analyzed with multivariable Cox-regression analysis. OxyHem was compared to known prognostic factors in SSc, which were identified by literature search. These included sex, type of SSc, age ≥ 60 years at baseline31, presence of ILD1,15, PVR ≥ 2 WU26, DLCO ≤ 65% predicted13. In addition, WBC count was included as it showed significant predictive power aside from OxyHem in the OxyHem discovery cohort7. For the multivariable Cox-regression all parameters listed above were considered.

Survival analysis was performed in subsets of patients with signs of cardiovascular disease with either increased PVR ≥ 2 WU, or with mPAP ≥ 25 mmHg. All analyses have been performed using IBM SPSS version 27 (SPSS Statistics V27, IBM Corporation, Somers, New York).

Ethics approval and consent to participate

All data were pseudonymized, and the study was approved by the ethics committee of the medical faculty of Heidelberg University Hospital (internal number S-126/2021). The study complies with the Declaration of Helsinki in its current version. The patients participating in this study provided written informed consent.