Abstract
Background Normative changes in right ventricular (RV) structure and function have not been characterized in the context of treatment-associated functional recovery (RVFnRec). The aim of this study is to assess the clinical relevance of a proposed RVFnRec definition.
Methods We evaluated 63 incident patients with PAH by right heart catheterization and cardiac MRI (CMR) at diagnosis and CMR and invasive cardiopulmonary exercise (CPET) following treatment (∼11 months). Sex, age, race/ethnicity matched healthy control subjects (n=62) with one-time CMR and non-invasive CPET were recruited from the PVDOMICS project. We examined therapeutic CMR changes relative to the evidence-based peak oxygen consumption (VO2peak)>15mL/kg/min to define RVFnRec by receiver operating curve analysis. Afterload was measured in the as mean pulmonary artery pressure, resistance, compliance, and elastance.
Results A drop in RV end-diastolic volume of -15 mL best defined RVFnRec (AUC 0.87, P=0.0001) and neared upper 95% CI RVEDV of controls. 22/63 (35%) of subjects met this cutoff which was reinforced by freedom from clinical worsening, RVFnRec 1/21 (5%) versus no RVFnRec 17/42, 40%, (log rank P=0.006). A therapy-associated increase of 0.8 mL/mmHg in compliance had the best predictive value of RVFnRec (AUC 0.76, CI 0.64-0.88, P=0.001). RVFnRec subjects had greater increases in stroke volume, and cardiac output at exercise.
Conclusions RVFnRec defined by RVEDV therapeutic decrease of -15mL predicts exercise capacity, freedom from clinical worsening, and nears normalization. A therapeutic improvement of compliance is superior to other measures of afterload in predicting RVFnRec. RVFnRec is also associated with increased RV output reserve at exercise.
What is new?Right ventricular functional recovery (RVFnRec) represents a novel endpoint of therapeutic success in PAH. We define RVFnRec as treatment associated normative RV changes related to function (peak oxygen consumption). Normative RV imaging changes are compared to a well phenotyped age, sex, and race/ethnicity matched healthy control cohort from the PVDOMICS project. Previous studies have focused on RV ejection fraction improvements. However, we show that changes in RVEDV are perhaps more important in that improvements in LV function also occur. Lastly, RVFnRec is best predicted by improvements in pulmonary artery compliance versus pulmonary vascular resistance, a more often cited metric of RV afterload.
What are the clinical implications?RVFnRec represents a potential non-invasive assessment of clinical improvement and therapeutic response. Clinicians with access to cardiac MRI can obtain a limited scan (i.e., ventricular volumes) before and after treatment. Future study should examine echocardiographic correlates of RVFnRec.
Introduction
Adaptation of the right ventricle (RV) to high afterload is the primary determinant of a patient’s symptoms and survival in PAH, making it a focus of research by expert working groups1,2. As afterload increases in PAH, the RV remodels with hypertrophy resulting in adaptive maintenance of cardiac output. Eventually, maladaptive processes occur such as myocardial interstitial fibrosis, detrimental myocyte intracellular and molecular changes3,4, and alterations in coronary flow5 leading to impaired systolic and diastolic function. When RV afterload is acutely reversed, for example with lung transplantation, there is often complete reversibility of these RV maladaptive changes6. In the pre-combination therapy era, the magnitude of afterload lowering was not large enough to appreciate substantial changes in RV structure7. However, combination therapy results in significant improvements in afterload8 making RV “reverse remodeling” a realistic therapeutic goal.
Reverse remodeling of the left ventricle (LV) is characterized by “normative” (returning toward normal structure and function) changes. But, also the left ventricular, cellular, and molecular processes associated with reverse remodeling have been well described9. Relative to the LV, characterization of RV reverse remodeling is in its infancy. Current definitions7 have relied upon the relationship to freedom from clinical worsening10 but moreover, they lack a substantive pathophysiological link of the consequences of pulmonary vascular load on RV function. In the absence of longitudinal RV cellular and molecular data, a formal definition of RV reverse remodeling is premature. However, RV functional recovery (RVFnRec) can be characterized by therapeutically associated “normative” imaging changes paired with peak exercise data (VO2peak) indicative of the RV stress response. To date, studies have not included a well-matched control group, so characterization of “normative” change is unknown. In addition, most studies thus far have used echocardiography to describe RVFnRec, whereas cardiac MRI (CMR) may be more sensitive to treatment-related changes11.
It is generally accepted that a large drop in afterload is required for RV reverse remodeling (or RVFnRec), but the most affected parameter of afterload is unknown7. Both steady (i.e., pulmonary vascular resistance, PVR) and pulsatile (i.e., pulmonary arterial compliance, Ca) components are important contributors to afterload in the pulmonary circulation12. Effective pulmonary arterial elastance (Ea), a composite measurement of total PVR and Ca, is also proposed as a more comprehensive way of representing afterload13. Lastly, parameters obtained during exercise such as alpha (α) distensibility, which describes the curvilinearity of the PA pressure-flow relationship, may give us insight into pulmonary vascular reserve/recruitment and the state of the distal vasculature14. No study has yet investigated these parameters longitudinally from PAH diagnosis to determine their impact on reverse remodeling.
In the present study, we aim to systematically define RVFnRec by CMR in a group of treatment naïve PAH patients relative to a VO2peak threshold and clinical outcome. We evaluate the degree of normative RV therapeutic changes relative to matched healthy controls. We further assess the best parameter of RV afterload for predicting RVFnRec. Lastly, we examine if afterload and RV output reserve at exercise add to our understanding of RVFnRec.
Methods
Subjects
Sixty-three incident treatment naive patients with World Symposium of PH (WSPH) Group I PAH gave an informed consent to participate in the study, which was approved by the institutional review board at the University of Arizona (IRB #1100000621). The subjects enrolled as part of a prospective protocolized study which is open to all incident patients (See supplemental methods). The diagnosis of PAH was established by dedicated PAH providers with invasive confirmation by right heart catheterization based on updated guidelines15. Patients were included if they had right heart catheterization and CMR at baseline (pre-treatment) and invasive cardiopulmonary exercise testing (iCPET) and CMR at > 6 months of follow-up. All patients were placed on PAH goal-directed therapy per current guidelines at the time. See supplemental methods for details on therapeutic strategy. Subjects were categorized by therapeutic strategy based on the European Respiratory Society/European Society of Cardiology (ERS/ESC) guidelines of 2015 at which goal-directed mono or sequential therapy was replaced by up-front combination therapy16. We selected a cohort of patients all placed on parenteral treprostinil (± oral therapy) to optimize the possibility of seeing large changes in afterload (and thus, RVFnRec)17. Ninety-six healthy control subjects were recruited for the NHLBI Redefining Pulmonary Hypertension through Pulmonary Vascular Disease Phenomics (PVDOMICS)18. These subjects underwent extensive testing including CMR and non-invasive CPET at a single time point. Sixty-two age, sex, and race/ethnicity matched controls from across all clinical sites were used for this analysis with steering committee approval. The United States Registry to Evaluate Early and Long-Term PAH Disease Management (REVEAL) score 2.0 score19 was calculated and presented at baseline and follow-up for the PAH subjects.
Cardiac MRI
All image analysis was performed using CVI42, Circle Cardiovascular Imaging® by TA (at Arizona) and DK (at the PVDOMICS Core, the Cleveland Clinic), expert cardiac MRI cardiologists with >10 years’ experience in CMR. Fifty subjects with overlap (between PVDOMICS and the UA registry) were analyzed by both TA and DK and used to assess inter-reader variability in RV volumes (see supplemental methods). Cardiac MRI imaging was performed as described previously20. Briefly, standard volumetric measurements were made from short axis cine projections for CMR. Contour smoothing was done to include trabeculations in end-systole (ESV) and diastole (EDV). Stroke volume (SV) was calculated as EDV – ESV and RVEF as [(SV/EDV) x 100]. Additional methods are available in the online supplement.
Right Heart Catheterization and Invasive Cardiopulmonary Exercise Testing
After study inclusion, a pulmonary artery catheter was advanced via the antecubital vein for measurements of pulmonary artery pressure (systolic, mean, and diastolic PAP), RV pressure, right atrial pressure (RAP), wedged PAP (PCWP) and cardiac output (direct Fick). PVR is calculated as mPAP-PCWP divided by direct Fick cardiac output (C.O.) and expressed as WU (mmHg*min*L-1) or mmHg*sec*mL-1 where indicated. Ca was calculated as stroke volume (SV) from direct Fick C.O. divided by PA pulse pressure and expressed as mL*mmHg-1. Effective pulmonary arterial elastance (Ea) was calculated as sRVP/SV 21 expressed as mmHg*mL-1. Distensibility (α), expressed as % change in vessel diameter/mmHg distending pressure, was calculated by methods previously described 14,22. All stages of a single subject’s mPAP, PCWP, and C.O. were used to calculate a single α value. RC time, the product of PVR and Ca, was expressed as seconds.
Our comprehensive resting and exercise catheterization protocol has been previously published23. Briefly, after obtaining supine resting measurements, the patient was placed in full upright position with an electronic fluoroscopy chair. Fluoroscopy was used to re-zero at left atrial level. A cycle ergometer was positioned below the patient. The patient then proceeded with exercise at predetermined workload based on their level of dyspnea 23 for steady-state two-minute stages until respiratory exchange ratio (RER) ∼ 1.1. PAP, PCWP, RAP, arterial and venous O2 content (CaO2, CvO2) were obtained during the last 30 seconds of each stage. Metabolic cart analysis (Vyaire Medical™, Mettawa, IL) was used for simultaneous collection gas exchange and lung volume. Hemodynamics presented were averaged over 3 respiratory cycles in accordance with current guidelines at rest and exercise24. The control patients underwent an identical protocol in the non-invasive exercise lab (without a catheter in place).
Definition of RV Functional Recovery
Since exercise VO2peak is physiologically meaningful (dependent on RV function) and prognostic, the guideline25 and evidence (clinical worsening) confirmed26 cutoff of >15 mL/kg/min was chosen to identify potential CMR cutoffs indicative of RVFnRec. Exercise data is presented only for the follow-up visit as many patients (N=55 without iCPET, N=8 with iCPET) were not exercised at baseline due to safety considerations. We examined therapeutically associated change in RV size (RVESV and EDV, % change RV volumes) and function (RVEF) relative to this VO2peak cutoff at follow-up. The Youden index was used to define the best cutoffs in RV volumes (absolute and relative change from baseline) and RVEF (change from baseline). Receiver operating characteristic (ROC) curves are used to display optimal cut-offs. Resulting RVFnRec groups were then examined relative to freedom from clinical worsening. Clinical worsening was defined as death, transplant, or hospitalization. Freedom from clinical worsening was evaluated using Kaplan-Meier Plots and the log-rank test. A similar procedure (Youden index and ROC analysis) was used to define the best cutoffs for mPAP, Ca, PVR, Ea, and α at follow-up and change from baseline to follow-up for predicting RVFnRec.
Statistical Analysis
Continuous data are expressed as mean ± standard deviation or median [25,75 percentile]. Categorical data are expressed as counts and percentages. The number of PAH subjects was chosen based primarily on data availability given strict inclusion/exclusion criteria. However, based on previously published data in prostacyclin treated patients, we would expect 100% power to detect a mean difference of 7±5.4 WU PVR difference between RVFnRec and no RVFnRec groups (R.B. personal communication, 9/25/22 based on ref.17. Kruskal-Wallis test for continuous and ordinal variables or Fisher exact test for categorical variables were used to test baseline differences. Normality was assessed with the Shapiro-Wilk test. Longitudinal differences were evaluated with Wilcoxon signed-rank (within group) or repeated measures analysis of variance (ANOVA)(between group). Variables were log-transformed in cases of non-normal distribution. Linear or non-linear regression was conducted where appropriate after testing assumptions. Binary logistic regression was done to assess the relationship of treatment type with RVFnRec. Statistical analyses were performed using SPSS software (version 28.0, IBM, Armonk, NY) and SigmaPlot (version 14.5, Systat©, San Jose, CA). Statistical tests were 2-sided, and a p-value <0.05 was considered statistically significant. P values and 95% CIs presented in this report have not been adjusted for multiplicity, and therefore, inferences drawn from these statistics may not be reproducible.
Results
Patient Characteristics
Sixty-three treatment naïve PAH and 62 control subjects were prospectively enrolled. Subjects were predominantly female and PAH subjects were majority idiopathic (Table 1). At enrollment, the subjects had very advanced PAH, characterized by low cardiac index and significantly elevated afterload. CMR demonstrated nearly twice the RVEDVI and depressed RVEF relative to controls. Therapeutic improvements in RVEF were accompanied by improvements in 6-minute walk, PVR, pulmonary compliance, cardiac output, and BNP (Table 1).
Values are mean ± SD, median [P25,P75], or N (%). BMI, body mass index; CHD, congenital heart disease; CTD, connective tissue disease; HIV, human immunodeficiency virus; iPAH, idiopathic pulmonary arterial hypertension; PoPH, portopulmonary hypertension; BNP, brain natriuretic peptide; MRI, magnetic resonance imaging; RVESVI, right ventricular end-systolic volume index; RVEDVI, right ventricular end-diastolic volume index; RVEF, right ventricular ejection fraction; PAP, pulmonary artery pressure; PCWP, pulmonary capillary wedge pressure; RAP, right atrial pressure; PVR, pulmonary vascular resistance. P values calculated as follows: Kruskal-Wallis test for continuous and ordinal variables or Wilcoxon signed-rank non-parametric test for related samples, Fisher exact test for categorical variables. *Significant difference at level <0.05 within PAH follow-up versus baseline; † P<0.05 control versus follow-up 1 PAH.
Right Ventricular Functional Recovery
As demonstrated in Figure 1A-B, VO2peak was correlated with RVEDV and RVEF. 22/63 (35%) subjects had VO2peak>15mL/kg/min with 20 (91% of the 22) having a change in RVEDV and RVEF meeting the cutoff defined by the highest Youden index (0.74 at - 15mL and 0.51 at 4% change, respectively) while 2 (9%) met only one of these cutoffs. Of the 41/63 (65%) of subjects with VO2peak≤15mL/kg/min, 22 (54% of the 41) were below both cutoffs, 15 (37%) had RVEF above the 4% change but below the RVEDV cutoff, whereas 4 (9%) met both the RVEF and RVEDV cutoffs. ROC analysis (Figure 1C) confirmed changes in RVEDV were superior to changes in RVEF in the prediction of VO2peak. Indexing RVEDV to BSA (RVEDVI) was not more predictive than RVEDV (Youden 0.65 at -10.4 mL/m2). Change in RVESV and relative change in RVESV and EDV were not more predictive than change in RVEDV (Table E1). All PAH subjects who achieved the high VO2peak cutoff demonstrated improvements in both RVESV and RVEDV, whereas some subjects below the VO2peak cutoff had improvements in RVESV alone (Figure E1).
A change in RV end-diastolic volume (RVEDV) (A) and RV ejection fraction (RVEF) (B) from baseline to follow-up both correlated with peak oxygen consumption, VO2peak, at follow-up. However, a change in RVEDV better discriminated high exercise capacity (VO2peak >15mL/kg/min) than did a change in RV ejection fraction (RVEF). Vertical solid lines indicate best Youden index of 0.72 at an RVEDV of -15mL or 0.51 for RVEF at +4%. Based on these cutoffs, the number of true positives (TP) and true (TN) and false (FN) negatives were similar in RVEDV and RVEF. However, the number of false positives (FP) was higher for RVEF. C) Receiver operating characteristic (ROC) analysis predicting VO2peak>15 mL/kg/min shows a change in RVEDV of -15mL is more predictive than RVEF +4mL as well. D) Kaplan-Meier plot of RVFnRec and no RVFnRec by clinical worsening shows higher freedom from clinical worsening in the RVFnRec group.
22/63 (35%) of subjects met criteria for RVFnRec using a drop of -15 mL in RVEDV. Using this definition of RVFnRec, the average length of follow up was 50±28 and 38±25 months for the RVFnRec and no RVFnRec groups, respectively. All-cause mortality was 1 (5%) and 10 (24%) and clinical worsening was 1 (5%) and 17 (41%) in the RVFnRec and no RVFnRec groups, respectively. As shown in Figure 1D, this definition of RVFnREc also resulted in a significant difference in time to clinical worsening, P=0.006. There were significantly more subjects in the RVFnRec group treated with up-front combination therapy versus the no RVFnRec group (50 versus 24%, P=0.01, Table 2 and Figure E2). Treatment approach differences were due in part to a difference between groups in time of enrollment (RVFnRec median enrollment date 4/2017 [3/2014-5/2021] versus no RVFnRec 5/2018 [11/2013-12/2021]). Up-front combination therapy was associated with odds-ratio of 10.4 (CI 1.9-56.6, P=0.007) of RVFnRec relative to goal-directed sequential therapy. There were significant improvements in 6-minute walk distance, REVEAL 2.0, RV volumes and EF in the RVFnRec group by follow-up compared to the no RVFnRec group (Table 2 and 3). Although both groups had an increase in LVEDV, only the RVFnRec group had an improvement in LVEF (Table 3).
Values are mean ± SD, median [P25,P75], or N (%). ADO, aldosterone; BMI, body mass index; CHD, congenital heart disease; CTD, connective tissue disease; HIV, human immunodeficiency virus; iPAH, idiopathic pulmonary arterial hypertension; PoPH, portopulmonary hypertension; REVEAL, United States Registry to Evaluate Early and Long-Term PAH Disease Management; ERS, European Respiratory Society; BNP, brain natriuretic peptide; MRI, magnetic resonance imaging; RVESVI, right ventricular end-systolic volume index; RVEDVI, right ventricular end-diastolic volume index; RVEF, right ventricular ejection fraction; PAP, pulmonary artery pressure; PCWP, pulmonary capillary wedge pressure; RAP, right atrial pressure; PVR, pulmonary vascular resistance; TRE, treprostinil. P values calculated as follows: Kruskal-Wallis test for continuous and ordinal variables or Wilcoxon signed-rank non-parametric test for related samples (within cohort) or repeated measures ANOVA (between cohort), Fisher exact test for categorical variables. * P<0.05 versus baseline within RVFnRec or no RVFnRec cohort; † P<0.05 between cohorts.
Values are mean ± SD, median [P25,P75], or N (%). RVEDVI, right ventricular end-diastolic volume index; RVESVI, right ventricular end-systolic volume index; RVEF, right ventricular ejection; RVMI, right ventricular mass index; RV M/V, right ventricular mass/volume; RV SV/ESV, right ventricular stroke volume/end-systolic volume; LVEDVI, left ventricular end-diastolic volume index; LVESVI, left ventricular end-systolic volume index; LVEF, left ventricular ejection-fraction; LVMI, left ventricular mass index; TR volume, tricuspid regurgitant volume; TR fraction, tricuspid regurgitant volume fraction (relative to stroke volume) fraction. P values calculated as follows: Wilcoxon signed-rank non-parametric test for related samples (within cohort) or repeated measures ANOVA (between cohort). * within cohort change from baseline to follow-up; † p value for between cohort change baseline to follow-up; ‡ p-value for between group relative change (as % baseline value) baseline to follow-up.
As shown in Figure 2A, only the RVFnRec group demonstrated a therapeutic change in RVEDV that nears “normalization” of RVEDVI at the upper 95% confidence interval (CI) of healthy controls (90 mL/m2). -15 mL decrease in RVEDV identifies most of the patients included in this confidence interval. Although most of both the RVFnRec(20/22, 91%) and no RVFnRec (32/41, 78%) groups met the SVI lower 95% CI of controls (30 mL/m2), only the RVFnRec group neared the upper and lower 95% CI of controls for RVEDVI (11/22, 50% versus 10/41, 24%, respectively) and RVEF (8/22, 35%) versus 2/41, 5%, respectively) at follow-up (Figure 2B). Despite near normalization of EDV and RVEF among the RVFnRec group, relative wall thickness (RV mass/volume) remained approximately twice that of controls, RVFnRec 0.43±0.13 gm/mL versus control 0.21±0.05 gm/mL (Figure 2C, Table 1 and 2). Table E2 compares RV volume and mass in PAH and control subjects with previously published control reference values.
A) The treatment related change in RV end-diastolic volume (RVEDV) relative to RVEDV index at follow-up in the functional recovery (RVFnRec) and no RVFnRec subjects. Only RVFnRec subjects achieve a therapeutic change in RVEDV that nears normalization. B) RVEF and RVEDVI at follow-up in the RVFnRec, no RVFnRec, and healthy controls. Solid vertical line indicates the upper 95% CI and the dark horizontal line the lower 95% CI for RVEDV and RVEF in the control cohort, respectively. C) Treatment related changes in short (first column) and long (second column) axis cardiac MRI images at Baseline and Follow-up. Substantial, improvements in RV volume and ejection fraction are seen in the RVFnRec group. However, relative wall thickness (mass/volume) remained almost twice the healthy control value in both the RVFnRec and no RVFnRec groups.
RV Afterload Parameters Associated with RV Functional Recovery
Patients with RVFnRec had larger absolute and relative (to baseline) drops in mPAP, and Ca, as well as relative changes in PVR, versus the no RVFnRec group by follow-up (Table 4). The absolute change in PVR and Ea and relative changes in Ea were not significant (Table 4). A greater proportion of RVFnRec than no RVFnRec subjects had moved to the “steep” portion of the PVR-Ca curve (Figure 3A) by follow-up. 9/63 (∼14%) of PAH subjects had near normalization of PVR (∼3WU) (Figure 3A). The accuracy of the change in the absolute value, change in relative value (as % baseline), and absolute value at follow-up were assessed as shown in Figure 3B and E3A-B. The value giving the best AUC 0.76 (95% CI 0.64-0.88, P<0.0001) was change in Ca at a cutoff of 0.8 mL/mmHg yielding a sensitivity of 85% and a specificity of 72% for RVFnRec. Changes in mPAP, PVR, and Ca were associated with changes in RVEDV (Figure 3 C-E) and RVEF (Figure E3C-E) with Ca being the most highly predictive.
Values are mean ± SD, median [P25,P75], or N (%). mPAP, mean pulmonary artery pressure; PVR, pulmonary vascular resistance; PA, pulmonary artery. P values calculated as follows: Wilcoxon signed-rank non-parametric test for related samples (within cohort) or repeated measures ANOVA (between cohort). * within cohort change from baseline to follow-up; † p value for between cohort mean change baseline to follow-up; ‡ p-value for between group relative change (as % baseline value) baseline to follow-up.
(A) non-linear relationship between pulmonary vascular resistance (PVR) and compliance (Ca) at baseline and follow-up. Most RVFnRec subjects had moved to the steep portion of the curve by follow-up. Solid line indicates an upper normal PVR of <3 WU. (B) receiver operating curve analysis of longitudinal changes between baseline and follow-up in mean pulmonary artery pressure (mPAP), PVR, Ca, and effective pulmonary elastance (Ea). Changes in Ca were the most closely related afterload parameter to changes in RV end-diastolic volume (RVEDV) (C-E). Vertical line indicates the highest (“best”) Youdin index cutoff for each afterload parameter.
The Relationship of RV Functional Recovery to Exercise Afterload and RV Output Reserve
The RVFnRec subjects had higher VO2peak predicted, mPAP, cardiac output, and stroke volume at exercise than did the no RVFnRec subjects at follow-up (Figure 4 and Table 5). The VO2peak was much lower than matched controls, however. Ventilatory efficiency index for carbon dioxide (Ve/VCO2) was lower (better) in the RVFnRec group. The mPAP-cardiac output relationship from rest to peak exercise, was similar in RVFnRec versus the no RVFnRec group (Figure 4B). α distensibility was also similar between RVFnRec and no RVFnRec groups (0.17±0.15 versus 0.19±0.3 %/mmHg, P=NS, respectively) and much lower than predicted normal values of 1.4 (range 0.77-2.34) %/mmHg14(Figure 4D). Using exercise afterload parameters did not improve the accuracy of resting values at follow-up to predict RVFnRec (Figure E4).
Values are mean ± SD, median [P25,P75]. P values calculated from the Kruskal-Wallis test. RER, respiratory exchange ratio; VO2, peak oxygen consumption; RAP, right atrial pressure; mPAP, mean pulmonary artery pressure; PCWP, pulmonary capillary wedge pressure; PAO2, pulmonary arterial oxygen saturation; PVR, pulmonary vascular resistance; Ve/VCO2, minute ventilation/exhaled carbon dioxide; PetCO2, end-tidal CO2; O2pulse, oxygen pulse (VO2/heart rate), CaO2, arterial oxygen content; CvO2, pulmonary arterial oxygen content; Ca-vO2, arterio-venous oxygen content difference; α distensibility, distensibility of the pulmonary arterial system (calculated at exercise). *P<0.05 for RVFnRec versus No RVFnRec; †=P<0.05 for Health Control versus RVFnRec: ‡P<0.05 for Healthy Control versus No RVFnRec.
Peak predicted oxygen consumption was higher in RVFnRec than no RVFnRec but both were well below normal matched controls (A). The linear pressure-flow relationship of RVFnRec subjects and no RVFnRec subjects was similar (B). The arrowheads represent peak and the tails resting mPAP/cardiac output. The thick arrow represents the median vector in each group. Despite similar mPAP/cardiac output, the cardiac output at peak was higher (B) and the change in stroke volume index (SVI) (C) was higher in RVFnRec than that of no RVFnRec subjects. Pulmonary vascular α distensibility was similar between groups (D). Taken together (B-D) these findings indicate similar RV afterload at exercise but improved RV output reserve in the RVFnRec cohort.
Discussion
This study is the first to demonstrate in treatment naïve PAH patients that 1) a -15 mL decrease in RVEDV is a useful CMR definition of RVFnRec associated with exercise capacity (VO2peak) and freedom from clinical worsening, 2) an increase in resting Ca of 0.8 mL/mmHg is superior to other afterload parameters in predicting RVFnRec, and 3) RVFnRec is associated with high exercise RV output reserve.
Although therapeutic improvements in RV size and function add independently prognostic information to well validated risk scores (such as REVEAL 2.0 and ERS risk score) in PAH10,27, the lack of a consensus definition of RV recovery hampers its use as a therapeutic goal and a clinical trial endpoint. As normalization of RV size and function would be an optimal goal, defining the limits of “normal” is challenging given reference populations have varied by age, sex, race/ethnicity, body mass index, and acquisition/reading methodology28-30. We have attempted to mitigate these challenges by enrolling a control group matched for age, sex, and race/ethnicity. Although many of our subjects met the CMR 95% confidence interval in the controls, normalization appears to be an unrealistic definition of RVFnRec. However, complete normalization while a laudable goal may be unnecessary for meaningful RV functional recovery. Unlike definitions of RV reverse remodeling which rely on the absence of clinical worsening alone10, we thought it important to examine RVFnRec by CMR cutoffs relating to clinical improvement and exercise capacity. The VO2peak chosen at 15 mL/kg/min is dependent on cardiac output (RV functional reserve) in PAH26 and is derived from guidelines25 and evidence26. We confirmed the findings of others by showing changes in RVEDV are related to freedom from clinical worsening in the RVFnRec group.
There is consensus that RVEF is consistently related to mortality 1 and therapeutic improvements in RVEF are related to a drop in afterload 10,31. Studies enrolling subjects on oral monotherapy or sequential therapy result in minimal afterload reduction. In these studies, improvement in RVEF occurs by reductions in RVESV but with minimal change in RVEDV27,31-33. Early therapeutic changes in RVESV result from improvements in afterload and RV diastolic function20. But RV work efficiency does not improve (vis a vis improved RV-PA coupling) until there is a drop in RVEDV with larger afterload reduction (∼7WU or >45% drop in PVR)34. In fact, we demonstrate that a drop in both RVESV and RVEDV are highly associated with an acceptable (VO2peak >15mL/kg/min) exercise capacity. We attempted to optimize the probability of large changes in afterload by employing an aggressive therapeutic strategy based on parenteral prostacyclin 7,17. A drop in RVEDV may signify a physiological shift from both heterometric (Frank-Starling) and homeometric (contractility) adaptation to maintain stroke volume to a more energetically optimized condition of homeometric adaptation alone35. Although improvements in RVEF from decreased RVESV represent a formidable short-term therapeutic goal, the increased myocardial efficiency and reduced wall stress accompanied by decreasing RVEDV should represent the ultimate long-term goal 35.
Clinical observation indicates that biventricular maladaptive remodeling in PH is reversible with rapid hemodynamic normalization such as lung transplantation6 and pulmonary thomboendarterectomy (PTE) for chronic thromboembolic PH (CTEPH)36. Even RV fibrosis appears reversible in experimental PAH when afterload is gradually but completely normalized37. We have shown for the first time that biventricular normative changes occur particularly in up-front prostacyclin with oral combination therapy resulting in a significant and gradual drop in afterload. Only patients with RVFnRec demonstrated large reductions in RV volume and improvements in RV and LVEF. Therefore, a reduction in RVEDV may be requisite to an improvement in the LV maladaptive changes, such as reduced LVEF36, seen with PAH and CTEPH. Although RVEDV shows near “normalization” relative to controls, normalization of RVEF and mass/volume was less common. RV relative wall thickness (mass/volume) was maintained more than double that of controls likely reflecting reduced but persistently high load. In this context, the present definition of RVFnRec is near the prognostically significant cut-point of 0.47 gm/mL38. The combination of RVFnRec and a high RV mass/volume appears physiologically significant during exercise where these subjects can maintain higher cardiac output and stroke volume (RV output reserve) as compared to those without RVFnRec despite similar PAP39,40.
As recently reviewed7, knowledge gaps exist on which afterload parameter predicts RVFnRec. Although the pulmonary impedance spectrum is the gold standard of afterload measurement, it requires simultaneous pressure and flow measurement in the frequency domain making it difficult to measure and interpret12. We have therefore chosen a comprehensive list of frequently measured afterload parameters in clinical medicine. Consistent with previously published data, we found that large changes in PVR (∼>50%) are necessary for RVFnRec 10. However, pulmonary compliance was significantly more predictive of RVFnRec than measures of steady load such as PVR and the composite measure of load, Ea. Studies have demonstrated that compliance has a more consistent relationship with outcome 41-43 and RV function 44,45 than does PVR. Also, experimental 46 and clinical 47 evidence indicates that compliance plays a significant role in maintaining optimal RV-PA coupling even at high PVR. Nevertheless, we found only a mild correlation between all standard load parameters and RV function. This finding may indicate that some afterload parameters not measured, like wave reflection 48, or load-independent factors, like possible direct prostacyclin-mediated effects on the RV such as reduced myocardial inflammation 49 or coronary microcirculatory effects 50, may play a role in RVFnRec. It is possible that RVFnRec may represent a more comprehensive therapeutic target, encompassing both load dependent and independent effects, than targeting load (i.e., PVR) alone.
We found that resting therapeutic improvements in afterload dissipate quickly during exercise as high PAP relative to cardiac output is seen even in subjects with RVFnRec. Likely, the drop in load seen at rest is negated at exercise as less affected pulmonary vessels meet maximal recruitment and distension as evidenced by very low α distensibility and high mPAP-C.O. slope in both groups. This observation suggests that pulmonary vascular remodeling is still likely advanced in RVFnRec despite such significant improvements in resting measurements of afterload. The demands on the RV by high afterload at exercise may explain why RVFnRec subjects maintain nearly twice the relative RV wall thickness (mass/volume) to that of controls despite near normalization of ventricular volume.
Our study has several limitations. It is a single center PAH cohort and therefore the results require validation amongst a larger population and preferably multiple centers. Validation is particularly important given that our ROC cutoff was assessed only at two time points (baseline and follow-up) and therefore does not account for intrinsic variability during the course of disease. For example, the RVEDV and Ca cutoffs may only be applicable to subjects with advanced PAH. Also, our study assessed VO2peak only at follow-up and not the change in VO2 from diagnosis. However, the change in 6-minute walk distance from baseline to follow-up confirms a substantial difference in functional capacity between RVfnRec groups. Also, the feasibility of MRI is difficult at some centers where echocardiography is preferable. Future studies should validate our findings against echocardiographic metrics such as those recently proposed 7.
Conclusions
Our data suggest that defining RVFnRec as a therapeutic normative -15 mL drop in RVEDV appears to be clinically useful. Among standard measures of afterload, an increase in Ca of 0.8 mL/mmHg by follow-up best predicts RVFnRec although factors outside afterload alone may be at play. We have shown that our definition of RVFnRec corresponds with near normalization of RV volume, but relative wall thickness remains high. The high relative wall thickness may be needed for output reserve in the face of high afterload at exercise.
Data Availability
All data is available upon reasonable request from the contact author.
Author Contributions
F.P.R., R.B., D.H.K., T.A., W.W.T.: study design; F.P.R., M.I., S.K.: patient recruitment, care and follow-up; MRI core interpretation; A.K., M.M.P., T.A., D.H.K.: rest and exercise hemodynamic core interpretation; F.P.R., M.I., S.K., J.D.W., W.W.T. .: data collection, maintenance, and analysis; R.J.B., F.P.R., R.P.F., P.M.H., A.R.H, E.M.H., J.A.L., E.B.: statistical analysis; F.P.R, A.B.L.:. drafted the original manuscript; F.P.R., R.N., R.R.V., R.B., R.J.B., A.B.L., D.H.K, T.A., critical revision of the manuscript for important intellectual content; F.P.R., M.I., S.K., G.J.B., N.S.H., R.P.F., P.M.H., A.R.H., E.M.H., J.G.G., J.A.L., E.B., S.C.E : principal investigator, had access to the study data and takes full responsibility for the integrity and accuracy of the data.
Acknowledgements
We would like to thank Andrew Swift (Department of Infection, Immunity and Cardiovascular Disease, University of Sheffield, Western Bank, Sheffield, UK) for sharing his CMR predicted reference equations with us. We would like to graciously acknowledge the assistance of the UAHD Biorepository at the University of Arizona. We would also like to acknowledge the many PH subjects of the UA registry/biorepository, the PVDOMICS study, and coordinators who have dedicated significant effort in the academic pursuit to further our field.
Footnotes
Sources of funding: The PVDOMICS study received grants U01 HL125218 (PI: E.B. Rosenzweig), U01 HL125205 (PI: R.P. Frantz), U01 HL125212 (PI: A.R. Hemnes), U01 HL125208 (PI: F.P. Rischard), U01 HL125175 (PI: P.M.Hassoun), U01 HL125215 (PI: J.A. Leopold), and U01 HL125177 (PI: G.J. Beck) and the Pulmonary Hypertension Association. The PAH cohort is a continuously enrolling (since 1/21) cohort at UA of treatment naïve patients. This registry/biorepository is funded directly by the UA.
Author Disclosures: Please see attached ICMJE disclosure forms.
Supplement: This article has an online data supplement, which is accessible from this issue’s table of content online.
Abbreviations
- ANOVA
- analysis of variance
- BNP
- brain natriuretic peptide
- Ca
- pulmonary vascular compliance
- CI
- confidence interval
- CMR
- cardiac MRI
- CTEPH
- chronic thromboembolic pulmonary hypertension
- EDV
- end-diastolic volume
- EF
- ejection fraction
- ERS/ESC
- European Respiratory Society/European Society of Cardiology
- ESV
- end-systolic volume
- Ea
- effective pulmonary arterial elastance
- LV
- Left ventricle
- PA
- pulmonary artery
- PAH
- pulmonary arterial hypertension
- PAP
- pulmonary artery pressure
- PCWP
- pulmonary capillary wedge pressure
- PVR
- pulmonary vascular resistance
- RAP
- right atrial pressure
- REVEAL
- United States Registry to Evaluate Early and Long-Term PAH Disease Management
- RV
- right ventricle
- RVFnRec
- RV functional recovery
- sRVP
- systolic RV pressure
- SV
- stroke volume
- PTE
- thromboendarterectomy
- WSPH
- World Symposium of Pulmonary Hypertension
- VO2peak
- peak oxygen consumption
References
