Abstract
Different Doppler echocardiography (DE) models have been proposed for estimation of mean pulmonary arterial pressures (PAMP) from tricuspid regurgitation (TR) jet velocity. We aimed to compare four TR-derived DE models in predicting the PAMP measured by right heart catheterization (RHC) in different groups of precapillary pulmonary hypertension (PH). A total of 287 patients with hemodynamically pre-capillary PH were enrolled (mean age = 51 ± 17.4 years, 59.9% female). All patients underwent DE before RHC (< 3 h) and four formulae (F) were used for TR-derived PAMP estimation (PAMP-DE). These were as follows: F1 = Chemla (0.61 × systolic pulmonary artery pressure [PASP] + 2); F2 = Friedberg (0.69 × PASP − 0.22), F3 = Aduen (0.70 × PASP); and F4 = Bech-Hanssen (0.65 × PASP − 1.2). The PASP and PAMP (mmHg) measured by RHC were 89.1 ± 30.4 and 55.8 ± 20.8, respectively. In the overall PH group, DE estimates for PASP (r = 0.59,
Keywords
Introduction
According to the current practice guidelines, pulmonary hypertension (PH) is defined as an increase in mean pulmonary arterial (PA) pressure (PAMP) ≥ 25 mmHg at rest as assessed by right heart catheterization (RHC). Although RHC is the method of choice for the diagnosis of PH, continuous wave Doppler echocardiography (DE) is frequently used as an initial screening method to estimate pulmonary pressures.1–3 Moreover, DE has been used as a non-invasive, inexpensive, and widely available method of monitoring disease progression over time in patients with PH. The most commonly used and studied method is the estimation of systolic PA pressure (PASP) by measuring the peak velocity of the tricuspid regurgitant (TR) jet using a modified Bernoulli equation.1–11 In this approach, right ventricular systolic pressure or PASP can then be estimated by adding the TR systolic pressure gradient to the estimated right atrial pressure (RAP) and this pressure is equal to PASP in the absence of pulmonic stenosis. However, the reliability of RAP estimation is considered to be an important limitation of the peak TR jet velocity method.2,4–6 Furthermore, accuracy of this non-invasive estimation of PAMP for initial screening and follow-up has remained an unmet need in this setting.12–21 The majority of DE methods proposed for the estimation of PAMP from TR jet velocity (PAMP-DE) have been based on the assumption of linearity between PASP and PAMP.12–23 However, utility of these DE methods across the subgroup of PH has not been evaluated.
The purpose of our study was to assess the reliability of the PASP and PAMP values estimated by TR jet-derived DE methods for corresponding PA pressures measured with RHC in an overall PH population and prespecified subgroups of PH.
Methods
In this study, a total of 287 patients with hemodynamically pre-capillary PH were enrolled (mean age = 51 ± 17.4 years, 59.9% female). Subgroups were as follows: idiopathic pulmonary arterial hypertension (IPAH) (n = 48); PAH associated with congenital heart disease and connective tissue disease (APAH-CHD, APAH-CTD) (n = 106 and 17); chronic thromboembolic PH (CTEPH) (n = 64); and other PH groups (n = 52). Hemodynamic definitions of pre-capillary PH and clinical diagnostic algorithms were initially performed according to the 2009 European Society of Cardiology / European Respiratory Society (ESC/ERS) PH Guidelines and retrospectively re-evaluated after the publication of the 2015 ESC/ERS PH Guidelines.1,2 Patients in whom PAMP was measured to be < 25 mmHg, which is the definitive threshold for PH, and patients with pulmonary arterial wedge pressure (PAWP) ≥ 15, indicating postcapillary PH, were excluded from the study.
The PAMP-DE was indirectly derived from PASP using the four formulae (F) as follows: F1 = Chemla (0.61 ×PASP + 2); 13 F2 = Friedberg (0.69 × PASP–0.22); 15 F3 =Aduen (0.70 × PASP); 16 and F4 = Bech-Hanssen (0.65 ×PASP – 1.2). 18
The study protocol was approved by the Institutional Review Board and written informed consent was obtained from all patients before the DE and RHC assessments.
Statistical methods
Continuous variables were presented as mean ± standard deviation (SD); categorical variables were presented as percentage. Linear association between variables was assessed using Pearson correlation coefficient. Paired-t test was performed between means (PAMP-RHC vs. PAMP estimated by four DE equations) to detect any possible systematic bias. The degree of agreement (DOA) between continuous measurements of PAMP-RHC and four DE formulae were quantified through intraclass correlation coefficient (ICC), Lin’s concordance correlation coefficient, and Bland–Altman (BA) analysis. ICC < 0.4 indicates poor agreement, ICC in the range of 0.4–0.75 indicates moderate agreement, and ICC values > 0.75 indicate excellent agreement. The DOA between PAMP-RHC and PAMP-DE by four equations categorized into deciles was quantified through Kappa statistics: Kappa value < 0.2 indicates poor agreement; 0.2–0.4 fair agreement; 0.4–0.6 moderate agreement; 0.6–0.8 substantial agreement; and > 0.8 almost perfect agreement.
24
Also, to assess the DOA between continuous variables, the Passing–Bablok non-parametric regression method was used. Passing–Bablok non-parametric regression provides intercept and slope with 95% confidence intervals (CI).
25
The intercept assesses systematic bias; a CI below zero denotes a constant underestimation trend whereas values over zero denote an overestimation trend. The slope assesses proportional bias; a good agreement is present when its CI includes 1. In BA analysis, PAMP-RHC and PAMP estimated by four DE formulae plotted with line represen mean difference ± 1.96 SD (95% limits of agreement).
26
A
Results
Baseline clinical, echocardiographic, and invasive characteristics.
IPAH, idiopathic pulmonary arterial hypertension; APAH, associated pulmonary arterial hypertension; CTD, connective tissue disorders; CHD, congenital heart disease; CTEPH, chronic thromboembolic pulmonary hypertension; NYHA, New York Heart Association; PASP, systolic pulmonary arterial pressure; TAPSE, tricuspid annular plane systolic excursion; LVEF, left ventricular ejection fraction; PADP, diastolic pulmonary arterial pressure; PAMP, mean pulmonary arterial pressure; PVR, pulmonary vascular resistance; SVR, systemic vascular resistance.
In the overall PH group, PASP-DE showed significant correlation to PASP-RHC (r = 0.593, Scatter plot between PASP measured by RHC and echocardiography (some spots are darker, which indicate that some patients have overlapping values).
The ICC, Lin’s concordance coefficient, and Pearson correlation coefficient between PAMP-RHC and PAMP-DE by four equations were presented in Table 2. The cumulative frequency and individual values of PAMP-DE obtained by four equations vs. PAMP-RHC were presented in Figs. 2 and 4, respectively. In the visual assessment of cumulative frequency plot (Fig. 2), the Chemla, Friedberg, Aduen, and Bech-Hanssen equations provided concordant results, respectively.
The cumulative frequency of PAMP obtained by four equations vs. RHC. The individual value of PAMP obtained by four equations vs. RHC (yellow colored areas indicate the mean value; red boxes indicate the median and interquartile range; and green range lines indicate standard error of the mean). Pearson correlation and absolute agreement (ICC and Lin’s coefficient) between PAMP_Catheter and different echocardiographic equations. PAMP, mean pulmonary arterial pressure; ICC, intraclass correlation coefficient.

Results of paired t test for PAMP_Catheter compared to each PAMP calculation by formula.
CI, confidence interval; SD, standard deviation.
Passing–Bablok non-parametic analysis for each PAMP estimating formula.
The cusum test was used for linearity.
PAMP, mean pulmonary arterial pressure; CI, confidence interval.

Passing-Bablok regression plot for PAMP obtained by four equations vs. RHC.
Bland–Altman analysis for measured PAMP by catheterization vs. different PAMP calculated by formulae using echocardiography in the overall population and in IPAH, APAH, and CTEPH populations.
IPAH, idiopathic pulmonary hypertension; APAH, associated pulmonary hypertension; CHD, congenital heart disease; CTEPH, chronic thromboembolic pulmonary hypertension; LOA, limit of agreement.

Bland–Altman plot for PAMP obtained by four equations vs. RHC in the overall population.
Discussion
In this study comprising patients with confirmed PH, DE estimates for PASP and PAMP from TR jet velocity showed good and significant correlations to invasively measured PASP and PAMP, respectively. Among the PAMP estimates by four DE formulae, concordant results were noted between the Chemla and Bech-Hanssen equations and the Friedberg and Auden equations, respectively. However, agreements among four DE formulae were found to be poor and the Chemla and Bech-Hanssen formulae underestimated PAMP-RHC while other two overestimated PAMP-RHC.
The estimation of PASP and PAMP from TR jet velocity suffers from several limitations.1–22 Therefore, PH can not be reliably defined by a specific threshold cut-off for TR jet peak velocity and estimation of PASP based solely on this measurement may not be appropriate to screen for mild, asymptomatic PH. In a large review including 29 studies, the correlation between PASP-DE and PASP-RHC was found to be modest, with a summary correlation coefficient of 0.70 (95% CI = 0.67–0.73), regardless of the disease severity. The diagnostic accuracy (DA) of DE for PH proven by RHC was also found to be modest, even adopting a cut-off value of 40 mmHg that is the most commonly used threshold for PASP by DE studies and some risks for false-positive and false-negative results also arose. 11 Although significant heterogeneity was evident in both correlation and DA analyses, any source for the heterogenity was not found by sensitivity analysis. 1 The possible source for the heterogenity for both analyses was associated with a combination of variation in study populations in terms of co-morbidity, study design, spectrum bias, disease progression bias, review and publication bias, and population bias. 11 In the review by Janda et al., spectrum bias and disease progression bias were noted in some studies. 15 However, because of the lack of information regarding blinding during testing, the possibility of review bias remains unclear whereas poor reporting of any co-morbid conditions within the PH group seemed to limit generalizability. 11 Furthermore, intra- and inter-observer variability during DE assessment and reading the recordings, physiological variability in PA hemodynamics among repeat measurements related with volume status, systemic blood pressure and changes in blood oxygenation should be considered as the confounding factors even by the same method on the same day.1–3,11
The mathematical and physiological links between PASP and PAMP are crucially important and a strong linearity between PASP and PAMP has been shown in several studies performed in adults and pediatric populations.12,13,22,23 This relationship was reported to be unchanged in several situations including exercise, pacing, treatment with inotropic agents or pulmonary vasodilators, following heart or double lung transplantation. 23 According to the classical physiologic approach based on empirical formula of the two-pressure model, PAMP is defined by the following equation: PAMP = 2/3 × PADP + 1/3 × PASP, where PADP refers to diastolic PA pressure. 27 The PAMP has been documented to be more sensitive to changes in PADP compared to those with PASP, and PADP may represent distal vascular resistance more accurately than PASP that is more dependent to the characteristics of right ventricular ejection dynamics, PA compliance, and wave reflections.14,22,27 Moreover, in a simpler single-pressure model using a high-fidelity manometer a new equation was proposed: PAMP = 0.61 ×PASP + 2 mmHg; PASP was reported to account for 98% of PAMP variability. 12 Chemla et al. analysed five studies in which high-fidelity manometers were utilized and reported that the most accurate estimate for PAMP was found by using PASP only, while combining the PASP with PADP provided the most precise PAMP estimate. 13
Several DE models have been developed for the estimation of the PAMP from TR or pulmonary regurgitant jet velocities.12–20 In addition to the Chemla formula, Syyed (PAMP = 0.65 × PASP + 0.55 mmHg), the Friedberg, Aduen and Bech-Hanssen equations have also been developed for TR-derived DE estimations for PAMP.12–19 A comparison analysis revealed that estimates of the Chemla or Syyed formulae reflect more precisely the PAMP measured by RHC than that measured by PASP-DE.12–14 As another method, PAMP-DE estimation by Aduen’s mean gradient method also showed a similar accuracy and precision compared to the methods of Chemla and Syyed. 17 The acceptable accuracy of these methods were considered to suggest their equal suitability for clinical use.12–19 Although all these DE formulae have been shown to provide high DA for PH, comparison among these equations in the pre-specified subsets of PH has not been reported.
Although our study revealed good and significant correlations between four PAMP-DE estimates and PAMP-RHC (r = 0.563,
Limitations
Because our analyses were performed in a pre-defined population in whom precapillary PH was already confirmed by RHC, the DA of any threshold cut-off for PASP or PAMP defined by DE for PH screening in a general population was not a concern in our study. Moreover, the exclusion of patients without PH as assessed by RHC might cause a spectrum bias carrying the risk for distortions in a diagnostic test’s performance and eventually overestimation of the sensitivity and specificity of the DE. These results represent drawbacks of TR jet-derived DE methods, but not those with other DE methods in which acceleration time of right ventricular outflow ejection or early diastolic pulmonary regurgitant jet velocity are utilized. As another probable limitation, intra- and inter-observer variability analysis in DE assessments were not included in the study. Finally, physiologic variability of the pulmonary vascular hemodynamics at sequential assessments might be a confounding factor for comparisons between DE and RHC measures in the absence of simultaneous evaluations.
Conclusions
The results of this study performed in a pre-defined population with confirmed PH showed a poor agreement between four TR jet-derived PAMP-DE estimates and PAMP-RHC measures, regardless of the underlying pathology. Therefore, these DE estimates solely might not be reliable for periodical non-invasive assessments of PAMP in patients with PH who were under targeted treatments. However, these results should not be extrapolated to DE screening procedures for PH performed in normal population. Whether implementation of other confirmatory echocardiographic findings suggestive for PH may eliminate these limitations of DE remains to be determined.
Footnotes
Conflict of interest
The author(s) declare that there is no conflict of interest.
Funding
Article processing procedure is financially supported by Actelion Pharmaceuticals – Turkey as an unrestricted scientific grant.
