Advertisement
Original Article|Articles in Press

Association between the pulmonary artery pulsatility index and prognosis in pulmonary arterial hypertension: A multicenter study

Open AccessPublished:April 25, 2023DOI:https://doi.org/10.1016/j.cjco.2023.04.005

      Abstract

      Background

      Risk stratification is fundamental in the management of pulmonary arterial hypertension (PAH). Pulmonary artery pulsatility index (PAPi), defined as pulmonary arterial pulse pressure divided by right atrial pressure (RAP), is a hemodynamic index shown to predict acute right ventricular (RV) dysfunction in several settings. Our objective was to test the prognostic utility of PAPi in a diverse multicentre cohort of patients with PAH.

      Methods

      A multi-centre retrospective cohort study of consecutive adult patients with a new diagnosis of pulmonary arterial hypertension (PAH) on right heart catheterization (RHC) between January 2016 and December 2020 was undertaken across four major centres in Canada. Hemodynamic data, clinical data, and outcomes were collected. The association between PAPi and other hemodynamic variables and mortality was assessed by receiver-operating characteristic (ROC) curves and Cox proportional hazards modeling.

      Results

      We identified 590 patients with a mean age 61.4±15.5 years with 66.3% being female. A low PAPi (defined as < 5.3) was associated with higher mortality at 1 year 10.2% vs. 5.2% (P=0.02). In a multivariable model including age, sex, body mass index and functional class, a low PAPi was associated with mortality at 1 year (AUC of 0.64, 95% CI 0.55-0.74). However, high RAP (> 8 mmHg) was similarly predictive of mortality with an AUC of 0.65.

      Conclusion

      PAPi was associated with mortality in a large incident PAH cohort. However, it’s discriminative value was not higher than RAP alone.

      Introduction

      Pulmonary arterial hypertension (PAH) is a progressive chronic disease characterized by pulmonary vascular remodeling, eventually culminating in right ventricular (RV) failure and death if untreated(
      • Parikh V.
      • Bhardwaj A.
      • Nair A.
      Pharmacotherapy for pulmonary arterial hypertension.
      ). Despite advances in treatment, PAH continues to be associated with high morbidity and mortality(
      • Boucly A.
      • Weatherald J.
      • Savale L.
      • de Groote P.
      • Cottin V.
      • Prevot G.
      • et al.
      External validation of a refined four-stratum risk assessment score from the French pulmonary hypertension registry.
      ,
      • Zelt J.G.E.
      • Sugarman J.
      • Weatherald J.
      • Partridge A.C.R.
      • Liang J.C.
      • Swiston J.
      • et al.
      Mortality trends in pulmonary arterial hypertension in Canada: a temporal analysis of survival per ESC/ERS guideline era.
      ). The most recent guidelines for the diagnosis and treatment of PAH propose a risk stratification tool based on multiple variables, including hemodynamic parameters such as mean right atrial pressure (mRAP), cardiac index (CI), and mixed venous oxygen saturation (SvO2)(
      • Galie N.
      • Channick R.N.
      • Frantz R.P.
      • Grunig E.
      • Jing Z.C.
      • Moiseeva O.
      • et al.
      Risk stratification and medical therapy of pulmonary arterial hypertension.
      ,
      • Galie N.
      • Humbert M.
      • Vachiery J.L.
      • Gibbs S.
      • Lang I.
      • Torbicki A.
      • et al.
      2015 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension: The Joint Task Force for the Diagnosis and Treatment of Pulmonary Hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS): Endorsed by: Association for European Paediatric and Congenital Cardiology (AEPC), International Society for Heart and Lung Transplantation (ISHLT).
      ). A comprehensive assessment is recommended as no single variable provides sufficient diagnostic and prognostic information. Additional markers may be beneficial for further refining assessment of prognosis, which is fundamental as it guides therapeutic decisions, including listing for transplantation.
      RV function and RV-pulmonary arterial (PA) coupling are strongly associated with prognosis in PAH(6). RV-PA coupling, defined by a ratio of end-systolic to arterial elastances, has been shown to detect developing RV failure(
      • Tello K.
      • Dalmer A.
      • Axmann J.
      • Vanderpool R.
      • Ghofrani H.A.
      • Naeije R.
      • et al.
      Reserve of Right Ventricular-Arterial Coupling in the Setting of Chronic Overload.
      ). Additionally, this parameter has been shown to be an independent predictor of survival in patients with PAH(8). The pulmonary artery pulsatility index (PAPi), defined as the ratio of the PA pulse pressure to the mRAP, has been proposed as an easily derivable marker of RV-PA coupling. Since its initial description in patients with acute inferior wall myocardial infarction, it has not yet received widespread evaluation in a large multicentre study of patients with PAH(9). The objectives of this study were to investigate the association between PAPi and mortality in patients with PAH, and to determine whether it provides superior prognostic value compared to other hemodynamic variables from right heart catheterization (RHC).

      Methods

      Study design and population

      This was a retrospective, multi-center cohort study across four pulmonary hypertension expert referral centres participating in the Canadian Pulmonary Hypertension Registry (CPHR) in Canada (Vancouver, Calgary, Ottawa, Hamilton). The study was coordinated through the University of British Columbia (H21-03574) and institutional ethics board approval was obtained at each participating center. Eligible patients were identified through a structured chart review by study authors (NB, AA, JS, JW, NH, LM). All incident patients with a clinical diagnosis of PAH confirmed on RHC between January 1, 2016 to Dec 31, 2020 were included. In this study, the hemodynamic criteria for PAH proposed by the 6th World Symposium on PH were used (pulmonary arterial pressure (mPAP) > 20 mmHg, a pulmonary arterial wedge pressure ≤ 15 mmHg and a pulmonary vascular resistance > 3 WU). Other causes of PH, including patients with group 2 (PH due to left heart disease), group 3 (PH due to lung disease and/or hypoxia), group 4 (PH due to pulmonary artery obstructions), and group 5 (PH with unclear and/or multifactorial mechanisms) were excluded. Patients with missing hemodynamic data were excluded from the analysis. Patient demographics including age, sex, body mass index (BMI), and etiology of PH as determined by the treating specialist were obtained from electronic medical records or clinical files. Hemodynamic variables were extracted from baseline RHCs. The primary outcome was all-cause mortality at 1 year following the index RHC, as determined from clinic records. We chose to analyze survival at 1 year as longer-term survival was expected to be more strongly influenced by response to PAH therapy as opposed to baseline hemodynamic parameters.

      Hemodynamic variables

      RHC were performed by experienced operators according to routine clinical care. During RHC, the zero line was set at the mid-thoracic height. Measurements were obtained from hemodynamic tracing according to current guidelines. Variables measured during the RHC included, but were not limited to, the mRAP; systolic (sPAP), diastolic (dPAP), and mean pulmonary artery pressures (mPAP); pulmonary arterial wedge pressure (PAWP); and SvO2. Cardiac output and index were determined by either thermodilution or Fick method, at the discretion of the operator. However, the thermodilution method was preferred in the absence of shunting. PAPi was calculated as pulmonary artery pulse pressure (sPAP – dPAP) divided by mRAP. PAPi was compared to other parameters shown to have prognostic importance in PAH, including the RV stroke work index (RVSWi), effective PA elastance (PAE), and PA capacitance (PAC).
      PAPi was analyzed as a continuous variable, and also by dichotomizing it into a binary variable. We used our data to find the cutoff value associated with the optimal c-statistic, and also investigated the use of cutoff values published in previous literature(
      • Lim Y.
      • Low T.T.
      • Chan S.P.
      • Lin W.
      • Teo T.W.
      • Jang J.J.
      • et al.
      Does pulmonary artery pulsatility index predict mortality in pulmonary arterial hypertension?.
      ). As mRAP and CI already have well known reference ranges, we dichotomized these values accordingly. A high mRAP was defined as greater than 8 mmHg and a low CI was defined as less than or equal to 2.4 L/min/m2(
      • Humbert M.
      • Kovacs G.
      • Hoeper M.M.
      • Badagliacca R.
      • Berger R.M.F.
      • Brida M.
      • et al.
      2022 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension.
      ).

      Statistical Analysis

      Categorical variables were summarized as numbers and percentages, whereas continuous variables were summarized as mean values ± standard deviations. Comparisons between groups were done using the chi-square test or the Kruskal-Wallis test, as appropriate.
      To assess the association between each hemodynamic parameter and 1-year mortality, we constructed a multivariate logistic regression model and its corresponding AUC curve; each model was adjusted for age, BMI, sex and advanced functional class.
      For each hemodynamic parameter, we also produced a Cox proportional hazards model adjusted for the same covariates as above to also assess the association between the hemodynamic parameter and 1-year mortality via hazards ratios. These last regressions were performed as a sensitivity analysis to confirm the findings from our logistic regressions despite Cox models being a time to event analysis.
      Statistical significance was determined as a P-value < 0.01.

      Results

      We identified 590 incident patients with PAH from the participating centres. Table 1 demonstrates the baseline characteristics for the total cohort, and for patients who had survived at one year vs. those who did not. The mean age was 61.4±15.5 years, and 66.3% were female. The majority had idiopathic PAH (51.5%). Only 2 patients underwent a lung transplant over the time course of the study. Table 2 shows all hemodynamic variables tested, analyzed as continuous variables. Neither PAPi nor any other of the variables tested were significantly different for survivors vs. non-survivors (P > 0.01 for each).
      Table 1Baseline characteristics of total cohort, and based on survival status at 1 year
      1-year Mortality
      0 (N=547)1 (N=43)Total (N=590)P-value
      Age61.2 (15.44)64.4 (16.19)61.4 (15.51)0.15481
      Female, n (%)365 (66.7%)26 (60.5%)391 (66.3%)0.40302
      PAH subtype, n (%)
      1)iPAH274 (51.3%)23 (53.5%)297 (51.5%)
      2)Heritable3 (0.6%)0 (0.0%)3 (0.5%)
      3)Drug and toxins20 (3.7%)0 (0.0%)20 (3.5%)
      4)CTD187 (35.0%)8 (18.6%)195 (33.8%)
      5)CHD30 (5.6%)3 (7.0%)33 (5.7%)
      6)Portal hypertension11 (2.1%)8 (18.6%)19 (3.3%)
      7)HIV4 (0.7%)0 (0.0%)4 (0.7%)
      Low Risk, functional class ( WHO FC 1,2) n (%)169 (30.9%)11 (25.6%)180 (30.5%)0.46622
      Height, Mean (SD)164.8 (9.81)165.0 (11.48)164.8 (9.92)0.94521
      Weight, Mean (SD)79.2 (22.33)75.8 (20.53)79.0 (22.20)0.30441
      Body mass index, Mean (SD)29.1 (7.73)27.8 (6.70)29.0 (7.66)0.36461
      Lung transplant, n (%)2 (0.4%)0 (0.0%)2 (0.3%)0.69122
      PAPi= Pulmonary artery pulsatility index, IPAH =idiopathic pulmonary arterial hypertension, CTD=connective tissue disease, CHD=congenital heart disease, HIV=human immunodeficiency virus, BMI = body mass index, WHO FC= world health organization functional class, SD=standard deviation
      Table 2Hemodynamic parameters for total cohort, and based on survival status at 1 year
      Death at 1-year
      Total (N=590)No (N=547)Yes (N=43)P-value
      PAWP
      Mean (SD)10.2 (4.75)10.2 (4.76)10.3 (4.76)0.67441
      Median (IQR)10.0 (7.0, 13.0)10.0 (7.0, 13.0)10.0 (6.0, 14.0)
      PVR
      Mean (SD)9.4 (5.50)9.3 (5.53)10.3 (5.10)0.09111
      Median (IQR)8.1 (5.4, 12.2)8.0 (5.3, 12.2)8.7 (6.6, 12.7)
      SvO2
      Mean (SD)61.7 (9.46)62.1 (9.39)57.3 (9.29)0.00141
      Median (IQR)62.8 (56.0, 68.5)63.0 (56.3, 68.9)58.0 (53.2, 65.0)
      mPAP
      Mean (SD)44.8 (11.91)44.5 (12.06)48.3 (9.27)0.01891
      Median (IQR)45.0 (36.0, 52.0)45.0 (35.0, 52.0)47.0 (43.0, 55.0)
      mRAP
      Mean (SD)8.3 (5.16)8.2 (5.07)9.7 (6.09)0.06051
      Median (IQR)7.0 (5.0, 10.0)7.0 (5.0, 10.0)9.0 (6.0, 12.0)
      sPAP
      Mean (SD)73.2 (19.26)72.9 (19.63)77.0 (13.39)0.08991
      Median (IQR)73.0 (60.0, 85.0)73.0 (59.0, 85.0)80.0 (67.0, 85.0)
      Stroke volume
      Mean (SD)57.9 (21.94)58.3 (22.10)53.1 (19.27)0.21881
      Median (IQR)54.8 (41.7, 69.8)55.1 (42.7, 70.0)50.0 (38.6, 65.2)
      Stroke Volume Index
      Mean (SD)30.7 (10.51)30.9 (10.58)28.7 (9.40)0.25151
      Median (IQR)29.6 (23.0, 36.9)29.7 (23.0, 37.1)26.5 (22.1, 34.8)
      Cardiac Index
      Mean (SD)2.3 (0.74)2.3 (0.74)2.2 (0.64)0.89631
      Median (IQR)2.1 (1.8, 2.7)2.1 (1.8, 2.7)2.1 (1.7, 2.7)
      PAPI
      Mean (SD)8.7 (8.63)8.6 (8.41)8.8 (11.22)0.10951
      Median (IQR)6.1 (4.2, 9.7)6.2 (4.2, 9.7)4.7 (3.4, 9.8)
      RV stroke work index
      Mean (SD)151.7 (122.87)149.2 (120.94)184.3 (143.46)0.05841
      Median (IQR)120.7 (70.3, 200.4)119.2 (69.6, 197.9)156.3 (110.1, 221.2)
      Effective PA Elastance
      Mean (SD)0.9 (0.51)0.9 (0.51)1.1 (0.51)0.02221
      Median (IQR)0.8 (0.5, 1.2)0.8 (0.5, 1.2)1.0 (0.7, 1.3)
      PA capacitance
      Mean (SD)1.6 (3.26)1.6 (3.38)1.2 (0.65)0.09071
      Median (IQR)1.2 (0.9, 1.8)1.2 (0.9, 1.8)1.0 (0.8, 1.5)
      1Kruskal-Wallis p-value;
      As the purpose of this work was to evaluate PAPi specifically, we evaluated whether dichotomizing PAPI into low PAPi vs. high PAPi gave a stronger association with mortality. We investigated the ideal value of PAPi that would perform best in a multivariate model. The ideal cutoff value for our data was 5.1, giving an area under the curve of 0.66. However, as this value was very close to the cuttoff value of 5.3 used in previous work(
      • Lim Y.
      • Low T.T.
      • Chan S.P.
      • Lin W.
      • Teo T.W.
      • Jang J.J.
      • et al.
      Does pulmonary artery pulsatility index predict mortality in pulmonary arterial hypertension?.
      ), we decided to use that same cutoff to maintain consistency in the literature.
      Table 3 shows the characteristics for low PAPi and high PAPi, respectively. At one year, 43 (7%) of the patients had died. Table 4 shows the number of deaths for low and high PAPi, CI and mRAP. In the univariate analysis, a PAPi of <5.3 was seen in 58% of patients that died at 1 year compared to 40% of surviving patients (P=0.02). Similarly, a mRAP > 8 mmHg was seen in 65% of patients who died and 47% of patients who survived (P=0.03). Neither baseline functional class nor cardiac index > 2.4 L/min/m2 were associated with mortality at 1 year, nor was age as a continuous variable.
      Table 3Characteristics of patients with high vs. low PAPI
      Low PAPi (N=246)High PAPi (N=345)P-value
      Age, Mean (SD)60.8 (16.34)61.9 (14.87)0.6326
      Female, n (%)161 (65.4%)231 (67.0%)0.7019
      Etiology of PAH, n (%)0.8547
      1)iPAH129 (52.4%)168 (48.7%)
      2)Heritable1 (0.4%)2 (0.6%)
      3)Drugs and toxins11 (4.5%)9 (2.6%)
      4)CTD78 (31.7%)118 (34.2%)
      5)CHD12 (4.9%)21 (6.1%)
      6)Portal hypertension7 (2.8%)12 (3.5%)
      7)HIV1 (0.4%)3 (0.9%)
      Low risk functional class (1,2) n (%)67 (27.2%)114 (33.0%)0.1311
      Height, Mean (SD)164.7 (9.95)164.8 (9.91)0.8999
      Weight, Mean (SD)83.9 (25.19)75.5 (19.10)0.0003
      BMI, Mean (SD)30.9 (8.89)27.7 (6.37)<.0001
      Lung transplant, n (%)2 (0.8%)0 (0.0%)0.0934
      PAPi= Pulmonary artery pulsatility index, IPAH =idiopathic pulmonary arterial hypertension, CTD=connective tissue disease, CHD=congenital heart disease, HIV=human immunodeficiency virus, BMI = body mass index
      Table 4Haemodynamic Parameters by 1-year mortality
      Death at 1 year
      No (N=548)Yes (N=43)TotalP-value
      PAPI< 5.3, n (%)221 (89.8%)25 (10.2%)2460.02
      PAPI>=5.3, n (%)327 (94.8%)18 (5.2%)345
      mRAP< 8, n (%)292 (95.1%)15 (4.9%)3070.02
      mRAP>=8, n (%)256 (90.1%)28 (9.9%)284
      CI< 2.4, n (%)350 (92.6%)28 (7.4%)3780.87
      CI>=2.4, n (%)198 (93.0%)15 (7%)213
      1Chi-Square p-value;
      PAPi = pulmonary artery pulsatility Index, mRAP = mean right atrial pressure, CI = cardiac index
      We constructed receiver operator curves to investigate the ability of PAPi and other parameters to discriminate 1-year mortality. In a multivariate model including age, sex, BMI and functional class, PAPi was associated with 1-year mortality with an AUC of 0.64, 95% CI 0.55-0.73. (Table 5, Figure 1). Similarly, mRAP > 8 mmHg was also able to predict mortality with a similar AUC of 0.65, 95% CI 0.56-0.74. A cardiac index of < 2.4 was not predictive of mortality. Kaplan-Meier curves were created for high vs low PAPi and mRAP (Figure 2). As shown, survival curves for PAPi < 5.3 and mRAP > 8 up to 1 year post catheterization were similar. Table 3 shows the adjusted estimates for the logistic regression and Cox Proportional Hazards Models for each of the three hemodynamic variables evaluated.
      Table 5Logistic and Cox proportional hazards model adjusted estimates for haemodynamic parameters
      Logistic RegressionCox PH
      HaemodynamicAdjusted estimatesAdjusted estimates
      Parameter:OR (95% CI)P-valueAUC (95% CI)HR (95% CI)P-value
      PAPi<5.32.4 (1.2 ,4.6)0.0090.64 (0.55, 0.73)2.5 (1.3 ,4.6)0.006
      mRAP>=82.5 (1.3 ,4.8)0.0090.65 (0.56 ,0.74)2.56 (1.32 ,4.98)0.006
      Cardiac index>=2.41.1 (0.6 ,2.2)0.7930.59 (0.51 ,0.68)0.99 (0.51 ,1.92)0.971
      PAPi = Pulmonary artery pulsatility index, mRAP=mean right atrial pressure, AUC = area under curve, Cox PH = Cox proportional hazards
      Figure thumbnail gr1
      Figure 1Receiver operating curves for PAPi, RAP and CI
      Figure thumbnail gr2
      Figure 2Kaplan-Meier curves for PAPI and RAP
      We performed a subgroup analysis, stratifying for sex. The results for PAPi are shown in Table 6. Given the low number of events seen for both females and males, the confidence intervals are wide, and P was > 0.01 for both. The model appeared to perform better for males compared to females, with a higher AUC (0.73 vs.0.66). However, this difference was not related to PAPi, but the other parameters in the model.
      Table 6Logistic and Cox proportional hazards model adjusted estimates for PAPi, stratified by sex
      Logistic RegressionCox PH
      Adjusted estimatesAdjusted estimates
      OR (95% CI)P-valueAUC (95% CI)HR (95% CI)P-value
      PAPi<5.3Total Cohort2.4 (1.2 ,4.6)0.0090.64 (0.55, 0.73)2.5 (1.3 ,4.6)0.006
      PAPi<5.3Females1.79 (0.79 ,4.06)0.160.66 (0.54 ,0.77)1.9 (0.9 ,4.3)0.11
      PAPi<5.3Males3.85 (1.26 ,11.78)0.0180.73 (0.61 ,0.85)3.5 (1.2 ,9.9)0.02
      As shown in Table 7, we did notice an association between BMI and both mRAP and PAPi. As BMI increased, we saw strong association with increasing mRAP, and a decrease in the median PAPi. The association between BMI appeared stronger with mRAP compared to PAPi. However, the mean PAPi still remained in the low risk range for the highest BMI category.
      Table 7
      BMI (kg/m2)
      < 18.5 (N=20)18.5-25 (N=181)25-30 (N=180)>30 (N=209)Total (N=590)P-value
      mRAP<.00011
      Mean (SD)6.2 (3.11)7.2 (4.74)8.2 (4.98)9.4 (5.57)8.3 (5.16)
      Median (IQR)6 (3.5, 7.5)6 (
      • Galie N.
      • Channick R.N.
      • Frantz R.P.
      • Grunig E.
      • Jing Z.C.
      • Moiseeva O.
      • et al.
      Risk stratification and medical therapy of pulmonary arterial hypertension.
      ,
      • Lim Y.
      • Low T.T.
      • Chan S.P.
      • Lin W.
      • Teo T.W.
      • Jang J.J.
      • et al.
      Does pulmonary artery pulsatility index predict mortality in pulmonary arterial hypertension?.
      )
      7.0 (5.0, 11.0)9 (
      • Sachdev A.
      • Villarraga H.R.
      • Frantz R.P.
      • McGoon M.D.
      • Hsiao J.F.
      • Maalouf J.F.
      • et al.
      Right ventricular strain for prediction of survival in patients with pulmonary arterial hypertension.
      ,
      • Kang G.
      • Ha R.
      • Banerjee D.
      Pulmonary artery pulsatility index predicts right ventricular failure after left ventricular assist device implantation.
      )
      7 (
      • Galie N.
      • Humbert M.
      • Vachiery J.L.
      • Gibbs S.
      • Lang I.
      • Torbicki A.
      • et al.
      2015 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension: The Joint Task Force for the Diagnosis and Treatment of Pulmonary Hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS): Endorsed by: Association for European Paediatric and Congenital Cardiology (AEPC), International Society for Heart and Lung Transplantation (ISHLT).
      ,
      • Lim Y.
      • Low T.T.
      • Chan S.P.
      • Lin W.
      • Teo T.W.
      • Jang J.J.
      • et al.
      Does pulmonary artery pulsatility index predict mortality in pulmonary arterial hypertension?.
      )
      PAPI0.00141
      Mean (SD)8.3 (4.16)10.0 (10.22)8.0 (6.67)8.1 (8.86)8.7 (8.63)
      Median (IQR)7.0 (5.7, 10.1)6.7 (4.6, 11.2)5.8 (4.2, 9.4)5.3 (3.5, 8.8)6.1 (4.2,9.7)
      1Kruskal-Wallis p-value;

      Discussion

      In this large, multicenter Canadian study of 590 patients with incident PAH, we found that PAPi was not associated with survival when analyzed as a continuous variable. However, when analyzed as low vs. high PAPi with a cutoff of 5.3, low PAPi was associated with a higher 1-year mortality. The ability of PAPi to discriminate 1-year mortality was no different than that of RAP alone.
      This is the largest study to evaluate the clinical significance and prognostic utility of PAPi in a PAH population. The PAPi is a novel hemodynamic parameter of RV-PA coupling that has recently been found to be useful in predicting outcomes in a variety of states influenced by RV function. It predicts right ventricular (RV) dysfunction in patients with acute inferior wall myocardial infarction or after left ventricular assist device implantation(
      • Korabathina R.
      • Heffernan K.S.
      • Paruchuri V.
      • Patel A.R.
      • Mudd J.O.
      • Prutkin J.M.
      • et al.
      The pulmonary artery pulsatility index identifies severe right ventricular dysfunction in acute inferior myocardial infarction.
      ,
      • Kang G.
      • Ha R.
      • Banerjee D.
      Pulmonary artery pulsatility index predicts right ventricular failure after left ventricular assist device implantation.
      ,
      • Morine K.J.
      • Kiernan M.S.
      • Pham D.T.
      • Paruchuri V.
      • Denofrio D.
      • Kapur N.K.
      Pulmonary Artery Pulsatility Index Is Associated With Right Ventricular Failure After Left Ventricular Assist Device Surgery.
      ). A more recent study has reported PAPi to be a significant predictor of mortality and hospitalization in individuals with advanced heart failure(
      • Kochav S.M.
      • Flores R.J.
      • Truby L.K.
      • Topkara V.K.
      Prognostic Impact of Pulmonary Artery Pulsatility Index (PAPi) in Patients With Advanced Heart Failure: Insights From the ESCAPE Trial.
      ). Additionally, a lower PAPi was associated with decreased left ventricular ejection fraction, more severe tricuspid regurgitation, inferior vena cava dilatation, and lower cardiac index. In PAH, the literature on PAPi is limited to small and predominantly single center studies. Lim et al investigated the association between PAPi and survival in a population of 102 patients with PAH from Singapore(
      • Lim Y.
      • Low T.T.
      • Chan S.P.
      • Lin W.
      • Teo T.W.
      • Jang J.J.
      • et al.
      Does pulmonary artery pulsatility index predict mortality in pulmonary arterial hypertension?.
      ). These patients were recruited over a large time range from 2003 to 2016, during which therapeutic options and use combination therapy changed significantly. Another study by Mazimba and colleagues evaluated the association between PAPi and survival in the historical National Institutes of Health (NIH) primary PH database in which baseline PAPi was predictive of survival(
      • Mazimba S.
      • Welch T.S.
      • Mwansa H.
      • Breathett K.K.
      • Kennedy J.L.W.
      • Mihalek A.D.
      • et al.
      Haemodynamically Derived Pulmonary Artery Pulsatility Index Predicts Mortality in Pulmonary Arterial Hypertension.
      ). However, it should be noted that the NIH study cohort predates the development of modern PAH therapy and is therefore of limited relevance in the modern care of PAH. We limited our population to the modern era of combination therapy use (after the publication of the AMBITION trial and the 2015 European Society of Cardiology/European Respiratory Society Guidelines(
      • Galie N.
      • Humbert M.
      • Vachiery J.L.
      • Gibbs S.
      • Lang I.
      • Torbicki A.
      • et al.
      2015 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension: The Joint Task Force for the Diagnosis and Treatment of Pulmonary Hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS): Endorsed by: Association for European Paediatric and Congenital Cardiology (AEPC), International Society for Heart and Lung Transplantation (ISHLT).
      ,
      • Galie N.
      • Barbera J.A.
      • Frost A.E.
      • Ghofrani H.A.
      • Hoeper M.M.
      • McLaughlin V.V.
      • et al.
      Initial Use of Ambrisentan plus Tadalafil in Pulmonary Arterial Hypertension.
      )), making these results more applicable to modern practice.
      In our study, PAPi was associated with 1-year mortality only when dichotomized. Moreover, many other prognostic parameters were not associated with morality at 1 year. The primary reason for this discrepancy is likely low statistical power, as only 7% mortality was noted in 1 year. When PAPi was dichotomized, it was associated with mortality, but the association was not stronger than mPAP alone. Similar findings were seen in the Lim et al. study. We chose a shorter term follow up of 1 year in our study, understanding that response to PAH therapy would likely be a greater predictor of survival over the long term compared to baseline hemodynamic data. Our analysis adjusted for age, sex, BMI and functional class at baseline, and found a similar AUC for PAPi as seen by Lim and colleagues. Interestingly, neither our study nor the Lim et al. study showed an association between baseline cardiac index and mortality despite the fact it is an established prognostic variable from older studies(
      • McLaughlin V.V.
      • Shillington A.
      • Rich S.
      Survival in primary pulmonary hypertension: the impact of epoprostenol therapy.
      ,
      • Sitbon O.
      • Humbert M.
      • Nunes H.
      • Parent F.
      • Garcia G.
      • Herve P.
      • et al.
      Long-term intravenous epoprostenol infusion in primary pulmonary hypertension: prognostic factors and survival.
      ). Data from the French PAH registry in the modern era also found that baseline hemodynamic variables were not associated with long-term survival whereas follow-up hemodynamic variables were(
      • Weatherald J.
      • Boucly A.
      • Chemla D.
      • Savale L.
      • Peng M.
      • Jevnikar M.
      • et al.
      Prognostic Value of Follow-Up Hemodynamic Variables After Initial Management in Pulmonary Arterial Hypertension.
      ,
      • Weatherald J.
      • Boucly A.
      • Launay D.
      • Cottin V.
      • Prevot G.
      • Bourlier D.
      • et al.
      Haemodynamics and serial risk assessment in systemic sclerosis associated pulmonary arterial hypertension.
      ). Together, these data support the concept that prognosis is not so much governed by baseline hemodynamics as it is response to therapy. We strongly suspect that differential changes in hemodynamics with initiation of PAH therapy also explains the limited prognostic utility of baseline PAPi and the other hemodynamic parameters tested. However, once decompensation and right sided heart failure with elevated RA pressures has developed, there may a signal of worse survival that persists despite initiation of therapy.
      This study has some limitations to acknowledge. The design was retrospective and observational based on registry data, which has inherent limitations and bias. Additionally, lack of follow up RHC is a limitation as it was not routinely done in each centre, making it impossible to evaluate associations between survival and follow up PAPi or change in PAPi after initial therapy. We were not powered to investigate other outcomes of interest, such as lung transplantation or hospital admission for PAH. We did not employ a competing risks model for lung transplantation and mortality, however the number of patients transplanted within 1 year was so low that omitting this is justifiable. The major strength of our study is the large multicenter nature, with well phenotyped patients who were diagnosed and treated in PAH expert centres. An additional strength is the limitation of inclusion to patients diagnosed in the modern era of initial combination therapy.

      Conclusions

      As in other studies, we found the association between PAPi and survival was modest, and not superior to RA pressure alone. Our findings cast further doubt as to whether PAPi is a useful measure in PAH, and on the long-term prognostic utility of baseline hemodynamic measures in the era of effective early combination therapy.

      Uncited reference

      • Nie L.
      • Li J.
      • Zhang S.
      • Dong Y.
      • Xu M.
      • Yan M.
      • et al.
      Correlation between right ventricular-pulmonary artery coupling and the prognosis of patients with pulmonary arterial hypertension.
      .

      Acknowledgements:

      The Canadian Pulmonary Hypertension registry was supported by grants from Janssen, United Therapeutics and Bayer Pharmaceuticals.

      References

        • Parikh V.
        • Bhardwaj A.
        • Nair A.
        Pharmacotherapy for pulmonary arterial hypertension.
        J Thorac Dis. 2019; 11: S1767-S1781
        • Boucly A.
        • Weatherald J.
        • Savale L.
        • de Groote P.
        • Cottin V.
        • Prevot G.
        • et al.
        External validation of a refined four-stratum risk assessment score from the French pulmonary hypertension registry.
        Eur Respir J. 2022; 59
        • Zelt J.G.E.
        • Sugarman J.
        • Weatherald J.
        • Partridge A.C.R.
        • Liang J.C.
        • Swiston J.
        • et al.
        Mortality trends in pulmonary arterial hypertension in Canada: a temporal analysis of survival per ESC/ERS guideline era.
        Eur Respir J. 2022; 59
        • Galie N.
        • Channick R.N.
        • Frantz R.P.
        • Grunig E.
        • Jing Z.C.
        • Moiseeva O.
        • et al.
        Risk stratification and medical therapy of pulmonary arterial hypertension.
        Eur Respir J. 2019; 53
        • Galie N.
        • Humbert M.
        • Vachiery J.L.
        • Gibbs S.
        • Lang I.
        • Torbicki A.
        • et al.
        2015 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension: The Joint Task Force for the Diagnosis and Treatment of Pulmonary Hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS): Endorsed by: Association for European Paediatric and Congenital Cardiology (AEPC), International Society for Heart and Lung Transplantation (ISHLT).
        Eur Respir J. 2015; 46: 903-975
        • Sachdev A.
        • Villarraga H.R.
        • Frantz R.P.
        • McGoon M.D.
        • Hsiao J.F.
        • Maalouf J.F.
        • et al.
        Right ventricular strain for prediction of survival in patients with pulmonary arterial hypertension.
        Chest. 2011; 139: 1299-1309
        • Tello K.
        • Dalmer A.
        • Axmann J.
        • Vanderpool R.
        • Ghofrani H.A.
        • Naeije R.
        • et al.
        Reserve of Right Ventricular-Arterial Coupling in the Setting of Chronic Overload.
        Circ Heart Fail. 2019; 12e005512
        • Nie L.
        • Li J.
        • Zhang S.
        • Dong Y.
        • Xu M.
        • Yan M.
        • et al.
        Correlation between right ventricular-pulmonary artery coupling and the prognosis of patients with pulmonary arterial hypertension.
        Medicine (Baltimore). 2019; 98e17369
        • Korabathina R.
        • Heffernan K.S.
        • Paruchuri V.
        • Patel A.R.
        • Mudd J.O.
        • Prutkin J.M.
        • et al.
        The pulmonary artery pulsatility index identifies severe right ventricular dysfunction in acute inferior myocardial infarction.
        Catheter Cardiovasc Interv. 2012; 80: 593-600
        • Lim Y.
        • Low T.T.
        • Chan S.P.
        • Lin W.
        • Teo T.W.
        • Jang J.J.
        • et al.
        Does pulmonary artery pulsatility index predict mortality in pulmonary arterial hypertension?.
        ESC Heart Fail. 2021; 8: 3835-3844
        • Humbert M.
        • Kovacs G.
        • Hoeper M.M.
        • Badagliacca R.
        • Berger R.M.F.
        • Brida M.
        • et al.
        2022 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension.
        Eur Heart J. 2022; 43: 3618-3731
        • Kang G.
        • Ha R.
        • Banerjee D.
        Pulmonary artery pulsatility index predicts right ventricular failure after left ventricular assist device implantation.
        J Heart Lung Transplant. 2016; 35: 67-73
        • Morine K.J.
        • Kiernan M.S.
        • Pham D.T.
        • Paruchuri V.
        • Denofrio D.
        • Kapur N.K.
        Pulmonary Artery Pulsatility Index Is Associated With Right Ventricular Failure After Left Ventricular Assist Device Surgery.
        J Card Fail. 2016; 22: 110-116
        • Kochav S.M.
        • Flores R.J.
        • Truby L.K.
        • Topkara V.K.
        Prognostic Impact of Pulmonary Artery Pulsatility Index (PAPi) in Patients With Advanced Heart Failure: Insights From the ESCAPE Trial.
        J Card Fail. 2018; 24: 453-459
        • Mazimba S.
        • Welch T.S.
        • Mwansa H.
        • Breathett K.K.
        • Kennedy J.L.W.
        • Mihalek A.D.
        • et al.
        Haemodynamically Derived Pulmonary Artery Pulsatility Index Predicts Mortality in Pulmonary Arterial Hypertension.
        Heart Lung Circ. 2019; 28: 752-760
        • Galie N.
        • Barbera J.A.
        • Frost A.E.
        • Ghofrani H.A.
        • Hoeper M.M.
        • McLaughlin V.V.
        • et al.
        Initial Use of Ambrisentan plus Tadalafil in Pulmonary Arterial Hypertension.
        N Engl J Med. 2015; 373: 834-844
        • McLaughlin V.V.
        • Shillington A.
        • Rich S.
        Survival in primary pulmonary hypertension: the impact of epoprostenol therapy.
        Circulation. 2002; 106: 1477-1482
        • Sitbon O.
        • Humbert M.
        • Nunes H.
        • Parent F.
        • Garcia G.
        • Herve P.
        • et al.
        Long-term intravenous epoprostenol infusion in primary pulmonary hypertension: prognostic factors and survival.
        J Am Coll Cardiol. 2002; 40: 780-788
        • Weatherald J.
        • Boucly A.
        • Chemla D.
        • Savale L.
        • Peng M.
        • Jevnikar M.
        • et al.
        Prognostic Value of Follow-Up Hemodynamic Variables After Initial Management in Pulmonary Arterial Hypertension.
        Circulation. 2018; 137: 693-704
        • Weatherald J.
        • Boucly A.
        • Launay D.
        • Cottin V.
        • Prevot G.
        • Bourlier D.
        • et al.
        Haemodynamics and serial risk assessment in systemic sclerosis associated pulmonary arterial hypertension.
        Eur Respir J. 2018; 52