Vira Tseluyko, MD, PhD, Professor, Honored leader of science and technics of Ukraine, Head of the Department of cardiology, laboratory and functional diagnostics, medical faculty, V.N. Karazin Kharkiv National University, Nezalezhnosti sq., 4, Kharkiv, 61022, Ukraine.
Tseluyko V.I, Yakovleva L.M., Mishchuk N.E., Kurinna M.V., Kharchenko L.V., Askierov R.N., Shylo N.G, (2024). Risk Stratification of In-Hospital Mortality in Patients with High-Risk Pulmonary Embolism. Clinical Cardiovascular Research. 3(2); DOI: 10.58489/2836-5917/025
© 2024 Vira Tseluyko, this is an open-access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
high-risk pulmonary embolism, additional risk, in-hospital mortality risk stratification formula.
Objective. To identify factors associated with in-hospital mortality and to develop and validate a formula for additional risk of early mortality in patients (pts) with high-risk pulmonary embolism (PE).
Materials and Methods. The study consisted of two phases. In Phase I, a retrospective analysis was conducted on the medical records of 113 pts with high-risk PE who were consecutively hospitalized at Kharkiv City Clinical Hospital No. 8 from 01.01. 2017 to 01.01.2023. The mean age of pts was 65.14±14.13 years, with 45 (38.8%) males and 68 (60.2%) females. Based on the course of the disease, the pts were divided into two groups: Group 1 – 63 pts (55.8%) who improved and were discharged, and Group 2 – 50 pts (44.2%) who died during hospitalization. Phase II of the study involved a multicenter cohort (9 centers in Ukraine that provided data on pts with high-risk pulmonary embolism hospitalized from January 2023 to December 2023). This phase included 75 pts with high-risk PE, with mean age of 60.60±13.40 years; 38 (50.7%) were male, and 37 (49.3%) were female. Group I consisted of 65 patients (80%) discharged with improvement, Group II comprised 15 patients (20%) who died during hospitalization. Clinical and anamnestic and laboratory-instrumental data were analyzed, statistical analysis of the resulting figures was performed.
Results and Discussion. In Phase I of the study, a formula for additional risk stratification of early mortality in pts with high-risk PE was developed. Using ROC analysis, the sensitivity of the resulting prognostic formula was found to be 64.0%, and the specificity was 92.1%. Based on this formula, we proposed a simplified scale called SBAFS (S - saturation, B - bifurcation, A - arterial hypotension, F - fraction of ejection, S - female sex). In Phase II of the study, the sensitivity and specificity of the obtained formula were validated using ROC analysis: the sensitivity was 93.3%, and the specificity was 74.6%, confirming the results of Phase I. When evaluating the proposed simplified SBAFS scale, it was also found that in patients of Group II, the mean score on the scale was (2.53±0.83) and was higher than that in Group I: (1.57±0.96), p=0.0006.
Conclusions. The main factors associated with the risk of in-hospital mortality in pts with high-risk PE are decreased oxygen saturation, thrombus localization in the bifurcation of the pulmonary artery as determined by СТPA angiography, arterial hypotension, female sex, and decreased left ventricular ejection fraction as assessed by echocardiography. The previously proposed SBAFS formula confirmed its significance upon validation involving data on PE pts from 9 regions of Ukraine, with a sensitivity of 93% and specificity of 75%.
Pulmonary embolism (PE) ranks third, following myocardial infarction and stroke, in the structure of cardiovascular mortality [1,2]. A significant number of risk factors and the presence of comorbidities contribute to substantial variability in both the course and outcomes of PE [3,6]. In light of this, current guidelines recommend conducting risk stratification for PE, considering the predicted early mortality [1,6]. It has been established that in pts with high-risk PE, the clinical course of the disease is the most severe, and the prognosis is extremely unfavorable; hospital mortality rates in various studies can range from 25.4% to 70% [3,4]. Currently, patients with high-risk PE are classified as those with hemodynamic instability, which is defined by criteria such as cardiac arrest requiring cardiopulmonary resuscitation, the presence of obstructive shock, or hypotension lasting more than 15 minutes with a systolic blood pressure (SBP) of less than 90 mm Hg, or a decrease in SBP of 40 mm Hg from baseline without obvious reasons [1]. The question of additional clinical and instrumental factors related to in-hospital mortality in pts with high-risk PE remains insufficiently studied.
The aim of the study was to identify factors associated with the risk of in-hospital mortality and to develop and validate a formula for additional risk of early death in pts with high-risk pulmonary embolism.
In Phase I of our study, a retrospective analysis of the medical records of 113 patients with high-risk pulmonary embolism (PE) was conducted. These pts were consecutively hospitalized at Kharkiv City Clinical Hospital No. 8 from January 1, 2017, to January 1, 2023 (Table 1). The mean age of the pts was 65.14±14.13 years, with 45 men (38.8%) and 68 women (60.2%).
Phase II of the study was conducted on a multicenter cohort involving centers from Kyiv, Sumy, Kremenchuk, Poltava, Mykolaiv, Dnipro, Lutsk and Kharkiv. These centers provided information on pts with high-risk PE who were treated from January 2023 to December 2023 (Table 4). In total, 75 pts with high-risk PE were included in Phase II, with a mean age of 60.60±13.40 years; among them, there were 38 men (50.7%) and 37 women (49.3%).
The diagnosis of PE was verified through computed tomography pulmonary angiography (CTPA) and/or post-mortem findings, while risk assessment was performed in accordance with the European guidelines for the treatment of pts with PE (2019) [1].
Statistical Analysis. The statistical analysis of the obtained data was performed using the Statistica 10.0 software package (StatSoft Inc, USA) and Microsoft Office Excel 2013. Since the preliminary analysis of quantitative characteristics did not reveal deviations from a normal distribution, they are presented as mean ± standard deviation (M ± σ). The Mann-Whitney U-test was used to compare the means of two samples. Differences were considered statistically significant at p < 0.05. Univariate and multivariate logistic regression analyses were employed to identify predictors of the combined clinical endpoint. The β coefficient, standard error, odds ratio (OR), and 95% confidence interval (CI) were calculated for each factor. Discrimination was assessed using the area under the receiver operating characteristic curve (AUC).
In the first phase pts with high-risk PE were divided into two groups: Group 1 – 63 (55.8%), who were discharged with improvement; Group 2 – 50 (44.2%), who died. (Table 1). Nearly 70% of pts from Group II died within 1,5 hours of arriving at the hospital.
Table 1: Comparative Characteristics of High-Risk Patients Based on Prognosis
Parameters |
All (n=113) |
Group 1 (n=63) |
Group 2 (n=50) |
Р |
Age (М ± σ), years |
65,14±14,13 |
64,92±14,08 |
65,44±14,33 |
0,768 |
Men, n (%) |
45 (38,8) |
32 (50,8) |
13 (26,0) |
0,007 |
Women, n (%) |
68 (60,2) |
31 (49,2) |
37 (74,0) |
|
Clinical signs at admission |
||||
HR (М ± σ), bpm |
105,92±20,77 |
103,58±20,28 |
109,04±21,23 |
0,142 |
SBP (М ± σ), mm Hg |
78,14±18,52 |
85,79±12,42 |
67,87±20,42 |
0,00001 |
Body temperature (М ± σ), °С |
36,39±0,52 |
36,42±0,59 |
36,35±0,43 |
0,276 |
SpO2 (М ± σ), % |
82,28±10,43 |
86,41±6,65 |
77,13±11,96 |
0,00004 |
Clinical characteristics |
||||
History of previous VTE, n (%) |
18 (16,2) (15,9) |
13 (20,6) |
5(10,0) |
0,202 |
Varicose veins, n (%) |
28 (25,2) |
13 (20,6) |
15 (31,3) (30) |
0,202 |
Obesity, n (%) |
40 (36,0) |
14(22,2) |
26 (54,2) (52) |
0,0005 |
DVT, n (%) |
42 (37,8) |
17(27,0) |
25 (52,1) (50) |
0,007 |
SVT, n (%) |
28 (25,2) |
13 (20,6) |
15 (31,3) (30) |
0,202 |
Congestive heart failure or respiratory failure, n (%) |
28 (25,2) |
16 (25,4) |
12 (25,0) (24) |
0,962 |
AF, n (%) |
27 (24,3) |
15 (23,8) |
12 (25,0) (24) |
0,885 |
DM, n (%) |
22 (19,8) |
12 (19,0) |
10 (20,8) (20) |
0,815 |
AH, n (%) |
88 (79,3) |
47 (74,6) |
41 (85,4) (82) |
0,248 |
History of cancer, n (%) |
11 (9,9) |
6 (9,5) |
5 (10,4) (10) |
0,56 |
Infectious disease, n (%) |
15 (13,5) |
8 (12,7) |
7 (14,6) (14) |
0,994 |
Level of lesion verified by CTPA |
||||
Bifurcation of the pulmonary trunk (saddle PE) n (%) |
26 (23,4) |
12 (19,0) |
14 (29,2) |
0,212 |
Lobar branches of PA, n (%) |
54 (48,6) |
32 (50,8) |
22 (45,8) |
0,605 |
Bilateral segmental branches of PA, n (%) |
10 (9,0) |
9 (14,3) |
1 (2,1) |
0,025 |
Unilateral segmental branches of PA, n (%) |
2 (1,8) |
2 (3,2) |
0 |
0,32 |
Parameters obtained by Echocardiography |
||||
LPA (М ± σ), mm |
25,49±2,68 |
25,44±2,75 |
25,73±2,41 |
0,426 |
RPA (М ± σ), mm |
26,00±3,91 |
25,90±3,90 |
26,45±4,13 |
0,8 |
PT (М ± σ), mm |
31,76±3,88 |
31,62±4,03 |
32,45±3,21 |
0,515 |
LA (М ± σ), mm |
39,63±5,81 |
39,60±6,11 |
39,80±3,61 |
0,948 |
RA (М ± σ), mm |
45,90±6,87 |
45,58±6,20 |
47,90±10,33 |
0,225 |
RV (М ± σ), mm |
36,23±10,28 |
35,56±9,93 |
41,38±12,19 |
0,216 |
LV EDV (М ± σ), mm |
45,56±9,62 |
45,69±10,14 |
44,70±5,76 |
0,864 |
LV ESV (М ± σ), мм |
32,19±6,64 |
32,00±6,82 |
33,40±5,54 |
0,449 |
Interventricular septum (М ± σ), mm |
10,58±1,75 |
10,42±1,60 |
11,45±2,30 |
0,157 |
Cardiac output (М ± σ), ml |
52,92±18,71 |
54,12±19,11 |
46,50±15,94 |
0,476 |
Diameter of aorta (М ± σ), mm |
33,96±4,55 |
34,30±4,31 |
31,67±5,74 |
0,38 |
MPAP (М ± σ), mm Hg |
51,70±16,90 |
50,78±16,86 |
58,00±16,85 |
0,163 |
LV EF (М ± σ), % |
55,52±9,45 |
56,43±9,25 |
49,90±9,16 |
0,008 |
Note. AF – atrial fibrillation, AH – arterial hypertension, CTPA – computed tomography of pulmonary arteries, DM – diabetes mellitus, DVT – deep vein thrombosis, HR – heart rate, LA – left atrium, LPA – left pulmonary artery, LV EDV – end diastolic volume of left ventricle, LV ESV – end systolic volume of left ventricle, LV EF – ejection fraction of left ventricle, MPAP – mean pressure in pulmonary artery, RA – right atrium, RPA – right pulmonary artery, RV – right ventricle , SBP – systolic blood pressure, SpO2 – blood oxygen saturation, SVT – superficial vein thrombosis.
There were more women in Group 2 (p=0,006). Also pts of Group 2 were more severe having significantly lower SpO2 values (p=0.000043) and SBP (p=0.00001). In Group 2, risk factors of PE, such as obesity (p=0.0005) and DVT (p=0.007), were more common than in Group 1. According to CTPA, pts of Group 2 had more bilateral segmental lesions (p=0.025) and (according to echocardiography) lower LVEF (p=0.008).
We used uni- and multivariate analyses to find out independent prognostic factors associated with in-hospital mortality in pts with high-risk PE (Table 4). The regression analysis revealed a significant independent association between in-hospital mortality and the presence of thrombus in the bifurcation of pulmonary trunk (p=0.015) and female sex (p=0.0048). Additionally, an inverse independent association was found between mortality and SBP (p=0.0008), SpO₂ (p=0.0014) and LVEF (p=0.011). There were more women in Group 2 (p=0,006). Also pts of Group 2 were more severe having significantly lower SpO2 values (p=0.000043) and SBP (p=0.00001). In Group 2, risk factors of PE, such as obesity (p=0.0005) and DVT (p=0.007), were more common than in Group 1. According to CTPA, pts of Group 2 had more bilateral segmental lesions (p=0.025) and (according to echocardiography) lower LVEF (p=0.008).
We used uni- and multivariate analyses to find out independent prognostic factors associated with in-hospital mortality in pts with high-risk PE (Table 4). The regression analysis revealed a significant independent association between in-hospital mortality and the presence of thrombus in the bifurcation of pulmonary trunk (p=0.015) and female sex (p=0.0048). Additionally, an inverse independent association was found between mortality and SBP (p=0.0008), SpO₂ (p=0.0014) and LVEF (p=0.011).
Table 2: Factors Affecting In-Hospital Mortality in High-Risk Patients with Pulmonary Embolism (Stage I)
Parameters |
Dependent variable: mortality. |
|||||||
univariate linear regression analysis (χ2=77,35; P <0,0001) AUC=0.922(0.856-0.964) |
multivariate linear regression analysis (χ2=63,11; P<0,0001) AUC=0.888(0.815-0.940) |
|||||||
β- coefficient |
OR |
95% СI |
Р |
β- coefficient |
OR |
95% CI |
Р |
|
SpO2, n (%) |
-0,118 |
0,889 |
0,82- 0,97 |
0,007 |
-0,114 |
0,893 |
0,83-0,96 |
0,0014 |
АH, n (%) |
3,483 |
32,558 |
2,37-446,41 |
0,009 |
|
|
|
|
Thrombus in bifurcation of PT (%) |
2,292 |
9,894 |
1,36-71,84 |
0,024 |
1,58 |
4,857 |
1,35-17,48 |
0,015 |
SBP, n (%) mm Hg |
-0,127 |
0,881 |
0,82- 0,95 |
0,001 |
-0,091 |
0,913 |
0,87-16,77 |
0,0008 |
Female sex, n (%) |
2,363 |
10,624 |
2,03-55,53 |
0,005 |
1,663 |
5,273 |
1,66-16,76 |
0,0048 |
LV EF, n (%) |
-0,097 |
0,908 |
0,83-0,99 |
0,044 |
-0,078 |
0,924 |
0,87-0,98 |
0,011 |
Note. AH – arterial hypertension, LV EF – ejection fraction of left ventricle, SBP – systolic blood pressure, SpO2 – blood oxygen saturation.
To determine the threshold levels of variable indicators associated with mortality, we also performed ROC analysis. The results indicate that the critical thresholds are: SpO₂ below 78% (sensitivity - 54.0%, specificity - 87.3%, p < 0.0001), SBP below 80 mmHg (sensitivity - 84.0%, specificity - 54.7%, p < 0.0001), and LVEF below 52% according to echocardiography (sensitivity - 66.0%, specificity - 82.5%, p < 0.0001) (Fig.1).
Fig 1: ROC Curve for the Predictive Formula of Logistic Regression to Determine (Positive or Negative) Prognosis.
Based on our results, we created a formula for risk stratification of in-hospital mortality for pts with high-risk PE:
Y=(exp(17,47-0,11*SpO2+1,58* Thrombus in bifurcation of PA-0,09*SBP
+1,66*female sex -0,08*LV EF))/(1+exp (17,47-0.11*SpO2+1.58* Thrombus in bifurcation of PA -0,09*SBP +1,66*female sex -0,08*LV EF)).
Calculating the individual prognosis for pts with high-risk PE, the additional risk of in-hospital mortality is higher if the value of Y > 0.5, and does not increase when Y < 0.5.
Using ROC analysis, it was established that the sensitivity of the obtained prognostic formula is 64.0%; specificity is 92.1% (area under the ROC curve is 0.7861; 95% confidence interval 0.784-0.93, p < 0.0001).
Based on the analysis and verification of the formula, we proposed the SBAFS scale (S – saturation, B – bifurcation, A – arterial hypotension, F – LV ejection fraction, S – female sex) for calculating the additional risk of in-hospital mortality in high-risk patients with PE (Table 3).
Table 3: Scale for Determining Additional Mortality Risk in High-Risk Pulmonary Embolism Patients (SBAFS)
Parameter |
Point |
SpO2<78% |
1 – yes 0 - no |
Bifurcation lesion |
1 – yes 0 - no |
Arterial Hypotension (SBP<80 mmHg) |
1 – yes 0 - no |
LV EF<52 % |
1 – yes 0 - no |
Sex |
|
Note. LV EF – ejection fraction of left ventricle, SpO2 – blood oxygen saturation.
In order to verify and check our formula we organized the second phase of our study based on a multicenter cohort. Patients were divided into two subgroups: I - 65 (80%) pts who were discharged with improved status; II - 15 (20%) pts who died during hospitalization (Table 4). The SBAFS scale was evaluated in Phase 2.
Table 4: Indicators of the SBAFS Scale in High-Risk Pulmonary Embolism Patients Based on Prognosis (Stage II of the Study).
Показник
|
All (n=75) |
І subgroup (n=60) |
ІІ subgroup (n=15) |
M-U, p |
Age (М ± σ), years |
60,60±13,40 |
59.67±13.61 |
64.44±12.64 |
0,162 |
Men, n (%) |
38 (50,7) |
31 (51,7) |
7 (46,7) |
0,954 |
Women, n (%) |
37 (49,3) |
29 (48,3) |
8 (53,3) |
|
Women, n (%) |
37 (49,3) |
29 (48,3) |
8 (53,3) |
0,954 |
SpO2 (М ± σ), % |
82,92±9,73 |
85,02±7.18 |
74,21±13,75 |
0,003 |
SpO2<78 %, n (%) |
15 (20) |
8 (13,3) |
7 (46,7) |
0,01 |
SBP (М ± σ), mm Hg. |
86,83±20,58 |
89,53±18.77 |
76,00±24,44 |
0,011 |
SBP<80 mm Hg, n (%) |
18 (24,0) |
12 (20,0) |
6 (40,0) |
0,199 |
LV EF by Echo (М ± σ), % |
54,27±11,02 |
56,24±9.71 |
43,73±12,08 |
0,002 |
LV EF<52% by Echo, n (%) |
22 (29,3) |
15 (25,0 |
7 (46,7) |
0,183 |
Bifurcation lesion by CTPA, n (%) |
17 (22,7) |
14 (23,3) |
3 (20,0) |
0,543 |
Mean value of Y according to the formula |
0,401±0,369 |
0,306±0,330 |
0,785±0,250 |
0,00002 |
Value of Y > 0,5, n (%) |
29 (38,7) |
16 (26,7) |
13 (86,7) |
0,00001 |
SBAFS score |
||||
0 |
5 (6,7) |
5 (8,3) |
0 |
0,316 |
1
|
28 (37,3) |
28 (46,7) |
0 |
0,0003 |
2 |
27 (36,0) |
18 (30,0) |
9 (60,0) |
0,088 |
3 |
11(14,7) |
6 (10,0) |
5 (33,3) |
0,037 |
4 |
3 (4,0) |
3 (5,0) |
0 |
0,507 |
5 |
1(1,3) |
0 |
1 (6,7) |
0,200 |
Mean score |
1,76±1,01 |
1,57±0,96 |
2,53±0,83 |
0,0006 |
|
|
|
|
|
Note. SBP – systolic blood pressure, SpO2 – blood oxygen saturation.
As it was shown in the first stage of our study, SBP, LVEF and SpO2 were significantly lower in subgroup II pts compared to subgroup I (Table 5).
When verifying the formula, it was established that in subgroup II pts, the average value of the indicator Y was significantly higher than in subgroup I: (0.785 ± 0.250) vs. (0.306 ± 0.330), p = 0.00002. Additionally, the proportion of pts with Y > 0.5 was also significantly higher: 13 (86.7%) vs. 16 (26.7%), respectively, p = 0.00001, confirming the presence of an additional risk of early death in subgroup II pts.
To assess the sensitivity and specificity of the proposed formula, we performed ROC analysis for the Phase 2 of our study (Figure 2). It was shown that in the multicenter cohort, the sensitivity of our prognostic formula was 93.3%, and the specificity was 74.6% (area under the ROC curve was 0.803; 95% confidence interval 0.690-0.889, p < 0.001). Thus, during the verification of the formula in the second phase of our study, the sensitivity and specificity not only reproduced the results from the first stage but even exceeded the previous findings. Additionally, the sensitivity and specificity of indicators such as SBP, LVEF and SpO2, as evaluated by ROC analysis conducted in both stages of the study, were significant and comparable.
Fig 2: ROC Curve for the Predictive Formula of Logistic Regression to Determine (Positive or Negative) Prognosis (Phase2).
While evaluating the proposed SBAFS scale, we also found that the average score for subgroup II on the scale was (2.53 ± 0.83) and was significantly higher than in subgroup I: (1.57 ± 0.96), (p = 0.0006). In subgroup II there was also a significantly greater proportion of pts who had more than 3 points according to SBAFS.
As we mentioned before lower levels of SpO2, SBP and bifurcational lesions verified by CTPA are associated with increased in-hospital mortality in pts with high-risk PE. This statement is logical as those factors are directly pathophysiologically linked to the hemodynamic overload and dysfunction of the right ventricle. Moreover, these factors are also used in various scales for predicting the risk of death in patients with acute PE and are consistent with the results of other studies [3, 5, 13].
Taking into account female sex according to our previous studies (as well as other scientific works), it is not an additional risk factor in patients with non-high-risk PE [7, 14]. However, in pts with high-risk PE, female sex may significantly serve as an additional factor for unfavorable short-term prognosis, warranting further investigation. Analysis of registries has shown that high-risk PE occurs more frequently in women than in men (14.6% vs. 9.2%, respectively; p = 0.0002) [4]. Possible explanations for these findings may include sex differences in risk factors for VTE and anatomical characteristics of the cardiovascular and respiratory systems, leading to more significant right ventricular decompensation in women during PE [15, 16].
Thus, the results of our study suggest that, in addition to known factors associated with unfavorable 30-day prognosis in high-risk PE patients, additional factors associated with in-hospital mortality in such patients are: SpO2 below 78%, SBP below 80 mmHg, LVEF below 52%, saddle PE verified by CTPA, female sex. The formula and SBAFS scale, which we developed, can be used for determining additional mortality risk in pts with high-risk PE based in Ukraine.
Limitations: In this study, invasive treatment methods for patients with pulmonary embolism were not applied.
Conflict of Interest: The authors declare that there are no conflicts of interest.
Author Contributions: Concept and study design, article editing, conclusion formulation – Vira Tseluyko; statistical data analysis – L. Yakovleva; article writing – L. Yakovleva, N. Mishchuk, M. Kurinna, R. Askierov; data collection – M. Kurinna, L. Kharchenko, R. Askierov.