Abstract
We tested the ability of exercise testing to predict not only survival, but also time to clinical worsening (TTCW) in idiopathic versus associated pulmonary arterial hypertension (PAH).
136 patients with PAH (85 idiopathic and 51 with associated conditions) underwent cardiopulmonary exercise testing and a 6-min walk test. Death or transplantation, and clinical worsening events were recorded.
32 patients died and four had lung transplantation. In a univariate analysis, PAH patients survival was associated with oxygen uptake (V′O2) at peak exercise and at the anaerobic threshold, ventilatory equivalent for carbon dioxide (minute ventilation (V′E)/carbon dioxide production (V′CO2) at the anaerobic threshold (at)), V′E/V′CO2 slope and distance walked. TTCW was associated with peak V′O2 and V′O2,at, V′E/V′CO2,at, end-tidal carbon dioxide tension measured at the anaerobic threshold, peak oxygen pulse, increase in oxygen pulse and distance walked. In a multivariable analysis, distance walked and V′E/V′CO2,at predicted survival, and only peak V′O2 predicted TTCW. The receiver operating characteristic curve-derived cut-off values were 305 m for the 6-min walk distance, 54 for V′E/V′CO2,at and 11.6 mL·kg−1·min for peak V′O2. In the subgroup with associated PAH, no variable independently predicted either survival or clinical worsening.
We conclude that several exercise variables predict survival and clinical stability in idiopathic PAH. Exercise variables are less accurate predictors of outcome in associated PAH.
- Cardiopulmonary exercise testing
- clinical worsening
- pulmonary arterial hypertension
- pulmonary hypertension
- 6-min walk test
- survival
The symptomatology of pulmonary arterial hypertension (PAH) is dominated by dyspnoea and fatigue induced by exercise until the final stages of the disease, when the patients are symptomatic at rest [1, 2]. Accordingly, a variety of variables measured during a cardiopulmonary exercise test (CPET) and, more simply the 6-min walk distance (6MWD) are used in clinical practice to estimate disease severity [3, 4].
The 6-min walk test (6MWT) has been shown to be an independent prognostic marker [5–10], and as such has served as a primary end-point in most randomised controlled trials of new therapies for PAH [11]. Peak oxygen uptake (V′O2) and ventilatory responses during CPET also relate to survival in pulmonary hypertension. This has been shown in idiopathic PAH (IPAH) [12] and in a cohort of pulmonary hypertension patients composed of IPAH, PAH with associated conditions (APAH) and chronic thromboembolic pulmonary hypertension (CTEPH) patients [7].
Because PAH is still an incurable disease with limited survival, clinical stability is a desirable therapeutic goal, especially in patients who are not too severely ill, or still in New York Heart Association (NYHA) functional class II or early NYHA functional class III [1, 2]. Time to clinical worsening (TTCW) has emerged as an improved primary end-point in newly designed event-driven randomised controlled trials in PAH [11]. However, no study to date has addressed the question whether exercise capacity predicts the TTCW. Another incompletely answered question is whether the predictive value of exercise testing is equivalent in IPAH or APAH.
The purpose of the present study is to evaluate the prognostic value of CPET variables and 6MWT in IPAH versus APAH, and to determine their ability not only to predict survival, but also TTCW.
METHODS
Protocol
Our study included the CPET and 6MWT data of 136 PAH patients from the pulmonary hypertension clinic of Erasmus University Hospital, Brussels, Belgium, between November 2001 and July 2010. The study was approved by the Erasmus Hospital Institutional Review Board. The mortality end-point was defined as all-cause mortality or lung transplantation, with the remaining cases designated as event-free survival. The clinical worsening (CW) end-point was defined as previously reported [11]: all causes of mortality; nonelective hospital stay for PAH (for initiation of prostanoids, lung transplantation or atrial septostomy); disease progression defined as a reduction from baseline in the 6MWD by 15% confirmed by two studies performed within 2 weeks plus worsening functional class.
Patients
The diagnosis of PAH rested on a right heart catheterisation with demonstration of an increase in mean pulmonary artery pressure (>25 mmHg), a normal pulmonary artery wedge pressure (<15 mmHg), no identifiable cardiac or pulmonary cause, and possibly associated with conditions such as appetite suppressant intake, connective tissue disease (CTD), liver cirrhosis, HIV infection and congenital left-to-right shunt (congenital heart disease; CHD) [1, 2]. 70 patients had IPAH and 66 patients had APAH that were previous to intake of anorexigen in 15 patients, CTD (all systemic sclerosis) without lung function impairment in 19, hepatic cirrhosis in 11, HIV infection in four, CHD in 16 patients or schistosomiasis in one patient. Patients with PAH associated with the intake of anorexigens were considered as IPAH patients, as recent studies have shown that anorexigens only trigger the disease, which is otherwise indiscernible [2, 13]. Thus, the study considered two subgroups of 85 patients with IPAH and 51 patients with APAH.
Exercise testing
CPET
Each patient underwent standard cycle ergometer incremental CPET until the symptom-limited maximum [14]. The CPET protocol consisted of pedalling at 0–20 W during the first 3 min and then an incremental increase in load from 5 to 15 W·min−1, aiming to obtain an exercise duration of between 8 and 12 min. Because of equipment renewal over the years, ventilation and gas analysis were performed using a CPX/D meter (Medical Graphics, St Paul, MN, USA) in 16 tests and a VMax meter (SensorMedics, Yorba Linda, CA, USA) in 120 tests. The gas analysers and pneumotachograph were calibrated prior to each test. Cardiac frequency (fc)and blood pressure were obtained via automatic standard ECG and sphygmomanometer.
Peak V′O2, V′O2 at the anaerobic threshold (V′O2,at), the ventilatory equivalent for carbon dioxide (minute ventilation (V′E)/carbon dioxide production (V′CO2)) measured at anaerobic threshold (V′E/V′CO2,at) or as a slope from 1 min after the beginning of loaded exercise to the end of the isocapnic buffering period (V′E/V′CO2 slope), the end-tidal carbon dioxide tension (PET,CO2) measured at the anaerobic threshold, the maximum–rest change in PET,CO2 (ΔPET,CO2), the oxygen pulse, calculated as the V′O2/fC ratio, at peak exercise (O2pulse), the difference in O2 pulse between rest and peak exercise (ΔO2pulse), the peak systolic blood pressures and the occurrence of a right-to-left exercise-induced shunt (EIS) through a patent foramen ovale following previously described criteria [15] were reviewed as potential prognostic markers. Peak V′O2 was defined as the highest V′O2 measured during a period of 20 s at the end of the CPET and anaerobic threshold was determined using the V-slope method [16]. In the case of uncertainty, the anaerobic threshold was counter checked using the nadir of ventilatory equivalents [14].
Statistics
Statistical analysis SPSS version 18.0.0 (SPSS Inc., Chicago, IL, USA) was performed, including the 136 PAH patients. The 85 IPAH and 51 APAH patients were also considered separately. PAH patients associated with CHD were removed from the analysis of the influence of EIS on survival and TTCW. Time of origin was the date of exercise to date of death or transplantation, the patient was censored at the end of the study if still alive.
Data are presented in mean±sd. PAH, and IPAH and APAH subgroups were compared by unpaired t-tests. Proportion differences were tested by either Chi-squared or Fisher’s exact tests depending of the number of patients in each group.
A Cox proportional hazards regression analysis was used to detect predictors associated with survival and with TTCW. Hazard ratio, 95% confidence intervals and p-values from the likelihood ratio test are given.
For the variables that were predictive of survival or TTCW, receiver operating characteristic (ROC) curves were plotted at 4 yrs for death and at 2 yrs for TTCW. The area under the curve (AUC) with 95% confidence interval and p-value was determined using the nonparametric method. When the lowest 95% confidence interval was >50% and the p-value <0.05, the optimal cut-off point for predicting survival was identified on the basis of the highest sum of sensitivity and specificity, and Kaplan–Meier cumulative survival plots constructed for pattern above and below the threshold were used to describe survival rates. The log rank test was used to compare survival curves (figs 1–4).
Kaplan–Meier cumulative survival curves for the 6-min walk distance (6MWD) in a) 136 patients with pulmonary arterial hypertension (PAH) and b) 85 patients with idiopathic PAH. Cut-off value determined by receiver-operating characteristic curve.
Kaplan–Meier cumulative survival curves for the ventilatory equivalent for carbon dioxide at the anaerobic threshold (V′E/V′CO2,at) in a) 136 patients with pulmonary arterial hypertension (PAH) and b) 85 patients with idiopathic PAH. Cut-off value determined by receiver-operating characteristic curve.
Kaplan–Meier cumulative survival curves for no, one or two additional risks factors being a 6-min walk distance <305 m in pulmonary arterial hypertension (PAH) and 307 m in idiopathic PAH (IPAH) and a ventilatory equivalent for carbon dioxide at the anaerobic threshold >54 in PAH and IPAH in a) 136 patients with PAH and b) 85 patients with IPAH.
Kaplan–Meier cumulative curves for clinical worsening for peak oxygen uptake (peak V′O2) in a) 136 patients with pulmonary arterial hypertension (PAH) and b) 85 patients with idiopathic PAH. Cut-off value determined by receiver-operating characteristic curve.
Multivariable Cox regression analysis with a forward selection procedure was used to determine independent predictors from the variables with p<0.10 in univariate analysis. In all analyses, a p-value <0.05 was considered significant.
RESULTS
Demographic, haemodynamic and clinical data, and pertinent exercise variables are given in tables 1 and 2. The diagnostic right heart catheterisations had been performed at the time of the exercise test in 82 patients, and otherwise 6±11 months before. At the time of the exercise test, 37 patients received targeted treatment for PAH, such as endothelin receptor antagonists, phosphodiesterase-5 inhibitors or prostacyclins.
Survival
Of the 136 PAH patients followed (44.2±28.3 months), 32 died (at between 1.1 and 77.2 months) and four had lung transplantation (at between 14.1 and 72.6 months). PAH patients had a Kaplan–Meier survival rate at 1, 2, 3 and 4 yrs of 94, 86.5, 75.1 and 72.7%, respectively.
In the IPAH patients (followed over 45.3±30.4 months), 18 died (at between 1.1 and 41.7 months) and four benefited from a lung transplantation (at between 14.1 and 72.6 months); survival at 1, 2, 3 and 4 yrs was 94.0, 88.7, 73.4 and 69.4%, respectively.
In the APAH patients (followed over 42.3±24.6 months), 14 died (at between 2.5 and 77.2 months) and survival at 1, 2, 3 and 4 yrs was 94.1, 82.9, 77.9 and 77.9%, respectively.
As can be seen in table 1, nonsurvivors compared with survivors had similar sex distribution and were of the same age, functional state and haemodynamic severity of pulmonary hypertension (with exception of a lower mixed venous oxygen saturation (Sv,O2) in PAH and APAH and a higher right atrial pressure (Pra) in the APAH subgroups). Exercise capacity variables were more altered in the nonsurvivors PAH patients and this appeared to be driven entirely by the IPAH subgroup.
Univariate Cox analyses for prediction of survival are shown in table 3.
Peak V′O2, V′O2,at, V′E/V′CO2,at, V′E/V′CO2 slope and 6MWD were predictive of survival in PAH and IPAH patients. No variable was associated with mortality in APAH patients. The 6MWD and V′E/V′CO2,at predicted death independently in PAH and in IPAH patients. Correcting the 6MWD for bodyweight did not affect its predictive value.
Results for receiver operating characteristics analysis are shown in table 4. In the PAH patients, optimal cut-off values to predict survival were: 11.5 mL·kg−1·min−1 for peak V′O2; 8.8 mL·kg−1·min−1 for V′O2,at; 54 for V′E/V′CO2,at (fig. 2a); 62 for V′E/V′CO2 slope; and 305 m for 6MWD (fig. 1a).
In the IPAH patients, optimal cut-off values to predict survival were: 10.6 mL·kg−1·min−1 for peak V′O2; 9.8 mL·kg−1·min for V′O2,at; 54 for V′E/V′CO2,at (fig. 2b); and 307 m for 6MWD (fig. 1b). No optimal cut-off point was determined for V′E/V′CO2 slope as the ROC curve p-value was >0.05 and the lowest 95% confidence interval of the area under the curve was <0.5. Kaplan–Meier cumulative survival curves for 0, 1 or 2 additional independant risks factors are shown in fig. 3.
Clinical worsening
Of the 136 PAH patients followed, 88 encountered a clinical worsening (at between 1.0 and 101 months). The rate of CW after 1, 2, 3 and 4 yrs for the PAH group was 33.6, 51.6, 59.8 and 64.6%, respectively. Of the IPAH patients, 56 had a CW (at between 1.0 and 101 months) and rate of CW at 1, 2, 3 and 4 yrs was 41.9, 58.2, 67.3 and 73.1%, respectively. Of the APAH patients, 24 had a CW (at between 1.0 and 79.3 months), and the rate of CW at 1, 2, 3 and 4 yrs was 19.8, 38, 47.5 and 51.6%, respectively.
As shown in table 2, the PAH patients who presented with a CW were as frequently male, and of similar age and functional state. Patients with subsequent CW had higher pulmonary vascular resistance. PAH had higher Pra, IPAH had higher mean pulmonary artery pressure, and PAH and IPAH had lower cardiac output, cardiac index and Sv,O2. Exercise capacity variables were more altered in the patients with CW, and this appeared to be driven by the IPAH subgroup. Univariate Cox analysis for prediction of CW is shown in table 5.
Peak V′O2, V′O2,at, V′E/V′CO2 slope, V′E/V′CO2,at, PET,CO2, O2pulse, ΔO2pulse and 6MWD predicted TTCW in both PAH and IPAH patients. Correcting the 6MWD for bodyweight did not affect its predictive value. No variables were associated with TTCW in the APAH patients. Peak V′O2 predicted TTCW independently in both PAH and IPAH patients, as shown in by the ROC analysis.
Results for receiver operating characteristics analysis are shown in table 6. Of the PAH patients, optimal cut-off values to predict TTCW were: 11.6 mL·kg−1·min−1 for peak V′O2 (fig. 4b); 9 mL·kg−1·min−1 for V′O2,at; 46 for V′E/V′CO2,at; 55 for V′E/V′CO2 slope; 23.5 for PET,CO2,at; 5.3 mL·beat−1 for O2pulse; 2.6 mL·beat−1 for ΔO2pulse; and 367 m for 6MWD.
Of the IPAH patients, optimal cut-off values to predict CW were: 11.8 mL·kg−1·min−1 for peak V′O2 (fig. 4b); 9.6 mL·kg−1·min−1 for V′O2,at; 51 for V′E/V′CO2,at; 59 for V′E/V′CO2 slope; 23.5 for PET,CO2,at; 2.1 mL·beat−1 for ΔO2pulse; and 367 m for 6MWD. No optimal cut-off point was determined for O2pulse IPAH as the ROC curve p-value was >0.05 and the lowest 95% confidence interval of the AUC was <0.5.
DISCUSSION
The present results show that exercise capacity predicts not only survival, but also clinical stability in PAH, with a better discrimination of exercise testing variables in IPAH than in APAH. In this study, independent predictors of survival were the 6MWD and V′E/V′CO2,at for IPAH, while no exercise variable independently predicted survival in APAH. As for TTCW, this was predicted only by peak V′O2 and only in IPAH.
PAH is a right heart failure syndrome [19]. Exercise capacity in heart failure is largely determined by maximal cardiac output, which, in pulmonary hypertension, is determined by right ventricular function. The afterloaded right ventricle relies on fC more than on stroke volume to increase flow output [19], which translates into a more decreased oxygen pulse during exercise for patient with a worse prognosis [7], as is also shown in the present results by univariate analysis. However, O2pulse is only an indirect measure of stroke volume derived from the Fick equation with the assumption of unchanged arterio-venous oxygen content difference. Maximal cardiac output is related to maximum V′O2 and also to maximal workload, which is in turn related to the maximal average running or walking speed [4]. According to this reasoning, the information content of peak V′O2 and the 6MWD in PAH are equivalent. Accordingly, in the present study, both peak V′O2 and the 6MWD predicted survival and TTCW in IPAH.
Previous studies in PAH have shown that the 6MWD was more sensitive than peak V′O2 to targeted therapies, such as beraprost [20] or sitaxsentan [21]. Part of this greater sensitivity of the 6MWD was ascribed to relatively less expertise in the practice of CPET in the centres that participated to these studies [22]. In the present study, the predictive capability of CPET variables in IPAH was found to be similar to those previously reported, making insufficient expertise unlikely. From another point of view, while the 6MWD and not peak V′O2 independently predicted survival, peak V′O2, and not the 6MWD, predicted CW. The reasons for this apparent contradiction are unclear, but this result underlines the interest of CPET when added to a 6MWT in the evaluation of PAH patients.
The correlation between a maximal average running speed and maximum V′O2 or peak V′O2 is generally significant, but rather loose, with correlation coefficients in the range of 0.5–0.7 in normal subjects [23], as well as in patients with IPAH [5, 22, 24]. This is explained by variable mechanical efficiency of running or walking, which is related in part to different body dimensions or weight. Correction of the 6MWD for bodyweight has indeed been shown to improve these correlations in chronic obstructive pulmonary disease [18], as well as in PAH [22, 24]. However, in the present study, correcting the 6MWD for bodyweight did not improve its prognostic value.
In the present study, the nonsurvivors had similar pulmonary haemodynamics compared with survivors, even though Sv,O2 was lower and Pra higher (in the APAH group only) in the nonsurvivors. Haemodynamic severity of PAH was more clearly associated with the occurrence of clinical deterioration. However, while these results are in keeping with the previously known poorer prognosis associated with more severe pulmonary hypertension and a lower cardiac output [1, 2, 5, 6], further analysis of the compared prognostic values of haemodynamic and exercise test variables was limited by the fact that the measurements were separated in time in too many of the patients.
The present results confirm the prognostic value of the ventilatory equivalents for carbon dioxide previously reported in chronic heart failure [25, 26], as well as in PAH [7, 12], but with significance only in the patients with IPAH and not those with APAH. This is probably related to the inhomogeneity of the APAH group, as it is known that survival is much better in CHD-APAH and worse survival in systemic sclerosis APAH [27]. Whether the V′E/V′CO2 slopes during CPET are different in APAH subcategories is not exactly known.
It has been recommended to measure V′E/V′CO2 below or at the anaerobic threshold, as the V′E/V′CO2 slope calculated on the entire CPET may be influenced by the maximal or peak level of exercise and variable associated acidosis [14]. In the present study, the predictive value of V′E/V′CO2,at was better than that of the V′E/V′CO2 slope, even though the latter was measured after the beginning of loaded exercise to the end of the isocapnic buffering period.
Another problem of the V′E/V′CO2 ratio or slope may be the sudden increase occurring with right-to-left shunting [15]. For that reason, Wensel et al. [12] excluded patients with a resting patent foramen ovale from the analysis of V′E/V′CO2 slope in relation to survival. The development or the persistence of such shunting during the course of therapy has been found to be associated with an altered prognosis in PAH [28]. In the present study, however, EIS did not predict survival or TTCW.
Groepenhoff et al. [7] found that V′E/V′CO2,at or V′E/V′CO2 slope predicted survival in their mixed PAH–CTEPH population with cut-off values in the same range as ours (i.e. 52 for V′E/V′CO2,at and 48 for the V′E/V′CO2 slope). It is of interest that survival cut-off values for V′E/V′CO2 are higher in PAH than in heart failure, with typical values >50 and of ∼35, respectively [25, 26, 29]. Both conditions are associated with increased ventilation at any given level of metabolic rate [30] in relation to increased chemosensitivity and physiological dead space [31, 32]; however, the respective contributions of both mechanisms are not exactly known.
Peak V′O2 values inferior to 11.5 and 10.6 mL·kg−1·min−1 in our PAH and IPAH subjects, respectively, were associated with a decreased survival. The cut-off of 10.6 mL·kg−1·min−1 is in agreement with a report by Wensel et al. [12] on similar severity IPAH patients, while Groepenhoff et al. [7] found a higher cut-off of 13.2 mL·kg−1·min−1 in a mixed group of PAH and CTEPH patients. The isolated impact of CTEPH on the predictive value of peak V′O2 remains to be investigated. Slight differences in “cut-off values” and predictions are attributable to differences in severity of disease in apparently similar patient populations.
A limitation to the present study is the heterogeneous nature of the APAH group, with too small numbers to identify specific profiles of subgroups such as CTD–PAH or CHD–PAH. Another limitation is that the walking ability may be impaired in CTD patients, so that the 6MWT as a measure of exercise capacity in these patients remains insufficiently validated [33]. Finally, because of the heterogeneous nature and smaller size of the APAH group, it happened that all conclusions for IPAH patients were transposable to the PAH patient population as a whole in the present study, but this of course may be affected the nature of associated conditions and sizes of APAH subgroups.
To our knowledge, this is the first report of prediction of TTCW from exercise testing variables. Our results clearly show that clinical stability is better predicted in IPAH than in APAH and that, for this purpose, peak V′O2 may be superior to the 6MWT. This is relevant to the definition of clinical deterioration and inclusion criteria for event-driven clinical trials of new therapeutic interventions in PAH. Moreover, our study documents the importance of V′E/V′CO2 as an outcome predictor independently associated with survival and also able to predict TTCW in PAH and IPAH.
Footnotes
Support Statement
G. Deboeck is recipient of a European Respiratory Society Fellowship (LTRF 21-2010).
Statement of Interest
None declared.
- Received December 12, 2011.
- Accepted February 16, 2012.
- ©ERS 2012