Bidimensional Measurements of Right Ventricular Function for Prediction of Survival in Patients with Pulmonary Hypertension: Comparison of Reproducibility and Time of Analysis with Volumetric Cardiac Magnetic Resonance Imaging Analysis

Abstract

We tested the hypothesis that bidimensional measurements of right ventricular (RV) function obtained by cardiac magnetic resonance imaging (CMR) in patients with pulmonary arterial hypertension (PAH) are faster than volumetric measures and highly reproducible, with comparable ability to predict patient survival. CMR-derived tricuspid annular plane systolic excursion (TAPSE), RV fractional shortening (RVFS), RV fractional area change (RVFAC), standard functional and volumetric measures, and ventricular mass index (VMI) were compared with right heart catheterization data. CMR analysis time was recorded. Receiver operating characteristic curves, Kaplan-Meier, Cox proportional hazard (CPH), and Bland-Altman test were used for analysis. Forty-nine subjects with PAH and 18 control subjects were included. TAPSE, RVFS, RVFAC, RV ejection fraction, and VMI correlated significantly with pulmonary vascular resistance and mean pulmonary artery pressure (all P < 0.05). Patients were followed up for a mean (± standard deviation) of 2.5 ± 1.6 years. Kaplan-Meier curves showed that death was strongly associated with TAPSE <18 mm, RVFS <16.7%, and RVFAC <18.8%. In CPH models with TAPSE as dichotomized at 18 mm, TAPSE was significantly associated with risk of death in both unadjusted and adjusted models (hazard ratio, 4.8; 95% confidence interval, 2.0–11.3; P = 0.005 for TAPSE <18 mm). There was high intra- and interobserver agreement. Bidimensional measurements were faster (1.5 ± 0.3 min) than volumetric measures (25 ± 6 min). In conclusion, TAPSE, RVFS, and RVFAC measures are efficient measures of RV function by CMR that demonstrate significant correlation with invasive measures of PAH severity. In patients with PAH, TAPSE, RVFS, and RVFAC have high intra- and interobserver reproducibility and are more rapidly obtained than volumetric measures. TAPSE <18 mm by CMR was strongly and independently associated with survival in PAH.

Keywords

pulmonary arterial hypertension tricuspid annular plane systolic excursion (TAPSE)cardiac MRI,survival

Pulmonary arterial hypertension (PAH) is characterized by pulmonary vascular remodeling that results in a progressive increase in pulmonary artery pressure and vascular resistance, eventually leading to right ventricular (RV) failure. Increased intracavitary pressure in the RV initially leads to elevated wall stress and adaptive myocardial hypertrophy, followed by contractile dysfunction and maladaptive cavity dilatation with disease progression. RV response to the increased afterload is a major determinant of survival in subjects with PAH.^1–3 Thus, serial assessment of RV function with indices such as RV ejection fraction (RVEF) may be important for patient management.⁴

Echocardiography is routinely used in the clinical assessment of RV function (RVF) in PAH. However, evaluation of RVF with echocardiography is challenging due to the complex shape of the RV and its trabeculated and thin myocardium.⁵ In addition, optimal imaging windows may not be attainable in up to 25% of patients. In spite of these limitations, echocardiography is extremely valuable as a screening modality for the diagnosis of PAH because of its wide availability and low cost. Earlier studies have shown that tricuspid annular plane systolic excursion (TAPSE) by echocardiography is a fast and easily obtainable measure of RVF⁶ that predicts survival in patients with PAH.^7–9

Reliable, noninvasive evaluation of the RV in patients with suspected PAH is feasible with cardiac magnetic resonance imaging (CMR). CMR accurately assesses RV structure and function, providing precise bidimensional (2D) and three-dimensional (3D) volumetric measurements of both ventricles. CMR has been validated for the RV and is considered the gold standard method for determining RVF.^10–12 CMR-derived ventricular mass index (VMI) has been shown to predict prognosis in PAH.^13,14 However, calculation of VMI requires time-consuming manual or semiautomated contouring of the left ventricle (LV) and RV. On the other hand, TAPSE by CMR is a 2D measurement that is simple and quick to obtain and has a high correlation with RVEF in healthy subjects.^15,16 Therefore, we hypothesized that 2D measures by CMR would have the same prognostic value as 3D measures, such as VMI, but would be faster to obtain and more reproducible. The objectives of this study were to compare the utility of 2D measurements of RVF in assessment of prognosis in PAH by evaluating 2D and 3D measures in healthy control subjects and patients with PAH; comparing 2D measurements of RVF (TAPSE and others) and standard 3D RVF measurements with invasive measurements obtained by right heart catheterization (RHC); determining whether 2D measures can be used for prediction of survival and whether they are equivalent to the 3D approach; and determining inter- and intraobserver agreement and time needed for analysis for 2D and 3D measures by CMR.

METHODS

This prospective study was Health Insurance Portability and Accountability Act compliant and was approved by our institutional review board. From April 2004 to February 2013, 49 consecutive patients with PAH were enrolled after referral for a clinically indicated RHC. Inclusion criteria were PAH by RHC, defined as a mean pulmonary artery pressure (mPAP) ≥25 mmHg with a pulmonary capillary wedge pressure (PCWP) ≤15 mmHg¹⁰ and the ability to undergo CMR. Exclusion criteria included chronic thromboembolic disease and obstructive or restrictive lung disease, as described elsewhere,⁸ among others. Eighteen healthy control subjects underwent the same CMR protocol but were not subjected to RHC. History of previous treatment for PAH was recorded.

CMR acquisition

CMR for this study was performed on a 3-T (MAGNETOM Trio, Siemens Healthcare, Erlangen, Germany) system using a body coil array. Retrospectively gated segmented gradient echo (FLASH) cine short-axis images including both ventricles from base to apex (repetition time [TR]/echo time [TE], 64.2/2.97 milliseconds; slice thickness, 6 mm; flip angle, 18°; acquisition matrix, 256 × 192; field of view [FOV], 360 mm) and cine four-chamber images (TR/TE, 52.2/2.82 milliseconds; slice thickness, 8 mm; flip angle, 12°; acquisition matrix, 256 × 192; FOV, 360 mm) were acquired.

RHC

Subjects with PAH underwent RHC within 24 hours of CMR. Standard hemodynamic parameters were measured. Mean (end-expiratory) right atrial pressure, pulmonary capillary wedge pressure (PCWP), and pulmonary arterial pressure were recorded, and thermodilution cardiac output was obtained and reported as the average of at least three values with <20% variation. Mixed venous oxygen saturation was determined while patients were breathing room air. Heart rate and noninvasive blood pressure were recorded during the procedure. Stroke volume (SV), defined as cardiac output/heart rate; cardiac index, defined as cardiac output/body surface area; stroke volume index (SVI), defined as cardiac index/heart rate; pulmonary vascular resistance (PVR), defined as (mPAP-PCWP)/cardiac output; and RV stroke work index, defined as (mPAP–right atrial pressure) × SVI, were calculated.

Image analysis

Quantification of the RVF was obtained using commercially available software (Qmass version 6.2.2 and Mass 5.1, Medis, Netherlands).

2D analysis

Assessment of the RVF was performed using cine four-chamber images. TAPSE was calculated by measuring the displacement of the lateral tricuspid valve annulus toward the apex between end-diastolic (ED) and end-systolic (ES) phases (Fig. 1). RV fractional shortening (RVFS) was calculated as follows:

RVFS = (\frac{TAPSE}{ED - annulus to apex length}) \times 100.

(1)

Figure 1.

Fifty-year-old male (A, B) and 53-year-old female (C, D) with pulmonary arterial hypertension. End-diastolic (A, C) and end-systolic (B, D) images from a cine four-chamber cardiac magnetic resonance imaging (CMR) acquisition. Tricuspid annulus plane systolic excursion (TAPSE) was determined as the difference between the distance from the right ventricular (RV) apex to tricuspid annulus in end diastole and end systole, measured in millimeters. In first patient (A, B), CMR shows a measurement of TAPSE >18 mm; in the second patient, TAPSE is <18 mm (C, D). LA: left aorta; LV: left ventricle; RA: right aorta.

Right ventricular fractional area change (RVFAC) was calculated by measuring percentage change as follows:

RVFAC = (\frac{ED - RV area - ES - RV area}{ED - RV area}) \times 100.

(2)

3D analysis

Assessment of the RVF was performed using contiguous cine short-axis images covering the entire RV and LV. The RV and LV epicardial and endocardial ventricular borders were manually contoured in end systole (smallest LV volume) and end diastole (largest LV volume) for functional analysis. Papillary muscles and trabeculations were included in the blood pool volume as described elsewhere.^17,18 Morphologic and functional parameters were normalized to body surface area. RVSV and RVEF were calculated as SV = RVED volume - RVES volume, and RVEF (%) = (RVSV/RVED volume) × 100. Ventricular mass index (VMI) was calculated as the ratio between RVED mass and LVED mass.

Reproducibility analysis

A subset of 36 patients was analyzed by two readers (with 2 years and 8 months of experience in CMR, respectively). Imaging analysis time was recorded for each case. Intraobserver agreement for 2D and 3D measures was assessed by blinded repeat analysis of the images 1 month after the initial analysis by reader 1.

Statistical analysis was performed with MedCalc (Ostend, Belgium). Variables were summarized by mean ± standard deviation and were compared using independent t test. Spearman correlation coefficient was used to compare the relationship between 2D and 3D measures and hemodynamic measurements. A univariate and repeated bivariate Cox regression analysis was used to look at the relative strength of association of continuous variables with survival. Survival analysis was performed for 2D, RHC, and 3D measures given previous reports demonstrating their prognostic significance^7,14 by receiver operating characteristic (ROC) curves. The proportional hazards assumption was examined for all covariates using a continuous time-varying predictor and generalized linear regression of scaled Schoenfeld residuals on functions of time and visually using graphical methods.^19,20 Kaplan-Meier curves were used to assess patient survival with respect to CMR analysis. Intra-and interobserver agreement was assessed by correlation coefficients and Bland-Altman plots.

RESULTS

Forty-nine patients with PAH were enrolled. Table 1 summarizes the demographic characteristics, RHC data, and differences in 2D and 3D CMR measures between control subjects and subjects with PAH. Table 2 shows the relationship between 2D and 3D measures by CMR and hemodynamic measures by RHC. Table 3 shows the differences between CMR and RHC findings in subjects with scleroderma-associated PAH (SSc-PAH) and idiopathic PAH (IPAH). Patients with IPAH had higher mPAP at RHC (49 ± 14 mmHg vs. 39 ± 11 mmHg; P = 0.01) and RV stroke work index (21.0 ± 6.2 mL/m² vs. 14.6 ± 4.6 mL/m²; P < 0.0001); all other RHC and CMR measures were not significantly different between the groups. To evaluate the accuracy of CMR volumetric measures, we compared RHC- and CMR-derived RVSVI (defined as (cardiac index/heart rate) × 1,000 body surface area), which were nearly identical (CMR, 38 ± 11; RHC, 36 ± 11; P = 0.33).

Table 1.

Demographic characteristics, hemodynamic test results, and cardiac magnetic resonance imaging (CMR) findings for patients with pulmonary arterial hypertension (PAH) and control subjects

Clinical characteristic	Control subjects (n = 18)	Patients with PAH (n = 49)	P
Age, years	51 ± 9	57 ± 12
Sex, no. (%)
Male	6 (33)	4 (8)
Female	12 (67)	45 (92)
Race, no. (%)
White	11 (61)	43 (88)
Black	5 (28)	5 (10)
Other	2 (11)	1 (2)
NYHA class, no. (%)
I	…	6 (12)
II	…	22 (45)
III	…	20 (41)
IV	…	1 (2)
Diagnosis, no. (%)
IPAH	…	17 (35)
SSc-PAH	…	32 (65)
Hemodynamic characteristics at right heart catheterization
HR, beats/min	NA	74 ± 9
Systolic BP, mmHg	NA	121 ± 21
Mean right atrial pressure, mmHg	NA	7 ± 4
Mean pulmonary arterial pressure, mmHg	NA	43 ± 13
Pulmonary capillary wedge pressure, mmHg	NA	10 ± 3
Cardiac index, L/min/m²	NA	2.7 ± 0.7
Pulmonary vascular resistance, WU	NA	7.9 ± 5.3
RV stroke work index, mL/m²	NA	16.8 ± 5.8
Stroke volume index, mL/m²	NA	35.9 ± 11
CMR measures
2D measures
TAPSE, mm	22 ± 5	15 ± 6	<0.01
RV fraction shortening, %	24 ± 4	16 ± 6	<0.01
RV fractional area change, %	24 ± 5	17 ± 6	<0.01
3D measures
Volume mass index	0.36 ± 0.1	0.58 ± 0.2	<0.01
RV end-diastolic mass index, g/m²	26.6 ± 6	37.5 ± 17	<0.01
RV ejection fraction, %	57.8 ± 6	45.2 ± 14	<0.01
RV end-diastolic volume index, mL/m²	76.0 ± 14	89.7 ± 33	0.09
RV stroke volume index	43.2 ± 5	38.1 ± 11	0.06
LV end-diastolic mass index, g/m²	74.4 ± 12	65.6 ± 18	0.05
LV ejection fraction, %	66.6 ± 6	64.2 ± 8	0.22
LV stroke volume index	43.6 ± 5	37 ± 11	<0.05
LV end-diastolic volume index, mL/m²	66.1 ± 10	56.9 ± 17	<0.05

Note: Data are mean value ± standard deviation (SD), unless otherwise indicated. BP: blood pressure; HR: heart rate; IPAH: idiopathic pulmonary arterial hypertension; LV: left ventricle; NA: not available; NYHA: New York Health Association; RV: right ventricle; SSc-PAH: scleroderma-associated PAH; TAPSE: tricuspid annular plane systolic excursion; WU: Wood units; 2D: bidiminsional; 3D: three diminsional.

Table 2.

Correlation between bidimensional (2D) and three-diminsional (3D) volumetric cardiac magnetic resonance imaging measures and invasive hemodynamic parameters

	2D measures			3D measures
Parameter	TAPSE	RVFS	RVFAC	RVEDVI	RVEF	VMI
PVRI	−0.52	−0.50	−0.54	0.17^a	−0.67	0.64
mPAP	−0.37	−0.42	−0.47	0.18^a	−0.52	0.53
RVSVI	0.49	0.40	0.40	−0.22^a	0.62	−0.50

Note: mPAP: mean pulmonary arterial pressure; PVRI: pulmonary vascular resistance index; RVEDVI: RV end-diastolic volume index; RVEF: RV ejection fraction; RVFAC: RV fraction area change; RVFS: RV fraction shortening; RVSVI: right ventricle stroke volume index; TAPSE: tricuspid annular plane systolic excursion; VMI: volume mass index.

P value was not significant.

Table 3.

Differences between cardiac magnetic resonance imaging (CMR) and right heart catheterization findings in subjects with scleroderma-related pulmonary arterial hypertension (SSc-PAH) and idiopathic pulmonary arterial hypertension (IPAH)

Variable	SSc-PAH (n = 32)	IPAH (n = 17)	P
CMR measure
TAPSE, mm	15.4 ± 5.6	14.6 ± 5.5	0.64
RVFS, %	17.2 ± 6.1	14.9 ± 5.6	0.22
RVFAC, %	17.8 ± 6.3	16.2 ± 5.5	0.36
RVEDVI, mL/m²	87.0 ± 34.7	94.7 ± 28.7	0.43
RVSVI, mL/m²	38.4 ± 11.0	37.5 ± 11.3	0.78
RVCO, L/min	4.8 ± 1.2	5.2 ± 1.6	0.28
RVEF, %	46.7 ± 13.6	42.4 ± 15	0.32
RVMI, g/m²	37.6 ± 17.8	37.3 ± 17.9	0.96
VMI, g/m²	0.55 ± 0.2	0.65 ± 0.3	0.18
Hemodynamic characteristic
mRAP, mmHg	7 ± 4	8 ± 5	0.36
mPAP, mmHg	39 ± 11	49 ± 14	0.01
PCWP, mmHg	9.8 ± 3.4	10.6 ± 3.4	0.40
CO, L/min	4.6 ± 1.6	5.3 ± 1.4	0.13
Cardiac index, L/min/m²	2.6 ± 0.8	2.7 ± 0.61	0.60
PVR, WU	7.9 ± 6.0	8.0 ± 3.7	0.98
RVSWI, mL/m²	14.6 ± 4.6	21.0 ± 6.2	<0.0001
RVSVI, mL/m²	34.2 ± 11.2	39.2 ± 10.1	0.13

Note: Data are mean value ± standard deviation. Boldface type indicates statistical significance. CO: cardiac output; mPAP: mean pulmonary artery pressure; mRAP: main right atrial pressure; PCWP: pulmonary capillary wedge pressure; PVR: pulmonary vascular resistance; RV: right ventricular; RVCO: RV cardiac output; RVEDVI: RV end-diastolic volume index; RVEF: RV ejection fraction; RVFAC: RV fraction area change; RVFS: RV fraction shortening; RVMI: RV mass index; RVSVI: RV stroke volume index; RVSWI: RV stroke work index; TAPSE: tricuspid annular plane systolic excursion; VMI: volume mass index; WU: Wood units.

CMR and survival

Individuals with PAH were followed up for a mean of 2.5 ± 1.6 years after CMR evaluation. Twenty-one subjects (43%) died during this follow-up period. Univariate Cox regression analyses were performed to evaluate the relative strength of association of demographic, hemodynamic, functional, and morphologic parameters and treatment status with survival. Table 4 shows the strength of association of mortality with study parameters expressed as continuous variables. There was a significant association between mortality and age, New York Heart Association class, RHC measures of PVR and cardiac index, and all 2D CMR measures of RVF (TAPSE, RVFAC, and RVFS). Interestingly, there was a trend but not a significant association between mortality and RVEF or VMI. ROC curves were used to generate thresholds with the optimal sensitivity and specificity for prediction of death. The optimal ROC-derived cutoff values for prediction of death were 18 mm for TAPSE, 16.7% for RVFS, 18.8% for RVFAC, 44.1% for RVEF, 0.55 for VMI, and 37.7 mL/m² for SVI. ROC curve analysis showed moderate-to-good prediction of mortality for TAPSE (area under the ROC curve [AUC], 0.77 [95% confidence interval (CI), 0.63–0.88]), RVFS (AUC, 0.75 [95% CI, 0.60–0.86]), RVFAC (AUC, 0.74 [95% CI, 0.60–0.86]), RVEF (AUC, 0.65 [95% CI, 0.50–0.79]), VMI (AUC, 0.59 [95% CI, 0.44–0.73]), PVR (AUC, 0.73 [95% CI, 0.58–0.84]), and SVI (AUC, 0.61 [95% CI, 0.46–0.74]). There were no significant differences between the AUC of TAPSE compared with RVEF (P = 0.14) and VMI (P = 0.06). The AUC of TAPSE was higher than that of SVI (P = 0.05).

Table 4.

Univariate Cox regression analyses to evaluate the relative strength of association of demographic, hemodynamic, functional, and morphologic parameters and treatment status with survival

Measure	Hazard ratio (95% confidence interval)	P
Age	1.02 (0.98–1.07)	0.04
Race	0.57 (0.18–1.80)	0.34
NYHA	2.27 (1.05–4.97)	0.04
TAPSE	0.91 (0.85–0.98)	0.02
RVFS	0.93 (0.87–0.99)	0.02
RVFAC	0.92 (0.86–0.98)	0.01
RVEF	0.98 (0.95–1.01)	0.20
LVEDVI	0.97 (0.94–1.00)	0.05
LVEF	1.03 (0.96–1.10)	0.40
LVSVI	0.98 (0.94–1.02)	0.34
VMI	1.98 (0.45–8.78)	0.37
Cardiac index	0.56 (0.29–1.09)	0.08
SVI	0.97 (0.93–1.01)	0.15
mPAP	1.02 (0.99–1.06)	0.14
PVR	1.06 (1.00–1.13)	0.03
Treatment status	1.3 (0.6–3.5)	0.45

Note: LVEDVI: left ventricular end-diastolic volume index; LVEF: left ventricular ejection fraction; LVSVI: left ventricular stroke volume index; mPAP: mean pulmonary artery pressure; PVR: pulmonary vascular resistance; NYHA: New York Health Association; RVEF: right ventricular ejection fraction; RVFAC: right ventricular fraction area change; RVFS: right ventricular fraction shortening; SVI: stroke volume index; TAPSE: tricuspid annular plane systolic excursion; VMI: volume mass index.

Kaplan-Meier survival curves demonstrated good survival prediction for TAPSE and SVI using ROC-derived cutoff values. However, there were no significant differences in survival between groups separated by ROC-derived cutoffs for RVEF and VMI (Fig. 2). Univariate and repeated bivariate Cox regression analysis of dichotomized variables are shown in Table 5 and Table 6. TAPSE remained a significant predictor of mortality after adjusting for potential confounders.

Table 5.

Univariate Cox regression analyses to evaluate the relative strength of association of demographic, hemodynamic, functional, and morphologic parameters and treatment status with survival with variables dichotomized according to the optimal receiver-operating characteristic curve threshold for prediction of death

Measure	Hazard ratio (95% confidence interval)	P
TAPSE ≤18 mm	4.8 (2.0–11.3)	0.005
RVFS ≤16.7%	3.6 (1.5–8.6)	0.004
RVFAC ≤18.8%	3.8 (1.6–9.2)	0.003
RVEF ≤40%	1.9 (0.8–4.9)	0.11
VMI ≥0.55	1.5 (0.6–3.6)	0.32
Cardiac index ≤2.5 L/min/m²	2.5 (1.0–6.1)	0.02
SVI ≤37.7 mL/m²	2.0 (0.9–4.9)	0.14

Note: RVEF: right ventricular ejection fraction; RVFAC: right ventricular fraction area change; RVFS: right ventricular fraction shortening; SVI: stroke volume index; TAPSE: tricuspid annular plane systolic excursion; VMI: volume mass index.

Table 6.

Repeated bivariate cox regression models for risk of death in the cohort

TAPSE ≤;18 mm	Hazard ratio (95% confidence interval)	P
Unadjusted	4.8 (2.0–11.3)	0.005
Adjusted for
Age	4.6 (1.4–15.7)	0.01
Race	5.0 (1.5–17.0)	0.01
NYHA	4.3 (1.3–14.7)	0.02
Treatment status	5.3 (1.5–18.3)	0.008
Disease type (IPAH vs. SSc-PAH)	5.8 (1.7–20.3)	0.006
RVFS	3.9 (0.9–17.9)	0.08
RVFAC	3.4 (0.8–14.4)	0.10
RVEF	4.5 (1.3–15.5)	0.02
LVEDVI	4.0 (1.1–14.5)	0.04
VMI	4.6 (1.3–15.8)	0.01
Cardiac index	4.2 (1.2–14.8)	0.03
SVI	4.4 (1.2–15.6)	0.02
PVR	4.1 (1.2–14.2)	0.03

Note: IPAH: ideopathic pulmonary arterial hypertension; LVEDVI: left ventricular end-diastolic volume index; NYHA: New York Health Association; PVR: pulmonary vascular resistance; RVEF: right ventricular ejection fraction; RVFAC: right ventricular fraction area change; RVFS: right ventricular fraction shortening; SSc-PAH: scleroderma-related pulmonary arterial hypertension; SVI: stroke volume index; TAPSE: tricuspid annular plane systolic excursion; VMI: volume mass index.

Figure 2.

Kaplan-Meier survival curves in a cohort of 49 patients with pulmonary arterial hypertension for both two-dimensional and three-dimensional cardiac magnetic resonance imaging (CMR) measures. A, Tricuspid annulus plane systolic excursion (TAPSE) of ≤18 mm (P < 0.01). B, Right ventricle (RV) ejection fraction (EF) of ≤40% (P = 0.11). C, Ventricular mass index (VMI) of ≥;0.55 (P = 0.32). D, Stroke volume index (SVI) of ≤;37.7 mL/m² (P = 0.14).

Treatment

Although all patients received PAH-specific therapy on the basis of our program treatment protocol during the follow-up observation period, 15 patients (31%) were treatment naive and 34 patients were receiving treatment within 3 months before baseline CMR. All treatment-naive patients had SSc-PAH, and 8 had TAPSE <18 mm. Of the 34 treated subjects, 19 received single-therapy (prostacyclin [PGI2], 1 patient; calcium channel blockers [CCB], 2 patients; endothelin receptor antagonist [ERA], 5 patients; phosphodiesterase inhibitor [PDEi], 11 patients); 10 patients received dual therapy (ERA + PDEi, 9 patients; PGI2 + PDEi, 1 patient); and 5 patients received triple therapy (CCB + ERA + PDEi, 2 patients; ERA + PGI2 + PDEi, 3 patients). Direct comparison of survival in treatment-naive versus treated patients showed no difference in survival for treated patients (median survival, 4.3 years) compared with treatment-naive patients (median survival, 3.1 years; P = 0.45, log-rank test). Significant differences in survival were seen in both treated and treatment-naive groups when subdivided by a TAPSE cutoff of 18 mm (Fig. 3).

Figure 3.

Kaplan-Meier survival curves for treated and treatment-naive patients with pulmonary arterial hypertension with respect to cardiac magnetic resonance imaging (CMR) categorized by tricuspid annulus plane systolic excursion (TAPSE) above and below 18 mm (P = 0.02). Patients who had TAPSE >18 mm demonstrate better survival than patients who had TAPSE ≤18 mm in both patient groups, irrespective of treatment status.

CMR reproducibility

Mean time imaging analysis for 2D measures was 90 ± 20 seconds per case compared with 25 ± 6 minutes for the 3D measures analysis. There was excellent intraobserver and interobserver correlation for TAPSE, RVFS, and RVFAC (intraobserver: r = 0.99, 0.97, and 0.96, respectively; interobserver: r = 0.96, 0.93, and 0.98, respectively). The intraobserver and interobserver correlation for RVED volume index, RVEF, and VMI were also excellent (intraobserver: r = 0.94, 0.97, and 0.96, respectively; interobserver: r = 0.92, 0.87, and 0.85, respectively). The interobserver Bland-Altman plot for TAPSE showed good agreement, with mean difference (bias) of 0.7 mm and 95% limits of agreement of 3.0 to −4.4 mm. The intraobserver Bland-Altman plot for TAPSE also showed good agreement, with mean difference of 0.1 mm and 95% limits of agreement of 2.3 to −2.0 mm (Fig. 4). RVFS and RVFAC had similarly good inter- and intraobserver reproducibility, although interobserver variability for RVFS was higher than that of TAPSE and RVFAC, and intraobserver variability was slightly higher with RVFAC. Reproducibility of 3D measures was assessed for VMI and RVEF using Bland-Altman plots. The intraobserver Bland-Altman plot for VMI showed a mean difference of 0.03 and 95% limits of agreement of 0.45 to −0.39. The interobserver Bland-Altman plot for VMI had a mean difference of 0.05 and 95% limit of agreement of 0.49 to −0.39 (Fig. 5). The intraobserver Bland-Altman plot for RVEF showed a mean difference of −0.4 and 95% limits of agreement of 5.9 to −6.7. The interobserver Bland-Altman plot for RVEF had a mean difference of −3.6 and 95% limits of agreement of 8.5 to −15.8 (Fig. 5).

Figure 4.

Bland-Altman plots demonstrate the intraobserver (tricuspid annulus plane systolic excursion [TAPSE] and TAPSE1) and interobserver (TAPSE and TAPSE2) agreement by two-dimensional assessment for TAPSE (mm), right ventricular (RV) fractional shortening (RVFS; %), and RV fractional area change (RVFAC; %). Horizontal solid line represents mean differences, and dashed line represents 95% limits of agreement. SD: standard deviation.

Figure 5.

Bland-Altman plots demonstrate the intra- and interobserver agreement by three-dimensional assessment for ventricular mass index (VMI) and right ventricular ejection fraction (RVEF; %). Horizontal solid line represents mean differences, and dashed line represents 95% limits of agreement. SD: standard deviation.

DISCUSSION

In this study, we found that TAPSE, RVFS, and RVFAC by CMR were associated with invasive measures of pulmonary hypertension severity and RVF in patients with PAH. These 2D measures are reproducible strong predictors of survival, with significant time saving over 3D analyses.

Our results confirm that RV dysfunction is reflected by both 2D and 3D measures and that each modality correlates well with invasive hemodynamic measures of PAH severity. In our cohort, there was an inverse relationship between TAPSE and RV afterload, demonstrated by higher PVR in subjects with reduced TAPSE. TAPSE was shown to correlate with reduced RVSVI, which is comparable to a previous study showing that longitudinal RV displacement by echocardiography predicted a reduced RVSVI in patients with PAH.⁷ This maladaptive RV response to elevated pulmonary arterial pressure is well described and is believed to be a major contributor to mortality in patients with IPAH and SSc-PAH.²

CMR and survival

TAPSE obtained by echocardiography and VMI obtained by CMR have previously been shown to be predictors of mortality in PAH.^7–9,13,14 We found that CMR-derived TAPSE measurements have similar prognostic value with improved survival prediction compared with VMI. An earlier study evaluating TAPSE by echocardiography has found an optimal cutoff point of 18 mm for survival prediction,⁷ which is identical to our cutoff point of 18 mm derived from ROC analysis. In addition, other studies that have evaluated the use of TAPSE by CMR to diagnose pulmonary hypertension and to monitor RVF in acquired heart diseases have used similar cutoff points, ranging from 15 to 18 mm.^15,21 In our cohort, TAPSE <18 mm was associated with a nearly 5-fold increased risk of death, whereas a 4-fold increased risk was associated with RVFS <16.7% and RVFAC <18.8%. This relationship persisted when adjusting for potential confounders such as age, treatment status, and disease type. Surprisingly, the association of TAPSE with survival remained even when controlling for invasive and noninvasive morphologic and hemodynamic measures previously demonstrated to predict outcome in PAH, such as RVEF, cardiac index, and SVI. Not surprisingly, TAPSE was also strongly associated with RVF, assessed by both RHC and CMR parameters.

Similar to earlier studies, we found a significant correlation of VMI with mPAP and PVR, both of which are invasive measures of PAH severity.^14,22 Our VMI cutoff point of 0.55 is similar to that identified in earlier studies by Swift et al.²³ (0.4), Hagger et al.¹⁴ (0.56), and Saba et al.²² (0.6). Interestingly, in contrast to the earlier studies by Rajaram et al.¹³ and Hagger et al.,¹⁴ we found only a trend toward increased VMI and survival, which did not reach statical significance (P = 0.32). This could be due to differences in patient populations between our studies and a greater number of deaths in our cohort. Rajaram et al.¹³ examined a cohort of patients with a variety of connective tissue diseases, whereas our patients had either IPAH or SSc-PAH. In the article by Hagger et al.,¹⁴ VMI was examined in a cohort at risk for PAH, not a cohort of patients with established PAH, as in our study. A strong association between VMI and mortality in the study by Hagger et al.¹⁴ would be expected, because increased VMI identified those subjects with PAH who are known to have higher mortality.¹⁴

Treatment

Even though some patients were treatment naive at the time of CMR, all patients received treatment during the observation period. Therefore, it is not surprising that there were no significant differences in survival between groups. As expected, when we grouped our cohort according to whether they were treatment naive or treated and then stratified them by whether they had TAPSE above or below 18 mm, TAPSE remained a powerful survival predictor in both groups. Furthermore, TAPSE remained significantly associated with outcome when controlling for treatment status.

Reproducibility

Analysis of 2D measures by CMR took significantly less time than did analysis of 3D measurements (one-tenth of the time). These findings suggest the unique value of using TAPSE for rapid prognostic assessment in PAH compared with more labor-intensive 3D measures. The interobserver reproducibility between the 2D CMR measures was high with excellent agreement, consistent with earlier reports.^12,15,21 Bland-Altman plots demonstrated low inter-and intraobserver variability, with limits of agreement (1.96 SDs) ±3–4 units from the mean for each measurement (±3.7 mm for TAPSE, ±4.2% for RVFS, ±3.5% for RVFAC). Intra- and interobserver agreement was slightly better for 2D measures than for 3D measures. Correlation coefficients for both RVEF and VMI were high; however, the 95% limits of agreement were broad (±0.44 from the mean for VMI, ±12.2% for RVEF), indicating a higher degree of variability for these measures. Although it is considered the gold standard for volumetric analysis for the RV, manual contouring of the RV cavity by CMR is more challenging than manual contouring of the LV cavity, as reflected by increased measurement variability.²⁴ This is in part attributable to differences in contouring of the highly trabeculated RV endocardial surface or differences in identification of RV chamber on short-axis views near the tricuspid annulus.²⁴

Our study has limitations. First, there was a relatively limited sample size; therefore, conclusions about the relative strengths of survival prediction for various measures are made with caution. Second, 71% of deaths occurred among patients with SSc-PAH; therefore, it is uncertain whether the ability of CMR TAPSE to predict survival is generalizable to all patients with PAH or has prognostic significance only for patients with SSc-PAH. It is not surprising that the deaths occurred in the SSc-PAH group, given the known higher mortality rates associated with SSc-PAH compared with IPAH.²⁵ However, an earlier echocardiographic study demonstrated the prognostic significance of TAPSE in a mixed cohort of patients with pulmonary hypertension,⁶ and therefore, we suspect that the CMR results would have similar generalizability to all patients with PAH.

We did not measure transverse shortening of the RV, a 2D marker that has recently been suggested to have greater value for monitoring response to therapy and survival prediction than longitudinal shortening assessed by TAPSE.²⁶ However, we would suspect that changes in transverse shortening in our cohort would be reflected by differences in RVFAC. In our cohort, RVFAC was a strong predictor of survival, and its performance was similar to that of TAPSE.

Given the similarities between magnetic resonance imaging and echocardiography for survival prediction, one might argue against the merit of performing expensive and time-consuming CMR examinations for patients with PAH. Echocardiograms were unavailable to perform a direct head-to-head comparison with CMR for efficacy of survival prediction in our cohort; therefore, it is uncertain whether there are any differences between the two modalities in this respect. However, we believe CMR has several advantages, including more reliable evaluation of the RV due to a lack of user dependence, need for adequate imaging windows, and the ability to assess additional novel imaging biomarkers of RV performance. These include RV strain, pulmonary transit time, myocardial perfusion, and myocardial delayed enhancement.^27–31 The relative importance of these biomarkers for treatment decisions and survival prediction remains an active area of research.

In this study, we evaluated outcomes in a cohort of patients that included patients receiving treatment and patients who were treatment naive. This is a potential limitation, because the imaging parameters discussed could have different prognostic ability in these two groups. Multivariate analysis was not performed for treated and treatment-naive subgroups separately because of the small sample size. However, we think it is probable that 2D measures of RVF are prognostically important in both groups, as indicated by Kaplan-Meier survival analysis performed separately for both groups (Fig. 3).

Conclusion

TAPSE by CMR is a simplified measure of RVF that can be obtained with high interobserver agreement and in significantly less time compared with labor-intensive 3D measures, allowing rapid assessment of RVF and patient prognosis. CMR measures of RVF and morphology are strong, independent predictors of survival in PAH and correlate with invasive measures of pulmonary hemodynamics.

Footnotes

ACKNOWLEDGMENTS

This work has been presented in part at the 98th Scientific Assembly and Annual Meeting of the Radiological Society of North America (Chicago, IL; November 2012), abstract 12032983; the annual meeting of the American Roentgen Ray Society (Washington, DC; April 2013), abstract E061; and the Royal Australian and New Zealand College of Radiologists’ 64th Annual Scientific Meeting (Auckland, New Zealand; October 2013), abstract R-0171.

Source of Support: This work was supported by National Institutes of Health/National Heart, Lung, and Blood Institute grants P01HL84946-01 and R01 HL114910.

Conflict of Interest: IRK has received grants from Bayer Healthcare, Siemens Medical, and Bracco Diagnostics. PMH has received honoraria for consultations and/or being on a scientific advisory board for Gilead, Pfizer, Novartis, and Merck, activities which are irrelevant to the current work.

References

Bogaard

Abe

VonkNoordegraaf

Voelkel

. The right ventricle under pressure: cellular and molecular mechanisms of right-heart failure in pulmonary hypertension. Chest 2009;135:794–804.

Campo

Mathai

Le Pavec

Zaiman

Hummers

Boyce

Housten

. Hemodynamic predictors of survival in scleroderma-related pulmonary arterial hypertension. Am J Respir Crit Care Med 2010;182:252–260.

van Wolferen

Marcus

Boonstra

Marques

Bronzwaer

Spreeuwenberg

Postmus

Vonk-Noordegraaf

. Prognostic value of right ventricular mass, volume, and function in idiopathic pulmonary arterial hypertension. Eur Heart J 2007;28:1250–1257.

van de Veerdonk

Kind

Marcus

Mauritz

Heymans

Bogaard

Boonstra

Marques

Westerhof

Vonk-Noordegraaf

. Progressive right ventricular dysfunction in patients with pulmonary arterial hypertension responding to therapy. J Am Coll Cardiol 2011;58:2511–2519.

Haddad

Hunt

Rosenthal

Murphy

. Right ventricular function in cardiovascular disease, part I: anatomy, physiology, aging, and functional assessment of the right ventricle. Circulation 2008;117: 1436–1448.

Kaul

Tei

Hopkins

Shah

. Assessment of right ventricular function using two-dimensional echocardiography. Am Heart J 1984;107: 526–531.

Forfia

Fisher

Mathai

Housten-Harris

Hemnes

Borlaug

Chamera

. Tricuspid annular displacement predicts survival in pulmonary hypertension. Am J Respir Crit Care Med 2006; 174:1034–1041.

Mathai

Sibley

Forfia

Mudd

Fisher

Tedford

Lechtzin

. Tricuspid annular plane systolic excursion is a robust outcome measure in systemic sclerosis–associated pulmonary arterial hypertension. J Rheumatol 2011;38:2410–2418.

Ghio

Klersy

Magrini

D'Armini

Scelsi

Raineri

Pasotti

Serio

Campana

Vigano

. Prognostic relevance of the echocardiographic assessment of right ventricular function in patients with idiopathic pulmonary arterial hypertension. Int J Cardiol 2010;140:272–278.

10.

McLaughlin

Archer

Badesch

Barst

Farber

Lindner

Mathier

; American College of Cardiology Foundation Task Force on Expert Consensus Documents, American Heart Association, American College of Chest Physicians, American Thoracic Society I, Pulmonary Hypertension Association . ACCF/AHA 2009 expert consensus document on pulmonary hypertension a report of the American College of Cardiology Foundation Task Force on Expert Consensus Documents and the American Heart Association developed in collaboration with the American College of Chest Physicians; American Thoracic Society, Inc.; and the Pulmonary Hypertension Association. J Am Coll Cardiol 2009;53:1573–1619.

11.

Kjaergaard

Petersen

Kjaer

Schaadt

Hassager

. Evaluation of right ventricular volume and function by 2D and 3D echocardiography compared to MRI. Eur J Echocardiogr 2006;7:430–438.

12.

van der Zwaan

Geleijnse

McGhie

Boersma

Helbing

Meijboom

Roos-Hesselink

. Right ventricular quantification in clinical practice: two-dimensional vs. three-dimensional echocardiography compared with cardiac magnetic resonance imaging. Eur J Echocardiogr 2011;12:656–664.

13.

Rajaram

Swift

Capener

Elliot

Condliffe

Davies

Hill

. Comparison of the diagnostic utility of cardiac magnetic resonance imaging, computed tomography, and echocardiography in assessment of suspected pulmonary arterial hypertension in patients with connective tissue disease. J Rheumatol 2012;39:1265–1274.

14.

Hagger

Condliffe

Woodhouse

Elliot

Armstrong

Davies

Hill

Akil

Wild

Kiely

. Ventricular mass index correlates with pulmonary artery pressure and predicts survival in suspected systemic sclerosis-associated pulmonary arterial hypertension. Rheumatology (Oxford) 2009;48:1137–1142.

15.

Nijveldt

Germans

McCann

Beek

van Rossum

. Semi-quantitative assessment of right ventricular function in comparison to a 3D volumetric approach: a cardiovascular magnetic resonance study. Eur Radiol 2008;18:2399–2405.

16.

Speiser

Hirschberger

Pilz

Heer

Sievers

Strasser

Schoen

. Tricuspid annular plane systolic excursion assessed using MRI for semi-quantification of right ventricular ejection fraction. Br J Radiol 2012;85:e716–e721.

17.

Alfakih

Plein

Thiele

Jones

Ridgway

Sivananthan

. Normal human left and right ventricular dimensions for MRI as assessed by turbo gradient echo and steady-state free precession imaging sequences. J Magn Reson Imaging 2003;17:323–329.

18.

Stevens

Garcia-Alvarez

Sahni

Garcia

Fuster

Sanz

. RV dysfunction in pulmonary hypertension is independently related to pulmonary artery stiffness. JACC Cardiovasc Imaging 2012;5:378–387.

19.

Szklo

Nieton

, eds. Epidemiology: beyond the basics. Gaithersburg, MD: Aspen, 2000.

20.

Hess

. Graphical methods for assessing violations of the proportional hazards assumption in Cox regression. Stat Med 1995;14:1707–1723.

21.

Caudron

Fares

Vivier

Lefebvre

Petitjean

Dacher

. Diagnostic accuracy and variability of three semi-quantitative methods for assessing right ventricular systolic function from cardiac MRI in patients with acquired heart disease. Eur Radiol 2011;21:2111–2120.

22.

Saba

Foster

Cockburn

Cowan

Peacock

. Ventricular mass index using magnetic resonance imaging accurately estimates pulmonary artery pressure. Eur Respir J 2002;20:1519–1524.

23.

Swift

Rajaram

Condliffe

Capener

Hurdman

Elliot

Wild

Kiely

. Diagnostic accuracy of cardiovascular magnetic resonance imaging of right ventricular morphology and function in the assessment of suspected pulmonary hypertension results from the ASPIRE registry. J Cardiovasc Magn Reson 2012;14:40-429X-14-40.

24.

Grothues

Moon

Bellenger

Smith

Klein

Pennell

. Interstudy reproducibility of right ventricular volumes, function, and mass with cardiovascular magnetic resonance. Am Heart J 2004;147: 218–223.

25.

Fisher

Mathai

Champion

Girgis

Housten-Harris

Hummers

Krishnan

Wigley

Hassoun

. Clinical differences between idiopathic and scleroderma-related pulmonary hypertension. Arthritis Rheum 2006;54:3043–3050.

26.

Mauritz

Kind

Marcus

Bogaard

van de Veerdonk

Postmus

Boonstra

Westerhof

Vonk-Noordegraaf

. Progressive changes in right ventricular geometric shortening and long-term survival in pulmonary arterial hypertension. Chest 2012;141:935–943.

27.

Shehata

Harouni

Skrok

Basha

Boyce

Lechtzin

Mathai

. Regional and global biventricular function in pulmonary arterial hypertension: a cardiac MR imaging study. Radiology 2013; 266:114–122.

28.

Skrok

Shehata

Mathai

Girgis

Zaiman

Mudd

Boyce

. Pulmonary arterial hypertension: MR imaging-derived first-pass bolus kinetic parameters are biomarkers for pulmonary hemodynamics, cardiac function, and ventricular remodeling. Radiology 2012; 263:678–687.

29.

Vogel-Claussen

Skrok

Shehata

Singh

Sibley

Boyce

Lechtzin

. Right and left ventricular myocardial perfusion reserves correlate with right ventricular function and pulmonary hemodynamics in patients with pulmonary arterial hypertension. Radiology 2011;258:119–127.

30.

Shehata

Lossnitzer

Skrok

Boyce

Lechtzin

Mathai

Girgis

. Myocardial delayed enhancement in pulmonary hypertension: pulmonary hemodynamics, right ventricular function, and remodeling. AJR Am J Roentgenol 2011;196:87–94.

31.

Ibrahim

White

. Cardiovascular magnetic resonance for the assessment of pulmonary arterial hypertension: toward a comprehensive CMR exam. Magn Reson Imaging 2012;30:1047–1058.