Despite the existence of several risk scores, the accurate prediction of the prognosis in pulmonary embolism (PE) remains a challenge. The Global Registry of Acute Coronary Events (GRACE) risk score has a high diagnostic performance for adverse outcomes in acute coronary syndrome. We aimed to assess the applicability and extend the use of the GRACE risk score to PE. A case-control study of 206 consecutive patients admitted with PE was performed. The GRACE, Geneva, Simplified Pulmonary Embolism Severity Index, Shock Index, and European Society of Cardiology risk scores were tested for the prediction of the primary end point: all-cause 30-day mortality. Comparisons between GRACE and the other risk scores were performed using receiver operating characteristic area under the curve and the integrated discrimination improvement index. All-cause 30-day mortality was observed in 18.9% of the patients. Unlike the other classifications, no adverse outcomes were observed in patients classified as low risk using the GRACE risk score (100% negative predictive value for GRACE risk score ≤113). The GRACE score showed greater discriminative performance than the Geneva score (area under the curve 0.623, 95% confidence interval [CI] 0.53 to 0.71), Shock Index (area under the curve 0.639, 95% CI 0.55 to 0.73), European Society of Cardiology (area under the curve 0.662, 95% CI 0.57 to 0.76), and Simplified Pulmonary Embolism Severity Index (area under the curve 0.705, 95% CI 0.61 to 0.80), although statistical significance was not reached. The integrated discrimination improvement index suggested a more appropriate risk classification with the GRACE score. In conclusion, our results have demonstrated that the GRACE risk score can accurately predict 30-day mortality in patients admitted for acute PE. Compared to previously proposed PE prediction rules, the GRACE risk score presented improved overall risk classification.
Acute pulmonary embolism (PE) risk stratification should reliably identify patients at high risk of death who could receive specific therapeutic interventions and those suitable for outpatient management. At present, the most used PE risk scores are the clinical prediction rule proposed by the European Society of Cardiology (ESC), Geneva prognostic score, Simplified Pulmonary Embolism Severity Index (sPESI), and the Shock Index. The Global Registry of Acute Coronary Events (GRACE) score is a prognostic assessment tool used in acute coronary syndromes, consisting of clinical, analytical, and electrocardiographic variables, making it an interesting assessment model for acute cardiovascular diseases. The GRACE score has never been tested for PE prognosis assessment. Our aim was to establish the appropriateness of the GRACE score in this setting and to compare it with established PE risk scores.
Methods
We performed a case-control study of 206 consecutive patients (mean age 70.3 ± 15.6 years, 41.3% men) with an “International Diseases Classification, 10th revision,” discharge code of PE from January 2007 to December 2010. A clinical review of each case was performed to ensure that the “International Diseases Classification, 10th revision,” coding referred to acute PE. Baseline data were collected at the PE diagnosis, and the prognostic scores (Geneva prognostic score, sPESI, ESC, Shock Index, and GRACE) were assessed retrospectively. The primary outcome measure was all-cause 30-day mortality assessed through hospital medical records revision or proxy interviews. No patient was lost to follow-up (30 days from presentation). The human subjects committee of Coimbra’s Hospital Centre and University approved the study protocol.
All patients were recruited from the emergency department of a tertiary care hospital, and all patients provided informed consent. Eligibility for the present study required that patients had acute symptomatic PE confirmed by multidetector computed tomography. Patients with PE during hospitalization for other clinical reasons were not included.
The creatinine clearance was estimated using the Modification of Diet in Renal Disease equation: glomerular filtration rate (ml/min/1.73 m 2 ) = 186 × (plasma creatinine) − 1.154 × (age) − 0.203 × (0.742 if female gender) × (1.210 if Afro-Caribbean). Troponin I was considered elevated if >0.034 ng/ml. Right ventricular failure was defined as an increase in the right ventricular diameter and/or the occurrence of right ventricular pressure overload in the echocardiographic assessment or N-terminal pro-B-type natriuretic peptide >4,500 pg/ml.
The patients were risk classified according to previously reported cutoff points. Regarding the Geneva prognostic score, the patients were classified as low (score ≤2) or high (score >2) risk. The sPESI dichotomized subjects into low (score <1) and high (score ≥1) risk. The ESC prediction rule stratified patients as low risk if no clinical or right ventricular dysfunction/myocardial injury markers were present, intermediate if either right ventricular dysfunction and/or myocardial injury markers were present, and high risk if shock or persistent hypotension were documented. The Shock Index, defined as the heart rate divided by the systolic blood pressure, divided patients into low risk (score <1) and high risk (score ≥1). For the GRACE score, we distributed patients into 3 risk strata using the criteria recommended for acute coronary syndromes (low risk, score <109; intermediate risk, score 109 to 140; and high risk, score >140) using the cutoff points of the receiver operating characteristic curve that ensured the most favorable compromise between sensitivity and specificity. This score was termed the “PE risk score.”
Statistical analysis was performed using SPSS, version 17.0 (SPSS, Chicago, Illinois) for the comparison of nominal variables with the chi-square test and continual variables with the Student’s t test. Univariate analysis was performed with the chi-square test, determining the respective relative risk with the 95% confidence interval (CI). Results with p <0.05 were regarded as statistically significant. We assessed the overall discriminatory power of each risk score (Geneva prognostic score, sPESI, ESC, Shock Index, and GRACE) by calculating the area under each receiver operating characteristic curve. Discrimination, measured in terms of the area under the curve, refers to a model’s ability to assign a greater probability to nonsurvivors than to survivors. However, it has been widely reported that the area under the curve cannot always measure a clinically meaningful quantity or demonstrate the value of new markers or models useful for prediction. Recently, other measures of incremental value have been proposed that examine the extent to which a new marker or model reclassifies subjects, such as the net reclassification improvement index and the integrated discrimination improvement. The net reclassification improvement method, described by Pencina and D’Agostino, states that a positive and significant net reclassification improvement translates a net overall successful reclassification of subjects into a more appropriate risk category. The integrated discrimination improvement, which can be seen as a continuous form of the net reclassification improvement method, assesses the improvement in risk discrimination by estimating the change in the difference of the average of the predicted probabilities of an event between those with and without the event under consideration. It is a more appropriate measure of risk reclassification when comparing scores with considerably different risk categorization (e.g., sPESI stratifies patients into 2 risk strata and GRACE into 3 risk strata). The GRACE score was also tested for calibration. Calibration refers to the agreement between the predicted and true probabilities and is most often measured using Hosmer-Lemeshow goodness-of-fit statistics, suggesting a good fit when associated with p >0.05.
Comparisons using receiver operating characteristics curve analysis and a measure of risk reclassification (the integrated discrimination improvement index) were performed between the GRACE model and the currently available PE risk scores (Geneva prognostic score, Shock Index, sPESI, and ESC). Furthermore, GRACE goodness of fit was assessed with the modified Hosmer-Lemeshow chi-square statistic.
Results
The patients’ baseline clinical and analytical and imaging characteristics are listed in Table 1 . In our population sample, 39 patients (18.9%) died within 30 days of presentation, with 33 corresponding to in-hospital mortality. No patient was lost to follow-up.
Variable | 30-d Mortality | p Value | ||
---|---|---|---|---|
Overall (n = 206) | No (n = 167) | Yes (n = 39) | ||
Age (yrs) | 70.27 ± 15.59 | 68.10 ± 16.25 | 74.74 ± 11.89 | 0.040 |
Age >80 yrs | 55 (26.7%) | 39 (23.4%) | 16 (41.0%) | 0.055 |
Men | 85 (41.3%) | 69 (41.3%) | 16 (41.0%) | 0.973 |
Hypertension | 91 (44.2%) | 69 (41.3%) | 22 (56.0%) | 0.084 |
Diabetes mellitus | 38 (18.4%) | 27 (16.2%) | 11 (28.0%) | 0.082 |
Heart failure and/or chronic lung disease | 58 (28.2%) | 41 (24.6%) | 17 (44.0%) | 0.041 |
Active or past malignancy | 42 (20.4%) | 29 (17.4%) | 13 (33.0%) | 0.050 |
Systolic arterial pressure (mm Hg) | 122.62 ± 23.36 | 123.61 ± 21.99 | 119.13 ± 26.19 | 0.286 |
Systolic blood pressure <100 mm Hg | 25 (12.1%) | 17 (10.2%) | 8 (21.0%) | 0.151 |
Heart rate (beats/min) | 101.29 ± 21.30 | 98.47 ± 20.75 | 110.76 ± 19.80 | 0.001 |
Heart rate ≥110 beats/min | 62 (30.1%) | 42 (25.1%) | 20 (51.0%) | 0.007 |
Respiratory frequency | 26.43 ± 8.67 | 24.67 ± 7.81 | 29.75 ± 8.53 | 0.017 |
Respiratory rate >30 breaths/min | 25 (12.1%) | 24.92 ± 8.0 | 29.75 ± 8.5 | 0.024 |
Arterial oxyhemoglobin saturation (%) | 90.49 ± 6.66 | 91.06 ± 6.5 | 88.35 ± 7.5 | 0.024 |
Arterial oxyhemoglobin saturation <90% | 59 (28.6%) | 41 (24,6%) | 18 (46.0%) | 0.020 |
Killip-Kimball class (admission) | 1.38 ± 0.65 | 1.39 ± 0.71 | 1.87 ± 1.11 | 0.007 |
Shock | 15 (7.3%) | 8 (4.8%) | 7 (18.0%) | 0.012 |
Thrombolysis | 19 (9.2%) | 13 (7.8%) | 6 (15.0%) | 0.008 |
Mechanical ventilation | 7 (3.4%) | 3 (1.8%) | 4 (10.0%) | 0.210 |
Inotropic therapy | 14 (6.8%) | 8 (4.8%) | 6 (15.0%) | 0.290 |
Admission hemoglobin (g/dl) | 12.83 ± 2.0 | 12.99 ± 1.99 | 12.08 ± 1.92 | 0.013 |
Admission C-reactive protein (mg/dl) | 7.65 ± 8.33 | 6.58 ± 7.33 | 10.54 ± 9.99 | 0.032 |
Admission glycemia (mmol/L) | 8.23 ± 4.99 | 7.89 ± 5.12 | 9.20 ± 5.03 | 0.042 |
D-dimers (μg/ml) | 4,758.69 ± 8,584.50 | 4,603 ± 8,231 | 5,920 ± 11,170 | 0.386 |
Maximum troponin I (ng/ml) | 0.56 ± 1.64 | 0.49 ± 1.63 | 0.86 ± 1.83 | 0.015 |
Troponin I >0.034 ng/ml | 113 (54.9%) | 83 (49.7%) | 30 (77.0%) | 0.023 |
Natriuretic pro-B-type natriuretic peptide (pg/ml) | 5,704.32 ± 11,592.98 | 3,937 ± 5,369 | 13,652 ± 21,942 | 0.001 |
Admission creatinine (μmol/L) | 101.67 ± 53.14 | 96.09 ± 51.98 | 123.59 ± 57.99 | <0.001 |
Creatinine clearance ∗ (ml/min) | 68.42 ± 28.30 | 72.59 ± 27.17 | 55.59 ± 29.31 | 0.001 |
Right ventricle dysfunction | 39 (18.9%) | 26 (15.6%) | 13 (33.0%) | 0.018 |
ST-segment elevation | 3 (1.5%) | 2 (1.2%) | 1 (3.0%) | 0.277 |
∗ Creatinine clearance determined according to Modified Diet in Renal Disease.
Hemodynamic shock was documented in 15 patients (7.3%), and no cases of cardiac arrest at admission were reported. Regarding specific medical interventions, thrombolysis was used in 19 patients (9.2%), mechanical ventilation in 7 (3.4%), and inotropic support in 14 (6.8%). Major hemorrhagic complications after thrombolysis were seen in 4 patients. In 2 of those patients, bleeding was considered related to the cause of death. A group comparison according to the occurrence of all-cause 30-day mortality is listed in Table 1 , and the predictors of death on univariate analysis (categorical variables) are listed in Table 2 .
Variable | OR | 95% CI | p Value |
---|---|---|---|
Age >80 yrs | 2.03 | 0.98–4.24 | 0.055 |
Heart failure and/or chronic lung disease | 2.11 | 1.02–4.37 | 0.041 |
Active or past malignancy | 2.14 | 1.00–4.66 | 0.050 |
Heart rate ≥110 beats/min | 2.70 | 1.30–5.61 | 0.007 |
Respiratory rate >30 breaths/min | 3.27 | 1.14–9.34 | 0.023 |
Arterial oxyhemoglobin saturation <90% | 2.38 | 1.14–4.99 | 0.020 |
Shock | 3.75 | 1.27–11.09 | 0.012 |
Troponin I >0.034 ng/ml | 3.11 | 1.13–8.59 | 0.023 |
Right ventricular dysfunction | 2.54 | 1.16–5.58 | 0.036 |
Using the previously defined criteria and the previously described scores, we risk stratified our cohort ( Table 3 ). The Geneva prognostic score and Shock Index were the scores with a greater proportion of patients in the low-risk group (score ≤ 2; Geneva prognostic score, 161 of 206 patients [78.2%]; Shock Index, 153 of 206 patients [74.3%]) compared to the sPESI (score <1; 46 of 206 patients [22.3%]) and ESC (low-risk class; 54 of 206 patients [26.2%]). The primary end point was reached in the Geneva prognostic score low-risk group (69.2% of deaths; 27 of 39 deaths), in Shock Index low-risk group (56.4% of deaths; 22 of 39 deaths), and in the sPESI and ESC low-risk groups (10.3%, 4 of 39 deaths, and 7.7%, 3 of 39 deaths, respectively).
Variable | Geneva ∗ | Shock Index † | sPESI ‡ | ESC § | ACS score ¶ (GRACE) | PE Risk Score || (GRACE) |
---|---|---|---|---|---|---|
Low risk | 161 (78.2%) | 153 (74.3%) | 46 (22.3%) | 54 (26.2%) | 39 (18.9%) | 45 (21.8%) |
Intermediate risk | — | — | — | 129 (62.6%) | 50 (24.3%) | 60 (29.1%) |
High risk | 45 (21.8%) | 53 (25.7%) | 160 (77.7%) | 23 (11.2%) | 117 (56.8%) | 101 (49.0%) |
Mortality | ||||||
Low-risk group | 27 (69.2%) | 22 (56.4%) | 4 (10.3%) | 3 (7.7%) | 0 (0%) | 0 (0%) |
Intermediate-risk group | — | — | — | 29 (74.4%) | 9 (23.1%) | 9 (23.1%) |
High-risk group | 12 (30.8%) | 17 (43.6%) | 35 (89.7%) | 7 (17.9%) | 30 (76.9%) | 30 (76.9%) |
Sensitivity | 30.8% (17.5–47.7) | 43.6% (28.2–60.2) | 89.7% (74.8–96.7) | 92.3% (78.0–98.0) | 100% (88.8–100.0) | 100% (88.8–100.0) |
Specificity | 80.2% (73.2–85.8) | 78.4% (71.3–84.3) | 25.1% (18.9–32.6) | 30.5% (23.8–38.2) | 23.4% (17.3–30.6) | 26.9% (20.5–34.4) |
Positive predictive value | 26.7% | 32.1% | 21.9% | 23.7% | 23.4% | 24.2% |
Negative predictive value | 83.2% | 85.6% | 91.3% | 94.4% | 100% | 100% |
∗ Geneva prognostic score >2 designates high-risk category, with lower values indicating low risk.
† Shock Index score of ≥1 designates high-risk category, with lower values indicating low risk.
‡ sPESI score of ≥1 designates high-risk category, with lower values indicating low risk.
§ ESC risk model stratifies into high risk (presence of shock or persistent hypotension), intermediate risk (presence of myocardial injury markers), and low risk (none of these markers).
¶ Criteria recommended for acute coronary syndromes (ACS score ) >140 designated as high risk; 109–140 as intermediate risk; and <109 as low risk.
|| PE risk score (GRACE): >153 designated high risk; 114–153 designated intermediate risk; and <114 designated low risk.
The receiver operating characteristic curve and the area under the curve for the 4 algorithms for all-cause 30-day mortality prediction are shown in the Figure 1 and listed in Table 4 . The receiver operating characteristic curve for the Geneva prognostic score (area under the curve 0.623, 95% CI 0.53 to 0.71; p = 0.016), receiver operating characteristic curve for Shock Index (area under the curve 0.639, 95% CI 0.55 to 0.73; p = 0.008), receiver operating characteristic curve for sPESI (area under the curve 0.705, 95% CI 0.61 to 0.80; p <0.001), and receiver operating characteristic curve for ESC (area under the curve 0.662, 95% CI 0.57 to 0.76; p = 0.002). No statistically significant differences were found between the area under the curve for the 4 scores (p >0.05; Table 4 ). The performance of all risk scores in our cohort is listed in Table 3 . We used the stratification criteria proposed by the original investigators: Geneva prognostic score >2, Shock Index ≥1, sPESI ≥1, and ESC (intermediate or more). ESC’s prediction rule had greater sensitivity and a greater negative predictive value than the remaining PE scores for predicting all-cause 30-day mortality.
ROC AUC for Prediction of All-cause 30-d Mortality | |||||
---|---|---|---|---|---|
Geneva | Shock Index | sPESI | ESC | GRACE | |
Area under curve | 0.623 | 0.639 | 0.705 | 0.662 | 0.715 |
95% Confidence interval | 0.53–0.71 | 0.55–0.73 | 0.61–0.80 | 0.57–0.76 | 0.63–0.80 |
p Value | 0.016 | 0.008 | <0.001 | 0.002 | <0.001 |
Cross-tabulation of area under receiver operating characteristic curve among different scores | |||||
Geneva | Geneva | Area differences = 0.023; p = 0.721 | Area differences = 0.066; p = 0.181 | Area differences = 0.002; p = 0.980 | Area differences = 0.076; p = 0.271 |
Shock Index | Area differences = 0.023; p = 0.721 | Shock Index | Area differences = 0.043; p = 0.454 | Area differences = 0.021; p = 0.752 | Area differences = 0.053; p = 0.351 |
Area differences = 0.066; p = 0.181 | Area differences = 0.043; p = 0.454 | sPESI | Area differences = 0.064; p = 0.353 | Area differences = 0.010; p = 0.860 | |
European Society of Cardiology | Area differences = 0.002; p = 0.980 | Area differences = 0.021 P value = 0.752 | Area differences = 0.064 P value = 0.353 | ESC | Area differences = 0.074 P value = 0.212 |
Global Registry of Acute Coronary Events | Area differences = 0.076; p = 0.271 | Area differences = 0.053; p = 0.351 | Area differences = 0.010; p = 0.860 | Area differences = 0.074; p = 0.212 | GRACE |