Abstract
Introduction:
Glomerular filtration rate (GFR) and blood urea nitrogen (BUN) are important prognostic indicators for cardiovascular disease. However, data on the relationship between renal dysfunction (RD) and prognosis in patients with acute pulmonary embolism (APE) are limited. The estimated-GFR (eGFR), based on the Modification of Diet in Renal Disease (MDRD) equation, has been suggested as a possible prognostic marker in patients with APE; however, at present, the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation is thought to be more accurate than the MDRD equation for the estimation of RD.
Objective:
We investigated whether eGFRCKD-EPI or BUN could predict adverse outcomes (AOs) better than eGFRMDRD in normotensive patients with APE.
Methods:
Ninety-nine normotensive patients with APE (aged 22-96, 56% male) were enrolled in the study retrospectively. Adverse outcomes were defined as the occurrence of any of the following: death, cardiopulmonary resuscitation, use of vasopressors, thrombolysis, or mechanical ventilation.
Results:
In univariate analyses, age, gender (male), heart rate (>110 bpm), serum creatinine, BUN, cardiac troponin (cTn) positivity, right ventricle–left ventricle ratio, eGFRMDRD, and eGFRCKD-EPI were found to be significantly different between those with and without AOs. Comparing area under the curves for AO, we found statistically significant differences between eGFRCKD-EPI and eGFRMDRD (P = .01) but not between BUN and eGFRCKD-EPI or BUN and eGFRMDRD. Furthermore, 30-day mortality was 36% versus 11% in cTn-positive patients with an eGFRCKD-EPI < and ≥ 60 mL/min, respectively.
Conclusion:
There is a close relationship between RD and APE prognosis. We conclude eGFRCKD-EPI is a potential prognostic marker for risk stratification in normotensive patients with APE.
Introduction
In acute pulmonary embolism (APE), hemodynamic instability occurring due to an abrupt increase in pulmonary vascular resistance (PVR) is closely associated with adverse events. Increased PVR results in backward right heart dysfunction, increased central venous pressure (CVP), decreased pulmonary blood flow, and hypotension. Both pulmonary and systemic hemodynamics are maintained at normal levels by increased neurohumoral activation (NHA). 1
Glomerular filtration rate (GFR) and blood urea nitrogen (BUN) are routinely used as indicators of renal function. Glomerular filtration rate is not only closely associated with preexisting renal pathology but also with hemodynamic alterations. Several formulas have been developed to estimate GFR because reference methods regarding its measurement are expensive, time consuming, and not readily available in many institutions. The most popular method for GFR estimation is the simplified Modification of Diet in Renal Disease (MDRD) equation, which takes into account serum creatinine (sCr) levels, age, race, and gender and is reported in mL/min/1.73 m2. However, this formula seems to systematically underestimate GFR in patients with normal renal function. Therefore, the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation has been developed, which is reported to be more accurate than the MDRD formula. 2 Blood urea nitrogen is a less specific measure of renal dysfunction (RD) compared to GFR. Blood urea nitrogen concentration is also affected by hemodynamics, and in particular neurohumoral mediators, and recent evidence suggests that it is a better prognostic marker for heart failure (HF) than is GFR. 3,4
Renal dysfunction is an important prognostic factor for cardiovascular disease (CVD) and is associated with increased mortality and morbidity, 5 but data on the relationship between RD and APE prognosis are limited to a few studies. The use of GFR to improve APE risk stratification has been suggested; however, previous studies have mostly used the MDRD equation to calculate the estimated GFR (eGFR). 6 –8 In this study, we assessed whether eGFRCKD-EPI or BUN is a better prognostic factor for APE compared to eGFRMDRD.
Materials and Methods
Patients and Study Design
This cross-sectional, retrospective observational study was conducted at the Department of Pulmonary Medicine at Bulent Ecevit University Hospital, Zonguldak, Turkey. The study was approved by the Ethics Review Board of Bulent Ecevit University. Patients with APE were identified from January 1, 2013, to December 31, 2013, using an electronic patient database and patient charts. One hundred and nine (n = 109) patients with symptomatic APE were diagnosed and treated during this period (consecutively). Of the 109 patients, 10 were excluded from the study due to high-risk APE (systolic blood pressure < 90 mm Hg on admission, n = 3) and uncertain diagnosis (intermediate probability ventilation/perfusion [V/Q] scan, n = 7). Ultimately, the study population consisted of 99 patients. The APE diagnosis was based on computed tomography pulmonary angiography (CTPA; n = 79), high probability V/Q scan (n = 18), and compression sonography of the leg veins with high clinical probability (n = 2). Chronic obstructive pulmonary disease, asthma, interstitial lung diseases, and pneumoconiosis were regarded as chronic lung diseases (CLDs). Heart failure was determined based on less than 50% ejection fraction of the left ventricle and/or clinically overt HF. Stroke (ischemic or hemorrhagic) or paresis was regarded as cerebrovascular accident. Cardiovascular disease was defined as the presence of at least one of the following conditions: ischemic heart disease, HF, stroke, or chronic kidney disease (CKD). Adverse clinical outcome was defined as the occurrence of any of the following: death, cardiopulmonary resuscitation, use of vasopressors, thrombolysis, or mechanical ventilation. Three senior physicians determined 30-day mortality and adverse outcome after carefully reviewing all available data. PE-related mortality was defined as death after cardiopulmonary resuscitation, use of vasopressors, thrombolysis, or need for mechanical ventilation.
Laboratory Analyses
On admission, venous blood samples were collected for creatinine, BUN, and troponin concentrations. All measurements were performed before the diagnostic procedures. Serum creatinine levels were analyzed using the Jaffe method and an Advia 2400 analyzer (Siemens Diagnostics, Tarrytown, New York). The BUN levels were measured using a biuret reagent and an Advia 2400 analyzer.
Troponin I and T were measured at different times (before and after June 2013, respectively). Therefore, statistical analysis was based on troponin positivity instead of continuous values. Troponin T levels were measured by electrochemiluminescence using the Elecsys e411 analyzer (Roche, Mannheim, Germany). Troponin I levels were measured by chemiluminescence using the Advia Centaur CP (Siemens, Munich, Germany).
Estimation of GFR
The equations used for GFR estimation are summarized in Table 1. The National Kidney Disease Education Program recommends reporting categorical (simply as ≥60 mL/min/1.73 m2) instead of exact numbers for eGFR values ≥ 60 mL/min/1.73 m2, based on known inaccuracies for GFR values > 60 mL/min/1.73 m2. 9 Therefore, 60 mL/min/1.73 m2 was used as the cutoff for normal eGFR in the statistical analysis.
Equations Utilized in Study.a
Abbreviations: min, the minimum of Scr/κ or 1; max, the maximum of Scr/κ or 1; Scr, serum creatinine (mg/dL), GFR, glomerular filtration rate, BUN, blood urea nitrogen, CKD-EPI, Chronic Kidney Disease Epidemiology Collaboration; MDRD, Modification of Diet in Renal Disease.
aκ = .7, α = −.329 for females, κ = .9, α = −.411 for males.
Computed Tomography Pulmonary Angiography
A multidetector CT system (Activision 16-row CT scanner; Toshiba Medical Systems, Otawara, Japan) was used for CT imaging. A total of 100 mL nonionic contrast agent (Ultravist 370; Bayer Schering Pharma, Berlin, Germany) were given at a rate of 3.0 mL/s via a peripheral venous line. Routine CTPA protocols were used for all patients, and the parameters were 120 kV, 144 effective mAs, a pitch factor of 0.938, a helical factor of 15.0, a rotation time of 0.75 s, and a reconstruction interval of 1 mm. Automatic bolus tracking was performed in the pulmonary trunk with a trigger of 120 Hounsfield units.
For each patient, CT imaging analysis was performed by the radiologist (I.O.) using a multimodality workstation (Infinitt PACS version 3.0.9.1; Infinitt Co, Ltd, Seoul, Korea). A 4-chamber heart (4-CH) view was constructed as described previously. 10 The maximum diameter of the right ventricle (RV) and left ventricle (LV) was measured on a reconstructed 4-CH view between the inner aspect of the interventricular septum and the ventricular endocardium, perpendicular to the long axis of the heart. In 2 patients, the pulmonary arteries were not visualized adequately, but the RV–LV ratio and superior vena cava (SVC) diameter could be measured. In those patients, diagnosis was confirmed by V/Q scintigraphy.
The pulmonary vascular obstruction index (PAOI) was calculated as described previously by Qanadli et al. 11 With the Qanadli index, each lung has 10 segmental artery branches (3 to the upper lobes, 2 to both the middle and lingual lobes, and 5 to the lower lobes). The number of obstructed segmental arteries is corrected by a factor of 1 for partial obstruction and by a factor of 2 for complete obstruction. The maximal obstruction index is 40 per patient, which is equivalent to 100% obstruction.
Statistical Analysis
Descriptive statistics of the categorical variables are given as numbers or percentages; continuous variables are given as means (standard deviations) or medians (min-max). A chi-square test was used to evaluate categorical variables. The Student t test or Mann-Whitney U test was used to compare the means/medians of variables where appropriate. Receiver–operating characteristic (ROC) curves were used to evaluate the predictive power of each eGFR calculation method on adverse clinical events. Statistical comparisons of ROC curves are based on the methods of DeLong et al. 12 Predictors of adverse clinical events among patients with APE were identified using univariate and multivariate logistic regression analyses. Statistical analysis was performed using SPSS version 18.0 for Windows (SPSS, IBM Inc, Chicago, Illinois) and MedCalc for Windows, version 12.2.1.0 (MedCalc Software, Ostend, Belgium). P values were 2-sided, and values less than .05 were considered statistically significant. The post hoc power of the study was found to be 0.86 with an effect size of (n1 = 22, p1 = 0.72; n2 = 77, p2 = 0.36) and α 0.05 by Power and Precision Trial 4.1.0.
Results
Baseline Characteristics and Clinical Course
The study population consisted of 44 females and 55 males with a median age of 68 (range: 22-96) years. Overall, 22 (22%) patients had an adverse outcome. In patients with intermediate-risk APE (n = 83), the mortality rate was 18%. 1 No deaths were observed in patients with low-risk APE (n = 16). The all-cause 30-day mortality rate was 15% (n = 14), and median survival time was 17 (range 2-30) days. APE-related mortality was 10% (n = 10); the remaining 4 patients died due to pneumonia (n = 1), COPD (n = 1), brain and hepatic abscesses (n = 1), and septic shock (n = 1). All patients were given standard anticoagulant therapy with intravenous unfractionated heparin (UFH) or a subcutaneous low-molecular-weight heparin (body mass adjusted). Massive hemorrhage (gross hematuria) caused by UFH was observed in 1 patient who died because of APE-related respiratory failure. No hemorrhage-associated deaths were observed. Eight patients received thrombolytic treatment. Hemodynamic deterioration in patients who presented as hemodynamically stable was considered an indication for thrombolytic treatment. The mortality rate was 38% (n = 3) in the group receiving thrombolytic treatment. The values of eGFRMDRD (<60 mL/min) and eGFRCKD- EPI (<60 mL/min) were noted in 42% and 49% of patients with intermediate risk and in 19% and 19% of patients with low-risk APE, respectively. Significant differences between patients with eGFRCKD-EPI < and ≥ 60 mL/min/1.73 m2 with respect to APE-related death (21% vs 4%, P = .008) and thrombolytic therapy (17% vs 2%, P = .008) were observed, but no differences were found with respect to all-cause 30-day mortality (21% vs 9%, P = .107). No significant differences were observed between patients with eGFRMDRD < and > 60 mL/min/1.73 m2 with respect to APE-related death (16% vs 8%, P = .242) and all-cause 30-day mortality (19% vs 12%, P = .335), but thrombolytic therapy was significantly different between the groups (19% vs 2%, P = .002). Clinical characteristics are presented in Table 2.
Baseline Characteristics.
Abbreviations: SBP, systolic blood pressure; BUN, blood urea nitrogen; eGFR, estimated glomerular filtration rate; CVD, cardiovascular disease; CVA, cerebrovascular accident; MDRD, Modification of Diet in Renal Disease; cTn, cardiac troponin; CKD-EPI, Chronic Kidney Disease Epidemiology Collaboration; RV/LV ratio, right ventricular to left ventricular diameter ratio; PAOI, pulmonary vascular obstruction index; SVC, superior vena cava.
a P < .05.
bcTn (n = 88).
cRV–LV ratio and PAOI (n = 81), SVC diameter (n = 79).
Renal Function Markers and Adverse Outcome
On admission, 3 patients were classified as having severe RD (<30 mL/min), 1 using eGFRMDRD and 2 using eGFRCKD-EPI. eGFRCKD-EPI < 60 mL/min was noted in 44 (44%) patients, while eGFRMDRD < 60 mL/min was noted in 38 (38%) patients. BUN > 20 mg/dL was noted in 47 (46%) patients. Areas under the adverse outcome ROC curves were compared for BUN, eGFRMDRD, and eGFRCKD-EPI (Figure 1). Statistically significant differences was present between eGFRCKD-EPI and eGFRMDRD (P = .008, 95% confidence interval [CI]: 0.010-0.067) in prediction of adverse outcome but not between BUN and eGFRCKD-EPI (P = .306, 95% CI: −0.055-0.174) or between BUN and eGFRMDRD (P = .116, 95% CI: −0.024-0.220).

Receiver–operating characteristic curves of indicators of renal dysfunction.
Multivariate Analysis
Variables significantly associated with adverse outcome in the univariate analyses were included in a multivariate logistic regression model (Table 3). Using multivariate analysis, we found that eGFRCKD-EPI (odds ratio [OR] 7.18, 95% CI 1.12-45.92, P = .037), cTn positivity (OR 9.84, 95% CI 1.24-77.93, P = .03), and heart rate (>110 beat/min; OR 8.47, 95% CI 1.44-49.77, P = .018) were independent predictors of adverse outcome. In the next analysis, eGFRMDRD was used instead of eGFRCKD-EPI, and eGFRMDRD (< 60 mL/min/1.73 m2) (OR 6.42, 95% CI 1.08-38.26, P = .041), cTn positivity (OR 9.24, 95% CI 1.18-72.18, P = .034), and heart rate (>110 beats/min; OR 7.20, 95% CI: 1.29-40.21, P = .025) remained significant.
Univariate and Multivariate Analyses of Factors Associated With Adverse Outcome.
Abbreviations: SBP, systolic blood pressure; CVD, cardiovascular disease; CVA, cerebrovascular accident; BUN, blood urea nitrogen; RV/LV ratio, right ventricular to left ventricular diameter ratio; eGFR, estimated glomerular filtration rate; MDRD, Modification of Diet in Renal Disease; cTn, cardiac troponin; CKD-EPI, Chronic Kidney Disease Epidemiology Collaboration; PAOI, pulmonary vascular obstruction index.
a P < .05.
beGFR categorical (<60 mL/min/1.73 m2).
Relationship Among CT Findings, Systolic Blood Pressure, and Renal Function Markers
eGFRCKD-EPI correlated significantly with RV/LV ratio, PAOI, and SVC diameter (r = −.31, P < .01; r = −.23, P < .05; and r = −.25, P < .05, respectively). eGFRMDRD was also correlated significantly with the above-mentioned variables (r = −.32, P < .01; r = −.24, P < .05; and r = −.31, P < .01, respectively). Blood urea nitrogen correlated significantly with SVC diameter only (r = .23, P < .05). eGFRCKD-EPI, eGFRMDRD, nor BUN was significantly correlated with SBP (P = .72, P = .91, and P = .12, respectively). In addition, SBP and RV–LV ratio were significantly correlated (r = −.22, P < .05).
Stratified Evaluation and Combined Approach
When 45 mL/min was used as the cutoff value for eGFRCKD-EPI, significant differences between patients above and below the cutoff were found for APE-related death (33% vs 6%, P = .001), all-cause mortality (40% vs 10%, P = .002), and thrombolytic therapy (27% vs 5%, P = .004). Similarly, BUN > 32 mg/dL was closely associated with APE-related death (30% vs 5%, P = .001), all-cause mortality (40% vs 8%, P = .000), and the need for thrombolytic therapy (30% vs 5%, P = .001).
Also, we compared combined cTn-positivity and eGFRCKD-EPI (<45 mL/min and <60 mL/min, respectively) and cTn-positivity alone for 30-day mortality in patients with APE. Combined approach improved positive predictive value (PPV) without significant change in negative predictive value (NPV) for 30-day mortality (Table 3). In addition, cTn-positive patients with eGFRCKD-EPI (<45 mL/min and <60 mL/min, respectively) showed 50%, 36% 30-day mortality, although cTn-positive patients with eGFRCKD-EPI (≥45 mL/min and ≥60 mL/min, respectively)) 14%, 11% mortality (Table 4).
Comparison of cTn Positivity Alone Versus cTn Positivity in Combination With eGFRCKD-EPI for 30-Day Mortality.
Abbreviations: cTn, cardiac troponin (n = 88), death (n = 11); PPV, positive predictive value; NPV, negative predictive value; eGFR, estimated glomerular filtration rate; CKD-EPI, Chronic Kidney Disease Epidemiology Collaboration.
Discussion
This cross-sectional study demonstrated that both eGFR and BUN are closely associated with APE prognosis in normotensive patients. The same median eGFR values (64 mL/min) were found using both equations. Prevalence of intermediate RD (<60 mL/min/1.73 m2) was found to be 38% and 44% with eGFRMDRD and eGFRCKD-EPI, respectively. These findings are consistent with previous reports. 6 –8,13 Both estimates of eGFR, together with troponin positivity and heart rate (>110 beats/min), were independent predictors for adverse outcome. However, ROC curve analysis showed that eGFRCKD-EPI was a better predictor of adverse outcome than was eGFRMDRD. Blood urea nitrogen was a significant predictor for adverse outcome in the univariate analysis but not in the multivariate analysis.
Renal Dysfunction
Renal arterial blood flow (RABF) is considered to be the most important determinant of GFR. However, increasing evidence suggests that increased CVP is what mainly determines renal function. 14 –16 In addition, a number of neurohumoral mediators take part in the regulation of GFR. 17 Decreased GFR in APE can be explained by several mechanisms. First, impaired RABF secondary to decreased cardiac output (CO) can lead to a decrease in GFR. Decreased CO is due to a decrease in preload that results from a number of mechanisms, including increased pulmonary vascular resistance, septal shift, interventricular asynchrony, and/or pericardium-mediated right ventricle–left ventricle interaction. 18 Reduction in RBAF is often considered to be synonymous with hypotension; however, impaired GFR in normotensive patients with HF is a well-known condition and is explained by the phenomenon of arterial underfilling. 18,19 Similarly, normotensive patients having APE with a low cardiac index have also been identified. 20 The presence of normotension (normal arterial blood pressure), despite decreased CO appears to result from NHA including sympathetic overactivity. As with HF, APE has been shown to be associated with elevated neuroendocrine markers. 21 –23 In addition to estimating renal function, BUN is accepted as a biomarker of neurohormonal activation and is reported to be a better predictor of survival than is eGFR in patients with HF. 24 In the present study, increased BUN and heart rate in patients with adverse outcomes indicates an existing NHA. However, BUN was not a better predictor of adverse outcome compared with eGFR, in the univariate or multivariate analyses. Second, accumulating evidence suggests that poor forward flow alone does not explain the decrease in GFR. Renal venous congestion seems to be secondary to increased CVP and could have a significant impact on GFR. This view is supported by previous studies performed in patients with HF. 14 –16 Similar observations were also made in studies investigating the relationship between APE and GFR. 6 –8 These studies suggest a significant correlation between GFR and various radiological and echocardiographic right ventricular functional parameters, since eGFRMDRD is significantly correlated with the RV–LV ratio, PAOI, and SVC diameter. These findings support the hypothesis that this may be a predominant mechanism.
Estimated Glomerular Filtration Rate Calculations
There are a number of methods used to estimate GFR in daily practice; however, the simplified MDRD equation is the most popular. This method was derived from individuals with renal disease and has been criticized due to its underestimation of GFR in healthy individuals. 2 The eGFRCKD-EPI equation was developed in a study involving healthy individuals, and it includes the same 4 variables as MDRD years but uses a 2-slope “spline” model. 25 Currently, the CKD-EPI equation is considered to be the most accurate method for estimating GFR. 2 Interestingly, in the present study, the eGFRCKD-EPI equation classified more patients as having intermediate RD. Similar findings have already been reported in patients with multiple myeloma, 26 but the advanced age of the patients in this study (>70) may be the reason for this discrepancy. 27
Renal Disease and Mortality
Chronic kidney disease and chronic renal failure are associated with increased APE risk. 28 However, there is conflicting evidence about the relationship between APE mortality and CKD, especially in patients undergoing dialysis. Autopsy studies suggest that APE-related mortality rates in patients with end-stage renal disease (ESRD) are lower than those in the general population. 29,30 A variety of reasons have been suggested to explain these observations, including platelet dysfunction, bleeding tendency, and anticoagulation therapy during dialysis. In a recent study, Fabbian et al found no relationship between APE-related mortality and CKD or ESRD, 31 yet other studies contradict these findings. 32 Several studies have shown an association between low eGFR levels and CKD; however, it is not clear whether the observed low eGFR levels were due to preexisting renal impairment or kidney injury. To clarify these findings, 2 studies by Kost et al were performed. In the first study, neutrophil gelatinase-associated lipocalin (N-GAL), an accepted marker of acute kidney injury, and cystatin C levels were evaluated in patients with APE. 7 Serum N-GAL levels were significantly higher in patients who died from APE; however, cystatin C levels were a significant predictor of all-cause mortality both in univariate and in multivariate analyses. In the second study, eGFR was measured at admission and 72 hours later, and it was found that eGFR levels improved in patients with a good prognosis but remained low in patients with a poor prognosis. 8 As a result, low GFR levels in patients with APE seem to be an outcome of multiple pathophysiological processes including CKD and hemodynamic dysfunction.
According to the current guidelines, risk stratification of APE is based on cardiac biomarkers. 1 Although the NPV of these biomarkers for APE mortality is high, the PPV is quite low. 1,33 In this study, combination of cTn and eGFRCKD-EPI has been shown to improve the risk stratification based on cTn only by increasing PPV.
Limitations
This study was retrospective, and the patient population was relatively small. Our laboratory switched from measuring troponin I to troponin T during the study; however, previous reports have shown that both troponins have a similar ability to predict the prognoses of patients with APE. In this study, we found a relatively higher adverse event rate (22%), although similar rates have been reported in the literature. 34 –37 Our study population was comprised of a relatively more advanced age and high comorbidity group including CLD and diabetes mellitus, which might have contributed to higher rate of adverse outcomes. Variation in creatinine assays may be a source of bias in eGFR measurement. Isotope dilution mass spectrometry is a standard method to measure reference levels of sCr; however, our laboratory does not have this capability. In this study, the prevalence of renal disease was lower relative to that in other studies, which may be result of recall bias. It may also indicate the presence of subclinical kidney disease, since patients with adverse outcomes had a higher incidence of comorbidities (HF, DM, and CVD), although this did not reach significance. However, in daily practice, it is difficult to determine whether low eGFR values in patients with APE are dependent on hemodynamic changes or CKD. Regardless of the reason, decreased eGFR seems to be closely associated with APE prognosis.
Conclusion
There is a close relationship between RD and the prognosis of APE. eGFRCKD-EPI, but not BUN, predicts adverse outcomes better than does eGFRMDRD. eGFRCKD-EPI seems to be a potential marker for risk stratification in normotensive patients with APE.
Footnotes
Authors’ Note
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
