Cardiovascular Risk Assessment of Pulmonary Embolism With the GRACE Risk Score




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.



Table 1

Sample characteristics and sub-analysis according to the presence of 30-day mortality






































































































































































































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 .



Table 2

Predictors of 30-day mortality on univariate analysis (categorical variables)






















































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

OR = odds ratio.


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).



Table 3

Patient classification and receiver operator characteristic test characteristics for prediction of 30-day mortality




































































































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.




Figure 1


PE scores and GRACE receiver operating characteristic curves. Gray line indicates ESC; blue line, Geneva risk score; red line, GRACE risk score; yellow line, Shock Index; and green line, sPESI.


Table 4

Receiver operating characteristic curve analysis for prediction of 30-day mortality using each risk score












































































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

Only gold members can continue reading. Log In or Register to continue

Stay updated, free articles. Join our Telegram channel

Dec 7, 2016 | Posted by in CARDIOLOGY | Comments Off on Cardiovascular Risk Assessment of Pulmonary Embolism With the GRACE Risk Score

Full access? Get Clinical Tree

Get Clinical Tree app for offline access