Next Article in Journal
Identification of Glycoxidative Lesion in Isolated Low-Density Lipoproteins from Diabetes Mellitus Subjects
Next Article in Special Issue
Enhanced Microvascular Adaptation to Acute Physical Stress and Reduced Oxidative Stress in Male Athletes Who Consumed Chicken Eggs Enriched with n-3 Polyunsaturated Fatty Acids and Antioxidants—Randomized Clinical Trial
Previous Article in Journal
Unraveling the Power of Topical Inhaled Steroids in Treating Laryngeal Granulomas: A Systematic Review
Previous Article in Special Issue
Recovery after Running an “Everesting” Mountain Ultramarathon
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cardiopulmonary Exercise Testing in Patients with Heart Failure: Impact of Gender in Predictive Value for Heart Transplantation Listing

by
Pedro Garcia Brás
1,*,
António Valentim Gonçalves
1,
João Ferreira Reis
1,
Rita Ilhão Moreira
1,
Tiago Pereira-da-Silva
1,
Pedro Rio
1,
Ana Teresa Timóteo
1,2,
Sofia Silva
1,
Rui M. Soares
1 and
Rui Cruz Ferreira
1
1
Cardiology Department, Santa Marta Hospital, Central Lisbon Hospital University Center, 1169-024 Lisbon, Portugal
2
NOVA Medical School, Faculdade de Ciências Médicas (NMS|FCM), 1169-056 Lisbon, Portugal
*
Author to whom correspondence should be addressed.
Life 2023, 13(10), 1985; https://doi.org/10.3390/life13101985
Submission received: 30 August 2023 / Revised: 22 September 2023 / Accepted: 27 September 2023 / Published: 29 September 2023
(This article belongs to the Special Issue New Insights into Cardiovascular and Exercise Physiology)

Abstract

:
Background: Exercise testing is key in the risk stratification of patients with heart failure (HF). There are scarce data on its prognostic power in women. Our aim was to assess the predictive value of the heart transplantation (HTx) thresholds in HF in women and in men. Methods: Prospective evaluation of HF patients who underwent cardiopulmonary exercise testing (CPET) from 2009 to 2018 for the composite endpoint of cardiovascular mortality and urgent HTx. Results: A total of 458 patients underwent CPET, with a composite endpoint frequency of 10.5% in females vs. 16.0% in males in 36-month follow-up. Peak VO2 (pVO2), VE/VCO2 slope and percent of predicted pVO2 were independent discriminators of the composite endpoint, particularly in women. The International Society for Heart Lung Transplantation recommended values of pVO2 ≤ 12 mL/kg/min or ≤14 if the patient is intolerant to β-blockers, VE/VCO2 slope > 35, and percent of predicted pVO2 ≤ 50% showed a higher diagnostic effectiveness in women. Specific pVO2, VE/VCO2 slope and percent of predicted pVO2 cut-offs in each sex group presented a higher prognostic power than the recommended thresholds. Conclusion: Individualized sex-specific thresholds may improve patient selection for HTx. More evidence is needed to address sex differences in HF risk stratification.

1. Introduction

Cardiopulmonary exercise testing (CPET) is a critical complementary test in the evaluation of patients with heart failure (HF) with reduced ejection fraction (HFrEF), particularly in selectin patients who may benefit from heart transplantation (HTx) [1,2]. Peak O2 consumption (pVO2) [3,4,5] and the VE/VCO2 slope (minute ventilation–CO2 production relationship) [3,5,6] are reliable indicators of heart failure events. A cut-off for pVO2 of ≤12 mL/kg/min is recommended to guide HTx listing for patients receiving β-blocker therapy, and a cut-off of 14 mL/kg/min may be used for patients intolerant to β-blockers, according to the 2016 International Society for Heart Lung Transplantation (ISHLT) listing criteria for HTx [7,8]. In female patients, alternative parameters such as a VE/VCO2 slope of >35 and a percent of predicted pVO2 ≤ 50% may be considered to guide HTx listing [7]. However, the data supporting these values come from studies that enrolled mostly male patients, with a sample that was between 80 and 90 percent male [2,5].
Indeed, female patients are underrepresented in HFrEF trials, although they account for around half of the adult HFrEF population [9]. Notably, in studies exploring CPET parameters in HFrEF, this gap in female representation is even larger [1,10,11,12]. Thus, the current evidence on female HFrEF pathophysiology and exercise testing prognostic power is scarce, and therapy changes, risk stratification, and recommendations for advanced HF therapies may be impacted by our insufficient comprehension of potential sex variations in HF [13,14]. Several trials evaluating pVO2 have reported lower values in female patients compared to male patients, which can be explained by anatomical and physiological differences [9,15]. Female patients exhibit lower left ventricular dimensions, with lower stroke volume and lower diastolic compliance [16,17]; women show a higher prevalence of iron deficiency, have lower hemoglobin levels [18], and have inferior lean mass compared to male patients [9,19].
The HF-ACTION trial [10] assessed the prognostic power of CPET variables to predict all-cause mortality in HFrEF and reported that the prognosis associated with a given pVO2 differed by sex. Female patients generally present a more favorable outcome, and have a lower pVO2 and a higher percent of predicted pVO2 [1,10]. Taking into account the sex-based variations in the pathophysiology and development of HFrEF, several authors proposed that prognostic values for pVO2 and VE/VCO2 slope should be tailored for different patient populations [12,14,20]. According to the ISHLT [7], different CPET variables can be used for risk stratification in women. Nevertheless, there is insufficient evidence to support these at this time, as unbiased data are not available [1].
This study’s objective was to assess the predictive power of the traditional HTx CPET cut-off values in HF patients, comparing women and men.

2. Materials and Methods

2.1. Study Population

From 2009 to 2018, we performed a retrospective study of a prospective database in our center. We assessed consecutive HfrEF patients who underwent CPET, were in New York Heart Association (NYHA) classes II or III, and presented left ventricular (LV) dysfunction (LV ejection fraction ≤ 40%). Patients were referred to the Heart Failure team for evaluation to determine whether HTx or mechanical circulatory support (MCS) were indicated.

2.2. Study Protocol

The patient’s comorbidities, HF etiology, medication, NYHA class, HFSS (Heart Failure Survival Score) [21], laboratory tests, CPET data, and electrocardiographic and echocardiographic results were evaluated.

2.3. Patients Were Excluded If One of the following Was Present

Age under 18 years; submaximal CPET (peak RER of ≤1.05 [7]); previous HTx or elective HTx during follow-up; coronary revascularization in the last six months; concomitant conditions limiting maximal exercise, including previous stroke, peripheral arterial disease, or musculoskeletal conditions.

2.4. Cardiorespiratory Exercise Testing

The modified Bruce protocol was employed to assess maximal exercise tolerance on a GE Marquette Series 2000 treadmill, with equipment calibration before each exercise exam. The VE, VO2, and VCO2 values were acquired with a Vmax 229 (SensorMedics, Yorba Linda, CA, USA) gas analyzer. Continuous ECG monitoring was used to assess the heart rate (HRt). Blood pressure (BP) was obtained with a sphygmomanometer, and O2 saturation was tracked with pulse oximetry. An exercise test was considered maximal if the RER (respiratory exchange ratio) was above 1.05 [7].
The pVO2 was defined as the highest achieved 30 s average in maximal exercise, which was then normalised for body mass. The standard methods (combining V-slope preferentially and ventilatory equivalents) were used to determine the gas exchange threshold (GET). The VE/VCO2 slope was determined with the least squares linear regression. The minimum ventilatory equivalent for oxygen (minimum VE/VO2) was employed to calculate the COP (cardiorespiratory optimal point). The partial pressure of end-tidal carbon dioxide (PetCO2) was recorded both before exercise and at GET. Peak O2 pulse, measured in millilitres per beat, was computed by dividing the derived pVO2 by the highest HRt during exercise. The peak systolic BP was divided by the VE/VCO2 slope to determine the ventilatory power. The circulatory power was estimated by multiplying the peak systolic BP by pVO2. The heart rate reserve was calculated using the difference between the highest HRt attained during maximal effort and the resting heart rate. The difference between the maximal heart rate attained with exercise and the heart rate one minute in recovery was used to determine the HRt recovery.

2.5. Follow-Up and Endpoint

All patients with HF were under follow-up for a 36-month period. The composite endpoint was defined as the combination of cardiovascular mortality or urgent HTx. Data were collected from medical records from inpatient and outpatient visits.

2.6. Statistical Analysis

All analytical tests compared patients according to female or male sex. Statistical analyses were performed with the Statistical Package for the Social Sciences (SPSS) v23.0.
Regarding categorical variables, results were reported as absolute frequency (number) and relative frequency (%). Continuous variables were presented as mean and standard deviation if normal distribution, or as median and interquartile range (IQR) if non-normal distribution. Normality assumptions were tested using the Kolmogorov–Smirnov test and a visual histogram analysis.
The comparison of categorical variables was performed using the Pearson’s X2 test. The Mann–Whitney U test was used to compare variables with non-normal distribution and the Student’s t-test was used to compare variables with normal distribution.
The correlation between the CPET parameters and the composite endpoint was evaluated using a Cox hazards regression analysis. Variables presenting a p-value < 0.200 in the univariate analysis were included in a multivariate analysis, adjusted for potential confounders, in order to identify independent predictors of the composite endpoint and calculate adjusted hazard ratios (HR) in each sex subgroup. The HR and the 95% confidence interval (CI) were used to report the results.
A receiver operating characteristic (ROC) curve analysis was used to examine the sensitivity and specificity of each CPET parameter in predicting the composite endpoint, in accordance with the thresholds defined by the ISHLT [7]: pVO2 ≤ 12 mL/Kg/min (pVO2 ≤ 14 in patients not tolerant to β-blockers), VE/VCO2 slope > 35 and percent of predicted pVO2 ≤ 50%.
The threshold with the highest combination of specificity and sensitivity was estimated using the Youden index (J). The DeLong test [22] was employed to evaluate the difference in area under the curve (AUC) between groups. Additionally, the Kaplan–Meier analysis was used to assess the event-free survival rate. A log-rank test was performed to compare the sex subgroups based on the different pVO2, VE/VCO2 slope, and percent of predicted pVO2 thresholds indicated by the ISHLT [7] and based on the proposed cut-offs. A significance threshold of α = 5% was considered whenever a statistical hypothesis was being tested.

3. Results

3.1. Patient Characteristics

Our study included 458 patients who underwent maximal exercise testing (Figure 1). Of these patients, 79% were men, 57% had ischemic etiology, 76% were in NYHA II and 24% in NYHA III, with a mean LVEF of 29.7 ± 8.0%, and 24% had atrial fibrillation (AF). In addition, 79% were taking either an ACEi (angiotensin-converting enzyme inhibitor) or an ARB (angiotensin receptor blockers), with 17% on an angiotensin receptor/neprilysin inhibitor. Mineralocorticoid receptor antagonists (MRAs) were being taken by 73% and β-blockers by 86%. Additionally, sodium-glucose cotransporter-2 inhibitors (SGLT2i) were being taken by 10% of the patients; 64% of patients had an ICD and 22% had a cardiac resynchronization device (CRT-D). Moreover, there was no difference in the mean Heart Failure Survival Score (HFSS). Compared to male patients, female patients had a similar pVO2 and a higher percent of predicted pVO2. The mean respiratory exchange ratio (RER) was 1.14 ± 0.07. Table 1 lists the baseline characteristics of both groups as well as the CPET values.

3.2. Composite Endpoint

The composite endpoint occurred in 68 (14.8%) patients in 36 months of follow-up, with cardiovascular death occurring in 54 individuals and urgent HTx occurring in 14 patients (Table 2). No urgent MCS was required; 10.5% of female patients and 16.0% of male patients experienced the composite endpoint, with no significant difference between groups.

3.3. Prognostic Power of CPET Parameters

The pVO2 (HR 0.856), the VE/VCO2 slope (HR 1.064), and the percent of predicted pVO2 (HR 0.955) were associated with the composite endpoint in a multivariable Cox regression analysis, regardless of the sex subgroup. Table 3 displays the results of the uni- and multivariable models. The correlations in the multivariable model were independent of potential confounders such as body mass index, LVEF, age, sex, smoking, diabetes mellitus, or estimated glomerular filtration rate. In the multivariable analysis, most of the other exercise testing variables were not linked with the primary endpoint. The peak O2 pulse was associated with the endpoint in both female and male patients. The ventilatory power, the circulatory power, and the PetCO2 at GET were linked with the primary endpoint in male patients, as shown in Table 3.
In an ROC curve analysis, the pVO2, the VE/VCO2 slope, and the percent of predicted pVO2 were linked to the composite endpoint, both in females and males. The predictive ability of these variables was significantly higher in women compared to males, including for pVO2, VE/VCO2 slope, and the percent of predicted pVO2, as presented in Table 4. The ROC curves for these subgroups are illustrated in Figure 2 and Supplementary Figure S1. In addition, the predictive power of the peak O2 pulse was also significantly higher in female patients compared to males (AUC 0.816 vs. AUC 0.616, p = 0.023).
The circulatory power presented a slightly higher prognostic power than the recommended exercise testing parameters in men (AUC 0.713 vs. AUC 0.701, p = 0.161), albeit with no statistically significant differences in predictive power. Despite being significant predictors of the composite endpoint, additional CPET variables such as peak O2 pulse, ventilatory power, COP, PetCO2 at rest, and PetCO2 at GET had an inferior predictive power than the traditional CPET parameters (Table 4).

3.4. ISHLT Recommended Thresholds for HTx Listing

A pVO2 of ≤ 12 mL/kg/min (≤14 if the patient is intolerant to β-blockers) was present in 49 (11%) patients. This threshold was linked with poor HF outcomes (HR 3.487, p < 0.001). This pVO2 cut-off showed a sensitivity of 40% and a specificity of 94% in women, presenting a higher Youden index compared to men (J 0.34 vs. J 0.12), with a sensitivity of 21% and a specificity of 91%, as shown in Table 5. This cut-off was shown to be a strong discriminator of HF outcomes for both sex subgroups in a Kaplan–Meier analysis (Figure 3a).
A total of 166 (36%) patients showed a VE/VCO2 slope value higher than 35. The composite endpoint occurred at a higher rate in individuals over this threshold as well (HR 3.587, 95% CI 2.194–5.864, p < 0.001). This threshold revealed a substantially higher Youden index in women (J 0.65 vs. J 0.23), with sensitivity of 90% and a specificity of 75%, in comparison with male patients, with a sensitivity of 57% and a specificity of 66%. In the survival analysis, this VE/VCO2 slope cut-off was a reliable indicator of the composite endpoint in both sex categories (Figure 3b).
In our cohort, a percent of predicted pVO2 of less than 50% was present in 120 (26%) patients. This cut-off was associated with the composite endpoint (HR 4.355, 95% CI 2.694–7.039, p < 0.001). This cut-off showed a sensitivity of 60% and a specificity of 89% in females, while it had a sensitivity of 48% and specificity of 78% in males. As a result, the Youden index in females was higher than in male patients (J 0.49 vs. J 0.26). This threshold was a reliable discriminator in both subgroups according to the survival curve analysis (Supplementary Figure S2a).

3.5. Alternative Thresholds for pVO2 and VE/VCO2 Slope

In an assessment of potential alternative thresholds, a pVO2 ≤ 14 mL/kg/min yielded a higher Youden index in female patients compared to the pVO2 ≤ 12 mL/kg/min cut-off (J 0.60 vs. J 0.34) (Table 5). Similarly, a pVO2 ≤ 15 mL/kg/min value showed a higher overall diagnostic effectiveness in male patients compared to the traditional cut-off (J 0.36 vs. J 0.12). The predictive value of this cut-off was supported by the Kaplan–Meier analysis (Figure 4a).
In males, a VE/VCO2 slope threshold of > 32 demonstrated sensitivity of 78% and a specificity of 57%, exhibiting a higher Youden index than the traditional VE/VCO2 slope cut-off (J 0.35 vs. J 0.23). Regarding female patients, the traditional VE/VCO2 slope > 35 threshold was associated with the highest overall diagnostic effectiveness (J 0.65). Additionally, it was demonstrated in the survival analysis that these cut-off values accurately predicted worse outcomes (Figure 4b).
A percent of predicted pVO2 of ≤ 55% yielded a significantly higher Youden index in female patients compared to the ≤ 50% threshold (J 0.76 vs. J 0.49) while a percent of predicted pVO2 of ≤ 58% showed a higher Youden index in comparison to the traditional thresholds (J 0.32 vs. J 0.26) (Table 5). These cut-offs were accurate discriminators of the composite endpoint in both sex subgroups (log-rank p < 0.001) (Supplementary Figure S2b).

4. Discussion

Our study’s key conclusion was that the traditional CPET variables had a considerably higher predictive power for HF outcomes in women compared to men. Furthermore, the ISHLT recommended thresholds for pVO2 (≤12 mL/kg/min, or ≤ 14 mL/kg/min if intolerant to β-blockers), VE/VCO2 slope (>35), and percent of predicted pVO2 ≤ 55% showed a significantly higher overall diagnostic effectiveness in women compared to men. Additionally, our study assessed the predictive capacity of various CPET variables and proposed sex-specific cut-offs for pVO2, VE/VCO2 slope, and the percent of predicted pVO2, which may assist in a more precise risk assessment in women and men with HFrEF. However, one of the main limitations of our study was that 79% of the enrolled patients were male; thus, further studies should include a higher proportion of female patients.
The current evidence on the predictive value of CPET in women with HFrEF was evaluated in a recent article by the Heart Failure Association’s Committee on Exercise Physiology and Training [1,14,23,24,25,26,27]. The mean age of female patients enrolled in these studies was slightly lower than that of male patients, and one of the explanations for female underrepresentation in HfrEF trials was a larger proportion of older women who were excluded due to the policy of non-inclusion of elderly patients [28].
pVO2 is influenced by gender, age, motivation, pulmonary status, and muscle mass [29], which raised concerns that this parameter’s role as a prognostic indicator in female patients may lead to premature cardiac transplantation in women [14].
However, several observational studies showed that pVO2 is a reliable discriminator for HF events in female patients [24,26] and a large trial [12] showed that predictive pVO2 cut-offs for men and women with HfrEF should be independent. Although thresholds such as the GET were described to provide incremental value in the assessment of cardiorespiratory fitness in healthy controls [30,31], VO2 measured at GET did not show a significant prognostic power compared to pVO2 in our HfrEF cohort.
Women generally exhibited a lower corrected pVO2 than male patients; however, female patients presented a lower rate of HF events. Notably, female patients showed a nearly 10% higher percent of predicted pVO2 compared to men [23,27]. However, a study by Corrà et al. [11] postulated that HF outcomes in women may not actually be better than in men, as the female prognostic advantage is lost when sex-specific variations are properly taken into account with propensity score matching. Therefore, adjusting for sex-related characteristics should be undertaken. Indeed, female patients in our cohort showed a significantly higher percent of predicted pVO2 despite having a similar absolute pVO2 value, with a numerically inferior frequency of the composite endpoint.
The VE/VCO2 slope is an alternative CPET parameter with proven prognostic power, and the HF event risk is constant throughout a large range of VE/VCO2 slope values [32,33,34]. A study by Guazzi et al. [25] demonstrated that in both men and women with HfrEF, the predictive power of pVO2 and the VE/VCO2 slope are similar. Notably, the discriminative power of the VE/VCO2 slope was greater than that of pVO2 in female patients. Our findings are in keeping with these results, as the VE/VCO2 slope also showed a slightly higher prognostic power compared to pVO2 in the ROC curve analysis in the female subgroup.
The percent of predicted pVO2, an age- and gender-adjusted parameter assessing exercise capacity, was shown to stratify the risk for HF events with a higher accuracy compared to pVO2 in women [27]. The role of CPET in pre-surgical risk stratification in women has also been studied. In a study by Rose et al. [35,36], sex-specific CPET thresholds improved surgical risk stratification and thus may contribute to optimise clinical decision-making.
There is a paucity of randomized clinical trials evaluating the value of CPET variables in women with HFrEF. The HF-ACTION [10], a randomized trial with 2100 patients, also concluded that women presented a better clinical outcome, showing a lower pVO2 and a higher percent of predicted pVO2 compared to men. The parameter with the highest predictive power in women was the percent of predicted pVO2. Our study had similar findings, as the percent of predicted pVO2 was the CPET parameter with the highest predictive power for HF outcomes in the female subgroup. This result is in keeping with the ISHLT guidelines [7], which recommend that alternative parameters such as percent of predicted pVO2 may be considered in conjunction with pVO2 to guide HTx listing in female patients. In contrast, in the male subgroup, pVO2 and percent of predicted pVO2 had a similar prognostic power for risk stratification of HF events. Moreover, our study showed that, in a cohort with similar pVO2 values between sexes, the predictive power of the traditional CPET parameters was notably lower in men than in women, which is in contrast with the results reported in a previous trial by Elmariah et al. [14].
The position paper by Corrà et al. [1] proposes three different threshold values of pVO2 for male HFrEF patients: <10 mL/kg/min, 10 to 18, or > 18 mL/kg/min. However, there are still limited data to define an accurate cut-off for other subgroups of patients, women or elderly patients in particular [9]. Extrapolating these three advocated thresholds of pVO2 in male patients to other subgroups may lead to misconceptions and inaccuracies of the objective pVO2 [9]. Consequently, further studies are necessary to define an accurate threshold to guide patient selection for HTx listing in women.
In a trial by Green et al. [26], the proposed pVO2 thresholds in females with HFrEF for high- (≤10), medium- (10.1 to 14), and low-risk (>14 mL/kg/min) showed a one-year event-free survival of 80%, 84%, and 93%, respectively. Elmariah et al. [14] reported that, in the current era, HTx may be deferred if the pVO2 is over 10 mL/kg/min. However, this study had several disparities between sexes in the baseline characteristics and it did not consider patients with CRT, which can affect pVO2 values [37,38].
In our cohort, the ISHLT recommended thresholds of pVO2 ≤ 12 or ≤ 14 mL/kg/min, VE/VCO2 slope > 35 and percent of predicted pVO2 ≤ 50% showed a higher overall diagnostic effectiveness in women compared to men, in keeping with the higher prognostic power these parameters showed in female patients. We proposed alternative thresholds that may improve risk discrimination among female patients. A threshold of pVO2 ≤ 14 mL/kg/min (including patients on β-blockers) and a percent of predicted pVO2 ≤ 55% showed a slightly lower specificity but a higher sensitivity, with an overall higher overall diagnostic effectiveness. The recommended cut-off of VE/VCO2 slope > 35 was the strongest predictor of HF events in women. Regarding male patients, a pVO2 threshold of ≤ 15 mL/kg/min, a VE/VCO2 slope of > 32, and a percent of predicted pVO2 of ≤ 58% may also provide an improved diagnostic effectiveness compared to the traditional thresholds for HTx listing.
Circulatory power is a surrogate of left ventricular stroke work index, incorporating pVO2 and peak systolic BP [39]. Circulatory power was a significant predictor of HF events in our cohort, especially in males. In a recent study by Martinez et al. [40] evaluating patients with advanced HF, circulatory power presented the highest discriminative power for HF outcomes and mortality, concluding that this parameter should also be considered for risk stratification in conjunction with the traditional CPET variables. However, further research is needed to determine whether circulatory power can contribute to the decision of the optimal timing for HTx in women.
In our study, peak O2 pulse presented a significantly higher predictive power for HF outcomes in women. Peak O2 pulse, a non-invasive measure of stroke volume and arteriovenous O2 differential, represents the pVO2 corrected for HRt [41]. Several CPET measures, including the pVO2, are corrected for total weight rather than lean body mass. There is a high variability in body fat as a percentage of total body weight [27,42,43] which can also contribute to the lower pVO2 reported in women [44]. The use of corrected pVO2 adjusted for lean body mass may be a more accurate measurement of exercise intolerance, particularly in groups with a greater body fat percentage such as women [42,45].
A trial by Lavie [45] et al. found that pVO2 lean and peak O2 pulse lean outperformed pVO2 as predictors of major HF events, including among obese patients and women. The authors noted that, when combined with conventional CPET variables, peak O2 pulse and lean body mass-adjusted O2 pulse were powerful predictors of HF outcomes in patients with HFrEF, particularly in populations with a higher percent of body fat.
Prognostic risk scores such as a high to medium risk HFSS [21] or a Seattle Heart Failure Model (SHFM) [46,47] <80% are also recommended as alternative parameters to consider HTx listing [7]. Although the SHFM was also an accurate predictor of HF outcomes in female patients, the Meta-Analysis Global Group in Chronic Heart Failure (MAGGIC) [48,49] score outperformed the SHFM owing to improved risk classification, presenting a similar discriminatory ability in both sexes, despite an overestimation of death in female patients at the 3-year follow-up [50].

Limitations

Firstly, since this is a retrospective study of a prospective database in our center, our findings require confirmation in larger, randomized studies. Additionally, the majority of patients enrolled were men (79%), which is a high proportion, particularly in a study evaluating disparities between female and male patients.
Secondly, our study enrolled unmatched patient subgroups. However, consecutive patient enrollment attenuated the lack of randomization. Furthermore, most baseline characteristics were comparable between sex groups. Men had a numerically higher proportion of ischemic HF, although there was no statistical difference among groups.
Most patients in each subgroup were receiving optimal disease-modifying pharmacological therapies for HF. However, only 10% of patients in our cohort were taking SGLT2i, as they were included in our study between 2009 and 2018, when this drug class was not yet considered as an optimised standard of care medication for patients without diabetes [51]. Future studies should include more patients taking SGLT2i, as they have shown to significantly reduce HF events [52,53]. Less than 25% of patients were taking angiotensin receptor/neprilysin inhibitors, as this therapy was not available for patients enrolled before 2016. Future trials should include more patients receiving sacubitril/valsartan. Moreover, new therapies such as selective cardiac myosin activators or guanylate cyclase stimulators were not available at the time of patient enrolment.
Thirdly, our research lacked the statistical power to infer a new pVO2 threshold for patients who were intolerant to β-blockers, as the majority of patients (86%) were taking β-blockers. As a result, the proposed pVO2 thresholds might not be reproducible in this subgroup of HFrEF patients.
Our study only included patients who had a maximal CPET. There is no current agreement on the best peak RER cut-off to determine maximal effort, especially in patients with HFrEF. A number of cut-offs ranging from 1.0 to 1.10 were suggested [1,54,55,56]. As our aim was to assess the recommendations for HTx, a peak RER of 1.05 was considered to determine a maximal CPET, as recommended by the ISHLT [7]. Consequently, our proposed cut-offs might not be applicable to an HF population with submaximal exercise testing, particularly considering the lower prognostic power of pVO2 in submaximal exercise capacity [57]. In patients with submaximal exercise capacity, VE/VCO2 slope and percent of predicted pVO2 may assist in the clinical stratification [7,34,55]. Indeed, the reliability of RER-based assessment of maximal exercise is suboptimal as there are methodological issues thwarting the accurate assessment of VO2max in submaximal exercise. Pool and Jones [58] caution against the acceptance of pVO2 measured during ramp incremental exercise as a maximum value in patients with submaximal exercise and proposed the inclusion of a second short constant work rate CPET, completed at a higher work rate than that previously achieved during the ramp test, in order to accurately verify the VO2max. Therefore, serial CPET may be more informative than a single cardiopulmonary exercise test and thus provide a more accurate assessment of the VO2max.
Our study evaluated the GET, as described by Beaver et al. [59]. However, it is now recognized that insufficient O2 is not the primary basis for lactatemia. Critical power likely represents the threshold above which there is a sustained glycolytic contribution with lactate accumulation. Although lactate is a key energy source, there is no evidence that the muscle becomes dysoxic or anoxic [60]. Thus, instead of the GET, critical power may potentially be a more accurate predictor of exercise capacity [60].
Lastly, our study had a lower rate of HF outcomes, especially urgent HTx, compared to other studies [32]. As all the recruited patients were reviewed by the specialized Heart Failure team for a possible indication for HTx, our results may not be applicable to the overall HFrEF population encompassing older patients or patients with significant comorbidities.

5. Conclusions

In an HFrEF cohort undergoing CPET, pVO2, VE/VCO2 slope, and the percent of predicted pVO2 were the variables with the highest discriminative power for HF events, with a higher predictive power in female patients compared to male patients. The ISHLT guideline thresholds for pVO2 and VE/VCO2 slope showed a higher diagnostic effectiveness in women. Sex-specific pVO2, VE/VCO2 slope, and percent of predicted pVO2 cut-offs presented a higher prognostic power than the recommended thresholds. Our results indicate that sex-specific cut-offs may assist in patient selection for HTx. However, more data are necessary to help close the gap in evidence between sexes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life13101985/s1, Figure S1: ROC curves for the composite endpoint in a 36-month follow up. (a) Percent of predicted peak oxygen consumption (pVO2) in female patients. (b) Percent of predicted pVO2 in male patients; Figure S2: Survival analysis for the composite endpoint in female patients and male patients according to (a) the International Society for Heart and Lung Transplantation (ISHLT) thresholds of percent of predicted peak O2 consumption (pVO2) ≤ 50% and (b) thresholds of percent of predicted pVO2 ≤ 55% in females and percent of predicted pVO2 ≤ 58% in males.

Author Contributions

Conceptualization, P.G.B. and A.V.G.; methodology, J.F.R. and P.G.B.; software, S.S. and A.V.G.; validation, P.R.; formal analysis, P.G.B.; investigation, R.I.M. and T.P.-d.-S.; resources, A.T.T. and J.F.R.; data curation, P.G.B.; writing—original draft preparation, P.G.B., J.F.R.; writing—review and editing, A.T.T., T.P.-d.-S., A.V.G., R.M.S. and P.G.B.; visualization, R.I.M.; supervision, R.M.S.; project administration, R.M.S. and R.C.F.; funding acquisition, none. All authors have read and agreed to the published version of the manuscript.

Funding

This study did not receive any funding.

Institutional Review Board Statement

The National Committee for Patient Information Protection and the Institutional Ethics Committee (authorization number 1232/2022 by the Comissão Nacional de Protecção de Dados) approved the study, which was conducted in conformity with the Declaration of Helsinki. The study did not involve any experiments on humans or animals.

Informed Consent Statement

All participants in the study provided written informed consent for publishing of this study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to patient consent regarding availability of individual patient data, applicable only to the local investigation team.

Acknowledgments

The authors would like to acknowledge the invaluable contribution of the Heart Failure team and the Cardiorespiratory Exercise Testing Laboratory team of Santa Marta Hospital in this study.

Conflicts of Interest

The authors declare no conflict of interest. The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.

References

  1. Corrà, U.; Agostoni, P.G.; Anker, S.D.; Coats, A.J.; Leiro, M.G.C.; de Boer, R.A.; Harjola, V.-P.; Hill, L.; Lainscak, M.; Lund, L.H.; et al. Role of cardiopulmonary exercise testing in clinical stratification in heart failure. A position paper from the Committee on Exercise Physiology and Training of the Heart Failure Association of the European Society of Cardiology. Eur. J. Heart Fail. 2018, 20, 3–15. [Google Scholar] [CrossRef]
  2. Chua, T.P.; Ponikowski, P.; Harrington, D.; Anker, S.D.; Webb-Peploe, K.; Clark, A.L.; A Poole-Wilson, P.; Coats, A.J. Clinical Correlates and Prognostic Significance of the Ventilatory Response to Exercise in Chronic Heart Failure. J. Am. Coll. Cardiol. 1997, 29, 1585–1590. [Google Scholar] [CrossRef]
  3. Francis, D.; Shamim, W.; Davies, L.; Piepoli, M.; Ponikowski, P.; Anker, S.; Coats, A. Cardiopulmonary exercise testing for prognosis in chronic heart failure: Continuous and independent prognostic value from VE/VCO2slope and peak VO2. Eur. Heart J. 2000, 21, 154–161. [Google Scholar] [CrossRef]
  4. Arena, R.; Humphrey, R. Comparison of ventilatory expired gas parameters used to predict hospitalization in patients with heart failure. Am. Heart J. 2002, 143, 427–432. [Google Scholar] [CrossRef]
  5. Arena, R.; Myers, J.; Aslam, S.S.; Varughese, E.B.; Peberdy, M.A. Peak VO2 and VE/VCO2 slope in patients with heart failure: A prognostic comparison. Am. Heart J. 2004, 147, 354–360. [Google Scholar] [CrossRef]
  6. Kleber, F.X.; Vietzke, G.; Wernecke, K.D.; Bauer, U.; Opitz, C.; Wensel, R.; Sperfeld, A.; Glaser, S. Impairment of Ventilatory Efficiency in Heart Failure. Circulation 2000, 101, 2803–2809. [Google Scholar] [CrossRef]
  7. Mehra, M.R.; Canter, C.E.; Hannan, M.M.; Semigran, M.J.; Uber, P.A.; Baran, D.A.; Danziger-Isakov, L.; Kirklin, J.K.; Kirk, R.; Kushwaha, S.S.; et al. The 2016 International Society for Heart Lung Transplantation listing criteria for heart transplantation: A 10-year update. J. Heart Lung Transpl. 2016, 35, 1–23. [Google Scholar] [CrossRef]
  8. Peterson, L.R.; Schechtman, K.B.; Ewald, G.A.; Geltman, E.M.; Fuentes, L.d.L.; Meyer, T.; Krekeler, P.; Moore, M.L.; Rogers, J.G. Timing of cardiac transplantation in patients with heart failure receiving β-adrenergic blockers. J. Heart Lung Transpl. 2003, 22, 1141–1148. [Google Scholar] [CrossRef]
  9. Palau, P.; Domínguez, E.; Núñez, J. Sex differences on peak oxygen uptake in heart failure. ESC Heart Fail. 2019, 6, 921–926. [Google Scholar] [CrossRef]
  10. Keteyian, S.J.; Patel, M.; Kraus, W.E.; Brawner, C.A.; McConnell, T.R.; Piña, I.L.; Leifer, E.S.; Fleg, J.L.; Blackburn, G.; Fonarow, G.C.; et al. Variables Measured During Cardiopulmonary Exercise Testing as Predictors of Mortality in Chronic Systolic Heart Failure. J. Am. Coll. Cardiol. 2016, 67, 780–789. [Google Scholar] [CrossRef]
  11. Corrà, U.; Agostoni, P.; Giordano, A.; Cattadori, G.; Battaia, E.; La Gioia, R.; Scardovi, A.B.; Emdin, M.; Metra, M.; Sinagra, G.; et al. Sex Profile and Risk Assessment with Cardiopulmonary Exercise Testing in Heart Failure: Propensity Score Matching for Sex Selection Bias. Can. J. Cardiol. 2015, 32, 754–759. [Google Scholar] [CrossRef]
  12. Ehrman, J.K.; Brawner, C.A.; Shafiq, A.; Lanfear, D.E.; Saval, M.; Keteyian, S.J. Cardiopulmonary Exercise Measures of Men and Women with HFrEF Differ in Their Relationship to Prognosis: The Henry Ford Hospital Cardiopulmonary Exercise Testing (FIT-CPX) Project. J. Card. Fail. 2018, 24, 227–233. [Google Scholar] [CrossRef]
  13. Piña, I.L.; Kokkinos, P.; Kao, A.; Bittner, V.; Saval, M.; Clare, B.; Goldberg, L.; Johnson, M.; Swank, A.; Ventura, H.; et al. Baseline differences in the HF-ACTION trial by sex. Am. Heart J. 2009, 158, S16–S23. [Google Scholar] [CrossRef]
  14. Elmariah, S.; Goldberg, L.R.; Allen, M.T.; Kao, A. Effects of Gender on Peak Oxygen Consumption and the Timing of Cardiac Transplantation. J. Am. Coll. Cardiol. 2006, 47, 2237–2242. [Google Scholar] [CrossRef]
  15. Bruce, R.A.; Kusumi, F.; Hosmer, D. Maximal oxygen intake and nomographic assessment of functional aerobic impairment in cardiovascular disease. Am. Heart J. 1973, 85, 546–562. [Google Scholar] [CrossRef]
  16. Beale, A.L.; Meyer, P.; Marwick, T.H.; Lam, C.S.; Kaye, D.M. Sex Differences in Cardiovascular Pathophysiology. Circulation 2018, 138, 198–205. [Google Scholar] [CrossRef]
  17. Redfield, M.M.; Jacobsen, S.J.; Borlaug, B.A.; Rodeheffer, R.J.; Kass, D.A. Age- and Gender-Related Ventricular-Vascular Stiffening. Circulation 2005, 112, 2254–2262. [Google Scholar] [CrossRef]
  18. Martens, P.; Nijst, P.; Verbrugge, F.H.; Smeets, K.; Dupont, M.; Mullens, W. Impact of iron deficiency on exercise capacity and outcome in heart failure with reduced, mid-range and preserved ejection fraction. Acta Cardiol. 2017, 73, 115–123. [Google Scholar] [CrossRef]
  19. Karastergiou, K.; Smith, S.R.; Greenberg, A.S.; Fried, S.K. Sex differences in human adipose tissues—The biology of pear shape. Biol. Sex Differ. 2012, 3, 13. [Google Scholar] [CrossRef]
  20. Jessup, M.; Piña, I.L. Is it important to examine gender differences in the epidemiology and outcome of severe heart failure? J. Thorac. Cardiovasc. Surg. 2004, 127, 1247–1252. [Google Scholar] [CrossRef]
  21. Aaronson, K.D.; Schwartz, J.S.; Chen, T.-M.; Wong, K.-L.; Goin, J.E.; Mancini, D.M. Development and Prospective Validation of a Clinical Index to Predict Survival in Ambulatory Patients Referred for Cardiac Transplant Evaluation. Circulation 1997, 95, 2660–2667. [Google Scholar] [CrossRef]
  22. DeLong, E.R.; DeLong, D.M.; Clarke-Pearson, D.L. Comparing the Areas under Two or More Correlated Receiver Operating Characteristic Curves: A Nonparametric Approach. Biometrics 1988, 44, 837–845. [Google Scholar] [CrossRef]
  23. Corrà, U.; Mezzani, A.; Giordano, A.; Pistono, M.; Gnemmi, M.; Caruso, R.; Giannuzzi, P. Peak oxygen consumption and prognosis in heart failure. Int. J. Cardiol. 2013, 167, 157–161. [Google Scholar] [CrossRef]
  24. Hsich, E.; Chadalavada, S.; Krishnaswamy, G.; Starling, R.C.; Pothier, C.E.; Blackstone, E.H.; Lauer, M.S. Long-Term Prognostic Value of Peak Oxygen Consumption in Women Versus Men with Heart Failure and Severely Impaired Left Ventricular Systolic Function. Am. J. Cardiol. 2007, 100, 291–295. [Google Scholar] [CrossRef]
  25. Guazzi, M.; Arena, R.; Myers, J. Comparison of the prognostic value of cardiopulmonary exercise testing between male and female patients with heart failure. Int. J. Cardiol. 2006, 113, 395–400. [Google Scholar] [CrossRef]
  26. Green, P.; Lund, L.H.; Mancini, D. Comparison of Peak Exercise Oxygen Consumption and the Heart Failure Survival Score for Predicting Prognosis in Women Versus Men. Am. J. Cardiol. 2007, 99, 399–403. [Google Scholar] [CrossRef]
  27. Richards, D.R.; Mehra, M.R.; O Ventura, H.; Lavie, C.J.; Smart, F.W.; Stapleton, D.D.; Milani, R.V. Usefulness of Peak Oxygen Consumption in Predicting Outcome of Heart Failure in Women Versus Men. Am. J. Cardiol. 1997, 80, 1236–1238. [Google Scholar] [CrossRef]
  28. Sciomer, S.; Moscucci, F.; Salvioni, E.; Marchese, G.; Bussotti, M.; Corrà, U.; Piepoli, M.F. Role of gender, age and BMI in prognosis of heart failure. Eur. J. Prev. Cardiol. 2020, 27, 46–51. [Google Scholar] [CrossRef]
  29. Stelken, A.M.; Younis, L.T.; Jennison, S.H.; Miller, D.D.; Miller, L.W.; Shaw, L.J.; Kargl, D.; Chaitman, B.R. Prognostic value of cardiopulmonary exercise testing using percent achieved of predicted peak oxygen uptake for patients with ischemic and dilated cardiomyopathy. J. Am. Coll. Cardiol. 1996, 27, 345–352. [Google Scholar] [CrossRef]
  30. Meyler, S.; Bottoms, L.; Wellsted, D.; Muniz-Pumares, D. Variability in exercise tolerance and physiological responses to exercise prescribed relative to physiological thresholds and to maximum oxygen uptake. Exp. Physiol. 2023, 108, 581–594. [Google Scholar] [CrossRef]
  31. Collins, J.; Leach, O.; Dorff, A.; Linde, J.; Kofoed, J.; Sherman, M.; Proffit, M.; Gifford, J.R. Critical power and work-prime account for variability in endurance training adaptations not captured by Vo2max. J. Appl. Physiol. 2022, 133, 986–1000. [Google Scholar] [CrossRef] [PubMed]
  32. Arena, R.; Myers, J.; Abella, J.; Peberdy, M.A.; Bensimhon, D.; Chase, P.; Guazzi, M. Development of a Ventilatory Classification System in Patients with Heart Failure. Circulation 2007, 115, 2410–2417. [Google Scholar] [CrossRef] [PubMed]
  33. Tabet, J.-Y.; Beauvais, F.; Thabut, G.; Tartière, J.-M.; Logeart, D.; Cohen-Solal, A. A critical appraisal of the prognostic value of the VE/VCO2 slope in chronic heart failure. Eur. J. Cardiovasc. Prev. Rehabil. 2003, 10, 267–272. [Google Scholar] [CrossRef] [PubMed]
  34. Ingle, L.; Goode, K.; Carroll, S.; Sloan, R.; Boyes, C.; Cleland, J.G.; Clark, A.L. Prognostic value of the VE/VCO2 slope calculated from different time intervals in patients with suspected heart failure. Int. J. Cardiol. 2007, 118, 350–355. [Google Scholar] [CrossRef] [PubMed]
  35. Rose, G.A.; Davies, R.G.; Torkington, J.; Berg, R.M.G.; Appadurai, I.R.; Poole, D.C.; Bailey, D.M. Assessing cardiorespiratory fitness relative to sex improves surgical risk stratification. Eur. J. Clin. Investig. 2023, 53, e13981. [Google Scholar] [CrossRef] [PubMed]
  36. Rose, G.A.; Davies, R.G.; Appadurai, I.R.; Williams, I.M.; Bashir, M.; Berg, R.M.G.; Poole, D.C.; Bailey, D.M. ‘Fit for surgery’: The relationship between cardiorespiratory fitness and postoperative outcomes. Exp. Physiol. 2022, 107, 787–799. [Google Scholar] [CrossRef]
  37. Anand, I.S.; Carson, P.; Galle, E.; Song, R.; Boehmer, J.; Ghali, J.K.; Jaski, B.; Lindenfeld, J.; O’Connor, C.; Steinberg, J.S.; et al. Cardiac Resynchronization Therapy Reduces the Risk of Hospitalizations in Patients with Advanced Heart Failure. Circulation 2009, 119, 969–977. [Google Scholar] [CrossRef]
  38. Abraham, W.T.; Fisher, W.G.; Smith, A.L.; Delurgio, D.B.; Leon, A.R.; Loh, E.; Kocovic, D.Z.; Packer, M.; Clavell, A.L.; Hayes, D.L.; et al. Cardiac Resynchronization in Chronic Heart Failure. N. Engl. J. Med. 2002, 346, 1845–1853. [Google Scholar] [CrossRef]
  39. Cohen-Solal, A.; Tabet, J.; Logeart, D.; Bourgoin, P.; Tokmakova, M.; Dahan, M. A non-invasively determined surrogate of cardiac power (‘circulatory power’) at peak exercise is a powerful prognostic factor in chronic heart failure. Eur. Heart J. 2002, 23, 806–814. [Google Scholar] [CrossRef]
  40. Martinez, L.; Cacciottolo, P.; Barnes, J.; Sylvester, K.; Oates, K.; Kydd, A.; Lewis, C.; Parameshwar, J.; Pettit, S.; Bhagra, S. Circulatory Power is Superior to Peak Oxygen Consumption in Predicting Adverse Outcomes in Ambulatory Patients Assessed for Heart Transplantation. J. Heart Lung Transplant. 2022, 41, S214. [Google Scholar] [CrossRef]
  41. Cohen-Solal, A.; Barnier, P.; Pessione, F.; Seknadji, P.; Logeart, D.; Laperche, T.; Gourgon, R. Comparison of the long term prognostic value of peak exercise oxygen pulse and peak oxygen uptake in patients with chronic heart failure. Heart 1997, 78, 572–576. [Google Scholar] [CrossRef]
  42. Osman, A.F.; Mehra, M.R.; Lavie, C.J.; Nunez, E.; Milani, R.V. The incremental prognostic importance of body fat adjusted peak oxygen consumption in chronic heart failure. J. Am. Coll. Cardiol. 2000, 36, 2126–2131. [Google Scholar] [CrossRef] [PubMed]
  43. Milani, R.V.; Lavie, C.J. The Effects of Body Composition Changes to Observed Improvements in Cardiopulmonary Parameters After Exercise Training with Cardiac Rehabilitation. Chest 1998, 113, 599–601. [Google Scholar] [CrossRef] [PubMed]
  44. Toth, M.J.; Gardner, A.W.; Ades, P.A.; Poehlman, E.T. Contribution of body composition and physical activity to age-related decline in peak VO2 in men and women. J. Appl. Physiol. 1994, 77, 647–652. [Google Scholar] [CrossRef] [PubMed]
  45. Lavie, C.J.; Milani, R.V.; Mehra, M.R. Peak exercise oxygen pulse and prognosis in chronic heart failure. Am. J. Cardiol. 2004, 93, 588–593. [Google Scholar] [CrossRef]
  46. Gorodeski, E.Z.; Chu, E.C.; Chow, C.H.; Levy, W.C.; Hsich, E.; Starling, R.C. Application of the Seattle Heart Failure Model in Ambulatory Patients Presented to an Advanced Heart Failure Therapeutics Committee. Circ. Heart Fail. 2010, 3, 706–714. [Google Scholar] [CrossRef]
  47. Kalogeropoulos, A.P.; Georgiopoulou, V.V.; Giamouzis, G.; Smith, A.L.; Agha, S.A.; Waheed, S.; Laskar, S.; Puskas, J.; Dunbar, S.; Vega, D.; et al. Utility of the Seattle Heart Failure Model in Patients with Advanced Heart Failure. J. Am. Coll. Cardiol. 2009, 53, 334–342. [Google Scholar] [CrossRef]
  48. Pocock, S.J.; Ariti, C.A.; McMurray, J.J.; Maggioni, A.; Køber, L.; Squire, I.B.; Swedberg, K.; Dobson, J.; Poppe, K.K.; Whalley, G.A.; et al. Predicting survival in heart failure: A risk score based on 39 372 patients from 30 studies. Eur. Heart J. 2012, 34, 1404–1413. [Google Scholar] [CrossRef]
  49. Sartipy, U.; Dahlström, U.; Edner, M.; Lund, L.H. Predicting survival in heart failure: Validation of the MAGGIC heart failure risk score in 51 043 patients from the Swedish Heart Failure Registry. Eur. J. Heart Fail. 2013, 16, 173–179. [Google Scholar] [CrossRef]
  50. Vishram-Nielsen, J.K.; Foroutan, F.; Ross, H.J.; Gustafsson, F.; Alba, A.C. Performance of Prognostic Risk Scores in Heart Failure Patients: Do Sex Differences Exist? Can. J. Cardiol. 2019, 36, 45–53. [Google Scholar] [CrossRef]
  51. Ponikowski, P.; Voors, A.A.; Anker, S.D.; Bueno, H.; Cleland, J.G.F.; Coats, A.J.S.; Falk, V.; González-Juanatey, J.R.; Harjola, V.-P.; Jankowska, E.A.; et al. 2016 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure. Eur. J. Heart Fail. 2016, 18, 891–975. [Google Scholar] [CrossRef] [PubMed]
  52. McMurray, J.J.V.; Solomon, S.D.; Inzucchi, S.E.; Køber, L.; Kosiborod, M.N.; Martinez, F.A.; Ponikowski, P.; Sabatine, M.S.; Anand, I.S.; Bělohlávek, J.; et al. Dapagliflozin in Patients with Heart Failure and Reduced Ejection Fraction. N. Engl. J. Med. 2019, 381, 1995–2008. [Google Scholar] [CrossRef] [PubMed]
  53. Packer, M.; Anker, S.D.; Butler, J.; Filippatos, G.; Pocock, S.J.; Carson, P.; Januzzi, J.; Verma, S.; Tsutsui, H.; Brueckmann, M.; et al. Cardiovascular and Renal Outcomes with Empagliflozin in Heart Failure. N. Engl. J. Med. 2020, 383, 1413–1424. [Google Scholar] [CrossRef] [PubMed]
  54. Guazzi, M.; Arena, R.; Halle, M.; Piepoli, M.F.; Myers, J.; Lavie, C.J. 2016 focused update: Clinical recommendations for cardiopulmonary exercise testing data assessment in specific patient populations. Eur. Heart J. 2016, 39, 1144–1161. [Google Scholar] [CrossRef]
  55. Malhotra, R.; Bakken, K.; D’elia, E.; Lewis, G.D. Cardiopulmonary Exercise Testing in Heart Failure. JACC Heart Fail. 2016, 4, 607–616. [Google Scholar] [CrossRef]
  56. Guazzi, M.; Bandera, F.; Ozemek, C.; Systrom, D.; Arena, R. Cardiopulmonary Exercise Testing: What is its value? J. Am. Coll. Cardiol. 2017, 70, 1618–1636. [Google Scholar] [CrossRef]
  57. Balady, G.J.; Arena, R.; Sietsema, K.; Myers, J.; Coke, L.; Fletcher, G.F.; Forman, D.; Franklin, B.; Guazzi, M.; Gulati, M.; et al. Clinician’s Guide to Cardiopulmonary Exercise Testing in Adults. Circulation 2010, 122, 191–225. [Google Scholar] [CrossRef]
  58. Poole, D.C.; Jones, A.M. Measurement of the maximum oxygen uptake Vo2max: Vo2peak is no longer acceptable. J. Appl. Physiol. 2017, 122, 997–1002. [Google Scholar] [CrossRef]
  59. Beaver, W.L.; Inglis, E.C.; Iannetta, D.; Murias, J.M.; Salvadego, D.; Sartorio, A.; Agosti, F.; Tringali, G.; Patrizi, A.; Isola, M.; et al. A new method for detecting anaerobic threshold by gas exchange. J. Appl. Physiol. 1986, 60, 2020–2027. [Google Scholar] [CrossRef]
  60. Poole, D.C.; Rossiter, H.B.; Brooks, G.A.; Gladden, L.B. The anaerobic threshold: 50+ years of controversy. J. Physiol. 2020, 599, 737–767. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the study population. * in patients intolerant to β-blockers. HFrEF: Heart failure with reduced ejection fraction; LVEF: Left ventricular ejection fraction; NYHA: New York Heart Association; CPET: Cardiopulmonary exercise test; HTx: Heart transplantation; ISHLT: International Society for Heart and Lung Transplantation; pVO2: Peak oxygen consumption; VE/VCO2 slope: Minute ventilation–carbon dioxide production relationship.
Figure 1. Flowchart of the study population. * in patients intolerant to β-blockers. HFrEF: Heart failure with reduced ejection fraction; LVEF: Left ventricular ejection fraction; NYHA: New York Heart Association; CPET: Cardiopulmonary exercise test; HTx: Heart transplantation; ISHLT: International Society for Heart and Lung Transplantation; pVO2: Peak oxygen consumption; VE/VCO2 slope: Minute ventilation–carbon dioxide production relationship.
Life 13 01985 g001
Figure 2. ROC curves for the composite endpoint in a 36-month follow up. (a) Peak oxygen consumption (pVO2) in female patients. (b) pVO2 in male patients. (c) Minute ventilation–carbon dioxide production relationship (VE/VCO2 slope) in female patients. (d) VE/VCO2 slope in male patients.
Figure 2. ROC curves for the composite endpoint in a 36-month follow up. (a) Peak oxygen consumption (pVO2) in female patients. (b) pVO2 in male patients. (c) Minute ventilation–carbon dioxide production relationship (VE/VCO2 slope) in female patients. (d) VE/VCO2 slope in male patients.
Life 13 01985 g002
Figure 3. Survival analysis for the composite endpoint according to the International Society for Heart and Lung Transplantation (ISHLT) thresholds in female patients and male patients. (a) Peak oxygen consumption (pVO2) ≤ 12 mL/Kg/min (≤14 mL/kg/min if intolerant to β-blockers [βB]). (b) Minute ventilation–carbon dioxide production ratio (VE/VCO2 slope) of >35.
Figure 3. Survival analysis for the composite endpoint according to the International Society for Heart and Lung Transplantation (ISHLT) thresholds in female patients and male patients. (a) Peak oxygen consumption (pVO2) ≤ 12 mL/Kg/min (≤14 mL/kg/min if intolerant to β-blockers [βB]). (b) Minute ventilation–carbon dioxide production ratio (VE/VCO2 slope) of >35.
Life 13 01985 g003
Figure 4. Survival analysis for the composite endpoint in female patients and male patients according to (a) Peak O2 consumption (pVO2) of ≤ 14 and ≤ 15 mL/Kg/min, respectively. (b) Minute ventilation–CO2 production ratio (VE/VCO2 slope) of > 35 and > 32, respectively.
Figure 4. Survival analysis for the composite endpoint in female patients and male patients according to (a) Peak O2 consumption (pVO2) of ≤ 14 and ≤ 15 mL/Kg/min, respectively. (b) Minute ventilation–CO2 production ratio (VE/VCO2 slope) of > 35 and > 32, respectively.
Life 13 01985 g004
Table 1. Baseline characteristics of the study population (n = 458).
Table 1. Baseline characteristics of the study population (n = 458).
Overall
(n = 458)
Female
(n = 95)
Male
(n = 363)
p-Value
Clinical and demographic data
Age (years)56 ± 1254 ± 1456 ± 120.328
Body mass index (kg/m2)27.1 ± 4.326.3 ± 4.627.4 ± 4.20.335
Ischemic etiology (n, %)261 (57)47 (49)214 (59)0.092
ACEi/ARB (n, %)361 (79)77 (81)284 (78)0.199
ARNI (n, %)80 (17)13 (14)67 (18)0.273
β-blocker (n, %)392 (86)81 (85)311 (86)0.726
MRA (n, %)336 (73)72 (76)264 (73)0.789
iSGLT2 (n, %)47 (10)8 (8)39 (11)0.164
Digoxin (n, %)129 (28)23 (24)106 (29)0.372
Diabetes104 (23)15 (16)89 (25)0.094
CKD (n, %)145 (32)25 (26)120 (34)0.138
AF (n, %)109 (24)14 (15)95 (26)0.021
ICD * (n, %)293 (64)59 (62)234 (64)0.617
Cardiac resynchronization therapy (n, %)102 (22)27 (28)75 (21)0.128
NYHA II347 (76)74 (78)273 (75)0.485
NYHA III111 (24)21 (22)90 (25)0.485
HFSS8.6 ± 1.18.8 ± 0.98.6 ± 1.20.109
Laboratory data
eGFR, ml/min/1.73 m275.3 ± 29.277.1 ± 30.974.8 ± 28.70.517
Na+, mEq/l138.0 ± 3.0138.4 ± 2.8137.9 ± 3.10.108
N-terminal pro b-type natriuretic
peptide, pg/mL
2196 ± 21012204 ± 17242193 ± 20990.979
Echocardiographic data
LVEDD, mm/m267.4 ± 10.363.8 ± 9.768.0 ± 10.30.064
LVEF, %29.7 ± 8.031.3 ± 7.929.0 ± 7.50.213
Mitral regurgitation severity III–IV, %67 (14)19 (20)48 (13)0.097
Right ventricular dysfunction (n, %)69 (15)9 (10)60 (16)0.630
Exercise testing data
Peak Respiratory Exchange Ratio1.14 ± 0.071.13 ± 0.081.14 ± 0.070.566
Delta heart rate during exercise51 (37–68)48 (34–67)52 (38–69)0.819
HHR117 (11–27)19 (14–29)16 (11–26)0.058
pVO2, mL/kg/min18.5 ± 5.818.0 ± 5.618.6 ± 5.90.363
Percent of predicted pVO2 (%)63.8 ± 18.767.4 ± 16.762.8 ± 19.10.021
VE/VCO2 slope33.9 ± 9.633.0 ± 8.934.2 ± 9.80.246
pVO2, mL/kg/min at GET13.6 ± 4.610.9 ± 2.814.2 ± 4.70.001
Peak O2 pulse0.14 ± 0.060.13 ± 0.030.14 ± 0.070.149
Circulatory power2883 ± 15432715 ± 10352927 ± 16490.235
Ventilatory power4.8 ± 1.74.8 ± 1.54.8 ± 1.70.739
COP28.9 ± 7.229.5 ± 7.928.8 ± 7.00.630
PetCO2 at rest, mmHg33.6 ± 4.833.9 ± 5.133.5 ± 4.70.558
PetCO2 at GET, mmHg36.8 ± 6.037.5 ± 5.936.6 ± 6.10.262
* including patients with a cardiac resynchronization therapy device. CPET: Cardiopulmonary exercise test; ACEi: Angiotensin-converting enzyme inhibitors; ARNI: Angiotensin receptor neprilysin inhibitors; ARB: Angiotensin receptor blockers; MRA: Mineralocorticoid receptor antagonists; CKD: Chronic kidney disease; AF: Atrial fibrillation; ACEi: Angiotensin-converting enzyme inhibitors; ICD: Implantable cardioverter-defibrillator; HFSS: Heart Failure Survival Score; eGFR: estimated glomerular filtration rate; LVEF: Left ventricular ejection fraction; LVEDD: Left ventricular end-diastolic diameter; pVO2: Peak O2 consumption; VE/VCO2 slope: Minute ventilation-carbon dioxide production relationship; GET: Gas exchange threshold; COP: Cardiorespiratory optimal point; HRR1: Heart rate recovery in the first minute after finishing CPET; PetCO2: Partial pressure of end-tidal carbon dioxide.
Table 2. Total adverse events during follow-up.
Table 2. Total adverse events during follow-up.
Total Cohort (n = 458)Female
(n = 95)
Male
(n = 363)
p-Value
Composite endpoint (n, %)68 (14.8%)10 (10.5%)58 (16.0%)0.199
Total mortality (n, %)67 (14.6%)13 (13.7%)54 (14.9%)0.597
Cardiac mortality (n, %)54 (11.8%)8 (8.4%)46 (12.7%)0.098
Sudden cardiac death (n, %)19 (4.1%)2 (2.1%)17 (4.7%)0.147
Death from worsening HF (n, %)35 (7.6%)6 (6.3%)29 (7.9%)0.638
Urgent HTx (n, %)14 (3.1%)2 (2.1%)12 (3.3%)0.744
HF: Heart failure; HTx: Heart transplantation.
Table 3. Univariable and multivariable analysis of the composite endpoint.
Table 3. Univariable and multivariable analysis of the composite endpoint.
Total Cohort
ModelUnivariable HR95% CIp-valueMultivariable HR95% CIp-value
Male sex1.5470.791–3.0260.203
Age1.0020.983–1.0210.829
BMI0.9530.897–1.0130.1210.9540.887–1.0270.210
LVEF0.9270.900–0.955<0.0010.9350.905–0.966<0.001
eGFR0.9790.969–0.989<0.0010.9860.976–0.9960.009
Diabetes1.1960.254–5.6320.821
Smoker1.7161.405–2.8200.0331.3950.835–2.3280.203
Peak VO20.8350.789–0.883<0.0010.8560.804–0.912<0.001
Percent of predicted pVO20.9480.934–0.963<0.0010.9550.939–0.971<0.001
VE/VCO2 slope1.0581.041–1.075<0.0011.0641.039–1.090<0.001
Peak VO2 at GET, mL/kg/min0.8540.737–0.9890.0350.8790.687–1.1240.305
O2 pulse, mL/kg/beat0.8580.791–0.932<0.0010.8650.780–0.9610.007
Circulatory power, mmHg.mL/kg/min0.9990.999–0.999<0.0010.9990.998–1.000<0.001
Ventilatory power, mmHg0.5750.483–0.684<0.0010.6320.521–0.768<0.001
COP1.1181.054–1.186<0.0011.0600.956–1.1740.268
PetCO2 at rest, mmHg0.8870.839–0.937<0.0010.9480.889–1.0110.102
PetCO2 at GET, mmHg0.8620.826–0.900<0.0010.8900.845–0.993<0.001
Female sex
ModelUnivariable HR95% CIp-valueMultivariable HR95% CIp-value
Age1.0030.960–1.0480.888
BMI0.8970.770–1.0450.1620.8610.694–1.0670.171
LVEF0.8930.820–0.9730.0100.9410.864–1.0160.160
eGFR0.9770.952–1.0030.0860.9910.966–1.0160.459
Diabetes1.1350.629–2.0530.674
Smoker0.9400.199–4.4360.9371.5650.178–13.6990.686
Peak VO20.7040.583–0.850<0.0010.7460.604–0.9220.007
Percent of predicted pVO20.9110.875–0.948<0.0010.9130.858–0.9720.004
VE/VCO2 slope1.0931.052–1.135<0.0011.1431.039–1.2570.006
Peak VO2 at GET, mL/kg/min0.2230.010–5.1590.350
O2 pulse, mL/kg/beat0.4930.346–0.703<0.0010.4580.261–0.8020.006
Circulatory power, mmHg.mL/kg/min0.9980.997–0.9990.0020.9990.998–1.0000.069
Ventilatory power, mmHg0.4050.240–0.6840.0010.5650.297–1.0720.080
COP1.7750.100–3.4500.903
PetCO2 at rest, mmHg0.9030.792–1.0280.1230.9810.841–1.1440.807
PetCO2 at GET, mmHg0.8140.715–0.9270.0020.8710.736 –1.0310.108
Male sex
ModelUnivariable HR95% CIp-valueMultivariable HR95% CIp-value
Age1.0010.979–1.0220.963
BMI0.9600.898–1.0270.240
LVEF0.9330.905–0.963<0.0010.9380.905–0.971<0.001
eGFR0.9800.969–0.991<0.0010.9870.976–0.9980.020
Diabetes1.2110.639–2.2300.558
Smoker1.7911.024–3.1330.0411.4250.805–2.5210.224
Peak VO20.8540.806–0.905<0.0010.8690.813–0.928<0.001
Percent of predicted pVO20.9560.941–0.971<0.0010.9600.943–0.977<0.001
VE/VCO2 slope1.0511.032–1.070<0.0011.0561.030–1.084<0.001
Peak VO2 at GET, mL/kg/min0.8620.746–0.9960.0440.8800.691–1.1210.302
O2 pulse, mL/kg/beat0.8730.802–0.9490.0010.8840.794–0.9850.026
Circulatory power, mmHg.mL/kg/min0.9990.999–0.999<0.0010.9990.999–1.000<0.001
Ventilatory power, mmHg0.6110.510–0.733<0.0010.6450.526–0.792<0.001
COP1.0951.027–1.1670.0051.0620.962–1.1730.230
PetCO2 at rest, mmHg0.8860.834–0.942<0.0010.9370.873–1.0050.070
PetCO2 at GET, mmHg0.8700.831–0.911<0.0010.8870.839–0.939<0.001
BMI: Body mass index; eGFR: Estimated glomerular filtration rate; LVEF: Left ventricular ejection fraction; pVO2: Peak oxygen consumption; VE/VCO2 slope: Minute ventilation–carbon dioxide production relationship; GET: Gas exchange threshold; COP: Cardiorespiratory optimal point; PetCO2: Partial pressure of end-tidal carbon dioxide.
Table 4. Receiver operating characteristic (ROC) curve analysis of the composite endpoint.
Table 4. Receiver operating characteristic (ROC) curve analysis of the composite endpoint.
Female (n = 95)Male (n = 363)
CPET ParametersAUC95% CIp-ValueAUC95% CIp-Valuep-Value
(Interaction)
pVO2, mL/kg/min0.8490.740–0.958<0.0010.7010.629–0.773<0.0010.031
Predicted pVO2 (%)0.9180.860–0.975<0.0010.7010.628–0.774<0.001<0.001
VE/VCO2 slope0.8940.803–0.986<0.0010.6880.615–0.761<0.001<0.001
pVO2, mL/kg/min at GET0.6480.464–0.8320.0960.6350.451–0.8200.1400.594
O2 pulse, mL/kg/beat0.8160.669–0.9620.0010.6160.537–0.6950.0050.023
Circulatory power, mmHg.ml/kg/min0.7880.642–0.9350.0030.7130.646–0.780<0.0010.444
Ventilatory power, mmHg0.7820.597–0.9670.0040.7110.641–0.780<0.0010.504
COP0.6260.482–0.7700.0950.7040.560–0.8480.0190.372
PetCO2 at rest, mmHg0.6060.390–0.8220.2750.6540.580–0.728<0.0010.694
PetCO2 at GET, mmHg0.7840.638–0.9300.0040.7190.644–0.794<0.0010.461
CPET: Cardiopulmonary exercise testing; pVO2: Peak oxygen consumption; VE/VCO2 slope: Minute ventilation–carbon dioxide production relationship; GET: Gas exchange threshold; COP: Cardiorespiratory optimal point; PetCO2: Partial pressure of end-tidal carbon dioxide.
Table 5. Evaluation of traditional and alternative thresholds cut-off values of the composite endpoint.
Table 5. Evaluation of traditional and alternative thresholds cut-off values of the composite endpoint.
Female (n = 95)Male (n = 363)
Exercise Testing ParametersSpecificitySensitivityYouden (J) IndexSpecificitySensitivityYouden (J) Index
pVO2 ≤ 12 mL/kg/min *94%40%0.3491%21%0.12
pVO2 ≤ 14 mL/kg/min80%80%0.6082%47%0.29
pVO2 ≤ 15 mL/kg/min67%80%0.4779%57%0.36
VE/VCO2 slope > 3575%90%0.6566%57%0.23
VE/VCO2 slope > 3268%90%0.5857%78%0.35
Percent of predicted pVO2 ≤ 50%89%60%0.4978%48%0.26
Percent of predicted pVO2 ≤ 55%86%90%0.7669%60%0.29
Percent of predicted pVO2 ≤ 58%81%90%0.7163%69%0.32
* pVO2 ≤ 12 mL/kg/min (≤14 if the patient is intolerant to β-blockers). The highest Youden index (J) of each CPET variable is highlighted in bold. pVO2: Peak O2 consumption; VE/VCO2 slope: Minute ventilation–CO2 production relationship.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Garcia Brás, P.; Gonçalves, A.V.; Reis, J.F.; Moreira, R.I.; Pereira-da-Silva, T.; Rio, P.; Timóteo, A.T.; Silva, S.; Soares, R.M.; Ferreira, R.C. Cardiopulmonary Exercise Testing in Patients with Heart Failure: Impact of Gender in Predictive Value for Heart Transplantation Listing. Life 2023, 13, 1985. https://doi.org/10.3390/life13101985

AMA Style

Garcia Brás P, Gonçalves AV, Reis JF, Moreira RI, Pereira-da-Silva T, Rio P, Timóteo AT, Silva S, Soares RM, Ferreira RC. Cardiopulmonary Exercise Testing in Patients with Heart Failure: Impact of Gender in Predictive Value for Heart Transplantation Listing. Life. 2023; 13(10):1985. https://doi.org/10.3390/life13101985

Chicago/Turabian Style

Garcia Brás, Pedro, António Valentim Gonçalves, João Ferreira Reis, Rita Ilhão Moreira, Tiago Pereira-da-Silva, Pedro Rio, Ana Teresa Timóteo, Sofia Silva, Rui M. Soares, and Rui Cruz Ferreira. 2023. "Cardiopulmonary Exercise Testing in Patients with Heart Failure: Impact of Gender in Predictive Value for Heart Transplantation Listing" Life 13, no. 10: 1985. https://doi.org/10.3390/life13101985

APA Style

Garcia Brás, P., Gonçalves, A. V., Reis, J. F., Moreira, R. I., Pereira-da-Silva, T., Rio, P., Timóteo, A. T., Silva, S., Soares, R. M., & Ferreira, R. C. (2023). Cardiopulmonary Exercise Testing in Patients with Heart Failure: Impact of Gender in Predictive Value for Heart Transplantation Listing. Life, 13(10), 1985. https://doi.org/10.3390/life13101985

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop